Xinxin You

The Sinocyclocheilus cavefish genome provides insights into cave adaptation

BackgroundAn emerging cavefish model, the cyprinid genus Sinocyclocheilus, is endemic to the massive southwestern karst area adjacent to the Qinghai-Tibetan Plateau of China. In order to understand whether orogeny influenced the evolution of these species, and how genomes change under isolation, especially in subterranean habitats, we performed whole-genome sequencing and comparative analyses of three species in this genus, S. grahami, S. rhinocerous and S. anshuiensis. These species are surface-dwelling, semi-cave-dwelling and cave-restricted, respectively.ResultsThe assembled genome sizes of S. grahami, S. rhinocerous and S. anshuiensis are 1.75 Gb, 1.73 Gb and 1.68 Gb, respectively. Divergence time and population history analyses of these species reveal that their speciation and population dynamics are correlated with the different stages of uplifting of the Qinghai-Tibetan Plateau. We carried out comparative analyses of these genomes and found that many genetic changes, such as gene loss (e.g. opsin genes), pseudogenes (e.g. crystallin genes), mutations (e.g. melanogenesis-related genes), deletions (e.g. scale-related genes) and down-regulation (e.g. circadian rhythm pathway genes), are possibly associated with the regressive features (such as eye degeneration, albinism, rudimentary scales and lack of circadian rhythms), and that some gene expansion (e.g. taste-related transcription factor gene) may point to the constructive features (such as enhanced taste buds) which evolved in these cave fishes.ConclusionAs the first report on cavefish genomes among distinct species in Sinocyclocheilus, our work provides not only insights into genetic mechanisms of cave adaptation, but also represents a fundamental resource for a better understanding of cavefish biology.

Mudskipper genomes provide insights into the terrestrial adaptation of amphibious fishes

Mudskippers are amphibious fishes that have developed morphological and physiological adaptations to match their unique lifestyles. Here we perform whole-genome sequencing of four representative mudskippers to elucidate the molecular mechanisms underlying these adaptations. We discover an expansion of innate immune system genes in the mudskippers that may provide defence against terrestrial pathogens. Several genes of the ammonia excretion pathway in the gills have experienced positive selection, suggesting their important roles in mudskippers’ tolerance to environmental ammonia. Some vision-related genes are differentially lost or mutated, illustrating genomic changes associated with aerial vision. Transcriptomic analyses of mudskippers exposed to air highlight regulatory pathways that are up- or down-regulated in response to hypoxia. The present study provides a valuable resource for understanding the molecular mechanisms underlying water-to-land transition of vertebrates.

The Asian arowana (Scleropages formosus) genome provides new insights into the evolution of an early lineage of teleosts

The Asian arowana (Scleropages formosus), one of the world’s most expensive cultivated ornamental fishes, is an endangered species. It represents an ancient lineage of teleosts: the Osteoglossomorpha. Here, we provide a high-quality chromosome-level reference genome of a female golden-variety arowana using a combination of deep shotgun sequencing and high-resolution linkage mapping. In addition, we have also generated two draft genome assemblies for the red and green varieties. Phylogenomic analysis supports a sister group relationship between Osteoglossomorpha (bonytongues) and Elopomorpha (eels and relatives), with the two clades together forming a sister group of Clupeocephala which includes all the remaining teleosts. The arowana genome retains the full complement of eight Hox clusters unlike the African butterfly fish (Pantodon buchholzi), another bonytongue fish, which possess only five Hox clusters. Differential gene expression among three varieties provides insights into the genetic basis of colour variation. A potential heterogametic sex chromosome is identified in the female arowana karyotype, suggesting that the sex is determined by a ZW/ZZ sex chromosomal system. The high-quality reference genome of the golden arowana and the draft assemblies of the red and green varieties are valuable resources for understanding the biology, adaptation and behaviour of Asian arowanas.

Draft genome of the protandrous Chinese black porgy, Acanthopagrus schlegelii

