Human Copy Number Variation Research Articles

Mosaic variegated aneuploidy in development, ageing and cancer.

Mosaic variegated aneuploidy (MVA) is a rare condition in which abnormal chromosome counts (that is, aneuploidies), affecting different chromosomes in each cell (making it variegated) are found only in a certain number of cells (making it mosaic). MVA is characterized by various developmental defects and, despite its rarity, presents a unique clinical scenario to understand the consequences of chromosomal instability and copy number variation in humans. Research from patients with MVA, genetically engineered mouse models and functional cellular studies have found the genetic causes to be mutations in components of the spindle-assembly checkpoint as well as in related proteins involved in centrosome dynamics during mitosis. MVA is accompanied by tumour susceptibility (depending on the genetic basis) as well as cellular and systemic stress, including chronic immune response and the associated clinical implications.

Parallel Genomic Analysis from Paired Bone Marrow and Peripheral Blood Samples of 200 Cytopenic Patients

Background: Analysis of bone marrow is the gold standard in patients with unexplainable cytopenia for diagnosing (or excluding) myelodysplastic neoplasm (MDS). Besides morphology also molecular genetics and cytogenetics are currently assessed from BM as both are required for diagnosis and risk stratification. Substitution of BM analysis with peripheral blood (PB) analysis would be a substantial progress especially in the elderly to allow for a less invasive strategy. Aims: Validate, whether PB reliably reflects molecular and cytogenetic aberrations of BM by evaluating somatic mutations and copy number variations (CNVs) in concurrent PB and BM samples. Methods: We analyzed paired BM and PB samples (time difference of sampling: <15 days) of 200 cytopenic patients. All samples were subjected to targeted panel sequencing analyzing mutations (assessed as positive at variant allelic frequencies/VAFs >1%) in 40 genes associated with myeloid neoplasms (coverage: 1500x). If available, BM chromosome analysis (CA) and BM FISH analysis were compared to NGS CNV analysis from PB using a human CNV backbone hybridization panel (coverage: 50x). Results: In 69 patients (35%) no mutations (MUTs) were detected in either BM or PB (Figure 1A). For the remaining 131 patients at least one MUT was found in BM (median MUTs/patient: 2; range 1-10). In 130/131 patients MUTs were detected in both BM and PB (same number of MUTs: n=121; higher number of MUTs in BM than in PB: n=9) while in only one patient harboring a DNMT3A MUT as the sole MUT in BM (VAF: 9%) no MUT was detected in PB. This resulted in a negative predictive value of 98.6% (69/70). Overall, we detected 364 MUTs in 29 of the analyzed genes (BM MUTs: n=362; PB MUTs: n=349) and a concordance between BM and PB of 95% (347/364; Figure 1A). The median VAF was lower in PB compared to BM across all MUTs (median VAF: 20% vs. 12%; p<0.001). All MUTs detected in BM at VAF ≥13% (n=229) were also detected in corresponding PB. Next, we compared CNVs detected by PB NGS to BM karyotype (KT). For 119/200 patients PB NGS including a CNV-backbone panel was available. BM CA revealed a normal KT in 71 and an aberrant KT in 48 patients. In total, 126 aberrations were detected in BM CA including 7 balanced aberrations that can technically not be detected using the CNV backbone panel. CNV analysis from PB revealed 68 aberrations (of 119 detectable BM CNVs; 54%) with 64 being concordant with BM CA. In total, in 88/119 patients (74%) concordant results were found in BM and PB (70/71 with normal and 18/48 with aberrant BM KT). In 3 patients (2.5%) different CNVs were detected (1/3: no CNV in BM CA), while in 28 (23.5%) either fewer CNVs (n=8; all complex KTs) or no CNVs (n=20) were detected in PB compared to BM. The majority (76%) of BM CNVs not detected by PB NGS with available BM FISH data were small BM clones (clone size ≤30%) and therefore not reliably captured by NGS PB with a technical sensitivity of 30%. Interestingly, 4 of the discordant patients with aberrant BM KT did not harbor any somatic MUT in BM and PB. Regarding the IPSS-R cytogenetic risk group, 95 (80%) patients fell into the same category, while 22 (18%; all with aberrant BM KT) were assigned to a better risk category by NGS PB than based on BM CA data and 2 patients (2%) to a worse risk category by NGS PB (Figure 1B). Conclusion: In cytopenic patients mutations can be reliably detected from PB showing a concordance with BM of individual mutations of 95% and a consistency regarding overall mutation detection between PB and BM in 99.5% of patients. Cytogenetic risk group assignment based on NGS CNV analysis from PB compared to chromosome analysis from BM resulted in the same risk group in 80% of patients. However, to further decrease the number of discrepant cases, the technical limitations (sensitivity, inability to detect balanced CNVs) of the CNV backbone method might be overcome by increasing the sequencing depth and/or using CD34-enriched PB samples and adding probe sets. Alternatively, FISH on PB smears could be used for CNV detection. Overall, our data indicate a high degree of overlap between PB and BM regarding mutational and cytogenetic analysis in patients with unexplained cytopenia and suggest PB as reliable surrogate especially in patients without any genetic aberrations.

