Echocardiographic heart ageing patterns predict cardiovascular and non-cardiovascular events and reflect biological age: the SardiNIA study

AIMS

Age is a crucial risk factor for cardiovascular (CV) and non-CV diseases. As people age at different rates, the concept of biological age has been introduced as a personalized measure of functional deterioration. Associations of age with echocardiographic quantitative traits were analysed to assess different heart ageing rates and their ability to predict outcomes and reflect biological age.

METHODS AND RESULTS

Associations of age with left ventricular mass, geometry, diastolic function, left atrial volume, and aortic root size were measured in 2614 healthy subjects. Based on the 95% two-sided tolerance intervals of each correlation, three discrete ageing trajectories were identified and categorized as 'slow', 'normal', and 'accelerated' heart ageing patterns. The primary endpoint included fatal and non-fatal CV events, and the secondary endpoint was a composite of CV and non-CV events and all-cause death. The phenotypic age of the heart (HeartPhAge) was estimated as a proxy of biological age. The slow ageing pattern was found in 8.7% of healthy participants, the normal pattern in 76.9%, and the accelerated pattern in 14.4%. Kaplan-Meier curves of the heart ageing patterns diverged significantly (P = 0.0001) for both primary and secondary endpoints, with the event rate being lowest in the slow, intermediate in the normal, and highest in the accelerated pattern. In the Cox proportional hazards model, heart ageing patterns predicted both primary (P = 0.01) and secondary (P = 0.03 to <0.0001) endpoints, independent of chronological age and risk factors. Compared with chronological age, HeartPhAge was 9 years younger in slow, 4 years older in accelerated (both P < 0.0001), and overlapping in normal ageing patterns.

CONCLUSION

Standard Doppler echocardiography detects slow, normal, and accelerated heart ageing patterns. They predict CV and non-CV events, reflect biological age, and provide a new tool to calibrate prevention timing and intensity.

Transcriptome organization of white blood cells through gene co-expression network analysis in a large RNA-seq dataset

Gene co-expression network analysis enables identification of biologically meaningful clusters of co-regulated genes (modules) in an unsupervised manner. We present here the largest study conducted thus far of co-expression networks in white blood cells (WBC) based on RNA-seq data from 624 individuals. We identify 41 modules, 13 of them related to specific immune-related functions and cell types (e.g. neutrophils, B and T cells, NK cells, and plasmacytoid dendritic cells); we highlight biologically relevant lncRNAs for each annotated module of co-expressed genes. We further characterize with unprecedented resolution the modules in T cell sub-types, through the availability of 95 immune phenotypes obtained by flow cytometry in the same individuals. This study provides novel insights into the transcriptional architecture of human leukocytes, showing how network analysis can advance our understanding of coding and non-coding gene interactions in immune system cells.

Multi-trait analysis characterizes the genetics of thyroid function and identifies causal associations with clinical implications

To date only a fraction of the genetic footprint of thyroid function has been clarified. We report a genome-wide association study meta-analysis of thyroid function in up to 271,040 individuals of European ancestry, including reference range thyrotropin (TSH), free thyroxine (FT4), free and total triiodothyronine (T3), proxies for metabolism (T3/FT4 ratio) as well as dichotomized high and low TSH levels. We revealed 259 independent significant associations for TSH (61% novel), 85 for FT4 (67% novel), and 62 novel signals for the T3 related traits. The loci explained 14.1%, 6.0%, 9.5% and 1.1% of the total variation in TSH, FT4, total T3 and free T3 concentrations, respectively. Genetic correlations indicate that TSH associated loci reflect the thyroid function determined by free T3, whereas the FT4 associations represent the thyroid hormone metabolism. Polygenic risk score and Mendelian randomization analyses showed the effects of genetically determined variation in thyroid function on various clinical outcomes, including cardiovascular risk factors and diseases, autoimmune diseases, and cancer. In conclusion, our results improve the understanding of thyroid hormone physiology and highlight the pleiotropic effects of thyroid function on various diseases.

Genetic and biomarker modulation of arterial stiffness change in the SardiNIA population cohort

Background and aims

Arterial stiffness (AS), quantified by pulse wave velocity (PWV), arises due to impaired arterial elastic tissue and smooth muscle dysfunction. We aimed to examine the longitudinal association of genetic, lipid and inflammation biomarkers with PWV and how these associations may change with aging.

Materials and methods

We utilized genotype and four time-point biomarker data from the SardiNIA cohort [n = 6,301; mean baseline age 43.3 (SD 17.3); 58% females]. To investigate the association of PWV with genetic variants, lipid, and inflammation biomarkers, we employed linear mixed modeling, using age as the time scale. Biomarkers exhibiting significant longitudinal associations were categorized into tertiles and individuals within the second tertile or those with heterozygous alleles were excluded, leaving a cohort of 2,000 individuals. This cohort was further divided into four risk groups: low genetic and low biomarker (L-L), low genetic and high biomarker (L-H), high genetic and low biomarker (H-L), and high genetic and high biomarker risk (H-H). Subsequent analyses focused on these risk groups to assess their association to PWV with time.

Results

Using the complete dataset, we found a significant longitudinal association of total cholesterol (TC), triglycerides (TG), fibrinogen (FGN), and total white blood cell count (TWBC) with PWV, all with p < 3.33 × 10-3. After grouping, individuals with homogeneous risk alleles of SNP rs3742207 and high baseline TG levels (H-H group) exhibited a 1.39-fold higher PWV (m/s) (95% CI, 1.17-1.64, p = 1.21 × 10-4) compared to the L-L group. Similarly, individuals in the H-H group of rs3742207-TWBC combination showed 1.75 times higher PWV (95% CI, 1.48-0.2.07, p = 1.01 × 10-10) compared to the L-L group. Similar patterns were observed for groups based on SNP rs7152623-TWBC risk. Furthermore, these associations became more pronounced with increasing age (p < 3.33 × 10-3).

Conclusion

The longitudinal association of TG and TWBC biomarkers with PWV varied by SNPs rs3742207 and rs7152623 genotype. Further studies are warranted to investigate the function of genetics, lipids, and inflammation biomarkers on PWV change.

Glucose-6-phosphate dehydrogenase deficiency accelerates arterial aging in diabetes

AIMS

High glucose levels and Glucose-6-Phosphate Dehydrogenase deficiency (G6PDd) have both tissue inflammatory effects. Here we determined whether G6PDd accelerates arterial aging (information linked stiffening) in diabetes.

METHODS

Plasma glucose, interleukin 6 (IL6), and arterial stiffness (indexed as carotid-femoral Pulse Wave Velocity, PWV) and red blood cell G6PD activity were assessed in a large (4448) Sardinian population.

RESULTS

Although high plasma glucose in diabetics, did not differ by G6DP status (178.2 ± 55.1 vs 169.0 ± 50.1 mg/dl) in G6DPd versus non-G6PDd subjects, respectively, IL6, and PWV (adjusted for age and glucose) were significantly increased in G6PDd vs non-G6PDd subjects (PWV, 8.0 ± 0.4 vs 7.2 ± 0.2 m/sec) and (IL6, 6.9 ± 5.0 vs 4.2 ± 3.0 pg/ml). In non-diabetics, neither fasting plasma glucose, nor IL6, nor PWV were impacted by G6PDd.

CONCLUSION

G6PDd in diabetics is associated with increased inflammatory markers and accelerated arterial aging.

Understanding the genetic complexity of puberty timing across the allele frequency spectrum

