Segmental copy number variation shapes tissue transcriptomes

Inference of Causative Genes for Alzheimer's Disease Due to Dosage Imbalance

Copy number variations (CNVs) have recently drawn attention as an important genetic factor for diseases, especially common neuropsychiatric disorders including Alzheimer's disease (AD). Because most of the pathogenic CNV regions overlap with multiple genes, it has been challenging to identify the true disease-causing genes amongst them. Notably, a recent study reported that CNV regions containing ohnologs, which are dosage-sensitive genes, are likely to be deleterious. Utilizing the unique feature of ohnologs could be useful for identifying causative genes with pathogenic CNVs, however its effectiveness is still unclear. Although it has been reported that AD is strongly affected by CNVs, most of AD-causing genes with pathogenic CNVs have not been identified yet. Here, we show that dosage-sensitive ohnologs within CNV regions reported in patients with AD are related to the nervous system and are highly expressed in the brain, similar to other known susceptible genes for AD. We found that CNV regions in patients with AD contained dosage-sensitive genes, which are ohnologs not overlapping with control CNV regions, frequently. Furthermore, these dosage-sensitive genes in pathogenic CNV regions had a strong enrichment in the nervous system for mouse knockout phenotype and high expression in the brain similar to the known susceptible genes for AD. Our results demonstrated that selecting dosage-sensitive ohnologs out of multiple genes with pathogenic CNVs is effective in identifying the causative genes for AD. This methodology can be applied to other diseases caused by dosage imbalance and might help to establish the medical diagnosis by analysis of CNVs.

Segmental duplications: evolution and impact among the current Lepidoptera genomes

Background

Structural variation among genomes is now viewed to be as important as single nucleoid polymorphisms in influencing the phenotype and evolution of a species. Segmental duplication (SD) is defined as segments of DNA with homologous sequence.

Results

Here, we performed a systematic analysis of segmental duplications (SDs) among five lepidopteran reference genomes (Plutella xylostella, Danaus plexippus, Bombyx mori, Manduca sexta and Heliconius melpomene) to understand their potential impact on the evolution of these species. We find that the SDs content differed substantially among species, ranging from 1.2% of the genome in B. mori to 15.2% in H. melpomene. Most SDs formed very high identity (similarity higher than 90%) blocks but had very few large blocks. Comparative analysis showed that most of the SDs arose after the divergence of each linage and we found that P. xylostella and H. melpomene showed more duplications than other species, suggesting they might be able to tolerate extensive levels of variation in their genomes. Conserved ancestral and species specific SD events were assessed, revealing multiple examples of the gain, loss or maintenance of SDs over time. SDs content analysis showed that most of the genes embedded in SDs regions belonged to species-specific SDs (“Unique” SDs). Functional analysis of these genes suggested their potential roles in the lineage-specific evolution. SDs and flanking regions often contained transposable elements (TEs) and this association suggested some involvement in SDs formation. Further studies on comparison of gene expression level between SDs and non-SDs showed that the expression level of genes embedded in SDs was significantly lower, suggesting that structure changes in the genomes are involved in gene expression differences in species.

Conclusions

The results showed that most of the SDs were “unique SDs”, which originated after species formation. Functional analysis suggested that SDs might play different roles in different species. Our results provide a valuable resource beyond the genetic mutation to explore the genome structure for future Lepidoptera research.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-1007-y) contains supplementary material, which is available to authorized users.

Tandem duplications lead to novel expression patterns through exon shuffling in Drosophila yakuba

One common hypothesis to explain the impacts of tandem duplications is that whole gene duplications commonly produce additive changes in gene expression due to copy number changes. Here, we use genome wide RNA-seq data from a population sample of Drosophila yakuba to test this ‘gene dosage’ hypothesis. We observe little evidence of expression changes in response to whole transcript duplication capturing 5′ and 3′ UTRs. Among whole gene duplications, we observe evidence that dosage sharing across copies is likely to be common. The lack of expression changes after whole gene duplication suggests that the majority of genes are subject to tight regulatory control and therefore not sensitive to changes in gene copy number. Rather, we observe changes in expression level due to both shuffling of regulatory elements and the creation of chimeric structures via tandem duplication. Additionally, we observe 30 de novo gene structures arising from tandem duplications, 23 of which form with expression in the testes. Thus, the value of tandem duplications is likely to be more intricate than simple changes in gene dosage. The common regulatory effects from chimeric gene formation after tandem duplication may explain their contribution to genome evolution.

The enclosed work shows that whole gene duplications rarely affect gene expression, in contrast to widely held views that the adaptive value of duplicate genes is related to additive changes in gene expression due to gene copy number. We further explain how tandem duplications that create shuffled gene structures can force upregulation of gene sequences, de novo gene creation, and multifold changes in transcript levels. These results show that tandem duplications can produce new genes that are a source of immediate novelty associated with more extreme expression changes than previously suggested by theory. Further, these gene expression changes are a potential source of both beneficial and pathogenic mutations, immediately relevant to clinical and medical genetics in humans and other metazoans.

Characterization of Copy Number Variation's Potential Role in Marek's Disease

Marek’s Disease (MD) is a highly contagious pathogenic and oncogenic disease primarily affecting chickens. Chicken Lines 63 and 72, as well as their recombinant congenic strains (RCS) with varied susceptibility to MD, are ideal models to study the complex mechanisms of genetic resistance to MD. In this study, we investigated copy number variation (CNV) in these inbred chicken lines using the Affymetrix Axiom HD 600 K SNP genotyping array. We detected 393 CNV segments across all ten chicken lines, of which 12 CNVs were specifically identified in Line 72. We then assessed genetic structure based on CNV and observed markedly different patterns. Finally, we validated two deletion events in Line 72 and correlated them with genes expression using qPCR and RNA-seq, respectively. Our combined results indicated that these two CNV deletions were likely to contribute to MD susceptibility.

Inferring gene regulatory relationships with a high-dimensional robust approach

Gene expression (GE) levels have important biological and clinical implications. They are regulated by copy number alterations (CNAs). Modeling the regulatory relationships between GEs and CNAs facilitates understanding disease biology and can also have values in translational medicine. The expression level of a gene can be regulated by its cis-acting as well as trans-acting CNAs, and the set of trans-acting CNAs is usually not known, which poses a high-dimensional selection and estimation problem. Most of the existing studies share a common limitation in that they cannot accommodate long-tailed distributions or contamination of GE data. In this study, we develop a high-dimensional robust regression approach to infer the regulatory relationships between GEs and CNAs. A high-dimensional regression model is used to accommodate the effects of both cis-acting and trans-acting CNAs. A density power divergence loss function is used to accommodate long-tailed GE distributions and contamination. Penalization is adopted for regularized estimation and selection of relevant CNAs. The proposed approach is effectively realized using a coordinate descent algorithm. Simulation shows that it has competitive performance compared to the nonrobust benchmark and the robust LAD (least absolute deviation) approach. We analyze TCGA (The Cancer Genome Atlas) data on cutaneous melanoma and study GE-CNA regulations in the RAP (regulation of apoptosis) pathway, which further demonstrates the satisfactory performance of the proposed approach.

Diversity and regulatory impact of copy number variation in the primate Macaca fascicularis

BACKGROUND

