Segregation analysis identifies specific alpha-defensin (DEFA1A3) SNP-CNV haplotypes in predisposition to IgA nephropathy.

BACKGROUND: Immunoglobulin A (IgA) nephropathy is a disorder of the immune system affecting kidney function, and genome-wide association studies (GWAS) have defined numerous loci with associated variation, all implicating components of innate or adaptive immunity. Among these, single nucleotide polymorphisms (SNPs) in a region including the multiallelic copy number variation (CNV) of DEFA1A3 are associated with IgA nephropathy in both European and Asian populations. At present, the precise factors underlying the observed associations at DEFA1A3 have not been defined, although the key alleles differ between Asian and European populations, and multiple independent factors may be involved even within a single population. METHODS: In this study, we measured DEFA1A3 copy number in UK family trios with an offspring affected by IgA nephropathy, used the population distributions of joint SNP-CNV haplotypes to infer the likely segregation in trios, and applied transmission disequilibrium tests (TDT) to examine joint SNP-CNV haplotypes for over- or undertransmission into affected offspring from heterozygous parents. RESULTS AND CONCLUSIONS: We observed overtransmission of 3-copy class 2 haplotypes (raw p = 0.029) and some evidence for under-transmission of 3-copy class 1 haplotypes (raw p = 0.051), although these apparent effects were not statistically significant after correction for testing of multiple haplotypes.

Skin microbiome alters attractiveness to Anopheles mosquitoes.

BACKGROUND: Some people produce specific body odours that make them more attractive than others to mosquitoes, and consequently are at higher risk of contracting vector-borne diseases. The skin microbiome can break down carbohydrates, fatty acids and peptides on the skin into volatiles that mosquitoes can differentiate. RESULTS: Here, we examined how skin microbiome composition of women differs in relation to level of attractiveness to Anopheles coluzzii mosquitoes, to identify volatiles in body odour and metabolic pathways associated with individuals that tend to be poorly-attractive to mosquitoes. We used behavioural assays to measure attractiveness of participants to An. coluzzii mosquitoes, 16S rRNA amplicon sequencing of the bacteria sampled from the skin and gas chromatography of volatiles in body odour. We found differences in skin microbiome composition between the poorly- and highly-attractive groups, particularly eight Amplicon Sequence Variants (ASVs) belonging to the Proteobacteria, Actinobacteria and Firmicutes phyla. Staphylococcus 2 ASVs are four times as abundant in the highly-attractive compared to poorly-attractive group. Associations were found between these ASVs and volatiles known to be attractive to Anopheles mosquitoes. Propanoic pathways are enriched in the poorly-attractive participants compared to those found to be highly-attractive. CONCLUSIONS: Our findings suggest that variation in attractiveness of people to mosquitoes is related to the composition of the skin microbiota, knowledge that could improve odour-baited traps or other next generation vector control tools.

No Evidence for Association of BMI with Salivary Amylase Gene Copy Number in the UK 1958 Birth Cohort.

OBJECTIVE: In a 2014 publication, evidence was presented supporting the association of BMI with the copy number of the salivary amylase 1 (AMY1) gene, with an unprecedented effect size of -0.15 kg/m2 (SE 0.02) per copy of AMY1. Most well-powered attempts to reproduce these findings have not been successful. However, because of different study designs, a significant association may still apply under restricted conditions such as in particular age groups. This study specifically tested the BMI-AMY1 association at different age points in the same individuals using longitudinal BMI information from participants in the UK 1958 Birth Cohort study. METHODS: This study measured the AMY1 copy number by paralogue ratio tests in genomic DNA and by using array comparative genomic hybridization data. BMI data from 1958 Birth Cohort participants were available from eight different age points between 7 and 50 years. RESULTS: No evidence, even at nominal significance, was found for association of the AMY1 copy number with BMI at any age point in approximately 1,400 members of the 1958 Birth Cohort or in 2,835 people from two disease cohorts from the Wellcome Trust Case Control Consortium. CONCLUSIONS: The results do not support an association between BMI and AMY1 copy number at any age point.

