The sequence and de novo assembly of Oxygymnocypris stewartii genome

Chromosome-scale assembly of Artemia tibetiana genome, first aquatic invertebrate genome from Tibet Plateau

Genomic-level studies on the adaptive evolution of animals in the Qinghai-Tibet Plateau have been rapidly increasing. However, most studies are concentrated on vertebrates, and there are few reports on invertebrates. Here, we report the chromosome-level genome assembly for the brine shrimp Artemia tibetiana from Kyêbxang Co, a high-altitude (4620 m above sea level) salt lake on the plateau, based on the combination of Illumina, Nanopore long-reads and Hi-C sequencing data. The assembled genome is 1.69 Gb, and 94.83% of the assembled sequences are anchored to 21 pseudo-chromosomes. Approximately 75% of the genome was identified as repetitive sequences, which is higher than most crustaceans documented so far. A total of 17,988 protein-coding genes were identified, among them 14,388 were functionally annotated. This genomic resource provides the foundation for whole-genome level investigation on the genetic adaptation of Artemia to the harsh conditions in the Qinghai-Tibet Plateau.

The First Genome-Wide Survey Analysis of the Tibetan Plateau Tetraploid Schizothorax curvilabiatus Reveals Its Microsatellite Characteristics and Phylogenetic Relationships

Background/Objectives: Schizothorax curvilabiatus, a typical highland polyploid species within the subfamily Schizothoracinae, holds economic value and ecological research significance. Currently, there are no related genomic studies. To obtain its genetic information and lay the foundation for subsequent whole-genome map construction, this study conducted a genome survey analysis, preliminary genome assembly, microsatellite identification, repeat sequence annotation, mitochondrial genome characterization, and phylogenetic relationship research. Methods: DNA was sequenced on a DNBSEQ-T7 platform to obtain paired-end genomic data. The genome was analyzed using GCE, and the draft genome was assembled with SOAPdenovo. Microsatellites were identified using MISA, and the mitochondrial genome was assembled with NOVOPlasty. Genome features were analyzed, and phylogenetic trees were constructed using PhyloSuite and MEGA. Results: The genome size was estimated at 2.53 Gb, with a heterozygosity of 6.55% and 47.66% repeat sequences. A 1.324 Gb preliminary genome draft was obtained, with repeat sequences comprising 47.17%, the majority being DNA transposons (24.64%). Dinucleotide repeats were most abundant (46.91%), followed by mononucleotide repeats (38.31%), with A/T and AC/GT being the most frequent. A complete mitochondrial genome of 16,589 bp was assembled, and a 939 bp D-loop was annotated. Phylogenetic relationships among genera in the Schizothoracinae subfamily were also clarified. Conclusions: This study provides the latest molecular data for analysis of the S. curvilabiatus genome and its related populations, and for the first time offers genomic resources for research on genomic adaptive evolution and polyploidization in high-altitude environments.

A chromosome-level genome assembly of big-barbel schizothorcin, Schizothorax macropogon

Big-barbel schizothorcin (Schizothorax macropogon), an endemic and vulnerable species to the mid-reaches of the Yarlung Zangbo River, epitomizes survival in harsh conditions yet suffers significant population contractions due to human activities. This species was the subject of our study in which we leveraged PacBio, MGI-Seq, and Hi-C data to assemble a chromosome-scale genome. This assembly comprises 25 pseudo-chromosomes, yielding a genome size of 1.42 Gb with a scaffold N50 length of 59.4 Mb, indicative of a highly contiguous assembly. A BUSCO assessment ascertained the comprehensiveness of the genome at 97.9%. Annotation efforts identified 46,246 putative protein-coding genes, with 49.61% of the assembled genome annotated as repetitive sequences. This genome assembly is instrumental for advancing conservation of the giant whiskered schizothoracines and related species, and for illuminating the evolution and ecology of schizothoracine fishes in the Qinghai-Tibet Plateau.

