Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “phylogeny”

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.

88 records · Page 2Linked to original sources
PubMed2026

In Situ Hybridization and RT-PCR Detection of Nervous Necrosis Virus in Fourfinger Threadfin, Eleutheronema tetradactylum, in Taiwan.

Between April and July 2020, suspected outbreaks of nervous necrosis virus (NNV) infection were observed in fourfinger threadfin (Eleutheronema tetradactylum) fingerling hatcheries in Pingtung County, southern Taiwan. Affected fish exhibited spiral swimming behaviour and abdominal distension, resulting in mortality rates between 50% and 100%. Histopathological examination showed severe vacuolation in the brain and ocular tissues, with large oval and/or irregular basophilic cytoplasmic inclusion bodies in the brain. Phylogenetic analysis of the viral replicase (RNA1) and capsid protein (RNA2) genes revealed high nucleotide sequence identities among the isolates in this study, with sequence similarity rates of 96.9%-99% for RNA1 and 98.2%-99.0% for RNA2 compared to RGNNV reference strains available in the NCBI GenBank database. This is the first detection of betanodavirus in fourfinger threadfin in Taiwan, using RT-PCR and ISH. A positive correlation between elevated water temperatures and disease severity indicates the need for year-round surveillance and genomic analysis to clarify the epidemiology of FTNNV. The data suggest that infected eggs may facilitate the vertical transmission of Betanodavirus. Crucially, utilising virus-free broodstock, alongside routine health screening and environmental control, is essential for mitigating NNV risks in fourfinger threadfin aquaculture.

Animals↗
PubMed2026

Genome-wide identification and characterization of ABC transporters and their expression in response to saline-alkaline stress and WSSV infection in Fenneropenaeus chinensis.

ATP-binding cassette (ABC) transporters play crucial roles in stress responses across organisms, yet their functions in Fenneropenaeus chinensis remain largely unknown. In this study, we identified 42 FcABC genes (FcABCs) in the F. chinensis genome and analyzed their phylogenetic relationships, gene structures, and chromosomal distributions. Phylogenetic analysis grouped the FcABCs into eight subfamilies (ABCA-ABCH), with conserved motif and domain compositions within each subfamily. Expression analysis showed that several FcABC genes, including FcABCG5, FcABCA1, and FcABCC3, were significantly induced under saline-alkaline stress in gill and hepatopancreas tissues. In contrast, most FcABCs were downregulated after WSSV challenge, though a subset (e.g., FcABCB1, FcABCC1) exhibited early upregulation. Functional validation via RNA interference demonstrated that knockdown of FcABCG5 increased shrimp mortality under saline-alkaline stress. Cis-regulatory element analysis revealed an enrichment of stress- and immune-related elements in FcABC promoters. Protein-protein interaction network predictions indicated potential roles for FcABCs in cholesterol metabolism and organic anion transport. Our findings provide insights into the roles of FcABC genes in stress adaptation and immune defense, offering candidate genes for the breeding of stress-resistant shrimp varieties.

Animals↗
PubMed2026

Identification of aquaporin (AQP) genes in the noble scallop Chlamys nobilis and characterization of their expression under low-temperature stress.

Aquaporins (AQPs) are transmembrane channel proteins essential for water homeostasis and cellular stress responses. In marine bivalves, their roles in cold tolerance remain poorly understood despite frequent winter mortality events in aquaculture. Here, we identified nine AQP genes in the genome of the economically important noble scallop Chlamys nobilis. Phylogenetic analysis revealed strong conservation with other bivalve AQPs, and structural features, including conserved NPA motifs and ar/R selectivity filters, support their canonical water/glycerol transport functions. Tissue-specific expression profiling showed predominant enrichment in osmoregulatory tissues (gills, intestine) and gonads. Under both chronic and acute low-temperature stress from 23 °C to 9 °C, most CnAQP genes exhibited transient upregulation followed by suppression. Notably, CnAQP4 displayed sustained upregulation, implicating it as a key mediator of long-term cold adaptation. Promoter analysis further revealed abundant cis-elements linked to growth and development as well as immune regulation. Our findings provide the first comprehensive characterization of the AQP family in C. nobilis, highlighting its critical role in maintaining cellular integrity during cold stress and offering molecular targets for selective breeding of cold-tolerant scallop strains.

Animals↗
PubMed2026

Beyond species trees: pervasive gene flow limits phylogenomic resolution in the diversification of Juniperus from the Qinghai-Tibet Plateau.

