Untangling the evolutionary history of a highly polymorphic species: introgressive hybridization and high genetic structure in the desert cichlid fish Herichtys minckleyi

Genetic and ecomorphological divergence between sympatric Astyanax morphs from Central America

Intraspecific ecological and morphological polymorphism can promote ecological speciation and the build-up of reproductive isolation. Here, we evaluate correlations among morphology, trophic ecology and genetic differentiation between two divergent morphs (elongate and deep-body) of the fish genus Astyanax in the San Juan River basin in Central America, to infer the putative evolutionary mechanism shaping this system. We collected the two morphs from three water bodies and analysed: (1) the correlation between body shape and the shape of the premaxilla, a relevant trophic morphological structure, (2) the trophic level and niche width of each morph, (3) the correspondence between trophic level and body and premaxillary shape, and (4) the genetic differentiation between morphs using mitochondrial and nuclear markers. We found a strong correlation between the body and premaxillary shape of the morphs. The elongate-body morph had a streamlined body, a premaxilla with acuter angles and a narrower ascending process, and a higher trophic level, characteristic of species with predatorial habits. By contrast, the deep-body morph had a higher body depth, a premaxilla with less acute angles and a broader trophic niche, suggesting generalist habits. Despite the strong correlation between morphological and ecological divergence, the morphs showed limited genetic differentiation, supporting the idea that morphs may be undergoing incipient ecological speciation, although alternative scenarios such as stable polymorphism or plasticity should also be considered. This study provides support for the role of ecological factors promoting diversification in both lake and stream-dwelling freshwater fish.

Pervasive admixture and the spread of a large-lipped form in a cichlid fish radiation

Adaptive radiations have proven important for understanding the mechanisms and processes underlying biological diversity. The convergence of form and function, as well as admixture and adaptive introgression, are common in adaptive radiations. However, distinguishing between these two scenarios remains a challenge for evolutionary research. The Midas cichlid species complex (Amphilophus spp.) is a prime example of adaptive radiation, with phenotypic diversification occurring at various stages of genetic differentiation. One species, A. labiatus, has large fleshy lips, is associated with rocky lake substrates, and occurs patchily within Lakes Nicaragua and Managua. By contrast, the similar, but thin-lipped, congener, A. citrinellus, is more common and widespread. We investigated the evolutionary history of the large-lipped form, specifically regarding whether the trait has evolved independently in both lakes from ancestral thin-lipped populations, or via dispersal and/or admixture events. We collected samples from distinct locations in both lakes, and assessed differences in morphology and ecology. Using RAD-seq, we genotyped thousands of SNPs to measure population structure and divergence, demographic history, and admixture. We found significant between-species differences in ecology and morphology, local intraspecific differences in body shape and trophic traits, but only limited intraspecific variation in lip shape. Despite clear ecological differences, our genomic approach uncovered pervasive admixture between the species and low genomic differentiation, with species within lakes being genetically more similar than species between lakes. Taken together, our results suggest a single origin of large-lips, followed by pervasive admixture and adaptive introgression, with morphology being driven by local ecological opportunities, despite ongoing gene-flow.

Species specificity and sexual dimorphism in tooth shape among the three sympatric haplochromine species in Lake Kivu cichlids

Tooth shape is used to differentiate between morphologically similar species of vertebrates, including fish. This study aimed to quantify tooth shape of three sympatric species: Haplochromis kamiranzovu, H. insidiae, and H. astatodon endemic to Lake Kivu, whose existing identification criteria are currently only qualitative. A quantitative tooth shape analysis was performed based on digitized tooth outline data with a subsequent elliptic Fourier analysis to test for differences among the three species. We looked at crown shape and size differences within H. kamiranzovu and H. insidiae at geographical, habitat, and gender levels. No comparison at habitat level was done for H. astatodon because it is found only in littoral zone. The analysis revealed significant tooth shape differences among the three species. Haplochromis astatodon had a significantly longer major cusp height and a longer and larger minor cusp than that of H. insidiae. It had also a longer major cusp height and a longer and larger minor cusp than that of H. kamiranzovu. Tooth shape differences of H. kamiranzovu and H. insidiae species were not significantly different between littoral and pelagic fish (p > .05) while differences were significant between southern and northern Lake Kivu populations (p < .05). Tooth sizes in H. kamiranzovu and H. insidiae were significantly different, both in height and width as well as in their ratios, and this was true at sex and geographic levels (p < .05), but not at habitat level (p > .05). Tooth shape was also significantly different with sharp teeth for males compared with females of southern populations versus northern ones. These shape- and size-related differences between sexes suggest differences in the foraging strategies toward available food resources in the lake habitat. Further research should explain the genetic basis of the observed pattern.

