Paternal genetic variants and risk of obstructive heart defects: A parent-of-origin approach

Previous research on risk factors for obstructive heart defects (OHDs) focused on maternal and infant genetic variants, prenatal environmental exposures, and their potential interaction effects. Less is known about the role of paternal genetic variants or environmental exposures and risk of OHDs. We examined parent-of-origin effects in transmission of alleles in the folate, homocysteine, or transsulfuration pathway genes on OHD occurrence in offspring. We used data on 569 families of liveborn infants with OHDs born between October 1997 and August 2008 from the National Birth Defects Prevention Study to conduct a family-based case-only study. Maternal, paternal, and infant DNA were genotyped using an Illumina Golden Gate custom single nucleotide polymorphism (SNP) panel. Relative risks (RR), 95% confidence interval (CI), and likelihood ratio tests from log-linear models were used to estimate the parent-of-origin effect of 877 SNPs in 60 candidate genes in the folate, homocysteine, and transsulfuration pathways on the risk of OHDs. Bonferroni correction was applied for multiple testing. We identified 3 SNPs in the transsulfuration pathway and 1 SNP in the folate pathway that were statistically significant after Bonferroni correction. Among infants who inherited paternally-derived copies of the G allele for rs6812588 in the RFC1 gene, the G allele for rs1762430 in the MGMT gene, and the A allele for rs9296695 and rs4712023 in the GSTA3 gene, RRs for OHD were 0.11 (95% CI: 0.04, 0.29, P = 9.16x10^-7), 0.30 (95% CI: 0.17, 0.53, P = 9.80x10^-6), 0.34 (95% CI: 0.20, 0.57, P = 2.28x10^-5), and 0.34 (95% CI: 0.20, 0.58, P = 3.77x10^-5), respectively, compared to infants who inherited maternally-derived copies of the same alleles. We observed statistically significant decreased risk of OHDs among infants who inherited paternal gene variants involved in folate and transsulfuration pathways.     

Genome-wide association study identifies genetic loci associated with iron deficiency

The existence of multiple inherited disorders of iron metabolism in man, rodents and other vertebrates suggests genetic contributions to iron deficiency. To identify new genomic locations associated with iron deficiency, a genome-wide association study (GWAS) was performed using DNA collected from white men aged≥25 y and women≥50 y in the Hemochromatosis and Iron Overload Screening (HEIRS) Study with serum ferritin (SF)≤12 μg/L (cases) and iron replete controls (SF>100 μg/L in men, SF>50 μg/L in women). Regression analysis was used to examine the association between case-control status (336 cases, 343 controls) and quantitative serum iron measures and 331,060 single nucleotide polymorphism (SNP) genotypes, with replication analyses performed in a sample of 71 cases and 161 controls from a population of white male and female veterans screened at a US Veterans Affairs (VA) medical center. Five SNPs identified in the GWAS met genome-wide statistical significance for association with at least one iron measure, rs2698530 on chr. 2p14; rs3811647 on chr. 3q22, a known SNP in the transferrin (TF) gene region; rs1800562 on chr. 6p22, the C282Y mutation in the HFE gene; rs7787204 on chr. 7p21; and rs987710 on chr. 22q11 (GWAS observed P<1.51×10(-7) for all). An association between total iron binding capacity and SNP rs3811647 in the TF gene (GWAS observed P=7.0×10(-9), corrected P=0.012) was replicated within the VA samples (observed P=0.012). Associations with the C282Y mutation in the HFE gene also were replicated. The joint analysis of the HEIRS and VA samples revealed strong associations between rs2698530 on chr. 2p14 and iron status outcomes. These results confirm a previously-described TF polymorphism and implicate one potential new locus as a target for gene identification.     

Exogenous visual orienting is associated with specific neurotransmitter genetic markers: a population-based genetic association study

Background

