Search PubMedSearch

SEARCH · Search PubMed

Results for “functional annotations”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7Linked to original sources

Whole-Genome Analysis of Bacillus Licheniformis Ali5 and Synthesis of Lichenysin via Genome Shuffling.

Whole-genome sequencing of Bacillus licheniformis Ali5 was performed via MGI-seq PE150 and Nanopore single-molecule real-time sequencing. The strain has a 4,114,664 bp circular genome encoding 4030 protein-coding genes. Functional annotation across NR, COG, GO, KEGG, CARD, BacMet, and CAZy databases identified 4025, 2812, 988, 1242, 72, 69, and 94 corresponding genes, respectively, and antiSMASH 6.0 revealed multiple antimicrobial biosynthetic gene clusters, including intact lichenysin and lichenicidin VK21 A1/A2 gene clusters. Three rounds of recursive protoplast fusion-based genome shuffling, paired with a dual-index screening system, significantly improved strain growth and lichenysin biosynthesis. Recombinants exhibited shortened lag phase, enhanced proliferation, improved stationary-phase stability, and higher diauxic peak biomass. PP3-176 and PP3-186 showed 4.6%-8.1% higher 12-h shake-flask titer and 3.1%-4.0% higher maximum titer than the parental average, with excellent fermentation stability. 1-L bioreactor validation confirmed strong scale-up potential. PP3-186 achieved 27.2% and 31.6% titer increases at 12 h and 20 h, while PP3-176 yielded 20.4% and 14.6% improvements with robust metabolic performance. This study validates genome shuffling as an effective strategy for enhancing lichenysin production, providing candidate strains and technical support for industrial application.

Bacillus licheniformis

Proteomics-Based Identification of the Pyroptosis-Related Biomarker PCSK9 and Its Association With the Pathogenesis of Rheumatoid Arthritis.

Rheumatoid arthritis (RA) is a common autoimmune disease, and early diagnosis is critical for effective treatment. This study aims to identify potential biomarkers related to pyroptosis through serum proteomics analysis, offering new insights for the early diagnosis of RA. We enrolled 100 participants, including 50 patients with RA and 50 healthy controls. Serum samples were collected and analyzed using high-resolution liquid chromatography-tandem mass spectrometry (LC-MS/MS) for proteomics profiling. Differential protein expression analysis and functional annotation revealed significant upregulation of pyroptosis-related proteins in the serum of patients with RA. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses, along with protein-protein interaction (PPI) network analysis, showed that these proteins are involved in inflammation and immune pathways, particularly the activation of the NOD-like receptor protein 3 (NLRP3) inflammasome. Enzyme-linked immunosorbent assay (ELISA) validation confirmed a significant increase in PCSK9 levels in patients with RA, suggesting that PCSK9 may play a key role in the pathogenesis of RA. This study provides new directions for biomarker research in RA, particularly regarding the potential involvement of the pyroptosis pathway, with significant clinical application prospects.

Humans

Network based approach identifies miR-145-3p as a central regulatory hub associated to the progression from localized to metastatic medullary thyroid carcinoma.

Medullary thyroid carcinoma (MTC) is a neuroendocrine tumor originating from calcitonin producing C-cells and accounts for 1-5% of thyroid cancers. Total thyroidectomy is curative in localized disease (N0), whereas lymph node metastases (N1) are associated with poorer prognosis. However, the molecular mechanisms driving the metastatic shift remain poorly understood. This study aimed to identify miRNA features linked to metastatic spread in MTC, focusing on the transition from N0 to N1. Co-expression networks were constructed for N0 and N1 tumors, and differential connectivity analysis was used to identify key miRNAs acting as regulatory hubs. Functional annotation of their target genes was performed using the Kyoto Encyclopedia of Genes and Genomes (KEGG), Gene Ontology (GO), and Reactome pathway analyses. Validation experiments were carried out in MTC cells to evaluate the effects of selected miRNAs on cell proliferation, survival, and MAPK pathway activation. Network analysis revealed distinct miRNA co-expression patterns between N0 and N1 tumors. Differential network analysis highlighted miR-145-3p as a central regulatory hub, exhibiting 29 altered co-expression changes and a marked loss of connectivity in N1. Target enrichment identified 59 validated genes, including key oncogenic drivers such as MYC, PTEN, BCL2, PIK3CA, AKT1, and MAPK7. In MTC cells, simultaneous inhibition of miR-145-3p together with its top co-expressed miRNAs increased proliferation and survival, and enhanced ERK phosphorylation, indicating MAPK pathway activation and a shift toward a more aggressive phenotype. In conclusion, this study identifies a miRNA regulatory hub centered on miR-145-3p that is associated with metastatic progression and highlights the value of network-based approaches in uncovering mechanisms of cancer dissemination. © 2026 The Author(s). The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

