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.

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.

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 F2-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 F2-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 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.

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. 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. 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. 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.

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.

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.

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.

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.

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.

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.

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).

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.

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.

Plants, as any other organisms, possess evolutionary old mechanisms to cope with the various stresses they are exposed to day by day. The management of stresses and their consequences requires substantial energy, which is frequently subtracted from biomass (in crops: yield). Therefore, a deeper understanding of stress biology has been, is, and will be of paramount importance for plant breeding. One goal of plant stress research centers around the transcriptome, the entirety of transcripts from expressed genes, and aims at identifying major genes in the stress management of the inflicted plant. The development of appropriate technologies to quantitatively study the transcriptomes (indeed the various sub-transcriptomes) in stressed plants and to extract biological meaning from the massive data will be demonstrated here. In particular, reduced complexity sequencing techniques such as deepSuperSAGE and MACE (massive analysis of cDNA ends) and their potential in stress biology are portrayed.

A 454-FLX low-coverage sequencing approach was used to assemble the mitochondrial genome of Radix balthica. The mtDNA sequence is 13,993nt long and contains 37 genes (13 protein coding genes, two rRNAs and 22 tRNAs). Four genes, the 12S RNA and seven tRNAs are transcribed in reverse order. The sequence is AT rich (71.3%), similar to other basommatophoran species. Comparison with the most closely related mt genomes available (Biomphalaria glabrata and Biomphalaria tenagophila) revealed identical gene orders except for five tRNAs. Next generation sequencing proved to be a fast and easy method for sequencing an entire mitochondrial genome.