A mutation creating an upstream initiation codon in the SOX9 5' UTR causes acampomelic campomelic dysplasia

Variants in the SOX9 transactivation middle domain induce axial skeleton dysplasia and scoliosis

SOX9 is an essential transcriptional regulator of cartilage development and homeostasis. In humans, dysregulation of SOX9 is associated with a wide spectrum of skeletal disorders, including campomelic and acampomelic dysplasia, and scoliosis. The mechanism of how SOX9 variants contribute to the spectrum of axial skeletal disorders is not well understood. Here, we report four novel pathogenic variants of SOX9 identified in a large cohort of patients with congenital vertebral malformations. Three of these heterozygous variants are in the HMG and DIM domains, and for the first time, we report a pathogenic variant within the transactivation middle (TAM) domain of SOX9 . Probands with these variants exhibit variable skeletal dysplasia, ranging from isolated vertebral malformation to acampomelic dysplasia. We also generated a Sox9 hypomorphic mutant mouse model bearing a microdeletion within the TAM domain ( Sox9 ^Asp272del ). We demonstrated that disturbance of the TAM domain with missense mutation or microdeletion results in reduced protein stability but does not affect the transcriptional activity of SOX9. Homozygous Sox9 ^Asp272del mice exhibited axial skeletal dysplasia including kinked tails, ribcage anomalies, and scoliosis, recapitulating phenotypes observed in human, while heterozygous mutants display a milder phenotype. Analysis of primary chondrocytes and the intervertebral discs in Sox9 ^Asp272del mutant mice revealed dysregulation of a panel of genes with major contributions of the extracellular matrix, angiogenesis, and ossification-related processes. In summary, our work identified the first pathologic variant of SOX9 within the TAM domain and demonstrated that this variant is associated with reduced SOX9 protein stability. Our finding suggests that reduced SOX9 stability caused by variants in the TAM domain may be responsible for the milder forms of axial skeleton dysplasia in humans.

Novel start codons introduce novel coding sequences in the human genomes

Start-gain mutations can introduce novel start codons and generate novel coding sequences that may affect the function of genes. In this study, we systematically investigated the novel start codons that were either polymorphic or fixed in the human genomes. 829 polymorphic start-gain SNVs were identified in the human populations, and the novel start codons introduced by these SNVs have significantly higher activity in translation initiation. Some of these start-gain SNVs were reported to be associated with phenotypes and diseases in previous studies. By comparative genomic analysis, we found 26 human-specific start codons that were fixed after the divergence between the human and chimpanzee, and high-level translation initiation activity was observed on them. The negative selection signal was detected in the novel coding sequences introduced by these human-specific start codons, indicating the important function of these novel coding sequences.

Annotation of uORFs in the OMIM genes allows to reveal pathogenic variants in 5'UTRs

An increasing number of studies emphasize the role of non-coding variants in the development of hereditary diseases. However, the interpretation of such variants in clinical genetic testing still remains a critical challenge due to poor knowledge of their pathogenicity mechanisms. It was previously shown that variants in 5'-untranslated regions (5'UTRs) can lead to hereditary diseases due to disruption of upstream open reading frames (uORFs). Here, we performed a manual annotation of upstream translation initiation sites (TISs) in human disease-associated genes from the OMIM database and revealed ∼4.7 thousand of TISs related to uORFs. We compared our TISs with the previous studies and provided a list of 'high confidence' uORFs. Using a luciferase assay, we experimentally validated the translation of uORFs in the ETFDH, PAX9, MAST1, HTT, TTN,GLI2 and COL2A1 genes, as well as existence of N-terminal CDS extension in the ZIC2 gene. Besides, we created a tool to annotate the effects of genetic variants located in uORFs. We revealed the variants from the HGMD and ClinVar databases that disrupt uORFs and thereby could lead to Mendelian disorders. We also showed that the distribution of uORFs-affecting variants differs between pathogenic and population variants. Finally, drawing on manually curated data, we developed a machine-learning algorithm that allows us to predict the TISs in other human genes.

Hypomorphic and dominant-negative impact of truncated SOX9 dysregulates Hedgehog-Wnt signaling, causing campomelia

