Optimized qRT-PCR Approach for the Detection of Intra- and Extra-Cellular SARS-CoV-2 RNAs.

The novel coronavirus SARS-CoV-2 is the causative agent of the acute respiratory disease COVID-19, which has become a global concern due to its rapid spread. Meanwhile, increased demand for testing has led to a shortage of reagents and supplies and compromised the performance of diagnostic laboratories in many countries. Both the World Health Organization (WHO) and the Center for Disease Control and Prevention (CDC) recommend multi-step RT-PCR assays using multiple primer and probe pairs, which might complicate the interpretation of the test results, especially for borderline cases. In this study, we describe an alternative RT-PCR approach for the detection of SARS-CoV-2 RNA that can be used for the probe-based detection of clinical isolates in diagnostics as well as in research labs using a low-cost SYBR green method. For the evaluation, we used samples from patients with confirmed SARS-CoV-2 infections and performed RT-PCR assays along with successive dilutions of RNA standards to determine the limit of detection. We identified an M-gene binding primer and probe pair highly suitable for the quantitative detection of SARS-CoV-2 RNA for diagnostic and research purposes.

Deep Super-SAGE transcriptomic analysis of cold acclimation in lentil (Lens culinaris Medik.).

BACKGROUND: Frost is one of the main abiotic stresses limiting plant distribution and crop production. To cope with the stress, plants evolved adaptations known as cold acclimation or chilling tolerance to maximize frost tolerance. Cold acclimation is a progressive acquisition of freezing tolerance by plants subjected to low non-freezing temperatures which subsequently allows them to survive exposure to frost. Lentil is a cool season grain legume that is challenged by winter frost in some areas of its cultivation. RESULTS: To better understand the genetic base of frost tolerance differential gene expression in response to cold acclimation was investigated. Recombinant inbred lines (RILs) from the cross Precoz x WA8649041 were first classified as cold tolerant or cold susceptible according to their response to temperatures between -3 to -15 °C. Then, RILs from both extremes of the response curve were cold acclimated and the leaf transcriptomes of two bulks each of eight frost tolerant and seven cold susceptible RILs were investigated by Deep Super-SAGE transcriptome profiling. Thus, four RNA bulks were analysed: the acclimated susceptible, the acclimated tolerant and the respective controls (non-acclimated susceptible and non-acclimated tolerant). Approximately 16.5 million 26 nucleotide long Super-SAGE tags were sequenced in the four sets (between ~3 and 5.4 millions). In total, 133,077 different unitags, each representing a particular transcript isoform, were identified in these four sets. Tags which showed a significantly different abundance in any of the bulks (fold change ≥4.0 and a significant p-value <0.001) were selected and used to identify the corresponding lentil gene sequence. Three hundred of such lentil sequences were identified. Most of their known homologs coded for glycine-rich, cold and drought-regulated proteins, dormancy-associated proteins, proline-rich proteins (PRPs) and other membrane proteins. These were generally but not exclusively over-expressed in the acclimated tolerant lines. CONCLUSIONS: This set of candidate genes implicated in the response to frost in lentil represents an useful base for deeper and more detailed investigations into this important agronomic trait in future.

Dual 3'Seq using deepSuperSAGE uncovers transcriptomes of interacting Salmonella enterica Typhimurium and human host cells.

