Genetics of circulating inflammatory proteins identifies drivers of immune-mediated disease risk and therapeutic targets.

Circulating proteins have important functions in inflammation and a broad range of diseases. To identify genetic influences on inflammation-related proteins, we conducted a genome-wide protein quantitative trait locus (pQTL) study of 91 plasma proteins measured using the Olink Target platform in 14,824 participants. We identified 180 pQTLs (59 cis, 121 trans). Integration of pQTL data with eQTL and disease genome-wide association studies provided insight into pathogenesis, implicating lymphotoxin-α in multiple sclerosis. Using Mendelian randomization (MR) to assess causality in disease etiology, we identified both shared and distinct effects of specific proteins across immune-mediated diseases, including directionally discordant effects of CD40 on risk of rheumatoid arthritis versus multiple sclerosis and inflammatory bowel disease. MR implicated CXCL5 in the etiology of ulcerative colitis (UC) and we show elevated gut CXCL5 transcript expression in patients with UC. These results identify targets of existing drugs and provide a powerful resource to facilitate future drug target prioritization.

A genetic association study of circulating coagulation factor VIII and von Willebrand factor levels.

ABSTRACT: Coagulation factor VIII (FVIII) and its carrier protein von Willebrand factor (VWF) are critical to coagulation and platelet aggregation. We leveraged whole-genome sequence data from the Trans-Omics for Precision Medicine (TOPMed) program along with TOPMed-based imputation of genotypes in additional samples to identify genetic associations with circulating FVIII and VWF levels in a single-variant meta-analysis, including up to 45 289 participants. Gene-based aggregate tests were implemented in TOPMed. We identified 3 candidate causal genes and tested their functional effect on FVIII release from human liver endothelial cells (HLECs) and VWF release from human umbilical vein endothelial cells. Mendelian randomization was also performed to provide evidence for causal associations of FVIII and VWF with thrombotic outcomes. We identified associations (P < 5 × 10-9) at 7 new loci for FVIII (ST3GAL4, CLEC4M, B3GNT2, ASGR1, F12, KNG1, and TREM1/NCR2) and 1 for VWF (B3GNT2). VWF, ABO, and STAB2 were associated with FVIII and VWF in gene-based analyses. Multiphenotype analysis of FVIII and VWF identified another 3 new loci, including PDIA3. Silencing of B3GNT2 and the previously reported CD36 gene decreased release of FVIII by HLECs, whereas silencing of B3GNT2, CD36, and PDIA3 decreased release of VWF by HVECs. Mendelian randomization supports causal association of higher FVIII and VWF with increased risk of thrombotic outcomes. Seven new loci were identified for FVIII and 1 for VWF, with evidence supporting causal associations of FVIII and VWF with thrombotic outcomes. B3GNT2, CD36, and PDIA3 modulate the release of FVIII and/or VWF in vitro.

Variability of human Alpha-1-acid glycoprotein N-glycome in a Caucasian population.

AIM: Alpha-1-acid glycoprotein (AGP) is a highly glycosylated protein in human plasma and one of the most abundant acute phase proteins in humans. Glycosylation plays a crucial role in its biological functions, and alterations in AGP N-glycome have been associated with various diseases and inflammatory conditions. However, large-scale studies of AGP N-glycosylation in the general population are lacking. METHODS: Using recently developed high-throughput glycoproteomic workflow for site-specific AGP N-glycosylation analysis, 803 individuals from the Croatian island of Korcula were analyzed and their AGP N-glycome data associated with biochemical and physiological traits, as well as different environmental factors. RESULTS: After regression analysis, we found that AGP N-glycosylation is strongly associated with sex, somewhat less with age, along with multiple biochemical and physiological traits (e.g. BMI, triglycerides, uric acid, glucose, smoking status, fibrinogen). CONCLUSION: For the first time we have extensively explored the inter-individual variability of AGP N-glycome in a general human population, demonstrating its changes with sex, age, biochemical, and physiological status of individuals, providing the baseline for future population and clinical studies.

