Transmission/disequilibrium tests for extended marker haplotypes

Estimation of Direct and Indirect Polygenic Effects and Gene-Environment Interactions using Polygenic Scores in Case-Parent Trio Studies

Family-based studies provide a unique opportunity to characterize genetic risks of diseases in the presence of population structure, assortative mating, and indirect genetic effects. We propose a novel framework, PGS-TRI, for the analysis of polygenic scores (PGS) in case-parent trio studies for estimation of the risk of an index condition associated with direct effects of inherited PGS, indirect effects of parental PGS, and gene-environment interactions. Extensive simulation studies demonstrate the robustness of PGS-TRI in the presence of complex population structure and assortative mating compared to alternative methods. We apply PGS-TRI to multi-ancestry trio studies of autism spectrum disorders (Ntrio = 1,517) and orofacial clefts (Ntrio = 1,904) to establish the first transmission-based estimates of risk associated with pre-defined PGS for these conditions and other related traits. For both conditions, we further explored offspring risk associated with polygenic gene-environment interactions, and direct and indirect effects of genetically predicted levels of gene expression and metabolite traits.

Powerful testing via hierarchical linkage disequilibrium in haplotype association studies

Marginal tests based on individual SNPs are routinely used in genetic association studies. Studies have shown that haplotype-based methods may provide more power in disease mapping than methods based on single markers when, for example, multiple disease-susceptibility variants occur within the same gene. A limitation of haplotype-based methods is that the number of parameters increases exponentially with the number of SNPs, inducing a commensurate increase in the degrees of freedom and weakening the power to detect associations. To address this limitation, we introduce a hierarchical linkage disequilibrium model for disease mapping, based on a reparametrization of the multinomial haplotype distribution, where every parameter corresponds to the cumulant of each possible subset of a set of loci. This hierarchy present in the parameters enables us to employ flexible testing strategies over a range of parameter sets: from standard single SNP analyses through the full haplotype distribution tests, reducing degrees of freedom and increasing the power to detect associations. We show via extensive simulations that our approach maintains the type I error at nominal level and has increased power under many realistic scenarios, as compared to single SNP and standard haplotype-based studies. To evaluate the performance of our proposed methodology in real data, we analyze genome-wide data from the Wellcome Trust Case-Control Consortium.

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.

Regionally Smoothed Meta-Analysis Methods for GWAS Datasets

Genome-wide association studies are proven tools for finding disease genes, but it is often necessary to combine many cohorts into a meta-analysis to detect statistically significant genetic effects. Often the component studies are performed by different investigators on different populations, using different chips with minimal SNPs overlap. In some cases, raw data are not available for imputation so that only the genotyped single nucleotide polymorphisms (SNPs) results can be used in meta-analysis. Even when SNP sets are comparable, different cohorts may have peak association signals at different SNPs within the same gene due to population differences in linkage disequilibrium or environmental interactions. We hypothesize that the power to detect statistical signals in these situations will improve by using a method that simultaneously meta-analyzes and smooths the signal over nearby markers. In this study, we propose regionally smoothed meta-analysis methods and compare their performance on real and simulated data.

Genes and environment in neonatal intraventricular hemorrhage

Emerging data suggest intraventricular hemorrhage (IVH) of the preterm neonate is a complex disorder with contributions from both the environment and the genome. Environmental analyses suggest factors mediating both cerebral blood flow and angiogenesis contribute to IVH, while candidate gene studies report variants in angiogenesis, inflammation, and vascular pathways. Gene-by-environment interactions demonstrate the interaction between the environment and the genome, and a non-replicated genome-wide association study suggests that both environmental and genetic factors contribute to the risk for severe IVH in very low-birth weight preterm neonates.

A multi-SNP association test for complex diseases incorporating an optimal P-value threshold algorithm in nuclear families

Background

Genome-wide association studies (GWAS) have become a common approach to identifying single nucleotide polymorphisms (SNPs) associated with complex diseases. As complex diseases are caused by the joint effects of multiple genes, while the effect of individual gene or SNP is modest, a method considering the joint effects of multiple SNPs can be more powerful than testing individual SNPs. The multi-SNP analysis aims to test association based on a SNP set, usually defined based on biological knowledge such as gene or pathway, which may contain only a portion of SNPs with effects on the disease. Therefore, a challenge for the multi-SNP analysis is how to effectively select a subset of SNPs with promising association signals from the SNP set.

Results

We developed the Optimal P-value Threshold Pedigree Disequilibrium Test (OPTPDT). The OPTPDT uses general nuclear families. A variable p-value threshold algorithm is used to determine an optimal p-value threshold for selecting a subset of SNPs. A permutation procedure is used to assess the significance of the test. We used simulations to verify that the OPTPDT has correct type I error rates. Our power studies showed that the OPTPDT can be more powerful than the set-based test in PLINK, the multi-SNP FBAT test, and the p-value based test GATES. We applied the OPTPDT to a family-based autism GWAS dataset for gene-based association analysis and identified MACROD2-AS1 with genome-wide significance (p-value= 2.5 × 10^− 6).

Conclusions

Our simulation results suggested that the OPTPDT is a valid and powerful test. The OPTPDT will be helpful for gene-based or pathway association analysis. The method is ideal for the secondary analysis of existing GWAS datasets, which may identify a set of SNPs with joint effects on the disease.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1620-3) contains supplementary material, which is available to authorized users.

Adaptive robust genetic association tests using case-parents triad families

Using the conditional likelihood, a model-specific score test can be derived under a given genetic model to test genetic association for case-parents triad family data. When the underlying genetic model is correctly specified, the score test is most powerful. However, it can lose substantial power when the model is misspecified. Several robust tests have been proposed to deal with the problem, such as the maximum test statistic, the maximin efficiency robust test, and the constrained likelihood ratio test. These tests have been shown to be robust against model misspecification compared with those model-specific score tests, but they are either time-consuming in computation or not sufficiently high in power robustness under some situations. In this study, a data-driven procedure is proposed to construct two adaptive robust genetic association tests WMERT and WMAX . The WMERT is simple in calculation and has fairly high power robustness. The empirical power of WMAX is quite stable and close to those of the model-specific score tests. The two proposed tests should be beneficial to practical genetic association studies. A real dataset consisting of neural tube defect triad families is used for illustration of the methods. R-scripts are also provided for numerical calculation of the proposed methods in practical studies.

Polymorphisms in immune mediators associate with risk of cervical cancer

OBJECTIVE

The immune system is critical for controlling the progression of HPV cervical disease and the development of cancer. This study aimed to identify cervical cancer susceptibility alleles in candidate immune-modulating genes.

METHODS

Our family-based study involved a cohort of 641 probands (women with ICC/CIN III) and their biologic parents or siblings (641 trios). In the discovery phase (stage 1), involving 288 of the trios, 80 tag single nucleotide polymorphisms (SNPs) in 11 immune-modulating genes (IFNG, IFNGR1, IFNGR2, JAK1, JAK2, STAT1, STAT6, IL12A, TNF, LTA and LTB) were evaluated on the GoldenGate platform. We used the combined dataset for a total of 641 trios (stage 2) and the Taqman platform to validate the SNPs that had proved significant in the discovery dataset. The transmission disequilibrium test was used to detect significant shifts in allelic transmissions in the datasets.

