Whole-genome resequencing from bulked-segregant analysis reveals gene set based association analyses for the Vibrio anguillarum resistance of turbot (Scophthalmus maximus)

Cloning of Toll3 and Toll4 and association analysis among their SNP haplotypes and disease resistance in red swamp crayfish (Procambarus clarkii)

With the expansion of the culture scale of red swamp crayfish (Procambarus clarkii), the high incidence of diseases has seriously threatened the development of its industry. In this study, PcToll3 and PcToll4 were respectively cloned and explored SNPs among the germplasm populations, which had been identified relating to disease resistance in crayfish based on our previous study. A total of 3036 bp and 2820 bp of the open reading frame of PcToll3 and PcToll4 encoded 1011 and 939 amino acids, respectively. They were specially expressed in haemolymph, and significantly up-regulated expression after stimulation by Vibrio parahaemolyticus, Aeromonas hydrophila and white spot syndrome virus. It was found that the expression of downstream genes PcALF, PcCru, PcIMD, PcMyD88, and PcNF-κB were repressed after interference of PcToll3 and/or PcToll4. Totally, 16 and 19 SNPs in the coding region of PcToll3 and PcToll4 were mined, and the favoured haplotypes and the combinations of them were classified according to the associated SNPs with the disease resistance in crayfish. The haplotypes of Toll3-Hap1, Toll4-Hap1 and the combination of Toll3+Toll4-Hap1 were further validated that they had the stronger disease resistance comparing to others haplotypes, and the related KASP markers were developed for further breeding application. This study will advance our understanding of the function of the two Toll genes in crayfish, and provide the markers for the molecular breeding.

Characterization of Myf6 and association with growth traits in swimming crab (Portunus trituberculatus)

BACKGROUND

Myogenic factor 6 (Myf6) plays an important role in muscle growth and differentiation. In aquatic animals and livestock, Myf6 contributes to improving meat quality and strengthening the accumulation of muscle flavor substances. However, studies on Myf6 gene polymorphisms in crustaceans have not been reported.

RESULTS

In the current study, we characterized the Myf6 gene for Portunus trituberculatus to better understand its biological function. The full-length cDNA of Myf6 was 4,101 bp, with a 915 bp open reading frame encoding 304 amino acids. In addition, Myf6 included a conservative bHLH domain. Homology analysis showed that Myf6 shared the highest identity with Penaeus vannamei. Expression pattern analysis of Myf6 in fast- and slow-growing groups revealed that the expression level of the latter was significantly higher than that of the former (P < 0.05). qPCR studies revealed that Myf6 was expressed in various tissues with the highest level in muscle. Nineteen single nucleotide polymorphisms (SNPs) of Myf6 were identified and five of them were significantly associated with growth-related traits of P. trituberculatus (P < 0.05), including full carapace width, carapace length, body height, and body weight. The AG and GG genotypes of g.1,187,834 A > G exhibited superior growth-related traits than the AA genotype. In the combined genotypes of g.1,187,324 C > T and g.1,187,834 A > G, the average body weight of diplotype D5 (CT-GG) was higher than that of diplotype D1 (CC-AA), D2 (CC-AG), and D3 (CC-GG) in a cultivated population. A haploblock was generated by three significant SNPs (g.1187834 A > G, g.1188616 A > G, and g.1189024 C > A), containing four haplotypes (AAA, AAC, AGC, and GGC), among which GGC haplotype exhibited superior growth traits (full carapace width and body weight) than the AAA haplotype.

CONCLUSIONS

To our knowledge, this is the first report on Myf6 in crustaceans. The results of this study would contribute to elucidating multiple functions of the Myf6 gene in crustaceans and exploring the potential as a candidate gene in selective breeding programs of P. trituberculatus.

Mapping and Screening of Candidate Gene Regulating the Biomass Yield of Sorghum (Sorghum bicolor L.)

