Gene-environment interaction in type 2 diabetes in Korean cohorts: Interaction of a type 2 diabetes polygenic risk score with triglyceride and cholesterol on fasting glucose levels.

Type 2 diabetes (T2D) is caused by genetic and environmental factors as well as gene-environment interactions. However, these interactions have not been systematically investigated. We analyzed these interactions for T2D and fasting glucose levels in three Korean cohorts, HEXA, KARE, and CAVAS, using the baseline data with a multiple regression model. Two polygenic risk scores for T2D (PRST2D ) and fasting glucose (PRSFG ) were calculated using 488 and 82 single nucleotide polymorphisms (SNP) for T2D and fasting glucose, respectively, which were extracted from large-scaled genome-wide association studies with multiethnic data. Both lifestyle risk factors and T2D-related biochemical measurements were assessed. The effect of interactions between PRST2D -triglyceride (TG) and PRST2D -total cholesterol (TC) on fasting glucose levels was observed as follows: ß ± SE = 0.0005 ± 0.0001, p = 1.06 × 10-19 in HEXA, ß ± SE = 0.0008 ± 0.0001, p = 2.08 × 10-8 in KARE for TG; ß ± SE = 0.0006 ± 0.0001, p = 2.00 × 10-6 in HEXA, ß ± SE = 0.0020 ± 0.0004, p = 2.11 × 10-6 in KARE, ß ± SE = 0.0007 ± 0.0004, p = 0.045 in CAVAS for TC. PRST2D -based classification of the participants into four groups showed that the fasting glucose levels in groups with higher PRST2D were more adversely affected by both the TG and TC. In conclusion, blood TG and TC levels may affect the fasting glucose level through interaction with T2D genetic factors, suggesting the importance of consideration of gene-environment interaction in the preventive medicine of T2D.

Ramlibacter paludis sp. nov., isolated from wetland.

A Gram-stain-negative, non-motile, rod-shaped, aerobic and white-coloured bacterium (designated XY19T) was isolated from a soil sample of wetland from Godeok Ecological Park, Gangdong-gu, Seoul, Republic of Korea. On the basis of 16S rRNA gene sequencing, strain XY19T clustered with species of the genus Ramlibacter and appeared closely related to R. ginsenosidimutans DSM 23480T (98.42 %), R. alkalitolerans JCM 32081T (97.68 %) and R. monticola JCM 31918T (97.66 %). The average nucleotide identity between strain XY19T and three strains (R. ginsenosidimutans DSM 23480T, R. alkalitolerans JCM 32081T and R. monticola JCM 31918T) were 80.7, 81.1 and 81.4 %. And the digital DNA-DNA hybridization (dDDH) calculated between strain XY19T and each of the three strains (R. ginsenosidimutans DSM 23480T, R. alkalitolerans JCM 32081T and R. monticola JCM 31918T) were 24.1, 24.4 and 24.5 %. ANI value and dDDH results were a novel species of the genus Ramlibacter. Growth occurs at 10-37 °C on R2A medium in the pressence of 0-1 % NaCl (w/v) and at pH 6.0-8.5. The DNA G+C content of the genomic DNA was 68.7 mol%, and ubiquinone-8 (Q-8) was the major respiratory quinone. The major cellular fatty acids (>5 %) were C16:1 ω7c and/or C16:1 ω6c (summed feature 3), C16 : 0, C17 : 0 cyclo and C18:1 ω7c and/or C18:1 ω6c (summed feature 8). The polar lipids consisted of diphosphatidylglycerol, phosphatidylglycerol, phosphatidylethanolamine, three unidentified lipids and unidentified aminophospholipid. Physiological and biochemical characteristics indicated that strain XY19T represents a novel species of the genus Ramlibacter, for which the name Ramlibacter paludis sp. nov. is proposed. The type strain is XY19T (= KACC 22220T = LMG 32190T).

Genome-wide gene and serum ferritin interaction in the development of type 2 diabetes in adults aged 40 years or older.

