Search

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., Author Summary Understanding how genetic variation can control phenotypic variation is a fundamental goal of modern biology. We combined genome-wide association mapping with metabolomics in the plant Arabidopsis thaliana to explore how species-wide genetic variation controls metabolism. We identified numerous naturally-variable genes that may influence plant metabolism, often clustering in “hotspots.” These hotspots were proximal to selective sweeps, regions of the genome showing decreased diversity possibly from a strong selective advantage of specific variants within the region. This suggests that metabolism may be connected to the selective advantage. Interestingly, metabolite variation in wild Arabidopsis is highly constrained despite the significant genetic variation, thus providing the plant un-sampled metabolic space if the environment shifts. The observed structuring of genetic and metabolic variation suggests individual convergence upon similar phenotypes via different genotypes, possibly intra-specific parallel evolution. This phenotypic convergence couples with a pattern of genotype—phenotype association consistent with metabolite variation largely controlled by numerous small effect genetic variants. This supports the supposition that large magnitude variation is likely unstable in a complex and interconnected metabolism. If this pattern proves generally applicable to other species, it could present a significant hurdle to identifying genes controlling metabolic trait variation via genome-wide association studies.

The functional structure of all biologically active molecules is dependent on intra- and inter-molecular interactions. This is especially evident for RNA molecules whose functionality, maturation, and regulation require formation of correct secondary structure through encoded base-pairing interactions. Unfortunately, intra- and inter-molecular base-pairing information is lacking for most RNAs. Here, we marry classical nuclease-based structure mapping techniques with high-throughput sequencing technology to interrogate all base-paired RNA in Arabidopsis thaliana and identify ∼200 new small (sm)RNA–producing substrates of RNA–DEPENDENT RNA POLYMERASE6. Our comprehensive analysis of paired RNAs reveals conserved functionality within introns and both 5′ and 3′ untranslated regions (UTRs) of mRNAs, as well as a novel population of functional RNAs, many of which are the precursors of smRNAs. Finally, we identify intra-molecular base-pairing interactions to produce a genome-wide collection of RNA secondary structure models. Although our methodology reveals the pairing status of RNA molecules in the absence of cellular proteins, previous studies have demonstrated that structural information obtained for RNAs in solution accurately reflects their structure in ribonucleoprotein complexes. Furthermore, our identification of RNA–DEPENDENT RNA POLYMERASE6 substrates and conserved functional RNA domains within introns and both 5′ and 3′ untranslated regions (UTRs) of mRNAs using this approach strongly suggests that RNA molecules are correctly folded into their secondary structure in solution. Overall, our findings highlight the importance of base-paired RNAs in eukaryotes and present an approach that should be widely applicable for the analysis of this key structural feature of RNA., Author Summary At the heart of RNA functionality, maturation, and regulation is the formation of intricate secondary structures that are dependent on specific nucleotide base-pairing interactions encoded within their sequences. These interactions can either be within (intra-molecular) or between (inter-molecular (heteroduplex)) RNA molecules. Although it is clear that secondary structure is abundantly important for the functionality and regulation of RNAs, comprehensive base-pairing interaction data are completely lacking for the majority of these molecules. To address this, we have developed a new approach for studying the base-pairing interactions of RNA molecules by marrying classical nuclease-based structure mapping techniques with high-throughput sequencing technology. We have used this approach to identify known and novel substrates of the base-paired RNA producing enzyme RNA–DEPENDENT RNA POLYMERASE6, reveal conserved functionality within introns and both 5′ and 3′ untranslated regions (UTRs) of mRNAs, uncover a novel population of functional RNAs, and produce a genome-wide collection of RNA secondary structure models by identifying the base-pairing interactions within each RNA molecule. Our findings demonstrate that our methodology should be widely applicable for the identification and analysis of base-paired RNAs in all biological organisms.

A set of experiments demonstrates the involvement of a chromosomal inversion in the adaptive transition between annual and perennial ecotypes of the yellow monkeyflower, Mimulus guttatus, The role of chromosomal inversions in adaptation and speciation is controversial. Historically, inversions were thought to contribute to these processes either by directly causing hybrid sterility or by facilitating the maintenance of co-adapted gene complexes. Because inversions suppress recombination when heterozygous, a recently proposed local adaptation mechanism predicts that they will spread if they capture alleles at multiple loci involved in divergent adaptation to contrasting environments. Many empirical studies have found inversion polymorphisms linked to putatively adaptive phenotypes or distributed along environmental clines. However, direct involvement of an inversion in local adaptation and consequent ecological reproductive isolation has not to our knowledge been demonstrated in nature. In this study, we discovered that a chromosomal inversion polymorphism is geographically widespread, and we test the extent to which it contributes to adaptation and reproductive isolation under natural field conditions. Replicated crosses between the prezygotically reproductively isolated annual and perennial ecotypes of the yellow monkeyflower, Mimulus guttatus, revealed that alternative chromosomal inversion arrangements are associated with life-history divergence over thousands of kilometers across North America. The inversion polymorphism affected adaptive flowering time divergence and other morphological traits in all replicated crosses between four pairs of annual and perennial populations. To determine if the inversion contributes to adaptation and reproductive isolation in natural populations, we conducted a novel reciprocal transplant experiment involving outbred lines, where alternative arrangements of the inversion were reciprocally introgressed into the genetic backgrounds of each ecotype. Our results demonstrate for the first time in nature the contribution of an inversion to adaptation, an annual/perennial life-history shift, and multiple reproductive isolating barriers. These results are consistent with the local adaptation mechanism being responsible for the distribution of the two inversion arrangements across the geographic range of M. guttatus and that locally adaptive inversion effects contribute directly to reproductive isolation. Such a mechanism may be partially responsible for the observation that closely related species often differ by multiple chromosomal rearrangements., Author Summary Genome rearrangements that change the order of genes along a chromosome are known as inversions and have long been hypothesized to be involved in the origin of species. Yet the way inversions contribute to adaptation and speciation remains mysterious. In this study, we identified a geographically widespread adaptive inversion polymorphism in the yellow monkeyflower, Mimulus guttatus. One arrangement of the inverted region is found in an annual ecotype of this species that lives in Mediterranean habitats characterized by reduced soil water availability in the summer. The other arrangement appears in a perennial ecotype that lives in habitats with high year-round soil moisture. The inversion was observed to influence morphological and flowering time differences between these ecotypes across most of western North America. To test whether the inversion polymorphism contributes to adaptation and reproductive isolation, we conducted a field experiment by breeding plants to reciprocally swap the alternative chromosomal arrangements between the annual and perennial genetic backgrounds. We demonstrated that this inversion polymorphism contributes to local adaptation, the annual/perennial life-history transition, and three reproductive isolating barriers. These results are consistent with the theory that adaptation to local environments can drive the spread of chromosomal inversions and promote speciation.

