Common SNPs explain a large proportion of the heritability for human height

SNPs discovered by genome-wide association studies (GWASs) account for only a small fraction of the genetic variation of complex traits in human populations. Where is the remaining heritability? We estimated the proportion of variance for human height explained by 294,831 SNPs genotyped on 3,925 unrelated individuals using a linear model analysis, and validated the estimation method with simulations based on the observed genotype data. We show that 45% of variance can be explained by considering all SNPs simultaneously. Thus, most of the heritability is not missing but has not previously been detected because the individual effects are too small to pass stringent significance tests. We provide evidence that the remaining heritability is due to incomplete linkage disequilibrium between causal variants and genotyped SNPs, exacerbated by causal variants having lower minor allele frequency than the SNPs explored to date.

Sequence variants in three loci influence monocyte counts and erythrocyte volume

Blood cells participate in vital physiological processes, and their numbers are tightly regulated so that homeostasis is maintained. Disruption of key regulatory mechanisms underlies many blood-related Mendelian diseases but also contributes to more common disorders, including atherosclerosis. We searched for quantitative trait loci (QTL) for hematology traits through a whole-genome association study, because these could provide new insights into both hemopoeitic and disease mechanisms. We tested 1.8 million variants for association with 13 hematology traits measured in 6015 individuals from the Australian and Dutch populations. These traits included hemoglobin composition, platelet counts, and red blood cell and white blood cell indices. We identified three regions of strong association that, to our knowledge, have not been previously reported in the literature. The first was located in an intergenic region of chromosome 9q31 near LPAR1, explaining 1.5% of the variation in monocyte counts (best SNP rs7023923, p=8.9x10(-14)). The second locus was located on chromosome 6p21 and associated with mean cell erythrocyte volume (rs12661667, p=1.2x10(-9), 0.7% variance explained) in a region that spanned five genes, including CCND3, a member of the D-cyclin gene family that is involved in hematopoietic stem cell expansion. The third region was also associated with erythrocyte volume and was located in an intergenic region on chromosome 6q24 (rs592423, p=5.3x10(-9), 0.6% variance explained). All three loci replicated in an independent panel of 1543 individuals (p values=0.001, 9.9x10(-5), and 7x10(-5), respectively). The identification of these QTL provides new opportunities for furthering our understanding of the mechanisms regulating hemopoietic cell fate.

Variants in EMX2 and PTEN do not contribute to risk of endometriosis

Endometriosis has a genetic component, and significant linkage has been found to a region on chromosome 10q. Two candidate genes, EMX2 and PTEN, implicated in both endometriosis and endometrial cancer, lie on chromosome 10q. We hypothesized that variation in EMX2 and/or PTEN could contribute to the risk of endometriosis and may account for some of the linkage signal on 10q. We genotyped single nucleotide polymorphisms (SNPs) in a case-control design to evaluate association between endometriosis and common variations in these two genes. The genotyping and statistical analysis were based on samples collected from Australian volunteers. The cases were 768 unrelated women with surgically confirmed endometriosis selected from affected sister pair (ASP) families participating in the Australian Genes behind Endometriosis Study. The controls were 768 female participants in twin studies who, based on screening questions, did not have a diagnosis of endometriosis. Genotypes of 22 SNPs in the EMX2 gene and 15 SNPs in the PTEN gene were the main outcome measures. Statistical analysis provided measures of linkage disequilibrium and association. Permutation testing showed no globally significant association between any SNPs or haplotypes and endometriosis for either gene. It is unlikely that the EMX2 or PTEN gene variants investigated contribute to risk for initiation and/or development of endometriosis.

Separate and interacting effects within the catechol-O-methyltransferase (COMT) are associated with schizophrenia

