Chromosomal‐level assembly of Takifugu obscurus (Abe, 1949) genome using third‐generation DNA sequencing and Hi‐C analysis

The Tetraodontidae family are known to have relatively small and compact genomes compared to other vertebrates. The obscure puffer fish Takifugu obscurus is an anadromous species that migrates to freshwater from the sea for spawning. Thus the euryhaline characteristics of T. obscurus have been investigated to gain understanding of their survival ability, osmoregulation, and other homeostatic mechanisms in both freshwater and seawater. In this study, a high quality chromosome‐level reference genome for T. obscurus was constructed using long‐read Pacific Biosciences (PacBio) Sequel sequencing and a Hi‐C‐based chromatin contact map platform. The final genome assembly of T. obscurus is 381 Mb, with a contig N50 length of 3,296 kb and longest length of 10.7 Mb, from a total of 62 Gb of raw reads generated using single‐molecule real‐time sequencing technology from a PacBio Sequel platform. The PacBio data were further clustered into chromosome‐scale scaffolds using a Hi‐C approach, resulting in a 373 Mb genome assembly with a contig N50 length of 15.2 Mb and and longest length of 28 Mb. When we directly compared the 22 longest scaffolds of T. obscurus to the 22 chromosomes of the tiger puffer Takifugu rubripes, a clear one‐to‐one orthologous relationship was observed between the two species, supporting the chromosome‐level assembly of T. obscurus. This genome assembly can serve as a valuable genetic resource for exploring fugu‐specific compact genome characteristics, and will provide essential genomic information for understanding molecular adaptations to salinity fluctuations and the evolution of osmoregulatory mechanisms.     

Plasma sex steroids and gonadosomatic index variations during ovarian development of female wild yellowfin seabream (Acanthopagrus latus)

The yellowfin seabream (Acanthopagrus latus) has a good potential for mariculture in Iran, so studies on various aspects of its culture has been carried out. In this study, the characteristic of reproduction in this species was investigated. The involvement of steroid hormones in the development of gonads was assayed by sampling blood and gonad tissue in 85 females of A. latus in fish caught from the Musa Creek, north-west Persian Gulf from October 2010 to May 2011. Gonadosomatic index (GSI) and hepatosomatic index (HSI) were recorded. Developmental stages were associated with changes in sex steroids, GSI and HSI. GSI and HSI began to increase in November (immature stage) and reached the highest value in March (ripe stage), decreasing sharply in April (spent stage). Plasma values of 17β-estradiol (E2) levels reached their highest levels in spawning females in March and then decreased in April. Progesterone levels were low through ovarian development, reached a maximum simultaneously with testosterone in February (mature stage) and then decreased in March (ripe stage). The seasonal changes observed in plasma steroid levels and GSI reflected the pattern of oocyte development and the spawning behavior of yellowfin seabream and were typical of the patterns described in most multiple spawners studied to date.     

Dipylidium caninum draft genome - a new resource for comparative genomic and genetic explorations of flatworms

Here, we present a draft genome of the tapeworm Dipylidium caninum (family Dipylidiidae) and compare it with other cestode genomes. This draft genome of D. caninum is 110 Mb in size, has a repeat content of ~13.4% and is predicted to encode ~10,000 protein-coding genes. We inferred excretory/secretory molecules (representing the secretome), other key groups of proteins (including peptidases, kinases, phosphatases, GTPases, receptors, transporters and ion-channels) and predicted potential intervention targets for future evaluation. Using 144 shared single-copy orthologous sequences, we investigated the genetic relationships of cestodes for which nuclear genomes are available. This study provides first insights into the molecular biology of D. caninum and a new resource for comparative genomic and genetic explorations of this and other flatworms.     

A genome-wide association study,supported by a new chromosome-level genome assembly,suggests sox2 as a main driver of the undifferentiatiated ZZ/ZW sex determination of turbot (Scophthalmus maximus)

