Blood RNA alternative splicing events as diagnostic biomarkers for infectious disease.

Assays detecting blood transcriptome changes are studied for infectious disease diagnosis. Blood-based RNA alternative splicing (AS) events, which have not been well characterized in pathogen infection, have potential normalization and assay platform stability advantages over gene expression for diagnosis. Here, we present a computational framework for developing AS diagnostic biomarkers. Leveraging a large prospective cohort of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection and whole-blood RNA sequencing (RNA-seq) data, we identify a major functional AS program switch upon viral infection. Using an independent cohort, we demonstrate the improved accuracy of AS biomarkers for SARS-CoV-2 diagnosis compared with six reported transcriptome signatures. We then optimize a subset of AS-based biomarkers to develop microfluidic PCR diagnostic assays. This assay achieves nearly perfect test accuracy (61/62 = 98.4%) using a naive principal component classifier, significantly more accurate than a gene expression PCR assay in the same cohort. Therefore, our RNA splicing computational framework enables a promising avenue for host-response diagnosis of infection.

The human salivary RNA transcriptome revealed by massively parallel sequencing.

BACKGROUND: Evaluation of the salivary transcriptome is an emerging diagnostic technology with discriminatory power for disease detection. This study explored massively parallel sequencing for providing nucleotide-level sequence information for each RNA in saliva. METHODS: Transcriptome profiling with the SOLiD™ system was applied to RNA isolated from unstimulated cell-free saliva (CFS) and whole saliva (WS) from healthy human volunteers. Sequenced reads were aligned to human genome build 18 and the Human Oral Microbiome Database (HOMD). RESULTS: Massively parallel sequencing enabled the acquisition of complete sequence information for each nucleotide position of the human salivary transcriptome through alignment to multiple sequence databases. Approximately 20%-25% of the sequenced reads from CFS aligned to the human genome, and approximately 30% of the sequenced reads aligned to the HOMD. We detected the expression of >4000 coding and noncoding genes in CFS and WS. Gene expression at different genome loci showed that the structural integrity of the transcripts for the annotated genes was preserved in saliva. CONCLUSIONS: A single measurement provided RNA sequence information of gene transcript abundance for both coding and noncoding RNAs and identified sequences from >400 different microbial species within a single sample. Contrary to previous data suggesting that salivary RNA is highly fragmented, in our study the structural integrity of RNA was preserved. The high degree of sequence alignment to annotated exons and introns for each of the respective reference genomes, with sequence coverage spanning the full length of the genes, provides strong evidence that the salivary transcriptome is a complex RNA network.

RNA-Seq mapping and detection of gene fusions with a suffix array algorithm.

High-throughput RNA sequencing enables quantification of transcripts (both known and novel), exon/exon junctions and fusions of exons from different genes. Discovery of gene fusions-particularly those expressed with low abundance- is a challenge with short- and medium-length sequencing reads. To address this challenge, we implemented an RNA-Seq mapping pipeline within the LifeScope software. We introduced new features including filter and junction mapping, annotation-aided pairing rescue and accurate mapping quality values. We combined this pipeline with a Suffix Array Spliced Read (SASR) aligner to detect chimeric transcripts. Performing paired-end RNA-Seq of the breast cancer cell line MCF-7 using the SOLiD system, we called 40 gene fusions among over 120,000 splicing junctions. We validated 36 of these 40 fusions with TaqMan assays, of which 25 were expressed in MCF-7 but not the Human Brain Reference. An intra-chromosomal gene fusion involving the estrogen receptor alpha gene ESR1, and another involving the RPS6KB1 (Ribosomal protein S6 kinase beta-1) were recurrently expressed in a number of breast tumor cell lines and a clinical tumor sample.

RecQ helicase stimulates both DNA catenation and changes in DNA topology by topoisomerase III.