BACKGROUND AND AIMS: Elevated serum ferritin is associated with incident Type 2 diabetes (T2D), but the interactions between serum ferritin and genetic factors which may improve understanding underlying mechanism in the development of T2D are still unclear. We determined the gene-ferritin interactions on the development of T2D by genome-wide gene-ferritin interaction analyses. METHODS AND RESULTS: A total of 3405 participants from two prospective cohorts of community living residents were included, and the median follow-time was 3.99 years. Genome-wide gene-ferritin interactions were analyzed using the joint test with two degrees of freedom and the interaction test with one degree of freedom. There were 18 SNPs selected in the joint test. Finally, four independent variants [rs355140 (LINC00312), rs4075576 (nearby PDGFA), rs1332202 (PTPRD), and rs713157 (nearby LINC00900)] with low pairwise linkage disequilibrium (r2<0.2) and located at least 1000 kb from the index SNP showed interactions with serum ferritin level. In the association analyses between serum ferritin levels (tertiles of ferritin and ferritin status) and the incidence of T2D according to genotype, the Incidence Rate Ratios (IRRs) in the highest tertile of ferritin level (vs. the lowest tertile) were greater for participants with heterozygotes of risk alleles of each of the four SNP than IRRs for those with wild type. Compared with the normal group, the elevated ferritin group also had a higher risk of T2D for all genetic variants of risk alleles, particularly its homozygotes. CONCLUSION: Serum ferritin level interacts with genetic variants (rs355140, rs4075576, rs1332202, and rs713157) in the development of T2D.

Genome-wide interaction study of single-nucleotide polymorphisms and alcohol consumption on blood pressure: The Ansan and Ansung study of the Korean Genome and Epidemiology Study (KoGES).

Hypertension is a common disease worldwide. Alcohol consumption is one of the risk factors for hypertension, however, it is unclear how alcohol consumption elevates blood pressure. Blood pressure could be affected by interactions between genetic variations and alcohol consumption. Thus, we performed a genome-wide interaction study (GWIS) to assess the effect of gene-alcohol consumption interaction on blood pressure among adults aged ≥40 years from the Ansan and Ansung cohort study (n = 6,176), a part of the Korean Genome Epidemiology Study (KoGES). As a result, rs1297184, single-nucleotide polymorphism (SNP) in locus LGR5 was significant (PGWIS = 8.78 × 10-9 ) in GWIS analysis on diastolic blood pressure, but not on systolic blood pressure. However, there was a heteroscedasticity of alcohol consumption. In the GWIS analysis, applying the inverse-variance weighting to correct the systematic inflation slightly attenuated the strength of interaction (PGWIS_IVW = 7.14 × 10-8 ). This interaction was replicated in the Health Examinees cohort (p = .026), a large-scale community-based cohort (n = 18,708). In conclusion, we identified a possible novel interaction between an SNP (rs1297184) and alcohol consumption on blood pressure.

Gene-environment interactions related to blood pressure traits in two community-based Korean cohorts.

Hypertension is a complex disorder caused by genetic and environmental risk factors. Recently, genome-wide association studies (GWASs) identified more than 100 genetic variants for blood pressure traits and hypertension. However, the interactions between these genetic variants and environmental factors have not been systematically investigated. Therefore, we examined the interaction between genetic and environmental risk factors in blood pressure traits using the genetic risk score (GRS). Two Korean community-based cohorts, Cohort I (KARE; N = 8,840) and Cohort II (CAVAS; N = 9,599), were used for this study, and GRSs were calculated from 42 GWAS single-nucleotide polymorphisms (SNPs) that were validated for their association in these cohorts. We calculated GRSs in both ways by considering the effect sizes of each SNP (weighted GRS) and not considering the effect sizes (unweighted GRS). The unweighted GRS was strongly associated with systolic blood pressure, diastolic blood pressure, and hypertension (p = 9.03 × 10 -47 , p = 9.41 × 10 -48 , and p = 3.22 × 10 -55 by meta-analysis, respectively) and the weighted GRS showed the similar results. The environmental factors of body mass index, waist circumference, and drinking status were significantly associated with blood pressure traits, and the interaction between these factors and GRSs were examined. However, no interactions were found with either the GRS or the individual SNPs considered for the GRS. Our findings show that it is challenging to find GRS-environment interactions regarding blood pressure traits.

