The Allelic Landscape of Human Blood Cell Trait Variation and Links to Common Complex Disease.

Many common variants have been associated with hematological traits, but identification of causal genes and pathways has proven challenging. We performed a genome-wide association analysis in the UK Biobank and INTERVAL studies, testing 29.5 million genetic variants for association with 36 red cell, white cell, and platelet properties in 173,480 European-ancestry participants. This effort yielded hundreds of low frequency (<5%) and rare (<1%) variants with a strong impact on blood cell phenotypes. Our data highlight general properties of the allelic architecture of complex traits, including the proportion of the heritable component of each blood trait explained by the polygenic signal across different genome regulatory domains. Finally, through Mendelian randomization, we provide evidence of shared genetic pathways linking blood cell indices with complex pathologies, including autoimmune diseases, schizophrenia, and coronary heart disease and evidence suggesting previously reported population associations between blood cell indices and cardiovascular disease may be non-causal.

Polygenic risk score adds to a clinical risk score in the prediction of cardiovascular disease in a clinical setting.

BACKGROUND AND AIMS: A cardiovascular disease polygenic risk score (CVD-PRS) can stratify individuals into different categories of cardiovascular risk, but whether the addition of a CVD-PRS to clinical risk scores improves the identification of individuals at increased risk in a real-world clinical setting is unknown. METHODS: The Genetics and the Vascular Health Check Study (GENVASC) was embedded within the UK National Health Service Health Check (NHSHC) programme which invites individuals between 40-74 years of age without known CVD to attend an assessment in a UK general practice where CVD risk factors are measured and a CVD risk score (QRISK2) is calculated. Between 2012-2020, 44,141 individuals (55.7% females, 15.8% non-white) who attended an NHSHC in 147 participating practices across two counties in England were recruited and followed. When 195 individuals (cases) had suffered a major CVD event (CVD death, myocardial infarction or acute coronary syndrome, coronary revascularisation, stroke), 396 propensity-matched controls with a similar risk profile were identified, and a nested case-control genetic study undertaken to see if the addition of a CVD-PRS to QRISK2 in the form of an integrated risk tool (IRT) combined with QRISK2 would have identified more individuals at the time of their NHSHC as at high risk (QRISK2 10-year CVD risk of ≥10%), compared with QRISK2 alone. RESULTS: The distribution of the standardised CVD-PRS was significantly different in cases compared with controls (cases mean score .32; controls, -.18, P = 8.28×10-9). QRISK2 identified 61.5% (95% confidence interval [CI]: 54.3%-68.4%) of individuals who subsequently developed a major CVD event as being at high risk at their NHSHC, while the combination of QRISK2 and IRT identified 68.7% (95% CI: 61.7%-75.2%), a relative increase of 11.7% (P = 1×10-4). The odds ratio (OR) of being up-classified was 2.41 (95% CI: 1.03-5.64, P = .031) for cases compared with controls. In individuals aged 40-54 years, QRISK2 identified 26.0% (95% CI: 16.5%-37.6%) of those who developed a major CVD event, while the combination of QRISK2 and IRT identified 38.4% (95% CI: 27.2%-50.5%), indicating a stronger relative increase of 47.7% in the younger age group (P = .001). The combination of QRISK2 and IRT increased the proportion of additional cases identified similarly in women as in men, and in non-white ethnicities compared with white ethnicity. The findings were similar when the CVD-PRS was added to the atherosclerotic cardiovascular disease pooled cohort equations (ASCVD-PCE) or SCORE2 clinical scores. CONCLUSIONS: In a clinical setting, the addition of genetic information to clinical risk assessment significantly improved the identification of individuals who went on to have a major CVD event as being at high risk, especially among younger individuals. The findings provide important real-world evidence of the potential value of implementing a CVD-PRS into health systems.

The influence of rare variants in circulating metabolic biomarkers.

