Natural variability of Kozak sequences correlates with function in a zebrafish model

Mutational cooperativity of RUNX1::RUNX1T1 isoform 9a and oncogenic NRAS in zebrafish myeloid leukaemia

RUNX1::RUNX1T1 (R::RT1) acute myeloid leukaemia (AML) remains a clinical challenge, and further research is required to model and understand leukaemogenesis. Previous zebrafish R::RT1 models were hampered by embryonic lethality and low penetrance of the malignant phenotype. Here, we overcome this by developing an adult zebrafish model in which the human R::RT1 isoform 9a is co-expressed with the frequently co-occurring oncogenic NRASG12D mutation in haematopoietic stem and progenitor cells (HSPCs), using the Runx1+23 enhancer. Approximately 50% of F0 9a+NRASG12D transgenic zebrafish developed signs of haematological disease between 5 and 14 months, with 27% exhibiting AML-like pathology: myeloid precursor expansion, erythrocyte reduction, kidney marrow hypercellularity and the presence of blasts. Moreover, only 9a+NRASG12D transplant recipients developed leukaemia with high rates of mortality within 40 days, inferring the presence of leukaemia stem cells. These leukaemic features were rare or not observed in animals expressing either the NRAS or 9a oncogenes alone, suggesting 9a and NRAS cooperation drives leukaemogenesis. This novel adult AML zebrafish model provides a powerful new tool for investigating the basis of R::RT1 - NRAS cooperativity with the potential to uncover new therapeutic targets.

U6 snRNA m6A modification is required for accurate and efficient splicing of C. elegans and human pre-mRNAs

pre-mRNA splicing is a critical feature of eukaryotic gene expression. Both cis- and trans-splicing rely on accurately recognising splice site sequences by spliceosomal U snRNAs and associated proteins. Spliceosomal snRNAs carry multiple RNA modifications with the potential to affect different stages of pre-mRNA splicing. Here, we show that the conserved U6 snRNA m6A methyltransferase METT-10 is required for accurate and efficient cis- and trans-splicing of C. elegans pre-mRNAs. The absence of METT-10 in C. elegans and METTL16 in humans primarily leads to alternative splicing at 5' splice sites with an adenosine at +4 position. In addition, METT-10 is required for splicing of weak 3' cis- and trans-splice sites. We identified a significant overlap between METT-10 and the conserved splicing factor SNRNP27K in regulating 5' splice sites with +4A. Finally, we show that editing endogenous 5' splice site +4A positions to +4U restores splicing to wild-type positions in a mett-10 mutant background, supporting a direct role for U6 snRNA m6A modification in 5' splice site recognition. We conclude that the U6 snRNA m6A modification is important for accurate and efficient pre-mRNA splicing.

UPF1 regulates mRNA stability by sensing poorly translated coding sequences

Post-transcriptional mRNA regulation shapes gene expression, yet how cis-elements and mRNA translation interface to regulate mRNA stability is poorly understood. We find that the strength of translation initiation, upstream open reading frame (uORF) content, codon optimality, AU-rich elements, microRNA binding sites, and open reading frame (ORF) length function combinatorially to regulate mRNA stability. Machine-learning analysis identifies ORF length as the most important conserved feature regulating mRNA decay. We find that Upf1 binds poorly translated and untranslated ORFs, which are associated with a higher decay rate, including mRNAs with uORFs and those with exposed ORFs after stop codons. Our study emphasizes Upf1's converging role in surveilling mRNAs with exposed ORFs that are poorly translated, such as mRNAs with long ORFs, ORF-like 3' UTRs, and mRNAs containing uORFs. We propose that Upf1 regulation of poorly/untranslated ORFs provides a unifying mechanism of surveillance in regulating mRNA stability and homeostasis in an exon-junction complex (EJC)-independent nonsense-mediated decay (NMD) pathway that we term ORF-mediated decay (OMD).

Plxnd1-mediated mechanosensing of blood flow controls the caliber of the Dorsal Aorta via the transcription factor Klf2

The cardiovascular system generates and responds to mechanical forces. The heartbeat pumps blood through a network of vascular tubes, which adjust their caliber in response to the hemodynamic environment. However, how endothelial cells in the developing vascular system integrate inputs from circulatory forces into signaling pathways to define vessel caliber is poorly understood. Using vertebrate embryos and in vitro-assembled microvascular networks of human endothelial cells as models, flow and genetic manipulations, and custom software, we reveal that Plexin-D1, an endothelial Semaphorin receptor critical for angiogenic guidance, employs its mechanosensing activity to serve as a crucial positive regulator of the Dorsal Aorta's (DA) caliber. We also uncover that the flow-responsive transcription factor KLF2 acts as a paramount mechanosensitive effector of Plexin-D1 that enlarges endothelial cells to widen the vessel. These findings illuminate the molecular and cellular mechanisms orchestrating the interplay between cardiovascular development and hemodynamic forces.

The regulatory landscape of 5' UTRs in translational control during zebrafish embryogenesis

The 5' UTRs of mRNAs are critical for translation regulation, but their in vivo regulatory features are poorly characterized. Here, we report the regulatory landscape of 5' UTRs during early zebrafish embryogenesis using a massively parallel reporter assay of 18,154 sequences coupled to polysome profiling. We found that the 5' UTR is sufficient to confer temporal dynamics to translation initiation, and identified 86 motifs enriched in 5' UTRs with distinct ribosome recruitment capabilities. A quantitative deep learning model, DaniO5P, revealed a combined role for 5' UTR length, translation initiation site context, upstream AUGs and sequence motifs on in vivo ribosome recruitment. DaniO5P predicts the activities of 5' UTR isoforms and indicates that modulating 5' UTR length and motif grammar contributes to translation initiation dynamics. This study provides a first quantitative model of 5' UTR-based translation regulation in early vertebrate development and lays the foundation for identifying the underlying molecular effectors.

