Identification of novel susceptibility loci for non‐syndromic cleft lip with or without cleft palate

Although several genome‐wide association studies (GWAS) of non‐syndromic cleft lip with or without cleft palate (NSCL/P) have been reported, more novel association signals are remained to be exploited. Here, we performed an in‐depth analysis of our previously published Chinese GWAS cohort study with replication in an extra dbGaP case‐parent trios and another in‐house Nanjing cohort, and finally identified five novel significant association signals (rs11119445: 3’ of SERTAD4, P = 6.44 × 10⁻¹⁴; rs227227 and rs12561877: intron of SYT14, P = 5.02 × 10⁻¹³ and 2.80 × 10⁻¹¹, respectively; rs643118: intron of TRAF3IP3, P = 4.45 × 10⁻⁶; rs2095293: intron of NR6A1, P = 2.98 × 10⁻⁵). The mean (standard deviation) of the weighted genetic risk score (wGRS) from these SNPs was 1.83 (0.65) for NSCL/P cases and 1.58 (0.68) for controls, respectively (P = 2.67 × 10⁻¹⁶). Rs643118 was identified as a shared susceptible factor of NSCL/P among Asians and Europeans, while rs227227 may contribute to the risk of NSCL/P as well as NSCPO. In addition, sertad4 knockdown zebrafish models resulted in down‐regulation of sox2 and caused oedema around the heart and mandibular deficiency, compared with control embryos. Taken together, this study has improved our understanding of the genetic susceptibility to NSCL/P and provided further clues to its aetiology in the Chinese population.     

Genetic associations for two biological age measures point to distinct aging phenotypes

Biological age measures outperform chronological age in predicting various aging outcomes, yet little is known regarding genetic predisposition. We performed genome‐wide association scans of two age‐adjusted biological age measures (PhenoAgeAcceleration and BioAgeAcceleration), estimated from clinical biochemistry markers (Levine et al., 2018; Levine, 2013) in European‐descent participants from UK Biobank. The strongest signals were found in the APOE gene, tagged by the two major protein‐coding SNPs, PhenoAgeAccel—rs429358 (APOE e4 determinant) (p = 1.50 × 10⁻⁷²); BioAgeAccel—rs7412 (APOE e2 determinant) (p = 3.16 × 10⁻⁶⁰). Interestingly, we observed inverse APOE e2 and e4 associations and unique pathway enrichments when comparing the two biological age measures. Genes associated with BioAgeAccel were enriched in lipid related pathways, while genes associated with PhenoAgeAccel showed enrichment for immune system, cell function, and carbohydrate homeostasis pathways, suggesting the two measures capture different aging domains. Our study reaffirms that aging patterns are heterogeneous across individuals, and the manner in which a person ages may be partly attributed to genetic predisposition.     

Integrative analysis of clinical and epigenetic biomarkers of mortality

DNA methylation (DNAm) has been reported to be associated with many diseases and with mortality. We hypothesized that the integration of DNAm with clinical risk factors would improve mortality prediction. We performed an epigenome‐wide association study of whole blood DNAm in relation to mortality in 15 cohorts (n = 15,013). During a mean follow‐up of 10 years, there were 4314 deaths from all causes including 1235 cardiovascular disease (CVD) deaths and 868 cancer deaths. Ancestry‐stratified meta‐analysis of all‐cause mortality identified 163 CpGs in European ancestry (EA) and 17 in African ancestry (AA) participants at p < 1 × 10⁻⁷, of which 41 (EA) and 16 (AA) were also associated with CVD death, and 15 (EA) and 9 (AA) with cancer death. We built DNAm‐based prediction models for all‐cause mortality that predicted mortality risk after adjusting for clinical risk factors. The mortality prediction model trained by integrating DNAm with clinical risk factors showed an improvement in prediction of cancer death with 5% increase in the C‐index in a replication cohort, compared with the model including clinical risk factors alone. Mendelian randomization identified 15 putatively causal CpGs in relation to longevity, CVD, or cancer risk. For example, cg06885782 (in KCNQ4) was positively associated with risk for prostate cancer (Beta = 1.2, P _MR = 4.1 × 10⁻⁴) and negatively associated with longevity (Beta = −1.9, P _MR = 0.02). Pathway analysis revealed that genes associated with mortality‐related CpGs are enriched for immune‐ and cancer‐related pathways. We identified replicable DNAm signatures of mortality and demonstrated the potential utility of CpGs as informative biomarkers for prediction of mortality risk.     

