aTRAM 2.0: An Improved, Flexible Locus Assembler for NGS Data

Characterization of the chemoreceptor repertoire of a highly specialized fly with comparisons to other Drosophila species

To explore the diversity of scenarios in nature, animals have evolved tools to interact with different environmental conditions. Chemoreceptors are an important interface component and among them, olfactory receptors (ORs) and gustatory receptors (GRs) can be used to find food and detect healthy resources. Drosophila is a model organism in many scientific fields, in part due to the diversity of species and niches they occupy. The contrast between generalists and specialists Drosophila species provides an important model for studying the evolution of chemoreception. Here, we compare the repertoire of chemoreceptors of different species of Drosophila with that of D. incompta, a highly specialized species whose ecology is restricted to Cestrum flowers, after reporting the preferences of D. incompta to the odor of Cestrum flowers in olfactory tests. We found evidence that the chemoreceptor repertoire in D. incompta is smaller than that presented by species in the Sophophora subgenus. Similar patterns were found in other non-Sophophora species, suggesting the presence of underlying phylogenetic trends. Nevertheless, we also found autapomorphic gene losses and detected some genes that appear to be under positive selection in D. incompta, suggesting that the specific lifestyle of these flies may have shaped the evolution of individual genes in each of these gene families.

Rapid Targeted Assembly of the Proteome Reveals Evolutionary Variation of GC Content in Avian Lice

Nucleotide base composition plays an influential role in the molecular mechanisms involved in gene function, phenotype, and amino acid composition. GC content (proportion of guanine and cytosine in DNA sequences) shows a high level of variation within and among species. Many studies measure GC content in a small number of genes, which may not be representative of genome-wide GC variation. One challenge when assembling extensive genomic data sets for these studies is the significant amount of resources (monetary and computational) associated with data processing, and many bioinformatic tools have not been optimized for resource efficiency. Using a high-performance computing (HPC) cluster, we manipulated resources provided to the targeted gene assembly program, automated target restricted assembly method (aTRAM), to determine an optimum way to run the program to maximize resource use. Using our optimum assembly approach, we assembled and measured GC content of all of the protein-coding genes of a diverse group of parasitic feather lice. Of the 499 426 genes assembled across 57 species, feather lice were GC-poor (mean GC = 42.96%) with a significant amount of variation within and between species (GC range = 19.57%-73.33%). We found a significant correlation between GC content and standard deviation per taxon for overall GC and GC3, which could indicate selection for G and C nucleotides in some species. Phylogenetic signal of GC content was detected in both GC and GC3. This research provides a large-scale investigation of GC content in parasitic lice laying the foundation for understanding the basis of variation in base composition across species.

Shifts in evolutionary lability underlie independent gains and losses of root-nodule symbiosis in a single clade of plants

Root nodule symbiosis (RNS) is a complex trait that enables plants to access atmospheric nitrogen converted into usable forms through a mutualistic relationship with soil bacteria. Pinpointing the evolutionary origins of RNS is critical for understanding its genetic basis, but building this evolutionary context is complicated by data limitations and the intermittent presence of RNS in a single clade of ca. 30,000 species of flowering plants, i.e., the nitrogen-fixing clade (NFC). We developed the most extensive de novo phylogeny for the NFC and an RNS trait database to reconstruct the evolution of RNS. Our analysis identifies evolutionary rate heterogeneity associated with a two-step process: An ancestral precursor state transitioned to a more labile state from which RNS was rapidly gained at multiple points in the NFC. We illustrate how a two-step process could explain multiple independent gains and losses of RNS, contrary to recent hypotheses suggesting one gain and numerous losses, and suggest a broader phylogenetic and genetic scope may be required for genome-phenome mapping.

GeneMiner: A tool for extracting phylogenetic markers from next-generation sequencing data

