Case-parent triads: estimating single- and double-dose effects of fetal and maternal disease gene haplotypes

Efficient inference of parent-of-origin effect using case-control mother-child genotype data

Parent-of-origin effect plays an important role in mammal development and disorder. Case-control mother-child pair genotype data can be used to detect parent-of-origin effect and is often convenient to collect in practice. Most existing methods for assessing parent-of-origin effect do not incorporate any covariates, which may be required to control for confounding factors. We propose to model the parent-of-origin effect through a logistic regression model, with predictors including maternal and child genotypes, parental origins, and covariates. The parental origins may not be fully inferred from genotypes of a target genetic marker, so we propose to use genotypes of markers tightly linked to the target marker to increase inference efficiency. A robust statistical inference procedure is developed based on a modified profile log-likelihood in a retrospective way. A computationally feasible expectation-maximization algorithm is devised to estimate all unknown parameters involved in the modified profile log-likelihood. This algorithm differs from the conventional expectation-maximization algorithm in the sense that it is based on a modified instead of the original profile log-likelihood function. The convergence of the algorithm is established under some mild regularity conditions. This expectation-maximization algorithm also allows convenient handling of missing child genotypes. Large sample properties, including weak consistency, asymptotic normality, and asymptotic efficiency, are established for the proposed estimator under some mild regularity conditions. Finite sample properties are evaluated through extensive simulation studies and the application to a real dataset.

Statistical methods to detect mother-father genetic interaction effects on risk of infertility: A genome-wide approach

Infertility is a heterogeneous phenotype, and for many couples, the causes of fertility problems remain unknown. One understudied hypothesis is that allelic interactions between the genotypes of the two parents may influence the risk of infertility. Our aim was, therefore, to investigate how allelic interactions can be modeled using parental genotype data linked to 15,789 pregnancies selected from the Norwegian Mother, Father, and Child Cohort Study. The newborns in 1304 of these pregnancies were conceived using assisted reproductive technologies (ART), and the remainder were conceived naturally. Treating the use of ART as a proxy for infertility, different parameterizations were implemented in a genome-wide screen for interaction effects between maternal and paternal alleles at the same locus. Some of the models were more similar in the way they were parameterized, and some produced similar results when implemented on a genome-wide scale. The results showed near-significant interaction effects in genes relevant to the phenotype under study, such as Dynein axonemal heavy chain 17 (DNAH17) with a recognized role in male infertility. More generally, the interaction models presented here are readily adaptable to the study of other phenotypes in which maternal and paternal allelic interactions are likely to be involved.

Benchmarking statistical methods for analyzing parent-child dyads in genetic association studies

Genetic association studies of child health outcomes often employ family-based study designs. One of the most popular family-based designs is the case-parent trio design that considers the smallest possible nuclear family consisting of two parents and their affected child. This trio design is particularly advantageous for studying relatively rare disorders because it is less prone to type 1 error inflation due to population stratification compared to population-based study designs (e.g., case-control studies). However, obtaining genetic data from both parents is difficult, from a practical perspective, and many large studies predominantly measure genetic variants in mother-child dyads. While some statistical methods for analyzing parent-child dyad data (most commonly involving mother-child pairs) exist, it is not clear if they provide the same advantage as trio methods in protecting against population stratification, or if a specific dyad design (e.g., case-mother dyads vs. case-mother/control-mother dyads) is more advantageous. In this article, we review existing statistical methods for analyzing genome-wide marker data on dyads and perform extensive simulation experiments to benchmark their type I errors and statistical power under different scenarios. We extend our evaluation to existing methods for analyzing a combination of case-parent trios and dyads together. We apply these methods on genotyped and imputed data from multiethnic mother-child pairs only, case-parent trios only or combinations of both dyads and trios from the Gene, Environment Association Studies consortium (GENEVA), where each family was ascertained through a child affected by nonsyndromic cleft lip with or without cleft palate. Results from the GENEVA study corroborate the findings from our simulation experiments. Finally, we provide recommendations for using statistical genetic association methods for dyads.

Methylome-wide association study of early life stressors and adult mental health

The environment and events that we are exposed to in utero, during birth and in early childhood influence our future physical and mental health. The underlying mechanisms that lead to these outcomes are unclear, but long-term changes in epigenetic marks, such as DNA methylation, could act as a mediating factor or biomarker. DNA methylation data were assayed at 713 522 CpG sites from 9537 participants of the Generation Scotland: Scottish Family Health Study, a family-based cohort with extensive genetic, medical, family history and lifestyle information. Methylome-wide association studies of eight early life environment phenotypes and two adult mental health phenotypes (major depressive disorder and brief resilience scale) were conducted using DNA methylation data collected from adult whole blood samples. Two genes involved with different developmental pathways (PRICKLE2, Prickle Planar Cell Polarity Protein 2 and ABI1, Abl-Interactor-1) were annotated to CpG sites associated with preterm birth (P < 1.27 × 10⁻⁹). A further two genes important to the development of sensory pathways (SOBP, Sine Oculis Binding Protein Homolog and RPGRIP1, Retinitis Pigmentosa GTPase Regulator Interacting Protein) were annotated to sites associated with low birth weight (P < 4.35 × 10⁻⁸). The examination of methylation profile scores and genes and gene-sets annotated from associated CpGs sites found no evidence of overlap between the early life environment and mental health conditions. Birth date was associated with a significant difference in estimated lymphocyte and neutrophil counts. Previous studies have shown that early life environments influence the risk of developing mental health disorders later in life; however, this study found no evidence that this is mediated by stable changes to the methylome detectable in peripheral blood.

Covariate adjusted inference of parent-of-origin effects using case-control mother-child paired multilocus genotype data

It is of great interest to identify parent-of-origin effects (POEs) since POEs play an important role in many human heritable disorders and human early life growth and development. POE is sometimes referred to as imprinting effect in the literature. Compared with the standard logistic regression analyses, retrospective likelihood-based statistical methods are more powerful in identifying POEs when data are collected from related individuals retrospectively. However, none of existing retrospective-based methods can appropriately incorporate covariates that should be adjusted for if they are confounding factors. In this paper, a novel semiparametric statistical method, M-HAP, is developed to detect POEs by fully exploring available information from multilocus genotypes of case-control mother-child pairs and covariates. Some large sample properties are established for M-HAP. Finite sample properties of M-HAP are illustrated by extensive simulation studies and real data applications to the Jerusalem Perinatal Study and the Danish National Birth Cohort study, which confirm the desired superiority of M-HAP over some existing methods. M-HAP has been implemented in the updated R package CCMO.

Aquaporin-1 and aquaporin-9 gene variations in sudden infant death syndrome

Several studies have indicated that a vulnerability in the development and regulation of brain function is involved in sudden infant death syndrome (SIDS). The aim of this study was to investigate the genes encoding the brain aquaporins (AQPs) AQP1 and AQP9 in SIDS. The hypothesis was that specific variants of these genes are part of the genetic vulnerability predisposing infants to sudden unexpected death. The study included 168 SIDS cases with a median age of 15.5 (range 2–52) weeks and 372 adolescent/adult deceased controls with a median age of 44 (range 11–91) years. In the AQP1 gene, the rs17159702 CC/CT genotypes were found to be associated with SIDS (p = 0.02). In the AQP9 gene, the combination of a TT genotype of rs8042354, rs2292711 and rs13329178 was more frequent in SIDS cases than in controls (p = 0.03). In the SIDS group, an association was found between genetic variations in the AQP1 gene and maternal smoking and between the 3xTT combination in the AQP9 gene and being found lifeless in a prone position. In conclusion, this study adds further evidence to the involvement of brain aquaporins in SIDS, suggesting that specific variants of AQP genes constitute a genetic predisposition, making the infant vulnerable to sudden death together with external risk factors and probably other genetic factors.

Altered DNA methylation profiles in blood from patients with sporadic Creutzfeldt-Jakob disease

Prion diseases are fatal and transmissible neurodegenerative disorders caused by the misfolding and aggregation of prion protein. Although recent studies have implicated epigenetic variation in common neurodegenerative disorders, no study has yet explored their role in human prion diseases. Here we profiled genome-wide blood DNA methylation in the most common human prion disease, sporadic Creutzfeldt–Jakob disease (sCJD). Our case–control study (n = 219), when accounting for differences in cell type composition between individuals, identified 38 probes at genome-wide significance (p < 1.24 × 10^–7). Nine of these sites were taken forward in a replication study, performed in an independent case–control (n = 186) cohort using pyrosequencing. Sites in or close to FKBP5, AIM2 (2 probes), UHRF1, KCNAB2 successfully replicated. The blood-based DNA methylation signal was tissue- and disease-specific, in that the replicated probe signals were unchanged in case–control studies using sCJD frontal-cortex (n = 84), blood samples from patients with Alzheimer’s disease, and from inherited and acquired prion diseases. Machine learning algorithms using blood DNA methylation array profiles accurately distinguished sCJD patients and controls. Finally, we identified sites whose methylation levels associated with prolonged survival in sCJD patients. Altogether, this study has identified a peripheral DNA methylation signature of sCJD with a variety of potential biomarker applications.

