Resolving the phylogeny of Thladiantha (Cucurbitaceae) with three different target capture pipelines.

BACKGROUND: Despite recent advances, reliable tools to simultaneously handle different types of sequencing data (e.g., target capture, genome skimming) for phylogenomics are still scarce. Here, we evaluate the performance of the recently developed pipeline Captus in comparison with the well-known target capture pipelines HybPiper and SECAPR. As test data, we analyzed newly generated sequences for the genus Thladiantha (Cucurbitaceae) for which no well-resolved phylogeny estimate has been available so far, as well as simulated reads derived from the genome of Arabidopsis thaliana. RESULTS: Our pipeline comparisons are based on (1) the time needed for data assembly and locus extraction, (2) locus recovery per sample, (3) the number of informative sites in nucleotide alignments, and (4) the topology of the nuclear and plastid phylogenies. Additionally, the simulated reads derived from the genome of Arabidopsis thaliana were used to evaluate the accuracy and completeness of the recovered loci. In terms of computation time, locus recovery per sample, and informative sites, Captus outperforms HybPiper and SECAPR. The resulting topologies of Captus and SECAPR are identical for coalescent trees but differ when trees are inferred from concatenated alignments. The HybPiper phylogeny is similar to Captus in both methods. The nuclear genes recover a deep split of Thladiantha in two clades, but this is not supported by the plastid data. CONCLUSIONS: Captus is the best choice among the three pipelines in terms of computation time and locus recovery. Even though there is no significant topological difference between the Thladiantha species trees produced by the three pipelines, Captus yields a higher number of gene trees in agreement with the topology of the species tree (i.e., fewer genes in conflict with the species tree topology).

Pan-European study of genotypes and phenotypes in the Arabidopsis relative Cardamine hirsuta reveals how adaptation, demography, and development shape diversity patterns.

We study natural DNA polymorphisms and associated phenotypes in the Arabidopsis relative Cardamine hirsuta. We observed strong genetic differentiation among several ancestry groups and broader distribution of Iberian relict strains in European C. hirsuta compared to Arabidopsis. We found synchronization between vegetative and reproductive development and a pervasive role for heterochronic pathways in shaping C. hirsuta natural variation. A single, fast-cycling ChFRIGIDA allele evolved adaptively allowing range expansion from glacial refugia, unlike Arabidopsis where multiple FRIGIDA haplotypes were involved. The Azores islands, where Arabidopsis is scarce, are a hotspot for C. hirsuta diversity. We identified a quantitative trait locus (QTL) in the heterochronic SPL9 transcription factor as a determinant of an Azorean morphotype. This QTL shows evidence for positive selection, and its distribution mirrors a climate gradient that broadly shaped the Azorean flora. Overall, we establish a framework to explore how the interplay of adaptation, demography, and development shaped diversity patterns of 2 related plant species.

Variability in mycorrhizal status of plant species is much larger within than between plots in grassland and coastal habitats.

Community-level studies linking plant mycorrhizal status to environment usually do not account for within-plot mycorrhizal status variability; thus, patterns of plant mycorrhizal status diversity are largely unknown. Here, we assessed the relative importance of within- and between-plot variability components in mycorrhizal status and examined how plant mycorrhizal status diversity is related to soil nutrient availability. We hypothesised larger between-plot variability in mycorrhizal status and higher plant mycorrhizal status diversity in P-poor soils. To test these hypotheses, we used plant phylogenies, vegetation, soil and plant mycorrhizal status data from Czech semi-natural grasslands and Scottish coastal habitats. We divided plant mycorrhizal status diversity into divergence and evenness and tested their relations to soil P, K, Ca and Mg. Within-plot variability component of mycorrhizal status was always, on average, at least 2.2 times larger than between-plot variability in our datasets. Plant mycorrhizal status divergence was positively related to Ca (in both datasets) and Mg (only in grasslands and when accounting for phylogeny). In grasslands, the relationship between Mg and plant mycorrhizal status evenness was negative when accounting for phylogeny, while it was positive when not accounting for phylogeny. Plant mycorrhizal status diversity was not linked to P and its relation to K was inconsistent. Our results suggest that high Ca in the soil can promote coexistence of mycorrhizal, facultatively mycorrhizal and non-mycorrhizal plant species. We encourage future studies to also focus on within-plot variability in mycorrhizal status, because it appears to be highly relevant in herbaceous systems.