Haploinsufficiency for SOX9, the master chondrogenesis transcription factor, can underlie campomelic dysplasia (CD), an autosomal dominant skeletal malformation syndrome, because heterozygous Sox9 null mice recapitulate the bent limb (campomelia) and some other phenotypes associated with CD. However, in vitro cell assays suggest haploinsufficiency may not apply for certain mutations, notably those that truncate the protein, but in these cases in vivo evidence is lacking and underlying mechanisms are unknown. Here, using conditional mouse mutants, we compared the impact of a heterozygous Sox9 null mutation (Sox9+/-) with the Sox9+/Y440X CD mutation that truncates the C-terminal transactivation domain but spares the DNA-binding domain. While some Sox9+/Y440X mice survived, all Sox9+/- mice died perinatally. However, the skeletal defects were more severe and IHH signaling in developing limb cartilage was significantly enhanced in Sox9+/Y440X compared with Sox9+/-. Activating Sox9Y440X specifically in the chondrocyte-osteoblast lineage caused milder campomelia, and revealed cell- and noncell autonomous mechanisms acting on chondrocyte differentiation and osteogenesis in the perichondrium. Transcriptome analyses of developing Sox9+/Y440X limbs revealed dysregulated expression of genes for the extracellular matrix, as well as changes consistent with aberrant WNT and HH signaling. SOX9Y440X failed to interact with β-catenin and was unable to suppress transactivation of Ihh in cell-based assays. We propose enhanced HH signaling in the adjacent perichondrium induces asymmetrically localized excessive perichondrial osteogenesis resulting in campomelia. Our study implicates combined haploinsufficiency/hypomorphic and dominant-negative actions of SOX9Y440X, cell-autonomous and noncell autonomous mechanisms, and dysregulated WNT and HH signaling, as the cause of human campomelia.

Autosomal Recessive ACTG2-Related Visceral Myopathy in Brothers

Pediatric intestinal pseudo-obstruction (PIPO) is a heterogeneous condition characterized by impaired gastrointestinal propulsion, a broad clinical spectrum, and variable severity. Several molecular bases underlying primary PIPO have been identified, of which autosomal dominant ACTG2-related visceral myopathy is the most common in both familial or sporadic primary PIPO cases. We present a family with autosomal recessive ACTG2-related disease in which both parents have mild gastrointestinal symptoms and sons have severe PIPO and bladder dysfunction.

Methods

Clinical genome sequencing was performed on the patients and the mother. Immunohistochemistry was performed on intestinal tissue from the patients to show expression levels of the ACTG2.

Results

Genome sequencing identified a 6.8 kb 2p13.1 loss that includes the ACTG2 gene and a maternally inherited missense variant p.Val10Met in the ACTG2 gene.

Discussion

This case demonstrates that monoallelic hypomorphic ACTG2 variants may underly mild primary gastrointestinal symptoms, while biallelic mild variants can cause severe diseases. The Deletions of the noncoding ACTG2 exon can be an under-recognized cause of mild gastrointestinal symptoms unidentifiable by exome sequencing, explaining some instances of interfamilial variability with an apparent autosomal dominant inheritance. Genome sequencing is recommended as a genetic work-up for primary or idiopathic PIPO because of genetic heterogeneity.

Exploring the Peptide Potential of Genomes

Recent studies attribute a central role to the noncoding genome in the emergence of novel genes. The widespread transcription of noncoding regions and the pervasive translation of the resulting RNAs offer to the organisms a vast reservoir of novel peptides. Although the majority of these peptides are anticipated as deleterious or neutral, and thereby expected to be degraded right away or short-lived in evolutionary history, some of them can confer an advantage to the organism. The latter can be further subjected to natural selection and be established as novel genes. In any case, characterizing the structural properties of these pervasively translated peptides is crucial to understand (1) their impact on the cell and (2) how some of these peptides, derived from presumed noncoding regions, can give rise to structured and functional de novo proteins. Therefore, we present a protocol that aims to explore the potential of a genome to produce novel peptides. It consists in annotating all the open reading frames (ORFs) of a genome (i.e., coding and noncoding ones) and characterizing the fold potential and other structural properties of their corresponding potential peptides. Here, we apply our protocol to a small genome and show how to apply it to very large genomes. Finally, we present a case study which aims to probe the fold potential of a set of 721 translated ORFs in mouse lncRNAs, identified with ribosome profiling experiments. Interestingly, we show that the distribution of their fold potential is different from that of the nontranslated lncRNAs and more generally from the other noncoding ORFs of the mouse.

