Clinical characteristics with inflammation profiling of long COVID and association with 1-year recovery following hospitalisation in the UK: a prospective observational study

BACKGROUND

No effective pharmacological or non-pharmacological interventions exist for patients with long COVID. We aimed to describe recovery 1 year after hospital discharge for COVID-19, identify factors associated with patient-perceived recovery, and identify potential therapeutic targets by describing the underlying inflammatory profiles of the previously described recovery clusters at 5 months after hospital discharge.

METHODS

The Post-hospitalisation COVID-19 study (PHOSP-COVID) is a prospective, longitudinal cohort study recruiting adults (aged ≥18 years) discharged from hospital with COVID-19 across the UK. Recovery was assessed using patient-reported outcome measures, physical performance, and organ function at 5 months and 1 year after hospital discharge, and stratified by both patient-perceived recovery and recovery cluster. Hierarchical logistic regression modelling was performed for patient-perceived recovery at 1 year. Cluster analysis was done using the clustering large applications k-medoids approach using clinical outcomes at 5 months. Inflammatory protein profiling was analysed from plasma at the 5-month visit. This study is registered on the ISRCTN Registry, ISRCTN10980107, and recruitment is ongoing.

FINDINGS

2320 participants discharged from hospital between March 7, 2020, and April 18, 2021, were assessed at 5 months after discharge and 807 (32·7%) participants completed both the 5-month and 1-year visits. 279 (35·6%) of these 807 patients were women and 505 (64·4%) were men, with a mean age of 58·7 (SD 12·5) years, and 224 (27·8%) had received invasive mechanical ventilation (WHO class 7-9). The proportion of patients reporting full recovery was unchanged between 5 months (501 [25·5%] of 1965) and 1 year (232 [28·9%] of 804). Factors associated with being less likely to report full recovery at 1 year were female sex (odds ratio 0·68 [95% CI 0·46-0·99]), obesity (0·50 [0·34-0·74]) and invasive mechanical ventilation (0·42 [0·23-0·76]). Cluster analysis (n=1636) corroborated the previously reported four clusters: very severe, severe, moderate with cognitive impairment, and mild, relating to the severity of physical health, mental health, and cognitive impairment at 5 months. We found increased inflammatory mediators of tissue damage and repair in both the very severe and the moderate with cognitive impairment clusters compared with the mild cluster, including IL-6 concentration, which was increased in both comparisons (n=626 participants). We found a substantial deficit in median EQ-5D-5L utility index from before COVID-19 (retrospective assessment; 0·88 [IQR 0·74-1·00]), at 5 months (0·74 [0·64-0·88]) to 1 year (0·75 [0·62-0·88]), with minimal improvements across all outcome measures at 1 year after discharge in the whole cohort and within each of the four clusters.

INTERPRETATION

The sequelae of a hospital admission with COVID-19 were substantial 1 year after discharge across a range of health domains, with the minority in our cohort feeling fully recovered. Patient-perceived health-related quality of life was reduced at 1 year compared with before hospital admission. Systematic inflammation and obesity are potential treatable traits that warrant further investigation in clinical trials.

FUNDING

UK Research and Innovation and National Institute for Health Research.

Spontaneous coronary artery dissection: novel pathophysiological insights from histological and ultrastructural tissue analysis

Introduction

Spontaneous coronary artery dissection (SCAD) is a cause of acute coronary syndromes and rarely sudden cardiac death (SCD). SCAD is characterised by medial false lumen haematoma formation and periadventitial inflammatory cell infiltrate. Although SCAD has been linked to connective tissue disorders, its pathophysiology remains poorly understood and the role of inflammation unknown.

Purpose

We sought to establish the definitive histopathological features of SCAD and explore pathophysiological mechanisms through assessment of dermal connective tissue ultrastructure.

Methods

