Open-source Tools for Training Resources - OTTR.

Data science and informatics tools are developing at a blistering rate, but their users often lack the educational background or resources to efficiently apply the methods to their research. Training resources and vignettes that accompany these tools often deprecate because their maintenance is not prioritized by funding, giving teams little time to devote to such endeavors. Our group has developed Open-source Tools for Training Resources (OTTR) to offer greater efficiency and flexibility for creating and maintaining these training resources. OTTR empowers creators to customize their work and allows for a simple workflow to publish using multiple platforms. OTTR allows content creators to publish training material to multiple massive online learner communities using familiar rendering mechanics. OTTR allows the incorporation of pedagogical practices like formative and summative assessments in the form of multiple choice questions and fill in the blank problems that are automatically graded. No local installation of any software is required to begin creating content with OTTR. Thus far, 15 training courses have been created with OTTR repository template. By using the OTTR system, the maintenance workload for updating these courses across platforms has been drastically reduced. For more information about OTTR and how to get started, go to ottrproject.org.

Endovascular image-guided sampling of tumor-draining veins provides an enriched source of oncological biomarkers.

Introduction: Circulating tumor-derived biomarkers can potentially impact cancer management throughout the continuum of care. This small exploratory study aimed to assess the relative levels of such biomarkers in the tumor-draining vascular beds in patients with solid tumors compared to levels in their peripheral veins. Methods: Using an endovascular image-guided approach, we obtained blood samples from peripheral veins and other vascular compartments-including the most proximal venous drainage from solid tumors-from a set of nine oncology patients with various primary and metastatic malignancies. We then interrogated these samples for a panel of oncological biomarkers, including circulating tumor cells (CTCs), exosome-derived microRNAs (miRNAs), circulating tumor DNA (ctDNA) mutations, and certain cancer-related proteins/biochemical markers. Results: We found substantially higher levels of CTCs, certain miRNAs, and specific ctDNA mutations in samples from vascular beds closer to the tumor compared with those from peripheral veins and also noted that some of these signals were altered by treatment procedures. Discussion: Our results indicate that tumor-proximal venous samples are highly enriched for some oncological biomarkers and may allow for more robust molecular analysis than peripheral vein samples.

Motivation, inclusivity, and realism should drive data science education.

Data science education provides tremendous opportunities but remains inaccessible to many communities. Increasing the accessibility of data science to these communities not only benefits the individuals entering data science, but also increases the field's innovation and potential impact as a whole. Education is the most scalable solution to meet these needs, but many data science educators lack formal training in education. Our group has led education efforts for a variety of audiences: from professional scientists to high school students to lay audiences. These experiences have helped form our teaching philosophy which we have summarized into three main ideals: 1) motivation, 2) inclusivity, and 3) realism. 20 we also aim to iteratively update our teaching approaches and curriculum as we find ways to better reach these ideals. In this manuscript we discuss these ideals as well practical ideas for how to implement these philosophies in the classroom.

Characterization of an RNA binding protein interactome reveals a context-specific post-transcriptional landscape of MYC-amplified medulloblastoma.

Pediatric medulloblastoma (MB) is the most common solid malignant brain neoplasm, with Group 3 (G3) MB representing the most aggressive subgroup. MYC amplification is an independent poor prognostic factor in G3 MB, however, therapeutic targeting of the MYC pathway remains limited and alternative therapies for G3 MB are urgently needed. Here we show that the RNA-binding protein, Musashi-1 (MSI1) is an essential mediator of G3 MB in both MYC-overexpressing mouse models and patient-derived xenografts. MSI1 inhibition abrogates tumor initiation and significantly prolongs survival in both models. We identify binding targets of MSI1 in normal neural and G3 MB stem cells and then cross referenced these data with unbiased large-scale screens at the transcriptomic, translatomic and proteomic levels to systematically dissect its functional role. Comparative integrative multi-omic analyses of these large datasets reveal cancer-selective MSI1-bound targets sharing multiple MYC associated pathways, providing a valuable resource for context-specific therapeutic targeting of G3 MB.

Discovery and functional interrogation of SARS-CoV-2 protein-RNA interactions.