Identification of deleterious single nucleotide polymorphism (SNP)s in the human TBX5 gene & prediction of their structural & functional consequences: An in silico approach

T-box transcription factor 5 gene (TBX5) encodes the transcription factor TBX5, which plays a crucial role in the development of heart and upper limbs. Damaging single nucleotide variants in this gene alter the protein structure, disturb the functions of TBX5, and ultimately cause Holt-Oram Syndrome (HOS). By analyzing the available single nucleotide polymorphism information in the dbSNP database, this study was designed to identify the most deleterious TBX5 SNPs through insilico approaches and predict their structural and functional consequences.

Fifty-eight missense substitutions were found damaging by sequence homology-based tools: SIFT and PROVEAN, and structure homology-based tool PolyPhen-2. Various disease association meta-predictors further scrutinized these SNPs. Additionally, conservation profile of the amino acid residues, their surface accessibility, stability, and structural integrity of the native protein upon mutations were assessed. From these analyses, finally 5 SNPs were detected as the most damaging ones: [rs1565941579 (P85S), rs1269970792 (W121R), rs772248871 (V153D), rs769113870 (E208D), and rs1318021626 (I222N)]. Analyses of stop-lost, nonsense, UTR, and splice site SNPs were also conducted.

Through integrative bioinformatics analyses, this study has identified the SNPs that are deleterious to the TBX5 protein structure and have the potential to cause HOS. Further wet-lab experiments can validate these findings.

•

Deleterious SNPs in the human TBX5 gene responsible for Holt-Oram Syndrome have been identified.

•

58 missense and 2 nonsense SNPs were identified as deleterious.

•

86 3′ UTR SNPs were predicted to be located on miRNA target sites.

•

Possible effects of missense SNPs on the TBX5 protein structure have been studied.

Upstream ORF frameshift variants in the PAX6 5'UTR cause congenital aniridia

Congenital aniridia (AN) is a severe autosomal dominant panocular disorder associated with pathogenic variants in the PAX6 gene. Previously, we performed a molecular genetic study of a large cohort of Russian patients with AN and revealed four noncoding nucleotide variants in the PAX6 5'UTR. 14 additional PAX6-5'UTR variants were also reported in the literature, but the mechanism of their pathogenicity remained unclear. In the present study, we experimentally analyze five patient-derived PAX6 5'UTR-variants: four variants that we identified in Russian patients (c.-128-2delA, c.-125dupG, c.-122dupG, c.-118_-117del) and one previously reported (c.-52+5G>C). We show that the variants lead to a decrease in the protein translation efficiency, while mRNA expression level is not significantly reduced. Two of these variants also affect splicing. Furthermore, we predict and experimentally validate the presence of an evolutionarily conserved small uORF in the PAX6 5'UTR. All studied variants lead to the frameshift of the uORF, resulting in its extension. This extended out-of-frame uORF overlaps with the downstream CDS and thereby reduces its translation efficiency. We conclude that the uORF frameshift may be the main mechanism of pathogenicity for at least 15 out of 18 known PAX6 5'UTR variants. Moreover, we predict additional uORFs in the PAX6 5'UTR.

Mutations in cis that affect mRNA synthesis, processing and translation

