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

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

Low hybridization temperatures improve target capture success of invertebrate loci: a case study of leaf-footed bugs (Hemiptera: Coreoidea).

Target capture is widely used in phylogenomic, ecological and functional genomic studies. Bait sets that allow capture from a diversity of species can be advantageous, but high-sequence divergence from baits can limit yields. Currently, only four experimental comparisons of a critical target capture parameter, hybridization temperature, have been published. These have been in vertebrates, where bait divergences are typically low, and none include invertebrates where bait-target divergences may be higher. Most invertebrate capture studies use a fixed, high hybridization temperature to maximize the proportion of on-target data, but many report low locus recovery. Using leaf-footed bugs (Hemiptera: Coreoidea), we investigate the effect of hybridization temperature on capture success of ultraconserved elements targeted by (i) baits developed from divergent hemipteran genomes and (ii) baits developed from less divergent coreoid transcriptomes. Lower temperatures generally resulted in more contigs and improved recovery of targets despite a lower proportion of on-target reads, lower read depth and more putative paralogues. Hybridization temperatures had less of an effect when using transcriptome-derived baits, which is probably due to lower bait-target divergences and greater bait tiling density. Thus, accommodating low hybridization temperatures during target capture can provide a cost-effective, widely applicable solution to improve invertebrate locus recovery.

New alignment-based sequence extraction software (ALiBaSeq) and its utility for deep level phylogenetics.

Despite many bioinformatic solutions for analyzing sequencing data, few options exist for targeted sequence retrieval from whole genomic sequencing (WGS) data with the ultimate goal of generating a phylogeny. Available tools especially struggle at deep phylogenetic levels and necessitate amino-acid space searches, which may increase rates of false positive results. Many tools are also difficult to install and may lack adequate user resources. Here, we describe a program that uses freely available similarity search tools to find homologs in assembled WGS data with unparalleled freedom to modify parameters. We evaluate its performance compared to other commonly used bioinformatics tools on two divergent insect species (>200 My) for which annotated genomes exist, and on one large set each of highly conserved and more variable loci. Our software is capable of retrieving orthologs from well-curated or unannotated, low or high depth shotgun, and target capture assemblies as well or better than other software as assessed by recovering the most genes with maximal coverage and with a low rate of false positives throughout all datasets. When assessing this combination of criteria, ALiBaSeq is frequently the best evaluated tool for gathering the most comprehensive and accurate phylogenetic alignments on all types of data tested. The software (implemented in Python), tutorials, and manual are freely available at https://github.com/AlexKnyshov/alibaseq.

Insight from an ultraconserved element bait set designed for hemipteran phylogenetics integrated with genomic resources.

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.

Sylvatic host associations of Triatominae and implications for Chagas disease reservoirs: a review and new host records based on archival specimens.

BACKGROUND: The 152 extant species of kissing bug include important vectors of the debilitating, chronic, and often fatal Chagas disease, which affects several million people mainly in Central and South America. An understanding of the natural hosts of this speciose group of blood-feeding insects has and will continue to aid ongoing efforts to impede the spread of Chagas disease. However, information on kissing bug biology is piecemeal and scattered, developed using methods with varying levels of accuracy over more than 100 years. Existing host records are heavily biased towards well-studied primary vector species and are derived from primarily three different types of observations, associational, immunological or DNA-based, with varying reliability. METHODS: We gather a comprehensive and unparalleled number of sources reporting host associations via rigorous targeted searches of publication databases to review all known natural, or sylvatic, host records including information on how each record was collected. We integrate this information with novel host records obtained via attempted amplification and sequencing of a ∼160 base pair (bp) region of the vertebrate 12S mitochondrial gene from the gastrointestinal tract of 64 archival specimens of Triatominae representing 19 species collected primarily in sylvatic habitats throughout the southern United States and Central and South America during the past 10 years. We show the utility of this method for uncovering novel and under-studied groups of Triatominae hosts, as well as detecting the presence of the Chagas disease pathogen via Polymerase Chain Reaction (PCR) of a ∼400 bp sequence of the trypanosome 18S gene. RESULTS: New host associations for several groups of arboreal mammals were determined including sloths, New World monkeys, coatis, arboreal porcupines and, for the first time as a host of any Triatominae, tayras. A thorough review of previously documented sylvatic hosts, organized by triatomine species and the type of observation (associational, antibody-based, or DNA-based), is presented in a phylogenetic context and highlights large gaps in our knowledge of Triatominae biology. CONCLUSION: The application of DNA-based methods of host identification towards additional species of Triatominae, including rarely collected species that may require use of archival specimens, is the most efficient and promising way to resolve recognized shortfalls.

