Discovering the hidden function in fungal genomes.

New molecular technologies have helped unveil previously unexplored facets of the genome beyond the canonical proteome, including microproteins and short ORFs, products of alternative splicing, regulatory non-coding RNAs, as well as transposable elements, cis-regulatory DNA, and other highly repetitive regions of DNA. In this Review, we highlight what is known about this 'hidden genome' within the fungal kingdom. Using well-established model systems as a contextual framework, we describe key elements of this hidden genome in diverse fungal species, and explore how these factors perform critical functions in regulating fungal metabolism, stress tolerance, and pathogenesis. Finally, we discuss new technologies that may be adapted to further characterize the hidden genome in fungi.

Epidemiology and molecular evolution of GI.1 sapovirus in the recent era.

Sapovirus (SaV) infection is increasing worldwide. Herein, we provided evidence of a significant increase in SaV infection in Japan during 2010-2022, primarily due to the considerable (p = 0.0003) rise of the GI.1 genotype. Furthermore, we found that all major and minor SaV outbreaks in Japan, including the largest SaV outbreak in 2021-2022, were caused by the GI.1 genotype. Therefore, to get insight into the underlying molecular mechanism behind this rising trend of the SaV GI.1 type, we selected 15 SaV GI.1 outbreak strains for complete genome analysis through next-generation sequencing. Phylogenetically, our strains remained clustered in different branches in lineages I and II among the GI.1 genotype. We showed all amino acid (aa) substitutions in different open reading frames (ORFs) in these strains. Importantly, we have demonstrated that the strains involved in the largest SaV outbreak in Japan in 2021-2022 belonged to lineage II and possessed the third ORF. We have identified some unique aa mutations in these major outbreak strains in the NS1 and NS6-NS7 regions that are thought to be associated with viral pathogenicity, cell tropism, and epidemiological competence. Thus, in addition to enriching the database of SaV's complete sequences, this study provides insights into its important mutations.

Comparative Analysis of the Mitochondrial Genome of Eggplant (Solanum melongena L.) to Identify Cytoplasmic Male Sterility Candidate Genes.

Cytoplasmic male sterility (CMS) is important for commercial hybrid seed production. However, it is still not used in eggplant (Solanum melongena L.), and corresponding regulatory genes and mechanisms of action have not been reported. We report CMS line 327A, which was derived from the hybridization between cultivated and wild eggplants. By looking at different stages of anther development under a microscope, we saw that the 327A anther's tapetum layer vacuolized during meiosis, which caused abortion. To investigate the 327A CMS regulatory genes, the mitochondrial genomes of 327A and its maintainer line 327B were assembled de novo. It was found that 15 unique ORFs (Open Reading Frame) were identified in 327A. RT-PCR and RT-QPCAR tests confirmed that orf312a and orf172a, 327A-specific ORFs with a transmembrane domain, were strongly expressed in sterile anthers of 327A. In addition, orf312a has a chimeric structure with the ribosomal protein subunit rpl16. Therefore, orf312a and orf172a can be considered strong candidate genes for CMS. Concurrently, we analyzed the characteristics of CMS to develop a functional molecular marker, CMS312, targeting a future theoretical basis for eggplant CMS three-line molecular breeding.

Polyomavirus ALTOs, but not MTs, downregulate viral early gene expression by activating the NF-κB pathway.

