Mitochondrial point heteroplasmy: insights from deep-sequencing of human replicate samples.

BACKGROUND: Human mitochondrial heteroplasmy is an extensively investigated phenomenon in the context of medical diagnostics, forensic identification and molecular evolution. However, technical limitations of high-throughput sequencing hinder reliable determination of point heteroplasmies (PHPs) with minor allele frequencies (MAFs) within the noise threshold. RESULTS: To investigate the PHP landscape at an MAF threshold down to 0.1%, we sequenced whole mitochondrial genomes at approximately 7.700x coverage, in multiple technical and biological replicates of longitudinal blood and buccal swab samples from 11 human donors (159 libraries in total). The results obtained by two independent sequencing platforms and bioinformatics pipelines indicate distinctive PHP patterns below and above the 1% MAF cut-off. We found a high inter-individual prevalence of low-level PHPs (MAF < 1%) at polymorphic positions of the mitochondrial DNA control region (CR), their tissue preference, and a tissue-specific minor allele linkage. We also established the position-dependent potential of minor allele expansion in PHPs, and short-term PHP instability in a mitotically active tissue. We demonstrate that the increase in sensitivity of PHP detection to minor allele frequencies below 1% within a robust experimental and analytical pipeline, provides new information with potential applicative value. CONCLUSIONS: Our findings reliably show different mutational loads between tissues at sub-1% allele frequencies, which may serve as an informative medical biomarker of time-dependent, tissue-specific mutational burden, or help discriminate forensically relevant tissues in a single person, close maternal relatives or unrelated individuals of similar phylogenetic background.

Transfection of Sponge Cells and Intracellular Localization of Cancer-Related MYC, RRAS2, and DRG1 Proteins.

The determination of the protein's intracellular localization is essential for understanding its biological function. Protein localization studies are mainly performed on primary and secondary vertebrate cell lines for which most protocols have been optimized. In spite of experimental difficulties, studies on invertebrate cells, including basal Metazoa, have greatly advanced. In recent years, the interest in studying human diseases from an evolutionary perspective has significantly increased. Sponges, placed at the base of the animal tree, are simple animals without true tissues and organs but with a complex genome containing many genes whose human homologs have been implicated in human diseases, including cancer. Therefore, sponges are an innovative model for elucidating the fundamental role of the proteins involved in cancer. In this study, we overexpressed human cancer-related proteins and their sponge homologs in human cancer cells, human fibroblasts, and sponge cells. We demonstrated that human and sponge MYC proteins localize in the nucleus, the RRAS2 in the plasma membrane, the membranes of the endolysosomal vesicles, and the DRG1 in the cell's cytosol. Despite the very low transfection efficiency of sponge cells, we observed an identical localization of human proteins and their sponge homologs, indicating their similar cellular functions.

Genomic profiling of HIV-1 integration in microglia cells links viral integration to the topologically associated domains.

HIV-1 encounters the hierarchically organized host chromatin to stably integrate and persist in anatomically distinct latent reservoirs. The contribution of genome organization in HIV-1 infection has been largely understudied across different HIV-1 targets. Here, we determine HIV-1 integration sites (ISs), associate them with chromatin and expression signatures at different genomic scales in a microglia cell model, and profile them together with the primary T cell reservoir. HIV-1 insertions into introns of actively transcribed genes with IS hotspots in genic and super-enhancers, characteristic of blood cells, are maintained in the microglia cell model. Genome organization analysis reveals dynamic CCCTC-binding factor (CTCF) clusters in cells with active and repressed HIV-1 transcription, whereas CTCF removal impairs viral integration. We identify CTCF-enriched topologically associated domain (TAD) boundaries with signatures of transcriptionally active chromatin as HIV-1 integration determinants in microglia and CD4+ T cells, highlighting the importance of host genome organization in HIV-1 infection.

Structure and function of cancer-related developmentally regulated GTP-binding protein 1 (DRG1) is conserved between sponges and humans.

