Evidence Of Copy Number Variation Research Articles

An Africa-wide genomic evolution of insecticide resistance in the malaria vector Anopheles funestus involves selective sweeps, copy number variations, gene conversion and transposons.

Insecticide resistance in malaria vectors threatens to reverse recent gains in malaria control. Deciphering patterns of gene flow and resistance evolution in malaria vectors is crucial to improving control strategies and preventing malaria resurgence. A genome-wide survey of Anopheles funestus genetic diversity Africa-wide revealed evidences of a major division between southern Africa and elsewhere, associated with different population histories. Three genomic regions exhibited strong signatures of selective sweeps, each spanning major resistance loci (CYP6P9a/b, GSTe2 and CYP9K1). However, a sharp regional contrast was observed between populations correlating with gene flow barriers. Signatures of complex molecular evolution of resistance were detected with evidence of copy number variation, transposon insertion and a gene conversion between CYP6P9a/b paralog genes. Temporal analyses of samples before and after bed net scale up suggest that these genomic changes are driven by this control intervention. Multiple independent selective sweeps at the same locus in different parts of Africa suggests that local evolution of resistance in malaria vectors may be a greater threat than trans-regional spread of resistance haplotypes.

Adding Complexity to Complexity: Gene Family Evolution in Polyploids

Comparative genomics of non-model organisms has resurrected whole genome duplication from being viewed as a somewhat obscure process that happens in plants to a primary driver of eukaryotic diversification. The shadow of past ploidy increases has left a strong signature of duplicated genes organized into gene families, even in small genomes that have undergone effectively complete rediploidization. Nevertheless, despite continually advancing technologies and bioinformatics pipelines, resolving the fate of duplicate genes remains a substantial challenge. For example, many important recognition processes are driven not only by allelic expansion through retention of duplicates but also by diversification and copy number variation. This creates technical difficulties with assembly to reference genomes and accurate interpretation of homology. Thus, relatively little is known about the impacts of recent polyploidization and hybridization on the evolution of gene families under selective forces that maintain diversity, such as balancing selection. Here we use a complex of species and ploidy levels in the genus Arabidopsis (A. lyrata and A. arenosa) as a model to investigate the evolutionary dynamics of a large and complicated gene family known to be under strong balancing selection: the S-related kinases, which include the female component of genetically controlled self-incompatibility. Specifically, we question: 1) How does diversity of S-related kinase (SRK) alleles in tetraploids compare to that in their close diploid relatives? 2) Is there increased trans-specific polymorphism (i.e. sharing of alleles that transcend speciation, characteristic of balancing selection) in tetraploids compared to diploids due to the higher number of copies they carry? 3) Do these highly variable loci show evidence of introgression among extant species/ploidy levels within or outside known zones of hybridization? 4) Is there evidence for copy number variation among paralogs? We use this example to highlight specific issues to consider when interpreting gene family evolution, particularly in relation to polyploids but also more generally in diploids. We conclude with recommendations for strategies to address the challenges of resolving such complex loci in the future, using advances in deep sequencing approaches.

No Evidence of Copy Number Variation in Acidic Mammalian Chitinase Genes (CHIA) in New World and Old World Monkeys

Copy number variation may be the most common form of structural genetic variation in the genome. Numerous studies have shown that gene copy number variation can correlate with phenotypic variation, where higher copy numbers correspond to increased expression of the protein and vice versa. Examples include some digestive enzyme genes, where variation in copy numbers and protein expression may be related to dietary differences. Increasing the expression of a digestive enzyme through higher gene copy numbers may thus be a potential mechanism for altering an organism’s digestive capabilities. I investigated copy number variation in genes coding for acidic mammalian chitinase, a chitinolytic digestive enzyme that may be used for the digestion of insect exoskeletons, in nonhuman primates with varying levels of insect consumption. I hypothesized that CHIA copy number correlates positively with level of insectivory, predicting higher copy numbers in more insectivorous primates. I assessed copy number variation with the QuantStudio 3D digital PCR platform, in a comparative sample of Old World and New World primate species (N = 10 species, one or two individuals each). Contrary to my prediction, no evidence of copy number variation was found and all species tested had two gene copies per diploid genome. These findings suggest that if acidic mammalian chitinase expression varies according to insect consumption in primates, it may be up- or downregulated through another mechanism.

Targeted capture enrichment and sequencing identifies extensive nucleotide variation in the turkey MHC-B.

