Genome-Wide Association Mapping and Gene Expression Analyses Reveal Genetic Mechanisms of Disease Resistance Variations in Cynoglossus semilaevis

Mapping of Candidate Genes in Response to Low Nitrogen in Rice Seedlings

Nitrogen is not only a macronutrient essential for crop growth and development, but also one of the most critical nutrients in farmland ecosystem. Insufficient nitrogen supply will lead to crop yield reduction, while excessive application of nitrogen fertilizer will cause agricultural and eco-environment damage. Therefore, mining low-nitrogen tolerant rice genes and improving nitrogen use efficiency are of great significance to the sustainable development of agriculture. This study was conducted by Genome-wide association study on a basis of two root morphological traits (root length and root diameter) and 788,396 SNPs of a natural population of 295 rice varieties. The transcriptome of low-nitrogen tolerant variety (Longjing 31) and low-nitrogen sensitive variety (Songjing 10) were sequenced between low and high nitrogen treatments. A total of 35 QTLs containing 493 genes were mapped. 3085 differential expressed genes were identified. Among these 493 genes, 174 genes showed different haplotype patterns. There were significant phenotype differences among different haplotypes of 58 genes with haplotype differences. These 58 genes were hypothesized as candidate genes for low nitrogen tolerance related to root morphology. Finally, six genes (Os07g0471300, Os11g0230400, Os11g0229300, Os11g0229400, Os11g0618300 and Os11g0229333) which expressed differentially in Longjing 31 were defined as more valuable candidate genes for low-nitrogen tolerance. The results revealed the response characteristics of rice to low-nitrogen, and provided insights into regulatory mechanisms of rice to nitrogen deficiency.

Genomics and transcriptomics reveal new molecular mechanism of vibriosis resistance in fish

Infectious diseases have caused dramatic production decline and economic loss for fish aquaculture. However, the poor understanding of fish disease resistance severely hampered disease prevention. Chinese tongue sole (Cynoglossus semilaevis) is an important economic flatfish suffering from vibriosis. Here we used genomic, transcriptomic and experimental approaches to investigate the molecular genetic mechanisms underlying fish vibriosis resistance. A genome-wide comparison revealed that the genes under selective sweeps were enriched for glycosaminoglycan (GAG) chondroitin sulfate (CS)/dermatan sulfate (DS) metabolism. Transcriptomic analyses prioritized synergic gene expression patterns in this pathway, which may lead to an increased CS/DS content in the resistant family. Further experimental evidence showed that carbohydrate sulfotransferases 12 (Chst12), a key enzyme for CS/DS biosynthesis, has a direct antibacterial activity. To the best of our knowledge, this is the first report that the chst12 gene has a bactericidal effect. In addition, CS/DS is a major component of the extracellular matrix (ECM) and the selection signatures and fine-tuned gene expressions of ECM-receptor interaction genes indicated a modification in the ECM structure with an enhancement of the barrier function. Furthermore, functional studies conducted on Col6a2, encoding a collagen gene which constitutes the ECM, pointed to that it may act as a cellular receptor for Vibrio pathogens, thus plays an important role for the Vibrio invasion. Taken together, these findings provide new insights into the molecular protective mechanism underlying vibriosis resistance in fish, which offers crucial genomic resources for the resistant germplasm breeding and infectious disease control in fish culturing.

Intestinal microbiome-mediated resistance against vibriosis for Cynoglossus semilaevis

BACKGROUND

Infectious diseases have caused huge economic loss and food security issues in fish aquaculture. Current management and breeding strategies heavily rely on the knowledge of regulative mechanisms underlying disease resistance. Though the intestinal microbial community was linked with disease infection, there is little knowledge about the roles of intestinal microbes in fish disease resistance. Cynoglossus semilaevis is an economically important and widely cultivated flatfish species in China. However, it suffers from outbreaks of vibriosis, which results in huge mortalities and economic loss.

RESULTS

