Gene and species trees reveal mitochondrial and nuclear discordance in the Drosophila cardini group (Diptera: Drosophilidae)

Abstract. The Drosophila cardini group includes 15 species, which are subdivided into the D. cardini and D. dunni subgroups. Although many phylogenetic hypotheses have been proposed for this group during the last five decades (based on patterns of reproductive isolation, morphology of male genitalia, chromosomal inversions, isozyme variation, or molecular sequence data), these are mostly discordant with each other. We aimed to clarify some of the evolutionary patterns related to the origin of this incongruence, while also attempting to provide a better-supported phylogenetic hypothesis for the D. cardini subgroup. For this purpose, sequences from three mitochondrial and three nuclear loci were gathered for at least eight species, and both individual gene trees and joint species tree estimates were evaluated. Although there was concordance among gene trees within each of the nuclear and mitochondrial sets, considerable incongruence was revealed in the comparisons between these two data sets. The branching position of D. neocardini was the main source of incongruence, and species trees reconstructed using different approaches with and without this species were particularly incongruent. In addition to providing a better approximation of the evolutionary history of the D. cardini group, this study suggests that incomplete lineage sorting or introgression may be biasing previous species tree estimates. More generally, the results also suggest that the use of supermatrix methods can lead to an overestimation of support for the inferred relationships, and highlight the potential effects of different taxon sampling strategies in phylogenetic reconstruction.     

Phylogeny of Mysis (Crustacea, Mysida): history of continental invasions inferred from molecular and morphological data

We studied the phylogenetic history of opossum shrimps of the genus Mysis Latreille, 1802 (Crustacea: Mysida) using parsimony analyses of morphological characters, DNA sequence data from mitochondrial (16S, COI and CytB) and nuclear genes (ITS2, 18S), and eight allozyme loci. With these data we aimed to resolve a long‐debated question of the origin of the non‐marine (continental) taxa in the genus, i.e., “glacial relicts” in circumpolar postglacial lakes and “arctic immigrants” in the Caspian Sea. A simultaneous analysis of the data sets gave a single tree supporting monophyly of all continental species, as well as monophyly of the taxa from circumpolar lakes and from the Caspian Sea. A clade of three circumarctic marine species was sister group to the continental taxa, whereas Atlantic species had more distant relationships to the others. Small molecular differentiation among the morphologically diverse endemic species from the Caspian Sea suggested their recent speciation, while the phenotypically more uniform “glacial relict” species from circumpolar lakes (Mysis relicta group) showed deep molecular divergences. For the length‐variable ITS2 region both direct optimization and a priori alignment procedures gave similar topologies, although the former approach provided a better overall resolution. In terms of partitioned Bremer support (PBS), mitochondrial protein coding genes provided the largest contribution (83%) to the total tree resolution. This estimate however, appears to be partly spurious, due to the concerted inheritance of mitochondrial characters and probable cases of introgression or ancestral polymorphism. © The Willi Hennig Society 2005.     

Multi‐locus phylogenetics,lineage sorting,and reticulation in Pinus subsection Australes

Premise of the Study

Both incomplete lineage sorting and reticulation have been proposed as causes of phylogenetic incongruence. Disentangling these factors may be most difficult in long‐lived, wind‐pollinated plants with large population sizes and weak reproductive barriers.

Methods

We used solution hybridization for targeted enrichment and massive parallel sequencing to characterize low‐copy‐number nuclear genes and high‐copy‐number plastomes (Hyb‐Seq) in 74 individuals of Pinus subsection Australes, a group of ~30 New World pine species of exceptional ecological and economic importance. We inferred relationships using methods that account for both incomplete lineage sorting and reticulation.

Key Results