BackgroundUnderstanding sex determination (SD) across taxa is a major challenge for evolutionary biology. The new genomic tools are paving the way to identify genomic features underlying SD in fish, a group frequently showing limited sex chromosome differentiation and high SD evolutionary turnover. Turbot (Scophthalmus maximus) is a commercially important flatfish with an undifferentiated ZW/ZZ SD system and remarkable sexual dimorphism. Here we describe a new long-read turbot genome assembly used to disentangle the genetic architecture of turbot SD by combining genomics and classical genetics approaches.ResultsThe new turbot genome assembly consists of 145 contigs (N50 = 22.9 Mb), 27 of them representing >95% of its estimated genome size. A genome wide association study (GWAS) identified a ~ 6.8 Mb region on chromosome 12 associated with sex in 69.4% of the 36 families analyzed. The highest associated markers flanked sox2, the only gene in the region showing differential expression between sexes before gonad differentiation. A single SNP showed consistent differences between Z and W chromosomes. The analysis of a broad sample of families suggested the presence of additional genetic and/or environmental factors on turbot SD.ConclusionsThe new chromosome-level turbot genome assembly, one of the most contiguous fish assemblies to date, facilitated the identification of sox2 as a consistent candidate gene putatively driving SD in this species. This chromosome SD system barely showed any signs of differentiation, and other factors beyond the main QTL seem to control SD in a certain proportion of families.     

Characterization of 11 polymorphic microsatellite loci in the yellowfin seabream Acanthopagrus latus

The development and characterization of polymorphic microsatellite loci in the yellowfin seabream Acanthopagrus latus are presented. Allelic diversity was estimated in 50 samples collected from southeastern China. Eleven loci displayed polymorphism that ranged from 10 to 25 alleles per locus and the observed heterozygosity ranged from 0.48 to 0.91. These markers could be very useful for the study of population genetics of the yellowfin seabream.     

Comprehensive whole genome survey analyses of male and female brown-spotted flathead fish Platycephalus sp.1

The flathead fish Platycephalus sp.1 is an ecologically and commercially important marine fish in the northwestern Pacific with notable sexual differences in growth and development. Yet the genomic data of this species is lacking. In the present study, whole genome sequencing of two individuals (one male and one female) of Platycephalus sp.1 were conducted to provide fundamental genomic information. The genome sizes were estimated to be 674.96 Mb (male) and 684.15 Mb (female) by using k-mer analyses. The heterozygosity and repeat ratios suggested possible male heterogamety of this species. The draft genome sequences were initially assembled and genome-wide microsatellite motifs were identified. Besides, the complete mitochondrial genome sequences were assembled and the phylogenetic analyses genetically supported the validation of Platycephalus sp.1. The reported genomic data and genetic markers in this study could be useful in future comparative genomics and evolutionary biology studies.     

Chromosome-level Genome Assembly of Acanthopagrus latus Provides Insights into Salinity Stress Adaptation of Sparidae

Marine Biotechnology - The yellowfin seabream, Acanthopagrus latus, is widely distributed throughout the Indo-West Pacific. This species, as a euryhaline Sparidae fish, inhabits in coastal...     

Dual in vitro effects of cortisol on cell turnover in the medaka esophagus via the glucocorticoid receptor

AimsCortisol is a glucocorticoid in mammals, but has both gluco- and mineralocorticoid activities in teleost fish. Our previous in vivo studies on osmoregulatory esophagi of euryhaline fish showed that epithelial apoptosis for the simple epithelium in seawater and cell proliferation for the stratified epithelium in fresh water are both induced by cortisol. The aim of the present study was to examine the mechanism of these dual cortisol effects on esophageal cell turnover.Main methodsWe developed a tissue culture method for the esophagus from euryhaline medaka (Oryzias latipes) and assessed cell proliferation and apoptosis in vitro in response to cortisol and 11-deoxycorticosterone (DOC), a recently identified agonist of the teleostean mineralocorticoid receptor.Key findingsEpithelial apoptosis, a well-established glucocorticoid function, was stimulated by treatment of the esophagus culture with 10 nM cortisol for 8 days, but no effects were seen at higher doses (100 and 1000 nM). In contrast, cell proliferation was induced by 1000 nM cortisol treatment for 8 days and this response was dose-dependent. Both effects were blocked by RU-486, a glucocorticoid receptor antagonist. DOC showed no significant effects at 10–1000 nM.SignificanceIn the esophageal epithelium in euryhaline fish, cortisol induces either apoptosis or cell proliferation via the glucocorticoid receptor, depending on the cortisol concentration. The glucocorticoid signaling may play a more important role than mineralocorticoid signaling in differentiation of the osmoregulatory esophagus in euryhaline fishes.     

The Chimonanthus salicifolius genome provides insight into magnoliid evolution and flavonoid biosynthesis