Pubertal timing varies considerably and has been associated with a range of health outcomes in later life. To elucidate the underlying biological mechanisms, we performed multi-ancestry genetic analyses in ~800,000 women, identifying 1,080 independent signals associated with age at menarche. Collectively these loci explained 11% of the trait variance in an independent sample, with women at the top and bottom 1% of polygenic risk exhibiting a ~11 and ~14-fold higher risk of delayed and precocious pubertal development, respectively. These common variant analyses were supported by exome sequence analysis of ~220,000 women, identifying several genes, including rare loss of function variants in ZNF483 which abolished the impact of polygenic risk. Next, we implicated 660 genes in pubertal development using a combination of in silico variant-to-gene mapping approaches and integration with dynamic gene expression data from mouse embryonic GnRH neurons. This included an uncharacterized G-protein coupled receptor GPR83, which we demonstrate amplifies signaling of MC3R, a key sensor of nutritional status. Finally, we identified several genes, including ovary-expressed genes involved in DNA damage response that co-localize with signals associated with menopause timing, leading us to hypothesize that the ovarian reserve might signal centrally to trigger puberty. Collectively these findings extend our understanding of the biological complexity of puberty timing and highlight body size dependent and independent mechanisms that potentially link reproductive timing to later life disease.

Immune resilience despite inflammatory stress promotes longevity and favorable health outcomes including resistance to infection

Some people remain healthier throughout life than others but the underlying reasons are poorly understood. Here we hypothesize this advantage is attributable in part to optimal immune resilience (IR), defined as the capacity to preserve and/or rapidly restore immune functions that promote disease resistance (immunocompetence) and control inflammation in infectious diseases as well as other causes of inflammatory stress. We gauge IR levels with two distinct peripheral blood metrics that quantify the balance between (i) CD8+ and CD4+ T-cell levels and (ii) gene expression signatures tracking longevity-associated immunocompetence and mortality-associated inflammation. Profiles of IR metrics in ~48,500 individuals collectively indicate that some persons resist degradation of IR both during aging and when challenged with varied inflammatory stressors. With this resistance, preservation of optimal IR tracked (i) a lower risk of HIV acquisition, AIDS development, symptomatic influenza infection, and recurrent skin cancer; (ii) survival during COVID-19 and sepsis; and (iii) longevity. IR degradation is potentially reversible by decreasing inflammatory stress. Overall, we show that optimal IR is a trait observed across the age spectrum, more common in females, and aligned with a specific immunocompetence-inflammation balance linked to favorable immunity-dependent health outcomes. IR metrics and mechanisms have utility both as biomarkers for measuring immune health and for improving health outcomes.

Effect of Genetic Factors, Age and Sex on Levels of Circulating Extracellular Vesicles and Platelets

Extracellular vesicles (EVs) mediate cell interactions in biological processes, such as receptor activation or molecule transfer. Estimates of variation by age and sex have been limited by small sample size, and no report has assessed the contribution of genetic factors to levels of EVs. Here, we evaluated blood levels of 25 EV and 3 platelet traits in 974 individuals (933 genotyped) and reported the first genome-wide association study (GWAS) on levels of these traits. EV levels all decreased with age, whereas the trend for their surface markers was more heterogeneous. Platelets and CD31dim platelet EVs significantly increased in females compared to males, although CD31 expression on both platelets and platelet EVs decreased in females. Levels of the other EV subsets were similar between sexes. GWAS revealed three statistically significant genetic signals associated with EV levels in the F10 and GBP1 genes and in the intergenic region between LRIG1 and KBTBD8. These add to a signal in the 3'UTR of RHOF associated with CD31 expression on platelets that was previously found to be associated with other platelet traits. These findings suggest that EV formation is not a simple, constant adjunct of metabolism but is under both age-related and genetic control that can be independent of the regulation of the levels of the cells from which the EVs derive.

A fast linkage method for population GWAS cohorts with related individuals

Linkage analysis, a class of methods for detecting co-segregation of genomic segments and traits in families, was used to map disease-causing genes for decades before genotyping arrays and dense SNP genotyping enabled genome-wide association studies in population samples. Population samples often contain related individuals, but the segregation of alleles within families is rarely used because traditional linkage methods are computationally inefficient for larger datasets. Here, we describe Population Linkage, a novel application of Haseman-Elston regression as a method of moments estimator of variance components and their standard errors. We achieve additional computational efficiency by using modern methods for detection of IBD segments and variance component estimation, efficient preprocessing of input data, and minimizing redundant numerical calculations. We also refined variance component models to account for the biases in population-scale methods for IBD segment detection. We ran Population Linkage on four blood lipid traits in over 70,000 individuals from the HUNT and SardiNIA studies, successfully detecting 25 known genetic signals. One notable linkage signal that appeared in both was for low-density lipoprotein (LDL) cholesterol levels in the region near the gene APOE (LOD = 29.3, variance explained = 4.1%). This is the region where the missense variants rs7412 and rs429358, which together make up the ε2, ε3, and ε4 alleles each account for 2.4% and 0.8% of variation in circulating LDL cholesterol. Our results show the potential for linkage analysis and other large-scale applications of method of moments variance components estimation.

Peripheral blood mononuclear cells reactivity in recent-onset type I diabetes patients is directed against the leader peptide of preproinsulin, GAD65271-285 and GAD65431-450

Introduction

T cell reactivity against pancreatic autoantigens is considered one of the main contributors to the destruction of insulin-producing cells in type 1 diabetes (T1D). Over the years, peptide epitopes derived from these autoantigens have been described in NOD mice and in both HLA class II transgenic mice and humans. However, which ones are involved in the early onset or in the progressive phases of the disease is still unclear.

Methods

In this work we have investigated, in early-onset T1D pediatric patients and HLA-matched controls from Sardinia, the potential of preproinsulin (PPI) and glutamate decarboxylase 65 (GAD65)-derived peptides to induce spontaneous T cell proliferation responses of peripheral blood mononuclear cells (PBMCs).

Results

Significant T cell responses against PPI1-18, PPI7-19 and PPI31-49, the first two belonging to the leader sequence of PPI, and GAD65271-285 and GAD65431-450, were found in HLA-DR4, -DQ8 and -DR3, -DQ2 T1D children.

Conclusions

These data show that cryptic epitopes from the leader sequence of the PPI and GAD65271-285 and GAD65431-450 peptides might be among the critical antigenic epitopes eliciting the primary autoreactive responses in the early phases of the disease. These results may have implications in the design of immunogenic PPI and GAD65 peptides for peptide-based immunotherapy.

Polygenic risk score and biochemical/environmental variables predict a low-risk profile of age-related macular degeneration in Sardinia

PURPOSE

To ascertain the prevalence and clinical and genetic features of age-related macular degeneration (AMD) in subjects living in the Lanusei valley, Central Sardinia, Italy, involved in a study on ageing (SardiNIA project).

METHODS

A total of 814 volunteers aged ≥ 50 years, randomly selected from the SardiNIA project dataset, were included. A color fundus (CF) photograph of the 30° central retina of each eye was obtained and graded according to the Age-Related Eye Disease Study system. Life-style choices were investigated using standardized questionnaires. The concentrations of several inflammatory biomarkers (i.e., complement component, fibrinogen, and C-reactive protein) were measured. Polygenic risk score (PRS) was calculated and compared with results obtained from a European cohort.

RESULTS

A total of 756 subjects had gradable CF photographs for AMD detection. In 91.3%, no signs of AMD were observed. The prevalence rates of early and late AMDs were 6.9% and 0.6%, respectively. A total of 85% of subjects were physically active; only 13.5% were current smokers. Low concentrations of complement component, fibrinogen, and C-reactive protein were found. We calculated the polygenic risk scores (PRS) using 40 AMD markers distributed on several candidate genes in Europeans and Sardinians. The mean PRS value was significantly lower in Sardinians than in the Europeans (0.21 vs. 0.248, respectively, p = 1.18 × 10-77).

CONCLUSIONS

In our cohort, most subjects showed no sign of any AMD type and late AMD was a condition rarely observed. Results of genetic, biochemical, and life-style investigation support the hypothesis that Sardinia population may present of a peculiar background with a protective effect against AMD development.