Concatenation‐ and coalescent‐based trees inferred from nuclear genes mainly agreed with one another, but they contradicted the plastid DNA tree in recovering the Attenuatae (the California closed‐cone pines) and Oocarpae (the egg‐cone pines of Mexico and Central America) as monophyletic and the Australes sensu stricto (the southern yellow pines) as paraphyletic to the Oocarpae. The plastid tree featured some relationships that were discordant with morphological and geographic evidence and species limits. Incorporating gene flow into the coalescent analyses better fit the data, but evidence supporting the hypothesis that hybridization explains the non‐monophyly of the Attenuatae in the plastid tree was equivocal.

Conclusions

Our analyses document cytonuclear discordance in Pinus subsection Australes. We attribute this discordance to ancient and recent introgression and present a phylogenetic hypothesis in which mostly hierarchical relationships are overlain by gene flow.     

Phylogenetic relationships in the Sceloporus variabilis (Squamata: Phrynosomatidae) complex based on three molecular markers,continuous characters and geometric morphometric data

The monophyly of the Sceloporus variabilis group is well established with five species and two species complexes, but phylogenetic relationships within species complexes are still uncertain. We studied 278 specimens in 20 terminals to sample all taxa in the “variabilis group,” including three subspecies in the “variabilis complex,” and two outgroups (Sceloporus grammicus and Sceloporus megalepidurus). We assembled an extensive morphological data set with discrete and continuous characters (distances and scale counts), including geometric morphometric data (landmark coordinates of three shapes), and a three‐marker molecular data set as well (ND4, 12S and RAG1). We conducted parsimony and Bayesian phylogenetic inferences on these data, including several partitioning and weighting schemes. We suggest elevating three subspecies to full species status. Therefore, we recommend recognition of nine species in the “variabilis group.” First, S. variabilis is sister to Sceloporus teapensis. In turn, Sceloporus cozumelae is sister to Sceloporus olloporus. These four species are a monophyletic group, which is sister to Sceloporus smithi. Finally, Sceloporus marmoratus is sister of the clade of five species. The other species in the “variabilis group” (Sceloporus chrysostictus, Sceloporus couchii and Sceloporus parvus) are a paraphyletic grade at the base of the tree. Our analyses reject the existence of the “variabilis complex.” We conducted a parsimony‐based ancestral reconstruction on body size (snout–vent length), femoral pores and dorsal scales and related morphological changes to geographic distribution of the species. Our phylogenetic hypothesis will allow best designs of comparative studies with species in the “variabilis group,” one of the earliest divergent lineages in the genus.     

The role of mitochondrial introgression in illuminating the evolutionary history of Nearctic treefrogs

Inferring the evolutionary history of lineages often becomes difficult when gene histories are in conflict with each other. Introgression, for example, can cause DNA sequences from one species to be more similar to sequences of a different species and lead to incongruence amongst gene trees. However, incorporating congruent and incongruent locus‐specific phylogenetic estimates with the geographical distribution of lineages may provide valuable insight into evolutionary processes important to speciation. In this study, we investigated mitochondrial introgression within the Hyla eximia group to better understand its role in illuminating the evolutionary history and phylogeography of these treefrogs. We reconstructed and compared the matrilineal history of the Hyla eximia group with estimates of evolutionary history inferred from nuclear genes. We tested for introgression within the mitochondrial and nuclear genes using a posterior predictive checking approach. Reconstructions of the species tree based on the mitochondrial DNA (mtDNA) and nuclear DNA data were strongly discordant. Introgression between lineages was widespread in the mtDNA data set (145 occurrences amongst 11 of the 16 lineages), but uncommon in the nuclear genes (12 occurrences amongst four of the 16 lineages). Nonetheless, the geographical structuring of mtDNA within species provides valuable information on biogeographical areas, ancient areas of hybridization, and unique histories of lineages within the H. eximia group. These results suggest that the combination of nuclear, mitochondrial, and spatial information can provide a more complete picture of ‘how evolutionary history played out’, particularly in cases where mitochondrial introgression is known to occur. © 2014 The Linnean Society of London     

