Skip to main content

ORIGINAL RESEARCH article

Front. Genet., 27 July 2023
Sec. Livestock Genomics
This article is part of the Research Topic Advances in genomic research on diseases affecting aquaculture View all 9 articles

A simple sequence repeats marker of disease resistance in shrimp Litopenaeus vannamei and its application in selective breeding

Bin Yin,,Bin Yin1,2,3Haiyang Wang,,Haiyang Wang1,2,3Shaoping Weng,,Shaoping Weng1,2,3Sedong LiSedong Li4Jianguo He,,,
Jianguo He1,2,3,5*Chaozheng Li,,,
Chaozheng Li1,2,3,5*
  • 1State Key Laboratory of Biocontrol, Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai), School of Marine Sciences, School of Life Sciences, Sun Yat-sen University, Guangzhou, China
  • 2Guangdong Provincial Key Laboratory of Marine Resources and Coastal Engineering, Guangdong Provincial Key Laboratory for Aquatic Economic Animals, Guangzhou, China
  • 3China-ASEAN Belt and Road Joint Laboratory on Marine Aquaculture Technology, Guangzhou, China
  • 4Guangdong Evergreen Feed Industry Co., Ltd., Zhanjiang, China
  • 5Maoming Branch, Guangdong Laboratory for Lingnan Modern Agriculture, Maoming, China

The polymorphism of the simple sequence repeat (SSR) in the 5′ untranslated coding region (5′-UTR) of the antiviral gene IRF (LvIRF) has been shown to be implicated in the resistance to viral pathogens in shrimp Litopenaeus vannamei (L. vannamei). In this study, we explored the potential of this (CT)n-SSR marker in disease resistance breeding and the hereditary property of disease resistance traits in offspring. From 2018 to 2021, eight populations were generated through crossbreeding by selecting individuals according to microsatellite genotyping. Our results demonstrated that shrimp with the shorter (CT)n repeat exhibited higher resistance to white spot syndrome virus (WSSV) or Decapod iridescent virus 1 (DIV1); meanwhile, these resistance traits could be inherited in offspring. Interestingly, we observed that the longer (CT)n repeats were associated with bacterial resistance traits. Accordingly, shrimp with longer (CT)n repeats exhibited higher tolerance to Vibrio parahaemolyticus infection. Taken together, these results indicate that the single (CT)n-SSR marker could be used to selective breeding for both resistance to virus and bacteria in shrimps.

1 Introduction

Microsatellites, also known as simple sequence repeats (SSRs) or short tandem repeats (STRs), were first coined by Litt and Luty in 1989 (Litt and Luty, 1989), which are simple repeated motifs consisting of 1–6 base pairs found in both coding and non-coding regions. The advantage of SSRs as genetic markers is that they are inherited as co-dominant markers in a Mendelian fashion. Furthermore, high polymorphism rates and a broad distribution throughout the genome have made SSRs one of the popular genetic markers for use in breeding programs (Wright and Bentzen, 1994; Morgante et al., 2002; Gous et al., 2013). In fact, genetic markers have been widely applied in the breeding of economically important crops such as rice (Xue et al., 2018). By analysing the agronomic traits of rice varieties throughout the world, a number of SSRs have been obtained for marketable traits, such as grain size (GLW7/OsSPL13) (Si et al., 2016), chilling tolerance (COLD1) (Ma et al., 2015), disease resistance (Pigm) (Deng et al., 2017) and nitrate-use efficiency (NITR1.1) (Hu et al., 2015). SSR-assisted breeding has also been reported in aquatic animals. In shrimp Penaeus monodon, two specific SSR DNA fingerprints were present in considerably higher frequencies in the white spot syndrome virus (WSSV) susceptible shrimp population (Chakrabarty et al., 2015). In addition, eight IFN system genes were identified as anti-disease molecular markers for resistance breeding of gibel carp by transcriptome sequencing analysis (Mou et al., 2018).

Shrimp Litopenaeus vannamei is a worldwide aquaculture species, and the total of yield exceeds 4 million tons per year. However, various shrimp diseases are responsible for huge losses, especially WSSV-caused white spot syndrome (WSS) (Flegel, 2012). To date, there is no effective method for preventing or controlling WSSV infection in shrimp, but selective breeding of the WSSV-resistant population could be effective. In shrimp, innate immunity plays a key role in the defence against a wide variety of pathogenic invasions, such as bacteria, fungi and viruses. Notably, the IRF-Vago-JAK/STAT pathway, which is similar to the IRF-IFN-JAK/STAT pathway in vertebrates, has been functionally identified as playing a significant role in defence against multiple viruses, including WSSV (Li et al., 2015a). In this cascade, the L. vannamei interferon regulatory factor (LvIRF) can be activated during viral infection and then translocates to the nucleus to initiate the expression of the L. vannamei Vago4 (LvVago4) gene (Li et al., 2015a). LvVago4 is an arthropod cytokine encoding a viral-activated secreted peptide that inhibits viral infection by activating the JAK-STAT axis (Chen et al., 2011; Wang and He, 2019). Previously, we have shown that a microsatellite with a variable number of CT repeats presents in the 5′ untranslated coding region (5′-UTR) of LvIRF, and the length of the (CT)n can influence the expressional of LvIRF (Yin et al., 2019). As a result, the SSR in LvIRF (LvIRF-SSR) can be used as a molecular marker to selectively breed a new generation of WSSV-resistant shrimp.

