Detecting haplotype-specific transcript variation in long reads with FLAIR2.

BACKGROUND: RNA-seq has brought forth significant discoveries regarding aberrations in RNA processing, implicating these RNA variants in a variety of diseases. Aberrant splicing and single nucleotide variants (SNVs) in RNA have been demonstrated to alter transcript stability, localization, and function. In particular, the upregulation of ADAR, an enzyme that mediates adenosine-to-inosine editing, has been previously linked to an increase in the invasiveness of lung adenocarcinoma cells and associated with splicing regulation. Despite the functional importance of studying splicing and SNVs, the use of short-read RNA-seq has limited the community's ability to interrogate both forms of RNA variation simultaneously. RESULTS: We employ long-read sequencing technology to obtain full-length transcript sequences, elucidating cis-effects of variants on splicing changes at a single molecule level. We develop a computational workflow that augments FLAIR, a tool that calls isoform models expressed in long-read data, to integrate RNA variant calls with the associated isoforms that bear them. We generate nanopore data with high sequence accuracy from H1975 lung adenocarcinoma cells with and without knockdown of ADAR. We apply our workflow to identify key inosine isoform associations to help clarify the prominence of ADAR in tumorigenesis. CONCLUSIONS: Ultimately, we find that a long-read approach provides valuable insight toward characterizing the relationship between RNA variants and splicing patterns.

Probing chromatin accessibility with small molecule DNA intercalation and nanopore sequencing.

Genome-wide identification of chromatin organization and structure has been generally probed by measuring accessibility of the underlying DNA to nucleases or methyltransferases. These methods either only observe the positioning of a single nucleosome or rely on large enzymes to modify or cleave the DNA. We developed adduct sequencing (Add-seq), a method to probe chromatin accessibility by treating chromatin with the small molecule angelicin, which preferentially intercalates into DNA not bound to core nucleosomes. We show that Nanopore sequencing of the angelicin-modified DNA is possible and allows visualization and analysis of long single molecules with distinct chromatin structure. The angelicin modification can be detected from the Nanopore current signal data using a neural network model trained on unmodified and modified chromatin-free DNA. Applying Add-seq to Saccharomyces cerevisiae nuclei, we identified expected patterns of accessibility around annotated gene loci in yeast. We also identify individual clusters of single molecule reads displaying different chromatin structure at specific yeast loci, which demonstrates heterogeneity in the chromatin structure of the yeast population. Thus, using Add-seq, we are able to profile DNA accessibility in the yeast genome across long molecules.

DeepSea is an efficient deep-learning model for single-cell segmentation and tracking in time-lapse microscopy.

Time-lapse microscopy is the only method that can directly capture the dynamics and heterogeneity of fundamental cellular processes at the single-cell level with high temporal resolution. Successful application of single-cell time-lapse microscopy requires automated segmentation and tracking of hundreds of individual cells over several time points. However, segmentation and tracking of single cells remain challenging for the analysis of time-lapse microscopy images, in particular for widely available and non-toxic imaging modalities such as phase-contrast imaging. This work presents a versatile and trainable deep-learning model, termed DeepSea, that allows for both segmentation and tracking of single cells in sequences of phase-contrast live microscopy images with higher precision than existing models. We showcase the application of DeepSea by analyzing cell size regulation in embryonic stem cells.

Detecting haplotype-specific transcript variation in long reads with FLAIR2.

Background: RNA-Seq has brought forth significant discoveries regarding aberrations in RNA processing, implicating these RNA variants in a variety of diseases. Aberrant splicing and single nucleotide variants in RNA have been demonstrated to alter transcript stability, localization, and function. In particular, the upregulation of ADAR, an enzyme which mediates adenosine-to-inosine editing, has been previously linked to an increase in the invasiveness of lung ADC cells and associated with splicing regulation. Despite the functional importance of studying splicing and SNVs, short read RNA-Seq has limited the community's ability to interrogate both forms of RNA variation simultaneously. Results: We employed long-read technology to obtain full-length transcript sequences, elucidating cis-effects of variants on splicing changes at a single molecule level. We have developed a computational workflow that augments FLAIR, a tool that calls isoform models expressed in long-read data, to integrate RNA variant calls with the associated isoforms that bear them. We generated nanopore data with high sequence accuracy of H1975 lung adenocarcinoma cells with and without knockdown of ADAR. We applied our workflow to identify key inosine-isoform associations to help clarify the prominence of ADAR in tumorigenesis. Conclusions: Ultimately, we find that a long-read approach provides valuable insight toward characterizing the relationship between RNA variants and splicing patterns.

Full-length transcript alterations in human bronchial epithelial cells with U2AF1 S34F mutations.

