The reliability of published scientific papers has been the topic of much recent discussion, notably in the biomedical sciences [1]. Although small sample size is regularly pointed as one of the culprits, big data can also be a concern. The advent of high-throughput sequencing, and the processing of sequence data by opaque bioinformatics workflows, mean that sequences with often high error rates are produced, and that exact but slow analyses are not feasible.
The troubles with bioinformatics arise from the increased complexity of the tools used by scientists, and from the lack of incentives and/or skills from authors (but also reviewers and editors) to make sure of the quality of those tools. As a much discussed example, a bug in the widely used PLINK software [2] has been pointed as the explanation [3] for incorrect inference of selection for increased height in European Human populations [4].
High-throughput sequencing often generates high rates of genotyping errors, so that the development of bioinformatics tools to assess the quality of data and correct them is a major issue. The work of Gillingham et al. [5] contributes to the latter goal. In this work, the authors propose a new bioinformatics workflow (ACACIA) for performing genotyping analysis of multigene complexes, such as self-incompatibility genes in plants, major histocompatibility genes (MHC) in vertebrates, and homeobox genes in animals, which are particularly challenging to genotype in non-model organisms. PCR and sequencing of multigene families generate artefacts, hence spurious alleles. A key to Gillingham et al.‘ s method is to call candidate genes based on Oligotyping, a software pipeline originally conceived for identifying variants from microbiome 16S rRNA amplicons [6]. This allows to reduce the number of false positives and the number of dropout alleles, compared to previous workflows.
This method is not based on an explicit probability model, and thus it is not conceived to provide a control of the rate of errors as, say, a valid confidence interval should (a confidence interval with coverage c for a parameter should contain the parameter with probability c, so the error rate 1- c is known and controlled by the user who selects the value of c). However, the authors suggest a method to adapt the settings of ACACIA to each application.
To compare and validate the new workflow, the authors have constructed new sets of genotypes representing different extents copy number variation, using already known genotypes from chicken MHC. In such conditions, it was possible to assess how many alleles are not detected and what is the rate of false positives. Gillingham et al. additionally investigated the effect of using non-optimal primers. They found better performance of ACACIA compared to a preexisting pipeline, AmpliSAS [7], for optimal settings of both methods. However, they do not claim that ACACIA will always be better than AmpliSAS. Rather, they warn against the common practice of using the default settings of the latter pipeline. Altogether, this work and the ACACIA workflow should allow for better ascertainment of genotypes from multigene families.

References

[1] Ioannidis, J. P. A, Greenland, S., Hlatky, M. A., Khoury, M. J., Macleod, M. R., Moher, D., Schulz, K. F. and Tibshirani, R. (2014) Increasing value and reducing waste in research design, conduct, and analysis. The Lancet, 383, 166-175. doi: 10.1016/S0140-6736(13)62227-8
[2] Chang, C. C., Chow, C. C., Tellier, L. C. A. M., Vattikuti, S., Purcell, S. M. and Lee, J. J. (2015) Second-generation PLINK: rising to the challenge of larger and richer datasets. GigaScience, 4, 7, s13742-015-0047-8. doi: 10.1186/s13742-015-0047-8
[3] Robinson, M. R. and Visscher, P. (2018) Corrected sibling GWAS data release from Robinson et al. http://cnsgenomics.com/data.html
[4] Field, Y., Boyle, E. A., Telis, N., Gao, Z., Gaulton, K. J., Golan, D., Yengo, L., Rocheleau, G., Froguel, P., McCarthy, M.I . and Pritchard J. K. (2016) Detection of human adaptation during the past 2000 years. Science, 354(6313), 760-764. doi: 10.1126/science.aag0776
[5] Gillingham, M. A. F., Montero, B. K., Wihelm, K., Grudzus, K., Sommer, S. and Santos P. S. C. (2020) A novel workflow to improve multi-locus genotyping of wildlife species: an experimental set-up with a known model system. bioRxiv 638288, ver. 3 peer-reviewed and recommended by Peer Community In Evolutionary Biology. doi: 10.1101/638288
[6] Eren, A. M., Maignien, L., Sul, W. J., Murphy, L. G., Grim, S. L., Morrison, H. G., and Sogin, M.L. (2013) Oligotyping: differentiating between closely related microbial taxa using 16S rRNA gene data. Methods in Ecology and Evolution 4(12), 1111-1119. doi: 10.1111/2041-210X.12114
[7] Sebastian, A., Herdegen, M., Migalska, M. and Radwan, J. (2016) AMPLISAS: a web server for multilocus genotyping using next‐generation amplicon sequencing data. Mol Ecol Resour, 16, 498-510. doi: 10.1111/1755-0998.12453

