Sequential development of several RT-qPCR tests using LNA nucleotides and dual probe technology to differentiate SARS-CoV-2 from influenza A and B.

Sensitive and accurate RT-qPCR tests are the primary diagnostic tools to identify SARS-CoV-2-infected patients. While many SARS-CoV-2 RT-qPCR tests are available, there are significant differences in test sensitivity, workflow (e.g. hands-on-time), gene targets and other functionalities that users must consider. Several publicly available protocols shared by reference labs and public health authorities provide useful tools for SARS-CoV-2 diagnosis, but many have shortcomings related to sensitivity and laborious workflows. Here, we describe a series of SARS-CoV-2 RT-qPCR tests that are originally based on the protocol targeting regions of the RNA-dependent RNA polymerase (RdRp) and envelope (E) coding genes developed by the Charité Berlin. We redesigned the primers/probes, utilized locked nucleic acid nucleotides, incorporated dual probe technology and conducted extensive optimizations of reaction conditions to enhance the sensitivity and specificity of these tests. By incorporating an RNase P internal control and developing multiplexed assays for distinguishing SARS-CoV-2 and influenza A and B, we streamlined the workflow to provide quicker results and reduced consumable costs. Some of these tests use modified enzymes enabling the formulation of a room temperature-stable master mix and lyophilized positive control, thus increasing the functionality of the test and eliminating cold chain shipping and storage. Moreover, a rapid, RNA extraction-free version enables high sensitivity detection of SARS-CoV-2 in about an hour using minimally invasive, self-collected gargle samples. These RT-qPCR assays can easily be implemented in any diagnostic laboratory and can provide a powerful tool to detect SARS-CoV-2 and the most common seasonal influenzas during the vaccination phase of the pandemic.

Plasma microRNA Interindividual Variability in Healthy Individuals, Pregnant Women, and an Individual with a Stably Altered Neuroendocrine Phenotype.

BACKGROUND: Extracellular RNAs (exRNAs) in biofluids are amenable to quantitative analysis and proposed as noninvasive biomarkers for monitoring organ function. Cell-lineage-specific microRNAs (miRNAs) are present in plasma as soluble ribonucleoproteins or enclosed in exRNA carriers and transported through the vasculature. However, more extensive studies of healthy individuals are needed to gain insights into the variability of plasma miRNA abundance and composition. METHODS: The exRNA composition of platelet-depleted plasma collected twice from 236 healthy individuals was characterized by small RNA sequencing. Plasma of pregnant women featuring dramatically increased placental miRNAs and samples from subject P12 with noticeably increased epithelial- and neuroendocrine-origin miRNAs were included for comparison. The miRNA content of 10 000g and 100 000g pellet fractions of plasma generated by ultracentrifugation was also determined. Data analysis methods included Pearson correlation, differential gene expression, and unsupervised clustering. RESULTS: The abundance changes for more variable miRNAs in plasma of normal individuals correlated between coexpressed cell-lineage-specific miRNAs of the liver, neuroendocrine organs, epithelial cells, and muscle. ExRNA of pellet fractions contained <2% of total plasma miRNA with modest enrichment of lineage-specific and variable miRNAs compared to supernatant. The abundance fold changes of miRNAs observed in pregnancy and P12 compared to normal exceeded interquartile variability of healthy individuals. The neuroendocrine miRNA signature of P12 persisted for more than 4 years and was absent in other individuals. CONCLUSIONS: This study defines the framework and effect size for screening of extensive plasma collections for miRNA phenotypes and biomarker discovery.

Open Problems in Extracellular RNA Data Analysis: Insights From an ERCC Online Workshop.