The advancement of next-generation sequencing (NGS) technologies has been revolutionary for the field of evolutionary biology. This technology has led to an abundance of available genomes and transcriptomes for researchers to mine. Specifically, researchers can mine for various types of molecular markers that are vital for phylogenetic, evolutionary and ecological studies. Numerous tools have been developed to extract these molecular markers from NGS data. However, due to an insufficient number of well-annotated reference genomes for non-model organisms, it remains challenging to obtain these markers accurately and efficiently. Here, we present GeneMiner, an improved and expanded version of our previous tool, Easy353. GeneMiner combines the reference-guided de Bruijn graph assembly with seed self-discovery and greedy extension. Additionally, it includes a verification step using a parameter-bootstrap method to reduce the pitfalls associated with using a relatively distant reference. Our results, using both experimental and simulation data, showed GeneMiner can accurately acquire phylogenetic molecular markers for plants using transcriptomic, genomic and other NGS data. GeneMiner is designed to be user-friendly, fast and memory-efficient. Further, it is compatible with Linux, Windows and macOS. All source codes are publicly available on GitHub (https://github.com/sculab/GeneMiner) and Gitee (https://gitee.com/sculab/GeneMiner) for easy accessibility and transparency.

Patchwork: Alignment-Based Retrieval and Concatenation of Phylogenetic Markers from Genomic Data

Low-coverage whole-genome sequencing (also known as "genome skimming") is becoming an increasingly affordable approach to large-scale phylogenetic analyses. While already routinely used to recover organellar genomes, genome skimming is rather rarely utilized for recovering single-copy nuclear markers. One reason might be that only few tools exist to work with this data type within a phylogenomic context, especially to deal with fragmented genome assemblies. We here present a new software tool called Patchwork for mining phylogenetic markers from highly fragmented short-read assemblies as well as directly from sequence reads. Patchwork is an alignment-based tool that utilizes the sequence aligner DIAMOND and is written in the programming language Julia. Homologous regions are obtained via a sequence similarity search, followed by a "hit stitching" phase, in which adjacent or overlapping regions are merged into a single unit. The novel sliding window algorithm trims away any noncoding regions from the resulting sequence. We demonstrate the utility of Patchwork by recovering near-universal single-copy orthologs within a benchmarking study, and we additionally assess the performance of Patchwork in comparison with other programs. We find that Patchwork allows for accurate retrieval of (putatively) single-copy genes from genome skimming data sets at different sequencing depths with high computational speed, outperforming existing software targeting similar tasks. Patchwork is released under the GNU General Public License version 3. Installation instructions, additional documentation, and the source code itself are all available via GitHub at https://github.com/fethalen/Patchwork.

Novel insights into symbiont population structure: Globe-trotting avian feather mites contradict the specialist-generalist variation hypothesis

Researchers often examine symbiont host specificity as a species-level pattern, but it can also be key to understanding processes occurring at the population level, which are not as well understood. The specialist-generalist variation hypothesis (SGVH) attempts to explain how host specificity influences population-level processes, stating that single-host symbionts (specialists) exhibit stronger population genetic structure than multi-host symbionts (generalists) because of fewer opportunities for dispersal and more restricted gene flow between populations. However, this hypothesis has not been tested in systems with highly mobile hosts, in which population connectivity may vary temporally and spatially. To address this gap, we tested the SGVH on proctophyllodid feather mites found on migratory warblers (family Parulidae) with contrasting host specificities, Amerodectes protonotaria (a host specialist of Protonotaria citrea) and A. ischyros (a host generalist of 17 parulid species). We used a pooled-sequencing approach and a novel workflow to analyse genetic variants obtained from whole genome data. Both mite species exhibited fairly weak population structure overall, and contrary to predictions of the SGVH, the generalist was more strongly structured than the specialist. These results may suggest that specialists disperse more freely among conspecifics, whereas generalists sort according to geography. Furthermore, our results may reflect an unexpected period for mite transmission - during the nonbreeding season of migratory hosts - as mite population structure more closely reflects the distributions of hosts during the nonbreeding season. Our findings alter our current understanding of feather mite biology and highlight the potential for studies to explore factors driving symbiont diversification at multiple evolutionary scales.

Phylogenomics of transcriptionally active AP2/ERF and bHLH transcription factors and study of their promoter regions in Nothapodytes nimmoniana (J.Graham) Mabb

Nothapodytes nimmoniana is a medicinally important plant producing anticancer monoterpene indole alkaloid (MIA), camptothecin (CPT). The CPT is synthesised through the strictosidine intermediate following the MIA pathway; however, transcriptional regulation of CPT pathway is still elusive in N. nimmoniana. Biosynthesis of MIA is regulated by various transcription factors (TFs) belonging to AP2/ERF, bHLH, MYB, and WRKY families. The present study identified transcriptionally active full-length 105 AP2/ERF and 68 bHLH family TFs from the N. nimmoniana. AP2/ERF TFs were divided into three subfamilies along with a soloist, while bHLH TFs were divided into 10 subfamilies according to their phylogenetic similarities. Three group IXa ERFs, Nn-ERF22, Nn-ERF29, and Nn-ERF41, one subfamily IVa TF Nn-bHLH7, and three subfamilies IIIe Nn-bHLH33, Nn-bHLH51, and Nn-bHLH52 clustered with the TFs regulating alkaloid biosynthesis in Catharanthus roseus, tomato, tobacco, and Artemisia annua. Expression of these TFs in N. nimmoniana was higher in roots, which is a primary CPT accumulating tissue. Moreover, genome skimming approach was used to reconstruct the promoter regions of candidate ERF genes to identify the cis-regulatory elements. The presence of G-boxes and other jasmonic acid-responsive elements in the promoter suggests the regulation of ERFs by bHLHs. The present study effectively generated and used genomics resource for characterisation of regulatory TFs from non-model medicinal plant.

Phylogenomic analysis provides diagnostic tools for the identification of Anastrepha fraterculus (Diptera: Tephritidae) species complex

Insect pests cause tremendous impact to agriculture worldwide. Species identification is crucial for implementing appropriate measures of pest control but can be challenging in closely related species. True fruit flies of the genus Anastrepha Schiner (Diptera: Tephritidae) include some of the most serious agricultural pests in the Americas, with the Anastrepha fraterculus (Wiedemann) complex being one of the most important due to its extreme polyphagy and wide distribution across most of the New World tropics and subtropics. The eight morphotypes described for this complex as well as other closely related species are classified in the fraterculus species group, whose evolutionary relationships are unresolved due to incomplete lineage sorting and introgression. We performed multifaceted phylogenomic approaches using thousands of genes to unravel the evolutionary relationships within the A. fraterculus complex to provide a baseline for molecular diagnosis of these pests. We used a methodology that accommodates variable sources of data (transcriptome, genome, and whole-genome shotgun sequencing) and developed a tool to align and filter orthologs, generating reliable datasets for phylogenetic studies. We inferred 3031 gene trees that displayed high levels of discordance. Nevertheless, the topologies of the inferred coalescent species trees were consistent across methods and datasets, except for one lineage in the A. fraterculus complex. Furthermore, network analysis indicated introgression across lineages in the fraterculus group. We present a robust phylogeny of the group that provides insights into the intricate patterns of evolution of the A. fraterculus complex supporting the hypothesis that this complex is an assemblage of closely related cryptic lineages that have evolved under interspecific gene flow. Despite this complex evolutionary scenario, our subsampling analysis revealed that a set of as few as 80 loci has a similar phylogenetic resolution as the genome-scale dataset, offering a foundation to develop more efficient diagnostic tools in this species group.

Dispersal-Limited Symbionts Exhibit Unexpectedly Wide Variation in Host Specificity

A fundamental aspect of symbiotic relationships is host specificity, ranging from extreme specialists associated with only a single host species to generalists associated with many different species. Although symbionts with limited dispersal capabilities are expected to be host specialists, some are able to associate with multiple hosts. Understanding the micro- and macro-evolutionary causes of variations in host specificity is often hindered by sampling biases and the limited power of traditional evolutionary markers. Here, we studied feather mites to address the barriers associated with estimates of host specificity for dispersal-limited symbionts. We sampled feather mites (Proctophyllodidae) from a nearly comprehensive set of North American breeding warblers (Parulidae) to study mite phylogenetic relationships and host-symbiont codiversification. We used pooled-sequencing (Pool-Seq) and short-read Illumina technology to interpret results derived from a traditional barcoding gene (cytochrome c oxidase subunit 1) versus 11 protein-coding mitochondrial genes using concatenated and multispecies coalescent approaches. Despite the statistically significant congruence between mite and host phylogenies, mite-host specificity varies widely, and host switching is common regardless of the genetic marker resolution (i.e., barcode vs. multilocus). However, the multilocus approach was more effective than the single barcode in detecting the presence of a heterogeneous Pool-Seq sample. These results suggest that presumed symbiont dispersal capabilities are not always strong indicators of host specificity or of historical host-symbiont coevolutionary events. A comprehensive sampling at fine phylogenetic scales may help to better elucidate the microevolutionary filters that impact macroevolutionary processes regulating symbioses, particularly for dispersal-limited symbionts. [Codiversification; cophylogenetics; feather mites; host switching; pooled sequencing; species delineation; symbiosis, warblers.].

Identifying Climatic Drivers of Hybridization with a New Ancestral Niche Reconstruction Method

Applications of molecular phylogenetic approaches have uncovered evidence of hybridization across numerous clades of life, yet the environmental factors responsible for driving opportunities for hybridization remain obscure. Verbal models implicating geographic range shifts that brought species together during the Pleistocene have often been invoked, but quantitative tests using paleoclimatic data are needed to validate these models. Here, we produce a phylogeny for Heuchereae, a clade of 15 genera and 83 species in Saxifragaceae, with complete sampling of recognized species, using 277 nuclear loci and nearly complete chloroplast genomes. We then employ an improved framework with a coalescent simulation approach to test and confirm previous hybridization hypotheses and identify one new intergeneric hybridization event. Focusing on the North American distribution of Heuchereae, we introduce and implement a newly developed approach to reconstruct potential past distributions for ancestral lineages across all species in the clade and across a paleoclimatic record extending from the late Pliocene. Time calibration based on both nuclear and chloroplast trees recovers a mid- to late-Pleistocene date for most inferred hybridization events, a timeframe concomitant with repeated geographic range restriction into overlapping refugia. Our results indicate an important role for past episodes of climate change, and the contrasting responses of species with differing ecological strategies, in generating novel patterns of range contact among plant communities and therefore new opportunities for hybridization. The new ancestral niche method flexibly models the shape of niche while incorporating diverse sources of uncertainty and will be an important addition to the current comparative methods toolkit. [Ancestral niche reconstruction; hybridization; paleoclimate; pleistocene.].

MTG-Link: leveraging barcode information from linked-reads to assemble specific loci

BACKGROUND

Local assembly with short and long reads has proven to be very useful in many applications: reconstruction of the sequence of a locus of interest, gap-filling in draft assemblies, as well as alternative allele reconstruction of large Structural Variants. Whereas linked-read technologies have a great potential to assemble specific loci as they provide long-range information while maintaining the power and accuracy of short-read sequencing, there is a lack of local assembly tools for linked-read data.

RESULTS

We present MTG-Link, a novel local assembly tool dedicated to linked-reads. The originality of the method lies in its read subsampling step which takes advantage of the barcode information contained in linked-reads mapped in flanking regions. We validated our approach on several datasets from different linked-read technologies. We show that MTG-Link is able to assemble successfully large sequences, up to dozens of Kb. We also demonstrate that the read subsampling step of MTG-Link considerably improves the local assembly of specific loci compared to other existing short-read local assembly tools. Furthermore, MTG-Link was able to fully characterize large insertion variants and deletion breakpoints in a human genome and to reconstruct dark regions in clinically-relevant human genes. It also improved the contiguity of a 1.3 Mb locus of biological interest in several individual genomes of the mimetic butterfly Heliconius numata.

CONCLUSIONS

MTG-Link is an efficient local assembly tool designed for different linked-read sequencing technologies. MTG-Link source code is available at https://github.com/anne-gcd/MTG-Link and as a Bioconda package.

Phylogenomics reveals widespread hybridization and polyploidization in Henckelia (Gesneriaceae)

BACKGROUND AND AIMS

Hybridization has long been recognized as an important process for plant evolution and is often accompanied by polyploidization, another prominent force in generating biodiversity. Despite its pivotal importance in evolution, the actual prevalence and distribution of hybridization across the tree of life remain unclear.

METHODS

We used whole-genome shotgun (WGS) sequencing and cytological data to investigate the evolutionary history of Henckelia, a large genus in the family Gesneriaceae with a high frequency of suspected hybridization and polyploidization events. We generated WGS sequencing data at about 10× coverage for 26 Chinese Henckelia species plus one Sri Lankan species. To untangle the hybridization history, we separately extracted whole plastomes and thousands of single-copy nuclear genes from the sequencing data, and reconstructed phylogenies based on both nuclear and plastid data. We also explored sources of both genealogical and cytonuclear conflicts and identified signals of hybridization and introgression within our phylogenomic dataset using several statistical methods. Additionally, to test the polyploidization history, we evaluated chromosome counts for 45 populations of the 27 Henckelia species studied.

KEY RESULTS

We obtained well-supported phylogenetic relationships using both concatenation- and coalescent-based methods. However, the nuclear phylogenies were highly inconsistent with the plastid phylogeny, and we observed intensive discordance among nuclear gene trees. Further analyses suggested that both incomplete lineage sorting and gene flow contributed to the observed cytonuclear and genealogical discordance. Our analyses of introgression and phylogenetic networks revealed a complex history of hybridization within the genus Henckelia. In addition, based on chromosome counts for 27 Henckelia species, we found independent polyploidization events occurred within Henckelia after different hybridization events.

CONCLUSIONS

Our findings demonstrated that hybridization and polyploidization are common in Henckelia. Furthermore, our results revealed that H. oblongifolia is not a member of the redefined Henckelia and they suggested several other taxonomic treatments in this genus.

Phylogenetic and Evolutionary Comparison of Mitogenomes Reveal Adaptive Radiation of Lampriform Fishes

Lampriform fishes (Lampriformes), which primarily inhabit deep-sea environments, are large marine fishes varying from the whole-body endothermic opah to the world's longest bony fish-giant oarfish, with species morphologies varying from long and thin to deep and compressed, making them an ideal model for studying the adaptive radiation of teleost fishes. Moreover, this group is important from a phylogenetic perspective owing to their ancient origins among teleosts. However, knowledge about the group is limited, which is, at least partially, due to the dearth of recorded molecular data. This study is the first to analyze the mitochondrial genomes of three lampriform species (Lampris incognitus, Trachipterus ishikawae, and Regalecus russelii) and infer a time-calibrated phylogeny, including 68 species among 29 orders. Our phylomitogenomic analyses support the classification of Lampriformes as monophyletic and sister to Acanthopterygii; hence, addressing the longstanding controversy regarding the phylogenetic status of Lampriformes among teleosts. Comparative mitogenomic analyses indicate that tRNA losses existed in at least five Lampriformes species, which may reveal the mitogenomic structure variation associated with adaptive radiation. However, codon usage in Lampriformes did not change significantly, and it is hypothesized that the nucleus transported the corresponding tRNA, which led to function substitutions. The positive selection analysis revealed that atp8 and cox3 were positively selected in opah, which might have co-evolved with the endothermic trait. This study provides important insights into the systematic taxonomy and adaptive evolution studies of Lampriformes species.

Late Cenozoic History And The Role Of Beringia In Assembling A Holarctic Cestode Species Complex

The dynamic climate history that drove sea level fluctuation during past glacial periods mediated the movement of organisms between Asia and North America via the Bering Land Bridge. Investigations of the biogeographic histories of small mammals and their parasites demonstrate facets of a complex history of episodic geographic colonization and refugial isolation that structured diversity across the Holarctic. We use a large multi-locus nuclear DNA sequence dataset to robustly resolve relationships within the cestode genus Arostrilepis (Cyclophyllidea: Hymenolepididae), a widespread parasite of predominantly arvicoline rodents (voles, lemmings). Using this phylogeny, we confirm that several Asian Arostrilepis lineages colonized North America during up to four distinct glacial periods in association with different rodent hosts, consistent with taxon-pulse dynamics. A previously inferred westward dispersal across the land bridge is rejected. We also refine interpretations of past host colonization, providing evidence for several distinct episodes of expanding host range, which probably contributed to diversification by Arostrilepis. Finally, Arostrilepis is shown to be paraphyletic with respect to Hymenandrya thomomyis, a parasite of pocket gophers, confirming that ancient Arostrilepis species colonized new host lineages upon arriving in North America.

Supergroup F Wolbachia with extremely reduced genome: transition to obligate insect symbionts

BACKGROUND

Wolbachia belong to highly abundant bacteria which are frequently found in invertebrate microbiomes and manifest by a broad spectrum of lifestyles from parasitism to mutualism. Wolbachia supergroup F is a particularly interesting clade as it gave rise to symbionts of both arthropods and nematodes, and some of its members are obligate mutualists. Investigations on evolutionary transitions among the different symbiotic stages have been hampered by a lack of the known diversity and genomic data for the supergroup F members.

RESULTS

Based on amplicon screening, short- and long-read WGS approaches, and laser confocal microscopy, we characterize five new supergroup F Wolbachia strains from four chewing lice species. These strains reached different evolutionary stages and represent two remarkably different types of symbiont genomes. Three of the genomes resemble other known members of Wolbachia F supergroup, while the other two show typical signs of ongoing gene inactivation and removal (genome size, coding density, low number of pseudogenes). Particularly, wMeur1, a symbiont fixed in microbiomes of Menacanthus eurysternus across four continents, possesses a highly reduced genome of 733,850 bp. The horizontally acquired capacity for pantothenate synthesis and localization in specialized bacteriocytes suggest its obligate nutritional role.

CONCLUSIONS

The genome of wMeur1 strain, from the M. eurysternus microbiome, represents the smallest currently known Wolbachia genome and the first example of Wolbachia which has completed genomic streamlining as known from the typical obligate symbionts. This points out that despite the large amount and great diversity of the known Wolbachia strains, evolutionary potential of these bacteria still remains underexplored. The diversity of the four chewing lice microbiomes indicates that this vast parasitic group may provide suitable models for further investigations. Video Abstract.

Phylogenomics and the flowering plant tree of life

The advances accelerated by next-generation sequencing and long-read sequencing technologies continue to provide an impetus for plant phylogenetic study. In the past decade, a large number of phylogenetic studies adopting hundreds to thousands of genes across a wealth of clades have emerged and ushered plant phylogenetics and evolution into a new era. In the meantime, a roadmap for researchers when making decisions across different approaches for their phylogenomic research design is imminent. This review focuses on the utility of genomic data (from organelle genomes, to both reduced representation sequencing and whole-genome sequencing) in phylogenetic and evolutionary investigations, describes the baseline methodology of experimental and analytical procedures, and summarizes recent progress in flowering plant phylogenomics at the ordinal, familial, tribal, and lower levels. We also discuss the challenges, such as the adverse impact on orthology inference and phylogenetic reconstruction raised from systematic errors, and underlying biological factors, such as whole-genome duplication, hybridization/introgression, and incomplete lineage sorting, together suggesting that a bifurcating tree may not be the best model for the tree of life. Finally, we discuss promising avenues for future plant phylogenomic studies.

Draft genome sequencing data of a feather mite, Amerodectes protonotaria Hernandes 2018 (Acariformes: Proctophyllodidae)

Feather mites are ubiquitous, permanent, obligate ectosymbionts of avian hosts and are a valuable natural system for studying host-symbiont evolutionary and ecological dynamics at multiple levels of biological organization. However, a lack of a sequenced genome impedes molecular studies using this system. Therefore, we present the first draft genome of a symbiotic feather mite, Amerodectes protonotaria Hernandes 2018. The genome sequence data presented here were derived from an individual female mite that was collected in the field from Protonotaria citrea, its only known host species. Short read sequence data were obtained using an Illumina NovaSeq 6000 platform. From these data, we assembled a 59,665,063 bp draft genome consisting of 2,399 contigs. Raw short reads and the assembled genome sequence are available at the National Center for Biotechnology Information (NCBI)'s Sequence Read Archive (SRA) under BioProject PRJNA884722. The data presented here are beneficial for future research on the biology and evolution of closely related mites and the genomics of host-symbiont interactions.

The mitochondrial genome of the red icefish (Channichthys rugosus) casts doubt on its species status

Antarctic notothenioid fishes are recognised as one of the rare examples of adaptive radiation in the marine system. Withstanding the freezing temperatures of Antarctic waters, these fishes have diversified into over 100 species within no more than 10–20 million years. However, the exact species richness of the radiation remains contested. In the genus Channichthys, between one and nine species are recognised by different authors. To resolve the number of Channichthys species, genetic information would be highly valuable; however, so far, only sequences of a single species, C. rhinoceratus, are available. Here, we present the nearly complete sequence of the mitochondrial genome of C. rugosus, obtained from a formalin-fixed museum specimen sampled in 1974. This sequence differs from the mitochondrial genome of C. rhinoceratus in no more than 27 positions, suggesting that the two species may be synonymous.

Genomic Approaches to Uncovering the Coevolutionary History of Parasitic Lice

Simple Summary

New sequencing technologies have now made it possible to sequence entire genomes for a diversity of life on earth. Parasites comprise nearly half of all species. Lice are one important group of parasites of birds and mammals, including humans. Genome sequencing approaches have been applied to this group of parasites to uncover patterns of diversification. These patterns can be compared to the patterns of diversification in their hosts. Key findings from these studies have revealed that parasitic lice likely originated on birds and then switched to mammals multiple times. Within groups of birds and mammals, the evolutionary trees of lice match those for mammal hosts more than those for birds. Genomic approaches have also revealed that individual birds and mammals harbor distinct populations of lice. Thus, these new techniques allow for the study of patterns of diversification at a wide variety of scales.

Abstract

Next-generation sequencing technologies are revolutionizing the fields of genomics, phylogenetics, and population genetics. These new genomic approaches have been extensively applied to a major group of parasites, the lice (Insecta: Phthiraptera) of birds and mammals. Two louse genomes have been assembled and annotated to date, and these have opened up new resources for the study of louse biology. Whole genome sequencing has been used to assemble large phylogenomic datasets for lice, incorporating sequences of thousands of genes. These datasets have provided highly supported trees at all taxonomic levels, ranging from relationships among the major groups of lice to those among closely related species. Such approaches have also been applied at the population scale in lice, revealing patterns of population subdivision and inbreeding. Finally, whole genome sequence datasets can also be used for additional study beyond that of the louse nuclear genome, such as in the study of mitochondrial genome fragmentation or endosymbiont function.

Independent evolution of highly variable, fragmented mitogenomes of parasitic lice

The mitochondrial genomes (mitogenomes) of bilaterian animals are highly conserved structures that usually consist of a single circular chromosome. However, several species of parasitic lice (Insecta: Phthiraptera) possess fragmented mitogenomes, where the mitochondrial genes are present on separate, circular chromosomes. Nevertheless, the extent, causes, and consequences of this structural variation remain poorly understood. Here, we combined new and existing data to better understand the evolution of mitogenome fragmentation in major groups of parasitic lice. We found strong evidence that fragmented mitogenomes evolved many times within parasitic lice and that the level of fragmentation is highly variable, including examples of heteroplasmic arrangements. We also found a significant association between mitochondrial fragmentation and signatures of relaxed selection. Mitochondrial fragmentation was also associated with changes to a lower AT%, possibly due to differences in mutation biases. Together, our results provide a significant advance in understanding the process of mitogenome fragmentation and provide an important perspective on mitochondrial evolution in eukaryotes.

Phylogenomics reveals the origin of mammal lice out of Afrotheria

Mammals host a wide diversity of parasites. Lice, comprising more than 5,000 species, are one group of ectoparasites whose major lineages have a somewhat patchwork distribution across the major groups of mammals. Here we explored patterns in the diversification of mammalian lice by reconstructing a higher-level phylogeny of these lice, leveraging whole genome sequence reads to assemble single-copy orthologue genes across the genome. The evolutionary tree of lice indicated that three of the major lineages of placental mammal lice had a single common ancestor. Comparisons of this parasite phylogeny with that for their mammalian hosts indicated that the common ancestor of elephants, elephant shrews and hyraxes (that is, Afrotheria) was the ancestral host of this group of lice. Other groups of placental mammals obtained their lice via host-switching out of these Afrotherian ancestors. In addition, reconstructions of the ancestral host group (bird versus mammal) for all parasitic lice supported an avian ancestral host, indicating that the ancestor of Afrotheria acquired these parasites via host-switching from an ancient avian host. These results shed new light on the long-standing question of why the major groups of parasitic lice are not uniformly distributed across mammals and reveal the origins of mammalian lice.

Receptor editing constrains development of phosphatidyl choline-specific B cells in VH12-transgenic mice

B1 B cells reactive to phosphatidyl choline (PtC) exhibit restricted immunoglobulin heavy chain (HC) and light chain (LC) combinations, exemplified by V_H12/Vκ4/5H. Two checkpoints are thought to focus PtC⁺ B cell maturation in V_H12-transgenic mice (VH12 mice): V-J rearrangements encoding a "permissive" LC capable of V_H12 HC pairing are selected first, followed by positive selection based on PtC binding, often requiring LC receptor editing to salvage PtC^- B cells and acquire PtC reactivity. However, evidence obtained from breeding VH12 mice to editing-defective dnRAG1 mice and analyzing LC sequences from PtC⁺ and PtC^- B cell subsets instead suggests that receptor editing functions after initial positive selection to remove PtC⁺ B cells in VH12 mice. This offers a mechanism to constrain natural, polyreactive B cells to limit their frequency. Sequencing also reveals occasional in-frame hybrid LC genes, reminiscent of type 2 gene replacement, that, testing suggests, arise via a recombination-activating gene (RAG)-independent mechanism.

Local Adaptation and the Evolution of Genome Architecture in Threespine Stickleback

Theory predicts that local adaptation should favor the evolution of a concentrated genetic architecture, where the alleles driving adaptive divergence are tightly clustered on chromosomes. Adaptation to marine versus freshwater environments in threespine stickleback has resulted in an architecture that seems consistent with this prediction: divergence among populations is mainly driven by a few genomic regions harboring multiple quantitative trait loci for environmentally adapted traits, as well as candidate genes with well-established phenotypic effects. One theory for the evolution of these "genomic islands" is that rearrangements remodel the genome to bring causal loci into tight proximity, but this has not been studied explicitly. We tested this theory using synteny analysis to identify micro- and macro-rearrangements in the stickleback genome and assess their potential involvement in the evolution of genomic islands. To identify rearrangements, we conducted a de novo assembly of the closely related tubesnout (Aulorhyncus flavidus) genome and compared this to the genomes of threespine stickleback and two other closely related species. We found that small rearrangements, within-chromosome duplications, and lineage-specific genes (LSGs) were enriched around genomic islands, and that all three chromosomes harboring large genomic islands have experienced macro-rearrangements. We also found that duplicates and micro-rearrangements are 9.9× and 2.9× more likely to involve genes differentially expressed between marine and freshwater genotypes. While not conclusive, these results are consistent with the explanation that strong divergent selection on candidate genes drove the recruitment of rearrangements to yield clusters of locally adaptive loci.

Protease and gag diversity and drug resistance mutations among treatment-naive Mexican people living with HIV

Introduction

In Mexico, HIV genotyping is performed in people living with HIV (PLWH) failing their first-line antiretroviral (ARV) regimen; it is not routinely done for all treatment-naive PLWH before ARV initiation. The first nationally representative survey published in 2016 reported that the prevalence of pretreatment drug mutations in treatment-naive Mexican PLWH was 15.5% to any antiretroviral drug and 10.6% to non-nucleoside reverse transcriptase inhibitors (NNRTIs) using conventional Sanger sequencing. Most reports in Mexico focus on HIV pol gene and nucleoside and non-nucleoside reverse transcriptase inhibitor (NRTI and NNRTI) drug resistance mutations (DRMs) prevalence, using Sanger sequencing, next-generation sequencing (NGS) or both. To our knowledge, NGS has not be used to detect pretreatment drug resistance mutations (DRMs) in the HIV protease (PR) gene and its substrate the Gag polyprotein.

Methods

Treatment-naive adult Mexican PLWH were recruited between 2016 and 2019. HIV Gag and protease sequences were obtained by NGS and DRMs were identified using the WHO surveillance drug resistance mutation (SDRM) list.

Results

One hundred PLWH attending a public national reference hospital were included. The median age was 28 years-old, and most were male. The median HIV viral load was 4.99 [4.39–5.40] log copies/mL and median CD4 cell count was 150 [68.0–355.78] cells/mm³. As expected, most sequences clustered with HIV-1 subtype B (97.9%). Major PI resistance mutations were detected: 8 (8.3%) of 96 patients at a detection threshold of 1% and 3 (3.1%) at a detection threshold of 20%. A total of 1184 mutations in Gag were detected, of which 51 have been associated with resistance to PI, most of them were detected at a threshold of 20%. Follow-up clinical data was available for 79 PLWH at 6 months post-ART initiation, seven PLWH failed their first ART regimen; however no major PI mutations were identified in these individuals at baseline.

Conclusions

The frequency of DRM in the HIV protease was 7.3% at a detection threshold of 1% and 3.1% at a detection threshold of 20%. NGS-based HIV drug resistance genotyping provide improved detection of DRMs. Viral load was used to monitor ARV response and treatment failure was 8.9%.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12879-022-07446-8.

Diversification of the phytophagous lineages of true bugs (Insecta: Hemiptera: Heteroptera) shortly after that of the flowering plants

More than 95% of phytophagous true bug (Hemiptera: Heteroptera) species belong to four superfamilies: Miroidea (Cimicomorpha), Pentatomoidea, Coreoidea, and Lygaeoidea (all Pentatomomorpha). These iconic groups of highly diverse, overwhelmingly phytophagous insects include several economically prominent agricultural and silvicultural pest species, though their evolutionary history has not yet been well resolved. In particular, superfamily- and family-level phylogenetic relationships of these four lineages have remained controversial, and the divergence times of some crucial nodes for phytophagous true bugs have hitherto been little known, which hampers a better understanding of the evolutionary processes and patterns of phytophagous insects. In the present study, we used 150 species and concatenated nuclear and mitochondrial protein-coding genes and rRNA genes to infer the phylogenetic relationships within the Terheteroptera (Cimicomorpha + Pentatomomorpha) and estimated their divergence times. Our results support the monophyly of Cimicomorpha, Pentatomomorpha, Miroidea, Pentatomoidea, Pyrrhocoroidea, Coreoidea, and Lygaeoidea. The phylogenetic relationships across phytophagous lineages are largely congruent at deep nodes across the analyses based on different datasets and tree-reconstructing methods with just a few exceptions. Estimated divergence times and ancestral state reconstructions for feeding habit indicate that phytophagous true bugs explosively radiated in the Early Cretaceous-shortly after the angiosperm radiation-with the subsequent diversification of the most speciose clades (Mirinae, Pentatomidae, Coreinae, and Rhyparochromidae) in the Late Cretaceous.

SLAG: A Program for Seeded Local Assembly of Genes in Complex Genomes

Although finished genomes have become more common, there is still a need for assemblies of individual genes or chromosomal regions when only unassembled reads are available. slag (Seeded Local Assembly of Genes) fulfils this need by performing iterative local assembly based on cycles of matching‐read retrieval with blast and assembly with cap3, phrap, spades, canu or unicycler. The target sequence can be nucleotide or protein. Read fragmentation allows slag to use phrap or cap3 to assemble long reads at lower coverage (e.g., 5×) than is possible with canu or unicycler. In simple, nonrepetitive genomes, a slag assembly can cover a whole chromosome, but in complex genomes the growth of target‐matching contigs is limited as additional reads are consumed by consensus contigs consisting of repetitive elements. Apart from genomic complexity, contig length and correctness depend on read length and accuracy. With pyrosequencing or Illumina reads, slag‐assembled contigs are accurate enough to allow design of PCR primers, while contigs assembled from Oxford Nanopore or pre‐HiFi Pacific Biosciences long reads are generally only accurate enough to design baiting sequences for further targeted sequencing. In an application with real reads, slag successfully extended sequences for four wheat genes, which were verified by cloning and Sanger sequencing of overlapping amplicons. slag is a robust alternative to atram2 for local assemblies, especially for read sets with less than 20× coverage. slag is freely available at https://github.com/cfcrane/SLAG.

The first true millipede-1306 legs long

The name “millipede” translates to a thousand feet (from mille “thousand” and pes “foot”). However, no millipede has ever been described with more than 750 legs. We discovered a new record-setting species of millipede with 1,306 legs, Eumillipes persephone, from Western Australia. This diminutive animal (0.95 mm wide, 95.7 mm long) has 330 segments, a cone-shaped head with enormous antennae, and a beak for feeding. A distant relative of the previous record holder, Illacme plenipes from California, it belongs to a different order, the Polyzoniida. Discovered 60 m below ground in a drill hole created for mineral exploration, E. persephone possesses troglomorphic features; it lacks eyes and pigmentation, and it has a greatly elongated body—features that stand in stark contrast to its closest surface-dwelling relatives in Australia and all other members of its order. Using phylogenomics, we found that super-elongation (> 180 segments) evolved repeatedly in the millipede class Diplopoda. The striking morphological similarity between E. persephone and I. plenipes is a result of convergent evolution, probably for locomotion in similar soil habitats. Discovered in the resource-rich Goldfields-Esperance region and threatened by encroaching surface mining, documentation of this species and conservation of its habitat are of critical importance.

The interplay between host biogeography and phylogeny in structuring diversification of the feather louse genus Penenirmus

Parasite diversification is influenced by many of the same factors that affect speciation of free-living organisms, such as biogeographic barriers. However, the ecology and evolution of the host lineage also has a major impact on parasite speciation. Here we explore the interplay between biogeography and host-association on the pattern of diversification in a group of ectoparasitic lice (Insecta: Phthiraptera: Penenirmus) that feeds on the feathers of woodpeckers, barbets, and honeyguides (Piciformes) and some songbirds (Passeriformes). We use whole genome sequencing of 41 ingroup and 12 outgroup samples to develop a phylogenomic dataset of DNA sequences from a reference set of 2395 single copy ortholog genes, for a total of nearly four million aligned base positions. The phylogenetic trees resulting from both concatenated and gene-tree/species-tree coalescent analyses were nearly identical and highly supported. These trees recovered the genus Penenirmus as monophyletic and identified several major clades, which tended to be associated with one major host group. However, cophylogenetic analysis revealed that host-switching was a prominent process in the diversification of this group. This host-switching generally occurred within single major biogeographic regions. We did, however, find one case in which it appears that a rare dispersal event by a woodpecker lineage from North America to Africa allowed its associated louse to colonize a woodpecker in Africa, even though the woodpecker lineage from North America never became established there.

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.

The Effects of Herbarium Specimen Characteristics on Short-Read NGS Sequencing Success in Nearly 8000 Specimens: Old, Degraded Samples Have Lower DNA Yields but Consistent Sequencing Success

Phylogenetic datasets are now commonly generated using short-read sequencing technologies unhampered by degraded DNA, such as that often extracted from herbarium specimens. The compatibility of these methods with herbarium specimens has precipitated an increase in broad sampling of herbarium specimens for inclusion in phylogenetic studies. Understanding which sample characteristics are predictive of sequencing success can guide researchers in the selection of tissues and specimens most likely to yield good results. Multiple recent studies have considered the relationship between sample characteristics and DNA yield and sequence capture success. Here we report an analysis of the relationship between sample characteristics and sequencing success for nearly 8,000 herbarium specimens. This study, the largest of its kind, is also the first to include a measure of specimen quality ("greenness") as a predictor of DNA sequencing success. We found that taxonomic group and source herbarium are strong predictors of both DNA yield and sequencing success and that the most important specimen characteristics for predicting success differ for DNA yield and sequencing: greenness was the strongest predictor of DNA yield, and age was the strongest predictor of proportion-on-target reads recovered. Surprisingly, the relationship between age and proportion-on-target reads is the inverse of expectations; older specimens performed slightly better in our capture-based protocols. We also found that DNA yield itself is not a strong predictor of sequencing success. Most literature on DNA sequencing from herbarium specimens considers specimen selection for optimal DNA extraction success, which we find to be an inappropriate metric for predicting success using next-generation sequencing technologies.

Continued Adaptation of C4 Photosynthesis After an Initial Burst of Changes in the Andropogoneae Grasses

C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{4}$\end{document} photosynthesis is a complex trait that sustains fast growth and high productivity in tropical and subtropical conditions and evolved repeatedly in flowering plants. One of the major C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{4}$\end{document} lineages is Andropogoneae, a group of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$\sim $\end{document}1200 grass species that includes some of the world’s most important crops and species dominating tropical and some temperate grasslands. Previous efforts to understand C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{4}$\end{document} evolution in the group have compared a few model C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{4}$\end{document} plants to distantly related C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{3}$\end{document} species so that changes directly responsible for the transition to C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{4}$\end{document} could not be distinguished from those that preceded or followed it. In this study, we analyze the genomes of 66 grass species, capturing the earliest diversification within Andropogoneae as well as their C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{3}$\end{document} relatives. Phylogenomics combined with molecular dating and analyses of protein evolution show that many changes linked to the evolution of C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{4}$\end{document} photosynthesis in Andropogoneae happened in the Early Miocene, between 21 and 18 Ma, after the split from its C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{3}$\end{document} sister lineage, and before the diversification of the group. This initial burst of changes was followed by an extended period of modifications to leaf anatomy and biochemistry during the diversification of Andropogoneae, so that a single C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{4}$\end{document} origin gave birth to a diversity of C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{4}$\end{document} phenotypes during 18 million years of speciation events and migration across geographic and ecological spaces. Our comprehensive approach and broad sampling of the diversity in the group reveals that one key transition can lead to a plethora of phenotypes following sustained adaptation of the ancestral state. [Adaptive evolution; complex traits; herbarium genomics; Jansenelleae; leaf anatomy; Poaceae; phylogenomics.]

The genomic timeline of cichlid fish diversification across continents

Cichlid fishes are celebrated for their vast taxonomic, phenotypic, and ecological diversity; however, a central aspect of their evolution - the timeline of their diversification - remains contentious. Here, we generate draft genome assemblies of 14 species representing the global cichlid diversity and integrate these into a new phylogenomic hypothesis of cichlid and teleost evolution that we time-calibrate with 58 re-evaluated fossil constraints and a new Bayesian model accounting for fossil-assignment uncertainty. Our results support cichlid diversification long after the breakup of the supercontinent Gondwana and lay the foundation for precise temporal reconstructions of the exceptional continental cichlid adaptive radiations.

Structure, gene order, and nucleotide composition of mitochondrial genomes in parasitic lice from Amblycera

Parasitic lice have unique mitochondrial (mt) genomes characterized by rearranged gene orders, variable genome structures, and less AT content compared to most other insects. However, relatively little is known about the mt genomes of Amblycera, the suborder sister to all other parasitic lice. Comparing among nine different genera (including representative of all seven families), we show that Amblycera have variable and highly rearranged mt genomes. Some genera have fragmented genomes that vary considerably in length, whereas others have a single mt chromosome. Notably, these genomes are more AT-biased than most other lice. We also recover genus-level phylogenetic relationships among Amblycera that are consistent with those reported from large nuclear datasets, indicating that mt sequences are reliable for reconstructing evolutionary relationships in Amblycera. However, gene order data cannot reliably recover these same relationships. Overall, our results suggest that the mt genomes of lice, already know to be distinctive, are even more variable than previously thought.

Disentangling lousy relationships: Comparative phylogenomics of two sucking louse lineages parasitizing chipmunks

The evolution of obligate parasites is often interpreted in light of their hosts' evolutionary history. An expanded approach is to examine the histories of multiple lineages of parasites that inhabit similar environments on a particular host lineage. Western North American chipmunks (genus Tamias) have a broad distribution, a history of divergence with gene flow, and host two species of sucking lice (Anoplura), Hoplopleura arboricola and Neohaematopinus pacificus. From total genomic sequencing, we obtained sequences of over 1100 loci sampled across the genomes of these lice to compare their evolutionary histories and examine the roles of host association in structuring louse relationships. Within each louse species, clades are largely associated with closely related chipmunk host species. Exceptions to this pattern appear to have a biogeographic component, but differ between the two louse species. Phylogenetic relationships among these major louse clades, in both species, are not congruent with chipmunk relationships. In the context of host associations, each louse lineage has a different evolutionary history, supporting the hypothesis that host-parasite assemblages vary both across the landscape and with the taxa under investigation. In addition, the louse Hoplopleura erratica (parasitizing the eastern Tamias striatus) is embedded within H. arboricola, rendering it paraphyletic. This phylogenetic result, together with comparable divergences within H. arboricola, indicate a need for taxonomic revision. Both host divergence and biogeographic components shape parasite diversification as demonstrated by the distinctive diversification patterns of these two independently evolving lineages that parasitize the same hosts.

Molecular Epidemiology of Acinetobacter calcoaceticus-Acinetobacter baumannii Complex Isolated From Children at the Hospital Infantil de México Federico Gómez

The Acinetobacter calcoaceticus-baumannii (Acb) complex is regarded as a group of phenotypically indistinguishable opportunistic pathogens responsible for mainly causing hospital-acquired pneumonia and bacteremia. The aim of this study was to determine the frequency of isolation of the species that constitute the Acb complex, as well as their susceptibility to antibiotics, and their distribution at the Hospital Infantil de Mexico Federico Gomez (HIMFG). A total of 88 strains previously identified by Vitek 2®, 40 as Acinetobacter baumannii and 48 as Acb complex were isolated from 52 children from 07, January 2015 to 28, September 2017. A. baumannii accounted for 89.77% (79/88) of the strains; Acinetobacter pittii, 6.82% (6/88); and Acinetobacter nosocomialis, 3.40% (3/88). Most strains were recovered mainly from patients in the intensive care unit (ICU) and emergency wards. Blood cultures (BC) provided 44.32% (39/88) of strains. The 13.63% (12/88) of strains were associated with primary bacteremia, 3.4% (3/88) with secondary bacteremia, and 2.3% (2/88) with pneumonia. In addition, 44.32% (39/88) were multidrug-resistant (MDR) strains and, 11.36% (10/88) were extensively drug-resistant (XDR). All strains amplified the bla_OXA-51 gene; 51.13% (45/88), the bla_OXA-23 gene; 4.54% (4/88), the bla_OXA-24 gene; and 2.27% (2/88), the bla_OXA-58 gene. Plasmid profiles showed that the strains had 1–6 plasmids. The strains were distributed in 52 pulsotypes, and 24 showed identical restriction patterns, with a correlation coefficient of 1.0. Notably, some strains with the same pulsotype were isolated from different patients, wards, or years, suggesting the persistence of more than one clone. Twenty-seven sequence types (STs) were determined for the strains based on a Pasteur multilocus sequence typing (MLST) scheme using massive sequencing; the most prevalent was ST 156 (27.27%, 24/88). The Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)-Cas I-Fb system provided amplification in A. baumannii and A. pittii strains (22.73%, 20/88). This study identified an increased number of MDR strains and the relationship among strains through molecular typing. The data suggest that more than one strain could be causing an infection in some patient. The implementation of molecular epidemiology allowed the characterization of a set of strains and identification of different attributes associated with its distribution in a specific environment.