A non-coding variant in the Kozak sequence of RARS2 strongly decreases protein levels and causes pontocerebellar hypoplasia

Bi-allelic variants in the mitochondrial arginyl-transfer RNA synthetase (RARS2) gene have been involved in early-onset encephalopathies classified as pontocerebellar hypoplasia (PCH) type 6 and in epileptic encephalopathy. A variant (NM_020320.3:c.-2A > G) in the promoter and 5'UTR of the RARS2 gene has been previously identified in a family with PCH. Only a mild impact of this variant on the mRNA level has been detected. As RARS2 is non-dosage-sensitive, this observation is not conclusive in regard of the pathogenicity of the variant.We report and describe here a new patient with the same variant in the RARS2 gene, at the homozygous state. This patient presents with a clinical phenotype consistent with PCH6 although in the absence of lactic acidosis. In agreement with the previous study, we measured RARS2 mRNA levels in patient's fibroblasts and detected a partially preserved gene expression compared to control. Importantly, this variant is located in the Kozak sequence that controls translation initiation. Therefore, we investigated the impact on protein translation using a bioinformatic approach and western blotting. We show here that this variant, additionally to its effect on the transcription, also disrupts the consensus Kozak sequence, and has a major impact on RARS2 protein translation. Through the identification of this additional case and the characterization of the molecular consequences, we clarified the involvement of this Kozak variant in PCH and on protein synthesis. This work also points to the current limitation in the pathogenicity prediction of variants located in the translation initiation region.

Viral proteogenomic and expression profiling during productive replication of a skin-tropic herpesvirus in the natural host

Efficient transmission of herpesviruses is essential for dissemination in host populations; however, little is known about the viral genes that mediate transmission, mostly due to a lack of natural virus-host model systems. Marek's disease is a devastating herpesviral disease of chickens caused by Marek's disease virus (MDV) and an excellent natural model to study skin-tropic herpesviruses and transmission. Like varicella zoster virus that causes chicken pox in humans, the only site where infectious cell-free MD virions are efficiently produced is in epithelial skin cells, a requirement for host-to-host transmission. Here, we enriched for heavily infected feather follicle epithelial skin cells of live chickens to measure both viral transcription and protein expression using combined short- and long-read RNA sequencing and LC/MS-MS bottom-up proteomics. Enrichment produced a previously unseen breadth and depth of viral peptide sequencing. We confirmed protein translation for 84 viral genes at high confidence (1% FDR) and correlated relative protein abundance with RNA expression levels. Using a proteogenomic approach, we confirmed translation of most well-characterized spliced viral transcripts and identified a novel, abundant isoform of the 14 kDa transcript family via IsoSeq transcripts, short-read intron-spanning sequencing reads, and a high-quality junction-spanning peptide identification. We identified peptides representing alternative start codon usage in several genes and putative novel microORFs at the 5' ends of two core herpesviral genes, pUL47 and ICP4, along with strong evidence of independent transcription and translation of the capsid scaffold protein pUL26.5. Using a natural animal host model system to examine viral gene expression provides a robust, efficient, and meaningful way of validating results gathered from cell culture systems.

TP53BP1, a dual-coding gene, uses promoter switching and translational reinitiation to express a smORF protein

The complexity of the metazoan proteome is significantly increased by the expression of small proteins (<100 aa) derived from smORFs within lncRNAs, uORFs, 3' UTRs and, reading frames overlapping the CDS. These smORF encoded proteins (SEPs) have diverse roles, ranging from the regulation of cellular physiological to essential developmental functions. We report the characterization of a new member of this protein family, SEP53BP1, derived from a small internal ORF that overlaps the CDS encoding 53BP1. Its expression is coupled to the utilization of an alternative, cell-type specific promoter coupled to translational reinitiation events mediated by a uORF in the alternative 5' TL of the mRNA. This uORF-mediated reinitiation at an internal ORF is also observed in zebrafish. Interactome studies indicate that the human SEP53BP1 associates with components of the protein turnover pathway including the proteasome, and the TRiC/CCT chaperonin complex, suggesting that it may play a role in cellular proteostasis.

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.

Translational enhancement by base editing of the Kozak sequence rescues haploinsufficiency

A variety of single-gene human diseases are caused by haploinsufficiency, a genetic condition by which mutational inactivation of one allele leads to reduced protein levels and functional impairment. Translational enhancement of the spare allele could exert a therapeutic effect. Here we developed BOOST, a novel gene-editing approach to rescue haploinsufficiency loci by the change of specific single nucleotides in the Kozak sequence, which controls translation by regulating start codon recognition. We evaluated for translational strength 230 Kozak sequences of annotated human haploinsufficient genes and 4621 derived variants, which can be installed by base editing, by a high-throughput reporter assay. Of these variants, 149 increased the translation of 47 Kozak sequences, demonstrating that a substantial proportion of haploinsufficient genes are controlled by suboptimal Kozak sequences. Validation of 18 variants for 8 genes produced an average enhancement in an expression window compatible with the rescue of the genetic imbalance. Base editing of the NCF1 gene, whose monoallelic loss causes chronic granulomatous disease, resulted in the desired increase of NCF1 (p47phox) protein levels in a relevant cell model. We propose BOOST as a fine-tuned approach to modulate translation, applicable to the correction of dozens of haploinsufficient monogenic disorders independently of the causing mutation.

The Oldest Co-opted gag Gene of a Human Endogenous Retrovirus Shows Placenta-Specific Expression and Is Upregulated in Diffuse Large B-Cell Lymphomas

