The diversity in our genome is crucial to understanding the demographic history of worldwide populations. However, we have yet to know whether subtle genetic differences within a population can be disentangled, or whether they have an impact on complex traits. Here we apply dimensionality reduction methods (PCA, t-SNE, PCA-t-SNE, UMAP, and PCA-UMAP) to biobank-derived genomic data of a Japanese population (n = 169,719). Dimensionality reduction reveals fine-scale population structure, conspicuously differentiating adjacent insular subpopulations. We further enluciate the demographic landscape of these Japanese subpopulations using population genetics analyses. Finally, we perform phenome-wide polygenic risk score (PRS) analyses on 67 complex traits. Differences in PRS between the deconvoluted subpopulations are not always concordant with those in the observed phenotypes, suggesting that the PRS differences might reflect biases from the uncorrected structure, in a trait-dependent manner. This study suggests that such an uncorrected structure can be a potential pitfall in the clinical application of PRS.
Groundwater provides critical freshwater supply, particularly in dry regions where surface water availability is limited. Climate change impacts on GWS (groundwater storage) could affect the sustainability of freshwater resources. Here, we used a fully-coupled climate model to investigate GWS changes over seven critical aquifers identified as significantly distressed by satellite observations. We assessed the potential climate-driven impacts on GWS changes throughout the 21st century under the business-as-usual scenario (RCP8.5). Results show that the climate-driven impacts on GWS changes do not necessarily reflect the long-term trend in precipitation; instead, the trend may result from enhancement of evapotranspiration, and reduction in snowmelt, which collectively lead to divergent responses of GWS changes across different aquifers. Finally, we compare the climate-driven and anthropogenic pumping impacts. The reduction in GWS is mainly due to the combined impacts of over-pumping and climate effects; however, the contribution of pumping could easily far exceed the natural replenishment.
The realization of operationally stable blue organic light-emitting diodes is a challenging issue across the field. While device optimization has been a focus to effectively prolong device lifetime, strategies based on molecular engineering of chemical structures, particularly at the subatomic level, remains little. Herein, we explore the effect of targeted deuteration on donor and/or acceptor units of thermally activated delayed fluorescence emitters and investigate the structure-property relationship between intrinsic molecular stability, based on isotopic effect, and device operational stability. We show that the deuteration of the acceptor unit is critical to enhance the photostability of thermally activated delayed fluorescence compounds and hence device lifetime in addition to that of the donor units, which is commonly neglected due to the limited availability and synthetic complexity of deuterated acceptors. Based on these isotopic analogues, we observe a gradual increase in the device operational stability and achieve the long-lifetime time to 90% of the initial luminance of 23.4 h at the luminance of 1000 cd m-2 for thermally activated delayed fluorescence-sensitized organic light-emitting diodes. We anticipate our strategic deuteration approach provides insights and demonstrates the importance on structural modification materials at a subatomic level towards prolonging the device operational stability.
Polygenic risk scores (PRS) have been shown to predict breast cancer risk in European women, but their utility in Asian women is unclear. Here we evaluate the best performing PRSs for European-ancestry women using data from 17,262 breast cancer cases and 17,695 controls of Asian ancestry from 13 case-control studies, and 10,255 Chinese women from a prospective cohort (413 incident breast cancers). Compared to women in the middle quintile of the risk distribution, women in the highest 1% of PRS distribution have a ~2.7-fold risk and women in the lowest 1% of PRS distribution has ~0.4-fold risk of developing breast cancer. There is no evidence of heterogeneity in PRS performance in Chinese, Malay and Indian women. A PRS developed for European-ancestry women is also predictive of breast cancer risk in Asian women and can help in developing risk-stratified screening programmes in Asia.
GWAS have identified a breast cancer susceptibility locus on 2q35. Here we report the fine mapping of this locus using data from 101,943 subjects from 50 case-control studies. We genotype 276 SNPs using the 'iCOGS' genotyping array and impute genotypes for a further 1,284 using 1000 Genomes Project data. All but two, strongly correlated SNPs (rs4442975 G/T and rs6721996 G/A) are excluded as candidate causal variants at odds against >100:1. The best functional candidate, rs4442975, is associated with oestrogen receptor positive (ER+) disease with an odds ratio (OR) in Europeans of 0.85 (95% confidence interval=0.84-0.87; P=1.7 × 10(-43)) per t-allele. This SNP flanks a transcriptional enhancer that physically interacts with the promoter of IGFBP5 (encoding insulin-like growth factor-binding protein 5) and displays allele-specific gene expression, FOXA1 binding and chromatin looping. Evidence suggests that the g-allele confers increased breast cancer susceptibility through relative downregulation of IGFBP5, a gene with known roles in breast cell biology.
Steroid-resistant nephrotic syndrome (SRNS) causes 15% of chronic kidney disease (CKD). Here we show that recessive mutations in FAT1 cause a distinct renal disease entity in four families with a combination of SRNS, tubular ectasia, haematuria and facultative neurological involvement. Loss of FAT1 results in decreased cell adhesion and migration in fibroblasts and podocytes and the decreased migration is partially reversed by a RAC1/CDC42 activator. Podocyte-specific deletion of Fat1 in mice induces abnormal glomerular filtration barrier development, leading to podocyte foot process effacement. Knockdown of Fat1 in renal tubular cells reduces migration, decreases active RAC1 and CDC42, and induces defects in lumen formation. Knockdown of fat1 in zebrafish causes pronephric cysts, which is partially rescued by RAC1/CDC42 activators, confirming a role of the two small GTPases in the pathogenesis. These findings provide new insights into the pathogenesis of SRNS and tubulopathy, linking FAT1 and RAC1/CDC42 to podocyte and tubular cell function.
Recent literature on marine fish farming brands it as potentially compatible with sustainable resource use, conservation, and human nutrition goals, and aligns with the emerging policy discourse of 'blue growth'. We advance a two-pronged critique. First, contemporary narratives tend to overstate marine finfish aquaculture's potential to deliver food security and environmental sustainability. Second, they often align with efforts to enclose maritime space that could facilitate its allocation to extractive industries and conservation interests and exclude fishers. Policies and investments that seek to increase the availability and accessibility of affordable and sustainable farmed aquatic foods should focus on freshwater aquaculture.
To investigate large structural clonal mosaicism of chromosome X, we analysed the SNP microarray intensity data of 38,303 women from cancer genome-wide association studies (20,878 cases and 17,425 controls) and detected 124 mosaic X events >2 Mb in 97 (0.25%) women. Here we show rates for X-chromosome mosaicism are four times higher than mean autosomal rates; X mosaic events more often include the entire chromosome and participants with X events more likely harbour autosomal mosaic events. X mosaicism frequency increases with age (0.11% in 50-year olds; 0.45% in 75-year olds), as reported for Y and autosomes. Methylation array analyses of 33 women with X mosaicism indicate events preferentially involve the inactive X chromosome. Our results provide further evidence that the sex chromosomes undergo mosaic events more frequently than autosomes, which could have implications for understanding the underlying mechanisms of mosaic events and their possible contribution to risk for chronic diseases.
Prostate cancer is a polygenic disease with a large heritable component. A number of common, low-penetrance prostate cancer risk loci have been identified through GWAS. Here we apply the Bayesian multivariate variable selection algorithm JAM to fine-map 84 prostate cancer susceptibility loci, using summary data from a large European ancestry meta-analysis. We observe evidence for multiple independent signals at 12 regions and 99 risk signals overall. Only 15 original GWAS tag SNPs remain among the catalogue of candidate variants identified; the remainder are replaced by more likely candidates. Biological annotation of our credible set of variants indicates significant enrichment within promoter and enhancer elements, and transcription factor-binding sites, including AR, ERG and FOXA1. In 40 regions at least one variant is colocalised with an eQTL in prostate cancer tissue. The refined set of candidate variants substantially increase the proportion of familial relative risk explained by these known susceptibility regions, which highlights the importance of fine-mapping studies and has implications for clinical risk profiling.
A locus at 19p13 is associated with breast cancer (BC) and ovarian cancer (OC) risk. Here we analyse 438 SNPs in this region in 46,451 BC and 15,438 OC cases, 15,252 BRCA1 mutation carriers and 73,444 controls and identify 13 candidate causal SNPs associated with serous OC (P=9.2 × 10(-20)), ER-negative BC (P=1.1 × 10(-13)), BRCA1-associated BC (P=7.7 × 10(-16)) and triple negative BC (P-diff=2 × 10(-5)). Genotype-gene expression associations are identified for candidate target genes ANKLE1 (P=2 × 10(-3)) and ABHD8 (P<2 × 10(-3)). Chromosome conformation capture identifies interactions between four candidate SNPs and ABHD8, and luciferase assays indicate six risk alleles increased transactivation of the ADHD8 promoter. Targeted deletion of a region containing risk SNP rs56069439 in a putative enhancer induces ANKLE1 downregulation; and mRNA stability assays indicate functional effects for an ANKLE1 3'-UTR SNP. Altogether, these data suggest that multiple SNPs at 19p13 regulate ABHD8 and perhaps ANKLE1 expression, and indicate common mechanisms underlying breast and ovarian cancer risk.
A panel of influenza virus-like sequences were recently documented in fish and amphibians. Of these, the Wuhan spiny eel influenza virus (WSEIV) was found to phylogenetically cluster with influenza B viruses as a sister clade. Influenza B viruses have been documented to circulate only in humans, with certain virus isolates found in harbor seals. It is therefore interesting that a similar virus was potentially found in fish. Here we characterize the putative hemagglutinin (HA) and neuraminidase (NA) surface glycoproteins of the WSEIV. Functionally, we show that the WSEIV NA-like protein has sialidase activity comparable to B/Malaysia/2506/2004 influenza B virus NA, making it a bona fide neuraminidase that is sensitive to NA inhibitors. We tested the functionality of the HA by addressing the receptor specificity, stability, preferential airway protease cleavage, and fusogenicity. We show highly specific binding to monosialic ganglioside 2 (GM2) and fusogenicity at a range of different pH conditions. In addition, we found limited antigenic conservation of the WSEIV HA and NA relative to the B/Malaysia/2506/2004 virus HA and NA. In summary, we perform a functional and antigenic characterization of the glycoproteins of WSEIV to assess if it is indeed a bona fide influenza virus potentially circulating in ray-finned fish.
Overhunting reduces important plant-animal interactions such as vertebrate seed dispersal and seed predation, thereby altering plant regeneration and even above-ground biomass. It remains unclear, however, if non-hunted species can compensate for lost vertebrates in defaunated ecosystems. We use a nested exclusion experiment to isolate the effects of different seed enemies in a Bornean rainforest. In four of five tree species, vertebrates kill many seeds (13-66%). Nonetheless, when large mammals are excluded, seed mortality from insects and fungi fully compensates for the lost vertebrate predation, such that defaunation has no effect on seedling establishment. The switch from seed predation by generalist vertebrates to specialist insects and fungi in defaunated systems may alter Janzen-Connell effects and density-dependence in plants. Previous work using simulation models to explore how lost seed dispersal will affect tree species composition and carbon storage may require reevaluation in the context of functional redundancy within complex species interactions networks.
In 2020, 146,063 deaths due to pancreatic cancer are estimated to occur in Europe and the United States combined. To identify common susceptibility alleles, we performed the largest pancreatic cancer GWAS to date, including 9040 patients and 12,496 controls of European ancestry from the Pancreatic Cancer Cohort Consortium (PanScan) and the Pancreatic Cancer Case-Control Consortium (PanC4). Here, we find significant evidence of a novel association at rs78417682 (7p12/TNS3, P = 4.35 × 10-8). Replication of 10 promising signals in up to 2737 patients and 4752 controls from the PANcreatic Disease ReseArch (PANDoRA) consortium yields new genome-wide significant loci: rs13303010 at 1p36.33 (NOC2L, P = 8.36 × 10-14), rs2941471 at 8q21.11 (HNF4G, P = 6.60 × 10-10), rs4795218 at 17q12 (HNF1B, P = 1.32 × 10-8), and rs1517037 at 18q21.32 (GRP, P = 3.28 × 10-8). rs78417682 is not statistically significantly associated with pancreatic cancer in PANDoRA. Expression quantitative trait locus analysis in three independent pancreatic data sets provides molecular support of NOC2L as a pancreatic cancer susceptibility gene.
The incidence of Plasmodium vivax infection has declined markedly in Malaysia over the past decade despite evidence of high-grade chloroquine resistance. Here we investigate the genetic changes in a P. vivax population approaching elimination in 51 isolates from Sabah, Malaysia and compare these with data from 104 isolates from Thailand and 104 isolates from Indonesia. Sabah displays extensive population structure, mirroring that previously seen with the emergence of artemisinin-resistant P. falciparum founder populations in Cambodia. Fifty-four percent of the Sabah isolates have identical genomes, consistent with a rapid clonal expansion. Across Sabah, there is a high prevalence of loci known to be associated with antimalarial drug resistance. Measures of differentiation between the three countries reveal several gene regions under putative selection in Sabah. Our findings highlight important factors pertinent to parasite resurgence and molecular cues that can be used to monitor low-endemic populations at the end stages of P. vivax elimination.
Small populations are often exposed to high inbreeding and mutational load that can increase the risk of extinction. The Sumatran rhinoceros was widespread in Southeast Asia, but is now restricted to small and isolated populations on Sumatra and Borneo, and most likely extinct on the Malay Peninsula. Here, we analyse 5 historical and 16 modern genomes from these populations to investigate the genomic consequences of the recent decline, such as increased inbreeding and mutational load. We find that the Malay Peninsula population experienced increased inbreeding shortly before extirpation, which possibly was accompanied by purging. The populations on Sumatra and Borneo instead show low inbreeding, but high mutational load. The currently small population sizes may thus in the near future lead to inbreeding depression. Moreover, we find little evidence for differences in local adaptation among populations, suggesting that future inbreeding depression could potentially be mitigated by assisted gene flow among populations.
As one of the great survivors of the plant kingdom, barnyard grasses (Echinochloa spp.) are the most noxious and common weeds in paddy ecosystems. Meanwhile, at least two Echinochloa species have been domesticated and cultivated as millets. In order to better understand the genomic forces driving the evolution of Echinochloa species toward weed and crop characteristics, we assemble genomes of three Echinochloa species (allohexaploid E. crus-galli and E. colona, and allotetraploid E. oryzicola) and re-sequence 737 accessions of barnyard grasses and millets from 16 rice-producing countries. Phylogenomic and comparative genomic analyses reveal the complex and reticulate evolution in the speciation of Echinochloa polyploids and provide evidence of constrained disease-related gene copy numbers in Echinochloa. A population-level investigation uncovers deep population differentiation for local adaptation, multiple target-site herbicide resistance mutations of barnyard grasses, and limited domestication of barnyard millets. Our results provide genomic insights into the dual roles of Echinochloa species as weeds and crops as well as essential resources for studying plant polyploidization, adaptation, precision weed control and millet improvements.
Chromosome 8q24 is a susceptibility locus for multiple cancers, including prostate cancer. Here we combine genetic data across the 8q24 susceptibility region from 71,535 prostate cancer cases and 52,935 controls of European ancestry to define the overall contribution of germline variation at 8q24 to prostate cancer risk. We identify 12 independent risk signals for prostate cancer (p
Protected Areas (PAs) are the cornerstone of biodiversity conservation. Here, we collated distributional data for >14,000 (~70% of) species of amphibians and reptiles (herpetofauna) to perform a global assessment of the conservation effectiveness of PAs using species distribution models. Our analyses reveal that >91% of herpetofauna species are currently distributed in PAs, and that this proportion will remain unaltered under future climate change. Indeed, loss of species' distributional ranges will be lower inside PAs than outside them. Therefore, the proportion of effectively protected species is predicted to increase. However, over 7.8% of species currently occur outside PAs, and large spatial conservation gaps remain, mainly across tropical and subtropical moist broadleaf forests, and across non-high-income countries. We also predict that more than 300 amphibian and 500 reptile species may go extinct under climate change over the course of the ongoing century. Our study highlights the importance of PAs in providing herpetofauna with refuge from climate change, and suggests ways to optimize PAs to better conserve biodiversity worldwide.
Comparative whole-genome sequencing enables the identification of specific mutations during adaptation of bacteria to new environments and allelic replacement can establish their causality. However, the mechanisms of action are hard to decipher and little has been achieved for epistatic mutations, especially at the metabolic level. Here we show that a strain of Escherichia coli carrying mutations in the rpoC and glpK genes, derived from adaptation in glycerol, uses two distinct metabolic strategies to gain growth advantage. A 27-bp deletion in the rpoC gene first increases metabolic efficiency. Then, a point mutation in the glpK gene promotes growth by improving glycerol utilization but results in increased carbon wasting as overflow metabolism. In a strain carrying both mutations, these contrasting carbon/energy saving and wasting mechanisms work together to give an 89% increase in growth rate. This study provides insight into metabolic reprogramming during adaptive laboratory evolution for fast cellular growth.
Antimicrobial resistance (AMR) is a serious threat to global public health, but obtaining representative data on AMR for healthy human populations is difficult. Here, we use metagenomic analysis of untreated sewage to characterize the bacterial resistome from 79 sites in 60 countries. We find systematic differences in abundance and diversity of AMR genes between Europe/North-America/Oceania and Africa/Asia/South-America. Antimicrobial use data and bacterial taxonomy only explains a minor part of the AMR variation that we observe. We find no evidence for cross-selection between antimicrobial classes, or for effect of air travel between sites. However, AMR gene abundance strongly correlates with socio-economic, health and environmental factors, which we use to predict AMR gene abundances in all countries in the world. Our findings suggest that global AMR gene diversity and abundance vary by region, and that improving sanitation and health could potentially limit the global burden of AMR. We propose metagenomic analysis of sewage as an ethically acceptable and economically feasible approach for continuous global surveillance and prediction of AMR.