Major Revisions in Pancrustacean Phylogeny and Evidence of Sensitivity to Taxon Sampling.

The clade Pancrustacea, comprising crustaceans and hexapods, is the most diverse group of animals on earth, containing over 80% of animal species and half of animal biomass. It has been the subject of several recent phylogenomic analyses, yet relationships within Pancrustacea show a notable lack of stability. Here, the phylogeny is estimated with expanded taxon sampling, particularly of malacostracans. We show small changes in taxon sampling have large impacts on phylogenetic estimation. By analyzing identical orthologs between two slightly different taxon sets, we show that the differences in the resulting topologies are due primarily to the effects of taxon sampling on the phylogenetic reconstruction method. We compare trees resulting from our phylogenomic analyses with those from the literature to explore the large tree space of pancrustacean phylogenetic hypotheses and find that statistical topology tests reject the previously published trees in favor of the maximum likelihood trees produced here. Our results reject several clades including Caridoida, Eucarida, Multicrustacea, Vericrustacea, and Syncarida. Notably, we find Copepoda nested within Allotriocarida with high support and recover a novel relationship between decapods, euphausiids, and syncarids that we refer to as the Syneucarida. With denser taxon sampling, we find Stomatopoda sister to this latter clade, which we collectively name Stomatocarida, dividing Malacostraca into three clades: Leptostraca, Peracarida, and Stomatocarida. A new Bayesian divergence time estimation is conducted using 13 vetted fossils. We review our results in the context of other pancrustacean phylogenetic hypotheses and highlight 15 key taxa to sample in future studies.

Detecting and Removing Sample Contamination in Phylogenomic Data: An Example and its Implications for Cicadidae Phylogeny (Insecta: Hemiptera).

Contamination of a genetic sample with DNA from one or more nontarget species is a continuing concern of molecular phylogenetic studies, both Sanger sequencing studies and next-generation sequencing studies. We developed an automated pipeline for identifying and excluding likely cross-contaminated loci based on the detection of bimodal distributions of patristic distances across gene trees. When contamination occurs between samples within a data set, a comparison between a contaminated sample and its contaminant taxon will yield bimodal distributions with one peak close to zero patristic distance. This new method does not rely on a priori knowledge of taxon relatedness nor does it determine the causes(s) of the contamination. Exclusion of putatively contaminated loci from a data set generated for the insect family Cicadidae showed that these sequences were affecting some topological patterns and branch supports, although the effects were sometimes subtle, with some contamination-influenced relationships exhibiting strong bootstrap support. Long tip branches and outlier values for one anchored phylogenomic pipeline statistic (AvgNHomologs) were correlated with the presence of contamination. While the anchored hybrid enrichment markers used here, which target hemipteroid taxa, proved effective in resolving deep and shallow level Cicadidae relationships in aggregate, individual markers contained inadequate phylogenetic signal, in part probably due to short length. The cleaned data set, consisting of 429 loci, from 90 genera representing 44 of 56 current Cicadidae tribes, supported three of the four sampled Cicadidae subfamilies in concatenated-matrix maximum likelihood (ML) and multispecies coalescent-based species tree analyses, with the fourth subfamily weakly supported in the ML trees. No well-supported patterns from previous family-level Sanger sequencing studies of Cicadidae phylogeny were contradicted. One taxon (Aragualna plenalinea) did not fall with its current subfamily in the genetic tree, and this genus and its tribe Aragualnini is reclassified to Tibicininae following morphological re-examination. Only subtle differences were observed in trees after the removal of loci for which divergent base frequencies were detected. Greater success may be achieved by increased taxon sampling and developing a probe set targeting a more recent common ancestor and longer loci. Searches for contamination are an essential step in phylogenomic analyses of all kinds and our pipeline is an effective solution. [Auchenorrhyncha; base-composition bias; Cicadidae; Cicadoidea; Hemiptera; phylogenetic conflict.].

Molecular dating and viral load growth rates suggested that the eclipse phase lasted about a week in HIV-1 infected adults in East Africa and Thailand.