Patterns of Dietary Blood Markers Are Related to Frailty Status in the FRAILOMIC Validation Phase

The influence of nutritional factors on frailty syndrome is still poorly understood. Thus, we aimed to confirm cross-sectional associations of diet-related blood biomarker patterns with frailty and pre-frailty statuses in 1271 older adults from four European cohorts. Principal component analysis (PCA) was performed based on plasma levels of α-carotene, β-carotene, lycopene, lutein + zeaxanthin, β-cryptoxanthin, α-tocopherol, γ-tocopherol and retinol. Cross-sectional associations between biomarker patterns and frailty status, according to Fried's frailty criteria, were assessed by using general linear models and multinomial logistic regression models as appropriate with adjustments for the main potential confounders. Robust subjects had higher concentrations of total carotenoids, β-carotene and β-cryptoxanthin than frail and pre-frail subjects and had higher lutein + zeaxanthin concentrations than frail subjects. No associations between 25-Hydroxyvitamin D3 and frailty status were observed. Two distinct biomarker patterns were identified in the PCA results. The principal component 1 (PC1) pattern was characterized by overall higher plasma levels of carotenoids, tocopherols and retinol, and the PC2 pattern was characterized by higher loadings for tocopherols, retinol and lycopene together and lower loadings for other carotenoids. Analyses revealed inverse associations between PC1 and prevalent frailty. Compared to participants in the lowest quartile of PC1, those in the highest quartile were less likely to be frail (odds ratio: 0.45, 95% CI: 0.25-0.80, p = 0.006). In addition, those in the highest quartile of PC2 showed higher odds for prevalent frailty (2.48, 1.28-4.80, p = 0.007) than those in the lowest quartile. Our findings strengthen the results from the first phase of the FRAILOMIC project, indicating carotenoids are suitable components for future biomarker-based frailty indices.

GWAS of genetic factors affecting white blood cell morphological parameters in Sardinians uncovers influence of chromosome 11 innate immunity gene cluster on eosinophil morphology

Few genome-wide association studies (GWAS) analyzing genetic regulation of morphological traits of white blood cells have been reported. We carried out a GWAS of 12 morphological traits in 869 individuals from the general population of Sardinia, Italy. These traits, included measures of cell volume, conductivity and light scatter in four white-cell populations (eosinophils, lymphocytes, monocytes, neutrophils). This analysis yielded seven statistically significant signals, four of which were novel (four novel, PRG2, P2RX3, two of CDK6). Five signals were replicated in the independent INTERVAL cohort of 11 822 individuals. The most interesting signal with large effect size on eosinophil scatter (P-value = 8.33 x 10-32, beta = -1.651, se = 0.1351) falls within the innate immunity cluster on chromosome 11, and is located in the PRG2 gene. Computational analyses revealed that a rare, Sardinian-specific PRG2:p.Ser148Pro mutation modifies PRG2 amino acid contacts and protein dynamics in a manner that could possibly explain the changes observed in eosinophil morphology. Our discoveries shed light on genetics of morphological traits. For the first time, we describe such large effect size on eosinophils morphology that is relatively frequent in Sardinian population.

708 PATTERNS OF BIOLOGICAL AGING OF THE HEART AND THEIR PROGNOSTIC VALUE: THE SARDI-NIA STUDY

Background

As people age at different rates, the concept of biological age has been introduced as a measure of functional deterioration. Associations of echocardiographic traits with age were analyzed to detect different biological aging rates of the heart and their prognostic value.

Methods

Echocardiographic left ventricular mass, geometry, left atrial volume, diastolic function, and aortic root size were measured in 3817 adults (58% females). Sex-specific correlations of age with six cardiac traits were assessed in 2614 healthy subjects. According to two-sided 95% tolerance intervals, three patterns of heart aging rates were identified in the whole cohort: slow aging pattern (n=287), normal aging pattern (n=2669), and accelerated aging pattern (n=796). The phenotypic (biological) age of the heart (HeartPhAge) was estimated from the cardiac traits by multiple linear regression. Cardiovascular and non-cardiovascular incident events at follow-up (46 months (IQR 34–64 months) were grouped into three endpoints and analyzed in survival analyses.

Results

Compared to chronological age, HeartPhAge was nine years younger (p&lt;0.0001) in the slow aging pattern, similar in the normal pattern, and three years older (p&lt;0.0001) in the accelerated pattern. The incidence of all endpoints was higher in the accelerated pattern than in the normal and slow aging patterns (p&lt;0.0001 to p=0.0003). In Cox proportional-hazards analysis, the heart aging patterns significantly predicted all endpoints in both the whole cohort and the healthy subset, independent of chronological age and risk factors.

Conclusion

Heart aging patterns reflect the rate of biological aging and independently predict cardiovascular and non-cardiovascular events.

Dilated hypertrophic phenotype of the carotid artery is associated with accelerated age-associated central arterial stiffening

Hypertrophic carotid geometric phenotypes (h-CGP) are predictors of incident cardiovascular disease (CVD). While arterial aging is hypothesized as a contributor to this associated risk, the association of CGPs with chronological age is not clear. In this manuscript we examine whether hypertrophic CGPs represent accelerated biological, rather than chronological, aging by examining their association with carotid-femoral pulse wave velocity (PWV), the hallmark of arterial aging. We analyzed data from 5516 participants of the SardiNIA study with a wide range of age at baseline (20-101 years), and a median follow-up time of 13 years (mean 11.5 years; maximum 17.9 years). Baseline CGPs were defined based on the common carotid lumen diameter, wall thickness, and their ratio. Subject-specific rates of change of PWV, blood pressure parameters, body mass index, glucose, and lipids were estimated using linear mixed effects models. Compared to those with typical(t-) CGP, those with dilated hypertrophy (dh-) CGP had a greater longitudinal increase in PWV; this increase was significantly greater among older individuals and men. The greater PWV longitudinal increase in dh-CGP remained significant after adjusting for baseline values and rates of change of covariates. Dilated hypertrophic CGP is independently associated with accelerated increase in age-associated arterial stiffening over time, with a strong association in men than in women. Future studies are needed to examine if this association mediates the increased risk for CVD observed in individuals with hypertrophic cardiac remodelling and the role of retarding it to reduce this risk. HIGHLIGHTS: • Individuals with dilated hypertrophic geometric phenotypes of the common carotid artery (increased age- and sex-specific wall thickness and lumen diameter) have greater future central arterial stiffening, independently of other determinants of arterial stiffening. • The dilated hypertrophic phenotype group has a greater age-specific arterial dilation, wall thickening, and stiffness (the arterial aging triad). This suggests that this phenotype is a form of accelerated aging that might explain the worse clinic outcomes observed in this group. • Understanding the natural history of the carotid geometric phenotype across the lifespan and the determinants of the deleterious progression towards the dilated hypertrophic phenotype are needed to develop interventions that reduce the adverse clinical outcomes associated with it.

Genetic diversity fuels gene discovery for tobacco and alcohol use

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 injury1-4. These substances are used across the globe, yet genome-wide association studies have focused largely on individuals of European ancestries5. 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.

A saturated map of common genetic variants associated with human height

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 sizes1. 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 panel2) 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.

Genetic analyses of the electrocardiographic QT interval and its components identify additional loci and pathways

