A comparison of software for analysis of rare and common short tandem repeat (STR) variation using human genome sequences from clinical and population-based samples

Short tandem repeat (STR) variation is an often overlooked source of variation between genomes. STRs comprise about 3% of the human genome and are highly polymorphic. Some cause Mendelian disease, and others affect gene expression. Their contribution to common disease is not well-understood, but recent software tools designed to genotype STRs using short read sequencing data will help address this. Here, we compare software that genotypes common STRs and rarer STR expansions genome-wide, with the aim of applying them to population-scale genomes. By using the Genome-In-A-Bottle (GIAB) consortium and 1000 Genomes Project short-read sequencing data, we compare performance in terms of sequence length, depth, computing resources needed, genotyping accuracy and number of STRs genotyped. To ensure broad applicability of our findings, we also measure genotyping performance against a set of genomes from clinical samples with known STR expansions, and a set of STRs commonly used for forensic identification. We find that HipSTR, ExpansionHunter and GangSTR perform well in genotyping common STRs, including the CODIS 13 core STRs used for forensic analysis. GangSTR and ExpansionHunter outperform HipSTR for genotyping call rate and memory usage. ExpansionHunter denovo (EHdn), STRling and GangSTR outperformed STRetch for detecting expanded STRs, and EHdn and STRling used considerably less processor time compared to GangSTR. Analysis on shared genomic sequence data provided by the GIAB consortium allows future performance comparisons of new software approaches on a common set of data, facilitating comparisons and allowing researchers to choose the best software that fulfils their needs.

Gene fusions during the early evolution of mesothelioma correlate with impaired DNA repair and Hippo pathways

Malignant pleural mesothelioma (MPM), a rare cancer a long latency period (up to 40 years) between asbestos exposure and disease presentation. The mechanisms coupling asbestos to recurrent somatic alterations are poorly defined. Gene fusions arising through genomic instability may create novel drivers during early MPM evolution. We explored the gene fusions that occurred early in the evolutionary history of the tumor. We conducted multiregional whole exome sequencing (WES) of 106 samples from 20 patients undergoing pleurectomy decortication and identified 24 clonal nonrecurrent gene fusions, three of which were novel (FMO9P-OR2W5, GBA3, and SP9). The number of early gene fusion events detected varied from zero to eight per tumor, and presence of gene fusions was associated with clonal losses involving the Hippo pathway genes and homologous recombination DNA repair genes. Fusions involved known tumor suppressors BAP1, MTAP, and LRP1B, and a clonal oncogenic fusion involving CACNA1D-ERC2, PARD3B-NT5DC2, and STAB2-NT5DC2 fusions were also identified as clonal fusions. Gene fusions events occur early during MPM evolution. Individual fusions are rare as no recurrent truncal fusions event were found. This suggests the importance of early disruption of these pathways in generating genomic rearrangements resulting in potentially oncogenic gene fusions.

Association study of human leukocyte antigen variants and idiopathic pulmonary fibrosis

Introduction

Idiopathic pulmonary fibrosis (IPF) is a chronic interstitial pneumonia marked by progressive lung fibrosis and a poor prognosis. Recent studies have highlighted the potential role of infection in the pathogenesis of IPF, and a prior association of the HLA-DQB1 gene with idiopathic fibrotic interstitial pneumonia (including IPF) has been reported. Owing to the important role that the human leukocyte antigen (HLA) region plays in the immune response, here we evaluated if HLA genetic variation was associated specifically with IPF risk.

Methods

We performed a meta-analysis of associations of the HLA region with IPF risk in individuals of European ancestry from seven independent case-control studies of IPF (comprising 5159 cases and 27 459 controls, including a prior study of fibrotic interstitial pneumonia). Single nucleotide polymorphisms, classical HLA alleles and amino acids were analysed and signals meeting a region-wide association threshold of p<4.5×10-4 and a posterior probability of replication >90% were considered significant. We sought to replicate the previously reported HLA-DQB1 association in the subset of studies independent of the original report.

Results

The meta-analysis of all seven studies identified four significant independent single nucleotide polymorphisms associated with IPF risk. However, none met the posterior probability for replication criterion. The HLA-DQB1 association was not replicated in the independent IPF studies.

Conclusion

Variation in the HLA region was not consistently associated with risk in studies of IPF. However, this does not preclude the possibility that other genomic regions linked to the immune response may be involved in the aetiology of IPF.

Genome-wide association study of thyroid-stimulating hormone highlights new genes, pathways and associations with thyroid disease

Thyroid hormones play a critical role in regulation of multiple physiological functions and thyroid dysfunction is associated with substantial morbidity. Here, we use electronic health records to undertake a genome-wide association study of thyroid-stimulating hormone (TSH) levels, with a total sample size of 247,107. We identify 158 novel genetic associations, more than doubling the number of known associations with TSH, and implicate 112 putative causal genes, of which 76 are not previously implicated. A polygenic score for TSH is associated with TSH levels in African, South Asian, East Asian, Middle Eastern and admixed American ancestries, and associated with hypothyroidism and other thyroid disease in South Asians. In Europeans, the TSH polygenic score is associated with thyroid disease, including thyroid cancer and age-of-onset of hypothyroidism and hyperthyroidism. We develop pathway-specific genetic risk scores for TSH levels and use these in phenome-wide association studies to identify potential consequences of pathway perturbation. Together, these findings demonstrate the potential utility of genetic associations to inform future therapeutics and risk prediction for thyroid diseases.

Identification and characterisation of a rare MTTP variant underlying hereditary non-alcoholic fatty liver disease

Background & Aims

Non-alcoholic fatty liver disease (NAFLD) is a complex trait with an estimated prevalence of 25% globally. We aimed to identify the genetic variant underlying a four-generation family with progressive NAFLD leading to cirrhosis, decompensation, and development of hepatocellular carcinoma in the absence of common risk factors such as obesity and type 2 diabetes.

Methods

Exome sequencing and genome comparisons were used to identify the likely causal variant. We extensively characterised the clinical phenotype and post-prandial metabolic responses of family members with the identified novel variant in comparison with healthy non-carriers and wild-type patients with NAFLD. Variant-expressing hepatocyte-like cells (HLCs) were derived from human-induced pluripotent stem cells generated from homozygous donor skin fibroblasts and restored to wild-type using CRISPR-Cas9. The phenotype was assessed using imaging, targeted RNA analysis, and molecular expression arrays.

