Differentiating between hypotheses of lineage sorting and introgression in New Zealand alpine cicadas (Maoricicada Dugdale)

Phylogeny of North Americas largest cicada radiation redefines Tibicinoides and Okanagana (Hemiptera: Auchenorrhyncha: Cicadidae: Tibicininae)

Tibicinoides, with three small endemic California cicada species, has a confusing, intertwined systematic history with Okanagana that we unravel here. An ingroup including all species of Tibicinoides and the majority (84.7%) of Okanagana species were sampled for six gene regions, polarized with Clidophleps, Okanagodes, Subpsaltria, and Tibicina outgroups, and subjected to Bayesian phylogenetic analysis. Although the ingroup was monophyletic from all outgroups including Tibicina, Tibicinoides rendered Okanagana paraphyletic among two major ingroup clades. To bring classification into agreement with phylogeny, we redescribe and redefine Tibicinoides to include all Okanagana species with a hooked uncus in the male genitalia, all of which grouped with the type T. cupreosparsa (Uhler, 1889) in the first of these clades: T. boweni (Chatfield-Taylor & Cole, 2020) comb. n., T. catalina (Davis, 1936) comb. n., T. hesperia (Uhler, 1876) comb. n., T. mercedita (Davis, 1915), T. minuta (Davis, 1915), T. pallidula (Davis, 1917a) comb. n., T. pernix (Bliven, 1964) comb. n., T. rubrovenosa (Davis, 1915) comb. n., T. simulata (Davis, 1921) comb. n., T. striatipes (Haldeman, 1852) comb. n., T. uncinata (Van Duzee, 1915) comb. n., T. utahensis (Davis, 1919) comb. n., and T. vanduzeei (Distant, 1914) comb. n. Okanagana is redescribed and restricted to the species of the second major clade which contained the type O. rimosa (Say, 1830). We describe two new genera for morphologically distinct orphan lineages: Chlorocanta gen. nov. for C. viridis (Davis, 1918) comb. n. and Hewlettia gen. nov. for H. nigriviridis (Davis, 1921) comb. n. We recognize O. rubrobasalis Davis, 1926 stat. rev. as a species and relegate two former species to junior subjective synonyms: O. noveboracensis (Emmons, 1854) = O. canadensis (Provancher, 1889) and O. occidentalis (Walker in Lord, 1866) = O. lurida Davis, 1919. Tibicinoides and Okanagana together represent a rapid radiation that presents challenges to phylogenetic analysis including suboptimal outgroups and short internodes.

The phylogenetic structure and coalescent species delimitation of an endemic trapdoor spider genus, Stasimopus (Araneae, Mygalomorphae, Stasimopidae) in the Karoo region of South Africa

The Karoo region of South Africa is a unique and sensitive ecosystem which is facing pressure for development due to economic incentives such as mining, farming and shale gas exploration. The species diversity of many taxa in the area is largely unknown. A phylogenetic analysis of the cork-lid trapdoor spider genus, Stasimopus (Stasimopidae) was undertaken in order to gain insight into the relationships between the species that may be present in the area. The species within Stasimopus are challenging to identify and define using traditional morphological methods due to a high degree of morphological conservatism within the genus. For this reason, multiple coalescent based species delimitation methods were used to attempt to determine the species present for Stasimopus in the region which was tested against the morphological identifications and genetic clades (based on CO1, 16S and EF-1ɣ). We tested single-locus methods Automatic Barcode Gap Discovery (ABGD), Bayesian implementation of Poisson Tree Processes (bPTP) and General Mixed Yule- Coalescent (GMYC), as well as multi-locus Brownie. The phylogenetic analysis of Stasimopus in the Karoo showed that there is a high degree of genetic diversity within the genus. The species delimitation results proved unfruitful for the genus, as they appear to delimit population structure rather than species for most methods. Alternative methods should be investigated to aid in the identification of the species in order truly understand the species diversity of the genus.

Phylogeny, biogeography, and character evolution of the genus Sophora s.l. (Fabaceae, Papilionoideae)

The papilionoid legume genus Sophora (Fabaceae) exhibits a worldwide distribution, but a phylogenetic framework to understand the evolution of this group is lacking to date. Previous studies have demonstrated that Sophora is not monophyletic and might include Ammodendron, Ammothamnus, and Echinosophora, but the relationships among these four genera (defined as Sophora s.l.) are unclear. Here we used a nuclear DNA dataset (ETS, ITS, SQD1) and a plastid DNA dataset (matK, rbcL, rpl32-trnL, trnL-F) of 654 accession sequences to reconstruct the phylogenetic relationships, estimate the divergence times and ancestral range of Sophora s.l., and infer the evolution of chromosome number and morphological characteristics. Our major aim was to reconstruct phylogenetic relationships to test monophyly and elucidate relationships within the genus. Our results indicated that Ammodendron, Ammothamnus, and Echinosophora are embedded within Sophora s.s. and that nine well-supported clades can be recognized within comprise Sophora s.l. Ancestral character state estimation revealed that the most recent common ancestor of Sophora s.l. was a deciduous shrub that lacks rhizome spines and has unwinged legumes. Divergence times estimation and ancestral area reconstruction showed that Sophora s.l. originated in Central Asia and/or adjacent Southeast China in the early Oligocene (ca. 31 Mya) and dispersed from these regions into East and South Asia's adjacent areas and North America via the Bering land bridge. The analyses also supported a South American origin for S. sect. Edwardsia, which experienced rapid radiation with its major lineages diversifying over a relatively narrow timescale (8 Mya).

Plastome phylogenomics provide new perspective into the phylogeny and evolution of Betulaceae (Fagales)

BACKGROUND

Betulaceae is a relatively small but morphologically diverse family, with many species having important economic and ecological values. Although plastome structure of Betulaceae has been reported sporadically, a comprehensive exploration for plastome evolution is still lacking. Besides, previous phylogenies had been constructed based on limited gene fragments, generating unrobust phylogenetic framework and hindering further studies on divergence ages, biogeography and character evolution. Here, 109 plastomes (sixteen newly assembled and 93 previously published) were subject to comparative genomic and phylogenomic analyses to reconstruct a robust phylogeny and trace the diversification history of Betulaceae.

