Lampreys are extant representatives of the jawless vertebrate lineage that diverged from jawed vertebrates around 500 million years ago. Lamprey genomes contain information crucial for understanding the evolution of gene families in vertebrates. The ATP-binding cassette (ABC) gene family is found from prokaryotes to eukaryotes. The recent availability of two lamprey draft genomes from sea lamprey Petromyzon marinus and Japanese lamprey Lethenteron japonicum presents an opportunity to infer early evolutionary events of ABC genes in vertebrates. We conducted a genome-wide survey of the ABC gene family in two lamprey draft genomes. A total of 37 ABC transporters were identified and classified into seven subfamilies; namely seven ABCA genes, 10 ABCB genes, 10 ABCC genes, three ABCD genes, one ABCE gene, three ABCF genes, and three ABCG genes. The ABCA subfamily has expanded from three genes in sea squirts, seven and nine in lampreys and zebrafish, to 13 and 16 in human and mouse. Conversely, the multiple copies of ABCB1-, ABCG1-, and ABCG2-like genes found in sea squirts have contracted in the other species examined. ABCB2 and ABCB3 seem to be new additions in gnathostomes (not in sea squirts or lampreys), which coincides with the emergence of the gnathostome-specific adaptive immune system. All the genes in the ABCD, ABCE and ABCF subfamilies were conserved and had undergone limited duplication and loss events. In the sea lamprey transcriptomes, the ABCE and ABCF gene subfamilies were ubiquitously and highly expressed in all tissues while the members in other gene subfamilies were differentially expressed. Thirteen more lamprey ABC transporter genes were identified in this study compared with a previous study. By concatenating the same gene sequences from the two lampreys, more full length sequences were obtained, which significantly improved both the assignment of gene names and the phylogenetic trees compared with a previous analysis using partial sequences. The ABC gene subfamilies in chordates have undergone obvious expansion or contraction. The ABCA subfamily showed the highest gene expansion rate during chordate evolution. The evolution of ABC transporters in lampreys requires further evaluation because the present results are based on a draft genome.

Goldfish is an ornamental fish with diverse phenotypes. However, the limited genomic resources of goldfish hamper our understanding of the genetic basis for its phenotypic diversity. To provide enriched genomic resources and infer possible mechanisms underlying skin pigmentation, we performed a large-scale transcriptomic sequencing on 13 adult goldfish tissues, larvae at one- and three-days post hatch, and skin tissues with four different color pigmentation. A total of 25.52 Gb and 149.80 Gb clean data were obtained using the PacBio and Illumina platforms, respectively. Onto the goldfish reference genome, we mapped 137,674 non-redundant transcripts, of which 5.54% was known isoforms and 78.53% was novel isoforms of the reference genes, and the remaining 21,926 isoforms are novel isoforms of additional new genes. Both skin-specific and color-specific transcriptomic analyses showed that several significantly enriched genes were known to be involved in melanogenesis, tyrosine metabolism, PPAR signaling pathway, folate biosynthesis metabolism and so on. Thirteen differentially expressed genes across different color skins were associated with melanogenesis and pteridine synthesis including mitf, ednrb, mc1r, tyr, mlph and gch1, and xanthophore differentiation such as pax7, slc2a11 and slc2a15. These transcriptomic data revealed pathways involved in goldfish pigmentation and improved the gene annotation of the reference genome.

Accumulation of excessive ammonia is a big threat to aquatic animals, which causes adverse effects on the health, production reduction, and even high mortality. The razor clam Sinonovacula constricta, a bivalve living in intertidal mudflat with a deep-burrowing lifestyle, often faces a high concentration of ambient ammonia. However, there is less available information concerning the toxic effects of ammonia on razor clam and its molecular mechanisms of adaptation to ammonia stress. The aim of this study was to investigate the effects of ammonia exposure on the gill and hepatopancreas of razor clam by transcriptome sequencing. The results showed that the median lethal concentration of ammonia was 244.55 mg/L for 96 h. A total of 1,415 and 306 differentially expressed genes (DEGs) were identified in the gill and hepatopancreas, respectively. The functional annotation showed that DEGs of the gill were mainly involved in the regulation of nitrogen compound metabolic process, nitrogen compound transport, and amide transport. The DEGs of the hepatopancreas were mostly enriched in oxidation-reduction process, response to stress, and amine metabolic process. The expression levels of NH₃/NH₄⁺ transporting channels and H⁺ excreting-related genes, including Rhesus glycoproteins (Rh), Na⁺/K⁺-ATPase (NKA), Na⁺/H⁺ exchanger, V-ATPase (VHA), and carbonic anhydrase (CA), were upregulated significantly in the gill (p < 0.05). In addition, the expression levels of glutamine and urea synthesis-related genes that played vital roles in ammonia detoxification, such as glutamine synthetase (GS), arginase (ARG), and argininosuccinate synthetase (ASS), were also increased obviously in the hepatopancreas (p < 0.05). Taken together, our results indicate that the synergistic action of ammonia excretion in the gill and ammonia metabolism in the hepatopancreas might be the mechanism through which the clams tolerate to environmental ammonia. This study provides a molecular basis for the better evaluation of the responding mechanism of ammonia tolerance.

