Escherichia albertii is a recently identified gastrointestinal bacterial pathogen of humans and animals which is typically misidentified as pathotypes of diarrhoeagenic Escherichia coli or Shigella species and is generally only detected during genomic surveillance of other Enterobacteriaceae. The incidence of E. albertii is likely underestimated, and its epidemiology and clinical relevance are poorly characterised. Here, we whole genome sequenced E. albertii isolates from humans (n = 83) and birds (n = 79) isolated in Great Britain between 2000 and 2021 and analysed these alongside a broader public dataset (n = 475) to address these gaps. We found human and avian isolates typically (90%; 148/164) belonged to host-associated monophyletic groups with distinct virulence and antimicrobial resistance profiles. Overlaid patient epidemiological data suggested that human infection was likely related to travel and possibly foodborne transmission. The Shiga toxin encoding stx2f gene was associated with clinical disease (OR = 10.27, 95% CI = 2.98-35.45 p = 0.0002) in finches. Our results suggest that improved future surveillance will further elucidate disease ecology and public and animal health risks associated with E. albertii.

The risks of cooling water shortages to thermo-electric power plants are increasingly studied as an important climate risk to the energy sector. Whilst electricity transmission networks reduce the risks during disruptions, more costly plants must provide alternative supplies. Here, we investigate the electricity price impacts of cooling water shortages on Britain's power supplies using a probabilistic spatial risk model of regional climate, hydrological droughts and cooling water shortages, coupled with an economic model of electricity supply, demand and prices. We find that on extreme days (p99), almost 50% (7GWe) of freshwater thermal capacity is unavailable. Annualized cumulative costs on electricity prices range from £29-66m.yr-1 GBP2018, whilst in 20% of cases from £66-95m.yr-1. With climate change, the median annualized impact exceeds £100m.yr-1. The single year impacts of a 1-in-25 year event exceed >£200m, indicating the additional investments justifiable to mitigate the 1st-order economic risks of cooling water shortage during droughts.

Humans survive and thrive through social exchange. Yet, social dependency also comes at a cost. Perceived social isolation, or loneliness, affects physical and mental health, cognitive performance, overall life expectancy, and increases vulnerability to Alzheimer's disease-related dementias. Despite severe consequences on behavior and health, the neural basis of loneliness remains elusive. Using the UK Biobank population imaging-genetics cohort (n = ~40,000, aged 40-69 years when recruited, mean age = 54.9), we test for signatures of loneliness in grey matter morphology, intrinsic functional coupling, and fiber tract microstructure. The loneliness-linked neurobiological profiles converge on a collection of brain regions known as the 'default network'. This higher associative network shows more consistent loneliness associations in grey matter volume than other cortical brain networks. Lonely individuals display stronger functional communication in the default network, and greater microstructural integrity of its fornix pathway. The findings fit with the possibility that the up-regulation of these neural circuits supports mentalizing, reminiscence and imagination to fill the social void.

The recent years have seen a growing number of studies investigating evolutionary questions using ancient DNA. To address these questions, one of the most frequently-used method is principal component analysis (PCA). When PCA is applied to temporal samples, the sample dates are, however, ignored during analysis, leading to imperfect representations of samples in PC plots. Here, we present a factor analysis (FA) method in which individual scores are corrected for the effect of allele frequency drift over time. We obtained exact solutions for the estimates of corrected factors, and we provided a fast algorithm for their computation. Using computer simulations and ancient European samples, we compared geometric representations obtained from FA with PCA and with ancestry estimation programs. In admixture analyses, FA estimates agreed with tree-based statistics, and they were more accurate than those obtained from PCA projections and from ancestry estimation programs. A great advantage of FA over existing approaches is to improve descriptive analyses of ancient DNA samples without requiring inclusion of outgroup or present-day samples.

The fossil record documents the history of life, but the reliability of that record has often been questioned. Spatiotemporal variability in sedimentary rock volume, sampling and research effort especially frustrates global-scale diversity reconstructions. Various proposals have been made to rectify palaeodiversity estimates using proxy measures for the availability and sampling of the rock record, but the validity of these approaches remains controversial. Targeting the rich fossil record of Great Britain as a highly detailed regional exemplar, our statistical analysis shows that marine outcrop area contains a signal useful for predicting changes in diversity, collections and formations, whereas terrestrial outcrop area contains a signal useful for predicting formations. In contrast, collection and formation counts are information redundant with fossil richness, characterized by symmetric, bidirectional information flow. If this is true, the widespread use of collection and formation counts as sampling proxies to correct the raw palaeodiversity data may be unwarranted.