Currently, there is a sense that the spatial orienting of attention is related to genotypic variations in cholinergic genes but not to variations in dopaminergic genes. However, reexamination of associations with both cholinergic and dopaminergic genes is warranted because previous studies used endogenous rather than exogenous cues and costs and benefits were not analyzed separately. Examining costs (increases in response time following an invalid pre-cue) and benefits (decreases in response time following a valid pre-cue) separately could be important if dopaminergic genes (implicated in disorders such as attention deficit disorder) independently influence the different processes of orienting (e.g., disengage, move, engage).

Methodology/Principal Findings

We tested normal subjects (N = 161) between 18 and 61 years. Participants completed a computer task in which pre-cues preceded the presence of a target. Subjects responded (with a key press) to the location of the target (right versus left of fixation). The cues could be valid (i.e., appear where the target would appear) or invalid (appear contralateral to where the target would appear). DNA sequencing assays were performed on buccal cells to genotype known genetic markers and these were examined for association with task scores. Here we show significant associations between visual orienting and genetic markers (on COMT, DAT1, and APOE; R²s from 4% to 9%).

Conclusions/Significance

One measure in particular – the response time cost of a single dim, invalid cue – was associated with dopaminergic markers on COMT and DAT1. Additionally, variations of APOE genotypes based on the ε2/ε3/ε4 alleles were also associated with response time differences produced by simultaneous cues with unequal luminances. We conclude that individual differences in visual orienting are related to several dopaminergic markers as well as to a cholinergic marker. These results challenge the view that orienting is not associated with genotypic variation in dopaminergic genes.     

Cystic fibrosis mutations and associated haplotypes in Bulgaria – a comparative population genetic study

We present data on the population genetics of cystic fibrosis (CF) in Bulgaria, obtained by comprehensive mutation analysis and the construction of intragenic microsatellite haplotypes. The sample of 262 CF alleles analysed is representative of the patients diagnosed during the period of referral and of the three main ethnic groups in the country. ΔF508 accounted for 100% of Gypsy CF alleles, which thus differed significantly from both Bulgarians and ethnic Turks. In Bulgarian and Turkish CF patients, 92% of the mutant alleles were identified, yielding a total of 25 different mutations, of which only 7 occurred at frequencies higher than 1%. The findings were compared to other European populations and to the distribution of phenylketonuria mutations. Genetic distances and population trees demonstrated that in the south-eastern tip of Europe, the overall distribution of CF mutations and polymorphic haplotypes is very close to that of Mediterranean populations, with a high frequency of N1303K and G542X, a large number of rare mutations and a prevalence of the 23 31 13 haplotype in association with ΔF508. These findings are consistent with a main role for the Neolithic expansion in the shaping of the CF mutation spectrum in Bulgaria and southern Europe. Received: 1 September 1996     

Polymorphisms and haplotypes in MyD88 are associated with the development of sarcoidosis: a candidate-gene association study

Sarcoidosis is considered as a disorder of protracted immune response to an as yet unidentified causative agent that leads to granuloma formation. Material from M. tuberculosis and P. acne has been repeatedly detected in the sarcoidosis lesions, implying the involvement of the Toll-like receptor2 (TLR2) gene that responds to these intracellular pathogens. Since TLR2 association studies have produced controversial results, we sought to investigate whether the downstream signalling molecule MyD88 could be linked to disease susceptibility. We analyzed a total of 93 cases with sarcoidosis and of 89 controls for the most common MyD88 SNPs: ?938C>A (rs4988453) and 1944C>G (rs4988457). There is evidence that the genotype distributions of both variants are associated with the development of sarcoidosis (p = 0.038 for ?938C>A and p = 0.026 for 1944C>G). In particular, ?938A and 1944G carriers were associated with risk of sarcoidosis [OR = 2.48 (1.23–5.02) and OR = 0.33 (0.14–0.76)], respectively, indicating dominance of the mutant alleles; however, the adjustment of the effect size for age and sex diminished the significance. The haplotype analysis showed association for the ?938A/1944G haplotype (p < 0.001). Since genetic association studies have linked MyD88 to Hodgkin’s lymphoma it is tempting to speculate that MyD88 may contribute to the granuloma formation that characterizes sarcoidosis.     