Understanding how lineages diversify despite persistent ancestral polymorphism and recurrent gene flow remains a central challenge in evolutionary biology. Juniperus distributed across the Qinghai-Tibet Plateau provide an ideal system for addressing this question because repeated geological uplift and climatic oscillations have likely promoted cycles of lineage divergence, range shifts, and secondary contact. Here, we combined approximately 1.08 million genome-wide SNPs from 164 individuals representing thirteen Juniperus lineages with phylogenomic datasets comprising 3,381 nuclear single-copy genes and nearly complete plastomes. We detected extensive phylogenomic discordance and cytonuclear incongruence across genomic datasets. Topology weighting, coalescent simulations, quartet-based tests, and analyses of gene flow and reticulation collectively support the interpretation that these patterns were shaped by the combined effects of prolonged incomplete lineage sorting and gene flow during lineage diversification. Ecological niche analyses further provide a spatial and climatic context in which environmentally similar lineages may have had greater opportunities for secondary contact during historical range shifts. Collectively, our results reveal that the evolutionary history of Qinghai-Tibet Plateau Juniperus is characterized by reticulate diversification rather than strictly bifurcating evolution, and demonstrate how genome-wide discordance can provide biological insights into the evolutionary processes underlying lineage diversification.

Gene Flow↗
PubMed2026

Genome-wide identification and expression analysis of the CREB/ATF family and its potential role in melanogenesis in the Manila clam (Ruditapes philippinarum).

Ruditapes philippinarum is an economically important bivalve species in China, and shell color is a trait of ecological and commercial significance. Melanin is a key determinant of shell color, and members of the CREB/ATF family have been reported to participate in melanogenesis in other organisms. In this study, members of the CREB/ATF family were systematically identified at the whole-genome level based on genomic and transcriptomic datasets, followed by analyses of their phylogenetic relationships, gene structures, and expression patterns. A total of six CREB/ATF family members were identified and classified into five subfamilies. Expression profiling and RT-qPCR validation revealed that most CREB/ATF genes were highly expressed in the mantle and displayed clear differences among shell-color phenotypes. Except for RpATF4 and RpCREBZF, most members exhibited relatively high expression levels in dark-colored shell strains, particularly in black and zebra-striped clams. Moreover, most genes showed low expression during early embryonic and larval stages but increased expression at the single-siphon spat and juvenile stages. These results suggest that the CREB/ATF family may be involved in melanin-associated shell-color regulation in R. philippinarum, providing important candidate genes and a theoretical basis for further elucidating the molecular mechanisms of shell-color formation in mollusks.

Animals↗
PubMed2026

Complete mitochondrial genomes of eight cyclophyllidean tapeworms: genome pattern and phylogenetic analysis.

Cyclophyllidean tapeworms are widespread parasites of significant medical and veterinary importance. However, mitochondrial (mt) genomic resources for cyclophyllideans from China, particularly those recovered from wildlife hosts, remain comparatively limited. In this study, we sequenced and characterized the complete mt genomes of eight cyclophyllidean isolates collected from diverse wild and domestic hosts in China, including two Hymenolepis sp. isolates and two Raillietina sp. isolates from China, and four additional isolates of previously sequenced Taenia species. The circular mt genomes ranged from 13,387 to 14,021 bp in length, encoding 36 typical genes with variable non-coding regions. Comparative analysis revealed highly conserved gene composition and mostly conserved mt architecture, with localized rearrangement patterns detected among the cyclophyllidean lineages examined. In particular, all sampled Taeniidae exhibited a consistent trnL1-trnS2 arrangement, whereas the examined non-Taeniidae families showed the trnS2-trnL1 arrangement, confirming and extending, across additional wildlife-associated isolates, a previously proposed family-associated gene-order marker within Cyclophyllidea. Phylogenetic analyses based on concatenated amino acid sequences of the 12 protein-coding genes placed the eight isolates within their expected families, in topologies broadly consistent with previous mitogenomic studies. These data provide additional Chinese mitogenomic references, especially for underrepresented wildlife-associated isolates, and support family-associated gene-order patterns in Cyclophyllidea.

Animals↗
PubMed2026

Genome-wide scans reveal candidate genes associated with wing morph differentiation in Tetrix japonica.

