Robert J. Hardwick

A worldwide analysis of beta-defensin copy number variation suggests recent selection of a high-expressing DEFB103 gene copy in East Asia.

Beta‐defensins are a family of multifunctional genes with roles in defense against pathogens, reproduction, and pigmentation. In humans, six beta‐defensin genes are clustered in a repeated region which is copy‐number variable (CNV) as a block, with a diploid copy number between 1 and 12. The role in host defense makes the evolutionary history of this CNV particularly interesting, because morbidity due to infectious disease is likely to have been an important selective force in human evolution, and to have varied between geographical locations. Here, we show CNV of the beta‐defensin region in chimpanzees, and identify a beta‐defensin block in the human lineage that contains rapidly evolving noncoding regulatory sequences. We also show that variation at one of these rapidly evolving sequences affects expression levels and cytokine responsiveness of DEFB103, a key inhibitor of influenza virus fusion at the cell surface. A worldwide analysis of beta‐defensin CNV in 67 populations shows an unusually high frequency of high‐DEFB103‐expressing copies in East Asia, the geographical origin of historical and modern influenza epidemics, possibly as a result of selection for increased resistance to influenza in this region. Hum Mutat 32:743–750, 2011.

Genetic Architecture of Ambulatory Blood Pressure in the General Population. Insights From Cardiovascular Gene-Centric Array

Genetic determinants of blood pressure are poorly defined. We undertook a large-scale, gene-centric analysis to identify loci and pathways associated with ambulatory systolic and diastolic blood pressure. We measured 24-hour ambulatory blood pressure in 2020 individuals from 520 white European nuclear families (the Genetic Regulation of Arterial Pressure of Humans in the Community Study) and genotyped their DNA using the Illumina HumanCVD BeadChip array, which contains ≈50 000 single nucleotide polymorphisms in >2000 cardiovascular candidate loci. We found a strong association between rs13306560 polymorphism in the promoter region of MTHFR and CLCN6 and mean 24-hour diastolic blood pressure; each minor allele copy of rs13306560 was associated with 2.6 mm Hg lower mean 24-hour diastolic blood pressure (P=1.2×10−8). rs13306560 was also associated with clinic diastolic blood pressure in a combined analysis of 8129 subjects from the Genetic Regulation of Arterial Pressure of Humans in the Community Study, the CoLaus Study, and the Silesian Cardiovascular Study (P=5.4×10−6). Additional analysis of associations between variants in gene ontology-defined pathways and mean 24-hour blood pressure in the Genetic Regulation of Arterial Pressure of Humans in the Community Study showed that cell survival control signaling cascades could play a role in blood pressure regulation. There was also a significant overrepresentation of rare variants (minor allele frequency: <0.05) among polymorphisms showing at least nominal association with mean 24-hour blood pressure indicating that a considerable proportion of its heritability may be explained by uncommon alleles. Through a large-scale gene-centric analysis of ambulatory blood pressure, we identified an association of a novel variant at the MTHFR/CLNC6 locus with diastolic blood pressure and provided new insights into the genetic architecture of blood pressure.

The effects of in utero irradiation on mutation induction and transgenerational instability in mice

Epidemiological evidence suggests that the deleterious effects of prenatal irradiation can manifest during childhood, resulting in an increased risk of leukaemia and solid cancers after birth. However, the mechanisms underlying the long-term effects of foetal irradiation remain poorly understood. This study was designed to analyse the impact of in utero irradiation on mutation rates at expanded simple tandem repeat (ESTR) DNA loci in directly exposed mice and their first-generation (F(1)) offspring. ESTR mutation frequencies in the germline and somatic tissues of male and female mice irradiated at 12 days of gestation remained highly elevated during adulthood, which was mainly attributed to a significant increase in the frequency of singleton mutations. The prevalence of singleton mutations in directly exposed mice suggests that foetal irradiation results in genomic instability manifested both in utero and during adulthood. The frequency of ESTR mutation in the F(1) offspring of prenatally irradiated male mice was equally elevated across all tissues, which suggests that foetal exposure results in transgenerational genomic instability. In contrast, maternal in utero exposure did not affect the F(1) stability. Our data imply that the passive erasure of epigenetic marks in the maternal genome can diminish the transgenerational effects of foetal irradiation and therefore provide important clues to the still unknown mechanisms of radiation-induced genomic instability. The results of this study offer a plausible explanation for the effects of in utero irradiation on the risk of leukaemia and solid cancers after birth.

Determination of Beta-Defensin Genomic Copy Number in Different Populations: A Comparison of Three Methods