The COVID-19 pandemic is caused by severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2). The betacoronvirus has a positive sense RNA genome which encodes for several RNA binding proteins. Here, we use enhanced crosslinking and immunoprecipitation to investigate SARS-CoV-2 protein interactions with viral and host RNAs in authentic virus-infected cells. SARS-CoV-2 proteins, NSP8, NSP12, and nucleocapsid display distinct preferences to specific regions in the RNA viral genome, providing evidence for their shared and separate roles in replication, transcription, and viral packaging. SARS-CoV-2 proteins expressed in human lung epithelial cells bind to 4773 unique host coding RNAs. Nine SARS-CoV-2 proteins upregulate target gene expression, including NSP12 and ORF9c, whose RNA substrates are associated with pathways in protein N-linked glycosylation ER processing and mitochondrial processes. Furthermore, siRNA knockdown of host genes targeted by viral proteins in human lung organoid cells identify potential antiviral host targets across different SARS-CoV-2 variants. Conversely, NSP9 inhibits host gene expression by blocking mRNA export and dampens cytokine productions, including interleukin-1α/ß. Our viral protein-RNA interactome provides a catalog of potential therapeutic targets and offers insight into the etiology of COVID-19 as a safeguard against future pandemics.

Discovery and functional interrogation of SARS-CoV-2 protein-RNA interactions.

The COVID-19 pandemic is caused by severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2). The betacoronvirus has a positive sense RNA genome which encodes for several RNA binding proteins. Here, we use enhanced crosslinking and immunoprecipitation to investigate SARS-CoV-2 protein interactions with viral and host RNAs in authentic virus-infected cells. SARS-CoV-2 proteins, NSP8, NSP12, and nucleocapsid display distinct preferences to specific regions in the RNA viral genome, providing evidence for their shared and separate roles in replication, transcription, and viral packaging. SARS-CoV-2 proteins expressed in human lung epithelial cells bind to 4773 unique host coding RNAs. Nine SARS-CoV-2 proteins upregulate target gene expression, including NSP12 and ORF9c, whose RNA substrates are associated with pathways in protein N-linked glycosylation ER processing and mitochondrial processes. Furthermore, siRNA knockdown of host genes targeted by viral proteins in human lung organoid cells identify potential antiviral host targets across different SARS-CoV-2 variants. Conversely, NSP9 inhibits host gene expression by blocking mRNA export and dampens cytokine productions, including interleukin-1α/ß. Our viral protein-RNA interactome provides a catalog of potential therapeutic targets and offers insight into the etiology of COVID-19 as a safeguard against future pandemics.

Inverting the model of genomics data sharing with the NHGRI Genomic Data Science Analysis, Visualization, and Informatics Lab-space.