Copy number variations (CNVs) are a significant source of genetic diversity and commonly found in mammalian genomes. We have generated a genome-wide CNV map for Cynomolgus monkeys (Macaca fascicularis). This crab-eating macaque is the closest animal model to humans that is used in biomedical research.

RESULTS

We show that Cynomolgus monkey CNVs are in general much smaller in size than gene loci and are specific to the population of origin. Genome-wide expression data from five vitally important organs demonstrates that CNVs in close proximity to transcription start sites associate strongly with expression changes. Among these eQTL genes we find an overrepresentation of genes involved in metabolism, receptor activity, and transcription.

CONCLUSION

These results provide evidence that CNVs shape tissue transcriptomes in monkey populations, potentially offering an adaptive advantage. We suggest that this genetic diversity should be taken into account when using Cynomolgus macaques as models.

Copy Number Variations in Tilapia Genomes

Discovering the nature and pattern of genome variation is fundamental in understanding phenotypic diversity among populations. Although several millions of single nucleotide polymorphisms (SNPs) have been discovered in tilapia, the genome-wide characterization of larger structural variants, such as copy number variation (CNV) regions has not been carried out yet. We conducted a genome-wide scan for CNVs in 47 individuals from three tilapia populations. Based on 254 Gb of high-quality paired-end sequencing reads, we identified 4642 distinct high-confidence CNVs. These CNVs account for 1.9% (12.411 Mb) of the used Nile tilapia reference genome. A total of 1100 predicted CNVs were found overlapping with exon regions of protein genes. Further association analysis based on linear model regression found 85 CNVs ranging between 300 and 27,000 base pairs significantly associated to population types (R ² > 0.9 and P > 0.001). Our study sheds first insights on genome-wide CNVs in tilapia. These CNVs among and within tilapia populations may have functional effects on phenotypes and specific adaptation to particular environments.

iGC-an integrated analysis package of gene expression and copy number alteration

BACKGROUND

With the advancement in high-throughput technologies, researchers can simultaneously investigate gene expression and copy number alteration (CNA) data from individual patients at a lower cost. Traditional analysis methods analyze each type of data individually and integrate their results using Venn diagrams. Challenges arise, however, when the results are irreproducible and inconsistent across multiple platforms. To address these issues, one possible approach is to concurrently analyze both gene expression profiling and CNAs in the same individual.

RESULTS

We have developed an open-source R/Bioconductor package (iGC). Multiple input formats are supported and users can define their own criteria for identifying differentially expressed genes driven by CNAs. The analysis of two real microarray datasets demonstrated that the CNA-driven genes identified by the iGC package showed significantly higher Pearson correlation coefficients with their gene expression levels and copy numbers than those genes located in a genomic region with CNA. Compared with the Venn diagram approach, the iGC package showed better performance.

CONCLUSION

The iGC package is effective and useful for identifying CNA-driven genes. By simultaneously considering both comparative genomic and transcriptomic data, it can provide better understanding of biological and medical questions. The iGC package's source code and manual are freely available at https://www.bioconductor.org/packages/release/bioc/html/iGC.html .

Evolutionary and Functional Features of Copy Number Variation in the Cattle Genome

Genomic structural variations are an important source of genetic diversity. Copy number variations (CNVs), gains and losses of large regions of genomic sequence between individuals of a species, have been associated with a wide variety of phenotypic traits. However, in cattle, as well as many other species, relatively little is understood about CNV, including frequency of CNVs in the genome, sizes, and locations, chromosomal properties, and evolutionary processes acting to shape CNV. In this work, we focused on copy number variation in the bovine genome, with the aim to detect CNVs in Bos taurus coding sequence and explore potential evolutionary mechanisms shaping these CNV. We identified and characterized CNV regions by utilizing exome sequence from 175 influential sires used in the Germplasm Evaluation project, representing 10 breeds. We examined various evolutionary and functional aspects of these CNVs, including selective constraint on CNV-overlapped genes, centrality of CNV genes in protein-protein interaction networks, and tissue-specific expression of CNV genes. Patterns of CNV in the Bos taurus genome reveal that reduced functional constraint and mutational bias may play a prominent role in shaping this type of structural variation.

Identification of a RAI1-associated disease network through integration of exome sequencing, transcriptomics, and 3D genomics

Background

Smith-Magenis syndrome (SMS) is a developmental disability/multiple congenital anomaly disorder resulting from haploinsufficiency of RAI1. It is characterized by distinctive facial features, brachydactyly, sleep disturbances, and stereotypic behaviors.

Methods

We investigated a cohort of 15 individuals with a clinical suspicion of SMS who showed neither deletion in the SMS critical region nor damaging variants in RAI1 using whole exome sequencing. A combination of network analysis (co-expression and biomedical text mining), transcriptomics, and circularized chromatin conformation capture (4C-seq) was applied to verify whether modified genes are part of the same disease network as known SMS-causing genes.

Results

Potentially deleterious variants were identified in nine of these individuals using whole-exome sequencing. Eight of these changes affect KMT2D, ZEB2, MAP2K2, GLDC, CASK, MECP2, KDM5C, and POGZ, known to be associated with Kabuki syndrome 1, Mowat-Wilson syndrome, cardiofaciocutaneous syndrome, glycine encephalopathy, mental retardation and microcephaly with pontine and cerebellar hypoplasia, X-linked mental retardation 13, X-linked mental retardation Claes-Jensen type, and White-Sutton syndrome, respectively. The ninth individual carries a de novo variant in JAKMIP1, a regulator of neuronal translation that was recently found deleted in a patient with autism spectrum disorder. Analyses of co-expression and biomedical text mining suggest that these pathologies and SMS are part of the same disease network. Further support for this hypothesis was obtained from transcriptome profiling that showed that the expression levels of both Zeb2 and Map2k2 are perturbed in Rai1^–/– mice. As an orthogonal approach to potentially contributory disease gene variants, we used chromatin conformation capture to reveal chromatin contacts between RAI1 and the loci flanking ZEB2 and GLDC, as well as between RAI1 and human orthologs of the genes that show perturbed expression in our Rai1^–/– mouse model.

Conclusions

These holistic studies of RAI1 and its interactions allow insights into SMS and other disorders associated with intellectual disability and behavioral abnormalities. Our findings support a pan-genomic approach to the molecular diagnosis of a distinctive disorder.

Electronic supplementary material

The online version of this article (doi:10.1186/s13073-016-0359-z) contains supplementary material, which is available to authorized users.

Evolutionary impact of transposable elements on genomic diversity and lineage-specific innovation in vertebrates

Since their discovery, a growing body of evidence has emerged demonstrating that transposable elements are important drivers of species diversity. These mobile elements exhibit a great variety in structure, size and mechanisms of transposition, making them important putative actors in organism evolution. The vertebrates represent a highly diverse and successful lineage that has adapted to a wide range of different environments. These animals also possess a rich repertoire of transposable elements, with highly diverse content between lineages and even between species. Here, we review how transposable elements are driving genomic diversity and lineage-specific innovation within vertebrates. We discuss the large differences in TE content between different vertebrate groups and then go on to look at how they affect organisms at a variety of levels: from the structure of chromosomes to their involvement in the regulation of gene expression, as well as in the formation and evolution of non-coding RNAs and protein-coding genes. In the process of doing this, we highlight how transposable elements have been involved in the evolution of some of the key innovations observed within the vertebrate lineage, driving the group's diversity and success.

