An intrauterine genomic classifier reliably delineates the location of nonviable pregnancies.

OBJECTIVE: To compare the intrauterine gene expression signatures of women with surgically confirmed ectopic pregnancy (ECT) and those of women with miscarriage to inform the development of a genomic classifier for the reliable delineation of pregnancy location in women with clinically nonviable pregnancies of unknown location (NV-PULs). DESIGN: Discovery-based prospective cohort study. SETTING: Academic medical center. PATIENT(S): Women with clinically nonviable early pregnancy to include abnormal intrauterine pregnancy (AIUP), ECT, or NV-PUL. INTERVENTION(S): Endometrial (EM) pipelle sampling of the uterus was conducted at the time of scheduled surgery for clinically nonviable early pregnancy (dilation and curettage, manual vacuum aspiration, or laparoscopy). All pregnancy locations were surgically and/or histologically confirmed as intrauterine or ectopic. MAIN OUTCOME MEASURE(S): Gene expression profiles as determined by array hybridization, quantitative real-time polymerase chain reaction, and nCounter technology. RESULT(S): Intrauterine samples were obtained by EM pipelle from 27 women undergoing surgery for a clinically nonviable early pregnancy. Comparison of array-based global gene expression signatures from women with histologically confirmed ECT versus AIUP revealed 61 differentially expressed genes from which the 5 most informative were included in the pregnancy location classifier. All 5 genes (C20orf85, LRRC46, RSPH4A, WDR49, and ZBBX) were cilia-associated and showed increased expression in pipelle samples from women with ECT relative to expression in samples from women with AIUP. The 5-gene classifier demonstrated an average area under the receiver operator characteristic curve of 0.97 for the detection of ECT. In an external test set composed of publicly available EM pipelle-based gene expression data from a study with similar ECT and AIUP cohorts (n = 19), the classifier revealed an average area under the receiver operator characteristic curve of 0.84. CONCLUSION(S): Consistently increased expression of cilia-associated genes in the uterine cavity of women with ECT provides a reliable molecular signal for the delineation of pregnancy location in women with clinically assessed NV-PUL. A classifier consisting of the 5 most informative cilia-associated genes demonstrated 91% (42/46) accuracy in predicting the pregnancy location.

Opposing immune and genetic mechanisms shape oncogenic programs in synovial sarcoma.

Synovial sarcoma (SyS) is an aggressive neoplasm driven by the SS18-SSX fusion, and is characterized by low T cell infiltration. Here, we studied the cancer-immune interplay in SyS using an integrative approach that combines single-cell RNA sequencing (scRNA-seq), spatial profiling and genetic and pharmacological perturbations. scRNA-seq of 16,872 cells from 12 human SyS tumors uncovered a malignant subpopulation that marks immune-deprived niches in situ and is predictive of poor clinical outcomes in two independent cohorts. Functional analyses revealed that this malignant cell state is controlled by the SS18-SSX fusion, is repressed by cytokines secreted by macrophages and T cells, and can be synergistically targeted with a combination of HDAC and CDK4/CDK6 inhibitors. This drug combination enhanced malignant-cell immunogenicity in SyS models, leading to induced T cell reactivity and T cell-mediated killing. Our study provides a blueprint for investigating heterogeneity in fusion-driven malignancies and demonstrates an interplay between immune evasion and oncogenic processes that can be co-targeted in SyS and potentially in other malignancies.

Reproducible Bioconductor workflows using browser-based interactive notebooks and containers.

