Inferring phylogenetic structure, hybridization and divergence times within Salmoninae (Teleostei: Salmonidae) using RAD-sequencing

Phylogenomic analysis of the Lake Kronotskoe species flock of Dolly Varden charr reveals genetic and developmental signatures of sympatric radiation

Recent adaptive radiations provide evolutionary case studies, which provide the context to parse the relationship between genomic variation and the origins of distinct phenotypes. Sympatric radiations of the charr complex (genus Salvelinus) present a trove for phylogenetics as charrs have repeatedly diversified into multiple morphs with distinct feeding specializations. However, species flocks normally comprise only two to three lineages. Dolly Varden charr inhabiting Lake Kronotske represent the most extensive radiation described for the charr genus, containing at least seven lineages, each with defining morphological and ecological traits. Here, we perform the first genome-wide analysis of this species flock to parse the foundations of adaptive change. Our data support distinct, reproductively isolated lineages with little evidence of hybridization. We also find that specific selection on thyroid signaling and craniofacial genes forms a genomic basis for the radiation. Thyroid hormone is further implicated in subsequent lineage partitioning events. These results delineate a clear genetic basis for the diversification of specialized lineages, and highlight the role of developmental mechanisms in shaping the forms generated during adaptive radiation.

A new approach to exploratory data analysis in hominin phylogenetic reconstruction

The phylogenetic relationships between fossil hominin taxa have been a contentious topic for decades. Recent discoveries of new taxa, rather than resolving the issue, have only further confused it. Compounding this problem are the limitations of some of the tools frequently used by paleoanthropologists to analyze these relationships. Most commonly, phylogenetic questions are investigated using analytical methods such as maximum parsimony and Bayesian analysis. While these are useful analytical tools, these tree-building methods can have limitations when investigating taxa that may have complex evolutionary histories. Exploratory data analysis can provide information about patterns in a dataset that are obscured by tree-based methods. These patterns include phylogenetic signal conflict, which is not depicted in tree-based methods. Signal conflict can have a number of sources, including methodological issues with character choice, taxonomic issues, homoplasy, and gene flow between taxa. In this study, an exploratory data analysis of fossil hominin morphological data is conducted using the tree-based analytical method neighbor-joining and the network-based analytical method neighbor-net with the goal of visualizing phylogenetic signal conflict within a hominin morphological data set. The data set is divided into cranial regions, and each cranial region is analyzed individually to investigate which regions of the skull contain the highest levels of signal conflict. Results of this analysis show that conflicting phylogenetic signals are present in the hominin fossil record during the relatively speciose period between 3 and 1 Ma, and they also indicate that levels of signal conflict vary by cranial region. Possible sources of these conflicting signals are then explored. Exploratory data analyses such as this can be a useful tool in generating phylogenetic hypotheses and in refining character choice. This study also highlights the value network-based approaches can bring to the hominin phylogenetic analysis toolkit.

Can cryptic female choice prevent invasive hybridization in external fertilizing fish?

Polyandrous mating systems result in females mating with multiple males, generating opportunities for strong pre-mating and post-mating sexual selection. Polyandry also creates the potential for unintended matings and subsequent sperm competition with hybridizing species. Cryptic female choice allows females to bias paternity towards preferred males under sperm competition and may include conspecific sperm preference when under hybridization risk. The potential for hybridization becomes particularly important in context of invasive species that can novelly hybridize with natives, and by definition, have evolved allopatrically. We provide the first examination of conspecific sperm preference in a system of three species with the potential to hybridize: North American native Atlantic salmon (Salmo salar) and brook char (Salvelinus fontinalis), and invasive brown trout (Salmo trutta) from Europe. Using naturalized populations on the island of Newfoundland, we measured changes in sperm swimming performance, a known predictor of paternity, to determine the degree of modification in sperm swimming to female cues related to conspecific sperm preference. Compared to water alone, female ovarian fluid in general had a pronounced effect and changed sperm motility (by a mean of 53%) and swimming velocity (mean 30%), but not linearity (mean 6%). However, patterns in the degree of modification suggest there is no conspecific sperm preference in the North American populations. Furthermore, female cues from both native species tended to boost the sperm of invasive males more than their own. We conclude that cryptic female choice via ovarian fluid mediated sperm swimming modification is too weak in this system to prevent invasive hybridization and is likely insufficient to promote or maintain reproductive isolation between the native North American species.

The systematics and evolution of the Sri Lankan rainforest land snail Corilla: New insights from RADseq-based phylogenetics

The stylommatophoran land-snail genus Corilla is endemic to Sri Lanka and India's Western Ghats. On the basis of habitat distribution and shell morphology, the 10 extant Sri Lankan species fall into two distinct groups, lowland and montane. Here, we use phylogenetic analyses of restriction-site-associated DNA sequencing (RADseq) data and ancestral-state reconstructions of habitat association and shell morphology to clarify the systematics and evolution of Sri Lankan Corilla. Our dataset consists of 9 species of Corilla. Phylogenetic analyses were based on 88 assemblies (9,604-4,132,850 bp) generated by the RADseq assembler ipyrad, using four parameter combinations and different levels of missing data. Trees were inferred using a maximum likelihood (ML) approach. Ancestral states were reconstructed using maximum parsimony (MP) and ML approaches, with 1 binary state character analysed for habitat association (lowland vs montane) and 6 binary state characters analysed for shell morphology (shape, colour, lip width, length of upper palatal folds, orientation of upper palatal folds and collabral sculpture). Over a wide range of missing data (40-87 % missing individuals per locus) and assembly sizes (62,279-4,132,850 bp), nearly all trees conformed to one of two topologies (A and B), most relationships were strongly supported and total branch support approached the maximal value. Apart from the position of Corilla odontophora 'south', topologies A and B showed similar, well-resolved relationships at and above the species level. Our study agrees with the shell-based taxonomy of C. adamsi, C. beddomeae, C. carabinata, C. colletti and C. humberti (all maximally supported as monophyletic species). It shows that C. erronea and C. fryae constitute a single relatively widespread species (for which the valid name is C. erronea) and that the names C. gudei and C. odontophora each apply to at least two distinct, yet conchologically-cryptic species. The MP and ML ancestral-state reconstructions yielded broadly similar results and provide firm evidence that diversification in Sri Lankan Corilla has involved evolutionary convergence in the shell morphology of lowland lineages, with a pale shell and wide lip having evolved on at least two separate occasions (in C. carabinata and C. colletti) from montane ancestors having a dark, narrow-lipped shell.