Cancer is a disease caused by errors within the multicellular system and it represents a major health issue in multicellular organisms. Although cancer research has advanced substantially, new approaches focusing on fundamental aspects of cancer origin and mechanisms of spreading are necessary. Comparative genomic studies have shown that most genes linked to human cancer emerged during the early evolution of Metazoa. Thus, basal animals without true tissues and organs, such as sponges (Porifera), might be an innovative model system for understanding the molecular mechanisms of proteins involved in cancer biology. One of these proteins is developmentally regulated GTP-binding protein 1 (DRG1), a GTPase stabilized by interaction with DRG family regulatory protein 1 (DFRP1). This study reveals a high evolutionary conservation of DRG1 gene/protein in metazoans. Our biochemical analysis and structural predictions show that both recombinant sponge and human DRG1 are predominantly monomers that form complexes with DFRP1 and bind non-specifically to RNA and DNA. We demonstrate the conservation of sponge and human DRG1 biological features, including intracellular localization and DRG1:DFRP1 binding, function of DRG1 in α-tubulin dynamics, and its role in cancer biology demonstrated by increased proliferation, migration and colonization in human cancer cells. These results suggest that the ancestor of all Metazoa already possessed DRG1 that is structurally and functionally similar to the human DRG1, even before the development of real tissues or tumors, indicating an important function of DRG1 in fundamental cellular pathways.

Embryo-Like Features in Developing Bacillus subtilis Biofilms.

Correspondence between evolution and development has been discussed for more than two centuries. Recent work reveals that phylogeny-ontogeny correlations are indeed present in developmental transcriptomes of eukaryotic clades with complex multicellularity. Nevertheless, it has been largely ignored that the pervasive presence of phylogeny-ontogeny correlations is a hallmark of development in eukaryotes. This perspective opens a possibility to look for similar parallelisms in biological settings where developmental logic and multicellular complexity are more obscure. For instance, it has been increasingly recognized that multicellular behavior underlies biofilm formation in bacteria. However, it remains unclear whether bacterial biofilm growth shares some basic principles with development in complex eukaryotes. Here we show that the ontogeny of growing Bacillus subtilis biofilms recapitulates phylogeny at the expression level. Using time-resolved transcriptome and proteome profiles, we found that biofilm ontogeny correlates with the evolutionary measures, in a way that evolutionary younger and more diverged genes were increasingly expressed toward later timepoints of biofilm growth. Molecular and morphological signatures also revealed that biofilm growth is highly regulated and organized into discrete ontogenetic stages, analogous to those of eukaryotic embryos. Together, this suggests that biofilm formation in Bacillus is a bona fide developmental process comparable to organismal development in animals, plants, and fungi. Given that most cells on Earth reside in the form of biofilms and that biofilms represent the oldest known fossils, we anticipate that the widely adopted vision of the first life as a single-cell and free-living organism needs rethinking.

Microbiota Alters Urinary Bladder Weight and Gene Expression.

We studied the effect of microbiota on the transcriptome and weight of the urinary bladder by comparing germ-free (GF) and specific pathogen-free (SPF) housed mice. In total, 97 genes were differently expressed (fold change > ±2; false discovery rate (FDR) p-value < 0.01) between the groups, including genes regulating circadian rhythm (Per1, Per2 and Per3), extracellular matrix (Spo1, Spon2), and neuromuscular synaptic transmission (Slc18a3, Slc5a7, Chrnb4, Chrna3, Snap25). The highest increase in expression was observed for immunoglobulin genes (Igkv1-122, Igkv4-68) of unknown function, but surprisingly the absence of microbiota did not change the expression of the genes responsible for recognizing microbes and their products. We found that urinary bladder weight was approximately 25% lighter in GF mice (p = 0.09 for males, p = 0.005 for females) and in mice treated with broad spectrum of antibiotics (p = 0.0002). In conclusion, our data indicate that microbiota is an important determinant of urinary bladder physiology controlling its gene expression and size.

Evolutionary Analysis of the Bacillus subtilis Genome Reveals New Genes Involved in Sporulation.

Bacilli can form dormant, highly resistant, and metabolically inactive spores to cope with extreme environmental challenges. In this study, we examined the evolutionary age of Bacillus subtilis sporulation genes using the approach known as genomic phylostratigraphy. We found that B. subtilis sporulation genes cluster in several groups that emerged at distant evolutionary time-points, suggesting that the sporulation process underwent several stages of expansion. Next, we asked whether such evolutionary stratification of the genome could be used to predict involvement in sporulation of presently uncharacterized genes (y-genes). We individually inactivated a representative sample of uncharacterized genes that arose during the same evolutionary periods as the known sporulation genes and tested the resulting strains for sporulation phenotypes. Sporulation was significantly affected in 16 out of 37 (43%) tested strains. In addition to expanding the knowledge base on B. subtilis sporulation, our findings suggest that evolutionary age could be used to help with genome mining.

