Transcriptome Data Research Articles

Starch and sucrose metabolism plays an important role in the stem development in Medicago sativa

The forage quality of alfalfa (Medicago sativa) stems is greater than the leaves. Sucrose hydrolysis provides energy for stem development, with starch being enzymatically converted into sucrose to maintain energy homeostasis. To understand the physiological and molecular networks controlling stem development, morphological characteristics and transcriptome profiles in the stems of two alfalfa cultivars (Zhungeer and WL168) were investigated. Based on transcriptome data, we analysed starch and sugar contents, and enzyme activity related to starch-sugar interconversion. Zhungeer stems were shorter and sturdier than WL168, resulting in significantly higher mechanical strength. Transcriptome analysis showed that starch and sucrose metabolism were significant enriched in the differentially expressed genes of stems development in both cultivars. Genes encoding INV, bglX, HK, TPS and glgC downregulated with the development of stems, while the gene encoding was AMY upregulated. Weighted gene co-expression network analysis revealed that the gene encoding glgC was pivotal in determining the variations in starch and sucrose contents between the two cultivars. Soluble carbohydrate, sucrose, and starch content of WL168 were higher than Zhungeer. Enzyme activities related to sucrose synthesis and hydrolysis (INV, bglX, HK, TPS) showed a downward trend. The change trend of enzyme activity was consistent with gene expression. WL168 stems had higher carbohydrate content than Zhungeer, which accounted for more rapid growth and taller plants. WL168 formed hollow stems were formed during rapid growth, which may be related to the redistribution of carbohydrates in the pith tissue. These results indicated that starch and sucrose metabolism play important roles in the stem development in alfalfa.

Nanoformulation of dasatinib cannot overcome therapy resistance of pancreatic cancer cells with low LYN kinase expression.

Pancreatic ductal adenocarcinoma (PDAC) is one of the most difficult to treat tumors. The Src (sarcoma) inhibitor dasatinib (DASA) has shown promising efficacy in preclinical studies of PDAC. However, clinical confirmation could not be achieved. Overall, our aim was to deliver arguments for the possible reinitiating clinical testing of this compound in a biomarker-stratifying therapy trial for PDAC patients. We tested if the nanofunctionalization of DASA can increase the drug efficacy and whether certain Src members can function as clinical predictive biomarkers. Methods include manufacturing of poly(vinyl alcohol) stabilized gold nanoparticles and their drug loading, dynamic light scattering, transmission electron microscopy, thermogravimetric analysis, Zeta potential measurement, sterile human cell culture, cell growth quantification, accessing and evaluating transcriptome and clinical data from molecular tumor dataset TCGA, as well as various statistical analyses. We generated homo-dispersed nanofunctionalized DASA as an AuNP@PVA-DASA conjugate. The composite did not enhance the anti-growth effect of DASA on PDAC cell lines. The cell model with high LYN expression showed the strongest response to the therapy. We confirm deregulated Src kinetome activity as a prevalent feature of PDAC by revealing mRNA levels associated with higher malignancy grade of tumors. BLK (B lymphocyte kinase) expression predicts shorter overall survival of diabetic PDAC patients. Nanofunctionalization of DASA needs further improvement to overcome the therapy resistance of PDAC. LYN mRNA is augmented in tumors with higher malignancy and can serve as a predictive biomarker for the therapy resistance of PDAC cells against DASA. Studying the biological roles of BLK might help to identify underlying molecular mechanisms associated with PDAC in diabetic patients.

Coevolution of Non-Homologous End Joining efficiency and encephalization.

Double strand breaks, the most difficult to repair DNA damage, are mainly repaired by Non-Homologous End-Joining (NHEJ) or Homologous Recombination (HR). Previous studies seem to indicate that primates, and particularly humans, have a better NHEJ system. A distinctive feature of the primate lineage (beside longevity) is encephalization, i.e. the expansion of brain relative to body mass. Using existing transcriptome data from 34 mammalian species, we investigated the possible correlations between the expression of genes involved in NHEJ and encephalization, body mass and longevity. The same was done also for genes involved in the HR pathway. We found that, while HR gene expression is better correlated with longevity, NHEJ gene expression is strongly (and better) correlated with encephalization. Since the brain is composed of post-mitotic cells, double-strand breaks repair should be mainly performed by NHEJ in this organ. Therefore, we interpret the correlation we found as an indication that NHEJ efficiency coevolved with encephalization.