Genome-wide association study identifies four genetic loci associated with thyroid volume and goiter risk

Thyroid disorders such as goiters represent important diseases, especially in iodine-deficient areas. Sibling studies have demonstrated that genetic factors substantially contribute to the interindividual variation of thyroid volume. We performed a genome-wide association study of this phenotype by analyzing a discovery cohort consisting of 3620 participants of the Study of Health in Pomerania (SHIP). Four genetic loci were associated with thyroid volume on a genome-wide level of significance. Of these, two independent loci are located upstream of and within CAPZB, which encodes the β subunit of the barbed-end F-actin binding protein that modulates actin polymerization, a process crucial in the colloid engulfment during thyroglobulin mobilization in the thyroid. The third locus marks FGF7, which encodes fibroblast growth factor 7. Members of this protein family have been discussed as putative signal molecules involved in the regulation of thyroid development. The fourth locus represents a "gene desert" on chromosome 16q23, located directly downstream of the predicted coding sequence LOC440389, which, however, had already been removed from the NCBI database as a result of the standard genome annotation processing at the time that this study was initiated. Experimental proof of the formerly predicted mature mRNA, however, demonstrates that LOC440389 indeed represents a real gene. All four associations were replicated in an independent sample of 1290 participants of the KORA study. These results increase the knowledge about genetic factors and physiological mechanisms influencing thyroid volume.     

Genome-wide association study to identify genetic loci associated with gastrointestinal nematode resistance in Katahdin sheep

Resistance to gastrointestinal nematodes has previously been shown to be a moderately heritable trait in some breeds of sheep, but the mechanisms of resistance are not well understood. Selection for resistance currently relies upon faecal egg counts (FEC), blood packed cell volumes and FAMACHA visual indicator scores of anaemia. Identifying genomic markers associated with disease resistance would potentially improve the selection process and provide a more reliable means of classifying and understanding the biology behind resistant and susceptible sheep. A GWAS was conducted to identify possible genetic loci associated with resistance to Haemonchus contortus in Katahdin sheep. Forty animals were selected from the top and bottom 10% of estimated breeding values for FEC from a total pool of 641 sires and ram lambs. Samples were genotyped using Applied Biosystems™ Axiom™ Ovine Genotyping Array (50K) consisting of 51 572 SNPs. Following quality control, 46 268 SNPs were included in subsequent analyses. Analyses were conducted using a linear regression model in plink v1.90 and a single-locus mixed model in snp and variation suite . Genome-wide significance was determined by a Bonferroni correction for multiple testing. Using linear regression, loci on chromosomes 2, 3, 16, 23 and 24 were significantly associated at the genome level with FEC estimated breeding values, and we identified a region on chromosome 2 that was significant using both statistical analyses. We suggest a potential role for the gene DIS3L2 for gastrointestinal nematode resistance in Katahdin sheep, although further research is needed to validate these findings.     

Selection of genetic markers for association analyses,using linkage disequilibrium and haplotypes

The genotyping of closely spaced single-nucleotide polymorphism (SNP) markers frequently yields highly correlated data, owing to extensive linkage disequilibrium (LD) between markers. The extent of LD varies widely across the genome and drives the number of frequent haplotypes observed in small regions. Several studies have illustrated the possibility that LD or haplotype data could be used to select a subset of SNPs that optimize the information retained in a genomic region while reducing the genotyping effort and simplifying the analysis. We propose a method based on the spectral decomposition of the matrices of pairwise LD between markers, and we select markers on the basis of their contributions to the total genetic variation. We also modify Clayton's "haplotype tagging SNP" selection method, which utilizes haplotype information. For both methods, we propose sliding window-based algorithms that allow the methods to be applied to large chromosomal regions. Our procedures require genotype information about a small number of individuals for an initial set of SNPs and selection of an optimum subset of SNPs that could be efficiently genotyped on larger numbers of samples while retaining most of the genetic variation in samples. We identify suitable parameter combinations for the procedures, and we show that a sample size of 50-100 individuals achieves consistent results in studies of simulated data sets in linkage equilibrium and LD. When applied to experimental data sets, both procedures were similarly effective at reducing the genotyping requirement while maintaining the genetic information content throughout the regions. We also show that haplotype-association results that Hosking et al. obtained near CYP2D6 were almost identical before and after marker selection.     