A-242 Integrating Respiratory Metagenomics and Metatranscriptomics for Diagnosis of Lung Cancer and Infection in Patients with Pulmonary Diseases

Abstract Background Advances in unbiased metagenomic next generation sequencing (mNGS) technologies have enabled the study of microbial and host genetic material (DNA and RNA) in one test. In this study, we aimed to develop machine learning-based differential diagnostic models (MLBDDMs) using the metagenomic and human transcriptomic data generated by an affordable bronchoalveolar lavage fluid (BLAF) mNGS assay and investigated their clinical utility for early differential diagnosis of lung cancer and pulmonary infection in patients with pulmonary diseases. Methods We recruited 775 patients with respiratory disease, including 160 pathologically diagnosed lung cancer and clinically diagnosed 615 infectious causes (131 tuberculosis, 172 fungal pneumonia and 312 bacterial pneumonia). An affordable mNGS assay on BALF samples collected from these patients on admission were performed. Using the generated mNGS data, we compared the differences in microbial diversity and host gene expression between lung cancer patients and pulmonary infection patients. The BLAF mNGS datasets of lung cancer group and each infection group were then randomly divided into a training dataset and a validation dataset at a ratio of approximately 3:1 for developing optimal MLBDDMs that can be used to distinguish lung cancer from various pulmonary infections. Results By comparing the BALF mNGS data of lung cancer (n = 160) and pulmonary infection (n = 615), we found that the infection group had higher microbial diversity than lung cancer group (P-value < 0.05). Respiratory colonizing microorganisms (e.g., Corynebacterium propinquum and Bacteroides uniformis) and pathogen (Mycobacterium tuberculosis and Cryptococcus neoformans) were found as differential microbes (adjusted p-value < 0.05, LDA score > 2). From BALF gene expression data, we detected 175 genes enriched in NOD-like receptor signaling pathway and chemokine signaling pathway differentially expressed between lung cancer and pulmonary infection groups (False Discovery Rate, FDR < 0.05). Cell composition analysis revealed that macrophage M1 was higher in lung infection group (P-value < 0.001), whereas mast cell activated and DCs activated were higher in lung cancer group (P-value < 0.001, P-value = 0.016). We integrated the metagenomic (microbial composition and human copy number variation) and transcriptomic data (host differentially expressed genes and cell composition) generated by the BALF mNGS assay with eleven machine learning classifiers to establish diagnosis models for distinguishing lung cancer from pulmonary infection (we named LC/PI model). The results showed that a Random Forest diagnostic model (the RF-LC/PI model) had optimal performance, with a sensitivity and specificity of 86.7% and 87.8%, respectively, in distinguishing lung cancer from pulmonary infection (area under the receiver operating characteristic curve [AUC] = 0.838 in the training dataset; AUC = 0.79 in a held-out validation dataset). Similar to the establishment of the LC/PI model, we further developed three diagnostic models for distinguishing lung cancer and tuberculosis (LC/TB model), lung cancer and fungal pneumonia (LC/FP model), and lung cancer and bacterial pneumonia (LC/BP model), respectively. The AUC of these three models were 0.91, 0.88, 0.91, respectively, showing a high differential diagnosis accuracy. Conclusions We have established MLBDDMs using BALF metagenomic and metatranscriptomic data and achieved superior accuracy for differentiating lung cancer and pulmonary infections, which could promote early diagnosis of pulmonary diseases and benefit more patients with one test.

