Towards understanding the dynamics of hybridization and apomixis in the evolution of the genus Boechera (Brassicaceae)

We have reviewed past and recent research on the radiation of the North American genus Boechera (Brassicaceae), and discuss future prospects for our understanding of the evolutionary processes and patterns of differentiation in this highly polymorphic genus. Boechera comprises about 80 biennial to perennial species of exclusively North American and Greenlandic distribution. Hybridization and apomixis, in conjunction with polyploidy, were found to play a major role in the origin and maintenance of intraspecific and interspecific polymorphism. However, taxonomic classification within the genus has mainly been based on morphological characters, and diagnostic features discriminating between taxa are often limited. Molecular markers are promising tools for: the assessment of infrageneric phylogen‐etic relationships, the isolation of basal sexual taxa, the detection of hybridization and estimation of its frequency among lineages in time and space, the identification of centres of genetic and species diversity, and the comparison of various reproductive modes. We present and discuss the major achievements obtained since the introduction of molecular methods to Boechera research with respect to these aims, and point out possible shortcomings of specific marker systems. The distribution and transmission of apomixis, polyphyletic origin of genotypes and cytotypes, variation in reproductive system, and the lack of consensus between current taxonomic concepts and evolutionary evidence are the major topics discussed here in the context of hybridization.     

Multiple hybrid formation in natural populations: concerted evolution of the internal transcribed spacer of nuclear ribosomal DNA (ITS) in North American Arabis divaricarpa (Brassicaceae)

DNA sequence variation of the internal transcribed spacer (ITS) region of nuclear ribosomal DNA from Arabis holboellii, A. drummondii, and its putative hybrid A. divaricarpa was analyzed to study hybrid speciation in a species system geographically covering nearly the entire North American continent. Based on molecular systematics the investigated species are better combined under the genus Boechera. Multiple intraindividual ITS copies were detected in numerous accessions of A. divaricarpa, and, to a minor extent, in the parental taxa. Comparative phylogenetic analysis demonstrates that reticulate evolution is common. Consequently, concerted evolution of ITS regions resulted in different types of ITS fragments not only in hybrid populations but also in one of the parental taxa, A. holboellii. Hybrid formation often occurred independently at different sites and at different times, which is reflected by ITS copies resampling the original parental sequence variation in different ways. Some biogeographic structuring of genetic diversity is apparent and mirrors postglacial migration routes. Hybridization, reticulation, and apomixis are assumed to be the major forces driving speciation processes in this species complex. Analysis of conserved regions and secondary structures of the ITS region provided no evidence that, in this system, hybrid ITS evolution is predominantly driven in a particular direction. However, two regions in the ITS1 and ITS2, respectively, show higher mutation rates than expected from outgroup comparisons. Strong evidence for the occurrence of apomixis in A. holboellii and A. divaricarpa has come from pollen size measurements and estimations of pollen quality, which favor the hypothesis that A. drummondii served as paternal hybridization partner more frequently than A. holboellii.     

Supernetwork identifies multiple events of plastid trnF(GAA) pseudogene evolution in the Brassicaceae

The occurrence of nonfunctional trnF pseudogenes has been rarely described in flowering plants. However, we describe the first large-scale supernetwork for the Brassiccaeae built from gene trees for 5 loci (adh, chs, matK, trnL-F, and ITS) and report multiple independent origins for trnF pseudogenes in crucifers. The duplicated regions of the original trnF gene are comprised of its anticodon domain and several other highly structured motifs not related to the original gene. Length variation of the trnL-F intergenic spacer region in different taxa ranges from 219 to 900 bp as a result of differences in pseudocopy number (1-14). It is speculated that functional constraints favor 2-3 or 5-6 copies, as found in Arabidopsis and Boechera. The phylogenetic distribution of microstructural changes for the trnL-F region supports ancient patterns of divergence in crucifer evolution for some but not all gene loci.     

Asexual reproduction in a close relative of Arabidopsis: a genetic investigation of apomixis in Boechera (Brassicaceae)

