Search PubMedSearch

SEARCH · Search PubMed

Results for “Phylogenetic relationship”

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.

303 records · Page 2Linked to original sources

Genome-wide characterisation of the myosin light chain gene family in Chinese perch (Siniperca chuatsi) and its expression patterns in muscle fibre types and injury response.

The Class II myosin light chain (myl) genes in Chinese perch (Siniperca chuatsi) have not yet been systematically characterised, and relationships with muscle fibre specification, development, and injury-associated remodelling remain unclear. In this study, fast and slow muscle fibres were initially distinguished using myofibrillar ATPase histochemistry. Subsequently, genome-wide mining identified 16 Class II myl genes, comprising eight essential and eight regulatory light-chain subunits. Their conserved-domain features, chromosomal distribution, phylogenetic relationships and expression profiles were analysed. Transcriptomic profiling showed that summed myl transcript abundance was higher in fast muscle than in slow muscle, accounting for 67.2% of the pooled myl transcript pool across the two muscle types (paired t-test, raw P = 0.036). mylpfa, myl1 and mylz3 were the major fast-muscle-associated genes, whereas myl10, myl2b and myl13 were preferentially expressed in slow muscle at the transcript level. These patterns support these genes as candidate fibre-type-associated expression markers. Developmental profiling identified stage-associated myl expression patterns, including a possible expression shift between mylpfb and mylpfa. In the descriptive injury-repair time course (d0-d7), FPKM profiles indicated that fast-muscle-associated genes (mylpfa, mylz3 and myl1) were lower at d1 and recovered by d3, whereas several slow-muscle-associated genes showed biphasic transcript-level increases. The slow-muscle-associated RLC gene mylpfb showed a delayed expression peak at d7. Notably, the embryonic isoform myl6l showed a modest increase from approximately 2 FPKM at d0 to 4-5 FPKM after injury, suggesting a possible injury-associated expression pattern that requires further validation. Together, these findings provide a genome-wide description of the Chinese perch myl gene family and identify candidate fibre-type-associated genes and descriptive injury-associated isoform expression patterns.

Animals

Genome-wide identification and expression profiling of CSP and OBP genes in Stictocephala bisonia reveals candidate genes potentially associated with insecticide response.

Stictocephala bisonia is an important invasive agricultural pest. Due to the frequent application of insecticides in its habitat, this species is under intense selection pressure. Chemosensory proteins (CSPs) and odorant-binding proteins (OBPs) are known to play key roles in insecticide resistance, but their specific functions in S. bisonia remain unclear. In this study, we identified a total of 22 SbisCSPs and 16 SbisOBPs based on the S. bisonia genome. To screen for candidate genes potentially linked to insecticide resistance, we adopted a multi-criteria screening strategy that integrated phylogenetic analysis, molecular docking with three insecticides, and tissue-specific expression profiling. Phylogenetic analysis identified several SbisCSPs and SbisOBPs clustering with genes known to be involved in insecticide resistance, serving as an initial evolutionary filter. Molecular docking results indicated that λ-Cyhalothrin exhibited the strong predicted binding affinity with most of SbisCSPs and SbisOBPs. Subsequent qPCR validation of seven prioritized candidates revealed distinct expression patterns: SbisCSP22 was highly expressed in adults and demonstrated strong binding affinity to all three insecticides tested, suggesting a potential role in mediating multi-insecticide response. Conversely, SbisCSP17 was significantly upregulated in larvae, clustered with genes known to mediate imidacloprid resistance, and exhibited strong binding affinity to imidacloprid. Given its larval-specific expression and the soil-dwelling behavior of larvae, we hypothesize that SbisCSP17 is a key candidate gene for larvae coping with soil-treated insecticides.

Animals

Mitochondrial DNA diversity in Ecuadorian populations: Recurrence of variant 16136 within haplogroup B2.

