Search

Given the 2,400-fold range of genome sizes (0.06–148.9 Gbp (gigabase pair)) of seed plants (angiosperms and gymnosperms) with a broadly similar gene content (amounting to approximately 0.03 Gbp), the repeat-sequence content of the genome might be expected to increase with genome size, resulting in the largest genomes consisting almost entirely of repetitive sequences. Here we test this prediction, using the same bioinformatic approach for 101 species to ensure consistency in what constitutes a repeat. We reveal a fundamental change in repeat turnover in genomes above around 10 Gbp, such that species with the largest genomes are only about 55% repetitive. Given that genome size influences many plant traits, habits and life strategies, this fundamental shift in repeat dynamics is likely to affect the evolutionary trajectory of species lineages., We thank Natural Environment Research Council (NE/G020256/1), the Czech Academy of Sciences (RVO:60077344) and Ramón y Cajal Fellowship (RYC-2017-2274) funded by the Ministerio de Ciencia y Tecnología (Gobierno de España) for support. We also thank Natural Environment Research Council for funding a studentship to S.D. and the China Scholarship Council for funding W.W., Main Methods Data availability Code availability References Acknowledgements Author information Ethics declarations Additional information Extended data Supplementary information Rights and permissions About this article Further reading

Societal Impact Statement For millennia, humans have used plants and fungi, as foods, fuels, fibers, and medicines; and have developed techniques for improving their usefulness to our species, mostly through selection of desirable traits. With human populations forecast to rise, the availability of arable land likely to fall amid climate change and increasing urbanization, and modern communications technologies accelerating the dispersal of pathogens, further improvement is urgently needed. However, ensuring long‐term resilience involves conservation of existing genetic diversity in addition to selection. New technologies, particularly those based on molecular biology, are increasingly driving conservation and improvement strategies. Summary Humans use plants and fungi for a wide range of purposes and, over millennia, have improved wild species by selecting for and combining genetic variation. Improvements in DNA sequencing technologies have enhanced our capacity to identify and manipulate genetic diversity, increasing the range of variation that can be utilized, and accelerating the breeding cycle to reduce the time taken to develop and put new varieties to use. Most recently, the CRISPR/Cas9 gene editing technology has greatly increased our capacity to directly introduce novel genetic variants without unwanted associated material. Moreover, increased knowledge of metabolic pathways resulting from genomic analysis can be used to design new varieties with desired properties with increased precision. Selecting for, or engineering, desirable variants has increased the usefulness of plants and fungi to humans, but at the cost of reducing their genetic diversity, decreasing their resilience and reducing the stock of variation available for future use. Conservation of genetic biodiversity is thus an essential counterpart of crop improvement and is essential to ensure that crop species retain resilience to emerging threats. Conservation efforts are focused on orphan crops, wild relatives of crop species, and landraces; in and exsitu efforts are complementary. Informatic approaches can inform use of these materials in breeding programmes even in the absence of genomic information. The application of some of these approaches may be restricted by ethical, legal, or organizational obstacles. If these can be overcome, there is great potential to unlock previously untapped reservoirs of biodiversity for human benefit.

Recent studies show that molecular convergence plays an unexpectedly common role in the evolution of convergent phenotypes. We exploited this phenomenon to find candidate loci underlying resistance to the emerald ash borer (EAB, Agrilus planipennis), the United States’ most costly invasive forest insect to date, within the pan-genome of ash trees (the genus Fraxinus). We show that EAB-resistant taxa occur within three independent phylogenetic lineages. In genomes from these resistant lineages, we detect 53 genes with evidence of convergent amino acid evolution. Gene-tree reconstruction indicates that, for 48 of these candidates, the convergent amino acids are more likely to have arisen via independent evolution than by another process such as hybridization or incomplete lineage sorting. Seven of the candidate genes have putative roles connected to the phenylpropanoid biosynthesis pathway and 17 relate to herbivore recognition, defence signalling or programmed cell death. Evidence for loss-of-function mutations among these candidates is more frequent in susceptible species than in resistant ones. Our results on evolutionary relationships, variability in resistance, and candidate genes for defence response within the ash genus could inform breeding for EAB resistance, facilitating ecological restoration in areas invaded by this beetle. By assembling the genomes of 22 Fraxinus species and conducting comparisons including further species, the authors identify candidate loci for emerald ash borer resistance that have evolved convergently.

Green ash (Fraxinus pennsylvanica) is the most widely distributed ash tree in North America. Once common, it has experienced high mortality from the non-native invasive emerald ash borer (EAB; Agrilus planipennis). A small percentage of native green ash trees that remain healthy in long-infested areas, termed "lingering ash," display partial resistance to the insect, indicating that breeding and propagating populations with higher resistance to EAB may be possible. To assist in ash breeding, ecology and evolution studies, we report the first chromosome-level assembly from the genus Fraxinus for F. pennsylvanica with over 99% of bases anchored to 23 haploid chromosomes, spanning 757 Mb in total, composed of 49.43% repetitive DNA, and containing 35,470 high-confidence gene models assigned to 22,976 Asterid orthogroups. We also present results of range-wide genetic variation studies, the identification of candidate genes for important traits including potential EAB-resistance genes, and an investigation of comparative genome organization among Asterids based on this reference genome platform. Residual duplicated regions within the genome probably resulting from a recent whole genome duplication event in Oleaceae were visualized in relation to wild olive (Olea europaea var. sylvestris). We used our F. pennsylvanica chromosome assembly to construct reference-guided assemblies of 27 previously sequenced Fraxinus taxa, including F. excelsior. Thus, we present a significant step forward in genomic resources for research and protection of Fraxinus species.

