Insights into the genome of the 'Loco' Concholepas concholepas (Gastropoda: Muricidae) from low-coverage short-read sequencing: genome size, ploidy, transposable elements, nuclear RNA gene operon, mitochondrial genome, and phylogenetic placement in the family Muricidae.

BACKGROUND: The Peruvian 'chanque' or Chilean 'loco' Concholepas concholepas is an economically, ecologically, and culturally important muricid gastropod heavily exploited by artisanal fisheries in the temperate southeastern Pacific Ocean. In this study, we have profited from a set of bioinformatics tools to recover important biological information of C. concholepas from low-coverage short-read NGS datasets. Specifically, we calculated the size of the nuclear genome, ploidy, and estimated transposable elements content using an in silico k-mer approach, we discovered, annotated, and quantified those transposable elements, we assembled and annotated the 45S rDNA RNA operon and mitochondrial genome, and we confirmed the phylogenetic position of C. concholepas within the muricid subfamily Rapaninae based on translated protein coding genes. RESULTS: Using a k-mer approach, the haploid genome size estimated for the predicted diploid genome of C. concholepas varied between 1.83 Gbp (with kmer = 24) and 2.32 Gbp (with kmer = 36). Between half and two thirds of the nuclear genome of C. concholepas was composed of transposable elements. The most common transposable elements were classified as Long Interspersed Nuclear Elements and Short Interspersed Nuclear Elements, which were more abundant than DNA transposons, simple repeats, and Long Terminal Repeats. Less abundant repeat elements included Helitron mobile elements, 45S rRNA DNA, and Satellite DNA, among a few others.The 45S rRNA DNA operon of C. concholepas that encodes for the ssrRNA, 5.8S rRNA, and lsrRNA genes was assembled into a single contig 8,090 bp long. The assembled mitochondrial genome of C. concholepas is 15,449 bp long and encodes 13 protein coding genes, two ribosomal genes, and 22 transfer RNAs. CONCLUSION: The information gained by this study will inform the assembly of a high quality nuclear genome for C. concholepas and will support bioprospecting and biomonitoring using environmental DNA to advance development of conservation and management plans in this overexploited marine snail.

Phylogenomic analyses sheds new light on the phylogeny and diversification of Corydalis DC. in Himalaya-Hengduan Mountains and adjacent regions.

The Himalaya-Hengduan Mountains (HHM), a renowned biodiversity hotspot of the world, harbors the most extensive habitats for alpine plants with extraordinary high levels of endemism. Although the general evolution pattern has been elucidated, the underlying processes driving spectacular radiations in many species-rich groups remain elusive. Corydalis DC. is widely distributed throughout the Northern Hemisphere containing more than 500 species, with high diversity in HHM and adjacent regions. Using 95 plastid genes, 3,258,640 nuclear single nucleotide polymorphisms (SNPs) and eight single-copy nuclear genes (SCNs) generated from genome skimming data, we reconstructed a robust time-calibrated phylogeny of Corydalis comprising more than 100 species that represented all subgenera and most sections. Molecular dating indicated that all main clades of Corydalis began to diverge in the Eocene, with the majority of extant species in HHM emerged from a diversification burst after the middle Miocene. Global pattern of mean divergence times indicated that species distributed in HHM were considerably younger than those in other regions, particularly for the two most species-rich clades (V and VI) of Corydalis. The early divergence and the recent diversification of Corydalis were most likely promoted by the continuous orogenesis and climate change associated with the uplift of the Qinghai-Tibetan Plateau (QTP). Our study demonstrates the effectivity of phylogenomic analyses with genome skimming data on the phylogeny of species-rich taxa, and sheds lights on how the uplift of QTP has triggered the evolutionary radiations of large plant genera in HHM and adjacent regions.

Imperiled wanderlust lichens in steppe habitats of western North America comprise geographically structured mycobiont lineages and a reversal to sexual reproduction within this asexual clade.