Most HIV-1 infected individuals do not know their infection dates. Precise infection timing is crucial information for studies that document transmission networks or drug levels at infection. To improve infection timing, we used the prospective RV217 cohort where the window when plasma viremia becomes detectable is narrow: the last negative visit occurred a median of four days before the first detectable HIV-1 viremia with an RNA test, referred below as diagnosis. We sequenced 1,280 HIV-1 genomes from 39 participants at a median of 4, 32 and 170 days post-diagnosis. HIV-1 infections were dated by using sequence-based methods and a viral load regression method. Bayesian coalescent and viral load regression estimated that infections occurred a median of 6 days prior to diagnosis (IQR: 9-3 and 11-4 days prior, respectively). Poisson-Fitter, which analyzes the distribution of hamming distances among sequences, estimated a median of 7 days prior to diagnosis (IQR: 15-4 days) based on sequences sampled 4 days post-diagnosis, but it did not yield plausible results using sequences sampled at 32 days. Fourteen participants reported a high-risk exposure event at a median of 8 days prior to diagnosis (IQR: 12 to 6 days prior). These different methods concurred that HIV-1 infection occurred about a week before detectable viremia, corresponding to 20 days (IQR: 34-15 days) before peak viral load. Together, our methods comparison helps define a framework for future dating studies in early HIV-1 infection.

Factors influencing estimates of HIV-1 infection timing using BEAST.

While large datasets of HIV-1 sequences are increasingly being generated, many studies rely on a single gene or fragment of the genome and few comparative studies across genes have been done. We performed genome-based and gene-specific Bayesian phylogenetic analyses to investigate how certain factors impact estimates of the infection dates in an acute HIV-1 infection cohort, RV217. In this cohort, HIV-1 diagnosis corresponded to the first RNA positive test and occurred a median of four days after the last negative test, allowing us to compare timing estimates using BEAST to a narrow window of infection. We analyzed HIV-1 sequences sampled one week, one month and six months after HIV-1 diagnosis in 39 individuals. We found that shared diversity and temporal signal was limited in acute infection, and insufficient to allow timing inferences in the shortest HIV-1 genes, thus dated phylogenies were primarily analyzed for env, gag, pol and near full-length genomes. There was no one best-fitting model across participants and genes, though relaxed molecular clocks (73% of best-fitting models) and the Bayesian skyline (49%) tended to be favored. For infections with single founders, the infection date was estimated to be around one week pre-diagnosis for env (IQR: 3-9 days) and gag (IQR: 5-9 days), whilst the genome placed it at a median of 10 days (IQR: 4-19). Multiply-founded infections proved problematic to date. Our ability to compare timing inferences to precise estimates of HIV-1 infection (within a week) highlights that molecular dating methods can be applied to within-host datasets from early infection. Nonetheless, our results also suggest caution when using uniform clock and population models or short genes with limited information content.

Discovery of Aphis ruborum (Hemiptera: Aphididae) and Aphelinus varipes (Hymenoptera: Aphelinidae) on Cultivated Strawberry in Mississippi, USA.

An adventive aphid and novel host-parasitoid association from cultivated strawberry (Fragaria × ananessa Duch. cv. Chandler; Fragaria × ananessa Duch. cv. Camarosa) in Mississippi, USA are reported herein. The aphid, first detected in high tunnel cultivation, was found predominately on newly emerged, not fully developed leaflets of daughter plants in the Fall of 2016. By 2017, aphids and their associated mummies were observed on fully developed leaflets on mother plants of both cultivars. The aphid was identified as Aphis ruborum (Börner & Schilder) using morphology and DNA barcoding studies. In addition, DNA barcoding identified parasitoid adults emerging from aphid mummies as two cryptic species, Aphelinus varipes (Foerster) and Aphelinus albipodus Hayat and Fatima. Occurrence of A. ruborum in Mississippi represents a new state record and the eastern-most established record in the United States. The A. ruborum - A. varipes or A. albipodus host-parasitoid association is reported for the first time anywhere in the world.

How the Aridification of Australia Structured the Biogeography and Influenced the Diversification of a Large Lineage of Australian Cicadas.