Analysis of microbiota-host communication mediated by butyrate in Atlantic Salmon

Butyrate is a microbiota-produced metabolite, sensed by host short-chain fatty acid receptors FFAR2 (Gpr43), FFAR3 (Gpr41), HCAR2 (Gpr109A), and Histone deacetylase (HDAC) that promotes microbiota-host crosstalk. Butyrate influences energy uptake, developmental and immune response in mammals. This microbial metabolite is produced by around 79 anaerobic genera present in the mammalian gut, yet little is known about the role of butyrate in the host-microbiota interaction in salmonid fish. To further our knowledge of this interaction, we analyzed the intestinal microbiota and genome of Atlantic salmon (Salmo salar), searching for butyrate-producing genera and host butyrate receptors. We identified Firmicutes, Proteobacteria, and Actinobacteria as the main butyrate-producing bacteria in the salmon gut microbiota. In the Atlantic salmon genome, we identified an expansion of genes orthologous to FFAR2 and HCAR2 receptors, and class I and IIa HDACs that are sensitive to butyrate. In addition, we determined the expression levels of orthologous of HCAR2 in the gut, spleen, and head-kidney, and FFAR2 in RTgutGC cells. The effect of butyrate on the Atlantic salmon immune response was evaluated by analyzing the pro and anti-inflammatory cytokines response in vitro in SHK-1 cells by RT-qPCR. Butyrate decreased the expression of the pro-inflammatory cytokine IL-1β and increased anti-inflammatory IL-10 and TGF-β cytokines. Butyrate also reduced the expression of interferon-alpha, Mx, and PKR, and decreased the viral load at a higher concentration (4 mM) in cells treated with this molecule before the infection with Infectious Pancreatic Necrosis Virus (IPNV) by mechanisms independent of FFAR2, FFAR3 and HCAR2 expression that probably inhibit HDAC. Moreover, butyrate modified phosphorylation of cytoplasmic proteins in RTgutGC cells. Our data allow us to infer that Atlantic salmon have the ability to sense butyrate produced by their gut microbiota via different specific targets, through which butyrate modulates the immune response of pro and anti-inflammatory cytokines and the antiviral response.

Evaluating a species phylogeny using ddRAD SNPs: Cyto-nuclear discordance and introgression in the salmonid genus Thymallus (Salmonidae)

Hybridization and introgression are very common among freshwater fishes due to the dynamic nature of hydrological landscapes. Cyclic patterns of allopatry and secondary contact provide numerous opportunities for interspecific gene flow, which can lead to discordant paths of evolution for mitochondrial and nuclear genomes. Here, we used double digest restriction-site associated DNA sequencing (ddRADseq) to obtain a genome-wide single nucleotide polymorphism (SNP) dataset comprehensive for allThymallus (Salmonidae)species to infer phylogenetic relationships and evaluate potential recent and historical gene flow among species. The newly obtained nuclear phylogeny was largely concordant with a previously published mitogenome-based topology but revealed a few cyto-nuclear discordances. These incongruencies primarily involved the placement of internal nodes rather than the resolution of species, except for one European species where anthropogenic stock transfers are thought to be responsible for the observed pattern. The analysis of four contact zones where multiple species are found revealed a few cases of mitochondrial capture and limited signals of nuclear introgression. Interestingly, the mechanisms restricting interspecific gene flow might be distinct; while in zones of secondary contact, small-scale physical habitat separation appeared as a limiting factor, biologically based reinforcement mechanisms are presumed to be operative in areas where species presumably evolved in sympatry. Signals of historical introgression were largely congruent with the routes of species dispersal previously inferred from mitogenome data. Overall, the ddRADseq dataset provided a robust phylogenetic reconstruction of the genus Thymallus including new insights into historical hybridization and introgression, opening up new questions concerning their evolutionary history.

Salmonidae Genome: Features, Evolutionary and Phylogenetic Characteristics

The salmon family is one of the most iconic and economically important fish families, primarily possessing meat of excellent taste as well as irreplaceable nutritional and biological value. One of the most common and, therefore, highly significant members of this family, the Atlantic salmon (Salmo salar L.), was not without reason one of the first fish species for which a high-quality reference genome assembly was produced and published. Genomic advancements are becoming increasingly essential in both the genetic enhancement of farmed salmon and the conservation of wild salmon stocks. The salmon genome has also played a significant role in influencing our comprehension of the evolutionary and functional ramifications of the ancestral whole-genome duplication event shared by all Salmonidae species. Here we provide an overview of the current state of research on the genomics and phylogeny of the various most studied subfamilies, genera, and individual salmonid species, focusing on those studies that aim to advance our understanding of salmonid ecology, physiology, and evolution, particularly for the purpose of improving aquaculture production. This review should make potential researchers pay attention to the current state of research on the salmonid genome, which should potentially attract interest in this important problem, and hence the application of new technologies (such as genome editing) in uncovering the genetic and evolutionary features of salmoniforms that underlie functional variation in traits of commercial and scientific importance.

Identification and quantification of chimeric sequencing reads in a highly multiplexed RAD-seq protocol

Highly multiplexed approaches have become common in genomic studies. They have improved the cost-effectiveness of genotyping hundreds of individuals using combinatorially barcoded adapters. These strategies, however, can potentially misassigned reads to incorrect samples. Here, we used a modified quaddRAD protocol to analyse the occurrence of index hopping and PCR chimeras in a series of experiments with up to 100 multiplexed samples per sequencing lane (639 samples in total). We created two types of sequencing libraries: four libraries of type A, where PCRs were run on individual samples before multiplexing, and three libraries of type B, where PCRs were run on pooled samples. We used fixed pairs of inner barcodes to identify chimeric reads. Type B libraries show a higher percentage of misassigned reads (1.15%) than type A libraries (0.65%). We also quantify the commonly undetectable chimeric sequences that occur whenever multiplexed groups of samples with different outer barcodes are sequenced together on a single flow cell. Our results suggest that these types of chimeric sequences represent up to 1.56% and 1.29% of reads in type A and B libraries, respectively. We also show that increasing the number of mismatches allowed for barcode rescue to above 2 dramatically increases the number of recovered chimeric reads. We provide recommendations for developing highly multiplexed RAD-seq protocols and analysing the resulting data to minimize the generation of chimeric sequences, allowing their quantification and a finer control on the number of PCR cycles necessary to generate enough input DNA for library preparation.

