Genomic patterns of nucleotide diversity in divergent populations of U.S. weedy rice

Background  

Weedy rice (red rice), a conspecific weed of cultivated rice (Oryza sativa L.), is a significant problem throughout the world and an emerging threat in regions where it was previously absent. Despite belonging to the same species complex as domesticated rice and its wild relatives, the evolutionary origins of weedy rice remain unclear. We use genome-wide patterns of single nucleotide polymorphism (SNP) variation in a broad geographic sample of weedy, domesticated, and wild Oryza samples to infer the origin and demographic processes influencing U.S. weedy rice evolution.     

Avian β-defensin variation in bottlenecked populations: the Seychelles warbler and other congeners

β-defensins are important components of the vertebrate innate immune system responsible for encoding a variety of anti-microbial peptides. Pathogen-mediated selection is thought to act on immune genes and potentially maintain allelic variation in the face of genetic drift. The Seychelles warbler, Acrocephalus sechellensis, is an endemic passerine that underwent a recent bottleneck in its last remaining population, resulting in a considerable reduction in genome-wide variation. We genotyped avian β-defensin (AvBD) genes in contemporary (2000–2008) and museum samples (1876–1940) of the Seychelles warbler to investigate whether immunogenetic variation was lost through this bottleneck, and examined AvBD variation across four other Acrocephalus species with varying demographic histories. No variation was detected at four of the six AvBD loci screened in the post-bottleneck population of Seychelles warbler, but two silent nucleotide polymorphisms were identified at AvBD8 and one potentially functional amino-acid variation was observed at AvBD11. Variation in the Seychelles warbler was significantly lower than in the mainland migratory congeneric species investigated, but it similar to that found in other bottlenecked species. In addition, screening AvBD7 in 15 museum specimens of Seychelles warblers sampled prior to the bottleneck (1877–1905) revealed that this locus possessed two alleles previously, compared to the single allele in the contemporary population. Overall, the results show that little AvBD variation remains in the Seychelles warbler, probably as a result of having low AvBD diversity historically rather than the loss of variation due to drift associated with past demographic history. Given the limited pathogen fauna, this lack of variation at the AvBD loci may currently not pose a problem for this isolate population of Seychelles warblers, but it may be detrimental to the species’ long-term survival if new pathogens reach the population in the future.     

Deleterious amino acid polymorphisms in Arabidopsis thaliana and rice

Plant genetic diversity has been mainly investigated with neutral markers, but large-scale DNA sequencing projects now enable the identification and analysis of different classes of genetic polymorphisms, such as non-synonymous single nucleotide polymorphisms (nsSNPs) in protein coding sequences. Using the SIFT and MAPP programs to predict whether nsSNPs are tolerated (i.e., effectively neutral) or deleterious for protein function, genome-wide nsSNP data from Arabidopsis thaliana and rice were analyzed. In both species, about 20% of polymorphic sites with nsSNPs were classified as deleterious; they segregate at lower allele frequencies than tolerated nsSNPs due to purifying selection. Furthermore, A. thaliana accessions from marginal populations show a higher relative proportion of deleterious nsSNPs, which likely reflects differential selection or demographic effects in subpopulations. To evaluate the sensitivity of predictions, genes from model and crop plants with known functional effects of nsSNPs were inferred with the algorithms. The programs predicted about 70% of nsSNPs correctly as tolerated or deleterious, i.e., as having a functional effect. Forward-in-time simulations of bottleneck and domestication models indicated a high power to detect demographic effects on nsSNP frequencies in sufficiently large datasets. The results indicate that nsSNPs are useful markers for analyzing genetic diversity in plant genetic resources and breeding populations to infer natural/artificial selection and genetic drift.     

Levels and patterns of nucleotide variation in domestication QTL regions on rice chromosome 3 suggest lineage-specific selection

