Functional Genomic Analysis Research Articles

Spatially resolved transcriptomics reveals gene expression characteristics in uveal melanoma

PurposeUveal melanoma (UM) is the most common intraocular malignancy in adults. Previous studies have examined the intra-tumoral heterogeneity. However, the spatial distribution of tumor cells within the tumor microenvironment and its relationship with tumor progression still remains largely unclear. Our study aimed to analyze the correlation between cell distribution patterns and the prognosis of UM.MethodsIn this paper, we performed spatial transcriptomics (ST) sequencing on two UM samples to describe the different cellular distribution patterns. Gene Ontology (GO) and Kyoto Encyclopedia of Genes, Genomes (KEGG) functional enrichment analysis, and protein–protein interaction (PPI) network were performed to define the biological function of each cluster. Differentially expressed genes (DEGs) and survival analysis based on datasets from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) database further confirmed the correlation between cellular distribution and clinical prognosis.ResultsWe found two different patterns of tumor cell distribution. The focal tumor cells have a distinct ribosome synthesis and rRNA pathway. In contrast, the subpopulation tented to distribute diffusely was related to fatty acids metabolism profile, presumably supporting tumor growth by providing energy. The scattered tumor cell cluster was associated with malignant biological behaviors and was involved in extensive cellular interactions, including COLLAGEN. Moreover, pseudo-time analysis showed that migration started from the basal region through cell differentiation. According to the TCGA and GEO database, genes expressed characteristically in the scattered tumor cell cluster were related to poor prognosis.ConclusionsOur study drew the ST maps for UM for the first time. These findings revealed the distribution patterns of tumor cells associated with different biological functions and pointed towards specific tumor subpopulations with higher invasiveness as potential therapeutic targets. Together, our study displayed an overview of UM transcriptome and explored the intra-tumoral heterogeneity of UM at the spatial level.

Genetic parameters and genome-wide association studies including the X chromosome for various reproduction and semen quality traits in Nellore cattle

BackgroundThe profitability of the beef industry is directly influenced by the fertility rate and reproductive performance of both males and females, which can be improved through selective breeding. When performing genomic analyses, genetic markers located on the X chromosome have been commonly ignored despite the X chromosome being one of the largest chromosomes in the cattle genome. Therefore, the primary objectives of this study were to: (1) estimate variance components and genetic parameters for eighteen male and five female fertility and reproductive traits in Nellore cattle including X chromosome markers in the analyses; and (2) perform genome-wide association studies and functional genomic analyses to better understand the genetic background of male and female fertility and reproductive performance traits in Nellore cattle.ResultsThe percentage of the total direct heritability (h2total) explained by the X chromosome markers (h2x) ranged from 3 to 32% (average: 16.4%) and from 9 to 67% (average: 25.61%) for female reproductive performance and male fertility traits, respectively. Among the traits related to breeding soundness evaluation, the overall bull and semen evaluation and semen quality traits accounted for the highest proportion of h2x relative to h2total with an average of 39.5% and 38.75%, respectively. The total number of significant genomic markers per trait ranged from 7 (seminal vesicle width) to 43 (total major defects). The number of significant markers located on the X chromosome ranged from zero to five. A total of 683, 252, 694, 382, 61, and 77 genes overlapped with the genomic regions identified for traits related to female reproductive performance, semen quality, semen morphology, semen defects, overall bulls’ fertility evaluation, and overall semen evaluation traits, respectively. The key candidate genes located on the X chromosome are PRR32, STK26, TMSB4X, TLR7, PRPS2, SMS, SMARCA1, UTP14A, and BCORL1. The main gene ontology terms identified are “Oocyte Meiosis”, “Progesterone Mediated Oocyte Maturation”, “Thermogenesis”, “Sperm Flagellum”, and “Innate Immune Response”.ConclusionsOur findings indicate the key role of genes located on the X chromosome on the phenotypic variability of male and female reproduction and fertility traits in Nellore cattle. Breeding programs aiming to improve these traits should consider adding the information from X chromosome markers in their genomic analyses.

Exploring caffeine as a disruptor of membrane integrity and genomic stability in Staphylococcus aureus: functional and in silico analysis.