BCL11A Haploinsufficiency Causes an Intellectual Disability Syndrome and Dysregulates Transcription

Intellectual disability (ID) is a common condition with considerable genetic heterogeneity. Next-generation sequencing of large cohorts has identified an increasing number of genes implicated in ID, but their roles in neurodevelopment remain largely unexplored. Here we report an ID syndrome caused by de novo heterozygous missense, nonsense, and frameshift mutations in BCL11A, encoding a transcription factor that is a putative member of the BAF swi/snf chromatin-remodeling complex. Using a comprehensive integrated approach to ID disease modeling, involving human cellular analyses coupled to mouse behavioral, neuroanatomical, and molecular phenotyping, we provide multiple lines of functional evidence for phenotypic effects. The etiological missense variants cluster in the amino-terminal region of human BCL11A, and we demonstrate that they all disrupt its localization, dimerization, and transcriptional regulatory activity, consistent with a loss of function. We show that Bcl11a haploinsufficiency in mice causes impaired cognition, abnormal social behavior, and microcephaly in accordance with the human phenotype. Furthermore, we identify shared aberrant transcriptional profiles in the cortex and hippocampus of these mouse models. Thus, our work implicates BCL11A haploinsufficiency in neurodevelopmental disorders and defines additional targets regulated by this gene, with broad relevance for our understanding of ID and related syndromes.

Copy number variations and stroke

Stroke is the third leading cause of death worldwide after heart disease and all forms of cancers. Monogenic disorders, genetic, and environmental risk factors contribute to damaging cerebral blood vessels and, consequently, cause stroke. Developments in genomic research led to the discovery of numerous copy number variants (CNVs) that have been recently identified as a new tool for understanding the genetic basis of many diseases. This review discusses the current understanding of the types of stroke, the existing knowledge on the involvement of specific CNVs in stroke as well as the limitations of the methods used for detecting CNVs like SNP-microarray. To confirm an unequivocally association between CNVs and stroke and extend the current findings, it would be desirable to use another methodology to detect smaller CNVs or CNVs in genomic regions poorly covered by this technique, for instance, CGH-array.

Targeted resequencing of regulatory regions at schizophrenia risk loci: Role of rare functional variants at chromatin repressive states

There is mounting evidence that regulatory variation plays an important role in genetic risk for schizophrenia. Here, we specifically search for regulatory variants at risk by sequencing promoter regions of twenty-three genes implied in schizophrenia by copy number variant or genome-wide association studies. After strict quality control, a total of 55,206bp per sample were analyzed in 526 schizophrenia cases and 516 controls from Galicia, NW Spain, using the Applied Biosystems SOLiD System. Variants were filtered based on frequency from public databases, chromatin states from the RoadMap Epigenomics Consortium at tissues relevant for schizophrenia, such as fetal brain, mid-frontal lobe, and angular gyrus, and prediction of functionality from RegulomeDB. The proportion of rare variants at polycomb repressive chromatin state at relevant tissues was higher in cases than in controls. The proportion of rare variants with predicted regulatory role was significantly higher in cases than in controls (P=0.0028, OR=1.93, 95% C.I.=1.23-3.04). Combination of information from both sources led to the identification of an excess of carriers of rare variants with predicted regulatory role located at polycomb repressive chromatin state at relevant tissues in cases versus controls (P=0.0016, OR=19.34, 95% C.I.=2.45-2495.26). The variants are located at two genes affected by the 17q12 copy number variant, LHX1 and HNF1B. These data strongly suggest that a specific epigenetic mechanism, chromatin remodeling by histone modification during early development, may be impaired in a subset of schizophrenia patients, in agreement with previous data.

Genome-Wide Detection of CNVs and Their Association with Meat Tenderness in Nelore Cattle

Brazil is one of the largest beef producers and exporters in the world with the Nelore breed representing the vast majority of Brazilian cattle (Bos taurus indicus). Despite the great adaptability of the Nelore breed to tropical climate, meat tenderness (MT) remains to be improved. Several factors including genetic composition can influence MT. In this article, we report a genome-wide analysis of copy number variation (CNV) inferred from Illumina® High Density SNP-chip data for a Nelore population of 723 males. We detected >2,600 CNV regions (CNVRs) representing ≈6.5% of the genome. Comparing our results with previous studies revealed an overlap in ≈1400 CNVRs (>50%). A total of 1,155 CNVRs (43.6%) overlapped 2,750 genes. They were enriched for processes involving guanosine triphosphate (GTP), previously reported to influence skeletal muscle physiology and morphology. Nelore CNVRs also overlapped QTLs for MT reported in other breeds (8.9%, 236 CNVRs) and from a previous study with this population (4.1%, 109 CNVRs). Two CNVRs were also proximal to glutathione metabolism genes that were previously associated with MT. Genome-wide association study of CN state with estimated breeding values derived from meat shear force identified 6 regions, including a region on BTA3 that contains genes of the cAMP and cGMP pathway. Ten CNVRs that overlapped regions associated with MT were successfully validated by qPCR. Our results represent the first comprehensive CNV study in Bos taurus indicus cattle and identify regions in which copy number changes are potentially of importance for the MT phenotype.

Association of a Chromosomal Rearrangement Event with Mouse Posterior Polymorphous Corneal Dystrophy and Alterations in Csrp2bp, Dzank1, and Ovol2 Gene Expression

We have previously described a mouse model of human posterior polymorphous corneal dystrophy (PPCD) and localized the causative mutation to a 6.2 Mbp region of chromosome 2, termed Ppcd1. We now show that the gene rearrangement linked to mouse Ppcd1 is a 3.9 Mbp chromosomal inversion flanked by 81 Kbp and 542 bp deletions. This recombination event leads to deletion of Csrp2bp Exons 8 through 11, Dzank1 Exons 20 and 21, and the pseudogene Znf133. In addition, we identified translocation of novel downstream sequences to positions adjacent to Csrp2bp Exon 7 and Dzank1 Exon 20. Twelve novel fusion transcripts involving Csrp2bp or Dzank1 linked to downstream sequences have been identified. Eight are expressed at detectable levels in PPCD1 but not wildtype eyes. Upregulation of two Csrp2bp fusion transcripts, as well as upregulation of the adjacent gene, Ovol2, was observed. Absence of the PPCD1 phenotype in animals haploinsufficient for Csrp2bp or both Csrp2bp and Dzank1 rules out haploinsufficiency of these genes as a cause of mouse PPCD1. Complementation experiments confirm that PPCD1 embryonic lethality is due to disruption of Csrp2bp expression. The ocular expression pattern of Csrp2bp is consistent with a role for this protein in corneal development and pathogenesis of PPCD1.

Genome-wide CNV analysis reveals variants associated with growth traits in Bos indicus

BACKGROUND

Apart from single nucleotide polymorphism (SNP), copy number variation (CNV) is another important type of genetic variation, which may affect growth traits and play key roles for the production of beef cattle. To date, no genome-wide association study (GWAS) for CNV and body traits in beef cattle has been reported, so the present study aimed to investigate this type of association in one of the most important cattle subspecies: Bos indicus (Nellore breed).

RESULTS

