Molecular evidence on the origin and evolution of glutinous rice

Glutinous rice is a major type of cultivated rice with long-standing cultural importance in Asia. A mutation in an intron 1 splice donor site of the Waxy gene is responsible for the change in endosperm starch leading to the glutinous phenotype. Here we examine an allele genealogy of the Waxy locus to trace the evolutionary and geographical origins of this phenotype. On the basis of 105 glutinous and nonglutinous landraces from across Asia, we find evidence that the splice donor mutation has a single evolutionary origin and that it probably arose in Southeast Asia. Nucleotide diversity measures indicate that the origin of glutinous rice is associated with reduced genetic variation characteristic of selection at the Waxy locus; comparison with an unlinked locus, RGRC2, confirms that this pattern is specific to Waxy. In addition, we find that many nonglutinous varieties in Northeast Asia also carry the splice donor site mutation, suggesting that partial suppression of this mutation may have played an important role in the development of Northeast Asian nonglutinous rice. This study demonstrates the utility of phylogeographic approaches for understanding trait diversification in crops, and it contributes to growing evidence on the importance of modifier loci in the evolution of domestication traits.     

Generation of new glutinous rice by CRISPR/Cas9‐targeted mutagenesis of the Waxy gene in elite rice varieties

In rice, amylose content(AC) is controlled by a single dominant Waxy gene. We used Clustered Regularly Interspaced Short Palindromic Repeats(CRISPR)/CRISPR-associated 9(Cas_9) to introduce a loss-of-function mutation into the Waxy gene in two widely cultivated elite japonica varieties. Our results show that mutations in the Waxy gene reduce AC and convert the rice into glutinous ones without affecting other desirable agronomic traits, offering an effective and easy strategy to improve glutinosity in elite varieties. Importantly, we successfully removed the transgenes from the progeny. Our study provides an example of generating improved crops with potential for commercialization, by editing a gene of interest directly in elite crop varieties.     

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.     

Molecular evidence for post-domestication selection in the <Emphasis Type="Italic">Waxy</Emphasis> gene of Chinese waxy maize

Waxy maize was first reported in China in 1909 and is mainly used in food production in Asia. The evidence for strong domestication selection in the Waxy locus of rice and a selective sweep around its genomic region make us to wonder whether there has been similar selection in Waxy in glutinous maize. To address this issue, DNA sequences of Waxy, three flanking genes and an unlinked gene (Adh1) of 30 accessions sampled from Chinese waxy maize accessions, including representative landraces and inbred lines, were determined in this study. Sharp reduction of nucleotide diversity and significant neutrality tests (Tajima’s D and Fu and Li’s F*) were observed in the Waxy locus in Chinese waxy maize but not in nonglutinous maize; comparison with the unlinked gene confirmed that this pattern was different to Waxy. Sequence analysis across a 143 kb genomic segment centered on the Waxy locus revealed patterns consistent with a selective sweep in the upstream region of Waxy. The selective sweep detected based on current limited genomic sequences exceeded over 50 kb, indicating strong selection in this or a bigger region. However, No sweep effect was detected in the repetitive downstream region of Waxy. Phylogenetic analysis indicated that Chinese waxy maize was domesticated from the cultivated flint maize (Zea mays ssp. mays) that was introduced from the new world. At least two independent deletions in exon 7 (30 bp) and 10 (15 bp) were identified in the Chinese accessions respectively. These findings demonstrate a similar pattern of domestication selection in the Waxy genomic region in both glutinous maize and rice, suggesting that this pattern in the rise of glutinous phenotype is likely in other cereal crops.     

Natural selection in gene-dense regions shapes the genomic pattern of polymorphism in wild and domesticated rice