Over the last 30 million years, Australia's landscape has undergone dramatic cooling and drying due to the establishment of the Antarctic Circumpolar Current and change in global CO$_{2}$ levels. Studies have shown that many Australian organisms went extinct during these major cooling events, while others experienced adaptive radiations and increases in diversification rates as a result of exploiting new niches in the arid zone. Despite the many studies on diversification and biogeography in Australia, few have been continent-wide and none have focused on a group of organisms adapted to feeding on plants. We studied 162 species of cicadas in the Australian Pauropsalta complex, a large generic lineage within the tribe Cicadettini. We asked whether there were changes in the diversification rate of Pauropsalta over time and if so: 1) which clades were associated with the rate change? 2) did timing of rate shifts correspond to known periods of dramatic historical climate change, 3) did increases in diversification rate along select lineages correspond to adaptive radiations with movement into the arid zone? To address these questions, we estimated a molecular phylogeny of the Pauropsalta complex using ${\sim}$5300 bp of nucleotide sequence data distributed among five loci (one mtDNA locus and four nDNA loci). We found that this large group of cicadas did not diversify at a constant rate as they spread through Australia; instead the signature of decreasing diversification rate changed roughly around the time of the expansion of the east Antarctic ice sheets ${\sim}$16 Ma and the glaciation of the northern hemisphere ${\sim}$3 Ma. Unlike other Australian taxa, the Pauropsalta complex did not explosively radiate in response to an early invasion of the arid zone. Instead multiple groups invaded the arid zone and experienced rates of diversification similar to mesic-distributed taxa. We found evidence for relictual groups, located in pre-Mesozoic habitat, that have not diversified and continue to reside on mesic hosts in isolated "habitat islands". Future work should focus on groups of similar ages with similar distribution patterns to determine whether this tempo and pattern of diversification and biogeography is consistent with evidence from other phytophagous insects.

Synthesis of phylogeny and taxonomy into a comprehensive tree of life.

Reconstructing the phylogenetic relationships that unite all lineages (the tree of life) is a grand challenge. The paucity of homologous character data across disparately related lineages currently renders direct phylogenetic inference untenable. To reconstruct a comprehensive tree of life, we therefore synthesized published phylogenies, together with taxonomic classifications for taxa never incorporated into a phylogeny. We present a draft tree containing 2.3 million tips-the Open Tree of Life. Realization of this tree required the assembly of two additional community resources: (i) a comprehensive global reference taxonomy and (ii) a database of published phylogenetic trees mapped to this taxonomy. Our open source framework facilitates community comment and contribution, enabling the tree to be continuously updated when new phylogenetic and taxonomic data become digitally available. Although data coverage and phylogenetic conflict across the Open Tree of Life illuminate gaps in both the underlying data available for phylogenetic reconstruction and the publication of trees as digital objects, the tree provides a compelling starting point for community contribution. This comprehensive tree will fuel fundamental research on the nature of biological diversity, ultimately providing up-to-date phylogenies for downstream applications in comparative biology, ecology, conservation biology, climate change, agriculture, and genomics.

The phylogenetic utility of acetyltransferase (ARD1) and glutaminyl tRNA synthetase (QtRNA) for reconstructing Cenozoic relationships as exemplified by the large Australian cicada Pauropsalta generic complex.

The Pauropsalta generic complex is a large group of cicadas (72 described spp.; >82 undescribed spp.) endemic to Australia. No previous molecular work on deep level relationships within this complex has been conducted, but a recent morphological revision and phylogenetic analysis proposed relationships among the 11 genera. We present here the first comprehensive molecular phylogeny of the complex using five loci (1 mtDNA, 4 nDNA), two of which are from nuclear genes new to cicada systematics. We compare the molecular phylogeny to the morphological phylogeny. We evaluate the phylogenetic informativeness of the new loci to traditional cicada systematics loci to generate a baseline of performance and behavior to aid in gene choice decisions in future systematic and phylogenomic studies. Our maximum likelihood and Bayesian inference phylogenies strongly support the monophyly of most of the newly described genera; however, relationships among genera differ from the morphological phylogeny. A comparison of phylogenetic informativeness among all loci revealed that COI 3rd positions dominate the informativeness profiles relative to all other loci but exhibit some among taxon nucleotide bias. After removing COI 3rd positions, COI 1st positions dominate near the terminals, while the period intron has the most phylogenetic informativeness near the root. Among the nuclear loci, ARD1 and QtRNA have lower phylogenetic informativeness than period intron and elongation factor 1 alpha intron, but the informativeness increases at you move from the tips to the root. The increase in phylogenetic informativeness deeper in the tree suggests these loci may be useful for resolving older relationships.

