Abstract
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
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.
TABLE 1
| Number of observed haplotypes | Y-filer (17Y-STR) | Y-filer Plus (27Y-STR) |
| 1 (unique) | 328 | 405 |
| 2 | 33 | 10 |
| 3 | 6 | 2 |
| 4 | 1 | – |
| 5 | 3 | – |
| Sample size | 431 | 431 |
| Number of unique haplotypes | 371 | 417 |
| Proportion of unique haplotypes | 0.8841 | 0.9712 |
| HD | 0.9994 | 0.9997 |
| DC | 0.8608 | 0.9675 |
| RMP | 0.0033 | 0.0027 |
Forensic features of Yfiler and Yfiler Plus loci in 431 Maoming Han.
HD, haplotype diversity; DC, discrimination capacity; RMP, random match probability.
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
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
FIGURE 4
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.
Statements
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).
Acknowledgments
We would like to thank the donors who contributed samples for this study.
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.
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
1
AlamS.AliM. E.FerdousA.HossainT.HasanM. M.AkhteruzzamanS. (2010). Haplotype diversity of 17 Y-chromosomal STR loci in the Bangladeshi population.Forensic Sci. Int. Genet.4e59–e60. 10.1016/j.fsigen.2009.05.005
2
Alonso MoralesL. A.Casas-VargasA.Rojas CastroM.ResqueR.Ribeiro-Dos-SantosÂK.SantosS.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 One13:e0207130. 10.1371/journal.pone.0207130
3
BallantyneK. N.KayserM. (2012). Additional Y-STRs in forensics: why, which, and when.Forensic Sci. Rev.2463–78.
4
BianY.ZhangS.ZhouW.ZhaoQ.Siqintuya, ZhuR.et al (2016). Analysis of genetic admixture in uyghur using the 26 Y-STR loci system.Sci. Rep.6:19998. 10.1038/srep19998
5
CalafellF.LarmuseauM. H. D. (2017). The Y chromosome as the most popular marker in genetic genealogy benefits interdisciplinary research.Hum. Genet.136559–573. 10.1007/s00439-016-1740-0
6
CarracedoA.ButlerJ. M.GusmaoL.LinacreA.ParsonW.RoewerL.et al (2013). New guidelines for the publication of genetic population data.Forensic Sci. Int. Genet.7217–220. 10.1016/j.fsigen.2013.01.001
7
ChenP.HeG.ZouX.ZhangX.LiJ.WangZ.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. 10.1038/s41598-018-33751-x
8
ChenY. (2013). Cultural Business Card in Local Maoming. Nan Fang Lun Kan(09). 87–89.
9
DingJ.FanH.ZhouY.WangZ.WangX.SongX.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.10.1080/20961790.2020.1810882[Epub ahead of print].
10
DingS. (2010). The dialects and cultures in maoming.Forward Position14164–166.
11
DuW.WuW.WuZ.GuoL.WangB.ChenL. (2019). Genetic polymorphisms of 32 Y-STR loci in Meizhou Hakka population.Int. J. Legal Med.133465–466. 10.1007/s00414-018-1845-1
12
FanH.WangX.ChenH.LiW.WangW.DengJ. (2019a). The Ong Be language-speaking population in Hainan island: genetic diversity, phylogenetic characteristics and reflections on ethnicity.Mol. Biol. Rep.464095–4103. 10.1007/s11033-019-04859-8
13
FanH.WangX.ChenH.LongR.LiangA.LiW.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.37e6–e11. 10.1016/j.fsigen.2018.09.008
14
FanH.WangX.ChenH.ZhangX.HuangP.LongR.et al (2018b). Population analysis of 27 Y-chromosomal STRs in the Li ethnic minority from Hainan province, southernmost China.Forensic Sci. Int. Genet.34e20–e22. 10.1016/j.fsigen.2018.01.007
15
FanH.WangX.RenZ.HeG.LongR.LiangA.et al (2019b). Population data of 19 autosomal STR loci in the Li population from Hainan province in southernmost China.Int. J. Legal Med.133429–431. 10.1007/s00414-018-1828-2
16
FanH.ZhangX.WangX.RenZ.LiW.LongR.et al (2018c). Genetic analysis of 27 Y-STR loci in Han population from Hainan province, southernmost China.Forensic Sci. Int. Genet.33e9–e10. 10.1016/j.fsigen.2017.12.009
17
FengR.ZhaoY.ChenS.LiQ.FuY.ZhaoL.et al (2020). Genetic analysis of 50 Y-STR loci in Dong, Miao, Tujia, and Yao populations from Hunan.Int. J. Legal Med.134981–983. 10.1007/s00414-019-02115-z
18
FuX.FuY.LiuY.GuoJ.LiuY.GuoY.et al (2016). Genetic polymorphisms of 26 Y-STR loci in the Mongolian minority from Horqin district. China.Int. J. Legal Med.130941–946. 10.1007/s00414-016-1387-3
19
GaoH. (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.
20
GaoT.YunL.GuY.HeW.WangZ.HouY. (2015). Phylogenetic analysis and forensic characteristics of 12 populations using 23 Y-STR loci.Forensic Sci. Int. Genet.19130–133. 10.1016/j.fsigen.2015.07.006
21
GopinathS.ZhongC.NguyenV.GeJ.LagaceR. E.ShortM. 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.24164–175. 10.1016/j.fsigen.2016.07.006
22
GuanT.SongX.XiaoC.SunH.YangX.LiuC.et al (2020). Analysis of 23 Y-STR loci in Chinese Jieyang Han population.Int. J. Legal Med.134505–507. 10.1007/s00414-019-02019-y
23
GuoF.ZhangL.JiangX. (2015). Population genetics of 17 Y-STR loci in Xibe ethnic minority from Liaoning province, Northeast China.Forensic Sci. Int. Genet.1686–87. 10.1016/j.fsigen.2014.12.007
24
GusmaoL.ButlerJ. M.CarracedoA.GillP.KayserM.MayrW. 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.120191–200. 10.1007/s00414-005-0026-1
25
HanY.LiL.LiuX.ChenW.YangS.WeiL.et al (2016). Genetic analysis of 17 Y-STR loci in Han and Korean populations from Jilin province, Northeast China.Forensic Sci. Int. Genet.228–10. 10.1016/j.fsigen.2016.01.003
26
HeH. (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.
27
HenryJ.DaoH.ScandrettL.TaylorD. (2019). Population genetic analysis of Yfiler(®) Plus haplotype data for three South Australian populations.Forensic Sci. Int. Genet.41e23–e25. 10.1016/j.fsigen.2019.03.021
28
HuL.GuT.FanX.YuanX.RaoM.PangJ. B.et al (2017). Genetic polymorphisms of 24 Y-STR loci in Hani ethnic minority from Yunnan province, Southwest China.Int. J. Legal Med.1311235–1237. 10.1007/s00414-017-1543-4
29
HughesJ. F.PageD. C. (2016). The history of the Y chromosome in man.Nat. Genet.48588–589. 10.1038/ng.3580
30
JiangW.GongZ.RongH.GuanH.ZhangT.ZhaoY.et al (2017). Population genetics of 26 Y-STR loci for the Han ethnic in Hunan province, China.Int. J. Legal Med.131115–117. 10.1007/s00414-016-1411-7
31
JoblingM. A.Tyler-SmithC. (2017). Human Y-chromosome variation in the genome-sequencing era.Nat. Rev. Genet.18485–497. 10.1038/nrg.2017.36
32
KayserM. (2017). Forensic use of Y-chromosome DNA: a general overview.Hum. Genet.136621–635. 10.1007/s00439-017-1776-9
33
KayserM.BrauerS.WeissG.UnderhillP. A.RoewerL.SchiefenhovelW.et al (2000). Melanesian origin of polynesian Y chromosomes.Curr. Biol.101237–1246. 10.1016/s0960-9822(00)00734-x
34
KayserM.VermeulenM.KnoblauchH.SchusterH.KrawczakM.RoewerL. (2007). Relating two deep-rooted pedigrees from central Germany by high-resolution Y-STR haplotyping.Forensic Sci. Int. Genet.1125–128. 10.1016/j.fsigen.2007.02.004
35
KimY. J.ShinD. J.KimJ. M.JinH. J.KwakK. D.HanM. S.et al (2001). Y-chromosome STR haplotype profiling in the Korean population.Forensic Sci. Int.115231–237. 10.1016/s0379-0738(00)00332-7
36
KumarS.StecherG.LiM.KnyazC.TamuraK. (2018). MEGA X: molecular evolutionary genetics analysis across computing platforms.Mol. Biol. Evol.351547–1549. 10.1093/molbev/msy096
37
LangM.LiuH.SongF.QiaoX.YeY.RenH.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.42e13–e20. 10.1016/j.fsigen.2019.07.011
38
LaouinaA.El HouateB.YahiaH.AzeddougH.BoulouizR.ChbelF. (2011). Allele frequencies and population data for 17 Y-STR loci (The AmpFlSTR® Y-filerTM) in Casablanca resident population.Forensic Sci. Int. Genet.5e1–e3. 10.1016/j.fsigen.2010.10.016
39
LiL.XuY.LuisJ. R.Alfonso-SanchezM. A.ZengZ.Garcia-BertrandR.et al (2019). Cebú, Thailand and Taiwanese aboriginal populations according to Y-STR loci.Gene X1:100001. 10.1016/j.gene.2018.100001
40
LiW.WangX.WangX.WangF.DuZ.FuF.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. 10.1002/mgg3.1462
41
LiuJ.WangR.ShiJ.ChengX.HaoT.GuoJ.et al (2020). The construction and application of a new 17-plex Y-STR system using universal fluorescent PCR.Int. J. Legal Med.1342015–2027. 10.1007/s00414-020-02291-3
42
MirandaJ. J.BeneckeM.HiddingM.SchmittC. (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.461250–1253.
43
MizunoN.NakaharaH.SekiguchiK.YoshidaK.NakanoM.KasaiK. (2008). 16 Y chromosomal STR haplotypes in Japanese.Forensic Sci Int17471–76. 10.1016/j.forsciint.2007.01.032
44
OppenheimerS. (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.367770–784. 10.1098/rstb.2011.0306
45
OuX.WangY.LiuC.YangD.ZhangC.DengS.et al (2015). Haplotype analysis of the polymorphic 40 Y-STR markers in Chinese populations.Forensic Sci. Int. Genet.19255–262. 10.1016/j.fsigen.2015.08.007
46
Ramos-LuisE.Blanco-VereaA.BriónM.Van HuffelV.Sánchez-DizP.CarracedoA. (2014). Y-chromosomal DNA analysis in French male lineages.Forensic Sci. Int. Genet.9162–168. 10.1016/j.fsigen.2013.12.008
47
RaponeC.D’AtanasioE.AgostinoA.MarianoM.PapalucaM. T.CrucianiF.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.21e1–e5. 10.1016/j.fsigen.2015.11.006
48
ReidK. M.HeathfieldL. J. (2020). Allele frequency data for 23 Y-chromosome short tandem repeats (STRs) for the South African population.Forensic Sci. Int. Genet.46:102270. 10.1016/j.fsigen.2020.102270
49
RoewerL.AndersenM. M.BallantyneJ.ButlerJ. M.CaliebeA.CorachD.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. 10.1016/j.fsigen.2020.102308
50
RoewerL.CroucherP. J.WilluweitS.LuT. T.KayserM.LessigR.et al (2005). Signature of recent historical events in the European Y-chromosomal STR haplotype distribution.Hum. Genet.116279–291. 10.1007/s00439-004-1201-z
51
SaitouN.NeiM. (1987). The neighbor-joining method: a new method for reconstructing phylogenetic trees.Mol. Biol. Evol.4406–425. 10.1093/oxfordjournals.molbev.a040454
52
SalvadorJ.RodriguezJ.CarandangL.AgmataA.HonradoM.DelfinF.et al (2019). Filipino DNA variation at 36 Y-chromosomal short tandem repeat (STR) marker units.Philipp. J. Sci.14843–52.
53
ShanW.AblimitA.ZhouW.ZhangF.MaZ.ZhengX. (2014). Genetic polymorphism of 17 Y chromosomal STRs in Kazakh and Uighur populations from Xinjiang. China.Int. J. Legal Med.128743–744. 10.1007/s00414-013-0948-y
54
ShiM.BaiR.BaiL.YuX. (2011). Population genetics for Y-chromosomal STRs haplotypes of Chinese Xibe ethnic group.Forensic Sci. Int. Genet.5e119–e121. 10.1016/j.fsigen.2010.08.004
55
ShuL.LiL.YuG.YuB.LiuY.LiS.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.19250–251. 10.1016/j.fsigen.2015.07.007
56
SinghM.SarkarA.NandineniM. R. (2018). A comprehensive portrait of Y-STR diversity of Indian populations and comparison with 129 worldwide populations.Sci. Rep.8:15421. 10.1038/s41598-018-33714-2
57
SongM.WangZ.ZhangY.ZhaoC.LangM.XieM.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.39e14–e20. 10.1016/j.fsigen.2018.11.016
58
SWGDAM (2010). Scientific Working Group on DNA Analysis Interpretation Guidelines for Autosomal STR Typing by Forensic DNA Testing Laboratories.
59
UnderhillP. A.KivisildT. (2007). Use of y chromosome and mitochondrial DNA population structure in tracing human migrations.Annu. Rev. Genet.41539–564. 10.1146/annurev.genet.41.110306.130407
60
WangC. Z.SuM. J.LiY.ChenL.JinX.WenS. 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.40e252–e255. 10.1016/j.fsigen.2019.02.003
61
WangM.WangZ.ZhangY.HeG.LiuJ.HouY. (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.31e17–e23. 10.1016/j.fsigen.2017.10.009
62
WangY.ZhangY. J.ZhangC. C.LiR.YangY.OuX. 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.215–9. 10.1016/j.fsigen.2015.09.013
63
WilluweitS.RoewerL. (2015). The new Y chromosome haplotype reference database.Forensic Sci. Int. Genet.1543–48. 10.1016/j.fsigen.2014.11.024
64
Wolfgramm EdeV.SilvaB. C.AguiarV. R.MaltaF. S.de CastroA. M.FerreiraA. 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.5e41–e43. 10.1016/j.fsigen.2010.05.001
65
XieM.SongF.LiJ.LangM.LuoH.WangZ.et al (2019). Genetic substructure and forensic characteristics of Chinese Hui populations using 157 Y-SNPs and 27 Y-STRs.Forensic Sci. Int. Genet.4111–18. 10.1016/j.fsigen.2019.03.022
66
XuH.WangC. C.ShresthaR.WangL. X.ZhangM.HeY.et al (2015). Inferring population structure and demographic history using Y-STR data from worldwide populations.Mol. Genet. Genomics290141–150. 10.1007/s00438-014-0903-8
67
XuJ.LiL.WeiL.NieZ.YangS.XiaM.et al (2016). Genetic analysis of 17 Y-STR loci in Han population from Shandong province in East China.Forensic Sci Int Genet22e15–e17. 10.1016/j.fsigen.2016.01.016
68
YangY.YuanW.GuoF.JiangX. (2014). Population data of 17 Y-STR loci in Nanyang Han population from Henan province, central China.Forensic Sci. Int. Genet.13145–146. 10.1016/j.fsigen.2014.07.013
69
YaoH. B.WangC. C.TaoX.ShangL.WenS. Q.ZhuB.et al (2016). Genetic evidence for an East Asian origin of Chinese Muslim populations Dongxiang and Hui.Sci. Rep.6:38656. 10.1038/srep38656
70
YinC.SuK.HeZ.ZhaiD.GuoK.ChenX.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. 10.3390/genes11070743
71
ZengZ.RowoldD. J.Garcia-BertrandR.CalderonS.RegueiroM.LiL.et al (2014). Taiwanese aborigines: genetic heterogeneity and paternal contribution to oceania.Gene542240–247. 10.1016/j.gene.2014.03.005
72
ZhangJ.WangJ.LiuY.ShiM.BaiR.MaS. (2017). Haplotype data for 27 Y-chromosomal STR loci in the Chaoshan Han population, South China.Forensic Sci. Int. Genet.31e54–e56. 10.1016/j.fsigen.2017.08.003
73
ZhaoQ.BianY.ZhangS.ZhuR.ZhouW.GaoY.et al (2017). Population genetics study using 26 Y-chromosomal STR loci in the Hui ethnic group in China.Forensic Sci. Int. Genet.28e26–e27. 10.1016/j.fsigen.2017.01.018
74
ZhouJ. (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).
75
ZhuB.LiX.WangZ.WuH.HeY.ZhaoJ.et al (2005). Y-STRs haplotypes of Chinese Mongol ethnic group using Y-PLEX 12.Forensic Sci. Int.153260–263. 10.1016/j.forsciint.2004.11.002
Summary
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
Volume
12 - 2021
Edited by
Chuan-Chao Wang, Xiamen University, China
Reviewed by
Cemal Gurkan, Turkish Cypriot DNA Laboratory, Cyprus; Pankaj Shrivastava, State Forensic Science Laboratory, Sagar, India
Updates
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-4720Pingming Qiu, qiupmfy@126.com; orcid.org/0000-0002-5579-1124Haoliang Fan, fanhaoliang198931@163.com; orcid.org/0000-0002-3214-0177
†These authors have contributed equally to this work
This article was submitted to Evolutionary and Population Genetics, a section of the journal Frontiers in Genetics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.