Genome-wide analysis of copy-number variation in humans with cleft lip and/or cleft palate identifies COBLL1, RIC1, and ARHGEF38 as clefting genes

Cleft lip with or without cleft palate (CL/P) is a common birth defect with a complex, heterogeneous etiology. It is well established that common and rare sequence variants contribute to the formation of CL/P, but the contribution of copy-number variants (CNVs) to cleft formation remains relatively understudied. To fill this knowledge gap, we conducted a large-scale comparative analysis of genome-wide CNV profiles of 869 individuals from the Philippines and 233 individuals of European ancestry with CL/P with three primary goals: first, to evaluate whether differences in CNV number, amount of genomic content, or amount of coding genomic content existed within clefting subtypes; second, to assess whether CNVs in our cohort overlapped with known Mendelian clefting loci; and third, to identify unestablished Mendelian clefting genes. Significant differences in CNVs across cleft types or in individuals with non-syndromic versus syndromic clefts were not observed; however, several CNVs in our cohort overlapped with known syndromic and non-syndromic Mendelian clefting loci. Moreover, employing a filtering strategy relying on population genetics data that rare variants are on the whole more deleterious than common variants, we identify several CNV-associated gene losses likely driving non-syndromic clefting phenotypes. By prioritizing genes deleted at a rare frequency across multiple individuals with clefts yet enriched in our cohort of individuals with clefts compared to control subjects, we identify COBLL1, RIC1, and ARHGEF38 as clefting genes. CRISPR-Cas9 mutagenesis of these genes in Xenopus laevis and Danio rerio yielded craniofacial dysmorphologies, including clefts analogous to those seen in human clefting disorders.

Segmental duplications and their variation in a complete human genome.

Despite their importance in disease and evolution, highly identical segmental duplications (SDs) are among the last regions of the human reference genome (GRCh38) to be fully sequenced. Using a complete telomere-to-telomere human genome (T2T-CHM13), we present a comprehensive view of human SD organization. SDs account for nearly one-third of the additional sequence, increasing the genome-wide estimate from 5.4 to 7.0% [218 million base pairs (Mbp)]. An analysis of 268 human genomes shows that 91% of the previously unresolved T2T-CHM13 SD sequence (68.3 Mbp) better represents human copy number variation. Comparing long-read assemblies from human (n = 12) and nonhuman primate (n = 5) genomes, we systematically reconstruct the evolution and structural haplotype diversity of biomedically relevant and duplicated genes. This analysis reveals patterns of structural heterozygosity and evolutionary differences in SD organization between humans and other primates.

CNVIntegrate: the first multi-ethnic database for identifying copy number variations associated with cancer.

Human copy number variations (CNVs) and copy number alterations (CNAs) are DNA segments (>1000 base pairs) of duplications or deletions with respect to the reference genome, potentially causing genomic imbalance leading to diseases such as cancer. CNVs further cause genetic diversity in healthy populations and are predominant drivers of gene/genome evolution. Initiatives have been taken by the research community to establish large-scale databases to comprehensively characterize CNVs in humans. Exome Aggregation Consortium (ExAC) is one such endeavor that catalogs CNVs, of nearly 60 000 healthy individuals across five demographic clusters. Furthermore, large projects such as the Catalogue of Somatic Mutations in Cancer (COSMIC) and the Cancer Cell Line Encyclopedia (CCLE) combine CNA data from cancer-affected individuals and large panels of human cancer cell lines, respectively. However, we lack a structured and comprehensive CNV/CNA resource including both healthy individuals and cancer patients across large populations. CNVIntegrate is the first web-based system that hosts CNV and CNA data from both healthy populations and cancer patients, respectively, and concomitantly provides statistical comparisons between copy number frequencies of multiple ethnic populations. It further includes, for the first time, well-cataloged CNV and CNA data from Taiwanese healthy individuals and Taiwan Breast Cancer data, respectively, along with imported resources from ExAC, COSMIC and CCLE. CNVIntegrate offers a CNV/CNA-data hub for structured information retrieval for clinicians and scientists towards important drug discoveries and precision treatments. Database URL: http://cnvintegrate.cgm.ntu.edu.tw/