The northern North American Cordillera is a globally significant center of endemism. In western North America, imperiled arid steppe habitats support a number of unique species, including several endemic lichens. However, processes driving diversification and endemism in this region remain unclear. In this study, we investigate diversity and phylogeography of the threatened wanderlust lichens (mycobiont = Rhizoplaca species) which occur unattached on calcareous soils in steppe habitats. Wanderlust lichens comprise three species of lichen-forming fungi (LFF) - Rhizoplaca arbuscula, R. haydenii, and R. idahoensis (endangered, IUCN Red List) - which occur in fragmented populations in Idaho and Wyoming, with more limited populations in southern Montana and northern Utah. These lichens reproduce almost exclusively via large, asexual vegetative propagules. Here, our aims were to (i) assess the evolutionary origin of this group and identify phylogeographic structure, (ii) infer ancestral geographic distributions for lineages within this clade, and (iii) use species distribution modeling to better understand the distribution of contemporary populations. Using a genome-skimming approach, we generated a 19.1Mb alignment, spanning ca. half of the complete LFF genome, from specimens collected throughout the entire range of wanderlust lichens. Based on this phylogeny, we investigated phylogeographic patterns using RASP. Finally, we used MaxEnt to estimate species distribution models for R. arbuscula and R. haydenii. We inferred a highly structured topology, with clades corresponding to distinct geographic regions and morphologies collected throughout the group's distribution. We found that R. robusta, a sexually reproducing taxon, is clearly nested within the vagrant Rhizoplaca clade. Phylogeographic analyses suggest that both dispersal and vicariance played significant roles throughout the evolutionary history of the vagrant Rhizoplaca clade, with most of the dispersal events originating from the Salmon Basin in eastern Idaho - the center of diversity for this group. Despite the fact that wanderlust lichens are dispersal limited due to large, unspecialized vegetative propagules, we inferred multiple dispersal events crossing the Continental Divide. Comparing herbarium records with species distribution models suggests that wanderlust lichens don't fully occupy the areas of highest distribution probability. In fact, documented records often occur in areas predicted to be only marginally suitable. These data suggest a potential mismatch between contemporary habitats outside of the center of diversity in eastern Idaho with the most suitable habitat, adding to the vulnerability of this imperiled complex of endemic lichens.

Population genomics of the Isoetes appalachiana (Isoetaceae) complex supports a 'diploids-first' approach to conservation.

BACKGROUND AND AIMS: Allopolyploidy is an important driver of diversification and a key contributor to genetic novelty across the tree of life. However, many studies have questioned the importance of extant polyploid lineages, suggesting that the vast majority may constitute evolutionary 'dead ends'. This has important implications for conservation efforts where polyploids and diploid progenitors often compete for wildlife management resources. Isoetes appalachiana is an allotetraploid that is broadly distributed throughout the eastern USA alongside its diploid progenitors, I. valida and I. engelmannii. As such, this species complex provides an excellent opportunity to investigate the processes that underpin the formation and survival of allopolyploid lineages. METHODS: Here we utilized RADseq and whole-chloroplast sequencing to unravel the demographic and evolutionary history of hybridization in this widespread species complex. We developed a modified protocol for phasing RADseq loci from an allopolyploid in order to examine each progenitor's genetic contribution independently in a phylogenetic context. Additionally, we conducted population-level analyses to examine genetic diversity and evidence of gene flow within species. KEY RESULTS: Isoetes appalachiana is the product of multiple phylogenetic origins, suggesting that formation and establishment of allopolyploids are common in this group. Hybridization appears to be unidirectional, with I. engelmannii consistently being the maternal progenitor. Additionally, we find that polyploid lineages are genetically isolated, rarely if ever experiencing gene flow between geographically distinct populations. CONCLUSIONS: Allopolyploid lineages of I. appalachiana appear to form frequently and experience a high degree of genetic isolation following formation. Thus, our results appear to corroborate the hypothesis that the vast majority of recently formed polyploids may represent evolutionary dead ends. However, this does not necessarily lessen the evolutionary importance or ecological impact of polyploidy per se. Accordingly, we propose a conservation strategy that prioritizes diploid taxa, thus preserving downstream processes that recurrently generate allopolyploid diversity.

Plastome phylogenomics and new fossil evidence from Dominican amber shed light on the evolutionary history of the Neotropical fern genus Pecluma.