Wing dimorphism is an important dispersal-related trait in insects, but its genomic basis remains poorly understood in pygmy grasshoppers. Here, we integrated genome-wide single-nucleotide polymorphism (SNP) analyses, population structure inference, selection scans, and functional annotation to investigate genomic differentiation between long- and short-winged Tetrix japonica. Principal component analysis (PCA), ADMIXTURE, and phylogenetic analyses revealed weak genome-wide separation between morphs, indicating differentiation on a largely shared genetic background. Genome-wide scans based on the fixation index (FST), nucleotide diversity ratios, and Tajima's D, using 50-kb non-overlapping windows and empirical top-5% outlier thresholds, identified multiple candidate regions across seven chromosomes. The broader long- and short-winged candidate sets spanned 9.35 Mb and 9.37 Mb and directly overlapped 82 and 77 genes, respectively. Candidate genes were associated with signaling/hormone regulation, membrane transport, metabolism, cytoskeletal organization, extracellular matrix structure, and development. Short-winged candidate genes were significantly enriched for ABC-type transporter activity and ATP hydrolysis activity. Because all individuals originated from a single laboratory-maintained population with weak genome-wide structure, these regions should be regarded as candidate loci from a screening-stage analysis that require validation in independent populations and by functional assays, rather than as confirmed targets of selection.

Animals↗
PubMed2026

Complex evolutionary history of Rosales mediated by extensive incomplete lineage sorting and hybridization.

The angiosperm order Rosales still represents a major challenge for phylogenetic reconstruction. Although its circumscription is now well-defined, phylogenetic relationships among families are still uncertain. Here, we used nuclear, plastid, and mitochondrial genomic data from 33 species representing all nine families to further clarify interfamilial relationships and the group's evolutionary history. We detected significant phylogenetic conflict among the three datasets. Further analyses at the nuclear level identified incomplete lineage sorting (ILS) as the main cause of unstable phylogenetic positions among families. The discordant placements of Rhamnaceae and Elaeagnaceae based on plastid and mitochondrial data are caused by ancient hybridization events, potentially involving differences in organellar inheritance. Our molecular dating confirms earlier suggestions that the ancient rapid diversification of the three Rosaceae subfamilies could be the main reason for the difficulties in resolving their phylogenetic relationships. Our findings provide new insights into the interfamilial relationships of Rosales and demonstrate that the evolutionary history of this order was shaped by ancient and rapid radiation as well as extensive ILS and reticulate evolution. They also suggest that previous attempts to clarify interfamilial relationships in this order were hampered by combining nuclear and organellar sequence data, leading to inconsistent topologies observed across earlier studies.

Phylogeny↗
PubMed2026

Viral community in Aspergillus spp. isolated from commercially available fermented dried bonito.

Katsuobushi is a traditional processed seafood product used in Japanese-style cooking, and when it is produced through fermentation by fungi, it is called karebushi. The fungi involved in katsuobushi fermentation are collectively referred to as katsuobushi molds. We previously discovered seven novel viruses from katsuobushi molds and determined their genome sequences. However, our previous explorations used only nine fungal strains available from culture collections, leaving the diversity of viruses infecting fungi involved in katsuobushi fermentation unclear. Therefore, in this study, we aimed to isolate fungi from commercially available karebushi and clarify the prevalence of viruses in the isolates. Karebushi produced by three manufacturers was obtained, and 30 fungal strains (including Aspergillus spp.) were isolated from each. Double-stranded RNA (dsRNA) fractions were prepared from the mycelia of the isolated strains. Electrophoresis suggested that a relatively high proportion of the isolates harbored dsRNA elements consistent with RNA virus infection (30-70% per manufacturer; 59% overall). Furthermore, dsRNA sequencing identified four novel viruses in isolates of Aspergillus chevalieri and Aspergillus montevidensis: a beny-like virus, a gammapartitivirus, a narnavirus, and a victorivirus, in addition to two previously reported viruses. Notably, this represents the first report of a beny-like virus in Aspergillus spp. This study provides insights into the diversity of viruses infecting fungi involved in katsuobushi fermentation.

Aspergillus↗
PubMed2026

Comparative genomics and full-length transcriptome profiling of wing morphs in Tetrix grossus (Orthoptera: Tetrigidae).

