LAFITE Reveals the Complexity of Transcript Isoforms in Subcellular Fractions.

Characterization of the subcellular distribution of RNA is essential for understanding the molecular basis of biological processes. Here, the subcellular nanopore direct RNA-sequencing (DRS) of four lung cancer cell lines (A549, H1975, H358, and HCC4006) is performed, coupled with a computational pipeline, Low-abundance Aware Full-length Isoform clusTEr (LAFITE), to comprehensively analyze the full-length cytoplasmic and nuclear transcriptome. Using additional DRS and orthogonal data sets, it is shown that LAFITE outperforms current methods for detecting full-length transcripts, particularly for low-abundance isoforms that are usually overlooked due to poor read coverage. Experimental validation of six novel isoforms exclusively identified by LAFITE further confirms the reliability of this pipeline. By applying LAFITE to subcellular DRS data, the complexity of the nuclear transcriptome is revealed in terms of isoform diversity, 3'-UTR usage, m6A modification patterns, and intron retention. Overall, LAFITE provides enhanced full-length isoform identification and enables a high-resolution view of the RNA landscape at the isoform level.

A Chromosome-level assembly of the Japanese eel genome, insights into gene duplication and chromosomal reorganization.

Japanese eels (Anguilla japonica) are commercially important species, harvested extensively for food. Currently, this and related species (American and European eels) are challenging to breed on a commercial basis. As a result, the wild stock is used for aquaculture. Moreover, climate change, habitat loss, water pollution, and altered ocean currents affect eel populations negatively. Accordingly, the International Union for Conservation of Nature lists Japanese eels as endangered and on its red list. Here we presented a high-quality genome assembly for Japanese eels and demonstrated that large chromosome reorganizations occurred in the events of third-round whole-genome duplications (3R-WRDs). Several chromosomal fusions and fissions have reduced the ancestral protochromosomal number of 25 to 19 in the Anguilla lineage. A phylogenetic analysis of the expanded gene families showed that the olfactory receptors (group δ and ζ genes) and voltage-gated Ca2+ channels expanded significantly. Both gene families are crucial for olfaction and neurophysiology. Additional tandem and proximal duplications occurred following 3R-WGD to acquire immune-related genes for an adaptive advantage against various pathogens. The Japanese eel assembly presented here can be used to study other Anguilla species relating to evolution and conservation.

Enhanced transcriptome-wide RNA G-quadruplex sequencing for low RNA input samples with rG4-seq 2.0.

BACKGROUND: RNA G-quadruplexes (rG4s) are non-canonical structural motifs that have diverse functional and regulatory roles, for instance in transcription termination, alternative splicing, mRNA localization and stabilization, and translational process. We recently developed the RNA G-quadruplex structure sequencing (rG4-seq) technique and described rG4s in both eukaryotic and prokaryotic transcriptomes. However, rG4-seq suffers from a complicated gel purification step and limited PCR product yield, thus requiring a high amount of RNA input, which limits its applicability in more physiologically or clinically relevant studies often characterized by the limited availability of biological material and low RNA abundance. Here, we redesign and enhance the workflow of rG4-seq to address this issue. RESULTS: We developed rG4-seq 2.0 by introducing a new ssDNA adapter containing deoxyuridine during library preparation to enhance library quality with no gel purification step, less PCR amplification cycles and higher yield of PCR products. We demonstrate that rG4-seq 2.0 produces high-quality cDNA libraries that support reliable and reproducible rG4 identification at varying RNA inputs, including RNA mounts as low as 10 ng. rG4-seq 2.0 also improved the rG4-seq calling outcome and nucleotide bias in rG4 detection persistent in rG4-seq 1.0. We further provide in vitro mapping of rG4 in the HEK293T cell line, and recommendations for assessing RNA input and sequencing depth for individual rG4 studies based on transcript abundance. CONCLUSIONS: rG4-seq 2.0 can improve the identification and study of rG4s in low abundance transcripts, and our findings can provide insights to optimize cDNA library preparation in other related methods.

An RNA G-Quadruplex Structure within the ADAR 5'UTR Interacts with DHX36 Helicase to Regulate Translation.