Combined transcriptomics and proteomics unveil the impact of vitamin C in modulating specific protein abundance in the mouse liver

BackgroundVitamin C (ascorbate) is a water-soluble antioxidant and an important cofactor for various biosynthetic and regulatory enzymes. Mice can synthesize vitamin C thanks to the key enzyme gulonolactone oxidase (Gulo) unlike humans. In the current investigation, we used Gulo−/− mice, which cannot synthesize their own ascorbate to determine the impact of this vitamin on both the transcriptomics and proteomics profiles in the whole liver. The study included Gulo−/− mouse groups treated with either sub-optimal or optimal ascorbate concentrations in drinking water. Liver tissues of females and males were collected at the age of four months and divided for transcriptomics and proteomics analysis. Immunoblotting, quantitative RT-PCR, and polysome profiling experiments were also conducted to complement our combined omics studies.ResultsPrincipal component analyses revealed distinctive differences in the mRNA and protein profiles as a function of sex between all the mouse cohorts. Despite such sexual dimorphism, Spearman analyses of transcriptomics data from females and males revealed correlations of hepatic ascorbate levels with transcripts encoding a wide array of biological processes involved in glucose and lipid metabolisms as well as in the acute-phase immune response. Moreover, integration of the proteomics data showed that ascorbate modulates the abundance of various enzymes involved in lipid, xenobiotic, organic acid, acetyl-CoA, and steroid metabolism mainly at the transcriptional level, especially in females. However, several proteins of the mitochondrial complex III significantly correlated with ascorbate concentrations in both males and females unlike their corresponding transcripts. Finally, poly(ribo)some profiling did not reveal significant enrichment difference for these mitochondrial complex III mRNAs between Gulo−/− mice treated with sub-optimal and optimal ascorbate levels.ConclusionsThus, the abundance of several subunits of the mitochondrial complex III are regulated by ascorbate at the post-transcriptional levels. Our extensive omics analyses provide a novel resource of altered gene expression patterns at the transcriptional and post-transcriptional levels under ascorbate deficiency.

A multilevel investigation to reveal the regulatory mechanism of lignin accumulation in juice sac granulation of pomelo

Granulation of juice sacs is a physiological disorder, which affects pomelo fruit quality. Here, the transcriptome and ubiquitinome of the granulated juice sacs were analyzed in Guanxi pomelo. We found that lignin accumulation in the granulated juice sacs was regulated at transcription and protein modification levels. In transcriptome data, we found that the genes in lignin biosynthesis pathway and antioxidant enzyme system of the granulated juice sacs were significantly upregulated. However, in ubiquitinome data, we found that ubiquitinated antioxidant enzymes increased in abundance but the enzyme activities decreased after the modification, which gave rise to reactive oxygen species (ROS) contents in granulated juice sacs. This finding suggests that ubiquitination level of the antioxidant enzymes is negatively correlated with the enzyme activities. Increased H2O2 is considered to be a signaling molecule to activate the key gene expressions in lignin biosynthesis pathway, which leads to the lignification in granulated juice sacs of pomelo. This regulatory mechanism in juice sac granulation of pomelo was further confirmed through the verification experiment using tissue culture by adding H2O2 or dimethylthiourea (DMTU). Our findings suggest that scavenging H2O2 and other ROS are important for reducing lignin accumulation, alleviating juice sac granulation and improving pomelo fruit quality.

A novel disulfide death-related genes prognostic signature identifies the role of IPO4 in glioma progression

