Search

Climate change reduces snowpack, advances snowmelt phenology, drives summer warming, alters growing season precipitation regimes, and consequently modifies vegetation phenology in mountain systems. Elevational migrants track spatial variation in seasonal plant growth by moving between ranges at different elevations during spring, so climate-driven vegetation change may disrupt historic benefits of migration. Elevational migrants can furthermore cope with short-term environmental variability by undertaking brief vertical movements to refugia when sudden adverse conditions arise. We uncover drivers of fine-scale vertical movement variation during upland migration in an endangered alpine specialist, Sierra Nevada bighorn sheep (Ovis canadensis sierrae) using a 20-year study of GPS collar data collected from 311 unique individuals. We used integrated step-selection analysis to determine factors that promote vertical movements and drive selection of destinations following vertical movements. Our results reveal that relatively high temperatures consistently drive uphill movements, while precipitation likely drives downhill movements. Furthermore, bighorn select destinations at their peak annual biomass and maximal time since snowmelt. These results indicate that although Sierra Nevada bighorn sheep seek out foraging opportunities related to landscape phenology, they compensate for short-term environmental stressors by undertaking brief up- and downslope vertical movements. Migrants may therefore be impacted by future warming and increased storm frequency or intensity, with shifts in annual migration timing, and fine-scale vertical movement responses to environmental variability., (© 2024. The Author(s).)

We construct non-linear machine learning (ML) prediction models for systolic and diastolic blood pressure (SBP, DBP) using demographic and clinical variables and polygenic risk scores (PRSs). We developed a two-model ensemble, consisting of a baseline model, where prediction is based on demographic and clinical variables only, and a genetic model, where we also include PRSs. We evaluate the use of a linear versus a non-linear model at both the baseline and the genetic model levels and assess the improvement in performance when incorporating multiple PRSs. We report the ensemble model's performance as percentage variance explained (PVE) on a held-out test dataset. A non-linear baseline model improved the PVEs from 28.1 to 30.1% (SBP) and 14.3% to 17.4% (DBP) compared with a linear baseline model. Including seven PRSs in the genetic model computed based on the largest available GWAS of SBP/DBP improved the genetic model PVE from 4.8 to 5.1% (SBP) and 4.7 to 5% (DBP) compared to using a single PRS. Adding additional 14 PRSs computed based on two independent GWASs further increased the genetic model PVE to 6.3% (SBP) and 5.7% (DBP). PVE differed across self-reported race/ethnicity groups, with primarily all non-White groups benefitting from the inclusion of additional PRSs. In summary, non-linear ML models improves BP prediction in models incorporating diverse populations., (© 2024. The Author(s).)

Type 2 diabetes (T2D) is a heterogeneous disease that develops through diverse pathophysiological processes 1,2 and molecular mechanisms that are often specific to cell type 3,4 . Here, to characterize the genetic contribution to these processes across ancestry groups, we aggregate genome-wide association study data from 2,535,601 individuals (39.7% not of European ancestry), including 428,452 cases of T2D. We identify 1,289 independent association signals at genome-wide significance (P < 5 × 10 -8 ) that map to 611 loci, of which 145 loci are, to our knowledge, previously unreported. We define eight non-overlapping clusters of T2D signals that are characterized by distinct profiles of cardiometabolic trait associations. These clusters are differentially enriched for cell-type-specific regions of open chromatin, including pancreatic islets, adipocytes, endothelial cells and enteroendocrine cells. We build cluster-specific partitioned polygenic scores 5 in a further 279,552 individuals of diverse ancestry, including 30,288 cases of T2D, and test their association with T2D-related vascular outcomes. Cluster-specific partitioned polygenic scores are associated with coronary artery disease, peripheral artery disease and end-stage diabetic nephropathy across ancestry groups, highlighting the importance of obesity-related processes in the development of vascular outcomes. Our findings show the value of integrating multi-ancestry genome-wide association study data with single-cell epigenomics to disentangle the aetiological heterogeneity that drives the development and progression of T2D. This might offer a route to optimize global access to genetically informed diabetes care., (© 2024. The Author(s).)

Background and Objective: Persistent prostatic specific antigen (PSA) represents a poor prognostic factor for recurrence after radical prostatectomy (RP). However, the impact of persistent PSA on oncologic outcomes in patients undergoing salvage RP is unknown. To investigate the impact of persistent PSA after salvage RP on long-term oncologic outcomes., Material and Methods: Patients who underwent salvage RP for recurrent prostate cancer between 2000 and 2021 were identified from twelve high-volume centers. Only patients with available PSA after salvage RP were included. Kaplan-Meier analyses and multivariable Cox regression models were used to test the effect of persistent PSA on biochemical recurrence (BCR), metastasis and any death after salvage RP. Persistent PSA was defined as a PSA-value ≥ 0.1 ng/ml, at first PSA-measurement after salvage RP., Results: Overall, 580 patients were identified. Of those, 42% (n = 242) harbored persistent PSA. Median follow-up after salvage RP was 38 months, median time to salvage RP was 64 months and median time to first PSA after salvage RP was 2.2 months. At 84 months after salvage RP, BCR-free, metastasis-free, and overall survival was 6.6 vs. 59%, 71 vs. 88% and 77 vs. 94% for patients with persistent vs. undetectable PSA after salvage RP (all p < 0.01). In multivariable Cox models persistent PSA was an independent predictor for BCR (HR: 5.47, p < 0.001) and death (HR: 3.07, p < 0.01)., Conclusion: Persistent PSA is common after salvage RP and represents an independent predictor for worse oncologic outcomes. Patients undergoing salvage RP should be closely monitored after surgery to identify those with persistent PSA., (© 2023. The Author(s).)

