Genetic analysis and natural history of Parkinson's disease due to the LRRK2 G2019S variant.

The LRRK2 G2019S variant is the most common cause of monogenic Parkinson's disease (PD); however, questions remain regarding the penetrance, clinical phenotype and natural history of carriers. We performed a 3.5-year prospective longitudinal online study in a large number of 1286 genotyped LRRK2 G2019S carriers and 109 154 controls, with and without PD, recruited from the 23andMe Research Cohort. We collected self-reported motor and non-motor symptoms every 6 months, as well as demographics, family histories and environmental risk factors. Incident cases of PD (phenoconverters) were identified at follow-up. We determined lifetime risk of PD using accelerated failure time modelling and explored the impact of polygenic risk on penetrance. We also computed the genetic ancestry of all LRRK2 G2019S carriers in the 23andMe database and identified regions of the world where carrier frequencies are highest. We observed that despite a 1 year longer disease duration (P = 0.016), LRRK2 G2019S carriers with PD had similar burden of motor symptoms, yet significantly fewer non-motor symptoms including cognitive difficulties, REM sleep behaviour disorder (RBD) and hyposmia (all P-values ≤ 0.0002). The cumulative incidence of PD in G2019S carriers by age 80 was 49%. G2019S carriers had a 10-fold risk of developing PD versus non-carriers. This rose to a 27-fold risk in G2019S carriers with a PD polygenic risk score in the top 25% versus non-carriers in the bottom 25%. In addition to identifying ancient founding events in people of North African and Ashkenazi descent, our genetic ancestry analyses infer that the G2019S variant was later introduced to Spanish colonial territories in the Americas. Our results suggest LRRK2 G2019S PD appears to be a slowly progressive predominantly motor subtype of PD with a lower prevalence of hyposmia, RBD and cognitive impairment. This suggests that the current prodromal criteria, which are based on idiopathic PD, may lack sensitivity to detect the early phases of LRRK2 PD in G2019S carriers. We show that polygenic burden may contribute to the development of PD in the LRRK2 G2019S carrier population. Collectively, the results should help support screening programmes and candidate enrichment strategies for upcoming trials of LRRK2 inhibitors in early-stage disease.

Exploring the relationship between condition severity and health-related quality of life in people with haemophilia A across Europe: a multivariable analysis of data from the CHESS II study.

BACKGROUND: Haemophilia A (HA; Factor VIII deficiency) is a congenital X-linked bleeding disorder characterized by trauma-related or spontaneous bleeding events, most notably arising within the intraarticular space and resulting in chronic inflammation and degeneration of affected joints. Endogenous clotting factor activity relative to normal levels determines the severity of HA symptoms, as mild (> 5-40%), moderate (1-5%), or severe (< 1%). Within the current environment of rapid evolution in HA management, we seek to understand the interplay of condition severity and health-related quality of life (HRQoL) to characterise and differentiate unmet needs among people with HA (PwHA). METHODS: A generalised linear regression model (GLM) was developed to explore the relationship between HA severity and EQ-5D-5 L index score from adult HA patients sampled in the "Cost of Haemophilia across Europe - a Socioeconomic Survey II" (CHESS II) cross-sectional, retrospective burden of illness study among adults with hereditary haemophilia A or B from eight European countries. HA patients of any severity with no active inhibitors during the 12 months prior to data capture and a completeEQ-5D-5 L response were included. A base GLM model was specified with covariates for demographic and clinical characteristics (age, body mass index, country, employment, HA severity, annual bleeding rate, problem joints, and chronic pain). RESULTS: Of 381 evaluable patients, 221 (58.0%) had severe HA, 96 (25.2%) had moderate HA, and 64 (16.8%) had mild HA. Among the covariates included in the GLM model and after controlling for haemophilia-related outcomes, a significant association was observed between mild HA and higher EQ-5D-5 L index score (average marginal effects, 0.084; p = 0.016) relative to severe HA. Patient country of residence and magnitude of HA-related chronic pain were also associated with significant differences in index scores, with the latter showing a negative relationship with HRQoL outcomes. CONCLUSIONS: Condition severity and chronic pain are significant predictors of HRQoL in PwHA. Durable bleeding protection and effective management of chronic pain have the potential to address unmet treatment needs in this population.