We have used intensity data from over 700,000 SNP probes across the bovine genome to detect common CNVs in a sample of 2230 Nellore cattle, and performed GWAS between the detected CNVs and nine growth traits. After filtering for frequency and length, a total of 231 CNVs ranging from 894 bp to 4,855,088 bp were kept and tested as predictors for each growth trait using linear regression analysis with principal components correction. There were 49 significant associations identified among 17 CNVs and seven body traits after false discovery rate correction (P < 0.05). Among the 17 CNVs, three were significant or marginally significant for all the traits. We have compared the locations of associated CNVs with quantitative trait locus and the RefGene database, and found two sets of 9 CNVs overlapping with either known QTLs or genes, respectively. The gene overlapping with CNV100, KCNJ12, is a functional candidate for muscle development and plays critical roles in muscling traits.

CONCLUSION

This study presents the first CNV-based GWAS of growth traits using high density SNP microarray data in cattle. We detected 17 CNVs significantly associated with seven growth traits and one of them (CNV100) may be involved in growth traits through KCNJ12.

Integrated Genomics for Pinpointing Survival Loci within Arm-Level Somatic Copy Number Alterations

The identification of driver loci underlying arm-level somatic copy number alterations (SCNAs) in cancer has remained challenging and incomplete. Here, we assess the relative impact and present a detailed landscape of arm-level SCNAs in 10,985 patient samples across 33 cancer types from The Cancer Genome Atlas (TCGA). Furthermore, using chromosome 9p loss in lower grade glioma (LGG) as a model, we employ a unique multi-tiered genomic dissection strategy using 540 patients from three independent LGG datasets to identify genetic loci that govern tumor aggressiveness and poor survival. This comprehensive approach uncovered several 9p loss-specific prognostic markers, validated existing ones, and redefined the impact of CDKN2A loss in LGG.

Mitochondrial Dysfunction may explain symptom variation in Phelan-McDermid Syndrome

Phelan-McDermid Syndrome (PMS), which is defined by a deletion within 22q13, demonstrates significant phenotypic variation. Given that six mitochondrial genes are located within 22q13, including complex I and IV genes, we hypothesize that mitochondrial complex activity abnormalities may explain phenotypic variation in PMS symptoms. Complex I, II, II + III and IV activity was measured in 51 PMS participants. Caretakers completed questionnaires and provided genetic information through the PMS foundation registry. Complex activity was abnormal in 59% of PMS participants. Abnormalities were found in complex I and IV but not complex II + III and II activity, consistent with disruption of genes within the 22q13 region. However, complex activity abnormalities were not related to specific gene deletions suggesting a "neighboring effect" of regional deletions on adjacent gene expression. A specific combination of symptoms (autism spectrum disorder, developmental regression, failure-to-thrive, exercise intolerance/fatigue) was associated with complex activity abnormalities. 64% of 106 individuals in the PMS foundation registry who did not have complex activity measured also endorsed this pattern of symptoms. These data suggest that mitochondrial abnormalities, specifically abnormalities in complex I and IV activity, may explain some phenotypic variation in PMS individuals. These results point to novel pathophysiology mechanisms and treatment targets for PMS patients.

Genome-wide identification of copy number variation using high-density single-nucleotide polymorphism array in Japanese Black cattle

Background

Copy number variation (CNV) is an important source of genetic variability associated with phenotypic variation and disease susceptibility. Comprehensive genome-wide CNV maps provide valuable information for genetic and functional studies. To identify CNV in Japanese Black cattle, we performed a genome-wide autosomal screen using genomic data from 1,481 animals analyzed with the Illumina Bovine High-Density (HD) BeadChip Array (735,293 single-nucleotide polymorphisms (SNPs) with an average marker interval of 3.4 kb on the autosomes).

Results

We identified a total of 861 CNV regions (CNVRs) across all autosomes, which covered 43.65 Mb of the UMD3.1 genome assembly and corresponded to 1.74 % of the 29 bovine autosomes. Overall, 35 % of the CNVRs were present at a frequency of > 1 % in 1,481 animals. The estimated lengths of CNVRs ranged from 1.1 kb to 1.4 Mb, with an average of 50.7 kb. The average number of CNVR events per animal was 35. Comparisons with previously reported cattle CNV showed that 72 % of the CNVR calls detected in this study were within or overlapped with known CNVRs. Experimentally, three CNVRs were validated using quantitative PCR, and one CNVR was validated using PCR with flanking primers for the deleted region. Out of the 861 CNVRs, 390 contained 717 Ensembl-annotated genes significantly enriched for stimulus response, cellular defense response, and immune response in the Gene Ontology (GO) database. To associate genes contained in CNVRs with phenotypes, we converted 560 bovine Ensembl gene IDs to their 438 orthologous associated mouse gene IDs, and 195 of these mouse orthologous genes were categorized into 1,627 phenotypes in the Mouse Genome Informatics (MGI) database.

Conclusions

We identified 861 CNVRs in 1,481 Japanese Black cattle using the Illumina BovineHD BeadChip Array. The genes contained in CNVRs were characterized using GO analysis and the mouse orthologous genes were characterized using the MGI database. The comprehensive genome-wide CNVRs map will facilitate identification of genetic variation and disease-susceptibility alleles in Japanese Black cattle.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-016-0335-z) contains supplementary material, which is available to authorized users.

Different Facets of Copy Number Changes: Permanent, Transient, and Adaptive

Chromosomal copy number changes are frequently associated with harmful consequences and are thought of as an underlying mechanism for the development of diseases. However, changes in copy number are observed during development and occur during normal biological processes. In this review, we highlight the causes and consequences of copy number changes in normal physiologic processes as well as cover their associations with cancer and acquired drug resistance. We discuss the permanent and transient nature of copy number gains and relate these observations to a new mechanism driving transient site-specific copy gains (TSSGs). Finally, we discuss implications of TSSGs in generating intratumoral heterogeneity and tumor evolution and how TSSGs can influence the therapeutic response in cancer.

The Mouse Universal Genotyping Array: From Substrains to Subspecies

Genotyping microarrays are an important resource for genetic mapping, population genetics, and monitoring of the genetic integrity of laboratory stocks. We have developed the third generation of the Mouse Universal Genotyping Array (MUGA) series, GigaMUGA, a 143,259-probe Illumina Infinium II array for the house mouse (Mus musculus). The bulk of the content of GigaMUGA is optimized for genetic mapping in the Collaborative Cross and Diversity Outbred populations, and for substrain-level identification of laboratory mice. In addition to 141,090 single nucleotide polymorphism probes, GigaMUGA contains 2006 probes for copy number concentrated in structurally polymorphic regions of the mouse genome. The performance of the array is characterized in a set of 500 high-quality reference samples spanning laboratory inbred strains, recombinant inbred lines, outbred stocks, and wild-caught mice. GigaMUGA is highly informative across a wide range of genetically diverse samples, from laboratory substrains to other Mus species. In addition to describing the content and performance of the array, we provide detailed probe-level annotation and recommendations for quality control.

A relay race on the evolutionary adaptation spectrum

Adaptation is the process in which organisms improve their fitness by changing their phenotype using genetic or non-genetic mechanisms. The adaptation toolbox consists of varied molecular and genetic means that we posit span an almost continuous "adaptation spectrum." Different adaptations are characterized by the time needed for organisms to attain them and by their duration. We suggest that organisms often adapt by progressing the adaptation spectrum, starting with rapidly attained physiological and epigenetic adaptations and culminating with slower long-lasting genetic ones. A tantalizing possibility is that earlier adaptations facilitate realization of later ones.