Chimonanthus salicifolius, a member of the Calycanthaceae of magnoliids, is one of the most famous medicinal plants in Eastern China. Here, we report a chromosome‐level genome assembly of C. salicifolius, comprising 820.1 Mb of genomic sequence with a contig N50 of 2.3 Mb and containing 36 651 annotated protein‐coding genes. Phylogenetic analyses revealed that magnoliids were sister to the eudicots. Two rounds of ancient whole‐genome duplication were inferred in the C. salicifolious genome. One is shared by Calycanthaceae after its divergence with Lauraceae, and the other is in the ancestry of Magnoliales and Laurales. Notably, long genes with > 20 kb in length were much more prevalent in the magnoliid genomes compared with other angiosperms, which could be caused by the length expansion of introns inserted by transposon elements. Homologous genes within the flavonoid pathway for C. salicifolius were identified, and correlation of the gene expression and the contents of flavonoid metabolites revealed potential critical genes involved in flavonoids biosynthesis. This study not only provides an additional whole‐genome sequence from the magnoliids, but also opens the door to functional genomic research and molecular breeding of C. salicifolius.     

A chromosome-level genome assembly of Ostrea denselamellosa provides initial insights into its evolution

The oyster Ostrea denselamellosa is a live-bearing species with a sharp decline in the natural population. Despite recent breakthroughs in long-read sequencing, high quality genomic data are very limited in O. denselamellosa. Here, we carried out the first whole genome sequencing at the chromosome-level in O. denselamellosa. Our studies yielded a 636 Mb assembly with scaffold N50 around 71.80 Mb. 608.3 Mb (95.6% of the assembly) were anchored to 10 chromosomes. A total of 26,412 protein-coding genes were predicted, of which 22,636 (85.7%) were functionally annotated. By comparative genomics, we found that long interspersed nuclear element (LINE) and short interspersed nuclear element (SINE) made up a larger proportion in O. denselamellosa genome than in other oysters'. Moreover, gene family analysis showed some initial insight into its evolution. This high-quality genome of O. denselamellosa provides a valuable genomic resource for studies of evolution, adaption and conservation in oysters.     

Stock discrimination and connectivity assessment of yellowfin seabream (Acanthopagrus latus) in northern South China Sea using otolith elemental fingerprints

Connectivity between fish stocks is fundamental to the understanding of population dynamics and the implementation of sustainable fisheries management. Otolith microchemistry is a promising tool as it can provide information on the continuous growth of otoliths and the environmental effects on otolith composition. Such elemental fingerprints can help distinguish different stocks or life history stages, identify the origins or nursery areas of fish, and assess population structure. In this study, we examined the stock discrimination and spatial connectivity of cage-cultured and wild stocks of yellowfin seabream (Acanthopagrus latus) from the coastal waters of Shantou, Yangjiang, and Zhanjiang in China southern province Guangdong during 2012–2014, based on otolith trace-elemental signatures using multivariate statistical analysis and machine learning approaches. The concentrations of 13 elements (⁷Li, ²³Na, ²⁴Mg, ⁴⁰Ca, ⁵⁵Mn, ⁵⁶Fe, ⁵⁹Co, ⁵⁹Ni, ⁶⁴Cu, ⁶⁵Zn, ⁸⁸Sr, ¹²²Sb, and ¹³⁷Ba) in the natal spot of fish otoliths, representing the embryonic and paralarval stages of fish, were analyzed using laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS). Stepwise discriminant analysis and random forests were used to distinguish the cultured and wild stocks of yellowfin seabream, and non-metric multidimensional scaling (NMDS) and cluster analysis were used to determine the spatial variation and connectivity of yellowfin seabream stocks. Overall, the cultured and wild stocks of yellowfin seabream could be identified with classification accuracy of 80.7% and 99.2% by using stepwise discriminant analysis and random forests respectively. When we compared site difference between cultured and wild stocks (site × stock interactions), the classification success was 60.4% for stepwise discriminant analysis and 85.7% for random forests. The misclassification of cultured and wild stocks within the three sites suggested the spatial connectivity between stocks and among sampling locations. Our findings suggested that the three wild stocks of yellowfin seabream from Guangdong coastal waters could be considered as one unit for management, and the difference between cultured and wild stocks was significant for yellowfin seabream from Shantou and Yangjiang, but less significant for yellowfin seabream from Zhanjiang. This study demonstrated that otolith elemental fingerprints can help improve our knowledge on the spatial connectivity, population structure, and life history of fish stocks, and random forests can be a useful tool for identifying cultured and wild stocks compared to the traditional stepwise discriminant analysis.     

The pomegranate (Punica granatum L.) draft genome dissects genetic divergence between soft‐ and hard‐seeded cultivars