Phylogeographic Structure in Anastrepha ludens (Diptera: Tephritidae) Populations Inferred With mtDNA Sequencing.

Anastrepha ludens (Loew) (Diptera: Tephritidae), the Mexican fruit fly, is a major pest of citrus and mango. It has a wide distribution in Mexico and Central America, with infestations occurring in Texas, California, and Florida with origins believed to have been centered in northeastern Mexico. This research evaluates the utility of a sequence-based approach for two mitochondrial (COI and ND6) gene regions. We use these markers to examine genetic diversity, estimate population structure, and identify diagnostic information for A. ludens populations. We analyzed 543 individuals from 67 geographic collections and found one predominant haplotype occurring in the majority of specimens. We observed 68 haplotypes in all and see differences among haplotypes belonging to northern and southern collections. Mexico haplotypes differ by few bases possibly as a result of a recent bottleneck event. In contrast to the hypothesis suggesting northeastern Mexico as the origin of this species, we see that specimens from two southern collections show high genetic variability delineating three mitochondrial groups. These data suggest that Central America is the origin for A. ludens. We show that COI and ND6 are useful for phylogeographic studies of A. ludens.

Phylogeography of Anastrepha obliqua inferred with mtDNA sequencing.

Anastrepha obliqua (Macquart) (Diptera: Tephritidae), the West Indian fruit fly, is a frugivorous pest that occasionally finds its way to commercial growing areas outside its native distribution. It inhabits areas in Mexico, Central and South America, and the Caribbean with occasional infestations having occurred in the southern tier states (California, Florida, and Texas) of the United States. This fly is associated with many plant species and is a major pest of mango and plum. We examine the genetic diversity of the West Indian fruit fly based on mitochondrial COI and ND6 DNA sequences. Our analysis of 349 individuals from 54 geographic collections from Mexico, Central America, the Caribbean, and South America detected 61 haplotypes that are structured into three phylogenetic clades. The distribution of these clades among populations is associated with geography. Six populations are identified in this analysis: Mesoamerica, Central America, Caribbean, western Mexico, Andean South America, and eastern Brazil. In addition, substantial differences exist among these genetic types that warrants further taxonomic review.

Reexamination of Rhopalosiphum (Hemiptera: Aphididae) using linear discriminant analysis to determine the validity of synonymized species, with some new synonymies and distribution data.

Although 17 species of Rhopalosiphum (Hemiptera: Aphididae) are currently recognized, 85 taxonomic names have been proposed historically. Some species are morphologically similar, especially alate individuals and most synonymies were proposed in catalogues without evidence. This has led to both confusion and difficulty in making accurate species-level identifications. In an attempt to address these issues, we developed a new approach to resolve synonymies based on linear discriminant analysis (LDA) and suggest that this approach may be useful for other taxonomic groups to reassess previously proposed synonymies. We compared 34 valid and synonymized species using 49 measurements and 20 ratios from 1,030 individual aphids. LDA was repeatedly applied to subsets of the data after removing clearly separated groups found in a previous iteration. We found our characters and technique worked well to distinguish among apterae. However, it separated well only those alatae with some distinctive traits, while those apterate which were morphologically similar were not well separated using LDA. Based on our morphological investigation, we transfer R. arundinariae (Tissot, 1933) to Melanaphis supported by details of the wing veination and other morphological traits and propose Melanaphis takahashii Skvarla and Miller as a replacement name for M. arundinariae (Takahashi, 1937); we also synonymize R. momo (Shinji, 1922) with R. nymphaeae (Linnaeus, 1761). Our analyses confirmed many of the proposed synonymies, which will help to stabilize the nomenclature and species concepts within Rhopalosiphum.

Hemiptera phylogenomic resources: Tree-based orthology prediction and conserved exon identification.