Domesticated organisms are excellent models for understanding ecology and evolution and they are important for our food production and safety. While less studied than plants and animals, micro-organisms have also been domesticated, in particular for food fermentation [1]. The most studied domesticated micro-organism is the yeast used to make wine, beer and bread, Saccharomyces cerevisiae [2, 3, 4].

Filamentous fungi used for cheese-making have recently gained interest, for example Penicillium roqueforti used to make blue cheeses and P. camemberti to make soft cheeses [5, 6, 7, 8]. As for plants and animals, domestication has led to beneficial traits for food production in fermenting fungi, but also to bottlenecks and degeneration [6, 7, 9]; P. camemberti for example does not produce enough spores any more for optimal culture and inoculation and P. roqueforti has lost sexual fertility [9]. The loss of genetic diversity and of species diversity in our food production system is concerning for multiple reasons : i) it jeopardizes future improvement in the face of global changes ; ii) it causes the loss of evolved diversity during centuries under human selection, and therefore of beneficial characteristics and specificities that we may never be able to recover ; iii) it leads to degeneration in the few cultivated strains; iv) it impoverishes the diversity of our food products and local adaptation of production practices. 

The study of domesticated fungi used for food fermentation has focused so far on the evolution of lineages and on their metabolic specificities. Microbiological assemblages and species diversity have been much less studied, while they likely also have a strong impact on the quality and safety of final products. This study by Elisa Michel and colleagues [10] addresses this question, using an interdisciplinary participatory research approach including bakers, psycho-sociologists and microbiologists to analyse bread-making practices and their impact on microbial communities in sourdough. Elisa Michel and colleagues [10] identified two distinct groups of bread-making practices based on interviews and surveys, with farmer-like practices (low bread production, use of ancient wheat populations, manual kneading, working at ambient temperature, long fermentation periods and no use of commercial baker’s yeast) versus more intensive, artisanal-like practices. Metabarcoding and microbial culture-based analyses showed that the well-known baker’s yeast, Saccharomyces cerevisiae, was, surprisingly, not the most common species in French organic sourdoughs. Kazachstania was the most represented yeast genus over all sourdoughs, both in terms of read abundance and of species diversity. Kazachstania species were also often dominant in individual sourdoughs, but Saccharomyces uvarum or Torulaspora delbrueckii could also be the dominant yeast species.

Metabarcoding analyses further revealed that the composition of the fungal communities differed between the farmer-like and more intensive practices, representing the first evidence of the influence of artisanal practices on microbial communities. The fungal communities were impacted by a combination of bread-making variables including the type of wheat varieties, the length of fermentation, the quantity of bread made per week and the use of commercial yeast. Maintaining on farm less intensive bread-making practices, may allow the preservation of typical species and phenotypic diversity in microbial communities in sourdough. Farmer-like practices did not lead to higher diversity within sourdoughs but, overall, the diversity of bread-making practices allow maintaining a larger diversity in sourdoughs. For example, different Kazachstania species were most abundant in sourdoughs from artisanal-like and farmer-like practices. Interviews with the bakers suggested the role of dispersal of Kazachstania species in shaping sourdough microbial communities, dispersal occurring by seed exchanges, sourdough mixing or gifts, bread-making training in common or working in one another’s bakery. Nikolai Vavilov [11] had already highlighted for crops the importance of isolated cultures and selection in different farms for generating and preserving crop diversity, but also the importance of seed exchange for fostering adaptation. 