Background There have been conflicting reports in the literature on association of gene copy number with disease, including CCL3L1 and HIV susceptibility, and β-defensins and Crohns disease. Quantification of precise gene copy numbers is important in order to define any association of gene copy number with disease. At present, real-time quantitative PCR (QPCR) is the most commonly used method to determine gene copy number, however the Paralogue Ratio Test (PRT) is being used in more and more laboratories. Findings In this study we compare a Pyrosequencing-based Paralogue Ratio Test (PPRT) for determining beta-defensin gene copy number with two currently used methods for gene copy number determination, QPCR and triplex PRT by typing five different cohorts (UK, Danish, Portuguese, Ghanaian and Czech) of DNA from a total of 576 healthy individuals. We found a systematic measurement bias between DNA cohorts revealed by QPCR, but not by the PRT-based methods. Using PRT, copy number ranged from 2 to 9 copies, with a modal copy number of 4 in all populations. Conclusions QPCR is very sensitive to quality of the template DNA, generating systematic biases that could produce false-positive or negative disease associations. Both triplex PRT and PPRT do not show this systematic bias, and type copy number within the correct range, although triplex PRT appears to be a more precise and accurate method to type beta-defensin copy number.

Large-Scale Candidate Gene Analysis of HDL Particle Features

Background HDL cholesterol (HDL-C) is an established marker of cardiovascular risk with significant genetic determination. However, HDL particles are not homogenous, and refined HDL phenotyping may improve insight into regulation of HDL metabolism. We therefore assessed HDL particles by NMR spectroscopy and conducted a large-scale candidate gene association analysis. Methodology/Principal Findings We measured plasma HDL-C and determined mean HDL particle size and particle number by NMR spectroscopy in 2024 individuals from 512 British Caucasian families. Genotypes were 49,094 SNPs in >2,100 cardiometabolic candidate genes/loci as represented on the HumanCVD BeadChip version 2. False discovery rates (FDR) were calculated to account for multiple testing. Analyses on classical HDL-C revealed significant associations (FDR<0.05) only for CETP (cholesteryl ester transfer protein; lead SNP rs3764261: p = 5.6*10−15) and SGCD (sarcoglycan delta; rs6877118: p = 8.6*10−6). In contrast, analysis with HDL mean particle size yielded additional associations in LIPC (hepatic lipase; rs261332: p = 6.1*10−9), PLTP (phospholipid transfer protein, rs4810479: p = 1.7*10−8) and FBLN5 (fibulin-5; rs2246416: p = 6.2*10−6). The associations of SGCD and Fibulin-5 with HDL particle size could not be replicated in PROCARDIS (n = 3,078) and/or the Womens Genome Health Study (n = 23,170). Conclusions We show that refined HDL phenotyping by NMR spectroscopy can detect known genes of HDL metabolism better than analyses on HDL-C.

β-defensin Genomic Copy Number Is Associated With HIV Load and Immune Reconstitution in Sub-Saharan Africans

AIDS, caused by the retrovirus human immunodeficiency virus (HIV), is the leading cause of death of economically active people (age, 15-59 years) in sub-Saharan Africa. The host genetic variability of immune response to HIV and immune reconstitution following initiation of highly active antiretroviral therapy (HAART) is poorly understood. Here we focused on copy number variation of the β-defensin genes, which have been shown to have anti-HIV activity, and are important chemoattractants for Th17 lymphocytes via the chemokine receptor CCR6. We determined β-defensin gene copy number for 1002 Ethiopian and Tanzanian patients. We show that higher β-defensin copy number variation is associated with increased HIV load prior to HAART (P=.005) and poor immune reconstitution following initiation of HAART (P=.003). We suggest a model where variable amounts of β-defensin expression by mucosal cells, due to gene copy number variation, alters the efficacy of recruitment of Th17 lymphocytes to the site of infection, altering the dynamics of infection.

Evolutionary History of Copy-Number-Variable Locus for the Low-Affinity Fcγ Receptor: Mutation Rate, Autoimmune Disease, and the Legacy of Helminth Infection

Both sequence variation and copy-number variation (CNV) of the genes encoding receptors for immunoglobulin G (Fcγ receptors) have been genetically and functionally associated with a number of autoimmune diseases. However, the molecular nature and evolutionary context of this variation is unknown. Here, we describe the structure of the CNV, estimate its mutation rate and diversity, and place it in the context of the known functional alloantigen variation of these genes. Deletion of Fcγ receptor IIIB, associated with systemic lupus erythematosus, is a result of independent nonallelic homologous recombination events with a frequency of approximately 0.1%. We also show that pathogen diversity, in particular helminth diversity, has played a critical role in shaping the functional variation at these genes both between mammalian species and between human populations. Positively selected amino acids are involved in the interaction with IgG and include some amino acids that are known polymorphic alloantigens in humans. This supports a genetic contribution to the hygiene hypothesis, which states that past evolution in the context of helminth diversity has left humans with an array of susceptibility alleles for autoimmune disease in the context of a helminth-free environment. This approach shows the link between pathogens and autoimmune disease at the genetic level and provides a strategy for interrogating the genetic variation underlying autoimmune-disease risk and infectious-disease susceptibility.

