Genomics clarifies taxonomic boundaries in a difficult species complex

PAPipe: A Pipeline for Comprehensive Population Genetic Analysis

Advancements in next-generation sequencing (NGS) technologies have led to a substantial increase in the availability of population genetic variant data, thus prompting the development of various population analysis tools to enhance our understanding of population structure and evolution. The tools that are currently used to analyze population genetic variant data generally require different environments, parameters, and formats of the input data, which can act as a barrier preventing the wide-spread usage of such tools by general researchers who may not be familiar with bioinformatics. To address this problem, we have developed an automated and comprehensive pipeline called PAPipe to perform nine widely used population genetic analyses using population NGS data. PAPipe seamlessly interconnects and serializes multiple steps, such as read trimming and mapping, genetic variant calling, data filtering, and format converting, along with nine population genetic analyses such as principal component analysis, phylogenetic analysis, population tree analysis, population structure analysis, linkage disequilibrium decay analysis, selective sweep analysis, population admixture analysis, sequentially Markovian coalescent analysis, and fixation index analysis. PAPipe also provides an easy-to-use web interface that allows for the parameters to be set and the analysis results to be browsed in intuitive manner. PAPipe can be used to generate extensive results that provide insights that can help enhance user convenience and data usability. PAPipe is freely available at https://github.com/jkimlab/PAPipe.

Plotting for change: an analytical framework to aid decisions on which lineages are candidate species in phylogenomic species discovery

A recent study argued that coalescent-based models of species delimitation mostly delineate population structure, not species, and called for the validation of candidate species using biological information additional to the genetic information, such as phenotypic or ecological data. Here, we introduce a framework to interrogate genomic datasets and coalescent-based species trees for the presence of candidate species in situations where additional biological data are unavailable, unobtainable or uninformative. For de novo genomic studies of species boundaries, we propose six steps: (1) visualize genetic affinities among individuals to identify both discrete and admixed genetic groups from first principles and to hold aside individuals involved in contemporary admixture for independent consideration; (2) apply phylogenetic techniques to identify lineages; (3) assess diagnosability of those lineages as potential candidate species; (4) interpret the diagnosable lineages in a geographical context (sympatry, parapatry, allopatry); (5) assess significance of difference or trends in the context of sampling intensity; and (6) adopt a holistic approach to available evidence to inform decisions on species status in the difficult cases of allopatry. We adopt this approach to distinguish candidate species from within-species lineages for a widespread species complex of Australian freshwater fishes (Retropinna spp.). Our framework addresses two cornerstone issues in systematics that are often not discussed explicitly in genomic species discovery: diagnosability and how to determine it, and what criteria should be used to decide whether diagnosable lineages are conspecific or represent different species.

Phylogeography of the iconic Australian pink cockatoo, Lophochroa leadbeateri

The pink cockatoo (Lophochroa leadbeateri; or Major Mitchell’s cockatoo) is one of Australia’s most iconic bird species. Two subspecies based on morphology are separated by a biogeographical divide, the Eyrean Barrier. Testing the genetic basis for this subspecies delineation, clarifying barriers to gene flow and identifying any cryptic genetic diversity will likely have important implications for conservation and management. Here, we used genome-wide single nucleotide polymorphisms (SNPs) and mitochondrial DNA data to conduct the first range-wide genetic assessment of the species. The aims were to investigate the phylogeography of the pink cockatoo, to characterize conservation units and to reassess subspecies boundaries. We found consistent but weak genetic structure between the two subspecies based on nuclear SNPs. However, phylogenetic analysis of nuclear SNPs and mitochondrial DNA sequence data did not recover reciprocally monophyletic groups, indicating incomplete evolutionary separation between the subspecies. Consequently, we have proposed that the two currently recognized subspecies be treated as separate management units rather than evolutionarily significant units. Given that poaching is suspected to be a threat to this species, we assessed the utility of our data for wildlife forensic applications. We demonstrated that a subspecies identification test could be designed using as few as 20 SNPs.

Genomic and morphological evidence of distinct populations in the endemic common (weedy) seadragon Phyllopteryx taeniolatus (Syngnathidae) along the east coast of Australia