"Parasite turnover zone" at secondary contact: A new pattern in host-parasite population genetics

We describe here a new pattern of population genetic structure in a host-parasite system that can arise after secondary contact of previously isolated populations. Due to different generation times, and therefore different tempos of molecular evolution, the host and parasite populations reach different degrees of genetic differentiation during their separation (e.g., in refugia). Consequently, upon secondary contact, the host populations are able to re-establish a single panmictic population across the area of contact, while the parasite populations stop their dispersal at the secondary contact zone and create a narrow hybrid zone. From the host's perspective, the parasite's hybrid zone functions on a microevolutionary scale as a "parasite turnover zone": while the hosts are passing from area A to area B, their parasites turn genetically from the area A genotypes to the area B genotypes. We demonstrate this novel pattern with a model composed of Apodemus mice and Polyplax lice by comparing maternally inherited markers (complete mitochondrial genomes, and complete genomes of the vertically transmitted symbiont Legionella polyplacis) with single nucleotide polymorphisms derived from louse genomic data. We discuss the circumstances that may lead to this pattern and possible reasons why it has been overlooked in studies of host-parasite population genetics.

Unraveling the Spiraling Radiation: A Phylogenomic Analysis of Neotropical Costus L

The family of pantropical spiral gingers (Costaceae Nakai; c. 125 spp.) can be used as a model to enhance our understanding of the mechanisms underlying Neotropical diversity. Costaceae has higher taxonomic diversity in South and Central America (c. 72 Neotropical species, c. 30 African, c. 23 Southeast Asian), particularly due to a radiation of Neotropical species of the genus Costus L. (c. 57 spp.). However, a well-supported phylogeny of the Neotropical spiral gingers including thorough sampling of proposed species encompassing their full morphologic and geographic variation is lacking, partly due to poor resolution recovered in previous analyses using a small sampling of loci. Here we use a phylogenomic approach to estimate the phylogeny of a sample of Neotropical Costus species using a targeted enrichment approach. Baits were designed to capture conserved elements' variable at the species level using available genomic sequences of Costus species and relatives. We obtained 832 loci (generating 791,954 aligned base pairs and 31,142 parsimony informative sites) for samples that encompassed the geographical and/or morphological diversity of some recognized species. Higher support values that improve the results of previous studies were obtained when including all the available loci, even those producing unresolved gene trees and having a low proportion of variable sites. Concatenation and coalescent-based species trees methods converge in almost the same topology suggesting a robust estimation of the relationships, even under the high levels of gene tree conflict presented here. The bait set design here presented made inferring a robust phylogeny to test taxonomic hypotheses possible and will improve our understanding of the origins of the charismatic diversity of the Neotropical spiral gingers.

