SNP high-throughput screening in grapevine using the SNPlex™ genotyping system

Background

Until recently, only a small number of low- and mid-throughput methods have been used for single nucleotide polymorphism (SNP) discovery and genotyping in grapevine (Vitis vinifera L.). However, following completion of the sequence of the highly heterozygous genome of Pinot Noir, it has been possible to identify millions of electronic SNPs (eSNPs) thus providing a valuable source for high-throughput genotyping methods.

Results

Herein we report the first application of the SNPlex? genotyping system in grapevine aiming at the anchoring of an eukaryotic genome. This approach combines robust SNP detection with automated assay readout and data analysis. 813 candidate eSNPs were developed from non-repetitive contigs of the assembled genome of Pinot Noir and tested in 90 progeny of Syrah × Pinot Noir cross. 563 new SNP-based markers were obtained and mapped. The efficiency rate of 69% was enhanced to 80% when multiple displacement amplification (MDA) methods were used for preparation of genomic DNA for the SNPlex assay.

Conclusion

Unlike other SNP genotyping methods used to investigate thousands of SNPs in a few genotypes, or a few SNPs in around a thousand genotypes, the SNPlex genotyping system represents a good compromise to investigate several hundred SNPs in a hundred or more samples simultaneously. Therefore, the use of the SNPlex assay, coupled with whole genome amplification (WGA), is a good solution for future applications in well-equipped laboratories.     

Analysis of TNFAIP3, a feedback inhibitor of nuclear factor-κB and the neighbor intergenic 6q23 region in rheumatoid arthritis susceptibility

Introduction

Genome-wide association studies of rheumatoid arthritis (RA) have identified an association of the disease with a 6q23 region devoid of genes. TNFAIP3, an RA candidate gene, flanks this region, and polymorphisms in both the TNFAIP3 gene and the intergenic region are associated with systemic lupus erythematosus. We hypothesized that there is a similar association with RA, including polymorphisms in TNFAIP3 and the intergenic region.

Methods

To test this hypothesis, we selected tag-single nucleotide polymorphisms (SNPs) in both loci. They were analyzed in 1,651 patients with RA and 1,619 control individuals of Spanish ancestry.

Results

Weak evidence of association was found both in the 6q23 intergenic region and in the TNFAIP3 locus. The rs582757 SNP and a common haplotype in the TNFAIP3 locus exhibited association with RA. In the intergenic region, two SNPs were associated, namely rs609438 and rs13207033. The latter was only associated in patients with anti-citrullinated peptide antibodies. Overall, statistical association was best explained by the interdependent contribution of SNPs from the two loci TNFAIP3 and the 6q23 intergenic region.

Conclusions

Our data are consistent with the hypothesis that several RA genetic factors exist in the 6q23 region, including polymorphisms in the TNFAIP3 gene, like that previously described for systemic lupus erythematosus.     

FitSNPs: highly differentially expressed genes are more likely to have variants associated with disease

Background

Candidate single nucleotide polymorphisms (SNPs) from genome-wide association studies (GWASs) were often selected for validation based on their functional annotation, which was inadequate and biased. We propose to use the more than 200,000 microarray studies in the Gene Expression Omnibus to systematically prioritize candidate SNPs from GWASs.

Results

We analyzed all human microarray studies from the Gene Expression Omnibus, and calculated the observed frequency of differential expression, which we called differential expression ratio, for every human gene. Analysis conducted in a comprehensive list of curated disease genes revealed a positive association between differential expression ratio values and the likelihood of harboring disease-associated variants. By considering highly differentially expressed genes, we were able to rediscover disease genes with 79% specificity and 37% sensitivity. We successfully distinguished true disease genes from false positives in multiple GWASs for multiple diseases. We then derived a list of functionally interpolating SNPs (fitSNPs) to analyze the top seven loci of Wellcome Trust Case Control Consortium type 1 diabetes mellitus GWASs, rediscovered all type 1 diabetes mellitus genes, and predicted a novel gene (KIAA1109) for an unexplained locus 4q27. We suggest that fitSNPs would work equally well for both Mendelian and complex diseases (being more effective for cancer) and proposed candidate genes to sequence for their association with 597 syndromes with unknown molecular basis.

Conclusions

Our study demonstrates that highly differentially expressed genes are more likely to harbor disease-associated DNA variants. FitSNPs can serve as an effective tool to systematically prioritize candidate SNPs from GWASs.     

Highly efficient genotyping of rice biparental populations by GoldenGate assays based on parental resequencing

Key message

A new time- and cost-effective strategy was developed for medium-density SNP genotyping of rice biparental populations, using GoldenGate assays based on parental resequencing.

