Copy number variants and selective sweeps in natural populations of the house mouse (Mus musculus domesticus)

Copy–number variants (CNVs) may play an important role in early adaptations, potentially facilitating rapid divergence of populations. We describe an approach to study this question by investigating CNVs present in natural populations of mice in the early stages of divergence and their involvement in selective sweeps. We have analyzed individuals from two recently diverged natural populations of the house mouse (Mus musculus domesticus) from Germany and France using custom, high–density, comparative genome hybridization arrays (CGH) that covered almost 164 Mb and 2444 genes. One thousand eight hundred and sixty one of those genes we previously identified as differentially expressed between these populations, while the expression of the remaining genes was invariant. In total, we identified 1868 CNVs across all 10 samples, 200 bp to 600 kb in size and affecting 424 genic regions. Roughly two thirds of all CNVs found were deletions. We found no enrichment of CNVs among the differentially expressed genes between the populations compared to the invariant ones, nor any meaningful correlation between CNVs and gene expression changes. Among the CNV genes, we found cellular component gene ontology categories of the synapse overrepresented among all the 2444 genes tested. To investigate potential adaptive significance of the CNV regions, we selected six that showed large differences in frequency of CNVs between the two populations and analyzed variation in at least two microsatellites surrounding the loci in a sample of 46 unrelated animals from the same populations collected in field trappings. We identified two loci with large differences in microsatellite heterozygosity (Sfi1 and Glo1/Dnahc8 regions) and one locus with low variation across the populations (Cmah), thus suggesting that these genomic regions might have recently undergone selective sweeps. Interestingly, the Glo1 CNV has previously been implicated in anxiety–like behavior in mice, suggesting a differential evolution of a behavioral trai

Identification of copy number variants from exome sequence data

Background With advances in next generation sequencing technologies and genomic capture techniques, exome sequencing has become a cost-effective approach for mutation detection in genetic diseases. However, computational prediction of copy number variants (CNVs) from exome sequence data is a challenging task. Whilst numerous programs are available, they have different sensitivities, and have low sensitivity to detect smaller CNVs (1–4 exons). Additionally, exonic CNV discovery using standard aCGH has limitations due to the low probe density over exonic regions. The goal of our study was to develop a protocol to detect exonic CNVs (including shorter CNVs that cover 1–4 exons), combining computational prediction algorithms and a high-resolution custom CGH array. Results We used six published CNV prediction programs (ExomeCNV, CONTRA, ExomeCopy, ExomeDepth, CoNIFER, XHMM) and an in-house modification to ExomeCopy and ExomeDepth (ExCopyDepth) for computational CNV prediction on 30 exomes from the 1000 genomes project and 9 exomes from primary immunodeficiency patients. CNV predictions were tested using a custom CGH array designed to capture all exons (exaCGH). After this validation, we next evaluated the computational prediction of shorter CNVs. ExomeCopy and the in-house modified algorithm, ExCopyDepth, showed the highest capability in detecting shorter CNVs. Finally, the performance of each computational program was assessed by calculating the sensitivity and false positive rate. Conclusions In this paper, we assessed the ability of 6 computational programs to predict CNVs, focussing on short (1–4 exon) CNVs. We also tested these predictions using a custom array targeting exons. Based on these results, we propose a protocol to identify and confirm shorter exonic CNVs combining computational prediction algorithms and custom aCGH experiments

A homozygous contiguous gene deletion in chromosome 16p13.3 leads to autosomal recessive osteopetrosis in a Jordanian patient

