Widespread stable noncanonical peptides identified by integrated analyses of ribosome profiling and ORF features.

Studies have revealed dozens of functional peptides in putative 'noncoding' regions and raised the question of how many proteins are encoded by noncanonical open reading frames (ORFs). Here, we comprehensively annotate genome-wide translated ORFs across five eukaryotes (human, mouse, zebrafish, worm, and yeast) by analyzing ribosome profiling data. We develop a logistic regression model named PepScore based on ORF features (expected length, encoded domain, and conservation) to calculate the probability that the encoded peptide is stable in humans. Systematic ectopic expression validates PepScore and shows that stable complex-associating microproteins can be encoded in 5'/3' untranslated regions and overlapping coding regions of mRNAs besides annotated noncoding RNAs. Stable noncanonical proteins follow conventional rules and localize to different subcellular compartments. Inhibition of proteasomal/lysosomal degradation pathways can stabilize some peptides especially those with moderate PepScores, but cannot rescue the expression of short ones with low PepScores suggesting they are directly degraded by cellular proteases. The majority of human noncanonical peptides with high PepScores show longer lengths but low conservation across species/mammals, and hundreds contain trait-associated genetic variants. Our study presents a statistical framework to identify stable noncanonical peptides in the genome and provides a valuable resource for functional characterization of noncanonical translation during development and disease.

Low-input RNase footprinting for simultaneous quantification of cytosolic and mitochondrial translation.

We describe a low-input RNase footprinting approach for the rapid quantification of ribosome-protected fragments with as few as 1000 cultured cells. The assay uses a simplified procedure to selectively capture ribosome footprints based on optimized RNase digestion. It simultaneously maps cytosolic and mitochondrial translation with single-nucleotide resolution. We applied it to reveal selective functions of the elongation factor TUFM in mitochondrial translation, as well as synchronized repression of cytosolic translation after TUFM perturbation. We show the assay is applicable to small amounts of primary tissue samples with low protein synthesis rates, including snap-frozen tissues and immune cells from an individual's blood draw. We showed its feasibility to characterize the personalized immuno-translatome. Our analyses revealed that thousands of genes show lower translation efficiency in monocytes compared with lymphocytes, and identified thousands of translated noncanonical open reading frames (ORFs). Altogether, our RNase footprinting approach opens an avenue to assay transcriptome-wide translation using low-input samples from a wide range of physiological conditions.

Dynamic sex chromosome expression in Drosophila male germ cells.

Given their copy number differences and unique modes of inheritance, the evolved gene content and expression of sex chromosomes is unusual. In many organisms the X and Y chromosomes are inactivated in spermatocytes, possibly as a defense mechanism against insertions into unpaired chromatin. In addition to current sex chromosomes, Drosophila has a small gene-poor X-chromosome relic (4th) that re-acquired autosomal status. Here we use single cell RNA-Seq on fly larvae to demonstrate that the single X and pair of 4th chromosomes are specifically inactivated in primary spermatocytes, based on measuring all genes or a set of broadly expressed genes in testis we identified. In contrast, genes on the single Y chromosome become maximally active in primary spermatocytes. Reduced X transcript levels are due to failed activation of RNA-Polymerase-II by phosphorylation of Serine 2 and 5.

Drosophila Heterochromatin Stabilization Requires the Zinc-Finger Protein Small Ovary.

Heterochromatin-mediated repression is essential for controlling the expression of transposons and for coordinated cell type-specific gene regulation. The small ovary (sov) locus was identified in a screen for female-sterile mutations in Drosophila melanogaster, and mutants show dramatic ovarian morphogenesis defects. We show that the null sov phenotype is lethal and map the locus to the uncharacterized gene CG14438, which encodes a nuclear zinc-finger protein that colocalizes with the essential Heterochromatin Protein 1 (HP1a). We demonstrate Sov functions to repress inappropriate gene expression in the ovary, silence transposons, and suppress position-effect variegation in the eye, suggesting a central role in heterochromatin stabilization.

Re-annotation of eight Drosophila genomes.

The sequenced genomes of the Drosophila phylogeny are a central resource for comparative work supporting the understanding of the Drosophila melanogaster non-mammalian model system. These have also facilitated evolutionary studies on the selected and random differences that distinguish the thousands of extant species of Drosophila. However, full utility has been hampered by uneven genome annotation. We have generated a large expression profile dataset for nine species of Drosophila and trained a transcriptome assembly approach on D. melanogaster that best matched the extensively curated annotation. We then applied this to the other species to add more than 10000 transcript models per species. We also developed new orthologs to facilitate cross-species comparisons. We validated the new annotation of the distantly related Drosophila grimshawi with an extensive collection of newly sequenced cDNAs. This re-annotation will facilitate understanding both the core commonalities and the species differences in this important group of model organisms, and suggests a strategy for annotating the many forthcoming genomes covering the tree of life.

Drosophila melanogaster positive transcriptional elongation factors regulate metabolic and sex-biased expression in adults.

