The evolutionary fate of the chloroplast and nuclear rps16 genes as revealed through the sequencing and comparative analyses of four novel legume chloroplast genomes from Lupinus.

The Fabaceae family is considered as a model system for understanding chloroplast genome evolution due to the presence of extensive structural rearrangements, gene losses and localized hypermutable regions. Here, we provide sequences of four chloroplast genomes from the Lupinus genus, belonging to the underinvestigated Genistoid clade. Notably, we found in Lupinus species the functional loss of the essential rps16 gene, which was most likely replaced by the nuclear rps16 gene that encodes chloroplast and mitochondrion targeted RPS16 proteins. To study the evolutionary fate of the rps16 gene, we explored all available plant chloroplast, mitochondrial and nuclear genomes. Whereas no plant mitochondrial genomes carry an rps16 gene, many plants still have a functional nuclear and chloroplast rps16 gene. Ka/Ks ratios revealed that both chloroplast and nuclear rps16 copies were under purifying selection. However, due to the dual targeting of the nuclear rps16 gene product and the absence of a mitochondrial copy, the chloroplast gene may be lost. We also performed comparative analyses of lupine plastomes (SNPs, indels and repeat elements), identified the most variable regions and examined their phylogenetic utility. The markers identified here will help to reveal the evolutionary history of lupines, Genistoids and closely related clades.

The chloroplast genome of the hexaploid Spartina maritima (Poaceae, Chloridoideae): Comparative analyses and molecular dating.

The history of many plant lineages is complicated by reticulate evolution with cases of hybridization often followed by genome duplication (allopolyploidy). In such a context, the inference of phylogenetic relationships and biogeographic scenarios based on molecular data is easier using haploid markers like chloroplast genome sequences. Hybridization and polyploidization occurred recurrently in the genus Spartina (Poaceae, Chloridoideae), as illustrated by the recent formation of the invasive allododecaploid S. anglica during the 19th century in Europe. Until now, only a few plastid markers were available to explore the history of this genus and their low variability limited the resolution of species relationships. We sequenced the complete chloroplast genome (plastome) of S. maritima, the native European parent of S. anglica, and compared it to the plastomes of other Poaceae. Our analysis revealed the presence of fast-evolving regions of potential taxonomic, phylogeographic and phylogenetic utility at various levels within the Poaceae family. Using secondary calibrations, we show that the tetraploid and hexaploid lineages of Spartina diverged 6-10 my ago, and that the two parents of the invasive allopolyploid S. anglica separated 2-4 my ago via long distance dispersal of the ancestor of S. maritima over the Atlantic Ocean. Finally, we discuss the meaning of divergence times between chloroplast genomes in the context of reticulate evolution.

Exploring the genome of the salt-marsh Spartina maritima (Poaceae, Chloridoideae) through BAC end sequence analysis.

Spartina species play an important ecological role on salt marshes. Spartina maritima is an Old-World species distributed along the European and North-African Atlantic coasts. This hexaploid species (2n = 6x = 60, 2C = 3,700 Mb) hybridized with different Spartina species introduced from the American coasts, which resulted in the formation of new invasive hybrids and allopolyploids. Thus, S. maritima raises evolutionary and ecological interests. However, genomic information is dramatically lacking in this genus. In an effort to develop genomic resources, we analysed 40,641 high-quality bacterial artificial chromosome-end sequences (BESs), representing 26.7 Mb of the S. maritima genome. BESs were searched for sequence homology against known databases. A fraction of 16.91% of the BESs represents known repeats including a majority of long terminal repeat (LTR) retrotransposons (13.67%). Non-LTR retrotransposons represent 0.75%, DNA transposons 0.99%, whereas small RNA, simple repeats and low-complexity sequences account for 1.38% of the analysed BESs. In addition, 4,285 simple sequence repeats were detected. Using the coding sequence database of Sorghum bicolor, 6,809 BESs found homology accounting for 17.1% of all BESs. Comparative genomics with related genera reveals that the microsynteny is better conserved with S. bicolor compared to other sequenced Poaceae, where 37.6% of the paired matching BESs are correctly orientated on the chromosomes. We did not observe large macrosyntenic rearrangements using the mapping strategy employed. However, some regions appeared to have experienced rearrangements when comparing Spartina to Sorghum and to Oryza. This work represents the first overview of S. maritima genome regarding the respective coding and repetitive components. The syntenic relationships with other grass genomes examined here help clarifying evolution in Poaceae, S. maritima being a part of the poorly-known Chloridoideae sub-family.