Furthermore, one of the yeast frequently found in artisanal sourdoughs, Kazachstania humilis, displayed phenotypic differences between sourdough and non-sourdough strains, suggesting domestication. The sourdough strains exhibited significantly higher CO2 production rate and a lower fermentation latency-phase time. 

The study by Elisa Michel and colleagues [10] is thus novel and inspiring in showing the importance of interdisciplinary studies, combining metabarcoding, microbiology and interviews for assessing the composition and diversity of microbial communities in human-made food, and in revealing the impact of artisanal-like bread-making practices in preserving microbial community diversity.

Interdisciplinary studies are still rare but have already shown the importance of combining ethno-ecology, biology and evolution to decipher the role of human practices on genetic diversity in crops, animals and food microorganisms and to help preserving genetic resources [12]. For example, in the case of the bread wheat Triticum aestivum, such interdisciplinary studies have shown that genetic diversity has been shaped by farmers’ seed diffusion and farming practices [13]. We need more of such interdisciplinary studies on the impact of farmer versus industrial agricultural and food-making practices as we urgently need to preserve the diversity of micro-organisms used in food production that we are losing at a rapid pace [6, 7, 14]. 

References

[1] Dupont J, Dequin S, Giraud T, Le Tacon F, Marsit S, Ropars J, Richard F, Selosse M-A (2017) Fungi as a Source of Food. Microbiology Spectrum, 5, 5.3.09. https://doi.org/10.1128/microbiolspec.FUNK-0030-2016

[2] Legras J-L, Galeote V, Bigey F, Camarasa C, Marsit S, Nidelet T, Sanchez I, Couloux A, Guy J, Franco-Duarte R, Marcet-Houben M, Gabaldon T, Schuller D, Sampaio JP, Dequin S (2018) Adaptation of S. cerevisiae to Fermented Food Environments Reveals Remarkable Genome Plasticity and the Footprints of Domestication. Molecular Biology and Evolution, 35, 1712–1727. https://doi.org/10.1093/molbev/msy066

[3] Bai F-Y, Han D-Y, Duan S-F, Wang Q-M (2022) The Ecology and Evolution of the Baker’s Yeast Saccharomyces cerevisiae. Genes, 13, 230. https://doi.org/10.3390/genes13020230

[4] Fay JC, Benavides JA (2005) Evidence for Domesticated and Wild Populations of Saccharomyces cerevisiae. PLOS Genetics, 1, e5. https://doi.org/10.1371/journal.pgen.0010005

[5] Ropars J, Rodríguez de la Vega RC, López-Villavicencio M, Gouzy J, Sallet E, Dumas É, Lacoste S, Debuchy R, Dupont J, Branca A, Giraud T (2015) Adaptive Horizontal Gene Transfers between Multiple Cheese-Associated Fungi. Current Biology, 25, 2562–2569. https://doi.org/10.1016/j.cub.2015.08.025

[6] Dumas E, Feurtey A, Rodríguez de la Vega RC, Le Prieur S, Snirc A, Coton M, Thierry A, Coton E, Le Piver M, Roueyre D, Ropars J, Branca A, Giraud T (2020) Independent domestication events in the blue-cheese fungus Penicillium roqueforti. Molecular Ecology, 29, 2639–2660. https://doi.org/10.1111/mec.15359

[7] Ropars J, Didiot E, Rodríguez de la Vega RC, Bennetot B, Coton M, Poirier E, Coton E, Snirc A, Le Prieur S, Giraud T (2020) Domestication of the Emblematic White Cheese-Making Fungus Penicillium camemberti and Its Diversification into Two Varieties. Current Biology, 30, 4441-4453.e4. https://doi.org/10.1016/j.cub.2020.08.082

[8] Caron T, Piver ML, Péron A-C, Lieben P, Lavigne R, Brunel S, Roueyre D, Place M, Bonnarme P, Giraud T, Branca A, Landaud S, Chassard C (2021) Strong effect of Penicillium roqueforti populations on volatile and metabolic compounds responsible for aromas, flavor and texture in blue cheeses. International Journal of Food Microbiology, 354, 109174. https://doi.org/10.1016/j.ijfoodmicro.2021.109174