N=36 SCD cases diagnosed as SCAD on autopsy were identified in pathology archives at four international centres. Their demographic and clinical characteristics were compared with n=359 survivors recruited in a SCAD survivors cohort. Haematoxylin & eosin sections were examined under light microscope. Immunohistochemistry (IHC) was employed for quantification of inflammatory cell infiltrate (CD68, CD3) and vasa vasorum density (CD31) of SCAD cases (n=20) compared to age- and sex-matched controls (n=10). Dermal extracellular matrix components (EMC) of n=32 SCAD survivors and n=16 healthy volunteers (HV) were compared using electron microscopy (EM).

Results

The autopsy series cases were more likely to be male (p=0.0256) and had higher incidence of left main stem (p=0.0475) and proximal left anterior descending (p<0.001) disease compared to SCAD survivors. N=24 (66%) of SCAD autopsy case showed no evidence of myocardial necrosis. N=17 (47%) showed mild-moderate atherosclerotic changes but no features of fibromuscular dysplasia. There were no differences in vasa vasorum density between SCAD and control cases (A). The degree of inflammatory cell infiltrate varied greatly but significantly higher than controls (B), comprising CD68+ macrophages, eosinophils and CD3+ positive T-cells. There was a statistically significant association (p=0.006) between the degree of inflammatory cell infiltrate and the length of time from onset of symptoms to death (Panel C), as well as significantly (p<0.001) denser inflammatory cell infiltrate adjacent to the dissection plane (D, exemplary sections E&F). EM revealed no differences between SCAD and HV in dermal fibroblast size & activity or elastin size & damage indicators, but possible changes in subgroups with more extreme clinical phenotype or pregnancy-related SCAD (G).

Conclusions

To our knowledge this is the largest SCAD pathology case series so far. We show for the first time that periadvential inflammation in SCAD appears to be time-dependent and localising to the dissected coronary segment, suggesting healing response to injury rather than causal contribution. We found no evidence to suggest increased vasa vasorum density is pathophysiologically important. Connective tissue changes were only linked to a small proportion of cases. These novel findings may have important implications for the management of SCAD patients.

Funding Acknowledgement

Type of funding source: Foundation. Main funding source(s): British Heart Foundation, Leicester NIHR Biomedical Research Centre

A Method to Exploit the Structure of Genetic Ancestry Space to Enhance Case-Control Studies

One goal of human genetics is to understand the genetic basis of disease, a challenge for diseases of complex inheritance because risk alleles are few relative to the vast set of benign variants. Risk variants are often sought by association studies in which allele frequencies in case subjects are contrasted with those from population-based samples used as control subjects. In an ideal world we would know population-level allele frequencies, releasing researchers to focus on case subjects. We argue this ideal is possible, at least theoretically, and we outline a path to achieving it in reality. If such a resource were to exist, it would yield ample savings and would facilitate the effective use of data repositories by removing administrative and technical barriers. We call this concept the Universal Control Repository Network (UNICORN), a means to perform association analyses without necessitating direct access to individual-level control data. Our approach to UNICORN uses existing genetic resources and various statistical tools to analyze these data, including hierarchical clustering with spectral analysis of ancestry; and empirical Bayesian analysis along with Gaussian spatial processes to estimate ancestry-specific allele frequencies. We demonstrate our approach using tens of thousands of control subjects from studies of Crohn disease, showing how it controls false positives, provides power similar to that achieved when all control data are directly accessible, and enhances power when control data are limiting or even imperfectly matched ancestrally. These results highlight how UNICORN can enable reliable, powerful, and convenient genetic association analyses without access to the individual-level data.

Effects of long-term averaging of quantitative blood pressure traits on the detection of genetic associations

Blood pressure (BP) is a heritable, quantitative trait with intraindividual variability and susceptibility to measurement error. Genetic studies of BP generally use single-visit measurements and thus cannot remove variability occurring over months or years. We leveraged the idea that averaging BP measured across time would improve phenotypic accuracy and thereby increase statistical power to detect genetic associations. We studied systolic BP (SBP), diastolic BP (DBP), mean arterial pressure (MAP), and pulse pressure (PP) averaged over multiple years in 46,629 individuals of European ancestry. We identified 39 trait-variant associations across 19 independent loci (p < 5 × 10(-8)); five associations (in four loci) uniquely identified by our LTA analyses included those of SBP and MAP at 2p23 (rs1275988, near KCNK3), DBP at 2q11.2 (rs7599598, in FER1L5), and PP at 6p21 (rs10948071, near CRIP3) and 7p13 (rs2949837, near IGFBP3). Replication analyses conducted in cohorts with single-visit BP data showed positive replication of associations and a nominal association (p < 0.05). We estimated a 20% gain in statistical power with long-term average (LTA) as compared to single-visit BP association studies. Using LTA analysis, we identified genetic loci influencing BP. LTA might be one way of increasing the power of genetic associations for continuous traits in extant samples for other phenotypes that are measured serially over time.