Vertebrate genomes contain endogenous retroviruses (ERVs) that represent remnants of past germline infections by ancient retroviruses. Despite comprising 8% of the human genome, the human ERVs (HERVs) do not encode a replication competent retrovirus. However, some HERV genes have been co-opted to serve host functions, most notably the viral envelope-derived syncytins involved in placentation. Here, we identify the oldest HERV intact gag gene with an open reading frame, gagV1. Its provirus contains an intact env, envV1, and the first open reading frame found in an HERV gag leader, pre-gagV1, which encodes a novel protein. This HERV is linked to a related gag gene, gagV3, and these three genes all show patterns of evolutionary conservation in primates. gagV1 and pre-gagV1 orthologs are present in all simian primate lineages indicating that this HERV entered the germline of the common simian primate ancestor at least 43 Ma, whereas gagV3 is found in Old and New World monkeys. gagV1 and gagV3 have undergone recurrent gene conversion events and positive selection. Expression of gagV1, gagV3, and pre-gagV1 is restricted to the placenta in humans and macaques suggesting co-option for placenta-specific host functions. Transcriptomic analysis of human tumors also found upregulated levels of gagV1 transcripts in diffuse large B-cell lymphomas. These findings suggest that these HERV-V genes may be useful markers for the most common type of non-Hodgkin's lymphoma and that they may have contributed to the successive domestications of env and gag genes in eutherians involved in the ongoing ERV-driven evolution of the placenta.

A Toolbox for Efficient Proximity-Dependent Biotinylation in Zebrafish Embryos

Understanding how proteins are organized in compartments is essential to elucidating their function. While proximity-dependent approaches such as BioID have enabled a massive increase in information about organelles, protein complexes, and other structures in cell culture, to date there have been only a few studies on living vertebrates. Here, we adapted proximity labeling for protein discovery in vivo in the vertebrate model organism, zebrafish. Using lamin A (LMNA) as bait and green fluorescent protein (GFP) as a negative control, we developed, optimized, and benchmarked in vivo TurboID and miniTurbo labeling in early zebrafish embryos. We developed both an mRNA injection protocol and a transgenic system in which transgene expression is controlled by a heat shock promoter. In both cases, biotin is provided directly in the egg water, and we demonstrate that 12 h of labeling are sufficient for biotinylation of prey proteins, which should permit time-resolved analysis of development. After statistical scoring, we found that the proximal partners of LMNA detected in each system were enriched for nuclear envelope and nuclear membrane proteins and included many orthologs of human proteins identified as proximity partners of lamin A in mammalian cell culture. The tools and protocols developed here will allow zebrafish researchers to complement genetic tools with powerful proteomics approaches.

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.

ORFik: a comprehensive R toolkit for the analysis of translation

BACKGROUND

With the rapid growth in the use of high-throughput methods for characterizing translation and the continued expansion of multi-omics, there is a need for back-end functions and streamlined tools for processing, analyzing, and characterizing data produced by these assays.

RESULTS

Here, we introduce ORFik, a user-friendly R/Bioconductor API and toolbox for studying translation and its regulation. It extends GenomicRanges from the genome to the transcriptome and implements a framework that integrates data from several sources. ORFik streamlines the steps to process, analyze, and visualize the different steps of translation with a particular focus on initiation and elongation. It accepts high-throughput sequencing data from ribosome profiling to quantify ribosome elongation or RCP-seq/TCP-seq to also quantify ribosome scanning. In addition, ORFik can use CAGE data to accurately determine 5'UTRs and RNA-seq for determining translation relative to RNA abundance. ORFik supports and calculates over 30 different translation-related features and metrics from the literature and can annotate translated regions such as proteins or upstream open reading frames (uORFs). As a use-case, we demonstrate using ORFik to rapidly annotate the dynamics of 5' UTRs across different tissues, detect their uORFs, and characterize their scanning and translation in the downstream protein-coding regions.

CONCLUSION

In summary, ORFik introduces hundreds of tested, documented and optimized methods. ORFik is designed to be easily customizable, enabling users to create complete workflows from raw data to publication-ready figures for several types of sequencing data. Finally, by improving speed and scope of many core Bioconductor functions, ORFik offers enhancement benefiting the entire Bioconductor environment.

AVAILABILITY

http://bioconductor.org/packages/ORFik .

Profiling of Small Ribosomal Subunits Reveals Modes and Regulation of Translation Initiation

Translation initiation is often attributed as the rate-determining step of eukaryotic protein synthesis and key to gene expression control. Despite this centrality, the series of steps involved in this process is poorly understood. Here, we capture the transcriptome-wide occupancy of ribosomes across all stages of translation initiation, enabling us to characterize the transcriptome-wide dynamics of ribosome recruitment to mRNAs, scanning across 5' UTRs and stop codon recognition, in a higher eukaryote. We provide mechanistic evidence for ribosomes attaching to the mRNA by threading the mRNA through the small subunit. Moreover, we identify features that regulate the recruitment and processivity of scanning ribosomes and redefine optimal initiation contexts. Our approach enables deconvoluting translation initiation into separate stages and identifying regulators at each step.

An analysis of the Kozak consensus in retinal genes and its relevance to gene therapy

PURPOSE

The classic Kozak consensus is a critical genetic element included in gene therapy transgenes to encourage the translation of the therapeutic coding sequence. Despite optimizations of other transgene elements, the Kozak consensus has not yet been considered for potential tissue-specific sequence refinement. We screened the -9 to -1 region relative to the AUG start codon of retina-specific genes to identify whether a Kozak consensus that is different from the classic sequence may be more appropriate for inclusion in gene therapy transgenes that treat inherited retinal disease.

METHODS