RESULTS

Two SNPs in JAK2 and one SNP in STAT6 showed significant allelic association with cervical cancer in the stage 1 discovery dataset and were replicated in the larger joint analysis stage 2 dataset (JAK2 rs10815144, P=0.0029 and rs12349785, P=0.0058; and STAT6 rs3024971, P=0.0127). An additional SNP in exon 19 of JAK2 (rs2230724) was also examined in the combined dataset due to its strong linkage disequilibrium (LD) with rs10815144. It was also significant (P=0.0335).

CONCLUSIONS

Our results suggest an association of SNPs in JAK2 and STAT6 with cervical cancer. This association should be investigated in additional cervical cancer populations.

Imputation without doing imputation: a new method for the detection of non-genotyped causal variants

Genome-wide association studies allow detection of non-genotyped disease-causing variants through testing of nearby genotyped SNPs. This approach may fail when there are no genotyped SNPs in strong LD with the causal variant. Several genotyped SNPs in weak LD with the causal variant may, however, considered together, provide equivalent information. This observation motivates popular but computationally intensive approaches based on imputation or haplotyping. Here we present a new method and accompanying software designed for this scenario. Our approach proceeds by selecting, for each genotyped "anchor" SNP, a nearby genotyped "partner" SNP, chosen via a specific algorithm we have developed. These two SNPs are used as predictors in linear or logistic regression analysis to generate a final significance test. In simulations, our method captures much of the signal captured by imputation, while taking a fraction of the time and disc space, and generating a smaller number of false-positives. We apply our method to a case/control study of severe malaria genotyped using the Affymetrix 500K array. Previous analysis showed that fine-scale sequencing of a Gambian reference panel in the region of the known causal locus, followed by imputation, increased the signal of association to genome-wide significance levels. Our method also increases the signal of association from P ≈ 2 × 10⁻⁶ to P ≈ 6 × 10⁻¹¹. Our method thus, in some cases, eliminates the need for more complex methods such as sequencing and imputation, and provides a useful additional test that may be used to identify genetic regions of interest.

Gene genealogies for genetic association mapping, with application to Crohn's disease

A gene genealogy describes relationships among haplotypes sampled from a population. Knowledge of the gene genealogy for a set of haplotypes is useful for estimation of population genetic parameters and it also has potential application in finding disease-predisposing genetic variants. As the true gene genealogy is unknown, Markov chain Monte Carlo (MCMC) approaches have been used to sample genealogies conditional on data at multiple genetic markers. We previously implemented an MCMC algorithm to sample from an approximation to the distribution of the gene genealogy conditional on haplotype data. Our approach samples ancestral trees, recombination and mutation rates at a genomic focal point. In this work, we describe how our sampler can be used to find disease-predisposing genetic variants in samples of cases and controls. We use a tree-based association statistic that quantifies the degree to which case haplotypes are more closely related to each other around the focal point than control haplotypes, without relying on a disease model. As the ancestral tree is a latent variable, so is the tree-based association statistic. We show how the sampler can be used to estimate the posterior distribution of the latent test statistic and corresponding latent p-values, which together comprise a fuzzy p-value. We illustrate the approach on a publicly-available dataset from a study of Crohn's disease that consists of genotypes at multiple SNP markers in a small genomic region. We estimate the posterior distribution of the tree-based association statistic and the recombination rate at multiple focal points in the region. Reassuringly, the posterior mean recombination rates estimated at the different focal points are consistent with previously published estimates. The tree-based association approach finds multiple sub-regions where the case haplotypes are more genetically related than the control haplotypes, and that there may be one or multiple disease-predisposing loci.

Association study in three different populations between the GPR88 gene and major psychoses

GPR88, coding for a G protein-coupled orphan receptor that is highly represented in the striatum, is a strong functional candidate gene for neuropsychiatric disorders and is located at 1p22-p21, a chromosomal region that we have previously linked to bipolar disorder (BD) in the Sardinian population. In order to ascertain the relevance of GPR88 as a risk factor for psychiatric diseases, we performed a genetic association analysis between GPR88 and BD in a sample of triads (patient and both parents) recruited in the Sardinian and the Palestinian population as well as between GPR88 and schizophrenia (SZ) in triads from the Xhosa population in South Africa. We found a positive association between GPR88 and BD in the Sardinian and Palestinian triads. Moreover, we found a positive association between GPR88 and SZ in triads from the Xhosa population in South Africa. When these results were corrected for multiple testing, the association between GPR88 and BD was maintained in the Palestinian population. Thus, these results suggest that GPR88 deserves consideration as a candidate gene for psychiatric diseases and requires to be further investigated in other populations.

Genetic analysis of common variants in the CMYA5 (cardiomyopathy-associated 5) gene with schizophrenia

Recently, CMYA5 was suggested as a susceptibility gene for schizophrenia based on two independent studies utilizing different ethnic samples. We designed a case-control study to examine whether 21 SNPs contained within CMYA5 were associated with the disorder in a western Han Chinese sample comprised of 488 schizophrenia patients and 516 healthy control subjects. The allele distribution of SNPs rs7714250, rs16877135 and rs13158477 showed significant association with schizophrenia (Puncorrected=0.008, Puncorrected=0.04, and Puncorrected=0.009, respectively) as well as the genotype distribution in the Cochran-Armitage trend test (Puncorrected=0.008, Puncorrected=0.037 and Puncorrected=0.011, respectively). After Bonferroni correction, rs7714250 showed a trend of association with schizophrenia both in allele distribution (Pcorrected=0.088) and genotype distribution (Pcorrected=0.088). Furthermore, significant associations were found in several two-, three-, four-, and five-SNP tests of haplotype analyses. Replications of the association of CMYA5 with schizophrenia across various studies suggest that it is very likely a potential common schizophrenia-related gene worldwide. Functional studies correlating CMYA5 with DTNBP1 and PKA warrant further investigation of the molecular basis of this gene in relationship to the signal transduction pathway(s) underlying the pathogenesis of schizophrenia.

Genetic variations in EGFR and ERBB4 increase susceptibility to cervical cancer

OBJECTIVES

Inherited genetic variability contributes to susceptibility to cervical cancer. We investigated the association of single nucleotide polymorphisms (SNPs) in the human epidermal growth factor receptor (ERBB) family with cervical cancer.

METHODS

We used the transmission disequilibrium test (TDT) to look for excessive transmission of tag single nucleotide polymorphisms (tSNPs) in ERBB family members EGFR, ERBB2, ERBB3, and ERBB4 in a large sample of women with invasive and in situ cervical cancer and their biological parents (628 trios). The study used a discovery set of trios (244) analyzed by Illumina GoldenGate in which SNPs reaching a P<.05 were re-tested by TaqMan in the combined set of 628. We also explored collaborative effects of different ERBB alleles.