Polyomaviruses are small, circular dsDNA viruses that can cause cancer. Alternative splicing of polyomavirus early transcripts generates large and small tumor antigens (LT, ST) that play essential roles in viral replication and tumorigenesis. Some polyomaviruses also express middle tumor antigens (MTs) or alternate LT open reading frames (ALTOs), which are evolutionarily related but have distinct gene structures. MTs are a splice variant of the early transcript whereas ALTOs are overprinted on the second exon of the LT transcript in an alternate reading frame and are translated via an alternative start codon. Merkel cell polyomavirus (MCPyV), the only human polyomavirus that causes cancer, encodes an ALTO but its role in the viral lifecycle and tumorigenesis has remained elusive. Here, we show MCPyV ALTO acts as a tumor suppressor and is silenced in Merkel cell carcinoma (MCC). Rescuing ALTO in MCC cells induces growth arrest and activates NF-κB signaling. ALTO activates NF-κB by binding SQSTM1 and TRAF2&3 via two N-Terminal Activating Regions (NTAR1+2), resembling Epstein-Barr virus (EBV) Latent Membrane Protein 1 (LMP1). Following activation, NF-κB dimers bind the MCPyV noncoding control region (NCCR) and downregulate early transcription. Beyond MCPyV, NTAR motifs are conserved in other polyomavirus ALTOs, which activate NF-κB signaling, but are lacking in MTs that do not. Furthermore, polyomavirus ALTOs downregulate their respective viral early transcription in an NF-κB- and NTAR-dependent manner. Our findings suggest that ALTOs evolved to suppress viral replication and promote viral latency and that MCPyV ALTO must be silenced for MCC to develop.

Rp3: Ribosome profiling-assisted proteogenomics improves coverage and confidence during microprotein discovery.

There has been a dramatic increase in the identification of non-canonical translation and a significant expansion of the protein-coding genome. Among the strategies used to identify unannotated small Open Reading Frames (smORFs) that encode microproteins, Ribosome profiling (Ribo-Seq) is the gold standard for the annotation of novel coding sequences by reporting on smORF translation. In Ribo-Seq, ribosome-protected footprints (RPFs) that map to multiple genomic sites are removed since they cannot be unambiguously assigned to a specific genomic location. Furthermore, RPFs necessarily result in short (25-34 nucleotides) reads, increasing the chance of multi-mapping alignments, such that smORFs residing in these regions cannot be identified by Ribo-Seq. Moreover, it has been challenging to identify protein evidence for Ribo-Seq. To solve this, we developed Rp3, a pipeline that integrates proteogenomics and Ribosome profiling to provide unambiguous evidence for a subset of microproteins missed by current Ribo-Seq pipelines. Here, we show that Rp3 maximizes proteomics detection and confidence of microprotein-encoding smORFs.

Ribosome Profiling and RNA Sequencing Reveal Translation and Transcription Regulation under Acute Heat Stress in Rainbow Trout (Oncorhynchus mykiss, Walbaum, 1792) Liver.

Rainbow trout (Oncorhynchus mykiss, Walbaum, 1792) is an important economic cold-water fish that is susceptible to heat stress. To date, the heat stress response in rainbow trout is more widely understood at the transcriptional level, while little research has been conducted at the translational level. To reveal the translational regulation of heat stress in rainbow trout, in this study, we performed a ribosome profiling assay of rainbow trout liver under normal and heat stress conditions. Comparative analysis of the RNA-seq data with the ribosome profiling data showed that the folding changes in gene expression at the transcriptional level are moderately correlated with those at the translational level. In total, 1213 genes were significantly altered at the translational level. However, only 32.8% of the genes were common between both levels, demonstrating that heat stress is coordinated across both transcriptional and translational levels. Moreover, 809 genes exhibited significant differences in translational efficiency (TE), with the TE of these genes being considerably affected by factors such as the GC content, coding sequence length, and upstream open reading frame (uORF) presence. In addition, 3468 potential uORFs in 2676 genes were identified, which can potentially affect the TE of the main open reading frames. In this study, Ribo-seq and RNA-seq were used for the first time to elucidate the coordinated regulation of transcription and translation in rainbow trout under heat stress. These findings are expected to contribute novel data and theoretical insights to the international literature on the thermal stress response in fish.

Assembly and Comparative Analysis of the Complete Mitochondrial Genome of Saussurea inversa (Asteraceae).

Saussurea inversa is a perennial herb used in traditional Chinese medicine and is effective against rheumatoid arthritis. In this study, we sequenced the complete mitochondrial (mt) genome of S. inversa (GenBank accession number: ON584565.1). The circular mt genome of S. inversa was 335,372 bp in length, containing 62 genes, including 33 mRNAs, 22 tRNAs, 6 rRNAs, and 1 pseudogene, along with 1626 open reading frames. The GC content was 45.14%. Predictive analysis revealed substantial RNA editing, with ccmFn being the most abundantly edited gene, showing 36 sites. Gene migration between the mt and chloroplast (cp) genomes of S. inversa was observed through the detection of homologous gene fragments. Phylogenetic analysis revealed that S. inversa was clustered with Arctium tomentosum (Asteraceae). Our findings provide extensive information regarding the mt genome of S. inversa and help lay the foundation for future studies on its genetic variations, phylogeny, and breeding via the analysis of the mt genome.

