Copy number variation modifies expression time courses

AAV1.NT-3 gene therapy in the SOD1KO mouse model of accelerated sarcopenia

BACKGROUND

Sarcopenia, an age-related loss of muscle mass, is a critical factor that affects the health of the older adults. The SOD1KO mouse is deficient of Cu/Zn superoxide dismutase, used as an accelerated aging model. We previously showed that NT-3 improves muscle fibre size by activating the mTOR pathway, suggesting a potential for attenuating age-related muscle loss. This study assessed the therapeutic efficacy of AAV1.NT-3 in this accelerated aging model.

METHODS

Twelve 6 months old SOD1KO mice were injected intramuscularly with a 1 × 1011 vg dose of AAV1.tMCK.NT-3, and 13 age-matched SOD1KO mice were used as controls. The treatment effect was evaluated using treadmill, rotarod and gait analyses as well as histological studies assessing changes in muscle fibre, and fibre type switch, in tibialis anterior, gastrocnemius, and triceps muscles, and myelin thickness by calculating G ratio in sciatic and tibial nerves. Molecular studies involved qPCR experiments to analyse the expression levels of mitochondrial and glycolysis markers and western blot experiments to assess the activity of mTORC1 pathway.

RESULTS

Treatment resulted in a 36% (154.9 vs. 114.1; P < 0.0001) and 76% increase (154.3 vs. 87.6; P < 0.0001) in meters ran, with treadmill test at 3 and 6 months post gene delivery. In addition, the treated cohort stayed on rotarod 30% (52.7 s vs. 40.4 s; P = 0.0095) and 54% (50.4 s vs. 32.7 s; P = 0.0007) longer, compared with untreated counterparts at 3 and 6 months post injection. Gait analysis, performed at endpoint, showed that stride width was normalized to wild type levels (29.3 mm) by an 11% decrease, compared with untreated cohort (28.6 mm vs. 32.1 mm; P = 0.0014). Compared with wild-type, SOD1KO mice showed 9.4% and 11.4% fibre size decrease in tibialis anterior and gastrocnemius muscles, respectively, which were normalized to wild type levels with treatment. Fibre diameter increase was observed prominently in FTG fibre type. G ratio analysis revealed hypomyelination in the tibial (0.721) and sciatic (0.676) nerves of SOD1KO model, which was reversed in the NT-3 cohort (0.646 and 0.634, respectively). Fibre size increase correlated with the increase in the p-S6 and p-4E-BP1 levels, and in the glycolysis markers in tibialis anterior. Alterations observed in the mitochondrial markers were not rescued with treatment. Overall, response to NT-3 was subdued in gastrocnemius muscle.

CONCLUSIONS

This study shows that AAV1.NT-3 gene therapy protected SOD1KO mouse from accelerated aging effects functionally and histologically. We further confirmed that NT-3 has potential to activate the mTOR and glycolytic pathways in muscle.

AAV1.NT-3 gene therapy prevents age-related sarcopenia

Sarcopenia is progressive loss of muscle mass and strength, occurring during normal aging with significant consequences on the quality of life for elderly. Neurotrophin 3 (NT-3) is an important autocrine factor supporting Schwann cell survival and differentiation and stimulating axon regeneration and myelination. NT-3 is involved in the maintenance of neuromuscular junction (NMJ) integrity, restoration of impaired radial growth of muscle fibers through activation of the Akt/mTOR pathway. We tested the efficacy of NT-3 gene transfer therapy in wild type (WT)-aged C57BL/6 mice, a model for natural aging and sarcopenia, via intramuscular injection 1 × 1011 vg AAV1.tMCK.NT-3, at 18 months of age. The treatment efficacy was assessed at 6 months post-injection using run to exhaustion and rotarod tests, in vivo muscle contractility assay, and histopathological studies of the peripheral nervous system, including NMJ connectivity and muscle. AAV1.NT-3 gene therapy in WT-aged C57BL/6 mice resulted in functional and in vivo muscle physiology improvements, supported by quantitative histology from muscle, peripheral nerves and NMJ. Hindlimb and forelimb muscles in the untreated cohort showed the presence of a muscle- and sex-dependent remodeling and fiber size decrease with aging, which was normalized toward values obtained from 10 months old WT mice with treatment. The molecular studies assessing the NT-3 effect on the oxidative state of distal hindlimb muscles, accompanied by western blot analyses for mTORC1 activation were in accordance with the histological findings. Considering the cost and quality of life to the individual, we believe our study has important implications for management of age-related sarcopenia.

An effect of large-scale deletions and duplications on transcript expression

Since copy number variants (CNVs) have been recognized as an important source of genetic and transcriptomic variation, we aimed to characterize the impact of CNVs located within coding, intergenic, upstream, and downstream gene regions on the expression of transcripts. Regions in which deletions occurred most often were introns, while duplications in coding regions. The transcript expression was lower for deleted coding (P = 0.008) and intronic regions (P = 1.355 × 10^-10), but it was not changed in the case of upstream and downstream gene regions (P = 0.085). Moreover, the expression was decreased if duplication occurred in the coding region (P = 8.318 × 10^-5). Furthermore, a negative correlation (r =  - 0.27) between transcript length and its expression was observed. The correlation between the percent of deleted/duplicated transcript and transcript expression level was not significant for all concerned genomic regions in five out of six animals. The exceptions were deletions in coding regions (P = 0.004) and duplications in introns (P = 0.01) in one individual. CNVs in coding (deletions, duplications) and intronic (deletions) regions are important modulators of transcripts by reducing their expression level. We hypothesize that deletions imply severe consequences by interrupting genes. The negative correlation between the size of the transcript and its expression level found in this study is consistent with the hypothesis that selection favours shorter introns and a moderate number of exons in highly expressed genes. This may explain the transcript expression reduction by duplications. We did not find the correlation between the size of deletions/duplications and transcript expression level suggesting that expression is modulated by CNVs regardless of their size.