Phenotypic Diversity of 15q11.2 BP1–BP2 Deletion in Three Korean Families with Development Delay and/or Intellectual Disability: A Case Series and Literature Review

The 15q11.2 breakpoint (BP) 1–BP2 deletion syndrome is emerging as the most frequent pathogenic copy number variation in humans related to neurodevelopmental diseases, with changes in cognition, behavior, and brain morphology. Previous publications have reported that patients with 15q11.2 BP1–BP2 deletion showed intellectual disability (ID), speech impairment, developmental delay (DD), and/or behavioral problems. We describe three new cases, aged 3 or 6 years old and belonging to three unrelated Korean families, with a 350-kb 15q11.2 BP1–BP2 deletion of four highly conserved genes, namely, the TUBGCP5, CYFIP1, NIPA2, and NIPA1 genes. All of our cases presented with global DD and/or ID, and the severity ranged from mild to severe, but common facial dysmorphism and congenital malformations in previous reports were not characteristic. The 15q11.2 BP1–BP2 deletion was inherited from an unaffected parent in all cases. Our three cases, together with previous findings from the literature review, confirm some of the features earlier reported to be associated with 15q11.2 BP1–BP2 deletion and help to further delineate the phenotype associated with 15q11.2 deletion. Identification of more cases with 15q11.2 BP1–BP2 deletion will allow us to obtain a better understanding of the clinical phenotypes. Further explanation of the functions of the genes within the 15q11.2 BP1–BP2 region is required to resolve the pathogenic effects on neurodevelopment.

The ELIXIR Human Copy Number Variations Community: building bioinformatics infrastructure for research

Copy number variations (CNVs) are major causative contributors both in the genesis of genetic diseases and human neoplasias. While “High-Throughput” sequencing technologies are increasingly becoming the primary choice for genomic screening analysis, their ability to efficiently detect CNVs is still heterogeneous and remains to be developed. The aim of this white paper is to provide a guiding framework for the future contributions of ELIXIR’s recently established h uman CNV Community, with implications beyond human disease diagnostics and population genomics. This white paper is the direct result of a strategy meeting that took place in September 2018 in Hinxton (UK) and involved representatives of 11 ELIXIR Nodes. The meeting led to the definition of priority objectives and tasks, to address a wide range of CNV-related challenges ranging from detection and interpretation to sharing and training. Here, we provide suggestions on how to align these tasks within the ELIXIR Platforms strategy, and on how to frame the activities of this new ELIXIR Community in the international context.

The ELIXIR Human Copy Number Variations Community: building bioinformatics infrastructure for research.

Copy number variations (CNVs) are major causative contributors both in the genesis of genetic diseases and human neoplasias. While "High-Throughput" sequencing technologies are increasingly becoming the primary choice for genomic screening analysis, their ability to efficiently detect CNVs is still heterogeneous and remains to be developed. The aim of this white paper is to provide a guiding framework for the future contributions of ELIXIR's recently established human CNV Community, with implications beyond human disease diagnostics and population genomics. This white paper is the direct result of a strategy meeting that took place in September 2018 in Hinxton (UK) and involved representatives of 11 ELIXIR Nodes. The meeting led to the definition of priority objectives and tasks, to address a wide range of CNV-related challenges ranging from detection and interpretation to sharing and training. Here, we provide suggestions on how to align these tasks within the ELIXIR Platforms strategy, and on how to frame the activities of this new ELIXIR Community in the international context.

Drosophila models of pathogenic copy-number variant genes show global and non-neuronal defects during development.

While rare pathogenic copy-number variants (CNVs) are associated with both neuronal and non-neuronal phenotypes, functional studies evaluating these regions have focused on the molecular basis of neuronal defects. We report a systematic functional analysis of non-neuronal defects for homologs of 59 genes within ten pathogenic CNVs and 20 neurodevelopmental genes in Drosophila melanogaster. Using wing-specific knockdown of 136 RNA interference lines, we identified qualitative and quantitative phenotypes in 72/79 homologs, including 21 lines with severe wing defects and six lines with lethality. In fact, we found that 10/31 homologs of CNV genes also showed complete or partial lethality at larval or pupal stages with ubiquitous knockdown. Comparisons between eye and wing-specific knockdown of 37/45 homologs showed both neuronal and non-neuronal defects, but with no correlation in the severity of defects. We further observed disruptions in cell proliferation and apoptosis in larval wing discs for 23/27 homologs, and altered Wnt, Hedgehog and Notch signaling for 9/14 homologs, including AATF/Aatf, PPP4C/Pp4-19C, and KIF11/Klp61F. These findings were further supported by tissue-specific differences in expression patterns of human CNV genes, as well as connectivity of CNV genes to signaling pathway genes in brain, heart and kidney-specific networks. Our findings suggest that multiple genes within each CNV differentially affect both global and tissue-specific developmental processes within conserved pathways, and that their roles are not restricted to neuronal functions.