Evolution of the assassin's arms: insights from a phylogeny of combined transcriptomic and ribosomal DNA data (Heteroptera: Reduvioidea).

Assassin bugs (Reduvioidea) are one of the most diverse (>7,000 spp.) lineages of predatory animals and have evolved an astounding diversity of raptorial leg modifications for handling prey. The evolution of these modifications is not well understood due to the lack of a robust phylogeny, especially at deeper nodes. We here utilize refined data from transcriptomes (370 loci) to stabilize the backbone phylogeny of Reduvioidea, revealing the position of major clades (e.g., the Chagas disease vectors Triatominae). Analyses combining transcriptomic and Sanger-sequencing datasets result in the first well-resolved phylogeny of Reduvioidea. Despite amounts of missing data, the transcriptomic loci resolve deeper nodes while the targeted ribosomal genes anchor taxa at shallower nodes, both with high support. This phylogeny reveals patterns of raptorial leg evolution across major leg types. Hairy attachment structures (fossula spongiosa), present in the ancestor of Reduvioidea, were lost multiple times within the clade. In contrast to prior hypotheses, this loss is not directly correlated with the evolution of alternative raptorial leg types. Our results suggest that prey type, predatory behavior, salivary toxicity, and morphological adaptations pose intricate and interrelated factors influencing the evolution of this diverse group of predators.

Efficient capture of natural history data reveals prey conservatism of cryptic termite predators.

Stenophagy, specialization of a clade on a narrow range of taxa, has not been well studied in speciose clades of predators, principally due to the difficulty of obtaining adequate natural history data. The pantropical Salyavatinae (Hemiptera: Reduviidae; 17 genera, 107 species) contains members with enigmatic morphology and specialized behavior for feeding on termites. All Salyavatinae are suspected specialist termite predators; however, existing observations are limited to seven species. Prior analyses indicate that Salyavatinae may be paraphyletic with respect to another subfamily, Sphaeridopinae, also hypothesized to feed on termites. A molecular phylogeny of these putative termite assassins is here constructed using seven loci from 28 species in nine genera and is used in a dating analysis to shed light on the timing of Neotropical colonization by this primarily Old World clade. DNA extracted from gut contents of 50 individuals was assayed using PCR with prey-specific primers.Molecular assays, along with recent photographs and observations, provide substantial evidence that this clade feeds specifically upon termites, documenting 28 new individual associations. Our phylogeny supports a sister group relationship of the Neotropical genus Salyavata with Sphaeridopinae. Termite association data combined with our phylogeny provide evidence of previously unknown prey conservatism among clades of one of the most diverse groups of specialist termite predators.

A rapid method for assaying thiaminase I activity in diverse biological samples.

Vitamin B1 (thiamine) deficiencies can lead to neurological disorders, reproductive failure and death in wild and domestic animal populations. In some cases, disease is brought about by the consumption of foods high in thiaminase I activity. Levels of thiaminase activity in these foods are highly variable and the factors leading to production of this enzyme are poorly understood. Here we describe improvements in a spectrophotometric thiaminase I activity assay that measures the disappearance of 4-nitrothiophenol, a favored nucleophile co-substrate that replaces the thiazole portion of thiamine during the inactivation of thiamine by the enzyme. Scalable sample processing protocols and a 96-well microtiter plate format are presented that allow the rapid evaluation of multiple, replicated samples in the course of only a few hours. Observed levels of activity in bacterial culture supernatant, fish, ferns and molluscs using this colorimetric assay were similar to previously published reports that employed a radiometric method. Organisms devoid of thiaminase I, based upon previous work, showed no activity with this assay. In addition, activity was found in a variety of fishes and one fern species from which this enzyme had not previously been reported. Overall, we demonstrate the suitability of this technique for measuring thiaminase I activity within small amounts of tissue and environmental samples with replication levels that were heretofore prohibitive. The assay provides a considerable improvement in the ability to examine and understand the properties of an enzyme that has a substantial influence on organism and ecosystem health.