Case-Parent Trio Studies in Cleft Lip and Palate

Cleft lip with or without cleft palate (CL/P) is one of the most common congenital malformations in humans involving various genetic and environmental risk factors. The prevalence of CL/P varies according to geographical location, ethnicity, race, gender, and socioeconomic status, affecting approximately 1 in 800 live births worldwide. Genetic studies aim to understand the mechanisms contributory to a phenotype by measuring the association between genetic variants and also between genetic variants and phenotype population. Genome-wide association studies are standard tools used to discover genetic loci related to a trait of interest. Genetic association studies are generally divided into two main design types: population-based studies and family-based studies. The epidemiological population-based studies comprise unrelated individuals that directly compare the frequency of genetic variants between (usually independent) cases and controls. The alternative to population-based studies (case-control designs) includes various family-based study designs that comprise related individuals. An example of such a study is a case-parent trio design study, which is commonly employed in genetics to identify the variants underlying complex human disease where transmission of alleles from parents to offspring is studied. This article describes the fundamentals of case-parent trio study, trio design and its significances, statistical methods, and limitations of the trio studies.

Gene-methylation interactions: discovering region-wise DNA methylation levels that modify SNP-associated disease risk

BACKGROUND

Current technology allows rapid assessment of DNA sequences and methylation levels at a single-site resolution for hundreds of thousands of sites in the human genome, in thousands of individuals simultaneously. This has led to an increase in epigenome-wide association studies (EWAS) of complex traits, particularly those that are poorly explained by previous genome-wide association studies (GWAS). However, the genome and epigenome are intertwined, e.g., DNA methylation is known to affect gene expression through, for example, genomic imprinting. There is thus a need to go beyond single-omics data analyses and develop interaction models that allow a meaningful combination of information from EWAS and GWAS.

RESULTS

We present two new methods for genetic association analyses that treat offspring DNA methylation levels as environmental exposure. Our approach searches for statistical interactions between SNP alleles and DNA methylation (G ×Me) and between parent-of-origin effects and DNA methylation (PoO ×Me), using case-parent triads or dyads. We use summarized methylation levels over nearby genomic region to ease biological interpretation. The methods were tested on a dataset of parent-offspring dyads, with EWAS data on the offspring. Our results showed that methylation levels around a SNP can significantly alter the estimated relative risk. Moreover, we show how a control dataset can identify false positives.

CONCLUSIONS