To explore the mechanistic underpinnings of caffeine as a potent antibacterial against Staphylococcus aureus ATCC 25923 via in vitro functional assays, whole-genome sequencing, and in silico docking studies. In vitro studies established that caffeine's minimum inhibitory concentration (MIC) against S. aureus ATCC 25923 is 0.01544mmol/mL. Functional assays along with Scanning Electron Microscopy confirmed that caffeine at 0.030089mmol/mL (2MIC) released nucleotide constituents (nucleotide leakage assay) and effluxed potassium ions (potassium efflux assay) thereby, further validating caffeine's role as a membrane-active antimicrobial agent. Whole genome sequencing of control versus caffeine treated samples revealed a significant drop in read mapping percentage from 99.96 to 23.68% and GC content from 30.69 to 6.93%. This massive reduction in the treated sample was a consequence of single nucleotide polymorphisms (SNPs, 50,303), along with insertions and deletions (InDels, 62). Several of these caffeine-induced mutations were found to be harbouring the coding regions of genes involved in processes such as cell membrane organization, bacterial virulence, and DNA repair processes. Thus, implying a caffeine-mediated genomic rearrangement and instability. In silico docking studies revealed a strong binding affinity of caffeine to key cell wall proteins ltaA (-6.9kcal/mol) and ltaS (-6.5kcal/mol) respectively. The dynamic simulation studies revealed caffeine's interaction with receptor ltaS remained stable, with low deviations and minimal fluctuations. Although caffeine has been widely investigated for its antibacterial properties, its specific mechanisms of action, notably its effects on the cell membrane and genomic integrity in S. aureus ATCC 25923, are little understood. This study thus offers a comprehensive functional genomic analysis of caffeine as an antibacterial against S. aureus.

Isolation and characterization of cellulose-mineralizing haloalkaliphilic bacteria from Siberian soda lakes

Soda lakes are unique double-extreme habitats characterized by high salinity and soluble carbonate alkalinity, yet harboring rich prokaryotic life. Despite intensive microbiology studies, little is known about the identity of the soda lake hydrolytic bacteria responsible for the primary degradation of the biomass organic matter, in particular cellulose. In this study, aerobic and anaerobic enrichment cultures with three forms of native insoluble cellulose inoculated with sediments from five soda lakes in south-western Siberia resulted in the isolation of four cellulotrophic haloalkaliphilic bacteria and their four saccharolytic satellites. The final aerobic enrichment included a cellulotrophic bacteroidetes (strain ABcell3) related to Sporocytophaga accompanied by a hemicellulolytic Marinimicrobium strain ABcell2. The anaerobic enrichments resolved in three primary cellulotrophic bacteria and their three saccharolytic bacteroidetes satellites. The culture selected on amorphous cellulose (ac) included a new cellulotrophic member of the Chitinispirillaceae (Fibrobacterota)—strain ANBcel5, and two different saccharolytic satellites from the Marinilabiliales and Balneolales orders. The final enrichment selected on Sigma 101 cellulose consisted of an endospore-forming cellulotrophic strain ANBcel31 belonging to the genus Herbivorax (Acetivibrionales) and its saccharolytic satellite from the Balneolales order. The anaerobic enrichment on a filter paper yielded a binary consortium with the cellulotrophic endospore-forming Halanaerobiales strain ANBcel28 in obligate syntrophy with a cellobiose-utilizing Natronincola. A functional genome analysis of the cellulotrophic isolates confirmed the presence of a large repertoire of genes encoding excreted cellulases, mostly from the GH9 and GH5 families, and indicated that in the endospore-forming anaerobic strains, ANBcel28 and ANBcel31 most of their endo-glucanases are assembled in cellulosomes. Overall, this study showed that cellulose can be mineralized in soda lakes at moderately saline and highly alkaline conditions either by aerobic or fermentative haloalkaliphilic bacteria.

Identification, Characterization, and Ultrastructure Analysis of the Phenol-Degrading Rhodococcus erythropolis 7Ba and Its Viable but Nonculturable Forms.

Phenol and its chlorinated derivatives are introduced into the environment with wastewater effluents from various industries, becoming toxic pollutants. Phenol-degrading bacteria are important objects of research; among them, representatives of the genus Rhodoccocus are often highlighted as promising. Strain 7Ba was isolated by enrichment culture. A new isolate was characterized using culturing, biochemistry, high-throughput sequencing, microscopy (including electron microscopy), and functional genome analysis. Rhodococcus erythropolis strain 7Ba is able to grow on phenol and chlorophenols without losing its properties during long-term storage. It was shown that strain 7Ba is able to form viable but nonculturable (VBNC) forms during long-term storage under nutrient limitation, preserving both cell viability and the ability to degrade phenols. The ultrastructural organization of the vegetative forms of cells and VBNC forms was characterized. The following distinctive features were found: modifications (thickening) of cell membranes, cell size reduction, nucleoid condensation. Functional analysis of the genome showed the presence of genes for the degradation of alkanes, and two branches of the β-ketoadipate pathway for the degradation of aromatic compounds. Also, the genome of strain 7Ba contains several copies of Rpf (resuscitation promoting factor) genes, a resuscitation factor of resting bacterial forms. The new isolate strain 7Ba is a promising biotechnological agent that can not only utilize toxic aromatic compounds but also remain viable during long-term storage. For this reason, its further application as an agent for bioremediation can be successful under changing conditions of climate and given the deficiency of nutrient compounds in nature. Minor biostimulation will allow the strain to recover its metabolic activity and effectively degrade pollution.

