Genetic diversity of the spotted seal (Phoca largha) in Liaodong Gulf following long-term conservation activities

doi:10.21203/rs.3.rs-1749397/v1

Download PDF

Research Article

Genetic diversity of the spotted seal (Phoca largha) in Liaodong Gulf following long-term conservation activities

https://doi.org/10.21203/rs.3.rs-1749397/v1

This work is licensed under a CC BY 4.0 License

Version 1

posted

You are reading this latest preprint version

The spotted seal (Phcoa largha) is the only pinniped species reproducing in the wild in China. Wild populations of this species have declined sharply in China due to anthropogenic impacts. Forty-six individuals died accidentally in Liaodong Gulf between 2018 to 2020 and were collected. The genetic diversity and population structure of this population were investigated using mitochondrial ND4L and D-loop sequences, and the results were compared with the same population in 2005. The haplotype diversity of Liaodong Gulf spotted seals was low compared to other spotted seal populations in other habitats. Fortunately, the genetic diversity of the current Liaodong Gulf spotted seal population was higher than 15 years ago, and the population expanded between 2005 and 2020. Our findings indicated that conservation activities for Liaodong Gulf spotted seals in recent decades have been valuable, and should persist in the future.

Spotted seal Phoca largha

Liaodong Gulf

ND4L gene

D-loop region

population expansion

The spotted seal (Phoca largha) is a small-bodied pinniped widespread in cold sea areas in the northern hemisphere (Jefferson et al., 2007). It is the only pinniped species found in China that reproduces in the wild (Rugh et al., 1997). The Liaodong Gulf in China, one of eight sanctuaries for spotted seals, is the southernmost protected area. Between November and May every year, spotted seals assemble in large groups in the Liaodong Gulf to mate and reproduce (Wang 1986). Due to the destruction of their habitat by anthropogenic impacts, P. largha has been listed as one of the most endangered species in China (Gao et al., 2013). Since the 1980s, local scientific and conservation organisations have worked to protect the spotted seal population in Liaodong Bay, through the establishment of the Dalian Spotted Seal Reserve of China in 1992. Since the independent colony of spotted seals in the Liaodong Gulf, they have been isolated from the Hokkaido and Okhotsk spotted seal populations and formed an independent subspecies (Han et al., 2007). The surviving population declined to about 1,000 individuals by 2005, even though some protection measures had been put in place and seal hunting was prohibited in this region for the past two decades (Gao et al., 2015).

The genetic diversity of a population is an important indicator for assessing the environmental adaptability of individuals within the population (Markert et al., 2010). Our previous study indicated that the population of Liaodong Gulf spotted seals had relatively low genetic variation compared with those from the coast of Hokkaido in Japan and from the Sea of Okhotsk in Russia (Han et al., 2007). Following continuous efforts in recent years, the P. largha population in the Liaodong Gulf is currently estimated to number ~ 2,000 animals. From 2005 to the present, Liaodong Gulf spotted seals have raised approximately three or four generations. Therefore, understanding the variation in genetic diversity of Liaodong Gulf spotted seals is critical for the rational conservation of this population.

Mitochondrial DNA (mtDNA) is convenient for assessing population genetics due to the characteristics of strict maternal inheritance, high uniformity in different tissues, and rapid evolution rate (Cui et al., 2010). The displacement loop (D-loop), a noncoding sequence, is the region with the greatest variation in mtDNA (Brown et al., 1986), and can be used to study the genetic diversity and population structure of spotted seals (Stanley et al., 1996). Another mtDNA region, the ND4L gene, has also been applied to investigate the genetic diversity of spotted seal populations (Ayako et al., 2003). In our previous studies, we sequenced the D-loop region and ND4L gene of the Liaodong Gulf spotted seal population to assess genetic diversity in 2005 (Han et al., 2006; 2007). In the present study, we collected several spotted seal samples from the Liaodong Gulf from 2018−2020. The genetic diversity of this population was investigated using the ND4L gene and D-loop region in mtDNA as molecular markers, and the results were compared with those of our previous studies. The findings illuminated variations in genetic diversity and will promote the conservation of this spotted seal population.