Testing for shared biogeographic history in the lower Central American freshwater fish assemblage using comparative phylogeography: concerted,independent, or multiple evolutionary responses?

A central goal of comparative phylogeography is determining whether codistributed species experienced (1) concerted evolutionary responses to past geological and climatic events, indicated by congruent spatial and temporal patterns (“concerted‐response hypothesis”); (2) independent responses, indicated by spatial incongruence (“independent‐response hypothesis”); or (3) multiple responses (“multiple‐response hypothesis”), indicated by spatial congruence but temporal incongruence (“pseudocongruence”) or spatial and temporal incongruence (“pseudoincongruence”). We tested these competing hypotheses using DNA sequence data from three livebearing fish species codistributed in the Nicaraguan depression of Central America (Alfaro cultratus, Poecilia gillii, and Xenophallus umbratilis) that we predicted might display congruent responses due to co‐occurrence in identical freshwater drainages. Spatial analyses recovered different subdivisions of genetic structure for each species, despite shared finer‐scale breaks in northwestern Costa Rica (also supported by phylogenetic results). Isolation‐with‐migration models estimated incongruent timelines of among‐region divergences, with A. cultratus and Xenophallus populations diverging over Miocene–mid‐Pleistocene while P. gillii populations diverged over mid‐late Pleistocene. Approximate Bayesian computation also lent substantial support to multiple discrete divergences over a model of simultaneous divergence across shared spatial breaks (e.g., Bayes factor [B₁₀] = 4.303 for Ψ [no. of divergences] > 1 vs. Ψ = 1). Thus, the data support phylogeographic pseudoincongruence consistent with the multiple‐response hypothesis. Model comparisons also indicated incongruence in historical demography, for example, support for intraspecific late Pleistocene population growth was unique to P. gillii, despite evidence for finer‐scale population expansions in the other taxa. Empirical tests for phylogeographic congruence indicate that multiple evolutionary responses to historical events have shaped the population structure of freshwater species codistributed within the complex landscapes in/around the Nicaraguan depression. Recent community assembly through different routes (i.e., different past distributions or colonization routes), and intrinsic ecological differences among species, has likely contributed to the unique phylogeographical patterns displayed by these Neotropical fishes.     

Molecular phylogeny of Osmanthus (Oleaceae) based on non‐coding chloroplast and nuclear ribosomal internal transcribed spacer regions

Abstract The phylogenetic relationships of Osmanthus Lour. were investigated using the nuclear ribosomal internal transcribed spacer (ITS) regions and non‐coding chloroplast regions (psbA‐trnH, trnL‐F). The two datasets support the conclusion that Osmanthus is polyphyletic, with some species of the subtribe Oleinae nested within Osmanthus. Osmanthus didymopetalus P. S. Green is nested within the clade formed by species of section Osmanthus in two trees. Osmanthus attenuatus P. S. Green, O. yunnanensis P. S. Green, and O. gracilinervis R. L. Lu of traditional section Osmanthus are clearly divergent from other accessions, and do not form a monophyletic group with other Osmanthus accessions. Osmanthus marginatus Hemsl. is embedded in the clade formed by species of section Osmanthus in the ITS tree. In cpDNA trees all species of section Osmanthus are placed in the large clade and all species of section Leiolea formed a group. The taxonomic incongruence among trees for ITS and cpDNA indicate hybridization, as introgression may have occurred among some species of sections Osmanthus and Leiolea. Phylogeny of Osmanthus is discussed in light of molecular and morphological data, and a revised infrageneric classification with three sections (Leiolea, Siphosmanthu, and Osmanthus) is presented. The section Linocieroides is abandoned and united with section Osmanthus.     

Incongruence among mitochondrial,chloroplast and nuclear gene trees in <Emphasis Type="Italic">Pinus</Emphasis> subgenus <Emphasis Type="Italic">Strobus</Emphasis> (Pinaceae)