Societal Impact Statement The current ash dieback epidemic in Europe caused by Hymenoscyphus fraxineus poses a key question to policy makers: whether or not to commit time and resources to the initiation of a breeding programme for the development of more resistant ash, as a long‐term policy of adaptation to the epidemic. Here we review current evidence on the potential viability of such a programme, from a biological perspective. We conclude that a breeding programme for ash aimed at resistance to current strains of H. fraxineus in the British Isles is biologically feasible. Summary To evaluate the viability and feasibility of a future breeding programme to produce trees resistant to an emerging pest or pathogen, it is helpful to ask the following questions: How much variation in resistance exists in tree populations? To what extent is this resistance heritable? How many genetic loci are involved? What level of resistance is found in other species of the same genus? Here, we survey current knowledge of these issues in relation to the degree of resistance of European ash (Fraxinus excelsior) to H. fraxineus, the fungus causing ash dieback (ADB). Several studies have found a low frequency of heritable resistance in F. excelsior populations, which seems to be determined by many genetic loci. This suggests that a breeding programme is viable and that natural selection may also increase the mean resistance of populations over time. More research is needed on the genetic basis of resistance to ADB to understand how quickly natural selection can operate in woodlands and what acceleration may be possible in breeding programmes, including via use of genetic markers. Hybrid breeding programmes may also be a possibility, as some ash species appear to be more resistant to ADB than is F. excelsior, but more research is needed on this issue. We do not yet know if it will be possible to breed F. excelsior to have resistance to both ADB and the emerging threat of emerald ash borer. We recommend short‐term mitigation measures for the ADB epidemic and future research directions.

Numerous plant genera have a history including frequent hybridisation and polyploidisation, which often means that their phylogenies are not yet fully resolved. The genus Betula, which contains many ecologically important allopolyploid tree species, is a case in point. We generated genome-wide sequence data for 27 diploid and 31 polyploid Betula species or subspecies using restriction site associated DNA (RAD) sequences assembled into contigs with a mean length of 675 bp. We reconstructed the evolutionary relationships among diploid Betula species using both supermatrix and species tree methods. We identified progenitors of the polyploids according to the relative rates at which their reads mapped to contigs from different diploid species. We sorted the polyploid reads into different putative sub-genomes and used the extracted contigs, along with the diploid sequences, to build new phylogenies that included the polyploid sub-genomes. This approach yielded a highly evidenced phylogenetic hypothesis for the genus Betula, including the complex reticulate origins of the majority of its polyploid taxa. The genus was split into two well supported clades, which differ in their seed-wing morphology. We propose a new taxonomy for Betula, splitting it into two subgenera. We have resolved the parentage of many widespread and economically important polyploid tree species, opening the way for their population genomic study.

Genome-wide discovery of candidate genes for functional traits within a species typically involves the sequencing of large samples of phenotyped individuals1, or linkage analysis through multiple generations2. When a trait occurs repeatedly among phylogenetically independent lineages within a genus, a more efficient approach may be to identify genes via detection of amino acid residues shared by species possessing that trait3,4. Here, by taking this approach, we identify candidate loci in the genus Fraxinus (ash trees) for resistance to the emerald ash borer beetle (EAB; Agrilus planipennis), a pest species that appears innocuous to otherwise healthy ash in its native East Asian range5 but is highly destructive in North America6 and poses a threat to ash trees in Europe7. Assembling whole genome sequences for 24 diploid species and subspecies of ash, and estimating resistance to EAB for 26 taxa from egg bioassays, we find 53 genes containing amino acid variants shared between two or more independent Fraxinus lineages with EAB-resistant species, that are unlikely to be due to chance or undetected paralogy. Of these, seven genes have putative roles relating to the phenylpropanoid biosynthesis pathway and 17 are potentially connected to herbivore recognition, defence signalling or programmed cell death. We also find that possible loss-of-function mutations among our 53 candidate genes are more frequent in susceptible species, than in resistant ones. Patterns of polymorphism for the EAB-associated amino acid variants in ash trees representing different European populations suggest that selection may be able to enhance their resistance to EAB.