Transcriptomic networks implicate neuronal energetic abnormalities in three mouse models harboring autism and schizophrenia-associated mutations.

Genetic risk for psychiatric illness is complex, so identification of shared molecular pathways where distinct forms of genetic risk might coincide is of substantial interest. A growing body of genetic and genomic studies suggest that such shared molecular pathways exist across disorders with different clinical presentations, such as schizophrenia and autism spectrum disorder (ASD). But how this relates to specific genetic risk factors is unknown. Further, whether some of the molecular changes identified in brain relate to potentially confounding antemortem or postmortem factors are difficult to prove. We analyzed the transcriptome from the cortex and hippocampus of three mouse lines modeling human copy number variants (CNVs) associated with schizophrenia and ASD: Df(h15q13)/+, Df(h22q11)/+, and Df(h1q21)/+ which carry the 15q13.3 deletion, 22q11.2 deletion, and 1q21.1 deletion, respectively. Although we found very little overlap of differential expression at the level of individual genes, gene network analysis identified two cortical and two hippocampal modules of co-expressed genes that were dysregulated across all three mouse models. One cortical module was associated with neuronal energetics and firing rate, and overlapped with changes identified in postmortem human brain from SCZ and ASD patients. These data highlight aspects of convergent gene expression in mouse models harboring major risk alleles, and strengthen the connection between changes inneuronal energetics and neuropsychiatric disorders in humans.

Dosage sensitivity intolerance of VIPR2 microduplication is disease causative to manifest schizophrenia-like phenotypes in a novel BAC transgenic mouse model.

Recent genome-wide association studies (GWAS) have identified copy number variations (CNVs) at chromosomal locus 7q36.3 that significantly contribute to the risk of schizophrenia, with all of the microduplications occurring within a single gene: vasoactive intestinal peptide receptor 2 (VIPR2). To confirm disease causality and translate such a genetic vulnerability into mechanistic and pathophysiological insights, we have developed a series of conditional VIPR2 bacterial artificial chromosome (BAC) transgenic mouse models of VIPR2 CNV. VIPR2 CNV mouse model recapitulates gene expression and signaling deficits seen in human CNV carriers. VIPR2 microduplication in mice elicits prominent dorsal striatal dopamine dysfunction, cognitive, sensorimotor gating, and social behavioral deficits preceded by an increase of striatal cAMP/PKA signaling and the disrupted early postnatal striatal development. Genetic removal of VIPR2 transgene expression via crossing with Drd1a-Cre BAC transgenic mice rescued the dopamine D2 receptor abnormality and multiple behavioral deficits, implicating a pathogenic role of VIPR2 overexpression in dopaminoceptive neurons. Thus, our results provide further evidence to support the GWAS studies that the dosage sensitivity intolerance of VIPR2 is disease causative to manifest schizophrenia-like dopamine, cognitive, and social behavioral deficits in mice. The conditional BAC transgenesis offers a novel strategy to model CNVs with a gain-of -copies and facilitate the genetic dissection of when/where/how the genetic vulnerabilities affect development, structure, and function of neural circuits. Our findings have important implications for therapeutic development, and the etiology-relevant mouse model provides a useful preclinical platform for drug discovery.

Origins of Intraindividual Genetic Variation in Human Fetuses