Introgression has been considered to be one of main factors leading to phylogenetic incongruence among different datasets at lower taxonomic levels. In the plants of Pinaceae, the mtDNA, cpDNA, and nuclear DNA (nrDNA) may have different evolutionary histories through introgression because they are inherited maternally, paternally and biparentally, respectively. We compared mtDNA, cpDNA, and two low-copy nrDNA phylogenetic trees in the genus Pinus subgenus Strobus, in order to detect unknown past introgression events in this group. nrDNA trees were mostly congruent with the cpDNA tree, and supported the recent sectional and subsectional classification system. In contrast, mtDNA trees split the members of sect. Quinquefoliae into two groups that were not observed in the other gene trees. The factors constituting incongruence may be divided into the following two categories: the different splits within subsect. Strobus, and the non-monophyly of subsect. Gerardianae. The former was hypothesized to have been caused by the past introgression of cpDNA, mtDNA or both between Eurasian and North American species through Beringia. The latter was likely caused by the chimeric structure of the mtDNA sequence of P. bungeana, which might have originated through past hybridization, or through a horizontal transfer event and subsequent recombination. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

A tale of worldwide success: Behind the scenes of Carex (Cyperaceae) biogeography and diversification

The megadiverse genus Carex (c. 2000 species, Cyperaceae) has a nearly cosmopolitan distribution, displaying an inverted latitudinal richness gradient with higher species diversity in cold‐temperate areas of the Northern Hemisphere. Despite great expansion in our knowledge of the phylogenetic history of the genus and many molecular studies focusing on the biogeography of particular groups during the last few decades, a global analysis of Carex biogeography and diversification is still lacking. For this purpose, we built the hitherto most comprehensive Carex‐dated phylogeny based on three markers (ETS–ITS–matK), using a previous phylogenomic Hyb‐Seq framework, and a sampling of two‐thirds of its species and all recognized sections. Ancestral area reconstruction, biogeographic stochastic mapping, and diversification rate analyses were conducted to elucidate macroevolutionary biogeographic and diversification patterns. Our results reveal that Carex originated in the late Eocene in E Asia, where it probably remained until the synchronous diversification of its main subgeneric lineages during the late Oligocene. E Asia is supported as the cradle of Carex diversification, as well as a “museum” of extant species diversity. Subsequent “out‐of‐Asia” colonization patterns feature multiple asymmetric dispersals clustered toward present times among the Northern Hemisphere regions, with major regions acting both as source and sink (especially Asia and North America), as well as several independent colonization events of the Southern Hemisphere. We detected 13 notable diversification rate shifts during the last 10 My, including remarkable radiations in North America and New Zealand, which occurred concurrently with the late Neogene global cooling, which suggests that diversification involved the colonization of new areas and expansion into novel areas of niche space.     

Evolutionary history of Trachylepis skinks in the Seychelles islands: introgressive hybridization,morphological evolution and geographic structure

The Seychelles is a remarkably interesting archipelago for evolutionary studies, but only recently have molecular markers been used to explore its biogeographic patterns. Here we used morphological and molecular data to examine diversity and phylogenetic relationships of two endemic skink sister‐species from this archipelago: Trachylepis sechellensis and Trachylepis wrightii. Mitochondrial DNA genealogy rendered a monophyletic T. wrightii nested within a paraphyletic T. sechellensis, whereas nuclear DNA sequences from five unlinked markers reflected the accepted taxonomy. Hybridization and massive mtDNA introgression leading to the complete replacement of the native mtDNA lineage of T. sechellensis in some of the islands were invoked to explain this result, and morphological variation also seemed to reflect this pattern of reticulation. A Mio‐Pliocene divergence between both species is suggested. Multilocus molecular data were used to uncover biogeographic patterns within the archipelago, which reflected shared patterns with other co‐distributed lizard taxa; specifically a north–south marked structure, a close relationship between populations from Fregate and the southern islands, and a detectable isolation within the southern group, between Mahé, and Silhouette and North Islands. Gene flow from these latter islands towards the northern group was also suggested. These results add to the growing body of evidence of the influence of geographic distance and sea‐level oscillations in shaping the genetic structure of Seychellois taxa and of the existence of common biogeographic patterns across the archipelago.     

