A catalog of genetic loci associated with kidney function from analyses of a million individuals
(2019)
Chronic kidney disease (CKD) is responsible for a public health burden with multi-systemic complications. Through transancestry meta-analysis of genome-wide association studies of estimated glomerular filtration rate (eGFR) and independent replication (n = 1,046,070), we identified 264 associated loci (166 new). Of these,147 were likely to be relevant for kidney function on the basis of associations with the alternative kidney function marker blood urea nitrogen (n = 416,178). Pathway and enrichment analyses, including mouse models with renal phenotypes, support the kidney as the main target organ. A genetic risk score for lower eGFR was associated with clinically diagnosed CKD in 452,264 independent individuals. Colocalization analyses of associations with eGFR among 783,978 European-ancestry individuals and gene expression across 46 human tissues, including tubulo-interstitial and glomerular kidney compartments, identified 17 genes differentially expressed in kidney. Fine-mapping highlighted missense driver variants in 11 genes and kidney-specific regulatory variants. These results provide a comprehensive priority list of molecular targets for translational research.
Rapid decline of glomerular filtration rate estimated from creatinine (eGFRcrea) is associated with severe clinical endpoints. In contrast to cross-sectionally assessed eGFRcrea, the genetic basis for rapid eGFRcrea decline is largely unknown. To help define this, we meta-analyzed 42 genome-wide association studies from the Chronic Kidney Diseases Genetics Consortium and United Kingdom Biobank to identify genetic loci for rapid eGFRcrea decline. Two definitions of eGFRcrea decline were used: 3 mL/min/1.73m(2)/year or more ("Rapid3"; encompassing 34,874 cases, 107,090 controls) and eGFRcrea decline 25% or more and eGFRcrea under 60 mL/min/1.73m(2) at follow-up among those with eGFRcrea 60 mL/min/1.73m(2) or more at baseline ("CKDi25"; encompassing 19,901 cases, 175,244 controls). Seven independent variants were identified across six loci for Rapid3 and/or CKDi25: consisting of five variants at four loci with genome-wide significance (near UMOD-PDILT (2), PRKAG2, WDR72, OR2S2) and two variants among 265 known eGFRcrea variants (near GATM, LARP4B). All these loci were novel for Rapid3 and/or CKDi25 and our bioinformatic follow-up prioritized variants and genes underneath these loci. The OR2S2 locus is novel for any eGFRcrea trait including interesting candidates. For the five genome-wide significant lead variants, we found supporting effects for annual change in blood urea nitrogen or cystatin-based eGFR, but not for GATM or (LARP4B). Individuals at high compared to those at low genetic risk (8-14 vs. 0-5 adverse alleles) had a 1.20-fold increased risk of acute kidney injury (95% confidence interval 1.08-1.33). Thus, our identified loci for rapid kidney function decline may help prioritize therapeutic targets and identify mechanisms and individuals at risk for sustained deterioration of kidney function.
Rapid decline of glomerular filtration rate estimated from creatinine (eGFRcrea) is associated with severe clinical endpoints. In contrast to cross-sectionally assessed eGFRcrea, the genetic basis for rapid eGFRcrea decline is largely unknown. To help define this, we meta-analyzed 42 genome-wide association studies from the Chronic Kidney Diseases Genetics Consortium and United Kingdom Biobank to identify genetic loci for rapid eGFRcrea decline. Two definitions of eGFRcrea decline were used: 3 mL/min/1.73m(2)/year or more ("Rapid3"; encompassing 34,874 cases, 107,090 controls) and eGFRcrea decline 25% or more and eGFRcrea under 60 mL/min/1.73m(2) at follow-up among those with eGFRcrea 60 mL/min/1.73m(2) or more at baseline ("CKDi25"; encompassing 19,901 cases, 175,244 controls). Seven independent variants were identified across six loci for Rapid3 and/or CKDi25: consisting of five variants at four loci with genome-wide significance (near UMOD-PDILT (2), PRKAG2, WDR72, OR2S2) and two variants among 265 known eGFRcrea variants (near GATM, LARP4B). All these loci were novel for Rapid3 and/or CKDi25 and our bioinformatic follow-up prioritized variants and genes underneath these loci. The OR2S2 locus is novel for any eGFRcrea trait including interesting candidates. For the five genome-wide significant lead variants, we found supporting effects for annual change in blood urea nitrogen or cystatin-based eGFR, but not for GATM or (LARP4B). Individuals at high compared to those at low genetic risk (8-14 vs. 0-5 adverse alleles) had a 1.20-fold increased risk of acute kidney injury (95% confidence interval 1.08-1.33). Thus, our identified loci for rapid kidney function decline may help prioritize therapeutic targets and identify mechanisms and individuals at risk for sustained deterioration of kidney function.
Distinct cellular roles for PDCD10 define a gut-brain axis in cerebral cavernous malformation
(2019)
Cerebral cavernous malformation (CCM) is a genetic, cerebrovascular disease. Familial CCM is caused by genetic mutations in KRIT1, CCM2, or PDCD10. Disease onset is earlier and more severe in individuals with PDCD10 mutations. Recent studies have shown that lesions arise from excess mitogen-activated protein kinase kinase kinase 3 (MEKK3) signaling downstream of Toll-like receptor 4 (TLR4) stimulation by lipopolysaccharide derived from the gut microbiome. These findings suggest a gut-brain CCM disease axis but fail to define it or explain the poor prognosis of patients with PDCD10 mutations. Here, we demonstrate that the gut barrier is a primary determinant of CCM disease course, independent of microbiome configuration, that explains the increased severity of CCM disease associated with PDCD10 deficiency. Chemical disruption of the gut barrier with dextran sulfate sodium augments CCM formation in a mouse model, as does genetic loss of Pdcd10, but not Krit1, in gut epithelial cells. Loss of gut epithelial Pdcd10 results in disruption of the colonic mucosal barrier. Accordingly, loss of Mucin-2 or exposure to dietary emulsifiers that reduce the mucus barrier increases CCM burden analogous to loss of Pdcd10 in the gut epithelium. Last, we show that treatment with dexamethasone potently inhibits CCM formation in mice because of the combined effect of action at both brain endothelial cells and gut epithelial cells. These studies define a gut-brain disease axis in an experimental model of CCM in which a single gene is required for two critical components: gut epithelial function and brain endothelial signaling.
VERITAS and Fermi-LAT Observations of TeV Gamma-Ray Sources Discovered by HAWC in the 2HWC Catalog
(2018)
The High Altitude Water Cherenkov (HAWC) collaboration recently published their 2HWC catalog, listing 39 very high energy (VHE; >100 GeV) gamma-ray sources based on 507 days of observation. Among these, 19 sources are not associated with previously known teraelectronvolt (TeV) gamma-ray sources. We have studied 14 of these sources without known counterparts with VERITAS and Fermi-LAT. VERITAS detected weak gamma-ray emission in the 1 TeV-30 TeV band in the region of DA 495, a pulsar wind nebula coinciding with 2HWC J1953+294, confirming the discovery of the source by HAWC. We did not find any counterpart for the selected 14 new HAWC sources from our analysis of Fermi-LAT data for energies higher than 10 GeV. During the search, we detected gigaelectronvolt (GeV) gamma-ray emission coincident with a known TeV pulsar wind nebula, SNR G54.1+0.3 (VER J1930+188), and a 2HWC source, 2HWC J1930+188. The fluxes for isolated, steady sources in the 2HWC catalog are generally in good agreement with those measured by imaging atmospheric Cherenkov telescopes. However, the VERITAS fluxes for SNR G54.1+0.3, DA 495, and TeV J2032+4130 are lower than those measured by HAWC, and several new HAWC sources are not detected by VERITAS. This is likely due to a change in spectral shape, source extension, or the influence of diffuse emission in the source region.
During the summer of 2013, a 4-month spectroscopic campaign took place to observe the variabilities in three Wolf-Rayet stars. The spectroscopic data have been analysed for WR 134 (WN6b), to better understand its behaviour and long-term periodicity, which we interpret as arising from corotating interaction regions (CIRs) in the wind. By analysing the variability of the He ii lambda 5411 emission line, the previously identified period was refined to P = 2.255 +/- 0.008 (s.d.) d. The coherency time of the variability, which we associate with the lifetime of the CIRs in the wind, was deduced to be 40 +/- 6 d, or similar to 18 cycles, by cross-correlating the variability patterns as a function of time. When comparing the phased observational grey-scale difference images with theoretical grey-scales previously calculated from models including CIRs in an optically thin stellar wind, we find that two CIRs were likely present. A separation in longitude of Delta I center dot a parts per thousand integral 90A degrees was determined between the two CIRs and we suggest that the different maximum velocities that they reach indicate that they emerge from different latitudes. We have also been able to detect observational signatures of the CIRs in other spectral lines (C iv lambda lambda 5802,5812 and He i lambda 5876). Furthermore, a DAC was found to be present simultaneously with the CIR signatures detected in the He i lambda 5876 emission line which is consistent with the proposed geometry of the large-scale structures in the wind. Small-scale structures also show a presence in the wind, simultaneously with the larger scale structures, showing that they do in fact co-exist.
We present the results of a multi-wavelength campaign targeting the blazar 1ES 1218+30.4 with observations with the 1.3 m McGraw-Hill optical telescope, the Rossi X-ray Timing Explorer (RXTE), the Fermi Gamma-Ray Space Telescope, and the Very Energetic Radiation Imaging Telescope Array System (VERITAS). The RXTE and VERITAS observations were spread over a 13 day period and revealed clear evidence for flux variability, and a strong X-ray and gamma-ray flare on 2009 February 26 (MJD 54888). The campaign delivered a well-sampled broadband energy spectrum with simultaneous RXTE and VERITAS very high energy (VHE, > 100 GeV) observations, as well as contemporaneous optical and Fermi observations. The 1ES 1218+30.4 broadband energy spectrum-the first with simultaneous X-ray and VHE gamma-ray energy spectra-is of particular interest as the source is located at a high cosmological redshift for a VHE source (z = 0.182), leading to strong absorption of VHE gamma rays by photons from the optical/infrared extragalactic background light (EBL) via gamma VHE +gamma EBL -> e(+) e(-)pair-creation processes. We model the data with a one-zone synchrotron self-Compton (SSC) emission model and with the extragalactic absorption predicted by several recent EBL models. We find that the observations are consistent with the SSC scenario and all the EBL models considered in this work. We discuss observational and theoretical avenues to improve on the EBL constraints.
The selaginella genome identifies genetic changes associated with the evolution of vascular plants
(2011)
Vascular plants appeared similar to 410 million years ago, then diverged into several lineages of which only two survive: the euphyllophytes (ferns and seed plants) and the lycophytes. We report here the genome sequence of the lycophyte Selaginella moellendorffii (Selaginella), the first nonseed vascular plant genome reported. By comparing gene content in evolutionarily diverse taxa, we found that the transition from a gametophyte- to a sporophyte-dominated life cycle required far fewer new genes than the transition from a nonseed vascular to a flowering plant, whereas secondary metabolic genes expanded extensively and in parallel in the lycophyte and angiosperm lineages. Selaginella differs in posttranscriptional gene regulation, including small RNA regulation of repetitive elements, an absence of the trans-acting small interfering RNA pathway, and extensive RNA editing of organellar genes.
Birth weight variation is influenced by fetal and maternal genetic and non-genetic factors, and has been reproducibly associated with future cardio-metabolic health outcomes. In expanded genome-wide association analyses of own birth weight (n = 321,223) and offspring birth weight (n = 230,069 mothers), we identified 190 independent association signals (129 of which are novel). We used structural equation modeling to decompose the contributions of direct fetal and indirect maternal genetic effects, then applied Mendelian randomization to illuminate causal pathways. For example, both indirect maternal and direct fetal genetic effects drive the observational relationship between lower birth weight and higher later blood pressure: maternal blood pressure-raising alleles reduce offspring birth weight, but only direct fetal effects of these alleles, once inherited, increase later offspring blood pressure. Using maternal birth weight-lowering genotypes to proxy for an adverse intrauterine environment provided no evidence that it causally raises offspring blood pressure, indicating that the inverse birth weight-blood pressure association is attributable to genetic effects, and not to intrauterine programming.
The Early Growth Genetics (EGG) and EArly Genetics and Lifecourse Epidemiology (EAGLE) consortia
(2019)
The impact of many unfavorable childhood traits or diseases, such as low birth weight and mental disorders, is not limited to childhood and adolescence, as they are also associated with poor outcomes in adulthood, such as cardiovascular disease. Insight into the genetic etiology of childhood and adolescent traits and disorders may therefore provide new perspectives, not only on how to improve wellbeing during childhood, but also how to prevent later adverse outcomes. To achieve the sample sizes required for genetic research, the Early Growth Genetics (EGG) and EArly Genetics and Lifecourse Epidemiology (EAGLE) consortia were established. The majority of the participating cohorts are longitudinal population-based samples, but other cohorts with data on early childhood phenotypes are also involved. Cohorts often have a broad focus and collect(ed) data on various somatic and psychiatric traits as well as environmental factors. Genetic variants have been successfully identified for multiple traits, for example, birth weight, atopic dermatitis, childhood BMI, allergic sensitization, and pubertal growth. Furthermore, the results have shown that genetic factors also partly underlie the association with adult traits. As sample sizes are still increasing, it is expected that future analyses will identify additional variants. This, in combination with the development of innovative statistical methods, will provide detailed insight on the mechanisms underlying the transition from childhood to adult disorders. Both consortia welcome new collaborations. Policies and contact details are available from the corresponding authors of this manuscript and/or the consortium websites.
This review presents recommended nomenclature for the biosynthesis of ribosomally synthesized and post-translationally modified peptides (RiPPs), a rapidly growing class of natural products. The current knowledge regarding the biosynthesis of the >20 distinct compound classes is also reviewed, and commonalities are discussed.
We employ P to S converted waveforms to investigate effects of the hot mantle plume on seismic discontinuities of the crust and upper mantle. We observe the Moho at depths between 13 and 17 km, regionally covered by a strong shallow intracrustal converted phase. Coherent phases on the transverse component indicate either dipping interfaces, 3- D heterogeneities or lower crustal anisotropy. We find anomalies related to discontinuities in the upper mantle down to the transition zone evidently related to the hot mantle plume. Lithospheric thinning is confirmed in greater detail than previously reported by Li et al., and we determine the dimensions of the low-velocity zone within the asthenosphere with greater accuracy. Our study mainly focuses on the temperature-pressure dependent discontinuities of the upper mantle transition zone. Effects of the hot diapir on the depths of mineral phase transitions are verified at both major interfaces at 410 and 660 km. We determine a plume radius of about 200 km at the 660 km discontinuity with a core zone of about 120 km radius. The plume conduit is located southwest of Big Island. A conduit tilted in the northeast direction is required in the upper mantle to explain the observations. The determined positions of deflections of the discontinuities support the hypothesis of decoupled upper and lower mantle convection
The photochemical ring-opening of 1,3-cyclohexadiene imaged by ultrafast electron diffraction
(2019)
The ultrafast photoinduced ring-opening of 1,3-cyclohexadiene constitutes a textbook example of electrocyclic reactions in organic chemistry and a model for photobiological reactions in vitamin D synthesis. Although the relaxation from the photoexcited electronic state during the ring-opening has been investigated in numerous studies, the accompanying changes in atomic distance have not been resolved. Here we present a direct and unambiguous observation of the ring-opening reaction path on the femtosecond timescale and subangstrom length scale using megaelectronvolt ultrafast electron diffraction. We followed the carbon-carbon bond dissociation and the structural opening of the 1,3-cyclohexadiene ring by the direct measurement of time-dependent changes in the distribution of interatomic distances. We observed a substantial acceleration of the ring-opening motion after internal conversion to the ground state due to a steepening of the electronic potential gradient towards the product minima. The ring-opening motion transforms into rotation of the terminal ethylene groups in the photoproduct 1,3,5-hexatriene on the subpicosecond timescale.
Polymeric matrices mimicking multiple functions of the ECM are expected to enable a material induced regeneration of tissues. Here, we investigated the adipogenic differentiation of human adipose derived mesenchymal stem cells (hADSCs) in a 3D architectured gelatin based hydrogel (ArcGel) prepared from gelatin and L-lysine diisocyanate ethyl ester (LDI) in an one-step process, in which the formation of an open porous morphology and the chemical network formation were integrated. The ArcGel was designed to support adipose tissue regeneration with its 3D porous structure, high cell biocompatibility, and mechanical properties compatible with human subcutaneous adipose tissue. The ArcGel could support initial cell adhesion and survival of hADSCs. Under static culture condition, the cells could migrate into the inner part of the scaffold with a depth of 840 +/- 120 mu m after 4 days, and distributed in the whole scaffold (2mm in thickness) within 14 days. The cells proliferated in the scaffold and the fold increase of cell number after 7 days of culture was 2.55 +/- 0.08. The apoptotic rate of hADSCs in the scaffold was similar to that of cells maintained on tissue culture plates. When cultured in adipogenic induction medium, the hADSCs in the scaffold differentiated into adipocytes with a high efficiency (93 +/- 1%). Conclusively, this gelatin based 3D scaffold presented high cell compatibility for hADSC cultivation and differentiation, which could serve as a potential implant material in clinical applications for adipose tissue reparation and regeneration.
For over a hundred years, the "river sharks" of the genus Glyphis were only known from the type specimens of species that had been collected in the 19th century. They were widely considered extinct until populations of Glyphis-like sharks were rediscovered in remote regions of Borneo and Northern Australia at the end of the 20th century. However, the genetic affinities between the newly discovered Glyphis-like populations and the poorly preserved, original museum-type specimens have never been established. Here, we present the first (to our knowledge) fully resolved, complete phylogeny of Glyphis that includes both archival-type specimens and modern material. We used a sensitive DNA hybridization capture method to obtain complete mitochondrial genomes from all of our samples and show that three of the five described river shark species are probably conspecific and widely distributed in Southeast Asia. Furthermore we show that there has been recent gene flow between locations that are separated by large oceanic expanses. Our data strongly suggest marine dispersal in these species, overturning the widely held notion that river sharks are restricted to freshwater. It seems that species in the genus Glyphis are euryhaline with an ecology similar to the bull shark, in which adult individuals live in the ocean while the young grow up in river habitats with reduced predation pressure. Finally, we discovered a previously unidentified species within the genus Glyphis that is deeply divergent from all other lineages, underscoring the current lack of knowledge about the biodiversity and ecology of these mysterious sharks.
Nanoparticles can improve topical drug delivery: size, surface properties and flexibility of polymer nanoparticles are defining its interaction with the skin. Only few studies have explored skin penetration for one series of structurally related polymer particles with systematic alteration of material composition. Here, a series of rigid poly[acrylonitrile-co-(N-vinyl pyrrolidone)] model nanoparticles stably loaded with Nile Red or Rhodamin B, respectively, was comprehensively studied for biocompatibility and functionality. Surface properties were altered by varying the molar content of hydrophilic NVP from 0 to 24.1% and particle size ranged from 35 to 244 nm. Whereas irritancy and genotoxicity were not revealed, lipophilic and hydrophilic nanoparticles taken up by keratinocytes affected cell viability. Skin absorption of the particles into viable skin ex vivo was studied using Nile Red as fluorescent probe. Whilst an intact stratum corneum efficiently prevented penetration, almost complete removal of the horny layer allowed nanoparticles of smaller size and hydrophilic particles to penetrate into viable epidermis and dermis. Hence, systematic variations of nanoparticle properties allows gaining insights into critical criteria for biocompatibility and functionality of novel nanocarriers for topical drug delivery and risks associated with environmental exposure.
Derivatization of fullerene (C60) with branched aliphatic chains softens C60-based materials and enables the formation of thermotropic liquid crystals and room temperature nonvolatile liquids. This work demonstrates that by carefully tuning parameters such as type, number and substituent position of the branched chains, liquid crystalline C60 materials with mesophase temperatures suited for photovoltaic cell fabrication and room temperature nonvolatile liquid fullerenes with tunable viscosity can be obtained. In particular, compound 1, with branched chains, exhibits a smectic liquid crystalline phase extending from 84 °C to room temperature. Analysis of bulk heterojunction (BHJ) organic solar cells with a ca. 100 nm active layer of compound 1 and poly(3-hexylthiophene) (P3HT) as an electron acceptor and an electron donor, respectively, reveals an improved performance (power conversion efficiency, PCE: 1.6 ± 0.1%) in comparison with another compound, 10 (PCE: 0.5 ± 0.1%). The latter, in contrast to 1, carries linear aliphatic chains and thus forms a highly ordered solid lamellar phase at room temperature. The solar cell performance of 1 blended with P3HT approaches that of PCBM/P3HT for the same active layer thickness. This indicates that C60 derivatives bearing branched tails are a promising class of electron acceptors in soft (flexible) photovoltaic devices.
Derivatization of fullerene (C60) with branched aliphatic chains softens C60-based materials and enables the formation of thermotropic liquid crystals and room temperature nonvolatile liquids. This work demonstrates that by carefully tuning parameters such as type, number and substituent position of the branched chains, liquid crystalline C60 materials with mesophase temperatures suited for photovoltaic cell fabrication and room temperature nonvolatile liquid fullerenes with tunable viscosity can be obtained. In particular, compound 1, with branched chains, exhibits a smectic liquid crystalline phase extending from 84°C to room temperature. Analysis of bulk heterojunction (BHJ) organic solar cells with a ca. 100 nm active layer of compound 1 and poly(3-hexylthiophene) (P3HT) as an electron acceptor and an electron donor, respectively, reveals an improved performance (power conversion efficiency, PCE: 1.6 ñ 0.1%) in comparison with another compound, 10 (PCE: 0.5 ñ 0.1%). The latter, in contrast to 1, carries linear aliphatic chains and thus forms a highly ordered solid lamellar phase at room temperature. The solar cell performance of 1 blended with P3HT approaches that of PCBM/P3HT for the same active layer thickness. This indicates that C60 derivatives bearing branched tails are a promising class of electron acceptors in soft (flexible) photovoltaic devices.
Derivatization of fullerene (C-60) with branched aliphatic chains softens C-60-based materials and enables the formation of thermotropic liquid crystals and room temperature nonvolatile liquids. This work demonstrates that by carefully tuning parameters such as type, number and substituent position of the branched chains, liquid crystalline C-60 materials with mesophase temperatures suited for photovoltaic cell fabrication and room temperature nonvolatile liquid fullerenes with tunable viscosity can be obtained. In particular, compound 1, with branched chains, exhibits a smectic liquid crystalline phase extending from 84 degrees C to room temperature. Analysis of bulk heterojunction (BHJ) organic solar cells with a ca. 100 nm active layer of compound 1 and poly(3-hexylthiophene) (P3HT) as an electron acceptor and an electron donor, respectively, reveals an improved performance (power conversion efficiency, PCE: 1.6 + 0.1%) in comparison with another compound, 10 (PCE: 0.5 + 0.1%). The latter, in contrast to 1, carries linear aliphatic chains and thus forms a highly ordered solid lamellar phase at room temperature. The solar cell performance of 1 blended with P3HT approaches that of PCBM/P3HT for the same active layer thickness. This indicates that C-60 derivatives bearing branched tails are a promising class of electron acceptors in soft (flexible) photovoltaic devices.
The largest uncertainty in projections of future sea-level change results from the potentially changing dynamical ice discharge from Antarctica. Basal ice-shelf melting induced by a warming ocean has been identified as a major cause for additional ice flow across the grounding line. Here we attempt to estimate the uncertainty range of future ice discharge from Antarctica by combining uncertainty in the climatic forcing, the oceanic response and the ice-sheet model response. The uncertainty in the global mean temperature increase is obtained from historically constrained emulations with the MAGICC-6.0 (Model for the Assessment of Greenhouse gas Induced Climate Change) model. The oceanic forcing is derived from scaling of the subsurface with the atmospheric warming from 19 comprehensive climate models of the Coupled Model Intercomparison Project (CMIP-5) and two ocean models from the EU-project Ice2Sea. The dynamic ice-sheet response is derived from linear response functions for basal ice-shelf melting for four different Antarctic drainage regions using experiments from the Sea-level Response to Ice Sheet Evolution (SeaRISE) intercomparison project with five different Antarctic ice-sheet models. The resulting uncertainty range for the historic Antarctic contribution to global sea-level rise from 1992 to 2011 agrees with the observed contribution for this period if we use the three ice-sheet models with an explicit representation of ice-shelf dynamics and account for the time-delayed warming of the oceanic subsurface compared to the surface air temperature. The median of the additional ice loss for the 21st century is computed to 0.07 m (66% range: 0.02-0.14 m; 90% range: 0.0-0.23 m) of global sea-level equivalent for the low-emission RCP-2.6 (Representative Concentration Pathway) scenario and 0.09 m (66% range: 0.04-0.21 m; 90% range: 0.01-0.37 m) for the strongest RCP-8.5. Assuming no time delay between the atmospheric warming and the oceanic subsurface, these values increase to 0.09 m (66% range: 0.04-0.17 m; 90% range: 0.02-0.25 m) for RCP-2.6 and 0.15 m (66% range: 0.07-0.28 m; 90% range: 0.04-0.43 m) for RCP-8.5. All probability distributions are highly skewed towards high values. The applied ice-sheet models are coarse resolution with limitations in the representation of grounding-line motion. Within the constraints of the applied methods, the uncertainty induced from different ice-sheet models is smaller than that induced by the external forcing to the ice sheets.