Assembly of Genome and Resequencing Provide Insights into Genetic Differentiation between Parents of Hulong Hybrid Grouper (Epinephelus fuscoguttatus ♀ × E. lanceolatus ♂).

The Hulong hybrid grouper was bred from the brown-marbled grouper (Epinephelus fuscoguttatus) ♀ and the giant grouper (E. lanceolatus) ♂, combining the advantageous traits of both parents. Possessing an excellent performance, this hybrid's cultivation promotes the development of the grouper industry. Its male parent, the giant grouper, possesses the fastest growth and the largest body size among all coral-reef-dwelling fish. This species is not only an economically important species in marine aquaculture, but it is also an ideal male parent in the interspecific crossing of grouper species. In the present study, a high-quality chromosome-level genome of the giant grouper was constructed with a total length of 1.06 Gb, consisting of 24 chromosomes and 69 scaffolds. To analyze the genetic differences between the parents of the Hulong hybrid grouper, the structural variations (SVs) between both parental genomes were detected, and a total of 46,643 SVs were obtained. High-quality SNPs were identified from resequencing data. There were significant differences between the two genomes, and the average FST reached 0.685. A total of 234 highly differentiated regions were detected with an FST > 0.9. The protein-coding genes involved in SVs and highly differentiated regions were significantly enriched in metabolic pathways, including fatty metabolism, carbohydrate metabolism, amino acid metabolism and the TCA cycle. These genes may be related to the differences in feeding preferences and the ability to digest carbohydrates between the two grouper species under natural conditions. In addition, protein-coding genes related to the cell cycle and p53-signaling pathway were also detected. These genes may play important roles in the regulation of body size and growth performance. This research provides genomic resources for further breeding works and evolutionary analyses.

First identification of two co-existing genome-wide significant sex quantitative trait loci (QTL) in red tilapia using integrative QTL mapping.

Red tilapia ( Oreochromis spp .) is one of the most popular fish in China due to its bright red appearance, fast growth rate, and strong adaptability. Understanding the sex determination mechanisms is of vital importance for the selection of all-male lines to increase aquacultural production of red tilapia. In this research, the genetic architecture for sex from four mapping populations ( n=1 090) of red tilapia was analyzed by quantitative trait loci (QTL)-seq, linkage-based QTL mapping, and linkage disequilibrium (LD)-based genome-wide association studies. Two genome-wide significant QTL intervals associated with sex were identified on ChrLG1 (22.4-23.9 Mb) and ChrLG23 (32.0-35.9 Mb), respectively. The QTL on ChrLG1 was detected in family 1 (FAM1), FAM2, and FAM4, and the other QTL on ChrLG23 was detected in FAM3 and FAM4. Four microsatellite markers located within the QTL were successfully developed for marker-assisted selection. Interestingly, three ( lpp, sox14, and amh) of the 12 candidate genes located near or on the two QTL intervals were abundantly expressed in males, while the remaining genes were more highly expressed in females. Seven genes ( scly, ube3a, lpp, gpr17, oca2, cog4, and atp10a) were significantly differentially expressed between the male and female groups. Furthermore, LD block analysis suggested that a cluster of genes on ChrLG23 may participate in regulating sex development in red tilapia. Our study provides important information on the genetic architecture of sex in red tilapia and should facilitate further exploration of sex determination mechanisms in this species.

Whole-genome sequencing of brown-marbled grouper (Epinephelus fuscoguttatus) provides insights into adaptive evolution and growth differences.

The brown-marbled grouper (Epinephelus fuscoguttatus) is an important species of fish in the coral reef ecosystem and marine aquaculture industry. In this study, a high-quality chromosome-level genome of brown-marbled grouper was assembled using Oxford Nanopore technology and Hi-C technology. The GC content and heterozygosity were approximately 42% and 0.35%, respectively. A total of 230 contigs with a total length of 1047 Mb and contig N50 of 13.8 Mb were assembled, and 228 contigs (99.13%) were anchored into 24 chromosomes. A total of 24,005 protein-coding genes were predicted, among which 23,862 (99.4%) predicted genes were annotated. Phylogenetic analysis showed that brown-marbled grouper and humpback grouper were clustered into one clade that separated approximately 11-23 million years ago. Collinearity analyses showed that there was no obvious duplication of large fragments between chromosomes in the brown-marbled grouper. Genomes of the humpback grouper and giant grouper showed a high collinearity with that of the brown-marbled grouper. A total of 305 expanded gene families were detected in the brown-marbled grouper genome, which is mainly involved in disease resistance. In addition, a genetic linkage map with 3061.88 cM was constructed. Based on the physical and genetic map, one growth-related quantitative trait loci was detected in 32,332,447 bp of chromosome 20, and meox1 and etv4 were considered candidate growth-related genes. This study provides pivotal genetic resources for further evolutionary analyses and artificial breeding of groupers.