The most abundant maternal lncRNA Sirena1 acts post-transcriptionally and impacts mitochondrial distribution.

Tens of thousands of rapidly evolving long non-coding RNA (lncRNA) genes have been identified, but functions were assigned to relatively few of them. The lncRNA contribution to the mouse oocyte physiology remains unknown. We report the evolutionary history and functional analysis of Sirena1, the most expressed lncRNA and the 10th most abundant poly(A) transcript in mouse oocytes. Sirena1 appeared in the common ancestor of mouse and rat and became engaged in two different post-transcriptional regulations. First, antisense oriented Elob pseudogene insertion into Sirena1 exon 1 is a source of small RNAs targeting Elob mRNA via RNA interference. Second, Sirena1 evolved functional cytoplasmic polyadenylation elements, an unexpected feature borrowed from translation control of specific maternal mRNAs. Sirena1 knock-out does not affect fertility, but causes minor dysregulation of the maternal transcriptome. This includes increased levels of Elob and mitochondrial mRNAs. Mitochondria in Sirena1-/- oocytes disperse from the perinuclear compartment, but do not change in number or ultrastructure. Taken together, Sirena1 contributes to RNA interference and mitochondrial aggregation in mouse oocytes. Sirena1 exemplifies how lncRNAs stochastically engage or even repurpose molecular mechanisms during evolution. Simultaneously, Sirena1 expression levels and unique functional features contrast with the lack of functional importance assessed under laboratory conditions.

Spatially clustered loci with multiple enhancers are frequent targets of HIV-1 integration.

HIV-1 recurrently targets active genes and integrates in the proximity of the nuclear pore compartment in CD4+ T cells. However, the genomic features of these genes and the relevance of their transcriptional activity for HIV-1 integration have so far remained unclear. Here we show that recurrently targeted genes are proximal to super-enhancer genomic elements and that they cluster in specific spatial compartments of the T cell nucleus. We further show that these gene clusters acquire their location during the activation of T cells. The clustering of these genes along with their transcriptional activity are the major determinants of HIV-1 integration in T cells. Our results provide evidence of the relevance of the spatial compartmentalization of the genome for HIV-1 integration, thus further strengthening the role of nuclear architecture in viral infection.

Coupling AAV-mediated promoterless gene targeting to SaCas9 nuclease to efficiently correct liver metabolic diseases.

Non-integrative AAV-mediated gene therapy in the liver is effective in adult patients, but faces limitations in pediatric settings due to episomal DNA loss during hepatocyte proliferation. Gene targeting is a promising approach by permanently modifying the genome. We previously rescued neonatal lethality in Crigler-Najjar mice by inserting a promoterless human uridine glucuronosyl transferase A1 (UGT1A1) cDNA in exon 14 of the albumin gene, without the use of nucleases. To increase recombination rate and therapeutic efficacy, here we used CRISPR/SaCas9. Neonatal mice were transduced with two AAVs: one expressing the SaCas9 and sgRNA, and one containing a promoterless cDNA flanked by albumin homology regions. Targeting efficiency increased ~26-fold with an eGFP reporter cDNA, reaching up to 24% of eGFP-positive hepatocytes. Next, we fully corrected the diseased phenotype of Crigler-Najjar mice by targeting the hUGT1A1 cDNA. Treated mice had normal plasma bilirubin up to 10 months after administration, hUGT1A1 protein levels were ~6-fold higher than in WT liver, with a 90-fold increase in recombination rate. Liver histology, inflammatory markers, and plasma albumin were normal in treated mice, with no off-targets in predicted sites. Thus, the improved efficacy and reassuring safety profile support the potential application of the proposed approach to other liver diseases.

Dynamics of exhaled breath temperature after smoking a cigarette and its association with lung function changes predictive of COPD risk in smokers: a cross-sectional study.

