Tracing the genetic footprints of vertebrate landing in non-teleost ray-finned fishes.

Rich fossil evidence suggests that many traits and functions related to terrestrial evolution were present long before the ancestor of lobe- and ray-finned fishes. Here, we present genome sequences of the bichir, paddlefish, bowfin, and alligator gar, covering all major early divergent lineages of ray-finned fishes. Our analyses show that these species exhibit many mosaic genomic features of lobe- and ray-finned fishes. In particular, many regulatory elements for limb development are present in these fishes, supporting the hypothesis that the relevant ancestral regulation networks emerged before the origin of tetrapods. Transcriptome analyses confirm the homology between the lung and swim bladder and reveal the presence of functional lung-related genes in early ray-finned fishes. Furthermore, we functionally validate the essential role of a jawed vertebrate highly conserved element for cardiovascular development. Our results imply the ancestors of jawed vertebrates already had the potential gene networks for cardio-respiratory systems supporting air breathing.

African lungfish genome sheds light on the vertebrate water-to-land transition.

Lungfishes are the closest extant relatives of tetrapods and preserve ancestral traits linked with the water-to-land transition. However, their huge genome sizes have hindered understanding of this key transition in evolution. Here, we report a 40-Gb chromosome-level assembly of the African lungfish (Protopterus annectens) genome, which is the largest genome assembly ever reported and has a contig and chromosome N50 of 1.60 Mb and 2.81 Gb, respectively. The large size of the lungfish genome is due mainly to retrotransposons. Genes with ultra-long length show similar expression levels to other genes, indicating that lungfishes have evolved high transcription efficacy to keep gene expression balanced. Together with transcriptome and experimental data, we identified potential genes and regulatory elements related to such terrestrial adaptation traits as pulmonary surfactant, anxiolytic ability, pentadactyl limbs, and pharyngeal remodeling. Our results provide insights and key resources for understanding the evolutionary pathway leading from fishes to humans.

Multi-omics Investigation of Freeze Tolerance in the Amur Sleeper, an Aquatic Ectothermic Vertebrate.

Freeze tolerance, the ability of an organism to survive internal ice formation, is a striking survival strategy employed by some ectotherms living in cold environments. However, the genetic bases of this remarkable adaptation are largely unknown. The Amur sleeper (Perccottus glenii), the only known freeze-tolerant fish species, can overwinter with its entire body frozen in ice. Here, we sequenced the chromosome-level genome of the Amur sleeper and performed comparative genomic, transcriptomic, and metabolomic analyses to investigate its strategies for surviving freezing. Evolutionary analysis suggested that the Amur sleeper diverged from its closest non-cold-hardy relative about 15.07 million years ago and has experienced a high rate of protein evolution. Transcriptomic and metabolomic data identified a coordinated and tissue-specific regulation of genes and metabolites involved in hypometabolism, cellular stress response, and cryoprotectant accumulation involved in freezing and thawing. Several genes show evidence of accelerated protein sequence evolution or family size expansion were found as adaptive responses to freezing-induced stresses. Specifically, genetic changes associated with cytoskeleton stability, cryoprotectant synthesis, transmembrane transport, and neuroprotective adaptations were identified as potentially key innovations that aid in freezing survival. Our work provides valuable resources and opportunities to unveil the molecular adaptations supporting freeze tolerance in ectothermic vertebrates.

The Role of Introgression During the Radiation of Endemic Fishes Adapted to Living at Extreme Altitudes in the Tibetan Plateau.