The common or weedy seadragon, Phyllopteryx taeniolatus, is an iconic and endemic fish found across temperate reefs of southern Australia. Despite its charismatic nature, few studies have been published, and the extent of population sub-structuring remains poorly resolved. Here we used 7462 single nucleotide polymorphisms (SNPs) to identify the extent of population structure in the weedy seadragon along the temperate southeast coast of Australia. We identified four populations, with strong genetic structure (F_ST = 0.562) between them. Both Discriminant Analysis of Principle Components (DAPC) and Bayesian clustering analyses support four distinct genetic clusters (north to south: central New South Wales, southern NSW, Victoria and Tasmania). In addition to these genetic differences, geographical variation in external morphology was recorded, with individuals from New South Wales shaped differently for a few measurements to those from the Mornington Peninsula (Victoria). We posit that these genetic and morphological differences suggest that the Victorian population of P. taeniolatus was historically isolated by the Bassian Isthmus during the last glacial maximum and should now be considered at least a distinct population. We also recorded high levels of genetic structure among the other locations. Based on the genomic and to a degree morphological evidence presented in this study, we recommend that the Victorian population be managed separately from the eastern populations (New South Wales and Tasmania).

Two New Aspergillus flavus Reference Genomes Reveal a Large Insertion Potentially Contributing to Isolate Stress Tolerance and Aflatoxin Production

Efforts in genome sequencing in the Aspergillus genus have led to the development of quality reference genomes for several important species including A. nidulans, A. fumigatus, and A. oryzae However, less progress has been made for A. flavus As part of the effort of the USDA-ARS Annual Aflatoxin Workshop Fungal Genome Project, the isolate NRRL3357 was sequenced and resulted in a scaffold-level genome released in 2005. Our goal has been biologically driven, focusing on two areas: isolate variation in aflatoxin production and drought stress exacerbating aflatoxin production by A. flavus Therefore, we developed two reference pseudomolecule genome assemblies derived from chromosome arms for two isolates: AF13, a MAT1-2, highly stress tolerant, and highly aflatoxigenic isolate; and NRRL3357, a MAT1-1, less stress tolerant, and moderate aflatoxin producer in comparison to AF13. Here, we report these two reference-grade assemblies for these isolates through a combination of PacBio long-read sequencing and optical mapping, and coupled them with comparative, functional, and phylogenetic analyses. This analysis resulted in the identification of 153 and 45 unique genes in AF13 and NRRL3357, respectively. We also confirmed the presence of a unique 310 Kb insertion in AF13 containing 60 genes. Analysis of this insertion revealed the presence of a bZIP transcription factor, named atfC, which may contribute to isolate pathogenicity and stress tolerance. Phylogenomic analyses comparing these and other available assemblies also suggest that the species complex of A. flavus is polyphyletic.

Phylogenomic Analyses Clarify True Species within the Butterfly Genus Speyeria despite Evidence of a Recent Adaptive Radiation

When confronted with an adaptive radiation, considerable evidence is needed to resolve the evolutionary relationships of these closely related lineages. The North American genus Speyeria is one especially challenging radiation of butterflies due to potential signs of incomplete lineage sorting, ongoing hybridization, and similar morphological characters between species. Previous studies have found species to be paraphyletic and have been unable to disentangle taxa, often due to a lack of data and/or incomplete sampling. As a result, Speyeria remains unresolved. To achieve phylogenetic resolution of the genus, we conducted phylogenomic and population genomic analyses of all currently recognized North American Speyeria species, as well as several subspecies, using restriction-site-associated DNA sequencing (RADseq). Together, these analyses confirm the 16 canonical species, and clarify many internal relationships. However, a few relationships within Speyeria were poorly supported depending on the evolutionary model applied. This lack of resolution among certain taxa corroborates Speyeria is experiencing an ongoing adaptive radiation, with incomplete lineage sorting and lack of postzygotic reproductive barriers contributing to hybridization and further ambiguity. Given that many Speyeria taxa are under duress from anthropogenic factors, their legal protection must be viewed cautiously and on a case by case basis in order to properly conserve the diversity being generated.

Fish diversification at the pace of geomorphological changes: evolutionary history of western Iberian Leuciscinae (Teleostei: Leuciscidae) inferred from multilocus sequence data