Results

We identified a rare causal variant c.1691T>C p.I564T (rs745447480) in MTTP, encoding microsomal triglyceride transfer protein (MTP), associated with progressive NAFLD, unrelated to metabolic syndrome and without characteristic features of abetalipoproteinaemia. HLCs derived from a homozygote donor had significantly lower MTP activity and lower lipoprotein ApoB secretion than wild-type cells, while having similar levels of MTP mRNA and protein. Cytoplasmic triglyceride accumulation in HLCs triggered endoplasmic reticulum stress, secretion of pro-inflammatory mediators, and production of reactive oxygen species.

Conclusions

We have identified and characterised a rare causal variant in MTTP, and homozygosity for MTTP p.I564T is associated with progressive NAFLD without any other manifestations of abetalipoproteinaemia. Our findings provide insights into mechanisms driving progressive NAFLD.

Impact and Implications

A rare genetic variant in the gene MTTP has been identified as responsible for the development of severe non-alcoholic fatty liver disease in a four-generation family with no typical disease risk factors. A cell line culture created harbouring this variant gene was characterised to understand how this genetic variation leads to a defect in liver cells, which results in accumulation of fat and processes that promote disease. This is now a useful model for studying the disease pathways and to discover new ways to treat common types of fatty liver disease.

Association study of human leukocyte antigen (HLA) variants and idiopathic pulmonary fibrosis

Introduction

Idiopathic pulmonary fibrosis (IPF) is a chronic interstitial pneumonia marked by progressive lung fibrosis and a poor prognosis. Recent studies have highlighted the potential role of infection in the pathogenesis of IPF and a prior association of the HLA-DQB1 gene with idiopathic fibrotic interstitial pneumonia (including IPF) has been reported. Due to the important role that the Human Leukocyte Antigen (HLA) region plays in the immune response, here we evaluated if HLA genetic variation was associated specifically with IPF risk.

Methods

We performed a meta-analysis of associations of the HLA region with IPF risk in individuals of European ancestry from seven independent case-control studies of IPF (comprising a total of 5,159 cases and 27,459 controls, including the prior study of fibrotic interstitial pneumonia). Single nucleotide polymorphisms, classical HLA alleles and amino acids were analysed and signals meeting a region-wide association threshold p<4.5×10-4 and a posterior probability of replication >90% were considered significant. We sought to replicate the previously reported HLA-DQB1 association in the subset of studies independent of the original report.

Results

The meta-analysis of all seven studies identified four significant independent single nucleotide polymorphisms associated with IPF risk. However, none met the posterior probability for replication criterion. The HLA-DQB1 association was not replicated in the independent IPF studies.

Conclusion

Variation in the HLA region was not consistently associated with risk in studies of IPF. However, this does not preclude the possibility that other genomic regions linked to the immune response may be involved in the aetiology of IPF.

DeepPheWAS: an R package for phenotype generation and association analysis for phenome-wide association studies

SUMMARY

DeepPheWAS is an R package for phenome wide association studies that creates clinically-curated composite phenotypes, and integrates quantitative phenotypes from primary care data, longitudinal trajectories of quantitative measures, disease progression, and drug response phenotypes. Tools are provided for efficient analysis of association with any genetic input, under any genetic model, with optional sex-stratified analysis, and for developing novel phenotypes.

AVAILABILITY AND IMPLEMENTATION

The DeepPheWAS R package is freely available under GNU general public licence v3.0 from at https://github.com/Richard-Packer/DeepPheWAS.

SUPPLEMENTARY INFORMATION

Supplementary methods and results are available at Bioinformatics online.

Exome-wide analysis of copy number variation shows association of the human leukocyte antigen region with asthma in UK Biobank

BACKGROUND

The role of copy number variants (CNVs) in susceptibility to asthma is not well understood. This is, in part, due to the difficulty of accurately measuring CNVs in large enough sample sizes to detect associations. The recent availability of whole-exome sequencing (WES) in large biobank studies provides an unprecedented opportunity to study the role of CNVs in asthma.

METHODS

We called common CNVs in 49,953 individuals in the first release of UK Biobank WES using ClinCNV software. CNVs were tested for association with asthma in a stage 1 analysis comprising 7098 asthma cases and 36,578 controls from the first release of sequencing data. Nominally-associated CNVs were then meta-analysed in stage 2 with an additional 17,280 asthma cases and 115,562 controls from the second release of UK Biobank exome sequencing, followed by validation and fine-mapping.

RESULTS

Five of 189 CNVs were associated with asthma in stage 2, including a deletion overlapping the HLA-DQA1 and HLA-DQB1 genes, a duplication of CHROMR/PRKRA, deletions within MUC22 and TAP2, and a duplication in FBRSL1. The HLA-DQA1, HLA-DQB1, MUC22 and TAP2 genes all reside within the human leukocyte antigen (HLA) region on chromosome 6. In silico analyses demonstrated that the deletion overlapping HLA-DQA1 and HLA-DQB1 is likely to be an artefact arising from under-mapping of reads from non-reference HLA haplotypes, and that the CHROMR/PRKRA and FBRSL1 duplications represent presence/absence of pseudogenes within the HLA region. Bayesian fine-mapping of the HLA region suggested that there are two independent asthma association signals. The variants with the largest posterior inclusion probability in the two credible sets were an amino acid change in HLA-DQB1 (glutamine to histidine at residue 253) and a multi-allelic amino acid change in HLA-DRB1 (presence/absence of serine, glycine or leucine at residue 11).

CONCLUSIONS

At least two independent loci characterised by amino acid changes in the HLA-DQA1, HLA-DQB1 and HLA-DRB1 genes are likely to account for association of SNPs and CNVs in this region with asthma. The high divergence of haplotypes in the HLA can give rise to spurious CNVs, providing an important, cautionary tale for future large-scale analyses of sequencing data.

Balancing selection at the human salivary agglutinin gene (DMBT1) driven by host-microbe interactions

