Unraveling the association between mRNA expressions and mutant phenotypes in a genome-wide assessment of mice

Reduced NET1 adversely affects early embryonic development in mice

Neuroepithelial transforming gene 1 (NET1) is a RhoA subfamily-specific guanine/nucleotide-exchange factor that exhibits critical roles in diverse biological processes. However, the functions in mouse preimplantation embryonic development have not yet determined. In the present study we demonstrated that NET1 is a key factor in the outcome of early mouse embryonic development. Immunofluorescence detection showed that NET1 is principally localized to the nucleus during mouse pre-implantation embryonic development. Silencing Net1 at the zygote stage using a specific siRNA impaired the developmental competence of early mouse embryos, and Net1-knockdown (Net1-KD) induced mitotic spindle-assembly defects and chromosomal alignment abnormalities at the first embryonic cleavage. In addition, reduced NET1 exacerbated reactive oxygen species production and DNA lesions in two-cell stage embryos, further augmenting cellular apoptosis in the preimplantation blastocyst. In summary, our data display key roles for NET1 in mitotic spindle assembly, oxidative stress, and DNA damage during early mouse embryonic development.

Degeneration of the Olfactory System in a Murid Rodent that Evolved Diurnalism

In mammalian research, it has been debated what can initiate an evolutionary tradeoff between different senses, and the phenomenon of sensory tradeoff in rodents, the most abundant mammalian clade, is not evident. The Nile rat (Arvicanthis niloticus), a murid rodent, recently adapted to a diurnal niche through an evolutionary acquisition of daylight vision with enhanced visual acuity. As such, this model provides an opportunity for a cross-species investigation where comparative morphological and multi-omic analyses of the Nile rat are made with its closely related nocturnal species, e.g. the mouse (Mus musculus) and the rat (Rattus norvegicus). Thus, morphological examinations were performed, and evolutionary reductions in relative sizes of turbinal bone surfaces, the cribriform plate, and the olfactory bulb were discovered in Nile rats. Subsequently, we compared multiple murid genomes, and profiled olfactory epithelium transcriptomes of mice and Nile rats at various ages with RNA sequencing. The results further demonstrate that, in comparison with mouse olfactory receptor (OR) genes, Nile rat OR genes have experienced less frequent gain, more frequent loss, and more frequent expression reduction during their evolution. Furthermore, functional degeneration of coding sequences in the Nile rat lineage was found in OR genes, yet not in other genes. Taken together, these results suggest that acquisition of improved vision in the Nile rat has been accompanied by degeneration of both olfaction-related anatomical structures and OR gene repertoires, consistent with the hypothesis of an olfaction-vision tradeoff initiated by the switch from a nocturnal to a diurnal lifestyle in mammals.

The Bacillus subtilis ywbD gene encodes RlmQ, the 23S rRNA methyltransferase forming m7G2574 in the A-site of the peptidyl transferase center

Ribosomal RNA contains many posttranscriptionally modified nucleosides, particularly in the functional parts of the ribosome. The distribution of these modifications varies from one organism to another. In Bacillus subtilis, the model organism for Gram-positive bacteria, mass spectrometry experiments revealed the presence of 7-methylguanosine (m7G) at position 2574 of the 23S rRNA, which lies in the A-site of the peptidyl transferase center of the large ribosomal subunit. Testing several m7G methyltransferase candidates allowed us to identify the RlmQ enzyme, encoded by the ywbD open reading frame, as the MTase responsible for this modification. The enzyme methylates free RNA and not ribosomal 50S or 70S particles, suggesting that modification occurs in the early steps of ribosome biogenesis.

Identification of Two Flip-Over Genes in Grass Family as Potential Signature of C4 Photosynthesis Evolution

In flowering plants, C4 photosynthesis is superior to C3 type in carbon fixation efficiency and adaptation to extreme environmental conditions, but the mechanisms behind the assembly of C4 machinery remain elusive. This study attempts to dissect the evolutionary divergence from C3 to C4 photosynthesis in five photosynthetic model plants from the grass family, using a combined comparative transcriptomics and deep learning technology. By examining and comparing gene expression levels in bundle sheath and mesophyll cells of five model plants, we identified 16 differentially expressed signature genes showing cell-specific expression patterns in C3 and C4 plants. Among them, two showed distinctively opposite cell-specific expression patterns in C3 vs. C4 plants (named as FOGs). The in silico physicochemical analysis of the two FOGs illustrated that C3 homologous proteins of LHCA6 had low and stable pI values of ~6, while the pI values of LHCA6 homologs increased drastically in C4 plants Setaria viridis (7), Zea mays (8), and Sorghum bicolor (over 9), suggesting this protein may have different functions in C3 and C4 plants. Interestingly, based on pairwise protein sequence/structure similarities between each homologous FOG protein, one FOG PGRL1A showed local inconsistency between sequence similarity and structure similarity. To find more examples of the evolutionary characteristics of FOG proteins, we investigated the protein sequence/structure similarities of other FOGs (transcription factors) and found that FOG proteins have diversified incompatibility between sequence and structure similarities during grass family evolution. This raised an interesting question as to whether the sequence similarity is related to structure similarity during C4 photosynthesis evolution.