Background - Posterior mapping is an increasingly popular hierarchical Bayesian based method used to infer character histories and reconstruct ancestral states at nodes of molecular phylogenies, notably of morphological characters. As for all Bayesian analyses specification of prior values is an integrative and important part of the analysis. He we provide an example of how alternative prior choices can seriously influence results and mislead interpretations. Methods/Principal Findings - For two contrasting discrete morphological characters, namely a slow and a fast evolving character found in the plant family Annonaceae, we specified a total of eight different prior distributions per character. We investigated how these prior settings affected important summary statistics. Our analyses showed that the different prior distributions had marked effects on the results in terms of average number of character state changes. These differences arise because priors play a crucial role in determining which areas of parameter space the values of the simulation will be drawn from, independent of the data at hand. However, priors seemed to fit the data better if they would result in a more even sampling of parameter space (normal posterior distribution), in which case alternative standard deviation values had little effect on the results. The most probable character history for each character was affected differently by the prior. For the slower evolving character, the same character history always had the highest posterior probability independent of the priors used. In contrast, the faster evolving character showed different most probable character histories depending on the prior. These differences could be related to the level of homoplasy exhibited by each character. Conclusions - Although our analyses were restricted to two morphological characters within a single family, our results underline the importance of carefully choosing prior values for posterior mapping. Prior specification will be of crucial importance when interpreting the results in a meaningful way. It is hard to suggest a statistically sound method for prior specification without more detailed studies. Meanwhile, we propose that the data could be used to estimate the prior value of the gamma distribution placed on the transformation rate in posterior mapping

Background The endosymbiotic birth of organelles is accompanied by massive transfer of endosymbiont genes to the eukaryotic host nucleus. In the centric diatom Thalassiosira pseudonana the Psb28 protein is encoded in the plastid genome while a second version is nuclear-encoded and possesses a bipartite N-terminal presequence necessary to target the protein into the diatom complex plastid. Thus it can represent a gene captured during endosymbiotic gene transfer. Methodology/Principal Findings To specify the origin of nuclear- and plastid-encoded Psb28 in T. pseudonana we have performed extensive phylogenetic analyses of both mentioned genes. We have also experimentally tested the intracellular location of the nuclear-encoded Psb28 protein (nuPsb28) through transformation of the diatom Phaeodactylum tricornutum with the gene in question fused to EYFP. Conclusions/Significance We show here that both versions of the psb28 gene in T. pseudonana are transcribed. We also provide experimental evidence for successful targeting of the nuPsb28 fused with EYFP to the diatom complex plastid. Extensive phylogenetic analyses demonstrate that nucleotide composition of the analyzed genes deeply influences the tree topology and that appropriate methods designed to deal with a compositional bias of the sequences and the long branch attraction artefact (LBA) need to be used to overcome this obstacle. We propose that nuclear psb28 in T. pseudonana is a duplicate of a plastid localized version, and that it has been transferred from its endosymbiont.

Much of the eukaryotic genome is known to be mobile, largely due to the movement of transposons and other parasitic elements. Recent work in plants and Drosophila suggests that mobility is also a feature of many nontransposon genes and gene families. Indeed, analysis of the Arabidopsis genome suggested that as many as half of all genes had moved to unlinked positions since Arabidopsis diverged from papaya roughly 72 million years ago, and that these mobile genes tend to fall into distinct gene families. However, the mechanism by which single gene transposition occurred was not deduced. By comparing two closely related species, Arabidopsis thaliana and Arabidopsis lyrata, we sought to determine the nature of gene transposition in Arabidopsis. We found that certain categories of genes are much more likely to have transposed than others, and that many of these transposed genes are flanked by direct repeat sequence that was homologous to sequence within the orthologous target site in A. lyrata and which was predominantly genic in identity. We suggest that intrachromosomal recombination between tandemly duplicated sequences, and subsequent insertion of the circular product, is the predominant mechanism of gene transposition., Author Summary Repetitive DNA, such as satellite repeats and transposons, is ubiquitous throughout the genome. Such repeats have been associated with DNA loss, circle formation, and gene transposition in plants and Drosophila. In this work we suggest that, in plants, one mechanism of gene mobility is intrachromosomal recombination via tandem repeats. In addition, we have demonstrated that the classes of genes that tend to form tandem duplications are more likely to have transposed than other gene classes. We conclude that tandem duplications may particularly facilitate gene excision and may also provide targets for gene insertion.

Background The concept of DNA barcoding for species identification has gained considerable momentum in animals because of fairly successful species identification using cytochrome oxidase I (COI). In plants, matK and rbcL have been proposed as standard barcodes. However, barcoding in complex genera is a challenging task. Methodology and Principal Findings We investigated the species discriminatory power of four reportedly most promising plant DNA barcoding loci (one from nuclear genome- ITS, and three from plastid genome- trnH-psbA, rbcL and matK) in species of Indian Berberis L. (Berberidaceae) and two other genera, Ficus L. (Moraceae) and Gossypium L. (Malvaceae). Berberis species were delineated using morphological characters. These characters resulted in a well resolved species tree. Applying both nucleotide distance and nucleotide character-based approaches, we found that none of the loci, either singly or in combinations, could discriminate the species of Berberis. ITS resolved all the tested species of Ficus and Gossypium and trnH-psbA resolved 82% of the tested species in Ficus. The highly regarded matK and rbcL could not resolve all the species. Finally, we employed amplified fragment length polymorphism test in species of Berberis to determine their relationships. Using ten primer pair combinations in AFLP, the data demonstrated incomplete species resolution. Further, AFLP analysis showed that there was a tendency of the Berberis accessions to cluster according to their geographic origin rather than species affiliation. Conclusions/Significance We reconfirm the earlier reports that the concept of universal barcode in plants may not work in a number of genera. Our results also suggest that the matK and rbcL, recommended as universal barcode loci for plants, may not work in all the genera of land plants. Morphological, geographical and molecular data analyses of Indian species of Berberis suggest probable reticulate evolution and thus barcode markers may not work in this case.