RESULTS

Based on single SNP TDT tests we identified 16 significant SNPs in the discover stage and six of 14 SNPs that could be assayed by TaqMan were significantly overtransmitted in women with cervical cancer in the combined replication set. Four SNPs were located in intron 1 of EGFR and two SNPs in intron 24 of ERBB4. The EGFR variants are located near multiple enhancers, silencers, and the previously identified functional common polymorphisms in intron 1.

CONCLUSIONS

Our data provide evidence for the involvement of intron 1 EGFR variants and intron 24 ERBB4 variants in modulating risk for the development of in situ and invasive cervical cancer. These variants should be examined in additional populations and functional studies would be needed to confirm this hypothesis.

A family-based association analysis and meta-analysis of the reading disabilities candidate gene DYX1C1

Reading disabilities (RD) have a significant genetic basis and have shown linkage to multiple regions including chromosome 15q. Dyslexia susceptibility 1 candidate gene 1 (DYX1C1) on chromosome 15q21 was originally proposed as a candidate gene with two potentially functional polymorphisms at the -3G/A and 1249G/T positions showing association with RD. However, subsequent studies have yielded mixed results. We performed a literature review and meta-analysis of the -3G/A and 1249G/T polymorphisms, including new unpublished data from two family-based samples. Ten markers in DYX1C1 were genotyped in the two independently ascertained samples. Single marker and -3G/A:1249G/T haplotype analyses were performed for RD in both samples, and quantitative trait analyses using standardized reading-related measures was performed in one of the samples. For the meta-analysis, we used a random-effects model to summarize studies that tested for association between -3G/A or 1249G/T and RD. No significant association was found between the DYX1C1 SNPs and RD or any of the reading-related measures tested after correction for the number of tests performed. The previously reported risk haplotype (-3A:1249T) was not biased in transmission. A total of 9 and 10 study samples were included in the meta-analysis of the -3G/A and 1249G/T polymorphisms, respectively. Neither polymorphism reached statistical significance, but the heterogeneity for the 1249G/T polymorphism was high. The results of this study do not provide evidence for association between the putatively functional SNPs -3G/A and 1249G/T and RD.

Further evidence for the association of genetic variants of ZNF804A with schizophrenia and a meta-analysis for genome-wide significance variant rs1344706

Recent accumulating evidence has indicated that ZNF804A (zinc finger protein 804A) may be one of the most robustly implicated genes in schizophrenia. In this report, we examined ZNF804A single nucleotide polymorphisms (SNPs) encompassing exon 4 by performing an association study that used a Han Chinese sample comprised of 492 schizophrenia patients and 516 healthy control subjects. A meta-analysis based on previous studies was also performed. For markers rs4667000 and rs1366842, significant differences in allele frequencies were found between cases and controls (Mantel-Haenszel corrected P=0.014 and P=0.025, respectively). Analysis of haplotype rs61739290-rs1366842 showed significant association with schizophrenia (global P=0.0018). Moreover, several other two-, three-, and four-SNP tests of haplotype association were also significant. A meta-analysis comprised of studies that utilized sample sets of either European and/or Han Chinese origin revealed statistically significant associations for two SNPs (rs1366842, P=0.002; and rs3731834, P=0.03) and schizophrenia. In addition, we observed a significant association between marker rsl344706 and schizophrenia (P<1.0×10(-5)) in combined populations. When we separately analyzed the studies by population, consistent and significant differences were found between cases and controls both in the European samples (P<1.0×10(-4)) and in the Chinese samples (P=0.03). In summary, we have added new evidence supporting the association between ZNF804A and schizophrenia in our Han Chinese sample. Further functional exploration of ZNF804A will greatly help us to elucidate the pathogenesis of schizophrenia and find promising new approaches for the treatment of this disorder.

Wound healing genes and susceptibility to cutaneous leishmaniasis in Brazil

Leishmania braziliensis causes cutaneous (CL) and mucosal (ML) leishmaniasis. In the mouse, Fli1 was identified as a gene influencing enhanced wound healing and resistance to CL caused by Leishmania major. Polymorphism at FLI1 is associated with CL caused by L. braziliensis in humans, with an inverse association observed for ML disease. Here we extend the analysis to look at other wound healing genes, including CTGF, TGFB1, TGFBR1/2, SMADS 2/3/4/7 and FLII, all functionally linked along with FLI1 in the TGF beta pathway. Haplotype tagging single nucleotide polymorphisms (tag-SNPs) were genotyped using Taqman technology in 325 nuclear families (652 CL cases; 126 ML cases) from Brazil. Robust case-pseudocontrol (CPC) conditional logistic regression analysis showed associations between CL and SNPs at CTGF (SNP rs6918698; CC genotype; OR 1.67; 95%CI 1.10-2.54; P=0.016), TGFBR2 (rs1962859; OR 1.50; 95%CI 1.12-1.99; P=0.005), SMAD2 (rs1792658; OR 1.57; 95%CI 1.04-2.38; P=0.03), SMAD7 (rs4464148; AA genotype; OR 2.80; 95%CI 1.00-7.87; P=0.05) and FLII (rs2071242; OR 1.60; 95%CI 1.14-2.24; P=0.005), and between ML and SNPs at SMAD3 (rs1465841; OR 2.15; 95%CI 1.13-4.07; P=0.018) and SMAD7 (rs2337107; TT genotype; OR 3.70; 95%CI 1.27-10.7; P=0.016). Stepwise logistic regression analysis showed that all SNPs associated with CL at FLI1, CTGF, TGFBR2, and FLII showed independent effects from each other, but SNPs at SMAD2 and SMAD7 did not add independent effects to SNPs from other genes. These results suggest that TGFβ signalling via SMAD2 is important in directing events that contribute to CL, whereas signalling via SMAD3 is important in ML. Both are modulated by the inhibitory SMAD7 that acts upstream of SMAD2 and SMAD3 in this signalling pathway. Along with the published FLI1 association, these data further contribute to the hypothesis that wound healing processes are important determinants of pathology associated with cutaneous forms of leishmaniasis.

A sparse transmission disequilibrium test for haplotypes based on Bradley-Terry graphs

BACKGROUND

Linkage and association analysis based on haplotype transmission disequilibrium can be more informative than single marker analysis. Several works have been proposed in recent years to extend the transmission disequilibrium test (TDT) to haplotypes. Among them, a powerful approach called the evolutionary tree TDT (ET-TDT) incorporates information about the evolutionary relationship among haplotypes using the cladogram of the locus.

METHODS

In this work we extend this approach by taking into consideration the sparsity of causal mutations in the evolutionary history. We first introduce the notion of a Bradley-Terry (BT) graph representation of a haplotype locus. The most important property of the BT graph is that sparsity of the edge set of the graph corresponds to small number of causal mutations in the evolution of the haplotypes. We then propose a method to test the null hypothesis of no linkage and association against sparse alternatives under which a small number of edges on the BT graph have non-nil effects.

