Search

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 (~150 kb) and indica (~75 kb). LD extends over a shorter distance in O. rufipogon ([much less than] 440 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.

Barley (Hordeum vulgare ssp. vulgare) was first cultivated 10,500 years ago in the Fertile Crescent and is one of the founder crops of Eurasian agriculture. Phylogeographic analysis of five nuclear loci and morphological assessment of two traits in >250 domesticated barley accessions reveal that landraces found in South and East Asia are genetically distinct from those in Europe and North Africa. A Bayesian population structure assessment method indicates that barley accessions are subdivided into six clusters and that barley landraces from 10 different geographical regions of Eurasia and North Africa show distinct patterns of distribution across these clusters. Using haplotype frequency data, it appears that the Europe/ North Africa landraces are most similar to the Near East population ([F.sub.ST] = 0.15) as well as to wild barley ([F.sub.ST] = 0.11) and are strongly differentiated from all other Asian populations ([F.sub.ST] = 0.34-0.74). A neighbor-joining analysis using these [F.sub.ST] estimates also supports a division between European, North African, and Near East barley types from more easterly Asian accessions. There is also differentiation in the presence of a naked caryopsis and spikelet row number between eastern and western barley accessions. The data support the differential migration of barley from two domestication events that led to the origin of barley--one in the Fertile Crescent and another farther east, possibly at the eastern edge of the Iranian Plateau--with European and North African barley largely originating from the former and much of Asian barley arising from the latter. This suggests that cultural diffusion or independent innovation is responsible for the expansion of agriculture to areas of South and East Asia during the Neolithic revolution.

Association mapping focused on 36 genes involved in branch development was used to identify candidate genes for variation in shoot branching in Arabidopsis thaliana. The associations between four branching traits and moderate-frequency haplogroups at the studied genes were tested in a panel of 96 accessions from a restricted geographic range in Central Europe. Using a mixed-model association-mapping method, we identified three loci--MORE AXILLARY GROWTH 2 (MAX2), MORE AXILLARY GROWFH 3 (MAX3), and SUPERSHOOT 1 (SPS1)--that were significantly associated with branching variation. On the basis of a more extensive examination of the MAX2 and MAX3 genomic regions, we find that linkage disequilibrium in these regions decays within ~10 kb and trait associations localize to the candidate genes in these regions. When the significant associations are compared to relevant quantitative trait loci (QTL) from previous Ler x Col and Cvi x Ler recombinant inbred line (RIL) mapping studies, no additive QTL overlapping these candidate genes are observed, although epistatic QTL for branching, including one that spans the SPS1, are found. These results suggest that epistasis is prevalent in determining branching variation in A. thaliana and may need to be considered in linkage disequilibrium mapping studies of genetically diverse accessions.

Rice (Oryza sativa) was cultivated by Asian Neolithic farmers > 11,000 years ago, and different cultures have selected for divergent starch qualities in the rice grain during and after the domestication process. An intron 1 splice donor site mutation of the Waxy gene is responsible for the absence of amylose in glutinous rice varieties. This mutation appears to have also played an important role in the origin of low amylose, nonglutinous temperate japonica rice varieties, which form a primary component of Northeast Asian cuisines. Waxy DNA sequence analyses indicate that the splice donor mutation is prevalent in temperate japonica rice varieties, but rare or absent in tropical japonica, indica, aus, and aromatic varieties. Sequence analysis across a 500-kb genomic region centered on Wax), reveals patterns consistent with a selective sweep in the temperate japonicas associated with the mutation. The size of the selective sweep (>250 kb) indicates very strong selection in this region, with an inferred selection coefficient that is higher than similar estimates from maize domestication genes or wild species. These findings demonstrate that selection pressures associated with crop domestication regimes can exceed by one to two orders of magnitude those observed for genes under even strong selection in natural systems.

High-diversity genes represent an important class of loci in organismal genomes. Since elevated levels of nucleotide variation are a key component of the molecular signature for balancing selection or local adaptation, high-diversity genes may represent loci whose alleles are selectively maintained as balanced polymorphisms. Comparison of 4300 random shotgun sequence fragments of the Arabidopsis thaliana L er ecotype genome with the whole genomic sequence of the Col-0 ecotype identified 60 genes with putatively high levels of intraspecific variability. Eleven of these genes were sequenced in multiple A. thaliana accessions, 3 of which were found to display elevated levels of nucleotide polymorphism. These genes encode the myb-like transcription factor MYB103, a putative soluble starch synthase I, and a homeodomainleucine zipper transcription factor. Analysis of these genes and 4-7 flanking genes in 14-20 A. thaliana ecotypes revealed that two of these loci show other characteristics of balanced polymorphisms, including broad peaks of nucleotide diversity spanning multiple linked genes and an excess of intermediate-frequency polymorphisms. Scanning genomes for high-diversity genomic regions may be useful in approaches to adaptive trait locus mapping for uncovering candidate balanced polymorphisms.

