Potato is among one of the most important food crops, yet maintaining plant productivity in this drought-sensitive crop has become a challenge. Competition for scarce water resources and the continued effects of global warming exacerbate current constraints on crop production. While plants’ response to drought in above-ground tissues has been well documented, the regulatory cascades in developing tubers have been largely unexplored. Using the commercial Canadian cultivar ‘Vigor’, plants were subjected to a drought treatment under high-tunnels causing a 4 ℃ increase in canopy temperature when compared to the well-watered control. Tubers were sampled for RNAseq and metabolite analysis. Approximately 2600 genes and 3898 transcripts were differentially expressed by at least four-fold in drought-stressed potato tubers, with 75 % and 69 % being down-regulated respectively. A further 229 small RNAs were implicated in gene regulation during drought. The comparison of protein homologues between Solanum tuberosum L. and Arabidopsis thaliana L. indicates that downregulated genes are associated with phenylpropanoid, carotenoid, and patatin biosynthesis. This suggests that there may be nutritive implications to drought stress occurring during the potato tuber bulking phase in sensitive cultivars.

The New Zealand brushtail possum (Trichosurus vulpecula), introduced from Australia in the 1850s, is an invasive species that is widespread throughout New Zealand and forms the largest self-sustained reservoir of bovine tuberculosis (Mycobacterium bovis) in the wild. Conservation and agricultural authorities regularly apply a series of population control measures to suppress brushtail possum populations. The evolutionary consequence of more than half a century of intensive population control operations on the species genomic diversity and population structure is hindered by a paucity of available genomic resources for the species. This study is the first to characterise the functional content and diversity of brushtail possum liver and brain cerebral cortex transcriptomes. Raw sequences from hepatic cells and cerebral cortex were assembled into 58,001 and 64,735 transcripts in hepatic cells and cerebral cortex, respectively. Functional annotation and polymorphism assignment of the assembled transcripts demonstrated a considerable level of variation in the core metabolic pathways that represent potential targets for selection pressure exerted by chemical toxicants. This study suggests that the brushtail possum population in New Zealand harbours considerable levels of variation in metabolic pathways that could potentially promote the development of tolerance against chemical toxicants.

Prebiotic oligosaccharides are widely used as human and animal feed additives for their beneficial effects on the gut microbiota. However, there are limited data to assess the direct effect of such functional foods on the transcriptome of intestinal epithelial cells. The purpose of this study is to describe the differential transcriptomes and cellular pathways of colonic cells directly exposed to galacto-oligosaccharides (GOS) and fructo-oligosaccharides (FOS). We have examined the differential gene expression of polarized Caco-2 cells treated with GOS or FOS and their respective mock-treated cells using mRNA sequencing (RNA-seq). A total of 89 significant differentially expressed genes were identified between GOS and mock-treated groups. For FOS treatment, a reduced number of 12 significant genes were observed to be differentially expressed relative to the control group. KEGG and Gene Ontology functional analysis revealed that genes up-regulated in the presence of GOS were involved in digestion and absorption processes, fatty acids and steroids metabolism, potential antimicrobial proteins, energy-dependent and -independent transmembrane trafficking of solutes and amino acids. Using our data, we have established complementary non-prebiotic modes of action for these frequently used dietary fibers.

Combined awareness about the power and limitations of bioinformatics and molecular biology enables advanced research based on high-throughput data. Despite an increasing demand for scientists with a combined background in both fields, the education in dry lab and wet lab is often separated. This work describes an example of integrated education with focus on genomics and transcriptomics. Participants learn computational and molecular biology methods in the same practical course. Peer-review is applied as a teaching method to foster cooperative learning of students with heterogeneous backgrounds. Evaluation results indicate acceptance and appreciation of this approach.