BackgroundIntraindividual copy number variation (CNV) origin is largely unknown. They might be due to aging and/or common genome instability at the preimplantation stage while contribution of preimplantation in human intraindividual CNVs occurrence is unknown. To address this question, we investigated mosaicism and its origin in the fetuses of natural conception.MethodsWe studied normal fetuses following therapeutic abortion due to maternal indications. We analyzed the genome of 22 tissues of each fetus by array comparative genomic hybridization for intraindividual CNVs. Each tissue was studied in 2 microarray experiments; the reciprocal aberrations larger than 40 Kb, identified by comparing tissues of each fetus, were subsequently validated using quantitative polymerase chain reaction.ResultsThrough intraindividual comparison, frequency of reciprocal events varied from 2 to 9. According to the distribution pattern of the frequent CNV in derivatives of different germ layers, we found that its origin is early development including preimplantation, whereas CNVs with low frequency have occurred in later stages. Shared CNVs in both fetuses were belonged to thymus and related to the functional role of genes located in these CNVs.ConclusionsThe origin of some of fetal CNVs is preimplantation stage. Each organ might inherit CNVs with an unpredictable pattern due to the extensive cell mixing/migration in embryonic development.

Genome plasticity in Candida albicans is driven by long repeat sequences.

Genome rearrangements resulting in copy number variation (CNV) and loss of heterozygosity (LOH) are frequently observed during the somatic evolution of cancer and promote rapid adaptation of fungi to novel environments. In the human fungal pathogen Candida albicans, CNV and LOH confer increased virulence and antifungal drug resistance, yet the mechanisms driving these rearrangements are not completely understood. Here, we unveil an extensive array of long repeat sequences (65-6499 bp) that are associated with CNV, LOH, and chromosomal inversions. Many of these long repeat sequences are uncharacterized and encompass one or more coding sequences that are actively transcribed. Repeats associated with genome rearrangements are predominantly inverted and separated by up to ~1.6 Mb, an extraordinary distance for homology-based DNA repair/recombination in yeast. These repeat sequences are a significant source of genome plasticity across diverse strain backgrounds including clinical, environmental, and experimentally evolved isolates, and represent previously uncharacterized variation in the reference genome.

Copy Number Analysis of Lung Tumor Progression in Mice

Background & Hypothesis 
 The xeroderma pigmentosum group C (XPC) gene encodes a protein that repairs DNA damage, like that caused by UV light. Its mutation in patients is associated with increased risk of melanoma. In our mouse model, the Xpc gene has a targeted deletion, leading to increased frequency of lung tumors in a mouse strain typically resistant to lung chemical carcinogenesis. We hypothesize that there will be enrichment of specific copy number variations (CNVs) in the genome throughout tumor progression that can distinguish the different stages of carcinogenesis in lung cancer. 
 Project Methods 
 Urethane treatments induced lung tumors in Xpc -/- (ko) mice. Laser capture microscopy was used to isolate tumor, hyperplasia, and normal tissue from lung sections. After DNA extraction, genomic DNA (gDNA) was quantified using Quantifluor assay. Whole genome amplification was completed using Repli-G kit. PCR, copy number analysis, and restriction enzyme digestion were quality control measures for amplified gDNA. Single nucleotide polymorphism (SNP) microarrays will assess CNVs in samples that are successfully amplified. 
 Results 
 DNA was successfully amplified, as confirmed by the Quantifluor assay. PCR, copy number analysis, and restriction enzyme digestion confirmed 3 out of the 4 amplified samples were comparable to matched non-amplified samples, suggesting mouse gDNA was maintained. 
 Potential Impact 
 The mouse model findings can be applied to publically available human genome databases of lung cancer CNVs. Overlapping findings in mice and human CNVs may give insight into novel pathways for effective treatment in lung cancer.

Human copy number variants are enriched in regions of low mappability.

Copy number variants (CNVs) are known to affect a large portion of the human genome and have been implicated in many diseases. Although whole-genome sequencing (WGS) can help identify CNVs, most analytical methods suffer from limited sensitivity and specificity, especially in regions of low mappability. To address this, we use PopSV, a CNV caller that relies on multiple samples to control for technical variation. We demonstrate that our calls are stable across different types of repeat-rich regions and validate the accuracy of our predictions using orthogonal approaches. Applying PopSV to 640 human genomes, we find that low-mappability regions are approximately 5 times more likely to harbor germline CNVs, in stark contrast to the nearly uniform distribution observed for somatic CNVs in 95 cancer genomes. In addition to known enrichments in segmental duplication and near centromeres and telomeres, we also report that CNVs are enriched in specific types of satellite and in some of the most recent families of transposable elements. Finally, using this comprehensive approach, we identify 3455 regions with recurrent CNVs that were missing from existing catalogs. In particular, we identify 347 genes with a novel exonic CNV in low-mappability regions, including 29 genes previously associated with disease.

