Comparison of single-nucleotide polymorphisms and microsatellites in detecting quantitative trait loci for alcoholism: The Collaborative Study on the Genetics of Alcoholism

Background

The feasibility of effectively analyzing high-density single nucleotide polymorphism (SNP) maps in whole genome scans of complex traits is not known. The purpose of this study was to compare variance components linkage results using different density marker maps in data from the Collaborative Study on the Genetics of Alcoholism (COGA). Marker maps having an average spacing of 10 cM (microsatellite), 0.78 cM (SNP1), and 0.31 cM (SNP2) were used to identify quantitative trait loci (QTLs) affecting maximum number of alcoholic drinks consumed in a 24-hour period (lnmaxalc).

Results

Heritability of lnmaxalc was estimated to be 15%. Multipoint variance components linkage analysis revealed similar linkage patterns among the three marker panels, with the SNP maps consistently yielding higher LOD scores. Robust LOD scores > 1.0 were observed on chromosomes 1 and 13 for all three marker maps. Additional LODs > 1.0 were observed on chromosome 4 with both SNP maps and on chromosomes 18 and 21 with the SNP2 map. Peak LOD scores for lnmaxalc were observed on chromosome 1, although none reached genome-wide statistical significance. Quantile-quantile plots revealed that the multipoint distribution of SNP results appeared to fit the asymptotic null distribution better than the twopoint results.

Conclusion

In conclusion, variance-components linkage analysis using high-density SNP maps provided higher LOD scores compared with the standard microsatellite map, similar to studies using nonparametric linkage methods. Widespread application of SNP maps will depend on further improvements in the computational methods implemented in current software packages.

     

Linkage analysis of complex diseases using microsatellites and single-nucleotide polymorphisms: application to alcoholism

The efficacy of linkage studies using microsatellites and single-nucleotide polymorphisms (SNPs) was evaluated. Analyzed data were supplied by the Collaborative Study on the Genetics of Alcoholism (COGA). Alcoholism was analyzed together with a simulated trait caused by a gene of known position, through a nonparametric linkage test (NPL). For the alcoholism trait, four densities of SNPs (1 SNP per 0.2 cM, 0.5 cM, 1 cM and 2 cM) showed higher peaks of NPL z scores and smaller significant p-values than the usual 10-cM density of microsatellites. However, the two highest densities of SNPs had unstable z score signals, and therefore were difficult to interpret. Analyzing a simulated trait with the same markers in the same pedigrees, we confirmed the higher power of all four densities of SNPs compared to the 10-cM microsatellites panel, although the existence of other confounding peaks was confirmed for maps that are denser than 1 SNP/cM. We further showed that estimating the gene position using SNPs is far less biased than using the usual panel of microsatellites (biases of 0-2 cM for SNPs vs. 8.9 cM for microsatellites). We conclude that using dense maps of SNPs in linkage analysis is more powerful and less biased than using the 10-cM maps of microsatellites. However, linkage signals can be unstable and difficult to interpret when several SNPs are genotyped per centimorgan. The power and accuracy of 1 SNP/cM or 1 SNP/2 cM may be sufficient in a genome-wide linkage scan while denser maps may be most useful in fine-gene mapping studies exploiting linkage disequilibrium.     

Comparing single-nucleotide polymorphism marker-based and microsatellite marker-based linkage analyses

We compared linkage analysis results for an alcoholism trait, ALDX1 (DSM-III-R and Feigner criteria) using a nonparametric linkage analysis method, which takes into account allele sharing among several affected persons, for both microsatellite and single-nucleotide polymorphism (SNP) markers (Affymetrix and Illumina) in the Collaborative Study on the Genetics of Alcoholism (COGA) dataset provided to participants at the Genetic Analysis Workshop 14 (GAW14). The two sets of linkage results from the dense Affymetrix SNP markers and less densely spaced Illumina SNP markers are very similar. The linkage analysis results from microsatellite and SNP markers are generally similar, but the match is not perfect. Strong linkage peaks were found on chromosome 7 in three sets of linkage analyses using both SNP and microsatellite marker data. We also observed that for SNP markers, using the given genetic map and using the map by converting 1 megabase pair (1 Mb) to 1 centimorgan (cM), did not change the linkage results. We recommend the use of the 1 Mb-to-1 cM converted map in a first round of linkage analysis with SNP markers in which map integration is an issue.     