Lousy grouse: Comparing evolutionary patterns in Alaska galliform lice to understand host evolution and host-parasite interactions

Understanding both sides of host-parasite relationships can provide more complete insights into host and parasite biology in natural systems. For example, phylogenetic and population genetic comparisons between a group of hosts and their closely associated parasites can reveal patterns of host dispersal, interspecies interactions, and population structure that might not be evident from host data alone. These comparisons are also useful for understanding factors that drive host-parasite coevolutionary patterns (e.g., codivergence or host switching) over different periods of time. However, few studies have compared the evolutionary histories between multiple groups of parasites from the same group of hosts at a regional geographic scale. Here, we used genomic data to compare phylogenomic and population genomic patterns of Alaska ptarmigan and grouse species (Aves: Tetraoninae) and two genera of their associated feather lice: Lagopoecus and Goniodes. We used whole-genome sequencing to obtain hundreds of genes and thousands of single-nucleotide polymorphisms (SNPs) for the lice and double-digest restriction-associated DNA sequences to obtain SNPs from Alaska populations of two species of ptarmigan. We found that both genera of lice have some codivergence with their galliform hosts, but these relationships are primarily characterized by host switching and phylogenetic incongruence. Population structure was also uncorrelated between the hosts and lice. These patterns suggest that grouse, and ptarmigan in particular, share habitats and have likely had historical and ongoing dispersal within Alaska. However, the two genera of lice also have sufficient dissimilarities in the relationships with their hosts to suggest there are other factors, such as differences in louse dispersal ability, that shape the evolutionary patterns with their hosts.