Genetic mutations that cause hereditary diseases usually affect the composition of the transcribed mRNA and its encoded protein, leading to instability of the mRNA and/or the protein. Sometimes, however, such mutations affect the synthesis, the processing or the translation of the mRNA, with similar disastrous effects. We here present an overview of mRNA synthesis, its posttranscriptional modification and its translation into protein. We then indicate which elements in these processes are known to be affected by pathogenic mutations, but we restrict our review to mutations in cis, in the DNA of the gene that encodes the affected protein. These mutations can be in enhancer or promoter regions of the gene, which act as binding sites for transcription factors involved in pre-mRNA synthesis. We also describe mutations in polyadenylation sequences and in splice site regions, exonic and intronic, involved in intron removal. Finally, we include mutations in the Kozak sequence in mRNA, which is involved in protein synthesis. We provide examples of genetic diseases caused by mutations in these DNA regions and refer to databases to help identify these regions. The over-all knowledge of mRNA synthesis, processing and translation is essential for improvement of the diagnosis of patients with genetic diseases.

Determinants of genome-wide distribution and evolution of uORFs in eukaryotes

Upstream open reading frames (uORFs) play widespread regulatory functions in modulating mRNA translation in eukaryotes, but the principles underlying the genomic distribution and evolution of uORFs remain poorly understood. Here, we analyze ~17 million putative canonical uORFs in 478 eukaryotic species that span most of the extant taxa of eukaryotes. We demonstrate how positive and purifying selection, coupled with differences in effective population size (Ne), has shaped the contents of uORFs in eukaryotes. Besides, gene expression level is important in influencing uORF occurrences across genes in a species. Our analyses suggest that most uORFs might play regulatory roles rather than encode functional peptides. We also show that the Kozak sequence context of uORFs has evolved across eukaryotic clades, and that noncanonical uORFs tend to have weaker suppressive effects than canonical uORFs in translation regulation. This study provides insights into the driving forces underlying uORF evolution in eukaryotes.

A novel rare c.-39C>T mutation in the PROS1 5'UTR causing PS deficiency by creating a new upstream translation initiation codon

Autosomal dominant inherited Protein S deficiency (PSD) (MIM 612336) is a rare disorder caused by rare mutations, mainly located in the coding sequence of the structural PROS1 gene, and associated with an increased risk of venous thromboembolism. To identify the molecular defect underlying PSD observed in an extended French pedigree with seven PSD affected members in whom no candidate deleterious PROS1 mutation was detected by Sanger sequencing of PROS1 exons and their flanking intronic regions or via an multiplex ligation-dependent probe amplification (MLPA) approach, a whole genome sequencing strategy was adopted. This led to the identification of a never reported C to T substitution at c.-39 from the natural ATG codon of the PROS1 gene that completely segregates with PSD in the whole family. This substitution ACG→ATG creates a new start codon upstream of the main ATG. We experimentally demonstrated in HeLa cells that the variant generates a novel overlapping upstream open reading frame (uORF) and inhibits the translation of the wild-type PS. This work describes the first example of 5'UTR PROS1 mutation causing PSD through the creation of an uORF, a mutation that is not predicted to be deleterious by standard annotation softwares, and emphasizes the need for better exploration of such type of non-coding variations in clinical genomics.

A Recurrent Rare SOX9 Variant (M469V) is Associated with Congenital Vertebral Malformations

OBJECTIVE

The genetic variations contributed to a substantial proportion of congenital vertebral malformations (CVM). SOX9 gene, a member of the SOX gene family, has been implicated in CVM. To study the SOX9 mutation in CVM patients is of great significance to explain the pathogenesis of scoliosis (the clinical manifestation of CVM) and to explore the pathogenesis of SOX9-related skeletal deformities.

METHODS

A total of 50 singleton patients with CVM were included in this study. Exome Sequencing (ES) was performed on all the patients. The recurrent candidate variant of SOX9 gene was validated by Sanger sequencing. Luciferase assay was performed to investigate the functional changes of this variant.

RESULTS

A recurrent rare heterozygous missense variant in SOX9 gene (NM_000346.3: c.1405A>G, p.M469V) which had not been reported previously was identified in three CVM patients who had the clinical findings of congenital scoliosis without deformities in other systems. This variant was absent from our in-house database and it was predicted to be deleterious (CADD = 24.5). The luciferase assay demonstrated that transactivation capacity of the mutated SOX9 protein was significantly lower than that of the wild-type for the two luciferase reporters (p = 0.0202, p = 0.0082, respectively).

CONCLUSION

