- 1School of Forensic Medicine, Southern Medical University, Guangzhou, China
- 2Institute of Archaeological Science, Fudan University, Shanghai, China
- 3School of Basic Medicine and Life Science, Hainan Medical University, Haikou, China
- 4School of Basic Medicine, Gannan Medical University, Ganzhou, China
Maoming is located in the southwest region of Guangdong Province and is the cradle of Gaoliang culture, which is the representative branch of Lingnan cultures. Historical records showed that the amalgamations between Gaoliang aborigines and distinct ethnic minorities had some influences on the shaping of Gaoliang culture, especially for the local Tai-kadai language-speaking Baiyue and Han Chinese from Central China. However, there is still no exact genetic evidence for the influences on the genetic pool of Maoming Han, and the genetic relationships between Maoming Han and other Chinese populations are still unclear. Hence, in order to get a better understanding of the paternal genetic structures and characterize the forensic features of 27 Y-chromosomal short tandem repeats (Y-STRs) in Han Chinese from Guangdong Maoming, we firstly applied the AmpFLSTR® Yfiler® Plus PCR Amplification Kit (Thermo Fisher Scientific, Waltham, MA, United States) to genotype the haplotypes in 431 Han males residing in Maoming. A total of 263 different alleles were determined across all 27 Y-STRs with the corresponding allelic frequencies from 0.0004 to 0.7401, and the range of genetic diversity (GD) was 0.4027 (DYS391) to 0.9596 (DYS385a/b). In the first batch of 27 Yfiler data in Maoming Han, 417 distinct haplotypes were discovered, and nine off-ladder alleles were identified at six Y-STRs; in addition, no copy number variant or null allele was detected. The overall haplotype diversity (HD) and discrimination capacity (DC) of 27 Yfiler were 0.9997 and 0.9675, respectively, which demonstrated that the 6-dye and 27-plex system has sufficient system effectiveness for forensic applications in Maoming Han. What is more, the phylogenetic analyses indicated that Maoming Han, which is a Southern Han Chinese population, has a close relationship with Meizhou Kejia, which uncovered that the role of the gene flows from surrounding Han populations in shaping the genetic pool of Maoming Han cannot be ignored. From the perspectives of genetics, linguistics, and geographies, the genetic structures of Han populations correspond to the patterns of the geographical-scale spatial distributions and the relationships of language families. Nevertheless, no exact genetic evidence supports the intimate relationships between Maoming Han and Tai-Kadai language-speaking populations and Han populations of Central Plains in the present study.
Introduction
Maoming, a city located in the southwest of Guangdong Province (Figure 1), is the cradle of Gaoliang culture (Zhou, 2019). Gaoliang culture, one of the representative Lingnan cultures, could be dated back to the Han Dynasty (111 B.C.) in Chinese history (He, 2012). The aborigines living in Gaoliang mountainous areas and the basins between Jian River and Moyang River are the inheritors of Gaoliang culture, which are represented by the customs of Nianli (a special celebration for New Year) and Piaose (a form of dramatic plastic arts on moving stages) (Chen, 2013). Since the Southern and Northern Dynasties (420–589 A.D.), the intermarriages accelerated national amalgamations between Gaoliang aborigines and other ethnic minorities in ancient Gaoliang District (He, 2012; Chen, 2013)12). Therefore, Gaoliang culture was influenced by the convergences between Gaoliang aborigines and different ethnic groups (Gao, 2007). Moreover, some archeological records also hinted that the population structures of Gaoliang aborigines might be affected by the local Baiyue (a Tai-kadai language-speaking population in ancient China) and Han Chinese from Central China with the increasingly social activities of mixed marriages, population migrations, and trade contacts in the long course of history (Gao, 2007). Maoming Hans, the descendants of Gaoliang aborigines, speak Cantonese (Gaoyang Pian), which is one branch of Sino-Tibetan language family (Ding, 2010). From the perspective of languages, the language of Maoming Han (Cantonese) did not seem to be impacted by Tai-Kadai groups (Baiyue). Hence, there is still no exact genetic evidence for the influences on the genetic pool of Maoming Han, and the genetic relationships between Maoming Han and other surrounding populations are still unclear.
Figure 1. Geographical locations of population distributions and sampling information. (A) The geographical distribution of Maoming Han and other Chinese populations (Han Chinese and ethnic minorities), which was analyzed in the present study. (B) The geographical distribution of Maoming Han and other five Han Chinese populations in South China (our studied group is marked as a red pentagram).
Y-chromosomal short tandem repeats (Y-STRs) have been regarded as a valuable tool in forensic genetics (Kayser, 2017), genealogy (Kayser et al., 2007), human evolution (Underhill and Kivisild, 2007), archeology (Calafell and Larmuseau, 2017), population history (Kayser et al., 2000; Jobling and Tyler-Smith, 2017), and male medical genetics (Hughes and Page, 2016; Jobling and Tyler-Smith, 2017). In addition, Y chromosomal variant analysis for determining the patterns of present and past flows of genes between populations is very helpful (Oppenheimer, 2012). The use of Y-STRs also allows the simultaneous analysis of closely related and distantly related populations (Ballantyne and Kayser, 2012). The 6-dye and 27-plex AmpFLSTR® Yfiler® Plus PCR Amplification Kit (Thermo Fisher Scientific, Waltham, MA, United States) includes 17 Yfiler loci (DYS19, DYS385a/b, DYS389I/II, DYS390, DYS391, DYS392, DYS393, DYS437, DYS438, DYS439, DYS448, DYS456, DYS458, DYS635, and Y GATA H4) plus three highly polymorphic Y-STR loci (DYS460, DYS481, and DYS533) and seven rapidly mutating Y-STR loci (DYF387S1a/b, DYS449, DYS518, DYS570, DYS576, and DYS627) in an effort to improve discrimination of related individuals (Gopinath et al., 2016). Y-chromosome STR haplotype reference database (YHRD)1 is an internet-accessible worldwide reference database of Y chromosome profiles, which contributed to provide a worldwide and high-quality Y-STR haplotype data from distinct human populations for forensic purposes and population genetics (Willuweit and Roewer, 2015). As yet, little is known about the genetic backgrounds of the aforementioned 27 Y-STRs in Maoming Han, and the forensic-related Y chromosome variation data in Guangdong Maoming still remains blank in YHRD.
Hence, in order to get a better understanding of the paternal genetic structure and characterize the forensic resolution of 27 Y-STRs in Han Chinese from Guangdong Maoming, we used the 6-dye and 27-plex Y-STR system to genotype the haplotypes in 431 Han males residing in Maoming city. Furthermore, we explored the genetic relationships between Maoming Han and Chinese populations of Southern and Northern China from the perspectives of geographies, linguistics, and genetics.
Materials and Methods
Sample Preparation
In this study, a total of 431 unrelated Han Chinese males were recruited from Maoming city, Guangdong Province, China (Figure 1). The inclusion criteria were as follows: (1) healthy individuals without any underlying diseases (including but not limited to cardiovascular diseases, metabolic diseases, chronic wasting diseases, immunologic diseases, etc.); (2) unrelated males and any two individuals who have no blood relationship for up to three generations; (3) the volunteers’ parents and grandparents are aboriginals and have non-consanguineous marriages of the same ethnic group for at least three generations, which was confirmed by the volunteers’ self-declared statements; and (4) the language Cantonese is the mother tongue of Maoming volunteers, and any self-declared Maoming Han who could not speak Cantonese would be excluded from our cohort. Blood samples of all Maoming volunteers were collected using FTA cards (WhatmanTM, GE Healthcare, Chicago, IL, United States) with written informed consents from participants. All the experimental procedures were performed following the standards of the Declaration of Helsinki. This study was approved by the Medical Ethics Committee of Hainan Medical University (no. HYLL-2020-012).
DNA Extraction, Amplification, and Genotyping
Genomic DNA was extracted using the TIANamp Blood Spot DNA Kit (TIANGEN BIOTECH, Beijing, China) according to the manufacturer’s protocol. The quantity of the DNA templates was determined using QubitTM dsDNA HS Assay Kit (Thermo Fisher Scientific, Waltham, MA, United States) on the Qubit 4.0 Fluorometer (Thermo Fisher Scientific, Waltham, MA, United States) according to the manufacturer’s instructions. Based on the quantitative results, DNA samples were normalized to 2.0 ng/μl and stored at −20°C until amplification.
The amplification of the 6-dye multiplex PCR-CE-based AmpFLSTR® Yfiler® Plus PCR Amplification Kit (Thermo Fisher Scientific, Waltham, MA, United States) was performed in a single multiplex PCR reaction (25 μl in total, containing 10 μl master mix, 5 μl primer mix, and 10 μl genomic DNA) on a Veriti® 96-Well Thermal Cycler System (Thermo Fisher Scientific, Waltham, MA, United States) following the manufacturer’s instructions. Amplified products were separated by capillary electrophoresis (CE) on a 3500xL Genetic Analyzer (Thermo Fisher Scientific, Waltham, MA, United States). The separation of CE-based amplified products was conducted according to our previous studies (Fan et al., 2019a; Liu et al., 2020).
Statistic and Population Genetic Analyses
Allele and haplotype frequencies as well as forensic parameters were calculated using direct counting. The relevant forensic parameters contained genetic diversity (GD), haplotype diversity (HD), discrimination capacity (DC), and random match probability (RMP). GD was calculated according to the following formula:
where n is the total sample size, and p_i indicates the frequency of i-th allele. HD was computed in the same formula as GD, except that p_i refers to the frequency of i-th haplotype. DC is equal to the ratio of different haplotypes to the total sample size. Computed with the formula , RMP is the probability that a particular DNA profile would appear in a population and that a “match” would occur by coincidence. In forensic statistics, a lower RMP value indicates higher strength of evidence provided by genetic analysis.
Population pairwise genetic distance (Rst) is commonly used for estimating the population differences and computing the genetic relationships among different populations (Fan et al., 2019b; Li et al., 2020). By using the “AMOVA&MDS tool” on YHRD, pairwise Rst and corresponding p values based on 17 Yfiler between Maoming Han and reference populations were estimated by analysis of molecular variance (AMOVA) and visualized in multidimensional scaling (MDS) plot, which were used show the reduced dimensionality spatial representation of the populations. Additionally, phylogenetic relationships among Han Chinese populations from Southern and Northern mainland China as well as those between 6 Han Chinese and 16 ethnic minorities were depicted in the Molecular Evolutionary Genetics Analysis-X (MEGA-X) software (Kumar et al., 2018) by a neighbor-joining (N-J) phylogenetic tree (Saitou and Nei, 1987) based upon the Rst genetic distance matrix, respectively.
Quality Control
The recommendations of the DNA Commission of the Chinese National Standards, the Scientific Working Group on DNA Analysis Methods (SWGDAM) (SWGDAM, 2010), and the DNA Commission of the International Society of Forensic Genetics (ISFG) (Gusmao et al., 2006; Carracedo et al., 2013; Roewer et al., 2020) for analysis of Y-STRs were strictly followed. Control DNA 007 was employed as a positive control, while ddH2O was used as a negative control for each batch of amplification and genotyping. Additionally, the laboratory has passed the proficiency testing for Y-STR typing organized by YHRD and has been accredited in accordance with ISO/IEC 17025:2005 and the China National Accreditation Service for Conformity Assessment (CNAS). The haplotype data of 431 unrelated male individuals from Guangdong Maoming Han population in the present study have been submitted to YHRD database and received the accession number YA004720 (Maoming Han, n = 431). The Y-STR profiles with off-ladders were re-amplified and re-genotyped by Goldeneye DNATM ID 27YB system (Goldeneye® Technology Ltd., Beijing, China).
Results and Discussion
In the present study, a total of 431 unrelated male individuals from Han Chinese in Guangdong Maoming were genotyped including 27 Y-STR loci using the AmpFLSTR® Yfiler® Plus PCR Amplification Kit (Thermo Fisher Scientific, Waltham, MA, United States). In order to evaluate the forensic features of Maoming Han population, we set up two datasets, Yfiler set and Yfiler Plus set, including 17 and 27 Y-STRs, respectively. In addition, a series of comprehensive population genetic analyses were conducted between Maoming Han and other southern and northern Chinese populations. In short, the aims of this study were to feature the forensic characteristics of 27 Y-STRs in Maoming Han, clarify the paternal genetic structures of Maoming Han, and get a better understanding of the genetic relationships between Maoming Han and other Chinese populations from the perspectives of geographics, linguistics, and genetics.
Forensic Characteristics
Forensic Features of Yfiler Set (17 Y-STRs)
As illustrated in Supplementary Table 1, a total of 147 distinct alleles were identified across all 17 Y-STRs in Maoming Han with the corresponding allelic frequencies from 0.0023 to 0.7401 (DYS391). Overall, 17 Yfiler loci were relatively highly polymorphic in Maoming Han. The range of allele numbers was 4 (DYS391, DYS437, and DYS438) ∼ 55 (DYS385a/b), and the lowest and highest estimates of GD corresponded to loci DYS391 (0.4027) and DYS385a/b (0.9596). Except for DYS391 (0.4027) and DYS438 (0.4049), the GD values for other 17 Yfiler loci were greater than 0.5. The haplotypes and haplotype frequencies of Yfiler in Maoming Han are shown in Supplementary Table 2. There were 371 different haplotypes observed in 431 Maoming Han individuals, of which 328 (88.41%) were unique, 33 occurred twice (S011–S043), 6 (S005–S010) were observed thrice, 1 (H004) was shared by 4 individuals, and 3 (S001–S003) were shared by 5 individuals. We observed four confirmed microvariants [18.2 (twice) and 19.2 at DYS448 and 18.2 at DYS458]. The overall HD was 0.9994 with a DC of 0.8608.
Forensic Features of Yfiler Plus Set (27 Y-STRs)
Allele frequency distributions and haplotype frequencies of Yfiler Plus for Maoming Han are presented in Supplementary Tables 3, 4. A total of 263 different alleles were observed, and the number of distinct alleles ranged from 4 for DYS391, DYS437, and DYS438 to 55 for DYS385a/b. Allele frequencies varied from 0.0004 to 0.7401. All 10 newly added loci got GD values higher than 0.5, especially for the added multi-copy DYF387S1a/b (0.9682). DYS385a/b (0.9596) on the one hand while DYS391 (0.4027) and DYS438 (0.4049) on the other marked the extremes of the GD distribution (with GD values less than 0.5). Genotyping with the 27 Y-STRs determined 417 distinct haplotypes in the population of Maoming Han, of which 405 (97.12%) were unique, 10 different haplotypes were identified twice (H003–H012), and 2 (H001-H002) appeared thrice. In addition to 18.2 and 19.2 at DYS448 and 18.2 at DYS458, intermediate alleles were also observed at the DYS449 (34.2), DYS518 (37.2), DYF387S1 (37.2), and DYS627 (17.2 and 18.2) loci. The overall HD and DC were calculated to be 0.9997 and 0.9675, respectively.
In this study, duplicated or triplicated alleles and null alleles were not detected in both Yfiler set and Yfiler Plus set. The analysis of genotype data revealed that DYS385a/b and DYF387S1a/b showed higher GD in Maoming Han, which were the same as other Chinese populations (Fan et al., 2018a). Forensic parameters based on different sets of Y-STR loci were calculated and listed in Table 1, indicating that as the number of Y-STR loci increased, more distinct haplotypes were identified, and HD and DC were also increased in the present study.
Genetic Differences Between Maoming Han and Han Chinese Populations From Southern and Northern Mainland China
From the intercontinental perspective, a MDS was performed between Maoming Han and 21 worldwide populations (Kim et al., 2001; Miranda et al., 2001; Roewer et al., 2005; Mizuno et al., 2008; Alam et al., 2010; Laouina et al., 2011; Wolfgramm Ede et al., 2011; Ramos-Luis et al., 2014; Xu et al., 2015; Han et al., 2016; Rapone et al., 2016; Wang et al., 2017; Alonso Morales et al., 2018; Fan et al., 2018c; Singh et al., 2018; Henry et al., 2019; Lang et al., 2019; Li et al., 2019; Salvador et al., 2019; Reid and Heathfield, 2020). As shown in Supplementary Figure 1, Maoming Han clustered with Han Chinese populations, while other populations got together in accordance with geographical patterns relatively. In addition, to further explore the genetic affinity among Han Chinese populations in mainland China, the degree of differentiation between Maoming Han and Han populations from different administrative divisions of China was assessed by AMOVA and visualized in an MDS plot. Pairwise Rst and corresponding p values based on the 17 Yfiler among Maoming Han and 11 other Han Chinese populations from North to South across the mainland China (Yang et al., 2014; Shu et al., 2015; Han et al., 2016; Wang et al., 2016, 2017; Xu et al., 2016; Yao et al., 2016; Jiang et al., 2017; Chen et al., 2018; Fan et al., 2018c; Lang et al., 2019; Yin et al., 2020) are listed in Supplementary Table 5. Significant genetic differences were observed between Maoming Han and all other Han populations (p = 0.05) except for Guangdong Han. However, after Bonferroni’s correction (p values = 0.05/431 ≈ 0.0001, n = 431), there were differences between Maoming Han and Northern Han Chinese from Gansu (Rst = 0.0156, p < 0.0001) and Shandong (Rst = 0.0123, p < 0.0001), while other Han Chinese populations had no significances with Maoming Han, especially for Guangdong Han (Rst = 0.0008, p = 0.1509), indicating the genetic affinity between Southern Han and Northern Han in mainland China. The MDS plot (Figure 2A) based on Rst values clearly demonstrated that the Northern Han populations, Gansu, Jilin, Hebei, Shandong, and Henan (Nanyang), were grouped in the bottom left side, while Hainan, Hubei (Wuhan), Hunan, Guizhou (Zunyi), and Guangdong isolated from the northern cluster and gathered into a Southern Han cluster in the upper right side with Maoming Han.
Figure 2. Genetic similarities and differences between Maoming Han and other 11 Han Chinese populations from Northern and Southern China based on Rst values. (A) Multidimensional scaling (MDS) plot between Maoming Han and other 11 Han Chinese populations; (B) neighbor-joining (N-J) phylogenetic tree among 12 Han Chinese populations from Southern and Northern mainland China.
Subsequently, to make further confirmation about the genetic relationships between Maoming Han and other Han Chinese populations, an N-J phylogenetic tree based on Rst values was constructed (Figure 2B). We found that two main branches could be clearly identified in the N-J phylogenetic tree. The upper branch was Southern Han cluster, which was composed of Maoming, Guangdong, Guangxi, Hainan, Guizhou (Zunyi), Hunan, and Hubei (Wuhan), while Henan (Nanyang), Hebei, Shandong, Gansu, and Jilin got together in the bottom clade as Northern Han cluster. From the perspective of genetics, the analyses above indicated that Maoming Han is a Southern Han population and has relatively close relationships with Guangdong Han, followed by Guangxi Han (Rst = 0.0031, p = 0.0155) and Hainan Han (Rst = 0.0033, p = 0.0183). The phylogenetic structures of Han Chinese populations from Southern and Northern mainland China in the phylogenetic dendrogram were in line with the results of the MDS. From the geographical scale, Guangdong, Guangxi, Hainan, Guizhou, and Hunan belong to southern administrative divisions of mainland China, while Jilin, Gansu, Shandong, Hebei, and Henan are subordinated to the northern administrative divisions of mainland China (Figure 1). Maoming, lying in the southwest of Guangdong Province, has close geographical distances with Guangxi Zhuang Autonomous Region and Hainan Province. Even though different Han Chinese populations from distinct administrative divisions of mainland China have genetic and linguistic homogeneousness, the genetic distances and population structures of Han Chinese are in accordance with the geographical-scale pattern to a certain extent in mainland China.
Genetic Affinities and Differentiations Among Maoming Han, Other Han Populations, and Ethnic Minorities From China
According to the history records, the population structures of Maoming Han were mainly affected by the intermarriages with local Tai-kadai language-speaking Baiyue population and the south migrations of Han Chinese from Central China (He, 2012; Chen, 2013)12), while the above population analyses between our studied population and other Southern and Northern Han populations did not hint the relatively intimate relationships between Maoming Han and Han populations of Central Plains. To reveal the genetic structures among Maoming Han, surrounding Han populations and other 16 Chinese ethnic groups (Zhu et al., 2005; Shi et al., 2011; Shan et al., 2014; Zeng et al., 2014; Gao et al., 2015; Guo et al., 2015; Ou et al., 2015; Shu et al., 2015; Bian et al., 2016; Fu et al., 2016; Hu et al., 2017; Wang et al., 2017, 2019; Zhang et al., 2017; Zhao et al., 2017; Chen et al., 2018; Fan et al., 2018a,b,c; Du et al., 2019; Lang et al., 2019; Song et al., 2019; Xie et al., 2019; Ding et al., 2020; Feng et al., 2020; Guan et al., 2020), pairwise Rst and corresponding p values were calculated based on 17 Yfiler. As presented in Supplementary Table 6, no difference was observed between Maoming Han and Meizhou Kejia (Rst 0.0007, p = 0.2899), while significant genetic differences were observed between Maoming Han and all other Han Chinese or ethnic groups (p < 0.05). However, after Bonferroni’s correction (p values ≈ 0.0001), there were no differences between Maoming Han and surrounding Han Chinese populations (Supplementary Table 7). Furthermore, we found that Inner Mongolia Daur (Rst = 0.1614), Xinjiang Uighur (Rst = 0.1407), and Lhasa U-Tsang Tibetan (Rst = 0.1294) in Northern China had the longest genetic distances with Maoming Han, while the closest genetic distance was seen in Meizhou Kejia (Rst = 0.0007), followed by Guangxi Han (Rst = 0.0031), Hainan Han (Rst = 0.0033), then by Jieyang Han (Rst = 0.0039), and Chaoshan Han (Rst = 0.0083).
On the basis of Rst values of 22 Chinese populations, an MDS plot (Figure 3) and an N-J phylogenetic tree (Figure 4) were performed to depict the forensic genetic landscape of Chinese Han and ethnic groups. As shown in Figure 3A, the Han Chinese populations were closely related to each other and therefore formed a Southern Han cluster, while other 16 minorities were relatively isolated from the Southern Han cluster and dispersed into four main clusters, which were in accord with the distributions of language families in some degree. The Hmong-Mien language-speaking groups, Miao, Yao, and She, clustered together at the upper left, and Dai, Lingao, and Li gathered in the bottom left as the Tai-Kadai language-speaking cluster, while the Tibeto-Burman-language speaking and Altaic-language speaking groups located together at the bottom right with relative separated positions. In addition, Figure 3B indicates the genetic relationships between Maoming Han and surrounding Han populations, which indicated that Maoming Han had a close relationship with Meizhou Kejia. Kejia, also known as Hakka, is a branch of Han Chinese that has a wide distribution in Guangdong Province. The genetic pool of Maoming Han was influenced by the surrounding Han populations, while no direct genetic evidence verified that the Tai-Kadai language-speaking populations contributed to the Maoming Han genetic pool. Furthermore, Meizhou Kejia was first clustered with the Maoming Han, followed by Guangxi Han and Hainan Han, then by Jieyang and Chaoshan in the phylogenetic tree (Figure 4). The tree also revealed that different populations were gathered into two cluster according to their geographical distributions and separated into two main branches: one represented the Altaic language-speaking populations; the other one stood for the Sino-Tibetan language-speaking populations (Han Chinese, Tibeto-Burman, Hmong-Mien, and Tai-Kadai), which was roughly congruent with the results of corresponding MDS (Figure 3).
Figure 3. Multidimensional scaling plots among Maoming Han, surrounding Han populations, and Chinese ethnic minorities based on pairwise Rstvalues. (A) MDS plot among Maoming Han, five other Chinese Han populations, as well as 16 ethnic minorities from China; (B) MDS plot between Maoming Han and surrounding Southern Han Chinese populations.
Figure 4. Phylogenetic analysis among the Maoming people, five other Chinese Han populations, as well as 16 reference minorities from China based on pairwise Rst values (the language families and geographic distributions of different populations are displayed in the inner circle and the outer layer, respectively).
From the perspective of linguistics, geographies, and genetics, the phylogenetic analyses (both the MDS plots and N-J phylogenetic tree) demonstrated that Maoming Han was isolated from Chinese ethnic minority groups relatively and had a relatively close genetic relationships with Southern Han populations, especially for those with the same dialects and intimate geographical distances (Meizhou Kejia, Guangxi Han, and Hainan Han), which indicated that there might be gene flows between Maoming Han and the surrounding Han populations. In addition, the genetic structures of Han populations correspond to the patterns of the geographical-scale spatial distributions and the relationships of language families. In total, the results of above population genetic analyses indicated that Maoming Han, which is a Southern Han Chinese population, has a relatively close relationship with Meizhou Kejia; therefore, the role of the gene flows from surrounding Han populations in shaping the genetic pool of Maoming Han cannot be ignored.
Conclusion
In the present study, a total of 431 unrelated Guangdong Maoming Han were investigated using the AmpFLSTR® Yfiler® Plus PCR Amplification Kit (Thermo Fisher Scientific, Waltham, MA, United States). The high-quality 27 Y-STR haplotype data of Maoming Han were obtained and submitted to YHRD with the accession number YA004720. Overall, 263 different alleles were identified across all 27 Y-STRs with the number of distinct alleles from 4 to 55. Allele frequencies varied from 0.0004 to 0.7401, and the lowest and highest estimates of GD corresponded to loci DYS391 (0.4027) and DYS385a/b (0.9596), respectively. Genotyping with the 27 Y-STRs determined 417 distinct haplotypes in the population of Maoming Han, of which 405 (97.12%) were unique. In the first batch of 27 Yfiler data for Maoming Han, nine intermediate alleles were detected at six Y-STR loci; in addition, duplicated or triplicated alleles and null alleles were not observed. Based on the comparisons of forensic parameters for different sets of Y-STRs (17 Yfiler set and 27 Yfiler Plus set), it demonstrated that the improvements of HD and DC are accompanied by the increasing numbers of Y-STRs. The overall HD and DC of 27 Yfiler in Maoming Han were calculated to be 0.9997 and 0.9675, respectively.
From the perspectives of genetics, linguistics, and geographies, different Han Chinese populations from distinct administrative divisions of mainland China have genetic and linguistic homogeneousness, and the genetic distances and population structures of Han Chinese are in accordance with the geographical-scale pattern to a certain extent in mainland China. Maoming Han, a Southern Han population, has a relatively close genetic relationship with Meizhou Kejia, which has the same language family and has intimate geographical distances with Maoming Han, while no exact genetic evidence supports that there are intimate relationships between Maoming Han and Tai-Kadai language-speaking populations and Han populations of Central Plains. At the same time, we found that the genetic structures of Han populations correspond to the patterns of the geographical-scale spatial distributions and the relationships of language families. As a whole, the sufficient systematic efficiencies of AmpFLSTR® Yfiler® Plus PCR Amplification Kit in Maoming Han demonstrated that it can be widely applied in the population of Guangdong Maoming Han for forensic purposes, and Maoming Han, which is a Southern Han Chinese population, and has a relatively close relationship with Meizhou Kejia; therefore, the role of the gene flows from surrounding Han populations in shaping the genetic pool of Maoming Han cannot be ignored.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics Statement
The studies involving human participants were reviewed and approved by Medical Ethics Committee of the Hainan Medical University. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
HF made significant contributions in the conceptualization, resources, software, formal analysis, and project administration. YL and QX made significant contributions in the investigation. LW performed the validation. QX performed the data curation. HF and QX performed the visualization, wrote and prepared the original draft, and reviewed and edited the manuscript. PQ and S-QW made significant contributions in the supervision of the study. PQ acquired funding for the study. All authors reviewed the manuscript.
Funding
This study was supported by grants from the National Natural Science Foundation of China (NSFC, No. 81971786).
Conflict of Interest
The 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.
The handling editor declared a past co-authorship with the authors LW and S-QW.
Acknowledgments
We would like to thank the donors who contributed samples for this study.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2021.690504/full#supplementary-material
Footnotes
References
Alam, S., Ali, M. E., Ferdous, A., Hossain, T., Hasan, M. M., and Akhteruzzaman, S. (2010). Haplotype diversity of 17 Y-chromosomal STR loci in the Bangladeshi population. Forensic Sci. Int. Genet. 4, e59–e60. doi: 10.1016/j.fsigen.2009.05.005
Alonso Morales, L. A., Casas-Vargas, A., Rojas Castro, M., Resque, R., Ribeiro-Dos-Santos, ÂK., Santos, S., et al. (2018). Paternal portrait of populations of the middle Magdalena river region (Tolima and Huila, Colombia): new insights on the peopling of central America and northernmost South America. PLoS One 13:e0207130. doi: 10.1371/journal.pone.0207130
Ballantyne, K. N., and Kayser, M. (2012). Additional Y-STRs in forensics: why, which, and when. Forensic Sci. Rev. 24, 63–78.
Bian, Y., Zhang, S., Zhou, W., Zhao, Q., Siqintuya, Zhu, R., et al. (2016). Analysis of genetic admixture in uyghur using the 26 Y-STR loci system. Sci. Rep. 6:19998. doi: 10.1038/srep19998
Calafell, F., and Larmuseau, M. H. D. (2017). The Y chromosome as the most popular marker in genetic genealogy benefits interdisciplinary research. Hum. Genet. 136, 559–573. doi: 10.1007/s00439-016-1740-0
Carracedo, A., Butler, J. M., Gusmao, L., Linacre, A., Parson, W., Roewer, L., et al. (2013). New guidelines for the publication of genetic population data. Forensic Sci. Int. Genet. 7, 217–220. doi: 10.1016/j.fsigen.2013.01.001
Chen, P., He, G., Zou, X., Zhang, X., Li, J., Wang, Z., et al. (2018). Genetic diversities and phylogenetic analyses of three Chinese main ethnic groups in southwest China: a Y-chromosomal STR study. Sci. Rep. 8:15339. doi: 10.1038/s41598-018-33751-x
Ding, J., Fan, H., Zhou, Y., Wang, Z., Wang, X., Song, X., et al. (2020). Genetic polymorphisms and phylogenetic analyses of the Ü-Tsang Tibetan from Lhasa based on 30 slowly and moderately mutated Y-STR loci. Forensic Sci. Res. doi: 10.1080/20961790.2020.1810882 [Epub ahead of print].
Du, W., Wu, W., Wu, Z., Guo, L., Wang, B., and Chen, L. (2019). Genetic polymorphisms of 32 Y-STR loci in Meizhou Hakka population. Int. J. Legal Med. 133, 465–466. doi: 10.1007/s00414-018-1845-1
Fan, H., Wang, X., Chen, H., Li, W., Wang, W., and Deng, J. (2019a). The Ong Be language-speaking population in Hainan island: genetic diversity, phylogenetic characteristics and reflections on ethnicity. Mol. Biol. Rep. 46, 4095–4103. doi: 10.1007/s11033-019-04859-8
Fan, H., Wang, X., Chen, H., Long, R., Liang, A., Li, W., et al. (2018a). The evaluation of forensic characteristics and the phylogenetic analysis of the Ong Be language-speaking population based on Y-STR. Forensic Sci. Int. Genet. 37, e6–e11. doi: 10.1016/j.fsigen.2018.09.008
Fan, H., Wang, X., Chen, H., Zhang, X., Huang, P., Long, R., et al. (2018b). Population analysis of 27 Y-chromosomal STRs in the Li ethnic minority from Hainan province, southernmost China. Forensic Sci. Int. Genet. 34, e20–e22. doi: 10.1016/j.fsigen.2018.01.007
Fan, H., Wang, X., Ren, Z., He, G., Long, R., Liang, A., et al. (2019b). Population data of 19 autosomal STR loci in the Li population from Hainan province in southernmost China. Int. J. Legal Med. 133, 429–431. doi: 10.1007/s00414-018-1828-2
Fan, H., Zhang, X., Wang, X., Ren, Z., Li, W., Long, R., et al. (2018c). Genetic analysis of 27 Y-STR loci in Han population from Hainan province, southernmost China. Forensic Sci. Int. Genet. 33, e9–e10. doi: 10.1016/j.fsigen.2017.12.009
Feng, R., Zhao, Y., Chen, S., Li, Q., Fu, Y., Zhao, L., et al. (2020). Genetic analysis of 50 Y-STR loci in Dong, Miao, Tujia, and Yao populations from Hunan. Int. J. Legal Med. 134, 981–983. doi: 10.1007/s00414-019-02115-z
Fu, X., Fu, Y., Liu, Y., Guo, J., Liu, Y., Guo, Y., et al. (2016). Genetic polymorphisms of 26 Y-STR loci in the Mongolian minority from Horqin district. China. Int. J. Legal Med. 130, 941–946. doi: 10.1007/s00414-016-1387-3
Gao, H. (2007). A Brief Analysis of the Integration of the Han and Li Nationalities in Ancient Gaoliang. Journal of Guangzhou Institute of Socialism(04). 38–41.
Gao, T., Yun, L., Gu, Y., He, W., Wang, Z., and Hou, Y. (2015). Phylogenetic analysis and forensic characteristics of 12 populations using 23 Y-STR loci. Forensic Sci. Int. Genet. 19, 130–133. doi: 10.1016/j.fsigen.2015.07.006
Gopinath, S., Zhong, C., Nguyen, V., Ge, J., Lagace, R. E., Short, M. L., et al. (2016). Developmental validation of the Yfiler((R)) plus PCR amplification kit: an enhanced Y-STR multiplex for casework and database applications. Forensic Sci. Int. Genet. 24, 164–175. doi: 10.1016/j.fsigen.2016.07.006
Guan, T., Song, X., Xiao, C., Sun, H., Yang, X., Liu, C., et al. (2020). Analysis of 23 Y-STR loci in Chinese Jieyang Han population. Int. J. Legal Med. 134, 505–507. doi: 10.1007/s00414-019-02019-y
Guo, F., Zhang, L., and Jiang, X. (2015). Population genetics of 17 Y-STR loci in Xibe ethnic minority from Liaoning province, Northeast China. Forensic Sci. Int. Genet. 16, 86–87. doi: 10.1016/j.fsigen.2014.12.007
Gusmao, L., Butler, J. M., Carracedo, A., Gill, P., Kayser, M., Mayr, W. R., et al. (2006). DNA commission of the international society of forensic genetics (ISFG): an update of the recommendations on the use of Y-STRs in forensic analysis. Int. J. Legal Med. 120, 191–200. doi: 10.1007/s00414-005-0026-1
Han, Y., Li, L., Liu, X., Chen, W., Yang, S., Wei, L., et al. (2016). Genetic analysis of 17 Y-STR loci in Han and Korean populations from Jilin province, Northeast China. Forensic Sci. Int. Genet. 22, 8–10. doi: 10.1016/j.fsigen.2016.01.003
He, H. (2012). Expanding the Research Field and Creating a New Cultural Field: A Review of Feng Guixiong and Wu Gang’s Madam Xian’s Family and Yangjiang in Sui and Tang dynasties. Nan Fang Lun Kan(03). 91–93.
Henry, J., Dao, H., Scandrett, L., and Taylor, D. (2019). Population genetic analysis of Yfiler(®) Plus haplotype data for three South Australian populations. Forensic Sci. Int. Genet. 41, e23–e25. doi: 10.1016/j.fsigen.2019.03.021
Hu, L., Gu, T., Fan, X., Yuan, X., Rao, M., Pang, J. B., et al. (2017). Genetic polymorphisms of 24 Y-STR loci in Hani ethnic minority from Yunnan province, Southwest China. Int. J. Legal Med. 131, 1235–1237. doi: 10.1007/s00414-017-1543-4
Hughes, J. F., and Page, D. C. (2016). The history of the Y chromosome in man. Nat. Genet. 48, 588–589. doi: 10.1038/ng.3580
Jiang, W., Gong, Z., Rong, H., Guan, H., Zhang, T., Zhao, Y., et al. (2017). Population genetics of 26 Y-STR loci for the Han ethnic in Hunan province, China. Int. J. Legal Med. 131, 115–117. doi: 10.1007/s00414-016-1411-7
Jobling, M. A., and Tyler-Smith, C. (2017). Human Y-chromosome variation in the genome-sequencing era. Nat. Rev. Genet. 18, 485–497. doi: 10.1038/nrg.2017.36
Kayser, M. (2017). Forensic use of Y-chromosome DNA: a general overview. Hum. Genet. 136, 621–635. doi: 10.1007/s00439-017-1776-9
Kayser, M., Brauer, S., Weiss, G., Underhill, P. A., Roewer, L., Schiefenhovel, W., et al. (2000). Melanesian origin of polynesian Y chromosomes. Curr. Biol. 10, 1237–1246. doi: 10.1016/s0960-9822(00)00734-x
Kayser, M., Vermeulen, M., Knoblauch, H., Schuster, H., Krawczak, M., and Roewer, L. (2007). Relating two deep-rooted pedigrees from central Germany by high-resolution Y-STR haplotyping. Forensic Sci. Int. Genet. 1, 125–128. doi: 10.1016/j.fsigen.2007.02.004
Kim, Y. J., Shin, D. J., Kim, J. M., Jin, H. J., Kwak, K. D., Han, M. S., et al. (2001). Y-chromosome STR haplotype profiling in the Korean population. Forensic Sci. Int. 115, 231–237. doi: 10.1016/s0379-0738(00)00332-7
Kumar, S., Stecher, G., Li, M., Knyaz, C., and Tamura, K. (2018). MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 35, 1547–1549. doi: 10.1093/molbev/msy096
Lang, M., Liu, H., Song, F., Qiao, X., Ye, Y., Ren, H., et al. (2019). Forensic characteristics and genetic analysis of both 27 Y-STRs and 143 Y-SNPs in eastern Han Chinese population. Forensic Sci. Int. Genet. 42, e13–e20. doi: 10.1016/j.fsigen.2019.07.011
Laouina, A., El Houate, B., Yahia, H., Azeddoug, H., Boulouiz, R., and Chbel, F. (2011). Allele frequencies and population data for 17 Y-STR loci (The AmpFlSTR® Y-filerTM) in Casablanca resident population. Forensic Sci. Int. Genet. 5, e1–e3. doi: 10.1016/j.fsigen.2010.10.016
Li, L., Xu, Y., Luis, J. R., Alfonso-Sanchez, M. A., Zeng, Z., Garcia-Bertrand, R., et al. (2019). Cebú, Thailand and Taiwanese aboriginal populations according to Y-STR loci. Gene X 1:100001. doi: 10.1016/j.gene.2018.100001
Li, W., Wang, X., Wang, X., Wang, F., Du, Z., Fu, F., et al. (2020). Forensic characteristics and phylogenetic analyses of one branch of Tai-Kadai language-speaking Hainan Hlai (Ha Hlai) via 23 autosomal STRs included in the Huaxia() platinum system. Mol. Genet. Genomic Med. 8:e1462. doi: 10.1002/mgg3.1462
Liu, J., Wang, R., Shi, J., Cheng, X., Hao, T., Guo, J., et al. (2020). The construction and application of a new 17-plex Y-STR system using universal fluorescent PCR. Int. J. Legal Med. 134, 2015–2027. doi: 10.1007/s00414-020-02291-3
Miranda, J. J., Benecke, M., Hidding, M., and Schmitt, C. (2001). Y-chromosomal short tandem repeat haplotypes at the loci DYS393, DYS19, DYS392, and DYS385-I/II, DYS390, DYS389-I/II, and DYS391 in a Filipino population sample. J. Forensic Sci. 46, 1250–1253.
Mizuno, N., Nakahara, H., Sekiguchi, K., Yoshida, K., Nakano, M., and Kasai, K. (2008). 16 Y chromosomal STR haplotypes in Japanese. Forensic Sci Int 174, 71–76. doi: 10.1016/j.forsciint.2007.01.032
Oppenheimer, S. (2012). Out-of-Africa, the peopling of continents and islands: tracing uniparental gene trees across the map. Philos. Trans. R. Soc. Lond. B Biol. Sci. 367, 770–784. doi: 10.1098/rstb.2011.0306
Ou, X., Wang, Y., Liu, C., Yang, D., Zhang, C., Deng, S., et al. (2015). Haplotype analysis of the polymorphic 40 Y-STR markers in Chinese populations. Forensic Sci. Int. Genet. 19, 255–262. doi: 10.1016/j.fsigen.2015.08.007
Ramos-Luis, E., Blanco-Verea, A., Brión, M., Van Huffel, V., Sánchez-Diz, P., and Carracedo, A. (2014). Y-chromosomal DNA analysis in French male lineages. Forensic Sci. Int. Genet. 9, 162–168. doi: 10.1016/j.fsigen.2013.12.008
Rapone, C., D’Atanasio, E., Agostino, A., Mariano, M., Papaluca, M. T., Cruciani, F., et al. (2016). Forensic genetic value of a 27 Y-STR loci multiplex (Yfiler(§) Plus kit) in an Italian population sample. Forensic Sci. Int. Genet. 21, e1–e5. doi: 10.1016/j.fsigen.2015.11.006
Reid, K. M., and Heathfield, L. J. (2020). Allele frequency data for 23 Y-chromosome short tandem repeats (STRs) for the South African population. Forensic Sci. Int. Genet. 46:102270. doi: 10.1016/j.fsigen.2020.102270
Roewer, L., Andersen, M. M., Ballantyne, J., Butler, J. M., Caliebe, A., Corach, D., et al. (2020). DNA commission of the international society of forensic genetics (ISFG): recommendations on the interpretation of Y-STR results in forensic analysis. Forensic Sci. Int. Genet. 48:102308. doi: 10.1016/j.fsigen.2020.102308
Roewer, L., Croucher, P. J., Willuweit, S., Lu, T. T., Kayser, M., Lessig, R., et al. (2005). Signature of recent historical events in the European Y-chromosomal STR haplotype distribution. Hum. Genet. 116, 279–291. doi: 10.1007/s00439-004-1201-z
Saitou, N., and Nei, M. (1987). The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol. Biol. Evol. 4, 406–425. doi: 10.1093/oxfordjournals.molbev.a040454
Salvador, J., Rodriguez, J., Carandang, L., Agmata, A., Honrado, M., Delfin, F., et al. (2019). Filipino DNA variation at 36 Y-chromosomal short tandem repeat (STR) marker units. Philipp. J. Sci. 148, 43–52.
Shan, W., Ablimit, A., Zhou, W., Zhang, F., Ma, Z., and Zheng, X. (2014). Genetic polymorphism of 17 Y chromosomal STRs in Kazakh and Uighur populations from Xinjiang. China. Int. J. Legal Med. 128, 743–744. doi: 10.1007/s00414-013-0948-y
Shi, M., Bai, R., Bai, L., and Yu, X. (2011). Population genetics for Y-chromosomal STRs haplotypes of Chinese Xibe ethnic group. Forensic Sci. Int. Genet. 5, e119–e121. doi: 10.1016/j.fsigen.2010.08.004
Shu, L., Li, L., Yu, G., Yu, B., Liu, Y., Li, S., et al. (2015). Genetic analysis of 17 Y-STR loci in Han, Dong, Miao and Tujia populations from Hunan province, central-southern China. Forensic Sci. Int. Genet. 19, 250–251. doi: 10.1016/j.fsigen.2015.07.007
Singh, M., Sarkar, A., and Nandineni, M. R. (2018). A comprehensive portrait of Y-STR diversity of Indian populations and comparison with 129 worldwide populations. Sci. Rep. 8:15421. doi: 10.1038/s41598-018-33714-2
Song, M., Wang, Z., Zhang, Y., Zhao, C., Lang, M., Xie, M., et al. (2019). Forensic characteristics and phylogenetic analysis of both Y-STR and Y-SNP in the Li and Han ethnic groups from Hainan island of China. Forensic Sci. Int. Genet. 39, e14–e20. doi: 10.1016/j.fsigen.2018.11.016
SWGDAM (2010). Scientific Working Group on DNA Analysis Interpretation Guidelines for Autosomal STR Typing by Forensic DNA Testing Laboratories.
Underhill, P. A., and Kivisild, T. (2007). Use of y chromosome and mitochondrial DNA population structure in tracing human migrations. Annu. Rev. Genet. 41, 539–564. doi: 10.1146/annurev.genet.41.110306.130407
Wang, C. Z., Su, M. J., Li, Y., Chen, L., Jin, X., Wen, S. Q., et al. (2019). Genetic polymorphisms of 27 Yfiler(®) plus loci in the daur and Mongolian ethnic minorities from Hulunbuir of inner Mongolia autonomous region, China. Forensic Sci. Int. Genet. 40, e252–e255. doi: 10.1016/j.fsigen.2019.02.003
Wang, M., Wang, Z., Zhang, Y., He, G., Liu, J., and Hou, Y. (2017). Forensic characteristics and phylogenetic analysis of two Han populations from the southern coastal regions of China using 27 Y-STR loci. Forensic Sci. Int. Genet. 31, e17–e23. doi: 10.1016/j.fsigen.2017.10.009
Wang, Y., Zhang, Y. J., Zhang, C. C., Li, R., Yang, Y., Ou, X. L., et al. (2016). Genetic polymorphisms and mutation rates of 27 Y-chromosomal STRs in a Han population from Guangdong province, Southern China. Forensic Sci. Int. Genet. 21, 5–9. doi: 10.1016/j.fsigen.2015.09.013
Willuweit, S., and Roewer, L. (2015). The new Y chromosome haplotype reference database. Forensic Sci. Int. Genet. 15, 43–48. doi: 10.1016/j.fsigen.2014.11.024
Wolfgramm Ede, V., Silva, B. C., Aguiar, V. R., Malta, F. S., de Castro, A. M., Ferreira, A. C., et al. (2011). Genetic analysis of 15 autosomal and 12 Y-STR loci in the Espirito Santo state population, Brazil. Forensic Sci. Int. Genet. 5, e41–e43. doi: 10.1016/j.fsigen.2010.05.001
Xie, M., Song, F., Li, J., Lang, M., Luo, H., Wang, Z., et al. (2019). Genetic substructure and forensic characteristics of Chinese Hui populations using 157 Y-SNPs and 27 Y-STRs. Forensic Sci. Int. Genet. 41, 11–18. doi: 10.1016/j.fsigen.2019.03.022
Xu, H., Wang, C. C., Shrestha, R., Wang, L. X., Zhang, M., He, Y., et al. (2015). Inferring population structure and demographic history using Y-STR data from worldwide populations. Mol. Genet. Genomics 290, 141–150. doi: 10.1007/s00438-014-0903-8
Xu, J., Li, L., Wei, L., Nie, Z., Yang, S., Xia, M., et al. (2016). Genetic analysis of 17 Y-STR loci in Han population from Shandong province in East China. Forensic Sci Int Genet 22, e15–e17. doi: 10.1016/j.fsigen.2016.01.016
Yang, Y., Yuan, W., Guo, F., and Jiang, X. (2014). Population data of 17 Y-STR loci in Nanyang Han population from Henan province, central China. Forensic Sci. Int. Genet. 13, 145–146. doi: 10.1016/j.fsigen.2014.07.013
Yao, H. B., Wang, C. C., Tao, X., Shang, L., Wen, S. Q., Zhu, B., et al. (2016). Genetic evidence for an East Asian origin of Chinese Muslim populations Dongxiang and Hui. Sci. Rep. 6:38656. doi: 10.1038/srep38656
Yin, C., Su, K., He, Z., Zhai, D., Guo, K., Chen, X., et al. (2020). Genetic reconstruction and forensic analysis of Chinese Shandong and Yunnan Han populations by Co-analyzing Y chromosomal STRs and SNPs. Genes (Basel) 11:743. doi: 10.3390/genes11070743
Zeng, Z., Rowold, D. J., Garcia-Bertrand, R., Calderon, S., Regueiro, M., Li, L., et al. (2014). Taiwanese aborigines: genetic heterogeneity and paternal contribution to oceania. Gene 542, 240–247. doi: 10.1016/j.gene.2014.03.005
Zhang, J., Wang, J., Liu, Y., Shi, M., Bai, R., and Ma, S. (2017). Haplotype data for 27 Y-chromosomal STR loci in the Chaoshan Han population, South China. Forensic Sci. Int. Genet. 31, e54–e56. doi: 10.1016/j.fsigen.2017.08.003
Zhao, Q., Bian, Y., Zhang, S., Zhu, R., Zhou, W., Gao, Y., et al. (2017). Population genetics study using 26 Y-chromosomal STR loci in the Hui ethnic group in China. Forensic Sci. Int. Genet. 28, e26–e27. doi: 10.1016/j.fsigen.2017.01.018
Zhou, J. (2019). A Preface to “Gaoliang Culture List” in Gaozhou City. Available online at: http://blog.sina.com.cn/s/blog_15029bea60102zec2.html (accessed December 07, 2019).
Keywords: Maoming Han, Gaoliang culture, Y-STR, forensic features, genetic structures
Citation: Fan H, Xie Q, Li Y, Wang L, Wen S-Q and Qiu P (2021) Insights Into Forensic Features and Genetic Structures of Guangdong Maoming Han Based on 27 Y-STRs. Front. Genet. 12:690504. doi: 10.3389/fgene.2021.690504
Received: 03 April 2021; Accepted: 14 May 2021;
Published: 18 June 2021.
Edited by:
Chuan-Chao Wang, Xiamen University, ChinaReviewed by:
Cemal Gurkan, Turkish Cypriot DNA Laboratory, CyprusPankaj Shrivastava, State Forensic Science Laboratory, Sagar, India
Copyright © 2021 Fan, Xie, Li, Wang, Wen and Qiu. 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: Shao-Qing Wen, wenshaoqing@fudan.edu.cn; orcid.org/0000-0003-1223-4720; Pingming Qiu, qiupmfy@126.com; orcid.org/0000-0002-5579-1124; Haoliang Fan, fanhaoliang198931@163.com; orcid.org/0000-0002-3214-0177
†These authors have contributed equally to this work