Teaching bioinformatics through the analysis of SARS-CoV-2: project-based training for computer science students.

MOTIVATION: We learn more effectively through experience and reflection than through passive reception of information. Bioinformatics offers an excellent opportunity for project-based learning. Molecular data are abundant and accessible in open repositories, and important concepts in biology can be rediscovered by reanalyzing the data. RESULTS: In the manuscript, we report on five hands-on assignments we designed for master's computer science students to train them in bioinformatics for genomics. These assignments are the cornerstones of our introductory bioinformatics course and are centered around the study of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). They assume no prior knowledge of molecular biology but do require programming skills. Through these assignments, students learn about genomes and genes, discover their composition and function, relate SARS-CoV-2 to other viruses, and learn about the body's response to infection. Student evaluation of the assignments confirms their usefulness and value, their appropriate mastery-level difficulty, and their interesting and motivating storyline. AVAILABILITY AND IMPLEMENTATION: The course materials are freely available on GitHub at https://github.com/IB-ULFRI.

SIMBA Method-Simultaneous Detection of Antimicrobial and Anti-biofilm Activity of New Compounds Using Salmonella Infantis.

The development of antimicrobial resistance and the formation of Salmonella biofilms are serious public health problems. For this reason, new natural compounds with antimicrobial and anti-biofilm activity are being sought, and wild fungi represent an untapped potential. Various extraction agents, including organic solvents and aqueous buffers, can be used to obtain bioactive compounds from natural sources. To evaluate their bioactivity, extensive screening studies are required to determine antimicrobial and anti-biofilm activity using methods such as broth microdilution or crystal violet assay, respectively, but none of these methods allow simultaneous evaluation of both activities against bacteria. Cold water extraction from wild fungi offers the advantage of extracting water-soluble compounds. The SIMultaneous detection of antiMicrobial and anti-Biofilm Activity (SIMBA) method combines the testing of both types of activity against bacteria with the evaluation of the 20 h growth curve of the Salmonella Infantis ZM9 strain determined with absorbance measurements at 600 nm in a 96-well plate. SIMBA method thus shortens the time to determine the bioactivity of extracts, reduces material consumption, and eliminates the need for additional reagents. SIMBA enables rapid selection of bioactive extracts for their fractionation and shortens the time to determine new natural products with antimicrobial and anti-biofilm activity. Graphical overview.

Identification of novel genes involved in neutral lipid storage by quantitative trait loci analysis of Saccharomyces cerevisiae.

BACKGROUND: The accumulation of intracellular fat depots is a polygenic trait. Therefore, the extent of lipid storage in the individuals of a species covers a broad range and is determined by many genetic factors. Quantitative trait loci analysis can be used to identify those genetic differences between two strains of the same species that are responsible for the differences in a given phenotype. We used this method and complementary approaches to identify genes in the yeast Saccharomyces cerevisiae that are involved in neutral lipid storage. RESULTS: We selected two yeast strains, the laboratory strain BY4741 and the wine yeast AWRI1631, with a more than two-fold difference in neutral lipid content. After crossing, sporulation and germination, we used fluorescence activated cell sorting to isolate a subpopulation of cells with the highest neutral lipid content from the pool of segregants. Whole genome sequencing of this subpopulation and of the unsorted pool of segregants implicated several loci that are involved in lipid accumulation. Three of the identified genes, PIG1, PHO23 and RML2, were investigated in more detail. Deletions of these genes and the exchange of the alleles between the two parental strains confirmed that the encoded proteins contribute to neutral lipid storage in S. cerevisiae and that PIG1, PHO23 and RML2 are the major causative genes. Backcrossing of one of the segregants with the parental strains for seven generations revealed additional regions in the genomes of both strains with potential causative genes for the high lipid accumulation phenotype. CONCLUSIONS: We identified several genes that contribute to the phenotype of lipid accumulation in an allele-specific manner. Surprisingly, no allelic variations of genes with known functions in lipid metabolism were found, indicating that the level of storage lipid accumulation is determined by many cellular processes that are not directly related to lipid metabolism.