Complete and highly accurate reference genomes and gene annotations are indispensable for basic biological research and trait improvement of woody tree species. In this study, we integrated single‐molecule sequencing and high‐throughput chromosome conformation capture techniques to produce a high‐quality and long‐range contiguity chromosome‐scale genome assembly of the soft‐seeded pomegranate cultivar ‘Tunisia’. The genome covers 320.31 Mb (scaffold N50 = 39.96 Mb; contig N50 = 4.49 Mb) and includes 33 594 protein‐coding genes. We also resequenced 26 pomegranate varieties that varied regarding seed hardness. Comparative genomic analyses revealed many genetic differences between soft‐ and hard‐seeded pomegranate varieties. A set of selective loci containing SUC8‐like, SUC6, FoxO and MAPK were identified by the selective sweep analysis between hard‐ and soft‐seeded populations. An exceptionally large selective region (26.2 Mb) was identified on chromosome 1. Our assembled pomegranate genome is more complete than other currently available genome assemblies. Our results indicate that genomic variations and selective genes may have contributed to the genetic divergence between soft‐ and hard‐seeded pomegranate varieties.     

The genome of the marine medaka Oryzias melastigma

Marine medaka (Oryzias melastigma) is considered to be a useful fish model for marine and estuarine ecotoxicology studies and has good potential for field‐based population genomics because of its geographical distribution in Asian estuarine and coastal areas. In this study, we present the first whole‐genome draft of O. melastigma. The genome assembly consists of 8,602 scaffolds (N50 = 23.737 Mb) and a total genome length of 779.4 Mb. A total of 23,528 genes were predicted, and 12,670 gene families shared with three teleost species (Japanese medaka, mangrove killifish and zebrafish) were identified. Genome analyses revealed that the O. melastigma genome is highly heterozygous and contains a large number of repeat sequences. This assembly represents a useful genomic resource for fish scientists.     

Finding Nemo’s Genes: A chromosome‐scale reference assembly of the genome of the orange clownfish Amphiprion percula

The iconic orange clownfish, Amphiprion percula, is a model organism for studying the ecology and evolution of reef fishes, including patterns of population connectivity, sex change, social organization, habitat selection and adaptation to climate change. Notably, the orange clownfish is the only reef fish for which a complete larval dispersal kernel has been established and was the first fish species for which it was demonstrated that antipredator responses of reef fishes could be impaired by ocean acidification. Despite its importance, molecular resources for this species remain scarce and until now it lacked a reference genome assembly. Here, we present a de novo chromosome‐scale assembly of the genome of the orange clownfish Amphiprion percula. We utilized single‐molecule real‐time sequencing technology from Pacific Biosciences to produce an initial polished assembly comprised of 1,414 contigs, with a contig N50 length of 1.86 Mb. Using Hi‐C‐based chromatin contact maps, 98% of the genome assembly were placed into 24 chromosomes, resulting in a final assembly of 908.8 Mb in length with contig and scaffold N50s of 3.12 and 38.4 Mb, respectively. This makes it one of the most contiguous and complete fish genome assemblies currently available. The genome was annotated with 26,597 protein‐coding genes and contains 96% of the core set of conserved actinopterygian orthologs. The availability of this reference genome assembly as a community resource will further strengthen the role of the orange clownfish as a model species for research on the ecology and evolution of reef fishes.     

Characterization of a Y‐specific duplication/insertion of the anti‐Mullerian hormone type II receptor gene based on a chromosome‐scale genome assembly of yellow perch,Perca flavescens

Yellow perch, Perca flavescens, is an ecologically and economically important species native to a large portion of the northern United States and southern Canada and is also a promising candidate species for aquaculture. However, no yellow perch reference genome has been available to facilitate improvements in both fisheries and aquaculture management practices. By combining Oxford Nanopore Technologies long‐reads, 10X Genomics Illumina short linked reads and a chromosome contact map produced with Hi‐C, we generated a high‐continuity chromosome‐scale yellow perch genome assembly of 877.4 Mb. It contains, in agreement with the known diploid chromosome yellow perch count, 24 chromosome‐size scaffolds covering 98.8% of the complete assembly (N50 = 37.4 Mb, L50 = 11). We also provide a first characterization of the yellow perch sex determination locus that contains a male‐specific duplicate of the anti‐Mullerian hormone type II receptor gene (amhr2by) inserted at the proximal end of the Y chromosome (chromosome 9). Using this sex‐specific information, we developed a simple PCR genotyping assay which accurately differentiates XY genetic males (amhr2by⁺) from XX genetic females (amhr2by^?). Our high‐quality genome assembly is an important genomic resource for future studies on yellow perch ecology, toxicology, fisheries and aquaculture research. In addition, characterization of the amhr2by gene as a candidate sex‐determining gene in yellow perch provides a new example of the recurrent implication of the transforming growth factor beta pathway in fish sex determination, and highlights gene duplication as an important genomic mechanism for the emergence of new master sex determination genes.     

