人类群体遗传结构的协方差阵主成分分析方法

目的:探讨基因频率矩阵的中心化(或均值化)协方差阵主成分分析方法在人类群体遗传结构研究中的适用性和合理性。方法:从基因频率矩阵的结构特征入手,分析中心化、均值化协方差阵主成分分析与标准化相关阵主成分分析在特征根、特征向量以及降维效果等方面的差异,并通过实例比较不同方法在解释群体遗传结构特征上合理性。结果:中心化(或均值化)协方差阵的主成分不仅反映了基因变异程度的“方差信息量权”,而且反映了基因间相互影响程度的“相关信息量权”;标准化相关阵的主成分反映的仅是“相关信息量权”,不包括“方差信息量权”。通过比较中国26个汉族人群HLA-A基因座中心化协方差阵和标准化相关阵2种主成分分析结果,证实中心化协方差阵主成分分析方法在特征根与特征向量、保留主成分的个数和对主成分的群体遗传学解释的合理性等方面均优于标准化相关阵主成分分析方法。结论:在对群体遗传结构进行主成分分析时,应使用中心化(或均值化)变换消除基因频率矩阵中量级的影响,然后在用其协方差阵提取主成分。     

人类群体遗传空间结构的"克立格"模型

通过将“克立格”技术应用于人类群体遗传学领域,构建了人类群体遗传空间结构的“克立格”模型,并论述了其原理和计算方法。以HLA-A基因座为例,应用“克立格”模型,定量分析了中国人群HLA-A基因座的空间遗传异质性;对HLA-A基因频率的空间数据矩阵进行了主成分分析,进而定义了人类群体遗传结构的综合遗传测度(SPC),绘制了综合遗传测度和主成分(PC)的“克立格”地图,分析了其群体遗传空间结构特性。与其他空间插值或平滑方法相比,人类群体遗传空间结构的“克立格”模型具有明显优点：1)“克立格”估计以空间遗传变异函数模型为基础,在绘制空间遗传结构地图之前,可利用变异函数模型定量分析所研究基因座(或多基因座)的空间遗传异质性;2)“克立格”插值方法是真正意义上的无偏估计模型,它利用待估区域周围的已知群体遗传调查点数据,并充分考虑调查点的空间影响范围,给出待估区域的最优估计值;3)“克立格”模型允许估计插值误差,这种插值误差既可用于评价空间估计效果,又可通过绘制误差地图指导在误差过高的地点增加新的群体遗传调查样本点,以优化估计效果。然而,人类群体遗传空间结构的“克立格”模型也存在一定缺点：1)若不能用任何理论遗传变异函数模型拟合观察遗传变异函数值,则不能建立“克立格”模型;2)若理论遗传变异函数的拟合优度很低,则据此建立的“克立格”模型的估计标准差在整个空间范围内会很大,此时“克立格”模型不适用于估计群体遗传空间结构。出现上述两种情形时,应选用不考虑空间相关性的空间随机插值方法绘制群体遗传结构地图,如基因绘图软件中的Cavalli-Sforza方法,反向距离加权法和条样函数插值法等。     

群体遗传学研究中STR数据的统计方法应用

STR作为遗传多态性较高的标记, 被广泛地运用于群体遗传学的研究。对于STR分型产生的基因型频率及等位基因频率数据, 文章总结了各种参数指标的计算及分析方法。其中参数指标包括杂合度、多态信息量、连锁不平衡系数、近交系数、遗传距离以及固定指数等; 分析方法包括主成分分析、系统发生树、分子方差分析、R矩阵、地理信息系统以及空间自相关分析。通过这些参数指标及分析方法的使用, 可以既直观又科学地揭示群体遗传结构、群体间遗传分化以及人类起源与进化等群体遗传学中研究的关键问题。     

等位基因多态性群体遗传结构的多元非线性分析方法

