PubMed HealthSearch

SEARCH · PubMed Health

Results for “evolutionary analysis”

Explore indexed PubMed citations for clinical trials, systematic reviews and public health research. Read source abstracts and follow each citation to its original 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

Comprehensive genomic and computational insights into Brucella suis: pan-genome analysis, evolutionary perspectives, and in-silico vaccine design.

BACKGROUND: Brucella suis is a zoonotic intracellular pathogen responsible for brucellosis, mainly in swine and humans. Although numerous genome sequences are publicly available, an integrative genomic analysis combining pan-genome architecture, structural organization, evolutionary relationships, and vaccine-associated targets remains limited. RESULTS: In this study, we analyzed 91 publicly available B.suis genomes to characterize their pan-genome composition and genomic structure. The pan-genome exhibited an open configuration, indicating continued genomic diversification. A total of 2,146 core genes were identified, representing conserved functions essential for species maintenance, while the accessory genome reflected strain-level variability. Phylogenetic reconstruction based on single-copy orthologs revealed distinct evolutionary clades among the strains. A complementary phylogenetic analysis of pan-genome gene presence-absence patterns further supported clade differentiation and highlighted variation in accessory gene repertoires. Comparative synteny and genome structural analyses demonstrated largely conserved chromosomal organization with localized rearrangements across strains. Screening of the core proteome identified 64 putative antigenic proteins with predicted surface localization and immunogenic properties. Additionally, resistance-associated determinants related to tetracycline and doxycycline were detected in one genome within the dataset. CONCLUSIONS: This comprehensive genomic analysis defines the pan-genome structure, evolutionary relationships, and genome organization of B.suis. The integration of core and pan-genome-based phylogenies provides complementary insights into strain diversification, while the identified conserved antigenic candidates offer a foundation for future experimental validation and rational vaccine development strategies.

Genome, Bacterial

Genome-wide identification, characterization, evolutionary analysis, and expression profiling of the FCS-like zinc finger (FLZ) gene family in soybean (Glycine max L.) under abiotic stresses.

Drought and salinity limit soybean yield. Despite their role in the SnRK1 energy-sensing complex, a systematic study of FCS-Like Zinc Finger (FLZ) proteins in soybean has not been reported. We performed a genome-wide identification of the GmFLZ gene family, identifying 40 members distributed across 18 of the 20 soybean chromosomes. Phylogenetic analysis of 87 FLZ proteins from Glycine max, Arabidopsis thaliana, and Oryza sativa revealed four major evolutionary clades, suggesting that diversification predates the separation of monocots and dicots. Structural analysis identified ten conserved motifs, with Motifs 1 and 2 present in all family members. Gene duplication analysis identified 304 paralogous pairs, most arising from segmental duplication. Ka/Ks analysis indicated localized positive selection in six gene pairs and purifying selection in 97.9% of pairs. Tissue-specific expression profiling across nine tissues showed that GmFLZ5, GmFLZ15, GmFLZ25, and GmFLZ34 had the highest expression levels detected across the GmFLZ family, with GmFLZ5 the most highly expressed member in leaves, nodules, and stem and showing moderate expression in pod, root, and root hairs, whereas GmFLZ18, GmFLZ23, and GmFLZ37 showed root-preferential expression. RT-qPCR validation under drought (20% PEG-6000) and salt (200 mM NaCl) treatments in the Giza 5 cultivar showed that 36 and 34 of the 40 GmFLZ genes, respectively, exhibited at least a two-fold change in expression, with GmFLZ21 and GmFLZ35 among the most strongly induced under salt stress. These findings provide an evolutionary and functional framework for the GmFLZ family and identify candidate genes for future functional studies in soybean stress tolerance.

Glycine max

Comprehensive identification and evolutionary analysis of the Wnt gene family in bivalves: Insights into the larval development of the noble scallop Chlamys nobilis.