Circulating metabolite levels are biomarkers for cardiovascular disease (CVD). Here we studied, association of rare variants and 226 serum lipoproteins, lipids and amino acids in 7,142 (discovery plus follow-up) healthy participants. We leveraged the information from multiple metabolite measurements on the same participants to improve discovery in rare variant association analyses for gene-based and gene-set tests by incorporating correlated metabolites as covariates in the validation stage. Gene-based analysis corrected for the effective number of tests performed, confirmed established associations at APOB, APOC3, PAH, HAL and PCSK (p<1.32x10-7) and identified novel gene-trait associations at a lower stringency threshold with ACSL1, MYCN, FBXO36 and B4GALNT3 (p<2.5x10-6). Regulation of the pyruvate dehydrogenase (PDH) complex was associated for the first time, in gene-set analyses also corrected for effective number of tests, with IDL and LDL parameters, as well as circulating cholesterol (pMETASKAT<2.41x10-6). In conclusion, using an approach that leverages metabolite measurements obtained in the same participants, we identified novel loci and pathways involved in the regulation of these important metabolic biomarkers. As large-scale biobanks continue to amass sequencing and phenotypic information, analytical approaches such as ours will be useful to fully exploit the copious amounts of biological data generated in these efforts.

Genetic architecture of human thinness compared to severe obesity.

The variation in weight within a shared environment is largely attributable to genetic factors. Whilst many genes/loci confer susceptibility to obesity, little is known about the genetic architecture of healthy thinness. Here, we characterise the heritability of thinness which we found was comparable to that of severe obesity (h2 = 28.07 vs 32.33% respectively), although with incomplete genetic overlap (r = -0.49, 95% CI [-0.17, -0.82], p = 0.003). In a genome-wide association analysis of thinness (n = 1,471) vs severe obesity (n = 1,456), we identified 10 loci previously associated with obesity, and demonstrate enrichment for established BMI-associated loci (pbinomial = 3.05x10-5). Simulation analyses showed that different association results between the extremes were likely in agreement with additive effects across the BMI distribution, suggesting different effects on thinness and obesity could be due to their different degrees of extremeness. In further analyses, we detected a novel obesity and BMI-associated locus at PKHD1 (rs2784243, obese vs. thin p = 5.99x10-6, obese vs. controls p = 2.13x10-6 pBMI = 2.3x10-13), associations at loci recently discovered with much larger sample sizes (e.g. FAM150B and PRDM6-CEP120), and novel variants driving associations at previously established signals (e.g. rs205262 at the SNRPC/C6orf106 locus and rs112446794 at the PRDM6-CEP120 locus). Our ability to replicate loci found with much larger sample sizes demonstrates the value of clinical extremes and suggest that characterisation of the genetics of thinness may provide a more nuanced understanding of the genetic architecture of body weight regulation and may inform the identification of potential anti-obesity targets.

Deletion of the 2-acyl-glycerophosphoethanolamine cycle improve glucose metabolism in Escherichia coli strains employed for overproduction of aromatic compounds.

