Chromosome-Aware Phylogenomics of Assassin Bugs (Hemiptera: Reduvioidea) Elucidates Ancient Gene Conflict.

Though the phylogenetic signal of loci on sex chromosomes can differ from those on autosomes, chromosomal-level genome assemblies for nonvertebrates are still relatively scarce and conservation of chromosomal gene content across deep phylogenetic scales has therefore remained largely unexplored. We here assemble a uniquely large and diverse set of samples (17 anchored hybrid enrichment, 24 RNA-seq, and 70 whole-genome sequencing samples of variable depth) for the medically important assassin bugs (Reduvioidea). We assess the performance of genes based on multiple features (e.g., nucleotide vs. amino acid, nuclear vs. mitochondrial, and autosomal vs. X chromosomal) and employ different methods (concatenation and coalescence analyses) to reconstruct the unresolved phylogeny of this diverse (∼7,000 spp.) and old (>180 Ma) group. Our results show that genes on the X chromosome are more likely to have discordant phylogenies than those on autosomes. We find that the X chromosome conflict is driven by high gene substitution rates that impact the accuracy of phylogenetic inference. However, gene tree clustering showed strong conflict even after discounting variable third codon positions. Alternative topologies were not particularly enriched for sex chromosome loci, but spread across the genome. We conclude that binning genes to autosomal or sex chromosomes may result in a more accurate picture of the complex evolutionary history of a clade.

Phylogenomics and biogeography of arid-adapted Chlamydogobius goby fishes.

The progressive aridification of the Australian continent from âˆ¼ 20 million years ago posed severe challenges for the persistence of its resident biota. A key question involves the role of refugial habitats - specifically, their ability to mediate the effects of habitat loss and fragmentation, and their potential to shape opportunities for allopatric speciation. With freshwater species, for example, the patchiness, or absence, of water will constrain distributions. However, aridity may not necessarily isolate populations if disjunct refugia experience frequent hydrological connections. To investigate this potential dichotomy, we explored the evolutionary history of the Chlamydogobius gobies (Gobiiformes: Gobiidae), an arid-adapted genus of six small, benthic fish species that exploit all types of waterbodies (i.e. desert springs, waterholes and bore-fed wetlands, coastal estuarine creeks and mangroves) across parts of central and northern Australia. We used Anchored Phylogenomics to generate a highly resolved phylogeny of the group from sequence data for 260 nuclear loci. Buttressed by companion allozyme and mtDNA datasets, our molecular findings infer the diversification of Chlamydogobius in arid Australia, and provide a phylogenetic structure that cannot be simply explained by invoking allopatric speciation events reflecting current geographic proximity. Our findings are generally consistent with the existing morphological delimitation of species, with one exception: at the shallowest nodes of phylogenetic reconstruction, the molecular data do not fully support the current dichotomous delineation of C. japalpa from C. eremius in Kati Thanda-Lake Eyre-associated waterbodies. Together these findings illustrate the ability of structural (hydrological) connections to generate patterns of connectivity and isolation for an ecologically moderate disperser in response to ongoing habitat aridification. Finally, we explore the implications of these results for the immediate management of threatened (C. gloveri) and critically endangered (C. micropterus, C. squamigenus) congeners.

Phylogenomics of Anguis and Pseudopus (Squamata, Anguidae) indicates Balkan-Apennine mitochondrial capture associated with the Messinian event.