Biomass yield is one of the important traits of sorghum, which is greatly affected by leaf morphology. In this study, a lobed-leaf mutant (sblob) was screened and identified, and its F2 inbred segregating line was constructed. Subsequently, MutMap and whole-genome sequencing were employed to identify the candidate gene (sblob1), the locus of which is Sobic.003G010300. Pfam and homologous analysis indicated that sblob1 encodes a Cytochrome P450 protein and plays a crucial role in the plant serotonin/melatonin biosynthesis pathway. Structural and functional changes in the sblob1 protein were elucidated. Hormone measurements revealed that sblob1 regulates both leaf morphology and sorghum biomass through regulation of the melatonin metabolic pathway. These findings provide valuable insights for further research and the enhancement of breeding programs, emphasizing the potential to optimize biomass yield in sorghum cultivation.

Validation of a QTL associated with resistance to Vibrio anguillarum in rainbow trout (Oncorhynchus mykiss)

Vibriosis is a bacterial disease in fish caused by the Gram negative bacterium Vibrio anguillarum with severe impact on rainbow trout (Oncorhynchus mykiss) farming. Sustainable control methods should be developed and we here show that marker assisted selective breeding of fish naturally resistant to the disease is feasible. We have validated the use of a single nucleotide polymorphism (SNP) marker SNP AX-89,945,921 (QTL on chromosome 21). The QTL was previously found associated with resistance to vibriosis and described following a genome wide association analysis (GWAS) of trout exposed to the bacterium. For this validation spawners were genotyped by use of the 57 K Axiom®Trout Microarray (Affymetrix) and homozygous male fish carrying the allele with the SNP AX-89,945,921 were then selected and used to fertilize eggs from outbred female trout resulting in fish all carrying the SNP (QTL-fish). Control fish (non-QTL fish) were produced by fertilizing the same batch of eggs by use of male parents negative for the SNP. The fish were exposed in freshwater to V. anguillarum (water bath infection) at 19 C°. A total of 900 fish were challenged in a common garden set-up in triplicate. A bacterial solution of V. anguillarum (serotype O1) was added to each of three freshwater fish tanks, each with 150 QTL and 150 non-QTL fish. Fish were tagged by tail fin cut (upper/lower) to discern the two groups, whereafter fish were monitored around the clock to detect disease signs and remove moribund fish. Clinical vibriosis developed within two days in non-QTL-fish (overall morbidity of 70%). QTL fish developed clinical signs later and the morbidity was significantly lower and did not reach 50%. Rainbow trout farming may benefit from using the QTL associated with higher resistance towards vibriosis. The effect may be optimized in the future by use of both male and female parents homozygous for the marker allele.

Efficacy of bivalent vaccine against Aeromonas salmonicida and Edwardsiella tarda infections in turbot

In recent years, more than one pathogenic organism has usually been isolated from diseased turbot Scophthalmus maximus, creating a pressing need for the development of combination vaccines to prevent fish diseases brought on simultaneously by various infections. In this study, the inactivated bivalent vaccine of Aeromonas salmonicida and Edwardsiella tarda was prepared by the formalin inactivation method. After challenge with A. salmonicida and E. tarda at 4 weeks post-vaccination in turbot, the relative percentage survival (RPS) of the inactivated bivalent vaccine was 77.1%. In addition, we assessed the effects of the inactivated bivalent vaccine and evaluated the immunological processes after immunization in a turbot model. Serum antibody titer and lysozyme activity of the vaccinated group were both upregulated and higher than that in control group after vaccination. The expression levels of genes (TLR2, IL-1β, CD4, MHCI, MHCⅡ) that related to antigen recognition, processing and presentation were also studied in the liver, spleen and kidney tissues of vaccinated turbot. All the detected genes in the vaccinated group had a significant upward trend, and most of them reached the maximum value at 3-4 weeks, which had significant differences from the control group, suggesting that antigen recognition, processing and presentation pathway was activated by the inactivated bivalent vaccine. Our study provides a basis for further application of the killed bivalent vaccine against A. salmonicida and E. tarda in turbot, making it good potential that can be applied in aquaculture.

Effects of berberine hydrochloride on immune response in the crab Charybdis japonica