MAPK signaling

Genomic and Structural Analysis of Gamete Recognition Proteins in a Broadcast Spawning Echinoderm Mesocentrotus franciscanus.

Gamete recognition proteins are expressed on the surfaces of sperm and eggs, where they mediate interactions between gametes. The genetic basis for gamete recognition proteins, as well as their structure and interactions, have yet to be fully resolved. Using a new high-quality de novo genome assembly for the sea urchin Mesocentrotus franciscanus, we investigated the genomic structure, expression, and protein forms of several gamete recognition proteins: sperm bindin, egg receptor for sperm (HSP110), and egg bindin receptor (EBR1), as well as the receptor for egg jelly (REJ) and its paralogs. To inform future population genetic and evolutionary studies, we resolve the genomic structure of the large EBR1 protein, identifying fewer tandem CUB-TSP1 repeats in EBR1 compared to the initial characterization of this protein. As expected for an egg receptor for sperm, EBR1 is highly expressed in female reproductive tissues (eggs and female gonad), compared to other tissues. In contrast, HSP110 shows similar levels of expression across male and female reproductive tissues, as well as across non-reproductive tissues and development stages. HSP110 might be a pleiotropic gene that in part influences fertilization. Using protein structural modeling and functional domain predictions, we propose hypotheses about potential interactions among EBR1, bindin, and HSP110 proteins that may provide insight into sperm-egg interactions in sea urchins. Resolving the genomic structure of genes encoding gamete recognition proteins, in combination with functional annotations and protein structural modeling, enables deeper investigation into the consequences of variation in gamete recognition proteins and the evolution of reproductive isolation.

Mesocentrotus franciscanus

Hybrid genome assembly and phenotypic assays reveal carbohydrate metabolism diversity in Lacticaseibacillus strains.

Investigation of carbohydrate metabolism in lactic acid bacteria is essential for the rational selection of strains for fermentation processes, particularly in emerging applications involving non-conventional substrates or building of synthetic microbial consortia. However, establishing robust genotype-phenotype relationships remains challenging, as gene presence alone often fails to explain observed metabolic traits without considering the genomic context and regulatory architecture. In the present study, we combined hybrid genome assembly (Illumina and Oxford Nanopore) with high-throughput phenotype profiling (Biolog GENIII and PM2A) to investigate carbohydrate utilization in five Lacticaseibacillus strains. Phenotypic assays revealed clear intra- and inter-specific variability in substrate utilization. We therefore investigated whether such differences could be attributed to the organization and regulatory context of carbohydrate-associated loci, rather than to gene presence alone. Functional annotation based on COG and CAZyme databases revealed candidate genomic regions potentially involved in carbohydrate metabolism. Comparative analysis between predicted and experimentally observed substrate usage highlighted specific loci associated with carbohydrate utilization profile. The trehalose (tre) operon was conserved across all strains, while at least two distinct cellobiose-associated loci were detected in each genome. Despite the presence of these loci, L. paracasei strains were unable to metabolize cellobiose, a phenotype likely linked to the presence of a downstream TetR-type transcriptional repressor within the cellobiose (cel) operon. Additionally, a genomic region uniquely found in L. rhamnosus strains was associated with gentiobiose utilization, consistent with phenotypic observations. Overall, these findings highlight the importance of integrating phenotypic validation with complete genome context to support the identification of candidate structural and regulatory determinants of carbohydrate utilization in lactic acid bacteria. KEY POINTS: • Phenotype microarrays reveal metabolic traits of interest in isolated strains. • Regulatory context is key to understanding carbohydrate metabolism differences. • Basis of subspecies-dependent cellobiose metabolism in L. paracasei is provided.