Sequences for 135 genes known to cause nonsyndromic inherited retinal disease were extracted from the NCBI database, and the -9 to -1 nucleotides were compared. This panel was then refined to 75 genes with specific retinal functions, for which the -9 to -1 nucleotides were placed in front of a GFP transcript sequence and RNAfold predictions performed. These were compared with a GFP sequence with the classic Kozak consensus (GCCGCCACC), and sequences from retinal genes with minimum free energy (MFE) predictions greater than the reference sequence were selected to generate an optimized Kozak consensus sequence. The original Kozak consensus and the refined retina Kozak consensus were placed upstream of the Renilla luciferase coding sequence, which were used to transfect retinoblastoma cell lines Y-79 and WERI-RB-1 and HEK 293T/17 cells.

RESULTS

The nucleotide frequencies of the original panel of genes were determined to be comparable to the classic Kozak consensus. RNAfold analysis of a GFP transcript with the classic Kozak sequence in the 5' untranslated region (UTR) generated an MFE prediction of -503.3 kcal/mol. RNAfold analysis was then performed with a GFP transcript containing each -9 to -1 Kozak sequence of 75 retinal genes. Thirty-eight of the 75 genes provided a greater MFE value than -503.3 kcal/mol and exhibited an absence of stable secondary structures before the AUG codon. The -9 to -1 nucleotide frequencies of these genes identified a Kozak consensus of ACCGAGACC, differing from the classic Kozak consensus at positions -9, -5, and -4. Applying this sequence to the GFP transcript increased the MFE prediction to -500.1 kcal/mol. The newly identified retina Kozak sequence was also applied to Renilla luciferase plus the REP1 and RPGR transcripts used in current clinical trials. In all examples, the predicted transcript MFE score increased when compared with the current transcript sequences containing classic Kozak consensus sequences. In vitro transfections identified a 7%-9% increase in Renilla activity when incorporating the optimized Kozak sequence.

CONCLUSIONS

The Kozak consensus is a critical element of eukaryotic genes; therefore, it is a required feature of gene therapy transgenes. To date, the classic sequence of GCCRCC (-6 to -1) has typically been incorporated in gene therapy transgenes, but the analysis described here suggests that, for vectors targeting the retina, using a Kozak consensus derived from retinal genes can provide increased expression of the target product.

TAEL 2.0: An Improved Optogenetic Expression System for Zebrafish

Inducible gene expression systems are valuable tools for studying biological processes. We previously developed an optogenetic gene expression system called TAEL that is optimized for use in zebrafish. When illuminated with blue light, TAEL transcription factors dimerize and activate gene expression downstream of the TAEL-responsive C120 promoter. By using light as the inducing agent, the TAEL/C120 system overcomes limitations of traditional inducible expression systems by enabling fine spatial and temporal regulation of gene expression. In this study, we describe ongoing efforts to improve the TAEL/C120 system. We made modifications to both the TAEL transcriptional activator and the C120 regulatory element, collectively referred to as TAEL 2.0. We demonstrate that TAEL 2.0 consistently induces higher levels of reporter gene expression and at a faster rate, but with comparable background and toxicity as the original TAEL system. With these improvements, we were able to create functional stable transgenic lines to express the TAEL 2.0 transcription factor either ubiquitously or with a tissue-specific promoter. We demonstrate that the ubiquitous line in particular can be used to induce expression at late embryonic and larval stages, addressing a major deficiency of the original TAEL system. This improved optogenetic expression system will be a broadly useful resource for the zebrafish community.

Cell-intrinsic Fgf signaling contributes to primordial germ cell homing in zebrafish

Primordial germ cells (PGCs) are specified before gastrulation and migrate toward the developing gonads. Previous in vitro studies have demonstrated a cell-intrinsic requirement of fibroblast growth factors (FGFs) by PGCs; however, no evidence suggests FGFs signal directly to PGCs in vivo. Here, using zebrafish as the animal model, we identified the mRNA expressions of Fgf receptors (Fgfrs) and determined the roles of Fgf signaling in migrating PGCs. To clarify the functions of Fgf signaling, we manipulated Fgf signaling specifically in PGCs using dominant-negative (dn) and constitutively-active (ca) Fgfrs and revealed a requirement of a basal Fgf signaling level for the robust arrival of PGCs. Repression of Fgf signaling in PGCs swayed the marginal positioning of PGCs as early as 6 h post-fertilization (6 hpf) and disrupted their arrival at the gonadal ridge at 24 hpf. On the other hand, the ectopic PGC phenotypes caused by the dn-Fgfrs could be alleviated by constitutive activation of Fgf signaling. In addition, we carefully ruled out the somatic effects in mosaic embryos by injecting RNA materials into one blastomere of the four- or eight-cell stage embryos. Injection of dn-Fgfrs into one of eight blastomeres hampered the arrival of only the treated PGCs, while the other PGCs remained unaffected. Furthermore, mosaic treatment of ca-Fgfrs rescued the ectopic rates of dn-Fgfr treated PGCs, while the other PGCs remained more ectopic within the same embryos. Interestingly, PGC-specific repression of Fgf signaling did not compromise the PGC number. To our knowledge, this is the first in vivo evidence to show that Fgf signaling plays a cell-intrinsic role in the migration of vertebrate PGCs.

Reconstruction of the birth of a male sex chromosome present in Atlantic herring

Understanding the evolution of sex determination mechanisms and sex chromosomes is of fundamental importance in biology. Here we have reconstructed the evolution of the sex-determining region in the Atlantic herring. The region is small and contains only three intact genes. The candidate sex-determining factor BMPR1BBY is an evolutionary innovation in the herring lineage. It encodes a truncated form of a BMP type I receptor, which originated by gene duplication and underwent rapid protein evolution. The receptor has maintained its kinase activity and has the potential to induce development of testis. The other two genes in the sex-determining region, CATSPERG and CATSPER3Y, are male beneficial genes because they encode proteins predicted to be essential for sperm to fertilize the egg.

