Evolutionary rates of commonly used nuclear and organelle markers of Arabidopsis relatives (Brassicaceae)

The complex polyploid genome architecture of sugarcane

Sugarcane, the world's most harvested crop by tonnage, has shaped global history, trade and geopolitics, and is currently responsible for 80% of sugar production worldwide1. While traditional sugarcane breeding methods have effectively generated cultivars adapted to new environments and pathogens, sugar yield improvements have recently plateaued2. The cessation of yield gains may be due to limited genetic diversity within breeding populations, long breeding cycles and the complexity of its genome, the latter preventing breeders from taking advantage of the recent explosion of whole-genome sequencing that has benefited many other crops. Thus, modern sugarcane hybrids are the last remaining major crop without a reference-quality genome. Here we take a major step towards advancing sugarcane biotechnology by generating a polyploid reference genome for R570, a typical modern cultivar derived from interspecific hybridization between the domesticated species (Saccharum officinarum) and the wild species (Saccharum spontaneum). In contrast to the existing single haplotype ('monoploid') representation of R570, our 8.7 billion base assembly contains a complete representation of unique DNA sequences across the approximately 12 chromosome copies in this polyploid genome. Using this highly contiguous genome assembly, we filled a previously unsized gap within an R570 physical genetic map to describe the likely causal genes underlying the single-copy Bru1 brown rust resistance locus. This polyploid genome assembly with fine-grain descriptions of genome architecture and molecular targets for biotechnology will help accelerate molecular and transgenic breeding and adaptation of sugarcane to future environmental conditions.

A taxonomic revision of the genus Angelica (Apiaceae) in Taiwan with a new species A. aliensis

BACKGROUND

Angelica L. sensu lato is a taxonomically complex genus, and many studies have utilized morphological and molecular features to resolve its classification issues. In Taiwan, there are six taxa within Angelica, and their taxonomic treatments have been a subject of controversy. In this study, we conducted a comprehensive analysis incorporating morphological and molecular (cpDNA and nrDNA) characteristics to revise the taxonomic treatments of Angelica in Taiwan.

RESULTS

As a result of our research, we have revised the classification between A. dahurica var. formosana and A. pubescens and merged two varieties of A. morrisonicola into a single taxon. A new taxon, A. aliensis, has been identified and found to share a close relationship with A. tarokoensis. Based on the morphological and molecular characteristics data, it has been determined that the former three taxa should be grouped into the Eurasian Angelica clade, while the remaining four taxa should belong to the littoral Angelica clade. Furthermore, Angelica species in Taiwan distributed at higher altitudes displayed higher genetic diversity, implying that the central mountain range of Taiwan serves as a significant reservoir of plant biodiversity. Genetic drift, such as bottlenecks, has been identified as a potential factor leading to the fixation or reduction of genetic diversity of populations in most Angelica species. We provide key to taxa, synopsis, phenology, and distribution for each taxon of Taiwan.

CONCLUSIONS

Our comprehensive analysis of morphological and molecular features has shed light on the taxonomic complexities within Angelica in Taiwan, resolving taxonomic issues and providing valuable insights into the phylogenetic relationships of Angelica in Taiwan.

Effect of geographic isolation on genetic variation and population structure of Euphrasia nankotaizanensis, a threatened endemic alpine herb in Taiwan

Euphrasia nankotaizanensis (Orobanchaceae) is a rare alpine herb that is endemic to Taiwan. Only four small populations remain in Xue, Nanhu, and Cilai Mountains of Taiwan. The distribution of alpine herbs is severely threatened by climate change, which influences genetic variation and population structure. In this study, we investigated the effects of the natural isolation of alpine habitats on the genetic diversity and geographic structure of populations of E. nankotaizanensis using chloroplast (cp) and nuclear DNA (nrDNA) markers. We found lower levels of genetic diversity in E. nankotaizanensis than in other alpine plants and little to no genetic variation within populations, which could be mainly attributed to the small population size and genetic drift. Only one nrDNA haplotype was present in each population. The lack of monophyly of the four populations in cpDNA probably resulted from lineage sorting or occasional long-distance seed dispersal. Phylogeographic analysis suggested that Nanhu Mountain was probably a refugium over the glacial maxima, agreeing with the potential refugia in central Taiwan. The STRUCTURE and AMOVA analyses revealed significant genetic differentiation in nrDNA among the mountains, which resulted from geographical isolation among these mountains. Estimates of the effective population size (Ne) and demography reflected lower Ne values and a recent population decline, probably implying a greater extinction risk for E. nankotaizanensis. We observed genetic depletion and considerable genetic differentiation among mountain populations, which should be considered in future conservation efforts for this species. In addition, this study provides important insights into the long-term potential of alpine herbs in Taiwan, which are useful for a better prediction of their responses to future climate change.