Carbohydrate Metabolism

A genetic signal at 8q12.3 modulates GGT levels via the Runx1-CYP7B1 axis in female ethnic minorities from Guizhou.

Gamma-glutamyl transferase (GGT) regarded as a biomarker of liver dysfunction or excessive alcohol consumption; however, existing genome-wide association studies (GWAS) have been conducted predominantly in European populations and East Asian populations from Japan and the Taiwan region, with limited investigation in ethnic minorities from Guizhou Province. Previous genetic studies have demonstrated that Guizhou ethnic minorities share an East Asian genetic background while exhibiting specific genetic structures, a pattern that is also confirmed by our principal component analysis (PCA) results. We therefore performed a GWAS in this population and identified a genome-wide significant signal at 8q12.3 in female ethnic minorities from Guizhou. Fine-mapping and functional annotation analyses suggest that a regulatory pathway involving Runt-related transcription factor 1 (Runx1)-Cytochrome P450 family 7 subfamily B member 1 (CYP7B1)-cholesterol-reactive oxygen species (ROS)-glutathione (GSH) may contribute to the regulation of GGT levels. Mendelian randomization (MR) analyses further supported a causal relationship between GGT levels and autoimmune hepatitis (AIH). These findings uncover a genetic mechanism underlying GGT variation at 8q12.3 in female ethnic minorities from Guizhou, implicating a pathway linked to cholesterol metabolism and oxidative stress, and providing potential targets and insights for precision prevention and treatment of related diseases.

Female

Genomic and polyphasic characterization of six novel Hymenobacter species isolated from soil in Korea.

Six bacterial strains (BT523T, BT559T, 15J16-1T3BT, BT730T, DG25AT, and DG25BT) were isolated from soil samples in Korea and assigned to the family Hymenobacteraceae (order Cytophagales, class Cytophagia). Phylogenetic analysis based on 16S rRNA gene sequences showed that the strains formed distinct lineages within the genus Hymenobacter. Strains BT523T, BT559T, and 15J16-1T3BT exhibited highest sequence similarities to Hymenobacter armeniacus BT189T (97.6%), Hymenobacter polaris RP-2-7 T (97.8%), and Hymenobacter paludis KBP-30 T (98.3%), respectively, while strains BT730T, DG25AT, and DG25BT were most closely related to Hymenobacter tibetensis XTM003T, with similarities of 96.3-96.6%. All strains were Gram-negative, aerobic, rod-shaped, and formed red to pink pigmented colonies. Whole-genome analysis revealed genome sizes ranging from 3.78 to 6.33 Mb with G + C contents of 55.5-65.0%. Average nucleotide identity (ANI) and digital DNA-DNA hybridization (dDDH) values between the strains and their closest relatives were below the accepted thresholds for species delineation, supporting their classification as novel species. Functional annotation indicated the presence of genes associated with core metabolism, stress response, and pigment biosynthesis, reflecting adaptation to soil environments. Secondary metabolite analysis further revealed the presence of biosynthetic gene clusters, including terpene and siderophore pathways. Based on polyphasic taxonomic evidence, the six strains are proposed to represent six novel species of the genus Hymenobacter, for which the names Hymenobacter miniatus sp. nov., Hymenobacter madidus sp. nov., Hymenobacter convexus sp. nov., Hymenobacter rubellus sp. nov., Hymenobacter erythromyxa sp. nov., and Hymenobacter radioresistens sp. nov. are proposed. The type strains are BT523T (= KCTC 72341 T = NBRC 114851 T), BT559T (= KACC 21821 T = NBRC 114852 T), 15J16-1T3BT (= KCTC 42995 T = NBRC 112818 T), BT730T (= KACC 22457 T = NBRC 116482 T), DG25AT (= KCTC 32451 T = TBRC 19856 T), and DG25BT (= KCTC 32452 T = JCM 19444 T).