Sample collection and DNA extraction

Forty-six spotted seals were obtained in the Liaodong Gulf from 2018 to 2020, all killed accidentally and found by fishermen. They were transported to the laboratory under cold conditions and abdominal muscles were collected and stored at -80°C until DNA extraction. Genomic DNA was isolated from the body wall using a TaKaRa MiniBEST Universal Genomic DNA Extraction Kit (TaKaRa, Dalian, China) according to the manufacturer’s instructions. Agarose gel electrophoresis (1%) was performed to confirm DNA extraction, and DNA concentration and purity were assessed by a NanoPhotometer Classic Launched instrument (IMPLEN, Germany). All DNA samples were stored at -20°C for further application.

PCR amplification and sequencing

The ND4L gene and D-loop region were amplified by PCR using specific primers (Table 1). PCR amplifications were carried out in an ABI2720 Thermal Cycler (Applied Biosystems, USA) in 20 µL reactions containing 10 µL of 2 × Taq Master Mix (Takara, Dalian, China), 0.5 µL of each primer, and ~ 50 ng template DNA. Thermal cycling parameters are listed in Table 2. A negative control was included in each round of PCR to check contamination, and all negative controls failed to yield products, as expected. PCR products were analysed by electrophoresis on a 1.5% agarose gel. The main band was purified and recovered using a QIAquick Gel Extraction Kit (Qiagen, Germany). Purified PCR products were sequenced with an ABI 3730 XL automatic sequencer (Perkin-Elmer, Waltham, MA, USA).

Table 1

The primer information and thermal cycling for *ND4L* gene and D-loop region.
Target	Primer information		Thermal cycling
N4DL gene	ND4L-F	CTTCCATGAGCATCGCACACAGA	94 ℃ 3 min; 35 cycles for 94 ℃ 15 s, 57 ℃ 20 s, and 72 ℃ 90 s; 72 ℃ 10 min.
N4DL gene	ND4L-R	GGCTTATGCAATTGTCACCGAGT
D-loop region	H34	CCAAATGCATGACACCACAG	94 ℃ 3 min; 35 cycles for 94 ℃ 15 s, 56 ℃ 45 s, and 72 ℃ 60 s; 72 ℃ 5 min.
D-loop region	L1624	TACACTGGTCTTGTAAACC

Table 2

Genetic diversity parameters for *P. largha* populations in Liaodong Gulf between 2005 and 2020 based on *ND4L* or D-loop sequences.
Population	Numbers	h	S	Hd	Pi	Tajima’s D	Fu’s Fs
2005-ND4L	15	2	2	0.248	0.00079	-0.51271	1.188
2020-ND4L	46	2	2	0.043	0.00008	-1.11103	-1.567
2005-D-loop	46	10	10	0.632	0.00217	-1.37483	-4.594
2020-D-loop	46	20	37	0.866	0.00561	-2.17225*	-10.349*
*Significant at level (p < 0.05).

Data analysis

Forty-six ND4L and D-loop sequences were obtained in the present study. All unique sequences of both markers have been deposited in GenBank under accession numbers MZ505010-MZ505033 Information about the genetic diversity of the Liaodong Gulf spotted seal population in 2005 was obtained from our previous studies (GenBank accession numbers AB244723-AB244729, DQ153234-DQ153242, and DQ244045-DQ244052) (Han et al., 2006; 2007). All sequences were aligned and manually corrected using ClustalX v1.83 software (Thompson et al., 1997) under default settings. The number of haplotypes (h), number of polymorphic sites (S), haplotype diversity (Hd), and nucleotide diversity (Pi) were calculated using DnaSP v6.12 (Librado and Rozas 2009). Meanwhile, Fu’s Fs statistics (Fu, 1997) and Tajima’s D tests (Tajima 1989) were also performed by DnaSP v6.12 to assess the neutrality of P. largha. Analysis of molecular variance (AMOVA) and the fixation index (Fst) were employed to estimate the genetic differentiation between the Liaodong Gulf spotted seal populations in 2005 and 2020 using Arlequin v3.5.2.2 software (Excoffier et al., 2005). Phylogenetic trees based on the haplotypes of D-loop sequences were constructed using the maximum likelihood method by MEGA v7.0 with an appropriate substitution model sequence chosen by Modeltest v3.7 (Posada and Crandall 1998). The robustness of the phylogenetic results was tested by bootstrap analysis with 1000 replicates.

