Ecological divergence in the presence of gene flow in two closely related Oryza species (Oryza rufipogon and O. nivara)

Wild rice: unlocking the future of rice breeding

Germplasm resources serve as the foundations of advancements in breeding and are crucial for maintaining food security. Wild rice species of the genus Oryza include rich sources of genetic diversity and high adaptability, making them a substantial resource for rice breeding. The discovery of wild-type cytoplasmic male sterility resources enabled the achievement of the 'three lines' goal in hybrid rice, significantly increasing rice yields. The application of resistance alleles from wild rice enables rice production to withstand losses caused by stress. Reduced genetic diversity due to rice breeding poses a significant limitation to further advances and can be alleviated through a systematic use of wild genetic resources that integrate geographic, climatic and environmental data of the original habitat, along with extensive germplasm collection and identification using advanced methods. Leveraging technological advancements in plant genomics, the understanding of genetic mechanisms and the application of artificial intelligence and gene editing can further enhance the efficiency and accuracy of this process. These advancements facilitate rapid isolation and functional studies of genes, and precise genome manipulation. This review systematically summarizes the utilization of superior genes and germplasm resources derived from wild rice sources, while also exploring the collection, conservation, identification and utilization of further wild rice germplasm resources. A focus on genome sequencing and biotechnology developments is leading to new breeding and biotechnology opportunities. These new opportunities will not only promote the development of rice varieties that exhibit high yields, superior stress resistance and high quality but also expand the genetic diversity among rice cultivars.

Genetic architecture of ecological divergence between Oryza rufipogon and Oryza nivara

Ecological divergence due to habitat difference plays a prominent role in the formation of new species, but the genetic architecture during ecological speciation and the mechanism underlying phenotypic divergence remain less understood. Two wild ancestors of rice (Oryza rufipogon and Oryza nivara) are a progenitor-derivative species pair with ecological divergence and provide a unique system for studying ecological adaptation/speciation. Here, we constructed a high-resolution linkage map and conducted a quantitative trait locus (QTL) analysis of 19 phenotypic traits using an F2 population generated from a cross between the two Oryza species. We identified 113 QTLs associated with interspecific divergence of 16 quantitative traits, with effect sizes ranging from 1.61% to 34.1% in terms of the percentage of variation explained (PVE). The distribution of effect sizes of QTLs followed a negative exponential, suggesting that a few genes of large effect and many genes of small effect were responsible for the phenotypic divergence. We observed 18 clusters of QTLs (QTL hotspots) on 11 chromosomes, significantly more than that expected by chance, demonstrating the importance of coinheritance of loci/genes in ecological adaptation/speciation. Analysis of effect direction and v-test statistics revealed that interspecific differentiation of most traits was driven by divergent natural selection, supporting the argument that ecological adaptation/speciation would proceed rapidly under coordinated selection on multiple traits. Our findings provide new insights into the understanding of genetic architecture of ecological adaptation and speciation in plants and help effective manipulation of specific genes or gene cluster in rice breeding.

Multiple domestications of Asian rice

The origin of domesticated Asian rice (Oryza sativa L.) has been controversial for more than half a century. The debates have focused on two leading hypotheses: a single domestication event in China or multiple domestication events in geographically separate areas. These two hypotheses differ in their predicted history of genes/alleles selected during domestication. Here we amassed a dataset of 1,578 resequenced genomes, including an expanded sample of wild rice from throughout its geographic range. We identified 993 selected genes that generated phylogenetic trees on which japonica and indica formed a monophyletic group, suggesting that the domestication alleles of these genes originated only once in either japonica or indica. Importantly, the domestication alleles of most selected genes (~80%) stemmed from wild rice in China, but the domestication alleles of a substantial minority of selected genes (~20%) originated from wild rice in South and Southeast Asia, demonstrating separate domestication events of Asian rice.

Gene Genealogy-Based Mutation Analysis Reveals Emergence of Aus, Tropical japonica, and Aromatic of Oryza sativa during the Later Stage of Rice Domestication

Asian rice (Oryza sativa L.) has become a model for understanding gene functions and domestication in recent decades; however, its own diversification is still controversial. Although the division of indica and japonica and five subgroups (aus, indica (sensu stricto), japonica (sensu stricto), tropical japonica, and aromatic) are broadly accepted, how they are phylogenetically related is not transparent. To clarify their relationships, a sample of 121 diverse genes was chosen here from 12 Oryza genomes (two parental and ten O. sativa (Os)) in parallel to allow gene genealogy-based mutation (GGM) analysis. From the sample, 361 Os mutations were shared by two or more subgroups (referred to here as trans mutations) from 549 mutations identified at 51 Os loci. The GGM analysis and related tests indicates that aus diverged from indica at a time significantly earlier than when tropical japonica split from japonica. The results also indicate that aromatic was selected from hybrid progeny of aus and tropical japonica and that all five subgroups share a significant number of the early mutations identified previously. The results suggest that aus, tropical japonica, and aromatic emerged sequentially within the most recent 4-5 millennia of rice domestication after the split of indica and japonica.

Genetic diversity and population dynamic of Ziziphus jujuba var. spinosa (Bunge) Hu ex H. F. Chow in Central China

Phylogeographic research concerning Central China has been rarely conducted. Population genetic and phylogeography of Ziziphus jujuba var. spinosa (also called sour jujube) were investigated to improve our understanding of plant phylogeographic patterns in Central China. Single-copy nuclear gene markers and complete chloroplast genome data were applied to 328 individuals collected from 21 natural populations of sour jujube in China. Nucleotide variation of sour jujube was relatively high (π = 0.00720, θ _w = 0.00925), which resulted from the mating system and complex population dynamics. Analysis of molecular variation analysis revealed that most of the total variation was attributed to variation within populations, and a high level of genetic differentiation among populations was detected (F _st = 0.197). Relatively low long-distance dispersal capability and vitality of pollen contributed to high genetic differentiation among populations. Differences in the environmental conditions and long distance among populations further restricted gene flow. Structure clustering analysis uncovered intraspecific divergence between central and marginal populations. Migrate analysis found a high level of gene flow between these two intraspecific groups. Bayesian skyline plot detected population expansion of these two intraspecific groups. Network and phylogeny analysis of chloroplast haplotypes also found intraspecific divergence, and the divergence time was estimated to occur at about 55.86 Ma. Haplotype native to the Loess Plateau was more ancient, and multiple glacial refugia of sour jujube were found to locate at the Loess Plateau, areas adjacent to the Qinling Mountains and Tianmu Mountains. Species distribution model analysis found a typical contraction-expansion model corresponding to the Quaternary climatic oscillations. In the future, the distribution of sour jujube may shift to high-latitude areas. This study provides new insights for phylogeographic research of temperate plant species distributed in Central China and sets a solid foundation for the application of the scientific management strategy of Z. jujuba var. spinosa.

Phenotypic Variation and the Impact of Admixture in the Oryza rufipogon Species Complex (ORSC)

