Highly structured nucleotide variation within and among Arabidopsis lyrata populations at the FAH1 and DFR gene regions

Nucleotide variation at the FAH1 and DFR gene regions was surveyed in four populations of Arabidopsis lyrata (two European A. l. petraea and two North American A. l. lyrata populations). In contrast to previous results, levels of variation were not consistently lower in A. l. lyrata than in A. l. petraea, and similar degrees of genetic differentiation were detected between and within subspecies. These observations and the significant genetic differentiation detected among populations suggest population substructure and no real subdivision between subspecies. For each gene studied, genotypic data were obtained, which allowed comparing nucleotide diversity within individuals (between sequences from the same individual) and within populations (between sequences from the same population). The generally lower level of variation within than among individuals detected in each population yielded a significant deviation from panmixia within populations. In three of the four populations studied, two highly divergent alleles were detected within populations at the highly variable DFR locus. This pattern and the significant excess of derived variants detected in most populations suggest that most variation segregating within populations results from rare migration events between relatively small and isolated populations exhibiting reduced panmixia.     

Population genetic structure of Arabidopsis lyrata in Europe

Population genetic theory predicts that the self-incompatible and perennial herb, Arabidopsis lyrata, will have a genetic structure that differs from the self-fertilizing, annual Arabidopsis thaliana. We quantified the genetic structure for eight populations of A. lyrata ssp. petraea in historically nonglaciated regions of central Europe. Analysis of 20 microsatellite loci for 344 individuals demonstrated that, in accordance with predictions, diploid populations had high genome-wide heterozygosity (H(O) = 0.48; H(E) = 0.52), high within-population diversity (83% of total) compatible with mutation-drift equilibrium, and moderate differentiation among populations (F(ST) = 0.17). Within a single population, the vast majority of genetic variability (92%) was found at the smallest spatial scale (< 3 m). Although there was no evidence of biparental inbreeding or clonal propagation at this scale (F(IS) = 0.003), significant fine-scale spatial autocorrelation indicated localized gene flow presumably due to gravity dispersed seeds (Sp = 0.018). Limited gene flow between isolated population clusters (regions) separated by hundreds of kilometres has given rise to an isolation by distance pattern of diversification, with low, but significant, differentiation among regions (F(ST) = 0.05). The maintenance of geographically widespread polymorphisms and uniformly high diversity throughout central Europe is consistent with periglacial survival of A. lyrata ssp. petraea north of the Alps in steppe-tundra habitats during the last glacial maximum. As expected of northern and previously glaciated localities, A. lyrata in Iceland was genetically less diverse and highly differentiated from central Europe (H(E) = 0.37; F(ST) = 0.27).     

Extremely low genetic variability and highly structured local populations of Arabidopsis thaliana at higher latitudes

The genetic diversity and population structure of Arabidopsis thaliana populations from Norway were studied and compared to a worldwide sample of A. thaliana to investigate the demographic history and elucidate possible colonization routes of populations at the northernmost species limit. We genotyped 282 individuals from 31 local populations using 149 single nucleotide polymorphism markers. A high level of population subdivision (F(ST) = 0.85 ± 0.007) was found indicating that A. thaliana is highly structured at the regional level. Significant relationships between genetic and geographical distances were found, suggesting an isolation by distance mode of evolution. Genetic diversity was much lower, and the level of linkage disequilibrium was higher in populations from the north (65-68°N) compared to populations from the south (59-62°N); this is consistent with a northward expansion pattern. A neighbour-joining tree showed that populations from northern Norway form a separate cluster, while the remaining populations are distributed over a few minor clusters. Minimal gene flow seems to have occurred between populations in different regions, especially between the geographically distant northern and southern populations. Our data suggest that northern populations represent a homogenous group that may have been established from a few founders during northward expansions, while populations in the central part of Norway constitute an admixed group established by founders of different origins, most probably as a result of human-mediated gene flow. Moreover, Norwegian populations appeared to be homogenous and isolated compared to a worldwide sample of A. thaliana, but they are still grouped with Swedish populations, which may indicate common colonization histories.     

Patterns of genetic diversity in outcrossing and selfing populations of Arabidopsis lyrata

