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.

Self-powered soft robot in the Mariana Trench.

The deep sea remains the largest unknown territory on Earth because it is so difficult to explore1-4. Owing to the extremely high pressure in the deep sea, rigid vessels5-7 and pressure-compensation systems8-10 are typically required to protect mechatronic systems. However, deep-sea creatures that lack bulky or heavy pressure-tolerant systems can thrive at extreme depths11-17. Here, inspired by the structure of a deep-sea snailfish15, we develop an untethered soft robot for deep-sea exploration, with onboard power, control and actuation protected from pressure by integrating electronics in a silicone matrix. This self-powered robot eliminates the requirement for any rigid vessel. To reduce shear stress at the interfaces between electronic components, we decentralize the electronics by increasing the distance between components or separating them from the printed circuit board. Careful design of the dielectric elastomer material used for the robot's flapping fins allowed the robot to be actuated successfully in a field test in the Mariana Trench down to a depth of 10,900 metres and to swim freely in the South China Sea at a depth of 3,224 metres. We validate the pressure resilience of the electronic components and soft actuators through systematic experiments and theoretical analyses. Our work highlights the potential of designing soft, lightweight devices for use in extreme conditions.

Genomic insights into the seawater adaptation in Cyprinidae.

BACKGROUND: Cyprinidae, the largest fish family, encompasses approximately 367 genera and 3006 species. While they exhibit remarkable adaptability to diverse aquatic environments, it is exceptionally rare to find them in seawater, with the Far Eastern daces being of few exceptions. Therefore, the Far Eastern daces serve as a valuable model for studying the genetic mechanisms underlying seawater adaptation in Cyprinidae. RESULTS: Here, we sequenced the chromosome-level genomes of two Far Eastern daces (Pseudaspius brandtii and P. hakonensis), the two known cyprinid fishes found in seawater, and performed comparative genomic analyses to investigate their genetic mechanism of seawater adaptation. Demographic history reconstruction of the two species reveals that their population dynamics are correlated with the glacial-interglacial cycles and sea level changes. Genomic analyses identified Pseudaspius-specific genetic innovations related to seawater adaptation, including positively selected genes, rapidly evolving genes, and conserved non-coding elements (CNEs). Functional assays of Pseudaspius-specific variants of the prolactin (prl) gene showed enhanced cell adaptation to greater osmolarity. Functional assays of Pseudaspius specific CNEs near atg7 and usp45 genes suggest that they exhibit higher promoter activity and significantly induced at high osmolarity. CONCLUSIONS: Our results reveal the genome-wide evidence for the evolutionary adaptation of cyprinid fishes to seawater, offering valuable insights into the molecular mechanisms supporting the survival of migratory fish in marine environments. These findings are significant as they contribute to our understanding of how cyprinid fishes navigate and thrive in diverse aquatic habitats, providing useful implications for the conservation and management of marine ecosystems.

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.

Proteomics and metabolomics profiling reveal panels of circulating diagnostic biomarkers and molecular subtypes in stable COPD.

BACKGROUND: Chronic obstructive pulmonary disease (COPD) is a complex and heterogeneous disease with high morbidity and mortality, especially in advanced patients. We aimed to develop multi-omics panels of biomarkers for the diagnosis and explore its molecular subtypes. METHODS: A total of 40 stable patients with advanced COPD and 40 controls were enrolled in the study. Proteomics and metabolomics techniques were applied to identify potential biomarkers. An additional 29 COPD and 31 controls were enrolled for validation of the obtained proteomic signatures. Information on demographic, clinical manifestation, and blood test were collected. The ROC analyses were carried out to evaluate the diagnostic performance, and experimentally validated the final biomarkers on mild-to-moderate COPD. Next, molecular subtyping was performed using proteomics data. RESULTS: Theophylline, palmitoylethanolamide, hypoxanthine, and cadherin 5 (CDH5) could effectively diagnose advanced COPD with high accuracy (auROC = 0.98, sensitivity of 0.94, and specificity of 0.95). The performance of the diagnostic panel was superior to that of other single/combined results and blood tests. Proteome based stratification of COPD revealed three subtypes (I-III) related to different clinical outcomes and molecular feature: simplex COPD, COPD co-existing with bronchiectasis, and COPD largely co-existing with metabolic syndrome, respectively. Two discriminant models were established using the auROC of 0.96 (Principal Component Analysis, PCA) and 0.95 (the combination of RRM1 + SUPV3L1 + KRT78) in differentiating COPD and COPD with co-morbidities. Theophylline and CDH5 were exclusively elevated in advanced COPD but not in its mild form. CONCLUSIONS: This integrative multi-omics analysis provides a more comprehensive understanding of the molecular landscape of advanced COPD, which may suggest molecular targets for specialized therapy.

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.

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.