A dated phylogenetic hypothesis on the evolutionary history of the extant taxa of the Western Palearctic lizards Anguis and Pseudopus is revised using genome-wide nuclear DNA and mitogenomes. We found overall concordance between nuclear and mitochondrial DNA phylogenies, with one significant exception - the Apennine A. veronensis. In mitochondrial DNA, this species forms a common clade with the earliest diverging lineage, the southern Balkan endemic A. cephallonica, while it clusters together with A. fragilis in nuclear DNA. The nuclear phylogeny conforms to the morphology, which is relatively similar between A. veronensis and A. fragilis. The most plausible explanation for the mitonuclear discordance is ancient mitochondrial capture from the Balkan ancestor of A. cephallonica to the Apennine population of the A. fragilis-veronensis ancestor. We hypothesize that this capture occurred only in a geographically restricted population. The dating of this presumed mitochondrial introgression and capture coincides with the Messinian event, when the Balkan and Apennine Peninsulas were presumably largely connected. The dated nuclear phylogenomic reconstruction estimated the divergence of A. cephallonica around 12 Mya, while the sister clade representing the A. fragilis species complex consisting of the sister species A. fragilis-A. veronensis and A. colchica-A. graeca further diversified around 7 Mya. The depth of nuclear divergence among the evolutionary lineages of Pseudopus (0.5-1.2 Mya) supports their subspecies status.

Sympatric and independently evolving lineages in the Thoropa miliaris - T. taophora species complex (Anura: Cycloramphidae).

Species delimitation can be challenging and affected by subjectivity. Sibling lineages that occur in sympatry constitute good candidates for species delimitation regardless of the adopted species concept. The Thoropa miliaris + T. taophora species complex exhibits high genetic diversity distributed in several lineages that occur sympatrically in the southeastern Atlantic Forest of Brazil. We used 414 loci obtained by anchored hybrid enrichment to characterize genetic variation in the Thoropa miliaris species group (T. saxatilis, T megatympanum, T. miliaris, and T. taophora), combining assignment analyses with traditional and coalescent phylogeny reconstruction. We also investigated evolutionary independence in co-occurring lineages by estimating gene flow, and validated lineages under the multispecies coalescent. We recovered most previously described lineages as unique populations in assignment analyses; exceptions include two lineages within T. miliaris that are further substructured, and the merging of all T. taophora lineages. We found very low probabilities of gene flow between sympatric lineages, suggesting independent evolution. Species tree inferences and species delimitation yielded resolved relationships and indicate that all lineages constitute putative species that diverged during the Pliocene and Pleistocene, later than previously estimated.

Phylogenomics of African radiation of Praomyini (Muridae: Murinae) rodents: First fully resolved phylogeny, evolutionary history and delimitation of extant genera.

The tribe Praomyini is a diversified group including 64 species and eight extant rodent genera. They live in a broad spectrum of habitats across whole sub-Saharan Africa. Members of this tribe are often very abundant, they have a key ecological role in ecosystems, they are hosts of many potentially pathogenic microorganisms and comprise numerous agricultural pests. Although this tribe is well supported by both molecular and morphological data, its intergeneric relationships and the species contents of several genera are not yet fully resolved. Recent molecular data suggest that at least three genera in current sense are paraphyletic. However, in these studies the species sampling was sparse and the resolution of relationships among genera was poor, probably due to a fast radiation of the tribe dated to the Miocene and insufficient amount of genetic data. Here we used genomic scale data (395 nuclear loci = 610,965 bp long alignment and mitogenomes = 14,745 bp) and produced the first fully resolved species tree containing most major lineages of the Praomyini tribe (i.e. all but one currently delimited genera and major intrageneric clades). Results of a fossil-based divergence dating analysis suggest that the radiation started during the Messinian stage (ca. 7 Ma) and was likely linked to a fragmentation of the pan-African Miocene forest. Some lineages remained in the rain forests, while many others adapted to a broad spectrum of new open lowland and montane habitats that appeared at the beginning of Pliocene. Our analyses clearly confirmed the presence of three polyphyletic genera (Praomys, Myomyscus and Mastomys). We review current knowledge of these three genera and suggest corresponding taxonomic changes. To keep genera monophyletic, we propose taxonomic re-arrangements and delimit four new genera. Furthermore, we discovered a new highly divergent genetic lineage of Praomyini in southwestern Ethiopia, which is described as a new species and genus.

Phylogeny and classification of Odonata using targeted genomics.