Multiplexed shotgun genotyping (MSG) data resolve phylogenetic relationships within and among archipelagos in Macaronesian Tolpis.

PREMISE: Plants endemic to oceanic archipelagos are suitable for studying evolution, being isolated on substrates of different ages. Evolution has been recent, rendering traditionally employed sequences insufficiently variable for resolving relationships. This study includes sampling in the genus Tolpis (Asteraceae) from the Azores, Madeira, and Cape Verde, and expands upon an earlier study demonstrating the efficacy of multiplexed shotgun genotyping (MSG) for resolving relationships in Canarian Tolpis. METHODS: Genomic libraries for 90 accessions of Tolpis and two from the outgroup were generated for genotyping individuals using MSG. Loci were de novo assembled with iPyrad, which clusters MSG loci within and between samples. A maximum likelihood phylogeny was generated with RAxML. Ancestral area reconstruction was inferred using R package BioGeoBEARS. RESULTS: MSG data recovered a highly resolved phylogeny from population to inter-archipelago levels. Ancestral area reconstruction provided biogeographic hypotheses for the radiation of Macaronesian Tolpis. CONCLUSIONS: Four major clades were resolved. The Madeiran endemic T. macrorhiza is sister to other Tolpis. Species from the Canaries, Cape Verdes, and the continent are sister to T. succulenta from Madeira, which has a sister subclade of Azorean populations composed of T. succulenta and T. azorica. Population-level resolution suggests unrecognized taxa on several archipelagos. Ancestral reconstruction suggests initial dispersal from the continent to Madeira, with dispersal to the Azores, then dispersal from Madeira to the Canary Islands, with both subsequent dispersal to the Cape Verdes and back-dispersal to the continent. Single-island radiations and inter-island dispersal are implicated in divergence in Macaronesian Tolpis.

Mycorrhizal status is a poor predictor of the distribution of herbaceous species along the gradient of soil nutrient availability in coastal and grassland habitats.

Plant mycorrhizal status (a trait indicating the ability to form mycorrhizas) can be a useful plant trait for predicting changes in vegetation influenced by increased fertility. Mycorrhizal fungi enhance nutrient uptake and are expected to provide a competitive advantage for plants growing in nutrient-poor soils; while in nutrient-rich soils, mycorrhizal symbiosis may be disadvantageous. Some studies in natural systems have shown that mycorrhizal plants can be more frequent in P and N-poor soils (low nutrient availability) or Ca and Mg-high (high pH) soils, but empirical support is still not clear. Using vegetation and soil data from Scottish coastal habitats, and Latvian and Czech grasslands, we examined whether there is a link between plant mycorrhizal status and plant-available P, N, Ca and Mg. We performed the max test analysis (to examine the central tendency) and a combination of quantile regression and meta-analysis (to examine tendencies in different quantiles) on both community and plant species data combined with plant phylogenies. We consistently found no changes in mycorrhizal status at the community and species levels along the gradients of plant-available P, N, Ca and Mg in the central tendency and in almost all quantiles across all datasets. Thus, we found no support for the hypotheses that herbaceous species which are able to form mycorrhizas are more frequent in nutrient-poor and high pH environments. Obligatory, facultatively and non-mycorrhizal herbaceous species appear to assemble randomly along the gradients of nutrient availability in several European herbaceous habitats, suggesting that all these strategies perform similarly under non-extreme soil nutrient conditions.

Enhancing photosynthesis at high light levels by adaptive laboratory evolution.

Photosynthesis is readily impaired by high light (HL) levels. Photosynthetic organisms have therefore evolved various mechanisms to cope with the problem. Here, we have dramatically enhanced the light tolerance of the cyanobacterium Synechocystis by adaptive laboratory evolution (ALE). By combining repeated mutagenesis and exposure to increasing light intensities, we generated strains that grow under extremely HL intensities. HL tolerance was associated with more than 100 mutations in proteins involved in various cellular functions, including gene expression, photosynthesis and metabolism. Co-evolved mutations were grouped into five haplotypes, and putative epistatic interactions were identified. Two representative mutations, introduced into non-adapted cells, each confer enhanced HL tolerance, but they affect photosynthesis and respiration in different ways. Mutations identified by ALE that allow photosynthetic microorganisms to cope with altered light conditions could be employed in assisted evolution approaches and could strengthen the robustness of photosynthesis in crop plants.