Human malignant autosomal recessive osteopetrosis (ARO) is a genetically heterogeneous disorder caused by reduced bone resorption by osteoclasts. Mutations in the CLCN7 gene are responsible not only for a substantial portion of ARO patients, but also for other forms of osteopetrosis characterized by different severity and inheritance. The lack of a clear genotype/phenotype correlation makes genetic counselling a tricky issue for CLCN7-dependent osteopetrosis. Here we characterize the first homozygous interstitial deletion in 16p13.3, detected by array Comparative Genomic Hybridization (a-CGH) in an ARO patient of Jordanian origin. The deletion involved other genes beside CLCN7, while the proband displayed a classic ARO phenotype; however her early death did not allow more extensive clinical investigations. The identification of this novel genomic deletion involving a large part of the CLCN7 gene is of clinical relevance, especially in prenatal diagnosis, and suggests the possibility that this kind of mutation has been underestimated so far. This data highlights the need for alternative approaches to genetic analysis also in other ARO-causative genes

Complex X chromosome rearrangement associated with multiorgan autoimmunity

BACKGROUND: Turner syndrome, a congenital condition that affects 1/2,500 births, results from absence or structural alteration of the second sex chromosome. Turner syndrome is usually associated with short stature, gonadal dysgenesis and variable dysmorphic features. The classical 45,X karyotype accounts approximately for half of all patients, the remainder exhibit mosaicism or structural abnormalities of the X chromosome. However, complex intra-X chromosomal rearrangements involving more than three breakpoints are extremely rare. RESULTS: We present a unique case of a novel complex X chromosome rearrangement in a young female patient presenting successively a wide range of autoimmune diseases including insulin dependent diabetes mellitus, Hashimoto's thyroiditis, celiac disease, anaemia perniciosa, possible inner ear disease and severe hair loss. For the genetic evaluation, conventional cytogenetic analysis and FISH with different X specific probes were initially performed. The complexity of these results and the variety of autoimmune problems of the patient prompted us to identify the exact composition and breakpoints of the rearranged X as well as methylation status of the X chromosomes. The high resolution array-CGH (assembly GRCh37/hg19) detected single copy for the whole chromosome X short arm. Two different sized segments of Xq arm were present in three copies: one large size of 80,3 Mb from Xq11.1 to Xq27.3 region and another smaller (11,1 Mb) from Xq27.3 to Xq28 region. An 1,6 Mb Xq27.3 region of the long arm was present in two copies. Southern blot analysis identified a skewed X inactivation with approximately 70:30 % ratios of methylated/unmethylated fragments. The G-band and FISH patterns of the rearranged X suggested the aspect of a restructured i(Xq) chromosome which was shattered and fortuitously repaired. The X-STR genotype analysis of the family detected that the patient inherited intact maternal X chromosome and a rearranged paternal X chromosome. The multiple Xq breakages and fusions as well as inverted duplication would have been expected to cause a severe Turner phenotype. However, the patient lacks many of the classic somatic features of Turner syndrome, instead she presented multiorgan autoimmune diseases. CONCLUSIONS: The clinical data of the presented patient suggest that fragmentation of the i(Xq) chromosome elevates the risk of autoimmune diseases

Rapid identification of homologous recombinants and determination of gene copy number with reference/query pyrosequencing (RQPS)

Manipulating the mouse genome is a widespread technology with important applications in many biological fields ranging from cancer research to developmental biology. Likewise, correlations between copy number variations (CNVs) and human diseases are emerging. We have developed the reference-query pyrosequencing (RQPS) method, which is based on quantitative pyrosequencing and uniquely designed probes containing single nucleotide variations (SNVs), to offer a simple and affordable genotyping solution capable of identifying homologous recombinants independent of the homology arm size, determining the micro-amplification status of endogenous human loci, and quantifying virus/transgene copy number in experimental or commercial species. In addition, we also present a simple pyrosequencing-based protocol that could be used for the enrichment of homologous recombinant embryonic stem (ES) cells

Portrait of Ependymoma Recurrence in Children: Biomarkers of Tumor Progression Identified by Dual-Color Microarray-Based Gene Expression Analysis