Dragonflies and damselflies are a charismatic, medium-sized insect order (~6300 species) with a unique potential to approach comparative research questions. Their taxonomy and many ecological traits for a large fraction of extant species are relatively well understood. However, until now, the lack of a large-scale phylogeny based on high throughput data with the potential to connect both perspectives has precluded comparative evolutionary questions for these insects. Here, we provide an ordinal hypothesis of classification based on anchored hybrid enrichment using a total of 136 species representing 46 of the 48 families or incertae sedis, and a total of 478 target loci. Our analyses recovered the monophyly for all three suborders: Anisoptera, Anisozygoptera and Zygoptera. Although the backbone of the topology was reinforced and showed the highest support values to date, our genomic data was unable to stronglyresolve portions of the topology. In addition, a quartet sampling approach highlights the potential evolutionary scenarios that may have shaped evolutionary phylogeny (e.g., incomplete lineage sorting and introgression) of this taxon. Finally, in light of our phylogenomic reconstruction and previous morphological and molecular information we proposed an updated odonate classification and define five new families (Amanipodagrionidae fam. nov., Mesagrionidae fam. nov., Mesopodagrionidae fam. nov., Priscagrionidae fam. nov., Protolestidae fam. nov.) and reinstate another two (Rhipidolestidae stat. res., Tatocnemididae stat. res.). Additionally, we feature the problematic taxonomic groupings for examination in future studies to improve our current phylogenetic hypothesis.

First worldwide molecular phylogeny of the morphologically and ecologically hyperdiversified snapping shrimp genus Alpheus (Malacostraca: Decapoda).

Hyperdiverse animal groups raise intriguing questions regarding the factors that generate and maintain their diversity. The snapping shrimp genus Alpheus (with >300 described species) is a spectacularly diversified group of decapod crustaceans that serves as an exemplary system for addressing evolutionary questions regarding morphological adaptations, symbiosis, cryptic diversity and molecular divergence. A lack of information regarding evolutionary relationships among species has limited investigations into the mechanisms that drive the diversification of Alpheus. Previous phylogenetic studies of Alpheus have been restricted in scope, while molecular datasets used for phylogenetic reconstructions have been based solely on mitochondrial and a handful of nuclear markers. Here we use an anchored hybrid enrichment (AHE) approach to resolve phylogenetic relationships among species of Alpheus. The AHE method generated sequence data for 240 loci (>72,000 bp) for 65 terminal species that span the geographic, ecological and taxonomic diversity of Alpheus. Our resulting, well-supported phylogeny demonstrates a lack of monophyly for five out of seven morphologically defined species groups that have traditionally been used as a framework in Alpheus taxonomy. Our results also suggest that symbiotic associations with a variety of other animals have evolved independently in at least seven lineages in this genus. Our AHE phylogeny represents the most comprehensive phylogenetic treatment of Alpheus to date and will provide a useful evolutionary framework to further investigate questions, such as various modifications of the snapping claw and the role of habitat specialization and symbiosis in promoting speciation. Running head: PHYLOGENY OF THE SNAPPING SHRIMP GENUS ALPHEUS.

Nuclear phylogenomics, but not mitogenomics, resolves the most successful Late Miocene radiation of African mammals (Rodentia: Muridae: Arvicanthini).

The tribe Arvicanthini (Muridae: Murinae) is a highly diversified group of rodents (ca. 100 species) and with 18 African genera (plus one Asiatic) represents probably the most successful adaptive radiation of extant mammals in Africa. They colonized a broad spectrum of habitats (from rainforests to semi-deserts) in whole sub-Saharan Africa and their members often belong to most abundant parts of mammal communities. Despite intensive efforts, the phylogenetic relationships among major lineages (i.e. genera) remained obscured, which was likely caused by the intensive radiation of the group, dated to the Late Miocene. Here we used genomic scale data (377 nuclear loci; 581,030 bp) and produced the first fully resolved species tree containing all currently delimited genera of the tribe. Mitogenomes were also extracted, and while the results were largely congruent, there was less resolution at basal nodes of the mitochondrial phylogeny. Results of a fossil-based divergence dating analysis suggest that the African radiation started early after the colonization of Africa by a single arvicanthine ancestor from Asia during the Messinian stage (ca. 7 Ma), and was likely linked with a fragmentation of the pan-African Miocene forest. Some lineages remained in the rain forest, while many others successfully colonized broad spectrum of new open habitats (e.g. savannas, wetlands or montane moorlands) that appeared at the beginning of Pliocene. One lineage even evolved partially arboricolous life style in savanna woodlands, which allowed them to re-colonize equatorial forests. We also discuss delimitation of genera in Arvicanthini and propose corresponding taxonomic changes.