Human placental stem cell-based therapies for prevention of abdominal adhesions: A prospective randomized preclinical trial.

Abdominal adhesions are networks of fibrotic tissues that form between organs postoperatively. Current prophylactic strategies do not reproducibly prevent adhesive small bowel obstruction across the entire abdomen. Human placental-derived stem cells produce an anti-inflammatory secretome that has been applied to multiple fibrosing diseases. The purpose of this project is to test human placental stem cell (hPSC)-based therapies for prevention of abdominal adhesions in a clinically relevant rat model. Fifty-four (n = 54, n = 6/group) male Sprague-Dawley rats (250-350 g) underwent model creation and treatment randomization under anesthesia. Experimental groups included human placental-derived stem cells (hPSC, 5 × 106 cells/10 mL Plasmalyte A), human placental-derived stem cells in a hyaluronic acid (HA-Mal-hPSC) hydrogel, the human placental-derived stem cell secretome from conditioned media in 10 mL Plasmalyte A, human placental-derived stem cells' conditioned media in a hyaluronic acid (HA-Mal-CM) hydrogel, Plasmalyte A (media alone, 10 mL), hyaluronic acid hydrogel alone (HA-Mal), Seprafilm (Baxter, Deerfield, IL), and the control groups, model with no treatment (MNT) and sham animals. Treatments were administered intraperitoneally, and the study period was 14 days postoperation. Adhesions were scored at necropsy and analyzed as the difference between means of an index statistic (Animal Index Score) versus MNT. Underlying molecular mechanisms were explored by functional genomic analysis and histology of peritoneal tissues. Hyaluronic acid hydrogel alone, HA-Mal-CM hydrogel, and Seprafilm significantly reduced the overall appearance of abdominal adhesions by mean Animal Index Score at 14 days versus MNT. Human placental stem cell, HA-Mal-hPSC hydrogel, HA-Mal-CM hydrogel, HA-Mal hydrogel alone, and Seprafilm significantly reduced the collagen content of injured peritoneal tissues. Human placental stem cell and HA-Mal-hPSC hydrogel suppressed expression of the most profibrotic genes. Conditioned media, HA-Mal hydrogel alone, and media alone significantly altered the expression of proteins associated with peritoneal fibrotic pathways. Human placental stem cell-based therapies reduce abdominal adhesions in a prospective randomized preclinical trial. This effect is supported by suppression of profibrotic genomic and proteomic pathways.

Study on the anti-colon cancer effect of Renshen Yangrong Decoction based on network pharmacology, molecular docking, and in vivo and in vitro experiments

Renshen Yangrong Decoction (RSYRD) is a traditional Chinese medicinal compound widely used as an adjuvant in colon cancer treatment. However, the specific role of RSYRD in anti-colon cancer remains unclear. This study aimed to investigate the anti-colon cancer effects of RSYRD using network pharmacology, molecular docking, and in vitro and in vivo experiments. The active ingredients and targets of RSYRD were identified utilizing databases, and RSYRD-ingredient-target and protein-protein interaction (PPI) networks were constructed. Gene Ontology (GO) functional analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis were performed on the anti-colon cancer targets of RSYRD, followed by molecular docking using AutoDockTools 1.5.6. In vitro and in vivo experiments were conducted to explore the mechanism of RSYRD and evaluate its inhibitory efficacy on colon cancer. A total of 233 active ingredients and 88 colon cancer-related targets were identified, including 9 core targets. GO functional analysis revealed that RSYRD exerts anti-tumor effects through processes such as oxidative reactions, closely associated with intracellular reactive oxygen species (ROS). KEGG enrichment analysis indicated that the anti-tumor effect of RSYRD involved the interaction of cancer-related and viral disease-associated pathways, as well as the IL-17, TNF, and PI3K-AKT signaling pathways. Molecular docking showed stable ligand-receptor binding, with binding energies below −5.0 kcal/mol. In vitro experiments demonstrated that RSYRD increased ROS levels, reduced mitochondrial membrane potential (MMP), and promoted apoptosis in SW620 cells, thereby inhibiting cell proliferation and migration. In vivo experiments showed that RSYRD significantly suppressed colon cancer proliferation. This study identified Quercetin, Kaempferol, 7-Methoxy-2-methyl-isoflavone, Licochalcone-a, and Isorhamnetin as core ingredients in RSYRD, likely exerting anti-colon cancer effects by inhibiting proteins such as ESR1, HSP90AA1, NCOA2, and MAPK14, and suppressing the PI3K-AKT signaling pathway, thereby regulating tumor cell proliferation, migration, and apoptosis.

Emerging functions of FMNL1 in myeloid neoplasms: insights from bioinformatics to biological and pharmacological landscapes.