RNA G-quadruplex (rG4) structures in the 5' untranslated region (5'UTR) play crucial roles in fundamental cellular processes. ADAR is an important enzyme that binds to double-strand RNA and accounts for the conversion of Adenosine to Inosine in RNA editing. However, so far there is no report on the formation and regulatory role of rG4 on ADAR expression. Here, we identify and characterize a thermostable rG4 structure within the 5'UTR of the ADAR1 mRNA and demonstrate its formation and inhibitory role on translation in reporter gene and native gene constructs. We reveal rG4-specific helicase DHX36 interacts with this rG4 in vitro and in cells under knockdown and knockout conditions by GTFH (G-quadruplex-triggered fluorogenic hybridization) probes and modulates translation in an rG4-dependent manner. Our results further substantiate the rG4 structure-DHX36 protein interaction in cells and highlight rG4 to be a key player in controlling ADAR1 translation.

G-quadruplex RNA motifs influence gene expression in the malaria parasite Plasmodium falciparum.

G-quadruplexes are non-helical secondary structures that can fold in vivo in both DNA and RNA. In human cells, they can influence replication, transcription and telomere maintenance in DNA, or translation, transcript processing and stability of RNA. We have previously showed that G-quadruplexes are detectable in the DNA of the malaria parasite Plasmodium falciparum, despite a very highly A/T-biased genome with unusually few guanine-rich sequences. Here, we show that RNA G-quadruplexes can also form in P. falciparum RNA, using rG4-seq for transcriptome-wide structure-specific RNA probing. Many of the motifs, detected here via the rG4seeker pipeline, have non-canonical forms and would not be predicted by standard in silico algorithms. However, in vitro biophysical assays verified formation of non-canonical motifs. The G-quadruplexes in the P. falciparum transcriptome are frequently clustered in certain genes and associated with regions encoding low-complexity peptide repeats. They are overrepresented in particular classes of genes, notably those that encode PfEMP1 virulence factors, stress response genes and DNA binding proteins. In vitro translation experiments and in vivo measures of translation efficiency showed that G-quadruplexes can influence the translation of P. falciparum mRNAs. Thus, the G-quadruplex is a novel player in post-transcriptional regulation of gene expression in this major human pathogen.

RNA G-quadruplexes (rG4s): genomics and biological functions.

G-quadruplexes (G4s) are non-classical DNA or RNA secondary structures that have been first observed decades ago. Over the years, these four-stranded structural motifs have been demonstrated to have significant regulatory roles in diverse biological processes, but challenges remain in detecting them globally and reliably. Compared to DNA G4s (dG4s), the study of RNA G4s (rG4s) has received less attention until recently. In this review, we will summarize the innovative high-throughput methods recently developed to detect rG4s on a transcriptome-wide scale, highlight the many novel and important functions of rG4 being discovered in vivo across the tree of life, and discuss the key biological questions to be addressed in the near future.

rG4-seeker enables high-confidence identification of novel and non-canonical rG4 motifs from rG4-seq experiments.