CRISPR/Cas9-mediated disruption of orf76 as an antiviral therapy against BmNPV in the transgenic silkworm.

Viral diseases pose a significant threat to livestock husbandry and plant cultivation. CRISPR/Cas9-mediated targeted editing of viral genes offers a promising approach to antiviral therapy. The silkworm, Bombyx mori, is an economically important insect susceptible to infection by B. mori nucleopolyhedrovirus (BmNPV), and viral outbreaks cause severe economic losses to the sericulture industry. Here, we identified BmNPV orf76 as a viral late gene that is highly similar to Autographa californica multiple nucleopolyhedrovirus Ac93. The deletion of orf76 abolished BmNPV proliferation and hindered the production of infectious budded viruses. We generated a transgenic line, Cas9(+)/sgorf76(+), that did not affect the growth or development of the silkworm and demonstrated that the transgenic line Cas9(+)/sgorf76(+) efficiently cleaved orf76 at the sgorf76 site, resulting in large deletions at 120 h post-infection, with no observed off-target effects. Survival analyses revealed that the transgenic line Cas9(+)/sgorf76(+) exhibited significantly higher survival rates than the control lines Cas9(-)/sgorf76(-), regardless of the BmNPV inoculation dose. Additionally, the number of BmNPV DNA copies and the expression levels of viral genes were markedly inhibited in the transgenic line Cas9(+)/sgorf76(+) compared with the control line Cas9(-)/sgorf76(-). The results provide a promising target for Cas9-mediated antiviral therapy against BmNPV, and the findings provide new insights for baculovirus gene function studies and lepidopteran pest control.

NMDtxDB: data-driven identification and annotation of human NMD target transcripts.

The nonsense-mediated RNA decay (NMD) pathway is a crucial mechanism of mRNA quality control. Current annotations of NMD substrate RNAs are rarely data-driven, but use generally established rules. We present a data set with four cell lines and combinations for SMG5, SMG6, and SMG7 knockdowns or SMG7 knockout. Based on this data set, we implemented a workflow that combines Nanopore and Illumina sequencing to assemble a transcriptome, which is enriched for NMD target transcripts. Moreover, we use coding sequence information (CDS) from Ensembl, Gencode consensus Ribo-seq ORFs, and OpenProt to enhance the CDS annotation of novel transcript isoforms. In summary, 302,889 transcripts were obtained from the transcriptome assembly process, out of which 24% are absent from Ensembl database annotations, 48,213 contain a premature stop codon, and 6433 are significantly upregulated in three or more comparisons of NMD active versus deficient cell lines. We present an in-depth view of these results through the NMDtxDB database, which is available at https://shiny.dieterichlab.org/app/NMDtxDB, and supports the study of NMD-sensitive transcripts. We open sourced our implementation of the respective web-application and analysis workflow at https://github.com/dieterich-lab/NMDtxDB and https://github.com/dieterich-lab/nmd-wf.

A catalog of small proteins from the global microbiome.

Small open reading frames (smORFs) shorter than 100 codons are widespread and perform essential roles in microorganisms, where they encode proteins active in several cell functions, including signal pathways, stress response, and antibacterial activities. However, the ecology, distribution and role of small proteins in the global microbiome remain unknown. Here, we construct a global microbial smORFs catalog (GMSC) derived from 63,410 publicly available metagenomes across 75 distinct habitats and 87,920 high-quality isolate genomes. GMSC contains 965 million non-redundant smORFs with comprehensive annotations. We find that archaea harbor more smORFs proportionally than bacteria. We moreover provide a tool called GMSC-mapper to identify and annotate small proteins from microbial (meta)genomes. Overall, this publicly-available resource demonstrates the immense and underexplored diversity of small proteins.

