Genetics and beyond--the transcriptome of human monocytes and disease susceptibility

Study on Allele Specific Expression of Long-Term Residents in High Altitude Areas

In diploid organisms, half of the chromosomes in each cell come from the father and half from the mother. Through previous studies, it was found that the paternal chromosome and the maternal chromosome can be regulated and expressed independently, leading to the emergence of allele specific expression (ASE). In this study, we analyzed the differential expression of alleles in the high-altitude population and the normal population based on the RNA sequencing data. Through gene cluster analysis and protein interaction network analysis, we found some changes occurred at the gene level, and some negative effects. During the study, we realized that the calmodulin homology domain may have a certain correlation with long-term survival at high altitude. The plateau environment is characterized by hypoxia, low air pressure, strong ultraviolet radiation, and low temperature. Accordingly, the genetic changes in the process of adaptation are mainly reflected in these characteristics. High altitude generation living is also highly related to cancer, immune disease, cardiovascular disease, neurological disease, endocrine disease, and other diseases. Therefore, the medical system in high altitude areas should pay more attention to these diseases.

Roles and mechanisms of optineurin in bone metabolism

Optineurin (OPTN) is a widely expressed multifunctional articulatory protein that participates in cellular or mitochondrial autophagy, vesicular transport, and endoplasmic reticulum (ER) stress via interactions with various proteins. Skeletal development is a complex biological process that requires the participation of various osteoblasts, such as bone marrow mesenchymal stem cells (BMSCs), and osteogenic, osteoclastic, and chondrogenic cells. OPTN was recently found to be involved in the regulation of osteoblast activity, which affects bone metabolism. OPTN inhibits osteoclastogenesis via signaling pathways, including NF-κB, IFN-β, and NRF2. OPTN can promote the differentiation of BMSCs toward osteogenesis and inhibit lipogenic differentiation by delaying BMSC senescence and autophagy. These effects are closely related to the development of bone metabolism disorders, such as Paget's disease of bone, rheumatoid arthritis, and osteoporosis. Therefore, this review aims to explore the role and mechanism of OPTN in the regulation of bone metabolism and related bone metabolic diseases. Our findings will provide new targets and strategies for the prevention and treatment of bone metabolic diseases.

Association of single nucleotide polymorphisms in the F2RL1 gene with clinical and inflammatory characteristics of patients with asthma

BACKGROUND

Proteinase-activated receptor 2 (PAR-2) is a G-protein coupled receptor associated with many inflammatory diseases, including asthma. We have shown an association between PAR-2 expression in peripheral blood monocytes and asthma severity as well as blood PAR-2 mRNA level and lung function. Since F2RL1 (the gene encoding PAR-2) polymorphisms affect PAR-2 expression, we hypothesize they may affect asthma severity.

METHODS

We recruited 76 subjects with asthma of varying severity and collected clinical (FEV1 [% predicted], FEV1/FVC, IgE) and immunological (PAR-2 mRNA, blood eosinophils) disease parameters. We also genotyped these individuals for 3 F2RL1 SNPs (-45C/T, -149C/G, c.621C/T).

RESULTS

We found that the F2RL1 SNP "C" allele of -45C/T (rs1529505) was associated with PAR-2 mRNA and blood eosinophils. F2RL1 SNP c.621C/T (rs631465) was associated with PAR-2 mRNA. The F2RL1 SNP -149C/G (rs2242991) had no association with any of the parameters studied. This study identified one F2RL1 SNP rs1529505 is associated with parameters of asthma, but not asthma severity.

CONCLUSION

Larger studies are needed to further elucidate the role of PAR-2 in the pathophysiology of asthma and the influence of genetic variation.

Association between genetic variants of transmembrane transporters and susceptibility to anthracycline-induced cardiotoxicity: Current understanding and existing evidence

Anthracyclines remain the cornerstone of numerous chemotherapeutic protocols, with beneficial effects against haematological malignancies and solid tumours. Unfortunately, the clinical usefulness of anthracyclines is compromised by the development of cardiotoxic side effects, leading to dose limitations or treatment discontinuation. There is no absolute linear correlation between the incidence of cardiotoxicity and the threshold dose, suggesting that genetic factors may modify the association between anthracyclines and cardiotoxicity risk. And the majority of single nucleotide polymorphisms (SNPs) associated with anthracycline pharmacogenomics were identified in the ATP-binding cassette (ABC) and solute carrier (SLC) transporters, generating increasing interest in the pharmacogenetic implications of their genetic variations for anthracycline-induced cardiotoxicity (AIC). This review focuses on the influence of SLC and ABC polymorphisms on AIC and highlights the prospects and clinical significance of pharmacogenetics for individualised preventive approaches.

Socioeconomic inequalities in early adulthood disrupt the immune transcriptomic landscape via upstream regulators

Disparities in socio-economic status (SES) predict many immune system-related diseases, and previous research documents relationships between SES and the immune cell transcriptome. Drawing on a bioinformatically-informed network approach, we situate these findings in a broader molecular framework by examining the upstream regulators of SES-associated transcriptional alterations. Data come from the National Longitudinal Study of Adolescent to Adult Health (Add Health), a nationally representative sample of 4543 adults in the United States. Results reveal a network-of differentially expressed genes, transcription factors, and protein neighbors of transcription factors-that shows widespread SES-related dysregulation of the immune system. Mediational models suggest that body mass index (BMI) plays a key role in accounting for many of these associations. Overall, the results reveal the central role of upstream regulators in socioeconomic differences in the molecular basis of immunity, which propagate to increase risk of chronic health conditions in later-life.

Genetic Association and Differential RNA Expression of Histone (De)Acetylation-Related Genes in Pemphigus Foliaceus-A Possible Epigenetic Effect in the Autoimmune Response

Pemphigus foliaceus (PF) is an autoimmune skin blistering disease characterized by antidesmoglein-1 IgG production, with an endemic form (EPF) in Brazil. Genetic and epigenetic factors have been associated with EPF, but its etiology is still not fully understood. To evaluate the genetic association of histone (de)acetylation-related genes with EPF susceptibility, we evaluated 785 polymorphisms from 144 genes, for 227 EPF patients and 194 controls. Carriers of HDAC4_rs4852054*A were more susceptible (OR = 1.79, p = 0.0038), whereas those with GSE1_rs13339618*A (OR = 0.57, p = 0.0011) and homozygotes for PHF21A_rs4756055*A (OR = 0.39, p = 0.0006) were less susceptible to EPF. These variants were not associated with sporadic PF (SPF) in German samples of 75 SPF patients and 150 controls, possibly reflecting differences in SPF and EPF pathophysiology. We further evaluated the expression of histone (de)acetylation-related genes in CD4+ T lymphocytes, using RNAseq. In these cells, we found a higher expression of KAT2B, PHF20, and ZEB2 and lower expression of KAT14 and JAD1 in patients with active EPF without treatment compared to controls from endemic regions. The encoded proteins cause epigenetic modifications related to immune cell differentiation and cell death, possibly affecting the immune response in patients with PF.

Genetics and epidemiology of mutational barcode-defined clonal hematopoiesis

Clonal hematopoiesis (CH) arises when a substantial proportion of mature blood cells is derived from a single hematopoietic stem cell lineage. Using whole-genome sequencing of 45,510 Icelandic and 130,709 UK Biobank participants combined with a mutational barcode method, we identified 16,306 people with CH. Prevalence approaches 50% in elderly participants. Smoking demonstrates a dosage-dependent impact on risk of CH. CH associates with several smoking-related diseases. Contrary to published claims, we find no evidence that CH is associated with cardiovascular disease. We provide evidence that CH is driven by genes that are commonly mutated in myeloid neoplasia and implicate several new driver genes. The presence and nature of a driver mutation alters the risk profile for hematological disorders. Nevertheless, most CH cases have no known driver mutations. A CH genome-wide association study identified 25 loci, including 19 not implicated previously in CH. Splicing, protein and expression quantitative trait loci were identified for CD164 and TCL1A.

Missense variants in SORT1 are associated with LDL-C in an Amish population

Common noncoding variants at the human 1p13.3 locus associated with SORT1 expression are among those most strongly associated with low-density lipoprotein cholesterol (LDL-C) in human genome-wide association studies. However, validation studies in mice and cell lines have produced variable results regarding the directionality of the effect of SORT1 on LDL-C. This, together with the fact that the 1p13.3 variants are associated with expression of several genes, has raised the question of whether SORT1 is the causal gene at this locus. Using whole exome sequencing in members of an Amish population, we identified coding variants in SORT1 that are associated with increased (rs141749679, K302E) and decreased (rs149456022, Q225H) LDL-C. Further, analysis of plasma lipoprotein particle subclasses by ion mobility in a subset of rs141749679 (K302E) carriers revealed higher levels of large LDL particles compared to noncarriers. In contrast to the effect of these variants in the Amish, the sortilin K302E mutation introduced into a C57BL/6J mouse via CRISPR/Cas9 resulted in decreased non-high-density lipoprotein cholesterol, and the sortilin Q225H mutation did not alter cholesterol levels in mice. This is indicative of different effects of these mutations on cholesterol metabolism in the two species. To our knowledge, this is the first evidence that naturally occurring coding variants in SORT1 are associated with LDL-C, thus supporting SORT1 as the gene responsible for the association of the 1p13.3 locus with LDL-C.

