Inbreeding evaluation using microsatellite confirmed inbreeding depression in growth in the Fenneropenaeus chinensis natural population

Abstract

Understanding inbreeding depressions (IBDs), the effect on the phenotypic performance of inbreeding, is of major importance for evolution and conservation genetics. Inbreeding depressions in aquatic animals were well documented in a domestic or captive population, while there is less evidence of inbreeding depression in natural populations. Chinese shrimp, Fenneropenaeus chinensis, is an important species in both aquaculture and fishery activities in China. To investigate inbreeding depression in natural populations, four Fenneropenaeus chinensis natural populations (Huanghua, Qinhuangdao, Qingdao, and Haiyang) were collected from the Bohai and Yellow seas. Microsatellite markers were used to evaluate individual inbreeding coefficients (F) of all samples. Furthermore, the effects of inbreeding on growth traits were investigated. The results showed marker-based F was continuous and ranged from 0 to 0.585, with an average of 0.191 ± 0.127, and there was no significant difference among the average F of the four populations. Regression analysis using the four populations showed inbreeding had a very significant (p < 0.01) effect on body weight. When analyzing a single population, regression coefficients were also all negative and those in Huanghua and in Qingdao were significant at the level of p < 0.05 and < 0.01, respectively. Inbreeding depressions, expressed as the percent change in body weight per 10% increase in F, were 2.75% in Huanghua, 2.22% in Qingdao, and 3.69% in all samples. This study provided a piece of rare evidence of inbreeding depression in natural populations and also guidance toward the conservation of wild Fenneropenaeus chinensis resources.

1 Introduction

Inbreeding is defined as the mating of individuals that are related by ancestry and results in the reduction of heterozygosity (Falconer and Mackay, 1996). The inbreeding coefficient (F) is a measure of an inbreeding level and can be defined as both the probability that two alleles at any given locus are identical by descent (alleles are descendants from a single ancestor) and the probable proportion of an individual’s loci containing genes that are identical by descent (Falconer and Mackay, 1996; Bourdon, 1997). Inbreeding depression (IBD) is the effect of inbreeding, measured as the reduction in mean phenotypic performance with increasing levels of inbreeding within a population (Falconer and Mackay, 1996; Lynch and Walsh, 1998). The existence of IBD has long been known, especially for fitness traits. Understanding IBD is of major importance for the evolution and conservation of genetics. These effects have been well documented in livestock species (reviewed by Leroy, 2014), and also in aquatic animals (Keys et al., 2004; Zheng et al., 2012; Luo et al., 2014; Gao et al., 2015). However, most of these study examples were carried out in the domestic or captive population, and relatively less evidence of IBD was illustrated in natural populations (Hoffman et al., 2014). The main reason is that F was very accessible in a captive population with pedigree information, in which F of any individual can be obtained by calculating its parents’ coancestry (Falconer and Mackay, 1996). However, investigating the inbreeding level and IBD in natural populations can also be of great significance, and IBD is one of the core research fields in conservation genetics (Frankham et al., 2009). An alternative approach to the pedigree method is calculating F based on molecular markers such as microsatellites (Ritland, 1996; Lynch and Ritland, 1999; Milligan, 2003; Wang, 2007). Compared to the traditional method, F based on a molecular marker can be obtained directly, without the need for pedigree information. As a result, it provided the possibility to study IBD in natural populations.

Chinese shrimp, F. chinensis, is an important species in both aquaculture and fishery activities in China. Over the past few decades, the natural population of F. chinensis was largely reliant on released shrimps to maintain its size, and their contribution to the total landings has been consistently > 90% (Wang et al., 2006). There is a general belief that a genetic threat of the loss of variation in wild populations was one of the main concerns about an artificial propagation release (Aho et al., 2006; Araki and Schmid, 2010). Also, it was inferred that continuous artificial propagation and release had lowered the level of genetic diversity of F. chinensis in Chinese stocks (Wang et al., 2006). Although a previous study demonstrated that inbreeding has a negative effect on economic traits (especially on growth) in the F. chinensis breeding population by comparing different levels of inbreeding (Luo et al., 2014), this phenomenon has not been demonstrated under natural conditions. Such a study is important both because inbreeding may affect the extinction risk in wild populations and because understanding IBD is of major importance for the conservation genetics of F. chinensis.

In the current study, multiple F. chinensis natural populations were collected from the Bohai and Yellow seas in northern China. Microsatellite markers were used to analyze these F. chinensis samples to calculate individual F. Furthermore, its effects on growth were investigated. The results obtained in this study should provide evidence of IBD in natural populations and also guidance to the conservation of F. chinensis resources.