Elucidating the Subcellular Localization of GLRaV-3 Proteins Encoded by the Unique Gene Block in N. benthamiana Suggests Implications on Plant Host Suppression.

Grapevine leafroll-associated virus 3 (GLRaV-3) is a formidable threat to the stability of the global grape and wine industries. It is the primary etiological agent of grapevine leafroll disease (GLD) and significantly impairs vine health, fruit quality, and yield. GLRaV-3 is a member of the genus Ampelovirus, Closteroviridae family. Viral genes within the 3' proximal unique gene blocks (UGB) remain highly variable and poorly understood. The UGBs of Closteroviridae viruses include diverse open reading frames (ORFs) that have been shown to contribute to viral functions such as the suppression of the host RNA silencing defense response and systemic viral spread. This study investigates the role of GLRaV-3 ORF8, ORF9, and ORF10, which encode the proteins p21, p20A, and p20B, respectively. These genes represent largely unexplored facets of the GLRaV-3 genome. Here, we visualize the subcellular localization of wildtype and mutagenized GLRaV-3 ORFs 8, 9, and 10, transiently expressed in Nicotiana benthamiana. Our results indicate that p21 localizes to the cytosol, p20A associates with microtubules, and p20B is trafficked into the nucleus to carry out the suppression of host RNA silencing. The findings presented herein provide a foundation for future research aimed at the characterization of the functions of these ORFs. In the long run, it would also facilitate the development of innovative strategies to understand GLRaV-3, mitigate its spread, and impacts on grapevines and the global wine industry.

Comprehensive analysis of the functional impact of single nucleotide variants of human CHEK2.

Loss of function mutations in the checkpoint kinase gene CHEK2 are associated with increased risk of breast and other cancers. Most of the 3,188 unique amino acid changes that can result from non-synonymous single nucleotide variants (SNVs) of CHEK2, however, have not been tested for their impact on the function of the CHEK2-enocded protein (CHK2). One successful approach to testing the function of variants has been to test for their ability to complement mutations in the yeast ortholog of CHEK2, RAD53. This approach has been used to provide functional information on over 100 CHEK2 SNVs and the results align with functional assays in human cells and known pathogenicity. Here we tested all but two of the 4,887 possible SNVs in the CHEK2 open reading frame for their ability to complement RAD53 mutants using a high throughput technique of deep mutational scanning (DMS). Among the non-synonymous changes, 770 were damaging to protein function while 2,417 were tolerated. The results correlate well with previous structure and function data and provide a first or additional functional assay for all the variants of uncertain significance identified in clinical databases. Combined, this approach can be used to help predict the pathogenicity of CHEK2 variants of uncertain significance that are found in susceptibility screening and could be applied to other cancer risk genes.

Regulation viral RNA transcription and replication by higher-order RNA structures within the nsp1 coding region of MERS coronavirus.

Coronavirus (CoV) possesses numerous functional cis-acting elements in its positive-strand genomic RNA. Although most of these RNA structures participate in viral replication, the functions of RNA structures in the genomic RNA of CoV in viral replication remain unclear. In this study, we investigated the functions of the higher-order RNA stem-loop (SL) structures SL5B, SL5C, and SL5D in the ORF1a coding region of Middle East respiratory syndrome coronavirus (MERS-CoV) in viral replication. Our approach, using reverse genetics of a bacterial artificial chromosome system, revealed that SL5B and SL5C play essential roles in the discontinuous transcription of MERS-CoV. In silico analyses predicted that SL5C interacts with a bulged stem-loop (BSL) in the 3' untranslated region, suggesting that the RNA structure of SL5C is important for viral RNA transcription. Conversely, SL5D did not affect transcription, but mediated the synthesis of positive-strand genomic RNA. Additionally, the RNA secondary structure of SL5 in the revertant virus of the SL5D mutant was similar to that of the wild-type, indicating that the RNA structure of SL5D can finely tune RNA replication in MERS-CoV. Our data indicate novel regulatory mechanisms of viral RNA transcription and replication by higher-order RNA structures in the MERS-CoV genomic RNA.