The new methods, G ×Me and PoO ×Me, integrate DNA methylation in the assessment of genetic relative risks and thus enable a more comprehensive biological interpretation of genome-wide scans. Moreover, our strategy of condensing DNA methylation levels within regions helps overcome specific disadvantages of using sparse chip-based measurements. The methods are implemented in the freely available R package Haplin ( https://cran.r-project.org/package=Haplin ), enabling fast scans of multi-omics datasets.

Design efficiency in genetic association studies

Selecting the best design for genetic association studies requires careful deliberation; different study designs can be used to scan for different genetic effects, and each design has its own set of strengths and limitations. A variety of family and unrelated control configurations are amenable to genetic association analyses, including the case-control design, case-parent triads, and case-parent triads in combination with unrelated controls or control-parent triads. Ultimately, the goal is to choose the design that achieves the highest statistical power using the lowest cost. For given parameter values and genotyped individuals, designs can be compared directly by computing the power. However, a more informative and general design comparison can be achieved by studying the relative efficiency, defined as the ratio of variances of two different parameter estimators, corresponding to two separate designs. Using log-linear modeling, we derive the relative efficiency from the asymptotic variance of the parameter estimators and relate it to the concept of Pitman efficiency. The relative efficiency takes into account the fact that different designs impose different costs relative to the number of genotyped individuals. We show that while optimal efficiency for analyses of regular autosomal effects is achieved using the standard case-control design, the case-parent triad design without unrelated controls is efficient when searching for parent-of-origin effects. Due to the potential loss of efficiency, maternal genes should generally not be adjusted for in an initial genome-wide association study scan of offspring genes but instead checked post hoc. The relative efficiency calculations are implemented in our R package Haplin.

A genome-wide scan of cleft lip triads identifies parent-of-origin interaction effects between ANK3 and maternal smoking, and between ARHGEF10 and alcohol consumption

Background: Although both genetic and environmental factors have been reported to influence the risk of isolated cleft lip with or without cleft palate (CL/P), the exact mechanisms behind CL/P are still largely unaccounted for. We recently developed new methods to identify parent-of-origin (PoO) interactions with environmental exposures (PoOxE) and now apply them to data from a genome-wide association study (GWAS) of families with children born with isolated CL/P. Methods: Genotypes from 1594 complete triads and 314 dyads (1908 nuclear families in total) with CL/P were available for the current analyses. Of these families, 1024 were Asian, 825 were European and 59 had other ancestries. After quality control, 341,191 SNPs remained from the original 569,244. The exposures were maternal cigarette smoking, use of alcohol, and use of vitamin supplements in the periconceptional period. Our new methodology detects if PoO effects are different across environmental strata and is implemented in the R-package Haplin. Results: Among Europeans, there was evidence of a PoOxSmoke effect for ANK3 with three SNPs (rs3793861, q=0.20, p=2.6e-6; rs7087489, q=0.20, p=3.1e-6; rs4310561, q=0.67, p=4.0e-5) and a PoOxAlcohol effect for ARHGEF10 with two SNPs (rs2294035, q=0.32, p=2.9e-6; rs4876274, q=0.76, p=1.3e-5). Conclusion: Our results indicate that the detected PoOxE effects have a plausible biological basis, and thus warrant replication in other independent cleft samples. Our demonstration of the feasibility of identifying complex interactions between relevant environmental exposures and PoO effects offers new avenues for future research aimed at unravelling the complex etiology of cleft lip defects.

A genome-wide scan of cleft lip triads identifies parent-of-origin interaction effects between ANK3 and maternal smoking, and between ARHGEF10 and alcohol consumption

Background: Although both genetic and environmental factors have been reported to influence the risk of isolated cleft lip with or without cleft palate (CL/P), the exact mechanisms behind CL/P are still largely unaccounted for. We recently developed new methods to identify parent-of-origin (PoO) interactions with environmental exposures (PoOxE) and applied them to families with children born with isolated cleft palate only. Here, we used the same genome-wide association study (GWAS) dataset and methodology to screen for PoOxE effects in the larger sample of CL/P triads. Methods: Genotypes from 1594 complete triads and 314 dyads (1908 nuclear families in total) with CL/P were available for the current analyses. Of these families, 1024 were Asian, 825 were European and 59 had other ancestries. After quality control, 341,191 SNPs remained from the original 569,244. The exposures were maternal cigarette smoking, use of alcohol, and use of vitamin supplements in the periconceptional period. The methodology applied in the analyses is implemented in the R-package Haplin. Results: Among Europeans, there was evidence of a PoOxSmoke effect for ANK3 with three SNPs (rs3793861, q=0.20, p=2.6e-6; rs7087489, q=0.20, p=3.1e-6; rs4310561, q=0.67, p=4.0e-5) and a PoOxAlcohol effect for ARHGEF10 with two SNPs (rs2294035, q=0.32, p=2.9e-6; rs4876274, q=0.76, p=1.3e-5). Conclusion: Our results indicate that the detected PoOxE effects have a plausible biological basis, and thus warrant replication in other independent cleft samples. Our demonstration of the feasibility of identifying complex interactions between relevant environmental exposures and PoO effects offers new avenues for future research aimed at unravelling  the complex etiology of cleft lip defects.

Haplin power analysis: a software module for power and sample size calculations in genetic association analyses of family triads and unrelated controls

Background

Log-linear and multinomial modeling offer a flexible framework for genetic association analyses of offspring (child), parent-of-origin and maternal effects, based on genotype data from a variety of child-parent configurations. Although the calculation of statistical power or sample size is an important first step in the planning of any scientific study, there is currently a lack of software for genetic power calculations in family-based study designs. Here, we address this shortcoming through new implementations of power calculations in the R package Haplin, which is a flexible and robust software for genetic epidemiological analyses. Power calculations in Haplin can be performed analytically using the asymptotic variance-covariance structure of the parameter estimator, or else by a straightforward simulation approach. Haplin performs power calculations for child, parent-of-origin and maternal effects, as well as for gene-environment interactions. The power can be calculated for both single SNPs and haplotypes, either autosomal or X-linked. Moreover, Haplin enables power calculations for different child-parent configurations, including (but not limited to) case-parent triads, case-mother dyads, and case-parent triads in combination with unrelated control-parent triads.

Results

We compared the asymptotic power approximations to the power of analysis attained with Haplin. For external validation, the results were further compared to the power of analysis attained by the EMIM software using data simulations from Haplin. Consistency observed between Haplin and EMIM across various genetic scenarios confirms the computational accuracy of the inference methods used in both programs. The results also demonstrate that power calculations in Haplin are applicable to genetic association studies using either log-linear or multinomial modeling approaches.

Conclusions

Haplin provides a robust and reliable framework for power calculations in genetic association analyses for a wide range of genetic effects and etiologic scenarios, based on genotype data from a variety of child-parent configurations.

Electronic supplementary material

The online version of this article (10.1186/s12859-019-2727-3) contains supplementary material, which is available to authorized users.

Comparison of different functional prediction scores using a gene-based permutation model for identifying cancer driver genes

Background

Identifying cancer driver genes (CDG) is a crucial step in cancer genomic toward the advancement of precision medicine. However, driver gene discovery is a very challenging task because we are not only dealing with huge amount of data; but we are also faced with the complexity of the disease including the heterogeneity of background somatic mutation rate in each cancer patient. It is generally accepted that CDG harbor variants conferring growth advantage in the malignant cell and they are positively selected, which are critical to cancer development; whereas, non-driver genes harbor random mutations with no functional consequence on cancer. Based on this fact, function prediction based approaches for identifying CDG have been proposed to interrogate the distribution of functional predictions among mutations in cancer genomes (eLS 1–16, 2016). Assuming most of the observed mutations are passenger mutations and given the quantitative predictions for the functional impact of the mutations, genes enriched of functional or deleterious mutations are more likely to be drivers. The promises of these methods have been continually refined and can therefore be applied to increase accuracy in detecting new candidate CDGs. However, current function prediction based approaches only focus on coding mutations and lack a systematic way to pick the best mutation deleteriousness prediction algorithms for usage.

Results

In this study, we propose a new function prediction based approach to discover CDGs through a gene-based permutation approach. Our method not only covers both coding and non-coding regions of the genes; but it also accounts for the heterogeneous mutational context in cohort of cancer patients. The permutation model was implemented independently using seven popular deleteriousness prediction scores covering splicing regions (SPIDEX), coding regions (MetaLR, and VEST3) and pan-genome (CADD, DANN, Fathmm-MKL coding and Fathmm-MKL noncoding). We applied this new approach to somatic single nucleotide variants (SNVs) from whole-genome sequences of 119 breast and 24 lung cancer patients and compared the seven deleteriousness prediction scores for their performance in this study.

Conclusion

The new function prediction based approach not only predicted known cancer genes listed in the Cancer Gene Census (CGC), but also new candidate CDGs that are worth further investigation. The results showed the advantage of utilizing pan-genome deleteriousness prediction scores in function prediction based methods. Although VEST3 score, a deleteriousness prediction score for missense mutations, has the best performance in breast cancer, it was topped by CADD and Fathmm-MKL coding, two pan-genome deleteriousness prediction scores, in lung cancer.

Electronic supplementary material

The online version of this article (10.1186/s12920-018-0452-9) contains supplementary material, which is available to authorized users.

A Genome-Wide Search for Gene-Environment Effects in Isolated Cleft Lip with or without Cleft Palate Triads Points to an Interaction between Maternal Periconceptional Vitamin Use and Variants in ESRRG

Background: It is widely accepted that cleft lip with or without cleft palate (CL/P) results from the complex interplay between multiple genetic and environmental factors. However, a robust investigation of these gene-environment (GxE) interactions at a genome-wide level is still lacking for isolated CL/P. Materials and Methods: We used our R-package Haplin to perform a genome-wide search for GxE effects in isolated CL/P. From a previously published GWAS, genotypes and information on maternal periconceptional cigarette smoking, alcohol intake, and vitamin use were available on 1908 isolated CL/P triads of predominantly European or Asian ancestry. A GxE effect is present if the relative risk estimates for gene-effects in the offspring are different across exposure strata. We tested this using the relative risk ratio (RRR). Besides analyzing all ethnicities combined ("pooled analysis"), separate analyses were conducted on Europeans and Asians to investigate ethnicity-specific effects. To control for multiple testing, q-values were calculated from the p-values. Results: We identified significant GxVitamin interactions with three SNPs in "Estrogen-related receptor gamma" (ESRRG) in the pooled analysis. The RRRs (95% confidence intervals) were 0.56 (0.45-0.69) with rs1339221 (q = 0.011), 0.57 (0.46-0.70) with rs11117745 (q = 0.011), and 0.62 (0.50-0.76) with rs2099557 (q = 0.037). The associations were stronger when these SNPs were analyzed as haplotypes composed of two-SNP and three-SNP combinations. The strongest effect was with the "t-t-t" haplotype of the rs1339221-rs11117745-rs2099557 combination [RRR = 0.50 (0.40-0.64)], suggesting that the effects observed with the other SNP combinations, including those in the single-SNP analyses, were mainly driven by this haplotype. Although there were potential GxVitamin effects with rs17734557 and rs1316471 and GxAlcohol effects with rs9653456 and rs921876 in the European sample, respectively, none of the SNPs was located in or near genes with strong links to orofacial clefts. GxAlcohol and GxSmoke effects were not assessed in the Asian sample because of a lack of observations for these exposures. Discussion/Conclusion: We identified significant interactions between vitamin use and variants in ESRRG in the pooled analysis. These GxE effects are novel and warrant further investigations to elucidate their roles in orofacial clefting. If validated, they could provide prospects for exploring the impact of estrogens and vitamins on clefting, with potential translational applications.

Analysis of Parent-of-Origin Effects on the X Chromosome in Asian and European Orofacial Cleft Triads Identifies Associations with DMD, FGF13, EGFL6, and Additional Loci at Xp22.2

Background: Although both the mother's and father's alleles are present in the offspring, they may not operate at the same level. These parent-of-origin (PoO) effects have not yet been explored on the X chromosome, which motivated us to develop new methods for detecting such effects. Orofacial clefts (OFCs) exhibit sex-specific differences in prevalence and are examples of traits where a search for various types of effects on the X chromosome might be relevant. Materials and Methods: We upgraded our R-package Haplin to enable genome-wide analyses of PoO effects, as well as power simulations for different statistical models. 14,486 X-chromosome SNPs in 1,291 Asian and 1,118 European case-parent triads of isolated OFCs were available from a previous GWAS. For each ethnicity, cleft lip with or without cleft palate (CL/P) and cleft palate only (CPO) were analyzed separately using two X-inactivation models and a sliding-window approach to haplotype analysis. In addition, we performed analyses restricted to female offspring. Results: Associations were identified in "Dystrophin" (DMD, Xp21.2-p21.1), "Fibroblast growth factor 13" (FGF13, Xq26.3-q27.1) and "EGF-like domain multiple 6" (EGFL6, Xp22.2), with biologically plausible links to OFCs. Unlike EGFL6, the other associations on chromosomal region Xp22.2 had no apparent connections to OFCs. However, the Xp22.2 region itself is of potential interest because it contains genes for clefting syndromes [for example, "Oral-facial-digital syndrome 1" (OFD1) and "Midline 1" (MID1)]. Overall, the identified associations were highly specific for ethnicity, cleft subtype and X-inactivation model, except for DMD in which associations were identified in both CPO and CL/P, in the model with X-inactivation and in Europeans only. Discussion/Conclusion: The specificity of the associations for ethnicity, cleft subtype and X-inactivation model underscores the utility of conducting subanalyses, despite the ensuing need to adjust for additional multiple testing. Further investigations are needed to confirm the associations with DMD, EGF16, and FGF13. Furthermore, chromosomal region Xp22.2 appears to be a hotspot for genes implicated in clefting syndromes and thus constitutes an exciting direction to pursue in future OFCs research. More generally, the new methods presented here are readily adaptable to the study of X-linked PoO effects in other outcomes that use a family-based design.

Influence of IL15 gene variations on the clinical features, treatment response and risk of developing childhood acute lymphoblastic leukemia in Latvian population

Acute lymphoblastic leukemia (ALL) is the most common childhood malignancy. Modern treatment protocols allow achievement of long-term event-free survival rates in up to 85% of cases, although the treatment response varies among different patient groups. It is hypothesized that treatment response is influenced by the IL15 gene variations, although research results are conflicting. To analyze IL15 gene variations influence treatment response, clinical course and the risk of developing ALL we performed a case-control and family-based study. The study included 81 patients with childhood ALL. DNA samples of both or one biological parent were available for 62 of ALL patients and 130 age and gender adjusted healthy samples were used as a control group. Analyzed IL15 gene variations: rs10519612, rs10519613 and rs17007695 were genotyped using PCR-RFLP assay. Our results shows that IL15 gene variations haplotypes are associated with the risk of developing childhood ALL (p < 0.05), although there is no such association for the variations separately. The variations rs10519612 and rs1059613 in a recessive pattern of inheritance were associated with hyperdiploidy (p = 0.048). Analyzed genetic variations had no impact on other clinical features and treatment response (assessed by the minimal residual disease) in our study.

Association of ARID5B Genetic Variants with Risk of Childhood B Cell Precursor Acute Lymphoblastic Leukaemia in Latvia

Background: Acute lymphoblastic leukaemia (ALL) is the most common malignancy in childhood. Despite numerous investigations very little is still known about its aetiology. However, in one genome wide association study conducted to identify the possible genetic risk factors, two allelic variations rs10821936 and rs10994982 in the 3rd intron of the ARID5B gene were identified as possible ALL risk alleles. Association between ARID5B gene variants and ALL risk was also been confirmed for different ethnic groups. Materials and Methods: Eight genetic variants in the gene ARID5B were genotyped - rs10994982, rs7908445, rs7923074, rs10821936, rs10821937, rs7896246, rs10821938 and rs7089424 in 77 ALL patients in remission and in 122 age and gender matched controls; parental samples were also genotyped in 50 cases. Results: Six out of the eight (rs7908445, rs7923074, rs10821936, rs10821937, rs7896246 and rs7089424) analysed allelic variations were identified in the case-control analysis as statistically significant risk alleles for ALL development. In the family study and using hybrid analysis, all allelic variations were significantly associated with ALL. During the study, risk haplotype was identified rs10994982/rs7908445/rs7923074/ rs10821936/ rs10821937/rs7896246/rs10821938/rs7089424 – ATACCAAG – with a frequency in cases of 0.17 and in the control group at 0.29 (chi square = 6.69, p value = 0.009). In the family association study the same haplotype showed statistical significance (chi squared = 10.3, p value = 0.001). Conclusions: Results of the study replicate and extend previously published findings for ARID5B localized allelic variants, but do not explain the mechanism of action related to the pathogenesis of ALL.

Parent-of-origin-environment interactions in case-parent triads with or without independent controls

With case–parent triad data, one can frequently deduce parent of origin of the child's alleles. This allows a parent‐of‐origin (PoO) effect to be estimated as the ratio of relative risks associated with the alleles inherited from the mother and the father, respectively. A possible cause of PoO effects is DNA methylation, leading to genomic imprinting. Because environmental exposures may influence methylation patterns, gene–environment interaction studies should be extended to allow for interactions between PoO effects and environmental exposures (i.e., PoOxE). One should thus search for loci where the environmental exposure modifies the PoO effect.

We have developed an extensive framework to analyze PoOxE effects in genome‐wide association studies (GWAS), based on complete or incomplete case–parent triads with or without independent control triads. The interaction approach is based on analyzing triads in each exposure stratum using maximum likelihood estimation in a log‐linear model. Interactions are then tested applying a Wald‐based posttest of parameters across strata. Our framework includes a complete setup for power calculations. We have implemented the models in the R software package Haplin.

To illustrate our PoOxE test, we applied the new methodology to top hits from our previous GWAS, assessing whether smoking during the periconceptional period modifies PoO effects on cleft palate only.

Genome-wide analysis of parent-of-origin interaction effects with environmental exposure (PoOxE): An application to European and Asian cleft palate trios

Cleft palate only is a common birth defect with high heritability. Only a small fraction of this heritability is explained by the genetic variants identified so far, underscoring the need to investigate other disease mechanisms, such as gene-environment (GxE) interactions and parent-of-origin (PoO) effects. Furthermore, PoO effects may vary across exposure levels (PoOxE effects). Such variation is the focus of this study. We upgraded the R-package Haplin to enable direct tests of PoOxE effects at the genome-wide level. From a previous GWAS, we had genotypes for 550 case-parent trios, of mainly European and Asian ancestry, and data on three maternal exposures (smoking, alcohol, and vitamins). Data were analyzed for Europeans and Asians separately, and also for all ethnicities combined. To account for multiple testing, a false discovery rate method was used, where q-values were generated from the p-values. In the Europeans-only analyses, interactions with maternal smoking yielded the lowest q-values. Two SNPs in the 'Interactor of little elongation complex ELL subunit 1' (ICE1) gene had a q-value of 0.14, and five of the 20 most significant SNPs were in the 'N-acetylated alpha-linked acidic dipeptidase-like 2' (NAALADL2) gene. No evidence of PoOxE effects was found in the other analyses. The connections to ICE1 and NAALADL2 are novel and warrant further investigation. More generally, the new methodology presented here is easily applicable to other traits and exposures in which a family-based study design has been implemented.

A new approach to chromosome-wide analysis of X-linked markers identifies new associations in Asian and European case-parent triads of orofacial clefts

Background

GWAS discoveries on the X-chromosome are underrepresented in the literature primarily because the analytical tools that have been applied were originally designed for autosomal markers. Our objective here is to employ a new robust and flexible tool for chromosome-wide analysis of X-linked markers in complex traits. Orofacial clefts are good candidates for such analysis because of the consistently observed excess of females with cleft palate only (CPO) and excess of males with cleft lip with or without cleft palate (CL/P).

Methods

Genotypes for 14,486 X-chromosome SNPs in 1,291 Asian and 1,118 European isolated cleft triads were available from a previously published GWAS. The R-package HAPLIN enables genome-wide–level analyses as well as statistical power simulations for a range of biologic scenarios. We analyzed isolated CL/P and isolated CPO for each ethnicity in HAPLIN, using a sliding-window approach to haplotype analysis and two different statistical models, with and without X-inactivation in females.

Results

There was a larger number of associations in the Asian versus the European sample, and similar to previous reports that have analyzed the same GWAS dataset using different methods, we identified associations with EFNB1/PJA1 and DMD. In addition, new associations were detected with several other genes, among which KLHL4, TBX22, CPXCR1 and BCOR were noteworthy because of their roles in clefting syndromes. A few of the associations were only detected by one particular X-inactivation model, whereas a few others were only detected in one sex.

Discussion/Conclusion

We found new support for the involvement of X-linked variants in isolated clefts. The associations were specific for ethnicity, sex and model parameterization, highlighting the need for flexible tools that are capable of detecting and estimating such effects. Further efforts are needed to verify and elucidate the potential roles of EFNB1/PJA1, KLHL4, TBX22, CPXCR1 and BCOR in isolated clefts.

The gene encoding the inwardly rectifying potassium channel Kir4.1 may be involved in sudden infant death syndrome

AIM

Disturbances in brain function and development may play a role in sudden infant death syndrome (SIDS). This Norwegian study aimed to test the hypothesis that specific variants of genes involved in water transport and potassium homeostasis would be predisposing factors for SIDS.

METHODS

Genetic variation in the genes encoding aquaporin-4 (AQP4), Kir4.1 (KCNJ10) and α-syntrophin was analysed in 171 SIDS cases (62.6% male) with a median age of 15.5 (2-52) weeks and 398 adult controls (70.6% male) with a median age of 44 (11-91) years. All the subjects were Caucasians who were autopsied from 1988 to 2013.

RESULTS

The CC genotype of rs72878794 in the AQP4 gene and a combination of the CC genotype in rs17375748, rs1130183, rs12133079 and rs1186688 in KCNJ10 (4xCC) were found to be associated with SIDS. The SIDS cases with the 4xCC SNP combination were younger than the SIDS cases with other genotype combinations (p = 0.006).

CONCLUSION

This study indicates that genetic variations in KCNJ10 and AQP4 may be predisposing factors for SIDS. Alterations in the expression of the AQP4/Kir4.1 complex can disrupt water and ion homeostasis, which may influence brain development and facilitate brain oedema formation This may be especially unfavourable during the first weeks of life.

A Population-Based Study of Effects of Genetic Loci on Orofacial Clefts

Prior genome-wide association studies for oral clefts have focused on clinic-based samples with unclear generalizability. Prior samples were also small for investigating effects by cleft type and exclusively studied isolated clefts (those occurring without other birth defects). We estimated the effects of 17 top loci on cleft types in both isolated and nonisolated cases in the largest consortium to date of European-descent population-based studies. Our analytic approach focused on a mother-child dyad case-control design, but it also allowed analyzing mother-only or child-only genotypes to maximize power. Our total sample included 1,875 cases with isolated clefts, 459 cases with nonisolated clefts, and 3,749 controls. After correcting for multiple testing, we observed significant associations between fetal single-nucleotide polymorphisms (SNPs) at IRF6, PAX7, 8q21.3, 8q24, KIAA1598-VAX1, and MAFB and isolated cleft lip only (CLO) and cleft lip and palate (CLP). Significant associations were observed between isolated CLO and fetal SNPs near TPM1 and NOG1 and between CLP and fetal SNPs at ABCA4-ARHGAP29, THADA, FOXE1, and SPRY2. Overall, effects were similar for isolated CLO and CLP, except for ABCA4-ARHGAP29. A protective effect was observed for the fetal NOG1 SNP on cleft palate only, opposite in direction to the effect on CLO. For most fetal SNPs, a dose-response allelic effect was observed. No evidence of parent-of-origin or maternal genome effects was observed. Overall, effect direction and magnitude were similar between isolated and nonisolated clefts, suggesting that several loci are modifiers of cleft risk in both isolated and nonisolated forms. Our results provide reliable estimates of the effects of top loci on risks of oral clefts in a population of European descent.

MSX1 gene polymorphisms in Mexican patients with non-syndromic cleft lip/palate

OBJECTIVE

Non-syndromic cleft lip/palate malformation (CL/P) is one of the most common birth defects in humans and has a complex etiology involving genetic and environmental factors. Mutations in the MSX1 gene are critical during craniofacial development. The purpose of this study was to investigate the contribution of MSX1 gene polymorphisms to the risk of developing CL/P in a sample of Mexican patients.

METHODS

The sample consisted of 282 subjects (69 cases and 213 relatives). Four single-nucleotide polymorphisms (SNP1, P147Q, SNP5 and P278S) were tested for association with CL/P in triad and case-pseudo-control analyses. Polymorphism typing was performed by restriction fragment length polymorphism and dot-blot techniques. Allele and genotype frequencies were calculated between patients and pseudo-controls and compared using the Chi square test with Yates correction. Odds ratios and 95% confidence intervals were obtained using SPSS software (v19). Triad analysis was also performed using the program HAPLIN (v5.3).

RESULTS

In the cases and pseudo-controls, an association was found between CL/P and the SNP1-G allele (P = 0.031) and the SNP1-G/G genotype (P = 0.032), a polymorphism located near MSX1. Triad analysis showed a tendency toward CL/P susceptibility for the genotype SNP1-G/G (P = 0.075) and an association between CL/P and the haplotype GCTC (P = 0.037). No associated haplotype was found in the cases and pseudo-controls. Two partial haplotypes, GT (SNP1-SNP5) (P = 0.032) and GC (SNP1-P278S) (P = 0.033), were associated with susceptibility in the heterozygous and homozygous types, respectively. In contrast, haplotype AT (SNP1-SNP5) was associated with protection (P = 0.012) in the homozygous type.

CONCLUSIONS

The results of this study suggest an association between CL/P susceptibility and SNP1, located near the MSX1 gene, in the Mexican population.

Analysis of Case-Parent Trios Using a Loglinear Model with Adjustment for Transmission Ratio Distortion

Transmission of the two parental alleles to offspring deviating from the Mendelian ratio is termed Transmission Ratio Distortion (TRD), occurs throughout gametic and embryonic development. TRD has been well-studied in animals, but remains largely unknown in humans. The Transmission Disequilibrium Test (TDT) was first proposed to test for association and linkage in case-trios (affected offspring and parents); adjusting for TRD using control-trios was recommended. However, the TDT does not provide risk parameter estimates for different genetic models. A loglinear model was later proposed to provide child and maternal relative risk (RR) estimates of disease, assuming Mendelian transmission. Results from our simulation study showed that case-trios RR estimates using this model are biased in the presence of TRD; power and Type 1 error are compromised. We propose an extended loglinear model adjusting for TRD. Under this extended model, RR estimates, power and Type 1 error are correctly restored. We applied this model to an intrauterine growth restriction dataset, and showed consistent results with a previous approach that adjusted for TRD using control-trios. Our findings suggested the need to adjust for TRD in avoiding spurious results. Documenting TRD in the population is therefore essential for the correct interpretation of genetic association studies.

Family-Based Multi-SNP X Chromosome Analysis Using Parent Information

We propose a method for association analysis of haplotypes on the X chromosome that offers both improved power and robustness to population stratification in studies of affected offspring and their parents if all three have been genotyped. The method makes use of assumed parental haplotype exchangeability (PHE), a weaker assumption than Hardy-Weinberg equilibrium (HWE). PHE requires that in the source population, of the three X chromosome haplotypes carried by the two parents, each is equally likely to be carried by the father. We propose a pseudo-sibling approach that exploits that exchangeability assumption. Our method extends the single-SNP PIX-LRT method to multiple SNPs in a high linkage block. We describe methods for testing the PHE assumption and also for determining how apparent violations can be distinguished from true fetal effects or maternally-mediated effects. We show results of simulations that demonstrate nominal type I error rate and good power. The methods are then applied to dbGaP data on the birth defect oral cleft, using both Asian and Caucasian families with cleft.

Expanding the toolbox of ADHD genetics. How can we make sense of parent of origin effects in ADHD and related behavioral phenotypes?

Genome-wide association (GWA) studies have shown that many different genetic variants cumulatively contribute to the risk of psychiatric disorders. It has also been demonstrated that various parent-of-origin effects (POE) may differentially influence the risk of these disorders. Together, these observations have provided important new possibilities to uncover the genetic underpinnings of such complex phenotypes. As POE so far have received little attention in neuropsychiatric disorders, there is still much progress to be made. Here, we mainly focus on the new and emerging role of POE in attention-deficit hyperactivity disorder (ADHD). We review the current evidence that POE play an imperative role in vulnerability to ADHD and related disorders. We also discuss how POE can be assessed using statistical genetics tools, expanding the resources of modern psychiatric genetics. We propose that better comprehension and inspection of POE may offer new insight into the molecular basis of ADHD and related phenotypes, as well as the potential for preventive and therapeutic interventions.

Increased Power for Detection of Parent-of-Origin Effects via the Use of Haplotype Estimation

Parent-of-origin (or imprinting) effects relate to the situation in which traits are influenced by the allele inherited from only one parent and the allele from the other parent has little or no effect. Given SNP genotype data from case-parent trios, the parent of origin of each allele in the offspring can often be deduced unambiguously; however, this is not true when all three individuals are heterozygous. Most existing methods for investigating parent-of-origin effects operate on a SNP-by-SNP basis and either perform some sort of averaging over the possible parental transmissions or else discard ambiguous trios. If the correct parent of origin at a SNP could be determined, this would provide extra information and increase the power for detecting the effects of imprinting. We propose making use of the surrounding SNP information, via haplotype estimation, to improve estimation of parent of origin at a test SNP for case-parent trios, case-mother duos, and case-father duos. This extra information is then used in a multinomial modeling approach for estimating parent-of-origin effects at the test SNP. We show through computer simulations that our approach has increased power over previous approaches, particularly when the data consist only of duos. We apply our method to two real datasets and find a decrease in significance of p values in genomic regions previously thought to possibly harbor imprinting effects, thus weakening the evidence that such effects actually exist in these regions, although some regions retain evidence of significant effects.

IRF6 polymorphisms in Mexican patients with non-syndromic cleft lip

Cleft lip with or without cleft palate (CL/P) is one of the most common birth defects; it is a multifactorial disease affecting > 1/1,000 live births in Europe, and its etiology is largely unknown, although it is very likely genetic and environmental factors contribute to this malformation. Orofacial development is a complex process involving many genes and signaling pathways. Mutations in the gene for the interferon regulatory factor 6 (IRF6) cause a hereditary dominant malformation syndrome including CL/P, and polymorphisms are associated with non-syndromic CL/P (MIM 119530). Five SNPs at the locus with high heterozygosity in Caucasian populations were chosen for the present research due to their very strong association with CL/P. A case–parent trio study was performed using 292 samples from Mexico. Association with the rs1319435-C/C genotype (P = 0.02) was found in patients (73) as compared to pseudocontrols (219), while the genotype rs1319435-T/C was related with protection (P = 0.041) in the triad design. Significant over-transmission of the G allele for marker rs2235375 (P = 0.049) was found. Only the TACGT haplotype was diminished in the affected child, either in single (P = 0.0208) or double (P = 0.0208) dose. The pairwise analysis showed rs2235543 and rs2235371 were in strong linkage disequilibrium. These results point to a substantial contribution of IRF6 in the etiology of non-syndromic CL/P in a sample of the Mexican population.

•

Cleft lip with or without cleft palate (CL/P) is one of the most common birth defects.

•

Association with the rs1319435-C/C genotype was found in patients as compared to pseudocontrols.

•

In the triad design, the genotype rs1319435-T/C was related with protection.

•

The IRF6 contributes to the etiology of non-syndromic CL/P in a sample of the Mexican population.

Learning about the X from our parents

The X chromosome is generally understudied in association studies, in part because the analyst has had limited methodological options. For nuclear-family-based association studies, most current methods extend the transmission disequilibrium test (TDT) to the X chromosome. We present a new method to study association in case-parent triads: the parent-informed likelihood ratio test for the X chromosome (PIX-LRT). Our method enables estimation of relative risks and takes advantage of parental genotype information and the sex of the affected offspring to increase statistical power to detect an effect. Under a parental exchangeability assumption for the X, if case-parent triads are complete, the parents of affected offspring provide an independent replication sample for estimates based on transmission distortion to their affected offspring. For each offspring sex we combine the parent-level and the offspring-level information to form a likelihood ratio test statistic; we then combine the two to form a combined test statistic. Our method can estimate relative risks under different modes of inheritance or a more general co-dominant model. In triads with missing parental genotypes, the method accounts for missingness with the Expectation-Maximization algorithm. We calculate non-centrality parameters to assess the power gain and robustness of our method compared to alternative methods. We apply PIX-LRT to publically available data from an international consortium of genotyped families affected by the birth defect oral cleft and find a strong, internally-replicated signal for a SNP marker related to cleft lip with or without cleft palate.

Review of statistical methodologies for the detection of parent-of-origin effects in family trio genome-wide association data with binary disease traits

The detection of parent-of-origin effects aims to identify whether the functionality of alleles, and in turn associated phenotypic traits, depends on the parental origin of the alleles. Different parent-of-origin effects have been identified through a variety of mechanisms and a number of statistical methodologies for their detection have been proposed, in particular for genome-wide association studies (GWAS). GWAS have had limited success in explaining the heritability of many complex disorders and traits, but successful identification of parent-of-origin effects using trio (mother, father and offspring) GWAS may help shed light on this missing heritability. However, it is important to choose the most appropriate parent-of-origin test or methodology, given knowledge of the phenotype, amount of available data and the type of parent-of-origin effect(s) being considered. This review brings together the parent-of-origin detection methodologies available, comparing them in terms of power and type I error for a number of different simulated data scenarios, and finally offering guidance as to the most appropriate choice for the different scenarios.

X-chromosomal maternal and fetal SNPs and the risk of spontaneous preterm delivery in a Danish/Norwegian genome-wide association study

BACKGROUND

Recent epidemiological studies suggest that the maternal genome is an important contributor to spontaneous preterm delivery (PTD). There is also a significant excess of males among preterm born infants, which may imply an X-linked mode of inheritance for a subset of cases. To explore this, we examined the effect of maternal and fetal X-chromosomal single nucleotide polymorphisms (SNPs) on the risk of PTD in two independent genome-wide association studies and one replication study.

METHODS

Participants were recruited from the Danish National Birth Cohort and the Norwegian Mother and Child cohort studies. Data from these two populations were first analyzed independently, and then combined in a meta-analysis. Overall, we evaluated 12,211 SNPs in 1,535 case-mother dyads and 1,487 control-mother dyads. Analyses were done using a hybrid design that combines case-mother dyads and control-mother dyads, as implemented in the Haplin statistical software package. A sex-stratified analysis was performed for the fetal SNPs. In the replication study, 10 maternal and 16 fetal SNPs were analyzed using case-parent triads from independent studies of PTD in the United States, Argentina and Denmark.

RESULTS

In the meta-analysis, the G allele at the maternal SNP rs2747022 in the FERM domain containing 7 gene (FRMD7) increased the risk of spontaneous PTD by 1.2 (95% confidence interval (CI): 1.1, 1.4). Although an association with this SNP was confirmed in the replication study, it was no longer statistically significant after a Bonferroni correction for multiple testing.

CONCLUSION

We did not find strong evidence in our data to implicate X-chromosomal SNPs in the etiology of spontaneous PTD. Although non-significant after correction for multiple testing, the mother's G allele at rs2747022 in FRMD7 increased the risk of spontaneous PTD across all populations in this study, thus warranting further investigation in other populations.

Maternal and offspring xenobiotic metabolism haplotypes and the risk of childhood acute lymphoblastic leukemia

Discovering genetic predictors of childhood acute lymphoblastic leukemia (ALL) necessitates the evaluation of novel factors including maternal genetic effects, which are a proxy for the intrauterine environment, and robust epidemiologic study designs. Therefore, we evaluated five maternal and offspring xenobiotic metabolism haplotypes and the risk of childhood ALL among 120 case-parent triads. Two of the five haplotypes were significantly associated with risk: GSTM3/GSTM4 (P=0.01) and GSTP1 (P=0.02). The EPHX1 haplotype was marginally associated with risk (P=0.05), whereas haplotypes in CYP1B1 and GSTA4 were not. Our results suggest genetic variation in xenobiotic metabolism is important in childhood ALL etiology.

A multi-locus likelihood method for assessing parent-of-origin effects using case-control mother-child pairs

Parent-of-origin effects have been pointed out to be one plausible source of the heritability that was unexplained by genome-wide association studies. Here, we consider a case-control mother-child pair design for studying parent-of-origin effects of offspring genes on neonatal/early-life disorders or pregnancy-related conditions. In contrast to the standard case-control design, the case-control mother-child pair design contains valuable parental information and therefore permits powerful assessment of parent-of-origin effects. Suppose the region under study is in Hardy-Weinberg equilibrium, inheritance is Mendelian at the diallelic locus under study, there is random mating in the source population, and the SNP under study is not related to risk for the phenotype under study because of linkage disequilibrium (LD) with other SNPs. Using a maximum likelihood method that simultaneously assesses likely parental sources and estimates effect sizes of the two offspring genotypes, we investigate the extent of power increase for testing parent-of-origin effects through the incorporation of genotype data for adjacent markers that are in LD with the test locus. Our method does not need to assume the outcome is rare because it exploits supplementary information on phenotype prevalence. Analysis with simulated SNP data indicates that incorporating genotype data for adjacent markers greatly help recover the parent-of-origin information. This recovery can sometimes substantially improve statistical power for detecting parent-of-origin effects. We demonstrate our method by examining parent-of-origin effects of the gene PPARGC1A on low birth weight using data from 636 mother-child pairs in the Jerusalem Perinatal Study.

A GCH1 haplotype and risk of neural tube defects in the National Birth Defects Prevention Study

Tetrahydrobiopterin (BH(4)) is an essential cofactor and an important cellular antioxidant. BH(4) deficiency has been associated with diseases whose etiologies stem from excessive oxidative stress. GTP cyclohydrolase I (GCH1) catalyzes the first and rate-limiting step of de novo BH(4) synthesis. A 3-SNP haplotype in GCH1 (rs8007267, rs3783641, and rs10483639) is known to modulate GCH1 gene expression levels and has been suggested as a major determinant of plasma BH(4) bioavailability. As plasma BH(4) bioavailability has been suggested as a mechanism of neural tube defect (NTD) teratogenesis, we evaluated the association between this GCH1 haplotype and the risk of NTDs. Samples were obtained from 760 NTD case-parent triads included in the National Birth Defects Prevention Study (NBDPS). The three SNPs were genotyped using TaqMan® SNP assays. An extension of the log-linear model was used to assess the association between NTDs and both offspring and maternal haplotypes. Offspring carrying two copies of haplotype C-T-C had a significantly increased NTD risk (risk ratio [RR]=3.40, 95% confidence interval [CI]: 1.02-11.50), after adjusting for the effect of the maternal haplotype. Additionally, mothers carrying two copies of haplotype C-T-C had a significantly increased risk of having an NTD-affected offspring (RR=3.46, 95% CI: 1.05-11.00), after adjusting for the effect of the offspring haplotype. These results suggest offspring and maternal variation in the GCH1 gene and altered BH(4) biosynthesis may contribute to NTD risk.

Association of TGFβ1 and SMAD4 variants in the etiology of keloid scar in the Malay population

Keloid is a complex condition with environmental and genetic risk-contributing factors. Two candidate genes, TGFβ1 and SMAD4, located in the same signaling pathway are highly expressed in the keloid fibroblast cells. In a case-control design, TGFβ1 haplotypes showed association with the risk of keloid in the present study. The CC haplotype, composed of both c.29C>T and -509T>C variants, was observed more frequently among cases (Corrected p = 0.037, OR = 2.07, 95 % CI = 0.87-4.93), showing a 4.5-fold increased risk for keloid. The AG genotype of the SMAD4 c.5131A>G variant showed a trend of significance (p = 0.0573, OR = 1.75, 95 % CI = 0.99-3.13). Taken together, either of these variants is most probably causative at the expression level or is in linkage disequilibrium with other causative variants in a complex pattern together with the environmental factors that contribute to the condition. To the best of our knowledge, there is only one documented report on a relationship between TGFβ1 and keloid with no association within the Caucasian population, while there have not been any reports for SMAD4. Therefore, the present study is likely the first research showing a significant association between TGFβ1 variants and keloids in the Malay population.

Application of a novel hybrid study design to explore gene-environment interactions in orofacial clefts

Orofacial clefts are common birth defects with strong evidence for both genetic and environmental causal factors. Candidate gene studies combined with exposures known to influence the outcome provide a highly targeted approach to detecting GxE interactions. We developed a new statistical approach that combines the case-control and offspring-parent triad designs into a "hybrid design" to search for GxE interactions among 334 autosomal cleft candidate genes and maternal first-trimester exposure to smoking, alcohol, coffee, folic acid supplements, dietary folate and vitamin A. The study population comprised 425 case-parent triads of isolated clefts and 562 control-parent triads derived from a nationwide study of orofacial clefts in Norway (1996-2001). A full maximum-likelihood model was used in combination with a Wald test statistic to screen for statistically significant GxE interaction between strata of exposed and unexposed mothers. In addition, we performed pathway-based analyses on 28 detoxification genes and 21 genes involved in folic acid metabolism. With the possible exception of the T-box 4 gene (TBX4) and dietary folate interaction in isolated CPO, there was little evidence overall of GxE interaction in our data. This study is the largest to date aimed at detecting interactions between orofacial clefts candidate genes and well-established risk exposures.

PREMIM and EMIM: tools for estimation of maternal, imprinting and interaction effects using multinomial modelling

Background

Here we present two new computer tools, PREMIM and EMIM, for the estimation of parental and child genetic effects, based on genotype data from a variety of different child-parent configurations. PREMIM allows the extraction of child-parent genotype data from standard-format pedigree data files, while EMIM uses the extracted genotype data to perform subsequent statistical analysis. The use of genotype data from the parents as well as from the child in question allows the estimation of complex genetic effects such as maternal genotype effects, maternal-foetal interactions and parent-of-origin (imprinting) effects. These effects are estimated by EMIM, incorporating chosen assumptions such as Hardy-Weinberg equilibrium or exchangeability of parental matings as required.

Results

In application to simulated data, we show that the inference provided by EMIM is essentially equivalent to that provided by alternative (competing) software packages such as MENDEL and LEM. However, PREMIM and EMIM (used in combination) considerably outperform MENDEL and LEM in terms of speed and ease of execution.

Conclusions

Together, EMIM and PREMIM provide easy-to-use command-line tools for the analysis of pedigree data, giving unbiased estimates of parental and child genotype relative risks.

X-linked genes and risk of orofacial clefts: evidence from two population-based studies in Scandinavia

BACKGROUND

Orofacial clefts are common birth defects of complex etiology, with an excess of males among babies with cleft lip and palate, and an excess of females among those with cleft palate only. Although genes on the X chromosome have been implicated in clefting, there has been no association analysis of X-linked markers.

METHODOLOGY/PRINCIPAL FINDINGS

We added new functionalities in the HAPLIN statistical software to enable association analysis of X-linked markers and an exploration of various causal scenarios relevant to orofacial clefts. Genotypes for 48 SNPs in 18 candidate genes on the X chromosome were analyzed in two population-based samples from Scandinavia (562 Norwegian and 235 Danish case-parent triads). For haplotype analysis, we used a sliding-window approach and assessed isolated cleft lip with or without cleft palate (iCL/P) separately from isolated cleft palate only (iCPO). We tested three statistical models in HAPLIN, allowing for: i) the same relative risk in males and females, ii) sex-specific relative risks, and iii) X-inactivation in females. We found weak but consistent associations with the oral-facial-digital syndrome 1 (OFD1) gene (formerly known as CXORF5) in the Danish iCL/P samples across all models, but not in the Norwegian iCL/P samples. In sex-specific analyses, the association with OFD1 was in male cases only. No analyses showed associations with iCPO in either the Norwegian or the Danish sample.