In this study, we uncovered that shrimp with shorter (CT)n repeats in the 5′-UTR of LvIRF exhibited high tolerance to viruses, such as WSSV and decapod iridescent virus 1 (DIV1), while shrimp with longer (CT)n repeats in the 5′-UTR of LvIRF exhibited high tolerance to Vibrio parahaemolyticus (VP). We used LvIRF-SSR as a genetic marker for the breeding programme of shrimp L. vannamei. Through selecting individuals for cross-matching through microsatellite genotyping and artificial challenge for comparing resistance to virus and bacteria in offspring, we demonstrated that the LvIRF-SSR could be used in L. vannamei breeding to target improvement of antiviral and antibacterial traits.

2 Materials and methods

2.1 Microsatellite genotyping

As previously described in our study, the CT microsatellite from the LvIRF 5′-UTR region was genotyped using FAM fluorogenic probes to discriminate allele size (Yin et al., 2019). Primers (Table 1) were designed to amplify a 202 bp fragment containing 18 CT repeats that were used as a reference, and deviation from this size allowed us to deduce the number of (CT)n repeats for different alleles. Samples were sequenced and analyzed using an ABI Genetic Analyser 3730 XL. The alleles were divided into two groups: short (S) alleles were ≤ 18 repeats, and long (L) alleles were > 18 repeats, based on the median as previously described (Yin et al., 2019).

TABLE 1
www.frontiersin.org

TABLE 1. Primers used in this study.

2.2 Experimental animals

The L. vannamei breeding programme was carried out at the Hengxing shrimp farm in Zhanjiang City, China. The HX150301 and HX1403 populations served as the founder stocks, and two populations were stocked in F1. One population, known as the CTS, had parents with short CT repeats (S) in the LvIRF 5′-UTR, while the other population had parents with only long CT repeats (L), and was known as the CTL population. On the basis of F1, the same method was used for breeding F2 and F3 (Figure 1). To investigate the differences in the disease-resistant traits between F2 and F3 populations, a total of 400 shrimps were collected from each population, and 100 shrimps were received the treatment of WSSV, DIV1 and V. parahaemolyticus, as well as PBS as a control. The shrimp’s body weight was 6.0 ± 1.0 g each, and each population of shrimp was cultured in filtered sea water with 2.5% salinity at 26°C in a recirculating water tank and fed with fodder at a rate of 5% of body weight per day.

FIGURE 1
www.frontiersin.org

FIGURE 1. Schematic representation of the selective population. The HX150301 and HX1403 populations as founder stocks, with LvIRF-SSR genotyping used to select individuals as parents for cross-matching to obtain F1 in 2019. Based on F1, the same method was used for breeding F2 and F3 in 2020 and 2021, respectively. The homozygous populations CTSS and CTSS02, as well as the heterozygous populations CTS02 and CTS03, were separated from CTS with LvIRF-SSR genotyping.

2.3 Pathogens

The WSSV isolate (accession No. AF440570.1) originated from a batch of WSSV-infected L. vannamei collected in Zhanjiang in 2021, and the DIV1 isolate (accession No. NC_055165.1) originated from a batch of DIV1-infected L. vannamei collected in Maoming in 2021. WSSV and DIV1 were prepared from shrimp muscle tissue previously infected with WSSV or DIV1 and stored at −80°C. The muscle tissue was homogenised to prepare a WSSV inoculum at a final concentration of about 1 × 106 virions/50 μL and a DIV1 inoculum with a final concentration of about 1 × 105 virions/50 μL, which was injected into each shrimp, respectively (Li et al., 2014; Liao et al., 2020). All tanks were checked every 4 h, and dead shrimp were collected and marked. All shrimp were collected and marked on the fifth day (120 h) after infection to calculate cumulative mortality. All experimental materials were stored at −20°C with DNA and RNA holders (TakaRa, Dalian, China) for downstream DNA and RNA extractions.

The bacteria V. parahaemolyticus (VP) is sourced from Guangdong Microbiological Collection Center (GDMCC No. 1.306). VP was cultured in Luria broth (LB) medium overnight at 37°C, and the bacteria were harvested by centrifugation (5,000 g, 5 min) and washed twice in phosphate buffer saline (PBS) to remove growth medium. In a previous study, a final injection density of VP was adjusted to yield approximately 1 × 105 CFU/50 μL (Li et al., 2017). On the 5th day (120 h) after infection, all shrimp were collected and marked to estimate cumulative mortality.

2.4 Quantitative RT-PCR

Total RNA was extracted using the RNeasy Mini Kit (Qiagen, Hilden, Germany) and reverse transcribed into cDNA using the PrimeScript RT reagent kit (TakaRa, Dalian, China). Quantitative RT-PCR was performed to analyze LvIRF gene expression, and the LvEF-1α gene was detected as an internal control.

2.5 Statistical analysis

All data are presented as means ± SD. Student t-test was used to calculate the comparisons between groups of numerical data. Using GraphPad Prism software, cumulative mortality data were subjected to statistical analysis to generate the Kaplan ± Meier plot (log-rank χ2 test).