Levels of nucleotide variability are frequently positively correlated with recombination rate and negatively associated with gene density due to the effects of selection on linked variation. These relationships are determined by properties that frequently differ among species, including the mating system, and aspects of genome organization such as how genes are distributed along chromosomes. In rice, genes are found at highest density in regions with frequent crossing-over. This association between gene density and recombination rate provides an opportunity to evaluate the effects of selection in a genomic context that differs from other model organisms. Using single-nucleotide polymorphism data from Asian domesticated rice Oryza sativa ssp. japonica and ssp. indica and their progenitor species O. rufipogon, we observe a significant negative association between levels of polymorphism and both gene and coding site density, but either no association, or a negative correlation, between nucleotide variability and recombination rate. We establish that these patterns are unlikely to be explained by neutral mutation rate biases and demonstrate that a model of background selection with variable rates of deleterious mutation is sufficient to account for the gene density effect in O. rufipogon. In O. sativa ssp. japonica, we report a strong negative correlation between polymorphism and recombination rate and greater losses of variation during domestication in the euchromatic chromosome arms than heterochromatin. This is consistent with Hill-Robertson interference in low-recombination regions, which may limit the efficacy of selection for domestication traits. Our results suggest that the physical distribution of selected mutations is a primary factor that determines the genomic pattern of polymorphism in wild and domesticated rice species.     

Global dissemination of a single mutation conferring white pericarp in rice

Here we report that the change from the red seeds of wild rice to the white seeds of cultivated rice (Oryza sativa) resulted from the strong selective sweep of a single mutation, a frame-shift deletion within the Rc gene that is found in 97.9% of white rice varieties today. A second mutation, also within Rc, is present in less than 3% of white accessions surveyed. Haplotype analysis revealed that the predominant mutation originated in the japonica subspecies and crossed both geographic and sterility barriers to move into the indica subspecies. A little less than one Mb of japonica DNA hitchhiked with the rc allele into most indica varieties, suggesting that other linked domestication alleles may have been transferred from japonica to indica along with white pericarp color. Our finding provides evidence of active cultural exchange among ancient farmers over the course of rice domestication coupled with very strong, positive selection for a single white allele in both subspecies of O. sativa.     

同源四倍体水稻籼型突变体Wx基因序列分析及特异识别分子标记的建立

同源四倍体水稻突变株D4063-1直链淀粉含量比来源二倍体明恢63下降一半,即其直链淀粉含量为5.23%。为研究其直链淀粉含量下降的原因, 根据普通水稻Wx基因设计引物, 扩增测序获得了D4063-1Wx基因的全序列, 并与已报道的Wx基因进行比对分析; 同源四倍体水稻D4063-1Wx基因最显著变化为在外显子序列中发生碱基缺失, 导致移码突变, 在第9外显子终止密码子提前出现。D4063-1Wx基因碱基位点的变化还导致其序列上酶切位点的变化,对常用限制性内切酶位点分析结果表明, 同源四倍体水稻相对于籼稻和粳稻多了2个sphⅠ酶切位点, 相对于粳稻减少了6个AccⅠ, 增加了4个XbaⅠ, 1个XhoⅠ, 1个PstⅠ和1个SalⅠ酶切位点。聚类分析表明D4063-1Wx基因序列与籼稻亲源关系较近, 由此推测D4063-1Wx基因来源于籼稻的Wxa基因型。另外, 根据D4063-1Wx基因的碱基差异, 推测D4063-1Wx基因外显子碱基变化导致的RNA加工障碍是其直链淀粉降低的主要原因, 并可能与其米饭较软等品质相关。本研究还根据D4063-1和籼稻、粳稻的序列差异及D4063-1在该片段上的特征序列位点设计了用于识别D4063-1的寡核苷酸片段,并作为PCR反应的引物命名为AUT4063-1,将该引物与作者设计的扩增普通籼稻、粳稻Wx基因的引物F5配合使用, 建立了识别D4063-1的显性和共显性两种检测方式的分子标记, 为快速、准确鉴别低直链淀粉含量突变体D4063-1创造了条件。     

Independent domestication of Asian rice followed by gene flow from japonica to indica