An important objective in genome research is to relate genome structure to gene function. Sequence comparisons among orthologous and paralogous genes and their allelic variants can reveal sequences of functional significance. Here, we describe a 379-kb region on chromosome 1 of maize that enables us to reconstruct chromosome breakage, transposition, non-homologous end-joining, and homologous recombination events. Such a high-density composition of various mechanisms in a small chromosomal interval exemplifies the evolution of gene regulation and allelic diversity in general. It also illustrates the evolutionary pace of changes in plants, where many of the above mechanisms are of somatic origin. In contrast to animals, somatic alterations can easily be transmitted through meiosis because the germline in plants is contiguous to somatic tissue, permitting the recovery of such chromosomal rearrangements. The analyzed region contains the P1-wr allele, a variant of the genetically well-defined p1 gene, which encodes a Myb-like transcriptional activator in maize. The P1-wr allele consists of eleven nearly perfect P1-wr 12-kb repeats that are arranged in a tandem head-to-tail array. Although a technical challenge to sequence such a structure by shotgun sequencing, we overcame this problem by subcloning each repeat and ordering them based on nucleotide variations. These polymorphisms were also critical for recombination and expression analysis in presence and absence of the trans-acting epigenetic factor Ufo1. Interestingly, chimeras of the p1 and p2 genes, p2/p1 and p1/p2, are framing the P1-wr cluster. Reconstruction of sequence amplification steps at the p locus showed the evolution from a single Myb-homolog to the multi-gene P1-wr cluster. It also demonstrates how non-homologous end-joining can create novel gene fusions. Comparisons to orthologous regions in sorghum and rice also indicate a greater instability of the maize genome, probably due to diploidization following allotetraploidization., Author Summary Plant genomes analyzed to date contain 15% or more genes that are arranged in tandem arrays. Tandem duplications are a source for allelic variability since their homologous sequences can serve in recombination events. For example, unequal crossing over between amplified genes can result in contraction and expansion of the array. Tandem gene multiplications are also subject to repeat induced gene silencing (RIGS). Most importantly, gene duplications create the evolutionary potential for genetic novelty (neo- or subfunctionalization). In addition to homologous recombination during meiosis, illegitimate recombination in somatic tissues of plants can create events that potentially can be transmitted through reproductive tissue to further enrich genetic diversity. Here we illustrate the evolution from a single Myb homolog to a multigene cluster that exemplifies the evolution of the maize genome. We used the p locus to demonstrate how plant genomes expand by polyploidization, gene duplication, and transposition. We characterized in detail the structural changes at the p cluster that resulted from genomic instability. Because structure determines function, we linked genomic rearrangements at the P1-wr cluster to functional consequences. At the P1-wr locus, structural changes caused regulatory/transcriptional modifications that in turn give rise to phenotypic alterations.

Background Genome level analyses have enhanced our view of phylogenetics in many areas of the tree of life. With the production of whole genome DNA sequences of hundreds of organisms and large-scale EST databases a large number of candidate genes for inclusion into phylogenetic analysis have become available. In this work, we exploit the burgeoning genomic data being generated for plant genomes to address one of the more important plant phylogenetic questions concerning the hierarchical relationships of the several major seed plant lineages (angiosperms, Cycadales, Gingkoales, Gnetales, and Coniferales), which continues to be a work in progress, despite numerous studies using single, few or several genes and morphology datasets. Although most recent studies support the notion that gymnosperms and angiosperms are monophyletic and sister groups, they differ on the topological arrangements within each major group. Methodology We exploited the EST database to construct a supermatrix of DNA sequences (over 1,200 concatenated orthologous gene partitions for 17 taxa) to examine non-flowering seed plant relationships. This analysis employed programs that offer rapid and robust orthology determination of novel, short sequences from plant ESTs based on reference seed plant genomes. Our phylogenetic analysis retrieved an unbiased (with respect to gene choice), well-resolved and highly supported phylogenetic hypothesis that was robust to various outgroup combinations. Conclusions We evaluated character support and the relative contribution of numerous variables (e.g. gene number, missing data, partitioning schemes, taxon sampling and outgroup choice) on tree topology, stability and support metrics. Our results indicate that while missing characters and order of addition of genes to an analysis do not influence branch support, inadequate taxon sampling and limited choice of outgroup(s) can lead to spurious inference of phylogeny when dealing with phylogenomic scale data sets. As expected, support and resolution increases significantly as more informative characters are added, until reaching a threshold, beyond which support metrics stabilize, and the effect of adding conflicting characters is minimized.

Background: The B3 DNA binding domain includes five families: auxin response factor (ARF), abscisic acid-insensitive3 (ABI3), high level expression of sugar inducible (HSI), related to ABI3/VP1 (RAV) and reproductive meristem (REM). The release of the complete genomes of the angiosperm eudicots Arabidopsis thaliana and Populus trichocarpa, the monocot Orysa sativa, the bryophyte Physcomitrella patens,the green algae Chlamydomonas reinhardtii and Volvox carteri and the red algae Cyanidioschyzon melorae provided an exceptional opportunity to study the evolution of this superfamily. Methodology: In order to better understand the origin and the diversification of B3 domains in plants, we combined comparative phylogenetic analysis with exon/intron structure and duplication events. In addition, we investigated the conservation and divergence of the B3 domain during the origin and evolution of each family. Conclusions: Our data indicate that showed that the B3 containing genes have undergone extensive duplication events, and that the REM family B3 domain has a highly diverged DNA binding. Our results also indicate that the founding member of the B3 gene family is likely to be similar to the ABI3/HSI genes found in C. reinhardtii and V. carteri. Among the B3 families, ABI3, HSI, RAV and ARF are most structurally conserved, whereas the REM family has experienced a rapid divergence. These results are discussed in light of their functional and evolutionary roles in plant development. © 2009 Romanel et al.

A common yet poorly understood evolutionary transition among flowering plants is a switch from outbreeding to an inbreeding mode of mating. The model plant Arabidopsis thaliana evolved to an inbreeding state through the loss of self-incompatibility, a pollen-rejection system in which pollen recognition by the stigma is determined by tightly linked and co-evolving alleles of the S-locus receptor kinase (SRK) and its S-locus cysteine-rich ligand (SCR). Transformation of A. thaliana, with a functional AlSRKb-SCRb gene pair from its outcrossing relative A. lyrata, demonstrated that A. thaliana accessions harbor different sets of cryptic self-fertility–promoting mutations, not only in S-locus genes, but also in other loci required for self-incompatibility. However, it is still not known how many times and in what manner the switch to self-fertility occurred in the A. thaliana lineage. Here, we report on our identification of four accessions that are reverted to full self-incompatibility by transformation with AlSRKb-SCRb, bringing to five the number of accessions in which self-fertility is due to, and was likely caused by, S-locus inactivation. Analysis of S-haplotype organization reveals that inter-haplotypic recombination events, rearrangements, and deletions have restructured the S locus and its genes in these accessions. We also perform a Quantitative Trait Loci (QTL) analysis to identify modifier loci associated with self-fertility in the Col-0 reference accession, which cannot be reverted to full self-incompatibility. Our results indicate that the transition to inbreeding occurred by at least two, and possibly more, independent S-locus mutations, and identify a novel unstable modifier locus that contributes to self-fertility in Col-0., Author Summary The mating system adopted by a species has a profound influence on extent of polymorphism, population structure, and evolutionary potential. In flowering plants, the switch from outbreeding to inbreeding has occurred repeatedly, yet little is known about the underlying genetic events. This is true even for the model species A. thaliana, a highly self-fertile member of the crucifer family. In this family, outbreeding is enforced by a self-incompatibility system controlled by the S locus, which involves the recognition of pollen by the stigma to prevent self-fertilization and familial inbreeding. We recently demonstrated that A. thaliana accessions may be reverted to full or partial self-incompatibility by transformation with S-locus genes isolated from its close self-incompatible relative A. lyrata. Despite much recent debate, however, we still do not know how A. thaliana became self-fertile. Here, we use our recently established A. thaliana transgenic self-incompatible experimental model to address these issues. Analysis of the S locus in accessions that can be reverted to full self-incompatibility demonstrates that self-fertility in A. thaliana arose by at least two independent S-locus mutations. Furthermore, analysis of an accession that expresses only partial self-incompatibility shows that self-fertility is associated with an unstable allele at a locus unlinked to the S locus.