Degradation of key photosynthetic genes in the critically endangered semi-aquatic flowering plant Saniculiphyllum guangxiense (Saxifragaceae)

BACKGROUND

Plastid gene loss and pseudogenization has been widely documented in parasitic and mycoheterotrophic plants, which have relaxed selective constraints on photosynthetic function. More enigmatic are sporadic reports of pseudogenization and loss of important photosynthesis genes in lineages thought to be fully photosynthetic. Here we report the complete plastid genome of Saniculiphyllum guangxiense, a critically endangered and phylogenetically isolated plant lineage, along with genomic evidence of reduced chloroplast function. We also report 22 additional plastid genomes representing the diversity of its containing clade Saxifragales, characterizing gene content and placing variation in a broader phylogenetic context.

RESULTS

We find that the plastid genome of Saniculiphyllum has experienced pseudogenization of five genes of the ndh complex (ndhA, ndhB, ndhD, ndhF, and ndhK), previously reported in flowering plants with an aquatic habit, as well as the surprising pseudogenization of two genes more central to photosynthesis (ccsA and cemA), contrasting with strong phylogenetic conservatism of plastid gene content in all other sampled Saxifragales. These genes participate in photooxidative protection, cytochrome synthesis, and carbon uptake. Nuclear paralogs exist for all seven plastid pseudogenes, yet these are also unlikely to be functional.