Cardiac-specific inactivation of Prdm16 effects cardiac conduction abnormalities and cardiomyopathy-associated phenotypes.

A variant in the PRDM16 locus has been correlated with QRS duration in an electrocardiogram genome-wide association study, and the deletion of PRDM16 has been implicated as a causal factor of the dilated cardiomyopathy that is linked to 1p36 deletion syndrome. We aimed to determine how a null mutation of Prdm16 affects cardiac function and study the underlying mechanism of the resulting phenotype in an appropriate mouse model. We used cardiac-specific Prdm16 conditional knockout mice to examine cardiac function by electrocardiography. QRS duration and QTc interval increased significantly in cardiac-specific Prdm16 knockout animals compared with wild-type mice. Further, we assessed cardiomyopathy-associated features by trichrome staining, densitometry, and hydroxyproline assay. Prdm16-null hearts showed greater fibrosis and cardiomyocyte hypertrophy. By quantitative real-time PCR, Prdm16-null hearts upregulated extracellular matrix-related genes (Ctgf, Timp1) and α-smooth muscle actin (Acta2), a myofibroblast marker. Moreover, TGF-ß signaling was activated in Prdm16-null hearts, as evidenced by increased Tgfb1-3 transcript levels and phosphorylated Smad2. However, the inhibition of TGF-ß receptor did not reverse the aberrations in conduction in cardiac-specific Prdm16 knockout mice. To determine the underlying mechanisms, we performed RNA-seq using mouse left ventricular tissue. By functional analysis, Prdm16-null hearts experienced dysregulated expression of ion channel genes, including Kcne1, Scn5a, Cacna1h, and Cacna2d2. Mice with Prdm16-null hearts develop abnormalities in cardiac conduction and cardiomyopathy-associated phenotypes, including fibrosis and cellular hypertrophy. Further, the RNA-seq findings suggest that impairments in ion homeostasis (Ca2+, K+, and Na+) may at least partially underlie the abnormal conduction in cardiac-specific Prdm16 knockout mice.NEW & NOTEWORTHY This is the first study that describes aberrant cardiac function and cardiomyopathy-associated phenotypes in an appropriate murine genetic model with cardiomyocyte-specific Prdm16-null mutation. It is noteworthy that the correlation of PRDM16 with QRS duration is replicated in a murine animal model and the potential underlying mechanism may be the impairment of ion homeostasis.

Identification of five novel genetic loci related to facial morphology by genome-wide association studies.

BACKGROUND: Face morphology is strongly determined by genetic factors. However, only a small number of genes related to face morphology have been identified to date. Here, we performed a two-stage genome-wide association study (GWAS) of 85 face morphological traits in 7569 Koreans (5643 in the discovery set and 1926 in the replication set). RESULTS: In this study, we analyzed 85 facial traits, including facial angles. After discovery GWAS, 128 single nucleotide polymorphisms (SNPs) showing an association of P < 5 × 10- 6 were selected to determine the replication of the associations, and meta-analysis of discovery GWAS and the replication analysis resulted in five genome-wide significant loci. The OSR1-WDR35 [rs7567283, G allele, beta (se) = -0.536 (0.096), P = 2.75 × 10- 8] locus was associated with the facial frontal contour; the HOXD1-MTX2 [rs970797, A allele, beta (se) = 0.015 (0.003), P = 3.97 × 10- 9] and WDR27 [rs3736712, C allele, beta (se) = 0.293 (0.048), P = 8.44 × 10- 10] loci were associated with eye shape; and the SOX9 [rs2193054, C allele, beta (se) (ln-transformed) = -0.007 (0.001), P = 6.17 × 10- 17] and DHX35 [rs2206437, A allele, beta (se) = -0.283 (0.047), P = 1.61 × 10- 9] loci were associated with nose shape. WDR35 and SOX9 were related to known craniofacial malformations, i.e., cranioectodermal dysplasia 2 and campomelic dysplasia, respectively. In addition, we found three independent association signals in the SOX9 locus, and six known loci for nose size and shape were replicated in this study population. Interestingly, four SNPs within these five face morphology-related loci showed discrepancies in allele frequencies among ethnic groups. CONCLUSIONS: We identified five novel face morphology loci that were associated with facial frontal contour, nose shape, and eye shape. Our findings provide useful genetic information for the determination of face morphology.