SummaryPopulations of European ash trees (Fraxinus excelsior) are being devastated by the invasive alien fungusHymenoscyphus fraxineus, which causes ash dieback (ADB). We sequenced whole genomic DNA from 1250 ash trees in 31 DNA pools, each pool containing trees with the same ADB damage status in a screening trial and from the same seed-source zone. A genome-wide association study (GWAS) identified 3,149 single nucleotide polymorphisms (SNPs) associated with low versus high ADB damage. Sixty-one of the 203 most significant SNPs were in, or close to, genes with putative homologs already known to be involved in pathogen responses in other plant species. We also used the pooled sequence data to train a genomic prediction (GP) model, cross-validated using individual whole genome sequence data generated for 75 healthy and 75 damaged trees from a single seed source. Using the top 30% of our genomic estimated breeding values from 200 SNPs, we could predict tree health with over 90% accuracy. We infer that ash dieback resistance inF. excelsioris a polygenic trait that should respond well to both natural selection and breeding, which could be accelerated using GP.

Populations of European ash trees (Fraxinus excelsior) are being devastated by the invasive alien fungus Hymenoscyphus fraxineus, which causes ash dieback. We sequenced whole genomic DNA from 1,250 ash trees in 31 DNA pools, each pool containing trees with the same ash dieback damage status in a screening trial and from the same seed-source zone. A genome-wide association study identified 3,149 single nucleotide polymorphisms (SNPs) associated with low versus high ash dieback damage. Sixty-one of the 192 most significant SNPs were in, or close to, genes with putative homologues already known to be involved in pathogen responses in other plant species. We also used the pooled sequence data to train a genomic prediction model, cross-validated using individual whole genome sequence data generated for 75 healthy and 75 damaged trees from a single seed source. The model's genomic estimated breeding values (GEBVs) allocated these 150 trees to their observed health statuses with 67% accuracy using 10,000 SNPs. Using the top 20% of GEBVs from just 200 SNPs, we could predict observed tree health with over 90% accuracy. We infer that ash dieback resistance in F. excelsior is a polygenic trait that should respond well to both natural selection and breeding, which could be accelerated using genomic prediction.

Numerous plant genera have a history including frequent hybridisation and polyploidisation (allopolyploidisation), which means that their phylogeny is a network of reticulate evolution that cannot be accurately depicted as a bifurcating tree with a single tip per species. The genus Betula, which contains many ecologically important tree species, is a case in point. We generated genome-wide sequence reads for 27 diploid and 36 polyploid Betula species or subspecies using restriction site associated DNA (RAD) sequences. These reads were assembled into contigs with a mean length of 675 bp. We reconstructed the evolutionary relationships among diploid Betula species using both supermatrix (concatenation) and species tree methods. We identified the closest diploid relatives of the polyploids according to the relative rates at which reads from polyploids mapped to contigs from different diploid species within a concatenated reference sequence. By mapping reads from allopolyploids to their different putative diploid relatives we assembled contigs from the putative sub-genomes of allopolyploid taxa. We used these to build new phylogenies that included allopolyploid sub-genomes as separate tips. This approach yielded a highly evidenced phylogenetic hypothesis for the genus Betula, including the complex reticulate origins of the majority of its polyploid taxa. Our phylogeny divides the genus into two well supported clades, which, interestingly, differ in their seed-wing morphology. We therefore propose to split Betula into two subgenera.

Plants exhibit an extraordinary range of genome sizes, varying by > 2000-fold between the smallest and largest recorded values. In the absence of polyploidy, changes in the amount of repetitive DNA (transposable elements and tandem repeats) are primarily responsible for genome size differences between species. However, there is ongoing debate regarding the relative importance of amplification of repetitive DNA versus its deletion in governing genome size. Using data from 454 sequencing, we analysed the most repetitive fraction of some of the largest known genomes for diploid plant species, from members of Fritillaria. We revealed that genomic expansion has not resulted from the recent massive amplification of just a handful of repeat families, as shown in species with smaller genomes. Instead, the bulk of these immense genomes is composed of highly heterogeneous, relatively low-abundance repeat-derived DNA, supporting a scenario where amplified repeats continually accumulate due to infrequent DNA removal. Our results indicate that a lack of deletion and low turnover of repetitive DNA are major contributors to the evolution of extremely large genomes and show that their size cannot simply be accounted for by the activity of a small number of high-abundance repeat families.

Gnetophytes are an enigmatic gymnosperm lineage comprising three genera, Gnetum, Welwitschia and Ephedra, which are morphologically distinct from all other seed plants. Their distinctiveness has triggered much debate as to their origin, evolution and phylogenetic placement among seed plants. To increase our understanding of the evolution of gnetophytes, and their relation to other seed plants, we report here a high-quality draft genome sequence for Gnetum montanum, the first for any gnetophyte. By using a novel genome assembly strategy to deal with high levels of heterozygosity, we assembled >4 Gb of sequence encoding 27,491 protein-coding genes. Comparative analysis of the G. montanum genome with other gymnosperm genomes unveiled some remarkable and distinctive genomic features, such as a diverse assemblage of retrotransposons with evidence for elevated frequencies of elimination rather than accumulation, considerable differences in intron architecture, including both length distribution and proportions of (retro) transposon elements, and distinctive patterns of proliferation of functional protein domains. Furthermore, a few gene families showed Gnetum-specific copy number expansions (for example, cellulose synthase) or contractions (for example, Late Embryogenesis Abundant protein), which could be connected with Gnetum’s distinctive morphological innovations associated with their adaptation to warm, mesic environments. Overall, the G. montanum genome enables a better resolution of ancestral genomic features within seed plants, and the identification of genomic characters that distinguish Gnetum from other gymnosperms.