PREMISE: Molecular studies based on chloroplast markers have questioned the monophyly of the fern genus Pecluma (Polypodioideae, Polypodiaceae), which has several species of Polypodium nested within it. We explored the delimitation of Pecluma and its biogeographic pattern by evaluating the phylogenetic position of four Polypodium species not sequenced thus far and integrating the first fossil evidence of Pecluma. METHODS: Using herbarium material, we applied a genome-skimming approach to obtain a phylogenetic hypothesis of Polypodioideae; assessed the combination of character states observed in the fossil from Miocene Dominican amber using a previously published phylogeny of Polypodioideae based on four plastid markers as framework; calculated divergence times; and conducted an ancestral area estimation. RESULTS: Within Polypodioideae, Pecluma was recovered as sister to Phlebodium. Three of the newly sequenced species-Polypodium otites, P. pinnatissimum, and P. ursipes-were recovered with maximum support within the Pecluma clade, whereas P. christensenii remained within Polypodium. The closest combination of character states of the fossil was found within Pecluma. Our biogeographic analyses suggest an Eocene origin of the genus in South America, with several subsequent Oligocene and Miocene colonization events to Mexico-Central America and to the West Indies. CONCLUSIONS: Although the circumscription of Pecluma is still challenging, our results elucidate the origin and age of the genus. The newly described fossil, Pecluma hispaniolae sp. nov., supports the hypothesis that the epiphytic communities of the Greater Antilles exhibit a constant generic composition since the Miocene. We propose new combinations (Pecluma otites, Pecluma pinnatissima, and Pecluma ursipes) to accommodate three species previously classified in Polypodium.

Complete Mitochondrial Genome and Phylogenetic Position of Nurudea zhengii Ren (Insecta, Hemiptera, Aphididae).

Nurudea zhengii Ren was identified by aphid morphological characteristics as well as the gall shape and host plant species, and placed in the tribe Fordini (Hemiptera, Aphididae, Eriosomatinae). Here, its whole genome was firstly sequenced by a genome-skimming method and its mitochondrial genome (mitogenome) was assembled to examine its genetic variation and phylogenetic position. The complete mitogenome of Nurudea zhengii is 15,392 bp in length, and consists of 13 protein-coding genes, 22 tRNAs, two rRNAs and one D-loop region. The gene order follows the mitogenomes of the other Rhus gall aphids, and similarly has an AT bias with the content of 83.9%. The majority strand is A-skewed and C-skewed, and shows opposite skewness for G-skewed in the minority strands. The ratios of nonsynonymous to synonymous substitution rates of protein-coding genes are lower than one except for ATP8, which indicated that ATP8 was undergoing positive selection. Phylogenetic analysis among the Rhus gall aphids based on 13 protein-coding genes and two rRNA genes showed that N. zhengii was sister to N. shiraii, and then clustered with N. yanoniella as a group with high support value. The two species, N. shiraii and N. yanoniella, share the same host plant Rhus chinensis, while the host of N. zhengii is R. hypoleuca. However, the phylogenetic relationship indicated that the taxa sharing the same host plant were not absolutely clustered as the closest taxa at least at species level.

Exploring Taxonomic and Genetic Relationships in the Pinus mugo Complex Using Genome Skimming Data.

Genome skimming is a novel approach that enables obtaining large-scale genomic information based on high-copy DNA fractions from shallow whole-genome sequencing. The simplicity of this method, low analysis costs, and large amounts of generated data have made it widely used in plant research, including species identification, especially in the case of protected or endangered taxa. This task is particularly difficult in the case of closely related taxa. The Pinus mugo complex includes several dozen closely related taxa occurring in the most important mountain ranges in Europe. The taxonomic rank, origin, or distribution of many of these taxa have been debated for years. In this study, we used genome skimming and multilocus DNA barcoding approaches to obtain different sequence data sets and also to determine their genetic diversity and suitability for distinguishing closely related taxa in the Pinus mugo complex. We generated seven different data sets, which were then analyzed using three discrimination methods, i.e., tree based, distance based, and assembling species by automatic partitioning. Genetic diversity among populations and taxa was also investigated using haplotype network analysis and principal coordinate analysis. The proposed data set based on divergence hotspots is even twenty-times more variable than the other analyzed sets and improves the phylogenetic resolution of the Pinus mugo complex. In light of the obtained results, Pinus × rhaetica does not belong to the Pinus mugo complex and should not be identified with either Pinus uliginosa or Pinus rotundata. It seems to represent a fixed hybrid or introgressant between Pinus sylvestris and Pinus mugo. In turn, Pinus mugo and Pinus uncinata apparently played an important role in the origins of Pinus uliginosa and Pinus rotundata.

Using low-coverage whole genome sequencing (genome skimming) to delineate three introgressed species of buffalofish (Ictiobus).