Thwaites Glacier is one of the fastest-changing ice-ocean systems in Antarctica 1-3 . Much of the ice sheet within the catchment of Thwaites Glacier is grounded below sea level on bedrock that deepens inland 4 , making it susceptible to rapid and irreversible ice loss that could raise the global sea level by more than half a metre 2,3,5 . The rate and extent of ice loss, and whether it proceeds irreversibly, are set by the ocean conditions and basal melting within the grounding-zone region where Thwaites Glacier first goes afloat 3,6 , both of which are largely unknown. Here we show-using observations from a hot-water-drilled access hole-that the grounding zone of Thwaites Eastern Ice Shelf (TEIS) is characterized by a warm and highly stable water column with temperatures substantially higher than the in situ freezing point. Despite these warm conditions, low current speeds and strong density stratification in the ice-ocean boundary layer actively restrict the vertical mixing of heat towards the ice base 7,8 , resulting in strongly suppressed basal melting. Our results demonstrate that the canonical model of ice-shelf basal melting used to generate sea-level projections cannot reproduce observed melt rates beneath this critically important glacier, and that rapid and possibly unstable grounding-line retreat may be associated with relatively modest basal melt rates., (© 2023. The Author(s).)

Thwaites Glacier represents 15% of the ice discharge from the West Antarctic Ice Sheet and influences a wider catchment 1-3 . Because it is grounded below sea level 4,5 , Thwaites Glacier is thought to be susceptible to runaway retreat triggered at the grounding line (GL) at which the glacier reaches the ocean 6,7 . Recent ice-flow acceleration 2,8 and retreat of the ice front 8-10 and GL 11,12 indicate that ice loss will continue. The relative impacts of mechanisms underlying recent retreat are however uncertain. Here we show sustained GL retreat from at least 2011 to 2020 and resolve mechanisms of ice-shelf melt at the submetre scale. Our conclusions are based on observations of the Thwaites Eastern Ice Shelf (TEIS) from an underwater vehicle, extending from the GL to 3 km oceanward and from the ice-ocean interface to the sea floor. These observations show a rough ice base above a sea floor sloping upward towards the GL and an ocean cavity in which the warmest water exceeds 2 °C above freezing. Data closest to the ice base show that enhanced melting occurs along sloped surfaces that initiate near the GL and evolve into steep-sided terraces. This pronounced melting along steep ice faces, including in crevasses, produces stratification that suppresses melt along flat interfaces. These data imply that slope-dependent melting sculpts the ice base and acts as an important response to ocean warming., (© 2023. The Author(s), under exclusive licence to Springer Nature Limited.)

The Harmonized Diagnostic Assessment of Dementia for the Longitudinal Aging Study in India (LASI-DAD) is a nationally representative in-depth study of cognitive aging and dementia. We present a publicly available dataset of harmonized cognitive measures of 4,096 adults 60 years of age and older in India, collected across 18 states and union territories. Blood samples were obtained to carry out whole blood and serum-based assays. Results are included in a venous blood specimen datafile that can be linked to the Harmonized LASI-DAD dataset. A global screening array of 960 LASI-DAD respondents is also publicly available for download, in addition to neuroimaging data on 137 LASI-DAD participants. Altogether, these datasets provide comprehensive information on older adults in India that allow researchers to further understand risk factors associated with cognitive impairment and dementia., (© 2023. The Author(s).)

Tobacco and alcohol use are heritable behaviours associated with 15% and 5.3% of worldwide deaths, respectively, due largely to broad increased risk for disease and injury 1-4 . These substances are used across the globe, yet genome-wide association studies have focused largely on individuals of European ancestries 5 . Here we leveraged global genetic diversity across 3.4 million individuals from four major clines of global ancestry (approximately 21% non-European) to power the discovery and fine-mapping of genomic loci associated with tobacco and alcohol use, to inform function of these loci via ancestry-aware transcriptome-wide association studies, and to evaluate the genetic architecture and predictive power of polygenic risk within and across populations. We found that increases in sample size and genetic diversity improved locus identification and fine-mapping resolution, and that a large majority of the 3,823 associated variants (from 2,143 loci) showed consistent effect sizes across ancestry dimensions. However, polygenic risk scores developed in one ancestry performed poorly in others, highlighting the continued need to increase sample sizes of diverse ancestries to realize any potential benefit of polygenic prediction., (© 2022. The Author(s).)