Background“Disulfide death,” a form of cellular demise, is triggered by the abnormal accumulation of intracellular disulfides under conditions of glucose deprivation. However, its role in the prognosis of glioma remains undetermined. Therefore, the main objective of this study is to establish prognostic signature based on disulfide death-related genes (DDRGs) and to provide new solutions in choosing the effective treatment of glioma.MethodsThe RNA transcriptome, clinical information, and mutation data of glioma samples were sourced from The Cancer Genome Atlas (TCGA) and the Chinese Glioma Genome Atlas (CGGA), while normal samples were obtained from the Genotype-Tissue Expression (GTEx). DDRGs were compiled from previous studies and selected through differential analysis and univariate Cox regression analysis. The molecular subtypes were determined through consensus clustering analysis. Further, LASSO analysis was employed to select characteristic genes, and subsequently, a risk model comprising seven DDRGs was constructed based on multivariable Cox analysis. Kaplan-Meier survival curves were employed to assess survival differences between high and low-risk groups. Additionally, functional analyses (GO, KEGG, GSEA) were conducted to explore the potential biological functions and signaling pathways of genes associated with the model. The study also explored immune checkpoint (ICP) genes, immune cell infiltration levels, and immune stromal scores. Finally, the effect of Importin-4(IPO4) on glioma has been further confirmed through RT-qPCR, Western blot, and cell functional experiments.Results7 genes associated with disulfide death were obtained and two subgroups of patients with different prognosis and clinical characteristics were identified. Risk signature was subsequently developed and proved to serve as an prognostic predictor. Notably, the high-risk group exhibited an immunosuppressive microenvironment characterized by a high concentration of M2 macrophages and regulatory T cells (Tregs). In contrast, the low-risk group showed lower half-maximal inhibitory concentration (IC50) values. Therefore, patients in the high-risk group may benefit more from immunotherapy, while patients in the low-risk group may benefit more from chemotherapy. In addition, in vitro experiments have shown that inhibition of the expression of IPO4 leads to a significant reduction in the proliferation, migration, and invasion of glioma cells.ConclusionThis study identified two glioma subtypes and constructed a prognostic signature based on DDRGs. The signature has the potential to optimize the selection of patients for immune- and chemotherapy and provided a potential therapeutic target for glioma.

Retraction Note: Identification of biomarkers of hepatocellular carcinoma gene prognosis based on the immune-related lncRNA signature of transcriptome data.

Retraction Note: Identification of biomarkers of hepatocellular carcinoma gene prognosis based on the immune-related lncRNA signature of transcriptome data.

Molecular evolution of transcription factors AF4/FMR2 family member (AFF) gene family and the role of lamprey AFF3 in cell proliferation.

AF4/FMR2 family member (AFF) proteins are a group of transcriptional regulators that can regulate gene transcription and play an important role in cellular physiological processes such as proliferation and differentiation. The transcriptome data of the lamprey spinal cord injury were analyzed in previous research. We then identified a hub gene, Lr-AFF3, from this dataset. Phylogenetic tree analysis determined the evolutionary relationships of the AFF gene family across different species. In addition, analysis of motifs, domains, and 3D structures further confirmed the conservatism of the AFF gene family. In particular, the gene structure of the AFF3 gene was not conserved, possibly because of intron insertion. It was also found that the neighboring genes of the Lr-AFF3 gene had a higher diversity than that in jawed vertebrates through synteny analysis. The results of the MTT and EdU experiments showed that the C-terminal homology domain (CHD) and N-terminal homology domain (NHD) of Lr-AFF3 promoted cell proliferation. In summary, our research will not only provide new insights into the origin and evolution of the AFF gene family in different species, but also provide new clues for the functions of Lr_AFF3.

Expression of Anoikis-Related Genes and Potential Biomarkers in Colon Cancer Based on RNA-seq and scRNA-seq.

Colon cancer (CC) is a malignant tumor in the colon. Despite some progress in the early detection and treatment of CC in recent years, some patients still experience recurrence and metastasis. Therefore, it is urgent to better predict the prognosis of CC patients and identify new biomarkers. Recent studies have shown that anoikis-related genes (ARGs) play a significant role in the progression of many tumors. Hence, it is essential to confirm the role of ARGs in the development and treatment of CC by integrating scRNA-seq and transcriptome data. This study integrated transcriptome and single-cell sequencing (scRNA-seq) data from CC samples to evaluate patient stratification, prognosis, and ARG expression in different cell types. Specifically, differential expression of ARGs was identified through consensus clustering to classify CC subtypes. Subsequently, a CC risk model composed of CDKN2A, NOX4, INHBB, CRYAB, TWIST1, CD36, SERPINE1, and MMP3 was constructed using prognosis-related ARGs. Finally, using scRNA-seq data of CC, the expression landscape of prognostic genes in different cell types and the relationship between important immune cells and other cells were explored. Through the above analysis, two CC subtypes were identified, showing significant differences in prognosis and clinical factors. Subsequently, a risk model comprising aforementioned genes successfully categorized all CC samples into two risk groups, which also exhibited significant differences in prognosis, clinical factors, involved pathways, immune landscape, and drug sensitivity. Multiple pathways (cell adhesion molecules (CAMs), and extracellular matrix (ECM) receptor interaction) and immune cells/immune functions (B cell naive, dendritic cell activate, plasma cells, and T cells CD4 memory activated) related to CC were identified. Furthermore, it was found that prognostic genes were highly expressed in various immune cells, and B cells exhibited more and stronger interaction pathways with other cells. The results of this study may provide references for personalized treatment and potential biomarker identification in CC.

