Characterization of Sf9 cell clones with differential susceptibilities to Sf-rhabdovirus X+3.7 and Sf-rhabdovirus X- replication

Our laboratory previously discovered a novel rhabdovirus in the Spodoptera frugiperda Sf9 insect cell line that was designated as Sf-rhabdovirus. Using limiting dilution, this cell line was found to be a mixed population of cells infected by Sf-rhabdovirus variants containing either the full length X accessory gene with a 3.7 kb internal duplication (designated as Sf-rhabdovirus X+3.7) or lacking the duplication and part of the X gene (designated as Sf-rhabdovirus X-), and cells that were negative for Sf-rhabdovirus. In this paper, we found that the Sf-rhabdovirus negative cell clones had sub-populations with different susceptibilities to the replication of Sf-rhabdovirus X+3.7 and X- variants: cell clone Sf9-13F12 was more sensitive to replication by both virus variants compared to Sf9-3003; moreover, Sf9-3003 showed more resistance to X+3.7 replication than to X- replication. RNA-Seq analysis indicated significant differentially expressed genes in the Sf9-13F12 and Sf9-3003 cell clones further supporting that distinct sub-populations of virus-negative cells co-exist in the parent Sf9 cell line.

A rapid return to normal: temporal gene expression patterns following cold exposure in the bumble bee Bombus impatiens

Bumble bees are common in cooler climates and many species likely experience periodic exposure to very cold temperatures, but little is known about the temporal dynamics of cold response mechanisms following chill exposure, especially how persistent effects of cold exposure may facilitate tolerance of future events. To investigate molecular processes involved in the temporal response by bumble bees to acute cold exposure, we compared mRNA transcript abundance in Bombus impatiens workers exposed to 0°C for 75 min (inducing chill coma) and control bees maintained at a constant ambient temperature (28°C). We sequenced the 3' end of mRNA transcripts (TagSeq) to quantify gene expression in thoracic tissue of bees at several time points (0, 10, 30, 120 and 720 min) following cold exposure. Significant differences from control bees were only detectable within 30 min after the treatment, with most occurring at the 10 min recovery time point. Genes associated with gluconeogenesis and glycolysis were most notably upregulated, while genes related to lipid and purine metabolism were downregulated. The observed patterns of expression indicate a rapid recovery after chill coma, suggesting an acute differential transcriptional response during recovery from chill coma and return to baseline expression levels within an hour, with no long-term gene expression markers of this cold exposure. Our work highlights the functions and pathways important for acute cold recovery, provides an estimated time frame for recovery from cold exposure in bumble bees, and suggests that cold hardening may be less important for these heterothermic insects.

Integrated bioinformatics analysis and experimental animal models identify a robust biomarker and its correlation with the immune microenvironment in pulmonary arterial hypertension

Background

Pulmonary arterial hypertension (PAH) represents a substantial global risk to human health. This study aims to identify diagnostic biomarkers for PAH and assess their association with the immune microenvironment through the utilization of sophisticated bioinformatics techniques.

Methods

Based on two microarray datasets, differentially expressed genes (DEGs) were detected, and hub genes underwent a sequence of machine learning analyses. After pathways associated with PAH were assessed by gene enrichment analysis, the identified genes were validated using external datasets and confirmed in a monocrotaline (MCT)-induced rat model. In addition, three algorithms were employed to estimate the proportions of various immune cell types, and the link between hub genes and immune cells was substantiated.

Results

Using SVM, LASSO, and WGCNA, we identified seven hub genes, including (BPIFA1, HBA2, HBB, LOC441081, PI15, S100A9, and WIF1), of which only BPIFA1 remained stable in the external datasets and was validated in an MCT-induced rat model. Furthermore, the results of the functional enrichment analysis established a link between PAH and both metabolism and the immune system. Correlation assessment showed that BPIFA1 expression in the MCP-counter algorithm was negatively associated with various immune cell types, positively correlated with macrophages in the ssGSEA algorithm, and correlated with M1 and M2 macrophages in the CIBERSORT algorithm.

Conclusion

BPIFA1 serves as a modulator of PAH, with the potential to impact the immune microenvironment and disease progression, possibly through its regulatory influence on both M1 and M2 macrophages.

Differential RNA Expression Between Metastatic and Primary Neuroblastoma Cells

INTRODUCTION

Neuroblastoma (NB) is the most common extra-cranial malignancy in children. Poor survival in high-risk NB is attributed to recurrent metastatic disease. To better study metastatic disease, we used a novel mouse model to investigate differential gene expression between primary tumor cells and metastatic cells. We hypothesized that metastatic NB cells have a different gene expression profile from primary tumor cells and cultured cells.

METHODS

Using three human NB cell lines (NGP, CHLA255, and SH-SY5Y), orthotopic xenografts were established in immunodeficient nod/scid gamma mice via subcapsular renal injection. Mice were sacrificed and NB cells were isolated from the primary tumor and from sites of metastasis (bone marrow, liver). RNA sequencing, gene set analysis, and pathway analysis were performed to identify differentially expressed genes and molecular pathways in the metastatic cells compared to primary tumor cells.

RESULTS

There were 266 differentially expressed genes in metastatic tumor cells (bone marrow and liver combined) compared to primary tumor cells. The top upregulated gene was KCNK1 and the top downregulated genes were PDE7B and NEBL. Top upregulated pathways in the metastatic cells were involved in ion transport, cell signaling, and cell proliferation. Top downregulated pathways were involved in DNA synthesis, transcription, and cellular metabolism.

CONCLUSIONS

In metastatic NB cells, our study identified the upregulation of biologic processes involved in cell cycle regulation, cell proliferation, migration, and invasion. Ongoing studies aim to validate downstream translation of these genomic alterations, as well as target these pathways to more effectively suppress and inhibit recurrent metastatic disease in NB.

Transcriptomic analysis of Anopheles gambiae from Benin reveals overexpression of salivary and cuticular proteins associated with cross-resistance to pyrethroids and organophosphates

BACKGROUND

Insecticide resistance (IR) is one of the major threats to malaria vector control programs in endemic countries. However, the mechanisms underlying IR are poorly understood. Thus, investigating gene expression patterns related to IR can offer important insights into the molecular basis of IR in mosquitoes. In this study, RNA-Seq was used to characterize gene expression in Anopheles gambiae surviving exposure to pyrethroids (deltamethrin, alphacypermethrin) and an organophosphate (pirimiphos-methyl).

RESULTS

Larvae of An. gambiae s.s. collected from Bassila and Djougou in Benin were reared to adulthood and phenotyped for IR using a modified CDC intensity bottle bioassay. The results showed that mosquitoes from Djougou were more resistant to pyrethroids (5X deltamethrin: 51.7% mortality; 2X alphacypermethrin: 47.4%) than Bassila (1X deltamethrin: 70.7%; 1X alphacypermethrin: 77.7%), while the latter were more resistant to pirimiphos-methyl (1.5X: 48.3% in Bassila and 1X: 21.5% in Djougou). RNA-seq was then conducted on resistant mosquitoes, non-exposed mosquitoes from the same locations and the laboratory-susceptible An. gambiae s.s. Kisumu strain. The results showed overexpression of detoxification genes, including cytochrome P450s (CYP12F2, CYP12F3, CYP4H15, CYP4H17, CYP6Z3, CYP9K1, CYP4G16, and CYP4D17), carboxylesterase genes (COEJHE5E, COE22933) and glutathione S-transferases (GSTE2 and GSTMS3) in all three resistant mosquito groups analyzed. Genes encoding cuticular proteins (CPR130, CPR10, CPR15, CPR16, CPR127, CPAP3-C, CPAP3-B, and CPR76) were also overexpressed in all the resistant groups, indicating their potential role in cross resistance in An. gambiae. Salivary gland protein genes related to 'salivary cysteine-rich peptide' and 'salivary secreted mucin 3' were also over-expressed and shared across all resistant groups.

CONCLUSION

Our results suggest that in addition to metabolic enzymes, cuticular and salivary gland proteins could play an important role in cross-resistance to multiple classes of insecticides in Benin. These genes warrant further investigation to validate their functional role in An. gambiae resistance to insecticides.

Long noncoding RNA expression in acute lymphoblastic leukemia: A systematic review

Long noncoding RNAs (lncRNAs), as gene expression modulators, are potential players in Acute Lymphoblastic Leukemia (ALL) pathogenesis. We systematically explored current literature on lncRNA expression in ALL to identify lncRNAs consistently reported as differentially expressed (DE) either in ALL versus controls or between ALL subtypes. By comparing articles that provided global expression data for DE lncRNAs in the ETV6::RUNX1-positive ALL subtype, we identified four DE lncRNAs in three independent studies (two versus other subtypes and one versus controls), showing concordant expression of LINC01013, CRNDE and lnc-KLF7-1. Additionally, LINC01503 was consistently downregulated on ALL versus controls. Within RT-qPCR studies, twelve lncRNA were DE in more than one source. Thus, several lncRNAs were supported as DE in ALL by multiple sources, highlighting their potential role as candidate biomarkers or therapeutic targets. Finally, as lncRNA annotation is rapidly expanding, standardization of reporting and nomenclature is urgently needed to improve data verifiability and compilation.

RNA-seq data exploration after trypanosome RNA-binding protein UBP1 expression is altered by CRISPR-Cas9 gene editing and overexpression

Previous studies have shown that overexpression of the Trypanosoma cruzi U-rich RNA-binding protein 1 (TcUBP1) in insect-dwelling epimastigotes results in a gene expression pattern resembling that of the infective form of the pathogen. Here, we used CRISPR-Cas9-induced edition of TcUBP1 and full-length protein overexpression in epimastigote cells to monitor transcriptomic changes during the epimastigote-to-metacyclic trypomastigote stage transition of T. cruzi. This dataset includes the bioinformatics analysis of three different RNA-seq samples, each with three biological replicates, showing differential mRNA abundances. The current transcriptome report has the potential to shed light on the quantitative variances in the expression of significant up- or down-regulated mRNAs as a consequence of the levels of the UBP1 protein. Raw data files were deposited at the NCBI Sequence Read Archive - SRA at http://ncbi.nlm.nih.gov/Traces/sra/sra.cgi with accession numbers PRJNA907231 and PRJNA949967.