Population Genomic Evidence Reveals Subtle Patterns of Differentiation in the Trophically Polymorphic Cuatro Ciénegas Cichlid, Herichthys minckleyi

In recent decades, an increased understanding of molecular ecology has led to a reinterpretation of the role of gene flow during the evolution of reproductive isolation and biological novelty. For example, even in the face of ongoing gene flow strong selection may maintain divergent polymorphisms, or gene flow may introduce novel biological diversity via hybridization and introgression from a divergent species. Herein, we elucidate the evolutionary history and genomic basis of a trophically polymorphic trait in a species of cichlid fish, Herichthys minckleyi. We explored genetic variation at 3 hierarchical levels; between H. minckleyi (n = 69) and a closely related species Herichthys cyanoguttatus (n = 10), between H. minckleyi individuals from 2 geographic locations, and finally between individuals with alternate morphotypes at both a genome-wide and locus-specific scale. We found limited support for the hypothesis that the H. minckleyi polymorphism is the result of ongoing hybridization between the 2 species. Within H. minckleyi we found evidence of geographic genetic structure, and using traditional population genetic analyses found that individuals of alternate morphotypes within a pool appear to be panmictic. However, when we used a locus-specific approach to examine the relationship between multi-locus genotype, tooth size, and geographic sampling, we found the first evidence for molecular genetic differences between the H. minckleyi morphotypes.

Trophic specialization and morphological divergence between two sympatric species in Lake Catemaco, Mexico

The association of morphological divergence with ecological segregation among closely related species could be considered as a signal of divergent selection in ecological speciation processes. Environmental signals such as diet can trigger phenotypic evolution, making polymorphic species valuable systems for studying the evolution of trophic-related traits. The main goal of this study was to analyze the association between morphological differences in trophic-related traits and ecological divergence in two sympatric species, Astyanax aeneus and A. caballeroi, inhabiting Lake Catemaco, Mexico. The trophic differences of a total of 70 individuals (35 A. aeneus and 35 A. caballeroi) were examined using stable isotopes and gut content analysis; a subset of the sample was used to characterize six trophic and six ecomorphological variables. In our results, we recovered significant differences between both species in the values of stable isotopes, with higher values of δ15N for A. caballeroi than for A. aeneus. Gut content results were consistent with the stable isotope data, with a higher proportion of invertebrates in A. caballeroi (a consumption of invertebrates ten times higher than that of A. aeneus, which in turn consumed three times more vegetal material than A. caballeroi). Finally, we found significant relationship between ecomorphology and stable isotopes (r = .24, p < .01), hence, head length, preorbital length, eye diameter, and δ15N were all positively correlated; these characteristics correspond to A. caballeroi. While longer gut and gill rakers, deeper bodies, and vegetal material consumption were positively correlated and corresponded to A. aeneus. Our results are consistent with the hypothesis that morphological divergence in trophic-related traits could be associated with niche partitioning, allowing the coexistence of closely related species and reducing interspecific competition.

Interspecies introgressive hybridization in spiny frogs Quasipaa (Family Dicroglossidae) revealed by analyses on multiple mitochondrial and nuclear genes

