• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Next-generation Sequencing of MHC Class I Genes Reveals Trans-species Polymorphism in Eutropis multifasciata and Other Species of Scincidae

    2024-01-02 06:49:22ShufangZHANGYoufuLINYingzhiCHENGHaiyunYANGXiamingZHUYuDULonghuiLINYanfuQULianCHENandHongLI
    Asian Herpetological Research 2023年4期

    Shufang ZHANG ,Youfu LIN ,Yingzhi CHENG ,Haiyun YANG ,Xiaming ZHU ,Yu DU,3 ,Longhui LIN,Yanfu QU,Lian CHEN* and Hong LI*

    1 College of Life Sciences,Nanjing Normal University,Nanjing 210023,Jiangsu,China

    2 College of Life Sciences,Nanjing Forestry University,Nanjing 210037,Jiangsu,China

    3 Hainan Key Laboratory of Herpetological Research,College of Fisheries and Life Science,Hainan Tropical Ocean University,Sanya 572022,Hainan,China

    4 Herpetological Research Center,College of Life and Environmental Sciences,Hangzhou Normal University,Hangzhou 311121,Zhejiang,China

    Abstract The genes of the major histocompatibility complex(MHC)encode cell surface proteins that are essential for adaptive immunity.MHC genes show the most prominent genetic diversity in vertebrates,reflecting the adaptation of populations to their evolving environment,population survival and reproduction.In the present study,we used nextgeneration sequencing (NGS) to study the loci polymorphism of exon 3 of the MHC class I genes in an ovoviviparous skink,the many-lined sun skink,Eutropis multifasciata and five other species of Scincidae,to quantify genetic variation.In addition,we genotyped the same MHC class I genes of E.multifasciata using clone sequencing,to directly compare the effectiveness of both analytical techniques for MHC genotyping.NGS detected 20 MHC class I alleles in E. multifasciata,and 2 to 15 alleles in the other five Scincidae species.However,clone sequencing detected only 15 of those MHC class I alleles in E. multifasciata.In addition,transspecies polymorphism of MHC class I genes was studied by constructing a phylogenetic tree using the gene sequences obtained by NGS.Phylogenetic analysis revealed that MHC class I alleles were shared among different species of Scincidae with trans-species polymorphism,and did not exhibit specific genealogical inheritance.These results have important implications for understanding polymorphism interspecies diversity in the MHC genes of Scincidae,and the evolution of the MHC more broadly.

    Keywords Eutropis multifasciata,major histocompatibility complex,next-generation sequencing,Scincidae,trans-species polymorphism

    1.Introduction

    The immune system is a complex network of organs,cells and proteins which play vital roles in protecting an organism from pathogens (Elbers and Taylor,2016).The major histocompatibility complex(MHC)is part of the sophisticated“adaptive immune system”;it encodes the molecules that are responsible for antigen recognition,and subsequent binding(Sommer,2005).The MHC is divided into three classes(MHC class I,II and III) based on the location and function of the genes within it.MHC class I genes primarily encode regions which identify and respond to polypeptides of endogenous antigens (i.e.,antigens that have been newly synthesized in antigen-presenting cells)and initiate the subsequent immune response (Neefjeset al.,2011).In order to bind to different antigenic peptides,the antigen-binding regions(exons 2 and 3)of MHC class I must be highly polymorphic.Polymorphisms in MHC class I genes have been studied across the Cercopithecinae (Caskeyet al.,2019),Lacertidae (Yuanet al.,2014) and Muscicapidae (Biedrzyckaet al.,2017),generally using exon 3 as the typing target.Perhaps not surprisingly,individuals with high polymorphisms in MHC class I genes have stronger immune function in response to a wide range of pathogens(Radwanet al.,2020).

    Within the jawed vertebrate genome,the evolutionary history of the MHC has varied across clades,resulting in the MHC consisting of a highly polymorphic,multigene family.Polymorphisms of MHC genes are created by individual adaptation in a population,resulting in the occurrence of multiple alleles for a single gene locus.As such,MHC genes are important markers for species adaptations,like developing resistance to parasites(Hackinget al.,2018)or acclimating to different geographical locations(Sagonaset al.,2019).

    Another form of polymorphism in MHC genes is transspecies polymorphism (TSP),which is the occurrence of identical or similar alleles among closely related species,that have not resulted from convergent evolution or gradual crossfertilization (Kleinet al.,2007).In general,MHC genes are homologous among species,have high sequence polymorphism,and linkage disequilibrium among different loci.For example,using exon 3 of the MHC class I genes researchers reconstructed the phylogenetic tree of three

    Eremiasskink species (E.przewalskii,E.brenchleyi,andE.multiocellata).The resulting tree showed that the MHC class I gene lineages were not clustered by species,rather,they exhibited TSP across the three species (Yuanet al.,2014).In such cases,MHC alleles could be older than the species and would thus have sufficient time to accumulate multiple interallelic variations,resulting in the development of the TSP(Kleinet al.,2007).

    The genes characteristics of the MHC are well described in mammals(especially humans and mice),birds,amphibians and fish,but little is known about their structural features in reptiles(Arnold and H?mmerling,1991;Jongsmaet al.,2019),leaving a large phylogenetic gap.To understand the wider evolution of the MHC,we must study the sequence characteristics of reptilian MHC genes.The reptilian family Scincidae is large and diverse,both morphologically and ecologically.Because of this,the Scincidae provide an ideal model system to obtain fundamental information about the evolution of MHC diversity and TSP across the wider family.The many-lined sun skink(Eutropis multifasciata)is a mediumsized,tropical lizard,native to China,India,Bangladesh and the Philippines (Duet al.,2012).It is ovoviviparous and a suitable species to investigate the proximate and ultimate reasons for variation in the reproductive patterns of reptiles(Sunet al.,2012).Populations of this skink are expected to decline dramatically in the future,due to the medicinal and culinary values of this species(Dunget al.,2011).In 2023,the State Forestry Administration of the People’s Republic of China recognized this threat and identified thatE.multifasciatawas of important economic and scientific value and recommended state protection as within China.Given this context,the study of the MHC inE.multifasciataalso provides information for the conservation of endangered species.

    When sequencing genes,researchers generally employ a protocol using the polymerase chain reaction(PCR)to amplify the corresponding alleles.Initially,the cloning and then sequencing technique uses the plasmid incompatibility principle to isolate heterozygous alleles and obtain haplotypic gene information,commonly used for the confirmation of new alleles,obtaining gene sequences and polymorphism analyses (Christet al.,1994).However,this may not be the most suitable method for sequencing MHC genes.Previous studies have found that in many species sequencing MHC genes with PCR can cause gene duplication,copy number variation and co-amplification of alleles from different loci,making the determination of individual genotypes impossible(Babik,2010).Fortunately,next-generation sequencing(NGS)is a massively parallel sequencing technology with high throughput,scalability and speed.This technology facilitates the accuracy and reproducibility of genotyping polymorphic gene families,such as the MHC,in non-model organisms(Groganet al.,2016;Shaheenuzzamnet al.,2020).

    In the present study,we use NGS to investigate exon 3 of the MHC class I genes in six skink species(E.multifasciata,P.chinensis,S.modesta,S.reevesii,S.indicusandS.incognitus).MHC genotyping is difficult due to the high degree of polymorphism between MHC gene sequences and the variation in allele numbers between species (Rekdalet al.,2018).In order to compare the ability of different techniques to detect MHC genes,we also used clone sequencing to detect exon 3 of the MHC class I genes inE.multifasciata.The results of this study have important implications for understanding polymorphism and interspecific diversity in the MHC genes of the Scincidae.

    2.Materials and Methods

    2.1.Sample collection and DNA extractionWe collected 14E.multifasciataand 8 other lizards in China(Table 1).The samples were kept in the Reptile Laboratory at Nanjing Normal University.DNA was extracted from the muscle tissues of 22 lizard samples using the DNeasy Tissue Kit(Qiagen,Germany)according to the manufacturer’s protocol.

    Table 1 Sample information of six species of Scincidae.

    2.2.MHC genotyping by an amplicon-based NGS approachThe universal primers(Murphyet al.,2009)2MHCF(5'-CAG CAG ATG TAT GGC TGT GA-3') and 2MHCR (5'-GCAGAT CTC CTC CAG GTA G-3')were used to amplify exon 3 of the MHC class I genes(corresponding to the α-2 region)in six species(E.multifasciata,Plestiodon chinensis,Scincella modesta,S.reevesii,S.indicus,andS.incognitus)and a unique barcode was added to the 5' end of the primers to identify different individuals.The PCR amplification reaction system was 50 μL,including 10 μL Premix Taq,1 μL of each primer(10 μmol/L),2 μL template DNA(100 ng)and ddH2O.The following cycle conditions were used for PCR amplification:94°C for 5 min;35 cycles of 94°C for 30 s,60°C for 30 s and 72°C for 1 min;72°C for 10 min.PCR products were mixed and purified with the QIAquick Gel Extraction Kit(Qiagen,Germany).The size and quality of the pooled amplicons were verified using the TruSeq? DNA PCR-Free Sample Preparation Kit (Illumina,USA)on the 2100 Bioanalyzer(Agilent,USA).The library was then sequenced on the Illumina NovaSeq and 250 bp pairedend sequences were obtained (Novogene Bioinformatics Technology,China).

    The raw sequencing data contained a variety of error messages that interfered with the subsequent analysis,such as sequencing junctions,low-quality bases and undetected bases(expressed as N).Therefore,the data were filtered by four criteria: (1) filtering out splice sequences and barcode sequences from the reads;(2) splicing pairs of reads with overlaps using FLASH VI.2.7(Mago?et al.,2011);(3)filtering the spliced data using QIIME(Caporasoet al.,2010)to remove chimeric sequences with high N or low quality;and (4)filtering out chimeric sequences from the spliced sequences.The Illumina data were then clustered using Vsearch(Rogneset al.,2016) and amplicons with a frequency of less than 1%were discarded.We only considered sequences that were present in more than one individual to be alleles(Zhaoet al.,2013).

    2.3.MHC genotyping by clone sequencingUniversal primers 2MHCF and 2MHCR were used to amplify exon 3 of MHC class I genes inE.multifasciata.PCR reaction procedures and systems were as described above.Amplified PCR products were then separated using agarose gel electrophoresis,and the amplified bands of 14 individuals ofE.multifasciatawere excised and cleaned (MiniElut RNA-Pure Kit,Genebase Biogene Technology).The purified PCR product was ligated to the pMD19-T vector(TaKaRa,Japan)and transformed into DH5α receptor cells.PCRs were then performed for positive clones.Each individual then had 10 colonies and the bacterial solution containing the target fragment was sent to Shanghai Sangon Biological Engineering Technology &Services Co.,Ltd.for sequencing.

    The initial quality and length filtering of raw data was performed with ChromasPro (v 2.1) to output FASTA files.The target sequences were obtained by removing redundant sequences using the MEGA X (Kumaret al.,2018).The sequences obtained from clone sequencing were imported into NCBI (https://www.ncbi.nlm.nih.gov/) for BLAST to identify MHC class I genes.The sequences were imported into MEGA X and aligned using Clustal W (Larkinet al.,2007).Haplotype information for MHC sequences was obtained by DnaSP 5.0.MHC class I alleles were named with reference to the naming principle:a new allele needed to be identified by its occurrence in at least two clones of the same individual or in different individuals (Kennedyet al.,2002) and named sequentially according to the frequency of allele occurrence.

    2.4.Phylogenetic analysisThe sequences from clone and next-generation sequencing were imported into DnaSP 5.0 to calculate the nucleotide diversity (π) and variable sites (S) of MHC class I genes.The TCS Network algorithm in Popart(Hayeset al.,2022) was used to construct separate network relationships between the diverse haplotypes of Scincidae.The exon 3 sequences of MHC class I derived from 6 species in the present study combined with other 4 species (Ctenophorus decresii,Eremias multiocellata,E.przewalskii,andE.brenchleyi)were used for phylogenetic analysis(Supplementary Table S1 for GenBank accession numbers).The sequences were aligned using MAFFT implemented in PhyloSuite(Zhanget al.,2020).The phylogenetic tree was created by the maximumlikelihood (ML) method implemented in IQ-TREE with a bootstrap of 1 000 replicates.The phylogenetic tree was visualized with the iTOL online tool(https://itol.embl.de/).

    3.Results

    In this study,the exon 3 sequences of MHC class I genes were successfully amplified in 22 samples fromE.multifasciata,P.chinensis,S.modesta,S.reevesii,S.indicus,andS.incognitus.

    3.1.Comparison of clone sequencing and NGSThe nucleotide sequences of all alleles identified by sequencing did not reveal a stop codon,indicating that these alleles all encode functional proteins.A total of 20 alleles were identified in 14 individuals.According to the naming principle of MHC class I gene sequences,20 MHC class I sequences were namedEumu-UA*01toEumu-UA*20(Table 2).Fifteen alleles were identified by clone sequencing and all of these alleles were also detected by NGS.Furthermore,clone sequencing identified samples with 3 to 7 alleles,with half of the samples having 3 alleles and only one sample having 7 alleles.The number of alleles detected by NGS ranged from 4 to 14,with four samples containing nine alleles,five samples containing more than ten alleles,and five samples containing fewer than nine alleles(Supplementary Table S2).

    Table 2 Allele frequency of MHC I gene exon 3 of E.multifasciata based on clone sequencing and NGS.

    3.2.Class MHC class I genes sequence polymorphismThere were some differences in nucleotide diversity among different species.The number of variable sites detected forE.multifasciatawas 90,with a nucleotide diversity of 0.101.No singleton variable sites were presented,but 90 parsimony informative sites were detected(Table 3).

    Table 3 Genetic diversity of MHC class I genes in 6 species of Scincidae.

    In the other five Scincidae species,the number of variable sites ranged from 39 to 97,and the nucleotide diversity was 0.103-0.232.There were no singleton variable sites in exon 3 of the MHC class I genes inS.modesta,and the value of the parsimony informative sites was 54.No parsimony informative sites were present inS.reevesii,and the value of the singleton variable sites was 39.Three skinks(P.chinensis,S.indicus,andS.incognitus)had both singleton variable sites and parsimony informative sites.A total of 15 alleles were present inS.incognitus,13 alleles were present inS.indicusand 10 alleles were present inP.chinensis.S.modestaandS.reevesiieach contained at least two allelic loci(Table 3).

    3.3.MHC class I genes trans-species polymorphismThere were eight alleles(Eumu-UA*01,Eumu-UA*02,Eumu-UA*05,Eumu-UA*06,Eumu-UA*07,Eumu-UA*08,Eumu-UA*11,andEumu-UA *14) shared with four species (E.multifasciata,S.indicus,S.incognitus,andP.chinensis).E.multifasciataandS.indicusshared 13 alleles.E.multifasciataandP.chinensisshared 10 alleles.No shared alleles were found amongE.multifasciata,S.modesta,andS.reevesii.TheE.multifasciatahad five unique alleles(Eumu-UA*09,Eumu-UA*10,Eumu-UA*15,Eumu-UA*16,andEumu-UA*18),and theS.indicusonly had one alleleSpinc-UA*01(Figure 1).

    Figure 1 Median-joining network of MHC class I exon 3 alleles found in Scincidae.Each circle represents a haplotype,and the size of the circle indicates the number of samples contained in that haplotype.The line between two circles indicates that the two haplotypes are related,and the short line on the line indicates the number of nucleobase substitutions required to change from one haplotype to another.

    The evolutionary relationship of MHC class I was explored by constructing a phylogenetic tree of the sequences of exon 3 from ten lizard species(Figure 2).It was found that the MHC class I exon 3 sequences of theEremiaslizards(E.multiocellata,E.przewalskii,andE.brenchleyi) clustered into one large branch,while those ofC.decresiiclustered into a single branch.However,the MHC class I exon 3 sequences of the six lizards in the family Scincidae were not clustered into one branch according to the family.Although the phylogenetic tree formed sub-branches of Scincidae,Lacertidae and Agamidae,it did not show obvious phylogenetic characteristics,and there was TSP in the exon 3 sequence of MHC class I genes.

    Figure 2 The Maximum Likelihood (ML) phylogenetic tree based on the MHC class I alleles from Eutropis multifasciata, Plestiodon chinensis,Scincella modesta, Scincella reevesii, Sphenomorphus indicus, Sphenomorphus incognitus, Eremias multiocellata, Eremias przewalskii, Eremias brenchleyi,and Ctenophorus decresii.

    4.Discussion

    4.1.Characteristics of MHC class I genes in E.multifasciataThe MHC plays an important role in the immune response of vertebrates.The inherent level of genetic variation and function within a species’ MHC underlies disease resistance,directly affecting the survival of individuals and even the fitness of entire populations(Miniaset al.,2016).Additionally,the level of genetic variation in MHC genes can be an important indicator for evaluating the evolutionary potential of adaptation in endangered wildlife(Akiyamaet al.,2017).Genetic diversity within the MHC is reflected in allelic and nucleotide diversity,as well as in the number of MHC genes.In the present study,amplicon-based NGS yielded a total of 20 MHC class I alleles inE.multifasciata.In comparison with other reptile families,the number of MHC class I alleles in the ornate dragon lizard(C.ornatus)was close to that ofE.multifasciata,ranging from 3-14 (Radwanet al.,2014).In contrast,the Tuatara(Sphenodon spp.)possessed 49 MHC class I alleles(Milleret al.,2010),which was significantly higher than theE.multifasciata.However,only six MHC class I alleles were found in a more closely related skink speciesC.decresii(Yuanet al.,2014),which was lower than that ofE.multifasciatain the present study.The main driver of genetic diversity within the MHC is pathogen-mediated selection,but other evolutionary processes(e.g.,genetic drift and gradual infiltration)may also increase genetic diversity within the MHC (Sagonaset al.,2019).

    4.2.Comparison of MHC typing using clone sequencing and next-generation sequencingTo understand the functional significance of allelic diversity in the MHC within the context of evolutionary ecology,pathogen resistance and conservation,it is imperative that the MHC composition of individuals is correctly assessed.Several methods have been used to type MHC genes,such as denaturing gradient gel electrophoresis (DGGE);singlestrand conformation polymorphism (SSCP);reference strand-mediated conformational analysis (RSCA);and direct sequencing after PCR purification (Lenzet al.,2008).The most widely used of these methods are clone sequencing and next-generation sequencing,which have their own advantages and limitations.The number of classical MHC genes may vary intra-or inter-species(de Grootet al.,2015;Plasilet al.,2022).The traditional method is to amplify individual samples using primers and pool the amplification products from the individual samples for sequencing.However,the MHC genes of many species exhibit gene duplications,copy number variations,and common alleles from different loci,making it impossible to determine individual genotypes using traditional methods (Babik,2010).In addition,clone sequencing can only be performed on small batches of samples,which is cumbersome and less sensitive,i.e.,susceptible to allelic deletions.Clone sequencing is also limited by the size of the cloning vector,the larger the amplified fragment,the more difficult it is to clone.NGS technology is easy to use and powerful,it detects large samples of MHC alleles concurrently during typing.However,it is difficult to obtain complete sequences and usually targets one or more exons(Chenget al.,2022).

    In the present study,the sequences of exon 3 of the MHC class I genes ofE.multifasciatawere sequenced using conventional clone and next-generation sequencing techniques,respectively.A total of 15 MHC alleles were detected by both clone sequencing,and NGS techniques,however NGS also detected an additional five new MHC alleles.This demonstrates that NGS is superior for sequencing highly polymorphic MHC genes,and can thus reduce the risk of allelic loss caused by preferential amplification of fewer alleles.For example,Eumu-UA*01was detected in only 6 individuals by clone sequencing,whereas it was detected in all 14 individuals by NGS.Therefore,NGS technology has clear advantages over traditional clone sequencing for genotyping of the MHC.This builds on previous studies utilizing NGS for MHC in other taxa.Genotyping of the MHC in Korean native chickens using long-range PCR in combination with NGS confirmed that NGS technology could be used effectively to detect MHC variants in avian fauna(Ediriweeraet al.,2022).Another study utilizing NGS to genotype the MHC of the Stokes skinkEgernia stokesiirevealed that 43 MHC I alleles were present in the offspring and two MHC I alleles were missing in its parents.This provides a useful reference for MHC genotyping in other taxa (Pearsonet al.,2017).In addition,clone sequencing and next-generation sequencing were used to type the MHC genes of 28 goat antelope,Rupicapra rupicapra,showing differences in MHC typing between the two methods in 25%of the individuals and again confirming the superior ability of next-generation sequencing to detect highly polymorphic genes(Stipoljevet al.,2020).

    However,NGS methods can be prone to sequencing errors and false positives,i.e.,false “alleles”.To minimize the occurrence of false positives,sequencing depth and coverage can be increased,thus avoiding errors generated by sequencing,and additional bioinformatics methods can be developed to address these issues(Chenget al.,2022;Sebastianet al.,2018).Therefore,the use of NGS technology to obtain MHC gene sequences should take the corresponding sequencing depth into account.In the present study,NGS was used based on the amplicon principle,i.e.,using a high sequencing depth to exclude sequencing errors (the average amplicon reads were 62 533 base pairs) to ensure genotype accuracy.In addition,the deeper the sequencing depth,the higher the probability of detecting alleles,which could compensate for the limitations of difficult allele amplification in clone sequencing (differences in amplification efficiency between alleles)(Biedrzyckaet al.,2017).

    4.3.Intraspecific and trans-species polymorphisms of MHC class I genesMHC genes form a highly polymorphic gene family with a higher rate of variable sites than other genes,and most animal MHC genes have genetic variability(Kohnet al.,2006).However,reptiles have low levels of genetic diversity in their MHC genes.The nucleotide diversity of different populations ofPodarcis erhardiiranged from 0 to 0.00564 (Hurstonet al.,2009).The nucleotide diversity of three populations in the southern Flinders Ranges ofE.stokesiiranged from 0.151 to 0.163(Pearsonet al.,2017).The π value of the MHC allele for 208Sphenodon punctatusfrom 12 island populations was 0.152±0.013.Compared to the above lizards,the π value ofE.multifasciatawas lower(0.101).The level of genetic diversity within a species can predict its evolutionary potential and viability(Guoet al.,2019).The higher the level of genetic diversity of the species,the smaller the impact of external environmental change,and the stronger the ability to survive(Paulset al.,2013).Multiple factors can influence the level of genetic diversity within a species,such as environmental variables,genetic drift within populations,bottleneck effects,inbreeding and life history traits,and diversity in a population.(Raposo Do Amaralet al.,2013).

    The six species of lizards in this study are listed as threatened under three kinds of animal protection laws for the state,andS.indicusis assessed as a near-threatened species.Because of these classifications,and the vulnerability of these species,our sample sizes were small.The sample size may have affected the nucleotide diversity of the alleles to some extent.However,results generated from this study could help us gain a simple understanding of the level of genetic diversity in a given population,with implications for the conservation and management of particular populations.

    Both the allele network and the phylogenetic tree showed that alleles were shared among the six lizard species in the family Scincidae,and they were not clustered by family in the phylogenetic tree.This indicated the existence of TSP in the MHC class I genes.This TSP could have arisen because the allelic lineage emerged before the species was formed and was retained after the species became established (Kleinet al.,1998).It could also be possible that the six species have overlapping ranges,with hybrid gradual infiltration between them(Guoet al.,2011),and that the evolutionary relationships exhibited by their MHC class I genes may be the result of hybridization.The phenomenon of TSP in MHC locus allele lineages has been reported in many vertebrates.For example,a phylogenetic tree constructed from exon 3 of the MHC class I genes of the oviparousE.brenchleyi,ovoviviparousE.multiocellataandE.przewalskiishowed that the MHC genes of the three species did not form their own genetic lineages according to species and did not show specific evolutionary lineage inheritance (Yuanet al.,2014).The phylogenetic tree showed that there was no interspecific mixing between the MHC class I alleles ofE.multiocellataand the three species of the generaEremias.It is therefore hypothesized that the divergence between the family Scincidae and Lacertidae occurred before the allelic lineage was formed.In addition,the phylogenetic relationships between the MHC I genes ofE.multiocellataand other lizards that are closely related were also studied,and alleles from each class did not fit neatly into several branches.This suggested that either the shared MHC class I genes ofE.multiocellatahad recently undergone too much gene duplication to exhibit homologous clustering of alleles,or that the alleles had not evolved independently after gene duplication(Winternitzet al.,2023).Some MHC class I alleles exhibit a greater degree of phylogenetic similarity between species than within species.This can be common for immune genes in general,which are subject to parasitemediated balancing selection and maintain allelic diversity over time(Tě?icky and Vinkler,2015).Convergent evolution,resulting from similar selective pressures,can also cause shared polymorphisms among independent evolutionary lineages(Krieneret al.,2000;Srithayakumaret al.,2012).

    5.Conclusions

    In most animals,the major histocompatibility complexes are polymorphic and polygenic,which makes genotyping problematic.This study demonstrated that when sequencing highly variable regions,such as the MHC,NGS technology performs better than clone sequencing.The level of genetic diversity of the MHC genes inE.multiocellatais low.In addition,the exon 3 sequences of the MHC gene of six species of lizards in the family Scincidae display TSP.This study contributes to our understanding of the MHC structure in lizards of the family Scincidae,and explores methods suitable for MHC genotyping in reptiles,while providing a reference for future studies on reptile genetic diversity.

    Data availabilityThe raw data presented in the study have been uploaded on the SRA database from NCBI under the BioProject PRJNA962605.

    AcknowledgementsWe would like to thank Wenjia SHEN for her lab assistance during the research.We thank Georgia WARD-FEAR and Zheng LI for their constructive comments on the revised version.This work was supported by grants from the National Natural Science Foundation of China(32171495 and 31971414),and the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD).

    国产亚洲精品一区二区www| 精品乱码久久久久久99久播| 五月伊人婷婷丁香| avwww免费| 免费电影在线观看免费观看| 欧美日韩中文字幕国产精品一区二区三区| 国产av麻豆久久久久久久| 国产 一区 欧美 日韩| 蜜桃久久精品国产亚洲av| 亚洲天堂国产精品一区在线| 99热这里只有是精品50| 99精品在免费线老司机午夜| 熟妇人妻久久中文字幕3abv| 又紧又爽又黄一区二区| 国产99白浆流出| 2021天堂中文幕一二区在线观| 男女午夜视频在线观看| 亚洲成人中文字幕在线播放| 他把我摸到了高潮在线观看| 色综合站精品国产| 母亲3免费完整高清在线观看| 高清毛片免费观看视频网站| 青草久久国产| 国内精品久久久久久久电影| 男女那种视频在线观看| 熟妇人妻久久中文字幕3abv| 日韩精品青青久久久久久| 在线观看一区二区三区| 90打野战视频偷拍视频| 亚洲五月天丁香| aaaaa片日本免费| 真实男女啪啪啪动态图| 亚洲国产日韩欧美精品在线观看 | 国产午夜福利久久久久久| 欧美又色又爽又黄视频| 精品国产三级普通话版| 久久精品影院6| 99在线人妻在线中文字幕| 国产精品爽爽va在线观看网站| 露出奶头的视频| 欧美黑人巨大hd| 可以在线观看的亚洲视频| 久久久国产成人精品二区| 制服人妻中文乱码| 亚洲无线观看免费| 网址你懂的国产日韩在线| 亚洲精品在线观看二区| 久久精品国产清高在天天线| 日韩欧美免费精品| 亚洲欧美日韩卡通动漫| 一进一出抽搐动态| 制服人妻中文乱码| 免费看光身美女| 亚洲中文字幕一区二区三区有码在线看| 亚洲色图av天堂| 老熟妇仑乱视频hdxx| 成年女人看的毛片在线观看| av天堂在线播放| 国产精品久久久久久精品电影| www.999成人在线观看| 舔av片在线| 亚洲五月天丁香| 少妇人妻精品综合一区二区 | 亚洲中文字幕日韩| 亚洲国产精品999在线| 哪里可以看免费的av片| 欧美成人一区二区免费高清观看| 久99久视频精品免费| 久久久久精品国产欧美久久久| 欧美日本视频| 亚洲精品在线美女| 亚洲乱码一区二区免费版| 搡女人真爽免费视频火全软件 | 一级黄色大片毛片| 日韩人妻高清精品专区| 色噜噜av男人的天堂激情| 精品国产三级普通话版| 男女床上黄色一级片免费看| 黄色视频,在线免费观看| 亚洲精品久久国产高清桃花| 狠狠狠狠99中文字幕| 国产亚洲av嫩草精品影院| 国产日本99.免费观看| 91麻豆av在线| 看免费av毛片| 日韩国内少妇激情av| 成年女人看的毛片在线观看| 人人妻人人澡欧美一区二区| 草草在线视频免费看| 成人特级av手机在线观看| 亚洲精品影视一区二区三区av| 18美女黄网站色大片免费观看| 久久久久国内视频| 亚洲人成伊人成综合网2020| 黄色丝袜av网址大全| 国产熟女xx| 亚洲成人免费电影在线观看| 欧美一区二区精品小视频在线| av天堂在线播放| 精品国产美女av久久久久小说| 淫秽高清视频在线观看| 久久精品91无色码中文字幕| 国产视频内射| 日韩欧美精品v在线| 亚洲七黄色美女视频| 91麻豆av在线| 国产99白浆流出| 少妇的逼好多水| 18禁黄网站禁片午夜丰满| 在线播放无遮挡| 亚洲美女黄片视频| 俺也久久电影网| 精品一区二区三区av网在线观看| 日韩中文字幕欧美一区二区| 国产一区二区在线av高清观看| 成人性生交大片免费视频hd| 看片在线看免费视频| 国产三级在线视频| 国模一区二区三区四区视频| 亚洲欧美日韩高清在线视频| 午夜影院日韩av| 国产伦人伦偷精品视频| 岛国在线观看网站| 精品电影一区二区在线| 人人妻,人人澡人人爽秒播| 国产色婷婷99| 久久国产精品影院| 亚洲激情在线av| 99国产综合亚洲精品| 无人区码免费观看不卡| 免费观看精品视频网站| 变态另类成人亚洲欧美熟女| 91av网一区二区| 天美传媒精品一区二区| 97超级碰碰碰精品色视频在线观看| 最近最新中文字幕大全电影3| 少妇高潮的动态图| 久久国产乱子伦精品免费另类| 长腿黑丝高跟| www.www免费av| 最新在线观看一区二区三区| 久久人妻av系列| netflix在线观看网站| 夜夜看夜夜爽夜夜摸| 色噜噜av男人的天堂激情| 色综合站精品国产| 男人和女人高潮做爰伦理| 欧美日韩亚洲国产一区二区在线观看| 精品国内亚洲2022精品成人| 亚洲av五月六月丁香网| 2021天堂中文幕一二区在线观| 久久精品国产清高在天天线| 亚洲七黄色美女视频| 中文字幕精品亚洲无线码一区| 精品人妻一区二区三区麻豆 | 91av网一区二区| 18美女黄网站色大片免费观看| 制服人妻中文乱码| 日韩免费av在线播放| 国产一区二区三区在线臀色熟女| 亚洲avbb在线观看| 一进一出抽搐gif免费好疼| 欧美一区二区亚洲| 久久久精品大字幕| 国产真人三级小视频在线观看| 亚洲精品影视一区二区三区av| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 精品一区二区三区视频在线观看免费| 国产一区二区在线av高清观看| 可以在线观看的亚洲视频| 成人永久免费在线观看视频| 久久久久久久久大av| 精品一区二区三区视频在线 | 亚洲七黄色美女视频| 精品人妻一区二区三区麻豆 | 欧美在线一区亚洲| 日日摸夜夜添夜夜添小说| 在线视频色国产色| 99riav亚洲国产免费| 真人做人爱边吃奶动态| 丰满的人妻完整版| xxxwww97欧美| 日韩中文字幕欧美一区二区| 亚洲av电影在线进入| 欧美极品一区二区三区四区| 午夜两性在线视频| 日本免费a在线| 中文字幕人妻熟人妻熟丝袜美 | 国产97色在线日韩免费| 亚洲中文字幕一区二区三区有码在线看| 久久香蕉国产精品| 蜜桃久久精品国产亚洲av| 悠悠久久av| 免费大片18禁| www.熟女人妻精品国产| 亚洲乱码一区二区免费版| 久久亚洲真实| 久久久久免费精品人妻一区二区| 免费av毛片视频| 午夜激情欧美在线| 精品福利观看| 国产精品自产拍在线观看55亚洲| 亚洲天堂国产精品一区在线| 国产午夜精品久久久久久一区二区三区 | 午夜福利免费观看在线| 亚洲国产精品久久男人天堂| 成人国产一区最新在线观看| 啪啪无遮挡十八禁网站| 国产亚洲精品av在线| 色综合亚洲欧美另类图片| 美女 人体艺术 gogo| 国产精品99久久99久久久不卡| 五月玫瑰六月丁香| 内地一区二区视频在线| 欧美日本亚洲视频在线播放| 久9热在线精品视频| 两性午夜刺激爽爽歪歪视频在线观看| www.色视频.com| 国产中年淑女户外野战色| 国产精品98久久久久久宅男小说| 久久久久久大精品| 亚洲久久久久久中文字幕| 91久久精品电影网| 日韩高清综合在线| 手机成人av网站| 久久香蕉国产精品| 国产综合懂色| 在线观看一区二区三区| 一进一出抽搐动态| 国产一级毛片七仙女欲春2| 69人妻影院| 美女免费视频网站| 午夜精品久久久久久毛片777| 国产精品亚洲美女久久久| 岛国视频午夜一区免费看| 国产69精品久久久久777片| 午夜福利在线观看免费完整高清在 | 一二三四社区在线视频社区8| 免费在线观看成人毛片| 欧美日韩综合久久久久久 | 欧美一区二区精品小视频在线| 我要搜黄色片| 国产精品野战在线观看| 久久中文看片网| 美女高潮的动态| 国内精品美女久久久久久| 日本a在线网址| 成人一区二区视频在线观看| 亚洲,欧美精品.| 蜜桃亚洲精品一区二区三区| 此物有八面人人有两片| 真人一进一出gif抽搐免费| 欧美激情久久久久久爽电影| 免费大片18禁| 亚洲中文日韩欧美视频| 国产国拍精品亚洲av在线观看 | 啦啦啦观看免费观看视频高清| 久久欧美精品欧美久久欧美| 麻豆国产av国片精品| 免费看光身美女| av国产免费在线观看| 久久久久国内视频| 内射极品少妇av片p| 搡老岳熟女国产| 午夜免费激情av| 亚洲va日本ⅴa欧美va伊人久久| 久久久久久久久久黄片| 亚洲精品乱码久久久v下载方式 | 熟女少妇亚洲综合色aaa.| 老熟妇仑乱视频hdxx| 久久天躁狠狠躁夜夜2o2o| 亚洲无线观看免费| 九九热线精品视视频播放| 两个人看的免费小视频| 国产精品久久久人人做人人爽| 国产麻豆成人av免费视频| 又紧又爽又黄一区二区| 丝袜美腿在线中文| 伊人久久精品亚洲午夜| 欧美一区二区亚洲| 午夜日韩欧美国产| 变态另类丝袜制服| 亚洲av成人精品一区久久| avwww免费| 国产亚洲欧美98| 欧美日韩黄片免| 免费搜索国产男女视频| 国产伦一二天堂av在线观看| 精品久久久久久久毛片微露脸| www.熟女人妻精品国产| 少妇人妻一区二区三区视频| 国产中年淑女户外野战色| 精品一区二区三区视频在线观看免费| 久久亚洲真实| 日韩精品中文字幕看吧| 日本a在线网址| 国产伦精品一区二区三区四那| 久久这里只有精品中国| 蜜桃亚洲精品一区二区三区| 亚洲av成人不卡在线观看播放网| 久久午夜亚洲精品久久| 日本黄色视频三级网站网址| 午夜老司机福利剧场| 757午夜福利合集在线观看| 久久久精品欧美日韩精品| 国产色婷婷99| 国产高清视频在线观看网站| 国产 一区 欧美 日韩| avwww免费| 美女高潮的动态| 好男人电影高清在线观看| 99riav亚洲国产免费| 色老头精品视频在线观看| 国产精品精品国产色婷婷| 可以在线观看毛片的网站| 国产精品,欧美在线| 国产成人啪精品午夜网站| 欧美av亚洲av综合av国产av| bbb黄色大片| 国产极品精品免费视频能看的| 日本a在线网址| 国产亚洲精品综合一区在线观看| 少妇的逼好多水| 观看免费一级毛片| 婷婷六月久久综合丁香| 亚洲国产精品久久男人天堂| 五月伊人婷婷丁香| 最新在线观看一区二区三区| 精品人妻偷拍中文字幕| 午夜精品久久久久久毛片777| 白带黄色成豆腐渣| 69av精品久久久久久| 搞女人的毛片| 18禁美女被吸乳视频| 国产亚洲精品综合一区在线观看| 九色成人免费人妻av| 国产aⅴ精品一区二区三区波| 亚洲精品成人久久久久久| xxxwww97欧美| 亚洲自拍偷在线| 国产日本99.免费观看| 亚洲精品亚洲一区二区| 三级毛片av免费| 一进一出抽搐gif免费好疼| 精品欧美国产一区二区三| 久久精品91蜜桃| 一个人观看的视频www高清免费观看| 真人做人爱边吃奶动态| 久久久久久久午夜电影| 一区二区三区国产精品乱码| 欧美一区二区国产精品久久精品| 欧美最黄视频在线播放免费| 亚洲内射少妇av| 日韩大尺度精品在线看网址| 变态另类丝袜制服| 亚洲av中文字字幕乱码综合| 丰满乱子伦码专区| 中文字幕高清在线视频| 国产黄a三级三级三级人| 国产色婷婷99| 床上黄色一级片| 午夜精品一区二区三区免费看| 他把我摸到了高潮在线观看| 制服人妻中文乱码| 欧美黑人欧美精品刺激| 午夜精品一区二区三区免费看| www.色视频.com| 国产亚洲欧美98| 全区人妻精品视频| 欧美成人性av电影在线观看| 91在线观看av| 高清日韩中文字幕在线| 亚洲精品456在线播放app | 一个人观看的视频www高清免费观看| 在线播放国产精品三级| 国产高潮美女av| 国产精品女同一区二区软件| 老师上课跳d突然被开到最大视频| 永久免费av网站大全| 久久精品久久精品一区二区三区| av线在线观看网站| av网站免费在线观看视频 | 51国产日韩欧美| 少妇高潮的动态图| 国产一区二区亚洲精品在线观看| 欧美丝袜亚洲另类| 免费无遮挡裸体视频| 色哟哟·www| 欧美成人午夜免费资源| 亚洲精品日韩在线中文字幕| 亚洲精品影视一区二区三区av| 国产在视频线在精品| 国产av在哪里看| 国产免费一级a男人的天堂| 啦啦啦中文免费视频观看日本| 亚洲图色成人| 91久久精品电影网| 国产精品久久久久久av不卡| 免费看不卡的av| 免费电影在线观看免费观看| 亚洲精品成人av观看孕妇| 国产男女超爽视频在线观看| 大片免费播放器 马上看| .国产精品久久| av国产久精品久网站免费入址| 国产白丝娇喘喷水9色精品| 成人午夜高清在线视频| 国产精品一区二区三区四区免费观看| 搡老妇女老女人老熟妇| 亚洲精华国产精华液的使用体验| 成人国产麻豆网| 中文字幕亚洲精品专区| 国产伦精品一区二区三区视频9| 韩国av在线不卡| 一区二区三区高清视频在线| 亚洲va在线va天堂va国产| 国产黄片美女视频| 午夜亚洲福利在线播放| 久久久欧美国产精品| 国产黄a三级三级三级人| 成人美女网站在线观看视频| 精品久久久久久电影网| 国产精品伦人一区二区| 亚洲av中文av极速乱| 在线 av 中文字幕| 黄色一级大片看看| 九色成人免费人妻av| 国产久久久一区二区三区| 久久久亚洲精品成人影院| 国产亚洲最大av| videos熟女内射| 国产一区二区三区av在线| 毛片女人毛片| 99久久精品一区二区三区| av免费在线看不卡| 99热这里只有精品一区| 自拍偷自拍亚洲精品老妇| 成人av在线播放网站| 亚洲电影在线观看av| 国产老妇伦熟女老妇高清| 亚洲精品乱码久久久久久按摩| 国产美女午夜福利| 国产 一区精品| 中文字幕亚洲精品专区| 最近最新中文字幕免费大全7| 亚洲成人中文字幕在线播放| 极品教师在线视频| 一级毛片 在线播放| 亚洲三级黄色毛片| 成人漫画全彩无遮挡| 日韩国内少妇激情av| 亚洲精品影视一区二区三区av| 亚洲,欧美,日韩| 夫妻午夜视频| 天堂网av新在线| 美女黄网站色视频| 别揉我奶头 嗯啊视频| 日韩精品有码人妻一区| 免费看日本二区| 免费观看a级毛片全部| 国产片特级美女逼逼视频| 建设人人有责人人尽责人人享有的 | 亚洲最大成人中文| 国产综合懂色| 三级经典国产精品| 又粗又硬又长又爽又黄的视频| 一级av片app| 免费看a级黄色片| 欧美xxxx性猛交bbbb| 我的老师免费观看完整版| 最后的刺客免费高清国语| 国产老妇女一区| 精品人妻熟女av久视频| 3wmmmm亚洲av在线观看| 国产日韩欧美在线精品| 欧美性感艳星| 大陆偷拍与自拍| 国产精品美女特级片免费视频播放器| 久久久亚洲精品成人影院| av女优亚洲男人天堂| 免费观看av网站的网址| 一区二区三区乱码不卡18| 亚洲成人av在线免费| 国产黄片视频在线免费观看| av在线亚洲专区| 国产精品爽爽va在线观看网站| 国产乱人偷精品视频| 久久国产乱子免费精品| 色视频www国产| 神马国产精品三级电影在线观看| 99久久精品国产国产毛片| 人体艺术视频欧美日本| 亚洲婷婷狠狠爱综合网| 久久久久久久久久久免费av| 欧美激情在线99| 日韩,欧美,国产一区二区三区| 淫秽高清视频在线观看| 美女被艹到高潮喷水动态| 国产黄a三级三级三级人| 青春草视频在线免费观看| 亚洲18禁久久av| 天堂网av新在线| 五月玫瑰六月丁香| 国产免费又黄又爽又色| 亚洲精品色激情综合| 男女那种视频在线观看| 在线免费观看的www视频| 免费黄网站久久成人精品| 最新中文字幕久久久久| 欧美成人午夜免费资源| 精品久久久久久电影网| 日韩欧美国产在线观看| 日韩大片免费观看网站| 国产麻豆成人av免费视频| 天堂av国产一区二区熟女人妻| 亚洲欧美成人综合另类久久久| 51国产日韩欧美| 国产黄频视频在线观看| 中文字幕亚洲精品专区| 亚洲国产精品专区欧美| 国产免费又黄又爽又色| 久久久久久国产a免费观看| a级一级毛片免费在线观看| av在线天堂中文字幕| 国内精品一区二区在线观看| 丰满少妇做爰视频| 亚洲av日韩在线播放| 国产成人aa在线观看| 亚洲精品,欧美精品| 免费看日本二区| 男人舔女人下体高潮全视频| 亚洲美女搞黄在线观看| 亚洲一级一片aⅴ在线观看| 欧美潮喷喷水| 嫩草影院精品99| 人妻制服诱惑在线中文字幕| 性色avwww在线观看| 在线天堂最新版资源| 欧美日韩国产mv在线观看视频 | 成人无遮挡网站| 午夜福利在线观看免费完整高清在| 高清欧美精品videossex| ponron亚洲| 三级毛片av免费| 黄片无遮挡物在线观看| 高清毛片免费看| 久久久久久久久久黄片| 午夜免费观看性视频| 深爱激情五月婷婷| www.av在线官网国产| 最后的刺客免费高清国语| av在线蜜桃| 国产乱来视频区| 国产 一区精品| 国产久久久一区二区三区| 最近最新中文字幕免费大全7| 久久久a久久爽久久v久久| 有码 亚洲区| 国产亚洲精品av在线| 一个人免费在线观看电影| 看十八女毛片水多多多| 18禁在线无遮挡免费观看视频| 国产成人精品福利久久| 国产精品久久久久久av不卡| 国产av码专区亚洲av| a级毛片免费高清观看在线播放| 中文字幕久久专区| 精品久久国产蜜桃| 免费av观看视频| 日日干狠狠操夜夜爽| 三级经典国产精品| 老司机影院成人| 国产黄色免费在线视频| 如何舔出高潮| 精品一区二区免费观看| 欧美97在线视频| 激情 狠狠 欧美| 久久久久国产网址| 一区二区三区乱码不卡18| 男的添女的下面高潮视频| 韩国高清视频一区二区三区| 在线观看一区二区三区| 99热这里只有是精品50| 国产精品爽爽va在线观看网站| 干丝袜人妻中文字幕| 亚洲美女视频黄频| 欧美97在线视频| 美女xxoo啪啪120秒动态图| h日本视频在线播放| 男人爽女人下面视频在线观看| 日韩精品青青久久久久久| 欧美 日韩 精品 国产| 久久久久久久久大av| 午夜精品一区二区三区免费看| 亚洲av免费在线观看| 亚洲无线观看免费| 亚洲国产色片| 免费av毛片视频| 内射极品少妇av片p| 亚洲成人久久爱视频| 国产精品一及| 免费大片黄手机在线观看| 午夜激情欧美在线| 黑人高潮一二区| 免费看光身美女| 禁无遮挡网站| 只有这里有精品99| 直男gayav资源| 国产精品99久久久久久久久| 日日干狠狠操夜夜爽| 成人亚洲欧美一区二区av| 久久久久精品性色| 国产三级在线视频| 亚洲精品乱码久久久v下载方式| 久久久久精品久久久久真实原创|