The endometrium, the inner mucosal lining of the uterus, undergoes complex molecular and cellular changes across the menstrual cycle in preparation for embryo implantation. Transcrip-tome-wide analyses have mainly been utilized to study the endometrial receptivity, the prerequisite for successful implantation, with most of the studies, so far, comparing the endometrial transcriptomes between (i) secretory versus proliferative endometrium or (ii) mid-secretory versus early-secretory endometrium. In the current study, we provide the complete transcriptome description of the endometrium across the entire menstrual cycle and, for the first time, comprehensively characterize the proliferative phase of the endometrium. Our temporal transcriptome analysis includes five time points: mid-proliferative, late-proliferative (peri-ovulatory phase), early-secretory, mid-secretory, and late-secretory phases. Thus, we exhaustively unveil the transitions between the consecutive proliferative and secretory phases, highlighting their unique gene expression profiles and possible distinct biological functions. Transcriptome analysis reveals many differentially expressed genes (DEGs) across the menstrual cycle, most of which are phase-specific. As an example of the coordinated gene activity, the expression profile of his-tone-encoding genes within the HIST cluster on chromosome 6 showed an increase in cluster activity during the late-proliferative and a decline during the mid-secretory phase. Moreover, numerous DEGs are shared among all phases. In conclusion, in the current study, we delineate the endometrial proliferative phase-centered view of transcriptome dynamics across the menstrual cycle. Data analysis highlights significant transcriptomic and functional changes occurring during the late-proliferative phase – an essential transition point from the proliferative to the secretory phase. Future studies should explore how the biology of the late-proliferative phase endometrium impacts the achievement of mid-secretory endometrial receptivity or contributes to molecular aberrations leading to embryo implantation failure.

Coffee Coffea arabica L. depends on abundantly distributed rainfall and drought impacts their development, fruit production, bean quality, and ultimately, beverage quality. Plant Biotechnology by means of genetic manipulation and plant regeneration by the somatic embryogenic process is an alternative technology to overcome these prob-lems. In the present work we used the molecular approach of the trehalose/trehalase role in plant surviving in extreme drought/salt environments. We used a cassette containing the antisense C. arabica L. trehalase gene under the RD29 promoter from A. thaliana and the NOS terminator to genetically modified an em-bryogenic coffee C. arabica L. cv typica line under osmotic stress supplemented with mannitol (0.3M) and sorbitol (0.3M) containing medium. Drought resistant SE lines were recovered and regenerate plants. Resistant SE lines showed a higher rate of competence to induce secondary SE capacity. A transcrip-tome analysis was performed from high competent SE lines to understand the acquired molecular mechanisms to developed SE. From the up-regulated genes, a PPI network made by STRING v12.0 with high confidence (0.700) revealed the presence of the 10 modules. The cell cycle, chromatin remodeling, somatic embryogenesis, oxidative stress, generic transcription pathway, carbon metabolism, phenylpropanoid biosynthesis, trehalose biosynthesis, proline bio-synthesis and glycerolipid metabolism.

In this study, a structural profiling of equid alphaherpesvirus 1 (EHV-1) transcriptome was carried out using next-generation (Illumina) and third-generation (Oxford Nanopore Technologies) sequencing platforms. We annotated the canonical mRNA molecules and their isoforms, including transcript start and end site isoforms, and splice variants. Additionally, a number of putative 5′-truncated mRNAs containing shorter in-frame ORFs were detected. We also demonstrated that EHV-1 produces a high number of non-coding transcripts, including antisense and intergenic RNAs. One of the most remarkable features of the EHV-1 is the generation of abundant fusion transcripts some of which encoding chimeric polypeptides. We observed a higher number of splicing and transcriptional overlaps than in related viruses. Additionally, we found that many upstream genes of tandem gene clusters have their own transcript end sites (TESs) besides the co-terminal TESs, which is rare in other alphaherpesviruses. We show here that the replication origins (OriS and OriL) of the virus are co-localized with promoter sequences and overlap with specific RNA molecules. Furthermore, we discovered a novel non-coding RNA (designated as NOIR) that overlaps the 5′-ends of the longer transcript variants encoded by the two main transactivator genes ORF64 and 65 bracketing the OriL. These all suggest the existence of a central regulatory system which controls the genome-wide transcription and the replication through a mechanism based on the interference between the machineries carrying out the synthesis of DNA and RNA.