The mechanisms underlying sex determination are astonishingly plastic. Particularly the triggers for the molecular machinery, which recalls either the male or female developmental program, are highly variable and have evolved independently and repeatedly. Fish show a huge variety of sex determination systems, including both genetic and environmental triggers. The advent of sex chromosomes is assumed to stabilize genetic sex determination. However, because sex chromosomes are notoriously cluttered with repetitive DNA and pseudogenes, the study of their evolution is hampered. Here we reconstruct the birth of a Y chromosome present in the Atlantic herring. The region is tiny (230 kb) and contains only three intact genes. The candidate male-determining gene BMPR1BBY encodes a truncated form of a BMP1B receptor, which originated by gene duplication and translocation and underwent rapid protein evolution. BMPR1BBY phosphorylates SMADs in the absence of ligand and thus has the potential to induce testis formation. The Y region also contains two genes encoding subunits of the sperm-specific Ca²⁺ channel CatSper required for male fertility. The herring Y chromosome conforms with a characteristic feature of many sex chromosomes, namely, suppressed recombination between a sex-determining factor and genes that are beneficial for the given sex. However, the herring Y differs from other sex chromosomes in that suppression of recombination is restricted to an ∼500-kb region harboring the male-specific and sex-associated regions. As a consequence, any degeneration on the herring Y chromosome is restricted to those genes located in the small region affected by suppressed recombination.

Knockout of family with sequence similarity 170 member A (Fam170a) causes male subfertility, while Fam170b is dispensable in mice†

Families with sequence similarity 170 members A and B (FAM170A and FAM170B) are testis-specific, paralogous proteins that share 31% amino acid identity and are conserved throughout mammals. While previous in vitro experiments suggested that FAM170B, an acrosome-localized protein, plays a role in the mouse sperm acrosome reaction and fertilization, the role of FAM170A in the testis has not been explored. In this study, we used CRISPR/Cas9 to generate null alleles for each gene, and homozygous null (-/-) male mice were mated to wild-type females for 6 months to assess fertility. Fam170b-/- males were found to produce normal litter sizes and had normal sperm counts, motility, and sperm morphology. In contrast, mating experiments revealed significantly reduced litter sizes and a reduced pregnancy rate from Fam170a-/- males compared with controls. Fam170a-/-;Fam170b-/- double knockout males also produced markedly reduced litter sizes, although not significantly different from Fam170a-/- alone, suggesting that Fam170b does not compensate for the absence of Fam170a. Fam170a-/- males exhibited abnormal spermiation, abnormal head morphology, and reduced progressive sperm motility. Thus, FAM170A has an important role in male fertility, as the loss of the protein leads to subfertility, while FAM170B is expendable. The molecular functions of FAM170A in spermatogenesis are as yet unknown; however, the protein localizes to the nucleus of elongating spermatids and may mediate its effects on spermatid head shaping and spermiation by regulating the expression of other genes. This work provides the first described role of FAM170A in reproduction and has implications for improving human male infertility diagnoses.

Fluorescent protein expression in the ectomycorrhizal fungus Laccaria bicolor: a plasmid toolkit for easy use of fluorescent markers in basidiomycetes

For long time, studies on ectomycorrhiza (ECM) have been limited by inefficient expression of fluorescent proteins (FPs) in the fungal partner. To convert this situation, we have evaluated the basic requirements of FP expression in the model ECM homobasidiomycete Laccaria bicolor and established eGFP and mCherry as functional FP markers. Comparison of intron-containing and intronless FP-expression cassettes confirmed that intron-processing is indispensable for efficient FP expression in Laccaria. Nuclear FP localization was obtained via in-frame fusion of FPs between the intron-containing genomic gene sequences of Laccaria histone H2B, while cytosolic FP expression was produced by incorporating the intron-containing 5' fragment of the glyceraldehyde-3-phosphate dehydrogenase encoding gene. In addition, we have characterized the consensus Kozak sequence of strongly expressed genes in Laccaria and demonstrated its boosting effect on transgene mRNA accumulation. Based on these results, an Agrobacterium-mediated transformation compatible plasmid set was designed for easy use of FPs in Laccaria. The four cloning plasmids presented here allow fast and highly flexible construction of C-terminal in-frame fusions between the sequences of interest and the two FPs, expressed either from the endogenous gene promoter, allowing thus evaluation of the native regulation modes of the gene under study, or alternatively, from the constitutive Agaricus bisporus gpdII promoter for enhanced cellular protein localization assays. The molecular tools described here for cell-biological studies in Laccaria can also be exploited in studies of other biotrophic or saprotrophic basidiomycete species susceptible to genetic transformation.

Rational Design of Bioavailable Photosensitizers for Manipulation and Imaging of Biological Systems

Light-mediated chemical reactions are powerful methods for manipulating and interrogating biological systems. Photosensitizers, compounds that generate reactive oxygen species upon excitation with light, can be utilized for numerous biological experiments, but the repertoire of bioavailable photosensitizers is limited. Here, we describe the synthesis, characterization, and utility of two photosensitizers based upon the widely used rhodamine scaffold and demonstrate their efficacy for chromophore-assisted light inactivation, cell ablation in culture and in vivo, and photopolymerization of diaminobenzidine for electron microscopy. These chemical tools will facilitate a broad range of applications spanning from targeted destruction of proteins to high-resolution imaging.

Structural Insights into the Mammalian Late-Stage Initiation Complexes

