鸡单核苷酸多态性与高清晰度QTL图谱的构建

作为一种重要的经济动物和模式动物, 鸡SNP多样性的研究以及鸡主要经济性状QTL定位的研究近年来成绩斐然。文章综述了上述研究成果, 并就SNP标记在鸡QTL精细定位方面的研究前景进行了展望。     

中国人群中STK15基因内两种非同义SNP的连锁不平衡和单体型分析

STK15基因编码一种丝氨酸苏氨酸蛋白激酶,哺乳动物细胞中其过量表达将导致中心体扩增、染色体不稳定和细胞癌变。STK15基因外显子3中有3种非同义单核苷酸多态(SNP),即：91A→T(131F)、169G→A(V571)和311C→T(S104L)。新近研究发现,91A→T与人类肿瘤遗传易感性相关。应用PCR-RFLP技术确定了91A→T(131F)和169G→A(V571)两种SNP在中国人群中的基因型和单体型。采用巢式PCR方法扩增了193例正常个体的DNA样品,通过错配正向巢式内引物引入EcoRⅠ酶切位点。巢式PCR扩增产物用限制性内切酶EcoRⅠ和ACCⅡ双酶切消化,其中EcoRⅠ能酶切91A,AccⅡ能切开169G,用聚丙烯酰胺凝胶电泳银染法鉴定双酶切结果,发现了4种可能的单体型中的3种,其单体型频率分别为：p(91A-169G)=68．65％,p(91T-169A)=10．88％,p(91T-169G)=20．47％,p(91A-169A)=0％;它们组成的6种基因型及频率分别为：91A-169G／91A-169G(46．11％),91A-169G／91T-169A(14．51％),91A-169G／91T-169G(30．57％)．91T-169G／91T-169G(3．11％),91T-169G／91T-169A(4．15％),91T-169A／91T-169A(1．55％)。等位基因及单体型数据分析结果表明,91A→T(131F)和169G→A(V571)之间存在连锁不平衡。     

Linkage disequilibrium fine mapping of quantitative trait loci: A simulation study

Recently, the use of linkage disequilibrium (LD) to locate genes which affect quantitative traits (QTL) has received an increasing interest, but the plausibility of fine mapping using linkage disequilibrium techniques for QTL has not been well studied. The main objectives of this work were to (1) measure the extent and pattern of LD between a putative QTL and nearby markers in finite populations and (2) investigate the usefulness of LD in fine mapping QTL in simulated populations using a dense map of multiallelic or biallelic marker loci. The test of association between a marker and QTL and the power of the test were calculated based on single-marker regression analysis. The results show the presence of substantial linkage disequilibrium with closely linked marker loci after 100 to 200 generations of random mating. Although the power to test the association with a frequent QTL of large effect was satisfactory, the power was low for the QTL with a small effect and/or low frequency. More powerful, multi-locus methods may be required to map low frequent QTL with small genetic effects, as well as combining both linkage and linkage disequilibrium information. The results also showed that multiallelic markers are more useful than biallelic markers to detect linkage disequilibrium and association at an equal distance.     

Allelic association and disease mapping.

The application of allelic association to map genes for complex traits, particularly using high-density maps of single nucleotide polymorphisms in candidate regions, is an area of very active research. Here we present some aspects of the methodology and applications to both major gene mapping, which illustrates the effectiveness of the method, and oligogenes, where methods are still in flux and for which there have been relatively few successes to date. Several important considerations emerge, including the selection of the optimal metric for measuring association and the importance of modelling the decline in association with distance given the variability in association in a candidate region. The Malecot model of association with distance is shown to have a resolution of greater than 50 kilobases but the available evidence suggests that considerably higher resolution might be achieved with dense single nucleotide polymorphism (SNP) maps.     

Ovine alpha-amylase genes: isolation, linkage mapping and association analysis with milk traits

On the basis of comparisons between cattle and sheep genome mapping information the ovine alpha-amylase gene was examined as a possible genetic marker for milk traits in sheep. The objective of the present study was to isolate, map and determine whether this gene is a candidate gene for milk traits. DNA fragments (832 and 2360 bp) corresponding to two different AMY genes were isolated, and one SNP in intron 3 and one GTG deletion in exon 3 of the 2360 bp DNA fragment were found. The 2360 bp ovine AMY DNA fragment was located on chromosome 1 by linkage mapping using the International Mapping Flock. No association was found between estimated breeding values for milk yield, protein and fat contents and AMY genotypes in a daughter design comprising 13 Manchega families with an average of 29 daughters (12-62) per sire.     