Oryza sativa or Asian cultivated rice is one of the major cereal grass species domesticated for human food use during the Neolithic. Domestication of this species from the wild grass Oryza rufipogon was accompanied by changes in several traits, including seed shattering, percent seed set, tillering, grain weight, and flowering time. Quantitative trait locus (QTL) mapping has identified three genomic regions in chromosome 3 that appear to be associated with these traits. We would like to study whether these regions show signatures of selection and whether the same genetic basis underlies the domestication of different rice varieties. Fragments of 88 genes spanning these three genomic regions were sequenced from multiple accessions of two major varietal groups in O. sativa--indica and tropical japonica--as well as the ancestral wild rice species O. rufipogon. In tropical japonica, the levels of nucleotide variation in these three QTL regions are significantly lower compared to genome-wide levels, and coalescent simulations based on a complex demographic model of rice domestication indicate that these patterns are consistent with selection. In contrast, there is no significant reduction in nucleotide diversity in the homologous regions in indica rice. These results suggest that there are differences in the genetic and selective basis for domestication between these two Asian rice varietal groups.     

Approximate Bayesian inference reveals evidence for a recent, severe bottleneck in a Netherlands population of Drosophila melanogaster

Genome-wide nucleotide variation in non-African populations of Drosophila melanogaster is a subset of variation found in East sub-Saharan African populations, suggesting a bottleneck in the history of the former. We implement an approximate Bayesian approach to infer the timing, duration, and severity of this putative bottleneck and ask whether this inferred model is sufficient to account for patterns of variability observed at 115 loci scattered across the X chromosome. We estimate a recent bottleneck 0.006N(e) generations ago, somewhat further in the past than suggested by biogeographical evidence. Using various proposed statistical tests, we find that this bottleneck model is able to predict the majority of observed features of diversity and linkage disequilibrium in the data. Thus, while precise estimates of bottleneck parameters (like inferences of selection) are sensitive to model assumptions, our results imply that it may be unnecessary to invoke frequent selective sweeps associated with the dispersal of D. melanogaster from Africa to explain patterns of variability in non-African populations.     

Genome-wide association mapping of virulence gene in rice blast fungus Magnaporthe oryzae using a genotyping by sequencing approach

Magnaporthe oryzae is a fungal pathogen causing blast disease in many plant species. In this study, seventy three isolates of M. oryzae collected from rice (Oryza sativa) in 1996–2014 were genotyped using a genotyping-by-sequencing approach to detect genetic variation. An association study was performed to identify single nucleotide polymorphisms (SNPs) associated with virulence genes using 831 selected SNP and infection phenotypes on local and improved rice varieties. Population structure analysis revealed eight subpopulations. The division into eight groups was not related to the degree of virulence. Association mapping showed five SNPs associated with fungal virulence on chromosome 1, 2, 3, 4 and 7. The SNP on chromosome 1 was associated with virulence against RD6-Pi7 and IRBL7-M which might be linked to the previously reported AvrPi7.     

Revisiting protein heterozygosity in plants��nucleotide diversity in allozyme coding genes of conifer Pinus sylvestris

Allozyme variation has been and continues to be a major source of information on the level of genetic variation among plant species. Deciphering the molecular basis of electrophoretic variation is essential for understanding the forces affecting the protein level variation. In this study, the relationship between allozyme heterozygosity and nucleotide diversity in plants is investigated among and within species. Allozyme and nucleotide diversity in 27 plant species was reviewed. At the multilocus level, the two methods are congruent: a clear correlation between the two measures of genetic diversity among plant species was observed, strengthening the view that effective population size is the major determinant of genome-wide diversity. Nucleotide diversity at six allozyme coding genes (6pgdB, aco, gdh, gotC, mdhA, and mdhB) in conifer Pinus sylvestris was investigated jointly with electrophoretic data. Single non-synonymous charge-changing mutations were found together with electrophoretic alleles that consequently were mutationally unique. Synonymous site nucleotide diversity (point estimate of θ _W—0.009 per bp) and silent site divergence from Pinus pinaster at allozyme coding loci were at comparable levels with other loci in the species. Linkage disequilibrium was extensive compared to earlier estimates from P. sylvestris and other trees, spanning several kilobases. Allozyme coding genes had an excess of closely related haplotypes whose frequency has recently increased possibly as a result of partial selective sweeps or balancing selection, but complex demographic effects cannot be excluded.     

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.     

