Whole-exome sequencing of 14 389 individuals from the ESP and CHARGE consortia identifies novel rare variation associated with hemostatic factors.

Plasma levels of fibrinogen, coagulation factors VII and VIII and von Willebrand factor (vWF) are four intermediate phenotypes that are heritable and have been associated with the risk of clinical thrombotic events. To identify rare and low-frequency variants associated with these hemostatic factors, we conducted whole-exome sequencing in 10 860 individuals of European ancestry (EA) and 3529 African Americans (AAs) from the Cohorts for Heart and Aging Research in Genomic Epidemiology Consortium and the National Heart, Lung and Blood Institute's Exome Sequencing Project. Gene-based tests demonstrated significant associations with rare variation (minor allele frequency < 5%) in fibrinogen gamma chain (FGG) (with fibrinogen, P = 9.1 × 10-13), coagulation factor VII (F7) (with factor VII, P = 1.3 × 10-72; seven novel variants) and VWF (with factor VIII and vWF; P = 3.2 × 10-14; one novel variant). These eight novel rare variant associations were independent of the known common variants at these loci and tended to have much larger effect sizes. In addition, one of the rare novel variants in F7 was significantly associated with an increased risk of venous thromboembolism in AAs (Ile200Ser; rs141219108; P = 4.2 × 10-5). After restricting gene-based analyses to only loss-of-function variants, a novel significant association was detected and replicated between factor VIII levels and a stop-gain mutation exclusive to AAs (rs3211938) in CD36 molecule (CD36). This variant has previously been linked to dyslipidemia but not with the levels of a hemostatic factor. These efforts represent the largest integration of whole-exome sequence data from two national projects to identify genetic variation associated with plasma hemostatic factors.

Epigenomic partitioning of a polygenic risk score for asthma reveals distinct genetically driven disease pathways.

BACKGROUND: Individual differences in susceptibility to developing asthma, a heterogeneous chronic inflammatory lung disease, are poorly understood. Whether genetics can predict asthma risk and how genetic variants modulate the complex pathophysiology of asthma are still debated. AIM: To build polygenic risk scores for asthma risk prediction and epigenomically link predictive genetic variants to pathophysiological mechanisms. METHODS: Restricted polygenic risk scores were constructed using single nucleotide variants derived from genome-wide association studies and validated using data generated in the Rotterdam Study, a Dutch prospective cohort of 14 926 individuals. Outcomes used were asthma, childhood-onset asthma, adulthood-onset asthma, eosinophilic asthma and asthma exacerbations. Genome-wide chromatin analysis data from 19 disease-relevant cell types were used for epigenomic polygenic risk score partitioning. RESULTS: The polygenic risk scores obtained predicted asthma and related outcomes, with the strongest associations observed for childhood-onset asthma (2.55 odds ratios per polygenic risk score standard deviation, area under the curve of 0.760). Polygenic risk scores allowed for the classification of individuals into high-risk and low-risk groups. Polygenic risk score partitioning using epigenomic profiles identified five clusters of variants within putative gene regulatory regions linked to specific asthma-relevant cells, genes and biological pathways. CONCLUSIONS: Polygenic risk scores were associated with asthma(-related traits) in a Dutch prospective cohort, with substantially higher predictive power observed for childhood-onset than adult-onset asthma. Importantly, polygenic risk score variants could be epigenomically partitioned into clusters of regulatory variants with different pathophysiological association patterns and effect estimates, which likely represent distinct genetically driven disease pathways. Our findings have potential implications for personalised risk mitigation and treatment strategies.

Effects and interaction of single nucleotide polymorphisms at the pharmacokinetic/pharmacodynamic site: insights from the Rotterdam study into metformin clinical response and dose titration.

Our study investigated the impact of genetic variations on metformin glycemic response in a cohort from the Rotterdam Study, comprising 14,926 individuals followed for up to 27 years. Among 1285 metformin users of European ancestry, using linear mixed models, we analyzed the association of single nucleotide polymorphisms (SNPs) and a Polygenic Risk Score (PRS) with glycemic response, measured by changes in metformin dosage or HbA1c levels. While individual genetic variants showed no significant association, rs622342 on SLC2A1 correlated with increased glycemic response only in metformin monotherapy patients (ß = -2.09, P-value < 0.001). The collective effect of variants, as represented by PRS, weakly correlated with changes in metformin dosage (ß = 0.023, P-value = 0.027). Synergistic interaction was observed between rs7124355 and rs8192675. Our findings suggest that while higher PRS correlates with increased metformin dosage, its modest effect size limits clinical utility, emphasizing the need for future research in diverse populations to refine genetic risk models.