The effect of linkage disequilibrium on linkage analysis of incomplete pedigrees

Dense SNP maps can be highly informative for linkage studies. But when parental genotypes are missing, multipoint linkage scores can be inflated in regions with substantial marker-marker linkage disequilibrium (LD). Such regions were observed in the Affymetrix SNP genotypes for the Genetic Analysis Workshop 14 (GAW14) Collaborative Study on the Genetics of Alcoholism (COGA) dataset, providing an opportunity to test a novel simulation strategy for studying this problem. First, an inheritance vector (with or without linkage present) is simulated for each replicate, i.e., locations of recombinations and transmission of parental chromosomes are determined for each meiosis. Then, two sets of founder haplotypes are superimposed onto the inheritance vector: one set that is inferred from the actual data and which contains the pattern of LD; and one set created by randomly selecting parental alleles based on the known allele frequencies, with no correlation (LD) between markers. Applying this strategy to a map of 176 SNPs (66 Mb of chromosome 7) for 100 replicates of 116 sibling pairs, significant inflation of multipoint linkage scores was observed in regions of high LD when parental genotypes were set to missing, with no linkage present. Similar inflation was observed in analyses of the COGA data for these affected sib pairs with parental genotypes set to missing, but not after reducing the marker map until r2 between any pair of markers was 相似文献   

Levels of linkage disequilibrium in a wild bird population

Population-based mapping approaches are attractive for tracing the genetic background to phenotypic traits in wild species, given that it is often difficult to gather extensive and well-defined pedigrees needed for quantitative trait locus analysis. However, the feasibility of association or hitch-hiking mapping is dependent on the degree of linkage disequilibrium (LD) in the population, on which there is yet limited information for wild species. Here we use single nucleotide polymorphism (SNP) markers from 23 genes in a recently established linkage map of the Z chromosome of the collared flycatcher, to study the extent of LD in a natural bird population. In most but not all cases we find SNPs within the same intron (less than 500 bp) to be in perfect LD. However, LD then decays to background level at a distance 1cM or 400-500 kb. Although LD seems more extensive than in other species, if the observed pattern is representative for other regions of the genome and turns out to be a general feature of natural bird populations, dense marker maps might be needed for genome scans aimed at identifying association between marker and trait loci.     

Comparison of marker types and map assumptions using Markov chain Monte Carlo-based linkage analysis of COGA data

We performed multipoint linkage analysis of the electrophysiological trait ECB21 on chromosome 4 in the full pedigrees provided by the Collaborative Study on the Genetics of Alcoholism (COGA). Three Markov chain Monte Carlo (MCMC)-based approaches were applied to the provided and re-estimated genetic maps and to five different marker panels consisting of microsatellite (STRP) and/or SNP markers at various densities. We found evidence of linkage near the GABRB1 STRP using all methods, maps, and marker panels. Difficulties encountered with SNP panels included convergence problems and demanding computations.     

Linkage analysis of the GAW14 simulated dataset with microsatellite and single-nucleotide polymorphism markers in large pedigrees

Recent studies have suggested that a high-density single nucleotide polymorphism (SNP) marker set could provide equivalent or even superior information compared with currently used microsatellite (STR) marker sets for gene mapping by linkage. The focus of this study was to compare results obtained from linkage analyses involving extended pedigrees with STR and single-nucleotide polymorphism (SNP) marker sets. We also wanted to compare the performance of current linkage programs in the presence of high marker density and extended pedigree structures. One replicate of the Genetic Analysis Workshop 14 (GAW14) simulated extended pedigrees (n = 50) from New York City was analyzed to identify the major gene D2. Four marker sets with varying information content and density on chromosome 3 (STR [7.5 cM]; SNP [3 cM, 1 cM, 0.3 cM]) were analyzed to detect two traits, the original affection status, and a redefined trait more closely correlated with D2. Multipoint parametric and nonparametric linkage analyses (NPL) were performed using programs GENEHUNTER, MERLIN, SIMWALK2, and S.A.G.E. SIBPAL. Our results suggested that the densest SNP map (0.3 cM) had the greatest power to detect linkage for the original trait (genetic heterogeneity), with the highest LOD score/NPL score and mapping precision. However, no significant improvement in linkage signals was observed with the densest SNP map compared with STR or SNP-1 cM maps for the redefined affection status (genetic homogeneity), possibly due to the extremely high information contents for all maps. Finally, our results suggested that each linkage program had limitations in handling the large, complex pedigrees as well as a high-density SNP marker set.     