AAV1.NT-3 gene therapy in a CMT2D model: phenotypic improvements in GarsP278KY/+ mice

Glycyl-tRNA synthetase mutations are associated to the Charcot-Marie-Tooth disease type-2D. The Gars^P278KY/+ model for Charcot-Marie-Tooth disease type-2D is known best for its early onset severe neuropathic phenotype with findings including reduced axon size, slow conduction velocities and abnormal neuromuscular junction. Muscle involvement remains largely unexamined. We tested the efficacy of neurotrophin 3 gene transfer therapy in two Gars mutants with severe (Gars^P278KY/+ ) and milder (Gars^ΔETAQ/+ ) phenotypes via intramuscular injection of adeno-associated virus setoype-1, triple tandem muscle creatine kinase promoter, neurotrophin 3 (AAV1.tMCK.NT-3) at 1 × 10¹¹ vg dose. In the Gars^P278KY/+ mice, the treatment efficacy was assessed at 12 weeks post-injection using rotarod test, electrophysiology and detailed quantitative histopathological studies of the peripheral nervous system including neuromuscular junction and muscle. Neurotrophin 3 gene transfer therapy in Gars^P278KY/+ mice resulted in significant functional and electrophysiological improvements, supported with increases in myelin thickness and improvements in the denervated status of neuromuscular junctions as well as increases in muscle fibre size along with attenuation of myopathic changes. Improvements in the milder phenotype Gars^ΔETAQ/+ was less pronounced. Furthermore, oxidative enzyme histochemistry in muscles from Gars mutants revealed alterations in the content and distribution of oxidative enzymes with increased expression levels of Pgc1a. Cox1, Cox3 and Atp5d transcripts were significantly decreased suggesting that the muscle phenotype might be related to mitochondrial dysfunction. Neurotrophin 3 gene therapy attenuated these abnormalities in the muscle. This study shows that neurotrophin 3 gene transfer therapy has disease modifying effect in a mouse model for Charcot-Marie-Tooth disease type-2D, leading to meaningful improvements in peripheral nerve myelination and neuromuscular junction integrity as well as in a unique myopathic process, associated with mitochondria dysfunction, all in combination contributing to functional outcome. Based on the multiple biological effects of this versatile molecule, we predict neurotrophin 3 has the potential to be beneficial in other aminoacyl-tRNA synthetase-linked Charcot-Marie-Tooth disease subtypes.

Genome-wide study of copy number variation implicates multiple novel loci for schizophrenia risk in Han Chinese family trios

Schizophrenia (SCZ) is a severe neuropsychiatric disorder that affects 1% of the global population. Copy number variations (CNVs) have been shown to play a critical role in its pathophysiology; however, only case-control studies on SCZ susceptibility CNVs have been conducted in Han Chinese. Here, we performed an array comparative genomic hybridization-based genome-wide CNV analysis in 100 Chinese family trios with SCZ. Burden test suggested that the SCZ probands carried more duplications than their healthy parents and unrelated healthy controls. Besides, five CNV loci were firstly reported to be associated with SCZ here, including both unbalanced transmitted CNVs and enriched de novo CNVs. Moreover, two genes (CTDSPL and MGAM) in these CNVs showed significant SCZ relevance in the expression level. Our findings support the crucial role of CNVs in the etiology of SCZ and provide new insights into the underlying mechanism of SCZ pathogenesis.

A de novo 13q31.3 microduplication encompassing the miR-17 ~ 92 cluster results in features mirroring those associated with Feingold syndrome 2

Hemizygosity of the MIR17HG gene encoding the miR-17 ~ 92 cluster is associated with Feingold syndrome 2 characterized by intellectual disability, skeletal abnormalities, short stature, and microcephaly. Here, we report on a female with a de novo 13q31.3 microduplication encompassing MIR17HG but excluding GPC5. She presented developmental delay, skeletal and digital abnormalities, and features such as tall stature and macrocephaly mirroring those of Feingold syndrome 2 patients. The limited extent of the proband's rearrangement to the miR cluster and the corresponding normal expression level of the neighboring GPC5 in her cells, together with previously described data on affected individuals of two families carrying overlapping duplications of the miR-17 ~ 92 cluster that comprise part of GPC5, who likewise presented macrocephaly, developmental delay, as well as skeletal, digital and stature abnormalities, allow to define a new syndrome due to independent microduplication of the miR-17 ~ 92 cluster.

Integrative analysis of genomic amplification-dependent expression and loss-of-function screen identifies ASAP1 as a driver gene in triple-negative breast cancer progression

The genetically heterogeneous triple-negative breast cancer (TNBC) continues to be an intractable disease, due to lack of effective targeted therapies. Gene amplification is a major event in tumorigenesis. Genes with amplification-dependent expression are being explored as therapeutic targets for cancer treatment. In this study, we have applied Analytical Multi-scale Identification of Recurring Events analysis and transcript quantification in the TNBC genome across 222 TNBC tumors and identified 138 candidate genes with positive correlation in copy number gain (CNG) and gene expression. siRNA-based loss-of-function screen of the candidate genes has validated EGFR, MYC, ASAP1, IRF2BP2, and CCT5 genes as drivers promoting proliferation in different TNBC cells. MYC, ASAP1, IRF2BP2, and CCT5 display frequent CNG and concurrent expression over 2173 breast cancer tumors (cBioPortal dataset). More frequently are MYC and ASAP1 amplified in TNBC tumors (>30%, n = 320). In particular, high expression of ASAP1, the ADP-ribosylation factor GTPase-activating protein, is significantly related to poor metastatic relapse-free survival of TNBC patients (n = 257, bc-GenExMiner). Furthermore, we have revealed that silencing of ASAP1 modulates numerous cytokine and apoptosis signaling components, such as IL1B, TRAF1, AIFM2, and MAP3K11 that are clinically relevant to survival outcomes of TNBC patients. ASAP1 has been reported to promote invasion and metastasis in various cancer cells. Our findings that ASAP1 is an amplification-dependent TNBC driver gene promoting TNBC cell proliferation, functioning upstream apoptosis components, and correlating to clinical outcomes of TNBC patients, support ASAP1 as a potential actionable target for TNBC treatment.