A test for goodness-of-fit analysis was performed to determine whether the alleles were inherited in a Mendelian fashion. The equation used with the respective degrees of freedom (df) was χ2=i=1kOE2/E, df = k-1, where O is the observed number of individuals in the genotypic class, E is the expected number of individuals in the same genotypic class and k is the number of genotypic classes (Eguiarte, 1998).

Associated analysis using 2 × 2 contingency tables and Fisher’s exact test. Bonferroni corrections were applied to account for the multiple testing of the LvIRF-SSR genotype groups [*p corr < 0.0167 (0.05/3)], following Fisher exact tests (Bush and Moore, 2012).

3 Results

3.1 Selective breeding of shrimp with LvIRF microsatellite marker

The L. vannamei breeding programme at Hengxing shrimp farm in Zhanjiang City, China, began in 2018, with a captive population in Hengxing serving as a founder stock. These populations were imported from abroad or domesticated in China for many years. The breeding programme was based on genotyping by the LvIRF-SSR marker to establish multiple families. As previously described in our study (Yin et al., 2019), the HX150301 and HX1403 populations were used as founder stocks, with LvIRF-SSR genotyping used to select individuals as parents for cross-matching. In the CTS population, where the parents contained allele S ((CT)n ≤ 18 repeats), the CTL population only had allele L ((CT)n > 18 repeats). Based on F1, the same method was used for breeding F2 and F3, and the homozygous population CTSS was separated from CTS in F2 (Figure 1).

3.2 LvIRF-SSR genotyping in offspring

Based on the breeding scheme illustrated in Figure 1, we obtained eight new shrimp populations in three generations from 2018 to 2021. LvIRF-SSR was used to genotype each population and count the allele frequencies to validate selection breeding and analysis of Mendelian inheritance. The statistical results of CTS and CTL were shown in previous studies: the CTS population contained the S allele, such as the (CT)14 and (CT)15 genotypes, while the CTL population did not, as expected (Yin et al., 2019). Similar to the genotyping results in F1, in F2 and F3, the homozygous populations CTSS and CTSS02 only contained the S allele of (CT)14 and (CT)15; CTL02 and CTL03 had the L allele of (CT)19, (CT)20, (CT)21, (CT)22 and (CT)23 (Figure 2). As expected, the heterozygous populations CTS02 and CTS03 contained S and L alleles.

FIGURE 2
www.frontiersin.org

FIGURE 2. Genotyped of LvIRF-SSR in F2 and F3. (A–C) The allelic distribution in F2 three populations: CTSS (A), CTS02 (B) and CTL02 (C). (D–F) The allelic distribution in F3 three populations: CTSS02 (D), CTS03 (E) and CTL03 (F).

The goodness-of-fit test in three heterozygous populations showed that the alleles of LvIRF-SSR did not segregate in a Mendelian ratio of 1:2:1. In the CTS population, the χ2 values of 21.90 was higher than 5.99 (df = 2), indicating that the heterozygous populations of CTS deviated from the Mendelian ratio of 1:2:1, manifested as homozygous of the LvIRF-SSR genotype SS deletion (Table 2). Homozygous LvIRF-SSR genotype SS deletion in the CTS02 and CTS03 populations also deviated from the Mendelian ratio of 1:2:1 (Table 2).

TABLE 2
www.frontiersin.org

TABLE 2. LvIRF-SSR genotypic frequencies among F1 and F2 populations and analysis of Mendelian inheritance.

3.3 Short CT microsatellite genotypes improve the resistance of shrimps to virus

WSSV and DIV1 were used to compare the resistance to virus in breeding populations. In F2, 120 h after WSSV infection, the cumulative mortality rate of the CTL02 population was 70%, significantly higher than those of CTSS and CTS02 populations (Figure 3A). Furthermore, the cumulative mortality of the CTL02 population was substantially higher than that of the CTS02 and CTL02 populations after 120 h of DIV1 infection (Figure 3C). The LvIRF gene expression in gills by quantitative RT-PCR revealed that after 24 h and 48 h of WSSV infection, the expression of LvIRF in CTSS and CTS02 populations was significantly higher than that in the CTL02 population, which was 1.85 times and 1.35 times after WSSV infection, respectively (Figure 3B). After DIV1 infection, the expression of LvIRF in the CTSS and CTS02 populations was 1.68-fold and 1.51-fold higher than that in the CTL02 population (Figure 3D). The results of the artificial challenge in the F2 generation showed that the resistance to virus of shrimps with the S allele (CTSS and CTS02 populations) were higher than those of the CTL02 population, and the increased expression of LvIRF could significantly improve the antiviral traits of the shrimp population. A similar result was shown in the F3 generation; the CTSS02 and CTS03 populations were resistant to the virus compared to CTL03 (Figures 3E, G), and the expression of LvIRF in CTSS02 and CTS03 populations was significantly higher than in CTL02 population after WSSV and DIV1 infection (Figures 3F, H). The above results suggest that the population of shrimp containing short CT genotypes in the 5′-UTR of LvIRF were resistant to WSSV and DIV1 infection, thereby improving their survival rate in this assay.

FIGURE 3
www.frontiersin.org