Crop wild relatives represent valuable reservoirs of variation for breeding, but their populations are threatened in natural habitats, are sparsely represented in genebanks, and most are poorly characterized. The focus of this study is the Oryza rufipogon species complex (ORSC), wild progenitor of Asian rice (Oryza sativa L.). The ORSC comprises perennial, annual and intermediate forms which were historically designated as O. rufipogon, O. nivara, and O. sativa f. spontanea (or Oryza spp., an annual form of mixed O. rufipogon/O. nivara and O. sativa ancestry), respectively, based on non-standardized morphological, geographical, and/or ecologically-based species definitions and boundaries. Here, a collection of 240 diverse ORSC accessions, characterized by genotyping-by-sequencing (113,739 SNPs), was phenotyped for 44 traits associated with plant, panicle, and seed morphology in the screenhouse at the International Rice Research Institute, Philippines. These traits included heritable phenotypes often recorded as characterization data by genebanks. Over 100 of these ORSC accessions were also phenotyped in the greenhouse for 18 traits in Stuttgart, Arkansas, and 16 traits in Ithaca, New York, United States. We implemented a Bayesian Gaussian mixture model to infer accession groups from a subset of these phenotypic data and ascertained three phenotype-based group assignments. We used concordance between the genotypic subpopulations and these phenotype-based groups to identify a suite of phenotypic traits that could reliably differentiate the ORSC populations, whether measured in tropical or temperate regions. The traits provide insight into plant morphology, life history (perenniality versus annuality) and mating habit (self- versus cross-pollinated), and are largely consistent with genebank species designations. One phenotypic group contains predominantly O. rufipogon accessions characterized as perennial and largely out-crossing and one contains predominantly O. nivara accessions characterized as annual and largely inbreeding. From these groups, 42 "core" O. rufipogon and 25 "core" O. nivara accessions were identified for domestication studies. The third group, comprising 20% of our collection, has the most accessions identified as Oryza spp. (51.2%) and levels of O. sativa admixture accounting for more than 50% of the genome. This third group is potentially useful as a "pre-breeding" pool for breeders attempting to incorporate novel variation into elite breeding lines.

A snapshot of progenitor-derivative speciation in Iberodes (Boraginaceae)

Traditional classification of speciation modes has focused on physical barriers to gene flow. Allopatric speciation with complete reproductive isolation is viewed as the most common mechanism of speciation. Parapatry and sympatry, by contrast, entail speciation in the face of ongoing gene flow, making them more difficult to detect. The genus Iberodes (Boraginaceae, NW Europe) comprises five species with contrasting morphological traits, habitats, and species distributions. Based on the predominance of narrow and geographically distant endemic species, we hypothesized that geographic barriers were responsible for most speciation events in Iberodes. We undertook an integrative study including: (i) phylogenomics through restriction-site associated DNA sequencing, (ii) genetic structure analyses, (iii) demographic modeling, (iv) morphometrics, and (v) climatic niche modeling and niche overlap analysis. Results revealed a history of recurrent progenitor-derivative speciation manifested by a paraphyletic pattern of nested species differentiation. Budding speciation mediated by ecological differentiation is suggested for the coastal lineage, deriving from the inland widespread I. linifolia during Late Pliocene. Meanwhile, geographic isolation followed by niche shifts are suggested for the more recent differentiation of the coastland taxa. Our work provides a model for distinguishing speciation via ecological differentiation of peripheral, narrowly endemic I. kuzinskyanae and I. littoralis from a widespread extant ancestor, I. linifolia. Ultimately, our results illustrate a case of Pliocene speciation in the probable absence of geographic barriers and get away from the traditional cladistic perspective of speciation as producing two species from an extinct ancestor, thus reminding us that phylogenetic trees tell only part of the story.

Study of Rhizosphere Microbial Community Structures of Asian Wild and Cultivated Rice Showed That Cultivated Rice Had Decreased and Enriched Some Functional Microorganisms in the Process of Domestication

Asian cultivated rice (Oryza sativa L.), domesticated from Asian wild rice, is a staple food crop for populations around the world. Asian cultivated rice has undergone physiological changes in the process of its evolution from Asian wild rice, and the closely related rhizosphere microorganisms may have changed in the process of plant domestication. However, the rhizosphere microorganisms of different Asian wild rice species and their related indica and japonica cultivated rice have not yet been illustrated clearly. This study aimed to illustrate the microbial community structures in the rhizosphere of Asian wild rice (common wild rice, nivara wild rice, medicinal wild rice, and spotted wild rice) and Asian cultivated rice (indica and japonica accessions) through the high-throughput sequencing of 16S rDNA, ITS amplifiers and metagenomic data. The results showed that there were significant differences between wild and cultivated rice in their rhizosphere microbial community structures. In view of the indica and japonica rice, the bacterial and fungal community structures of indica rice with the nivara wild rice and medicinal wild rice were more similar than the japonica rice species. The indica and japonica rice had the lowest proportion of Actinobacteria than the wild rice species, and indica rice has the highest relative abundance of Nitrospira. As for the microbial functions, methane metabolism and pyruvate metabolism were found to be the common pathway enriched in the rhizosphere of common and nivara wild rice in comparison with the indica and japonica rice; in addition, though it was found that the relative abundances of the pathogenic fungi in the rhizosphere soil of indica and japonica rice were significantly lower than that of the wild rice, the relative abundances of Magnaporthales and Ustilaginales were significantly higher in indica and japonica rice than that of the wild rice. This study is expected to provide a theoretical basis for the development and utilization of rhizosphere microbial resources for wild and cultivated rice.

Species diversification in the lowlands of mid-latitude South America: Turnera sidoides subsp. carnea as a case study

The lowlands of mid-latitude South America comprise complex temperate ecoregions characterized by a unique biodiversity. However, the processes responsible for shaping its species diversity are still largely unknown. Turnera sidoides subsp. carnea is a variable subspecies occurring in the lowlands of northeastern Argentina and Uruguay, extending to southern Paraguay and Brazil. It constitutes a good model to perform evolutionary studies. Here we used an integrative approach to understand the process of diversification within this subspecies and to increase the knowledge concerning patterns and processes responsible for shaping the species diversity in the temperate lowlands of South America. The results provided strong evidences that this subspecies is an autopolyploid complex per se, being in an active process of intrasubspecific diversification. Morphological and genetic data show that the diversity of T. sidoides subsp. carnea is in congruence with the great past and present abiotic and biotic variability of the mid-latitude South American lowlands. The evolutionary history of this subspecies is consistent with past fragmentation and allopatric differentiation at diploid level. Geographic isolation and local adaptation would have promoted strong morphological, ecological, and genetic differentiation, resulting in two morphotypes and different genetic groups indicative of incipient speciation.

Rapid adaptive radiation in a hillstream cyprinid fish in the East African White Nile River basin