Here, we used C. semilaevis as a research model to investigate the host-microbiome interactions in regulating vibriosis resistance. The resistance to vibriosis was reflected in intestinal microbiome on both taxonomic and functional levels. Such differences also influenced the host gene expressions in the resistant family. Moreover, the intestinal microbiome might control the host immunological homeostasis and inflammation to enhance vibriosis resistance through the microbe-intestine-immunity axis. For example, Phaeobacter regulated its hdhA gene and host cyp27a1 gene up-expressed in bile acid biosynthesis pathways, but regulated its trxA gene and host akt gene down-expressed in proinflammatory cytokines biosynthesis pathways, to reduce inflammation and resist disease infection in the resistant family. Furthermore, the combination of intestinal microbes and host genes as biomarkers could accurately differentiate resistant family from susceptible family.

CONCLUSION

Our study uncovered the regulatory patterns of the microbe-intestine-immunity axis that may contribute to vibriosis resistance in C. semilaevis. These findings could facilitate the disease control and selective breeding of superior germplasm with high disease resistance in fish aquaculture. Video Abstract.

Identification of crucial factors involved in Cynoglossus semilaevis sexual size dimorphism by GWAS and demonstration of zbed1 regulatory network by DAP-seq

Sexual size dimorphism (SSD), whereby females and males exhibit different body sizes, are widely documented in animals. To explore crucial regulators implicated in female-biased SSD of Chinese tongue sole (Cynoglossus semilaevis), GWAS was conducted on 350 females and 59 males. Twenty SNPs and 25 genes including zbed1, nsd3, cdc45, klhl29, and smad4 with -log(p) > 7 were screened, mainly mapping to sex chromosome. The chromosome W-linked gene zbed1 attracted particular attention because it is a master key for cell proliferation. Thus, the regulatory network of zbed1 in C. semilaevis was explored by DAP-seq and 1352 peaks were discovered in the female brain. Moreover, zbed1 potentially regulated hippo signaling pathway, cell cycle, translation, and PI3k-Akt signaling pathway in C. semilaevis. These findings identify crucial SNPs and genes associated with female-biased SSD in C. semilaevis, also provide the first genome-wide survey for the zbed1 regulatory network in fish species.

Genome-wide association study reveals genetic variations associated with ocean acidification resilience in Yesso scallop Patinopecten yessoensis

Ocean acidification (OA), which refers to a gradual decrease in seawater pH due to the absorption of atmospheric carbon dioxide, profoundly affects the growth, development and survival of bivalves. Relatively limited studies have assessed the resilience of bivalve to OA. In the present study, Patinopecten yessoensis, an economically and ecologically significant species, were exposed to low pH (pH = 7.5) for 4 weeks. Forty-seven scallops that died in the first week were considered pH-sensitive population, and 20 that were alive at the end of the experiment were considered pH-tolerant population. A genome-wide association study was conducted to identify the genomic loci associated the resilience of P. yessoensis to OA. Twenty-one single nucleotide polymorphisms were significantly associated with resilience, which were distributed in 11 linkage groups. Within the linkage disequilibrium block region (± 300 kb) surrounding the 21 SNPs, 193 candidate genes were successfully identified. Particularly, five associated SNPs were directly located on five genes, including SP24, CFDH, 5HTR3, HSDL1 and ZFP346. The GO enrichment and KEGG pathway analyses showed that the molecular response of P. yessoensis to OA mainly involved neural signal transmission, energy metabolism and redox reaction. Candidate genes were expressed during larval development and in adult tissues. Furthermore, the expression of 30 candidate genes changed significantly under low pH stress in the mantle. Our results reveal certain SNPs and candidate genes that could elucidate the different responses of P. yessoensis to OA. The genetic variations indicated molecular resilience in P. yessoensis populations, which may enable adaptation to future acidification stress.

Identification and functional analysis of the perforin-1 like gene in disease resistance in half smooth tongue sole (Cynoglossus semilaevis)