Recent genomic analyses of evolutionary radiations suggest that ancient introgression may facilitate rapid diversification and adaptive radiation. The loach genus Triplophysa, a genus with most species endemic to Tibetan Plateau, shows ecological diversity and rapid evolution and represents a potential example of adaptive radiation linked to the uplift of the Tibetan Plateau. Here, we interrogate the complex evolutionary history of Triplophysa fishes through the analysis of whole-genome sequences. By reconstructing the phylogeny of Triplophysa, quantifying introgression across this clade, and simulating speciation and migration processes, we confirm that extensive gene flow events occurred across disparate Triplophysa species. Our results suggest that introgression plays a more substantial role than incomplete lineage sorting in underpinning phylogenetic discordance in Triplophysa. The results also indicate that genomic regions affected by ancient gene flow exhibit characteristics of lower recombination rates and nucleotide diversity and may associate with selection. Simulation analysis of Triplophysa tibetana suggests that the species may have been affected by the Gonghe Movement in the third uplift of the Tibetan Plateau, resulting in founder effects and a subsequent reduction in Ne.

Comparative iTRAQ proteomics revealed proteins associated with lobed fin regeneration in Bichirs.

BACKGROUND: Polypterus senegalus can fully regenerate its pectoral lobed fins, including a complex endoskeleton, with remarkable precision. However, despite the enormous potential of this species for use in medical research, its regeneration mechanisms remain largely unknown. METHODS: To identify the differentially expressed proteins (DEPs) during the early stages of lobed fin regeneration in P. senegalus, we performed a differential proteomic analysis using isobaric tag for relative and absolute quantitation (iTRAQ) approach based quantitative proteome from the pectoral lobed fins at 3 time points. Furthermore, we validated the changes in protein expression with multiple-reaction monitoring (MRM) analysis. RESULTS: The experiment yielded a total of 3177 proteins and 15,091 unique peptides including 1006 non-redundant (nr) DEPs. Of these, 592 were upregulated while 349 were downregulated after lobed fin amputation when compared to the original tissue. Bioinformatics analyses showed that the DEPs were mainly associated with Ribosome and RNA transport, metabolic, ECM-receptor interaction, Golgi and endoplasmic reticulum, DNA replication, and Regulation of actin cytoskeleton. CONCLUSIONS: To our knowledge, this is the first proteomic research to investigate alterations in protein levels and affected pathways in bichirs' lobe-fin/limb regeneration. In addition, our study demonstrated a highly dynamic regulation during lobed fin regeneration in P. senegalus. These results not only provide a comprehensive dataset on differentially expressed proteins during the early stages of lobe-fin/limb regeneration but also advance our understanding of the molecular mechanisms underlying lobe-fin/limb regeneration.

Pseudo-chromosome-length genome assembly for a deep-sea eel Ilyophis brunneus sheds light on the deep-sea adaptation.

High hydrostatic pressure, low temperature, and scarce food supply are the major factors that limit the survival of vertebrates in extreme deep-sea environments. Here, we constructed a high-quality genome of the deep-sea Muddy arrowtooth eel (MAE, Ilyophis brunneus, captured below a depth of 3,500 m) by using Illumina, PacBio, and Hi-C sequencing. We compare it against those of shallow-water eel and other outgroups to explore the genetic basis that underlies the adaptive evolution to deep-sea biomes. The MAE genome was estimated to be 1.47 Gb and assembled into 14 pseudo-chromosomes. Phylogenetic analyses indicated that MAE diverged from its closely related shallow-sea species, European eel, ∼111.9 Mya and experienced a rapid evolution. The genome evolutionary analyses primarily revealed the following: (i) under high hydrostatic pressure, the positively selected gene TUBGCP3 and the expanded family MLC1 may improve the cytoskeleton stability; ACOX1 may enhance the fluidity of cell membrane and maintain transport activity; the expansion of ABCC12 gene family may enhance the integrity of DNA; (ii) positively selected HARS likely maintain the transcription ability at low temperatures; and (iii) energy metabolism under a food-limited environment may be increased by expanded and positively selected genes in AMPK and mTOR signaling pathways.

Monsoon boosted radiation of the endemic East Asian carps.