Inherited genetic variants associated with glucocorticoid sensitivity in leukaemia cells

Identification of genetic variants associated with glucocorticoids (GC) sensitivity of leukaemia cells may provide insight into potential drug targets and tailored therapy. In the present study, within 72 leukaemic cell lines derived from Japanese patients with B‐cell precursor acute lymphoblastic leukaemia (ALL), we conducted genome‐wide genotyping of single nucleotide polymorphisms (SNP) and attempted to identify genetic variants associated with GC sensitivity and NR3C1 (GC receptor) gene expression. IC50 measures for prednisolone (Pred) and dexamethasone (Dex) were available using an alamarBlue cell viability assay. IC50 values of Pred showed the strongest association with rs904419 (P = 4.34 × 10⁻⁸), located between the FRMD4B and MITF genes. The median IC50 values of prednisolone for cell lines with rs904419 AA (n = 13), AG (n = 31) and GG (n = 28) genotypes were 0.089, 0.139 and 297 µmol/L, respectively. For dexamethasone sensitivity, suggestive association was observed for SNP rs2306888 (P = 1.43 × 10⁻⁶), a synonymous SNP of the TGFBR3 gene. For NR3C1 gene expression, suggestive association was observed for SNP rs11982167 (P = 6.44 × 10⁻⁸), located in the PLEKHA8 gene. These genetic variants may affect GC sensitivity of ALL cells and may give rise to opportunities in personalized medicine for effective and safe chemotherapy in ALL patients.     

Genetic trade‐offs between complex diseases and longevity

Longevity was influenced by many complex diseases and traits. However, the relationships between human longevity and genetic risks of complex diseases were not broadly studied. Here, we constructed polygenic risk scores (PRSs) for 225 complex diseases/traits and evaluated their relationships with human longevity in a cohort with 2178 centenarians and 2299 middle‐aged individuals. Lower genetic risks of stroke and hypotension were observed in centenarians, while higher genetic risks of schizophrenia (SCZ) and type 2 diabetes (T2D) were detected in long‐lived individuals. We further stratified PRSs into cell‐type groups and significance‐level groups. The results showed that the immune component of SCZ genetic risk was positively linked to longevity, and the renal component of T2D genetic risk was the most deleterious. Additionally, SNPs with very small p‐values (p ≤ 1x10^‐5) for SCZ and T2D were negatively correlated with longevity. While for the less significant SNPs (1x10^‐5 < p ≤ 0.05), their effects on disease and longevity were positively correlated. Overall, we identified genetically informed positive and negative factors for human longevity, gained more insights on the accumulation of disease risk alleles during evolution, and provided evidence for the theory of genetic trade‐offs between complex diseases and longevity.     

Identification of Multiple Genetic Susceptibility Loci in Takayasu Arteritis

Takayasu arteritis is a rare inflammatory disease of large arteries. The etiology of Takayasu arteritis remains poorly understood, but genetic contribution to the disease pathogenesis is supported by the genetic association with HLA-B^∗52. We genotyped ∼200,000 genetic variants in two ethnically divergent Takayasu arteritis cohorts from Turkey and North America by using a custom-designed genotyping platform (Immunochip). Additional genetic variants and the classical HLA alleles were imputed and analyzed. We identified and confirmed two independent susceptibility loci within the HLA region (r² < 0.2): HLA-B/MICA (rs12524487, OR = 3.29, p = 5.57 × 10⁻¹⁶) and HLA-DQB1/HLA-DRB1 (rs113452171, OR = 2.34, p = 3.74 × 10⁻⁹; and rs189754752, OR = 2.47, p = 4.22 × 10⁻⁹). In addition, we identified and confirmed a genetic association between Takayasu arteritis and the FCGR2A/FCGR3A locus on chromosome 1 (rs10919543, OR = 1.81, p = 5.89 × 10⁻¹²). The risk allele in this locus results in increased mRNA expression of FCGR2A. We also established the genetic association between IL12B and Takayasu arteritis (rs56167332, OR = 1.54, p = 2.18 × 10⁻⁸).     