长期以来,对于多维基因多态性数据的多元统计分析,如计算遗传距离时昕用的聚类分析、分析群体遗传结构时所用的主成分分析、因子分析和典型相关分析等,一直应用为无约束条件数据而设计的经典多元线性分析方法,并没有注意基因多态性数据的“闭合效应”所带来的问题。从分析基因多态性数据的分布和结构特征入手,文中指出了基因多态性分布具有“闭合数据”的特点,分析了由于“闭合效应”的影响,经典多元线性方法用于群体遗传结构分析昕面临的困难。根据成分数据统计分析的理论和方法,提出了基因多态性群体遗传结构的多元非线性分析基本方法。并以主成分分析为例,通过实例比较和分析了经典线性主成分分析和“对数比”非线性主成分分析的结果,证明“对数比”非线性主成分分析方法是研究基因多态性群体遗传结构的良好方法,具有特异、灵敏等优点,其结果符合群体遗传学规律。     

辽宁满族11个Y-STR基因座多态性及遗传关系的分析

调查了辽宁满族群体11个Y-STR基因座的多态性分布, 探讨其群体遗传学及法医学应用价值。应用Powerplex Y System荧光标记复合扩增系统检测203名满族无关男性个体的11个Y-STR基因座, 用ABI310遗传分析仪进行基因分型, 计算等位基因和单倍型频率, 并结合已公开发表的国内其他16个群体相应基因座的遗传学资料, 分析其遗传距离和聚类关系。满族个体中共检出189种单倍型, 单倍型频率多样性0.9991, 基因多样性GD值在0.4594(DYS391)～0.9258(DYS385a/b); 从遗传距离分析发现, 满族和东北汉族的遗传距离最小(0.0015), 与维吾尔族的遗传距离最大(0.1485)。结果表明, 11个Y-STR基因座在满族群体中具有较好的遗传多态性, 适用于当地区的法医学应用。与其他民族群体遗传多样性的研究, 对了解满族的起源、迁移以及相互关系有重要意义。     

信息

“东亚现代人史前迁徙的遗传学研究”取得新进展东亚现代人的史前迁徙一直是人类学界颇多争议的问题之一。以前的遗传学研究揭示了东亚南方和北方人群在遗传背景上存在一定的差异,但如何解释这种遗传差异以及它同东亚现代人史前迁徙路线的关系是大家争论的焦点。昆明动物研究所“引进国外杰出人才”计划获得者、副所长宿兵研究员领导的研究小组采用父系遗传的Y染色体的遗传标记,对40个东亚南北方代表群体,共2332个男性个体进行了系统的比较分析。他们的研究结果表明,南方群体的Y染色体单倍型较北方群体更为丰富,并存在南方群体特有的单倍型…     

山西汉族17个Y-STR基因座遗传多态性及遗传关系

为了调查山西汉族群体17个Y-STR基因座的多态性分布,探讨其群体遗传学及法医学应用价值,文章应用Y-filer TM试剂盒检测222名山西汉族无关男性个体的17个Y-STR基因座,用ABI3130遗传分析仪进行基因分型,计算等位基因频率及单倍型多样性,并结合已公开发表的国内其他13个群体相关数据资料,分析山西汉族群体遗传距离和聚类关系。结果:山西汉族个体中共检出219种单倍型,单倍型多样性为0.9999;基因多样性GD值在0.3894(DYS391)~0.9755(DYS385a/b)。从遗传距离分析发现,山西汉族与吉黑汉族之间的遗传距离最近(?0.0001),与台湾群体(0.0152)之间的遗传距离相对较远。结果表明该17个Y-STR基因座在山西汉族群体中具有丰富的遗传多态性,对建立Y染色体STR数据库、研究群体遗传学和进行法医学应用有重要意义。     

赣南地区汉族人群29个Y-STR基因座的遗传多样性及系统进化分析