RESULTS

All Betulaceae plastomes were highly conserved in genome size, gene order, and structure, although specific variations such as gene loss and IR boundary shifts were revealed. Ten divergent hotspots, including five coding regions (P_i > 0.02) and five noncoding regions (P_i > 0.035), were identified as candidate DNA barcodes for phylogenetic analysis and species delimitation. Phylogenomic analyses yielded high-resolution topology that supported reciprocal monophyly between Betula and Alnus within Betuloideae, and successive divergence of Corylus, Ostryopsis, and Carpinus-Ostrya within Coryloideae. Incomplete lineage sorting and hybridization may be responsible for the mutual paraphyly between Ostrya and Carpinus. Betulaceae ancestors originated from East Asia during the upper Cretaceous; dispersals and subsequent vicariance accompanied by historical environment changes contributed to its diversification and intercontinental disjunction. Ancestral state reconstruction indicated the acquisition of many taxonomic characters was actually the results of parallel or reversal evolution.

CONCLUSIONS

Our research represents the most comprehensive taxon-sampled and plastome-level phylogenetic inference for Betulaceae to date. The results clearly document global patterns of plastome structural evolution, and established a well-supported phylogeny of Betulaceae. The robust phylogenetic framework not only provides new insights into the intergeneric relationships, but also contributes to a perspective on the diversification history and evolution of the family.

Lack of host phylogenetic structure in the gut bacterial communities of New Zealand cicadas and their interspecific hybrids

Host-microbe interactions are intimately linked to eukaryotic evolution, particularly in sap-sucking insects that often rely on obligate microbial symbionts for nutrient provisioning. Cicadas (Cicadidae: Auchenorrhyncha) specialize on xylem fluid and derive many essential amino acids and vitamins from intracellular bacteria or fungi (Hodgkinia, Sulcia, and Ophiocordyceps) that are propagated via transmission from mothers to offspring. Despite the beneficial role of these non-gut symbionts in nutrient provisioning, the role of beneficial microbiota within the gut remains unclear. Here, we investigate the relative abundance and impact of host phylogeny and ecology on gut microbial diversity in cicadas using 16S ribosomal RNA gene amplicon sequencing data from 197 wild-collected cicadas and new mitochondrial genomes across 38 New Zealand cicada species, including natural hybrids between one pair of two species. We find low abundance and a lack of phylogenetic structure and hybrid effects but a significant role of elevation in explaining variation in gut microbiota.

Out of chaos: Phylogenomics of Asian Sonerileae

Sonerileae is a diverse Melastomataceae lineage comprising ca. 1000 species in 44 genera, with >70% of genera and species distributed in Asia. Asian Sonerileae are taxonomically intractable with obscure generic circumscriptions. The backbone phylogeny of this group remains poorly resolved, possibly due to complexity caused by rapid species radiation in early and middle Miocene, which hampers further systematic study. Here, we used genome resequencing data to reconstruct the phylogeny of Asian Sonerileae. Three parallel datasets, viz. single-copy ortholog (SCO), genomic SNPs, and whole plastome, were assembled from genome resequencing data of 205 species for this purpose. Based on these genome-scale data, we provided the first well resolved phylogeny of Asian Sonerileae, with 34 major clades identified and 74% of the interclade relationships consistently resolved by both SCO and genomic data. Meanwhile, widespread phylogenetic discordance was detected among SCO gene trees as well as species trees reconstructed using different tree estimation methods (concatenation/site-based coalescent method/summary method) or different datasets (SCO/genomic/plastome). We explored sources of discordance using multiple approaches and found that the observed discordance in Asian Sonerileae was mainly caused by a combination of biased distribution of missing data, random noise from uninformative genes, incomplete lineage sorting, and hybridization/introgression. Exploration of these sources can enable us to generate hypotheses for future testing, which is the first step towards understanding the evolution of Asian Sonerileae. We also detected high levels of homoplasy for some characters traditionally used in taxonomy, which explains current chaotic generic delimitations. The backbone phylogeny of Asian Sonerileae revealed in this study offers a solid basis for future taxonomic revision at the generic level.

Cryptic Species Diversification of the Pedicularis siphonantha Complex (Orobanchaceae) in the Mountains of Southwest China Since the Pliocene

Morphological approaches often fail to delimit species in recently derived species complexes. This can be exacerbated in historical collections which may have lost key features in specimen preparation and preservation. Here, we examine the Pedicularis siphonantha complex, endemic to the Mountains of Southwest China. This complex is characterized by its red/purple/pink and long-tubular corolla, and twisted, beaked galea. However, herbarium specimens are often difficult to identify to species. Molecular approaches using nrITS or nuclear ribosomal internal transcribed spacer (nrITS) + plastid DNA (ptDNA) have been successfully used for species identification in Pedicularis. To resolve taxonomic confusion in the Pedicularis siphonantha complex, we reconstructed phylogenies of the complex using nrITS and four plastid DNA loci (matK, rbcL, trnH-psbA, and trnL-F). To recover as much of the phylogenetic history as possible, we sampled individuals at the population level. Topological incongruence between the nrITS and ptDNA datasets was recovered in clades including two widely distributed species, Pedicularis milliana and Pedicularis tenuituba. Based on morphological, geographical, and genetic evidence, we suggest that hybridization/introgression has occurred between P. milliana and Pedicularis sigmoidea/Pedicularis sp. 1 in the Yulong Snow Mountain of Lijiang, northwest Yunnan, and between P. tenuituba and Pedicularis leptosiphon in Ninglang, northwest Yunnan. After removing conflicting DNA regions in Pedicularis dolichosiphon (nrITS) and P. milliana (ptDNA), the concatenated nrITS and ptDNA phylogenies distinguish 11 species in the P. siphonantha complex, including two undescribed species, from the Jiaozi and Yulong Snow Mountains, respectively. Phylogeographical analyses indicate that the P. siponantha complex originated from south of the Hengduan Mountains, expanding north to the Himalayas and the Yunnan-Guizhou Plateau. Moreover, the uplift of the Qinghai-Tibet Plateau and climate oscillations may have driven further diversification in the complex.

Utilizing museomics to trace the complex history and species boundaries in an avian-study system of conservation concern