BACKGROUND: The interaction of eukaryotic host and prokaryotic pathogen cells is linked to specific changes in the cellular proteome, and consequently to infection-related gene expression patterns of the involved cells. To simultaneously assess the transcriptomes of both organisms during their interaction we developed dual 3'Seq, a tag-based sequencing protocol that allows for exact quantification of differentially expressed transcripts in interacting pro- and eukaryotic cells without prior fixation or physical disruption of the interaction. RESULTS: Human epithelial cells were infected with Salmonella enterica Typhimurium as a model system for invasion of the intestinal epithelium, and the transcriptional response of the infected host cells together with the differential expression of invading and intracellular pathogen cells was determined by dual 3'Seq coupled with the next-generation sequencing-based transcriptome profiling technique deepSuperSAGE (deep Serial Analysis of Gene Expression). Annotation to reference transcriptomes comprising the operon structure of the employed S. enterica Typhimurium strain allowed for in silico separation of the interacting cells including quantification of polycistronic RNAs. Eighty-nine percent of the known loci are found to be transcribed in prokaryotic cells prior or subsequent to infection of the host, while 75% of all protein-coding loci are represented in the polyadenylated transcriptomes of human host cells. CONCLUSIONS: Dual 3'Seq was alternatively coupled to MACE (Massive Analysis of cDNA ends) to assess the advantages and drawbacks of a library preparation procedure that allows for sequencing of longer fragments. Additionally, the identified expression patterns of both organisms were validated by qRT-PCR using three independent biological replicates, which confirmed that RELB along with NFKB1 and NFKB2 are involved in the initial immune response of epithelial cells after infection with S. enterica Typhimurium.

Identification of novel small ncRNAs in pollen of tomato.

BACKGROUND: The unprecedented role of sncRNAs in the regulation of pollen biogenesis on both transcriptional and epigenetic levels has been experimentally proven. However, little is known about their global regulation, especially under stress conditions. We used tomato pollen in order to identify pollen stage-specific sncRNAs and their target mRNAs. We further deployed elevated temperatures to discern stress responsive sncRNAs. For this purpose high throughput sncRNA-sequencing as well as Massive Analysis of cDNA Ends (MACE) were performed for three-replicated sncRNAs libraries derived from tomato tetrad, post-meiotic, and mature pollen under control and heat stress conditions. RESULTS: Using the omiRas analysis pipeline we identified known and predicted novel miRNAs as well as sncRNAs from other classes, responsive or not to heat. Differential expression analysis revealed that post-meiotic and mature pollen react most strongly by regulation of the expression of coding and non-coding genomic regions in response to heat. To gain insight to the function of these miRNAs, we predicted targets and annotated them to Gene Ontology terms. This approach revealed that most of them belong to protein binding, transcription, and Serine/Threonine kinase activity GO categories. Beside miRNAs, we observed differential expression of both tRNAs and snoRNAs in tetrad, post-meiotic, and mature pollen when comparing normal and heat stress conditions. CONCLUSIONS: Thus, we describe a global spectrum of sncRNAs expressed in pollen as well as unveiled those which are regulated at specific time-points during pollen biogenesis. We integrated the small RNAs into the regulatory network of tomato heat stress response in pollen.

Next-generation sequencing reveals novel differentially regulated mRNAs, lncRNAs, miRNAs, sdRNAs and a piRNA in pancreatic cancer.

BACKGROUND: Previous studies identified microRNAs (miRNAs) and messenger RNAs with significantly different expression between normal pancreas and pancreatic cancer (PDAC) tissues. Due to technological limitations of microarrays and real-time PCR systems these studies focused on a fixed set of targets. Expression of other RNA classes such as long intergenic non-coding RNAs or sno-derived RNAs has rarely been examined in pancreatic cancer. Here, we analysed the coding and non-coding transcriptome of six PDAC and five control tissues using next-generation sequencing. RESULTS: Besides the confirmation of several deregulated mRNAs and miRNAs, miRNAs without previous implication in PDAC were detected: miR-802, miR-2114 or miR-561. SnoRNA-derived RNAs (e.g. sno-HBII-296B) and piR-017061, a piwi-interacting RNA, were found to be differentially expressed between PDAC and control tissues. In silico target analysis of miR-802 revealed potential binding sites in the 3' UTR of TCF4, encoding a transcription factor that controls Wnt signalling genes. Overexpression of miR-802 in MiaPaCa pancreatic cancer cells reduced TCF4 protein levels. Using Massive Analysis of cDNA Ends (MACE) we identified differential expression of 43 lincRNAs, long intergenic non-coding RNAs, e.g. LINC00261 and LINC00152 as well as several natural antisense transcripts like HNF1A-AS1 and AFAP1-AS1. Differential expression was confirmed by qPCR on the mRNA/miRNA/lincRNA level and by immunohistochemistry on the protein level. CONCLUSIONS: Here, we report a novel lncRNA, sncRNA and mRNA signature of PDAC. In silico prediction of ncRNA targets allowed for assigning potential functions to differentially regulated RNAs.