Variations in the genetic diversity of Liaodong Gulf spotted seals

For the ND4L gene, one polymorphic site and two haplotypes were detected in both populations of Liaodong Gulf spotted seals in 2005 and 2020 (Table 2). Two ND4L haplotypes detected in past and present populations were completely consistent (Fig. 1). The ND4L_1 haplotype was present in 13 of 15 and 45 of 46 samples in 2005 and 2020 populations, respectively. In contrast, the ND4L_2 haplotype was only detected in 2 and one individual in the populations in 2005 and 2020, respectively. Moreover, lower values of Hd and Pi based on the ND4L gene were observed for the 2020 population (0.043 and 0.00008 for Hd and Pi, respectively) compared with 2005 (0.248 and 0.00079 for Hd and Pi, respectively; Table 2).

For the D-loop region, 10 polymorphic sites and 10 haplotypes were identified in the 2005 population, compared with 37 polymorphic sites and 20 haplotypes in the 2020 population (Table 2). Among these haplotypes, five and 14 haplotypes were only detected in the 2005 and 2020 populations, respectively (Fig. 2). DL6 and DL3 were the most abundant haplotypes of the 2005 and 2020 populations (27 and 15 individuals, respectively; Fig. 2). Based on the D-loop sequences, the Hd and Pi of total samples were 0.632 and 0.00217 for 2005 and 0.866 and 0.00561 for 2020 populations, respectively (Table 2).

Variations in population genetic structure

AMOVA results indicated that 91.64% and 91.82% of total genetic variance based on the ND4L and D-loop sequences could be attributed to variations within populations (Tables 3 and 4). The fixation index Fst values between Liaodong Gulf spotted seal populations in 2005 and 2020 based on the ND4L gene and D-loop region were 0.08364 and 0.08179, respectively, and differences were significant for the D-loop region.

Tajima’s D and Fu’s Fs tests were performed based on the ND4L and D-loop sequences. For all samples, both ND4L and D-loop neutrality tests yielded negative Tajima’s D values, among which the results for the 2020 population based on the D-loop region were significant (Table 2). Similarly, for Fu’s Fs tests, significantly negative values were also only obtained for the 2020 population based on the D-loop region (Table 2).

In the present study, two and 20 haplotypes of ND4L genes and D-loop regions were respectively identified in the Liaodong Gulf spotted seal population in 2020. In the present study, the genetic diversity indices based on the ND4L gene of spotted seal population in 2020 showed slight lower than those in 2005. However, we speculated that the lower genetic diversity based on the ND4L gene could be a mask due to the increase of detected individuals in 2020 (46 vs 17). The ND4L gene is a conserved region in mtDNA that can be applied to evaluate the phylogeny of different species belonging to the same family or order (Yu et al., 1999). Meanwhile, the ND4L genes were completely consistent between the Liaodong Gulf spotted seal populations in 2005 and 2020. Thus, we believed that the genetic diversity of Liaodong Gulf spotted seal population based on the ND4L gene remained stable between 2005 and 2020.

