Further improvements to linear mixed models for genome-wide association studies

We examine improvements to the linear mixed model (LMM) that better correct for population structure and family relatedness in genome-wide association studies (GWAS). LMMs rely on the estimation of a genetic similarity matrix (GSM), which encodes the pairwise similarity between every two individuals in a cohort. These similarities are estimated from single nucleotide polymorphisms (SNPs) or other genetic variants. Traditionally, all available SNPs are used to estimate the GSM. In empirical studies across a wide range of synthetic and real data, we find that modifications to this approach improve GWAS performance as measured by type I error control and power. Specifically, when only population structure is present, a GSM constructed from SNPs that well predict the phenotype in combination with principal components as covariates controls type I error and yields more power than the traditional LMM. In any setting, with or without population structure or family relatedness, a GSM consisting of a mixture of two component GSMs, one constructed from all SNPs and another constructed from SNPs that well predict the phenotype again controls type I error and yields more power than the traditional LMM. Software implementing these improvements and the experimental comparisons are available at http://microsoft.com/science

Comparative Treatment Outcomes for Patients With Idiopathic Subglottic Stenosis.

To access publisher's full text version of this article, please click on the hyperlink in Additional Links field or click on the hyperlink at the top of the page marked DownloadImportance: Surgical treatment comparisons in rare diseases are difficult secondary to the geographic distribution of patients. Fortunately, emerging technologies offer promise to reduce these barriers for research. Objective: To prospectively compare the outcomes of the 3 most common surgical approaches for idiopathic subglottic stenosis (iSGS), a rare airway disease. Design, setting, and participants: In this international, prospective, 3-year multicenter cohort study, 810 patients with untreated, newly diagnosed, or previously treated iSGS were enrolled after undergoing a surgical procedure (endoscopic dilation [ED], endoscopic resection with adjuvant medical therapy [ERMT], or cricotracheal resection [CTR]). Patients were recruited from clinician practices in the North American Airway Collaborative and an online iSGS community on Facebook. Main outcomes and measures: The primary end point was days from initial surgical procedure to recurrent surgical procedure. Secondary end points included quality of life using the Clinical COPD (chronic obstructive pulmonary disease) Questionnaire (CCQ), Voice Handicap Index-10 (VHI-10), Eating Assessment Test-10 (EAT-10), the 12-Item Short-Form Version 2 (SF-12v2), and postoperative complications. Results: Of 810 patients in this cohort, 798 (98.5%) were female and 787 (97.2%) were white, with a median age of 50 years (interquartile range, 43-58 years). Index surgical procedures were ED (n = 603; 74.4%), ERMT (n = 121; 14.9%), and CTR (n = 86; 10.6%). Overall, 185 patients (22.8%) had a recurrent surgical procedure during the 3-year study, but recurrence differed by modality (CTR, 1 patient [1.2%]; ERMT, 15 [12.4%]; and ED, 169 [28.0%]). Weighted, propensity score-matched, Cox proportional hazards regression models showed ED was inferior to ERMT (hazard ratio [HR], 3.16; 95% CI, 1.8-5.5). Among successfully treated patients without recurrence, those treated with CTR had the best CCQ (0.75 points) and SF-12v2 (54 points) scores and worst VHI-10 score (13 points) 360 days after enrollment as well as the greatest perioperative risk. Conclusions and relevance: In this cohort study of 810 patients with iSGS, endoscopic dilation, the most popular surgical approach for iSGS, was associated with a higher recurrence rate compared with other procedures. Cricotracheal resection offered the most durable results but showed the greatest perioperative risk and the worst long-term voice outcomes. Endoscopic resection with medical therapy was associated with better disease control compared with ED and had minimal association with vocal function. These results may be used to inform individual patient treatment decision-making.Patient-Centered Outcomes Research Institute - PCOR

The detrital zircon U-Pb-Hf fingerprint of the northern Arabian-Nubian Shield as reflected by a Late Ediacaran arkosic wedge (Zenifim Formation; subsurface Israel)

