PhyloMap: an algorithm for visualizing relationships of large sequence data sets and its application to the influenza A virus genome

Genetic Diversity in Natural Populations of the Near-Threatened Species Lignobrycon myersi (Characiformes, Triportheidae): Implications for Species Conservation

The river basins of Brazil contain a highly diverse ichthyofauna of remarkable endemism, including several threatened species. Accordingly, Lignobrycon myersi is a fish species distributed only in a few rivers from the state of Bahia, northeastern Brazil. Since this species is classified as Near Threatened and is poorly studied, efforts to understand the genetic structure of populations and putative cryptic forms should help define efficient strategies of management and conservation. Herein, the molecular identification and the population genetic diversity of specimens of L. myersi across their range (Almada, Contas, and Cachoeira river basins) were assessed using mitochondrial markers (16S rDNA and D-Loop, respectively). The inferences based on phylogenetics, genetic distance, and species delimitation methods invariably identified all samples as L. myersi. In addition, sequencing of D-loop fragments revealed significant haplotype diversity and a considerable level of population genetic structure. Despite their geographic isolation, these data suggested that populations from Almada and Contas rivers represent a single evolutionary lineage that could be managed as a whole. In contrast, the population from Cachoeira River was highly differentiated from the others and should be managed separately as a unique and endemic unit, particularly focused on the conservation of native habitats.

Taxonomic Delineation of the Old World Species Stomphastis thraustica (Lepidoptera: Gracillariidae) Feeding on Jatropha gossypiifolia (Euphorbiaceae) that Was Collected in the New World and Imported as a Biocontrol Agent to Australia

We provide the identification and species delineation of this biocontrol agent as Stomphastis thraustica (Meyrick in Trans Ent Soc Lond 80(1):107-120, 1908) belonging to the family Gracillariidae. We clarify the distribution pattern of S. thraustica, its host plant preferences, and present taxonomic and molecular diagnoses based on original morphological and genetic data as well as data retrieved from historic literature and genetic databases. Following our own collecting efforts in three continents Africa, South America, and Australia as well as our study of historic museum collection material, we present many new distribution records of S. thraustica for countries and territories in the world including the new discovery of this species in the Neotropical region and we report its introduction in Australia as a biocontrol agent. Using mitogenomic and COI gene data, we clarified that the closest relative of S. thraustica is Stomphastis sp. that occurs in Madagascar and Australia and feeds on the same host plant as S. thraustica - Jatropha gossypiifolia L. (Euphorbiaceae). The molecular sequence divergence in the mitochondrial DNA barcode fragment between these two closely related species S. thraustica and Stomphastis sp. is over 5.7% supporting that they are different species.

Environmental genomics of Late Pleistocene black bears and giant short-faced bears

Analysis of ancient environmental DNA (eDNA) has revolutionized our ability to describe biological communities in space and time,1-3 by allowing for parallel sequencing of DNA from all trophic levels.4-8 However, because environmental samples contain sparse and fragmented data from multiple individuals, and often contain closely related species,9 the field of ancient eDNA has so far been limited to organellar genomes in its contribution to population and phylogenetic studies.5,6,10,11 This is in contrast to data from fossils12,13 where full-genome studies are routine, despite these being rare and their destruction for sequencing undesirable.14-16 Here, we report the retrieval of three low-coverage (0.03×) environmental genomes from American black bear (Ursus americanus) and a 0.04× environmental genome of the extinct giant short-faced bear (Arctodus simus) from cave sediment samples from northern Mexico dated to 16-14 thousand calibrated years before present (cal kyr BP), which we contextualize with a new high-coverage (26×) and two lower-coverage giant short-faced bear genomes obtained from fossils recovered from Yukon Territory, Canada, which date to ∼22-50 cal kyr BP. We show that the Late Pleistocene black bear population in Mexico is ancestrally related to the present-day Eastern American black bear population, and that the extinct giant short-faced bears present in Mexico were deeply divergent from the earlier Beringian population. Our findings demonstrate the ability to separately analyze genomic-scale DNA sequences of closely related species co-preserved in environmental samples, which brings the use of ancient eDNA into the era of population genomics and phylogenetics.

DISSEQT-DIStribution-based modeling of SEQuence space Time dynamics