Wing polymorphism represents a paradigmatic dispersal-reproduction trade-off, yet its molecular basis remains uncharacterised in the phylogenetically distant pygmy grasshoppers (Tetrigidae). Here we integrate comparative genomics across ten orthopteran species with full-length transcriptomics of long-winged (FL) and short-winged (FS) Tetrix grossus. OrthoFinder recovered 118 orthogroups specific to T. grossus. Against a backdrop of pronounced gene-family contraction (36 expansions versus 222 contractions; net -186, mirrored at the ancestral Tetrix node, +37/-140), we identified an ancestral, Tetrix-specific expansion of hormone-regulation (12 genes; fold enrichment 7.93) and lipid/carbohydrate-metabolic families organised into syntenic clusters, alongside 513 positively selected genes enriched for integrin-mediated cell adhesion (6 genes), a process relevant to epithelial and appendage morphogenesis. Full-length transcriptomics of one long-winged (FL) and one short-winged (FS) adult female detected 7530 (FL) and 7515 (FS) expressed genes, with 794 FL- and 776 FS-restricted transcriptome-derived SNP-associated genes. The FL morph was enriched for an EGFR/Ras-Rho developmental-patterning axis and neuromuscular flight genes, whereas the FS morph was enriched for insulin/peptide-hormone response and growth-regulatory loci. Overall, we present genomic resources and testable hypotheses concerning the evolution and regulation of wing morphs in Tetrigidae rather than a validated genetic architecture of wing-morph determination.

Animals↗
PubMed2026

Comparative genomic and proteomic analysis reveals orthogroup structured evolution of tick protease inhibitors.

Protease inhibitors (PIs) play central roles in regulating endogenous proteolysis and host-parasite interactions in ticks. However, the evolutionary architecture underlying their diversification across tick lineages remains insufficiently resolved. Here, we performed a genome-wide comparative analysis of predicted proteomes from 14 tick species to systematically characterize PI repertoires. In total, 4931 putative PIs were identified and grouped into 20 families using the MEROPS classification system. Further, PI families such as Antistasin, WAP-type, and Pacifastin, which have not previously been systematically reported in tick genomes, were classified. Orthogroup inference demonstrated that PI expansion is structured at the level of evolutionary lineages rather than uniformly across families. By stratifying orthogroups according to duplication burden and taxonomic conservation, we identified a broadly conserved single-copy core under strong purifying selection. Motif level analysis of serpin reactive center loops further revealed conservation of inhibitory specificity within single copy orthogroups and diversification of key functional residues in duplication-associated lineages. Integration of secretion prediction and tissue-resolved proteomics from Hyalomma anatolicum and Rhipicephalus microplus demonstrated that evolutionary stratification is reflected at the protein level. Together, these findings provide an orthogroup-resolved evolutionary framework linking duplication dynamics, molecular evolution, and tissue-level protein deployment. This integrative approach offers a systematic basis for prioritizing conserved and diversified PI lineages for future functional and anti-tick intervention studies.

Animals↗
PubMed2026

Fructophilic lactic acid bacteria as a window into multi-scale convergent evolution.

Fructophilic lactic acid bacteria (FLAB) are a group of lactic acid bacteria with unique growth characteristics, that is, poor growth on glucose. Their growth is enhanced in the presence of fructose or external electron acceptors. These organisms inhabit fructose-rich environments such as flowers, fruits, and pollinating insects, particularly honey bees. Apilactobacillus spp. and Fructobacillus spp. are representatives of FLAB, although they belong to phylogenetically distant clades. These organisms commonly possess markedly small genomes with a low number of coding DNA sequences. Furthermore, their genomes are characterized by a markedly reduced number of genes involved in carbohydrate transport and metabolism. Genome reduction in FLAB reflects convergent adaptation to fructose-rich environments rather than general genome streamlining. The two distinct FLAB genera, Fructobacillus and Apilactobacillus, independently lost more than 100 genes in statistically similar orders. In contrast, genes involved in carbohydrate and amino acid metabolism exhibited reversed orders of loss between the two genera. Furthermore, FLAB genomes lack an intact bifunctional alcohol/aldehyde dehydrogenase gene (adhE), which causes their poor growth on glucose. A comparative genomic study suggested the evolutionary process underlying adhE gene decay during adaptation to the fructose-rich environments, including pollinating insects. In conclusion, FLAB represent a unique example of habitat-driven convergent reductive evolution that can be investigated across multiple biological scales - from individual genes to whole genomes - in the diverse LAB group with a wide range of habitats, and partially share the fructophilic evolution with eukaryotic yeasts found in fructose-rich habitats.

Fructose↗
PubMed2026

Structural and tissue-specific organisation of endocrine Fgf19 and Fgf21 signalling in rainbow trout.