The identification of lineage-defining variants, frequently found in the coding region of mitochondrial DNA (mtDNA), is essential for refining haplogroup classification. Most mtDNA studies in South American populations have focused on the control region (CR), which has provided important insights into population structure and maternal lineage origins, although information needed for more robust phylogenetic resolution has been neglected. This study investigates the maternal genetic structure of Ecuadorian populations by combining CR and whole mitogenome analyses. Sequences from the mtDNA CR were obtained from 461 individuals (253 Mestizos and 208 Native Americans), while complete mitogenomes were sequenced for 127 individuals to improve phylogenetic resolution by identifying lineage-defining variants present in coding region. Most mtDNA haplogroups in the two population groups analyzed were of Native American origin (A2, B2, B4, C1, D1, D4), with significant differences in the distribution of specific lineages between them. Among Mestizos, African haplogroups (all within the L branches) and Eurasian haplogroups (H, K, R, U) were detected at low frequencies, whereas no African lineages were observed among Native Americans. The results obtained highlighted a heterogeneity within Ecuadorian populations that must be considered when developing mtDNA haplotype databases for forensic purposes. Whole mitogenome sequences enabled the identification of variants that refined haplogroup classifications, provided a more accurate reconstruction of the maternal genetic diversity, and improve the discrimination between Native American and Asian maternal lineages within haplogroup B4b.

Humans

Genome-wide SNP data support species boundaries in sympatric Polylepis Ruiz & Pav. (Rosaceae) species from Bolivia and Ecuador.

Species delimitation in the South American genus Polylepis is notoriously challenging due to high morphological similarity and phenotypic plasticity, likely driven by hybridization and gene flow. Previous phylogenetic studies suggested that genetic structure aligns more strongly with geography than with taxonomy, questioning existing species concepts and hampering conservation efforts. We used double-digest RAD sequencing (ddRADseq) to generate genome-wide SNP data for 11 Polylepis species sampled across multiple localities in Bolivia and Ecuador. Population genetic analyses, phylogenetic inference, and network approaches were combined to assess whether genetic structure aligns more closely with taxonomy or geography. Morphologically defined species formed largely cohesive genetic lineages across regions, with species identity explaining substantially more genetic variation than locality. While localized admixture and reticulation were detected among closely related taxa, widespread species showed strong genetic cohesion and clear separation from congeners. Our results indicate that the sampled Polylepis species from Bolivia and Ecuador maintain distinct genetic identities despite localized signals consistent with gene flow. This genome-wide support for current taxonomy highlights Polylepis as a valuable model for studying speciation under gene flow and indicates that multiple geographic sampling will be essential in reconstructing a robust phylogeny of the genus, with important implications for conservation planning in Andean montane forests.

Bolivia

Cytonuclear conflict and reticulate evolution in the Morelloid clade (Solanum, Solanaceae): Insights from genome skimming and network Phylogenomics.

The Morelloid clade (black nightshades) is one of the most strongly supported clades within the megadiverse Solanum genus. It comprises 76 globally distributed, non-spiny herbaceous and suffrutescent species. While often erroneously considered poisonous weeds, several species are economically important as orphan crops. The clade is closely related to tomato and potato but, due to a lack of focused breeding efforts, remains a putative reservoir of genetic diversity for crop improvement. Despite this potential, we lack fundamental knowledge on the evolution of the Morelloid clade. The group includes polyploid species with unknown parental origins-likely reflecting reticulate processes such as hybridization, introgression, and associated backcrossing events. Prior analyses have been unable to disentangle these processes, leaving the mechanisms underlying reticulate evolution in the Morelloid clade poorly understood. Here, we use genome skimming to produce a well-supported maximum likelihood plastid phylogeny from complete circularized plastomes and a coalescent-based species tree from combined Angiosperms353 and conserved ortholog set nuclear markers. Our dataset, composed of previously published data and deep genome skimming from herbarium samples, spans 26 Morelloid species. To investigate phylogenetic discordance, we used a nuclear phylogenetic network, multispecies coalescent simulations, a fused rooted nuclear chloroplast tree, and quantification of nuclear gene tree concordance. We show that incongruence between nuclear and plastid trees is pervasive and cannot be explained by incomplete lineage sorting alone. Instead, our results demonstrate that events consistent with repeated chloroplast capture have shaped the reticulate evolutionary history of the clade, especially among African polyploid and Pan-American diploid lineages.

Phylogeny

First identification and molecular subtyping of Blastocystis spp. in donkeys in Aksaray province, Türkiye.