Dissecting gene expression networks in the developing hippocampus through the lens of NEIL3 depletion

Gene regulation in the hippocampus is fundamental for its development, synaptic plasticity, memory formation, and adaptability. Comparisons of gene expression among different developmental stages, distinct cell types, and specific experimental conditions have identified differentially expressed genes contributing to the organization and functionality of hippocampal circuits. The NEIL3 DNA glycosylase, one of the DNA repair enzymes, plays an important role in hippocampal maturation and neuron functionality by shaping transcription. While differential gene expression (DGE) analysis has identified key genes involved, broader gene expression patterns crucial for high-order hippocampal functions remain uncharted. By utilizing the weighted gene co-expression network analysis (WGCNA), we mapped gene expression networks in immature (p8-neonatal) and mature (3 m-adult) hippocampal circuits in wild-type and NEIL3-deficient mice. Our study unveiled intricate gene network structures underlying hippocampal maturation, delineated modules of co-expressed genes, and pinpointed highly interconnected hub genes specific to the maturity of hippocampal subregions. We investigated variations within distinct gene network modules following NEIL3 depletion, uncovering NEIL3-targeted hub genes that influence module connectivity and specificity. By integrating WGCNA with DGE, we delve deeper into the NEIL3-dependent molecular intricacies of hippocampal maturation. This study provides a comprehensive systems-level analysis for assessing the potential correlation between gene connectivity and functional connectivity within the hippocampal network, thus shaping hippocampal function throughout development.

Curare and GenExVis: a versatile toolkit for analyzing and visualizing RNA-Seq data

Even though high-throughput transcriptome sequencing is routinely performed in many laboratories, computational analysis of such data remains a cumbersome process often executed manually, hence error-prone and lacking reproducibility. For corresponding data processing, we introduce Curare, an easy-to-use yet versatile workflow builder for analyzing high-throughput RNA-Seq data focusing on differential gene expression experiments. Data analysis with Curare is customizable and subdivided into preprocessing, quality control, mapping, and downstream analysis stages, providing multiple options for each step while ensuring the reproducibility of the workflow. For a fast and straightforward exploration and visualization of differential gene expression results, we provide the gene expression visualizer software GenExVis. GenExVis can create various charts and tables from simple gene expression tables and DESeq2 results without the requirement to upload data or install software packages. In combination, Curare and GenExVis provide a comprehensive software environment that supports the entire data analysis process, from the initial handling of raw RNA-Seq data to the final DGE analyses and result visualizations, thereby significantly easing data processing and subsequent interpretation.

Transcriptome profiling and metabolic pathway analysis towards reliable biomarker discovery in early-stage lung cancer

Earlier diagnosis of lung cancer is crucial for reducing mortality and morbidity in high-risk patients. Liquid biopsy is a critical technique for detecting the cancer earlier and tracking the treatment outcomes. However, noninvasive biomarkers are desperately needed due to the lack of therapeutic sensitivity and early-stage diagnosis. Therefore, we have utilized transcriptomic profiling of early-stage lung cancer patients to discover promising biomarkers and their associated metabolic functions. Initially, PCA highlights the diversity level of gene expression in three stages of lung cancer samples. We have identified two major clusters consisting of highly variant genes among the three stages. Further, a total of 7742, 6611, and 643 genes were identified as DGE for stages I-III respectively. Topological analysis of the protein-protein interaction network resulted in seven candidate biomarkers such as JUN, LYN, PTK2, UBC, HSP90AA1, TP53, and UBB cumulatively for the three stages of lung cancers. Gene enrichment and KEGG pathway analyses aid in the comprehension of pathway mechanisms and regulation of identified hub genes in lung cancer. Importantly, the medial survival rates up to ~ 70 months were identified for hub genes during the Kaplan-Meier survival analysis. Moreover, the hub genes displayed the significance of risk factors during gene expression analysis using TIMER2.0 analysis. Therefore, we have reason that these biomarkers may serve as a prospective targeting candidate with higher treatment efficacy in early-stage lung cancer patients.

Discovering the regulators of heat stress tolerance in Ziziphus nummularia (Burm.f) wight and walk.-arn

Ziziphus nummularia an elite heat-stress tolerant shrub, grows in arid regions of desert. However, its molecular mechanism responsible for heat stress tolerance is unexplored. Therefore, we analysed whole transcriptome of Jaisalmer (heat tolerant) and Godhra (heat sensitive) genotypes of Z. nummularia to understand its molecular mechanism responsible for heat stress tolerance. De novo assembly of 16,22,25,052 clean reads yielded 276,029 transcripts. A total of 208,506 unigenes were identified which contains 4290 and 1043 differentially expressed genes (DEG) in TGO (treated Godhra at 42 °C) vs. CGO (control Godhra) and TJR (treated Jaisalmer at 42 °C) vs. CJR (control Jaisalmer), respectively. A total of 987 (67 highly enriched) and 754 (34 highly enriched) pathways were obsorved in CGO vs. TGO and CJR vs. TJR, respectively. Antioxidant pathways and TFs like Homeobox, HBP, ARR, PHD, GRAS, CPP, and E2FA were uniquely observed in Godhra genotype and SET domains were uniquely observed in Jaisalmer genotype. Further transposable elements were highly up-regulated in Godhra genotype but no activation in Jaisalmer genotype. A total of 43,093 and 39,278 simple sequence repeats were identified in the Godhra and Jaisalmer genotypes, respectively. A total of 10 DEGs linked to heat stress were validated in both genotypes for their expression under different heat stresses using quantitative real-time PCR. Comparing expression patterns of the selected DEGs identified ClpB1 as a potential candidate gene for heat tolerance in Z. nummularia. Here we present first characterized transcriptome of Z. nummularia in response to heat stress for the identification and characterization of heat stress-responsive genes.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-024-01431-y.

Transcriptomic analysis identifies candidate genes for Aphanomyces root rot disease resistance in pea

BACKGROUND

Aphanomyces euteiches is a soil-borne oomycete that causes root rot in pea and other legume species. Symptoms of Aphanomyces root rot (ARR) include root discoloration and wilting, leading to significant yield losses in pea production. Resistance to ARR is known to be polygenic but the roles of single genes in the pea immune response are still poorly understood. This study uses transcriptomics to elucidate the immune response of two pea genotypes varying in their levels of resistance to A. euteiches.

RESULTS

In this study, we inoculated roots of the pea (P. sativum L.) genotypes 'Linnea' (susceptible) and 'PI180693' (resistant) with two different A. euteiches strains varying in levels of virulence. The roots were harvested at 6 h post-inoculation (hpi), 20 hpi and 48 hpi, followed by differential gene expression analysis. Our results showed a time- and genotype-dependent immune response towards A. euteiches infection, involving several WRKY and MYB-like transcription factors, along with genes associated with jasmonic acid (JA) and abscisic acid (ABA) signaling. By cross-referencing with genes segregating with partial resistance to ARR, we identified 39 candidate disease resistance genes at the later stage of infection. Among the genes solely upregulated in the resistant genotype 'PI180693', Psat7g091800.1 was polymorphic between the pea genotypes and encoded a Leucine-rich repeat receptor-like kinase reminiscent of the Arabidopsis thaliana FLAGELLIN-SENSITIVE 2 receptor.

CONCLUSIONS

This study provides new insights into the gene expression dynamics controlling the immune response of resistant and susceptible pea genotypes to A. euteiches infection. We present a set of 39 candidate disease resistance genes for ARR in pea, including the putative immune receptor Psat7g091800.1, for future functional validation.

A regulatory element associated to NAFLD in the promoter of DIO1 controls LDL-C, HDL-C and triglycerides in hepatic cells

BACKGROUND

Genome-wide association studies (GWAS) have identified genetic variants linked to fat metabolism and related traits, but rarely pinpoint causative variants. This limitation arises from GWAS not considering functional implications of noncoding variants that can affect transcription factor binding and potentially regulate gene expression. The aim of this study is to investigate a candidate noncoding functional variant within a genetic locus flagged by a GWAS SNP associated with non-alcoholic fatty liver disease (NAFLD), a condition characterized by liver fat accumulation in non-alcohol consumers.

METHODS

CRISPR-Cas9 gene editing in HepG2 cells was used to modify the regulatory element containing the candidate functional variant linked to NAFLD. Global gene expression in mutant cells was assessed through RT-qPCR and targeted transcriptomics. A phenotypic assay measured lipid droplet accumulation in the CRISPR-Cas9 mutants.

RESULTS

The candidate functional variant, rs2294510, closely linked to the NAFLD-associated GWAS SNP rs11206226, resided in a regulatory element within the DIO1 gene's promoter region. Altering this element resulted in changes in transcription factor binding sites and differential expression of candidate target genes like DIO1, TMEM59, DHCR24, and LDLRAD1, potentially influencing the NAFLD phenotype. Mutant HepG2 cells exhibited increased lipid accumulation, a hallmark of NAFLD, along with reduced LDL-C, HDL-C and elevated triglycerides.

CONCLUSIONS

This comprehensive approach, that combines genome editing, transcriptomics, and phenotypic assays identified the DIO1 promoter region as a potential enhancer. Its activity could regulate multiple genes involved in the NAFLD phenotype or contribute to defining a polygenic risk score for enhanced risk assessment in NAFLD patients.

Differential Expression of Nuclear-Encoded Mitochondrial Protein Genes of ATP Synthase Across Different Tissues of Female Buffalo