Transcriptome de novo assembly from next-generation sequencing and comparative analyses in the hexaploid salt marsh species Spartina maritima and Spartina alterniflora (Poaceae).

Spartina species have a critical ecological role in salt marshes and represent an excellent system to investigate recurrent polyploid speciation. Using the 454 GS-FLX pyrosequencer, we assembled and annotated the first reference transcriptome (from roots and leaves) for two related hexaploid Spartina species that hybridize in Western Europe, the East American invasive Spartina alterniflora and the Euro-African S. maritima. The de novo read assembly generated 38 478 consensus sequences and 99% found an annotation using Poaceae databases, representing a total of 16 753 non-redundant genes. Spartina expressed sequence tags were mapped onto the Sorghum bicolor genome, where they were distributed among the subtelomeric arms of the 10 S. bicolor chromosomes, with high gene density correlation. Normalization of the complementary DNA library improved the number of annotated genes. Ecologically relevant genes were identified among GO biological function categories in salt and heavy metal stress response, C4 photosynthesis and in lignin and cellulose metabolism. Expression of some of these genes had been found to be altered by hybridization and genome duplication in a previous microarray-based study in Spartina. As these species are hexaploid, up to three duplicated homoeologs may be expected per locus. When analyzing sequence polymorphism at four different loci in S. maritima and S. alterniflora, we found up to four haplotypes per locus, suggesting the presence of two expressed homoeologous sequences with one or two allelic variants each. This reference transcriptome will allow analysis of specific Spartina genes of ecological or evolutionary interest, estimation of homoeologous gene expression variation using RNA-seq and further gene expression evolution analyses in natural populations.

Transcriptome divergence between the hexaploid salt-marsh sister species Spartina maritima and Spartina alterniflora (Poaceae).

Invasive species are ideal model systems to investigate the evolutionary processes associated with their ecological success by comparison with closely related species. In this article, we explore transcriptome evolution following divergence between two closely related salt-marsh species, the invasive Spartina alterniflora (native to the East-American Atlantic coast, introduced in several continents) and the declining Spartina maritima (native to the Euro-African Atlantic coast). We have explored the utility of cross-species hybridization microarrays using rice (Oryza sativa) oligo-microarrays to compare leaf expression patterns between these species. Coding sequence comparisons from 10 nuclear genes (2256 bp) revealed that nucleotide divergence between Spartina and Oryza range from 8% to 14%. More than 70% of the 60-mer oligonucleotide sequences spotted on the rice microarray exhibited stable and repeatable patterns when hybridized against Spartina RNA. In total, 9353 (44.5%) genes on the array hybridized with both species S. maritima and S. alterniflora. Among these genes, 1247 genes were found to be differentially expressed between the two Spartina species, most of them (957) being up-regulated in S. alterniflora. In particular, developmental and cellular growth genes (gene ontology, biological process) were highly up-regulated in S. alterniflora and down-regulated in S. maritima, whereas genes involved in stress response were up-regulated in S. maritima. Our findings indicate the suitability of cross-species microarray hybridization between Spartina and O. sativa and reveal the extent of leaf transcriptome evolution that took place during the divergence between S. alterniflora and S. maritima. Expression patterns are consistent with the morphological differentiation and differential expansion of the two species.

Making a functional diploid: from polysomic to disomic inheritance.