Adaptive radiation of freshwater fishes was long thought to be possible only in lacustrine environments. Recently, several studies have shown that riverine and stream environments also provide the ecological opportunity for adaptive radiation. In this study, we report on a riverine adaptive radiation of six ecomorphs of cyprinid hillstream fishes of the genus Garra in a river located in the Ethiopian Highlands in East Africa. Garra are predominantly highly specialized algae-scrapers with a wide distribution ranging from Southeast Asia to West Africa. However, adaptive phenotypic diversification in mouth type, sucking disc morphology, gut length and body shape have probably been found among these ecomorphs in a single Ethiopian river. Moreover, we found two novel phenotypes of Garra ("thick-lipped" and "predatory") that had not been discovered before in this species-rich genus (>160 species). Mitochondrial and genome-wide data suggest monophyletic, intrabasin evolution of Garra phenotypic diversity with signatures of gene flow from other local populations. Although sympatric ecomorphs are genetically distinct and can be considered to being young species as suggested by genome-wide single nucleotide polymorphism data, mitochondrial DNA was unable to identify any genetic structure suggesting recent and rapid speciation events. Some data suggest a hybrid origin of the novel "thick-lipped" ecomorph. Here we highlight how, driven by ecological opportunity, an ancestral trophically highly specialized lineage is likely to have rapidly radiated in a riverine environment promoted by the evolution of novel feeding strategies.

Boundaries and hybridization in a secondary contact zone between freshwater mussel species (Family:Unionidae)

Correct species identification and delineation are crucial for effective conservation and management. However, species delineation can be problematic in the presence of morphological ambiguities due to phenotypic plasticity, convergence, and/or interspecific hybridization. Here, we investigated the degree of hybridization between two closely related freshwater mussel species [Bivalvia: Unionidae; Lampsilis siliquoidea (Barnes) and L. radiata (Gmelin)] that present intermediate forms in areas of sympatry. Unionids have a distinct form of mitochondrial DNA (mtDNA) inheritance, termed doubly uniparental inheritance (DUI) where female mtDNA (F-type) is transmitted to all progeny but male mtDNA (M-type) is mostly inherited by the males resulting in mostly homoplasmic females and heteroplasmic males. An individual was identified as hybrid when F-type and M-type mtDNA of the two different species were found in the same individual. Twelve out of 116 sequenced males were identified as hybrids indicating that these species hybridize where their geographic range overlaps in the lower Great Lakes and St. Lawrence basins. Microsatellite analyses further support the occurrence of hybridization but at a larger spatial scale than indicated by the mitochondrial analyses. We also found that strong within-species population genetic structure affects the detection of purebred individuals overestimating the number of hybrids. Given the large geographic scale and proportion of hybrids found in this study, natural hybridization and introgression need to be considered when implementing local biodiversity inventories, identifying waterbodies as source of organisms for relocation and restoration projects and when setting appropriate conservation policies.

Genome-Wide Analysis of CCT Transcript Factors to Identify Genes Contributing to Photoperiodic Flowering in Oryza rufipogon

Photoperiod sensitivity is a dominant determinant for the phase transition in cereal crops. CCT (CONSTANS, CO-like, and TOC1) transcription factors (TFs) are involved in many physiological functions including the regulation of the photoperiodic flowering. However, the functional roles of CCT TFs have not been elucidated in the wild progenitors of crops. In this study, we identified 41 CCT TFs, including 19 CMF, 17 COL, and five PRR TFs in Oryza rufipogon, the presumed wild ancestor of Asian cultivated rice. There are thirty-eight orthologous CCT genes in Oryza sativa, of which ten pairs of duplicated CCT TFs are shared with O. rufipogon. We investigated daily expression patterns, showing that 36 OrCCT genes exhibited circadian rhythmic expression. A total of thirteen OrCCT genes were identified as putative flowering suppressors in O. rufipogon based on rhythmic and developmental expression patterns and transgenic phenotypes. We propose that OrCCT08, OrCCT24, and OrCCT26 are the strong functional alleles of rice DTH2, Ghd7, and OsPRR37, respectively. The SD treatment at 80 DAG stimulated flowering of the LD-grown O. rufipogon plants. Our results further showed that the nine OrCCT genes were significantly downregulated under the treatment. Our findings would provide valuable information for the construction of photoperiodic flowering regulatory network and functional characterization of the CCT TFs in both O. rufipogon and O. sativa.

Divergence in flowering time is a major component contributing to reproductive isolation between two wild rice species (Oryza rufipogon and O. nivara)

It is of critical importance for our understanding of speciation process to determine the forms of reproductive isolation and their relative importance in species divergence. Oryza nivara and O. rufipogon are direct ancestors of Asian cultivated rice and a progenitor-daughter species pair. Investigating the reproductive isolation between them provides insights into plant speciation and helps understanding of the rice domestication. Here, we quantitatively measured the major components of reproductive isolation between the two species based on common garden and crossing experiments for three pairs of sympatric populations in Nepal, Cambodia and Laos. We revealed significant differences in the flowering times between species pairs, with O. nivara flowering much earlier than O. rufipogon. A very weak reduction in seed set but no reduction in F1 viability and fertility were detected for the crosses between species relative to those within species. Moreover, we detected asymmetrical compatibility between species and found that emasculation significantly decreased pollination success in O. nivara but not in O. rufipogon. Our study demonstrates that the divergence between O. nivara and O. rufipogon is maintained almost entirely by the difference in flowering times and suggests that differential flowering times contribute to both habitat preferences and reproductive isolation between species.

Conservation recommendations for Oryza rufipogon Griff. in China based on genetic diversity analysis

Over the past 30 years, human disturbance and habitat fragmentation have severely endangered the survival of common wild rice (Oryza rufipogon Griff.) in China. A better understanding of the genetic structure of O. rufipogon populations will therefore be useful for the development of conservation strategies. We examined the diversity and genetic structure of natural O. rufipogon populations at the national, provincial, and local levels using simple sequence repeat (SSR) markers. Twenty representative populations from sites across China showed high levels of genetic variability, and approximately 44% of the total genetic variation was among populations. At the local level, we studied fourteen populations in Guangxi Province and four populations in Jiangxi Province. Populations from similar ecosystems showed less genetic differentiation, and local environmental conditions rather than geographic distance appeared to have influenced gene flow during population genetic evolution. We identified a triangular area, including northern Hainan, southern Guangdong, and southwestern Guangxi, as the genetic diversity center of O. rufipogon in China, and we proposed that this area should be given priority during the development of ex situ and in situ conservation strategies. Populations from less common ecosystem types should also be given priority for in situ conservation.

Discrepancies between genetic and ecological divergence patterns suggest a complex biogeographic history in a Neotropical genus

Phylogenetic patterns and the underlying speciation processes can be deduced from morphological, functional, and ecological patterns of species similarity and divergence. In some cases, though, species retain multiple similarities and remain almost indistinguishable; in other cases, evolutionary convergence can make such patterns misleading; very often in such cases, the "true" picture only emerges from carefully built molecular phylogenies, which may come with major surprises. In addition, closely related species may experience gene flow after divergence, thus potentially blurring species delimitation. By means of advanced inferential methods, we studied molecular divergence between species of the Virola genus (Myristicaceae): widespread Virola michelii and recently described, endemic V. kwatae, using widespread V. surinamensis as a more distantly related outgroup with different ecology and morphology-although with overlapping range. Contrary to expectations, we found that the latter, and not V. michelii, was sister to V. kwatae. Therefore, V. kwatae probably diverged from V. surinamensis through a recent morphological and ecological shift, which brought it close to distantly related V. michelii. Through the modeling of the divergence process, we inferred that gene flow between V. surinamensis and V. kwatae stopped soon after their divergence and resumed later, in a classical secondary contact event which did not erase their ecological and morphological differences. While we cannot exclude that initial divergence occurred in allopatry, current species distribution and the absence of geographical barriers make complete isolation during speciation unlikely. We tentatively conclude that (a) it is possible that divergence occurred in allopatry/parapatry and (b) secondary contact did not suppress divergence.