Discovering loci under balancing selection in humans can identify loci with alleles that affect response to the environment and disease. Genome variation data have identified the 5′ region of the DMBT1 gene as undergoing balancing selection in humans. DMBT1 encodes the pattern-recognition glycoprotein DMBT1, also known as SALSA, gp340, or salivary agglutinin. DMBT1 binds to a variety of pathogens through a tandemly arranged scavenger receptor cysteine-rich (SRCR) domain, with the number of domains polymorphic in humans. We show that the signal of balancing selection is driven by one haplotype usually carrying a shorter SRCR repeat and another usually carrying a longer SRCR repeat. DMBT1 encoded by a shorter SRCR repeat allele does not bind a cariogenic and invasive Streptococcus mutans strain, in contrast to the long SRCR allele that shows binding. Our results suggest that balancing selection at DMBT1 is due to host-microbe interactions of encoded SRCR tandem repeat alleles.

•

Clear evidence from many analyses show balancing selection at DMBT1

•

Scavenger-receptor cysteine-rich domain array associated with balancing selection

•

Genetic variation, not alternative splicing, responsible for protein isoforms

•

Long, but not short, DMBT1 isoforms bind a cariogenic strain of Streptococcus mutans

Human intelectin-2 (ITLN2) is selectively expressed by secretory Paneth cells

Intelectins (intestinal lectins) are highly conserved across chordate evolution and have been implicated in various human diseases, including Crohn's disease (CD). The human genome encodes two intelectin genes, intelectin-1 (ITLN1) and intelectin-2 (ITLN2). Other than its high sequence similarity with ITLN1, little is known about ITLN2. To address this void in knowledge, we report that ITLN2 exhibits discrete, yet notable differences from ITLN1 in primary structure, including a unique amino terminus, as well as changes in amino acid residues associated with the glycan-binding activity of ITLN1. We identified that ITLN2 is a highly abundant Paneth cell-specific product, which localizes to secretory granules, and is expressed as a multimeric protein in the small intestine. In surgical specimens of ileal CD, ITLN2 mRNA levels were reduced approximately five-fold compared to control specimens. The ileal expression of ITLN2 was unaffected by previously reported disease-associated variants in ITLN2 and CD-associated variants in neighboring ITLN1 as well as NOD2 and ATG16L1. ITLN2 mRNA expression was undetectable in control colon tissue; however, in both ulcerative colitis (UC) and colonic CD, metaplastic Paneth cells were found to express ITLN2. Together, the data reported establish the groundwork for understanding ITLN2 function(s) in the intestine, including its possible role in CD.

Extended Cohort for E-health, Environment and DNA (EXCEED) COVID-19 focus

Background: New data collection in established longitudinal population studies provides an opportunity for studying the risk factors and sequelae of the novel coronavirus disease 2019 (COVID-19), plus the indirect impacts of the COVID-19 pandemic on wellbeing. The Extended Cohort for E-health, Environment and DNA (EXCEED) cohort is a population-based cohort (N>11,000), recruited from 2013 in Leicester, Leicestershire and Rutland. EXCEED includes consent for electronic healthcare record (EHR) linkage, spirometry, genomic data, and questionnaire data. Methods: Between May 2020 and July 2021, a new questionnaire was deployed in EXCEED, which captured COVID-19 symptoms, general physical and mental health, plus socioeconomic and environmental factors during the pandemic. An online system was developed to invite new participants to join EXCEED, with informed consent being provided online. New and existing participants then completed the COVID-19 questionnaire online. A subset of the new questionnaire respondents were invited to participate in COVID-19 serology substudies, using home antibody testing kits. Results: In total, 3,693 participants provided COVID-19 infection status (median age 62.9 (IQR 54.7-69.2), 58.9% female). Trends of monthly incidence proportions of COVID-19 in EXCEED (self-report or symptom-predicted) approximated local and national figures. Regression analysis of 2,768 participants with linked EHR data showed no obvious monotonic relationship between number of chronic diseases (of 16 pre-specified diseases) and COVID-19 infection. There were 2,144 participants with valid results from a kit allowing differentiation between antibodies due to vaccination or infection. Of these, 8.5% had results consistent with previous COVID-19 infection, and 85.9% had evidence of COVID-19 vaccination, but without evidence of infection. Conclusions: Enriching EXCEED with a new COVID-19 questionnaire and serology data may improve understanding of the risk factors, clinical sequelae and broader impacts of the COVID-19 pandemic in the general population. Controlled access to these data for bona fide researchers is via application to the EXCEED study.

Genome structural variation in human evolution

Structural variation (SV) is a large difference (typically >100 bp) in the genomic structure of two genomes and includes both copy number variation and variation that does not change copy number of a genomic region, such as an inversion. Improved reference genomes, combined with widespread genome sequencing using short-read sequencing technology, and increasingly using long-read sequencing, have reignited interest in SV. Recent large-scale studies and functional focused analyses have highlighted the role of SV in human evolution. In this review, we highlight human-specific SVs involved in changes in the brain, population-specific SVs that affect response to the environment, including adaptation to diet and infectious diseases, and summarise the contribution of archaic hominin admixture to present-day human SV.

Human intelectin-1 (ITLN1) genetic variation and intestinal expression

Intelectins are ancient carbohydrate binding proteins, spanning chordate evolution and implicated in multiple human diseases. Previous GWAS have linked SNPs in ITLN1 (also known as omentin) with susceptibility to Crohn's disease (CD); however, analysis of possible functional significance of SNPs at this locus is lacking. Using the Ensembl database, pairwise linkage disequilibrium (LD) analyses indicated that several disease-associated SNPs at the ITLN1 locus, including SNPs in CD244 and Ly9, were in LD. The alleles comprising the risk haplotype are the major alleles in European (67%), but minor alleles in African superpopulations. Neither ITLN1 mRNA nor protein abundance in intestinal tissue, which we confirm as goblet-cell derived, was altered in the CD samples overall nor when samples were analyzed according to genotype. Moreover, the missense variant V109D does not influence ITLN1 glycan binding to the glycan β-D-galactofuranose or protein-protein oligomerization. Taken together, our data are an important step in defining the role(s) of the CD-risk haplotype by determining that risk is unlikely to be due to changes in ITLN1 carbohydrate recognition, protein oligomerization, or expression levels in intestinal mucosa. Our findings suggest that the relationship between the genomic data and disease arises from changes in CD244 or Ly9 biology, differences in ITLN1 expression in other tissues, or an alteration in ITLN1 interaction with other proteins.