GluD1 knockout mice with a pure C57BL/6N background show impaired fear memory, social interaction, and enhanced depressive-like behavior

The GluD1 gene is associated with susceptibility for schizophrenia, autism, depression, and bipolar disorder. However, the function of GluD1 and how it is involved in these conditions remain elusive. In this study, we generated a Grid1 gene-knockout (GluD1-KO) mouse line with a pure C57BL/6N genetic background and performed several behavioral analyses. Compared to a control group, GluD1-KO mice showed no significant anxiety-related behavioral differences, evaluated using behavior in an open field, elevated plus maze, a light-dark transition test, the resident-intruder test of aggression and sensorimotor gating evaluated by the prepulse inhibition test. However, GluD1-KO mice showed (1) higher locomotor activity in the open field, (2) decreased sociability and social novelty preference in the three-chambered social interaction test, (3) impaired memory in contextual, but not cued fear conditioning tests, and (4) enhanced depressive-like behavior in a forced swim test. Pharmacological studies revealed that enhanced depressive-like behavior in GluD1-KO mice was restored by the serotonin reuptake inhibitors imipramine and fluoxetine, but not the norepinephrine transporter inhibitor desipramine. In addition, biochemical analysis revealed no significant difference in protein expression levels, such as other glutamate receptors in the synaptosome and postsynaptic densities prepared from the frontal cortex and the hippocampus. These results suggest that GluD1 plays critical roles in fear memory, sociability, and depressive-like behavior.

C57BL/6 and 129 inbred mouse strains differ in Gbp2 and Gbp2b expression in response to inflammatory stimuli in vivo

Background: Infections cause the production of inflammatory cytokines such as Interferon gamma (IFNγ). IFNγ in turn prompts the upregulation of a range of host defence proteins including members of the family of guanylate binding proteins (Gbps). In humans and mice alike, GBPs restrict the intracellular replication of invasive microbes and promote inflammation. To study the physiological functions of Gbp family members, the most commonly chosen in vivo models are mice harbouring loss-of-function mutations in either individual Gbp genes or the entire Gbp gene cluster on mouse chromosome 3. Individual Gbp deletion strains differ in their design, as some strains exist on a pure C57BL/6 genetic background, while other strains contain a 129-derived genetic interval encompassing the Gbp gene cluster on an otherwise C57BL/6 genetic background. Methods: To determine whether the presence of 129 alleles of paralogous Gbps could influence the phenotypes of 129-congenic Gbp-deficient strains, we studied the expression of Gbps in both C57BL/6J and 129/Sv mice following in vivo stimulation with adjuvants and after infection with either Toxoplasma gondii or Shigella flexneri. Results: We show that C57BL/6J relative to 129/Sv mice display moderately elevated expression of Gbp2, but more prominently, are also defective for Gbp2b (formerly Gbp1) mRNA induction upon immune priming. Notably, Toxoplasma infections induce robust Gbp2b protein expression in both strains of mice, suggestive of a Toxoplasma-activated mechanism driving Gbp2b protein translation. We further find that the higher expression of Gbp2b mRNA in 129/Sv mice correlates with a gene duplication event at the Gbp2b locus resulting in two copies of the Gbp2b gene on the haploid genome of the 129/Sv strain. Conclusions: Our findings demonstrate functional differences between 129 and C57BL/6 Gbp alleles which need to be considered in the design and interpretation of studies utilizing mouse models, particularly for phenotypes influenced by Gbp2 or Gbp2b expression.

An investigation of the role of gene copy number variations in sorafenib sensitivity in metastatic hepatocellular carcinoma patients

Background: Metastatic hepatocellular carcinoma (HCC) is a highly aggressive tumor with limited treatment options. While sorafenib has recently been shown to provide a survival advantage in patients with advanced HCC, the overall outcomes such as time to progression (TTP) and overall survival (OS) ought to be further improved. To that end, several targeted agents aimed at amplified oncogenes such as HER2 and FGFR2 have recently been developed. In this study, we aimed to identify genetic markers in the form of copy number variations (CNVs) that influence clinical outcomes post-sorafenib treatment in advanced HCC patients.

Methods: We surveyed 38 metastatic HCC patients who were treated with sorafenib for the presence of CNVs using the NanoString nCounter assay.

Results: The median TTP and OS for all patients were 2.7 months (95% confidence interval [CI]: 2.0-3.3 months) and 13.4 months (95% CI: 8.4-18.4 months), respectively. Several genes previously implicated in liver cancer were amplified, including CCND1 (n = 4), CDKN1A (n = 2), KRAS (n = 2), MDM2 (n = 1), and ERBB2 (n = 1). However, we found no correlations between CNVs and survival in our sorafenib-treated patients.

Conclusions: The clinical features and biomarkers that account for sensitivity to sorafenib in HCC are complicated and remain unclear. Further investigation to identify predictive biomarkers and therapeutic strategies, including combining sorafenib with other target agents, are warranted.

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.

Identifying anti-cancer drug response related genes using an integrative analysis of transcriptomic and genomic variations with cell line-based drug perturbations

Background

Clinical responses to anti-cancer therapies often only benefit a defined subset of patients. Predicting the best treatment strategy hinges on our ability to effectively translate genomic data into actionable information on drug responses.

Results