Several lines of evidence have implicated the catechol-O-methyltransferase (COMT) gene as a candidate for schizophrenia (SZ) susceptibility, not only because it encodes a key dopamine catabolic enzyme but also because it maps to the velocardiofacial syndrome region of chromosome 22q11 which has long been associated with SZ predisposition. The interest in COMT as a candidate SZ risk factor has led to numerous case-control and family-based studies, with the majority placing emphasis on examining a functional Val/Met polymorphism within this enzyme. Unfortunately, these studies have continually produced conflicting results. To assess the genetic contribution of other COMT variants to SZ susceptibility, we investigated three single-nucleotide polymorphisms (SNPs) (rs737865, rs4633, rs165599) in addition to the Val/Met variant (rs4680) in a highly selected sample of Australian Caucasian families containing 107 patients with SZ. The Val/Met and rs4633 variants showed nominally significant associations with SZ (P<0.05), although neither of the individual SNPs remained significant after adjusting for multiple testing (most significant P=0.1174). However, haplotype analyses showed strong evidence of an association; the most significant being the three-marker haplotype rs737865-rs4680-rs165599 (global P=0.0022), which spans more than 26 kb. Importantly, conditional analyses indicated the presence of two separate and interacting effects within this haplotype, irrespective of gender. In addition, our results indicate the Val/Met polymorphism is not disease-causing and is simply in strong linkage disequilibrium with a causative effect, which interacts with another as yet unidentified variant approximately 20 kb away. These results may help explain the inconsistent results reported on the Val/Met polymorphism and have important implications for future investigations into the role of COMT in SZ susceptibility.

GENEHUNTER: your 'one-stop shop' for statistical genetic analysis?

The past decade has brought a proliferation of statistical genetic (linkage) analysis techniques, incorporating new methodology and/or improvement of existing methodology in gene mapping, specifically targeted towards the localization of genes underlying complex disorders. Most of these techniques have been implemented in user-friendly programs and made freely available to the genetics community. Although certain packages may be more 'popular' than others, a common question asked by genetic researchers is 'which program is best for me?'. To help researchers answer this question, the following software review aims to summarize the main advantages and disadvantages of the popular GENEHUNTER package.

Further evidence for linkage of Gilles de la Tourette syndrome (GTS) susceptibility loci on chromosomes 2p11, 8q22 and 11q23-24 in South African Afrikaners

Utilizing DNA samples from 91 Afrikaner nuclear families with one or more affected children, five genomic regions on chromosomes 2p, 8q, 11q, 20q, and 21q that gave evidence for association with GTS in previous case-control association studies were investigated for linkage and association with GTS. Highly polymorphic markers with mean heterozygosity of 0.77 were typed and resulting genotypes evaluated using single marker transmission disequilibrium (TDT), single marker haplotype relative risk (HRR), and multi-marker "extended" TDT and HRR methods. Single marker TDT analysis showed evidence for linkage or association, with p-values near 0.05, for markers D2S139, GATA28F12, and D11S1377 on chromosomes 2p11, 8q22 and 11q23-24, respectively. Extended, two-locus TDT and HRR analysis provided further evidence for linkage or association on chromosome 2 with p-values of 0.007 and 0.025, and chromosome 8 with p-values of 0.059 and 0.013, respectively. These results provide important additional evidence for the location of GTS susceptibility loci.

Genetic Structure of Cochliobolus sativus Populations Sampled from Roots and Leaves of Barley and Wheat in North Dakota

Common root rot (CRR) and spot blotch, caused by Cochliobolus sativus (Ito and Kurib.) Drechsl. ex Dast., are important diseases of barley (Hordeum vulgare L.) and wheat (Triticum aestivum L.) worldwide. However, the population biology of C. sativus is still poorly understood. In this study, the genetic structure of three C. sativus populations, consisting of isolates sampled respectively from barley leaves (BL), barley roots (BR) and wheat roots (WR) in North Dakota, was analysed with amplified fragment length polymorphism (AFLP) markers. A total of 127 AFLP loci were generated among 208 C. sativus isolates analysed with three primer combinations. Gene diversity (H = 0.277-0.335) were high in all three populations. Genetic variation among C. sativus individuals within population accounted for 74%, whereas 26% of the genetic variation was explained among populations. Genetic differentiation was high (empty set PT = 0.261, corrected G ''(st)= 0.39), whereas gene flow (Nm) ranged from 1.27 to 1.56 among the three populations analysed. The multilocus linkage disequilibrium (LD) ((r) over bard = 0.0760.117) was moderate in C. sativus populations. Cluster analyses indicate that C. sativus populations differentiated according to the hosts (barley and wheat) and tissues (root and leaf) although generalists also exist in North Dakota. Crop breeding may benefit from combining genes for resistance against both specialists and generalists of C. sativus.