RESULTS AND CONCLUSION

We compare the performance of our approach to that of the ET-TDT through a power study, and show that incorporating sparsity of causal mutations can significantly improve the power of a haplotype-based TDT.

Sample reproducibility of genetic association using different multimarker TDTs in genome-wide association studies: characterization and a new approach

Multimarker Transmission/Disequilibrium Tests (TDTs) are very robust association tests to population admixture and structure which may be used to identify susceptibility loci in genome-wide association studies. Multimarker TDTs using several markers may increase power by capturing high-degree associations. However, there is also a risk of spurious associations and power reduction due to the increase in degrees of freedom. In this study we show that associations found by tests built on simple null hypotheses are highly reproducible in a second independent data set regardless the number of markers. As a test exhibiting this feature to its maximum, we introduce the multimarker 2-Groups TDT (mTDT(2G)), a test which under the hypothesis of no linkage, asymptotically follows a χ2 distribution with 1 degree of freedom regardless the number of markers. The statistic requires the division of parental haplotypes into two groups: disease susceptibility and disease protective haplotype groups. We assessed the test behavior by performing an extensive simulation study as well as a real-data study using several data sets of two complex diseases. We show that mTDT(2G) test is highly efficient and it achieves the highest power among all the tests used, even when the null hypothesis is tested in a second independent data set. Therefore, mTDT(2G) turns out to be a very promising multimarker TDT to perform genome-wide searches for disease susceptibility loci that may be used as a preprocessing step in the construction of more accurate genetic models to predict individual susceptibility to complex diseases.

Genetic and functional evaluation of the role of CXCR1 and CXCR2 in susceptibility to visceral leishmaniasis in north-east India

Background

IL8RA and IL8RB, encoded by CXCR1 and CXCR2, are receptors for interleukin (IL)-8 and other CXC chemokines involved in chemotaxis and activation of polymorphonuclear neutrophils (PMN). Variants at CXCR1 and CXCR2 have been associated with susceptibility to cutaneous and mucocutaneous leishmaniasis in Brazil. Here we investigate the role of CXCR1/CXCR2 in visceral leishmaniasis (VL) in India.

Methods

Three single nucleotide polymorphisms (SNPs) (rs4674259, rs2234671, rs3138060) that tag linkage disequilibrium blocks across CXCR1/CXCR2 were genotyped in primary family-based (313 cases; 176 nuclear families; 836 individuals) and replication (941 cases; 992 controls) samples. Family- and population-based analyses were performed to look for association between CXCR1/CXCR2 variants and VL. Quantitative RT/PCR was used to compare CXCR1/CXCR2 expression in mRNA from paired splenic aspirates taken before and after treatment from 19 VL patients.

Results

Family-based analysis using FBAT showed association between VL and SNPs CXCR1_rs2234671 (Z-score = 2.935, P = 0.003) and CXCR1_rs3138060 (Z-score = 2.22, P = 0.026), but not with CXCR2_rs4674259. Logistic regression analysis of the case-control data under an additive model of inheritance showed association between VL and SNPs CXCR2_rs4674259 (OR = 1.15, 95%CI = 1.01-1.31, P = 0.027) and CXCR1_rs3138060 (OR = 1.25, 95%CI = 1.02-1.53, P = 0.028), but not with CXCR1_rs2234671. The 3-locus haplotype T_G_C across these SNPs was shown to be the risk haplotype in both family- (TRANSMIT; P = 0.014) and population- (OR = 1.16, P = 0.028) samples (combined P = 0.002). CXCR2, but not CXCR1, expression was down regulated in pre-treatment compared to post-treatment splenic aspirates (P = 0.021).

Conclusions

This well-powered primary and replication genetic study, together with functional analysis of gene expression, implicate CXCR2 in determining outcome of VL in India.

Simple strategies for haplotype analysis of the X chromosome with application to age-related macular degeneration

For haplotype analysis of the X chromosome, haplotype-sharing (HS) statistics with sliding windows are defined for males and females separately, which are then combined to a single HS test for the X chromosome. When independent replication samples are not available, the training-testing sets approach is used to validate this procedure and a permutation method is used to obtain its P-value. We applied this method to the X chromosome (with 1804 SNPs) for age-related macular degeneration (AMD). We found a window of five SNPs over a 272 kb region associated with AMD after Bonferroni correction. An examination of the odds ratio and the population attributable risks revealed a disease-preventive haplotype, ATGAC, on these five SNPs. For elderly females without this haplotype, the likelihood of AMD is increased by a factor of 4.75 with a 95% confidence interval (1.43, 15.82). The frequency of ATGAC in HapMap CEU is 0.276. These five SNPs are covered by the gene DIAPH2, which is known to cause premature ovarian failure (POF) in females. Our results indicated that DIAPH2 may be a polygenic pleiotropy for POF and AMD.

Weak association of the platelet-derived growth factor beta (PDGFB) and PDGF receptor beta (PDGFRB) genes with schizophrenia and schizoaffective disorder

Schizophrenia is a severe neuropsychiatric disorder with diverse characterization of symptoms. Extensive research has been performed to elucidate the etiology of schizophrenia. One of the most convincing hypotheses comes from the dopaminergic system although none of the core genes has been consistently positive in association studies.

OBJECTIVE

In this investigation, we explored the possibility that the genes for platelet-derived growth factor beta (PDGFB) and its receptor (PDGFRB) might play an important role in the development of schizophrenia based on previous reports pointing to their ability to interact with the dopamine D(2)/D(4) and NMDA receptors as well as their role in neurite outgrowth.

METHODS

We investigated the association of variants around these genes with schizophrenia and schizoaffective disorder in 104 small nuclear families using the Sib-Transmission Disequilibrium Test (TDT-STDT). Furthermore, quantitative trait analysis using family-based association test was applied to determine possible association of age at onset (AAO).

RESULTS

Allele G in PDGFRB(rs758588) was associated with AAO (P=0.019). An over-transmission of allele T in PDGFB(rs130650) polymorphism (P=0.043) and an over-transmission of allele A in PDGFRB(rs6865659) polymorphism (P=0.046) were observed. Furthermore, the combined TDT-STDT yielded consistent results.

CONCLUSION

Overall, PDGFB and PDGFRB genes might play a role in the etiology of schizophrenia.

The NADH-ubiquinone oxidoreductase 1 alpha subcomplex 5 (NDUFA5) gene variants are associated with autism

OBJECTIVE

Autism appears to have a strong genetic component. The product of the NADH-ubiquinone oxidoreductase 1 alpha subcomplex 5 (NDUFA5) gene is included in the mitochondrial electron transport chain.

METHOD

We performed a case-control study of 235 patients with autism and 214 controls and examined three single-nucleotide polymorphisms (SNPs) within this gene in a Japanese population. We then conducted a transmission disequilibrium test (TDT) analysis in 148 autistic trios.

RESULTS