To achieve this goal, we compiled a comprehensive collection of baseline cancer genome data and drug response information derived from a large panel of cancer cell lines. This data set was applied to identify the signature genes relevant to drug sensitivity and their resistance by integrating CNVs and the gene expression of cell lines with in vitro drug responses. We presented an efficient in-silico pipeline for integrating heterogeneous cell line data sources with the simultaneous modeling of drug response values across all the drugs and cell lines. Potential signature genes correlated with drug response (sensitive or resistant) in different cancer types were identified. Using signature genes, our collaborative filtering-based drug response prediction model outperformed the 44 algorithms submitted to the DREAM competition on breast cancer cells. The functions of the identified drug response related signature genes were carefully analyzed at the pathway level and the synthetic lethality level. Furthermore, we validated these signature genes by applying them to the classification of the different subtypes of the TCGA tumor samples, and further uncovered their in vivo implications using clinical patient data.

Conclusions

Our work may have promise in translating genomic data into customized marker genes relevant to the response of specific drugs for a specific cancer type of individual patients.

Duplication of chicken defensin7 gene generated by gene conversion and homologous recombination

Defensins constitute an evolutionary conserved family of cationic antimicrobial peptides that play a key role in host innate immune responses to infection. Defensin genes generally reside in complex genomic regions that are prone to structural variation, and defensin genes exhibit extensive copy number variation in humans and in other species. Copy number variation of defensin genes was examined in inbred lines of Leghorn and Fayoumi chickens, and a duplication of defensin7 was discovered in the Fayoumi breed. Analysis of junction sequences confirmed the occurrence of a simple tandem duplication of defensin7 with sequence identity at the junction, suggesting nonallelic homologous recombination between defensin7 and defensin6 The duplication event generated two chimeric promoters that are best explained by gene conversion followed by homologous recombination. Expression of defensin7 was not elevated in animals with two genes despite both genes being transcribed in the tissues examined. Computational prediction of promoter regions revealed the presence of several putative transcription factor binding sites generated by the duplication event. These data provide insight into the evolution and possible function of large gene families and specifically, the defensins.

Mechanisms underlying structural variant formation in genomic disorders

With the recent burst of technological developments in genomics, and the clinical implementation of genome-wide assays, our understanding of the molecular basis of genomic disorders, specifically the contribution of structural variation to disease burden, is evolving quickly. Ongoing studies have revealed a ubiquitous role for genome architecture in the formation of structural variants at a given locus, both in DNA recombination-based processes and in replication-based processes. These reports showcase the influence of repeat sequences on genomic stability and structural variant complexity and also highlight the tremendous plasticity and dynamic nature of our genome in evolution, health and disease susceptibility.

Clinical and Genetic Heterogeneity of the 15q13.3 Microdeletion Syndrome

The 15q13.3 microdeletion is a recurrent CNV, presumably mediated by NAHR between segmental duplications in chromosome 15. The 15q13.3 deletion and duplication are associated with a wide range of clinical manifestations, such as intellectual deficits, seizures, autism, language and developmental delay, neuropsychiatric impairments, and behavioral problems illustrating incomplete penetrance and expressivity. This study comprises an evaluation of 106 symptomatic patients carrying the heterozygous deletion, as well as of 21 patients carrying the duplication, who have been described in previous studies. The analysis shows considerable heterogeneity for the manifestation of different key symptoms and familiar occurrence. Furthermore, 8 new patients are introduced. Convoluted familiar connections give new insights into the complexity of symptomatic manifestation. In previous studies, different opinions have been expressed as to the nature and precise location of the deletion breakpoints. Here, we show that not CHRNA7 and CHRFAM7A, but rather FAM7A or GOLGA8, serve as breakpoint regions concerning our patients. The deletion is described as heterogeneous in size. However, we assume that not only different breakpoints but also the imprecision of aCGH analysis on chromosome 15 due to segmental duplications accounts for the variability in size.

Genomic Profiling Reveals Unique Molecular Alterations in Hepatoblastomas and Adjacent Hepatocellular Carcinomas in B6C3F1 Mice

The cell of origin of hepatoblastoma (HB) in humans and mice is unknown; it is hypothesized to be a transformed hepatocyte, oval cell, or hepatic progenitor cell. In mice, current dogma is that HBs arise from preexisting hepatocellular neoplasms as a result of further neoplastic transformation. However, there is little evidence supporting this direct relationship. To better understand the relationship between hepatocellular carcinoma (HCC) and HB and determine molecular similarities between mouse and human HB, global gene expression analysis and targeted mutation analysis were performed using HB, HCC, and adjacent liver from the same animals in a recent National Toxicology Program bioassay. There were significant differences in Hras and Ctnnb1 mutation spectra, and by microarray, HBs showed dysregulation of embryonic development, stem cell pluripotency, and genomic imprinting compared to HCC. Meta-analysis showed similarities between HB, early mouse embryonic liver, and hepatocyte-derived stem/progenitor cells compared to HCC. Our data show that there are striking differences between HB and HCC and suggest that HB is a significantly different entity that may arise from a hepatic precursor cell. Furthermore, mouse HB is similar to the human disease at the pathway level and therefore is likely a relevant model for evaluating human cancer hazard.

Deciphering the Correlation between Breast Tumor Samples and Cell Lines by Integrating Copy Number Changes and Gene Expression Profiles

Breast cancer is one of the most common cancers with high incident rate and high mortality rate worldwide. Although different breast cancer cell lines were widely used in laboratory investigations, accumulated evidences have indicated that genomic differences exist between cancer cell lines and tissue samples in the past decades. The abundant molecular profiles of cancer cell lines and tumor samples deposited in the Cancer Cell Line Encyclopedia and The Cancer Genome Atlas now allow a systematical comparison of the breast cancer cell lines with breast tumors. We depicted the genomic characteristics of breast primary tumors based on the copy number variation and gene expression profiles and the breast cancer cell lines were compared to different subgroups of breast tumors. We identified that some of the breast cancer cell lines show high correlation with the tumor group that agrees with previous knowledge, while a big part of them do not, including the most used MCF7, MDA-MB-231, and T-47D. We presented a computational framework to identify cell lines that mostly resemble a certain tumor group for the breast tumor study. Our investigation presents a useful guide to bridge the gap between cell lines and tumors and helps to select the most suitable cell line models for personalized cancer studies.

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.

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.

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.