Numerous dynamic and complicated processes characterize development from the oocyte to the embryo. However, given the importance of functional transcriptome profiles, long non-coding RNAs, single-nucleotide polymorphisms, and alternative splicing during embryonic development, the effect that these features have on the blastomeres of 2-, 4-, 8-, 16-cell, and morula stages of development have not been studied. Here, we conducted a scRNA-seq survey of cells from sheep from the oocyte to the blastocyst developmental stages. We then carried out experiments to identify and functionally analyze the transcriptome profiles, long non-coding RNAs, single-nucleotide polymorphisms (SNPs), and alternative splicing (AS). We founded that between the oocyte and zygote groups significantly down-regulated genes and the second-largest change in gene expression occurred between the 8- and 16-cell stages. We used various methods to construct a profile to characterize cellular and molecular features and systematically analyze the related GO and KEGG profile of cells of all stages from the oocyte to the blastocyst. This large-scale, single-cell atlas provides key cellular information and will likely assist clinical studies in improving preimplantation genetic diagnosis.

Ophelia syndrome is characterized by the coincidence of severe neuropsychiatric symptoms, classical Hodgkin lymphoma, and the presence of antibodies to the metabotropic glutamate 5 receptor (mGluR5). Little is known about the pathogenetic link between these symptoms and the role that anti-mGluR5-antibodies play. We investigated lymphoma tissue from patients with Ophelia syndrome and with isolated classical Hodgkin lymphoma by quantitative immunocytochemistry for mGluR5-expression. Further, we studied the L-1236, L-428, L-540, SUP-HD1, KM-H2, and HDLM-2 classical Hodgkin lymphoma cell lines by FACS and Western blot for mGluR5-expression, and by transcriptome analysis. mGluR5 surface expression differed significantly in terms of receptor density, distribution pattern, and percentage of positive cells. Highest expression levels were found in the L-1236 line. RNA-sequencing revealed more than 800 genes that were higher expressed in the L-1236 line in comparison to the other classical Hodgkin lymphoma cell lines. High mGluR5-expression was associated with upregulation of PI3K/AKT and MAPK pathways and of downstream targets (e.g. EGR1) known to be involved in classical Hodgkin lymphoma progression. Finally, mGluR5 expression was increased in the classical Hodgkin lymphoma-tissue of our Ophelia syndrome patient in contrast to five classical Hodgkin lymphoma-patients without autoimmune encephalitis. Given the association of encephalitis and classical Hodgkin lymphoma in Ophelia syndrome, it is possible that mGluR5-expression on classical Hodgkin lymphoma cells not only drives tumor progression but also triggers anti-mGluR5 encephalitis even before classical Hodgkin lymphoma becomes manifest.

Leishmania major is the main causative agent of cutaneous leishmaniasis in humans. The Friedlin strain of this species (LmjF) was chosen when a multi-laboratory consortium undertook the objective of deciphering the first genome sequence for a parasite of the genus Leishmania. The objective was successfully attained in 2005, and this represented a milestone for Leishmania molecular biology studies around the world. Although the LmjF genome sequence was done following a shotgun strategy and using classical Sanger sequencing, the results were excellent and this genome assembly served as the reference for subsequent genome assemblies in other Leishmania species. Here, we present a new assembly for the genome of this strain (named LMJFC for clarity), generated by the combination of two high throughput sequencing platforms, Illumina short-read sequencing and PacBio Single Molecular Real-Time (SMRT) sequencing, which provides long-read sequences. Apart from resolving uncertain nucleotide positions, several genomic regions have been reorganized and a more precise composition of tandemly repeated gene loci was attained. Additionally, the genome annotation has been improved by adding 542 genes and more accurate coding-sequences defined for around two hundred genes, based on the transcriptome delimitation also carried out in this work. As a result, we are providing gene models (including untranslated regions and introns) for 11,238 genes. Genomic information ultimately determines the biology of every organism; therefore, our understanding of molecular mechanisms will depend on the availability of precise genome sequences and accurate gene annotations. In this regards, this work is providing an improved genome sequence and updated transcriptome annotations for the reference L. major Friedlin strain.