Genetic counselling legislation and practice in cancer in EU Member States.

BACKGROUND: Somatic and germline genetic alterations are significant drivers of cancer. Increasing integration of new technologies which profile these alterations requires timely, equitable and high-quality genetic counselling to facilitate accurate diagnoses and informed decision-making by patients and their families in preventive and clinical settings. This article aims to provide an overview of genetic counselling legislation and practice across European Union (EU) Member States to serve as a foundation for future European recommendations and action. METHODS: National legislative databases of all 27 Member States were searched using terms relevant to genetic counselling, translated as appropriate. Interviews with relevant experts from each Member State were conducted to validate legislative search results and provide detailed insights into genetic counselling practice in each country. RESULTS: Genetic counselling is included in national legislative documents of 22 of 27 Member States, with substantial variation in legal mechanisms and prescribed details (i.e. the 'who, what, when and where' of counselling). Practice is similarly varied. Workforce capacity (25 of 27 Member States) and genetic literacy (all Member States) were common reported barriers. Recognition and/or better integration of genetic counsellors and updated legislation and were most commonly noted as the 'most important change' which would improve practice. CONCLUSIONS: This review highlights substantial variability in genetic counselling across EU Member States, as well as common barriers notwithstanding this variation. Future recommendations and action should focus on addressing literacy and capacity challenges through legislative, regulatory and/or strategic approaches at EU, national, regional and/or local levels.

Multi-source data approach for personalized outcome prediction in lung cancer screening: update from the NELSON trial.

Trials show that low-dose computed tomography (CT) lung cancer screening in long-term (ex-)smokers reduces lung cancer mortality. However, many individuals were exposed to unnecessary diagnostic procedures. This project aims to improve the efficiency of lung cancer screening by identifying high-risk participants, and improving risk discrimination for nodules. This study is an extension of the Dutch-Belgian Randomized Lung Cancer Screening Trial, with a focus on personalized outcome prediction (NELSON-POP). New data will be added on genetics, air pollution, malignancy risk for lung nodules, and CT biomarkers beyond lung nodules (emphysema, coronary calcification, bone density, vertebral height and body composition). The roles of polygenic risk scores and air pollution in screen-detected lung cancer diagnosis and survival will be established. The association between the AI-based nodule malignancy score and lung cancer will be evaluated at baseline and incident screening rounds. The association of chest CT imaging biomarkers with outcomes will be established. Based on these results, multisource prediction models for pre-screening and post-baseline-screening participant selection and nodule management will be developed. The new models will be externally validated. We hypothesize that we can identify 15-20% participants with low-risk of lung cancer or short life expectancy and thus prevent ~140,000 Dutch individuals from being screened unnecessarily. We hypothesize that our models will improve the specificity of nodule management by 10% without loss of sensitivity as compared to assessment of nodule size/growth alone, and reduce unnecessary work-up by 40-50%.

Donor lamella thickness after ultrathin Descemet stripping automated endothelial keratoplasty and its relation to postoperative visual acuity and pre-operative lamella measures.

BACKGROUND: To accurately measure ultrathin Descemet stripping automated endothelial keratoplasty (DSAEK) donor lamella thickness during the first postoperative year and to correlate this with pre-operative and other postoperative measurements. METHODS: Donor lamella thickness in 41 eyes undergoing DSAEK for Fuchs endothelial dystrophy (FED) was measured using the Tomey Casia OCT directly after graft preparation and at 1 week and 1, 3, 6 and 12 months postoperatively. Visual acuity and endothelial cell density were measured as the secondary parameters. RESULTS: Individual graft thickness profiles were shown to be fairly regular within the optically relevant area. There was a strong and highly significant correlation between the pre- and postoperative lamellar thicknesses at all time points (p < 0.0001). Compared with the measurements directly after preparation at the cornea bank, the lamella thickness decreased by 12% after 12 months. Between 1 and 12 months postoperatively, the lamella thickness (mean ± SD) changed from 112 ± 27 µm to 101 ± 21 µm. Best spectacle-corrected visual acuity (BSCVA) changed from 0.46 ± 0.30 logMAR pre-operatively through 0.36 ± 0.33 at 1 month to 0.13 ± 0.16 at 1 year postoperatively. The endothelial cell counts were comparable to those reported in previous studies. CONCLUSIONS: Thickness profiles of individual grafts were fairly regular within the optically relevant area. A strong relationship between pre- and postoperative graft thicknesses was detected, and ultrathin DSAEK grafts prepared using methods similar to that applied in this study are expected to show a deswelling of around 12% during the first postoperative year. No correlation was detected between graft thickness and BSCVA.