Eukaryotes bearing red alga-derived plastids — photosynthetic alveolates (dinoflagellates plus the apicomplexan Toxoplasma gondii plus the chromerid Chromera velia), photosynthetic stramenopiles, haptophytes, and cryptophytes — possess unique plastid-targeted glyceraldehyde-3-phosphate dehydrogenases (henceforth designated as “GapC1”). Pioneering phylogenetic studies have indicated a single origin of the GapC1 enzymes in eukaryotic evolution, but there are two potential idiosyncrasies in the GapC1 phylogeny: Firstly, the GapC1 tree topology is apparently inconsistent with the organismal relationship among the “GapC1-containing” groups. Secondly, four stramenopile GapC1 homologues are consistently paraphyletic in previously published studies, although these organisms have been widely accepted as monophyletic. For a closer examination of the above issues, in this study GapC1 gene sampling was improved by determining/identifying nine stramenopile and two cryptophyte genes. Phylogenetic analyses of our GapC1 dataset, which is particularly rich in the stramenopile homologues, prompt us to propose a new scenario that assumes multiple, lateral GapC1 gene transfer events to explain the incongruity between the GapC1 phylogeny and the organismal relationships amongst the “GapC1-containing” groups. Under our new scenario, GapC1 genes uniquely found in photosynthetic alveolates, photosynthetic stramenopiles, haptophytes, and cryptopyhytes are not necessarily a character vertically inherited from a common ancestor.

Background\ud \ud Archaeobotanical samples of barley (Hordeum vulgare L.) found at Qasr Ibrim display a two-row phenotype that is unique to the region of archaeological sites upriver of the first cataract of the Nile, characterised by the development of distinctive lateral bracts. The phenotype occurs throughout all strata at Qasr Ibrim, which range in age from 3000 to a few hundred years.\ud Methodology and Findings\ud \ud We extracted ancient DNA from barley samples from the entire range of occupancy of the site, and studied the Vrs1 gene responsible for row number in extant barley. Surprisingly, we found a discord between the genotype and phenotype in all samples; all the barley had a genotype consistent with the six-row condition. These results indicate a six-row ancestry for the Qasr Ibrim barley, followed by a reassertion of the two-row condition. Modelling demonstrates that this sequence of evolutionary events requires a strong selection pressure.\ud Conclusions\ud \ud The two-row phenotype at Qasr Ibrim is caused by a different mechanism to that in extant barley. The strength of selection required for this mechanism to prevail indicates that the barley became locally adapted in the region in response to a local selection pressure. The consistency of the genotype/phenotype discord over time supports a scenario of adoption of this barley type by successive cultures, rather than the importation of new barley varieties associated with individual cultures.

BACKGROUND: Sexual conflict theory predicts sexually antagonistic coevolution of reproductive traits driven by conflicting evolutionary interests of two reproducing individuals. Most studies of the evolutionary consequences of sexual conflicts have, however, to date collectively investigated only a few species. In this study we used the annual herb Collinsia heterophylla to experimentally test the existence and evolutionary consequences of a potential sexual conflict over onset of stigma receptivity. METHODOLOGY/PRINCIPAL FINDINGS: We conducted crosses within and between four greenhouse-grown populations originating from two regions. Our experimental setup allowed us to investigate male-female interactions at three levels of geographic distances between interacting individuals. Both recipient and pollen donor identity affected onset of stigma receptivity within populations, confirming previous results that some pollen donors can induce stigma receptivity. We also found that donors were generally better at inducing stigma receptivity following pollen deposition on stigmas of recipients from another population than their own, especially within a region. On the other hand, we found that donors did worse at inducing stigma receptivity in crosses between regions. Interestingly, recipient costs in terms of lowered seed number after early fertilisation followed the same pattern: the cost was apparent only if the pollen donor belonged to the same region as the recipient. CONCLUSION/SIGNIFICANCE: Our results indicate that recipients are released from the cost of interacting with local pollen donors when crossed with donors from a more distant location, a pattern consistent with a history of sexually antagonistic coevolution within populations. Accordingly, sexual conflicts may have important evolutionary consequences also in plants. (Less)

Most phenotypic variation present in natural populations is under polygenic control, largely determined by genetic variation at quantitative trait loci (QTLs). These genetic loci frequently interact with the environment, development, and each other, yet the importance of these interactions on the underlying genetic architecture of quantitative traits is not well characterized. To better study how epistasis and development may influence quantitative traits, we studied genetic variation in Arabidopsis glucosinolate activation using the moderately sized Bayreuth×Shahdara recombinant inbred population, in terms of number of lines. We identified QTLs for glucosinolate activation at three different developmental stages. Numerous QTLs showed developmental dependency, as well as a large epistatic network, centered on the previously cloned large-effect glucosinolate activation QTL, ESP. Analysis of Heterogeneous Inbred Families validated seven loci and all of the QTL×DPG (days post-germination) interactions tested, but was complicated by the extensive epistasis. A comparison of transcript accumulation data within 211 of these RILs showed an extensive overlap of gene expression QTLs for structural specifiers and their homologs with the identified glucosinolate activation loci. Finally, we were able to show that two of the QTLs are the result of whole-genome duplications of a glucosinolate activation gene cluster. These data reveal complex age-dependent regulation of structural outcomes and suggest that transcriptional regulation is associated with a significant portion of the underlying ontogenic variation and epistatic interactions in glucosinolate activation., Author Summary A principal interest in biology is to understand how natural genetic variation translates into phenotypic variation. A key component of this connection is how the genetic variation interacts with other sources of variation, such as environment (G×E), development (G×D), or other genetic loci (G×G or epistasis). To analyze the molecular underpinnings of these quantitative genetics interaction terms, we investigated the genetic architecture of an adaptive trait, glucosinolate activation, in Arabidopsis thaliana during the development of what is considered a static mature rosette. Variation in glucosinolate activation was principally controlled by epistatic and G×D interactions. Epistatic interactions identified both Mendelian epistasis, where regulatory loci controlled enzymatic loci, and quantitative interactions between regulatory loci. G×D appeared to involve master regulatory loci as determined by trans-eQTL hotspot analysis. Finally, two common glucosinolate activation QTLs appear to have evolved via gene loss and sub-functionalization following quadruplication of an ancestral genomic fragment, potentially by two whole-genome duplications. Thus, genomic duplication events may facilitate the formation of quantitative genetic variation. This study provides insights into the molecular basis of the link between genetic and phenotypic variation in a potentially adaptive trait.