In silico genome wide identification of long non-coding RNAs differentially expressed during Candida auris host pathogenesis.

Candida auris is an invasive fungal pathogen of high concern due to acquired drug tolerance against antifungals used in clinics. The prolonged persistence on biotic and abiotic surfaces can result in onset of hospital outbreaks causing serious health threat. An in depth understanding of pathology of C. auris is highly desirable for development of efficient therapeutics. Non-coding RNAs play crucial role in fungal pathology. However, the information about ncRNAs is scanty to be utilized. Herein our aim is to identify long noncoding RNAs with potent role in pathobiology of C. auris. Thereby, we analyzed the transcriptomics data of C. auris infection in blood for identification of potential lncRNAs with regulatory role in determining invasion, survival or drug tolerance under infection conditions. Interestingly, we found 275 lncRNAs, out of which 253 matched with lncRNAs reported in Candidamine, corroborating for our accurate data analysis pipeline. Nevertheless, we obtained 23 novel lncRNAs not reported earlier. Three lncRNAs were found to be under expressed throughout the course of infection, in the transcriptomics data. 16 of potent lncRNAs were found to be coexpressed with coding genes, emphasizing for their functional role. Noteworthy, these ncRNAs are expressed from intergenic regions of the genes associated with transporters, metabolism, cell wall biogenesis. This study recommends for possible association between lncRNA expression and C. auris pathogenesis.

Glutamine Metabolism Heterogeneity in Glioblastoma Unveils an Innovative Combination Therapy Strategy.

Treatment of glioblastoma multiforme (GBM) remains challenging. Unraveling the orchestration of glutamine metabolism may provide a novel viewpoint on GBM therapy. The study presented a full and comprehensive comprehending of the glutamine metabolism atlas and heterogeneity in GBM for facilitating the development of a more effective therapeutic choice. Transcriptome data from large GBM cohorts were integrated in this study. A glutamine metabolism-based classification was established through consensus clustering approach, and a classifier by LASSO analysis was defined for differentiating the classification. Prognosis, signaling pathway activity, tumor microenvironment, and responses to immune checkpoint blockade (ICB) and small molecular drugs were characterized in each cluster. A combinational therapy of glutaminase inhibitor CB839 with dihydroartemisinin (DHA) was proposed, and the influence on glutamine metabolism, apoptosis, reactive oxygen species (ROS), and migration was measured in U251 and U373 cells. We discovered that GBM presented heterogeneous glutamine metabolism-based clusters, with unique survival outcomes, activity of signaling pathways, tumor microenvironment, and responses to ICB and small molecular compounds. In addition, the classifier could accurately differentiate the two clusters. Strikingly, the combinational therapy of CB839 with DHA synergistically attenuated glutamine metabolism, triggered apoptosis and ROS accumulation, and impaired migrative capacity in GBM cells, demonstrating the excellent preclinical efficacy. Altogether, our findings unveil the glutamine metabolism heterogeneity in GBM and propose an innovative combination therapy of CB839 with DHA for this malignant disease.

Identification of potential hub genes and regulatory networks of smoking-related endothelial dysfunction in atherosclerosis using bioinformatics analysis.