Correction: The Swine Plasma Metabolome Chronicles "Many Days" Biological Timing and Functions Linked to Growth.

[This corrects the article DOI: 10.1371/journal.pone.0145919.].

Two variants of the assembly factor Surf1 target specific terminal oxidases in Paracoccus denitrificans.

Biogenesis of cytochrome c oxidase (COX) relies on a large number of assembly proteins, one of them being Surf1. In humans, the loss of Surf1 function is associated with Leigh syndrome, a fatal neurodegenerative disorder. In the soil bacterium Paracoccus denitrificans, homologous genes specifying Surf1 have been identified and located in two operons of terminal oxidases: surf1q is the last gene of the qox operon (coding for a ba(3)-type ubiquinol oxidase), and surf1c is found at the end of the cta operon (encoding subunits of the aa(3)-type cytochrome c oxidase). We introduced chromosomal single and double deletions for both surf1 genes, leading to significantly reduced oxidase activities in membrane. Our experiments on P. denitrificans surf1 single deletion strains show that both Surf1c and Surf1q are functional and act independently for the aa(3)-type cytochrome c oxidase and the ba(3)-type quinol oxidase, respectively. This is the first direct experimental evidence for the involvement of a Surf1 protein in the assembly of a quinol oxidase. Analyzing the heme content of purified cytochrome c oxidase, we conclude that Surf1, though not indispensable for oxidase assembly, is involved in an early step of cofactor insertion into subunit I.

Mapping-by-sequencing using NGS-based 3'-MACE-Seq reveals a new mutant allele of the essential nodulation gene Sym33 (IPD3) in pea (Pisum sativum L.).

Large collections of pea symbiotic mutants were accumulated in the 1990s, but the causal genes for a large portion of the mutations are still not identified due to the complexity of the task. We applied a Mapping-by-Sequencing approach including Bulk Segregant Analysis and Massive Analysis of cDNA Ends (MACE-Seq) sequencing technology for genetic mapping the Sym11 gene of pea which controls the formation of symbioses with both nodule bacteria and arbuscular-mycorrhizal fungi. For mapping we developed an F 2-population from the cross between pea line N24 carrying the mutant allele of sym11 and the wild type NGB1238 (=JI0073) line. Sequencing libraries were prepared from bulks of 20 plants with mutant and 12 with wild-type phenotype. MACE-Seq differential gene expression analysis between mutant-phenotype and wild-type-phenotype bulks revealed 2,235 genes, of which 514 (23%) were up-regulated and 1,721 (77%) were down-regulated in plant roots inoculated with rhizobia as a consequence of sym11 mutation. MACE-Seq also detected single nucleotide variants between bulks in 217 pea genes. Using a novel mathematical model we calculated the recombination frequency (RF) between the Sym11 gene and these 217 polymorphic genes. Six genes with the lowest RF were converted into CAPS or dCAPS markers and genetically mapped on the complete mapping population of 108 F 2-plants which confirmed their tight linkage to Sym11 and to each other. The Medicago truncatula Gaertn. (Mt) homologs of these genes are located in a distinct region of Mt chromosome 5, which corresponds to linkage group I of pea. Among 94 candidate genes from this region only one was down-regulated-the pea Sym33 homolog of the Mt IPD3 gene which is essential for nodulation. Sequencing of the Sym33 allele of the N24 (sym11) mutant revealed a single nucleotide deletion (c.C319del) in its third exon resulting in a codon shift in the open reading frame and premature translation termination. Thus, we identified a novel mutant allele sym33-4 most probably responsible for the mutant phenotype of the N24 (sym11) line, thereby demonstrating that mapping by MACE-Seq can be successfully used for genetic mapping of mutations and identification of candidate genes in pea.