We recently developed the rG4-seq method to detect and map in vitro RNA G-quadruplex (rG4s) structures on a transcriptome-wide scale. rG4-seq of purified human HeLa RNA has revealed many non-canonical rG4s and the effects adjacent sequences have on rG4 formation. In this study, we aimed to improve the outcomes and false-positive discrimination in rG4-seq experiments using a bioinformatic approach. By establishing connections between rG4-seq library preparation chemistry and the underlying properties of sequencing data, we identified how to mitigate indigenous sampling errors and background noise in rG4-seq. We applied these findings to develop a novel bioinformatics pipeline named rG4-seeker (https://github.com/TF-Chan-Lab/rG4-seeker), which uses tailored noise models to autonomously assess and optimize rG4 detections in a replicate-independent manner. Compared with previous methods, rG4-seeker exhibited better false-positive discrimination and improved sensitivity for non-canonical rG4s. Using rG4-seeker, we identified novel features in rG4 formation that were missed previously. rG4-seeker provides a reliable and sensitive approach for rG4-seq investigations, laying the foundations for further elucidation of rG4 biology.

Author Correction: Systematic evaluation and optimization of the experimental steps in RNA G-quadruplex structure sequencing.

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

Characterization of Hepatocellular Carcinoma Cell Lines Using a Fractionation-Then-Sequencing Approach Reveals Nuclear-Enriched HCC-Associated lncRNAs.

Background: Advances in sequencing technologies have greatly improved our understanding of long noncoding RNA (lncRNA). These transcripts with lengths of >200 nucleotides may play significant regulatory roles in various biological processes. Importantly, the dysregulation of better characterized lncRNAs has been associated with multiple types of cancers, including hepatocellular carcinoma (HCC). There are many studies on altered lncRNA expression levels, very few, however, have focused on their subcellular localizations, from which accumulating evidences have indicated their close relationships to lncRNA functions. A transcriptome-wide investigation of the subcellular distributions of lncRNAs might thus provide new insights into their roles and functions in cancers. Results: In this study, we subjected eight patient-derived HCC cell lines to subcellular fractionation and independently sequenced RNAs from the nuclear and cytoplasmic compartments. With the integration of tumor and tumor-adjacent RNA-seq datasets of liver hepatocellular carcinoma (LIHC) from The Cancer Genome Atlas (TCGA), de novo transcriptome assembly and differential expression analysis were conducted successively and identified 26 nuclear-enriched HCC-associated lncRNAs shared between the HCC samples and the TCGA datasets, including the reported cancer driver PXN-AS1. The majority of nuclear-enriched HCC-associated lncRNAs were associated with the survival outcomes of HCC patients, exhibited characteristics similar to those of many experimentally supported HCC prognostic lncRNAs, and were co-expressed with protein-coding genes that have been linked to disease progression in various cancer types. Conclusion: We adopted a fractionation-then-sequencing approach on multiple patient-derived HCC samples and identified nuclear-enriched, HCC-associated lncRNAs that could serve as important targets for HCC diagnosis and therapeutic development. This approach could be widely applicable to other studies into the disease etiologies of lncRNA.

Systematic evaluation and optimization of the experimental steps in RNA G-quadruplex structure sequencing.

cDNA library preparation is important for many high-throughput sequencing applications, such as RNA G-quadruplex structure sequencing (rG4-seq). A systematic evaluation of the procedures of the experimental pipeline, however, is lacking. Herein, we perform a comprehensive assessment of the 5 key experimental steps involved in the cDNA library preparation of rG4-seq, and identify better reaction conditions and/or enzymes to carry out each of these key steps. Notably, we apply the improved methods to fragmented cellular RNA, and show reduced RNA input requirement, lower transcript abundance variations between biological replicates, as well as lower transcript coverage bias when compared to prior arts. In addition, the time to perform these steps is substantially reduced to hours. Our method and results can be directly applied in protocols that require cDNA library preparation, and provide insights to the further development of simple and efficient cDNA library preparation for different biological applications.

Genome maps across 26 human populations reveal population-specific patterns of structural variation.

Large structural variants (SVs) in the human genome are difficult to detect and study by conventional sequencing technologies. With long-range genome analysis platforms, such as optical mapping, one can identify large SVs (>2 kb) across the genome in one experiment. Analyzing optical genome maps of 154 individuals from the 26 populations sequenced in the 1000 Genomes Project, we find that phylogenetic population patterns of large SVs are similar to those of single nucleotide variations in 86% of the human genome, while ~2% of the genome has high structural complexity. We are able to characterize SVs in many intractable regions of the genome, including segmental duplications and subtelomeric, pericentromeric, and acrocentric areas. In addition, we discover ~60 Mb of non-redundant genome content missing in the reference genome sequence assembly. Our results highlight the need for a comprehensive set of alternate haplotypes from different populations to represent SV patterns in the genome.

A microarray MEMS device for biolistic delivery of vaccine and drug powders.

We report a biolistic technology platform for physical delivery of particle formulations of drugs or vaccines using parallel arrays of microchannels, which generate highly collimated jets of particles with high spatial resolution. Our approach allows for effective delivery of therapeutics sequentially or concurrently (in mixture) at a specified target location or treatment area. We show this new platform enables the delivery of a broad range of particles with various densities and sizes into both in vitro and ex vivo skin models. Penetration depths of ∼1 mm have been achieved following a single ejection of 200 µg high-density gold particles, as well as 13.6 µg low-density polystyrene-based particles into gelatin-based skin simulants at 70 psi inlet gas pressure. Ejection of multiple shots at one treatment site enabled deeper penetration of ∼3 mm in vitro, and delivery of a higher dose of 1 mg gold particles at similar inlet gas pressure. We demonstrate that particle penetration depths can be optimized in vitro by adjusting the inlet pressure of the carrier gas, and dosing is controlled by drug reservoirs that hold precise quantities of the payload, which can be ejected continuously or in pulses. Future investigations include comparison between continuous versus pulsatile payload deliveries. We have successfully delivered plasmid DNA (pDNA)-coated gold particles (1.15 µm diameter) into ex vivo murine and porcine skin at low inlet pressures of ∼30 psi. Integrity analysis of these pDNA-coated gold particles confirmed the preservation of full-length pDNA after each particle preparation and jetting procedures. This technology platform provides distinct capabilities to effectively deliver a broad range of particle formulations into skin with specially designed high-speed microarray ejector nozzles.

Epitope switching as a novel escape mechanism of HIV to CCR5 monoclonal antibodies.

In passaging experiments, we isolated HIV strains resistant to MAb3952, a chemokine (C-C motif) receptor 5 (CCR5) monoclonal antibody (MAb) that binds to the second extracellular domain (extracellular loop 2 [ECL-2]) of CCR5. MAb3952-resistant viruses remain CCR5-tropic and are cross-resistant to a second ECL-2-specific antibody. Surprisingly, MAb3952-resistant viruses were more susceptible to RoAb13, a CCR5 antibody binding to the N terminus of CCR5. Using CCR5 receptor mutants, we show that MAb3952-resistant virus strains preferentially use the N terminus of CCR5, while the wild-type viruses preferentially use ECL-2. We propose this switch in the CCR5 binding site as a novel mechanism of HIV resistance.

CD4-BFFI: a novel, bifunctional HIV-1 entry inhibitor with high and broad antiviral potency.

Resistance to antiretroviral drugs is a common problem in the treatment of HIV-1-infected patients. To overcome resistance, we generated a novel, bifunctional HIV-1 entry inhibitor by combining the anti-CD4 monoclonal antibody (mAb) 6314 with a fusion inhibitor similar to T-651 (anti-CD4 mAb based BiFunctional Fusion Inhibitor, CD4-BFFI). CD4-BFFI has potent antiviral activity against a multitude of HIV-1 isolates independent of their co-receptor usage and genetic background. It has higher antiviral potency compared to the fusion inhibitor T-651 or the anti-CD4 mAb 6314 used independently. More importantly, every HIV-1 strain tested was fully inhibited by CD4-BFFI while many strains were only partially inhibited by 6314. CD4-BFFI also retained antiviral potency against virus strains resistant to two fusion inhibitors, a CCR5 antagonist and an anti-CCR5 mAb. Pre-incubation of cells with a saturating concentration of anti-CD4 mAbs reduced the antiviral potency of CD4-BFFI, suggesting that binding of CD4-BFFI to the cell surface via its CD4 mAb portion is required for the antiviral potency of its fusion inhibitor moiety. Collectively, we present a novel HIV-1 inhibitor with a dual mode of action and excellent antiviral potency against wildtype and entry-inhibitor resistant virus strains suggesting that CD4-BFFI may have a high barrier to resistance.

Intracellular calcium waves in bone cell networks under single cell nanoindentation.

In this study, bone cells were successfully cultured into a micropatterned network with dimensions close to that of in vivo osteocyte networks using microcontact printing and self-assembled monolyers (SAMs). The optimal geometric parameters for the formation of these networks were determined in terms of circle diameters and line widths. Bone cells patterned in these networks were also able to form gap junctions with each other, shown by immunofluorescent staining for the gap junction protein connexin 43, as well as the transfer of gap-junction permeable calcein-AM dye. We have demonstrated for the first time, that the intracellular calcium response of a single bone cell indented in this bone cell network, can be transmitted to neighboring bone cells through multiple calcium waves. Furthermore, the propagation of these calcium waves was diminished with increased cell separation distance. Thus, this study provides new experimental data that support the idea of osteocyte network memory of mechanical loading similar to memory in neural networks.