Evolutionary history of Hemerocallis in Japan inferred from chloroplast and nuclear phylogenies and levels of interspecific gene flow

The perennial herb genus Hemerocallis (Asphodelaceae) shows four flowering types: diurnal half-day, diurnal one-day, nocturnal half-day, and nocturnal one-day flowering. These flowering types are corresponding to their main pollinators, and probably act as a primary mechanism of reproductive isolation. To examine how the four flowering types diverged, we reconstructed the phylogeny of the Japanese species of Hemerocallis using 1615 loci of nuclear genome-wide SNPs and 2078 bp sequences of four cpDNA regions. We also examined interspecific gene flows among taxa by an Isolation-with-Migration model and a population structure analysis. Our study revealed an inconsistency between chloroplast and nuclear genome phylogenies, which may have resulted from chloroplast capture. Each of the following five clusters is monophyletic and clearly separated on the nuclear genome-wide phylogenetic tree: (I) two nocturnal flowering species with lemon-yellow flowers, H. citrina (half-day flowering) and H. lilioasphodelus (one-day flowering); (II) a diurnal one-day flowering species with yellow-orange flowers, H. middendorffii; (III) a variety of a diurnal half-day flowering species with reddish orange flowers, H. fulva var. disticha; (IV) another variety of a diurnal half-day flowering species with reddish orange flowers, H. fulva var. aurantiaca, and a diurnal one-day flowering species with yellow-orange flowers, H. major; (V) a diurnal half-day flowering species with yellow-orange flowers, H. hakuunensis. The five clusters are consistent with traditional phenotype-based taxonomy (cluster I, cluster II, and clusters III-V correspond to Hemerocallis sect. Hemerocallis, Capitatae, and Fulvae, respectively). These findings could indicate that three flowering types (nocturnal flowering, diurnal one-day flowering, and diurnal half-day flowering) diverged in early evolutionary stages of Hemerocallis and subsequently a change from diurnal half-day flowering to diurnal one-day flowering occurred in a lineage of H. major. While genetic differentiation among the five clusters was well maintained, significant gene flow was detected between most pairs of taxa, suggesting that repeated hybridization played a role in the evolution of those taxa.

Subspecies divergence and pronounced phylogenetic incongruence in the East-Asia-endemic shrub Magnolia sieboldii

BACKGROUND AND AIMS

The biogeographic patterns of the East-Asia-endemic shrub Magnolia sieboldii, in which the range of the subsp. sieboldii is interposed with the disjunct distribution of subsp. japonica, implies a complex evolutionary history, involving rapid speciation and hybridization. Here, we aim to reveal the evolutionary and phylogeographic histories of the species with a particular focus on the time of subspecies divergence, the hypothesis of secondary hybridization and the Pleistocene survival of each subspecies, using a combination of genetic analyses and ecological niche modelling.

METHODS

Genetic variation, genetic structures and phylogenetic relationships were elucidated based on nuclear low-copy genes, chloroplast DNA, and nuclear simple sequence repeats (SSRs). A scenario selection analysis and divergence time estimation were performed using coalescent simulation in DIYABC and *BEAST. Ecological niche modelling and a test of niche differentiation were performed using Maxent and ENMTools.

KEY RESULTS

All marker types showed deep, but pronouncedly incongruent, west-east genetic divergences, with the subspecies being delineated only by the nuclear low-copy genes. Phylogenetic tree topologies suggested that ancient hybridization and introgression were likely to have occurred; however, this scenario did not receive significant support in the DIYABC analysis. The subspecies differentiated their niches, but both showed a dependence on high humidity and were predicted to have persisted during the last glacial cycle by maintaining a stable latitudinal distribution via migration to lower altitudes.

CONCLUSIONS

We found a deep genetic divergence and a pronounced phylogenetic incongruence among the two subspecies of M. sieboldii, which may have been driven by major paleogeographic and paleoclimatic events that have occurred since the Neogene in East Asia, including global cooling, climate oscillations and the formation of land bridges. Both subspecies were, however, considered to persist in situ in stable climatic conditions during the late Pleistocene.

Turnip Time Travels: Age Estimates in Brassicaceae

Results of research in life sciences acquire a deeper meaning if they can also be discussed in temporal contexts of evolution. Despite the importance of the mustard family (Brassicaceae) as a prominent angiosperm model family, a robust, generally accepted hypothesis for a family-wide temporal framework does not yet exist. The main cause for this situation is a poor fossil record of the family. We suggest that the few known fossils require a critical re-evaluation of phylogenetic and temporal assignments as a prerequisite for appropriate molecular dating analyses within the family. In addition, (palaeo)biogeographical calibrations, not explored so far in the family, should be integrated in a synthesis of various dating approaches, with each contributing their specific possibilities and limitations.