In the case-control study, two SNPs (rs12666974 and rs3779262) showed a significant association with autism (P=0.00064 and 0.00046 respectively). Furthermore, a haplotype containing these two SNPs showed a significant association (P-global=0.0013, individual haplotype A-A: P=0.010). In TDT analysis, the global and A-A haplotype P-values also indicated significant associations. Minor allele and genotype frequencies were decreased in the autistic subjects.

CONCLUSION

We found significant association between the NDFA5 gene and autism.

Genetic variation in the KIAA0319 5' region as a possible contributor to dyslexia

Reading disabilities (RD) have been linked and associated with markers on chromosome 6p with results from multiple independent samples pointing to KIAA0319 as a risk gene and specifically, the 5' region of this gene. Here we focus genetic studies on a 2.3 kb region spanning the predicted promoter, the first untranslated exon, and part of the first intron, a region we identified as a region of open chromatin. Using DNA from probands with RD, we screened for genetic variants and tested select variants for association. We identified 17 DNA variants in this sample of probands, 16 of which were previously reported in public databases and one previously identified in a screen of this region. Based on the allele frequencies in the probands compared to public databases, and on possible functional consequences of the variation, we selected seven variants to test for association in a sample of families with RD, in addition to four variants which had been tested previously. We also tested two markers 5' of this region that were previously reported as associated. The strongest evidence for association was observed with alleles of the microsatellite marker located in the first untranslated exon and haplotypes of that marker. These results support previous studies indicating the 5' region of the KIAA0319 gene as the location of risk alleles contributing to RD.

Distribution of the number of false discoveries in large-scale family-based association testing with application to the association between PTPN1 and hypertension and obesity

We present a model-free approach to the study of the number of false discoveries for large-scale simultaneous family-based association tests (FBATs) in which the set of discoveries is decided by applying a threshold to the test statistics. When the association between a set of markers in a candidate gene and a group of phenotypes is studied by a class of FBATs, we indicate that a joint null hypothesis distribution for these statistics can be obtained by the fundamental statistical method of conditioning on sufficient statistics for the null hypothesis. Based on the joint null distribution of these statistics, we can obtain the distribution of the number of false discoveries for the set of discoveries defined by a threshold; the size of this set is referred to as its tail count. Simulation studies are presented to demonstrate that the conditional, not the unconditional, distribution of the tail count is appropriate for the study of false discoveries. The usefulness of this approach is illustrated by re-examining the association between PTPN1 and a group of blood-pressure-related phenotypes reported by Olivier et al. (Hum Mol Genet 13:1885-1892, 2004); our results refine and reinforce this association.

HLA-Cw group 1 ligands for KIR increase susceptibility to invasive cervical cancer

Inherited genetic polymorphisms within immune response genes have been shown to associate with risk of invasive cervical cancer (ICC) and its immediate precursor, cervical intraepithelial neoplasia grade 3. Here, we used the transmission/disequilibrium test to detect disease-liability alleles and investigate haplotype transmission of KIR and HLA class I polymorphisms in a large family-based population of women with cervical cancer and their biological parents (359 trios). The effect of distinct human papillomavirus types was also explored. HLA-Cw group 1 (HLA-Cw alleles with asparagine at position 80), which serves as ligand for certain killer immunoglobulin-like receptors (KIR), was significantly overtransmitted in women with ICC (P = 0.04), and particularly in the subgroup of women infected with high risk HPV16 or 18 subtypes (P = 0.008). These data support the involvement of the HLA-C locus in modulating the risk of cervical neoplasia perhaps through its function as ligands for KIR, but functional studies are essential to confirm this hypothesis.

Allergic rhinitis and genetic components: focus on Toll-like receptors (TLRs) gene polymorphism

Allergic rhinitis represents a global health issue affecting 10% to 25% of the population worldwide. Over the years, studies have found that allergic diseases, including allergic rhinitis, are associated with immunological responses to antigens driven by a Th2-mediated immune response. Because Toll-like receptors (TLRs) are involved in both innate and adaptive immune responses to a broad variety of antigens, the association between polymorphisms of TLRs and allergic diseases has been the focus in many animal and human studies. Although the etiology of allergic rhinitis is still unknown, extensive research over the years has confirmed that the underlying causes of allergic diseases are due to many genetic and environmental factors, along with the interactions among them, which include gene-environment, gene-gene, and environment-environment interactions. Currently, there is great inconsistency among studies mainly due to differences in genetic background and unique gene-environment interactions. This paper reviews studies focusing on the association between TLR polymorphisms and allergic diseases, including allergic rhinitis, which would help researchers better understand the role of TLR polymorphisms in the development of allergic rhinitis, and ultimately lead to more efficient therapeutic interventions being developed.

Members 6B and 14 of the TNF receptor superfamily in multiple sclerosis predisposition

TNFRSF6B and TNFRSF14 genes were recently associated with Crohn's disease and rheumatoid arthritis. TNFRSF14 is known as herpes virus entry mediator (HVEM), and herpes viruses have been involved in the aetiology of multiple sclerosis (MS). MS patients present human herpes virus 6 (HHV6) in active plaques and increased antibody responses to HHV6. We aimed to ascertain the role of these genes in MS susceptibility and to investigate the relationship of the gene encoding the widely expressed HVEM receptor with the active replication of HHV6 found in some MS patients. Genotyping of 1370 Spanish MS patients and 1715 ethnically matched controls was performed. HHV6A DNA levels (surrogate of active viral replication) were analysed in serum of MS patients during a 2-year follow-up. Both polymorphisms were associated with MS predisposition, with stronger effect in patients with HHV6 active replication-TNFRSF6B-rs4809330(*)A: P=0.028, OR=1.13; TNFRSF14-rs6684865(*)A: overall P=0.0008, OR=1.2; and HHV6-positive patients vs controls: P=0.017, OR=1.69.

Comparison of multimarker logistic regression models, with application to a genomewide scan of schizophrenia

Background

Genome-wide association studies (GWAS) are a widely used study design for detecting genetic causes of complex diseases. Current studies provide good coverage of common causal SNPs, but not rare ones. A popular method to detect rare causal variants is haplotype testing. A disadvantage of this approach is that many parameters are estimated simultaneously, which can mean a loss of power and slower fitting to large datasets.

Haplotype testing effectively tests both the allele frequencies and the linkage disequilibrium (LD) structure of the data. LD has previously been shown to be mostly attributable to LD between adjacent SNPs. We propose a generalised linear model (GLM) which models the effects of each SNP in a region as well as the statistical interactions between adjacent pairs. This is compared to two other commonly used multimarker GLMs: one with a main-effect parameter for each SNP; one with a parameter for each haplotype.

Results

We show the haplotype model has higher power for rare untyped causal SNPs, the main-effects model has higher power for common untyped causal SNPs, and the proposed model generally has power in between the two others. We show that the relative power of the three methods is dependent on the number of marker haplotypes the causal allele is present on, which depends on the age of the mutation. Except in the case of a common causal variant in high LD with markers, all three multimarker models are superior in power to single-SNP tests.