Arabidopsis lyrata is normally considered an obligately outcrossing species with a strong self-incompatibility system, but a shift in mating system towards inbreeding has been found in some North American populations (subspecies A. lyrata ssp. lyrata). This study provides a survey of the Great Lakes region of Canada to determine the extent of this mating system variation and how outcrossing rates are related to current population density, geographical distribution, and genetic diversity. Based on variation at microsatellite markers (progeny arrays to estimate multilocus outcrossing rates and population samples to estimate diversity measures) and controlled greenhouse pollinations, populations can be divided into two groups: (i) group A, consisting of individuals capable of setting selfed seed (including autogamous fruit set in the absence of pollinators), showing depressed outcrossing rates (T(m) = 0.2-0.6), heterozygosity (H(O) = 0.02-0.06) and genetic diversity (H(E) = 0.08-0.10); and (ii) group B, consisting of individuals that are predominantly self-incompatible (T(m) > 0.8), require pollinators for seeds set, and showing higher levels of heterozygosity (H(O) = 0.13-0.31) and diversity (H(E) = 0.19-0.410). Current population density is not related to the shift in mating system but does vary with latitude. Restricted gene flow among populations was evident among all but two populations (F(ST) = 0.11-0.8). Group A populations were more differentiated from one another (F(ST) = 0.78) than they were from group B populations (F(ST) = 0.59), with 41% of the variation partitioned within populations, 47% between populations, and 12% between groups. No significant relationship was found between genetic and geographical distance. Results are discussed in the context of possible postglacial expansion scenarios in relation to loss of self-incompatibility.     

Effects of gene expression on molecular evolution in Arabidopsis thaliana and Arabidopsis lyrata

We analyzed the complete genome sequence of Arabidopsis thaliana and sequence data from 83 genes in the outcrossing A. lyrata, to better understand the role of gene expression on the strength of natural selection on synonymous and replacement sites in Arabidopsis. From data on tRNA gene abundance, we find a good concordance between codon preferences and the relative abundance of isoaccepting tRNAs in the complete A. thaliana genome, consistent with models of translational selection. Both EST-based and new quantitative measures of gene expression (MPSS) suggest that codon preferences derived from information on tRNA abundance are more strongly associated with gene expression than those obtained from multivariate analysis, which provides further support for the hypothesis that codon bias in Arabidopsis is under selection mediated by tRNA abundance. Consistent with previous results, analysis of protein evolution reveals a significant correlation between gene expression level and amino acid substitution rate. Analysis by MPSS estimates of gene expression suggests that this effect is primarily the result of a correlation between the number of tissues in which a gene is expressed and the rate of amino acid substitution, which indicates that the degree of tissue specialization may be an important determinant of the rate of protein evolution in Arabidopsis.     

Evidence for a large-scale population structure among accessions of Arabidopsis thaliana: possible causes and consequences for the distribution of linkage disequilibrium

The existence of a large-scale population structure was investigated in Arabidopsis thaliana by studying patterns of polymorphism in a set of 71 European accessions. We used sequence polymorphism surveyed in 10 fragments of approximately 600 nucleotides and a set of nine microsatellite markers. Population structure was investigated using a model-based inference framework. Among the accessions studied, the presence of four groups was inferred using genetic data, without using prior information on the geographical origin of the accessions. Significant genetic isolation by geographical distance was detected at the group level, together with a geographical gradient in allelic richness across groups. These results are discussed with respect to the previously proposed scenario of postglacial colonization of Europe from putative glacial refugia. Finally, the contribution of the inferred structure to linkage disequilibrium among 171 pairs of essentially unlinked markers was also investigated. Linkage disequilibrium analysis revealed that significant associations detected in the whole sample were mainly due to genetic differentiation among the inferred groups. We discuss the implication of this finding for future association studies in A. thaliana.     

Gene, phenotype and function: GLABROUS1 and resistance to herbivory in natural populations of Arabidopsis lyrata