Novel allelic variation in the Phospholipase D alpha1 gene (OsPLDα1) of wild Oryza species implies to its low expression in rice bran

Rice bran, a by-product after milling, is a rich source of phytonutrients like oryzanols, tocopherols, tocotrienols, phytosterols, and dietary fibers. Moreover, exceptional properties of the rice bran oil make it unparalleled to other vegetable oils. However, a lipolytic enzyme Phospholipase D alpha1 (OsPLDα1) causes rancidity and ‘stale flavor’ in the oil, and thus limits the rice bran usage for human consumption. To improve the rice bran quality, sequence based allele mining at OsPLDα1 locus (3.6 Kb) was performed across 48 accessions representing 11 wild Oryza species, 8 accessions of African cultivated rice, and 7 Oryza sativa cultivars. From comparative sequence analysis, 216 SNPs and 30 InDels were detected at the OsPLDα1 locus. Phylogenetic analysis revealed 20 OsPLDα1 cDNA variants which further translated into 12 protein variants. The O. officinalis protein variant, when compared to Nipponbare, showed maximum variability comprising 22 amino acid substitutions and absence of two peptides and two β-sheets. Further, expression profiling indicated significant differences in transcript abundance within as well as between the OsPLDα1 variants. Also, a new OsPLDα1 transcript variant having third exon missing in it, Os01t0172400-06, has been revealed. An O. officinalis accession (IRGC101152) had lowest gene expression which suggests the presence of novel allele, named as OsPLDα1-1a (GenBank accession no. MF966931). The identified novel allele could be further deployed in the breeding programs to overcome rice bran rancidity in elite cultivars.

Potential of Platinum Standard Reference Genomes to Exploit Natural Variation in the Wild Relatives of Rice

As the world's population expands from 7.6 billion to 10 billion over the next 30 years, scientists and farmers across the globe must explore every angle necessary to provide a safe, stable and sustainable food supply for generations to come. Rice, and its wild relatives in the genus Oryza, will play a significant role in helping to solve this 10 billion people question due to its place as a staple food for billions. The genus Oryza is composed of 27 species that span 15 million years of evolutionary diversification and have been shown to contain a plethora of untapped adaptive traits, e.g., biotic and abiotic resistances, which can be used to improve cultivated rice. Such traits can be introduced into cultivated rice, in some cases by conventional crossing, and others via genetic transformation and gene editing methods. In cases where traits are too complex to easily transfer to cultivated rice [e.g., quantitative trait loci (QTL)], an alternative strategy is to domesticate the wild relative that already contains the desired adaptive traits - i.e., "neodomestication". To utilize the Oryza genus for crop improvement and neodomestication, we first need a set of genomic resources that can be used to efficiently identify, capture, and guide molecular crop improvement. Here, we introduce the concept of platinum standard reference genome sequences (PSRefSeq) - a new standard by which contiguous near-gap free reference genomes can now be produced. By having a set of PSRefSeqs for every Oryza species we set a new bar for how crop wild relatives can be integrated into crop improvement programs.

A Combination of Long-Day Suppressor Genes Contributes to the Northward Expansion of Rice

Growing cultivated rice with a moderate heading date is the key to expanding its cultivation area and maintaining stable yields. The genes that regulate heading date are largely cloned; however, it remains unclear how genetic mutations and their combinations affect the heading date and adaptability of cultivated rice. Here, we report the analysis of genetic variation in eight long-day flowering suppressor genes (Hd1, DTH8, Ghd7, OsCOL4, DTH7, Hd6, Se5, and PhyB) and the phylogenetic relationship of eight genes. Genetic variations in DTH8, Ghd7, Hd1, DTH7, PhyB, and OsCOL4 are correlated with differences in heading date and the correlation between the genetic diversity of Hd6 and Se5 and rice heading data are weak. One group of haplotypes of DTH8, Ghd7, Hd1, DTH7, PhyB, and OsCOL4 are associated with earlier heading dates and appear to have accumulated during the northward expansion of rice cultivation. A minimum of four group A alleles of DTH8, Ghd7, Hd1, DTH7, PhyB, and OsCOL4 are required for the growth of cultivated rice at latitudes above 30°N. This study presents a preliminary investigation of the genetic patterns and adaptation mechanisms of long-day flowering suppressor genes and provides a useful reference for the molecular breeding of rice cultivars for various environments and farming systems.

Genome-wide analyses reveal the role of noncoding variation in complex traits during rice domestication

Genomes carry millions of noncoding variants, and identifying the tiny fraction with functional consequences is a major challenge for genomics. We assessed the role of selection on long noncoding RNAs (lncRNAs) for domestication-related changes in rice grains. Among 3363 lncRNA transcripts identified in early developing panicles, 95% of those with differential expression (329 lncRNAs) between Oryza sativa ssp. japonica and wild rice were significantly down-regulated in the domestication event. Joint genome and transcriptome analyses reveal that directional selection on lncRNAs altered the expression of energy metabolism genes during domestication. Transgenic experiments and population analyses with three focal lncRNAs illustrate that selection on these loci led to increased starch content and grain weight. Together, our findings indicate that genome-wide selection for lncRNA down-regulation was an important mechanism for the emergence of rice domestication traits.

Parallel Speciation of Wild Rice Associated with Habitat Shifts

The occurrence of parallel speciation strongly implies the action of natural selection. However, it is unclear how general a phenomena parallel speciation is since it was only shown in a small number of animal species. In particular, the adaptive process and mechanisms underlying the process of parallel speciation remain elusive. Here, we used an integrative approach incorporating population genomics, common garden, and crossing experiments to investigate parallel speciation of the wild rice species Oryza nivara from O. rufipogon. We demonstrated that O. nivara originated multiple times from different O. rufipogon populations and revealed that different O. nivara populations have evolved similar phenotypes under divergent selection, a reflection of recurrent local adaptation of ancient O. rufipogon populations to dry habitats. Almost completed premating isolation was detected between O. nivara and O. rufipogon in the absence of any postmating barriers between and within these species. These results suggest that flowering time is a “magic” trait that contributes to both local adaptation and reproductive isolation in the origin of wild rice species. Our study thus demonstrates a convincing case of parallel ecological speciation as a consequence of adaptation to new environments.

An overview on reproductive isolation in Oryza sativa complex

Reproductive isolation is generally regarded as the essence of the speciation process. Studying closely related species is convenient for understanding the genetic basis of reproductive isolation. Therefore, the present review is restricted to the species and subspecies of the Oryza sativa complex, which includes the two domestic rice cultivars and six wild species. Although closely related, these rice species are separated from each other by a range reproductive barriers. This review presents a comprehensive understanding of the forces that shaped the formation of reproductive barriers among and between the species of the O. sativa complex. We suggest the possibility that domestication and artificial breeding in these rice species can lead to the early stages of speciation. Understanding the evolutionary and molecular mechanisms underlying reproductive isolation in rice will increase our knowledge in speciation and would also offer practical significance for the implementation of crop improvement strategies.