An Integrative Breakage Model of genome architecture, reshuffling and evolution: The Integrative Breakage Model of genome evolution, a novel multidisciplinary hypothesis for the study of genome plasticity

Our understanding of genomic reorganization, the mechanics of genomic transmission to offspring during germ line formation, and how these structural changes contribute to the speciation process, and genetic disease is far from complete. Earlier attempts to understand the mechanism(s) and constraints that govern genome remodeling suffered from being too narrowly focused, and failed to provide a unified and encompassing view of how genomes are organized and regulated inside cells. Here, we propose a new multidisciplinary Integrative Breakage Model for the study of genome evolution. The analysis of the high-level structural organization of genomes (nucleome), together with the functional constrains that accompany genome reshuffling, provide insights into the origin and plasticity of genome organization that may assist with the detection and isolation of therapeutic targets for the treatment of complex human disorders.

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.

Adaptive genomic structural variation in the grape powdery mildew pathogen, Erysiphe necator

BACKGROUND

Powdery mildew, caused by the obligate biotrophic fungus Erysiphe necator, is an economically important disease of grapevines worldwide. Large quantities of fungicides are used for its control, accelerating the incidence of fungicide-resistance. Copy number variations (CNVs) are unbalanced changes in the structure of the genome that have been associated with complex traits. In addition to providing the first description of the large and highly repetitive genome of E. necator, this study describes the impact of genomic structural variation on fungicide resistance in Erysiphe necator.

RESULTS

A shotgun approach was applied to sequence and assemble the genome of five E. necator isolates, and RNA-seq and comparative genomics were used to predict and annotate protein-coding genes. Our results show that the E. necator genome is exceptionally large and repetitive and suggest that transposable elements are responsible for genome expansion. Frequent structural variations were found between isolates and included copy number variation in EnCYP51, the target of the commonly used sterol demethylase inhibitor (DMI) fungicides. A panel of 89 additional E. necator isolates collected from diverse vineyard sites was screened for copy number variation in the EnCYP51 gene and for presence/absence of a point mutation (Y136F) known to result in higher fungicide tolerance. We show that an increase in EnCYP51 copy number is significantly more likely to be detected in isolates collected from fungicide-treated vineyards. Increased EnCYP51 copy numbers were detected with the Y136F allele, suggesting that an increase in copy number becomes advantageous only after the fungicide-tolerant allele is acquired. We also show that EnCYP51 copy number influences expression in a gene-dose dependent manner and correlates with fungal growth in the presence of a DMI fungicide.

CONCLUSIONS

Taken together our results show that CNV can be adaptive in the development of resistance to fungicides by providing increasing quantitative protection in a gene-dosage dependent manner. The results of this work not only demonstrate the effectiveness of using genomics to dissect complex traits in organisms with very limited molecular information, but also may have broader implications for understanding genomic dynamics in response to strong selective pressure in other pathogens with similar genome architectures.

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.

The effect of homozygous deletion of the BBOX1 and Fibin genes on carnitine level and acyl carnitine profile

Background

Carnitine is a key molecule in energy metabolism that helps transport activated fatty acids into the mitochondria. Its homeostasis is achieved through oral intake, renal reabsorption and de novo biosynthesis. Unlike dietary intake and renal reabsorption, the importance of de novo biosynthesis pathway in carnitine homeostasis remains unclear, due to lack of animal models and description of a single patient defective in this pathway.

Case presentation

We identified by array comparative genomic hybridization a 42 months-old girl homozygote for a 221 Kb interstitial deletions at 11p14.2, that overlaps the genes encoding Fibin and butyrobetaine-gamma 2-oxoglutarate dioxygenase 1 (BBOX1), an enzyme essential for the biosynthesis of carnitine de novo. She presented microcephaly, speech delay, growth retardation and minor facial anomalies. The levels of almost all evaluated metabolites were normal. Her serum level of free carnitine was at the lower limit of the reference range, while her acylcarnitine to free carnitine ratio was normal.

Conclusions

We present an individual with a completely defective carnitine de novo biosynthesis. This condition results in mildly decreased free carnitine level, but not in clinical manifestations characteristic of carnitine deficiency disorders, suggesting that dietary carnitine intake and renal reabsorption are sufficient to carnitine homeostasis. Our results also demonstrate that haploinsufficiency of BBOX1 and/or Fibin is not associated with Primrose syndrome as previously suggested.

The impact of chromosomal rearrangements on regulation of gene expression

The effects that coding region single-nucleotide polymorphisms or mutations have on gene expression have been well documented, predominantly owing to their association with disease. The effects of structural chromosomal rearrangements are also receiving increasing attention with the development of new techniques that allow accurate, high-resolution data, whether genomic interaction or transcriptome data, to be generated right down to the single-cell level. Over the past 18 months, these advances in experimental techniques have been used to further confirm and delineate the substantial effects that chromosome rearrangements can have on the regulation of gene expression and provide evidence of direct links between the two.

Integrated analysis of copy number variation and genome-wide expression profiling in colorectal cancer tissues