Identification of genetic loci associated with growth traits at weaning in yak through a genome-wide association study

A multilocus GWAS was performed to explore the genetic architecture of four growth traits in yak. In total, 354 female yaks for which measurements of body weight (BW), withers height (WH), body length (BL) and chest girth (CG) at weaning were available underwent genotyping with the Illumina BovineHD BeadChip (770K). After quality control, we retained 98 688 SNPs and 354 animals for GWAS analysis. We identified seven, 18, seven and nine SNPs (corresponding to seven, 17, seven and eight candidate genes) associated with BW, WH, BL and CG at weaning respectively. Interestingly, most of these candidate genes were reported to be involved in growth-related processes such as muscle formation, lipid deposition, feed efficiency, carcass composition and development of the central and peripheral nervous system. Our results offer novel insight into the molecular architecture underpinning yak growth traits. Further functional analyses are thus warranted to explore the molecular mechanisms whereby these genes affect these traits of interest.     

Complement C7 deficiency: seven further molecular defects and their associated marker haplotypes

Seven further molecular bases of C7 deficiency are described. All these new molecular defects involve single-nucleotide events, deletions and substitutions, some of which alter splice sites, and others codons. They are distributed along the C7 gene, but predominantly towards the 3′ end. All were found in compound heterozygous individuals. The C6/C7 marker haplotypes associated with most C7 defects are tabulated. Received: 30 March 1998 / Accepted: 3 July 1998     

Regression-based association analysis with clustered haplotypes through use of genotypes

Haplotype-based association analysis has been recognized as a tool with high resolution and potentially great power for identifying modest etiological effects of genes. However, in practice, its efficacy has not been as successfully reproduced as expected in theory. One primary cause is that such analysis tends to require a large number of parameters to capture the abundant haplotype varieties, and many of those are expended on rare haplotypes for which studies would have insufficient power to detect association even if it existed. To concentrate statistical power on more-relevant inferences, in this study, we developed a regression-based approach using clustered haplotypes to assess haplotype-phenotype association. Specifically, we generalized the probabilistic clustering methods of Tzeng to the generalized linear model (GLM) framework established by Schaid et al. The proposed method uses unphased genotypes and incorporates both phase uncertainty and clustering uncertainty. Its GLM framework allows adjustment of covariates and can model qualitative and quantitative traits. It can also evaluate the overall haplotype association or the individual haplotype effects. We applied the proposed approach to study the association between hypertriglyceridemia and the apolipoprotein A5 gene. Through simulation studies, we assessed the performance of the proposed approach and demonstrate its validity and power in testing for haplotype-trait association.     

'Smoking genes': a genetic association study

