A Novel, Improved, Application for the Normalization of RNA-seq Based Expression Data in Complex Polyploids.

Much of the work on the normalization of RNA-seq data has been performed on human, notably cancer tissue. Little work has been done in plants, particularly polyploids and those species with incomplete or no genomes. We present a novel implementation of GeTMM (Gene Length Corrected TMM) that accounts for GC bias and works at the transcript level. The algorithm also employs transcript length as a factor, allowing for incomplete transcripts and alternate transcripts. This significantly improves overall normalization. The GCGeTMM methodology also allows for simultaneous determination of differentially expressed transcripts (and by extension genes) and stably expressed genes to act as references for qRT-PCR and microarray analyses.

Metastatic and recurrent adrenocortical cancer is not defined by its genomic landscape

Background Adrenocortical carcinoma (ACC) is a rare, often-aggressive neoplasm of the adrenal cortex, with a 14–17 month median overall survival. We asked whether tumors from patients with advanced or metastatic ACC would offer clues as to putative genes that might have critical roles in disease progression or in more aggressive disease biology. Methods We conducted comprehensive genomic and expression analyses of ACCs from 43 patients, 30 female, and 42 from metastatic sites, including deep sequencing, copy number analysis, mRNA expression and microRNA arrays. Results Copy number gains and losses were similar to that previously reported for ACC. We identified a median mutation rate of 3.38 per megabase (Mb). The mutational signature was characterized by a predominance of C > T, C > A and T > C transitions. Only cancer genes TP53 (26%) and beta-catenin (CTNNB1, 14%) were mutated in more than 10% of samples. The TCGA-identified putative cancer genes MEN1 and PRKAR1A were found in low frequency—4.7 and 2.3%, respectively. The majority of the mutations were in genes not implicated in the etiology or maintenance of cancer. Specifically, amongst the 38 genes that were mutated in more than 9% of samples, only four were represented in Tier 1 of the 576 COSMIC Cancer Gene Census (CCGC). Thus, 82% of genes found to have mutations likely have no role in the etiology or biology of ACC; while the role of the other 18%, if any, remains to be proven. Finally, the transcript length for the 38 most frequently mutated genes in ACC is statistically longer than the average of all coding genes, raising the question of whether transcript length in part determined mutation probability. Conclusions We conclude that the mutational and expression profiles of advanced and metastatic tumors are very similar to those from newly diagnosed patients—with very little in the way of genomic aberration to explain differences in biology. With relatively low mutation rates, few major oncogenic drivers, and loss of function mutations in several epigenetic regulators, an epigenetic basis for ACC may be postulated and serve as the basis for future studies.

Single-Cell RNA Sequencing Reveals mRNA Splice Isoform Switching during Kidney Development

BackgroundDuring mammalian kidney development, nephron progenitors undergo a mesenchymal-to-epithelial transition and eventually differentiate into the various tubular segments of the nephron. Recently, Drop-seq single-cell RNA sequencing technology for measuring gene expression from thousands of individual cells identified the different cell types in the developing kidney. However, that analysis did not include the additional layer of heterogeneity that alternative mRNA splicing creates.MethodsFull transcript length single-cell RNA sequencing characterized the transcriptomes of 544 individual cells from mouse embryonic kidneys.ResultsGene expression levels measured with full transcript length single-cell RNA sequencing identified each cell type. Further analysis comprehensively characterized splice isoform switching during the transition between mesenchymal and epithelial cellular states, which is a key transitional process in kidney development. The study also identified several putative splicing regulators, including the genes Esrp1/2 and Rbfox1/2.ConclusionsDiscovery of the sets of genes that are alternatively spliced as the fetal kidney mesenchyme differentiates into tubular epithelium will improve our understanding of the molecular mechanisms that drive kidney development.

Transgenic Testing Does Not Support a Role for Additional Candidate Genes in Wolbachia Male Killing or Cytoplasmic Incompatibility