The Effect of Mediterranean Diet on Thyroid Gland Activity.

The main goal of this research was to determine whether there is a correlation between adherence to the Mediterranean diet (assessed by the Mediterranean Diet Serving Score (MDSS)) and parameters indicating thyroid gland activity, such as concentration of thyroid-stimulating hormone (TSH), thyroid hormones (free triiodothyronine (fT3), free thyroxine (fT4)), thyroglobulin (Tg), antibodies to thyroid proteins (thyroglobulin antibodies (TgAb) and thyroid peroxidase antibodies (TPOAb)), and calcitonin (CT) in plasma and serum samples. An additional objective was to investigate whether there are differences in the values of the MDSS among clinical groups (euthyroid individuals, euthyroid individuals with positive TgAb and/or TPOAb, and hypothyroid and hyperthyroid participants). This cross-sectional study included 4620 participants over 18 years of age from the islands of Korcula and Vis, and the mainland city of Split. The MDSS was assessed from a food frequency questionnaire (FFQ). MDSS values were significantly higher in females compared to males and showed a positive association with the age of the participants. There was no significant difference in the MDSS values among the examined clinical groups. In the group of subjects with euthyroidism, a significant positive association was found between fT3 and the MDSS, while in the group of subjects with subclinical hypothyroidism, a significant positive association was observed between the MDSS and both fT3 and fT4. CT levels were also positively associated with the MDSS. Considering the significant positive association of the MDSS and both fT3 and fT4 levels in patients with subclinical hypothyroidism, the results of this study could be used to create guidelines for selecting an appropriate, potentially protective diet for these patients.

Genetic Landscape of the ACE2 Coronavirus Receptor.

BACKGROUND: SARS-CoV-2, the causal agent of COVID-19, enters human cells using the ACE2 (angiotensin-converting enzyme 2) protein as a receptor. ACE2 is thus key to the infection and treatment of the coronavirus. ACE2 is highly expressed in the heart and respiratory and gastrointestinal tracts, playing important regulatory roles in the cardiovascular and other biological systems. However, the genetic basis of the ACE2 protein levels is not well understood. METHODS: We have conducted the largest genome-wide association meta-analysis of plasma ACE2 levels in >28 000 individuals of the SCALLOP Consortium (Systematic and Combined Analysis of Olink Proteins). We summarize the cross-sectional epidemiological correlates of circulating ACE2. Using the summary statistics-based high-definition likelihood method, we estimate relevant genetic correlations with cardiometabolic phenotypes, COVID-19, and other human complex traits and diseases. We perform causal inference of soluble ACE2 on vascular disease outcomes and COVID-19 severity using mendelian randomization. We also perform in silico functional analysis by integrating with other types of omics data. RESULTS: We identified 10 loci, including 8 novel, capturing 30% of the heritability of the protein. We detected that plasma ACE2 was genetically correlated with vascular diseases, severe COVID-19, and a wide range of human complex diseases and medications. An X-chromosome cis-protein quantitative trait loci-based mendelian randomization analysis suggested a causal effect of elevated ACE2 levels on COVID-19 severity (odds ratio, 1.63 [95% CI, 1.10-2.42]; P=0.01), hospitalization (odds ratio, 1.52 [95% CI, 1.05-2.21]; P=0.03), and infection (odds ratio, 1.60 [95% CI, 1.08-2.37]; P=0.02). Tissue- and cell type-specific transcriptomic and epigenomic analysis revealed that the ACE2 regulatory variants were enriched for DNA methylation sites in blood immune cells. CONCLUSIONS: Human plasma ACE2 shares a genetic basis with cardiovascular disease, COVID-19, and other related diseases. The genetic architecture of the ACE2 protein is mapped, providing a useful resource for further biological and clinical studies on this coronavirus receptor.

Meta-GWAS Reveals Novel Genetic Variants Associated with Urinary Excretion of Uromodulin.

