X-APL: an improved family-based test of association in the presence of linkage for the X chromosome

The X-factor in ART: does the use of assisted reproductive technologies influence DNA methylation on the X chromosome?

BACKGROUND

Assisted reproductive technologies (ART) may perturb DNA methylation (DNAm) in early embryonic development. Although a handful of epigenome-wide association studies of ART have been published, none have investigated CpGs on the X chromosome. To bridge this knowledge gap, we leveraged one of the largest collections of mother-father-newborn trios of ART and non-ART (natural) conceptions to date to investigate sex-specific DNAm differences on the X chromosome. The discovery cohort consisted of 982 ART and 963 non-ART trios from the Norwegian Mother, Father, and Child Cohort Study (MoBa). To verify our results from the MoBa cohort, we used an external cohort of 149 ART and 58 non-ART neonates from the Australian 'Clinical review of the Health of adults conceived following Assisted Reproductive Technologies' (CHART) study. The Illumina EPIC array was used to measure DNAm in both datasets. In the MoBa cohort, we performed a set of X-chromosome-wide association studies ('XWASs' hereafter) to search for sex-specific DNAm differences between ART and non-ART newborns. We tested several models to investigate the influence of various confounders, including parental DNAm. We also searched for differentially methylated regions (DMRs) and regions of co-methylation flanking the most significant CpGs. Additionally, we ran an analogous model to our main model on the external CHART dataset.

RESULTS

In the MoBa cohort, we found more differentially methylated CpGs and DMRs in girls than boys. Most of the associations persisted after controlling for parental DNAm and other confounders. Many of the significant CpGs and DMRs were in gene-promoter regions, and several of the genes linked to these CpGs are expressed in tissues relevant for both ART and sex (testis, placenta, and fallopian tube). We found no support for parental DNAm-dependent features as an explanation for the observed associations in the newborns. The most significant CpG in the boys-only analysis was in UBE2DNL, which is expressed in testes but with unknown function. The most significant CpGs in the girls-only analysis were in EIF2S3 and AMOT. These three loci also displayed differential DNAm in the CHART cohort.

CONCLUSIONS

Genes that co-localized with the significant CpGs and DMRs associated with ART are implicated in several key biological processes (e.g., neurodevelopment) and disorders (e.g., intellectual disability and autism). These connections are particularly compelling in light of previous findings indicating that neurodevelopmental outcomes differ in ART-conceived children compared to those naturally conceived.

Robust association tests for quantitative traits on the X chromosome

The genome-wide association study is an elementary tool to assess the genetic contribution to complex human traits. However, such association tests are mainly proposed for autosomes, and less attention has been given to methods for identifying loci on the X chromosome due to their distinct biological features. In addition, the existing association tests for quantitative traits on the X chromosome either fail to incorporate the information of males or only detect variance heterogeneity. Therefore, we propose four novel methods, which are denoted as QXcat, QZmax, QMVXcat and QMVZmax. When using these methods, it is assumed that the risk alleles for females and males are the same and that the locus being studied satisfies the generalized genetic model for females. The first two methods are based on comparing the means of the trait value across different genotypes, while the latter two methods test for the difference of both means and variances. All four methods effectively incorporate the information of X chromosome inactivation. Simulation studies demonstrate that the proposed methods control the type I error rates well. Under the simulated scenarios, the proposed methods are generally more powerful than the existing methods. We also apply our proposed methods to data from the Minnesota Center for Twin and Family Research and find 10 single nucleotide polymorphisms that are statistically significantly associated with at least two traits at the significance level of 1 × 10-3.

Detection of imprinting effects for qualitative traits on X chromosome based on nuclear families