In contrast to ND4L gene, the D-loop region, which exhibits the greatest variation among mtDNA regions, is often used to assess differences in genetic diversity among different populations of a given species (Zhang et al., 2021). In 2003, a study investigated the genetic diversity of spotted seal population in the seas of Japan and Okhotsk based on the D-loop region, and 57 haplotypes were found in 66 individuals (Ayako et al., 2003). This implied that the current genetic diversity in the Liaodong Gulf spotted seal population remains at a relatively low level (20 of 46 vs. 57 of 66). There are some potential reasons that may explain the lower genetic diversity of Liaodong Gulf spotted seals. Particularly, the Liaodong Gulf is a narrow habitat located at the margin of the spotted seal living area, resulting in reproductive isolation from other populations. In addition, severe environmental pollution and intensive fishing activities in the Liaodong Gulf have resulted in too few successfully breeding adults in the existing population.

Although the genetic diversity of D-loop region in the Liaodong Gulf spotted seal population was relatively low, the obvious increase of it in 2020 compared to 2005 was observed in the present study. For the obtained haplotypes of D-loop region, the unique haplotypes emerged in 2020 almost detected in only one individual. In contrast, the disappeared haplotypes in 2020 were detected in several individuals in 2005. This could be the reason for the higher genetic diversity indices based on the D-loop region in the Liaodong Gulf spotted seal population in 2020 compared to 2005. The Fst value is an effective index for assessing differentiation among populations (Wright 1972). The significant Fst value for D-loop sequences between the populations investigated in the present study indicated appreciable genetic differentiation between the Liaodong Gulf spotted seal populations in 2005 and 2020. In addition, neutrality tests, such as Tajima’s D and Fu’s Fs tests, were used to examine recent population expansion when the null hypothesis was rejected (Zhang et al., 2017). Negative and significant neutrality test values indicated that sequences contained more nucleotide changes than would be predicted based on the neutral evolution model, which may suggest a population expansion event in recent history (Zhang et al., 2017). Herein, both Tajima’s D and Fu’s Fs values based on D-loop sequences in the 2020 population were significantly negative. In addition, the actual individual number of the Liaodong Gulf spotted seal population was increased from less than 1000 in 2005 to more than 2000 in 2020. These results further indicated that a population expansion event could have occurred recently.

It is well-known that the spotted seal population in the Liaodong Gulf has rapidly declined in the past few decades. Systematic conservation activities have been conducted, and the population of Liaodong Gulf spotted seals has increased, genetic diversity is lower than in other spotted seal populations. Fortunately, the genetic diversity of the Liaodong Gulf spotted seal population based on the D-loop region in now higher than it was in 2005. Our findings confirmed that conservation activities for Liaodong Gulf spotted seals are valuable and should be continued in the future.

Acknowledgements This work was supported by the China Environment and Zoology Protection for Offshore Oil and Ocean Foundation(CF-MEEC/ER/2021-15 and CF-MEEC/TR/2020-04), and the Foundation of Liaoning Province Department of Ocean and Fisheries (201812 and 201822).

Ethic approval Samples collection for spotted seals collected in Liaodong Bay, China, were authorized under the Ministry of Agriculture and Rural Affairs of the People's Republic of China, permit number: 1376. This study was conducted under a permit issued by Liaoning Fisheries Administration Bureau, Liaoning Province, China (approval number: LSYXFZ20111105).

Data Availability Statement All sequences have been deposited in GenBank at https://www.ncbi.nlm.nih.gov/nuccore/ (accession numbers MZ505010-MZ505033). Information about the genetic diversity of the Liaodong Gulf spotted seal population in 2005 was obtained from our previous studies (GenBank, https://www.ncbi.nlm.nih.gov/nuccore/, accession numbers AB244723-AB244729, DQ153234-DQ153242, and DQ244045-DQ244052).