BACKGROUND: Uromodulin, the most abundant protein excreted in normal urine, plays major roles in kidney physiology and disease. The mechanisms regulating the urinary excretion of uromodulin remain essentially unknown. METHODS: We conducted a meta-analysis of genome-wide association studies for raw (uUMOD) and indexed to creatinine (uUCR) urinary levels of uromodulin in 29,315 individuals of European ancestry from 13 cohorts. We tested the distribution of candidate genes in kidney segments and investigated the effects of keratin-40 (KRT40) on uromodulin processing. RESULTS: Two genome-wide significant signals were identified for uUMOD: a novel locus (P 1.24E-08) over the KRT40 gene coding for KRT40, a type 1 keratin expressed in the kidney, and the UMOD-PDILT locus (P 2.17E-88), with two independent sets of single nucleotide polymorphisms spread over UMOD and PDILT. Two genome-wide significant signals for uUCR were identified at the UMOD-PDILT locus and at the novel WDR72 locus previously associated with kidney function. The effect sizes for rs8067385, the index single nucleotide polymorphism in the KRT40 locus, were similar for both uUMOD and uUCR. KRT40 colocalized with uromodulin and modulating its expression in thick ascending limb (TAL) cells affected uromodulin processing and excretion. CONCLUSIONS: Common variants in KRT40, WDR72, UMOD, and PDILT associate with the levels of uromodulin in urine. The expression of KRT40 affects uromodulin processing in TAL cells. These results, although limited by lack of replication, provide insights into the biology of uromodulin, the role of keratins in the kidney, and the influence of the UMOD-PDILT locus on kidney function.

Genome-Wide Meta-Analysis Identifies Multiple Novel Rare Variants to Predict Common Human Infectious Diseases Risk.

Infectious diseases still threaten global human health, and host genetic factors have been indicated as determining risk factors for observed variations in disease susceptibility, severity, and outcome. We performed a genome-wide meta-analysis on 4624 subjects from the 10,001 Dalmatians cohort, with 14 infection-related traits. Despite a rather small number of cases in some instances, we detected 29 infection-related genetic associations, mostly belonging to rare variants. Notably, the list included the genes CD28, INPP5D, ITPKB, MACROD2, and RSF1, all of which have known roles in the immune response. Expanding our knowledge on rare variants could contribute to the development of genetic panels that could assist in predicting an individual's life-long susceptibility to major infectious diseases. In addition, longitudinal biobanks are an interesting source of information for identifying the host genetic variants involved in infectious disease susceptibility and severity. Since infectious diseases continue to act as a selective pressure on our genomes, there is a constant need for a large consortium of biobanks with access to genetic and environmental data to further elucidate the complex mechanisms behind host-pathogen interactions and infectious disease susceptibility.

Role of ST6GAL1 in Thyroid Cancers: Insights from Tissue Analysis and Genomic Datasets.

Thyroid cancer is the predominant endocrine-related malignancy. ST6 ß-galactoside α2,6-sialyltransferase 1 (ST6GAL1) has been studied in various types of cancers; however, the expression and function of ST6GAL1 in thyroid cancer has not been investigated so far. Previously, we conducted two genome-wide association studies and have identified the association of the ST6GAL1 gene with plasma thyroglobulin (Tg) levels. Since Tg levels are altered in thyroid pathologies, in the current study, we wanted to evaluate the expression of ST6GAL1 in thyroid cancer tissues. We performed an immunohistochemical analysis using human thyroid tissue from 89 patients and analyzed ST6GAL1 protein expression in papillary thyroid cancer (including follicular variant and microcarcinoma) and follicular thyroid cancer in comparison to normal thyroid tissue. Additionally, ST6GAL1 mRNA levels from The Cancer Genome Atlas (TCGA, n = 572) and the Genotype-Tissue Expression (GTEx) project (n = 279) were examined. The immunohistochemical analysis revealed higher ST6GAL1 protein expression in all thyroid tumors compared to normal thyroid tissue. TCGA data revealed increased ST6GAL1 mRNA levels in both primary and metastatic tumors versus controls. Notably, the follicular variant of papillary thyroid cancer exhibited significantly higher ST6GAL1 mRNA levels than classic papillary thyroid cancer. High ST6GAL1 mRNA levels significantly correlated with lymph node metastasis status, clinical stage, and reduced survival rate. ST6GAL1 emerges as a potential cancer-associated glycosyltransferase in thyroid malignancies, offering valuable insights into its diagnostic and prognostic significance.