Abstract

Since the advent of molecular markers, crop researchers and breeders have dedicated huge amounts of effort to detecting quantitative trait loci (QTL) in biparental populations for genetic analysis and marker-assisted selection (MAS). In this study, we developed a new time- and cost-effective strategy for genotyping a population of progeny from a rice cross using medium-density single nucleotide polymorphisms (SNPs). Using this strategy, 728,362 “high quality” SNPs were identified by resequencing Teqing and Lemont, the parents of the population. We selected 384 informative SNPs that were evenly distributed across the genome for genotyping the biparental population using the Illumina GoldenGate assay. 335 (87.2 %) validated SNPs were used for further genetic analyses. After removing segregation distortion markers, 321 SNPs were used for linkage map construction and QTL mapping. This strategy generated SNP markers distributed more evenly across the genome than previous SSR assays. Taking the GW5 gene that controls grain shape as an example, our strategy provided higher accuracy (0.8 Mb) and significance (LOD 5.5 and 10.1) in QTL mapping than SSR analysis. Our study thus provides a rapid and efficient strategy for genetic studies and QTL mapping using SNP genotyping assays.     

Investigating the utility of combining Φ29 whole genome amplification and highly multiplexed single nucleotide polymorphism BeadArray™ genotyping

Background

Sustainable DNA resources and reliable high-throughput genotyping methods are required for large-scale, long-term genetic association studies. In the genetic dissection of common disease it is now recognised that thousands of samples and hundreds of thousands of markers, mostly single nucleotide polymorphisms (SNPs), will have to be analysed. In order to achieve these aims, both an ability to boost quantities of archived DNA and to genotype at low costs are highly desirable. We have investigated Φ29 polymerase Multiple Displacement Amplification (MDA)-generated DNA product (MDA product), in combination with highly multiplexed BeadArray? genotyping technology. As part of a large-scale BeadArray genotyping experiment we made a direct comparison of genotyping data generated from MDA product with that from genomic DNA (gDNA) templates.

Results

Eighty-six MDA product and the corresponding 86 gDNA samples were genotyped at 345 SNPs and a concordance rate of 98.8% was achieved. The BeadArray sample exclusion rate, blind to sample type, was 10.5% for MDA product compared to 5.8% for gDNA.

Conclusions

We conclude that the BeadArray technology successfully produces high quality genotyping data from MDA product. The combination of these technologies improves the feasibility and efficiency of mapping common disease susceptibility genes despite limited stocks of gDNA samples.     

Genetic analyses of bolting in bulb onion (Allium cepa L.)

Key message

We present the first evidence for a QTL conditioning an adaptive trait in bulb onion, and the first linkage and population genetics analyses of candidate genes involved in photoperiod and vernalization physiology.

Abstract

Economic production of bulb onion (Allium cepa L.) requires adaptation to photoperiod and temperature such that a bulb is formed in the first year and a flowering umbel in the second. ‘Bolting’, or premature flowering before bulb maturation, is an undesirable trait strongly selected against by breeders during adaptation of germplasm. To identify genome regions associated with adaptive traits we conducted linkage mapping and population genetic analyses of candidate genes, and QTL analysis of bolting using a low-density linkage map. We performed tagged amplicon sequencing of ten candidate genes, including the FT-like gene family, in eight diverse populations to identify polymorphisms and seek evidence of differentiation. Low nucleotide diversity and negative estimates of Tajima’s D were observed for most genes, consistent with purifying selection. Significant population differentiation was observed only in AcFT2 and AcSOC1. Selective genotyping in a large ‘Nasik Red × CUDH2150’ F₂ family revealed genome regions on chromosomes 1, 3 and 6 associated (LOD > 3) with bolting. Validation genotyping of two F₂ families grown in two environments confirmed that a QTL on chromosome 1, which we designate AcBlt1, consistently conditions bolting susceptibility in this cross. The chromosome 3 region, which coincides with a functionally characterised acid invertase, was not associated with bolting in other environments, but showed significant association with bulb sucrose content in this and other mapping pedigrees. These putative QTL and candidate genes were placed on the onion map, enabling future comparative studies of adaptive traits.     

Candidate gene association mapping of Sclerotinia stalk rot resistance in sunflower (Helianthus annuus L.) uncovers the importance of COI1 homologs

Key message

Functional markers for Sclerotinia basal stalk rot resistance in sunflower were obtained using gene-level information from the model species Arabidopsis thaliana.

Abstract