Genome-wide association identifies a common variant in the reelin gene that increases the risk of schizophrenia only in women

Sex differences in schizophrenia are well known, but their genetic basis has not been identified. We performed a genome-wide association scan for schizophrenia in an Ashkenazi Jewish population using DNA pooling. We found a female-specific association with rs7341475, a SNP in the fourth intron of the reelin (RELN) gene (p = 2.9 × 10⁻⁵ in women), with a significant gene-sex effect (p = 1.8 × 10⁻⁴). We studied rs7341475 in four additional populations, totaling 2,274 cases and 4,401 controls. A significant effect was observed only in women, replicating the initial result (p = 2.1 × 10⁻³ in women; p = 4.2 × 10⁻³ for gene-sex interaction). Based on all populations the estimated relative risk of women carrying the common genotype is 1.58 (p = 8.8 × 10⁻⁷; p = 1.6 × 10⁻⁵ for gene-sex interaction). The female-specific association between RELN and schizophrenia is one of the few examples of a replicated sex-specific genetic association in any disease.     

Association of Genetic Variants Influencing Lipid Levels with Coronary Artery Disease in Japanese Individuals

Background/Objective

In Japanese populations, we performed a replication study of genetic loci previously identified in European-descent populations as being associated with lipid levels and risk of coronary artery disease (CAD).

Methods

We genotyped 48 single nucleotide polymorphisms (SNPs) from 22 candidate loci that had previously been identified by genome-wide association (GWA) meta-analyses for low-density lipoprotein cholesterol (LDL-C), high-density lipoprotein cholesterol (HDL-C), and/or triglycerides in Europeans. We selected 22 loci with 2 parallel tracks from 95 reported loci: 16 significant loci (p<1×10⁻³⁰ in Europeans) and 6 other loci including those with suggestive evidence of lipid associations in 1292 GWA-scanned Japanese samples. Genotyping was done in 4990 general population samples, and 1347 CAD cases and 1337 controls. For 9 SNPs, we further examined CAD associations in an additional panel of 3052 CAD cases and 6335 controls.

Principal Findings

Significant lipid associations (one-tailed p<0.05) were replicated for 18 of 22 loci in Japanese samples, with significant inter-ethnic heterogeneity at 4 loci–APOB, APOE-C1, CETP, and APOA5–and allelic heterogeneity. The strongest association was detected at APOE rs7412 for LDL-C (p = 1.3×10⁻⁴¹), CETP rs3764261 for HDL-C (p = 5.2×10⁻²⁴), and APOA5 rs662799 for triglycerides (p = 5.8×10⁻⁵⁴). CAD association was replicated and/or verified for 4 loci: SORT1 rs611917 (p = 1.7×10⁻⁸), APOA5 rs662799 (p = 0.0014), LDLR rs1433099 (p = 2.1×10⁻⁷), and APOE rs7412 (p = 6.1×10⁻¹³).

Conclusions

Our results confirm that most of the tested lipid loci are associated with lipid traits in the Japanese, further indicating that in genetic susceptibility to lipid levels and CAD, the related metabolic pathways are largely common across the populations, while causal variants at individual loci can be population-specific.     

Genome-wide Association Study Identifies TNFSF15 and POU2AF1 as Susceptibility Loci for Primary Biliary Cirrhosis in the Japanese Population

For the identification of susceptibility loci for primary biliary cirrhosis (PBC), a genome-wide association study (GWAS) was performed in 963 Japanese individuals (487 PBC cases and 476 healthy controls) and in a subsequent replication study that included 1,402 other Japanese individuals (787 cases and 615 controls). In addition to the most significant susceptibility region, human leukocyte antigen (HLA), we identified two significant susceptibility loci, TNFSF15 (rs4979462) and POU2AF1 (rs4938534) (combined odds ratio [OR] = 1.56, p = 2.84 × 10⁻¹⁴ for rs4979462, and combined OR = 1.39, p = 2.38 × 10⁻⁸ for rs4938534). Among 21 non-HLA susceptibility loci for PBC identified in GWASs of individuals of European descent, three loci (IL7R, IKZF3, and CD80) showed significant associations (combined p = 3.66 × 10⁻⁸, 3.66 × 10⁻⁹, and 3.04 × 10⁻⁹, respectively) and STAT4 and NFKB1 loci showed suggestive association with PBC (combined p = 1.11 × 10⁻⁶ and 1.42 × 10⁻⁷, respectively) in the Japanese population. These observations indicated the existence of ethnic differences in genetic susceptibility loci to PBC and the importance of TNF signaling and B cell differentiation for the development of PBC in individuals of European descent and Japanese individuals.     