Soil Microbiology

Whole-genome analysis of Brevibacterium sanguinis AZMABM HM27: a bacterial isolate from the sea anemone Radianthus magnifica and exhibiting promising multi-therapeutic properties.

BACKGROUND: The marine anemone Radianthus magnifica harbors symbiotic microbes with promising biomedical potential, yet their diversity and therapeutic properties remain underexplored. This study aimed to characterize a symbiotic bacterium isolated from R. magnifica collected from Samalona Island, Indonesia, and to evaluate its multi-therapeutic potential. METHODS: Strain AZMABM HM27 was characterized using whole-genome sequencing, functional annotation, biosynthetic gene cluster prediction, molecular docking, and in vitro bioactivity assays. RESULTS: Phylogenetic and genome-based analyses confirmed AZMABM HM27 as Brevibacterium sanguinis, with an OrthoANI value of 97.37% and a dDDH value of 76.50% against the type strain. The genome comprises a 3,834,082 bp chromosome encoding 3,362 protein-coding genes, including 95 genes involved in secondary metabolite biosynthesis. Five biosynthetic gene clusters were predicted, including those associated with ectoine, terpene, and siderophore production. The crude extract demonstrated antioxidant activity (IC₅₀ = 0.87 mg/mL), anti-inflammatory activity (up to 60% inhibition), antidiabetic activity through α-glucosidase inhibition (up to 40% inhibition), and dose-dependent antiproliferative activity against MCF-7 breast cancer cells (74.10% viability at 1 mg/mL). Molecular docking identified a lead compound, 8,9,9,10,10,11-hexafluoro-4,4-dimethyl-3,5-dioxatetracyclo [5.4.1.0(2,6)0.0(8,11)] dodecane, with strong binding affinities to selected therapeutic targets. CONCLUSIONS: B. sanguinis AZMABM HM27 represents a marine symbiotic strain associated with R. magnifica and a promising source of bioactive compounds with antioxidant, anti-inflammatory, antidiabetic, and antiproliferative potential. Further purification, structural elucidation, and in vivo studies are warranted to validate its therapeutic potential.

Animals

Integrated Genome Mining and Bioactivity-Guided Isolation of Antimicrobial Peptides from Bacillus amyloliquefaciens BS4.

Bacterial resistance remains a critical global health challenge, driving the continuous search for novel antimicrobial agents. Bacillus amyloliquefaciens is a recognized repository of bioactive metabolites; however, its full biosynthetic potential requires integrated genomic and experimental validation. This study characterized the antimicrobial profile of B. amyloliquefaciens BS4 through a hybrid pipeline. Genome sequencing and de novo assembly revealed a 3.9 Mb chromosome with a G + C content of 46.14%. Functional annotation identified 3,887 coding sequences, including pathways for siderophore biosynthesis and a complete bacilysin biosynthetic cluster. BGC analysis using antiSMASH v7.1.0 and BAGEL4 identified 18 biosynthetic gene clusters, while similarity network analysis via BiG-SCAPE highlighted unique singleton BGCs, indicating untapped biosynthetic diversity. Although in silico screening via Macrel predicted two putative cationic antimicrobial peptides (AMPs), bioactivity-guided purification utilizing sequential RP-HPLC, and de novo sequencing revealed a distinct set of four active peptides. Notably, three of these sequences were identified as fragments derived from the BclA exosporium protein family, highlighting the structural proteome as a non-canonical source of antimicrobials. The purified fractions exhibited activity against M. luteus and E. coli, while displaying no significant hemolytic activity or cytotoxicity, even above the MIC values. Molecular docking further supported the interaction of these candidates with bacterial targets. Overall, this hybrid strategy effectively uncovers the antimicrobial complexity of BS4, revealing 'cryptic' peptide candidates with therapeutic potential.

Bacillus amyloliquefaciens BS4

Genomic diversity, inbreeding, and selection signatures in duroc, landrace, and yorkshire pigs from a long-term closed breeding system.