BACKGROUND: The genus Oryza is being used as a model in plant genomic studies although there are several issues still to be resolved regarding the spatio-temporal evolution of this ancient genus. Particularly contentious is whether undated transoceanic natural dispersal or recent human interference has been the principal agent determining its present distribution and differentiation. In this context, we studied the origin and distribution history of the allotetraploid CD rice genome. It is endemic to the Neotropics but the genus is thought to have originated in the Paleotropics, and there is relatively little genetic divergence between some orthologous sequences of the C genome component and their Old World counterparts. METHODOLOGY/PRINCIPAL FINDINGS: Because of its allotetraploidy, there are several potential pitfalls in trying to date the formation of the CD genome using molecular data and this could lead to erroneous estimates. Therefore, we rather chose to rely on historical evidence to determine whether or not the CD genome was present in the Neotropics before the arrival of Columbus. We searched early collections of herbarium specimens and studied the reports of explorers of the tropical Americas for references to rice. In spite of numerous collectors traveling inland and collecting Oryza, plants determined as CD genome species were not observed away from cultivated rice fields until 1869. Various arguments suggest that they only consisted of weedy forms until that time. CONCLUSIONS/SIGNIFICANCE: The spatio-temporal distribution of herbarium collections fits a simple biogeographical scenario for the emergence in cultivated rice fields followed by radiation in the wild of the CD genome in the Neotropics during the last four centuries. This probably occurred from species introduced to the Americas by humans and we found no evidence that the CD genome pre-existed in the Old World. We therefore propose a new evolutionary hypothesis for such a recent origin of the CD genome. Moreover, we exemplify how an historical approach can provide potentially important information and help to disentangle the timing of evolutionary events in the history of the Oryza genomes.

BackgroundMany of the processes affecting genetic diversity act on local populations. However, studies of plant nucleotide diversity have largely ignored local sampling, making it difficult to infer the demographic history of populations and to assess the importance of local adaptation. Arabidopsis lyrata, a self-incompatible, perennial species with a circumpolar distribution, is an excellent model system in which to study the roles of demographic history and local adaptation in patterning genetic variation.Principal findingsWe studied nucleotide diversity in six natural populations of Arabidopsis lyrata, using 77 loci sampled from 140 chromosomes. The six populations were highly differentiated, with a median FST of 0.52, and structure analysis revealed no evidence of admixed individuals. Average within-population diversity varied among populations, with the highest diversity found in a German population; this population harbors 3-fold higher levels of silent diversity than worldwide samples of A. thaliana. All A. lyrata populations also yielded positive values of Tajima's D. We estimated a demographic model for these populations, finding evidence of population divergence over the past 19,000 to 47,000 years involving non-equilibrium demographic events that reduced the effective size of most populations. Finally, we used the inferred demographic model to perform an initial test for local adaptation and identified several genes, including the flowering time gene FCA and a disease resistance locus, as candidates for local adaptation events.ConclusionsOur results underscore the importance of population-specific, non-equilibrium demographic processes in patterning diversity within A. lyrata. Moreover, our extensive dataset provides an important resource for future molecular population genetic studies of local adaptation in A. lyrata.

The latitudinal biodiversity gradient remains one of the most widely recognized yet puzzling patterns in nature [1]. Presently, the high level of extinction of tropical species, referred to as the “tropical biodiversity crisis”, has the potential to erode this pattern. While the connection between species richness, extinction, and speciation has long intrigued biologists [2], [3], these interactions have experienced increased poignancy due to their relevancy to where we should concentrate our conservation efforts. Natural extinction is a phenomenon thought to have its own latitudinal gradient, with lower extinction rates in the tropics being reported in beetles, birds, mammals, and bivalves [4]–[7]. Processes that have buffered ecosystems from high extinction rates in the past may also buffer ecosystems against disturbance of anthropogenic origin. While potential parallels between historical and present-day extinction patterns have been acknowledged, they remain only superficially explored and plant extinction patterns have been particularly neglected. Studies on the disappearances of animal species have reached conflicting conclusions, with the rate of extinction appearing either higher [8] or lower [9] in species richness hotspots. Our global study of extinction risk in vascular plants finds disproportionately higher extinction risk in tropical countries, even when indicators of human pressure (GDP, population density, forest cover change) are taken into account. Our results are at odds with the notion that the tropics represent a museum of plant biodiversity (places of historically lowered extinction) and we discuss mechanisms that may reconcile this apparent contradiction.

BackgroundIn eukaryotes the photosynthetic antenna system is composed of subunits encoded by the light harvesting complex (Lhc) multigene family. These proteins play a key role in photosynthesis and are involved in both light harvesting and photoprotection. The moss Physcomitrella patens is a member of a lineage that diverged from seed plants early after land colonization and therefore by studying this organism, we may gain insight into adaptations to the aerial environment.Principal findingsIn this study, we characterized the antenna protein multigene family in Physcomitrella patens, by sequence analysis as well as biochemical and functional investigations. Sequence identification and analysis showed that some antenna polypeptides, such as Lhcb3 and Lhcb6, are present only in land organisms, suggesting they play a role in adaptation to the sub-aerial environment. Our functional analysis which showed that photo-protective mechanisms in Physcomitrella patens are very similar to those in seed plants fits with this hypothesis. In particular, Physcomitrella patens also activates Non Photochemical Quenching upon illumination, consistent with the detection of an ortholog of the PsbS protein. As a further adaptation to terrestrial conditions, the content of Photosystem I low energy absorbing chlorophylls also increased, as demonstrated by differences in Lhca3 and Lhca4 polypeptide sequences, in vitro reconstitution experiments and low temperature fluorescence spectra.ConclusionsThis study highlights the role of Lhc family members in environmental adaptation and allowed proteins associated with mechanisms of stress resistance to be identified within this large family.

BackgroundTranscription factors of the basic leucine zipper (bZIP) family control important processes in all eukaryotes. In plants, bZIPs are regulators of many central developmental and physiological processes including photomorphogenesis, leaf and seed formation, energy homeostasis, and abiotic and biotic stress responses. Here we performed a comprehensive phylogenetic analysis of bZIP genes from algae, mosses, ferns, gymnosperms and angiosperms.Methodology/principal findingsWe identified 13 groups of bZIP homologues in angiosperms, three more than known before, that represent 34 Possible Groups of Orthologues (PoGOs). The 34 PoGOs may correspond to the complete set of ancestral angiosperm bZIP genes that participated in the diversification of flowering plants. Homologous genes dedicated to seed-related processes and ABA-mediated stress responses originated in the common ancestor of seed plants, and three groups of homologues emerged in the angiosperm lineage, of which one group plays a role in optimizing the use of energy.Conclusions/significanceOur data suggest that the ancestor of green plants possessed four bZIP genes functionally involved in oxidative stress and unfolded protein responses that are bZIP-mediated processes in all eukaryotes, but also in light-dependent regulations. The four founder genes amplified and diverged significantly, generating traits that benefited the colonization of new environments.