Phylogenomic analysis of the bowfin (Amia calva) reveals unrecognized species diversity in a living fossil lineage

The Bowfin (Amia calva), as currently recognized, represents the sole living member of the family Amiidae, which dates back to approximately 150 Ma. Prior to 1896, 13 species of extant Bowfins had been described, but these were all placed into a single species with no rationale or analysis given. This situation has persisted until the present day, with little attention given to re-evaluation of those previously described nominal forms. Here, we present a phylogenomic analysis based on over 21,000 single nucleotide polymorphisms (SNPs) from 94 individuals that unambiguously demonstrates the presence of at least two independent evolutionary lineages within extant Amia populations that merit species-level standing, as well as the possibility of two more. These findings not only expand the recognizable species diversity in an iconic, ancient lineage, but also demonstrate the utility of such methods in addressing previously intractable questions of molecular systematics and phylogeography in slowly evolving groups of ancient fishes.

The Italian Alpine and Subalpine trouts: Taxonomy, Evolution, and Conservation

Simple Summary

In a great part of the world, trout fishing has long inspired human spiritual ideals of immersion into nature and recreation, far removed from the fast-encroaching urbanization. Concurrently, these values and emotions fueled a white-hot business, establishing a florid market of outdoor recreation. Since the 20th century, the trout-culture industry strived to provide anglers with fishing entertainment by stocking massive amounts of non-native trouts in dozens of countries, irrespective of the lakes’ and rivers’ carrying capacity. This had dire consequences on the structural and functional diversity of these ecosystems. “Trout wars” sparked throughout the world between the promoters of stocking activities and the promoters of “wild trout management” and ethics. The “Italian trout war” has been fought on the harsh battleground of trout taxonomy, ecology, distribution, and native vs. non-native interfertile species. Northern Italy, home to the Italian Alpine and subalpine trouts and economic center of the national trout-culture and stocking industry, was particularly affected by this clash. We review here the state of art of this ongoing debate, outlining our scientific view of the taxonomy, evolution, distribution, and sustainable management of the native Italian trouts of northern Italy.

Abstract

During the last 150 years, the trout-culture industry focused on enhancing trout populations by stocking, in response to the growing anglers’ demand and the habitat degradation associated to the rapid urbanization and hydropower development. The industrialized north of Italy, home to the Italian Alpine and subalpine trout populations, is the source of most of the revenues of the national trout-culture industry. Its rapid growth, and the massive introduction of non-native interfertile trouts eroded the genetic diversity of native lineages, leading to harsh confrontations between scientists, institutions, and sportfishing associations. We review here the state of the art of the taxonomy and distribution of the northern Italian native trouts, presenting both scientific results and historical documentation. We think the only native trouts in this region are Salmo marmoratus, widespread in this region, plus small and fragmented populations of S. ghigii, present only in the South-western Alps. We strongly recommend the interruption of stocking of domesticated interfertile non-native trouts in this area, and recommend the adoption of Evolutionary Significant Units for salmonid fishery management. We further propose future research directions for a sustainable approach to the conservation and ecosystem management of the fishery resources and inland waters of northern Italy.

Genetic and Phenotypic Characteristics of the Salmo trutta Complex in Italy

Salmonid fish have become ecological and research models of study in the field of conservation genetics and genomics. Over the last decade, brown trout have received a high level of interest in research and publications. The term Salmo trutta complex is used to indicate the large number of geographic forms present in the species Salmo trutta. In Europe, the S. trutta complex consists (based on mitochondrial DNA control region analysis) of seven major evolutionary lineages: Atlantic (AT), Mediterranean (ME), Adriatic (AD), Danubian (DA), Marmoratus (MA), Duero (DU) and Tigris (TI). In several nations, the difficulty of identifying some lineages derives from their wide phenotypic and geographic plasticity and the presence of mixed lineages (due to introgressive hybridization with domestic AT populations). In Italy, the S. trutta complex populations living in the Tyrrhenian area and on the main islands (Sicily, Sardinia and Corsica) showed high genetic diversity. Currently, on the Italian Red List, the protected (near threatened) populations are the AD and ME lineages. Recent studies based on traditional (mitochondrial and nuclear markers) and NGS (next-generation sequencing) analyses have clarified some genetic differences between the populations of the Tyrrhenian region, Sicily, Sardinia and Corsica. Native populations in Sardinia belong to the AD lineage, while those living in Corsica are mainly characterized by the AD, MA and ME haplotypes. In Sicily, in the area of the Iblei mountains, an AT lineage (North African) exists. According to some authors, the term Salmo macrostigma should only be used for populations in North Africa. The use of genotyping methods based on mtDNA and nuclear markers and the latest generation sequencing techniques can improve the study of populations and evolutionary lineages in areas where there are overlaps and hybridization phenomena.

Genome-Wide Reconstruction of Rediploidization Following Autopolyploidization across One Hundred Million Years of Salmonid Evolution