BACKGROUND: Transcriptional elongation is a generic function, but is also regulated to allow rapid transcription responses. Following relatively long initiation and promoter clearance, RNA polymerase II can pause and then rapidly elongate following recruitment of positive elongation factors. Multiple elongation complexes exist, but the role of specific components in adult Drosophila is underexplored. RESULTS: We conducted RNA-seq experiments to analyze the effect of RNAi knockdown of Suppressor of Triplolethal and lilliputian. We similarly analyzed the effect of expressing a dominant negative Cyclin-dependent kinase 9 allele. We observed that almost half of the genes expressed in adults showed reduced expression, supporting a broad role for the three tested genes in steady-state transcript abundance. Expression profiles following lilliputian and Suppressor of Triplolethal RNAi were nearly identical raising the possibility that they are obligatory co-factors. Genes showing reduced expression due to these RNAi treatments were short and enriched for genes encoding metabolic or enzymatic functions. The dominant-negative Cyclin-dependent kinase 9 profiles showed both overlapping and specific differential expression, suggesting involvement in multiple complexes. We also observed hundreds of genes with sex-biased differential expression following treatment. CONCLUSION: Transcriptional profiles suggest that Lilliputian and Suppressor of Triplolethal are obligatory cofactors in the adult and that they can also function with Cyclin-dependent kinase 9 at a subset of loci. Our results suggest that transcriptional elongation control is especially important for rapidly expressed genes to support digestion and metabolism, many of which have sex-biased function.

Sxl-Dependent, tra/tra2-Independent Alternative Splicing of the Drosophila melanogaster X-Linked Gene found in neurons.

Somatic sexual determination and behavior in Drosophila melanogaster are under the control of a genetic cascade initiated by Sex lethal (Sxl). In the female soma, SXL RNA-binding protein regulates the splicing of transformer (tra) transcripts into a female-specific form. The RNA-binding protein TRA and its cofactor TRA2 function in concert in females, whereas SXL, TRA, and TRA2 are thought to not function in males. To better understand sex-specific regulation of gene expression, we analyzed male and female head transcriptome datasets for expression levels and splicing, quantifying sex-biased gene expression via RNA-Seq and qPCR. Our data uncouple the effects of Sxl and tra/tra2 in females in the-sex-biased alternative splicing of head transcripts from the X-linked locus found in neurons (fne), encoding a pan-neuronal RNA-binding protein of the ELAV family. We show that FNE protein levels are downregulated by Sxl in female heads, also independently of tra/tra2. We argue that this regulation may have important sexually dimorphic consequences for the regulation of nervous system development or function.

Expression profile and gene age jointly shaped the genome-wide distribution of premature termination codons in a Drosophila melanogaster population.

Widespread premature termination codon mutations (PTCs) were recently observed in human and fly populations. We took advantage of the population resequencing data in the Drosophila Genetic Reference Panel to investigate how the expression profile and the evolutionary age of genes shaped the allele frequency distribution of PTCs. After generating a high-quality data set of PTCs, we clustered genes harboring PTCs into three categories: genes encoding low-frequency PTCs (≤ 1.5%), moderate-frequency PTCs (1.5-10%), and high-frequency PTCs (>10%). All three groups show narrow transcription compared with PTC-free genes, with the moderate- and high-PTC frequency groups showing a pronounced pattern. Moreover, nearly half (42%) of the PTC-encoding genes are not expressed in any tissue. Interestingly, the moderate-frequency PTC group is strongly enriched for genes expressed in midgut, whereas genes harboring high-frequency PTCs tend to have sex-specific expression. We further find that although young genes born in the last 60 My compose a mere 9% of the genome, they represent 16%, 30%, and 50% of the genes containing low-, moderate-, and high-frequency PTCs, respectively. Among DNA-based and RNA-based duplicated genes, the child copy is approximately twice as likely to contain PTCs as the parent copy, whereas young de novo genes are as likely to encode PTCs as DNA-based duplicated new genes. Based on these results, we conclude that expression profile and gene age jointly shaped the landscape of PTC-mediated gene loss. Therefore, we propose that new genes may need a long time to become stably maintained after the origination.

Roles of young serine-endopeptidase genes in survival and reproduction revealed rapid evolution of phenotypic effects at adult stages.

Our recent study found that 30% of young genes were essential for viability that determines development through stages from embryo to pupae in Drosophila melanogaster, revealing rapidly evolving genetic components involved in the evolution of development. Meanwhile, many young genes did not produce complete lethal phenotype upon constitutive knockdown, suggesting that they may not be essential for viability. These genes, nevertheless, were fixed by natural selection, and might play an important functional role in their adult stage. Here we present a detailed demonstration that a newly duplicated serine-type endopeptidase gene that originated in the common ancestor in the D. melanogaster subgroup 6~11 million years ago, named Slfc, revealing a strong effect in post-eclosion. Although animals survived constitutive knockdown of Slfc to adult stage, however, their life span reduced significantly by two-thirds compared to wildtype. Furthermore, the Slfc-RNAi males dropped their fertility to less than 10% of the wildtype level, with over 80% of these males being sterile. The Slfc-RNAi females, on the other hand, showed a slight reduction in fertility. This case study demonstrates that a young gene can contribute to fitness on the three important traits of life history in adults, including the life expectancy, male fertility and female fertility, suggesting that new genes can quickly evolve and impact multiple phenotypes.