BACKGROUND: As a metabolic engineering tool, an adaptive laboratory evolution (ALE) experiment was performed to increase the specific growth rate (µ) in an Escherichia coli strain lacking PTS, originally engineered to increase the availability of intracellular phosphoenolpyruvate and redirect to the aromatic biosynthesis pathway. As result, several evolved strains increased their growth fitness on glucose as the only carbon source. Two of these clones isolated at 120 and 200 h during the experiment, increased their µ by 338 and 373 %, respectively, compared to the predecessor PB11 strain. The genome sequence and analysis of the genetic changes of these two strains (PB12 and PB13) allowed for the identification of a novel strategy to enhance carbon utilization to overcome the absence of the major glucose transport system. RESULTS: Genome sequencing data of evolved strains revealed the deletion of chromosomal region of 10,328 pb and two punctual non-synonymous mutations in the dhaM and glpT genes, which occurred prior to their divergence during the early stages of the evolutionary process. Deleted genes related to increased fitness in the evolved strains are rppH, aas, lplT and galR. Furthermore, the loss of mutH, which was also lost during the deletion event, caused a 200-fold increase in the mutation rate. CONCLUSIONS: During the ALE experiment, both PB12 and PB13 strains lost the galR and rppH genes, allowing the utilization of an alternative glucose transport system and allowed enhanced mRNA half-life of many genes involved in the glycolytic pathway resulting in an increment in the µ of these derivatives. Finally, we demonstrated the deletion of the aas-lplT operon, which codes for the main components of the phosphatidylethanolamine turnover metabolism increased the further fitness and glucose uptake in these evolved strains by stimulating the phospholipid degradation pathway. This is an alternative mechanism to its regeneration from 2-acyl-glycerophosphoethanolamine, whose utilization improved carbon metabolism likely by the elimination of a futile cycle under certain metabolic conditions. The origin and widespread occurrence of a mutated population during the ALE indicates a strong stress condition present in strains lacking PTS and the plasticity of this bacterium that allows it to overcome hostile conditions.

A polygenic risk score added to a QRISK®2 cardiovascular disease risk calculator demonstrated robust clinical acceptance and clinical utility in the primary care setting.

AIMS: The aim of the study was to assess the real-world feasibility, acceptability, and impact of an integrated risk tool for cardiovascular disease (CVD IRT, combining the standard QRISK®2 risk algorithm with a polygenic risk score), implemented within routine primary practice in the UK National Health Service. METHODS AND RESULTS: The Healthcare Evaluation of Absolute Risk Testing Study (NCT05294419) evaluated participants undergoing primary care health checks. Both QRISK2 and CVD IRT scores were returned to the healthcare providers (HCPs), who then communicated the results to participants. The primary outcome of the study was feasibility of CVD IRT implementation. Secondary outcomes included changes in CVD risk (QRISK2 vs. CVD IRT) and impact of the CVD IRT on clinical decision-making. A total of 832 eligible participants (median age 55 years, 62% females, 97.5% White ethnicity) were enrolled across 12 UK primary care practices. Cardiovascular disease IRT scores were obtained on 100% of the blood samples. Healthcare providers stated that the CVD IRT could be incorporated into routine primary care in a straightforward manner in 90.7% of reports. Participants stated they were 'likely' or 'very likely' to recommend the use of this test to their family or friends in 86.9% of reports. Participants stated that the test was personally useful (98.8%) and that the results were easy to understand (94.6%). When CVD IRT exceeded QRISK2, HCPs planned changes in management for 108/388 (27.8%) of participants and 47% (62/132) of participants with absolute risk score changes of >2%. CONCLUSION: Amongst HCPs and participants who agreed to the trial of genetic data for refinement of clinical risk prediction in primary care, we observed that CVD IRT implementation was feasible and well accepted. The CVD IRT results were associated with planned changes in prevention strategies.

When a standard cardiovascular risk tool, as currently used in National Health Service Health Checks, was expanded to include genetic risk information, it was well accepted by both participants and healthcare providers and generated impactful changes in planned clinical decision-making.Most participants found the test useful and easy to understand, and healthcare providers found it straightforward to use in most cases.When risk was increased by the addition of genetic information, this influenced planned management decisions.

A systematic evaluation of the performance and properties of the UK Biobank Polygenic Risk Score (PRS) Release.

We assess the UK Biobank (UKB) Polygenic Risk Score (PRS) Release, a set of PRSs for 28 diseases and 25 quantitative traits that has been made available on the individuals in UKB, using a unified pipeline for PRS evaluation. We also release a benchmarking software tool to enable like-for-like performance evaluation for different PRSs for the same disease or trait. Extensive benchmarking shows the PRSs in the UKB Release to outperform a broad set of 76 published PRSs. For many of the diseases and traits we also validate the PRS algorithms in a separate cohort (100,000 Genomes Project). The availability of PRSs for 53 traits on the same set of individuals also allows a systematic assessment of their properties, and the increased power of these PRSs increases the evidence for their potential clinical benefit.