The Bacillus subtilis open reading frame ysgA encodes the SPOUT methyltransferase RlmP forming 2'-O-methylguanosine at position 2553 in the A-loop of 23S rRNA

A previous bioinformatic analysis predicted that the ysgA open reading frame of Bacillus subtilis encodes an RNA methyltransferase of the SPOUT superfamily. Here we show that YsgA is the 2'-O-methyltransferase that targets position G2553 (Escherichia coli numbering) of the A-loop of 23S rRNA. This was shown by a combination of biochemical and mass spectrometry approaches using both rRNA extracted from B. subtilis wild-type or ΔysgA cells and in vitro synthesized rRNA. When the target G2553 is mutated, YsgA is able to methylate the ribose of adenosine. However, it cannot methylate cytidine nor uridine. The enzyme modifies free 23S rRNA but not the fully assembled ribosome nor the 50S subunit, suggesting that the modification occurs early during ribosome biogenesis. Nevertheless, ribosome subunits assembly is unaffected in a B. subtilis ΔysgA mutant strain. The crystal structure of the recombinant YsgA protein, combined with mutagenesis data, outlined in this article highlights a typical SPOUT fold preceded by an L7Ae/L30 (eL8/eL30 in a new nomenclature) amino-terminal domain.

Transformation to ischaemia tolerance of frog brain function corresponds to dynamic changes in mRNA co-expression across metabolic pathways

Neural activity is costly and requires continuous ATP from aerobic metabolism. Brainstem motor function of American bullfrogs normally collapses after minutes of ischaemia, but following hibernation, it becomes ischaemia-tolerant, generating output for up to 2 h without oxygen or glucose delivery. Transforming the brainstem to function during ischaemia involves a switch to anaerobic glycolysis and brain glycogen. We hypothesized that improving neural performance during ischaemia involves a transcriptional program for glycogen and glucose metabolism. Here we measured mRNA copy number of genes along the path from glycogen metabolism to lactate production using real-time quantitative PCR. The expression of individual genes did not reflect enhanced glucose metabolism. However, the number of co-expressed gene pairs increased early into hibernation, and by the end, most genes involved in glycogen metabolism, glucose transport and glycolysis exhibited striking linear co-expression. By contrast, co-expression of genes in the Krebs cycle and electron transport chain decreased throughout hibernation. Our results uncover reorganization of the metabolic transcriptional network associated with a shift to ischaemia tolerance in brain function. We conclude that modifying gene co-expression may be a critical step in synchronizing storage and use of glucose to achieve ischaemia tolerance in active neural circuits.

Plasticity and conditional essentiality of modification enzymes for domain V of Escherichia coli 23S ribosomal RNA

Escherichia coli rRNAs are post-transcriptionally modified at 36 positions but their modification enzymes are dispensable individually for growth, bringing into question their significance. However, a major growth defect was reported for deletion of the RlmE enzyme, which abolished a 2'O methylation near the peptidyl transferase center (PTC) of the 23S rRNA. Additionally, an adjacent 80-nt "critical region" around the PTC had to be modified to yield significant peptidyl transferase activity in vitro. Surprisingly, we discovered that an absence of just two rRNA modification enzymes is conditionally lethal (at 20°C): RlmE and RluC. At a permissive temperature (37°C), this double knockout was shown to abolish four modifications and be defective in ribosome assembly, though not more so than the RlmE single knockout. However, the double knockout exhibited an even lower rate of tripeptide synthesis than did the single knockout, suggesting an even more defective ribosomal translocation. A combination knockout of the five critical-region-modifying enzymes RluC, RlmKL, RlmN, RlmM, and RluE (not RlmE), which synthesize five of the seven critical-region modifications and 14 rRNA and tRNA modifications altogether, was viable (minor growth defect at 37°C, major at 20°C). This was surprising based on prior in vitro studies. This five-knockout combination had minimal effects on ribosome assembly and frameshifting at 37°C, but greater effects on ribosome assembly and in vitro peptidyl transferase activity at cooler temperatures. These results establish the conditional essentiality of bacterial rRNA modification enzymes and also reveal unexpected plasticity of modification of the PTC region in vivo.