Endothelial dysfunction, the earliest stage of atherosclerosis, can be caused by smoking, but its molecular mechanism requires further investigation. This study aimed to use bioinformatics analysis to identify potential mechanisms involved in smoking-related atherosclerotic endothelial dysfunction. The transcriptome data used for this bioinformatics analysis were obtained from the Gene Expression Omnibus (GEO) database. The GSE137578 and GSE141136 datasets were used to identify common differentially expressed genes (co-DEGs) in endothelial cells treated with oxidized low-density lipoprotein (ox-LDL) and tobacco. The co-DEGs were annotated using Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomics (KEGG) databases. Additionally, a protein-protein interaction (PPI) network was constructed to visualize their interactions and screen for hub genes. GSE120521 dataset was used to verify the expression of hub genes in unstable plaques. The miRNA expression profile GSE137580 and online databases (starBase 2.0, TargetScan 8.0 and DGIdb v4.2.0) were used to predict the related non-coding RNAs and drugs. A total of 232 co-DEGs were identified, including 113 up-regulated genes and 119 down-regulated genes. These DEGs were primarily enriched in detrimental autophagy, cell death, transcription factors, and cytokines, and were implicated in ferroptosis, abnormal lipid metabolism, inflammation, and oxidative stress pathways. Ten hub genes were screened from the constructed PPI network, including up-regulated genes such as FOS, HMOX1, SQSTM1, PTGS2, ATF3, DDIT3, and down-regulated genes MCM4, KIF15, UHRF1, and CCL2. Importantly, HMOX1 was further up-regulated in unstable plaques (p= 0.034). Finally, a regulatory network involving lncRNA/circRNA-miRNA-hub genes and drug-hub genes was established. Atherosclerotic endothelial dysfunction is associated with smoking-induced injury. Through bioinformatics analysis, we identified potential mechanisms and provided potential therapeutic targets.

Diversity and host interaction of the gut microbiota in specific pathogen-free pigs

Pigs are widely used as animal models in various studies related to humans. The interaction between the gut microbiota and the host has significant effects on the host’s health and disease status. However, although there have been many studies investigating the pig gut microbiota, the findings have been inconsistent due to variations in rearing conditions. Interactions between the gut microbiota and host have not been fully explored in pigs. Specific pathogen-free (SPF) pigs are ideal non-primate large animals to study the interactions between the gut microbiota and the host. In this study, we performed high-throughput sequencing analysis of the gut microbiota and the gut tissue transcriptome of six SPF pigs to provide a systematic understanding of the composition, function, and spatial distribution of gut microbiota in SPF pigs. We identified significant differences in microbial diversity and functionality among different gastrointestinal tract sites. Metagenomics data analysis revealed significant differences in alpha diversity and beta diversity of microbiota in different gastrointestinal sites of SPF pigs. Additionally, transcriptomic data indicated significant differences in gene expression as well as KEGG and GO functional enrichment between the small intestine and large intestine. Furthermore, by combining microbial metagenomics and host transcriptomics analyses, specific correlations were found between gut microbiota and host genes. These included a negative correlation between the TCN1 gene and Prevotella dentalis, possibly related to bacterial metabolic pathways involving vitamin B12, and a positive correlation between the BDH1 gene and Roseburia hominis, possibly because both are involved in fatty acid metabolism. These findings lay the groundwork for further exploration of the co-evolution between the microbiota and the host, specifically in relation to nutrition, metabolism, and immunity. In conclusion, we have elucidated the diversity of the gut microbiota in SPF pigs and conducted a detailed investigation into the interactions between the gut microbiota and host gene expression. These results contribute to our understanding of the intricate dynamics between the gut microbiota and the host, offering important references for advancements in life science research, bioproduct production, and sustainable development in animal husbandry.

SARS-CoV-2-Induced Type I Interferon Signaling Dysregulation in Olfactory Networks Implications for Alzheimer’s Disease

Type I interferon signaling (IFN-I) perturbations are major drivers of COVID-19. Dysregulated IFN-I in the brain, however, has been linked to both reduced cognitive resilience and neurodegenerative diseases such as Alzheimer’s. Previous works from our group have proposed a model where peripheral induction of IFN-I may be relayed to the CNS, even in the absence of fulminant infection. The aim of our study was to identify significantly enriched IFN-I signatures and genes along the transolfactory route, utilizing published datasets of the nasal mucosa and olfactory bulb amygdala transcriptomes of COVID-19 patients. We furthermore sought to identify these IFN-I signature gene networks associated with Alzheimer’s disease pathology and risk. Gene expression data involving the nasal epithelium, olfactory bulb, and amygdala of COVID-19 patients and transcriptomic data from Alzheimer’s disease patients were scrutinized for enriched Type I interferon pathways. Gene set enrichment analyses and gene–Venn approaches were used to determine genes in IFN-I enriched signatures. The Agora web resource was used to identify genes in IFN-I signatures associated with Alzheimer’s disease risk based on its aggregated multi-omic data. For all analyses, false discovery rates (FDR) <0.05 were considered statistically significant. Pathways associated with type I interferon signaling were found in all samples tested. Each type I interferon signature was enriched by IFITM and OAS family genes. A 14-gene signature was associated with COVID-19 CNS and the response to Alzheimer’s disease pathology, whereas nine genes were associated with increased risk for Alzheimer’s disease based on Agora. Our study provides further support to a type I interferon signaling dysregulation along the extended olfactory network as reconstructed herein, ranging from the nasal epithelium and extending to the amygdala. We furthermore identify the 14 genes implicated in this dysregulated pathway with Alzheimer’s disease pathology, among which HLA-C, HLA-B, HLA-A, PSMB8, IFITM3, HLA-E, IFITM1, OAS2, and MX1 as genes with associated conferring increased risk for the latter. Further research into its druggability by IFNb therapeutics may be warranted.