Genetic Diversity and Population Structure of Rhodeus uyekii in the Republic of Korea Revealed by Microsatellite Markers from Whole Genome Assembly

This study is the first report to characterize the Rhodus uyekii genome and study the development of microsatellite markers and their markers applied to the genetic structure of the wild population. Genome assembly was based on PacBio HiFi and Illumina HiSeq paired-end sequencing, resulting in a draft genome assembly of R. uyekii. The draft genome was assembled into 2652 contigs. The integrity assessment of the assemblies indicates that the quality of the draft assemblies is high, with 3259 complete BUSCOs (97.2%) in the database of Verbrata. A total of 31,166 predicted protein-coding genes were annotated in the protein database. The phylogenetic tree showed that R. uyekii is a close but distinct relative of Onychostoma macrolepis. Among the 10 fish genomes, there were significant gene family expansions (8-2387) and contractions (16-2886). The average number of alleles amplified by the 21 polymorphic markers ranged from 6 to 23, and the average PIC value was 0.753, which will be useful for evolutionary and genetic analysis. Using population genetic analysis, we analyzed genetic diversity and the genetic structures of 120 individuals from 6 populations. The average number of alleles per population ranged from 7.6 to 9.9, observed heterozygosity ranged from 0.496 to 0.642, and expected heterozygosity ranged from 0.587 to 0.783. Discriminant analysis of principal components According to the analysis method, the population was divided into three populations (BS vs. DC vs. GG, GC, MS, DC). In conclusion, our study provides a useful resource for comparative genomics, phylogeny, and future population studies of R. uyekii.

Evolutionary divergence of subgenomes in common carp provides insights into speciation and allopolyploid success

Hybridization and polyploidization have made great contributions to speciation, heterosis, and agricultural production within plants, but there is still limited understanding and utilization in animals. Subgenome structure and expression reorganization and cooperation post hybridization and polyploidization are essential for speciation and allopolyploid success. However, the mechanisms have not yet been comprehensively assessed in animals. Here, we produced a high-fidelity reference genome sequence for common carp, a typical allotetraploid fish species cultured worldwide. This genome enabled in-depth analysis of the evolution of subgenome architecture and expression responses. Most genes were expressed with subgenome biases, with a trend of transition from the expression of subgenome A during the early stages to that of subgenome B during the late stages of embryonic development. While subgenome A evolved more rapidly, subgenome B contributed to a greater level of expression during development and under stressful conditions. Stable dominant patterns for homoeologous gene pairs both during development and under thermal stress suggest a potential fixed heterosis in the allotetraploid genome. Preferentially expressing either copy of a homoeologous gene at higher levels to confer development and response to stress indicates the dominant effect of heterosis. The plasticity of subgenomes and their shifting of dominant expression during early development, and in response to stressful conditions, provide novel insights into the molecular basis of the successful speciation, evolution, and heterosis of the allotetraploid common carp.

Chromosome-level genome of butterflyfish unveils genomic features of unique colour patterns and morphological traits

Chaetodontidae, known as butterflyfishes, are typical fish in coral ecosystems, exhibiting remarkable interspecific differences including body colour patterns and feeding ecology. In this study, we report genomes of three butterflyfish species (Chelmon rostratus, Chaetodon trifasciatus and Chaetodon auriga) and a closely related species from the Pomacanthidae family, Centropyge bicolour, with an average genome size of 65,611 Mb. Chelmon rostratus, comprising 24 chromosomes assembled to the chromosome level, could be served as a reference genome for butterflyfish. By conducting a collinearity analysis between butterflyfishes and several fishes, we elucidated the specific and conserved genomic features of butterflyfish, with particular emphasis on novel genes arising from tandem duplications and their potential functions. In addition to the two melanocyte-specific tyr genes commonly found in fish, we found the gene tyrp3, a new tyrosinase-related proteins gene in the reef fish, including butterflyfish and clownfish, implicating their involvement in the pigmentation diversity of fish. Additionally, we observed a tandem duplication expansion of three copies of nell1 gene in C. rostratus genome, which likely contribute to its unique jaw development and distinctive morphology of its sharp mouth. These results provided valuable genomic resources for further investigations into the genetic diversity and evolutionary adaptations of reef fish.

