Development,validation and genetic analysis of a large soybean SNP genotyping array

Cultivated soybean (Glycine max) suffers from a narrow germplasm relative to other crop species, probably because of under‐use of wild soybean (Glycine soja) as a breeding resource. Use of a single nucleotide polymorphism (SNP) genotyping array is a promising method for dissecting cultivated and wild germplasms to identify important adaptive genes through high‐density genetic mapping and genome‐wide association studies. Here we describe a large soybean SNP array for use in diversity analyses, linkage mapping and genome‐wide association analyses. More than four million high‐quality SNPs identified from high‐depth genome re‐sequencing of 16 soybean accessions and low‐depth genome re‐sequencing of 31 soybean accessions were used to select 180 961 SNPs for creation of the Axiom^® SoyaSNP array. Validation analysis for a set of 222 diverse soybean lines showed that 170 223 markers were of good quality for genotyping. Phylogenetic and allele frequency analyses of the validation set data indicated that accessions showing an intermediate morphology between cultivated and wild soybeans collected in Korea were natural hybrids. More than 90 unanchored scaffolds in the current soybean reference sequence were assigned to chromosomes using this array. Finally, dense average spacing and preferential distribution of the SNPs in gene‐rich chromosomal regions suggest that this array may be suitable for genome‐wide association studies of soybean germplasm. Taken together, these results suggest that use of this array may be a powerful method for soybean genetic analyses relating to many aspects of soybean breeding.     

Evaluation of the metabochip genotyping array in African Americans and implications for fine mapping of GWAS-identified loci: the PAGE study

The Metabochip is a custom genotyping array designed for replication and fine mapping of metabolic, cardiovascular, and anthropometric trait loci and includes low frequency variation content identified from the 1000 Genomes Project. It has 196,725 SNPs concentrated in 257 genomic regions. We evaluated the Metabochip in 5,863 African Americans; 89% of all SNPs passed rigorous quality control with a call rate of 99.9%. Two examples illustrate the value of fine mapping with the Metabochip in African-ancestry populations. At CELSR2/PSRC1/SORT1, we found the strongest associated SNP for LDL-C to be rs12740374 (p = 3.5 × 10(-11)), a SNP indistinguishable from multiple SNPs in European ancestry samples due to high correlation. Its distinct signal supports functional studies elsewhere suggesting a causal role in LDL-C. At CETP we found rs17231520, with risk allele frequency 0.07 in African Americans, to be associated with HDL-C (p = 7.2 × 10(-36)). This variant is very rare in Europeans and not tagged in common GWAS arrays, but was identified as associated with HDL-C in African Americans in a single-gene study. Our results, one narrowing the risk interval and the other revealing an associated variant not found in Europeans, demonstrate the advantages of high-density genotyping of common and rare variation for fine mapping of trait loci in African American samples.     

Tm-shift基因分型方法在遗传学中的应用

Melting Temperature shift(Tm-shift)是一种新的基因分型方法,主要通过在两条特异性引物5′端加入不同长度的GC序列,PCR扩增后根据熔解曲线中产物Tm值的差异来完成分型。文章建立了Tm-shift法对2 048份样品的29个SNP进行分型,通过分型成功率、重复检测一致率、测序验证准确度综合评价分型效果。结果显示,29个SNP中有27个可以采用本方法分型,分型成功率为93.1%。测序验证准确性达到100%。3种基因型阳性标准对照重复检测一致率为100%;100个随机样品重复检测,重复性为97%。因此,Tm-shift基因分型法是一种成本低廉、准确灵敏、稳定可靠、通量灵活、操作简便的基因分型方法,可在遗传学研究中推广应用。     

A new method for quantifying genotyping errors for noninvasive genetic studies

More and more noninvasive genetic data are being produced but a general methodology to quantify genotyping error rates from non-pilot data remains lacking. Here we propose a mathematical approach to estimate genotyping error rates by exploring the relationship between errors and PCR replicates. This method can be used to quantify the error rates for either the multi-tubes approach designed by Taberlet et al. (Nucleic Acids Res 24: 3189–3194, 1996) or the pilot method by Prugh et al. (Mol Ecol 14: 1585–1596, 2005).     

Heritability and genetic correlations explained by common SNPs for metabolic syndrome traits