Ayako WM, Manabu O, Manami T, Ohtaishi N (2003) Population genetic structure of the spotted seals Phoca largha along the coast of Hokkaido, based on Mitochondrial DNA sequences. Zool Sci 20:783–788
Brown GG, Gadaleta G, Pepe G, Saccone C, Sbisa E (1986) Structural conservation and variation in the D-loop-containing region of vertebrate mitochondrial DNA. J Mol Biol 192(3):503–511
Cui Z, Liu Y, Chu KH (2010) Broader pattern of tandem repeats in the mitochondrial control region of Perciformes. Chin J Ocean Limnol 28(4):785–794
Excoffier L, Laval G, Schneider S (2005) Arlequin (version 3.0): an integrated software package for population genetics data analysis. Evol Bioinform 1:47–50
Fu YX (1997) Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics 147(2):915–925
Gao XG, Han JB, Lu ZC, Zhang PJ, He CB (2013) Denovo assembly and characterization of spotted seal Phoca largha transcriptome using Illumina paired-end sequencing. Comp Biochem Phys D 8(2):103–110
Gao XG, Han JB, Lu ZC, Zhang PJ, He CB (2015) Short communication sequence variation and gene duplication at the MHC DRB loci of the spotted seal Phoca largha. Genet Mol Res 14:2055–2062
Han JB, He CB, Wang Q, Ma ZQ, Xu XH (2006) Sequence Analysis of Mitochondrial ND4, tRNAArg, ND4L and ND3 from spotted seal (Phoca largha) in Liaodong Gulf. Fish Sci 25:500–504
Han JB, He CB, Wang XM, Wang Q, Wang PL (2007) Sequence analysis of mitochondrial tRNAThr, tRNAPro and control region from spotted seals (Phoca largha) in Liaodong Gulf. Fish Sci 26:74–78
Jefferson TA, Webber MA, Pitman RL (2007) Marine mammals of the world: a comprehensive guide to their identification. Academic Press, UK
Librado P, Rozas J (2009) DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics 25(11):1451–1452
Markert JA, Champlin DM, Gutjahrgobell RE, Grear JS, Kuhn A, Mcgreevy TG, Roth AC, Bagley MJ, Nacci DE (2010) Population genetic diversity and fitness in multiple environments. BMC Evol Biol 10:205–213
Posada D, Crandall KA (1998) Modeltest: testing the model of DNA substitution. Bioinformatics 14:817–818
Rugh DJ, Shelden KE, Withrow DE (1997) Spotted seals, Phoca largha, in Alaska. Mar Fisheries Rev 59:1–18
Stanley HF, Casey S, Comahan JM, Goodman S, Harwood J, Wayne RK (1996) Worldwide patterns of mitochondrial DNA differentiation in the Harbor Seal (Phoca vitulina). Mol Biol Evol 13:368–382
Tajima F (1989) Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics 123:585–595
Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG (1997) The Clustal_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res 25:4876–4882
Wang P (1986) Distribution, ecology and resource conservation of the spotted seal in the Huanghai and Bohai Seas. Acta Oceanol Sin 1:1–14
Wright S (1972) Evolution and the genetics of populations. J Bio Sci 4:253–256
Yu H, Wang W, Fang S, Zhang YP, Lin FJ, Geng ZC (1999) Phylogeny and evolution of the Drosophila nasuta subgroup based on mitochondrial ND4 and ND4L gene sequences. Mol Phylogenet Evol 13(3):556–565
Zhang R, Leng P, Wang X, Zhang Z (2017) Molecular analysis and genetic diversity of Aedes albopictus (Diptera, Culicidae) from China. Mitochondr DNA A 13:1–6
Zhang S, Li M, Sun Y, Shang H, Wang L, Yang T, Ma L, Chen Y, Zhang B, Liu T, Chen W (2021) Analysis of the relationship between geography and body color with the genetic diversity in the Echiura worm Urechis unicinctus based on the mitochondrial COI and D-loop sequences. Mitochondr DNA B 6:1380–1386

No competing interests reported.

Download PDF

Version 1

posted

You are reading this latest preprint version

Genetic diversity of the spotted seal (Phoca largha) in Liaodong Gulf following long-term conservation activities

Status:

Version 1

Abstract

Figures

Introduction

Materials And Methods

Sample collection and DNA extraction

PCR amplification and sequencing

Data analysis

Results

Variations in the genetic diversity of Liaodong Gulf spotted seals

Variations in population genetic structure

Discussion

Conclusion

Declarations

References

Additional Declarations

Status:

Version 1