Robber or assassin flies (Asilidae) are a diverse family of venomous predators. The most recent classification organizes Asilidae into 14 subfamilies based on a morphological phylogeny, but many of these are not supported by molecular data. To test the monophyly of various clades in Asilidae, we used the recently developed Diptera-wide ultraconserved element bait set to compile seven datasets comprising 151 robber flies and 145–2496 loci. We also compared the performance of various nodal support metrics. Our Maximum Likelihood phylogeny was fully resolved and well supported, but partially incongruent with the coalescent phylogeny. Further examination of datasets suggested that GC bias had influenced gene tree inference and subsequent species tree analyses. The subfamilies Brachyrhopalinae, Dasypogoninae, Dioctriinae, Stenopogoninae, Tillobromatinae, Trigonomiminae and Willistonininae were not recovered as monophyletic. The inter-subfamily relationships are summarized as follows: Laphriinae and Dioctriinae (in part) successively sister to the remaining subfamilies, which form two clades: the first consists of a grade of Stenopogoninae (in part), Willistonininae (in part), Bathypogoninae+Phellinae, Stichopogoninae, Leptogastrinae, Ommatiinae and Asilinae; the second clade consists of a paraphyletic assemblage of genera from Dioctriinae (in part), Trigonomiminae, Stenopogoninae (in part), Tillobromatinae, Brachyrhopalinae and Dasypogoninae. This phylogeny demonstrates that the higher classification of Asilidae is far from settled, but does provide a much-needed foundation for a thorough revision of the subfamily classification.
One of the most important non-Apis groups of bees for agriculture is the mason bee subgenus Osmia Panzer (Osmia), or Osmia s.s. (Hymenoptera: Megachilidae). Out of the 29 known species, four have been developed as managed pollinators of orchards. In addition, the group is important as a source of non-native pollinators, given that several species have been introduced into new areas. Osmia s.s. occurs naturally throughout the northern temperate zone with greatest species richness in Europe and Asia. Here, we integrate phylogenomic data from ultraconserved elements (UCEs), near complete taxon sampling, and a diversity of analytical approaches to infer the phylogeny, divergence times and biogeographic history of Osmia s.s. We also demonstrate how mitochondrial sequence data can be extracted from ultraconserved element data and combined with sequences from public repositories in order to test the phylogeny, examine species boundaries and identify specimen-associated, non-bee DNA. We resolve the phylogeny of Osmia s.s. and show strong support that Nearctic Osmia ribifloris is the sister group to the rest of the subgenus. Biogeographic analyses indicate that the group originated during the Late Miocene in the West Nearctic plus East Palearctic region following dispersal from the East Palearctic to the West Nearctic across the Bering land bridge prior to its closure 5.5–4.8 Ma. The mitochondrial DNA results reveal potential taxonomic synonymies involving Osmia yanbianensis and Osmia opima, and Osmia rufina, Osmia rufinoides and Osmia taurus.
PREMISE Recent, rapid radiations present a challenge for phylogenetic reconstruction. Fast successive speciation events typically lead to low sequence divergence and poorly resolved relationships with standard phylogenetic markers. Target sequence capture of many independent nuclear loci has the potential to improve phylogenetic resolution for rapid radiations. METHODS Here we applied target sequence capture with 353 protein-coding genes (Angiosperms353 bait kit) to Veronica sect. Hebe (common name hebe) to determine its utility for improving the phylogenetic resolution of rapid radiations. Veronica section Hebe originated 5–10 million years ago in New Zealand, forming a monophyletic radiation of ca 130 extant species. RESULTS We obtained approximately 150 kbp of 353 protein-coding exons and an additional 200 kbp of flanking noncoding sequences for each of 77 hebe and two outgroup species. When comparing coding, noncoding, and combined data sets, we found that the latter provided the best overall phylogenetic resolution. While some deep nodes in the radiation remained unresolved, our phylogeny provided broad and often improved support for subclades identified by both morphology and standard markers in previous studies. Gene-tree discordance was nonetheless widespread, indicating that additional methods are needed to disentangle fully the history of the radiation. CONCLUSIONS Phylogenomic target capture data sets both increase phylogenetic signal and deliver new insights into the complex evolutionary history of rapid radiations as compared with traditional markers. Improving methods to resolve remaining discordance among loci from target sequence capture is now important to facilitate the further study of rapid radiations.
The genetic consequences of species-wide declines are rarely quantified because the timing and extent of the decline varies across the species’ range. The sea otter (Enhydra lutris) is a unique model in this regard. Their dramatic decline from thousands to fewer than 100 individuals per population occurred range-wide and nearly simultaneously due to the 18th-19th century fur trade. Consequently, each sea otter population represents an independent natural experiment of recovery after extreme population decline. We designed sequence capture probes for 50 megabases of sea otter exonic and neutral genomic regions. We sequenced 107 sea otters from five populations that span the species range to high coverage (18-76x) and three historic Californian samples from 1500 and 200 years ago to low coverage (1.5-3.5X). We observe distinct population structure and find that sea otters in California are the last survivors of a divergent lineage isolated for thousands of years and therefore warrant special conservation concern. We detect signals of extreme population decline in every surviving sea otter population and use this demographic history to design forward-in-time simulations of coding sequence. Our simulations indicate that this decline could lower the fitness of recovering populations for generations. However, the simulations also demonstrate how historically low effective population sizes prior to the fur trade may have mitigated the effects of population decline on genetic health. Our comprehensive approach shows how demographic inference from genomic data, coupled with simulations, allows assessment of extinction risk and different models of recovery.
Until recently many historical museum specimens were largely inaccessible to genomic inquiry, but high-throughput sequencing (HTS) approaches have allowed researchers to successfully sequence genomic DNA from dried and fluid-preserved museum specimens. In addition to preserved specimens, many museums contain large series of allozyme supernatant samples but the amenability of these samples to HTS has not yet been assessed. Here, we compared the performance of a target-capture approach using alternative sources of genomic DNA from ten specimens of spring salamanders (Plethodontidae: Gyrinophilus porphyriticus) collected between 1985 and 1990: allozyme supernatants, allozyme homogenate pellets, and formalin-fixed tissues. We designed capture probes based on double-digest restriction-site associated sequencing (RADseq) derived loci from frozen blood samples available for seven of the specimens and assessed the success and consistency of capture and RADseq approaches. This study design enabled direct comparisons of data quality and potential biases among the different datasets for phylogenomic and population genomic analyses. We found that in phylogenetic analyses, all enrichment types for a given specimen clustered together. In principal component space all capture-based samples clustered together, but RADseq samples did not cluster with corresponding capture-based samples. SNP calls were on average 18.3% different between enrichment types for a given individual, but these discrepancies were primarily due to differences in heterozygous/homozygous SNP calls. We demonstrate that both allozyme supernatant and formalin-fixed samples can be successfully used for population genomic analyses and we discuss ways to identify and reduce biases associated with combining capture and RADseq data.
The octocoral genus Chrysogorgia Duchassaing & Michelotti, 1864 contains about 80 nominal species that are ecologically important components of benthic communities. A review of morphological variation within Chrysogorgia identified nine groups that may represent distinct lineages; which would help clarify their taxonomy and guide genus revision. Here, we applied a recently developed universal target enrichment bait method for octocoral exons and ultraconserved elements (UCEs) on 96 specimens varying in morphology, collection ages and DNA quality and quantity to determine whether there was genetic support for morphologically defined groups within Chrysogorgia. Following Illumina sequencing and SPAdes assembly we recovered 1682 of 1700 targeted exon loci and 1333 of 1340 targeted UCE loci. Loci recovery per sample was highly variable and significantly correlated with time since specimen collection (2-60 years) and DNA quantity and quality. Phylogenetically informative sites in UCE and exon loci ranged from 34.75 – 36.10% for 50% and 75% taxon-occupancy matrices, respectively. Maximum likelihood analyses recovered highly resolved trees with topologies supporting the recognition of 11 candidate genera, nine of which are novel. Our results also demonstrate that this target-enrichment approach can be applied to degraded museum specimens of up to 60 years old. This study shows that an integrative approach consisting of molecular and morphological methods is essential in a proper revision of Chrysogorgiath taxonomy and regional diversity of these ecologically important corals.
Crop brassicas include three diploid [ Brassica rapa (AA; 2 n = 2 x = 16), B. nigra (BB; 2 n = 2 x = 18), and B. oleracea (CC; 2 n = 2 x = 20)] and three derived allotetraploid species. It is difficult to distinguish Brassica chromosomes as they are small and morphologically similar. We aimed to develop a genome-sequence based cytogenetic toolkit for reproducible identification of Brassica chromosomes and their structural variations. A bioinformatic pipeline was used to extract repeat-free sequences from the whole genome assembly of B. rapa . Identified sequences were subsequently used to develop four c. 47-mer oligonucleotide libraries comprising 27,100, 11,084, 9,291, and 16,312 oligonucleotides. We selected these oligonucleotides after removing repeats from 18 identified sites (500–1,000 kb) with 1,997–5,420 oligonucleotides localized at each site in B. rapa . For one set of probes, a new method for amplification or immortalization of the library is described. oligonucleotide probes produced specific and reproducible in situ hybridization patterns for all chromosomes belonging to A, B, C, and R ( Raphanus sativu s) genomes. The probes were able to identify structural changes between the genomes, including translocations, fusions, and deletions. Furthermore, the probes were able to identify a structural translocation between a pak choi and turnip cultivar of B. rapa. Overall, the comparative chromosomal mapping helps understand the role of chromosome structural changes during genome evolution and speciation in the family Brassicaceae. The probes can also be used to identify chromosomes in aneuploids such as addition lines used for gene mapping, and to track transfer of chromosomes in hybridization and breeding programs.
Ann Arbor, MI 48103
(d/b/a Daicel Arbor Biosciences)
All Rights Reserved.