Some controversy exists on the specific genetic variants that are associated with nicotine dependence and smoking-related phenotypes. The purpose of this study was to analyse the association of smoking status and smoking-related phenotypes (included nicotine dependence) with 17 candidate genetic variants: CYP2A6*1×2, CYP2A6*2 (1799T>A) [rs1801272], CYP2A6*9 (-48T>G) [rs28399433], CYP2A6*12, CYP2A13*2 (3375C>T) [rs8192789], CYP2A13*3 (7520C>G), CYP2A13*4 (579G>A), CYP2A13*7 (578C>T) [rs72552266], CYP2B6*4 (785A>G), CYP2B6*9 (516G>T), CHRNA3 546C>T [rs578776], CHRNA5 1192G>A [rs16969968], CNR1 3764C>G [rs6928499], DRD2-ANKK1 2137G>A (Taq1A) [rs1800497], 5HTT LPR, HTR2A -1438A>G [rs6311] and OPRM1 118A>G [rs1799971]. We studied the genotypes of the aforementioned polymorphisms in a cohort of Spanish smokers (cases, N = 126) and ethnically matched never smokers (controls, N = 80). The results showed significant between-group differences for CYP2A6*2 and CYP2A6*12 (both P<0.001). Compared with carriers of variant alleles, the odds ratio (OR) for being a non-smoker in individuals with the wild-type genotype of CYP2A6*12 and DRD2-ANKK1 2137G>A (Taq1A) polymorphisms was 3.60 (95%CI: 1.75, 7.44) and 2.63 (95%CI: 1.41, 4.89) respectively. Compared with the wild-type genotype, the OR for being a non-smoker in carriers of the minor CYP2A6*2 allele was 1.80 (95%CI: 1.24, 2.65). We found a significant genotype effect (all P≤0.017) for the following smoking-related phenotypes: (i) cigarettes smoked per day and CYP2A13*3; (ii) pack years smoked and CYP2A6*2, CYP2A6*1×2, CYP2A13*7, CYP2B6*4 and DRD2-ANKK1 2137G>A (Taq1A); (iii) nicotine dependence (assessed with the Fagestrom test) and CYP2A6*9. Overall, our results suggest that genetic variants potentially involved in nicotine metabolization (mainly, CYP2A6 polymorphisms) are those showing the strongest association with smoking-related phenotypes, as opposed to genetic variants influencing the brain effects of nicotine, e.g., through nicotinic acetylcholine (CHRNA5), serotoninergic (HTR2A), opioid (OPRM1) or cannabinoid receptors (CNR1).     

HaploBuild: an algorithm to construct non-contiguous associated haplotypes in family based genetic studies

We have created a program that searches densely genotyped regions for associated non-contiguous haplotypes using a standard family based haplotype association test. This program was designed to expand upon the 'sliding window' methodologies commonly used for haplotype construction by allowing the association of subsets of single nucleotide polymorphisms (SNPs) to drive the construction of the haplotype. This strategy permits HaploBuild to construct more biologically relevant haplotypes that are not constrained by arbitrary length and contiguous orientation. Availability: http://snp.bumc.bu.edu.     

A study of genetic association with electrophysiological measures related to alcoholism: GAW14 data

Recently, alcohol-related traits have been shown to have a genetic component. Here, we study the association of specific genetic measures in one of the three sets of electrophysiological measures in families with alcoholism distributed as part of the Genetic Analysis Workshop 14 data, the NTTH (non-target case of Visual Oddball experiment for 4 electrode placements) phenotypes: ntth1, ntth2, ntth3, and ntth4. We focused on the analysis of the 786 Affymetrix markers on chromosome 4. Our desire was to find at least a partial answer to the question of whether ntth1, ntth2, ntth3, and ntth4 are separately or jointly genetically controlled, so we studied the principal components that explain most of the covariation of the four quantitative traits. The first principal component, which explains 70% of the covariation, showed association but not genetic linkage to two markers: tsc0272102 and tsc0560854. On the other hand, ntth1 appeared to be the trait driving the variation in the second principal component, which showed association and genetic linkage at markers in four regions: tsc0045058, tsc1213381, tsc0055068, and tsc0051777 at map distances 53.26, 85.42, 89.31, and 172.86, respectively. These results show that the partial answer to our starting question for this brief analysis is that the NTTH phenotypes are not jointly genetically controlled. The component ntth1 displays marked genetic linkage.     

Economic implications of using Japanese Black sires carrying recessive genes associated with genetic defects

