Joint association analysis of bivariate quantitative and qualitative traits

Univariate genome-wide association analysis of quantitative and qualitative traits has been investigated extensively in the literature. In the presence of correlated phenotypes, it is more intuitive to analyze all phenotypes simultaneously. We describe an efficient likelihood-based approach for the joint association analysis of quantitative and qualitative traits in unrelated individuals. We assume a probit model for the qualitative trait, under which an unobserved latent variable and a prespecified threshold determine the value of the qualitative trait. To jointly model the quantitative and qualitative traits, we assume that the quantitative trait and the latent variable follow a bivariate normal distribution. The latent variable is allowed to be correlated with the quantitative phenotype. Simultaneous modeling of the quantitative and qualitative traits allows us to make more precise inference on the pleiotropic genetic effects. We derive likelihood ratio tests for the testing of genetic effects. An application to the Genetic Analysis Workshop 17 data is provided. The new method yields reasonable power and meaningful results for the joint association analysis of the quantitative trait Q1 and the qualitative trait disease status at SNPs with not too small MAF

A question of scale: Human migrations writ large and small

Several recent papers illustrate the importance of migration and gene flow in molding the patterns of genetic variation observed in humans today. We place the varied demographic processes covered by these terms into a more general framework, and discuss some of the challenges facing attempts to reconstruct past human mobility and determine its influence on our genetic heritage

Identity-by-descent estimation with population- and pedigree-based imputation in admixed family data

BACKGROUND: In the past few years, imputation approaches have been mainly used in population-based designs of genome-wide association studies, although both family- and population-based imputation methods have been proposed. With the recent surge of family-based designs, family-based imputation has become more important. Imputation methods for both designs are based on identity-by-descent (IBD) information. Apart from imputation, the use of IBD information is also common for several types of genetic analysis, including pedigree-based linkage analysis. METHODS: We compared the performance of several family- and population-based imputation methods in large pedigrees provided by Genetic Analysis Workshop 19 (GAW19). We also evaluated the performance of a new IBD mapping approach that we propose, which combines IBD information from known pedigrees with information from unrelated individuals. RESULTS: Different combinations of the imputation methods have varied imputation accuracies. Moreover, we showed gains from the use of both known pedigrees and unrelated individuals with our IBD mapping approach over the use of known pedigrees only. CONCLUSIONS: Our results represent accuracies of different combinations of imputation methods that may be useful for data sets similar to the GAW19 pedigree data. Our IBD mapping approach, which uses both known pedigree and unrelated individuals, performed better than classical linkage analysis

Genome-wide linkage analysis of 972 bipolar pedigrees using single-nucleotide polymorphisms.

Because of the high costs associated with ascertainment of families, most linkage studies of Bipolar I disorder (BPI) have used relatively small samples. Moreover, the genetic information content reported in most studies has been less than 0.6. Although microsatellite markers spaced every 10 cM typically extract most of the genetic information content for larger multiplex families, they can be less informative for smaller pedigrees especially for affected sib pair kindreds. For these reasons we collaborated to pool family resources and carried out higher density genotyping. Approximately 1100 pedigrees of European ancestry were initially selected for study and were genotyped by the Center for Inherited Disease Research using the Illumina Linkage Panel 12 set of 6090 single-nucleotide polymorphisms. Of the ~1100 families, 972 were informative for further analyses, and mean information content was 0.86 after pruning for linkage disequilibrium. The 972 kindreds include 2284 cases of BPI disorder, 498 individuals with bipolar II disorder (BPII) and 702 subjects with recurrent major depression. Three affection status models (ASMs) were considered: ASM1 (BPI and schizoaffective disorder, BP cases (SABP) only), ASM2 (ASM1 cases plus BPII) and ASM3 (ASM2 cases plus recurrent major depression). Both parametric and non-parametric linkage methods were carried out. The strongest findings occurred at 6q21 (non-parametric pairs LOD 3.4 for rs1046943 at 119 cM) and 9q21 (non-parametric pairs logarithm of odds (LOD) 3.4 for rs722642 at 78 cM) using only BPI and schizoaffective (SA), BP cases. Both results met genome-wide significant criteria, although neither was significant after correction for multiple analyses. We also inspected parametric scores for the larger multiplex families to identify possible rare susceptibility loci. In this analysis, we observed 59 parametric LODs of 2 or greater, many of which are likely to be close to maximum possible scores. Although some linkage findings may be false positives, the results could help prioritize the search for rare variants using whole exome or genome sequencing

Clinical and Biomarker Changes in Dominantly Inherited Alzheimer's Disease

BACKGROUND: The order and magnitude of pathologic processes in Alzheimer's disease are not well understood, partly because the disease develops over many years. Autosomal dominant Alzheimer's disease has a predictable age at onset and provides an opportunity to determine the sequence and magnitude of pathologic changes that culminate in symptomatic disease. METHODS: In this prospective, longitudinal study, we analyzed data from 128 participants who underwent baseline clinical and cognitive assessments, brain imaging, and cerebrospinal fluid (CSF) and blood tests. We used the participant's age at baseline assessment and the parent's age at the onset of symptoms of Alzheimer's disease to calculate the estimated years from expected symptom onset (age of the participant minus parent's age at symptom onset). We conducted cross-sectional analyses of baseline data in relation to estimated years from expected symptom onset in order to determine the relative order and magnitude of pathophysiological changes. RESULTS: Concentrations of amyloid-beta (Aβ)(42) in the CSF appeared to decline 25 years before expected symptom onset. Aβ deposition, as measured by positron-emission tomography with the use of Pittsburgh compound B, was detected 15 years before expected symptom onset. Increased concentrations of tau protein in the CSF and an increase in brain atrophy were detected 15 years before expected symptom onset. Cerebral hypometabolism and impaired episodic memory were observed 10 years before expected symptom onset. Global cognitive impairment, as measured by the Mini-Mental State Examination and the Clinical Dementia Rating scale, was detected 5 years before expected symptom onset, and patients met diagnostic criteria for dementia at an average of 3 years after expected symptom onset. CONCLUSIONS: We found that autosomal dominant Alzheimer's disease was associated with a series of pathophysiological changes over decades in CSF biochemical markers of Alzheimer's disease, brain amyloid deposition, and brain metabolism as well as progressive cognitive impairment. Our results require confirmation with the use of longitudinal data and may not apply to patients with sporadic Alzheimer's disease. (Funded by the National Institute on Aging and others; DIAN ClinicalTrials.gov number, NCT00869817.)