Climate change-induced temperature fluctuations pose a significant threat to crop production, particularly in the Southern Hemisphere. This study investigates the transcriptome and physiological responses of rapeseed to post-flowering temperature increases, providing valuable insights into the molecular mechanisms underlying rapeseed tolerance to heat stress. Two rapeseed genotypes, Lumen and Solar, were assessed under control and heat stress conditions in field experiments conducted in Valdivia, Chile. Results showed that seed yield and seed number were negatively affected by heat stress, with genotype-specific responses. Lumen exhibited a 9.3% average seed yield reduction, while Solar showed a 28.7% reduction. RNA-seq analysis of siliques and seeds revealed tissue-specific responses to heat stress, with siliques being more sensitive to temperature stress. Hierarchical clustering analysis identified distinct gene clusters reflecting different aspects of heat stress adaptation in siliques, with a role for protein folding in maintaining silique development and seed quality under high temperature conditions. In seeds, three distinct patterns of heat-responsive gene expression were observed, with genes involved in protein folding and response to heat showing genotype-specific expression. Gene coexpression network analysis revealed major modules for rapeseed yield and quality, as well as the trade-off between seed number and seed weight. Overall, this study contributes to understanding the molecular mechanisms underlying rapeseed tolerance to heat stress and can inform crop improvement strategies targeting yield optimization under changing environmental conditions.

Truffles are ascomycete hypogeous fungi belonging to the Tuberaceae family of the Pezizales order that grow in ectomycorrhizal symbiosis with tree roots and are known for their peculiar aromas and flavors. Axenic culture of truffle mycelium is problematic because it is not possible in many cases, and the growth rate is meager when it is possible. This limitation prompts searching and characterizing new strains that can be handled in laboratory conditions for basic and applied studies. In this work, a new strain of Tuber borchii (strain SP1) has been isolated and cultured, and its transcriptome has been analyzed under different in vitro culture conditions. The results show that the best T. borchii SP1 growth was obtained using maltose-enriched cultures made with soft-agar and in static submerged cultures made at 22ºC. The transcriptome analysis of this strain cultured in different media indicated that most of the gene transcription effort is due to a limited number of genes (20% of genes account for 80% of the transcription), that the transcription profile of the central metabolism genes was similar in the different conditions analyzed with a transcription signal detected for around 80% of the annotated genes. The gene expression profile suggests that T. borchii uses a fermentative rather than respiratory metabolism, even in aerobic conditions. Finally, there is a reduced expression of genes belonging to secondary metabolite clusters, whereas there is a significative transcription of those involved in producing volatile aromatic compounds.

The transcriptional response of human blood leukocytes to SARS-CoV-2 infection was investigated focusing on the differences between mild and severe cases and between age subgroups. Weighted gene co-expression network analysis and comparative gene expression analysis were used. Three transcriptional modules positively associated with the traits of interest and their respective high hierarchy genes were identified. Enrichment analyses showed that the yellow module, associated with severe cases and older patients, had an overrepresentation of genes involved in inflammatory and innate immune responses, and neutrophil activation. The magenta and black modules, associated with disease severity and younger patients, contained genes related to innate immunity and inflammation and genes that regulate those responses. Subnetworks for these modules were constructed using genes enriched for innate immunity, inflammation, immunoregulation and differentially expressed genes (severe vs. mild). Their analysis evidenced that immunoregulatory functions are more activated in the modules associated with younger patients, what may help to explain the better disease course and faster recovery observed in younger COVID-19 patients. Comparative gene expression analysis between severe and mild groups, followed by gene enrichment and normalized gene expression analyses, revelated a set of 23 potential biomarkers for COVID-19 severity, of which 13 are newly described.

Alternative splicing (AS) is a gene regulatory mechanism modulating gene expression in multiple ways. AS is prevalent in all eukaryotes including plants. AS generates two or more mRNAs from the precursor mRNA (pre-mRNA) to regulate transcriptome complexity and proteome diversity. Advances in next-generation sequencing, omics technology and bioinformatics tools, and computational methods provide new opportunities to quantify and visualize AS-based quantitative trait variation associated with plant growth, development, reproduction, and stress tolerance. Domestication, polyploidization and environmental perturbation may evolve novel splicing variants associated with agronomically beneficial traits. To date, pre-mRNAs from many genes are spliced into multiple transcripts that cause phenotypic variation for complex traits, both in model plant Arabidopsis and field crops. Cataloguing and exploiting such variation may provide new paths to enhance climate resilience, resource-use efficiency, productivity, and nutritional quality of staple food crops. This review provides insights into AS variation alongside gene expression analysis to select for novel phenotypic diversity for use in breeding programs. AS contributes to heterosis, enhances plant symbiosis (mycorrhiza and rhizobium), and provides a mechanistic link between the core clock genes and diverse environmental clues.