Including the adjacent statistical interactions results in lower inflation in test statistics when a realistic level of population stratification is present in a dataset.

Using the multimarker models, we analyse data from the Molecular Genetics of Schizophrenia study. The multimarker models find potential associations that are not found by single-SNP tests. However, multimarker models also require stricter control of data quality since biases can have a larger inflationary effect on multimarker test statistics than on single-SNP test statistics.

Conclusions

Analysing a GWAS with multimarker models can yield candidate regions which may contain rare untyped causal variants. This is useful for increasing prior odds of association in future whole-genome sequence analyses.

Linkage of angiotensinogen gene polymorphisms with hypertension in a sibling study of Hong Kong Chinese

OBJECTIVE

The angiotensinogen gene has been linked with human essential hypertension in whites but the relationship in Asian populations has been less consistent. This study aimed to examine genetic associations between hypertension and the M235T, T174M, and G-217A polymorphisms of the angiotensinogen gene in Chinese siblings.

METHODS

We studied members of 126 families with a hypertensive proband, including 434 siblings, of which 178 were hypertensive. Parental history of hypertension was recorded. The M235T, T174M, and G-217A polymorphisms were examined using a microarray method, validated by sequencing. The transmission disequilibrium test was applied to identify whether the genetic polymorphism loci were related to hypertension. Haplotype analysis of the combined polymorphisms was applied using the TRANSMIT program. Linkage study was conducted by applying the affected pedigree member method.

RESULTS

A significant overtransmission was observed for the T235 allele at the M235T polymorphism and hypertension (chi2 = 4.41, P = 0.036) but not for the T174M and G-217A polymorphisms. The haplotype analysis showed a significant association with the haplotypes of paired markers (T174 and T235) with chi2 value of 8.131 (P = 0.004; global test chi2 = 9.131, P = 0.028). Linkage between M235T and hypertension was detected (T = -2.25, P = 0.019), and a tendency for linkage with central obesity-related hypertension was found for the M235T and T174M polymorphisms (P = 0.0087 and P = 0.01).

CONCLUSION

The M235T and T174M variants, especially the T235 allele, contribute to an increased risk of hypertension in these Chinese patients.

Association of the GABRD gene and childhood-onset mood disorders

The chromosome 1p36 region was previously indicated as a locus for susceptibility to recurrent major depressive disorder based on a linkage study in a sample of 497 sib pairs. We investigated the gamma-aminobutyric acid A (GABA(A)) delta receptor subunit gene, GABRD, as a susceptibility gene to childhood-onset mood disorders (COMD) because of substantial evidence implicating GABAergic dysfunction in mood disorders and the position of this gene near the 1p36 linkage region. Using a sample consisting of 645 Hungarian families with a child/adolescent proband diagnosed with a mood disorder with the onset of the first episode before age 15, we found some evidence for the association of two polymorphisms located within the gene, rs2376805 and rs2376803, as well as significant evidence for biased transmission of the haplotypes of these two markers (global chi(2) test for haplotypes = 12.746, 3 df, P = 0.0052). Furthermore, significant evidence of association was only observed in male subjects (n = 438) when the results were analyzed by sex (chi(2) = 9.000 1 df, P = 0.003 for rs2376805). This was in contrast with the previous linkage findings, as LOD scores exceeding 3 were only in female-female pairs in that study. These findings point to the GABRD gene as a susceptibility gene for COMD; however, this gene may not explain the previous linkage finding.

Evidence for associations between the purinergic receptor P2X(7) (P2RX7) and toxoplasmosis

Congenital Toxoplasma gondii infection can result in intracranial calcification, hydrocephalus, and retinochoroiditis. Acquired infection is commonly associated with ocular disease. Pathology is characterized by strong pro-inflammatory responses. Ligation of ATP by purinergic receptor P2X₇, encoded by P2RX7, stimulates pro-inflammatory cytokines and can lead directly to killing of intracellular pathogens. To determine whether P2X₇ plays a role in susceptibility to congenital toxoplasmosis, we examined polymorphisms at P2RX7 in 149 child/parent trios from North America. We found association (FBAT Z scores ±2.429; P= 0.015) between the derived C(+)G(−) allele (f= 0.68; OR= 2.06; 95% CI: 1.14–3.75) at SNP rs1718119 (1068T>C; Thr-348-Ala), and a second synonymous variant rs1621388 in linkage disequilibrium with it, and clinical signs of disease per se. Analysis of clinical sub-groups showed no association with hydrocephalus, with effect sizes for associations with retinal disease and brain calcifications enhanced (OR=3.0 to 4.25; 0.004<P<0.009) when hydrocephalus was removed from the analysis. Association with toxoplasmic retinochoroiditis was replicated (FBAT Z scores ±3.089; P= 0.002) in a small family-based study (60 families; 68 affected offspring) of acquired infection in Brazil, where the ancestral T(+) allele (f= 0.296) at SNP rs1718119 was strongly protective (OR= 0.27; 95% CI: 0.09–0.80). (Words 194)

Polymorphisms in leucine-rich repeat genes are associated with autism spectrum disorder susceptibility in populations of European ancestry

BACKGROUND

Autism spectrum disorders (ASDs) are a group of highly heritable neurodevelopmental disorders which are characteristically comprised of impairments in social interaction, communication and restricted interests/behaviours. Several cell adhesion transmembrane leucine-rich repeat (LRR) proteins are highly expressed in the nervous system and are thought to be key regulators of its development. Here we present an association study analysing the roles of four promising candidate genes - LRRTM1 (2p), LRRTM3 (10q), LRRN1 (3p) and LRRN3 (7q) - in order to identify common genetic risk factors underlying ASDs.

METHODS

In order to gain a better understanding of how the genetic variation within these four gene regions may influence susceptibility to ASDs, a family-based association study was undertaken in 661 families of European ancestry selected from four different ASD cohorts. In addition, a case-control study was undertaken across the four LRR genes, using logistic regression in probands with ASD of each population against 295 ECACC controls.

RESULTS

Significant results were found for LRRN3 and LRRTM3 (P < 0.005), using both single locus and haplotype approaches. These results were further supported by a case-control analysis, which also highlighted additional SNPs in LRRTM3.

CONCLUSIONS

Overall, our findings implicate the neuronal leucine-rich genes LRRN3 and LRRTM3 in ASD susceptibility.

Association of reading disabilities with regions marked by acetylated H3 histones in KIAA0319

Reading disabilities (RDs) have been associated with chromosome 6p with recent studies pointing to two genes, DCDC2 and KIAA0319. In this study, markers across the 6p region were tested for association with RD. Our strongest findings were for association with markers in KIAA0319, although with the opposite alleles compared with a previous study. We also found association with markers in VMP, but not with DCDC2. Current evidence indicates that differential regulation of KIAA0319 and DCDC2 contributes to RD, thus we used chromatin immunoprecipitation coupled with genomic tiling arrays (ChIP-chip) to map acetylated histones, a molecular marker for regulatory elements, across a 500 kb genomic region covering the RD locus on 6p. This approach identified several regions marked by acetylated histones that mapped near associated markers, including intron 7 of DCDC2 and the 5' region of KIAA0319. The latter is located within the 70 kb region previously associated with differential expression of KIAA0319. Interestingly, five markers associated with RD in independent studies were also located within the 2.7 kb acetylated region, and six additional associated markers, including the most significant one in this study, were located within a 22 kb haplotype block that encompassed this region. Our data indicates that this putative regulatory region is a likely site of genetic variation contributing to RD in our sample, further narrowing the candidate region.