Barcodes are short segments of DNA that can be used to uniquely identify an unknown specimen to species, particularly when diagnostic morphological features are absent. These sequences could offer a new forensic tool in plant and animal conservation-especially for endangered species such as members of the Cycadales. Ideally, barcodes could be used to positively identify illegally obtained material even in cases where diagnostic features have been purposefully removed or to release confiscated organisms into the proper breeding population. In order to be useful, a DNA barcode sequence must not only easily PCR amplify with universal or near-universal reaction conditions and primers, but also contain enough variation to generate unique identifiers at either the species or population levels. Chloroplast regions suggested by the Plant Working Group of the Consortium for the Barcode of Life (CBoL), and two alternatives, the chloroplast psbA-trnH intergenic spacer and the nuclear ribosomal internal transcribed spacer (nrITS), were tested for their utility in generating unique identifiers for members of the Cycadales. Ease of amplification and sequence generation with universal primers and reaction conditions was determined for each of the seven proposed markers. While none of the proposed markers provided unique identifiers for all species tested, nrITS showed the most promise in terms of variability, although sequencing difficulties remain a drawback. We suggest a workflow for DNA barcoding, including database generation and management, which will ultimately be necessary if we are to succeed in establishing a universal DNA barcode for plants.

Population structure can confound the identification of correlations in biological data. Such confounding has been recognized in multiple biological disciplines, resulting in a disparate collection of proposed solutions. We examine several methods that correct for confounding on discrete data with hierarchical population structure and identify two distinct confounding processes, which we call coevolution and conditional influence. We describe these processes in terms of generative models and show that these generative models can be used to correct for the confounding effects. Finally, we apply the models to three applications: identification of escape mutations in HIV-1 in response to specific HLA-mediated immune pressure, prediction of coevolving residues in an HIV-1 peptide, and a search for genotypes that are associated with bacterial resistance traits in Arabidopsis thaliana. We show that coevolution is a better description of confounding in some applications and conditional influence is better in others. That is, we show that no single method is best for addressing all forms of confounding. Analysis tools based on these models are available on the internet as both web based applications and downloadable source code at http://atom.research.microsoft.com/bio/phylod.aspx.

In plants, microRNAs (miRNAs) comprise one of two classes of small RNAs that function primarily as negative regulators at the posttranscriptional level. Several MIRNA genes in the plant kingdom are ancient, with conservation extending between angiosperms and the mosses, whereas many others are more recently evolved. Here, we use deep sequencing and computational methods to identify, profile and analyze non-conserved MIRNA genes in Arabidopsis thaliana. 48 non-conserved MIRNA families, nearly all of which were represented by single genes, were identified. Sequence similarity analyses of miRNA precursor foldback arms revealed evidence for recent evolutionary origin of 16 MIRNA loci through inverted duplication events from protein-coding gene sequences. Interestingly, these recently evolved MIRNA genes have taken distinct paths. Whereas some non-conserved miRNAs interact with and regulate target transcripts from gene families that donated parental sequences, others have drifted to the point of non-interaction with parental gene family transcripts. Some young MIRNA loci clearly originated from one gene family but form miRNAs that target transcripts in another family. We suggest that MIRNA genes are undergoing relatively frequent birth and death, with only a subset being stabilized by integration into regulatory networks.

International audience; Background: Host-pathogen interactions that may lead to a competitive co-evolution of virulence and resistance mechanisms present an attractive system to study molecular evolution because strong, recent (or even current) selective pressure is expected at many genomic loci. However, it is unclear whether these selective forces would act to preserve existing diversity, promote novel diversity, or reduce linked neutral diversity during rapid fixation of advantageous alleles. In plants, the lack of adaptive immunity places a larger burden on genetic diversity to ensure survival of plant populations. This burden is even greater if the generation time of the plant is much longer than the generation time of the pathogen. Methodology/Principal Findings: Here, we present nucleotide polymorphism and substitution data for 41 candidate genes from the long-lived forest tree loblolly pine, selected primarily for their prospective influences on host-pathogen interactions. This dataset is analyzed together with 15 drought-tolerance and 13 wood-quality genes from previous studies. A wide range of neutrality tests were performed and tested against expectations from realistic demographic models. Conclusions/Significance: Collectively, our analyses found that axr (auxin response factor), caf1 (chromatin assembly factor) and gatabp1 (gata binding protein 1) candidate genes carry patterns consistent with directional selection and erd3 (early response to drought 3) displays patterns suggestive of a selective sweep, both of which are consistent with the arm-race model of disease response evolution. Furthermore, we have identified patterns consistent with diversifying selection at erf1-like (ethylene responsive factor 1), ccoaoemt (caffeoyl-CoA-O-methyltransferase), cyp450-like (cytochrome p450-like) and pr4.3 (pathogen response 4.3), expected under the trench-warfare evolution model. Finally, a drought-tolerance candidate related to the plant cell wall, lp5, displayed patterns consistent with balancing selection. In conclusion, both arms-race and trench-warfare models seem compatible with patterns of polymorphism found in different disease-response candidate genes, indicating a mixed strategy of disease tolerance evolution for loblolly pine, a major tree crop in southeastern United States.

Alfred Sturtevant, who invented genetic mapping while still an undergraduate, published the first evidence of a chromosomal inversion in 1921 [1]. He suggested then, and later proved, that they have a dramatic effect on transmission: when heterozygous, inversions suppress recombination. Over the next half century, inspired largely by Dobzhansky and his coworkers, much of empirical population genetics devoted itself to studying the abundant polymorphisms within and fixed differences of inversions between species of Drosophila [2]. Starting in the 1970s, this rich literature largely sank from view with the rise of biochemical and then molecular genetics. But inversions are ascendant again. Comparative genomics is now revealing that chromosomes are far more structurally fluid than even Dobzhansky dared to suppose. Where classic cytogenetics identified only nine inversions that distinguish humans and chimpanzees, comparison of their genomic sequences reveals on the order of 1,500 [3] (Figure 1). Despite the importance of inversions as a major mechanism for reorganizing the genome, we are still struggling to understand how and why they evolve almost a century after Sturtevant's discovery. Figure 1 Chromosome inversions that distinguish humans and chimpanzees inferred from a comparison of their genomic sequences [3]. An inversion occurs when a chromosome breaks at two points and the segment bounded by the breakpoints is reinserted in the reversed orientation. Several molecular mechanisms can mediate this event [4]. Box 1 gives an overview of some basic properties of inversions and the ways that they are detected. Box 1. What are chromosome inversions? Inversions are a diverse class of chromsomal mutation. The majority are small (