Major historical events often trigger the rapid flourishing of a few lineages, which in turn shape established biodiversity patterns. How did this process occur and develop? This study provides a window into this issue. The endemic East Asian carps (EEAC) dominated the ichthyofauna of East Asia and exhibited a high degree of adaptation to monsoonal river-lake ecosystems. A series of evidence, including ecogeography, phylogenetics, and macroevolution, suggests that the EEAC is a lineage that arose with the East Asian monsoon and thrived intimately with subsequent monsoon activities. We further deduce the evolution of the EEAC and find that a range of historical events in the monsoon setting (e.g., marine transgression and regression and glacial-interglacial cycle) have further reshaped the distribution patterns of EEAC's members. Comparative genomics analyses reveal that introgressions during the initial period of EEAC radiation and innovations in the regulation of the brain and nervous system may have aided their adaptation to river-lake ecosystems in a monsoon setting, which boosted radiation. Overall, this study strengthens knowledge of the evolutionary patterns of freshwater fishes in East Asia and provides a model case for understanding the impact of major historical events on the evolution of biota.

Chromosome-level genome assembly of hadal snailfish reveals mechanisms of deep-sea adaptation in vertebrates.

As the deepest vertebrate in the ocean, the hadal snailfish (Pseudoliparis swirei), which lives at a depth of 6,000-8,000 m, is a representative case for studying adaptation to extreme environments. Despite some preliminary studies on this species in recent years, including their loss of pigmentation, visual and skeletal calcification genes, and the role of trimethylamine N-oxide in adaptation to high-hydrostatic pressure, it is still unknown how they evolved and why they are among the few vertebrate species that have successfully adapted to the deep-sea environment. Using genomic data from different trenches, we found that the hadal snailfish may have entered and fully adapted to such extreme environments only in the last few million years. Meanwhile, phylogenetic relationships show that they spread into different trenches in the Pacific Ocean within a million years. Comparative genomic analysis has also revealed that the genes associated with perception, circadian rhythms, and metabolism have been extensively modified in the hadal snailfish to adapt to its unique environment. More importantly, the tandem duplication of a gene encoding ferritin significantly increased their tolerance to reactive oxygen species, which may be one of the important factors in their adaptation to high-hydrostatic pressure.

Sympatric or micro-allopatric speciation in a glacial lake? Genomic islands support neither.

Apparent cases of sympatric speciation may actually be due to micro-allopatric or micro-parapatric speciation. One way to distinguish between these models is to examine the existence and nature of genomic islands of divergence, wherein divergent DNA segments are interspersed with low-divergence segments. Such islands should be rare or absent under micro-allopatric speciation but common in cases of speciation with gene flow. Sympatric divergence of endemic fishes is known from isolated saline, crater, postglacial, and ancient lakes. Two morphologically distinct cyprinid fishes, Gymnocypris eckloni scoliostomus (GS) and G. eckloni eckloni (GE), in a small glacial lake on the Qinghai-Tibet Plateau, Lake Sunmcuo, match the biogeographic criteria of sympatric speciation. In this study, we examined genome-wide variation in 46 individuals from these two groups. The divergence time between the GS and GE lineages was estimated to be 20-60 Kya. We identified 54 large genomic islands (≥100 kb) of speciation, which accounted for 89.4% of the total length of all genomic islands. These islands harboured divergent genes related to olfactory receptors and olfaction signals that may play important roles in food selection and assortative mating in fishes. Although the genomic islands clearly indicated speciation with gene flow and rejected micro-allopatric speciation, they were too large to support the hypothesis of sympatric speciation. Theoretical and recent empirical studies suggested that continual gene flow in sympatry should give rise to many small genomic islands (as small as a few kilobases in size). Thus, the observed pattern is consistent with the extensive evidence on parapatric speciation, in which adjacent habitats facilitate divergent selection but also permit gene flow during speciation. We suggest that many, if not most, of the reported cases of sympatric speciation are likely to be micro-parapatric speciation.

A phylogeographic, demographic and historical analysis of the short-tailed pit viper (Gloydius brevicaudus): evidence for early divergence and late expansion during the Pleistocene.