The physiological well-being of buffaloes, encompassing phenotypic traits, reproductive health, and productivity, depends on their energy status. Mitochondria, the architects of energy production, orchestrate a nuanced interplay between nuclear and mitochondrial domains. Oxidative phosphorylation complexes and associated proteins wield significant influence over metabolic functions, energy synthesis, and organelle dynamics, often linked to tissue-specific pathologies. The unexplored role of ATP synthase in buffalo tissues prompted a hypothesis: in-depth exploration of nuclear-derived mitochondrial genes, notably ATP synthase, reveals distinctive tissue-specific diversity. RNA extraction and sequencing of buffalo tissues (kidney, heart, brain, and ovary) enabled precise quantification of nuclear-derived mitochondrial protein gene expression. The analysis unveiled 24 ATP synthase transcript variants, each with unique tissue-specific patterns. Kidney, brain, and heart exhibited elevated gene expression compared to ovaries, with 10, 8, and 19 up-regulated genes, respectively. The kidney showed 3 and 12 down-regulated genes compared to the brain and heart. The heart-brain comparison highlighted ten highly expressed genes in ATP synthase functions. Gene ontology and pathway analyses revealed enriched functions linked to ATP synthesis and oxidative phosphorylation, offering a comprehensive understanding of energy production in buffalo tissues. This analysis enhances understanding of tissue-specific gene expression, emphasizing the influence of energy demands. Revealing intricate links between mitochondrial gene expression and tissue specialization in buffaloes, it provides nuanced insights into tissue-specific expression of nuclear-encoded mitochondrial protein genes, notably ATP synthase, advancing the comprehension of buffalo tissue biology.

Peripheral inflammation-induced changes in songbird brain gene expression: 3' mRNA transcriptomic approach

Species-specific neural inflammation can be induced by profound immune signalling from periphery to brain. Recent advances in transcriptomics offer cost-effective approaches to study this regulation. In a population of captive zebra finch (Taeniopygia guttata), we compare the differential gene expression patterns in lipopolysaccharide (LPS)-triggered peripheral inflammation revealed by RNA-seq and QuantSeq. The RNA-seq approach identified more differentially expressed genes but failed to detect any inflammatory markers. In contrast, QuantSeq results identified specific expression changes in the genes regulating inflammation. Next, we adopted QuantSeq to relate peripheral and brain transcriptomes. We identified subtle changes in the brain gene expression during the peripheral inflammation (e.g. up-regulation in AVD-like and ACOD1 expression) and detected co-structure between the peripheral and brain inflammation. Our results suggest benefits of the 3'end transcriptomics for association studies between peripheral and neural inflammation in genetically heterogeneous models and identify potential targets for the future brain research in birds.

CrustyBase v.2.0: new features and enhanced utilities to support open science

BACKGROUND

Transcriptomes present a rich, multi-dimensional subset of genomics data. They provide broad insights into genetic sequence, and more significantly gene expression, across biological samples. This technology is frequently employed for describing the genetic response to experimental conditions and has created vast libraries of datasets which shed light on gene function across different tissues, diseases, diets and developmental stages in many species. However, public accessibility of these data is impeded by a lack of suitable software interfaces and databases with which to locate and analyse them. BODY: Here we present an update on the status of CrustyBase.org, an online resource for analysing and sharing crustacean transcriptome datasets. Since its release in October 2020, the resource has provided many thousands of transcriptome sequences and expression profiles to its users and received 19 new dataset imports from researchers across the globe. In this article we discuss user analytics which point towards the utilization of this resource. The architecture of the application has proven robust with over 99.5% uptime and effective reporting of bugs through both user engagement and the error logging mechanism. We also introduce several new features that have been developed as part of a new release of CrustyBase.org. Two significant features are described in detail, which allow users to navigate through transcripts directly by submission of transcript identifiers, and then more broadly by searching for encoded protein domains by keyword. The latter is a novel and experimental feature, and grants users the ability to curate gene families from any dataset hosted on CrustyBase in a matter of minutes. We present case studies to demonstrate the utility of these features.

CONCLUSION

Community engagement with this resource has been very positive, and we hope that improvements to the service will further enable the research of users of the platform. Web-based platforms such as CrustyBase have many potential applications across life science domains, including the health sector, which are yet to be realised. This leads to a wider discussion around the role of web-based resources in facilitating an open and collaborative research community.

A chromosome-level genome for the nudibranch gastropod Berghia stephanieae helps parse clade-specific gene expression in novel and conserved phenotypes

BACKGROUND

How novel phenotypes originate from conserved genes, processes, and tissues remains a major question in biology. Research that sets out to answer this question often focuses on the conserved genes and processes involved, an approach that explicitly excludes the impact of genetic elements that may be classified as clade-specific, even though many of these genes are known to be important for many novel, or clade-restricted, phenotypes. This is especially true for understudied phyla such as mollusks, where limited genomic and functional biology resources for members of this phylum have long hindered assessments of genetic homology and function. To address this gap, we constructed a chromosome-level genome for the gastropod Berghia stephanieae (Valdés, 2005) to investigate the expression of clade-specific genes across both novel and conserved tissue types in this species.

RESULTS

The final assembled and filtered Berghia genome is comparable to other high-quality mollusk genomes in terms of size (1.05 Gb) and number of predicted genes (24,960 genes) and is highly contiguous. The proportion of upregulated, clade-specific genes varied across tissues, but with no clear trend between the proportion of clade-specific genes and the novelty of the tissue. However, more complex tissue like the brain had the highest total number of upregulated, clade-specific genes, though the ratio of upregulated clade-specific genes to the total number of upregulated genes was low.

CONCLUSIONS

Our results, when combined with previous research on the impact of novel genes on phenotypic evolution, highlight the fact that the complexity of the novel tissue or behavior, the type of novelty, and the developmental timing of evolutionary modifications will all influence how novel and conserved genes interact to generate diversity.

Immunomodulatory differences between mesenchymal stem cells from different oral tissues

Mesenchymal stem cells (MSCs) have recently been identified as having potentially therapeutic immunomodulatory properties. MSCs isolated from different oral tissues have similar morphology and immunophenotypes, however, direct comparisons of their gene expression and immunomodulatory properties have not been conducted. We isolated alveolar bone-derived MSCs (aBMSCs), dental pulp stem cells (DPSCs) and gingiva-derived MSCs (GMSCs) from the same patients and compared their immunophenotypes and transcriptomes. Additionally, we compared their production of soluble immunomodulatory cytokines as well as their immunoregulatory properties in coculture with THP-1 human monocytic cells. RNA sequencing revealed distinct gene expression in DPSCs while aBMSCs and GMSCs had less differentially expressed genes. DPSCs also had significantly less secretion of osteopontin compared to aBMSCs and GMSCs. Finally, DPSCs did not exhibit an immunosuppresive effect on THP-1 cells to the same degree as aBMSCs and GMSCs. These findings demonstrate that MSCs from different oral tissues have distinct transcriptomes and immunoregulatory properties.

Epigenetic variation mediated by lncRNAs accounts for adaptive genomic differentiation of the endemic blue mussel Mytiluschilensis

Epigenetic variation affects gene expression without altering the underlying DNA sequence of genes controlling ecologically relevant phenotypes through different mechanisms, one of which is long non-coding RNAs (lncRNAs). This study identified and evaluated the gene expression of lncRNAs in the gill and mantle tissues of Mytilus chilensis individuals from two ecologically different sites: Cochamó (41°S) and Yaldad (43°S), southern Chile, both impacted by climatic-related conditions and by mussel farming given their use as seedbeds. Sequences identified as lncRNAs exhibited tissue-specific differences, mapping to 3.54 % of the gill transcriptome and 1.96 % of the mantle transcriptome, representing an average of 2.76 % of the whole transcriptome. Using a high fold change value (≥|100|), we identified 43 and 47 differentially expressed lncRNAs (DE-lncRNAs) in the gill and mantle tissue of individuals sampled from Cochamó and 21 and 17 in the gill and mantle tissue of individuals sampled from Yaldad. Location-specific DE-lncRNAs were also detected in Cochamó (65) and Yaldad (94) samples. Via analysis of the differential expression of neighboring protein-coding genes, we identified enriched GO terms related to metabolic, genetic, and environmental information processing and immune system functions, reflecting how the impact of local ecological conditions may influence the M. chilensis (epi)genome expression. These DE-lncRNAs represent complementary biomarkers to DNA sequence variation for maintaining adaptive differences and phenotypic plasticity to cope with natural and human-driven perturbations.

Deciphering the correlation between metabolic activity through 18F-FDG-PET/CT and immune landscape in soft-tissue sarcomas: an insight from the NEOSARCOMICS study

Metabolic elevation in soft-tissue sarcomas (STS), as documented with 18F-Fluorodeoxyglucose positron emission tomography (18F-FDG-PET/CT) has been linked with cell proliferation, higher grade, and lower survivals. However, the recent diagnostic innovations (CINSARC gene-expression signature and tertiary lymphoid structure [TLS]) and therapeutic innovations (immune checkpoint inhibitors [ICIs]) for STS patients underscore the need to re-assess the role of 18F-FDG-PET/CT. Thus, in this correspondence, our objective was to investigate the correlations between STS metabolism as assessed by nuclear imaging, and the immune landscape as estimated by transcriptomics analysis, immunohistochemistry panels, and TLS assessment. Based on a prospective cohort of 85 adult patients with high-grade STS recruited in the NEOSARCOMICS trial (NCT02789384), we identified 3 metabolic groups according to 18F-FDG-PET/CT metrics (metabolic-low [60%], -intermediate [15.3%] and high [24.7%]). We found that T-cells CD8 pathway was significantly enriched in metabolic-high STS. Conversely, several pathways involved in antitumor immune response, cell differentiation and cell cycle, were downregulated in extreme metabolic-low STS. Next, multiplex immunofluorescence showed that densities of CD8+, CD14+, CD45+, CD68+, and c-MAF cells were significantly higher in the metabolic-high group compared to the metabolic-low group. Lastly, no association was found between metabolic group and TLS status. Overall, these results suggest that (i) rapidly proliferating and metabolically active STS can instigate a more robust immune response, thereby attracting immune cells such as T cells and macrophages, and (ii) metabolic activity and TLS could independently influence immune responses.