2 Materials and methods

2.1 Experimental materials

The F. chinensis samples were collected from four locations in the Bohai and Yellow seas: Huanghua and Qinhuangdao populations in the Bohai Sea and Qingdao and Haiyang populations in the Yellow Sea (Figure 1). The sampling time in Bohai Sea was in the autumn of 2021 and that in the Yellow Sea was in the spring of 2022. The sample size in each population ranged from 75–194, and the total was 564 (Table 1). All the samples in Qingdao and Haiyang were females because males died after mating in October or November. Those in Huanghua and Qinhuangdao contained both sexes. The body weight of each individual was measured and all samples were transported to the laboratory in liquid nitrogen and stored at −80 °C until analysis.

FIGURE 1

Genomic DNA was extracted from swimming legs in all individuals using standard phenol-chloroform procedures (Sambrook et al., 1989). The primer sequence, fluorescent dyes and anneal of genotyping microsatellite are shown in Table 2. PCR thermal cycling was performed as follows: an initial denaturing at 94°C for 5 min, followed by 30 cycles including 30 s denaturing at 94°C, 30 s annealing at locus-specific temperatures, and 40 s extension at 72°C, and then with a final extension at 72°C for 5 min. The PCR products were separated by an ABI-3130 automated Genetic Analyzer (Applied Biosystems). Alleles from the microsatellite loci were sized with a GeneScanTM-500 LIZ Size Standard (Applied Biosystems) and scored using GeneMapper™ V4.1 (Applied Biosystems).

2.2 Inbreeding coefficient calculation

The genetic diversity parameters, including the number of alleles (N), observed heterozygosity (Ho), expected heterozygosity (He), and polymorphism information content (PIC) at 11 genotyping microsatellite loci were obtained using CERVUS software (Kalinowski et al., 2007). A fixation index () was used to measure genetic differentiation among populations and calculated using Arlequin software (Excoffier et al., 2005). There are multiple ways of calculating F from microsatellite genotyping data (Ritland, 1996; Lynch and Ritland, 1999; Milligan, 2003; Wang, 2007), and a triadic likelihood method (Wang, 2007) was adopted in the current study. This was because it was also proven to be the most accurate method to calculate the paired coancestry, having either the lowest root mean error or close to the smallest one in different population structures and sizes, number of loci, and alleles (Wang, 2007; Hammerly et al., 2013). The triadic likelihood method was carried out in coancestry software (Wang, 2011) to obtain F from microsatellite genotyping data.

2.3 Inbreeding effect evaluation on body weight

A linear-regression model was used to evaluate the effects of F on body weight:where y is the measured value for body weight; a is the y-intercept; b is the regression coefficient of F; and e is the residual term. The residual term in this model includes individual genetic effects, day-age effects, and error terms, all of which are unrelated to inbreeding. In particular, since the interaction effect between sex and F was not significant (p > 0.05), the analysis of sex effects in Huanghua and Qinhuangdao populations were also subsumed into the residual term.

, expressed as the percent change in phenotype per 10% increase in F, was calculated with the following equation:where a and b are y-intercept from the regression of y on F and the regression coefficient of F, respectively (see Eq. 1). All of the aforementioned statistical analysis processes were performed with corresponding functions in R software (R Development Core Team, 2013).

3 Results

3.1 Inbreeding coefficient calculation

The N, Ho, He, and PIC at each microsatellite locus are shown in Table 3, and the average N was 25.18 ± 14.68. The ranged from 0.004–0.012 (Table 4), with an average of 0.008 ± 0.003, and only that between Qinhuangdao and Qingdao was significant (p < 0.05). The marker-based F was continuous and ranged from 0 to 0.585, with an average of 0.191 ± 0.127. The average F of Qingdao and Haiyang were 0.180 and 0.185, respectively, lower than those of Qinhuangdao and Huanghua, 0.202 and 0.194, respectively (Figure 2). However, there was no significant difference among the average F of the four populations.

FIGURE 2

3.2 Inbreeding effects on body weight