Background Elymus (Poaceae) is a large genus of polyploid species in the wheat tribe Triticeae. It is polyphyletic, exhibiting many distinct allopolyploid genome combinations, and its history might be further complicated by introgression and lineage sorting. We focus on a subset of Elymus species with a tetraploid genome complement derived from Pseudoroegneria (genome St) and Hordeum (H). We confirm the species' allopolyploidy, identify possible genome donors, and pinpoint instances of apparent introgression or incomplete lineage sorting. Methodology/Principal Findings We sequenced portions of three unlinked nuclear genes—phosphoenolpyruvate carboxylase, β-amylase, and granule-bound starch synthase I—from 27 individuals, representing 14 Eurasian and North American StStHH Elymus species. Elymus sequences were combined with existing data from monogenomic representatives of the tribe, and gene trees were estimated separately for each data set using maximum likelihood. Trees were examined for evidence of allopolyploidy and additional reticulate patterns. All trees confirm the StStHH genome configuration of the Elymus species. They suggest that the StStHH group originated in North America, and do not support separate North American and European origins. Our results point to North American Pseudoroegneria and Hordeum species as potential genome donors to Elymus. Diploid P. spicata is a prospective St-genome donor, though conflict among trees involving P. spicata and the Eurasian P. strigosa suggests either introgression of GBSSI sequences from P. strigosa into North American Elymus and Pseudoroegneria, or incomplete lineage sorting of ancestral GBSSI polymorphism. Diploid H. californicum and/or allotetraploid H. jubatum are possible H-genome donors; direct involvement of an allotetraploid Hordeum species would simultaneously introduce two distinct H genomes to Elymus, consistent with some of the relationships among H-genome sequences in Hordeum and Elymus. Conclusions/Significance Comparisons among molecular phylogenetic trees confirm allopolyploidy, identify potential genome donors, and highlight cases of apparent introgression or incomplete lineage sorting. The complicated history of this group emphasizes an inherent problem with interpreting conflicts among bifurcating trees—identifying introgression and determining its direction depend on which tree is chosen as a starting point of comparison. In spite of difficulties with interpretation, differences among gene trees allow us to identify reticulate species and develop hypotheses about underlying evolutionary processes.

The population structure of an organism reflects its evolutionary history and influences its evolutionary trajectory. It constrains the combination of genetic diversity and reveals patterns of past gene flow. Understanding it is a prerequisite for detecting genomic regions under selection, predicting the effect of population disturbances, or modeling gene flow. This paper examines the detailed global population structure of Arabidopsis thaliana. Using a set of 5,707 plants collected from around the globe and genotyped at 149 SNPs, we show that while A. thaliana as a species self-fertilizes 97% of the time, there is considerable variation among local groups. This level of outcrossing greatly limits observed heterozygosity but is sufficient to generate considerable local haplotypic diversity. We also find that in its native Eurasian range A. thaliana exhibits continuous isolation by distance at every geographic scale without natural breaks corresponding to classical notions of populations. By contrast, in North America, where it exists as an exotic species, A. thaliana exhibits little or no population structure at a continental scale but local isolation by distance that extends hundreds of km. This suggests a pattern for the development of isolation by distance that can establish itself shortly after an organism fills a new habitat range. It also raises questions about the general applicability of many standard population genetics models. Any model based on discrete clusters of interchangeable individuals will be an uneasy fit to organisms like A. thaliana which exhibit continuous isolation by distance on many scales., Author Summary Much of the modern field of population genetics is premised on particular models of what an organism's population structure is and how it behaves. The classic models generally start with the idea of a single randomly mating population that has reached an evolutionary equilibrium. Many models relax some of these assumptions, allowing for phenomena such as assortative mating, discrete sub-populations with migration, self-fertilization, and sex-ratio distortion. Virtually all models, however, have as their core premise the notion that there exist classes of exchangeable individuals each of which represents an identical, independent sample from that class' distribution. For certain organisms, such as Drosophila melanogaster, these models do an excellent job of describing how populations work. For other organisms, such as humans, these models can be reasonable approximations but require a great deal of care in assembling samples and can begin to break down as sampling becomes locally dense. For the vast majority of organisms the applicability of these models has never been investigated.

We describe a comprehensive and general approach for mapping centromeres and present a detailed characterization of two maize centromeres. Centromeres are difficult to map and analyze because they consist primarily of repetitive DNA sequences, which in maize are the tandem satellite repeat CentC and interspersed centromeric retrotransposons of maize (CRM). Centromeres are defined epigenetically by the centromeric histone H3 variant, CENH3. Using novel markers derived from centromere repeats, we have mapped all ten centromeres onto the physical and genetic maps of maize. We were able to completely traverse centromeres 2 and 5, confirm physical maps by fluorescence in situ hybridization (FISH), and delineate their functional regions by chromatin immunoprecipitation (ChIP) with anti-CENH3 antibody followed by pyrosequencing. These two centromeres differ substantially in size, apparent CENH3 density, and arrangement of centromeric repeats; and they are larger than the rice centromeres characterized to date. Furthermore, centromere 5 consists of two distinct CENH3 domains that are separated by several megabases. Succession of centromere repeat classes is evidenced by the fact that elements belonging to the recently active recombinant subgroups of CRM1 colonize the present day centromeres, while elements of the ancestral subgroups are also found in the flanking regions. Using abundant CRM and non-CRM retrotransposons that inserted in and near these two centromeres to create a historical record of centromere location, we show that maize centromeres are fluid genomic regions whose borders are heavily influenced by the interplay of retrotransposons and epigenetic marks. Furthermore, we propose that CRMs may be involved in removal of centromeric DNA (specifically CentC), invasion of centromeres by non-CRM retrotransposons, and local repositioning of the CENH3., Author Summary Centromeres tend to be the last regions to be assembled in genome projects, as their mapping is hampered by their characteristically high repeat DNA content and lack of genetic recombination. Using unique markers derived from these repeat-rich regions, we were able to generate and annotate physical maps of two maize centromeres. Functional centromeres are defined not so much by their primary DNA sequence as by the presence of CENH3, a special histone that replaces canonical histone H3 in centromeric nucleosomes. Little is known about how deposition of CENH3 is regulated, or about the interplay between centromeric repeats and CENH3. By graphing the density of CENH3 nucleosomes onto the physical map, we delineated the functional centromeres in today's maize genome. We then used the large number of LTR retrotransposon insertions, for which the corn genome is well known, as “archeological evidence” to reconstruct the historic centromere boundaries. This was possible because i) some retrotransposon families of maize (CRM) appear to possess a unique ability to preferentially target centromeres during integration and ii) insertion times of individual retrotransposons can be calculated. Here we show that the centromere boundaries in maize have changed over time and are heavily influenced by centromeric and non-centromeric repeats.

Background The Cocoseae is one of 13 tribes of Arecaceae subfam. Arecoideae, and contains a number of palms with significant economic importance, including the monotypic and pantropical Cocos nucifera L., the coconut, the origins of which have been one of the “abominable mysteries” of palm systematics for decades. Previous studies with predominantly plastid genes weakly supported American ancestry for the coconut but ambiguous sister relationships. In this paper, we use multiple single copy nuclear loci to address the phylogeny of the Cocoseae subtribe Attaleinae, and resolve the closest extant relative of the coconut. Methodology/Principal Findings We present the results of combined analysis of DNA sequences of seven WRKY transcription factor loci across 72 samples of Arecaceae tribe Cocoseae subtribe Attaleinae, representing all genera classified within the subtribe, and three outgroup taxa with maximum parsimony, maximum likelihood, and Bayesian approaches, producing highly congruent and well-resolved trees that robustly identify the genus Syagrus as sister to Cocos and resolve novel and well-supported relationships among the other genera of the Attaleinae. We also address incongruence among the gene trees with gene tree reconciliation analysis, and assign estimated ages to the nodes of our tree. Conclusions/Significance This study represents the as yet most extensive phylogenetic analyses of Cocoseae subtribe Attaleinae. We present a well-resolved and supported phylogeny of the subtribe that robustly indicates a sister relationship between Cocos and Syagrus. This is not only of biogeographic interest, but will also open fruitful avenues of inquiry regarding evolution of functional genes useful for crop improvement. Establishment of two major clades of American Attaleinae occurred in the Oligocene (ca. 37 MYBP) in Eastern Brazil. The divergence of Cocos from Syagrus is estimated at 35 MYBP. The biogeographic and morphological congruence that we see for clades resolved in the Attaleinae suggests that WRKY loci are informative markers for investigating the phylogenetic relationships of the palm family.