TP53, MDM2, NQO1, and susceptibility to cervical cancer

Host genetic variability modifies the risk of cervical cancer in women infected with oncogenic human papillomavirus (HPV). Studies have reported an association of the TP53 codon 72 arginine and cervical cancer, but the results are inconsistent. We examined the association of this single nucleotide polymorphism (SNP) in women with cervical cancer and cervical intraepithelial neoplasia grade 3, using a family-based association test. We further explored SNPs in two genes that regulate p53 stability: MDM2 (SNP309) and NQO1 (SNP609, SNP465). We also examined the relationship between host genotype and tumor HPV type. We genotyped 577 patients and their biological parents and/or siblings, using PCR-RFLP or Taqman assays. HPVs were typed by sequence-based methods. The transmission/disequilibrium test was used to detect disease-susceptibility alleles. The arginine peptide of TP53 codon 72 was overtransmitted in Caucasian families (P = 0.043), and the significance of this finding was enhanced in a subgroup of women infected with HPV16- and/or HPV18-related HPVs (P = 0.026). Allele C of NQO1 SNP609 was also overtransmitted in all cases (P = 0.026). We found no association between MDM2 SNP309 or NQO1 SNP465 and cervical cancer. Our results indicate that functional polymorphisms in TP53 codon 72 and NQO1 SNP609 associate with the risk of cervical cancer especially in women infected with type 16- and/or type 18-related HPVs.

CXCR1 and SLC11A1 polymorphisms affect susceptibility to cutaneous leishmaniasis in Brazil: a case-control and family-based study

Background

L. braziliensis causes cutaneous (CL) and mucosal (ML) leishmaniasis. Wound healing neutrophil (PMN) and macrophage responses made following the bite of the vector sand fly contribute to disease progression in mice. To look at the interplay between PMN and macrophages in disease progression in humans we asked whether polymorphisms at genes that regulate their infiltration or function are associated with different clinical phenotypes. Specifically, CXCR1 (IL8RA) and CXCR2 (IL8RB) are receptors for chemokines that attract PMN to inflammatory sites. They lie 30-260 kb upstream of SLC11A1, a gene known primarily for its role in regulating macrophage activation, resistance to leishmaniasis, and wound healing responses in mice, but also known to be expressed in PMN, macrophages and dendritic cells.

Methods

Polymorphic variants at CXCR1, CXCR2 and SLC11A1 were analysed using Taqman or ABI fragment separation technologies in cases (60 CL; 60 ML), unrelated controls (n = 120), and multicase families (104 nuclear families; 88 ML, 250 CL cases) from Brazil. Logistic regression analysis, family-based association testing (FBAT) and haplotype analysis (TRANSMIT) were performed.

Results

Case-control analysis showed association between the common C allele (OR 2.38; 95% CI 1.23-4.57; P = 0.009) of CXCR1_rs2854386 and CL, supported by family-based (FBAT; Z score 2.002; P = 0.045) analysis (104 nuclear families; 88 ML, 250 CL cases). ML associated with the rarer G allele (Z score 1.999; P = 0.046). CL associated with a 3' insertion/deletion polymorphism at SLC11A1 (Z score 2.549; P = 0.011).

Conclusions

The study supports roles for CXCR1 and SLC11A1 in the outcome of L. braziliensis infection in humans. Slc11a1 does not influence cutaneous lesion development following needle injection of Leishmania in mice, suggesting that its role here might relate to the action of PMN, macrophage and/or dendritic cells in the wound healing response to the sand fly bite. Together with the CXCR1 association, the data are consistent with hypotheses relating to the possible role of PMN in initiation of a lesion following the delivery of parasites via the sand fly bite. Association of ML with the rare derived G allele suggests that PMN also have an important positive role to play in preventing this form of the disease.

Two-stage testing strategies for genome-wide association studies in family-based designs

The analysis of genome-wide association studies (GWAS) poses statistical hurdles that have to be handled efficiently in order for the study to be successful. The two largest impediments in the analysis phase of the study are the multiple comparisons problem and maintaining robustness against confounding due to population admixture and stratification. For quantitative traits in family-based designs, Van Steen (1) proposed a two-stage testing strategy that can be considered a hybrid approach between family-based and population-based analysis. By including the population-based component into the family-based analysis, the Van Steen algorithm maximizes the statistical power, while at the same time, maintains the original robustness of family-based association tests (FBATs) (2-4). The Van Steen approach consists of two statistically independent steps, a screening step and a testing step. For all genotyped single nucleotide polymorphisms (SNPs), the screening step examines the evidence for association at a population-based level. Based on support for a potential genetic association from the screening step, the SNPs are prioritized for testing in the next step, where they are analyzed with a FBAT (3). By exploiting population-based information in the screening step that is not utilized in family-based association testing step, the two steps are statistically independent. Therefore, the use of the population-based data for the purposes of screening does not bias the FBAT statistic calculated in the testing step. Depending on the trait type and the ascertainment conditions, Van Steen-type testing strategies can achieve statistical power levels that are comparable to those of population-based studies with the same number of probands. In this chapter, we review the original Van Steen algorithm, its numerous extensions, and discuss its advantages and disadvantages.

Association of the alpha4 integrin subunit gene (ITGA4) with autism

In the present work, we provide further evidence for the involvement of the integrin alpha-4 precursor gene (ITGA4) in the etiology of autism, by replicating previous findings of a genetic association with autism in various independent populations. The ITGA4 gene maps to the autism linkage region on 2q31-33 and is therefore a plausible positional candidate. We tested eight single nucleotide polymorphisms (SNPs) in the ITGA4 gene region for association with autism in a sample of 164 nuclear families. Evidence for association was found for the rs155100 marker (P = 0.019) and for a number of specific marker haplotypes containing this SNP (0.00053 < P < 0.022). alpha4 integrins are known to play a key role in neuroinflammatory processes, which are hypothesized to contribute to autism. In this study, an association was found between the ITGA4 rs1449263 marker and levels of a serum autoantibody directed to brain tissue, which was previously shown to be significantly more frequent in autistic patients than in age-matched controls in our population. This result suggests that the ITGA4 gene could be involved in a neuroimmune process thought to occur in autistic patients and, together with previous findings, offers a new perspective on the role of integrins in the etiology of autism to which little attention has been paid so far.

A data-driven weighting scheme for family-based genome-wide association studies

