Family-based association analysis with ordered categorical phenotypes, covariates and interactions

Genetic association analyses of family-based studies with ordered categorical phenotypes are often conducted using methods either for quantitative or for binary traits, which can lead to suboptimal analyses. Here we present an alternative likelihood-based method of analysis for single nucleotide polymorphism (SNP) genotypes and ordered categorical phenotypes in nuclear families of any size. Our approach, which extends our previous work for binary phenotypes, permits straightforward inclusion of covariate, gene-gene and gene-covariate interaction terms in the likelihood, incorporates a simple model for ascertainment and allows for family-specific effects in the hypothesis test. Additionally, our method produces interpretable parameter estimates and valid confidence intervals. We assess the proposed method using simulated data, and apply it to a polymorphism in the c-reactive protein (CRP) gene typed in families collected to investigate human systemic lupus erythematosus. By including sex interactions in the analysis, we show that the polymorphism is associated with anti-nuclear autoantibody (ANA) production in females, while there appears to be no effect in males.

Sequentially testing for a gene-drug interaction in a genomewide analysis

Assaying a large number of genetic markers from patients in clinical trials is now possible in order to tailor drugs with respect to efficacy. The statistical methodology for analysing such massive data sets is challenging. The most popular type of statistical analysis is to use a univariate test for each genetic marker, once all the data from a clinical study have been collected. This paper presents a sequential method for conducting an omnibus test for detecting gene-drug interactions across the genome, thus allowing informed decisions at the earliest opportunity and overcoming the multiple testing problems from conducting many univariate tests. We first propose an omnibus test for a fixed sample size. This test is based on combining F-statistics that test for an interaction between treatment and the individual single nucleotide polymorphism (SNP). As SNPs tend to be correlated, we use permutations to calculate a global p-value. We extend our omnibus test to the sequential case. In order to control the type I error rate, we propose a sequential method that uses permutations to obtain the stopping boundaries. The results of a simulation study show that the sequential permutation method is more powerful than alternative sequential methods that control the type I error rate, such as the inverse-normal method. The proposed method is flexible as we do not need to assume a mode of inheritance and can also adjust for confounding factors. An application to real clinical data illustrates that the method is computationally feasible for a large number of SNPs. Copyright (c) 2007 John Wiley & Sons, Ltd.

Recent developments in genetic data analysis: what can they tell us about human demographic history?

Over the last decade, a number of new methods of population genetic analysis based on likelihood have been introduced. This review describes and explains the general statistical techniques that have recently been used, and discusses the underlying population genetic models. Experimental papers that use these methods to infer human demographic and phylogeographic history are reviewed. It appears that the use of likelihood has hitherto had little impact in the field of human population genetics, which is still primarily driven by more traditional approaches. However, with the current uncertainty about the effects of natural selection, population structure and ascertainment of single-nucleotide polymorphism markers, it is suggested that likelihood-based methods may have a greater impact in the future.

Linkage disequilibrium assessment via log-linear modeling of SNP haplotype frequencies

Analyses of high-density single-nucleotide polymorphism (SNP) data, such as genetic mapping and linkage disequilibrium (LD) studies, require phase-known haplotypes to allow for the correlation between tightly linked loci. However, current SNP genotyping technology cannot determine phase, which must be inferred statistically. In this paper, we present a new Bayesian Markov chain Monte Carlo (MCMC) algorithm for population haplotype frequency estimation, particulary in the context of LD assessment. The novel feature of the method is the incorporation of a log-linear prior model for population haplotype frequencies. We present simulations to suggest that 1) the log-linear prior model is more appropriate than the standard coalescent process in the presence of recombination (>0.02cM between adjacent loci), and 2) there is substantial inflation in measures of LD obtained by a "two-stage" approach to the analysis by treating the "best" haplotype configuration as correct, without regard to uncertainty in the recombination process. Genet Epidemiol 25:106-114, 2003. (C) 2003 Wiley-Liss, Inc.

Influence of apoA-V gene variants on postprandial triglyceride metabolism: impact of gender