Wild GWAS—association mapping in natural populations

The increasing affordability of sequencing and genotyping technologies has transformed the field of molecular ecology in recent decades. By correlating marker variants with trait variation using association analysis, large‐scale genotyping and phenotyping of individuals from wild populations has enabled the identification of genomic regions that contribute to phenotypic differences among individuals. Such “gene mapping” studies are enabling us to better predict evolutionary potential and the ability of populations to adapt to challenges, such as changing environment. These studies are also allowing us to gain insight into the evolutionary processes maintaining variation in natural populations, to better understand genotype‐by‐environment and epistatic interactions and to track the dynamics of allele frequency change at loci contributing to traits under selection. Gene mapping in the wild using genomewide association scans (GWAS) do, however, come with a number of methodological challenges, not least the population structure in space and time inherent to natural populations. We here provide an overview of these challenges, summarize the exciting methodological advances and applications of association mapping in natural populations reported in this special issue and provide some guidelines for future “wild GWAS” research.     

Porcine skeletal muscle differentially expressed gene CMYA1: isolation, characterization, mapping, expression and association analysis with carcass traits

To investigate the differences in gene expression between some obese and lean pig breeds, differential display of mRNA was employed in our previous research. One differentially expressed EST ( BI596262 ) was further identified as the porcine cardiomyopathy associated 1 ( CMYA1 ) gene because of its homology to the human CMYA1 gene. The full-length DNA of the porcine CMYA1 gene encompasses 9379 bp, including a complete open reading frame encoding 1839 amino acid residues, a 158-bp 5'-untranslated region and a 630-bp 3'-untranslated region. The porcine CMYA1 gene was assigned to chromosome 13 by the radiation hybrid panel (IMpRH). The porcine CMYA1 gene was expressed only in the striated muscle. Single nucleotide polymorphism (SNP) scanning in the coding region identified one synonymous mutation (c.1053C>T) and three missense mutations, c.1394A>G (p.His465Arg), c.1751A>G (p.Asp582Gly) and c.3290C>A (p.Thr1097Asp). The allele frequencies were tested among about 200 unrelated pigs from several pig breeds. Linkage mapping was further conducted with the SNP c.1751A>G (p.Asp582Gly) in a Berkshire × Yorkshire resource family and this confirmed that porcine CMYA1 is closely linked with Sw344 (distance  =  2 cM, LOD score is 129.47), an interesting region harbouring a QTL for back fat thickness. Association analysis in our experimental pig population showed that different genotypes of CMYA1 gene were associated with different back fat thicknesses ( P  <   0.05). Our results suggest that the porcine CMYA1 gene has effects on porcine back fat deposition and further investigation will be necessary to illustrate the underlying mechanisms.     

SHEsis, a powerful software platform for analyses of linkage disequilibrium, haplotype construction, and genetic association at polymorphism loci

In multiloci-based genetic association studies of complex diseases, a powerful and high efficient tool for analyses of linkage disequilibrium (LD) between markers, haplotype distributions and many chi-square/p values with a large number of samples has been sought for long. In order to achieve the goal of obtaining meaningful results directly from raw data,we developed a robust and user-friendly software platform with a series of tools for analysis in association study with high efficiency. The platform has been well evaluated by several sets of real data.     

Using multilocus sequence data to assess population structure, natural selection, and linkage disequilibrium in wild tomatoes

We employed a multilocus approach to examine the effects of population subdivision and natural selection on DNA polymorphism in 2 closely related wild tomato species (Solanum peruvianum and Solanum chilense), using sequence data for 8 nuclear loci from populations across much of the species' range. Both species exhibit substantial levels of nucleotide variation. The species-wide level of silent nucleotide diversity is 18% higher in S. peruvianum (pi(sil) approximately 2.50%) than in S. chilense (pi(sil) approximately 2.12%). One of the loci deviates from neutral expectations, showing a clinal pattern of nucleotide diversity and haplotype structure in S. chilense. This geographic pattern of variation is suggestive of an incomplete (ongoing) selective sweep, but neutral explanations cannot be entirely dismissed. Both wild tomato species exhibit moderate levels of population differentiation (average F(ST) approximately 0.20). Interestingly, the pooled samples (across different demes) exhibit more negative Tajima's D and Fu and Li's D values; this marked excess of low-frequency polymorphism can only be explained by population (or range) expansion and is unlikely to be due to population structure per se. We thus propose that population structure and population/range expansion are among the most important evolutionary forces shaping patterns of nucleotide diversity within and among demes in these wild tomatoes. Patterns of population differentiation may also be impacted by soil seed banks and historical associations mediated by climatic cycles. Intragenic linkage disequilibrium (LD) decays very rapidly with physical distance, suggesting high recombination rates and effective population sizes in both species. The rapid decline of LD seems very promising for future association studies with the purpose of mapping functional variation in wild tomatoes.     