Understanding apomixis (asexual reproduction through seeds) is of great interest to both plant breeders and evolutionary biologists. The genus Boechera is an excellent system for studying apomixis because of its close relationship to Arabidopsis, the occurrence of apomixis at the diploid level, and its potentially simple inheritance by transmission of a heterochromatic (Het) chromosome. Diploid sexual Boechera stricta and diploid apomictic Boechera divaricarpa (carrying a Het chromosome) were crossed. Flow cytometry, karyotype analysis, genomic in situ hybridization, pollen staining and seed-production measurements were used to analyse the parents and resulting F1, F2 and selected F3 and test-cross (TC) generations. The F1 plant was a low-fertility triploid that produced a swarm of aneuploid and polyploid F2 progeny. Two of the F2 plants were fertile near-tetraploids, and analysis of their F3 and TC progeny revealed that they were sexual and genomically stabilized. The apomictic phenotype was not transmitted by genetic crossing as a single dominant locus on the Het chromosome, suggesting a complex genetic control of apomixis that has implications for future genetic and evolutionary analyses in this group.     

Genome structure and apomixis in Phoenicaulis (Brassicaceae; Boechereae)

Apomixis in crucifer (Brassicaceae) species is rare, reported in just four genera (Boechera, Draba, Erysimum, and Parrya), and one intergeneric hybrid (Raphanobrassica). It is well studied only in Boechera, where it is widespread among 100+ recognized species. However, its occurrence in related genera of the tribe Boechereae has not been documented previously. Here we analyzed genome evolution, mode of reproduction, and fertility of the monospecific Boechereae genus Phoenicaulis (P. cheiranthoides), endemic to the northwestern United States. We discovered that the species encompasses diploid (2n = 2x = 14), triploid (2n = 3x = 21), and tetraploid (2n = 4x = 28) populations. Comparative chromosome painting analyses revealed that the three karyotypes are essentially structurally identical, differing only in the presence of a largely heterochromatic chromosome (Het) in the triploid and tetraploid cytotypes. The genome structure of Phoenicaulis appeared identical to that of Boechera species previously analyzed, suggesting genomic stasis despite the morphological and molecular divergence of the two genera. This genome colinearity extended to the presence and structure of the Het chromosomes, which are closely associated with apomictic reproduction in Boechera. Interestingly, all three cytotypes of Phoenicaulis proved to be apomictic, regardless of the presence or absence of a Het chromosome, and sexual populations have yet to be identified.     

Genome structure and evolution in the cruciferous tribe Thlaspideae (Brassicaceae)

Whole-genome duplications (WGDs) and chromosome rearrangements (CRs) play the key role in driving the diversification and evolution of plant lineages. Although the direct link between WGDs and plant diversification is well documented, relatively few studies focus on the evolutionary significance of CRs. The cruciferous tribe Thlaspideae represents an ideal model system to address the role of large-scale chromosome alterations in genome evolution, as most Thlaspideae species share the same diploid chromosome number (2n = 2x = 14). Here we constructed the genome structure in 12 Thlaspideae species, including field pennycress (Thlaspi arvense) and garlic mustard (Alliaria petiolata). We detected and precisely characterized genus- and species-specific CRs, mostly pericentric inversions, as the main genome-diversifying drivers in the tribe. We reconstructed the structure of seven chromosomes of an ancestral Thlaspideae genome, identified evolutionary stable chromosomes versus chromosomes prone to CRs, estimated the rate of CRs, and uncovered an allohexaploid origin of garlic mustard from diploid taxa closely related to A. petiolata and Parlatoria cakiloidea. Furthermore, we performed detailed bioinformatic analysis of the Thlaspideae repeatomes, and identified repetitive elements applicable as unique species- and genus-specific barcodes and chromosome landmarks. This study deepens our general understanding of the evolutionary role of CRs, particularly pericentric inversions, in plant genome diversification, and provides a robust base for follow-up whole-genome sequencing efforts.     

Transposable element-medicated evolution of sex: A population genetic model