CONCLUSIONS

Saniculiphyllum appears to represent the greatest degree of plastid gene loss observed to date in any fully photosynthetic lineage, perhaps related to its extreme habitat specialization, yet plastid genome length, structure, and substitution rate are within the variation previously reported for photosynthetic plants. These results highlight the increasingly appreciated dynamism of plastid genomes, otherwise highly conserved across a billion years of green plant evolution, in plants with highly specialized life history traits.

Mitochondrial genomes of Columbicola feather lice are highly fragmented, indicating repeated evolution of minicircle-type genomes in parasitic lice

Most animals have a conserved mitochondrial genome structure composed of a single chromosome. However, some organisms have their mitochondrial genes separated on several smaller circular or linear chromosomes. Highly fragmented circular chromosomes (“minicircles”) are especially prevalent in parasitic lice (Insecta: Phthiraptera), with 16 species known to have between nine and 20 mitochondrial minicircles per genome. All of these species belong to the same clade (mammalian lice), suggesting a single origin of drastic fragmentation. Nevertheless, other work indicates a lesser degree of fragmentation (2–3 chromosomes/genome) is present in some avian feather lice (Ischnocera: Philopteridae). In this study, we tested for minicircles in four species of the feather louse genus Columbicola (Philopteridae). Using whole genome shotgun sequence data, we applied three different bioinformatic approaches for assembling the Columbicola mitochondrial genome. We further confirmed these approaches by assembling the mitochondrial genome of Pediculus humanus from shotgun sequencing reads, a species known to have minicircles. Columbicola spp. genomes are highly fragmented into 15–17 minicircles between ∼1,100 and ∼3,100 bp in length, with 1–4 genes per minicircle. Subsequent annotation of the minicircles indicated that tRNA arrangements of minicircles varied substantially between species. These mitochondrial minicircles for species of Columbicola represent the first feather lice (Philopteridae) for which minicircles have been found in a full mitochondrial genome assembly. Combined with recent phylogenetic studies of parasitic lice, our results provide strong evidence that highly fragmented mitochondrial genomes, which are otherwise rare across the Tree of Life, evolved multiple times within parasitic lice.