Here we report the results of multiple analytical techniques on sub-mm particulate material derived from Unit 1 of the Fukushima Daiichi Nuclear Power Plant to provide a better understanding of the events that occurred and the environmental legacy. Through combined x-ray fluorescence and absorption contrast micro-focused x-ray tomography, entrapped U particulate are observed to exist around the exterior circumference of the highly porous Si-based particle. Further synchrotron radiation analysis of a number of these entrapped particles shows them to exist as UO2-identical to reactor fuel, with confirmation of their nuclear origin shown via mass spectrometry analysis. While unlikely to represent an environmental or health hazard, such assertions would likely change should break-up of the Si-containing bulk particle occur. However, more important to the long-term decommissioning of the reactors at the FDNPP (and environmental clean-upon), is the knowledge that core integrity of reactor Unit 1 was compromised with nuclear material existing outside of the reactors primary containment.

Socioeconomic position (SEP) is a multi-dimensional construct reflecting (and influencing) multiple socio-cultural, physical, and environmental factors. In a sample of 286,301 participants from UK Biobank, we identify 30 (29 previously unreported) independent-loci associated with income. Using a method to meta-analyze data from genetically-correlated traits, we identify an additional 120 income-associated loci. These loci show clear evidence of functionality, with transcriptional differences identified across multiple cortical tissues, and links to GABAergic and serotonergic neurotransmission. By combining our genome wide association study on income with data from eQTL studies and chromatin interactions, 24 genes are prioritized for follow up, 18 of which were previously associated with intelligence. We identify intelligence as one of the likely causal, partly-heritable phenotypes that might bridge the gap between molecular genetic inheritance and phenotypic consequence in terms of income differences. These results indicate that, in modern era Great Britain, genetic effects contribute towards some of the observed socioeconomic inequalities.

Despite improvements in medical and surgical therapies, a significant portion of patients with critical limb ischemia (CLI) are considered as "no option" for revascularization. In this work, a nitric oxide (NO)-boosted and activated nanovesicle regeneration kit (n-BANK) is constructed by decorating stem cell-derived nanoscale extracellular vesicles with NO nanocages. Our results demonstrate that n-BANKs could store NO in endothelial cells for subsequent release upon pericyte recruitment for CLI revascularization. Notably, n-BANKs enable endothelial cells to trigger eNOS activation and form tube-like structures. Subsequently, eNOS-derived NO robustly recruits pericytes to invest nascent endothelial cell tubes, giving rise to mature blood vessels. Consequently, n-BANKs confer complete revascularization in female mice following CLI, and thereby achieve limb preservation and restore the motor function. In light of n-BANK evoking pericyte-endothelial interactions to create functional vascular networks, it features promising therapeutic potential in revascularization to reduce CLI-related amputations, which potentially impact regeneration medicine.

Rapid-eye movement (REM) sleep behavior disorder (RBD), enactment of dreams during REM sleep, is an early clinical symptom of alpha-synucleinopathies and defines a more severe subtype. The genetic background of RBD and its underlying mechanisms are not well understood. Here, we perform a genome-wide association study of RBD, identifying five RBD risk loci near SNCA, GBA, TMEM175, INPP5F, and SCARB2. Expression analyses highlight SNCA-AS1 and potentially SCARB2 differential expression in different brain regions in RBD, with SNCA-AS1 further supported by colocalization analyses. Polygenic risk score, pathway analysis, and genetic correlations provide further insights into RBD genetics, highlighting RBD as a unique alpha-synucleinopathy subpopulation that will allow future early intervention.

There is increasing evidence for a strong inherited genetic basis of susceptibility to acute lymphoblastic leukaemia (ALL) in children. To identify new risk variants for B-cell ALL (B-ALL) we conducted a meta-analysis with four GWAS (genome-wide association studies), totalling 5321 cases and 16,666 controls of European descent. We herein describe novel risk loci for B-ALL at 9q21.31 (rs76925697, P = 2.11 × 10-8), for high-hyperdiploid ALL at 5q31.1 (rs886285, P = 1.56 × 10-8) and 6p21.31 (rs210143 in BAK1, P = 2.21 × 10-8), and ETV6-RUNX1 ALL at 17q21.32 (rs10853104 in IGF2BP1, P = 1.82 × 10-8). Particularly notable are the pleiotropic effects of the BAK1 variant on multiple haematological malignancies and specific effects of IGF2BP1 on ETV6-RUNX1 ALL evidenced by both germline and somatic genomic analyses. Integration of GWAS signals with transcriptomic/epigenomic profiling and 3D chromatin interaction data for these leukaemia risk loci suggests deregulation of B-cell development and the cell cycle as central mechanisms governing genetic susceptibility to ALL.