Sparse data embedding and prediction by tropical matrix factorization.

BACKGROUND: Matrix factorization methods are linear models, with limited capability to model complex relations. In our work, we use tropical semiring to introduce non-linearity into matrix factorization models. We propose a method called Sparse Tropical Matrix Factorization (STMF) for the estimation of missing (unknown) values in sparse data. RESULTS: We evaluate the efficiency of the STMF method on both synthetic data and biological data in the form of gene expression measurements downloaded from The Cancer Genome Atlas (TCGA) database. Tests on unique synthetic data showed that STMF approximation achieves a higher correlation than non-negative matrix factorization (NMF), which is unable to recover patterns effectively. On real data, STMF outperforms NMF on six out of nine gene expression datasets. While NMF assumes normal distribution and tends toward the mean value, STMF can better fit to extreme values and distributions. CONCLUSION: STMF is the first work that uses tropical semiring on sparse data. We show that in certain cases semirings are useful because they consider the structure, which is different and simpler to understand than it is with standard linear algebra.

Author Correction: A systems view of spliceosomal assembly and branchpoints with iCLIP.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

A systems view of spliceosomal assembly and branchpoints with iCLIP.

Studies of spliceosomal interactions are challenging due to their dynamic nature. Here we used spliceosome iCLIP, which immunoprecipitates SmB along with small nuclear ribonucleoprotein particles and auxiliary RNA binding proteins, to map spliceosome engagement with pre-messenger RNAs in human cell lines. This revealed seven peaks of spliceosomal crosslinking around branchpoints (BPs) and splice sites. We identified RNA binding proteins that crosslink to each peak, including known and candidate splicing factors. Moreover, we detected the use of over 40,000 BPs with strong sequence consensus and structural accessibility, which align well to nearby crosslinking peaks. We show how the position and strength of BPs affect the crosslinking patterns of spliceosomal factors, which bind more efficiently upstream of strong or proximally located BPs and downstream of weak or distally located BPs. These insights exemplify spliceosome iCLIP as a broadly applicable method for transcriptomic studies of splicing mechanisms.

Metabolomic profiling of CHO fed-batch growth phases at 10, 100, and 1,000 L.

Established bioprocess monitoring is based on quick and reliable methods, including cell count and viability measurement, extracellular metabolite measurement, and the measurement of physicochemical qualities of the cultivation medium. These methods are sufficient for monitoring of process performance, but rarely give insight into the actual physiological states of the cell culture. However, understanding of the latter is essential for optimization of bioprocess development. Our study used LC-MS metabolomics as a tool for additional resolution of bioprocess monitoring and was designed at three bioreactors scales (10 L, 100 L, and 1,000 L) to gain insight into the basal metabolic states of the Chinese hamster ovary (CHO) cell culture during fed-batch. Metabolites characteristics of the four growth stages (early and late exponential phase, stationary phase, and the phase of decline) were identified by multivariate analysis. Enriched metabolic pathways were then established for each growth phase using the CHO metabolic network model. Biomass generation and nucleotide synthesis were enriched in early exponential phase, followed by increased protein production and imbalanced glutathione metabolism in late exponential phase. Glycolysis became downregulated in stationary phase and amino-acid metabolism increased. Phase of culture decline resulted in rise of oxidized glutathione and fatty acid concentrations. Intracellular metabolic profiles of the CHO fed-batch culture were also shown to be consistent with scale and thus demonstrate metabolomic profiling as an informative method to gain physiological insight into the cell culture states during bioprocess regardless of scale.

An ancient germ cell-specific RNA-binding protein protects the germline from cryptic splice site poisoning.

Male germ cells of all placental mammals express an ancient nuclear RNA binding protein of unknown function called RBMXL2. Here we find that deletion of the retrogene encoding RBMXL2 blocks spermatogenesis. Transcriptome analyses of age-matched deletion mice show that RBMXL2 controls splicing patterns during meiosis. In particular, RBMXL2 represses the selection of aberrant splice sites and the insertion of cryptic and premature terminal exons. Our data suggest a Rbmxl2 retrogene has been conserved across mammals as part of a splicing control mechanism that is fundamentally important to germ cell biology. We propose that this mechanism is essential to meiosis because it buffers the high ambient concentrations of splicing activators, thereby preventing poisoning of key transcripts and disruption to gene expression by aberrant splice site selection.