Genetic Variation in Schizothorax kozlovi Nikolsky in the Upper Reaches of the Chinese Yangtze River Based on Genotyping for Simplified Genome Sequencing

Simple Summary

Schizothorax kozlovi Nikolsky is a unique cold−water fish in the upper reaches of the Yangtze River in China and has high economic value. In our study, genetic diversity and population structure analyses were performed on seven wild populations in the upper reaches of the Yangtze River by GBS. The above results indicate that the populations of S. kozlovi have different degrees of tolerance and selection pressure in response to temperature and altitude. The Wujiang population was genetically differentiated from the Jinsha River and Yalong River populations. The Wujiang intrapopulation has greater genetic diversity and differentiation than the Jinsha River and Yalong River populations, which demonstrates that the Jinsha and Yalong populations require more attention and resources for their protection. The results of this study will increase our understanding of the diversity of S. kozlovi in the upper reaches of the Yangtze River and provide a basis for the conservation and utilization of wild resources.

Abstract

Schizothorax kozlovi Nikolsky is a unique cold−water fish in the upper reaches of the Yangtze River in China and has high economic value. In our study, genetic diversity and population structure analyses were performed on seven wild populations (originating from the Jinsha River, Yalong River, and Wujiang River) in the upper reaches of the Yangtze River by genotyping by sequencing (GBS). The results indicated that a total of 303,970 single−nucleotide polymorphisms (SNPs) were identified from the seven wild populations. Lower genetic diversity was exhibited among the intrapopulations of the three tributaries, and the Wujiang River population had significant genetic differentiation when compared to the Jinsha River and Yalong River populations. Furthermore, the selected SNPs were enriched in cellular processes, environmental adaptation, signal transduction, and related metabolic processes between the Wujiang population and the other two populations. The above results indicate that the populations of S. kozlovi have different degrees of tolerance and selection pressure in response to temperature and altitude. The Wujiang intrapopulation has greater genetic diversity and differentiation than the Jinsha River and Yalong River populations, which demonstrates that the Jinsha and Yalong populations require more attention and resources for their protection. The results of this study will increase our understanding of the diversity of S. kozlovi in the upper reaches of the Yangtze River and provide a basis for the conservation and utilization of wild resources.

Chromosome-level genome of Tibetan naked carp (Gymnocypris przewalskii) provides insights into Tibetan highland adaptation

Gymnocypris przewalskii, a cyprinid fish endemic to the Qinghai-Tibetan Plateau, has evolved unique morphological, physiological and genetic characteristics to adapt to the highland environment. Herein, we assembled a high-quality G. przewalskii tetraploid genome with a size of 2.03 Gb and scaffold N50 of 44.93 Mb, which was anchored onto 46 chromosomes. The comparative analysis suggested that gene families related to highland adaptation were significantly expanded in G. przewalskii. According to the G. przewalskii genome, we evaluated the phylogenetic relationship of 13 schizothoracine fishes, and inferred that the demographic history of G. przewalskii was strongly associated with geographic and eco-environmental alterations. We noticed that G. przewalskii experienced whole-genome duplication, and genes preserved post duplication were functionally associated with adaptation to high salinity and alkalinity. In conclusion, a chromosome-scale G. przewalskii genome provides an important genomic resource for teleost fish, and will particularly promote our understanding of the molecular evolution and speciation of fish in the highland environment.

Temperature Acclimation Alters the Thermal Tolerance and Intestinal Heat Stress Response in a Tibetan Fish Oxygymnocypris stewarti