Differential gene expression (DGE) analysis in persons with a history of giardiasis

Giardiasis, which is caused by Giardia duodenalis, has clinical symptoms such as steatorrhea and can be very dangerous in children. In addition, some documents reported that this parasite is present inside the tissue of patients with cancer. In this study, we analyzed the gene expression profiles of some main genes important to apoptosis and anti-apoptosis in humans.Expression profile arrays of Genomic Spatial Event (GSE) 113666, GSE113667, and GSE113679 obtained from Gene Expression Omnibus were used for meta-analysis using R commands. Cytoscape and STRING databases used the protein-protein Interaction network. Then, the Kyoto Encyclopedia of Genes and Genomes and Gene Ontology analysis was performed. Similar genes in Homo sapiens were identified using Basic Local Alignment Search Tool analysis. The validation was performed on eight people using real-time Polymerase chain reaction. In addition to the candidate genes, the gene expression of some other genes, including Serine/Threonine Kinase 1 (AKT1), Cyclin Dependent Kinase Inhibitor 2A (CDKN2A), Kirsten Rat Sarcoma (KRAS), and Phosphatidylinositol-4,5-Bisphosphate 3-Kinase Catalytic Subunit Alpha (PIK3CA) were also examined. Analysis of the expression of serum amyloid A1 (SAA1), Regenerating Islet-Derived 3 Gamma (REG3G), and REG3A genes did not show any difference between the two groups of healthy and diseased people. Examining the mean expression of the four genes AKT1, CDKN2A, KRAS, and PIK3CA showed that three genes of AKT1, CDKN2A, and KRAS had increased expression in people with a history of giardiasis compared to healthy people. We showed that the gene expression pattern differs in apoptosis and anti-apoptosis signaling in people with a history of giardiasis. Giardia duodenalis seems to induce post-non-infectious symptoms with stimulation of human gene expression.

Inference of Dynamic Growth Regulatory Network in Cancer Using High-Throughput Transcriptomic Data

Growth is regulated by gene expression variation at different developmental stages of biological processes such as cell differentiation, disease progression, or drug response. In cancer, a stage-specific regulatory model constructed to infer the dynamic expression changes in genes contributing to tissue growth or proliferation is referred as a dynamic growth regulatory network (dGRN). Over the past decade, gene expression data has been widely used for reconstructing dGRN by computing correlations between the differentially expressed genes (DEGs). A wide variety of pipelines are available to construct the GRNs using DEGs and the choice of a particular method or tool depends on the nature of the study. In this protocol, we have outlined a step-by-step guide for the analysis of DEGs using RNA-Seq data, beginning from data acquisition, pre-processing, mapping to reference genome, and construction of a correlation-based co-expression network to further downstream analysis. We have also outlined the steps for the inclusion of publicly available interaction/regulation information into the dGRN followed by relevant topological inferences. This tutorial has been designed in a way that early researchers can refer to for an easy and comprehensive glimpse of methodologies used in the inference of dGRN using transcriptomics data.

Alteration of PBMC transcriptome profile after interaction with multipotent mesenchymal stromal cells under "physiological" hypoxia

Multipotent mesenchymal stromal cells (MSCs) have demonstrated a pronounced immunosuppressive activity, the manifestation of which depends on the microenvironmental factors, including O2 level. Here we examined the effects of MSCs on transcriptomic profile of allogeneic phytohemagglutinin-stimulated peripheral blood mononuclear cells (PBMCs) after interaction at ambient (20%) or "physiological" hypoxia (5%) O2. As revealed with microarray analysis, PBMC transcriptome at 20% O2 was more affected, which was manifested as differential expression of more than 300 genes, whereas under 5% O2 220 genes were changed. Most of genes at 20% O2 were downregulated, while at hypoxia most of genes were upregulated. Altered gene patterns were only partly overlapped at different O2 levels. A set of altered genes at hypoxia only was of particular interest. According to Gene Ontology a part of above genes was responsible for adhesion, cell communication, and immune response. At both oxygen concentrations, MSCs demonstrated effective immunosuppression manifested as attenuation of T cell activation and proliferation as well as anti-inflammatory shift of cytokine profile. Thus, MSC-mediated immunosuppression is executed with greater efficacy at a "physiological" hypoxia, since the same result has been achieved through a change in the expression of a fewer genes in target PBMCs.

Recent findings in molecular reactions of E. coli as exposed to alkylated, nano- and ordinary chitosans

The antibacterial effects of chitosan have been widely studied, but the underlying molecular mechanisms are not fully understood. We investigated the molecular responses of Escherichia coli MG1655 cell, a model gram-negative bacterium, upon exposure to chitosan (Cs), alkylated Cs (AlkCs), and chitosan nanoparticles (CsNPs). Nine target genes involved in relevant signaling pathways (ompF, ompC, ompA, mrcA, mrcB, mgtA, glnA, kdpA, lptA) were selected for analysis. A significant reduction in the expression of mrcA, mgtA, glnA, and lptA genes was observed in the cells treated with Cs. Those treated with Cs, AlkCs, and CsNPs revealed an increase in ompF gene expression, but the expression level was lower in the cells treated with AlkCs and CsNPs compared to Cs. This increase in porin expression suggests compromised membrane integrity and disrupted nutrient transport. In addition, the changes in the expression of mgtA, kdpA, and glnA are related to different effects on membrane permeability. The higher expression in the genes mrcA and mrcB is associated with morphological changes of cells treated with AlkCs and CsNPs. These findings contribute to our understanding of the molecular mechanisms underlying chitosan-induced stress responses and provide insights for the development of safer antimicrobial compounds in the future.

The forgotten variable? Does the euthanasia method and sample storage condition influence an organisms transcriptome - a gene expression analysis on multiple tissues in pigs

BACKGROUND

Transcriptomic studies often require collection of fresh tissues post euthanasia. The chosen euthanasia method might have the potential to induce variations in gene expressions that are unlinked with the experimental design. The present study compared the suitability of 'nitrogen gas in foam' (ANOXIA) in comparison to a non-barbiturate anaesthetic, T-61® (T61), for euthanizing piglets used in transcriptome research. Further, the effect of common tissue storage conditions, RNAlater™ (RL) and snap freezing in liquid nitrogen (LN2), on gene expression profiles were also analysed.

RESULTS

On comparison of the 3'mRNA-Seq data generated from pituitary, hypothalamus, liver and lung tissues, no significant differential expression in the protein coding genes were detected between the euthanasia methods. This implies that the nitrogen anoxia method could be a suitable alternative for euthanasia of piglets used in transcriptomic research. However, small nuclear RNAs (snRNAs) that constitute the eukaryotic spliceosomal machinery were found to be significantly higher (log2fold change ≥ 2.0, and adjusted p value ≤ 0.1) in pituitary samples collected using ANOXIA. Non-protein coding genes like snRNAs that play an important role in pre-mRNA splicing can subsequently modify gene expression. Storage in RL was found to be superior in preserving RNA compared to LN2 storage, as evidenced by the significantly higher RIN values in representative samples. However, storage in RL as opposed to LN2, also influenced differential gene expression in multiple tissues, perhaps as a result of its inability to inhibit biological activity during storage. Hence such external sources of variations should be carefully considered before arriving at research conclusions.

CONCLUSIONS

Source of biological variations like euthanasia method and storage condition can confound research findings. Even if we are unable to prevent the effect of these external factors, it will be useful to identify the impact of these variables on the parameter under observation and thereby prevent misinterpretation of our results.

A computational biology approach to identify potential protein biomarkers and drug targets for sporadic amyotrophic lateral sclerosis

Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disease caused by the loss of upper and lower motor neurons. The sporadic ALS (sALS) is a multigenic disorder and the complex mechanisms underlying its onset are still not fully delineated. Despite the recent scientific advancements, certain aspects of ALS pathogenic targets need to be yet clarified. The aim of the presented study is to identify potential genetic biomarkers and drug targets for sALS, by analysing gene expression profiles, presented in the publicly available GSE68605 dataset, of motor neurons cells obtained from sALS patients. We used different computational approaches including differential expression analysis, protein network mapping, candidate protein biomarker (CPB) identification, elucidation of the role of functional modules, and molecular docking analysis. The resultant top ten up- and downregulated genes were further used to construct protein-protein interaction network (PPIN). The PPIN analysis resulted in identifying four CPBs (namely RIOK2, AKT1, CTNNB1, and TNF) that commonly overlapped with one another in network parameters (degree, bottleneck and maximum neighbourhood component). The RIOK2 protein emerged as a potential mediator of top five functional modules that are associated with RNA binding, lipoprotein particle receptor binding in pre-ribosome, and interferon, cytokine-mediated signaling pathway. Furthermore, molecular docking analysis revealed that cyclosporine exhibited the highest binding affinity (-8.6 kJ/mol) with RIOK2, and surpassed the FDA-approved ALS drugs, such as riluzole and edaravone. This suggested that cyclosporine may serve as a promising candidate for targeting RIOK2 downregulation observed in sALS patients. In order to validate our computational results, it is suggested that in vitro and in vivo studies may be conducted in future to provide a more detailed understanding of ALS diagnosis, prognosis, and therapeutic intervention.

[Bioinformatics analysis of primary biliary cholangitis key genes and molecular mechanisms]