Differential expression of microRNAs and other small RNAs in muscle tissue of patients with ALS and healthy age-matched controls.

Amyotrophic lateral sclerosis is a late-onset disorder primarily affecting motor neurons and leading to progressive and lethal skeletal muscle atrophy. Small RNAs, including microRNAs (miRNAs), can serve as important regulators of gene expression and can act both globally and in a tissue-/cell-type-specific manner. In muscle, miRNAs called myomiRs govern important processes and are deregulated in various disorders. Several myomiRs have shown promise for therapeutic use in cellular and animal models of ALS; however, the exact miRNA species differentially expressed in muscle tissue of ALS patients remain unknown. Following small RNA-Seq, we compared the expression of small RNAs in muscle tissue of ALS patients and healthy age-matched controls. The identified snoRNAs, mtRNAs and other small RNAs provide possible molecular links between insulin signaling and ALS. Furthermore, the identified miRNAs are predicted to target proteins that are involved in both normal processes and various muscle disorders and indicate muscle tissue is undergoing active reinnervation/compensatory attempts thus providing targets for further research and therapy development in ALS.

Genome Sequence of a Lethal Strain of Xylem-Invading Verticillium nonalfalfae.

Verticillium nonalfalfae, a soilborne vascular phytopathogenic fungus, causes wilt disease in several crop species. Of great concern are outbreaks of highly aggressive V. nonalfalfae strains, which cause a devastating wilt disease in European hops. We report here the genome sequence and annotation of V. nonalfalfae strain T2, providing genomic information that will allow better understanding of the molecular mechanisms underlying the development of highly aggressive strains.

Identification of RNA-binding domains of RNA-binding proteins in cultured cells on a system-wide scale with RBDmap.

This protocol is an extension to: Nat. Protoc. 8, 491-500 (2013); doi:10.1038/nprot.2013.020; published online 14 February 2013RBDmap is a method for identifying, in a proteome-wide manner, the regions of RNA-binding proteins (RBPs) engaged in native interactions with RNA. In brief, cells are irradiated with UV light to induce protein-RNA cross-links. Following stringent denaturing washes, the resulting covalently linked protein-RNA complexes are purified with oligo(dT) magnetic beads. After elution, RBPs are subjected to partial proteolysis, in which the protein regions still bound to the RNA and those released to the supernatant are separated by a second oligo(dT) selection. After sample preparation and mass-spectrometric analysis, peptide intensity ratios between the RNA-bound and released fractions are used to determine the RNA-binding regions. As a Protocol Extension, this article describes an adaptation of an existing Protocol and offers additional applications. The earlier protocol (for the RNA interactome capture method) describes how to identify the active RBPs in cultured cells, whereas this Protocol Extension also enables the identification of the RNA-binding domains of RBPs. The experimental workflow takes 1 week plus 2 additional weeks for proteomics and data analysis. Notably, RBDmap presents numerous advantages over classic methods for determining RNA-binding domains: it produces proteome-wide, high-resolution maps of the protein regions contacting the RNA in a physiological context and can be adapted to different biological systems and conditions. Because RBDmap relies on the isolation of polyadenylated RNA via oligo(dT), it will not provide RNA-binding information on proteins interacting exclusively with nonpolyadenylated transcripts. Applied to HeLa cells, RBDmap uncovered 1,174 RNA-binding sites in 529 proteins, many of which were previously unknown.

Tia1 dependent regulation of mRNA subcellular location and translation controls p53 expression in B cells.