Morphological and transcriptomic effects of endocrine modulators on the gonadal differentiation of chicken embryos: The case of tributyltin (TBT).

Morphological malformations induced by tributyltin (TBT) exposure during embryonic development have already been characterized in various taxonomic groups, but, nonetheless, the molecular processes underlying these changes remain obscure. The present study provides the first genome-wide screening for differentially expressed genes that are linked to morphological alterations of gonadal tissue from chicken embryos after exposure to TBT. We applied a single injection of TBT (between 0.5 and 30 pg as Sn/g egg) into incubated fertile eggs to simulate maternal transfer of the endocrine disruptive compound. Methyltestosterone (MT) served as a positive control (30 pg/g egg). After 19 days of incubation, structural features of the gonads as well as genome-wide gene expression profiles were assessed simultaneously. TBT induced significant morphological and histological malformations of gonadal tissue from female embryos that show a virilization of the ovaries. This phenotypical virilization was mirrored by altered expression profiles of sex-dependent genes. Among these are several transcription and growth factors (e.g. FGF12, CTCF, NFIB), whose altered expression might serve as a set of markers for early identification of endocrine active chemicals that affect embryonic development by transcriptome profiling without the need of elaborate histological analyses.

Increased expression of hypoxia-inducible factor-1 alpha and its impact on transcriptional changes and prognosis in malignant tumours of the ocular adnexa.

PURPOSE: To investigate the expression profile of the hypoxia-inducible transcription factor-1α (HIF-1α) and its downstream targets in malignancies of the ocular adnexa and to determine its relevance as a prognostic factor for clinical outcome. METHODS: We included 49 subjects with malignant tumours (25 squamous cell carcinomas (SCC), 15 non-Hodgkin lymphomas, 9 melanomas) and 30 patients with benign tumours of the ocular adnexa (13 papillomas, 7 reactive lymphoid hyperplasias (RLHs) and 10 nevi) as controls. We quantified HIF-1α protein expression by immunohistochemistry and assessed the association between HIF-1α and clinical outcome via Kaplan-Meier analysis. Furthermore, we assessed the expression of HIF-1α downstream factors by transcriptional sequencing using the MACE (massive analysis of cDNA ends) technology. RESULTS: SCCs revealed a strong HIF-1α expression in 61% of tumour cells in comparison with only 22% in papillomas (p < 0.0001). In contrast, malignant melanomas and lymphomas revealed a similar HIF-1α expression compared with nevi and RLHs. Transcriptional sequencing and Gene Ontology Cluster analysis demonstrated 37 hypoxia-associated factors, including HIF-1α, VEGF, SFRP1 and LOXL2 that are significantly increased in SCC and may contribute to tumour proliferation, angiogenesis, and metastasis. Association analysis between HIF-1α immunoreactivity and clinical outcome revealed a trend towards an unfavourable prognosis in malignant tumours with increased HIF-1α expression. CONCLUSIONS: HIF-1α protein is increased in malignant tumours of the ocular adnexa, which is associated with an increase in multiple HIF-1α-downstream factors and a trend towards an unfavourable clinical outcome.

Precision medicine for hepatocelluar carcinoma using molecular pattern diagnostics: results from a preclinical pilot study.