The impact of quaternary glaciation in eastern China on local fanua and flora has been a topic of considerable interest. We used mitochondrial DNA (mtDNA) sequence data and coalescent simulations to test two general biogeographic hypothesis related to the effects of Pleistocene climatic fluctuations for a widespread ophidian species (Gloydius brevicaudus) in eastern China and Korean Peninsula. The phylogenetic analysis revealed three major lineages, the southeast Coastal, Yangtze and North Lineages. The latter two are closely related and jointly form a continental lineage. Divergence dating and coalescent simulations indicate a Late Pliocene to Early Pleistocene divergence between lineages from the southeast coast and continental interior, followed by a mid-to-late Pleistocene divergence between lineages from the north and the middle-lower Yangtze Valley across East China, suggesting that all these lineages predated the last glacial maximum. An overlapping range between the two lineages within the continental lineage and a secondary contact associated with ecological transition zones on the margins of the North China Plain were also observed. These results show that vicariance patterns dominated the history of G. brevicaudus. Though the climatic events of the Pleistocene have had a marked effect on the historical distribution and intra-specific divergence of reptiles in China, coalescent and non-coalescent demographic analyses indicate that all lineages of G. brevicaudus seem not to have been adversely affected by glacial cycles during the Late Pleistocene, presumably because of an increase in the amount of climatically mild habitat in East Asia due to a decline in elevation and the development of monsoons since the Mid-End Pleistocene.

Phylogeography of the endemic Gymnocypris chilianensis (Cyprinidae): sequential westward colonization followed by allopatric evolution in response to cyclical Pleistocene glaciations on the Tibetan Plateau.

The schizothoracine Gymnocypris chilianensis is restricted to the Shiyang, Ruoshui and Shule Rivers, listed from east to west, along the northeast edge of the Tibetan Plateau. This distribution provides a valuable system to test hypotheses about postglacial colonization. We used mitochondrial DNA sequence data (a control region and the cytochrome b gene; 1894 bp) to assess the phylogeographic structure of this species based on 278 specimens sampled from throughout the species' entire geographical range. We found three lineages corresponding geographically to the three rivers, suggesting three independent glacial differentiation centers within the northeast edge of the Tibetan Plateau. The phylogenetic analysis suggested that the Shiyang River population forms a lineage that separated from the other populations of G. chilianensis at the basal phylogenetic split within this species. The molecular data further demonstrated a clear pattern of decreasing genetic diversity from the eastern Shiyang River towards the central Ruoshui River and western Shule River lineages, a pattern consistent with sequential western colonization. We therefore propose a phylogeographic scenario for G. chilianensis of a gradual westerly expansion from the Shiyang River population along the northeast edge of the Tibetan Plateau, with subsequent allopatric evolution at approximately 0.37 and 0.05 million years ago (Ma), through at least two glacial maxima. Together with the genetic evidence reported in other species, our findings suggest that this common biogeographic pattern emphasizes the importance of the northeastern edge region of the Tibetan Plateau as a hotspot of genetic diversity for some taxa.

The new chimeric chiron genes evolved essential roles in zebrafish embryonic development by regulating NAD+ levels.

The origination of new genes is important for generating genetic novelties for adaptive evolution and biological diversity. However, their potential roles in embryonic development, evolutionary processes into ancient networks, and contributions to adaptive evolution remain poorly investigated. Here, we identified a novel chimeric gene family, the chiron family, and explored its genetic basis and functional evolution underlying the adaptive evolution of Danioninae fishes. The ancestral chiron gene originated through retroposition of nampt in Danioninae 48-54 million years ago (Mya) and expanded into five duplicates (chiron1-5) in zebrafish 1-4 Mya. The chiron genes (chirons) likely originated in embryonic development and gradually extended their expression in the testis. Functional experiments showed that chirons were essential for zebrafish embryo development. By integrating into the NAD+ synthesis pathway, chirons could directly catalyze the NAD+ rate-limiting reaction and probably impact two energy metabolism genes (nmnat1 and naprt) to be under positive selection in Danioninae fishes. Together, these results mainly demonstrated that the origin of new chimeric chiron genes may be involved in adaptive evolution by integrating and impacting the NAD+ biosynthetic pathway. This coevolution may contribute to the physiological adaptation of Danioninae fishes to widespread and varied biomes in Southeast Asian.

