skip to main content


Title: Target gene enrichment in the cyclophyllidean cestodes, the most diverse group of tapeworms
Abstract

The Cyclophyllidea is the most diverse order of tapeworms, encompassing species that infect all classes of terrestrial tetrapods including humans and domesticated animals. Available phylogenetic reconstructions based either on morphology or molecular data lack the resolution to allow scientists to either propose a solid taxonomy or infer evolutionary associations. Molecular markers available for the Cyclophyllidea mostly include ribosomalDNAand mitochondrial loci. In this study, we identified 3641 single‐copy nuclear coding loci by comparing the genomes ofHymenolepis microstoma,Echinococcus granulosusandTaenia solium. We designedRNAbaits based on the sequence ofH. microstoma, and applied target enrichment and Illumina sequencing to test the utility of those baits to recover loci useful for phylogenetic analyses. We capturedDNAfrom five species of tapeworms representing two families of cyclophyllideans. We obtained an average of 3284 (90%) of the targets from the test samples and then used captured sequences (2 181 361 bp in total; fragment size ranging from 301 to 6969 bp) to reconstruct a phylogeny for the five test species plus the three species for which genomic data are available. The results were consistent with the current consensus regarding cyclophyllidean relationships. To assess the potential for our method to yield informative genetic variation at intraspecific scales, we extracted 14 074 single nucleotide polymorphisms (SNPs) from alignments of fourArostrilepis macrocirrosaand twoA. cookiand successfully inferred their relationships. The results showed that our target gene tools yield data sets that provide robust inferences at a range of taxonomic scales in the Cyclophyllidea.

 
more » « less
NSF-PAR ID:
10243958
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  
Publisher / Repository:
Wiley-Blackwell
Date Published:
Journal Name:
Molecular Ecology Resources
Volume:
16
Issue:
5
ISSN:
1755-098X
Page Range / eLocation ID:
p. 1095-1106
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Molecular ecologists seek to genotype hundreds to thousands of loci from hundreds to thousands of individuals at minimal cost per sample. Current methods, such as restriction‐site‐associatedDNAsequencing (RADseq) and sequence capture, are constrained by costs associated with inefficient use of sequencing data and sample preparation. Here, we introduceRADcap, an approach that combines the major benefits ofRADseq (low cost with specific start positions) with those of sequence capture (repeatable sequencing of specific loci) to significantly increase efficiency and reduce costs relative to current approaches.RADcap uses a new version of dual‐digestRADseq (3RAD) to identify candidateSNPloci for capture bait design and subsequently uses custom sequence capture baits to consistently enrich candidateSNPloci across many individuals. We combined this approach with a new library preparation method for identifying and removingPCRduplicates from 3RADlibraries, which allows researchers to processRADseq data using traditional pipelines, and we tested theRADcap method by genotyping sets of 96–384Wisteriaplants. Our results demonstrate that ourRADcap method: (i) methodologically reduces (to <5%) and allows computational removal ofPCRduplicate reads from data, (ii) achieves 80–90% reads on target in 11 of 12 enrichments, (iii) returns consistent coverage (≥4×) across >90% of individuals at up to 99.8% of the targeted loci, (iv) produces consistently high occupancy matrices of genotypes across hundreds of individuals and (v) costs significantly less than current approaches.

     
    more » « less
  2. Abstract Aim

    The Lesser Sunda Islands are situated between the Sunda and Sahul Shelves, with a linear arrangement that has functioned as a two‐way filter for taxa dispersing between the Asian and Australo‐Papuan biogeographical realms. Distributional patterns of many terrestrial vertebrates suggest a stepping‐stone model of island colonization. Here we investigate the timing and sequence of island colonization in Asian‐origin fanged frogs from the volcanic Sunda Arc islands with the goal of testing the stepping‐stone model of island colonization.

    Location

    The Indonesian islands of Java, Lombok, Sumbawa, Flores and Lembata.

    Taxon

    Limnonectes dammermaniandL. kadarsani(Family: Dicroglossidae)

    Methods

    MitochondrialDNAwas sequenced from 153 frogs to identify major lineages and to select samples for an exon‐capture experiment. We designed probes to capture sequence data from 974 exonic loci (1,235,981 bp) from 48 frogs including the outgroup species,L. microdiscus. The resulting data were analysed using phylogenetic, population genetic and biogeographical model testing methods.

    Results

    The mtDNAphylogeny findsL. kadarsaniparaphyletic with respect toL. dammermani, with a pectinate topology consistent with the stepping‐stone model. Phylogenomic analyses of 974 exons recovered the two species as monophyletic sister taxa that diverged ~7.6 Ma with no detectable contemporary gene flow, suggesting introgression of theL. dammermanimitochondrion intoL. kadarsanion Lombok resulting from an isolated ancient hybridization event ~4 Ma. WithinL. kadarsani,the Lombok lineage diverged first while the Sumbawa and Lembata lineages are nested within a Flores assemblage composed of two parapatrically distributed lineages meeting in central Flores. Biogeographical model comparison found strict stepping‐stone dispersal to be less likely than models involving leap‐frog dispersal events.

    Main conclusions

    These results suggest that the currently accepted stepping‐stone model of island colonization might not best explain the current patterns of diversity in the archipelago. The high degree of genetic structure, large divergence times, and absent or low levels of migration between lineages suggests thatL. kadarsanirepresents five distinct species.

     
    more » « less
  3. Abstract

    Custom sequence capture experiments are becoming an efficient approach for gathering large sets of orthologous markers in nonmodel organisms. Transcriptome‐based exon capture utilizes transcript sequences to design capture probes, typically using a reference genome to identify intron–exon boundaries to exclude shorter exons (<200 bp). Here, we test directly using transcript sequences for probe design, which are often composed of multiple exons of varying lengths. Using 1260 orthologous transcripts, we conducted sequence captures across multiple phylogenetic scales for frogs, including outgroups ~100 Myr divergent from the ingroup. We recovered a large phylogenomic data set consisting of sequence alignments for 1047 of the 1260 transcriptome‐based loci (~561 000 bp) and a large quantity of highly variable regions flanking the exons in transcripts (~70 000 bp), the latter improving substantially by only including ingroup species (~797 000 bp). We recovered both shorter (<100 bp) and longer exons (>200 bp), with no major reduction in coverage towards the ends of exons. We observed significant differences in the performance of blocking oligos for target enrichment and nontarget depletion during captures, and differences inPCRduplication rates resulting from the number of individuals pooled for capture reactions. We explicitly tested the effects of phylogenetic distance on capture sensitivity, specificity, and missing data, and provide a baseline estimate of expectations for these metrics based on a priori knowledge of nuclear pairwise differences among samples. We provide recommendations for transcriptome‐based exon capture design based on our results, cost estimates and offer multiple pipelines for data assembly and analysis.

     
    more » « less
  4. Abstract Aim

    To test the importance of alternative diversification drivers and biogeographical processes for the evolution of Amazonian upland forest birds through a densely sampled analysis of diversification of the endemic Amazonian genusRhegmatorhinaat multiple taxonomic and temporal scales.

    Location

    Amazonia.

    Taxon

    Antbirds (Thamnophilidae).

    Methods

    We sequenced four mtDNAand nuclear gene regions of 120 individuals from 50 localities representing all recognized species and subspecies of the genus. We performed molecular phylogenetic analyses using both gene tree and species tree methods, molecular dating analysis and estimated population demographic history and gene flow.

    Results

    Dense sampling throughout the distribution ofRhegmatorhinarevealed that the main Amazonian rivers delimit the geographic distribution of taxa as inferred from mtDNAlineages. Molecular phylogenetic analyses resulted in a strongly supported phylogenetic hypothesis for the genus, with two main clades currently separated by the Madeira River. Molecular dating analysis indicated diversification during the Quaternary. Reconstruction of recent demographic history of populations revealed a trend for population expansion in eastern Amazonia and stability in the west. Estimates of gene flow corroborate the possibility that migration after divergence had some influence on the current patterns of diversity.

    Main Conclusions

    Based on broad‐scale sampling, a clarification of taxonomic boundaries, and strongly supported phylogenetic relationships, we confirm that, first, mitochondrial lineages within this upland forest Amazonian bird genus agree with spatial patterns known for decades based on phenotypes, and second, that most lineages are geographically delimited by the large Amazonian rivers. The association between past demographic changes related to palaeoclimatic cycles and the historically varying strength and size of rivers as barriers to dispersal may be the path to the answer to the long‐standing question of identifying the main drivers of Amazonian diversification.

     
    more » « less
  5. Abstract

    Evolutionary biology often seeks to decipher the drivers of speciation, and much debate persists over the relative importance of isolation and gene flow in the formation of new species. Genetic studies of closely related species can assess if gene flow was present during speciation, because signatures of past introgression often persist in the genome. We test hypotheses on which mechanisms of speciation drove diversity among three distinct lineages of desert tortoise in the genusGopherus. These lineages offer a powerful system to study speciation, because different biogeographic patterns (physical vs. ecological segregation) are observed at opposing ends of their distributions. We use 82 samples collected from 38 sites, representing the entire species' distribution and generate sequence data for mtDNAand four nuclear loci. A multilocus phylogenetic analysis in *BEAST estimates the species tree.RNA‐seq data yield 20,126 synonymous variants from 7665 contigs from two individuals of each of the three lineages. Analyses of these data using the demographic inference package ∂a∂i serve to test the null hypothesis of no gene flow during divergence. The best‐fit demographic model for the three taxa is concordant with the *BEAST species tree, and the ∂a∂i analysis does not indicate gene flow among any of the three lineages during their divergence. These analyses suggest that divergence among the lineages occurred in the absence of gene flow and in this scenario the genetic signature of ecological isolation (parapatric model) cannot be differentiated from geographic isolation (allopatric model).

     
    more » « less