This SOX9 mutation (p.M469V) may contribute to CVM without other systematic deformity, which provides important implications and better understanding of phenotypic variability in SOX9-related skeletal deformities.

Some like it translated: small ORFs in the 5'UTR

The 5' untranslated region (5'UTR) is critical in determining post-transcriptional control, which is partly mediated by short upstream open reading frames (uORFs) present in half of mammalian transcripts. uORFs are generally considered to provide functionally important repression of the main-ORF by engaging initiating ribosomes, but under specific environmental conditions such as cellular stress, uORFs can become essential to activate the translation of the main coding sequence. In addition, a growing number of uORF-encoded bioactive microproteins have been described, which have the potential to significantly increase cellular protein diversity. Here we review the diverse cellular contexts in which uORFs play a critical role and discuss the molecular mechanisms underlying their function and regulation. The progress over the last decades in dissecting uORF function suggests that the 5'UTR remains an exciting frontier towards understanding how the cellular proteome is shaped in health and disease.

Characterising the loss-of-function impact of 5' untranslated region variants in 15,708 individuals

Upstream open reading frames (uORFs) are tissue-specific cis-regulators of protein translation. Isolated reports have shown that variants that create or disrupt uORFs can cause disease. Here, in a systematic genome-wide study using 15,708 whole genome sequences, we show that variants that create new upstream start codons, and variants disrupting stop sites of existing uORFs, are under strong negative selection. This selection signal is significantly stronger for variants arising upstream of genes intolerant to loss-of-function variants. Furthermore, variants creating uORFs that overlap the coding sequence show signals of selection equivalent to coding missense variants. Finally, we identify specific genes where modification of uORFs likely represents an important disease mechanism, and report a novel uORF frameshift variant upstream of NF2 in neurofibromatosis. Our results highlight uORF-perturbing variants as an under-recognised functional class that contribute to penetrant human disease, and demonstrate the power of large-scale population sequencing data in studying non-coding variant classes.

Detection of a 4 bp Mutation in the 3'UTR Region of Goat Sox9 Gene and Its Effect on the Growth Traits

Simple Summary

The sex determining region Y (SRY)-type high mobility group (HMG) box 9 (Sox9) gene is critically important in the formation and development of cartilage and is considered the “main regulator” of chondrogenesis. Additionally, a large number of studies have shown that mutations in a single allele of human Sox9 can lead to campomelic dysplasia syndrome. Therefore, the mutations of Sox9 have been the subject of increasing interest among researchers. However, no studies to date have examined the association between Sox9 gene variants and growth traits in goats. Here, we detected a 4 bp indel in the 3′Untranslated Regions (3′UTR) region of Sox9 in Shaanbei white cashmere (SBWC) goats (n = 1109) and studied the association between this indel and growth traits. The 4 bp indel of Sox9 was significantly associated with body length, heart girth, hip width, and all body measurement indexes (p < 0.05) in SBWC goats. Thus, this deletion could be used as an effective molecular marker for maximizing the growth traits of goats in breeding programs.

Abstract

The SRY-type HMG box 9 (Sox9) gene plays an important role in chondrocyte development as well as changes in hypertrophic chondrocytes, indicating that Sox9 can regulate growth in animals. However, no studies to date have examined the correlation between variations in Sox9 and growth traits in goats. Here, we found a 4 bp indel in the 3′UTR of Sox9 and verified its association with growth traits in Shaanbei white cashmere goats (n = 1109). The frequencies of two genotypes (ID and II) were 0.397 and 0.603, respectively, and polymorphic information content (PIC) values showed that the indel had a medium PIC (PIC > 0.25). The 4 bp indel was significantly correlated with body length (p = 0.006), heart girth (p = 0.001), and hip width (p = 4.37 × 10 ⁻⁴). Notably, individuals with the ID genotype had significantly superior phenotypic traits compared with individuals bearing the II genotype. Hence, we speculated that the 4 bp indel is an important mutation affecting growth traits in goat, and may serve as an effective DNA molecular marker for marker-assisted selection in goat breeding programs.

N-Terminal Proteoforms in Human Disease