Copy number variation of human AMY1 is a minor contributor to variation in salivary amylase expression and activity.

BACKGROUND: Salivary amylase in humans is encoded by the copy variable gene AMY1 in the amylase gene cluster on chromosome 1. Although the role of salivary amylase is well established, the consequences of the copy number variation (CNV) at AMY1 on salivary amylase protein production are less well understood. The amylase gene cluster is highly structured with a fundamental difference between odd and even AMY1 copy number haplotypes. In this study, we aimed to explore, in samples from 119 unrelated individuals, not only the effects of AMY1 CNV on salivary amylase protein expression and amylase enzyme activity but also whether there is any evidence for underlying difference between the common haplotypes containing odd numbers of AMY1 and even copy number haplotypes. RESULTS: AMY1 copy number was significantly correlated with the variation observed in salivary amylase production (11.7% of variance, P < 0.0005) and enzyme activity (13.6% of variance, P < 0.0005) but did not explain the majority of observed variation between individuals. AMY1-odd and AMY1-even haplotypes showed a different relationship between copy number and expression levels, but the difference was not statistically significant (P = 0.052). CONCLUSIONS: Production of salivary amylase is correlated with AMY1 CNV, but the majority of interindividual variation comes from other sources. Long-range haplotype structure may affect expression, but this was not significant in our data.

Recurrent Rearrangements of Human Amylase Genes Create Multiple Independent CNV Series.

The human amylase gene cluster includes the human salivary (AMY1) and pancreatic amylase genes (AMY2A and AMY2B), and is a highly variable and dynamic region of the genome. Copy number variation (CNV) of AMY1 has been implicated in human dietary adaptation, and in population association with obesity, but neither of these findings has been independently replicated. Despite these functional implications, the structural genomic basis of CNV has only been defined in detail very recently. In this work, we use high-resolution analysis of copy number, and analysis of segregation in trios, to define new, independent allelic series of amylase CNVs in sub-Saharan Africans, including a series of higher-order expansions of a unit consisting of one copy each of AMY1, AMY2A, and AMY2B. We use fiber-FISH (fluorescence in situ hybridization) to define unexpected complexity in the accompanying rearrangements. These findings demonstrate recurrent involvement of the amylase gene region in genomic instability, involving at least five independent rearrangements of the pancreatic amylase genes (AMY2A and AMY2B). Structural features shared by fundamentally distinct lineages strongly suggest that the common ancestral state for the human amylase cluster contained more than one, and probably three, copies of AMY1.

Analysis of Multiallelic CNVs by Emulsion Haplotype Fusion PCR.

Emulsion-fusion PCR recovers long-range sequence information by combining products in cis from individual genomic DNA molecules. Emulsion droplets act as very numerous small reaction chambers in which different PCR products from a single genomic DNA molecule are condensed into short joint products, to unite sequences in cis from widely separated genomic sites. These products can therefore provide information about the arrangement of sequences and variants at a larger scale than established long-read sequencing methods. The method has been useful in defining the phase of variants in haplotypes, the typing of inversions, and determining the configuration of sequence variants in multiallelic CNVs. In this description we outline the rationale for the application of emulsion-fusion PCR methods to the analysis of multiallelic CNVs, and give practical details for our own implementation of the method in that context.

Selective sweep on human amylase genes postdates the split with Neanderthals.

Humans have more copies of amylase genes than other primates. It is still poorly understood, however, when the copy number expansion occurred and whether its spread was enhanced by selection. Here we assess amylase copy numbers in a global sample of 480 high coverage genomes and find that regions flanking the amylase locus show notable depression of genetic diversity both in African and non-African populations. Analysis of genetic variation in these regions supports the model of an early selective sweep in the human lineage after the split of humans from Neanderthals which led to the fixation of multiple copies of AMY1 in place of a single copy. We find evidence of multiple secondary losses of copy number with the highest frequency (52%) of a deletion of AMY2A and associated low copy number of AMY1 in Northeast Siberian populations whose diet has been low in starch content.

Low α-defensin gene copy number increases the risk for IgA nephropathy and renal dysfunction.