Transcriptomics illuminate the phylogenetic backbone of tiger beetles

Phylogenomics is progressing rapidly, allowing large strides forward into our understanding of the tree of life. In this study, we generated transcriptomes from ethanol-preserved specimens of 13 tiger beetle species (Coleoptera: Cicindelinae) and one Scaritinae outgroup. From these 14 transcriptomes and seven publicly available transcriptomes, we recovered an average of 2538 loci for phylogenetic analysis. We constructed an evolutionary tree of tiger beetles to examine deep-level relationships and examined the extent to which the composition of the dataset, missing data, gene tree inconsistency and codon position saturation impacted phylogenetic accuracy. Ethanol-preserved specimens yielded similar numbers of loci to specimens originally preserved in costly reagents, showcasing more flexibility in transcriptomics than anticipated. The number of loci and gene tree inconsistency had less impact on downstream results than third codon position saturation and missing data. Our results recovered tiger beetles as sister to Carabidae with strong support, confirming their taxonomic status as an independent family within Adephaga. Within tiger beetles, phylogenetic relationships were robust across all nodes. This new phylogenomic backbone represents a useful framework for future endeavours in tiger beetle systematics and serves as a starting point for the development of less costly target capture toolkits to expand the taxonomic breadth of the future tiger beetle tree of life.