ABSTRACT Endosymbiotic bacteria in the genus Wolbachia remarkably infect nearly half of all arthropod species. They spread in part because of manipulations of host sexual reproduction that enhance the maternal transmission of the bacteria, including male killing (death of infected males) and unidirectional cytoplasmic incompatibility (CI; death of offspring from infected fathers and uninfected mothers). Recent discoveries identified several genes in prophage WO of Wolbachia (wmk, cifA, and cifB) that fully or partially recapitulate male killing or CI when transgenically expressed in Drosophila melanogaster. However, it is not yet fully resolved if other gene candidates contribute to these phenotypes. Here, we transgenically tested 10 additional gene candidates for their involvement in male killing and/or CI. The results show that despite sequence and protein architecture similarities or comparative associations with reproductive parasitism, transgenic expression of the candidates does not recapitulate male killing or CI. Sequence analysis across Wmk and its closest relatives reveals amino acids that may be important to its function. In addition, evidence is presented to propose new hypotheses regarding the relationship between wmk transcript length and its ability to kill a given host, as well as copy number of wmk homologs within a bacterial strain, which may be predictive of host resistance. Together, these analyses continue to build the evidence for identification of wmk, cifA, and cifB as the major genes that have thus far been shown to cause reproductive parasitism in Wolbachia, and the transgenic resources provide a basis for further functional study of phage WO genes. IMPORTANCE Wolbachia are widespread bacterial endosymbionts that manipulate the reproduction of diverse arthropods to spread through a population and can substantially shape host evolution. Recently, reports identified three prophage WO genes (wmk, cifA, and cifB) that transgenically recapitulate many aspects of reproductive manipulation in Drosophila melanogaster. Here, we transgenically tested 10 additional gene candidates for CI and/or male killing in flies. The results yield no evidence for the involvement of these gene candidates in reproductive parasitism, bolstering the evidence for identification of the cif and wmk genes as the major factors involved in their phenotypes. In addition, evidence supports new hypotheses for prediction of male-killing phenotypes or lack thereof based on wmk transcript length and copy number. These experiments inform efforts to understand the full basis of reproductive parasitism for basic and applied purposes and lay the foundation for future work on the function of an interesting group of Wolbachia and phage WO genes.

Gene length as a regulator for ribosome recruitment and protein synthesis: theoretical insights

Protein synthesis rates are determined, at the translational level, by properties of the transcript’s sequence. The efficiency of an mRNA can be tuned by varying the ribosome binding sites controlling the recruitment of the ribosomes, or the codon usage establishing the speed of protein elongation. In this work we propose transcript length as a further key determinant of translation efficiency. Based on a physical model that considers the kinetics of ribosomes advancing on the mRNA and diffusing in its surrounding, as well as mRNA circularisation and ribosome drop-off, we explain how the transcript length may play a central role in establishing ribosome recruitment and the overall translation rate of an mRNA. We also demonstrate how this process may be involved in shaping the experimental ribosome density-gene length dependence. Finally, we argue that cells could exploit this mechanism to adjust and balance the usage of its ribosomal resources.

Ribosome reinitiation can explain length-dependent translation of messenger RNA

Models of mRNA translation usually presume that transcripts are linear; upon reaching the end of a transcript each terminating ribosome returns to the cytoplasmic pool before initiating anew on a different transcript. A consequence of linear models is that faster translation of a given mRNA is unlikely to generate more of the encoded protein, particularly at low ribosome availability. Recent evidence indicates that eukaryotic mRNAs are circularized, potentially allowing terminating ribosomes to preferentially reinitiate on the same transcript. Here we model the effect of ribosome reinitiation on translation and show that, at high levels of reinitiation, protein synthesis rates are dominated by the time required to translate a given transcript. Our model provides a simple mechanistic explanation for many previously enigmatic features of eukaryotic translation, including the negative correlation of both ribosome densities and protein abundance on transcript length, the importance of codon usage in determining protein synthesis rates, and the negative correlation between transcript length and both codon adaptation and 5' mRNA folding energies. In contrast to linear models where translation is largely limited by initiation rates, our model reveals that all three stages of translation - initiation, elongation, and termination/reinitiation - determine protein synthesis rates even at low ribosome availability.