The Wnt gene family regulates fundamental developmental processes in metazoans, but its evolutionary composition and developmental deployment in bivalves remain largely unresolved. Here, we performed a comparative genomic analysis of Wnt genes in 19 bivalve species and examined developmental expression profiles in the noble scallop Chlamys nobilis, with Crassostrea gigas and Chlamys farreri used for cross-species comparison. A total of 235 Wnt genes were identified and assigned to 12 subfamilies. No reliable Wnt3 ortholog was detected in any analyzed bivalve, supporting the view that Wnt3 loss occurred early during lophotrochozoan evolution rather than representing a lineage-specific absence. Most Wnt proteins retained the conserved WNT domain, indicating strong structural conservation, whereas lineage-specific copy-number variation and gene loss were observed among species. C. farreri and C. gigas each retained 12 Wnt genes and lacked Wnt3, whereas C. nobilis lacked Wnt3, Wnt7, and Wnt16. Developmental transcriptome analysis and RT-qPCR revealed clear stage-specific expression patterns. In C. gigas, Wnt2/10/A were highly expressed during earlydevelopment and peaked around the D-shaped larval stage, while Wnt8 and Wnt11 showed distinct stage-specific peaks. By contrast, Wnt1/5/6/9 were more active during later larval development or juvenile formation. These results provide a comparative framework for bivalve Wnt evolution and identify candidate Wnt genes potentially involved in larval development and aquaculture-relevant developmental transitions.

Animals

Genome-wide identification and evolutionary analysis of the ERF-VII gene family in the tea plant (Camellia sinensis) and functional characterization of CsRAP2.2 in response to cold stress.

The ERF-VII gene family, a critical branch of the AP2/ERF superfamily, is central to plant stress adaptation. However, its evolutionary history and function in tea plant (Camellia sinensis) remain unclear. Here, we performed integrated evolutionary, genomic, and functional analyses of ERF-VII genes across 14 plant lineages and 20 tea plant cultivars. The phylogenetic analysis revealed that ERF-VII proteins originated after vascular plant divergence, coinciding with the emergence of the N-terminal MCGGA/I motif linked to the oxygen-dependent N-degron pathway. Gymnosperms retained few conserved members, whereas angiosperms exhibited lineage-specific expansion-extensive in monocots via whole-genome duplication, moderate in eudicots with functional diversification. Pan-genome analysis across 20 tea plant cultivars further revealed varietal differences in ERF-VII gene distribution. Transcriptome profiling via the Tea Plant Information Archive identified CsRAP2.2 as a cold-inducible ERF-VII member with sustained expression under low-temperature stress. Functional assays demonstrated that silencing CsRAP2.2 reduced cold tolerance, while overexpression in tea leaves and heterologous expression in Arabidopsis thaliana enhanced cold tolerance by maintaining photosystem II efficiency, reducing membrane lipid peroxidation, and improving antioxidant capacity. Weighted gene co-expression network analysis positioned CsRAP2.2 as a regulatory hub integrating cold, hormone, and oxygen-sensing pathways. These results clarify the evolutionary trajectory of ERF-VII genes and establish CsRAP2.2 as a core cold-tolerance regulator in tea plant. These findings may inform future breeding of cold-resilient tea cultivars.

Camellia sinensis

Genomic and evolutionary analysis reveals dynamic variations of MKK3 gene, a key regulator for seed dormancy in barley.

Barley (Hordeum vulgare L.) is an important crop in the world, and its seed dormancy is primarily controlled by a mitogen-activated protein kinase kinase 3 (MKK3) gene. Although kinase activity of MKK3 and its roles in barley post-domestication have been widely studied, the pre-domestication evolution of MKK3 and the spread of nondormant alleles among global barley varieties remain largely unexplored. In this study, we analyzed MKK3 sequences in barley and its wild progenitor (Hordeum spontaneum K. Koch) and identified two polymorphic miniature inverted-repeat transposable elements (MITEs). Comparative analyses indicated that the insertions/excision of the MITEs predated the current estimates of barley domestication. Examination of the barley pangenomes coupled with droplet digital polymerase chain reaction revealed extensive copy number variation of MKK3 and suggested that transposons likely contributed to tandem amplification of the MKK3 gene on chromosome 5H. Additionally, approximately 1-Kb MKK3 sequences were found on chromosomes 1H and 6H. Further analysis indicated that these short MKK3 sequences were captured by a CACTA transposon that also contained fragments from four other expressed genes. The acquisition of MKK3 was estimated to be between 1.9 and 2.5 million years ago. Together, these findings illuminate the dynamic pre-domestication evolution of the MKK3 gene and identify three divergent MKK3 haplotype groups including a unique lineage predominant in Ethiopian germplasm. This study highlights the contribution of transposons to structural diversification and evolutionary differentiation of the MKK3 locus and provides helpful information for understanding the complex history of MKK3 gene in barley and also for improving preharvest sprouting tolerant varieties under distinct natural conditions.