Duroc (DD), Landrace (LL), and Yorkshire (YY) are among the most widely used commercial pig breeds, having undergone intense long-term selection within closed breeding systems. This study presents a comprehensive genomic analysis of genetic diversity, inbreeding patterns, and selection signatures in DD, LL, and YY populations that have been subject to close breeding for over 15 years. Genomic and pedigree data were available for 1,088 animals (DD = 348, LL = 276, YY = 464), genotyped using the GenoBaits® Porcine 100 K SNP panel. Principal component analysis and genetic diversity metrics revealed distinct population structures among the three breeds. Pairwise genetic differentiation supported this pattern, with DD showing the greatest divergence from LL (0.34 ± 0.24) and YY (0.33 ± 0.24), while LL and YY were more closely related (FST = 0.22 ± 0.19). Linkage disequilibrium (LD) analysis further confirmed these differences, as DD exhibited the highest average r² (0.34), followed by LL (0.28) and YY (0.25). Within-breed genetic diversity metrics, including observed heterozygosity (HO: 0.37 in DD, 0.39 in LL, 0.38 in YY), expected heterozygosity (HE: 0.36 in DD, 0.37 in LL, 0.38 in YY), and minor allele frequency (MAF: 0.27 in DD, 0.28 in LL, 0.29 in YY), indicated greater genetic variability in LL and YY compared to DD. Runs of homozygosity (ROH) analyses revealed different patterns of autozygosity, with DD exhibiting more long ROH indicative of recent inbreeding, while YY harbored a higher number of short ROH, suggestive of more ancient demographic events. ROH-based inbreeding coefficients (FROH) consistently exceeded pedigree-based estimates (FPED) across all breeds, highlighting the presence of recent or unrecorded inbreeding that pedigree data may not fully capture. According to Generation Proxy Selection Mapping (GPSM), 17, 1, and 12 significant SNPs were detected in DD, LL, and YY, respectively. Functional annotation of ROH islands and GPSM-significant loci revealed both breed-specific and overlapping QTLs related to traits such as growth, reproduction, and carcass. In general, the findings of this study contribute to a deeper understanding of the genomic consequences of long-term closed breeding and provide reference information to support consideration of breeding strategies that balance continued selection for productivity with the maintenance of genetic diversity in modern commercial pig populations.

Animals

Mining metagenomes from extremophiles as a resource for novel glycoside hydrolases for industrial applications.

The exploration of metagenomes from extremophiles has emerged as a promising approach for discovering novel glycoside hydrolases (GHs) with potential industrial applications. Extremophiles, which thrive in harsh conditions such as high salinity, extreme temperatures, and acidic or alkaline environments, produce enzymes naturally adapted to function under these conditions. This unique adaptability makes them highly desirable for industrial processes requiring robust and efficient biocatalysts. These biocatalysts reduce reliance on harsh chemicals and energy-intensive processes, contributing to greener industrial operations. This review underscores the power of metagenomics in bypassing the need to culture large libraries of extremophiles in the lab. High-throughput sequencing and bioinformatics enable the identification of novel GH-encoding genes directly from environmental DNA. While metagenomic mining has yielded promising results, challenges such as the expression of extremophile-derived genes in mesophilic hosts, low activity yields, and scalability remain. Advances in synthetic biology and protein engineering could address these bottlenecks, enabling more efficient utilization of GHs. Additionally, integrating machine learning for predictive functional annotation may accelerate the identification of high-value candidates.

Glycoside Hydrolases

Novel HLA class I and II insights into the pathogenesis of systemic sclerosis-associated interstitial lung disease.

