Search PubMed⌕ Search

SEARCH · Search PubMed

Search Search PubMed

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 19 recordsLinked to original sources

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↗

Genome-wide Identification and Expression Profiling Reveal the Galectin Gene Family Diversity and their Possible Role in Antibacterial Mucosal Immunity in Japanese Flounder (Paralichthys olivaceus).

Galectins are a family of proteins that bind specifically to β-galactosides. Their importance in innate immunity of mammals has been well-documented. However, the systematic identification and characterization of galectin gene family remain limited in teleost. In this study, we identified 13 galectin genes (lgals2, lgals2a, lgals2b, lgals3, lgals3a, lgals3b, lgals4, lgals8, lgals8a, lgals9, grp, grp-b, grp-c) from Paralichthys olivaceus genome and analyzed their tissue expressions and expressions in response to Gram-negative and Gram-positive bacterial infections in mucosal tissues (gills, intestine and skin). The P. olivaceus galections were classified into three distinct types based on carbohydrate recognition domains (CRDs). Phylogenetic and syntenic analyses revealed that these galectins are closely related to their counterparts in turbot and zebrafish. Moreover, the transcripts of the 13 galectins were widespread across all tested tissues of healthy fish and regulated following challenge with Vibrio anguillarum or Streptococcus iniae in mucosal tissues, indicating their involvement in P. olivaceus immune response to bacterial infections. The lgals2a was significantly upregulated in the three mucosal tissues by either bacterial infection, whereas lgals9 and grp were basically downregulated in these tissues by either infection. On the other hand, the lgals3b and lgals4 exhibited a bacteria-specific responsive expression as they were upregulated by V. anguillarum whereas remained stable upon S. iniae infection in the gills. We also observed a positive correlation between expression level and bacterial load for the upregulated galectin genes and a negative correlation for the downregulated galectin genes. These results suggest a functional divergence among galectin members in mucosal immunity against bacterial infection in P. olivaceus.

Animals↗

Viral surveillance beyond detection: JMTV and the need for ensemble approaches in emerging virus discovery.