Coronary Artery Calcium Score and Polygenic Risk Score for the Prediction of Coronary Heart Disease Events.

Importance: Coronary artery calcium score and polygenic risk score have each separately been proposed as novel markers to identify risk of coronary heart disease (CHD), but no prior studies have directly compared these markers in the same cohorts. Objective: To evaluate change in CHD risk prediction when a coronary artery calcium score, a polygenic risk score, or both are added to a traditional risk factor-based model. Design, Setting, and Participants: Two observational population-based studies involving individuals aged 45 years through 79 years of European ancestry and free of clinical CHD at baseline: the Multi-Ethnic Study of Atherosclerosis (MESA) study involved 1991 participants at 6 US centers and the Rotterdam Study (RS) involved 1217 in Rotterdam, the Netherlands. Exposure: Traditional risk factors were used to calculate CHD risk (eg, pooled cohort equations [PCEs]), computed tomography for the coronary artery calcium score, and genotyped samples for a validated polygenic risk score. Main Outcomes and Measures: Model discrimination, calibration, and net reclassification improvement (at the recommended risk threshold of 7.5%) for prediction of incident CHD events were assessed. Results: The median age was 61 years in MESA and 67 years in RS. Both log (coronary artery calcium+1) and polygenic risk score were significantly associated with 10-year risk of incident CHD (hazards ratio per SD, 2.60; 95% CI, 2.08-3.26 and 1.43; 95% CI, 1.20-1.71, respectively), in MESA. The C statistic for the coronary artery calcium score was 0.76 (95% CI, 0.71-0.79) and for the polygenic risk score, 0.69 (95% CI, 0.63-0.71). The change in the C statistic when each was added to the PCEs was 0.09 (95% CI, 0.06-0.13) for the coronary artery calcium score, 0.02 (95% CI, 0.00-0.04) for the polygenic risk score, and 0.10 (95% CI, 0.07-0.14) for both. Overall categorical net reclassification improvement was significant when the coronary artery calcium score (0.19; 95% CI, 0.06-0.28) but was not significant when the polygenic risk score (0.04; 95% CI, -0.05 to 0.10) was added to the PCEs. Calibration of the PCEs and models with coronary artery calcium and/or polygenic risk scores was adequate (all χ2<20). Subgroup analysis stratified by the median age demonstrated similar findings. Similar findings were observed for 10-year risk in RS and in longer-term follow-up in MESA (median, 16.0 years). Conclusions and Relevance: In 2 cohorts of middle-aged to older adults from the US and the Netherlands, the coronary artery calcium score had better discrimination than the polygenic risk score for risk prediction of CHD. In addition, the coronary artery calcium score but not the polygenic risk score significantly improved risk discrimination and risk reclassification for CHD when added to traditional risk factors.

Identification of Side Chain Oxidized Sterols as Novel Liver X Receptor Agonists with Therapeutic Potential in the Treatment of Cardiovascular and Neurodegenerative Diseases.