Genetic changes during a laboratory adaptive evolution process that allowed fast growth in glucose to an Escherichia coli strain lacking the major glucose transport system.

BACKGROUND: Escherichia coli strains lacking the phosphoenolpyruvate: carbohydrate phosphotransferase system (PTS), which is the major bacterial component involved in glucose transport and its phosphorylation, accumulate high amounts of phosphoenolpyruvate that can be diverted to the synthesis of commercially relevant products. However, these strains grow slowly in glucose as sole carbon source due to its inefficient transport and metabolism. Strain PB12, with 400% increased growth rate, was isolated after a 120 hours adaptive laboratory evolution process for the selection of faster growing derivatives in glucose. Analysis of the genetic changes that occurred in the PB12 strain that lacks PTS will allow a better understanding of the basis of its growth adaptation and, therefore, in the design of improved metabolic engineering strategies for enhancing carbon diversion into the aromatic pathways. RESULTS: Whole genome analyses using two different sequencing methodologies: the Roche NimbleGen Inc. comparative genome sequencing technique, and high throughput sequencing with Illumina Inc. GAIIx, allowed the identification of the genetic changes that occurred in the PB12 strain. Both methods detected 23 non-synonymous and 22 synonymous point mutations. Several non-synonymous mutations mapped in regulatory genes (arcB, barA, rpoD, rna) and in other putative regulatory loci (yjjU, rssA and ypdA). In addition, a chromosomal deletion of 10,328 bp was detected that removed 12 genes, among them, the rppH, mutH and galR genes. Characterization of some of these mutated and deleted genes with their functions and possible functions, are presented. CONCLUSIONS: The deletion of the contiguous rppH, mutH and galR genes that occurred simultaneously, is apparently the main reason for the faster growth of the evolved PB12 strain. In support of this interpretation is the fact that inactivation of the rppH gene in the parental PB11 strain substantially increased its growth rate, very likely by increasing glycolytic mRNA genes stability. Furthermore, galR inactivation allowed glucose transport by GalP into the cell. The deletion of mutH in an already stressed strain that lacks PTS is apparently responsible for the very high mutation rate observed.

An Expanded Genome-Wide Association Study of Fructosamine Levels Identifies RCN3 as a Replicating Locus and Implicates FCGRT as the Effector Transcript.

Fructosamine is a measure of short-term glycemic control, which has been suggested as a useful complement to glycated hemoglobin (HbA1c) for the diagnosis and monitoring of diabetes. To date, a single genome-wide association study (GWAS) including 8,951 U.S. White and 2,712 U.S. Black individuals without a diabetes diagnosis has been published. Results in Whites and Blacks yielded different association loci, near RCN3 and CNTN5, respectively. In this study, we performed a GWAS on 20,731 European-ancestry blood donors and meta-analyzed our results with previous data from U.S. White participants from the Atherosclerosis Risk in Communities (ARIC) study (Nmeta = 29,685). We identified a novel association near GCK (rs3757840, ßmeta = 0.0062; minor allele frequency [MAF] = 0.49; Pmeta = 3.66 × 10-8) and confirmed the association near RCN3 (rs113886122, ßmeta = 0.0134; MAF = 0.17; Pmeta = 5.71 × 10-18). Colocalization analysis with whole-blood expression quantitative trait loci data suggested FCGRT as the effector transcript at the RCN3 locus. We further showed that fructosamine has low heritability (h2 = 7.7%), has no significant genetic correlation with HbA1c and other glycemic traits in individuals without a diabetes diagnosis (P > 0.05), but has evidence of shared genetic etiology with some anthropometric traits (Bonferroni-corrected P < 0.0012). Our results broaden knowledge of the genetic architecture of fructosamine and prioritize FCGRT for downstream functional studies at the established RCN3 locus.