We now know RNA can survive the harsh environment of biofluids when encapsulated in vesicles or by associating with lipoproteins or RNA binding proteins. These extracellular RNA (exRNA) play a role in intercellular signaling, serve as biomarkers of disease, and form the basis of new strategies for disease treatment. The Extracellular RNA Communication Consortium (ERCC) hosted a two-day online workshop (April 19-20, 2021) on the unique challenges of exRNA data analysis. The goal was to foster an open dialog about best practices and discuss open problems in the field, focusing initially on small exRNA sequencing data. Video recordings of workshop presentations and discussions are available (https://exRNA.org/exRNAdata2021-videos/). There were three target audiences: experimentalists who generate exRNA sequencing data, computational and data scientists who work with those groups to analyze their data, and experimental and data scientists new to the field. Here we summarize issues explored during the workshop, including progress on an effort to develop an exRNA data analysis challenge to engage the community in solving some of these open problems.

Characterizing and classifying neuroendocrine neoplasms through microRNA sequencing and data mining.

Neuroendocrine neoplasms (NENs) are clinically diverse and incompletely characterized cancers that are challenging to classify. MicroRNAs (miRNAs) are small regulatory RNAs that can be used to classify cancers. Recently, a morphology-based classification framework for evaluating NENs from different anatomical sites was proposed by experts, with the requirement of improved molecular data integration. Here, we compiled 378 miRNA expression profiles to examine NEN classification through comprehensive miRNA profiling and data mining. Following data preprocessing, our final study cohort included 221 NEN and 114 non-NEN samples, representing 15 NEN pathological types and 5 site-matched non-NEN control groups. Unsupervised hierarchical clustering of miRNA expression profiles clearly separated NENs from non-NENs. Comparative analyses showed that miR-375 and miR-7 expression is substantially higher in NEN cases than non-NEN controls. Correlation analyses showed that NENs from diverse anatomical sites have convergent miRNA expression programs, likely reflecting morphological and functional similarities. Using machine learning approaches, we identified 17 miRNAs to discriminate 15 NEN pathological types and subsequently constructed a multilayer classifier, correctly identifying 217 (98%) of 221 samples and overturning one histological diagnosis. Through our research, we have identified common and type-specific miRNA tissue markers and constructed an accurate miRNA-based classifier, advancing our understanding of NEN diversity.

A single-cell survey of the human first-trimester placenta and decidua.

The placenta and decidua interact dynamically to enable embryonic and fetal development. Here, we report single-cell RNA sequencing of 14,341 and 6754 cells from first-trimester human placental villous and decidual tissues, respectively. Bioinformatic analysis identified major cell types, many known and some subtypes previously unknown in placental villi and decidual context. Further detailed analysis revealed proliferating subpopulations, enrichment of cell type-specific transcription factors, and putative intercellular communication in the fetomaternal microenvironment. This study provides a blueprint to further the understanding of the roles of these cells in the placenta and decidua for maintenance of early gestation as well as pathogenesis in pregnancy-related disorders.

Human plasma and serum extracellular small RNA reference profiles and their clinical utility.

Circulating extracellular RNAs (exRNAs) have the potential to serve as biomarkers for a wide range of medical conditions. However, limitations in existing exRNA isolation methods and a lack of knowledge on parameters affecting exRNA variability in human samples may hinder their successful discovery and clinical implementation. Using combinations of denaturants, reducing agents, proteolysis, and revised organic extraction, we developed an automated, high-throughput approach for recovery of exRNAs and exDNA from the same biofluid sample. We applied this method to characterize exRNAs from 312 plasma and serum samples collected from 13 healthy volunteers at 12 time points over a 2-month period. Small RNA cDNA library sequencing identified nearly twofold increased epithelial-, muscle-, and neuroendocrine-cell-specific miRNAs in females, while fasting and hormonal cycle showed little effect. External standardization helped to detect quantitative differences in erythrocyte and platelet-specific miRNA contributions and in miRNA concentrations between biofluids. It also helped to identify a study participant with a unique exRNA phenotype featuring a miRNA signature of up to 20-fold elevated endocrine-cell-specific miRNAs and twofold elevated total miRNA concentrations stable for over 1 year. Collectively, these results demonstrate an efficient and quantitative method to discern exRNA phenotypes and suggest that plasma and serum RNA profiles are stable over months and can be routinely monitored in long-term clinical studies.

Urine cell-free microRNA as biomarkers for transitional cell carcinoma.

OBJECTIVE: MicroRNA (miRNA) are short nucleotide strands with a regulatory function in the cell. Several miRNAs have been shown to be useful as biomarkers for different neoplasms. The aim of this project was to assess whether levels of miRNA in cell free urine could be used as a biomarker in transitional cell carcinoma (TCC). RESULTS: cDNA libraries were produced based on small RNAs in urine samples of fourteen TCC patients and twenty healthy volunteers. Resulting reads were deep sequenced on Illumina HiSeq sequencer with the intent of characterizing cell free urine miRNA profiles. A statistically significant difference was found for a single miRNA; miR-210 was > sixfold higher in the TCC group compared to the control group. Furthermore, we were able to produce a diagnostic score by summing of standardized levels of overexpressed miRNA. This score was considerably higher in TCC patients with a sensitivity of 0.93, specificity of 0.76 and negative predictive value > 0.97.

Cell and Microvesicle Urine microRNA Deep Sequencing Profiles from Healthy Individuals: Observations with Potential Impact on Biomarker Studies.

BACKGROUND: Urine is a potential source of biomarkers for diseases of the kidneys and urinary tract. RNA, including microRNA, is present in the urine enclosed in detached cells or in extracellular vesicles (EVs) or bound and protected by extracellular proteins. Detection of cell- and disease-specific microRNA in urine may aid early diagnosis of organ-specific pathology. In this study, we applied barcoded deep sequencing to profile microRNAs in urine of healthy volunteers, and characterized the effects of sex, urine fraction (cells vs. EVs) and repeated voids by the same individuals. RESULTS: Compared to urine-cell-derived small RNA libraries, urine-EV-derived libraries were relatively enriched with miRNA, and accordingly had lesser content of other small RNA such as rRNA, tRNA and sn/snoRNA. Unsupervised clustering of specimens in relation to miRNA expression levels showed prominent bundling by specimen type (urine cells or EVs) and by sex, as well as a tendency of repeated (first and second void) samples to neighbor closely. Likewise, miRNA profile correlations between void repeats, as well as fraction counterparts (cells and EVs from the same specimen) were distinctly higher than correlations between miRNA profiles overall. Differential miRNA expression by sex was similar in cells and EVs. CONCLUSIONS: miRNA profiling of both urine EVs and sediment cells can convey biologically important differences between individuals. However, to be useful as urine biomarkers, careful consideration is needed for biofluid fractionation and sex-specific analysis, while the time of voiding appears to be less important.

Identification of the RNA recognition element of the RBPMS family of RNA-binding proteins and their transcriptome-wide mRNA targets.

Recent studies implicated the RNA-binding protein with multiple splicing (RBPMS) family of proteins in oocyte, retinal ganglion cell, heart, and gastrointestinal smooth muscle development. These RNA-binding proteins contain a single RNA recognition motif (RRM), and their targets and molecular function have not yet been identified. We defined transcriptome-wide RNA targets using photoactivatable-ribonucleoside-enhanced crosslinking and immunoprecipitation (PAR-CLIP) in HEK293 cells, revealing exonic mature and intronic pre-mRNA binding sites, in agreement with the nuclear and cytoplasmic localization of the proteins. Computational and biochemical approaches defined the RNA recognition element (RRE) as a tandem CAC trinucleotide motif separated by a variable spacer region. Similar to other mRNA-binding proteins, RBPMS family of proteins relocalized to cytoplasmic stress granules under oxidative stress conditions suggestive of a support function for mRNA localization in large and/or multinucleated cells where it is preferentially expressed.

Identification of mRNAs bound and regulated by human LIN28 proteins and molecular requirements for RNA recognition.

Human LIN28A and LIN28B are RNA-binding proteins (RBPs) conserved in animals with important roles during development and stem cell reprogramming. We used Photoactivatable-Ribonucleoside-Enhanced Crosslinking and Immunoprecipitation (PAR-CLIP) in HEK293 cells and identified a largely overlapping set of ∼3000 mRNAs at ∼9500 sites located in the 3' UTR and CDS. In vitro and in vivo, LIN28 preferentially bound single-stranded RNA containing a uridine-rich element and one or more flanking guanosines and appeared to be able to disrupt base-pairing to access these elements when embedded in predicted secondary structure. In HEK293 cells, LIN28 protein binding mildly stabilized target mRNAs and increased protein abundance. The top targets were its own mRNAs and those of other RBPs and cell cycle regulators. Alteration of LIN28 protein levels also negatively regulated the abundance of some but not all let-7 miRNA family members, indicating sequence-specific binding of let-7 precursors to LIN28 proteins and competition with cytoplasmic miRNA biogenesis factors.

RNA single strands bind to a conserved surface of the major cold shock protein in crystals and solution.

Bacterial cold shock proteins (CSPs) regulate the cellular response to temperature downshift. Their general principle of function involves RNA chaperoning and transcriptional antitermination. Here we present two crystal structures of cold shock protein B from Bacillus subtilis (Bs-CspB) in complex with either a hexanucleotide (5'-UUUUUU-3') or heptanucleotide (5'-GUCUUUA-3') single-stranded RNA (ssRNA). Hydrogen bonds and stacking interactions between RNA bases and aromatic sidechains characterize individual binding subsites. Additional binding subsites which are not occupied by the ligand in the crystal structure were revealed by NMR spectroscopy in solution on Bs-CspB·RNA complexes. Binding studies demonstrate that Bs-CspB associates with ssDNA as well as ssRNA with moderate sequence specificity. Varying affinities of oligonucleotides are reflected mainly in changes of the dissociation rates. The generally lower binding affinity of ssRNA compared to its ssDNA analog is attributed solely to the substitution of thymine by uracil bases in RNA.

The coxsackievirus-adenovirus receptor reveals complex homophilic and heterophilic interactions on neural cells.

The coxsackievirus-adenovirus receptor (CAR) is a member of the Ig superfamily strongly expressed in the developing nervous system. Our histological investigations during development reveal an initial uniform distribution of CAR on all neural cells with a concentration on membranes that face the margins of the nervous system (e.g., the basal laminae and the ventricular side). At more advanced stages, CAR becomes downregulated and restricted to specific regions including areas rich in axonal and dendritic surfaces. To study the function of CAR on neural cells, we used the fiber knob of the adenovirus, extracellular CAR domains, blocking antibodies to CAR, as well as CAR-deficient neural cells. Blocking antibodies were found to inhibit neurite extension in retina organ and retinal explant cultures, whereas the application of the recombinant fiber knob of the adenovirus subtype Ad2 or extracellular CAR domains promoted neurite extension and adhesion to extracellular matrices. We observed a promiscuous interaction of CAR with extracellular matrix glycoproteins, which was deduced from analytical ultracentrifugation experiments, affinity chromatography, and adhesion assays. The membrane proximal Ig domain of CAR, termed D2, was found to bind to a fibronectin fragment, including the heparin-binding domain 2, which promotes neurite extension of wild type, but not of CAR-deficient neural cells. In contrast to heterophilic interactions, homophilic association of CAR involves both Ig domains, as was revealed by ultracentrifugation, chemical cross-linking, and adhesion studies. The results of these functional and binding studies are correlated to a U-shaped homodimer of the complete extracellular domains of CAR detected by x-ray crystallography.

Dimer formation of a stabilized Gbeta1 variant: a structural and energetic analysis.

In previous work, a strongly stabilized variant of the beta1 domain of streptococcal protein G (Gbeta1) was obtained by an in vitro selection method. This variant, termed Gbeta1-M2, contains the four substitutions E15V, T16L, T18I, and N37L. Here we elucidated the molecular basis of the observed strong stabilizations. The contributions of these four residues were analyzed individually and in various combinations, additional selections with focused Gbeta1 gene libraries were performed, and the crystal structure of Gbeta1-M2 was determined. All single substitutions (E15V, T16L, T18I, and N37L) stabilize wild-type Gbeta1 by contributions of between 1.6 and 6.0 kJ mol(-1) (at 70 degrees C). Hydrophobic residues at positions 16 and 37 provide the major contribution to stabilization by enlarging the hydrophobic core of Gbeta1. They also increase the tendency to form dimers, as shown by dependence on the concentration of apparent molecular mass in analytical ultracentrifugation, by concentration-dependent stability, and by a strongly increased van't Hoff enthalpy of unfolding. The 0.88-A crystal structure of Gbeta1-M2 and NMR measurements in solution provide the explanation for the observed dimer formation. It involves a head-to-head arrangement of two Gbeta1-M2 molecules via six intermolecular hydrogen bonds between the two beta strands 2 and 2' and an adjacent self-complementary hydrophobic surface area, which is created by the T16L and N37L substitutions and a large 120 degrees rotation of the Tyr33 side chain. This removal of hydrophilic groups and the malleability of the created hydrophobic surface provide the basis for the dimer formation of stabilized Gbeta1 variants.

Optimization of the gbeta1 domain by computational design and by in vitro evolution: structural and energetic basis of stabilization.

Computational design and in vitro evolution are major strategies for stabilizing proteins. For the four critical positions 16, 18, 25, and 29 of the B domain of the streptococcal protein G (Gbeta1), they identified the same optimal residues at positions 16 and 25, but not at 18 and 29. Here we analyzed the energetic contributions of the residues from these two approaches by single and double mutant analyses and determined crystal structures for a variant from the calculation (I16/L18/E25/K29) and from the selection (I16/I18/E25/F29). The structural analysis explains the observed differences in stabilization. Residues 16, 18, and 29 line an invagination, which results from a packing defect between the helix and the beta-sheet of Gbeta1. In all stabilized variants, residues with larger side-chains occur at these positions and packing is improved. In the selected variant, packing is better optimized than in the computed variant. Such differences in side-chain packing strongly affect stability but are difficult to evaluate by computation.

Optimized variants of the cold shock protein from in vitro selection: structural basis of their high thermostability.

The bacterial cold shock proteins (Csp) are widely used as models for the experimental and computational analysis of protein stability. In a previous study, in vitro evolution was employed to identify strongly stabilizing mutations in Bs-CspB from Bacillus subtilis. The best variant found by this approach contained the mutations M1R, E3K and K65I, which raised the midpoint of thermal unfolding of Bs-CspB from 53.8 degrees C to 83.7 degrees C, and increased the Gibbs free energy of stabilization by 20.9 kJ mol(-1). Another selected variant with the two mutations A46K and S48R was stabilized by 11.1 kJ mol(-1). To elucidate the molecular basis of these stabilizations, we determined the crystal structures of these two Bs-CspB variants. The mutated residues are generally well ordered and provide additional stabilizing interactions, such as charge interactions, additional hydrogen bonds and improved side-chain packing. Several mutations improve the electrostatic interactions, either by the removal of unfavorable charges (E3K) or by compensating their destabilizing interactions (A46K, S48R). The stabilizing mutations are clustered at a contiguous surface area of Bs-CspB, which apparently is critically important for the stability of the beta-barrel structure but not well optimized in the wild-type protein.

Sequence specificity of single-stranded DNA-binding proteins: a novel DNA microarray approach.

We have developed a novel DNA microarray-based approach for identification of the sequence-specificity of single-stranded nucleic-acid-binding proteins (SNABPs). For verification, we have shown that the major cold shock protein (CspB) from Bacillus subtilis binds with high affinity to pyrimidine-rich sequences, with a binding preference for the consensus sequence, 5'-GTCTTTG/T-3'. The sequence was modelled onto the known structure of CspB and a cytosine-binding pocket was identified, which explains the strong preference for a cytosine base at position 3. This microarray method offers a rapid high-throughput approach for determining the specificity and strength of ss DNA-protein interactions. Further screening of this newly emerging family of transcription factors will help provide an insight into their cellular function.

Common mode of DNA binding to cold shock domains. Crystal structure of hexathymidine bound to the domain-swapped form of a major cold shock protein from Bacillus caldolyticus.

Bacterial cold shock proteins (CSPs) regulate cellular adaptation to cold stress. Functions ascribed to CSP include roles as RNA chaperones and in transcription antitermination. We present the crystal structure of the Bacillus caldolyticus CSP (Bc-Csp) in complex with hexathymidine (dT(6)) at a resolution of 1.29 A. Bound to dT(6), crystalline Bc-Csp forms a domain-swapped dimer in which beta strands 1-3 associate with strands 4 and 5 from the other subunit to form a closed beta barrel and vice versa. The globular units of dimeric Bc-Csp closely resemble the well-known structure of monomeric CSP. Structural reorganization from the monomer to the domain-swapped dimer involves a strictly localized change in the peptide bond linking Glu36 and Gly37 of Bc-Csp. Similar structural reorganizations have not been found in any other CSP or oligonucleotide/oligosaccharide-binding fold structures. Each dT(6) ligand is bound to one globular unit of Bc-Csp via an amphipathic protein surface. Individual binding subsites interact with the DNA bases through stacking and hydrogen bonding. The sugar-phosphate backbone remains solvent exposed. Based on crystallographic and biochemical studies of deoxyoligonucleotide binding to CSP, we suggest a common mode of binding of single-stranded heptanucleotide motifs to proteins containing cold shock domains, including the eukaryotic Y-box factors.

Recognition of T-rich single-stranded DNA by the cold shock protein Bs-CspB in solution.

Cold shock proteins (CSP) belong to the family of single-stranded nucleic acid binding proteins with OB-fold. CSP are believed to function as 'RNA chaperones' and during anti-termination. We determined the solution structure of Bs-CspB bound to the single-stranded DNA (ssDNA) fragment heptathymidine (dT7) by NMR spectroscopy. Bs-CspB reveals an almost invariant conformation when bound to dT7 with only minor reorientations in loop beta1-beta2 and beta3-beta4 and of few aromatic side chains involved in base stacking. Binding studies of protein variants and mutated ssDNA demonstrated that Bs-CspB associates with ssDNA at almost diffusion controlled rates and low sequence specificity consistent with its biological function. A variation of the ssDNA affinity is accomplished solely by changes of the dissociation rate. 15N NMR relaxation and H/D exchange experiments revealed that binding of dT7 increases the stability of Bs-CspB and reduces the sub-nanosecond dynamics of the entire protein and especially of loop beta3-beta4.

T-rich DNA single strands bind to a preformed site on the bacterial cold shock protein Bs-CspB.

Bacterial cold shock proteins (CSPs) are involved in cellular adaptation to cold stress. They bind to single-stranded nucleic acids with a KD value in the micro- to nanomolar range. Here we present the structure of the Bacillus subtilis CspB (Bs-CspB) in complex with hexathymidine (dT6) at a resolution of 1.78 A. Bs-CspB binds to dT6 with nanomolar affinity via an amphipathic interface on the protein surface. Individual binding subsites interact with single nucleobases through stacking interactions and hydrogen bonding. The sugar-phosphate backbone and the methyl groups of the thymine nucleobases remain solvent exposed and are not contacted by protein groups. Fluorescence titration experiments monitoring the binding of oligopyrimidines to Bs-CspB reveal binding preferences at individual subsites and allow the design of an optimised heptapyrimidine ligand, which is bound with sub-nanomolar affinity. This study reveals the stoichiometry and sequence determinants of the binding of single-stranded nucleic acids to a preformed site on Bs-CspB and thus provides the structural basis of the RNA chaperone and transcription antitermination activities of the CSP.