A Single nucleotide polymorphism in the ALDH2 gene modifies the risk of esophageal squamous cell carcinoma in BRCA2 p.K3326* carriers

Esophageal squamous cell carcinoma (ESCC) has a very high incidence rate in northeastern Iran. Our team previously reported the BReast CAncer gene 2 (BRCA2) p.K3326* mutation as a moderately penetrant ESCC susceptibility variant in northern Iran (odds ratio (OR) = 3.64, 95% confidence interval (CI) = 1.74-7.59, P = 0.0003). Recently, it has been reported that aldehydes can induce BRCA2 haploinsufficiency in cells with a heterozygous pathogenic BRCA2 mutation and predispose them to carcinogenic effects. Based on this observation, we speculate that dysfunctional variants in Aldehyde Dehydrogenase 2 Family Member (ALDH2) may result in aldehyde-induced BRCA2 haploinsufficiency and increase cancer risk in BRCA2 mutation carriers. In support of this hypothesis, our team recently reported the breast cancer risk modifying effect of an ALDH2 common polymorphism, rs10744777, among Polish carriers of the BRCA2 p.K3326* mutation. In the current case-control study, we aimed to investigate the ESCC risk modifying effect of this ALDH2 polymorphism among BRCA2 p.K3326* mutation carriers. We assessed the interaction between the ALDH2 rs10744777 polymorphism and BRCA2 p.K3326* mutation in ESCC risk by genotyping this ALDH2 variant in the germline DNA of 746 ESCC cases and 1,373 controls from northern Iran who were previously genotyped for the BRCA2 p.K3326* mutation. Among a total of 464 individuals with TT genotype of the ALDH2 rs10744777 polymorphism, which is associated with lower ALDH2 expression, we found 9 of 164 cases versus 3 of 300 controls who carried the BRCA2 p.K3326* variant (OR = 5.66, 95% CI = 1.22-26.2, P = 0.018). This finding supports our hypothesis that the ALDH2-rs10744777 TT genotype may be a significant risk modifier of ESCC in individuals with a BRCA2 p.K3326* mutation.

TNS1 and NRXN1 Genes Interacting With Early-Life Smoking Exposure in Asthma-Plus-Eczema Susceptibility

PURPOSE

Numerous genes have been associated with allergic diseases (asthma, allergic rhinitis, and eczema), but they explain only part of their heritability. This is partly because most previous studies ignored complex mechanisms such as gene-environment (G-E) interactions and complex phenotypes such as co-morbidity. However, it was recently evidenced that the co-morbidity of asthma-plus-eczema appears as a sub-entity depending on specific genetic factors. Besides, evidence also suggest that gene-by-early life environmental tobacco smoke (ETS) exposure interactions play a role in asthma, but were never investigated for asthma-plus-eczema. To identify genetic variants interacting with ETS exposure that influence asthma-plus-eczema susceptibility.

METHODS

To conduct a genome-wide interaction study (GWIS) of asthma-plus-eczema according to ETS exposure, we applied a 2-stage strategy with a first selection of single nucleotide polymorphisms (SNPs) from genome-wide association meta-analysis to be tested at a second stage by interaction meta-analysis. All meta-analyses were conducted across 4 studies including a total of 5,516 European-ancestry individuals, of whom 1,164 had both asthma and eczema.

RESULTS

Two SNPs showed significant interactions with ETS exposure. They were located in 2 genes, NRXN1 (2p16) and TNS1 (2q35), never reported associated and/or interacting with ETS exposure for asthma, eczema or more generally for allergic diseases. TNS1 is a promising candidate gene because of its link to lung and skin diseases with possible interactive effect with tobacco smoke exposure.

CONCLUSIONS

This first GWIS of asthma-plus-eczema with ETS exposure underlines the importance of studying sub-phenotypes such as co-morbidities as well as G-E interactions to detect new susceptibility genes.

Socioeconomic inequalities in young adulthood disrupt the immune transcriptomic landscape via upstream regulators

Disparities in socio-economic status (SES) predict many immune system-related diseases, and previous research documents relationships between SES and the immune cell transcriptome. Drawing on a bioinformatically-informed network approach, we situate these findings in a broader molecular framework by examining the upstream regulators of SES-associated transcriptional alterations. Data come from the National Longitudinal Study of Adolescent to Adult Health (Add Health), a nationally representative sample of 4,543 adults in the United States. Results reveal a network-of differentially-expressed genes, transcription factors, and protein neighbors of transcription factors- that shows widespread SES-related dysregulation of the immune system. Mediational models suggest that body mass index plays a key role in accounting for many of these associations. Overall, the results reveal the central role of upstream regulators in socioeconomic differences in the molecular basis of immunity, which propagate to increase risk of chronic health conditions in later-life.

Sequence variant affects GCSAML splicing, mast cell specific proteins, and risk of urticaria

Urticaria is a skin disorder characterized by outbreaks of raised pruritic wheals. In order to identify sequence variants associated with urticaria, we performed a meta-analysis of genome-wide association studies for urticaria with a total of 40,694 cases and 1,230,001 controls from Iceland, the UK, Finland, and Japan. We also performed transcriptome- and proteome-wide analyses in Iceland and the UK. We found nine sequence variants at nine loci associating with urticaria. The variants are at genes participating in type 2 immune responses and/or mast cell biology (CBLB, FCER1A, GCSAML, STAT6, TPSD1, ZFPM1), the innate immunity (C4), and NF-κB signaling. The most significant association was observed for the splice-donor variant rs56043070[A] (hg38: chr1:247556467) in GCSAML (MAF = 6.6%, OR = 1.24 (95%CI: 1.20-1.28), P-value = 3.6 × 10^-44). We assessed the effects of the variants on transcripts, and levels of proteins relevant to urticaria pathophysiology. Our results emphasize the role of type 2 immune response and mast cell activation in the pathogenesis of urticaria. Our findings may point to an IgE-independent urticaria pathway that could help address unmet clinical need.

Polymorphisms within Autophagy-Related Genes as Susceptibility Biomarkers for Multiple Myeloma: A Meta-Analysis of Three Large Cohorts and Functional Characterization

Multiple myeloma (MM) arises following malignant proliferation of plasma cells in the bone marrow, that secrete high amounts of specific monoclonal immunoglobulins or light chains, resulting in the massive production of unfolded or misfolded proteins. Autophagy can have a dual role in tumorigenesis, by eliminating these abnormal proteins to avoid cancer development, but also ensuring MM cell survival and promoting resistance to treatments. To date no studies have determined the impact of genetic variation in autophagy-related genes on MM risk. We performed meta-analysis of germline genetic data on 234 autophagy-related genes from three independent study populations including 13,387 subjects of European ancestry (6863 MM patients and 6524 controls) and examined correlations of statistically significant single nucleotide polymorphisms (SNPs; p < 1 × 10^-9) with immune responses in whole blood, peripheral blood mononuclear cells (PBMCs), and monocyte-derived macrophages (MDM) from a large population of healthy donors from the Human Functional Genomic Project (HFGP). We identified SNPs in six loci, CD46, IKBKE, PARK2, ULK4, ATG5, and CDKN2A associated with MM risk (p = 4.47 × 10^-4-5.79 × 10^-14). Mechanistically, we found that the ULK4_rs6599175 SNP correlated with circulating concentrations of vitamin D3 (p = 4.0 × 10^-4), whereas the IKBKE_rs17433804 SNP correlated with the number of transitional CD24⁺CD38⁺ B cells (p = 4.8 × 10^-4) and circulating serum concentrations of Monocyte Chemoattractant Protein (MCP)-2 (p = 3.6 × 10^-4). We also found that the CD46_rs1142469 SNP correlated with numbers of CD19⁺ B cells, CD19⁺CD3^- B cells, CD5⁺IgD^- cells, IgM^- cells, IgD^-IgM^- cells, and CD4^-CD8^- PBMCs (p = 4.9 × 10^-4-8.6 × 10^-4) and circulating concentrations of interleukin (IL)-20 (p = 0.00082). Finally, we observed that the CDKN2A_rs2811710 SNP correlated with levels of CD4⁺EMCD45RO⁺CD27^- cells (p = 9.3 × 10^-4). These results suggest that genetic variants within these six loci influence MM risk through the modulation of specific subsets of immune cells, as well as vitamin D3^-, MCP-2^-, and IL20-dependent pathways.

Gene expression variability in long-term survivors of childhood cancer and cancer-free controls in response to ionizing irradiation

BACKGROUND

Differential expression analysis is usually adjusted for variation. However, most studies that examined the expression variability (EV) have used computations affected by low expression levels and did not examine healthy tissue. This study aims to calculate and characterize an unbiased EV in primary fibroblasts of childhood cancer survivors and cancer-free controls (N0) in response to ionizing radiation.

METHODS