The molecular genetic basis of adaptive variation is of fundamental importance for evolutionary dynamics, but is still poorly known. Only in very few cases has the relationship between genetic variation at the molecular level, phenotype and function been established in natural populations. We examined the functional significance and genetic basis of a polymorphism in production of leaf hairs, trichomes, in the perennial herb Arabidopsis lyrata. Earlier studies suggested that trichome production is subject to divergent selection. Here we show that the production of trichomes is correlated with reduced damage from insect herbivores in natural populations, and using statistical methods developed for medical genetics we document an association between loss of trichome production and mutations in the regulatory gene GLABROUS1. Sequence data suggest that independent mutations in this regulatory gene have provided the basis for parallel evolution of reduced resistance to insect herbivores in different populations of A. lyrata and in the closely related Arabidopsis thaliana. The results show that candidate genes identified in model organisms provide a valuable starting point for analysis of the genetic basis of phenotypic variation in natural populations.     

新疆北部拟南芥自然居群表型变异与协变

拟南芥(Arabidopsis thaliana)自然居群的表型特征代表其在自然环境下的适应状况, 不同居群间特征的对比可以为了解拟南芥表型变化规律, 进而分析其形成过程和机制提供重要线索。本研究以分布于新疆北部天山、塔尔巴哈台山和阿尔泰山的10个种群的9个表型性状为基础, 对比分析了小尺度、局域尺度和区域尺度环境下原生境拟南芥种群表型性状的变化。结果发现, 不同性状对环境变化的反应不同, 其中株高、株重、根重、根长、单个果实重、果实开裂力度在3种环境尺度下种群间的差异均达到极显著水平, 而分枝数、果实长度的种群间变化不显著, 种群间的表型分化系数较低。不同环境尺度下株重、根重、单株果数均表现出一致的协变格局, 反映了生理功能性状之间整合对拟南芥适应环境的重要性。同时, 各种群间整体的性状协变差异性明显, 根长、单个果实重、分枝数、果实长度、果实开裂力度等特征与其他特征协变具有明显的局部性, 局域尺度和区域尺度环境之间的变化较大。聚类分析发现区域尺度上的不同种群聚合在一起的现象非常突出, 进一步表明拟南芥的表型特征受微环境的强烈影响。Mantel检验表明, 小尺度上10个种群株高、株重、根重、单个果实重、果实长度、果实开裂力度6个性状变化存在显著的空间相关性, 而分枝数、根长的相关性却不显著。因此, 我们认为拟南芥表型变化受小尺度环境的影响强烈, 但在表型层面并非所有性状都与原生境气候存在遗传关联。     

Cross-species microsatellite markers for elucidating population genetic structure in Arabidopsis and Arabis (Brassicaeae)

Species closely related to model organisms present the opportunity to efficiently apply molecular and functional tools developed by a large research community to taxa with different ecological and evolutionary histories. We complied 42 microsatellite loci that amplify under common conditions in four closely related Arabidopsis: A. thaliana; A. halleri; A. lyrata ssp. lyrata; and A. lyrata ssp. petraea, as well as in one more distantly related crucifer; Arabis drummondii. Variation at these loci is amenable to a diversity of applications including population genetics, phylogeographical analyses, mapping of inter and intraspecific crosses, and recombination mapping. Our analysis of microsatellite variation illustrates significant differences in population genetic parameters among three Arabidopsis species. A population of A. thaliana, an inbreeding annual plant associated with disturbed habitats, was highly monomorphic (P = 8% percent polymorphic loci) and only 0.2% heterozygous for 648 locus-by-individual combinations. A population of the self-incompatible perennial herb, A. halleri, was more genetically variable (P = 71%) and had an excess of heterozygosity that may reflect a recent population bottleneck associated with human-mediated founder events. A population of the self-incompatible perennial herb, A. lyrata ssp. petraea, was even more genetically variable (P = 86%) and appeared to be at mutation-drift equilibrium. Population structure estimated from neutrally evolving loci provides an empirical expectation against which hypotheses of adaptive evolution at functional loci can be tested.     

拟南芥抗旱转录因子CBF4基因区域的核苷酸多样性及其分子进化分析