Construction of a genetic linkage map in Fenneropenaeus chinensis using SNP markers

Genetic linkage maps of Fenneropenaeus chinensis were constructed using a “double pseudo-testcross” strategy with 200 single nucleotide polymorphisms (SNPs) markers. This study represents the first SNP genetic linkage map for F. chinensis. The parents and F ₁ progeny of 100 individuals were used as mapping populations. 21 genetic linkage groups in the male and female maps were identified. The male linkage map was composed of 115 loci and spanned 879.7 cM, with an average intermarker spacing of 9.4 cM, while the female map was composed of 119 loci and spanned 876.2 cM, with an average intermarker spacing of 8.9 cM. The estimated coverage of the linkage maps was 51.94% for the male and 53.77% for the female, based on two estimates of genome length. The integrated map contains 180 markers distributed in 16 linkage groups, and spans 899.3 cM with an average marker interval of 5.2 cM. This SNP genetic map lays the foundation for future shrimp genomics and genetic breeding studies, especially the discovery of gene or regions for economically important traits in Chinese shrimp.     

An apricot (Prunus armeniaca L.) F2 progeny linkage map based on SSR and AFLP markers,mapping plum pox virus resistance and self-incompatibility traits

A genetic linkage map of apricot ( Prunus armeniaca L.) was constructed using AFLP and SSR markers. The map is based on an F(2) population (76 individuals) derived from self-pollination of an F(1) individual ('Lito') originated from a cross between 'Stark Early Orange' and 'Tyrinthos'. This family, designated as 'Lito' x 'Lito', segregated for two important agronomical traits: plum pox virus resistance (PPV) and self-incompatibility. A total of 211 markers (180 AFLPs, 29 SSRs and two agronomic traits) were assigned to 11 linkage groups covering 602 cM of the apricot genome. The average distance (cM/marker) between adjacent markers is 3.84 cM. The PPV resistance trait was mapped on linkage group G1 and the self-incompatibility trait was mapped on linkage group G6. Twenty two loci held in common with other Prunus maps allowed us to compare and establish homologies among the respective linkage groups.     

Genetic mapping of local adaptation along the altitudinal gradient in <Emphasis Type="Italic">Abies sachalinensis</Emphasis>

Understanding the genetic bases of local adaptation in dominant conifer species is critical in predicting the impacts of rapid climate change on forest ecosystems. However, the genetic basis of adaptation is not yet fully understood due to the huge and complex genomes of conifers and the unavailability to date of suitable crossing material. In this study, we constructed a linkage map for Abies sachalinensis (2n = 24) and investigated quantitative trait loci (QTLs) associated with local adaptation along an altitudinal gradient. A segregating population of 239 seedlings was produced from a cross between two F₁ hybrids (high-altitude × low-altitude genotypes). QTL mapping of phenological and growth traits was performed using a pseudo-testcross strategy with linkage maps based on 1251 single-nucleotide polymorphism (SNP) and three simple sequence repeat (SSR) markers. Two maps consisting of 12 linkage groups with an average marker interval of ca. 3 cM were constructed for each parent. The total lengths of the maps were 1861 and 1949 cM. A permutation test identified four significant QTLs and 11 additional suggestive QTLs, with high logarithm of odds (LOD) scores (> 3.0). This is the first highly saturated linkage map produced for Abies taxa. Our results suggest that spring bud phenology is controlled by several QTLs with moderate effects. The use of the mapping population created by crossing two hybrids (high × low altitude genotypes) and numerous SNP markers enabled us to investigate the genetic basis of adaptive traits in conifer species.     

Linkage disequilibrium across two different single-nucleotide polymorphism genome scans