An Adaptive Three-Arm Comparative Clinical Endpoint Bioequivalence Study Design With Unblinded Sample Size Re-Estimation and Optimized Allocation Ratio.

A three-arm comparative clinical endpoint bioequivalence (BE) study is often used to establish bioequivalence (BE) between a locally acting generic drug (T) and reference drug (R), where superiority needs to be established for T and R over Placebo (P) and equivalence needs to be established for T vs. R. Sometimes, when study design parameters are uncertain, a fixed design study may be under- or over-powered and result in study failure or unnecessary cost. In this paper, we propose a two-stage adaptive clinical endpoint BE study with unblinded sample size re-estimation, standard or maximum combination method, optimized allocation ratio, optional re-estimation of the effect size based on likelihood estimation, and optional re-estimation of the R and P treatment means at interim analysis, which have not been done previously. Our proposed method guarantees control of Type 1 error rate analytically. It helps to reduce the average sample size when the original fixed design is overpowered and increases the sample size and power when the original study and group sequential design are under-powered. Our proposed adaptive design can help generic drug sponsors cut cost and improve success rate, making clinical study endpoint BE studies more affordable and more generic drugs accessible to the public.

Matching-adjusted indirect comparison of bleeding outcomes in severe haemophilia A: Comparing valoctocogene roxaparvovec gene therapy, emicizumab prophylaxis, and FVIII replacement prophylaxis.

INTRODUCTION: Head-to-head evaluation of valoctocogene roxaparvovec, the first gene therapy approved for haemophilia A, with emicizumab is not available. Therefore, phase 3 trial data were indirectly compared. AIM: To compare bleeding rates in trials evaluating 6 × 1013  vg/kg valoctocogene roxaparvovec (GENEr8-1; NCT03370913), 1.5 mg/kg emicizumab dosed every week (HAVEN 3; NCT02847637), and FVIII prophylaxis (270-902) in participants with severe haemophilia A (FVIII ≤1 IU/dL). METHODS: Valoctocogene roxaparvovec versus emicizumab and FVIII prophylaxis as used in 270-902 versus emicizumab cross-trial comparisons were performed using matching-adjusted indirect comparison (MAIC). Individual participant data from GENEr8-1 and 270-902 were weighted to equalise aggregate participant baseline characteristics from HAVEN 3. After MAIC weighting, annualised bleeding rates (ABR) and proportions of participants without bleeds were compared for treated bleeds, all bleeds, treated joint bleeds, and treated spontaneous bleeds. RESULTS: After MAIC weighting, ABR for all bleeds was statistically significantly lower with valoctocogene roxaparvovec than emicizumab (rate ratio [95% CI], .55 [.33-.93]). Additionally, significantly higher proportions of participants had no treated joint bleeds (odds ratio [95% CI], 2.75 [1.20-6.31]) and no treated bleeds (3.25 [1.53-6.90]) with valoctocogene roxaparvovec versus emicizumab. When compared with the mainly standard half-life FVIII prophylaxis regimens in 270-902, mean ABRs (except for all bleeds) were significantly lower, and significantly higher proportions reported 0 bleeds for all outcomes with emicizumab. CONCLUSION: Valoctocogene roxaparvovec provided generally lower bleeding rates and higher probability of no bleeds, including treated joint bleeds, than emicizumab. Emicizumab was more effective than FVIII prophylaxis regimens used in 270-902.

pH and Proton Sensor GPR65 Determine Susceptibility to Atopic Dermatitis.

