Genomic Footprints of Historical Introgression Between Ancient Lineages of Wild Oryza AA-Genome Species With Widely Separated Contemporary Distributions.

Kanako O Koyanagi; Yuta Kotoku; Yuji Kishima

doi:10.1093/gbe/evag191

Home/Archive/PubMed

PubMed · 42698189

Genomic Footprints of Historical Introgression Between Ancient Lineages of Wild Oryza AA-Genome Species With Widely Separated Contemporary Distributions.

Kanako O Koyanagi·Yuta Kotoku·Yuji Kishima

View original source ↗Export metadata ↓

Abstract

Phylogenetic incongruence is increasingly recognized as pervasive, yet the extent to which reticulate evolution occurs between groups separated by substantial geographical distances and deep phylogenetic divergence remains poorly characterized. In the Oryza AA-genome group-a model for plant speciation and domestication-the traditional bifurcation model posits that Australian Oryza meridionalis and African Oryza longistaminata occupy basal branches, distinct from the more recently diversified monophyletic clade comprising Asian and other African lineages, including major cultivars. However, recent evidence from endogenous viral sequences has hinted at unexpected genetic relatedness between African O. longistaminata and Asian Oryza sativa, which are geographically and phylogenetically distant. Here, we conducted a genome-wide survey across 11 Oryza species to systematically identify genomic regions exhibiting phylogenetic incongruence. Widespread phylogenetic discordance was observed, notably involving genomic segments in which O. longistaminata showed phylogenetic proximity to Asian species, contradicting their established deep divergence. To distinguish between introgression and incomplete lineage sorting, we performed four-taxon ABBA-BABA tests, which provided statistical support for introgression. Furthermore, divergence time estimates for these incongruent regions were younger than the species divergence times, suggesting historical introgression between the ancestors of lineages that are currently separated by vast geographical distances. Systematic assessments indicated that potential analytical artifacts, such as compositional bias and substitution saturation, were unlikely to explain the observations. These convergent lines of evidence suggest that ancient introgression had occurred between currently geographically separated and evolutionarily divergent Oryza lineages, leaving detectable footprints across their modern genomes.

Explore related subjects

Oryza Genome, Plant Phylogeny Genetic Introgression Evolution, Molecular Genetic Speciation ancient introgression molecular dating phylogenetic incongruence rice Journal Article

Keep this discovery

Explore connections, maps & timelines

BibTeX RIS

Kanako O Koyanagi, Yuta Kotoku, Yuji Kishima. 2026-09-02. Genomic Footprints of Historical Introgression Between Ancient Lineages of Wild Oryza AA-Genome Species With Widely Separated Contemporary Distributions.. https://doi.org/10.1093/gbe%2Fevag191

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related citations

PubMed2026

Circadian- and light-regulated oscillatory expression of CSA in rice leaves is required for pollen fertility.

The oscillatory expression of CSA in rice leaves is regulated by the circadian clock and red/far-red light signals, mediated through DOF5 and PIL11, and is required for normal pollen fertility. Photoperiod-sensitive male-sterile lines represent a pivotal innovation in the development of hybrid rice. However, the underlying mechanisms governing photoperiod-sensitive male reproductive development remain poorly understood. Our previous studies demonstrated that the carbon starved anther (csa) mutant exhibits male sterility under short-day (SD) conditions but partial fertility under long-day (LD) conditions. In this study, we report that CSA expression follows an oscillatory rhythm in rice leaves under both SD and LD conditions, a pattern regulated by both circadian clock and light signals. Tissue-specific RNA interference knockdown of CSA in leaves was associated with reduced pollen viability, suggesting that CSA expression in leaves contributes to normal male fertility. Promoter truncation assay results indicate that distinct regions of the CSA promoter contribute differentially to the regulation of CSA expression in leaves versus anthers, and that both the CSA expression level in anthers and the rhythmic expression pattern of CSA in leaves are associated with the restoration of male fertility. Using dual-luciferase, yeast one-hybrid, and electrophoretic mobility shift assays, we identified two proteins, PIL11 and DOF5, which directly bind to specific motifs (an E-box and T/AAAAG motif) within the CSA promoter truncation, thereby regulating its transcription. These findings elucidate novel mechanisms linking light sensing to the expression of circadian-controlled genes, thus connecting photoperiod with male reproductive development in rice.