Objective: To extract the differentially expressed key genes of primary biliary cholangitis (PBC) using bioinformatics methods, so as to provide information for further study into the mechanism. Methods: The GSE119600 dataset was downloaded from the GEO database to obtain differentially expressed genes. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed for differentially expressed genes. Protein-protein interaction (PPI) network reconstruction, Cytoscape software visualization, and core gene screening were performed. The area under the receiver operating characteristic curve (ROC AUC) was used to assess the diagnostic effectiveness of genes and plot the pROC software package. The x-Cell software was used to calculate the enrichment score of 34 immune cells in each sample. Finally, four key genes (PSMA4, PSMA1, PSMB1, and PSMA3) were selected. Blood samples were analyzed using the qPCR method. Results:: A total of 373 immune-related differentially expressed genes were identified. Eight genes (PSMC6, PSMB2, PSMB1, PSMA3, PSMA4, PSMA1, PSMD7, and PSMB5) were screened from the 178 nodes and 596 edges as hub genes of the PPI network, which were significantly related to amino acid metabolism, hematopoietic stem cell differentiation, cell cycle, and immune processes. PSMA4, PSMA1, PSMB1, and PSMA3 were defined as immunological biomarkers for PBC with an AUC value of the ROC curve > 0.7. Immunoinfiltrating cell analysis showed that the proportion of eosinophils was significantly higher in PBC patients compared to the control group, whereas the proportion of CD4+ memory T cells, plasma cells, Th2 cells, and cDC cells was significantly lower in PBC patients than the control group. Plasma cells were associated with all four immunological biomarkers. Seven PBC patients and seven healthy subjects were selected for peripheral blood qPCR validation, which demonstrates that PSMB1, PSMA3, PSMA1, and PSMA4 levels were significantly lower in PBC patients than healthy subjects, with a statistically significant difference. Conclusion:: Bioinformatics screened eight key genes, of which four were key immunological markers and may serve as a basis for clinical diagnosis and mechanism exploration.

Key pathways and genes that are altered during treatment with hyperbaric oxygen in patients with sepsis due to necrotizing soft tissue infection (HBOmic study)

BACKGROUND

For decades, the basic treatment strategies of necrotizing soft tissue infections (NSTI) have remained unchanged, primarily relying on aggressive surgical removal of infected tissue, broad-spectrum antibiotics, and supportive intensive care. One treatment strategy that has been proposed as an adjunctive measure to improve patient outcomes is hyperbaric oxygen (HBO2) treatment. HBO2 treatment has been linked to several immune modulatory effects; however, investigating these effects is complicated due to the disease's acute life-threatening nature, metabolic and cell homeostasis dependent variability in treatment effects, and heterogeneity with respect to both patient characteristics and involved pathogens. To embrace this complexity, we aimed to explore the underlying biological mechanisms of HBO2 treatment in patients with NSTI on the gene expression level.

METHODS

We conducted an observational cohort study on prospective collected data, including 85 patients admitted to the intensive care unit (ICU) for NSTI. All patients were treated with one or two HBO2 treatments and had one blood sample taken before and after the intervention. Total RNAs from blood samples were extracted and mRNA purified with rRNA depletion, followed by whole-transcriptome RNA sequencing with a targeted sequencing depth of 20 million reads. A model for differentially expressed genes (DEGs) was fitted, and the functional aspects of the obtained set of genes was predicted with GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of genes and Genomes) enrichment analyses. All analyses were corrected for multiple testing with FDR.

RESULTS

After sequential steps of quality control, a final of 160 biological replicates were included in the present study. We found 394 protein coding genes that were significantly DEGs between the two conditions with FDR < 0.01, of which 205 were upregulated and 189 were downregulated. The enrichment analysis of these DEGs revealed 20 GO terms in biological processes and 12 KEGG pathways that were significantly overrepresented in the upregulated DEGs, of which the term; "adaptive immune response" (GO:0002250) (FDR = 9.88E-13) and "T cell receptor signaling pathway" (hsa04660) (FDR = 1.20E-07) were the most significant. Among the downregulated DEGs two biological processes were significantly enriched, of which the GO term "apoptotic process" (GO:0006915) was the most significant (FDR = 0.001), followed by "Positive regulation of T helper 1 cell cytokine production" (GO:2000556), and "NF-kappa B signaling pathway" (hsa04064) was the only KEGG pathway that was significantly overrepresented (FDR = 0.001).

CONCLUSIONS

When one or two sessions of HBO2 treatment were administered to patients with a dysregulated immune response and systemic inflammation due to NSTI, the important genes that were regulated during the intervention were involved in activation of T helper cells and downregulation of the disease-induced highly inflammatory pathway NF-κB, which was associated with a decrease in the mRNA level of pro-inflammatory factors.

TRIAL REGISTRATION

Biological material was collected during the INFECT study, registered at ClinicalTrials.gov (NCT01790698).

De novo transcriptional analysis of the response to starvation stress in the white ridgetail prawn, Exopalaemon carinicauda

To study the mechanism of the biomolecular response in Exopalaemon carinicauda to starvation stress, we subjected muscle tissue RNA samples from four stress points, including 0 d(control group), 10 d, 20 d, and 30 d, to starvation stress on white ridgetail prawn with a body weight of 1.41 + 0.42 g, aquaculture water temperature of 23-25 °C, salinity of 26, dissolved oxygen ≥5 mg/L, and pH 8-8.5, Then performed de novo transcriptome assembly and gene expression analysis using BGISEQ-500 with a tag-based digital gene expression (DGE) system. By de novo assembling at the four times, we obtained 28,167, 21,115, 24,497, and 27,080 reads, respectively. The results showed that the stress at 10 d led to no significant difference in the expressed genes, while the stress at 20 d and 30 d showed a significant increase (or decrease) in the expression of 97 (276) and 143 (410) genes, respectively, which were involved in 8 different metabolic pathways. In addition, we detected 2647 unigene transcription factors. Eleven upregulated and sixteen downregulated genes from the different starvation stress groups were choose to verify the reliability of the transcriptome data, and the results showed that the expression trends of these genes were consistent with the results shown by the transcriptome. The analysis of the experimental data and our discussion of the response mechanism of white ridgetail prawn under starvation stress provides a foundation for further screening of the key genes of starvation stress and may help to elucidate their functions.

Differential gene expression profiling implicates altered network development in rat postnatal day 4 cortex following 4-Methylimidazole (4-MeI) induced maternal seizures

Compromised maternal health leading to maternal seizures can have adverse effects on the healthy development of offspring. This may be the result of inflammation, hypoxia-ischemia, and altered GABA signaling. The current study examined cortical tissue from F2b (2nd litter of the 2nd generation) postnatal day 4 (PND4) offspring of female Harlan SD rats chronically exposed to the seizuregenic compound, 4-Methylimidazole (0, 750, or 2500 ppm 4-MeI). Maternal seizures were evident only at 2500 ppm 4-MeI. GABA related gene expression as examined by qRT-PCR and whole genome microarray showed no indication of disrupted GABA or glutamatergic signaling. Canonical pathway hierarchical clustering and multi-omics combinatory genomic (CNet) plots of differentially expressed genes (DEG) showed alterations in genes associated with regulatory processes of cell development including neuronal differentiation and synaptogenesis. Functional enrichment analysis showed a similarity of cellular processes across the two exposure groups however, the genes comprising each cluster were primarily unique rather than shared and often showed different directionality. A dose-related induction of cytokine signaling was indicated however, pathways associated with individual cytokine signaling were not elevated, suggesting an alternative involvement of cytokine signaling. Pathways related to growth process and cell signaling showed a negative activation supporting an interpretation of disruption or delay in developmental processes at the 2500 ppm 4-MeI exposure level with maternal seizures. Thus, while GABA signaling was not altered as has been observed with maternal seizures, the pattern of DEG suggested a potential for alteration in neuronal network formation.

Differential effects of follicle-stimulating hormone glycoforms on the transcriptome profile of cultured rat granulosa cells as disclosed by RNA-seq

It has been documented that variations in glycosylation on glycoprotein hormones, confer distinctly different biological features to the corresponding glycoforms when multiple in vitro biochemical readings are analyzed. We here applied next generation RNA sequencing to explore changes in the transcriptome of rat granulosa cells exposed for 0, 6, and 12 h to 100 ng/ml of four highly purified follicle-stimulating hormone (FSH) glycoforms, each exhibiting different glycosylation patterns: human pituitary FSH18/21 and equine FSH (eqFSH) (hypo-glycosylated), and human FSH24 and chinese-hamster ovary cell-derived human recombinant FSH (recFSH) (fully-glycosylated). Total RNA from triplicate incubations was prepared from FSH glycoform-exposed cultured granulosa cells obtained from DES-pretreated immature female rats, and RNA libraries were sequenced in a HighSeq 2500 sequencer (2 × 125 bp paired-end format, 10-15 × 106 reads/sample). The computational workflow focused on investigating differences among the four FSH glycoforms at three levels: gene expression, enriched biological processes, and perturbed pathways. Among the top 200 differentially expressed genes, only 4 (0.6%) were shared by all 4 glycoforms at 6 h, whereas 118 genes (40%) were shared at 12 h. Follicle-stimulating hormone glycocoforms stimulated different patterns of exclusive and associated up regulated biological processes in a glycoform and time-dependent fashion with more shared biological processes after 12 h of exposure and fewer treatment-specific ones, except for recFSH, which exhibited stronger responses with more specifically associated processes at this time. Similar results were found for down-regulated processes, with a greater number of processes at 6 h or 12 h, depending on the particular glycoform. In general, there were fewer downregulated than upregulated processes at both 6 h and 12 h, with FSH18/21 exhibiting the largest number of down-regulated associated processes at 6 h while eqFSH exhibited the greatest number at 12 h. Signaling cascades, largely linked to cAMP-PKA, MAPK, and PI3/AKT pathways were detected as differentially activated by the glycoforms, with each glycoform exhibiting its own molecular signature. These data extend previous observations demonstrating glycosylation-dependent differential regulation of gene expression and intracellular signaling pathways triggered by FSH in granulosa cells. The results also suggest the importance of individual FSH glycoform glycosylation for the conformation of the ligand-receptor complex and induced signalling pathways.

Integrated bioinformatic analysis reveals NOS2 as a novel ferroptosis-related biomarker for pre-eclampsia

BACKGROUND