BACKGROUND: Children with ependymoma may experience a relapse in up to 50% of cases depending on the extent of resection. Key biological events associated with recurrence are unknown. METHODOLOGY/PRINCIPAL FINDINGS: To discover the biology behind the recurrence of ependymomas, we performed CGHarray and a dual-color gene expression microarray analysis of 17 tumors at diagnosis co-hybridized with the corresponding 27 first or subsequent relapses from the same patient. As treatment and location had only limited influence on specific gene expression changes at relapse, we established a common signature for relapse. Eighty-seven genes showed an absolute fold change ≥2 in at least 50% of relapses and were defined as the gene expression signature of ependymoma recurrence. The most frequently upregulated genes are involved in the kinetochore (ASPM, KIF11) or in neural development (CD133, Wnt and Notch pathways). Metallothionein (MT) genes were downregulated in up to 80% of the recurrences. Quantitative PCR for ASPM, KIF11 and MT3 plus immunohistochemistry for ASPM and MT3 confirmed the microarray results. Immunohistochemistry on an independent series of 24 tumor pairs at diagnosis and at relapse confirmed the decrease of MT3 expression at recurrence in 17/24 tumor pairs (p = 0.002). Conversely, ASPM expression was more frequently positive at relapse (87.5% vs 37.5%, p = 0.03). Loss or deletion of the MT genes cluster was never observed at relapse. Promoter sequencing after bisulfite treatment of DNA from primary tumors and recurrences as well as treatment of short-term ependymoma cells cultures with a demethylating agent showed that methylation was not involved in MT3 downregulation. However, in vitro treatment with a histone deacetylase inhibitor or zinc restored MT3 expression. CONCLUSIONS/SIGNIFICANCE: The most frequent molecular events associated with ependymoma recurrence were over-expression of kinetochore proteins and down-regulation of metallothioneins. Metallothionein-3 expression is epigenetically controlled and can be restored in vitro by histone deacetylase inhibitors

Identification and validation of copy number variants using SNP genotyping arrays from a large clinical cohort.

BACKGROUND: Genotypes obtained with commercial SNP arrays have been extensively used in many large case-control or population-based cohorts for SNP-based genome-wide association studies for a multitude of traits. Yet, these genotypes capture only a small fraction of the variance of the studied traits. Genomic structural variants (GSV) such as Copy Number Variation (CNV) may account for part of the missing heritability, but their comprehensive detection requires either next-generation arrays or sequencing. Sophisticated algorithms that infer CNVs by combining the intensities from SNP-probes for the two alleles can already be used to extract a partial view of such GSV from existing data sets. RESULTS: Here we present several advances to facilitate the latter approach. First, we introduce a novel CNV detection method based on a Gaussian Mixture Model. Second, we propose a new algorithm, PCA merge, for combining copy-number profiles from many individuals into consensus regions. We applied both our new methods as well as existing ones to data from 5612 individuals from the CoLaus study who were genotyped on Affymetrix 500K arrays. We developed a number of procedures in order to evaluate the performance of the different methods. This includes comparison with previously published CNVs as well as using a replication sample of 239 individuals, genotyped with Illumina 550K arrays. We also established a new evaluation procedure that employs the fact that related individuals are expected to share their CNVs more frequently than randomly selected individuals. The ability to detect both rare and common CNVs provides a valuable resource that will facilitate association studies exploring potential phenotypic associations with CNVs. CONCLUSION: Our new methodologies for CNV detection and their evaluation will help in extracting additional information from the large amount of SNP-genotyping data on various cohorts and use this to explore structural variants and their impact on complex traits

Evidence of novel finescale structural variation at autism spectrum disorder candidate loci