pH sensing by GPR65 regulates various inflammatory conditions, but its role in skin remains unknown. In this study, we performed a phenome-wide association study and report that the T allele of GPR65-intronic single-nucleotide polymorphism rs8005161, which reduces GPR65 signaling, showed a significant association with atopic dermatitis, in addition to inflammatory bowel diseases and asthma, as previously reported. Consistent with this genetic association in humans, we show that deficiency of GPR65 in mice resulted in markedly exacerbated disease in the MC903 experimental model of atopic dermatitis. Deficiency of GPR65 also increased neutrophil migration in vitro. Moreover, GPR65 deficiency in mice resulted in higher expression of the inflammatory cytokine TNF-α by T cells. In humans, CD4+ T cells from rs8005161 heterozygous individuals expressed higher levels of TNF-α after PMA/ionomycin stimulation, particularly under pH 6 conditions. pH sensing by GPR65 appears to be important for regulating the pathogenesis of atopic dermatitis.

Healthcare utilization of patients with hereditary angioedema treated with lanadelumab and subcutaneous C1-inhibitor concentrate.

Background: New hereditary angioedema (HAE) treatments have become available in recent years for the treatment of HAE due to C1-inhibitor (C1-INH) deficiency, including two subcutaneous (SC) options: a monoclonal antibody (lanadelumab) and a plasma-derived C1-INH concentrate (SC-C1-INH). Limited real-world data on these therapies have been reported. Objective: The objective was to describe new users of lanadelumab and SC-C1-INH, including demographics, healthcare resource utilization (HCRU), costs, and treatment patterns before and after beginning treatment. Methods: This was a retrospective cohort study that used an administrative claims data base. Two mutually exclusive cohorts of adult (ages ≥18 years) new users of lanadelumab or SC-C1-INH with ≥180 days of continuous use were identified. HCRU, costs, and treatment patterns were assessed in the 180-day period before the index date (new treatment use) and up to 365 days after the index date. HCRU and costs were calculated as annualized rates. Results: Forty-seven patients who used lanadelumab and 38 patients who used SC-C1-INH were identified. The most frequently used on-demand HAE treatments at baseline were the same for both cohorts: bradykinin B2 antagonists (48.9% of the patients on lanadelumab, 52.6% of the patients on SC-C1-INH) and C1-INHs (40.4% of the patients on lanadelumab, 57.9% of the patients on SC-C1-INH). More than 33% of the patients continued to fill on-demand medications after treatment initiation. Annualized angioedema-associated emergency department visits and hospitalizations decreased after initiation of treatment, from 1.8 to 0.6 for the patients on lanadelumab and from 1.3 to 0.5 for the patients on SC-C1-INH. Annualized total healthcare costs after treatment initiation in the database were $866,639 and $734,460 for the lanadelumab and SC-C1-INH cohorts, respectively. Pharmacy costs accounted for >95% of these total costs. Conclusion: Although HCRU decreased after the initiation of treatment, angioedema-associated emergency department visits and hospitalizations and on-demand treatment fills were not completely eliminated. This indicates ongoing disease and treatment burden despite use of modern HAE medicines.

Age-of-onset information helps identify 76 genetic variants associated with allergic disease.

Risk factors that contribute to inter-individual differences in the age-of-onset of allergic diseases are poorly understood. The aim of this study was to identify genetic risk variants associated with the age at which symptoms of allergic disease first develop, considering information from asthma, hay fever and eczema. Self-reported age-of-onset information was available for 117,130 genotyped individuals of European ancestry from the UK Biobank study. For each individual, we identified the earliest age at which asthma, hay fever and/or eczema was first diagnosed and performed a genome-wide association study (GWAS) of this combined age-of-onset phenotype. We identified 50 variants with a significant independent association (P<3x10-8) with age-of-onset. Forty-five variants had comparable effects on the onset of the three individual diseases and 38 were also associated with allergic disease case-control status in an independent study (n = 222,484). We observed a strong negative genetic correlation between age-of-onset and case-control status of allergic disease (rg = -0.63, P = 4.5x10-61), indicating that cases with early disease onset have a greater burden of allergy risk alleles than those with late disease onset. Subsequently, a multivariate GWAS of age-of-onset and case-control status identified a further 26 associations that were missed by the univariate analyses of age-of-onset or case-control status only. Collectively, of the 76 variants identified, 18 represent novel associations for allergic disease. We identified 81 likely target genes of the 76 associated variants based on information from expression quantitative trait loci (eQTL) and non-synonymous variants, of which we highlight ADAM15, FOSL2, TRIM8, BMPR2, CD200R1, PRKCQ, NOD2, SMAD4, ABCA7 and UBE2L3. Our results support the notion that early and late onset allergic disease have partly distinct genetic architectures, potentially explaining known differences in pathophysiology between individuals.