GmKIX8-1 regulates organ size in soybean and is the causative gene for the major seed weight QTL qSw17-1

Seed weight is one of the most important agronomic traits in soybean for yield improvement and food production. Several quantitative trait loci (QTLs) associated with the trait have been identified in soybean. However, the genes underlying the QTLs and their functions remain largely unknown. Using forward genetic methods and CRISPR/Cas9 gene editing, we identified and characterized the role of GmKIX8-1 in the control of organ size in soybean. GmKIX8-1 belongs to a family of KIX domain-containing proteins that negatively regulate cell proliferation in plants. Consistent with this predicted function, we found that loss-of-function GmKIX8-1 mutants showed a significant increase in the size of aerial plant organs, such as seeds and leaves. Likewise, the increase in organ size is due to increased cell proliferation, rather than cell expansion, and increased expression of CYCLIN D3;1-10. Lastly, molecular analysis of soybean germplasms harboring the qSw17-1 QTL for the big-seeded phenotype indicated that reduced expression of GmKIX8-1 is the genetic basis of the qSw17-1 phenotype.

High expression in maize pollen correlates with genetic contributions to pollen fitness as well as with coordinated transcription from neighboring transposable elements

In flowering plants, gene expression in the haploid male gametophyte (pollen) is essential for sperm delivery and double fertilization. Pollen also undergoes dynamic epigenetic regulation of expression from transposable elements (TEs), but how this process interacts with gene expression is not clearly understood. To explore relationships among these processes, we quantified transcript levels in four male reproductive stages of maize (tassel primordia, microspores, mature pollen, and sperm cells) via RNA-seq. We found that, in contrast with vegetative cell-limited TE expression in Arabidopsis pollen, TE transcripts in maize accumulate as early as the microspore stage and are also present in sperm cells. Intriguingly, coordinate expression was observed between highly expressed protein-coding genes and their neighboring TEs, specifically in mature pollen and sperm cells. To investigate a potential relationship between elevated gene transcript level and pollen function, we measured the fitness cost (male-specific transmission defect) of GFP-tagged coding sequence insertion mutations in over 50 genes identified as highly expressed in the pollen vegetative cell, sperm cell, or seedling (as a sporophytic control). Insertions in seedling genes or sperm cell genes (with one exception) exhibited no difference from the expected 1:1 transmission ratio. In contrast, insertions in over 20% of vegetative cell genes were associated with significant reductions in fitness, showing a positive correlation of transcript level with non-Mendelian segregation when mutant. Insertions in maize gamete expressed2 (Zm gex2), the sole sperm cell gene with measured contributions to fitness, also triggered seed defects when crossed as a male, indicating a conserved role in double fertilization, given the similar phenotype previously demonstrated for the Arabidopsis ortholog GEX2. Overall, our study demonstrates a developmentally programmed and coordinated transcriptional activation of TEs and genes in pollen, and further identifies maize pollen as a model in which transcriptomic data have predictive value for quantitative phenotypes.

Novel eye genes systematically discovered through an integrated analysis of mouse transcriptomes and phenome

In the last few decades, reverse genetic and high throughput approaches have been frequently applied to the mouse (Mus musculus) to understand how genes function in tissues/organs and during development in a mammalian system. Despite these efforts, the associated phenotypes for the majority of mouse genes remained to be fully characterized. Here, we performed an integrated transcriptome-phenome analysis by identifying coexpressed gene modules based on tissue transcriptomes profiled with each of various platforms and functionally interpreting these modules using the mouse phenotypic data. Consequently, >15,000 mouse genes were linked with at least one of the 47 tissue functions that were examined. Specifically, our approach predicted >50 genes previously unknown to be involved in mice (Mus musculus) visual functions. Fifteen genes were selected for further analysis based on their potential biomedical relevance and compatibility with further experimental validation. Gene-specific morpholinos were introduced into zebrafish (Danio rerio) to target their corresponding orthologs. Quantitative assessments of phenotypes of developing eyes confirmed predicted eye-related functions of 13 out of the 15 genes examined. These novel eye genes include: Adal, Ankrd33, Car14, Ccdc126, Dhx32, Dkk3, Fam169a, Grifin, Kcnj14, Lrit2, Ppef2, Ppm1n, and Wdr17. The results highlighted the potential for this phenome-based approach to assist the experimental design of mutating and phenotyping mouse genes that aims to fully reveal the functional landscape of mammalian genomes.