We used a bivariate (multivariate) linear mixed-effects model to estimate the narrow-sense heritability (h(2)) and heritability explained by the common SNPs (h(g)(2)) for several metabolic syndrome (MetS) traits and the genetic correlation between pairs of traits for the Atherosclerosis Risk in Communities (ARIC) genome-wide association study (GWAS) population. MetS traits included body-mass index (BMI), waist-to-hip ratio (WHR), systolic blood pressure (SBP), fasting glucose (GLU), fasting insulin (INS), fasting trigylcerides (TG), and fasting high-density lipoprotein (HDL). We found the percentage of h(2) accounted for by common SNPs to be 58% of h(2) for height, 41% for BMI, 46% for WHR, 30% for GLU, 39% for INS, 34% for TG, 25% for HDL, and 80% for SBP. We confirmed prior reports for height and BMI using the ARIC population and independently in the Framingham Heart Study (FHS) population. We demonstrated that the multivariate model supported large genetic correlations between BMI and WHR and between TG and HDL. We also showed that the genetic correlations between the MetS traits are directly proportional to the phenotypic correlations.     

Concept, design and implementation of a cardiovascular gene-centric 50 k SNP array for large-scale genomic association studies

A wealth of genetic associations for cardiovascular and metabolic phenotypes in humans has been accumulating over the last decade, in particular a large number of loci derived from recent genome wide association studies (GWAS). True complex disease-associated loci often exert modest effects, so their delineation currently requires integration of diverse phenotypic data from large studies to ensure robust meta-analyses. We have designed a gene-centric 50 K single nucleotide polymorphism (SNP) array to assess potentially relevant loci across a range of cardiovascular, metabolic and inflammatory syndromes. The array utilizes a "cosmopolitan" tagging approach to capture the genetic diversity across approximately 2,000 loci in populations represented in the HapMap and SeattleSNPs projects. The array content is informed by GWAS of vascular and inflammatory disease, expression quantitative trait loci implicated in atherosclerosis, pathway based approaches and comprehensive literature searching. The custom flexibility of the array platform facilitated interrogation of loci at differing stringencies, according to a gene prioritization strategy that allows saturation of high priority loci with a greater density of markers than the existing GWAS tools, particularly in African HapMap samples. We also demonstrate that the IBC array can be used to complement GWAS, increasing coverage in high priority CVD-related loci across all major HapMap populations. DNA from over 200,000 extensively phenotyped individuals will be genotyped with this array with a significant portion of the generated data being released into the academic domain facilitating in silico replication attempts, analyses of rare variants and cross-cohort meta-analyses in diverse populations. These datasets will also facilitate more robust secondary analyses, such as explorations with alternative genetic models, epistasis and gene-environment interactions.     

Archival,demographic and genetic studies define a Sardinian sub-isolate as a suitable model for mapping complex traits

Genetic isolates represent exceptional resources for the mapping of complex traits but not all isolates are similar. We have selected a genetic and cultural isolate, the village of Talana from an isolated area of Sardinia, and propose that this population is suitable for the mapping of complex traits. A wealth of historical and archive data allowed the reconstruction of the demographic and genealogical history of the village. Key features of the population, which has grown slowly with no significant immigration, were defined by using a combination of historical, demographic and genetic studies. The genealogy of each Talana inhabitant was reconstructed and the main maternal and paternal lineages of the village were defined. Haplotype and phylogenetic analyses of the Y chromosome and characterisation of mitochondrial DNA haplogroups were used to determine the number of ancestral village founders. The extent of linkage disequilibrium (LD) was evaluated by the analysis of several microsatellites in chromosomal region Xq13.3, which was previously used to asses the extension of LD. Genealogical reconstructions were confirmed and reinforced by the genetic analyses, since some lineages were found to have merged prior to the beginning of the archival records, suggesting an even smaller number of founders than initially predicted. About 80% of the present-day population appears to derive from eight paternal and eleven maternal ancestral lineages. LD was found to span, on average, a 5-Mb region in Xq13.3. This suggests the possibility of identifying identical-by-descent regions associated with complex traits in a genome-wide search by using a low-density marker map. The present study emphasises the importance of combining genetic studies with genealogical and historical information.     

Efficiency of selective genotyping for genetic analysis of complex traits and potential applications in crop improvement