Human skin fibroblasts of 52 donors with a first primary neoplasm in childhood (N1), 52 donors with at least one second primary neoplasm (N2 +), as well as 52 N0 were obtained from the KiKme case-control study and exposed to a high (2 Gray) and a low dose (0.05 Gray) of X-rays and sham- irradiation (0 Gray). Genes were then classified as hypo-, non-, or hyper-variable per donor group and radiation treatment, and then examined for over-represented functional signatures.

RESULTS

We found 22 genes with considerable EV differences between donor groups, of which 11 genes were associated with response to ionizing radiation, stress, and DNA repair. The largest number of genes exclusive to one donor group and variability classification combination were all detected in N0: hypo-variable genes after 0 Gray (n = 49), 0.05 Gray (n = 41), and 2 Gray (n = 38), as well as hyper-variable genes after any dose (n = 43). While after 2 Gray positive regulation of cell cycle was hypo-variable in N0, (regulation of) fibroblast proliferation was over-represented in hyper-variable genes of N1 and N2+. In N2+, 30 genes were uniquely classified as hyper-variable after the low dose and were associated with the ERK1/ERK2 cascade. For N1, no exclusive gene sets with functions related to the radiation response were detected in our data.

CONCLUSION

N2+ showed high degrees of variability in pathways for the cell fate decision after genotoxic insults that may lead to the transfer and multiplication of DNA-damage via proliferation, where apoptosis and removal of the damaged genome would have been appropriate. Such a deficiency could potentially lead to a higher vulnerability towards side effects of exposure to high doses of ionizing radiation, but following low-dose applications employed in diagnostics, as well.

Identification of MKNK1 and TOP3A as ovarian endometriosis risk-associated genes using integrative genomic analyses and functional experiments

The risk of endometriosis (EM), which is a common complex gynaecological disease, is related to genetic predisposition. However, it is unclear how genetic variants confer the risk of EM. Here, via Sherlock integrative analysis, we combined large-scale genome-wide association studies (GWAS) summary statistics on EM (N = 245,494) with a blood-based eQTL dataset (N = 1490) to identify EM risk-related genes. For validation, we leveraged two independent eQTL datasets (N = 769) for integration with the GWAS data. Thus, we prioritised 14 genes, including GIMAP4, TOP3A, and NMNAT3, which showed significant association with susceptibility to EM. We also utilised two independent methods, Multi-marker Analysis of GenoMic Annotation and S-PrediXcan, to further validate the EM risk-associated genes. Moreover, protein-protein interaction network analysis showed the 14 genes were functionally connected. Functional enrichment analyses further demonstrated that these genes were significantly enriched in metabolic and immune-related pathways. Differential gene expression analysis showed that in peripheral blood samples from patients with ovarian EM, TOP3A, MKNK1, SIPA1L2, and NUCB1 were significantly upregulated, while HOXB2, GIMAP5, and MGMT were significantly downregulated compared with their expression levels in samples from the controls. Immunohistochemistry further confirmed the increased expression levels of MKNK1 and TOP3A in the ectopic and eutopic endometrium compared to normal endometrium, while HOBX2 was downregulated in the endometrium of women with ovarian EM. Finally, in ex vivo functional experiments, MKNK1 knockdown inhibited ectopic endometrial stromal cells (EESCs) migration and invasion. TOP3A knockdown inhibited EESCs proliferation, migration, and invasion, while promoting their apoptosis. Convergent lines of evidence suggested that MKNK1 and TOP3A are novel EM risk-related genes.

Genetic Regulation of SMC Gene Expression and Splicing Predict Causal CAD Genes

BACKGROUND

Coronary artery disease (CAD) is the leading cause of death worldwide. Recent meta-analyses of genome-wide association studies have identified over 175 loci associated with CAD. The majority of these loci are in noncoding regions and are predicted to regulate gene expression. Given that vascular smooth muscle cells (SMCs) play critical roles in the development and progression of CAD, we aimed to identify the subset of the CAD loci associated with the regulation of transcription in distinct SMC phenotypes.

METHODS

We measured gene expression in SMCs isolated from the ascending aortas of 151 heart transplant donors of various genetic ancestries in quiescent or proliferative conditions and calculated the association of their expression and splicing with ~6.3 million imputed single-nucleotide polymorphism markers across the genome.

RESULTS

We identified 4910 expression and 4412 splicing quantitative trait loci (sQTLs) representing regions of the genome associated with transcript abundance and splicing. A total of 3660 expression quantitative trait loci (eQTLs) had not been observed in the publicly available Genotype-Tissue Expression dataset. Further, 29 and 880 eQTLs were SMC-specific and sex-biased, respectively. We made these results available for public query on a user-friendly website. To identify the effector transcript(s) regulated by CAD loci, we used 4 distinct colocalization approaches. We identified 84 eQTL and 164 sQTL that colocalized with CAD loci, highlighting the importance of genetic regulation of mRNA splicing as a molecular mechanism for CAD genetic risk. Notably, 20% and 35% of the eQTLs were unique to quiescent or proliferative SMCs, respectively. One CAD locus colocalized with a sex-specific eQTL (TERF2IP), and another locus colocalized with SMC-specific eQTL (ALKBH8). The most significantly associated CAD locus, 9p21, was an sQTL for the long noncoding RNA CDKN2B-AS1, also known as ANRIL, in proliferative SMCs.

CONCLUSIONS

Collectively, our results provide evidence for the molecular mechanisms of genetic susceptibility to CAD in distinct SMC phenotypes.

Development of a Dose-Adjusted Polygenic Risk Model for Anthracycline-Induced Cardiotoxicity

BACKGROUND

Anthracyclines, which are effective chemotherapeutic agents, cause cardiac dysfunction in up to 57% of patients. The cumulative anthracycline dose is a crucial predictor of cardiotoxicity; however, the cumulative dose alone cannot explain all cardiotoxic events. Strongly associated genetic variants in SLC28A3, UGT1A6, and RARG contribute to anthracycline-induced cardiotoxicity in pediatric patients and may help identify those most susceptible. This study aimed to examine how these pharmacogenetic effects are modulated by cumulative anthracycline doses in the development of cardiotoxicity.

METHODS

A total of 595 anthracycline-treated children were genotyped and cardiotoxicity cases were identified. A dose-stratified analysis was performed to compare the contributions of SLC28A3 rs7853758, UGT1A6 rs17863783, and RARG rs2229774 variants to the development of cardiotoxicity in low-dose (<150 mg/m2 cumulative dose) and high-dose (>250 mg/m2 cumulative dose) patient groups. Logistic regression was used to model the relationships between the cumulative anthracycline dose, genetic variants, and cardiotoxicity in the full cohort.

RESULTS

At < 150 mg/m2 cumulative anthracycline dose, the SLC28A3 protective variant did not reach statistical significance [odds ratio (OR) 0.46 (95% confidence interval (CI) 0.10-1.45), P = 0.23], but it was statistically significant at doses >250 mg/m2 [OR 0.43 (95% CI 0.22-0.78), P = 0.0093]. Conversely, the UGT1A6 and RARG risk variants were either statistically significant or approaching significance at doses <150 mg/m2 [OR 7.18 (95% CI 1.78-28.4), P = 0.0045 for UGT1A6 and OR 2.76 (95% CI 0.89-7.63), P = 0.057 for RARG], but not at doses >250 mg/m2 [OR 2.91 (95% CI 0.80-11.0), P = 0.10; OR 1.56 (95% CI 0.89-2.75), P = 0.12].

CONCLUSIONS

These findings suggest that the SLC28A3 variant imparts more significant protection for patients receiving higher anthracycline doses, whereas the UGT1A6 and RARG risk variants significantly increased the risk of cardiotoxicity at low anthracycline doses.

G Protein-Coupled Receptor 15 Expression Is Associated with Myocardial Infarction

Beyond the influence of lifestyle-related risk factors for myocardial infarction (MI), the mechanisms of genetic predispositions for MI remain unclear. We sought to identify and characterize differentially expressed genes in early-onset MI in a translational approach. In an observational case−control study, transcriptomes from 112 early-onset MI individuals showed upregulated G protein-coupled receptor 15 (GPR15) expression in peripheral blood mononuclear cells compared to controls (fold change = 1.4, p = 1.87 × 10−7). GPR15 expression correlated with intima-media thickness (β = 0.8498, p = 0.111), C-reactive protein (β = 0.2238, p = 0.0052), ejection fraction (β = −0.9991, p = 0.0281) and smoking (β = 0.7259, p = 2.79 × 10−10). The relation between smoking and MI was diminished after the inclusion of GPR15 expression as mediator in mediation analysis (from 1.27 (p = 1.9 × 10−5) to 0.46 (p = 0.21)). The DNA methylation of two GPR15 sites was 1%/5% lower in early-onset MI individuals versus controls (p = 2.37 × 10−6/p = 0.0123), with site CpG3.98251219 significantly predicting risk for incident MI (hazard ratio = 0.992, p = 0.0177). The nucleotide polymorphism rs2230344 (C/T) within GPR15 was associated with early-onset MI (odds ratio = 3.61, p = 0.044). Experimental validation showed 6.3-fold increased Gpr15 expression in an ischemic mouse model (p < 0.05) and 4-fold increased Gpr15 expression in cardiomyocytes under ischemic stress (p < 0.001). After the induction of MI, Gpr15gfp/gfp mice showed lower survival (p = 0.042) and deregulated gene expression for response to hypoxia and signaling pathways. Using a translational approach, our data provide evidence that GPR15 is linked to cardiovascular diseases, mediating the adverse effects of smoking.