The nuclear receptors-liver X receptors (LXR α and ß) are potential therapeutic targets in cardiovascular and neurodegenerative diseases because of their key role in the regulation of lipid homeostasis and inflammatory processes. Specific oxy(phyto)sterols differentially modulate the transcriptional activity of LXRs providing opportunities to develop compounds with improved therapeutic characteristics. We isolated oxyphytosterols from Sargassum fusiforme and synthesized sidechain oxidized sterol derivatives. Five 24-oxidized sterols demonstrated a high potency for LXRα/ß activation in luciferase reporter assays and induction of LXR-target genes APOE, ABCA1 and ABCG1 involved in cellular cholesterol turnover in cultured cells: methyl 3ß-hydroxychol-5-en-24-oate (S1), methyl (3ß)-3-aldehydeoxychol-5-en-24-oate (S2), 24-ketocholesterol (S6), (3ß,22E)-3-hydroxycholesta-5,22-dien-24-one (N10) and fucosterol-24,28 epoxide (N12). These compounds induced SREBF1 but not SREBP1c-mediated lipogenic genes such as SCD1, ACACA and FASN in HepG2 cells or astrocytoma cells. Moreover, S2 and S6 enhanced cholesterol efflux from HepG2 cells. All five oxysterols induced production of the endogenous LXR agonists 24(S)-hydroxycholesterol by upregulating the CYP46A1, encoding the enzyme converting cholesterol into 24(S)-hydroxycholesterol; S1 and S6 may also act via the upregulation of desmosterol production. Thus, we identified five novel LXR-activating 24-oxidized sterols with a potential for therapeutic applications in neurodegenerative and cardiovascular diseases.

Somatic TARDBP variants as a cause of semantic dementia.

The aetiology of late-onset neurodegenerative diseases is largely unknown. Here we investigated whether de novo somatic variants for semantic dementia can be detected, thereby arguing for a more general role of somatic variants in neurodegenerative disease. Semantic dementia is characterized by a non-familial occurrence, early onset (<65 years), focal temporal atrophy and TDP-43 pathology. To test whether somatic variants in neural progenitor cells during brain development might lead to semantic dementia, we compared deep exome sequencing data of DNA derived from brain and blood of 16 semantic dementia cases. Somatic variants observed in brain tissue and absent in blood were validated using amplicon sequencing and digital PCR. We identified two variants in exon one of the TARDBP gene (L41F and R42H) at low level (1-3%) in cortical regions and in dentate gyrus in two semantic dementia brains, respectively. The pathogenicity of both variants is supported by demonstrating impaired splicing regulation of TDP-43 and by altered subcellular localization of the mutant TDP-43 protein. These findings indicate that somatic variants may cause semantic dementia as a non-hereditary neurodegenerative disease, which might be exemplary for other late-onset neurodegenerative disorders.

The Genetics of Atypical Femur Fractures-a Systematic Review.

PURPOSE OF REVIEW: Atypical femur fractures (AFFs) are rare subtrochanteric or diaphyseal fractures regarded as side effects of bisphosphonates (BPs), possibly with a genetic background. Here, we summarize the most recent knowledge about genetics of AFFs. RECENT FINDINGS: AFF has been reported in 57 patients with seven different monogenic bone disorders including hypophosphatasia and osteogenesis imperfecta; 56.1% had never used BPs, while 17.5% were diagnosed with the disorder only after the AFF. Gene mutation finding in familial and sporadic cases identified possible AFF-related variants in the GGPS1 and ATRAID genes respectively. Functional follow-up studies of mutant proteins showed possible roles in AFF. A recent small genome-wide association study on 51 AFF cases did not identify significant hits associated with AFF. Recent findings have strengthened the hypothesis that AFFs have underlying genetic components but more studies are needed in AFF families and larger cohorts of sporadic cases to confirm previous results and/or find novel gene variants involved in the pathogenesis of AFFs.

Molecular Pathways Involved in Frontotemporal Lobar Degeneration with TDP-43 Proteinopathy: What Can We Learn from Proteomics?

Frontotemporal lobar degeneration (FTLD) is a neurodegenerative disorder clinically characterized by behavioral, language, and motor symptoms, with major impact on the lives of patients and their families. TDP-43 proteinopathy is the underlying neuropathological substrate in the majority of cases, referred to as FTLD-TDP. Several genetic causes have been identified, which have revealed some components of its pathophysiology. However, the exact mechanisms driving FTLD-TDP remain largely unknown, forestalling the development of therapies. Proteomic approaches, in particular high-throughput mass spectrometry, hold promise to help elucidate the pathogenic molecular and cellular alterations. In this review, we describe the main findings of the proteomic profiling studies performed on human FTLD-TDP brain tissue. Subsequently, we address the major biological pathways implicated in FTLD-TDP, by reviewing these data together with knowledge derived from genomic and transcriptomic literature. We illustrate that an integrated perspective, encompassing both proteomic, genetic, and transcriptomic discoveries, is vital to unravel core disease processes, and to enable the identification of disease biomarkers and therapeutic targets for this devastating disorder.