Consumption of buffalofish has been sporadically associated with Haff disease-like illnesses involving sudden onset muscle pain and weakness due to skeletal muscle rhabdomyolysis, but determination of precisely which species are associated with these illnesses has been impeded by a lack of species-specific DNA-based markers. Here, three closely related species of buffalofish native to the Mississippi River Basin (Ictiobus bubalus, Ictiobus cyprinellus and Ictiobus niger) that have previously proven genetically indistinguishable using both mitochondrial and nuclear single-locus sequencing were reliably discriminated using low-coverage whole genome sequencing ('genome skimming'). Using 44 specimens representing the three species collected from the mid/upper (Missouri) and lower (Louisiana) regions of the species' native ranges, the SISRS (Site Identification from Short Read Sequences) bioinformatics pipeline was adapted to (1) identify over 620Mbp of putatively homologous nuclear sequence data and (2) isolate over 140,000 single-nucleotide polymorphisms (SNPs) that supported accurate species delimitation, all without the use of a reference genome or annotation data. These sites were used to classify Ictiobus spp. samples with genome-skim data, along with a larger set (n = 67) where ultraconserved elements (UCEs) were sequenced. Analyses of whole mitochondrial data revealed more limited signal. Nearly all samples matched their purported species based on morphologic identification, but two Missouri samples morphologically identified as I. niger grouped with samples of I. bubalus, albeit with significant enrichment of I. niger SNPs. To our knowledge this is the first report of a DNA-based tool to reliably discriminate these three morphologically distinct species.

Phylogenomics and morphological evolution of the mega-diverse genus Artemisia (Asteraceae: Anthemideae): implications for its circumscription and infrageneric taxonomy.

BACKGROUND AND AIMS: Artemisia is a mega-diverse genus consisting of ~400 species. Despite its medicinal importance and ecological significance, a well-resolved phylogeny for global Artemisia, a natural generic delimitation and infrageneric taxonomy remain missing, owing to the obstructions from limited taxon sampling and insufficient information on DNA markers. Its morphological characters, such as capitulum, life form and leaf, show marked variations and are widely used in its infrageneric taxonomy. However, their evolution within Artemisia is poorly understood. Here, we aimed to reconstruct a well-resolved phylogeny for global Artemisia via a phylogenomic approach, to infer the evolutionary patterns of its key morphological characters and to update its circumscription and infrageneric taxonomy. METHODS: We sampled 228 species (258 samples) of Artemisia and its allies from both fresh and herbarium collections, covering all the subgenera and its main geographical areas, and conducted a phylogenomic analysis based on nuclear single nucleotide polymorphisms (SNPs) obtained from genome skimming data. Based on the phylogenetic framework, we inferred the possible evolutionary patterns of six key morphological characters widely used in its previous taxonomy. KEY RESULTS: The genus Kaschgaria was revealed to be nested in Artemisia with strong support. A well-resolved phylogeny of Artemisia consisting of eight highly supported clades was recovered, two of which were identified for the first time. Most of the previously recognized subgenera were not supported as monophyletic. Evolutionary inferences based on the six morphological characters showed that different states of these characters originated independently more than once. CONCLUSIONS: The circumscription of Artemisia is enlarged to include the genus Kaschgaria. The morphological characters traditionally used for the infrageneric taxonomy of Artemisia do not match the new phylogenetic tree. They experienced a more complex evolutionary history than previously thought. We propose a revised infrageneric taxonomy of the newly circumscribed Artemisia, with eight recognized subgenera to accommodate the new results.

Phylogenetic Relationships and Next-Generation Barcodes in the Genus Torreya Reveal a High Proportion of Misidentified Cultivated Plants.

Accurate species identification is key to conservation and phylogenetic inference. Living plant collections from botanical gardens/arboretum are important resources for the purpose of scientific research, but the proportion of cultivated plant misidentification are un-tested using DNA barcodes. Here, we assembled the next-generation barcode (complete plastid genome and complete nrDNA cistron) and mitochondrial genes from genome skimming data of Torreya species with multiple accessions for each species to test the species discrimination and the misidentification proportion of cultivated plants used in Torreya studies. A total of 38 accessions were included for analyses, representing all nine recognized species of genus Torreya. The plastid phylogeny showed that all 21 wild samples formed species-specific clades, except T. jiulongshanensis. Disregarding this putative hybrid, seven recognized species sampled here were successfully discriminated by the plastid genome. Only the T. nucifera accessions grouped into two grades. The species identification rate of the nrDNA cistron was 62.5%. The Skmer analysis based on nuclear reads from genome skims showed promise for species identification with seven species discriminated. The proportion of misidentified cultivated plants from arboreta/botanical gardens was relatively high with four accessions (23.5%) representing three species. Interspecific relationships within Torreya were fully resolved with maximum support by plastomes, where Torreya jackii was on the earliest diverging branch, though sister to T. grandis in the nrDNA cistron tree, suggesting that this is likely a hybrid species between T. grandis and an extinct Torreya ancestor lineage. The findings here provide quantitative insights into the usage of cultivated samples for phylogenetic study.