以生长于不同气候条件下的17个拟南芥核心生态型为材料,分析了它们的抗旱转录因子CBF4基因区域的序列多态性。结果表明：拟南芥CBF4基因区域具有高密度的单核苷酸多态性(SNP)和插入缺失(Indel),多态性频率为每35．8bp一个SNP,每143bp一个Indel,基因非编码区的多态性是编码区的4倍;在编码区,SNP的频率为每96．4bp一个SNP,其中发现25av、203av和244av 3个生态型CBF4基因区域1034位(以Gen—Bank登录号AB015478序列第19696位的核苷酸为1)碱基变化：G←→T,引起第205位氨基酸变化：gly←→val。核苷酸多样性统计分析显示,该基因内部大范围内存在连锁不平衡(linkage disequilibrium,LD),5’端非编码区有一个重组。与拟南芥等的研究结果类似,选择压力对不同的区域作用不同。3’端非编码区核苷酸多样性程度最高,是平衡性选择的结果,编码区核苷酸变化符合中性突变假说,而5’端非编码区是自然选择作用的靶位点。     

天山北部拟南芥生存群落特征及其与环境的关系

为了解拟南芥在天山北部的分布状况及环境依赖特点, 分析拟南芥的自然选择特征, 本文对天山北部分布的13个拟南芥(Arabidopsis thaliana)生存群落结构、组成及其与环境关系进行了研究, 并分析了拟南芥与群落主要物种的种间联结性。结果表明: 拟南芥生存的群落结构简单, 其中天山北坡中段的石河子、一四三团、沙湾、独山子地区的8个群落均为草本类型, 优势种相似, 而与伊犁果子沟、额敏和阿勒泰的5个群落差别较大。属的区系成分分析表明世界分布、北温带分布以及地中海、西亚至中亚分布型成分占大多数, 具有典型的地中海旱生植物区系分布特征, 体现了本地拟南芥分布及演化的干旱、半干旱的地理环境特点。采用双向指示种分析(TWINSPAN)将13个群落分为新疆绢蒿–猪毛菜–角果毛茛(Seriphidium kaschgaricum–Salsola collina–Ceratocephalus testiculatus)、新疆绢蒿–猪毛菜(S. kaschgaricum–S. collina)、新疆绢蒿–狭果鹤虱(S. kaschgaricum–Lappula semiglabra)、新疆绢蒿–旱麦草(S. kaschgaricum–Eremopyrum triticeum)、勿忘草–草原苔草(Myosotis sylvatica–Carex liparocarpos)5个群落类型。去势典范对应分析(DCCA)表明纬度、坡向、土壤有机质及pH值是决定天山北部拟南芥种群分布的主导因子。拟南芥分布与群落内其他物种有极强的依赖关系, 与13个群落62个主要物种的种间联结性分析表明, 共有119个正关联性种对, 明显高于72个负关联性种对, 与各群落优势种呈显著正关联。拟南芥种群分布数量在群落间差异较大, 分布于降雨较少的天山中部浅山地带拟南芥种群数量均高于降雨较丰富的天山西部伊犁果子沟地区, 是否发生适应性分化需要深入研究。     

Genetic isolation by distance in Arabidopsis thaliana: biogeography and postglacial colonization of Europe

Arabidopsis thaliana provides a useful model system for functional, evolutionary and ecological studies in plant biology. We have analysed natural genetic variation in A. thaliana in order to infer its biogeographical and historical distribution across Eurasia. We analysed 79 amplified fragment length polymorphism (AFLP) markers in 142 accessions from the species' native range, and found highly significant genetic isolation by distance among A. thaliana accessions from Eurasia and southern Europe. These spatial patterns of genetic variation suggest that A. thaliana colonized central and northern Europe from Asia and from Mediterranean Pleistocene refugia, a trend which has been identified in other species. Statistically significant levels of multilocus linkage disequilibrium suggest intermediate levels of disequilibrium among subsets of loci, and analysis of genetic relationships among accessions reveal a star or bush-like dendrogram with low bootstrap support. Taken together, it appears that there has been sufficient historical recombination in the A. thaliana genome such that accessions do not conform to a tree-like, bifurcating pattern of evolution - there is no 'ecotype phylogeny.' Nonetheless, significant isolation by distance provides a framework upon which studies of natural variation in A. thaliana may be designed and interpreted.     