Integrating genetics and epigenetics in breast cancer: biological insights, experimental, computational methods and therapeutic potential

BACKGROUND

Development of human cancer can proceed through the accumulation of different genetic changes affecting the structure and function of the genome. Combined analyses of molecular data at multiple levels, such as DNA copy-number alteration, mRNA and miRNA expression, can clarify biological functions and pathways deregulated in cancer. The integrative methods that are used to investigate these data involve different fields, including biology, bioinformatics, and statistics.

RESULTS

These methodologies are presented in this review, and their implementation in breast cancer is discussed with a focus on integration strategies. We report current applications, recent studies and interesting results leading to the identification of candidate biomarkers for diagnosis, prognosis, and therapy in breast cancer by using both individual and combined analyses.

CONCLUSION

This review presents a state of art of the role of different technologies in breast cancer based on the integration of genetics and epigenetics, and shares some issues related to the new opportunities and challenges offered by the application of such integrative approaches.

Mechanisms of Origin, Phenotypic Effects and Diagnostic Implications of Complex Chromosome Rearrangements

Complex chromosome rearrangements (CCRs) are currently defined as structural genome variations that involve more than 2 chromosome breaks and result in exchanges of chromosomal segments. They are thought to be extremely rare, but their detection rate is rising because of improvements in molecular cytogenetic technology. Their population frequency is also underestimated, since many CCRs may not elicit a phenotypic effect. CCRs may be the result of fork stalling and template switching, microhomology-mediated break-induced repair, breakage-fusion-bridge cycles, or chromothripsis. Patients with chromosomal instability syndromes show elevated rates of CCRs due to impaired DNA double-strand break responses during meiosis. Therefore, the putative functions of the proteins encoded by ATM, BLM, WRN, ATR, MRE11, NBS1, and RAD51 in preventing CCRs are discussed. CCRs may exert a pathogenic effect by either (1) gene dosage-dependent mechanisms, e.g. haploinsufficiency, (2) mechanisms based on disruption of the genomic architecture, such that genes, parts of genes or regulatory elements are truncated, fused or relocated and thus their interactions disturbed - these mechanisms will predominantly affect gene expression - or (3) mixed mutation mechanisms in which a CCR on one chromosome is combined with a different type of mutation on the other chromosome. Such inferred mechanisms of pathogenicity need corroboration by mRNA sequencing. Also, future studies with in vitro models, such as inducible pluripotent stem cells from patients with CCRs, and transgenic model organisms should substantiate current inferences regarding putative pathogenic effects of CCRs. The ramifications of the growing body of information on CCRs for clinical and experimental genetics and future treatment modalities are briefly illustrated with 2 cases, one of which suggests KDM4C (JMJD2C) as a novel candidate gene for mental retardation.

Divergence patterns of genic copy number variation in natural populations of the house mouse (Mus musculus domesticus) reveal three conserved genes with major population-specific expansions

Copy number variation represents a major source of genetic divergence, yet the evolutionary dynamics of genic copy number variation in natural populations during differentiation and adaptation remain unclear. We applied a read depth approach to genome resequencing data to detect copy number variants (CNVs) ≥1 kb in wild-caught mice belonging to four populations of Mus musculus domesticus. We complemented the bioinformatics analyses with experimental validation using droplet digital PCR. The specific focus of our analysis is CNVs that include complete genes, as these CNVs could be expected to contribute most directly to evolutionary divergence. In total, 1863 transcription units appear to be completely encompassed within CNVs in at least one individual when compared to the reference assembly. Further, 179 of these CNVs show population-specific copy number differences, and 325 are subject to complete deletion in multiple individuals. Among the most copy-number variable genes are three highly conserved genes that encode the splicing factor CWC22, the spindle protein SFI1, and the Holliday junction recognition protein HJURP. These genes exhibit population-specific expansion patterns that suggest involvement in local adaptations. We found that genes that overlap with large segmental duplications are generally more copy-number variable. These genes encode proteins that are relevant for environmental and behavioral interactions, such as vomeronasal and olfactory receptors, as well as major urinary proteins and several proteins of unknown function. The overall analysis shows that genic CNVs contribute more to population differentiation in mice than in humans and may promote and speed up population divergence.

Genomic copy number variation in Mus musculus

BACKGROUND

Copy number variation is an important dimension of genetic diversity and has implications in development and disease. As an important model organism, the mouse is a prime candidate for copy number variant (CNV) characterization, but this has yet to be completed for a large sample size. Here we report CNV analysis of publicly available, high-density microarray data files for 351 mouse tail samples, including 290 mice that had not been characterized for CNVs previously.

RESULTS

We found 9634 putative autosomal CNVs across the samples affecting 6.87% of the mouse reference genome. We find significant differences in the degree of CNV uniqueness (single sample occurrence) and the nature of CNV-gene overlap between wild-caught mice and classical laboratory strains. CNV-gene overlap was associated with lipid metabolism, pheromone response and olfaction compared to immunity, carbohydrate metabolism and amino-acid metabolism for wild-caught mice and classical laboratory strains, respectively. Using two subspecies of wild-caught Mus musculus, we identified putative CNVs unique to those subspecies and show this diversity is better captured by wild-derived laboratory strains than by the classical laboratory strains. A total of 9 genic copy number variable regions (CNVRs) were selected for experimental confirmation by droplet digital PCR (ddPCR).

CONCLUSION

The analysis we present is a comprehensive, genome-wide analysis of CNVs in Mus musculus, which increases the number of known variants in the species and will accelerate the identification of novel variants in future studies.

Inside the 8p23.1 duplication syndrome; eight microduplications of likely or uncertain clinical significance

The 8p23.1 duplication syndrome (8p23.1 DS) is a recurrent genomic condition with an estimated prevalence of 1 in 58,000. The core 3.68 Mb duplication contains 32 genes of which five are currently candidates for the phenotypic features. Here we describe four patients and five families with eight microduplications of 8p23.1 ranging from 187 to 1082 kb in size and one atypical duplication of 4 Mb. These indicate that a minimal region of overlap (MRO) in medial 8p23.1 can give rise to features of 8p23.1 DS including developmental delay, dysmorphism, macrocephaly and otitis media, but not congenital heart disease (CHD). This MRO spans 776 kb (chr8:10,167,881-10,943,836 hg19) and contains SOX7 and seven of the other 32 core 8p23.1 DS genes. In centromeric 8p23.1, microduplications including GATA4 can give rise to non-syndromic CHD but the clinical significance of two smaller centromeric microduplications without GATA4 was uncertain due to severe neurological profiles not usually found in 8p23.1 DS. The clinical significance of three further 8p23.1 microduplications was uncertain due to additional genetic factors without which the probands might not have come to medical attention. Variable expressivity was indicated by the almost entirely unaffected parents in all five families and the mildly affected sibling in one. Intronic interruptions of six genes by microduplication breakpoint intervals had no apparent additional clinical consequences. Our results suggest that 8p23.1 DS is an oligogenetic condition largely caused by the duplication and interactions of the SOX7 and GATA4 transcription factors.

Copy number variations and cognitive phenotypes in unselected populations

IMPORTANCE