Phylogeny and Historical Biogeography of the East Asian Clematis Group, Sect. Tubulosae, Inferred from Phylogenomic Data.

The evolutionary history of Clematis section Tubulosae, an East Asian endemic lineage, has not been comprehensively studied. In this study, we reconstruct the phylogeny of this section with a complete sampling using a phylogenomic approach. The genome skimming method was applied to obtain the complete plastome sequence, the nuclear ribosomal DNA (nrDNA), and the nuclear SNPs data for phylogenetic reconstruction. Using a Bayesian molecular clock approach and ancestral range reconstruction, we reconstruct biogeographical history and discuss the biotic and abiotic factors that may have shaped the distribution patterns of the section. Both nuclear datasets better resolved the phylogeny of the sect. Tubulosae than the plastome sequence. Sect. Tubulosae was resolved as a monophyletic group sister to a clade mainly containing species from the sect. Clematis and sect. Aspidanthera. Within sect. Tubulosae, two major clades were resolved by both nuclear datasets. Two continental taxa, C. heracleifolia and C. tubulosa var. ichangensis, formed one clade. One continental taxon, C. tubulosa, and all the other species from Taiwan island, the Korean peninsula, and the Japanese archipelago formed the other clade. Molecular dating results showed that sect. Tubulosae diverged from its sister clade in the Pliocene, and all the current species diversified during the Pleistocene. Our biogeographical reconstruction suggested that sect. Tubulosae evolved and began species diversification, most likely in mainland China, then dispersed to the Korean peninsula, and then expanded its range through the Japanese archipelago to Taiwan island. Island species diversity may arise through allopatric speciation by vicariance events following the range fragmentation triggered by the climatic oscillation and sea level change during the Pleistocene epoch. Our results highlight the importance of climatic oscillation during the Pleistocene to the spatial-temporal diversification patterns of the sect. Tubulosae.

An updated classification for the hyper-diverse genus Corydalis (Papaveraceae: Fumarioideae) based on phylogenomic and morphological evidence.

The genus Corydalis, with ca. 530 species, has long been considered taxonomically challenging because of its great variability. Previous molecular analyses, based on a few molecular markers and incomplete taxonomic sampling, were clearly inadequate to delimit sections and subgenera. We have performed phylogenetic analyses of Corydalis and related taxa, using 65 shared protein-coding plastid genes from 313 accessions (including 280 samples of ca. 226 species of Corydalis) and 152 universal low-copy nuclear genes from 296 accessions (including 271 samples of Corydalis) covering all 42 previously recognized sections and five independent "series". Phylogenetic trees were inferred using Bayesian Inference and Maximum Likelihood. Eight selected morphological characters were estimated using ancestral state reconstructions. Results include: (i) of the three subgenera of Corydalis, two are fully supported by both the plastid and nuclear data; the third, subg. Cremnocapnos, is weakly supported by plastid DNA only, whereas in the nuclear data the two included sections form successive outgroups to the rest of the genus; (ii) among all 42 sections and five "series", 25 sections and one "series" are resolved as monophyletic in both data sets; (iii) the common ancestor of Corydalis is likely to be a perennial plant with a taproot, yellow flowers with a short saccate spur, linear fruits with recurved fruiting pedicels, and seeds with elaiosomes; (iv) we provide a new classification of Corydalis with four subgenera (of which subg. Bipapillatae is here newly described), 39 sections, 16 of which are consistent with the previous classification, 16 sections have been recircumscribed, one section has been reinstated and six new sections are established. Characters associated with lifespan, underground structures, floral spur, fruit and elaiosomes are important for the recognition of subgenera and sections. These new phylogenetic analyses combined with ancestral character reconstructions uncovered previously unrecognized relationships, and greatly improved our understanding of the evolution of the genus.

Comparative analyses of three complete Primula mitogenomes with insights into mitogenome size variation in Ericales.