The QT interval is an electrocardiographic measure representing the sum of ventricular depolarization and repolarization, estimated by QRS duration and JT interval, respectively. QT interval abnormalities are associated with potentially fatal ventricular arrhythmia. Using genome-wide multi-ancestry analyses (>250,000 individuals) we identify 177, 156 and 121 independent loci for QT, JT and QRS, respectively, including a male-specific X-chromosome locus. Using gene-based rare-variant methods, we identify associations with Mendelian disease genes. Enrichments are observed in established pathways for QT and JT, and previously unreported genes indicated in insulin-receptor signalling and cardiac energy metabolism. In contrast for QRS, connective tissue components and processes for cell growth and extracellular matrix interactions are significantly enriched. We demonstrate polygenic risk score associations with atrial fibrillation, conduction disease and sudden cardiac death. Prioritization of druggable genes highlight potential therapeutic targets for arrhythmia. Together, these results substantially advance our understanding of the genetic architecture of ventricular depolarization and repolarization.

Cross-sectional analysis of the humoral response after SARS-CoV-2 vaccination in Sardinian multiple sclerosis patients, a follow-up study

Monitoring immune responses to SARS-CoV-2 vaccination and its clinical efficacy over time in Multiple Sclerosis (MS) patients treated with disease-modifying therapies (DMTs) help to establish the optimal strategies to ensure adequate COVID-19 protection without compromising disease control offered by DMTs. Following our previous observations on the humoral response one month after two doses of BNT162b2 vaccine (T1) in MS patients differently treated, here we present a cross-sectional and longitudinal follow-up analysis six months following vaccination (T2, n=662) and one month following the first booster (T3, n=185). Consistent with results at T1, humoral responses were decreased in MS patients treated with fingolimod and anti-CD20 therapies compared with untreated patients also at the time points considered here (T2 and T3). Interestingly, a strong upregulation one month after the booster was observed in patients under every DMTs analyzed, including those treated with fingolimod and anti-CD20 therapies. Although patients taking these latter therapies had a higher rate of COVID-19 infection five months after the first booster, only mild symptoms that did not require hospitalization were reported for all the DMTs analyzed here. Based on these findings we anticipate that additional vaccine booster shots will likely further improve immune responses and COVID-19 protection in MS patients treated with any DMT.

Facets of Personality and Risk of Cognitive Impairment: Longitudinal Findings in a Rural Community from Sardinia

Background: Few studies have examined the associations between personality facets and dementia risk and rarely included individuals from rural settings or with low education. Objective: To examine the association between personality and the risk of cognitive impairment. Methods: Participants (N = 1,668; age 50 to 94 at baseline; 56.4% women; 86.5% less than high school diploma) were from a rural region of Sardinia (Italy) who completed the Revised NEO Personality Inventory (NEO-PI-R) during the first wave (2001–2004) and the Mini-Mental State Examination (MMSE) at waves two to five (2005–2021). Cox regression was used to test personality and covariates as predictors of cognitive impairment based on MMSE education-adjusted cutoffs. Results: During the up to 18-year follow-up (M = 10.38; SD = 4.76), 187 individuals (11.2%) scored as cognitively impaired. Participants with higher neuroticism (particularly the depression facet [HR = 1.22, 95% CI = 1.06–1.40]), and lower agreeableness (particularly the modesty facet [HR = 0.83, 95% CI = 0.71–0.97]) and lower conscientiousness (particularly the dutifulness facet [HR = 0.78, 95% CI = 0.67–0.92]) were at higher risk of cognitive impairment. Lower warmth ([HR = 0.75, 95% CI = 0.65–0.87], facet of extraversion) and ideas ([HR = 0.76, 95% CI = 0.65–0.89], facet of openness) were also associated with increased risk of impairment. These associations were virtually unchanged in models that accounted for other risk factors, including smoking, depression, obesity, hypertension, diabetes, and apolipoprotein E (APOE) ɛ4 carrier status. Across the five domains, sex and the APOE variant did not moderate the associations. Conclusion: In a sample with demographic characteristics underrepresented in dementia research, this study identifies personality domains and facets most relevant to the risk of cognitive impairment.

Autoimmunity Is a Significant Feature of Idiopathic Pulmonary Arterial Hypertension

Rationale: Autoimmunity is believed to play a role in idiopathic pulmonary arterial hypertension (IPAH). It is not clear whether this is causative or a bystander of disease and if it carries any prognostic or treatment significance. Objectives: To study autoimmunity in IPAH using a large cross-sectional cohort. Methods: Assessment of the circulating immune cell phenotype was undertaken using flow cytometry, and the profile of serum immunoglobulins was generated using a standardized multiplex array of 19 clinically validated autoantibodies in 473 cases and 946 control subjects. Additional glutathione S-transferase fusion array and ELISA data were used to identify a serum autoantibody to BMPR2 (bone morphogenetic protein receptor type 2). Clustering analyses and clinical correlations were used to determine associations between immunogenicity and clinical outcomes. Measurements and Main Results: Flow cytometric immune profiling demonstrates that IPAH is associated with an altered humoral immune response in addition to raised IgG3. Multiplexed autoantibodies were significantly raised in IPAH, and clustering demonstrated three distinct clusters: "high autoantibody," "low autoantibody," and a small "intermediate" cluster exhibiting high concentrations of ribonucleic protein complex. The high-autoantibody cluster had worse hemodynamics but improved survival. A small subset of patients demonstrated immunoglobulin reactivity to BMPR2. Conclusions: This study establishes aberrant immune regulation and presence of autoantibodies as key features in the profile of a significant proportion of patients with IPAH and is associated with clinical outcomes.

The longitudinal dynamics and natural history of clonal haematopoiesis

Clonal expansions driven by somatic mutations become pervasive across human tissues with age, including in the haematopoietic system, where the phenomenon is termed clonal haematopoiesis1-4. The understanding of how and when clonal haematopoiesis develops, the factors that govern its behaviour, how it interacts with ageing and how these variables relate to malignant progression remains limited5,6. Here we track 697 clonal haematopoiesis clones from 385 individuals 55 years of age or older over a median of 13 years. We find that 92.4% of clones expanded at a stable exponential rate over the study period, with different mutations driving substantially different growth rates, ranging from 5% (DNMT3A and TP53) to more than 50% per year (SRSF2P95H). Growth rates of clones with the same mutation differed by approximately ±5% per year, proportionately affecting slow drivers more substantially. By combining our time-series data with phylogenetic analysis of 1,731 whole-genome sequences of haematopoietic colonies from 7 individuals from an older age group, we reveal distinct patterns of lifelong clonal behaviour. DNMT3A-mutant clones preferentially expanded early in life and displayed slower growth in old age, in the context of an increasingly competitive oligoclonal landscape. By contrast, splicing gene mutations drove expansion only later in life, whereas TET2-mutant clones emerged across all ages. Finally, we show that mutations driving faster clonal growth carry a higher risk of malignant progression. Our findings characterize the lifelong natural history of clonal haematopoiesis and give fundamental insights into the interactions between somatic mutation, ageing and clonal selection.

Comparison of Whole Blood Cryopreservation Methods for Extensive Flow Cytometry Immunophenotyping

Fresh blood immunophenotyping by flow cytometry, based on the reliable simultaneous detection of several markers in a cell, is the method of choice to study the circulating human immune system. Especially in large and multicenter studies, high sample quality is difficult to achieve, and adequate collection and storage of samples with fine-tuned whole blood cryopreservation is mandatory. Here, we compared the quality of immunophenotypic data obtained from fresh blood with those obtained after five cryopreservation methods by quantifying the levels of 41 immune cell populations. They comprised B and T lymphocyte subsets and their maturation stages, as well as monocytes and granulocytes. Three methods used fixative solutions and two other methods used dimethyl sulfoxide solutions to preserve cell viability. The fixative methods prevented detection of markers critical for identification of B and T cell subsets, including CD27, CXCR3, and CCR6. The other two methods permitted reliable discrimination of most immune-cell populations in thawed samples, though some cell frequencies varied compared to the corresponding fresh sample. Of those two methods, the one preserving blood in media containing dimethyl sulfoxide produced results that were most similar to those with fresh samples.

Non-coding RNAs in malaria infection