Numerous studies have shown that thermal tolerance and intestinal heat resistance are strongly associated with temperature acclimation. However, few reports have successfully conducted similar research on fishes from the Qinghai–Tibetan Plateau, an area that is facing the threat of climate warming. Therefore, the present study determined the growth, thermal tolerance, and intestinal heat stress (exposure to 30°C) responses in juveniles of a Tibetan fish, Oxygymnocypris stewarti, acclimated to three temperature levels (10°C, 15°C, and 20°C, named as T10, T15, and T20, respectively) for 30 days. The fastest growth was recorded in the T15 group. At 1°C/30 min heating rate, the critical thermal maximum (CT_Max) ranged from 31.3°C to 32.3°C, and the lethal thermal maximum (LT_Max) ranged from 31.8°C to 32.6°C among the three acclimation temperatures. According to the results of thermal tolerance tests, the heat stress temperature was set to 30°C. When the water temperature reached 30°C, the expression of the intestinal heat shock protein 70 (HSP70) gene as well as the intestinal microbiome and histology of experimental fish were monitored at 0, 2, 6, and 12 h. The expression of HSP70 reached the highest level at 2 h in all three temperature treatments. The histological analysis showed damage to intestinal cells, including diffuse infiltration of lymphocytes, villi epithelial cell swelling, decrease of intestinal villi length, and cytoplasmic light staining at 2 h in all three temperature treatments. In terms of the intestinal microbiome, phyla Proteobacteria and Firmicutes dominated the treatments at each monitored time in the T10 and T15 groups and at 0 h in T20 group, while phyla Fusobacteria, Proteobacteria, and Cyanobacteria were dominant in treatments at 2, 6, and 12 h in the T20 group. The overall results indicated that acclimation temperature could affect the growth, thermal tolerance, and intestinal heat stress response of O. stewarti juveniles. As the first report on intestinal heat stress response associated with temperature acclimation in a Tibetan fish, this study will help to understand the potential effects of climate change on highland fishes.

Whole-Genome Sequencing of Sinocyclocheilus maitianheensis Reveals Phylogenetic Evolution and Immunological Variances in Various Sinocyclocheilus Fishes

An adult Sinocyclocheilus maitianheensis, a surface-dwelling golden-line barbel fish, was collected from Maitian river (Kunming City, Yunnan Province, China) for whole-genome sequencing, assembly, and annotation. We obtained a genome assembly of 1.7 Gb with a scaffold N50 of 1.4 Mb and a contig N50 of 24.7 kb. A total of 39,977 protein-coding genes were annotated. Based on a comparative phylogenetic analysis of five Sinocyclocheilus species and other five representative vertebrates with published genome sequences, we found that S. maitianheensis is close to Sinocyclocheilus anophthalmus (a cave-restricted species with similar locality). Moreover, the assembled genomes of S. maitianheensis and other four Sinocyclocheilus counterparts were used for a fourfold degenerative third-codon transversion (4dTv) analysis. The recent whole-genome duplication (WGD) event was therefore estimated to occur about 18.1 million years ago. Our results also revealed a decreased tendency of copy number in many important genes related to immunity and apoptosis in cave-restricted Sinocyclocheilus species. In summary, we report the first genome assembly of S. maitianheensis, which provides a valuable genetic resource for comparative studies on cavefish biology, species protection, and practical aquaculture of this potentially economical fish.

Genomic and functional evidence reveals convergent evolution in fishes on the Tibetan Plateau

High-altitude environments are strong drivers of adaptive evolution in endemic organisms. However, little is known about the genetic mechanisms of convergent adaptation among different lineages, especially in fishes. There are three independent fish groups on the Tibetan Plateau: Tibetan Loaches, Schizothoracine fishes and Glyptosternoid fishes; all are well adapted to the harsh environmental conditions. They represent an excellent example of convergent evolution but with an unclear genetic basis. We used comparative genomic analyses between Tibetan fishes and fishes from low altitudes and detected genomic signatures of convergent evolution in fishes on the Tibetan Plateau. The Tibetan fishes exhibited genome-wide accelerated evolution in comparison with a control set of fishes from low altitudes. A total of 368 positively selected genes were identified in Tibetan fishes, which were enriched in functional categories related to energy metabolism and hypoxia response. Widespread parallel amino acid substitutions were detected among the Tibetan fishes and a subset of these substitutions occurred in positively selected genes associated with high-altitude adaptation. Functional assays suggested that von Hippel-Lindau (VHL) tumour suppressor genes from Tibetan fishes enhance hypoxia-inducible factor (HIF) activity convergently under hypoxia compared to low-altitude fishes. The results provide genomic and functional evidence supporting convergent genetic mechanisms for high-altitude adaptation in fishes on the Tibetan Plateau.