Transcriptome Analysis Reveals Cr(VI) Adaptation Mechanisms in Klebsiella sp. Strain AqSCr.

Klebsiella sp. strain AqSCr, isolated from Cr(VI)-polluted groundwater, reduces Cr(VI) both aerobically and anaerobically and resists up 34 mM Cr(VI); this resistance is independent of the ChrA efflux transporter. In this study, we report the whole genome sequence and the transcriptional profile by RNA-Seq of strain AqSCr under Cr(VI)-adapted conditions and found 255 upregulated and 240 downregulated genes compared to controls without Cr(VI) supplementation. Genes differentially transcribed were mostly associated with oxidative stress response, DNA repair and replication, sulfur starvation response, envelope-osmotic stress response, fatty acid (FA) metabolism, ribosomal subunits, and energy metabolism. Among them, genes not previously associated with chromium resistance, for example, cybB, encoding a putative superoxide oxidase (SOO), gltA2, encoding an alternative citrate synthase, and des, encoding a FA desaturase, were upregulated. The sodA gene encoding a manganese superoxide dismutase was upregulated in the presence of Cr(VI), whereas sodB encoding an iron superoxide dismutase was downregulated. Cr(VI) resistance mechanisms in strain AqSCr seem to be orchestrated by the alternative sigma factors fecl, rpoE, and rpoS (all of them upregulated). Membrane lipid analysis of the Cr(IV)-adapted strain showed a lower proportion of unsaturated lipids with respect to the control, which we hypothesized could result from unsaturated lipid peroxidation followed by degradation, together with de novo synthesis mediated by the upregulated FA desaturase-encoding gene, des. This report helps to elucidate both Cr(VI) toxicity targets and global bacterial response to Cr(VI).

Integrated Polygenic Tool Substantially Enhances Coronary Artery Disease Prediction.

BACKGROUND: There is considerable interest in whether genetic data can be used to improve standard cardiovascular disease risk calculators, as the latter are routinely used in clinical practice to manage preventative treatment. METHODS: Using the UK Biobank resource, we developed our own polygenic risk score for coronary artery disease (CAD). We used an additional 60 000 UK Biobank individuals to develop an integrated risk tool (IRT) that combined our polygenic risk score with established risk tools (either the American Heart Association/American College of Cardiology pooled cohort equations [PCE] or UK QRISK3), and we tested our IRT in an additional, independent set of 186 451 UK Biobank individuals. RESULTS: The novel CAD polygenic risk score shows superior predictive power for CAD events, compared with other published polygenic risk scores, and is largely uncorrelated with PCE and QRISK3. When combined with PCE into an IRT, it has superior predictive accuracy. Overall, 10.4% of incident CAD cases were misclassified as low risk by PCE and correctly classified as high risk by the IRT, compared with 4.4% misclassified by the IRT and correctly classified by PCE. The overall net reclassification improvement for the IRT was 5.9% (95% CI, 4.7-7.0). When individuals were stratified into age-by-sex subgroups, the improvement was larger for all subgroups (range, 8.3%-15.4%), with the best performance in 40- to 54-year-old men (15.4% [95% CI, 11.6-19.3]). Comparable results were found using a different risk tool (QRISK3) and also a broader definition of cardiovascular disease. Use of the IRT is estimated to avoid up to 12 000 deaths in the United States over a 5-year period. CONCLUSIONS: An IRT that includes polygenic risk outperforms current risk stratification tools and offers greater opportunity for early interventions. Given the plummeting costs of genetic tests, future iterations of CAD risk tools would be enhanced with the addition of a person's polygenic risk.

Validation of an Integrated Risk Tool, Including Polygenic Risk Score, for Atherosclerotic Cardiovascular Disease in Multiple Ethnicities and Ancestries.