Methods for detecting imprinting effects have been developed primarily for autosomal markers. However, no method is available in the literature to test for imprinting effects on X chromosome. Therefore, it is necessary to suggest methods for detecting such imprinting effects. In this article, the parental-asymmetry test on X chromosome (XPAT) is first developed to test for imprinting for qualitative traits in the presence of association, based on family trios each with both parents and their affected daughter. Then, we propose 1-XPAT to deal with parent-daughter pairs, each with one parent and his/her affected daughter. By simultaneously considering family trios and parent-daughter pairs, C-XPAT (the combined test statistic of XPAT and 1-XPAT) is constructed to test for imprinting. Further, we extend the proposed methods to accommodate complete (with both parents) and incomplete (with one parent) nuclear families having multiple daughters of which at least one is affected. Simulation results demonstrate that the proposed methods control the size well, irrespective of the inbreeding coefficient in females being zero or non-zero. By incorporating incomplete nuclear families, C-XPAT is more powerful than XPAT using only complete nuclear families. For practical use, these proposed methods are applied to analyse the rheumatoid arthritis data and Turner's syndrome data.

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.

Likelihood Ratio Test for Excess Homozygosity at Marker Loci on X Chromosome

The assumption of Hardy-Weinberg equilibrium (HWE) is generally required for association analysis using case-control design on autosomes; otherwise, the size may be inflated. There has been an increasing interest of exploring the association between diseases and markers on X chromosome and the effect of the departure from HWE on association analysis on X chromosome. Note that there are two hypotheses of interest regarding the X chromosome: (i) the frequencies of the same allele at a locus in males and females are equal and (ii) the inbreeding coefficient in females is zero (without excess homozygosity). Thus, excess homozygosity and significantly different minor allele frequencies between males and females are used to filter X-linked variants. There are two existing methods to test for (i) and (ii), respectively. However, their size and powers have not been studied yet. Further, there is no existing method to simultaneously detect both hypotheses till now. Therefore, in this article, we propose a novel likelihood ratio test for both (i) and (ii) on X chromosome. To further investigate the underlying reason why the null hypothesis is statistically rejected, we also develop two likelihood ratio tests for detecting (i) and (ii), respectively. Moreover, we explore the effect of population stratification on the proposed tests. From our simulation study, the size of the test for (i) is close to the nominal significance level. However, the size of the excess homozygosity test and the test for both (i) and (ii) is conservative. So, we propose parametric bootstrap techniques to evaluate their validity and performance. Simulation results show that the proposed methods with bootstrap techniques control the size well under the respective null hypothesis. Power comparison demonstrates that the methods with bootstrap techniques are more powerful than those without bootstrap procedure and the existing methods. The application of the proposed methods to a rheumatoid arthritis dataset indicates their utility.

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.

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.

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.

An X chromosome-wide association study in autism families identifies TBL1X as a novel autism spectrum disorder candidate gene in males

BACKGROUND

Autism spectrum disorder (ASD) is a complex neurodevelopmental disorder with a strong genetic component. The skewed prevalence toward males and evidence suggestive of linkage to the X chromosome in some studies suggest the presence of X-linked susceptibility genes in people with ASD.

METHODS

We analyzed genome-wide association study (GWAS) data on the X chromosome in three independent autism GWAS data sets: two family data sets and one case-control data set. We performed meta- and joint analyses on the combined family and case-control data sets. In addition to the meta- and joint analyses, we performed replication analysis by using the two family data sets as a discovery data set and the case-control data set as a validation data set.

RESULTS

One SNP, rs17321050, in the transducin β-like 1X-linked (TBL1X) gene [OMIM:300196] showed chromosome-wide significance in the meta-analysis (P value = 4.86 × 10-6) and joint analysis (P value = 4.53 × 10-6) in males. The SNP was also close to the replication threshold of 0.0025 in the discovery data set (P = 5.89 × 10-3) and passed the replication threshold in the validation data set (P = 2.56 × 10-4). Two other SNPs in the same gene in linkage disequilibrium with rs17321050 also showed significance close to the chromosome-wide threshold in the meta-analysis.

CONCLUSIONS

TBL1X is in the Wnt signaling pathway, which has previously been implicated as having a role in autism. Deletions in the Xp22.2 to Xp22.3 region containing TBL1X and surrounding genes are associated with several genetic syndromes that include intellectual disability and autistic features. Our results, based on meta-analysis, joint analysis and replication analysis, suggest that TBL1X may play a role in ASD risk.