OBJECTIVES: Systemic sclerosis-associated interstitial lung disease (SSc-ILD) is the leading cause of mortality in systemic sclerosis (SSc), yet its genetic architecture remains incompletely understood. Therefore, given the key role of the major histocompatibility complex (MHC) in SSc, we aimed to perform a comprehensive MHC-wide association study in the largest SSc-ILD cohort to date. METHODS: We analysed 2412 patients with SSc-ILD⁺, 3550 patients with SSc-ILD⁻, and 15,076 controls of European ancestry from 10 international cohorts. After quality control, the MHC region was imputed, and inverse variance weighted meta-analysis was performed. Subsequently, conditional stepwise analyses, adjustment for antitopoisomerase autoantibody (ATA) status, and functional annotation of significant single-nucleotide polymorphisms were performed. Finally, we constructed a composite score combining genetic, clinical, and demographic variables to predict SSc-ILD. RESULTS: After conditional analysis, we detected 12 significant associations within class I and class II human leukocyte antigen (HLA) genes. ATA adjustment reduced the significance of class II HLA variants, whereas class I HLA variants remained unaffected. Finally, the built composite score had an area under the curve of 0.754, significantly outperforming the models including any of the variables alone. CONCLUSIONS: In this study, we identify genetic mechanisms underlying SSc-ILD that support the potential implication of CD8+ T cells and ATAs in its pathogenesis. Moreover, we also demonstrate the enhanced efficacy of integrating genetic information into predictive models to detect patients at high risk of SSc-ILD. These findings provide new insights into disease pathogenesis and suggest potential biomarkers and therapeutic targets for improved patient management.

Humans

Shared genetic architecture and therapeutic targets across paediatric immune-mediated diseases.

OBJECTIVES: Paediatric-onset immune-mediated inflammatory diseases (IMIDs), including juvenile idiopathic arthritis and related rheumatic diseases, remain genetically undercharacterised. We aimed to define shared and category-specific genetic architecture across paediatric IMIDs, compare signals with adult IMIDs, and identify therapeutic opportunities. METHODS: We analysed 24 paediatric IMIDs classified as autoimmune, polygenic-autoinflammatory, mixed-pattern, or allergic. Genome-wide association analyses included 18,086 cases and 131,019 controls of European ancestry. We estimated single nucleotide polymorphism (SNP)-based heritability, genetic correlations, and polygenic overlap; performed subset-based meta-analysis; and conducted functional annotation, gene prioritisation, pathway and protein network analyses, adult-IMID comparison, and drug-target prioritisation. RESULTS: SNP-based heritability ranged from 28.9% for allergic IMIDs to 61.9% for autoimmune IMIDs. Genetic correlation and polygenic modelling supported partial sharing across categories with category-specific components. Meta-analysis identified 39 genome-wide significant loci outside the Major Histocompatibility Complex (MHC) region, including 15 previously unreported loci; 19 loci were shared between categories. Gene-prioritisation and protein interaction analyses identified a core MHC-centred antigen-presentation network, with category-enriched modules involving complement, innate/barrier pathways, epithelial biology, and type 2 immunity. Enriched pathways included nuclear factor κB signalling, T helper 17 related pathways, Janus kinase-signal transducer and activator of transcription signalling, programmed cell death protein 1/programmed death‑ligand 1, cytotoxic T‑lymphocyte associated protein 4 regulation, and osteoclast differentiation, several of which are relevant to rheumatic diseases. Paediatric IMIDs shared broad polygenic architecture with adult IMIDs, whereas top-ranked genes converged strongly with adult rheumatic diseases. Priority Index analysis identified 178 high-scoring genes, including 43 approved or investigational IMID drug targets. CONCLUSIONS: Paediatric-onset IMIDs share core pathways with adult forms but exhibit distinct genetic architecture shaped by age-specific immune and neurodevelopmental biology. These findings provide a genomic framework for paediatric precision medicine, guiding classification, risk prediction, and therapeutic development.

Humans

Can Psychiatric Genetics Advance Without Incorporating a Life Course Perspective?