Phylogeography of a widespread species: pre-glacial vicariance, refugia, occasional blocking straits and long-distance migrations

Phylogeographic studies give us the opportunity to reconstruct the historical migrations of species and link them with climatic and geographic variation. They are, therefore, a key tool to understanding the relationships among biology, geology and history. One of the most interesting biogeographical areas of the world is the Mediterranean region. However, in this area, the description of concordant phylogeographic patterns is quite scarce, which limits the understanding of evolutionary patterns related to climate. Species with one-dimensional distribution ranges, such as the strawberry tree (Arbutus unedo), are particularly useful to unravel these patterns. Here, we describe its phylogeographic structure and check for concordance with patterns seen in other Mediterranean plants: longitudinal/latitudinal clines of diversity, evidence for glacial refugia and the role of sea straits in dispersal. We also identify the most likely source for the disjunct Irish population. With this aim, we sequenced four chloroplast non-coding fragments of A. unedo from 23 populations covering its whole distribution. We determined the genetic diversity, population structure, haplotype genealogy and time to the most recent common ancestor. The genealogy revealed two clades that separated during the last 700 ky but before the last glacial maximum. One clade occupies Atlantic Iberia and North Africa, while the other occurs in the Western Mediterranean. The Eastern Mediterranean is inhabited by newer haplotypes derived from both clades, while the Irish population is closely related to Iberian demes. The straits of Sicily and Gibraltar partially restricted the gene flow. We concluded that a vicariance event during the Late Quaternary in the western end of the species' range followed by eastward migration seems a likely explanation for the observed phylogeographic pattern. The role of straits indicates an occasional communication between Europe and North Africa, suggesting that the latter was a novel refugia. The East-West genetic split in Iberia is consistent with the refugia-within-refugia model. Finally, the strawberry tree possibly reached Ireland from Iberia instead of throughout the maritime fringe of France as previously thought.

Phylogeography of a widespread species: pre-glacial vicariance, refugia, occasional blocking straits and long-distance migrations

Phylogeographic studies give us the opportunity to reconstruct the historical migrations of species and link them with climatic and geographic variation. They are, therefore, a key tool to understanding the relationships among biology, geology and history. One of the most interesting biogeographical areas of the world is the Mediterranean region. However, in this area, the description of concordant phylogeographic patterns is quite scarce, which limits the understanding of evolutionary patterns related to climate. Species with one-dimensional distribution ranges, such as the strawberry tree (Arbutus unedo), are particularly useful to unravel these patterns. Here, we describe its phylogeographic structure and check for concordance with patterns seen in other Mediterranean plants: longitudinal/latitudinal clines of diversity, evidence for glacial refugia and the role of sea straits in dispersal. We also identify the most likely source for the disjunct Irish population. With this aim, we sequenced four chloroplast non-coding fragments of A. unedo from 23 populations covering its whole distribution. We determined the genetic diversity, population structure, haplotype genealogy and time to the most recent common ancestor. The genealogy revealed two clades that separated during the last 700 ky but before the last glacial maximum. One clade occupies Atlantic Iberia and North Africa, while the other occurs in the Western Mediterranean. The Eastern Mediterranean is inhabited by newer haplotypes derived from both clades, while the Irish population is closely related to Iberian demes. The straits of Sicily and Gibraltar partially restricted the gene flow. We concluded that a vicariance event during the Late Quaternary in the western end of the species' range followed by eastward migration seems a likely explanation for the observed phylogeographic pattern. The role of straits indicates an occasional communication between Europe and North Africa, suggesting that the latter was a novel refugia. The East-West genetic split in Iberia is consistent with the refugia-within-refugia model. Finally, the strawberry tree possibly reached Ireland from Iberia instead of throughout the maritime fringe of France as previously thought.

Phylogeny of Allium L. subgenus Anguinum (G. Don. ex W.D.J. Koch) N. Friesen (Amaryllidaceae)