Oryza↗

PubMed2026

Quantitative trait locus qLDC5 regulates primary root branching in an auxin-dependent manner.

L-type lateral root (LLR) density determines root system architecture, affecting nutrient acquisition in rice (Oryza sativa L.), particularly under low-phosphorus conditions. Previous studies identified genotypic differences in LLR density and a quantitative trait locus (QTL) enhancing LLR density on crown roots (qLDC5). We showed that LLR densities on crown and primary roots were closely correlated and confirmed higher LLR density on primary roots in qLDC5 donor DJ123 compared with the African variety NERICA4 using X-ray micro-computed tomography. We confirmed the qLDC5 effect in a field experiment for LLR density on primary roots. LLR densities on primary and crown roots, therefore, appear under similar genetic control. Developmental analyses revealed that DJ123 and NDJ188-a derivative line harboring qLDC5-initiate more lateral root primordia than NERICA4, with a higher proportion progressing to elongation, but that exogenous auxin application reversed this ranking. Within qLDC5, auxin biosynthesis gene OsYUCCA2 and auxin response factor OsARF15 were up-regulated in DJ123. Transcriptome analysis revealed an indirect auxin-mediated regulatory network underlying LLR variation. Differentially expressed genes in DJ123 and NDJ188 were enriched for ent-kaurene and gibberellin metabolism, including the robust induction of OsGA2ox5. These findings suggest qLDC5 increases lateral root density by coordinating gibberellin, auxin, and terpene pathways.

Oryza↗

PubMed2026

Circular RNAs orchestrate integrated post-transcriptional responses to combined heat and drought stress in rice.

Circular RNAs (circRNAs) are emerging post-transcriptional regulators, yet their landscape and functional roles in rice under combined abiotic stress remain largely unexplored. Here, we systematically reanalyzed strand-specific RNA-seq data to characterize circRNAs responsive to simultaneous heat and drought stress. Following quality control, read mapping, and dual-algorithm prediction using CIRI2 and CIRCexplorer2, we identified 208 high-confidence circRNAs distributed across all 12 chromosomes. Comparative profiling revealed 83 circRNAs uniquely expressed in control samples, 51 in stressed samples, and 74 shared between conditions, indicating stress-dependent circularization. Junction-read analysis highlighted a spectrum of circularization strength, ranging from highly abundant circRNAs with dominant junction reads to low-confidence candidates masked by linear transcript background. Genomic annotation showed that circRNAs primarily originated from exonic and intergenic regions, with a pronounced negative-strand bias; several genes generated multiple circRNA isoforms via alternative back-splicing. Functional enrichment of host genes suggested involvement in protein folding, nutrient reservoir activity, RNA degradation, and branched-chain amino acid catabolism, implicating roles in stress adaptation and metabolic regulation. Differential expression analysis identified seven circRNAs specifically induced under combined stress conditions. Network topology analysis pinpointed key miRNAs-including osa-miR414, osa-miR1439, and osa-miR2919-as candidate topological hubs within the predicted network. Their predicted target genes, such as those encoding stress-responsive transcription factors and signaling proteins, suggest potential roles in coordinating post-transcriptional responses to combined stress. Network topology analysis pinpointed key miRNAs-including osa-miR414, osa-miR1439, and osa-miR2919-as candidate topological hubs within the predicted network. Their predicted target genes, such as those encoding stress-responsive transcription factors and signaling proteins, suggest potential roles in coordinating post-transcriptional responses to combined stress. Overall, this study provides a comprehensive map of circRNAs in rice under combined heat and drought stress, suggests their potential as ceRNAs based on predictive analysis, and lays a foundation for future experimental validation of circRNA-mediated regulation.

Oryza↗