Reduced penetrance of pathogenic ACMG variants in a deeply phenotyped cohort study and evaluation of ClinVar classification over time.

PURPOSE: We studied the penetrance of pathogenically classified variants in an elderly Dutch population from the Rotterdam Study, for which deep phenotyping is available. We screened the 59 actionable genes for which reporting of known pathogenic variants was recommended by the American College of Medical Genetics and Genomics (ACMG), and demonstrate that determining what constitutes a known pathogenic variant can be quite challenging. METHODS: We defined "known pathogenic" as classified pathogenic by both ClinVar and the Human Gene Mutation Database (HGMD). In 2628 individuals, we performed exome sequencing and identified known pathogenic variants. We investigated the clinical records of carriers and evaluated clinical events during 25 years of follow-up for evidence of variant pathogenicity. RESULTS: Of 3815 variants detected in the 59 ACMG genes, 17 variants were considered known pathogenic. For 14/17 variants the ClinVar classification had changed over time. Of 24 confirmed carriers of these variants, we observed at least one clinical event possibly caused by the variant in only three participants (13%). CONCLUSION: We show that the definition of "known pathogenic" is often unclear and should be approached carefully. Additionally variants marked as known pathogenic do not always have clinical impact on their carriers. Definition and classification of true (individual) expected pathogenic impact should be defined carefully.

Validation of the BOADICEA model and a 313-variant polygenic risk score for breast cancer risk prediction in a Dutch prospective cohort.

PURPOSE: We evaluated the performance of the recently extended Breast and Ovarian Analysis of Disease Incidence and Carrier Estimation Algorithm (BOADICEA version 5) in a Dutch prospective cohort, using a polygenic risk score (PRS) based on 313 breast cancer (BC)-associated variants (PRS313) and other, nongenetic risk factors. METHODS: Since 1989, 6522 women without BC aged 45 or older of European descent have been included in the Rotterdam Study. The PRS313 was calculated per 1 SD in controls from the Breast Cancer Association Consortium (BCAC). Cox regression analysis was performed to estimate the association between the PRS313 and incident BC risk. Cumulative 10-year risks were calculated with BOADICEA including different sets of variables (age, risk factors and PRS313). C-statistics were used to evaluate discriminative ability. RESULTS: In total, 320 women developed BC. The PRS313 was significantly associated with BC (hazard ratio [HR] per SD of 1.56, 95% confidence interval [CI] [1.40-1.73]). Using 10-year risk estimates including age and the PRS313, other risk factors improved the discriminatory ability of the BOADICEA model marginally, from a C-statistic of 0.636 to 0.653. CONCLUSIONS: The effect size of the PRS313 is highly reproducible in the Dutch population. Our results validate the BOADICEA v5 model for BC risk assessment in the Dutch general population.

Genome-wide analyses as part of the international FTLD-TDP whole-genome sequencing consortium reveals novel disease risk factors and increases support for immune dysfunction in FTLD.

Frontotemporal lobar degeneration with neuronal inclusions of the TAR DNA-binding protein 43 (FTLD-TDP) represents the most common pathological subtype of FTLD. We established the international FTLD-TDP whole-genome sequencing consortium to thoroughly characterize the known genetic causes of FTLD-TDP and identify novel genetic risk factors. Through the study of 1131 unrelated Caucasian patients, we estimated that C9orf72 repeat expansions and GRN loss-of-function mutations account for 25.5% and 13.9% of FTLD-TDP patients, respectively. Mutations in TBK1 (1.5%) and other known FTLD genes (1.4%) were rare, and the disease in 57.7% of FTLD-TDP patients was unexplained by the known FTLD genes. To unravel the contribution of common genetic factors to the FTLD-TDP etiology in these patients, we conducted a two-stage association study comprising the analysis of whole-genome sequencing data from 517 FTLD-TDP patients and 838 controls, followed by targeted genotyping of the most associated genomic loci in 119 additional FTLD-TDP patients and 1653 controls. We identified three genome-wide significant FTLD-TDP risk loci: one new locus at chromosome 7q36 within the DPP6 gene led by rs118113626 (p value = 4.82e - 08, OR = 2.12), and two known loci: UNC13A, led by rs1297319 (p value = 1.27e - 08, OR = 1.50) and HLA-DQA2 led by rs17219281 (p value = 3.22e - 08, OR = 1.98). While HLA represents a locus previously implicated in clinical FTLD and related neurodegenerative disorders, the association signal in our study is independent from previously reported associations. Through inspection of our whole-genome sequence data for genes with an excess of rare loss-of-function variants in FTLD-TDP patients (n ≥ 3) as compared to controls (n = 0), we further discovered a possible role for genes functioning within the TBK1-related immune pathway (e.g., DHX58, TRIM21, IRF7) in the genetic etiology of FTLD-TDP. Together, our study based on the largest cohort of unrelated FTLD-TDP patients assembled to date provides a comprehensive view of the genetic landscape of FTLD-TDP, nominates novel FTLD-TDP risk loci, and strongly implicates the immune pathway in FTLD-TDP pathogenesis.

