基于DNA变异的中国汉族人群脱发表型推断及预测模型评估

目的 男性型脱发（male pattern baldness,MPB）,又称为雄激素性脱发（AGA）,是一种常见的男性脱发类型,大约80%的表型差异可以用遗传因素解释。目前的MPB遗传推断研究主要基于欧洲人群,东亚人群相关研究较少。本研究在中国人群中对欧洲人群MPB关联位点进行验证分析,并建立遗传推断模型。方法 本研究调查了486个与欧洲人群MPB相关单核苷酸多态性（SNP）位点在312名中国汉族男性中的关联性,分别使用逐步回归和Lasso回归方法对关联出的位点进行筛选。使用逻辑回归算法构建预测模型,通过十折交叉验证的方法评估。之后进一步比较了逻辑回归、k近邻分类器、随机森林、支持向量机4种常用分类器模型对MPB的预测准确性。结果 有174个SNP位点与中国汉族男性的MPB显著相关（P＜0.05）。通过不同的筛选方法,分别得到了22个SNP和25个SNP的位点集合。基于上述位点集合建立了22-SNP和 25-SNP两种逻辑回归预测模型。以AUC（ROC曲线下方的面积大小,area under curve）来衡量,两种模型对MPB预测的准确性分别为0.85和0.84;经十折交叉验证后预测准确性分别下降至0.81和0.77。当加入年龄作为预测因子后,两种模型的AUC均达到最大值0.89。从运行结果来看,逻辑回归预测模型较本研究中的其他分类器模型具有明显优势。结论 总体而言,虽然预测模型的准确性尚未达到临床期望水平,但SNP在MPB的遗传预测方面仍具备很大的潜力,可以为MPB的早期诊断、临床干预和法庭科学应用提供参考。     

30个祖先信息位点的筛选及应用

摘要：目的 筛选一组祖先信息SNPs位点（AIMs,Ancestry Informative Markers）,构建复合检测体系,用于东亚、欧洲和非洲人群遗传成分描述及个体种族来源推断。方法 以HapMap数据库9个人群的658份样本的分型数据为基础,从30个表型相关基因总共282个SNPs位点中筛选出30个AIMs位点,基于微测序-通用芯片技术构建复合检测体系,并建立人群等位基因频率数据库。使用这组位点分析HapMap数据库中658份人群样本,初步验证位点的区分效能;然后,使用研究构建的体系检验收集的5个人群194份无关个体的DNA样本。最后,通过Structure软件分析获取人群的成分构成以及个体的遗传成分,对个体样本进行种族来源推断。 结果 筛选的30个AIMs位点符合哈迪温伯格平衡（p>0.01）,位点之间没有连锁（r2<0.1）, 658份HapMap数据库样本和194份实验样本的祖先成分分析结果与已知结果完全一致。 结论 本文筛选并建立的30个AIMs位点复合检测体系,能够有效实现东亚、欧洲、非洲人群及混合人群的成分构成和个体遗传成分的分析,有效控制遗传连锁分析中由于人群分层现象带来的误差,也可以用于法医DNA检验中个体祖先来源推断。     

中国人群常见恶性肿瘤遗传风险预测模型研究进展

肿瘤的发生发展是一个涉及多基因、多阶段、多步骤的过程。由于肿瘤表型的高度异质性以及肿瘤病因的复杂性,其风险预测模型构建一直是研究的热点和难点之一。全基因组关联研究(genome wide association studies,GWAS)是分子流行病学研究的重要策略,迄今发现了大量与肿瘤风险相关的易感基因与遗传位点。在传统肿瘤风险预测模型基础上纳入遗传信息,对肿瘤高危人群早期识别、精准预防以及个体化干预具有重大的公共卫生转化意义。文章主要介绍中国人群常见恶性肿瘤遗传风险预测模型的研究现状,以及所面临的机遇和挑战。     

东亚三族群SNP高分辨推断模型构建与效能评估

