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.

Chromosome-level genome provides insight into the evolution and conservation of the threatened goral (Naemorhedus goral).

BACKGROUND: Gorals Naemorhedus resemble both goats and antelopes, which prompts much debate about the intragenus species delimitation and phylogenetic status of the genus Naemorhedus within the subfamily Caprinae. Their evolution is believed to be linked to the uplift of the Qinghai-Tibet Plateau (QTP). To better understand its phylogenetics, the genetic information is worth being resolved. RESULTS: Based on a sample from the eastern margin of QTP, we constructed the first reference genome for Himalayan goral Naemorhedus goral, using PacBio long-read sequencing and Hi-C technology. The 2.59 Gb assembled genome had a contig N50 of 3.70 Mb and scaffold N50 of 106.66 Mb, which anchored onto 28 pseudo chromosomes. A total of 20,145 protein-coding genes were predicted in the assembled genome, of which 99.93% were functionally annotated. Phylogenetically, the goral was closely related to muskox on the mitochondrial genome level and nested into the takin-muskox clade on the genome tree, rather than other so-called goat-antelopes. The cladogenetic event among muskox, takin and goral occurred sequentially during the late Miocene (~ 11 - 5 Mya), when the QTP experienced a third dramatic uplift with consequent profound changes in climate and environment. Several chromosome fusions and translocations were observed between goral and takin/muskox. The expanded gene families in the goral genome were mainly related to the metabolism of drugs and diseases, so as the positive selected genes. The Ne of goral continued to decrease since ~ 1 Mya during the Pleistocene with active glaciations. CONCLUSION: The high-quality goral genome provides insights into the evolution and valuable information for the conservation of this threatened group.

The Genomes of Two Billfishes Provide Insights into the Evolution of Endothermy in Teleosts.

Endothermy is a typical convergent phenomenon which has evolved independently at least eight times in vertebrates, and is of significant advantage to organisms in extending their niches. However, how vertebrates other than mammals or birds, especially teleosts, achieve endothermy has not previously been fully understood. In this study, we sequenced the genomes of two billfishes (swordfish and sailfish), members of a representative lineage of endothermic teleosts. Convergent amino acid replacements were observed in proteins related to heat production and the visual system in two endothermic teleost lineages, billfishes and tunas. The billfish-specific genetic innovations were found to be associated with heat exchange, thermoregulation, and the specialized morphology, including elongated bill, enlarged dorsal fin in sailfish and loss of the pelvic fin in swordfish.

Mitochondrial genome and phylogenetic relationship of Gymnocypris eckloni (Schizothoracinae) in Qaidam river basin.

Genomic data can improve our understanding on the phylogenetic relationship among Tibetan highland fishes. The whole mitochondrial genome of Gymnocypris eckloni generated in this study is 16,784 bp in length, containing 22 transfer RNA genes, 13 protein-coding genes, two ribosomal RNA genes, and one non-coding control region (D-Loop). Phylogenetic analysis recovers a non-monophyetic population of G. eckloni in the Qaidam basin, representing two distinct lineages designated 'Qiadam A' and 'Qaidam B', within Qaidam A clustering with Chuanchia labiosa and Schizopygopisis pylzovi captured in the Yellow River and Qaidam B grouping with G. eckloni from the Yellow River. Our research may helpful to further reconsideration of clearer taxonomy and improvement of biodiversity conservation strategy of Tibetan highland fishes.

Multiple convergent events created a nominal widespread species: Triplophysa stoliczkae (Steindachner, 1866) (Cobitoidea: Nemacheilidae).