A C6orf10/LOC101929163 locus is associated with age of onset in C9orf72 carriers.

The G4C2-repeat expansion in C9orf72 is the most common known cause of amyotrophic lateral sclerosis and frontotemporal dementia. The high phenotypic heterogeneity of C9orf72 patients includes a wide range in age of onset, modifiers of which are largely unknown. Age of onset could be influenced by environmental and genetic factors both of which may trigger DNA methylation changes at CpG sites. We tested the hypothesis that age of onset in C9orf72 patients is associated with some common single nucleotide polymorphisms causing a gain or loss of CpG sites and thus resulting in DNA methylation alterations. Combined analyses of epigenetic and genetic data have the advantage of detecting functional variants with reduced likelihood of false negative results due to excessive correction for multiple testing in genome-wide association studies. First, we estimated the association between age of onset in C9orf72 patients (n = 46) and the DNA methylation levels at all 7603 CpG sites available on the 450 k BeadChip that are mapped to common single nucleotide polymorphisms. This was followed by a genetic association study of the discovery (n = 144) and replication (n = 187) C9orf72 cohorts. We found that age of onset was reproducibly associated with polymorphisms within a 124.7 kb linkage disequilibrium block tagged by top-significant variation, rs9357140, and containing two overlapping genes (LOC101929163 and C6orf10). A meta-analysis of all 331 C9orf72 carriers revealed that every A-allele of rs9357140 reduced hazard by 30% (P = 0.0002); and the median age of onset in AA-carriers was 6 years later than GG-carriers. In addition, we investigated a cohort of C9orf72 negative patients (n = 2634) affected by frontotemporal dementia and/or amyotrophic lateral sclerosis; and also found that the AA-genotype of rs9357140 was associated with a later age of onset (adjusted P = 0.007 for recessive model). Phenotype analyses detected significant association only in the largest subgroup of patients with frontotemporal dementia (n = 2142, adjusted P = 0.01 for recessive model). Gene expression studies of frontal cortex tissues from 25 autopsy cases affected by amyotrophic lateral sclerosis revealed that the G-allele of rs9357140 is associated with increased brain expression of LOC101929163 (a non-coding RNA) and HLA-DRB1 (involved in initiating immune responses), while the A-allele is associated with their reduced expression. Our findings suggest that carriers of the rs9357140 GG-genotype (linked to an earlier age of onset) might be more prone to be in a pro-inflammatory state (e.g. by microglia) than AA-carriers. Further, investigating the functional links within the C6orf10/LOC101929163/HLA-DRB1 pathway will be critical to better define age-dependent pathogenesis of frontotemporal dementia and amyotrophic lateral sclerosis.

Novel Genetic Variants for Cartilage Thickness and Hip Osteoarthritis.

Osteoarthritis is one of the most frequent and disabling diseases of the elderly. Only few genetic variants have been identified for osteoarthritis, which is partly due to large phenotype heterogeneity. To reduce heterogeneity, we here examined cartilage thickness, one of the structural components of joint health. We conducted a genome-wide association study of minimal joint space width (mJSW), a proxy for cartilage thickness, in a discovery set of 13,013 participants from five different cohorts and replication in 8,227 individuals from seven independent cohorts. We identified five genome-wide significant (GWS, P≤5·0×10-8) SNPs annotated to four distinct loci. In addition, we found two additional loci that were significantly replicated, but results of combined meta-analysis fell just below the genome wide significance threshold. The four novel associated genetic loci were located in/near TGFA (rs2862851), PIK3R1 (rs10471753), SLBP/FGFR3 (rs2236995), and TREH/DDX6 (rs496547), while the other two (DOT1L and SUPT3H/RUNX2) were previously identified. A systematic prioritization for underlying causal genes was performed using diverse lines of evidence. Exome sequencing data (n = 2,050 individuals) indicated that there were no rare exonic variants that could explain the identified associations. In addition, TGFA, FGFR3 and PIK3R1 were differentially expressed in OA cartilage lesions versus non-lesioned cartilage in the same individuals. In conclusion, we identified four novel loci (TGFA, PIK3R1, FGFR3 and TREH) and confirmed two loci known to be associated with cartilage thickness.The identified associations were not caused by rare exonic variants. This is the first report linking TGFA to human OA, which may serve as a new target for future therapies.