Validation of the new AJCC eighth edition of the TNM classification for breast cancer with a single-center breast cancer cohort.

PURPOSE: The new eighth edition TNM classification by the AJCC for breast cancer (BC) incorporates biologic factors and gene expression prognostic panels, in addition to traditional anatomic factors. In this study, we evaluated the prognostic value of this new staging system compared to the previous AJCC 7th edition staging system. METHODS: We conducted a retrospective analysis of women with stage I, II, or III BC who underwent curative surgery with/without adjuvant systemic therapy at Samsung Medical Center between July 2004 and December 2008. RESULTS: Of 3,208 BCs, this study was analyzed using the information of 2,790 BC patients. Hormone receptor-positive (HR+) and human epidermal growth factor 2 (HER2)- BCs were observed in 62.9% of BCs, HR+/ HER2+ in 9.3%, HR-/HER2- in 17.0%, and HR-/HER2+ in 10.8%. In survival analysis, we observed 245 distant recurrences and 198 deaths caused by BC progression. The median follow-up duration was 116.2 months. 10-year disease-specific survival (DSS) rates according to the AJCC 7th edition criteria were 97.2% of stage IA, 100% of IB, 94.9% of IIA, 87.9% of IIB, 86.4% of IIIA, 95.7% of IIIB, and 65.7% of IIIC (p < 0.001). After applying 8th edition criteria, the 10-year DSS rates were 98.1% of stage IA, 97.7% of IB, 93.8% of IIA, 92.7% of IIB, 88.2% of IIIA, 80.8% of IIIB, and 70.3% of IIIC (p < 0.001). CONCLUSIONS: The AJCC 8th edition clinical staging system provides a good prognostic value and addresses the weakness of the AJCC 7th edition, which uses only anatomical pathologic staging.

GAREM1 regulates the PR interval on electrocardiograms.

PR interval is the period from the onset of P wave to the start of the QRS complex on electrocardiograms. A recent genomewide association study (GWAS) suggested that GAREM1 was linked to the PR interval on electrocardiograms. This study was designed to validate this correlation using additional subjects and examined the function of Garem1 in a mouse model. We analyzed the association of rs17744182, a variant in the GAREM1 locus, with the PR interval in 5646 subjects who were recruited from 2 Korean replication sets, Yangpyeong (n = 2471) and Yonsei (n = 3175), and noted a significant genomewide association by meta-analysis (P = 2.39 × 10-8). To confirm the function of Garem1 in mice, Garem1 siRNA was injected into mouse tail veins to reduce the expression of Garem1. Garem1 transcript levels declined by 53% in the atrium of the heart (P = 0.029), and Garem1-siRNA injected mice experienced a significant decrease in PR interval (43.27 ms vs. 44.89 ms in control, P = 0.007). We analyzed the expression pattern of Garem1 in the heart by immunohistology and observed specific expression of Garem1 in intracardiac ganglia. Garem1 was expressed in most neurons of the ganglion, including cholinergic and adrenergic cells. We have provided evidence that GAREM1 is involved in the PR interval of ECGs. These findings increase our understanding of the regulatory signals of heart rhythm through intracardiac ganglia of the autonomic nervous system and can be used to guide the development of a therapeutic target for heart conditions, such as atrial fibrillation.

Csk Regulates Blood Pressure by Controlling the Synthetic Pathways of Aldosterone.