BACKGROUND: Triplophysa stoliczkae is the most widespread species in the genus Triplophysa and may have originated from morphological convergence. To understand the evolutionary history of T. stoliczkae, we employed a multilocus approach to investigate the phylogenetics and the morphological evolution of T. stoliczkae on the Qinghai-Tibetan Plateau. RESULTS: All phylogenetic analyses (two mitochondrial and five nuclear loci), a genealogical sorting index and species tree inferences suggested that T. stoliczkae consists of distinct lineages that were not closest relatives. The time estimation indicated that the divergence events between "T. stoliczkae" and other Triplophysa species occurred from approximately 0.10 to 4.51 Ma. The ancestral state analyses supported the independent evolution of T. stoliczkae morphology in distinct lineages. The morphometric analysis and convergence estimates demonstrated significant phenotypic convergence among "T. stoliczkae" lineages. CONCLUSIONS: Triplophysa stoliczkae includes 4 different lineages with similar morphologies. The increasingly harsh environments that have occurred since the Pliocene have driven the occurrences of scrape-feeding fish in the genus Triplophysa. Morphological adaptations associated with scrape-feeding behavior resulted in convergences and the artificial lumping of four different species in the nominal taxon T. stoliczkae. A taxonomic revision for T. stoliczkae is needed.

Convergent evolution misled taxonomy in schizothoracine fishes (Cypriniformes: Cyprinidae).

Highly specialized grade (HSG; genera Gymnocypris, Oxygymnocypris, Schizopygopsis, Platypharodon and Chuanchia) of the Schizothoracinae (Cypriniformes: Cyprinidae) are endemic to the Qinghai-Tibet Plateau (QTP). Previously, two distinct ecomorphs were recognized according to trophic traits. One was a limnetic omnivore with normal lower jaw morphology, terminal mouth, and moderate or dense gill rakers, mostly inhabiting in open water of lakes, including Gymnocypris and Oxygymnocypris. Another was a benthic feeder with inferior mouth, sparse gill rakers and sharp horny sheath on the lower jaw for scraping of attached prey off hard substrates, including Schizopygopsis, Platypharodon and Chuanchia. However, traditional taxonomy of HSG based on these trophic traits presented extensive conflicts with the molecular studies in recent years. The possible cause could be convergent evolution in morphology, retention of ancestral polymorphisms or mitochondrial introgression, but these hypotheses could not be assessed due to incomplete taxon sampling and only mitochondrial data employed in previous works. Here, we conducted the most comprehensive molecular analysis on HSG fishes to date, using four mitochondrial loci and 152,464 genome-wide SNPs, and including 21 of 24 putative species and one undescribed Schizopygopsis species. Both SNP and mtDNA trees confirmed extensive paraphyly of genera Gymnocypris and Schizopygopsis, where species often were clustered together by watershed instead of by genus. Basal split into the north clade B and the south clade C (ca. 3.03 Ma) approximately by the Tanggula-Tanitawen Mountains in SLAF tree coincided with a violent uplift of the QTP during the phase A of 'Qingzang movement' (ca. 3.6 Ma). Ancestral state reconstruction of the trophic ecomorph showed that the limnetic omnivore ecomorph had evolved repeatedly in clade B and C. Furthermore, we presented a striking case of convergent evolution between two 'subspecies' Gymnocypris chui chui and G. chui longimandibularis, which had diverged as early as two million years ago (ca. 2.42 Ma). Ecological analyses revealed that similar food utilization, particularly in zooplankton, was the main underlying driving force. This work showed an example of taxonomy with the most extensive errors at the genus/species levels due to convergent evolution and suggested that trophic traits could be misleading in fish taxonomy. Therefore, we propose a major generic revision for HSG species.

Molecular systematics of the Triplophysa robusta (Cobitoidea) complex: Extensive gene flow in a depauperate lineage.

Gene flow between populations assumed to be isolated frequently leads to incorrect inferences of evolutionary history. Understanding gene flow and its causes has long been a key topic in evolutionary biology. In this study, we explored the evolutionary history of the Triplophysa robusta complex, using a combination of multilocus analyses and coalescent simulation. Our multilocus approach detected conspicuous mitonuclear discordances in the T. robusta complex. Mitochondrial results showed reticular clades, whereas the nuclear results corresponded with the morphological data. Coalescent simulation indicated that gene flow was the source of these discordances. Molecular clock analysis combined with geological processes suggest that intense geological upheavals have shaped a complicated evolutionary history for the T. robusta complex since the late Miocene, causing extensive gene flow which has distorted the molecular systematics of the T. robusta complex. We suggest that frequent gene flow may restrict speciation in the T. robusta complex, leading to such a depauperate lineage. Based on this comprehensive understanding, we provide our proposals for taxonomic revision of the T. robusta complex.