Common single-nucleotide polymorphisms (SNPs) are predicted to collectively explain 40-50% of phenotypic variation in human height, but identifying the specific variants and associated regions requires huge sample sizes 1 . Here, using data from a genome-wide association study of 5.4 million individuals of diverse ancestries, we show that 12,111 independent SNPs that are significantly associated with height account for nearly all of the common SNP-based heritability. These SNPs are clustered within 7,209 non-overlapping genomic segments with a mean size of around 90 kb, covering about 21% of the genome. The density of independent associations varies across the genome and the regions of increased density are enriched for biologically relevant genes. In out-of-sample estimation and prediction, the 12,111 SNPs (or all SNPs in the HapMap 3 panel 2 ) account for 40% (45%) of phenotypic variance in populations of European ancestry but only around 10-20% (14-24%) in populations of other ancestries. Effect sizes, associated regions and gene prioritization are similar across ancestries, indicating that reduced prediction accuracy is likely to be explained by linkage disequilibrium and differences in allele frequency within associated regions. Finally, we show that the relevant biological pathways are detectable with smaller sample sizes than are needed to implicate causal genes and variants. Overall, this study provides a comprehensive map of specific genomic regions that contain the vast majority of common height-associated variants. Although this map is saturated for populations of European ancestry, further research is needed to achieve equivalent saturation in other ancestries., (© 2022. The Author(s).)

The molecular evolution of cytochromes P450 and associated redox-driven oxidative catalysis remains a mystery in biology. It is widely believed that sterol 14α-demethylase (CYP51), an essential enzyme of sterol biosynthesis, is the ancestor of the whole P450 superfamily given its conservation across species in different biological kingdoms. Herein we have utilized X-ray crystallography, molecular dynamics simulations, phylogenetics and electron transfer measurements to interrogate the nature of P450-redox partner binding using the naturally occurring fusion protein, CYP51-ferredoxin found in the sterol-producing bacterium Methylococcus capsulatus. Our data advocates that the electron transfer mechanics in the M. capsulatus CYP51-ferredoxin fusion protein involves an ensemble of ferredoxin molecules in various orientations and the interactions are transient. Close proximity of ferredoxin, however, is required to complete the substrate-induced large-scale structural switch in the P450 domain that enables proton-coupled electron transfer and subsequent oxygen scission and catalysis. These results have fundamental implications regarding the early evolution of electron transfer proteins and for the redox reactions in the early steps of sterol biosynthesis. They also shed new light on redox protein mechanics and the subsequent diversification of the P450 electron transfer machinery in nature., (© 2022. The Author(s).)

Chronic cough is globally prevalent across all age groups. This disorder is challenging to treat because many pulmonary and extrapulmonary conditions can present with chronic cough, and cough can also be present without any identifiable underlying cause or be refractory to therapies that improve associated conditions. Most patients with chronic cough have cough hypersensitivity, which is characterized by increased neural responsivity to a range of stimuli that affect the airways and lungs, and other tissues innervated by common nerve supplies. Cough hypersensitivity presents as excessive coughing often in response to relatively innocuous stimuli, causing significant psychophysical morbidity and affecting patients' quality of life. Understanding of the mechanisms that contribute to cough hypersensitivity and excessive coughing in different patient populations and across the lifespan is advancing and has contributed to the development of new therapies for chronic cough in adults. Owing to differences in the pathology, the organs involved and individual patient factors, treatment of chronic cough is progressing towards a personalized approach, and, in the future, novel ways to endotype patients with cough may prove valuable in management., (© 2022. Springer Nature Limited.)

The oncogenic potential of the latent transcription factor signal transducer and activator of transcription (STAT)3 in many human cancers, including lung cancer, has been largely attributed to its nuclear activity as a tyrosine-phosphorylated (pY 705 site) transcription factor. By contrast, an alternate mitochondrial pool of serine phosphorylated (pS 727 site) STAT3 has been shown to promote tumourigenesis by regulating metabolic processes, although this has been reported in only a restricted number of mutant RAS-addicted neoplasms. Therefore, the involvement of STAT3 serine phosphorylation in the pathogenesis of most cancer types, including mutant KRAS lung adenocarcinoma (LAC), is unknown. Here, we demonstrate that LAC is suppressed in oncogenic Kras G12D -driven mouse models engineered for pS 727 -STAT3 deficiency. The proliferative potential of the transformed Kras G12D lung epithelium, and mutant KRAS human LAC cells, was significantly reduced upon pS 727 -STAT3 deficiency. Notably, we uncover the multifaceted capacity of constitutive pS 727 -STAT3 to metabolically reprogramme LAC cells towards a hyper-proliferative state by regulating nuclear and mitochondrial (mt) gene transcription, the latter via the mtDNA transcription factor, TFAM. Collectively, our findings reveal an obligate requirement for the transcriptional activity of pS 727 -STAT3 in mutant KRAS-driven LAC with potential to guide future therapeutic targeting approaches., (© 2021. The Author(s), under exclusive licence to Springer Nature Limited.)

Increased blood lipid levels are heritable risk factors of cardiovascular disease with varied prevalence worldwide owing to different dietary patterns and medication use 1 . Despite advances in prevention and treatment, in particular through reducing low-density lipoprotein cholesterol levels 2 , heart disease remains the leading cause of death worldwide 3 . Genome-wideassociation studies (GWAS) of blood lipid levels have led to important biological and clinical insights, as well as new drug targets, for cardiovascular disease. However, most previous GWAS 4-23 have been conducted in European ancestry populations and may have missed genetic variants that contribute to lipid-level variation in other ancestry groups. These include differences in allele frequencies, effect sizes and linkage-disequilibrium patterns 24 . Here we conduct a multi-ancestry, genome-wide genetic discovery meta-analysis of lipid levels in approximately 1.65 million individuals, including 350,000 of non-European ancestries. We quantify the gain in studying non-European ancestries and provide evidence to support the expansion of recruitment of additional ancestries, even with relatively small sample sizes. We find that increasing diversity rather than studying additional individuals of European ancestry results in substantial improvements in fine-mapping functional variants and portability of polygenic prediction (evaluated in approximately 295,000 individuals from 7 ancestry groupings). Modest gains in the number of discovered loci and ancestry-specific variants were also achieved. As GWAS expand emphasis beyond the identification of genes and fundamental biology towards the use of genetic variants for preventive and precision medicine 25 , we anticipate that increased diversity of participants will lead to more accurate and equitable 26 application of polygenic scores in clinical practice., (© 2021. The Author(s), under exclusive licence to Springer Nature Limited.)