Although apolipoprotein AN (apoA-V) polymorphisms have been consistently associated with fasting triglyceride (TG) levels, their impact on postprandial lipemia remains relatively unknown. In this study, we investigate the impact of two common apoA-V polymorphisms (-1131 T>C and S19W) and apoA-V haplotypes on fasting and postprandial lipid metabolism in adults in the United Kingdom (n = 259). Compared with the wild-type TT, apoA-V -1131 TC heterozygotes had 15% (P = 0.057) and 21% (P = 0.002) higher fasting TG and postprandial TG area under the curve (AUC), respectively. Significant (P = 0.038) and nearly significant (P = 0.057) gender X genotype interactions were observed for fasting TG and TG AUC, with a greater impact of genotype in males. Lower HDL-cholesterol was associated with the rare TC genotype (P = 0.047). Significant linkage disequilibrium was found between the apoA-V -1131 T>C and the apoC-III 3238 C>G variants, with univariate analysis indicating an impact of this apoC-III single nucleotide polymorphism (SNP) on TG AUC (P = 0.015). However, in linear regression analysis, a significant independent association with TG AUC (P = 0.007) was only evident for the apoA-V -1131 T>C SNP, indicating a greater relative importance of the apoA-V genotype.

ApoE genotype, cardiovascular risk and responsiveness to dietary fat manipulation

Cardiovascular risk is determined by the complex interactions between genetic and environmental factors. The apoE genotype represents the most-widely-studied single nucleotide polymorphism in relation to CVD risk, with >3600 publications cited in PubMed. Although originally described as a mediator of lipoprotein metabolism, the lipoprotein-independent functions of apoE are being increasingly recognised, with limited data available on the potential impact of genotype on these metabolic processes. Furthermore, although meta-analyses suggest that apoE4 carriers may have a 40-50% increased CVD risk, the associations reported in individual studies are highly heterogeneous and it is recognised that environmental factors such as smoking status and dietary fat composition influence genotype-phenotype associations. However, information is often derived from observational studies or small intervention trials in which retrospective genotyping of the cohort results in small group sizes in the rarer E2 and E4 subgroups. Either larger well-standardised intervention trials or smaller trials with prospective recruitment according to apoE genotype are needed to fully establish the impact of diet on genotype-CVD associations and to establish the potential of dietary strategies such as reduced total fat, saturated fat, or increased antioxidant intakes to counteract the increased CVD burden in apoE4 carriers.

NOS3 gene polymorphisms are associated with risk markers of cardiovascular disease, and interact with omega-3 polyunsaturated fatty acids

Objective Omega-3 polyunsaturated fatty acids (n-3 PUFA) may protect against the development of cardiovascular disease (CVD). Genotype at key genes such as nitric oxide synthase (NOS3) may determine responsiveness to fatty acids. Gene–nutrient interactions may be important in modulating the development of CVD, particularly in high-risk individuals with the metabolic syndrome (MetS). Methods Biomarkers of CVD risk, plasma fatty acid composition, and NOS3 single nucleotide polymorphism (SNP) genotype (rs11771443, rs1800783, rs1800779, rs1799983, rs3918227, and rs743507) were determined in 450 individuals with the MetS from the LIPGENE dietary intervention cohort. The effect of dietary fat modification for 12 weeks on metabolic indices of the MetS was determined to understand potential NOS3 gene–nutrient interactions. Results Several markers of inflammation and dyslipidaemia were significantly different between the genotype groups. A significant gene–nutrient interaction was observed between the NOS3 rs1799983 SNP and plasma n-3 PUFA status on plasma triacylglycerol (TAG) concentrations. Minor allele carriers (AC + AA) showed an inverse association with significantly higher plasma TAG concentrations in those with low plasma n-3 PUFA status and vice versa but the major allele homozygotes (CC) did not. Following n-3 PUFA supplementation, plasma TAG concentrations of minor allele carriers of rs1799983 were considerably more responsive to changes in plasma n-3 PUFA, than major allele homozygotes. Conclusions Carriers of the minor allele at rs1799983 in NOS3 have plasma TAG concentrations which are more responsive to n-3 PUFA. This suggests that these individuals might show greater beneficial effects of n-3 PUFA consumption to reduce plasma TAG concentrations.

Detection of somaclonal variation during cocoa somatic embryogenesis characterised using cleaved amplified polymorphic sequence and the new freeware Artbio