CONCLUSIONS

The association of OFD1 with iCL/P is plausible given the biological relevance of this gene. However, the lack of replication in the Norwegian samples highlights the need to verify these preliminary findings in other large datasets. More generally, the novel analytic methods presented here are widely applicable to investigations of the role of X-linked genes in complex traits.

Assessing the impact of nicotine dependence genes on the risk of facial clefts: An example of the use of national registry and biobank data

BACKGROUND

Maternal smoking during pregnancy has been associated with risk of facial clefts in offspring, but causation has not yet been established. It is possible that the effect of maternal smoking on facial clefts is mediated through genes that are involved in nicotine dependence. Gamma-aminobutyric acid B receptor 2 (GABBR2), dopa decarboxylase (DDC), and cholinergic receptor nicotinic alpha 4 (CHRNA4) are three examples of genes that have previously shown strong associations with nicotine dependence.

METHODS

We used a population-based sample of 377 case-parent trios of cleft lip with or without cleft palate (CL/P) and 762 control-parent trios from Norway (1996-2001) to investigate whether variants in GABBR2, DDC and CHRNA4 are associated with maternal first-trimester smoking and with clefting risk. We used HAPLIN (Gjessing et al. 2006), a statistical software tailored for family-based association tests, to perform haplotype-based analyses on 12 SNPs in these genes (rs10985765, rs1435252, rs3780422, rs2779562, and rs3750344 in GABBR2; rs2060762, rs3757472, rs1451371, rs3735273, and rs921451 in DDC; rs4522666 and rs1044393 in CHRNA4).