Acute myeloid leukemia (AML) is a hematological malignancy with an undefined heritable risk. Here we perform a meta-analysis of three genome-wide association studies, with replication in a fourth study, incorporating a total of 4018 AML cases and 10488 controls. We identify a genome-wide significant risk locus for AML at 11q13.2 (rs4930561; P = 2.15 × 10-8; KMT5B). We also identify a genome-wide significant risk locus for the cytogenetically normal AML sub-group (N = 1287) at 6p21.32 (rs3916765; P = 1.51 × 10-10; HLA). Our results inform on AML etiology and identify putative functional genes operating in histone methylation (KMT5B) and immune function (HLA).

Schizophrenia and bipolar disorder are major psychiatric disorders with high heritability and overlapping genetic variance. Here we perform a genome-wide association study in an ethnically homogeneous cohort of 904 schizophrenia cases and 1,640 controls drawn from the Ashkenazi Jewish population. We identify a novel genome-wide significant risk locus at chromosome 4q26, demonstrating the potential advantages of this founder population for gene discovery. The top single-nucleotide polymorphism (SNP; rs11098403) demonstrates consistent effects across 11 replication and extension cohorts, totalling 23, 191 samples across multiple ethnicities, regardless of diagnosis (schizophrenia or bipolar disorder), resulting in Pmeta=9.49 × 10(-12) (odds ratio (OR)=1.13, 95% confidence interval (CI): 1.08-1.17) across both disorders and Pmeta=2.67 × 10(-8) (OR=1.15, 95% CI: 1.08-1.21) for schizophrenia alone. In addition, this intergenic SNP significantly predicts postmortem cerebellar gene expression of NDST3, which encodes an enzyme critical to heparan sulphate metabolism. Heparan sulphate binding is critical to neurite outgrowth, axon formation and synaptic processes thought to be aberrant in these disorders.

Genetic interactions have been reported to underlie phenotypes in a variety of systems, but the extent to which they contribute to complex disease in humans remains unclear. In principle, genome-wide association studies (GWAS) provide a platform for detecting genetic interactions, but existing methods for identifying them from GWAS data tend to focus on testing individual locus pairs, which undermines statistical power. Importantly, a global genetic network mapped for a model eukaryotic organism revealed that genetic interactions often connect genes between compensatory functional modules in a highly coherent manner. Taking advantage of this expected structure, we developed a computational approach called BridGE that identifies pathways connected by genetic interactions from GWAS data. Applying BridGE broadly, we discover significant interactions in Parkinson's disease, schizophrenia, hypertension, prostate cancer, breast cancer, and type 2 diabetes. Our novel approach provides a general framework for mapping complex genetic networks underlying human disease from genome-wide genotype data.

Genetic correlation analysis has quickly gained popularity in the past few years and provided insights into the genetic etiology of numerous complex diseases. However, existing approaches oversimplify the shared genetic architecture between different phenotypes and cannot effectively identify precise genetic regions contributing to the genetic correlation. In this work, we introduce LOGODetect, a powerful and efficient statistical method to identify small genome segments harboring local genetic correlation signals. LOGODetect automatically identifies genetic regions showing consistent associations with multiple phenotypes through a scan statistic approach. It uses summary association statistics from genome-wide association studies (GWAS) as input and is robust to sample overlap between studies. Applied to seven phenotypically distinct but genetically correlated neuropsychiatric traits, we identify 227 non-overlapping genome regions associated with multiple traits, including multiple hub regions showing concordant effects on five or more traits. Our method addresses critical limitations in existing analytic strategies and may have wide applications in post-GWAS analysis.

Genome-wide association studies (GWAS) have cataloged many significant associations between genetic variants and complex traits. However, most of these findings have unclear biological significance, because they often have small effects and occur in non-coding regions. Integration of GWAS with gene regulatory networks addresses both issues by aggregating weak genetic signals within regulatory programs. Here we develop a Bayesian framework that integrates GWAS summary statistics with regulatory networks to infer genetic enrichments and associations simultaneously. Our method improves upon existing approaches by explicitly modeling network topology to assess enrichments, and by automatically leveraging enrichments to identify associations. Applying this method to 18 human traits and 38 regulatory networks shows that genetic signals of complex traits are often enriched in interconnections specific to trait-relevant cell types or tissues. Prioritizing variants within enriched networks identifies known and previously undescribed trait-associated genes revealing biological and therapeutic insights.