Association genetics in Corymbia citriodora subsp variegata identifies single nucleotide polymorphisms affecting wood growth and cellulosic pulp yield

Wood is an important biological resource which contributes to nutrient and hydrology cycles through ecosystems, and provides structural support at the plant level. Thousands of genes are involved in wood development, yet their effects on phenotype are not well understood. We have exploited the low genomic linkage disequilibrium (LD) and abundant phenotypic variation of forest trees to explore allelic diversity underlying wood traits in an association study. Candidate gene allelic diversity was modelled against quantitative variation to identify SNPs influencing wood properties, growth and disease resistance across three populations of Corymbia citriodora subsp. variegata, a forest tree of eastern Australia. Nine single nucleotide polymorphism (SNP) associations from six genes were identified in a discovery population (833 individuals). Associations were subsequently tested in two smaller populations (130160 individuals), validating our findings in three cases for actin 7 (ACT7) and COP1 interacting protein 7 (CIP7). The results imply a functional role for these genes in mediating wood chemical composition and growth, respectively. A flip in the effect of ACT7 on pulp yield between populations suggests gene by environment interactions are at play. Existing evidence of gene function lends strength to the observed associations, and in the case of CIP7 supports a role in cortical photosynthesis.

Genetic Structure of Cochliobolus sativus Populations Sampled from Roots and Leaves of Barley and Wheat in North Dakota

Common root rot (CRR) and spot blotch, caused by Cochliobolus sativus (Ito and Kurib.) Drechsl. ex Dast., are important diseases of barley (Hordeum vulgare L.) and wheat (Triticum aestivum L.) worldwide. However, the population biology of C. sativus is still poorly understood. In this study, the genetic structure of three C. sativus populations, consisting of isolates sampled respectively from barley leaves (BL), barley roots (BR) and wheat roots (WR) in North Dakota, was analysed with amplified fragment length polymorphism (AFLP) markers. A total of 127 AFLP loci were generated among 208 C. sativus isolates analysed with three primer combinations. Gene diversity (H = 0.277-0.335) were high in all three populations. Genetic variation among C. sativus individuals within population accounted for 74%, whereas 26% of the genetic variation was explained among populations. Genetic differentiation was high (empty set PT = 0.261, corrected G ''(st)= 0.39), whereas gene flow (Nm) ranged from 1.27 to 1.56 among the three populations analysed. The multilocus linkage disequilibrium (LD) ((r) over bard = 0.0760.117) was moderate in C. sativus populations. Cluster analyses indicate that C. sativus populations differentiated according to the hosts (barley and wheat) and tissues (root and leaf) although generalists also exist in North Dakota. Crop breeding may benefit from combining genes for resistance against both specialists and generalists of C. sativus.

Genetics of T cell co-stimulatory receptors -CD28, CTLA4, ICOS and PDCD1 in immunity and transplantation