Comparative mitochondrial genome analysis of Neodontobutis hainanensis and Perccottus glenii reveals conserved genome organization and phylogeny.

To investigate the molecular evolution of mitochondrial genomes among the family Odontobutidae, the complete mitochondrial genomes of Neodontobutis hainanensis and Perccottus glenii were sequenced and compared with seven odontobutids species. The genome organization, base composition, codon usage, and gene arrangement of N. hainanensis exhibited high similarity to P. glenii compared to those of other Odontobutidae species. Reconstructed phylogenetic analyses of Odontobutidae strongly supported that Neodontobutis and Perccottus formed a unifying group sister to Odontobutis. Our molecular dating time revealed that the two species diverged approximately 21.7 Ma during Miocene, later than that of Odontobutis. Selection analyses showed stronger selective constraints in mitochondrial genes for P. glenii. However, two positively selected sites in NADH4 and NADH6 genes were respectively detected in N. hainanensis and P. glenii, indicating that they might evolve different metabolic performance in response to the contrasting environments.

FishDB: an integrated functional genomics database for fishes.

BACKGROUND: Hundreds of genomes and transcriptomes of fish species have been sequenced in recent years. However, fish scholarship currently lacks a comprehensive, integrated, and up-to-date collection of fish genomic data. RESULTS: Here we present FishDB, the first database for fish multi-level omics data, available online at http://fishdb.ihb.ac.cn . The database contains 233 fish genomes, 201 fish transcriptomes, 5841 fish mitochondrial genomes, 88 fish gene sets, 16,239 miRNAs of 65 fishes, 1,330,692 piRNAs and 4852 lncRNAs of Danio rerio, 59,040 Mb untranslated regions (UTR) of 230 fishes, and 31,918 Mb coding sequences (CDS) of 230 fishes. Among these, we newly generated a total of 11 fish genomes and 53 fish transcriptomes. CONCLUSIONS: This release contains over 410,721.67 Mb sequences and provides search functionality, a BLAST server, JBrowse, and PrimerServer modules.

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.

Evidence of high-altitude adaptation in the glyptosternoid fish, Creteuchiloglanis macropterus from the Nujiang River obtained through transcriptome analysis.

BACKGROUND: Organisms living at high altitudes face low oxygen and temperature conditions; thus, the genetic mechanisms underlying the adaptations in these organisms merit investigation. The glyptosternoid fish, Creteuchiloglanis macropterus mainly inhabits regions with gradual increases in altitudes along the Nujiang River and might serve as an appropriate evolutionary model for detecting adaptation processes in environments with altitude changes. RESULTS: We constructed eleven RNA-sequencing (RNA-seq) libraries of C. macropterus collected from five locations at different altitudes to identify the genetic signatures of high-altitude adaptation. The comparative genomic analysis indicated that C. macropterus has an accelerated evolutionary rate compared with that of fishes in the lowland, and fishes at higher altitudes might evolve faster. Functional enrichment analysis of the fast-evolving and positively selected genes, differentially expressed genes and highly expressed genes, showed that these genes were involved in many functions related to energy metabolism and hypoxia. CONCLUSIONS: Our study provides evidence of high-altitude adaptation in C. macropterus, and the detected adaptive genes might be a resource for future investigations of adaptations to high-altitude environments in other fishes.

Phylogeographic structure, cryptic speciation and demographic history of the sharpbelly (Hemiculter leucisculus), a freshwater habitat generalist from southern China.

