Multiple testing correction in linear mixed models.

BackgroundMultiple hypothesis testing is a major issue in genome-wide association studies (GWAS), which often analyze millions of markers. The permutation test is considered to be the gold standard in multiple testing correction as it accurately takes into account the correlation structure of the genome. Recently, the linear mixed model (LMM) has become the standard practice in GWAS, addressing issues of population structure and insufficient power. However, none of the current multiple testing approaches are applicable to LMM.ResultsWe were able to estimate per-marker thresholds as accurately as the gold standard approach in real and simulated datasets, while reducing the time required from months to hours. We applied our approach to mouse, yeast, and human datasets to demonstrate the accuracy and efficiency of our approach.ConclusionsWe provide an efficient and accurate multiple testing correction approach for linear mixed models. We further provide an intuition about the relationships between per-marker threshold, genetic relatedness, and heritability, based on our observations in real data

Linkage and linkage disequilibrium analysis of the lipoprotein lipase gene with lipid profiles in Chinese hypertensive families

A B S T R A C T Elevated TG [triacylglycerol (triglyceride)] is a significant independent risk factor for cardiovascular disease. LPL (lipoprotein lipase) is one of the key enzymes in the metabolism of the TG-rich lipoproteins which hydrolyses TG from the chylomicrons and very-LDL (low-density lipoprotein). To investigate the relationship between the LPL gene and lipid profiles, especially TG, in 148 hypertensive families, we have chosen seven flanking microsatellite markers and four internal markers of the LPL gene and conducted linkage analysis by SOLAR and S.A.G.E. (statistical analysis for genetic epidemiology)/SIBPAL 2 programs, and linkage disequilibrium analysis by QTDT (quantitative transmission/disequilibrium test) and GOLD (graphical overview of linkage disequilibrium). There were statistically significant differences in lipid levels between subjects without and with hypertension within families. A maximum LOD score of 1.3 with TG at the marker D8S261 was observed by SOLAR. Using S.A.G.E./SIBPAL 2, we identified a linkage with TG at the marker 'ATTT' located within intron 6 of the LPL gene (P = 0.0095). Two SNPs (single nucleotide polymorphisms), HindIII and HinfI, were found in linkage disequilibrium with LDLcholesterol levels (P = 0.0178 and P = 0.0088 respectively). A strong linkage disequilibrium was observed between the HindIII in intron 8 and HinfI in the exon 9 (P < 0.00001, D = 0.895). Linkage disequilibrium was also found between the 'ATTT' polymorphism in intron 6 and two SNPs (P = 0.0021 and D = 0.611 for HindIII; and P = 0.00004, D = 0.459 for HinfI). The present study in the Chinese families with hypertension suggested that the LPL gene might influence lipid levels, especially TG metabolism. Replication studies both in Chinese and other populations are warranted to confirm these results

Role of estrogen related receptor beta (ESRRB) in DFN35B hearing impairment and dental decay

BACKGROUND: Congenital forms of hearing impairment can be caused by mutations in the estrogen related receptor beta (ESRRB) gene. Our initial linkage studies suggested the ESRRB locus is linked to high caries experience in humans. METHODS: We tested for association between the ESRRB locus and dental caries in 1,731 subjects, if ESRRB was expressed in whole saliva, if ESRRB was associated with the microhardness of the dental enamel, and if ESRRB was expressed during enamel development of mice. RESULTS: Two families with recessive ESRRB mutations and DFNB35 hearing impairment showed more extensive dental destruction by caries. Expression levels of ESRRB in whole saliva samples showed differences depending on sex and dental caries experience. CONCLUSIONS: The common etiology of dental caries and hearing impairment provides a venue to assist in the identification of individuals at risk to either condition and provides options for the development of new caries prevention strategies, if the associated ESRRB genetic variants are correlated with efficacy.Fil: Weber, Megan L.. University of Pittsburgh; Estados UnidosFil: Hsin, Hong Yuan. University of Pittsburgh; Estados UnidosFil: Kalay, Ersan. Karadeniz Technical University; TurquíaFil: Brožková, Dana Š. Charles University; República Checa. University Hospital Motol; República ChecaFil: Shimizu, Takehiko. Nihon University. School of Dentistry; JapónFil: Bayram, Merve. Medipol Istanbul University; TurquíaFil: Deeley, Kathleen. University of Pittsburgh; Estados UnidosFil: Küchler, Erika C.. University of Pittsburgh; Estados UnidosFil: Forella, Jessalyn. University of Pittsburgh; Estados UnidosFil: Ruff, Timothy D.. University of Pittsburgh; Estados UnidosFil: Trombetta, Vanessa M.. University of Pittsburgh; Estados UnidosFil: Sencak, Regina C.. University of Pittsburgh; Estados UnidosFil: Hummel, Michael. University of Pittsburgh; Estados UnidosFil: Briseño Ruiz, Jessica. University of Pittsburgh; Estados UnidosFil: Revu, Shankar K.. University of Pittsburgh; Estados UnidosFil: Granjeiro, José M.. Universidade Federal Fluminense; BrasilFil: Antunes, Leonardo S.. Universidade Federal Fluminense; BrasilFil: Antunes, Livia A.. Universidade Federal Fluminense; BrasilFil: Abreu, Fernanda V.. Universidade Federal Fluminense; BrasilFil: Costabel, Marcelo C.. Universidade Federal do Rio de Janeiro; BrasilFil: Tannure, Patricia N.. Veiga de Almeida University; Brasil. Salgado de Oliveira University; BrasilFil: Koruyucu, Mine. Istanbul University; TurquíaFil: Patir, Asli. Medipol Istanbul University; TurquíaFil: Poletta, Fernando Adrián. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Parque Centenario. CEMIC-CONICET. Centro de Educaciones Médicas e Investigaciones Clínicas ; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas; ArgentinaFil: Mereb, Juan C.. Estudio Colaborativo Latino Americano de Malformaciones Congénitas; ArgentinaFil: Castilla, Eduardo Enrique. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Parque Centenario. CEMIC-CONICET. Centro de Educaciones Médicas e Investigaciones Clínicas ; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas; ArgentinaFil: Orioli, Iêda M.. Universidade Federal do Rio de Janeiro; BrasilFil: Marazita, Mary L.. University of Pittsburgh; Estados UnidosFil: Ouyang, Hongjiao. University of Pittsburgh; Estados UnidosFil: Jayaraman, Thottala. University of Pittsburgh; Estados UnidosFil: Seymen, Figen. Istanbul University; TurquíaFil: Vieira, Alexandre R.. University of Pittsburgh; Estados Unido