A taxonomic classification that accurately captures evolutionary history is essential for conservation. Genomics provides powerful tools for delimiting species and understanding their evolutionary relationships. This allows for a more accurate and detailed view on conservation status compared with other, traditionally used, methods. However, from a practical and ethical perspective, gathering sufficient samples for endangered taxa may be difficult. Here, we use museum specimens to trace the evolutionary history and species boundaries in an Asian oriole clade. The endangered silver oriole has long been recognized as a distinct species based on its unique coloration, but a recent study suggested that it might be nested within the maroon oriole-species complex. To evaluate species designation, population connectivity, and the corresponding conservation implications, we assembled a de novo genome and used whole-genome resequencing of historical specimens. Our results show that the silver orioles form a monophyletic lineage within the maroon oriole complex and that maroon and silver forms continued to interbreed after initial divergence, but do not show signs of recent gene flow. Using a genome scan, we identified genes that may form the basis for color divergence and act as reproductive barriers. Taken together, our results confirm the species status of the silver oriole and highlight that taxonomic revision of the maroon forms is urgently needed. Our study demonstrates how genomics and Natural History Collections (NHC) can be utilized to shed light on the taxonomy and evolutionary history of natural populations and how such insights can directly benefit conservation practitioners when assessing wild populations.

Thirteen new species of Chilecicada Sanborn, 2014 (Hemiptera: Auchenorrhyncha: Cicadidae: Tibicininae) expand the highly endemic cicada fauna of Chile

The genus Chilecicada Sanborn, 2014 is shown to be a complex of closely related species rather than a monospecific genus. Chilecicada citatatemporaria Sanborn & Cole n. sp., C. culenesensis Sanborn & Cole n. sp., C. curacaviensis Sanborn & Cole n. sp., C. impartemporaria Sanborn & Cole n. sp., C. magna Sanborn & Cole n. sp., C. mapuchensis Sanborn n. sp., C. oraria Sanborn & Cole n. sp., C. parrajaraorum Sanborn n. sp., C. partemporaria Sanborn & Cole n. sp., C. pehuenchesensis Sanborn & Cole n. sp., C. trifascia Sanborn n. sp., C. trifasciunca Sanborn & Cole n. sp., and C. viridicitata Sanborn & Cole n. sp. are described as new. Chilecicada occidentis Walker, 1850 is re-described to facilitate separation of the new species from the only previously known species. Song and cytochrome oxidase I analysis available for most species support the separation of the new taxa from the type species of the genus. Known species distributions and a key to the species of the genus are also provided. The new species increases the known cicada diversity 61.9% to 34 species, 91.2% of which are endemic to Chile.

Phylogeography of the endemic red-tailed cicadas of New Zealand (Hemiptera: Cicadidae: Rhodopsalta), and molecular, morphological and bioacoustical confirmation of the existence of Hudson’s Rhodopsalta microdora