In higher eukaryotes, the mRNA sequence in the direct vicinity of the start codon, called the Kozak sequence (CRCCaugG, where R is a purine), is known to influence the rate of the initiation process. However, the molecular basis underlying its role remains poorly understood. Here, we present the cryoelectron microscopy (cryo-EM) structures of mammalian late-stage 48S initiation complexes (LS48S ICs) in the presence of two different native mRNA sequences, β-globin and histone 4, at overall resolution of 3 and 3.5 Å, respectively. Our high-resolution structures unravel key interactions from the mRNA to eukaryotic initiation factors (eIFs): 1A, 2, 3, 18S rRNA, and several 40S ribosomal proteins. In addition, we are able to study the structural role of ABCE1 in the formation of native 48S ICs. Our results reveal a comprehensive map of ribosome/eIF-mRNA and ribosome/eIF-tRNA interactions and suggest the impact of mRNA sequence on the structure of the LS48S IC.

Conservation and Variability of the AUG Initiation Codon Context in Eukaryotes

Selection of the translation initiation site (TIS) is a crucial step during translation. In the 1980s Marylin Kozak performed key studies on vertebrate mRNAs to characterize the optimal TIS consensus sequence, the Kozak motif. Within this motif, conservation of nucleotides in crucial positions, namely a purine at -3 and a G at +4 (where the A of the AUG is numbered +1), is essential for TIS recognition. Ever since its characterization the Kozak motif has been regarded as the optimal sequence to initiate translation in all eukaryotes. We revisit here published in silico data on TIS consensus sequences, as well as experimental studies from diverse eukaryotic lineages, and propose that, while the -3A/G position is universally conserved, the remaining variability of the consensus sequences enables their classification as optimal, strong, and moderate TIS sequences.

uORF-mediated translational control: recently elucidated mechanisms and implications in cancer

Protein synthesis is tightly regulated, and its dysregulation can contribute to the pathology of various diseases, including cancer. Increased or selective translation of mRNAs can promote cancer cell proliferation, metastasis and tumor expansion. Translational control is one of the most important means for cells to quickly adapt to environmental stresses. Adaptive translation involves various alternative mechanisms of translation initiation. Upstream open reading frames (uORFs) serve as a major regulator of stress-responsive translational control. Since recent advances in omics technologies including ribo-seq have expanded our knowledge of translation, we discuss emerging mechanisms for uORF-mediated translation regulation and its impact on cancer cell biology. A better understanding of dysregulated translational control of uORFs in cancer would facilitate the development of new strategies for cancer therapy.

Combining near-infrared fluorescence with Brainbow to visualize expression of specific genes within a multicolor context

Fluorescent proteins are a powerful experimental tool, allowing the visualization of gene expression and cellular behaviors in a variety of systems. Multicolor combinations of fluorescent proteins, such as Brainbow, have expanded the range of possible research questions and are useful for distinguishing and tracking cells. The addition of a separately driven color, however, would allow researchers to report expression of a manipulated gene within the multicolor context to investigate mechanistic effects. A far-red or near-infrared protein could be particularly suitable in this context, as these can be distinguished spectrally from Brainbow. We investigated five far-red/near-infrared proteins in zebrafish: TagRFP657, mCardinal, miRFP670, iRFP670, and mIFP. Our results show that both mCardinal and iRFP670 are useful fluorescent proteins for zebrafish expression. We also introduce a new transgenic zebrafish line that expresses Brainbow under the control of the neuroD promoter. We demonstrate that mCardinal can be used to track the expression of a manipulated bone morphogenetic protein receptor within the Brainbow context. The overlay of near-infrared fluorescence onto a Brainbow background defines a clear strategy for future research questions that aim to manipulate or track the effects of specific genes within a population of cells that are delineated using multicolor approaches.

Redox Signaling by Reactive Electrophiles and Oxidants

The concept of cell signaling in the context of nonenzyme-assisted protein modifications by reactive electrophilic and oxidative species, broadly known as redox signaling, is a uniquely complex topic that has been approached from numerous different and multidisciplinary angles. Our Review reflects on five aspects critical for understanding how nature harnesses these noncanonical post-translational modifications to coordinate distinct cellular activities: (1) specific players and their generation, (2) physicochemical properties, (3) mechanisms of action, (4) methods of interrogation, and (5) functional roles in health and disease. Emphasis is primarily placed on the latest progress in the field, but several aspects of classical work likely forgotten/lost are also recollected. For researchers with interests in getting into the field, our Review is anticipated to function as a primer. For the expert, we aim to stimulate thought and discussion about fundamentals of redox signaling mechanisms and nuances of specificity/selectivity and timing in this sophisticated yet fascinating arena at the crossroads of chemistry and biology.

Adenine Enrichment at the Fourth CDS Residue in Bacterial Genes Is Consistent with Error Proofing for +1 Frameshifts

Beyond selection for optimal protein functioning, coding sequences (CDSs) are under selection at the RNA and DNA levels. Here, we identify a possible signature of “dual-coding,” namely extensive adenine (A) enrichment at bacterial CDS fourth sites. In 99.07% of studied bacterial genomes, fourth site A use is greater than expected given genomic A-starting codon use. Arguing for nucleotide level selection, A-starting serine and arginine second codons are heavily utilized when compared with their non-A starting synonyms. Several models have the ability to explain some of this trend. In part, A-enrichment likely reduces 5′ mRNA stability, promoting translation initiation. However T/U, which may also reduce stability, is avoided. Further, +1 frameshifts on the initiating ATG encode a stop codon (TGA) provided A is the fourth residue, acting either as a frameshift “catch and destroy” or a frameshift stop and adjust mechanism and hence implicated in translation initiation. Consistent with both, genomes lacking TGA stop codons exhibit weaker fourth site A-enrichment. Sequences lacking a Shine–Dalgarno sequence and those without upstream leader genes, that may be more error prone during initiation, have greater utilization of A, again suggesting a role in initiation. The frameshift correction model is consistent with the notion that many genomic features are error-mitigation factors and provides the first evidence for site-specific out of frame stop codon selection. We conjecture that the NTG universal start codon may have evolved as a consequence of TGA being a stop codon and the ability of NTGA to rapidly terminate or adjust a ribosome.