Recently, Steen et al proposed a novel two-stage approach for family-based genome-wide association studies. In the first stage, a test based on between-family information is used to rank SNPs according to their P-values or conditional power of the test. In the second stage, the R most promising SNPs are tested using a family-based association test. We call this two-stage approach top R method. Ionita-Laza et al proposed an exponential weighting method within a two-stage framework. In the second stage of this approach, instead of testing top R SNPs, it tests all SNPs and weights the P-values of association test according to the information of the first stage. However, both of the top R and exponential weighting methods only use the information from the first stage to rank SNPs. It seems that the two methods do not use information from the first stage efficiently. Furthermore, it may be unreasonable for the exponential weighting method to use the same weight for all SNPs within a group when only one or a few SNPs are related with a disease. In this article, we propose a data-driven weighting scheme within a two-stage framework. In this method, we use the information from the first stage to determine a SNP-specific weight for each SNP. We use simulation studies to evaluate the performance of our method. The simulation results showed that our proposed method is consistently more powerful than the top R method and the exponential weighting method, regardless of the LD structure, population structure, and family structure.

A variable-sized sliding-window approach for genetic association studies via principal component analysis

Recently with the rapid improvements in high-throughout genotyping techniques, researchers are facing the very challenging task of analysing large-scale genetic associations, especially at the whole-genome level, without an optimal solution. In this study, we propose a new approach for genetic association analysis that is based on a variable-sized sliding-window framework and employs principal component analysis to find the optimum window size. With the help of the bisection algorithm in window-size searching, our method is more computationally efficient than available approaches. We evaluate the performance of the proposed method by comparing it with two other methods-a single-marker method and a variable-length Markov chain method. We demonstrate that, in most cases, the proposed method out-performs the other two methods. Furthermore, since the proposed method is based on genotype data, it does not require any computationally intensive phasing program to account for uncertain haplotype phase.

Association of attention-deficit/hyperactivity disorder with a candidate region for reading disabilities on chromosome 6p

BACKGROUND

Reading disabilities (RD) and attention-deficit hyperactivity/disorder (ADHD) are two common childhood disorders that co-occur by chance more often than expected. Twin studies and overlapping genetic linkage findings indicate that shared genetic factors partially contribute to this comorbidity. Linkage of ADHD to 6p, an identified RD candidate locus, has previously been reported, suggesting the possibility of a pleiotropic gene at this locus. RD has been previously associated with five genes in the region, particularly DCDC2 and KIAA0319.

METHODS

To test whether these genes also contribute to ADHD, we investigated markers previously associated with RD for association with ADHD and ADHD symptoms in a sample of families with ADHD (n = 264). Markers were located in two subregions, VMP/DCDC2 and KIAA0319/TTRAP.

RESULTS

Across all analyses conducted, strong evidence for association was observed in the VMP/DCDC2 region. Association was equally strong with symptoms of both inattention and hyperactivity/impulsivity, suggesting that this locus contributes to both symptom dimensions. Markers were also tested for association with measures of reading skills (word identification, decoding); however, there was virtually no overlap in the markers associated with ADHD and those associated with reading skills in this sample.

CONCLUSIONS

Overall this study supports a previous linkage study of ADHD indicating a risk gene for ADHD on 6p and points to VMP or DCDC2 as the most likely candidates.

HLA haplotypes in recurrent aphthous stomatitis: a mode of inheritance?

The aim of this study was to investigate the genetic association between recurrent aphthous stomatitis (RAS) and human leucocyte antigen (HLA) class I and II alleles and HLA haplotypes. Families selected had at least one child suffering from recurrent aphthous stomatitis in addition to one or both of the parents. HLA-A, -B and -DR alleles were typed in 29 families, 27 nuclear and two extended (121 subjects). HLA haplotypes of all family members with RAS were compared with those who were RAS negative. Although major histocompatibility complex class I and II gene analysis failed to demonstrate any significant association between RAS and HLA antigens, the study of HLA haplotypes revealed a significant association between HLA haplotypes and susceptibility to RAS. The results indicate that susceptibility to RAS segregates in families in association with HLA haplotypes.

The KIAA0319-like (KIAA0319L) gene on chromosome 1p34 as a candidate for reading disabilities

A locus on chromosome 1p34-36 (DYX8) has been linked to developmental dyslexia or reading disabilities (RD) in three independent samples. In the current study, we investigated a candidate gene KIAA0319-Like (KIAA0319L) within DYX8, as it is homologous to KIAA0319, a strong RD candidate gene on chromosome 6p (DYX2). Association was assessed by using five tagging single nucleotide polymorphisms in a sample of 291 nuclear families ascertained through a proband with reading difficulties. Evidence of association was found for a single marker (rs7523017; P=0.042) and a haplotype (P=0.031), with RD defined as a categorical trait in a subset of the sample (n=156 families) with a proband that made our criteria for RD. The same haplotype also showed evidence for association with quantitative measures of word-reading efficiency (i.e., a composite score of word identification and decoding; P=0.032) and rapid naming of objects and colors (P=0.047) when analyzed using the entire sample. Although the results from the current study are modestly significant and would not withstand a correction for multiple testing, KIAA0319L remains an intriguing positional and functional candidate for RD, especially when considered alongside the supporting evidence for its homolog KIAA0319 on chromosome 6p. Additional studies in independent samples are now required to confirm these findings.

Application of homozygosity haplotype analysis to genetic mapping with high-density SNP genotype data

BACKGROUND

In families segregating a monogenic genetic disorder with a single disease gene introduction, patients share a mutation-carrying chromosomal interval with identity-by-descent (IBD). Such a shared chromosomal interval or haplotype, surrounding the actual pathogenic mutation, is typically detected and defined by multipoint linkage and phased haplotype analysis using microsatellite or SNP genotype data. High-density SNP genotype data presents a computational challenge for conventional genetic analyses. A novel non-parametric method termed Homozygosity Haplotype (HH) was recently proposed for the genome-wide search of the autosomal segments shared among patients using high density SNP genotype data.

METHODOLOGY/PRINCIPAL FINDINGS

The applicability and the effectiveness of HH in identifying the potential linkage of disease causative gene with high-density SNP genotype data were studied with a series of monogenic disorders ascertained in eastern Canadian populations. The HH approach was validated using the genotypes of patients from a family affected with a rare autosomal dominant disease Schnyder crystalline corneal dystrophy. HH accurately detected the approximately 1 Mb genomic interval encompassing the causative gene UBIAD1 using the genotypes of only four affected subjects. The successful application of HH to identify the potential linkage for a family with pericentral retinal disorder indicates that HH can be applied to perform family-based association analysis by treating affected and unaffected family members as cases and controls respectively. A new strategy for the genome-wide screening of known causative genes or loci with HH was proposed, as shown the applications to a myoclonus dystonia and a renal failure cohort.

CONCLUSIONS/SIGNIFICANCE

Our study of the HH approach demonstrates that HH is very efficient and effective in identifying potential disease linked region. HH has the potential to be used as an efficient alternative approach to sequencing or microsatellite-based fine mapping for screening the known causative genes in genetic disease study.