Clonal architecture in mesothelioma is prognostic and shapes the tumour microenvironment

Malignant Pleural Mesothelioma (MPM) is typically diagnosed 20-50 years after exposure to asbestos and evolves along an unknown evolutionary trajectory. To elucidate this path, we conducted multi-regional exome sequencing of 90 tumour samples from 22 MPMs acquired at surgery. Here we show that exomic intratumour heterogeneity varies widely across the cohort. Phylogenetic tree topology ranges from linear to highly branched, reflecting a steep gradient of genomic instability. Using transfer learning, we detect repeated evolution, resolving 5 clusters that are prognostic, with temporally ordered clonal drivers. BAP1/-3p21 and FBXW7/-chr4 events are always early clonal. In contrast, NF2/-22q events, leading to Hippo pathway inactivation are predominantly late clonal, positively selected, and when subclonal, exhibit parallel evolution indicating an evolutionary constraint. Very late somatic alteration of NF2/22q occurred in one patient 12 years after surgery. Clonal architecture and evolutionary clusters dictate MPM inflammation and immune evasion. These results reveal potentially drugable evolutionary bottlenecking in MPM, and an impact of clonal architecture on shaping the immune landscape, with potential to dictate the clinical response to immune checkpoint inhibition.

Genotyping complex structural variation at the malaria-associated human glycophorin locus using a PCR-based strategy

Structural variation in the human genome can affect risk of disease. An example is a complex structural variant of the human glycophorin gene cluster, called DUP4, which is associated with a clinically significant level of protection against severe malaria. The human glycophorin gene cluster harbours at least 23 distinct structural variants, and accurate genotyping of this complex structural variation remains a challenge. Here, we use a polymerase chain reaction-based strategy to genotype structural variation at the human glycophorin gene cluster, including the alleles responsible for the U- blood group. We validate our approach, based on a triplex paralogue ratio test, on publically available samples from the 1000 Genomes project. We then genotype 574 individuals from a longitudinal birth cohort (Tori-Bossito cohort) using small amounts of DNA at low cost. Our approach readily identifies known deletions and duplications, and can potentially identify novel variants for further analysis. It will allow exploration of genetic variation at the glycophorin locus, and investigation of its relationship with malaria, in large sample sets at minimal cost, using standard molecular biology equipment.

Structural variation of the malaria-associated human glycophorin A-B-E region

BACKGROUND

Approximately 5% of the human genome shows common structural variation, which is enriched for genes involved in the immune response and cell-cell interactions. A well-established region of extensive structural variation is the glycophorin gene cluster, comprising three tandemly-repeated regions about 120 kb in length and carrying the highly homologous genes GYPA, GYPB and GYPE. Glycophorin A (encoded by GYPA) and glycophorin B (encoded by GYPB) are glycoproteins present at high levels on the surface of erythrocytes, and they have been suggested to act as decoy receptors for viral pathogens. They are receptors for the invasion of the protist parasite Plasmodium falciparum, a causative agent of malaria. A particular complex structural variant, called DUP4, creates a GYPB-GYPA fusion gene known to confer resistance to malaria. Many other structural variants exist across the glycophorin gene cluster, and they remain poorly characterised.

RESULTS

Here, we analyse sequences from 3234 diploid genomes from across the world for structural variation at the glycophorin locus, confirming 15 variants in the 1000 Genomes project cohort, discovering 9 new variants, and characterising a selection of these variants using fibre-FISH and breakpoint mapping at the sequence level. We identify variants predicted to create novel fusion genes and a common inversion duplication variant at appreciable frequencies in West Africans. We show that almost all variants can be explained by non-allelic homologous recombination and by comparing the structural variant breakpoints with recombination hotspot maps, confirm the importance of a particular meiotic recombination hotspot on structural variant formation in this region.

CONCLUSIONS

We identify and validate large structural variants in the human glycophorin A-B-E gene cluster which may be associated with different clinical aspects of malaria.

The Malaria-Protective Human Glycophorin Structural Variant DUP4 Shows Somatic Mosaicism and Association with Hemoglobin Levels

Glycophorin A and glycophorin B are red blood cell surface proteins and are both receptors for the parasite Plasmodium falciparum, which is the principal cause of malaria in sub-Saharan Africa. DUP4 is a complex structural genomic variant that carries extra copies of a glycophorin A-glycophorin B fusion gene and has a dramatic effect on malaria risk by reducing the risk of severe malaria by up to 40%. Using fiber-FISH and Illumina sequencing, we validate the structural arrangement of the glycophorin locus in the DUP4 variant and reveal somatic variation in copy number of the glycophorin B-glycophorin A fusion gene. By developing a simple, specific, PCR-based assay for DUP4, we show that the DUP4 variant reaches a frequency of 13% in the population of a malaria-endemic village in south-eastern Tanzania. We genotype a substantial proportion of that village and demonstrate an association of DUP4 genotype with hemoglobin levels, a phenotype related to malaria, using a family-based association test. Taken together, we show that DUP4 is a complex structural variant that may be susceptible to somatic variation and show that DUP4 is associated with a malarial-related phenotype in a longitudinally followed population.

Human CCL3L1 copy number variation, gene expression, and the role of the CCL3L1-CCR5 axis in lung function

