The multi-faceted regulation of nuclear tRNA gene transcription

Locus-specific proteome decoding reveals Fpt1 as a chromatin-associated negative regulator of RNA polymerase III assembly

Transcription of tRNA genes by RNA polymerase III (RNAPIII) is tuned by signaling cascades. The emerging notion of differential tRNA gene regulation implies the existence of additional regulatory mechanisms. However, tRNA gene-specific regulators have not been described. Decoding the local chromatin proteome of a native tRNA gene in yeast revealed reprogramming of the RNAPIII transcription machinery upon nutrient perturbation. Among the dynamic proteins, we identified Fpt1, a protein of unknown function that uniquely occupied RNAPIII-regulated genes. Fpt1 binding at tRNA genes correlated with the efficiency of RNAPIII eviction upon nutrient perturbation and required the transcription factors TFIIIB and TFIIIC but not RNAPIII. In the absence of Fpt1, eviction of RNAPIII was reduced, and the shutdown of ribosome biogenesis genes was impaired upon nutrient perturbation. Our findings provide support for a chromatin-associated mechanism required for RNAPIII eviction from tRNA genes and tuning the physiological response to changing metabolic demands.

Epigenetic regulation of functional candidate genes for milk production traits in dairy sheep subjected to protein restriction in the prepubertal stage

BACKGROUND

As the prepubertal stage is a crucial point for the proper development of the mammary gland and milk production, this study aims to evaluate how protein restriction at this stage can affect methylation marks in milk somatic cells. Here, 28 Assaf ewes were subjected to 42.3% nutritional protein restriction (14 animals, NPR) or fed standard diets (14 animals, C) during the prepubertal stage. During the second lactation, the milk somatic cells of these ewes were sampled, and the extracted DNA was subjected to whole-genome bisulfite sequencing.

RESULTS

A total of 1154 differentially methylated regions (DMRs) were identified between the NPR and C groups. Indeed, the results of functional enrichment analyses of the genes harboring these DMRs suggested their relevant effects on the development of the mammary gland and lipid metabolism in sheep. The additional analysis of the correlations of the mean methylation levels within these DMRs with fat, protein, and dry extract percentages in the milk and milk somatic cell counts suggested associations between several DMRs and milk production traits. However, there were no phenotypic differences in these traits between the NPR and C groups.

CONCLUSION

In light of the above, the results obtained in the current study might suggest potential candidate genes for the regulation of milk production traits in the sheep mammary gland. Further studies focusing on elucidating the genetic mechanisms affected by the identified DMRs may help to better understand the biological mechanisms modified in the mammary gland of dairy sheep as a response to nutritional challenges and their potential effects on milk production.

Deciphering D4Z4 CpG methylation gradients in fascioscapulohumeral muscular dystrophy using nanopore sequencing

Fascioscapulohumeral muscular dystrophy (FSHD) is caused by a unique genetic mechanism that relies on contraction and hypomethylation of the D4Z4 macrosatellite array on the Chromosome 4q telomere allowing ectopic expression of the DUX4 gene in skeletal muscle. Genetic analysis is difficult because of the large size and repetitive nature of the array, a nearly identical array on the 10q telomere, and the presence of divergent D4Z4 arrays scattered throughout the genome. Here, we combine nanopore long-read sequencing with Cas9-targeted enrichment of 4q and 10q D4Z4 arrays for comprehensive genetic analysis including determination of the length of the 4q and 10q D4Z4 arrays with base-pair resolution. In the same assay, we differentiate 4q from 10q telomeric sequences, determine A/B haplotype, identify paralogous D4Z4 sequences elsewhere in the genome, and estimate methylation for all CpGs in the array. Asymmetric, length-dependent methylation gradients were observed in the 4q and 10q D4Z4 arrays that reach a hypermethylation point at approximately 10 D4Z4 repeat units, consistent with the known threshold of pathogenic D4Z4 contractions. High resolution analysis of individual D4Z4 repeat methylation revealed areas of low methylation near the CTCF/insulator region and areas of high methylation immediately preceding the DUX4 transcriptional start site. Within the DUX4 exons, we observed a waxing/waning methylation pattern with a 180-nucleotide periodicity, consistent with phased nucleosomes. Targeted nanopore sequencing complements recently developed molecular combing and optical mapping approaches to genetic analysis for FSHD by adding precision of the length measurement, base-pair resolution sequencing, and quantitative methylation analysis.

Deciphering D4Z4 CpG methylation gradients in fascioscapulohumeral muscular dystrophy using nanopore sequencing