The collection of chemically different protein variants, or proteoforms, by far exceeds the number of protein-coding genes in the human genome. Major contributors are alternative splicing and protein modifications. In this review, we focus on those proteoforms that differ at their N termini with a molecular link to disease. We describe the main underlying mechanisms that give rise to such N-terminal proteoforms, these being splicing, initiation of protein translation, and protein modifications. Given their role in several human diseases, it is becoming increasingly clear that several of these N-terminal proteoforms may have potential as therapeutic interventions and/or for diagnosing and prognosing their associated disease.

A mutation creating an upstream translation initiation codon in SLC22A5 5'UTR is a frequent cause of primary carnitine deficiency

Primary carnitine deficiency is caused by a defect in the active cellular uptake of carnitine by Na⁺ -dependent organic cation transporter novel 2 (OCTN2). Genetic diagnostic yield for this metabolic disorder has been relatively low, suggesting that disease-causing variants are missed. We Sanger sequenced the 5' untranslated region (UTR) of SLC22A5 in individuals with possible primary carnitine deficiency in whom no or only one mutant allele had been found. We identified a novel 5'-UTR c.-149G>A variant which we characterized by expression studies with reporter constructs in HeLa cells and by carnitine-transport measurements in fibroblasts using a newly developed sensitive assay based on tandem mass spectrometry. This variant, which we identified in 57 of 236 individuals of our cohort, introduces a functional upstream out-of-frame translation initiation codon. We show that the codon suppresses translation from the wild-type ATG of SLC22A5, resulting in reduced OCTN2 protein levels and concomitantly lower transport activity. With an allele frequency of 24.2% the c.-149G>A variant is the most frequent cause of primary carnitine deficiency in our cohort and may explain other reported cases with an incomplete genetic diagnosis. Individuals carrying this variant should be clinically re-evaluated and monitored to determine if this variant has clinical consequences.

Function and Evolution of Upstream ORFs in Eukaryotes

There is growing interest in the role of translational regulation in cellular homeostasis during organismal development. Translation initiation is the rate-limiting step in mRNA translation and is central to translational regulation. Upstream open reading frames (uORFs) are regulatory elements that are prevalent in eukaryotic mRNAs. uORFs modulate the translation initiation rate of downstream coding sequences (CDSs) by sequestering ribosomes. Over the past several years, genome-wide studies have revealed the widespread regulatory functions of uORFs in different species in different biological contexts. Here, we review the current understanding of uORF-mediated translational regulation from the perspective of functional and evolutionary genomics and address remaining gaps that deserve further study.

Mining for missed sORF-encoded peptides

INTRODUCTION

Small open reading frames (sORFs) with potential protein-coding capacity have been disclosed in various transcripts, including long noncoding RNAs (LncRNAs), mRNAs (5'-upstream, coding domain, and 3'-downstream), circular RNAs, pri-miRNAs, and ribosomal RNAs (rRNAs). Recent characterization of several sORF-encoded peptides (SEPs or micropeptides) revealed their important roles in many fundamental biological processes in a broad range of species from yeast to human. The success in the mining of micropeptides attributes to the advanced bioinformatics and high-throughput sequencing techniques. Areas covered: sORFs and SEPs were overlooked for their tiny size and the difficulty of identification by bioinformatics analyses. With more and more sORFs and SEPs have been identified, this field has attracted more attention. This review covers recent advances in the strategies for the detection and identification of sORFs and SEPs. Expert commentary: The advantages and drawbacks of the strategies for detection and identification of sORFs and SEPs are discussed, as well as the techniques that are used to decipher the roles of micropeptides in organisms are described.

Characterization of a family mutation in the 5' untranslated region of the endoglin gene causative of hereditary hemorrhagic telangiectasia