Meta-analysis of gene-level associations for rare variants based on single-variant statistics

Meta-analysis of genome-wide association studies (GWASs) has led to the discoveries of many common variants associated with complex human diseases. There is a growing recognition that identifying "causal" rare variants also requires large-scale meta-analysis. The fact that association tests with rare variants are performed at the gene level rather than at the variant level poses unprecedented challenges in the meta-analysis. First, different studies may adopt different gene-level tests, so the results are not compatible. Second, gene-level tests require multivariate statistics (i.e., components of the test statistic and their covariance matrix), which are difficult to obtain. To overcome these challenges, we propose to perform gene-level tests for rare variants by combining the results of single-variant analysis (i.e., p values of association tests and effect estimates) from participating studies. This simple strategy is possible because of an insight that multivariate statistics can be recovered from single-variant statistics, together with the correlation matrix of the single-variant test statistics, which can be estimated from one of the participating studies or from a publicly available database. We show both theoretically and numerically that the proposed meta-analysis approach provides accurate control of the type I error and is as powerful as joint analysis of individual participant data. This approach accommodates any disease phenotype and any study design and produces all commonly used gene-level tests. An application to the GWAS summary results of the Genetic Investigation of ANthropometric Traits (GIANT) consortium reveals rare and low-frequency variants associated with human height. The relevant software is freely available.

Large-scale gene-centric meta-analysis across 32 studies identifies multiple lipid loci

Genome-wide association studies (GWASs) have identified many SNPs underlying variations in plasma-lipid levels. We explore whether additional loci associated with plasma-lipid phenotypes, such as high-density lipoprotein cholesterol (HDL-C), low-density lipoprotein cholesterol (LDL-C), total cholesterol (TC), and triglycerides (TGs), can be identified by a dense gene-centric approach. Our meta-analysis of 32 studies in 66,240 individuals of European ancestry was based on the custom ∼50,000 SNP genotyping array (the ITMAT-Broad-CARe array) covering ∼2,000 candidate genes. SNP-lipid associations were replicated either in a cohort comprising an additional 24,736 samples or within the Global Lipid Genetic Consortium. We identified four, six, ten, and four unreported SNPs in established lipid genes for HDL-C, LDL-C, TC, and TGs, respectively. We also identified several lipid-related SNPs in previously unreported genes: DGAT2, HCAR2, GPIHBP1, PPARG, and FTO for HDL-C; SOCS3, APOH, SPTY2D1, BRCA2, and VLDLR for LDL-C; SOCS3, UGT1A1, BRCA2, UBE3B, FCGR2A, CHUK, and INSIG2 for TC; and SERPINF2, C4B, GCK, GATA4, INSR, and LPAL2 for TGs. The proportion of explained phenotypic variance in the subset of studies providing individual-level data was 9.9% for HDL-C, 9.5% for LDL-C, 10.3% for TC, and 8.0% for TGs. This large meta-analysis of lipid phenotypes with the use of a dense gene-centric approach identified multiple SNPs not previously described in established lipid genes and several previously unknown loci. The explained phenotypic variance from this approach was comparable to that from a meta-analysis of GWAS data, suggesting that a focused genotyping approach can further increase the understanding of heritability of plasma lipids.

Genetic associations for activated partial thromboplastin time and prothrombin time, their gene expression profiles, and risk of coronary artery disease