The scarcity and stochastic nature of genetic mutations presents a significant challenge for scientists seeking to characterise de novo mutation frequency at specific loci. Such mutations can be particularly numerous during regeneration of plants from in vitro culture and can undermine the value of germplasm conservation efforts. We used cleaved amplified polymorphic sequence (CAPS) analysis to characterise new mutations amongst a clonal population of cocoa plants regenerated via a somatic embryogenesis protocol used previously for cocoa cryopreservation. Efficacy of the CAPS system for mutation detection was greatly improved after an ‘a priori’ in silico screen of reference target sequences for actual and potential restriction enzyme recognition sites using a new freely available software called Artbio. Artbio surveys known sequences for existing restriction enzyme recognition sites but also identifies all single nucleotide polymorphism (SNP) deviations from such motifs. Using this software, we performed an in silico screen of seven loci for restriction sites and their potential mutant SNP variants that were possible from 21 restriction enzymes. The four most informative locus-enzyme combinations were then used to survey the regenerant populations for de novo mutants. We characterised the pattern of point mutations and, using the outputs of Artbio, calculated the ratio of base substitution in 114 somatic embryo-derived cocoa regenerants originating from two explant genotypes. We found 49 polymorphisms, comprising 26.3% of the samples screened, with an inferred rate of 2.8 × 10−3 substitutions/screened base. This elevated rate is of a similar order of magnitude to previous reports of de novo microsatellite length mutations arising in the crop and suggests caution should be exercised when applying somatic embryogenesis for the conservation of plant germplasm.

Brodifacoum is effective against Norway Rats (Rattus norvegicus) in a tyrosine139cysteine focus of anticoagulant resistance, Westphalia, Germany

BACKGROUND: The single nucleotide polymorphism (SNP), and consequent amino acid exchange from tyrosine to cysteine at location 139 of the vkorc1 gene (i.e. tyrosine139cysteine or Y139C), is the most widespread anticoagulant resistance mutation in Norway rats (Rattus norvegicus Berk.) in Europe. Field trials were conducted to determine incidence of the Y139C SNP at two rat infested farms in Westphalia, Germany, and to estimate the practical efficacy against them of applications, using a pulsed baiting treatment regime, of a proprietary bait (KleratTM) containing 50 ppm brodifacoum. RESULTS: DNA analysis for the Y139C mutation showed that resistant rats were prevalent at the two farms, with an incidence of 80.0% and 78.6% respectively. Applications of brodifacoum bait achieved results of 99.2% and 100.0% control at the two farms, when measured by census baiting, although the treatment was somewhat prolonged at one site due to the abundance of attractive alternative food. CONCLUSION: The study showed that 50 ppm brodifacoum bait is fully effective against the Y139C SNP at the Münsterland focus and is likely to be so elsewhere in Europe where this mutation is found. The pulsed baiting regime reduced to relatively low levels the quantity of bait required to control these two substantial resistant Norway rat infestations. Previous studies had shown much larger quantities of bromadiolone and difenacoum baits used in ineffective treatments against Y139C resistant rats in the Münsterland. These results should be considered when making decisions about the use of anticoagulant against resistant Norway rats and their potential environmental impacts.

Isolation and characterization of an AINTEGUMENTA-like gene in different coconut (Cocos nucifera L.) varieties from Sri Lanka

Development of an efficient tissue culture protocol in coconut is hampered by numerous technical constraints. Thus a greater understanding of the fundamental aspects of embryogenesis is essential. The role of AINTEGUMENTA-like genes in embryogenesis has been elucidated not only in model plants but also in economically important crops. A coconut gene, CnANT, that encodes two APETALA2 (AP2) domains and a conserved linker region similar to those of the BABY BOOM transcription factor was cloned, characterized, and its tissue specific expression was examined. The full-length cDNA of 1,780 bp contains a 1,425-bp open reading frame that encodes a putative peptide of 474 amino acids. The genomic DNA sequence includes 2,317 bp and consists of nine exons interrupted by eight introns. The exon/intron organization of CnANT is similar to that of homologous genes in other plant species. Analysis of differential tissue expression by real-time polymerase chain reaction indicated that CnANT is expressed more highly in in vitro grown tissues than in other vegetative tissues. Sequence comparison of the genomic sequence of CnANT in different coconut varieties revealed one single nucleotide polymorphism and one indel in the first exon and first intron, respectively, which differentiate the Tall group of trees from Dwarfs. The indel sequence, which can be considered a simple sequence repeats marker, was successfully used to distinguish the Tall and Dwarf groups as well as to develop a marker system, which may be of value in the identification of parental varieties that are used in coconut breeding programs in Sri Lanka.