A systematic analysis highlights multiple long non-coding RNAs associated with cardiometabolic disorders.

Genome-wide association studies (GWAS) have identified many susceptibility loci for cardiometabolic disorders. Most of the associated variants reside in non-coding regions of the genome including long non-coding RNAs (lncRNAs), which are thought to play critical roles in diverse biological processes. Here, we leveraged data from the available GWAS meta-analyses on lipid and obesity-related traits, blood pressure, type 2 diabetes, and coronary artery disease and identified 179 associated single-nucleotide polymorphisms (SNPs) in 102 lncRNAs (p-value < 2.3 × 10-7). Of these, 55 SNPs, either the lead SNP or in strong linkage disequilibrium with the lead SNP in the related loci, were selected for further investigations. Our in silico predictions and functional annotations of the SNPs as well as expression and DNA methylation analysis of their lncRNAs demonstrated several lncRNAs that fulfilled predefined criteria for being potential functional targets. In particular, we found evidence suggesting that LOC157273 (at 8p23.1) is involved in regulating serum lipid-cholesterol. Our results showed that rs4841132 in the second exon and cg17371580 in the promoter region of LOC157273 are associated with lipids; the lncRNA is expressed in liver and associates with the expression of its nearby coding gene, PPP1R3B. Collectively, we highlight a number of loci associated with cardiometabolic disorders for which the association may act through lncRNAs.

Whole-Exome Sequencing in Age-Related Macular Degeneration Identifies Rare Variants in COL8A1, a Component of Bruch's Membrane.

PURPOSE: Genome-wide association studies and targeted sequencing studies of candidate genes have identified common and rare variants that are associated with age-related macular degeneration (AMD). Whole-exome sequencing (WES) studies allow a more comprehensive analysis of rare coding variants across all genes of the genome and will contribute to a better understanding of the underlying disease mechanisms. To date, the number of WES studies in AMD case-control cohorts remains scarce and sample sizes are limited. To scrutinize the role of rare protein-altering variants in AMD cause, we performed the largest WES study in AMD to date in a large European cohort consisting of 1125 AMD patients and 1361 control participants. DESIGN: Genome-wide case-control association study of WES data. PARTICIPANTS: One thousand one hundred twenty-five AMD patients and 1361 control participants. METHODS: A single variant association test of WES data was performed to detect variants that are associated individually with AMD. The cumulative effect of multiple rare variants with 1 gene was analyzed using a gene-based CMC burden test. Immunohistochemistry was performed to determine the localization of the Col8a1 protein in mouse eyes. MAIN OUTCOME MEASURES: Genetic variants associated with AMD. RESULTS: We detected significantly more rare protein-altering variants in the COL8A1 gene in patients (22/2250 alleles [1.0%]) than in control participants (11/2722 alleles [0.4%]; P = 7.07×10-5). The association of rare variants in the COL8A1 gene is independent of the common intergenic variant (rs140647181) near the COL8A1 gene previously associated with AMD. We demonstrated that the Col8a1 protein localizes at Bruch's membrane. CONCLUSIONS: This study supported a role for protein-altering variants in the COL8A1 gene in AMD pathogenesis. We demonstrated the presence of Col8a1 in Bruch's membrane, further supporting the role of COL8A1 variants in AMD pathogenesis. Protein-altering variants in COL8A1 may alter the integrity of Bruch's membrane, contributing to the accumulation of drusen and the development of AMD.

Excessive burden of lysosomal storage disorder gene variants in Parkinson's disease.