Sequencing an F1 hybrid of Silurus asotus and S. meridionalis enabled the assembly of high-quality parental genomes

Genome complexity such as heterozygosity may heavily influence its de novo assembly. Sequencing somatic cells of the F1 hybrids harboring two sets of genetic materials from both of the paternal and maternal species may avoid alleles discrimination during assembly. However, the feasibility of this strategy needs further assessments. We sequenced and assembled the genome of an F1 hybrid between Silurus asotus and S. meridionalis using the SequelII platform and Hi-C scaffolding technologies. More than 300 Gb raw data were generated, and the final assembly obtained 2344 scaffolds composed of 3017 contigs. The N50 length of scaffolds and contigs was 28.55 Mb and 7.49 Mb, respectively. Based on the mapping results of short reads generated for the paternal and maternal species, each of the 29 chromosomes originating from S. asotus and S. meridionalis was recognized. We recovered nearly 94% and 96% of the total length of S. asotus and S. meridionalis. BUSCO assessments and mapping analyses suggested that both genomes had high completeness and accuracy. Further analyses demonstrated the high collinearity between S. asotus, S. meridionalis, and the related Pelteobagrus fulvidraco. Comparison of the two genomes with that assembled only using the short reads from non-hybrid parental species detected a small portion of sequences that may be incorrectly assigned to the different species. We supposed that at least part of these situations may have resulted from mitotic recombination. The strategy of sequencing the F1 hybrid genome can recover the vast majority of the parental genomes and may improve the assembly of complex genomes.

Hb adaptation to hypoxia in high-altitude fishes: Fresh evidence from schizothoracinae fishes in the Qinghai-Tibetan Plateau

Uncovering the genetic basis of hypoxic adaptation is one of the most active research areas in evolutionary biology. Among air-breathing vertebrates, modifications of hemoglobin (Hb) play a pivotal role in mediating an adaptive response to high-altitude hypoxia. However, the relative contributions in water-breathing organisms are still unclear. Here, we tested the Hb concentration of fish at different altitudes. All species showed species-specific Hb concentration, which has a non-positive correlation with altitude. Moreover, we investigated the expression of Hb genes by the RNA-seq and quantitative real-time PCR (qRT-PCR), and Hb composition by two-dimensional electrophoresis (2-DE). The results showed that the multiple Hb genes and isoforms are co-expressed in schizothoracinae fishes endemic to the Qinghai-Tibetan Plateau (QTP). Phylogenetic analyses of Hb genes indicated that the evolutionary relationships are not easily reconciled with the organismal phylogeny. Furthermore, evidence of positive selection was found in the Hb genes of schizothoracinae fishes through the selection pressure analysis. We demonstrated that positively selected sites likely facilitated the functional divergence of Hb isoforms. Taken together, this study indicated that the long-term maintenance of high Hb concentration may be a disadvantage for physiologically acclimating to high altitude hypoxia. Meanwhile, the genetically based modification of Hb-O₂ affinity in schizothoracinae fishes might facilitate the evolutionary adaptation to Tibetan aqueous environments.

Genomic Signature of Shifts in Selection and Alkaline Adaptation in Highland Fish