The Neoproterozoic Arabian-Nubian Shield (ANS), which comprises the northern sector of the East African Orogen (EAO), has been deeply denuded in Ediacaran times and its erosional products were deposited across NE Africa and Arabia; they are also found i

RL-SKAT: An exact and efficient score Test for heritability and set tests.

Testing for the existence of variance components in linear mixed models is a fundamental task in many applicative fields. In statistical genetics, the score test has recently become instrumental in the task of testing an association between a set of genetic markers and a phenotype. With few markers, this amounts to set-based variance component tests, which attempt to increase power in association studies by aggregating weak individual effects. When the entire genome is considered, it allows testing for the heritability of a phenotype, defined as the proportion of phenotypic variance explained by genetics. In the popular score-based Sequence Kernel Association Test (SKAT) method, the assumed distribution of the score test statistic is uncalibrated in small samples, with a correction being computationally expensive. This may cause severe inflation or deflation of p-values, even when the null hypothesis is true. Here, we characterize the conditions under which this discrepancy holds, and show it may occur also in large real datasets, such as a dataset from the Wellcome Trust Case Control Consortium 2 (n=13,950) study, and in particular when the individuals in the sample are unrelated. In these cases the SKAT approximation tends to be highly over-conservative and therefore underpowered. To address this limitation, we suggest an efficient method to calculate exact p-values for the score test in the case of a single variance component and a continuous response vector, which can speed up the analysis by orders of magnitude. Our results enable fast and accurate application of the score test in heritability and in set-based association tests. Our method is available in http://github.com/cozygene/RL-SKAT

Functional disease architectures reveal unique biological role of transposable elements

Transposable elements (TE) comprise roughly half of the human genome. Though initially derided as junk DNA, they have been widely hypothesized to contribute to the evolution of gene regulation. However, the contribution of TE to the genetic architecture of diseases remains unknown. Here, we analyze data from 41 independent diseases and complex traits to draw three conclusions. First, TE are uniquely informative for disease heritability. Despite overall depletion for heritability (54% of SNPs, 39 ± 2% of heritability), TE explain substantially more heritability than expected based on their depletion for known functional annotations. This implies that TE acquire function in ways that differ from known functional annotations. Second, older TE contribute more to disease heritability, consistent with acquiring biological function. Third, Short Interspersed Nuclear Elements (SINE) are far more enriched for blood traits than for other traits. Our results can help elucidate the biological roles that TE play in the genetic architecture of diseases

Titin mutation in familial restrictive cardiomyopathy

BACKGROUND: Familial restrictive cardiomyopathy (RCM) caused by a single gene mutation is the least common of the inherited cardiomyopathies. Only a few RCM-causing mutations have been described. Most mutations causing RCM are located in sarcomere protein genes which also cause hypertrophic cardiomyopathy (HCM). Other genes associated with RCM include the desmin and familial amyloidosis genes. In the present study we describe familial RCM with severe heart failure triggered by a de novo mutation in TTN, encoding the huge muscle filament protein titin. METHODS AND RESULTS: Family members underwent physical examination, ECG and Doppler echocardiogram studies. The family comprised 6 affected individuals aged 12-35years. Linkage to candidate loci was performed, followed by gene sequencing. Candidate loci/gene analysis excluded 18 candidate genes but showed segregation with a common haplotype surrounding the TTN locus. Sequence analysis identified a de novo mutation within exon 266 of the TTN gene, resulting in the replacement of tyrosine by cysteine. p.Y7621C affects a highly conserved region in the protein within a fibronectin-3 domain, belonging to the A/I junction region of titin. No other disease-causing mutation was identified in cardiomyopathy genes by whole exome sequencing. CONCLUSIONS: Our study shows, for the first time, that mutations in TTN can cause restrictive cardiomyopathy. The giant filament titin is considered to be a determinant of a resting tension of the sarcomere and this report provides genetic evidence of its crucial role in diastolic function