The signature of positive selection on standing genetic variation

Considerable interest is focused on the use of polymorphism data to identify regions of the genome that underlie recent adaptations. These searches are guided by a simple model of positive selection, in which a mutation is favored as soon as it arises. This assumption may not be realistic, as environmental changes and range expansions may lead previously neutral or deleterious alleles to become beneficial. We examine what effect this mode of selection has on patterns of variation at linked neutral sites by implementing a new coalescent model of positive directional selection on standing variation. In this model, a neutral allele arises and drifts in the population, then at frequency f becomes beneficial, and eventually reaches fixation. Depending on the value of f, this scenario can lead to a large variance in allele frequency spectra and in levels of linkage disequilibrium at linked, neutral sites. In particular, for intermediate f, the beneficial substitution often leads to a loss of rare alleles--a pattern that differs markedly from the signature of directional selection currently relied on by researchers. These findings highlight the importance of an accurate characterization of the effects of positive selection, if we are to reliably identify recent adaptations from polymorphism data.     

Identification of X-linked quantitative trait loci affecting cold tolerance in Drosophila melanogaster and fine mapping by selective sweep analysis

Drosophila melanogaster is a cosmopolitan species that colonizes a great variety of environments. One trait that shows abundant evidence for naturally segregating genetic variance in different populations of D. melanogaster is cold tolerance. Previous work has found quantitative trait loci (QTL) exclusively on the second and the third chromosomes. To gain insight into the genetic architecture of cold tolerance on the X chromosome and to compare the results with our analyses of selective sweeps, a mapping population was derived from a cross between substitution lines that solely differed in the origin of their X chromosome: one originates from a European inbred line and the other one from an African inbred line. We found a total of six QTL for cold tolerance factors on the X chromosome of D. melanogaster. Although the composite interval mapping revealed slightly different QTL profiles between sexes, a coherent model suggests that most QTL overlapped between sexes, and each explained around 5–14% of the genetic variance (which may be slightly overestimated). The allelic effects were largely additive, but we also detected two significant interactions. Taken together, this provides evidence for multiple QTL that are spread along the entire X chromosome and whose effects range from low to intermediate. One detected transgressive QTL influences cold tolerance in different ways for the two sexes. While females benefit from the European allele increasing their cold tolerance, males tend to do better with the African allele. Finally, using selective sweep mapping, the candidate gene CG16700 for cold tolerance colocalizing with a QTL was identified.     

Molecular diversity and association mapping of fiber quality traits in exotic G. hirsutum L. germplasm

The narrow genetic base of cultivated cotton germplasm is hindering the cotton productivity worldwide. Although potential genetic diversity exists in Gossypium genus, it is largely ‘underutilized’ due to photoperiodism and the lack of innovative tools to overcome such challenges. The application of linkage disequilibrium (LD)-based association mapping is an alternative powerful molecular tool to dissect and exploit the natural genetic diversity conserved within cotton germplasm collections, greatly accelerating still ‘lagging’ cotton marker-assisted selection (MAS) programs. However, the extent of genome-wide linkage disequilibrium (LD) has not been determined in cotton. We report the extent of genome-wide LD and association mapping of fiber quality traits by using a 95 core set of microsatellite markers in a total of 285 exotic Gossypium hirsutum accessions, comprising of 208 landrace stocks and 77 photoperiodic variety accessions. We demonstrated the existence of useful genetic diversity within exotic cotton germplasm. In this germplasm set, 11–12% of SSR loci pairs revealed a significant LD. At the significance threshold (r² ≥ 0.1), a genome-wide average of LD declines within the genetic distance at < 10 cM in the landrace stocks germplasm and > 30 cM in variety germplasm. Genome wide LD at r² ≥ 0.2 was reduced on average to  1–2 cM in the landrace stock germplasm and 6–8 cM in variety germplasm, providing evidence of the potential for association mapping of agronomically important traits in cotton. We observed significant population structure and relatedness in assayed germplasm. Consequently, the application of the mixed liner model (MLM), considering both kinship (K) and population structure (Q) detected between 6% and 13% of SSR markers associated with the main fiber quality traits in cotton. Our results highlight for the first time the feasibility and potential of association mapping, with consideration of the population structure and stratification existing in cotton germplasm resources. The number of SSR markers associated with fiber quality traits in diverse cotton germplasm, which broadly covered many historical meiotic events, should be useful to effectively exploit potentially new genetic variation by using MAS programs.     