The association of copy number variations (CNVs), differing numbers of copies of genetic sequence at locations in the genome, with phenotypes such as intellectual disability has been almost exclusively evaluated using clinically ascertained cohorts. The contribution of these genetic variants to cognitive phenotypes in the general population remains unclear.

OBJECTIVE

To investigate the clinical features conferred by CNVs associated with known syndromes in adult carriers without clinical preselection and to assess the genome-wide consequences of rare CNVs (frequency ≤0.05%; size ≥250 kilobase pairs [kb]) on carriers' educational attainment and intellectual disability prevalence in the general population.

DESIGN, SETTING, AND PARTICIPANTS

The population biobank of Estonia contains 52,000 participants enrolled from 2002 through 2010. General practitioners examined participants and filled out a questionnaire of health- and lifestyle-related questions, as well as reported diagnoses. Copy number variant analysis was conducted on a random sample of 7877 individuals and genotype-phenotype associations with education and disease traits were evaluated. Our results were replicated on a high-functioning group of 993 Estonians and 3 geographically distinct populations in the United Kingdom, the United States, and Italy.

MAIN OUTCOMES AND MEASURES

Phenotypes of genomic disorders in the general population, prevalence of autosomal CNVs, and association of these variants with educational attainment (from less than primary school through scientific degree) and prevalence of intellectual disability.

RESULTS

Of the 7877 in the Estonian cohort, we identified 56 carriers of CNVs associated with known syndromes. Their phenotypes, including cognitive and psychiatric problems, epilepsy, neuropathies, obesity, and congenital malformations are similar to those described for carriers of identical rearrangements ascertained in clinical cohorts. A genome-wide evaluation of rare autosomal CNVs (frequency, ≤0.05%; ≥250 kb) identified 831 carriers (10.5%) of the screened general population. Eleven of 216 (5.1%) carriers of a deletion of at least 250 kb (odds ratio [OR], 3.16; 95% CI, 1.51-5.98; P = 1.5e-03) and 6 of 102 (5.9%) carriers of a duplication of at least 1 Mb (OR, 3.67; 95% CI, 1.29-8.54; P = .008) had an intellectual disability compared with 114 of 6819 (1.7%) in the Estonian cohort. The mean education attainment was 3.81 (P = 1.06e-04) among 248 (≥250 kb) deletion carriers and 3.69 (P = 5.024e-05) among 115 duplication carriers (≥1 Mb). Of the deletion carriers, 33.5% did not graduate from high school (OR, 1.48; 95% CI, 1.12-1.95; P = .005) and 39.1% of duplication carriers did not graduate high school (OR, 1.89; 95% CI, 1.27-2.8; P = 1.6e-03). Evidence for an association between rare CNVs and lower educational attainment was supported by analyses of cohorts of adults from Italy and the United States and adolescents from the United Kingdom.

CONCLUSIONS AND RELEVANCE

Known pathogenic CNVs in unselected, but assumed to be healthy, adult populations may be associated with unrecognized clinical sequelae. Additionally, individually rare but collectively common intermediate-size CNVs may be negatively associated with educational attainment. Replication of these findings in additional population groups is warranted given the potential implications of this observation for genomics research, clinical care, and public health.

Connecting Anxiety and Genomic Copy Number Variation: A Genome-Wide Analysis in CD-1 Mice

Genomic copy number variants (CNVs) have been implicated in multiple psychiatric disorders, but not much is known about their influence on anxiety disorders specifically. Using next-generation sequencing (NGS) and two additional array-based genotyping approaches, we detected CNVs in a mouse model consisting of two inbred mouse lines showing high (HAB) and low (LAB) anxiety-related behavior, respectively. An influence of CNVs on gene expression in the central (CeA) and basolateral (BLA) amygdala, paraventricular nucleus (PVN), and cingulate cortex (Cg) was shown by a two-proportion Z-test (p = 1.6 x 10^-31), with a positive correlation in the CeA (p = 0.0062), PVN (p = 0.0046) and Cg (p = 0.0114), indicating a contribution of CNVs to the genetic predisposition to trait anxiety in the specific context of HAB/LAB mice. In order to confirm anxiety-relevant CNVs and corresponding genes in a second mouse model, we further examined CD-1 outbred mice. We revealed the distribution of CNVs by genotyping 64 CD 1 individuals using a high-density genotyping array (Jackson Laboratory). 78 genes within those CNVs were identified to show nominally significant association (48 genes), or a statistical trend in their association (30 genes) with the time animals spent on the open arms of the elevated plus-maze (EPM). Fifteen of them were considered promising candidate genes of anxiety-related behavior as we could show a significant overlap (permutation test, p = 0.0051) with genes within HAB/LAB CNVs. Thus, here we provide what is to our knowledge the first extensive catalogue of CNVs in CD-1 mice and potential corresponding candidate genes linked to anxiety-related behavior in mice.

Adaptive potential of genomic structural variation in human and mammalian evolution

Because phenotypic innovations must be genetically heritable for biological evolution to proceed, it is natural to consider new mutation events as well as standing genetic variation as sources for their birth. Previous research has identified a number of single-nucleotide polymorphisms that underlie a subset of adaptive traits in organisms. However, another well-known class of variation, genomic structural variation, could have even greater potential to produce adaptive phenotypes, due to the variety of possible types of alterations (deletions, insertions, duplications, among others) at different genomic positions and with variable lengths. It is from these dramatic genomic alterations, and selection on their phenotypic consequences, that adaptations leading to biological diversification could be derived. In this review, using studies in humans and other mammals, we highlight examples of how phenotypic variation from structural variants might become adaptive in populations and potentially enable biological diversification. Phenotypic change arising from structural variants will be described according to their immediate effect on organismal metabolic processes, immunological response and physical features. Study of population dynamics of segregating structural variation can therefore provide a window into understanding current and historical biological diversification.

Human genomics. Effect of predicted protein-truncating genetic variants on the human transcriptome

Accurate prediction of the functional effect of genetic variation is critical for clinical genome interpretation. We systematically characterized the transcriptome effects of protein-truncating variants, a class of variants expected to have profound effects on gene function, using data from the Genotype-Tissue Expression (GTEx) and Geuvadis projects. We quantitated tissue-specific and positional effects on nonsense-mediated transcript decay and present an improved predictive model for this decay. We directly measured the effect of variants both proximal and distal to splice junctions. Furthermore, we found that robustness to heterozygous gene inactivation is not due to dosage compensation. Our results illustrate the value of transcriptome data in the functional interpretation of genetic variants.

A Potential Contributory Role for Ciliary Dysfunction in the 16p11.2 600 kb BP4-BP5 Pathology

The 16p11.2 600 kb copy-number variants (CNVs) are associated with mirror phenotypes on BMI, head circumference, and brain volume and represent frequent genetic lesions in autism spectrum disorders (ASDs) and schizophrenia. Here we interrogated the transcriptome of individuals carrying reciprocal 16p11.2 CNVs. Transcript perturbations correlated with clinical endophenotypes and were enriched for genes associated with ASDs, abnormalities of head size, and ciliopathies. Ciliary gene expression was also perturbed in orthologous mouse models, raising the possibility that ciliary dysfunction contributes to 16p11.2 pathologies. In support of this hypothesis, we found structural ciliary defects in the CA1 hippocampal region of 16p11.2 duplication mice. Moreover, by using an established zebrafish model, we show genetic interaction between KCTD13, a key driver of the mirrored neuroanatomical phenotypes of the 16p11.2 CNV, and ciliopathy-associated genes. Overexpression of BBS7 rescues head size and neuroanatomical defects of kctd13 morphants, whereas suppression or overexpression of CEP290 rescues phenotypes induced by KCTD13 under- or overexpression, respectively. Our data suggest that dysregulation of ciliopathy genes contributes to the clinical phenotypes of these CNVs.