Background: The CCL3L1-CCR5 signaling axis is important in a number of inflammatory responses, including macrophage function, and T-cell-dependent immune responses. Small molecule CCR5 antagonists exist, including the approved antiretroviral drug maraviroc, and therapeutic monoclonal antibodies are in development. Repositioning of drugs and targets into new disease areas can accelerate the availability of new therapies and substantially reduce costs. As it has been shown that drug targets with genetic evidence supporting their involvement in the disease are more likely to be successful in clinical development, using genetic association studies to identify new target repurposing opportunities could be fruitful. Here we investigate the potential of perturbation of the CCL3L1-CCR5 axis as treatment for respiratory disease. Europeans typically carry between 0 and 5 copies of CCL3L1 and this multi-allelic variation is not detected by widely used genome-wide single nucleotide polymorphism studies.  Methods: We directly measured the complex structural variation of CCL3L1 using the Paralogue Ratio Test and imputed (with validation) CCR5del32 genotypes in 5,000 individuals from UK Biobank, selected from the extremes of the lung function distribution, and analysed DNA and RNAseq data for CCL3L1 from the 1000 Genomes Project. Results: We confirmed the gene dosage effect of CCL3L1 copy number on CCL3L1 mRNA expression levels.  We found no evidence for association of CCL3L1 copy number or CCR5del32 genotype with lung function. Conclusions: These results suggest that repositioning CCR5 antagonists is unlikely to be successful for the treatment of airflow obstruction.

Human CCL3L1 copy number variation, gene expression, and the role of the CCL3L1-CCR5 axis in lung function

Background: The CCL3L1-CCR5 signaling axis is important in a number of inflammatory responses, including macrophage function, and T-cell-dependent immune responses. Small molecule CCR5 antagonists exist, including the approved antiretroviral drug maraviroc, and therapeutic monoclonal antibodies are in development. Repositioning of drugs and targets into new disease areas can accelerate the availability of new therapies and substantially reduce costs. As it has been shown that drug targets with genetic evidence supporting their involvement in the disease are more likely to be successful in clinical development, using genetic association studies to identify new target repurposing opportunities could be fruitful. Here we investigate the potential of perturbation of the CCL3L1-CCR5 axis as treatment for respiratory disease. Europeans typically carry between 0 and 5 copies of CCL3L1 and this multi-allelic variation is not detected by widely used genome-wide single nucleotide polymorphism studies.  Methods: We directly measured the complex structural variation of CCL3L1 using the Paralogue Ratio Test and imputed (with validation) CCR5del32 genotypes in 5,000 individuals from UK Biobank, selected from the extremes of the lung function distribution, and analysed DNA and RNAseq data for CCL3L1 from the 1000 Genomes Project. Results: We confirmed the gene dosage effect of CCL3L1 copy number on CCL3L1 mRNA expression levels.  We found no evidence for association of CCL3L1 copy number or CCR5del32 genotype with lung function. Conclusions: These results suggest that repositioning CCR5 antagonists is unlikely to be successful for the treatment of airflow obstruction.

Population genetics of immune-related multilocus copy number variation in Native Americans

While multiallelic copy number variation (mCNV) loci are a major component of genomic variation, quantifying the individual copy number of a locus and defining genotypes is challenging. Few methods exist to study how mCNV genetic diversity is apportioned within and between populations (i.e. to define the population genetic structure of mCNV). These inferences are critical in populations with a small effective size, such as Amerindians, that may not fit the Hardy-Weinberg model due to inbreeding, assortative mating, population subdivision, natural selection or a combination of these evolutionary factors. We propose a likelihood-based method that simultaneously infers mCNV allele frequencies and the population structure parameter f, which quantifies the departure of homozygosity from the Hardy-Weinberg expectation. This method is implemented in the freely available software CNVice, which also infers individual genotypes using information from both the population and from trios, if available. We studied the population genetics of five immune-related mCNV loci associated with complex diseases (beta-defensins, CCL3L1/CCL4L1, FCGR3A, FCGR3B and FCGR2C) in 12 traditional Native American populations and found that the population structure parameters inferred for these mCNVs are comparable to but lower than those for single nucleotide polymorphisms studied in the same populations.

No Evidence for Association of β-Defensin Genomic Copy Number with HIV Susceptibility, HIV Load during Clinical Latency, or Progression to AIDS

Common single-nucleotide variation in the host accounts for 25% of the variability in the plasma levels of HIV during the clinical latency stage (viral load set point). However, the role of rare variants and copy number variants remains relatively unexplored. Previous work has suggested copy number variation of a cluster of β-defensin genes affects HIV load in treatment-naïve sub-Saharan Africans and rate of response to antiretroviral treatment. Here we analyse a total of 1827 individuals from two cohorts of HIV-infected individuals from Europe and sub-Saharan Africa to investigate the role of β-defensin copy number variation on HIV load at set point. We find no evidence for association of copy number with viral load. We also compare distribution of β-defensin copy number between European cases and controls and find no differences, arguing against a role of β-defensin copy number in HIV acquisition. Taken together, our data argue against an effect of copy number variation of the β-defensin region in the spontaneous control of HIV infection.

Understanding the Genomic Structure of Copy-Number Variation of the Low-Affinity Fcγ Receptor Region Allows Confirmation of the Association of FCGR3B Deletion with Rheumatoid Arthritis

Fcγ receptors are a family of cell–surface receptors that are expressed by a host of different innate and adaptive immune cells, and mediate inflammatory responses by binding the Fc portion of immunoglobulin G. In humans, five low‐affinity receptors are encoded by the genes FCGR2A, FCGR2B, FCGR2C, FCGR3A, and FCGR3B, which are located in an 82.5‐kb segmental tandem duplication on chromosome 1q23.3, which shows extensive copy‐number variation (CNV). Deletions of FCGR3B have been suggested to increase the risk of inflammatory diseases such as systemic lupus erythematosus and rheumatoid arthritis (RA). In this study, we identify the deletion breakpoints of FCGR3B deletion alleles in the UK population and endogamous native American population, and show that some but not all alleles are likely to be identical‐by‐descent. We also localize a duplication breakpoint, confirming that the mechanism of CNV generation is nonallelic homologous recombination, and identify several alleles with gene conversion events using fosmid sequencing data. We use information on the structure of the deletion alleles to distinguish FCGR3B deletions from FCGR3A deletions in whole‐genome array comparative genomic hybridization (aCGH) data. Reanalysis of published aCGH data using this approach supports association of FCGR3B deletion with increased risk of RA in a large cohort of 1,982 cases and 3,271 controls (odds ratio 1.61, P = 2.9×10⁻³).

Analysis of copy number variation at DMBT1 and age-related macular degeneration

BACKGROUND