Sclerotinia stalk rot, caused by Sclerotinia sclerotiorum, is one of the most destructive diseases of sunflower (Helianthus annuus L.) worldwide. Markers for genes controlling resistance to S. sclerotiorum will enable efficient marker-assisted selection (MAS). We sequenced eight candidate genes homologous to Arabidopsis thaliana defense genes known to be associated with Sclerotinia disease resistance in a sunflower association mapping population evaluated for Sclerotinia stalk rot resistance. The total candidate gene sequence regions covered a concatenated length of 3,791 bp per individual. A total of 187 polymorphic sites were detected for all candidate gene sequences, 149 of which were single nucleotide polymorphisms (SNPs) and 38 were insertions/deletions. Eight SNPs in the coding regions led to changes in amino acid codons. Linkage disequilibrium decay throughout the candidate gene regions declined on average to an r ² = 0.2 for genetic intervals of 120 bp, but extended up to 350 bp with r ² = 0.1. A general linear model with modification to account for population structure was found the best fitting model for this population and was used for association mapping. Both HaCOI1-1 and HaCOI1-2 were found to be strongly associated with Sclerotinia stalk rot resistance and explained 7.4 % of phenotypic variation in this population. These SNP markers associated with Sclerotinia stalk rot resistance can potentially be applied to the selection of favorable genotypes, which will significantly improve the efficiency of MAS during the development of stalk rot resistant cultivars.     

Antimicrobial resistance of Escherichia coli isolates from broiler chickens and humans

Background

Inguinal hernias are usually caused by a congenital defect, which occurs as a weakness of the inguinal canal. Porcine β-glucuronidase gene (GUSB) was chosen as functional candidate gene because of its involvement in degradation of hyaluronan within gubernacular tissue during descent of testes. Since a genome-wide linkage analysis approach has shown evidence that two regions on porcine chromosome 3 (SSC 3) are involved in the inheritance of hernia inguinalis/scrotalis in German pig breeds, GUSB also attained status as a positional candidate gene by its localization within a hernia-associated chromosomal region.

Results

A contig spanning 17,157 bp, which contains the entire GUSB, was assembled. Comparative sequence analyses were conducted for the GUSB gene locus. Single nucleotide polymorphisms (SNPs) located within the coding region of GUSB were genotyped in 512 animals. Results of transmission disequilibrium test (TDT) for two out of a total of five detected SNPs gave no significant association with the outcome of hernia in pigs.

Conclusion

On the basis of our studies we are able to exclude the two analyzed SNPs within the porcine GUSB gene as causative for the transmission of inguinal hernia.     

Deletion of the V2 vasopressin receptor gene in two Chinese patients with nephrogenic diabetes insipidus

Background

The aetiology of the autoimmune disease type 1 diabetes (T1D) involves many genetic and environmental factors. Evidence suggests that innate immune responses, including the action of interferons, may also play a role in the initiation and/or pathogenic process of autoimmunity. In the present report, we have adopted a linkage disequilibrium (LD) mapping approach to test for an association between T1D and three regions encompassing 13 interferon alpha (IFNA) genes, interferon omega-1 (IFNW1), interferon beta-1 (IFNB1), interferon gamma (IFNG) and the interferon consensus-sequence binding protein 1 (ICSBP1).

Results

We identified 238 variants, most, single nucleotide polymorphisms (SNPs), by sequencing IFNA, IFNB1, IFNW1 and ICSBP1, 98 of which where novel when compared to dbSNP build 124. We used polymorphisms identified in the SeattleSNP database for INFG. A set of tag SNPs was selected for each of the interferon and interferon-related genes to test for an association between T1D and this complex gene family. A total of 45 tag SNPs were selected and genotyped in a collection of 472 multiplex families.

Conclusion

We have developed informative sets of SNPs for the interferon and interferon related genes. No statistical evidence of a major association between T1D and any of the interferon and interferon related genes tested was found.     

Rare allele of HvLox-1 associated with lipoxygenase activity in barley (Hordeum vulgare L.)

Key message

Identification and allele-specific marker development of a functional SNP of HvLox - 1 which associated with barley lipoxygenase activity.

Abstract

Improving the stability of the flavor of beer is one of the main objectives in breeding barley for malting, and lipoxygenase-1 (LOX-1) is a key enzyme controlling this trait. In this study, a modified LOX activity assay was used for null LOX-1 mutant screening. Four barley landraces with no detected level of LOX-1 activity were screened from 1,083 barley germplasm accessions from China. The genomic sequence diversity of the HvLox-1 gene of the four null LOX-1 Chinese landraces was compared with that of a further 76 accessions. A total of 104 nucleotide polymorphisms were found, which contained 83 single-nucleotide polymorphisms (SNPs), 7 multiple-nucleotide polymorphisms, and 14 insertions and deletions. Most notably, we found a rare C/G mutation (SNP-61) in the second intron which led to null LOX-1 activity through an altered splicing acceptor site. In addition, an allele-specific polymerase chain reaction marker was developed for the genotyping of SNP-61, which could be used in breeding programs for barley to be used for malting. The objective was to improve beer quality.     