Endogenous pararetroviral sequences are the most commonly found virus sequences integrated into angiosperm genomes. We describe an endogenous pararetrovirus (EPRV) repeat in Fritillaria imperialis, a species that is under study as a result of its exceptionally large genome (1C = 42 096 Mbp, approximately 240 times bigger than Arabidopsis thaliana). The repeat (FriEPRV) was identified from Illumina reads using the RepeatExplorer pipeline, and exists in a complex genomic organization at the centromere of most, or all, chromosomes. The repeat was reconstructed into three consensus sequences that formed three interconnected loops, one of which carries sequence motifs expected of an EPRV (including the gag and pol domains). FriEPRV shows sequence similarity to members of the Caulimoviridae pararetrovirus family, with phylogenetic analysis indicating a close relationship to Petuvirus. It is possible that no complete EPRV sequence exists, although our data suggest an abundance that exceeds the genome size of Arabidopsis. Analysis of single nucleotide polymorphisms revealed elevated levels of C→T and G→A transitions, consistent with deamination of methylated cytosine. Bisulphite sequencing revealed high levels of methylation at CG and CHG motifs (up to 100%), and 15-20% methylation, on average, at CHH motifs. FriEPRV's centromeric location may suggest targeted insertion, perhaps associated with meiotic drive. We observed an abundance of 24 nt small RNAs that specifically target FriEPRV, potentially providing a signature of RNA-dependent DNA methylation. Such signatures of epigenetic regulation suggest that the huge genome of F. imperialis has not arisen as a consequence of a catastrophic breakdown in the regulation of repeat amplification.

A large proportion of genomic information, particularly repetitive elements, is usually ignored when researchers are using next-generation sequencing. Here we demonstrate the usefulness of this repetitive fraction in phylogenetic analyses, utilizing comparative graph-based clustering of next-generation sequence reads, which results in abundance estimates of different classes of genomic repeats. Phylogenetic trees are then inferred based on the genome-wide abundance of different repeat types treated as continuously varying characters; such repeats are scattered across chromosomes and in angiosperms can constitute a majority of nuclear genomic DNA. In six diverse examples, five angiosperms and one insect, this method provides generally well-supported relationships at interspecific and intergeneric levels that agree with results from more standard phylogenetic analyses of commonly used markers. We propose that this methodology may prove especially useful in groups where there is little genetic differentiation in standard phylogenetic markers. At the same time as providing data for phylogenetic inference, this method additionally yields a wealth of data for comparative studies of genome evolution.(Repetitive DNA; continuous characters; genomics; next-generation sequencing; phylogenetics; molecular systematics.)

• Since the occurrence of giant genomes in angiosperms is restricted to just a few lineages, identifying where shifts towards genome obesity have occurred is essential for understanding the evolutionary mechanisms triggering this process. • Genome sizes were assessed using flow cytometry in 79 species and new chromosome numbers were obtained. Phylogenetically based statistical methods were applied to infer ancestral character reconstructions of chromosome numbers and nuclear DNA contents. • Melanthiaceae are the most diverse family in terms of genome size, with C-values ranging more than 230-fold. Our data confirmed that giant genomes are restricted to tribe Parideae, with most extant species in the family characterized by small genomes. Ancestral genome size reconstruction revealed that the most recent common ancestor (MRCA) for the family had a relatively small genome (1C = 5.37 pg). Chromosome losses and polyploidy are recovered as the main evolutionary mechanisms generating chromosome number change. • Genome evolution in Melanthiaceae has been characterized by a trend towards genome size reduction, with just one episode of dramatic DNA accumulation in Parideae. Such extreme contrasting profiles of genome size evolution illustrate the key role of transposable elements and chromosome rearrangements in driving the evolution of plant genomes.

Recent advances have highlighted the ubiquity of whole-genome duplication (polyploidy) in angiosperms, although subsequent genome size change and diploidization (returning to a diploid-like condition) are poorly understood. An excellent system to assess these processes is provided by Nicotiana section Repandae, which arose via allopolyploidy (approximately 5 million years ago) involving relatives of Nicotiana sylvestris and Nicotiana obtusifolia. Subsequent speciation in Repandae has resulted in allotetraploids with divergent genome sizes, including Nicotiana repanda and Nicotiana nudicaulis studied here, which have an estimated 23.6% genome expansion and 19.2% genome contraction from the early polyploid, respectively. Graph-based clustering of next-generation sequence data enabled assessment of the global genome composition of these allotetraploids and their diploid progenitors. Unexpectedly, in both allotetraploids, over 85% of sequence clusters (repetitive DNA families) had a lower abundance than predicted from their diploid relatives; a trend seen particularly in low-copy repeats. The loss of high-copy sequences predominantly accounts for the genome downsizing in N. nudicaulis. In contrast, N. repanda shows expansion of clusters already inherited in high copy number (mostly chromovirus-like Ty3/Gypsy retroelements and some low-complexity sequences), leading to much of the genome upsizing predicted. We suggest that the differential dynamics of low- and high-copy sequences reveal two genomic processes that occur subsequent to allopolyploidy. The loss of low-copy sequences, common to both allopolyploids, may reflect genome diploidization, a process that also involves loss of duplicate copies of genes and upstream regulators. In contrast, genome size divergence between allopolyploids is manifested through differential accumulation and/or deletion of high-copy-number sequences.