Using genotype data to perform accurate genetic prediction of complex traits can facilitate genomic selection in animal and plant breeding programs, and can aid in the development of personalized medicine in humans. Because most complex traits have a polygenic architecture, accurate genetic prediction often requires modeling all genetic variants together via polygenic methods. Here, we develop such a polygenic method, which we refer to as the latent Dirichlet process regression model. Dirichlet process regression is non-parametric in nature, relies on the Dirichlet process to flexibly and adaptively model the effect size distribution, and thus enjoys robust prediction performance across a broad spectrum of genetic architectures. We compare Dirichlet process regression with several commonly used prediction methods with simulations. We further apply Dirichlet process regression to predict gene expressions, to conduct PrediXcan based gene set test, to perform genomic selection of four traits in two species, and to predict eight complex traits in a human cohort.Genetic prediction of complex traits with polygenic architecture has wide application from animal breeding to disease prevention. Here, Zeng and Zhou develop a non-parametric genetic prediction method based on latent Dirichlet Process regression models.

Epithelial ovarian cancer (EOC) has a heritable component that remains to be fully characterized. Most identified common susceptibility variants lie in non-protein-coding sequences. We hypothesized that variants in the 3' untranslated region at putative microRNA (miRNA)-binding sites represent functional targets that influence EOC susceptibility. Here, we evaluate the association between 767 miRNA-related single-nucleotide polymorphisms (miRSNPs) and EOC risk in 18,174 EOC cases and 26,134 controls from 43 studies genotyped through the Collaborative Oncological Gene-environment Study. We identify several miRSNPs associated with invasive serous EOC risk (odds ratio=1.12, P=10(-8)) mapping to an inversion polymorphism at 17q21.31. Additional genotyping of non-miRSNPs at 17q21.31 reveals stronger signals outside the inversion (P=10(-10)). Variation at 17q21.31 is associated with neurological diseases, and our collaboration is the first to report an association with EOC susceptibility. An integrated molecular analysis in this region provides evidence for ARHGAP27 and PLEKHM1 as candidate EOC susceptibility genes.

Integrating results from genome-wide association studies (GWASs) and gene expression studies through transcriptome-wide association study (TWAS) has the potential to shed light on the causal molecular mechanisms underlying disease etiology. Here, we present a probabilistic Mendelian randomization (MR) method, PMR-Egger, for TWAS applications. PMR-Egger relies on a MR likelihood framework that unifies many existing TWAS and MR methods, accommodates multiple correlated instruments, tests the causal effect of gene on trait in the presence of horizontal pleiotropy, and is scalable to hundreds of thousands of individuals. In simulations, PMR-Egger provides calibrated type I error control for causal effect testing in the presence of horizontal pleiotropic effects, is reasonably robust under various types of model misspecifications, is more powerful than existing TWAS/MR approaches, and can directly test for horizontal pleiotropy. We illustrate the benefits of PMR-Egger in applications to 39 diseases and complex traits obtained from three GWASs including the UK Biobank.

HNF1B is overexpressed in clear cell epithelial ovarian cancer, and we observed epigenetic silencing in serous epithelial ovarian cancer, leading us to hypothesize that variation in this gene differentially associates with epithelial ovarian cancer risk according to histological subtype. Here we comprehensively map variation in HNF1B with respect to epithelial ovarian cancer risk and analyse DNA methylation and expression profiles across histological subtypes. Different single-nucleotide polymorphisms associate with invasive serous (rs7405776 odds ratio (OR)=1.13, P=3.1 × 10(-10)) and clear cell (rs11651755 OR=0.77, P=1.6 × 10(-8)) epithelial ovarian cancer. Risk alleles for the serous subtype associate with higher HNF1B-promoter methylation in these tumours. Unmethylated, expressed HNF1B, primarily present in clear cell tumours, coincides with a CpG island methylator phenotype affecting numerous other promoters throughout the genome. Different variants in HNF1B associate with risk of serous and clear cell epithelial ovarian cancer; DNA methylation and expression patterns are also notably distinct between these subtypes. These findings underscore distinct mechanisms driving different epithelial ovarian cancer histological subtypes.