单核苷酸多态性（single nucleotide polymorphism,SNP）是法医遗传学个体识别和族群推断常用的遗传标记. 本研究集合文献和公共库中祖先信息SNP位点（ancestry informative SNPs,AISNPs）,应用softmax回归、支持向量机和随机森林3种算法,研究东亚北方的3个主体人群（中国北方汉族人、日本人和韩国人）的族群推断效果. 我们分析了来自千人基因组计划的103份中国北方汉族人样本、104份日本人样本和亚洲多样性计划的100份韩国人样本的428个AISNP位点分型,采用多元线性回归共线性诊断筛选出67个高信息量的AISNPs位点组合,构建了softmax回归和支持向量机算法的两种族群推断模型,采用随机森林平均降准分析筛选出42个高信息量的AISNPs位点组合,并构建了随机森林算法的族群推断模型,将softmax回归、支持向量机与随机森林3种模型用于北方汉族人、日本人、韩国人的族群推断,五次十折交叉验证（training∶testing=9∶1）测试3种模型的平均准确率分别为95.19%、95.77%、94.53%. 本研究建立的3种族群推断模型均可用于东亚北方三大人群的遗传推断,42 AISNPs组合的位点数目较少,更适于构建法医检测体系,具有较高的实际应用价值.     

基于毛干蛋白质组的族群推断技术的建立与验证

毛干是一种案件现场常见的生物物证,由于核DNA含量极少且高度降解,难以采用现有的短串联重复序列(short tandem repeat,STR)检验方法进行个人识别鉴定,目前仅使用线粒体DNA检验进行母系亲缘关系的判定,利用率非常低.毛干中蛋白质非常稳定,而且具有遗传多态性,表现为基因组中的非同义单核苷酸多态性(non-synonymous single nucleotide polymorphisms,ns SNPs),转录翻译后形成蛋白质序列中的单氨基酸多态性(single amino acid polymorphisms,SAPs).充分利用毛干蛋白质中蕴含的遗传信息,为案件提供线索和证据,是实际公安业务的迫切需求,具有重要的应用价值.本文选取了104份中国汉族的毛干样本进行蛋白质组的检测,共获得了703个SAP位点,位于460个蛋白质上,共推导出552个nsSNP位点.进一步筛选在所有样本中检出率超过15%的位点,获得了88个nsSNP位点,使用毛干样本对应的口腔拭子DNA对88个ns SNP位点进行一代测序验证.为评估发现的nsSNP位点对于人群的区分能力,以千人数据库(1 000 Genome Project)为参考数据库,采用聚类分析和群体匹配概率等方法对检测的19份毛干样本进行人群来源推断.结果显示,通过检测毛干蛋白质组中的ns SNP可以实现东亚、欧洲、非洲三大洲际人群的区分.     

11q21区段与中国汉族人群Graves病相关性研究