Activated partial thromboplastin time (aPTT) and prothrombin time (PT) are clinical tests commonly used to screen for coagulation-factor deficiencies. One genome-wide association study (GWAS) has been reported previously for aPTT, but no GWAS has been reported for PT. We conducted a GWAS and meta-analysis to identify genetic loci for aPTT and PT. The GWAS for aPTT was conducted in 9,240 individuals of European ancestry from the Atherosclerosis Risk in Communities (ARIC) study, and the GWAS for PT was conducted in 2,583 participants from the Genetic Study of Three Population Microisolates in South Tyrol (MICROS) and the Lothian Birth Cohorts (LBC) of 1921 and 1936. Replication was assessed in 1,041 to 3,467 individuals. For aPTT, previously reported associations with KNG1, HRG, F11, F12, and ABO were confirmed. A second independent association in ABO was identified and replicated (rs8176704, p = 4.26 × 10(-24)). Pooling the ARIC and replication data yielded two additional loci in F5 (rs6028, p = 3.22 × 10(-9)) and AGBL1 (rs2469184, p = 3.61 × 10(-8)). For PT, significant associations were identified and confirmed in F7 (rs561241, p = 3.71 × 10(-56)) and PROCR/EDEM2 (rs2295888, p = 5.25 × 10(-13)). Assessment of existing gene expression and coronary artery disease (CAD) databases identified associations of five of the GWAS loci with altered gene expression and two with CAD. In summary, eight genetic loci that account for ∼29% of the variance in aPTT and two loci that account for ∼14% of the variance in PT were detected and supported by functional data.

Large-scale gene-centric meta-analysis across 39 studies identifies type 2 diabetes loci

To identify genetic factors contributing to type 2 diabetes (T2D), we performed large-scale meta-analyses by using a custom ∼50,000 SNP genotyping array (the ITMAT-Broad-CARe array) with ∼2000 candidate genes in 39 multiethnic population-based studies, case-control studies, and clinical trials totaling 17,418 cases and 70,298 controls. First, meta-analysis of 25 studies comprising 14,073 cases and 57,489 controls of European descent confirmed eight established T2D loci at genome-wide significance. In silico follow-up analysis of putative association signals found in independent genome-wide association studies (including 8,130 cases and 38,987 controls) performed by the DIAGRAM consortium identified a T2D locus at genome-wide significance (GATAD2A/CILP2/PBX4; p = 5.7 × 10(-9)) and two loci exceeding study-wide significance (SREBF1, and TH/INS; p < 2.4 × 10(-6)). Second, meta-analyses of 1,986 cases and 7,695 controls from eight African-American studies identified study-wide-significant (p = 2.4 × 10(-7)) variants in HMGA2 and replicated variants in TCF7L2 (p = 5.1 × 10(-15)). Third, conditional analysis revealed multiple known and novel independent signals within five T2D-associated genes in samples of European ancestry and within HMGA2 in African-American samples. Fourth, a multiethnic meta-analysis of all 39 studies identified T2D-associated variants in BCL2 (p = 2.1 × 10(-8)). Finally, a composite genetic score of SNPs from new and established T2D signals was significantly associated with increased risk of diabetes in African-American, Hispanic, and Asian populations. In summary, large-scale meta-analysis involving a dense gene-centric approach has uncovered additional loci and variants that contribute to T2D risk and suggests substantial overlap of T2D association signals across multiple ethnic groups.

Blood pressure loci identified with a gene-centric array

Raised blood pressure (BP) is a major risk factor for cardiovascular disease. Previous studies have identified 47 distinct genetic variants robustly associated with BP, but collectively these explain only a few percent of the heritability for BP phenotypes. To find additional BP loci, we used a bespoke gene-centric array to genotype an independent discovery sample of 25,118 individuals that combined hypertensive case-control and general population samples. We followed up four SNPs associated with BP at our p < 8.56 × 10(-7) study-specific significance threshold and six suggestively associated SNPs in a further 59,349 individuals. We identified and replicated a SNP at LSP1/TNNT3, a SNP at MTHFR-NPPB independent (r(2) = 0.33) of previous reports, and replicated SNPs at AGT and ATP2B1 reported previously. An analysis of combined discovery and follow-up data identified SNPs significantly associated with BP at p < 8.56 × 10(-7) at four further loci (NPR3, HFE, NOS3, and SOX6). The high number of discoveries made with modest genotyping effort can be attributed to using a large-scale yet targeted genotyping array and to the development of a weighting scheme that maximized power when meta-analyzing results from samples ascertained with extreme phenotypes, in combination with results from nonascertained or population samples. Chromatin immunoprecipitation and transcript expression data highlight potential gene regulatory mechanisms at the MTHFR and NOS3 loci. These results provide candidates for further study to help dissect mechanisms affecting BP and highlight the utility of studying SNPs and samples that are independent of those studied previously even when the sample size is smaller than that in previous studies.