Phylogenetic informativeness analyses to clarify past diversification processes in Cucurbitaceae.

Phylogenomic studies have so far mostly relied on genome skimming or target sequence capture, which suffer from representation bias and can fail to resolve relationships even with hundreds of loci. Here, we explored the potential of phylogenetic informativeness and tree confidence analyses to interpret phylogenomic datasets. We studied Cucurbitaceae because their small genome size allows cost-efficient genome skimming, and many relationships in the family remain controversial, preventing inferences on the evolution of characters such as sexual system or floral morphology. Genome skimming and PCR allowed us to retrieve the plastome, 57 single copy nuclear genes, and the nuclear ribosomal ITS from 29 species representing all but one tribe of Cucurbitaceae. Node support analyses revealed few inter-locus conflicts but a pervasive lack of phylogenetic signal among plastid loci, suggesting a fast divergence of Cucurbitaceae tribes. Data filtering based on phylogenetic informativeness and risk of homoplasy clarified tribe-level relationships, which support two independent evolutions of fringed petals in the family. Our study illustrates how formal analysis of phylogenomic data can increase our understanding of past diversification processes. Our data and results will facilitate the design of well-sampled phylogenomic studies in Cucurbitaceae and related families.

Origin and domestication of Cucurbitaceae crops: insights from phylogenies, genomics and archaeology.

Some of the World's most valuable crops, including watermelon, honey melon, cucumber, squash, zucchini and pumpkin, belong to the family Cucurbitaceae. We review insights on their domestication from new phylogenies, archaeology and genomic studies. Ancestral state estimation on the most complete Cucurbitaceae phylogeny to date suggests that an annual life cycle may have contributed to domestication. Domestication started c. 11 000 years ago in the New World and Asia, and apparently more recently in Africa. Some cucurbit crops were domesticated only once, others multiple times (e.g. melon from different Asian and African populations). Most wild cucurbit fruits are bitter and nonpalatable to humans, and nonbitterness of the pulp apparently was a trait favoured early during domestication, with genomic data showing how bitterness loss was achieved convergently. The genetic pathways underlying lycopene accumulation, red or orange pulp colour, and fruit size and shape are only just beginning to be understood. The study of cucurbit domestication in recent years has benefitted from the increasing integration of archaeological and genomic data with insights from herbarium collections, the most efficient way to understand species' natural geographic ranges and climate adaptations.

How rapidly do self-compatible populations evolve selfing? Mating system estimation within recently evolved self-compatible populations of Azorean Tolpis succulenta (Asteraceae).

Genome-wide genotyping and Bayesian inference method (BORICE) were employed to estimate outcrossing rates and paternity in two small plant populations of Tolpis succulenta (Asteraceae) on Graciosa island in the Azores. These two known extant populations of T. succulenta on Graciosa have recently evolved self-compatibility. Despite the expectation that selfing would occur at an appreciable rate (self-incompatible populations of the same species show low but nonzero selfing), high outcrossing was found in progeny arrays from maternal plants in both populations. This is inconsistent with an immediate transition to high selfing following the breakdown of a genetic incompatibility system. This finding is surprising given the small population sizes and the recent colonization of an island from self-incompatible colonists of T. succulenta from another island in the Azores, and a potential paucity of pollinators, all factors selecting for selfing through reproductive assurance. The self-compatible lineage(s) likely have high inbreeding depression (ID) that effectively halts the evolution of increased selfing, but this remains to be determined. Like their progeny, all maternal plants in both populations are fully outbred, which is consistent with but not proof of high ID. High multiple paternity was found in both populations, which may be due in part to the abundant pollinators observed during the flowering season.

A comprehensive genome variation map of melon identifies multiple domestication events and loci influencing agronomic traits.