Blastocystis is a common intestinal protist worldwide that can infect humans and animals. Although its molecular epidemiology in Türkiye is mostly focused primarily on humans and livestock, equids have received limited attention despite their traditional roles and frequent contact with humans and other animals in rural environments. This study aimed to determine the molecular prevalence and subtype (ST) distribution of Blastocystis spp. in donkeys in Aksaray Province, providing the first molecular data on donkeys in Türkiye. A total of 182 fresh fecal samples were collected from donkeys in nine villages within Aksaray province. Genomic DNA was extracted, and the small subunit ribosomal RNA (SSU rRNA) gene fragment of Blastocystis spp. was amplified via PCR analysis. Positive isolates were sequenced bidirectionally for identification and subsequent phylogenetic analysis of Blastocystis in donkeys. The overall molecular prevalence of Blastocystis spp. in donkeys was 4.4% (8/182). The infection rate was higher in young donkeys (under 3 years old; 8.33%) than in adults (3 years or older; 2.46%). However, this difference was not statistically significant. Sequence analysis of the positive PCR products revealed the presence of one known livestock-specific subtype, ST10. Phylogenetic analysis showed that the ST10 isolates characterized in this study clustered with isolates identified from different hosts. This study provides the first molecular data on Blastocystis presence in donkeys in Türkiye. The exclusive detection of ST10 suggests potential cross-species transmission, likely facilitated by the traditional practice of co-housing donkeys with other animals in confined barns. These findings indicate that donkeys may contribute to Blastocystis transmission, underscoring the importance of a "One Health" approach in future epidemiological surveillance.

Animals

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

Algae-to-host horizontal gene transfer in Paramecium bursaria is associated with host adaptation during endosymbiosis.

Paramecium bursaria maintains a stable endosymbiosis with green algae, yet the evolutionary consequences of this association remain unclear. Here, we screened the host genome for algal-derived horizontally transferred genes (HTGs) using a lineage-aware workflow designed to detect horizontal gene transfer (HGT) between two defined lineages. We identified 16 candidate HTGs, including four putative newly transferred genes and 12 homologous transferred genes, most of which were functionally associated with redox homeostasis and metabolism. Five HTGs showed symbiosis-dependent expression. RNAi knockdown of GH32s and SATs reduced host proliferation, total cell area, and motility, while GH32s knockdown also reduced endosymbiont load. Duplication patterns suggest that most transfers may have occurred after the P. bursaria lineage diverged from the sampled Paramecium species but before its lineage-specific whole-genome duplication (WGD). The HTGs also showed host-associated shifts in GC content and gene length, while representative HTGs retained conserved domains and functional motifs. Together, our results support algae-to-host HGT in P. bursaria and suggest that some transferred genes may contribute to metabolic integration during endosymbiosis.

Gene Transfer, Horizontal

Genetic Diversity of BK Polyomavirus Among Renal Transplant Recipients in Yunnan, China.

BK polyomavirus (BKV) infection, a common complication following kidney transplantation, can lead to BKV-associated nephropathy (BKVN). Molecular genetic studies have classified BKV into four genotypes (I-IV); however, comprehensive molecular characterization of BKV strains circulating in China remains limited. This study aimed to elucidate the predominant subtypes and clinical infection characteristics of BKV strains among kidney transplant recipients in Yunnan, a province in southwestern China. PCR-amplified BKV DNA sequences from kidney transplant recipients were aligned with reference strains and subjected to phylogenetic analysis. The viral VP1 gene was successfully amplified from 180 participants, spanning 16 ethnic groups. Genotype I was the predominant viral strain (56.66%, 102/180), followed by genotype IV (43.33%, 78/180), while genotypes II and III were not detected. Among genotypic subtypes, IVc-1 was most prevalent (40.0%, 72/180), followed by Ic (38.3%, 69/180) and Ib-1 (18.3%, 33/180). IVa-1 and IVa-2 were rare, identified in only 0.6% (n = 1) and 2.2% (n = 4) of cases, respectively. No significant differences in sex, age, BKVN incidence, BK viremia, or viruria were observed between patients with BKV-I and BKV-IV infections. Among the five confirmed BKVN cases, two were genotyped as subtype Ic, one as Ib-1, and two as IVc-1. Clinical phenotypes were also comparable between patients with BKV-I and BKV-IV infections. This study represents the largest single-center sequencing analysis of BKV in kidney transplant recipients in China, offering a valuable genomic resource for future research.

Humans

Pervasive hybridization and introgression in Diervilleae (Caprifoliaceae).

