REFMAKER: Make your own reference to target nuclear loci in low coverage genome skimming libraries. Phylogenomic application in Sapotaceae.

Genome skimming approach is widely used in plant systematics to infer phylogenies mostly from organelle genomes. However, organelles represent only 10 % of the produced libraries, and the low coverage associated with these libraries (<3X) prevents the capture of nuclear sequences, which are not always available in non-model organisms or limited to the ribosomal regions. We developed REFMAKER, a user-friendly pipeline, to create specific sets of nuclear loci that can be extracted directly from the genome skimming libraries. For this, a catalogue is built from the meta-assembly of each library contigs, and cleaned by selecting the nuclear regions and removing duplicates from clustering steps. Libraries are next mapped onto this catalogue and consensus sequences are generated to produce a ready-to-use phylogenetic matrix following different filtering parameters aiming at removing putative errors and paralogous sequences. REFMAKER allowed us to infer a well resolved phylogeny in Capurodendron (Sapotaceae) on 67 nuclear loci from low-coverage libraries (<1X). The resulting phylogeny is concomitant with one previously inferred on 638 nuclear genes from target enrichment libraries. While it remains preliminary because of this low sequencing depth, REFMAKER therefore opens perspectives in phylogenomics by allowing nuclear phylogeny reconstructions with genome skimming datasets.

Species Delimitation and Conservation in Taxonomically Challenging Lineages: The Case of Two Clades of Capurodendron (Sapotaceae) in Madagascar.

Capurodendron is the largest endemic genus of plants from Madagascar, with around 76% of its species threatened by deforestation and illegal logging. However, some species are not well circumscribed and many of them remain undescribed, impeding a confident evaluation of their conservation status. Here we focus on taxa delimitation and conservation of two species complexes within Capurodendron: the Arid and Western complexes, each containing undescribed morphologies as well as intermediate specimens alongside well-delimited taxa. To solve these taxonomic issues, we studied 381 specimens morphologically and selected 85 of them to obtain intergenic, intronic, and exonic protein-coding sequences of 794 nuclear genes and 227 microsatellite loci. These data were used to test species limits and putative hybrid patterns using different approaches such as phylogenies, PCA, structure analyses, heterozygosity level, FST, and ABBA-BABA tests. The potential distributions were furthermore estimated for each inferred species. The results show that the Capurodendron Western Complex contains three well-delimited species, C. oblongifolium, C. perrieri, and C. pervillei, the first two hybridizing sporadically with the last and producing morphologies similar to, but genetically distinct from C. pervillei. The Arid Complex shows a more intricate situation, as it contains three species morphologically well-delimited but genetically intermixed. Capurodendron mikeorum nom. prov. is shown to be an undescribed species with a restricted distribution, while C. androyense and C. mandrarense have wider and mostly sympatric distributions. Each of the latter two species contains two major genetic pools, one showing interspecific admixture in areas where both taxa coexist, and the other being less admixed and comprising allopatric populations having fewer contacts with the other species. Only two specimens out of 172 showed clear genetic and morphological signals of recent hybridization, while all the others were morphologically well-delimited, independent of their degree of genetic admixture. Hybridization between Capurodendron androyense and C. microphyllum, the sister species of the Arid Complex, was additionally detected in areas where both species coexist, producing intermediate morphologies. Among the two complexes, species are well-defined morphologically with the exception of seven specimens (1.8%) displaying intermediate patterns and genetic signals compatible with a F1 hybridization. A provisional conservation assessment for each species is provided.

Phylogeographic reconstructions can be biased by ancestral shared alleles: The case of the polymorphic lichen Bryoria fuscescens in Europe and North Africa.

Large phylogeographic studies on lichens are scarce, and none involves a single species within which different lineages show fixed alternative dispersal strategies. We investigated Bryoria fuscescens (including B. capillaris) in Europe and western North Africa by phenotypically characterizing 1400 specimens from 64 populations and genotyping them with 14 microsatellites. We studied population structure and genetic diversity at the local and continental scales, discussed the post-glacial phylogeography, and compared dispersal capacities of phenotypes with and without soralia. Our main hypothesis is that the estimated phylogeography, migration routes, and dispersal capacities may be strongly biased by ancestral shared alleles. Scandinavia is genetically the richest area, followed by the Iberian Peninsula, the Carpathians, and the Alps. Three gene pools were detected: two partially linked to phenotypic characteristics, and the third one genetically related to the American sister species B. pseudofuscescens. The comparison of one gene pool producing soredia and one not, suggested both as panmictic, with similar levels of isolation by distance (IBD). The migration routes were estimated to span from north to south, in disagreement with the assessed glacial refugia. The presence of ancestral shared alleles in distant populations can explain the similar IBD levels found in both gene pools while producing a false signal of panmixia, and also biasing the phylogeographic reconstruction. The incomplete lineage sorting recorded for DNA sequence loci also supports this hypothesis. Consequently, the high diversity in Scandinavia may rather come from recent immigration into northern populations than from an in situ diversification. Similar patterns of ancestral shared polymorphism may bias the phylogeographical reconstruction of other lichen species.

