Whole-genome resequencing of large yellow croaker (Larimichthys crocea) reveals the population structure and signatures of environmental adaptation

Population structure and adaptability analysis of Schizothorax o'connori based on whole-genome resequencing

BACKGROUND

Schizothorax o'connori is an endemic fish distributed in the upper and lower reaches of the Yarlung Zangbo River in China. It has experienced a fourth round of whole gene replication events and is a good model for exploring the genetic differentiation and environmental adaptability of fish in the Qinghai-Tibet Plateau. The uplift of the Qinghai-Tibet Plateau has led to changes in the river system, thereby affecting gene exchange and population differentiation between fish populations. With the release of fish whole genome data, whole genome resequencing has been widely used in genetic evolutionary analysis and screening of selected genes in fish, which can better elucidate the genetic basis and molecular environmental adaptation mechanisms of fish. Therefore, our purpose of this study was to understand the population structure and adaptive characteristics of S. o'connori using the whole-genome resequencing method.

RESULTS

The results showed that 23,602,746 SNPs were identified from seven populations, mostly distributed on chromosomes 2 and 23. There was no significant genetic differentiation between the populations, and the genetic diversity was relatively low. However, the Zangga population could be separated from the Bomi, Linzhi, and Milin populations in the cluster analysis. Based on historical dynamics analysis of the population, the size of the ancestral population of S. o'connori was affected by the late accelerated uplift of the Qinghai Tibet Plateau and the Fourth Glacial Age. The selected sites were mostly enriched in pathways related to DNA repair and energy metabolism.

CONCLUSION

Overall, the whole-genome resequencing analysis provides valuable insights into the population structure and adaptive characteristics of S. o'connori. There was no obvious genetic differentiation at the genome level between the S. o'connori populations upstream and downstream of the Yarlung Zangbo River. The current distribution pattern and genetic diversity are influenced by the late accelerated uplift of the Qinghai Tibet Plateau and the Fourth Ice Age. The selected sites of S. o'connori are enriched in the energy metabolism and DNA repair pathways to adapt to the low temperature and strong ultraviolet radiation environment at high altitude.

Transcriptomic analysis of large yellow croaker (Larimichthys crocea) reveals the suppression of the inflammatory response from Cryptocaryon irritans infection

Large yellow croaker (Larimichthys crocea) is the most productive marine fish in China. Cryptocaryon irritans is an extremely destructive parasite that causes great economic losses in large yellow croaker aquaculture industry. Therefore, it is very necessary to study the immune response of large yellow croaker in response to C. irritans infection. In this study, the transcriptomic profiles of large yellow croaker were sequenced and analyzed in the brain and head kidney at 72 h after C. irritans infection. Cytokines and chemokines related terms were significantly enriched based on the GO enrichment of down-regulated differentially expressed genes (DEGs) from the head kidney. Meanwhile, cytokine-cytokine receptor interaction was significantly enriched based on the KEGG enrichment of up-regulated DEGs from the brain and down-regulated DEGs from the head kidney, respectively. Moreover, the majority of inflammation-related DEGs were significantly up-regulated in the brain, but distinctly down-regulated in the head kidney. These results showed that the brain and head kidney might play different roles against C. irritans infection, and the inflammatory response of large yellow croaker may be restrained during C. irritans infection. Taken together, the transcriptomic analyses will be helpful to more comprehensively understand the immune mechanism of teleost against C. irritans infection, and provide a theoretical basis for the prevention and treatment of Cryptosporidiosis.

Different cold tolerances among three strains of large yellow croaker: related to antioxidant defense and energy metabolism

The aim of this study was to compare low-temperature tolerances in different strains of large yellow croaker. Dai Qu (DQ), Min-Yue Dong (MY), and Quan Zhou (NZ) strains of large yellow croaker were subjected to cold stress (8.6 °C) for 12 h, 24 h, 48 h, and 96 h. Survival rate, histological observation, and antioxidant and energy metabolism indicators were determined. The results showed that compared with the DQ group and MY group, NZ group aggravated hepatic structure, enhanced ROS, lactate, and anaerobic metabolism (PK gene expression and activity), while inhibited ATP, GSH, antioxidant enzymes (mRNA levels and activities of SOD, GPx, and CAT), and aerobic metabolism enzymes (mRNA levels and activities of F-ATPase, SDH, and MDH), indicating the reduction of cold tolerance in the NZ group was closely correlated with the decrement of antioxidative capacity and energy metabolism efficiency. Nrf2 and AMPK gene expressions were correlated with antioxidant and energy metabolism mRNA levels, respectively, suggesting Nrf2 and AMPK might participate in the modulation of target genes during the cold-stress adaptation. In conclusion, low temperature tolerance of fish depended on the antioxidant defense and energy metabolism efficiency, which contributes to understanding the underlying mechanisms of cold adaptation in large yellow croaker.