Deficiency of AtGFAT1 activity impairs growth, pollen germination and tolerance to tunicamycin in Arabidopsis

The hexosamine biosynthetic pathway (HBP) plays essential roles in growth and development in plants. However, insight into the biological function of glutamine:fructose-6-phosphate amidotransferase 1 (GFAT1), mediating the first regulatory step of the HBP, remains unclear in plants. Here, we report the molecular characterization of Arabidopsis AtGFAT1 gene. AtGFAT1 was highly expressed in mature pollen grains, but its expression was not detectable in the rest of the organs. Pollen grains bearing the gfat1-2 knockout allele displayed defects in a polar deposition of pectin and callose in the pollen cell wall, leading to no genetic transmission of the gfat1-2 allele through the male gametophyte. AtGFAT1 overexpression increased glucosamine (GlcN) content and enhanced resistance to tunicamycin (Tm) treatment, while RNAi-mediated suppression reduced GlcN content and resistance to Tm treatment. However, the decrease in Tm resistance by RNAi suppression of AtGFAT1 was recovered by a GlcN supplement. The exogenous GlcN supplement also rescued gfat1-2/gaft1-2 mutant plants, which were otherwise not viable. The gfat1-2/gfat1-2 plants stopped growing at the germination stage on GlcN-free medium, but GlcN supplement allowed wild-type growth of gfat1-2/gfat1-2 plants. In addition, reactive oxygen species production, cell death and a decrease in protein N-glycosylation were observed in gfat1-2/gaft1-2 mutant plants grown on GlcN-free medium, whereas these aberrant defects were not detectable on GlcN-sufficient medium. Taken together, these results show that the reduction of protein N-glycosylation was at least partially responsible for many aberrant phenotypes in growth and development as well as the response to Tm treatment caused by AtGFAT1 deficiency in Arabidopsis.

Lineage-Specific Expression Divergence in Grasses Is Associated with Male Reproduction, Host-Pathogen Defense, and Domestication

Poaceae (grasses) is an agriculturally important and widely distributed family of plants with extraordinary phenotypic diversity, much of which was generated under recent lineage-specific evolution. Yet, little is known about the genes and functional modules involved in the lineage-specific divergence of grasses. Here, I address this question on a genome-wide scale by applying a novel branch-based statistic of lineage-specific expression divergence, LED, to RNA-seq data from nine tissues of the wild grass Brachypodium distachyon and its domesticated relatives Oryza sativa japonica (rice) and Sorghum bicolor (sorghum). I find that LED is generally smallest in B. distachyon and largest in O. sativa japonica, which underwent domestication earlier than S. bicolor, supporting the hypothesis that domestication may increase the rate of lineage-specific expression divergence in grasses. Moreover, in all three species, LED is positively correlated with protein-coding sequence divergence and tissue specificity, and negatively correlated with network connectivity. Further analysis reveals that genes with large LED are often primarily expressed in anther, implicating lineage-specific expression divergence in the evolution of male reproductive phenotypes. Gene ontology enrichment analysis also identifies an overrepresentation of terms related to male reproduction in the two domesticated grasses, as well as to those involved in host-pathogen defense in all three species. Last, examinations of genes with the largest LED reveal that their lineage-specific expression divergence may have contributed to antimicrobial functions in B. distachyon, to enhanced adaptation and yield during domestication in O. sativa japonica, and to defense against a widespread and devastating fungal pathogen in S. bicolor. Together, these findings suggest that lineage-specific expression divergence in grasses may increase under domestication and preferentially target rapidly evolving genes involved in male reproduction, host-pathogen defense, and the origin of domesticated phenotypes.

Recruitment of histone modifications to assist mRNA dosage maintenance after degeneration of cytosine DNA methylation during animal evolution

Following gene duplication, mRNA expression of the duplicated gene is reduced to maintain mRNA dosage. In mammals, this process is achieved with increased cytosine DNA methylation of the promoters of duplicated genes to suppress transcriptional initiation. However, not all animal species possess a full apparatus for cytosine DNA methylation. For such species, such as the roundworm (Caenorhabditis elegans, "worm" hereafter) or fruit fly (Drosophila melanogaster, "fly" hereafter), it is unclear how reduced expression of duplicated genes has been achieved evolutionarily. Here, we hypothesize that in the absence of a classical cytosine DNA methylation pathway, histone modifications play an increasing role in maintaining mRNA dosage following gene duplication. We initially verified that reduced gene expression of duplicated genes had occurred in the worm, fly, and mouse (Mus musculus). Next, several histone marks, with the capacity to control mRNA abundance in the models studied, were examined. In the worm and fly, but not in the mouse, multiple histone modifications were found to assist mRNA dosage maintenance following gene duplication events and the possible involvement of adenine DNA methylation in this process was excluded. Furthermore, the histone marks and acting regions that mediated the reduction in duplicated gene expression were found to be largely organism specific. Thus, it appears that many of the histone marks that maintain mRNA dosage were independently recruited during the evolution of worms and flies to compensate for the loss of cytosine DNA methylation machinery from their genomes.