For a population made up of individuals capable of sexual as well as asexual modes of reproduction, conditions for the spread of a transposable element are explored using a one-locus, two-haplotype model. The analysis is then extended to include the possibility that the transposable element can modulate the probability of sexual reproduction, thus casting Hickey’s (1982,Genetics 101: 519–531) suggestion in a population genetics framework. The model explicitly includes the cost of sexual reproduction, fitness disadvantage to the transposable element, probability of transposition, and the predisposition for sexual reproduction in the presence and absence of the transposable element. The model predicts several kinds of outcome, including initial frequency dependence and stable polymorphism. More importantly, it is seen that for a wide range of parameter values, the transposable element can go to fixation. Therefore it is able to convert the population from a predominantly asexual to a predominantly sexual mode of reproduction. Viewed in conjunction with recent results implicating short stretches of apparently non-coding DNA in sex determination (McCoubreyet al. 1988,Science 242: 1146–1151), the model hints at the important role this mechanism could have played in the evolution of sexuality.     

Evolution of the Eastern Asian and Eastern North American disjunct genus Symplocarpus (Araceae): insights from chloroplast DNA restriction site data

Chloroplast DNA (cpDNA) variation was surveyed with 20 restriction endonucleases for the eastern Asian and eastern North American disjunct genus Symplocarpus (Araceae). The cpDNA phylogeny reveals a sister group relationship between S. foetidus from eastern North America and S. renifolius from eastern Asia. The cpDNA divergence between the two intercontinental sister species is 0.61%, which suggests an estimated divergence time of 6.1 million years ago during the late Miocene. The Bering land bridge hypothesis is compatible with the estimated time of divergence for the migration of Symplocarpus between eastern Asia and North America. Furthermore, a single origin of the exothermic spadices in Symplocarpus is suggested by the phylogeny. The cpDNA data also provide independent support for the recognition of three species within the genus.     

Evolution of the plastid ribosomal RNA operon in a nongreen parasitic plant: Accelerated sequence evolution,altered promoter structure,and tRNA pseudogenes

The nucleotide sequence of a 7.4 kb region containing the entire plastid ribosomal RNA operon of the nongreen parasitic plant Epifagus virginiana has been determined. Analysis of the sequence indicates that all four rRNA genes are intact and almost certainly functional. In contrast, the split genes for tRNA^Ile and tRNA^Ala present in the 16S-23S rRNA spacer region have become pseudogenes, and deletion upstream of the 16S rRNA gene has removed a tRNA^Val gene and most of the promoter region for the rRNA operon. The rate of nucleotide substitution in 16S and 23S rRNAs is several times higher in Epifagus than in tobacco, a related photosynthetic plant. Possible reasons for this, including relaxed translational constraints, are discussed.     

Molecular phylogeny and evolution of the freshwater eel, genus Anguilla

A molecular phylogenetic analysis was conducted on all of the known catadromous eel species of the genus Anguilla to assess their relationships and evolutionary history. The analyses of a total of 1427 bp of the mitochondrial 16S ribosomal RNA and 1140 bp of the complete cytochrome b gene sequences suggested that the genus Anguilla was monophyletic in origin, with A. borneensis as the most basal species. Four clades/species groups that correspond to their geographical ranges were indicated, Indo-Atlantic (three species), Oceania (two species), tropical Pacific (two species), and Indo-Pacific (five species), with ambiguous positions for A. japonica and A. reinhardti. This grouping conflicts with that of a previous morphological study, since the broad undivided maxillary and short-fin type, which were thought to be phylogenetically important, were paraphyletic in the molecular analysis. However, the molecular phylogeny and the present geographic distribution of species suggested historical dispersion of the genus Anguilla according to the Tethys corridor hypothesis, which proposed that anguillid eels originated near present-day Indonesia and dispersed westward along paleo-circumglobal equatorial currents. The westward-moving strain entered the paleo-Atlantic through the Tethys Sea and was ancestral to present-day European and American species.     

Speciation in North American black basses,Micropterus (Actinopterygii: Centrarchidae)