RESULTS

When analyzed one at a time, there was little evidence of association between any of the 12 SNPs and maternal first-trimester smoking. In haplotype analyses, however, one copy of the maternal G-G-c-G-c haplotype in DDC was linked with smoking prevalence (odds ratio: 1.5; 95% confidence interval: 1.0-2.1). This same haplotype also increased the risk of isolated CL/P in offspring by 1.5-fold with one copy and 2.4-fold with two copies (P_trend = 0.06). No statistically significant associations were detected with GABBR2 and CHRNA4.

CONCLUSIONS

Despite strong associations previously reported between nicotine dependence and variants in GABBR2, DDC and CHRNA4, these genes were poor predictors of maternal first-trimester smoking in our data. The direct association of the DDC haplotype with CL/P suggests that this haplotype may either have direct effects on clefts or it may influence clefting risks through other yet unexplored risk behavior(s).

Candidate gene analysis of spontaneous preterm delivery: new insights from re-analysis of a case-control study using case-parent triads and control-mother dyads

Background

Spontaneous preterm delivery (PTD) has a multifactorial etiology with evidence of a genetic contribution to its pathogenesis. A number of candidate gene case-control studies have been performed on spontaneous PTD, but the results have been inconsistent, and do not fully assess the role of how two genotypes can impact outcome. To elucidate this latter point we re-analyzed data from a previously published case-control candidate gene study, using a case-parent triad design and a hybrid design combining case-parent triads and control-mother dyads. These methods offer a robust approach to genetic association studies for PTD compared to traditional case-control designs.