U2AF1 is one of the most recurrently mutated splicing factors in lung adenocarcinoma and has been shown to cause transcriptome-wide pre-mRNA splicing alterations; however, the full-length altered mRNA isoforms associated with the mutation are largely unknown. To better understand the impact U2AF1 has on full-length isoform fate and function, we conducted high-throughput long-read cDNA sequencing from isogenic human bronchial epithelial cells with and without a U2AF1 S34F mutation. We identified 49,366 multi-exon transcript isoforms, more than half of which did not match GENCODE or short-read-assembled isoforms. We found 198 transcript isoforms with significant expression and usage changes relative to WT, only 68% of which were assembled by short reads. Expression of isoforms from immune-related genes is largely down-regulated in mutant cells and without observed splicing changes. Finally, we reveal that isoforms likely targeted by nonsense-mediated decay are down-regulated in U2AF1 S34F cells, suggesting that isoform changes may alter the translational output of those affected genes. Altogether, our work provides a resource of full-length isoforms associated with U2AF1 S34F in lung cells.

Mutant KRAS regulates transposable element RNA and innate immunity via KRAB zinc-finger genes.

RAS genes are the most frequently mutated oncogenes in cancer, yet the effects of oncogenic RAS signaling on the noncoding transcriptome remain unclear. We analyzed the transcriptomes of human airway and bronchial epithelial cells transformed with mutant KRAS to define the landscape of KRAS-regulated noncoding RNAs. We find that oncogenic KRAS signaling upregulates noncoding transcripts throughout the genome, many of which arise from transposable elements (TEs). These TE RNAs exhibit differential expression, are preferentially released in extracellular vesicles, and are regulated by KRAB zinc-finger (KZNF) genes, which are broadly downregulated in mutant KRAS cells and lung adenocarcinomas in vivo. Moreover, mutant KRAS induces an intrinsic IFN-stimulated gene (ISG) signature that is often seen across many different cancers. Our results indicate that mutant KRAS remodels the repetitive noncoding transcriptome, demonstrating the broad scope of intracellular and extracellular RNAs regulated by this oncogenic signaling pathway.

Full-length transcript characterization of SF3B1 mutation in chronic lymphocytic leukemia reveals downregulation of retained introns.

While splicing changes caused by somatic mutations in SF3B1 are known, identifying full-length isoform changes may better elucidate the functional consequences of these mutations. We report nanopore sequencing of full-length cDNA from CLL samples with and without SF3B1 mutation, as well as normal B cell samples, giving a total of 149 million pass reads. We present FLAIR (Full-Length Alternative Isoform analysis of RNA), a computational workflow to identify high-confidence transcripts, perform differential splicing event analysis, and differential isoform analysis. Using nanopore reads, we demonstrate differential 3' splice site changes associated with SF3B1 mutation, agreeing with previous studies. We also observe a strong downregulation of intron retention events associated with SF3B1 mutation. Full-length transcript analysis links multiple alternative splicing events together and allows for better estimates of the abundance of productive versus unproductive isoforms. Our work demonstrates the potential utility of nanopore sequencing for cancer and splicing research.

A two-step probing method to compare lysine accessibility across macromolecular complex conformations.

Structural models of large and dynamic molecular complexes are appearing in increasing numbers, in large part because of recent technical advances in cryo-electron microscopy. However, the inherent complexity of such biological assemblies comprising dozens of moving parts often limits the resolution of structural models and leaves the puzzle as to how each functional configuration transitions to the next. Orthogonal biochemical information is crucial to understanding the molecular interactions that drive those rearrangements. We present a two-step method for chemical probing detected by tandem mass-spectrometry to globally assess the reactivity of lysine residues within purified macromolecular complexes. Because lysine side chains often balance the negative charge of RNA in ribonucleoprotein complexes, the method is especially useful for detecting changes in protein-RNA interactions. By probing the E. coli 30S ribosome subunit, we established that the reactivity pattern of lysine residues quantitatively reflects structure models derived from X-ray crystallography. We also used the strategy to assess differences in three conformations of purified human spliceosomes in the context of recent cryo-electron microscopy models. Our results demonstrate that the probing method yields powerful biochemical information that helps contextualize architectural rearrangements of intermediate resolution structures of macromolecular complexes, often solved in multiple conformations.

Prp8 positioning of U5 snRNA is linked to 5' splice site recognition.

Prp8 is an essential protein that regulates spliceosome assembly and conformation during pre-mRNA splicing. Recent cryo-EM structures of the spliceosome model Prp8 as a scaffold for the spliceosome's catalytic U snRNA components. Using a new amino acid probing strategy, we identified a dynamic region in human Prp8 that is positioned to stabilize the pre-mRNA in the spliceosome active site through interactions with U5 snRNA. Mutagenesis of the identified Prp8 residues in yeast indicates a role in 5' splice site recognition. Genetic interactions with spliceosome proteins Isy1, which buttresses the intron branch point, and Snu114, a regulatory GTPase that directly contacts Prp8, further corroborate a role for the same Prp8 residues in substrate positioning and activation. Together the data suggest that adjustments in interactions between Prp8 and U5 snRNA help establish proper positioning of the pre-mRNA into the active site to enhance 5' splice site fidelity.

High-Throughput Small RNA Sequencing Enhanced by AlkB-Facilitated RNA de-Methylation (ARM-Seq).