Genome sequencing has been restricted to species with a small genome size. With the advent of second- and third-generation sequencing technologies, the potential to sequence genomes of all sizes is becoming a reality. As the field of whole genome sequencing has developed, there has been a growing appreciation of the need to better represent the major lineages of the plant tree of life, rather than just those that contain economically important taxa. We argue that as well as accounting for phylogenetic diversity when selecting species to analyse, in order to gain a comprehensive understanding of genome evolution, large-scale sequencing projects also need to reflect the diversity of genome sizes in plants. In this article we briefly outline evidence from the literature to support this view.

Nicotiana (Solanaceae) provides an ideal system for understanding polyploidization, a pervasive and powerful evolutionary force in plants, as this genus contains several groups of allotetraploids that formed at different times from different diploid progenitors. However, the parental lineages of the largest group of allotetraploids, Nicotiana section Suaveolentes, have been problematic to identify. Using data from four regions of three low-copy nuclear genes, nuclear ribosomal DNA, and regions of the plastid genome, we have reconstructed the evolutionary origin of sect. Suaveolentes and identified the most likely diploid progenitors by using a combination of gene trees and network approaches to uncover the most strongly supported evidence of species relationships. Our analyses best support a scenario where a member of the sect. Sylvestres lineage acted as the paternal progenitor and a member of either sect. Petunioides or sect. Noctiflorae that also contained introgressed DNA from the other, or a hypothetical hybrid species between these two sections, was the maternal progenitor. Nicotiana exemplifies many of the factors that can complicate the reconstruction of polyploid evolutionary history and highlights how reticulate evolution at the diploid level can add even greater complexity to allopolyploid genomes.

Se evaluaron seis regiones de ADN como posibles códigos de barras de ADN para Hongos, el segundo reino más grande de vida eucariota, por un consorcio multinacional y de laboratorios múltiples. La región de la mitocondria se utilizó la subunidad 1 de la citocromo c oxidasa como código de barras animal excluido como un marcador potencial, porque es difícil de amplificar en hongos, a menudo incluye intrones grandes y puede ser insuficientemente variable. Tres subunidades del cistrón de ARN ribosómico nuclear fueron en comparación con las regiones de tres genes de codificación proteica representativos (la subunidad más grande de la ARN polimerasa II, la segunda más grande subunidad de la ARN polimerasa II y mantenimiento de minicromosomas proteína). Aunque las regiones de genes que codifican proteínas a menudo tenían un porcentaje más alto de identificación correcta en comparación con ribosomal marcadores, baja amplificación por PCR y éxito de secuenciación eliminados ellos como candidatos para un código de barras universal para hongos. Entre la regiones del cistrón ribosómico, el espaciador transcrito interno (ITS) tiene la mayor probabilidad de identificación exitosa para la gama más amplia de hongos, con la brecha de código de barras más claramente definida entre la variación inter e intraespecífica. El nuclear subunidad ribosómica grande, un marcador filogenético popular en ciertos grupos, tuvieron una resolución de especies superior en algunos grupos taxonómicos, como los primeros linajes divergentes y las levaduras ascomicetas, pero por lo demás era ligeramente inferior al ITS. El ribosomal nuclear La subunidad pequeña tiene una resolución pobre a nivel de especie en los hongos. SU será propuesto formalmente para su adopción como el código de barras fúngico principal marcador al Consorcio para el Código de Barras de la Vida, con la posibilidad de que se desarrollen códigos de barras suplementarios para grupos taxonómicos estrechamente circunscritos., Six DNA regions were evaluated as potential DNA barcodes for Fungi, the second largest kingdom of eukaryotic life, by a multinational, multilaboratory consortium. The region of the mitochondrial cytochrome c oxidase subunit 1 used as the animal barcode was excluded as a potential marker, because it is difficult to amplify in fungi, often includes large introns, and can be insufficiently variable. Three subunits from the nuclear ribosomal RNA cistron were compared together with regions of three representative proteincoding genes (largest subunit of RNA polymerase II, second largest subunit of RNA polymerase II, and minichromosome maintenance protein). Although the protein-coding gene regions often had a higher percent of correct identification compared with ribosomal markers, low PCR amplification and sequencing success eliminated them as candidates for a universal fungal barcode. Among the regions of the ribosomal cistron, the internal transcribed spacer (ITS) region has the highest probability of successful identification for the broadest range of fungi, with the most clearly defined barcode gap between inter- and intraspecific variation. The nuclear ribosomal large subunit, a popular phylogenetic marker in certain groups, had superior species resolution in some taxonomic groups, such as the early diverging lineages and the ascomycete yeasts, but was otherwise slightly inferior to the ITS. The nuclear ribosomal small subunit has poor species-level resolution in fungi. ITS will be formally proposed for adoption as the primary fungal barcode marker to the Consortium for the Barcode of Life, with the possibility that supplementary barcodes may be developed for particular narrowly circumscribed taxonomic groups.