Why do some genera radiate, whereas others do not? The genetic structure of present-day populations can provide clues for developing hypotheses. In New Zealand, three Cicadidae genera are depauperate [Amphipsalta (three species), Notopsalta (one species) and Rhodopsalta (three species)], whereas two have speciated extensively [Kikihia (~30 species/subspecies) and Maoricicada (~20 species/subspecies). Here, we examine the evolution of Rhodopsalta, the last New Zealand genus to be studied phylogenetically and phylogeographically. We use Bayesian and maximum-likelihood analyses of mitochondrial cox1 and nuclear EF1α gene sequences. Concatenated and single-gene phylogenies for 70 specimens (58 localities) support its monophyly and three described species: Rhodopsalta cruentata, Rhodopsalta leptomera and Rhodopsalta microdora, the last taxon previously regarded as uncertain. We provide distribution maps, biological notes and the first descriptions of diagnostic songs. We show that both R. cruentata and R. microdora exhibit northern and southern genetic subclades. Subclades of the dry-adapted R. microdora clade show geographical structure, whereas those of the mesic R. cruentata and sand-dune specialist R. leptomera have few discernible patterns. Genetic, bioacoustical and detailed distributional evidence for R. microdora add to the known biodiversity of New Zealand. We designate a lectotype for Tettigonia cruentata  Fabricius, 1775, the type species of Rhodopsalta.

Repeat proliferation and partial endoreplication jointly shape the patterns of genome size evolution in orchids

Although the evolutionary drivers of genome size change are known, the general patterns and mechanisms of plant genome size evolution are yet to be established. Here we aim to assess the relative importance of proliferation of repetitive DNA, chromosomal variation (including polyploidy), and the type of endoreplication for genome size evolution of the Pleurothallidinae, the most species-rich orchid lineage. Phylogenetic relationships between 341 Pleurothallidinae representatives were refined using a target enrichment hybrid capture combined with high-throughput sequencing approach. Genome size and the type of endoreplication were assessed using flow cytometry supplemented with karyological analysis and low-coverage Illumina sequencing for repeatome analysis on a subset of samples. Data were analyzed using phylogeny-based models. Genome size diversity (0.2-5.1 Gbp) was mostly independent of profound chromosome count variation (2n = 12-90) but tightly linked with the overall content of repetitive DNA elements. Species with partial endoreplication (PE) had significantly greater genome sizes, and genomic repeat content was tightly correlated with the size of the non-endoreplicated part of the genome. In PE species, repetitive DNA is preferentially accumulated in the non-endoreplicated parts of their genomes. Our results demonstrate that proliferation of repetitive DNA elements and PE together shape the patterns of genome size diversity in orchids.

Genetic Introgression and Morphological Variation in Naked-Back Bats (Chiroptera: Mormoopidae: Pteronotus Species) along Their Contact Zone in Central America

Two sibling bare-backed bat species (Pteronotus fulvus and P. gymnonotus) have been traditionally differentiated by their size. However, intermediate specimens between the two species have been found in sympatric populations along southern Mexico and it has been suggested that they may be the outcome of a hybridization process between the two species. We used one mitochondrial (COI), three nuclear markers (PRKCL, STAT5A and RAG2) and 13 microsatellites to explore the evolutionary relationships between these two species and elucidate whether the intermediate morphotypes correspond to hybrid individuals. These markers have been analyzed in sympatric and allopatric populations of the two species plus the closely related species Pteronotus davyi. We confirmed the species-level differentiation of the three lineages (P. fulvus, P. davyi and P. gymnonotus), but the phylogenetic hypotheses suggested by the nuclear and mitochondrial markers were discordant. We confirm that the discordance between markers is due to genetic introgression through the mitochondrial capture of P. fulvus in P. gymnonotus populations. Such introgression was found in all P. gymnonotus specimens across its sympatric distribution range (Mexico to Costa Rica) and is related to expansion/retraction species distribution pulses associated with changes in forest distribution during the Quaternary climate cycles. Microsatellite analyses showed contemporary genetic contact between the two sympatric species and 3.0% of the samples studied were identified as hybrids. In conclusion, we found a historical and asymmetric genetic introgression (through mitochondrial capture) of P. fulvus into P. gymnonotus in Mexico and Central America and a limited contemporary gene exchange between the two species. However, no relationship was found between hybridization and the intermediate-sized specimens from southern Mexico, which might likely result from a clinal variation with latitude. These results confirm the need for caution when using forearm size to identify these species in the field and when differentiating them in the laboratory based on mitochondrial DNA alone.

Hidden in plain sight: Morphological and phylogenetic evidence for Bouteloua arizonica, a species distinct from Bouteloua aristidoides (Poaceae, Chloridoideae)

Two varieties of Bouteloua aristidoides have been recognized, the widespread var. aristidoides and the more narrowly distributed var. arizonica. The two varieties differ in inflorescence form even more than that seen between many other closely related species of Bouteloua. We therefore asked whether these taxa might be better regarded as distinct species. A total of 93 vouchers were studied by using morphometry (principal components analysis and statistical tests), leaf micromorphology, ancestral state reconstruction, and/or molecular (ITS, trnC-rpoB and trnT-L-F) phylogenetic analyses. Except from the ITS tree, all results supported elevation of B. aristidoides var. arizonica to the rank of species, thus the new combination Bouteloua arizonica (M.E. Jones) L.F. Cuellar & Columbus comb. nov. et stat. nov., is proposed. Chloroplast and combined chloroplast-nuclear molecular trees depicted var. arizonica as monophyletic (even in sympatric populations with var. aristidoides) and reveals phylogenetic structure within var. aristidoides for which the presence of new undescribed varieties of B. aristidoides (different from B. arizonica) is addressed. B. arizonica differs from B. aristidoides in having fewer branches per inflorescence, a bigger branch with more spikelets, and a shorter branch extension. Scanning electron microscopy revealed the presence of papillae on leaves of B. arizonica as a clear synapomorphy. Growing mature plants of B. arizonica from seeds in a greenhouse revealed a strong cleistogamous nature for this species for which gene flow in sympatric populations with B. aristidoides seems unlikely. A taxonomic treatment and distribution map for identification of B. arizonica is provided.

Parallel ddRAD and Genome Skimming Analyses Reveal a Radiative and Reticulate Evolutionary History of the Temperate Bamboos

Rapid evolutionary radiations are among the most challenging phylogenetic problems, wherein different types of data (e.g., morphology and molecular) or genetic markers (e.g., nuclear and organelle) often yield inconsistent results. The tribe Arundinarieae, that is, the temperate bamboos, is a clade of tetraploid originated 22 Ma and subsequently radiated in East Asia. Previous studies of Arundinarieae have found conflicting relationships and/or low support. Here, we obtain nuclear markers from ddRAD data for 213 Arundinarieae taxa and parallel sampling of chloroplast genomes from genome skimming for 147 taxa. We first assess the feasibility of using ddRAD-seq data for phylogenetic estimates of paleopolyploid and rapidly radiated lineages, optimize clustering thresholds, and analysis workflow for orthology identification. Reference-based ddRAD data assembly approaches perform well and yield strongly supported relationships that are generally concordant with morphology-based taxonomy. We recover five major lineages, two of which are notable (the pachymorph and leptomorph lineages), in that they correspond with distinct rhizome morphologies. By contrast, the phylogeny from chloroplast genomes differed significantly. Based on multiple lines of evidence, the ddRAD tree is favored as the best species tree estimation for temperate bamboos. Using a time-calibrated ddRAD tree, we find that Arundinarieae diversified rapidly around the mid-Miocene corresponding with intensification of the East Asian monsoon and the evolution of key innovations including the leptomorph rhizomes. Our results provide a highly resolved phylogeny of Arundinarieae, shed new light on the radiation and reticulate evolutionary history of this tribe, and provide an empirical example for the study of recalcitrant plant radiations. [Arundinarieae; ddRAD; paleopolyploid; genome skimming; rapid diversification; incongruence.]

Repeated intercontinental migrations and recurring hybridizations characterise the evolutionary history of yew (Taxus L.)

The genus Taxus (Taxaceae) consists of 16 genetically well-defined lineages that are predominantly distributed across the Northern hemisphere. We investigated its biogeographic origin and evolutionary history by sampling 13 chloroplast gene sequences, the nuclear internal transcribed spacers (ITS) and NEEDLY sequences for all 16 lineages. We applied Maximum Parsimony and Bayesian Inference analyses to infer their phylogenetic relationships, time-calibrated phylogenies using BEAST and inferred the ancestral area of occupancy with BioGeoBEARS. We found strong evidence for the hybrid origin of three lineages and dated these events to a rather narrow time window of 6.8-4.9 million years ago (Mya). The dated phylogenies inferred an Upper Cretaceous origin of the genus, with the extant lineages diversifying in North America much later during the Oligocene/early Miocene. Repeated migrations via the Bering land bridge to Eurasia and back were further inferred, with the return to North America as a possible result of vicariance. The diversification in Eurasia (from ~8 Mya onwards) coincided with the orogeny of the Hengduan Mountains, the intensification of the East Asian summer monsoon and the occupancy of ecological niches by lineages that experienced secondary contacts and hybridizations in the Hengduan Mountains and Qinling Mountain, especially around the Sichuan basin. We provide a hypothesis for the evolution of extant lineages of Taxus, a genus with an old and complex evolutionary history. The study highlights that the history of complex species can be unravelled with a careful dissection of phylogenetic signals.

Disentangling Sources of Gene Tree Discordance in Phylogenomic Data Sets: Testing Ancient Hybridizations in Amaranthaceae s.l

Gene tree discordance in large genomic data sets can be caused by evolutionary processes such as incomplete lineage sorting and hybridization, as well as model violation, and errors in data processing, orthology inference, and gene tree estimation. Species tree methods that identify and accommodate all sources of conflict are not available, but a combination of multiple approaches can help tease apart alternative sources of conflict. Here, using a phylotranscriptomic analysis in combination with reference genomes, we test a hypothesis of ancient hybridization events within the plant family Amaranthaceae s.l. that was previously supported by morphological, ecological, and Sanger-based molecular data. The data set included seven genomes and 88 transcriptomes, 17 generated for this study. We examined gene-tree discordance using coalescent-based species trees and network inference, gene tree discordance analyses, site pattern tests of introgression, topology tests, synteny analyses, and simulations. We found that a combination of processes might have generated the high levels of gene tree discordance in the backbone of Amaranthaceae s.l. Furthermore, we found evidence that three consecutive short internal branches produce anomalous trees contributing to the discordance. Overall, our results suggest that Amaranthaceae s.l. might be a product of an ancient and rapid lineage diversification, and remains, and probably will remain, unresolved. This work highlights the potential problems of identifiability associated with the sources of gene tree discordance including, in particular, phylogenetic network methods. Our results also demonstrate the importance of thoroughly testing for multiple sources of conflict in phylogenomic analyses, especially in the context of ancient, rapid radiations. We provide several recommendations for exploring conflicting signals in such situations. [Amaranthaceae; gene tree discordance; hybridization; incomplete lineage sorting; phylogenomics; species network; species tree; transcriptomics.]

Evolution of Termite Symbiosis Informed by Transcriptome-Based Phylogenies

Termitidae comprises ∼80% of all termite species [1] that play dominant decomposer roles in tropical ecosystems [2, 3]. Two major events during termite evolution were the loss of cellulolytic gut protozoans in the ancestor of Termitidae and the subsequent gain in the termitid subfamily Macrotermitinae of fungal symbionts cultivated externally in "combs" constructed within the nest [4, 5]. How these symbiotic transitions occurred remains unresolved. Phylogenetic analyses of mitochondrial data previously suggested that Macrotermitinae is the earliest branching termitid lineage, followed soon after by Sphaerotermitinae [6], which cultivates bacterial symbionts on combs inside its nests [7]. This has led to the hypothesis that comb building was an important evolutionary step in the loss of gut protozoa in ancestral termitids [8]. We sequenced genomes and transcriptomes of 55 termite species and reconstructed phylogenetic trees from up to 4,065 orthologous genes of 68 species. We found strong support for a novel sister-group relationship between the bacterial comb-building Sphaerotermitinae and fungus comb-building Macrotermitinae. This key finding indicates that comb building is a derived trait within Termitidae and that the creation of a comb-like "external rumen" involving bacteria or fungi may not have driven the loss of protozoa from ancestral termitids, as previously hypothesized. Instead, associations with gut prokaryotic symbionts, combined with dietary shifts from wood to other plant-based substrates, may have played a more important role in this symbiotic transition. Our phylogenetic tree provides a platform for future studies of comparative termite evolution and the evolution of symbiosis in this taxon.

Incomplete lineage sorting and introgression in the diversification of Chinese spot-billed ducks and mallards

Incomplete lineage sorting and introgression are 2 major and nonexclusive causes of species-level non-monophyly. Distinguishing between these 2 processes is notoriously difficult because they can generate similar genetic signatures. Previous studies have suggested that 2 closely related duck species, the Chinese spot-billed duck Anas zonorhyncha and the mallard A. platyrhynchos were polyphyletically intermixed. Here, we utilized a wide geographical sampling, multilocus data and a coalescent-based model to revisit this system. Our study confirms the finding that Chinese spot-billed ducks and Mallards are not monophyletic. There was no apparent interspecific differentiation across loci except those at the mitochondrial DNA (mtDNA) control region and the Z chromosome (CHD1Z). Based on an isolation-with-migration model and the geographical distribution of lineages, we suggest that both introgression and incomplete lineage sorting might contribute to the observed non-monophyly of the 2 closely related duck species. The mtDNA introgression was asymmetric, with high gene flow from Chinese spot-billed ducks to Mallards and negligible gene flow in the opposite direction. Given that the 2 duck species are phenotypically distinctive but weakly genetically differentiated, future work based on genome-scale data is necessary to uncover genomic regions that are involved in divergence, and this work may provide further insights into the evolutionary histories of the 2 species and other waterfowls.

Phylotranscriptomic Insights into the Diversification of Endothermic Thunnus Tunas

Birds, mammals, and certain fishes, including tunas, opahs and lamnid sharks, are endothermic, conserving internally generated, metabolic heat to maintain body or tissue temperatures above that of the environment. Bluefin tunas are commercially important fishes worldwide, and some populations are threatened. They are renowned for their endothermy, maintaining elevated temperatures of the oxidative locomotor muscle, viscera, brain and eyes, and occupying cold, productive high-latitude waters. Less cold-tolerant tunas, such as yellowfin tuna, by contrast, remain in warm-temperate to tropical waters year-round, reproducing more rapidly than most temperate bluefin tuna populations, providing resiliency in the face of large-scale industrial fisheries. Despite the importance of these traits to not only fisheries but also habitat utilization and responses to climate change, little is known of the genetic processes underlying the diversification of tunas. In collecting and analyzing sequence data across 29,556 genes, we found that parallel selection on standing genetic variation is associated with the evolution of endothermy in bluefin tunas. This includes two shared substitutions in genes encoding glycerol-3 phosphate dehydrogenase, an enzyme that contributes to thermogenesis in bumblebees and mammals, as well as four genes involved in the Krebs cycle, oxidative phosphorylation, β-oxidation, and superoxide removal. Using phylogenetic techniques, we further illustrate that the eight Thunnus species are genetically distinct, but found evidence of mitochondrial genome introgression across two species. Phylogeny-based metrics highlight conservation needs for some of these species.

Influence of a climatic gradient on genetic exchange between two oak species

PREMISE

In plant groups with limited intrinsic barriers to gene flow, it is thought that environmental conditions can modulate interspecific genetic exchange. Oaks are known for limited barriers to gene flow among closely related species. Here, we use Quercus as a living laboratory in which to pursue a fundamental question in plant evolution: Do environmental gradients restrict or promote genetic exchange between species?

METHODS

We focused on two North American oaks, the rare Quercus dumosa and the widespread Q. berberidifolia. We sampled intensively along a contact zone in California, USA. We sequenced restriction site-associated DNA markers and measured vegetative phenotype. We tested for genetic exchange, the association with climate, and the effect on phenotype.

RESULTS

There is evidence for genetic exchange between the species. Admixed plants are found in areas of intermediate climate, while less admixed plants are found at the extremes of the climatic gradient. Genetic and phenotypic patterns are out of phase in the contact zone; some plants display the phenotype of one species but are genetically associated with another.

CONCLUSIONS

Our results support the hypothesis that a strong climatic gradient can promote genetic exchange between species. The overall weak correlation between genotype and phenotype in the contact zone between the species suggests that genetic exchange can lead to the breakdown of trait combinations used to define species. This incongruency predicts ongoing problems for conservation of Q. dumosa, with implications for conservation of other oaks.

Lima A, Nakamura AM, Fernandes F, Sobrinho I, de Brito RA. Evidence for Introgression Among Three Species of the Anastrepha fraterculus Group, a Radiating Species Complex of Fruit Flies

Introgression should no longer be considered as rare a phenomenon as once thought, since several studies have recently documented gene flow between closely related and radiating species. Here, we investigated evolutionary relationships among three closely related species of fruit flies of the Anastrepha fraterculus group (Anastrepha fraterculus, A. obliqua and A. sororcula). We sequenced a set of 20 genes and implemented a combined populational and phylogenetic inference with a model selection approach by an ABC framework in order to elucidate the demographic history of these species. The phylogenetic histories inferred from most genes showed a great deal of discordance and substantial shared polymorphic variation. The analysis of several population and speciation models reveal that this shared variation is better explained by introgression rather than convergence by parallel mutation or incomplete lineage sorting. Our results consistently showed these species evolving under an isolation with migration model experiencing a continuous and asymmetrical pattern of gene flow involving all species pairs, even though still showed a more closely related relationship between A. fraterculus and A. sororcula when compared with A. obliqua. This suggests that these species have been exchanging genes since they split from their common ancestor ∼2.6 MYA ago. We also found strong evidence for recent population expansion that appears to be consequence of anthropic activities affecting host crops of fruit flies. These findings point that the introgression here found may have been driven by genetic drift and not necessary by selection, which has implications for tracking and managing fruit flies.

Phylogenomic analyses reveal a deep history of hybridization and polyploidy in the Neotropical genus Lachemilla (Rosaceae)

Hybridization, incomplete lineage sorting, and phylogenetic error produce similar incongruence patterns, representing a great challenge for phylogenetic reconstruction. Here, we use sequence capture data and multiple species tree and species network approaches to resolve the backbone phylogeny of the Neotropical genus Lachemilla, while distinguishing among sources of incongruence. We used 396 nuclear loci and nearly complete plastome sequences from 27 species to clarify the relationships among the major groups of Lachemilla, and explored multiple sources of conflict between gene trees and species trees inferred with a plurality of approaches. All phylogenetic methods recovered the four major groups previously proposed for Lachemilla, but species tree methods recovered different topologies for relationships between these four clades. Species network analyses revealed that one major clade, Orbiculate, is likely of ancient hybrid origin, representing one of the main sources of incongruence among the species trees. Additionally, we found evidence for a potential whole genome duplication event shared by Lachemilla and allied genera. Lachemilla shows clear evidence of ancient and recent hybridization throughout the evolutionary history of the group. Also, we show the necessity to use phylogenetic network approaches that can simultaneously accommodate incomplete lineage sorting and gene flow when studying groups that show patterns of reticulation.

Using Genomic Location and Coalescent Simulation to Investigate Gene Tree Discordance in Medicago L

Several well-documented evolutionary processes are known to cause conflict between species-level phylogenies and gene-level phylogenies. Three of the most challenging processes for species tree inference are incomplete lineage sorting, hybridization and gene duplication, which may result in unwarranted comparisons of paralogous genes. Several existing methods have dealt with these processes but none has yet been able to untangle all three at once. Here, we propose a stepwise method by which these processes can be discerned using information on genomic location coupled with coalescent simulations. In the first step, highly discordant genes within genomic blocks (putative paralogs) are identified and excluded from the data set and, in the second step, blocks of linked genes are grouped according to their hybrid history. Existing multispecies coalescent software can then be applied to recover the principal tree(s) that make up the species tree/network without violating the underlying model. The potential of the approach is evaluated on simulated data derived from a species network composed of nine species, of which one is of hybrid origin, and displaying a single-gene duplication that leads to paralogous comparisons. We apply our method to an empirical set of 12 genes from 7 species sampled in the plant genus Medicago that display phylogenetic discordance. We identify the causes of the discordance and demonstrate that the Medicago orbicularis lineage experienced an episode of ancient hybridization. Our results show promise as a new way to explore phylogenetic sequence data that can significantly improve species tree inference in presence of hybridization and undetected paralogy or other causes leading to extremely discordant gene trees. [Coalescent simulation; gene tree; genomic location; hybridization; incomplete lineage sorting; paralogy; phylogenetic incongruence; principal tree; species tree.].

Deep reticulation and incomplete lineage sorting obscure the diploid phylogeny of rain-lilies and allies (Amaryllidaceae tribe Hippeastreae)

Hybridization is a frequent and important force in plant evolution. Next-generation sequencing (NGS) methods offer new possibilities for clade resolution and ambitious sampling of gene genealogies, yet difficulty remains in detecting deep reticulation events using currently available methods. We reconstructed the phylogeny of diploid representatives of Amaryllidaceae tribe Hippeastreae to test the hypothesis of ancient hybridizations preceding the radiation of its major subclade, Hippeastrinae. Through hybrid enrichment of DNA libraries and NGS, we obtained data for 18 nuclear loci through a curated assembly approach and nearly complete plastid genomes for 35 ingroup taxa plus 5 outgroups. Additionally, we obtained alignments for 39 loci through an automated assembly algorithm. These data were analyzed with diverse phylogenetic methods, including concatenation, coalescence-based species tree estimation, Bayesian concordance analysis, and network reconstructions, to provide insights into the evolutionary relationships of Hippeastreae. Causes for gene tree heterogeneity and cytonuclear discordance were examined through a Bayesian posterior predictive approach (JML) and coalescent simulations. Two major clades were found, Hippeastrinae and Traubiinae, as previously reported. Our results suggest the presence of two major nuclear lineages in Hippeastrinae characterized by different chromosome numbers: (1) Tocantinia and Hippeastrum with 2n=22, and (2) Eithea, Habranthus, Rhodophiala, and Zephyranthes mostly with 2n=12, 14, and 18. Strong cytonuclear discordance was confirmed in Hippeastrinae, and a network scenario with at least six hybridization events is proposed to reconcile nuclear and plastid signals, along a backbone that may also have been affected by incomplete lineage sorting at the base of each major subclade.

Complex within a Complex: Integrative Taxonomy Reveals Hidden Diversity in Cicadetta brevipennis (Hemiptera: Cicadidae) and Unexpected Relationships with a Song Divergent Relative

Multiple sources of data in combination are essential for species delimitation and classification of difficult taxonomic groups. Here we investigate a cicada taxon with unusual cryptic diversity and we attempt to resolve seemingly contradictory data sets. Cicada songs act as species-specific premating barriers and have been used extensively to reveal hidden taxonomic diversity in morphologically similar species. The Palaearctic Cicadetta montana species complex is an excellent example where distinct song patterns have disclosed multiple recently described species. Indeed, two taxa turned out to be especially diverse in that they form a “complex within the complex”: the Cicadetta cerdaniensis song group (four species studied previously) and Cicadetta brevipennis (examined in details here). Based on acoustic, morphological, molecular, ecological and spatial data sampled throughout their broad European distribution, we find that Cicadetta brevipennis s. l. comprises five lineages. The most distinct lineage is identified as Cicadetta petryi Schumacher, 1924, which we re-assign to the species level. Cicadetta brevipennis litoralis Puissant & Hertach ssp. n. and Cicadetta brevipennis hippolaidica Hertach ssp. n. are new to science. The latter hybridizes with Cicadetta brevipennis brevipennis Fieber, 1876 at a zone inferred from intermediate song patterns. The fifth lineage requires additional investigation. The C. cerdaniensis and the C. brevipennis song groups exhibit characteristic, clearly distinct basic song patterns that act as reproductive barriers. However, they remain completely intermixed in the Bayesian and maximum likelihood COI and COII mitochondrial DNA phylogenies. The closest relative of each of the four cerdaniensis group species is a brevipennis group taxon. In our favoured scenario the phylogenetic pairs originated in common Pleistocene glacial refuges where the taxa speciated and experienced sporadic inter-group hybridization leading to extensive introgression and mitochondrial capture.

Hybridization among distantly related species: Examples from the polyploid genus Curcuma (Zingiberaceae)

Discerning relationships among species evolved by reticulate and/or polyploid evolution is not an easy task, although it is widely discussed. The economically important genus Curcuma (ca. 120 spp.; Zingiberaceae), broadly distributed in tropical SE Asia, is a particularly interesting example of a group of palaeopolyploid origin whose evolution is driven mainly by hybridization and polyploidization. Although a phylogeny and a new infrageneric classification of Curcuma, based on commonly used molecular markers (ITS and cpDNA), have recently been proposed, significant evolutionary questions remain unresolved. We applied a multilocus approach and a combination of modern analytical methods to this genus to distinguish causes of gene tree incongruence and to identify hybrids and their parental species. Five independent regions of nuclear DNA (DCS, GAPDH, GLOBOSA3, LEAFY, ITS) and four non-coding cpDNA regions (trnL-trnF, trnT-trnL, psbA-trnH and matK), analysed as a single locus, were employed to construct a species tree and hybrid species trees using (*)BEAST and STEM-hy. Detection of hybridogenous species in the dataset was also conducted using the posterior predictive checking approach as implemented in JML. The resulting species tree outlines the relationships among major evolutionary lineages within Curcuma, which were previously unresolved or which conflicted depending upon whether they were based on ITS or cpDNA markers. Moreover, by using the additional markers in tests of plausible topologies of hybrid species trees for C. vamana, C. candida, C. roscoeana and C. myanmarensis suggested by previous molecular and morphological evidence, we found strong evidence that all the species except C. candida are of subgeneric hybrid origin.

Phylogeny and biogeography of Primula sect. Armerina: implications for plant evolution under climate change and the uplift of the Qinghai-Tibet Plateau

BACKGROUND

The historical orogenesis and associated climatic changes of mountain areas have been suggested to partly account for the occurrence of high levels of biodiversity and endemism. However, their effects on dispersal, differentiation and evolution of many groups of plants are still unknown. In this study, we examined the detailed diversification history of Primula sect. Armerina, and used biogeographic analysis and macro-evolutionary modeling to investigate a series of different questions concerning the evolution of the geographical and ecological distribution of the species in this section.

RESULTS

We sequenced five chloroplast and one nuclear genes for species of Primula sect. Armerina. Neither chloroplast nor nuclear trees support the monophyly of the section. The major incongruences between the two trees occur among closely related species and may be explained by hybridization. Our dating analyses based on the chloroplast dataset suggest that this section began to diverge from its relatives around 3.55 million years ago, largely coinciding with the last major uplift of the Qinghai-Tibet Plateau (QTP). Biogeographic analysis supports the origin of the section in the Himalayan Mountains and dispersal from the Himalayas to Northeastern QTP, Western QTP and Hengduan Mountains. Furthermore, evolutionary models of ecological niches show that the two P. fasciculata clades have significantly different climatic niche optima and rates of niche evolution, indicating niche evolution under climatic changes and further providing evidence for explaining their biogeographic patterns.

CONCLUSION

Our results support the hypothesis that geologic and climatic events play important roles in driving biological diversification of organisms in the QTP area. The Pliocene uplift of the QTP and following climatic changes most likely promoted both the inter- and intraspecific divergence of Primula sect. Armerina. This study also illustrates how niche evolution under climatic changes influences biogeographic patterns.

Phylogeny and classification of the East Asian Amitostigma alliance (Orchidaceae: Orchideae) based on six DNA markers

BACKGROUND

Tribe Orchideae dominates the orchid flora of the temperate Northern Hemisphere but its representatives in East Asia had been subject to less intensive phylogenetic study than those in Eurasia and North America. Although this situation was improved recently by the molecular phylogenetic study of Jin et al., comparatively few species were analyzed from the species-rich and taxonomically controversial East Asian Amitostigma alliance. Here, we present a framework nrITS tree of 235 accessions of Orchideae plus an in-depth analysis of 110 representative accessions, encompassing most widely recognized species within the alliance, to elucidate their relationships.

RESULTS

We used parsimony, likelihood and Bayesian approaches to generate trees from data for two nuclear (nrITS, low-copy Xdh) and four chloroplast (matK, psbA-trnH, trnL-F, trnS-trnG) markers. Nuclear and plastid data were analyzed separately due to a few hard incongruences that most likely reflect chloroplast capture. Our results suggest key phylogenetic placements for Sirindhornia and Brachycorythis, and confirm previous assertions that the Amitostigma alliance is monophyletic and sister to the Eurasian plus European clades of subtribe Orchidinae. Seven robust clades are evident within the alliance, but none corresponds precisely with any of the traditional genera; the smaller and more morphologically distinct genera Tsaiorchis, Hemipilia, Neottianthe and Hemipiliopsis are monophyletic but each is nested within a polyphyletic plexus of species attributed to either Ponerorchis or the most plesiomorphic genus, Amitostigma. Two early-divergent clades that escaped analysis by Jin et al. undermine their attempt to circumscribe an expanded monophyletic genus Ponerorchis.

CONCLUSIONS

We provide a new framework on the complex phylogenetic relationships between Amitostigma and other genera traditionally included in its alliance; based on which, we combine the entire Amitostigma alliance into a morphologically and molecularly circumscribed Amitostigma sensu latissimo that also contains seven molecularly circumscribed sections. Our molecular trees imply unusually high levels of morphological homoplasy, but these will need to be quantified via a future group-wide review of the alliance based on living plants if morphology is to be fully integrated into our classification.

Distinguishing between incomplete lineage sorting and genomic introgressions: complete fixation of allospecific mitochondrial DNA in a sexually reproducing fish (Cobitis; Teleostei), despite clonal reproduction of hybrids

Distinguishing between hybrid introgression and incomplete lineage sorting causing incongruence among gene trees in that they exhibit topological differences requires application of statistical approaches that are based on biologically relevant models. Such study is especially challenging in hybrid systems, where usual vectors mediating interspecific gene transfers--hybrids with Mendelian heredity--are absent or unknown. Here we study a complex of hybridizing species, which are known to produce clonal hybrids, to discover how one of the species, Cobitis tanaitica, has achieved a pattern of mito-nuclear mosaic genome over the whole geographic range. We appplied three distinct methods, including the method using solely the information on gene tree topologies, and found that the contrasting mito-nuclear signal might not have resulted from the retention of ancestral polymorphism. Instead, we found two signs of hybridization events related to C. tanaitica; one concerning nuclear gene flow and the other suggested mitochondrial capture. Interestingly, clonal inheritance (gynogenesis) of contemporary hybrids prevents genomic introgressions and non-clonal hybrids are either absent or too rare to be detected among European Cobitis. Our analyses therefore suggest that introgressive hybridizations are rather old episodes, mediated by previously existing hybrids whose inheritance was not entirely clonal. Cobitis complex thus supports the view that the type of resulting hybrids depends on a level of genomic divergence between sexual species.

Assortative mating between two sympatric closely-related specialists: inferred from molecular phylogenetic analysis and behavioral data

Host plant shifting of phytophagous insects can lead to the formation of host associated differentiation and ultimately speciation. In some cases, host plant specificity alone acts as a nearly complete pre-mating isolating barrier among insect populations. We here test whether effective pre-mating isolation and host-independent behavioral isolation have evolved under the condition of extreme host specilization using two sympatric flea beetles with incomplete post-mating isolation under laboratory conditions. Phylogenetic analysis and coalescent simulation results showed that there is a limited interspecific gene flow, indicating effctive isolation between these species. Three types of mating tests in the absence of host plant cues showed that strong host-independent behavioral isolation has evolved between them. We conclude that almost perfect assortative mating between these two extreme host specialists results from a combination of reduced encounter rates due to differential host preference and strong sexual isolation.

The phylogeography of Eurasian Fraxinus species reveals ancient transcontinental reticulation

To investigate the biogeographical history of ashes species of the Eurasian section Fraxinus and to test the hypothesis of ancient reticulations, we sequenced nuclear DNA (nETS and nITS, 1075 bp) for 533 samples and scored AFLP for 63 samples of Eurasian ashes within the section Fraxinus. The nITS phylogeny retrieved the classical view of the evolution of the section, whereas nETS phylogeny indicated an unexpected separation of F. angustifolia in two paraphyletic groups, respectively found in southeastern Europe and in the other parts of the Mediterranean basin. In the nETS phylogeny, the former group was closely related to F. excelsior, whereas the later was closely related to F. mandshurica, a species which is restricted nowadays to northeastern Asia. This topological incongruence between the two loci indicated the occurrence of an ancient reticulation between European and Asian ash species. Several other ancient reticulation events between the two European species and the other species of the section were supported by the posterior predictive checking method. Some of these reticulation events would have occurred during the Miocene, when climatic variations may have lead these species to expand their distribution range and come into contact. The recurrent reticulations observed among Eurasian ash species indicate that they should be considered as conspecific taxa, with subspecific status for some groups. Altogether, the results of the present study provide a rare documented evidence for the occurrence of multiple ancient reticulations within a group of temperate tree taxa with modern disjunct distributions in Eurasia. These ancient reticulation events indicate that the speciation process is slow in ashes, necessitating long periods of geographical isolation. The implications for speciation processes in temperate trees with similar life history and reproductive biology are discussed.