BACKGROUND: Although knowledge of the sizes, contents, and forms of plant mitochondrial genomes (mitogenomes) is increasing, little is known about the mechanisms underlying their structural diversity. Evolutionary information on the mitogenomes of Primula, an important ornamental taxon, is more limited than the information on their nuclear and plastid counterparts, which has hindered the comprehensive understanding of Primula mitogenomic diversity and evolution. The present study reported and compared three Primula mitogenomes and discussed the size expansion of mitogenomes in Ericales. RESULTS: Mitogenome master circles were sequenced and successfully assembled for three Primula taxa and were compared with publicly available Ericales mitogenomes. The three mitogenomes contained similar gene contents and varied primarily in their structures. The Primula mitogenomes possessed relatively high nucleotide diversity among all examined plant lineages. In addition, high nucleotide diversity was found among Primula species between the Mediterranean and Himalaya-Hengduan Mountains. Most predicted RNA editing sites appeared in the second amino acid codon, increasing the hydrophobic character of the protein. An early stop in atp6 caused by RNA editing was conserved across all examined Ericales species. The interfamilial relationships within Ericales and interspecific relationships within Primula could be well resolved based on mitochondrial data. Transfer of the two longest mitochondrial plastid sequences (MTPTs) occurred before the divergence of Primula and its close relatives, and multiple independent transfers could also occur in a single MTPT sequence. Foreign sequence [MTPTs and mitochondrial nuclear DNA sequences (NUMTs)] uptake and repeats were to some extent associated with changes in Ericales mitogenome size, although none of these relationships were significant overall. CONCLUSIONS: The present study revealed relatively conserved gene contents, gene clusters, RNA editing, and MTPTs but considerable structural variation in Primula mitogenomes. Relatively high nucleotide diversity was found in the Primula mitogenomes. In addition, mitogenomic genes, collinear gene clusters, and locally collinear blocks (LCBs) all showed phylogenetic signals. The evolutionary history of MTPTs in Primula was complicated, even in a single MTPT sequence. Various reasons for the size variation observed in Ericales mitogenomes were found.

Plastid phylogenomics of the Sansevieria Clade of Dracaena (Asparagaceae) resolves a recent radiation.

Best known as low maintenance houseplants, sansevierias are a diverse group of flowering plants native to Africa, Madagascar, the Arabian Peninsula, and the Indian subcontinent. Traditionally recognised as a distinct genus, Sansevieria was recently merged with the larger genus Dracaena based on molecular phylogenetic data. Within the Sansevieria Clade of Dracaena, taxonomic uncertainties remain despite attempts to unravel the relationships between the species. To investigate the evolutionary relationships, morphological evolution and biogeographical history in the group, we aim to reconstruct a robust dated phylogenetic hypothesis. Using genome skimming, a chloroplast genome (cpDNA) dataset and a nuclear ribosomal (nrDNA) dataset were generated. The sampling included representatives of all sections and informal groups previously described in Sansevieria based on morphology. Analysis of the cpDNA dataset using a maximum likelihood approach resulted in a well-supported phylogeny. The time-calibrated phylogeny indicated a recent radiation with five main clades emerging in the Pliocene. Two strongly supported clades align with previously defined groups, i.e., Sansevieria section Dracomima, characterised by the Dracomima-type inflorescence, and the Zeylanica informal group, native to the Indian subcontinent. Other previously defined groups were shown to be polyphyletic; a result of convergent evolution of the identifying characters. Switches between flat and cylindrical leaves occurred multiple times in the evolution of the Sansevieria Clade. Similarly, the Cephalantha-type inflorescence has originated multiple times from an ancestor with a Sansevieria-type inflorescence. Analysis of the nrDNA dataset resulted in a phylogenetic hypothesis with low resolution, yet it supported the same two groups confirmed by the cpDNA dataset. This study furthers our understanding of the evolution of the Sansevieria Clade, which will benefit taxonomic and applied research, and aid conservation efforts.

High-throughput degraded DNA sequencing of subfossil shells of a critically endangered stenoendemic land snail in the Aegean.