The disjunct distribution of the subgenus Anguinum of the genus Allium makes it a good candidate to test models of Northern hemisphere biogeography. Here we conduct phylogenetic analysis with the nuclear marker ITS and three different chloroplast markers (rps16 intron, rbcL-atpB spacer, rpl32-trnL spacer). Divergence time estimations (Beast) relying on published ITS substitution rates and ancestral range reconstructions were calculated to elucidate the biogeographical history of the subgenus. Additionally we compiled distribution maps for all species with data taken from the literature, herbariums and data from field observations. The main radiation of the subgenus took place in the last one million years and is still going on. They have their origins in the mid Miocene in East Asia and were highly influenced by the climate fluctuations in the Pliocene/Pleistocene period. Conflicting tree topologies between nuclear and cpDNA markers of Allium tricoccum Solander indicate that the species is of hybridogenous origin. Cloning the ITS sequence revealed the parental copies and confirmed our conclusion. One originated from the Eurasian and the other from the East Asian clade. We were able to show that it reached North America most likely via the Beringia around 2.5mya (95% HPD of 0.35-5.26mya). Our data suggest that Allium victorialis L. is only distributed in mountain pastures in Europe as it forms a well-supported clade in the ITS tree. In the analysis of the molecular markers we found two distinct types of Allium ochotense Prokh. and we suggest splitting the species based on Prokhanov's (1930) proposal. Taxonomical remarks and an identification key to all species of the subgenus Anguinum is provided.

Phylogeography of Thlaspi arvense (Brassicaceae) in China Inferred from Chloroplast and Nuclear DNA Sequences and Ecological Niche Modeling

Thlaspi arvense is a well-known annual farmland weed with worldwide distribution, which can be found from sea level to above 4000 m high on the Qinghai-Tibetan Plateau (QTP). In this paper, a phylogeographic history of T. arvense including 19 populations from China was inferred by using three chloroplast (cp) DNA segments (trnL-trnF, rpl32-trnL and rps16) and one nuclear (n) DNA segment (Fe-regulated transporter-like protein, ZIP). A total of 11 chloroplast haplotypes and six nuclear alleles were identified, and haplotypes unique to the QTP were recognized (C4, C5, C7 and N4). On the basis of molecular dating, haplotypes C4, C5 and C7 have separated from others around 1.58 Ma for cpDNA, which corresponds to the QTP uplift. In addition, this article suggests that the T. arvense populations in China are a mixture of diverged subpopulations as inferred by hT/vT test (hT ≤ vT, cpDNA) and positive Tajima's D values (1.87, 0.05 < p < 0.10 for cpDNA and 3.37, p < 0.01 for nDNA). Multimodality mismatch distribution curves and a relatively large shared area of suitable environmental conditions between the Last Glacial Maximum (LGM) as well as the present time recognized by MaxEnt software reject the sudden expansion population model.

Efficient Detection of Novel Nuclear Markers for Brassicaceae by Transcriptome Sequencing

The lack of DNA sequence information for most non-model organisms impairs the design of primers that are universally applicable for the study of molecular polymorphisms in nuclear markers. Next-generation sequencing (NGS) techniques nowadays provide a powerful approach to overcome this limitation. We present a flexible and inexpensive method to identify large numbers of nuclear primer pairs that amplify in most Brassicaceae species. We first obtained and mapped NGS transcriptome sequencing reads from two of the distantly related Brassicaceae species, Cardamine hirsuta and Arabis alpina, onto the Arabidopsis thaliana reference genome, and then identified short conserved sequence motifs among the three species bioinformatically. From these, primer pairs to amplify coding regions (nuclear protein coding loci, NPCL) and exon-primed intron-crossing sequences (EPIC) were developed. We identified 2,334 universally applicable primer pairs, targeting 1,164 genes, which provide a large pool of markers as readily usable genomic resource that will help addressing novel questions in the Brassicaceae family. Testing a subset of the newly designed nuclear primer pairs revealed that a great majority yielded a single amplicon in all of the 30 investigated Brassicaceae taxa. Sequence analysis and phylogenetic reconstruction with a subset of these markers on different levels of phylogenetic divergence in the mustard family were compared with previous studies. The results corroborate the usefulness of the newly developed primer pairs, e.g., for phylogenetic analyses or population genetic studies. Thus, our method provides a cost-effective approach for designing nuclear loci across a broad range of taxa and is compatible with current NGS technologies.

The genomes of many yam species contain transcriptionally active endogenous geminiviral sequences that may be functionally expressed