Introgression may lead to discordant patterns of variation among loci and traits. For example, previous phylogeographic studies on the genus Quasipaa detected signs of genetic introgression from genetically and morphologically divergent Quasipaa shini or Quasipaa spinosa. In this study, we used mitochondrial and nuclear DNA sequence data to verify the widespread introgressive hybridization in the closely related species of the genus Quasipaa, evaluate the level of genetic diversity, and reveal the formation mechanism of introgressive hybridization. In Longsheng, Guangxi Province, signs of asymmetrical nuclear introgression were detected between Quasipaa boulengeri and Q. shini. Unidirectional mitochondrial introgression was revealed from Q. spinosa to Q. shini. By contrast, bidirectional mitochondrial gene introgression was detected between Q. spinosa and Q. shini in Lushan, Jiangxi Province. Our study also detected ancient hybridizations between a female Q. spinosa and a male Q. jiulongensis in Zhejiang Province. Analyses on mitochondrial and nuclear genes verified three candidate cryptic species in Q. spinosa, and a cryptic species may also exist in Q. boulengeri. However, no evidence of introgressive hybridization was found between Q. spinosa and Q. boulengeri. Quasipaa exilispinosa from all the sampling localities appeared to be deeply divergent from other communities. Our results suggest widespread introgressive hybridization in closely related species of Quasipaa and provide a fundamental basis for illumination of the forming mechanism of introgressive hybridization, classification of species, and biodiversity assessment in Quasipaa.

Phenotypes in phylogeography: Species' traits, environmental variation, and vertebrate diversification

Almost 30 y ago, the field of intraspecific phylogeography laid the foundation for spatially explicit and genealogically informed studies of population divergence. With new methods and markers, the focus in phylogeography shifted to previously unrecognized geographic genetic variation, thus reducing the attention paid to phenotypic variation in those same diverging lineages. Although phenotypic differences among lineages once provided the main data for studies of evolutionary change, the mechanisms shaping phenotypic differentiation and their integration with intraspecific genetic structure have been underexplored in phylogeographic studies. However, phenotypes are targets of selection and play important roles in species performance, recognition, and diversification. Here, we focus on three questions. First, how can phenotypes elucidate mechanisms underlying concordant or idiosyncratic responses of vertebrate species evolving in shared landscapes? Second, what mechanisms underlie the concordance or discordance of phenotypic and phylogeographic differentiation? Third, how can phylogeography contribute to our understanding of functional phenotypic evolution? We demonstrate that the integration of phenotypic data extends the reach of phylogeography to explain the origin and maintenance of biodiversity. Finally, we stress the importance of natural history collections as sources of high-quality phenotypic data that span temporal and spatial axes.

Functional modularity in lake-dwelling characin fishes of Mexico

Modular evolution promotes evolutionary change, allowing independent variation across morphological units. Recent studies have shown that under contrasting ecological pressures, patterns of modularity could be related to divergent evolution. The main goal of the present study was to evaluate the presence of modular evolution in two sister lacustrine species, Astyanax aeneus and A. caballeroi, which are differentiated by their trophic habits. Two different datasets were analyzed: (1) skull X-rays from 73 specimens (35 A. aeneus and 38 A. caballeroi) to characterize skull variation patterns, considering both species and sex effects. For this dataset, three different modularity hypotheses were tested, previously supported in other lacustrine divergent species; (2) a complete body shape dataset was also tested for four modularity hypotheses, which included a total of 196 individuals (110 Astyanax aeneus and 86 A. caballeroi). Skull shape showed significant differences among species and sex (P < 0.001), where Astyanax caballeroi species showed an upwardly projected mandible and larger preorbital region. For the skull dataset, the modularity hypothesis ranked first included three partitioning modules. While for the complete body dataset the best ranked hypothesis included two modules (head vs the rest of the body), being significant only for A. caballeroi.

New insight into hybridization and unidirectional introgression between Ammodytes japonicus and Ammodytes heian (Trachiniformes, Ammodytidae)