Diervilleae (Caprifoliaceae) is a horticulturally important lineage with striking floral diversity and a long history of interspecific crossing, suggesting reticulate evolution. We integrated nuclear SNPs and whole plastome data to reconstruct a phylogenomic backbone for the tribe and to identify hybrids, cultivated accessions, and introgression among lineages. Nuclear and plastid phylogenies consistently recover Weigela and Diervilla as reciprocally monophyletic and resolve four major lineages within Weigela, providing a reproducible framework for revising sectional limits and species boundaries. Cultivated accessions form a well supported clade sister to W. florida and show predominantly W. florida ancestry while retaining contributions from multiple wild lineages, consistent with recurrent crossing, backcrossing, and selection. Analyses of wild populations reveal recurrent hybrids and enable plausible parental combinations to be inferred. Tests across the genome further indicate strong evidence for historical introgression across Diervilleae, with the strongest signals involving W. middendorffiana, W. maximowiczii, and Diervilla. Fossil evidence, divergence time estimation, and paleodistribution modelling together suggest range expansion during the Miocene and Pliocene followed by climate driven contraction, providing a spatiotemporal context for episodic contact, introgression, and the East Asia-North America disjunction.

Hybridization, Genetic

Evolutionary diversification of invertase paralogs couples carbon metabolism and sexual reproduction in fission yeasts.

Dynamic patterns of gene gain and loss play a major role in the diversification of eukaryotes, reflecting adaptation to a broad range of ecological contexts. Reconstructing the evolutionary history of genes provides a powerful framework for understanding how functional innovation shapes life-history traits. Here we report a comprehensive analysis of gene gain and loss across the fission yeast clade, whose evolutionary trajectory remains elusive. Reductive evolution of metabolic genes is a major contributor to species diversification, as observed in other fungal taxa. Notably, we uncovered an evolutionary scenario in which an ancestral gene duplication was followed by lineage-specific loss of one or the other paralog, except in S. pombe, which retained both. We demonstrate that these paralogs encode catalytically-active invertases, named Inv1 and Inv2, with distinct enzymatic properties, localization, regulation, and physiological roles. Inv1 is a secreted enzyme subject to glucose catabolite repression and is the sole invertase required for sucrose assimilation, resembling canonical yeast invertases. In contrast, Inv2 is intracellular, constitutively expressed, and required for inducing sexual differentiation in response to nutrient availability. Overall, these findings reveal an unexpected role for carbon metabolism in modulating the haploid-diploid cycle of fission yeasts, suggesting that diversification of core metabolic functions may contribute to adaptation to environments with distinct sugar compositions.

Evolution

Uce-based phylogeny and classification of Megachilini.

The generic-level classification of the bee tribe Megachilini (Megachilidae) has remained controversial due to poor phylogenetic resolution at the base of the group, particularly among the brood parasitic genera and the numerous dauber ("Chalicodoma s. l.") lineages. We present a phylogenomic analysis of Megachilini based on ultraconserved elements (UCEs), sampling 52 ingroup taxa with emphasis on the dauber lineages. We also present a combined UCE + six-gene analysis to improve taxon coverage, resulting in a dataset with 127 ingroup taxa. Maximum likelihood, coalescent, and Bayesian analyses of multiple UCE matrices recover largely congruent topologies with substantially improved support relative to previous studies. Our results strongly support the monophyly of Megachilini, the early divergence of Noteriades and Gronoceras, and a single origin of brood parasitism. All remaining non-parasitic Megachilini form a moderately supported clade sister to the brood parasitic lineage. The leafcutter bees are monophyletic and nested within dauber lineages. Several major dauber clades are consistently recovered, including an exclusively Australian clade corresponding to the Hackeriapis group of subgenera, while several recognized subgenera are paraphyletic. The lineage known as Morphella, previously placed in synonymy with the subgenus Callomegachile, was not closely related to that subgenus and is here treated as a valid subgenus. Divergence-time analyses place the crown age of Megachilini in the late Eocene to early Oligocene, with major extant lineages diversifying during the Miocene. Limited morphological diagnosability of several clades indicates that splitting non-parasitic lineages into numerous genera would result in an impractical classification that would widen the gap between taxonomists and non-specialists and exacerbate the taxonomic impediment in bees. We therefore advocate retaining a single genus Megachile for non-parasitic Megachilini (excluding Noteriades and Gronoceras), as the classification best supported by phylogenomic evidence and most robust to future taxon sampling.

Animals

Evolutionary characterization and expression profiling of ACC and FASN genes in Chinese mitten crab Eriocheir sinensis.