Phylogenomics of Ichneumonoidea (Hymenoptera) and implications for evolution of mode of parasitism and viral endogenization.

Ichneumonoidea is one of the most diverse lineages of animals on the planet with >48,000 described species and many more undescribed. Parasitoid wasps of this superfamily are mostly beneficial insects that attack and kill other arthropods and are important for understanding diversification and the evolution of life history strategies related to parasitoidism. Further, some lineages of parasitoids within Ichneumonoidea have acquired endogenous virus elements (EVEs) that are permanently a part of the wasp's genome and benefit the wasp through host immune disruption and behavioral control. Unfortunately, understanding the evolution of viral acquisition, parasitism strategies, diversification, and host immune disruption mechanisms, is deeply limited by the lack of a robust phylogenetic framework for Ichneumonoidea. Here we design probes targeting 541 genes across 91 taxa to test phylogenetic relationships, the evolution of parasitoid strategies, and the utility of probes to capture polydnavirus genes across a diverse array of taxa. Phylogenetic relationships among Ichneumonoidea were largely well resolved with most higher-level relationships maximally supported. We noted codon use biases between the outgroups, Braconidae, and Ichneumonidae and within Pimplinae, which were largely solved through analyses of amino acids rather than nucleotide data. These biases may impact phylogenetic reconstruction and caution for outgroup selection is recommended. Ancestral state reconstructions were variable for Braconidae across analyses, but consistent for reconstruction of idiobiosis/koinobiosis in Ichneumonidae. The data suggest many transitions between parasitoid life history traits across the whole superfamily. The two subfamilies within Ichneumonidae that have polydnaviruses are supported as distantly related, providing strong evidence for two independent acquisitions of ichnoviruses. Polydnavirus capture using our designed probes was only partially successful and suggests that more targeted approaches would be needed for this strategy to be effective for surveying taxa for these viral genes. In total, these data provide a robust framework for the evolution of Ichneumonoidea.

Drift, selection and adaptive variation in small populations of a threatened rattlesnake.

An important goal of conservation genetics is to determine if the viability of small populations is reduced by a loss of adaptive variation due to genetic drift. Here, we assessed the impact of drift and selection on direct measures of adaptive variation (toxin loci encoding venom proteins) in the eastern massasauga rattlesnake (Sistrurus catenatus), a threatened reptile that exists in small isolated populations. We estimated levels of individual polymorphism in 46 toxin loci and 1,467 control loci across 12 populations of this species, and compared the results with patterns of selection on the same loci following speciation of S. catenatus and its closest relative, the western massasauga (S. tergeminus). Multiple lines of evidence suggest that both drift and selection have had observable impacts on standing adaptive variation. In support of drift effects, we found little evidence for selection on toxin variation within populations and a significant positive relationship between current levels of adaptive variation and long- and short-term estimates of effective population size. However, we also observed levels of directional selection on toxin loci among populations that are broadly similar to patterns predicted from interspecific selection analyses that pre-date the effects of recent drift, and that functional variation in these loci persists despite small short-term effective sizes. This suggests that much of the adaptive variation present in populations may represent an example of "drift debt," a nonequilibrium state where present-day levels of variation overestimate the amount of functional genetic diversity present in future populations.

Polyploidy breaks speciation barriers in Australian burrowing frogs Neobatrachus.