Malaria is one of the most severe infectious diseases affecting humans and it is caused by protozoan pathogens of the species Plasmodium (spp.). The malaria parasite Plasmodium is characterized by a complex, multistage life cycle that requires tight gene regulation which allows for host invasion and defense against host immune responses. Unfortunately, the mechanisms regulating gene expression during Plasmodium infection remain largely elusive, though several lines of evidence implicate a major involvement of non-coding RNAs (ncRNAs). The ncRNAs have been found to play a key role in regulating transcriptional and post-transcriptional events in a broad range of organisms including Plasmodium. In Plasmodium ncRNAs have been shown to regulate key events in the multistage life cycle and virulence ability. Here we review recent progress involving ncRNAs (microRNAs, long non-coding RNAs, and circular RNAs) and their role as regulators of gene expression during Plasmodium infection in human hosts with focus on the possibility of using these molecules as biomarkers for monitoring disease status. We also discuss the surprising function of ncRNAs in mediating the complex interplay between parasite and human host and future perspectives of the field. This article is categorized under: RNA in Disease and Development > RNA in Disease.

Systematic identification of NF90 target RNAs by iCLIP analysis

RNA-binding proteins (RBPs) interact with and determine the fate of many cellular RNAs directing numerous essential roles in cellular physiology. Nuclear Factor 90 (NF90) is an RBP encoded by the interleukin enhancer-binding factor 3 (ILF3) gene that has been found to influence RNA metabolism at several levels, including pre-RNA splicing, mRNA turnover, and translation. To systematically identify the RNAs that interact with NF90, we carried out iCLIP (individual-nucleotide resolution UV crosslinking and immunoprecipitation) analysis in the human embryonic fibroblast cell line HEK-293. Interestingly, many of the identified RNAs encoded proteins involved in the response to viral infection and RNA metabolism. We validated a subset of targets and investigated the impact of NF90 on their expression levels. Two of the top targets, IRF3 and IRF9 mRNAs, encode the proteins IRF3 and IRF9, crucial regulators of the interferon pathway involved in the SARS-CoV-2 immune response. Our results support a role for NF90 in modulating key genes implicated in the immune response and offer insight into the immunological response to the SARS-CoV-2 infection.

Personality traits are consistently associated with blood mitochondrial DNA copy number estimated from genome sequences in two genetic cohort studies

Background

Mitochondrial DNA copy number (mtDNAcn) in tissues and blood can be altered in conditions like diabetes and major depression and may play a role in aging and longevity. However, little is known about the association between mtDNAcn and personality traits linked to emotional states, metabolic health, and longevity. This study tests the hypothesis that blood mtDNAcn is related to personality traits and mediates the association between personality and mortality.

Methods

We assessed the big five personality domains and facets using the Revised NEO Personality Inventory (NEO-PI-R), assessed depressive symptoms with the Center for Epidemiologic Studies Depression Scale (CES-D), estimated mtDNAcn levels from whole-genome sequencing, and tracked mortality in participants from the Baltimore Longitudinal Study of Aging. Results were replicated in the SardiNIA Project.

Results

We found that mtDNAcn was negatively associated with the Neuroticism domain and its facets and positively associated with facets from the other four domains. The direction and size of the effects were replicated in the SardiNIA cohort and were robust to adjustment for potential confounders in both samples. Consistent with the Neuroticism finding, higher depressive symptoms were associated with lower mtDNAcn. Finally, mtDNAcn mediated the association between personality and mortality risk.

Conclusions

To our knowledge, this is the first study to show a replicable association between mtDNAcn and personality. Furthermore, the results support our hypothesis that mtDNAcn is a biomarker of the biological process that explains part of the association between personality and mortality.

Funding

Support for this work was provided by the Intramural Research Program of the National Institute on Aging (Z01-AG000693, Z01-AG000970, and Z01-AG000949) and the National Institute of Neurological Disorders and Stroke of the National Institutes of Health. AT was also supported by the National Institute on Aging of the National Institutes of Health Grant R01AG068093.

Effect of Different Disease-Modifying Therapies on Humoral Response to BNT162b2 Vaccine in Sardinian Multiple Sclerosis Patients

Objectives

Vaccination against COVID-19 is highly recommended to patients affected by multiple sclerosis (MS); however, the impact of MS disease-modifying therapies (DMTs) on the immune response following vaccination has been only partially investigated. Here, we aimed to elucidate the effect of DMTs on the humoral immune response to mRNA-based anti-SARS-CoV-2 vaccines in MS patients.

Methods

We obtained sera from 912 Sardinian MS patients and 63 healthy controls 30 days after the second dose of BNT162b2 vaccine and tested them for SARS-CoV-2 response using anti-Spike (S) protein-based serology. Previous SARS-CoV-2 infection was assessed by anti-Nucleocapsid (N) serology. Patients were either untreated or undergoing treatment with a total of 13 different DMTs. Differences between treatment groups comprised of at least 10 patients were assessed by generalized linear mixed-effects model. Demographic and clinical data and smoking status were analyzed as additional factors potentially influencing humoral immunity from COVID-19 vaccine.

Results

MS patients treated with natalizumab, teriflunomide, azathioprine, fingolimod, ocrelizumab, and rituximab showed significantly lower humoral responses compared to untreated patients. We did not observe a statistically significant difference in response between patients treated with the other drugs (dimethyl fumarate, interferon, alemtuzumab and glatiramer acetate) and untreated patients. In addition, older age, male sex and active smoking were significantly associated with lower antibody titers against SARS-CoV-2. MS patients previously infected with SARS-CoV-2 had significantly higher humoral responses to vaccine than uninfected patients.

Conclusion

Humoral response to BNT162b2 is significantly influenced by the specific DMTs followed by patients, as well as by other factors such as previous SARS-CoV-2 infection, age, sex, and smoking status. These results are important to inform targeted strategies to prevent clinically relevant COVID-19 in MS patients.

The power of genetic diversity in genome-wide association studies of lipids

Increased blood lipid levels are heritable risk factors of cardiovascular disease with varied prevalence worldwide owing to different dietary patterns and medication use1. Despite advances in prevention and treatment, in particular through reducing low-density lipoprotein cholesterol levels2, heart disease remains the leading cause of death worldwide3. 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 GWAS4-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 patterns24. 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 medicine25, we anticipate that increased diversity of participants will lead to more accurate and equitable26 application of polygenic scores in clinical practice.

Immunologic resilience and COVID-19 survival advantage

BACKGROUND

The risk of severe coronavirus disease 2019 (COVID-19) varies significantly among persons of similar age and is higher in males. Age-independent, sex-biased differences in susceptibility to severe COVID-19 may be ascribable to deficits in a sexually dimorphic protective attribute that we termed immunologic resilience (IR).

OBJECTIVE

We sought to examine whether deficits in IR that antedate or are induced by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection independently predict COVID-19 mortality.

METHODS

IR levels were quantified with 2 novel metrics: immune health grades (IHG-I [best] to IHG-IV) to gauge CD8⁺ and CD4⁺ T-cell count equilibrium, and blood gene expression signatures. IR metrics were examined in a prospective COVID-19 cohort (n = 522); primary outcome was 30-day mortality. Associations of IR metrics with outcomes in non-COVID-19 cohorts (n = 13,461) provided the framework for linking pre-COVID-19 IR status to IR during COVID-19, as well as to COVID-19 outcomes.

RESULTS

IHG-I, tracking high-grade equilibrium between CD8⁺ and CD4⁺ T-cell counts, was the most common grade (73%) among healthy adults, particularly in females. SARS-CoV-2 infection was associated with underrepresentation of IHG-I (21%) versus overrepresentation (77%) of IHG-II or IHG-IV, especially in males versus females (P < .01). Presentation with IHG-I was associated with 88% lower mortality, after controlling for age and sex; reduced risk of hospitalization and respiratory failure; lower plasma IL-6 levels; rapid clearance of nasopharyngeal SARS-CoV-2 burden; and gene expression signatures correlating with survival that signify immunocompetence and controlled inflammation. In non-COVID-19 cohorts, IR-preserving metrics were associated with resistance to progressive influenza or HIV infection, as well as lower 9-year mortality in the Framingham Heart Study, especially in females.