Co-stimulatory signals are essential for the activation of naïve T cells and productive immune response. Naïve T cells receive first, antigen-specific signal through T cell receptor. Co-stimulatory receptors provide the second signal which can be either activating or inhibitory. The balance between signals determines the outcome of an immune response. CD28 is crucial for T cell activation; whereas cytotoxic T lymphocyte associated antigen 4 (CTLA4) mediates critical inhibitory signal. Inducible co-stimulator (ICOS) augments cytokine expression and plays role in immunoglobulin class switching. Programmed cell death 1 (PDCD1) acts as negative regulator of T cell proliferation and cytokine responses. The co-stimulatory receptor pathways are potentially involved in self-tolerance and thus, they provide a promising therapeutic strategy for autoimmune diseases and transplantation. The genes encoding CD28, CTLA4 and ICOS are located adjacently in the chromosome region 2q33. The PDCD1 gene maps further, to the region 2q37. CTLA4 and PDCD1 are associated with the risk of a few autoimmune diseases. There is strong linkage disequilibrium (LD) on the 2q33 region; the whole gene of CD28 exists in its own LD block but CTLA4 and the 5' part of ICOS are within a same LD block. The 3' part of ICOS and PDCD1 are in their own separate LD blocks. Extended haplotypes covering the 2q33 region can be identified. This study focuses on immune related conditions like coeliac disease (CD) which is a chronic inflammatory disease with autoimmune features. Immunoglobulin A deficiency (IgAD) belongs to the group of primary antibody deficiencies characterised by reduced levels of immunoglobulins. IgAD co-occurs often with coeliac disease. Renal transplantation is needed in the end stage kidney diseases. Transplantation causes strong immune response which is tried to suppress with drugs. All these conditions are multifactorial with complex genetic background and multiple environmental factors affecting the outcome. We have screened ICOS for polymorphisms by sequencing the exon regions. We detected 11 new variants and determined their frequencies in Finnish population. We have measured linkage disequilibrium on the 2q33 region in Finnish as well as other European populations and observed conserved haplotypes. We analysed genetic association and linkage of the co-stimulatory receptor gene region aiming to study if it is a common risk locus for immune diseases. The 2q33 region was replicated to be linked to coeliac disease in Finnish population and CTLA4-ICOS haplotypes were found to be associated with CD and IgAD being the first non-HLA risk locus common for CD and immunodeficiencies. We also showed association between ICOS and the outcome of kidney transplantation. Our results suggest new evidence for CTLA4-ICOS gene region to be involved in susceptibility of coeliac disease. The earlier published contradictory association results can be explained by involvement of both CTLA4 and ICOS in disease susceptibility. The pattern of variants acting together rather than a single polymorphism may confer the disease risk. These genes may predispose also to immunodeficiencies as well as decreased graft survival and delayed graft function. Consequently, the present study indicates that like the well established HLA locus, the co-stimulatory receptor genes predispose to variety of immune disorders.

Genetic analysis of chromosomal regions 2q33, 7q32 and 19q13 in multiple sclerosis susceptibility

Multiple sclerosis (MS) is an immune-mediated demyelinating disorder of the central nervous system (CNS) affecting 0.1-0.2% of Northern European descent population. MS is considered to be a multifactorial disease, both environment and genetics play a role in its pathogenesis. Despite several decades of intense research, the etiological and pathogenic mechanisms underlying MS remain still largely unknown and no curative treatment exists. The genetic architecture underlying MS is complex with multiple genes involved. The strongest and the best characterized predisposing genetic factors for MS are located, as in other immune-mediated diseases, in the major histocompatibility complex (MHC) on chromosome 6. In humans MHC is called human leukocyte antigen (HLA). Alleles of the HLA locus have been found to associate strongly with MS and remained for many years the only consistently replicable genetic associations. However, recently other genes located outside the MHC region have been proposed as strong candidates for susceptibility to MS in several studies. In this thesis a new genetic locus located on chromosome 7q32, interferon regulatory factor 5 (IRF5), was identified in the susceptibility to MS. In particular, we found that common variation of the gene was associated with the disease in three different populations, Spanish, Swedish and Finnish. We also suggested a possible functional role for one of the risk alleles with impact on the expression of the IRF5 locus. Previous studies have pointed out a possible role played by chromosome 2q33 in the susceptibility to MS and other autoimmune disorders. The work described here also investigated the involvement of this chromosomal region in MS predisposition. After the detection of genetic association with 2q33 (article-1), we extended our analysis through fine-scale single nucleotide polymorphism (SNP) mapping to define further the contribution of this genomic area to disease pathogenesis (article-4). We found a trend (p=0.04) for association to MS with an intronic SNP located in the inducible T-cell co-stimulator (ICOS) gene, an important player in the co-stimulatory pathway of the immune system. Expression analysis of ICOS revealed a novel, previously uncharacterized, alternatively spliced isoform, lacking the extracellular domain that is needed for ligand binding. The stability of the newly-identified transcript variant and its subcellular localization were analyzed. These studies indicated that the novel isoform is stable and shows different subcellular localization as compared to full-length ICOS. The novel isoform might have a regulatory function, but further studies are required to elucidate its function. Chromosome 19q13 has been previously suggested as one of the genomic areas involved in MS predisposition. In several populations, suggestive linkage signals between MS predisposition and 19q13 have been obtained. Here, we analysed the role of allelic variation in 19q13 by family based association analysis in 782 MS families collected from Finland. In this dataset, we were not able to detect any statistically significant associations, although several previously suggested markers were included to the analysis. Replication of the previous findings on the basis of linkage disequilibrium between marker allele and disease/risk allele appears notoriously difficult because of limitations such as allelic heterogeneity. Re-sequencing based approaches may be required for elucidating the role of chromosome 19q13 with MS. This thesis has resulted in the identification of a new MS susceptibility locus (IRF5) previously associated with other inflammatory or autoimmune disorders, such as SLE. IRF5 is one of the mediators of interferons biological function. In addition to providing new insight in the possible pathogenetic pathway of the disease, this finding suggests that there might be common mechanisms between different immune-mediated disorders. Furthermore the work presented here has uncovered a novel isoform of ICOS, which may play a role in regulatory mechanisms of ICOS, an important mediator of lymphocyte activation. Further work is required to uncover its functions and possible involvement of the ICOS locus in MS susceptibility.