Polyploidy has played an important role in evolution across the tree of life but it is still unclear how polyploid lineages may persist after their initial formation. While both common and well-studied in plants, polyploidy is rare in animals and generally less understood. The Australian burrowing frog genus Neobatrachus is comprised of six diploid and three polyploid species and offers a powerful animal polyploid model system. We generated exome-capture sequence data from 87 individuals representing all nine species of Neobatrachus to investigate species-level relationships, the origin and inheritance mode of polyploid species, and the population genomic effects of polyploidy on genus-wide demography. We describe rapid speciation of diploid Neobatrachus species and show that the three independently originated polyploid species have tetrasomic or mixed inheritance. We document higher genetic diversity in tetraploids, resulting from widespread gene flow between the tetraploids, asymmetric inter-ploidy gene flow directed from sympatric diploids to tetraploids, and isolation of diploid species from each other. We also constructed models of ecologically suitable areas for each species to investigate the impact of climate on differing ploidy levels. These models suggest substantial change in suitable areas compared to past climate, which correspond to population genomic estimates of demographic histories. We propose that Neobatrachus diploids may be suffering the early genomic impacts of climate-induced habitat loss, while tetraploids appear to be avoiding this fate, possibly due to widespread gene flow. Finally, we demonstrate that Neobatrachus is an attractive model to study the effects of ploidy on the evolution of adaptation in animals.

Phylogenomics, Biogeography, and Morphometrics Reveal Rapid Phenotypic Evolution in Pythons After Crossing Wallace's Line.

Ecological opportunities can be provided to organisms that cross stringent biogeographic barriers towards environments with new ecological niches. Wallace's and Lyddeker's lines are arguably the most famous biogeographic barriers, separating the Asian and Australo-Papuan biotas. One of the most ecomorphologically diverse groups of reptiles, the pythons, is distributed across these lines, and are remarkably more diverse in phenotype and ecology east of Lydekker's line in Australo-Papua. We used an anchored hybrid enrichment approach, with near complete taxon sampling, to extract mitochondrial genomes and 376 nuclear loci to resolve and date their phylogenetic history. Biogeographic reconstruction demonstrates that they originated in Asia around 38-45 Ma and then invaded Australo-Papua around 23 Ma. Australo-Papuan pythons display a sizeable expansion in morphological space, with shifts towards numerous new adaptive optima in head and body shape, coupled with the evolution of new micro-habitat preferences. We provide an updated taxonomy of pythons and our study also demonstrates how ecological opportunity following colonization of novel environments can promote morphological diversification in a formerly ecomorphologically conservative group. [Adaptive radiation; anchored hybrid enrichment; biogeography; morphometrics; snakes.].

Species delimitation and systematics of the green pythons (Morelia viridis complex) of melanesia and Australia.

Molecular data sets and the increasing use of integrative systematics is revealing cryptic diversity in a range of taxa - particularly in remote and poorly sampled landscapes like the island of New Guinea. Green pythons (Morelia viridis complex) are one of the most conspicuous elements of this island's fauna, with large numbers taken from the wild to supply international demand for exotic pets. We test hypotheses about species boundaries in green pythons from across New Guinea and Australia with mitochondrial genomes, 389 nuclear exons, and comprehensive assessment of morphological variation. Strong genetic structuring of green python populations and species delimitation methods confirm the presence of two species, broadly occurring north and south of New Guinea's central mountains. Our data also support three subspecies within the northern species. Subtle but consistent morphological divergence among the putative taxa is concordant with patterns of molecular divergence. Our extensive sampling identifies several zones of hitherto unknown biogeographical significance on the island of New Guinea. We revise the taxonomy of the group, discuss the relevance of our findings in the context of Papuan biogeography and the implications of our systematic changes for the conservation management of these taxa.

Environmental temperatures shape thermal physiology as well as diversification and genome-wide substitution rates in lizards.