Nuclear progesterone receptor (nPR) is an evolutionary innovation in vertebrates that mediates genomic responses to progesterone. Vertebrates also respond to progesterone via membrane progesterone receptors (mPRs) or membrane associated progesterone receptors (MAPRs) through rapid nongenomic mechanisms. Lampreys are extant agnathan vertebrates, residing at the evolutionary juncture where vertebrates diverged from invertebrates. A survey of the progesterone receptor (PR) gene sequences in lamprey genomes would inform PR gene evolutionary events during the transition from invertebrates to vertebrates. In this study, we annotated sequences of one nPR, four mPR (β, γ, δ and ε) and four MAPR genes from genomes of two lamprey species (Petromyzon marinus and Lethenteron japonicum). To infer the origin and evolutionary history of PR genes, we constructed phylogenetic trees of PR homologous sequences across representative species of metazoans. Phylogenetic analyses revealed that the mPRγ gene first appeared in non-bilaterians, and the mPRβ gene likely arose from a duplication of mPRγ. On the other hand, the mPRγ gene gave rise to the mPRδ and ε genes much later in the vertebrate lineage. In addition, the mPRα gene first appeared in cartilaginous fishes, likely derived from duplication of mPRβ after the agnathan-gnathostome divergence. All known MAPR genes were present in the lamprey genomes. Progesterone receptor membrane component 1 (PGRMC1), neudesin and neuferricin genes probably evolved in parallel in non-bilaterians, whereas two copies of PGRMC genes probably derived from duplication of ancestral PGRMC1 sequence and appeared before the speciation of lampreys. Non-classical mPR and MAPR genes first evolved in non-bilaterians and classical nPR genes evolved later in basal vertebrates. Sequence repertoires for membrane progesterone receptor genes in vertebrates likely originated from an ancestral metazoan sequence and expanded via several duplication events.

Geological events and landscape features, as well as changes in the climate during the Pliocene period, have shaped the distribution of genetic diversity and demographic history of freshwater fish in mainland China. In this study, we investigated the phylogeny and population genetic structure of Acrossocheilus species (A. fasciatus, A. kreyenbergii and A. wenchowensis) in the Zhejiang region by the mitochondrial cytochrome b (cyt-b) and control (D-loop) region s. Mitochondrial phylogenetic analysis revealed three major lineages (lineages A, B and C), which may represent A. fasciatus, A. wenchowensis and A. kreyenbergii, respectively. Our results revealed that A. fasciatus and A. wenchowensis diverged from A. kreyenbergii in the Zhejiang region. The uplift of the Wuyi Mountains and the Xianxia Mountains served as an important geographic barrier in the diversification of the three Acrossocheilus species. The most recent common ancestors (TMRCAs) of the three lineages and lineages A + B were dated to 1.859 and 1.614 myr, respectively. Our results indicate that the effective population sizes of A. fasciatus and A. wenchowensis remained constant from the past to the present, as supported by ABC analysis. For conservation and protection, a strategy is required because of their genetic uniqueness, and we suggest that the two regions divided by the Xianxia Mountains be regarded as different management units (Mus), conforming to the major zoological regions of the Zhejiang region.

MicroRNAs (miRNAs) are regarded as key regulators in gonadal development and sex determination in diverse organisms. However, the functions of miRNAs in gonads of Acrossocheilus fasciatus, an economically important freshwater species in the south of China, are still unclear. Here, high-throughput sequencing was performed to investigate the mRNA and miRNAs on gonads of A. fasciatus. In total, 49,447 unigenes were obtained, including 11,635 differentially expressed genes (DEGs), among which 4147 upregulated genes and 7488 downregulated genes in the testis compared to the ovary, while 300 (237 known, and 63 novel) miRNAs with 36 differentially expressed miRNAs (DEMs) were identified, from which 17 upregulated and 19 downregulated DEMs. GO and KEGG enrichment analysis were performed to analyze the potential biological functions of DEGs and DEMs. Using qRT-PCR, 9 sex-related genes and 9 miRNAs were selected to verify the sequencing data. By dual-luciferase reporter assay, miR-22a-5p and miR-22b-5p interaction with piwil1, and miR-10d-5p interaction with piwil2 were identified. These findings could provide a reference for miRNA-regulated sex control of A. fasciatus and may reveal new insights into aquaculture and breeding concepts.