Although a major source of genetic variation, copy number variations (CNVs) and their involvement in disease development have not been well studied. Immunoglobulin A nephropathy (IgAN) is the most common primary glomerulonephritis worldwide. We performed association analysis of the DEFA1A3 CNV locus in two independent IgAN cohorts of southern Chinese Han (total of 1189 cases and 1187 controls). We discovered three independent copy number associations within the locus: DEFA1A3 [P = 3.99 × 10(-9); odds ratio (OR), 0.88], DEFA3 (P = 6.55 × 10(-5); OR, 0.82), and a noncoding deletion variant (211bp) (P = 3.50 × 10(-16); OR, 0.75) (OR per copy, fixed-effects meta-analysis). While showing strong association with an increased risk for IgAN (P = 9.56 × 10(-20)), low total copy numbers of the three variants also showed significant association with renal dysfunction in patients with IgAN (P = 0.03; hazards ratio, 3.69; after controlling for the effects of known prognostic factors) and also with increased serum IgA1 (P = 0.02) and galactose-deficient IgA1 (P = 0.03). For replication, we confirmed the associations of DEFA1A3 (P = 4.42 × 10(-4); OR, 0.82) and DEFA3 copy numbers (P = 4.30 × 10(-3); OR, 0.74) with IgAN in a Caucasian cohort (531 cases and 198 controls) and found the 211bp variant to be much rarer in Caucasians. We also observed an association of the 211bp copy number with membranous nephropathy (P = 1.11 × 10(-7); OR, 0.74; in 493 Chinese cases and 500 matched controls), but not with diabetic kidney disease (in 806 Chinese cases and 786 matched controls). By explaining 4.96% of disease risk and influencing renal dysfunction in patients with IgAN, the DEFA1A3 CNV locus may be a potential therapeutic target for developing treatments for this disease.

Heritability of attractiveness to mosquitoes.

Female mosquitoes display preferences for certain individuals over others, which is determined by differences in volatile chemicals produced by the human body and detected by mosquitoes. Body odour can be controlled genetically but the existence of a genetic basis for differential attraction to insects has never been formally demonstrated. This study investigated heritability of attractiveness to mosquitoes by evaluating the response of Aedes aegypti (=Stegomyia aegypti) mosquitoes to odours from the hands of identical and non-identical twins in a dual-choice assay. Volatiles from individuals in an identical twin pair showed a high correlation in attractiveness to mosquitoes, while non-identical twin pairs showed a significantly lower correlation. Overall, there was a strong narrow-sense heritability of 0.62 (SE 0.124) for relative attraction and 0.67 (0.354) for flight activity based on the average of ten measurements. The results demonstrate an underlying genetic component detectable by mosquitoes through olfaction. Understanding the genetic basis for attractiveness could create a more informed approach to repellent development.

Obesity, starch digestion and amylase: association between copy number variants at human salivary (AMY1) and pancreatic (AMY2) amylase genes.

The human salivary amylase genes display extensive copy number variation (CNV), and recent work has implicated this variation in adaptation to starch-rich diets, and in association with body mass index. In this work, we use paralogue ratio tests, microsatellite analysis, read depth and fibre-FISH to demonstrate that human amylase CNV is not a smooth continuum, but is instead partitioned into distinct haplotype classes. There is a fundamental structural distinction between haplotypes containing odd or even numbers of AMY1 gene units, in turn coupled to CNV in pancreatic amylase genes AMY2A and AMY2B. Most haplotypes have one copy each of AMY2A and AMY2B and contain an odd number of copies of AMY1; consequently, most individuals have an even total number of AMY1. In contrast, haplotypes carrying an even number of AMY1 genes have rearrangements leading to CNVs of AMY2A/AMY2B. Read-depth and experimental data show that different populations harbour different proportions of these basic haplotype classes. In Europeans, the copy numbers of AMY1 and AMY2A are correlated, so that phenotypic associations caused by variation in pancreatic amylase copy number could be detected indirectly as weak association with AMY1 copy number. We show that the quantitative polymerase chain reaction (qPCR) assay previously applied to the high-throughput measurement of AMY1 copy number is less accurate than the measures we use and that qPCR data in other studies have been further compromised by systematic miscalibration. Our results uncover new patterns in human amylase variation and imply a potential role for AMY2 CNV in functional associations.