The long-term evolutionary impacts of whole-genome duplication (WGD) are strongly influenced by the ensuing rediploidization process. Following autopolyploidization, rediploidization involves a transition from tetraploid to diploid meiotic pairing, allowing duplicated genes (ohnologs) to diverge genetically and functionally. Our understanding of autopolyploid rediploidization has been informed by a WGD event ancestral to salmonid fishes, where large genomic regions are characterized by temporally delayed rediploidization, allowing lineage-specific ohnolog sequence divergence in the major salmonid clades. Here, we investigate the long-term outcomes of autopolyploid rediploidization at genome-wide resolution, exploiting a recent "explosion" of salmonid genome assemblies, including a new genome sequence for the huchen (Hucho hucho). We developed a genome alignment approach to capture duplicated regions across multiple species, allowing us to create 121,864 phylogenetic trees describing genome-wide ohnolog divergence across salmonid evolution. Using molecular clock analysis, we show that 61% of the ancestral salmonid genome experienced an initial "wave" of rediploidization in the late Cretaceous (85-106 Ma). This was followed by a period of relative genomic stasis lasting 17-39 My, where much of the genome remained tetraploid. A second rediploidization wave began in the early Eocene and proceeded alongside species diversification, generating predictable patterns of lineage-specific ohnolog divergence, scaling in complexity with the number of speciation events. Using gene set enrichment, gene expression, and codon-based selection analyses, we provide insights into potential functional outcomes of delayed rediploidization. This study enhances our understanding of delayed autopolyploid rediploidization and has broad implications for future studies of WGD events.

Phylogenomics and loci dropout patterns of deeply diverged Zodarion ant-eating spiders suggest a high potential of RAD-seq for genus-level spider phylogenetics

RAD sequencing yields large amounts of genome-wide data at a relatively low cost and without requiring previous taxon-specific information, making it ideal for evolutionary studies of highly diversified and neglected organisms. However, concerns about information decay with phylogenetic distance have discouraged its use for assessing supraspecific relationships. Here, using Double Digest Restriction Associated DNA (ddRAD) data, we perform the first deep-level approach to the phylogeny of Zodarion, a highly diversified spider genus. We explore the impact of loci and taxon filtering across concatenated and multispecies coalescent reconstruction methods and investigate the patterns of information dropout in reference to both the time of divergence and the mitochondrial divergence between taxa. We found that relaxed loci-filtering and nested taxon-filtering strategies maximized the amount of molecular information and improved phylogenetic inference. As expected, there was a clear pattern of allele dropout towards deeper time and mitochondrial divergences, but the phylogenetic signal remained strong throughout the phylogeny. Therefore, we inferred topologies that were almost fully resolved, highly supported, and noticeably congruent between setups and inference methods, which highlights overall inconsistency in the taxonomy of Zodarion. Because Zodarion appears to be among the oldest and most mitochondrially diversified spider genera, our results suggest that ddRAD data show high potential for inferring intra-generic relationships across spiders and probably also in other taxonomic groups.

Rapid adaptive radiation in a hillstream cyprinid fish in the East African White Nile River basin

Adaptive radiation of freshwater fishes was long thought to be possible only in lacustrine environments. Recently, several studies have shown that riverine and stream environments also provide the ecological opportunity for adaptive radiation. In this study, we report on a riverine adaptive radiation of six ecomorphs of cyprinid hillstream fishes of the genus Garra in a river located in the Ethiopian Highlands in East Africa. Garra are predominantly highly specialized algae-scrapers with a wide distribution ranging from Southeast Asia to West Africa. However, adaptive phenotypic diversification in mouth type, sucking disc morphology, gut length and body shape have probably been found among these ecomorphs in a single Ethiopian river. Moreover, we found two novel phenotypes of Garra ("thick-lipped" and "predatory") that had not been discovered before in this species-rich genus (>160 species). Mitochondrial and genome-wide data suggest monophyletic, intrabasin evolution of Garra phenotypic diversity with signatures of gene flow from other local populations. Although sympatric ecomorphs are genetically distinct and can be considered to being young species as suggested by genome-wide single nucleotide polymorphism data, mitochondrial DNA was unable to identify any genetic structure suggesting recent and rapid speciation events. Some data suggest a hybrid origin of the novel "thick-lipped" ecomorph. Here we highlight how, driven by ecological opportunity, an ancestral trophically highly specialized lineage is likely to have rapidly radiated in a riverine environment promoted by the evolution of novel feeding strategies.

Responses of microRNA and predicted mRNA and enzymatic targets in liver of two salmonids (Oncorhynchus mykiss and Salvelinus fontinalis) following air exposure

The acute stress response is well-characterized, with rainbow trout as a teleost model for physiological and molecular responses. Air exposure, which stimulates an acute stress response, modulates liver microRNAs in rainbow trout; however, these highly conserved non-coding RNAs that bind to mRNA and repress translation, have never been measured in brook trout and it is unknown how miRNA expression responds following air exposure in this less studied salmonid. Our objective was to characterize the effects of air exposure on rainbow and brook trout liver miRNA expression, as well as the mRNA expression and enzyme activity that the miRNAs are predicted to target. Brook and rainbow trout were sampled pre- and 1-, 3-, and 24-h post- a three-minute air exposure. Plasma cortisol, glucose, and lactate were measured. Relative expression of miR-21a-5p, miR-143-3p, let-7a-5p and relative expression and enzyme activities of five predicted targets (pyruvate kinase, glucokinase, citrate synthase, cytochrome c oxidase, and catalase) were measured in liver. Rainbow and brook trout both had increases in plasma cortisol and lactate, while only rainbow trout had significant post-stress increases in plasma glucose. Furthermore, both trout species had increased miR-143-3p and miR-21a-5p relative expression 24-h post-stress. Four of the five enzymes measured had altered activity following stress. Brook trout miRNAs had inverse relative expression with relative catalase mRNA expression and cytochrome c oxidase enzyme activity, but no relationship was found in rainbow trout. Therefore, we have further characterized the transcriptional and enzymatic response to air exposure in two salmonids.

Some Examples of the Use of Molecular Markers for Needs of Basic Biology and Modern Society