Conserved microRNA targeting reveals preexisting gene dosage sensitivities that shaped amniote sex chromosome evolution.

Mammalian X and Y Chromosomes evolved from an ordinary autosomal pair. Genetic decay of the Y led to X Chromosome inactivation (XCI) in females, but some Y-linked genes were retained during the course of sex chromosome evolution, and many X-linked genes did not become subject to XCI. We reconstructed gene-by-gene dosage sensitivities on the ancestral autosomes through phylogenetic analysis of microRNA (miRNA) target sites and compared these preexisting characteristics to the current status of Y-linked and X-linked genes in mammals. Preexisting heterogeneities in dosage sensitivity, manifesting as differences in the extent of miRNA-mediated repression, predicted either the retention of a Y homolog or the acquisition of XCI following Y gene decay. Analogous heterogeneities among avian Z-linked genes predicted either the retention of a W homolog or gene-specific dosage compensation following W gene decay. Genome-wide analyses of human copy number variation indicate that these heterogeneities consisted of sensitivity to both increases and decreases in dosage. We propose a model of XY/ZW evolution incorporating such preexisting dosage sensitivities in determining the evolutionary fates of individual genes. Our findings thus provide a more complete view of the role of dosage sensitivity in shaping the mammalian and avian sex chromosomes and reveal an important role for post-transcriptional regulatory sequences (miRNA target sites) in sex chromosome evolution.

Dosage sensitivity is a major determinant of human copy number variant pathogenicity

Human copy number variants (CNVs) account for genome variation an order of magnitude larger than single-nucleotide polymorphisms. Although much of this variation has no phenotypic consequences, some variants have been associated with disease, in particular neurodevelopmental disorders. Pathogenic CNVs are typically very large and contain multiple genes, and understanding the cause of the pathogenicity remains a major challenge. Here we show that pathogenic CNVs are significantly enriched for genes involved in development and genes that have greater evolutionary copy number conservation across mammals, indicative of functional constraints. Conversely, genes found in benign CNV regions have more variable copy number. These evolutionary constraints are characteristic of genes in pathogenic CNVs and can only be explained by dosage sensitivity of those genes. These results implicate dosage sensitivity of individual genes as a common cause of CNV pathogenicity. These evolutionary metrics suggest a path to identifying disease genes in pathogenic CNVs.

Association of common copy number variations with diseases

Genomic polymorphisms come in various forms including single nucleotide variations, translocations, insertions and copy number variations (CNVs). As a form of structural variation, the CNVs comprise common and rare forms based on their populational frequencies. Studies have demonstrated that certain CNVs are associated with risks for neuro-developmental diseases, viral infections, chronic inflammations, and cancers. With the development of high-resolution genome typing technologies such as microarrays and whole genome sequencing, the human genomic CNVs map has been continuously improved and refined. In-depth study of CNVs not only can provide comprehensive understanding for their structural variations and genetic evolution, but also provide new insights into genetic factors contributing to such diseases. In this paper, the general characteristics, pathogenesis and detection methods for the CNVs, as well as their association with human diseases are reviewed.

Identifying rhesus macaque gene orthologs using heterospecific human CNV probes.

We used the Affymetrix® Genome-Wide Human SNP Array 6.0 to identify heterospecific markers and compare copy number and structural genomic variation between humans and rhesus macaques. Over 200,000 human copy number variation (CNV) probes were mapped to a Chinese and an Indian rhesus macaque sample. Observed genomic rearrangements and synteny were in agreement with the results of a previously published genomic comparison between humans and rhesus macaques. Comparisons between each of the two rhesus macaques and humans yielded 206 regions with copy numbers that differed by at least two fold in the Indian rhesus macaque and human, 32 in the Chinese rhesus macaque and human, and 147 in both rhesus macaques. The detailed genomic map and preliminary CNV data are useful for better understanding genetic variation in rhesus macaques, identifying derived changes in human CNVs that may have evolved by selection, and determining the suitability of rhesus macaques as human models for particular biomedical studies.