Analysis of Starch Gene Diversity in the Wild Relatives of Oryza sativa

Genetic loci influencing traits important to humans have been selected during crop domestication. The starch properties of rice influence the ease of cooking and attractiveness of rice as a human food. Starch biosynthesis genes likely to influence starch properties in the grain were compared in wild and domesticated rice genotypes. Sequence variation was investigated in starch biosynthesis gene exons that have been reported to have a direct influence on rice amylose content, gelatinization temperature, and amylopectin chain length. Exons 6 and 10 of GBSSI, exon 8 of SSIIa and exons 11, 13, 14 and 16 of SBEIIb were amplified and sequenced from 13 wild Oryza species encompassing genome types AA to HHJJ. Thirty two single nucleotide polymorphisms (SNPs) were identified in the exons of GBSSI; 176 in exon 8 of SSIIa, and 43 in SBEIIb, giving a total of 251 SNPs among the species. Eighty six of these SNP caused changes in the encoded amino acid, of which 28 were missense mutations that affected highly conserved amino acids within the protein sequence of GBSSI, SSIIa or SBEIIb. Two indels were identified in Potamophila parviflora, a close relative of Zizania palustris, a North American native wild rice. Most of the nucleotide variations and non-conservative changes were observed in the genomes other than the AA genome species. This represents a genetic resource for use in rice starch manipulation. The impact of human selection at these loci can be deduced by comparison of modern cultivated genotypes with their wild progenitors.     

Sequence variation of avirulence gene AVR-Pita1 in rice blast fungus, Magnaporthe oryzae

The interaction between rice, Oryza sativa, and rice blast fungus, Magnaporthe oryzae, is triggered by an interaction between the protein products of the host resistant gene, and the pathogen avirulence gene. This interaction follows the ‘gene-for-gene' concept. The resistant gene has effectively protected rice plants from rice blast infection. However, the resistant genes usually break down several years after the release of the resistant rice varieties because the fungus has evolved to new races. The objective of this study is to investigate the nucleotide sequence variation of the AVR-Pita1 gene that influences the adaption of rice blast fungus to overcome the resistant gene, Pi-ta. Thirty rice blast fungus isolates were collected in 2005 and 2010 from infected rice plants in northern and northeastern Thailand. The nucleotide sequences of AVR-Pita1 were amplified and analyzed. Phylogenetic analysis was conducted using the MEGA 5.0 program. The results showed a high level of nucleotide sequence polymorphisms and the positive genetic selection pressure in Thai rice blast isolates. The details of sequence variation analysis were described in this article. The information from this study can be used for rice blast resistant breeding program in the future.     

Seeing red: the origin of grain pigmentation in US weedy rice

Weedy forms of crop species infest agricultural fields worldwide and are a leading cause of crop losses, yet little is known about how these weeds evolve. Red rice (Oryza sativa), a major weed of cultivated rice fields in the US, is recognized by the dark‐pigmented grain that gives it its common name. Studies using neutral molecular markers have indicated a close relationship between US red rice and domesticated rice, suggesting that the weed may have originated through reversion of domesticated rice to a feral form. We have tested this reversion hypothesis by examining molecular variation at Rc, the regulatory gene responsible for grain pigmentation differences between domesticated and wild rice. Loss‐of‐function mutations at Rc account for the absence of proanthocyanidin pigments in cultivated rice grains, and the major rc domestication allele has been shown to be capable of spontaneous reversion to a functional form through additional mutations at the Rc locus. Using a diverse sample of 156 weedy, domesticated and wild Oryzas, we analysed DNA sequence variation at Rc and its surrounding 4 Mb genomic region. We find that reversion of domestication alleles does not account for the pigmented grains of weed accessions; moreover, we find that haplotypes characterizing the weed are either absent or very rare in cultivated rice. Sequences from genomic regions flanking Rc are consistent with a genomic footprint of the rc selective sweep in cultivated rice, and they are compatible with a close relationship of red rice to Asian Oryzas that have never been cultivated in the US.     