The NHGRI Genomic Data Science Analysis, Visualization, and Informatics Lab-space (AnVIL; https://anvilproject.org) was developed to address a widespread community need for a unified computing environment for genomics data storage, management, and analysis. In this perspective, we present AnVIL, describe its ecosystem and interoperability with other platforms, and highlight how this platform and associated initiatives contribute to improved genomic data sharing efforts. The AnVIL is a federated cloud platform designed to manage and store genomics and related data, enable population-scale analysis, and facilitate collaboration through the sharing of data, code, and analysis results. By inverting the traditional model of data sharing, the AnVIL eliminates the need for data movement while also adding security measures for active threat detection and monitoring and provides scalable, shared computing resources for any researcher. We describe the core data management and analysis components of the AnVIL, which currently consists of Terra, Gen3, Galaxy, RStudio/Bioconductor, Dockstore, and Jupyter, and describe several flagship genomics datasets available within the AnVIL. We continue to extend and innovate the AnVIL ecosystem by implementing new capabilities, including mechanisms for interoperability and responsible data sharing, while streamlining access management. The AnVIL opens many new opportunities for analysis, collaboration, and data sharing that are needed to drive research and to make discoveries through the joint analysis of hundreds of thousands to millions of genomes along with associated clinical and molecular data types.

Non-microRNA binding competitively inhibits LIN28 regulation.

RNA binding protein (RBP) expression is finite. For RBPs that are vastly outnumbered by their potential target sites, a simple competition for binding can set the magnitude of post-transcriptional control. Here, we show that LIN28, best known for its direct regulation of let-7 miRNA biogenesis, is also indirectly regulated by its widespread binding of non-miRNA transcripts. Approximately 99% of LIN28 binding sites are found on non-miRNA transcripts, like protein coding and ribosomal RNAs. These sites are bound specifically and strongly, but they do not appear to mediate direct post-transcriptional regulation. Instead, non-miRNA sites act to sequester LIN28 protein and effectively change its functional availability, thus impeding the regulation of let-7 in cells. Together, these data show that the binding properties of the transcriptome broadly influence the ability of an RBP to mediate changes in RNA metabolism and gene expression.

Large-scale tethered function assays identify factors that regulate mRNA stability and translation.

The molecular functions of the majority of RNA-binding proteins (RBPs) remain unclear, highlighting a major bottleneck to a full understanding of gene expression regulation. Here, we develop a plasmid resource of 690 human RBPs that we subject to luciferase-based 3'-untranslated-region tethered function assays to pinpoint RBPs that regulate RNA stability or translation. Enhanced UV-cross-linking and immunoprecipitation of these RBPs identifies thousands of endogenous mRNA targets that respond to changes in RBP level, recapitulating effects observed in tethered function assays. Among these RBPs, the ubiquitin-associated protein 2-like (UBAP2L) protein interacts with RNA via its RGG domain and cross-links to mRNA and rRNA. Fusion of UBAP2L to RNA-targeting CRISPR-Cas9 demonstrates programmable translational enhancement. Polysome profiling indicates that UBAP2L promotes translation of target mRNAs, particularly global regulators of translation. Our tethering survey allows rapid assignment of the molecular activity of proteins, such as UBAP2L, to specific steps of mRNA metabolism.

A point mutation decouples the lipid transfer activities of microsomal triglyceride transfer protein.

Apolipoprotein B-containing lipoproteins (B-lps) are essential for the transport of hydrophobic dietary and endogenous lipids through the circulation in vertebrates. Zebrafish embryos produce large numbers of B-lps in the yolk syncytial layer (YSL) to move lipids from yolk to growing tissues. Disruptions in B-lp production perturb yolk morphology, readily allowing for visual identification of mutants with altered B-lp metabolism. Here we report the discovery of a missense mutation in microsomal triglyceride transfer protein (Mtp), a protein that is essential for B-lp production. This mutation of a conserved glycine residue to valine (zebrafish G863V, human G865V) reduces B-lp production and results in yolk opacity due to aberrant accumulation of cytoplasmic lipid droplets in the YSL. However, this phenotype is milder than that of the previously reported L475P stalactite (stl) mutation. MTP transfers lipids, including triglycerides and phospholipids, to apolipoprotein B in the ER for B-lp assembly. In vitro lipid transfer assays reveal that while both MTP mutations eliminate triglyceride transfer activity, the G863V mutant protein unexpectedly retains ~80% of phospholipid transfer activity. This residual phospholipid transfer activity of the G863V mttp mutant protein is sufficient to support the secretion of small B-lps, which prevents intestinal fat malabsorption and growth defects observed in the mttpstl/stl mutant zebrafish. Modeling based on the recent crystal structure of the heterodimeric human MTP complex suggests the G865V mutation may block triglyceride entry into the lipid-binding cavity. Together, these data argue that selective inhibition of MTP triglyceride transfer activity may be a feasible therapeutic approach to treat dyslipidemia and provide structural insight for drug design. These data also highlight the power of yolk transport studies to identify proteins critical for B-lp biology.

Muscle stem cell renewal suppressed by Gas1 can be reversed by GDNF in mice.

Muscle undergoes progressive weakening and regenerative dysfunction with age due in part to the functional decline of skeletal muscle stem cells (MuSCs). MuSCs are heterogeneous but whether their gene expression changes with age and the implication of such changes are unclear. Here we show that in mice, Growth arrest-specific gene 1 (Gas1) is expressed in a small subset of young MuSCs with its expression progressively increasing in larger fractions of MuSCs later in life. Over-expression of Gas1 in young MuSCs and inactivation of Gas1 in aged MuSCs support that Gas1 reduces the quiescence and self-renewal capacity of MuSCs. Gas1 reduces Ret signaling, which is required for MuSC quiescence and self-renewal. Indeed, we show that the Ret ligand, Glial Cell-Derived Neurotrophic Factor (GDNF), can counteract Gas1 by stimulating Ret signaling and enhancing MuSC self-renewal and regeneration, thus improving muscle function. We propose that strategies aimed to target this pathway can be exploited to improve the regenerative decline of muscle stem cells.

A Transcriptome-wide Translational Program Defined by LIN28B Expression Level.

LIN28 RNA binding proteins are dynamically expressed throughout mammalian development and during disease. However, it remains unclear how changes in LIN28 expression define patterns of post-transcriptional gene regulation. Here we show that LIN28 expression level is a key variable that sets the magnitude of protein translation. By systematically varying LIN28B protein levels in human cells, we discovered a dose-dependent divergence in transcriptome-wide ribosome occupancy that enabled the formation of two discrete translational subpopulations composed of nearly all expressed genes. This bifurcation in gene expression was mediated by a redistribution in Argonaute association, from let-7 to non-let-7 microRNA families, resulting in a global shift in cellular miRNA activity. Post-transcriptional effects were scaled across the physiological LIN28 expression range. Together, these data highlight the central importance of RBP expression level and its ability to encode regulation.

Hijacking Oogenesis Enables Massive Propagation of LINE and Retroviral Transposons.

Although animals have evolved multiple mechanisms to suppress transposons, "leaky" mobilizations that cause mutations and diseases still occur. This suggests that transposons employ specific tactics to accomplish robust propagation. By directly tracking mobilization, we show that, during a short and specific time window of oogenesis, retrotransposons achieve massive amplification via a cell-type-specific targeting strategy. Retrotransposons rarely mobilize in undifferentiated germline stem cells. However, as oogenesis proceeds, they utilize supporting nurse cells-which are highly polyploid and eventually undergo apoptosis-as factories to massively manufacture invading products. Moreover, retrotransposons rarely integrate into nurse cells themselves but, instead, via microtubule-mediated transport, they preferentially target the DNA of the interconnected oocytes. Blocking microtubule-dependent intercellular transport from nurse cells significantly alleviates damage to the oocyte genome. Our data reveal that parasitic genomic elements can efficiently hijack a host developmental process to propagate robustly, thereby driving evolutionary change and causing disease.

The Hippo pathway effector proteins YAP and TAZ have both distinct and overlapping functions in the cell.

The Hippo pathway plays an important role in regulating tissue homeostasis, and its effectors, the transcriptional co-activators Yes-associated protein (YAP) and WW domain-containing transcription regulator 1 (WWTR1 or TAZ), are responsible for mediating the vast majority of its physiological functions. Although YAP and TAZ are thought to be largely redundant and similarly regulated by Hippo signaling, they have developmental, structural, and physiological differences that suggest they may differ in their regulation and downstream functions. To better understand the functions of YAP and TAZ in the Hippo pathway, using CRISPR/Cas9, we generated YAP KO, TAZ KO, and YAP/TAZ KO cell lines in HEK293A cells. We evaluated them in response to many environmental conditions and stimuli and used RNA-Seq to compare their transcriptional profiles. We found that YAP inactivation has a greater effect on cellular physiology (namely, cell spreading, volume, granularity, glucose uptake, proliferation, and migration) than TAZ inactivation. However, functional redundancy between YAP and TAZ was also observed. In summary, our findings confirm that the Hippo pathway effectors YAP and TAZ are master regulators for multiple cellular processes but also reveal that YAP has a stronger influence than TAZ.

mRNA processing in mutant zebrafish lines generated by chemical and CRISPR-mediated mutagenesis produces unexpected transcripts that escape nonsense-mediated decay.

As model organism-based research shifts from forward to reverse genetics approaches, largely due to the ease of genome editing technology, a low frequency of abnormal phenotypes is being observed in lines with mutations predicted to lead to deleterious effects on the encoded protein. In zebrafish, this low frequency is in part explained by compensation by genes of redundant or similar function, often resulting from the additional round of teleost-specific whole genome duplication within vertebrates. Here we offer additional explanations for the low frequency of mutant phenotypes. We analyzed mRNA processing in seven zebrafish lines with mutations expected to disrupt gene function, generated by CRISPR/Cas9 or ENU mutagenesis methods. Five of the seven lines showed evidence of altered mRNA processing: one through a skipped exon that did not lead to a frame shift, one through nonsense-associated splicing that did not lead to a frame shift, and three through the use of cryptic splice sites. These results highlight the need for a methodical analysis of the mRNA produced in mutant lines before making conclusions or embarking on studies that assume loss of function as a result of a given genomic change. Furthermore, recognition of the types of adaptations that can occur may inform the strategies of mutant generation.

Importance of anatomical dominance in the evaluation of coronary dilatation in Kawasaki disease.

Introduction In Kawasaki disease, although coronary dilatation is attributed to vasculitis, the effect of myocardial inflammation is underestimated. Coronary dilatations are determined by Z-scores, which do not take into account dominance. The aim of the present study was to describe the impact of coronary dominance on dilatation in Kawasaki disease. METHODS: We performed a retrospective analysis of coronary dilatations according to angiography categorisation of dominance. RESULTS: Of 28 patients (2.6 [0.2-10.1] years), right dominance was present in 15 patients and left in 13. Early dilatation was present in all patients, of whom 11 were ipsilateral to the dominant segment and 17 contralateral. Ipsilateral dilatations were present at diagnosis (9/11 versus 6/17, p=0.02) compared with contralateral dilatations, which developed 2 weeks after diagnosis (9/11 versus 16/17, p=0.29). Coronary artery Z-scores of patients with contralateral dilatation increased at 2 weeks, before returning to baseline values (2.0±2.2 at diagnosis, 4.1±1.8 at 2 weeks, 1.8±1.2 at 3-6 months, p=0.001), compared with patients with ipsilateral dilatation in whom Z-scores were maximal at diagnosis and remained stable (3.0±0.9, 2.7±1.1 and 2.6±1.5, respectively, p=0.13). Dominant coronary artery Z-scores were higher compared with non-dominant segments at diagnosis (3.0±0.9 versus 1.0±0.8, p<0.001) and at late follow-up (2.6±1.5 versus 0.4±1.4, p=0.002) in patients with ipsilateral dilatation. CONCLUSION: Progression of coronary dilatation after diagnosis may be a sign of dilatation secondary to vasculitis, as opposed to regression of Z-scores in ipsilateral dilatations, probably related to physiological vasodilatation in response to carditis. This needs to be validated in larger studies against vasculitic and myocardial inflammatory markers.

S1-DRIP-seq identifies high expression and polyA tracts as major contributors to R-loop formation.

R loops form when transcripts hybridize to homologous DNA on chromosomes, yielding a DNA:RNA hybrid and a displaced DNA single strand. R loops impact the genome of many organisms, regulating chromosome stability, gene expression, and DNA repair. Understanding the parameters dictating R-loop formation in vivo has been hampered by the limited quantitative and spatial resolution of current genomic strategies for mapping R loops. We report a novel whole-genome method, S1-DRIP-seq (S1 nuclease DNA:RNA immunoprecipitation with deep sequencing), for mapping hybrid-prone regions in budding yeast Saccharomyces cerevisiae Using this methodology, we identified ∼800 hybrid-prone regions covering 8% of the genome. Given the pervasive transcription of the yeast genome, this result suggests that R-loop formation is dictated by characteristics of the DNA, RNA, and/or chromatin. We successfully identified two features highly predictive of hybrid formation: high transcription and long homopolymeric dA:dT tracts. These accounted for >60% of the hybrid regions found in the genome. We demonstrated that these two factors play a causal role in hybrid formation by genetic manipulation. Thus, the hybrid map generated by S1-DRIP-seq led to the identification of the first global genomic features causal for R-loop formation in yeast.

Blurred Boundaries: The RNA Binding Protein Lin28A Is Also an Epigenetic Regulator.

Lin28A is best known as a post-transcriptional regulator of gene expression. In this issue, Zeng et al. (2016) show that Lin28A has an unexpected role as an epigenetic regulator of DNA.

Brf1 posttranscriptionally regulates pluripotency and differentiation responses downstream of Erk MAP kinase.

AU-rich element mRNA-binding proteins (AUBPs) are key regulators of development, but how they are controlled and what functional roles they play depends on cellular context. Here, we show that Brf1 (zfp36l1), an AUBP from the Zfp36 protein family, operates downstream of FGF/Erk MAP kinase signaling to regulate pluripotency and cell fate decision making in mouse embryonic stem cells (mESCs). FGF/Erk MAP kinase signaling up-regulates Brf1, which disrupts the expression of core pluripotency-associated genes and attenuates mESC self-renewal without inducing differentiation. These regulatory effects are mediated by rapid and direct destabilization of Brf1 targets, such as Nanog mRNA. Enhancing Brf1 expression does not compromise mESC pluripotency but does preferentially regulate mesendoderm commitment during differentiation, accelerating the expression of primitive streak markers. Together, these studies demonstrate that FGF signals use targeted mRNA degradation by Brf1 to enable rapid posttranscriptional control of gene expression in mESCs.