Climatic conditions changing over time and space shape the evolution of organisms at multiple levels, including temperate lizards in the family Lacertidae. Here we reconstruct a dated phylogenetic tree of 262 lacertid species based on a supermatrix relying on novel phylogenomic datasets and fossil calibrations. Diversification of lacertids was accompanied by an increasing disparity among occupied bioclimatic niches, especially in the last 10 Ma, during a period of progressive global cooling. Temperate species also underwent a genome-wide slowdown in molecular substitution rates compared to tropical and desert-adapted lacertids. Evaporative water loss and preferred temperature are correlated with bioclimatic parameters, indicating physiological adaptations to climate. Tropical, but also some populations of cool-adapted species experience maximum temperatures close to their preferred temperatures. We hypothesize these species-specific physiological preferences may constitute a handicap to prevail under rapid global warming, and contribute to explaining local lizard extinctions in cool and humid climates.

Phylogenomic Analysis of a Putative Missing Link Sparks Reinterpretation of Leech Evolution.

Leeches (Hirudinida) comprise a charismatic, yet often maligned group of organisms. Despite their ecological, economic, and medical importance, a general consensus on the phylogenetic relationships of major hirudinidan lineages is lacking. This absence of a consistent, robust phylogeny of early-diverging lineages has hindered our understanding of the underlying processes that enabled evolutionary diversification of this clade. Here, we used an anchored hybrid enrichment-based phylogenomic approach, capturing hundreds of loci to investigate phylogenetic relationships among major hirudinidan lineages and their closest living relatives. Our results suggest that a dramatic reinterpretation of early leech evolution is warranted. We recovered Branchiobdellida as sister to a clade that includes all major lineages of hirudinidans, but found Acanthobdella to be nested within Oceanobdelliformes. These results cast doubt on the utility of Acanthobdella as a "missing link" used to explain the origin of blood-feeding in hirudineans. Further, our results support a deep divergence between predominantly marine and freshwater lineages, while not supporting the reciprocal monophyly of jawed and proboscis-bearing leeches. To sum up, our phylogenomic resolution of early-diverging leeches provides a necessary foundation for illuminating the evolution of host-symbiont associations and key adaptations that have allowed leeches to colonize a wide diversity of habitats worldwide.

Phylogenomic Evidence Overturns Current Conceptions of Social Evolution in Wasps (Vespidae).

The hypothesis that eusociality originated once in Vespidae has shaped interpretation of social evolution for decades and has driven the supposition that preimaginal morphophysiological differences between castes were absent at the outset of eusociality. Many researchers also consider casteless nest-sharing an antecedent to eusociality. Together, these ideas endorse a stepwise progression of social evolution in wasps (solitary → casteless nest-sharing → eusociality with rudimentary behavioral castes → eusociality with preimaginal caste-biasing (PCB) → morphologically differentiated castes). Here, we infer the phylogeny of Vespidae using sequence data generated via anchored hybrid enrichment from 378 loci across 136 vespid species and perform ancestral state reconstructions to test whether rudimentary and monomorphic castes characterized the initial stages of eusocial evolution. Our results reject the single origin of eusociality hypothesis, contest the supposition that eusociality emerged from a casteless nest-sharing ancestor, and suggest that eusociality in Polistinae + Vespinae began with castes having morphological differences. An abrupt appearance of castes with ontogenetically established morphophysiological differences conflicts with the current conception of stepwise social evolution and suggests that the climb up the ladder of sociality does not occur through sequential mutation. Phenotypic plasticity and standing genetic variation could explain how cooperative brood care evolved in concert with nest-sharing and how morphologically dissimilar castes arose without a rudimentary intermediate. Furthermore, PCB at the outset of eusociality implicates a subsocial route to eusociality in Polistinae + Vespinae, emphasizing the role of mother-daughter interactions and subfertility (i.e. the cost component of kin selection) in the origin of workers.

Phylogenomic Data Yield New and Robust Insights into the Phylogeny and Evolution of Weevils.

