ABSTRACT
Pubertal timing varies considerably and is associated with later health outcomes. We performed multi-ancestry genetic analyses on ~800,000 women, identifying 1,080 signals for age at menarche. Collectively, these explained 11% of trait variance in an independent sample. Women at the top and bottom 1% of polygenic risk exhibited ~11 and ~14-fold higher risks of delayed and precocious puberty, respectively. We identified several genes harboring rare loss-of-function variants in ~200,000 women, including variants in ZNF483, which abolished the impact of polygenic risk. Variant-to-gene mapping approaches and mouse gonadotropin-releasing hormone neuron RNA sequencing implicated 665 genes, including an uncharacterized G-protein-coupled receptor, GPR83, which amplified the signaling of MC3R, a key nutritional sensor. Shared signals with menopause timing at genes involved in DNA damage response suggest that the ovarian reserve might signal centrally to trigger puberty. We also highlight body size-dependent and independent mechanisms that potentially link reproductive timing to later life disease.
Subject(s)
Gene Frequency , Menarche , Puberty , Humans , Female , Menarche/genetics , Puberty/genetics , Animals , Multifactorial Inheritance/genetics , Mice , Genome-Wide Association Study , Adolescent , Puberty, Precocious/genetics , Polymorphism, Single Nucleotide , Receptors, G-Protein-Coupled/genetics , Puberty, Delayed/genetics , ChildABSTRACT
Calcitonin (CT), a hormone secreted by thyroid parafollicular C cells, plays a role in calcium homeostasis and bone health. Understanding the relationship between CT levels and dietary, sociodemographic, and lifestyle factors is essential for public health and hormonal balance studies. This study encompassed 3323 healthy participants from the Croatian biobank. We utilized principal component analysis (PCA) to reduce food items into dietary patterns. Regression analysis was used to investigate the relationship between CT levels and data collected through questionnaires, accounting for age and sex. CT levels exhibited sex-specific differences, with higher values observed in males. Positive associations were found between CT levels and age, body mass index (BMI), as well as weekly consumption of white and red wine mixed with water. While height and sternal notch-finger length initially correlated positively with CT levels, this relationship reversed upon adjusting for age and sex. Regarding sport activities, CT levels were significantly increased in non-participants compared to occasional sport participants (p = 0.043). Dietary factors yielded intriguing findings, with frequent consumption of butter, animal fat and veal associated with lower CT levels, while higher CT levels were associated with the frequent consumption of white fish, blue fish, pasta, and rice. However, no significant correlation was found between CT levels and bone mineral density (BMD), weight, or body surface area (BSA). This study highlights the complex interplay of dietary, lifestyle, and sociodemographic factors influencing CT levels. These findings suggest that a broad range of factors should be considered in hormonal balance studies, underlining their potential implications for public health.
Subject(s)
Calcitonin , Diet , Life Style , Humans , Male , Female , Middle Aged , Calcitonin/blood , Adult , Aged , Sociodemographic Factors , Bone Density/physiologyABSTRACT
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.
Subject(s)
Diet, Mediterranean , Thyroglobulin , Thyroid Gland , Humans , Female , Male , Thyroid Gland/metabolism , Middle Aged , Adult , Cross-Sectional Studies , Thyroglobulin/blood , Autoantibodies/blood , Autoantibodies/immunology , Aged , Thyrotropin/blood , Triiodothyronine/blood , Hypothyroidism/blood , Thyroid Hormones/blood , Thyroxine/bloodABSTRACT
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.
Subject(s)
Orosomucoid , White People , Adult , Aged , Female , Humans , Male , Middle Aged , Croatia , Glycosylation , Orosomucoid/metabolismABSTRACT
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.
Subject(s)
Cell Adhesion Molecules , Factor VIII , Kininogens , Lectins, C-Type , Receptors, Cell Surface , von Willebrand Factor , Humans , von Willebrand Factor/genetics , von Willebrand Factor/metabolism , Factor VIII/genetics , Factor VIII/metabolism , Polymorphism, Single Nucleotide , Human Umbilical Vein Endothelial Cells/metabolism , Mendelian Randomization Analysis , Genome-Wide Association Study , Thrombosis/genetics , Thrombosis/blood , Genetic Association Studies , Male , Endothelial Cells/metabolism , FemaleABSTRACT
Background: In this priority-setting exercise, we sought to identify leading research priorities needed for strengthening future pandemic preparedness and response across countries. Methods: The International Society of Global Health (ISoGH) used the Child Health and Nutrition Research Initiative (CHNRI) method to identify research priorities for future pandemic preparedness. Eighty experts in global health, translational and clinical research identified 163 research ideas, of which 42 experts then scored based on five pre-defined criteria. We calculated intermediate criterion-specific scores and overall research priority scores from the mean of individual scores for each research idea. We used a bootstrap (n = 1000) to compute the 95% confidence intervals. Results: Key priorities included strengthening health systems, rapid vaccine and treatment production, improving international cooperation, and enhancing surveillance efficiency. Other priorities included learning from the coronavirus disease 2019 (COVID-19) pandemic, managing supply chains, identifying planning gaps, and promoting equitable interventions. We compared this CHNRI-based outcome with the 14 research priorities generated and ranked by ChatGPT, encountering both striking similarities and clear differences. Conclusions: Priority setting processes based on human crowdsourcing - such as the CHNRI method - and the output provided by ChatGPT are both valuable, as they complement and strengthen each other. The priorities identified by ChatGPT were more grounded in theory, while those identified by CHNRI were guided by recent practical experiences. Addressing these priorities, along with improvements in health planning, equitable community-based interventions, and the capacity of primary health care, is vital for better pandemic preparedness and response in many settings.
Subject(s)
COVID-19 , Pandemic Preparedness , Child , Humans , Consensus , Research Design , COVID-19/epidemiology , COVID-19/prevention & control , Child HealthABSTRACT
BACKGROUND: Uncovering the functional relevance underlying verbal declarative memory (VDM) genome-wide association study (GWAS) results may facilitate the development of interventions to reduce age-related memory decline and dementia. METHODS: We performed multi-omics and pathway enrichment analyses of paragraph (PAR-dr) and word list (WL-dr) delayed recall GWAS from 29,076 older non-demented individuals of European descent. We assessed the relationship between single-variant associations and expression quantitative trait loci (eQTLs) in 44 tissues and methylation quantitative trait loci (meQTLs) in the hippocampus. We determined the relationship between gene associations and transcript levels in 53 tissues, annotation as immune genes, and regulation by transcription factors (TFs) and microRNAs. To identify significant pathways, gene set enrichment was tested in each cohort and meta-analyzed across cohorts. Analyses of differential expression in brain tissues were conducted for pathway component genes. RESULTS: The single-variant associations of VDM showed significant linkage disequilibrium (LD) with eQTLs across all tissues and meQTLs within the hippocampus. Stronger WL-dr gene associations correlated with reduced expression in four brain tissues, including the hippocampus. More robust PAR-dr and/or WL-dr gene associations were intricately linked with immunity and were influenced by 31 TFs and 2 microRNAs. Six pathways, including type I diabetes, exhibited significant associations with both PAR-dr and WL-dr. These pathways included fifteen MHC genes intricately linked to VDM performance, showing diverse expression patterns based on cognitive status in brain tissues. CONCLUSIONS: VDM genetic associations influence expression regulation via eQTLs and meQTLs. The involvement of TFs, microRNAs, MHC genes, and immune-related pathways contributes to VDM performance in older individuals.
Subject(s)
Genome-Wide Association Study , MicroRNAs , Humans , Aged , Genome-Wide Association Study/methods , Multiomics , Memory , Cognition , Polymorphism, Single Nucleotide/geneticsABSTRACT
Glycans are an essential structural component of immunoglobulin G (IgG) that modulate its structure and function. However, regulatory mechanisms behind this complex posttranslational modification are not well known. Previous genome-wide association studies (GWAS) identified 29 genomic regions involved in regulation of IgG glycosylation, but only a few were functionally validated. One of the key functional features of IgG glycosylation is the addition of galactose (galactosylation), a trait which was shown to be associated with ageing. We performed GWAS of IgG galactosylation (N=13,705) and identified 16 significantly associated loci, indicating that IgG galactosylation is regulated by a complex network of genes that extends beyond the galactosyltransferase enzyme that adds galactose to IgG glycans. Gene prioritization identified 37 candidate genes. Using a recently developed CRISPR/dCas9 system we manipulated gene expression of candidate genes in the in vitro IgG expression system. Upregulation of three genes, EEF1A1, MANBA and TNFRSF13B, changed the IgG glycome composition, which confirmed that these three genes are involved in IgG galactosylation in this in vitro expression system.
Subject(s)
Galactose , Genome-Wide Association Study , Gene Regulatory Networks , Immunoglobulin G/genetics , Polysaccharides/metabolismABSTRACT
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.
Subject(s)
Genome-Wide Association Study , Thyroid Neoplasms , Humans , Thyroid Cancer, Papillary , Thyroid Neoplasms/genetics , Thyroid Neoplasms/pathology , Genomics , RNA, Messenger/genetics , beta-D-Galactoside alpha 2-6-Sialyltransferase , Antigens, CD/metabolismABSTRACT
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.
Subject(s)
Colitis, Ulcerative , Inflammatory Bowel Diseases , Multiple Sclerosis , Humans , Genome-Wide Association Study , Inflammatory Bowel Diseases/genetics , Quantitative Trait Loci , Colitis, Ulcerative/drug therapy , Colitis, Ulcerative/genetics , Inflammation/genetics , Multiple Sclerosis/genetics , Polymorphism, Single NucleotideABSTRACT
Resting heart rate is associated with cardiovascular diseases and mortality in observational and Mendelian randomization studies. The aims of this study are to extend the number of resting heart rate associated genetic variants and to obtain further insights in resting heart rate biology and its clinical consequences. A genome-wide meta-analysis of 100 studies in up to 835,465 individuals reveals 493 independent genetic variants in 352 loci, including 68 genetic variants outside previously identified resting heart rate associated loci. We prioritize 670 genes and in silico annotations point to their enrichment in cardiomyocytes and provide insights in their ECG signature. Two-sample Mendelian randomization analyses indicate that higher genetically predicted resting heart rate increases risk of dilated cardiomyopathy, but decreases risk of developing atrial fibrillation, ischemic stroke, and cardio-embolic stroke. We do not find evidence for a linear or non-linear genetic association between resting heart rate and all-cause mortality in contrast to our previous Mendelian randomization study. Systematic alteration of key differences between the current and previous Mendelian randomization study indicates that the most likely cause of the discrepancy between these studies arises from false positive findings in previous one-sample MR analyses caused by weak-instrument bias at lower P-value thresholds. The results extend our understanding of resting heart rate biology and give additional insights in its role in cardiovascular disease development.
Subject(s)
Atrial Fibrillation , Cardiovascular Diseases , Humans , Cardiovascular Diseases/genetics , Risk Factors , Heart Rate/genetics , Genetic Predisposition to Disease , Mendelian Randomization Analysis/methods , Genome-Wide Association Study/methods , Polymorphism, Single NucleotideABSTRACT
Genetic studies have identified numerous regions associated with plasma fibrinogen levels in Europeans, yet missing heritability and limited inclusion of non-Europeans necessitates further studies with improved power and sensitivity. Compared with array-based genotyping, whole genome sequencing (WGS) data provides better coverage of the genome and better representation of non-European variants. To better understand the genetic landscape regulating plasma fibrinogen levels, we meta-analyzed WGS data from the NHLBI's Trans-Omics for Precision Medicine (TOPMed) program (n=32,572), with array-based genotype data from the Cohorts for Heart and Aging Research in Genomic Epidemiology (CHARGE) Consortium (n=131,340) imputed to the TOPMed or Haplotype Reference Consortium panel. We identified 18 loci that have not been identified in prior genetic studies of fibrinogen. Of these, four are driven by common variants of small effect with reported MAF at least 10% higher in African populations. Three ( SERPINA1, ZFP36L2 , and TLR10) signals contain predicted deleterious missense variants. Two loci, SOCS3 and HPN , each harbor two conditionally distinct, non-coding variants. The gene region encoding the protein chain subunits ( FGG;FGB;FGA ), contains 7 distinct signals, including one novel signal driven by rs28577061, a variant common (MAF=0.180) in African reference panels but extremely rare (MAF=0.008) in Europeans. Through phenome-wide association studies in the VA Million Veteran Program, we found associations between fibrinogen polygenic risk scores and thrombotic and inflammatory disease phenotypes, including an association with gout. Our findings demonstrate the utility of WGS to augment genetic discovery in diverse populations and offer new insights for putative mechanisms of fibrinogen regulation. Key Points: Largest and most diverse genetic study of plasma fibrinogen identifies 54 regions (18 novel), housing 69 conditionally distinct variants (20 novel).Sufficient power achieved to identify signal driven by African population variant.Links to (1) liver enzyme, blood cell and lipid genetic signals, (2) liver regulatory elements, and (3) thrombotic and inflammatory disease.
ABSTRACT
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.
ABSTRACT
Obstructive sleep apnea (OSA) is a prevalent disease associated with increased risk for cardiovascular and metabolic diseases and shortened lifespan. The aim of this study was to explore the possibility of using N-glycome as a biomarker for the severe form of OSA. Seventy subjects who underwent a whole-night polysomnography/polygraphy and had apnea-hypopnea index (AHI) over 30 were compared to 23 controls (AHI under 5). Plasma samples were used to extract 39 glycan peaks using ultra-high-performance liquid chromatography (UPLC) and 27 IgG peaks using capillary gel electrophoresis (CGE). We also measured glycan age, a molecular proxy for biological aging. Three plasma and one IgG peaks were significant in a multivariate model controlling for the effects of age, sex, and body mass index. These included decreased GP24 (disialylated triantennary glycans as major structure) and GP28 (trigalactosylated, triantennary, disialylated, and trisialylated glycans), and increased GP32 (trisialylated triantennary glycan). Only one IgG glycan peak was significantly increased (P26), which contains biantennary digalactosylated glycans with core fucose. Patients with severe OSA exhibited accelerated biological aging, with a median of 6.9 years more than their chronological age (p < 0.001). Plasma N-glycome can be used as a biomarker for severe OSA.
Subject(s)
Sleep Apnea, Obstructive , Humans , Sleep Apnea, Obstructive/etiology , Polysomnography/adverse effects , Aging , Biomarkers , Immunoglobulin GABSTRACT
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.
Subject(s)
Communicable Diseases , Genetic Predisposition to Disease , Humans , Phenotype , Risk Factors , Genome-Wide Association Study , Communicable Diseases/genetics , Nuclear Proteins/genetics , Trans-Activators/geneticsABSTRACT
Human plasma transferrin (Tf) N-glycosylation has been mostly studied as a marker for congenital disorders of glycosylation, alcohol abuse, and hepatocellular carcinoma. However, inter-individual variability of Tf N-glycosylation is not known, mainly due to technical limitations of Tf isolation in large-scale studies. Here, we present a highly specific robust high-throughput approach for Tf purification from human blood plasma and detailed characterization of Tf N-glycosylation on the level of released glycans by ultra-high-performance liquid chromatography based on hydrophilic interactions and fluorescence detection (HILIC-UHPLC-FLD), exoglycosidase sequencing, and matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF-MS). We perform a large-scale comparative study of Tf and immunoglobulin G (IgG) N-glycosylation analysis in two human populations and demonstrate that Tf N-glycosylation is associated with age and sex, along with multiple biochemical and physiological traits. Observed association patterns differ compared to the IgG N-glycome corroborating tissue-specific N-glycosylation and specific N-glycans' role in their distinct physiological functions.
Subject(s)
Immunoglobulin G , Protein Processing, Post-Translational , Transferrin , Humans , Glycosylation , High-Throughput Screening Assays , Immunoglobulin G/blood , Immunoglobulin G/chemistry , Transferrin/chemistry , Transferrin/isolation & purification , Polysaccharides/analysisABSTRACT
During development, thyroid hormones play an important role in eye development, while in adults, some pathological thyroid conditions can affect the normal functioning of the eyes. Thyroid eye disease is the most well-known eye pathology caused by a pathological thyroid condition. Few studies have investigated the association between ocular parameters and thyroid function. Thus, in this study, we aimed to examine whether thyroid activity affects ocular parameters. This cross-sectional study included 4633 healthy adults recruited within the 10,001 Dalmatians project of the Croatian Biobank. The plasma levels of thyroid-stimulating hormone (TSH), free triiodothyronine (fT3), free thyroxine (fT4), thyroglobulin (Tg), thyroglobulin antibodies (TgAb), and thyroid peroxidase antibodies (TPOAb) were measured by an immunoassay. We determined 20 ocular parameters for each participant (10 for each eye, including corneal radius, corneal thickness, anterior chamber depth, anterior chamber angle, lens thickness, posterior chamber length, axial length, intraocular lens power (IOL), spherical power, and cylinder power). Patients with hyperthyroidism had thicker corneas compared to euthyroid individuals. Corneal thickness was also negatively associated with plasma TSH levels. Intra-ocular lens power was higher in patients with clinical hypothyroidism, while spherical power was higher in euthyroid individuals with positive antibodies compared to euthyroid individuals. Intra-ocular lens power negatively correlated with fT4 levels, while spherical power positively correlated with TgAb, TPOAb, and Tg levels and negatively correlated with TSH levels. The anterior chamber angle was positively associated with plasma TSH levels and TPOAb levels and negatively associated with plasma fT4 levels. These findings suggest an interesting interplay between ophthalmic measures and thyroid status, detectable even in the general adult population.