研究赣南地区汉族人群29个Y-STR基因座的遗传多态性以及与国内多个民族群体的遗传关系,探讨其在群体遗传学和法医学中的实际应用价值。采用DNATyperTM Y29试剂盒对1532例赣南地区健康男性无关个体进行Y-STR基因座扩增,3730型遗传分析仪进行毛细管电泳检测,运用GeneMapper IDX1.4软件对数据结果进行Y-STR分型,计算29个Y-STR基因座的等位基因频率及单倍型频率等遗传学参数。应用Mega-X软件对选取的群体构建进化树,用YHRD在线工具软件的分子方差分析(AMOVA),计算群体间遗传距离,同时构建多维尺度图(MDS)。经分析,29个Y-STR基因座在赣南汉族人群中的基因多样性(GD)值范围为0.3815~0.8766,除了DYS508、DYS437、DYS391和DYS438基因座外,其余25个基因座GD值均高于0.5,且单倍型多样性为0.999924,表明29个Y-STR基因座在赣南汉族人群中有较高的遗传多态性。与其他地区汉族人群比较,赣南汉族与福建汉族遗传距离最近(遗传距离Rst值是0.0002),与黑龙江汉族遗传距离最远(Rst值是0.0249);与其他地区少数民族人群比较,赣南汉族与云南白族遗传距离最近(Rst值是0.0059),与甘孜藏族遗传距离最远(Rst值是0.4689)。研究表明,这29个Y-STR基因座在赣南汉族人群中具有较好的遗传多态性分布,能够满足家系排查及法医学检案的要求,所得的数据可为该地区的群体遗传学和法医学研究与应用提供基础数据支持。     

贵州回族、苗族和彝族29个Y-STR基因座遗传多态性及与其他7个群体遗传结构分析

对贵州回族、苗族、彝族开展29个Y-STR基因座遗传多态性及群体遗传结构分析,并与其他7个群体进行遗传关系比较研究,探讨其在法医学和群体遗传学中的实际应用价值.应用DNATyperTMY29试剂盒对309名贵州苗族、331名贵州彝族和291名贵州回族无关个体进行检测,统计29个Y-STR基因座的等位基因频率及基因多样性...     

陕西渭南地区汉族人群17个Y-STR基因座多态性及遗传关系分析

调查陕西渭南地区汉族群体17个Y-STR基因座的多态性,探讨其群体遗传学及法医学应用价值。应用Y-fi ler荧光标记复合扩增系统,对413名陕西渭南地区汉族无关男性个体17个Y-STR基因座进行复合扩增,用ABI3130遗传分析仪进行基因分型,计算各基因座的群体遗传学参数,并结合已经发表的其他10个群体相应基因座的单倍型资料,分析各群体间的遗传距离。413名陕西渭南汉族个体共检出405种单倍型,其中397种单倍型仅出现1次,单倍型多样性达0.9999,基因多样性(GD)为0.4130(DYS391)~0.9734(DYS385a/b),累计GD值为0.9999。遗传距离分析提示,陕西渭南汉族与辽宁满族的遗传距离最小(0.00110),与青海藏族的遗传距离最大(0.22333)。结果表明,17个Y-STR基因座在陕西渭南汉族群体中具有丰富的遗传多态性和较高的非父排除能力,在法医学和人类群体遗传学研究中具有重要价值。     

Population Genetic Studies Revealed Local Adaptation in a High Gene-Flow Marine Fish,the Small Yellow Croaker (Larimichthys polyactis)

The genetic differentiation of many marine fish species is low. Yet local adaptation may be common in marine fish species as the vast and changing marine environment provides more chances for natural selection. Here, we used anonymous as well as known protein gene linked microsatellites and mitochondrial DNA to detect the population structure of the small yellow croaker (Larimichthys polyactis) in the Northwest Pacific marginal seas. Among these loci, we detected at least two microsatellites, anonymous H16 and HSP27 to be clearly under diversifying selection in outlier tests. Sequence cloning and analysis revealed that H16 was located in the intron of BAHCC1 gene. Landscape genetic analysis showed that H16 mutations were significantly associated with temperature, which further supported the diversifying selection at this locus. These marker types presented different patterns of population structure: (i) mitochondrial DNA phylogeny showed no evidence of genetic divergence and demonstrated only one glacial linage; (ii) population differentiation using putatively neutral microsatellites presented a pattern of high gene flow in the L. polyactis. In addition, several genetic barriers were identified; (iii) the population differentiation pattern revealed by loci under diversifying selection was rather different from that revealed by putatively neutral loci. The results above suggest local adaptation in the small yellow croaker. In summary, population genetic studies based on different marker types disentangle the effects of demographic history, migration, genetic drift and local adaptation on population structure and also provide valuable new insights for the design of management strategies in L. polyactis.     

关于最大信息熵原理与群体遗传平衡一致性的探讨