Understanding the population structure and movements of the invasive fall armyworm (FAW, Spodoptera frugiperda) is important as it can help mitigate crop damage, and highlight areas at risk of outbreaks or evolving insecticide resistance. Determining population structure in invasive FAW has been a challenge due to genetic mutations affecting the markers traditionally used for strain and haplotype identification; mitochondrial cytochrome oxidase I (COIB) and the Z-chromosome-linked Triosephosphate isomerase (Tpi). Here, we compare the results from COIB and Tpi markers with highly variable repeat regions (microsatellites) to improve our understanding of FAW population structure in Africa. There was very limited genetic diversity using the COIB marker, whereas using the TpiI4 marker there was greater diversity that showed very little evidence of genetic structuring between FAW populations across Africa. There was greater genetic diversity identified using microsatellites, and this revealed a largely panmictic population of FAW alongside some evidence of genetic structuring between countries. It is hypothesised here that FAW are using long-distance flight and prevailing winds to frequently move throughout Africa leading to population mixing. These approaches combined provide important evidence that genetic mixing between invasive FAW populations may be more common than previously reported., (© 2021. The Author(s).)

Reproductive longevity is essential for fertility and influences healthy ageing in women 1,2 , but insights into its underlying biological mechanisms and treatments to preserve it are limited. Here we identify 290 genetic determinants of ovarian ageing, assessed using normal variation in age at natural menopause (ANM) in about 200,000 women of European ancestry. These common alleles were associated with clinical extremes of ANM; women in the top 1% of genetic susceptibility have an equivalent risk of premature ovarian insufficiency to those carrying monogenic FMR1 premutations 3 . The identified loci implicate a broad range of DNA damage response (DDR) processes and include loss-of-function variants in key DDR-associated genes. Integration with experimental models demonstrates that these DDR processes act across the life-course to shape the ovarian reserve and its rate of depletion. Furthermore, we demonstrate that experimental manipulation of DDR pathways highlighted by human genetics increases fertility and extends reproductive life in mice. Causal inference analyses using the identified genetic variants indicate that extending reproductive life in women improves bone health and reduces risk of type 2 diabetes, but increases the risk of hormone-sensitive cancers. These findings provide insight into the mechanisms that govern ovarian ageing, when they act, and how they might be targeted by therapeutic approaches to extend fertility and prevent disease., (© 2021. The Author(s), under exclusive licence to Springer Nature Limited.)

The Trans-Omics for Precision Medicine (TOPMed) programme seeks to elucidate the genetic architecture and biology of heart, lung, blood and sleep disorders, with the ultimate goal of improving diagnosis, treatment and prevention of these diseases. The initial phases of the programme focused on whole-genome sequencing of individuals with rich phenotypic data and diverse backgrounds. Here we describe the TOPMed goals and design as well as the available resources and early insights obtained from the sequence data. The resources include a variant browser, a genotype imputation server, and genomic and phenotypic data that are available through dbGaP (Database of Genotypes and Phenotypes) 1 . In the first 53,831 TOPMed samples, we detected more than 400 million single-nucleotide and insertion or deletion variants after alignment with the reference genome. Additional previously undescribed variants were detected through assembly of unmapped reads and customized analysis in highly variable loci. Among the more than 400 million detected variants, 97% have frequencies of less than 1% and 46% are singletons that are present in only one individual (53% among unrelated individuals). These rare variants provide insights into mutational processes and recent human evolutionary history. The extensive catalogue of genetic variation in TOPMed studies provides unique opportunities for exploring the contributions of rare and noncoding sequence variants to phenotypic variation. Furthermore, combining TOPMed haplotypes with modern imputation methods improves the power and reach of genome-wide association studies to include variants down to a frequency of approximately 0.01%.

Age is the dominant risk factor for most chronic human diseases, but the mechanisms through which ageing confers this risk are largely unknown 1 . The age-related acquisition of somatic mutations that lead to clonal expansion in regenerating haematopoietic stem cell populations has recently been associated with both haematological cancer 2-4 and coronary heart disease 5 -this phenomenon is termed clonal haematopoiesis of indeterminate potential (CHIP) 6 . Simultaneous analyses of germline and somatic whole-genome sequences provide the opportunity to identify root causes of CHIP. Here we analyse high-coverage whole-genome sequences from 97,691 participants of diverse ancestries in the National Heart, Lung, and Blood Institute Trans-omics for Precision Medicine (TOPMed) programme, and identify 4,229 individuals with CHIP. We identify associations with blood cell, lipid and inflammatory traits that are specific to different CHIP driver genes. Association of a genome-wide set of germline genetic variants enabled the identification of three genetic loci associated with CHIP status, including one locus at TET2 that was specific to individuals of African ancestry. In silico-informed in vitro evaluation of the TET2 germline locus enabled the identification of a causal variant that disrupts a TET2 distal enhancer, resulting in increased self-renewal of haematopoietic stem cells. Overall, we observe that germline genetic variation shapes haematopoietic stem cell function, leading to CHIP through mechanisms that are specific to clonal haematopoiesis as well as shared mechanisms that lead to somatic mutations across tissues.