BACKGROUND: Blood pressure is regulated by a network of diverse physiological pathways. The C-terminal Src kinase (CSK) locus (15q24) is associated with blood pressure in various ethnic groups. It was recently reported thatCskinsufficiency increases blood pressure through Src. The mechanisms of hypertension inCsk+/-mice are examined further in this study.Methods and Results:To identify a causal component responsible for hypertension inCsk+/-, the heart rate was measured by electrocardiogram and plasma volume by Evans blue dilution. Plasma volume increased inCsk+/-compared with wild-types, while the heart rate did not change. Plasma sodium and aldosterone levels rose consistently inCsk+/-vs. wild-types, and spironolactone, a mineralocorticoid receptor antagonist, reduced blood pressure. The amounts of Sgk1 and Na+/K+-ATPase (NKA) increased in the kidney ofCsk+/-compared with wild-types. It was also found that Cyp11b2 (aldosterone synthase) was upregulated in the adrenal glands ofCsk+/-, and that Csk was enriched in the zona glomerulosa of adrenals, the major site of aldosterone production in the normal mouse. CONCLUSIONS: The results of the present study identify a physiological pathway by which blood pressure is regulated, in which the insufficiency ofCskinduces aldosterone production with zonal specificity in the adrenal glands, increasing sodium reabsorption and plasma volume and thus resulting in hypertension.

The Mechanism of p53 Rescue by SUSP4.

p53 is an important tumor-suppressor protein deactivation of which by mdm2 results in cancers. A SUMO-specific protease 4 (SUSP4) was shown to rescue p53 from mdm2-mediated deactivation, but the mechanism is unknown. The discovery by NMR spectroscopy of a "p53 rescue motif" in SUSP4 that disrupts p53-mdm2 binding is presented. This 29-residue motif is pre-populated with two transient helices connected by a hydrophobic linker. The helix at the C-terminus binds to the well-known p53-binding pocket in mdm2 whereas the N-terminal helix serves as an affinity enhancer. The hydrophobic linker binds to a previously unidentified hydrophobic crevice in mdm2. Overall, SUSP4 appears to use two synergizing modules, the p53 rescue motif described here and a globular-structured SUMO-binding catalytic domain, to stabilize p53. A p53 rescue motif peptide exhibits an anti-tumor activity in cancer cell lines expressing wild-type p53. A pre-structures motif in the intrinsically disordered proteins is thus important for target recognition.

Identification of three novel genetic variations associated with electrocardiographic traits (QRS duration and PR interval) in East Asians.

The electrocardiogram has several advantages in detecting cardiac arrhythmia-it is readily available, noninvasive and cost-efficient. Recent genome-wide association studies have identified single-nucleotide polymorphisms that are associated with electrocardiogram measures. We performed a genome-wide association study using Korea Association Resource data for the discovery phase (Phase 1, n = 6805) and two consecutive replication studies in Japanese populations (Phase 2, n = 2285; Phase 3, n = 5010) for QRS duration and PR interval. Three novel loci were identified: rs2483280 (PRDM16 locus) and rs335206 (PRDM6 locus) were associated with QRS duration, and rs17026156 (SLC8A1 locus) correlated with PR interval. PRDM16 was recently identified as a causative gene of left ventricular non-compaction and dilated cardiomyopathy in 1p36 deletion syndrome, which is characterized by heart failure, arrhythmia and sudden cardiac death. Thus, our finding that a PRDM16 SNP is linked to QRS duration strongly implicates PRDM16 in cardiac function. In addition, C allele of rs17026156 increases PR interval (beta ± SE, 2.39 ± 0.40 ms) and exists far more frequently in East Asians (0.46) than in Europeans and Africans (0.05 and 0.08, respectively).

Default Mode Network Functional Connectivity in Early and Late Mild Cognitive Impairment: Results From the Alzheimer's Disease Neuroimaging Initiative.