Integrative analyses of multiple genomic datasets for selected samples can provide better insight into the overall data and can enhance our knowledge of cancer. The objective of this study was to elucidate the association between copy number variation (CNV) and gene expression in colorectal cancer (CRC) samples and their corresponding non-cancerous tissues. Sixty-four paired CRC samples from the same patients were subjected to CNV profiling using the Illumina HumanOmni1-Quad assay, and validation was performed using multiplex ligation probe amplification method. Genome-wide expression profiling was performed on 15 paired samples from the same group of patients using the Affymetrix Human Gene 1.0 ST array. Significant genes obtained from both array results were then overlapped. To identify molecular pathways, the data were mapped to the KEGG database. Whole genome CNV analysis that compared primary tumor and non-cancerous epithelium revealed gains in 1638 genes and losses in 36 genes. Significant gains were mostly found in chromosome 20 at position 20q12 with a frequency of 45.31% in tumor samples. Examples of genes that were associated at this cytoband were PTPRT, EMILIN3 and CHD6. The highest number of losses was detected at chromosome 8, position 8p23.2 with 17.19% occurrence in all tumor samples. Among the genes found at this cytoband were CSMD1 and DLC1. Genome-wide expression profiling showed 709 genes to be up-regulated and 699 genes to be down-regulated in CRC compared to non-cancerous samples. Integration of these two datasets identified 56 overlapping genes, which were located in chromosomes 8, 20 and 22. MLPA confirmed that the CRC samples had the highest gains in chromosome 20 compared to the reference samples. Interpretation of the CNV data in the context of the transcriptome via integrative analyses may provide more in-depth knowledge of the genomic landscape of CRC.

Transcriptomic analysis of genetically defined autism candidate genes reveals common mechanisms of action

Background

Austism spectrum disorder (ASD) is a heterogeneous behavioral disorder or condition characterized by severe impairment of social engagement and the presence of repetitive activities. The molecular etiology of ASD is still largely unknown despite a strong genetic component. Part of the difficulty in turning genetics into disease mechanisms and potentially new therapeutics is the sheer number and diversity of the genes that have been associated with ASD and ASD symptoms. The goal of this work is to use shRNA-generated models of genetic defects proposed as causative for ASD to identify the common pathways that might explain how they produce a core clinical disability.

Methods

Transcript levels of Mecp2, Mef2a, Mef2d, Fmr1, Nlgn1, Nlgn3, Pten, and Shank3 were knocked-down in mouse primary neuron cultures using shRNA constructs. Whole genome expression analysis was conducted for each of the knockdown cultures as well as a mock-transduced culture and a culture exposed to a lentivirus expressing an anti-luciferase shRNA. Gene set enrichment and a causal reasoning engine was employed to identify pathway level perturbations generated by the transcript knockdown.

Results

Quantification of the shRNA targets confirmed the successful knockdown at the transcript and protein levels of at least 75% for each of the genes. After subtracting out potential artifacts caused by viral infection, gene set enrichment and causal reasoning engine analysis showed that a significant number of gene expression changes mapped to pathways associated with neurogenesis, long-term potentiation, and synaptic activity.

Conclusions

This work demonstrates that despite the complex genetic nature of ASD, there are common molecular mechanisms that connect many of the best established autism candidate genes. By identifying the key regulatory checkpoints in the interlinking transcriptional networks underlying autism, we are better able to discover the ideal points of intervention that provide the broadest efficacy across the diverse population of autism patients.

Structural variation-associated expression changes are paralleled by chromatin architecture modifications

Copy number variants (CNVs) influence the expression of genes that map not only within the rearrangement, but also to its flanks. To assess the possible mechanism(s) underlying this “neighboring effect”, we compared intrachromosomal interactions and histone modifications in cell lines of patients affected by genomic disorders and control individuals. Using chromosome conformation capture (4C-seq), we observed that a set of genes flanking the Williams-Beuren Syndrome critical region (WBSCR) were often looping together. The newly identified interacting genes include AUTS2, mutations of which are associated with autism and intellectual disabilities. Deletion of the WBSCR disrupts the expression of this group of flanking genes, as well as long-range interactions between them and the rearranged interval. We also pinpointed concomitant changes in histone modifications between samples.

We conclude that large genomic rearrangements can lead to chromatin conformation changes that extend far away from the structural variant, thereby possibly modulating expression globally and modifying the phenotype.

GEO Series accession number: GSE33784, GSE33867.

Identification of both copy number variation-type and constant-type core elements in a large segmental duplication region of the mouse genome

Background

Copy number variation (CNV), an important source of diversity in genomic structure, is frequently found in clusters called CNV regions (CNVRs). CNVRs are strongly associated with segmental duplications (SDs), but the composition of these complex repetitive structures remains unclear.

Results

We conducted self-comparative-plot analysis of all mouse chromosomes using the high-speed and large-scale-homology search algorithm SHEAP. For eight chromosomes, we identified various types of large SD as tartan-checked patterns within the self-comparative plots. A complex arrangement of diagonal split lines in the self-comparative-plots indicated the presence of large homologous repetitive sequences. We focused on one SD on chromosome 13 (SD13M), and developed SHEPHERD, a stepwise ab initio method, to extract longer repetitive elements and to characterize repetitive structures in this region. Analysis using SHEPHERD showed the existence of 60 core elements, which were expected to be the basic units that form SDs within the repetitive structure of SD13M. The demonstration that sequences homologous to the core elements (>70% homology) covered approximately 90% of the SD13M region indicated that our method can characterize the repetitive structure of SD13M effectively. Core elements were composed largely of fragmented repeats of a previously identified type, such as long interspersed nuclear elements (LINEs), together with partial genic regions. Comparative genome hybridization array analysis showed that whereas 42 core elements were components of CNVR that varied among mouse strains, 8 did not vary among strains (constant type), and the status of the others could not be determined. The CNV-type core elements contained significantly larger proportions of long terminal repeat (LTR) types of retrotransposon than the constant-type core elements, which had no CNV. The higher divergence rates observed in the CNV-type core elements than in the constant type indicate that the CNV-type core elements have a longer evolutionary history than constant-type core elements in SD13M.

Conclusions

Our methodology for the identification of repetitive core sequences simplifies characterization of the structures of large SDs and detailed analysis of CNV. The results of detailed structural and quantitative analyses in this study might help to elucidate the biological role of one of the SDs on chromosome 13.