Myeloid neoplasms encompass disorders characterized by abnormal myeloid cell proliferation and differentiation, including myelodysplastic syndromes (MDS), myeloproliferative neoplasms, acute myeloid leukemia (AML), and chronic myeloid leukemia (CML). Formin-like protein 1 (FMNL1) is involved in the regulation of the actin cytoskeleton and is predominantly expressed in hematopoietic cells. Given its role in leukemia cell proliferation, survival, migration, and invasion, this study investigates FMNL1 expression in normal hematopoiesis and myeloid neoplasms and explores associations with clinical-laboratory characteristics, mutational status, and survival outcomes in AML. Transcript levels of FMNL1 from several blood-forming cell populations and myeloid neoplasms were extracted from publicly available databases. Myeloid neoplasm cell lines were used for gene/protein expression and cell differentiation studies. Functional genomics analysis was performed using RNA-seq data from The Cancer Genome Atlas (TCGA) AML study, and drug sensitivity predictions were investigated using Beat AML and Genomics of Drug Sensitivity in Cancer (GDSC) datasets. Statistical analyses assessed the impact of FMNL1 expression on clinical outcomes. FMNL1 was highly expressed in metamyelocytes, neutrophils, and monocytes compared to hematopoietic stem cells, and its expression increased with granulocytic differentiation. FMNL1 expression was elevated in AML and CML patients compared to healthy donors. FMNL1 expression was not significantly associated with clinical-laboratory characteristics or survival outcomes but showed a higher frequency of WT1 transcription factor (WT1) mutations with low FMNL1 expression in AML patients. High FMNL1 expression in AML correlated with immune response and inflammatory activity pathways. FMNL1 mRNA levels influenced drug sensitivity in AML models, with correlations observed for specific antineoplastic agents. FMNL1 plays a potential role in granulocyte differentiation and function, and its differential expression is linked to critical signaling pathways in leukemogenesis and inflammation. These findings highlight FMNL1's potential therapeutic implications in myeloid neoplasia, warranting further investigation.

Analysis of Characteristics of Bovine-Derived Non-Enterotoxigenic Bacteroides fragilis and Validation of Potential Probiotic Effects.

Bacteroides fragilis is a new generation of probiotics, and its probiotic effects on humans and some animals have been verified. However, research on B. fragilis in cattle is still lacking. In this study, 24 stool samples were collected from two large-scale cattle farms in Wuzhong, Ningxia, including 12 diarrheal and 12 normal stools. A non-toxigenic Bacteroides fragilis (NTBF) was isolated and identified by 16S rRNA high-throughput sequencing and named BF-1153; genome composition and genome functional analyses were carried out to reflect the biological characteristics of the BF-1153 strain. A cluster analysis of BF-1153 was performed using Mega X to explore its genetic relationship. In addition, Cell Counting Kit-8 (CCK8) was used to determine the toxic effects of the strain on human ileocecal colorectal adenocarcinoma cell line cells (HCT-8), Madin-Darby bovine kidney cells (MDBK), and intestinal porcine epithelial cells (IPECs). The results showed that BF-1153 conformed to the biological characteristics of B. fragilis. BF-1153 had no toxic effects on HCT-8, MDBK, and IPEC. Animal experiments have shown that BF-1153 has no toxic effects on healthy SPF Kunming mice. Notably, the supernatant of BF-1153 enhanced cell activity and promoted cell growth in all three cell lines. At the same time, a cluster analysis of the isolated strains showed that the BF-1153 strain belonged to the same branch as the B. fragilis strain 23212, and B. fragilis strain 22998. The results of the animal experiments showed that BF-1153 had a certain preventive effect on diarrhea symptoms in SPF Kunming mice caused by a bovine rotavirus (BRV). In summary, the strain BF-1153 isolated in this experiment is NTBF, which has no toxic effect on MDBK, HCT-8, and IPEC, and has obvious cell growth-promoting effects, especially on MDBK. BF-1153 promotes the growth and development of SPF Kunming mice when compared with the control group. At the same time, BF-1153 alleviated the diarrhea symptoms caused by BRV in SPF Kunming mice. Therefore, BF-1153 has the potential to be a probiotic for cattle.

Comprehensive analysis of mRNA and miRNA differential expression profiles in the hypothalamus-pituitary-gonadal axis in laying and broodiness period of Wanxi white geese