High-throughput sequencing of transcriptomes and targeted genomic regions are advancing our knowledge of The Tree of Life. Building phylogenies with regions of the genome requires 1-to-1 orthologue resources of genes and noncoding loci. One organismal group that has received little attention in this area is the Hemiptera, the fifth largest insect order represented by ~103,590 named species. Here, we present a set of 3,872 Hemiptera 1-to-1 orthogroups based on tree-based orthology inference of eight Hemiptera species with publicly available genome sequences. We also estimate a set of 406 orthologous exons with similar mRNA splice sites that can be used for Sanger sequencing and develop enrichment probes for targeted genome sequencing for phylogenomic inference. We show this novel set of orthologues is informative at the protein, coding sequence and exon molecular levels and provides robust branch support in both gene tree-species tree methods and concatenated sequence phylogenies. In addition, we demonstrate the utility of these loci to resolve relationships in whiteflies, Bemisia tabaci, a large species complex with few phylogenomic resources. Last, we compare our Hemiptera phylogeny with previously published phylogenies and other orthologue databases, while providing suggestions on further improvement to this phylogenomic resource.

Towards a barnacle tree of life: integrating diverse phylogenetic efforts into a comprehensive hypothesis of thecostracan evolution.

Barnacles and their allies (Thecostraca) are a biologically diverse, monophyletic crustacean group, which includes both intensely studied taxa, such as the acorn and stalked barnacles, as well as cryptic taxa, for example, Facetotecta. Recent efforts have clarified phylogenetic relationships in many different parts of the barnacle tree, but the outcomes of these phylogenetic studies have not yet been combined into a single hypothesis for all barnacles. In the present study, we applied a new "synthesis" tree approach to estimate the first working Barnacle Tree of Life. Using this approach, we integrated phylogenetic hypotheses from 27 studies, which did not necessarily include the same taxa or used the same characters, with hierarchical taxonomic information for all recognized species. This first synthesis tree contains 2,070 barnacle species and subspecies, including 239 barnacle species with phylogenetic information and 198 undescribed or unidentified species. The tree had 442 bifurcating nodes, indicating that 79.3% of all nodes are still unresolved. We found that the acorn and stalked barnacles, the Thoracica, and the parasitic Rhizocephala have the largest amount of published phylogenetic information. About half of the thecostracan families for which phylogenetic information was available were polyphyletic. We queried publicly available geographic occurrence databases for the group, gaining a sense of geographic gaps and hotspots in our phylogenetic knowledge. Phylogenetic information is especially lacking for deep sea and Arctic taxa, but even coastal species are not fully incorporated into phylogenetic studies.

Transcriptomic analysis of changes in gene expression of immune proteins of gill tissue in response to low environmental temperature in fathead minnows (Pimephales promelas).

In the face of ongoing climate change, it is imperative to understand better the effects of temperature on immune function in freshwater teleosts. It is unclear whether previously observed changes were caused by temperature per se. We studied changes in the gill transcriptome of fathead minnows (Pimephales promelas) at low temperature to understand better the effects of temperature on immune function. De novo assembly of the transcriptome using Trinity software resulted in 73,378 assembled contigs. Annotation using the Trinotate package yielded 58,952 Blastx hits (accessions). Expression of 194 unique mRNA transcripts changed in gill tissue of fathead minnows acclimatized to 5° compared to controls at 22 °C. At 5 °C mRNAs coding for proteins involved in innate immune responses were up-regulated. Those included proteins that block early-stage viral replication and macrophage activation. Expression of mRNAs coding for pro-inflammatory molecules and mucus secretion were also enhanced. Messenger RNAs coding for proteins associated with adaptive immune responses were down-regulated at 5 °C. Those included antigen-presenting proteins and proteins involved in immunoglobin production. Messenger RNAs coding for proteins that stimulate the cell cycle were also down-regulated at 5 °C. Histological comparison revealed that gills of cold acclimated fish had fewer mucus cells but cells contained larger mucus droplets. We conclude that decreased temperature modifies the immune systems of freshwater teleosts, leading to genome-wide upregulation of innate immunity and down regulation of adaptive immunity. Such acclimation likely evolved as an adaptive strategy against seasonal changes in infectious insults.

