Identification of human mitochondrial RNA cleavage sites and candidate RNA processing factors.

BACKGROUND: The human mitochondrial genome is transcribed as long strands of RNA containing multiple genes, which require post-transcriptional cleavage and processing to release functional gene products that play vital roles in cellular energy production. Despite knowledge implicating mitochondrial post-transcriptional processes in pathologies such as cancer, cardiovascular disease and diabetes, very little is known about the way their function varies on a human population level and what drives changes in these processes to ultimately influence disease risk. Here, we develop a method to detect and quantify mitochondrial RNA cleavage events from standard RNA sequencing data and apply this approach to human whole blood data from > 1000 samples across independent cohorts. RESULTS: We detect 54 putative mitochondrial RNA cleavage sites that not only map to known gene boundaries, short RNA ends and RNA modification sites, but also occur at internal gene positions, suggesting novel mitochondrial RNA cleavage junctions. Inferred RNA cleavage rates correlate with mitochondrial-encoded gene expression across individuals, suggesting an impact on downstream processes. Furthermore, by comparing inferred cleavage rates to nuclear genetic variation and gene expression, we implicate multiple genes in modulating mitochondrial RNA cleavage (e.g. MRPP3, TBRG4 and FASTKD5), including a potentially novel role for RPS19 in influencing cleavage rates at a site near to the MTATP6-COX3 junction that we validate using shRNA knock down data. CONCLUSIONS: We identify novel cleavage junctions associated with mitochondrial RNA processing, as well as genes newly implicated in these processes, and detect the potential impact of variation in cleavage rates on downstream phenotypes and disease processes. These results highlight the complexity of the mitochondrial transcriptome and point to novel mechanisms through which nuclear-encoded genes can potentially influence key mitochondrial processes.

Development of a Novel NGS Methodology for Ultrasensitive Circulating Tumor DNA Detection as a Tool for Early-Stage Breast Cancer Diagnosis.

Breast cancer (BC) is the most prevalent cancer in women. While usually detected when localized, invasive procedures are still required for diagnosis. Herein, we developed a novel ultrasensitive pipeline to detect circulating tumor DNA (ctDNA) in a series of 75 plasma samples from localized BC patients prior to any medical intervention. We first performed a tumor-informed analysis to correlate the mutations found in tumor tissue and plasma. Disregarding the tumor data next, we developed an approach to detect tumor mutations in plasma. We observed a mutation concordance between the tumor and plasma of 29.50% with a sensitivity down to 0.03% in mutant variant allele frequency (VAF). We detected mutations in 33.78% of the samples, identifying eight patients with plasma-only mutations. Altogether, we determined a specificity of 86.36% and a positive predictive value of 88.46% for BC detection. We demonstrated an association between higher ctDNA median VAF and higher tumor grade, multiple plasma mutations with a likelihood of relapse and more frequent TP53 plasma mutations in hormone receptor-negative tumors. Overall, we have developed a unique ultra-sensitive sequencing workflow with a technology not previously employed in early BC, paving the way for its application in BC screening.

Protection against de novo methylation is instrumental in maintaining parent-of-origin methylation inherited from the gametes.

Identifying loci with parental differences in DNA methylation is key to unraveling parent-of-origin phenotypes. By conducting a MeDIP-Seq screen in maternal-methylation free postimplantation mouse embryos (Dnmt3L-/+), we demonstrate that maternal-specific methylation exists very scarcely at midgestation. We reveal two forms of oocyte-specific methylation inheritance: limited to preimplantation, or with longer duration, i.e. maternally imprinted loci. Transient and imprinted maternal germline DMRs (gDMRs) are indistinguishable in gametes and preimplantation embryos, however, de novo methylation of paternal alleles at implantation delineates their fates and acts as a major leveling factor of parent-inherited differences. We characterize two new imprinted gDMRs, at the Cdh15 and AK008011 loci, with tissue-specific imprinting loss, again by paternal methylation gain. Protection against demethylation after fertilization has been emphasized as instrumental in maintaining parent-of-origin methylation inherited from the gametes. Here we provide evidence that protection against de novo methylation acts as an equal major pivot, at implantation and throughout life.

Comparative methylomics reveals gene-body H3K36me3 in Drosophila predicts DNA methylation and CpG landscapes in other invertebrates.

In invertebrates that harbor functional DNA methylation enzymatic machinery, gene-bodies are the primary targets for CpG methylation. However, virtually all other aspects of invertebrate DNA methylation have remained a mystery until now. Here, using a comparative methylomics approach, we demonstrate that Nematostella vectensis, Ciona intestinalis, Apis mellifera, and Bombyx mori show two distinct populations of genes differentiated by gene-body CpG density. Genome-scale DNA methylation profiles for A. mellifera spermatozoa reveal CpG-poor genes are methylated in the germline, as predicted by the depletion of CpGs. We find an evolutionarily conserved distinction between CpG-poor and GpC-rich genes: The former are associated with basic biological processes, the latter with more specialized functions. This distinction is strikingly similar to that recently observed between euchromatin-associated genes in Drosophila that contain intragenic histone 3 lysine 36 trimethylation (H3K36me3) and those that do not, even though Drosophila does not display CpG density bimodality or methylation. We confirm that a significant number of CpG-poor genes in N. vectensis, C. intestinalis, A. mellifera, and B. mori are orthologs of H3K36me3-rich genes in Drosophila. We propose that over evolutionary time, gene-body H3K36me3 has influenced gene-body DNA methylation levels and, consequently, the gene-body CpG density bimodality characteristic of invertebrates that harbor CpG methylation.