Berberine hydrochloride is the main effective component of Coptis spp. used in Chinese herbal medicine and its underlying molecular mechanisms, responsible for inducing effects in crustacean species, are not fully understood. In this study, the molecular response of the crab Charybdis japonica to berberine hydrochloride exposure was studied using transcriptome sequencing. The survival rate, gene expression and activities of several immune enzymes were measured after berberine hydrochloride treatments, with or without injection of the pathogenic bacterium Aeromonas hydrophila. A total of 962 differentially expressed genes (464 up-regulated and 498 down-regulated) were observed during exposure to 100 mg/L of berberine hydrochloride and in the control group after 48 h. Enrichment analysis revealed that these genes are involved in metabolism, cellular processes, signal transduction and immune functions, indicating that exposure to berberine hydrochloride activated the immune complement system. This bioactive compound simultaneously activated fibrinogen beta (FGB), fibrinogen alpha (FGA), alpha-2-macroglobulin (A2M), kininogen (KNG), fibrinogen gamma chain (FGB), alpha-2-HS-glycoprotein (AHSG), caspase-8 (CASP8), cathepsin L (CTSL), adenylate cyclase 3 (Adcy3) and MMP1. Its action could significantly increase the survival rate of the crabs injected with A. hydrophila and promote the activity of LZM, Caspas8, FGA, ACP and AKP in the hepatopancreas. When A. hydrophila was added, the neutralization of 300 mg/L berberine hydrochloride maximized the activities of Caspas8, LZM, ACP and AKP. Our results provide a new understanding of the potential effects of berberine hydrochloride on the immune system mechanisms in crustaceans.

High Genetic Connectivity Inferred from Whole-Genome Resequencing Provides Insight into the Phylogeographic Pattern of Larimichthys polyactis

Compared with terrestrial biota, marine fishes usually present lower genetic differentiation among different geographical populations because of high-level gene flow and lack of physical barriers. Understanding the genetic structure of marine fishes is essential for dividing management unit and making reasonable protection measures. The small yellow croaker (Larimichthys polyactis) belongs to the family Sciaenidae, which is an economic fish and widely distributed in the Western Pacific. To delineate genetic diversity and phylogeographic pattern, whole-genome resequencing was used to evaluate genetic connectivity, genetic diversity, and spatial pattern of L. polyactis for the first time. We obtained 6,645,711 high-quality single nucleotide polymorphisms (SNPs) markers from 40 L. polyactis individuals. The phylogenetic analysis, STRUCTURE, principal component analysis, and F_st results all indicated that no genetic structure consistent with the distribution pattern was found. This result revealed high genetic connectivity of L. polyactis in different sampling sites. High genetic diversity was also detected, indicating that there was sufficient evolutionary potential to maintain its effective population size. Besides, a similar result of high genetic connectivity and genetic diversity was also detected by mitochondria DNA marker. Our study demonstrated the persistence of high levels of genetic connectivity and a lack of population structure across L. polyactis in different sea areas. This study aimed to analyze the division of population structure and the reason for the decline and not exhaustion of L. polyactis resource on a genetic level.

Bulk segregation analysis in the NGS era: a review of its teenage years

Bulk segregation analysis (BSA) utilizes a strategy of pooling individuals with extreme phenotypes to conduct economical and rapidly linked marker screening or quantitative trait locus (QTL) mapping. With the development of next-generation sequencing (NGS) technology in the past 10 years, BSA methods and technical systems have been gradually developed and improved. At the same time, the ever-decreasing costs of sequencing accelerate NGS-based BSA application in different species, including eukaryotic yeast, grain crops, economic crops, horticultural crops, trees, aquatic animals, and insects. This paper provides a landscape of BSA methods and reviews the BSA development process in the past decade, including the sequencing method for BSA, different populations, different mapping algorithms, associated region threshold determination, and factors affecting BSA mapping. Finally, we summarize related strategies in QTL fine mapping combining BSA.

Selection of functional EPHB2 genotypes from ENU mutated grass carp treated with GCRV

Background