Understanding how organisms adapt to aquatic life at high altitude is fundamental in evolutionary biology. This objective has been addressed primarily related to hypoxia adaptation by recent comparative studies, whereas highland fish has also long suffered extreme alkaline environment, insight into the genomic basis of alkaline adaptation has rarely been provided. Here, we compared the genomes or transcriptomes of 15 fish species, including two alkaline tolerant highland fish species and their six alkaline intolerant relatives, three alkaline tolerant lowland fish species, and four alkaline intolerant species. We found putatively consistent patterns of molecular evolution in alkaline tolerant species in a large number of shared orthologs within highland and lowland fish taxa. Remarkably, we identified consistent signatures of accelerated evolution and positive selection in a set of shared genes associated with ion transport, apoptosis, immune response, and energy metabolisms in alkaline tolerant species within both highland and lowland fish taxa. This is one of the first comparative studies that began to elucidate the consistent genomic signature of alkaline adaptation shared by highland and lowland fish. This finding also highlights the adaptive molecular evolution changes that support fish adapting to extreme environments at high altitude.

Genome sequence and transcriptome profiles of pathogenic fungus Paecilomyces penicillatus reveal its interactions with edible fungus Morchella importuna

Paecilomyces penicillatus is one of the pathogens of morels, which greatly affects the yield and quality of Morchella spp.. In the present study, we de novo assembled the genome sequence of the fungus P. penicillatus SAAS_ppe1. We analyzed the transcriptional profile of P. penicillatus SAAS_ppe1 infection of Morchella importuna at different stages (3 days and 6 days after infection) and the response of M. importuna using the transcriptome. The assembled genome sequence of P. penicillatus SAAS_ppe1 was 39.78 Mb in length (11 scaffolds; scaffold N50, 6.50 Mb), in which 99.7% of the expected genes were detected. A total of 7.48% and 19.83% clean transcriptional reads from the infected sites were mapped to the P. penicillatus genome at the early and late stages of infection, respectively. There were 3,943 genes differently expressed in P. penicillatus at different stages of infection, of which 24 genes had increased expression with the infection and infection stage, including diphthamide biosynthesis, aldehyde reductase, and NAD (P)H-hydrate epimerase (P < 0.05). Several genes had variable expression trends at different stages of infection, indicating P. penicillatus had diverse regulation patterns to infect M. importuna. GO function, involving cellular components, and KEGG pathways, involving glycerolipid metabolism, and plant-pathogen interaction were significantly enriched during infection by P. penicillatus. The expression of ten genes in M. importuna increased during the infection and infection stage, and these may regulate the response of M. importuna to P. penicillatus infection. This is the first comprehensive study on P. penicillatus infection mechanism and M. importuna response mechanism, which will lay a foundation for understanding the fungus-fungus interactions, gene functions, and variety breeding of pathogenic and edible fungi.

Chromosome-level genome assembly of the Arctic fox (Vulpes lagopus) using PacBio sequencing and Hi-C technology

The Arctic fox (Vulpes lagopus) is the only fox species occurring in the Arctic and has adapted to its extreme climatic conditions. Currently, the molecular basis of its adaptation to the extreme climate has not been characterized. Here, we applied PacBio sequencing and chromosome structure capture technique to assemble the first V. lagopus genome assembly, which is assembled into chromosome fragments. The genome assembly has a total length of 2.345 Gb with a contig N50 of 31.848 Mb and a scaffold N50 of 131.537 Mb, consisting of 25 pseudochromosomal scaffolds. The V. lagopus genome had approximately 32.33% repeat sequences. In total, 21,278 protein-coding genes were predicted, of which 99.14% were functionally annotated. Compared with 12 other mammals, V. lagopus was most closely related to V. Vulpes with an estimated divergence time of ~7.1 Ma. The expanded gene families and positively selected genes potentially play roles in the adaptation of V. lagopus to Arctic extreme environment. This high-quality assembled genome will not only promote future studies of genetic diversity and evolution in foxes and other canids but also provide important resources for conservation of Arctic species.

Chromosome-level genome of Poropuntius huangchuchieni provides a diploid progenitor-like reference genome for the allotetraploid Cyprinus carpio