Variation in the major histocompatibility complex (MHC) is increasingly associated with disease susceptibility and resistance in avian species of agricultural importance. This variation includes sequence polymorphisms but also structural differences (gene rearrangement) and copy number variation (CNV). The MHC has now been described for multiple galliform species including the best defined assemblies of the chicken (Gallus gallus) and domestic turkey (Meleagris gallopavo). Using this sequence resource, this study applied high-throughput sequencing to investigate MHC variation in turkeys of North America (NA turkeys). An MHC-specific SureSelect (Agilent) capture array was developed, and libraries were created for 14 turkeys representing domestic (commercial bred), heritage breed, and wild turkeys. In addition, a representative of the Ocellated turkey (M. ocellata) and chicken (G. gallus) was included to test cross-species applicability of the capture array allowing for identification of new species-specific polymorphisms. Libraries were hybridized to ∼12 K cRNA baits and the resulting pools were sequenced. On average, 98% of processed reads mapped to the turkey whole genome sequence and 53% to the MHC target. In addition to the MHC, capture hybridization recovered sequences corresponding to other MHC regions. Sequence alignment and de novo assembly indicated the presence of several additional BG genes in the turkey with evidence for CNV. Variant detection identified an average of 2245 polymorphisms per individual for the NA turkeys, 3012 for the Ocellated turkey, and 462 variants in the chicken (RJF-256). This study provides an extensive sequence resource for examining MHC variation and its relation to health of this agriculturally important group of birds.

Ultra‐deep Illumina sequencing accurately identifies MHC class IIb alleles and provides evidence for copy number variation in the guppy (Poecilia reticulata)

We address the bioinformatic issue of accurately separating amplified genes of the major histocompatibility complex (MHC) from artefacts generated during high-throughput sequencing workflows. We fit observed ultra-deep sequencing depths (hundreds to thousands of sequences per amplicon) of allelic variants to expectations from genetic models of copy number variation (CNV). We provide a simple, accurate and repeatable method for genotyping multigene families, evaluating our method via analyses of 209 b of MHC class IIb exon 2 in guppies (Poecilia reticulata). Genotype repeatability for resequenced individuals (N = 49) was high (100%) within the same sequencing run. However, repeatability dropped to 83.7% between independent runs, either because of lower mean amplicon sequencing depth in the initial run or random PCR effects. This highlights the importance of fully independent replicates. Significant improvements in genotyping accuracy were made by greatly reducing type I genotyping error (i.e. accepting an artefact as a true allele), which may occur when using low-depth allele validation thresholds used by previous methods. Only a small amount (4.9%) of type II error (i.e. rejecting a genuine allele as an artefact) was detected through fully independent sequencing runs. We observed 1-6 alleles per individual, and evidence of sharing of alleles across loci. Variation in the total number of MHC class II loci among individuals, both among and within populations was also observed, and some genotypes appeared to be partially hemizygous; total allelic dosage added up to an odd number of allelic copies. Collectively, observations provide evidence of MHC CNV and its complex basis in natural populations.

Integrative genetic, epigenetic and pathological analysis of paraganglioma reveals complex dysregulation of NOTCH signaling

Head and neck paragangliomas, rare neoplasms of the paraganglia composed of nests of neurosecretory and glial cells embedded in vascular stroma, provide a remarkable example of organoid tumor architecture. To identify genes and pathways commonly deregulated in head and neck paraganglioma, we integrated high-density genome-wide copy number variation (CNV) analysis with microRNA and immunomorphological studies. Gene-centric CNV analysis of 24 cases identified a list of 104 genes most significantly targeted by tumor-associated alterations. The “NOTCH signaling pathway” was the most significantly enriched term in the list (P = 0.002 after Bonferroni or Benjamini correction). Expression of the relevant NOTCH pathway proteins in sustentacular (glial), chief (neuroendocrine) and endothelial cells was confirmed by immunohistochemistry in 47 head and neck paraganglioma cases. There were no relationships between level and pattern of NOTCH1/JAG2 protein expression and germline mutation status in the SDH genes, implicated in paraganglioma predisposition, or the presence/absence of immunostaining for SDHB, a surrogate marker of SDH mutations. Interestingly, NOTCH upregulation was observed also in cases with no evidence of CNVs at NOTCH signaling genes, suggesting altered epigenetic modulation of this pathway. To address this issue we performed microarray-based microRNA expression analyses. Notably 5 microRNAs (miR-200a,b,c and miR-34b,c), including those most downregulated in the tumors, correlated to NOTCH signaling and directly targeted NOTCH1 in in vitro experiments using SH-SY5Y neuroblastoma cells. Furthermore, lentiviral transduction of miR-200s and miR-34s in patient-derived primary tympano-jugular paraganglioma cell cultures was associated with NOTCH1 downregulation and increased levels of markers of cell toxicity and cell death. Taken together, our results provide an integrated view of common molecular alterations associated with head and neck paraganglioma and reveal an essential role of NOTCH pathway deregulation in this tumor type.Electronic supplementary materialThe online version of this article (doi:10.1007/s00401-013-1165-y) contains supplementary material, which is available to authorized users.

