Current status of global conservation and characterisation of wild and cultivated Brassicaceae genetic resources.

BACKGROUND: The economic importance of the globally distributed Brassicaceae family resides in the large diversity of crops within the family and the substantial variety of agronomic and functional traits they possess. We reviewed the current classifications of crop wild relatives (CWRs) in the Brassicaceae family with the aim of identifying new potential cross-compatible species from a total of 1,242 species using phylogenetic approaches. RESULTS: In general, cross-compatibility data between wild species and crops, as well as phenotype and genotype characterisation data, were available for major crops but very limited for minor crops, restricting the identification of new potential CWRs. Around 70% of wild Brassicaceae did not have genetic sequence data available in public repositories, and only 40% had chromosome counts published. Using phylogenetic distances, we propose 103 new potential CWRs for this family, which we recommend as priorities for cross-compatibility tests with crops and for phenotypic characterisation, including 71 newly identified CWRs for 10 minor crops. From the total species used in this study, more than half had no records of being in ex situ conservation, and 80% were not assessed for their conservation status or were data deficient (IUCN Red List Assessments). CONCLUSIONS: Great efforts are needed on ex situ conservation to have accessible material for characterising and evaluating the species for future breeding programmes. We identified the Mediterranean region as one key conservation area for wild Brassicaceae species, with great numbers of endemic and threatened species. Conservation assessments are urgently needed to evaluate most of these wild Brassicaceae.

Phylogenomics and the rise of the angiosperms.

Angiosperms are the cornerstone of most terrestrial ecosystems and human livelihoods1,2. A robust understanding of angiosperm evolution is required to explain their rise to ecological dominance. So far, the angiosperm tree of life has been determined primarily by means of analyses of the plastid genome3,4. Many studies have drawn on this foundational work, such as classification and first insights into angiosperm diversification since their Mesozoic origins5-7. However, the limited and biased sampling of both taxa and genomes undermines confidence in the tree and its implications. Here, we build the tree of life for almost 8,000 (about 60%) angiosperm genera using a standardized set of 353 nuclear genes8. This 15-fold increase in genus-level sampling relative to comparable nuclear studies9 provides a critical test of earlier results and brings notable change to key groups, especially in rosids, while substantiating many previously predicted relationships. Scaling this tree to time using 200 fossils, we discovered that early angiosperm evolution was characterized by high gene tree conflict and explosive diversification, giving rise to more than 80% of extant angiosperm orders. Steady diversification ensued through the remaining Mesozoic Era until rates resurged in the Cenozoic Era, concurrent with decreasing global temperatures and tightly linked with gene tree conflict. Taken together, our extensive sampling combined with advanced phylogenomic methods shows the deep history and full complexity in the evolution of a megadiverse clade.

Strategies of diaspore dispersal investment in Compositae: the case of the Andean highlands.

BACKGROUND AND AIMS: Understanding diaspore morphology and how much a species invests on dispersal appendages is key for improving our knowledge of dispersal in fragmented habitats. We investigate diaspore morphological traits in high-Andean Compositae and their main abiotic and biotic drivers and test whether they play a role in species distribution patterns across the naturally fragmented high-Andean grasslands. METHODS: We collected diaspore trait data for 125 Compositae species across 47 tropical high-Andean summits, focusing on achene length and pappus-to-achene length ratio, with the latter as a proxy of dispersal investment. We analysed the role of abiotic (temperature, elevation and latitude) and biotic factors (phylogenetic signal and differences between tribes) on diaspore traits and whether they are related to distribution patterns across the Andes, using phylogenomics, distribution modelling and community ecology analyses. KEY RESULTS: Seventy-five percent of the studied species show small achenes (length <3.3 mm) and 67% have high dispersal investment (pappus length at least two times the achene length). Dispersal investment increases with elevation, possibly to compensate for lower air density, and achene length increases towards the equator, where non-seasonal climate prevails. Diaspore traits show significant phylogenetic signal, and higher dispersal investment is observed in Gnaphalieae, Astereae and Senecioneae, which together represent 72% of our species. High-Andean-restricted species found across the tropical Andes have, on average, the pappus four times longer than the achene, a significantly higher dispersal investment than species present only in the northern Andes or only in the central Andes. CONCLUSIONS: Small achenes and high diaspore dispersal investment dominate among high-Andean Compositae, traits typical of mostly three tribes of African origin; but traits are also correlated with the environmental gradients within the high-Andean grasslands. Our results also suggest that diaspore dispersal investment is likely to shape species distribution patterns in naturally fragmented habitats.

A Bioinformatic Pipeline to Estimate Ploidy Level from Target Capture Sequence Data Obtained from Herbarium Specimens.