Melon is an economically important fruit crop that has been cultivated for thousands of years; however, the genetic basis and history of its domestication still remain largely unknown. Here we report a comprehensive map of the genomic variation in melon derived from the resequencing of 1,175 accessions, which represent the global diversity of the species. Our results suggest that three independent domestication events occurred in melon, two in India and one in Africa. We detected two independent sets of domestication sweeps, resulting in diverse characteristics of the two subspecies melo and agrestis during melon breeding. Genome-wide association studies for 16 agronomic traits identified 208 loci significantly associated with fruit mass, quality and morphological characters. This study sheds light on the domestication history of melon and provides a valuable resource for genomics-assisted breeding of this important crop.

Morphological and Molecular Characterization of Zanthoxylum zanthoxyloides (Rutaceae) from Burkina Faso.

Zanthoxylum zanthoxyloides is a West African forest tree that is used for example against malaria and sickle cell anemia in Burkina Faso. The goal of this study was to analyze the genetic and morphological diversity of the species within wild populations in Burkina Faso, where it is potentially under threat due to the uncontrolled harvesting of its roots. Seventy-two trees from three different sites in Southwestern Burkina Faso were analyzed. Each tree was characterized by 12 traits specifying the period of flowering and maturity as well as morphological characteristics of the stem, leaves, and seeds. The molecular analysis was performed using two plastid DNA regions (psbA-trnH and trnL-trnF) and two nuclear regions (GBSSI and ITS) to identify the genetic diversity of the species for further development of a management plan for ex situ reproduction and in situ conservation. We found variability in morphological traits correlating with the geographic distance of the study sites. The molecular analysis, in contrast, revealed hardly any genetic variability among the tested trees and no population structure. Whether the differences in morphological traits are caused by different environmental conditions or by genetic variability in genes linked to morphological traits needs further testing. The apparent lack of genetic differentiation suggests that germplasm throughout the study region is suitable for planting in conservation actions. Efficient conservation management should involve local communities, especially those interested in traditional medicine.

European mushroom assemblages are darker in cold climates.

Thermal melanism theory states that dark-colored ectotherm organisms are at an advantage at low temperature due to increased warming. This theory is generally supported for ectotherm animals, however, the function of colors in the fungal kingdom is largely unknown. Here, we test whether the color lightness of mushroom assemblages is related to climate using a dataset of 3.2 million observations of 3,054 species across Europe. Consistent with the thermal melanism theory, mushroom assemblages are significantly darker in areas with cold climates. We further show differences in color phenotype between fungal lifestyles and a lifestyle differentiated response to seasonality. These results indicate a more complex ecological role of mushroom colors and suggest functions beyond thermal adaption. Because fungi play a crucial role in terrestrial carbon and nutrient cycles, understanding the links between the thermal environment, functional coloration and species' geographical distributions will be critical in predicting ecosystem responses to global warming.

An Ancient Peptide Family Buried within Vicilin Precursors.

New proteins can evolve by duplication and divergence or de novo, from previously noncoding DNA. A recently observed mechanism is for peptides to evolve within a "host" protein and emerge by proteolytic processing. The first examples of such interstitial peptides were ones hosted by precursors for seed storage albumin. Interstitial peptides have also been observed in precursors for seed vicilins, but current evidence for vicilin-buried peptides (VBPs) is limited to seeds of the broadleaf plants pumpkin and macadamia. Here, an extensive sequence analysis of vicilin precursors suggested that peptides buried within the N-terminal region of preprovicilins are widespread and truly ancient. Gene sequences indicative of interstitial peptides were found in species from Amborellales to eudicots and include important grass and legume crop species. We show the first protein evidence for a monocot VBP in date palm seeds as well as protein evidence from other crops including the common tomato, sesame and pumpkin relatives, cucumber, and the sponge loofah ( Luffa aegyptiaca). Their excision was consistent with asparaginyl endopeptidase-mediated maturation, and sequences were confirmed by tandem mass spectrometry. Our findings suggest that the family is large and ancient and that based on the NMR solution structures for loofah Luffin P1 and tomato VBP-8, VBPs adopt a helical hairpin fold stapled by two internal disulfide bonds. The first VBPs characterized were a protease inhibitor, antimicrobials, and a ribosome inactivator. The age and evolutionary retention of this peptide family suggest its members play important roles in plant biology.

Repeated domestication of melon (Cucumis melo) in Africa and Asia and a new close relative from India.