High-throughput sequencing has enabled the comprehensive genetic exploration of biological diversity, especially by using natural history collections to study hard-to-find, threatened or even extinct-in-the-wild taxa. Mollusk shells are under-exploited as a source for DNA-based approaches, despite their apparent advantages in the field of conservation genetics. More particularly, degraded DNA techniques combined with high-throughput sequencing have never been used to gain insights about the DNA preservation in land snail subfossil or historical shells. Here, we applied degraded DNA analysis on two historical shells of Levantina rechingeri, a stenoendemic Critically Endangered species that has never been found alive, in order to explore the patterns of DNA preservation on land snail shells originating from the eastern Mediterranean, as well as to infer its molecular phylogenetic placement. Our results showed that centuries to decades-old DNA from an empty shell collected in an Aegean island exhibits characteristic post-mortem damage patterns similar to those observed in ancient DNA from eastern Mediterranean terrestrial animals, setting a precedent for future museomics studies on taxa distributed in areas with similar climate. Finally, genome skimming of the empty shell allowed high coverage of multiple nuclear and mitochondrial loci, enabling the phylogenetic placement of the focal taxon, the re-evaluation of its taxonomic classification, and the revealing of a new Aegean land snail lineage, Aristena genus novum. This approach is a non-invasive way to sample DNA from threatened land snail species and suitable to study the evolutionary history of taxa with cryptic ecology, stenoendemics, or extinct-in-the-wild, as well as old museum specimens.

Resolution, conflict and rate shifts: insights from a densely sampled plastome phylogeny for Rhododendron (Ericaceae).

BACKGROUND AND AIMS: Rhododendron is a species-rich and taxonomically challenging genus due to recent adaptive radiation and frequent hybridization. A well-resolved phylogenetic tree would help to understand the diverse history of Rhododendron in the Himalaya-Hengduan Mountains where the genus is most diverse. METHODS: We reconstructed the phylogeny based on plastid genomes with broad taxon sampling, covering 161 species representing all eight subgenera and all 12 sections, including ~45 % of the Rhododendron species native to the Himalaya-Hengduan Mountains. We compared this phylogeny with nuclear phylogenies to elucidate reticulate evolutionary events and clarify relationships at all levels within the genus. We also estimated the timing and diversification history of Rhododendron, especially the two species-rich subgenera Rhododendron and Hymenanthes that comprise >90 % of Rhododendron species in the Himalaya-Hengduan Mountains. KEY RESULTS: The full plastid dataset produced a well-resolved and supported phylogeny of Rhododendron. We identified 13 clades that were almost always monophyletic across all published phylogenies. The conflicts between nuclear and plastid phylogenies suggested strongly that reticulation events may have occurred in the deep lineage history of the genus. Within Rhododendron, subgenus Therorhodion diverged first at 56 Mya, then a burst of diversification occurred from 23.8 to 17.6 Mya, generating ten lineages among the component 12 clades of core Rhododendron. Diversification in subgenus Rhododendron accelerated c. 16.6 Mya and then became fairly continuous. Conversely, Hymenanthes diversification was slow at first, then accelerated very rapidly around 5 Mya. In the Himalaya-Hengduan Mountains, subgenus Rhododendron contained one major clade adapted to high altitudes and another to low altitudes, whereas most clades in Hymenanthes contained both low- and high-altitude species, indicating greater ecological plasticity during its diversification. CONCLUSIONS: The 13 clades proposed here may help to identify specific ancient hybridization events. This study will help to establish a stable and reliable taxonomic framework for Rhododendron, and provides insight into what drove its diversification and ecological adaption. Denser sampling of taxa, examining both organelle and nuclear genomes, is needed to better understand the divergence and diversification history of Rhododendron.

A Pipeline for the Development of Microsatellite Markers using Next Generation Sequencing Data.

Background: Also known as Simple Sequence Repetitions (SSRs), microsatellites are profoundly informative molecular markers and powerful tools in genetics and ecology studies on plants. Objective: This research presents a workflow for developing microsatellite markers using genome skimming. Methods: The pipeline was proposed in several stages that must be performed sequentially: obtaining DNA sequences, identifying microsatellite regions, designing primers, and selecting candidate microsatellite regions to develop the markers. Our pipeline efficiency was analyzed using Illumina sequencing data from the non-model tree species Pterodon emarginatus Vog. Results: The pipeline revealed 4,382 microsatellite regions and drew 7,411 pairs of primers for P. emarginatus. However, a much larger number of microsatellite regions with the potential to develop markers were discovered from our pipeline. We selected 50 microsatellite regions with high potential for developing markers and organized 29 microsatellite regions in sets for multiplex PCR. Conclusion: The proposed pipeline is a powerful tool for fast and efficient development of microsatellite markers on a large scale in several species, especially nonmodel plant species.