Methods

The study participants were obtained from the Norwegian Mother and Child Cohort Study (MoBa). A total of 196 case triads and 211 control dyads were selected for the analysis. A case-parent triad design as well as a hybrid design was used to analyze 1,326 SNPs from 159 candidate genes. We compared our results to those from a previous case-control study on the same samples. Haplotypes were analyzed using a sliding window of three SNPs and a pathway analysis was performed to gain biological insight into the pathophysiology of preterm delivery.

Results

The most consistent significant fetal gene across all analyses was COL5A2. The functionally similar COL5A1 was significant when combining fetal and maternal genotypes. PON1 was significant with analytical approaches for single locus association of fetal genes alone, but was possibly confounded by maternal effects. Focal adhesion (hsa04510), Cell Communication (hsa01430) and ECM receptor interaction (hsa04512) were the most constant significant pathways.

Conclusion

This study suggests a fetal association of COL5A2 and a combined fetal-maternal association of COL5A1 with spontaneous PTD. In addition, the pathway analysis implied interactions of genes affecting cell communication and extracellular matrix.

Detection of intergenerational genetic effects with application to HLA-B matching as a risk factor for schizophrenia

BACKGROUND AND METHODS

Association studies using unrelated individuals cannot detect intergenerational genetic effects contributing to disease. To detect these effects, we improve the extended maternal-fetal genotype (EMFG) incompatibility test to estimate any combination of maternal effects, offspring effects, and their interactions at polymorphic loci or multiple SNPs, using any size pedigrees. We explore the advantages of using extended pedigrees rather than nuclear families. We apply our methods to schizophrenia pedigrees to investigate whether the previously associated mother-daughter HLA-B matching is a genuine risk or the result of bias.