New genetic markers for Sapotaceae phylogenomics: More than 600 nuclear genes applicable from family to population levels.

Some tropical plant families, such as the Sapotaceae, have a complex taxonomy, which can be resolved using Next Generation Sequencing (NGS). For most groups however, methodological protocols are still missing. Here we identified 531 monocopy genes and 227 Short Tandem Repeats (STR) markers and tested them on Sapotaceae using target capture and NGS. The probes were designed using two genome skimming samples from Capurodendron delphinense and Bemangidia lowryi, both from the Tseboneae tribe, as well as the published Manilkara zapota transcriptome from the Sapotoideae tribe. We combined our probes with 261 additional ones previously published and designed for the entire angiosperm group. On a total of 792 low-copy genes, 638 showed no signs of paralogy and were used to build a phylogeny of the family with 231 individuals from all main lineages. A highly supported topology was obtained at high taxonomic ranks but also at the species level. This phylogeny revealed the existence of more than 20 putative new species. Single nucleotide polymorphisms (SNPs) extracted from the 638 genes were able to distinguish lineages within a species complex and to highlight geographical structuration. STR were recovered efficiently for the species used as reference (C. delphinense) but the recovery rate decreased dramatically with the phylogenetic distance to the focal species. Altogether, the new loci will help reaching a sound taxonomic understanding of the family Sapotaceae for which many circumscriptions and relationships are still debated, at the species, genus and tribe levels.

Evolutionary Trends in the Mitochondrial Genome of Archaeplastida: How Does the GC Bias Affect the Transition from Water to Land?

Among the most intriguing mysteries in the evolutionary biology of photosynthetic organisms are the genesis and consequences of the dramatic increase in the mitochondrial and nuclear genome sizes, together with the concomitant evolution of the three genetic compartments, particularly during the transition from water to land. To clarify the evolutionary trends in the mitochondrial genome of Archaeplastida, we analyzed the sequences from 37 complete genomes. Therefore, we utilized mitochondrial, plastidial and nuclear ribosomal DNA molecular markers on 100 species of Streptophyta for each subunit. Hierarchical models of sequence evolution were fitted to test the heterogeneity in the base composition. The best resulting phylogenies were used for reconstructing the ancestral Guanine-Cytosine (GC) content and equilibrium GC frequency (GC*) using non-homogeneous and non-stationary models fitted with a maximum likelihood approach. The mitochondrial genome length was strongly related to repetitive sequences across Archaeplastida evolution; however, the length seemed not to be linked to the other studied variables, as different lineages showed diverse evolutionary patterns. In contrast, Streptophyta exhibited a powerful positive relationship between the GC content, non-coding DNA, and repetitive sequences, while the evolution of Chlorophyta reflected a strong positive linear relationship between the genome length and the number of genes.

Evolution of complex symbiotic relationships in a morphologically derived family of lichen-forming fungi.

We studied the evolutionary history of the Parmeliaceae (Lecanoromycetes, Ascomycota), one of the largest families of lichen-forming fungi with complex and variable morphologies, also including several lichenicolous fungi. We assembled a six-locus data set including nuclear, mitochondrial and low-copy protein-coding genes from 293 operational taxonomic units (OTUs). The lichenicolous lifestyle originated independently three times in lichenized ancestors within Parmeliaceae, and a new generic name is introduced for one of these fungi. In all cases, the independent origins occurred c. 24 million yr ago. Further, we show that the Paleocene, Eocene and Oligocene were key periods when diversification of major lineages within Parmeliaceae occurred, with subsequent radiations occurring primarily during the Oligocene and Miocene. Our phylogenetic hypothesis supports the independent origin of lichenicolous fungi associated with climatic shifts at the Oligocene-Miocene boundary. Moreover, diversification bursts at different times may be crucial factors driving the diversification of Parmeliaceae. Additionally, our study provides novel insight into evolutionary relationships in this large and diverse family of lichen-forming ascomycetes.

Characterization of microsatellite loci in lichen-forming fungi of Bryoria section Implexae (Parmeliaceae).

PREMISE OF THE STUDY: The locally rare, haploid, lichen-forming fungi Bryoria capillaris, B. fuscescens, and B. implexa are associated with boreal forests and belong to Bryoria sect. Implexae. Recent phylogenetic studies consider them to be conspecific. Microsatellite loci were developed to study population structure in Bryoria sect. Implexae and its response to ecosystem disturbances. • METHODS AND RESULTS: We developed 18 polymorphic microsatellite markers using 454 pyrosequencing data assessed in 82 individuals. The number of alleles per locus ranged from two to 13 with an average of 4.6. Nei's unbiased gene diversity, averaged over loci, ranged from 0.38 to 0.52. The markers amplified with all three species, except for markers Bi05, Bi15, and Bi18. • CONCLUSIONS: The new markers will allow the study of population subdivision, levels of gene introgression, and levels of clonal spread of Bryoria sect. Implexae. They will also facilitate an understanding of the effects of forest disturbance on genetic diversity of these lichen species.