N-ethyl-N-nitrosourea (ENU) mutagenesis is a useful method for the genetic engineering of plants, and the production of functional mutants in animal models including mice and zebrafish. Grass carp reovirus (GCRV) is a haemorrhagic disease of grass carp which has caused noteworthy losses in fingerlings over the last few years. To overcome this problem, we used ENU mutant grass carp in an attempt to identify functional resistance genes for future hereditary rearing projects in grass carp.

Results

This study used ENU-mutated grass carp to identify genetic markers associated with resistance to the haemorrhagic disease caused by GCRV. Bulked segregant analysis (BSA) was performed on two homozygous gynogenetic ENU grass carp groups who were susceptible or resistant to GCRV. This analysis identified 466,162 SNPs and 197,644 InDels within the genomes of these mixed pools with a total of 170 genes annotated in the associated region, including 49 genes with non-synonymous mutations at SNP sites and 25 genes with frame shift mutations at InDel sites. Of these 170 mutated genes, 5 randomly selected immune-related genes were shown to be more strongly expressed in the resistant group as compared to the susceptible animals. In addition, we found that one immune-related gene, EPHB2, presented with two heterozygous SNP mutations which altered the animal’s responded to GCRV disease. These SNPs were found in the intron region of EPHB2 at positions 5859 (5859^G > A) and 5968 (5968^G > A) and were significantly (p = 0.002, 0.003) associated with resistance to GCRV. These SNP sites were also shown to correlate with the GCRV-resistant phenotype in these ENU grass carp. We also evaluated the mortality of the different ENU fish genotypes in response to GCRV and the SNPs in EPHB2. The outcomes of these evaluations will be useful in future selections of GCRV-resistant genes for genetic breeding in grass carp.

Conclusion

Our results provide a proof of concept for the application of BSA-sequence analysis in detecting genes responsible for specific functional genotypes and may help to develop better methods for marker-assisted selection, especially for disease resistance in response to GCRV.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07858-x.

A Major QTL for Resistance to Vibrio anguillarum in Rainbow Trout

Genetic selection of disease resistant fish is a major strategy to improve health, welfare and sustainability in aquaculture. Mapping of single nucleotide polymorphisms (SNP) in the fish genome may be a fruitful tool to define relevant quantitative trait loci (QTL) and we here show its use for characterization of Vibrio anguillarum resistant rainbow trout (Oncorhynchus mykiss). Fingerlings were exposed to the pathogen V. anguillarum serotype O1 in a solution of 1.5 × 10⁷ cfu/ml and observed for 14 days. Disease signs appeared 3 days post exposure (dpe) whereafter mortality progressed exponentially until 6 dpe reaching a total mortality of 55% within 11 days. DNA was sampled from all fish – including survivors – and analyzed on a 57 k Affymetrix SNP platform whereby it was shown that disease resistance was associated with a major QTL on chromosome 21 (Omy 21). Gene expression analyses showed that diseased fish activated genes associated with innate and adaptive immune responses. The possible genes associated with resistance are discussed.

Conjunctive Analyses of BSA-Seq and BSR-Seq to Reveal the Molecular Pathway of Leafy Head Formation in Chinese Cabbage

As the storage organ of Chinese cabbage, the leafy head was harvested as a commercial product due to its edible value. In this study, the bulked segregant analysis (BSA) and bulked segregant RNA-Seq (BSR) were performed with F2 separation progeny to study the molecular mechanism of leafy head formation in Chinese cabbage. BSA-Seq analysis located four candidate regions containing 40 candidate genes, while BSR-Seq analysis revealed eight candidate regions containing 607 candidate genes. The conjunctive analyses of these two methods identified that Casein kinase gene BrCKL8 (Bra035974) is the common candidate gene related with leafy head formation in Chinese cabbage, and it showed high expression levels at the three segments of heading type plant leaves. The differentially expressed genes (DEGs) between two set pairs of cDNA sequencing bulks were divided into two categories: one category was related with five hormone pathways (Auxin, Ethylene, Abscisic acid, Jasmonic acid and Gibberellin), the other category was composed of genes that associate with the calcium signaling pathway. Moreover, a series of upregulated transcriptional factors (TFs) were also identified by the association analysis of BSR-Seq analysis. The leafy head development was regulated by various biological processes and effected by diverse external environment factors, so our research will contribute to the breeding of perfect leaf-heading types of Chinese cabbage.