PREMISE OF THE STUDY: The domestication history of melon is still unclear. An African or Asian origin has been suggested, but its closest wild relative was recently revealed to be an Australian species. The complicated taxonomic history of melon has resulted in additional confusion, with a high number of misidentified germplasm collections currently used by breeders and in genomics research. METHODS: Using seven DNA regions sequenced for 90% of the genus and the major cultivar groups, we sort out described names and infer evolutionary origins and domestication centers. KEY RESULTS: We found that modern melon cultivars go back to two lineages, which diverged ca. 2 million years ago. One is restricted to Asia (Cucumis melo subsp. melo), and the second, here described as C. melo subsp. meloides, is restricted to Africa. The Asian lineage has given rise to the widely commercialized cultivar groups and their market types, while the African lineage gave rise to cultivars still grown in the Sudanian region. We show that C. trigonus, an overlooked perennial and drought-tolerant species from India is among the closest living relatives of C. melo. CONCLUSIONS: Melon was domesticated at least twice: in Africa and Asia. The African lineage and the Indian C. trigonus are exciting new resources for breeding of melons tolerant to climate change.

Evolutionary dynamics of host specialization in wood-decay fungi.

BACKGROUND: The majority of wood decomposing fungi are mushroom-forming Agaricomycetes, which exhibit two main modes of plant cell wall decomposition: white rot, in which all plant cell wall components are degraded, including lignin, and brown rot, in which lignin is modified but not appreciably removed. Previous studies suggested that brown rot fungi tend to be specialists of gymnosperm hosts and that brown rot promotes gymnosperm specialization. However, these hypotheses were based on analyses of limited datasets of Agaricomycetes. Overcoming this limitation, we used a phylogeny with 1157 species integrating available sequences, assembled decay mode characters from the literature, and coded host specialization using the newly developed R package, rusda. RESULTS: We found that most brown rot fungi are generalists or gymnosperm specialists, whereas most white rot fungi are angiosperm specialists. A six-state model of the evolution of host specialization revealed high transition rates between generalism and specialization in both decay modes. However, while white rot lineages switched most frequently to angiosperm specialists, brown rot lineages switched most frequently to generalism. A time-calibrated phylogeny revealed that Agaricomycetes is older than the flowering plants but many of the large clades originated after the diversification of the angiosperms in the Cretaceous. CONCLUSIONS: Our results challenge the current view that brown rot fungi are primarily gymnosperm specialists and reveal intensive white rot specialization to angiosperm hosts. We thus suggest that brown rot associated convergent loss of lignocellulose degrading enzymes was correlated with host generalism, rather than gymnosperm specialism. A likelihood model of host specialization evolution together with a time-calibrated phylogeny further suggests that the rise of the angiosperms opened a new mega-niche for wood-decay fungi, which was exploited particularly well by white rot lineages.

Phylogeography of the Macaronesian Lettuce Species Lactuca watsoniana and L. palmensis (Asteraceae).

The phylogenetic relationships and phylogeography of two relatively rare Macaronesian Lactuca species, Lactuca watsoniana (Azores) and L. palmensis (Canary Islands), were, until this date, unclear. Karyological information of the Azorean species was also unknown. For this study, a chromosome count was performed and L. watsoniana showed 2n = 34. A phylogenetic approach was used to clarify the relationships of the Azorean endemic L. watsoniana and the La Palma endemic L. palmensis within the subtribe Lactucinae. Maximum parsimony, Maximum likelihood and Bayesian analysis of a combined molecular dataset (ITS and four chloroplast DNA regions) and molecular clock analyses were performed with the Macaronesian Lactuca species, as well as a TCS haplotype network. The analyses revealed that L. watsoniana and L. palmensis belong to different subclades of the Lactuca clade. Lactuca watsoniana showed a strongly supported phylogenetic relationship with North American species, while L. palmensis was closely related to L. tenerrima and L. inermis, from Europe and Africa. Lactuca watsoniana showed four single-island haplotypes. A divergence time estimation of the Macaronesian lineages was used to examine island colonization pathways. Results obtained with BEAST suggest a divergence of L. palmensis and L. watsoniana clades c. 11 million years ago, L. watsoniana diverged from its North American sister species c. 3.8 million years ago and L. palmensis diverged from its sister L. tenerrima, c. 1.3 million years ago, probably originating from an African ancestral lineage which colonized the Canary Islands. Divergence analyses with *BEAST indicate a more recent divergence of the L. watsoniana crown, c. 0.9 million years ago. In the Azores colonization, in a stepping stone, east-to-west dispersal pattern, associated with geological events might explain the current distribution range of L. watsoniana.