Relationships between F and measured body weight in each population and all four populations are shown in Figures 3, 4, respectively. Regression analysis’ results are shown in Table 5. Regression analysis using four populations showed that inbreeding had a very significant (p < 0.01) effect on body weight. When analyzing a single population, regression coefficients were also all negative and those in Huanghua and Qingdao were significant at the level of p < 0.05 and < 0.01, respectively. However, those in Qinhuangdao and Haiyang were not significantly different from zero due to their high standard errors. In particular, when regression analysis was based on four populations, the interaction effect between the population and F was not significant. For analyses with significant regression coefficients, estimates of IBD on body weight were further calculated. Those in Huanghua and Qingdao were 2.75% and 2.22%, respectively, and when all individuals were fitted, the estimate of IBD was 3.39% (Table 5).

FIGURE 3

FIGURE 4

4 Discussion

We estimated individual F within four F. chinensis populations directly, with no requirement of pedigree records or pedigree reconstruction. This method provided an alternative approach to detecting the inbreeding level of F. chinensis natural populations. The individual F estimate results showed that there was no significant difference in average F among different populations. All the genetic differentiation between different populations was very low, and only that between Qinhuangdao and Qingdao was significant (p < 0.05). This suggested that there was no or very little genetic differentiation among these populations, which was consistent with the expected average F. Most previous studies thought of F. chinensis distributed along northern China as one population and suggested a lack of genetic structure in this species (Hwang et al., 1997; Quan et al., 2001; Cui et al., 2007). Furthermore, there was no difference in average F before (Huanghua and Qinhuangdao) and after (Qingdao and Haiyang) the overwintering migration, which suggested that there was no significant relationship between F and the survival rate during overwintering migration. However, the average individual F obtained in this study was consistently high. There are several possible reasons responsible for such a high F. The first is that successive years of artificial propagation and release may have reduced the effective population size of F. chinensis and increased the possibility of inbreeding; another possible reason is the occurrence of null alleles that might have also led to Ho in most loci to be lower than that of He.

Based on individual F calculated from microsatellites, we provided a rare piece of evidence of IBD in a natural aquatic animal population as opposed to that in a captive population. The results of the current study showed that F had a significant negative effect on body weight in F. chinensis in the wild (p < 0.01). In addition, when the samples were divided into four populations according to where they were collected, regression analysis results in two of them supporting this conclusion. Insignificant regression coefficients in another two can likely be attributed to sampling errors which are associated with both a small sample size and a small number of microsatellite loci. Moreover, the interaction effect between the population and F was not significant, which also suggested that IBD should be universal in different populations. Similar to microsatellites, single nucleotide polymorphism (SNP) is another available molecular marker to estimate individual F (Ritland, 1996; Lynch and Ritland, 1999), and F based on whole-genome high-density SNPs should obviously be more accurate. However, the high cost restricts its routine usage in inbreeding level monitoring. If there are no applicable SNP data, using a small number of microsatellites is also an alternative method because microsatellite can provide more marker information contents due to its high polymorphism (Wang, 2016).

Although IBD typically was detected in fitness-related traits (e.g., survival and reproductive traits) (Falconer and Mackay, 1996; Lynch and Walsh, 1998), there was also numerous research on IBD of growth traits in aquaculture species. In another study on F. chinensis, Luo et al. (2014) estimated the IBD of body weight as 4.16%–4.74%, which was slightly larger than those in the current study. However, it should be noted that the estimated IBD in Luo et al. (2014) was obtained by comparing several designed high inbreeding levels (F = 0.25, 0.375, and 0.50) in a breeding population. Some research has indicated fast inbreeding (mating between close relatives, such as full sibs) caused more harm than slow inbreeding (Day et al., 2003; Thodesen et al., 2005). This view might explain the difference in results between these two studies. In addition, the occurrence of null alleles may affect the accuracy of the estimated F, and the IBD scales were also affected to some extent. Thodesen et al. (2005) summarized previous research on IBD in aquaculture species and found the average IBD of slow inbreeding on growth traits ranging from 0% to 13% with an average of 2.4%, and the results in the current study were fairly consistent with previous research.

In this study, we demonstrated the phenomenon of IBD in a natural population of F. chinensis for the first time (to the best of our knowledge). Although a study has demonstrated the existence of IBD in F. chinensis breeding population (Luo et al., 2014), the finding in the current study is still of great significance. This is because IBD under laboratory and captive conditions may not be representative of that under natural conditions (Crnokrak and Roff, 1999), and IBD may vary with environmental conditions (Armbruster and Reed, 2005). Over the past few decades, the natural population of F. chinensis in China was largely reliant on released shrimps to maintain its size (Deng and Zhuang, 2000; Wang et al., 2006). Due to its strong reproduction ability, often only a few gravid females are needed to produce enough released offspring for a hatchery. Some scholars have long been concerned about the effective population size decreasing in the F. chinensis natural population (Wang et al., 2006; Song et al., 2019). Therefore, to avoid IBD, it is very important to maintain the effective population size of the releasing population for F. chinensis.