Genetic analysis reveals a novel locus in Vicia faba decoupling pigmentation in the flower from that in the extra-floral nectaries

The genetics of the stipule spot pigmentation (SSP) in faba bean (Vicia faba L.) was studied using four inbred lines, of which Disco/2 was zero-tannin (zt2) with colourless stipule spots, ILB938/2 was normal-tannin (ZT2) with colourless stipule spots, and both Aurora/2 and Mélodie/2 were ZT2 with coloured stipule spots. Crosses Mélodie/2 × ILB 938/2, Mélodie/2 × Disco/2, ILB 938/2 × Aurora/2 and ILB 938/2 × Disco/2 (A, B, C and D, respectively) were prepared, along with reciprocals and backcrosses, and advanced through single-seed descent. All F1 hybrid plants had pigmented stipule spots, and in the F2 generation, the segregation ratio fit 3 coloured:1 colourless in crosses A, B and C and 9:7 in cross D. In the F3 generation, the ratio fit 5:3 in crosses A and C and 25:39 in cross D, and in the F4 generation, 9:7 in cross A. SSP was linked to the zero-tannin characteristics (white flower) only in cross B. The results show that coloured stipule spot is dominant to colourless and that colouration is determined by two unlinked complementary recessive genes. We propose the symbols ssp2 for the gene associated with zt2 in Disco/2 and ssp1 for the gene not associated with tannin content in ILB938/2. The novel ssp1 locus was mapped at F5 in cross ‘A’ using Medicago truncatula-derived single-nucleotide polymorphism and was on chromosome 1 of faba bean, in a well-conserved region of M. truncatula chromosome 5 containing some candidate Myb and basic helix–loop–helix transcription factor genes.

An eight-parent multiparent advanced generation inter-cross population for winter-sown wheat: creation, properties, and validation

MAGIC populations represent one of a new generation of crop genetic mapping resources combining high genetic recombination and diversity. We describe the creation and validation of an eight-parent MAGIC population consisting of 1091 F7 lines of winter-sown wheat (Triticum aestivum L.). Analyses based on genotypes from a 90,000-single nucleotide polymorphism (SNP) array find the population to be well-suited as a platform for fine-mapping quantitative trait loci (QTL) and gene isolation. Patterns of linkage disequilibrium (LD) show the population to be highly recombined; genetic marker diversity among the founders was 74% of that captured in a larger set of 64 wheat varieties, and 54% of SNPs segregating among the 64 lines also segregated among the eight founder lines. In contrast, a commonly used reference bi-parental population had only 54% of the diversity of the 64 varieties with 27% of SNPs segregating. We demonstrate the potential of this MAGIC resource by identifying a highly diagnostic marker for the morphological character "awn presence/absence" and independently validate it in an association-mapping panel. These analyses show this large, diverse, and highly recombined MAGIC population to be a powerful resource for the genetic dissection of target traits in wheat, and it is well-placed to efficiently exploit ongoing advances in phenomics and genomics. Genetic marker and trait data, together with instructions for access to seed, are available at http://www.niab.com/MAGIC/.

Changes in physical activity following a genetic-based internet-delivered personalized intervention: randomized controlled trial (Food4Me)

Background: There is evidence that physical activity (PA) can attenuate the influence of the fat mass- and obesity-associated (FTO) genotype on the risk to develop obesity. However, whether providing personalized information on FTO genotype leads to changes in PA is unknown. Objective: The purpose of this study was to determine if disclosing FTO risk had an impact on change in PA following a 6-month intervention. Methods: The single nucleotide polymorphism (SNP) rs9939609 in the FTO gene was genotyped in 1279 participants of the Food4Me study, a four-arm, Web-based randomized controlled trial (RCT) in 7 European countries on the effects of personalized advice on nutrition and PA. PA was measured objectively using a TracmorD accelerometer and was self-reported using the Baecke questionnaire at baseline and 6 months. Differences in baseline PA variables between risk (AA and AT genotypes) and nonrisk (TT genotype) carriers were tested using multiple linear regression. Impact of FTO risk disclosure on PA change at 6 months was assessed among participants with inadequate PA, by including an interaction term in the model: disclosure (yes/no) × FTO risk (yes/no). Results: At baseline, data on PA were available for 874 and 405 participants with the risk and nonrisk FTO genotypes, respectively. There were no significant differences in objectively measured or self-reported baseline PA between risk and nonrisk carriers. A total of 807 (72.05%) of the participants out of 1120 in the personalized groups were encouraged to increase PA at baseline. Knowledge of FTO risk had no impact on PA in either risk or nonrisk carriers after the 6-month intervention. Attrition was higher in nonrisk participants for whom genotype was disclosed (P=.01) compared with their at-risk counterparts. Conclusions: No association between baseline PA and FTO risk genotype was observed. There was no added benefit of disclosing FTO risk on changes in PA in this personalized intervention. Further RCT studies are warranted to confirm whether disclosure of nonrisk genetic test results has adverse effects on engagement in behavior change.