Genome-wide Association Study Identifies Five Susceptibility Loci for Follicular Lymphoma outside the HLA Region

Genome-wide association studies (GWASs) of follicular lymphoma (FL) have previously identified human leukocyte antigen (HLA) gene variants. To identify additional FL susceptibility loci, we conducted a large-scale two-stage GWAS in 4,523 case subjects and 13,344 control subjects of European ancestry. Five non-HLA loci were associated with FL risk: 11q23.3 (rs4938573, p = 5.79 × 10⁻²⁰) near CXCR5; 11q24.3 (rs4937362, p = 6.76 × 10⁻¹¹) near ETS1; 3q28 (rs6444305, p = 1.10 × 10⁻¹⁰) in LPP; 18q21.33 (rs17749561, p = 8.28 × 10⁻¹⁰) near BCL2; and 8q24.21 (rs13254990, p = 1.06 × 10⁻⁸) near PVT1. In an analysis of the HLA region, we identified four linked HLA-DRβ1 multiallelic amino acids at positions 11, 13, 28, and 30 that were associated with FL risk (p_omnibus = 4.20 × 10⁻⁶⁷ to 2.67 × 10⁻⁷⁰). Additional independent signals included rs17203612 in HLA class II (odds ratio [OR_per-allele] = 1.44; p = 4.59 × 10⁻¹⁶) and rs3130437 in HLA class I (OR_per-allele = 1.23; p = 8.23 × 10⁻⁹). Our findings further expand the number of loci associated with FL and provide evidence that multiple common variants outside the HLA region make a significant contribution to FL risk.     

Polymorphism in miR-31 and miR-584 binding site in the angiotensinogen gene differentially influences body fat distribution in both sexes

Angiotensinogen (AGT), its active fragments and microRNA-31 (miR-31) play an important role in adipocyte differentiation. AGT contains a miR-31 polymorphic binding site. We hypothesize that the rs7079 polymorphism in the miR-31/584 binding site of the AGT gene could influence body fat distribution. A total of 751 subjects (195 men, 556 women) were enrolled in the study. The rs7079 genotypes were determined by qRT-PCR. Anthropometric measurements were taken on all subjects, who were subsequently divided into two groups: obese (>30 kg m⁻²) and non-obese (<30 kg m⁻²). Linear regression models were created to determine the contributions of sex, obesity status and rs7079 to all measured parameters. Adding the rs7079 genotype significantly contributed to the linear regression model for waist circumference (p = 0.013), hip circumference (p = 0.018) and supraspinal skin-fold thickness (p = 1 × 10⁻³). Differences between sexes and between the obese and non-obese groups were observed. Waist circumference was lower in men carrying the A allele (p = 0.022); hip circumference was higher only in obese women carrying the A allele (p = 0.015). While men carrying the A allele had lower supraspinal skin-fold thickness (p = 0.022), this parameter was found to be higher in A allele carrying women (p = 3 × 10⁻³). The higher total sum of skin-fold thickness in A allele carrying women was restricted to obese individuals (p = 0.028). The presence of the A allele was associated with both lower tricipital skin-fold thickness in non-obese women (p = 0.023) and a trend of higher thickness in non-obese men (p = 0.065). Significant associations of rs7079 in the AGT gene and body fat distribution were observed. The distribution followed opposing patterns in both sexes.     

Modeling of Environmental Effects in Genome-Wide Association Studies Identifies SLC2A2 and HP as Novel Loci Influencing Serum Cholesterol Levels