Based on northern (NOL) and southern (SOL) mitochondrial lineages, recently, it proposed the new species Ammodytes heian and revived the species name Ammodytes japonicus to describe sand lances from the northwestern Pacific Ocean. This study used molecular methods to investigate genetic relationships between the two sand lance species in Korea and Japan. In total, 154 specimens were collected from four locations in Korea (Baengnyeongdo in the Yellow Sea, Tongyeong in the Korean Strait, and Jumunjin and Gijang in the East Sea), and 50 specimens were collected from a single location in Japan (Wakkanai in the Okhotsk Sea). Mitochondrial DNA analysis demonstrated that the individuals from Baengnyeongdo and Tongyeong all belonged to the SOL, whereas those from Gijang, Jumunjin, and Wakkanai included individuals from both the NOL and SOL (over 75% NOL). Population structure analyses were performed on the same individuals using seven microsatellite DNA markers. The population structure analysis based on 201 specimens identified two clusters (named as northern group and southern group), with the admixture proportion (q) of < 0.1 for the northern group in the Backyeongdo and Tongyeong sand lances and < 0.1 for the southern group in the Wakkanai sand lances. The high heterogeneity indicated that the former was probably A. japonicus and the latter probably A. heian. However, the admixture proportion in the Jumunjin and Gijang sand lances was 0.71–0.75 for the southern group, indicating that hybridization and unidirectional introgression from SOL to NOL occurs in southwestern margin of the East Sea. Our findings illustrate the speciation process based on different patterns of gene flow between Korean and Japanese sand lance, which is strongly influenced by both the paleo-climatic change and the contemporary local oceanic current pattern.

Molecular ecology studies of species radiations: current research gaps, opportunities and challenges

Understanding the drivers and limits of species radiations is a crucial goal of evolutionary genetics and molecular ecology, yet research on this topic has been hampered by the notorious difficulty of connecting micro- and macroevolutionary approaches to studying the drivers of diversification. To chart the current research gaps, opportunities and challenges of molecular ecology approaches to studying radiations, we examine the literature in the journal Molecular Ecology and revisit recent high-profile examples of evolutionary genomic research on radiations. We find that available studies of radiations are highly unevenly distributed among taxa, with many ecologically important and species-rich organismal groups remaining severely understudied, including arthropods, plants and fungi. Most studies employed molecular methods suitable over either short or long evolutionary time scales, such as microsatellites or restriction site-associated DNA sequencing (RAD-seq) in the former case and conventional amplicon sequencing of organellar DNA in the latter. The potential of molecular ecology studies to address and resolve patterns and processes around the species level in radiating groups of taxa is currently limited primarily by sample size and a dearth of information on radiating nuclear genomes as opposed to organellar ones. Based on our literature survey and personal experience, we suggest possible ways forward in the coming years. We touch on the potential and current limitations of whole-genome sequencing (WGS) in studies of radiations. We suggest that WGS and targeted ('capture') resequencing emerge as the methods of choice for scaling up the sampling of populations, species and genomes, including currently understudied organismal groups and the genes or regulatory elements expected to matter most to species radiations.

Do relaxed selection and habitat temperature facilitate biased mitogenomic introgression in a narrowly endemic fish?

Introgression might be exceptionally common during the evolution of narrowly endemic species. For instance, in the springs of the small and isolated Cuatro Ciénegas Valley, the mitogenome of the cichlid fish Herichthys cyanoguttatus could be rapidly introgressing into populations of the trophically polymorphic H. minckleyi. We used a combination of genetic and environmental data to examine the factors associated with this mitochondrial introgression. A reduced representation library of over 6220 single nucleotide polymorphisms (SNPs) from the nuclear genome showed that mitochondrial introgression into H. minckleyi is biased relative to the amount of nuclear introgression. SNP assignment probabilities also indicated that cichlids with more hybrid ancestry are not more commonly female providing no support for asymmetric backcrossing or hybrid‐induced sex‐ratio distortion in generating the bias in mitochondrial introgression. Smaller effective population size in H. minckleyi inferred from the SNPs coupled with sequences of all 13 mitochondrial proteins suggests that relaxed selection on the mitogenome could be facilitating the introgression of “H. cyanoguttatus” haplotypes. Additionally, we showed that springs with colder temperatures had greater amounts of mitochondrial introgression from H. cyanoguttatus. Relaxed selection in H. minckleyi coupled with temperature‐related molecular adaptation could be facilitating mitogenomic introgression into H. minckleyi.