Hordeum

Genome-wide identification, structural characterization, and evolutionary analysis of growth-related gene families in African catfish (Clarias gariepinus).

The somatotropic axis encompassing growth hormone (GH), insulin-like growth factor (IGF), myostatin (MSTN), and prolactin (PRL) signalling cascades is the master regulator of somatic growth, metabolism, and development in vertebrates. African catfish (Clarias gariepinus), a commercially pivotal aquaculture species, now possesses a chromosome-level reference genome (CGAR_prim_01v2); however, a systematic, genome-wide characterization spanning all five interconnected growth-related gene families has not previously been undertaken in this species. Here, we identified and characterized 15 growth-related genes spanning gh1, ghra, ghrb, Igf1, Igf2a, Igf2b, igf1ra, Igf1rb, Igf2r, Mstna, Mstnb, prl, prlra, prlrb, and smtlb distributed across 13 chromosomes. Complete one-to-one orthology with zebrafish confirmed strong dosage-balance conservation across >120 million years of teleost divergence. Physicochemical analysis resolved a clear biochemical dichotomy between compact, basic secreted ligands (19.88-45.81 kDa; pI up to 10.02) and large, acidic, heavily glycosylated membrane receptors (56.82-270.80 kDa; pI 4.85-5.97). Phylogenetic analysis confirmed 3R whole-genome duplication origins for all paralog pairs, while synteny analysis revealed a disruption of the ancestral gh1-prl chromosomal block in C. gariepinus, a finding that warrants further comparative and functional investigation. This genomic atlas provides the sequence and structural information including exon-intron boundaries, domain architecture, and chromosomal coordinates needed as a prerequisite for future marker-assisted selection and CRISPR-based myostatin-editing efforts in African catfish aquaculture, though translation into applied breeding outcomes will require subsequent functional and expression studies.

Animals

Meta-evolutionary exome analysis identifies novel type 2 diabetes mellitus genes in the UK Biobank and all of us.