Genome-wide association studies (GWAS) have identified 38 larger genetic regions affecting classical blood lipid levels without adjusting for important environmental influences. We modeled diet and physical activity in a GWAS in order to identify novel loci affecting total cholesterol, LDL cholesterol, HDL cholesterol, and triglyceride levels. The Swedish (SE) EUROSPAN cohort (N _SE = 656) was screened for candidate genes and the non-Swedish (NS) EUROSPAN cohorts (N _NS = 3,282) were used for replication. In total, 3 SNPs were associated in the Swedish sample and were replicated in the non-Swedish cohorts. While SNP rs1532624 was a replication of the previously published association between CETP and HDL cholesterol, the other two were novel findings. For the latter SNPs, the p-value for association was substantially improved by inclusion of environmental covariates: SNP rs5400 (p _{SE,unadjusted} = 3.6×10⁻⁵, p _SE,adjusted = 2.2×10⁻⁶, p _{NS,unadjusted} = 0.047) in the SLC2A2 (Glucose transporter type 2) and rs2000999 (p _{SE,unadjusted} = 1.1×10⁻³, p _SE,adjusted = 3.8×10⁻⁴, p _{NS,unadjusted} = 0.035) in the HP gene (Haptoglobin-related protein precursor). Both showed evidence of association with total cholesterol. These results demonstrate that inclusion of important environmental factors in the analysis model can reveal new genetic susceptibility loci.     

A Genome-Wide Association Study of Total Serum and Mite-Specific IgEs in Asthma Patients

Immunoglobulin E (IgE) is one of the central players in asthma and allergic diseases. Although the serum IgE level, a useful endophenotype, is generally increased in patients with asthma, genetic factors influencing IgE regulation in asthma are still not fully understood. To identify the genetic variations associated with total serum and mite-specific IgEs in asthmatics, a genome-wide association study (GWAS) of 657,366 single nucleotide polymorphisms (SNPs) was performed in 877 Korean asthmatics. This study found that several new genes might be associated with total IgE in asthmatics, such as CRIM1 (rs848512, P = 1.18×10⁻⁶; rs711254, P = 6.73×10⁻⁶), ZNF71 (rs10404342, P = 7.60×10⁻⁶), TLN1 (rs4879926, P = 7.74×10⁻⁶), and SYNPO2 (rs1472066, P = 8.36×10⁻⁶; rs1038770, P = 8.66×10⁻⁶). Regarding the association of specific IgE to house dust mites, it was observed that intergenic SNPs nearby to OPRK1 and LOC730217 might be associated with Dermatophagoides pteronyssinus (D.p.) and Dermatophagoides farinae (D.f.) in asthmatics, respectively. In further pathway analysis, the phosphatidylinositol signaling system and adherens junction pathways were estimated to play a role in the regulation of total IgE levels in asthma. Although functional evaluations and replications of these results in other populations are needed, this GWAS of serum IgE in asthmatics could facilitate improved understanding of the role of the newly identified genetic variants in asthma and its related phenotypes.     

Genome-Wide Scan Identifies Variant in TNFSF13 Associated with Serum IgM in a Healthy Chinese Male Population

IgM provides a first line of defense during microbial infections. Serum IgM levels are detected routinely in clinical practice. And IgM is a genetically complex trait. We conducted a two-stage genome-wide association study (GWAS) to identify genetic variants affecting serum IgM levels in a Chinese population of 3495, including 1999 unrelated subjects in the first stage and 1496 independent individuals in the second stage. Our data show that a common single nucleotide polymorphism (SNP), rs11552708 located in the TNFSF13 gene was significantly associated with IgM levels (p = 5.00×10⁻⁷ in first stage, p = 1.34×10⁻³ in second stage, and p = 4.22×10⁻⁹ when combined). Besides, smoking was identified to be associated with IgM levels in both stages (P<0.05), but there was no significant interaction between smoking and the identified SNP (P>0.05). It is suggested that TNFSF13 may be a susceptibility gene affecting serum IgM levels in Chinese male population.     

Low base‐substitution mutation rate and predominance of insertion‐deletion events in the acidophilic bacterium Acidobacterium capsulatum