汪小龙等建立了用最大信息熵原理推导一个基因座上群体遗传平衡的统一数学模型,并给出了模型的最大值解,此解正是Hardy-Weinberg平衡定律所给出的基因型频率。这说明当群体基因型信息熵最大时,群体基因型频率不再变化,达到平衡状态,从而证明了最大信息熵原理与Hardy-Weinberg平衡定律具有一致性,同时指出这一结论可以推广至有迁移、突变、选择、遗传漂变、近亲交配的群体以及多个基因座情形。概括地说就是：最大信息熵原理与群体遗传平衡具有一致性。但是,他们仅仅证明了最大信息熵原理与一个基因座上Hardy-Weinberg平衡定律具有一致性,本文在这个范围内将其推广至多个基因座,且每一个基因座均为复等位基因情形。至于最大信息熵原理是否与其它的群体遗传平衡具有一致性,他们的结论仅仅是猜想,并未严格推导。事实上,要想将这种一致性推广到迁移、突变、随机漂变和近亲交配等群体,则不见得正确。      

Population Genetics of Euphydryas Butterflies. I. Genetic Variation and the Neutrality Hypothesis

Twenty-one populations of the checkerspot butterfly, Euphydryas editha, and ten populations of Euphydryas chalcedona were sampled for genetic variation at eight polymorphic enzyme loci. Both species possessed loci that were highly variable from population to population and loci that were virtually identical across all populations sampled. Our data indicate that the neutrality hypothesis is untenable for the loci studied, and therefore selection is indicated as the major factor responsible for producing these patterns. Thorough ecological work allowed gene flow to be ruled out (in almost all instances) as a factor maintaining similar gene frequencies across populations. The Lewontin-Krakauer test indicated magnitudes of heterogeneity among standardized variances of gene frequencies inconsistent with the neutrality hypothesis. The question of whether or not to correct this statistic for sample size is discussed. Observed equitability of gene frequencies of multiple allelic loci was found to be greater than that predicted under the neutrality hypothesis. Genetic differentiation persisting through two generations was found between the one pair of populations known to exchange significant numbers of individuals per generation. Two matrices of genetic distance between populations, based on the eight loci sampled, were found to be significantly correlated with a matrix of environmental distance, based on measures of fourteen environmental parameters. Correlations between gene frequencies and environmental parameters, results of multiple regression analysis, and results of principle component analysis showed strong patterns of association and of "explained" variation. The correlation analyses suggest which factors might be further investigated as proximate selective agents.     

Genetic Differentiation of the Tuva Population with Respect to the AluInsertions

Polymorphism of three rural populations of the Tuva Republic was examined using a set of five autosomal Aluinsertions at the ACE, PLAT, PV92, APOA1,and F13Bloci. The allele frequency distribution patterns revealed in Tuvinians were typical to Mongoloid populations of Asia and were characterized by relatively high frequency of the Alu-repeat insertion at the PV92and F13Bloci along with relatively low insertion frequency at the APOA1locus. With respect to the test systems used, Tuvinian populations examined displayed high levels of genetic diversity. The mean expected heterozygosity values in the populations of Kungurtug, Toora-Khem, and Teeli were 0.433, 0.407, and 0.437, respectively. The level of genetic diversity in the pooled Tuvinian sample was 0.432. The coefficient of genetic differentiation in the three populations studied was 1.4% pointing to relatively low level of genetic subdivision of the indigenous Tuvinian populations. However, estimates of genetic differentiation of the Tuvinian gene pool made by use of the Alu-repeat system were higher compared to those performed using classical protein systems, mtDNA, or Y-chromosomal haplotypes. Even though Tuvinian populations were characterized by common gene pool, some features specific to Western Tuvinian population could be distinguished. These features could be associated with higher contribution of the Caucasian component to the gene pool of this population. Phylogenetic analysis demonstrated close genetic relationships between the Tuvinian and Altaic ethnic populations.     

Change in Genetic Variance under Selection in a Self-Fertilizing Population