Coding, or non-coding, that is the question.

The advent of high-throughput sequencing uncovered that our genome is pervasively transcribed into RNAs that are seemingly not translated into proteins. It was also found that non-coding RNA transcripts outnumber canonical protein-coding genes. This mindboggling discovery prompted a surge in non-coding RNA research that started unraveling the functional relevance of these new genetic units, shaking the classic definition of "gene". While the non-coding RNA revolution was still taking place, polysome/ribosome profiling and mass spectrometry analyses revealed that peptides can be translated from non-canonical open reading frames. Therefore, it is becoming evident that the coding vs non-coding dichotomy is way blurrier than anticipated. In this review, we focus on several examples in which the binary classification of coding vs non-coding genes is outdated, since the same bifunctional gene expresses both coding and non-coding products. We discuss the implications of this intricate usage of transcripts in terms of molecular mechanisms of gene expression and biological outputs, which are often concordant, but can also surprisingly be discordant. Finally, we discuss the methodological caveats that are associated with the study of bifunctional genes, and we highlight the opportunities and challenges of therapeutic exploitation of this intricacy towards the development of anticancer therapies.

Molecular Cloning, Characterization, and Function of Insulin-Related Peptide 1 (IRP1) in the Haliotis discus hanna.

Abalone is a popular mollusk in the marine aquaculture industry of China. However, existing challenges, like slow growth, individual miniaturization, and the absence of abundant abalone, have emerged as significant obstacles impeding its long-term progress in aquaculture. Studies have demonstrated that insulin-related peptide (IRP) is a crucial factor in the growth of marine organisms. However, limited studies have been conducted on IRP in abalone. This study indicated that the hdh-MIRP1 open reading frame (ORF) was composed of 456 base pairs, which encoded 151 amino acids. Based on the gene expression and immunofluorescence analyses, the cerebral ganglion of Haliotis discus hannai (H. discus hannai) was the primary site of hdh-MIRP1 mRNA expression. Moreover, hdh-MIRP1 expression was observed to be higher in the larger group than in the smaller group abalones. Only single nucleotide polymorphism (SNP) was related to their growth characteristics. However, approximately 82 proteins that may interact with hdh-MIRP1 were identified. The functional enrichment analysis of the 82 genes indicated that hdh-MIRP1 may be involved in the regulation of glucose metabolism and the process of growth. This study established a benchwork for further investigating the role of IRP in the growth of abalone.

Genetic variations and gene expression profiles of Rice Black-streaked dwarf virus (RBSDV) in different host plants and insect vectors: insights from RNA-Seq analysis.

Rice black-streaked dwarf virus (RBSDV) is an etiological agent of a destructive disease infecting some economically important crops from the Gramineae family in Asia. While RBSDV causes high yield losses, genetic characteristics of replicative viral populations have not been investigated within different host plants and insect vectors. Herein, eleven publicly available RNA-Seq datasets from Chinese RBSDV-infected rice, maize, and viruliferous planthopper (Laodelphax striatellus) were obtained from the NCBI database. The patterns of SNP and RNA expression profiles of expected RBSDV populations were analyzed by CLC Workbench 20 and Geneious Prime software. These analyses discovered 2,646 mutations with codon changes in RBSDV whole transcriptome and forty-seven co-mutated hotspots with high variant frequency within the crucial regions of S5-1, S5-2, S6, S7-1, S7-2, S9, and S10 open reading frames (ORFs) which are responsible for some virulence and host range functions. Moreover, three joint mutations are located on the three-dimensional protein of P9-1. The infected RBSDV-susceptible rice cultivar KTWYJ3 and indigenous planthopper datasets showed more co-mutated hotspot numbers than others. Our analyses showed the expression patterns of viral genomic fragments varied depending on the host type. Unlike planthopper, S5-1, S2, S6, and S9-1 ORFs, respectively had the greatest read numbers in host plants; and S5-2, S9-2, and S7-2 were expressed in the lowest level. These findings underscore virus/host complexes are effective in the genetic variations and gene expression profiles of plant viruses. Our analysis revealed no evidence of recombination events. Interestingly, the negative selection was observed at 12 RBSDV ORFs, except for position 1015 in the P1 protein, where a positive selection was detected. The research highlights the potential of SRA datasets for analysis of the virus cycle and enhances our understanding of RBSDV's genetic diversity and host specificity.