EMMA: An Extensible Mammalian Modular Assembly Toolkit for the Rapid Design and Production of Diverse Expression Vectors

Mammalian plasmid expression vectors are critical reagents underpinning many facets of research across biology, biomedical research, and the biotechnology industry. Traditional cloning methods often require laborious manual design and assembly of plasmids using tailored sequential cloning steps. This process can be protracted, complicated, expensive, and error-prone. New tools and strategies that facilitate the efficient design and production of bespoke vectors would help relieve a current bottleneck for researchers. To address this, we have developed an extensible mammalian modular assembly kit (EMMA). This enables rapid and efficient modular assembly of mammalian expression vectors in a one-tube, one-step golden-gate cloning reaction, using a standardized library of compatible genetic parts. The high modularity, flexibility, and extensibility of EMMA provide a simple method for the production of functionally diverse mammalian expression vectors. We demonstrate the value of this toolkit by constructing and validating a range of representative vectors, such as transient and stable expression vectors (transposon based vectors), targeting vectors, inducible systems, polycistronic expression cassettes, fusion proteins, and fluorescent reporters. The method also supports simple assembly combinatorial libraries and hierarchical assembly for production of larger multigenetic cargos. In summary, EMMA is compatible with automated production, and novel genetic parts can be easily incorporated, providing new opportunities for mammalian synthetic biology.

nr3c1 null mutant zebrafish are viable and reveal DNA-binding-independent activities of the glucocorticoid receptor

Glucocorticoids (GCs) play important roles in developmental and physiological processes through the transcriptional activity of their cognate receptor (Gr). Using CRISPR/Cas9 technology, we established a zebrafish null Gr mutant line and compared its phenotypes with wild type and a zebrafish line with partially silenced gr (gr^s357/s357). Homozygous gr^−/− larvae are morphologically inconspicuous and, in contrast to GR^−/− knockout mice, viable through adulthood, although with reduced fitness and early life survival. Mutants gr^−/− are fertile, but their reproductive capabilities fall at around 10 months of age, when, together with cardiac and intestinal abnormalities already visible at earlier stages, increased fat deposits are also observed. Mutants show higher levels of whole-body cortisol associated with overstimulated basal levels of crh and pomca transcripts along the HPI axis, which is unresponsive to a mechanical stressor. Transcriptional activity linked to immune response is also hampered in the gr^−/− line: after intestinal damage by dextran sodium sulphate exposure, there are neither inflammatory nor anti-inflammatory cytokine gene responses, substantiating the hypothesis of a dual-action of the GC-GR complex on the immune system. Hence, the zebrafish gr mutant line appears as a useful tool to investigate Gr functions in an integrated in vivo model.

N-terminal Proteomics Assisted Profiling of the Unexplored Translation Initiation Landscape in Arabidopsis thaliana

Proteogenomics is an emerging research field yet lacking a uniform method of analysis. Proteogenomic studies in which N-terminal proteomics and ribosome profiling are combined, suggest that a high number of protein start sites are currently missing in genome annotations. We constructed a proteogenomic pipeline specific for the analysis of N-terminal proteomics data, with the aim of discovering novel translational start sites outside annotated protein coding regions. In summary, unidentified MS/MS spectra were matched to a specific N-terminal peptide library encompassing protein N termini encoded in the Arabidopsis thaliana genome. After a stringent false discovery rate filtering, 117 protein N termini compliant with N-terminal methionine excision specificity and indicative of translation initiation were found. These include N-terminal protein extensions and translation from transposable elements and pseudogenes. Gene prediction provided supporting protein-coding models for approximately half of the protein N termini. Besides the prediction of functional domains (partially) contained within the newly predicted ORFs, further supporting evidence of translation was found in the recently released Araport11 genome re-annotation of Arabidopsis and computational translations of sequences stored in public repositories. Most interestingly, complementary evidence by ribosome profiling was found for 23 protein N termini. Finally, by analyzing protein N-terminal peptides, an in silico analysis demonstrates the applicability of our N-terminal proteogenomics strategy in revealing protein-coding potential in species with well- and poorly-annotated genomes.

Identification of long non-coding transcripts with feature selection: a comparative study

Background

The unveiling of long non-coding RNAs as important gene regulators in many biological contexts has increased the demand for efficient and robust computational methods to identify novel long non-coding RNAs from transcripts assembled with high throughput RNA-seq data. Several classes of sequence-based features have been proposed to distinguish between coding and non-coding transcripts. Among them, open reading frame, conservation scores, nucleotide arrangements, and RNA secondary structure have been used with success in literature to recognize intergenic long non-coding RNAs, a particular subclass of non-coding RNAs.

Results

In this paper we perform a systematic assessment of a wide collection of features extracted from sequence data. We use most of the features proposed in the literature, and we include, as a novel set of features, the occurrence of repeats contained in transposable elements. The aim is to detect signatures (groups of features) able to distinguish long non-coding transcripts from other classes, both protein-coding and non-coding. We evaluate different feature selection algorithms, test for signature stability, and evaluate the prediction ability of a signature with a machine learning algorithm. The study reveals different signatures in human, mouse, and zebrafish, highlighting that some features are shared among species, while others tend to be species-specific. Compared to coding potential tools and similar supervised approaches, including novel signatures, such as those identified here, in a machine learning algorithm improves the prediction performance, in terms of area under precision and recall curve, by 1 to 24%, depending on the species and on the signature.