Modeling the COVID-19 epidemic in Croatia: a comparison of three analytic approaches.

AIM: To facilitate the development of a COVID-19 predictive model in Croatia by analyzing three different methodological approaches. METHOD: We used the historical data to explore the fit of the extended SEIRD compartmental model, the Heidler function, an exponential approximation in analyzing electromagnetic phenomena related to lightning strikes, and the Holt-Winters smoothing (HWS) for short-term epidemic predictions. We also compared various methods for the estimation of R0. RESULTS: The R0 estimates for Croatia varied from 2.09 (95% CI 1.77-2.40) obtained by using an empirical post-hoc method to 2.28 (95% CI 2.27-2.28) when we assumed an exponential outbreak at the very beginning of the COVID-19 epidemic in Croatia. Although the SEIRD model provided a good fit for the early epidemic stages, it was outperformed by the Heidler function fit. HWS achieved accurate short-term predictions and depended the least on model entry parameters. Neither model performed well across the entire observed period, which was characterized by multiple wave-form events, influenced by the re-opening for the tourist season during the summer, mandatory masks use in closed spaces, and numerous measures introduced in retail stores and public places. However, an extension of the Heidler function achieved the best overall fit. CONCLUSIONS: Predicting future epidemic events remains difficult because modeling relies on the accuracy of the information on population structure and micro-environmental exposures, constant changes of the input parameters, varying societal adherence to anti-epidemic measures, and changes in the biological interactions of the virus and hosts.

Genome-Wide Association Analysis and Genomic Prediction of Thyroglobulin Plasma Levels.

Thyroglobulin (Tg) is an iodoglycoprotein produced by thyroid follicular cells which acts as an essential substrate for thyroid hormone synthesis. To date, only one genome-wide association study (GWAS) of plasma Tg levels has been performed by our research group. Utilizing recent advancements in computation and modeling, we apply a Bayesian approach to the probabilistic inference of the genetic architecture of Tg. We fitted a Bayesian sparse linear mixed model (BSLMM) and a frequentist linear mixed model (LMM) of 7,289,083 variants in 1096 healthy European-ancestry participants of the Croatian Biobank. Meta-analysis with two independent cohorts (total n = 2109) identified 83 genome-wide significant single nucleotide polymorphisms (SNPs) within the ST6GAL1 gene (p<5×10-8). BSLMM revealed additional association signals on chromosomes 1, 8, 10, and 14. For ST6GAL1 and the newly uncovered genes, we provide physiological and pathophysiological explanations of how their expression could be associated with variations in plasma Tg levels. We found that the SNP-heritability of Tg is 17% and that 52% of this variation is due to a small number of 16 variants that have a major effect on Tg levels. Our results suggest that the genetic architecture of plasma Tg is not polygenic, but influenced by a few genes with major effects.

N-Glycosylation Patterns across the Age-Related Macular Degeneration Spectrum.