Next-generation sequencing (NGS) is widely used in genetic testing for the highly sensitive detection of single nucleotide changes and small insertions or deletions. However, detection and phasing of structural variants, especially in repetitive or homologous regions, can be problematic due to uneven read coverage or genome reference bias, resulting in false calls. To circumvent this challenge, a computational approach utilizing customized scaffolds as supplementary reference sequences for read alignment was developed, and its effectiveness demonstrated with two CBS gene variants: NM&#95;000071.2:c.833T>C and NM&#95;000071.2:c.[833T>C; 844&#95;845ins68]. Variant c.833T>C is a known causative mutation for homocystinuria, but is not pathogenic when in cis with the insertion, c.844&#95;845ins68, because of alternative splicing. Using simulated reads, the custom scaffolds method resolved all possible combinations with 100% accuracy and, based on > 60,000 clinical specimens, exceeded the performance of current approaches that only align reads to GRCh37/hg19 for the detection of c.833T>C alone or in cis with c.844&#95;845ins68. Furthermore, analysis of two 1000 Genomes Project trios revealed that the c.[833T>C; 844&#95;845ins68] complex variant had previously been undetected in these datasets, likely due to the alignment method used. This approach can be configured for existing workflows to detect other challenging and potentially underrepresented variants, thereby augmenting accurate variant calling in clinical NGS testing.

The hydrogen isotopes deuterium (D) and tritium (T) have become essential tools in chemistry, biology and medicine 1 . Beyond their widespread use in spectroscopy, mass spectrometry and mechanistic and pharmacokinetic studies, there has been considerable interest in incorporating deuterium into drug molecules 1 . Deutetrabenazine, a deuterated drug that is promising for the treatment of Huntington's disease 2 , was recently approved by the United States' Food and Drug Administration. The deuterium kinetic isotope effect, which compares the rate of a chemical reaction for a compound with that for its deuterated counterpart, can be substantial 1,3,4 . The strategic replacement of hydrogen with deuterium can affect both the rate of metabolism and the distribution of metabolites for a compound 5 , improving the efficacy and safety of a drug. The pharmacokinetics of a deuterated compound depends on the location(s) of deuterium. Although methods are available for deuterium incorporation at both early and late stages of the synthesis of a drug 6,7 , these processes are often unselective and the stereoisotopic purity can be difficult to measure 7,8 . Here we describe the preparation of stereoselectively deuterated building blocks for pharmaceutical research. As a proof of concept, we demonstrate a four-step conversion of benzene to cyclohexene with varying degrees of deuterium incorporation, via binding to a tungsten complex. Using different combinations of deuterated and proteated acid and hydride reagents, the deuterated positions on the cyclohexene ring can be controlled precisely. In total, 52 unique stereoisotopomers of cyclohexene are available, in the form of ten different isotopologues. This concept can be extended to prepare discrete stereoisotopomers of functionalized cyclohexenes. Such systematic methods for the preparation of pharmacologically active compounds as discrete stereoisotopomers could improve the pharmacological and toxicological properties of drugs and provide mechanistic information related to their distribution and metabolism in the body.

We use ultrasonic spray-coating to sequentially deposit thin films of tin oxide, a triple-cation perovskite and spiro-OMeTAD, allowing us fabricate perovskite solar cells (PSCs) with a champion reverse scan power conversion efficiency (PCE) of 19.4% on small-area substrates. We show that the use of spray-deposition permits us to rapidly (>80 mm s -1 ) coat 25 mm × 75 mm substrates that were divided into a series of devices each with an active area of 15.4 mm 2 , yielding an average PCE of 10.3% and a peak PCE of 16.3%. By connecting seven 15.4 mm 2 devices in parallel on a single substrate, we create a device having an effective active area of 1.08 cm 2 and a PCE of 12.7%. This work demonstrates the possibility for spray-coating to fabricate high efficiency and low-cost perovskite solar cells at speed.

The mid-Cretaceous period was one of the warmest intervals of the past 140 million years 1-5 , driven by atmospheric carbon dioxide levels of around 1,000 parts per million by volume 6 . In the near absence of proximal geological records from south of the Antarctic Circle, it is disputed whether polar ice could exist under such environmental conditions. Here we use a sedimentary sequence recovered from the West Antarctic shelf-the southernmost Cretaceous record reported so far-and show that a temperate lowland rainforest environment existed at a palaeolatitude of about 82° S during the Turonian-Santonian age (92 to 83 million years ago). This record contains an intact 3-metre-long network of in situ fossil roots embedded in a mudstone matrix containing diverse pollen and spores. A climate model simulation shows that the reconstructed temperate climate at this high latitude requires a combination of both atmospheric carbon dioxide concentrations of 1,120-1,680 parts per million by volume and a vegetated land surface without major Antarctic glaciation, highlighting the important cooling effect exerted by ice albedo under high levels of atmospheric carbon dioxide.