The hypothalamus-pituitary-gonadal (HPG) axis is an important neuroendocrine regulatory center involved in egg-laying process in poultry. However, its mechanism of regulating broodiness behavior and laying performance in geese remains unclear. This study explored the molecular mechanism by which the HPG axis regulates brooding behavior in Wanxi white geese (WWG). The hypothalamus, pituitary, and ovarian tissues of Wanxi white geese were collected at laying and brooding periods for transcriptome sequencing analysis. A total of 240 (BH vs. LH), 319 (BP vs. LP), and 445 (BO vs. LO) differentially expressed genes, and 56 (BH vs. LH), 82 (BP vs. LP), and 48 (BO vs. LO) differentially expressed miRNAs were identified. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment analysis showed that differentially expressed genes (DEGs) and differentially expressed miRNAs (DEMs) were significantly enriched in hormone level regulation, cell communication, calcium signaling pathway, GnRH signaling pathway, MAPK signaling pathway, Wnt signaling pathway, and other processes. Six DEGs and four DEMs were randomly selected for real-time fluorescence quantitative reverse transcription PCR (RT-qPCR). The results showed that the transcriptome sequencing data were accurate and reliable. In addition, 22 potential hub miRNAs were screened. Dual luciferase reporter assays confirmed the targeting relationship between miR-144-y and DIO3. The results showed that the miRNAs mainly regulated the laying performance and brooding behavior of WWG by mediating the expression of target genes. In this study, we systematically elucidated the mechanisms by which the HPG axis regulates the broodiness behavior and laying performance of WWG at the post-transcriptional level. Several miRNAs and mRNAs associated with the reproductive performance of WWG were identified, providing a crucial reference for the subsequent use of gene editing technologies to breed new varieties and advance the development of WWG breeding industry.

Integrative Analyses of Genes of Pediatric Non-alcoholic Fatty Liver Disease Associated with Energy Metabolism.

Pediatric non-alcoholic fatty liver disease (NAFLD) is a chronic steatosis of the liver associated with energy metabolism in children and adolescents, failure to intervene promptly can elevate the risk of developing hepatocellular carcinoma. Therefore, this study aimed to understand the underlying mechanism of pediatric NAFLD and investigate potential biomarkers and therapeutic targets. We investigated genes using the GSE185051 data set related to energy metabolism from the GeneCards database, constructed protein-protein interaction network, identified hub genes and established networks representing interactions between these hub genes and miRNA, RNA-binding proteins, transcription factors, and drugs. Subsequently, we performed Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment analysis, Gene Set Enrichment Analysis (GSEA), and immune infiltration analysis. Our analysis identified 9 hub genes through the PPI network. The target molecules were identified through the interaction network between hub genes and miRNAs, RNA-binding proteins, transcription factors, and drugs. GO analysis revealed that hub genes were associated with oxidative stress responses and other pathways. KEGG analysis highlighted their involvement in pathways such as insulin resistance, among others. GSEA revealed that hub genes were highly enriched in pathways related to Omega-9 fatty acid synthesis, among others. Immune infiltration analysis suggested that mast cells and T follicular helper cells play significant roles in the pathogenesis of NAFLD. We identified the hub genes in pediatric NAFLD closely related to energy metabolism. These findings offer the potential for identifying potential novel diagnostic biomarkers, and establishing therapeutic targets for pediatric NAFLD.

Genome sequencing and functional analysis of potential Trichoderma species for controlling Pythium schmitthenneri-Induced root rot in olive trees

Root rot disease caused by Pythium schmitthenneri is a significant challenge in olive (Olea europaea L.) cultivation. Due to the limitations of chemical control, this study explores the potential of Trichoderma species as a sustainable alternative with an overview in the genetic properties involved in it. Eight Trichoderma isolates were obtained from olive roots and identified using rDNA internal transcribed spacer (ITS) sequences, revealing four species: T. harzianum, T. longibrachiatum, T. virens, and T. viride. Among these, T. harzianum T-BM6 and T. virens T-KH4, showed 81.4 and 78.77 % inhibition rates of pathogen growth, respectively, in vitro. Greenhouse trials demonstrated that these isolates reduced disease severity to 18.75 % for T-BM6 and 31.25 % for T-KH4. All isolates produced cellulase, with T-BM6 exhibiting high cellulolytic activity. Genomic analysis highlighted key genes related to the MAPK signaling pathway, sulfur metabolism, and inositol phosphate metabolism in T-BM6, and phenylalanine, tyrosine, and tryptophan metabolism in T-KH4. These results suggest that T. harzianum T-BM6 and T. virens T-KH4 are promising candidates for biological control, offering a sustainable alternative to chemical treatments for managing olive root rot.

The Functional Transcriptomic Landscape Informs Therapeutic Strategies in Multiple Myeloma.