The aim of this study was to design a road map for personalizing cancer therapy in hepatocellular carcinoma (HCC) by using molecular pattern diagnostics. As an exploratory study, we investigated molecular patterns of tissues of two tumors from individual HCC patients, which in previous experiments had shown contrasting reactions to the phase 2 transforming growth factor beta receptor 1 inhibitor galunisertib. Cancer-driving molecular patterns encompass - inter alias - altered transcription profiles and somatic mutations in coding regions differentiating tumors from their respective peritumoral tissues and from each other. Massive analysis of cDNA ends and all-exome sequencing demonstrate a highly divergent transcriptional and mutational landscape, respectively, for the two tumors, that offers potential explanations for the tumors contrasting responses to galunisertib. Molecular pattern diagnostics (MPDs) suggest alternative, individual-tumor-specific therapies, which in both cases deviate from the standard sorafenib treatment and from each other. Suggested personalized therapies use kinase inhibitors and immune-focused drugs as well as low-toxicity natural compounds identified using an advanced bioinformatics routine included in the MPD protocol. The MPD pipeline we describe here for the prediction of suitable drugs for treatment of two contrasting HCCs may serve as a blueprint for the design of therapies for various types of cancer.

MicroRNA Profiling in Aqueous Humor of Individual Human Eyes by Next-Generation Sequencing.

PURPOSE: Extracellular microRNAs (miRNAs) in aqueous humor were suggested to have a role in transcellular signaling and may serve as disease biomarkers. The authors adopted next-generation sequencing (NGS) techniques to further characterize the miRNA profile in single samples of 60 to 80 µL human aqueous humor. METHODS: Samples were obtained at the outset of cataract surgery in nine independent, otherwise healthy eyes. Four samples were used to extract RNA and generate sequencing libraries, followed by an adapter-driven amplification step, electrophoretic size selection, sequencing, and data analysis. Five samples were used for quantitative PCR (qPCR) validation of NGS results. Published NGS data on circulating miRNAs in blood were analyzed in comparison. RESULTS: One hundred fifty-eight miRNAs were consistently detected by NGS in all four samples; an additional 59 miRNAs were present in at least three samples. The aqueous humor miRNA profile shows some overlap with published NGS-derived inventories of circulating miRNAs in blood plasma with high prevalence of human miR-451a, -21, and -16. In contrast to blood, miR-184, -4448, -30a, -29a, -29c, -19a, -30d, -205, -24, -22, and -3074 were detected among the 20 most prevalent miRNAs in aqueous humor. Relative expression patterns of miR-451a, -202, and -144 suggested by NGS were confirmed by qPCR. CONCLUSIONS: Our data illustrate the feasibility of miRNA analysis by NGS in small individual aqueous humor samples. Intraocular cells as well as blood plasma contribute to the extracellular aqueous humor miRNome. The data suggest possible roles of miRNA in intraocular cell adhesion and signaling by TGF-ß and Wnt, which are important in intraocular pressure regulation and glaucoma.

The Swine Plasma Metabolome Chronicles "Many Days" Biological Timing and Functions Linked to Growth.

The paradigm of chronobiology is based almost wholly upon the daily biological clock, or circadian rhythm, which has been the focus of intense molecular, cellular, pharmacological, and behavioral, research. However, the circadian rhythm does not explain biological timings related to fundamental aspects of life history such as rates of tissue/organ/body size development and control of the timing of life stages such as gestation length, age at maturity, and lifespan. This suggests that another biological timing mechanism is at work. Here we focus on a "many days" (multidien) chronobiological period first observed as enigmatic recurring growth lines in developing mammalian tooth enamel that is strongly associate with all adult tissue, organ, and body masses as well as life history attributes such as gestation length, age at maturity, weaning, and lifespan, particularly among the well studied primates. Yet, knowledge of the biological factors regulating the patterning of mammalian life, such as the development of body size and life history structure, does not exist. To identify underlying molecular mechanisms we performed metabolome and genome analyses from blood plasma in domestic pigs. We show that blood plasma metabolites and small non-coding RNA (sncRNA) drawn from 33 domestic pigs over a two-week period strongly oscillate on a 5-day multidien rhythm, as does the pig enamel rhythm. Metabolomics and genomics pathway analyses actually reveal two 5-day rhythms, one related to growth in which biological functions include cell proliferation, apoptosis, and transcription regulation/protein synthesis, and another 5-day rhythm related to degradative pathways that follows three days later. Our results provide experimental confirmation of a 5-day multidien rhythm in the domestic pig linking the periodic growth of enamel with oscillations of the metabolome and genome. This association reveals a new class of chronobiological rhythm and a snapshot of the biological bases that regulate mammalian growth, body size, and life history.