Psychiatric disorders unfold over the life course; however, genomic studies of these conditions overwhelmingly rely on phenotypes collected at a single time point, often in adulthood. Therefore, genome-wide association studies (GWASs) of psychiatric conditions may miss genetic variants with time-varying relevance to etiology, prevention, and treatment, such as those that influence trajectories of symptoms and behaviors, age at onset, course of treatment response, and the co-evolution of comorbidities. With recent advances in longitudinal biobanks and analytic tools, we posit that incorporating a life course perspective in psychiatric genetics will enable critically relevant insights into each of these areas of investigation. We propose that the current inconsistent portability of polygenic scores across age groups can be reconciled through the design of carefully considered longitudinal GWASs in age-diverse samples. Pioneering longitudinal GWASs in psychiatry have revealed novel genomic signals associated with time-dependent phenotypes that are distinct from those influencing lifetime diagnosis, suggesting that the study of longitudinal phenotypes will complement cross-sectional approaches and empower biological and therapeutic discoveries. Advances in post-GWAS functional annotation resources and analytic approaches now enable us to contextualize the genetic contributions to psychiatric disorders as dynamic age- and exposure-dependent processes. Although longitudinal GWASs pose unique challenges with regard to data availability, selection bias, and missing data, integrating temporality into psychiatric genetics at scale is now attainable and promises to reveal novel biology and therapeutic opportunities for psychiatric conditions.

Cohort study

Sex- and development-specific transcriptomic profiling of venom and silk genes in the wolf spider Pardosa astrigera provides insights into ecological adaptation and predatory strategies.

Spider venom and silk glands are two major secretory systems that contribute to prey capture, defense, and reproduction, but their sex- and development-specific molecular regulation in wandering wolf spiders remains poorly understood. Here, the transcriptome of Pardosa astrigera, an important agricultural natural enemy in China, revealed significant sex- and development-associated molecular differentiation among adult females, adult males, and spiderlings. A total of 100,025 unigenes were obtained, of which 23,852 were functionally annotated, providing a comprehensive transcriptomic resource for this species. Differential expression patterns showed marked variation among groups, with 531, 1792, and 832 DEGs detected in PAF vs PAS, PAM vs PAS, and PAF vs PAM, respectively. These genes were mainly associated with metabolic, oxidation-reduction, cuticle development, MAPK signaling, and lysosome pathways. Fifteen co-expression modules revealed distinct expression patterns. The turquoise, pink, yellow, and red modules were development-related, whereas the blue module was male-biased. Venom- and spidroin-related genes were distributed across multiple modules, suggesting coordinated regulation. Overall, 42 venom peptides, 21 venom proteins, and 11 spidroins were identified. Representative genes showed strongly biased expression, including spiderling-biased U3_Pp1a and U5_Pp1e, female-biased U4_Pp1a, and male-biased SMase D_108750 and PaTuSp_108466. These findings reveal sex- and development-biased expression patterns of venom- and silk-related candidate genes in P. astrigera and may provide molecular insights into ecological adaptation and predatory strategies in wandering wolf spiders.

Animals

A deep metagenomic atlas of Qinghai-Xizang Plateau lakes reveals their microbial diversity and salinity adaptation mechanisms.

The Qinghai-Xizang Plateau (QXP), harboring the planet's highest density of plateau lakes, offers an exceptional biogeographic environment for studying extremophilic microbial communities and their adaptation to salinity. Through deep metagenomic sequencing, we construct the Qinghai-Xizang Lake Sediment Genome (QXLSG) catalog, a high-resolution genomic catalog comprising 5,866 metagenome-assembled genomes (MAGs), 58.16 million non-redundant protein encoding genes, and 19,008 biosynthetic gene clusters. Notably, 80.78% of the 2,742 species-level MAGs represent undescribed taxa, significantly expanding the known microbial diversity. Salinity emerges as the primary environmental factor influencing microbial community. Functional annotation highlights that the "salt-out" strategy, particularly the uptake of glycine betaine, is the main mechanism for salinity tolerance. This strategy is prevalent in both hypersaline lake communities and the dominant microbial phyla. Overall, this study provides a crucial genetic resource for future bioprospecting and deepens our understanding of the fundamental mechanisms of microbial adaptation to extreme saline environments.

Lakes

Comparative genomic analysis of Streptococcus parasuis and Streptococcus suis reveals mobile element-associated enrichment of antimicrobial resistance and lack of detectable same-MGE colocalization with virulence-associated genes within stable species boundaries.