DMBT1 is a gene that shows extensive copy number variation (CNV) that alters the number of bacteria-binding domains in the protein and has been shown to activate the complement pathway. It lies next to the ARMS2/HTRA1 genes in a region of chromosome 10q26, where single nucleotide variants have been strongly associated with age-related macular degeneration (AMD), the commonest cause of blindness in Western populations. Complement activation is thought to be a key factor in the pathogenesis of this condition. We sought to investigate whether DMBT1 CNV plays any role in the susceptibility to AMD.

METHODS

We analysed long-range linkage disequilibrium of DMBT1 CNV1 and CNV2 with flanking single nucleotide polymorphisms (SNPs) using our previously published CNV and HapMap Phase 3 SNP data in the CEPH Europeans from Utah (CEU). We then typed a large cohort of 860 AMD patients and 419 examined age-matched controls for copy number at DMBT1 CNV1 and CNV2 and combined these data with copy numbers from a further 480 unexamined controls.

RESULTS

We found weak linkage disequilibrium between DMBT1 CNV1 and CNV2 with the SNPs rs1474526 and rs714816 in the HTRA1/ARMS2 region. By directly analysing copy number variation, we found no evidence of association of CNV1 or CNV2 with AMD.

CONCLUSIONS

We have shown that copy number variation at DMBT1 does not affect risk of developing age-related macular degeneration and can therefore be ruled out from future studies investigating the association of structural variation at 10q26 with AMD.

Evaluation of High-Throughput Genomic Assays for the Fc Gamma Receptor Locus

Cancer immunotherapy has been revolutionised by the use monoclonal antibodies (mAb) that function through their interaction with Fc gamma receptors (FcγRs). The low-affinity FcγR genes are highly homologous, map to a complex locus at 1p23 and harbour single nucleotide polymorphisms (SNPs) and copy number variation (CNV) that can impact on receptor function and response to therapeutic mAbs. This complexity can hinder accurate characterisation of the locus. We therefore evaluated and optimised a suite of assays for the genomic analysis of the FcγR locus amenable to peripheral blood mononuclear cells and formalin-fixed paraffin-embedded (FFPE) material that can be employed in a high-throughput manner. Assessment of TaqMan genotyping for FCGR2A-131H/R, FCGR3A-158F/V and FCGR2B-232I/T SNPs demonstrated the need for additional methods to discriminate genotypes for the FCGR3A-158F/V and FCGR2B-232I/T SNPs due to sequence homology and CNV in the region. A multiplex ligation-dependent probe amplification assay provided high quality SNP and CNV data in PBMC cases, but there was greater data variability in FFPE material in a manner that was predicted by the BIOMED-2 multiplex PCR protocol. In conclusion, we have evaluated a suite of assays for the genomic analysis of the FcγR locus that are scalable for application in large clinical trials of mAb therapy. These assays will ultimately help establish the importance of FcγR genetics in predicting response to antibody therapeutics.

Determining multiallelic complex copy number and sequence variation from high coverage exome sequencing data

BACKGROUND

Copy number variation (CNV) is a major component of genomic variation, yet methods to accurately type genomic CNV lag behind methods that type single nucleotide variation. High-throughput sequencing can contribute to these methods by using sequence read depth, which takes the number of reads that map to a given part of the reference genome as a proxy for copy number of that region, and compares across samples. Furthermore, high-throughput sequencing also provides information on the sequence differences between copies within and between individuals.

METHODS

In this study we use high-coverage phase 3 exome sequences of the 1000 Genomes project to infer diploid copy number of the beta-defensin genomic region, a well-studied CNV that carries several beta-defensin genes involved in the antimicrobial response, signalling, and fertility. We also use these data to call sequence variants, a particular challenge given the multicopy nature of the region.

RESULTS

We confidently call copy number and sequence variation of the beta-defensin genes on 1285 samples from 26 global populations, validate copy number using Nanostring nCounter and triplex paralogue ratio test data. We use the copy number calls to verify the genomic extent of the CNV and validate sequence calls using analysis of cloned PCR products. We identify novel variation, mostly individually rare, predicted to alter amino-acid sequence in the beta-defensin genes. Such novel variants may alter antimicrobial properties or have off-target receptor interactions, and may contribute to individuality in immunological response and fertility.

CONCLUSIONS

Given that 81% of identified sequence variants were not previously in dbSNP, we show that sequence variation in multiallelic CNVs represent an unappreciated source of genomic diversity.

A comparison of assays for accurate copy number measurement of the low-affinity Fc gamma receptor genes FCGR3A and FCGR3B

The FCGR3 locus encoding the low affinity activating receptor FcγRIII, plays a vital role in immunity triggered by cellular effector and regulatory functions. Copy number of the genes FCGR3A and FCGR3B has previously been reported to affect susceptibility to several autoimmune diseases and chronic inflammatory conditions. However, such genetic association studies often yield inconsistent results; hence require assays that are robust with low error rate. We investigated the accuracy and efficiency in estimating FCGR3 CNV by comparing Sequenom MassARRAY and paralogue ratio test-restriction enzyme digest variant ratio (PRT-REDVR). In addition, since many genetic association studies of FCGR3B CNV were carried out using real-time quantitative PCR, we have also included the evaluation of that method’s performance in estimating the multi-allelic CNV of FCGR3B. The qPCR assay exhibited a considerably broader distribution of signal intensity, potentially introducing error in estimation of copy number and higher false positive rates. Both Sequenom and PRT-REDVR showed lesser systematic bias, but Sequenom skewed towards copy number normal (CN = 2). The discrepancy between Sequenom and PRT-REDVR might be attributed either to batch effects noise in individual measurements. Our study suggests that PRT-REDVR is more robust and accurate in genotyping the CNV of FCGR3, but highlights the needs of multiple independent assays for extensive validation when performing a genetic association study with multi-allelic CNVs.

Evidence of convergent evolution in humans and macaques supports an adaptive role for copy number variation of the β-defensin-2 gene