CONCLUSIONS

Preservation of immunocompetence with controlled inflammation during antigenic challenges is a hallmark of IR and associates with longevity and AIDS resistance. Independent of age, a male-biased proclivity to degrade IR before and/or during SARS-CoV-2 infection predisposes to severe COVID-19.

Ten millennia of hepatitis B virus evolution

[Figure: see text].

Application of Genetic Studies to Flow Cytometry Data and Its Impact on Therapeutic Intervention for Autoimmune Disease

In recent years, systematic genome-wide association studies of quantitative immune cell traits, represented by circulating levels of cell subtypes established by flow cytometry, have revealed numerous association signals, a large fraction of which overlap perfectly with genetic signals associated with autoimmune diseases. By identifying further overlaps with association signals influencing gene expression and cell surface protein levels, it has also been possible, in several cases, to identify causal genes and infer candidate proteins affecting immune cell traits linked to autoimmune disease risk. Overall, these results provide a more detailed picture of how genetic variation affects the human immune system and autoimmune disease risk. They also highlight druggable proteins in the pathogenesis of autoimmune diseases; predict the efficacy and side effects of existing therapies; provide new indications for use for some of them; and optimize the research and development of new, more effective and safer treatments for autoimmune diseases. Here we review the genetic-driven approach that couples systematic multi-parametric flow cytometry with high-resolution genetics and transcriptomics to identify endophenotypes of autoimmune diseases for the development of new therapies.

Mitochondrial genetic variation is enriched in G-quadruplex regions that stall DNA synthesis in vitro

As the powerhouses of the eukaryotic cell, mitochondria must maintain their genomes which encode proteins essential for energy production. Mitochondria are characterized by guanine-rich DNA sequences that spontaneously form unusual three-dimensional structures known as G-quadruplexes (G4). G4 structures can be problematic for the essential processes of DNA replication and transcription because they deter normal progression of the enzymatic-driven processes. In this study, we addressed the hypothesis that mitochondrial G4 is a source of mutagenesis leading to base-pair substitutions. Our computational analysis of 2757 individual genomes from two Italian population cohorts (SardiNIA and InCHIANTI) revealed a statistically significant enrichment of mitochondrial mutations within sequences corresponding to stable G4 DNA structures. Guided by the computational analysis results, we designed biochemical reconstitution experiments and demonstrated that DNA synthesis by two known mitochondrial DNA polymerases (Pol γ, PrimPol) in vitro was strongly blocked by representative stable G4 mitochondrial DNA structures, which could be overcome in a specific manner by the ATP-dependent G4-resolving helicase Pif1. However, error-prone DNA synthesis by PrimPol using the G4 template sequence persisted even in the presence of Pif1. Altogether, our results suggest that genetic variation is enriched in G-quadruplex regions that impede mitochondrial DNA replication.

Genetic insights into biological mechanisms governing human ovarian ageing

Reproductive longevity is essential for fertility and influences healthy ageing in women1,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 premutations3. 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.

Genetic variant within CDK6 regulates immune response to palbociclib treatment

Everyone carries a set of genetic variants that contribute to regulation of the levels of blood cells, with unknown clinical impact. One of them, rs445 within the cell-cycle checkpoint gene CDK6, reduces the levels of myeloid cell types including granulocytes. We treated CD3+ T cells and whole blood with palbociclib in 41 individuals, who were stratified by genotype for analyses. In T cells we assessed cell cycle and apoptosis, whereas in whole blood, apoptosis in activated (CD11b+), unactivated (CD11b-) granulocytes, cytotoxic (CD8 + CD4-), and helper (CD8-CD4+) T cells. We find that rs445 modulates the immune response of CD8+ T cells. It also increases the level of apoptotic CD11b + activated granulocytes after palbociclib treatment, which, in synergy with neutropenia, may affect drug related adverse events. These results suggest that the effect of palbociclib treatment may depend on underlying genetically encoded individual immune response as well as the direct response to the drug.

A robust machine learning framework to identify signatures for frailty: a nested case-control study in four aging European cohorts

Phenotype-specific omic expression patterns in people with frailty could provide invaluable insight into the underlying multi-systemic pathological processes and targets for intervention. Classical approaches to frailty have not considered the potential for different frailty phenotypes. We characterized associations between frailty (with/without disability) and sets of omic factors (genomic, proteomic, and metabolomic) plus markers measured in routine geriatric care. This study was a prevalent case control using stored biospecimens (urine, whole blood, cells, plasma, and serum) from 1522 individuals (identified as robust (R), pre-frail (P), or frail (F)] from the Toledo Study of Healthy Aging (R=178/P=184/F=109), 3 City Bordeaux (111/269/100), Aging Multidisciplinary Investigation (157/79/54) and InCHIANTI (106/98/77) cohorts. The analysis included over 35,000 omic and routine laboratory variables from robust and frail or pre-frail (with/without disability) individuals using a machine learning framework. We identified three protective biomarkers, vitamin D3 (OR: 0.81 [95% CI: 0.68-0.98]), lutein zeaxanthin (OR: 0.82 [95% CI: 0.70-0.97]), and miRNA125b-5p (OR: 0.73, [95% CI: 0.56-0.97]) and one risk biomarker, cardiac troponin T (OR: 1.25 [95% CI: 1.23-1.27]). Excluding individuals with a disability, one protective biomarker was identified, miR125b-5p (OR: 0.85, [95% CI: 0.81-0.88]). Three risks of frailty biomarkers were detected: pro-BNP (OR: 1.47 [95% CI: 1.27-1.7]), cardiac troponin T (OR: 1.29 [95% CI: 1.21-1.38]), and sRAGE (OR: 1.26 [95% CI: 1.01-1.57]). Three key frailty biomarkers demonstrated a statistical association with frailty (oxidative stress, vitamin D, and cardiovascular system) with relationship patterns differing depending on the presence or absence of a disability.

Analysis of overlapping genetic association in type 1 and type 2 diabetes

AIMS/HYPOTHESIS

Given the potential shared aetiology between type 1 and type 2 diabetes, we aimed to identify any genetic regions associated with both diseases. For associations where there is a shared signal and the allele that increases risk to one disease also increases risk to the other, inference about shared aetiology could be made, with the potential to develop therapeutic strategies to treat or prevent both diseases simultaneously. Alternatively, if a genetic signal co-localises with divergent effect directions, it could provide valuable biological insight into how the association affects the two diseases differently.

METHODS

Using publicly available type 2 diabetes summary statistics from a genome-wide association study (GWAS) meta-analysis of European ancestry individuals (74,124 cases and 824,006 controls) and type 1 diabetes GWAS summary statistics from a meta-analysis of studies on individuals from the UK and Sardinia (7467 cases and 10,218 controls), we identified all regions of 0.5 Mb that contained variants associated with both diseases (false discovery rate <0.01). In each region, we performed forward stepwise logistic regression to identify independent association signals, then examined co-localisation of each type 1 diabetes signal with each type 2 diabetes signal using coloc. Any association with a co-localisation posterior probability of ≥0.9 was considered a genuine shared association with both diseases.

RESULTS

Of the 81 association signals from 42 genetic regions that showed association with both type 1 and type 2 diabetes, four association signals co-localised between both diseases (posterior probability ≥0.9): (1) chromosome 16q23.1, near CTRB1/BCAR1, which has been previously identified; (2) chromosome 11p15.5, near the INS gene; (3) chromosome 4p16.3, near TMEM129 and (4) chromosome 1p31.3, near PGM1. In each of these regions, the effect of genetic variants on type 1 diabetes was in the opposite direction to the effect on type 2 diabetes. Use of additional datasets also supported the previously identified co-localisation on chromosome 9p24.2, near the GLIS3 gene, in this case with a concordant direction of effect.