关联分析及其在真菌遗传学研究中的应用

近年来,基于连锁不平衡的关联分析被逐渐应用于解析真菌数量性状的遗传基础。本文在介绍关联分析方法的基础上,综述了关联分析在真菌遗传学中的研究进展,并对其未来的应用前景进行了展望。     

A unified DNA sequence and non‐DNA sequence mapping model of complex traits

Increasing evidence shows that quantitative inheritance is based on both DNA sequence and non‐DNA sequence variants. However, how to simultaneously detect these variants from a mapping study has been unexplored, hampering our effort to illustrate the detailed genetic architecture of complex traits. We address this issue by developing a unified model of quantitative trait locus (QTL) mapping based on an open‐pollinated design composed of randomly sampling maternal plants from a natural population and their half‐sib seeds. This design forms a two‐level hierarchical platform for a joint linkage‐linkage disequilibrium analysis of population structure. The EM algorithm was implemented to estimate and test DNA sequence‐based effects and non‐DNA sequence‐based effects of QTLs. We applied this model to analyze genetic mapping data from the OP design of a gymnosperm coniferous species, Torreya grandis, identifying 25 significant DNA sequence and non‐DNA sequence QTLs for seedling height and diameter growth in different years. Results from computer simulation show that the unified model has good statistical properties and is powerful for QTL detection. Our model enables the tests of how a complex trait is affected differently by DNA‐based effects and non‐DNA sequence‐based transgenerational effects, thus allowing a more comprehensive picture of genetic architecture to be charted and quantified.     

A general method for controlling the genome-wide type I error rate in linkage and association mapping experiments in plants

Control of the genome-wide type I error rate (GWER) is an important issue in association mapping and linkage mapping experiments. For the latter, different approaches, such as permutation procedures or Bonferroni correction, were proposed. The permutation test, however, cannot account for population structure present in most association mapping populations. This can lead to false positive associations. The Bonferroni correction is applicable, but usually on the conservative side, because correlation of tests cannot be exploited. Therefore, a new approach is proposed, which controls the genome-wide error rate, while accounting for population structure. This approach is based on a simulation procedure that is equally applicable in a linkage and an association-mapping context. Using the parameter settings of three real data sets, it is shown that the procedure provides control of the GWER and the generalized genome-wide type I error rate (GWER(k)).     

Detecting signatures of selection in nine distinct lines of broiler chickens

Modern commercial chickens have been bred for one of two specific purposes: meat production (broilers) or egg production (layers). This has led to large phenotypic changes, so that the genomic signatures of selection may be detectable using statistical techniques. Genetic differentiation between nine distinct broiler lines was calculated using Weir and Cockerham's pairwise F_ST estimator for 11 003 genome‐wide markers to identify regions showing evidence of differential selection across lines. Differentiation measures were averaged into overlapping sliding windows for each line, and a permutation approach was used to determine the significance of each window. A total of 51 regions were found to show significant differentiation between the lines. Several lines were consistently found to share significant regions, suggesting that the pattern of line divergence is related to selection for broiler traits. The majority of the 51 regions contain QTL relating to broiler traits, but only five of them were found to be significantly enriched for broiler QTL, including a region on chromosome 27 containing 39 broiler QTL and 114 genes. Additionally, a number of these regions have been identified by other selection mapping studies. This study has identified a large number of potential selection signatures, and further tests with higher‐density marker data may narrow these regions down to individual genes.     

Recombination, balancing selection and adaptive evolution in the aflatoxin gene cluster of Aspergillus parasiticus