Acetyl-CoA carboxylase (ACC) and fatty acid synthase (FASN) are rate-limiting enzymes in the fatty acid biosynthetic pathway, yet their evolutionary relationships, sequence features, and expression profiles remain poorly understood in crustaceans, particularly in the economically important Chinese mitten crab (Eriocheir sinensis). Here, we identified and systematically analyzed ACC and FASN genes in E. sinensis using comparative genomic analyses across 43 species. ACC was highly conserved as a single-copy gene in invertebrates, in contrast to the multiple paralogs observed in vertebrates. Similarly, FASN was generally maintained as a single-copy gene across most taxa but exhibited lineage-specific expansion in certain insect groups. Phylogenetic and structural analysis revealed strong conservation of both genes within crustaceans, supported by multiple conserved motifs and canonical functional domains. Expression profiling showed predominant expression in the hepatopancreas and midgut, suggesting their potential involvement in crustacean lipid metabolism. During the molting cycle, ACC and FASN exhibited higher expression levels during stages C and D, suggesting an increased capacity for fatty acid biosynthesis before molting. In addition, dietary lipid levels experiment revealed that ACC and FASN expression responded dynamically to dietary lipid availability, with increased expression at moderate lipid levels but reduced expression under excessive lipid supplementation, indicating a possible adaptive transcriptional response to lipid status. Collectively, this study provides insights into the evolutionary conservation and expression dynamics of ACC and FASN and improves our understanding of lipid metabolic adaptation in crustaceans.

Animals

Longitudinal whole-genome analysis of bluetongue virus identifies conserved serotype-specific genomes and distinct genomic constellations within a Colorado sheep flock (2021-2023).

Bluetongue virus (BTV) is a segmented double-stranded RNA virus of ruminants transmitted by Culicoides spp. biting midges. Although the genome consists of ten segments, classification into serotypes is primarily based on genome segment 2. However, reassortment among genomic segments is a major driver of BTV evolution and diversity. This study used longitudinal whole-genome sequencing to characterize BTV genomes collected from 2021 to 2023 within a single sheep flock in Colorado, where multiple serotypes co-circulate. Whole-genome sequences were generated from fourteen blood samples representing four serotypes: BTV-6, -11, -13, and -17. Longitudinal sampling identified multiple BTV serotypes within individual sheep across consecutive years. Tanglegram analysis comparing segment phylogenies to the segment 2 tree demonstrated incongruent topologies across all genomic segments, suggestive of reassortment or the circulation of distinct genomic constellations. Nucleotide-level comparisons revealed high sequence homology among same-serotype samples from the same year, while the greatest genetic divergence was observed among BTV-17 genomes collected in different years. Additionally, all BTV-13 genomes contained a previously undescribed nonsynonymous substitution in segment 10 predicted to extend the encoded protein by three amino acids. Together, these findings demonstrate that highly conserved BTV genomes and distinct genomic constellations can be detected at the flock level across multiple years. This longitudinal whole-genome approach reveals the genetic complexity of endemic BTV populations, including novel variants and genomic patterns consistent with reassortment that are lost with conventional serotyped-based approaches, highlighting the need to integrate whole-genome characterization into endemic BTV monitoring programs.

Animals

Genome-wide insights into the evolutionary and demographic history of the red alga Mazzaella laminarioides: Evidence for speciation with ancient migration along the southeast Pacific coast.

The mechanisms driving lineage divergence in red algae remain unexplored, despite the group's remarkable diversity and ancient evolutionary history. The red alga Mazzaella laminarioides, a Chilean intertidal species complex composed of three parapatric cryptic lineages (North, Center, South), offers a valuable system to evaluate these processes, as its life history combines severe dispersal limitation with a haploid-diploid cycle that may influence the emergence of reproductive barriers. We reconstructed its evolutionary history using whole-genome sequencing and nuclear genome assembly of representative individuals from each lineage. Phylogenomic analyses based on 1,507 single-copy orthologs recovered three deeply divergent lineages with limited nuclear discordance consistent with incomplete lineage sorting. For both splits, demographic modelling was most consistent with an Ancient Migration scenario, although support over strict isolation was moderate, suggesting that divergence may have begun with low asymmetric ancestral gene flow followed by subsequent loss of connectivity, demographic bottlenecks, and later population expansion. Coding sequence analyses revealed lineage-specific dN/dS heterogeneity; only one South-lineage locus passed FDR correction (metaxin-1, mitochondrial protein import), with two further South-lineage candidates in chlorophyll and heme biosynthesis falling below the FDR threshold. Together, these signals suggest that divergent selective pressures on energy acquisition may have contributed to divergence at the southern end of the distribution. These results add to the small but growing body of whole-genome data for red algae and, alongside recent macroalgal studies, suggest that ancestral connectivity could be a recurrent feature of lineage divergence even in marine organisms with extremely restricted dispersal.

Rhodophyta

Introgression shapes the genomic conflict landscape of Malus, providing evidence for a reticulate backbone in a woody crop lineage.

