Subgroup-specific structural variation across 1,000 medulloblastoma genomes.

Medulloblastoma, the most common malignant paediatric brain tumour, is currently treated with nonspecific cytotoxic therapies including surgery, whole-brain radiation, and aggressive chemotherapy. As medulloblastoma exhibits marked intertumoural heterogeneity, with at least four distinct molecular variants, previous attempts to identify targets for therapy have been underpowered because of small samples sizes. Here we report somatic copy number aberrations (SCNAs) in 1,087 unique medulloblastomas. SCNAs are common in medulloblastoma, and are predominantly subgroup-enriched. The most common region of focal copy number gain is a tandem duplication of SNCAIP, a gene associated with Parkinson's disease, which is exquisitely restricted to Group 4α. Recurrent translocations of PVT1, including PVT1-MYC and PVT1-NDRG1, that arise through chromothripsis are restricted to Group 3. Numerous targetable SCNAs, including recurrent events targeting TGF-β signalling in Group 3, and NF-κB signalling in Group 4, suggest future avenues for rational, targeted therapy.

Understanding mutational and selective processes that govern gene duplication in mammals

Abstract : Gene duplication is an essential source of material for the origin of genetic novelties. The reverse transcription of source gene mRNA followed by the genomic insertion of the resulting cDNA - retroposition - has provided the human genome with at least ~3600 detectable retrocopies. We find that ~30% of these retrocopies are transcribed, generally in testes. Their transcription often relies on preexisting regulatory elements (or open chromatin) close to their insertion site, which is illustrated by mRNA molecules containing retrocopies fused to their neighboring genes. Retrocopies appear to have been profoundly shaped by selection. Consistently, human retrocopies with an intact open reading (ORF) are more often transcribed than retropseudogenes, which leads to a minimal estimate of 120 functional retrogenes present in our genome. We also performed an analysis of Ka/Ks for human retrocopies. This analysis demonstrates that several intact retrocopies evolved under purifying selection and yields an estimated formation rate of ~1 retrogene per million year in the primate lineage. Using DNA sequencing and evolutionary simulations, we have identified 7 such primate-specific retrogenes that emerged on the lineage leading to humans In therian genomes, we found an excess of retrogenes with X-linked parents. Expression analyses support the idea that this "out of X" movement was driven by natural selection to produce autosomal functional counterparts for X-linked genes, which are silenced during male meiosis. Phylogenetic dating of this "out of X" movement suggests that our sex chromosomes arose about 180 MYA ago and are thus much younger than previously thought. Finally, we have also analyzed young gene duplications (and deletions) that arose by non allelic-homologous recombination and are not fixed in species. Using wild-caught and laboratory animals, we detected thousands of DNA segments that are polymorphic in copy number in mice. These copy number variants were found to profoundly alter the transcriptome of several mouse tissues. Strikingly, their influence on gene expression is not limited to the gene they contain but seems to extend to genes located up to 1.5 million bases away.

Origin and evolution of homologous repeated sequences in the mitochondrial DNA control region of shrews.

The complete mitochondrial DNA (mtDNA) control region was amplified and directly sequenced in two species of shrew, Crocidura russula and Sorex araneus (Insectivora, Mammalia). The general organization is similar to that found in other mammals: a central conserved region surrounded by two more variable domains. However, we have found in shrews the simultaneous presence of arrays of tandem repeats in potential locations where repeats tend to occur separately in other mammalian species. These locations correspond to regions which are associated with a possible interruption of the replication processes, either at the end of the three-stranded D-loop structure or toward the end of the heavy-strand replication. In the left domain the repeated sequences (R1 repeats) are 78 bp long, whereas in the right domain the repeats are 12 bp long in C. russula and 14 bp long in S. araneus (R2 repeats). Variation in the copy number of these repeated sequences results in mtDNA control region length differences. Southern blot analysis indicates that level of heteroplasmy (more than one mtDNA form within an individual) differs between species. A comparative study of the R2 repeats in 12 additional species representing three shrew subfamilies provides useful indications for the understanding of the origin and the evolution of these homologous tandemly repeated sequences. An asymmetry in the distribution of variants within the arrays, as well as the constant occurrence of shorter repeated sequences flanking only one side of the R2 arrays, could be related to asymmetry in the replication of each strand of the mtDNA molecule. The pattern of sequence and length variation within and between species, together with the capability of the arrays to form stable secondary structures, suggests that the dominant mechanism involved in the evolution of these arrays in unidirectional replication slippage.

Calcium phosphate transfection generates mammalian recombinant cell lines with higher specific productivity than polyfection.