Elevated concentrations of fluoride ions (F - ) in natural groundwater are a worldwide problem. Discarded oyster shells were ground to ≤100 µm particle size to produce oyster shell powder (OS). A subset of the OS was heated to produce calcined oyster shell (COS). A subset of the COS was further treated with 1 M phosphoric acid to produce phosphoric-acid-treated oyster shell (POS). OS and COS were combined with phosphoric acid (1.6 mM and 3.2 mM) to produce OS + P (oyster shell with phosphoric acid) and COS + P (calcined oyster shell with phosphoric acid). OS and COS removed 46% and 50% (10 g/L of sorbent dose) but POS, OS + P and COS + P removed 96%, 100% and 76% (1 g/L of sorbent dose) when the initial concentration of fluoride was 10 mg/L. The sorption kinetics of POS, OS + P and COS + P followed second-order reaction rates, and sorption isotherms of all sorbents were well-described by the Freundlich sorption isotherm. These results indicate that oyster shells can be an effective sorbent for fluoride removal, with the added benefit of re-use of a waste product.

The Antarctic Peninsula Ice Sheet is currently experiencing sustained and accelerating loss of ice. Determining when these changes were initiated and identifying the main drivers is hampered by the short instrumental record (1992 to present). Here we present a 6,250 year record of glacial discharge based on the oxygen isotope composition of diatoms (δ 18 O diatom ) from a marine core located at the north-eastern tip of the Antarctic Peninsula. We find that glacial discharge - sourced primarily from ice shelf and iceberg melting along the eastern Antarctic Peninsula - remained largely stable between ~6,250 to 1,620 cal. yr BP, with a slight increase in variability until ~720 cal. yr. BP. An increasing trend in glacial discharge occurs after 550 cal. yr BP (A.D. 1400), reaching levels unprecedented during the past 6,250 years after 244 cal. yr BP (A.D. 1706). A marked acceleration in the rate of glacial discharge is also observed in the early part of twentieth century (after A.D. 1912). Enhanced glacial discharge, particularly after the 1700s is linked to a positive Southern Annular Mode (SAM). We argue that a positive SAM drove stronger westerly winds, atmospheric warming and surface ablation on the eastern Antarctic Peninsula whilst simultaneously entraining more warm water into the Weddell Gyre, potentially increasing melting on the undersides of ice shelves. A possible implication of our data is that ice shelves in this region have been thinning for at least ~300 years, potentially predisposing them to collapse under intensified anthropogenic warming.

Domestic dogs are the most abundant large carnivore on the planet, and their ubiquity has led to concern regarding the impacts of dogs as predators of and competitors with native wildlife. If native large carnivores perceive dogs as threatening, impacts could extend to the community level by altering interactions between large carnivores and their prey. Dog impacts may be further exacerbated if these human-associated predators are also perceived as indicators of risk from humans. However, observational approaches used to date have led to ambiguity regarding the effects of dog presence on wildlife. We experimentally quantified dog impacts on the behavior of a native large carnivore, presenting playbacks of dog vocalizations to pumas in central California. We show that the perceived presence of dogs has minimal impacts on puma behavior at their kill sites, and is no more likely to affect total feeding time at kills than non-threatening controls. We previously demonstrated that pumas exhibit strong responses to human cues, and here show that perceived risk from human presence far exceeds that from dogs. Our results suggest that protected areas management policies that restrict dogs but permit human access may in some cases be of limited value for large carnivores.

Background: The first gene therapy for an inherited retinal dystrophy recently received market approval in the United States; multiple other gene therapies are in the clinical pipeline. Thus far, gene therapy has commanded prices in the range of $500,000 to over $1,000,000 for the one-time doses and have been indicated for highly orphan diseases where there is no other viable treatment option. To be adopted by healthcare systems, gene therapy will need to show clinical benefit in line with its increased costs. Before longitudinal patient studies are available, model-based estimations will be necessary to project the full clinical benefit of gene therapy., Methods: To investigate the lifetime benefit of gene therapy for the retinal dystrophy choroideremia, we have built a Markov model of disease progression informed by clinical data of AAV.REP1 and voretigene neparvovec (Luxturna, Spark Therapeutics). Gene therapy patient benefit was estimated by quality-adjusted life years (QALYs) in three hypothetical disease severity patient groups. The severity of disease was defined by the combined effect of remaining retinal area and visual acuity and assigned corresponding health utility values., Results: Early-stage patients treated with gene therapy were estimated to gain, in average, 14.30 QALYs over standard-of-care, mid-stage patients 6.22 QALYs, and late-stage patients 1.48 QALYs over untreated patients during their lifetime owing to treatment. Cost-effectiveness was not assessed as AAV.REP1 is still in clinical trials., Conclusions: In young adults in the earlier stages of choroideremia, successful gene therapy is expected to provide a significant increase in health-related quality of life.