Several therapeutic agents have been approved for treating multiple myeloma (MM), a cancer of bone marrow resident plasma cells. Predictive biomarkers for drug response could help guide clinical strategies to optimize outcomes. Here, we present an integrated functional genomic analysis of tumor samples from MM patients that were assessed for their ex vivo drug sensitivity to 37 drugs, clinical variables, cytogenetics, mutational profiles, and transcriptomes. This analysis revealed a MM transcriptomic topology that generates "footprints" in association with ex vivo drug sensitivity that have both predictive and mechanistic applications. Validation of the transcriptomic footprints for the anti-CD38 monoclonal antibody daratumumab and the nuclear export inhibitor selinexor demonstrated that these footprints can accurately classify clinical responses. The analysis further revealed that daratumumab and selinexor have anti-correlated mechanisms of resistance, and treatment with a selinexor-based regimen immediately after a daratumumab-containing regimen was associated with improved survival in three independent clinical trials, supporting an evolutionary-based strategy involving sequential therapy. These findings suggest that this unique repository and computational framework can be leveraged to inform underlying biology and to identify therapeutic strategies to improve treatment of MM.

First De Novo genome assembly and characterization of Gaultheria prostrata.

Gaultheria Kalm ex L. (Ericaceae), a type of evergreen shrub, known as a natural source of methyl salicylate, possesses rich germplasm resources, strong habitat adaptability, significant ornamental value, and noteworthy pharmacological activities. However, due to the paucity of whole genomic information, genetically deep research in these areas remains limited. Consequently, we intend to obtain genome data through high-throughput sequencing, gene annotation, flow cytometry, transcription factors prediction and genetic marker analysis for a representative species of this genus, with Gaultheria prostrata selected for our study. In this study, we preliminarily obtained the genome of G. prostrata through next-generation sequencing methods. Utilizing 47.94 Gb of high-quality sequence data (108.95× coverage), assembled into 114,436 scaffolds, with an N50 length of 33,667 bp. The genome size assembled by SOAPdenovo, approximately 417 Mb, corresponded closely to predictions by flow cytometry (440 Mb) and k-mer analysis (447 Mb). The genome integrity was evaluated using BUSCO with 91%. The heterozygosity ratio was 0.159%, the GC content was 38.85%, and the repetitive regions encompassed over 34.6% of the genome. A total of 26,497 protein-coding genes have been predicted and annotated across Nr, Swissprot, GO, KEGG, and Pfam databases. Among these, 14,377 and 2,387 genes received functional annotation in Nr and Swissprot, respectively; 21,895, 24,424, and 22,330 genes were similarly annotated in GO, KEGG, and Pfam. Moreover, A total of 279,785 SSRs were identified and 345,270 primers for these SSRs were designed. Within the various nucleotide types of SSRs, AG/CT and AAG/CTT constituted the predominant dinucleotide and trinucleotide repeat types in G. prostrata. In addition, 1,395 transcription factors (TFs) from 75 TF families, 462 transcription regulators (TRs) from 33 TR families and 840 protein kinase (PKs) from 118 PK families were identified in this genome. We also performed phylogenetic analyses of G. prostrata and related species, including estimation of divergence times and expansion and contraction analyses, followed by positive selection analyses of orthologous gene pairs of G. prostrata and its close relative Vaccinium corymbosum. These results provide a reference for in-depth study of genus Gaultheria, contributing to future functional and comparative genomics analyses and providing supporting data for the development of molecular markers.

Genome-wide association study uncovers pea candidate genes and metabolic pathways involved in rust resistance.

Pea (Pisum sativum L.) is an important temperate legume crop providing plant-based proteins for food and feed worldwide. Pea yield can be limited by several biotic stresses, among which rust represents a major limiting factor in many temperate and subtropical regions. Some efforts have been made to assess the natural variation in pea resistance to rust, but its efficient exploitation in breeding is limited since the resistance loci identified so far are scarce and their responsible gene(s) unknown. To overcome this knowledge gap, a comprehensive genome-wide association study (GWAS) has been performed on pea rust, caused by Uromyces pisi, to uncover genetic loci associated with resistance. Utilizing a diverse collection of 320 pea accessions, we evaluated phenotypic responses to two rust isolates using both traditional methods and advanced image-based phenotyping. We detected 95 significant trait-marker associations using a set of 26,045 Diversity Arrays Technology-sequencing polymorphic markers. Our in silico analysis identified 62 candidate genes putatively involved in rust resistance, grouped into different functional categories such as gene expression regulation, vesicle trafficking, cell wall biosynthesis, and hormonal signaling. This research highlights the potential of GWAS to identify molecular markers associated with resistance and candidate genes against pea rust, offering new targets for precision breeding. By integrating our findings into current breeding programs, we can facilitate the development of pea varieties with improved resistance to rust, contributing to sustainable agricultural practices and food security. This study sets the stage for future functional genomic analyses and the application of genomic selection approaches to enhance disease resistance in peas.

Analysis of the composition and function of rhizosphere microbial communities in plants with tobacco bacterial wilt disease and healthy plants.