Biological mechanisms behind clozapine-induced agranulocytosis

Clozapine is the most effective drug in treating therapy-resistant schizophrenia and may even be superior to all other antipsychotics. However, its use is limited by a high incidence (approximately 0.8%) of a severe hematological side effect, agranulocytosis. The exact molecular mechanism(s) of clozapine-induced agranulocytosis is still unknown. We investigated the mechanisms behind responsiveness to clozapine therapy and the risk of developing agranulocytosis by performing an HLA (human leukocyte antigens) association study in patients with schizophrenia. The first group comprised patients defined by responsiveness to first-generation antipsychotics (FGAs) (n= 19). The second group was defined by a lack of response to FGAs but responsiveness to clozapine (n=19). The third group of patients had a history of clozapine-induced granulocytopenia or agranulocytosis (n=26). Finnish healthy blood donors served as controls (n= 120). We found a significantly increased frequency of HLA-A1 among patients who were refractory to FGAs but responsive to clozapine. We also found that the frequency of HLA-A1 was low in patients with clozapine-induced neutropenia or agranulocytosis. These results suggest that HLA-A1 may predict a good therapeutic outcome and a low risk of agranulocytosis and therefore HLA typing may aid in the selection of patients for clozapine therapy. Furthermore, in a subgroup of schizophrenia, HLA-A1 may be in linkage disequilibrium with some vulnerability genes in the MHC (major histocompatibility complex) region on chromosome 6. These genes could be involved in antipsychotic drug response and clozapine-induced agranulocytosis. In addition, we investigated the effect of clozapine on gene expression in granulocytes by performing a microarray analysis on blood leukocytes of 8 schizophrenic patients who had started clozapine therapy for the first time. We identified an altered expression in 4 genes implicated in the maturation or apoptosis of granulocytes: MPO (myeloperoxidase precursor), MNDA (myeloid cell nuclear differentiation antigen), FLT3LG (Fms-related tyrosine kinase 3 ligand) and ITGAL (antigen CD11A, lymphocyte function-associated antigen 1). The altered expression of these genes following clozapine administration may suggest their involvement in clozapine-induced agranulocytosis. Finally, we investigated whether or not normal human bone marrow mesenchymal stromal cells (MSC) are sensitive to clozapine. We treated cultures of human MSCs and human skin fibroblasts with 10 µM of unmodified clozapine and with clozapine bioactivated by oxidation. We found that, independent of bioactivation, clozapine was cytotoxic to MSCs in primary culture, whereas clozapine at the same concentration stimulated the growth of human fibroblasts. This suggests that direct cytotoxicity to MSCs is one possible mechanism by which clozapine induces agranulocytosis.

Evidence for an association with type 2 diabetes mellitus at the PPARG locus in a South Indian population