Endogenous viral sequences are essentially 'fossil records' that can sometimes reveal the genomic features of long extinct virus species. Although numerous known instances exist of single-stranded DNA (ssDNA) genomes becoming stably integrated within the genomes of bacteria and animals, there remain very few examples of such integration events in plants. The best studied of these events are those which yielded the geminivirus-related DNA elements found within the nuclear genomes of various Nicotiana species. Although other ssDNA virus-like sequences are included within the draft genomes of various plant species, it is not entirely certain that these are not contaminants. The Nicotiana geminivirus-related DNA elements therefore remain the only definitively proven instances of endogenous plant ssDNA virus sequences. Here, we characterize two new classes of endogenous plant virus sequence that are also apparently derived from ancient geminiviruses in the genus Begomovirus. These two endogenous geminivirus-like elements (EGV1 and EGV2) are present in the Dioscorea spp. of the Enantiophyllum clade. We used fluorescence in situ hybridization to confirm that the EGV1 sequences are integrated in the D. alata genome and showed that one or two ancestral EGV sequences likely became integrated more than 1.4 million years ago during or before the diversification of the Asian and African Enantiophyllum Dioscorea spp. Unexpectedly, we found evidence of natural selection actively favouring the maintenance of EGV-expressed replication-associated protein (Rep) amino acid sequences, which clearly indicates that functional EGV Rep proteins were probably expressed for prolonged periods following endogenization. Further, the detection in D. alata of EGV gene transcripts, small 21-24 nt RNAs that are apparently derived from these transcripts, and expressed Rep proteins, provides evidence that some EGV genes are possibly still functionally expressed in at least some of the Enantiophyllum clade species.

One or three species in Megadenia (Brassicaceae): insight from molecular studies

Megadenia Maxim. is a small genus of the Brassicaceae endemic to East Asia with three disjunct areas of distribution: the eastern edge of the Qinghai-Tibetan Plateau, the Eastern Sayan Mountains in southern Siberia, and Chandalaz Ridge in the southern Sikhote-Alin Mountains. Although distinct species (M. pygmaea Maxim., M. bardunovii Popov, and M. speluncarum Vorob., Vorosch. and Gorovoj) have been described from each area, they have lately been reduced to synonymy with M. pygmaea due to high morphological similarity. Here, we present the first molecular study of Megadenia. Using the sequences of 11 noncoding regions from the cytoplasmic (chloroplast and mitochondrial) and nuclear genomes, we assessed divergence within the genus and explored the relationships between Megadenia and Biscutella L. Although M. bardunovii, M. speluncarum, and M. pygmaea were found to be indiscernible with regard to the nuclear and mitochondrial markers studied, our data on the plastid genome revealed their distinctness and a clear subdivision of the genus into three lineages matching the three described species. All of the phylogenetic analyses of the chloroplast DNA sequences provide strong support for the inclusion of Megadenia and Biscutella in the tribe Biscutelleae. A dating analysis shows that the genus Megadenia is of Miocene origin and diversification within the genus, which has led to the three extant lineages, most likely occurred during the Early-Middle Pleistocene, in agreement with the vicariance pattern. Given the present-day distribution, differences in habitat preferences and in some anatomical traits, and lack of a direct genealogical relationship, M. pygmaea, M. bardunovii, and M. speluncarum should be treated as distinct species or at least subspecies.

Plant mitochondrial genome evolution can be explained by DNA repair mechanisms

Plant mitochondrial genomes are notorious for their large and variable size, nonconserved open reading frames of unknown function, and high rates of rearrangement. Paradoxically, the mutation rates are very low. However, mutation rates can only be measured in sequences that can be aligned—a very small part of plant mitochondrial genomes. Comparison of the complete mitochondrial genome sequences of two ecotypes of Arabidopsis thaliana allows the alignment of noncoding as well as coding DNA and estimation of the mutation rates in both. A recent chimeric duplication is also analyzed. A hypothesis is proposed that the mechanisms of plant mitochondrial DNA repair account for these features and includes different mechanisms in transcribed and nontranscribed regions. Within genes, a bias toward gene conversion would keep measured mutation rates low, whereas in noncoding regions, break-induced replication (BIR) explains the expansion and rearrangements. Both processes are types of double-strand break repair, but enhanced second-strand capture in transcribed regions versus BIR in nontranscribed regions can explain the two seemingly contradictory features of plant mitochondrial genome evolution—the low mutation rates in genes and the striking expansions of noncoding sequences.

Plant sex chromosome evolution

It is now well established that plants have an important place in studies of sex chromosome evolution because of the repeated independent evolution of separate sexes and sex chromosomes. There has been considerable recent progress in studying plant sex chromosomes. In this review, I focus on how these recent studies have helped clarify or answer several important questions about sex chromosome evolution, and I shall also try to clarify some common misconceptions. I also outline future work that will be needed to make further progress, including testing some important ideas by genetic, molecular, and developmental approaches. Systems with different ages can clearly help show the time course of events during changes from an ancestral co-sexual state (hermaphroditism or monoecy), and I will also explain how different questions can be studied in lineages whose dioecy or sex chromosomes evolved at different times in the past.