Demographic history of a common pioneer tree,Zanthoxylum ailanthoides,reconstructed using isolation-with-migration model

The Japanese prickly ash, Zanthoxylum ailanthoides Siebold & Zucc. (Rutaceae), is one of the most common pioneer tree species that grow in disturbed areas, such as canopy gaps, in Japanese warm-temperate evergreen oak forests. The strong genetic structure of its current population might suggest long-term isolation of subpopulations by demographic events, in addition to ecological features of this species. Analyses of genome-wide nucleotide variation using model-based approaches are necessary to achieve a good understanding of the demographic history of a species. In this study, we analyzed the nucleotide variation in natural populations of Z. ailanthoides using a computer program that applied the isolation-with-migration model to quantitatively infer the demographic history. Nucleotide variations at 10 or 26 nuclear loci in six populations from a wide range of the species distribution were analyzed. The maximum likelihood estimate of the divergence time between the current populations in the two hypothetical refugia of Z. ailanthoides was approximately 24,000 years. Unexpectedly, the estimated size of the ancestral population was larger than the sizes of the two current populations. The results suggested a relatively recent divergence of these two populations and rapid formation of strong genetic structure among subpopulations. The large ancestral population may indicate a more complex demographic history during or after the last glacial period than the simple isolation-with-migration model implies.     

Allele Distributions at Hybrid Incompatibility Loci Facilitate the Potential for Gene Flow between Cultivated and Weedy Rice in the US

The accumulation of independent mutations over time in two populations often leads to reproductive isolation. Reproductive isolation between diverging populations may be reinforced by barriers that occur either pre- or postzygotically. Hybrid sterility is the most common form of postzygotic isolation in plants. Four postzygotic sterility loci, comprising three hybrid sterility systems (Sa, s5, DPL), have been recently identified in Oryza sativa. These loci explain, in part, the limited hybridization that occurs between the domesticated cultivated rice varieties, O. sativa spp. japonica and O. sativa spp. indica. In the United States, cultivated fields of japonica rice are often invaded by conspecific weeds that have been shown to be of indica origin. Crop-weed hybrids have been identified in crop fields, but at low frequencies. Here we examined the possible role of these hybrid incompatibility loci in the interaction between cultivated and weedy rice. We identified a novel allele at Sa that seemingly prevents loss of fertility in hybrids. Additionally, we found wide-compatibility type alleles at strikingly high frequencies at the Sa and s5 loci in weed groups, and a general lack of incompatible alleles between crops and weeds at the DPL loci. Our results suggest that weedy individuals, particularly those of the SH and BRH groups, should be able to freely hybridize with the local japonica crop, and that prezygotic factors, such as differences in flowering time, have been more important in limiting weed-crop gene flow in the past. As the selective landscape for weedy rice changes due to increased use of herbicide resistant strains of cultivated rice, the genetic barriers that hinder indica-japonica hybridization cannot be counted on to limit the flow of favorable crop genes into weeds.     

Inferring the demographic history and rate of adaptive substitution in Drosophila

An important goal of population genetics is to determine the forces that have shaped the pattern of genetic variation in natural populations. We developed a maximum likelihood method that allows us to infer demographic changes and detect recent positive selection (selective sweeps) in populations of varying size from DNA polymorphism data. Applying this approach to single nucleotide polymorphism data at more than 250 noncoding loci on the X chromosome of Drosophila melanogaster from an (ancestral) African population and a (derived) European, we found that the African population expanded about 60,000 y ago and that the European population split off from the African lineage about 15,800 y ago, thereby suffering a severe population size bottleneck. We estimated that about 160 beneficial mutations (with selection coefficients s between 0.05% and 0.5%) were fixed in the euchromatic portion of the X in the African population since population size expansion, and about 60 mutations (with s around 0.5%) in the diverging European lineage.     