Adaptive Evolution of the Eda Gene and Scales Loss in Schizothoracine Fishes in Response to Uplift of the Tibetan Plateau.

Schizothoracine is the predominant wild fish subfamily of the Tibetan plateau (TP). Their scales, pharyngeal teeth and barbels have gradually regressed with increasing altitude. Schizothoracine have been divided into three groups: primitive, specialized and highly specialized. Ectodysplasin-A (Eda) has been considered as a major gene that contributes to the development of skin appendages. The present study cloned the Eda genes of 51 Schizothoracine fish species which represent the three groups and five Barbinae species. Phylogenetic analyses indicated that Eda may have acted as the genetic trigger for scale loss in the Schizothoracine. Furthermore, 14 single nucleotide polymorphisms (SNPs) and two deletions (18 bp and 6 bp in size), were also detected in the Eda coding sequence of the highly specialized group compared to the primitive group. The same SNPs and two indels result in four non-synonymous and two G-X-Y and 1 XY motif indels, which possibly contribute to significant structure changes in the Eda gene. The domain including (G-X-Y)n motif in the Eda gene is relatively conserved amongst teleosts. Based on the above results, we hypothesize that the evolution of Eda gene might be associated with the scale loss in Schizothoracine fishes in response to the phased uplift of the TP.

Transcriptomic profiling of Tibetan highland fish (Gymnocypris przewalskii) in response to the infection of parasite ciliate Ichthyophthirius multifiliis.

Gymnocypris przewalskii is a native cyprinid in the Lake Qinghai of the Qinghai-Tibetan Plateau. G. przewalskii is highly susceptible to the infection of a parasite, Ichthyophthirius multifiliis, in the artificial propagation and breeding. To better understand the host immune reaction to I. multifiliis infection, we characterize the gene expression profiles in the spleen of healthy and I. multifiliis infected G. przewalskii by RNA-seq. Totally, the transcriptomic analysis produces 463,031,110 high quality reads, which are assembled to 213,538 genes with N50 of 1918 bp and the average length of 1205 bp. Of assembled genes, 90.52% are annotated by public databases. The expression analysis shows 744 genes are significantly changed by the infection of I. multifiliis, which are validated by qRT-PCR with the correlation coefficient of 0.896. The differentially expressed genes are classified into 689 GO terms and 230 KEGG pathways, highlighting the promoted innate immunity in I. multifiliis infected G. przewalskii at 2 days post infection. Our results pinpoint that the up-regulated genes are enriched in TLR signaling pathway, inflammatory response and activation of immune cell migration. On the contrary, complement genes are down-regulated, indicating the evasion of host complement cascades by I. multifiliis. The repressed genes are also enriched in the pathways related to metabolism and endocrine, suggesting the metabolic disturbance in I. multifiliis treated G. przewalskii. In summary, the present study profiles the gene expression signature of G. przewalskii in the responses to I. multifiliis infection, and improves our understanding on molecular mechanisms of host-parasite interaction in G. przewalskii, which focuses the crucial function of TLRs, cytokines and complement components in the host defense against I. multifiliis.

Positive Darwinian selection within interferon regulatory factor genes of Gymnocypris przewalskii (Cyprinidae) on the Tibetan Plateau.

Tibetan Plateau (TP) had experienced phased uplift, resulting in inhospitable environment of low temperature, hypoxia and high ultraviolet radiation for Tibetan wildlife. Many organisms can well adapt to TP, it is of ecological and evolutionary interest to untangle how organisms adapt to extreme environment on TP through evolution. Previous studies mainly focused on hypoxia and metabolism related genes, but we know little about the evolutionary history of immune genes in Tibetan wildlife. In this study, we first identified 10 interferon regulatory factor (IRF) genes from Tibetan naked carp Gymnocypris przewalskii. Within this gene family, IRF3, IRF5, IRF7 and IRF8 contained positive selection sites. Evidences indicated that positive selection may lead to IRF genes functional alternations, presumably driving genes towards adaptation to the environmental changes. Taken together, our results suggested 4 candidate genes as interesting targets for further experimental confirmation of their functional variations and contributions to high altitude adaptation in Tibet fish.