Application of molecular genetic markers appeared to be very fruitful in achieving many goals, including (i) proving the theoretic basements of general biology and (ii) assessment of worldwide biodiversity. Both are provided in the present meta-analysis and a review as the main signal. One of the basic current challenges in modern biology in the face of new demands in the 21st century is the validation of its paradigms such as the synthetic theory of evolution (STE) and biological species concept (BSC). Another of most valuable goals is the biodiversity assessment for a variety of social needs including free web-based information resources about any living being, renovation of museum collections, nature conservation that recognized as a global project, iBOL, as well as resolving global trading problems such as false labeling of species specimens used as food, drug components, entertainment, etc. The main issues of the review are focused on animals and combine four items. (1) A combination of nDNA and mtDNA markers best suits the identification of hybrids and estimation of genetic introgression. (2) The available facts on nDNA and mtDNA diversity seemingly make introgression among many taxa obvious, although it is evident, that introgression may be quite restricted or asymmetric, thus, leaving at least the "source" taxon (taxa) intact. (3) If we consider sexually reproducing species in marine and terrestrial realms introgressed, as it is still evident in many cases, then we should recognize that the BSC, in view of the complete lack of gene flow among species, is inadequate because many zoological species are not biological ones yet. However, vast modern molecular data have proven that sooner or later they definitely become biological species. (4) An investigation into the fish taxa divergence using the BOLD database shows that most gene trees are basically monophyletic and interspecies reticulations are quite rare.

A genomic perspective on an old question: Salmo trouts or Salmo trutta (Teleostei: Salmonidae)?

There are particular challenges in defining the taxonomic status of recently radiated groups due to the low level of phylogenetic signal. Members of the Salmo trutta species-complex, which mostly evolved during and following the Pleistocene, show high morphological and ecological diversity that, along with their very wide geographic distribution, have led to morphological description of 47 extant nominal species. However, many of these species have not been supported by previous phylogenetic studies, which could be partly due to lack of significant genetic differences among them, the limited resolution offered by molecular methods previously used, as well as the often local scale of these studies. The development of next-generation sequencing (NGS) and related analytical tools have enhanced our ability to address such challenging questions. In this study, Genotyping-by-Sequencing (GBS) of 15,169 filtered SNPs and mitochondrial DNA (mtDNA) D-loop sequences were combined to assess the phylogenetic relationships among 166 brown trouts representing 21 described species and three undescribed groups collected from 84 localities throughout their natural distribution in Europe, west Asia, and North Africa. The data were analysed using different clustering algorithms (admixture analysis and discriminant analysis of principal components-DAPC), a Bayes Factor Delimitation (BFD) test, species tree reconstruction, gene flow tests (three- and four-population tests), and Rogue taxa identification tests. Genomic contributions of the Atlantic lineage brown trout were found in all major sea basins excluding the North African and Aral Sea basins, suggesting introgressive hybridization of native brown trouts driven by stocking using strains of the Atlantic lineage. After removing the phylogenetic noise caused by the Atlantic brown trout, admixture clusters and DAPC clustering based on GBS data, respectively, resolved 11 and 13 clusters among the previously described brown trout species, which were also supported by BFD test results. Our results suggest that natural hybridization between different brown trout lineages has probably played an important role in the origin of several of the putative species, including S. marmoratus, S. carpio, S. farioides, S. pellegrini, S. caspius (in the Kura River drainage) and Salmo sp. in the Danube River basin. Overall, our results support a multi-species taxonomy for brown trouts. They also resolve some species in the Adriatic-Mediterranean and Black Sea drainages as members of very closely related genomic clusters that may need taxonomic revision. However, any final conclusions pertaining to the taxonomy of the brown trout complex should be based on an integrative approach combining genomic, morphological, and ecological data. To avoid challenges in taxonomy and conservation of species complexes like brown trouts, it is suggested to describe species based on genomic clusters of populations instead of describing species based only on morphologically differentiated single type populations.

Global systematic diversity, range distributions, conservation and taxonomic assessments of graylings (Teleostei: Salmonidae; Thymallus spp.)

Graylings (Thymallus) are among the less well-studied groups of salmonid fishes, especially across their Asian distribution range. Here we perform a comprehensive global review of their phylogeography, systematic diversity and range distributions, including biogeographic reconstruction and assessment of both conservation and taxonomic status of each species. Based on a mitogenomic phylogenetic analysis, three approaches to the delineation of molecular operational units, and evaluation of 15 a-priori defined species, we provide biological support for the recognition of 13 grayling species, plus two additional species tentatively. Several instances of paraphyly and its potential effect on systematic inferences are discussed. Overall, the genus displays increasing species diversity and decreasing range size from higher to lower latitudes and ancestral trait reconstruction supports an East Asian origin for extant diversity, most likely centred in the Amur River drainage. Europe’s colonization by Thymallus took place as early as the late Miocene, at least two colonisations of North America are supported, and multiple dispersal events likely took place into Western Siberia. The conservation status for the 15 taxa was estimated to be: 6 least concern, 1 near-threatened, 2 vulnerable, 3 endangered and 3 data deficient.

Genome-wide RAD sequencing resolves the evolutionary history of serrate leaf Juniperus and reveals discordance with chloroplast phylogeny

Juniper (Juniperus) is an ecologically important conifer genus of the Northern Hemisphere, the members of which are often foundational tree species of arid regions. The serrate leaf margin clade is native to topologically variable regions in North America, where hybridization has likely played a prominent role in their diversification. Here we use a reduced-representation sequencing approach (ddRADseq) to generate a phylogenomic data set for 68 accessions representing all 22 species in the serrate leaf margin clade, as well as a number of close and distant relatives, to improve understanding of diversification in this group. Phylogenetic analyses using three methods (SVDquartets, maximum likelihood, and Bayesian) yielded highly congruent and well-resolved topologies. These phylogenies provided improved resolution relative to past analyses based on Sanger sequencing of nuclear and chloroplast DNA, and were largely consistent with taxonomic expectations based on geography and morphology. Calibration of a Bayesian phylogeny with fossil evidence produced divergence time estimates for the clade consistent with a late Oligocene origin in North America, followed by a period of elevated diversification between 12 and 5 Mya. Comparison of the ddRADseq phylogenies with a phylogeny based on Sanger-sequenced chloroplast DNA revealed five instances of pronounced discordance, illustrating the potential for chloroplast introgression, chloroplast transfer, or incomplete lineage sorting to influence organellar phylogeny. Our results improve understanding of the pattern and tempo of diversification in Juniperus, and highlight the utility of reduced-representation sequencing for resolving phylogenetic relationships in non-model organisms with reticulation and recent divergence.