References

• 1

  AhoT.RönnJ.PiironenJ.BjörklundM. (2006). Impacts of effective population size on genetic diversity in hatchery reared Brown trout (Salmo trutta L.) populations ective population size on genetic diversity in hatchery reared Brown trout (Salmo trutta L.) populations. Aquaculture253, 244–248. 10.1016/j.aquaculture.2005.09.013

  • CrossRef
  • Google Scholar

• 2

  ArakiH.SchmidC. (2010). Is hatchery stocking a help or harm?: Evidence, limitations and future directions in ecological and genetic surveys. Aquaculture308, S2–S11. 10.1016/j.aquaculture.2010.05.036

  • CrossRef
  • Google Scholar

• 3

  ArmbrusterP.ReedD. H. (2005). Inbreeding depression in benign and stressful environments. Heredity95 (3), 235–242. 10.1038/sj.hdy.6800721

  • CrossRef
  • Google Scholar

• 4

  BourdonR. M. (1997). Understanding animal breeding. 2nd ed. New Jersey, USA: Prentice-Hall.

  • Google Scholar

• 5

  CrnokrakP.RoffD. A. (1999). Inbreeding depression in the wild. Heredity83 (3), 260–270. 10.1038/sj.hdy.6885530

  • CrossRef
  • Google Scholar

• 6

  CuiZ.LiC.JangI.ChuK. (2007). Lack of genetic differentiation in the shrimp Penaeus chinensis in the Northwestern Pacific. Biochem. Genet.45, 579–588.

  • Google Scholar

• 7

  DayS. B.BryantE. H.MeffertL. M. (2003). The influence of variable rates of inbreeding on fitness, environmental responsiveness, and evolutionary potential. Evolution57 (6), 1314–1324. 10.1111/j.0014-3820.2003.tb00339.x

  • CrossRef
  • Google Scholar

• 8

  DengJ. Y.ZhuangZ. M. (2000). The cause of recruitment variation of Penaeus chinensis in the Bohai Sea. J. Fish. Sci. China.7, 125–128.

  • Google Scholar

• 9

  ExcoffierL.LavalL. G.SchneiderS. (2005). Arlequin (version 3.0): An integrated software package for population genetics data analysis. Evol. Bioinform.1, 117693430500100–117693430500150. 10.1177/117693430500100003

  • CrossRef
  • Google Scholar

• 10

  FalconerD. S.MackayT. F. C. (1996). Introduction to Quantitative genetics. 4th ed. Harlow, Essex, England: Longman.

  • Google Scholar

• 11

  FrankhamR.BallouJ. D.BriscoeD. A. (2009). Introduction to conservation genetics. 2nd edn. Cambridge, UK: Cambridge University Press.

  • Google Scholar

• 12

  GaoB.LiuP.LiJ.WangQ.HanZ. (2015). Effect of inbreeding on growth and genetic diversity of Portunus trituberculatus based on the full-sibling inbreeding families. Aquac. Int.23, 1401–1410. 10.1007/s10499-015-9892-9

  • CrossRef
  • Google Scholar

• 13

  HammerlyS. C.MorrowM. E.JohnsonJ. A. (2013). A comparison of pedigree- and DNA-based measures for identifying inbreeding depression in the critically endangered Attwater’s Prairie-chicken. Mol. Ecol.22, 5313–5328. 10.1111/mec.12482

  • CrossRef
  • Google Scholar

• 14

  HoffmanJ. I.SimpsonF.DavidP.RijksJ. M.KuikenT.ThorneM. A. S.et al (2014). High-throughput sequencing reveals inbreeding depression in a natural population. Proc. Natl. Acad. Sci.111 (10), 3775–3780. 10.1073/pnas.1318945111

  • CrossRef
  • Google Scholar

• 15

  HwangG. L.LeeY. C.ChangC. S. (1997). Mitochondrial DNA analysis of the fleshy prawn (Penaeus chinensis) for stock discrimination in the Yellow Sea. J. Korean Fish. Soc.30, 88–94.

  • Google Scholar

• 16

  KalinowskiS. T.TaperM. L.MarshallT. C. (2007). Revising how the computer program CERVUS accommodates genotyping error increases success in paternity assignment. Mol. Ecol.16 (5), 1099–1106. 10.1111/j.1365-294X.2007.03089.x

  • CrossRef
  • Google Scholar