Together, RecQ helicase and topoisomerase III (Topo III) of Escherichia coli comprise a potent DNA strand passage activity that can catenate covalently closed DNA (Harmon, F. G., DiGate, R. J., and Kowalczykowski, S. C. (1999) Mol. Cell 3, 611-620). Here we directly assessed the structure of the catenated DNA species formed by RecQ helicase and Topo III using atomic force microscopy. The images show complex catenated DNA species involving crossovers between multiple double-stranded DNA molecules that are consistent with full catenanes. E. coli single-stranded DNA-binding protein significantly stimulated both the topoisomerase activity of Topo III alone and the DNA strand passage activity of RecQ helicase and Topo III. Titration data suggest that an intermediate of the RecQ helicase unwinding process, perhaps a RecQ helicase-DNA fork, is the target for Topo III action. Catenated DNA is the predominant product under conditions of molecular crowding; however, we also discovered that RecQ helicase and single-stranded DNA-binding protein greatly stimulated the intramolecular strand passage ("supercoiling") activity of Topo III, as revealed by changes in the linking number of uncatenated DNA. Together our results demonstrate that RecQ helicase and Topo III function together to comprise a potent and concerted single-strand DNA passage activity that can mediate both catenation-decatenation processes and changes in DNA topology.

High accuracy genotyping directly from genomic DNA using a rolling circle amplification based assay.

BACKGROUND: Rolling circle amplification of ligated probes is a simple and sensitive means for genotyping directly from genomic DNA. SNPs and mutations are interrogated with open circle probes (OCP) that can be circularized by DNA ligase when the probe matches the genotype. An amplified detection signal is generated by exponential rolling circle amplification (ERCA) of the circularized probe. The low cost and scalability of ligation/ERCA genotyping makes it ideally suited for automated, high throughput methods. RESULTS: A retrospective study using human genomic DNA samples of known genotype was performed for four different clinically relevant mutations: Factor V Leiden, Factor II prothrombin, and two hemochromatosis mutations, C282Y and H63D. Greater than 99% accuracy was obtained genotyping genomic DNA samples from hundreds of different individuals. The combined process of ligation/ERCA was performed in a single tube and produced fluorescent signal directly from genomic DNA in less than an hour. In each assay, the probes for both normal and mutant alleles were combined in a single reaction. Multiple ERCA primers combined with a quenched-peptide nucleic acid (Q-PNA) fluorescent detection system greatly accellerated the appearance of signal. Probes designed with hairpin structures reduced misamplification. Genotyping accuracy was identical from either purified genomic DNA or genomic DNA generated using whole genome amplification (WGA). Fluorescent signal output was measured in real time and as an end point. CONCLUSIONS: Combining the optimal elements for ligation/ERCA genotyping has resulted in a highly accurate single tube assay for genotyping directly from genomic DNA samples. Accuracy exceeded 99 % for four probe sets targeting clinically relevant mutations. No genotypes were called incorrectly using either genomic DNA or whole genome amplified sample.

Escherichia coli RecX inhibits RecA recombinase and coprotease activities in vitro and in vivo.

In Escherichia coli the RecA protein plays a pivotal role in homologous recombination, DNA repair, and SOS repair and mutagenesis. A gene designated recX (or oraA) is present directly downstream of recA in E. coli; however, the function of RecX is unknown. In this work we demonstrated interaction of RecX and RecA in a yeast two-hybrid assay. In vitro, substoichiometric amounts of RecX strongly inhibited both RecA-mediated DNA strand exchange and RecA ATPase activity. In vivo, we showed that recX is under control of the LexA repressor and is up-regulated in response to DNA damage. A loss-of-function mutation in recX resulted in decreased resistance to UV irradiation; however, overexpression of RecX in trans resulted in a greater decrease in UV resistance. Overexpression of RecX inhibited induction of two din (damage-inducible) genes and cleavage of the UmuD and LexA repressor proteins; however, recX inactivation had no effect on any of these processes. Cells overexpressing RecX showed decreased levels of P1 transduction, whereas recX mutation had no effect on P1 transduction frequency. Our combined in vitro and in vivo data indicate that RecX can inhibit both RecA recombinase and coprotease activities.