CCL3L1 copy number, CCR5 genotype and susceptibility to tuberculosis.

BACKGROUND: Tuberculosis is a major infectious disease and functional studies have provided evidence that both the chemokine MIP-1α and its receptor CCR5 play a role in susceptibility to TB. Thus by measuring copy number variation of CCL3L1, one of the genes that encode MIP-1α, and genotyping a functional promoter polymorphism -2459A > G in CCR5 (rs1799987) we investigate the influence of MIP-1α and CCR5, independently and combined, in susceptibility to clinically active TB in three populations, a Peruvian population (n = 1132), a !Xhosa population (n = 605) and a South African Coloured population (n = 221). The three populations include patients with clinically diagnosed pulmonary TB, as well as other, less prevalent forms of extrapulmonary TB. METHODS AND RESULTS: Copy number of CCL3L1 was measured using the paralogue ratio test and exhibited ranges between 0-6 copies per diploid genome (pdg) in Peru, between 0-12 pdg in !Xhosa samples and between 0-10 pdg in South African Coloured samples. The CCR5 promoter polymorphism was observed to differ significantly in allele frequency between populations (*A; Peru f = 0.67, !Xhosa f = 0.38, Coloured f = 0.48). CONCLUSIONS: The case-control association studies performed however find, surprisingly, no evidence for an influence of variation in genes coding for MIP-1α or CCR5 individually or together in susceptibility to clinically active TB in these populations.

Accurate measurement of gene copy number for human alpha-defensin DEFA1A3.

BACKGROUND: Multi-allelic copy number variants include examples of extensive variation between individuals in the copy number of important genes, most notably genes involved in immune function. The definition of this variation, and analysis of its impact on function, has been hampered by the technical difficulty of large-scale but accurate typing of genomic copy number. The copy-variable alpha-defensin locus DEFA1A3 on human chromosome 8 commonly varies between 4 and 10 copies per diploid genome, and presents considerable challenges for accurate high-throughput typing. RESULTS: In this study, we developed two paralogue ratio tests and three allelic ratio measurements that, in combination, provide an accurate and scalable method for measurement of DEFA1A3 gene number. We combined information from different measurements in a maximum-likelihood framework which suggests that most samples can be assigned to an integer copy number with high confidence, and applied it to typing 589 unrelated European DNA samples. Typing the members of three-generation pedigrees provided further reassurance that correct integer copy numbers had been assigned. Our results have allowed us to discover that the SNP rs4300027 is strongly associated with DEFA1A3 gene copy number in European samples. CONCLUSIONS: We have developed an accurate and robust method for measurement of DEFA1A3 copy number. Interrogation of rs4300027 and associated SNPs in Genome-Wide Association Study SNP data provides no evidence that alpha-defensin copy number is a strong risk factor for phenotypes such as Crohn's disease, type I diabetes, HIV progression and multiple sclerosis.

Multilocus genetic models of handedness closely resemble single-locus models in explaining family data and are compatible with genome-wide association studies.

Right- and left-handedness run in families, show greater concordance in monozygotic than dizygotic twins, and are well described by single-locus Mendelian models. Here we summarize a large genome-wide association study (GWAS) that finds no significant associations with handedness and is consistent with a meta-analysis of GWASs. The GWAS had 99% power to detect a single locus using the conventional criterion of P < 5 × 10(-8) for the single locus models of McManus and Annett. The strong conclusion is that handedness is not controlled by a single genetic locus. A consideration of the genetic architecture of height, primary ciliary dyskinesia, and intelligence suggests that handedness inheritance can be explained by a multilocus variant of the McManus DC model, classical effects on family and twins being barely distinguishable from the single locus model. Based on the ENGAGE meta-analysis of GWASs, we estimate at least 40 loci are involved in determining handedness.