p22phox C242T SNP inhibits inflammatory oxidative damage to endothelial cells and vessels

Background— T NADPH oxidase, by generating reactive oxygen species, is involved in the pathophysiology of many cardiovascular diseases and represents a therapeutic target for the development of novel drugs. A single-nucleotide polymorphism (SNP) C242T of the p22phox subunit of NADPH oxidase has been reported to be negatively associated with coronary heart disease (CHD) and may predict disease prevalence. However, the underlying mechanisms remain unknown. Methods and Results— Using computer molecular modelling we discovered that C242T SNP causes significant structural changes in the extracellular loop of p22phox and reduces its interaction stability with Nox2 subunit. Gene transfection of human pulmonary microvascular endothelial cells showed that C242T p22phox reduced significantly Nox2 expression but had no significant effect on basal endothelial O2.- production or the expression of Nox1 and Nox4. When cells were stimulated with TNFα (or high glucose), C242T p22phox inhibited significantly TNFα-induced Nox2 maturation, O2.- production, MAPK and NFκB activation and inflammation (all p<0.05). These C242T effects were further confirmed using p22phox shRNA engineered HeLa cells and Nox2-/- coronary microvascular endothelial cells. Clinical significance was investigated using saphenous vein segments from non CHD subjects after phlebectomies. TT (C242T) allele was common (prevalence of ~22%) and compared to CC, veins bearing TT allele had significantly lower levels of Nox2 expression and O2.- generation in response to high glucose challenge. Conclusions— C242T SNP causes p22phox structural changes that inhibit endothelial Nox2 activation and oxidative response to TNFα or high glucose stimulation. C242T SNP may represent a natural protective mechanism against inflammatory cardiovascular diseases.

Analysis of DNA polymorphism in ancient barley herbarium material: validation of the KASP SNP genotyping platform

The availability of crop specimens archived in herbaria and old seed collections represent valuable resources for the analysis of plant genetic diversity and crop domestication. The ability to extract ancient DNA (aDNA) from such samples has recently allowed molecular genetic investigations to be undertaken in ancient materials. While analyses of aDNA initially focused on the use of markers which occur in multiple copies such as the internal transcribed spacer region (ITS) within ribosomal DNA and those requiring amplification of short DNA regions of variable length such as simple sequence repeats (SSRs), emphasis is now moving towards the genotyping of single nucleotide polymorphisms (SNPs), traditionally undertaken in aDNA by Sanger sequencing. Here, using a panel of barley aDNA samples previously surveyed by Sanger sequencing for putative causative SNPs within the flowering-time gene PPD-H1, we assess the utility of the Kompetitive Allele Specific PCR (KASP) genotyping platform for aDNA analysis. We find KASP to out-perform Sanger sequencing in the genotyping of aDNA samples (78% versus 61% success, respectively), as well as being robust to contamination. The small template size (≥46 bp) and one-step, closed-tube amplification/genotyping process make this platform ideally suited to the genotypic analysis of aDNA, a process which is often hampered by template DNA degradation and sample cross-contamination. Such attributes, as well as its flexibility of use and relatively low cost, make KASP particularly relevant to the genetic analysis of aDNA samples. Furthermore, KASP provides a common platform for the genotyping and analysis of corresponding SNPs in ancient, landrace and modern plant materials. The extended haplotype analysis of PPD-H1 undertaken here (allelic variation at which is thought to be important for the spread of domestication and local adaptation) provides further resolution to the previously identified geographic cline of flowering-time allele distribution, illustrating how KASP can be used to aid genetic analyses of aDNA from plant species. We further demonstrate the utility of KASP by genotyping ten additional genetic markers diagnostic for morphological traits in barley, shedding light on the phenotypic traits, alleles and allele combinations present in these unviable ancient specimens, as well as their geographic distributions.