Mutations in the glucocerebrosidase gene (GBA), which cause Gaucher disease, are also potent risk factors for Parkinson's disease. We examined whether a genetic burden of variants in other lysosomal storage disorder genes is more broadly associated with Parkinson's disease susceptibility. The sequence kernel association test was used to interrogate variant burden among 54 lysosomal storage disorder genes, leveraging whole exome sequencing data from 1156 Parkinson's disease cases and 1679 control subjects. We discovered a significant burden of rare, likely damaging lysosomal storage disorder gene variants in association with Parkinson's disease risk. The association signal was robust to the exclusion of GBA, and consistent results were obtained in two independent replication cohorts, including 436 cases and 169 controls with whole exome sequencing and an additional 6713 cases and 5964 controls with exome-wide genotyping. In secondary analyses designed to highlight the specific genes driving the aggregate signal, we confirmed associations at the GBA and SMPD1 loci and newly implicate CTSD, SLC17A5, and ASAH1 as candidate Parkinson's disease susceptibility genes. In our discovery cohort, the majority of Parkinson's disease cases (56%) have at least one putative damaging variant in a lysosomal storage disorder gene, and 21% carry multiple alleles. Our results highlight several promising new susceptibility loci and reinforce the importance of lysosomal mechanisms in Parkinson's disease pathogenesis. We suggest that multiple genetic hits may act in combination to degrade lysosomal function, enhancing Parkinson's disease susceptibility.

A Randomized Multicenter Clinical Trial of Ultrathin Descemet Stripping Automated Endothelial Keratoplasty (DSAEK) versus DSAEK.

OBJECTIVE: To compare visual acuity, refraction, endothelial cell density (ECD), and complications after Descemet stripping automated endothelial keratoplasty (DSAEK) and ultrathin DSAEK (UT-DSAEK). DESIGN: A multicenter, prospective, double-masked, randomized, controlled clinical trial. PARTICIPANTS: From 66 patients with irreversible corneal endothelial dysfunction dues to Fuchs' dystrophy who enrolled from 4 tertiary medical centers in the Netherlands, 66 eyes were studied. METHODS: Participants were centrally randomized to undergo either UT-DSAEK or DSAEK, based on preoperative best spectacle-corrected visual acuity (BSCVA), recipient central corneal thickness, patient age, and recruitment center. Donor corneas were precut by a single cornea bank. PARTICIPANTS: Participants underwent ophthalmic examinations preoperatively and 3, 6, and 12 months after the operation, including manifest refraction, BSCVA using an Early Treatment Diabetic Retinopathy Study chart, and endothelium imaging. MAIN OUTCOME MEASURES: BSCVA 12 months postoperatively. RESULTS: Preoperative BSCVA did not differ between patients undergoing DSAEK (0.35 logarithm of the minimum angle of resolution [logMAR] [95% confidence interval {CI} 0.27-0.43]; n = 32) and UT-DSAEK (0.37 logMAR [95% CI 0.31-0.43]; n = 34; P = 0.8). BSCVA was significantly better after UT-DSAEK compared with that after DSAEK at 3 months (0.17 logMAR [95% CI 0.13-0.21], n = 31 vs. 0.28 logMAR [95% CI 0.23-0.33], n = 31; P = 0.001), 6 months (0.14 logMAR [95% CI 0.10-0.18], n = 30 vs. 0.24 logMAR [95% CI 0.20-0.28], n = 30; P = 0.002), and 12 months (0.13 logMAR [95% CI 0.09-0.17], n = 33 vs. 0.20 logMAR [95% CI 0.15-0.25], n = 29; P = 0.03). Refraction, ECD loss (40% at 3 months; P < 0.001), donor loss (DSAEK n = 2 vs. UT-DSAEK n = 3 [relative risk {RR} 1.4 {95% CI 0.24-7.5}; P = 0.7]), and graft dislocation (DSAEK n = 5 vs. UT-DSAEK n = 5 [RR 1.0 {95% CI 0.34-3.33}; P = 0.9]) did not differ between UT-DSAEK and DSAEK. Donor thickness was significantly thinner for UT-DSAEK (101 µm [95% CI 93-110 µm]; range 50-145 µm) than for DSAEK (209 µm [95% CI 196-222 µm]; range 147-289 µm; P < 0.001). CONCLUSIONS: This study indicates that compared with DSAEK, UT-DSAEK results in faster and better recovery of BSCVA with similar refractive outcomes, endothelial cell loss, and incidence of complications.