N 1-methyladenosine (m1A), N 3-methylcytidine (m3C), and N 1-methylguanosine (m1G) are common in transfer RNA (tRNA) and tRNA-derived fragments. These modifications alter Watson-Crick base-pairing, and cause pauses or stops during reverse transcription required for most high-throughput RNA sequencing protocols, resulting in inefficient detection of methyl-modified RNAs. Here, we describe a procedure to demethylate RNAs containing m1A, m3C, or m1G using the Escherichia coli dealkylating enzyme AlkB, along with instructions for subsequent processing with widely used protocols for small RNA sequencing.

ARM-seq: AlkB-facilitated RNA methylation sequencing reveals a complex landscape of modified tRNA fragments.

High-throughput RNA sequencing has accelerated discovery of the complex regulatory roles of small RNAs, but RNAs containing modified nucleosides may escape detection when those modifications interfere with reverse transcription during RNA-seq library preparation. Here we describe AlkB-facilitated RNA methylation sequencing (ARM-seq), which uses pretreatment with Escherichia coli AlkB to demethylate N(1)-methyladenosine (m(1)A), N(3)-methylcytidine (m(3)C) and N(1)-methylguanosine (m(1)G), all commonly found in tRNAs. Comparative methylation analysis using ARM-seq provides the first detailed, transcriptome-scale map of these modifications and reveals an abundance of previously undetected, methylated small RNAs derived from tRNAs. ARM-seq demonstrates that tRNA fragments accurately recapitulate the m(1)A modification state for well-characterized yeast tRNAs and generates new predictions for a large number of human tRNAs, including tRNA precursors and mitochondrial tRNAs. Thus, ARM-seq provides broad utility for identifying previously overlooked methyl-modified RNAs, can efficiently monitor methylation state and may reveal new roles for tRNA fragments as biomarkers or signaling molecules.

Quantitative analysis of the transcription control mechanism.

Gene transcription requires a sequence of promoter state transitions, including chromatin remodeling, assembly of the transcription machinery, and clearance of the promoter by RNA polymerase. The rate-limiting steps in this sequence are regulated by transcriptional activators that bind at specific promoter elements. As the transition kinetics of individual promoters cannot be observed, the identity of the activator-controlled steps has remained a matter of speculation. In this study, we investigated promoter chromatin structure, and the intrinsic noise of expression over a wide range of expression values for the PHO5 gene of yeast. Interpretation of our results with regard to a stochastic model of promoter chromatin remodeling and gene expression suggests that the regulatory architecture of the gene expression process is measurably reflected in its intrinsic noise profile. Our chromatin structure and noise analyses indicate that the activator of PHO5 transcription stimulates the rates of promoter nucleosome disassembly, and assembly of the transcription machinery after nucleosome removal, but no other rates of the expression process.

Standard operating procedure for the collection and preparation of voucher plant specimens for use in the nutraceutical industry.

A vital part of the development of any standardized protocol for the extraction of plant-derived crude extracts to be used in herbal medicine or nutritional supplementation is proper documentation of the original botanical source of the extract via acquisition of a voucher specimen. The purpose of this document is to serve as an accepted protocol for voucher specimen collection, handling, and storage, with specific guidelines to address commercial and research uses.

Estrogen-triggered delays in mammary gland gene expression during the estrous cycle: evidence for a novel timing system.

During the estrous cycle and beginning in estrus, the mammary gland undergoes pregnancy-like development that depends on transcriptional regulation by the estrogen and progesterone receptors (ER, PR) and Pax-2 as well as the action of the growth factors Wnt-4 and RANKL. In this report, we first describe the decay and delayed expression of ERalpha, PR, and Pax-2 proteins as well as depression of Wnt-4 and RANKL mRNA coincident with the strong estrogen surge in proestrus. In time-course studies using ovari-ectomized mice, a single estrogen injection replicated these delays and caused an 18 h delay in Wnt-4 expression. Molecular time-delay systems are at the core of cellular cycles, most notably the circadian clock, and depend on proteasome degradation of transcriptional regulators that exhibit dedicated timing functions. The cytoplasmic dynamics of these regulators govern delay duration through negative transcription/translation feedback loops. A proteasome inhibitor, PS-341, blocked estrogen-stimulated ERalpha, PR, and Pax-2 decay and proteasome chymotryptic activity, assayed using a fluorogenic substrate, was elevated in proestrus correlating with the depletion of the transcription factors. The 18-h delay in Wnt-4 induction corresponded to the turnover time of Pax-2 protein in the cytoplasm and was eliminated in Pax-2 knockout mammary tissue, demonstrating that Pax-2 has a unique timing function. The patterns of estrogen-triggered ERalpha, PR, and Pax-2 turnover were consistent with a negative transcriptional feedback. Retarding the expression of ERalpha, PR, and Pax-2 may optimize preparations for pregnancy by coordinating expression of critical receptors and transcription factors with rising estrogen and progesterone levels in estrus. The estrogen surge in proestrus has no defined mammotropic function. This study provides the first evidence that it is a synchronizing signal triggering proteasome-dependent turnover of mammary gland ERalpha, PR, and Pax-2. We hypothesize that the delays reflect a previously unrecognized timing system, which is present in all ovarian target tissues.