DNA barcoding supports reclassification of Japanese Chironomus species (Diptera: Chironomidae)

Non‐biting midges (Diptera: Chironomidae) adapt to species‐specific environmental conditions and hence are promising bioindicators for aquatic and ecotoxicological monitoring. Although their utility for these purposes was historically limited by difficulties in their morphological identification, DNA barcoding offers a possible solution. Here, eight Japanese species of the genus Chironomus, which is characterized by its worldwide distribution and abundance among Chironomidae, were subjected to DNA barcoding using cytochromec oxidase subunit I (COI). To examine whether this DNA barcode is a useful indicator for Japanese species of Chironomus, we calculated genetic distances within and between the COI sequences of Chironomus species both from this study and worldwide and constructed phylogenetic trees. Based on 415 bp COI sequences and the Kimura two‐parameter model, the average genetic distances within 37 species and between 72 species were 2.6% and 17.2%, respectively. Although the ranges of genetic distances within and between species overlapped from 0.8% to 17.3%, 99.7% of average genetic distances between species were >3.0%. Some of this overlap is attributable to distances within species that were “too large” as well as those between species that were “too small”. Of eight Japanese species examined, two showed genetic distances between species that were below a 3.0% threshold, and four had distances within species that were greater than 3.0%. These results suggest a possible reclassification of these species and the need for further sampling to unveil biogeographic variations among different countries and regions.     

Sources of character conflict in a clade of water striders (Heteroptera: Gerridae)

Incongruence among trees reconstructed with different data may stem from historical (gene tree‐species tree conflict) or process (character change biases) phenomena. Regardless of the source, incongruent data, as determined with “global” measures of homoplasy, have often been excluded from parsimony analysis of the combined data. Recent studies suggest that these homoplasy measures do not predict the contribution of each character to overall tree structure. Branch support measures identify, on a character to node basis, sources of support and conflict resulting from a simultaneous analysis of the data. We implement these branch support measures to identify sources of character conflict in a clade of water striders consisting of Gerris Fabricius, Aquarius Schellenberg, and Limnoporus Stål species. Separate analyses of morphology, mitochondrial cytochrome oxidase I (COI), large mitochondrial ribosomal subunit (16SrRNA), and elongation factor‐1α (EF‐1α) data resulted in cladograms that varied in resolution and topological concordance. Simultaneous analysis of the data resulted in two trees that were unresolved for one node in a strict consensus. The topology agreed with current classification except for the placements of Aquarius chilensis and the Aquarius remigis species group closer to Gerris than to congeneric species. Branch support measures indicated that support derived from each data set varied among nodes, but COI had an overall negative effect on branch support. However, Spearman rank correlation of partitioned branch support values indicated no negative associations of branch support between any data sets and a positive association between EF‐1α and 16SrRNA. Thus incongruence among data sets was not drastic and the gene‐tree versus species tree phenomenon was not implicated. Biases in character change were a more likely reason for incongruence, although saturation curves and incongruence length difference for COI indicated little potential for homoplasy. However, a posteriori inspection of COI nucleotide change with reference to the simultaneous analysis tree revealed AT and codon biases. These biases were not associated with branch support measures. Therefore, it is difficult to predict incongruence or identify its cause. Exclusion of data is ill advised because every character is potentially parsimony informative.     

Phylogenetic relationships of toads of the Rhinella granulosa group (Anura: Bufonidae): a molecular perspective with comments on hybridization and introgression