• 17

  KeysS. J.CrocosP. J.BurridgeC. Y.ComanG. J.DavisG. P.PrestonN. P. (2004). Comparative growth and survival of inbred and outbred Penaeus (marsupenaeus) japonicus, reared under controlled environment conditions: Indications of inbreeding depression. Aquaculture241, 151–168. 10.1016/j.aquaculture.2004.08.039

  • CrossRef
  • Google Scholar

• 18

  LeroyG. (2014). Inbreeding depression in livestock species: Review and meta-analysis. Anim. Genet.45 (5), 618–628. 10.1111/age.12178

  • CrossRef
  • Google Scholar

• 19

  LuoK.KongJ.LuanS.MengX. H.ZhangT. S.WangQ. Y. (2014). Effect of inbreeding on survival, WSSV tolerance and growth at the postlarval stage of experimental full-sibling inbred populations of the Chinese shrimp Fenneropenaeus chinensis. Aquaculture420, 32–37. 10.1016/j.aquaculture.2013.10.030

  • CrossRef
  • Google Scholar

• 20

  LynchM.RitlandK. (1999). Estimation of pairwise relatedness with molecular markers. Genetics152, 1753–1766. 10.1093/genetics/152.4.1753

  • CrossRef
  • Google Scholar

• 21

  LynchM.WalshB. (1998). Genetics and analysis of Quantitative traits. Sunderland, MA: Sinauer Associates, 360–361.

  • Google Scholar

• 22

  MilliganB. G. (2003). Maximum-likelihood estimation of relatedness. Genetics163, 1153–1167. 10.1093/genetics/163.3.1153

  • CrossRef
  • Google Scholar

• 23

  QuanJ. X.LuX. M.ZhuangZ. M.DaiJ. X.DengJ. Y.ZhangY. P. (2001). Low genetic variation of Penaeus chinensis as revealed by mitochondrial COI and 16S rRNA gene sequences. Biochem. Genet.39, 279–284. 10.1023/a:1010234816514

  • CrossRef
  • Google Scholar

• 24

  R Development Core Team (2013). R: A Language and environment for statistical Computing. Vienna, Austria: R Foundation for Statistical Computing.

  • Google Scholar

• 25

  RitlandK. (1996). Estimators for pairwise relatedness and individual inbreeding coefficients. Genet. Res.67, 175–185. 10.1017/s0016672300033620

  • CrossRef
  • Google Scholar

• 26

  SambrookJ.FristchE. E.ManiatisT. (1989). Molecular Cloning: A laboratory manual. New York: Cold Spring Harbor Laboratory.

  • Google Scholar

• 27

  SongN.LiP.ZhangX.GaoT. (2019). Genetic diversity status and conservation priority of hatchery-produced offspring populations of Fenneropenaeus chinensis by microsatellite DNA. Mar. Freshw. Res.71 (12), 1592–1601. 10.1071/mf19235

  • CrossRef
  • Google Scholar

• 28

  ThodesenJ.HuH. L.KongJ. (2005). Inbreeding and its impact on aquaculture. J. Fish. China29 (6), 849–856.

  • Google Scholar

• 29

  WangJ. (2011). Coancestry: A program for simulating, estimating and analysing relatedness and inbreeding coefficients. Mol. Ecol. Resour.11, 141–145. 10.1111/j.1755-0998.2010.02885.x

  • CrossRef
  • Google Scholar

• 30

  WangJ. (2016). Pedigrees or markers: Which are better in estimating relatedness and inbreeding coefficient?Theor. Popul. Biol.107, 4–13. 10.1016/j.tpb.2015.08.006

  • CrossRef
  • Google Scholar

• 31

  WangJ. (2007). Triadic IBD coefficients and applications to estimating pairwise relatedness. Genet. Res.89, 135–153. 10.1017/S0016672307008798

  • CrossRef
  • Google Scholar

• 32

  WangQ.ZhuangZ.DengJ.YeY. (2006). Stock enhancement and translocation of the shrimp Penaeus chinensis in China. Fish. Res.80 (1), 67–79. 10.1016/j.fishres.2006.03.015

  • CrossRef
  • Google Scholar

• 33

  ZhengH.LiL.ZhangG. (2012). Inbreeding depression for fitness-related traits and purging the genetic load in the hermaphroditic bay scallop Argopecten irradians irradians (Mollusca: Bivalvia). Aquaculture366-367, 27–33. 10.1016/j.aquaculture.2012.08.029

  • CrossRef
  • Google Scholar