RESULTS

Simulations demonstrate that using the EMFG test with extended pedigrees increases power and precision, while partitioning extended pedigrees into nuclear families can underestimate intergenerational effects. Application to actual data demonstrates that mother-daughter HLA-B matching remains a schizophrenia risk factor. Furthermore, ascertainment and mate selection biases cannot by themselves explain the observed HLA-B matching and schizophrenia association.

CONCLUSIONS

Our results demonstrate the power of the EMFG test to examine intergenerational genetic effects, highlight the importance of pedigree rather than case/control or case-mother/control-mother designs, illustrate that pedigrees provide a means to examine alternative, non-causal mechanisms, and they strongly support the hypothesis that HLA-B matching is causally involved in the etiology of schizophrenia in females.

Contribution of MSX1 variants to the risk of non-syndromic cleft lip and palate in a Malay population

Oral clefts are clinically and genetically heterogeneous disorders that are influenced by both genetic and environmental factors. The present family-based association study investigated the role of the MSX1 and TGFB3 genes in the etiology of non-syndromic oral cleft in a Malay population. No transmission distortion was found in the transmission disequilibrium analysis for either MSX1-CA or TGFB3-CA intragenic markers, whereas TGFB3-CA exhibited a trend to excess maternal transmission. In sequencing the MSX1 coding regions in 124 patients with oral cleft, five variants were found, including three known variants (A34G, G110G and P147Q) and two novel variants (M37L and G267A). The P147Q and M37L variants were not observed in 200 control chromosomes, whereas G267A was found in one control sample, indicating a very rare polymorphic variant. Furthermore, the G110G variant displayed a significant association between patients with non-syndromic cleft lip, with or without cleft palate, and normal controls (P=0.001, odds ratio=2.241, 95% confidence interval, 1.357-3.700). Therefore, these genetic variants may contribute, along with other genetic and environmental factors, to this condition.