Phenotypic impact of genomic structural variation: insights from and for human disease

Genomic structural variants have long been implicated in phenotypic diversity and human disease, but dissecting the mechanisms by which they exert their functional impact has proven elusive. Recently however, developments in high-throughput DNA sequencing and chromosomal engineering technology have facilitated the analysis of structural variants in human populations and model systems in unprecedented detail. In this Review, we describe how structural variants can affect molecular and cellular processes, leading to complex organismal phenotypes, including human disease. We further present advances in delineating disease-causing elements that are affected by structural variants, and we discuss future directions for research on the functional consequences of structural variants.

8p23.1 duplication syndrome; common, confirmed, and novel features in six further patients

The 8p23.1 duplication syndrome is a relatively rare genomic condition that has been confirmed with molecular cytogenetic methods in only 11 probands and five family members. Here, we describe another prenatal and five postnatal patients with de novo 8p23.1 duplications analyzed with oligonucleotide array comparative genomic hybridization (oaCGH). Of the common features, mild or moderate developmental delays and/or learning difficulties have been found in 11/12 postnatal probands, a variable degree of mild dysmorphism in 8/12 and congenital heart disease (CHD) in 4/5 prenatal and 3/12 postnatal probands. Behavioral problems, cleft lip and/or palate, macrocephaly, and seizures were confirmed as additional features among the new patients, and novel features included neonatal respiratory distress, attention deficit hyperactivity disorder (ADHD), ocular anomalies, balance problems, hypotonia, and hydrocele. The core duplication of 3.68 Mb contains 31 genes and microRNAs of which only GATA4, TNKS, SOX7, and XKR6 are likely to be dosage sensitive genes and MIR124-1 and MIR598 have been implicated in neurocognitive phenotypes. A combination of the duplication of GATA4, SOX7, and related genes may account for the variable penetrance of CHD. Two of the duplications were maternal and intrachromosomal in origin with maternal heterozygosity for the common inversion between the repeats in 8p23.1. These additional patients and the absence of the 8p23.1 duplications in published controls, indicate that the 8p23.1 duplication syndrome may now be considered a pathogenic copy number variation (pCNV) with an estimated population prevalence of 1 in 58,000.

Evaluating dynamic effects of copy number alterations on gene expression using a single transcription model

DNA copy number alterations (CNAs) are known to be related to genetic diseases, including cancer. The unlimited transcription (UT) model, in which transcription occurs permissively with a simple activation probability, has been proposed to investigate long-term effects of CNAs on gene expression values. Queueing theory was applied, and the copy-number-gene-expression relationship has been shown to be generally nonlinear in the UT model. However, the dynamic effects of CNAs on transcription and the underlying disorders related to diseases remain greatly unknown. Since most genes in a single cell are permissively transcribed in short periods of time interspersed by long periods of limited transcription, an alternative model for transcription in the restrictive state is needed for unraveling the effects of CNAs on gene expression levels with time. To address these issues, herein a single transcription (ST) model is proposed, in which bound TFs are assumed to be unloaded immediately after stimulating a transcription. Using the Laplace-Stieltjes transform and numerical analysis, the relationship between DNA copy number and gene expression level is evaluated. Dynamic modeling reveals that CNAs would potentially alter, or even reverse, the burst-like gene expression modifications while shifting from the ST model to the UT model. Moreover, functional disorders in transcriptional oscillation due to CNAs are shown via simulation. This paper demonstrates how mathematical theories could be helpful to interpret statistical findings from real data and achieve a better understanding of cancer biology.

In silico QTL mapping of maternal nurturing ability with the mouse diversity panel

Significant variation exists for maternal nurturing ability in inbred mice. Although classical mapping approaches have identified quantitative trait loci (QTL) that may account for this variation, the underlying genes are unknown. In this study, lactation performance data among the mouse diversity panel were used to map genomic regions associated with this variation. Females from each of 32 inbred strains (n = 8-19 dams/strain) were studied during the first 8 days of lactation by allowing them to raise weight- and size-normalized cross-foster litters (10 pups/litter). Average daily weight gain (ADG) of litters served as the primary indicator of milk production. The number of pups successfully reared to 8 days (PNUM8) also served as a related indicator of maternal performance. Initial haplotype association analysis using a Bonferroni-corrected, genome-wide threshold revealed 10 and 15 associations encompassing 11 and 13 genes for ADG and PNUM8, respectively. The most significant of these associated haplotype blocks were found on MMU 8, 11, and 19 and contained the genes Nr3c2, Egfr, Sec61g, and Gnaq. Lastly, two haplotype blocks on MMU9 were detected in association with PNUM8. These overlapped with the previously described maternal performance QTL, Neogq1. These results suggest that the application of in silico QTL mapping is a useful tool in discovering the presence of novel candidate genes involved in determining lactation capacity in mice.

RNA-Seq and human complex diseases: recent accomplishments and future perspectives

The availability of the human genome sequence has allowed identification of disease-causing mutations in many Mendelian disorders, and detection of significant associations of nucleotide polymorphisms to complex diseases and traits. Despite these progresses, finding the causative variations for most of the common diseases remains a complex task. Several studies have shown gene expression analyses provide a quite unbiased way to investigate complex traits and common disorders' pathogenesis. Therefore, whole-transcriptome analysis is increasingly acquiring a key role in the knowledge of mechanisms responsible for complex diseases. Hybridization- and tag-based technologies have elucidated the involvement of multiple genes and pathways in pathological conditions, providing insights into the expression of thousand of coding and noncoding RNAs, such as microRNAs. However, the introduction of Next-Generation Sequencing, particularly of RNA-Seq, has overcome some drawbacks of previously used technologies. Identifying, in a single experiment, potentially novel genes/exons and splice isoforms, RNA editing, fusion transcripts and allele-specific expression are some of its advantages. RNA-Seq has been fruitfully applied to study cancer and host-pathogens interactions, and it is taking first steps for studying neurodegenerative diseases (ND) as well as neuropsychiatric diseases. In addition, it is emerging as a very powerful tool to study quantitative trait loci associated with gene expression in complex diseases. This paper provides an overview on gene expression profiling of complex diseases, with emphasis on RNA-Seq, its advantages over conventional technologies for studying cancer and ND, and for linking nucleotide variations to gene expression changes, also discussing its limitations.