FIGURE 3. Differences in antiviral traits between shrimp populations in F2 and F3. (A,C) Cumulative mortality of F2 populations infected with WSSV and DIV1. (B,D) The expression of LvIRF in the gill after WSSV and DIV1 infection. (E,G) Cumulative mortality of F3 populations infected with WSSV and DIV1. (F,H) The expression of LvIRF in the gill after WSSV and DIV1 infection. Cumulative mortality data were analyzed statistically by the Kaplan ± Meier plot (log-rank χ2 test). *p < 0.05, **p < 0.01, ***p < 0.001, ns, not significant. Quantitative RT-PCR was performed for each sample using the LvEF-1α gene as an internal control by the Livak (2−ΔΔCT) method. Expression levels in the gill were used to determine the mean fold change (means ± SD, n = 9. *p < 0.05, **p < 0.01, ***p < 0.001, ns, not significant.

3.4 Long CT microsatellite genotypes are associated with VP resistance

VP was used to compare the resistance to bacteria in breeding populations. In F2, after 120 h of VP infection, the cumulative mortality rate of the CTSS population was close to 100%, significantly higher than the CTS02 and CTL02 populations (Figure 4A), while the CTL02 and CTS02 populations showed no significant difference in the anti-VP trait (p = 0.078). A similar result was shown in F3, and the anti-VP trait of the CTL03 population were substantially higher than those of the CTSS02 population (Figure 4A). The CTSS02 and CTS03 populations showed no significant difference in the anti-VP trait (p = 0.066), while the CTS03 and CTL03 populations showed no significant difference in the anti-VP trait (p = 0.177). After 120 h of VP infection, we divided each shrimp population into two groups based on survival status: dead (D) and surviving (S). Through microsatellite genotyping, we could attribute one of three bi-allelic genotypes to each individual: SS, SL or LL. Fisher’s exact test and Bonferroni corrections showed that SS and LL genotypes were associated with VP resistance. The CTL02 population showed a significant difference in the LL genotype (p corr < 0.0003), while the CTSS population showed a significant difference in the SS genotype (p corr < 0.0003). Similar to the CTL02 population, the LL genotype associated with VP resistance showed a substantial difference in the CTL03 population (p corr < 0.0003), while the SS genotype exhibited a significant difference in the CTSS02 population (p corr < 0.0003) (Table 3).

FIGURE 4
www.frontiersin.org

FIGURE 4. Differences in Vibrio parahaemolyticus resistance between shrimp populations in F2 and F3. (A) Analysis of Vibrio parahaemolyticus resistance differences among the three populations in F2. (B) Analysis of VP resistance differences among the three populations in F2. Cumulative mortality data were analyzed statistically by the Kaplan ± Meier plot (log-rank χ2 test). **p < 0.01, ***p < 0.001, ns, not significant.

TABLE 3
www.frontiersin.org

TABLE 3. Specific (CT)n genotypes associated with VP resistance.

We compared the bi-allelic genotype distribution between the survive and dead groups and shrimp with the LL genotype in the 5′-UTR of the LvIRF gene were resistant to VP, while those with SS were susceptible to VP. The above results suggested that shrimp with a long (CT)n repeat in the 5′-UTR of LvIRF were more resistant to VP than those with shorter (CT)n repeats.

3.5 Short CT microsatellite genotypes inhibited the growth traits of shrimps

On comparing body length and body weight among the three populations of F2, we found that the growth rate of the CTSS population was considerably lower than those of the CTS02 and CTL02 populations (Figure 5). In the CTSS population, 135 shrimps had an average body length of 7.70 cm and an average weight of 3.55 g. In the CTS02 and CTL02 populations, 136 and 142 shrimps were measured and counted, with an average body length of 8.22 cm and 8.37 cm and an average body weight of 4.24 g and 4.34 g, respectively. T-text analysis results show that the growth rate of the CTSS population was considerably lower than that of the CTS02 and CTL02 populations, but the growth rates of the CTS02 and CTL02 populations were not significantly different (Figure 5). The above results suggest that the shrimp with the homozygous LvIRF-SSR SS genotype may inhibit shrimp growth.

FIGURE 5
www.frontiersin.org

FIGURE 5. Differences in growth traits between shrimp populations in F2. (A) Analysis of body length differences among the three populations in F2. (B) Analysis of body weight differences among the three populations in F2. ***p < 0.001, ns, not significant.

4 Discussion

Our study examined the role of an SSR with a (CT)n motif in the 5′-UTR, its repeat length variation and the association of its genetic variations with virus and bacteria resistance. Several studies have suggested that SSR repeats exist in abundance in eukaryotic genes; untranslated coding regions (UTRs) have more SSRs than coding regions and exhibit a strong bias towards di- and tri-nucleotides repeats (Li et al., 2002; Fujimori et al., 2003; Zhao et al., 2014). The SSRs are present in the UTRs, indicating that they are likely involved in regulating target gene transcription (Simpson and Ayyar, 2008; Colak et al., 2014). Previous research has shown that the length of the SSR in LvIRF can influence the expressional level, and polymorphisms of LvIRF-SSR are associated with WSSV resistance traits in shrimp (Yin et al., 2019). In this study, our goal was to investigate the application of LvIRF-SSR in improving the antiviral traits in shrimp. Genotyping selects individuals for cross-matching, compared the differences in resistance to virus of different populations and found that shrimp with smaller numbers of (CT)n repeats exhibited enhanced tolerance to virus infection.

In aquatic breeding, molecular markers have been used to associate with and improve biological traits. Applications of marker-assisted selection for disease resistance included the breeding of a Japanese flounder (Paralichthys olivaceus) strain with resistance to the viral disease lymphocystis (Fuji et al., 2006; Fuji et al., 2007), and based on molecular markers affecting resistance to bacterial cold-water disease in rainbow trout (Oncorhynchus mykiss) (Liu et al., 2018). Another well-known example is the major molecular markers affecting resistance to the infectious pancreatic necrosis (IPN) virus in Atlantic salmon by QTL; marker-assisted selection significantly prevented IPN outbreaks (Houston et al., 2008; Moen et al., 2009). In addition, several molecular markers associated with resistance to pancreas disease (Gonen et al., 2015) and cardiomyopathy syndrome have been used in salmon breeding (Boison et al., 2019; Hillestad and Moghadam, 2019). In this study, we used the LvIRF-SSR as a molecular marker to establish shrimp breeding programmes with diverse genotypes and tolerances to viral infection. The population of shrimps containing short CT genotypes had a stronger antiviral immune response, as manifested by the observation of increased LvIRF gene expression and prolonged survival time after WSSV and DIV1 infection. In addition, the resistance to virus could be stably passed on to offspring by marker-assisted selection. Thus, LvIRF-SSR could be a candidate marker for shrimp breeding for virus resistance.

Furthermore, we noted that the CTSS population is resistant to viruses while susceptible to bacteria (VP). Previous studies have shown that activating the antibacterial NF-κB pathway in L. vannamei benefits the transcription and replication of WSSV genes (Qiu et al., 2017; Li et al., 2019a). Host transcription factors—Dorsal, Relish, c-Jun and c-Fos—have all been shown to regulate a bulk of viral genes, such as wsv069, wsv079 and wsv249 (Huang et al., 2010; Li et al., 2015b; Li et al., 2016). Meanwhile, the antiviral STING-IRF-Vago-JAK/STAT pathway can negatively regulate the activity of the NF-κB pathway (Li et al., 2019b). The transcriptional activity of microRNA-1 (miR-1) molecule from L. vannamei could be regulated by the JAK-STAT pathway, with miR-1 serving as the mediator, and STAT could exert a negative regulatory effect on the Dorsal pathway (Zuo et al., 2020). In other words, resistance to virus and bacteria in L. vannamei may have antagonistic effects; shrimps with strong resistance to virus are susceptible to bacteria, while those with strong resistance to bacteria are susceptible to virus. In our study, we investigated the bi-allelic genotype distribution and phenotypic variations in two generations of offspring. Through associated analysis and artificial challenge, we found that shrimps with the LL genotype in LvIRF-SSR were resistant to VP and susceptible to WSSV, while those with the SS genotype in LvIRF-SSR were resistant to WSSV and susceptible to VP. The LvIRF-SSR marker in selective breeding can be used to improve antiviral or antibacterial traits.

SSRs are co-dominant Mendelian markers, and they can inherit in accordance with Mendelian inheritance (Powell et al., 1996). In sugarcane (Lu et al., 2015), Pseudotsuga menziesii (Slavov et al., 2004) and peanut (Zhou et al., 2016), the segregation ratios of the SSR markers were consistent with the Mendelian inheritance law for single loci. However, our investigation discovered that the LvIRF-SSR of L. vannamei did not segregate in a Mendelian ratio of 1:2:1, but rather as a phenomena of homozygous deletion in offspring. We speculated that homozygous deletion might relate to bacterial susceptibility and slow growth of the genotype SS in shrimps. The CTSS population was susceptible to bacteria and grew more slowly compared to the CTS and CTL populations, indicating that shrimps with the SS allele are more likely to be eliminated in natural environments. This also explains the inability to segregate the SS allele of shrimp in the natural population when breeding the F1 generation.

In summary, we identified that the LvIRF-SSR marker could be used as a molecular marker in shrimp breeding. By cross-matching, enriching the dominant allele in offspring can improve the antiviral or antibacterial traits. Our results provide some insights into how this SSR could be used as a molecular marker in shrimp breeding.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding authors.

Author contributions

CL and JH conceived and designed the experiments. BY, HW, SW, and SL performed the experiments and analyzed data. CL and BY wrote the draft manuscript. CL was responsible for forming the hypothesis, project development, data coordination, and writing, finalizing, and submitting the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This research was supported by National Key Research and Development Program of China (2022YFD2400204), National Natural Science Foundation of China (32102777/32022085/31930113), the open competition program of top ten critical priorities of Agricultural Science and Technology Innovation for the 14th Five-Year Plan of Guangdong Province (2022SDZG01), Independent Research and Development Projects of Maoming Laboratory (2021TDQD004), Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai) (SML2021SP301), and the Fundamental Research Funds for the Central Universities, Sun Yat-sen University (22lglj05).