BACKGROUND: Default mode network (DMN) functional connectivity is one of the neuroimaging candidate biomarkers of Alzheimer disease. However, no studies have investigated DMN connectivity at different stages of mild cognitive impairment (MCI). The aim of this study was to investigate patterns of DMN connectivity and its breakdown among cognitively normal (CN), early MCI (EMCI), and late MCI (LMCI) subjects. METHODS: Magnetic resonance imaging data and neuropsychological test scores from 130 subjects (CN=43, EMCI=47, LMCI=40) were obtained from the Alzheimer's Disease Neuroimaging Initiative. DMN functional connectivity was extracted using independent components analysis and compared between groups. RESULTS: Functional connectivity in the precuneus, bilateral medial frontal, parahippocampal, middle temporal, right superior temporal, and left angular gyri was decreased in EMCI subjects compared with CN subjects. When the 2 MCI groups were directly compared, LMCI subjects exhibited decreased functional connectivity in the precuneus, bilateral medial frontal gyri, and left angular gyrus. There was no significant difference in gray matter volume among the 3 groups. Amyloid-positive EMCI subjects revealed more widespread breakdown of DMN connectivity than amyloid-negative EMCI subjects. A quantitative index of DMN connectivity correlated well with measures of cognitive performance. CONCLUSIONS: Our results suggest that the breakdown of DMN connectivity may occur in the early stage of MCI.

Influence of blood pressure polygenic risk scores and environmental factors on coronary artery disease in the Korean Genome and Epidemiology Study.

The present study aimed to investigate the association of blood pressure polygenic risk scores (BP PRSs) with coronary artery disease (CAD) in a Korean population and the interaction effects between PRSs and environmental factors on CAD. Data were derived from the Cardiovascular Disease Association Study (CAVAS; N = 5100) and the Health Examinee Study (HEXA; N = 58,623) within the Korean Genome and Epidemiology Study. PRSs for systolic and diastolic BP were calculated with the weighted allele sum of >200 single-nucleotide polymorphisms. Multivariable logistic regression models were used. BP PRSs were strongly associated with systolic BP (SBP), diastolic BP (DBP), and hypertension in both CAVAS and HEXA (p < 0.0001). PRSSBP was significantly associated with CAD in CAVAS, while PRSSBP and PRSDBP were significantly associated with CAD in HEXA. There was an interaction effect between the BP PRSs and environmental factors on CAD. The odds ratios (ORs) for CAD were 1.036 (95% confidence interval [CI], 1.016-1.055) for obesity, 1.028 (95% CI, 1.011-1.045) for abdominal obesity, 1.030 (95% CI, 1.009-1.050) for triglyceride, 1.024 (95% CI, 1.008-1.041) for high-density lipoprotein cholesterol, and 1.039 for smoking (95% CI, 1.003-1.077) in CAVAS. There was no significant interaction in HEXA, except between PRSDBP and triglyceride (OR, 1.012; 95% CI, 1.001-1.024). BP PRS was associated with an increased risk of hypertension and CAD. The interactions among PRSs and environmental risk factors increased the risk of CAD. Multi-component interventions to lower BP in the population via healthy behaviors are needed to prevent CAD regardless of genetic predisposition.

Integration of risk factor polygenic risk score with disease polygenic risk score for disease prediction.

Polygenic risk score (PRS) is useful for capturing an individual's genetic susceptibility. However, previous studies have not fully exploited the potential of the risk factor PRS (RFPRS) for disease prediction. We explored the potential of integrating disease-related RFPRSs with disease PRS to enhance disease prediction performance. We constructed 112 RFPRSs and analyzed the association of RFPRSs with diseases to identify disease-related RFPRSs in 700 diseases, using the UK Biobank dataset. We uncovered 6157 statistically significant associations between 247 diseases and 109 RFPRSs. We estimated the disease PRSs of 70 diseases that exhibited statistically significant heritability, to generate RFDiseasemetaPRS-a combined PRS integrating RFPRSs and disease PRS-and compare the prediction performance metrics between RFDiseasemetaPRS and disease PRS. RFDiseasemetaPRS showed better performance for Nagelkerke's pseudo-R2, odds ratio (OR) per 1 SD, net reclassification improvement (NRI) values and difference of R2 considered by variance of R2 in 31 out of 70 diseases. Additionally, we assessed risk classification between two models by examining OR between the top 10% and remaining 90% individuals for the 31 diseases; RFDiseasemetaPRS exhibited better R2, NRI and OR than disease PRS. These findings highlight the importance of utilizing RFDiseasemetaPRS, which can provide personalized healthcare and tailored prevention strategies.