Whole genome sequencing of silver carp (Hypophthalmichthys molitrix) and bighead carp (Hypophthalmichthys nobilis) provide novel insights into their evolution and speciation.

The edible silver carp (Hypophthalmichthys molitrix) and bighead carp (H. nobilis), which are two of the "Four Domesticated Fish" of China, are cultivated intensively worldwide. Here, we constructed 837- and 845-Mb draft genome assemblies for the silver carp and the bighead carp, respectively, including 24,571 and 24,229 annotated protein-coding genes. Genetic maps, anchoring 71.7% and 83.8% of all scaffolds, were obtained for the silver and bighead carp, respectively. Phylogenetic analysis showed that the bighead carp formed a clade with the silver carp, with an estimated divergence time of 3.6 million years ago; the time of divergence between the silver carp and zebrafish was 50.7 million years ago. An East Asian cyprinid genome-specific chromosome fusion took place ~9.2 million years after this clade diverged from the clade containing the common carp and Sinocyclocheilus. KEGG and GO analyses indicated that the expanded gene families in the silver and bighead carp were associated with diseases, the immune system and environmental adaptations. Genomic regions differentiating the silver and bighead carp populations were detected based on the whole-genome sequences of 42 individuals. Genes associated with the divergent regions were associated with reproductive system development and the development of primary female sexual characteristics. Thus, our results provided a novel systematic genomic analysis of the East Asian cyprinids, as well as the evolution and speciation of the silver carp and bighead carp.

Multiple source genes of HAmo SINE actively expanded and ongoing retroposition in cyprinid genomes relying on its partner LINE.

BACKGROUND: We recently characterized HAmo SINE and its partner LINE in silver carp and bighead carp based on hybridization capture of repetitive elements from digested genomic DNA in solution using a bead-probe 1. To reveal the distribution and evolutionary history of SINEs and LINEs in cyprinid genomes, we performed a multi-species search for HAmo SINE and its partner LINE using the bead-probe capture and internal-primer-SINE polymerase chain reaction (PCR) techniques. RESULTS: Sixty-seven full-size and 125 internal-SINE sequences (as well as 34 full-size and 9 internal sequences previously reported in bighead carp and silver carp) from 17 species of the family Cyprinidae were aligned as well as 14 new isolated HAmoL2 sequences. Four subfamilies (type I, II, III and IV), which were divided based on diagnostic nucleotides in the tRNA-unrelated region, expanded preferentially within a certain lineage or within the whole family of Cyprinidae as multiple active source genes. The copy numbers of HAmo SINEs were estimated to vary from 104 to 106 in cyprinid genomes by quantitative RT-PCR. Over one hundred type IV members were identified and characterized in the primitive cyprinid Danio rerio genome but only tens of sequences were found to be similar with type I, II and III since the type IV was the oldest subfamily and its members dispersed in almost all investigated cyprinid fishes. For determining the taxonomic distribution of HAmo SINE, inter-primer SINE PCR was conducted in other non-cyprinid fishes, the results shows that HAmo SINE- related sequences may disperse in other families of order Cypriniforms but absent in other orders of bony fishes: Siluriformes, Polypteriformes, Lepidosteiformes, Acipenseriformes and Osteoglossiforms. CONCLUSIONS: Depending on HAmo LINE2, multiple source genes (subfamilies) of HAmo SINE actively expanded and underwent retroposition in a certain lineage or within the whole family of Cyprinidae. From this perspective, HAmo SINE should provide useful phylogenetic makers for future analyses of the evolutionary relationships among species in the family Cyprinidae.

Genome size evolution in pufferfish: an insight from BAC clone-based Diodon holocanthus genome sequencing.