Determination of haplotypes at structurally complex regions using emulsion haplotype fusion PCR.

BACKGROUND: Genotyping and massively-parallel sequencing projects result in a vast amount of diploid data that is only rarely resolved into its constituent haplotypes. It is nevertheless this phased information that is transmitted from one generation to the next and is most directly associated with biological function and the genetic causes of biological effects. Despite progress made in genome-wide sequencing and phasing algorithms and methods, problems assembling (and reconstructing linear haplotypes in) regions of repetitive DNA and structural variation remain. These dynamic and structurally complex regions are often poorly understood from a sequence point of view. Regions such as these that are highly similar in their sequence tend to be collapsed onto the genome assembly. This is turn means downstream determination of the true sequence haplotype in these regions poses a particular challenge. For structurally complex regions, a more focussed approach to assembling haplotypes may be required. RESULTS: In order to investigate reconstruction of spatial information at structurally complex regions, we have used an emulsion haplotype fusion PCR approach to reproducibly link sequences of up to 1kb in length to allow phasing of multiple variants from neighbouring loci, using allele-specific PCR and sequencing to detect the phase. By using emulsion systems linking flanking regions to amplicons within the CNV, this led to the reconstruction of a 59kb haplotype across the DEFA1A3 CNV in HapMap individuals. CONCLUSION: This study has demonstrated a novel use for emulsion haplotype fusion PCR in addressing the issue of reconstructing structural haplotypes at multiallelic copy variable regions, using the DEFA1A3 locus as an example.

Association of ß-defensin copy number and psoriasis in three cohorts of European origin.

A single previous study has demonstrated significant association of psoriasis with copy number of ß-defensin genes, using DNA from psoriasis cases and controls from Nijmegen and Erlangen. In this study, we attempted to replicate that finding in larger new cohorts from Erlangen (N=2,017) and Michigan (N=5,412), using improved methods for ß-defensin copy number determination based on the paralog ratio test, and enhanced methods of analysis and association testing implemented in the CNVtools resource. We demonstrate that the association with psoriasis found in the discovery sample is maintained after applying improved typing and analysis methods (P=5.5 × 10(-4), odds ratio (OR)=1.25). We also find that the association is replicated in 2,616 cases and 2,526 controls from Michigan, although at reduced significance (P=0.014), but not in new samples from Erlangen (1,396 cases and 621 controls, P=0.38). Meta-analysis across all cohorts suggests a nominally significant association (P=6.6 × 10(-3)/2 × 10(-4)) with an effect size (OR=1.081) much lower than found in the discovery study (OR=1.32). This reduced effect size and significance on replication is consistent with a genuine but weak association.

CCL3L1 copy number and susceptibility to malaria.

Copy number variation can contribute to the variation observed in susceptibility to complex diseases. Here we present the first study to investigate copy number variation of the chemokine gene CCL3L1 with susceptibility to malaria. We present a family-based genetic analysis of a Tanzanian population (n=922), using parasite load, mean number of clinical infections of malaria and haemoglobin levels as phenotypes. Copy number of CCL3L1 was measured using the paralogue ratio test (PRT) and the dataset exhibited copy numbers ranging between 1 and 10 copies per diploid genome (pdg). Association between copy number and phenotypes was assessed. Furthermore, we were able to identify copy number haplotypes in some families, using microsatellites within the copy variable region, for transmission disequilibrium testing. We identified a high level of copy number haplotype diversity and find some evidence for an association of low CCL3L1 copy number with protection from anaemia.

Copy number variation of ß-defensin genes in Behçet's disease.