Abstract.— The Pleistocene Epoch has been frequently cited as a period of intense speciation for a significant portion of temperate continental biotas. To critically assess the role of Pleistocene glaciations on the evolution of the freshwater fish clade Micropterus , we use a phylogenetic analysis of complete gene sequences from two mitochondrial genes (cytochrome b and ND2), and a fossil calibration of the molecular clock to estimate ages of speciation events and rates of diversification. The absence of substantial morphological and ecological divergence together with endemism of five of the eight species in North American tributaries of the Gulf of Mexico may be interpreted as the result of a recent Pleistocene origin for these species. Speciation dates in Micropterus range from 1.01 ± 0.32 to 11.17 ± 1.02 million years ago. Only one speciation event is dated to the Pleistocene, and rates of diversification are not significantly variable in Micropterus. The premise that the Pleistocene was an exceptional period of speciation in Micropterus is not supported. Instead, a Gulf Coast allopatric speciation model is proposed, and predicts periods of dynamic speciation driven by sea level fluctuations in the Late Miocene and Pliocene. The Pleistocene, however, was a period of significant intraspecific mitochondrial lineage diversification. The application of the Gulf Coast allopatric speciation model to the remaining aquatic fauna of the Gulf of Mexico coast in North America will rely on robust phylogenetic hypotheses and accurate age estimations of speciation events.     

Transposable elements reveal the impact of introgression,rather than transposition,in Pisum diversity,evolution, and domestication

The genetic structure and evolutionary history of the genus Pisum were studied exploiting our germplasm collection to compare the contribution of different mechanisms to the generation of diversity. We used sequence-specific amplification polymorphism (SSAP) markers to assess insertion site polymorphism generated by a representative of each of the two major groups of LTR-containing retrotransposons, PDR1 (Ty1/copia-like) and Cyclops (Ty3/gypsy-like), together with Pis1, a member of the En/Spm transposon superfamily. The analysis of extended sets of the four main Pisum species, P. fulvum, P. elatius, P. abyssinicum, and P. sativum, together with the reference set, revealed a distinct pattern of the NJ (Neighbor-Joining) tree for each basic lineage, which reflects the different evolutionary history of each species. The SSAP markers showed that Pisum is exceptionally polymorphic for an inbreeding species. The patterns of phylogenetic relationships deduced from different transposable elements were in general agreement. The retrotransposon-derived markers gave a clearer separation of the main lineages than the Pis1 markers and were able to distinguish the truly wild form of P. elatius from the antecedents of P. sativum. There were more species-specific and unique PDR1 markers than Pis1 markers in P. fulvum and P. elatius, pointing to PDR1 activity during speciation and diversification, but the proportion of these markers is low. The overall genetic diversity of Pisum and the extreme polymorphism in all species, except P. abyssinicum, indicate a high contribution of recombination between multiple ancestral lineages compared to transposition within lineages. The two independently domesticated pea species, P. abyssinicum and P. sativum, arose in contrasting ways from the common processes of hybridization, introgression, and selection without associated transpositional activity.     

Molecular evolution of cdc2 pseudogenes in spruce (Picea)

The p34^cdc2 protein and other cyclin-dependent protein kinases (CDK) are important regulators of eukaryotic cell cycle progression. We have previously cloned a functional cdc2 gene from Picea abies and found it to be part of a family of related sequences, largely consisting of pseudogenes. We now report on the isolation of partial cdc2 pseudogenes from Picea engelmannii and Picea sitchensis, as well as partial functional cdc2 sequences from P. engelmannii, P. sitchensis and Pinus contorta. A high level of conservation between species was detected for these sequences. Phylogenetic analyses of pseudogene and functional cdc2 sequences, as well as the presence of shared insertions or deletions, support the division of most of the cdc2 pseudogenes into two subfamilies. New cdc2 pseudogenes appear to have been formed in Picea at a much higher rate than they have been obliterated by neutral mutations. The pattern of nucleotide changes in the cdc2 pseudogenes, as compared to a presumed ancestral functional cdc2 gene, was similar to that previously found in mammalian pseudogenes, with a strong bias for the transitions C to T and G to A, and the transversions C to A and G to T.     

Plant transposable elements: Their role in evolution

Transposable elements (TE) are natural constituents of plant genomes. However, their presence only becomes apparent if they become dislodged from their resident positions in the genome and transpore into another gene, thereby inducing a mutation. Such TE-induced mutations are somatically unstable because they revert to wild type and hence reconstitute the expression of the mutated gene. The frequent somatic excision of the TE results in a variegated phenotype. Since this instability is inherited in a Mendelian manner the variegated phenotype is nuclear determined. By this criterion TE have been shown to occur in more than 30 species belonging to different families and genera. Many questions arise when dealing with TE: their structure and functions, and the biological significance of the activity of elements in the differentiation of a normal plant or in the evolution of plant genes.     