Background: Autism spectrum disorders (ASD) represent a group of neurodevelopmental disorders characterized by a core set of social-communicative and behavioral impairments. Gamma-aminobutyric acid (GABA) is the major inhibitory neurotransmitter in the brain, acting primarily via the GABA receptors (GABR). Multiple lines of evidence, including altered GABA and GABA receptor expression in autistic patients, indicate that the GABAergic system may be involved in the etiology of autism. Methods: As copy number variations (CNVs), particularly rare and de novo CNVs, have now been implicated in ASD risk, we examined the GABA receptors and genes in related pathways for structural variation that may be associated with autism. We further extended our candidate gene set to include 19 genes and regions that had either been directly implicated in the autism literature or were directly related (via function or ancestry) to these primary candidates. For the high resolution CNV screen we employed custom-designed 244 k comparative genomic hybridization (CGH) arrays. Collectively, our probes spanned a total of 11 Mb of GABA-related and additional candidate regions with a density of approximately one probe every 200 nucleotides, allowing a theoretical resolution for detection of CNVs of approximately 1 kb or greater on average. One hundred and sixty-eight autism cases and 149 control individuals were screened for structural variants. Prioritized CNV events were confirmed using quantitative PCR, and confirmed loci were evaluated on an additional set of 170 cases and 170 control individuals that were not included in the original discovery set. Loci that remained interesting were subsequently screened via quantitative PCR on an additional set of 755 cases and 1,809 unaffected family members. Results: Results include rare deletions in autistic individuals at JAKMIP1, NRXN1, Neuroligin4Y, OXTR, and ABAT. Common insertion/deletion polymorphisms were detected at several loci, including GABBR2 and NRXN3. Overall, statistically significant enrichment in affected vs. unaffected individuals was observed for NRXN1 deletions. Conclusions: These results provide additional support for the role of rare structural variation in ASD

Large-scale Identification of Chemically Induced Mutations in Drosophila melanogaster.

Forward genetic screens using chemical mutagens have been successful in defining the function of thousands of genes in eukaryotic model organisms. The main drawback of this strategy is the time-consuming identification of the molecular lesions causative of the phenotypes of interest. With whole-genome sequencing (WGS), it is now possible to sequence hundreds of strains, but determining which mutations are causative among thousands of polymorphisms remains challenging. We have sequenced 394 mutant strains, generated in a chemical mutagenesis screen, for essential genes on the Drosophila X chromosome and describe strategies to reduce the number of candidate mutations from an average of -3500 to 35 single-nucleotide variants per chromosome. By combining WGS with a rough mapping method based on large duplications, we were able to map 274 (-70%) mutations. We show that these mutations are causative, using small 80-kb duplications that rescue lethality. Hence, our findings demonstrate that combining rough mapping with WGS dramatically expands the toolkit necessary for assigning function to genes

Mesenchymal Transition and PDGFRA Amplification/Mutation Are Key Distinct Oncogenic Events in Pediatric Diffuse Intrinsic Pontine Gliomas

Diffuse intrinsic pontine glioma (DIPG) is one of the most frequent malignant pediatric brain tumor and its prognosis is universaly fatal. No significant improvement has been made in last thirty years over the standard treatment with radiotherapy. To address the paucity of understanding of DIPGs, we have carried out integrated molecular profiling of a large series of samples obtained with stereotactic biopsy at diagnosis. While chromosomal imbalances did not distinguish DIPG and supratentorial tumors on CGHarrays, gene expression profiling revealed clear differences between them, with brainstem gliomas resembling midline/thalamic tumours, indicating a closely-related origin. Two distinct subgroups of DIPG were identified. The first subgroup displayed mesenchymal and pro-angiogenic characteristics, with stem cell markers enrichment consistent with the possibility to grow tumor stem cells from these biopsies. The other subgroup displayed oligodendroglial features, and appeared largely driven by PDGFRA, in particular through amplification and/or novel missense mutations in the extracellular domain. Patients in this later group had a significantly worse outcome with an hazard ratio for early deaths, ie before 10 months, 8 fold greater that the ones in the other subgroup (p = 0.041, Cox regression model). The worse outcome of patients with the oligodendroglial type of tumors was confirmed on a series of 55 paraffin-embedded biopsy samples at diagnosis (median OS of 7.73 versus 12.37 months, p = 0.045, log-rank test). Two distinct transcriptional subclasses of DIPG with specific genomic alterations can be defined at diagnosis by oligodendroglial differentiation or mesenchymal transition, respectively. Classifying these tumors by signal transduction pathway activation and by mutation in pathway member genes may be particularily valuable for the development of targeted therapies