The selfing plant Arabidopsis thaliana has been proposed to be well suited for linkage disequilibrium (LD) mapping as a means of identifying genes underlying natural trait variation. Here we apply LD mapping to examine haplotype variation in the genomic region of the photoperiod receptor CRYPTOCHROME2 and associated flowering time variation. CRY2DNA sequences reveal strong LD and the existence of two highly differentiated haplogroups (A and B) across the gene; in addition, a haplotype possessing a radical glutamineto-serine replacement ([A.sup.s]) occurs within the more common haplogroup. Growth chamber and field experiments using an unstratified population of 95 ecotypes indicate that tinder short-day photoperiod, the [A.sup.s] and B haplogroups are both highly significantly associated with early flowering. Data from six genes flanking CRY2 indicate that these haplogroups are limited to an ~65-kb genomic region around CRY2. Whereas the B haplogroup cannot be delimited to < 16 kb around CRY2, the [A.sup.s] haplogroup is characterized almost exclusively by the nucleotide polymorphisms directly associated with the serine replacement in CRY2; this finding strongly suggests that the serine substitution is directly responsible for the [A.sup.s] early flowering phenotype. This study demonstrates the utility of LD mapping for elucidating the genetic basis of natural, ecologically relevant variation in Arabidopsis.

Phenotypic plasticity and genotype-environment interactions (GEI) play a prominent role in plant morphological diversity and in the potential functional capacities of plant life-history traits. The genetic basis of plasticity and GEI, however, is poorly understood in most organisms. In this report, inflorescence development patterns in Arabidopsis thaliana were examined under different, ecologically relevant photoperiod environments for two recombinant inbred mapping populations (Ler X Col and Cvi X Ler) using a combination of quantitative genetics and quantitative trait locus (QTL) mapping. Plasticity and GEI were regularly observed for the majority of 13 inflorescence traits. These observations can be attributable (at least partly) to variable effects of specific QTL. Pooled across traits, 12/44 (27.3%) and 32/62 (51.6%) of QTL exhibited significant QTL X environment interactions in the Ler X Col and Cvi X Ler lines, respectively. These interactions were attributable to changes in magnitude of effect of QTL more often than to changes in rank order (sign) of effect. Multiple QTL X environment interactions (in Cvi X Ler) clustered in two genomic regions on chromosomes 1 and 5, indicating a disproportionate contribution of these regions to the phenotypic patterns observed. High-resolution mapping will be necessary to distinguish between the alternative explanations of pleiotropy and tight linkage among multiple genes.

Genetic variation for quantitative traits is often greater than that expected to be maintained by mutation in the face of purifying natural selection. One possible explanation for this observed variation is the action of heterogeneous natural selection in the wild. Here we report that selection on quantitative trait loci (QTL) for fitness traits in the model plant species Arabidopsis thaliana differs among natural ecological settings and genetic backgrounds. At one QTL, the allele that enhanced the viability of fall-germinating seedlings in North Carolina reduced the fecundity of spring-germinating seedlings in Rhode Island. Several other QTL experienced strong directional selection, but only in one site and seasonal cohort. Thus, different loci were exposed to selection in different natural environments. Selection on allelic variation also depended upon the genetic background. The allelic fitness effects of two QTL reversed direction depending on the genotype at the other locus. Moreover, alternative alleles at each of these loci caused reversals in the allelic fitness effects of a QTL closely linked to TFL1, a candidate developmental gene displaying nucleotide sequence polymorphism consistent with balancing selection. Thus, both environmental heterogeneity and epistatic selection may maintain genetic variation for fitness in wild plant species.

The Arabidopsis thaliana CLAVATA2 (CLV2) gene encodes a leucine-rich repeat protein that regulates the development of the shoot meristem. The levels and patterns of nucleotide variation were assessed for CLV2 and 10 flanking genes that together span a 40-kb region of chromosome I. A total of 296 out of 7959 sequenced nucleotide sites were polymorphic. The mean levels of sequence diversity of the contiguous genes in this region are approximately twofold higher than those of other typical Arabidopsis nuclear loci. There is, however, wide variation in the levels and patterns of sequence variation among the 11 linked genes in this region, and adjacent genes appear to be subject to contrasting evolutionary forces. CLV2 has the highest levels of nucleotide variation in this region, a significant excess of intermediate frequency polymorphisms, and significant levels of intragenic linkage disequilibrium. Most alleles at CLV2 are found in one of three haplotype groups of moderate (>15%) frequency. These features suggest that CLV2 may harbor a balanced polymorphism.

Genes that have undergone positive or diversifying selection are likely to be associated with adaptive divergence between species. One indicator of adaptive selection at the molecular level is an excess of amino acid replacement fixed differences per replacement site relative to the number of synonymous fixed differences per synonymous site (to [omega] = [K.sub.a]/[K.sub.s]). We used an evolutionary expressed sequence tag (EST) approach to estimate the distribution of [omega] among 304 orthologous loci between Arabidopsis thaliana and A. lyrata to identify genes potentially involved in the adaptive divergence between these two Brassicaceae species. We find that 14 of 304 genes (~5%) have an estimated [omega] > 1 and are candidates for genes with increased selection intensities. Molecular population genetic analyses of 6 of these rapidly evolving protein loci indicate that, despite their high levels of between-species nonsynonymous divergence, these genes do not have elevated levels of intraspecific replacement polymorphisms compared to previously studied genes. A hierarchical Bayesian analysis of protein-coding region evolution within and between species also indicates that the selection intensities of these genes are elevated compared to previously studied A. thaliana nuclear loci.