The variation game: Cracking complex genetic disorders with NGS and omics data

Tremendous advances in Next Generation Sequencing (NGS) and high-throughput omics methods have brought us one step closer towards mechanistic understanding of the complex disease at the molecular level. In this review, we discuss four basic regulatory mechanisms implicated in complex genetic diseases, such as cancer, neurological disorders, heart disease, diabetes, and many others. The mechanisms, including genetic variations, copy-number variations, posttranscriptional variations, and epigenetic variations, can be detected using a variety of NGS methods. We propose that malfunctions detected in these mechanisms are not necessarily independent, since these malfunctions are often found associated with the same disease and targeting the same gene, group of genes, or functional pathway. As an example, we discuss possible rewiring effects of the cancer-associated genetic, structural, and posttranscriptional variations on the protein-protein interaction (PPI) network centered around P53 protein. The review highlights multi-layered complexity of common genetic disorders and suggests that integration of NGS and omics data is a critical step in developing new computational methods capable of deciphering this complexity.

The impact of human copy number variation on gene expression

Recent years have witnessed a flurry of important technological and methodological developments in the discovery and analysis of copy number variations (CNVs), which are increasingly enabling the systematic evaluation of their impact on a broad range of phenotypes from molecular-level (intermediate) traits to higher-order clinical phenotypes. Like single nucleotide variants in the human genome, CNVs have been linked to complex traits in humans, including disease and drug response. These recent developments underscore the importance of incorporating complex forms of genetic variation into disease mapping studies and promise to transform our understanding of genome function and the genetic basis of disease. Here we review some of the findings that have emerged from transcriptome studies of CNVs facilitated by the rapid advances in -omics technologies and corresponding methodologies.

Haplotype phasing and inheritance of copy number variants in nuclear families

DNA copy number variants (CNVs) that alter the copy number of a particular DNA segment in the genome play an important role in human phenotypic variability and disease susceptibility. A number of CNVs overlapping with genes have been shown to confer risk to a variety of human diseases thus highlighting the relevance of addressing the variability of CNVs at a higher resolution. So far, it has not been possible to deterministically infer the allelic composition of different haplotypes present within the CNV regions. We have developed a novel computational method, called PiCNV, which enables to resolve the haplotype sequence composition within CNV regions in nuclear families based on SNP genotyping microarray data. The algorithm allows to i) phase normal and CNV-carrying haplotypes in the copy number variable regions, ii) resolve the allelic copies of rearranged DNA sequence within the haplotypes and iii) infer the heritability of identified haplotypes in trios or larger nuclear families. To our knowledge this is the first program available that can deterministically phase null, mono-, di-, tri- and tetraploid genotypes in CNV loci. We applied our method to study the composition and inheritance of haplotypes in CNV regions of 30 HapMap Yoruban trios and 34 Estonian families. For 93.6% of the CNV loci, PiCNV enabled to unambiguously phase normal and CNV-carrying haplotypes and follow their transmission in the corresponding families. Furthermore, allelic composition analysis identified the co-occurrence of alternative allelic copies within 66.7% of haplotypes carrying copy number gains. We also observed less frequent transmission of CNV-carrying haplotypes from parents to children compared to normal haplotypes and identified an emergence of several de novo deletions and duplications in the offspring.

Quantitative effect of a CNV on a morphological trait in chickens

Copy Number Variation has been associated with morphological traits, developmental defects or disease susceptibility. The autosomal dominant Pea-comb mutation in chickens is due to the massive amplification of a CNV in intron 1 of SOX5 and provides a unique opportunity to assess the effect of variation in the number of repeats on quantitative traits such as comb size and comb mass in Pea-comb chickens. The quantitative variation of comb size was estimated by 2D morphometry and the number of repeats (RQ) was estimated by qPCR, in a total of 178 chickens from 3 experimental lines, two of them showing segregation for the Pea-comb mutation. This study included only Pea-comb chickens. Analysis of variance showed highly significant effects of line and sex on comb measurements. Adult body weight (BW) and RQ were handled as covariates. BW significantly influenced comb mass but not comb size. RQ values significantly influenced comb size, and the linear regression coefficient was highest for heterozygous carriers: the higher the number of repeats, the smaller the comb size. A similar trend was observed for comb mass. The CNV contributed to 3.4% of the phenotypic variance of comb size in heterozygous carriers of the CNV, an order of magnitude frequently encountered for QTLs. Surprisingly, there was no such relationship between RQ values and comb size in the homozygous line. It may be concluded that heterozygosity for a CNV in a non-coding region may contribute to phenotypic plasticity.

Multiple system atrophy: the application of genetics in understanding etiology

Classically defined phenotypically by a triad of cerebellar ataxia, parkinsonism, and autonomic dysfunction in conjunction with pyramidal signs, multiple system atrophy (MSA) is a rare and progressive neurodegenerative disease affecting an estimated 3-4 per every 100,000 individuals among adults 50-99 years of age. With a pathological hallmark of alpha-synuclein-immunoreactive glial cytoplasmic inclusions (GCIs; Papp-Lantos inclusions), MSA patients exhibit marked neurodegenerative changes in the striatonigral and/or olivopontocerebellar structures of the brain. As a member of the alpha-synucleinopathy family, which is defined by its well-demarcated alpha-synuclein-immunoreactive inclusions and aggregation, MSA's clinical presentation exhibits several overlapping features with other members including Parkinson's disease (PD) and dementia with Lewy bodies (DLB). Given the extensive fund of knowledge regarding the genetic etiology of PD revealed within the past several years, a genetic investigation of MSA is warranted. While a current genome-wide association study is underway for MSA to further clarify the role of associated genetic loci and single-nucleotide polymorphisms, several cases have presented solid preliminary evidence of a genetic etiology. Naturally, genes and variants manifesting known associations with PD (and other phenotypically similar neurodegenerative disorders), including SNCA and MAPT, have been comprehensively investigated in MSA patient cohorts. More recently variants in COQ2 have been linked to MSA in the Japanese population although this finding awaits replication. Nonetheless, significant positive associations with subsequent independent replication studies have been scarce. With very limited information regarding genetic mutations or alterations in gene dosage as a cause of MSA, the search for novel risk genes, which may be in the form of common variants or rare variants, is the logical nexus for MSA research. We believe that the application of next generation genetic methods to MSA will provide valuable insight into the underlying causes of this disease, and will be central to the identification of etiologic-based therapies.

Mechanisms of x chromosome dosage compensation

In many animals, males have one X and females have two X chromosomes. The difference in X chromosome dosage between the two sexes is compensated by mechanisms that regulate X chromosome transcription. Recent advances in genomic techniques have provided new insights into the molecular mechanisms of X chromosome dosage compensation. In this review, I summarize our current understanding of dosage imbalance in general, and then review the molecular mechanisms of X chromosome dosage compensation with an emphasis on the parallels and differences between the three well-studied model systems, M. musculus, D. melanogaster and C. elegans.