Addressing incomplete lineage sorting and paralogy in the inference of uncertain salmonid phylogenetic relationships

Recent and continued progress in the scale and sophistication of phylogenetic research has yielded substantial advances in knowledge of the tree of life; however, segments of that tree remain unresolved and continue to produce contradicting or unstable results. These poorly resolved relationships may be the product of methodological shortcomings or of an evolutionary history that did not generate the signal traits needed for its eventual reconstruction. Relationships within the euteleost fish family Salmonidae have proven challenging to resolve in molecular phylogenetics studies in part due to ancestral autopolyploidy contributing to conflicting gene trees. We examine a sequence capture dataset from salmonids and use alternative strategies to accommodate the effects of gene tree conflict based on aspects of salmonid genome history and the multispecies coalescent. We investigate in detail three uncertain relationships: (1) subfamily branching, (2) monophyly of Coregonus and (3) placement of Parahucho. Coregoninae and Thymallinae are resolved as sister taxa, although conflicting topologies are found across analytical strategies. We find inconsistent and generally low support for the monophyly of Coregonus, including in results of analyses with the most extensive dataset and complex model. The most consistent placement of Parahucho is as sister lineage of Salmo.

Evaluating the effect of reference genome divergence on the analysis of empirical RADseq datasets

The advent of high-throughput sequencing (HTS) has made genomic-level analyses feasible for nonmodel organisms. A critical step of many HTS pipelines involves aligning reads to a reference genome to identify variants. Despite recent initiatives, only a fraction of species has publically available reference genomes. Therefore, a common practice is to align reads to the genome of an organism related to the target species; however, this could affect read alignment and bias genotyping. In this study, I conducted an experiment using empirical RADseq datasets generated for two species of salmonids (Actinopterygii; Teleostei; Salmonidae) to address these questions. There are currently reference genomes for six salmonids of varying phylogenetic distance. I aligned the RADseq data to all six genomes and identified variants with several different genotypers, which were then fed into population genetic analyses. Increasing phylogenetic distance between target species and reference genome reduced the proportion of reads that successfully aligned and mapping quality. Reference genome also influenced the number of SNPs that were generated and depth at those SNPs, although the affect varied by genotyper. Inferences of population structure were mixed: increasing reference genome divergence reduced estimates of differentiation but similar patterns of population relationships were found across scenarios. These findings reveal how the choice of reference genome can influence the output of bioinformatic pipelines. It also emphasizes the need to identify best practices and guidelines for the burgeoning field of biodiversity genomics.

Parallelism in eco-morphology and gene expression despite variable evolutionary and genomic backgrounds in a Holarctic fish

Understanding the extent to which ecological divergence is repeatable is essential for predicting responses of biodiversity to environmental change. Here we test the predictability of evolution, from genotype to phenotype, by studying parallel evolution in a salmonid fish, Arctic charr (Salvelinus alpinus), across eleven replicate sympatric ecotype pairs (benthivorous-planktivorous and planktivorous-piscivorous) and two evolutionary lineages. We found considerable variability in eco-morphological divergence, with several traits related to foraging (eye diameter, pectoral fin length) being highly parallel even across lineages. This suggests repeated and predictable adaptation to environment. Consistent with ancestral genetic variation, hundreds of loci were associated with ecotype divergence within lineages of which eight were shared across lineages. This shared genetic variation was maintained despite variation in evolutionary histories, ranging from postglacial divergence in sympatry (ca. 10-15kya) to pre-glacial divergence (ca. 20-40kya) with postglacial secondary contact. Transcriptome-wide gene expression (44,102 genes) was highly parallel across replicates, involved biological processes characteristic of ecotype morphology and physiology, and revealed parallelism at the level of regulatory networks. This expression divergence was not only plastic but in part genetically controlled by parallel cis-eQTL. Lastly, we found that the magnitude of phenotypic divergence was largely correlated with the genetic differentiation and gene expression divergence. In contrast, the direction of phenotypic change was mostly determined by the interplay of adaptive genetic variation, gene expression, and ecosystem size. Ecosystem size further explained variation in putatively adaptive, ecotype-associated genomic patterns within and across lineages, highlighting the role of environmental variation and stochasticity in parallel evolution. Together, our findings demonstrate the parallel evolution of eco-morphology and gene expression within and across evolutionary lineages, which is controlled by the interplay of environmental stochasticity and evolutionary contingencies, largely overcoming variable evolutionary histories and genomic backgrounds.

Maternal Inheritance of U's Triangle and Evolutionary Process of Brassica Mitochondrial Genomes

The sequences and genomic structures of plant mitochondrial (mt) genomes provide unique material for phylogenetic studies. The nature of uniparental inheritance renders an advantage when utilizing mt genomes for determining the parental sources of hybridized taxa. In this study, a concatenated matrix of mt genes was used to infer the phylogenetic relationships of six cultivated Brassica taxa and explore the maternal origins of three allotetraploids. The well-resolved sister relationships between two pairs of diploid and allotetraploid taxa suggest that Brassica carinata (car) possessed a maternal origin from Brassica nigra, while Brassica juncea (jun) was maternally derived from Brassica rapa (cam). Another allotetraploid taxon, Brassica napus (cv. Wester) may have been maternally derived from the common ancestor of B. rapa and Brassica oleracea (ole), and/or have undergone (an) extra hybridization event(s) along its evolutionary history. The characteristics of Brassica mt genomic structures also supported the phylogenetic results. Sinapis arvensis was nested inside the Brassica species, sister to the B. nigra-B. carinata lineage, and possessed an mt genome structure that mostly resembled B. nigra. Collectively, the evidence supported a systematic revision that placed S. arvensis within Brassica. Finally, ancestral mt genomes at each evolutionary node of Brassica were reconstructed, and the detailed and dynamic evolution of Brassica mt genomes was successfully reproduced. The mt genome of B. nigra structurally resembled that of the Brassica ancestor the most, with only one reversion of a block, and the Brassica oleracea underwent the most drastic changes. These findings suggested that repeat-mediated recombinations were largely responsible for the observed structural variations in the evolutionary history of Brassica mt genomes.