Results from studies on the domestication process of Asian rice Oryza sativa have been controversial because of its complicated evolutionary history. Previous studies have yielded two alternative hypotheses about the origin(s) of the two major groups of O. sativa: japonica and indica. One study proposes a single common wild ancestor, whereas the other suggests that there were multiple domestication events of different types of wild rice. Here, we provide clear evidence of the independent domestication of japonica and indica obtained via high-throughput sequencing and a large-scale comparative analysis of two wild rice accessions (W1943 and W0106) and two cultivars (a japonica cultivar called "Nipponbare" and an indica cultivar called "Guangluai-4"). The different domestication processes of the two cultivar groups appear to have led to distinct patterns of molecular evolution in protein-coding regions. The intensity of purifying selection was relaxed only in the japonica group, possibly because of a bottleneck effect. Moreover, a genome-wide comparison between Nipponbare, Guangluai-4, and another indica cultivar (93-11) suggests multiple hybridization events between japonica and indica, both before and after the divergence of the indica cultivars. We found that a large amount of genomic DNA, including domestication-related genes, was transferred from japonica to indica, which might have been important in the development of modern rice. Our study provides an overview of the dynamic process of Asian rice domestication, including independent domestication events and subsequent gene flow.     

Evidence of selection on fatty acid biosynthetic genes during the evolution of cultivated sunflower

The identification of genes underlying the phenotypic transitions that took place during crop evolution, as well as the genomic extent of resultant selective sweeps, is of great interest to both evolutionary biologists and applied plant scientists. In this study, we report the results of a molecular evolutionary analysis of 11 genes that underlie fatty acid biosynthesis and metabolism in wild and cultivated sunflower (Helianthus annuus). Seven of these 11 genes showed evidence of selection at the nucleotide level, with 1 (FAD7) having experienced selection prior to domestication, 2 (FAD2-3 and FAD3) having experienced selection during domestication, and 4 (FAB1, FAD2-1, FAD6, and FATB) having experienced selection during the subsequent period of improvement. Sequencing of a subset of these genes from an extended panel of sunflower cultivars revealed little additional variation, and an analysis of the genomic region surrounding one of these genes (FAD2-1) revealed the occurrence of an extensive selective sweep affecting a region spanning at least ca. 100 kb. Given that previous population genetic analyses have revealed a relatively rapid decay of linkage disequilibrium in sunflower, this finding indicates the occurrence of strong selection and a rapid sweep.     

The extent of linkage disequilibrium in rice (Oryza sativa L.)

Despite its status as one of the world's major crops, linkage disequilibrium (LD) patterns have not been systematically characterized across the genome of Asian rice (Oryza sativa). Such information is critical to fully exploit the genome sequence for mapping complex traits using association techniques. Here we characterize LD in five 500-kb regions of the rice genome in three major cultivated rice varieties (indica, tropical japonica, and temperate japonica) and in the wild ancestor of Asian rice, Oryza rufipogon. Using unlinked SNPs to determine the amount of background linkage disequilibrium in each population, we find that the extent of LD is greatest in temperate japonica (probably >500 kb), followed by tropical japonica (approximately 150 kb) and indica (approximately 75 kb). LD extends over a shorter distance in O. rufipogon (<40 kb) than in any of the O. sativa groups assayed here. The differences in the extent of LD among these groups are consistent with differences in outcrossing and recombination rate estimates. As well as heterogeneity between groups, our results suggest variation in LD patterns among genomic regions. We demonstrate the feasibility of genomewide association mapping in cultivated Asian rice using a modest number of SNPs.     

2012年第11期《遗传》封面说明

作物的驯化是人类从开始种植和储存的野生作物中选择优良性状,使之形态特征适应于农业生产方向进化的过程,因此,大部分种子作物驯化后在落粒性、种子休眠和植株形态等方面都出现了相似的变化。水稻是研究谷类作物驯化的良好模式生物。稻属包含2种栽培稻,分别为亚洲栽培稻(Oryza sativa L.)和非洲栽培稻(O. glaberrima Steud.),其中亚洲栽培稻遍布全世界,包含两个主要亚种,粳稻亚种(O. sativa L. ssp. japonica)和籼稻亚种(O. sativa L. ssp. indica)。稻属丰富的近缘种和广泛的地域分布非常有利于研究确定现代栽培稻的驯化地域。此外,水稻基因组较小、具高质量精细图谱,加上功能基因研究上的进展,也为深入开展水稻驯化进程研究奠定了基础。详见本期第XX-XX页区树俊,汪鸿儒,储成才“亚洲栽培稻主要驯化性状研究进展”,对水稻关键驯化性状研究进行的比较全面的综述。封面图中央是选取23株AA基因组的亚洲栽培稻及其近缘野生稻,利用水稻驯化过程中受到选择的控制稻壳颜色基因Bh4上下游各50 kb中的SNP位点所构建的进化树;图外从左下至右下沿顺时针方向,反映的是水稻驯化过程中稻壳颜色、谷粒形状、穗型的变化趋势。 区树俊,汪鸿儒,储成才（绘图：区树俊）     

