The natural sex reversal severely affects the sex ratio and thus decreases the productivity of the rice field eel (Monopterus albus). How to understand and manipulate this process is one of the major issues for the rice field eel stocking. So far the genomics and transcriptomics data available for this species are still scarce. Here we provide a comprehensive study of transcriptomes of brain and gonad tissue in three sex stages (female, intersex and male) from the rice field eel to investigate changes in transcriptional level during the sex reversal process. Approximately 195 thousand unigenes were generated and over 44.4 thousand were functionally annotated. Comparative study between stages provided multiple differentially expressed genes in brain and gonad tissue. Overall 4668 genes were found to be of unequal abundance between gonad tissues, far more than that of the brain tissues (59 genes). These genes were enriched in several different signaling pathways. A number of 231 genes were found with different levels in gonad in each stage, with several reproduction-related genes included. A total of 19 candidate genes that could be most related to sex reversal were screened out, part of these genes’ expression patterns were validated by RT-qPCR. The expression of spef2, maats1, spag6 and dmc1 were abundant in testis, but was barely detected in females, while the 17β-hsd12, zpsbp3, gal3 and foxn5 were only expressed in ovary. This study investigated the complexity of brain and gonad transcriptomes in three sex stages of the rice field eel. Integrated analysis of different gene expression and changes in signaling pathways, such as PI3K-Akt pathway, provided crucial data for further study of sex transformation mechanisms.

MicroRNAs (miRNAs) are endogenous small RNAs that can regulate target mRNAs by binding to their sequences in the 3′ untranslated region. The expression of miRNAs and their biogenetic pathway are involved in sexual differentiation and in the regulation of the development of germ cells and gonadal somatic cells. The rice field eel (Monopterus albus) undergoes a natural sexual transformation from female to male via an intersex stage during its life cycle. To investigate the molecular mechanisms of this sexual transformation, miRNAs present in the different sexual stages of the rice field eel were identified by high-throughput sequencing technology. A significantly differential expression among the 3 genders (p < 0.001) was observed for 48 unique miRNAs and 3 miRNAs*. Only 9 unique miRNAs showed a more than 8-fold change in their expression among the 3 genders, including mal-miR-430a and mal-miR-430c which were higher in females than in males. However, mal-miR-430b was only detected in males. Several potential miRNA target genes (cyp19a, cyp19b, nr5a1b, foxl2 amh, and vasa) were also investigated. Real-time RT-PCR demonstrated highly specific expression patterns of these genes in the 3 genders of the rice field eel. Many of these genes are targets of mal-miR-430b according to the TargetScan and miRTarBase. These results suggest that the miR-430 family may be involved in the sexual transformation of the rice field eel.

Hypoxia‐inducible factor (HIF) is a crucial regulator of cellular and systemic responses to low oxygen levels. Here we firstly cloned three HIF‐α isoforms from the basal branches of Osteichthyes and used computational tools to characterise the molecular change underlying the functional divergence of HIF‐α isoforms in different lineages. Only the HIF‐1α and HIF‐2α in African lungfish and amphibians were found under positive selection. HIF‐1α and ‐2α were less functionally divergent in basal ray‐finned fish than in teleosts, and showed conserved but different transcriptional activity towards specific target genes.

Hypoxia-inducible factor 1 (HIF-1) functions as a master regulator of the cellular response to hypoxic stress. Two HIF-1α paralogs, HIF-1αA and HIF-1αB, were generated in euteleosts by the specific, third round of genome duplication, but one paralog was later lost in most families with the exception of cyprinid fish. How these duplicates function in mitochondrial regulation and whether their preservation contributes to the hypoxia tolerance demonstrated by cyprinid fish in freshwater environments is not clear. Here we demonstrated the divergent function of these two zebrafish Hif-1a paralogs through cellular approaches. The results showed that Hif-1aa played a role in tricarboxylic acid cycle by increasing the expression of Citrate synthase and the activity of mitochondrial complex II, and it also enhanced mitochondrial membrane potential and ROS production by reducing free Ca²⁺ in the cytosol. Hif-1ab promoted intracellular ATP content by up-regulating the activity of mitochondrial complexes I, III and IV and the expression of related genes. Furthermore, both the two zebrafish Hif-1a paralogs promoted mitochondrial mass and the expression level of mtDNA, contributing to mitochondrial biogenesis. Our study reveals the divergent functions of Hif-1aa and Hif-1ab in cellular mitochondrial regulation.

Molecular adaptation to the severe environments present during the uplift of the Qinghai-Tibet Plateau has attracted the attention of researchers. The divergence of the three specialization groups of schizothoracins (Primitive, Specialized and Highly Specialized) may correspond to the three phases of plateau uplift. Based on the transcripts of representative species of the three specialized groups and an outgroup, genes in schizothoracins that may have played important roles during the adaptation to new environments were investigated. The contigs of Gymnodiptychus dybowskii and Schizothorax pseudaksaiensis were compared with those of Gymnocypris przewalskii ganzihonensis and the outgroup Sinocyclocheilus angustiporus, and 5,894 ortholog groups with an alignment length longer than 90 nt after deleting gaps were retained. Evolutionary analyses indicated that the average evolutionary rate of the branch leading to the Specialized group was faster than that of the branch leading to the Highly Specialized group. Moreover, the numbers of gene categories in which more than half of the genes evolved faster than the average values of the genome were 117 and 15 along the branches leading to the Specialized and Highly Specialized groups, respectively. A total of 40, 36, and 55 genes were likely subject to positive selection along the branches leading to the Primitive, Specialized and Highly Specialized groups, respectively, and many of these genes are likely relevant to adaptation to the cold temperatures, low oxygen concentrations, and strong ultraviolet radiation that result from elevation. By selecting representative species of the three groups of schizothoracins and applying next-generation sequencing technology, several candidate genes corresponding to adaptation to the three phases of plateau uplift were identified. Some of the genes identified in this report that were likely subject to positive selection are good candidates for subsequent evolutionary and functional analyses of adaptation to high altitude.