Endocrine fibroblast growth factors (FGF19 subfamily) play a key role in regulating metabolic homeostasis in vertebrates. However, their functional diversification in salmonids remains poorly understood. In this study, we conducted an integrative characterisation of Fgf19 and Fgf21 signalling in rainbow trout (Oncorhynchus mykiss) by combining phylogenetic, structural and expression analyses. Phylogenetic analyses revealed the conservation of single fgf19 and fgf21 genes, despite the extensive expansion of receptors post-Ss4R (salmonid-specific fourth-round whole genome duplication). Structural modelling and molecular dynamics simulations demonstrated the stable interactions of both ligands to multiple Fgfr isoforms, with receptor-specific energetic profiles and conserved core interaction residues. Tissue expression profiling revealed clear differences from mammalian models, such as predominant hepatic fgf19 expression and the absence of hepatic fgf21 under basal conditions. In addition, there were complex and tissue-dependent distributions of fgfr and klotho transcripts. These findings support a receptor-driven diversification model of endocrine Fgf signalling in salmonids, suggesting enhanced endocrine plasticity associated with the retention of receptors following post-genomic duplication. Taken together, our findings provide new insights into the structural and regulatory organisation of endocrine Fgf signalling, as well as its potential role in metabolic regulation in rainbow trout.

Animals↗
PubMed2026

Identification of the BrSK gene family in flowering Chinese cabbage and functional characterization of BrSK2 subfamily involvement in heat stress.

Glycogen synthase kinase 3 (GSK3) kinases are evolutionarily conserved regulators of plant development and stress signaling, yet their contributions to thermotolerance in cool-adapted Brassica crops remain poorly understood. Here, we identified 16 BrSK genes in the Caixin (Brassica rapa ssp. chinensis var. parachinensis) genome, all harboring intact catalytic motifs indicative of functional kinase activity. Spatiotemporal expression profiling revealed preferential accumulation of BrSK transcripts in stem apices and floral organs during reproductive transition, while promoter analysis identified abundant heat- and abiotic stress-responsive cis-elements. Under heat stress, BrSK21, BrSK22, and BrSK23 displayed striking genotype-specific expression dynamics. BrSK21/22/23 transcripts were stably suppressed in the heat-tolerant cultivar '49-19' but transiently declined before rapidly rebounding in the heat-sensitive 'Liuye 50', mirroring RNA-seq profiles. Protein-protein interaction assays (Y2H, BiFC, and LCI) demonstrated specific associations between BrSK kinases and BrHSFA1. Functional validation via VIGS revealed that silencing of BrSK21 significantly enhanced thermotolerance, with triple silencing of BrSK21/22/23 conferring additive protection, indicating functional redundancy within the BrSK2 subfamily. Collectively, these findings establish the BrSK2 subfamily as negative regulators of heat tolerance in Caixin, likely via modulation of BrHSFA1 expression. This work identifies high-priority targets for molecular breeding of climate-resilient Brassica vegetables.

Plant Proteins↗
PubMed2026

Pangenome-wide identification and expression analysis of the chalcone synthase (CHS) gene family in five yellowhorn spp.

Chalcone synthase (CHS) is a pivotal enzyme in flavonoid biosynthesis involved in plant development, defense, and secondary metabolism. Xanthoceras sorbifolium (yellowhorn) is a medicinal and ornamental species with high resistance to environmental stresses, but its CHS gene family remains uncharacterized. We performed a pangenome-wide identification of CHS genes across five yellowhorn genomes (Xzs4, Xwf8, Xjg, Xg11, and Xzg2). Across the five yellowhorn genomes, 27 CHS genes were identified and classified into four core pangenes, present in all five genomes, and two dispensable genes, present only in a subset of genomes. Phylogenetic analysis grouped these genes into three major clades, and chromosomal mapping and duplication analyses identified four tandemly duplicated gene pairs under purifying selection. The analyses of conserved structural features, including protein motifs and exon-intron organization, together with promoter cis-regulatory elements and gene ontology annotation, further indicated the potential involvement of CHS genes in flavonoid biosynthesis and stress-responsive mechanisms. Gene expression profiling identified significant upregulation of Xg11_CHS1 and Xg11_CHS3 under cold and drought stress, with tissue-specific expression patterns. These findings provide valuable insights into the evolution, functional diversification, and stress-responsive roles of the CHS gene family, identifying candidate genes for future studies targeting stress tolerance and flavonoid biosynthesis in yellowhorn.

Acyltransferases↗
PubMed2026

The Complete Mitochondrial Genome of a Newly Recorded Chinese Species of Diglyphus sabulosus (Hymenoptera: Eulophidae) and Insights into Its Phylogenetic Position.