Nuclear genomes of dinoflagellates reveal evolutionarily conserved pattern of RNA editing relative to stress response

Dinoflagellates are a group of diverse protists with complex genomes whose gene expression regulation mechanisms remain little known. RNA editing is a post-transcriptional regulatory mechanism of gene expression utilized by diverse species, and has been described primarily in the plastid and mitochondrial genomes of dinoflagellates. Its role in post-transcriptional regulation in the nuclear genomes of dinoflagellates remains largely unexplored. Here, integrating genome and transcriptome data from two dinoflagellate taxa in a comparative analysis, we identified 10,486 and 69,953 putative RNA editing sites in the nuclear genomes of the coral symbiont, Durusdinium trenchii CCMP2556 and the free-living bloom-forming taxon, Prorocentrum cordatum CCMP1329. We recovered all 12 possible types of RNA edits, with more edits representing transitions than transversions. In contrast to other eukaryotes, we found a dominance of A-to-T transversion in non-coding regions, many of which were condition-specific. Overall, the RNA editing sites implicate 7.5% of D. trenchii genes and 13.2% of P. cordatum genes. Some sites (1.5% in D. trenchii and more-substantially 62.3% in P. cordatum) were edited at significantly different frequencies in distinct growth conditions. The distribution of editing types and locations exhibited conserved patterns between the two phylogenetically distant species. Interestingly, A-to-T editing within the untranslated regions appear to be associated with upregulation of the edited genes in response to heat stress. These results lend support to the hypothesis that RNA editing is a key molecular mechanism that underpins regulation of gene expression in dinoflagellates.

Aging clocks based on accumulating stochastic variation.

Aging clocks have provided one of the most important recent breakthroughs in the biology of aging, and may provide indicators for the effectiveness of interventions in the aging process and preventive treatments for age-related diseases. The reproducibility of accurate aging clocks has reinvigorated the debate on whether a programmed process underlies aging. Here we show that accumulating stochastic variation in purely simulated data is sufficient to build aging clocks, and that first-generation and second-generation aging clocks are compatible with the accumulation of stochastic variation in DNA methylation or transcriptomic data. We find that accumulating stochastic variation is sufficient to predict chronological and biological age, indicated by significant prediction differences in smoking, calorie restriction, heterochronic parabiosis and partial reprogramming. Although our simulations may not explicitly rule out a programmed aging process, our results suggest that stochastically accumulating changes in any set of data that have a ground state at age zero are sufficient for generating aging clocks.

The Search for and Functional Analysis of Genetic Variants in microRNA-Binding Sites using Massively Parallel Reporter Assay.

A large-scale search for the genetic variants with a bias in the representation of alleles in transcriptome data (AE SNPs) and the binding sites in microRNA 3'-UTRs was performed and their functional significance was assessed using massively parallel reporter assay (MPRA). Of the 629,559 associated "SNP-gene" pairs (eQTLs) discovered in the human liver tissue according to the GTEx Analysis V8 data, 4394 polymorphic positions in the 3'-UTRs of the genes, which represent the eQTLs for these genes were selected. The TargetScanHuman 7.0 algorithm and PolymiRTS database were searched for the potential microRNA-binding sites. Of the predicted microRNA sites affected by eQTL-SNPs, we selected 51 sites with the best evidence of functionality according to Ago2-CLIP-seq, CLEAR-CLIP, and eCLIP-seq for RNA-binding proteins. For MPRA, a library of the plasmids carrying the main and alternative alleles for each AE SNP (in total, 102 constructs) was created. Allele-specific expression for 6 SNPs was detected by transfection of the HepG2 cell line with the constructed plasmid library and sequencing of target DNA and RNA sequences using the Illumina (MiSeq) platform.