A molecular phylogeny of the cicadas (Hemiptera: Cicadidae) with a review of tribe and subfamily classification.

A molecular phylogeny and a review of family-group classification are presented for 137 species (ca. 125 genera) of the insect family Cicadidae, the true cicadas, plus two species of hairy cicadas (Tettigarctidae) and two outgroup species from Cercopidae. Five genes, two of them mitochondrial, comprise the 4992 base-pair molecular dataset. Maximum-likelihood and Bayesian phylogenetic results are shown, including analyses to address potential base composition bias. Tettigarcta is confirmed as the sister-clade of the Cicadidae and support is found for three subfamilies identified in an earlier morphological cladistic analysis. A set of paraphyletic deep-level clades formed by African genera are together named as Tettigomyiinae n. stat. Taxonomic reassignments of genera and tribes are made where morphological examination confirms incorrect placements suggested by the molecular tree, and 11 new tribes are defined (Arenopsaltriini n. tribe, Durangonini n. tribe, Katoini n. tribe, Lacetasini n. tribe, Macrotristriini n. tribe, Malagasiini n. tribe, Nelcyndanini n. tribe, Pagiphorini n. tribe, Pictilini n. tribe, Psaltodini n. tribe, and Selymbriini n. tribe). Tribe Tacuini n. syn. is synonymized with Cryptotympanini, and Tryellina n. syn. is synonymized with an expanded Tribe Lamotialnini. Tribe Hyantiini n. syn. is synonymized with Fidicinini. Tribe Sinosenini is transferred to Cicadinae from Cicadettinae, Cicadatrini is moved to Cicadettinae from Cicadinae, and Ydiellini and Tettigomyiini are transferred to Tettigomyiinae n. stat from Cicadettinae. While the subfamily Cicadinae, historically defined by the presence of timbal covers, is weakly supported in the molecular tree, high taxonomic rank is not supported for several earlier clades based on unique morphology associated with sound production.

Phylogenetic evidence from freshwater crayfishes that cave adaptation is not an evolutionary dead-end.

Caves are perceived as isolated, extreme habitats with a uniquely specialized biota, which long ago led to the idea that caves are "evolutionary dead-ends." This implies that cave-adapted taxa may be doomed for extinction before they can diversify or transition to a more stable state. However, this hypothesis has not been explicitly tested in a phylogenetic framework with multiple independently evolved cave-dwelling groups. Here, we use the freshwater crayfish, a group with dozens of cave-dwelling species in multiple lineages, as a system to test this hypothesis. We consider historical patterns of lineage diversification and habitat transition as well as current patterns of geographic range size. We find that while cave-dwelling lineages have small relative range sizes and rarely transition back to the surface, they exhibit remarkably similar diversification patterns to those of other habitat types and appear to be able to maintain a diversity of lineages through time. This suggests that cave adaptation is not a "dead-end" for freshwater crayfish, which has positive implications for our understanding of biodiversity and conservation in cave habitats.

A synthetic phylogeny of freshwater crayfish: insights for conservation.

Phylogenetic systematics is heading for a renaissance where we shift from considering our phylogenetic estimates as a static image in a published paper and taxonomies as a hardcopy checklist to treating both the phylogenetic estimate and dynamic taxonomies as metadata for further analyses. The Open Tree of Life project (opentreeoflife.org) is developing synthesis tools for harnessing the power of phylogenetic inference and robust taxonomy to develop a synthetic tree of life. We capitalize on this approach to estimate a synthesis tree for the freshwater crayfish. The crayfish make an exceptional group to demonstrate the utility of the synthesis approach, as there recently have been a number of phylogenetic studies on the crayfishes along with a robust underlying taxonomic framework. Importantly, the crayfish have also been extensively assessed by an IUCN Red List team and therefore have accurate and up-to-date area and conservation status data available for analysis within a phylogenetic context. Here, we develop a synthesis phylogeny for the world's freshwater crayfish and examine the phylogenetic distribution of threat. We also estimate a molecular phylogeny based on all available GenBank crayfish sequences and use this tree to estimate divergence times and test for divergence rate variation. Finally, we conduct EDGE and HEDGE analyses and identify a number of species of freshwater crayfish of highest priority in conservation efforts.