Analyses of spontaneous mutation have shown that total genome‐wide mutation rates are quantitatively similar for most prokaryotic organisms. However, this view is mainly based on organisms that grow best around neutral pH values (6.0–8.0). In particular, the whole‐genome mutation rate has not been determined for an acidophilic organism. Here, we have determined the genome‐wide rate of spontaneous mutation in the acidophilic Acidobacterium capsulatum using a direct and unbiased method: a mutation‐accumulation experiment followed by whole‐genome sequencing. Evaluation of 69 mutation accumulation lines of A. capsulatum after an average of ~2900 cell divisions yielded a base‐substitution mutation rate of 1.22 × 10⁻¹⁰ per site per generation or 4 × 10⁻⁴ per genome per generation, which is significantly lower than the consensus value (2.5−4.6 × 10⁻³) of mesothermophilic (~15–40°C) and neutrophilic (pH 6–8) prokaryotic organisms. However, the insertion‐deletion rate (0.43 × 10⁻¹⁰ per site per generation) is high relative to the base‐substitution mutation rate. Organisms with a similar effective population size and a similar expected effect of genetic drift should have similar mutation rates. Because selection operates on the total mutation rate, it is suggested that the relatively high insertion‐deletion rate may be balanced by a low base‐substitution rate in A. capsulatum, with selection operating on the total mutation rate.     

Clonal hematopoiesis associated with epigenetic aging and clinical outcomes

Clonal hematopoiesis of indeterminate potential (CHIP) is a common precursor state for blood cancers that most frequently occurs due to mutations in the DNA‐methylation modifying enzymes DNMT3A or TET2. We used DNA‐methylation array and whole‐genome sequencing data from four cohorts together comprising 5522 persons to study the association between CHIP, epigenetic clocks, and health outcomes. CHIP was strongly associated with epigenetic age acceleration, defined as the residual after regressing epigenetic clock age on chronological age, in several clocks, ranging from 1.31 years (GrimAge, p < 8.6 × 10⁻⁷) to 3.08 years (EEAA, p < 3.7 × 10⁻¹⁸). Mutations in most CHIP genes except DNA‐damage response genes were associated with increases in several measures of age acceleration. CHIP carriers with mutations in multiple genes had the largest increases in age acceleration and decrease in estimated telomere length. Finally, we found that ~40% of CHIP carriers had acceleration >0 in both Hannum and GrimAge (referred to as AgeAccelHG+). This group was at high risk of all‐cause mortality (hazard ratio 2.90, p < 4.1 × 10⁻⁸) and coronary heart disease (CHD) (hazard ratio 3.24, p < 9.3 × 10⁻⁶) compared to those who were CHIP−/AgeAccelHG−. In contrast, the other ~60% of CHIP carriers who were AgeAccelHG− were not at increased risk of these outcomes. In summary, CHIP is strongly linked to age acceleration in multiple clocks, and the combination of CHIP and epigenetic aging may be used to identify a population at high risk for adverse outcomes and who may be a target for clinical interventions.     

Telomere length and metabolic syndrome traits: A Mendelian randomisation study

Observational studies have revealed associations between short leucocyte telomere length (LTL), a TL marker in somatic tissues and multiple Metabolic Syndrome (MetS) traits. Animal studies have supported these findings by showing that increased telomere attrition leads to adipose tissue dysfunction and insulin resistance. We investigated the associations between genetically instrumented LTL and MetS traits using Mendelian Randomisation (MR). Fifty‐two independent variants identified at FDR<0.05 from a genome‐wide association study (GWAS) including 78,592 Europeans and collectively accounting for 2.93% of LTL variance were selected as genetic instruments for LTL. Summary‐level data for MetS traits and for the MetS as a binary phenotype were obtained from the largest publicly available GWAS and two‐sample MR analyses were used to estimate the associations of LTL with these traits. The combined effect of the genetic instruments was modelled using inverse variance weighted regression and sensitivity analyses with MR‐Egger, weighted‐median and MR‐PRESSO were performed to test for and correct horizonal pleiotropy. Genetically instrumented longer LTL was associated with higher waist‐to‐hip ratio adjusted for body mass index (β = 0.045 SD, SE = 0.018, p = 0.01), raised systolic (β = 1.529 mmHg, SE = 0.332, p = 4x10⁻⁶) and diastolic (β = 0.633 mmHg, SE = 0.222, p = 0.004) blood pressure, and increased MetS risk (OR = 1.133, 95% CI 1.057–1.215). Consistent results were obtained in sensitivity analyses, which provided no evidence of unbalanced horizontal pleiotropy. Telomere shortening might not be a major driver of cellular senescence and dysfunction in human adipose tissue. Future experimental studies should examine the mechanistic bases for the links between longer LTL and increased upper‐body fat distribution and raised blood pressure.