Multitrait genetic association analysis identifies 50 new risk loci for gastro-oesophageal reflux, seven new loci for Barrett's oesophagus and provides insights into clinical heterogeneity in reflux diagnosis.

OBJECTIVE: Gastro-oesophageal reflux disease (GERD) has heterogeneous aetiology primarily attributable to its symptom-based definitions. GERD genome-wide association studies (GWASs) have shown strong genetic overlaps with established risk factors such as obesity and depression. We hypothesised that the shared genetic architecture between GERD and these risk factors can be leveraged to (1) identify new GERD and Barrett's oesophagus (BE) risk loci and (2) explore potentially heterogeneous pathways leading to GERD and oesophageal complications. DESIGN: We applied multitrait GWAS models combining GERD (78 707 cases; 288 734 controls) and genetically correlated traits including education attainment, depression and body mass index. We also used multitrait analysis to identify BE risk loci. Top hits were replicated in 23andMe (462 753 GERD cases, 24 099 BE cases, 1 484 025 controls). We additionally dissected the GERD loci into obesity-driven and depression-driven subgroups. These subgroups were investigated to determine how they relate to tissue-specific gene expression and to risk of serious oesophageal disease (BE and/or oesophageal adenocarcinoma, EA). RESULTS: We identified 88 loci associated with GERD, with 59 replicating in 23andMe after multiple testing corrections. Our BE analysis identified seven novel loci. Additionally we showed that only the obesity-driven GERD loci (but not the depression-driven loci) were associated with genes enriched in oesophageal tissues and successfully predicted BE/EA. CONCLUSION: Our multitrait model identified many novel risk loci for GERD and BE. We present strong evidence for a genetic underpinning of disease heterogeneity in GERD and show that GERD loci associated with depressive symptoms are not strong predictors of BE/EA relative to obesity-driven GERD loci.

Insights into the genetic basis of retinal detachment.

Retinal detachment (RD) is a serious and common condition, but genetic studies to date have been hampered by the small size of the assembled cohorts. In the UK Biobank data set, where RD was ascertained by self-report or hospital records, genetic correlations between RD and high myopia or cataract operation were, respectively, 0.46 (SE = 0.08) and 0.44 (SE = 0.07). These correlations are consistent with known epidemiological associations. Through meta-analysis of genome-wide association studies using UK Biobank RD cases (N = 3 977) and two cohorts, each comprising ~1 000 clinically ascertained rhegmatogenous RD patients, we uncovered 11 genome-wide significant association signals. These are near or within ZC3H11B, BMP3, COL22A1, DLG5, PLCE1, EFEMP2, TYR, FAT3, TRIM29, COL2A1 and LOXL1. Replication in the 23andMe data set, where RD is self-reported by participants, firmly establishes six RD risk loci: FAT3, COL22A1, TYR, BMP3, ZC3H11B and PLCE1. Based on the genetic associations with eye traits described to date, the first two specifically impact risk of a RD, whereas the last four point to shared aetiologies with macular condition, myopia and glaucoma. Fine-mapping prioritized the lead common missense variant (TYR S192Y) as causal variant at the TYR locus and a small set of credible causal variants at the FAT3 locus. The larger study size presented here, enabled by resources linked to health records or self-report, provides novel insights into RD aetiology and underlying pathological pathways.

Genetic Architectures of Childhood- and Adult-Onset Asthma Are Partly Distinct.