Increased blood draws for ultrasensitive ctDNA and CTCs detection in early breast cancer patients.

Early breast cancer patients often experience relapse due to residual disease after treatment. Liquid biopsy is a methodology capable of detecting tumor components in blood, but low concentrations at early stages pose challenges. To detect them, next-generation sequencing has promise but entails complex processes. Exploring larger blood volumes could overcome detection limitations. Herein, a total of 282 high-volume plasma and blood-cell samples were collected for dual ctDNA/CTCs detection using a single droplet-digital PCR assay per patient. ctDNA and/or CTCs were detected in 100% of pre-treatment samples. On the other hand, post-treatment positive samples exhibited a minimum variant allele frequency of 0.003% for ctDNA and minimum cell number of 0.069 CTCs/mL of blood, surpassing previous investigations. Accurate prediction of residual disease before surgery was achieved in patients without a complete pathological response. A model utilizing ctDNA dynamics achieved an area under the ROC curve of 0.92 for predicting response. We detected disease recurrence in blood in the three patients who experienced a relapse, anticipating clinical relapse by 34.61, 9.10, and 7.59 months. This methodology provides an easily implemented alternative for ultrasensitive residual disease detection in early breast cancer patients.

Bionemo: molecular information on biodegradation metabolism.

Bionemo (http://bionemo.bioinfo.cnio.es) stores manually curated information about proteins and genes directly implicated in the Biodegradation metabolism. When possible, the database includes information on sequence, domains and structures for proteins; and sequence, regulatory elements and transcription units for genes. Thus, Bionemo is a unique resource that complements other biodegradation databases such as the University of Minessota Biocatalysis/Biodegradation Database, or Metarouter, which focus more on the biochemical aspects of biodegradation than in the nature of the biomolecules carrying out the reactions. Bionemo has been built by manually associating sequences database entries to biodegradation reactions, using the information extracted from published articles. Information on transcription units and their regulation was also extracted from the literature for biodegradation genes, and linked to the underlying biochemical network. In its current version, Bionemo contains sequence information for 324 reactions and transcription regulation information for more than 100 promoters and 100 transcription factors. The information in the Bionemo database is available via a web server and the full database is also downloadable as a PostgresSQL dump. To facilitate the programmatic use of the information contained in the database, an object-oriented Perl API is also provided.

Nuclear genetic regulation of the human mitochondrial transcriptome.

Mitochondria play important roles in cellular processes and disease, yet little is known about how the transcriptional regime of the mitochondrial genome varies across individuals and tissues. By analyzing >11,000 RNA-sequencing libraries across 36 tissue/cell types, we find considerable variation in mitochondrial-encoded gene expression along the mitochondrial transcriptome, across tissues and between individuals, highlighting the importance of cell-type specific and post-transcriptional processes in shaping mitochondrial-encoded RNA levels. Using whole-genome genetic data we identify 64 nuclear loci associated with expression levels of 14 genes encoded in the mitochondrial genome, including missense variants within genes involved in mitochondrial function (TBRG4, MTPAP and LONP1), implicating genetic mechanisms that act in trans across the two genomes. We replicate ~21% of associations with independent tissue-matched datasets and find genetic variants linked to these nuclear loci that are associated with cardio-metabolic phenotypes and Vitiligo, supporting a potential role for variable mitochondrial-encoded gene expression in complex disease.

Loss of Trem2 in microglia leads to widespread disruption of cell coexpression networks in mouse brain.

Rare heterozygous coding variants in the triggering receptor expressed in myeloid cells 2 (TREM2) gene, conferring increased risk of developing late-onset Alzheimer's disease, have been identified. We examined the transcriptional consequences of the loss of Trem2 in mouse brain to better understand its role in disease using differential expression and coexpression network analysis of Trem2 knockout and wild-type mice. We generated RNA-Seq data from cortex and hippocampus sampled at 4 and 8 months. Using brain cell-type markers and ontology enrichment, we found subnetworks with cell type and/or functional identity. We primarily discovered changes in an endothelial gene-enriched subnetwork at 4 months, including a shift toward a more central role for the amyloid precursor protein gene, coupled with widespread disruption of other cell-type subnetworks, including a subnetwork with neuronal identity. We reveal an unexpected potential role of Trem2 in the homeostasis of endothelial cells that goes beyond its known functions as a microglial receptor and signaling hub, suggesting an underlying link between immune response and vascular disease in dementia.

Early-life nutrition modulates the epigenetic state of specific rDNA genetic variants in mice.

A suboptimal early-life environment, due to poor nutrition or stress during pregnancy, can influence lifelong phenotypes in the progeny. Epigenetic factors are thought to be key mediators of these effects. We show that protein restriction in mice from conception until weaning induces a linear correlation between growth restriction and DNA methylation at ribosomal DNA (rDNA). This epigenetic response remains into adulthood and is restricted to rDNA copies associated with a specific genetic variant within the promoter. Related effects are also found in models of maternal high-fat or obesogenic diets. Our work identifies environmentally induced epigenetic dynamics that are dependent on underlying genetic variation and establishes rDNA as a genomic target of nutritional insults.