Whole-genome resequencing reveals recent signatures of selection in five populations of largemouth bass ( Micropterus salmoides)

Largemouth bass ( Micropterus salmoides) is an economically important fish species in North America, Europe, and China. Various genetic improvement programs and domestication processes have modified its genome sequence through selective pressure, leaving nucleotide signals that can be detected at the genomic level. In this study, we sequenced 149 largemouth bass fish, including protospecies (imported from the US) and improved breeds (four domestic breeding populations from China). We detected genomic regions harboring certain genes associated with improved traits, which may be useful molecular markers for practical domestication, breeding, and selection. Subsequent analyses of genetic diversity and population structure revealed that the improved breeds have undergone more rigorous genetic changes. Through selective signal analysis, we identified hundreds of putative selective sweep regions in each largemouth bass line. Interestingly, we predicted 103 putative candidate genes potentially subjected to selection, including several associated with growth (p sst1 and grb10), early development ( klf9, sp4, and sp8), and immune traits ( pkn2, sept2, bcl6, and ripk2). These candidate genes represent potential genomic landmarks that could be used to improve important traits of biological and commercial interest. In summary, this study provides a genome-wide map of genetic variations and selection footprints in largemouth bass, which may benefit genetic studies and accelerate genetic improvement of this economically important fish.

Genomic signatures of artificial selection in fecundity of Pacific white shrimp, Penaeus vannamei

Penaeusvannamei is the most important economic shrimp in the world. Many selective breeding programs are carried out to improve its production and performance traits. Although significant differences in the reproductive ability of female P. vannamei under artificial breeding conditions have been reported, the genome-wide adaption of the reproductive ability of domesticated female P. vannamei is less investigated. In this study, whole-genome analysis was performed along with pooled DNA sequencing on two fecundity separated bulks, high fecundity bulk (HB), and low fecundity bulk (LB). Each bulk contained 30 individuals from 3 commercial populations. A sequencing depth of >30× was achieved for each bulk, leading to the identification of 625,181 and 629,748 single nucleotide polymorphisms (SNPs) in HB and LB, respectively. Fixation index (Fst) combined with p ratio allowed for the identification of 145 selective sweep regions, with a sequence length of 14.5 Mb, accounting for 0.59% of the genome. Among the 145 selective sweep regions, a total of 64,046 SNPs were identified, and further verification was performed by genotyping 50 candidate SNPs on 60 samples from the offspring of the three populations. Furthermore, 121 genes were screened from the sweep regions. GO annotation and KEGG enrichment analyses showed that partial genes were essential for fecundity regulation. This study provides important information for in-depth investigation of genomic characteristics for long-term selective breeding on the fecundity of female P. vannamei and will also be important for genome-assisted breeding of P. vannamei in the future.

Genomic signatures of artificial selection in the Pacific oyster, Crassostrea gigas

The Pacific oyster, Crassostrea gigas, is an important aquaculture shellfish around the world with great economic and ecological value. Selective breeding programs have been carried out globally to improve production and performance traits, while genomic signatures of artificial selection remain largely unexplored. In China, we performed selective breeding of C. gigas for over a decade, leading to production of several fast‐growing strains. In the present study, we conducted whole‐genome resequencing of 20 oysters from two fast‐growing strains that have been successively selected for 10 generations, and 20 oysters from the two corresponding wild populations. Sequencing depth of >10× was achieved for each sample, leading to identification of over 12.20 million SNPs. The population structures investigated with three independent methods (principal component analysis, phylogenetic tree, and structure) suggested distinct patterns among selected and wild oyster populations. Assessment of the linkage disequilibrium (LD) decay clearly indicated the changes in genetic diversity during selection. Fixation index (F_st) combined with cross‐population composite likelihood ratio (XP‐CLR) allowed for identification of 768 and 664 selective sweeps (encompassing 1042 and 872 genes) tightly linked to selection in the two fast‐growing strains. KEGG enrichment and functional analyses revealed that 33 genes are important for growth regulation, which act as key components of various signaling pathways with close connection and further take part in regulating the process of cell cycle. This work provides valuable information for the understanding of genomic signatures for long‐term selective breeding and will also be important for growth study and genome‐assisted breeding of the Pacific oyster in the future.