Abstract Background As one of the most popular and valuable commercial marine fishes in China and East Asian countries, the Chinese black porgy (Acanthopagrus schlegelii), also known as the blackhead seabream, has some attractive characteristics such as fast growth rate, good meat quality, resistance to diseases, and excellent adaptability to various environments. Furthermore, the black porgy is a good model for investigating sex changes in fish due to its protandrous hermaphroditism. Here, we obtained a high-quality genome assembly of this interesting teleost species and performed a genomic survey on potential genes associated with the sex-change phenomenon. Findings We generated 175.4 gigabases (Gb) of clean sequence reads using a whole-genome shotgun sequencing strategy. The final genome assembly is approximately 688.1 megabases (Mb), accounting for 93% of the estimated genome size (739.6 Mb). The achieved scaffold N50 is 7.6 Mb, reaching a relatively high level among sequenced fish species. We identified 19 465 protein-coding genes, which had an average transcript length of 17.3 kb. By performing a comparative genomic analysis, we found 3 types of genes potentially associated with sex change, which are useful for studying the genetic basis of the protandrous hermaphroditism. Conclusions We provide a draft genome assembly of the Chinese black porgy and discuss the potential genetic mechanisms of sex change. These data are also an important resource for studying the biology and for facilitating breeding of this economically important fish.

High-throughput identification of novel conotoxins from the Chinese tubular cone snail (Conus betulinus) by multi-transcriptome sequencing

BackgroundThe venom of predatory marine cone snails mainly contains a diverse array of unique bioactive peptides commonly referred to as conopeptides or conotoxins. These peptides have proven to be valuable pharmacological probes and potential drugs because of their high specificity and affinity to important ion channels, receptors and transporters of the nervous system. Most previous studies have focused specifically on the conopeptides from piscivorous and molluscivorous cone snails, but little attention has been devoted to the dominant vermivorous species.ResultsThe vermivorous Chinese tubular cone snail, Conus betulinus, is the dominant Conus species inhabiting the South China Sea. The transcriptomes of venom ducts and venom bulbs from a variety of specimens of this species were sequenced using both next-generation sequencing and traditional Sanger sequencing technologies, resulting in the identification of a total of 215 distinct conopeptides. Among these, 183 were novel conopeptides, including nine new superfamilies. It appeared that most of the identified conopeptides were synthesized in the venom duct, while a handful of conopeptides were identified only in the venom bulb and at very low levels.ConclusionsWe identified 215 unique putative conopeptide transcripts from the combination of five transcriptomes and one EST sequencing dataset. Variation in conopeptides from different specimens of C. betulinus was observed, which suggested the presence of intraspecific variability in toxin production at the genetic level. These novel conopeptides provide a potentially fertile resource for the development of new pharmaceuticals, and a pathway for the discovery of new conotoxins.

Construction of the High-Density Genetic Linkage Map and Chromosome Map of Large Yellow Croaker (Larimichthys crocea)

High-density genetic maps are essential for genome assembly, comparative genomic analysis and fine mapping of complex traits. In this study, 31,191 single nucleotide polymorphisms (SNPs) evenly distributed across the large yellow croaker (Larimichthys crocea) genome were identified using restriction-site associated DNA sequencing (RAD-seq). Among them, 10,150 high-confidence SNPs were assigned to 24 consensus linkage groups (LGs). The total length of the genetic linkage map was 5451.3 cM with an average distance of 0.54 cM between loci. This represents the densest genetic map currently reported for large yellow croaker. Using 2889 SNPs to target specific scaffolds, we assigned 533 scaffolds, comprising 421.44 Mb (62.04%) of the large yellow croaker assembled sequence, to the 24 linkage groups. The mapped assembly scaffolds in large yellow croaker were used for genome synteny analyses against the stickleback (Gasterosteus aculeatus) and medaka (Oryzias latipes). Greater synteny was observed between large yellow croaker and stickleback. This supports the hypothesis that large yellow croaker is more closely related to stickleback than to medaka. Moreover, 1274 immunity-related genes and 195 hypoxia-related genes were mapped to the 24 chromosomes of large yellow croaker. The integration of the high-resolution genetic map and the assembled sequence provides a valuable resource for fine mapping and positional cloning of quantitative trait loci associated with economically important traits in large yellow croaker.

Molecular Evolution of Aralkylamine N-Acetyltransferase in Fish: A Genomic Survey