Phylogenomic conflict analyses in the apple genus Malus s.l. reveal widespread hybridization and allopolyploidy driving diversification, with insights into the complex biogeographic history in the Northern Hemisphere.

Phylogenomic evidence from an increasing number of studies has demonstrated that different data sets and analytical approaches often reconstruct strongly supported but conflicting relationships. In this study, 785 single-copy nuclear genes and 75 complete plastomes were used to infer the phylogenetic relationships and estimate the historical biogeography of the apple genus Malus sensu lato, an economically important lineage disjunctly distributed in the Northern Hemisphere and involved in known and suspected hybridization and allopolyploidy events. The nuclear phylogeny recovered the monophyly of Malus s.l. (including Docynia); however, the genus was supported to be biphyletic in the plastid phylogeny. An ancient chloroplast capture event in the Eocene in western North America best explains the cytonuclear discordance. Our conflict analysis demonstrated that ILS, hybridization, and allopolyploidy could explain the widespread nuclear gene tree discordance. One deep hybridization event (Malus doumeri) and one recent event (Malus coronaria) were detected in Malus s.l. Furthermore, our historical biogeographic analysis integrating living and fossil data supported a widespread East Asian-western North American origin of Malus s.l. in the Eocene, followed by several extinction and dispersal events in the Northern Hemisphere. We also propose a general workflow for assessing phylogenomic discordance and biogeographic analysis using deep genome skimming data sets.

Chloroplast genome sequencing based on genome skimming for identification of Eriobotryae Folium.

BACKGROUND: Whole chloroplast genome (cpDNA) sequence is becoming widely used in the phylogenetic studies of plant and species identification, but in most cases the cpDNA were acquired from silica gel dried fresh leaves. So far few reports have been available to describe cpDNA acquisition from crude drugs derived from plant materials, the DNA of which usually was seriously damaged during their processing. In this study, we retrieved cpDNA from the commonly used crude drug Eriobotryae Folium (Pipaye in Chinese, which is the dried leaves of Eriobotrya japonica, PPY) using genome skimming technique. RESULTS: We successfully recovered cpDNA sequences and rDNA sequences from the crude drug PPY, and bioinformatics analysis showed a high overall consistency between the cpDNA obtained from the crude drugs and fresh samples. In the ML tree, each species formed distinct monophyletic clades based on cpDNA sequence data, while the phylogenetic relationships between Eriobotrya species were poorly resolved based on ITS and ITS2. CONCLUSION: Our results demonstrate that both cpDNA and ITS/ITS2 are effective for identifying PPY and its counterfeits derived from distantly related species (i.e. Dillenia turbinata and Magnolia grandiflora), but cpDNA is more effective for distinguishing the counterfeits derived from the close relatives of Eriobotrya japonica, suggesting the potential of genome skimming for retrieving cpDNA from crude drugs used in Traditional Chinese Medicine for their identification.

Large vs small genomes in Passiflora: the influence of the mobilome and the satellitome.

MAIN CONCLUSIONS: While two lineages of retrotransposons were more abundant in larger Passiflora genomes, the satellitome was more diverse and abundant in the smallest genome analysed. Repetitive sequences are ubiquitous and fast-evolving elements responsible for size variation and large-scale organization of plant genomes. Within Passiflora genus, a tenfold variation in genome size, not attributed to polyploidy, is known. Here, we applied a combined in silico and cytological approach to study the organization and diversification of repetitive elements in three species of this genus representing its known range in genome size variation. Sequences were classified in terms of type and repetitiveness and the most abundant were mapped to chromosomes. We identified long terminal repeat (LTR) retrotransposons as the most abundant elements in the three genomes, showing a considerable variation among species. Satellite DNAs (satDNAs) were less representative, but highly diverse between subgenera. Our results clearly confirm that the largest genome species (Passiflora quadrangularis) presents a higher accumulation of repetitive DNA sequences, specially Angela and Tekay elements, making up most of its genome. Passiflora cincinnata, with intermediate genome and from the same subgenus, showed similarity with P. quadrangularis regarding the families of repetitive DNA sequences, but in different proportions. On the other hand, Passiflora organensis, the smallest genome, from a different subgenus, presented greater diversity and the highest proportion of satDNA. Altogether, our data indicates that while large genomes evolved by an accumulation of retrotransposons, the smallest genome known for the genus has evolved by diversification of different repeat types, particularly satDNAs.