Population genetics and evolutionary history of the wild rice species Oryza rufipogon and O. nivara in Sri Lanka

Genetic diversity and population genetic structure of the wild rice species Oryza rufipogon and O. nivara in Sri Lanka were studied using 33 microsatellite markers. A total of 315 individuals of 11 natural populations collected from the wet, intermediate, and dry zones of the country were used in the study. We found a moderate to high level of genetic diversity at the population level, with the polymorphic loci (P) ranging from 60.6% to 100% (average 81.8%) and the expected heterozygosity (H E) varying from 0.294 to 0.481 (average 0.369). A significant genetic differentiation between species and strong genetic structure within species were also observed. Based on species distribution modeling, we detected the dynamics of the preferred habitats for the two species in Sri Lanka and demonstrated that both O. rufipogon and O. nivara populations have expanded substantially since the last internal glacial. In addition, we showed that the geographical distribution of the two species corresponded to the climate zones and identified a few of key environmental variables that contribute to the distribution of the two species, implying the potential mechanism for ecological adaptation of these two species in Sri Lanka. These studies provided important insights into the population genetics and evolution of these wild species in Sri Lanka and are of great significance to the in situ conservation and utilization of these wild resources in genetic improvement of rice.

Ancient Endogenous Pararetroviruses in Oryza Genomes Provide Insights into the Heterogeneity of Viral Gene Macroevolution

Endogenous viral sequences in eukaryotic genomes, such as those derived from plant pararetroviruses (PRVs), can serve as genomic fossils to study viral macroevolution. Many aspects of viral evolutionary rates are heterogeneous, including substitution rate differences between genes. However, the evolutionary dynamics of this viral gene rate heterogeneity (GRH) have been rarely examined. Characterizing such GRH may help to elucidate viral adaptive evolution. In this study, based on robust phylogenetic analysis, we determined an ancient endogenous PRV group in Oryza genomes in the range of being 2.41-15.00 Myr old. We subsequently used this ancient endogenous PRV group and three younger groups to estimate the GRH of PRVs. Long-term substitution rates for the most conserved gene and a divergent gene were 2.69 × 10-8 to 8.07 × 10-8 and 4.72 × 10-8 to 1.42 × 10-7 substitutions/site/year, respectively. On the basis of a direct comparison, a long-term GRH of 1.83-fold was identified between these two genes, which is unexpectedly low and lower than the short-term GRH (>3.40-fold) of PRVs calculated using published data. The lower long-term GRH of PRVs was due to the slightly faster rate decay of divergent genes than of conserved genes during evolution. To the best of our knowledge, we quantified for the first time the long-term GRH of viral genes using paleovirological analyses, and proposed that the GRH of PRVs might be heterogeneous on time scales (time-dependent GRH). Our findings provide special insights into viral gene macroevolution and should encourage a more detailed examination of the viral GRH.

Morphological and molecular dissection of wild rices from eastern India suggests distinct speciation between O. rufipogon and O. nivara populations

The inter relationships between the two progenitors is interesting as both wild relatives are known to be the great untapped gene reservoirs. The debate continues on granting a separate species status to Oryza nivara. The present study was conducted on populations of Oryza rufipogon and Oryza nivara from Eastern India employing morphological and molecular characteristics. The cluster analysis of the data on morphological traits could clearly classify the two wild forms into two separate discrete groups without any overlaps i.e. lack of intermediate forms, suggesting the non-sympatric existence of the wild forms. Amplification of hyper variable regions of the genome could reveal 144 alleles suggesting high genetic diversity values (average He = 0.566). Moreover, with 42.37% of uncommon alleles between the two wild relatives, the molecular variance analysis (AMOVA) could detect only 21% of total variation (p < 0.001) among them and rest 59% was within them. The population structure analysis clearly classified these two wild populations into two distinct sub-populations (K = 2) without any overlaps i.e. lack of intermediate forms, suggesting the non-sympatric existence of the wild forms. Clear differentiation into two distinct groups indicates that O. rufipogon and O. nivara could be treated as two different species.

Natural Variation in the Promoter of GSE5 Contributes to Grain Size Diversity in Rice

The utilization of natural genetic variation greatly contributes to improvement of important agronomic traits in crops. Understanding the genetic basis for natural variation of grain size can help breeders develop high-yield rice varieties. In this study, we identify a previously unrecognized gene, named GSE5, in the qSW5/GW5 locus controlling rice grain size by combining the genome-wide association study with functional analyses. GSE5 encodes a plasma membrane-associated protein with IQ domains, which interacts with the rice calmodulin protein, OsCaM1-1. We found that loss of GSE5 function caused wide and heavy grains, while overexpression of GSE5 resulted in narrow grains. We showed that GSE5 regulates grain size predominantly by influencing cell proliferation in spikelet hulls. Three major haplotypes of GSE5 (GSE5, GSE5^DEL1+IN1, and GSE5^DEL2) in cultivated rice were identified based on the deletion/insertion type in its promoter region. We demonstrated that a 950-bp deletion (DEL1) in indica varieties carrying the GSE5^DEL1+IN1 haplotype and a 1212-bp deletion (DEL2) in japonica varieties carrying the GSE5^DEL2 haplotype associated with decreased expression of GSE5, resulting in wide grains. Further analyses indicate that wild rice accessions contain all three haplotypes of GSE5, suggesting that the GSE5 haplotypes present in cultivated rice are likely to have originated from different wild rice accessions during rice domestication. Taken together, our results indicate that the previously unrecognized GSE5 gene in the qSW5/GW5 locus, which is widely utilized by rice breeders, controls grain size, and reveal that natural variation in the promoter region of GSE5 contributes to grain size diversity in rice.

Genetic structure and isolation by altitude in rice landraces of Yunnan, China revealed by nucleotide and microsatellite marker polymorphisms

Rice landraces, a genetic reservoir for varietal improvement, are developed by farmers through artificial selection during the long-term domestication process. To efficiently conserve, manage, and use such germplasm resources, an understanding of the genetic structure and differentiation of local rice landraces is required. In this study, we analyzed 188 accessions of rice landraces collected from localities across an altitudinal gradient from 425 to 2, 274 m above sea level in Yunnan Province, China using ten target genes and 48 SSR markers. We detected clear differentiation of the rice landraces into indica and japonica groups and further separation of the accessions in each group into two subgroups according to altitude, including a lower altitude subgroup and higher altitude subgroup. The AMOVA results showed significant genetic differentiation among altitude zones at SSRs and most genes, except Os1977 and STS22. We further determined that differentiation among landrace populations followed a model of isolation by altitude, in which gene flow was higher among populations at similar altitude levels than across different altitude levels. Our findings demonstrated that both adaptation to altitude and altitude-dependent gene flow played key roles in the genetic differentiation of rice landraces in Yunnan, China.

Nonfunctional alleles of long-day suppressor genes independently regulate flowering time