To explore the factors influencing the occurrence of bacterial wilt, the differences in the physicochemical properties, microbial community composition and function between rhizosphere soil of tobacco plants with bacterial wilt and healthy plants in the tobacco planting area of Fuzhou City, Jiangxi Province were analyzed and compared. The results showed that the rhizosphere soil of diseased tobacco exhibited significantly reduced levels of exchangeable potassium, water-soluble potassium, nitrate nitrogen, total nitrogen and pH, in comparison to the rhizosphere soil of healthy plants. Conversely, the available phosphorus content of the rhizosphere soil of diseased tobacco was significantly increased. The amount of Ralstonia solanacearum in soil was negatively correlated with pH, nitrate nitrogen and total nitrogen, and positively correlated with exchangeable potassium and water-soluble potassium. A total of 43 genera were significantly different between the two groups of rhizosphere soil, of which 24 genera were enriched in the rhizosphere of healthy plants, including Ideonella, Rhizophagus, Rhizobacter, Altererythrobacter and Ignavibacterium associated with plant disease resistance, Thermodesulfovibrio, Syntrophorhabdus, Syntrophus, Chlorobium, Hydrogenophaga and Limnohabitans associated with soil sulfur metabolism, as well as Ignavibacterium, Ideonella, Derxia and Azohydromonas associated with soil nitrogen cycling. Kyoto Encyclopedia of Genes and Genomes functional analysis of the unigenes obtained by metagenomic sequencing also showed that the differential unigenes were significantly enriched in the sulfur metabolism pathway. In addition, the rhizosphere soil of diseased tobacco plants exhibited a higher abundance of antibiotic-producing actinomycetes and an increased load of antibiotic resistance genes compared to that of healthy plants. In general, lower pH value, less content of nitrate nitrogen and total nitrogen, and more content of exchangeable potassium and water-soluble potassium could contribute to onset of bacterial wilt. Twenty-four genera, including Ideonella and Rhizophagus, may construct a healthy microecological network in the rhizosphere of tobacco plants. All these factors may interact with each other to control the development of bacterial wilt. This complicated interaction network needs to be explored further.IMPORTANCEPrevious studies have mainly focused on the differences in microbial species composition between healthy and diseased soils, but the differences in microbial community functions between two types of soil have not been well characterized. In this study, soil samples in diseased and healthy plant rhizospheres were collected for physicochemical property testing and metagenomic sequencing. We focused on analyzing the differences in physicochemical properties and microbial community functions between these soils, as well as the correlation between these factors and pathogen content. The results of this study provide a theoretical basis for further understanding the occurrence of tobacco bacterial wilt in the field.

Integrative analysis of bulk and single-cell transcriptomic data reveals novel insights into lipid metabolism and prognostic factors in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is associated with high mortality rate. This study investigated the status of lipid metabolism-related genes in HCC. Bulk transcriptomic and single-cell sequencing data for HCC were retrieved from public databases. The single-cell sequencing data was subjected to dimensionality reduction, which facilitated the annotation of distinct cell subpopulations and marker gene expression analysis within each subpopulation. Genes associated with lipid metabolism in liver cells were identified, and a machine-learning model was developed using the bulk transcriptomic data randomly partitioned into training and validation sets. The efficacy of the model was validated using these two sets. A multifactorial Cox analysis on the model genes combined with clinical features, led to the identification of age, HMGCS2, HNRNPU, and RAN as independent prognostic factors, which were included in the nomogram model construction and validation. A weighted gene co-expression analysis of all genes of the bulk transcriptome samples revealed the correlation between gene modules and risk score. Genes with cor > 0.4 in the highest-expressing module were selected for Gene Ontology and Kyoto Encyclopedia of Genes and Genomes functional enrichment analysis. Immune-related analysis was conducted based on seven algorithms for immune cell infiltration prediction. For the genes in the nomogram model, the expression in clinical pathological factors was also analyzed. The drug sensitivity analysis offered a reference for the selection of targeting drugs. This investigation provides novel insights and a theoretical basis for the prognosis, treatment, and pharmaceutical advancements for patients diagnosed with HCC.

PneumoBrowse 2: an integrated visual platform for curated genome annotation and multiomics data analysis of Streptococcus pneumoniae.