All living organisms synchronize biological functions with environmental changes; melatonin plays a vital role in regulating daily and seasonal variations. Due to rhythmic activity of the timezyme aralkylamine N-acetyltransferase (AANAT), the blood level of melatonin increases at night and decreases during daytime. Whereas other vertebrates have a single form of AANAT, bony fishes possess various isoforms of aanat genes, though the reasons are still unclear. Here, we have taken advantage of multiple unpublished teleost aanat sequences to explore and expand our understanding of the molecular evolution of aanat in fish. Our results confirm that two rounds of whole-genome duplication (WGD) led to the existence of three fish isoforms of aanat, i.e., aanat1a, aanat1b, and aanat2; in addition, gene loss led to the absence of some forms from certain special fish species. Furthermore, we suggest the different roles of two aanat1s in amphibious mudskippers, and speculate that the loss of aanat1a, may be related to terrestrial vision change. Several important sites of AANAT proteins and regulatory elements of aanat genes were analyzed for structural comparison and functional forecasting, respectively, which provides insights into the molecular evolution of the differences between AANAT1 and AANAT2.

Genome-Wide Mapping of Growth-Related Quantitative Trait Loci in Orange-Spotted Grouper (Epinephelus coioides) Using Double Digest Restriction-Site Associated DNA Sequencing (ddRADseq).

Mapping of quantitative trait loci (QTL) is essential for the discovery of genetic structures that related to complex quantitative traits. In this study, we identified 264,072 raw SNPs (single-nucleotide polymorphisms) by double digest restriction site associated DNA sequencing (ddRADseq), and utilized 3029 of these SNPs to construct a genetic linkage map in orange-spotted grouper (Epinephelus coioides) using a regression mapping algorithm. The genetic map contained 24 linkage groups (LGs) spanning a total genetic distance of 1231.98 cM. Twenty-seven significant growth-related QTLs were identified. Furthermore, we identified 17 genes (fez2, alg3, ece2, arvcf, sla27a4, sgk223, camk2, prrc2b, mchr1, sardh, pappa, syk, tert, wdrcp91, ftz-f1, mate1 and notch1) including three (tert, ftz-f1 and notch1) that have been reported to be involved in fish growth. To summarize, we mapped growth-related QTLs in the orange-spotted grouper. These QTLs will be useful in marker-assisted selection (MAS) efforts to improve growth-related traits in this economically important fish.

GHRH, PRP-PACAP and GHRHR Target Sequencing via an Ion Torrent Personal Genome Machine Reveals an Association with Growth in Orange-Spotted Grouper (Epinephelus coioides).

Growth hormone-releasing hormone (GHRH) and the receptor, GHRHR, constitute important components of the hypothalamus-pituitary growth axis and act on the downstream growth hormone (GH). PACAP-related peptide/pituitary adenylate cyclase activating polypeptide (PRP-PACAP) is a paralog of GHRH. These genes all play key roles in development and growth patterns. To improve the quality of cultured fish strains, natural genetic variation must be examined and understood. A mixed linear model has been widely used in association mapping, taking the population structures and pairwise kinship patterns into consideration. In this study, a mass cross population of orange-spotted grouper (Epinephelus coioides) was examined. These candidate genes were found to harbor low nucleotide diversity (θw from 0.00154 to 0.00388) and linkage disequilibrium levels (delay of 50% within 2 kbp). Association mapping was employed, and two single-nucleotide polymorphisms (KR269823.1:g.475A>C and KR269823.1:g.2143T>C) were found to be associated with growth (false discovery rate Q < 0.05), explaining 9.0%–17.0% of the phenotypic variance. The association of KR269823.1:g.2143T>C was also found via haplotype-based association (p < 0.05). The identified associations offer new insights into gene functions, and the associated single-nucleotide polymorphisms (SNPs) may be used for breeding purposes.

A chromosome-level genome assembly of the Asian arowana, Scleropages formosus

Asian arowana (Scleropages formosus), an ancient teleost belonging to the Order Osteoglossomorpha, has been a valuable ornamental fish with some varieties. However, its biological studies and breeding germplasm have been remarkably limited by the lack of a reference genome. To solve these problems, here we report high-quality genome sequences of three common varieties of Asian arowana (the golden, red and green arowana). We firstly generated a chromosome-level genome assembly of the golden arowana, on basis of the genetic linkage map constructed with the restriction site-associated DNA sequencing (RAD-seq). In addition, we obtained draft genome assemblies of the red and green varieties. Finally, we annotated 22,016, 21,256 and 21,524 protein-coding genes in the genome assemblies of golden, red and green varieties respectively. Our data were deposited in publicly accessible repositories to promote biological research and molecular breeding of Asian arowana.