Linkage disequilibrium (LD) content was calculated for the Genetic Analysis Workshop 14 Affymetrix and Illumina single-nucleotide polymorphism (SNP) genome scans of the Collaborative Study on the Genetics of Alcoholism samples. Pair-wise LD was measured as both D' and r2 on 505 pedigree founder individuals. The r2 estimates were then used to correct the multipoint identity by descent matrix (MIBD) calculation to account for LD and LOD scores on chromosomes 3 and 18 were calculated for COGA's ttdt3 electrophysiological trait using those MIBDs. Extensive LD was observed throughout both marker sets, and it was higher in Affymetrix's more dense SNP map. However, SNP density did not solely account for Affymetrix's higher LD. MIBD estimation procedures assume linkage equilibrium to construct genotypes of non-genotyped pedigree founder individuals, and dense SNP genotyping maps are likely to contain moderate to high LD between markers. LOD score plots calculated after correction for LD followed the same general pattern as uncorrected ones. Since in our study almost half of the pedigree founders were genotyped, it is possible that LD had a minor impact on the LOD scores. Caution should probably be taken when using high density SNP maps when many non-genotyped founders are present in the study pedigrees.     

Construction of High-Density Linkage Maps of Populus deltoides × P. simonii Using Restriction-Site Associated DNA Sequencing

Although numerous linkage maps have been constructed in the genus Populus, they are typically sparse and thus have limited applications due to low throughput of traditional molecular markers. Restriction-site associated DNA sequencing (RADSeq) technology allows us to identify a large number of single nucleotide polymorphisms (SNP) across genomes of many individuals in a fast and cost-effective way, and makes it possible to construct high-density genetic linkage maps. We performed RADSeq for 299 progeny and their two parents in an F₁ hybrid population generated by crossing the female Populus deltoides ‘I-69’ and male Populus simonii ‘L3’. A total of 2,545 high quality SNP markers were obtained and two parent-specific linkage maps were constructed. The female genetic map contained 1601 SNPs and 20 linkage groups, spanning 4,249.12 cM of the genome with an average distance of 2.69 cM between adjacent markers, while the male map consisted of 940 SNPs and also 20 linkage groups with a total length of 3,816.24 cM and an average marker interval distance of 4.15 cM. Finally, our analysis revealed that synteny and collinearity are highly conserved between the parental linkage maps and the reference genome of P. trichocarpa. We demonstrated that RAD sequencing is a powerful technique capable of rapidly generating a large number of SNPs for constructing genetic maps in outbred forest trees. The high-quality linkage maps constructed here provided reliable genetic resources to facilitate locating quantitative trait loci (QTLs) that control growth and wood quality traits in the hybrid population.     

Genetic mapping of a new set of microsatellite markers in a reference common bean (Phaseolus vulgaris) population BAT93 x Jalo EEP558

The present study describes a new set of 61 polymorphic microsatellite markers for beans and the construction of a genetic map using the BAT93 x Jalo EEP558 (BJ) population for the purpose of developing a reference linkage map for common bean (Phaseolus vulgaris). The main objectives were to integrate new microsatellites on the existing framework map of the BJ population, and to develop the first linkage map for the BJ population based exclusively on microsatellites. Of the total of 264 microsatellites evaluated for polymorphism, 42.8% showed polymorphism between the genitors. An integrated map was created totaling 199 mapped markers in 13 linkage groups, with an observed length of 1358 cM and a mean distance between markers of 7.23 cM. For the map constructed exclusively with microsatellites, 106 markers were placed in 12 groups with a total length of 606.8 cM and average distance of 6.8 cM. Linkage group designation and marker order for BM microsatellites generally agreed with previous mapping, while the new microsatellites were well distributed across the genome, corroborating the utility of the BJ population for a reference map. The extensive use of the microsatellites and the availability of a reference map can help in the development of other genetic maps for common bean through the transfer of information of marker order and linkage, which will allow comparative analysis and map integration, especially for future quantitative trait loci and association mapping studies.     

Construction of a SNP-based genetic linkage map for kiwifruit using next-generation restriction-site-associated DNA sequencing (RADseq)