The diploid Poropuntius huangchuchieni in the cyprinid family, which is widely distributed in the Mekong and Red River basins, is one of the most closely related diploid progenitor-like species of allotetraploid common carp, which was generated by merging of two diploid genomes during evolution. Therefore, the P. huangchuchieni genome is essential for polyploid evolution studies in Cyprinidae. Here, we report a high-quality chromosome-level genome assembly of P. huangchuchieni by integrating Oxford Nanopore and Hi-C technologies. The assembled genome size was 1,021.38 Mb, 895.66 Mb of which was anchored onto 25 chromosomes with a N50 of 32.93 Mb. The genome contained 486.28 Mb repetitive elements and 24,099 protein-coding genes. Approximately 95.9% of the complete BUSCOs were detected, suggesting a high completeness of the genome. Evolutionary analysis revealed that P. huangchuchieni diverged from Cyprinus carpio at approximately 12 Mya. Genome comparison between P. huangchuchieni and the B subgenome of C. carpio provided insights into chromosomal rearrangements during the allotetraploid speciation. With the complete gene set, 17,474 orthologous genes were identified between P. huangchuchieni and C. carpio, providing a broad view of the gene component in the allotetraploid genome, which is critical for future genetic analyses. The high-quality genomic data set created for P. huangchuchieni provides a diploid progenitor-like reference for the evolution and adaptation of allotetraploid carps.

Comparative transcriptome analysis of scaled and scaleless skins in Gymnocypris eckloni provides insights into the molecular mechanism of scale degeneration

Background

The scale degeneration is thought to be related to the adaptation to the extreme environment with cold climate and high-altitude in schizothoracine fishes. Gymnocypris eckloni, a schizothoracine fish living in plateau waters with the elevation above 2500 m, is nearly esquamate and only covered with shoulder scales and anal scales, making it a good model species to study the molecular mechanism of scale degeneration.

Results

The transcriptomes of shoulder scaled skins (SSS), anal scaled skins (ASS) and scaleless skins (NSS) were sequenced and analyzed in G. eckloni at the age of 1 year. Histological examination showed that shoulder scale had completed its differentiation and anal scale just initiated the differentiation. A total of 578,046 unigenes were obtained from the transcriptomes, with 407,799 unigenes annotated in public databases. A total of 428 and 142 differentially expressed unigenes (DEUs) were identified between SSS and NSS, and between ASS and NSS, respectively, with 45 DEUs that were overlapped. Annotation analysis indicated that these DEUs were mainly enriched in Gene Ontology (GO) terms and KEGG pathways associated with bone and muscle formation, such as myofibril, contractile fiber, cytoskeletal protein binding, muscle structure development, cardiac muscle contraction, hypertrophic cardiomyopathy (HCM) and calcium signaling pathway.

Conclusions

Our results would provide insights into the molecular mechanisms of scale degeneration in G. eckloni and other congeneric fishes. In addition, the transcriptome data provides candidate genes and markers for future studies.

Chromosomal genome of Triplophysa bleekeri provides insights into its evolution and environmental adaptation

BACKGROUND

Intense stresses caused by high-altitude environments may result in noticeable genetic adaptions in native species. Studies of genetic adaptations to high elevations have been largely limited to terrestrial animals. How fish adapt to high-elevation environments is largely unknown. Triplophysa bleekeri, an endemic fish inhabiting high-altitude regions, is an excellent model to investigate the genetic mechanisms of adaptation to the local environment. Here, we assembled a chromosomal genome sequence of T. bleekeri, with a size of ∼628 Mb (contig and scaffold N50 of 3.1 and 22.9 Mb, respectively). We investigated the origin and environmental adaptation of T. bleekeri based on 21,198 protein-coding genes in the genome.

RESULTS

Compared with fish species living at low altitudes, gene families associated with lipid metabolism and immune response were significantly expanded in the T. bleekeri genome. Genes involved in DNA repair exhibit positive selection for T. bleekeri, Triplophysa siluroides, and Triplophysa tibetana, indicating that adaptive convergence in Triplophysa species occurred at the positively selected genes. We also analyzed whole-genome variants among samples from 3 populations. The results showed that populations separated by geological and artificial barriers exhibited obvious differences in genetic structures, indicating that gene flow is restricted between populations.