One little understood feature of polyploid speciation is the transition from polysomic to disomic inheritance, and much recent attention has focused on the role of pairing genes in this process. Using computer simulations we studied the effects of mutations, chromosomal inversions, chiasma, neofunctionalization, subfunctionalization and selection on the evolution of disomic inheritance in a polyploid over 10 000 generations. We show that: the evolution of pairing genes is not essential for the establishment of disomic inheritance, as genetic drift, coupled with a threshold for homologue pairing fidelity, is sufficient to explain the transition from polysomic to disomic inheritance; high rates of recombination increase the number of generations required for disomic inheritance to become established; both neofunctionalization and subfunctionalization speed up the transition to disomic inheritance. The data suggest that during polyploid species establishment, selection will favour reduced chiasma number and/or more focused distribution. The data also suggest a new role for subfunctionalization in that it can drive disomic inheritance. The evolution of subfunctionalization in genes across the genome will then act to maintain genes in syntenic blocks and may explain why such regions are so highly conserved.

Tracking the evolutionary history of polyploidy in Fragaria L. (strawberry): new insights from phylogenetic analyses of low-copy nuclear genes.

Phylogenetic utility of two nuclear genes (GBSSI-2 and DHAR) was explored in genus Fragaria in order to clarify phylogenetic relationships among taxa and to elucidate the origin of the polyploid species. Orthology of the amplified products was assessed by several methods. Our results strongly suggest the loss of one GBSSI duplicated copy (GBSSI-1) in the Fragariinae subtribe. Phylogenetic analyses provided new insights into the evolutionary history of Fragaria, such as evidence supporting the presence of three main diploid genomic pools in the genus and demonstrating the occurrence of independent events of polyploidisation. In addition, the data provide evidence supporting an allopolyploid origin of the hexaploid F. moschata, and the octoploids F. chiloensis, F. iturupensis and F. virginiana. Accordingly, a new pattern summarizing our present knowledge on the Fragaria evolutionary history is proposed. Additionally, sequence analyses also revealed relaxed constraints on homoeologous copies at high ploidy level, as demonstrated by deletion events within DHAR coding sequences of some allo-octoploid haplotypes.

The enigmatic invasive Spartina densiflora: a history of hybridizations in a polyploidy context.

The aim of this study was to explore the origin of the invasive Spartina densiflora by analysing samples from the native region (South America) and from a recently colonized area (California). A combination of various molecular data (chloroplast and nuclear sequences, molecular fingerprint) and ploidy level estimations was used to answer the question whether the reticulate phylogenetic pattern previously detected in S. densiflora was restricted to California, or alternatively, whether a more ancient hybrid origin preceded formation of this species in its native area. We found that this species is heptaploid in both its native and introduced range. Identification of nuclear homeologous sequences indicate that this species has a reticulate origin in its native range, involving a lineage related to the hexaploid clade formed by S. alterniflora, S.foliosa, and S. maritima, and another lineage related to the sub-Antarctic endemic S. arundinacea that provided the chloroplast genome. The samples from California displayed similar multilocus patterns to the samples from Chile, supporting the hypothesis that this species originated on the southeast American coast (Argentina), from where it eventually spread to the west coast of South America (Chile) first and to the Northern Hemisphere (California) later.

Evolutionary dynamics of Waxy and the origin of hexaploid Spartina species (Poaceae).

We investigated the evolutionary dynamics of duplicated copies of the granule-bound starch synthase I gene (GBSSI or Waxy) within polyploid Spartina species. Molecular cloning, sequencing, and phylogenetic analyses revealed incongruences between the expected species phylogeny and the inferred gene trees. Some genes within species were more divergent than expected from ploidy level alone, suggesting the existence of paralogous sets of Waxy loci in Spartina. Phylogenetic analyses indicate that this paralogy originated from a duplication that occurred prior to the divergence of Spartina from other Chloridoideae. Gene tree topologies revealed three divergent homoeologous sequences in the hexaploid S. alterniflora that are consistent with the proposal of an allopolyploid origin of the hexaploid clade. Waxy sequences differ in insertion-deletion events in introns, which may be used to diagnose gene copies. Both paralogous and homoeologous coding regions appear to evolving under selective constraints.

Uniformity of the nuclear and chloroplast genomes of Spartina maritima (Poaceae), a salt-marsh species in decline along the Western European Coast.

