Search

Fungus-farming ants cultivate multiple lineages of fungi for food, but, because fungal cultivar relationships are largely unresolved, the history of fungus-ant coevolution remains poorly known. We designed probes targeting >2000 gene regions to generate a dated evolutionary tree for 475 fungi and combined it with a similarly generated tree for 276 ants. We found that fungus-ant agriculture originated ~66 million years ago when the end-of-Cretaceous asteroid impact temporarily interrupted photosynthesis, causing global mass extinctions but favoring the proliferation of fungi. Subsequently, ~27 million years ago, one ancestral fungal cultivar population became domesticated, i.e., obligately mutualistic, when seasonally dry habitats expanded in South America, likely isolating the cultivar population from its free-living, wet forest-dwelling conspecifics. By revealing these and other major transitions in fungus-ant coevolution, our results clarify the historical processes that shaped a model system for nonhuman agriculture.

Phylogenomic approaches have recently helped elucidate various insect relationships, but large-scale comprehensive analyses on relationships within sawflies and woodwasps are still lacking. Here, we infer the relationships and long-term biogeographic history of these hymenopteran groups using a large dataset of 354 UCE loci collected from 385 species that represent all major lineages. Early Hymenoptera started diversifying during the Early Triassic ∼249 Ma and spread all over the ancient supercontinent Pangaea. We recovered Xyeloidea as a monophyletic sister group to other Hymenoptera and Pamphilioidea as sister to Unicalcarida. Within the diverse family Tenthredinidae, our taxonomically and geographically expanded taxon sampling highlights the non-monophyly of several traditionally defined subfamilies. In addition, the recent removal of Athalia and related genera from the Tenthredinidae into the separate family Athaliidae is supported. The deep historical biogeography of the group is characterised by independent dispersals and re-colonisations between the northern (Laurasia) and southern (Gondwana) palaeocontinents. The breakup of these landmasses led to ancient vicariance in several Gondwanan lineages, while interchange across the Northern Hemisphere has continued until the Recent. The little-studied African sawfly fauna is likewise a diverse mixture of groups with varying routes of colonization. Our results reveal interesting parallels in the evolution and biogeography of early hymenopterans and other ancient insect groups., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 The Author(s). Published by Elsevier Inc. All rights reserved.)

Despite tremendous efforts in the past decades, relationships among main avian lineages remain heavily debated without a clear resolution. Discrepancies have been attributed to diversity of species sampled, phylogenetic method and the choice of genomic regions 1-3 . Here we address these issues by analysing the genomes of 363 bird species 4 (218 taxonomic families, 92% of total). Using intergenic regions and coalescent methods, we present a well-supported tree but also a marked degree of discordance. The tree confirms that Neoaves experienced rapid radiation at or near the Cretaceous-Palaeogene boundary. Sufficient loci rather than extensive taxon sampling were more effective in resolving difficult nodes. Remaining recalcitrant nodes involve species that are a challenge to model due to either extreme DNA composition, variable substitution rates, incomplete lineage sorting or complex evolutionary events such as ancient hybridization. Assessment of the effects of different genomic partitions showed high heterogeneity across the genome. We discovered sharp increases in effective population size, substitution rates and relative brain size following the Cretaceous-Palaeogene extinction event, supporting the hypothesis that emerging ecological opportunities catalysed the diversification of modern birds. The resulting phylogenetic estimate offers fresh insights into the rapid radiation of modern birds and provides a taxon-rich backbone tree for future comparative studies., (© 2024. This is a U.S. Government work and not under copyright protection in the US; foreign copyright protection may apply.)

Humans have profoundly impacted the distribution of plant and animal species over thousands of years. The most direct example of these effects is human-mediated movement of individuals, either through translocation of individuals within their range or through the introduction of species to new habitats. While human involvement may be suspected in species with obvious range disjunctions, it can be difficult to detect natural versus human-mediated dispersal events for populations at the edge of a species' range, and this uncertainty muddles how we understand the evolutionary history of populations and broad biogeographical patterns. Studies combining genetic data with archaeological, linguistic and historical evidence have confirmed prehistoric examples of human-mediated dispersal; however, it is unclear whether these methods can disentangle recent dispersal events, such as species translocated by European colonizers during the past 500 years. We use genomic DNA from historical museum specimens and historical records to evaluate three hypotheses regarding the timing and origin of Northern Bobwhites (Colinus virginianus) in Cuba, whose status as an endemic or introduced population has long been debated. We discovered that bobwhites from southern Mexico arrived in Cuba between the 12th and 16th centuries, followed by the subsequent introduction of bobwhites from the southeastern USA to Cuba between the 18th and 20th centuries. These dates suggest the introduction of bobwhites to Cuba was human-mediated and concomitant with Spanish colonial shipping routes between Veracruz, Mexico and Havana, Cuba during this period. Our results identify endemic Cuban bobwhites as a genetically distinct population born of hybridization between divergent, introduced lineages., (© 2023 John Wiley & Sons Ltd.)