Empirical vs Bayesian approach for estimating haplotypes from genotypes of unrelated individuals

BACKGROUND: The completion of the HapMap project has stimulated further development of haplotype-based methodologies for disease associations. A key aspect of such development is the statistical inference of individual diplotypes from unphased genotypes. Several methodologies for inferring haplotypes have been developed, but they have not been evaluated extensively to determine which method not only performs well, but also can be easily incorporated in downstream haplotype-based association analyses. In this paper, we attempt to do so. Our evaluation was carried out by comparing the two leading Bayesian methods, implemented in PHASE and HAPLOTYPER, and the two leading empirical methods, implemented in PL-EM and HPlus. We used these methods to analyze real data, namely the dense genotypes on X-chromosome of 30 European and 30 African trios provided by the International HapMap Project, and simulated genotype data. Our conclusions are based on these analyses. RESULTS: All programs performed very well on X-chromosome data, with an average similarity index of 0.99 and an average prediction rate of 0.99 for both European and African trios. On simulated data with approximation of coalescence, PHASE implementing the Bayesian method based on the coalescence approximation outperformed other programs on small sample sizes. When the sample size increased, other programs performed as well as PHASE. PL-EM and HPlus implementing empirical methods required much less running time than the programs implementing the Bayesian methods. They required only one hundredth or thousandth of the running time required by PHASE, particularly when analyzing large sample sizes and large umber of SNPs. CONCLUSION: For large sample sizes (hundreds or more), which most association studies require, the two empirical methods might be used since they infer the haplotypes as accurately as any Bayesian methods and can be incorporated easily into downstream haplotype-based analyses such as haplotype-association analyses

LDLR Expression and Localization Are Altered in Mouse and Human Cell Culture Models of Alzheimer's Disease

Alzheimer's disease (AD) is a chronic neurodegenerative disorder and the most common form of dementia. The major molecular risk factor for late-onset AD is expression of the ε-4 allele of apolipoprotein E (apoE), the major cholesterol transporter in the brain. The low-density lipoprotein receptor (LDLR) has the highest affinity for apoE and plays an important role in brain cholesterol metabolism.Using RT-PCR and western blotting techniques we found that over-expression of APP caused increases in both LDLR mRNA and protein levels in APP transfected H4 neuroglioma cells compared to H4 controls. Furthermore, immunohistochemical experiments showed aberrant localization of LDLR in H4-APP neuroglioma cells, Aβ-treated primary neurons, and in the PSAPP transgenic mouse model of AD. Finally, immunofluorescent staining of LDLR and of γ- and α-tubulin showed a change in LDLR localization preferentially away from the plasma membrane that was paralleled by and likely the result of a disruption of the microtubule-organizing center and associated microtubule network.These data suggest that increased APP expression and Aβ exposure alters microtubule function, leading to reduced transport of LDLR to the plasma membrane. Consequent deleterious effects on apoE uptake and function will have implications for AD pathogenesis and/or progression

Haplotype association analyses in resources of mixed structure using Monte Carlo testing

<p>Abstract</p> <p>Background</p> <p>Genomewide association studies have resulted in a great many genomic regions that are likely to harbor disease genes. Thorough interrogation of these specific regions is the logical next step, including regional haplotype studies to identify risk haplotypes upon which the underlying critical variants lie. Pedigrees ascertained for disease can be powerful for genetic analysis due to the cases being enriched for genetic disease. Here we present a Monte Carlo based method to perform haplotype association analysis. Our method, hapMC, allows for the analysis of full-length and sub-haplotypes, including imputation of missing data, in resources of nuclear families, general pedigrees, case-control data or mixtures thereof. Both traditional association statistics and transmission/disequilibrium statistics can be performed. The method includes a phasing algorithm that can be used in large pedigrees and optional use of pseudocontrols.</p> <p>Results</p> <p>Our new phasing algorithm substantially outperformed the standard expectation-maximization algorithm that is ignorant of pedigree structure, and hence is preferable for resources that include pedigree structure. Through simulation we show that our Monte Carlo procedure maintains the correct type 1 error rates for all resource types. Power comparisons suggest that transmission-disequilibrium statistics are superior for performing association in resources of only nuclear families. For mixed structure resources, however, the newly implemented pseudocontrol approach appears to be the best choice. Results also indicated the value of large high-risk pedigrees for association analysis, which, in the simulations considered, were comparable in power to case-control resources of the same sample size.</p> <p>Conclusions</p> <p>We propose hapMC as a valuable new tool to perform haplotype association analyses, particularly for resources of mixed structure. The availability of meta-association and haplotype-mining modules in our suite of Monte Carlo haplotype procedures adds further value to the approach.</p