SEMPER: Stoichiometric expression of mRNA polycistrons by eukaryotic ribosomes for compact, ratio-tunable multi-gene expression.

Here, we present a method for expressing multiple open reading frames (ORFs) from single transcripts using the leaky scanning model of translation initiation. In this approach termed "stoichiometric expression of mRNA polycistrons by eukaryotic ribosomes" (SEMPER), adjacent ORFs are translated from a single mRNA at tunable ratios determined by their order in the sequence and the strength of their translation initiation sites. We validate this approach by expressing up to three fluorescent proteins from one plasmid in two different cell lines. We then use it to encode a stoichiometrically tuned polycistronic construct encoding gas vesicle acoustic reporter genes that enables efficient formation of the multi-protein complex while minimizing cellular toxicity. We also demonstrate that SEMPER enables polycistronic expression of recombinant monoclonal antibodies from plasmid DNA and of two fluorescent proteins from single mRNAs made through in vitro transcription. Finally, we provide a probabilistic model to elucidate the mechanisms underlying SEMPER. A record of this paper's transparent peer review process is included in the supplemental information.

How antisense transcripts can evolve to encode novel proteins.

Protein coding features can emerge de novo in non coding transcripts, resulting in emergence of new protein coding genes. Studies across many species show that a large fraction of evolutionarily novel non-coding RNAs have an antisense overlap with protein coding genes. The open reading frames (ORFs) in these antisense RNAs could also overlap with existing ORFs. In this study, we investigate how the evolution an ORF could be constrained by its overlap with an existing ORF in three different reading frames. Using a combination of mathematical modeling and genome/transcriptome data analysis in two different model organisms, we show that antisense overlap can increase the likelihood of ORF emergence and reduce the likelihood of ORF loss, especially in one of the three reading frames. In addition to rationalising the repeatedly reported prevalence of de novo emerged genes in antisense transcripts, our work also provides a generic modeling and an analytical framework that can be used to understand evolution of antisense genes.

The signal peptide of yeast killer toxin K2 confers producer self-protection and allows conversion into a modular toxin-immunity system.

Some microbial toxins also target the producer species itself, necessitating a means of self-protection. The M2 double-stranded RNA (dsRNA) killer virus in Saccharomyces cerevisiae contains a single open reading frame (ORF) encoding both the secreted pore-forming toxin K2 as well as a cognate immunity factor. Here, we show that expression of a 49-amino acid N-terminal peptide from the K2 precursor is both necessary and sufficient for immunity. This immunity peptide simultaneously functions as a signal peptide for toxin secretion and protects the cell against the cytotoxic K2 α subunit. The K2 toxin and immunity factor can be functionally separated into two ORFs, yielding a modular toxin-immunity system. This case further shows how a (signal) peptide can carry the potential for providing cellular protection against an antimicrobial toxin.

Evidence for widespread translation of 5' untranslated regions.

Ribosome profiling experiments support the translation of a range of novel human open reading frames. By contrast, most peptides from large-scale proteomics experiments derive from just one source, 5' untranslated regions. Across the human genome we find evidence for 192 translated upstream regions, most of which would produce protein isoforms with extended N-terminal ends. Almost all of these N-terminal extensions are from highly abundant genes, which suggests that the novel regions we detect are just the tip of the iceberg. These upstream regions have characteristics that are not typical of coding exons. Their GC-content is remarkably high, even higher than 5' regions in other genes, and a large majority have non-canonical start codons. Although some novel upstream regions have cross-species conservation - five have orthologues in invertebrates for example - the reading frames of two thirds are not conserved beyond simians. These non-conserved regions also have no evidence of purifying selection, which suggests that much of this translation is not functional. In addition, non-conserved upstream regions have significantly more peptides in cancer cell lines than would be expected, a strong indication that an aberrant or noisy translation initiation process may play an important role in translation from upstream regions.