Genome size in plants is characterised by its extraordinary range. Although it appears that the majority of plants have small genomes, in several lineages genome size has reached giant proportions. The recent advent of next-generation sequencing (NGS) methods has for the first time made detailed analysis of even the largest of plant genomes a possibility. In this review, we highlight investigations that have utilised NGS for the study of plants with large genomes, as well as describing ongoing work that aims to harness the power of these technologies to gain insights into their evolution. In addition, we emphasise some areas of research where the use of NGS has the potential to generate significant advances in our current understanding of how plant genomes evolve. Finally, we discuss some of the future developments in sequencing technology that may further improve our ability to explore the content and evolutionary dynamics of the very largest genomes.

DNA barcoding has the possibility for application in many fields, ranging from biological inventory and species discov ery to forensics and quality control. Barcoding could have great utility within the aquatic angiosperm family Podostemaceae, as these plants can be extremely difficult to identify on the basis of morphological characters alone, but are important components of tropical rivers and face threats from increased water pollution. Several potential barcodes for land plants have been suggested, including two 3-locus combinations proposed by the Plant Working Group of the Consortium for the Barcode of Life. These loci, matK, rpoB, rpoCl, and trnH-psbA, were tested for their potential utility in African members of the Podostemaceae. Our results show that although the rpoCl region performs best in terms of universality of amplification and sequencing, the superior resolving power of matK makes it the most promising of the regions for use in barcoding. In the light of these results, we discuss the potential to apply barcoding to Podostemaceae of Africa and make suggestions for the directions of future work to exploit the full potential of this new technology in this group of plants.

Reticulate evolution may function both at the species level, through homoploid and polyploid hybridization, and below the species level, through inter and intragenic recombination. These processes represent challenges for the reconstruction of evolutionary relationships between species, because they cannot be represented adequately with bifurcating trees. We use data from low-copy nuclear genes to evaluate long-standing hypotheses of homoploid (interspecific) hybrid speciation in Nicotiana (Solanaceae) and reconstruct a complex series of reticulation events that have been important in the evolutionary history of this genus. Hybrid origins for three diploid species (Nicotiana glauca, N. linearis, and N. spegazzinii) are inferred on the basis of gene tree incongruence, evidence for interallelic recombination between likely parental alleles, and support for incompatible splits in Lento plots. Phylogenetic analysis of recombinant gene sequences illustrates that recombinants may be resolved with one of their progenitor lineages with a high posterior probability under Bayesian inference, and thus there is no indication of the conflict between phylogenetic signals that results from reticulation. Our results illustrate the importance of hybridization in shaping evolution in Nicotiana and also show that intragenic recombination may be relatively common. This finding demonstrates that it is important to investigate the possibility of recombination when aiming to detect hybrids from DNA-sequence data and reconstruct patterns of reticulate evolution between species.

Ford, C. S., Ayres, K. L., Toomey, N., Haider, N., Stahl, J. V., Kelly, L. J., Wikstrom, N., Hollingsworth, P. M., Duff, R. J., Hoot, S. B., Cowan, R. S., Chase, M. W., Wilkinson, M. J. (2009). Selection of candidate coding DNA barcoding regions for use on land plants. Botanical Journal of the Linnean Society, 159, (1), 1-11. Sponsorship: Alfred P. Sloan Foundation Gordon and Betty Moore Foundation IMPF: 00.98 RONO: 00

Chase, M. W., Cowan, R. S., Hollingsworth, P. M., van den Berg, C., Madrinan, S., Petersen, G., Seberg, O., Jorgsensen, T., Cameron, K. M., Carine, M., Pedersen, N., Hedderson, T. A. J., Conrad, F., Salazar, G. A., Richardson, J. E., Hollingsworth, M. L., Barraclough, T. G., Kelly, L., Wilkinson, M. (2007). A proposal for a standardised protocol to barcode all land plants. Taxon, 56, (2), 295-299.