BACKGROUND: Species with broad distributions frequently divide into multiple genetic forms and may therefore be viewed as "cryptic species". Here, we used the mitochondrial cytochrome b (Cytb) and 12 nuclear DNA loci to investigate phylogeographic structures of the sharpbelly (Hemiculter leucisculus) in rivers in southern China and explored how the geological and climatic factors have shaped the genetic diversity and evolutionary history of this species. RESULTS: Our mitochondrial phylogenetic analysis identified three major lineages (lineages A, B, and C). Lineages B and C showed a relatively narrower geographic distribution, whereas lineage A was widely distributed in numerous drainages. Divergence dates suggested that H. leucisculus populations diverged between 1.61-2.38 Ma. Bayesian species delimitation analysis using 12 nuclear DNA loci indicated the three lineages probably represented three valid taxa. Isolation-with-migration (IM) analysis found substantial gene flow has occurred among the three lineages. Demographic analyses showed that lineages B and C have experienced rapid demographic expansion at 0.03 Ma and 0.08 Ma, respectively. CONCLUSIONS: Hemiculter leucisculus populations in drainages in southern China comprise three mtDNA lineages, and each of which may represent a separate species. Intense uplift of the Qinghai-Tibetan Plateau, evolution of Asian monsoons, changes in paleo-drainages, and poor dispersal ability may have driven the divergence of the three putative species. However, gene flow occurs among the three lineages. Climatic fluctuations have a prominent impact on the populations from the lineages B and C, but exerted little influence on the lineage A.

Correlated expression of retrocopies and parental genes in zebrafish.

Previous studies of the function and evolution of retrocopies in plants, Drosophila and non-mammalian chordates provided new insights into the origin of novel genes. However, little is known about retrocopies and their parental genes in teleosts, and it remains obscure whether there is any correlation between them. The present study aimed to characterize the spatial and temporal expression profiles of retrogenes and their parental genes based on RNA-Seq data from Danio rerio embryos and tissues from adult. Using a modified pipeline, 306 retrocopies were identified in the zebrafish genome, most of which exhibited ancient retroposition, and 76 of these showed a Ks < 2.0. Expression of a retrocopy is generally expected to present no correlation with its parental gene, as regulatory regions are not part of the retroposition event. Here, this assumption was tested based on RNA-Seq data from eight stages and thirteen tissue types of zebrafish. However, the result suggested that retrocopies displayed correlated expression with their parental genes. The level of correlation was found to decrease during embryogenesis, but to increase slightly within a tissue using Ks as the proxy for the divergence time. Tissue specificity was also observed: retrocopies were found to be expressed at a more specific level compared with their parental genes. Unlike Drosophila, which has sex chromosomes, zebrafish do not show testis-biased expression. Our study elaborated temporal and spatial patterns of expression of retrocopies in zebrafish, examined the correlation between retrocopies and parental genes and analyzed potential source of regulated elements of retrocopies, which lay a foundation for further functional study of retrocopies.

Genome-wide identification of lineage-specific genes within Caenorhabditis elegans.

With the rapid growth of sequencing technology, a number of genomes and transcriptomes of various species have been sequenced, contributing to the study of lineage-specific genes (LSGs). We identified two sets of LSGs using BLAST: one included Caenorhabditis elegans species-specific genes (1423, SSGs), and the other consisted of Caenorhabditis genus-specific genes (4539, GSGs). The subsequent characterization and analysis of the SSGs and GSGs showed that they have significant differences in evolution and that most LSGs were generated by gene duplication and integration of transposable elements (TEs). We then performed temporal expression profiling and protein function prediction and observed that many SSGs and GSGs are expressed and that genes involved with sex determination, specific stress, immune response, and morphogenesis are over-represented, suggesting that these specific genes may be related to the Caenorhabditis nematodes' special ability to survive in severe and extreme environments.

[Activity Study of Icarlin and Epimedin C Monomer under the Same Molarity in Osteoporosis Ze- brafish Model].