The pore-forming protein perforin is one of the effectors of cell-mediated killing via the granule exocytosis pathway. In this study, a genome-wide association study was conducted in Vibrio harveyi disease-resistant and disease-susceptible families of half smooth tongue sole (Cynoglossus semilaevis) to determine the genes accounting for host resistance, and a perforin homologue was identified, designated perforin-1 like (CsPRF1l). The full-length cDNA of CsPRF1l is 1835 bp, and encodes 514 amino acids. The CsPRF1l gene consists of 10 exons and 9 introns, spanning approximately 7 kb. The amino acid sequence of CsPRF1l shows 60.35, 54.03, 41.92, and 34.17% identities to Morone saxatilis PRF1l, Oryzias melastigma PRF1l, Danio rerio PRF1.5 and Homo sapiens PRF, respectively. Sequence analysis revealed the presence of membrane attack complex/perforin (MACPF) and C2 domains in CsPRF1l. Quantitative real-time PCR showed that CsPRF1l presented a higher intestinal expression level in disease-resistant families than in susceptible families. Tissue expression pattern analysis showed that CsPRF1l is present in most of the tested tissues and highly expressed in the intestine, brain, stomach and gills. After challenge with V. harveyi, CsPRF1l mRNA was markedly upregulated in the liver, spleen, kidney, intestine, gills and skin. In addition, the recombinant CsPRF1l protein exhibited obvious antimicrobial activity against V. harveyi in vitro and in a zebrafish model. Collectively, these data indicate that CsPRF1l modulates host immune defense against V. harveyi invasion and provide clues about the efficacy of rCsPRF1l in fish that will give rise to useful therapeutic applications for V. harveyi infection in C. semilaevis.

Gene expression associated with disease resistance and long-term growth in a reef-building coral

Rampant coral disease, exacerbated by climate change and other anthropogenic stressors, threatens reefs worldwide, especially in the Caribbean. Physically isolated yet genetically connected reefs such as Flower Garden Banks National Marine Sanctuary (FGBNMS) may serve as potential refugia for degraded Caribbean reefs. However, little is known about the mechanisms and trade-offs of pathogen resistance in reef-building corals. Here, we measure pathogen resistance in Montastraea cavernosa from FGBNMS. We identified individual colonies that demonstrated resistance or susceptibility to Vibrio spp. in a controlled laboratory environment. Long-term growth patterns suggest no trade-off between disease resistance and calcification. Predictive (pre-exposure) gene expression highlights subtle differences between resistant and susceptible genets, encouraging future coral disease studies to investigate associations between resistance and replicative age and immune cell populations. Predictive gene expression associated with long-term growth underscores the role of transmembrane proteins involved in cell adhesion and cell-cell interactions, contributing to the growing body of knowledge surrounding genes that influence calcification in reef-building corals. Together these results demonstrate that coral genets from isolated sanctuaries such as FGBNMS can withstand pathogen challenges and potentially aid restoration efforts in degraded reefs. Furthermore, gene expression signatures associated with resistance and long-term growth help inform strategic assessment of coral health parameters.

Development of a 38 K single nucleotide polymorphism array and application in genomic selection for resistance against Vibrio harveyi in Chinese tongue sole, Cynoglossus semilaevis

Based on 1572 re-sequenced Chinese tongue sole (Cynoglossus semilaevis), we investigated the accuracy of four genomic methods at predicting genomic estimated breeding values (GEBVs) of Vibrio harveyi resistance in C. semilaevis when SNPs varying from 500 to 500 k. All methods outperformed the pedigree-based best linear unbiased prediction when SNPs reached 50 k or more. Then, we developed an SNP array "Solechip No.1" for C. semilaevis breeding using the Affymetrix Axiom technology. This array contains 38,295 SNPs with an average of 10.5 kb inter-spacing between two adjacent SNPs. We selected 44 candidates as the parents of 23 families and genotyped them by the array. The challenge survival rates of offspring families had a correlation of 0.706 with the mid-parental GEBVs. This SNP array is a convenient and reliable tool in genotyping, which could be used for improving V. harveyi resistance in C. semilaevis coupled with the genomic selection methods.

Transcriptome Profiling Reveals the Sexual Dimorphism of Gene Expression Patterns during Gonad Differentiation in the Half-Smooth Tongue Sole (Cynoglossus semilaevis)