Common variants in RYR1 are associated with left ventricular hypertrophy assessed by electrocardiogram.

AIMS: To identify the genetic risk factors that influence the development of electrocardiographic (ECG) left ventricular hypertrophy (LVH), a major risk factor for cardiovascular (CV) morbidity and mortality. METHODS AND RESULTS: We performed a genomewide association study (GWAS) of ECG-LVH, in which the community-based Korea Association REsource (KARE) study (8432 controls and 398 cases) was analysed by Affymetrix SNP array 5.0. The GWAS results were validated in hospital-based samples (597 controls and 207 cases). Fourteen single-nucleotide polymorphisms (SNPs) in eight genetic loci (5q35.1, 6p22.3-22.1, 8q24.2, 11p15, 11q21-22.1, 14q12, 17q11.2, and 19q13.1) were associated with ECG-LVH in the original GWAS study (P < 1 × 10(-5)). Of these SNPs, 12 were genotyped in the hospital sample. There was consistent association with the 19q13.1 region which contains RYR1 gene. The most significant SNP in the region was rs10500279, which had genomewide significance in the combined GWAS/replication sample [odds ratio = 1.58 (confidence interval: 1.35-1.85), P = 1.0 × 10(-8)]. Mutations in RYR1, which encodes a major Ca(2+) channel in the skeletal muscle, have been reported to correlate with CV diseases. CONCLUSION: We performed the first GWAS for ECG-LVH, implicating the skeletal muscle Ca(2+) channel protein RYR1 as a genetic risk factor. These results might increase our understanding of the development of ECG-LVH.

Association between dyslipidemia and asthma in children: a systematic review and multicenter cohort study using a common data model.

BACKGROUND: The association between dyslipidemia and asthma in children remains unclear. PURPOSE: This study investigated the association between dyslipidemia and cholesterol levels in children. METHODS: A systematic literature review was performed to identify studies investigating the association between dyslipidemia and asthma in children. The PubMed database was searched for articles published from January 2000-March 2022. Data from a cohort study using electronic health records from 5 hospitals, converted to the Observational Medical Outcomes Partnership Common Data Model (OMOP-CDM), were used to identify the association between total cholesterol (TC) levels and asthma in children. This cohort study used the Cox proportional hazards model to examine hazard ratio (HR) of asthma after propensity score matching, and included an aggregate meta-analysis of HR. RESULTS: We examined 11 studies reporting an association between dyslipidemia and asthma in children. Most were cross-sectional; however, their results were inconsistent. In OMOP-CDM multicenter analysis, the high TC (>170 mg/dL) group included 29,038 children, while the normal TC (≤170 mg/dL) group included 88,823 children including all hospital datasets. In a meta-analysis of this multicenter cohort, a significant association was found between high TC levels and later development of asthma in children <15 years of age (pooled HR, 1.30; 95% confidence interval, 1.12-1.52). CONCLUSION: Elevated TC levels in children may be associated with asthma.

Investigation of heteroscedasticity in polygenic risk scores across 15 quantitative traits.

The polygenic risk score (PRS) could be used to stratify individuals with high risk of diseases and predict complex trait of individual in a population. Previous studies developed a PRS-based prediction model using linear regression and evaluated the predictive performance of the model using the R 2 value. One of the key assumptions of linear regression is that the variance of the residual should be constant at each level of the predictor variables, called homoscedasticity. However, some studies show that PRS models exhibit heteroscedasticity between PRS and traits. This study analyzes whether heteroscedasticity exists in PRS models of diverse disease-related traits and, if any, it affects the accuracy of PRS-based prediction in 354,761 Europeans from the UK Biobank. We constructed PRSs for 15 quantitative traits using LDpred2 and estimated the existence of heteroscedasticity between PRSs and 15 traits using three different tests of the Breusch-Pagan (BP) test, score test, and F test. Thirteen out of fifteen traits show significant heteroscedasticity. Further replication using new PRSs from the PGS catalog and independent samples (N = 23,620) from the UK Biobank confirmed the heteroscedasticity in ten traits. As a result, ten out of fifteen quantitative traits show statistically significant heteroscedasticity between the PRS and each trait. There was a greater variance of residuals as PRS increased, and the prediction accuracy at each level of PRS tended to decrease as the variance of residuals increased. In conclusion, heteroscedasticity was frequently observed in the PRS-based prediction models of quantitative traits, and the accuracy of the predictive model may differ according to PRS values. Therefore, prediction models using the PRS should be constructed by considering heteroscedasticity.