The pathogenesis of age-related macular degeneration (AMD) remains elusive, despite numerous research studies. Therefore, we aimed to investigate the changes of plasma and IgG-specific N-glycosylation across the disease severity spectrum. We examined 2835 subjects from the 10.001 Dalmatians project, originating from the isolated Croatian islands of Vis and Korcula. All subjects were classified into four groups, namely (i) bilateral AMD, (ii) unilateral AMD, (iii) early-onset drusen, and (iv) controls. We analysed plasma and IgG N-glycans measured by HPLC and their association with retinal fundus photographs. There were 106 (3.7%) detected cases of AMD; 66 of them were bilateral. In addition, 45 (0.9%) subjects were recorded as having early-onset retinal drusen. We detected several interesting differences across the analysed groups, suggesting that N-glycans can be used as a biomarker for AMD. Multivariate analysis suggested a significant decrease in the immunomodulatory bi-antennary glycan structures in unilateral AMD (adjusted odds ratio 0.43 (95% confidence interval 0.22-0.79)). We also detected a substantial increase in the pro-inflammatory tetra-antennary plasma glycans in bilateral AMD (7.90 (2.94-20.95)). Notably, some of these associations were not identified in the aggregated analysis, where all three disease stages were collapsed into a single category, suggesting the need for better-refined phenotypes and the use of disease severity stages in the analysis of more complex diseases. Age-related macular degeneration progression is characterised by the complex interplay of various mechanisms, some of which can be detected by measuring plasma and IgG N-glycans. As opposed to a simple case-control study, more advanced and refined study designs are needed to understand the pathogenesis of complex diseases.

ACEF performed better than other risk scores in non-ST-elevation acute coronary syndrome during long term follow-up.

BACKGROUND: Risk stratification of patients with non-ST-elevation acute coronary syndrome (NSTE-ACS) is an important clinical method, but long-term studies on patients subjected to all-treatment strategies are lacking. Therefore, the aim was to compare several established risk scores in the all-treatment NSTE-ACS cohort during long-term follow-up. METHODS: Consecutive patients (n = 276) with NSTE-ACS undergoing coronary angiography were recruited between September 2012 and May 2015. Six risk scores for all patients were calculated, namely GRACE 2.0, ACEF, SYNTAX, Clinical SYNTAX, SYNTAX II PCI and SYNTAX II CABG. The primary end-point was Major Adverse Cardiovascular Events (MACE) which was a composite of cardiac death, nonfatal myocardial infarction, ischemic stroke or urgent coronary revascularization. RESULTS: During a median follow-up of 33 months, 64 MACE outcomes were recorded (23.2%). There was no difference between risk score categories, except in the highest risk group of ACEF and SYNTAX II PCI scores which exhibited significantly more MACE (51.6%, N = 33 and 45.3%, N = 29, P = 0.024, respectively). In the multivariate Cox regression analysis of individual variables, only age and atrial fibrillation were significant predictors for MACE (HR 1.03, 95% CI 1.00-1.05, P = 0.023 and HR 2.02, 95% CI 1.04-3.89, P = 0.037, respectively). Furthermore, multivariate analysis of the risk scores showed significant prediction of MACE only with ACEF score (HR 2.16, 95% CI 1.36-3.44, P = 0.001). The overall performance of GRACE, SYNTAX, Clinical SYNTAX and SYNTAX II CABG was poor with AUC values of 0.596, 0.507, 0.530 and 0.582, respectively, while ACEF and SYNTAX II PCI showed the best absolute AUC values for MACE (0.630 and 0.626, respectively). CONCLUSIONS: ACEF risk score showed better discrimination than other risk scores in NSTE-ACS patients undergoing all-treatment strategies over long-term follow-up and it could represent a fast and user-friendly tool to stratify NSTE-ACS patients.

Genome-wide meta-analysis of 158,000 individuals of European ancestry identifies three loci associated with chronic back pain.