[9] Ropars J, Lo Y-C, Dumas E, Snirc A, Begerow D, Rollnik T, Lacoste S, Dupont J, Giraud T, López-Villavicencio M (2016) Fertility depression among cheese-making Penicillium roqueforti strains suggests degeneration during domestication. Evolution, 70, 2099–2109. https://doi.org/10.1111/evo.13015

[10] Michel E, Masson E, Bubbendorf S, Lapicque L, Nidelet T, Segond D, Guézenec S, Marlin T, Devillers H, Rué O, Onno B, Legrand J, Sicard D, Bakers TP (2022) Artisanal and farmer bread making practices differently shape fungal species community composition in French sourdoughs. bioRxiv, 679472, ver. 6 peer-reviewed and recommended by Peer Community in Evolutionary Biology. https://doi.org/10.1101/679472

[11] Vavilov NI, Vavylov MI, Dorofeev VF (1992) Origin and Geography of Cultivated Plants. Cambridge University Press.

[12] Saslis-Lagoudakis CH, Clarke AC (2013) Ethnobiology: the missing link in ecology and evolution. Trends in Ecology & Evolution, 28, 67–68. https://doi.org/10.1016/j.tree.2012.10.017

[13] Thomas M, Demeulenaere E, Dawson JC, Khan AR, Galic N, Jouanne-Pin S, Remoue C, Bonneuil C, Goldringer I (2012) On-farm dynamic management of genetic diversity: the impact of seed diffusions and seed saving practices on a population-variety of bread wheat. Evolutionary Applications, 5, 779–795. https://doi.org/10.1111/j.1752-4571.2012.00257.x

[14] Demeulenaere É, Lagrola M (2021) Des indicateurs pour accompagner “ les éleveurs de microbes” : Une communauté épistémique face au problème des laits “ paucimicrobiens ” dans la production fromagère au lait cru (1995-2015). Revue d’anthropologie des connaissances, 15. http://journals.openedition.org/rac/24953

Nobody can ignore that a full understanding of evolution requires an integrated approach from both conceptual and methodological viewpoints. Although some life-history traits, e.g. flowering time, have long been receiving more attention than others, in many cases because the former are more workable than the latter, we must acknowledge that our comprehension about how evolution works is strongly biased and limited. In the Arabidopsis community, such an integration is making good progress as an increasing number of research groups worldwide are changing the way in which evolution is put to the test.

This manuscript [1] is a good example of that as the authors raise an important issue in evolutionary biology by combining gene expression and flowering time data from different sources. In particular, the authors explore how variation in flowering time, which determines lifespan, and host immunity defenses co-vary, which is interpreted in terms of co-evolution between the two traits. Interestingly, the authors go beyond that pattern by separating lifespan-dependent from lifespan–independent defense genes, and by showing that defense genes with variants known to impact fitness in the field are among the genes whose expression co-varies most strongly with flowering time. Finally, these results are supported by a simple mathematical model indicating that such a relationship can also be expected theoretically.

Overall, the readers will find many conceptual and methodological elements of interest in this manuscript. The idea that evolution is better understood under the scope of life history variation is really exciting and challenging, and in my opinion on the right track for disentangling the inherent complexities of evolutionary research. However, only when we face complexity, we also face its costs and burdens. In this particular case, the well-known co-variation between seed dormancy and flowering time is a missing piece, as well as the identification of (variation in) putative selective pressures accounting for the co-evolution between defense mechanisms and life history (seed dormancy vs. flowering time) along environmental gradients. More intellectual, technical and methodological challenges that with no doubt are totally worth it.

Reference

[1] Glander S, He F, Schmitz G, Witten A, Telschow A, de Meaux J. 2017. Assortment of flowering time and defense alleles in natural Arabidopsis thaliana populations suggests co-evolution between defense and vegetative lifespan strategies. bioRxiv ver.1 of June 19, 2017. doi: 10.1101/131136

The present-day distribution of extant species is the result of the interplay between their past population demography (e.g., expansion, contraction, isolation, and migration) and adaptation to the environment. Shedding light on the timing and magnitude of key demographic events helps identify potential drivers of such events and interaction of those drivers, such as life history traits and past episodes of environmental shifts.