Bees of the Azores: an annotated checklist (Apidae, Hymenoptera).

We report 18 species of wild bees plus the domesticated honeybee from the Azores, which adds nine species to earlier lists. One species, Hylaeus azorae, seems to be a single island endemic, and three species are possibly native (Colletes eous, Halictus villosulus, and Hylaeus pictipes). All the remaining bee species are most likely accidental introductions that arrived after human colonization of the archipelago in the 15th century. Bee diversity in the Azores is similar to bee diversity of Madeira and Cape Verde but nearly ten times lower than it is in the Canary Islands.

Are there any widespread endemic flowering plant species in Macaronesia? Phylogeography of Ranunculus cortusifolius.

PREMISE OF THE STUDY: Oceanic island endemics typically exhibit very restricted distributions. In Macaronesia, only one endemic angiosperm species, Ranunculus cortusifolius, has a distribution spanning the archipelagos of the Azores, Madeira, and Canaries. Earlier work suggested possible differences between archipelagos and the multiple origins of the species. This paper tests the hypothesis that R. cortusifolius is a single widespread Macaronesian endemic species with a single origin. METHODS: Chloroplast (matK-trnK, psbJ-petA) and ITS sequences were generated from across the distribution of R. cortusifolius. Relationships were investigated using Bayesian inference and divergence times estimated using BEAST. Infraspecific variation was investigated using statistical parsimony. The general mixed Yule-coalescent model (GMYC) was further used to identify putative species boundaries based on maternally inherited plastid data. KEY RESULTS: The hypothesis of multiple independent origins of R. cortusifolius is rejected. Divergence of the R. cortusifolius lineage from a western Mediterranean sister group in the late Miocene is inferred. Distinct genotypes were resolved within R. cortusifolius that are endemic to the Azores, Madeira, and the Canaries. Four to five putative species were delimited by different versions of the GMYC model. CONCLUSION: Ranunculus cortusifolius is the result of a single colonization of Macaronesia. The large distances between archipelagos have been effective barriers to dispersal, promoting allopatric diversification at the molecular level with diversification also evident within the Canaries. Isolation has not been accompanied by marked morphological diversification, which may be explained by the typical association of R. cortusifolius with stable and climatically buffered laurel forest communities.

Towards a more holistic research approach to plant conservation: the case of rare plants on oceanic islands.

Research dedicated to rare endemic plants is usually focused on one given aspect. However, holistic studies, addressing several key issues, might be more useful, supporting management programmes while unravelling basic knowledge about ecological and population-level processes. A more comprehensive approach to research is proposed, encompassing: phylogenetics/systematics, pollination biology and seed dispersal, propagation, population genetics, species distribution models (SDMs), threats and monitoring. We present a holistic study dedicated to Veronica dabneyi Hochst. ex Seub., an endangered chamaephyte endemic to the Azores. Veronica dabneyi was mainly found associated with other endemic taxa; however, invasive plants were also present and together with introduced cattle, goats and rabbits are a major threat. Most populations grow at somewhat rocky and steep locations that appeared to work as refuges. Seed set in the wild was generally high and recruitment of young plants from seed seemed to be frequent. In the laboratory, it was possible to germinate and fully develop V. dabneyi seedlings, which were planted at their site of origin. No dormancy was detected and time for 50 % germination was affected by incubation temperature. Eight new microsatellite markers were applied to 72 individuals from 7 sites. A considerable degree of admixture was found between samples from the two islands Flores and Corvo, with 98 % of the genetic variability allocated within populations. Levels of heterozygosity were high and no evidence of inbreeding was found. Species distribution models based on climatic and topographic variables allowed the estimation of the potential distribution of V. dabneyi on Flores and Corvo using ecological niche factor analysis and Maxent. The inclusion of land-use variables only slightly increased the information explained by the models. Projection of the expected habitat in Faial largely coincided with the only historic record of V. dabneyi on that island. This research could be the basis for the design of a recovery plan, showing the pertinence of more holistic research approaches to plant conservation.