Contribution of 6p24 to non-syndromic cleft lip and palate in a Malay population: association of variants in OFC1

Non-syndromic cleft lip, with or without cleft palate, is a heterogeneous, complex disease with a high incidence in the Asian population. Several association studies have been done on cleft candidate genes, but no reports have been published thus far on the Orofacial Cleft 1 (OFC1) genomic region in an Asian population. This study investigated the association between the OFC1 genomic region and non-syndromic cleft lip with or without cleft palate in 90 Malay father-mother-offspring trios. Results showed a preferential over-transmission of a 101-bp allele of marker D6S470 in the allele- and haplotype-based transmission disequilibrium test (TDT), as well as an excess of maternal transmission. However, no significant p-value was found for a maternal genotype effect in a log-linear model, although single and double doses of the 101-bp allele showed a slightly increased cleft risk (RR = 1.37, 95% CI, 0.527-3.4, p-value = 0.516). Carrying two copies of the 101-bp allele was significantly associated with an increased cleft risk (RR = 2.53, 95% CI, 1.06-6.12, p-value = 0.035). In conclusion, we report evidence of the contribution of the OFC1 genomic region to the etiology of clefts in a Malay population.

Fetal genetic risk of isolated cleft lip only versus isolated cleft lip and palate: a subphenotype analysis using two population-based studies of orofacial clefts in Scandinavia

Background

Cleft lip only (CLO) and cleft lip and palate (CLP) are commonly regarded as variants of the same defect and traditionally combined to form the single group of cleft lip with or without cleft palate (CL/P) prior to analysis. However, recent data have suggested that at least a subgroup of isolated CLO may be etiologically distinct from isolated CLP.

Methods

To explore fetal genetic risk of isolated CLO separately from isolated CLP, we performed a sub-phenotype analysis using two population-based studies of clefts in Scandinavia. One hundred and twenty-one isolated CLO, 190 isolated CLP, and 592 control triads were available from Norway(1996–2001), and a further 76 isolated CLO and 107 isolated CLP triads were available from Denmark (1991–2001). Genotypes for 1315 SNPs in 334 autosomal cleft candidate genes were analyzed using two complementary statistical methods, TRIMM and HAPLIN, to look for genetic associations across the two national samples.

Results

Both TRIMM and HAPLIN identified strong associations between FGF12 and isolated CLO in both populations. In addition, only TRIMM identified associations with IRF6 and VCL , and only HAPLIN found an association with CX43 . When analyses were repeated on the larger sample of isolated CLP, no significant associations were found with FGF12 , IRF6, VCL or CX43 .

Conclusions

Despite some inconsistency in the pattern of associations across the two populations, the associations themselves were phenotype-specific. While both IRF6 and FGF12 have previously shown strong associations with isolated CL/P, the associations with VCL and CX43 are novel and warrant further investigation in other isolated CLO samples.

Maternal genes and facial clefts in offspring: a comprehensive search for genetic associations in two population-based cleft studies from Scandinavia

BACKGROUND

Fetal conditions can in principle be affected by the mother's genotype working through the prenatal environment.

METHODOLOGY/PRINCIPAL FINDINGS

Genotypes for 1536 SNPs in 357 cleft candidate genes were available from a previous analysis in which we focused on fetal gene effects. After data-cleaning, genotypes for 1315 SNPs in 334 autosomal genes were available for the current analysis of maternal gene effects. Two complementary statistical methods, TRIMM and HAPLIN, were used to detect multi-marker effects in population-based samples from Norway (562 case-parent and 592 control-parent triads) and Denmark (235 case-parent triads). We analyzed isolated cleft lip with or without cleft palate (iCL/P) and isolated cleft palate only (iCP) separately and assessed replication by looking for genes detected in both populations by both methods. In iCL/P, neither TRIMM nor HAPLIN detected more genes than expected by chance alone; furthermore, the selected genes were not replicated across the two methods. In iCP, however, FLNB was identified by both methods in both populations. Although HIC1 and ZNF189 did not fully satisfy our stringency criterion for replication, they were strongly associated with iCP in TRIMM analyses of the Norwegian triads.

CONCLUSION/SIGNIFICANCE

Except for FLNB, HIC1 and ZNF189, maternal genes did not appear to influence the risk of clefting in our data. This is consistent with recent epidemiological findings showing no apparent difference between mother-to-offspring and father-to-offspring recurrence of clefts in these two populations. It is likely that fetal genes make the major genetic contribution to clefting risk in these populations, but we cannot rule out the possibility that maternal genes can affect risk through interactions with specific teratogens or fetal genes.

Maternal angiotensinogen (AGT) haplotypes, fetal renin (REN) haplotypes and risk of preeclampsia; estimation of gene-gene interaction from family-triad data

Background

Preeclampsia is a debilitating disorder affecting approximately 3% of pregnant women in the Western world. Although inconclusive, current evidence suggests that the renin-angiotensin system may be involved in hypertension. Therefore, our objective was to determine whether the genes for placental renin (REN) and maternal angiotensinogen (AGT) interact to influence the risk of preeclampsia.

Methods

Three haplotype-tagging SNPs (htSNPs) covering REN (rs5705, rs1464818, and rs3795575) and another three covering AGT (rs2148582, rs2478545 and rs943580) were genotyped in 99 mother-father-child triads of preeclampsia pregnancies. We estimated relative risks (RR) conferred by maternal AGT and fetal REN haplotypes using HAPLIN, a statistical software designed to detect multi-marker transmission distortion among triads. To assess a combined effect of maternal AGT and fetal REN haplotypes, the preeclamptic triads were first stratified by presence/absence of maternal AGT haplotype C-T-A and tested for an effect of fetal REN across these strata.

Results

We found evidence that mothers carrying the most frequent AGT haplotype, C-T-A, had a reduced risk of preeclampsia (RR of 0.4, 95% CI = 0.2-0.8 for heterozygotes and 0.6, 95% CI = 0.2-1.5 for homozygotes). Mothers homozygous for AGT haplotypes t-c-g and C-c-g appeared to have a higher risk, but only the former was statistically significant. We found only weak evidence of an overall effect of fetal REN haplotypes and no support for our hypothesis that an effect of REN depended on whether the mother carried the C-T-A haplotype of AGT (p = 0.33).

Conclusion

Our findings indicate that the mother's AGT haplotypes affect her risk for developing preeclampsia. However, this risk is not influenced by fetal REN haplotypes.