The understanding of the key factors driving species evolution gives important insights into how the species may respond to changing conditions, which can be particularly relevant for the management of harmful species, such as agricultural pests (e.g. [1]). Meaningful demographic inferences present major challenges. These include formulating evolutionary scenarios fitting species biology and the eco-geographical context and choosing informative molecular markers and accurate quantitative approaches to statistically compare multiple demographic scenarios and estimate the parameters of interest. A further issue comes with result interpretation. Accurately dating the inferred events is far from straightforward since reliable calibration points are necessary to translate the molecular estimates of the evolutionary time into absolute time units (i.e. years). This can be attempted in different ways, such as by using fossil and archaeological records, heterochronous samples (e.g. ancient DNA), and/or mutation rate estimated from independent data (e.g. [2], [3] for review). Nonetheless, most experimental systems rarely meet these conditions, hindering the comprehensive interpretation of results.

The contribution of Chapuis et al. [4] addresses these issues to investigate the recent history of the African insect pest Schistocerca gregaria (desert locust). They apply Approximate Bayesian Computation-Random Forest (ABC-RF) approaches to microsatellite markers. Owing to their fast mutation rate microsatellite markers offer at least two advantages: i) suitability for analyzing recently diverged populations, and ii) direct estimate of the germline mutation rate in pedigree samples. The work of Chapuis et al. [4] benefits of both these advantages, since they have estimates of mutation rate and allele size constraints derived from germline mutations in the species [5].

The main aim of the study is to infer the history of divergence of the two subspecies of the desert locust, which have spatially disjoint distribution corresponding to the dry regions of North and West-South Africa. They first use paleo-vegetation maps to formulate hypotheses about changes in species range since the last glacial maximum. Based on them, they generate 12 divergence models. For the selection of the demographic model and parameter estimation, they apply the recently developed ABC-RF approach, a powerful inferential tool that allows optimizing the use of summary statistics information content, among other advantages [6]. Some methodological novelties are also introduced in this work, such as the computation of the error associated with the posterior parameter estimates under the best scenario. The accuracy of timing estimate is assured in two ways: i) by the use of microsatellite markers with known evolutionary dynamics, as underlined above, and ii) by assessing the divergence time threshold above which posterior estimates are likely to be biased by size homoplasy and limits in allele size range [7]. The best-supported model suggests a recent divergence event of the subspecies of S. gregaria (around 2.6 kya) and a reduction of populations size in one of the subspecies (S. g. flaviventris) that colonized the southern distribution area. As such, results did not support the hypothesis that the southward colonization was driven by the expansion of African dry environments associated with the last glacial maximum, as it has been postulated for other arid-adapted species with similar African disjoint distributions [8]. The estimated time of divergence points at a much more recent origin for the two subspecies, during the late Holocene, in a period corresponding to fairly stable arid conditions similar to current ones [9,10].

Although the authors cannot exclude that their microsatellite data bear limited information on older colonization events than the last one, they bring arguments in favour of alternative explanations. The hypothesis privileged does not involve climatic drivers, but the particularly efficient dispersal behaviour of the species, whose individuals are able to fly over long distances (up to thousands of kilometers) under favourable windy conditions. A single long-distance dispersal event by a few individuals would explain the genetic signature of the bottleneck. There is a growing number of studies in phylogeography in arid regions in the Southern hemisphere, but the impact of past climate changes on the species distribution in this region remains understudied relative to the Northern hemisphere [11,12].

The study presented by Chapuis et al. [4] offers several important insights into demographic changes and the evolutionary history of an agriculturally important pest species in Africa, which could also mirror the history of other organisms in the continent. As the authors point out, there are necessarily some uncertainties associated with the models of past ecosystems and climate, especially for Africa. Interestingly, the authors argue that the information on paleo-vegetation turnover was more informative than climatic niche modeling for the purpose of their study since it made them consider a wider range of bio-geographical changes and in turn a wider range of evolutionary scenarios (see discussion in Supplementary Material). Microsatellite markers have been offering a useful tool in population genetics and phylogeography for decades, but their popularity is perhaps being taken over by single nucleotide polymorphism (SNP) genotyping and whole-genome sequencing (WGS) (the peak year of the number of the publication with “microsatellite” is in 2012 according to PubMed).