Objective: Bioinformatics publications typically include complex software workflows that are difficult to describe in a manuscript. We describe and demonstrate the use of interactive software notebooks to document and distribute bioinformatics research. We provide a user-friendly tool, BiocImageBuilder, that allows users to easily distribute their bioinformatics protocols through interactive notebooks uploaded to either a GitHub repository or a private server. Materials and methods: We present four different interactive Jupyter notebooks using R and Bioconductor workflows to infer differential gene expression, analyze cross-platform datasets, process RNA-seq data and KinomeScan data. These interactive notebooks are available on GitHub. The analytical results can be viewed in a browser. Most importantly, the software contents can be executed and modified. This is accomplished using Binder, which runs the notebook inside software containers, thus avoiding the need to install any software and ensuring reproducibility. All the notebooks were produced using custom files generated by BiocImageBuilder. Results: BiocImageBuilder facilitates the publication of workflows with a point-and-click user interface. We demonstrate that interactive notebooks can be used to disseminate a wide range of bioinformatics analyses. The use of software containers to mirror the original software environment ensures reproducibility of results. Parameters and code can be dynamically modified, allowing for robust verification of published results and encouraging rapid adoption of new methods. Conclusion: Given the increasing complexity of bioinformatics workflows, we anticipate that these interactive software notebooks will become as necessary for documenting software methods as traditional laboratory notebooks have been for documenting bench protocols, and as ubiquitous.

SpermBase: A Database for Sperm-Borne RNA Contents.

Since their discovery ~three decades ago, sperm-borne RNAs, both large/small and coding/noncoding, have been reported in multiple organisms, and some have been implicated in spermatogenesis, early development, and epigenetic inheritance. Despite these advances, isolation, quantification and annotation of sperm-borne RNAs remain nontrivial. The yields and subspecies of sperm-borne RNAs isolated from sperm can vary drastically depending on the methods used, and no cross-species analyses of sperm RNA contents have ever been conducted using a standardized sperm RNA isolation protocol. To address these issues, we developed a simple RNA isolation method that is applicable to sperm of various species, thus allowing for reliable interspecies comparisons. Based on RNA-Seq analyses, we established SpermBase (www.spermbase.org), a database dedicated to sperm-borne RNA profiling of multiple species. Currently, SpermBase contains large and small RNA expression data for mouse, rat, rabbit and human total sperm and sperm heads. By analyzing large and small RNAs for conserved features, we found that many sperm-borne RNA species were conserved across all four species analyzed, and among the conserved small RNAs, sperm-borne tsRNAs and miRNAs can target a large number of genes known to be critical for early development.

Sperm-borne miRNAs and endo-siRNAs are important for fertilization and preimplantation embryonic development.

Although it is believed that mammalian sperm carry small noncoding RNAs (sncRNAs) into oocytes during fertilization, it remains unknown whether these sperm-borne sncRNAs truly have any function during fertilization and preimplantation embryonic development. Germline-specific Dicer and Drosha conditional knockout (cKO) mice produce gametes (i.e. sperm and oocytes) partially deficient in miRNAs and/or endo-siRNAs, thus providing a unique opportunity for testing whether normal sperm (paternal) or oocyte (maternal) miRNA and endo-siRNA contents are required for fertilization and preimplantation development. Using the outcome of intracytoplasmic sperm injection (ICSI) as a readout, we found that sperm with altered miRNA and endo-siRNA profiles could fertilize wild-type (WT) eggs, but embryos derived from these partially sncRNA-deficient sperm displayed a significant reduction in developmental potential, which could be rescued by injecting WT sperm-derived total or small RNAs into ICSI embryos. Disrupted maternal transcript turnover and failure in early zygotic gene activation appeared to associate with the aberrant miRNA profiles in Dicer and Drosha cKO spermatozoa. Overall, our data support a crucial function of paternal miRNAs and/or endo-siRNAs in the control of the transcriptomic homeostasis in fertilized eggs, zygotes and two-cell embryos. Given that supplementation of sperm RNAs enhances both the developmental potential of preimplantation embryos and the live birth rate, it might represent a novel means to improve the success rate of assisted reproductive technologies in fertility clinics.

In silico identification of novel endo-siRNAs.