Use of synteny to identify candidate genes underlying QTL controlling stomatal traits in faba bean (Vicia faba L.)

Key message

We have identified QTLs for stomatal characteristics on chromosome II of faba bean by applying SNPs derived from M. truncatula , and have identified candidate genes within these QTLs using synteny between the two species.

Abstract

Faba bean (Vicia faba L.) is a valuable food and feed crop worldwide, but drought often limits its production, and its genome is large and poorly mapped. No information is available on the effects of genomic regions and genes on drought adaptation characters such as stomatal characteristics in this species, but the synteny between the sequenced model legume, Medicago truncatula, and faba bean can be used to identify candidate genes. A mapping population of 211 F₅ recombinant inbred lines (Mélodie/2 × ILB 938/2) were phenotyped to identify quantitative trait loci (QTL) affecting stomatal morphology and function, along with seed weight, under well-watered conditions in a climate-controlled glasshouse in 2013 and 2014. Canopy temperature (CT) was evaluated in 2013 under water-deficit (CTd). In total, 188 polymorphic single nucleotide polymorphisms (SNPs), developed from M. truncatula genome data, were assigned to nine linkage groups that covered ~928 cM of the faba bean genome with an average inter-marker distance of 5.8 cM. 15 putative QTLs were detected, of which eight (affecting stomatal density, length and conductance and CT) co-located on chromosome II, in the vicinity of a possible candidate gene—a receptor-like protein kinase found in the syntenic interval of M. truncatula chromosome IV. A ribose-phosphate pyrophosphokinase from M. truncatula chromosome V, postulated as a possible candidate gene for the QTL for CTd, was found some distance away in the same chromosome. These results demonstrate that genomic information from M. truncatula can successfully be translated to the faba bean genome.     

Genetical genomic determinants of alcohol consumption in rats and humans

Background

We have used a genetical genomic approach, in conjunction with phenotypic analysis of alcohol consumption, to identify candidate genes that predispose to varying levels of alcohol intake by HXB/BXH recombinant inbred rat strains. In addition, in two populations of humans, we assessed genetic polymorphisms associated with alcohol consumption using a custom genotyping array for 1,350 single nucleotide polymorphisms (SNPs). Our goal was to ascertain whether our approach, which relies on statistical and informatics techniques, and non-human animal models of alcohol drinking behavior, could inform interpretation of genetic association studies with human populations.

Results

In the HXB/BXH recombinant inbred (RI) rats, correlation analysis of brain gene expression levels with alcohol consumption in a two-bottle choice paradigm, and filtering based on behavioral and gene expression quantitative trait locus (QTL) analyses, generated a list of candidate genes. A literature-based, functional analysis of the interactions of the products of these candidate genes defined pathways linked to presynaptic GABA release, activation of dopamine neurons, and postsynaptic GABA receptor trafficking, in brain regions including the hypothalamus, ventral tegmentum and amygdala. The analysis also implicated energy metabolism and caloric intake control as potential influences on alcohol consumption by the recombinant inbred rats. In the human populations, polymorphisms in genes associated with GABA synthesis and GABA receptors, as well as genes related to dopaminergic transmission, were associated with alcohol consumption.

Conclusion

Our results emphasize the importance of the signaling pathways identified using the non-human animal models, rather than single gene products, in identifying factors responsible for complex traits such as alcohol consumption. The results suggest cross-species similarities in pathways that influence predisposition to consume alcohol by rats and humans. The importance of a well-defined phenotype is also illustrated. Our results also suggest that different genetic factors predispose alcohol dependence versus the phenotype of alcohol consumption.     

Identification of novel genetic susceptibility loci for Behçet's disease using a genome-wide association study

Introduction

Behçet's disease is a chronic systemic inflammatory disease that remains incompletely understood. Herein, we perform the first genome-wide association study in Behçet's disease.

Methods

Using DNA pooling technology and the Affymetrix 500K arrays, we identified possible candidate gene associations with Behçet's disease in a cohort of 152 Behçet's disease patients and 172 healthy ethnically matched controls. Genetic loci that were identified in the pooling study were genotyped in patients and controls using TaqMan genotyping technology.