The phylogeny and evolution of weevils (the beetle superfamily Curculionoidea) has been extensively studied, but many relationships, especially in the large family Curculionidae (true weevils; > 50,000 species), remain uncertain. We used phylogenomic methods to obtain DNA sequences from 522 protein-coding genes for representatives of all families of weevils and all subfamilies of Curculionidae. Most of our phylogenomic results had strong statistical support, and the inferred relationships were generally congruent with those reported in previous studies, but with some interesting exceptions. Notably, the backbone relationships of the weevil phylogeny were consistently strongly supported, and the former Nemonychidae (pine flower snout beetles) were polyphyletic, with the subfamily Cimberidinae (here elevated to Cimberididae) placed as sister group of all other weevils. The clade comprising the sister families Brentidae (straight-snouted weevils) and Curculionidae was maximally supported and the composition of both families was firmly established. The contributions of substitution modeling, codon usage and/or mutational bias to differences between trees reconstructed from amino acid and nucleotide sequences were explored. A reconstructed timetree for weevils is consistent with a Mesozoic radiation of gymnosperm-associated taxa to form most extant families and diversification of Curculionidae alongside flowering plants-first monocots, then other groups-beginning in the Cretaceous.

Resolving Rapid Radiations within Angiosperm Families Using Anchored Phylogenomics.

Despite the promise that molecular data would provide a seemingly unlimited source of independent characters, many plant phylogenetic studies are still based on only two regions, the plastid genome and nuclear ribosomal DNA (nrDNA). Their popularity can be explained by high-copy numbers and universal polymerase chain reaction (PCR) primers that make their sequences easily amplified and converted into parallel datasets. Unfortunately, their utility is limited by linked loci and limited characters resulting in low confidence in the accuracy of phylogenetic estimates, especially when rapid radiations occur. In another contribution on anchored phylogenomics in angiosperms, we presented flowering plant-specific anchored enrichment probes for hundreds of conserved nuclear genes and demonstrated their use at the level of all angiosperms. In this contribution, we focus on a common problem in phylogenetic reconstructions below the family level: Weak or unresolved backbone due to rapid radiations ($\leqslant $10 million years) followed by long divergence, using the Cariceae-Dulichieae-Scirpeae (CDS, Cyperaceae) clade as a test case. By comparing our nuclear matrix of 461 genes to a typical Sanger-sequence dataset consisting of a few plastid genes (matK, ndhF) and an nrDNA marker (ETS), we demonstrate that our nuclear data is fully compatible with the Sanger dataset and resolves short backbone internodes with high support in both concatenated and coalescence-based analyses. In addition, we show that nuclear gene tree incongruence is inversely proportional to phylogenetic information content, indicating that incongruence is mostly due to gene tree estimation error. This suggests that large numbers of conserved nuclear loci could produce more accurate trees than sampling rapidly evolving regions prone to saturation and long-branch attraction. The robust phylogenetic estimates obtained here, and high congruence with previous morphological and molecular analyses, are strong evidence for a complete tribal revision of CDS clade. The anchored hybrid enrichment probes used in this study should be similarly effective in other flowering plant groups.

Quantity, Not Quality: Rapid Adaptation in a Polygenic Trait Proceeded Exclusively through Expression Differentiation.

A trait's genomic architecture can affect the rate and mechanism of adaptation, and although many ecologically-important traits are polygenic, most studies connecting genotype, phenotype, and fitness in natural populations have focused on traits with relatively simple genetic bases. To understand the genetic basis of polygenic adaptation, we must integrate genomics, phenotypic data, ecology, and fitness effects for a genetically tractable, polygenic trait; snake venoms provide such a system for studying polygenic adaptation because of their genetic tractability and vital ecological role in feeding and defense. We used a venom transcriptome-proteome map, quantitative proteomics, genomics, and fitness assays in sympatric prey to construct a genotype-phenotype-fitness map for the venoms of an island-mainland pair of rattlesnake populations. Reciprocal fitness experiments demonstrated that each population was locally adapted to sympatric prey. We identified significant expression differentiation with little to no coding-sequence variation across populations, demonstrating that expression differentiation was exclusively the genetic basis of polygenic adaptation. Previous research on the genetics of adaptation, however, has largely been biased toward investigating protein-coding regions because of the complexity of gene regulation. Our results showed that biases at the molecular level can be in the opposite direction, highlighting the need for more systematic comparisons of different molecular mechanisms underlying rapid, adaptive evolution in polygenic traits.