Molecular biologists are rapidly characterizing the genetic basis of flowering in model species such as Arabidopsis thaliana. However, it is not clear how the developmental pathways identified in controlled environments contribute to variation in reproductive timing in natural ecological settings. Here we report the first study of quantitative trait loci (QTL) for date of bolting (the transition from vegetative to reproductive growth) in A. thaliana in natural seasonal field environments and compare the results with those obtained under typical growth-chamber conditions. Two QTL specific to long days in the chamber were expressed only in spring-germinating cohorts in the field, and two loci specific to short days in the chamber were expressed only in fall-germinating cohorts, suggesting differential involvement of the photoperiod pathway in different seasonal environments. However, several other photoperiod-specific QTL with large effects in controlled conditions were undetectable in natural environments, indicating that expression of allelic variation at these loci was overridden by environmental factors specific to the field. Moreover, a substantial number of QTL with major effects on bolting date in one or more field environments were undetectable under controlled environment conditions. These novel loci suggest the involvement of additional genes in the transition to flowering under ecologically relevant conditions.

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.

The floral developmental pathway in Arabidopsis thaliana is composed of several interacting regulatory genes, including the inflorescence architecture gene TERMINAL FLOWER1 (TFL1), the floral meristem identity genes LEAFY (LFY), APETALA1 (AP1), and CAULIFLOWER (CAL), and the floral organ identity genes APETALA3 (AP3) and PISTILLATA (PI). Molecular population genetic analyses of these different genes indicate that the coding regions of AP3 and PI, as well as AP1 and CAL, share similar levels and patterns of nucleotide diversity. In contrast, the coding regions of TFL1 and LFY display a significant reduction in nucleotide variation, suggesting that these sequences have been subjected to a recent adaptive sweep. Moreover, the promoter of TFL1, unlike its coding region, displays high levels of diversity organized into two distinct haplogroups that appear to be maintained by selection. These results suggest that patterns of molecular evoution differ among regulatory genes in this developmental pathway, with the earlier acting genes exhibiting evidence of adaptive evolution.

Variation in inflorescence development patterns is a central factor in the evolutionary ecology of plants. The genetic architectures of 13 traits associated with inflorescence developmental timing, architecture, rosette morphology, and fitness were investigated in Arabidopsis thaliana, a model plant system. There is substantial naturally occurring genetic variation for inflorescence development traits, with broad sense heritabilities computed from 21 Arabidopsis ecotypes ranging from 0.134 to 0.772. Genetic correlations are significant for most (64/78) pairs of traits, suggesting either pleiotropy or tight linkage among loci. Quantitative trait locus (QTL) mapping indicates 47 and 63 QTL for inflorescence developmental traits in Ler x Col and Cvi x Ler recombinant inbred mapping populations, respectively. Several QTL associated with different developmental traits map to the same Arabidopsis chromosomal regions, in agreement with the strong genetic correlations observed. Epistasis among QTL was observed only in the Cvi x Ler population, and only between regions on chromosomes 1 and 5. Examination of the completed Arabidopsis genome sequence in three QTL regions revealed between 375 and 783 genes per region. Previously identified flowering time, inflorescence architecture, floral meristem identity, and hormone signaling genes represent some of the many candidate genes in these regions.

The evolution of plant morphologies during domestication events provides clues to the origin of crop species and the evolutionary genetics of structural diversification. The CAULIFLOWER gene, a floral regulatory locus, has been implicated in the cauliflower phenotype in both Arabidopsis thaliana and Brassica oleracea. Molecular population genetic analysis indicates that alleles carrying a nonsense mutation in exon 5 of the B. oleracea CAULIFLOWER (BoCAL) gene are segregating in both wild and domesticated B. oleracea subspecies. Alleles carrying this nonsense mutation are nearly fixed in B. oleracea ssp. botrytis (domestic cauliflower) and B. oleracea ssp. italica (broccoli), both of which show evolutionary modifications of inflorescence structures. Tests for selection indicate that the pattern of variation at this locus is consistent with positive selection at BoCAL in these two subspecies. This nonsense polymorphism, however, is also present in both B. oleracea ssp. acephala (kale) and B. oleracea ssp. oleracea (wild cabbage). These results indicate that specific alleles of BoCAL were selected by early farmers during the domestication of modified inflorescence structures in B. oleracea.

Investigates the molecular evolution of the myc-like R regulatory gene family in seven plant species. Conservation between monocot and dicot R genes; Mean replacement rates for the conserved regions; Analysis of R homologs within the grasses; Factors contributing to the diversification of the variable regions within the regulatory genes.