Results

We identified genetic associations between Behçet's disease and single-nucleotide polymorphisms (SNPs) in KIAA1529, CPVL, LOC100129342, UBASH3B, and UBAC2 (odds ratio = 2.04, 2.26, 1.84, 1.71, and 1.61, respectively; P value = 4.2 × 10^-5, 1.0 × 10^-4, 3.0 × 10^-4, 1.5 × 10^-3, and 5.8 × 10^-3, respectively). Among the associated SNPs, the Behçet's disease-risk allele in rs2061634 leads to substitution of serine to cysteine at amino acid position 995 (S995C) in the KIAA1529 protein.

Conclusions

Using an unbiased whole-genome genetic association approach, we identified novel candidate genetic loci that are associated with increased susceptibility for Behçet's disease. These findings will help to better understand the pathogenesis of Behçet's disease and identify novel targets for therapeutic intervention.     

Midkine is a NF-κB-inducible gene that supports prostate cancer cell survival

Background

Atherosclerotic peripheral arterial disease (PAD) affects 8–10 million people in the United States and is associated with a marked impairment in quality of life and an increased risk of cardiovascular events. Noninvasive assessment of PAD is performed by measuring the ankle-brachial index (ABI). Complex traits, such as ABI, are influenced by a large array of genetic and environmental factors and their interactions. We attempted to characterize the genetic architecture of ABI by examining the main and interactive effects of individual single nucleotide polymorphisms (SNPs) and conventional risk factors.

Methods

We applied linear regression analysis to investigate the association of 435 SNPs in 112 positional and biological candidate genes with ABI and related physiological and biochemical traits in 1046 non-Hispanic white, hypertensive participants from the Genetic Epidemiology Network of Arteriopathy (GENOA) study. The main effects of each SNP, as well as SNP-covariate and SNP-SNP interactions, were assessed to investigate how they contribute to the inter-individual variation in ABI. Multivariable linear regression models were then used to assess the joint contributions of the top SNP associations and interactions to ABI after adjustment for covariates. We reduced the chance of false positives by 1) correcting for multiple testing using the false discovery rate, 2) internal replication, and 3) four-fold cross-validation.

Results

When the results from these three procedures were combined, only two SNP main effects in NOS3, three SNP-covariate interactions (ADRB2 Gly 16 – lipoprotein(a) and SLC4A5 – diabetes interactions), and 25 SNP-SNP interactions (involving SNPs from 29 different genes) were significant, replicated, and cross-validated. Combining the top SNPs, risk factors, and their interactions into a model explained nearly 18% of variation in ABI in the sample. SNPs in six genes (ADD2, ATP6V1B1, PRKAR2B, SLC17A2, SLC22A3, and TGFB3) were also influencing triglycerides, C-reactive protein, homocysteine, and lipoprotein(a) levels.

Conclusion

We found that candidate gene SNP main effects, SNP-covariate and SNP-SNP interactions contribute to the inter-individual variation in ABI, a marker of PAD. Our findings underscore the importance of conducting systematic investigations that consider context-dependent frameworks for developing a deeper understanding of the multidimensional genetic and environmental factors that contribute to complex diseases.     

Association of acid phosphatase locus 1*C allele with the risk of cardiovascular events in rheumatoid arthritis patients

Introduction

Acid phosphatase locus 1 (ACP1) encodes a low molecular weight phosphotyrosine phosphatase implicated in a number of different biological functions in the cell. The aim of this study was to determine the contribution of ACP1 polymorphisms to susceptibility to rheumatoid arthritis (RA), as well as the potential contribution of these polymorphisms to the increased risk of cardiovascular disease (CV) observed in RA patients.

Methods

A set of 1,603 Spanish RA patients and 1,877 healthy controls were included in the study. Information related to the presence/absence of CV events was obtained from 1,284 of these participants. All individuals were genotyped for four ACP1 single-nucleotide polymorphisms (SNPs), rs10167992, rs11553742, rs7576247, and rs3828329, using a predesigned TaqMan SNP genotyping assay. Classical ACP1 alleles (*A, *B and *C) were imputed with SNP data.

Results

No association between ACP1 gene polymorphisms and susceptibility to RA was observed. However, when RA patients were stratified according to the presence or absence of CV events, an association between rs11553742*T and CV events was found (P = 0.012, odds ratio (OR) = 2.62 (1.24 to 5.53)). Likewise, the ACP1*C allele showed evidence of association with CV events in patients with RA (P = 0.024, OR = 2.43).

Conclusions

Our data show that the ACP1*C allele influences the risk of CV events in patients with RA.