Molecular evolution of the sheep prion protein gene

Transmissible spongiform encephalopathies (TSEs) are infectious, fatal neurodegenerative diseases characterized by aggregates of modified forms of the prion protein (PrP) in the central nervous system. Well known examples include variant Creutzfeldt-Jakob Disease (vCJD) in humans, BSE in cattle, chronic wasting disease in deer and scrapie in sheep and goats. In humans, sheep and deer, disease susceptibility is determined by host genotype at the prion protein gene (PRNP). Here I examine the molecular evolution of PRNP in ruminants and show that variation in sheep appears to have been maintained by balancing selection, a profoundly different process from that seen in other ruminants. Scrapie eradication programs such as those recently implemented in the UK, USA and elsewhere are based on the assumption that PRNP is under positive selection in response to scrapie. If, as these data suggest, that assumption is wrong, eradication programs will disrupt this balancing selection, and may have a negative impact on the fitness or scrapie resistance of national flocks

JLIN: A java based linkage disequilibrium plotter

BACKGROUND: A great deal of effort and expense are being expended internationally in attempts to detect genetic polymorphisms contributing to susceptibility to complex human disease. Techniques such as Linkage Disequilibrium mapping are being increasingly used to examine and compare markers across increasingly large datasets. Visualisation techniques are becoming essential to analyse the ever-growing volume of data and results available with any given analysis. RESULTS: JLIN (Java LINkage disequilibrium plotter) is a software package designed for customisable, intuitive visualisation of Linkage Disequilibrium (LD) across all common computing platforms. Customisation allows the user to choose particular visualisations, statistical measures and measurement ranges. JLIN also allows the user to export images of the LD visualisation in several common document formats. CONCLUSION: JLIN allows the user to visually compare and contrast the results of a range of statistical measures on the input dataset(s). These measures include the commonly used D' and r(2 )statistics and empirical p-values. JLIN has a number of unique and novel features that improve on existing LD visualisation tools

SNPAnalyzer: a web-based integrated workbench for single-nucleotide polymorphism analysis

SNPAnalyzer is a software that performs four essential statistical analyses of SNPs in a common computational environment. It is composed of three main modules: (i) data manipulation, (ii) analysis and (iii) visualization. The data manipulation module is responsible for data input and output, and handles genotype, phenotype and genetic distance data. To ensure user convenience, the data format is simple. The analysis module performs statistical calculations and consists of four subcomponents: (i) Hardy–Weinberg equilibrium, (ii) Haplotype Estimation, (iii) linkage disequilibrium (LD) and (iv) quantitative trait locus analysis. The main feature of the analysis module is multiple implementations of different algorithms and indices for haplotype estimation and for LD analysis. This enables users to compare separate results generated by different algorithms, which help to avoid biased results acquired by applying a single statistical algorithm. The performance of all implemented algorithms has been validated using experimentally proven datasets. The visualization module presents most of the analyzed results as figures, rather than as simple text, which aids in the intuitive understanding of complex data. The SNPAnalyzer has been developed using C and C++ and is available at

Haplotype reconstruction error as a classical misclassification problem