Chromosome Genome Assembly of Cromileptes altivelis Reveals Loss of Genome Fragment in Cromileptes Compared with Epinephelus Species.

The humpback grouper (Cromileptes altivelis), an Epinephelidae species, is patchily distributed in the reef habitats of Western Pacific water. This grouper possesses a remarkably different body shape and notably low growth rate compared with closely related grouper species. For promoting further research of the grouper, in the present study, a high-quality chromosome-level genome of humpback grouper was assembled using PacBio sequencing and high-throughput chromatin conformation capture (Hi-C) technology. The assembled genome was 1.013 Gb in size with 283 contigs, of which, a total of 143 contigs with 1.011 Gb in size were correctly anchored into 24 chromosomes. Moreover, a total of 26,037 protein-coding genes were predicted, of them, 25,243 (96.95%) genes could be functionally annotated. The high-quality chromosome-level genome assembly will provide pivotal genomic information for future research of the speciation, evolution and molecular-assisted breeding in humpback groupers. In addition, phylogenetic analysis based on shared single-copy orthologues of the grouper species showed that the humpback grouper is included in the Epinephelus genus and clustered with the giant grouper in one clade with a divergence time of 9.86 Myr. In addition, based on the results of collinearity analysis, a gap in chromosome 6 of the humpback grouper was detected; the missed genes were mainly associated with immunity, substance metabolism and the MAPK signal pathway. The loss of the parts of genes involved in these biological processes might affect the disease resistance, stress tolerance and growth traits in humpback groupers. The present research will provide new insight into the evolution and origin of the humpback grouper.

[Biosynthesis and metabolism regulation of terpenoids in Pogostemon cablin: a review].

The terpenoids in Pogostemon cablin have complex structures and abundant pharmacological effects. Patchouli alcohol(PA) and pogostone(PO) have a high medicinal value by virtue of anti-tumor, anti-inflammatory, antibacterial, antioxidant, and other biological activities. Due to the low content of terpenoid metabolites in P. cablin, the study of biosynthesis and metabolism regulation can provide a biosynthetic basis for obtaining high-content terpenoids. In this study, key enzyme genes in biosynthesis, transcription factors in metabolism regulation, spatio-temporal expression of terpene synthase were reviewed, aiming to provide a reference for the development, protection, and utilization of P. cablin resources.

Whole-genome sequencing of leopard coral grouper ( Plectropomus leopardus) and exploration of regulation mechanism of skin color and adaptive evolution.

Leopard coral groupers belong to the Plectropomus genus of the Epinephelidae family and are important fish for coral reef ecosystems and the marine aquaculture industry. To promote future research of this species, a high-quality chromosome-level genome was assembled using PacBio sequencing and Hi-C technology. A 787.06 Mb genome was assembled, with 99.7% (784.57 Mb) of bases anchored to 24 chromosomes. The leopard coral grouper genome size was smaller than that of other groupers, which may be related to its ancient status among grouper species. A total of 22 317 protein-coding genes were predicted. This high-quality genome of the leopard coral grouper is the first genomic resource for Plectropomus and should provide a pivotal genetic foundation for further research. Phylogenetic analysis of the leopard coral grouper and 12 other fish species showed that this fish is closely related to the brown-marbled grouper. Expanded genes in the leopard coral grouper genome were mainly associated with immune response and movement ability, which may be related to the adaptive evolution of this species to its habitat. In addition, we also identified differentially expressed genes (DEGs) associated with carotenoid metabolism between red and brown-colored leopard coral groupers. These genes may play roles in skin color decision by regulating carotenoid content in these groupers.