Spartina maritima is a salt-marsh species from European and African Atlantic coasts. In the northern range of the species (including north-west France), a rapid decline of the populations has been observed during the 20th century. In this paper, the molecular diversity of 10 populations of S. maritima from France has been investigated using nuclear and chloroplast DNA markers: inter-simple sequence polymorphism (ISSR), randomly amplified polymorphic DNA (RAPD), inter-retrotransposon amplified polymorphism (IRAP), and PCR-RFLP of a 5 kb long portion of chloroplast DNA. The results reveal an extremely low level of genetic variation in this species: only one nuclear marker (out of 98) was polymorphic, with the presence of two genotypes randomly distributed among the populations. The lack of genetic diversity is interpreted in light of the almost exclusive vegetative propagation of the species in its northern range, the colonization history of the populations, and the origin of S. maritima (2n = 60) in the hexaploid lineage of the genus and in the context of the management of S. maritima populations.

Molecular phylogeny of hybridizing species from the genus Spartina Schreb. (Poaceae).

Interspecific hybridization events have been reported in the genus Spartina Schreb. (Poaceae), involving the east American species Spartina alterniflora, and including either introgression (e.g., with the western American Spartina foliosa) or allopolyploid speciation (e.g., with the Euro-African Spartina maritima). Molecular phylogenetic analysis of the genus has been undertaken in order to understand phylogenetic relationships and genetic divergence among these hybridizing species. Twelve Spartina species have been sequenced for two nuclear DNA regions (ITS of ribosomal DNA, and part of the Waxy gene) and one chloroplast DNA spacer (trnT-trnL). Separate and conditional combined phylogenetic analyses using Cynodon dactylon as the outgroup have been conducted. Spartina is composed of two lineages. The first clade includes all hexaploid species: the Euro-African S. maritima (2n = 60), the East-American S. alterniflora (2n = 62) and the West-American S. foliosa (2n = 60). Spartina alterniflora appears as a closely related sister species to S. foliosa. Although belonging to the same lineage, Spartina maritima appears consistently more genetically differentiated from S. alterniflora than S. foliosa. The tetraploid species S. argentinensis (2n = 40) is placed at the base of this first clade according to the Waxy data, but its position is not well resolved by the other sequences. The second well-supported main lineage within genus Spartina includes the other tetraploid American species. Significant incongruence has been encountered between the waxy based tree and both the ITS and trnT-trnL trees concerning the position of S. densiflora, suggesting a possible reticulate evolution for this species. The results agree with hybridization patterns occurring in Spartina: introgression involving closely related species (S. alterniflora and S. foliosa) on one hand, and alloploid speciation involving more differentiated species (S. alterniflora and S. maritima) on the other hand.

Molecular investigations in populations of Spartina anglica C.E. Hubbard (Poaceae) invading coastal Brittany (France).

Spartina anglica is a classical example of recent alloploid speciation. It arose during the end of the nineteenth century in England by hybridization between the indigenous Spartina maritima and the introduced East-American Spartina alterniflora. Duplication of the hybrid genome (Spartina x townsendii) gave rise to a vigorous allopolyploid involved in natural and artificial invasions on different continents. Spartina anglica was first recorded in France in 1906, and since then, it has spread all along the western French coast. Earlier studies revealed that native British populations display consistent morphological plasticity and lack of isozyme variation. In this paper, we use different molecular markers (randomly amplified polymorphic DNA, intersimple sequence repeats and restriction patterns from nuclear and chloroplast DNA sequences) to analyse the genetic patterns of the French populations of S. anglica. Our results show that French populations are mainly composed of one "major" multilocus genotype. This genotype is identical to the first-generation hybrid S. x townsendii from England. Losses of few markers from this genotype are observed but are restricted to a few populations from Brittany; it is likely that they appeared independently, subsequent to their introduction. In southern Brittany, no hybrids between S. anglica and S. maritima have been found where the two species co-occur. All French populations of S. anglica display the same chloroplast DNA sequences as S. alterniflora, the maternal genome donor. These findings are consistent with a severe genetic bottleneck at the time of the species formation, as a consequence of a unique origin of the species. Both parental nuclear sequences are present in the allopolyploid populations, revealing that for the markers investigated, no extensive changes have occurred in this young species.