Many classes of small noncoding RNAs (sncRNAs), such as microRNAs (miRNAs) and endogenous small interfering RNAs (endo-siRNAs), have been identified as important regulators of gene expression. Endo-siRNAs represent an integral part of the endogenous RNAi pathway and have been identified in multiple organisms and cell types. Wide adoption of the next-generation deep sequencing (NGS)-based sncRNA profiling has made the identification of novel sncRNA species more accessible. However, it remains a challenge to identify novel endo-siRNAs that are not collected in the current endo-siRNA databases. We have developed an in silico method for identification of novel endo-siRNAs using small RNA NGS data. Here, we describe our protocol in detail.

Computer-assisted annotation of small RNA transcriptomes.

Small noncoding RNAs (sncRNAs) are widely expressed in the cell of almost all known species. Most sncRNAs appear to have regulatory roles, ranging from facilitating RNA production and modifications (e.g., snoRNAs) to control of mRNA stability and translational efficiency (e.g., miRNAs and endo-siRNA) and to transposon silencing (e.g., piRNAs). The affordability and efficiency of next-generation RNA deep sequencing (RNA-Seq) technologies have made sncRNA deep sequencing (sncRNA-Seq) analyses a routine in biomedical research. SncRNA-Seq analyses generate millions of reads and gigabytes of data; annotation of sncRNA-Seq data remains challenging due to a lack of comprehensive sncRNA annotation pipelines. To solve this problem, we have developed a computer-assisted sncRNA annotation pipeline, which uses open-source software and allows for not only proper classification of known sncRNAs, but also discovery of novel sncRNA species. In this chapter, we describe our sncRNA annotation protocol in detail.

A novel class of somatic small RNAs similar to germ cell pachytene PIWI-interacting small RNAs.

PIWI-interacting RNAs (piRNAs) are small noncoding RNAs that bind PIWI family proteins exclusively expressed in the germ cells of mammalian gonads. MIWI2-associated piRNAs are essential for silencing transposons during primordial germ cell development, and MIWI-bound piRNAs are required for normal spermatogenesis during adulthood in mice. Although piRNAs have long been regarded as germ cell-specific, increasing lines of evidence suggest that somatic cells also express piRNA-like RNAs (pilRNAs). Here, we report the detection of abundant pilRNAs in somatic cells, which are similar to MIWI-associated piRNAs mainly expressed in pachytene spermatocytes and round spermatids in the testis. Based on small RNA deep sequencing and quantitative PCR analyses, pilRNA expression is dynamic and displays tissue specificity. Although pilRNAs are similar to pachytene piRNAs in both size and genomic origins, they have a distinct ping-pong signature. Furthermore, pilRNA biogenesis appears to utilize a yet to be identified pathway, which is different from all currently known small RNA biogenetic pathways. In addition, pilRNAs appear to preferentially target the 3'-UTRs of mRNAs in a partially complementary manner. Our data suggest that pilRNAs, as an integral component of the small RNA transcriptome in somatic cell lineages, represent a distinct population of small RNAs that may have functions similar to germ cell piRNAs.

Murine follicular development requires oocyte DICER, but not DROSHA.

Both DICER and DROSHA are RNase III enzymes involved in the biogenesis of small noncoding RNAs. DROSHA cleaves the stem-loop portion of the primary miRNAs and produces precursor miRNAs in the nucleus, whereas DICER processes double-stranded RNA precursors into mature miRNAs and endogenous small interference RNAs in the cytoplasm. Selective inactivation of Dicer in growing oocytes of primary follicles leads to female infertility due to oocyte spindle defects. However, it remains unknown if oocyte Dicer expression in the fetal ovary is required for proper follicular development in the postnatal ovary. Moreover, the role of Drosha in folliculogenesis has never been investigated. Here, we report that conditional knockout of Dicer in prophase I oocytes of the fetal ovary led to compromised folliculogenesis, premature ovarian failure, and female infertility in the adult ovary, whereas selective inactivation of Drosha in oocytes of either the fetal or the developing ovary had no effects on normal folliculogenesis and female fertility in adulthood. Our data indicate that oocyte DICER expression in the fetal ovary is required, and oocyte DROSHA is dispensable, for postnatal follicular development and female fertility in adulthood.