Development and characterization of novel EST-SSR markers for Gentiana straminea Maxim., a traditional Tibetan herb in China and cross-amplification in related species

Abstract Gentiana straminea Maxim. (Gentianaceae) is an important traditional Tibetan herb that is mainly distributed on the Qinghai-Tibetan Plateau. Despite its agricultural and pharmacological importance, there remains a paucity of microsatellite markers, particularly expressed sequence tag-simple sequence repeat (EST-SSR) markers, available for this local endemic species. In this study, based on previous Illumina transcriptome data of G. straminea, a total of 96 EST-SSR markers were initially designed and tested. Thirty-two of 96 loci (33.33%) were successfully amplified and verified for validation. Among them, 10 were polymorphic and had clear bands. The polymorphism information content values were 0.09–0.799, the number of alleles per locus ranged from 3 to 14, and the levels of observed and expected heterozygosity were 0.078–0.722 and 0.238–0.884, respectively, which suggested a high level of information. Moreover, cross-amplification was successful for 10 loci in two other related species, Gentiana macrophylla Pallas and Gentiana dahurica Fischer. These EST-SSR markers provide a valuable tool for investigating the genetic diversity related to quantitative traits and population genetic studies on G. straminea and related species in sect. Cruciata Gaudin.

Transcriptome dynamics of Gossypium purpurascens in response to abiotic stresses by Iso-seq and RNA-seq data

Gossypium purpurascens is a member of the Malvaceae family, holds immense economic significance as a fiber crop worldwide. Abiotic stresses harm cotton crops, reduce yields, and cause economic losses. Generating high-quality reference genomes and large-scale transcriptomic datasets across diverse conditions can offer valuable insights into identifying preferred agronomic traits for crop breeding. The present research used leaf tissues to conduct PacBio Iso-seq and RNA-seq analysis. We carried out an in-depth analysis of DEGs using both correlations with cluster analysis and principal component analysis. Additionally, the study also involved the identification of both lncRNAs and CDS. We have prepared RNA-seq libraries from 75 RNA samples to study the effects of drought, salinity, alkali, and saline-alkali stress, as well as control conditions. A total of 454.06 Gigabytes of transcriptome data were effectively validated through the identification of differentially expressed genes and KEGG and GO analysis. Overwhelmingly, gene expression profiles and full-length transcripts from cotton tissues will aid in understanding the genetic mechanism of abiotic stress tolerance in G. purpurascens.

VD3/VDR attenuates bisphenol A-induced impairment in mouse Leydig cells via regulation of autophagy.

Bisphenol A (BPA) is an endocrine disruptor with reproductive toxicity. Further, 1,25-dihydroxyvitamin D3 (VD3) plays an important role in male reproduction by binding vitamin D receptor (VDR) and mediating the pleiotropic biological actions that include spermatogenesis. However, whether VD3/VDR regulates the effect of BPA on Leydig cells (LCs) injury remains unknown. This study aimed to explore the effects of VD on BPA-induced cytotoxicity in mouse LCs. Hereby, LCs treated with BPA, VD3, or both were subjected to the assays of cell apoptosis, proliferation, autophagy, and levels of target proteins. This study unveiled that cell viability was dose-dependently reduced after exposure to BPA. BPA treatment significantly inhibited LC proliferation, induced apoptosis, and also downregulated VDR expression. By jointly analyzing transcriptome data and Comparative Toxicogenomics Database (CTD) data, autophagy signaling pathways related to testicular development and male reproduction were screened out. Therefore, the autophagy phenomenon of cells was further detected. The results showed that BPA treatment could activate cell autophagy, Vdr-/- inhibits cell autophagy, and active VD3 does not have a significant effect on the autophagy of normal LCs. After VD3 and BPA were used in combination, the autophagy of cells was further enhanced, and VD3 could alleviate BPA-induced damage of LCs. In conclusion, this study found that supplementing VD3 could eliminate the inhibition of BPA on VDR expression, further enhance LCs autophagy effect, and alleviate the inhibition of LCs proliferation and induction of apoptosis by BPA, playing a protective role in cells. The research results will provide valuable strategies to alleviate BPA-induced reproductive toxicity.