Back pain is the #1 cause of years lived with disability worldwide, yet surprisingly little is known regarding the biology underlying this symptom. We conducted a genome-wide association study (GWAS) meta-analysis of chronic back pain (CBP). Adults of European ancestry were included from 15 cohorts in the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) consortium, and from the UK Biobank interim data release. CBP cases were defined as those reporting back pain present for ≥3-6 months; non-cases were included as comparisons ("controls"). Each cohort conducted genotyping using commercially available arrays followed by imputation. GWAS used logistic regression models with additive genetic effects, adjusting for age, sex, study-specific covariates, and population substructure. The threshold for genome-wide significance in the fixed-effect inverse-variance weighted meta-analysis was p<5×10(-8). Suggestive (p<5×10(-7)) and genome-wide significant (p<5×10(-8)) variants were carried forward for replication or further investigation in the remaining UK Biobank participants not included in the discovery sample. The discovery sample comprised 158,025 individuals, including 29,531 CBP cases. A genome-wide significant association was found for the intronic variant rs12310519 in SOX5 (OR 1.08, p = 7.2×10(-10)). This was subsequently replicated in 283,752 UK Biobank participants not included in the discovery sample, including 50,915 cases (OR 1.06, p = 5.3×10(-11)), and exceeded genome-wide significance in joint meta-analysis (OR 1.07, p = 4.5×10(-19)). We found suggestive associations at three other loci in the discovery sample, two of which exceeded genome-wide significance in joint meta-analysis: an intergenic variant, rs7833174, located between CCDC26 and GSDMC (OR 1.05, p = 4.4×10(-13)), and an intronic variant, rs4384683, in DCC (OR 0.97, p = 2.4×10(-10)). In this first reported meta-analysis of GWAS for CBP, we identified and replicated a genetic locus associated with CBP (SOX5). We also identified 2 other loci that reached genome-wide significance in a 2-stage joint meta-analysis (CCDC26/GSDMC and DCC).

Supplementation With the Sialic Acid Precursor N-Acetyl-D-Mannosamine Breaks the Link Between Obesity and Hypertension.

BACKGROUND: Obesity-related hypertension is a common disorder, and attempts to combat the underlying obesity are often unsuccessful. We previously revealed that mice globally deficient in the inhibitory immunoglobulin G (IgG) receptor FcγRIIB are protected from obesity-induced hypertension. However, how FcγRIIB participates is unknown. Studies were designed to determine if alterations in IgG contribute to the pathogenesis of obesity-induced hypertension. METHODS: Involvement of IgG was studied using IgG µ heavy chain-null mice deficient in mature B cells and by IgG transfer. Participation of FcγRIIB was interrogated in mice with global or endothelial cell-specific deletion of the receptor. Obesity was induced by high-fat diet (HFD), and blood pressure (BP) was measured by radiotelemetry or tail cuff. The relative sialylation of the Fc glycan on mouse IgG, which influences IgG activation of Fc receptors, was evaluated by Sambucus nigra lectin blotting. Effects of IgG on endothelial NO synthase were assessed in human aortic endothelial cells. IgG Fc glycan sialylation was interrogated in 3442 human participants by mass spectrometry, and the relationship between sialylation and BP was evaluated. Effects of normalizing IgG sialylation were determined in HFD-fed mice administered the sialic acid precursor N-acetyl-D-mannosamine (ManNAc). RESULTS: Mice deficient in B cells were protected from obesity-induced hypertension. Compared with IgG from control chow-fed mice, IgG from HFD-fed mice was hyposialylated, and it raised BP when transferred to recipients lacking IgG; the hypertensive response was absent if recipients were FcγRIIB-deficient. Neuraminidase-treated IgG lacking the Fc glycan terminal sialic acid also raised BP. In cultured endothelial cells, via FcγRIIB, IgG from HFD-fed mice and neuraminidase-treated IgG inhibited vascular endothelial growth factor activation of endothelial NO synthase by altering endothelial NO synthase phosphorylation. In humans, obesity was associated with lower IgG sialylation, and systolic BP was inversely related to IgG sialylation. Mice deficient in FcγRIIB in endothelium were protected from obesity-induced hypertension. Furthermore, in HFD-fed mice, ManNAc normalized IgG sialylation and prevented obesity-induced hypertension. CONCLUSIONS: Hyposialylated IgG and FcγRIIB in endothelium are critically involved in obesity-induced hypertension in mice, and supportive evidence was obtained in humans. Interventions targeting these mechanisms, such as ManNAc supplementation, may provide novel means to break the link between obesity and hypertension.