The half-smooth tongue sole (Cynoglossus semilaevis), one of the most economically-important fish species in China, exhibits sexually dimorphic growth. An understanding of sex-related gene expression patterns in the tongue sole may inform sex regulation and breeding processes that increase fish production. However, the gene expression patterns during gonad development in the tongue sole remain unknown. In this study, transcriptome sequencing analyses were performed during gonad differentiation in the tongue sole, namely, at 62 days post-hatching (dph), 100 dph, 120 dph, and 150 dph. Differentially expressed genes associated with sex differentiation and gonad development were identified at each time point. Trend analysis showed that gene expression patterns varied over time. These expression patterns either explained common, non-sexually-dimorphic features or indicated significant sexual dimorphism. Transcript structure analyses identified both sex and time differences among samples. This study investigated the time-dependent expression patterns of several sex-related genes, including Dmrt1, Amh, Foxl2, aromatase encoding gene, Esr, and the Sox gene family, during gonad differentiation in the tongue sole. These results might clarify the significant sexual differences during early development in the tongue sole and might provide insight into the mechanisms controlling sex differentiation and development.

High-throughput metabolomics method based on liquid chromatography-mass spectrometry: Insights into the underlying mechanisms of salinity-alkalinity exposure-induced metabolites changes in Barbus capito

It is critical to investigate the adaptive development and the physiological mechanism of fish in external stimulation. In this study, the response of Barbus capito to salinity-alkalinity exposure was explored by high-throughput nontargeted and liquid chromatography-mass spectrometry-based metabolomics to investigate metabolic biomarker and pathway changes. Meanwhile, the biochemical indexes of Barbus capito were measured to discover the chronic impairment response to salinity-alkalinity exposures. A total of 29 tissue metabolites were determined to deciphering the endogenous metabolic changes of fishes during the different concentration salinity-alkalinity exposures environment, which were mainly involved in the key metabolism including the phenylalanine, tyrosine, and tryptophan biosynthesis, arachidonic acid metabolism, pyruvate metabolism, citrate cycle, and glycerophospholipid metabolism. Finally, we found the amino acid metabolism as key target was associated with the endogenous metabolites and metabolic pathways of Barbus capito to salinity-alkalinity exposures. In conclusion, metabolomics is a potentially powerful tool to reveal the mechanism information of fish in various exposure environments.

Combining Multiple Approaches and Models to Dissect the Genetic Architecture of Resistance to Infections in Fish

Infectious diseases represent a major threat for the sustainable development of fish farming. Efficient vaccines are not available against all diseases, and growing antibiotics resistance limits the use of antimicrobial drugs in aquaculture. It is therefore important to understand the basis of fish natural resistance to infections to help genetic selection and to develop new approaches against infectious diseases. However, the identification of the main mechanisms determining the resistance or susceptibility of a host to a pathogenic microbe is challenging, integrating the complexity of the variation of host genetics, the variability of pathogens, and their capacity of fast evolution and adaptation. Multiple approaches have been used for this purpose: (i) genetic approaches, QTL (quantitative trait loci) mapping or GWAS (genome-wide association study) analysis, to dissect the genetic architecture of disease resistance, and (ii) transcriptomics and functional assays to link the genetic constitution of a fish to the molecular mechanisms involved in its interactions with pathogens. To date, many studies in a wide range of fish species have investigated the genetic determinism of resistance to many diseases using QTL mapping or GWAS analyses. A few of these studies pointed mainly toward adaptive mechanisms of resistance/susceptibility to infections; others pointed toward innate or intrinsic mechanisms. However, in the majority of studies, underlying mechanisms remain unknown. By comparing gene expression profiles between resistant and susceptible genetic backgrounds, transcriptomics studies have contributed to build a framework of gene pathways determining fish responsiveness to a number of pathogens. Adding functional assays to expression and genetic approaches has led to a better understanding of resistance mechanisms in some cases. The development of knock-out approaches will complement these analyses and help to validate putative candidate genes critical for resistance to infections. In this review, we highlight fish isogenic lines as a unique biological material to unravel the complexity of host response to different pathogens. In the future, combining multiple approaches will lead to a better understanding of the dynamics of interaction between the pathogen and the host immune response, and contribute to the identification of potential targets of selection for improved resistance.