The extent to which genetic risk factors are shared between childhood-onset (COA) and adult-onset (AOA) asthma has not been estimated. On the basis of data from the UK Biobank study (n = 447,628), we found that the variance in disease liability explained by common variants is higher for COA (onset at ages between 0 and 19 years; h2g = 25.6%) than for AOA (onset at ages between 20 and 60 years; h2g = 10.6%). The genetic correlation (rg) between COA and AOA was 0.67. Variation in age of onset among COA-affected individuals had a low heritability (h2g = 5%), which we confirmed in independent studies and also among AOA-affected individuals. To identify subtype-specific genetic associations, we performed a genome-wide association study (GWAS) in the UK Biobank for COA (13,962 affected individuals) and a separate GWAS for AOA (26,582 affected individuals) by using a common set of 300,671 controls for both studies. We identified 123 independent associations for COA and 56 for AOA (37 overlapped); of these, 98 and 34, respectively, were reproducible in an independent study (n = 262,767). Collectively, 28 associations were not previously reported. For 96 COA-associated variants, including five variants that represent COA-specific risk factors, the risk allele was more common in COA- than in AOA-affected individuals. Conversely, we identified three variants that are stronger risk factors for AOA. Variants associated with obesity and smoking had a stronger contribution to the risk of AOA than to the risk of COA. Lastly, we identified 109 likely target genes of the associated variants, primarily on the basis of correlated expression quantitative trait loci (up to n = 31,684). GWAS informed by age of onset can identify subtype-specific risk variants, which can help us understand differences in pathophysiology between COA and AOA and so can be informative for drug development.

Findings from the Morquio A Registry Study (MARS) after 6 years: Long-term outcomes of MPS IVA patients treated with elosulfase alfa.

BACKGROUND: The Morquio A Registry Study (MARS) is an ongoing, multinational, observational study of patients with MPS IVA. Key objectives of MARS are to characterize the heterogeneity and natural history of disease and to evaluate long-term effectiveness and safety of elosulfase alfa enzyme replacement therapy (ERT). Enrollment began in September 2014; data on medical history, clinical outcomes, and safety assessments are collected as part of routine care. RESULTS: As of February 2021, 381 subjects from 17 countries had enrolled in MARS: 58 ERT-naïve subjects and 323 ERT-treated subjects (≥1 infusion), with a mean ERT exposure of 5.5 years (SD 2.8) and median age at first ERT treatment of 9.8 years. ERT-treated subjects were younger at diagnosis (median 3.4 vs 6.5 years) relative to ERT-naïve subjects. Among ERT-treated subjects, urinary keratan sulfate (uKS) levels declined from pre-ERT baseline to last follow-up on treatment (mean % change [95% confidence interval]: -52.5% [-57.5%, -47.4%]; n = 115) and 6-min walk test distance remained stable (mean change: -6.1 [-27.6, 15.5] m; n = 131) over a mean follow-up of 5.5 years. Forced expiratory volume in 1 s (FEV1) and forced vital capacity (FVC) increased in subjects who were < 18 years of age at ERT initiation (mean change: +0.3 [0.1, 0.4] L and + 0.4 [0.3, 0.5] L; mean follow-up: ∼6 years; n = 82) and were stable in subjects ≥18 years (mean change: 0.0 [-0.0, 0.1] L and 0.0 [-0.1, 0.1] L; mean follow-up: 4.6 years; n = 38). Overall, 148 (47.1%) ERT-treated subjects experienced ≥1 adverse event (AE) and 110 subjects (35%) reported ≥1 serious AE. Drug-related AEs were reported in 39 (12.4%) subjects; the most common were hypersensitivity (9 subjects [2.9%]), urticaria (8 subjects [2.5%]), and pyrexia (7 subjects [2.2%]). CONCLUSIONS: MARS is the longest and largest observational study of MPS IVA patients to date, with a heterogenous population that is representative of the MPS IVA population overall. Data collected over the first 6 years of MARS provide real-world evidence for long-term stabilization of endurance and respiratory function among ERT-treated patients, with no new safety concerns identified.

SARS-CoV-2 (COVID-19) structural and evolutionary dynamicome: Insights into functional evolution and human genomics.

The pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has challenged the speed at which laboratories can discover the viral composition and study health outcomes. The small ∼30-kb ssRNA genome of coronaviruses makes them adept at cross-species spread while enabling a robust understanding of all of the proteins the viral genome encodes. We have employed protein modeling, molecular dynamics simulations, evolutionary mapping, and 3D printing to gain a full proteome- and dynamicome-level understanding of SARS-CoV-2. We established the Viral Integrated Structural Evolution Dynamic Database (VIStEDD at RRID:SCR_018793) to facilitate future discoveries and educational use. Here, we highlight the use of VIStEDD for nsp6, nucleocapsid (N), and spike (S) surface glycoprotein. For both nsp6 and N, we found highly conserved surface amino acids that likely drive protein-protein interactions. In characterizing viral S protein, we developed a quantitative dynamics cross-correlation matrix to gain insights into its interactions with the angiotensin I-converting enzyme 2 (ACE2)-solute carrier family 6 member 19 (SLC6A19) dimer. Using this quantitative matrix, we elucidated 47 potential functional missense variants from genomic databases within ACE2/SLC6A19/transmembrane serine protease 2 (TMPRSS2), warranting genomic enrichment analyses in SARS-CoV-2 patients. These variants had ultralow frequency but existed in males hemizygous for ACE2. Two ACE2 noncoding variants (rs4646118 and rs143185769) present in ∼9% of individuals of African descent may regulate ACE2 expression and may be associated with increased susceptibility of African Americans to SARS-CoV-2. We propose that this SARS-CoV-2 database may aid research into the ongoing pandemic.

Genomic, transcriptomic, and protein landscape profile of CFTR and cystic fibrosis.

Cystic Fibrosis (CF) is caused most often by removal of amino acid 508 (Phe508del, deltaF508) within CFTR, yet dozens of additional CFTR variants are known to give rise to CF and many variants in the genome are known to contribute to CF pathology. To address CFTR coding variants, we developed a sequence-to-structure-to-dynamic matrix for all amino acids of CFTR using 233 vertebrate species, CFTR structure within a lipid membrane, and 20 ns of molecular dynamic simulation to assess known variants from the CFTR1, CFTR2, ClinVar, TOPmed, gnomAD, and COSMIC databases. Surprisingly, we identify 18 variants of uncertain significance within CFTR from diverse populations that are heritable and a likely cause of CF that have been understudied due to nonexistence in Caucasian populations. In addition, 15 sites within the genome are known to modulate CF pathology, where we have identified one genome region (chr11:34754985-34836401) that contributes to CF through modulation of expression of a noncoding RNA in epithelial cells. These 15 sites are just the beginning of understanding comodifiers of CF, where utilization of eQTLs suggests many additional genomics of CFTR expressing cells that can be influenced by genomic background of CFTR variants. This work highlights that many additional insights of CF genetics are needed, particularly as pharmaceutical interventions increase in the coming years.

Evaluating Lipid-Lowering Drug Targets for Parkinson's Disease Prevention with Mendelian Randomization.

Long-term exposure to lipid-lowering drugs might affect Parkinson's disease (PD) risk. We conducted Mendelian randomization analyses where genetic variants indexed expected effects of modulating lipid-lowering drug targets on PD. Statin exposure was not predicted to increase PD risk, although results were not precise enough to support benefits for prevention clearly (odds ratio [OR] = 0.83; 95% confidence interval [CI] = 0.65, 1.07). Other target results were null, except for variants indicating Apolipoprotein-A5 or Apolipoprotein-C3 inhibition might confer protection. These findings suggest peripheral lipid variation may not have a prominent role in PD etiology, but some related drug targets could influence PD via alternate pathways. ANN NEUROL 2020;88:1043-1047.

Treatment patterns among non-active users of maintenance asthma medication in the United Kingdom: a retrospective cohort study in the Clinical Practice Research Datalink.