The Rhinella granulosa group consists of 13 species of toads distributed throughout open areas of South America and Panama. In this paper we perform a phylogenetic analysis considering all but one species of the group, employing five nuclear and four mitochondrial genes, for up to 7910 bp per specimen. Separate phylogenetic analyses under direct optimization (DO) of nuclear and mitochondrial sequences recovered the R. granulosa group as monophyletic and revealed topological incongruence that can be explained mainly by multiple events of hybridization and introgression, both mitochondrial and nuclear. The DO combined analysis, after the exclusion of putatively introgressed or heterozygous genomes, resulted in a phylogenetic hypothesis for the R. granulosa group in which most of the species are recovered as monophyletic, but with interspecific relationships poorly supported. The optimization of morphological (adult and larval), chromosomal, and behavioural characters resulted in 12 putative phenotypic synapomorphies for this species group and some other synapomorphies for internal clades. Our results indicate the need for additional population genetic studies on R. dorbignyi and R. fernandezae to corroborate the taxonomic status of both taxa. Finally, we discuss biological and genetic characteristics of Bufonidae, as possible explanations for the common occurrence of hybridization and introgression observed in some lineages of this family.     

Chloroplast DNA inheritance and variation in <Emphasis Type="Italic">Leucadendron</Emphasis> species (Proteaceae) as revealed by PCR-RFLP

The inheritance of chloroplast DNA (cpDNA) in Leucadendron species was studied by polymerase chain reaction–restriction fragment length polymorphism (PCR–RFLP) analysis. A total of 100 progeny from five interspecific crosses involving seven parental species were tested, and all progeny exhibited the cpDNA restriction fragment pattern of the female parent, indicating that cpDNA in Leucadendron is maternally inherited. PCR-RFLP was also employed to study cpDNA variation among 21 Leucadendron species. Parsimony analysis using a heuristic search resulted in a phylogenetic tree that showed limited agreement to the taxonomic classification of Leucadendron species, based on morphological characteristics. The incongruence between cpDNA phylogenetic and taxonomic groupings in Leucadendron may be due to reticulate evolution involving a combination of hybridization and introgression, convergent evolution and/or lineage sorting at the interspecific, intersubsectional and intersectional levels.     

Testing hypotheses of mitochondrial gene‐tree paraphyly: unravelling mitochondrial capture of the Streak‐breasted Scimitar Babbler (Pomatorhinus ruficollis) by the Taiwan Scimitar Babbler (Pomatorhinus musicus)

Species‐level paraphyly inferred from mitochondrial gene trees is a prevalent phenomenon in taxonomy and systematics, but there are several potential causes that are not easily explained by currently used methods. This study aimed to test the underlying causes behind the observed paraphyly of Streak‐breasted Scimitar Babbler (Pomatorhinus ruficollis) via statistical analyses of four mitochondrial (mtDNA) and nine nuclear (nuDNA) genes. Mitochondrial gene trees show paraphyly of P. ruficollis with respect to the Taiwan Scimitar Babbler (Pomatorhinus musicus), but nuclear genealogies support a sister‐group relationship. Predictive coalescent simulations imply several hypothetical explanations, the most likely being mitochondrial capture of P. ruficollis by P. musicus for the observed cyto‐nuclear incongruence. Further approximate Bayesian computation suggests a unidirectional introgression model with substantial level of gene flow from P. ruficollis to P. musicus during their initial divergence during the Late Pleistocene. This specific observation frames several potential causes for incongruent outcomes of mitochondrial and nuclear introgression in general, and on the whole, our results underscore the strength of multiple independent loci for species delimitation and importance of testing hypotheses that explain disparate causes of mitochondrial gene‐tree paraphyly.     

Morphological Identities of Two Different Marine Stramenopile Environmental Sequence Clades: Bicosoeca kenaiensis (Hilliard, 1971) and Cantina marsupialis (Larsen and Patterson, 1990) gen. nov., comb. nov.