CONCLUSIONS

These results will help us expand our understanding of environmental adaptation and genetic diversity of T. bleekeri and provide valuable genetic resources for future studies on the evolution and conservation of high-altitude fish species such as T. bleekeri.

Recent genome duplications facilitate the phenotypic diversity of Hb repertoire in the Cyprinidae

Whole-genome duplications (WGDs) are an important contributor to phenotypic innovations in evolutionary history. The diversity of blood oxygen transport traits is the perfect reflection of physiological versatility for evolutionary success among vertebrates. In this study, the evolutionary changes of hemoglobin (Hb) repertoire driven by the recent genome duplications were detected in representative Cyprinidae fish, including eight diploid and four tetraploid species. Comparative genomic analysis revealed a substantial variation in both membership composition and intragenomic organization of Hb genes in these species. Phylogenetic reconstruction analyses were conducted to characterize the evolutionary history of these genes. Data were integrated with the expression profiles of the genes during ontogeny. Our results indicated that genome duplications facilitated the phenotypic diversity of the Hb gene family; each was associated with species-specific changes in gene content via gene loss and fusion after genome duplications. This led to repeated evolutionary transitions in the ontogenic regulation of Hb gene expression. Our results revealed that genome duplications helped to generate phenotypic changes in Cyprinidae Hb systems.

Sequencing of the black rockfish chromosomal genome provides insight into sperm storage in the female ovary

Black rockfish (Sebastes schlegelii) is an economically important viviparous marine teleost in Japan, Korea, and China. It is characterized by internal fertilization, long-term sperm storage in the female ovary, and a high abortion rate. For better understanding the mechanism of fertilization and gestation, it is essential to establish a reference genome for viviparous teleosts. Herein, we used a combination of Pacific Biosciences sequel, Illumina sequencing platforms, 10× Genomics, and Hi-C technology to obtain a genome assembly size of 848.31 Mb comprising 24 chromosomes, and contig and scaffold N50 lengths of 2.96 and 35.63 Mb, respectively. We predicted 39.98% repetitive elements, and 26,979 protein-coding genes. S. schlegelii diverged from Gasterosteus aculeatus ∼32.1-56.8 million years ago. Furthermore, sperm remained viable within the ovary for up to 6 months. The glucose transporter SLC2 showed significantly positive genomic selection, and carbohydrate metabolism-related KEGG pathways were significantly up-regulated in ovaries after copulation. In vitro suppression of glycolysis with sodium iodoacetate reduced sperm longevity significantly. The results indicated the importance of carbohydrates in maintaining sperm survivability. Decoding the S. schlegelii genome not only provides new insights into sperm storage; additionally, it is highly valuable for marine researchers and reproduction biologists.

Comprehensive transcriptome data for endemic Schizothoracinae fish in the Tibetan Plateau

The Schizothoracinae fishes, endemic species in the Tibetan Plateau, are considered as ideal models for highland adaptation and speciation investigation. Despite several transcriptome studies for highland fishes have been reported before, the transcriptome information of Schizothoracinae is still lacking. To obtain comprehensive transcriptome data for Schizothoracinae, the transcriptome of a total of 183 samples from 14 representative Schizothoracinae species, were sequenced and de novo assembled. As a result, about 1,363 Gb transcriptome clean data was obtained. After the assembly, we obtain 76,602-154,860 unigenes for each species with sequence N50 length of 1,564-2,143 bp. More than half of the unigenes were functionally annotated by public databases. The Schizothoracinae fishes in this work exhibited diversified ecological distributions, phenotype characters and feeding habits; therefore, the comprehensive transcriptome data of those species provided valuable information for the environmental adaptation and speciation of Schizothoracinae in the Tibetan Plateau.