Rapidly evolving microbes are a challenge to model because of the volatile, complex, and dynamic nature of their populations. We developed the DISSEQT pipeline (DIStribution-based SEQuence space Time dynamics) for analyzing, visualizing, and predicting the evolution of heterogeneous biological populations in multidimensional genetic space, suited for population-based modeling of deep sequencing and high-throughput data. The pipeline is openly available on GitHub (https://github.com/rasmushenningsson/DISSEQT.jl, accessed 23 June 2019) and Synapse (https://www.synapse.org/#!Synapse: syn11425758, accessed 23 June 2019), covering the entire workflow from read alignment to visualization of results. Our pipeline is centered around robust dimension and model reduction algorithms for analysis of genotypic data with additional capabilities for including phenotypic features to explore dynamic genotype-phenotype maps. We illustrate its utility and capacity with examples from evolving RNA virus populations, which present one of the highest degrees of genetic heterogeneity within a given population found in nature. Using our pipeline, we empirically reconstruct the evolutionary trajectories of evolving populations in sequence space and genotype-phenotype fitness landscapes. We show that while sequence space is vastly multidimensional, the relevant genetic space of evolving microbial populations is of intrinsically low dimension. In addition, evolutionary trajectories of these populations can be faithfully monitored to identify the key minority genotypes contributing most to evolution. Finally, we show that empirical fitness landscapes, when reconstructed to include minority variants, can predict phenotype from genotype with high accuracy.

Molecular Data Reveal Multiple Lineages in Piranhas of the Genus Pygocentrus (Teleostei, Characiformes)

Carnivorous piranhas are distributed in four serrasalmid genera including Pygocentrus, which inhabit major river basins of South America. While P. cariba and P. piraya are endemics of the Orinoco and São Francisco basins, respectively, P. nattereri is widely distributed across the Amazonas, Essequibo, lower Paraná, Paraguay, and coastal rivers of northeastern Brazil, with recent records of introductions in Asia. Few studies have focused on the genetic diversity and systematics of Pygocentrus and the putative presence of additional species within P. nattereri has never been the subject of a detailed molecular study. Here we aimed to delimit species of Pygocentrus, test the phylogeographic structure of P. nattereri, and access the origin of introduced specimens of P. nattereri in Asia. Phylogenetic analyses based on a mitochondrial dataset involving maximum-likelihood tree reconstruction, genetic distances, Bayesian analysis, three delimitation approaches, and haplotype analysis corroborate the morphological hypothesis of the occurrence of three species of Pygocentrus. However, we provide here strong evidence that P. nattereri contains at least five phylogeographically-structured lineages in the Amazonas, Guaporé (type locality), Itapecuru, Paraná/Paraguay, and Tocantins/Araguaia river basins. We finally found that the introduced specimens in Asia consistently descend from the lineage of P. nattereri from the main Rio Amazonas. These results contribute to future research aimed to detect morphological variation that may occur in those genetic lineages of Pygocentrus.

Vector Quantized Spectral Clustering Applied to Whole Genome Sequences of Plants

We develop a Vector Quantized Spectral Clustering (VQSC) algorithm that is a combination of spectral clustering (SC) and vector quantization (VQ) sampling for grouping genome sequences of plants. The inspiration here is to use SC for its accuracy and VQ to make the algorithm computationally cheap (the complexity of SC is cubic in terms of the input size). Although the combination of SC and VQ is not new, the novelty of our work is in developing the crucial similarity matrix in SC as well as use of k-medoids in VQ, both adapted for the plant genome data. For Soybean, we compare our approach with commonly used techniques like Un-weighted Pair Graph Method with Arithmetic mean (UPGMA) and Neighbor Joining (NJ). Experimental results show that our VQSC outperforms both these techniques significantly in terms of cluster quality (average improvement of 21% over UPGMA and 24% over NJ) as well as time complexity (order of magnitude faster than both UPGMA and NJ).

An exhaustive phylogeny of the combtooth blenny genus Salaria (Pisces, Blenniidae) shows introgressive hybridization and lack of reciprocal mtDNA monophyly between the marine species Salaria basilisca and Salaria pavo

A comprehensive phylogeny of the genus Salaria based on mitochondrial and nuclear markers grouped the extant species of the genus in well-characterised marine and freshwater clades, thus rejecting the hypothesis of a polytypic origin of the freshwater Salaria populations and supporting the occurrence of a single invasion event of the inland waters by the genus. Based on both mitochondrial and nuclear DNA datasets, the Salaria species of the freshwater clade proved to be vicariant taxa originating from a common ancestor which could possibly spread throughout the circum-Mediterranean inland waters during the late Miocene Messinian salinity crisis, then experiencing a process of allopatric differentiation after the re-flooding of the Mediterranean basin. Within the marine clade, although the nuDNA datasets showed the existence of well-supported subclades in accordance to the morphological identification of the studied specimens, one of the two subclades obtained in the phylogenetic tree based on the mtDNA dataset included both S. basilisca and S. pavo specimens, thus failing to find the two species as reciprocally monophyletic. Such a mito-nuclear discordance is here ascribed to multiple mtDNA unidirectional introgression events from S. basilisca to S. pavo, and the molecular diversity pattern of the marine Salaria species is here ascribed to a Pleistocene speciation event nowadays partly concealed by the occurrence of introgressive hybridization phenomena between the two taxa. Our results urge for prudence when implementing DNA barcoding approaches since, in the presence of mito-nuclear discordance phenomena, single-marker mtDNA-only analyses might lead to significant misidentifications.

Cytogenetic analyses in Trinomys (Echimyidae, Rodentia), with description of new karyotypes

Trinomys Thomas (1921) is a terrestrial genus of spiny rats endemic to the Brazilian areas of Atlantic Forest and the transitional areas of Cerrado and Caatinga. Although most species have been already karyotyped, the available cytogenetic information is mostly restricted to diploid and fundamental numbers. We analyzed the chromosomes of two Trinomys species: Trinomys moojeni (2n = 56, FN = 106) and Trinomys setosus setosus (2n = 56, FN = 106 and 2n = 56, FN = 108). Our analyses included GTG- and CBG-banding, silver-staining of the nucleolar organizer regions, and chromosome mapping of telomeres and 45S rDNA by fluorescent in situ hybridization (FISH). Comparative GTG- and CBG-banding suggested that the interspecific variation may be due to rearrangements such as pericentric inversions, centromere repositioning, and heterochromatin variation. We report two new karyotypes for T. s. setosus and describe for the first time the banding patterns of the two Trinomys species.

Topo-phylogeny: Visualizing evolutionary relationships on a topographic landscape

Phylogenetic trees are the de facto standard for visualizing evolutionary relationships, but large trees can be difficult to interpret because they require a high cognitive load to identify relationships between multiple operational taxonomic units (OTUs). We present a new tool for displaying phylogenetic relationships as a topographic map in which OTUs autonomously attract or repel one another based on their individual branch lengths and distance to a common ancestor. This data visualization paradigm makes it possible to preattentively identify the nature of the relationship between items without having to trace a complex network of branches back to the root. This tool was developed for exploring phylogenetic data, but the technique could be extended for visualizing other hierarchical structures as well.

Cryptic Species or Inadequate Taxonomy? Implementation of 2D Geometric Morphometrics Based on Integumental Organs as Landmarks for Delimitation and Description of Copepod Taxa

Discovery of cryptic species using molecular tools has become common in many animal groups but it is rarely accompanied by morphological revision, creating ongoing problems in taxonomy and conservation. In copepods, cryptic species have been discovered in most groups where fast-evolving molecular markers were employed. In this study at Yeelirrie in Western Australia we investigate a subterranean species complex belonging to the harpacticoid genus Schizopera Sars, 1905, using both the barcoding mitochondrial COI gene and landmark-based two-dimensional geometric morphometrics. Integumental organs (sensilla and pores) are used as landmarks for the first time in any crustacean group. Complete congruence between DNA-based species delimitation and relative position of integumental organs in two independent morphological structures suggests the existence of three distinct evolutionary units. We describe two of them as new species, employing a condensed taxonomic format appropriate for cryptic species. We argue that many supposedly cryptic species might not be cryptic if researchers focus on analyzing morphological structures with multivariate tools that explicitly take into account geometry of the phenotype. A perceived supremacy of molecular methods in detecting cryptic species is in our view a consequence of disparity of investment and unexploited recent advancements in morphometrics among taxonomists. Our study shows that morphometric data alone could be used to find diagnostic morphological traits and gives hope to anyone studying small animals with a hard integument or shell, especially opening the door to assessing fossil diversity and rich museum collections. We expect that simultaneous use of molecular tools with geometry-oriented morphometrics may yield faster formal description of species. Decrypted species in this study are a good example for urgency of formal descriptions, as they display short-range endemism in small groundwater calcrete aquifers in a paleochannel, where their conservation may be threatened by proposed mining.

A new look at an old virus: patterns of mutation accumulation in the human H1N1 influenza virus since 1918

BACKGROUND

The H1N1 influenza A virus has been circulating in the human population for over 95 years, first manifesting itself in the pandemic of 1917-1918. Initial mortality was extremely high, but dropped exponentially over time. Influenza viruses have high mutation rates, and H1N1 has undergone significant genetic changes since 1918. The exact nature of H1N1 mutation accumulation over time has not been fully explored.

METHODS

We have made a comprehensive historical analysis of mutational changes within H1N1 by examining over 4100 fully-sequenced H1N1 genomes. This has allowed us to examine the genetic changes arising within H1N1 from 1918 to the present.

RESULTS

We document multiple extinction events, including the previously known extinction of the human H1N1 lineage in the 1950s, and an apparent second extinction of the human H1N1 lineage in 2009. These extinctions appear to be due to a continuous accumulation of mutations. At the time of its disappearance in 2009, the human H1N1 lineage had accumulated over 1400 point mutations (more than 10% of the genome), including approximately 330 non-synonymous changes (7.4% of all codons). The accumulation of both point mutations and non-synonymous amino acid changes occurred at constant rates (μ = 14.4 and 2.4 new mutations/year, respectively), and mutations accumulated uniformly across the entire influenza genome. We observed a continuous erosion over time of codon-specificity in H1N1, including a shift away from host (human, swine, and bird [duck]) codon preference patterns.

CONCLUSIONS

While there have been numerous adaptations within the H1N1 genome, most of the genetic changes we document here appear to be non-adaptive, and much of the change appears to be degenerative. We suggest H1N1 has been undergoing natural genetic attenuation, and that significant attenuation may even occur during a single pandemic. This process may play a role in natural pandemic cessation and has apparently contributed to the exponential decline in mortality rates over time, as seen in all major human influenza strains. These findings may be relevant to the development of strategies for managing influenza pandemics and strain evolution.