BACKGROUND: Variations in genome size within and between species have been observed since the 1950 s in diverse taxonomic groups. Serving as model organisms, smooth pufferfish possess the smallest vertebrate genomes. Interestingly, spiny pufferfish from its sister family have genome twice as large as smooth pufferfish. Therefore, comparative genomic analysis between smooth pufferfish and spiny pufferfish is useful for our understanding of genome size evolution in pufferfish. RESULTS: Ten BAC clones of a spiny pufferfish Diodon holocanthus were randomly selected and shotgun sequenced. In total, 776 kb of non-redundant sequences without gap representing 0.1% of the D. holocanthus genome were identified, and 77 distinct genes were predicted. In the sequenced D. holocanthus genome, 364 kb is homologous with 265 kb of the Takifugu rubripes genome, and 223 kb is homologous with 148 kb of the Tetraodon nigroviridis genome. The repetitive DNA accounts for 8% of the sequenced D. holocanthus genome, which is higher than that in the T. rubripes genome (6.89%) and that in the Te. nigroviridis genome (4.66%). In the repetitive DNA, 76% is retroelements which account for 6% of the sequenced D. holocanthus genome and belong to known families of transposable elements. More than half of retroelements were distributed within genes. In the non-homologous regions, repeat element proportion in D. holocanthus genome increased to 10.6% compared with T. rubripes and increased to 9.19% compared with Te. nigroviridis. A comparison of 10 well-defined orthologous genes showed that the average intron size (566 bp) in D. holocanthus genome is significantly longer than that in the smooth pufferfish genome (435 bp). CONCLUSION: Compared with the smooth pufferfish, D. holocanthus has a low gene density and repeat elements rich genome. Genome size variation between D. holocanthus and the smooth pufferfish exhibits as length variation between homologous region and different accumulation of non-homologous sequences. The length difference of intron is consistent with the genome size variation between D. holocanthus and the smooth pufferfish. Different transposable element accumulation is responsible for genome size variation between D. holocanthus and the smooth pufferfish.

Hox genes of the Japanese eel Anguilla japonica and Hox cluster evolution in teleosts.

Compared with other diploid teleosts (2n=48), anguilloid fish have a specialized karyotype (2n=38) and remarkable morphological variation, and represent one basal group species of teleosts. To investigate the Hox gene/cluster inventory in basal teleosts, a PCR-based survey of Hox genes in the Japanese eel (Anguilla japonica) was conducted with both gene-specific and homeobox-targeted degenerate primers. Our data provide evidence that at least 34 distinct Hox genes exist in the Japanese eel genome and that they represent eight Hox clusters. Duplication of Hox genes in the Japanese eel appears to be the result of the fish-specific genome duplication (FSGD) event. The Japanese eel shared the FSGD event with other teleosts such as zebrafish and pufferfish. A member of Hox paralog group one (HoxA1b) was preserved in the Japanese eel but was lost in other teleosts. Available Hox data revealed that the Hox cluster evolved distinctly in different teleost lineages. All duplicated Hox clusters were retained after the FSGD event in basal teleosts like in the Japanese eel, whereas crown teleosts lost one cluster (HoxCb or HoxDb). Based on current teleostean phylogeny, the HoxDb cluster was lost independently in the teleost lineages Otocephala and Euteleostei.

Distinct evolutionary patterns between two duplicated color vision genes within cyprinid fishes.