目的:在中国人群中验证Cooper团队在欧洲人群中鉴定的新的Graves病(Graves'disease,GD)易感区段11q21与中国汉族人群GD的相关性。方法:在前期1442例GD患者和1468例正常对照的GWAS数据基础上通过11q21区段精细定位选取主效位点,利用Taqman探针技术进行等位基因分析,在1594例GD患者和1679例正常对照中进行扩大样本验证,然后将两个阶段的结果合并分析并得出11q21区段与GD的相关性结果。结果:验证阶段11q21的rs12575636与GD相关的P值为2.98×10~(-5)(OR=1.42,95%CI=1.20-1.67),GWAS阶段和验证阶段合并后11q21的rs12575636与GD相关的P值达到1.26×10~(-6)(OR=1.35,95%CI=1.19-1.51)。结论:11q21的rs12575636与中国汉族人群GD显著相关,11q21的rs12575636是中国汉族人GD的易感位点。     

中国南北方汉族人群DNA甲基化表观遗传差异研究

目的 族群地域、体貌特征等表型是基因型与环境共同作用的结果。大量基因组学研究表明,汉族人群具有混合特征,内部存在明显的南北遗传差异。本研究旨在探索研究表观基因组在中国南北方汉族人群之间是否存在差异,并筛选差异遗传位点。方法 使用GLINT软件对483份汉族样本的全基因组甲基化芯片数据进行EWAS分析,使用Lasso回归方法筛选位点。使用多元逻辑回归算法构建南北方汉族人群预测模型,通过十折交叉验证的方法评估。结果 筛选出一组南北方汉族之间差异显著的CpG位点,准确性为99.03%,Kappa系数为0.979 6。结论 本研究表明南北方汉族人群之间存在表观遗传差异,本研究为进一步开展不同地域汉族人群之间的表观遗传差异研究奠定了基础。     

骨骼基因分化特征研究

为了揭开人类表型多样性的奥妙,科学家系统分析了骨骼发育相关基因群的群体多态性,并与其他基因进行比较,发现骨骼基因中的SNP位点在非洲人群和非非洲人群（东亚和欧洲人群）之间表现出很高的群体分化特征,但东亚和欧洲人群之间并没有此现象。骨骼基因中的高分化SNP在东亚人群和欧洲人群具有很高的衍生等位基因频率,但非洲人群中却没有。     

NBS1基因SNPs与中国汉族人群原发性肝癌的遗传易感性

为分析DNA损伤修复相关基因NBS1单核苷酸多态性(SNPs)与原发性肝癌遗传易感性的关系,并对高分辨率单链构象多态性(SSCP)检测技术在SNPs分型中的适用性进行评估,本研究对来自中国汉族人群的327例原发性肝癌以及295例阴性对照中NBS1基因常见SNPs的稀有等位基因频率进行检测和分析．此外,对NBS1基因6个常见SNPs分别选择部分样本同时进行直接序列测定,以比较2种方法的检测效果．119例原发性肝癌以及95例肝硬化/慢性肝炎组织标本的SSCP分析结果表明,6个常见NBS1基因SNPs位点(102G>A, 320+208G/A, 553G>C, 1197T>C, 2016A>G和2071-30A>T)中,SNP 1197T>C的稀有等位基因频率为68.1%,显著高于肝硬化/慢性肝炎对照的57.9% (P = 0.0298)．对该SNP位点另外采用208份肝细胞癌和200份健康人群血液标本进一步分析, 肝细胞癌SNP 1197T>C的稀有等位基因频率为66.8%,显著高于健康人群对照的58.8% (P = 0.0170)．其他5个SNPs的稀有等位基因频率在原发性肝癌与肝硬化/慢性肝炎之间均无显著性差异．高分辨率SSCP分析法与直接序列测定法对所选样本的SNPs基因分型结果完全一致,而且直接测序法对PCR扩增产物质量的要求相对高分辨率SSCP分析更高．研究表明,中国汉族人群NBS1基因SNP 1197T>C可能与原发性肝癌的发生相关,高分辨率SSCP技术准确度与直接测序法相当,且操作更加简便易行,非常适用于大量样本多个已知SNPs的基因分型．     

基于机器学习方法对人类基因组DNA双链断裂位点进行识别

DNA双链断裂(double-strand break, DSB)是细胞中一种严重的DNA损伤形式，与包括癌症、重组异常、神经元发育异常在内的多种基因组不稳定性疾病密切相关。由于成本和技术门槛的限制，高通量测序技术绘制的高分辨率DSB图谱十分有限，这阻碍了我们对不同物种基因组中DSB情况的认知。据此，我们建立了以随机森林(RF)、支持向量机(SVM)和逻辑回归(LR)三种分类器为基础算法的分类预测模型，对人类上皮细胞基因组DSB位点进行预测。除了之前预测研究中常用到的表观特征和DNA形状特征外，我们发现DNA序列特征(k-mer频数、GC含量、GC-偏移和互信息)也能表征DSB位点。同时，在考虑DNA物理性质、化学位移和自相关信息后，预测结果得到有效提高。将上述所有特征合并后进行预测，得到了较好的分类预测结果，其中逻辑回归(LR)的分类预测性能是最佳(AUC=0.97),与以往的预测结果相当(AUC=0.964)。另外，通过特征递增搜索方法，得到由294个特征组成的最优特征集，对应的AUC值达到0.974。     

Evaluation of DNA Variants Associated with Androgenetic Alopecia and Their Potential to Predict Male Pattern Baldness

Androgenetic alopecia, known in men as male pattern baldness (MPB), is a very conspicuous condition that is particularly frequent among European men and thus contributes markedly to variation in physical appearance traits amongst Europeans. Recent studies have revealed multiple genes and polymorphisms to be associated with susceptibility to MPB. In this study, 50 candidate SNPs for androgenetic alopecia were analyzed in order to verify their potential to predict MPB. Significant associations were confirmed for 29 SNPs from chromosomes X, 1, 5, 7, 18 and 20. A simple 5-SNP prediction model and an extended 20-SNP model were developed based on a discovery panel of 305 males from various European populations fitting one of two distinct phenotype categories. The first category consisted of men below 50 years of age with significant baldness and the second; men aged 50 years or older lacking baldness. The simple model comprised the five best predictors: rs5919324 near AR, rs1998076 in the 20p11 region, rs929626 in EBF1, rs12565727 in TARDBP and rs756853 in HDAC9. The extended prediction model added 15 SNPs from five genomic regions that improved overall prevalence-adjusted predictive accuracy measured by area under the receiver characteristic operating curve (AUC). Both models were evaluated for predictive accuracy using a test set of 300 males reflecting the general European population. Applying a 65% probability threshold, high prediction sensitivity of 87.1% but low specificity of 42.4% was obtained in men aged <50 years. In men aged ≥50, prediction sensitivity was slightly lower at 67.7% while specificity reached 90%. Overall, the AUC=0.761 calculated for men at or above 50 years of age indicates these SNPs offer considerable potential for the application of genetic tests to predict MPB patterns, adding a highly informative predictive system to the emerging field of forensic analysis of externally visible characteristics.     

基于高密度SNP数据的东亚人群遗传结构研究

目的 东亚疆域辽阔,民族众多,有着广泛多样的语言。中国34个省级行政区可划分为7个地理分区,人群主要分属世界七大语系。已有研究主要集中在东亚人群的起源、迁徙、融合等遗传历史。本文基于5 147份世界人群个体的高密度单核苷酸多态性（SNP）数据,从地域及语言两个角度研究东亚人群尤其是中国人群与世界其他人群的遗传关系,研究中国人群的遗传关系和遗传结构。方法 收集了5 147份世界人群个体的高密度SNP数据,并对其进行质控、合并。通过频率差异分析方法对最终获得的32 789个SNP进行统计学检验,并进一步使用主成分分析、系统发育树、祖先成分分析和D检验统计等方法,对东亚人群与世界其他人群的遗传关系,以及中国人群的遗传关系和遗传结构进行研究。结果 研究发现东亚人群与非洲、美洲和欧洲人群存在显著差异。中国人群可分为7个亚群,不同人群间的遗传聚类与其地理分布、语系语族和族源历史有很强的相关性。结论 本文研究了中国人群与世界人群的遗传关系和差异,并系统研究了中国人群的遗传亚结构。这将丰富东亚人群的群体遗传学、法医遗传学等研究基础,为个体化医疗等工作提供数据支撑。     

Distribution of Two Asian-Related Coding SNPs in the MC1R and OCA2 Genes

Very little is known about the genes and mechanisms affecting skin lightening in Asian populations. In this study, two coding SNPs, c.G1129A (R163Q) at the MC1R (melanocortin 1 receptor) gene and c.A1962G (H615R) at the OCA2 (oculocutaneous albinism type II) gene, were investigated in a total of 1,809 individuals in 16 populations from various areas. The Q163 and R615 alleles prevailed almost exclusively in East and Southeast Asian populations. Wright’s F _ST was 0.445 for R163Q and 0.385 for H615R among the 16 populations. The frequency of the Q163 allele was higher in Northeast Asians than in Southeast Asians. The frequency of the R615 allele was highest in South China and unlikely to be associated with levels of ultraviolet radiation. This allele may be a good marker to study the genetic affinity among East Asians because of its restricted distribution and marked difference in allele frequency.     

Genetic risk factors for BMI and obesity in an ethnically diverse population: Results from the population architecture using genomics and epidemiology (PAGE) study

Objective:

Several genome–wide association studies (GWAS) have demonstrated that common genetic variants contribute to obesity. However, studies of this complex trait have focused on ancestrally European populations, despite the high prevalence of obesity in some minority groups.

Design and Methods:

As part of the “Population Architecture using Genomics and Epidemiology (PAGE)” Consortium, we investigated the association between 13 GWAS‐identified single‐nucleotide polymorphisms (SNPs) and BMI and obesity in 69,775 subjects, including 6,149 American Indians, 15,415 African‐Americans, 2,438 East Asians, 7,346 Hispanics, 604 Pacific Islanders, and 37,823 European Americans. For the BMI‐increasing allele of each SNP, we calculated β coefficients using linear regression (for BMI) and risk estimates using logistic regression (for obesity defined as BMI ≥ 30) followed by fixed‐effects meta‐analysis to combine results across PAGE sites. Analyses stratified by racial/ethnic group assumed an additive genetic model and were adjusted for age, sex, and current smoking. We defined “replicating SNPs” (in European Americans) and “generalizing SNPs” (in other racial/ethnic groups) as those associated with an allele frequency‐specific increase in BMI.

Results:

By this definition, we replicated 9/13 SNP associations (5 out of 8 loci) in European Americans. We also generalized 8/13 SNP associations (5/8 loci) in East Asians, 7/13 (5/8 loci) in African Americans, 6/13 (4/8 loci) in Hispanics, 5/8 in Pacific Islanders (5/8 loci), and 5/9 (4/8 loci) in American Indians.

Conclusion:

Linkage disequilibrium patterns suggest that tagSNPs selected for European Americans may not adequately tag causal variants in other ancestry groups. Accordingly, fine‐mapping in large samples is needed to comprehensively explore these loci in diverse populations.     

The population genomic landscape of human genetic structure,admixture history and local adaptation in Peninsular Malaysia

Peninsular Malaysia is a strategic region which might have played an important role in the initial peopling and subsequent human migrations in Asia. However, the genetic diversity and history of human populations—especially indigenous populations—inhabiting this area remain poorly understood. Here, we conducted a genome-wide study using over 900,000 single nucleotide polymorphisms (SNPs) in four major Malaysian ethnic groups (MEGs; Malay, Proto-Malay, Senoi and Negrito), and made comparisons of 17 world-wide populations. Our data revealed that Peninsular Malaysia has greater genetic diversity corresponding to its role as a contact zone of both early and recent human migrations in Asia. However, each single Orang Asli (indigenous) group was less diverse with a smaller effective population size (N _e) than a European or an East Asian population, indicating a substantial isolation of some duration for these groups. All four MEGs were genetically more similar to Asian populations than to other continental groups, and the divergence time between MEGs and East Asian populations (12,000—6,000 years ago) was also much shorter than that between East Asians and Europeans. Thus, Malaysian Orang Asli groups, despite their significantly different features, may share a common origin with the other Asian groups. Nevertheless, we identified traces of recent gene flow from non-Asians to MEGs. Finally, natural selection signatures were detected in a batch of genes associated with immune response, human height, skin pigmentation, hair and facial morphology and blood pressure in MEGs. Notable examples include SYN3 which is associated with human height in all Orang Asli groups, a height-related gene (PNPT1) and two blood pressure-related genes (CDH13 and PAX5) in Negritos. We conclude that a long isolation period, subsequent gene flow and local adaptations have jointly shaped the genetic architectures of MEGs, and this study provides insight into the peopling and human migration history in Southeast Asia.     

Association of single-nucleotide polymorphisms in RHOB and TXNDC3 with knee osteoarthritis susceptibility: two case-control studies in East Asian populations and a meta-analysis

Introduction

Conflicting findings on the association of single nucleotide polymorphisms (SNPs) in RHOB and TXNDC3 with susceptibility to knee osteoarthritis (OA) have been reported in European Caucasians. To examine the associations of these SNPs with OA in East Asian populations and to evaluate their global significance, we conducted two case-control studies in 955 Chinese and 750 Japanese patients.

Methods

We genotyped the previously implicated SNPs rs585017 (in RHOB) and rs4720262 (in TXNDC3) in patients with primary symptomatic knee OA with radiographic confirmation and in matched control individuals, and analyzed their associations. We further conducted a meta-analysis of the study findings together with those of previously reported European studies using the DerSimonian-Laird procedure.

Results

A significant association of RHOB with knee OA was observed in male Chinese patients (P = 0.02). No significant associations were found for RHOB in any other comparisons in the East Asian populations. The association of TXNDC3 was replicated in Chinese female (P = 0.04) and Japanese (P = 0.03) patients, although none of these associations persisted after Bonferroni correction. Significant association (P = 0.02 for the allelic frequency) with nonsignificant heterogeneity was found in the East Asian replication study. No significant association was found in any comparison in the meta-analysis for all studies.

Conclusion

Our study replicates the association, previously reported in European Caucasians, of TXNDC3 with knee OA susceptibility in an East Asian population.     

Global diversity and genetic contributions of chicken populations from African,Asian and European regions

Genetic diversity and population structure of 113 chicken populations from Africa, Asia and Europe were studied using 29 microsatellite markers. Among these, three populations of wild chickens and nine commercial purebreds were used as reference populations for comparison. Compared to commercial lines and chickens sampled from the European region, high mean numbers of alleles and a high degree of heterozygosity were found in Asian and African chickens as well as in Red Junglefowl. Population differentiation (F_ST) was higher among European breeds and commercial lines than among African, Asian and Red Junglefowl populations. Neighbour‐Net genetic clustering and structure analysis revealed two main groups of Asian and north‐west European breeds, whereas African populations overlap with other breeds from Eastern Europe and the Mediterranean region. Broilers and brown egg layers were situated between the Asian and north‐west European clusters. structure analysis confirmed a lower degree of population stratification in African and Asian chickens than in European breeds. High genetic differentiation and low genetic contributions to global diversity have been observed for single European breeds. Populations with low genetic variability have also shown a low genetic contribution to a core set of diversity in attaining maximum genetic variation present from the total populations. This may indicate that conservation measures in Europe should pay special attention to preserving as many single chicken breeds as possible to maintain maximum genetic diversity given that higher genetic variations come from differentiation between breeds.     

Prediction of Fatness by Standing 8‐Electrode Bioimpedance: A Multiethnic Adolescent Population

The objective of this study was to validate an 8‐electrode bioimpedance analysis (BIA₈) device (BC‐418; Tanita, Tokyo, Japan) for use in populations of European, Maori, Pacific Island, and Asian adolescents. Healthy adolescents (215 M, 216 F; 129 Pacific Island, 120 Asian, 91 Maori, and 91 European; age range 12–19 years) were recruited by purposive sampling of high schools in Auckland, New Zealand. Weight, height, sitting height, leg length, waist circumference, and whole‐body impedance were measured. Fat mass (FM) and fat‐free mass (FFM) derived from the BIA₈ manufacturer's equations were compared with measurements by dual‐energy X‐ray absorptiometry (DXA). DXA‐measured FFM was used as the reference to develop prediction equations based on impedance. A double cross‐validation technique was applied. BIA₈ underestimated FM by 2.06 kg (P < 0.0001) and percent body fat (%BF) by 2.84% (P < 0.0001), on average. However, BIA₈ tended to overestimate FM and %BF in lean and underestimate FM and %BF in fat individuals. Sex‐specific equations developed showed acceptable accuracy on cross‐validation. In the total sample, the best prediction equations were, for boys: FFM (kg) = 0.607 height (cm)²/impedance (Ω) + 1.542 age (y) + 0.220 height (cm) + 0.096 weight (kg) + 1.836 ethnicity (0 = European or Asian, 1 = Maori or Pacific) ? 47.547, R² = 0.93, standard error of estimate (SEE) = 3.09 kg; and, for girls: FFM (kg) = 0.531 height (cm)²/impedance (Ω) + 0.182 height (cm) + 0.096 weight (kg) + 1.562 ethnicity (0 = non‐Pacific, 1 = Pacific) ? 15.782, R² = 0.91, SEE = 2.19 kg. In conclusion, equations for fatness estimation using BIA₈ developed for our sample perform better than reliance on the manufacturer's estimates. The relationship between BIA and body composition in adolescents is ethnicity dependent.