Genetic differences according to onset age and lung function in asthma: A cluster analysis.

BACKGROUND: The extent of differences between genetic risks associated with various asthma subtypes is still unknown. To better understand the heterogeneity of asthma, we employed an unsupervised method to identify genetic variants specifically associated with asthma subtypes. Our goal was to gain insight into the genetic basis of asthma. METHODS: In this study, we utilized the UK Biobank dataset to select asthma patients (All asthma, n = 50,517) and controls (n = 283,410). We excluded 14,431 individuals who had no information on predicted values of forced expiratory volume in one second percent (FEV1%) and onset age, resulting in a final total of 36,086 asthma cases. We conducted k-means clustering based on asthma onset age and predicted FEV1% using these samples (n = 36,086). Cluster-specific genome-wide association studies were then performed, and heritability was estimated via linkage disequilibrium score regression. To further investigate the pathophysiology, we conducted eQTL analysis with GTEx and gene-set enrichment analysis with FUMA. RESULTS: Clustering resulted in four distinct clusters: early onset asthmanormalLF (early onset with normal lung function, n = 8172), early onset asthmareducedLF (early onset with reduced lung function, n = 8925), late-onset asthmanormalLF (late-onset with normal lung function, n = 12,481), and late-onset asthmareducedLF (late-onset with reduced lung function, n = 6508). Our GWASs in four clusters and in All asthma sample identified 5 novel loci, 14 novel signals, and 51 cluster-specific signals. Among clusters, early onset asthmanormalLF and late-onset asthmareducedLF were the least correlated (rg  = 0.37). Early onset asthmareducedLF showed the highest heritability explained by common variants (h2  = 0.212) and was associated with the largest number of variants (71 single nucleotide polymorphisms). Further, the pathway analysis conducted through eQTL and gene-set enrichment analysis showed that the worsening of symptoms in early onset asthma correlated with lymphocyte activation, pathogen recognition, cytokine receptor activation, and lymphocyte differentiation. CONCLUSIONS: Our findings suggest that early onset asthmareducedLF was the most genetically predisposed cluster, and that asthma clusters with reduced lung function were genetically distinct from clusters with normal lung function. Our study revealed the genetic variation between clusters that were segmented based on onset age and lung function, providing an important clue for the genetic mechanism of asthma heterogeneity.

Gene-environment interaction explains a part of missing heritability in human body mass index.

Gene-environment (G×E) interaction could partially explain missing heritability in traits; however, the magnitudes of G×E interaction effects remain unclear. Here, we estimate the heritability of G×E interaction for body mass index (BMI) by subjecting genome-wide interaction study data of 331,282 participants in the UK Biobank to linkage disequilibrium score regression (LDSC) and linkage disequilibrium adjusted kinships-software for estimating SNP heritability from summary statistics (LDAK-SumHer) analyses. Among 14 obesity-related lifestyle factors, MET score, pack years of smoking, and alcohol intake frequency significantly interact with genetic factors in both analyses, accounting for the partial variance of BMI. The G×E interaction heritability (%) and standard error of these factors by LDSC and LDAK-SumHer are as follows: MET score, 0.45% (0.12) and 0.65% (0.24); pack years of smoking, 0.52% (0.13) and 0.93% (0.26); and alcohol intake frequency, 0.32% (0.10) and 0.80% (0.17), respectively. Moreover, these three factors are partially validated for their interactions with genetic factors in other obesity-related traits, including waist circumference, hip circumference, waist-to-hip ratio adjusted with BMI, and body fat percentage. Our results suggest that G×E interaction may partly explain the missing heritability in BMI, and two G×E interaction loci identified could help in understanding the genetic architecture of obesity.