Historical and current introgression in a Mesoamerican hummingbird species complex: a biogeographic perspective

The influence of geologic and Pleistocene glacial cycles might result in morphological and genetic complex scenarios in the biota of the Mesoamerican region. We tested whether berylline, blue-tailed and steely-blue hummingbirds, Amazilia beryllina, Amazilia cyanura and Amazilia saucerottei, show evidence of historical or current introgression as their plumage colour variation might suggest. We also analysed the role of past and present climatic events in promoting genetic introgression and species diversification. We collected mitochondrial DNA (mtDNA) sequence data and microsatellite loci scores for populations throughout the range of the three Amazilia species, as well as morphological and ecological data. Haplotype network, Bayesian phylogenetic and divergence time inference, historical demography, palaeodistribution modelling, and niche divergence tests were used to reconstruct the evolutionary history of this Amazilia species complex. An isolation-with-migration coalescent model and Bayesian assignment analysis were assessed to determine historical introgression and current genetic admixture. mtDNA haplotypes were geographically unstructured, with haplotypes from disparate areas interdispersed on a shallow tree and an unresolved haplotype network. Assignment analysis of the nuclear genome (nuDNA) supported three genetic groups with signs of genetic admixture, corresponding to: (1) A. beryllina populations located west of the Isthmus of Tehuantepec; (2) A. cyanura populations between the Isthmus of Tehuantepec and the Nicaraguan Depression (Nuclear Central America); and (3) A. saucerottei populations southeast of the Nicaraguan Depression. Gene flow and divergence time estimates, and demographic and palaeodistribution patterns suggest an evolutionary history of introgression mediated by Quaternary climatic fluctuations. High levels of gene flow were indicated by mtDNA and asymmetrical isolation-with-migration, whereas the microsatellite analyses found evidence for three genetic clusters with distributions corresponding to isolation by the Isthmus of Tehuantepec and the Nicaraguan Depression and signs of admixture. Historical levels of migration between genetically distinct groups estimated using microsatellites were higher than contemporary levels of migration. These results support the scenario of secondary contact and range contact during the glacial periods of the Pleistocene and strongly imply that the high levels of structure currently observed are a consequence of the limited dispersal of these hummingbirds across the isthmus and depression barriers.

Response paper: Morphometric article by Mejía et al. 2015 alluding genera Herichthys and Nosferatu displays serious inconsistencies

ABSTRACT In the present response paper, the article entitled "Morphometric variation of the Herichthys bartoni (Bean, 1892) species group (Teleostei: Cichlidae): How many species comprise H. labridens (Pellegrin, 1903)?" by Mejía et al . 2015 is critically reviewed. The current review pinpoints some of the more conspicuous conceptual inconsistencies and fundamental errors found in the study by Mejía et al . (2015), It is contended that the authors fail to provide any new insights into the complex biogeography and evolutionary history of the Nosferatu and Herichthys genus groups, and that while results of their Cox1 molecular analysis are comparable to those by De la Maza-Benignos et al . (2015), the conclusions of the two studies are not comparable. In addition, it is contested that, whereas the designation of genus Nosferatu by De la Maza-Benignos et al . (2015) was found on the principles of the biological and phylogenetic species concepts, the rejection of the genus by Mejía et al . (2015) is solely based "on the presence of (overlapping) morphometric characters" between genera. The assumption by Mejía et al . (2015),that because their geometric morphometrics analysis failed to provide separation of species, then Nosferatu genus does not correspond to a valid taxon; and their suggesting geometric morphometrics "as useful tool to discriminate species, because it allows to propose diagnostic characters" were not supported by their results. While Mejía et al . present some interesting thoughts on the systematics of Nosferatu , they unfortunately fail to provide any data that can be objectively assessed as relevant to motivate any changes in the current taxonomy.