Due to the remarkable adaptability to various environments, rice varieties with diverse flowering times have been domesticated or improved from Oryza rufipogon. Detailed knowledge of the genetic factors controlling flowering time will facilitate understanding the adaptation mechanism in cultivated rice and enable breeders to design appropriate genotypes for distinct preferences. In this study, four genes (Hd1, DTH8, Ghd7 and OsPRR37) in a rice long-day suppression pathway were collected and sequenced in 154, 74, 69 and 62 varieties of cultivated rice (Oryza sativa) respectively. Under long-day conditions, varieties with nonfunctional alleles flowered significantly earlier than those with functional alleles. However, the four genes have different genetic effects in the regulation of flowering time: Hd1 and OsPRR37 are major genes that generally regulate rice flowering time for all varieties, while DTH8 and Ghd7 only regulate regional rice varieties. Geographic analysis and network studies suggested that the nonfunctional alleles of these suppression loci with regional adaptability were derived recently and independently. Alleles with regional adaptability should be taken into consideration for genetic improvement. The rich genetic variations in these four genes, which adapt rice to different environments, provide the flexibility needed for breeding rice varieties with diverse flowering times.

Diachronic analysis of genetic diversity in rice landraces under on-farm conservation in Yunnan, China

Diachronic analysis showed no significant changes in the level of genetic diversity occurred over the past 27 years' domestication, which indicated genetic diversity was successfully maintained under on-farm conservation. Rice (Oryza sativa L.) is one of the earliest domesticated crop species. Its genetic diversity has been declining as a result of natural and artificial selection. In this study, we performed the first analysis of the levels and patterns of nucleotide variation in rice genomes under on-farm conservation in Yunnan during a 27-year period of domestication. We performed large-scale sequencing of 600 rice accessions with high diversity, which were collected in 1980 and 2007, using ten unlinked nuclear loci. Diachronic analysis showed no significant changes in the level of genetic diversity occurring over the past 27 years' domestication, which indicated genetic diversity was successfully maintained under on-farm conservation. Population structure revealed that the rice landraces could be grouped into two subpopulations, namely the indica and japonica groups. Interestingly, the alternate distribution of indica and japonica rice landraces could be found in each ecological zone. The results of AMOVA showed that on-farm conservation provides opportunities for continued differentiation and variation of landraces. Therefore, dynamic conservation measures such as on-farm conservation (which is a backup, complementary strategy to ex situ conservation) should be encouraged and enhanced, especially in crop genetic diversity centers. The results of this study offered accurate insights into short-term evolutionary processes and provided a scientific basis for on-farm management practices.

Population genetic structure of Oryza rufipogon and Oryza nivara: implications for the origin of O. nivara

Ecological speciation plays a primary role in driving species divergence and adaptation. Oryza rufipogon and Oryza nivara are two incipient species at the early stage of speciation with distinct differences in morphology, life history traits and habitat preference, and therefore provide a unique model for the study of ecological speciation. However, the population genetic structure of the ancestral O. rufipogon has been controversial despite substantial study, and the origin of the derivative O. nivara remains unclear. Here, based on sequences of 10 nuclear and two chloroplast loci from 26 wild populations across the entire geographic ranges of the two species, we conducted comprehensive analyses using population genetics, phylogeography and species distribution modelling (SDM) approaches. In addition to supporting the two previously reported major subdivisions, we detected four genetically distinct groups within O. rufipogon and found no correlation between the genetic groups and either species identity or geographical regions. The SDM clearly showed substantial change in the distribution range of O. rufipogon in history, demonstrating that the repeated extinction and colonization of local populations due to multiple glacial-interglacial cycles during the Quaternary was most likely the main factor shaping the confounding population genetic structure of O. rufipogon. Moreover, we found significant differences between the two species in climate preferences, suggestive of an important role for climatic factors in the adaptation, persistence and expansion of O. nivara. Finally, based on the genetic pattern and dynamics of the O. nivara populations, we hypothesize that O. nivara might have independently originated multiple times from different O. rufipogon populations.

Widespread and Adaptive Alterations in Genome-Wide Gene Expression Associated with Ecological Divergence of Two Oryza Species

Ecological speciation is a common mechanism by which new species arise. Despite great efforts, the role of gene expression in ecological divergence and speciation is poorly understood. Here, we conducted a genome-wide gene expression investigation of two Oryza species that are evolutionarily young and distinct in ecology and morphology. Using digital gene expression technology and the paired-end RNA sequencing method, we obtained 21,415 expressed genes across three reproduction-related tissues. Of them, approximately 8% (1,717) differed significantly in expression levels between the two species and these differentially expressed genes are randomly distributed across the genome. Moreover, 62% (1,064) of the differentially expressed genes exhibited a signature of directional selection in at least one species. Importantly, the genes with differential expression between species evolved more rapidly at the 5' flanking sequences than the genes without differential expression relative to coding sequences, suggesting that cis-regulatory changes are likely adaptive and play an important role in the ecological divergence of the two species. Finally, we showed evidence of significant differentiation between species in phenotype traits and observed that genes with differential expression were overrepresented with functional terms involving phenotypic and ecological differentiation between the two species, including reproduction- and stress-related characteristics. Our findings demonstrate that ecological speciation is associated with widespread and adaptive alterations in genome-wide gene expression and provide new insights into the importance of regulatory evolution in ecological speciation in plants.

Diversification in the South American Pampas: the genetic and morphological variation of the widespread Petunia axillaris complex (Solanaceae)

Understanding the spatiotemporal distribution of genetic variation and the ways in which this distribution is connected to the ecological context of natural populations is fundamental for understanding the nature and mode of intraspecific and, ultimately, interspecific differentiation. The Petunia axillaris complex is endemic to the grasslands of southern South America and includes three subspecies: P. a. axillaris, P. a. parodii and P. a. subandina. These subspecies are traditionally delimited based on both geography and floral morphology, although the latter is highly variable. Here, we determined the patterns of genetic (nuclear and cpDNA), morphological and ecological (bioclimatic) variation of a large number of P. axillaris populations and found that they are mostly coincident with subspecies delimitation. The nuclear data suggest that the subspecies are likely independent evolutionary units, and their morphological differences may be associated with local adaptations to diverse climatic and/or edaphic conditions and population isolation. The demographic dynamics over time estimated by skyline plot analyses showed different patterns for each subspecies in the last 100 000 years, which is compatible with a divergence time between 35 000 and 107 000 years ago between P. a. axillaris and P. a. parodii, as estimated with the IMa program. Coalescent simulation tests using Approximate Bayesian Computation do not support previous suggestions of extensive gene flow between P. a. axillaris and P. a. parodii in their contact zone.

Genetic divergence and phylogeographic history of two closely related species (Leucomeris decora and Nouelia insignis) across the 'Tanaka Line' in Southwest China

BACKGROUND

Leucomeris decora and Nouelia insignis (Asteraceae) are narrowly and allopatrically distributed species, separated by the important biogeographic boundary Tanaka Line in Southwest China. Previous morphological, cytogenetic and molecular studies suggested that L. decora is sister to N. insignis. However, it is less clear how the two species diverged, whether in full isolation or occurring gene flow across the Tanaka Line. Here, we performed a molecular study at the population level to characterize genetic differentiation and decipher phylogeographic history in two closely related species based on variation examined in plastid and nuclear DNAs using a coalescent-based approach.