Conflict of interest

Author SL was employed by Guangdong Evergreen Feed Industry Co. Ltd.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

Boison, S., Ding, J., Leder, E., Gjerde, B., Bergtun, P. H., Norris, A., et al. (2019). QTLs associated with resistance to cardiomyopathy syndrome in atlantic salmon. J. Hered. 110, 727–737. PMID 31287894. doi:10.1093/jhered/esz042

PubMed Abstract | CrossRef Full Text | Google Scholar

Bush, W. S., and Moore, J. H. (2012). Chapter 11: Genome-wide association studies. PLoS Comput. Biol. 8 (12), e1002822. PMID 23300413. doi:10.1371/journal.pcbi.1002822

PubMed Abstract | CrossRef Full Text | Google Scholar

Chakrabarty, U., Dutta, S., Mallik, A., Mondal, D., and Mandal, N. (2015). Identification and characterisation of microsatellite DNA markers in order to recognise the WSSV susceptible populations of marine giant black tiger shrimp, Penaeus monodon. Vet. Res. 46, 110. PMID 26407974. doi:10.1186/s13567-015-0248-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, Y. H., Jia, X. T., Zhao, L., Li, C. Z., Zhang, S., Chen, Y. G., et al. (2011). Identification and functional characterization of Dicer2 and five single VWC domain proteins of Litopenaeus vannamei. Dev. Comp. Immunol. 35, 661–671. PMID 21256150. doi:10.1016/j.dci.2011.01.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Colak, D., Zaninovic, N., Cohen, M. S., Rosenwaks, Z., Yang, W. Y., Gerhardt, J., et al. (2014). Promoter-bound trinucleotide repeat mRNA drives epigenetic silencing in fragile X syndrome. Science 343, 1002–1005. PMID 24578575. doi:10.1126/science.1245831