Infantile hypertrophic pyloric stenosis: evaluation of three positional candidate genes, TRPC1, TRPC5 and TRPC6, by association analysis and re-sequencing

Infantile hypertrophic pyloric stenosis (IHPS) is the most common inherited form of gastrointestinal obstruction in infancy with a striking male preponderance. Infants present with vomiting due to gastric outlet obstruction caused by hypertrophy of the smooth muscle of the pylorus. Two loci specific to extended pedigrees displaying autosomal dominant inheritance have been identified. A genome scan identified loci on chromosomes 11q14-q22 and Xq23-q24 which are predicted to be responsible for a subset of smaller families with IHPS demonstrating non-Mendelian inheritance. The two linked chromosomal regions both harbour functional candidate genes which are members of the canonical transient receptor potential (TRPC) family of ion channels. Both TRPC5 (Xq23-q24) and TRPC6 (11q14-q22) have a potential role in smooth muscle control and hypertrophy. Here, we report suggestive evidence for a third locus on chromosome 3q12-q25 (Zmax = 2.7, p < 0.004), a region which harbours a third TRPC gene, TRPC1. Fine mapping of all three genes using a tagSNP approach and re-sequencing identified a SNP in the promoter region of TRPC6 and a missense variant in exon 4 of TRPC6 which may be putative causal variants.

Variants in several genomic regions associated with asperger disorder

Asperger disorder (ASP) is one of the autism spectrum disorders (ASD) and is differentiated from autism largely on the absence of clinically significant cognitive and language delays. Analysis of a homogenous subset of families with ASP may help to address the corresponding effect of genetic heterogeneity on identifying ASD genetic risk factors. To examine the hypothesis that common variation is important in ASD, we performed a genome-wide association study (GWAS) in 124 ASP families in a discovery data set and 110 ASP families in a validation data set. We prioritized the top 100 association results from both cohorts by employing a ranking strategy. Novel regions on 5q21.1 (P = 9.7 × 10(-7) ) and 15q22.1-q22.2 (P = 7.3 × 10(-6) ) were our most significant findings in the combined data set. Three chromosomal regions showing association, 3p14.2 (P = 3.6 × 10(-6) ), 3q25-26 (P = 6.0 × 10(-5) ) and 3p23 (P = 3.3 × 10(-4) ) overlapped linkage regions reported in Finnish ASP families, and eight association regions overlapped ASD linkage areas. Our findings suggest that ASP shares both ASD-related genetic risk factors, as well as has genetic risk factors unique to the ASP phenotype.

Sex chromosomes and genetic association studies

Although the literature concerning statistical testing for genotype-phenotype association in family-based and population-based studies is very extensive, until recently the sex chromosomes have received little attention. Here it is shown that the X chromosome in particular presents special problems with respect to efficient analysis of mixed-sex population studies, and as a result of X inactivation. This paper reviews recent developments in approaching these problems.

Common genetic variants on 5p14.1 associate with autism spectrum disorders

Autism spectrum disorders (ASDs) represent a group of childhood neurodevelopmental and neuropsychiatric disorders characterized by deficits in verbal communication, impairment of social interaction, and restricted and repetitive patterns of interests and behaviour. To identify common genetic risk factors underlying ASDs, here we present the results of genome-wide association studies on a cohort of 780 families (3,101 subjects) with affected children, and a second cohort of 1,204 affected subjects and 6,491 control subjects, all of whom were of European ancestry. Six single nucleotide polymorphisms between cadherin 10 (CDH10) and cadherin 9 (CDH9)-two genes encoding neuronal cell-adhesion molecules-revealed strong association signals, with the most significant SNP being rs4307059 (P = 3.4 x 10(-8), odds ratio = 1.19). These signals were replicated in two independent cohorts, with combined P values ranging from 7.4 x 10(-8) to 2.1 x 10(-10). Our results implicate neuronal cell-adhesion molecules in the pathogenesis of ASDs, and represent, to our knowledge, the first demonstration of genome-wide significant association of common variants with susceptibility to ASDs.