'Black silicon' (bSi) samples with surfaces covered in nanoneedles of length ~5 µm were fabricated using a plasma etching process and then coated with a conformal uniform layer of diamond using hot filament chemical vapour deposition to produce 'black diamond' (bD) nanostructures. The diamond needles were then chemically terminated with H, O, NH 2 or F using plasma treatment, and the hydrophilicity of the resulting surfaces were assessed using water droplet contact-angle measurements, and scaled in the order O > H ≈NH 2 >F, with the F-terminated surface being superhydrophobic. The effectiveness of these differently terminated bD needles in killing the Gram-negative bacterium E. coli was semi-quantified by Live/Dead staining and fluorescence microscopy, and visualised by environmental scanning electron microscopy. The total number of adhered bacteria was consistent for all the nanostructured bD surfaces at around 50% of the value for the flat diamond control. This, combined with a chemical bactericidal effect of 20-30%, shows that the nanostructured bD surfaces supported significantly fewer viable E. coli than flat surfaces. Moreover, the bD surfaces were particularly effective at preventing the establishment of bacterial aggregates - a precursor to biofilm formation. The percentage of dead bacteria also decreased as a function of hydrophilicity. These results are consistent with a predominantly mechanical mechanism for bacteria death based on the stretching and disruption of the cell membrane, combined with an additional effect from the chemical nature of the surface.

High blood pressure (BP) is a major risk factor for cardiovascular disease (CVD) and is more prevalent in African Americans as compared to other US groups. Although large, population-based genome-wide association studies (GWAS) have identified over 300 common polymorphisms modulating inter-individual BP variation, largely in European ancestry subjects, most of them do not localize to regions previously identified through family-based linkage studies. This discrepancy has remained unexplained despite the statistical power differences between current GWAS and prior linkage studies. To address this issue, we performed genome-wide linkage analysis of BP traits in African-American families from the Family Blood Pressure Program (FBPP) and genotyped on the Illumina Human Exome BeadChip v1.1. We identified a genomic region on chromosome 1q31 with LOD score 3.8 for pulse pressure (PP), a region we previously implicated in DBP studies of European ancestry families. Although no reported GWAS variants map to this region, combined linkage and association analysis of PP identified 81 rare and low frequency exonic variants accounting for the linkage evidence. Replication analysis in eight independent African ancestry cohorts (N = 16,968) supports this specific association with PP (P = 0.0509). Additional association and network analyses identified multiple potential candidate genes in this region expressed in multiple tissues and with a strong biological support for a role in BP. In conclusion, multiple genes and rare variants on 1q31 contribute to PP variation. Beyond producing new insights into PP, we demonstrate how family-based linkage and association studies can implicate specific rare and low frequency variants for complex traits.

Objective: To examine low birthweight and preterm birth of second children born to home-visited first-time mothers., Subjects: Women were previously recruited for a randomized controlled trial of the home visiting model disseminated as Nurse-Family Partnership. 512 of these women had second children within 18 years of the first child's birth, and were included in our sample., Results: The intervention was associated with a lower likelihood of low birthweight for second children (odds ratio: 0.51, 95% CI: 0.27, 0.97), an effect apparent only if the first-born had low birthweight and mediated by close birth spacing. These moderation and mediation patterns were similar in the preterm birth outcome., Conclusion: A home visiting program provided for first-born children reduced low birthweight for second-born children, if the first-born had low birthweight. This finding implies a broader impact than previously documented, because few studies have included these second children.

Category 4 landfalling hurricane Harvey poured more than a metre of rainfall across the heavily populated Houston area, leading to unprecedented flooding and damage. Although studies have focused on the contribution of anthropogenic climate change to this extreme rainfall event 1-3 , limited attention has been paid to the potential effects of urbanization on the hydrometeorology associated with hurricane Harvey. Here we find that urbanization exacerbated not only the flood response but also the storm total rainfall. Using the Weather Research and Forecast model-a numerical model for simulating weather and climate at regional scales-and statistical models, we quantify the contribution of urbanization to rainfall and flooding. Overall, we find that the probability of such extreme flood events across the studied basins increased on average by about 21 times in the period 25-30 August 2017 because of urbanization. The effect of urbanization on storm-induced extreme precipitation and flooding should be more explicitly included in global climate models, and this study highlights its importance when assessing the future risk of such extreme events in highly urbanized coastal areas.

Larval fishes are a useful metric of marine ecosystem state and change, as well as species-specific patterns in phenology. The high level of taxonomic expertise required to identify larval fishes to species level, and the considerable effort required to collect samples, make these data very valuable. Here we collate 3178 samples of larval fish assemblages, from 12 research projects from 1983-present, from temperate and subtropical Australian pelagic waters. This forms a benchmark for the larval fish assemblage for the region, and includes recent monitoring of larval fishes at coastal oceanographic reference stations. Comparing larval fishes among projects can be problematic due to differences in taxonomic resolution, and identifying all taxa to species is challenging, so this study reports a standard taxonomic resolution (of 218 taxa) for this region to help guide future research. This larval fish database serves as a data repository for surveys of larval fish assemblages in the region, and can contribute to analysis of climate-driven changes in the location and timing of the spawning of marine fishes.

Consumption rates are the foundation of trophic ecology, yet bioenergetics models used to estimate these rates can lack realism by not incorporating the ontogeny of diet. We constructed a bioenergetics model of a marine predatory fish (tailor, Pomatomus saltatrix) that incorporated high-resolution ontogenetic diet variation, and compared consumption estimates to those derived from typical bioenergetics models that do not consider ontogenetic diet variation. We found tailor consumption was over- or under-estimated by ~5-25% when only including the most common prey item. This error was due to a positive relationship between mean prey energy density and predator body size. Since high-resolution diet data isn't always available, we also simulated how increasing dietary information progressively influenced consumption rate estimates. The greatest improvement in consumption rate estimates occurred when diet variation of 2-3 stanzas (1-2 juvenile stanzas, and adults) was included, with at least 5-6 most common prey types per stanza. We recommend increased emphasis on incorporating the ontogeny of diet and prey energy density in consumption rate estimates, especially for species with spatially segregated life stages or variable diets. A small-moderate increase in the resolution of dietary information can greatly benefit the accuracy of estimated consumption rates. We present a method of incorporating variable prey energy density into bioenergetics models.