PubMed Abstract | CrossRef Full Text | Google Scholar

Deng, Y., Zhai, K., Xie, Z., Yang, D., Zhu, X., Liu, J., et al. (2017). Epigenetic regulation of antagonistic receptors confers rice blast resistance with yield balance. Science 355, 962–965. PMID 28154240. doi:10.1126/science.aai8898

PubMed Abstract | CrossRef Full Text | Google Scholar

Eguiarte, L. E. (1998). Principles of population genetics by D. L. Hard and A. G. Clark, sinauer associates inc. Publishers, sunderland, MA, 1997. $58.95 (casebound), xiii + 542 pp. ISBN 0-87893-306-9. Bull. Math. Biol. 60, 1201–1205. doi:10.1016/s0092-8240(98)90008-8

CrossRef Full Text | Google Scholar

Flegel, T. W. (2012). Historic emergence, impact and current status of shrimp pathogens in Asia. J. Invertebr. Pathol. 110, 166–173. PMID 22429834. doi:10.1016/j.jip.2012.03.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Fuji, K., Hasegawa, O., Honda, K., Kumasaka, K., Sakamoto, T., and Okamoto, N. (2007). Marker-assisted breeding of a lymphocystis disease-resistant Japanese flounder (Paralichthys olivaceus). Aquaculture 272, 291–295. doi:10.1016/j.aquaculture.2007.07.210

CrossRef Full Text | Google Scholar

Fuji, K., Kobayashi, K., Hasegawa, O., Coimbra, M. R. M., Sakamoto, T., and Okamoto, N. (2006). Identification of a single major genetic locus controlling the resistance to lymphocystis disease in Japanese flounder (Paralichthys olivaceus). Aquaculture 254, 203–210. doi:10.1016/j.aquaculture.2005.11.024

CrossRef Full Text | Google Scholar

Fujimori, S., Washio, T., Higo, K., Ohtomo, Y., Murakami, K., Matsubara, K., et al. (2003). A novel feature of microsatellites in plants: A distribution gradient along the direction of transcription. FEBS Lett. 554, 17–22. PMID 14596907. doi:10.1016/s0014-5793(03)01041-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Gonen, S., Baranski, M., Thorland, I., Norris, A., Grove, H., Arnesen, P., et al. (2015). Mapping and validation of a major QTL affecting resistance to pancreas disease (salmonid alphavirus) in Atlantic salmon (Salmo salar). Hered. (Edinb). 115, 405–414. PMID 25990876. doi:10.1038/hdy.2015.37

CrossRef Full Text | Google Scholar

Hillestad, B., and Moghadam, H. K. (2019). Genome-wide association study of piscine myocarditis virus (PMCV) resistance in atlantic salmon (Salmo salar). J. Hered. 110, 720–726. PMID 31287547. doi:10.1093/jhered/esz040

PubMed Abstract | CrossRef Full Text | Google Scholar