The differential accumulation and elimination of repetitive DNA are key drivers of genome size variation in flowering plants, yet there have been few studies which have analysed how different types of repeats in related species contribute to genome size evolution within a phylogenetic context. This question is addressed here by conducting large-scale comparative analysis of repeats in 23 species from four genera of the monophyletic legume tribe Fabeae, representing a 7.6-fold variation in genome size. Phylogenetic analysis and genome size reconstruction revealed that this diversity arose from genome size expansions and contractions in different lineages during the evolution of Fabeae. Employing a combination of low-pass genome sequencing with novel bioinformatic approaches resulted in identification and quantification of repeats making up 55?83% of the investigated genomes. In turn, this enabled an analysis of how each major repeat type contributed to the genome size variation encountered. Differential accumulation of repetitive DNA was found to account for 85% of the genome size differences between the species, and most (57%) of this variation was found to be driven by a single lineage of Ty3/gypsy LTR-retrotransposons, the Ogre elements. Although the amounts of several other lineages of LTR-retrotransposons and the total amount of satellite DNA were also positively correlated with genome size, their contributions to genome size variation were much smaller (up to 6%). Repeat analysis within a phylogenetic framework also revealed profound differences in the extent of sequence conservation between different repeat types across Fabeae. In addition to these findings, the study has provided a proof of concept for the approach combining recent developments in sequencing and bioinformatics to perform comparative analyses of repetitive DNAs in a large number of non-model species without the need to assemble their genomes.

The RNA-directed DNA methylation (RdDM) pathway can be divided into three phases: 1) small interfering RNA biogenesis, 2) de novo methylation, and 3) chromatin modification. To determine the degree of conservation of this pathway we searched for key genes among land plants. We used OrthoMCL and the OrthoMCL Viridiplantae database to analyze proteomes of species in bryophytes, lycophytes, monilophytes, gymnosperms, and angiosperms. We also analyzed small RNA size categories and, in two gymnosperms, cytosine methylation in ribosomal DNA. Six proteins were restricted to angiosperms, these being NRPD4/NRPE4, RDM1, DMS3 (defective in meristem silencing 3), SHH1 (SAWADEE homeodomain homolog 1), KTF1, and SUVR2, although we failed to find the latter three proteins in Fritillaria persica, a species with a giant genome. Small RNAs of 24 nt in length were abundant only in angiosperms. Phylogenetic analyses of Dicer-like (DCL) proteins showed that DCL2 was restricted to seed plants, although it was absent in Gnetum gnemon and Welwitschia mirabilis. The data suggest that phases (1) and (2) of the RdDM pathway, described for model angiosperms, evolved with angiosperms. The absence of some features of RdDM in F. persica may be associated with its large genome. Phase (3) is probably the most conserved part of the pathway across land plants. DCL2, involved in virus defense and interaction with the canonical RdDM pathway to facilitate methylation of CHH, is absent outside seed plants. Its absence in G. gnemon, and W. mirabilis coupled with distinctive patterns of CHH methylation, suggest a secondary loss of DCL2 following the divergence of Gnetales.

1095-8290McCarthy, Elizabeth WArnold, Sarah E JChittka, LarsLe Comber, Steven CVerity, RobertDodsworth, StevenKnapp, SandraKelly, Laura JChase, Mark WBaldwin, Ian TKovarik, AlesMhiri, CorinneTaylor, LinLeitch, Andrew RJournal ArticleEnglandAnn Bot. 2015 Jun;115(7):1117-31. doi: 10.1093/aob/mcv048. Epub 2015 May 15.; BACKGROUND AND AIMS: Speciation in angiosperms can be accompanied by changes in floral colour that may influence pollinator preference and reproductive isolation. This study investigates whether changes in floral colour can accompany polyploid and homoploid hybridization, important processes in angiosperm evolution. METHODS: Spectral reflectance of corolla tissue was examined for 60 Nicotiana (Solanaceae) accessions (41 taxa) based on spectral shape (corresponding to pigmentation) as well as bee and hummingbird colour perception in order to assess patterns of floral colour evolution. Polyploid and homoploid hybrid spectra were compared with those of their progenitors to evaluate whether hybridization has resulted in floral colour shifts. KEY RESULTS: Floral colour categories in Nicotiana seem to have arisen multiple times independently during the evolution of the genus. Most younger polyploids displayed an unexpected floral colour, considering those of their progenitors, in the colour perception of at least one pollinator type, whereas older polyploids tended to resemble one or both of their progenitors. CONCLUSIONS: Floral colour evolution in Nicotiana is weakly constrained by phylogeny, and colour shifts do occur in association with both polyploid and homoploid hybrid divergence. Transgressive floral colour in N. tabacum has arisen by inheritance of anthocyanin pigmentation from its paternal progenitor while having a plastid phenotype like its maternal progenitor. Potentially, floral colour evolution has been driven by, or resulted in, pollinator shifts. However, those polyploids that are not sympatric (on a regional scale) with their progenitor lineages are typically not divergent in floral colour from them, perhaps because of a lack of competition for pollinators.