CONCLUSIONS/INTERPRETATION

Four of five association signals that co-localise between type 1 diabetes and type 2 diabetes are in opposite directions, suggesting a complex genetic relationship between the two diseases.

The trans-ancestral genomic architecture of glycemic traits

Glycemic traits are used to diagnose and monitor type 2 diabetes and cardiometabolic health. To date, most genetic studies of glycemic traits have focused on individuals of European ancestry. Here we aggregated genome-wide association studies comprising up to 281,416 individuals without diabetes (30% non-European ancestry) for whom fasting glucose, 2-h glucose after an oral glucose challenge, glycated hemoglobin and fasting insulin data were available. Trans-ancestry and single-ancestry meta-analyses identified 242 loci (99 novel; P < 5 × 10-8), 80% of which had no significant evidence of between-ancestry heterogeneity. Analyses restricted to individuals of European ancestry with equivalent sample size would have led to 24 fewer new loci. Compared with single-ancestry analyses, equivalent-sized trans-ancestry fine-mapping reduced the number of estimated variants in 99% credible sets by a median of 37.5%. Genomic-feature, gene-expression and gene-set analyses revealed distinct biological signatures for each trait, highlighting different underlying biological pathways. Our results increase our understanding of diabetes pathophysiology by using trans-ancestry studies for improved power and resolution.

Low-Risk Human Leukocyte Antigen Genes and Mild Villous Atrophy Typify Celiac Disease With Immunoglobulin A Deficiency

OBJECTIVES

We aimed to establish if in celiac disease (CD) with immunoglobulin A deficiency (IgAD) duodenal histopathology is influenced by human leukocyte antigen (HLA)-DQB1∗02 alleles dosage. Clinical differences between patients with CD and patients with CD and IgAD (CD-IgAD) were also evaluated.

METHODS

Five hundred and sixteen CD and 16 patients with CD-IgAD, enrolled over the time of 8 years, took part in this study. The severity of duodenal histopathology and frequency of CD at-risk HLA class II genes were compared in patients with CD versus patients with CD-IgAD. HLA class II genotypes were subdivided into two categories of genetic risk: high: HLA-DR3/DR7, -DR3/DR3, -DR4/DR4 -DR3/DR4 and low: HLA-DR5/DR7, -DR3/X, -DR4/X and X/X, where X means neither -DR3 nor -DR4. Then, they were compared with two types of duodenal histopathology: 0, 1, 2 and 3a of mild villous atrophy (MVA) and 3b and 3c of severe villous atrophy (SVA) according to the Marsh-Oberhuber classification. Clinical data concerning gender, number of esophagogastroduodenoscopies (EGDs) and association with other autoimmune diseases were obtained from medical records.

RESULTS

In comparison with CD, CD-IgAD showed an increased frequency of MVA (P < 0.0001). Furthermore, CD-IgAD with MVA showed an increase of HLA low-risk genotypes (P = 0.036) and half HLA-DQ2 heterodimers (P = 0.0443). Interestingly, CD-IgAD demanded an increased number of EGDs to reach the diagnosis of CD (P = 0.0104) and autoimmune liver diseases were more frequent compared to CD (P = 0.0049).

CONCLUSIONS

CD-IgAD is associated with MVA, low-risk HLA class II genes, an increased number of EGDs and autoimmune liver diseases.

Genome-wide search for genes affecting the age at diagnosis of type 1 diabetes

BACKGROUND

Type 1 diabetes (T1D) is an autoimmune disease affecting individuals in the early years of life. Although previous studies have identified genetic loci influencing T1D diagnosis age, these studies did not investigate the genome with high resolution.

OBJECTIVE AND METHODS

We performed a genome-wide meta-analysis for age at diagnosis with cohorts from Finland (Finnish Diabetic Nephropathy Study), the United Kingdom (UK Genetic Resource Investigating Diabetes) and Sardinia. Through SNP associations, transcriptome-wide association analysis linked T1D diagnosis age and gene expression.

RESULTS

We identified two chromosomal regions associated with T1D diagnosis age: multiple independent variants in the HLA region on chromosome 6 and a locus on chromosome 17q12. We performed gene-level association tests with transcriptome prediction models from two whole blood datasets, lymphocyte cell line, spleen, pancreas and small intestine tissues. Of the non-HLA genes, lower PNMT expression in whole blood, and higher IKZF3 and ZPBP2, and lower ORMDL3 and GSDMB transcription levels in multiple tissues were associated with lower T1D diagnosis age (FDR = 0.05). These genes lie on chr17q12 which is associated with T1D, other autoimmune diseases, and childhood asthma. Additionally, higher expression of PHF20L1, a gene not previously implicated in T1D, was associated with lower diagnosis age in lymphocytes, pancreas, and spleen. Altogether, the non-HLA associations were enriched in open chromatin in various blood cells, blood vessel tissues and foetal thymus tissue.

CONCLUSION

Multiple genes on chr17q12 and PHF20L1 on chr8 were associated with T1D diagnosis age and only further studies may elucidate the role of these genes for immunity and T1D onset.

PRF1 mutation alters immune system activation, inflammation, and risk of autoimmunity

BACKGROUND

Defective alleles within the PRF1 gene, encoding the pore-forming protein perforin, in combination with environmental factors, cause familial type 2 hemophagocytic lymphohistiocytosis (FHL2), a rare, severe autosomal recessive childhood disorder characterized by massive release of cytokines-cytokine storm.

OBJECTIVE

The aim of this study was to determine the function of hypomorph PRF1:p.A91V g.72360387 G > A on multiple sclerosis (MS) and type 1 diabetes (T1D).

METHODS

We cross-compare the association data for PRF1:p.A91V mutation derived from GWAS on adult MS and pediatric T1D in Sardinians. The novel association with T1D was replicated in metanalysis in 12,584 cases and 17,692 controls from Sardinia, the United Kingdom, and Scotland. To dissect this mutation function, we searched through the coincident association immunophenotypes in additional set of general population Sardinians.

RESULTS

We report that PRF1:p.A91V, is associated with increase of lymphocyte levels, especially within the cytotoxic memory T-cells, at general population level with reduced interleukin 7 receptor expression on these cells. The minor allele increased risk of MS, in 2903 cases and 2880 controls from Sardinia p = 2.06 × 10^-4, odds ratio OR = 1.29, replicating a previous finding, whereas it protects from T1D p = 1.04 × 10^-5, OR = 0.82.

CONCLUSION

Our results indicate opposing contributions of the cytotoxic T-cell compartment to MS and T1D pathogenesis.

Facing the urgency of therapies for progressive MS - a Progressive MS Alliance proposal

Therapies for infiltrative inflammation in multiple sclerosis (MS) have advanced greatly, but neurodegeneration and compartmentalized inflammation remain virtually untargeted as in other diseases of the nervous system. Consequently, many therapies are available for the relapsing-remitting form of MS, but the progressive forms remain essentially untreated. The objective of the International Progressive MS Alliance is to expedite the development of effective therapies for progressive MS through new initiatives that foster innovative thinking and concrete advancements. Based on these principles, the Alliance is developing a new funding programme that will focus on experimental medicine trials. Here, we discuss the reasons behind the focus on experimental medicine trials, the strengths and weaknesses of these approaches and of the programme, and why we hope to advance therapies while improving the understanding of progression in MS. We are soliciting public and academic feedback, which will help shape the programme and future strategies of the Alliance.

Loss-of-function genomic variants highlight potential therapeutic targets for cardiovascular disease