Putting the genome in insect phylogenomics

Next-generation sequencing technologies provide a substantial increase in the size of molecular phylogenetic datasets that can be obtained for studies of insect systematics. Several new genome reduction approaches are leveraging these technologies to generate large phylogenomic datasets: targeted amplicon sequencing, target capture, and transcriptome sequencing. Although cost effective, these approaches provide limited data for questions outside of phylogenetics. For many groups of insects, sequencing the entire genome at modest coverage is feasible. Using these genomic reads, an automated Target Restricted Assembly Method (aTRAM) can use the results of blast searches to assemble thousands of single copy ortholog genes across a group of interest. These locally assembled genes can then be compiled into very large phylogenomic datasets. These genomic libraries have the advantage in that they also contain reads from the mitochondrial genome and symbiont genomes, as well the entire insect genome, and can be leveraged for additional studies beyond phylogenetics.

Evolutionary genetics in insect phenotypic radiations: the value of a comparative genomic approach

Expanding genome sequencing and transgenic technologies are enabling the discovery of genes driving phenotypic diversity across insect taxa. Limitations in downstream functional genetic approaches, however, have been an obstacle for developing non-model systems for evolutionary genetics. Phenotypically diverse radiations, such as those exhibiting convergence and divergence as a result of mimicry, are ideal for evolutionary genetics as they can lead to insights using comparative genomic approaches alone. The varied and repeated instances of phenotypes in highly polymorphic systems allow assessment of whether similar loci are repeatedly targeted by selection and can inform how alleles sort across lineages. Comparative genomics of these taxa can be used to decipher components of gene regulatory networks, dissect regulatory regions, and validate genes.

Rates of niche and phenotype evolution lag behind diversification in a temperate radiation

Alternative models of evolutionary processes suggest different associations between species diversification and trait evolution, but limited empirical evidence is available to test these models across large clades at global extents. Here we investigate the relative timing of species diversification and niche and phenotypic evolution across a global plant radiation (Saxifragales) with enormous phenotypic and habitat variation. We demonstrate strong temporal lags among rates, with increased diversification occurring first, followed by niche and phenotype. Accelerated diversification rates are coincident with mid-Miocene expansion of temperate biomes. Later increases in niche and phenotypic evolutionary rates argue against density-dependent diversification alone, indicating a major role for ecological opportunity. These results have broad implications for understanding diversification processes and the origin of present-day temperate biotas.

Environmental change can create opportunities for increased rates of lineage diversification, but continued species accumulation has been hypothesized to lead to slowdowns via competitive exclusion and niche partitioning. Such density-dependent models imply tight linkages between diversification and trait evolution, but there are plausible alternative models. Little is known about the association between diversification and key ecological and phenotypic traits at broad phylogenetic and spatial scales. Do trait evolutionary rates coincide with rates of diversification, are there lags among these rates, or is diversification niche-neutral? To address these questions, we combine a deeply sampled phylogeny for a major flowering plant clade—Saxifragales—with phenotype and niche data to examine temporal patterns of evolutionary rates. The considerable phenotypic and habitat diversity of Saxifragales is greatest in temperate biomes. Global expansion of these habitats since the mid-Miocene provided ecological opportunities that, with density-dependent adaptive radiation, should result in simultaneous rate increases for diversification, niche, and phenotype, followed by decreases with habitat saturation. Instead, we find that these rates have significantly different timings, with increases in diversification occurring at the mid-Miocene Climatic Optimum (∼15 Mya), followed by increases in niche and phenotypic evolutionary rates by ∼5 Mya; all rates increase exponentially to the present. We attribute this surprising lack of temporal coincidence to initial niche-neutral diversification followed by ecological and phenotypic divergence coincident with more extreme cold and dry habitats that proliferated into the Pleistocene. A lack of density-dependence contrasts with investigations of other cosmopolitan lineages, suggesting alternative patterns may be common in the diversification of temperate lineages.

Genomic sequence capture of haemosporidian parasites: Methods and prospects for enhanced study of host-parasite evolution

Avian malaria and related haemosporidians (Plasmodium, [Para]Haemoproteus and Leucocytoozoon) represent an exciting multihost, multiparasite system in ecology and evolution. Global research in this field accelerated after the publication in 2000 of PCR protocols to sequence a haemosporidian mitochondrial (mtDNA) barcode and the development in 2009 of an open-access database to document the geographic and host ranges of parasite mtDNA haplotypes. Isolating haemosporidian nuclear DNA from bird hosts, however, has been technically challenging, slowing the transition to genomic-scale sequencing techniques. We extend a recently developed sequence capture method to obtain hundreds of haemosporidian nuclear loci from wild bird samples, which typically have low levels of infection, or parasitemia. We tested 51 infected birds from Peru and New Mexico and evaluated locus recovery in light of variation in parasitemia, divergence from reference sequences and pooling strategies. Our method was successful for samples with parasitemia as low as ~0.02% (2 of 10,000 blood cells infected) and mtDNA divergence as high as 15.9% (one Leucocytozoonsample), and using the most cost-effective pooling strategy tested. Phylogenetic relationships estimated with >300 nuclear loci were well resolved, providing substantial improvement over the mtDNA barcode. We provide protocols for sample preparation and sequence capture including custom probe sequences and describe our bioinformatics pipeline using atram 2.0, phyluce and custom Perl/Python scripts. This approach can be applied to thousands of avian samples that have already been found to have haemosporidian infections of at least moderate intensity, greatly improving our understanding of parasite speciation, biogeography and evolutionary dynamics.