Genome-Wide Association Transethnic Meta-Analyses Identifies Novel Associations Regulating Coagulation Factor VIII and von Willebrand Factor Plasma Levels.

BACKGROUND: Factor VIII (FVIII) and its carrier protein von Willebrand factor (VWF) are associated with risk of arterial and venous thrombosis and with hemorrhagic disorders. We aimed to identify and functionally test novel genetic associations regulating plasma FVIII and VWF. METHODS: We meta-analyzed genome-wide association results from 46 354 individuals of European, African, East Asian, and Hispanic ancestry. All studies performed linear regression analysis using an additive genetic model and associated ≈35 million imputed variants with natural log-transformed phenotype levels. In vitro gene silencing in cultured endothelial cells was performed for candidate genes to provide additional evidence on association and function. Two-sample Mendelian randomization analyses were applied to test the causal role of FVIII and VWF plasma levels on the risk of arterial and venous thrombotic events. RESULTS: We identified 13 novel genome-wide significant ( P≤2.5×10-8) associations, 7 with FVIII levels ( FCHO2/TMEM171/TNPO1, HLA, SOX17/RP1, LINC00583/NFIB, RAB5C-KAT2A, RPL3/TAB1/SYNGR1, and ARSA) and 11 with VWF levels ( PDHB/PXK/KCTD6, SLC39A8, FCHO2/TMEM171/TNPO1, HLA, GIMAP7/GIMAP4, OR13C5/NIPSNAP, DAB2IP, C2CD4B, RAB5C-KAT2A, TAB1/SYNGR1, and ARSA), beyond 10 previously reported associations with these phenotypes. Functional validation provided further evidence of association for all loci on VWF except ARSA and DAB2IP. Mendelian randomization suggested causal effects of plasma FVIII activity levels on venous thrombosis and coronary artery disease risk and plasma VWF levels on ischemic stroke risk. CONCLUSIONS: The meta-analysis identified 13 novel genetic loci regulating FVIII and VWF plasma levels, 10 of which we validated functionally. We provide some evidence for a causal role of these proteins in thrombotic events.

Novel Genetic Locus Influencing Retinal Venular Tortuosity Is Also Associated With Risk of Coronary Artery Disease.

OBJECTIVE: The retina may provide readily accessible imaging biomarkers of global cardiovascular health. Increasing evidence suggests variation in retinal vascular traits is highly heritable. This study aimed to identify the genetic determinants of retinal vascular traits. Approach and Results: We conducted a meta-analysis of genome-wide association studies for quantitative retinal vascular traits derived using semi-automatic image analysis of digital retinal photographs from the GoDARTS (Genetics of Diabetes Audit and Research in Tayside; N=1736) and ORCADES (Orkney Complex Disease Study; N=1358) cohorts. We identified a novel genome-wide significant locus at 19q13 (ACTN4/CAPN12) for retinal venular tortuosity (TortV), and one at 13q34 (COL4A2) for retinal arteriolar tortuosity (TortA); these 2 loci were subsequently confirmed in 3 independent cohorts (Ntotal=1413). In the combined analysis of discovery and replication cohorts, the lead single-nucleotide polymorphism in ACTN4/CAPN12 was rs1808382 (ßs.d.=-0.109; SE=0.015; P=2.39×10-13) and in COL4A2 was rs7991229 (ßs.d.=0.103; SE=0.015; P=4.66×10-12). Notably, the ACTN4/CAPN12 locus associated with TortV is also associated with coronary artery disease, heart rate, and atrial fibrillation. CONCLUSIONS: Genetic determinants of retinal vascular tortuosity are also linked to cardiovascular health. These findings provide a molecular pathophysiological foundation for the use of retinal vascular traits as biomarkers for cardiovascular diseases.