Although environmental DNA surveys improve our understanding of biodiversity, interpretation of unidentified lineages is limited by the absence of associated morphological traits and living cultures. Unidentified lineages of marine stramenopiles are called “MAST clades”. Twenty‐five MAST clades have been recognized: MAST‐1 through MAST‐25; seven of these have been subsequently discarded because the sequences representing those clades were found to either (1) be chimeric or (2) affiliate within previously described taxonomic groups. Eighteen MAST clades remain without a cellular identity. Moreover, the discarded “MAST‐13” has been used in different studies to refer to two different environmental sequence clades. After establishing four cultures representing two different species of heterotrophic stramenopiles and then characterizing their morphology and molecular phylogenetic positions, we determined that the two different species represented the two different MAST‐13 clades: (1) a lorica‐bearing Bicosoeca kenaiensis and (2) a microaerophilic flagellate previously named “Cafeteria marsupialis”. Both species were previously described with only light microscopy; no cultures, ultrastructural data or DNA sequences were available from these species prior to this study. The molecular phylogenetic position of three different “C. marsupialis” isolates was not closely related to the type species of Cafeteria; therefore, we established a new genus for these isolates, Cantina gen. nov.     

TESTING MITOCHONDRIAL CAPTURE AND DEEP COALESCENCE IN AMAZONIAN CICHLID FISHES (CICHLIDAE: CICHLA)

Hybridization and introgression have important consequences in evolution, such as increasing the genetic diversity and adaptive potential of a species. One of their most conspicuous footprints is discordance among gene trees or between genes and phenotypes. However, most studies that report introgression fail to disprove the null hypothesis that genetic incongruence may result from stochastic sorting of ancestral allelic polymorphisms. In the case of ancient introgression, these two processes may be especially difficult to distinguish topologically, but they make different predictions about the patterns of coalescence among loci. Here we apply three methods, molecular dating, multispecies coalescent models, and gene tree simulation under coalescence, to compare these two hypotheses that explain the polyphyletic mtDNA of the butterfly peacock bass, Cichla orinocensis. In comparison with a species tree based on 20 unlinked nuclear loci, we determined that mtDNA divergences were too recent to be explained by ancestral polymorphism. Similarly, coalescent species tree branches were significantly shorter when putative introgressed mtDNA was incorporated, and simulations showed the mtDNA topology to be unlikely under lineage sorting only. We conclude that introgression approximately 1.5 million years ago resulted in capture by C. orinocensis of an mtDNA lineage ancestral to the modern subspecies C. oc. monoculus.     

Diversification of Hawaiian Cyrtandra (Gesneriaceae) under the influence of incomplete lineage sorting and hybridization

Cyrtandra (Gesneriaceae) is a genus of flowering plants with over 800 species distributed throughout Southeast Asia and the Pacific Islands. On the Hawaiian Islands, 60 named species and over 89 putative hybrids exist, most of which are identified on the basis of morphology. Despite many previous studies on the Hawaiian Cyrtandra lineage, questions regarding the reconciliation of morphology and genetics remain, many of which can be attributed to the relatively young age and evidence of hybridization between species. We utilized targeted enrichment, high‐throughput sequencing, and modern phylogenomics tools to test 31 Hawaiian Cyrtandra samples (22 species, two putative hybrids, four species with two samples each, one species with four samples) and two outgroups for species relationships and hybridization in the presence of incomplete lineage sorting (ILS). Both concatenated and species‐tree methods were used to reconstruct species relationships, and network analyses were conducted to test for hybridization. We expected to see high levels of ILS and putative hybrids intermediate to their parent species. Phylogenies reconstructed from the concatenated and species‐tree methods were highly incongruent, most likely due to high levels of incomplete lineage sorting. Network analyses inferred gene flow within this lineage, but not always between taxa that we expected. Multiple hybridizations were inferred, but many were on deeper branches of the island lineages suggesting a long history of hybridization. We demonstrated the utility of high‐throughput sequencing and a phylogenomic approach using 569 loci to understanding species relationships and gene flow in the presence of ILS.