Houston, R. D., Haley, C. S., Hamilton, A., Guy, D. R., Tinch, A. E., Taggart, J. B., et al. (2008). Major quantitative trait loci affect resistance to infectious pancreatic necrosis in Atlantic salmon (Salmo salar). Genetics 178, 1109–1115. PMID 18245341. doi:10.1534/genetics.107.082974

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, B., Wang, W., Ou, S., Tang, J., Li, H., Che, R., et al. (2015). Variation in NRT1.1B contributes to nitrate-use divergence between rice subspecies. Nat. Genet. 47, 834–838. PMID 26053497. doi:10.1038/ng.3337

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, X. D., Zhao, L., Zhang, H. Q., Xu, X. P., Jia, X. T., Chen, Y. H., et al. (2010). Shrimp NF-κB binds to the immediate-early gene ie1 promoter of white spot syndrome virus and upregulates its activity. Virology 406, 176–180. PMID 20684968. doi:10.1016/j.virol.2010.06.046

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, C., Li, H., Chen, Y., Chen, Y., Wang, S., Weng, S. P., et al. (2015a). Activation of Vago by interferon regulatory factor (IRF) suggests an interferon system-like antiviral mechanism in shrimp. Sci. Rep. 5, 15078. PMID 26459861. doi:10.1038/srep15078

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, C., Li, H., Wang, S., Song, X., Zhang, Z., Qian, Z., et al. (2015b). The c-Fos and c-Jun from Litopenaeus vannamei play opposite roles in Vibrio parahaemolyticus and white spot syndrome virus infection. Dev. Comp. Immunol. 52, 26–36. PMID 25912357. doi:10.1016/j.dci.2015.04.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, C., Wang, S., and He, J. (2019a). The two NF-κB pathways regulating bacterial and WSSV infection of shrimp. Front. Immunol. 10, 1785. PMID 31417561. doi:10.3389/fimmu.2019.01785

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, C., Weng, S., and He, J. (2019b). WSSV-host interaction: Host response and immune evasion. Fish. Shellfish Immunol. 84, 558–571. PMID 30352263. doi:10.1016/j.fsi.2018.10.043

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, C. Z., Chen, Y. X., Weng, S. P., Li, S. D., Zuo, H. L., Yu, X. Q., et al. (2014). Presence of Tube isoforms in Litopenaeus vannamei suggests various regulatory patterns of signal transduction in invertebrate NF-kappa B pathway. Dev. Comp. Immunol. 42, 174–185. PMID 24012725. doi:10.1016/j.dci.2013.08.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, H. Y., Wang, S., Chen, Y. G., Lu, K., Yin, B., Li, S. D., et al. (2017). Identification of two p53 isoforms from Litopenaeus vannamei and their interaction with NF-kappa B to induce distinct immune response. Sci. Rep. 7, 45821. PMID 28361937. doi:10.1038/srep45821

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, H. Y., Wang, S., Qian, Z., Wu, Z. Z., Lu, K., Weng, S. P., et al. (2016). MKK6 from Pacific white shrimp Litopenaeus vannamei is responsive to bacterial and WSSV infection. Mol. Immunol. 70, 72–83. PMID 26745320. doi:10.1016/j.molimm.2015.12.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Y. C., Korol, A. B., Fahima, T., Beiles, A., and Nevo, E. (2002). Microsatellites: Genomic distribution, putative functions and mutational mechanisms: A review. Mol. Ecol. 11, 2453–2465. PMID 12453231. doi:10.1046/j.1365-294x.2002.01643.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Liao, X., Wang, C., Wang, B., Qin, H., Hu, S., Wang, P., et al. (2020). Comparative transcriptome analysis of Litopenaeus vannamei reveals that triosephosphate isomerase-like genes play an important role during decapod iridescent virus 1 infection. Front. Immunol. 11, 1904. PMID 32983114. doi:10.3389/fimmu.2020.01904

PubMed Abstract | CrossRef Full Text | Google Scholar

Litt, M., and Luty, J. A. (1989). A hypervariable microsatellite revealed by in vitro amplification of a dinucleotide repeat within the cardiac muscle actin gene. Am. J. Hum. Genet. 44, 397–401. PMID 2563634.

PubMed Abstract | Google Scholar

Liu, S., Vallejo, R. L., Evenhuis, J. P., Martin, K. E., Hamilton, A., Gao, G., et al. (2018). Retrospective evaluation of marker-assisted selection for resistance to bacterial cold water disease in three generations of a commercial rainbow trout breeding population. Front. Genet. 9, 286. PMID 30123238. doi:10.3389/fgene.2018.00286

PubMed Abstract | CrossRef Full Text | Google Scholar

Lu, X., Zhou, H., Pan, Y. B., Chen, C. Y., Zhu, J. R., Chen, P. H., et al. (2015). Segregation analysis of microsatellite (SSR) markers in sugarcane polyploids. Genet. Mol. Res. 14, 18384–18395. PMID 26782486. doi:10.4238/2015.December.23.26

PubMed Abstract | CrossRef Full Text | Google Scholar

Ma, Y., Dai, X., Xu, Y., Luo, W., Zheng, X., Zeng, D., et al. (2015). COLD1 confers chilling tolerance in rice. Cell. 160, 1209–1221. PMID 25728666. doi:10.1016/j.cell.2015.01.046

PubMed Abstract | CrossRef Full Text | Google Scholar

Miah, G., Rafii, M. Y., Ismail, M. R., Puteh, A. B., Rahim, H. A., Islam, K. h. N., et al. (2013). A review of microsatellite markers and their applications in rice breeding programs to improve blast disease resistance. Int. J. Mol. Sci. 14, 22499–22528. PMID 24240810. doi:10.3390/ijms141122499

PubMed Abstract | CrossRef Full Text | Google Scholar

Moen, T., Baranski, M., Sonesson, A. K., and Kjoglum, S. (2009). Confirmation and fine-mapping of a major QTL for resistance to infectious pancreatic necrosis in atlantic salmon (Salmo salar): Population-level associations between markers and trait. BMC Genomics 10, 368. PMID 19664221. doi:10.1186/1471-2164-10-368