A compact cartilaginous fish model genome

The genomes of several vertebrates, including six mammals, the chicken, Xenopus and four ray-finned fishes have been sequenced or are currently being sequenced to provide a better understanding of the human genome through comparative analysis. However, this list does not include cartilaginous fishes, which are the most basal living jawed vertebrates [1]. The genomes of the current ‘popular’ cartilaginous fishes such as the nurse shark, dogfish, and horn shark are larger than the human genome (∼3800 Mb to 7000 Mb) [2], and are not attractive for whole-genome sequencing. Here, we report the characterization of the relatively small genome (1200 Mb) of a cartilaginous fish, the elephant fish (Callorhinchus milii), and propose it as a model for whole-genome sequencing.     

A chromosome‐scale genome assembly of Antheraea pernyi (Saturniidae,Lepidoptera)

Antheraea pernyi is a semi‐domesticated lepidopteran insect species valuable to the silk industry, human health, and ecological tourism. Owing to its economic influence and developmental properties, it serves as an ideal model for investigating divergence of the Bombycoidea super family. However, studies on the karyotype evolution and functional genomics of A. pernyi are limited by scarce genomic resource. Here, we applied PacBio sequencing and chromosome structure capture technique to assemble the first high‐quality A. pernyi genome from a single male individual. The genome is 720.67 Mb long with 49 chromosomes and a 13.77‐Mb scaffold N50. Approximately 441.75 Mb, accounting for 60.74% of the genome, was identified as repeats. The genome comprises 21,431 protein‐coding genes, 85.22% of which were functionally annotated. Comparative genomics analysis suggested that A. pernyi diverged from its common ancestor with A. yamamai ~30.3 million years ago, and that chromosome fission contributed to the increased chromosome number. The genome assembled in this work will not only facilitate future research on A. pernyi and related species but also help to progress comparative genomics analyses in Lepidoptera.     

A Chromosome-level Genome Assembly of Wild Castor Provides New Insights into its Adaptive Evolution in Tropical Desert

Wild castor grows in the high-altitude tropical desert of the African Plateau, a region known for high ultraviolet radiation, strong light, and extremely dry condition. To investigate the potential genetic basis of adaptation to both highland and tropical deserts, we generated a chromosome-level genome sequence assembly of the wild castor accession WT05, with a genome size of 316 Mb, a scaffold N50 of 31.93 Mb, and a contig N50 of 8.96 Mb, respectively. Compared with cultivated castor and other Euphorbiaceae species, the wild castor exhibits positive selection and gene family expansion for genes involved in DNA repair, photosynthesis, and abiotic stress responses. Genetic variations associated with positive selection were identified in several key genes, such as LIG1, DDB2, and RECG1, involved in nucleotide excision repair. Moreover, a study of genomic diversity among wild and cultivated accessions revealed genomic regions containing selection signatures associated with the adaptation to extreme environments. The identification of the genes and alleles with selection signatures provides insights into the genetic mechanisms underlying the adaptation of wild castor to the high-altitude tropical desert and would facilitate direct improvement of modern castor varieties.     

Genome assembly of Scorias spongiosa and comparative genomics provide insights into ecological adaptation of honeydew-dependent sooty mould fungi

Sooty moulds are fungi of economic importance and with unique lifestyle mainly growing on insect honeydew. However, the lack of genomic data hinders investigation of genetic mechanisms underlying their ecological adaptation. With long-read sequencing technology, we generated the genome of Scorias spongiosa, an extraordinary sooty mould fungus associated with honeydew of colony aphids and producing large fruiting bodies. A 24.21 Mb high-quality genome assembly with a N50 length of 3.37 Mb was obtained. The genome contained 7758 protein coding genes, 97.13% of which were homologous to known genes, and approximately 0.29 Mb repeat sequences. Comparative genomics showed S. spongiosa lost relatively more gene families and contained fewer species-specific genes and gene families, with many CAZyme families and sugar transporters reduced or absent. This study not only promotes understanding of the ecological adaptation of sooty moulds, but also provides valuable genomic data resource for future comparative genomic and genetic studies.