Whole genome duplications (WGD) are frequent in many plant lineages; however, ploidy level variation is unknown in most species. The most widely used methods to estimate ploidy levels in plants are chromosome counts, which require living specimens, and flow cytometry estimates, which necessitate living or relatively recently collected samples. Newly described bioinformatic methods have been developed to estimate ploidy levels using high-throughput sequencing data, and these have been optimized in plants by calculating allelic ratio values from target capture data. This method relies on the maintenance of allelic ratios from the genome to the sequence data. For example, diploid organisms will generate allelic data in a 1:1 proportion, with an increasing number of possible allelic ratio combinations occurring in individuals with higher ploidy levels. In this chapter, we explain step-by-step this bioinformatic approach for the estimation of ploidy level.

Effective population size in a partially clonal plant is not predicted by the number of genetic individuals.

Estimating effective population size (N e) is important for theoretical and practical applications in evolutionary biology and conservation. Nevertheless, estimates of N e in organisms with complex life-history traits remain scarce because of the challenges associated with estimation methods. Partially clonal plants capable of both vegetative (clonal) growth and sexual reproduction are a common group of organisms for which the discrepancy between the apparent number of individuals (ramets) and the number of genetic individuals (genets) can be striking, and it is unclear how this discrepancy relates to N e. In this study, we analysed two populations of the orchid Cypripedium calceolus to understand how the rate of clonal versus sexual reproduction affected N e. We genotyped >1000 ramets at microsatellite and SNP loci, and estimated contemporary N e with the linkage disequilibrium method, starting from the theoretical expectation that variance in reproductive success among individuals caused by clonal reproduction and by constraints on sexual reproduction would lower N e. We considered factors potentially affecting our estimates, including different marker types and sampling strategies, and the influence of pseudoreplication in genomic data sets on N e confidence intervals. The magnitude of N e/N ramets and N e/N genets ratios we provide may be used as reference points for other species with similar life-history traits. Our findings demonstrate that N e in partially clonal plants cannot be predicted based on the number of genets generated by sexual reproduction, because demographic changes over time can strongly influence N e. This is especially relevant in species of conservation concern in which population declines may not be detected by only ascertaining the number of genets.

Plastid phylogenomics and molecular evolution of Thismiaceae (Dioscoreales).

PREMISE: Species in Thismiaceae can no longer photosynthesize and instead obtain carbon from soil fungi. Here we infer Thismiaceae phylogeny using plastid genome data and characterize the molecular evolution of this genome. METHODS: We assembled five Thismiaceae plastid genomes from genome skimming data, adding to previously published data for phylogenomic inference. We investigated plastid-genome structural changes, considering locally colinear blocks (LCBs). We also characterized possible shifts in selection pressure in retained genes by considering changes in the ratio of nonsynonymous to synonymous changes (ω). RESULTS: Thismiaceae experienced two major pulses of gene loss around the early diversification of the family, with subsequent scattered gene losses across descendent lineages. In addition to massive size reduction, Thismiaceae plastid genomes experienced occasional inversions, and there were likely two independent losses of the plastid inverted repeat (IR) region. Retained plastid genes remain under generally strong purifying selection (ω << 1), with significant and sporadic weakening or strengthening in several instances. The bifunctional trnE-UUC gene of Thismia huangii may retain a secondary role in heme biosynthesis, despite a probable loss of functionality in protein translation. Several cis-spliced group IIA introns have been retained, despite the loss of the plastid intron maturase, matK. CONCLUSIONS: We infer that most gene losses in Thismiaceae occurred early and rapidly, following the initial loss of photosynthesis in its stem lineage. As a species-rich, fully mycoheterotrophic lineage, Thismiaceae provide a model system for uncovering the unique and divergent ways in which plastid genomes evolve in heterotrophic plants.

Multiple pre- and postzygotic components of reproductive isolation between two co-occurring Lysimachia species.

Genetic divergence between species depends on reproductive isolation (RI) due to traits that reduce interspecific mating (prezygotic isolation) or are due to reduced hybrid fitness (postzygotic isolation). Previous research found that prezygotic barriers tend to be stronger than postzygotic barriers, but most studies are based on the evaluation of F1 hybrid fitness in early life cycle stages. We combined field and experimental data to determine the strength of 17 prezygotic and postzygotic reproductive barriers between two Lysimachia species that often co-occur and share pollinators. We assessed postzygotic barriers up to F2 hybrids and backcrosses. The two species showed near complete RI due to the cumulative effect of multiple barriers, with an uneven and asymmetric contribution to isolation. In allopatry, prezygotic barriers contributed more to reduce gene flow than postzygotic barriers, but their contributions were more similar in sympatry. The strength of postzygotic RI was up to three times lower for F1 progeny than for F2 or backcrossed progenies, and RI was only complete when late F1 stages and either F2 or backcrosses were accounted for. Our results thus suggest that the relative strength of postzygotic RI may be underestimated when its effects on late stages of the life cycle are disregarded.

Genomic, spatial and morphometric data for discrimination of four species in the Mediterranean Tamus clade of yams (Dioscorea, Dioscoreaceae).