The rapid expansion of transcriptomics from increased affordability of next-generation sequencing (NGS) technologies generates rocketing amounts of gene expression data across biology and medicine, and notably in cancer research. Concomitantly, many bioinformatics tools were developed to streamline gene expression analysis and quantification. We tested the concordance of NGS RNA sequencing (RNA-seq) analysis outcomes between the two predominant programs for reads alignment, HISAT2 and STAR, and the two most popular programs for quantifying gene expression in NGS experiments, edgeR and DESeq2, using RNA-seq data from a series of breast cancer progression specimens, which include histologically confirmed normal, early neoplasia, ductal carcinoma in situ and infiltrating ductal carcinoma samples microdissected from formalin fixed, paraffin embedded (FFPE) breast tissue blocks. We identified significant differences in aligners’ performance: HISAT2 was prone to misalign reads to retrogene genomic loci, STAR generated more precise alignments, especially for early neoplasia samples. edgeR and DESeq2 produced similar lists of differentially expressed genes in stage comparisons, with edgeR producing more conservative, though shorter, lists of genes. Albeit, Gene Ontology (GO) enrichment analysis revealed no skewness in significant GO categories identified among differentially expressed genes by edgeR vs DESeq2. As transcriptome analysis of archived FFPE samples becomes a vanguard of precision medicine, identification and fine-tuning of bioinformatics tools becomes critical for clinical research. Our results indicate that STAR and edgeR are well-suited tools for differential gene expression analysis from FFPE samples.

Coccomyxa subellipsoidea KJ (C-KJ) is a green alga with unique immunoregulatory characteristics. Here, we investigated the mechanism underlying the modification of T cell function by C-KJ components. The water-soluble extract of C-KJ was fractionated into protein (P) and sugar (S) fractions [acidic (AS), basic (BS), and neutral (NS)]. These fractions were used for the treatment of peripheral blood mononuclear cells stimulated with toxic shock syndrome toxin-1. Transcriptome analysis revealed that both P and AS enhanced the expression of the genes encoding metallothionein (MT) family proteins, inflammatory factors, and T helper (Th) 17 cytokine and suppressed that of those encoding Th2 cytokines in stimulated T cells. The kinetics of MT1 and MT2A gene expression showed a transient increase in MT1 and maintenance of MT2A mRNA after T cell stimulation in the presence of AS. The kinetics of Th17-related cytokine secretion in the early period were comparable to those of MT2A mRNA. Furthermore, our findings revealed that static, a STAT-3 inhibitor, significantly suppressed MT2A gene expression. These findings suggest that the expression of MTs is involved in the immune regulatory function of C-KJ components, which is partially regulated by Th17 responses, and may help develop innovative immunoregulatory drugs or functional foods.

Stem/progenitor cell transplantation is a potential novel therapeutic strategy to induce angiogenesis in ischemic tissue, which can prevent major amputation in patients with advanced peripheral artery disease (PAD). Thus, clinicians can use cell therapies worldwide to treat PAD. However, some cell therapy studies did not report beneficial outcomes. Clinical researchers suggested that classical risk factors and comorbidities may adversely affect the efficacy of cell therapy. Some studies have indicated that the response to stem cell therapy varies among patients even in those harboring limited risk factors. This suggested the role of undetermined risk factors, including genetic alterations, somatic mutations, and clonal hematopoiesis. Personalized stem cell-based therapy can be developed by analyzing individual risk factors. These approaches must consider several clinical biomarkers and perform studies (such as genome-wide association studies (GWAS)) on disease-related genetic traits and integrate the findings with those of transcriptome-wide association studies (TWAS) and whole-genome sequencing in PAD. Additional unbiased analyses with state-of-the-art computational methods, such as machine learning-based patient stratification, are suited for predictions in clinical investigations. The integration of these complex approaches into a unified analysis procedure for the identification of responders and non-responders before stem cell therapy, which can decrease treatment expenditure, is a major challenge to increase the efficacy of therapies.