Transfection with polyethylenimine (PEI) was evaluated as a method for the generation of recombinant Chinese hamster ovary (CHO DG44) cell lines by direct comparison with calcium phosphate-DNA coprecipitation (CaPO4) using both green fluorescent protein (GFP) and a monoclonal antibody as reporter proteins. Following transfection with a GFP expression vector, the proportion of GFP-positive cells as determined by flow cytometry was fourfold higher for the PEI transfection as compared to the CaPO4 transfection. However, the mean level of transient GFP expression for the cells with the highest level of fluorescence was twofold greater for the CaPO4 transfection. Fluorescence in situ hybridization on metaphase chromosomes from pools of cells grown under selective pressure demonstrated that plasmid integration always occurred at a single site regardless of the transfection method. Importantly, the copy number of integrated plasmids was measurably higher in cells transfected with CaPO4. The efficiency of recombinant cell line recovery under selective pressure was fivefold higher following PEI transfection, but the average specific productivity of a recombinant antibody was about twofold higher for the CaPO4-derived cell lines. Nevertheless, no difference between the two transfection methods was observed in terms of the stability of protein production. These results demonstrated the feasibility of generating recombinant CHO-derived cell lines by PEI transfection. However, this method appeared inferior to CaPO4 transfection with regard to the specific productivity of the recovered cell lines.

Islands of euchromatin-like sequence and expressed polymorphic sequences within the short arm of human chromosome 21.

The goals of the human genome project did not include sequencing of the heterochromatic regions. We describe here an initial sequence of 1.1 Mb of the short arm of human chromosome 21 (HSA21p), estimated to be 10% of 21p. This region contains extensive euchromatic-like sequence and includes on average one transcript every 100 kb. These transcripts show multiple inter- and intrachromosomal copies, and extensive copy number and sequence variability. The sequencing of the "heterochromatic" regions of the human genome is likely to reveal many additional functional elements and provide important evolutionary information.

Validation of real-time methylation-specific PCR to determine O6-methylguanine-DNA methyltransferase gene promoter methylation in glioma.