Association test for X-linked QTL in family-based designs

Family-based association methods for detecting quantitative trait loci (QTL) have been developed primarily for autosomes, and comparable methods for X-linked QTL have received less attention. We have developed a family-based association test for quantitative traits, named XQTL, which uses X-linked markers in a nuclear family design. XQTL adopts the framework of the orthogonal model implemented in the QTDT program, modifying the sex-specific score for X-linked genotypes. XQTL also takes into account the dosage effect due to female X chromosome inactivation. Restricted maximum likelihood (REML) and Fisher's scoring method are used to estimate variance components of random effects. Fixed effects, derived from the phenotypic differences among and within families, are estimated by the least-squares method. Our proposed XQTL can perform allelic and two-locus haplotypic association tests and can provide estimates of additive genetic effects and variance components. Simulation studies show correct type I error rates under the null hypothesis and robust statistical power under alternative scenarios. The loss of power observed when parental genotypes are missing can be compensated by an increase of offspring number. By treating age at onset of Parkinson disease as a quantitative trait, we illustrate our method, using MAO polymorphisms in 780 families.

X-LRT: a likelihood approach to estimate genetic risks and test association with X-linked markers using a case-parents design

Recently, there has been interest in family-based tests of association to identify X-chromosome genes. However, none of the approaches allow for estimation of genetic risks. We propose a likelihood approach to estimate disease-related marker relative risks and test genotype association using a case-parents design. The test uses nuclear families with a single affected proband and allows additional siblings and missing parental genotypes. Extension to a haplotype test is based on assumptions of random mating and multiplicative penetrance. We investigate power and type I error rate of the likelihood-based test, using simulated data and apply our method to marker data from the monoamine oxidase A&B genes in families with Parkinson disease. We show how efficiency with missing parental information can be improved with additional sibling genotype information. Our likelihood approach offers great flexibility for testing different penetrance relationships within and between sexes. In addition, estimation of disease-related marker relative risks provides a measure of the magnitude of X-linked genetic effects on complex disorders.

Family-based methods for linkage and association analysis

Traditional epidemiological study concepts such as case-control or cohort designs can be used in the design of genetic association studies, giving them a prominent role in genetic association analysis. A different class of designs based on related individuals, typically families, uses the concept of Mendelian transmission to achieve design-independent randomization, which permits the testing of linkage and association. Family-based designs require specialized analytic methods but they have distinct advantages: They are robust to confounding and variance inflation, which can arise in standard designs in the presence of population substructure; they test for both linkage and association; and they offer a natural solution to the multiple comparison problem. This chapter focuses on family-based designs. We describe some basic study designs as well as general approaches to analysis for qualitative, quantitative, and complex traits. Finally, we review available software.

Substance dependence low-density whole genome association study in two distinct American populations

Cocaine and opioid dependence are common, complex disorders with high heritability that commonly co-occur with other substance dependence disorders. Improved insight into the genetic basis of substance dependence would help elucidate its etiology and could inform its prevention and treatment. To generate new hypotheses about the genetics of substance dependence, we genotyped 5633 tagging single nucleotide polymorphism (SNP) markers in 1699 subjects from 339 African American (AA) families and 334 European American (EA) families ascertained through a sib pair meeting DSM-IV criteria for either cocaine or opioid dependence. The associations between genetic markers and five substance dependence traits (cocaine dependence, opioid dependence, cocaine-induced paranoia, alcohol dependence, and nicotine dependence) were assessed by family based association tests (FBAT). Results were ranked according to several criteria including statistical significance, concordance of results across population samples, and potential biological relevance of the implicated gene. The top-ranked result was an association of SNP rs1133503 in the MANEA gene with cocaine-induced paranoia (CIP). Our study provides an initial substance dependence trait-specific blueprint of associated regions for future candidate gene studies.