RESULTS

These morphologically distinct species share plastid DNA (cpDNA) haplotypes. In contrast, Bayesian analysis of nuclear DNA (nDNA) uncovered two distinct clusters corresponding to L. decora and N. insignis. Based on the IMa analysis, no strong indication of migration was detected based on both cpDNA and nDNA sequences. The molecular data pointed to a major west-east split in nuclear DNA between the two species corresponding with the Tanaka Line. The coalescent time estimate for all cpDNA haplotypes dated to the Mid-Late Pleistocene. The estimated demographic parameters showed that the population size of L. decora was similar to that of N. insignis and both experienced limited demographic fluctuations recently.

CONCLUSIONS

The study revealed comprehensive species divergence and phylogeographic histories of N. insignis and L. decora divided by the Tanaka Line. The phylogeographic pattern inferred from cpDNA reflected ancestrally shared polymorphisms without post-divergence gene flow between species. The marked genealogical lineage divergence in nDNA provided some indication of Tanaka Line for its role as a barrier to plant dispersal, and lent support to its importance in promoting strong population structure and allopatric divergence.

Rice genomes recorded ancient pararetrovirus activities: Virus genealogy and multiple origins of endogenization during rice speciation

Viral fossils in rice genomes are a best entity to understand ancient pararetrovirus activities through host plant history because of our advanced knowledge of the genomes and evolutionary history with rice and its related species. Here, we explored organization, geographic origins and genealogy of rice pararetroviruses, which were turned into endogenous rice tungro bacilliform virus-like (eRTBVL) sequences. About 300 eRTBVL sequences from three representative rice genomes were clearly classified into six families. Most of the endogenization events of the eRTBVLs were initiated before differentiation of the rice progenitor (> 160,000 years ago). We successfully followed the genealogy of old relic viruses during rice speciation, and inferred the geographical origins for these viruses. Possible virus genomic sequences were explained mostly by recombinations between different virus families. Interestingly, we discovered that only a few recombination events among the numerous occasions had determined the virus genealogy.

Pliocene intraspecific divergence and Plio-Pleistocene range expansions within Picea likiangensis (Lijiang spruce), a dominant forest tree of the Qinghai-Tibet Plateau

A knowledge of intraspecific divergence and range dynamics of dominant forest trees in response to past geological and climate change is of major importance to an understanding of their recent evolution and demography. Such knowledge is informative of how forests were affected by environmental factors in the past and may provide pointers to their response to future environmental change. However, genetic signatures of such historical events are often weak at individual loci due to large effective population sizes and long generation times of forest trees. This problem can be overcome by analysing genetic variation across multiple loci. We used this approach to examine intraspecific divergence and past range dynamics in the conifer Picea likiangensis, a dominant tree of forests occurring in eastern and southern areas of the Qinghai-Tibet Plateau (QTP). We sequenced 13 nuclear loci, two mitochondrial DNA regions and three plastid (chloroplast) DNA regions in 177 individuals sampled from 22 natural populations of this species, and tested the hypothesis that its evolutionary history was markedly affected by Pliocene QTP uplifts and Quaternary climatic oscillations. Consistent with the taxonomic delimitation of the three morphologically divergent varieties examined, all individuals clustered into three genetic groups with intervariety admixture detected in regions of geographical overlap. Divergence between varieties was estimated to have occurred within the Pliocene and ecological niche modelling based on 20 ecological variables suggested that niche differentiation was high. Furthermore, modelling of population-genetic data indicated that two of the varieties (var. rubescens and var. linzhiensis) expanded their population sizes after the largest Quaternary glaciation in the QTP, while expansion of the third variety (var. likiangensis) began prior to this, probably following the Pliocene QTP uplift. These findings point to the importance of geological and climatic changes during the Pliocene and Pleistocene as causes of intraspecific diversification and range shifts of dominant tree species in the QTP biodiversity hot spot region.

Local differentiation amidst extensive allele sharing in Oryza nivara and O. rufipogon

Genetic variation patterns within and between species may change along geographic gradients and at different spatial scales. This was revealed by microsatellite data at 29 loci obtained from 119 accessions of three Oryza series Sativae species in Asia Pacific: Oryza nivara Sharma and Shastry, O. rufipogon Griff., and O. meridionalis Ng. Genetic similarities between O. nivara and O. rufipogon across their distribution are evident in the clustering and ordination results and in the large proportion of shared alleles between these taxa. However, local-level species separation is recognized by Bayesian clustering and neighbor-joining analyses. At the regional scale, the two species seem more differentiated in South Asia than in Southeast Asia as revealed by F ST analysis. The presence of strong gene flow barriers in smaller spatial units is also suggested in the analysis of molecular variance (AMOVA) results where 64% of the genetic variation is contained among populations (as compared to 26% within populations and 10% among species). Oryza nivara (H E = 0.67) exhibits slightly lower diversity and greater population differentiation than O. rufipogon (H E = 0.70). Bayesian inference identified four, and at a finer structural level eight, genetically distinct population groups that correspond to geographic populations within the three taxa. Oryza meridionalis and the Nepalese O. nivara seemed diverged from all the population groups of the series, whereas the Australasian O. rufipogon appeared distinct from the rest of the species.

Pleistocene climate change and the origin of two desert plant species, Pugionium cornutum and Pugionium dolabratum (Brassicaceae), in northwest China

Pleistocene climate change has had an important effect in shaping intraspecific genetic variation in many species; however, its role in driving speciation is less clear. We examined the possibility of a Pleistocene origin of the only two representatives of the genus Pugionium (Brassicaceae), Pugionium cornutum and Pugionium dolabratum, which occupy different desert habitats in northwest China. We surveyed sequence variation for internal transcribed spacer (ITS), three chloroplast (cp) DNA fragments, and eight low-copy nuclear genes among individuals sampled from 11 populations of each species across their geographic ranges. One ITS mutation distinguished the two species, whereas mutations in cpDNA and the eight low-copy nuclear gene sequences were not species-specific. Although interspecific divergence varied greatly among nuclear gene sequences, in each case divergence was estimated to have occurred within the Pleistocene when deserts expanded in northwest China. Our findings point to the importance of Pleistocene climate change, in this case an increase in aridity, as a cause of speciation in Pugionium as a result of divergence in different habitats that formed in association with the expansion of deserts in China.

Understanding rice domestication and implications for cultivar improvement

Considerable insights were recently gained into the history and process of rice domestication. It becomes increasingly clear that artificial and natural selections coupled with extensive introgression have shaped the genomes of cultivated rice. The interplay of these evolutionary forces gave rise to the cultivated species, Oryza sativa, with divergent genomic backgrounds from two wild species, O. rufipogon and O. nivara, governed by a set of domestication alleles which had originated primarily at one location of initial cultivation. The mechanistic understanding of domestication suggests that the combination of quantitative trait locus mapping, genome-wide association study, and genome scan will be effective means for discovering potentially valuable alleles from the cultivated and wild species. The accumulation and appropriate sampling of germplasm collections for these analyses should effectively grow the useful allele pool, which combined with molecular breeding may get to a point literally triggering the re-domestication of rice varieties for sustainable food production.