The objective of this study was to calculate cumulative discounted expressions (CDE) for Japanese Black sires carrying a single defective allele in a herd by applying the gene-flow method to investigate the expression pattern of homozygous recessive genotype and to evaluate the monetary loss of using these sires. A single biallelic locus was considered with A representing the dominant allele and a representing the recessive allele. The gene-flow method was modified to consider the fitness of homozygous recessive genotype. Input parameters representing a typical situation in a Japanese Black cattle herd were used to calculate the CDE and the loss of using carrier sires. The effects of initial allele frequency and fitness on the CDE were determined for Aa and AA sires. The CDE of Aa sires were larger than those of AA sires under all initial allele frequencies and fitness. The difference in the CDE between using Aa and AA sires was largest when fitness was 0 (lethal recessive condition). The differences in the loss between Aa and AA sires were larger with increasing initial allele frequencies in lethal recessive condition. Applying the method used in this study to defects reported in Japanese Black cattle and with a population size of 628 000, the difference in the loss between using Aa and AA sires was US$48 575 800 and US$74 418 000 in the case of Band-3 and Claudin-16 deficiencies, respectively. The approach used in this study could be applied to other genetic defects in different breeds and species.     

A genome-wide association study identifies novel alleles associated with hair color and skin pigmentation

We conducted a multi-stage genome-wide association study of natural hair color in more than 10,000 men and women of European ancestry from the United States and Australia. An initial analysis of 528,173 single nucleotide polymorphisms (SNPs) genotyped on 2,287 women identified IRF4 and SLC24A4 as loci highly associated with hair color, along with three other regions encompassing known pigmentation genes. We confirmed these associations in 7,028 individuals from three additional studies. Across these four studies, SLC24A4 rs12896399 and IRF4 rs12203592 showed strong associations with hair color, with p=6.0×10⁻⁶² and p=7.46×10⁻¹²⁷, respectively. The IRF4 SNP was also associated with skin color (p=6.2×10⁻¹⁴), eye color (p=6.1×10⁻¹³), and skin tanning response to sunlight (p=3.9×10⁻⁸⁹). A multivariable analysis pooling data from the initial GWAS and an additional 1,440 individuals suggested that the association between rs12203592 and hair color was independent of rs1540771, a SNP between the IRF4 and EXOC2 genes previously found to be associated with hair color. After adjustment for rs12203592, the association between rs1540771 and hair color was not significant (p=0.52). One variant in the MATP gene was associated with hair color. A variant in the HERC2 gene upstream of the OCA2 gene showed the strongest and independent association with hair color compared with other SNPs in this region, including three previously reported SNPs. The signals detected in a region around the MC1R gene were explained by MC1R red hair color alleles. Our results suggest that the IRF4 and SLC24A4 loci are associated with human hair color and skin pigmentation.     

A genome-wide association study identified new variants associated with mathematical abilities in Chinese children

Mathematical ability is moderately heritable, and it is a complex trait which can be evaluated in several different categories. A few genetic studies have been published on general mathematical ability. However, no genetic study focused on specific mathematical ability categories. In this study, we separately performed genome-wide association studies on 11 mathematical ability categories in 1146 students from Chinese elementary schools. We identified seven genome-wide significant single nucleotide polymorphisms (SNPs) with strong linkage disequilibrium among each other (all r² > 0.8) associated with mathematical reasoning ability (top SNP: rs34034296, p = 2.01 × 10⁻⁸, nearest gene: CUB and Sushi multiple domains 3, CSMD3). We replicated one SNP (rs133885) from 585 SNPs previously reported to be associated with general mathematical ability associated with division ability in our data (p = 1.053 × 10⁻⁵). In the gene- and gene-set enrichment analysis by MAGMA, we found three significant enrichments of associations with three mathematical ability categories for three genes (LINGO2, OAS1 and HECTD1). We also observed four significant enrichments of associations with four mathematical ability categories for three gene sets. Our results suggest new candidate genetic loci for the genetics of mathematical ability.