High-fat diet aggravates the severity of the in vitro posttraumatic osteoarthritis model through macrophagic FBW7.

Osteoarthritis (OA) is a prevalent and distressing chronic degenerative joint disorder characterized by damaged articular cartilage and inflamed joints. Among risk factors, obesity has emerged as the second-leading contributor to OA after age. Obesity is believed to play a key role in the development and progression of OA. This study aimed to investigate the role and underlying mechanisms of high-fat diet (HFD)-induced obesity in the development of OA. Our findings revealed that HFD could aggravate the destabilization of the medial meniscus (DMM)-induced damage in the mouse model of obesity. Similar results were observed when macrophages obtained from HFD-fed mice were cocultured with cartilage and subsequently stimulated with interleukin-1ß (IL-1ß). Mechanistically, we observed a decrease in the expression of intraarticular macrophagic FBW7, which was implicated in the aggravation of OA in the HFD-fed animal. Furthermore, by modulating the immune status of macrophages, we found that reversing the macrophagic expression of FBW7 in these cells can alleviate the chondrocyte damage. In conclusion, this study provides novel insights into the pathological mechanisms underlying HFD-related OA development by identifying the role of FBW7 in synovial macrophages. These findings open up new avenues for research and therapeutic interventions targeting HFD-related OA.

Genetic adaptations of sea anemone to hydrothermal environment.

Hydrothermal vent habitats are characterized by high hydrostatic pressure, darkness, and the continuous release of toxic metal ions into the surrounding environment where sea anemones and other invertebrates thrive. Nevertheless, the understanding of metazoan metal ion tolerances and environmental adaptations remains limited. We assembled a chromosome-level genome for the vent sea anemone, Alvinactis idsseensis sp. nov. Comparative genomic analyses revealed gene family expansions and gene innovations in A. idsseensis sp. nov. as a response to high concentrations of metal ions. Impressively, the metal tolerance proteins MTPs is a unique evolutionary response to the high concentrations of Fe2+ and Mn2+ present in the environments of these anemones. We also found genes associated with high concentrations of polyunsaturated fatty acids that may respond to high hydrostatic pressure and found sensory and circadian rhythm-regulated genes that were essential for adaptations to darkness. Overall, our results provide insights into metazoan adaptation to metal ions, high pressure, and darkness in hydrothermal vents.

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.

Distinct and shared endothermic strategies in the heat producing tissues of tuna and other teleosts.

Although most fishes are ectothermic, some, including tuna and billfish, achieve endothermy through specialized heat producing tissues that are modified muscles. How these heat producing tissues evolved, and whether they share convergent molecular mechanisms, remain unresolved. Here, we generated a high-quality genome from the mackerel tuna (Euthynnus affinis) and investigated the heat producing tissues of this fish by single-nucleus and bulk RNA sequencing. Compared with other teleosts, tuna-specific genetic variation is strongly associated with muscle differentiation. Single-nucleus RNA-seq revealed a high proportion of specific slow skeletal muscle cell subtypes in the heat producing tissues of tuna. Marker genes of this cell subtype are associated with the relative sliding of actin and myosin, suggesting that tuna endothermy is mainly based on shivering thermogenesis. In contrast, cross-species transcriptome analysis indicated that endothermy in billfish relies mainly on non-shivering thermogenesis. Nevertheless, the heat producing tissues of the different species do share some tissue-specific genes, including vascular-related and mitochondrial genes. Overall, although tunas and billfishes differ in their thermogenic strategies, they share similar expression patterns in some respects, highlighting the complexity of convergent evolution.

Genomes of Two Flying Squid Species Provide Novel Insights into Adaptations of Cephalopods to Pelagic Life.