Identification of Candidate Growth-Related SNPs and Genes Using GWAS in Brown-Marbled Grouper (Epinephelus fuscoguttatus).

Brown-marbled grouper, Epinephelus fuscoguttatus, is not only an important commercial fish species, but also an important crossbreeding parent in grouper industry. Improvement of growth traits of this species contributes to the development of grouper breeding. Currently, the development of molecular marker associated with growth of brown-marbled grouper is rare. Thus, we performed the first genome-wide association study (GWAS) for five growth traits in 172 brown-marbled groupers with 43,688 SNPs detected by ddRAD-seq. We identified a total of 5 significant and 18 suggestive QTLs located in multiple chromosomes associated with growth traits. In the 20 kb window of the significant SNPs and suggestive SNPs, 5 and 14 potential candidate genes affecting growth were detected, respectively. Five potential candidate genes near the significantly associated SNPs were selected for expression analysis. Among of which, bmp2k, wasf1, and acyp2 involved in bone development, maintenance of mitochondrion structure, and metabolism were differentially expressed. Interestingly, the SNP 23:29601315 located in the intron of bmp2k was significantly associated with body weight, body length, body height, and body thickness and suggestively associated with total length. We verified the locus using another new group including 123 individuals. The results showed that individuals with CC genotype have better growth traits comparing other individuals. Our findings not only contribute to understanding the molecular mechanism of growth regulation, but also promote the advance of marker-assisted selection in brown-marbled grouper.

First Genome-wide Association Analysis for Growth Traits in the Largest Coral Reef-Dwelling Bony Fishes, the Giant Grouper (Epinephelus lanceolatus).

The giant grouper, Epinephelus lanceolatus, is the largest coral reef-dwelling bony fish species. However, despite extremely fast growth performance and the considerable economic importance in this species, its genetic regulation of growth remains unknown. Here, we performed the first genome-wide association study (GWAS) for five growth traits in 289 giant groupers using 42,323 single nucleotide polymorphisms (SNPs) obtained by genotyping-by-sequencing (GBS). We identified a total of 36 growth-related SNPs, of which 11 SNPs reached a genome-wide significance level. The phenotypic variance explained by these SNPs varied from 7.09% for body height to 18.42% for body length. Moreover, 22 quantitative trait loci (QTLs) for growth traits, including nine significant QTLs and 13 suggestive QTLs, were found on multiple chromosomes. Interestingly, the QTL (LG17: 6934451) was shared between body weight and body height, while two significant QTLs (LG7: 22596399 and LG15: 11877836) for body length were consistent with the associated regions of total length at the genome-wide suggestive level. Eight potential candidate genes close to the associated SNPs were selected for expression analysis, of which four genes (phosphatidylinositol transfer protein cytoplasmic 1, protein tyrosine phosphatase receptor type E, alpha/beta hydrolase domain-containing protein 17C, and vascular endothelial growth factor A-A) were differentially expressed and involved in metabolism, development, response stress, etc. This study improves our understanding of the complex genetic architecture of growth in the giant grouper. The results contribute to the selective breeding of grouper species and the conservation of coral reef fishes.

Comparative transcriptome analyses reveal changes of gene expression in fresh and cryopreserved yellow catfish (Pelteobagrus fulvidraco) sperm and the effects of Cryoprotectant Me2SO.

This study, for the first time in fish, compared the transcriptome of fresh and frozen-thawed sperm, and would help to better understand the effect of cryopreservation on fish sperm and then better preserve the aquatic germplasm resources. Here, we employed high-throughput sequencing technology to obtain the transcriptome of yellow catfish from fresh sperm, cryopreserved sperm with and without cryoprotectant. When cryoprotectant (Me2SO) was excluded, down-regulated genes were significantly enriched into calcium ion binding, cytoskeletal protein binding, microfilament motor activity, calmodulin binding and carnitine O-acyltransferase activity, which affected Ca2+ regulation, cellular morphology, motility and metabolism. Moreover, heat shock proteins and genes associated with regulation of cholesterol, HCO3- and protein tyrosine phosphorylation (PTP) were down-regulated, and thus would impair ability against stress, membrane rigidity, pH regulation and signal transduction of cryopreserved sperm. After Me2SO was added, the amounts of DEGs decreased significantly and down-regulation of genes were found mainly in cytoskeleton and heat shock proteins, thereby suggesting that Me2SO effectively reduced the impact caused by low temperature on gene expression. Whether adding Me2SO or not, the up-regulated genes were mainly found in ribosomal proteins genes. However, when Me2SO was added, over-expression of some genes might contribute to maintain normal function of cryopreserved sperm.