Streptococcus suis is a major porcine pathogen and a zoonotic agent that causes meningitis and septicemia in humans. Streptococcus parasuis, a recently recognized close relative, remains poorly characterized with regard to its clinical significance and genomic features. In this study, we generated a single-contig closed genome assembly with genome-wide DNA methylation profiles for S. parasuis strain A1, isolated from a diseased pig in Xinjiang, China, and complemented in silico genomic predictions with isolate-level experimental validation of antimicrobial resistance (AMR) genotypes, virulence genotypes, and phenotypic susceptibility for this reference strain. Using this high-quality genome as a reference anchor, we performed comparative genomic analyses across 195 streptococcal genomes, comprising 15 S. parasuis and 180 S. suis strains, to distinguish genome-level co-occurrence of resistance and virulence determinants from their physical colocalization on the same mobile genetic element (MGE).Species boundaries remained clearly delineated at the genomic level, with a median interspecies average nucleotide identity (ANI) of approximately 86.0%, compared with intraspecies ANI medians of 97.5% for S. parasuis and 96.2% for S. suis. Pangenome analysis identified 12,693 gene clusters, of which 1086 were core clusters, and functional annotation revealed significant differences in accessory gene repertoires between the two species. Within this stable genomic framework, S. parasuis genomes carried a higher AMR gene burden; strain A1 harbored 10 AMR genes, multiple virulence-associated genes, three genomic islands, and eight prophage regions. For strain A1, PCR validation confirmed six AMR genes and six virulence genes, and disk diffusion testing demonstrated a multidrug-resistant phenotype consistent with the genotypic profile.Among 235 predicted mobile elements, 19 harbored AMR genes and seven carried Virulence Factor Database (VFDB) homologs, but none carried both categories simultaneously. This finding reflects a lack of detectable same-MGE colocalization under the applied annotation and assembly framework; it should not be interpreted as evidence of biological physical decoupling. Under a random-placement model, the expected number of co-carrying regions was only 0.57, and the probability of observing zero co-carrying regions was P = 0.55. This negative result should be interpreted with caution, given the limited number of cargo-bearing regions and the predominantly draft status of most genomes. Furthermore, the A1 genome contained multiple restriction-modification systems, showed depletion of several methylation motif families in mobile regions, and had limited CRISPR spacer matching evidence, suggesting prior exposure to the relevant sequence space. None of the genomes met our predefined criteria for whole-genome convergence.Collectively, our results support a model in which S. parasuis accumulates AMR-related genes in a modular fashion via mobile elements within stable species boundaries, with no detectable same-MGE colocalization of AMR and virulence determinants under our analytical pipeline. These findings imply that AMR surveillance strategies for this species should prioritize tracking mobile genetic elements rather than inferring wholesale genomic convergence toward S. suis.

Streptococcus suis

Uncovering hub genes and key pathways responsive to drought stress in rice via meta-analysis of transcriptomic data.

Drought stress presents a formidable threat to global rice cultivation, triggering complex molecular responses that impact plant growth and productivity. To decipher the underlying gene expression dynamics, we performed a comprehensive meta-analysis of transcriptomic datasets derived from drought-tolerant rice genotypes. Via microarray data from three independent studies, we identified a set of consistently expressed differentially expressed genes (DEGs) under drought conditions. Integration of functional annotation tools, including GO and KEGG pathway enrichment, revealed key biological processes and signaling cascades involved in stress mitigation, such as ABA signaling, protein folding, and photosynthesis suppression. Protein-protein interaction (PPI) network construction, followed by hub gene identification via maximal clique centrality (MCC), highlighted pivotal regulators including LEA proteins, dehydrins, HSP70, and several transcription factors. Machine learning approaches further prioritize potential biomarkers, with Random Forest models achieving high classification accuracy and pinpointing key predictive genes. Chromosomal localization analysis provided spatial insights into the distribution of these hub genes, whose expression patterns were further compared against qRT-PCR data from previously published studies. This integrative approach identifies candidate genomic markers and mechanistic insights that may support future breeding strategies for drought-tolerant rice, pending experimental validation.

Cytoscape