In this study we show how the genetic variance of a quantitative trait changes in a self-fertilizing population under repeated cycles of truncation selection, with the analysis based on the infinitesimal model in which it is assumed that the trait is determined by an infinite number of unlinked loci without epistasis. The genetic variance is reduced not as a consequence of the genotypic frequency change but due to the build-up of linkage disequilibrium under truncation selection in this model. We assume that the order of the genotypic contribution from each locus is n(-1/2), where n is the number of loci involved, and investigate the change in linkage disequilibrium resulting from selection and self-fertilization using genotypic frequency dynamics in order to analyze the change in the genetic variance. Our analysis gives recurrence relations of genetic variance among the succeeding generations for the three cases of gene action, i.e., purely additive action, pure dominance without additive effect and the presence of both additive effect and dominance, respectively. Numerical examples are also given as a check on the recurrence formulas.     

Genetic variants influencing lipid levels and risk of dyslipidemia in Chinese population

Recently, several human genetic and genomewide association studies (GWAS) have discovered many genetic loci that are associated with the concentration of the blood lipids. To confirm the reported loci in Chinese population, we conducted a cross-section study to analyse the association of 25 reported SNPs, genotyped by the ABI SNaPshot method, with the blood levels of total cholesterol (TC), low-density lipoprotein cholesterol (LDL-C), high-density lipoprotein cholesterol (HDL-C) and triglycerides (TG) in 1900 individuals by multivariate analysis. Logistic regression was applied to assess the association of the genetic loci with the risk of different types of dyslipidemia. Our study has convincingly identified that 12 of 25 studied SNPs were strongly associated with one or more blood lipid parameters (TC, LDL, HDL and TG). Among the 12 associated SNPs, 10 significantly influence the risk of one or more types of dyslipidemia. We firstly found four SNPs (rs12654264 in HMGCR; rs2479409 in PCSK9; rs16996148 in CILP2, PBX4; rs4420638 in APOE-C1-C4-C2) robustly and independently associate with four types of dyslipidemia (MHL, mixed hyperlipidemia; IHTC, isolated hypercholesterolemia; ILH, isolated low HDL-C; IHTG, isolated hypertriglyceridemia). Our results suggest that genetic susceptibility is different on the same candidate locus for the different populations. Meanwhile, most of the reported genetic variants strongly influence one or more plasma lipid levels and the risk of dyslipidemia in Chinese population.     

Effect of balancing selection on spatial genetic structure within populations: theoretical investigations on the self-incompatibility locus and empirical studies in Arabidopsis halleri

The effect of selection on patterns of genetic structure within and between populations may be studied by contrasting observed patterns at the genes targeted by selection with those of unlinked neutral marker loci. Local directional selection on target genes will produce stronger population genetic structure than at neutral loci, whereas the reverse is expected for balancing selection. However, theoretical predictions on the intensity of this signal under precise models of balancing selection are still lacking. Using negative frequency-dependent selection acting on self-incompatibility systems in plants as a model of balancing selection, we investigated the effect of such selection on patterns of spatial genetic structure within a continuous population. Using numerical simulations, we tested the effect of the type of self-incompatibility system, the number of alleles at the self-incompatibility locus and the dominance interactions among them, the extent of gene dispersal, and the immigration rate on spatial genetic structure at the selected locus and at unlinked neutral loci. We confirm that frequency-dependent selection is expected to reduce the extent of spatial genetic structure as compared to neutral loci, particularly in situations with low number of alleles at the self-incompatibility locus, high frequency of codominant interactions among alleles, restricted gene dispersal and restricted immigration from outside populations. Hence the signature of selection on spatial genetic structure is expected to vary across species and populations, and we show that empirical data from the literature as well as data reported here on three natural populations of the herb Arabidopsis halleri confirm these theoretical results.     

人类混血群体连锁不平衡的遗传学模型与基因定位

人类混血群体可以说是混合群体的一种特例．在无选择、无突变、无限随机交配群体的假定前提下,研究了亲本群体的基因频率对混血群体及其衍生后代群体连锁不平衡结构的影响,导出了各群体连锁不平衡值的表达式,建立了一个估计基因间重组率的简便方法;同时, 采用估算分子标记与QTL之间连锁不平衡系数的统计分析方法,分析了人类混血群体及其衍生后代群体QTL检测与估计的关系,建立了该关系的系列理论公式．研究结果表明,本方法不仅适用于人类疾病(包括复杂遗传疾病)基因定位,而且适合于人类正常基因的定位,同时也适用于人类普通多基因性状的QTL分析．