Genome-wide meta-analysis identifies novel gender specific loci associated with thyroid antibodies level in Croatians.

Autoimmune thyroid diseases (AITD) are multifactorial endocrine diseases most frequently accompanied by Tg and TPO autoantibodies. Both antibodies have a higher prevalence in females and act under a strong genetic influence. To identify novel variants underlying thyroid antibody levels, we performed GWAS meta-analysis on the plasma levels of TgAb and TPOAb in three Croatian cohorts, as well as gender specific GWAS and a bivariate analysis. No significant association was detected with the level of TgAb and TPOAb in the meta-analysis of GWAS or bivariate results for all individuals. The bivariate analysis in females only revealed a genome-wide significant association for the locus near GRIN3A (rs4457391, P = 7.76 × 10-9). The same locus had borderline association with TPOAb levels in females (rs1935377, P = 8.58 × 10-8). In conclusion, we identified a novel gender specific locus associated with TgAb and TPOAb levels. Our findings provide a novel insight into genetic and gender differences associated with thyroid antibodies.

Coronavirus epidemic in Croatia: case fatality decline during summer?

AIM: To describe the SARS-CoV-2 epidemic pattern in Croatia during February-September 2020 and compare the case fatality ratio (CFR) between spring and summer. METHODS: National data were used to calculate the weekly and monthly CFRs, stratified by three age groups: 0-64, 65-79, and 80+ years. We also calculated the standardized mortality ratios (SMR) to offset the differences in age composition. RESULTS: The epidemic consisted of the initial wave, a trough in June, and two conjoined summer waves, yielding 17206 coronavirus disease 2019 cases and 290 deaths. While the number of confirmed cases nearly quadrupled during summer, case fatality estimates decreased; CFR in spring was 4.81 (95% confidence interval 3.91-5.71), compared with 1.24 (1.06-1.42) in summer. The SMR for summer was 0.45 (0.37-0.55), suggesting that the case fatality risk halved compared with spring. Cardiovascular comorbidity was an important risk factor for case fatality (SMR 2.63 [2.20-3.13] during spring and 1.28 [1.02-1.59] during summer). The risk of death in ventilated patients remained unchanged (SMR 0.98 [0.77-1.24]). CONCLUSIONS: The epidemic dynamics suggests summer decline in case fatality, except in ventilated patients. While the effect of comorbidity also decreased, cardiovascular comorbidity remained an important risk factor for death even during summer. A plethora of possible confounders and an ever-changing landscape of SARS-CoV-2 epidemic in Croatia require constant monitoring and evaluation, with an aim to prevent the uncontrolled spread of the virus and a disruption of health care functioning.

Thyroid hormone levels are associated with metabolic components: a cross-sectional study.

AIM: To analyze the association of thyroid function and hormone levels with metabolic syndrome (MetS) and its components. METHODS: This cross-sectional population-based study involved 2183 Croatian individuals with no history of thyroid disease, hypertension, diabetes, and hyperlipidemia. MetS was diagnosed according to the National Cholesterol Education Program's Adult Treatment Panel III criteria. RESULTS: We found no association between thyroid function groups and the prevalence of MetS and its components. Clinically hypothyroid participants showed significantly higher triceps skinfold measurements than subclinically hypothyroid and euthyroid participants. Furthermore, clinically hypothyroid participants had higher abdominal skinfold thickness than subclinically hypothyroid participants. Otherwise, suprailiac and abdominal skinfold measurements were higher in the subclinically and clinically hyperthyroid group of participants compared with euthyroid and subclinically hypothyroid participants. A strong positive association of thyroid-stimulating hormone (TSH) and strong negative association of free triiodothyronine (fT3) and free thyroxine (fT4) levels with HOMA-IR and cholesterol levels were found. Furthermore, the fT4 level also showed a strong negative association with HDL and triceps skinfold thickness. CONCLUSIONS: This study supports the standing that TSH, fT3, and fT4 levels are important variables to determine the association of thyroid function with MetS.