Post-transcriptional regulation of cellular mRNA is essential for protein synthesis. Here we describe the importance of mRNA translational repression and mRNA subcellular location for protein expression during B lymphocyte activation and the DNA damage response. Cytoplasmic RNA granules are formed upon cell activation with mitogens, including stress granules that contain the RNA binding protein Tia1. Tia1 binds to a subset of transcripts involved in cell stress, including p53 mRNA, and controls translational silencing and RNA granule localization. DNA damage promotes mRNA relocation and translation in part due to dissociation of Tia1 from its mRNA targets. Upon DNA damage, p53 mRNA is released from stress granules and associates with polyribosomes to increase protein synthesis in a CAP-independent manner. Global analysis of cellular mRNA abundance and translation indicates that this is an extended ATM-dependent mechanism to increase protein expression of key modulators of the DNA damage response.Sequestering mRNA in cytoplasmic stress granules is a mechanism for translational repression. Here the authors find that p53 mRNA, present in stress granules in activated B lymphocytes, is released upon DNA damage and is translated in a CAP-independent manner.

High-Resolution RNA Maps Suggest Common Principles of Splicing and Polyadenylation Regulation by TDP-43.

Many RNA-binding proteins (RBPs) regulate both alternative exons and poly(A) site selection. To understand their regulatory principles, we developed expressRNA, a web platform encompassing computational tools for integration of iCLIP and RNA motif analyses with RNA-seq and 3' mRNA sequencing. This reveals at nucleotide resolution the "RNA maps" describing how the RNA binding positions of RBPs relate to their regulatory functions. We use this approach to examine how TDP-43, an RBP involved in several neurodegenerative diseases, binds around its regulated poly(A) sites. Binding close to the poly(A) site generally represses, whereas binding further downstream enhances use of the site, which is similar to TDP-43 binding around regulated exons. Our RNAmotifs2 software also identifies sequence motifs that cluster together with the binding motifs of TDP-43. We conclude that TDP-43 directly regulates diverse types of pre-mRNA processing according to common position-dependent principles.

Insights into the design and interpretation of iCLIP experiments.

BACKGROUND: Ultraviolet (UV) crosslinking and immunoprecipitation (CLIP) identifies the sites on RNAs that are in direct contact with RNA-binding proteins (RBPs). Several variants of CLIP exist, which require different computational approaches for analysis. This variety of approaches can create challenges for a novice user and can hamper insights from multi-study comparisons. Here, we produce data with multiple variants of CLIP and evaluate the data with various computational methods to better understand their suitability. RESULTS: We perform experiments for PTBP1 and eIF4A3 using individual-nucleotide resolution CLIP (iCLIP), employing either UV-C or photoactivatable 4-thiouridine (4SU) combined with UV-A crosslinking and compare the results with published data. As previously noted, the positions of complementary DNA (cDNA)-starts depend on cDNA length in several iCLIP experiments and we now find that this is caused by constrained cDNA-ends, which can result from the sequence and structure constraints of RNA fragmentation. These constraints are overcome when fragmentation by RNase I is efficient and when a broad cDNA size range is obtained. Our study also shows that if RNase does not efficiently cut within the binding sites, the original CLIP method is less capable of identifying the longer binding sites of RBPs. In contrast, we show that a broad size range of cDNAs in iCLIP allows the cDNA-starts to efficiently delineate the complete RNA-binding sites. CONCLUSIONS: We demonstrate the advantage of iCLIP and related methods that can amplify cDNAs that truncate at crosslink sites and we show that computational analyses based on cDNAs-starts are appropriate for such methods.

Exon Junction Complexes Show a Distributional Bias toward Alternatively Spliced mRNAs and against mRNAs Coding for Ribosomal Proteins.

The exon junction complex (EJC) connects spliced mRNAs to posttranscriptional processes including RNA localization, transport, and regulated degradation. Here, we provide a comprehensive analysis of bona fide EJC binding sites across the transcriptome including all four RNA binding EJC components eIF4A3, BTZ, UPF3B, and RNPS1. Integration of these data sets permits definition of high-confidence EJC deposition sites as well as assessment of whether EJC heterogeneity drives alternative nonsense-mediated mRNA decay pathways. Notably, BTZ (MLN51 or CASC3) emerges as the EJC subunit that is almost exclusively bound to sites 20-24 nucleotides upstream of exon-exon junctions, hence defining EJC positions. By contrast, eIF4A3, UPF3B, and RNPS1 display additional RNA binding sites suggesting accompanying non-EJC functions. Finally, our data show that EJCs are largely distributed across spliced RNAs in an orthodox fashion, with two notable exceptions: an EJC deposition bias in favor of alternatively spliced transcripts and against the mRNAs that encode ribosomal proteins.