Fascioscapulohumeral muscular dystrophy (FSHD) is caused by a unique genetic mechanism that relies on contraction and hypomethylation of the D4Z4 macrosatellite array on the chromosome 4q telomere allowing ectopic expression of the DUX4 gene in skeletal muscle. Genetic analysis is difficult due to the large size and repetitive nature of the array, a nearly identical array on the 10q telomere, and the presence of divergent D4Z4 arrays scattered throughout the genome. Here, we combine nanopore long-read sequencing with Cas9-targeted enrichment of 4q and 10q D4Z4 arrays for comprehensive genetic analysis including determination of the length of the 4q and 10q D4Z4 arrays with base-pair resolution. In the same assay, we differentiate 4q from 10q telomeric sequences, determine A/B haplotype, identify paralogous D4Z4 sequences elsewhere in the genome, and estimate methylation for all CpGs in the array. Asymmetric, length-dependent methylation gradients were observed in the 4q and 10q D4Z4 arrays that reach a hypermethylation point at approximately 10 D4Z4 repeat units, consistent with the known threshold of pathogenic D4Z4 contractions. High resolution analysis of individual D4Z4 repeat methylation revealed areas of low methylation near the CTCF/insulator region and areas of high methylation immediately preceding the DUX4 transcriptional start site. Within the DUX4 exons, we observed a waxing/waning methylation pattern with a 180-nucleotide periodicity, consistent with phased nucleosomes. Targeted nanopore sequencing complements recently developed molecular combing and optical mapping approaches to genetic analysis for FSHD by adding precision of the length measurement, base-pair resolution sequencing, and quantitative methylation analysis.

PlantRNA 2.0: an updated database dedicated to tRNAs of photosynthetic eukaryotes