Uneven segregation of sporophytic self-incompatibility alleles in Arabidopsis lyrata

Self-incompatibility in Arabidopsis lyrata is sporophytically controlled by the multi-allelic S-locus. Self-incompatibility alleles (S-alleles) are under strong negative frequency dependent selection because pollen carrying common S-alleles have fewer mating opportunities. Population genetics theory predicts that deleterious alleles can accumulate if linked to the S-locus. This was tested by studying segregation of S-alleles in 11 large full sib families in A. lyrata. Significant segregation distortion leading to an up to fourfold difference in transmission rates was found in six families. Differences in transmission rates were not significantly different in reciprocal crosses and the distortions observed were compatible with selection acting at the gametic stage alone. The S-allele with the largest segregation advantage is also the most recessive, and is very common in natural populations concordant with its apparent segregation advantage. These results imply that frequencies of S-alleles in populations of A. lyrata cannot be predicted based on simple models of frequency-dependent selection alone.     

Conditional neutrality at two adjacent NBS-LRR disease resistance loci in natural populations of Arabidopsis lyrata

We examined patterns of nucleotide diversity at a genomic region containing two linked candidate disease resistance (NBS-LRR) genes in seven populations of the outcrossing plant Arabidopsis lyrata. In comparison with two adjacent control genes and neutral reference genes across the genome, the NBS-LRR genes exhibited elevated nonsynonymous variation and a large number of major-effect polymorphisms causing early stop codons and/or frameshift mutations. In contrast, analysis of synonymous diversity provided no evidence that the region was subject to long-term balancing selection or recent selective sweeps in any of the seven populations surveyed. Also in contrast with earlier surveys of one of these R genes, there was no evidence that the resistance genes or the major-effect mutations were subject to elevated differentiation between populations. We suggest that conditional neutrality in the absence of the corresponding pathogen, rather than long-term balancing selection or local adaptation, may in some circumstances be a significant cause of elevated functional polymorphism at R genes. In contrast with the R genes, analysis of diversity and differentiation at the flanking FERONIA locus showed high population divergence, suggesting local adaptation on this locus controlling male-female signalling during fertilization.     

Development of EST‐derived microsatellite markers for Arabidopsis lyrata subspecies petraea (L.)

A set of expressed sequence tag–simple sequence repeat (EST‐SSR) loci has been developed for Arabidopsis lyrata ssp. petraea. From 768 root cDNA clones, 126 microsatellites, including di‐, tri‐, tetra‐ and pentanucleotide repeat motifs were identified and primers were designed to 24 EST‐SSRs. Eleven loci were subsequently screened on 150 individuals sampled from five natural populations, which revealed three to nine alleles per locus (mean 5.36) and expected heterozygosity (H_E) estimates ranging from 0.046 to 0.698. Significant deviations from random mating were observed at 10 EST‐SSR loci, likely due to inbreeding (global F_IS = 0.151) and population structure (global F_ST = 0.246).     

Genome-wide population structure and evolutionary history of the Frizarta dairy sheep

In the present study, we used genomic data, generated with a medium density single nucleotide polymorphisms (SNP) array, to acquire more information on the population structure and evolutionary history of the synthetic Frizarta dairy sheep. First, two typical measures of linkage disequilibrium (LD) were estimated at various physical distances that were then used to make inferences on the effective population size at key past time points. Population structure was also assessed by both multidimensional scaling analysis and k-means clustering on the distance matrix obtained from the animals’ genomic relationships. The Wright’s fixation F_ST index was also employed to assess herds’ genetic homogeneity and to indirectly estimate past migration rates. The Wright’s fixation F_IS index and genomic inbreeding coefficients based on the genomic relationship matrix as well as on runs of homozygosity were also estimated. The Frizarta breed displays relatively low LD levels with r² and |Dʹ| equal to 0.18 and 0.50, respectively, at an average inter-marker distance of 31 kb. Linkage disequilibrium decayed rapidly by distance and persisted over just a few thousand base pairs. Rate of LD decay (β) varied widely among the 26 autosomes with larger values estimated for shorter chromosomes (e.g. β=0.057, for OAR6) and smaller values for longer ones (e.g. β=0.022, for OAR2). The inferred effective population size at the beginning of the breed’s formation was as high as 549, was then reduced to 463 in 1981 (end of the breed’s formation) and further declined to 187, one generation ago. Multidimensional scaling analysis and k-means clustering suggested a genetically homogenous population, F_ST estimates indicated relatively low genetic differentiation between herds, whereas a heat map of the animals’ genomic kinship relationships revealed a stratified population, at a herd level. Estimates of genomic inbreeding coefficients suggested that most recent parental relatedness may have been a major determinant of the current effective population size. A denser than the 50k SNP panel may be more beneficial when performing genome wide association studies in the breed.     