The American College of Cardiology / American Heart Association pooled cohort equations tool (ASCVD-PCE) is currently recommended to assess 10-year risk for atherosclerotic cardiovascular disease (ASCVD). ASCVD-PCE does not currently include genetic risk factors. Polygenic risk scores (PRSs) have been shown to offer a powerful new approach to measuring genetic risk for common diseases, including ASCVD, and to enhance risk prediction when combined with ASCVD-PCE. Most work to date, including the assessment of tools, has focused on performance in individuals of European ancestries. Here we present evidence for the clinical validation of a new integrated risk tool (IRT), ASCVD-IRT, which combines ASCVD-PCE with PRS to predict 10-year risk of ASCVD across diverse ethnicity and ancestry groups. We demonstrate improved predictive performance of ASCVD-IRT over ASCVD-PCE, not only in individuals of self-reported White ethnicities (net reclassification improvement [NRI]; with 95% confidence interval = 2.7% [1.1 to 4.2]) but also Black / African American / Black Caribbean / Black African (NRI = 2.5% [0.6-4.3]) and South Asian (Indian, Bangladeshi or Pakistani) ethnicities (NRI = 8.7% [3.1 to 14.4]). NRI confidence intervals were wider and included zero for ethnicities with smaller sample sizes, including Hispanic (NRI = 7.5% [-1.4 to 16.5]), but PRS effect sizes in these ethnicities were significant and of comparable size to those seen in individuals of White ethnicities. Comparable results were obtained when individuals were analyzed by genetically inferred ancestry. Together, these results validate the performance of ASCVD-IRT in multiple ethnicities and ancestries, and favor their generalization to all ethnicities and ancestries.

A transcriptomic atlas of mammalian olfactory mucosae reveals an evolutionary influence on food odor detection in humans.

The mammalian olfactory system displays species-specific adaptations to different ecological niches. To investigate the evolutionary dynamics of olfactory sensory neuron (OSN) subtypes across mammalian evolution, we applied RNA sequencing of whole olfactory mucosa samples from mouse, rat, dog, marmoset, macaque, and human. We find that OSN subtypes, representative of all known mouse chemosensory receptor gene families, are present in all analyzed species. Further, we show that OSN subtypes expressing canonical olfactory receptors are distributed across a large dynamic range and that homologous subtypes can be either highly abundant across all species or species/order specific. Highly abundant mouse and human OSN subtypes detect odorants with similar sensory profiles and sense ecologically relevant odorants, such as mouse semiochemicals or human key food odorants. Together, our results allow for a better understanding of the evolution of mammalian olfaction in mammals and provide insights into the possible functions of highly abundant OSN subtypes.

Demographic history and biologically relevant genetic variation of Native Mexicans inferred from whole-genome sequencing.

Understanding the genetic structure of Native American populations is important to clarify their diversity, demographic history, and to identify genetic factors relevant for biomedical traits. Here, we show a demographic history reconstruction from 12 Native American whole genomes belonging to six distinct ethnic groups representing the three main described genetic clusters of Mexico (Northern, Southern, and Maya). Effective population size estimates of all Native American groups remained below 2,000 individuals for up to 10,000 years ago. The proportion of missense variants predicted as damaging is higher for undescribed (~ 30%) than for previously reported variants (~ 15%). Several variants previously associated with biological traits are highly frequent in the Native American genomes. These findings suggest that the demographic and adaptive processes that occurred in these groups shaped their genetic architecture and could have implications in biological processes of the Native Americans and Mestizos of today.

Draft Genome Sequence of Leuconostoc mesenteroides P45 Isolated from Pulque, a Traditional Mexican Alcoholic Fermented Beverage.

Leuconostoc mesenteroides P45 was isolated from the traditional Mexican pulque beverage. We report its draft genome sequence, assembled in 6 contigs consisting of 1,874,188 bp and no plasmids. Genome annotation predicted a total of 1,800 genes, 1,687 coding sequences, 52 pseudogenes, 9 rRNAs, 51 tRNAs, 1 noncoding RNA, and 44 frameshifted genes.