Rivers frequently delimit the geographic ranges of species in the Amazon Basin. These rivers also define the boundaries between genetic clusters within many species, yet river boundaries have been documented to break down in headwater regions where rivers are narrower. To explore the evolutionary implications of headwater contact zones in Amazonia, we examined genetic variation in the Blue-capped Manakin (Lepidothrix coronata), a species previously shown to contain several genetically and phenotypically distinct populations across the western Amazon Basin. We collected restriction site-associated DNA sequence data (RADcap) for 706 individuals and found that spatial patterns of genetic structure indicate several rivers, particularly the Amazon and Ucayali, are dispersal barriers for L. coronata. We also found evidence that genetic connectivity is elevated across several headwater regions, highlighting the importance of headwater gene flow for models of Amazonian diversification. The headwater region of the Ucayali River provided a notable exception to findings of headwater gene flow by harboring non-admixed populations of L. coronata on opposite sides of a < 1-km-wide river channel with a known dynamic history, suggesting that additional prezygotic barriers may be limiting gene flow in this region., (© The Author(s) 2023. Published by Oxford University Press on behalf of The Society for the Study of Evolution (SSE). All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

A universal paradigm describing patterns of speciation across the tree of life has been debated for decades. In marine organisms, inferring patterns of speciation using contemporary and historical patterns of biogeography is challenging due to the deficiency of species-level phylogenies and information on species' distributions, as well as conflicting relationships between species' dispersal, range size and co-occurrence. Most research on global patterns of marine fish speciation and biogeography has focused on coral reef or pelagic species. Carangoidei is an ecologically important clade of marine fishes that use coral reef and pelagic environments. We used sequence capture of 1314 ultraconserved elements (UCEs) from 154 taxa to generate a time-calibrated phylogeny of Carangoidei and its parent clade, Carangiformes. Age-range correlation analyses of the geographical distributions and divergence times of sister species pairs reveal widespread sympatry, with 73% of sister species pairs exhibiting sympatric geographical distributions, regardless of node age. Most species pairs coexist across large portions of their ranges. We also observe greater disparity in body length and maximum depth between sympatric relative to allopatric sister species. These and other ecological or behavioural attributes probably facilitate sympatry among the most closely related carangoids.

The clapper rail (Rallus crepitans), of the family Rallidae, is a secretive marsh bird species that is adapted for high salinity habitats. They are very similar in appearance to the closely related king rail (R. elegans), but while king rails are limited primarily to freshwater marshes, clapper rails are highly adapted to tolerate salt marshes. Both species can be found in brackish marshes where they freely hybridize, but the distribution of their respective habitats precludes the formation of a continuous hybrid zone and secondary contact can occur repeatedly. This system, thus, provides unique opportunities to investigate the underlying mechanisms driving their differential salinity tolerance as well as the maintenance of the species boundary between the 2 species. To facilitate these studies, we assembled a de novo reference genome assembly for a female clapper rail. Chicago and HiC libraries were prepared as input for the Dovetail HiRise pipeline to scaffold the genome. The pipeline, however, did not recover the Z chromosome so a custom script was used to assemble the Z chromosome. We generated a near chromosome level assembly with a total length of 994.8 Mb comprising 13,226 scaffolds. The assembly had a scaffold N50 was 82.7 Mb, L50 of four, and had a BUSCO completeness score of 92%. This assembly is among the most contiguous genomes among the species in the family Rallidae. It will serve as an important tool in future studies on avian salinity tolerance, interspecific hybridization, and speciation., Competing Interests: Conflicts of interest The authors declare no conflict of interest., (© The Author(s) 2023. Published by Oxford University Press on behalf of The Genetics Society of America.)

Bacterioplankton of the SAR11 clade are the most abundant marine microorganisms and consist of numerous subclades spanning order-level divergence ( Pelagibacterales ). The assignment of the earliest diverging subclade V (a.k.a. HIMB59) to the Pelagibacterales is highly controversial, with multiple recent phylogenetic studies placing them completely separate from SAR11. Other than through phylogenomics, subclade V has not received detailed examination due to limited genomes from this group. Here, we assessed the ecogenomic characteristics of subclade V to better understand the role of this group in comparison to the Pelagibacterales . We used a new isolate genome, recently released single-amplified genomes and metagenome-assembled genomes, and previously established SAR11 genomes to perform a comprehensive comparative genomics analysis. We paired this analysis with the recruitment of metagenomes spanning the open ocean, coastal, and brackish systems. Phylogenomics, average amino acid identity, and 16S rRNA gene phylogeny indicate that SAR11 subclade V is synonymous with the ubiquitous AEGEAN-169 clade and support the contention that this group represents a taxonomic family. AEGEAN-169 shared many bulk genome qualities with SAR11, such as streamlining and low GC content, but genomes were generally larger. AEGEAN-169 had overlapping distributions with SAR11 but was metabolically distinct from SAR11 in its potential to transport and utilize a broader range of sugars as well as in the transport of trace metals and thiamin. Thus, regardless of the ultimate phylogenetic placement of AEGEAN-169, these organisms have distinct metabolic capacities that likely allow them to differentiate their niche from canonical SAR11 taxa. IMPORTANCE One goal of marine microbiologists is to uncover the roles various microorganisms are playing in biogeochemical cycles. Success in this endeavor relies on differentiating groups of microbes and circumscribing their relationships. An early-diverging group (subclade V) of the most abundant bacterioplankton, SAR11, has recently been proposed as a separate lineage that does not share a most recent common ancestor. But beyond phylogenetics, little has been done to evaluate how these organisms compare with SAR11. Our work leverages dozens of new genomes to demonstrate the similarities and differences between subclade V and SAR11. In our analysis, we also establish that subclade V is synonymous with a group of bacteria established from 16S rRNA gene sequences, AEGEAN-169. Subclade V/AEGEAN-169 has clear metabolic distinctions from SAR11 and their shared traits point to remarkable convergent evolution if they do not share a most recent common ancestor., Competing Interests: The authors declare no conflict of interest.

The Order Pelagibacterales (SAR11) is the most abundant group of heterotrophic bacterioplankton in global oceans and comprises multiple subclades with unique spatiotemporal distributions. Subclade IIIa is the primary SAR11 group in brackish waters and shares a common ancestor with the dominant freshwater IIIb (LD12) subclade. Despite its dominance in brackish environments, subclade IIIa lacks systematic genomic or ecological studies. Here, we combine closed genomes from new IIIa isolates, new IIIa MAGS from San Francisco Bay (SFB), and 460 highly complete publicly available SAR11 genomes for the most comprehensive pangenomic study of subclade IIIa to date. Subclade IIIa represents a taxonomic family containing three genera (denoted as subgroups IIIa.1, IIIa.2, and IIIa.3) that had distinct ecological distributions related to salinity. The expansion of taxon selection within subclade IIIa also established previously noted metabolic differentiation in subclade IIIa compared to other SAR11 subclades such as glycine/serine prototrophy, mosaic glyoxylate shunt presence, and polyhydroxyalkanoate synthesis potential. Our analysis further shows metabolic flexibility among subgroups within IIIa. Additionally, we find that subclade IIIa.3 bridges the marine and freshwater clades based on its potential for compatible solute transport, iron utilization, and bicarbonate management potential. Pure culture experimentation validated differential salinity ranges in IIIa.1 and IIIa.3 and provided detailed IIIa cell size and volume data. This study is an important step forward for understanding the genomic, ecological, and physiological differentiation of subclade IIIa and the overall evolutionary history of SAR11., (© 2023. The Author(s).)

North American sunfishes (Family Centrarchidae) are among the most popular sportfish throughout the United States and Canada. Despite the popularity of sunfishes, their ecological importance, and their extensive stocking and aquacultural history, few molecular studies have examined the evolutionary relationships and species boundaries among members of this group, many of which are known to hybridize. Here, we describe a chromosome-scale genome assembly representing Bluegill (Lepomis macrochirus), one of the most widespread centrarchid species. By combining long-read, Oxford Nanopore sequencing data with short-insert, whole-genome and HiC sequence reads, we produced an assembly (Lm&#95;LA&#95;1.1) having a total length of 889 Mb including 1,841 scaffolds and having a scaffold N50 of 36 Mb, L50 of 12, N90 of 29 Mb, and L90 of 22. We detected 99% (eukaryota&#95;odb10) and 98% (actinopterygii&#95;odb10) universal single-copy orthologs (BUSCOs), and ab initio gene prediction performed using this new assembly identified a set of 17,233 genes that were supported by external (OrthoDB v10) data. This new assembly provides an important addition to the growing set of assemblies already available for spiny-rayed fishes (Acanthomorpha), and it will serve as a resource for future studies that focus on the complex evolutionary history of centrarchids., Competing Interests: Conflicts of interest None declared., (© The Author(s) 2023. Published by Oxford University Press on behalf of the Genetics Society of America.)

The Gulf Coast tick, Amblyomma maculatum, inhabits the Southeastern states of the USA bordering the Gulf of Mexico, Mexico, and other Central and South American countries. More recently, its U.S. range has extended West to Arizona and Northeast to New York state and Connecticut. It is a vector of Rickettsia parkeri and Hepatozoon americanum. This tick species has become a model to study tick/Rickettsia interactions. To increase our knowledge of the basic biology of A. maculatum we report here a draft genome of this tick and an extensive functional classification of its proteome. The DNA from a single male tick was used as a genomic source, and a 10X genomics protocol determined 28,460 scaffolds having equal or more than 10 Kb, totaling 1.98 Gb. The N50 scaffold size was 19,849 Kb. The BRAKER pipeline was used to find the protein-coding gene boundaries on the assembled A. maculatum genome, discovering 237,921 CDS. After trimming and classifying the transposable elements, bacterial contaminants, and truncated genes, a set of 25,702 were annotated and classified as the core gene products. A BUSCO analysis revealed 83.4% complete BUSCOs. A hyperlinked spreadsheet is provided, allowing browsing of the individual gene products and their matches to several databases., (Copyright © 2022. Published by Elsevier GmbH.)

As phylogenomics focuses on comprehensive taxon sampling at the species and population/subspecies levels, incorporating genomic data from historical specimens has become increasingly common. While historical samples can fill critical gaps in our understanding of the evolutionary history of diverse groups, they also introduce additional sources of phylogenomic uncertainty, making it difficult to discern novel evolutionary relationships from artifacts caused by sample quality issues. These problems highlight the need for improved strategies to disentangle artifactual patterns from true biological signal as historical specimens become more prevalent in phylogenomic datasets. Here, we tested the limits of historical specimen-driven phylogenomics to resolve subspecies-level relationships within a highly polytypic family, the New World quails (Odontophoridae), using thousands of ultraconserved elements (UCEs). We found that relationships at and above the species-level were well-resolved and highly supported across all analyses, with the exception of discordant relationships within the two most polytypic genera which included many historical specimens. We examined the causes of discordance and found that inferring phylogenies from subsets of taxa resolved the disagreements, suggesting that analyzing subclades can help remove artifactual causes of discordance in datasets that include historical samples. At the subspecies-level, we found well-resolved geographic structure within the two most polytypic genera, including the most polytypic species in this family, Northern Bobwhites (Colinus virginianus), demonstrating that variable sites within UCEs are capable of resolving phylogenetic structure below the species level. Our results highlight the importance of complete taxonomic sampling for resolving relationships among polytypic species, often through the inclusion of historical specimens, and we propose an integrative strategy for understanding and addressing the uncertainty that historical samples sometimes introduce to phylogenetic analyses., (Copyright © 2022 Elsevier Inc. All rights reserved.)

Spiny-rayed fishes (Acanthomorpha) dominate modern marine habitats and account for more than a quarter of all living vertebrate species. Previous time-calibrated phylogenies and patterns from the fossil record explain this dominance by correlating the origin of major acanthomorph lineages with the Cretaceous-Palaeogene mass extinction. Here we infer a time-calibrated phylogeny using ultraconserved elements that samples 91.4% of all acanthomorph families and investigate patterns of body shape disparity. Our results show that acanthomorph lineages steadily accumulated throughout the Cenozoic and underwent a significant expansion of among-clade morphological disparity several million years after the end-Cretaceous. These acanthomorph lineages radiated into and diversified within distinct regions of morphospace that characterize iconic lineages, including fast-swimming open-ocean predators, laterally compressed reef fishes, bottom-dwelling flatfishes, seahorses and pufferfishes. The evolutionary success of spiny-rayed fishes is the culmination of multiple species-rich and phenotypically disparate lineages independently diversifying across the globe under a wide range of ecological conditions., (© 2022. The Author(s), under exclusive licence to Springer Nature Limited.)

Although recent molecular phylogenetic analyses of Lepidothrix manakins (family Pipridae) have helped clarify their evolutionary relationships, the placement of several lineages remains in question because of low or conflicting branch support. In particular, the relationship of L. coronata to other members of the genus and relationships within the L. nattereri + L. vilasboasi + L. iris clade have been difficult to resolve. We used RADcap to collect restriction site-associated DNA sequence data and estimate the first subspecies-level phylogeny of the genus Lepidothrix (17 of 18 currently recognized subspecies), and we included extensive geographic representation of the widespread and phenotypically variable L. coronata. We found strong support for the phylogenetic position and monophyly of L. coronata, and we resolved two clades separated by the Andes that, along with previous divergence time estimates and our assessment of morphological and vocal evidence, suggest the presence of two biological species: Velvety Manakin (L. velutina) west of the Andes and Blue-capped Manakin (L. coronata) east of the Andes. Species-level relationships within the L. nattereri + L. vilasboasi + L. iris clade remained poorly resolved in concatenated and coalescent-based analyses, with SNAPP analyses suggesting that the lack of reciprocal monophyly is due to extensive allele sharing among these taxa. Finally, we confirmed a previously documented hybrid between L. coronata and L. suavissima as an F1 individual, consistent with the view that hybridization between these two species is a rare event and that postmating reproductive barriers prevent successful backcrossing., (Copyright © 2022 Elsevier Inc. All rights reserved.)

Secondary contact between species often results in the formation of a hybrid zone, with the eventual fates of the hybridizing species dependent on evolutionary and ecological forces. We examine this process in the Amazon Basin by conducting the first genomic and phenotypic characterization of the hybrid zone formed after secondary contact between two obligate army-ant-followers: the White-breasted Antbird (Rhegmatorhina hoffmannsi) and the Harlequin Antbird (Rhegmatorhina berlepschi). We found a major geographic displacement (∼120 km) between the mitochondrial and nuclear clines, and we explore potential hypotheses for the displacement, including sampling error, genetic drift, and asymmetric cytonuclear incompatibilities. We cannot exclude roles for sampling error and genetic drift in contributing to the discordance; however, the data suggest expansion and unidirectional introgression of hoffmannsi into the distribution of berlepschi., (© 2021 The Authors. Evolution © 2021 The Society for the Study of Evolution.)

The Neotropics harbor the most species-rich freshwater fish fauna on the planet, but the timing of that exceptional diversification remains unclear. Did the Neotropics accumulate species steadily throughout their long history, or attain their remarkable diversity recently? Biologists have long debated the relative support for these museum and cradle hypotheses, but few phylogenies of megadiverse tropical clades have included sufficient taxa to distinguish between them. We used 1288 ultraconserved element loci spanning 293 species, 211 genera, and 21 families of characoid fishes to reconstruct a new, fossil-calibrated phylogeny and infer the most likely diversification scenario for a clade that includes a third of Neotropical fish diversity. This phylogeny implies paraphyly of the traditional delimitation of Characiformes because it resolves the largely Neotropical Characoidei as the sister lineage of Siluriformes (catfishes), rather than the African Citharinodei. Time-calibrated phylogenies indicate an ancient origin of major characoid lineages and reveal a much more recent emergence of most characoid species. Diversification rate analyses infer increased speciation and decreased extinction rates during the Oligocene at around 30 Ma during a period of mega-wetland formation in the proto-Orinoco-Amazonas. Three species-rich and ecomorphologically diverse lineages (Anostomidae, Serrasalmidae, and Characidae) that originated more than 60 Ma in the Paleocene experienced particularly notable bursts of Oligocene diversification and now account collectively for 68% of the approximately 2150 species of Characoidei. In addition to paleogeographic changes, we discuss potential accelerants of diversification in these three lineages. While the Neotropics accumulated a museum of ecomorphologically diverse characoid lineages long ago, this geologically dynamic region also cradled a much more recent birth of remarkable species-level diversity. [Biodiversity; Characiformes; macroevolution; Neotropics; phylogenomics; ultraconserved elements.]., (© The Author(s) 2021. Published by Oxford University Press, on behalf of the Society of Systematic Biologists. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

Black-throated Flowerpiercers (Diglossa brunneiventris) are one species representing a phenotypically specialized group of tanagers (Thraupidae) that have hooked bills which allow them to feed by stealing nectar from the base of flowers. Members of the genus are widely distributed in montane regions from Mexico to northern Argentina, and previous studies of Diglossa have focused on their systematics, phylogenetics, and interesting natural history. Despite numerous studies of species within the genus, no genome assembly exists to represent these nectivorous tanagers. We described the assembly of a genome sequence representing a museum-vouchered, wild, female D. brunneiventris collected in Peru. By combining Pacific Biosciences Sequel long-read technology with 10× linked-read and reference-based scaffolding, we produced a 1.08 Gbp pseudochromosomal assembly including 600 scaffolds with a scaffold N50 of 67.3 Mbp, a scaffold L50 of 6, and a BUSCO completeness score of 95%. This new assembly improves representation of the diverse species that comprise the tanagers, improves on scaffold lengths and contiguity when compared to existing genomic resources for tanagers, and provides another avenue of research into the genetic basis of adaptations common to a nectivorous lifestyle among vertebrates., (© The Author(s) 2021. Published by Oxford University Press on behalf of Genetics Society of America.)

Several bird taxa have been recently described or elevated to full species and almost twice as many bird species than are currently recognized may exist. Defining species is one of the most basic and important issues in biological science because unknown or poorly defined species hamper subsequent studies. Here, we evaluate the species limits and evolutionary history of Tunchiornis ochraceiceps-a widespread forest songbird that occurs in the lowlands of Central America, Chocó and Amazonia-using an integrative approach that includes plumage coloration, morphometrics, vocalization and genomic data. The species has a relatively old crown age (~9 Ma) and comprises several lineages with little, if any, evidence of gene flow among them. We propose a taxonomic arrangement composed of four species, three with a plumage coloration diagnosis and one deeply divergent cryptic species. Most of the remaining lineages have variable but unfixed phenotypic characters despite their relatively old origin. This decoupling of genomic and phenotypic differentiation reveals a remarkable case of phenotypic conservatism, possibly due to strict habitat association. Lineages are geographically delimited by the main Amazonian rivers and the Andes, a pattern observed in studies of other understory upland forest Neotropical birds, although phylogenetic relationships and divergence times among populations are idiosyncratic., (Copyright © 2021 Elsevier Inc. All rights reserved.)

A voucher is a permanently preserved specimen that is maintained in an accessible collection. In genomics, vouchers serve as the physical evidence for the taxonomic identification of genome assemblies. Unfortunately, the vast majority of vertebrate genomes stored in the GenBank database do not refer to voucher specimens. Here, we urge researchers generating new genome assemblies to deposit voucher specimens in accessible, permanent research collections, and to link these vouchers to publications, public databases, and repositories. We also encourage scientists to deposit voucher specimens in order to recognize the work of local field biologists and promote a diverse and inclusive knowledge base, and we recommend best practices for voucher deposition to prevent taxonomic errors and ensure reproducibility and legality in genetic studies., Competing Interests: JB, RS, BF, PC No competing interests declared, (© 2021, Buckner et al.)

During the Last Glacial Maximum (LGM), global sea levels were 120-130 m lower than today, resulting in the emergence of most continental shelves and extirpation of subtidal organisms from these areas. During the interglacial periods, rapid inundation of shelf regions created a dynamic environment for coastal organisms, such as the charismatic leafy seadragon (Phycodurus eques, Syngnathidae), a brooder with low dispersal ability inhabiting kelp beds in temperate Australia. Reconstructions of the palaeoshoreline revealed that the increase of shallow areas since the LGM was not uniform across the species' range and we investigated the effects of these asymmetries on genetic diversity and structuring. Using targeted capture of 857 variable ultraconserved elements (UCEs, 2,845 single nucleotide polymorphisms) in 68 individuals, we found that the regionally different shelf topographies were paralleled by contrasting population genetic patterns. In the west, populations may not have persisted through sea-level lows because shallow seabed was very limited. Shallow genetic structure, weak expansion signals and a westward cline in genetic diversity indicate a postglacial recolonization of the western part of the range from a more eastern location following sea-level rise. In the east, shallow seabed persisted during the LGM and increased considerably after the flooding of large bays, which resulted in strong demographic expansions, deeper genetic structure and higher genetic diversity. This study suggests that postglacial flooding with rising sea levels produced locally variable signatures in colonizing populations., (© 2020 John Wiley & Sons Ltd.)

The lack of genomic resources for tropical canopy trees is impeding several research avenues in tropical forest biology. We present genome assemblies for two Neotropical hardwood species, Jacaranda copaia and Handroanthus (formerly Tabebuia) guayacan, that are model systems for research on tropical tree demography and flowering phenology. For each species, we combined Illumina short-read data with in vitro proximity-ligation (Chicago) libraries to generate an assembly. For Jacaranda copaia, we obtained 104X physical coverage and produced an assembly with N50/N90 scaffold lengths of 1.020/0.277 Mbp. For H. guayacan, we obtained 129X coverage and produced an assembly with N50/N90 scaffold lengths of 0.795/0.165 Mbp. J. copaia and H. guayacan assemblies contained 95.8% and 87.9% of benchmarking orthologs, although they constituted only 77.1% and 66.7% of the estimated genome sizes of 799 and 512 Mbp, respectively. These differences were potentially due to high repetitive sequence content (>59.31% and 45.59%) and high heterozygosity (0.5% and 0.8%) in each species. Finally, we compared each new assembly to a previously sequenced genome for Handroanthus impetiginosus using whole-genome alignment. This analysis indicated extensive gene duplication in H. impetiginosus since its divergence from H. guayacan., (© The Author(s) 2021. Published by Oxford University Press on behalf of Genetics Society of America.)

The tropics are the source of most biodiversity yet inadequate sampling obscures answers to fundamental questions about how this diversity evolves. We leveraged samples assembled over decades of fieldwork to study diversification of the largest tropical bird radiation, the suboscine passerines. Our phylogeny, estimated using data from 2389 genomic regions in 1940 individuals of 1283 species, reveals that peak suboscine species diversity in the Neotropics is not associated with high recent speciation rates but rather with the gradual accumulation of species over time. Paradoxically, the highest speciation rates are in lineages from regions with low species diversity, which are generally cold, dry, unstable environments. Our results reveal a model in which species are forming faster in environmental extremes but have accumulated in moderate environments to form tropical biodiversity hotspots., (Copyright © 2020 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

Whole-genome sequencing projects are increasingly populating the tree of life and characterizing biodiversity 1-4 . Sparse taxon sampling has previously been proposed to confound phylogenetic inference 5 , and captures only a fraction of the genomic diversity. Here we report a substantial step towards the dense representation of avian phylogenetic and molecular diversity, by analysing 363 genomes from 92.4% of bird families-including 267 newly sequenced genomes produced for phase II of the Bird 10,000 Genomes (B10K) Project. We use this comparative genome dataset in combination with a pipeline that leverages a reference-free whole-genome alignment to identify orthologous regions in greater numbers than has previously been possible and to recognize genomic novelties in particular bird lineages. The densely sampled alignment provides a single-base-pair map of selection, has more than doubled the fraction of bases that are confidently predicted to be under conservation and reveals extensive patterns of weak selection in predominantly non-coding DNA. Our results demonstrate that increasing the diversity of genomes used in comparative studies can reveal more shared and lineage-specific variation, and improve the investigation of genomic characteristics. We anticipate that this genomic resource will offer new perspectives on evolutionary processes in cross-species comparative analyses and assist in efforts to conserve species.

Identifying how past environmental conditions shaped the evolution of corals and their skeletal traits provides a framework for predicting their persistence and that of their non-calcifying relatives under impending global warming and ocean acidification. Here we show that ocean geochemistry, particularly aragonite-calcite seas, drives patterns of morphological evolution in anthozoans (corals, sea anemones) by examining skeletal traits in the context of a robust, time-calibrated phylogeny. The lability of skeletal composition among octocorals suggests a greater ability to adapt to changes in ocean chemistry compared with the homogeneity of the aragonitic skeleton of scleractinian corals. Pulses of diversification in anthozoans follow mass extinctions and reef crises, with sea anemones and proteinaceous corals filling empty niches as tropical reef builders went extinct. Changing environmental conditions will likely diminish aragonitic reef-building scleractinians, but the evolutionary history of the Anthozoa suggests other groups will persist and diversify in their wake.

Determining how genetic diversity is structured between populations that span the divergence continuum from populations to biological species is key to understanding the generation and maintenance of biodiversity. We investigated genetic divergence and gene flow in eight lineages of birds with a trans-Beringian distribution, where Asian and North American populations have likely been split and reunited through multiple Pleistocene glacial cycles. Our study transects the speciation process, including eight pairwise comparisons in three orders (ducks, shorebirds and passerines) at population, subspecies and species levels. Using ultraconserved elements (UCEs), we found that these lineages represent conditions from slightly differentiated populations to full biological species. Although allopatric speciation is considered the predominant mode of divergence in birds, all of our best divergence models included gene flow, supporting speciation with gene flow as the predominant mode in Beringia. In our eight lineages, three were best described by a split-migration model (divergence with gene flow), three best fit a secondary contact scenario (isolation followed by gene flow), and two showed support for both models. The lineages were not evenly distributed across a divergence space defined by gene flow (M) and differentiation (F ST ), instead forming two discontinuous groups: one with relatively shallow divergence, no fixed single nucleotide polymorphisms (SNPs), and high rates of gene flow between populations; and the second with relatively deeply divergent lineages, multiple fixed SNPs, and low gene flow. Our results highlight the important role that gene flow plays in avian divergence in Beringia., (© 2020 John Wiley & Sons Ltd.)

Northern bobwhites ( Colinus virginianus ) are small quails in the New World Quail family (Odontophoridae) and are one of the most phenotypically diverse avian species. Despite extensive research on bobwhite ecology, genomic studies investigating the evolution of phenotypic diversity in this species are lacking. Here, we present a new, highly contiguous assembly for bobwhites using tissue samples from a vouchered, wild, female bird collected in Louisiana. By performing a de novo assembly and scaffolding the assembly with Dovetail Chicago and HiC libraries and the HiRise pipeline, we produced an 866.8 Mb assembly including 1,512 scaffolds with a scaffold N50 of 66.8 Mb, a scaffold L90 of 17, and a BUSCO completeness score of 90.8%. This new assembly represents approximately 96% of the non-repetitive and 84% of the entire bobwhite genome size, greatly improves scaffold lengths and contiguity compared to an existing draft bobwhite genome, and provides an important tool for future studies of evolutionary and functional genomics in bobwhites., (Copyright © 2019 Salter et al.)

The advent of next-generation sequencing allows researchers to use large-scale datasets for species delimitation analyses, yet one can envision an inflection point where the added accuracy of including more loci does not offset the increased computational burden. One alternative to including all loci could be to prioritize the analysis of loci for which there is an expectation of high informativeness. Here, we explore the issue of species delimitation and locus selection with montane species from two anuran genera that have been isolated in sky islands across the southern Brazilian Atlantic Forest: Melanophryniscus (Bufonidae) and Brachycephalus (Brachycephalidae). To delimit species, we obtained genetic data using target enrichment of ultraconserved elements from 32 populations (13 for Melanophryniscus and 19 for Brachycephalus), and we were able to create datasets that included over 800 loci with no missing data. We ranked loci according to their number of parsimony-informative sites, and we performed species delimitation analyses using BPP with the most informative 10, 20, 40, 80, 160, 320, and 640 loci. We identified three types of phylogenetic node: nodes with either consistently high or low support regardless of the number of loci or their informativeness and nodes that were initially poorly supported where support became stronger as we included more data. When viewed across all sensitivity analyses, our results suggest that the current species richness in both genera is likely underestimated. In addition, our results show the effects of different sampling strategies on species delimitation using phylogenomic datasets., (Copyright © 2019 Elsevier Inc. All rights reserved.)

Massively parallel DNA sequencing offers many benefits, but major inhibitory cost factors include: (1) start-up (i.e., purchasing initial reagents and equipment); (2) buy-in (i.e., getting the smallest possible amount of data from a run); and (3) sample preparation. Reducing sample preparation costs is commonly addressed, but start-up and buy-in costs are rarely addressed. We present dual-indexing systems to address all three of these issues. By breaking the library construction process into universal, re-usable, combinatorial components, we reduce all costs, while increasing the number of samples and the variety of library types that can be combined within runs. We accomplish this by extending the Illumina TruSeq dual-indexing approach to 768 (384 + 384) indexed primers that produce 384 unique dual-indexes or 147,456 (384 × 384) unique combinations. We maintain eight nucleotide indexes, with many that are compatible with Illumina index sequences. We synthesized these indexing primers, purifying them with only standard desalting and placing small aliquots in replicate plates. In qPCR validation tests, 206 of 208 primers tested passed (99% success). We then created hundreds of libraries in various scenarios. Our approach reduces start-up and per-sample costs by requiring only one universal adapter that works with indexed PCR primers to uniquely identify samples. Our approach reduces buy-in costs because: (1) relatively few oligonucleotides are needed to produce a large number of indexed libraries; and (2) the large number of possible primers allows researchers to use unique primer sets for different projects, which facilitates pooling of samples during sequencing. Our libraries make use of standard Illumina sequencing primers and index sequence length and are demultiplexed with standard Illumina software, thereby minimizing customization headaches. In subsequent Adapterama papers, we use these same primers with different adapter stubs to construct amplicon and restriction-site associated DNA libraries, but their use can be expanded to any type of library sequenced on Illumina platforms., Competing Interests: Brant C. Faircloth is an Academic Editor for PeerJ. The EHS DNA lab (baddna.uga.edu) and Georgia Genomics Bioinformatics Core (dna.uga.edu) provide oligonucleotide aliquots and library preparation services at cost, including some oligonucleotides and services that make use of the adapters and primers presented in this manuscript., (©2019 Glenn et al.)

Next-generation sequencing (NGS) of amplicons is used in a wide variety of contexts. In many cases, NGS amplicon sequencing remains overly expensive and inflexible, with library preparation strategies relying upon the fusion of locus-specific primers to full-length adapter sequences with a single identifying sequence or ligating adapters onto PCR products. In Adapterama I , we presented universal stubs and primers to produce thousands of unique index combinations and a modifiable system for incorporating them into Illumina libraries. Here, we describe multiple ways to use the Adapterama system and other approaches for amplicon sequencing on Illumina instruments. In the variant we use most frequently for large-scale projects, we fuse partial adapter sequences (TruSeq or Nextera) onto the 5' end of locus-specific PCR primers with variable-length tag sequences between the adapter and locus-specific sequences. These fusion primers can be used combinatorially to amplify samples within a 96-well plate (8 forward primers + 12 reverse primers yield 8 × 12 = 96 combinations), and the resulting amplicons can be pooled. The initial PCR products then serve as template for a second round of PCR with dual-indexed iTru or iNext primers (also used combinatorially) to make full-length libraries. The resulting quadruple-indexed amplicons have diversity at most base positions and can be pooled with any standard Illumina library for sequencing. The number of sequencing reads from the amplicon pools can be adjusted, facilitating deep sequencing when required or reducing sequencing costs per sample to an economically trivial amount when deep coverage is not needed. We demonstrate the utility and versatility of our approaches with results from six projects using different implementations of our protocols. Thus, we show that these methods facilitate amplicon library construction for Illumina instruments at reduced cost with increased flexibility. A simple web page to design fusion primers compatible with iTru primers is available at: http://baddna.uga.edu/tools-taggi.html. A fast and easy to use program to demultiplex amplicon pools with internal indexes is available at: https://github.com/lefeverde/Mr&#95;Demuxy., Competing Interests: Brant C. Faircloth is an Academic Editor for PeerJ. The EHS DNA lab provides oligonucleotide aliquots and library preparation services at cost, including some oligonucleotides and services that make use of the adapters and primers presented in this manuscript (baddna.uga.edu). This document has been reviewed in accordance with U.S Environmental Protection Agency policy and approved for publication. Any mention of trade names, manufacturer or products does not imply an endorsement by the United States Government or the U.S. Environmental Protection Agency. EPA and its employees do not endorse any commercial products., (©2019 Glenn et al.)

Molecular ecologists frequently use genome reduction strategies that rely upon restriction enzyme digestion of genomic DNA to sample consistent portions of the genome from many individuals (e.g., RADseq, GBS). However, researchers often find the existing methods expensive to initiate and/or difficult to implement consistently, especially because it is difficult to multiplex sufficient numbers of samples to fill entire sequencing lanes. Here, we introduce a low-cost and highly robust approach for the construction of dual-digest RADseq libraries that build on adapters and primers designed in Adapterama I . Major features of our method include: (1) minimizing the number of processing steps; (2) focusing on a single strand of sample DNA for library construction, allowing the use of a non-phosphorylated adapter on one end; (3) ligating adapters in the presence of active restriction enzymes, thereby reducing chimeras; (4) including an optional third restriction enzyme to cut apart adapter-dimers formed by the phosphorylated adapter, thus increasing the efficiency of adapter ligation to sample DNA, which is particularly effective when only low quantity/quality DNA samples are available; (5) interchangeable adapter designs; (6) incorporating variable-length internal indexes within the adapters to increase the scope of sample indexing, facilitate pooling, and increase sequence diversity; (7) maintaining compatibility with universal dual-indexed primers and thus, Illumina sequencing reagents and libraries; and, (8) easy modification for the identification of PCR duplicates. We present eight adapter designs that work with 72 restriction enzyme combinations. We demonstrate the efficiency of our approach by comparing it with existing methods, and we validate its utility through the discovery of many variable loci in a variety of non-model organisms. Our 2RAD/3RAD method is easy to perform, has low startup costs, has increased utility with low-concentration input DNA, and produces libraries that can be highly-multiplexed and pooled with other Illumina libraries., Competing Interests: Brant C. Faircloth is an Academic Editor for PeerJ. Kerin E. Bentley is currently affiliated with LeafWorks, Inc. The EHS DNA lab (http://baddna.uga.edu) provides oligonucleotide aliquots and library preparation services at cost, including some oligonucleotides and services that make use of the adapters and primers., (©2019 Bayona-Vásquez et al.)

The fish clade Pelagiaria, which includes tunas as its most famous members, evolved remarkable morphological and ecological variety in a setting not generally considered conducive to diversification: the open ocean. Relationships within Pelagiaria have proven elusive due to short internodes subtending major lineages suggestive of rapid early divergences. Using a novel sequence dataset of over 1000 ultraconserved DNA elements (UCEs) for 94 of the 286 species of Pelagiaria (more than 70% of genera), we provide a time-calibrated phylogeny for this widely distributed clade. Some inferred relationships have clear precedents (e.g. the monophyly of 'core' Stromateoidei, and a clade comprising 'Gempylidae' and Trichiuridae), but others are unexpected despite strong support (e.g. Chiasmodontidae + Tetragonurus ). Relaxed molecular clock analysis using node-based fossil calibrations estimates a latest Cretaceous origin for Pelagiaria, with crown-group families restricted to the Cenozoic. Estimated mean speciation rates decline from the origin of the group in the latest Cretaceous, although credible intervals for root and tip rates are broad and overlap in most cases, and there is higher-than-expected partitioning of body shape diversity (measured as fineness ratio) between clades concentrated during the Palaeocene-Eocene. By contrast, more direct measures of ecology show either no substantial deviation from a null model of diversification (diet) or patterns consistent with evolutionary constraint or high rates of recent change (depth habitat). Collectively, these results indicate a mosaic model of diversification. Pelagiarians show high morphological disparity and modest species richness compared to better-studied fish radiations in contrasting environments. However, this pattern is also apparent in other clades in open-ocean or deep-sea habitats, and suggests that comparative study of such groups might provide a more inclusive model of the evolution of diversity in fishes.

Objective: Hybrid zones are geographic regions where genetically distinct taxa interbreed, resulting in offspring of mixed ancestry. California Scrub-Jays (Aphelocoma californica) and Woodhouse's Scrub-Jays (A. woodhouseii) come into secondary contact and hybridize in western Nevada. Although previous work investigated divergence and gene flow between these species using a handful of microsatellite markers, the hybrid zone has not been studied using genome-scale markers, providing an opportunity to assess genome-wide introgression, test for a genetic basis for ecomorphological traits, and compare these estimates to those derived from microsatellites., Results: Using variant sites flanking ultraconserved elements (UCEs), we performed population assignment and quantified hybrid ancestry for 16 individuals across the zone of secondary contact. Our study included 2468 SNPs distributed throughout the genome, allowing discrimination of genetic affinities of hybrid individuals that were similar to estimates from microsatellites. We show a relationship between bill and wing length and the genetic composition of individuals that was not found in prior work using microsatellites, suggesting a genetic basis for these traits. Our analyses demonstrate the utility of UCEs for the analysis of hybrid zones and provide a basis for future studies to identify the genomic architecture of speciation and phenotypic differences between these incipient species.

Resolving patterns of ancient and rapid diversifications is one of the most challenging tasks in evolutionary biology. These difficulties arise from confusing phylogenetic signals that are associated with the interplay of incomplete lineage sorting (ILS) and homoplasy. Phylogenomic analyses of hundreds, or even thousands, of loci offer the potential to resolve such contentious relationships. Yet, how much useful phylogenetic information these large data sets contain remains uncertain and often goes untested. Here, we assess the utility of different data filtering approaches to maximize phylogenetic information and minimize noise when reconstructing an ancient radiation of Neotropical electric knifefishes (Order Gymnotiformes) using ultraconserved elements. We found two contrasting hypotheses of gymnotiform evolutionary relationships depending on whether phylogenetic inferences were based on concatenation or coalescent methods. In the first case, all analyses inferred a previously-and commonly-proposed hypothesis, where the family Apteronotidae was found as the sister group to all other gymnotiform families. In contrast, coalescent-based analyses suggested a novel hypothesis where families producing pulse-type (viz., Gymnotidae, Hypopomidae, and Rhamphichthyidae) and wave-type electric signals (viz., Apteronotidae, Sternopygidae) were reciprocally monophyletic. Nodal support for this second hypothesis increased when analyzing loci with the highest phylogenetic information content and further increased when data were pruned using targeted filtering methods that maximized phylogenetic informativeness at the deepest nodes of the Gymnotiformes. Bayesian concordance analyses and topology tests of individual gene genealogies demonstrated that the difficulty of resolving this radiation was likely due to high gene-tree incongruences that resulted from ILS. We show that data filtering reduces gene-tree heterogeneity and increases nodal support and consistency of species trees using coalescent methods; however, we failed to observe the same effect when using concatenation methods. Furthermore, the targeted filtering strategies applied here support the use of "gene data interrogation" rather than "gene genealogy interrogation" approaches in phylogenomic analyses, to extract phylogenetic signal from intractable portions of the Tree of Life., (© The Author(s) 2018. Published by Oxford University Press, on behalf of the Society of Systematic Biologists. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

Avian diversification has been influenced by global climate change, plate tectonic movements, and mass extinction events. However, the impact of these factors on the diversification of the hyperdiverse perching birds (passerines) is unclear because family level relationships are unresolved and the timing of splitting events among lineages is uncertain. We analyzed DNA data from 4,060 nuclear loci and 137 passerine families using concatenation and coalescent approaches to infer a comprehensive phylogenetic hypothesis that clarifies relationships among all passerine families. Then, we calibrated this phylogeny using 13 fossils to examine the effects of different events in Earth history on the timing and rate of passerine diversification. Our analyses reconcile passerine diversification with the fossil and geological records; suggest that passerines originated on the Australian landmass ∼47 Ma; and show that subsequent dispersal and diversification of passerines was affected by a number of climatological and geological events, such as Oligocene glaciation and inundation of the New Zealand landmass. Although passerine diversification rates fluctuated throughout the Cenozoic, we find no link between the rate of passerine diversification and Cenozoic global temperature, and our analyses show that the increases in passerine diversification rate we observe are disconnected from the colonization of new continents. Taken together, these results suggest more complex mechanisms than temperature change or ecological opportunity have controlled macroscale patterns of passerine speciation., Competing Interests: The authors declare no conflict of interest., (Copyright © 2019 the Author(s). Published by PNAS.)

Target enrichment of conserved genomic regions facilitates collecting sequences of many orthologous loci from non-model organisms to address phylogenetic, phylogeographic, population genetic, and molecular evolution questions. Bait sets for sequence capture can simultaneously target thousands of loci, which opens new avenues of research on speciose groups. Current phylogenetic hypotheses on the >103,000 species of Hemiptera have failed to unambiguously resolve major nodes, suggesting that alternative datasets and more thorough taxon sampling may be required to resolve relationships. We use a recently designed ultraconserved element (UCE) bait set for Hemiptera, with a focus on the suborder Heteroptera, or the true bugs, to test previously proposed relationships. We present newly generated UCE data for 36 samples representing three suborders, all seven heteropteran infraorders, 23 families, and 34 genera of Hemiptera and one thysanopteran outgroup. To improve taxon sampling, we also mined additional UCE loci in silico from published hemipteran genomic and transcriptomic data. We obtained 2271 UCE loci for newly sequenced hemipteran taxa, ranging from 265 to 1696 (average 904) per sample. These were similar in number to the data mined from transcriptomes and genomes, but with fewer loci overall. The amount of missing data correlates with greater phylogenetic divergence from taxa used to design the baits. This bait set hybridizes to a wide range of hemipteran taxa and specimens of varying quality, including dried specimens as old as 1973. Our estimated phylogeny yielded topologies consistent with other studies for most nodes and was strongly-supported. We also demonstrate that UCE loci are almost exclusively from the transcribed portion of the genome, thus data can be successfully integrated with existing genomic and transcriptomic resources for more comprehensive phylogenetic sampling, an important feature in the era of phylogenomics. UCE approaches can be used by other researchers for additional studies on hemipteran evolution and other research that requires well resolved phylogenies., (Copyright © 2018 Elsevier Inc. All rights reserved.)

Advances in high-throughput sequencing techniques now allow relatively easy and affordable sequencing of large portions of the genome, even for nonmodel organisms. Many phylogenetic studies reduce costs by focusing their sequencing efforts on a selected set of targeted loci, commonly enriched using sequence capture. The advantage of this approach is that it recovers a consistent set of loci, each with high sequencing depth, which leads to more confidence in the assembly of target sequences. High sequencing depth can also be used to identify phylogenetically informative allelic variation within sequenced individuals, but allele sequences are infrequently assembled in phylogenetic studies. Instead, many scientists perform their phylogenetic analyses using contig sequences which result from the de novo assembly of sequencing reads into contigs containing only canonical nucleobases, and this may reduce both statistical power and phylogenetic accuracy. Here, we develop an easy-to-use pipeline to recover allele sequences from sequence capture data, and we use simulated and empirical data to demonstrate the utility of integrating these allele sequences to analyses performed under the multispecies coalescent model. Our empirical analyses of ultraconserved element locus data collected from the South American hummingbird genus Topaza demonstrate that phased allele sequences carry sufficient phylogenetic information to infer the genetic structure, lineage divergence, and biogeographic history of a genus that diversified during the last 3 myr. The phylogenetic results support the recognition of two species and suggest a high rate of gene flow across large distances of rainforest habitats but rare admixture across the Amazon River. Our simulations provide evidence that analyzing allele sequences leads to more accurate estimates of tree topology and divergence times than the more common approach of using contig sequences.

Conflicts between nuclear and mitochondrial phylogenies have led to uncertainty for some relationships within the tree of life. These conflicts have led some to question the value of mitochondrial DNA in phylogenetics now that genome-scale nuclear data can be readily obtained. However, since mitochondrial DNA is maternally inherited and does not recombine, its phylogeny should be closer to the species tree. Additionally, its rapid evolutionary rate may drive accumulation of mutations along short internodes where relevant information from nuclear loci may be limited. In this study, we examine the mitochondrial phylogeny of Cavitaves to elucidate its congruence with recently published nuclear phylogenies of this group of birds. Cavitaves includes the orders Trogoniformes (trogons), Bucerotiformes (hornbills), Coraciiformes (kingfishers and allies), and Piciformes (woodpeckers and allies). We hypothesized that sparse taxon sampling in previously published mitochondrial trees was responsible for apparent cyto-nuclear discordance. To test this hypothesis, we assembled 27 additional Cavitaves mitogenomes and estimated phylogenies using seven different taxon sampling schemes ranging from five to 42 ingroup species. We also tested the role that partitioning and model choice played in the observed discordance. Our analyses demonstrated that improved taxon sampling could resolve many of the disagreements. Similarly, partitioning was valuable in improving congruence with the topology from nuclear phylogenies, though the model used to generate the mitochondrial phylogenies had less influence. Overall, our results suggest that the mitochondrial tree is trustworthy when partitioning is used with suitable taxon sampling., (Copyright © 2018 Elsevier Inc. All rights reserved.)

The formation of independent evolutionary lineages involves neutral and selective factors, and understanding their relative roles in population divergence is a fundamental goal of speciation research. Correlations between allele frequencies and environmental variability can reveal the role of selection, yet the relative contribution of drift can be difficult to establish. Recently diversified taxa like the Oregon junco (Aves, Passerellidae, Junco hyemalis oreganus) of western North America provide ideal scenarios to apply genetic-environment association analyses (GEA) while controlling for population structure. Analysis of genome-wide SNP loci revealed marked genetic structure consisting of differentiated populations in isolated, dry southern mountain ranges, and less divergent, recently expanded populations in humid northern latitudes. We used correlations between genomic and environmental variance to test for three specific modes of evolutionary divergence: (a) drift in geographic isolation, (b) differentiation along continuous selective gradients and (c) isolation-by-adaptation. We found evidence of strong drift in southern mountains, but also signals of local adaptation driven by temperature, precipitation, elevation and vegetation, especially when controlling for population history. We identified numerous variants under selection scattered across the genome, suggesting that local adaptation can promote rapid differentiation when acting over multiple independent loci., (© 2018 John Wiley & Sons Ltd.)

Using a large, consistent set of loci shared by descent (orthologous) to study relationships among taxa would revolutionize among-lineage comparisons of divergence and speciation processes. Ultraconserved elements (UCEs), highly conserved regions of the genome, offer such genomic markers. The utility of UCEs for deep phylogenetics is clearly established and there are mature analytical frameworks available, but fewer studies apply UCEs to recent evolutionary events, creating a need for additional example datasets and analytical approaches. We used UCEs to study population genomics in snow and McKay's buntings ( Plectrophenax nivalis and P. hyperboreus ). Prior work suggested divergence of these sister species during the last glacial maximum (∼18-74 Kya). With a sequencing depth of ∼30× from four individuals of each species, we used a series of analysis tools to genotype both alleles, obtaining a complete dataset of 2,635 variable loci (∼3.6 single nucleotide polymorphisms/locus) and 796 invariable loci. We found no fixed allelic differences between the lineages, and few loci had large allele frequency differences. Nevertheless, individuals were 100% diagnosable to species, and the two taxa were different genetically ( F ST = 0.034; P = 0.03). The demographic model best fitting the data was one of divergence with gene flow. Estimates of demographic parameters differed from published mtDNA research, with UCE data suggesting lower effective population sizes (∼92,500-240,500 individuals), a deeper divergence time (∼241,000 years), and lower gene flow (2.8-5.2 individuals per generation). Our methods provide a framework for future population studies using UCEs, and our results provide additional evidence that UCEs are useful for answering questions at shallow evolutionary depths., Competing Interests: Brant C. Faircloth is an Academic Editor for PeerJ. The authors declare that they have no other competing interests.

[This corrects the article DOI: 10.1371/journal.pone.0188044.].

Evolutionary transitions between fresh and salt water happen infrequently among bacterioplankton. Within the ubiquitous and highly abundant heterotrophic Alphaproteobacteria order Pelagibacterales (SAR11), most members live in marine habitats, but the LD12 subclade has evolved as a unique freshwater lineage. LD12 cells occur as some of the most dominant freshwater bacterioplankton, yet this group has remained elusive to cultivation, hampering a more thorough understanding of its biology. Here, we report the first successful isolation of an LD12 representative, strain LSUCC0530, using high-throughput dilution-to-extinction cultivation methods, and its complete genome sequence. Growth experiments corroborate ecological data suggesting active populations of LD12 in brackish water up to salinities of ~5. LSUCC0530 has the smallest closed genome thus far reported for a SAR11 strain (1.16 Mbp). The genome affirms many previous metabolic predictions from cultivation-independent analyses, like a complete Embden-Meyerhof-Parnas glycolysis pathway, but also provides novel insights, such as the first isocitrate dehydrogenase in LD12, a likely homologous recombination of malate synthase from outside of the SAR11 clade, and analogous substitutions of ion transporters with others that occur throughout the rest of the SAR11 clade. Growth data support metagenomic recruitment results suggesting temperature-based ecotype diversification within LD12. Key gene losses for osmolyte uptake provide a succinct hypothesis for the evolutionary transition of LD12 from salt to freshwater. For strain LSUCC0530, we propose the provisional nomenclature Candidatus fonsibacter ubiquis.

The Cretaceous-Palaeogene (K-Pg) mass extinction is linked to the rapid emergence of ecologically divergent higher taxa (for example, families and orders) across terrestrial vertebrates, but its impact on the diversification of marine vertebrates is less clear. Spiny-rayed fishes (Acanthomorpha) provide an ideal system for exploring the effects of the K-Pg on fish diversification, yet despite decades of morphological and molecular phylogenetic efforts, resolution of both early diverging lineages and enormously diverse subclades remains problematic. Recent multilocus studies have provided the first resolved phylogenetic backbone for acanthomorphs and suggested novel relationships among major lineages. However, these new relationships and associated timescales have not been interrogated using phylogenomic approaches. Here, we use targeted enrichment of >1,000 ultraconserved elements in conjunction with a divergence time analysis to resolve relationships among 120 major acanthomorph lineages and provide a new timescale for acanthomorph radiation. Our results include a well-supported topology that strongly resolves relationships along the acanthomorph backbone and the recovery of several new relationships within six major percomorph subclades. Divergence time analyses also reveal that crown ages for five of these subclades, and for the bulk of the species diversity in the sixth, coincide with the K-Pg boundary, with divergences between anatomically and ecologically distinctive suprafamilial clades concentrated in the first 10 million years of the Cenozoic.

Anthozoans (e.g., corals, anemones) are an ecologically important and diverse group of marine metazoans that occur from shallow to deep waters worldwide. However, our understanding of the evolutionary relationships among the ~7,500 species within this class is hindered by the lack of phylogenetically informative markers that can be reliably sequenced across a diversity of taxa. We designed and tested 16,306 RNA baits to capture 720 ultraconserved element loci and 1,071 exon loci. Library preparation and target enrichment were performed on 33 taxa from all orders within the class Anthozoa. Following Illumina sequencing and Trinity assembly, we recovered 1,774 of 1,791 targeted loci. The mean number of loci recovered from each species was 638 ± 222, with more loci recovered from octocorals (783 ± 138 loci) than hexacorals (475 ± 187 loci). Parsimony informative sites ranged from 26 to 49% for alignments at differing hierarchical taxonomic levels (e.g., Anthozoa, Octocorallia, Hexacorallia). The per cent of variable sites within each of three genera (Acropora, Alcyonium, and Sinularia) for which multiple species were sequenced ranged from 4.7% to 30%. Maximum-likelihood analyses recovered highly resolved trees with topologies matching those supported by other studies, including the monophyly of the order Scleractinia. Our results demonstrate the utility of this target-enrichment approach to resolve phylogenetic relationships from relatively old to recent divergences. Redesigning the baits with improved affinities to capture loci within each subclass will provide a valuable toolset to address systematic questions, further our understanding of the timing of diversifications and help resolve long-standing controversial relationships in the class Anthozoa., (© 2017 John Wiley & Sons Ltd.)

The rapid diversification of Myotis bats into more than 100 species is one of the most extensive mammalian radiations available for study. Efforts to understand relationships within Myotis have primarily utilized mitochondrial markers and trees inferred from nuclear markers lacked resolution. Our current understanding of relationships within Myotis is therefore biased towards a set of phylogenetic markers that may not reflect the history of the nuclear genome. To resolve this, we sequenced the full mitochondrial genomes of 37 representative Myotis, primarily from the New World, in conjunction with targeted sequencing of 3648 ultraconserved elements (UCEs). We inferred the phylogeny and explored the effects of concatenation and summary phylogenetic methods, as well as combinations of markers based on informativeness or levels of missing data, on our results. Of the 294 phylogenies generated from the nuclear UCE data, all are significantly different from phylogenies inferred using mitochondrial genomes. Even within the nuclear data, quartet frequencies indicate that around half of all UCE loci conflict with the estimated species tree. Several factors can drive such conflict, including incomplete lineage sorting, introgressive hybridization, or even phylogenetic error. Despite the degree of discordance between nuclear UCE loci and the mitochondrial genome and among UCE loci themselves, the most common nuclear topology is recovered in one quarter of all analyses with strong nodal support. Based on these results, we re-examine the evolutionary history of Myotis to better understand the phenomena driving their unique nuclear, mitochondrial, and biogeographic histories.

Weevils (Curculionoidea) comprise one of the most diverse groups of organisms on earth. There is hardly a vascular plant or plant part without its own species of weevil feeding on it and weevil species diversity is greater than the number of fishes, birds, reptiles, amphibians and mammals combined. Here, we employ ultraconserved elements (UCEs) designed for beetles and a novel partitioning strategy of loci to help resolve phylogenetic relationships within the radiation of Australasian smurf-weevils (Eupholini). Despite being emblematic of the New Guinea fauna, no previous phylogenetic studies have been conducted on the Eupholini. In addition to a comprehensive collection of fresh specimens, we supplement our taxon sampling with museum specimens, and this study is the first target enrichment phylogenomic dataset incorporating beetle specimens from museum collections. We use both concatenated and species tree analyses to examine the relationships and taxonomy of this group. For species tree analyses we present a novel partitioning strategy to better model the molecular evolutionary process in UCEs. We found that the current taxonomy is problematic, largely grouping species on the basis of similar color patterns. Finally, our results show that most loci required multiple partitions for nucleotide rate substitution, suggesting that single partitions may not be the optimal partitioning strategy to accommodate rate heterogeneity for UCE loci.

Ostariophysi is a superorder of bony fishes including more than 10,300 species in 1100 genera and 70 families. This superorder is traditionally divided into five major groups (orders): Gonorynchiformes (milkfishes and sandfishes), Cypriniformes (carps and minnows), Characiformes (tetras and their allies), Siluriformes (catfishes), and Gymnotiformes (electric knifefishes). Unambiguous resolution of the relationships among these lineages remains elusive, with previous molecular and morphological analyses failing to produce a consensus phylogeny. In this study, we use over 350 ultraconserved element (UCEs) loci comprising 5 million base pairs collected across 35 representative ostariophysan species to compile one of the most data-rich phylogenies of fishes to date. We use these data to infer higher level (interordinal) relationships among ostariophysan fishes, focusing on the monophyly of the Characiformes-one of the most contentiously debated groups in fish systematics. As with most previous molecular studies, we recover a non-monophyletic Characiformes with the two monophyletic suborders, Citharinoidei and Characoidei, more closely related to other ostariophysan clades than to each other. We also explore incongruence between results from different UCE data sets, issues of orthology, and the use of morphological characters in combination with our molecular data. [Conserved sequence; ichthyology; massively parallel sequencing; morphology; next-generation sequencing; UCEs.]., (© The Author(s) 2017. Published by Oxford University Press, on behalf of the Society of Systematic Biologists. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

Melanophryniscus is a bufonid frog genus with a broad geographic distribution over southeastern South America. In recent years, several new species of Melanophryniscus have been discovered in southern Brazil showing a distinctive life-history strategy for the genus - breeding in phytotelmata - as well as a strong association with high-altitude regions. In this study, we use mitogenomic data to infer the phylogenetic relationships among diurnal, phytotelm-breeding Melanophryniscus and to determine the timing of their divergence. We obtained the mitochondrial genomes (not including the control region) for eight individuals of Melanophryniscus representing all three described species (M. alipioi, M. milanoi, and M. xanthostomus), as well as some recently-discovered and potentially new species. Gene order was conserved in all species and corresponded to the general order found in bufonids. Although the phylogenetic relationships among the studied species was poorly supported, dating confirmed that they diverged during the Pleistocene, suggesting that phytotelm breeding could have arisen during drier periods in the glacial/interglacial cycles due to a decrease in the availability of permanent streams or ephemeral/temporary streams or ponds in which Melanophryniscus species commonly breed., (Copyright © 2017 Elsevier B.V. All rights reserved.)

Adaptive radiations could often occur in discrete stages. For instance, the species flock of ∼1000 species of Lake Malawi cichlid fishes might have only diverged once between rocky and sandy environments during the initial stage of their diversification. All further diversification within the rock-dwelling (mbuna) or sand-dwelling (utaka) cichlids would have occurred during a subsequent second stage of extensive trophic evolution that was followed by a third stage of sexual trait divergence. We provide an improved phylogenetic framework for Malawi cichlids to test this three-stage hypothesis based on newly reconstructed phylogenetic relationships among 32 taxonomically disparate Malawi cichlids species. Using several reconstruction methods and 1037 ultra-conserved element (UCE) markers, we recovered a molecular phylogeny that confidently resolved relationships among most of the Malawi lineages sampled when a bifurcating framework was enforced. These bifurcating reconstructions also indicated that the sand-dwelling species Cyathochromis obliquidens was well-nested within the primarily rock-dwelling radiation known as the mbuna. In contrast to predictions from the three-stage model of vertebrate diversification, the recovered phylogeny reveals an initial colonization of rocky reefs, followed by substantial diversification of rock-dwelling lineages, and then at least one instance of subsequent evolution back into sandy habitats. This repeated evolution into major habitat types provides further evidence that the three-stage model of Malawi cichlid diversification has numerous exceptions., (Copyright © 2017 Elsevier Inc. All rights reserved.)