Peroxisome proliferator activated receptor-gamma 2 (PPARG2) is a nuclear hormone receptor of ligand-dependent ranscription factor involved in adipogenesis and a molecular target of the insulin sensitizers thiazolidinediones. We addressed the question of whether the 3 variants (-1279G/A, Pro12Ala, and His478His) in the PPARG2 gene are associated with type 2 diabetes mellitus and its related traits in a South Indian population. The study subjects (1000 type 2 diabetes mellitus and 1000 normal glucose-tolerant subjects) were chosen randomly from the Chennai Urban Rural Epidemiology Study, an ongoing population-based study in southern India. The variants were screened by single-stranded conformational variant, direct sequencing, and restriction fragment length polymorphism. Linkage disequilibrium was estimated from the estimates of haplotypic frequencies. The -1279G/A, Pro12Ala, and His478His variants of the PPARG2 gene were not associated with type 2 diabetes mellitus. However, the 2-loci analyses showed that, in the presence of Pro/Pro genotype of the Pro12Ala variant, the -1279G/A promoter variant showed increased susceptibility to type 2 diabetes mellitus (odds ratio, 2.092; 95% confidence interval, 1.22-3.59; P = .008), whereas in the presence of 12Ala allele, the -1279G/A showed a protective effect against type 2 diabetes mellitus (odds ratio, 0.270; 95% confidence interval, 0.15-0.49; P < .0001). The 3-loci haplotype analysis showed that the A-Ala-T (-1279G/A-Pro12Ala-His478His) haplotype was associated with a reduced risk of type 2 diabetes mellitus (P < .0001). Although our data indicate that the PPARG2 gene variants, independently, have no association with type 2 diabetes mellitus, the 2-loci genotype analysis involving -1279G/A and Pro12Ala variants and the 3-loci haplotype analysis have shown a significant association with type 2 diabetes mellitus in this South Indian population. (C) 2010 Elsevier Inc. All rights reserved.

Genetic Mechanisms Underlying Autism Spectrum Disorders

Autism is a childhood-onset developmental disorder characterized by deficits in reciprocal social interaction, verbal and non-verbal communication, and dependence on routines and rituals. It belongs to a spectrum of disorders (autism spectrum disorders, ASDs) which share core symptoms but show considerable variation in severity. The whole spectrum affects 0.6-0.7% of children worldwide, inducing a substantial public health burden and causing suffering to the affected families. Despite having a very high heritability, ASDs have shown exceptional genetic heterogeneity, which has complicated the identification of risk variants and left the etiology largely unknown. However, recent studies suggest that rare, family-specific factors contribute significantly to the genetic basis of ASDs. In this study, we investigated the role of DISC1 (Disrupted-in-schizophrenia-1) in ASDs, and identified association with markers and haplotypes previously associated with psychiatric phenotypes. We identified four polymorphic micro-RNA target sites in the 3 UTR of DISC1, and showed that hsa-miR-559 regulates DISC1 expression in vitro in an allele-specific manner. We also analyzed an extended autism pedigree with genealogical roots in Central Finland reaching back to the 17th century. To take advantage of the beneficial characteristics of population isolates to gene mapping and reduced genetic heterogeneity observed in distantly related individuals, we performed a microsatellite-based genome-wide screen for linkage and linkage disequilibrium in this pedigree. We identified a putative autism susceptibility locus on chromosome 19p13.3 and obtained further support for previously reported loci at 1q23 and 15q11-q13. To follow-up these findings, we extended our study sample from the same sub-isolate and initiated a genome-wide analysis of homozygosity and allelic sharing using high-density SNP markers. We identified a small number of haplotypes shared by different subsets of the genealogically connected cases, along with convergent biological pathways from SNP and gene expression data, which highlighted axon guidance molecules in the pathogenesis of ASDs. In conclusion, the results obtained in this thesis show that multiple distinct genetic variants are responsible for the ASD phenotype even within single pedigrees from an isolated population. We suggest that targeted resequencing of the shared haplotypes, linkage regions, and other susceptibility loci is essential to identify the causal variants. We also report a possible micro-RNA mediated regulatory mechanism, which might partially explain the wide-range neurobiological effects of the DISC1 gene.

高山松的群体历史及其适应性进化研究