Exhaled breath temperature (EBT) is a biomarker of inflammation and vascularity of the airways already shown to predict incident COPD. This cross-sectional study was aimed to assess the potential of EBT in identifying "healthy" smokers susceptible to cigarette smoke toxicity of the airways and to the risk of developing COPD by analysing the dynamics of EBT after smoking a cigarette and its associations with their demographics (age, smoking burden) and lung function. The study included 55 current smokers of both sexes, 29-62 years of age, with median smoking exposure of 15 (10-71.8) pack-years. EBT was measured at baseline and 5, 15, 30, 45, and 60 min after smoking a single cigarette. Lung function was measured with spirometry followed by a bronchodilator test. To compare changes in EBT between repeated measurements we used the analysis of variance and the area under the curve (EBTAUC) as a dependent variable. Multivariate regression analysis was used to look for associations with patient characteristics and lung function in particular. The average (±SD) baseline EBT was 33.42±1.50 °C. The highest significant increase to 33.84 (1.25) °C was recorded 5 min after the cigarette was smoked (p=0.003), and it took one hour for it to return to the baseline. EBTAUC showed significant repeatability (ICC=0.85, p<0.001) and was significantly associated with age, body mass index, number of cigarettes smoked a day, baseline EBT, and baseline FEF75 (R2=0.39, p<0.001 for the model). Our results suggest that EBT after smoking a single cigarette could be used as early risk predictor of changes associated with chronic cigarette smoke exposure.

Role of Cnot6l in maternal mRNA turnover.

Removal of poly(A) tail is an important mechanism controlling eukaryotic mRNA turnover. The major eukaryotic deadenylase complex CCR4-NOT contains two deadenylase components, CCR4 and CAF1, for which mammalian CCR4 is encoded by Cnot6 or Cnot6l paralogs. We show that Cnot6l apparently supplies the majority of CCR4 in the maternal CCR4-NOT in mouse, hamster, and bovine oocytes. Deletion of Cnot6l yielded viable mice, but Cnot6l -/- females exhibited ∼40% smaller litter size. The main onset of the phenotype was post-zygotic: fertilized Cnot6l -/- eggs developed slower and arrested more frequently than Cnot6l +/- eggs, suggesting that maternal CNOT6L is necessary for accurate oocyte-to-embryo transition. Transcriptome analysis revealed major transcriptome changes in Cnot6l -/- ovulated eggs and one-cell zygotes. In contrast, minimal transcriptome changes in preovulatory Cnot6l -/- oocytes were consistent with reported Cnot6l mRNA dormancy. A minimal overlap between transcripts sensitive to decapping inhibition and Cnot6l loss suggests that decapping and CNOT6L-mediated deadenylation selectively target distinct subsets of mRNAs during oocyte-to-embryo transition in mouse.

Heat-induced longevity in budding yeast requires respiratory metabolism and glutathione recycling.

Heat-induced hormesis is a well-known conserved phenomenon in aging, traditionally attributed to the benefits conferred by increased amounts of heat shock (HS) proteins. Here we find that the key event for the HS-induced lifespan extension in budding yeast is the switch from glycolysis to respiratory metabolism. The resulting increase in reactive oxygen species activates the antioxidant response, supported by the redirection of glucose from glycolysis to the pentose phosphate pathway, increasing the production of NADPH. This sequence of events culminates in replicative lifespan (RLS) extension, implying decreased mortality per generation that persists even after the HS has finished. We found that switching to respiratory metabolism, and particularly the consequent increase in glutathione levels, were essential for the observed RLS extension. These results draw the focus away solely from the HS response and demonstrate that the antioxidant response has a key role in heat-induced hormesis. Our findings underscore the importance of the changes in cellular metabolic activity for heat-induced longevity in budding yeast.

TORC1-mediated sensing of chaperone activity alters glucose metabolism and extends lifespan.

Protein quality control mechanisms, required for normal cellular functioning, encompass multiple functions related to protein production and maintenance. However, the existence of communication between proteostasis and metabolic networks and its underlying mechanisms remain elusive. Here, we report that enhanced chaperone activity and consequent improved proteostasis are sensed by TORC1 via the activity of Hsp82. Chaperone enrichment decreases the level of Hsp82, which deactivates TORC1 and leads to activation of Snf1/AMPK, regardless of glucose availability. This mechanism culminates in the extension of yeast replicative lifespan (RLS) that is fully reliant on both TORC1 deactivation and Snf1/AMPK activation. Specifically, we identify oxygen consumption increase as the downstream effect of Snf1 activation responsible for the entire RLS extension. Our results set a novel paradigm for the role of proteostasis in aging: modulation of the misfolded protein level can affect cellular metabolic features as well as mitochondrial activity and consequently modify lifespan. The described mechanism is expected to open new avenues for research of aging and age-related diseases.