BACKGROUND AND AIMS: Among the numerous pantropical species of the yam genus, Dioscorea, only a small group occurs in the Mediterranean basin, including two narrow Pyrenean endemics (Borderea clade) and two Mediterranean-wide species (D. communis and D. orientalis, Tamus clade). However, several currently unrecognized species and infraspecific taxa have been described in the Tamus clade due to significant morphological variation associated with D. communis. Our overarching aim was to investigate taxon delimitation in the Tamus clade using an integrative approach combining phylogenomic, spatial and morphological data. METHODS: We analysed 76 herbarium samples using Hyb-Seq genomic capture to sequence 260 low-copy nuclear genes and plastomes, together with morphometric and environmental modelling approaches. KEY RESULTS: Phylogenomic reconstructions confirmed that the two previously accepted species of the Tamus clade, D. communis and D. orientalis, are monophyletic and form sister clades. Three subclades showing distinctive geographic patterns were identified within D. communis. These subclades were also identifiable from morphometric and climatic data, and introgression patterns were inferred between subclades in the eastern part of the distribution of D. communis. CONCLUSIONS: We propose a taxonomy that maintains D. orientalis, endemic to the eastern Mediterranean region, and splits D. communis sensu lato into three species: D. edulis, endemic to Macaronesia (Canary Islands and Madeira); D. cretica, endemic to the eastern Mediterranean region; and D. communis sensu stricto, widespread across western and central Europe. Introgression inferred between D. communis s.s. and D. cretica is likely to be explained by their relatively recent speciation at the end of the Miocene, disjunct isolation in eastern and western Mediterranean glacial refugia and a subsequent westward recolonization of D. communis s.s. Our study shows that the use of integrated genomic, spatial and morphological approaches allows a more robust definition of species boundaries and the identification of species that previous systematic studies failed to uncover.

Genome-wide footprints in the carob tree (Ceratonia siliqua) unveil a new domestication pattern of a fruit tree in the Mediterranean.

Intense research efforts over the last two decades have renewed our understanding of plant phylogeography and domestication in the Mediterranean basin. Here we aim to investigate the evolutionary history and the origin of domestication of the carob tree (Ceratonia siliqua), which has been cultivated for millennia for food and fodder. We used >1000 microsatellite genotypes to delimit seven carob evolutionary units (CEUs). We investigated genome-wide diversity and evolutionary patterns of the CEUs with 3557 single nucleotide polymorphisms generated by restriction-site associated DNA sequencing (RADseq). To address the complex wild vs. cultivated status of sampled trees, we classified 56 sampled populations across the Mediterranean basin as wild, seminatural or cultivated. Nuclear and cytoplasmic loci were identified from RADseq data and separated for analyses. Phylogenetic analyses of these genomic-wide data allowed us to resolve west-to-east expansions from a single long-term refugium probably located in the foothills of the High Atlas Mountains near the Atlantic coast. Our findings support multiple origins of domestication with a low impact on the genetic diversity at range-wide level. The carob was mostly domesticated from locally selected wild genotypes and scattered long-distance westward dispersals of domesticated varieties by humans, concomitant with major historical migrations by Romans, Greeks and Arabs. Ex situ efforts to preserve carob genetic resources should prioritize accessions from both western and eastern populations, with emphasis on the most differentiated CEUs situated in southwest Morocco, south Spain and eastern Mediterranean. Our study highlights the relevance of wild and seminatural habitats in the conservation of genetic resources for cultivated trees.

Genome Insights into Autopolyploid Evolution: A Case Study in Senecio doronicum (Asteraceae) from the Southern Alps.

Polyploidy is a widespread phenomenon across angiosperms, and one of the main drivers of diversification. Whilst it frequently involves hybridisation, autopolyploidy is also an important feature of plant evolution. Minority cytotypes are frequently overlooked due to their lower frequency in populations, but the development of techniques such as flow cytometry, which enable the rapid screening of cytotype diversity across large numbers of individuals, is now providing a more comprehensive understanding of cytotype diversity within species. Senecio doronicum is a relatively common daisy found throughout European mountain grasslands from subalpine to almost nival elevations. We have carried out a population-level cytotype screening of 500 individuals from Tête Grosse (Alpes-de-Haute-Provence, France), confirming the coexistence of tetraploid (28.2%) and octoploid cytotypes (71.2%), but also uncovering a small number of hexaploid individuals (0.6%). The analysis of repetitive elements from short-read genome-skimming data combined with nuclear (ITS) and whole plastid DNA sequences support an autopolyploid origin of the polyploid S. doronicum individuals and provide molecular evidence regarding the sole contribution of tetraploids in the formation of hexaploid individuals. The evolutionary impact and resilience of the new cytotype have yet to be determined, although the coexistence of different cytotypes may indicate nascent speciation.