PlantRNA (http://plantrna.ibmp.cnrs.fr/) is a comprehensive database of transfer RNA (tRNA) gene sequences retrieved from fully annotated nuclear, plastidial and mitochondrial genomes of photosynthetic organisms. In the first release (PlantRNA 1.0), tRNA genes from 11 organisms were annotated. In this second version, the annotation was implemented to 51 photosynthetic species covering the whole phylogenetic tree of photosynthetic organisms, from the most basal group of Archeplastida, the glaucophyte Cyanophora paradoxa, to various land plants. tRNA genes from lower photosynthetic organisms such as streptophyte algae or lycophytes as well as extremophile photosynthetic species such as Eutrema parvulum were incorporated in the database. As a whole, about 37 000 tRNA genes were accurately annotated. In the frame of the tRNA genes annotation from the genome of the Rhodophyte Chondrus crispus, non-canonical splicing sites in the D- or T-regions of tRNA molecules were identified and experimentally validated. As for PlantRNA 1.0, comprehensive biological information including 5'- and 3'-flanking sequences, A and B box sequences, region of transcription initiation and poly(T) transcription termination stretches, tRNA intron sequences and tRNA mitochondrial import are included.

Organization and epigenomic control of RNA polymerase III-transcribed genes in plants

The genetic information linearly scripted in chromosomes is wrapped in a ribonucleoprotein complex called chromatin. The adaptation of its compaction level and spatiotemporal organization refines gene expression in response to developmental and environmental cues. RNA polymerase III (RNAPIII) is responsible for the biogenesis of elementary non-coding RNAs. Their genes are subjected to high duplication and mutational rates, and invade nuclear genomes. Their insertion into different epigenomic environments raises the question of how chromatin packing affects their individual transcription. In this review, we provide a unique perspective to this issue in plants. In addition, we discuss how the genomic organization of RNAPIII-transcribed loci, combined with epigenetic differences, might participate to plant trait variations.

PltRNAdb: Plant transfer RNA database

Transfer RNAs (tRNAs) are intermediate-sized non-coding RNAs found in all organisms that help translate messenger RNA into protein. Recently, the number of sequenced plant genomes has increased dramatically. The availability of this extensive data greatly accelerates the study of tRNAs on a large scale. Here, 8,768,261 scaffolds/chromosomes containing 229,093 giga-base pairs representing whole-genome sequences of 256 plant species were analyzed to identify tRNA genes. As a result, 331,242 nuclear, 3,216 chloroplast, and 1,467 mitochondrial tRNA genes were identified. The nuclear tRNA genes include 275,134 tRNAs decoding 20 standard amino acids, 1,325 suppressor tRNAs, 6,273 tRNAs with unknown isotypes, 48,475 predicted pseudogenes, and 37,873 tRNAs with introns. Efforts also extended to the creation of PltRNAdb (https://bioinformatics.um6p.ma/PltRNAdb/index.php), a data source for tRNA genes from 256 plant species. PltRNAdb website allows researchers to search, browse, visualize, BLAST, and download predicted tRNA genes. PltRNAdb will help improve our understanding of plant tRNAs and open the door to discovering the unknown regulatory roles of tRNAs in plant genomes.

Targeted genome editing of plants and plant cells for biomanufacturing

Plants have provided humans with useful products since antiquity, but in the last 30 years they have also been developed as production platforms for small molecules and recombinant proteins. This initially niche area has blossomed with the growth of the global bioeconomy, and now includes chemical building blocks, polymers and renewable energy. All these applications can be described as “plant molecular farming” (PMF). Despite its potential to increase the sustainability of biologics manufacturing, PMF has yet to be embraced broadly by industry. This reflects a combination of regulatory uncertainty, limited information on process cost structures, and the absence of trained staff and suitable manufacturing capacity. However, the limited adaptation of plants and plant cells to the requirements of industry-scale manufacturing is an equally important hurdle. For example, the targeted genetic manipulation of yeast has been common practice since the 1980s, whereas reliable site-directed mutagenesis in most plants has only become available with the advent of CRISPR/Cas9 and similar genome editing technologies since around 2010. Here we summarize the applications of new genetic engineering technologies to improve plants as biomanufacturing platforms. We start by identifying current bottlenecks in manufacturing, then illustrate the progress that has already been made and discuss the potential for improvement at the molecular, cellular and organism levels. We discuss the effects of metabolic optimization, adaptation of the endomembrane system, modified glycosylation profiles, programmable growth and senescence, protease inactivation, and the expression of enzymes that promote biodegradation. We outline strategies to achieve these modifications by targeted gene modification, considering case-by-case examples of individual improvements and the combined modifications needed to generate a new general-purpose “chassis” for PMF.

RNA polymerase III and antiviral innate immune response

The innate immune system has numerous signal transduction pathways that lead to the production of type I interferons in response to exposure of cells to external stimuli. One of these pathways comprises RNA polymerase (Pol) III that senses common DNA viruses, such as cytomegalovirus, vaccinia, herpes simplex virus-1 and varicella zoster virus. This polymerase detects and transcribes viral genomic regions to generate AU-rich transcripts that bring to the induction of type I interferons. Remarkably, Pol III is also stimulated by foreign non-viral DNAs and expression of one of its subunits is induced by an RNA virus, the Sindbis virus. Moreover, a protein subunit of RNase P, which is known to associate with Pol III in initiation complexes, is induced by viral infection. Accordingly, alliance of the two tRNA enzymes in innate immunity merits a consideration.

Combining tRNA sequencing methods to characterize plant tRNA expression and post-transcriptional modification

Differences in tRNA expression have been implicated in a remarkable number of biological processes. There is growing evidence that tRNA genes can play dramatically different roles depending on both expression and post-transcriptional modification, yet sequencing tRNAs to measure abundance and detect modifications remains challenging. Their secondary structure and extensive post-transcriptional modifications interfere with RNA-seq library preparation methods and have limited the utility of high-throughput sequencing technologies. Here, we combine two modifications to standard RNA-seq methods by treating with the demethylating enzyme AlkB and ligating with tRNA-specific adapters in order to sequence tRNAs from four species of flowering plants, a group that has been shown to have some of the most extensive rates of post-transcriptional tRNA modifications. This protocol has the advantage of detecting full-length tRNAs and sequence variants that can be used to infer many post-transcriptional modifications. We used the resulting data to produce a modification index of almost all unique reference tRNAs in Arabidopsis thaliana, which exhibited many anciently conserved similarities with humans but also positions that appear to be 'hot spots' for modifications in angiosperm tRNAs. We also found evidence based on northern blot analysis and droplet digital PCR that, even after demethylation treatment, tRNA-seq can produce highly biased estimates of absolute expression levels most likely due to biased reverse transcription. Nevertheless, the generation of full-length tRNA sequences with modification data is still promising for assessing differences in relative tRNA expression across treatments, tissues or subcellular fractions and help elucidate the functional roles of tRNA modifications.

Epigenetic silencing of clustered tRNA genes in Arabidopsis

Beyond their key role in translation, cytosolic transfer RNAs (tRNAs) are involved in a wide range of other biological processes. Nuclear tRNA genes (tDNAs) are transcribed by the RNA polymerase III (RNAP III) and cis-elements, trans-factors as well as genomic features are known to influence their expression. In Arabidopsis, besides a predominant population of dispersed tDNAs spread along the 5 chromosomes, some clustered tDNAs have been identified. Here, we demonstrate that these tDNA clusters are transcriptionally silent and that pathways involved in the maintenance of DNA methylation play a predominant role in their repression. Moreover, we show that clustered tDNAs exhibit repressive chromatin features whilst their dispersed counterparts contain permissive euchromatic marks. This work demonstrates that both genomic and epigenomic contexts are key players in the regulation of tDNAs transcription. The conservation of most of these regulatory processes suggests that this pioneering work in Arabidopsis can provide new insights into the regulation of RNA Pol III transcription in other organisms, including vertebrates.

The Target of Rapamycin Signalling Pathway in Ageing and Lifespan Regulation

Ageing is a complex trait controlled by genes and the environment. The highly conserved mechanistic target of rapamycin signalling pathway (mTOR) is a major regulator of lifespan in all eukaryotes and is thought to be mediating some of the effects of dietary restriction. mTOR is a rheostat of energy sensing diverse inputs such as amino acids, oxygen, hormones, and stress and regulates lifespan by tuning cellular functions such as gene expression, ribosome biogenesis, proteostasis, and mitochondrial metabolism. Deregulation of the mTOR signalling pathway is implicated in multiple age-related diseases such as cancer, neurodegeneration, and auto-immunity. In this review, we briefly summarise some of the workings of mTOR in lifespan and ageing through the processes of transcription, translation, autophagy, and metabolism. A good understanding of the pathway's outputs and connectivity is paramount towards our ability for genetic and pharmacological interventions for healthy ageing and amelioration of age-related disease.

Interchangeable parts: The evolutionarily dynamic tRNA population in plant mitochondria

Transfer RNAs (tRNAs) remain one of the very few classes of genes still encoded in the mitochondrial genome. These key components of the protein translation system must interact with a large enzymatic network of nuclear-encoded gene products to maintain mitochondrial function. Plants have an evolutionarily dynamic mitochondrial tRNA population, including ongoing tRNA gene loss and replacement by both horizontal gene transfer from diverse sources and import of nuclear-expressed tRNAs from the cytosol. Thus, plant mitochondria represent an excellent model for understanding how anciently divergent genes can act as "interchangeable parts" during the evolution of complex molecular systems. In particular, understanding the integration of the mitochondrial translation system with elements of the corresponding machinery used in cytosolic protein synthesis is a key area for eukaryotic cellular evolution. Here, we review the increasingly detailed phylogenetic data about the evolutionary history of mitochondrial tRNA gene loss, transfer, and functional replacement that has created extreme variation in mitochondrial tRNA populations across plant species. We describe emerging tRNA-seq methods with promise for refining our understanding of the expression and subcellular localization of tRNAs. Finally, we summarize current evidence and identify open questions related to coevolutionary changes in nuclear-encoded enzymes that have accompanied turnover in mitochondrial tRNA populations.

The GATA Transcription Factor Gaf1 Represses tRNAs, Inhibits Growth, and Extends Chronological Lifespan Downstream of Fission Yeast TORC1

Target of Rapamycin Complex 1 (TORC1) signaling promotes growth and aging. Inhibition of TORC1 leads to reduced protein translation, which promotes longevity. TORC1-dependent post-transcriptional regulation of protein translation has been well studied, while analogous transcriptional regulation is less understood. Here we screen fission yeast mutants for resistance to Torin1, which inhibits TORC1 and cell growth. Cells lacking the GATA factor Gaf1 (gaf1Δ) grow normally even in high doses of Torin1. The gaf1Δ mutation shortens the chronological lifespan of non-dividing cells and diminishes Torin1-mediated longevity. Expression profiling and genome-wide binding experiments show that upon TORC1 inhibition, Gaf1 directly upregulates genes for small-molecule metabolic pathways and indirectly represses genes for protein translation. Surprisingly, Gaf1 binds to and downregulates the tRNA genes, so it also functions as a transcription factor for RNA polymerase III. Thus, Gaf1 controls the transcription of both protein-coding and tRNA genes to inhibit translation and growth downstream of TORC1.

Predicting transfer RNA gene activity from sequence and genome context

Transfer RNA (tRNA) genes are among the most highly transcribed genes in the genome owing to their central role in protein synthesis. However, there is evidence for a broad range of gene expression across tRNA loci. This complexity, combined with difficulty in measuring transcript abundance and high sequence identity across transcripts, has severely limited our collective understanding of tRNA gene expression regulation and evolution. We establish sequence-based correlates to tRNA gene expression and develop a tRNA gene classification method that does not require, but benefits from, comparative genomic information and achieves accuracy comparable to molecular assays. We observe that guanine + cytosine (G + C) content and CpG density surrounding tRNA loci is exceptionally well correlated with tRNA gene activity, supporting a prominent regulatory role of the local genomic context in combination with internal sequence features. We use our tRNA gene activity predictions in conjunction with a comprehensive tRNA gene ortholog set spanning 29 placental mammals to estimate the evolutionary rate of functional changes among orthologs. Our method adds a new dimension to large-scale tRNA functional prediction and will help prioritize characterization of functional tRNA variants. Its simplicity and robustness should enable development of similar approaches for other clades, as well as exploration of functional diversification of members of large gene families.