Streptococcus pneumoniae is an opportunistic human pathogen responsible for high morbidity and mortality rates. Extensive genome sequencing revealed its large pangenome, serotype diversity, and provided insight into genome dynamics. However, functional genome analysis has lagged behind, as that requires detailed and time-consuming manual curation of genome annotations and integration of genomic and phenotypic data. To remedy this, PneumoBrowse was presented in 2018, a user-friendly interactive online platform, which provided the detailed annotation of the S. pneumoniae D39V genome, alongside transcriptomic data. Since 2018, many new studies on S. pneumoniae genome biology and protein functioning have been performed. Here, we present PneumoBrowse 2 (https://veeninglab.com/pneumobrowse), fully rebuilt in JBrowse 2. We updated annotations for transcribed and transcriptional regulatory features in the D39V genome. We added genome-wide data tracks for high-resolution chromosome conformation capture (Hi-C) data, chromatin immunoprecipitation coupled to high-throughput sequencing (ChIP-Seq), ribosome profiling, CRISPRi-seq gene essentiality data and more. Additionally, we included 18 phylogenetically diverse S. pneumoniae genomes and their annotations. By providing easy access to diverse high-quality genome annotations and links to other databases (including UniProt and AlphaFold), PneumoBrowse 2 will further accelerate research and development into preventive and treatment strategies, through increased understanding of the pneumococcal genome.

To Reveal the Potential Mechanism of Quercetin against NSCLC Based on Network Pharmacology and Experimental Validation.

This study aimed to initially clarify the potential mechanism of quercetin in the treatment of non-small cell lung cancer (NSCLC) based on network pharmacology, molecular docking and in vitro experiments. TCMSP, SwissTargetPrediction, TCMIP, STITCH, and ETCM databases were applied to obtain the targets of quercetin. NSCLC-related targets were retrieved from GeneCards, OMIM, PharmGKB, TTD, and NCBI databases. Their intersection targets were imported into the STRING database to construct a protein-protein interaction (PPI) network and core targets were identified through the Cytoscape 3.10.0 soft and the CytoHubba tool. Furthermore, Gene Ontology (GO) functional analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were performed on the intersection targets. A compound-targetspathways network was subsequently constructed to screen for key targets and pathways. Molecular docking was performed with Discovery Studio software to verify the interactions between quercetin and core targets. In vitro validations were conducted employing CCK-8 assays, flow cytometry, and Western blotting (WB). 193 potential targets of quercetin for treating NSCLC were obtained. The top ten core targets identified within the PPI network included TP53, HSP90AA1, AKT1, JUN, SRC, EGFR, ACTB, TNF, MAPK1, and VEGFA. GO analysis yielded 2319 items, and KEGG analysis resulted in 211 enriched pathways. Molecular docking results demonstrated a high affinity of quercetin towards the core targets. Based on the compound-targets-pathways network and molecular docking, the PI3K/AKT/P53 pathway and its key-related proteins (PIK3R1, AKT1, and TP53) were selected for further validation. Quercetin(20 and 40 μg/mL) significantly decreased the viability of A549 NSCLC cells but not BEAS-2B normal cells via CCK-8 assays. Flow cytometry and WB analyses further corroborated that quercetin could promote apoptosis of A549 cells by downregulating and upregulating the expression of Bcl-2 and Bax (P<0.05), respectively. Notably, quercetin did not significantly alter the total protein levels of PI3K, AKT, and P53 but downregulated the phosphorylation levels of PI3K and AKT (P<0.05) and upregulated the phosphorylation level of P53 (P<0.05). Quercetin exhibits therapeutic potential in NSCLC by regulating the PI3K/AKT/P53 pathway to promote cell apoptosis.

Functional genomic analysis of the 68-1 RhCMV-Mycobacteria tuberculosis vaccine reveals an IL-15 response signature that is conserved with vector attenuation.

Tuberculosis (TB), caused by Mycobacterium tuberculosis (Mtb) is a deadly infectious disease having a major impact on global health. Using the CMV vector for development of novel vaccines is a promising new strategy that elicits strong and durable, high frequency memory T cell responses against heterologous immunogens. We conducted functional transcriptomic analysis of whole blood samples collected from cohorts of rhesus (Rh) macaques that were administered RhCMV/TB vector using a prime-boost strategy. Two modified CMV vectors were used in this study, including 68-1 RhCMV/TB-6Ag (encoding 6 Mtb protein immunogens, including Ag85A, ESAT-6, Rv3407, Rv2626, Rpf A, and Rpf D) and its attenuated variant, 68-1 RhCMV/Δpp71-TB-6Ag (a cell-to-cell spread-deficient vaccine vector lacking the Rh110 gene encoding the pp71 tegument protein). Bulk mRNA sequencing, differential gene expression, and functional enrichment analyses showed that these RhCMV/TB vaccines induce the innate and adaptive immune responses with specific transcriptomic signatures, including the IL-15-induced protective gene signature previously defined to be linked with protection against simian immunodeficiency virus (SIV) by the 68-1 RhCMV/SIV vaccine. While both vectors exhibited a transcriptomic response of the IL-15 protective signature in whole blood, we show that lack of pp71 does not maintain induction of the protective signature for the full duration of the study compared to the parental non-attenuated vector. Our observations indicate that RhCMV vector vaccines induce a transcriptomic response in whole blood that include a conserved IL-15 signature of which vector-encoded pp71 is an important component of response durability that upon future Mtb challenge may define specific vaccine protection outcomes against Mtb infection.