Pre-eclampsia (PE) is a common condition in pregnancy; however, methods for early diagnosis and effective treatment options are lacking. Ferroptosis is a newly identified iron-dependent cell death pathway. The aim of this study was to investigate the role of ferroptosis-related genes in PE, the underlying mechanism, and their potential diagnostic value using a bioinformatics approach.

METHODS

We downloaded the GSE48424 and GSE98224 datasets from the Gene Expression Omnibus database. Differentially expressed genes (DEGs) between PE and healthy pregnancy samples were identified in the GSE48424 dataset and subjected to weighted gene co-expression network analysis; the most relevant modules were intersected with known ferroptosis-related genes to distinctly identify the role of ferroptosis in PE. We further searched transcription factors and microRNAs that are predicted to regulate these ferroptosis-related genes, and patients in the GSE48424 dataset were divided into two groups according to high or low expression of the key ferroptosis-related genes associated with PE. To obtain robust key ferroptosis-related genes in PE, we validated their expression levels in the external dataset GSE98224. Finally, the reverse transcription-quantitative polymerase chain reaction (RT-qPCR) assay was utilized to access the expression of these genes in the PE and normal blood samples.

RESULTS

Six ferroptosis-related genes involved in PE were obtained by overlapping 3661 genes most associated with PE, 565 DEGs between PE and normal samples, and 259 known ferroptosis-related genes. Among these genes, patients with PE displaying lower expression levels of NOS2 and higher expression levels of PTGS2 had a higher ferroptosis potential index. The expression pattern of NOS2 was consistent in the GSE48424 and GSE98224 datasets. RT-qPCR data confirmed that NOS2 expression was more significantly elevated in patients with PE than in those with a normal pregnancy.

CONCLUSIONS

Our study explored the diagnostic value of ferroptosis-related genes in PE, and identified NOS2 as the key gene linking ferroptosis and PE, suggesting a new candidate biomarker for early PE diagnosis.

Machine learning analysis of lung squamous cell carcinoma gene expression datasets reveals novel prognostic signatures

BACKGROUND

Lung squamous cell carcinoma (LUSC) patients are often diagnosed at an advanced stage and have poor prognoses. Thus, identifying novel biomarkers for the LUSC is of utmost importance.

METHODS

Multiple datasets from the NCBI-GEO repository were obtained and merged to construct the complete dataset. We also constructed a subset from this complete dataset with only known cancer driver genes. Further, machine learning classifiers were employed to obtain the best features from both datasets. Simultaneously, we perform differential gene expression analysis. Furthermore, survival and enrichment analyses were performed.

RESULTS

The kNN classifier performed comparatively better on the complete and driver datasets' top 40 and 50 gene features, respectively. Out of these 90 gene features, 35 were found to be differentially regulated. Lasso-penalized Cox regression further reduced the number of genes to eight. The median risk score of these eight genes significantly stratified the patients, and low-risk patients have significantly better overall survival. We validated the robust performance of these eight genes on the TCGA dataset. Pathway enrichment analysis identified that these genes are associated with cell cycle, cell proliferation, and migration.

CONCLUSION

This study demonstrates that an integrated approach involving machine learning and system biology may effectively identify novel biomarkers for LUSC.

Expanding the genetic toolbox of Rhodotorula toruloides by identification and validation of six novel promoters induced or repressed under nitrogen starvation

BACKGROUND

The non-conventional yeast Rhodotorula toruloides is an emerging host organism in biotechnology by merit of its natural capacity to accumulate high levels of carotenoids and intracellular storage lipids from a variety of carbon sources. While the number of genetic engineering strategies that employ R. toruloides is increasing, the lack of genetic tools available for modification of this yeast is still limiting strain development. For instance, several strong, constitutive R. toruloides promoters have been characterized, but to date, only five inducible promoters have been identified. Although nitrogen-limited cultivation conditions are commonly used to induce lipid accumulation in this yeast, no promoters regulated by nitrogen starvation have been described for R. toruloides.

RESULTS

In this study, we used a combination of genomics and transcriptomics methods to identify novel R. toruloides promoter sequences that are either inducible or repressible by nitrogen starvation. RNA sequencing was used to assess gene expression in the recently isolated strain R. toruloides BOT-A2 during exponential growth and during nitrogen starvation, when cultivated with either glucose or xylose as the carbon source. The genome of BOT-A2 was sequenced using a combination of long- and short-read sequencing and annotated with support of the RNAseq data. Differential expression analysis was used to identify genes with a |log2 fold change|≥ 2 when comparing their expression during nitrogen depletion to that during exponential growth. The promoter regions from 16 of these genes were evaluated for their ability to drive the expression of a fluorescent reporter gene. Three promoters that were clearly upregulated under nitrogen starvation and three that were downregulated were selected and further characterized. One promoter, derived from gene RTBOTA2_003877, was found to function like an on-off switch, as it was only upregulated under full nitrogen depletion and downregulated in the presence of the nitrogen source.

CONCLUSIONS

Six new R. toruloides promoters that were either upregulated or downregulated under nitrogen-starvation were identified. These substantially contribute to the available promoters when engineering this organism and are foreseen to be particularly useful for future engineering strategies requiring specific regulation of target genes in accordance with nitrogen availability.

Spatial transcriptomics of human placentas reveal distinct RNA patterns associated with morphology and preeclampsia

Spatial transcriptomics (ST) maps RNA level patterns within a tissue. This technology has not been previously applied to human placental tissue. We demonstrate analysis of human placental samples with ST. Unsupervised clustering revealed that distinct RNA patterns were found corresponding to different morphological structures. Additionally, when focusing upon terminal villi and hemoglobin associated structures, RNA levels differed between placentas from full term healthy pregnancies and those complicated by preeclampsia. The results from this study can provide a benchmark for future ST studies in placenta.

Differential expression of genes during recovery of Nicotiana tabacum from tomato leaf curl Gujarat virus infection

MAIN CONCLUSION

Nicotiana tabacum exhibits recovery response towards tomato leaf curl Gujarat virus. Transcriptome analysis revealed the differential expression of defense-related genes. Genes encoding for cysteine protease inhibitor, hormonal- and stress-related to DNA repair mechanism are found to be involved in the recovery process. Elucidating the role of host factors in response to viral infection is crucial in understanding the plant host-virus interaction. Begomovirus, a genus in the family Geminiviridae, is reported throughout the globe and is known to cause serious crop diseases. Tomato leaf curl Gujarat virus (ToLCGV) infection in Nicotiana tabacum resulted in initial symptom expression followed by a quick recovery in the systemic leaves. Transcriptome analysis using next-generation sequencing (NGS) revealed a large number of differentially expressed genes both in symptomatic as well as recovered leaves when compared to mock-inoculated plants. The virus infected N. tabacum results in alteration of various metabolic pathways, phytohormone signaling pathway, defense related protein, protease inhibitor, and DNA repair pathway. RT-qPCR results indicated that Germin-like protein subfamily T member 2 (NtGLPST), Cysteine protease inhibitor 1-like (NtCPI), Thaumatin-like protein (NtTLP), Kirola-like (NtKL), and Ethylene-responsive transcription factor ERF109-like (NtERTFL) were down-regulated in symptomatic leaves when compared to recovered leaves of ToLCGV-infected plants. In contrast, the Auxin-responsive protein SAUR71-like (NtARPSL) was found to be differentially down-regulated in recovered leaves when compared to symptomatic leaves and the mock-inoculated plants. Lastly, Histone 2X protein like (NtHH2L) gene was found to be down-regulated, whereas Uncharacterized (NtUNCD) was up-regulated in both symptomatic as well as recovered leaves compared to the mock-inoculated plants. Taken together, the present study suggests potential roles of the differentially expressed genes that might govern tobacco's susceptibility and/or recovery response towards ToLCGV infection.

Study of the cross-talk between Fasciola hepatica juveniles and the intestinal epithelial cells of the host by transcriptomics in an in vitro model

Fasciolosis is a globally widespread trematodiasis with a major economic and veterinary impact. Therefore, this disease is responsible for millions of dollars in losses to the livestock industry, and also constitutes an emerging human health problem in endemic areas. The ubiquitous nature of Fasciola hepatica, the main causative agent, is one of the key factors for the success of fasciolosis. Accordingly, this parasite is able to subsist in a wide variety of ecosystems and hosts, thanks to the development of a plethora of strategies for adaption and immune evasion. Fasciolosis comprises a growing concern due to its high prevalence rates, together with the emergence of strains of the parasite resistant to the treatment of choice (triclabendazole). These facts highlight the importance of developing novel control measures which allow for an effective protection against the disease before F. hepatica settles in a niche inaccessible to the immune system. However, knowledge about the initial phases of the infection, including the migration mechanisms of the parasite and the early innate host response, is still scarce. Recently, our group developed an in vitro host-parasite interaction model that allowed the early events to be unveiled after the first contact between the both actors. This occurs shortly upon ingestion of F. hepatica metacercariae and the emergence of the newly excysted juveniles (FhNEJ) in the host duodenum. Here, we present a transcriptomic analysis of such model using an approach based on RNA sequencing (RNA-Seq), which reveals changes in gene expression related to proteolysis and uptake of metabolites in FhNEJ. Additionally, contact with the parasite triggered changes in host intestinal cells related to pseudogenes expression and host defence mechanisms, including immune response, among others. In sum, these results provide a better understanding of the early stages of fasciolosis at molecular level, and a pool of targets that could be used in future therapeutic strategies against the disease.

Preliminary transcriptomic analysis of peripheral blood from German Shepherd dogs with degenerative joint disease for the identification of diagnostic biomarkers

OBJECTIVE

Degenerative joint disease (DJD) includes a group of disorders characterised by the deterioration of the articular cartilage. In this study, we investigated the transcriptomic profile of peripheral blood in German Shepherd dogs with DJD to identify putative diagnostic biomarkers.

METHODS

Differential gene expression (DGE) and gene ontology (GO) analyses of the bulk RNA-seq experiment were performed in a cohort of 12 adult dogs (five cases and seven controls, classified by clinical and radiographic analyses).