Parallel selection on ecologically relevant gene functions in the transcriptomes of highly diversifying salmonids

Background

Salmonid fishes are characterised by a very high level of variation in trophic, ecological, physiological, and life history adaptations. Some salmonid taxa show exceptional potential for fast, within-lake diversification into morphologically and ecologically distinct variants, often in parallel; these are the lake-resident charr and whitefish (several species in the genera Salvelinus and Coregonus). To identify selection on genes and gene categories associated with such predictable diversifications, we analysed 2702 orthogroups (4.82 Mbp total; average 4.77 genes/orthogroup; average 1783 bp/orthogroup). We did so in two charr and two whitefish species and compared to five other salmonid lineages, which do not evolve in such ecologically predictable ways, and one non-salmonid outgroup.

Results

All selection analyses are based on Coregonus and Salvelinus compared to non-diversifying taxa. We found more orthogroups were affected by relaxed selection than intensified selection. Of those, 122 were under significant relaxed selection, with trends of an overrepresentation of serine family amino acid metabolism and transcriptional regulation, and significant enrichment of behaviour-associated gene functions. Seventy-eight orthogroups were under significant intensified selection and were enriched for signalling process and transcriptional regulation gene ontology terms and actin filament and lipid metabolism gene sets. Ninety-two orthogroups were under diversifying/positive selection. These were enriched for signal transduction, transmembrane transport, and pyruvate metabolism gene ontology terms and often contained genes involved in transcriptional regulation and development. Several orthogroups showed signs of multiple types of selection. For example, orthogroups under relaxed and diversifying selection contained genes such as ap1m2, involved in immunity and development, and slc6a8, playing an important role in muscle and brain creatine uptake. Orthogroups under intensified and diversifying selection were also found, such as genes syn3, with a role in neural processes, and ctsk, involved in bone remodelling.

Conclusions

Our approach pinpointed relevant genomic targets by distinguishing among different kinds of selection. We found that relaxed, intensified, and diversifying selection affect orthogroups and gene functions of ecological relevance in salmonids. Because they were found consistently and robustly across charr and whitefish and not other salmonid lineages, we propose these genes have a potential role in the replicated ecological diversifications.

Carotenoid pigmentation in salmon: variation in expression at BCO2-l locus controls a key fitness trait affecting red coloration

Carotenoids are primarily responsible for the characteristic red flesh coloration of salmon. Flesh coloration is an economically and evolutionarily significant trait that varies inter- and intra-specifically, yet the underlying genetic mechanism is unknown. Chinook salmon (Oncorhynchus tshawytscha) represents an ideal system to study carotenoid variation as, unlike other salmonids, they exhibit extreme differences in carotenoid utilization due to genetic polymorphisms. Here, we crossed populations of Chinook salmon with fixed differences in flesh coloration (red versus white) for a genome-wide association study to identify loci associated with pigmentation. Here, the beta-carotene oxygenase 2-like (BCO2-l) gene was significantly associated with flesh colour, with the most significant single nucleotide polymorphism explaining 66% of the variation in colour. BCO2 gene disruption is linked to carotenoid accumulation in other taxa, therefore we hypothesize that an ancestral mutation partially disrupting BCO2-l activity (i.e. hypomorphic mutation) allowed the deposition and accumulation of carotenoids within Salmonidae. Indeed, we found elevated transcript levels of BCO2-l in white Chinook salmon relative to red. The long-standing mystery of why salmon are red, while no other fishes are, is thus probably explained by a hypomorphic mutation in the proto-salmonid at the time of divergence of red-fleshed salmonid genera (approx. 30 Ma).

Have Niche, Will Travel. New Means of Linking Diet and Ecomorphology Reveals Niche Conservatism in Freshwater Cottoid Fishes

Evolutionary transitions between habitats have been catalysts for some of the most stunning examples of adaptive diversification, with novel niches and new resources providing ecological opportunity for such radiations. In aquatic animals, transitions from saltwater to freshwater habitats are rare, but occur often enough that in the Neotropics for example, marine-derived fishes contribute noticeably to regional ichthyofaunal diversity. Here, we investigate how morphology has evolved in a group of temperate fishes that contain a marine to freshwater transition: the sculpins (Percomorpha; Cottoidea). We devised a novel method for classifying dietary niche and relating functional aspects of prey to their predators. Coupled with functional measurements of the jaw apparatus in cottoids, we explored whether freshwater sculpins have fundamentally changed their niche after invading freshwater (niche lability) or if they retain a niche similar to their marine cousins (niche conservatism). Freshwater sculpins exhibit both phylogeographical and ecological signals of phylogenetic niche conservatism, meaning that regardless of habitat, sculpins fill similar niche roles in either saltwater or freshwater. Rather than competition guiding niche conservatism in freshwater cottoids, we argue that strong intrinsic constraints on morphological and ecological evolution are at play, contra to other studies of diversification in marine-derived freshwater fishes. However, several intertidal and subtidal sculpins as well as several pelagic freshwater species from Lake Baikal show remarkable departures from the typical sculpin bauplan. Our method of prey categorization provides an explicit, quantitative means of classifying dietary niche for macroevolutionary studies, rather than relying on somewhat arbitrary means used in previous literature.

A hybrid swarm of manatees along the Guianas coastline, a peculiar environment under the influence of the Amazon River plume

The West Indian (Trichechus manatus) and Amazonian (T. inunguis) manatees have a sympatric occurrence at the mouth of the Amazon River. A result of this interspecific encounter is the occurrence of hybrids, which are frequently found along the coasts of Amapá state in Brazil, French Guiana and Guyana. Here we present new genetic evidence indicating the occurrence of a hybrid swarm along the Guianas Shield coastline, which is an interspecific hybrid zone that also separates T. manatus populations located east (Brazil) and west (Caribbean, Gulf of Mexico, Florida and Antilles). In addition, we suggest that this hybrid population occupies a peculiar mangrove-rich environment under strong influence of the Amazon River plume, which requires an independent management and should be considered a special conservation area.

Riverscape genetic variation, migration patterns, and morphological variation of the threatened Round Rocksnail, Leptoxis ampla