X‐linked inhibitor of apoptosis single nucleotide polymorphisms and copy number variation are not risk factors for asthma

Aberrant apoptosis in asthma contributes to airway inflammation. Early apoptosis and fragility of airway epithelial cells and delayed apoptosis of inflammatory lymphocytes can cooperate to increase airway inflammation. In this study, single nucleotide polymorphisms (SNPs) and copy number variation (CNV) in the Baculoviral inhibitor of apoptosis protein repeat-containing 4 (BIRC4) gene (which encodes X-linked inhibitor of apoptosis protein) were evaluated for associations with asthma. Asthma cases (n = 203) were identified from Caucasian cohort participants in the North West Adelaide Health Study and matched with 198 controls. Asthma status was defined using self-report of doctor-diagnosed asthma, in conjunction with spirometry and bronchodilator response. Seven SNPs, which spanned the entire BIRC4 gene, were selected for the study on the basis of a haplotype tagging approach. SNPs genotyping was performed on the SEQUENOM MassARRAY iPLEX Gold platform, and genotyping success rate was > 98%. BIRC4 gene CNV was measured using a duplex Taqman qPCR assay, with RNAseP as the reference gene. Alleles and haplotype associations were analysed by logistic regression, assuming an additive genetic model, and adjusted for gender and atopy. BIRC4 gene copy number was determined entirely by gender. All SNPs were in Hardy-Weinberg equilibrium for both case and control females. BIRC4 allele and haplotype frequencies were comparable between asthma cases and controls. There is no evidence of CNV in BIRC4, and BIRC4 is not a susceptibility gene for asthma.

Insights into the genome structure and copy-number variation of Eimeria tenella

BackgroundEimeria is a genus of parasites in the same phylum (Apicomplexa) as human parasites such as Toxoplasma, Cryptosporidium and the malaria parasite Plasmodium. As an apicomplexan whose life-cycle involves a single host, Eimeria is a convenient model for understanding this group of organisms. Although the genomes of the Apicomplexa are diverse, that of Eimeria is unique in being composed of large alternating blocks of sequence with very different characteristics - an arrangement seen in no other organism. This arrangement has impeded efforts to fully sequence the genome of Eimeria, which remains the last of the major apicomplexans to be fully analyzed. In order to increase the value of the genome sequence data and aid in the effort to gain a better understanding of the Eimeria tenella genome, we constructed a whole genome map for the parasite.ResultsA total of 1245 contigs representing 70.0% of the whole genome assembly sequences (Wellcome Trust Sanger Institute) were selected and subjected to marker selection. Subsequently, 2482 HAPPY markers were developed and typed. Of these, 795 were considered as usable markers, and utilized in the construction of a HAPPY map. Markers developed from chromosomally-assigned genes were then integrated into the HAPPY map and this aided the assignment of a number of linkage groups to their respective chromosomes. BAC-end sequences and contigs from whole genome sequencing were also integrated to improve and validate the HAPPY map. This resulted in an integrated HAPPY map consisting of 60 linkage groups that covers approximately half of the estimated 60 Mb genome. Further analysis suggests that the segmental organization first seen in Chromosome 1 is present throughout the genome, with repeat-poor (P) regions alternating with repeat-rich (R) regions. Evidence of copy-number variation between strains was also uncovered.ConclusionsThis paper describes the application of a whole genome mapping method to improve the assembly of the genome of E. tenella from shotgun data, and to help reveal its overall structure. A preliminary assessment of copy-number variation (extra or missing copies of genomic segments) between strains of E. tenella was also carried out. The emerging picture is of a very unusual genome architecture displaying inter-strain copy-number variation. We suggest that these features may be related to the known ability of this parasite to rapidly develop drug resistance.

Androgen Receptor Copy Number Variation and Androgenetic Alopecia: A Case-Control Study

BackgroundThe functional polymorphism that explains the established association of the androgen receptor (AR) with androgenetic alopecia (AGA) remains unidentified, but Copy Number Variation (CNV) might be relevant. CNV involves changes in copy number of large segments of DNA, leading to the altered dosage of gene regulators or genes themselves. Two recent reports indicate regions of CNV in and around AR, and these have not been studied in relation to AGA. The aim of this preliminary case-control study was to determine if AR CNV is associated with AGA, with the hypothesis that CNV is the functional AR variant contributing to this condition.Methodology/Principal FindingsMultiplex Ligation-dependent Probe Amplification was used to screen for CNV in five AR exons and a conserved, non-coding region upstream of AR in 85 men carefully selected as cases and controls for maximal phenotypic contrast. There was no evidence of CNV in AR in any of the cases or controls, and thus no evidence of significant association between AGA and AR CNV.Conclusions/SignificanceThe results suggest this form of genomic variation at the AR locus is unlikely to predispose to AGA.