Single nucleotide polymorphisms generated by genotyping by sequencing to characterize genome-wide diversity,linkage disequilibrium,and selective sweeps in cultivated watermelon

Background

A large single nucleotide polymorphism (SNP) dataset was used to analyze genome-wide diversity in a diverse collection of watermelon cultivars representing globally cultivated, watermelon genetic diversity. The marker density required for conducting successful association mapping depends on the extent of linkage disequilibrium (LD) within a population. Use of genotyping by sequencing reveals large numbers of SNPs that in turn generate opportunities in genome-wide association mapping and marker-assisted selection, even in crops such as watermelon for which few genomic resources are available. In this paper, we used genome-wide genetic diversity to study LD, selective sweeps, and pairwise F_ST distributions among worldwide cultivated watermelons to track signals of domestication.

Results

We examined 183 Citrullus lanatus var. lanatus accessions representing domesticated watermelon and generated a set of 11,485 SNP markers using genotyping by sequencing. With a diverse panel of worldwide cultivated watermelons, we identified a set of 5,254 SNPs with a minor allele frequency of ≥ 0.05, distributed across the genome. All ancestries were traced to Africa and an admixture of various ancestries constituted secondary gene pools across various continents. A sliding window analysis using pairwise F_ST values was used to resolve selective sweeps. We identified strong selection on chromosomes 3 and 9 that might have contributed to the domestication process. Pairwise analysis of adjacent SNPs within a chromosome as well as within a haplotype allowed us to estimate genome-wide LD decay. LD was also detected within individual genes on various chromosomes. Principal component and ancestry analyses were used to account for population structure in a genome-wide association study. We further mapped important genes for soluble solid content using a mixed linear model.

Conclusions

Information concerning the SNP resources, population structure, and LD developed in this study will help in identifying agronomically important candidate genes from the genomic regions underlying selection and for mapping quantitative trait loci using a genome-wide association study in sweet watermelon.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-767) contains supplementary material, which is available to authorized users.     

Comparative phylogeography among hydrothermal vent species along the East Pacific Rise reveals vicariant processes and population expansion in the South

The use of sequence polymorphism from individual mitochondrial genes to infer past demography has recently proved controversial because of the recurrence of selective sweeps acting over genes and the need for unlinked multilocus data sets. However, comparative analyses using several species for one gene and/or multiple genes for one species can serve as a test for potential selective effects and clarify our understanding of historical demographic effects. This study compares nucleotide polymorphisms in mitochondrial cytochrome oxidase I across seven deep-sea hydrothermal vent species that live along the volcanically active East Pacific Rise. Approximate Bayesian Computation (ABC) method, developed to trace shared vicariant events across species pairs, indicates the occurrence of two across species divergence times, and suggests that the present geographical patterns of genetic differentiation may be explained by two periods of significant population isolation. The oldest period dates back 11.6 Ma and is associated with the vent limpet Lepetodrilus elevatus , while the most recent period of isolation is 1.3 Ma, which apparently affected all species examined and coincides with a transition zone across the equator. Moreover, significant negative Tajima's D and star-like networks were observed for all southern lineages, suggesting that these lineages experienced a concomitant demographic and geographical expansion about 100 000–300 000 generations ago. This expansion may have initiated from a wave of range expansions during the secondary colonization of new sites along the Southern East Pacific Rise (founder effects below the equator) or recurrent bottleneck events because of the increase of eruptive phases associated with the higher spreading rates of the ridge in this region.     

Unraveling the Complex Hybrid Ancestry and Domestication History of Cultivated Strawberry