Within riverine systems, headwater populations are hypothesized to harbour higher amounts of genetic distinctiveness than populations in the main stem of a river and display increased genetic diversity in large, downstream habitats. However, these hypotheses were mostly developed with insects and fish, and they have not been tested on many invertebrate lineages. Pleuroceridae gastropods are of particular ecological importance to rivers of eastern North America, sometimes comprising over 90% of macroinvertebrate biomass. Yet, virtually nothing is known of pleurocerid landscape genetics, including whether genetic diversity follows predictions made by hypotheses developed on more mobile species. Moreover, the commonly repeated hypothesis that intraspecific morphological variation in gastropods results from ecophenotypic plasticity has not been well tested on pleurocerids. Using 2bRAD-seq to discover single nucleotide polymorphisms, we show that the threatened, Cahaba River endemic pleurocerid, Leptoxis ampla, has limited gene flow among populations and that migration is downstream-biased, conflicting with previous hypotheses. Both tributary and main stem populations harbour unique genomic profiles, and genetic diversity was highest in downstream populations. Furthermore, L. ampla shell morphology was more correlated with genetic differences among individuals and populations than habitat characteristics. We anticipate similar genetic and demographic patterns to be seen in other pleurocerids, and hypotheses about gene flow and population demographics that were based on more mobile taxa often, but not always, apply to freshwater gastropods. From a conservation standpoint, genetic structure of L. ampla populations suggests distinctive genetic diversity is lost with localized extirpation, a phenomenon common across the range of Pleuroceridae.

Atlantic salmon (Salmo salar L.) genetics in the 21st century: taking leaps forward in aquaculture and biological understanding

Atlantic salmon (Salmo salar L.) is among the most iconic and economically important fish species and was the first member of Salmonidae to have a high-quality reference genome assembly published. Advances in genomics have become increasingly central to the genetic improvement of farmed Atlantic salmon as well as conservation of wild salmon stocks. The salmon genome has also been pivotal in shaping our understanding of the evolutionary and functional consequences arising from an ancestral whole-genome duplication event characterising all Salmonidae members. Here, we provide a review of the current status of Atlantic salmon genetics and genomics, focussed on progress made from genome-wide research aimed at improving aquaculture production and enhancing understanding of salmonid ecology, physiology and evolution. We present our views on the future direction of salmon genomics, including the role of emerging technologies (e.g. genome editing) in elucidating genetic features that underpin functional variation in traits of commercial and evolutionary importance.

Phylogeny and species delimitation of near Eastern Neurergus newts (Salamandridae) based on genome-wide RADseq data analysis

We reconstruct the molecular phylogeny of Near Eastern mountain brook newts of the genus Neurergus (family Salamandridae) based on newly determined RADseq data, and compare the outcomes of concatenation-based phylogenetic reconstruction with species-tree inference. Furthermore, we test the current taxonomy of Neurergus (with four species: Neurergus strauchii, N. crocatus, N. kaiseri, and N. derjugini) against coalescent-based species-delimitation approaches of our genome-wide genetic data set. While the position of N. strauchii as sister species to all other Neurergus species was consistent in all of our analyses, the phylogenetic relationships between the three remaining species changed depending on the applied method. The concatenation approach, as well as quartet-based species-tree inference, supported a topology with N. kaiseri as the closest relative to N. derjugini, while full-coalescent species-tree inference approaches supported N. crocatus as sister species of N. derjugini. Investigating the individual signal of gene trees highlighted an extensive variation among gene histories, most likely resulting from incomplete lineage sorting. Coalescent-based species-delimitation models suggest that the current taxonomy might underestimate the species richness within Neurergus and supports seven species. Based on the current sampling, our analysis suggests that N. strauchii, N. derjugini and N. kaiseri might each be subdivided into further species. However, as amphibian species are known to be composed of deep conspecific lineages that do not always warrant species status, these results need to be cautiously interpreted in an integrative taxonomic framework. We hypothesize that the rather shallow divergences detected within N. kaiseri and N. derjugini likely reflect an ongoing speciation process and thus require further investigation. On the contrary, the much deeper genetic divergence found between the two morphologically and geographically differentiated subspecies of N. strauchii leads us to propose that N. s. barani should be considered a distinct species, Neurergus barani Öz, 1994.

Small and genetically highly structured populations in a long-legged bee, Rediviva longimanus, as inferred by pooled RAD-seq

Adaptation to local host plants may impact a pollinator's population genetic structure by reducing gene flow and driving population genetic differentiation, representing an early stage of ecological speciation. South African Rediviva longimanus bees exhibit elongated forelegs, a bizarre adaptation for collecting oil from floral spurs of their Diascia hosts. Furthermore, R. longimanus foreleg length (FLL) differs significantly among populations, which has been hypothesised to result from selection imposed by inter-population variation in Diascia floral spur length. Here, we used a pooled restriction site-associated DNA sequencing (pooled RAD-seq) approach to investigate the population genetic structure of R. longimanus and to test if phenotypic differences in FLL translate into increased genetic differentiation (i) between R. longimanus populations and (ii) between phenotypes across populations. We also inferred the effects of demographic processes on population genetic structure and tested for genetic markers underpinning local adaptation. RESULTS: Populations showed marked genetic differentiation (average FST = 0.165), though differentiation was not statistically associated with differences between populations in FLL. All populations exhibited very low genetic diversity and were inferred to have gone through recent bottleneck events, suggesting extremely low effective population sizes. Genetic differentiation between samples pooled by leg length (short versus long) rather than by population of origin was even higher (FST = 0.260) than between populations, suggesting reduced interbreeding between long and short-legged individuals. Signatures of selection were detected in 1119 (3.8%) of a total of 29,721 SNP markers, CONCLUSIONS: Populations of R. longimanus appear to be small, bottlenecked and isolated. Though we could not detect the effect of local adaptation (FLL in response to floral spurs of host plants) on population genetic differentiation, short and long legged bees appeared to be partially differentiated, suggesting incipient ecological speciation. To test this hypothesis, greater resolution through the use of individual-based whole-genome analyses is now needed to quantify the degree of reproductive isolation between long and short legged bees between and even within populations.