Conclusions

Understanding which features are best suited for the prediction of long non-coding RNAs allows for the development of more effective automatic annotation pipelines especially relevant for poorly annotated genomes, such as zebrafish. We provide a web tool that recognizes novel long non-coding RNAs with the obtained signatures from fasta and gtf formats. The tool is available at the following url: http://www.bioinformatics-sannio.org/software/.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-017-1594-z) contains supplementary material, which is available to authorized users.

Genetics of Combined Pituitary Hormone Deficiency: Roadmap into the Genome Era

The genetic basis for combined pituitary hormone deficiency (CPHD) is complex, involving 30 genes in a variety of syndromic and nonsyndromic presentations. Molecular diagnosis of this disorder is valuable for predicting disease progression, avoiding unnecessary surgery, and family planning. We expect that the application of high throughput sequencing will uncover additional contributing genes and eventually become a valuable tool for molecular diagnosis. For example, in the last 3 years, six new genes have been implicated in CPHD using whole-exome sequencing. In this review, we present a historical perspective on gene discovery for CPHD and predict approaches that may facilitate future gene identification projects conducted by clinicians and basic scientists. Guidelines for systematic reporting of genetic variants and assigning causality are emerging. We apply these guidelines retrospectively to reports of the genetic basis of CPHD and summarize modes of inheritance and penetrance for each of the known genes. In recent years, there have been great improvements in databases of genetic information for diverse populations. Some issues remain that make molecular diagnosis challenging in some cases. These include the inherent genetic complexity of this disorder, technical challenges like uneven coverage, differing results from variant calling and interpretation pipelines, the number of tolerated genetic alterations, and imperfect methods for predicting pathogenicity. We discuss approaches for future research in the genetics of CPHD.

Contemporary zebrafish transgenesis with Tol2 and application for Cre/lox recombination experiments

Spatiotemporal transgene regulation by transgenic DNA recombinases is a central tool for reverse genetics in multicellular organisms, with excellent applications for misexpression and lineage tracing experiments. One of the most widespread technologies for this purpose is Cre recombinase-controlled lox site recombination that is attracting increasing interest in the zebrafish field. Tol2-mediated zebrafish transgenesis provides a stable platform to integrate lox cassette transgenes, while the amenability of the zebrafish embryo to drug treatments makes the model an ideal candidate for tamoxifen-inducible CreERT2 experiments. In addition, advanced transgenesis technologies such as phiC31 or CRISPR-Cas9-based knock-ins are even further promoting zebrafish transgenesis for Cre/lox applications. In this chapter, we will first introduce the basics of Cre/lox methodology, CreERT2 regulation by tamoxifen, as well as the utility of Tol2 and other contemporary transgenesis techniques for Cre/lox experiments. We will then outline in detail practical experimental steps for efficient transgenesis toward the creation of single-insertion transgenes and will introduce protocols for 4-hydroxytamoxifen-mediated CreERT2 induction to perform spatiotemporal lox transgene regulation experiments in zebrafish embryos. Last, we will discuss advanced experimental applications of Cre/lox beyond traditional lineage tracing approaches.

Upstream ORFs are prevalent translational repressors in vertebrates

Regulation of gene expression is fundamental in establishing cellular diversity and a target of natural selection. Untranslated mRNA regions (UTRs) are key mediators of post-transcriptional regulation. Previous studies have predicted thousands of ORFs in 5'UTRs, the vast majority of which have unknown function. Here, we present a systematic analysis of the translation and function of upstream open reading frames (uORFs) across vertebrates. Using high-resolution ribosome footprinting, we find that (i)uORFs are prevalent within vertebrate transcriptomes, (ii) the majority show signatures of active translation, and (iii)uORFs act as potent regulators of translation and RNA levels, with a similar magnitude to miRNAs. Reporter experiments reveal clear repression of downstream translation by uORFs/oORFs. uORF number, intercistronic distance, overlap with the CDS, and initiation context most strongly influence translation. Evolution has targeted these features to favor uORFs amenable to regulation over constitutively repressive uORFs/oORFs. Finally, we observe that the regulatory potential of uORFs on individual genes is conserved across species. These results provide insight into the regulatory code within mRNA leader sequences and their capacity to modulate translation across vertebrates.

A Comprehensive Analysis of Codon Usage Patterns in Blunt Snout Bream (Megalobrama amblycephala) Based on RNA-Seq Data

Blunt snout bream (Megalobrama amblycephala) is an important fish species for its delicacy and high economic value in China. Codon usage analysis could be helpful to understand its codon biology, mRNA translation and vertebrate evolution. Based on RNA-Seq data for M. amblycephala, high-frequency codons (CUG, AGA, GUG, CAG and GAG), as well as low-frequency ones (NUA and NCG codons) were identified. A total of 724 high-frequency codon pairs were observed. Meanwhile, 14 preferred and 199 avoided neighboring codon pairs were also identified, but bias was almost not shown with one or more intervening codons inserted between the same pairs. Codon usage bias in the regions close to start and stop codons indicated apparent heterogeneity, which even occurs in the flanking nucleotide sequence. Codon usage bias (RSCU and SCUO) was related to GC₃ (GC content of 3rd nucleotide in codon) bias. Six GO (Gene ontology) categories and the number of methylation targets were influenced by GC₃. Codon usage patterns comparison among 23 vertebrates showed species specificities by using GC contents, codon usage and codon context analysis. This work provided new insights into fish biology and new information for breeding projects.