A functional polymorphism in the promoter of RhoB is associated with susceptibility to Vibrio anguillarum in turbot (Scophthalmus maximus)

As an isoform of Rho family GTPases, RhoB plays a pivotal role in cytoskeletal organization, cell proliferation, apoptosis and immune response. However, the regulatory mechanisms of RhoB expression in aquatic animals are still unknown. In the present study, we first construct Vibrio anguillarum infection model in S. maximus, including susceptible and resistant individuals. Then the temporal expression of RhoB was detected after V. anguillarum challenge using qRT-PCR and found that RhoB transcripts were significantly induced in the liver, gill and blood despite of differential expression levels and responsive time points. In addition, the mRNA levels of RhoB in resistant individuals were significantly higher than in susceptible ones. The length of 2083 bp sequences of RhoB promoter was cloned and characterized. Moreover, DNA methylation of the RhoB promoter was measured by bisulfite sequencing (BSP) and hypo-methylated was detected in the CpG islands. Three SNPs (-1590, -1575 and -1449) and two haplotypes in the promoter region of RhoB were identified to be associated with V. anguillarum resistance in turbot by association analysis in group 17-R and 17-S. Deletion analysis indicated that these SNPs could negatively mediate the activity of RhoB promoter. Site-directed mutagenesis and qRT-PCR of individuals with different genotypes demonstrated that -1575 T/A polymorphism affected promoter activity. Further study showed that this mutation altered the binding site of the transcription factor CREB. Co-transfection of SmCREB and RhoB promoter was performed in HEK293T cells which confirmed the -1575 allelic differences on transcriptional activity, with the susceptibility allele showing reduced activity. Taken together, our findings implicate that losing of binding of CREB to SmRhoB promoter due to -1575T/A polymorphisms enhances SmRhoB expression in resistant turbot, which provide insights into the effect of SmRhoB expression in response to V. anguillarum infection.

Functional differentiation of three phosphatidylinositol 3-kinase catalytic subunit alpha (PIK3CA) in response to Vibrio anguillarum infection in turbot (Scophthalmus maximus)

PIK3CA has been extensively investigated from its molecular mechanism perspective and association with its mutations in different types of cancers. However, little has been reported regarding the pathological significance of PIK3CA expression in teleost. Here, in our present study, three PIK3CA genes termed SmPIK3CAa, SmPIK3CAb and SmPIK3CA-like were firstly identified in the genome of turbot S. maximus. Although these three genes located in different chromosomes, all of them share the same five domains. Phylogenetic and synteny analysis indicated that SmPIK3CAa, SmPIK3CAb and SmPIK3CA-like were three paralogs that may originate from duplication of the same ancestral PIK3CA gene. Subcellular localization analysis confirmed the cytoplasm distribution of these three paralogs. All three SmPIK3CA were ubiquitously expressed in examined tissues in turbot, with the higher expression levels in immune-related tissues such as blood, spleen, kidney, gills and intestines. Upon Vibrio anguillarum challenge, SmPIK3CAa and SmPIK3CA-like transcripts were significantly induced in spleen, intestine and blood despite of differential expression levels and responsive time points. Additionally, individuals in resistant group showed significantly higher expression level of both two genes than in the susceptible group. Moreover, four SNPs (102, 2530, 3027 and 3060) and one haplotype (Hap2) located in exon region of SmPIK3CA-like were identified and confirmed to be associated with V. anguillarum resistance in turbot by association analysis in different populations. Taken together, these results suggested that functional differentiation occurred in three SmPIK3CA paralogs with Vibrio anguillarum resistance and SmPIK3CAa and SmPIK3CA-like probable play potential roles in innate immune response to pathogenic invasions in turbot.