Selective genotyping of individuals from the two tails of the phenotypic distribution of a population provides a cost efficient alternative to analysis of the entire population for genetic mapping. Past applications of this approach have been confounded by the small size of entire and tail populations, and insufficient marker density, which result in a high probability of false positives in the detection of quantitative trait loci (QTL). We studied the effect of these factors on the power of QTL detection by simulation of mapping experiments using population sizes of up to 3,000 individuals and tail population sizes of various proportions, and marker densities up to one marker per centiMorgan using complex genetic models including QTL linkage and epistasis. The results indicate that QTL mapping based on selective genotyping is more powerful than simple interval mapping but less powerful than inclusive composite interval mapping. Selective genotyping can be used, along with pooled DNA analysis, to replace genotyping the entire population, for mapping QTL with relatively small effects, as well as linked and interacting QTL. Using diverse germplasm including all available genetics and breeding materials, it is theoretically possible to develop an “all-in-one plate” approach where one 384-well plate could be designed to map almost all agronomic traits of importance in a crop species. Selective genotyping can also be used for genomewide association mapping where it can be integrated with selective phenotyping approaches. We also propose a breeding-to-genetics approach, which starts with identification of extreme phenotypes from segregating populations generated from multiple parental lines and is followed by rapid discovery of individual genes and combinations of gene effects together with simultaneous manipulation in breeding programs.     

Generalized indicator plate for genetic, metabolic, and taxonomic studies with microorganisms.

We have developed an indicator plate that works well for diverse types of substrates and microorganisms. The plates are inexpensive and easy to prepare. The essential components are agar, buffer, growth-supporting nutrients, a test substrate, and 2,3,5-triphenyl tetrazolium chloride (TTC). Using various strains of Salmonella typhimurium and Escherichia coli, we have studied and defined the contribution of each component to the satisfactory function of the plate. Colonies capable of catabolizing the test substrate reduce TTC and produce a deep red formazan, whereas colonies failing to catabolize the substrate remain uncoloured. Those with intermediate rates of catabolism differ in rate and/or extent of color formation. In all cases the color is stable because TTC reduction is essentially irreversible. Since the mode of action of these plates is fairly well understood, alternative formulations can be devised to meet specific needs. The general applicability of this TTC indicator system makes it an extremely useful tool in microbial genetics, metabolism, and taxonomy.     

Generalized indicator plate for genetic, metabolic, and taxonomic studies with microorganisms.

We have developed an indicator plate that works well for diverse types of substrates and microorganisms. The plates are inexpensive and easy to prepare. The essential components are agar, buffer, growth-supporting nutrients, a test substrate, and 2,3,5-triphenyl tetrazolium chloride (TTC). Using various strains of Salmonella typhimurium and Escherichia coli, we have studied and defined the contribution of each component to the satisfactory function of the plate. Colonies capable of catabolizing the test substrate reduce TTC and produce a deep red formazan, whereas colonies failing to catabolize the substrate remain uncoloured. Those with intermediate rates of catabolism differ in rate and/or extent of color formation. In all cases the color is stable because TTC reduction is essentially irreversible. Since the mode of action of these plates is fairly well understood, alternative formulations can be devised to meet specific needs. The general applicability of this TTC indicator system makes it an extremely useful tool in microbial genetics, metabolism, and taxonomy.     

Development of a 135K SNP genotyping array for Actinidia arguta and its applications for genetic mapping and QTL analysis in kiwifruit

Kiwifruit (Actinidia spp) is a woody, perennial and deciduous vine. In this genus, there are multiple ploidy levels but the main cultivated cultivars are polyploid. Despite the availability of many genomic resources in kiwifruit, SNP genotyping is still a challenge given these different levels of polyploidy. Recent advances in SNP array technologies have offered a high-throughput genotyping platform for genome-wide DNA polymorphisms. In this study, we developed a high-density SNP genotyping array to facilitate genetic studies and breeding applications in kiwifruit. SNP discovery was performed by genome-wide DNA sequencing of 40 kiwifruit genotypes. The identified SNPs were stringently filtered for sequence quality, predicted conversion performance and distribution over the available Actinidia chinensis genome. A total of 134 729 unique SNPs were put on the array. The array was evaluated by genotyping 400 kiwifruit individuals. We performed a multidimensional scaling analysis to assess the diversity of kiwifruit germplasm, showing that the array was effective to distinguish kiwifruit accessions. Using a tetraploid F1 population, we constructed an integrated linkage map covering 3060.9 cM across 29 linkage groups and performed QTL analysis for the sex locus that has been identified on Linkage Group 3 (LG3) in Actinidia arguta. Finally, our dataset presented evidence of tetrasomic inheritance with partial preferential pairing in A. arguta. In conclusion, we developed and evaluated a 135K SNP genotyping array for kiwifruit. It has the advantage of a comprehensive design that can be an effective tool in genetic studies and breeding applications in this high-value crop.     