Statistically reconstructing haplotypes from single nucleotide polymorphism (SNP) genotypes, can lead to falsely classified haplotypes. This can be an issue when interpreting haplotype association results or when selecting subjects with certain haplotypes for subsequent functional studies. It was our aim to quantify haplotype reconstruction error and to provide tools for it. By numerous simulation scenarios, we systematically investigated several error measures, including discrepancy, error rate, and R(2), and introduced the sensitivity and specificity to this context. We exemplified several measures in the KORA study, a large population-based study from Southern Germany. We find that the specificity is slightly reduced only for common haplotypes, while the sensitivity was decreased for some, but not all rare haplotypes. The overall error rate was generally increasing with increasing number of loci, increasing minor allele frequency of SNPs, decreasing correlation between the alleles and increasing ambiguity. We conclude that, with the analytical approach presented here, haplotype-specific error measures can be computed to gain insight into the haplotype uncertainty. This method provides the information, if a specific risk haplotype can be expected to be reconstructed with rather no or high misclassification and thus on the magnitude of expected bias in association estimates. We also illustrate that sensitivity and specificity separate two dimensions of the haplotype reconstruction error, which completely describe the misclassification matrix and thus provide the prerequisite for methods accounting for misclassification

SNPTransformer: A Lightweight Toolkit for Genome-Wide Association Studies

AbstractHigh-throughput genotyping chips have produced huge datasets for genome-wide association studies (GWAS) that have contributed greatly to discovering susceptibility genes for complex diseases. There are two strategies for performing data analysis for GWAS. One strategy is to use open-source or commercial packages that are designed for GWAS. The other is to take advantage of classic genetic programs with specific functions, such as linkage disequilibrium mapping, haplotype inference and transmission disequilibrium tests. However, most classic programs that are available are not suitable for analyzing chip data directly and require custom-made input, which results in the inconvenience of converting raw genotyping files into various data formats. We developed a powerful, user-friendly, lightweight program named SNPTransformer for GWAS that includes five major modules (Transformer, Operator, Previewer, Coder and Simulator). The toolkit not only works for transforming the genotyping files into ten input formats for use with classic genetics packages, but also carries out useful functions such as relational operations on IDs, previewing data files, recoding data formats and simulating marker files, among other functions. It bridges upstream raw genotyping data with downstream genetic programs, and can act as an in-hand toolkit for human geneticists, especially for non-programmers. SNPTransformer is freely available at http://snptransformer.sourceforge.net

Interaction between IRF6 and TGFA Genes Contribute to the Risk of Nonsyndromic Cleft Lip/Palate

Previous evidence from tooth agenesis studies suggested IRF6 and TGFA interact. Since tooth agenesis is commonly found in individuals with cleft lip/palate (CL/P), we used four large cohorts to evaluate if IRF6 and TGFA interaction contributes to CL/P. Markers within and flanking IRF6 and TGFA genes were tested using Taqman or SYBR green chemistries for case-control analyses in 1,000 Brazilian individuals. We looked for evidence of gene-gene interaction between IRF6 and TGFA by testing if markers associated with CL/P were overtransmitted together in the case-control Brazilian dataset and in the additional family datasets. Genotypes for an additional 142 case-parent trios from South America drawn from the Latin American Collaborative Study of Congenital Malformations (ECLAMC), 154 cases from Latvia, and 8,717 individuals from several cohorts were available for replication of tests for interaction. Tgfa and Irf6 expression at critical stages during palatogenesis was analyzed in wild type and Irf6 knockout mice. Markers in and near IRF6 and TGFA were associated with CL/P in the Brazilian cohort (p<10-6). IRF6 was also associated with cleft palate (CP) with impaction of permanent teeth (p<10-6). Statistical evidence of interaction between IRF6 and TGFA was found in all data sets (p = 0.013 for Brazilians; p = 0.046 for ECLAMC; p = 10-6 for Latvians, and p = 0.003 for the 8,717 individuals). Tgfa was not expressed in the palatal tissues of Irf6 knockout mice. IRF6 and TGFA contribute to subsets of CL/P with specific dental anomalies. Moreover, this potential IRF6-TGFA interaction may account for as much as 1% to 10% of CL/P cases. The Irf6-knockout model further supports the evidence of IRF6-TGFA interaction found in humans. © 2012 Letra et al

Sequence Variations and Haplotypes of the GJB2 Gene Revealed by Resequencing of 192 Chromosomes from the General Population in Korea

ObjectivesHearing impairment (HI) is the most common sensory deficit in human. The Gap Junction Protein, Beta-2 (GJB2) gene encodes the protein connexin 26, and this gene accounts for up to half of the cases of autosomal recessive nonsyndromic HI. This study was conducted to obtain a set of sequence variations (SVs) of the GJB2 gene among Koreans from the general population for making molecular genetic diagnoses and performing genetic counseling.MethodsWe resequenced the GJB2 gene in 192 chromosomes from 96 adult individuals of Korean descent and who were without a history of hearing difficulty. The data of the SVs was obtained and the haplotypes were reconstructed from the data.ResultsFive SVs were observed, including a novel one (c.558G>A; p.T186T), with the allele frequencies ranging from 0.5% (1/192) to 41% (79/192). The linkage disequilibrium study and haplotype construction showed that some of the SVs are in tight linkage, resulting in a limited number of haplotypes.ConclusionWe observed SVs of the GJB2 gene with different allele frequencies, and a limited number of haplotypes were constructed. The data from this study can be used as reference data for GJB2-related hearing genetic studies, including studies on the founder effect and population genetics, and this data is particularly relevant to people of East Asian decent