Long terminal repeats power evolution of genes and gene expression programs in mammalian oocytes and zygotes.

Retrotransposons are "copy-and-paste" insertional mutagens that substantially contribute to mammalian genome content. Retrotransposons often carry long terminal repeats (LTRs) for retrovirus-like reverse transcription and integration into the genome. We report an extraordinary impact of a group of LTRs from the mammalian endogenous retrovirus-related ERVL retrotransposon class on gene expression in the germline and beyond. In mouse, we identified more than 800 LTRs from ORR1, MT, MT2, and MLT families, which resemble mobile gene-remodeling platforms that supply promoters and first exons. The LTR-mediated gene remodeling also extends to hamster, human, and bovine oocytes. The LTRs function in a stage-specific manner during the oocyte-to-embryo transition by activating transcription, altering protein-coding sequences, producing noncoding RNAs, and even supporting evolution of new protein-coding genes. These functions result, for example, in recycling processed pseudogenes into mRNAs or lncRNAs with regulatory roles. The functional potential of the studied LTRs is even higher, because we show that dormant LTR promoter activity can rescue loss of an essential upstream promoter. We also report a novel protein-coding gene evolution-D6Ertd527e-in which an MT LTR provided a promoter and the 5' exon with a functional start codon while the bulk of the protein-coding sequence evolved through a CAG repeat expansion. Altogether, ERVL LTRs provide molecular mechanisms for stochastically scanning, rewiring, and recycling genetic information on an extraordinary scale. ERVL LTRs thus offer means for a comprehensive survey of the genome's expression potential, tightly intertwining with gene expression and evolution in the germline.

Long non-coding RNA exchange during the oocyte-to-embryo transition in mice.

The oocyte-to-embryo transition (OET) transforms a differentiated gamete into pluripotent blastomeres. The accompanying maternal-zygotic RNA exchange involves remodeling of the long non-coding RNA (lncRNA) pool. Here, we used next generation sequencing and de novo transcript assembly to define the core population of 1,600 lncRNAs expressed during the OET (lncRNAs). Relative to mRNAs, OET lncRNAs were less expressed and had shorter transcripts, mainly due to fewer exons and shorter 5' terminal exons. Approximately half of OET lncRNA promoters originated in retrotransposons suggesting their recent emergence. Except for a small group of ubiquitous lncRNAs, maternal and zygotic lncRNAs formed two distinct populations. The bulk of maternal lncRNAs was degraded before the zygotic genome activation. Interestingly, maternal lncRNAs seemed to undergo cytoplasmic polyadenylation observed for dormant mRNAs. We also identified lncRNAs giving rise to trans-acting short interfering RNAs, which represent a novel lncRNA category. Altogether, we defined the core OET lncRNA transcriptome and characterized its remodeling during early development. Our results are consistent with the notion that rapidly evolving lncRNAs constitute signatures of cells-of-origin while a minority plays an active role in control of gene expression across OET. Our data presented here provide an excellent source for further OET lncRNA studies.

Big Data, Evolution, and Metagenomes: Predicting Disease from Gut Microbiota Codon Usage Profiles.

Metagenomics projects use next-generation sequencing to unravel genetic potential in microbial communities from a wealth of environmental niches, including those associated with human body and relevant to human health. In order to understand large datasets collected in metagenomics surveys and interpret them in context of how a community metabolism as a whole adapts and interacts with the environment, it is necessary to extend beyond the conventional approaches of decomposing metagenomes into microbial species' constituents and performing analysis on separate components. By applying concepts of translational optimization through codon usage adaptation on entire metagenomic datasets, we demonstrate that a bias in codon usage present throughout the entire microbial community can be used as a powerful analytical tool to predict for community lifestyle-specific metabolism. Here we demonstrate this approach combined with machine learning, to classify human gut microbiome samples according to the pathological condition diagnosed in the human host.