Complete genome sequence of the gram-negative probiotic Escherichia coli strain Nissle 1917.

Escherichia coli strain Nissle 1917 (EcN) is the active principle of a probiotic preparation (trade name Mutaflor(®)) used for the treatment of patients with intestinal diseases such as ulcerative colitis and diarrhea. It has GRAS (generally recognized as save) status and has been shown to be a therapeutically effective drug (Sonnenborn and Schulze, 2009). The complete genomic DNA sequence will help in identifying genes and their products which are essential for the strains probiotic nature. Genbank/EMBL/DDBJ accession number: CP007799 (chromosome).

Gene expression of chicken gonads is sex- and side-specific.

In chicken, the left and right female gonads undergo a completely different program during development. To learn more about the molecular factors underlying side-specific development and to identify potential sex- and side-specific genes in developing gonads, we separately performed next-generation sequencing-based deepSuperSAGE transcription profiling from left and right, female and male gonads of 19-day-old chicken embryos. A total of 836 transcript variants were significantly differentially expressed (p < 10(-5)) between combined male and female gonads. Left-right comparison revealed 1,056 and 822 differentially (p < 10(-5)) expressed transcript variants for male and female gonads, respectively, of which 72 are side-specific in both sexes. At least some of these may represent key players for lateral development in birds. Additionally, several genes with laterally differential expression in the ovaries seem to determine female gonads for growth or regression, whereas right-left differences in testes are mostly limited to the differentially expressed genes present in both sexes. With a few exceptions, side-specific genes are not located on the sex chromosomes. The large differences in lateral gene expression in the ovaries in almost all metabolic pathways suggest that the regressing right gonad might have undergone a change of function during evolution.

Genome-based analysis of the transcriptome from mature chickpea root nodules.

Symbiotic nitrogen fixation (SNF) in root nodules of grain legumes such as chickpea is a highly complex process that drastically affects the gene expression patterns of both the prokaryotic as well as eukaryotic interacting cells. A successfully established symbiotic relationship requires mutual signaling mechanisms and a continuous adaptation of the metabolism of the involved cells to varying environmental conditions. Although some of these processes are well understood today many of the molecular mechanisms underlying SNF, especially in chickpea, remain unclear. Here, we reannotated our previously published transcriptome data generated by deepSuperSAGE (Serial Analysis of Gene Expression) to the recently published draft genome of chickpea to assess the root- and nodule-specific transcriptomes of the eukaryotic host cells. The identified gene expression patterns comprise up to 71 significantly differentially expressed genes and the expression of twenty of these was validated by quantitative real-time PCR with the tissues from five independent biological replicates. Many of the differentially expressed transcripts were found to encode proteins implicated in sugar metabolism, antioxidant defense as well as biotic and abiotic stress responses of the host cells, and some of them were already known to contribute to SNF in other legumes. The differentially expressed genes identified in this study represent candidates that can be used for further characterization of the complex molecular mechanisms underlying SNF in chickpea.

Introducing a novel human mtDNA mutation into the Paracoccus denitrificans COX I gene explains functional deficits in a patient.

We identified a novel mutation (S142F) in the human mtDNA CO I gene in a patient with a clinical phenotype resembling mitochondrial cardioencephalomyopathy. To substantiate pathogenicity, we modeled the identified mutation in the homologous gene in Paracoccus denitrificans and analyzed the biochemical consequences. We observed a deleterious effect on enzyme activity, with a lack of heme a3. Taking advantage of the extensive structural homology between the bacterial enzyme and the mammalian core complex, we conclude that the novel S142F mutation is disease-related. This approach can be used in other cases to support the pathogenicity of novel variants in the mitochondrial genome.