Recombination diversifies chloroplast trnF pseudogenes in Arabidopsis lyrata

Extensive intraspecific variation in the chloroplast trnL(UAA)-trnF(GAA) spacer of model plant Arabidopsis lyrata is caused by multiple copies of a tandemly repeated trnF pseudogene undergoing parallel independent changes in copy number. Linkage disequilibrium and secondary structure analyses indicate that the diversification of pseudogene copies is driven by complex processes of structurally mediated illegitimate recombination. Disperse repeats sharing similar secondary structures interact, facilitating reciprocal exchange of structural motifs between copies via intramolecular and intermolecular recombinations, forming chimeric sequences and iterative expansion and contraction in pseudogene copy numbers. Widely held assumptions that chloroplast sequence evolution is simple and structural changes are informative are violated. Our findings have important implications for the use of this highly variable region in Brassicaceae studies. The reticulate evolution and nonindependent nucleotide substitution render the pseudogene inappropriate for standard phylogenetic reconstruction, but over short evolutionary timescales they may be useful for assessing gene flow, hybridization and introgression.     

On the correlation between heterozygosity and fitness in natural populations

Three primary hypotheses currently prevail for correlations between heterozygosity at a set of molecular markers and fitness in natural populations. First, multilocus heterozygosity-fitness correlations might result from selection acting directly on the scored loci, such as at particular allozyme loci. Second, significant levels of linkage disequilibrium, as in recently bottlenecked-and-expanded populations, might cause associations between the markers and fitness loci in the local chromosomal vicinity. Third, in partially inbred populations, heterozygosity at the markers might reflect variation in the inbreeding coefficient and might associate with fitness as a result of effects of homozygosity at genome-wide distributed loci. Despite years of research, the relative importance of these hypotheses remains unclear. The screening of heterozygosity at polymorphic DNA markers offers an opportunity to resolve this issue, and relevant empirical studies have now emerged. We provide an account of the recent progress on the subject, and give suggestions on how to distinguish between the three hypotheses in future studies.     

拟南芥和琴叶拟南芥中MADS-box基因的比较进化分析

MADS-box基因编码一类转录因子。在被子植物中,MADS-box基因对于营养生长和生殖发育都有重要的调控作用,是植物体(特别是花序、花和果实)的正常发育所不可或缺的。为了理解近缘物种在遗传基础上的异同,我们对拟南芥(Arabidopsis thaliana)和琴叶拟南芥(A.lyrata)基因组中MADS-box基因的拷贝数目和进化式样进行了比较分析。通过搜索公共数据库,我们在拟南芥和琴叶拟南芥中分别鉴定出了106和115个基因。系统发育分析的结果表明,这些基因属于I型和II型MADS-box基因。在两个物种分化之后,II型基因的拷贝数目变化不大,I型基因则经历了多次独立的基因丢失和获得事件。通过比较这些基因在染色体上的排列,我们不但鉴定出了存在微共线性的基因组区段,而且发现新基因产生的主要机制是串联重复和散在重复。分子进化的研究进一步表明,I型和II型基因在进化式样上存在着显著差异:II型基因在进化中一般都受到了较强的选择压力,而I型基因大多受到的选择压力较弱。本研究将为深入理解近缘物种在基因和基因组层面上的异同、探讨物种分化和生物多样性形成的机制等问题提供新思路。