Morphometric and population genomic evidence for species divergence in the Chimarrichthys fish complex of the Tibetan Plateau

Abstract

The uplift of the Tibetan Plateau altered the environmental conditions of the local area substantially. Here, we conducted a comprehensive investigation based on morphometrics, population genomics, and climatic factors to evaluate phenotypic and genome-level variations in a radiation of Chimarrichthys catfish endemic to the Plateau. Discriminant function analysis showed phenotypic differences of Chimarrichthys between rivers with respect to elevation. Genetic structure analysis based on 6606 single nucleotide polymorphisms (SNPs) deduced genetic differences between rivers, and species delimitation indicated that the Chimarrichthys fish complex could be divided into three species. Restriction site-associated DNA tags were mapped to the gene sets of Glyptosternon maculatum, and matches were searched against databases for Gene Ontology annotation. Genomic regions exhibiting marked differences among localities represented a range of biological functions, including growth (gdf11), bone development (bmp8a), cellular response to light stimulus (opn3), regulation of the rhodopsin-mediated signalling pathway (grk1), immune response (rag1 and ung), reproductive process (antxr2), and regulation of intracellular iron levels (ireb2). The tag44126, where gene gdf11 is located, was identified as an outlier exhibiting divergence between rivers with altitude differences, and the SNP is thymine (T) in Dadu and Yalong River (~2700 m), but guanine (G) in Jinsha and Qingyi rivers (~2200 and ~ 684 m), suggesting a possible effect of altitude on its differentiation.