Type 2 diabetes mellitus (T2DM) risk is heavily influenced by genetics, yet current association tests have explained only parts of its heritability. We developed MEVA (Meta-Evolutionary Action), a meta-analytic framework that integrates three complementary methods-EAML, Sigma-Diff, and GeneEMBED-to assess the functional burden of protein-coding variants using evolutionary data. MEVA was applied to exome data from 28,115 T2DM cases and 28,115 controls in the UK Biobank (UKB), identifying 101 genes (p&#x2009;<&#x2009;1e-5). MEVA outperformed its component methods, each of which substantially outperformed a conventional burden test (MAGMA), in recovering known T2DM genes (AUROC&#x2009;=&#x2009;0.925) and maintaining robustness in progressively smaller cohorts (AUROC&#x2009;=&#x2009;0.917). MEVA showed significant enrichment for T2DM-related loci (p&#x2009;=&#x2009;6.8e-10, p&#x2009;=&#x2009;2.0e-34), protein interactions (z&#x2009;=&#x2009;4.6, z&#x2009;=&#x2009;4.2), pathways (p&#x2009;=&#x2009;1.3e-6, z&#x2009;=&#x2009;2.0), phenotypes (p&#x2009;=&#x2009;1.3e-21, z&#x2009;=&#x2009;9.1), and literature mentions (z&#x2009;=&#x2009;7.2). Replication in 16,915 T2DM cases and 16,915 controls from All of Us (AoU) yielded 99 genes (p&#x2009;<&#x2009;1e-5), 23 of which were also recovered in the UKB cohort - far exceeding random chance. These included established genes (SLC30A8, WFS1, HNF1A) and less-characterized candidates (NRIP1, ADAM30, CALCOCO2, TUBB1, ZFP36L2, WDR90). Notably, NRIP1 loss-of-function variants were associated with increased T2DM risk in both the UKB (OR = 1.09, FDR&#x2009;=&#x2009;5.4e-4) and AoU (OR = 1.09, FDR&#x2009;=&#x2009;0.046), and TUBB1 and CALCOCO2 gain-of-function variants showed consistent risk effects (FDR&#x2009;<&#x2009;0.05). Pathway analyses revealed convergence on endoplasmic reticulum chaperone complexes (FDR&#x2009;=&#x2009;0.02) and Hippo signaling (FDR&#x2009;=&#x2009;8.5e-4). Finally, all 177 candidate genes were functionally prioritized using ten orthogonal criteria to guide experimental follow-up. These results demonstrate that combining complementary, impact-aware association tests increases sensitivity, improves replication, and expands the catalog of genetic risk factors for T2DM.

Humans

Genome-wide identification of WOX transcription factors and functional characterization of WOX4 and WOX13 involved in cold stress response in Malus baccata.

INTRODUCTION: Cold stress is a major abiotic threat to apple production. Malus baccata has exceptional cold hardiness and is widely used as a superior cold-resistant rootstock. The WUSCHEL-related homeobox (WOX) transcription factor family regulates plant growth, development and stress adaptation, whereas the functions of WOX genes in cold tolerance of M. baccata remain elusive. METHODS: In the present work, 19 MbWOX family members were identified and characterized at the genome-wide level. Evolutionary analysis, cis-element prediction, transcriptome profiling and real-time quantitative PCR (RT-qPCR) were performed to screen core cold-responsive genes. Overexpression vectors were constructed and transformed into Arabidopsis seedlings for functional verification. RESULTS: Evolutionary analysis revealed that segmental duplication drove the expansion of the MbWOX family, and these genes contained a variety of stress-responsive cis-elements. Combined transcriptome and RT-qPCR analyses confirmed that MbWOX4 and MbWOX13 were core cold-responsive genes with distinct expression patterns. The two genes participated in cold signal transduction by interacting with different transcription factor networks. Functional tests revealed that MbWOX4 and MbWOX13 isoforms differentially modulated seedling cold tolerance under low-temperature stress.

Malus baccata

Natural variation in the cytokinin oxidase gene ZmCKX6 influences leaf morphology and yield-related traits in maize.

Leaf width (LW) is a critical determinant of maize architecture and yield. To uncover its genetic basis, we performed a genome-wide association study (GWAS) on 348 maize inbred lines and identified ZmCKX6, encoding cytokinin oxidase/dehydrogenase, as a key gene associated with LW. Natural variation in the ZmCKX6 promoter significantly influenced its expression levels, leading to differences in LW across various haplotypes. Functional validation using CRISPR/Cas9 revealed that ZmCKX6 knockout results in pleiotropic effects, including narrower leaves, reduced plant height, and decreased grain yield components. These phenotypes were accompanied by elevated levels of active cytokinins but reduced levels of auxin, gibberellins, and salicylic acid. Transcriptome analysis revealed a significant downregulation of photosynthesis-related genes, corresponding to reduced photosynthetic rates in knockout lines. Evolutionary analysis demonstrated that the allele associated with narrower leaves were preferentially selected during maize domestication and breeding. This study highlights the role of ZmCKX6 in modulating cytokinin homeostasis and its subsequent impact on multiple agronomic traits in maize, providing insights into the complex genetic control of plant architecture and yield. The identified natural variations could be valuable for marker-assisted selection aimed at optimizing plant architecture and improving yield.

Zea mays

A high-quality draft genome assembly of Johnsongrass illuminates relationships between polyploidization, crop-wild hybridization, and reproductive biology.

Johnsongrass [Sorghum halepense (L.) Pers.] is an allopolyploid, rhizomatous, perennial grass species and one of the most troublesome weeds in global agriculture. We assembled the first Johnsongrass genome to clarify poorly understood genetic factors influencing variable rates of crop-wild hybridization with cultivated sorghum [S. bicolor (L.) Moench]. The draft genome assembly has a total size of 3.26 Gb and BUSCO completeness of 95.3%. We also report the first evolutionary analysis of INHIBITION OF ALIEN POLLEN (IAP), the only known cross-(in)compatibility locus in the genus. Our results reveal an evolutionary history of genome instability, including the loss of distinct parental subgenomes, and suggest that Nebraska accession 'J-37,' the genome donor, is a segmental allotetraploid that may function as a diploid or aneuploid during meiosis. Genome instability could explain observations of variable ploidies in Johnsongrass and facilitate ongoing hybridization with sorghum where gamete ploidies and IAP alleles match. Given this information, we provide a suggested research framework for studying evolution and gene expression in the Sorghum genus where crop-wild hybridization occurs and for predicting the potential for hybridization between specific crossing partners. Collectively, this work will bolster efforts to study and manage reproductive biology in other crop-wild polyploid complexes.

Sorghum

Expanding the human proteome with microproteins and peptideins.

A major scientific drive is to characterize the protein-coding genome, which is a primary basis for studying human health. But the fundamental question remains of what has been missed in previous analyses. Over the past decade, the translation of non-canonical open reading frames (ncORFs) has been observed across human cell types and disease states1-3, with major implications for biomedical science. However, a key gap in knowledge has been which ncORFs produce small microproteins or alternative protein molecules that contribute to the human proteome. Here we report the collaborative efforts of the TransCODE Consortium4 to produce a consensus landscape of protein-level evidence for ncORFs. We show that about 25% of a set of 7,264 ncORFs gives rise to detectable peptides in a large-scale analysis of 95,520 proteomics experiments. We develop an annotation framework for ncORF-encoded microproteins as human proteins and codify the new conceptual model of 'peptideins' as microproteins that have indeterminate potential as functional proteins. To probe the biological implications of peptideins, we create an evolutionary analysis approach, termed ORF relative branch length (ORBL), and determine that evolutionary constraint is common and associates with observation of ncORF-derived peptides. We then characterize a pan-essential cellular phenotype for one peptidein from the OLMALINC long non-coding RNA. Overall, we generate public research tools supported by GENCODE and PeptideAtlas and advance biomedical discovery for understudied components of the human proteome.

Humans

Comparative genomic analysis of Artemisia argyi reveals asymmetric expansion of terpene synthases and conservation of artemisinin biosynthesis.

Artemisia argyi, a perennial herb of the Asteraceae family, possesses significant therapeutic and economic value. We present a 7.88&#x2009;Gb chromosome-level haplotype-resolved genome assembly, revealing its unique evolutionary trajectory. The karyotype (2n&#x2009;=&#x2009;34) of A. argyi is that of an autotetraploid, which underwent gametic chromosome fusion prior to species-specific whole-genome duplication (WGD-3). The genome exhibits pronounced multivalent chromosome pairing and frequent recombination among homologous groups. Asymmetrical evolution following WGD-3 is a hallmark feature, evidenced by imbalanced allelic gene loss and widespread neofunctionalization. The terpene synthase (TPS) gene family exemplifies this pattern, having expanded through four duplication events in A. argyi. Recent tandem duplications and allelic functional differentiation have generated substantial gene functional diversity. Notably, we identified a tandem-duplicated six-copy ADS homolog (AarADS)-a key TPS gene in the artemisinin biosynthetic pathway of Artemisia annua (AanADS)-localized exclusively to a single chromosome in A. argyi. Unlike AanADS, which converts farnesyl pyrophosphate (FPP) to amorpha-4,11-diene, AarADS catalyzes FPP to &#x3b1;-bisabolol. Evolutionary analysis suggested that AanADS acquired its specialized function via a derived mutation in the A. annua lineage. This study elucidates the genomic evolution underpinning A. argyi's distinctive medicinal properties.

Alkyl and Aryl Transferases

Understanding Mycobacterium tuberculosis through its genomic diversity and evolution.

Pathogen evolution and genomic diversity are shaped by specific host immune pressures and therapeutic interventions. Analysis of the extant genomes of circulating strains of Mycobacterium tuberculosis, a leading cause of infectious mortality that has co-evolved with humans for thousands of years, can provide new insights into host-pathogen interactions that underlie specific aspects of pathogenesis and onward transmission. With the explosion in the number of fully sequenced M. tuberculosis strains that are now paired with detailed clinical data, there are new opportunities to understand the evolutionary basis for and consequences of M. tuberculosis strain diversity. This review examines mechanistic findings that have emerged from pairing whole genome sequencing data and evolutionary analysis with functional dissection of specific bacterial variants. These include improved understanding of secreted effectors that modulate the properties and migratory behavior of infected macrophages as well as bacterial genetic alterations important for survival within hypoxic microenvironments. Genomic, evolutionary, and functional analyses across diverse M. tuberculosis strains will identify prominent bacterial adaptations to their human hosts and shape our understanding of TB disease biology and the host immune response.

Mycobacterium tuberculosis

Fast, accurate construction of multiple sequence alignments from protein language embeddings.

Multiple sequence alignment (MSA) is a foundational task in computational biology, underpinning protein structure prediction, evolutionary analysis, and domain annotation. Traditional MSA algorithms rely on pairwise amino acid substitution matrices derived from conserved protein families. While effective for aligning closely related sequences, these scoring schemes struggle in the low-identity "twilight zone." Here, we present a new approach for constructing MSAs leveraging amino acid embeddings generated by protein language models (PLMs), which capture rich evolutionary and contextual information from massive and diverse sequence datasets. We introduce a windowed reciprocal-weighted embedding similarity metric that is surprisingly effective in identifying corresponding amino acids across sequences. Building on this metric, we develop ARIES (Alignment via RecIprocal Embedding Similarity), an algorithm that constructs a PLM-generated template embedding and aligns each sequence to this template via dynamic time warping in order to build a global MSA. Across diverse benchmark datasets, ARIES achieves higher accuracies than existing state-of-the-art approaches, especially in low-identity regimes where traditional methods degrade, while scaling almost linearly with the number of sequences to be aligned. Together, these results provide the first large-scale demonstration of the power of PLMs for accurate and scalable MSA construction across protein families of varying sizes and levels of similarity, highlighting the potential of PLMs to transform comparative sequence analysis.

Deep Learning

Dynamic and non-additive gene regulation shapes maize responses to simultaneous salt and cold stress.

Salt and cold stresses often occur together in nature and severely impact crop productivity, yet their transcriptional regulation remains poorly understood. Here, we conducted a time-series transcriptomic analysis of maize under salt, cold, and their combination at 0, 6, 12, and 24&#xa0;h. Differential expression analysis revealed dynamic, condition-specific gene responses grouped into eight distinct temporal patterns. Promoter motif analysis of genes within each pattern identified 5-39 significantly enriched motifs, with over 40% lacking known counterparts, suggesting the involvement of previously uncharacterized cis-regulatory elements in stress-responsive transcriptional regulation. By comparing combined stress responses to the sum of single-stress effects, we found that about 74% of DEGs showed non-additive patterns, suggesting that combined stress triggers a distinct transcriptional program. Evolutionary analysis showed that additive DEGs tend to be more recently evolved, subject to weaker purifying selection, and enriched in transposed duplications, contrasting with the stronger constraint observed in non-additive DEGs. WGCNA identified 24 co-expression modules, among which 65 hub DEGs were detected in modules significantly correlated with specific stress conditions. Furthermore, we reconstructed 228, 20, and 200 sequential transcription factor cascades spanning 6&#xa0;h, 12&#xa0;h, and 24&#xa0;h under cold, salt, and combined stress, respectively, with no cascade shared across all three conditions. Together, these results reveal that maize responses to combined salt and cold stress are largely non-additive and temporally dynamic, with distinct evolutionary patterns underlying different response types, offering insights and candidate regulators for enhancing crop stress resilience.

Zea mays

A single-nucleus transcriptome atlas of soybean anthers.

Anther development is crucial for plant sexual reproduction. However, a high-resolution, cell-type-specific transcriptomic atlas of this process is lacking for the legume crop soybean (Glycine max). Here, we construct a comprehensive transcriptional atlas of developing soybean anthers using single-nucleus RNA sequencing (snRNA-seq). We identify and characterize nine distinct cell types spanning both somatic and reproductive lineages. Our analysis reveals robust transcriptional continuity across anther developmental stages and dynamic reprogramming during key transitions. Notably, the shift from diploid meiocytes to haploid unicellular microspores is marked by the induction of previously inactive genes, despite an overall reduction in transcript abundance. Subsequently, within bicellular microspores, generative and vegetative cell lineages exhibit sharply divergent transcriptional programs: generative cells specialize in mRNA export and turnover, whereas vegetative cells up-regulate translational machinery. Evolutionary analysis further indicates that generative-cell-specific genes are subject to more relaxed purifying selection compared to those specific to vegetative cells. Functional validation using mutants generated by CRISPR/Cas9-mediated genome editing and EMS mutagenesis reveals the essential roles of OSD1A and PKSA in pollen development and fertility. This high-resolution atlas provides fundamental insights into the transcriptional regulation of soybean anther development and serves as a valuable resource for manipulating male fertility to advance hybrid breeding programs. The data are available at https://databases.genedenovo.com/pollen.

Glycine max

Identification of rice DUF1719 gene family and analysis of alkaline tolerance function of OsDUF1719.8.

Alkaline stress severely constrains the physiological metabolism and growth and development of rice through high pH and ionic toxicity. Domains of unknown function (DUF) play significant roles in plant stress responses. However, the function of the DUF1719 family (PF08224) in rice has not been reported and further research is needed. This study systematically identified the OsDUF1719 gene family in rice and investigated the function of OsDUF1719.8 under alkaline stress. The results demonstrate that the rice DUF1719 family comprises 13 protein members, all containing the PF08224 domain. It is predicted that this domain may play a role in ATPase activation. Evolutionary analysis divided DUF1719 proteins from eight grass species into six subgroups, with highly conserved gene structures, motifs, and tertiary architectures within each subgroup. Promoter analysis indicated enrichment of stress- and hormone-responsive elements, implying broad involvement in stress regulation. Expression analysis revealed that several genes, including OsDUF1719.5 and OsDUF1719.8, were upregulated under multiple abiotic stresses. Notably, OsDUF1719.8 was strongly induced during early alkaline stress. Consequently, we further analyzed the function of OsDUF1719.8 in the rice alkaline stress response. The results demonstrate that overexpression of OsDUF1719.8 enhanced rice alkaline tolerance, whereas knockout mutants exhibited stress sensitivity. OsDUF1719.8 enhances rice tolerance to alkaline stress by coordinately regulating reactive oxygen species metabolism, promoting the accumulation of osmotic adjustment compounds, and modulating ion homeostasis. This study provides the first systematic identification of the DUF1719 family and elucidates the function of OsDUF1719.8 in positively regulating rice alkaline tolerance, offering a novel gene for alkali-tolerant molecular breeding of rice.

Oryza

Comparison of phylogenetic metrics of transmission in symptomatic and asymptomatic tuberculosis.

BACKGROUND: Understanding drivers of Mycobacterium tuberculosis (Mtb) transmission remains a critical challenge in high-burden settings. Tuberculosis control efforts traditionally target symptomatic individuals, yet the role of asymptomatic cases in sustaining transmission is increasing recognized. METHODS: We conducted a genomic and epidemiological analysis of Mtb isolates collected in Mato Grosso do Sul, Brazil, between 2008 and 2024. From 2017 to 2022, active case finding was performed in three of the state's largest prisons, whereby sputum was collected from individuals irrespective of symptoms and tested by GeneXpert and culture. We evaluated several metrics of recent transmission from symptomatic and asymptomatic individuals, including phylogenetic clustering, Time-scaled Haplotype Density (THD), Local Branching Index (LBI), and transmission probabilities inferred using the Bayesian Reconstruction and Evolutionary Analysis of Transmission Histories (BREATH). FINDINGS: We sequenced 2,362 Mtb strains, of which 3.5% (115/2,362) were resistant to at least one drug, and 0.6% (16/2,362) were multi-drug resistant. Most strains were lineage 4, and 78.2% of all isolates were part of a genomic cluster. Among 2,362 individuals with tuberculosis, 1,137 were incarcerated at the time of diagnosis. Among these, 505 were identified through active case finding: 277 had symptomatic disease and 228 had asymptomatic tuberculosis. There was no significant difference in phylogenetic clustering proportion (77% vs. 85%; p= 0.816), THD (median 0.50 vs. 0.39; p = 0.120), or LBI (median 0.00863 vs. 0.00871; p = 0.086) between symptomatic and asymptomatic individuals. Bayesian transmission trees revealed no significant difference in the number of secondary infections inferred from symptomatic compared with asymptomatic individuals (p = 0.56). These findings were consistent across genomic clusters and robust to model assumptions. INTERPRETATION: We identified no differences in transmission from symptomatic compared with asymptomatic individuals, using several genomic measures of transmission, underscoring the substantial contribution that asymptomatic tuberculosis makes to transmission at the population level.

Asymptomatic