Age-Related Accumulation of Mutations Supports a Replication-Dependent Mechanism of Spontaneous Mutation at Tandem Repeat DNA Loci in Mice

Expanded simple tandem repeat (ESTR) loci belong to the class of highly unstable loci in the mouse genome. The mechanisms underlying the very high spontaneous instability at these loci still remain poorly understood. Using single-molecule polymerase chain reaction, here we have compared the pattern of mutation accumulation in tissues with different proliferation capacities in male mice of age 12, 26, 48, and 96 weeks. In the nonproliferating brain, we did not observe any measurable age-related accumulation of ESTR mutations. In contrast, a highly elevated frequency of ESTR mutation was detected in the sperm samples taken from old mice; similar changes were also observed in the bone marrow tissue. The spectra of ESTR mutations accumulated in all tissues of young and old mice did not significantly differ. Taken together, these data clearly imply that spontaneous ESTR mutations occur almost exclusively in replication-proficient cells. To gain further insights into the mechanisms of ESTR mutation, we developed a stochastic model of age-related mutation accumulation. The observed spectra of ESTR mutants accumulated in the brain and sperm were fairly accurately approximated assuming the values of ESTR mutation rate, ranging from 0.01 to 0.04 per cell division. As these estimates dramatically exceed those for protein-coding genes and microsatellite loci, our data therefore suggest that ESTRs represent one of the most unstable loci in the mammalian genome. The results of our study also imply that ESTR loci can be regarded as a class of expanded microsatellites, with the mechanism of spontaneous mutation most probably attributed to replication slippage.

CCL3L1 copy number, HIV load, and immune reconstitution in sub-Saharan Africans

BackgroundThe role of copy number variation of the CCL3L1 gene, encoding MIP1α, in contributing to the host variation in susceptibility and response to HIV infection is controversial. Here we analyse a sub-Saharan African cohort from Tanzania and Ethiopia, two countries with a high prevalence of HIV-1 and a high co-morbidity of HIV with tuberculosis.MethodsWe use a form of quantitative PCR called the paralogue ratio test to determine CCL3L1 gene copy number in 1134 individuals and validate our copy number typing using array comparative genomic hybridisation and fiber-FISH.ResultsWe find no significant association of CCL3L1 gene copy number with HIV load in antiretroviral-naïve patients prior to initiation of combination highly active anti-retroviral therapy. However, we find a significant association of low CCL3L1 gene copy number with improved immune reconstitution following initiation of highly active anti-retroviral therapy (p = 0.012), replicating a previous study.ConclusionsOur work supports a role for CCL3L1 copy number in immune reconstitution following antiretroviral therapy in HIV, and suggests that the MIP1α -CCR5 axis might be targeted to aid immune reconstitution.

Automated design of paralogue ratio test assays for the accurate and rapid typing of copy number variation

Motivation: Genomic copy number variation (CNV) can influence susceptibility to common diseases. High-throughput measurement of gene copy number on large numbers of samples is a challenging, yet critical, stage in confirming observations from sequencing or array Comparative Genome Hybridization (CGH). The paralogue ratio test (PRT) is a simple, cost-effective method of accurately determining copy number by quantifying the amplification ratio between a target and reference amplicon. PRT has been successfully applied to several studies analyzing common CNV. However, its use has not been widespread because of difficulties in assay design. Results: We present PRTPrimer (www.prtprimer.org) software for automated PRT assay design. In addition to stand-alone software, the web site includes a database of pre-designed assays for the human genome at an average spacing of 6 kb and a web interface for custom assay design. Other reference genomes can also be analyzed through local installation of the software. The usefulness of PRTPrimer was tested within known CNV, and showed reproducible quantification. This software and database provide assays that can rapidly genotype CNV, cost-effectively, on a large number of samples and will enable the widespread adoption of PRT. Availability: PRTPrimer is available in two forms: a Perl script (version 5.14 and higher) that can be run from the command line on Linux systems and as a service on the PRTPrimer web site (www.prtprimer.org). Contact: cjt14@le.ac.uk Supplementary Information: Supplementary data are available at Bioinformatics online.