Increased complexity of gene structure and base composition in vertebrates.

How the structure and base composition of genes changed with the evolution of vertebrates remains a puzzling question. Here we analyzed 895 orthologous protein-coding genes in six multicellular animals: human, chicken, zebrafish, sea squirt, fruit fly, and worm. Our analyses reveal that many gene regions, particularly intron and 3' UTR, gradually expanded throughout the evolution of vertebrates from their invertebrate ancestors, and that the number of exons per gene increased. Studies based on all protein-coding genes in each genome provide consistent results. We also find that GC-content increased in many gene regions (especially 5' UTR) in the evolution of endotherms, except in coding-exons. Analysis of individual genomes shows that 3' UTR demonstrated stronger length and GC-content correlation with intron than 5' UTR, and gene with large intron in all six species demonstrated relatively similar GC-content. Our data indicates a great increase in complexity in vertebrate genes and we propose that the requirement for morphological and functional changes is probably the driving force behind the evolution of structure and base composition complexity in multicellular animal genes.

Important role of indels in somatic mutations of human cancer genes.

BACKGROUND: Cancer is clonal proliferation that arises owing to mutations in a subset of genes that confer growth advantage. More and more cancer related genes are found to have accumulated somatic mutations. However, little has been reported about mutational patterns of insertions/deletions (indels) in these genes. RESULTS: We analyzed indels' abundance and distribution, the relative ratio between indels and somatic base substitutions and the association between those two forms of mutations in a large number of somatic mutations in the Catalogue of Somatic Mutations in Cancer database. We found a strong correlation between indels and base substitutions in cancer-related genes and showed that they tend to concentrate at the same locus in the coding sequences within the same samples. More importantly, a much higher proportion of indels were observed in somatic mutations, as compared to meiotic ones. Furthermore, our analysis demonstrated a great diversity of indels at some loci of cancer-related genes. Particularly in the genes with abundant mutations, the proportion of 3n indels in oncogenes is 7.9 times higher than that in tumor suppressor genes. CONCLUSIONS: There are three distinct patterns of indel distribution in somatic mutations: high proportion, great abundance and non-random distribution. Because of the great influence of indels on gene function (e.g., the effect of frameshift mutation), these patterns indicate that indels are frequently under positive selection and can often be the 'driver mutations' in oncogenesis. Such driver forces can better explain why much less frameshift mutations are in oncogenes while much more in tumor suppressor genes, because of their different function in oncogenesis. These findings contribute to our understanding of mutational patterns and the relationship between indels and cancer.

Variation in the ratio of nucleotide substitution and indel rates across genomes in mammals and bacteria.

Rates of nucleotide substitution and insertion/deletion (indel) are known to vary across the functional components of a genome. Little attention has been paid, however, to the quantitative relationship between the two. Here we investigate the ratio of nucleotide substitutions to indels (S/I) in different regions of 4 primates, 70 bacteria, and 8 other genomes. We find that the ratio differs at 5.4-times between coding and noncoding, 3.3-times between conserved and less conserved coding sequences, and 1.46-times between nonrepeat and repeat regions. The S/I ratio is also positively correlated with the level of divergence between the genomes compared. Our results suggest that the S/I ratio may reflect differences in the efficacy of selection against indels. Due to the sensitivity of indel density in different regions, this ratio varies over a much larger range. With the recent discovery suggesting that indels act as local enhancers of mutation in surrounding sequences, nucleotide substitution rates are expected to be accelerated in regions of low constraint, where indels tend to accumulate, but will otherwise be modulated in proportion to the level of a sequence's functional constraint. Indels, therefore, may play a nontrivial role in controlling differences in genetic variation and divergence across functional regions of a genome.

Evolutionary pattern of protein architecture in mammal and fruit fly genomes.

Mutations, which can alter amino acid constitution, contribute greatly to protein evolution. However, little is reported of their pattern during protein structural evolution. We investigated the distribution of non-synonymous single nucleotide polymorphisms (nsSNPs) and insertions/deletions (indels) along mammal and fruit fly proteins. We found the nsSNPs (and d(N)) and indels increased in protein boundary regions, and this pattern is inversely correlated with the distribution of protein domain density. Additionally, synonymous substitutions (and d(S)) are reduced in 5' and 3' regions, indicating more variable protein boundaries, compared with central interior. All evidence suggests that the inner part of coding sequences (CDSs) is comparatively conserved, whereas the 5' and 3' regions, with higher evolution rates, are more variable. We assumed that due to greater frequencies of nsSNPs and indels in adaptive regions of CDSs it could be easier to ultimately alter, gain, or lose amino acids, thus becoming the front line of protein evolution.