Marker-assisted selection of YY supermales from a genetically improved farmed tilapia-derived strain.

Genetically improved farmed tilapia (GIFT) and GIFT-derived strains account for the majority of farmed tilapia worldwide. As male tilapias grow much faster than females, they are often considered more desirable in the aquacultural industry. Sex reversal of females to males using the male sex hormone 17-α-methyltestosterone (MT) is generally used to induce phenotypic males during large-scale production of all male fingerlings. However, the widespread use of large quantities of sex reversal hormone in hatcheries may pose a health risk to workers and ecological threats to surrounding environments. Breeding procedures to produce genetically all-male tilapia with limited or no use of sex hormones are therefore urgently needed. In this study, by applying marker-assisted selection (MAS) for the selection of YY supermales from a GIFT-derived strain, we identified 24 XY pseudofemale and 431 YY supermale tilapias. Further performance evaluation on the progenies of the YY supermales resulted in male rates of 94.1%, 99.5% and 99.6%, respectively, in three populations, and a daily increase in body weight of 1.4 g at 3 months (n=997). Our study established a highly effective MAS procedure in the selection of YY supermales from a GIFT-derived strain. Furthermore, the development of MAS-selected YY supermales will help reduce the utilization of hormones for controlling sex in the tilapia aquaculture.

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.

Complete mitochondrial DNA sequence of the yellowfin seabream Acanthopagrus latus and a genomic comparison among closely related sparid species.

The complete mitochondrial genome of the yellowfin seabream Acanthopagrus latus was determined in the present study. The genome was 16,609 bp in length and contained 37 genes (2 ribosomal RNA, 22 transfer RNA and 13 protein-coding genes) and the control region (CR), with the content and order of genes being similar to those in typical teleosts. Comparisons of the 37 genes and CR among species indicate the CR was the highest divergent (0.3341), but tRNA(Gly) possesses the lowest genetic variation (0.0542). Much greater p-genetic distances [mean = 0.1559, standard deviation (SD) = 0.0235; n = 1653] for the interspecies level with high frequency (99.4%) than those of the intraspecies level (mean = 0.0098, SD = 0.0090; n = 20) were inferred from 212 Cyt b sequence data, suggesting the Cyt b gene is conserved within Sparidae species and supporting the barcoding validity of Cyt b sequence data for Sparidae species identification. Phylogenetic analysis using amino acid sequences of 13 protein-coding genes supported that the genus Pagrus was not monophyletic, showing the need to re-evaluate the morphological characteristics of Pagrus fishes.

Population genetic structure of the Yangtze finless porpoise (Neophocaena phocaenoides asiaeorientalis): implications for management and conservation.

Understanding the population genetic structure is a prerequisite for conservation of a species. The degree of genetic variability characteristic of the mitochondrial DNA control region has been widely exploited in studies of population genetic structure and can be useful in identifying meaningful population subdivisions. To estimate the genetic profile of the Yangtze finless porpoise (Neophocaena phocaenoides asiaeorientalis), an endangered freshwater population endemic to China, the complete mtDNA control region was examined in 39 individuals belonging to seven different stocks inhabiting the middle and lower reaches of the Yangtze River. Very low genetic diversity was found (nucleotide diversity 0.0011 +/- 0.0002 and haplotypic diversity 0.65 +/- 0.05). The mtDNA genetic pattern of the Yangtze population appears to indicate a founder event in its evolutionary history and to support the marine origin for this population. Analyses by Fst and phi(st) yielded statistically significant population genetic structure (Fst = 0.44, P < 0.05; phi(st) = 0.36, P < 0.05). These results may have significant implications for the management and conservation of the Yangtze finless porpoise in the future.