Development of an integrated 200K SNP genotyping array and application for genetic mapping,genome assembly improvement and genome wide association studies in pear (Pyrus)

Pear (Pyrus; 2n = 34), the third most important temperate fruit crop, has great nutritional and economic value. Despite the availability of many genomic resources in pear, it is challenging to genotype novel germplasm resources and breeding progeny in a timely and cost‐effective manner. Genotyping arrays can provide fast, efficient and high‐throughput genetic characterization of diverse germplasm, genetic mapping and breeding populations. We present here 200K AXIOM^® PyrSNP, a large‐scale single nucleotide polymorphism (SNP) genotyping array to facilitate genotyping of Pyrus species. A diverse panel of 113 re‐sequenced pear genotypes was used to discover SNPs to promote increased adoption of the array. A set of 188 diverse accessions and an F₁ population of 98 individuals from ‘Cuiguan’ × ‘Starkrimson’ was genotyped with the array to assess its effectiveness. A large majority of SNPs (166 335 or 83%) are of high quality. The high density and uniform distribution of the array SNPs facilitated prediction of centromeric regions on 17 pear chromosomes, and significantly improved the genome assembly from 75.5% to 81.4% based on genetic mapping. Identification of a gene associated with flowering time and candidate genes linked to size of fruit core via genome wide association studies showed the usefulness of the array in pear genetic research. The newly developed high‐density SNP array presents an important tool for rapid and high‐throughput genotyping in pear for genetic map construction, QTL identification and genomic selection.     

Ascertainment adjustment in genetic studies of ordinal traits

Most genetic studies recruit high risk families and the discoveries are based on non-random selected groups. We must consider the consequences of this ascertainment process in order to apply the results of genetic research to the general population. In previous reports, we developed a latent variable model to assess the familial aggregation and inheritability of ordinal-scaled diseases, and found a major gene component of alcoholism after applying the model to the data from the Yale family study of comorbidity of alcoholism and anxiety (YFSCAA). In this report, we examine the ascertainment effects on parameter estimates and correct potential bias in the latent variable model. The simulation studies for various ascertainment schemes suggest that our ascertainment adjustment is necessary and effective. We also find that the estimated effects are relatively unbiased for the particular ascertainment scheme used in the YFSCAA, which assures the validity of our earlier conclusion.     

Mutation and copy number detection in human cancers using a custom genotyping assay

Identification of biomarkers for positive and negative predictors of response to cancer therapeutics can help direct clinical strategies. However, challenges with tissue availability and costs are significant limiting factors for diagnostic assays. To address these challenges, we have customized a high-throughput single nucleotide polymorphism genotyping assay with the objective of simultaneously surveying known somatic mutations and copy number alterations for translational studies in cancer. As constructed, this assay can interrogate 376 known somatic mutations and quantify copy number alterations of genes commonly implicated in tumorigenesis or progression. Validation of this assay on a panel of 321 cell lines demonstrates sensitivity to accurately detect mutations, robust accuracy in the presence of infiltrating normal tissue, and the ability to detect both DNA copy number amplifications and deletions. This technology, with its high sensitivity, small DNA requirements, and low costs is an attractive platform for biomarker exploration in cancer.     

Familial resemblance for anthropometric traits in dizygotic twins and siblings

Familial resemblance for fifty anthropometric traits was studied on a sample of 45 MZ, 101 DZ twin pairs and their 125 singleton siblings. Intraclass correlation coefficients were significant for all the traits. However, resemblance within DZ twin pairs was significantly greater than within sibs for 22 variables, showing that the former had a more correlated environment than the latter. The study also showed that head and facial traits were relatively more stable to the environmental factors than the body traits and hence more suitable for cross-cultural comparisons. The study listed measures of girth and skinfold, thickness as the most labile traits.