Pelagic cephalopods have evolved a series of fascinating traits, such as excellent visual acuity, high-speed agility, and photophores for adaptation to open pelagic oceans. However, the genetic mechanisms underpinning these traits are not well understood. Thus, in this study, we obtained high-quality genomes of two purpleback flying squid species (Sthenoteuthis oualaniensis and Sthenoteuthis sp.), with sizes of 5450 Mb and 5651 Mb, respectively. Comparative genomic analyses revealed that the S-crystallin subfamily SL20-1 associated with visual acuity in the purpleback flying squid lineage was significantly expanded, and the evolution of high-speed agility for the species was accompanied by significant positive selection pressure on genes related to energy metabolism. These molecular signals might have contributed to the evolution of their adaptative predatory and anti-predatory traits. In addition, the transcriptomic analysis provided clear indications of the evolution of the photophores of purpleback flying squids, especially the recruitment of new genes and energy metabolism-related genes which may have played key functional roles in the process.

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.

Comparative genome anatomy reveals evolutionary insights into a unique amphitriploid fish.

Triploids are rare in nature because of difficulties in meiotic and gametogenic processes, especially in vertebrates. The Carassius complex of cyprinid teleosts contains sexual tetraploid crucian carp/goldfish (C. auratus) and unisexual hexaploid gibel carp/Prussian carp (C. gibelio) lineages, providing a valuable model for studying the evolution and maintenance mechanism of unisexual polyploids in vertebrates. Here we sequence the genomes of the two species and assemble their haplotypes, which contain two subgenomes (A and B), to the chromosome level. Sequencing coverage analysis reveals that C. gibelio is an amphitriploid (AAABBB) with two triploid sets of chromosomes; each set is derived from a different ancestor. Resequencing data from different strains of C. gibelio show that unisexual reproduction has been maintained for over 0.82 million years. Comparative genomics show intensive expansion and alterations of meiotic cell cycle-related genes and an oocyte-specific histone variant. Cytological assays indicate that C. gibelio produces unreduced oocytes by an alternative ameiotic pathway; however, sporadic homologous recombination and a high rate of gene conversion also exist in C. gibelio. These genomic changes might have facilitated purging deleterious mutations and maintaining genome stability in this unisexual amphitriploid fish. Overall, the current results provide novel insights into the evolutionary mechanisms of the reproductive success in unisexual polyploid vertebrates.

The First Draft Genome of a Cold-Water Coral Trachythela sp. (Alcyonacea: Stolonifera: Clavulariidae).

Cold-water corals (CWCs) are important habitats for creatures in the deep-sea environment, but they have been degraded by anthropogenic activity. So far, no genome for any CWC has been reported. Here, we report a draft genome of Trachythela sp., which represents the first genome of CWCs to date. In total, 56 and 65 Gb of raw reads were generated from Illumina and Nanopore sequencing platforms, respectively. The final assembled genome was 578.26 Mb, which consisted of 396 contigs with a contig N50 of 3.56 Mb, and the genome captured 90.1% of the metazoan Benchmarking Universal Single-Copy Orthologs. We identified 335 Mb (57.88% of the genome) of repetitive elements, which is a higher proportion compared with others in the Cnidarians, along with 35,305 protein-coding genes. We also detected 483 expanded and 51 contracted gene families, and many of them were associated with longevity, ion transposase, heme-binding nicotinamide adenine dinucleotide, and metabolic regulators of transcription. Overall, we believe this genome will serve as an important resource for studies on community protection for CWCs.

A new species of Triplophysa (Cypriniformes, Nemacheilidae) from Weihe River in Gansu Province, China.

A new species of Tibetan loach, Triplophysa weiheensis sp. nov., is described from the Weihe River in Gansu Province, China, based on morphological and molecular analyses. The new species can be distinguished from all known congeners by a unique combination of the following characters: scaleless; snout abruptly sloping downward, anterior to anterior nostril; lower jaw crescentic, not sharp; body without obvious mottling; lateral line interrupted on posterior trunk at pelvic-fin distal extremity; caudal-peduncle length 2.0-2.7 times its depth; branched rays of pectoral fin 10-11; branched rays of pelvic fin 5-6; inner gill rakers on 1 st gill arch 14-16; vertebrae 4+34-36; intestine with 6-7 loops, length ca. 1.8 times SL ( n=3); bony capsule of air bladder small and thin; posterior chamber of air bladder absent.