PubMed Abstract | CrossRef Full Text | Google Scholar

Morgante, M., Hanafey, M., and Powell, W. (2002). Microsatellites are preferentially associated with nonrepetitive DNA in plant genomes. Nat. Genet. 30, 194–200. PMID 11799393. doi:10.1038/ng822

PubMed Abstract | CrossRef Full Text | Google Scholar

Mou, C. Y., Wang, Y., Zhang, Q. Y., Gao, F. X., Li, Z., Tong, J. F., et al. (2018). Differential interferon system gene expression profiles in susceptible and resistant gynogenetic clones of gibel carp challenged with herpesvirus CaHV. Dev. Comp. Immunol. 86, 52–64. PMID 29727627. doi:10.1016/j.dci.2018.04.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Powell, W., Machray, G. C., and Provan, J. (1996). Polymorphism revealed by simple sequence repeats. Trends Plant Sci. 1, 215–222. doi:10.1016/1360-1385(96)86898-1

CrossRef Full Text | Google Scholar

Qiu, W., He, J. H., Zuo, H., Niu, S., Li, C., Zhang, S., et al. (2017). Identification, characterization, and function analysis of the NF-κB repressing factor (NKRF) gene from Litopenaeus vannamei. Dev. Comp. Immunol. 76, 83–92. PMID 28564581. doi:10.1016/j.dci.2017.05.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Si, L., Chen, J., Huang, X., Gong, H., Luo, J., Hou, Q., et al. (2016). OsSPL13 controls grain size in cultivated rice. Nat. Genet. 48, 447–456. PMID 26950093. doi:10.1038/ng.3518

PubMed Abstract | CrossRef Full Text | Google Scholar

Simpson, P., and Ayyar, S. (2008). Evolution of cis-regulatory sequences in Drosophila. Adv. Genet. 61, 67–106. PMID 18282503. doi:10.1016/S0065-2660(07)00003-X

PubMed Abstract | CrossRef Full Text | Google Scholar

Slavov, G. T., Howe, G. T., Yakovlev, I., Edwards, K. J., Krutovskii, K. V., Tuskan, G. A., et al. (2004). Highly variable SSR markers in Douglas-fir: Mendelian inheritance and map locations. Theor. Appl. Genet. 108, 873–880. PMID 14625671. doi:10.1007/s00122-003-1490-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, P. H., and He, J. G. (2019). Nucleic acid sensing in invertebrate antiviral immunity. Int. Rev. Cell. Mol. Biol. 345, 287–360. PMID 30904195. doi:10.1016/bs.ircmb.2018.11.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Wright, J. M., and Bentzen, P. (1994). Microsatellites: Genetic markers for the future. Rev. Fish Biol. Fish. 4, 384–388. doi:10.1007/BF00042912

CrossRef Full Text | Google Scholar

Xue, Y., Zhong, K., Han, B., Gui, J., and Jing, H. (2018). Innovation and achievements of designer breeding by molecular modules in China. Bull. Chin. Acad. Sci. 33, 893–899. doi:10.16418/j.issn.1000-3045.2018.09.001

CrossRef Full Text | Google Scholar

Yin, B., Wang, H., Zhu, P., Weng, S., He, J., and Li, C. (2019). A polymorphic (CT)n-SSR influences the activity of the Litopenaeus vannamei IRF gene implicated in viral resistance. Front. Genet. 10, 1257. PMID 31921300. doi:10.3389/fgene.2019.01257

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, Z., Guo, C., Sutharzan, S., Li, P., Echt, C. S., Zhang, J., et al. (2014). Genome-wide analysis of tandem repeats in plants and green algae. G3 (Bethesda) 4, 67–78. PMID 24192840. doi:10.1534/g3.113.008524

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, X., Dong, Y., Zhao, J., Huang, L., Ren, X., Chen, Y., et al. (2016). Genomic survey sequencing for development and validation of single-locus SSR markers in peanut (Arachis hypogaea L). BMC Genomics 17, 420. PMID 27251557. doi:10.1186/s12864-016-2743-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Zuo, H., Weng, K., Luo, M., Yang, L., Weng, S., He, J., et al. (2020). A MicroRNA-1-mediated inhibition of the NF-κB pathway by the JAK-STAT pathway in the invertebrate Litopenaeus vannamei. J. Immunol. 204, 2918–2930. PMID 32303554. doi:10.4049/jimmunol.2000071

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: shrimp, simple sequence repeat (SSR), selective breeding, disease resistance, IRF

Citation: Yin B, Wang H, Weng S, Li S, He J and Li C (2023) A simple sequence repeats marker of disease resistance in shrimp Litopenaeus vannamei and its application in selective breeding. Front. Genet. 14:1144361. doi: 10.3389/fgene.2023.1144361

Received: 14 January 2023; Accepted: 18 April 2023;
Published: 27 July 2023.

Edited by:

Jordan Poley, Scientist at the Center for Aquaculture Technologies in Canada, Canada

Reviewed by:

Yueling Zhang, Shantou University, China
Yang Yu, Chinese Academy of Sciences (CAS), China

Copyright © 2023 Yin, Wang, Weng, Li, He and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Chaozheng Li, lichaozh@mail2.sysu.edu.cn; Jianguo He, lsshjg@mail.sysu.edu.cn

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.