Deciphering the pathogenic consequences of chromosomal aberrations in human genetic disease

Chromosomal aberrations include translocations, deletions, duplications, inversions, aneuploidies and complex rearrangements. They underlie genetic disease in roughly 15% of patients with multiple congenital abnormalities and/or mental retardation (MCA/MR). In genetic diagnostics, the pathogenicity of chromosomal aberrations in these patients is typically assessed based on criteria such as phenotypic similarity to other patients with the same or overlapping aberration, absence in healthy individuals, de novo occurrence, and protein coding gene content. However, a thorough understanding of the molecular mechanisms that lead to MCA/MR as a result of chromosome aberrations is often lacking. Chromosome aberrations can affect one or more genes in a complex manner, such as by changing the regulation of gene expression, by disrupting exons, and by creating fusion genes. The precise delineation of breakpoints by whole-genome sequencing enables the construction of local genomic architecture and facilitates the prediction of the molecular determinants of the patient's phenotype. Here, we review current methods for breakpoint identification and their impact on the interpretation of chromosome aberrations in patients with MCA/MR. In addition, we discuss opportunities to dissect disease mechanisms based on large-scale genomic technologies and studies in model organisms.

The synapse in schizophrenia

It has been several decades since synaptic dysfunction was first suggested to play a role in schizophrenia, but only in the last few years has convincing evidence been obtained as progress has been made in elucidating the genetic underpinnings of the disorder. In the intervening years much has been learned concerning the complex macromolecular structure of the synapse itself, and genetic studies are now beginning to draw upon these advances. Here we outline our current understanding of the genetic architecture of schizophrenia and examine the evidence for synaptic involvement. A strong case can now be made that disruption of glutamatergic signalling pathways regulating synaptic plasticity contributes to the aetiology of schizophrenia.

Expansion of stochastic expression repertoire by tandem duplication in mouse Protocadherin-α cluster

Tandem duplications are concentrated within the Pcdh cluster throughout vertebrate evolution and as copy number variations (CNVs) in human populations, but the effects of tandem duplication in the Pcdh cluster remain elusive. To investigate the effects of tandem duplication in the Pcdh cluster, here we generated and analyzed a new line of the Pcdh cluster mutant mice. In the mutant allele, a 218-kb region containing the Pcdh-α2 to Pcdh-αc2 variable exons with their promoters was duplicated and the individual duplicated Pcdh isoforms can be disctinguished. The individual duplicated Pcdh-α isoforms showed diverse expression level with stochastic expression manner, even though those have an identical promoter sequence. Interestingly, the 5'-located duplicated Pcdh-αc2, which is constitutively expressed in the wild-type brain, shifted to stochastic expression accompanied by increased DNA methylation. These results demonstrate that tandem duplication in the Pcdh cluster expands the stochastic expression repertoire irrespective of sequence divergence.

Copy number deletion has little impact on gene expression levels in racehorses

Copy number variations (CNVs), important genetic factors for study of human diseases, may have as large of an effect on phenotype as do single nucleotide polymorphisms. Indeed, it is widely accepted that CNVs are associated with differential disease susceptibility. However, the relationships between CNVs and gene expression have not been characterized in the horse. In this study, we investigated the effects of copy number deletion in the blood and muscle transcriptomes of Thoroughbred racing horses. We identified a total of 1,246 CNVs of deletion polymorphisms using DNA re-sequencing data from 18 Thoroughbred racing horses. To discover the tendencies between CNV status and gene expression levels, we extracted CNVs of four Thoroughbred racing horses of which RNA sequencing was available. We found that 252 pairs of CNVs and genes were associated in the four horse samples. We did not observe a clear and consistent relationship between the deletion status of CNVs and gene expression levels before and after exercise in blood and muscle. However, we found some pairs of CNVs and associated genes that indicated relationships with gene expression levels: a positive relationship with genes responsible for membrane structure or cytoskeleton and a negative relationship with genes involved in disease. This study will lead to conceptual advances in understanding the relationship between CNVs and global gene expression in the horse.

Genomic instability, driver genes and cell selection: Projections from cancer to stem cells

Cancer cells and stem cells share many traits, including a tendency towards genomic instability. Human cancers exhibit tumor-specific genomic aberrations, which often affect their malignancy and drug response. During their culture propagation, human pluripotent stem cells (hPSCs) also acquire characteristic genomic aberrations, which may have significant impact on their molecular and cellular phenotypes. These aberrations vary in size from single nucleotide alterations to copy number alterations to whole chromosome gains. A prominent challenge in both cancer and stem cell research is to identify "driver aberrations" that confer a selection advantage, and "driver genes" that underlie the recurrence of these aberrations. Following principles that are already well-established in cancer research, candidate driver genes have also been suggested in hPSCs. Experimental validation of the functional role of such candidates can uncover whether these are bona fide driver genes. The identification of driver genes may bring us closer to a mechanistic understanding of the genomic instability of stem cells. Guided by terminologies and methodologies commonly applied in cancer research, such understanding may have important ramifications for both stem cell and cancer biology. This article is part of a Special Issue entitled: Stress as a fundamental theme in cell plasticity.

Genome wide CNV analysis reveals additional variants associated with milk production traits in Holsteins

Background

Milk production is an economically important sector of global agriculture. Much attention has been paid to the identification of quantitative trait loci (QTL) associated with milk, fat, and protein yield and the genetic and molecular mechanisms underlying them. Copy number variation (CNV) is an emerging class of variants which may be associated with complex traits.

Results

In this study, we performed a genome-wide association between CNVs and milk production traits in 26,362 Holstein bulls and cows. A total of 99 candidate CNVs were identified using Illumina BovineSNP50 array data, and association tests for each production trait were performed using a linear regression analysis with PCA correlation. A total of 34 CNVs on 22 chromosomes were significantly associated with at least one milk production trait after false discovery rate (FDR) correction. Some of those CNVs were located within or near known QTL for milk production traits. We further investigated the relationship between associated CNVs with neighboring SNPs. For all 82 combinations of traits and CNVs (less than 400 kb in length), we found 17 cases where CNVs directly overlapped with tag SNPs and 40 cases where CNVs were adjacent to tag SNPs. In 5 cases, CNVs located were in strong linkage disequilibrium with tag SNPs, either within or adjacent to the same haplotype block. There were an additional 20 cases where CNVs did not have a significant association with SNPs, suggesting that the effects of those CNVs were probably not captured by tag SNPs.

Conclusion

We conclude that combining CNV with SNP analyses reveals more genetic variations underlying milk production traits than those revealed by SNPs alone.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-683) contains supplementary material, which is available to authorized users.

Identification of structural variation in mouse genomes

Structural variation is variation in structure of DNA regions affecting DNA sequence length and/or orientation. It generally includes deletions, insertions, copy-number gains, inversions, and transposable elements. Traditionally, the identification of structural variation in genomes has been challenging. However, with the recent advances in high-throughput DNA sequencing and paired-end mapping (PEM) methods, the ability to identify structural variation and their respective association to human diseases has improved considerably. In this review, we describe our current knowledge of structural variation in the mouse, one of the prime model systems for studying human diseases and mammalian biology. We further present the evolutionary implications of structural variation on transposable elements. We conclude with future directions on the study of structural variation in mouse genomes that will increase our understanding of molecular architecture and functional consequences of structural variation.