Evidence for asymmetrical divergence-gene flow of nuclear loci, but not mitochondrial loci, between seabird sister species: blue-footed (Sula nebouxii) and Peruvian (S. variegata) boobies

Understanding the process of speciation requires understanding how gene flow influences divergence. Recent analyses indicate that divergence can take place despite gene flow and that the sex chromosomes can exhibit different levels of gene flow than autosomes and mitochondrial DNA. Using an eight marker dataset including autosomal, z-linked, and mitochondrial loci we tested the hypothesis that blue-footed (Sula nebouxii) and Peruvian (S. variegata) boobies diverged from their common ancestor with gene flow, paying specific attention to the differences in gene flow estimates from nuclear and mitochondrial markers. We found no gene flow at mitochondrial markers, but found evidence from the combined autosomal and z-linked dataset that blue-footed and Peruvian boobies experienced asymmetrical gene flow during or after their initial divergence, predominantly from Peruvian boobies into blue-footed boobies. This gene exchange may have occurred either sporadically between periods of allopatry, or regularly throughout the divergence process. Our results add to growing evidence that diverging species can remain distinct but exchange genes.

Origin of Oryza sativa in China inferred by nucleotide polymorphisms of organelle DNA

China is rich of germplasm resources of common wild rice (Oryza rufipogon Griff.) and Asian cultivated rice (O. sativa L.) which consists of two subspecies, indica and japonica. Previous studies have shown that China is one of the domestication centers of O. sativa. However, the geographic origin and the domestication times of O. sativa in China are still under debate. To settle these disputes, six chloroplast loci and four mitochondrial loci were selected to examine the relationships between 50 accessions of Asian cultivated rice and 119 accessions of common wild rice from China based on DNA sequence analysis in the present study. The results indicated that Southern China is the genetic diversity center of O. rufipogon and it might be the primary domestication region of O. sativa. Molecular dating suggested that the two subspecies had diverged 0.1 million years ago, much earlier than the beginning of rice domestication. Genetic differentiations and phylogeography analyses indicated that indica was domesticated from tropical O. rufipogon while japonica was domesticated from O. rufipogon which located in higher latitude. These results provided molecular evidences for the hypotheses of (i) Southern China is the origin center of O. sativa in China and (ii) the two subspecies of O. sativa were domesticated multiple times.

Domestication and geographic origin of Oryza sativa in China: insights from multilocus analysis of nucleotide variation of O. sativa and O. rufipogon

Previous studies have indicated that China is one of the domestication centres of Asian cultivated rice (Oryza sativa), and common wild rice (O. rufipogon) is the progenitor of O. sativa. However, the number of domestication times and the geographic origin of Asian cultivated rice in China are still under debate. In this study, 100 accessions of Asian cultivated rice and 111 accessions of common wild rice in China were selected to examine the relationship between O. sativa and O. rufipogon and thereby infer the domestication and evolution of O. sativa in China through sequence analyses of six gene regions, trnC-ycf6 in chloroplast genomes, cox3 in mitochondrial genomes and ITS, Ehd1, Waxy, Hd1 in nuclear genomes. The results indicated that the two subspecies of O. sativa (indica and japonica) were domesticated independently from different populations of O. rufipogon with gene flow occurring later from japonica to indica; Southern China was the genetic diversity centre of O. rufipogon, and the Pearl River basin near the Tropic of Cancer was the domestication centre of O. sativa in China.

Chloroplast genome sequence confirms distinctness of Australian and Asian wild rice

Cultivated rice (Oryza sativa) is an AA genome Oryza species that was most likely domesticated from wild populations of O. rufipogon in Asia. O. rufipogon and O. meridionalis are the only AA genome species found within Australia and occur as widespread populations across northern Australia. The chloroplast genome sequence of O. rufipogon from Asia and Australia and O. meridionalis and O. australiensis (an Australian member of the genus very distant from O. sativa) was obtained by massively parallel sequencing and compared with the chloroplast genome sequence of domesticated O. sativa. Oryza australiensis differed in more than 850 sites single nucleotide polymorphism or indel from each of the other samples. The other wild rice species had only around 100 differences relative to cultivated rice. The chloroplast genomes of Australian O. rufipogon and O. meridionalis were closely related with only 32 differences. The Asian O. rufipogon chloroplast genome (with only 68 differences) was closer to O. sativa than the Australian taxa (both with more than 100 differences). The chloroplast sequences emphasize the genetic distinctness of the Australian populations and their potential as a source of novel rice germplasm. The Australian O. rufipogon may be a perennial form of O. meridionalis.

Genetic divergence, range expansion and possible homoploid hybrid speciation among pine species in Northeast China

Although homoploid hybrid speciation in plants is probably more common than previously realized, there are few well-documented cases of homoploid hybrid origin in conifers. We examined genetic divergence between two currently widespread pines in Northeast China, Pinus sylvestris var. mongolica and Pinus densiflora, and also whether two narrowly distributed pines in the same region, Pinus funebris and Pinus takahasii, might have originated from the two widespread species by homoploid hybrid speciation. Our results, based on population genetic analysis of chloroplast (cp), mitochondrial (mt) DNA, and nuclear gene sequence variation, showed that the two widespread species were divergent for both cp- and mtDNA variation, and also for haplotype variation at two of eight nuclear gene loci surveyed. Our analysis further indicated that P. sylvestris var. mongolica and P. densiflora remained allopatric during the most severe Quaternary glacial period that occurred in Northeast China, but subsequently exhibited rapid range expansions. P. funebris and P. takahasii, were found to contain a mixture of chlorotypes and nuclear haplotypes that distinguish P. sylvestris var. mongolica and P. densiflora, in support of the hypothesis that they possibly originated via homoploid hybrid speciation following secondary contact and hybridization between P. sylvestris var. mongolica and P. densiflora.

Pleistocene climatic oscillations and the speciation history of an alpine endemic and a widespread arctic-alpine plant

Even in cases in which geographic isolation appears to have driven the speciation of regional endemics, range shifts during the Pleistocene climatic oscillations may also have influenced their evolutionary history. Elucidating speciation history can provide novel insights into evolutionary dynamics following climatic oscillations. We demonstrated a sister relationship between the Japanese alpine endemic Cardamine nipponica and the currently allopatric, widespread arctic-alpine Cardamine bellidifolia (Brassicaceae) based on internal transcribed spacer (ITS) sequences and 10 other nuclear genes. Speciation history was inferred using demographic parameters under the isolation with migration model. The estimated demographic parameters showed that the population size of C. nipponica was similar to that of C. bellidifolia and that gene flow occurred exclusively from C. nipponica to C. bellidifolia after speciation. The inferred speciation history, which included gene flow, suggests that geographic barriers between the peripheral C. nipponica and the widespread C. bellidifolia were reduced during the Pleistocene. The asymmetric introgression implies that genetic isolation may have been involved in the speciation of C. nipponica. Our results suggest that even currently allopatric species may not have diverged solely under geographic isolation, and that their evolutionary history may have been influenced by Pleistocene range dynamics.