An updated genus‐wide phylogenetic analysis of Arisaema (Araceae) with reference to sections

Arisaema has a long and complicated taxonomic history regarding its infrageneric classification. In the latest system, 14 sections were recognized, based on an unpublished, tentative phylogenetic analysis; in addition the type species and nomenclatural priority for each section were confirmed. Here, we present an updated, genus‐wide phylogenetic analysis, based on four plastid non‐coding regions (3′trnL–trnF, rpl20–5′rps12, psbB–psbH and rpoC2–rps2) for > 150 accessions. The maximum parsimony and maximum likelihood phylogenetic analyses identified eight major clades and one branch with unique sequence variation, although the relationships were unclear due to a polytomy and weak support. In the phylogenetic trees, most of the sections proposed in the latest system were distinct and corresponded to the major clades, but some sections are not monophyletic. On the basis of the phylogenetic relationship: (1) A. schimperianum is treated as a member of section Arisaema, rather than section Tenuipistillata or section Sinarisaema; and (2) section Fimbriata, which was synonymised into section Attenuata in the latest system but is morphologically distinct from the other species, is redefined as a monotypic section. In conclusion, we recognize 15 sections of Arisaema and species‐level classifications are discussed in a phylogenetic context.     

Phylogenetic analysis of the genus Hexachlamys (Myrtaceae) based on plastid and nuclear DNA sequences and their taxonomic implications

Myrtaceae are one of the most species‐rich families of flowering plants in the Neotropics. They include several complex genera and species; Hexachlamys is one of the complex genera. It has not been recognized as a distinct genus and has been included in Eugenia, based on morphological grounds. Therefore, molecular systematic studies may be useful to understand and to help to solve these relationships. Here, we performed a molecular phylogenetic analysis using plastid and nuclear data in order to check the inclusion of Hexachlamys in Eugenia. Plastid (accD, rpoB, rpoC1, trnH‐psbA) and nuclear (ITS2) sequence data were analysed using Bayesian and maximum parsimony methods. The trees constructed using ITS2 and trnH‐psbA were the best able to resolve the relationships between species and genera, revealing the non‐monophyly of Hexachlamys. The molecular phylogenetic analyses were in agreement with previous morphological revisions that have included Hexachlamys in Eugenia. These results reinforce the importance of uniting knowledge and strategies to understand better issues of delimitation of genera and species in groups of plants with taxonomic problems. © 2013 The Linnean Society of London, Botanical Journal of the Linnean Society, 2013, 172 , 532–543.     

Genetic discontinuity revealed by chloroplast microsatellites in eastern North American Abies (Pinaceae)

Development of conservation strategies for Fraser fir (Abies fraseri) in the southern Appalachian Mountains depends in part on recognition of the extent to which Fraser fir is genetically distinct from the closely related balsam (A. balsamea) and intermediate (A. balsamea var. phanerolepis) fir. These sibling species have exhibited intergrading, clinal variation in morphological, chemical, and genetic characteristics in prior research. Chloroplast microsatellite markers were polymerase chain reaction amplified from genomic DNA samples of 78 individuals representing the geographic ranges of Fraser, balsam, and intermediate fir. Gene diversity levels at two loci ranged among taxa from 0.65 to 0.84. Allele frequencies demonstrated significant differentiation among taxa, with R(ST) values of 0.36 and 0.10. Haplotype diversity and D(SH) were highest for balsam fir and lowest for intermediate fir. A haplotype network analysis based on allele size distribution for the two loci revealed two distinct clusters of haplotypes and population-specific haplotypes. Ninety-two percent of the haplotypes in one cluster were from balsam fir and intermediate fir, and 84% of the haplotypes in the other cluster were from Fraser fir and intermediate fir. The genetic differentiation of chloroplast DNA markers provides justification for the recognition of Fraser fir as a distinct Management Unit (MU) for conservation purposes, regardless of its taxonomic classification.