Objective: To describe patient characteristics, treatment patterns and healthcare utilization (HCU) of non-active users of maintenance asthma medications in the United Kingdom.Methods: Retrospective, cohort analysis of patients with asthma, aged ≥ 6 years who were non-active users of maintenance therapy (no prescription for inhaled corticosteroids (ICS), combined ICS/long-acting beta agonists (ICS/LABA) or 'other' bronchodilatory therapies in last 12 months) were identified in the Clinical Practice Research Datalink (2012-2015) and followed-up for 2 years after a new prescription for an asthma maintenance medication (index date). Patient characteristics, most common maintenance treatment sequences and HCU were described.Results: 55,293 patients were identified (ICS: 46,297, ICS/LABA: 8,367; Other: 629). Mean age was 37 years and 56% were female. During follow-up, the most common treatment sequences across groups implied intermittent use, comprising periods of maintenance therapy interspersed with maintenance-free periods. During year 1 and year 2 of follow-up, the proportion of patients prescribed OCS was 19% and 13%, prescribed ≥ 4 short-acting bronchodilators (SABD) was 24% and 19%, having ≥ 3 asthma-related primary care consultations/year was 59% and 36% and experiencing ≥ 1 exacerbation/year was 15% and 11%, respectively.Conclusions: In previously non-active users of asthma maintenance medication subsequently commenced on maintenance therapy, intermittent use was common during the 2-year follow-up despite the potential need for regular use as evidenced by patient HCU and SABD usage patterns. This highlights the need for regular patient assessment and education on medication adherence to ensure appropriateness of prescribing to maintain asthma control.

Whole exome sequencing reveals HSPA1L as a genetic risk factor for spontaneous preterm birth.

Preterm birth is a leading cause of morbidity and mortality in infants. Genetic and environmental factors play a role in the susceptibility to preterm birth, but despite many investigations, the genetic basis for preterm birth remain largely unknown. Our objective was to identify rare, possibly damaging, nucleotide variants in mothers from families with recurrent spontaneous preterm births (SPTB). DNA samples from 17 Finnish mothers who delivered at least one infant preterm were subjected to whole exome sequencing. All mothers were of northern Finnish origin and were from seven multiplex families. Additional replication samples of European origin consisted of 93 Danish sister pairs (and two sister triads), all with a history of a preterm delivery. Rare exonic variants (frequency <1%) were analyzed to identify genes and pathways likely to affect SPTB susceptibility. We identified rare, possibly damaging, variants in genes that were common to multiple affected individuals. The glucocorticoid receptor signaling pathway was the most significant (p<1.7e-8) with genes containing these variants in a subgroup of ten Finnish mothers, each having had 2-4 SPTBs. This pathway was replicated among the Danish sister pairs. A gene in this pathway, heat shock protein family A (Hsp70) member 1 like (HSPA1L), contains two likely damaging missense alleles that were found in four different Finnish families. One of the variants (rs34620296) had a higher frequency in cases compared to controls (0.0025 vs. 0.0010, p = 0.002) in a large preterm birth genome-wide association study (GWAS) consisting of mothers of general European ancestry. Sister pairs in replication samples also shared rare, likely damaging HSPA1L variants. Furthermore, in silico analysis predicted an additional phosphorylation site generated by rs34620296 that could potentially affect chaperone activity or HSPA1L protein stability. Finally, in vitro functional experiment showed a link between HSPA1L activity and decidualization. In conclusion, rare, likely damaging, variants in HSPA1L were observed in multiple families with recurrent SPTB.

Correction: Whole exome sequencing reveals HSPA1L as a genetic risk factor for spontaneous preterm birth.

[This corrects the article DOI: 10.1371/journal.pgen.1007394.].

Assessment of rosacea symptom severity by genome-wide association study and expression analysis highlights immuno-inflammatory and skin pigmentation genes.