Pharmaceutical drugs targeting dyslipidemia and cardiovascular disease (CVD) may increase the risk of fatty liver disease and other metabolic disorders. To identify potential novel CVD drug targets without these adverse effects, we perform genome-wide analyses of participants in the HUNT Study in Norway (n = 69,479) to search for protein-altering variants with beneficial impact on quantitative blood traits related to cardiovascular disease, but without detrimental impact on liver function. We identify 76 (11 previously unreported) presumed causal protein-altering variants associated with one or more CVD- or liver-related blood traits. Nine of the variants are predicted to result in loss-of-function of the protein. This includes ZNF529:p.K405X, which is associated with decreased low-density-lipoprotein (LDL) cholesterol (P = 1.3 × 10-8) without being associated with liver enzymes or non-fasting blood glucose. Silencing of ZNF529 in human hepatoma cells results in upregulation of LDL receptor and increased LDL uptake in the cells. This suggests that inhibition of ZNF529 or its gene product should be prioritized as a novel candidate drug target for treating dyslipidemia and associated CVD.

Induction of therapeutic levels of HbF in genome-edited primary β0 39-thalassaemia haematopoietic stem and progenitor cells

Hereditary persistence of fetal haemoglobin (HPFH) is the major modifier of the clinical severity of β-thalassaemia. The homozygous mutation c.-196 C>T in the Aγ-globin (HBG1) promoter, which causes Sardinian δβ⁰ -thalassaemia, is able to completely rescue the β-major thalassaemia phenotype caused by the β⁰ 39-thalassaemia mutation, ensuring high levels of fetal haemoglobin synthesis during adulthood. Here, we describe a CRISPR/Cas9 genome-editing approach, combined with the non-homologous end joining (NHEJ) pathway repair, aimed at reproducing the effects of this naturally occurring HPFH mutation in both HBG promoters. After selecting the most efficient guide RNA in K562 cells, we edited the HBG promoters in human umbilical cord blood-derived erythroid progenitor 2 cells (HUDEP-2) and in haematopoietic stem and progenitor cells (HSPCs) from β⁰ -thalassaemia patients to assess the therapeutic potential of HbF induction. Our results indicate that small deletions targeting the -196-promoter region restore high levels of fetal haemoglobin (HbF) synthesis in all cell types tested. In pools of HSPCs derived from homozygous β⁰ 39-thalassaemia patients, a 20% editing determined a parallel 20% increase of HbF compared to unedited pools. These results suggest that editing the region of HBG promoters around the -196 position has the potential to induce therapeutic levels of HbF in patients with most types of β-thalassaemia irrespective of the β-globin gene (HBB) mutations.

Kidney size in relation to ageing, gender, renal function, birthweight and chronic kidney disease risk factors in a general population

BACKGROUND

The relationship of kidney size to ageing, kidney function and kidney disease risk factors is not fully understood.

METHODS

Ultrasound length and parenchymal kidney volume were determined from a population-based sample of 3972 Sardinians (age range 18-100 years). We then identified the subset of 2256 'healthy' subjects to define age- and sex-specific reference ranges (2.5-97.5 percentile) of kidney volume. Logistic regression (accounting for family clustering) was used to identify the clinical characteristics associated with abnormally large kidneys or abnormally small kidneys.

RESULTS

In the healthy subset, kidney volume and length increased up to the fourth to fifth decade of life followed by a progressive decrease in men, whereas there was a gradual kidney volume decrease throughout the lifespan of women. In the whole sample, independent predictors of lower kidney volume (<2.5 percentile for age and sex) were male sex, low body mass index, short height, low waist:hip ratio and high serum creatinine (SCr); the independent predictors of larger kidney volume (>97.5 percentile for age and sex) were younger age, female sex, diabetes, obesity, high height, high waist:hip ratio and lower SCr. Estimated heritability for kidney volume was 15%, and for length 27%; kidney volume correlated strongly with birthweight.

CONCLUSIONS

Overall, in a general healthy population, kidney measures declined with age differently in men and women. The determinants of kidney parenchymal volume include genetic factors and modifiable clinical factors.

Evolutionarily Selected Overexpression of the Cytokine BAFF Enhances Mucosal Immune Response Against P. falciparum

We have previously shown that a variant of the TNFSF13B gene that we called BAFF-var increases the production of the cytokine BAFF, upregulating humoral immunity and increasing the risk for certain autoimmune diseases. In addition, genetic population signatures revealed that BAFF-var was evolutionarily advantageous, most likely by increasing resistance to malaria infection, which is a prime candidate for selective pressure. To evaluate whether the increased soluble BAFF (sBAFF) production confers protection, we experimentally assessed the role of BAFF-var in response to malaria antigens. Lysates of erythrocytes infected with Plasmodium falciparum (iRBCs) or left uninfected (uRBCs, control) were used to treat peripheral blood mononuclear cells (PBMCs) with distinct BAFF genotypes. The PBMCs purified from BAFF-var donors and treated with iRBCs showed different levels of specific cells, immunoglobulins, and cytokines as compared with BAFF-WT. In particular, a relevant differential effect on mucosal immunity B subpopulations have been observed. These findings point to specific immune cells and molecules through which the evolutionary selected BAFF-var may have improved fitness during P. falciparum infection.

Meta-analysis of up to 622,409 individuals identifies 40 novel smoking behaviour associated genetic loci

Smoking is a major heritable and modifiable risk factor for many diseases, including cancer, common respiratory disorders and cardiovascular diseases. Fourteen genetic loci have previously been associated with smoking behaviour-related traits. We tested up to 235,116 single nucleotide variants (SNVs) on the exome-array for association with smoking initiation, cigarettes per day, pack-years, and smoking cessation in a fixed effects meta-analysis of up to 61 studies (up to 346,813 participants). In a subset of 112,811 participants, a further one million SNVs were also genotyped and tested for association with the four smoking behaviour traits. SNV-trait associations with P < 5 × 10-8 in either analysis were taken forward for replication in up to 275,596 independent participants from UK Biobank. Lastly, a meta-analysis of the discovery and replication studies was performed. Sixteen SNVs were associated with at least one of the smoking behaviour traits (P < 5 × 10-8) in the discovery samples. Ten novel SNVs, including rs12616219 near TMEM182, were followed-up and five of them (rs462779 in REV3L, rs12780116 in CNNM2, rs1190736 in GPR101, rs11539157 in PJA1, and rs12616219 near TMEM182) replicated at a Bonferroni significance threshold (P < 4.5 × 10-3) with consistent direction of effect. A further 35 SNVs were associated with smoking behaviour traits in the discovery plus replication meta-analysis (up to 622,409 participants) including a rare SNV, rs150493199, in CCDC141 and two low-frequency SNVs in CEP350 and HDGFRP2. Functional follow-up implied that decreased expression of REV3L may lower the probability of smoking initiation. The novel loci will facilitate understanding the genetic aetiology of smoking behaviour and may lead to the identification of potential drug targets for smoking prevention and/or cessation.

Author Correction: Complex genetic signatures in immune cells underlie autoimmunity and inform therapy

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Complex genetic signatures in immune cells underlie autoimmunity and inform therapy

We report on the influence of ~22 million variants on 731 immune cell traits in a cohort of 3,757 Sardinians. We detected 122 significant (P < 1.28 × 10^-11) independent association signals for 459 cell traits at 70 loci (53 of them novel) identifying several molecules and mechanisms involved in cell regulation. Furthermore, 53 signals at 36 loci overlapped with previously reported disease-associated signals, predominantly for autoimmune disorders, highlighting intermediate phenotypes in pathogenesis. Collectively, our findings illustrate complex genetic regulation of immune cells with highly selective effects on autoimmune disease risk at the cell-subtype level. These results identify drug-targetable pathways informing the design of more specific treatments for autoimmune diseases.