Diglyphus Walker, 1844 is an economically important genus which many species acting as biocontrol agents against agromyzid leafminer pests, but there is a lack of mitogenomic data on the evolutionary relationships within this genus, hindering a comprehensive understanding of its evolutionary history. We used traditional morphological methods to identify species, and present the first complete mitochondrial genome sequence and characterization of features of Diglyphus sabulosus and further infer its phylogenetic position based on the amino acid sequences of 13 protein-coding genes (PCGs). The complete mitochondrial genome of D. sabulosus is 15,690 bp in length, including 13 PCGs, 22 transfer RNA genes, 2 ribosomal RNA genes and a control region. The AT content of the whole genome sequence was 81.0%, indicating a significant AT bias. All protein-coding genes have the typical ATN as the start codon and TAA as the stop codon. Phylogenetic analysis inferred from the amino acid sequences of 13 PCGs revealed that all species within the family Eulophidae constituted a monophyletic clade, supporting the monophyly of this family. D. sabulosus and D. poppoea form a well-supported sister group, representing the species with the closest phylogenetic relationship within the analyzed taxa. In this study, the mitogenome structure was analyzed and the taxonomic status of D. sabulosus was clarified, thus providing a theoretical basis for understanding the phylogenetic relationships of Diglyphus.

Animals↗
PubMed2026

Tigecycline-resistant Staphylococcus in waiting pens of a pig slaughterhouse: genomic insights into a food safety alert.

BACKGROUND: The waiting pens of slaughterhouses represent a critical control point in the 'farm-to-fork' continuum, yet their role in the emergence and dissemination of antimicrobial resistance remains understudied. This study investigated tigecycline-resistant Staphylococcus (TRS) in these high-risk zones to assess their prevalence, resistance mechanisms, and transmission dynamics. METHODS: 400 samples were collected from the waiting pens of a pig slaughterhouse in Guangzhou, China. Antimicrobial susceptibility testing, whole-genome sequencing, phylogenetic analysis, and molecular cloning were employed to characterize resistance mechanisms and transmission patterns. RESULTS: 78 TRS strains were isolated and classified into three species, including S. borealis, S. ureilyticus, and S. pasteuri. These isolates exhibited multidrug-resistant phenotypes and carried new mutations in rpsJ and tet(M), which were functionally confirmed to reduce tigecycline susceptibility. Phylogenetic evidence demonstrated clonal transmission between pig farms and the slaughterhouse. The tet(M) gene was located within Staphylococcal cassette chromosome mec elements mediated by IS257, while tet(L) was carried by plasmids formed through IS256/IS257-mediated recombination. CONCLUSIONS: Waiting pens serve as crucial reservoirs for the amplification and dissemination of antimicrobial resistance. Our findings underscore the urgent need for enhanced biosecurity measures, improved waste management, and routine molecular surveillance in these high-risk zones to mitigate the spread of resistance along the food production chain.

Animals↗
PubMed2026

Genome-wide characterization of the TGF-β superfamily identifies bmp15, gdf9, and gsdf as sex-biased candidate regulators of gonadal differentiation in the synchronous hermaphrodite Plectropomus leopardus.

The transforming growth factor-β (TGF-β) superfamily plays conserved roles in vertebrate reproduction and gonadal sex differentiation. However, its genomic repertoire and sex-biased expression patterns remain unclear in the leopard coral grouper (Plectropomus leopardus), a species with synchronous hermaphroditism. Here, we performed a genome-wide identification of the TGF-β superfamily, identifying 42 genes from the chromosome-level genome. Phylogenetic and synteny analyses indicated that segmental duplication under purifying selection contributed to family expansion. Expression profiling across multiple tissues and four gonadal developmental stages (undifferentiated, 120 dph; early differentiated, 15 months; mature testis, 3 years; mature ovary, 3 years) identified eight gonad-enriched genes, among which bmp15 and gdf9 exhibited pronounced female-biased expression, with transcripts localized exclusively to the oocyte cytoplasm, particularly in stage II-III oocytes. In contrast, gsdf showed male-biased expression and was localized in spermatogenic cells of the testis. These reciprocal expression patterns indicate that bmp15/gdf9 and gsdf are candidate factors associated with gonadal sex differentiation. Our study provides the first comprehensive characterization of the TGF-β superfamily in P. leopardus and highlights bmp15, gdf9, and gsdf as candidate sex-differentiation factors in this hermaphroditic species.

Animals↗