Rosacea is a common, chronic skin disease of variable severity with limited treatment options. The cause of rosacea is unknown, but it is believed to be due to a combination of hereditary and environmental factors. Little is known about the genetics of the disease. We performed a genome-wide association study (GWAS) of rosacea symptom severity with data from 73 265 research participants of European ancestry from the 23andMe customer base. Seven loci had variants associated with rosacea at the genome-wide significance level (P < 5 × 10-8). Further analyses highlighted likely gene regions or effector genes including IRF4 (P = 1.5 × 10-17), a human leukocyte antigen (HLA) region flanked by PSMB9 and HLA-DMB (P = 2.2 × 10-15), HERC2-OCA2 (P = 4.2 × 10-12), SLC45A2 (P = 1.7 × 10-10), IL13 (P = 2.8 × 10-9), a region flanked by NRXN3 and DIO2 (P = 4.1 × 10-9), and a region flanked by OVOL1and SNX32 (P = 1.2 × 10-8). All associations with rosacea were novel except for the HLA locus. Two of these loci (HERC-OCA2 and SLC45A2) and another precedented variant (rs1805007 in melanocortin 1 receptor) with an association P value just below the significance threshold (P = 1.3 × 10-7) have been previously associated with skin phenotypes and pigmentation, two of these loci are linked to immuno-inflammation phenotypes (IL13 and PSMB9-HLA-DMA) and one has been associated with both categories (IRF4). Genes within three loci (PSMB9-HLA-DMA, HERC-OCA2 and NRX3-DIO2) were differentially expressed in a previously published clinical rosacea transcriptomics study that compared lesional to non-lesional samples. The identified loci provide specificity of inflammatory mechanisms in rosacea, and identify potential pathways for therapeutic intervention.

Multiethnic GWAS Reveals Polygenic Architecture of Earlobe Attachment.

The genetic basis of earlobe attachment has been a matter of debate since the early 20th century, such that geneticists argue both for and against polygenic inheritance. Recent genetic studies have identified a few loci associated with the trait, but large-scale analyses are still lacking. Here, we performed a genome-wide association study of lobe attachment in a multiethnic sample of 74,660 individuals from four cohorts (three with the trait scored by an expert rater and one with the trait self-reported). Meta-analysis of the three expert-rater-scored cohorts revealed six associated loci harboring numerous candidate genes, including EDAR, SP5, MRPS22, ADGRG6 (GPR126), KIAA1217, and PAX9. The large self-reported 23andMe cohort recapitulated each of these six loci. Moreover, meta-analysis across all four cohorts revealed a total of 49 significant (p < 5 × 10-8) loci. Annotation and enrichment analyses of these 49 loci showed strong evidence of genes involved in ear development and syndromes with auricular phenotypes. RNA sequencing data from both human fetal ear and mouse second branchial arch tissue confirmed that genes located among associated loci showed evidence of expression. These results provide strong evidence for the polygenic nature of earlobe attachment and offer insights into the biological basis of normal and abnormal ear development.

Genetic Associations with Gestational Duration and Spontaneous Preterm Birth.

BACKGROUND: Despite evidence that genetic factors contribute to the duration of gestation and the risk of preterm birth, robust associations with genetic variants have not been identified. We used large data sets that included the gestational duration to determine possible genetic associations. METHODS: We performed a genomewide association study in a discovery set of samples obtained from 43,568 women of European ancestry using gestational duration as a continuous trait and term or preterm (<37 weeks) birth as a dichotomous outcome. We used samples from three Nordic data sets (involving a total of 8643 women) to test for replication of genomic loci that had significant genomewide association (P<5.0×10-8) or an association with suggestive significance (P<1.0×10-6) in the discovery set. RESULTS: In the discovery and replication data sets, four loci (EBF1, EEFSEC, AGTR2, and WNT4) were significantly associated with gestational duration. Functional analysis showed that an implicated variant in WNT4 alters the binding of the estrogen receptor. The association between variants in ADCY5 and RAP2C and gestational duration had suggestive significance in the discovery set and significant evidence of association in the replication sets; these variants also showed genomewide significance in a joint analysis. Common variants in EBF1, EEFSEC, and AGTR2 showed association with preterm birth with genomewide significance. An analysis of mother-infant dyads suggested that these variants act at the level of the maternal genome. CONCLUSIONS: In this genomewide association study, we found that variants at the EBF1, EEFSEC, AGTR2, WNT4, ADCY5, and RAP2C loci were associated with gestational duration and variants at the EBF1, EEFSEC, and AGTR2 loci with preterm birth. Previously established roles of these genes in uterine development, maternal nutrition, and vascular control support their mechanistic involvement. (Funded by the March of Dimes and others.).