This corrects the article DOI: 10.1038/nature20136.

Glaciological and oceanographic observations coupled with numerical models show that warm Circumpolar Deep Water (CDW) incursions onto the West Antarctic continental shelf cause melting of the undersides of floating ice shelves. Because these ice shelves buttress glaciers feeding into them, their ocean-induced thinning is driving Antarctic ice-sheet retreat today. Here we present a multi-proxy data based reconstruction of variability in CDW inflow to the Amundsen Sea sector, the most vulnerable part of the West Antarctic Ice Sheet, during the Holocene epoch (from 11.7 thousand years ago to the present). The chemical compositions of foraminifer shells and benthic foraminifer assemblages in marine sediments indicate that enhanced CDW upwelling, controlled by the latitudinal position of the Southern Hemisphere westerly winds, forced deglaciation of this sector from at least 10,400 years ago until 7,500 years ago-when an ice-shelf collapse may have caused rapid ice-sheet thinning further upstream-and since the 1940s. These results increase confidence in the predictive capability of current ice-sheet models.

Accumulating evidence indicates that Toll-like receptor (TLR) signaling adapter protein interactions with Toll/Interleukin-1 Receptor (TIR) domains present in sensory neurons may modulate neuropathic pain states. Following ligand interaction with TLRs, TIR serves to both initiate intracellular signaling and facilitate recruitment of signaling adapter proteins to the intracytoplasmic domain. Although TLR TIR is central to a number of TLR signaling cascades, its role in sensory neurons is poorly understood. In this study we investigated the degree to which TLR TIR decoy peptide modified to include a TAT sequence (Trans-Activator of Transcription gene in HIV; TAT-4BB) affected LPS-induced intracellular calcium flux and excitation in sensory neurons, and behavioral changes due to TLR4 active metabolite, morphine-3-glucuronide (M3G) exposure in vivo. TAT-4BB inhibited LPS-induced calcium changes in a majority of sensory neurons and decreased LPS-dependent neuronal excitability in small diameter neurons. Acute systemic administration of the TAT-4BB reversed M3G-induced tactile allodynia in a dose-dependent manner but did not affect motor activity, anxiety or responses to noxious thermal stimulus. These data suggest that targeting TLR TIR domains may provide novel pharmacological targets to reduce or reverse TLR4-dependent pain behavior in the rodent.

Interferon regulatory factor 5 (IRF5) is a key transcription factor of innate immunity, which plays an important role in host restriction to viral infection and inflammation. Genome-wide association studies have implied the association of IRF5 with several autoimmune diseases, including systemic lupus erythematosus (SLE), Sjogren's syndrome, inflammatory bowel disease and multiple sclerosis. However, the regulation of IRF5-mediated immunity is not well understood. To uncover new regulators in IRF5 pathway, we used two "omics" approaches: affinity purification coupled with mass spectrometry and a high throughput RNAi screen. Proteomics identified 16 new IRF5 interactors while RNAi-mediated knockdown found 43 regulators of the TLR7-dependent IRF5 signaling pathway. NXF1 was identified in both screens. Stimulation with TLR7 ligand enhances formation of IRF5-NXF1 protein complexes. Gain or loss-of-function experiments revealed NXF1 selectively regulates TLR7-driven IRF5 transcriptional activity, suggesting a new role for NXF1 in the IRF5 signaling pathway.

This week, Scientific Data published a collection of eight papers that describe datasets from high-throughput functional genomics screens, primarily utilizing RNA interference (RNAi). The publications explore host-pathogen dependencies, innate immune response, disease pathways, and cell morphology and motility at the genome-level. All data, including raw images from the high content screens, are publically available in PubChem BioAssay, figshare, Harvard Dataverse or the Image Data Resource (IDR). Detailed data descriptors enable use of these data for analysis algorithm design, machine learning, data comparisons, as well as generating new scientific hypotheses.

Height is a highly heritable, classic polygenic trait with approximately 700 common associated variants identified through genome-wide association studies so far. Here, we report 83 height-associated coding variants with lower minor-allele frequencies (in the range of 0.1-4.8%) and effects of up to 2 centimetres per allele (such as those in IHH, STC2, AR and CRISPLD2), greater than ten times the average effect of common variants. In functional follow-up studies, rare height-increasing alleles of STC2 (giving an increase of 1-2 centimetres per allele) compromised proteolytic inhibition of PAPP-A and increased cleavage of IGFBP-4 in vitro, resulting in higher bioavailability of insulin-like growth factors. These 83 height-associated variants overlap genes that are mutated in monogenic growth disorders and highlight new biological candidates (such as ADAMTS3, IL11RA and NOX4) and pathways (such as proteoglycan and glycosaminoglycan synthesis) involved in growth. Our results demonstrate that sufficiently large sample sizes can uncover rare and low-frequency variants of moderate-to-large effect associated with polygenic human phenotypes, and that these variants implicate relevant genes and pathways.