Gene expression analysis of matched ovarian primary tumors and peritoneal metastasis

Background

Ovarian cancer is the most deadly gynecological cancer due to late diagnosis at advanced stage with major peritoneal involvement. To date most research has focused on primary tumor. However the prognosis is directly related to residual disease at the end of the treatment. Therefore it is mandatory to focus and study the biology of meatastatic disease that is most frequently localized to the peritoneal caivty in ovarian cancer.

Methods

We used high-density gene expression arrays to investigate gene expression changes between matched primary and metastatic (peritoneal) lesions.

Results

Here we show that gene expression profiles in peritoneal metastasis are significantly different than their matched primary tumor and these changes are affected by underlying copy number variation differences among other causes. We show that differentially expressed genes are enriched in specific pathways including JAK/STAT pathway, cytokine signaling and other immune related pathways. We show that underlying copy number variations significantly affect gene expression. Indeed patients with important differences in copy number variation displayed greater gene expression differences between their primary and matched metastatic lesions.

Conclusions

Our analysis shows a very specific targeting at both the genomic and transcriptomic level to upregulate certain pathways in the peritoneal metastasis of ovarian cancer. Moreover, while primary tumors use certain pathways we identify distinct differences with metastatic lesions. The variation between primary and metastatic lesions should be considered in personalized treatment of ovarian cancer.

Mediation of Drosophila autosomal dosage effects and compensation by network interactions

BACKGROUND

Gene dosage change is a mild perturbation that is a valuable tool for pathway reconstruction in Drosophila. While it is often assumed that reducing gene dose by half leads to two-fold less expression, there is partial autosomal dosage compensation in Drosophila, which may be mediated by feedback or buffering in expression networks.

RESULTS

We profiled expression in engineered flies where gene dose was reduced from two to one. While expression of most one-dose genes was reduced, the gene-specific dose responses were heterogeneous. Expression of two-dose genes that are first-degree neighbors of one-dose genes in novel network models also changed, and the directionality of change depended on the response of one-dose genes.

CONCLUSIONS

Our data indicate that expression perturbation propagates in network space. Autosomal compensation, or the lack thereof, is a gene-specific response, largely mediated by interactions with the rest of the transcriptome.

Relating CNVs to transcriptome data at fine resolution: assessment of the effect of variant size, type, and overlap with functional regions

Copy-number variants (CNVs) form an abundant class of genetic variation with a presumed widespread impact on individual traits. While recent advances, such as the population-scale sequencing of human genomes, facilitated the fine-scale mapping of CNVs, the phenotypic impact of most of these CNVs remains unclear. By relating copy-number genotypes to transcriptome sequencing data, we have evaluated the impact of CNVs, mapped at fine scale, on gene expression. Based on data from 129 individuals with ancestry from two populations, we identified CNVs associated with the expression of 110 genes, with 13% of the associations involving complex, multiallelic CNVs. Categorization of CNVs according to variant type, size, and gene overlap enabled us to examine the impact of different CNV classes on expression variation. While many small (<4 kb) CNVs were associated with expression variation, overall we observed an enrichment of large duplications and deletions, including large intergenic CNVs, relative to the entire set of expression-associated CNVs. Furthermore, the copy number of genes intersecting with CNVs typically correlated positively with the genes' expression, and also was more strongly correlated with expression than nearby single nucleotide polymorphisms, suggesting a frequent causal role of CNVs in expression quantitative trait loci (eQTLs). We also elucidated unexpected cases of negative correlations between copy number and expression by assessing the CNVs' effects on the structure and regulation of genes. Finally, we examined dosage compensation of transcript levels. Our results suggest that association studies can gain in resolution and power by including fine-scale CNV information, such as those obtained from population-scale sequencing.

Copy number variants in candidate genes are genetic modifiers of Hirschsprung disease

Hirschsprung disease (HSCR) is a neurocristopathy characterized by absence of intramural ganglion cells along variable lengths of the gastrointestinal tract. The HSCR phenotype is highly variable with respect to gender, length of aganglionosis, familiality and the presence of additional anomalies. By molecular genetic analysis, a minimum of 11 neuro-developmental genes (RET, GDNF, NRTN, SOX10, EDNRB, EDN3, ECE1, ZFHX1B, PHOX2B, KIAA1279, TCF4) are known to harbor rare, high-penetrance mutations that confer a large risk to the bearer. In addition, two other genes (RET, NRG1) harbor common, low-penetrance polymorphisms that contribute only partially to risk and can act as genetic modifiers. To broaden this search, we examined whether a set of 67 proven and candidate HSCR genes harbored additional modifier alleles. In this pilot study, we utilized a custom-designed array CGH with ∼33,000 test probes at an average resolution of ∼185 bp to detect gene-sized or smaller copy number variants (CNVs) within these 67 genes in 18 heterogeneous HSCR patients. Using stringent criteria, we identified CNVs at three loci (MAPK10, ZFHX1B, SOX2) that are novel, involve regulatory and coding sequences of neuro-developmental genes, and show association with HSCR in combination with other congenital anomalies. Additional CNVs are observed under relaxed criteria. Our research suggests a role for CNVs in HSCR and, importantly, emphasizes the role of variation in regulatory sequences. A much larger study will be necessary both for replication and for identifying the full spectrum of small CNV effects.