OBJECTIVES: Behçet's disease (BD) may be triggered by infectious agents in genetically susceptible persons. Human ß-defensin 2 is an inducible antimicrobial peptide, the level of which can be influenced by copy number (CN) of the DEFB4. We investigated the relationship between copy number variation (CNV) of DEFB4 and BD. METHODS: One hundred and ninety-seven patients with BD and 197 healthy controls were enrolled. After measuring CN of DEFB4 with a paralogue ratio test, the CNV was compared between patients and controls. CNV was also analysed in comparison with the clinical manifestations of BD. RESULTS: The CN of DEFB4 was unimodally distributed among the study subjects with mean CN of 4.57 and standard deviation of 1.28. BD samples had numerically lower CN than controls, but the difference was not statistically significant (4.49 ± 1.21 vs. 4.65 ± 1.36, p=0.245). Regarding the relationship between CN of DEFB4 and clinical manifestations, there was no difference of CNV depending on the clinical manifestations. CONCLUSIONS: We found no significant difference in CNV of DEFB4 between patients with BD and controls. Our results suggest that CNV of DEFB4 may not contribute to the pathogenesis of BD.

Meta-analysis confirms the LCE3C_LCE3B deletion as a risk factor for psoriasis in several ethnic groups and finds interaction with HLA-Cw6.

A multicenter meta-analysis including data from 9,389 psoriasis patients and 9,477 control subjects was performed to investigate the contribution of the deletion of genes LCE3C and LCE3B, involved in skin barrier defense, to psoriasis susceptibility in different populations. The study confirms that the deletion of LCE3C and LCE3B is a common genetic factor for susceptibility to psoriasis in the European populations (OR(Overall) = 1.21 (1.15-1.27)), and for the first time directly demonstrates the deletion's association with psoriasis in the Chinese (OR = 1.27 (1.16-1.34)) and Mongolian (OR = 2.08 (1.44-2.99)) populations. The analysis of the HLA-Cw6 locus showed significant differences in the epistatic interaction with the LCE3C and LCE3B deletion in at least some European populations, indicating epistatic effects between these two major genetic contributors to psoriasis. The study highlights the value of examining genetic risk factors in multiple populations to identify genetic interactions, and indicates the need of further studies to understand the interaction of the skin barrier and the immune system in susceptibility to psoriasis.

Screening for common copy-number variants in cancer genes.

For most cases of colorectal cancer that arise without a family history of the disease, it is proposed that an appreciable heritable component of predisposition is the result of contributions from many loci. Although progress has been made in identifying single nucleotide variants associated with colorectal cancer risk, the involvement of low-penetrance copy number variants is relatively unexplored. We have used multiplex amplifiable probe hybridization (MAPH) in a fourfold multiplex (QuadMAPH), positioned at an average resolution of one probe per 2 kb, to screen a total of 1.56 Mb of genomic DNA for copy number variants around the genes APC, AXIN1, BRCA1, BRCA2, CTNNB1, HRAS, MLH1, MSH2, and TP53. Two deletion events were detected, one upstream of MLH1 in a control individual and the other in APC in a colorectal cancer patient, but these do not seem to correspond to copy number polymorphisms with measurably high population frequencies. In summary, by means of our QuadMAPH assay, copy number measurement data were of sufficient resolution and accuracy to detect any copy number variants with high probability. However, this study has demonstrated a very low incidence of deletion and duplication variants within intronic and flanking regions of these nine genes, in both control individuals and colorectal cancer patients.

Measurement methods and accuracy in copy number variation: failure to replicate associations of beta-defensin copy number with Crohn's disease.

The copy number variation in beta-defensin genes on human chromosome 8 has been proposed to underlie susceptibility to inflammatory disorders, but presents considerable challenges for accurate typing on the scale required for adequately powered case-control studies. In this work, we have used accurate methods of copy number typing based on the paralogue ratio test (PRT) to assess beta-defensin copy number in more than 1500 UK DNA samples including more than 1000 cases of Crohn's disease. A subset of 625 samples was typed using both PRT-based methods and standard real-time PCR methods, from which direct comparisons highlight potentially serious shortcomings of a real-time PCR assay for typing this variant. Comparing our PRT-based results with two previous studies based only on real-time PCR, we find no evidence to support the reported association of Crohn's disease with either low or high beta-defensin copy number; furthermore, it is noteworthy that there are disagreements between different studies on the observed frequency distribution of copy number states among European controls. We suggest safeguards to be adopted in assessing and reporting the accuracy of copy number measurement, with particular emphasis on integer clustering of results, to avoid reporting of spurious associations in future case-control studies.