Objective To evaluate anti-osteoporotic activity of icariin and Epimedin C monomer under the same molarity in predinsolone-induced osteoporosis zebrafish. Methods Zebrafish larvae after 4-day fertilization were divided into group S [0. 5% dimethyl sulfoxide (DMSO) , A (25 µmol/L prednisolone, 0. 5% DMSO), B (2 IU/L salmon calcitonin, 25 µmol/L prednisolone,0. 5% DMSO), C (1. 5 1,mol/L icariin, 25 µmol/L prednisolone, 0. 5% DMSO) , D (15 µLmol/L icariin,25 µmol/L prednisolone, 0. 5% DM- SO), E (150 µmol/L icariin, 25 µmol/L prednisolone, 0. 5% DMSO), F (1. 5 µmol/L Epimediri C, 25 µmol/L prednisolone, 0. 5% DMSO) , G (15 µmol/L Epimedin C, 25 µmol/L prednisolone, 0.5% DM- SO) , H (150 µmol/L Epimedin C, 25 µmol/L prednisolone, 0. 5% DMSO). All culture solution contained 0. 5% DMSO. All the young fishes were grown in a 24-well plate. The culture medium was changed every day. They were cultured in a incubator box at 28. 5 °C and killed at day 9. Zebrafish skeleton was stained with alizarin red. The stained Zebrafish ventral skull was observed using microscope, and mineralized area was quantitatively analyzed. Results Compared with group S, accumulative integrated optical densi- ty(IOD)of the mineralized area significantly decreased in group A (P <0. 01) ; accumulative IOD of the mineralized area significantly increased in group B (P <0. 01). The accumulative IOD of the mineralized area showed weakly increasing tendency in group C, D, and E along with increased concentration (P < 0. 05). Compared with group A, accumulative IOD obviously increased in group B with statistical difference (P <0. 01) , but with no statistical difference as compared with group C or group D (P >0. 05). Statistical difference existed in accumulative IOD between group A and group E (P <0. 05). The mineralized area showed increasing tendency in group F and group G along with increased concentration (P <0. 05), and accumulative IOD obviously increased as well (P <0. 05). No Zebrafish embryo survived in group H. There was no statistical difference in Zebrafish embryo survival among group E, F, or G (P >0. 05). The staining of Zebrafish skull was clearly seen in group S, with vertebrae and bilateral branchial skeleton clearly seen. The intensity of staining in the same area was obviously attenuated in group A. The osteo- genesis was speeded up under the same condition in group B, with obviously enlarged mineralized area and more darkly stained bone tissue. The mineralization of skull was gradually increasing during the stai- ning process in group C, D, E, F, and G. The mineralized area and the intensity of staining were gradually enhanced, and changes of vertebrae were most obviously seen in group C, D, E, F, and G, but they were not arrived at the stained intensity level in group B. Conclusions Osteoporosis Zebrafish model is a simple and efficient model for screening bioactive ingredients of Chinese herbs. The activity of Epimedin C at low concentration was better than icariin in this model. But possible toxicity of Epimedin C at high concentration needs to be further studied.

Phylogenomic analyses reveal a molecular signature linked to subterranean adaptation in rodents.

BACKGROUND: Genome-wide signatures of convergent evolution are widely expected but rarely revealed in animals. Subterranean rodent genome and transcriptome data produced by next-generation sequencing facilitate the use of phylogenetic methods to infer non-synonymous and synonymous substitution rates within coding regions, which can reveal changes at the molecular level that are correlated with the dramatic shift from a terrestrial to subterranean habitat. RESULTS: Our study used previously sequenced genome or transcriptome data of two subterranean rodents, the blind mole rat and naked mole rat, and their terrestrial relatives, the mouse and guinea pig, to investigate the genetic basis of rodent subterranean adaptation. An analysis of 4996 orthologous genes revealed that the substitution pace of coding sequences was significantly slower in the blind mole rat than in the mouse, and slower in the naked mole rat than in the guinea pig. The dN/dS ratio was significantly higher in the blind mole rat than in the mouse and in the naked mole rat than in the guinea pig. These patterns are most likely related to the longer generation time and lower effective population size of subterranean rodents caused by subterranean ecological constraints. We also identified some genes and gene ontology (GO) categories that might be candidates for adaptation to subterranean life. CONCLUSIONS: Our study reveals a case of subterranean convergent evolution in rodents that is correlated with change in the pace and mode of molecular evolution observed at the genome scale. We believe that this genomic signature could have also evolved in other cases of subterranean convergence. Additionally, the genes that displayed the most radical changes in their patterns of evolution and their associated GO categories provide a strong basis for further comparative and functional studies, and potentially reveal molecular signatures of adaptation to subterranean life.

Biogeographic history and high-elevation adaptations inferred from the mitochondrial genome of Glyptosternoid fishes (Sisoridae, Siluriformes) from the southeastern Tibetan Plateau.

BACKGROUND: The distribution of the Chinese Glyptosternoid catfish is limited to the rivers of the Tibetan Plateau and peripheral regions, especially the drainage areas of southeastern Tibet. Therefore, Glyptosternoid fishes are ideal for reconstructing the geological history of the southeastern Tibet drainage patterns and mitochondrial genetic adaptions to high elevations. RESULTS: Our phylogenetic results support the monophyly of the Sisoridae and the Glyptosternoid fishes. The reconstructed ancestral geographical distribution suggests that the ancestral Glyptosternoids was widely distributed throughout the Brahmaputra drainage in the eastern Himalayas and Tibetan area during the Late Miocene (c. 5.5 Ma). We found that the Glyptosternoid fishes lineage had a higher ratio of nonsynonymous to synonymous substitutions than those found in non-Glyptosternoids. In addition, ωpss was estimated to be 10.73, which is significantly higher than 1 (p-value 0.0002), in COX1, which indicates positive selection in the common ancestral branch of Glyptosternoid fishes in China. We also found other signatures of positive selection in the branch of specialized species. These results imply mitochondrial genetic adaptation to high elevations in the Glyptosternoids. CONCLUSIONS: We reconstructed a possible scenario for the southeastern Tibetan drainage patterns based on the adaptive geographical distribution of the Chinese Glyptosternoids in this drainage. The Glyptosternoids may have experienced accelerated evolutionary rates in mitochondrial genes that were driven by positive selection to better adapt to the high-elevation environment of the Tibetan Plateau.

Evolutionary fate and implications of retrocopies in the African coelacanth genome.

BACKGROUND: The coelacanth is known as a "living fossil" because of its morphological resemblance to its fossil ancestors. Thus, it serves as a useful model that provides insight into the fish that first walked on land. Retrocopies are a type of novel genetic element that are likely to contribute to genome or phenotype innovations. Thus, investigating retrocopies in the coelacanth genome can determine the role of retrocopies in coelacanth genome innovations and perhaps even water-to-land adaptations. RESULTS: We determined the dS values, dN/dS ratios, expression patterns, and enrichment of functional categories for 472 retrocopies in the African coelacanth genome. Of the retrocopies, 85-355 were shown to be potentially functional (i.e., retrogenes). The distribution of retrocopies based on their dS values revealed a burst pattern of young retrocopies in the genome. The retrocopy birth pattern was shown to be more similar to that in tetrapods than ray-finned fish, which indicates a genomic transformation that accompanied vertebrate evolution from water to land. Among these retrocopies, retrogenes were more prevalent in old than young retrocopies, which indicates that most retrocopies may have been eliminated during evolution, even though some retrocopies survived, attained biological function as retrogenes, and became old. Transcriptome data revealed that many retrocopies showed a biased expression pattern in the testis, although the expression was not specifically associated with a particular retrocopy age range. We identified 225 Ensembl genes that overlapped with the coelacanth genome retrocopies. GO enrichment analysis revealed different overrepresented GO (gene ontology) terms between these "retrocopy-overlapped genes" and the retrocopy parent genes, which indicates potential genomic functional organization produced by retrotranspositions. Among the 225 retrocopy-overlapped genes, we also identified 46 that were coelacanth-specific, which could represent a potential molecular basis for coelacanth evolution. CONCLUSIONS: Our study identified 472 retrocopies in the coelacanth genome. Sequence analysis of these retrocopies and their parent genes, transcriptome data, and GO annotation information revealed novel insight about the potential role of genomic retrocopies in coelacanth evolution and vertebrate adaptations during the evolutionary transition from water to land.