物种形成一直以来是进化生物学的一个中心议题。因为杂交能够快速创造出新的遗传变异和提高基因组重组率，加速物种形成和适应性进化，所以杂交物种形成是其重要的组成部分。但是直到最近几年，人们才逐渐意识到同倍体杂种形成的重要性和普遍性。此外，人工杂交种群结果证实了杂种的适应性起源于亲本基因的超亲分离现象，这些基因之间通过加和效应或上位效应相互连锁，并在杂种后代中逐渐固定，因此，探讨维持这种适应机制的杂种的连锁不平衡模式就显得尤为重要。到目前为止，关于自然杂种的核苷酸多态性和连锁不平衡结构的报道非常稀见，对杂种的进化历史知之更少。高山松能生长于亲本种不能正常繁殖生长的青藏高原上，这种极端环境为杂种的适应性进化和与亲本的生态隔离造就了绝佳的条件，其生态和生殖的稳定性极大地方便了我们研究植物物种形成和适应性进化的遗传机理，对高山松开展的生态学、遗传学的系列分析已使高山松成为同倍体杂种形成的经典范例。 该论文通过随机挑选7个核基因位点对5个高山松群体、3个油松群体和3个云南松群体的164个大配子体单倍体基因组进行了核苷酸多态性分析。研究发现：所有基因座位的单倍型组成和基因谱系的拓扑结构都支持高山松的多次杂交起源。高山松的平均核苷酸多态性与油松相当，θW达到0.0107;比云南松高出一倍;高于已报道的裸子植物类群遗传变异水平。高山松如此高的遗传多态性与其杂交特性和有效群体大有关。我们运用分子钟和共祖模拟分析，推测出高山松的有效群体为7.32 x 105。等位基因共祖时间模拟分析发现，杂交过程早于青藏高原的隆升，也即是说，在杂种稳定成种以前，两个亲本种之间存在广泛的渐渗杂交。此外，高山松杂种的群体间分化严重，并且在不同的群体中，我们找到了多个偏离中性的基因座位。由此说明，高山松复杂的进化历史和适应的策略区域化。 在显著偏离中性的基因座位上，选取群体历史清晰的高山松群体，对其全基因序列多态性和连锁不平衡结构进行了深入调查，以期找到选择作用位点或区域。Ara-like和Dhn1两个基因在高山松与亲本之间都没有发现固定变异。在高山松中，Ara-like基因和Dhn1基因核苷酸多态性式样和连锁不平衡结构存在明显差异，这说明在两个基因上选择作用的程度和方向不尽相同。 我们对Dhn1基因的多态性分布、LD结构和中性偏离水平进行了分析，结果表明，该基因可能在PdNX和PdLZ群体中受到平衡选择的作用。从Dhn1基因谱系结构可以看出有两种来源的等位基因造成这种平衡多态性，一种属于祖先类型，另一种是从云南松继承衍生而来。这两种等位基因之间的分化很大，极有可能造成编码蛋白的亚功能化。事实上，Dhn1编码的脱水蛋白在植物对环境的抗逆过程中发挥着重要作用。由此我们推测高山松在高海拔极端环境中的适应性进化与Dhn1基因的平衡作用有关。 对Ara-like和Dhn1基因进行HKA检测，结果表明，与Dhn1基因相比，在Ara-like基因上，双维管束亚属与单维管束亚属之间存在显著分化，沉默突变位点的分化Ksil达到了0.1392，远远高于平均水平0.0508;在双维管束亚属共发现了43个固定突变位点，其中有6个能导致氨基酸突变，它们有可能导致了Ara-like基因功能的分化。结合裸子植物近缘物种间共享多态性的普遍性，我们推测Ara-like基因可能在单维管束亚属中的进化速率加快，暗示其在单维管束亚属中的适应性进化。 最后，基因内连锁不平衡分析结果显示，随机筛选的基因座位之间不存在连锁不平衡。高山松平均的基因内LD程度非常低，仅有18%的信息位点之间显著连锁。平均基因内LD在油松和高山松中的衰减速率很快，尤其在杂种中下降最快，在不到200bp以内就降到0.1以下。 LD的结果印证了高山松的有效群体大和多次起源特性。另外，我们也怀疑杂交物种形成过程中染色体组的重组和重排频繁发生，也是造成自然杂种现有群体的LD水平低的一个主要原因。 通过LD衰减曲线估计高山松的平均单位重组率,比亲本油松高17倍，比云南松高45倍。这个结果表明，杂种共适应的基因间要维持超亲分离，需要强烈的自然选择压力，才能保持它们之间的连锁不平衡。