Fritillaria (Liliaceae) is a genus of approximately 140 species of bulbous perennial plants that includes taxa of both horticultural and medicinal importance. As well as being commercially valuable, Fritillaria species have attracted attention because of their exceptionally large genome sizes, with all values recorded to date in excess of 30 Gb. Despite such interest in the genus, phylogenetic relationships between the majority of species have remained untested. Here we present the first phylogenetic reconstruction of relationships to encompass most of the currently recognised species diversity in the genus. Three regions of the plastid genome were sequenced in 117 individuals of Fritillaria, representing 92 species (c. 66% of the genus) and in representatives of nine other genera of Liliaceae. Eleven low-copy nuclear gene regions were also screened in selected species for their potential utility. Phylogenetic analysis of a combined plastid dataset using maximum parsimony and Bayesian inference provided support for the monophyly of the majority of currently recognised subgenera. However, subgenus Fritillaria, which is by far the largest of the subgenera and includes the most important species used in traditional Chinese medicine, is found to be polyphyletic. Moreover, several taxa that were represented by multiple individuals show evidence of species non-monophyly. The Japanese endemic subgenus Japonica, which contains the species with the largest recorded genome size for any diploid plant, is resolved as sister to the predominantly Middle Eastern and Central Asian subgenus Rhinopetalum. Whilst relationships between most of the major Fritillaria lineages can now be resolved, our results also highlight the need for data from additional independently evolving loci; an endeavour that may be particularly challenging in light of the huge nuclear genomes found in these plants.

Nymphaeales are the most species-rich lineage of the earliest diverging angiosperms known as the ANA grade (Amborellales, Nymphaeales, Austrobaileyales), and they have received considerable attention from morphological, physiological, and ecological perspectives. Although phylogenetic relationships between these three lineages of angiosperms are mainly well resolved, insights at the whole genome level are still limited because of a dearth of information. To address this, genome sizes and chromosome numbers in 34 taxa, comprising 28 species were estimated and analysed together with previously published data to provide an overview of genome size and chromosome diversity in Nymphaeales. Overall, genome sizes were shown to vary 10-fold and chromosome numbers and ploidy levels ranged from 2n = 2x = 18 to 2n = 16x = ∼224. Distinct patterns of genome diversity were apparent, reflecting the differential incidence of polyploidy, changes in repetitive DNA content, and chromosome rearrangements within and between genera. Using model-based approaches, ancestral genome size and basic chromosome numbers were reconstructed to provide insights into the dynamics of genome size and chromosome number evolution. Finally, by combining additional data from Amborellales and Austrobaileyales, a comprehensive overview of genome sizes and chromosome numbers in these early diverging angiosperms is presented.

Interspecies relationships in Nicotiana (Solanaceae) are complex because 40 species are diploid (two sets of chromosomes) and 35 species are allotetraploid (four sets of chromosomes, two from each progenitor diploid species). We sequenced a fragment (containing four introns) of the nuclear gene ‘chloroplast-expressed glutamine synthetase’ (ncpGS) in 65 species of Nicotiana . Here we present the first phylogenetic analysis based on a low-copy nuclear gene for this well studied and important genus. Diploid species have a single-copy of ncpGS, and allotetraploids as expected have two homeologous copies, each derived from their progenitor diploid. Results were particularly useful for determining the paternal lineage of previously enigmatic taxa (for which our previous analyses had revealed only the maternal progenitors). In particular, we were able to shed light on the origins of the two oldest and largest allotetraploid sections, N . sects. Suaveolentes and Repandae . All homeologues have an intact reading frame and apparently similar rates of divergence, suggesting both remain functional. Difficulties in fitting certain diploid species into the sectional classification of Nicotiana on morphological grounds, coupled with discordance between the ncpGS data and previous trees (i.e. plastid, nuclear ribosomal DNA), indicate a number of homoploid (diploid) hybrids in the genus. We have evidence for Nicotiana glutinosa and Nicotiana linearis being of hybrid origin and patterns of intra-allelic recombination also indicate the possibility of reticulate origins for other diploid species.

DNA barcoding involves sequencing a standard region of DNA as a tool for species identification. However, there has been no agreement on which region(s) should be used for barcoding land plants. To provide a community recommendation on a standard plant barcode, we have compared the performance of 7 leading candidate plastid DNA regions ( atpF–atpH spacer, matK gene, rbcL gene, rpoB gene, rpoC1 gene, psbK–psbI spacer, and trnH–psbA spacer). Based on assessments of recoverability, sequence quality, and levels of species discrimination, we recommend the 2-locus combination of rbcL + matK as the plant barcode. This core 2-locus barcode will provide a universal framework for the routine use of DNA sequence data to identify specimens and contribute toward the discovery of overlooked species of land plants.

Clam shell morphology was related to lake water chemistry and sediment properties for Elliptio complanata from 40 south central Ontario lakes. Principal component analyses described the major trends in the environmental data. Relationships between shell morphology and the environmental trends were ascertained through a canonical correlation analysis. Thicker shells were related to increasing alkalinity, conductivity, and pH, suggesting that shell formation is favoured in well-buffered habitats that are not affected by acidic deposition. We also determined that narrow shells were related to coarse sediment in deeper water. This pattern of variation in shell morphology may be adaptive because streamlined shells would better enable clams to maintain their position in turbulent habitats. The water chemistry and sediment property correlates with shell morphology were independent of each other. Any attempt to use unionid shells as long-term indicators of changes in environmental water chemistry must be preceded by a consideration of other habitat variables such as water turbulence, sediment type, and depth.