Kiwifruit is a perennial horticultural crop species of the Actinidiaceae family and has high nutritional value. For a species with a long generation time, traditional breeding and genetic improvement is predicted to take more than 20 years to obtain superior cultivars. Thus, marker-assisted selection (MAS) should be used to accelerate the breeding process. Development of a genetic linkage map and molecular markers are pre-requisites for MAS of crop species. Here, we report a genome-wide SNP-based genetic map of kiwifruit by analysing next-generation restriction-site-associated DNA sequencing (RADseq) reads. To construct a genetic linkage map, a 102 F1 line mapping population of Actinidia chinensis (2n = 58) was derived by combining parents that had contrasting phenotypic traits. The maternal map contained 4112 SNP loci and spanned a distance of 3821 cM, with an average adjacent-marker interval length of 0.929 cM. The map length of the 29 linkage groups ranged from 78.3 to 169.9 cM, with an average length of 131.8 cM. High levels of collinearity between the 29 genetic maps with the kiwifruit reference genome were found. The genetic map developed in this study can serve as an important platform to improve kiwifruit research, including anchoring unmapped scaffolds of the kiwifruit genome sequence and mapping QTLs (quantitative trait loci) that control economically important traits.     

A genome-wide linkage analysis of alcoholism on microsatellite and single-nucleotide polymorphism data,using alcohol dependence phenotypes and electroencephalogram measures

The Collaborative Study on the Genetics of Alcoholism (COGA) is a large-scale family study designed to identify genes that affect the risk for alcoholism and alcohol-related phenotypes. We performed genome-wide linkage analyses on the COGA data made available to participants in the Genetic Analysis Workshop 14 (GAW 14). The dataset comprised 1,350 participants from 143 families. The samples were analyzed on three technologies: microsatellites spaced at 10 cM, Affymetrix GeneChip Human Mapping 10 K Array (HMA10K) and Illumina SNP-based Linkage III Panel. We used ALDX1 and ALDX2, the COGA definitions of alcohol dependence, as well as electrophysiological measures TTTH1 and ECB21 to detect alcoholism susceptibility loci. Many chromosomal regions were found to be significant for each of the phenotypes at a p-value of 0.05. The most significant region for ALDX1 is on chromosome 7, with a maximum LOD score of 2.25 for Affymetrix SNPs, 1.97 for Illumina SNPs, and 1.72 for microsatellites. The same regions on chromosome 7 (96-106 cM) and 10 (149-176 cM) were found to be significant for both ALDX1 and ALDX2. A region on chromosome 7 (112-153 cM) and a region on chromosome 6 (169-185 cM) were identified as the most significant regions for TTTH1 and ECB21, respectively. We also performed linkage analysis on denser maps of markers by combining the SNPs datasets from Affymetrix and Illumina. Adding the microsatellite data to the combined SNP dataset improved the results only marginally. The results indicated that SNPs outperform microsatellites with the densest marker sets performing the best.     

RAPD和SSR两种标记构建的中国对虾遗传连锁图谱

利用RAPD和SSR分子标记结合拟测交策略,对中国对虾(Fenneropenaeuschinensis)“黄海1号”雌虾与野生雄虾作为亲本进行单对杂交产生的F1代,采用RAPD和SSR两种分子标记技术初步构建了中国对虾雌、雄遗传连锁图谱。对460个RAPD引物和44对SSR引物进行筛选,共选出61个RAPD引物和20对SSR引物,用于对父母本和82个F1个体进行遗传分析。共得到母本分离标记146个(RAPD标记128个,微卫星标记18个)和父本分离标记127个(RAPD标记109个,微卫星标记18个)。雌性图谱包括8个连锁群、9个三联体和14个连锁对,标记间平均间隔为11·28cM,图谱共覆盖1173cM,覆盖率为59·36%;雄性图谱包括10个连锁群、12个三联体和7个连锁对,标记间平均间隔为12·05cM,图谱共覆盖1144·6cM,覆盖率为62·01%。中国对虾遗传图谱的构建为其分子标记辅助育种、比较基因组作图及数量性状位点的定位与克隆奠定了基础。     

The extent and distribution of linkage disequilibrium in a multi-hierarchic outbred canine pedigree