Integrated Transcriptome Map Highlights Structural and Functional Aspects of the Normal Human Heart

A systematic meta-analysis of the available gene expression profiling datasets for the whole normal human heart generated a quantitative transcriptome reference map of this organ. Transcriptome Mapper (TRAM) software integrated 32 gene expression profile datasets from different sources returning a reference value of expression for each of the 43,360 known, mapped transcripts assayed by any of the experimental platforms used in this regard. Main findings include the visualization at the gene and chromosomal levels of the classical description of the basic histology and physiology of the heart, the identification of suitable housekeeping reference genes, the analysis of stoichiometry of gene products, and the focusing on chromosome 21 genes, which are present in one excess copy in Down syndrome subjects, presenting cardiovascular defects in 30-40% of cases. Independent in vitro validation showed an excellent correlation coefficient (r = 0.98) with the in silico data. Remarkably, heart/non-cardiac tissue expression ratio may also be used to anticipate that effects of mutations will most probably affect or not the heart. The quantitative reference global portrait of gene expression in the whole normal human heart illustrates the structural and functional aspects of the whole organ and is a general model to understand the mechanisms underlying heart pathophysiology. J. Cell. Physiol. 232: 759-770, 2017. © 2016 Wiley Periodicals, Inc.

Development and validation of a custom microarray for global transcriptome profiling of the fungus Aspergillus nidulans

Transcriptome profiling is a powerful tool for identifying gene networks from whole genome expression analysis in many living species. Here is described the first extensively characterized platform using Agilent microarray technology for transcriptome analysis in the filamentous fungus Aspergillus (Emericella) nidulans. We developed and validated a reliable gene expression microarray in 8 × 15 K format, with predictive and experimental data establishing its specificity and sensitivity. Either one or two 60-mer oligonucleotide probes were selected for each of 10,550 nuclear as well as 20 mitochondrial coding sequences. More than 99 % of probes were predicted to hybridize with 100 % identity to their aimed specific A. nidulans target only. Probe sensitivity was supported by a highly narrow distribution of melting temperatures together with thermodynamic features, which strongly favored probe-target perfect match hybridization, in comparison with predicted secondary structures. Array quality was evaluated through transcriptome comparison of two A. nidulans strains, differing by the presence or not of Escherichia coli LacZ transgene. High signal-to-noise ratios were measured, and signal reproducibility was established at intra-probe and inter-probe levels. Reproducibility of microarray performances was assessed by high correlation between two-color dye signals and between technical replicates. Results were confirmed by RT-qPCR analysis on five genes. Though it covers 100 % of the A. nidulans targeted coding sequences, this low density array allows limited experimental costs and simplified data analysis process, making it suitable for studying gene expression in this model organism through large numbers of experimental conditions, in basic, biomedical or industrial microbiology research fields.

Functionalities of expressed messenger RNAs revealed from mutant phenotypes

Total messenger RNAs mRNAs that are produced from a given gene under a certain set of conditions include both functional and nonfunctional transcripts. The high prevalence of nonfunctional mRNAs that have been detected in cells has raised questions regarding the functional implications of mRNA expression patterns and divergences. Phenotypes that result from the mutagenesis of protein-coding genes have provided the most straightforward descriptions of gene functions, and such data obtained from model organisms have facilitated investigations of the functionalities of expressed mRNAs. Mutant phenotype data from mouse tissues have revealed various attributes of functional mRNAs, including tissue-specificity, strength of expression, and evolutionary conservation. In addition, the role that mRNA expression evolution plays in driving morphological evolution has been revealed from studies designed to exploit morphological and physiological phenotypes of mouse mutants. Investigations into yeast essential genes (defined by an absence of colony growth after gene deletion) have further described gene regulatory strategies that reduce protein expression noise by mediating the rates of transcription and translation. In addition to the functional significance of expressed mRNAs as described in the abovementioned findings, the functionalities of other type of RNAs (i.e., noncoding RNAs) remain to be characterized with systematic mutations and phenotyping of the DNA regions that encode these RNA molecules. WIREs RNA 2016, 7:416-427. doi: 10.1002/wrna.1329 For further resources related to this article, please visit the WIREs website.