β-defensins are a family of important peptides of innate immunity, involved in host defense, immunomodulation, reproduction, and pigmentation. Genes encoding β-defensins show evidence of birth-and-death evolution, adaptation by amino acid sequence changes, and extensive copy number variation (CNV) within humans and other species. The role of CNV in the adaptation of β-defensins to new functions remains unclear, as does the adaptive role of CNV in general. Here, we fine-map CNV of a cluster of β-defensins in humans and rhesus macaques. Remarkably, we found that the structure of the CNV is different between primates, with distinct mutational origins and CNV boundaries defined by retroviral long terminal repeat elements. Although the human β-defensin CNV region is 322 kb and encompasses several genes, including β-defensins, a long noncoding RNA gene, and testes-specific zinc-finger transcription factors, the orthologous region in the rhesus macaque shows CNV of a 20-kb region, containing only a single gene, the ortholog of the human β-defensin-2 gene. Despite its independent origins, the range of gene copy numbers in the rhesus macaque is similar to humans. In addition, the rhesus macaque gene has been subject to divergent positive selection at the amino acid level following its initial duplication event between 3 and 9.5 Ma, suggesting adaptation of this gene as the macaque successfully colonized novel environments outside Africa. Therefore, the molecular phenotype of β-defensin-2 CNV has undergone convergent evolution, and this gene shows evidence of adaptation at the amino acid level in rhesus macaques.

Copy number variation of the beta defensin gene cluster on chromosome 8p influences the bacterial microbiota within the nasopharynx of otitis-prone children