Cultivated strawberry (Fragaria × ananassa) is one of our youngest domesticates, originating in early eighteenth-century Europe from spontaneous hybrids between wild allo-octoploid species (Fragaria chiloensis and Fragaria virginiana). The improvement of horticultural traits by 300 years of breeding has enabled the global expansion of strawberry production. Here, we describe the genomic history of strawberry domestication from the earliest hybrids to modern cultivars. We observed a significant increase in heterozygosity among interspecific hybrids and a decrease in heterozygosity among domesticated descendants of those hybrids. Selective sweeps were found across the genome in early and modern phases of domestication—59–76% of the selectively swept genes originated in the three less dominant ancestral subgenomes. Contrary to the tenet that genetic diversity is limited in cultivated strawberry, we found that the octoploid species harbor massive allelic diversity and that F. × ananassa harbors as much allelic diversity as either wild founder. We identified 41.8 M subgenome-specific DNA variants among resequenced wild and domesticated individuals. Strikingly, 98% of common alleles and 73% of total alleles were shared between wild and domesticated populations. Moreover, genome-wide estimates of nucleotide diversity were virtually identical in F. chiloensis,F. virginiana, and F. × ananassa (π = 0.0059–0.0060). We found, however, that nucleotide diversity and heterozygosity were significantly lower in modern F. × ananassa populations that have experienced significant genetic gains and have produced numerous agriculturally important cultivars.     

The complex genetic architecture of the metabolome

Discovering links between the genotype of an organism and its metabolite levels can increase our understanding of metabolism, its controls, and the indirect effects of metabolism on other quantitative traits. Recent technological advances in both DNA sequencing and metabolite profiling allow the use of broad-spectrum, untargeted metabolite profiling to generate phenotypic data for genome-wide association studies that investigate quantitative genetic control of metabolism within species. We conducted a genome-wide association study of natural variation in plant metabolism using the results of untargeted metabolite analyses performed on a collection of wild Arabidopsis thaliana accessions. Testing 327 metabolites against >200,000 single nucleotide polymorphisms identified numerous genotype-metabolite associations distributed non-randomly within the genome. These clusters of genotype-metabolite associations (hotspots) included regions of the A. thaliana genome previously identified as subject to recent strong positive selection (selective sweeps) and regions showing trans-linkage to these putative sweeps, suggesting that these selective forces have impacted genome-wide control of A. thaliana metabolism. Comparing the metabolic variation detected within this collection of wild accessions to a laboratory-derived population of recombinant inbred lines (derived from two of the accessions used in this study) showed that the higher level of genetic variation present within the wild accessions did not correspond to higher variance in metabolic phenotypes, suggesting that evolutionary constraints limit metabolic variation. While a major goal of genome-wide association studies is to develop catalogues of intraspecific variation, the results of multiple independent experiments performed for this study showed that the genotype-metabolite associations identified are sensitive to environmental fluctuations. Thus, studies of intraspecific variation conducted via genome-wide association will require analyses of genotype by environment interaction. Interestingly, the network structure of metabolite linkages was also sensitive to environmental differences, suggesting that key aspects of network architecture are malleable.     

Natural Selection and Recombination Rate Variation Shape Nucleotide Polymorphism Across the Genomes of Three Related Populus Species

A central aim of evolutionary genomics is to identify the relative roles that various evolutionary forces have played in generating and shaping genetic variation within and among species. Here we use whole-genome resequencing data to characterize and compare genome-wide patterns of nucleotide polymorphism, site frequency spectrum, and population-scaled recombination rates in three species of Populus: Populus tremula, P. tremuloides, and P. trichocarpa. We find that P. tremuloides has the highest level of genome-wide variation, skewed allele frequencies, and population-scaled recombination rates, whereas P. trichocarpa harbors the lowest. Our findings highlight multiple lines of evidence suggesting that natural selection, due to both purifying and positive selection, has widely shaped patterns of nucleotide polymorphism at linked neutral sites in all three species. Differences in effective population sizes and rates of recombination largely explain the disparate magnitudes and signatures of linked selection that we observe among species. The present work provides the first phylogenetic comparative study on a genome-wide scale in forest trees. This information will also improve our ability to understand how various evolutionary forces have interacted to influence genome evolution among related species.