The molecular basis of white pericarps in African domesticated rice: novel mutations at the Rc gene

Repeated phenotypic evolution can occur at both the inter- and intraspecific level and is especially prominent in domesticated plants, where artificial selection has favoured the same traits in many different species and varieties. The question of whether repeated evolution reflects changes at the same or different genes in each lineage can now be addressed using the domestication and improvement genes that have been identified in a variety of crops. Here, we document the genetic basis of nonpigmented ('white') pericarps in domesticated African rice (Oryza glaberrima) and compare it with the known genetic basis of the same trait in domesticated Asian rice (Oryza sativa). In some cases, white pericarps in African rice are apparently caused by unique mutations at the Rc gene, which also controls pericarp colour variation in Asian rice. In one case, white pericarps appear to reflect changes at a different gene or potentially a cis-regulatory region.     

The accumulation of deleterious mutations in rice genomes: a hypothesis on the cost of domestication

The extent of molecular differentiation between domesticated animals or plants and their wild relatives is postulated to be small. The availability of the complete genome sequences of two subspecies of the Asian rice, Oryza sativa (indica and japonica) and their wild relatives have provided an unprecedented opportunity to study divergence following domestication. We observed significantly more amino acid substitutions during rice domestication than can be expected from a comparison among wild species. This excess is disproportionately larger for the more radical kinds of amino acid changes (e.g. Cys<-->Tyr). We estimate that approximately a quarter of the amino acid differences between rice cultivars are deleterious, not accountable by the relaxation of selective constraints. This excess is negatively correlated with the rate of recombination, suggesting that 'hitchhiking' has occurred. We hypothesize that during domestication artificial selection increased the frequency of many deleterious mutations.     

Positive Selection in the Chromosome 16 VKORC1 Genomic Region Has Contributed to the Variability of Anticoagulant Response in Humans

VKORC1 (vitamin K epoxide reductase complex subunit 1, 16p11.2) is the main genetic determinant of human response to oral anticoagulants of antivitamin K type (AVK). This gene was recently suggested to be a putative target of positive selection in East Asian populations. In this study, we genotyped the HGDP-CEPH Panel for six VKORC1 SNPs and downloaded chromosome 16 genotypes from the HGDP-CEPH database in order to characterize the geographic distribution of footprints of positive selection within and around this locus. A unique VKORC1 haplotype carrying the promoter mutation associated with AVK sensitivity showed especially high frequencies in all the 17 HGDP-CEPH East Asian population samples. VKORC1 and 24 neighboring genes were found to lie in a 505 kb region of strong linkage disequilibrium in these populations. Patterns of allele frequency differentiation and haplotype structure suggest that this genomic region has been submitted to a near complete selective sweep in all East Asian populations and only in this geographic area. The most extreme scores of the different selection tests are found within a smaller 45 kb region that contains VKORC1 and three other genes (BCKDK, MYST1 (KAT8), and PRSS8) with different functions. Because of the strong linkage disequilibrium, it is not possible to determine if VKORC1 or one of the three other genes is the target of this strong positive selection that could explain present-day differences among human populations in AVK dose requirement. Our results show that the extended region surrounding a presumable single target of positive selection should be analyzed for genetic variation in a wide range of genetically diverse populations in order to account for other neighboring and confounding selective events and the hitchhiking effect.     

Introduction of Wx transgene into rice wx mutants leads to both high- and low-amylose rice