This study reaffirms the usefulness of these classic molecular markers to estimate past demographic events, especially when species- and locus-specific microsatellite mutation features are available and a powerful inferential approach is adopted. Nonetheless, there are still hurdles to overcome, such as the limitations in scenario choice associated with the simulation software used (e.g. not allowing for continuous gene flow in this particular case), which calls for further improvement of simulation tools allowing for more flexible modeling of demographic events and mutation patterns. In sum, this work not only contributes to our understanding of the makeup of the African biodiversity but also offers a useful statistical framework, which can be applied to a wide array of species and molecular markers (microsatellites, SNPs, and WGS).

References

[1] Lehmann, P. et al. (2018). Complex responses of global insect pests to climate change. bioRxiv, 425488. doi: https://dx.doi.org/10.1101/425488

[2] Donoghue, P. C., & Benton, M. J. (2007). Rocks and clocks: calibrating the Tree of Life using fossils and molecules. Trends in Ecology & Evolution, 22(8), 424-431. doi: https://dx.doi.org/10.1016/j.tree.2007.05.005

[3] Ho, S. Y., Lanfear, R., Bromham, L., Phillips, M. J., Soubrier, J., Rodrigo, A. G., & Cooper, A. (2011). Time‐dependent rates of molecular evolution. Molecular ecology, 20(15), 3087-3101. doi: https://dx.doi.org/10.1111/j.1365-294X.2011.05178.x

[4] Chapuis, M.-P., Raynal, L., Plantamp, C., Meynard, C. N., Blondin, L., Marin, J.-M. and Estoup, A. (2020). A young age of subspecific divergence in the desert locust Schistocerca gregaria, inferred by ABC Random Forest. bioRxiv, 671867, ver. 4 peer-reviewed and recommended by PCI Evolutionary Biology. doi: https://dx.doi.org/10.1101/671867

5] Chapuis, M.-P., Plantamp, C., Streiff, R., Blondin, L., & Piou, C. (2015). Microsatellite evolutionary rate and pattern in Schistocerca gregaria inferred from direct observation of germline mutations. Molecular ecology, 24(24), 6107-6119. doi: https://dx.doi.org/10.1111/mec.13465

[6] Raynal, L., Marin, J. M., Pudlo, P., Ribatet, M., Robert, C. P., & Estoup, A. (2018). ABC random forests for Bayesian parameter inference. Bioinformatics, 35(10), 1720-1728. doi: https://dx.doi.org/10.1093/bioinformatics/bty867

[7] Estoup, A., Jarne, P., & Cornuet, J. M. (2002). Homoplasy and mutation model at microsatellite loci and their consequences for population genetics analysis. Molecular ecology, 11(9), 1591-1604. doi: https://dx.doi.org/10.1046/j.1365-294X.2002.01576.x

[8] Moodley, Y. et al. (2018). Contrasting evolutionary history, anthropogenic declines and genetic contact in the northern and southern white rhinoceros (Ceratotherium simum). Proceedings of the Royal Society B, 285(1890), 20181567. doi: https://dx.doi.org/10.1098/rspb.2018.1567

[9] Kröpelin, S. et al. (2008). Climate-driven ecosystem succession in the Sahara: the past 6000 years. science, 320(5877), 765-768. doi: https://dx.doi.org/10.1126/science.1154913

[10] Maley, J. et al. (2018). Late Holocene forest contraction and fragmentation in central Africa. Quaternary Research, 89(1), 43-59. doi: https://dx.doi.org/10.1017/qua.2017.97

[11] Beheregaray, L. B. (2008). Twenty years of phylogeography: the state of the field and the challenges for the Southern Hemisphere. Molecular Ecology, 17(17), 3754-3774. doi: https://dx.doi.org/10.1111/j.1365-294X.2008.03857.x

[12] Dubey, S., & Shine, R. (2012). Are reptile and amphibian species younger in the Northern Hemisphere than in the Southern Hemisphere?. Journal of evolutionary biology, 25(1), 220-226. doi: https://dx.doi.org/10.1111/j.1420-9101.2011.02417.x