The FES Gene at the 15q26 Coronary-Artery-Disease Locus Inhibits Atherosclerosis

BACKGROUND

Genome-wide association studies have discovered a link between genetic variants on human chromosome 15q26.1 and increased coronary artery disease (CAD) susceptibility; however, the underlying pathobiological mechanism is unclear. This genetic locus contains the FES (FES proto-oncogene, tyrosine kinase) gene encoding a cytoplasmic protein-tyrosine kinase involved in the regulation of cell behavior. We investigated the effect of the 15q26.1 variants on FES expression and whether FES plays a role in atherosclerosis.

METHODS AND RESULTS

Analyses of isogenic monocytic cell lines generated by CRISPR (clustered regularly interspaced short palindromic repeats)-mediated genome editing showed that monocytes with an engineered 15q26.1 CAD risk genotype had reduced FES expression. Small-interfering-RNA-mediated knockdown of FES promoted migration of monocytes and vascular smooth muscle cells. A phosphoproteomics analysis showed that FES knockdown altered phosphorylation of a number of proteins known to regulate cell migration. Single-cell RNA-sequencing revealed that in human atherosclerotic plaques, cells that expressed FES were predominately monocytes/macrophages, although several other cell types including smooth muscle cells also expressed FES. There was an association between the 15q26.1 CAD risk genotype and greater numbers of monocytes/macrophage in human atherosclerotic plaques. An animal model study demonstrated that Fes knockout increased atherosclerotic plaque size and within-plaque content of monocytes/macrophages and smooth muscle cells, in apolipoprotein E-deficient mice fed a high fat diet.

CONCLUSIONS

We provide substantial evidence that the CAD risk variants at the 15q26.1 locus reduce FES expression in monocytes and that FES depletion results in larger atherosclerotic plaques with more monocytes/macrophages and smooth muscle cells. This study is the first demonstration that FES plays a protective role against atherosclerosis and suggests that enhancing FES activity could be a potentially novel therapeutic approach for CAD intervention.

Socioeconomic inequalities in molecular risk for chronic diseases observed in young adulthood

Many common chronic diseases of aging are negatively associated with socioeconomic status (SES). This study examines whether inequalities can already be observed in the molecular underpinnings of such diseases in the 30s, before many of them become prevalent. Data come from the National Longitudinal Study of Adolescent to Adult Health (Add Health), a large, nationally representative sample of US subjects who were followed for over two decades beginning in adolescence. We now have transcriptomic data (mRNA-seq) from a random subset of 4,543 of these young adults. SES in the household-of-origin and in young adulthood were examined as covariates of a priori-defined mRNA-based disease signatures and of specific gene transcripts identified de novo. An SES composite from young adulthood predicted many disease signatures, as did income and subjective status. Analyses highlighted SES-based inequalities in immune, inflammatory, ribosomal, and metabolic pathways, several of which play central roles in senescence. Many genes are also involved in transcription, translation, and diverse signaling mechanisms. Average causal-mediated effect models suggest that body mass index plays a key role in accounting for these relationships. Overall, the results reveal inequalities in molecular risk factors for chronic diseases often decades before diagnoses and suggest future directions for social signal transduction models that trace how social circumstances regulate the human genome.

Sorting through the extensive and confusing roles of sortilin in metabolic disease

Sortilin is a post-Golgi trafficking receptor homologous to the yeast vacuolar protein sorting receptor 10 (VPS10). The VPS10 motif on sortilin is a 10-bladed β-propeller structure capable of binding more than 50 proteins, covering a wide range of biological functions including lipid and lipoprotein metabolism, neuronal growth and death, inflammation, and lysosomal degradation. Sortilin has a complex cellular trafficking itinerary, where it functions as a receptor in the trans-Golgi network, endosomes, secretory vesicles, multivesicular bodies, and at the cell surface. In addition, sortilin is associated with hypercholesterolemia, Alzheimer's disease, prion diseases, Parkinson's disease, and inflammation syndromes. The 1p13.3 locus containing SORT1, the gene encoding sortilin, carries the strongest association with LDL-C of all loci in human genome-wide association studies. However, the mechanism by which sortilin influences LDL-C is unclear. Here, we review the role sortilin plays in cardiovascular and metabolic diseases and describe in detail the large and often contradictory literature on the role of sortilin in the regulation of LDL-C levels.

Arsenic Exposure, Blood DNA Methylation, and Cardiovascular Disease

BACKGROUND

Epigenetic dysregulation has been proposed as a key mechanism for arsenic-related cardiovascular disease (CVD). We evaluated differentially methylated positions (DMPs) as potential mediators on the association between arsenic and CVD.

METHODS

Blood DNA methylation was measured in 2321 participants (mean age 56.2, 58.6% women) of the Strong Heart Study, a prospective cohort of American Indians. Urinary arsenic species were measured using high-performance liquid chromatography coupled to inductively coupled plasma mass spectrometry. We identified DMPs that are potential mediators between arsenic and CVD. In a cross-species analysis, we compared those DMPs with differential liver DNA methylation following early-life arsenic exposure in the apoE knockout (apoE-/-) mouse model of atherosclerosis.

RESULTS