Aflatoxins are toxic and carcinogenic polyketides produced by several Aspergillus species that are known to contaminate agricultural commodities, posing a serious threat to animal and human health. Aflatoxin (AF) biosynthesis is almost fully characterized and involves the coordinated expression of approximately 25 genes clustered in a 70-kb DNA region. Aspergillus parasiticus is an economically important and common agent of AF contamination. Naturally occurring nonaflatoxigenic strains of A. parasiticus are rarely found and generally produce O-methylsterigmatocystin (OMST), the immediate precursor of AF. To elucidate the evolutionary forces acting to retain AF and OMST pathway extrolites (chemotypes), we sequenced 21 intergenic regions spanning the entire cluster in 24 A. parasiticus isolates chosen to represent the genetic diversity within a single Georgia field population. Linkage disequilibrium analyses revealed five distinct recombination blocks in the A. parasiticus cluster. Phylogenetic network analyses showed a history of recombination between chemotype-specific haplotypes, as well as evidence of contemporary recombination. We performed coalescent simulations of variation in recombination blocks and found an approximately twofold deeper coalescence for cluster genealogies compared to noncluster genealogies, our internal standard of neutral evolution. Significantly deeper cluster genealogies are indicative of balancing selection in the AF cluster of A. parasiticus and are further corroborated by the existence of trans-species polymorphisms and common haplotypes in the cluster for several closely related species. Estimates of Ka/Ks for representative cluster genes provide evidence of selection for OMST and AF chemotypes, and indicate a possible role of chemotypes in ecological adaptation and speciation.     

Haplotype structure and copy number polymorphism of the beta‐defensin 7 genes in diverse chicken breeds

Beta‐defensins is a family of avian peptides related to the innate immune system. Copy number variation was recently reported for the avian beta‐defensin 7 gene (AvBD7) between the highly inbred Leghorn and Fayoumi lines. Here, we examined copy number variants in 35 different chicken breeds and found that 31 of them have at least the same representation of the duplicated AvBD7 allele. We also found haplotypes upstream of the AvBD6 regions that are strongly linked to the AvBD7 duplication. We observed a strong linkage disequilibrium spanning of the upstream region of the AvBD6 gene, with two SNPs being flanking markers to detect duplication of the AvBD7.     

SNP-based association mapping of Arachnomelia in Fleckvieh cattle

Bovine arachnomelia is an inherited congenital disorder with malformation mainly of the limbs, the vertebral column and the skull, following a monogenic autosomal recessive heredity. Despite almost identical pathological findings, arachnomelia has previously been mapped to bovine chromosome 23 and 5 in Fleckvieh and Braunvieh respectively. Therefore, this disorder may be an example of locus heterogeneity in cattle. This study aimed to refine the candidate region to allow positional cloning and sequence analyses of candidate genes in Fleckvieh cattle. For that purpose, a case-control association mapping design was set up with a case group of 16 pre-selected affected individuals and a control group of 50 unrelated animals. The subset of affected animals was selected from a total of 129 pathologically confirmed cases due to the occurrence of recombination(s) within a 14.5 cM candidate interval previously mapped to chromosome 23. Six linked microsatellites currently used for indirect gene testing in Fleckvieh were analysed for this purpose. In all selected cases, a genome-wide scan using 44 473 informative SNPs revealed shared segments of homozygosity at 15 adjacent SNPs on chromosome 23. Additional haplotype analysis of 37 carrier bulls confirmed the localization of the arachnomelia locus to a region of 927 kb (13.622-14.549 Mb) containing molybdenum cofactor biosynthesis protein 1 gene, the most likely candidate gene for arachnomelia in Fleckvieh. The number of recombinant haplotypes observed in cases was more than doubled compared with the number of expected recombinations. This remarkably increased mapping resolution and thus illustrates the benefit of pre-selection in association studies.     

Genetic dissection of drought and heat‐responsive agronomic traits in wheat

High yield and wide adaptation are principal targets of wheat breeding but are hindered by limited knowledge on genetic basis of agronomic traits and abiotic stress tolerances. In this study, 277 wheat accessions were phenotyped across 30 environments with non‐stress, drought‐stressed, heat‐stressed, and drought‐heat‐stressed treatments and were subjected to genome‐wide association study using 395 681 single nucleotide polymorphisms. We detected 295 associated loci including consistent loci for agronomic traits across different treatments and eurytopic loci for multiple abiotic stress tolerances. A total of 22 loci overlapped with quantitative trait loci identified by biparental quantitative trait loci mapping. Six loci were simultaneously associated with agronomic traits and abiotic stress tolerance, four of which fell within selective sweep regions. Selection in Chinese wheat has increased the frequency of superior marker alleles controlling yield‐related traits in the four loci during past decades, which conversely diminished favourable genetic variation controlling abiotic stress tolerance in the same loci; two promising candidate paralogous genes colocalized with such loci, thereby providing potential targets for studying the molecular mechanism of stress tolerance–productivity trade‐off. These results uncovering promising alleles controlling agronomic traits and/or multiple abiotic stress tolerances, providing insights into heritable covariation between yield and abiotic stress tolerance, will accelerate future efforts for wheat improvement.