Association analysis of beta2 adrenoceptor polymorphisms with hypertension in a Black African population

OBJECTIVE

To determine whether or not beta2 adrenoceptor polymorphism is a risk factor for the development of hypertension in a Black South African population.

BACKGROUND

Attenuated vasodilator responses to endogenous catecholamines may contribute to the aetiology of hypertension. Downregulation of beta2 adrenoreceptors (beta2AR) following stimulation with agonists is determined in part by variation at the beta2AR gene locus. The Glu27 beta2AR genotype results in attenuated downregulation compared with the wild-type Gln27 receptor, whereas Gly16 exhibits enhanced down-regulation compared to Arg16. Possible racial differences in the prevalence of the beta2AR polymorphisms may be an explanation for the blunted responses to isoprenaline and the increased prevalence of hypertension in Black African populations.

METHODS

One hundred and ninety-two unrelated hypertensives and 123 normotensives of Black South African origin were studied. Hypertensives were recruited from hospital hypertension clinics in the province of Gauteng and if on treatment, had a 2-4 week washout period before 24-h ambulatory blood pressure assessment Normotensive controls were recruited from the same community.

RESULTS

There was no significant association between either the Arg-Gly16 polymorphism or the Gln-Glu27 polymorphism and hypertension status. Furthermore, in the hypertensives, no significant association was seen between beta2AR genotype at either site and clinical blood pressure, 24-h blood pressure or left ventricular mass. A significant association was seen between Arg16 homozygotes and lower body mass index in hypertensives (P = 0.007) although this was not a primary end point. Interestingly, the Glu27 polymorphism was much rarer in this population (allelic frequency 17%) compared to a Caucasian population.

CONCLUSION

These data suggest that beta2AR polymorphism is not a risk factor for hypertension per se in this defined population. The possibility that the decreased prevalence of Glu27 in black South African populations explains blunted vasodilator responses to isoprenaline requires further study.

The effect of sodium salicylate on the rat embryo in culture: an in vitro model for the morphological assessment of teratogenicity

The optimum teratogenic dose of subcutaneously administered sodium salicylate was determined in vivo at 9.5 days of gestation. Fetuses from rats injected with this dose at 9.5 days were examined at 11.5 days and its propensities for producing resorption as well as deformities noted. Next, maternal serum levels 3-3.5 hours and 18.5 hours after salicylate injection were determined. Having established the feasibility of determining an effect due to salicylate at 11.5 days, after injection at 9.5 days, the following experiments were performed in vitro in animals cultured between 9.5 and 11.5 days by the method of New et al. (1976a): (1) Sodium salicylate was added to the culture serum, at levels equivalent to those obtained 3-3.5 hours after maternal injection of the optimum teratogenic dose, for 24 of the 48 hours culture period. (2) Rats were cultured for the first 24 hours of the 48 hours culture period in serum taken from rats injected 3.5 hours previously with the optimum teratogenic dose. (3) Rats were cultured for 24 out of 48 hours in serum from animals which had been injected with the optimum teratogenic dose 18.5 hours before bleeding. (4) Rats were cultured for 24 out of 48 hours in serum containing salicylate added to make up levels normally associated with (3). (5) A control culture was performed. The experiment indicates great similarity between the results obtained from animals cultured in serum with salicylate added and results from culture in serum of salicylate treated rats. It appears therefore that sodium salicylate - acting directly on the feto-placental unit - is the active teratogen rather than any of its metabolites.