A total of 20 and 13 DMPs were potential mediators for CVD incidence and mortality, respectively, several of them annotated to genes related to diabetes. Eleven of these DMPs were similarly associated with incident CVD in 3 diverse prospective cohorts (Framingham Heart Study, Women's Health Initiative, and Multi-Ethnic Study of Atherosclerosis). In the mouse model, differentially methylated regions in 20 of those genes and DMPs in 10 genes were associated with arsenic.

CONCLUSIONS

Differential DNA methylation might be part of the biological link between arsenic and CVD. The gene functions suggest that diabetes might represent a relevant mechanism for arsenic-related cardiovascular risk in populations with a high burden of diabetes.

Patterns and distribution of de novo mutations in multiplex Middle Eastern families

While de novo mutations (DNMs) are key to genetic diversity, they are also responsible for a high number of rare disorders. To date, no study has systematically examined the rate and distribution of DNMs in multiplex families in highly consanguineous populations. Leveraging WGS profiles of 645 individuals in 146 families, we implemented a combinatorial approach using 3 complementary tools for DNM discovery in 353 unique trio combinations. We found a total of 27,168 DNMs (median: 70 single-nucleotide and 6 insertion-deletions per individual). Phasing revealed around 80% of DNMs were paternal in origin. Notably, using whole-genome methylation data of spermatogonial stem cells, these DNMs were significantly more likely to occur at highly methylated CpGs (OR: 2.03; p value = 6.62 × 10⁻¹¹). We then examined the effects of consanguinity and ethnicity on DNMs, and found that consanguinity does not seem to correlate with DNM rate, and special attention has to be considered while measuring such a correlation. Additionally, we found that Middle-Eastern families with Arab ancestry had fewer DNMs than African families, although not significant (p value = 0.16). Finally, for families with diseased probands, we examined the difference in DNM counts and putative impact across affected and unaffected siblings, but did not find significant differences between disease groups, likely owing to the enrichment for recessive disorders in this part of the world, or the small sample size per clinical condition. This study serves as a reference for DNM discovery in multiplex families from the globally under-represented populations of the Middle-East.

A de novo start-loss in EFTUD2 associated with mandibulofacial dysostosis with microcephaly: case report

Mandibulofacial dysostosis with microcephaly (MFDM) is a rare genetic disorder inherited in an autosomal dominant pattern. Major characteristics include developmental delay, craniofacial malformations such as malar and mandibular hypoplasia, and ear anomalies. Here, we report a 4.5-yr-old female patient with symptoms fitting MFDM. Using whole-genome sequencing, we identified a de novo start-codon loss (c.3G > T) in the EFTUD2. We examined EFTUD2 expression in the patient by RNA sequencing and observed a notable functional consequence of the variant on gene expression in the patient. We identified a novel variant for the development of MFDM in humans. To the best of our knowledge, this is the first report of a start-codon loss in EFTUD2 associated with MFDM.

Genetic architecture of band neutrophil fraction in Iceland

The characteristic lobulated nuclear morphology of granulocytes is partially determined by composition of nuclear envelope proteins. Abnormal nuclear morphology is primarily observed as an increased number of hypolobulated immature neutrophils, called band cells, during infection or in rare envelopathies like Pelger-Huët anomaly. To search for sequence variants affecting nuclear morphology of granulocytes, we performed a genome-wide association study using band neutrophil fraction from 88,101 Icelanders. We describe 13 sequence variants affecting band neutrophil fraction at nine loci. Five of the variants are at the Lamin B receptor (LBR) locus, encoding an inner nuclear membrane protein. Mutations in LBR are linked to Pelger-Huët anomaly. In addition, we identify cosegregation of a rare stop-gain sequence variant in LBR and Pelger Huët anomaly in an Icelandic eight generation pedigree, initially reported in 1963. Two of the other loci include genes which, like LBR, play a role in the nuclear membrane function and integrity. These GWAS results highlight the role proteins of the inner nuclear membrane have as important for neutrophil nuclear morphology.

A GWAS in Iceland reveals that variants in inner nuclear membrane proteins are associated with nuclear morphology of granulocytes and band neutrophil fraction.

The hematopoietic compartment is sufficient for lupus development resulting from the POLB-Y265C mutation

Systemic lupus erythematosus is a chronic disease characterized by autoantibodies, renal and cutaneous disease, and immune complex formation. Emerging evidence suggests that aberrant DNA repair is an underlying mechanism of lupus development. We previously showed that the POLBY265C/C mutation, which results in development of an aberrant immune repertoire, leads to lupus-like disease in mice. To address whether the hematopoietic compartment is sufficient for lupus development, we transplanted bone marrow cells from POLBY265C/C and POLB+/+ into wild-type congenic mice. Only mice transplanted with the POLBY265C/C bone marrow develop high levels of antinuclear antibodies and renal disease. In conclusion, we show that the hematopoietic compartment harvested from the POLBY265C/C mice is sufficient for development of autoimmune disease.

Systematic Review of Pharmacogenetics of ABC and SLC Transporter Genes in Acute Myeloid Leukemia

Antineoplastic uptake by blast cells in acute myeloid leukemia (AML) could be influenced by influx and efflux transporters, especially solute carriers (SLCs) and ATP-binding cassette family (ABC) pumps. Genetic variability in SLC and ABC could produce interindividual differences in clinical outcomes. A systematic review was performed to evaluate the influence of SLC and ABC polymorphisms and their combinations on efficacy and safety in AML cohorts. Anthracycline intake was especially influenced by SLCO1B1 polymorphisms, associated with lower hepatic uptake, showing higher survival rates and toxicity in AML studies. The variant alleles of ABCB1 were related to anthracycline intracellular accumulation, increasing complete remission, survival and toxicity. Similar findings have been suggested with ABCC1 and ABCG2 polymorphisms. Polymorphisms of SLC29A1, responsible for cytarabine uptake, demonstrated significant associations with survival and response in Asian populations. Promising results were observed with SLC and ABC combinations regarding anthracycline toxicities. Knowledge of the role of transporter pharmacogenetics could explain the differences observed in drug disposition in the blast. Further studies including novel targeted therapies should be performed to determine the influence of genetic variability to individualize chemotherapy schemes.

Common Variant in ALDH2 Modifies the Risk of Breast Cancer Among Carriers of the p.K3326* Variant in BRCA2

PURPOSE

The BRCA2 p.K3326* variant is considered a low-penetrance variant for breast cancer. Aldehydes that accumulate in cells under insufficient aldehyde oxidation were most recently shown to trigger carcinogenesis by promoting depletion of BRCA2 protein. Allele T of the common variant rs10744777 in the ALDH2 gene was associated with reduced expression of aldehyde dehydrogenase, the main enzyme in aldehyde oxidation. We hypothesized that this allele could modify breast cancer risk in women with the BRCA2 p.K3326* low-penetrance variant through reduced function of ALDH2, increased accumulation of cellular aldehydes, and depletion of BRCA2 protein.

MATERIALS AND METHODS

We genotyped 11,873 Polish women diagnosed with breast cancer and 7,615 ethnically matched controls for these two variants. Next, we extended our analysis of rs10744777 to 231 carriers of pathogenic BRCA2 mutations.

RESULTS

BRCA2 p.K3326* variant was associated with significant increase in breast cancer risk only in those who were homozygous for the T allele of the ALDH2 rs10744777 variant (odds ratio = 1.72; 95% CI, 1.19 to 2.48; P = .003). The BRCA2 p.K3326* variant did not increase the risk of breast cancer among those who were heterozygous or homozygous for the C allele of the ALDH2 rs10744777 variant (odds ratio = 1.05; 95% CI, 0.73 to 1.51; P = .81). In the carriers of high-risk BRCA2 mutations, the TT genotype of rs10744777 conferred a modest (18%) and not significant increase in breast cancer risk.

CONCLUSION

Our results suggest that BRCA2 p.K3326* variant, which is low-penetrance by itself, confers increased breast cancer risk on the background of the TT genotype of the ALDH2 rs10744777 variant in the Polish population.

Identification of a shared genetic risk locus for Kawasaki disease and immunoglobulin A vasculitis by a cross-phenotype meta-analysis

OBJECTIVES

Combining of genomic data of different pathologies as a single phenotype has emerged as a useful strategy to identify genetic risk loci shared among immune-mediated diseases. Our study aimed to increase our knowledge of the genetic contribution to Kawasaki disease (KD) and IgA vasculitis (IgAV) by performing the first comprehensive large-scale analysis on the genetic overlap between them.

METHODS

A total of 1190 vasculitis patients and 11 302 healthy controls were analysed. First, in the discovery phase, genome-wide data of 405 KD patients and 6252 controls and 215 IgAV patients and 1324 controls, all of European origin, were combined using an inverse variance meta-analysis. Second, the top associated polymorphisms were selected for replication in additional independent cohorts (570 cases and 3726 controls). Polymorphisms with P-values ≤5 × 10-8 in the global IgAV-KD meta-analysis were considered as shared genetic risk loci.

RESULTS

A genetic variant, rs3743841, located in an intron of the NAGPA gene, reached genome-wide significance in the cross-disease meta-analysis (P = 8.06 × 10-10). Additionally, when IgAV was individually analysed, a strong association between rs3743841 and this vasculitis was also evident [P = 1.25 × 10-7; odds ratio = 1.47 (95% CI 1.27, 1.69)]. In silico functional annotation showed that this polymorphism acts as a regulatory variant modulating the expression levels of the NAGPA and SEC14L5 genes.

CONCLUSION

We identified a new risk locus with pleiotropic effects on the two childhood vasculitides analysed. This locus represents the strongest non-HLA signal described for IgAV to date.

Periodontitis Risk Variants at SIGLEC5 Impair ERG and MAFB Binding

Periodontitis is a common complex inflammatory disease of the oral cavity. It is characterized by inflammation of gingival tissues and alveolar bone loss. Recently, a genome-wide association study and 2 genome-wide association study meta-analyses found 2 associated regions (haplotype blocks) at the inhibitory immune receptor gene SIGLEC5 to increase the risk for periodontitis. The aims of the current study were the identification of the putative causal variants underlying these associations, characterization of their molecular biological effects, and validation of SIGLEC5 as the target gene. We mapped the associated single-nucleotide polymorphisms to DNA elements with predictive features of regulatory functions and screened the associated alleles for transcription factor (TF) binding sites. Antibody electrophoretic mobility shift assays (EMSAs) with allele-specific probes were used to identify TF binding and to quantify allele-specific effects on binding affinities. Luciferase reporter assays were used to quantify the effect directions and allele-specific strength of the associated regulatory elements. We used CRISPR-dCas9 gene activation to validate SIGLEC5 as a target of the association. EMSA in peripheral blood mononuclear cells showed that E-26 transformation-specific TF-related gene (ERG) binds at rs11084095, with almost complete loss of binding at the minor A-allele. Allele-specific reporter genes showed enhancer function of the DNA sequence at rs11084095, which was abrogated in the background of the A-allele. EMSA in B lymphocytes showed that TF MAF bZIP (MAFB) binds at the common G-allele of rs4284742, whereas the minor A-allele reduced TF binding by 69%, corresponding to 9-fold reduction of luciferase reporter gene activity by the A-allele. Using CRISPR-dCas9, we showed that the enhancer at rs4284742 strongly activated SIGLEC5 expression, validating this gene as the target gene of the association. We conclude that rs11084095 and rs4284742 are putatively causal for the genome-wide significant associations with periodontitis at SIGLEC5 that impair ERG and MAFB binding, respectively.

Hygiene Hypothesis as the Etiology of Kawasaki Disease: Dysregulation of Early B Cell Development

Kawasaki disease (KD) is an acute systemic vasculitis that occurs predominantly in children under 5 years of age. Despite much study, the etiology of KD remains unknown. However, epidemiological and immunological data support the hygiene hypothesis as a possible etiology. It is thought that more sterile or clean modern living environments due to increased use of sanitizing agents, antibiotics, and formula feeding result in a lack of immunological challenges, leading to defective or dysregulated B cell development, accompanied by low IgG and high IgE levels. A lack of B cell immunity may increase sensitivity to unknown environmental triggers that are nonpathogenic in healthy individuals. Genetic studies of KD show that all of the KD susceptibility genes identified by genome-wide association studies are involved in B cell development and function, particularly in early B cell development (from the pro-B to pre-B cell stage). The fact that intravenous immunoglobulin is an effective therapy for KD supports this hypothesis. In this review, I discuss clinical, epidemiological, immunological, and genetic studies showing that the etiopathogenesis of KD in infants and toddlers can be explained by the hygiene hypothesis, and particularly by defects or dysregulation during early B cell development.

Seed sequence polymorphism rs2168518 and allele-specific target gene regulation of hsa-miR-4513

MicroRNAs (miRNAs) are small post-transcriptional regulators that offer promising targets for treating complex diseases. To this end, hsa-miR-4513 is an excellent candidate as this gene harbors within its conserved heptametrical seed sequence a frequent polymorphism (rs2168518), which has previously been associated with several complex phenotypes. So far, little is known about the biological mechanism(s) underlying these associations. In an initial step, we now aimed to identify allele-specific target genes of hsa-miR-4513. We performed RNA sequencing in a miRNA overexpression model in human umbilical vein endothelial cells transfected with separated hsa-miR-4513 alleles at rs2168518, namely hsa-miR-4513-G and hsa-miR-4513-A. Genes specifically regulated by the rs2168518 alleles were independently verified by quantitative reverse transcription PCR (qRT-PCR), western blot analysis and allele-specific miRNA binding via a luciferase reporter assay. By a text-based search publicly available databases such as Online Mendelian Inheritance in Man and Mouse Genome Informatics were utilized to link target genes of hsa-miR-4513 to previously described phenotypes. Overall, we identified 23 allele-specific hsa-miR-4513 target genes and replicated 19 of those independently via qRT-PCR. Western blot analysis and luciferase reporter assays conducted for an exemplary subsample further confirmed the allele-specific regulation of these genes by hsa-miR-4513. Remarkably, multiple allele-specific target genes identified are linked via text retrieval to several phenotypes previously reported to be associated with hsa-miR-4513. These genes offer promising candidates for ongoing research on the functional pathobiological impact of hsa-miR-4513 and its seed polymorphism rs2168518. This could give rise to therapeutic applications targeting this miRNA.

Leveraging auxiliary data from arbitrary distributions to boost GWAS discovery with Flexible cFDR

Genome-wide association studies (GWAS) have identified thousands of genetic variants that are associated with complex traits. However, a stringent significance threshold is required to identify robust genetic associations. Leveraging relevant auxiliary covariates has the potential to boost statistical power to exceed the significance threshold. Particularly, abundant pleiotropy and the non-random distribution of SNPs across various functional categories suggests that leveraging GWAS test statistics from related traits and/or functional genomic data may boost GWAS discovery. While type 1 error rate control has become standard in GWAS, control of the false discovery rate can be a more powerful approach. The conditional false discovery rate (cFDR) extends the standard FDR framework by conditioning on auxiliary data to call significant associations, but current implementations are restricted to auxiliary data satisfying specific parametric distributions, typically GWAS p-values for related traits. We relax these distributional assumptions, enabling an extension of the cFDR framework that supports auxiliary covariates from arbitrary continuous distributions ("Flexible cFDR"). Our method can be applied iteratively, thereby supporting multi-dimensional covariate data. Through simulations we show that Flexible cFDR increases sensitivity whilst controlling FDR after one or several iterations. We further demonstrate its practical potential through application to an asthma GWAS, leveraging various functional genomic data to find additional genetic associations for asthma, which we validate in the larger, independent, UK Biobank data resource.

Inflammatory Protein Profiles in Plasma of Candidaemia Patients and the Contribution of Host Genetics to Their Variability

Circulatory inflammatory proteins play a significant role in anti-Candida host immune defence. However, little is known about the genetic variation that contributes to the variability of inflammatory responses in response to C. albicans. To systematically characterize inflammatory responses in Candida infection, we profiled 91 circulatory inflammatory proteins in peripheral blood mononuclear cells (PBMCs) stimulated with C. albicans yeast isolated from 378 individuals of European origin from the 500 Functional Genomics (500FG) cohort of the Human Functional Genomics Project (HFGP) and Lifelines Deep cohort. To identify the genetic factors that determine variation in inflammatory protein responses, we correlated genome-wide single nucleotide polymorphism (SNP) genotypes with protein abundance (protein quantitative trait loci, pQTLs) produced by the Candida-stimulated PBMCs. Furthermore, we investigated whether differences in survival of candidaemia patients can be explained by modulating levels of inflammatory proteins. We identified five genome-wide significant pQTLs that modulate IL-8, MCP-2, MMP-1, and CCL3 in response to C. albicans. In addition, our genetic analysis suggested that GADD45G from rs10114707 locus that reached genome-wide significance could be a potential core gene that regulates a cytokine network upon Candida infection. Last but not least, we observed that a trans-pQTL marked from SNP rs7651677 at chromosome 3 that influences urokinase plasminogen activator (uPA) is strongly associated with patient survival (P_survival = 3.52 x 10^-5, OR 3). Overall, our genetic analysis showed that genetic variation determines the abundance of circulatory proteins in response to Candida infection.

Proteins Associated with Risk of Kidney Function Decline in the General Population

BACKGROUND

Proteomic profiling may allow identification of plasma proteins that associate with subsequent changesin kidney function, elucidating biologic processes underlying the development and progression of CKD.

METHODS

We quantified the association between 4877 plasma proteins and a composite outcome of ESKD or decline in eGFR by ≥50% among 9406 participants in the Atherosclerosis Risk in Communities (ARIC) Study (visit 3; mean age, 60 years) who were followed for a median of 14.4 years. We performed separate analyses for these proteins in a subset of 4378 participants (visit 5), who were followed at a later time point, for a median of 4.4 years. For validation, we evaluated proteins with significant associations (false discovery rate <5%) in both time periods in 3249 participants in the Chronic Renal Insufficiency Cohort (CRIC) and 703 participants in the African American Study of Kidney Disease and Hypertension (AASK). We also compared the genetic determinants of protein levels with those from a meta-analysis genome-wide association study of eGFR.

RESULTS

In models adjusted for multiple covariates, including baseline eGFR and albuminuria, we identified 13 distinct proteins that were significantly associated with the composite end point in both time periods, including TNF receptor superfamily members 1A and 1B, trefoil factor 3, and β-trace protein. Of these proteins, 12 were also significantly associated in CRIC, and nine were significantly associated in AASK. Higher levels of each protein associated with higher risk of 50% eGFR decline or ESKD. We found genetic evidence for a causal role for one protein, lectin mannose-binding 2 protein (LMAN2).

CONCLUSIONS

Large-scale proteomic analysis identified both known and novel proteomic risk factors for eGFR decline.

Transcriptional Regulation of RUNX1: An Informatics Analysis

The RUNX1/AML1 gene encodes a developmental transcription factor that is an important regulator of haematopoiesis in vertebrates. Genetic disruptions to the RUNX1 gene are frequently associated with acute myeloid leukaemia. Gene regulatory elements (REs), such as enhancers located in non-coding DNA, are likely to be important for Runx1 transcription. Non-coding elements that modulate Runx1 expression have been investigated over several decades, but how and when these REs function remains poorly understood. Here we used bioinformatic methods and functional data to characterise the regulatory landscape of vertebrate Runx1. We identified REs that are conserved between human and mouse, many of which produce enhancer RNAs in diverse tissues. Genome-wide association studies detected single nucleotide polymorphisms in REs, some of which correlate with gene expression quantitative trait loci in tissues in which the RE is active. Our analyses also suggest that REs can be variant in haematological malignancies. In summary, our analysis identifies features of the RUNX1 regulatory landscape that are likely to be important for the regulation of this gene in normal and malignant haematopoiesis.

Comprehensive transcriptome analysis of peripheral blood unravels key lncRNAs implicated in ABPA and asthma

Allergic bronchopulmonary aspergillosis (ABPA) is a complex hypersensitivity lung disease caused by a fungus known as Aspergillus fumigatus. It complicates and aggravates asthma. Despite their potential associations, the underlying mechanisms of asthma developing into ABPA remain obscure. Here we performed an integrative transcriptome analysis based on three types of human peripheral blood, which derived from ABPA patients, asthmatic patients and health controls, aiming to identify crucial lncRNAs implicated in ABPA and asthma. Initially, a high-confidence dataset of lncRNAs was identified using a stringent filtering pipeline. A comparative mutational analysis revealed no significant difference among these samples. Differential expression analysis disclosed several immune-related mRNAs and lncRNAs differentially expressed in ABPA and asthma. For each disease, three sub-networks were established using differential network analysis. Many key lncRNAs implicated in ABPA and asthma were identified, respectively, i.e., AL139423.1-201, AC106028.4-201, HNRNPUL1-210, PUF60-218 and SREBF1-208. Our analysis indicated that these lncRNAs exhibits in the loss-of-function networks, and the expression of which were repressed in the occurrences of both diseases, implying their important roles in the immune-related processes in response to the occurrence of both diseases. Above all, our analysis proposed a new point of view to explore the relationship between ABPA and asthma, which might provide new clues to unveil the pathogenic mechanisms for both diseases.

Identification of OCA2 as a novel locus for the co-morbidity of asthma-plus-eczema

BACKGROUND

Numerous genes have been associated with the three most common allergic diseases (asthma, allergic rhinitis or eczema) but these genes explain only a part of the heritability. In the vast majority of genetic studies, complex phenotypes such as co-morbidity of two of these diseases, have not been considered. This may partly explain missing heritability.

OBJECTIVE

To identify genetic variants specifically associated with the co-morbidity of asthma-plus-eczema.

METHODS

We first conducted a meta-analysis of four GWAS (Genome-Wide Association Study) of the combined asthma-plus-eczema phenotype (total of 8807 European-ancestry subjects of whom 1208 subjects had both asthma and eczema). To assess whether the association with SNP(s) was specific to the co-morbidity, we also conducted a meta-analysis of homogeneity test of association according to disease status ("asthma-plus-eczema" vs. the presence of only one disease "asthma only or eczema only"). We then used a joint test by combining the two test statistics from the co-morbidity-SNP association and the phenotypic heterogeneity of SNP effect meta-analyses.

RESULTS

Seven SNPs were detected for specific association to the asthma-plus-eczema co-morbidity, two with significant and five with suggestive evidence using the joint test after correction for multiple testing. The two significant SNPs are located in the OCA2 gene (Oculocutaneous Albinism II), a new locus never detected for significant evidence of association with any allergic disease. This gene is a promising candidate gene, because of its link to skin and lung diseases, and to epithelial barrier and immune mechanisms.

CONCLUSION

Our study underlines the importance of studying sub-phenotypes as co-morbidities to detect new susceptibility genes.

Genetic Determinants of Paget's Disease of Bone

PURPOSE OF REVIEW

To provide an overview of the role of genes and loci that predispose to Paget's disease of bone and related disorders.

RECENT FINDINGS

Studies over the past ten years have seen major advances in knowledge on the role of genetic factors in Paget's disease of bone (PDB). Genome wide association studies have identified six loci that predispose to the disease whereas family based studies have identified a further eight genes that cause PDB. This brings the total number of genes and loci implicated in PDB to fourteen. Emerging evidence has shown that a number of these genes also predispose to multisystem proteinopathy syndromes where PDB is accompanied by neurodegeneration and myopathy due to the accumulation of abnormal protein aggregates, emphasising the importance of defects in autophagy in the pathogenesis of PDB. Genetic factors play a key role in the pathogenesis of PDB and the studies in this area have identified several genes previously not suspected to play a role in bone metabolism. Genetic testing coupled to targeted therapeutic intervention is being explored as a way of halting disease progression and improving outcome before irreversible skeletal damage has occurred.

DNA Damage Response in the Adaptive Arm of the Immune System: Implications for Autoimmunity

In complex environments, cells have developed molecular responses to confront threats against the genome and achieve the maintenance of genomic stability assuring the transfer of undamaged DNA to their progeny. DNA damage response (DDR) mechanisms may be activated upon genotoxic or environmental agents, such as cytotoxic drugs or ultraviolet (UV) light, and during physiological processes requiring DNA transactions, to restore DNA alterations that may cause cellular malfunction and affect viability. In addition to the DDR, multicellular organisms have evolved specialized immune cells to respond and defend against infections. Both adaptive and innate immune cells are subjected to DDR processes, either as a prerequisite to the immune response, or as a result of random endogenous and exogenous insults. Aberrant DDR activities have been extensively studied in the immune cells of the innate arm, but not in adaptive immune cells. Here, we discuss how the aberrant DDR may lead to autoimmunity, with emphasis on the adaptive immune cells and the potential of therapeutic targeting.

Defining functional variants associated with Alzheimer's disease in the induced immune response

Defining the mechanisms involved in the aetiology of Alzheimer’s disease from genome-wide association studies alone is challenging since Alzheimer’s disease is polygenic and most genetic variants are non-coding. Non-coding Alzheimer’s disease risk variants can influence gene expression by affecting miRNA binding and those located within enhancers and within CTCF sites may influence gene expression through alterations in chromatin states. In addition, their function can be cell-type specific. They can function specifically in microglial enhancers thus affecting gene expression in the brain. Hence, transcriptome-wide association studies have been applied to test the genetic association between disease risk and cell-/tissue-specific gene expression. Many Alzheimer’s disease-associated loci are involved in the pathways of the innate immune system. Both microglia, the primary immune cells of the brain, and monocytes which can infiltrate the brain and differentiate into activated macrophages, have roles in neuroinflammation and β‐amyloid clearance through phagocytosis. In monocytes the function of regulatory variants can be context-specific after immune stimulation. To dissect the variants associated with Alzheimer’s disease in the context of monocytes, we utilized data from naïve monocytes and following immune stimulation in vitro, in combination with genome-wide association studies of Alzheimer’s disease in transcriptome-wide association studies. Of the nine genes with statistically independent transcriptome-wide association signals, seven are located in known Alzheimer’s disease risk loci: BIN1, PTK2B, SPI1, MS4A4A, MS4A6E, APOE and PVR. The transcriptome-wide association signal for MS4A6E, PTK2B and PVR and the direction of effect replicated in an independent genome-wide association studies. Our analysis identified two novel candidate genes for Alzheimer’s disease risk, LACTB2 and PLIN2/ADRP. LACTB2 replicated in a transcriptome-wide association study using independent expression weights. LACTB2 and PLIN2/ADRP are involved in mitochondrial function and lipid metabolism, respectively. Comparison of transcriptome-wide association study results from monocytes, whole blood and brain showed that the signal for PTK2B is specific to blood and MS4A6E is specific to LPS stimulated monocytes.

Integration of summary data from GWAS and eQTL studies identified novel risk genes for coronary artery disease

Several genetic loci have been reported to be significantly associated with coronary artery disease (CAD) by multiple genome-wide association studies (GWAS). Nevertheless, the biological and functional effects of these genetic variants on CAD remain largely equivocal. In the current study, we performed an integrative genomics analysis by integrating large-scale GWAS data (N = 459,534) and 2 independent expression quantitative trait loci (eQTL) datasets (N = 1890) to determine whether CAD-associated risk single nucleotide polymorphisms (SNPs) exert regulatory effects on gene expression. By using Sherlock Bayesian, MAGMA gene-based, multidimensional scaling (MDS), functional enrichment, and in silico permutation analyses for independent technical and biological replications, we highlighted 4 susceptible genes (CHCHD1, TUBG1, LY6G6C, and MRPS17) associated with CAD risk. Based on the protein-protein interaction (PPI) network analysis, these 4 genes were found to interact with each other. We detected a remarkably altered co-expression pattern among these 4 genes between CAD patients and controls. In addition, 3 genes of CHCHD1 (P = .0013), TUBG1 (P = .004), and LY6G6C (P = .038) showed significantly different expressions between CAD patients and controls. Together, we provide evidence to support that these identified genes such as CHCHD1 and TUBG1 are indicative factors of CAD.

Variation in genes implicated in B-cell development and antibody production affects susceptibility to pemphigus

Pemphigus foliaceus (PF) is an autoimmune blistering skin disease characterized by the presence of pathogenic autoantibodies against desmoglein 1, a component of intercellular desmosome junctions. PF occurs sporadically across the globe and is endemic in some Brazilian regions. Because PF is a B-cell-mediated disease, we aimed to study the impact of variants within genes encoding molecules involved in the different steps of B-cell development and antibody production on the susceptibility of endemic PF. We analysed 3,336 single nucleotide polymorphisms (SNPs) from 167 candidate genes genotyped with Illumina microarray in a cohort of 227 PF patients and 193 controls. After quality control and exclusion of non-informative and redundant SNPs, 607 variants in 149 genes remained in the logistic regression analysis, in which sex and ancestry were included as covariates. Our results revealed 10 SNPs within or nearby 11 genes that were associated with susceptibility to endemic PF (OR >1.56; p < 0.005): rs6657275*G (TGFB2); rs1818545*A (RAG1/RAG2/IFTAP);rs10781530*A (PAXX), rs10870140*G and rs10781522*A (TRAF2); rs535068*A (TNFRSF1B); rs324011*A (STAT6);rs6432018*C (YWHAQ); rs17149161*C (YWHAG); and rs2070729*C (IRF1). Interestingly, these SNPs have been previously associated with differential gene expression, mostly in peripheral blood, in publicly available databases. For the first time, we show that polymorphisms in genes involved in B-cell development and antibody production confer differential susceptibility to endemic PF, and therefore are candidates for possible functional studies to understand immunoglobulin gene rearrangement and its impact on diseases.

Promising genes and variants to reduce chemotherapy adverse effects in acute lymphoblastic leukemia

Almost two decades ago, the sequencing of the human genome and high throughput technologies came to revolutionize the clinical and therapeutic approaches of patients with complex human diseases. In acute lymphoblastic leukemia (ALL), the most frequent childhood malignancy, these technologies have enabled to characterize the genomic landscape of the disease and have significantly improved the survival rates of ALL patients. Despite this, adverse reactions from treatment such as toxicity, drug resistance and secondary tumors formation are still serious consequences of chemotherapy, and the main obstacles to reduce ALL-related mortality. It is well known that germline variants and somatic mutations in genes involved in drug metabolism impact the efficacy of drugs used in oncohematological diseases therapy. So far, a broader spectrum of clinically actionable alterations that seems to be crucial for the progression and treatment response have been identified. Although these results are promising, it is necessary to put this knowledge into the clinics to help physician make medical decisions and generate an impact in patients' health. This review summarizes the gene variants and clinically actionable mutations that modify the efficacy of antileukemic drugs. Therefore, knowing their genetic status before treatment is critical to reduce severe adverse effects, toxicities and life-threatening consequences in ALL patients.

Host genetics and infectious disease: new tools, insights and translational opportunities

Understanding how human genetics influence infectious disease susceptibility offers the opportunity for new insights into pathogenesis, potential drug targets, risk stratification, response to therapy and vaccination. As new infectious diseases continue to emerge, together with growing levels of antimicrobial resistance and an increasing awareness of substantial differences between populations in genetic associations, the need for such work is expanding. In this Review, we illustrate how our understanding of the host–pathogen relationship is advancing through holistic approaches, describing current strategies to investigate the role of host genetic variation in established and emerging infections, including COVID-19, the need for wider application to diverse global populations mirroring the burden of disease, the impact of pathogen and vector genetic diversity and a broad array of immune and inflammation phenotypes that can be mapped as traits in health and disease. Insights from study of inborn errors of immunity and multi-omics profiling together with developments in analytical methods are further advancing our knowledge of this important area.

Infectious diseases are an ever-present global threat. In this Review, Kwok, Mentzer and Knight discuss our latest understanding of how human genetics influence susceptibility to disease. Furthermore, they discuss emerging progress in the interplay between host and pathogen genetics, molecular responses to infection and vaccination, and opportunities to bring these aspects together for rapid responses to emerging diseases such as COVID-19.

Complex genetic dependencies among growth and neurological phenotypes in healthy children: Towards deciphering developmental mechanisms

The genetic mechanisms of childhood development in its many facets remain largely undeciphered. In the population of healthy infants studied in the Growing Up in Singapore Towards Healthy Outcomes (GUSTO) program, we have identified a range of dependencies among the observed phenotypes of fetal and early childhood growth, neurological development, and a number of genetic variants. We have quantified these dependencies using our information theory-based methods. The genetic variants show dependencies with single phenotypes as well as pleiotropic effects on more than one phenotype and thereby point to a large number of brain-specific and brain-expressed gene candidates. These dependencies provide a basis for connecting a range of variants with a spectrum of phenotypes (pleiotropy) as well as with each other. A broad survey of known regulatory expression characteristics, and other function-related information from the literature for these sets of candidate genes allowed us to assemble an integrated body of evidence, including a partial regulatory network, that points towards the biological basis of these general dependencies. Notable among the implicated loci are RAB11FIP4 (next to NF1), MTMR7 and PLD5, all highly expressed in the brain; DNMT1 (DNA methyl transferase), highly expressed in the placenta; and PPP1R12B and DMD (dystrophin), known to be important growth and development genes. While we cannot specify and decipher the mechanisms responsible for the phenotypes in this study, a number of connections for further investigation of fetal and early childhood growth and neurological development are indicated. These results and this approach open the door to new explorations of early human development.

Association of UBE3C Variants with Reduced Kidney Function in Patients with Diabetic Kidney Disease

Diabetic kidney disease (DKD) is the leading cause of morbidity and mortality in patients with diabetes mellitus (DM) and the most common variant of end-stage renal disease (ESRD) globally. The economic burden of ESRD treatment with dialysis is substantial. The incidence and prevalence of ESRD in Taiwan remain the highest worldwide. Therefore, identifying genetic factors affecting kidney function would have valuable clinical implications. We performed microarray experiments and identified that ubiquitin protein ligase E3C (UBE3C) is differentially expressed in two DKD patient groups with extreme (low and high) urine protein-to-creatinine ratios. A follow-up genotyping study was performed in a larger group to investigate any specific variants of UBE3C associated with DKD. A total of 263 patients were included in the study, comprising 172 patients with DKD and 91 control subjects (patients with DM without chronic kidney disease (CKD)). Two UBE3C variants (rs3802129(AA) and rs7807(CC)) were determined to be associated with reduced kidney function. The haplotype analysis revealed that rs3802129/rs3815217 (block 1) with A/G haplotype and rs8101/rs7807 (block 2) with T/C haplotype were associated with higher risks of CKD phenotypes. These findings suggest a clinical role of UBE3C variants in DKD risk.

Integrated genomics analysis highlights important SNPs and genes implicated in moderate-to-severe asthma based on GWAS and eQTL datasets

Background

Severe asthma is a chronic disease contributing to disproportionate disease morbidity and mortality. From the year of 2007, many genome-wide association studies (GWAS) have documented a large number of asthma-associated genetic variants and related genes. Nevertheless, the molecular mechanism of these identified variants involved in asthma or severe asthma risk remains largely unknown.

Methods

In the current study, we systematically integrated 3 independent expression quantitative trait loci (eQTL) data (N = 1977) and a large-scale GWAS summary data of moderate-to-severe asthma (N = 30,810) by using the Sherlock Bayesian analysis to identify whether expression-related variants contribute risk to severe asthma. Furthermore, we performed various bioinformatics analyses, including pathway enrichment analysis, PPI network enrichment analysis, in silico permutation analysis, DEG analysis and co-expression analysis, to prioritize important genes associated with severe asthma.

Results

In the discovery stage, we identified 1129 significant genes associated with moderate-to-severe asthma by using the Sherlock Bayesian analysis. Two hundred twenty-eight genes were prominently replicated by using MAGMA gene-based analysis. These 228 replicated genes were enriched in 17 biological pathways including antigen processing and presentation (Corrected P = 4.30 × 10^− 6), type I diabetes mellitus (Corrected P = 7.09 × 10^− 5), and asthma (Corrected P = 1.72 × 10^− 3). With the use of a series of bioinformatics analyses, we highlighted 11 important genes such as GNGT2, TLR6, and TTC19 as authentic risk genes associated with moderate-to-severe/severe asthma. With respect to GNGT2, there were 3 eSNPs of rs17637472 (P_eQTL = 2.98 × 10^− 8 and P_GWAS = 3.40 × 10^− 8), rs11265180 (P_eQTL = 6.0 × 10^− 6 and P_GWAS = 1.99 × 10^− 3), and rs1867087 (P_eQTL = 1.0 × 10^− 4 and P_GWAS = 1.84 × 10^− 5) identified. In addition, GNGT2 is significantly expressed in severe asthma compared with mild-moderate asthma (P = 0.045), and Gngt2 shows significantly distinct expression patterns between vehicle and various glucocorticoids (Anova P = 1.55 × 10^− 6).

Conclusions

Our current study provides multiple lines of evidence to support that these 11 identified genes as important candidates implicated in the pathogenesis of severe asthma.

Integrative genomics analysis identifies promising SNPs and genes implicated in tuberculosis risk based on multiple omics datasets

More than 10 GWASs have reported numerous genetic loci associated with tuberculosis (TB). However, the functional effects of genetic variants on TB remains largely unknown. In the present study, by combining a reported GWAS summary dataset (N = 452,264) with 3 independent eQTL datasets (N = 2,242) and other omics datasets downloaded from public databases, we conducted an integrative genomics analysis to highlight SNPs and genes implicated in TB risk. Based on independent biological and technical validations, we prioritized 26 candidate genes with eSNPs significantly associated with gene expression and TB susceptibility simultaneously; such as, CDC16 (rs7987202, rs9590408, and rs948182) and RCN3 (rs2946863, rs2878342, and rs3810194). Based on the network-based enrichment analysis, we found these 26 highlighted genes were jointly connected to exert effects on TB susceptibility. The co-expression patterns among these 26 genes were remarkably changed according to Mycobacterium tuberculosis (MTB) infection status. Based on 4 independent gene expression datasets, 21 of 26 genes (80.77%) showed significantly differential expressions between TB group and control group in mesenchymal stem cells, mice blood and lung tissues, as well as human alveolar macrophages. Together, we provide robust evidence to support 26 highlighted genes as important candidates for TB.

Prevention of adverse drug effects: a pharmacogenomic approach

PURPOSE OF REVIEW

Adverse drug reactions (ADRs) are a serious burden and can negatively impact patient quality of life. One of these ADRs, anthracycline-induced cardiotoxicity (ACT), occurs in up to 65% of treated patients and can lead to congestive heart failure. Pharmacogenetic studies have helped to reveal the mechanisms of ACT and, consequently, inform current strategies to prevent ACT in the clinic.

RECENT FINDINGS

Many pharmacogenetic studies have been conducted for ACT, but few have led to the development of clinical practice guidelines and clinical genetic testing for ACT. This is, in part, because of lack of replication in independent patient cohorts and/or validation of an affected biological pathway. Recent advances in pharmacogenetic studies have been made through the use of novel methods that directly implicate dysregulated genes and perturbed biological pathways in response to anthracycline treatment.

SUMMARY

Furthering the understanding of the genetics and altered biological pathways of ACT through these novel methods can inform clinical treatment strategies and enable refinement of current clinical practice guidelines. This can therefore lead to improvement in clinical pharmacogenetic testing for further reduction of the incidence of ACT in pediatric cancer patients taking anthracyclines.