*****

A video about this preprint is available here:

It has long been known that most multicellular eukaryotes rely on microbial partners for a variety of functions including nutrition, immune reactions and defence against enemies. Lichens are probably the most popular example of a symbiosis involving a photosynthetic microorganism (an algae, a cyanobacteria or both) living embedded within the filaments of a fungus (usually an ascomycete). The latter is the backbone structure of the lichen, whereas the former provides photosynthetic products. Lichens are unique among symbioses because the structures the fungus and the photosynthetic microorganism form together do not resemble any of the two species living in isolation. Classic textbook examples like lichens are not often challenged and this is what Toby Spribille and his co-authors did with their paper published in July 2016 in Science [1]. This story started with the study of two species of macrolichens from the class of Lecanoromycetes that are commonly found in the mountains of Montana (US): Bryoria fremontii and B. tortuosa. For more than 90 years, these species have been known to differ in their chemical composition and colour, but studies performed so far failed in finding differences at the molecular level in both the mycobiont and the photobiont. These two species were therefore considered as nomenclatural synonyms, and the origin of their differences remained elusive. To solve this mystery, the authors of this work performed a transcriptome-wide analysis that, relative to previous studies, expanded the taxonomic range to all Fungi. This analysis revealed higher abundances of a previously unknown basidiomycete yeast from the genus Cyphobasidium in one of the lichen species, a pattern that was further confirmed by combining microscopy imaging and the fluorescent in situ hybridisation technique (FISH).

Finding out that a previously unknown micro-organism changes the colour and the chemical composition of an organism is surprising but not new. For instance, bacterial symbionts are able to trigger colour changes in some insect species [2], and endophyte fungi are responsible for the production of defensive compounds in the leaves of several grasses [3]. The study by Spribille and his co-authors is fascinating because it demonstrates that Cyphobasidium yeasts have played a key role in the evolution and diversification of Lecanoromycetes, one of the most diverse classes of macrolichens. Indeed these basidiomycete yeasts were not only found in Bryoria but in 52 other lichen genera from all six continents, and these included 42 out of 56 genera in the family Parmeliaceae. Most of these sequences formed a highly supported monophyletic group, and a molecular clock revealed that the origin of many macrolichen groups occurred around the same time Cyphobasidium yeasts split from Cystobasidium, their nearest relatives. This newly discovered passenger is therefore an ancient inhabitant of lichens and has driven the evolution of this emblematic group of organisms.

This study raises an important question on the stability of complex symbiotic partnerships. In intimate obligatory symbioses the evolutionary interests of both partners are often identical and what is good for one is also good for the other. This is the case of several insects that feed on poor diets like phloem and xylem sap, and which carry vertically-transmitted symbionts that provide essential nutrients. Molecular phylogenetic studies have repeatedly shown that in several insect groups transition to phloem or xylem feeding occurred at the same time these nutritional symbionts were acquired [4]. In lichens, an outstanding question is to know what was the key feature Cyphobasidium yeasts brought to the symbiosis. As suggested by the authors, these yeasts are likely to be involved in the production of secondary defensive metabolites and architectural structures, but, are these services enough to explain the diversity found in macrolichens? This paper is an appealing example of a multipartite symbiosis where the different partners share an ancient evolutionary history.

References

[1] Spribille T, Tuovinen V, Resl P, et al. 2016. Basidiomycete yeasts in the cortex of ascomycete macrolichens. Science 353:488–92. doi: 10.1126/science.aaf8287

[2] Tsuchida T, Koga R, Horikawa M, et al. 2010. Symbiotic Bacterium Modifies Aphid Body Color. Science 330:1102–1104. doi: 10.1126/science.1195463

[3] Clay K. 1988. Fungal Endophytes of Grasses: A Defensive Mutualism between Plants and Fungi. Ecology 69:10–16. doi: 10.2307/1943155

[4] Moran NA. 2007. Symbiosis as an adaptive process and source of phenotypic complexity. Proceeding of the National Academy of Science USA 104:8627–8633. doi: 10.1073/pnas.0611659104