The recent report by T. Murillo, L. E. Enrique Chaves-González, S. Temmam, S. Bermúdez, et al. (Microbiol Spectr 14:e04078-25, 2026, https://doi.org/10.1128/spectrum.04078-25) expands the known geographic and ecological range of Jingmen tick virus (JMTV) by detecting the virus in Amblyomma mixtum ticks collected from horses in Costa Rica. This is an important finding because A. mixtum can feed on wildlife, domestic animals, and humans, creating a possible interface for virus movement across various hosts. The study also places the Costa Rican virus in a wider phylogenetic context, linking it to JMTV diversity reported from other regions. However, the detection of viral RNA in ticks should not be interpreted as proof of local disease, human infection, or active transmission, especially in the absence of supporting results. Instead, it reflects an important signal for careful viral surveillance. Here, I discuss how JMTV illustrates the need for ensemble approaches that combine field sampling, phylogeny, segment-level genome analysis, serology, experimental validation, and data-driven virus discovery tools.

emerging viruses↗

The genomic origins and evolutionary path to a key innovation in the world's most venomous snakes.

Evolutionary innovation is a catalyst for the colonization of new environments and the adaptive radiations of major groups. Novel traits typically evolve through the modification of preexisting characters, but the genetic paths underlying their origin have been challenging to trace, and the general requirements for and relative order of different kinds of gene mutations have been difficult to assess. Here, we trace the genomic origins of four procoagulant venom toxins (factor X, factor V, group I phospholipase A2, and Kunitz-type toxins) that collectively underlie a novel, especially potent blood-clotting venom type in the recently evolved Australian brown snake and taipan clade. We find evidence for a previously unknown fifth toxin, coagulation factor VII, and show that the toxins evolved through two distinct genetic paths. The factor X and factor V toxins evolved through the sequential de novo co-option of ancestral clotting factor proteins that entailed their heterotopic expression in the venom gland, the fixation of segmental duplications containing each locus, and subsequent gain-of-function mutations that rendered factor X and factor V constitutively active. In contrast, the phospholipase A2 and Kunitz-type toxins evolved by modifying the functions of neurotoxins that were part of the venom arsenal. Our findings support models in which innovative mutations in single-copy genes precede gene duplication in the evolution of novel proteins and offer a rare view into the genesis of a complex trait that has played a central role in a major adaptive radiation.

Animals↗

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↗

Genetic diversity and recombination of NA-PRRSV field strains in Vietnam: Implications for vaccine efficacy.

Porcine reproductive and respiratory syndrome (PRRS) causes severe reproductive losses in pregnant sows and piglets, resulting in substantial economic impact on the swine industry worldwide. However, due to the significant genetic diversity and rapid evolutionary changes of the pathogen, continuous surveillance and detailed genetic analysis of circulating strains are essential. The current study aimed to evaluate the genetic diversity of the hypervariable (HV) region of non-structural protein 2 (nsp2) among North American PRRSV strains isolated from swine farms in Vietnam. Phylogenetic analysis and multiple sequence alignment were conducted to determine subtype classification and assess genetic variability. A total of 48 field isolates were obtained, of which 12.5% belonged to classical NA-PRRSV, 16.6% to NADC30-like and 70.9% to HP-PRRSV, primarily distributed across sublineages 1.4, 5.1, 8.7 and 8.9. Amino acid comparisons found multiple insertions, deletions and substitutions at various positions within the hypervariable region of nsp2. The study revealed substantial genetic variation in the HV region of nsp2 among NA-PRRSV field strains, largely associated with recombination and immune escape. These findings highlight epidemiological risks to vaccine efficacy and underscore the need for continuous molecular surveillance to support effective PRRSV control in Vietnam.

PRRSV↗

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↗

Comparative genomic epidemiology of food- and patient-derived diarrheagenic Escherichia coli from sentinel surveillance in Southeast China.

Diarrheagenic Escherichia coli (DEC) remains an important foodborne pathogen, yet long-term comparative genomic surveillance data jointly characterizing food-derived and patient-derived isolates remain limited. This surveillance-based comparative study integrated antimicrobial susceptibility testing and whole-genome sequencing to characterize diarrheagenic Escherichia coli isolates recovered from food and patient sources in Lishui, Southeast China, during 2018-2025, with emphasis on occurrence, resistance profiles, genomic backgrounds, and plasmid replicon-associated features. Antimicrobial susceptibility testing was performed for 258 selected isolates, and whole-genome sequencing was conducted for a curated analytical subset of 204 isolates. The sequenced subset was used for diversity-oriented comparative genomic analysis rather than for unbiased prevalence estimation of the entire DEC collection. EAEC predominated in both sources, although food-associated occurrence was heterogeneous across categories, with the highest recovery rate observed in raw meat. Patient-derived isolates showed a broader overall resistance burden, whereas food-derived isolates retained substantial resistance to tetracycline, chloramphenicol, and florfenicol. Phylogenetic analysis showed partial overlap in genomic backgrounds between food-derived and patient-derived isolates, while representative resistance determinants displayed both broadly distributed and lineage-enriched patterns. Replicon-based plasmid profiling identified 42 plasmid types, including 12 detected in both sources, with IncF-related replicons predominating among these shared profiles. Several food-derived isolates carried multiple plasmid replicon types that were also observed in patient-derived isolates. Overall, food-derived and patient-derived DEC showed partial overlap in genomic backgrounds, resistance determinants, and replicon-defined plasmid profiles within this surveillance setting, while retaining source-associated heterogeneity. These findings should be interpreted as surveillance-based comparative evidence rather than as evidence of direct source attribution or transmission.

Humans↗

Contrasting redox-related physiological responses associated with HaGATA23 and HaGATA36 during Orobanche cumana parasitism in sunflower (Helianthus annuus L.).

Helianthus annuus is an economically important Asteraceae species used for seed oil production and ornamental purposes, but its production is seriously affected by the root-parasitic plant Orobanche cumana. GATA transcription factors are zinc-finger DNA-binding regulators involved in plant development and stress adaptation. However, the molecular characteristics of GATA transcription factors in Helianthus annuus and their contribution to Helianthus annuus -Orobanche cumana interaction remain poorly understood. Here, 36 HaGATA members were retrieved from the Helianthus annuus genome and classified into four phylogenetic clades. Chromosomal placement, collinearity, gene structure, motif composition, and promoter elements varied among the 36 HaGATA members, indicating evolutionary conservation coupled with functional diversification. Expression analysis and RT-qPCR analyses revealed differential expression patterns among HaGATA genes under O. cumana stress, with HaGATA23 markedly downregulated and HaGATA36 strongly upregulated. Overexpression of HaGATA23 was associated with increased malondialdehyde (MDA) accumulation and unfavorable changes in antioxidant enzyme activities, whereas its silencing showed the opposite physiological tendency. In contrast, overexpression of HaGATA36 reduced malondialdehyde accumulation, increased peroxidase (POD), catalase (CAT), and superoxide dismutase (SOD) activities, while its silencing showed the reverse tendency. These results indicate that HaGATA23 and HaGATA36 are candidate genes associated with contrasting redox-related physiological responses during O. cumana stress. This work provides evidence that GATA transcription factors are associated with redox-related physiological responses in sunflower under O. cumana treatment and identifies HaGATA23 and HaGATA36 as functionally divergent candidate genes for further validation.

Helianthus↗

Gloeotrichia echinulata genomes from the United States are nontoxigenic and likely geosmin producers.

Six Gloeotrichia echinulata genomes derived from planktonic harmful algal blooms (HABs) with similar colonial morphology have been sequenced from lakes in the west and northeast regions of USA, four of them to completion. The c. 7 Mbp genomes exhibit a high level of conservation, with 98-99% pairwise genome-wide average nucleotide identity and high levels of synteny, representing a single species cluster. We observed strong conservation of gene clusters responsible for the synthesis of the secondary metabolites and bioactive peptides that are characteristic of HAB-forming cyanobacteria. All six G. echinulata genomes lack genes for the synthesis of classic cyanotoxins, including microcystin, but possess genes responsible for the synthesis of the taste and odor compound geosmin. Interestingly, the geoA geosmin synthase gene in three genomes is homologous to other cyanobacterial geoA genes, while the other three geoA genes are related to actinomyces geoA. Phylogenomic analysis places the G. echinulata genomes within a clade of benthic Nostocales, reflecting an ecological niche featuring extensive growth on the sediment surface before colonies disperse into the epilimnion for planktonic growth. We identify genes conserved in all six genomes that could represent physiological adaptations supporting active growth on sediments and pelagic recruitment independent of wind-driven mixing: phycoerythrin light harvesting complexes for optimal photosynthesis at depth; gliding motility to access patchy nutrient distributions; and gas vesicles with relatively small GvpC proteins that predict resistance to higher hydrostatic pressure. The strong genomic similarity across geographically distant populations suggests that G. echinulata in the United States is a tightly related non-toxigenic species group with predictable properties relevant to public health and drinking water management.

Cyanobacteria↗

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↗

Genome-wide identification of CXE gene family in soybean and functional characterization of GmCXE31 in lipid biosynthesis and salt tolerance.

GmCXE31 negatively regulates salt tolerance and lipid synthesis in soybean, and the cxe31-edited lines improve soybean yield and seed quality. Carboxylesterases (CXEs), as essential lipid hydrolases of the α/β-hydrolase fold superfamily, are critical for plant stress responses, hormone signaling and secondary metabolism. The key candidate gene GmCXE31 was previously identified in our laboratory through a genome‑wide association study (GWAS) of soybean lipid‑related traits. In the present study, we further identified 60 GmCXE family genes in soybean. Phylogenetic analysis clustered them into 11 conserved subfamilies. Cis-acting element analysis showed their promoters are enriched with elements related to abiotic stress, growth and hormone signaling, suggesting potential roles in soybean development and stress adaptation. GmCXE31 is highly expressed in seedling roots and responsive to strigolactones (SLs) and salt stress. Functional assays revealed that GmCXE31 negatively regulates soybean salt tolerance: its overexpression reduced salt tolerance in Arabidopsis and soybean under 150 mM NaCl stress, while its knockout enhanced this trait. Lipid profiling revealed GmCXE31-edited lines had higher seed oil content, elevated oleic/linoleic acid ratio and lower saturated fatty acid proportion, which was achieved by regulating lipid synthesis-related genes like GmNFYA. Agronomic trait analysis showed GmCXE31-edited lines had increased nodule number, plant height and single-plant yield at maturity, with opposite phenotypes in overexpression lines. In conclusion, this study elucidates the multifaceted roles of GmCXE31 in coordinating soybean salt tolerance, lipid metabolism and agronomic traits, providing theoretical and genetic resources for salt-tolerant and high-quality soybean molecular breeding.

Glycine max↗

Genome-wide identification of the peanut HD-Zip gene family and AhHDZ15 positively regulating salt and drought stress in heterologously overexpressed Arabidopsis.

Homeodomain-leucine zipper (HD-Zip) transcription factors play important roles in plant growth, development, and abiotic stress responses. However, bioinformatic analyses and functional studies of HD-Zip family in peanut are scarce. In this study, 128 AhHDZ genes were identified and classified into four subfamilies in the phylogenetic analysis. Transcriptomic data and RT-qPCR analysis indicated the expression levels of AhHDZ4 and AhHDZ15 were significantly elevated in response to 12 h of salt stress, while AhHDZ4/15/60/69/126 all showed a progressive increase over time in response to drought stress. AhHDZ15 protein was localized in the nucleus. Under salt and drought stress, the germination rates of AhHDZ15-overexpressing in Arabidopsis were significantly higher than wild-type (WT), and root lengths were also significantly longer than WT. In addition, the SOD, CAT, chlorophyll content, and Relative Leaf Water Content (RLWC) value of leaves in AhHDZ15-overexpressing lines were significantly higher than WT, while the MDA content was significantly lower than WT. The above results indicate that heterologous overexpression of AhHDZ15 enhanced salt and drought tolerance in Arabidopsis. Furthermore, AhHDZ15 could bind to the L1-box element of the AhVNI2 promoter, thereby activating AhVNI2 transcription and enhancing the expression of downstream salt stress-responsive genes. These findings implies a potential function of AhHDZ15 in peanut that requires further validation.

Arabidopsis↗

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↗

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↗

Transcriptomic and RNAi analyses reveal chloride channel 3-associated osmoregulation in Litopenaeus vannamei under low-salinity stress.

Chloride channels and transporters are important for cellular volume regulation and salinity adaptation in euryhaline crustaceans, yet the intestinal transcriptional relationship between plasma-membrane and intracellular chloride pathways remains unclear in Litopenaeus vannamei. In this study, RNA interference of anoctamin 1 (ANO1) was combined with intestinal transcriptome sequencing under the production-relevant low-salinity condition of salinity 3. ANO1 silencing produced a focused transcriptional response, with 16 differentially expressed genes (DEGs) identified (11 upregulated and 5 downregulated). Functional enrichment indicated that these genes were associated with transporter activity, cytoskeletal organization, extracellular matrix-receptor interaction, membrane lipid metabolism, and vesicular processes. Notably, a transcript encoding chloride channel protein 3 (CLC-3) was significantly upregulated following ANO1 knockdown, suggesting a potential transcriptional relationship between ANO1 and CLC-3 in chloride homeostasis. Based on this finding, CLC-3 was selected for full-length cDNA cloning, sequence characterization, salinity-gradient expression analysis, and RNAi-based functional assessment. The cloned CLC-3 cDNA was 2883 bp in length and encoded an 850 amino acid protein containing a conserved voltage-gated chloride channel (Voltage-CLC) domain and two cystathionine β-synthase domains. Phylogenetic analysis placed LvCLC-3 within the intracellular CLC-c clade, and tissue distribution analysis showed the highest CLC-3 expression in the intestine. Intestinal CLC-3 expression responded nonlinearly to salinity variation, peaking at salinity 20. Under salinity 3, CLC-3 knockdown reduced ANO1, Na+/K+-ATPase alpha subunit, and Na+-K+-2Cl- cotransporter transcript levels, whereas glutamate-gated chloride channel expression increased. Mild hepatopancreatic structural alterations were also observed after CLC-3 knockdown. These findings suggest that CLC-3 is a salinity-responsive intracellular chloride-transporter candidate associated with intestinal ion-transport-related transcriptional responses after ANO1 suppression in L. vannamei, although the underlying physiological mechanism requires further validation.

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↗

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↗