As there is increasing evidence that aberrant defensin expression is related to susceptibility for infectious disease and inflammatory disorders, we sought to determine if copy number of the beta-defensin gene cluster located on chromosome 8p23.1 (DEFB107, 106, 105, 104, 103, DEFB4 and SPAG11), that shows copy number variation as a block, was associated with susceptibility to otitis media (OM). The gene DEFB103 within this complex encodes human beta defensin-3 (hBD-3), an antimicrobial peptide (AP) expressed by epithelial cells that line the mammalian airway, important for defense of mucosal surfaces and previously shown to have bactericidal activity in vitro against multiple human pathogens, including the three that predominate in OM. To this end, we conducted a retrospective case-control study of 113 OM prone children and 267 controls aged five to sixty months. We identified the copy number of the above defined beta-defensin gene cluster (DEFB-CN) in each study subject by paralogue ratio assays. The mean DEFB-CN was indistinguishable between subjects classified as OM prone based on a recent history of multiple episodes of OM and control subjects who had no history of OM (4.4 ± 0.96 versus 4.4 ± 1.08, respectively: Odds Ratio [OR]: 1.16 (95% CI: 0.61, 2.20). Despite a lack of direct association, we observed a statistically significant correlation between DEFB-CN and nasopharyngeal bacterial colonization patterns. Collectively, our findings suggested that susceptibility to OM might be mediated by genetic variation among individuals, wherein a DEFB-CN less than 4 exerts a marked influence on the microbiota of the nasopharynx, specifically with regard to colonization by the three predominant bacterial pathogens of OM.

Whole exome re-sequencing implicates CCDC38 and cilia structure and function in resistance to smoking related airflow obstruction

Chronic obstructive pulmonary disease (COPD) is a leading cause of global morbidity and mortality and, whilst smoking remains the single most important risk factor, COPD risk is heritable. Of 26 independent genomic regions showing association with lung function in genome-wide association studies, eleven have been reported to show association with airflow obstruction. Although the main risk factor for COPD is smoking, some individuals are observed to have a high forced expired volume in 1 second (FEV1) despite many years of heavy smoking. We hypothesised that these "resistant smokers" may harbour variants which protect against lung function decline caused by smoking and provide insight into the genetic determinants of lung health. We undertook whole exome re-sequencing of 100 heavy smokers who had healthy lung function given their age, sex, height and smoking history and applied three complementary approaches to explore the genetic architecture of smoking resistance. Firstly, we identified novel functional variants in the "resistant smokers" and looked for enrichment of these novel variants within biological pathways. Secondly, we undertook association testing of all exonic variants individually with two independent control sets. Thirdly, we undertook gene-based association testing of all exonic variants. Our strongest signal of association with smoking resistance for a non-synonymous SNP was for rs10859974 (P = 2.34 × 10(-4)) in CCDC38, a gene which has previously been reported to show association with FEV1/FVC, and we demonstrate moderate expression of CCDC38 in bronchial epithelial cells. We identified an enrichment of novel putatively functional variants in genes related to cilia structure and function in resistant smokers. Ciliary function abnormalities are known to be associated with both smoking and reduced mucociliary clearance in patients with COPD. We suggest that genetic influences on the development or function of cilia in the bronchial epithelium may affect growth of cilia or the extent of damage caused by tobacco smoke.

Copy-number variation of the neuronal glucose transporter gene SLC2A3 and age of onset in Huntington's disease

Huntington's disease (HD) is a devastating neurodegenerative disorder which is inherited in an autosomal dominant manner. HD is caused by a trinucleotide CAG repeat expansion that encodes a polyglutamine stretch in the huntingtin (HTT) protein. Mutant HTT expression leads to a myriad of cellular dysfunctions culminating in neuronal loss and consequent motor, cognitive and psychiatric disturbances in HD patients. The length of the CAG repeat is inversely correlated with age of onset (AO) in HD patients, while environmental and genetic factors can further modulate this parameter. Here, we explored whether the recently described copy-number variation (CNV) of the gene SLC2A3-which encodes the neuronal glucose transporter GLUT3-could modulate AO in HD. Strikingly, we found that increased dosage of SLC2A3 delayed AO in an HD cohort of 987 individuals, and that this correlated with increased levels of GLUT3 in HD patient cells. To our knowledge this is the first time that CNV of a candidate gene has been found to modulate HD pathogenesis. Furthermore, we found that increasing dosage of Glut1-the Drosophila melanogaster homologue of this glucose transporter-ameliorated HD-relevant phenotypes in fruit flies, including neurodegeneration and life expectancy. As alterations in glucose metabolism have been implicated in HD pathogenesis, this study may have important therapeutic relevance for HD.

Copy number variation of the beta-defensin genes in europeans: no supporting evidence for association with lung function, chronic obstructive pulmonary disease or asthma

Lung function measures are heritable, predict mortality and are relevant in diagnosis of chronic obstructive pulmonary disease (COPD). COPD and asthma are diseases of the airways with major public health impacts and each have a heritable component. Genome-wide association studies of SNPs have revealed novel genetic associations with both diseases but only account for a small proportion of the heritability. Complex copy number variation may account for some of the missing heritability. A well-characterised genomic region of complex copy number variation contains beta-defensin genes (DEFB103, DEFB104 and DEFB4), which have a role in the innate immune response. Previous studies have implicated these and related genes as being associated with asthma or COPD. We hypothesised that copy number variation of these genes may play a role in lung function in the general population and in COPD and asthma risk. We undertook copy number typing of this locus in 1149 adult and 689 children using a paralogue ratio test and investigated association with COPD, asthma and lung function. Replication of findings was assessed in a larger independent sample of COPD cases and smoking controls. We found evidence for an association of beta-defensin copy number with COPD in the adult cohort (OR = 1.4, 95%CI:1.02-1.92, P = 0.039) but this finding, and findings from a previous study, were not replicated in a larger follow-up sample(OR = 0.89, 95%CI:0.72-1.07, P = 0.217). No robust evidence of association with asthma in children was observed. We found no evidence for association between beta-defensin copy number and lung function in the general populations. Our findings suggest that previous reports of association of beta-defensin copy number with COPD should be viewed with caution. Suboptimal measurement of copy number can lead to spurious associations. Further beta-defensin copy number measurement in larger sample sizes of COPD cases and children with asthma are needed.

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

BACKGROUND

The 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.

METHODS

We 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.

RESULTS

We 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.

CONCLUSIONS

Our 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.

Copy number variation of Fc gamma receptor genes in HIV-infected and HIV-tuberculosis co-infected individuals in sub-Saharan Africa

AIDS, caused by the retrovirus HIV, remains the largest cause of morbidity in sub-Saharan Africa yet almost all genetic studies have focused on cohorts from Western countries. HIV shows high co-morbidity with tuberculosis (TB), as HIV stimulates the reactivation of latent tuberculosis (TB). Recent clinical trials suggest that an effective anti-HIV response correlates with non-neutralising antibodies. Given that Fcγ receptors are critical in mediating the non-neutralising effects of antibodies, analysis of the extensive variation at Fcγ receptor genes is important. Single nucleotide variation and copy number variation (CNV) of Fcγ receptor genes affects the expression profile, activatory/inhibitory balance, and IgG affinity of the Fcγ receptor repertoire of each individual. In this study we investigated whether CNV of FCGR2C, FCGR3A and FCGR3B as well as the HNA1 allotype of FCGR3B is associated with HIV load, response to highly-active antiretroviral therapy (HAART) and co-infection with TB. We confirmed an effect of TB-co-infection status on HIV load and response to HAART, but no conclusive effect of the genetic variants we tested. We observed a small effect, in Ethiopians, of FCGR3B copy number, where deletion was more frequent in HIV-TB co-infected patients than those infected with HIV alone.

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.

β-Defensin Genomic Copy Number Does Not Influence the Age of Onset in Huntington's Disease

BACKGROUND

Huntington's disease (HD) is an autosomal dominant neurodegenerative disorder caused by the abnormal expansion of a CAG triplet repeat tract in the huntingtin gene. While the length of this CAG expansion is the major determinant of the age of onset (AO), other genetic factors have also been shown to play a modulatory role. Recent evidence suggests that neuroinflammation is a pivotal factor in the pathogenesis of HD, and that targeting this process may have important therapeutic ramifications. The human β-defensin 2 (hBD2) - encoded by DEFB4 - is an antimicrobial peptide that exhibits inducible expression in astrocytes during inflammation and is an important regulator of innate and adaptive immune response. Therefore, DEFB4 may contribute to the neuroinflammatory processes observed in HD.

OBJECTIVE

In this study we tested the hypothesis that copy number variation (CNV) of the β-defensin region, including DEFB4, modifies the AO in HD.

METHODS AND RESULTS

We genotyped β-defensin CNV in 490 HD individuals using the paralogue ratio test and found no association between β-defensin CNV and onset of HD.

CONCLUSIONS

We conclude that it is unlikely that DEFB4 plays a role in HD pathogenesis.

β-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.

The challenges of studying complex and dynamic regions of the human genome

Recent work has emphasised that the human genome is not simple and static, but complex and dynamic. This review focuses on the regions that are particularly hard to dissect and analyse, yet hold clues to how the genome changes during evolution and disease. I begin by summarising recent key advances in the understanding of the variable structure of our genome, and then I discuss a medley of methods that may allow us to analyse this structure in fine detail. In the final part, I describe potential future developments in this field, and make an argument that, just as we routinely genotype single-nucleotide polymorphisms now and will routinely re-sequence genomes in the near future, we should be aiming to physically re-map the individual human genome for each individual we study.

A common mutation in the defensin DEFB126 causes impaired sperm function and subfertility

A glycosylated polypeptide, β-defensin 126 (DEFB126), derived from the epididymis and adsorbed onto the sperm surface, has been implicated in immunoprotection and efficient movement of sperm in mucosal fluids of the female reproductive tract. Here, we report a sequence variant in DEFB126 that has a two-nucleotide deletion in the open reading frame, which generates an abnormal mRNA. The allele frequency of this variant sequence was high in both a European (0.47) and a Chinese (0.45) population cohort. Binding of the Agaricus bisporus lectin to the sperm surface glycocalyx was significantly lower in men with the homozygous variant (del/del) genotype than in those with either a del/wt or a wt/wt genotype, suggesting an altered sperm glycocalyx with fewer O-linked oligosaccharides in del/del men. Moreover, sperm from del/del carriers exhibited an 84% reduction in the rate of penetration of a hyaluronic acid gel, a surrogate for cervical mucus, compared to the other genotypes. This reduction in sperm performance in hyaluronic acid gels was not a result of decreased progressive motility (average curvilinear velocity) or morphological deficits. Nevertheless, DEFB126 genotype and lectin binding were correlated with sperm performance in the penetration assays. In a prospective cohort study of newly married couples who were trying to conceive by natural means, couples were less likely to become pregnant and took longer to achieve a live birth if the male partner was homozygous for the variant sequence. This common sequence variation in DEFB126, and its apparent effect of impaired reproductive function, will allow a better understanding, clinical evaluation, and possibly treatment of human infertility.