The differential expression of alternatively polyadenylated transcripts is a common stress-induced response mechanism that modulates mammalian mRNA expression in a quantitative and qualitative fashion.

Stress adaptation plays a pivotal role in biological processes and requires tight regulation of gene expression. In this study, we explored the effect of cellular stress on mRNA polyadenylation and investigated the implications of regulated polyadenylation site usage on mammalian gene expression. High-confidence polyadenylation site mapping combined with global pre-mRNA and mRNA expression profiling revealed that stress induces an accumulation of genes with differentially expressed polyadenylated mRNA isoforms in human cells. Specifically, stress provokes a global trend in polyadenylation site usage toward decreased utilization of promoter-proximal poly(A) sites in introns or ORFs and increased utilization of promoter-distal polyadenylation sites in intergenic regions. This extensively affects gene expression beyond regulating mRNA abundance by changing mRNA length and by altering the configuration of open reading frames. Our study highlights the impact of post-transcriptional mechanisms on stress-dependent gene regulation and reveals the differential expression of alternatively polyadenylated transcripts as a common stress-induced mechanism in mammalian cells.

Comprehensive Identification of RNA-Binding Domains in Human Cells.

Mammalian cells harbor more than a thousand RNA-binding proteins (RBPs), with half of these employing unknown modes of RNA binding. We developed RBDmap to determine the RNA-binding sites of native RBPs on a proteome-wide scale. We identified 1,174 binding sites within 529 HeLa cell RBPs, discovering numerous RNA-binding domains (RBDs). Catalytic centers or protein-protein interaction domains are in close relationship with RNA-binding sites, invoking possible effector roles of RNA in the control of protein function. Nearly half of the RNA-binding sites map to intrinsically disordered regions, uncovering unstructured domains as prevalent partners in protein-RNA interactions. RNA-binding sites represent hot spots for defined posttranslational modifications such as lysine acetylation and tyrosine phosphorylation, suggesting metabolic and signal-dependent regulation of RBP function. RBDs display a high degree of evolutionary conservation and incidence of Mendelian mutations, suggestive of important functional roles. RBDmap thus yields profound insights into native protein-RNA interactions in living cells.

SR proteins are NXF1 adaptors that link alternative RNA processing to mRNA export.

Nuclear export factor 1 (NXF1) exports mRNA to the cytoplasm after recruitment to mRNA by specific adaptor proteins. How and why cells use numerous different export adaptors is poorly understood. Here we critically evaluate members of the SR protein family (SRSF1-7) for their potential to act as NXF1 adaptors that couple pre-mRNA processing to mRNA export. Consistent with this proposal, >1000 endogenous mRNAs required individual SR proteins for nuclear export in vivo. To address the mechanism, transcriptome-wide RNA-binding profiles of NXF1 and SRSF1-7 were determined in parallel by individual-nucleotide-resolution UV cross-linking and immunoprecipitation (iCLIP). Quantitative comparisons of RNA-binding sites showed that NXF1 and SR proteins bind mRNA targets at adjacent sites, indicative of cobinding. SRSF3 emerged as the most potent NXF1 adaptor, conferring sequence specificity to RNA binding by NXF1 in last exons. Interestingly, SRSF3 and SRSF7 were shown to bind different sites in last exons and regulate 3' untranslated region length in an opposing manner. Both SRSF3 and SRSF7 promoted NXF1 recruitment to mRNA. Thus, SRSF3 and SRSF7 couple alternative splicing and polyadenylation to NXF1-mediated mRNA export, thereby controlling the cytoplasmic abundance of transcripts with alternative 3' ends.