The mitochondrial genome encodes abundant small noncoding RNAs.

Small noncoding RNAs identified thus far are all encoded by the nuclear genome. Here, we report that the murine and human mitochondrial genomes encode thousands of small noncoding RNAs, which are predominantly derived from the sense transcripts of the mitochondrial genes (host genes), and we termed these small RNAs mitochondrial genome-encoded small RNAs (mitosRNAs). DICER inactivation affected, but did not completely abolish mitosRNA production. MitosRNAs appear to be products of currently unidentified mitochondrial ribonucleases. Overexpression of mitosRNAs enhanced expression levels of their host genes in vitro, and dysregulated mitosRNA expression was generally associated with aberrant mitochondrial gene expression in vivo. Our data demonstrate that in addition to 37 known mitochondrial genes, the mammalian mitochondrial genome also encodes abundant mitosRNAs, which may play an important regulatory role in the control of mitochondrial gene expression in the cell.

Computer-assisted annotation of murine Sertoli cell small RNA transcriptome.

Mammalian genomes encode a large number of small noncoding RNAs (sncRNAs) that play regulatory roles during development and adulthood by affecting gene expression. Several sncRNA species, including microRNAs (miRNAs), piwi-interacting RNAs (piRNAs), endogenous small interfering RNAs (endo-siRNAs), and small nucleolar RNAs (snoRNAs), are abundantly expressed in the testis and required for normal testicular development and spermatogenesis. To evaluate global changes in sncRNA expression, the next-generation sequencing (NGS)-based sncRNA transcriptomic analysis has become routine, because it allows rapid determination of the small RNA transcriptome of a particular testicular cell type. However, annotation of small RNA NGS reads can be challenging due to the volume of reads obtained, which is usually in the millions. Therefore, we developed a computer-assisted sncRNA annotation protocol that could identify not only known sncRNAs but also previously uncharacterized ones. Using this protocol, we annotated NGS reads of a Sertoli cell sncRNA library, and we report to our knowledge the first comprehensive annotation of the sncRNA transcriptome of immature murine Sertoli cells. Moreover, the computer-assisted sncRNA annotation pipeline that we report is applicable for annotating NGS reads derived from other cell types and/or sequencing platforms.

The RNase III enzyme DROSHA is essential for microRNA production and spermatogenesis.

DROSHA is a nuclear RNase III enzyme responsible for cleaving primary microRNAs (miRNAs) into precursor miRNAs and thus is essential for the biogenesis of canonical miRNAs. DICER is a cytoplasmic RNase III enzyme that not only cleaves precursor miRNAs to produce mature miRNAs but also dissects naturally formed/synthetic double-stranded RNAs to generate small interfering RNAs (siRNAs). To investigate the role of canonical miRNA and/or endogenous siRNA production in spermatogenesis, we generated Drosha or Dicer conditional knock-out (cKO) mouse lines by inactivating Drosha or Dicer exclusively in spermatogenic cells in postnatal testes using the Cre-loxp strategy. Both Drosha and Dicer cKO males were infertile due to disrupted spermatogenesis characterized by depletion of spermatocytes and spermatids leading to oligoteratozoospermia or azoospermia. The developmental course of spermatogenic disruptions was similar at morphological levels between Drosha and Dicer cKO males, but Drosha cKO testes appeared to be more severe in spermatogenic disruptions than Dicer cKO testes. Microarray analyses revealed transcriptomic differences between Drosha- and Dicer-null pachytene spermatocytes or round spermatids. Although levels of sex-linked mRNAs were mildly elevated, meiotic sex chromosome inactivation appeared to have occurred normally. Our data demonstrate that unlike DICER, which is required for the biogenesis of several small RNA species, DROSHA is essential mainly for the canonical miRNA production, and DROSHA-mediated miRNA production is essential for normal spermatogenesis and male fertility.