We investigated the molecular evolution of duplicated color vision genes (LWS-1 and SWS2) within cyprinid fish, focusing on the most cavefish-rich genus--Sinocyclocheilus. Maximum likelihood-based codon substitution approaches were used to analyze the evolution of vision genes. We found that the duplicated color vision genes had unequal evolutionary rates, which may lead to a related function divergence. Divergence of LWS-1 was strongly influenced by positive selection causing an accelerated rate of substitution in the proportion of pocket-forming residues. The SWS2 pigment experienced divergent selection between lineages, and no positively selected site was found. A duplicate copy of LWS-1 of some cyprinine species had become a pseudogene, but all SWS2 sequences remained intact in the regions examined in the cyprinid fishes examined in this study. The pseudogenization events did not occur randomly in the two copies of LWS-1 within Sinocyclocheilus species. Some cave species of Sinocyclocheilus with numerous morphological specializations that seem to be highly adapted for caves, retain both intact copies of color vision genes in their genome. We found some novel amino acid substitutions at key sites, which might represent interesting target sites for future mutagenesis experiments. Our data add to the increasing evidence that duplicate genes experience lower selective constraints and in some cases positive selection following gene duplication. Some of these observations are unexpected and may provide insights into the effect of caves on the evolution of color vision genes in fishes.

Massive Loss of Olfactory Receptors But Not Trace Amine-Associated Receptors in the World's Deepest-Living Fish (Pseudoliparis swirei).

Olfactory receptor repertoires show highly dynamic evolution associated with ecological adaptations in different species. The Mariana snailfish (Pseudoliparis swirei) living below a depth of 6000 m in the Mariana Trench evolved degraded vision and occupies a specific feeding habitat in a dark, low-food environment. However, whether such adaptations involve adaptive changes in the chemosensory receptor repertoire is not known. Here, we conducted a comparative analysis of the olfactory receptor (OR) and trace amine-associated receptor (TAAR) gene repertoires in nine teleosts with a focus on the evolutionary divergence between the Mariana snailfish and its shallow-sea relative, Tanaka's snailfish (Liparis tanakae). We found many fewer functional OR genes and a significantly higher fraction of pseudogenes in the Mariana snailfish, but the numbers of functional TAAR genes in the two species were comparable. Phylogenetic analysis showed that the expansion patterns of the gene families were shared by the two species, but that Mariana snailfish underwent massive gene losses in its OR repertoire. Despite an overall decreased size in OR subfamilies and a reduced number of TAAR subfamilies in the Mariana snailfish, expansion of certain subfamilies was observed. Selective pressure analysis indicated greatly relaxed selective strength in ORs but a slightly enhanced selective strength in TAARs of Mariana snailfish. Overall, our study reveals simplified but specific OR and TAAR repertoires in the Mariana snailfish shaped by natural selection with respect to ecological adaptations in the hadal environment. This is the first study on the chemosensation evolution in vertebrates living in the hadal zone, which could provide new insights into evolutionary adaptation to the hadal environment.

DNA barcoding the ichthyofauna of the Yangtze River: Insights from the molecular inventory of a mega-diverse temperate fauna.

Intensification of inland fisheries and aquatic landscape conversion led to a drastic decline of fish populations in the Yangtze River (YR) during the last decades. This situation urges for the development of a large-scale molecular assessment of YR ichthyofauna to further develop standardized methods of molecular identification for conservation and fisheries management purposes. We present here the results of a large-scale campaign to DNA barcode YR freshwater fishes that succeeded in producing 1,424 new DNA barcodes for 123 species. Together with 1,406 sequences mined from BOLD and GenBank, a reference library including 2,830 DNA barcodes for 238 species was compiled. By using four DNA-based species delimitation methods, RESL, ABGD, mPTP and mGMYC, 230 operational taxonomic units (OTUs) were identified and 195 species displayed OTUs that tightly match species boundaries. No barcoding gap was observed; however, and conflicting cases of species and OTU delimitation were identified. A total of 23 species with maximum intraspecific distances above 2% were detected and null genetic distances to the nearest phylogenetic relatives were detected in 11 species. Among those 23 species, 16 were represented by multiple OTUs amounting to 40 OTUs delineated. Several cases of multiple OTUs confined to species boundaries were detected suggesting the presence of overlooked species. A total of 18 OTUs, however, were shared by several species and particularly so for the Qinghai-Tibet plateau endemic species. These results are discussed with reference to previous large-scale DNA barcoding campaign and compared to previous phylogeographic studies in the YR.