Background: Genome evolution in the gymnosperm lineage of seed plants has given rise to many of the most complex and largest plant genomes, however the elements involved are poorly understood. Methodology/Principal Findings: Gymny is a previously undescribed retrotransposon family in Pinus that is related to Athila elements in Arabidopsis. Gymny elements are dispersed throughout the modern Pinus genome and occupy a physical space at least the size of the Arabidopsis thaliana genome. In contrast to previously described retroelements in Pinus, the Gymny family was amplified or introduced after the divergence of pine and spruce (Picea). If retrotransposon expansions are responsible for genome size differences within the Pinaceae, as they are in angiosperms, then they have yet to be identified. In contrast, molecular divergence of Gymny retrotransposons together with other families of retrotransposons can account for the large genome complexity of pines along with protein-coding genic DNA, as revealed by massively parallel DNA sequence analysis of Cot fractionated genomic DNA. Conclusions/Significance: Most of the enormous genome complexity of pines can be explained by divergence of retrotransposons, however the elements responsible for genome size variation are yet to be identified. Genomic resources for Pinus including those reported here should assist in further defining whether and how the roles of retrotransposons differ in the evolution of angiosperm and gymnosperm genomes.

Serpentine soil, which is naturally high in heavy metal content and has low calcium to magnesium ratios, comprises a difficult environment for most plants. An impressive number of species are endemic to serpentine, and a wide range of non-endemic plant taxa have been shown to be locally adapted to these soils. Locating genomic polymorphisms which are differentiated between serpentine and non-serpentine populations would provide candidate loci for serpentine adaptation. We have used the Arabidopsis thaliana tiling array, which has 2.85 million probes throughout the genome, to measure genetic differentiation between populations of Arabidopsis lyrata growing on granitic soils and those growing on serpentinic soils. The significant overrepresentation of genes involved in ion transport and other functions provides a starting point for investigating the molecular basis of adaptation to soil ion content, water retention, and other ecologically and economically important variables. One gene in particular, calcium-exchanger 7, appears to be an excellent candidate gene for adaptation to low CaratioMg ratio in A. lyrata.

Recently diverged species typically have incomplete reproductive barriers, allowing introgression of genetic material from one species into the genomic background of the other. The role of natural selection in preventing or promoting introgression remains contentious. Because of genomic co-adaptation, some chromosomal fragments are expected to be selected against in the new background and resist introgression. In contrast, natural selection should favor introgression for alleles at genes evolving under multi-allelic balancing selection, such as the MHC in vertebrates, disease resistance, or self-incompatibility genes in plants. Here, we test the prediction that negative, frequency-dependent selection on alleles at the multi-allelic gene controlling pistil self-incompatibility specificity in two closely related species, Arabidopsis halleri and A. lyrata, caused introgression at this locus at a higher rate than the genomic background. Polymorphism at this gene is largely shared, and we have identified 18 pairs of S-alleles that are only slightly divergent between the two species. For these pairs of S-alleles, divergence at four-fold degenerate sites (K = 0.0193) is about four times lower than the genomic background (K = 0.0743). We demonstrate that this difference cannot be explained by differences in effective population size between the two types of loci. Rather, our data are most consistent with a five-fold increase of introgression rates for S-alleles as compared to the genomic background, making this study the first documented example of adaptive introgression facilitated by balancing selection. We suggest that this process plays an important role in the maintenance of high allelic diversity and divergence at the S-locus in flowering plant families. Because genes under balancing selection are expected to be among the last to stop introgressing, their comparison in closely related species provides a lower-bound estimate of the time since the species stopped forming fertile hybrids, thereby complementing the average portrait of divergence between species provided by genomic data., Author Summary The role of natural selection in promoting or preventing genomic divergence between nascent species remains highly debated. As long as reproductive barriers remain incomplete, genetic material from one species is indeed exposed to natural selection into the genomic background of the other species. In some cases, genomic co-adaptations developing independently in each species are believed to select against such transfers. Yet, theory predicts that the transfer of some chromosomal fragments may be favored by natural selection. In particular, this should occur for alleles at genes evolving under a particular form of natural selection, i.e., multi-allelic balancing selection. We test this prediction using two closely related Arabidopsis species, and find a four-fold lower divergence at alleles at the gene controlling pistil self-incompatibility specificity than at the genomic background. We conclude that alleles at this gene have been transferred more readily between the two species than the genomic background. We suggest that natural selection may efficiently allow the maintenance of high allelic diversity and divergence across many species at S-loci as well as at all other loci under multi-allelic balancing selection, such as the MHC in vertebrates or disease resistance genes in plants.

BackgroundSynonymous DNA substitution rates in the plant chloroplast genome are generally relatively slow and lineage dependent. Non-synonymous rates are usually even slower due to purifying selection acting on the genes. Positive selection is expected to speed up non-synonymous substitution rates, whereas synonymous rates are expected to be unaffected. Until recently, positive selection has seldom been observed in chloroplast genes, and large-scale structural rearrangements leading to gene duplications are hitherto supposed to be rare.Methodology/principle findingsWe found high substitution rates in the exons of the plastid clpP1 gene in Oenothera (the Evening Primrose family) and three separate lineages in the tribe Sileneae (Caryophyllaceae, the Carnation family). Introns have been lost in some of the lineages, but where present, the intron sequences have substitution rates similar to those found in other introns of their genomes. The elevated substitution rates of clpP1 are associated with statistically significant whole-gene positive selection in three branches of the phylogeny. In two of the lineages we found multiple copies of the gene. Neighboring genes present in the duplicated fragments do not show signs of elevated substitution rates or positive selection. Although non-synonymous substitutions account for most of the increase in substitution rates, synonymous rates are also markedly elevated in some lineages. Whereas plant clpP1 genes experiencing negative (purifying) selection are characterized by having very conserved lengths, genes under positive selection often have large insertions of more or less repetitive amino acid sequence motifs.Conclusions/significanceWe found positive selection of the clpP1 gene in various plant lineages to correlated with repeated duplication of the clpP1 gene and surrounding regions, repetitive amino acid sequences, and increase in synonymous substitution rates. The present study sheds light on the controversial issue of whether negative or positive selection is to be expected after gene duplications by providing evidence for the latter alternative. The observed increase in synonymous substitution rates in some of the lineages indicates that the detection of positive selection may be obscured under such circumstances. Future studies are required to explore the functional significance of the large inserted repeated amino acid motifs, as well as the possibility that synonymous substitution rates may be affected by positive selection.