Hereditary hemorrhagic telangiectasia (HHT) is a vascular disease characterized by nose and gastrointestinal bleeding, telangiectases in skin and mucosa, and arteriovenous malformations in major internal organs. Most patients carry a mutation in the coding region of the endoglin (ENG) or activin A receptor type II-1 (ACVRL1) gene. Nonetheless, in around 15% of patients, sequencing analysis and duplication/deletion tests fail to pinpoint mutations in the coding regions of these genes. In these cases, it has been shown that sequencing of the 5’-untranslated region (5’UTR) of ENG may be useful to identify novel mutations in the ENG non-coding region. Here we report the genetic characterization and functional analysis of the heterozygous mutation c.-142A>T in the 5’UTR region of ENG found in a family with several members affected by HHT. This variant gives rise to a new initiation codon of the protein that involves the change in its open reading frame. Transfection studies in monkey cells using endoglin expression vectors demonstrated that c-142A>T mutation results in a clear reduction in the levels of the endoglin protein. These results support the inclusion of the 5’UTR of ENG in the standard genetic testing for HHT to increase its sensitivity.

Genome-wide maps of ribosomal occupancy provide insights into adaptive evolution and regulatory roles of uORFs during Drosophila development

Upstream open reading frames (uORFs) play important roles in regulating the main coding DNA sequences (CDSs) via translational repression. Despite their prevalence in the genomes, uORFs are overall discriminated against by natural selection. However, it remains unclear why in the genomes there are so many uORFs more conserved than expected under the assumption of neutral evolution. Here, we generated genome-wide maps of translational efficiency (TE) at the codon level throughout the life cycle of Drosophila melanogaster. We identified 35,735 uORFs that were expressed, and 32,224 (90.2%) of them showed evidence of ribosome occupancy during Drosophila development. The ribosome occupancy of uORFs is determined by genomic features, such as optimized sequence contexts around their start codons, a shorter distance to CDSs, and higher coding potentials. Our population genomic analysis suggests the segregating mutations that create or disrupt uORFs are overall deleterious in D. melanogaster. However, we found for the first time that many (68.3% of) newly fixed uORFs that are associated with ribosomes in D. melanogaster are driven by positive Darwinian selection. Our findings also suggest that uORFs play a vital role in controlling the translational program in Drosophila. Moreover, we found that many uORFs are transcribed or translated in a developmental stage-, sex-, or tissue-specific manner, suggesting that selective transcription or translation of uORFs could potentially modulate the TE of the downstream CDSs during Drosophila development.

Upstream open reading frames (uORFs) in the 5′ untranslated regions (UTRs) of messenger RNAs can potentially inhibit translation of the downstream regions that encode proteins by sequestering protein-making machinery the ribosome. Moreover, mutations that destroy existing uORFs or create new ones are known to cause human disease. Although mutations that create new uORFs are generally deleterious and are selected against, many uORFs are evolutionarily conserved across eukaryotic species. To resolve this dilemma, we used extensive mRNA-Seq and ribosome profiling to generate high-resolution genome-wide maps of ribosome occupancy and translational efficiency (TE) during the life cycle of the fruit fly D. melanogaster. This allowed us to identify the sequence features of uORFs that influence their ability to associate with ribosomes. We demonstrate for the first time that the majority of the newly fixed uORFs in D. melanogaster, especially the translated ones, are under positive Darwinian selection. We also show that uORFs exert widespread repressive effects on the translation of the downstream protein-coding region. We find that many uORFs are transcribed or translated in a developmental stage-, sex-, or tissue-specific manner. Our results suggest that during Drosophila development, changes in the TE of uORFs, as well as the inclusion/exclusion of uORFs, are frequently exploited to inversely influence the translation of the downstream protein-coding regions. Our study provides novel insights into the molecular mechanisms and functional consequences of uORF-mediated regulation.

In Search of Lost Small Peptides

A large body of evidence indicates that genome annotation pipelines have biased our view of coding sequences because they generally undersample small proteins and peptides. The recent development of genome-wide translation profiling reveals the prevalence of small/short open reading frames (smORFs or sORFs), which are scattered over all classes of transcripts, including both mRNAs and presumptive long noncoding RNAs. Proteomic approaches further confirm an unexpected variety of smORF-encoded peptides (SEPs), representing an overlooked reservoir of bioactive molecules. Indeed, functional studies in a broad range of species from yeast to humans demonstrate that SEPs can harbor key activities for the control of development, differentiation, and physiology. Here we summarize recent advances in the discovery and functional characterization of smORF/SEPs and discuss why these small players can no longer be ignored with regard to genome function.