A canine integrated linkage-radiation map has been recently constructed by using microsatellite markers. This map, with a good coverage of the canine genome, allows for a genome-wide search for the extent and distribution of linkage disequilibrium derived from linkage and evolutionary forces. In this study, we genotyped an outbred pedigree between Labrador retriever and Greyhound breeds with a set of microsatellite markers (240) from the canine linkage map. Linkage disequilibrium was measured between all syntenic and nonsyntenic marker pairs. Analysis of syntenic pairs revealed a significant correlation (–0.229, P < 0.001) between linkage disequilibrium and genetic distance (log transformed). Significant linkage disequilibria were observed more frequently between syntenic pairs spaced <40 cM than those paced >40 cM. There is a clear trend for linkage disequilibrium to decline with marker distance. From our results, a genome-wide screen with markers at low to moderate density (1–2 per 10 cM) should take full advantage of linkage disequilibrium for quantitative trait locus mapping in dogs. This study supports the appropriateness of linkage disequilibrium analysis to detect and map quantitative trait loci underlying complex traits in dogs.     

Genetic analysis of floral anthocyanin pigmentation traits in Asiatic hybrid lily using molecular linkage maps

To understand the genetic background of two floral anthocyanin pigmentation traits, anthocyanin pigmentation in the flower tepals and spot formation, in the Asiatic hybrid lily (2n = 24), segregation of the two traits among 96 F₁ plants derived from a cross between commercial cultivars 'Montreux' and 'Connecticut King' were investigated. 'Montreux' has anthocyanin pigmentation in the tepals with many spots, and 'Connecticut King' has flowers with carotenoid pigmentation without spots. The F₁ plants with or without anthocyanin pigment in the tepals segregated with a 1:1 segregation ratio, indicating that a single gene controls anthocyanin pigmentation in the tepals. The number of spots per square centimeter of all tepals showed continuous distribution in the F₁ plants. To map the loci for the two anthocyanin pigmentation traits, molecular linkage maps in the Asiatic hybrid lily were constructed using a double pseudo-testcross strategy, with the same F₁ plants used for phenotypic evaluation, and 212 PCR-based DNA markers. The trait for anthocyanin pigmentation in tepals was used as a trait marker. The map of 'Montreux' comprised 95 markers in 26 linkage groups, and the map of 'Connecticut King' used 119 markers in 24 linkage groups. The total map lengths were 867.5 and 1,114.8 cM, respectively. The trait locus for anthocyanin pigmentation in the tepals was between markers ASR35-180 and P506-40 in linkage group 1 of the 'Montreux' map with a map distance of 1.2 cM and 2.6 cM, respectively. A single-point analysis of quantitative trait loci (QTLs) for tepal spot number identified two putative QTLs in linkage groups 1 and 19 of the 'Connecticut King' map. One putative QTL in linkage group 19 explained 64% of the total phenotypic variation. Because both putative QTLs were mapped on the linkage map of 'Connecticut King' that has no spots, dominant alleles of them might suppress spot formation.     

A genetic linkage map based on AFLP and SSR markers and mapping of QTL for dry-matter content in sweetpotato

We developed a genetic linkage map of sweetpotato using amplified fragment length polymorphism (AFLP) and simple sequence repeat (SSR) markers and a mapping population consisting of 202 individuals derived from a broad cross between Xushu 18 and Xu 781, and mapped quantitative trait loci (QTL) for the storage root dry-matter content. The linkage map for Xushu 18 included 90 linkage groups with 2077 markers (1936 AFLP and 141 SSR) and covered 8,184.5 cM with an average marker distance of 3.9 cM, and the map for Xu 781 contained 90 linkage groups with 1954 markers (1824 AFLP and 130 SSR) and covered 8,151.7 cM with an average marker distance of 4.2 cM. The maps described herein have the best coverage of the sweetpotato genome and the highest marker density reported to date. These are the first maps developed that have 90 complete linkage groups, which is in agreement with the actual number of chromosomes. Duplex and triplex markers were used to detect the homologous groups, and 13 and 14 homologous groups were identified in Xushu 18 and Xu 781 maps, respectively. Interval mapping was performed first and, subsequently, a multiple QTL model was used to refine the position and magnitude of the QTL. A total of 27 QTL for dry-matter content were mapped, explaining 9.0–45.1 % of the variation; 77.8 % of the QTL had a positive effect on the variation. This work represents an important step forward in genomics and marker-assisted breeding of sweetpotato.