RESULTS

Radiographs of cases revealed severe signs of progressive DJD. Two up-regulated (LOC106559672 and THBS4) and one down-regulated (LOC106559235) differentially expressed genes (adjusted p value < 0.05) were identified. The DGE with log2 fold change < -1.5 and > 1.5 and non-adjusted p < 0.01 were selected for GO analysis. No significant enrichment terms were observed in the selected threshold.

CONCLUSION

The gene-encoding protein THBS4 is correlated with DJD severity and long noncoding RNA LOC106559235 is probably involved in the DJD process. The THBS4 gene should be considered a good biomarker for DJD in dogs. Future studies using independent cohorts will be necessary to validate the present results.

Physiological response and transcriptome analyses of leguminous Indigofera bungeana Walp. to drought stress

Objective

Indigofera bungeana is a shrub with high quality protein that has been widely utilized for forage grass in the semi-arid regions of China. This study aimed to enrich the currently available knowledge and clarify the detailed drought stress regulatory mechanisms in I. bungeana, and provide a theoretical foundation for the cultivation and resistance breeding of forage crops.

Methods

This study evaluates the response mechanism to drought stress by exploiting multiple parameters and transcriptomic analyses of a 1-year-old seedlings of I. bungeana in a pot experiment.

Results

Drought stress significantly caused physiological changes in I. bungeana. The antioxidant enzyme activities and osmoregulation substance content of I. bungeana showed an increase under drought. Moreover, 3,978 and 6,923 differentially expressed genes were approved by transcriptome in leaves and roots. The transcription factors, hormone signal transduction, carbohydrate metabolism of regulatory network were observed to have increased. In both tissues, genes related to plant hormone signaling transduction pathway might play a more pivotal role in drought tolerance. Transcription factors families like basic helix-loop-helix (bHLH), vian myeloblastosis viral oncogene homolog (MYB), basic leucine zipper (bZIP) and the metabolic pathway related-genes like serine/threonine-phosphatase 2C (PP2C), SNF1-related protein kinase 2 (SnRK2), indole-3-acetic acid (IAA), auxin (AUX28), small auxin up-regulated rna (SAUR), sucrose synthase (SUS), sucrosecarriers (SUC) were highlighted for future research about drought stress resistance in Indigofera bungeana.

Conclusion

Our study posited I. bungeana mainly participate in various physiological and metabolic activities to response severe drought stress, by regulating the expression of the related genes in hormone signal transduction. These findings, which may be valuable for drought resistance breeding, and to clarify the drought stress regulatory mechanisms of I. bungeana and other plants.

The metabolism and detoxification effects of lead exposure on Aleurolyphus ovatus (Acari: Acaridae) via transcriptome analysis

Aleurolyphus ovatus Troupeau is one of the most predominant species of the Acaridae family worldwide. Recent reports have demonstrated that the accumulation of lead in stored grains and dietary items exceeds the required standards. However, the molecular mechanism of heavy metal stress on mites has not been reported. To understand the mechanism underlying the heavy metal response of A. ovatus, comparative transcriptome analysis was performed in this study using an Illumina high throughput mRNA sequencing (RNA-seq) platform. A. ovatus was fed on artificial diets containing two different concentrations of lead, namely, a low concentration of 12.5 mg/kg (LAO) and a high concentration of 100 mg/kg (HAO), while the mites in the control (NAO) group were not exposed to lead. A total of 44,362 unigenes, with an average length of 1547 bp, were identified. Of these, 996 unigenes were successfully annotated in seven functional databases. The number of differentially expressed genes (DEGs) in A. ovatus under different lead concentrations was compared. In NAO versus LAO group, including 310 up-regulated and 1580 down-regulated DEGs. In NAO versus HAO group, including 3928 up-regulated and 1761 down-regulated DEGs. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment indicated that detoxification enzyme genes were significantly expressed in pathways, such as cytochrome P450 foreign body metabolism, glutathione metabolism and drug metabolism-cytochrome pathway. The results of gene annotation and quantitative real-time PCR showed that high concentration of lead significantly stimulated the expression of metabolic detoxification enzyme genes such as glutathione S transferase (GST) and superoxide dismutase (SOD), while low concentration inhibited their expression. This study will provide a basis for the molecular mechanism of A. ovatus in response to heavy metal lead stimulation in stored grain.

Robust identification of regulatory variants (eQTLs) using a differential expression framework developed for RNA-sequencing

BACKGROUND

A gap currently exists between genetic variants and the underlying cell and tissue biology of a trait, and expression quantitative trait loci (eQTL) studies provide important information to help close that gap. However, two concerns that arise with eQTL analyses using RNA-sequencing data are normalization of data across samples and the data not following a normal distribution. Multiple pipelines have been suggested to address this. For instance, the most recent analysis of the human and farm Genotype-Tissue Expression (GTEx) project proposes using trimmed means of M-values (TMM) to normalize the data followed by an inverse normal transformation.

RESULTS

In this study, we reasoned that eQTL analysis could be carried out using the same framework used for differential gene expression (DGE), which uses a negative binomial model, a statistical test feasible for count data. Using the GTEx framework, we identified 35 significant eQTLs (P < 5 × 10^-8) following the ANOVA model and 39 significant eQTLs (P < 5 × 10^-8) following the additive model. Using a differential gene expression framework, we identified 930 and six significant eQTLs (P < 5 × 10^-8) following an analytical framework equivalent to the ANOVA and additive model, respectively. When we compared the two approaches, there was no overlap of significant eQTLs between the two frameworks. Because we defined specific contrasts, we identified trans eQTLs that more closely resembled what we expect from genetic variants showing complete dominance between alleles. Yet, these were not identified by the GTEx framework.

CONCLUSIONS

Our results show that transforming RNA-sequencing data to fit a normal distribution prior to eQTL analysis is not required when the DGE framework is employed. Our proposed approach detected biologically relevant variants that otherwise would not have been identified due to data transformation to fit a normal distribution.

Transcriptomic profiling of near-isogenic lines reveals candidate genes for a significant locus conferring metribuzin resistance in wheat

BACKGROUND

Weeds reduce wheat yields in dryland farming systems. Herbicides such as metribuzin are commonly used to control weeds. However, wheat has a narrow safety margin against metribuzin. Standing crops such as wheat with weeds in the same field can also be killed by the same dose of metribuzin. Therefore, it is important to identify metribuzin resistance genes and understand the resistance mechanism in wheat for sustainable crop production. A previous study identified a significant metribuzin resistance wheat QTL, Qsns.uwa.4 A.2, explaining 69% of the phenotypic variance for metribuzin resistance.

RESULTS

Two NIL pairs with the most contrasting performance in the metribuzin treatment and different in genetic backgrounds were compared using RNA sequence analysis, identifying nine candidate genes underlying Qsns.uwa.4 A.2 responsible for metribuzin resistance. Quantitative RT-qPCR further validated the candidate genes, with TraesCS4A03G1099000 (nitrate excretion transporter), TraesCS4A03G1181300 (aspartyl protease), and TraesCS4A03G0741300 (glycine-rich proteins) identified as key factors for metribuzin resistance.

CONCLUSION

Identified markers and key candidate genes can be used for selecting metribuzin resistance in wheat.

Sex differences in gene expression and alternative splicing in the Chinese horseshoe bat

Sexually dimorphic traits are common in sexually reproducing organisms and can be encoded by differential gene regulation between males and females. Although alternative splicing is common mechanism in generating transcriptional diversity, its role in generating sex differences relative to differential gene expression is less clear. Here, we investigate the relative roles of differential gene expression and alternative splicing between male and female the horseshoe bat species, Rhinolophus sinicus. Horseshoe bats are an excellent model to study acoustic differences between sexes. Using RNA-seq analyses of two somatic tissues (brain and liver) from males and females of the same population, we identified 3,471 and 2,208 differentially expressed genes between the sexes (DEGs) in the brain and liver, respectively. DEGs were enriched with functional categories associated with physiological difference of the sexes (e.g.,gamete generation and energy production for reproduction in females). In addition, we also detected many differentially spliced genes between the sexes (DSGs, 2,231 and 1,027 in the brain and liver, respectively) which were mainly involved in regulation of RNA splicing and mRNA metabolic process. Interestingly, we found a significant enrichment of DEGs on the X chromosome, but not for DSGs. As for the extent of overlap between the two sets of genes, more than expected overlap of DEGs and DSGs was observed in the brain but not in the liver. This suggests that more complex tissues, such as the brain, may require the intricate and simultaneous interplay of both differential gene expression and splicing of genes to govern sex-specific functions. Overall, our results support that variation in gene expression and alternative splicing are important and complementary mechanisms governing sex differences.

Novel evolved Yarrowia lipolytica strains for enhanced growth and lipid content under high concentrations of crude glycerol

BACKGROUND

Yarrowia lipolytica is a well-studied oleaginous yeast known for its ability to accumulate and store intracellular lipids, while growing on diverse, non-conventional substrates. Amongst them, crude glycerol, a low-cost by-product of the biodiesel industry, appears to be an interesting option for scaling up a sustainable single-cell oil production process. Adaptive laboratory evolution (ALE) is a powerful tool to force metabolic adaptations endowing tolerance to stressful environmental conditions, generating superior phenotypes with industrial relevance.

RESULTS