Phylogenomic discordance is widespread across plants, but its evolutionary significance is often obscured when conflict is treated primarily as analytical noise rather than as evidence of underlying processes. In woody lineages in particular, incomplete lineage sorting, introgression, and genome duplication can interact over long timescales to produce complex genomic histories that are not adequately summarized by a strictly bifurcating tree. Here, we use Malus as a model woody genus to investigate how these processes structure conflict across a genus-scale, accession-based phylogenomic framework. Using broad taxon sampling, hundreds of nuclear loci, plastid genomes, and genome-wide SNP summaries, we reconstruct a robust nuclear backbone for sampled Malus lineages and evaluate where discordance is concentrated and which processes best explain it. Nuclear analyses resolve eight major clades, whereas conflict is non-random and localized to recurrent hotspots rather than evenly distributed across the tree. Cytonuclear discordance is similarly concentrated, especially around Clade H, represented by sampled accessions of M. tschonoskii, where localized plastid-nuclear disagreement is consistent with candidate plastid capture or organellar introgression. Multiple complementary analyses further indicate that the strongest conflict is not explained by ILS alone, but instead reflects lineage-structured introgression, while polyploid complexes represent additional localized sources of evolutionary complexity. Together, these results provide evidence for a reticulate genomic backbone in Malus and show how integrating nuclear, plastid, and genome-wide conflict analyses can help distinguish background discordance from process-specific signals in woody plant radiations. Several lineage-level reticulation hypotheses identified here should now be tested with broader population-level sampling and curated reference accessions.

Malus

First insights into the role of evolutionary history in shaping venom composition of Vipera ammodytes.

Understanding intraspecific venom variation requires distinguishing the contributions of neutral population history from natural selection. This study aims to determine whether venom variation in Vipera ammodytes species complex is structured across eight phylogenetic groups. Despite a complex evolutionary history, venom composition did not differ among phylogenetic groups within the analytical framework used, suggesting that shared ancestry alone does not explain venom variation. Whether local adaptation to environmental conditions explains the observed variation remains an open question for future studies.

Animals

Genome-wide identification and expression profiling of HSD3B and SDR42E1 genes in the Pacific oyster (Crassostrea gigas): potential associations with gonadal development.

Sex steroids are lipid-soluble signaling molecules that regulate sex differentiation, reproductive development and physiological homeostasis in animals. 3β-Hydroxysteroid dehydrogenase/Δ5-Δ4 isomerase (3β-HSD) is a key steroidogenic enzyme, whereas SDR42E1, an extended short-chain dehydrogenase/reductase, has been implicated in sterol- and steroid-related metabolism. However, the composition, evolutionary relationships and expression patterns of the HSD3B- and SDR42E1-related genes in bivalve gonadal development remain poorly characterized. In this study, five PF01073-containing genes, comprising three CgHsd3b and two CgSdr42e1 genes, were identified in the Pacific oyster Crassostrea gigas. Phylogenetic analysis separated the proteins into HSD3B-related and SDR42E1-related groups, and gene-structure and motif analyses indicated subfamily-level divergence. All five proteins retained the SDR domain but differed in exon-intron structure and motif composition. Each contained the extended-SDR TGxxGxxG motif, whereas exact classical [ST]GxxxGxG and NNAG motifs were absent. Tyr- and Lys-equivalent residues were conserved, while the HSD3B1 Ser-equivalent position contained Thr in two C. gigas proteins and Ser in one. These features support their classification as extended-SDR proteins but do not establish enzymatic activity or substrate specificity. The three CgHsd3b genes were dispersed on one chromosome, whereas CgSdr42e1-1 and CgSdr42e1-2 were adjacent on another chromosome, suggesting a possible local duplication event for the CgSdr42e1 pair. Public RNA-seq data showed distinct tissue- and gonadal-stage expression patterns, with several genes displaying gonad-biased or female-stage-associated expression. Independent RT-qPCR profiling of the representative genes CgHsd3b-3 and CgSdr42e1-1 detected stage-dependent expression, although tissue rankings differed from those in the public RNA-seq datasets. These differences may reflect the use of independent biological samples, tissue composition, normalization procedures, and platform-specific measurements. Because enzymatic assays, metabolite measurements, cellular localization, and functional perturbation were not performed, the results identify candidate genes whose expression is associated with gonadal development rather than demonstrating regulatory roles. This study provides a comparative framework for future functional investigation of sterol- and steroid-related metabolism in bivalves.

Animals