The Waxy (Wx) gene encodes a granule-bound starch synthase (GBSS) that plays a key role in the amylose synthesis of rice and other plant species. Two functional Wx alleles of rice exist: Wx(a), which produces a large amount of amylose, and Wx(b), which produces a smaller amount of amylose because of the mutation at the 5' splice site of intron 1. Wx(b) is largely distributed in Japonica cultivars, and high amylose cultivars do not exist in Japonica cultivars. We introduced the cloned Wx(a) cDNA into null-mutant Japonica rice (wx). The amylose contents of these transgenic plants were 6-11% higher than that of the original cultivar, Labelle, which carries the Wx(a) allele, although the levels of the Wx protein in the transgenic rice were equal to those of cv. Labelle. We also observed a gene-dosage effect of the Wx(a) transgene on Wx protein expression, but a smaller dosage effect was observed in amylose production with over 40% of amylose content in transgenic rice. Moreover, one transgenic line carrying eleven copies of the transgene showed low levels of Wx expression and amylose in the endosperm. This suggested that the integration of excessive copies of the transgene might lead to gene silencing.     

Molecular aspect of good eating quality formation in Japonica rice

The composition of amylopectin is the determinant of rice eating quality under certain threshold of protein content and the ratio of amylose and amylopectin. In molecular biology level, the fine structure of amylopectin is determined by relative activities of starch branching enzyme (SBE), granule-bound starch synthase (GBSS), and soluble starch synthase (SSS) in rice grain under the same ADP-Glucose level. But the underlying mechanism of eating quality in molecular biology level remains unclear. This paper reports the differences on major parameters such as SNP and insertion-deletion sites, RNA expressions, and enzyme activities associated with eating quality of japonica varieties. Eight japonica rice varieties with significant differences in various eating quality parameters such as palatability and protein content were used in this experiment. Association analysis between nucleotide polymorphism and eating quality showed that S12 and S13 loci in SBE1, S55 in SSS1, S58 in SSS2A were significantly associated with apparent amylose content, alkali digestion value, setback viscosity, consistency viscosity, pasting temperature, which explained most of the variation in apparent amylose content, setback viscosity, and consistency viscosity; and explained almost all variations in alkali digestion value and pasting temperature. Thirty-five SNPs and insertion-deletions from SBE1, SBE3, GBSS1, SSS1, and SSS2A differentiated high or intermediate palatability rice varieties from low palatability rice varieties. Correlation analysis between enzyme activities and eating quality properties revealed that SBE25 and SSS15/W15 were positively correlated with palatability, whereas GBSS10 and GBSS15 were negatively correlated. Gene expressions showed that SBE1 and SBE3 expressions in high palatability varieties tended to be higher than middle and low palatability varieties. Collectively, SBE1, SBE3, SSS1, and SSS2A, especially SBE1 and SBE3 could improve eating quality, but GBSS1 decreased eating quality. The results indicated the possibility of developing high palatability cultivars through modification of key genes related to japonica rice eating quality formation in starch biosynthesis.     

A genome-wide scan for signatures of directional selection in domesticated pigs

Background

Animal domestication involved drastic phenotypic changes driven by strong artificial selection and also resulted in new populations of breeds, established by humans. This study aims to identify genes that show evidence of recent artificial selection during pig domestication.

Results

Whole-genome resequencing of 30 individual pigs from domesticated breeds, Landrace and Yorkshire, and 10 Asian wild boars at ~16-fold coverage was performed resulting in over 4.3 million SNPs for 19,990 genes. We constructed a comprehensive genome map of directional selection by detecting selective sweeps using an F_ST-based approach that detects directional selection in lineages leading to the domesticated breeds and using a haplotype-based test that detects ongoing selective sweeps within the breeds. We show that candidate genes under selection are significantly enriched for loci implicated in quantitative traits important to pig reproduction and production. The candidate gene with the strongest signals of directional selection belongs to group III of the metabolomics glutamate receptors, known to affect brain functions associated with eating behavior, suggesting that loci under strong selection include loci involved in behaviorial traits in domesticated pigs including tameness.

Conclusions

We show that a significant proportion of selection signatures coincide with loci that were previously inferred to affect phenotypic variation in pigs. We further identify functional enrichment related to behavior, such as signal transduction and neuronal activities, for those targets of selection during domestication in pigs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1330-x) contains supplementary material, which is available to authorized users.