The evolutionary history of western Iberian Leuciscinae, obligatory freshwater fish, is directly linked to the evolution of the hydrographic network of the Iberian Peninsula after its isolation from the rest of Europe, which involved dramatic rearrangements such as the transition from endorheic lakes to open basins draining to the Atlantic. Previous phylogenetic research on western Iberian leuciscines, using mainly mitochondrial DNA and more recently one or two nuclear genes, has found contradictory results and there remain many unresolved issues regarding species relationships, taxonomy, and evolutionary history. Moreover, there is a lack of integration between phylogenetic and divergence time estimates and information on the timing of geomorphological changes and paleobasin rearrangements in the Iberian Peninsula. This study presents the first comprehensive fossil-calibrated multilocus coalescent species tree of western Iberian Leuciscinae (including 14 species of Achondrostoma, Iberochondrostoma, Pseudochondrostoma and Squalius endemic to the Iberian Peninsula, seven of which endemic to Portugal) based on seven nuclear genes, and from which we infer their biogeographic history by comparing divergence time estimates to known dated geological events. The phylogenetic pattern suggests slow-paced evolution of leuciscines during the Early-Middle Miocene endorheic phase of the main Iberian river basins, with the shift to exorheism in the late Neogene-Quaternary allowing westward dispersals that resulted in many cladogenetic events and a high rate of endemism in western Iberia. The results of this study also: (i) confirm the paraphyly of S. pyrenaicus with respect to S. carolitertii, and thus the possible presence of a new taxon in the Portuguese Tagus currently assigned to S. pyrenaicus; (ii) support the taxonomic separation of the Guadiana and Sado populations of S. pyrenaicus; (iii) show the need for further population sampling and taxonomic research to clarify the phylogenetic status of A. arcasii from the Minho basin and of the I. lusitanicum populations in the Sado and Tagus basins; and (iv) indicate that A. occidentale, I. olisiponensis and P. duriensis are the most ancient lineages within their respective genera.

Exploring malaria vector diversity on the Amazon Frontier

BACKGROUND

Deforestation in the Amazon and the social vulnerability of its settler communities has been associated with increased malaria incidence. The feeding biology of the most important malaria vectors in the region, notably Nyssorhynchus darlingi, compounds efforts to control vectors and reduce transmission of what has become known as "Frontier Malaria". Exploring Anophelinae mosquito diversity is fundamental to understanding the species responsible for transmission and developing appropriate management and intervention strategies for malaria control in the Amazon River basin.

METHODS

This study describes Anophelinae mosquito diversity from settler communities affected by Frontier Malaria in the states of Acre, Amazonas and Rondônia by analysing COI gene data using cluster and tree-based species delimitation approaches.

RESULTS

In total, 270 specimens from collection sites were sequenced and these were combined with 151 reference (GenBank) sequences in the analysis to assist in species identification. Conservative estimates found that the number of species collected at these sites was between 23 (mPTP partition) and 27 (strict ABGD partition) species, up to 13 of which appeared to be new. Nyssorhynchus triannulatus and Nyssorhynchus braziliensis displayed exceptional levels of intraspecific genetic diversity but there was little to no support for putative species complex status.

CONCLUSIONS

This study demonstrates that Anophelinae mosquito diversity continues to be underestimated in poorly sampled areas where frontier malaria is a major public health concern. The findings will help shape future studies of vector incrimination and transmission dynamics in these areas and support efforts to develop more effective vector control and transmission reduction strategies in settler communities in the Amazon River basin.

Salmo macrostigma (Teleostei, Salmonidae): Nothing more than a brown trout (S. trutta) lineage?

We examined specimens of the macrostigma trout Salmo macrostigma, which refers to big black spots on the flanks, to assess whether it is an example of taxonomic inflation within the brown trout Salmo trutta complex. Using new specimens, publicly available data and a mitogenomic protocol to amplify the control and cytochrome b regions of the mitochondrial genome from degraded museum samples, including one syntype specimen, the present study shows that the macrostigma trout is not a valid species. Our results suggest the occurrence of a distinct evolutionary lineage of S. trutta in North Africa and Sicily. The name of the North African lineage is proposed for this lineage, which was found to be sister to the Atlantic lineage of brown trout, S. trutta.