On the origins of the tetraploid Bromus species (section Bromus, Poaceae): insights from internal transcribed spacer sequences of nuclear ribosomal DNA.

The internal transcribed spacer (ITS) region of nuclear ribosomal DNA from 22 diploid and tetraploid annual Bromus species of section Bromus (Poaceae) and three species belonging to other Bromus sections, Bromus catharticus (section Ceratochloa), Bromus anomalus (section Pnigma), and Bromus sterilis (section Genea), were investigated by PCR amplification and direct sequencing. The length of the ITS-1 region varied from 215 to 218 bp, and that of the ITS-2 region from 215 to 216 bp, in the species analyzed. ITS-1 was more variable and provided more informative sites (49) than ITS-2 (32). No variation was encountered within species. In pairwise comparison among species of section Bromus, sequence divergence ranged from 0.0 to 8.0% for the combined ITS-1 and ITS-2 regions. Parsimony analysis using Avena longiglumis and Hordeum vulgare as outgroups resulted in well-resolved phylogenetic trees and showed that section Bromus is monophyletic according to the species analyzed outside of the section. The analysis clarified the phylogenetic relationships among monogenomic (diploid) species. Introduction of the allotetraploid species did not change the general topology of the trees obtained using only the diploid species. Although some tetraploid-diploid species relationships will have to be clarified with faster evolving markers, the ITS sequences are shown to be useful for assessing evolutionary relationships among closely related Bromus species, as well as for clarifying taxonomic problems in previously controversial cases (e.g., Bromus alopecuros and Bromus caroli-henrici). New hypotheses are proposed concerning the origin of several allotetraploid species. For example, it is shown that the tetraploid Bromus hordeaceus diverged earlier than all other species of section Bromus, excluding the diploid B. caroli-henrici, which is found to be basal in this group. The tetraploid Bromus arenarius, which was considered a hybrid between sections Bromus and Genea, and the tetraploid Bromus adoensis are sister taxa within section Bromus; they belong in a weakly differentiated clade with the diploids Bromus brachystachys, Bromus japonicus, Bromus squarrosus, Bromus arvensis, and Bromus intermedius.

Genetic diversity in Mediterranean diploid and tetraploid Bromus L. (section Bromus Sm.) populations.

The levels of genetic diversity assessed from allozyme data were investigated in 25 populations of Mediterranean Bromus intermedius, B. squarrosus, B. lanceolatus, and B. hordeaceus from Algeria. The geographically restricted diploids B. intermedius and B. squarrosus displayed less genetic diversity (the mean population gene diversity of Nei (Hu) ranged from 0.03 to 0.12) than the widespread tetraploid colonizers B. lanceolatus and B. hordeaceus (Hu = 0.07-0.27). Deviations from Hardy-Weinberg expectations in diploid populations of B. intermedius and B. squarrosus were observed owing to heterozygote excess at several loci and suggested that these self-fertilizing species may have substantial amounts of allogamy. Tetraploid populations of B. lanceolatus and B. hordeaceus were largely homozygous at homologous loci and frequently exhibited intergenomic fixed heterozygosity in accordance with their alloploid origin. Genetic variation at the infraspecific level was mostly distributed within populations in the four species, B. hordeaceus showing the lowest level of interpopulation differentiation (Gst = 0.06) and the highest level of gene flow (Nm = 3.75). Consistent gene flows are in agreement with the strongest intercontinental invasive behaviour of B. hordeaceus. Less differentiation was reported in the literature among later introduced B. hordeaceus populations from England and Australia, indicating reduced differentiation under the process of colonization. Moderate divergence occured among the four taxa, with interspecific genetic identities ranging from 0.87 to 0.93. In spite of substantial genetic similarity, species were clearly differentiated, with each tetraploid being more closely related to a diploid: B. hordeaceus to B. squarrosus and B. lanceolatus to B. intermedius.