Y. lipolytica MUCL 28849 underwent ALE in a synthetic medium with increasing concentration of pure or crude glycerol as a stressing factor (9-20% v/v) for 520 generations. In one case of pure glycerol, chemical mutagenesis with ethyl methanesulfonate (EMS) was applied prior to ALE. Growth profile, biomass production and lipid content of 660 evolved strains (EVS), revealed 5 superior isolates; exhibiting from 1.9 to 3.6-fold increase of dry biomass and from 1.1 to 1.6-fold increase of lipid concentration compared to the parental strain, when grown in 15% v/v crude glycerol. NGS for differential gene expression analysis, showed induced expression in all EVS affecting nucleosomal structure and regulation of transcription. As strains differentiated, further changes accumulated in membrane transport and protein transport processes. Genes involved in glycerol catabolism and triacylglycerol biosynthesis were overexpressed in two EVS. Mismatches and gaps in the expressed sequences identified altered splicing and mutations in the EVS, with most of them, affecting different components of septin ring formation in the budding process. The selected YLE155 EVS, used for scale-up cultivation in a 3L benchtop bioreactor with 20% v/v crude glycerol, achieved extended exponential phase, twofold increase of dry biomass and lipid yields at 48 h, while citric acid secretion and glycerol consumption rates were 40% and 50% lower, respectively, compared to the parental strain, after 24 h of cultivation.

CONCLUSION

ALE and EMS-ALE under increasing concentrations of pure or crude glycerol generated novel Y. lipolytica strains with enhanced biomass and lipid content. Differential gene expression analysis and scale-up of YLE155, illustrated the potential of the evolved strains to serve as suitable "chassis" for rational engineering approaches towards both increased lipid accumulation, and production of high-added value compounds, through efficient utilization of crude glycerol.

Functional genome annotation and transcriptome analysis of Pseudozyma hubeiensis BOT-O, an oleaginous yeast that utilizes glucose and xylose at equal rates

Pseudozyma hubeiensis is a basidiomycete yeast that has the highly desirable traits for lignocellulose valorisation of being equally efficient at utilization of glucose and xylose, and capable of their co-utilization. The species has previously mainly been studied for its capacity to produce secreted biosurfactants in the form of mannosylerythritol lipids, but it is also an oleaginous species capable of accumulating high levels of triacylglycerol storage lipids during nutrient starvation. In this study, we aimed to further characterize the oleaginous nature of P. hubeiensis by evaluating metabolism and gene expression responses during storage lipid formation conditions with glucose or xylose as a carbon source. The genome of the recently isolated P. hubeiensis BOT-O strain was sequenced using MinION long-read sequencing and resulted in the most contiguous P. hubeiensis assembly to date with 18.95 Mb in 31 contigs. Using transcriptome data as experimental support, we generated the first mRNA-supported P. hubeiensis genome annotation and identified 6540 genes. 80% of the predicted genes were assigned functional annotations based on protein homology to other yeasts. Based on the annotation, key metabolic pathways in BOT-O were reconstructed, including pathways for storage lipids, mannosylerythritol lipids and xylose assimilation. BOT-O was confirmed to consume glucose and xylose at equal rates, but during mixed glucose-xylose cultivation glucose was found to be taken up faster. Differential expression analysis revealed that only a total of 122 genes were significantly differentially expressed at a cut-off of |log₂ fold change| ≥ 2 when comparing cultivation on xylose with glucose, during exponential growth and during nitrogen-starvation. Of these 122 genes, a core-set of 24 genes was identified that were differentially expressed at all time points. Nitrogen-starvation resulted in a larger transcriptional effect, with a total of 1179 genes with significant expression changes at the designated fold change cut-off compared with exponential growth on either glucose or xylose.

Differential gene expression associated with flower development of mango (Mangifera indica L.) varieties with different shelf-life

Mango (Mangifera indica L.) is one of the most important commercial fruit crop grown in many parts of the world. Major challenges affecting mango trade are short shelf-life, high susceptibility to chilling injury, post-harvest diseases and consumer demand for improved fruit quality. The objective of the present study was to reveal the key regulators present in bud and flower tissues during flower development stage, associated with fruit development and affect the shelf-life of the mango fruit. RNA-sequencing of contrasting genotypes having short and long shelf-life, was carried out. Comparative differential expression pathway studies of long shelf-life (Totapuri) and short shelf-life (Bombay Green) mango genotypes revealed a total of 177 highly differentially expressed genes. Out of 177 total genes, 101 genes from endoplasmic reticulum pathway and very few from gibberellins (3) and jasmonic acid (1) pathway were identified. Genes from endoplasmic reticulum pathway like hsp 90, SRC2, DFRA, CHS, BG3 and ASPG1 mainly up regulated in Bombay Green. Uniprotein B9R8D3 also shows up regulation in Bombay Green. Ethylene insensitive pathway gene EIL1 up regulated in Bombay Green. Gene CAD1 from phenylpropanoid pathway mainly up regulated in Bombay Green. A total of 4 SSRs and 227 SNPs were mined from these pathways specific to the shelf-life. Molecular studies of endoplasmic reticulum, phenylpropanoid, ethylene, polygalacturonase and hormone pathways at the time of bud and flower formation revealed key regulators that determine the shelf-life of mango fruit.

Transcriptome Analysis Revealed Osmoregulation Related Regulatory Networks and Hub Genes in the Gills of Hilsa shad, Tenualosa ilisha, during the Migratory Osmotic Stress

Tenualosa ilisha (Hilsa shad), an anadromous fish, usually inhabits coastal and estuarine waters, and migrates to freshwater for spawning. In this study, large-scale gill transcriptome analyses from three salinity regions, i.e., fresh, brackish and marine water, revealed 3277 differentially expressed genes (DEGs), out of which 232 were found to be common between marine vs freshwater and brackish vs freshwater. These genes were mapped into 54 KEGG Pathways, and the most significant of these were focal adhesion, adherens junction, tight junction, and PI3K-Akt signaling pathways. A total of 24 osmoregulatory genes were found to be differentially expressed in different habitats. The gene members of slc16 and slc2 families showed a dissimilar pattern of expressions, while two claudin genes (cldn11 & cldn10), transmembrane tm56b, and voltage-gated potassium channel gene kcna10 were downregulated in freshwater samples, as compared to that of brackish and marine environment. Protein-protein interaction (PPI) network analysis of 232 DEGs showed 101 genes to be involved in PPI, while fn1 gene was found to be interacting with the highest number of genes (36). Twenty-five hub genes belonged to 12 functional groups, with muscle structure development with seven genes, forming the major group. These results provided valuable information about the genes, potentially involved in the molecular mechanisms regulating water homeostasis in gills, during migration for spawning and low-salinity adaptation in Hilsa shad. These genes may form the basis for the bio-marker development for adaptation to the stress levied by major environmental changes, due to hatchery/culture conditions.

Start codon targeted (SCoT) polymorphism marker in plant genome analysis: current status and prospects

The present review illustrates a comprehensive overview of the start codon targeted (SCoT) polymorphism marker and their utilization in various applications related to genetic and genomic studies. Start codon targeted (SCoT) polymorphism marker, a targeted fingerprinting marker technique, has gained considerable importance in plant genetics, genomics, and molecular breeding due to its many desirable features. SCoT marker targets the region flanking the start codon, a highly conserved region in plant genes. Therefore, it can distinguish genetic variations in a specific gene that link to a specific trait. It is a simple, novel, cost-effective, highly polymorphic, and reproducible molecular marker for which there is no need for prior sequence information. In the recent past, SCoT markers have been employed in many commercially important and underutilized plant species for a variety of applications, including genetic diversity analysis, interspecific/generic genetic relationships, cultivar/hybrid/species identification, sex determination, construction of linkage map, association mapping/analysis, differential gene expression, and genetic fidelity analysis of tissue culture-raised plants. The main aim of this review is to provide up-to-date information on SCoT markers and their application in many commercially important and underutilized plant species, mainly progress made in the last 8-10 years.

Differential Expression Feature Extraction (DEFE): A Case Study in Wheat FHB RNA-Seq Data Analysis

In differential gene expression data analysis, one objective is to identify groups of co-expressed genes from a large dataset in order to detect the association between such a group of genes and an experimental condition. This is often done through a clustering approach, such as k-means or bipartition hierarchical clustering, based on particular similarity measures in the grouping process. In such a dataset, the gene differential expression itself is an innate attribute that can be used in the feature extraction process. For example, in a dataset consisting of multiple treatments versus their controls, the expression of a gene in each treatment would have three possible behaviors, upregulated, downregulated, or unchanged. We present in this chapter, a differential expression feature extraction (DEFE) method by using a string consisting of three numerical values at each character to denote such behavior, i.e., 1 = up, 2 = down, and 0 = unchanged, which results in up to 3^B differential expression patterns across all B comparisons. This approach has been successfully applied in many research projects, and among these, we demonstrate the strength of DEFE in a case study on RNA-sequencing (RNA-seq) data analysis of wheat challenged with the phytopathogenic fungus, Fusarium graminearum. Combinations of multiple schemes of DEFE patterns revealed groups of genes putatively associated with resistance or susceptibility to FHB.

Possible Contribution of Alternative Transcript Isoforms in Mature Biofilm Growth Phase of Candida glabrata

Expression of genome-wide alternative transcript isoforms and differential transcript isoform usage in different biological conditions (isoform switching) are responsible for the varied proteomic functional diversity in higher eukaryotic organisms. However, these mechanisms have not been studied in Candida glabrata, which is a potent eukaryotic opportunistic pathogen. Biofilm formation is an important virulence factor of C. glabrata that withstands antifungal drug stress and overcomes the host-immune response. Here, we present the genome-wide differential transcript isoform expression (DTE) and differential transcript isoform usage (DTU) in a mature biofilm growth phase of C. glabrata (clinical isolate; NCCPF 100,037) using the RNA sequencing approach. The DTE analysis generated 7837 transcript isoforms from the C. glabrata genome (5293 genes in total), and revealed that transcript isoforms generated from 292 genes showed significant DTU in the mature biofilm cells. Gene ontology, pathway analysis and protein-protein interactions of significant transcript isoforms, further substantiated that their specific expression and differential usage is required for transitioning the planktonic cells to biofilm in C. glabrata. The present study reported the possible role of expression of alternative transcript isoforms and differential transcript isoform usage in the mature biofilms of C. glabrata. The observation derived from the study may prove to be beneficial for making future antifungal therapeutic strategies.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12088-022-01036-7.