Epigenetic silencing of the DNA repair protein O(6)-methylguanine-DNA methyltransferase (MGMT) by promoter methylation predicts successful alkylating agent therapy, such as with temozolomide, in glioblastoma patients. Stratified therapy assignment of patients in prospective clinical trials according to tumor MGMT status requires a standardized diagnostic test, suitable for high-throughput analysis of small amounts of formalin-fixed, paraffin-embedded tumor tissue. A direct, real-time methylation-specific PCR (MSP) assay was developed to determine methylation status of the MGMT gene promoter. Assay specificity was obtained by selective amplification of methylated DNA sequences of sodium bisulfite-modified DNA. The copy number of the methylated MGMT promoter, normalized to the beta-actin gene, provides a quantitative test result. We analyzed 134 clinical glioma samples, comparing the new test with the previously validated nested gel-based MSP assay, which yields a binary readout. A cut-off value for the MGMT methylation status was suggested by fitting a bimodal normal mixture model to the real-time results, supporting the hypothesis that there are two distinct populations within the test samples. Comparison of the tests showed high concordance of the results (82/91 [90%]; Cohen's kappa = 0.80; 95% confidence interval, 0.82-0.95). The direct, real-time MSP assay was highly reproducible (Pearson correlation 0.996) and showed valid test results for 93% (125/134) of samples compared with 75% (94/125) for the nested, gel-based MSP assay. This high-throughput test provides an important pharmacogenomic tool for individualized management of alkylating agent chemotherapy.

The Waddlia genome: a window into chlamydial biology.

Growing evidence suggests that a novel member of the Chlamydiales order, Waddlia chondrophila, is a potential agent of miscarriage in humans and abortion in ruminants. Due to the lack of genetic tools to manipulate chlamydia, genomic analysis is proving to be the most incisive tool in stimulating investigations into the biology of these obligate intracellular bacteria. 454/Roche and Solexa/Illumina technologies were thus used to sequence and assemble de novo the full genome of the first representative of the Waddliaceae family, W. chondrophila. The bacteria possesses a 2'116'312 bp chromosome and a 15'593 bp low-copy number plasmid that might integrate into the bacterial chromosome. The Waddlia genome displays numerous repeated sequences indicating different genome dynamics from classical chlamydia which almost completely lack repetitive elements. Moreover, W. chondrophila exhibits many virulence factors also present in classical chlamydia, including a functional type III secretion system, but also a large complement of specific factors for resistance to host or environmental stresses. Large families of outer membrane proteins were identified indicating that these highly immunogenic proteins are not Chlamydiaceae specific and might have been present in their last common ancestor. Enhanced metabolic capability for the synthesis of nucleotides, amino acids, lipids and other co-factors suggests that the common ancestor of the modern Chlamydiales may have been less dependent on their eukaryotic host. The fine-detailed analysis of biosynthetic pathways brings us closer to possibly developing a synthetic medium to grow W. chondrophila, a critical step in the development of genetic tools. As a whole, the availability of the W. chondrophila genome opens new possibilities in Chlamydiales research, providing new insights into the evolution of members of the order Chlamydiales and the biology of the Waddliaceae.

Identifying the parent-of-origin of de novo SNVs in schizophrenia

Several studies over the last few years have shown that newly arising (de novo) mutations contribute to the genetics of schizophrenia (SZ), autism (ASD) and other developmental disorders. The strongest evidence comes from studies of de novo Copy Number Variation (CNV), where the rate of new mutations is shown to be increased in cases when compared to controls [23, 24]. Research on de novo point mutations and small insertion-deletions (indels) has been more limited, but with the development of next-generation sequencing (NGS) technology, such studies are beginning to provide preliminary evidence that de novo single-nucleotide mutations (SNVs) might also increase risk of SZ and ASD [25, 26] Advanced paternal age is a major source of new mutations in human beings [27] and could thus be associated with increased risk for developing SZ, ASD or other developmental disorders. Indeed, advanced paternal age is found to be a risk factor for developing SZ and ASD in the offspring [28, 29] and new mutations related to advanced paternal age have been implicated as a cause of sporadic cases in several autosomal dominant diseases, some neurodevelopmental diseases, including SZ and ASD, and social functioning. New single-base substitutions occur at higher rates at males compared to females and this difference increases with paternal age. This is due to the fact that sperm cells go through a much higher number of cell divisions (~840 by the age of 50), which increases the risk for DNA copy errors in the male germ line [30] . By contrast, the female eggs (oocytes) undergo only 24 cell divisions and all but the last occur during foetal life. The aim of my project is to determine the parent-of-origin of de novo SNVs, using large samples of parent-offspring trios affected with schizophrenia (SZ). From whole exome sequencing of 618 Bulgarian proband-offspring trios affected, nearly 1000 de novo (SNVs or small indels) have been identified and from these, the parent-of-origin of at least 60% of the mutations (N=600) can be established. This project is contained in a main one that consists on the determination of the parental origin of different types of de novo mutations (SNVs, small indels and large CNVs).

Structural variation and its effect on expression.

Structural variation, whether it is caused by copy number variants or present in a balanced form, such as reciprocal translocations and inversions, can have a profound and dramatic effect on the expression of genes mapping within and close to the rearrangement, as well as affecting others genome wide. These effects can be caused by altering the copy number of one or more genes or regulatory elements (dosage effect) or from physical disruption of links between regulatory elements and their associated gene or genes, resulting in perturbation of expression. Similarly, large-scale structural variants can result in genome-wide expression changes by altering the positions that chromosomes occupy within the nucleus, potentially disrupting not only local cis interactions, but also trans interactions that occur throughout the genome. Structural variation is, therefore, a significant factor in the study of gene expression and is discussed here in more detail.

Towards the identification of a genetic basis for Landau-Kleffner syndrome.

OBJECTIVE: To establish the genetic basis of Landau-Kleffner syndrome (LKS) in a cohort of two discordant monozygotic (MZ) twin pairs and 11 isolated cases. METHODS: We used a multifaceted approach to identify genetic risk factors for LKS. Array comparative genomic hybridization (CGH) was performed using the Agilent 180K array. Whole genome methylation profiling was undertaken in the two discordant twin pairs, three isolated LKS cases, and 12 control samples using the Illumina 27K array. Exome sequencing was undertaken in 13 patients with LKS including two sets of discordant MZ twins. Data were analyzed with respect to novel and rare variants, overlapping genes, variants in reported epilepsy genes, and pathway enrichment. RESULTS: A variant (cG1553A) was found in a single patient in the GRIN2A gene, causing an arginine to histidine change at site 518, a predicted glutamate binding site. Following copy number variation (CNV), methylation, and exome sequencing analysis, no single candidate gene was identified to cause LKS in the remaining cohort. However, a number of interesting additional candidate variants were identified including variants in RELN, BSN, EPHB2, and NID2. SIGNIFICANCE: A single mutation was identified in the GRIN2A gene. This study has identified a number of additional candidate genes including RELN, BSN, EPHB2, and NID2. A PowerPoint slide summarizing this article is available for download in the Supporting Information section here.

Gene expression patterns unveil a new level of molecular heterogeneity in colorectal cancer.

The recognition that colorectal cancer (CRC) is a heterogeneous disease in terms of clinical behaviour and response to therapy translates into an urgent need for robust molecular disease subclassifiers that can explain this heterogeneity beyond current parameters (MSI, KRAS, BRAF). Attempts to fill this gap are emerging. The Cancer Genome Atlas (TGCA) reported two main CRC groups, based on the incidence and spectrum of mutated genes, and another paper reported an EMT expression signature defined subgroup. We performed a prior free analysis of CRC heterogeneity on 1113 CRC gene expression profiles and confronted our findings to established molecular determinants and clinical, histopathological and survival data. Unsupervised clustering based on gene modules allowed us to distinguish at least five different gene expression CRC subtypes, which we call surface crypt-like, lower crypt-like, CIMP-H-like, mesenchymal and mixed. A gene set enrichment analysis combined with literature search of gene module members identified distinct biological motifs in different subtypes. The subtypes, which were not derived based on outcome, nonetheless showed differences in prognosis. Known gene copy number variations and mutations in key cancer-associated genes differed between subtypes, but the subtypes provided molecular information beyond that contained in these variables. Morphological features significantly differed between subtypes. The objective existence of the subtypes and their clinical and molecular characteristics were validated in an independent set of 720 CRC expression profiles. Our subtypes provide a novel perspective on the heterogeneity of CRC. The proposed subtypes should be further explored retrospectively on existing clinical trial datasets and, when sufficiently robust, be prospectively assessed for clinical relevance in terms of prognosis and treatment response predictive capacity. Original microarray data were uploaded to the ArrayExpress database (http://www.ebi.ac.uk/arrayexpress/) under Accession Nos E-MTAB-990 and E-MTAB-1026. © 2013 Swiss Institute of Bioinformatics. Journal of Pathology published by John Wiley & Sons Ltd on behalf of Pathological Society of Great Britain and Ireland.

Combining chromosomal arm status and significantly aberrant genomic locations reveals new cancer subtypes.

Many types of tumors exhibit characteristic chromosomal losses or gains, as well as local amplifications and deletions. Within any given tumor type, sample specific amplifications and deletions are also observed. Typically, a region that is aberrant in more tumors, or whose copy number change is stronger, would be considered as a more promising candidate to be biologically relevant to cancer. We sought for an intuitive method to define such aberrations and prioritize them. We define V, the "volume" associated with an aberration, as the product of three factors: (a) fraction of patients with the aberration, (b) the aberration's length and (c) its amplitude. Our algorithm compares the values of V derived from the real data to a null distribution obtained by permutations, and yields the statistical significance (p-value) of the measured value of V. We detected genetic locations that were significantly aberrant, and combine them with chromosomal arm status (gain/loss) to create a succinct fingerprint of the tumor genome. This genomic fingerprint is used to visualize the tumors, highlighting events that are co-occurring or mutually exclusive. We apply the method on three different public array CGH datasets of Medulloblastoma and Neuroblastoma, and demonstrate its ability to detect chromosomal regions that were known to be altered in the tested cancer types, as well as to suggest new genomic locations to be tested. We identified a potential new subtype of Medulloblastoma, which is analogous to Neuroblastoma type 1.

Asymmetric stochastic switching driven by intrinsic molecular noise

Low-copy-number molecules are involved in many functions in cells. The intrinsic fluctuations of these numbers can enable stochastic switching between multiple steady states, inducing phenotypic variability. Herein we present a theoretical and computational study based on Master Equations and Fokker-Planck and Langevin descriptions of stochastic switching for a genetic circuit of autoactivation. We show that in this circuit the intrinsic fluctuations arising from low-copy numbers, which are inherently state-dependent, drive asymmetric switching. These theoretical results are consistent with experimental data that have been reported for the bistable system of the gallactose signaling network in yeast. Our study unravels that intrinsic fluctuations, while not required to describe bistability, are fundamental to understand stochastic switching and the dynamical relative stability of multiple states.

Association analyses of 249,796 individuals reveal 18 new loci associated with body mass index.

Obesity is globally prevalent and highly heritable, but its underlying genetic factors remain largely elusive. To identify genetic loci for obesity susceptibility, we examined associations between body mass index and ∼ 2.8 million SNPs in up to 123,865 individuals with targeted follow up of 42 SNPs in up to 125,931 additional individuals. We confirmed 14 known obesity susceptibility loci and identified 18 new loci associated with body mass index (P < 5 × 10⁻⁸), one of which includes a copy number variant near GPRC5B. Some loci (at MC4R, POMC, SH2B1 and BDNF) map near key hypothalamic regulators of energy balance, and one of these loci is near GIPR, an incretin receptor. Furthermore, genes in other newly associated loci may provide new insights into human body weight regulation.