Interplay Of Serum Bilirubin and Tobacco Smoking with Lung and Head and Neck Cancers in a Diverse, EHR-linked Los Angeles Biobank

Background

Bilirubin is a potent antioxidant with a protective role in many diseases. We examined the relationships between serum bilirubin (SB) levels, tobacco smoking (a known cause of low SB), and aerodigestive cancers, grouped as lung cancers (LC) and head and neck cancers (HNC).

Methods

We examined the associations between SB, LC, and HNC using data from 393,210 participants from a real-world, diverse, de-identified data repository and biobank linked to the UCLA Health system. We employed regression models, propensity score matching, and polygenic scores to investigate the associations and interactions between SB, tobacco smoking, LC, and HNC.

Results

Current tobacco smokers showed lower SB (-0.04mg/dL, 95% CI: [-0.04, -0.03]), compared to never-smokers. Lower SB levels were observed in HNC and LC cases (-0.10 mg/dL, [-0.13, -0.09] and - 0.09 mg/dL, CI [-0.1, -0.07] respectively) compared to cancer-free controls with the effect persisting after adjusting for smoking. SB levels were inversely associated with HNC and LC risk (ORs per SD change in SB: 0.64, CI [0.59,0.69] and 0.57, CI [0.43,0.75], respectively). Lastly, a polygenic score (PGS) for SB was associated with LC (OR per SD change of SB-PGS: 0.71, CI [0.67, 0.76]).

Conclusions

Low SB levels are associated with an increased risk of both HNC and LC, independent of the effect of tobacco smoking. Additionally, tobacco smoking demonstrated a strong interaction with SB on LC risk. Lastly, genetically predicted low SB (using a polygenic score) is negatively associated with LC. These findings suggest that SB could serve as a potential early and low-cost biomarker for LC and HNC. The interaction with tobacco smoking suggests that smokers with lower bilirubin could likely be at higher risk for LC compared to never smokers, suggesting the utility of SB in risk stratification for patients at risk for LC. Lastly, the results of the polygenic score analyses suggest potential shared biological pathways between the genetic control of SB and the risk of LC development.

MRBEE: A novel bias-corrected multivariable Mendelian Randomization method

Mendelian randomization (MR) is an instrumental variable approach used to infer causal relationships between exposures and outcomes and can apply to summary data from genome-wide association studies (GWAS). Since GWAS summary statistics are subject to estimation errors, most existing MR approaches suffer from measurement error bias, whose scale and direction are influenced by weak instrumental variables and GWAS sample overlap, respectively. We introduce MRBEE (MR using Bias-corrected Estimating Equation), a novel multivariable MR method capable of simultaneously removing measurement error bias and identifying horizontal pleiotropy. In simulations, we showed that MRBEE is capable of effectively removing measurement error bias in the presence of weak instrumental variables and sample overlap. In two independent real data analyses, we discovered that the causal effect of BMI on coronary artery disease risk is entirely mediated by blood pressure, and that existing MR methods may underestimate the causal effect of cannabis use disorder on schizophrenia risk compared to MRBEE. MRBEE possesses significant potential for advancing genetic research by providing a valuable tool to study causality between multiple risk factors and disease outcomes, particularly as a large number of GWAS summary statistics become publicly available.

Early-life farm exposures and eczema among adults in the Agricultural Lung Health Study

Background

Several studies conducted in Europe have suggested a protective association between early-life farming exposures and childhood eczema or atopic dermatitis; few studies have examined associations in adults.

Objectives

To investigate associations between early-life exposures and eczema among 3217 adult farmers and farm spouses (mean age 62.8 years) in a case-control study nested within an US agricultural cohort.

Methods

We used sampling-weighted logistic regression to estimate odds ratios (ORs) and 95% confidence intervals (95%CIs) for associations between early-life exposures and self-reported doctor-diagnosed eczema (273 cases) and polytomous logistic regression to estimate ORs (95%CIs) for a 4-level outcome combining information on eczema and atopy (specific IgE≥0.35). Additionally, we explored genetic and gene-environment associations with eczema.

Results

Although early-life farming exposures were not associated with eczema overall, several early-life exposures were associated with a reduced risk of having both eczema and atopy. Notably, results suggest stronger protective associations among individuals with both eczema and atopy than among those with either atopy alone or eczema alone. For example, ORs (95%CIs) for having a mother who did farm work while pregnant were 1.01 (0.60-1.69) for eczema alone and 0.80 (0.65-0.99) for atopy alone, but 0.54 (0.33-0.80) for having both eczema and atopy. A genetic risk score based on previously identified atopic dermatitis variants was strongly positively associated with eczema, and interaction testing suggested protective effects of several early-life farming exposures only in individuals at lower genetic risk.

Conclusions

In utero and childhood farming exposures are associated with decreased odds of having eczema with atopy in adults.

Arsenic-gene interactions and beta-cell function in the Strong Heart Family Study

We explored arsenic-gene interactions influencing pancreatic beta-cell activity in the Strong Heart Family Study (SHFS). We considered 42 variants selected for associations with either beta-cell function (31 variants) or arsenic metabolism (11 variants) in the SHFS. Beta-cell function was calculated as homeostatic model - beta corrected for insulin resistance (cHOMA-B) by regressing homeostatic model - insulin resistance (HOMA-IR) on HOMA-B and adding mean HOMA-B. Arsenic exposure was dichotomized at the median of the sum of creatinine-corrected inorganic and organic arsenic species measured by high performance liquid chromatography-inductively coupled plasma mass spectrometry (HPLC-ICPMS). Additive GxE models for cHOMA-B were adjusted for age and ancestry, and accounted for family relationships. Models were stratified by center (Arizona, Oklahoma, North Dakota and South Dakota) and meta-analyzed. The two interactions between higher vs. lower arsenic and SNPs for cHOMA-B that were nominally significant at P < 0.05 were with rs10738708 (SNP overall effect -3.91, P = 0.56; interaction effect with arsenic -31.14, P = 0.02) and rs4607517 (SNP overall effect +16.61, P = 0.03; interaction effect with arsenic +27.02, P = 0.03). The corresponding genes GCK and TUSC1 suggest oxidative stress and apoptosis as possible mechanisms for arsenic impacts on beta-cell function. No interactions were Bonferroni-significant (1.16 × 10^-3). Our findings are suggestive of oligogenic moderation of arsenic impacts on pancreatic β-cell endocrine function, but were not Bonferroni-significant.

Epidemiology of Esophageal Squamous Cell Carcinoma

Esophageal squamous cell carcinoma (ESCC) accounts for about 90% of the 456,000 incident esophageal cancers each year. Regions of high incidence include Eastern to Central Asia, along the Rift Valley in East Africa, and into South Africa. There are many causes of ESCC, which vary among regions. Early studies in France associated smoking cigarettes and heavy alcohol consumption with high rates of ESCC, but these factors cannot explain the high incidence in other regions. We discuss other risk factors for ESCC, including polycyclic aromatic hydrocarbons from a variety of sources, high-temperature foods, diet, and oral health and the microbiome-all require further research. A growing list of defined genomic regions affects susceptibility, but large genome-wide association studies have been conducted with ethnic Chinese subjects only; more studies are called for in the rest of Asia and Africa. ESCC has been understudied, but growing infrastructure in more high-incidence countries will allow rapid progress in our understanding.

Rare, protein-truncating variants in ATM, CHEK2 and PALB2, but not XRCC2, are associated with increased breast cancer risks

BACKGROUND

Breast cancer (BC) is the most common malignancy in women and has a major heritable component. The risks associated with most rare susceptibility variants are not well estimated. To better characterise the contribution of variants in ATM, CHEK2, PALB2 and XRCC2, we sequenced their coding regions in 13 087 BC cases and 5488 controls from East Anglia, UK.

METHODS

Gene coding regions were enriched via PCR, sequenced, variant called and filtered for quality. ORs for BC risk were estimated separately for carriers of truncating variants and of rare missense variants, which were further subdivided by functional domain and pathogenicity as predicted by four in silico algorithms.

RESULTS

Truncating variants in PALB2 (OR=4.69, 95% CI 2.27 to 9.68), ATM (OR=3.26; 95% CI 1.82 to 6.46) and CHEK2 (OR=3.11; 95% CI 2.15 to 4.69), but not XRCC2 (OR=0.94; 95% CI 0.26 to 4.19) were associated with increased BC risk. Truncating variants in ATM and CHEK2 were more strongly associated with risk of oestrogen receptor (ER)-positive than ER-negative disease, while those in PALB2 were associated with similar risks for both subtypes. There was also some evidence that missense variants in ATM, CHEK2 and PALB2 may contribute to BC risk, but larger studies are necessary to quantify the magnitude of this effect.

CONCLUSIONS

Truncating variants in PALB2 are associated with a higher risk of BC than those in ATM or CHEK2. A substantial risk of BC due to truncating XRCC2 variants can be excluded.

Association study of the SLC6A3 VNTR (DAT) and DRD2/ANKK1 Taq1A polymorphisms with alcohol dependence in a population from northeastern Brazil

BACKGROUND

Alcohol dependence (AD) is a complex psychiatric disorder, affecting 5.4% of the general population lifetime, characterized by excessive alcohol consumption influenced by environmental risk factors and genetic factors. Genetic alterations in dopaminergic system are involved in the treatment and etiology of AD. The aim of this search was to test the association of the SLC6A3 40 bp-VNTR and DRD2/ANKK1 Taq1A single nucleotide polymorphism (SNP), a transporter and receptor of the dopaminergic system, with AD through a study in a population of northeastern Brazil.

METHODS

The study design was a case-control that included 227 males of northeastern Brazil (113 alcoholics and 114 controls). Alcoholics were classified according to the DSM-IV criteria for AD and controls were subjects who had nonalcohol problems or who never drank. Genotyping was detected through polymerase chain reaction (PCR) for SLC6A3 40 bp-VNTR and RFLP-PCR for DRD2/ANKK1 Taq1A, and subsequent electrophoresis on a 2% agarose gel. The distribution of allele and genotype frequencies and association of polymorphisms with AD were assessed by chi-square, Fisher's exact test, and odds ratio (OR) with a confidence interval of 95% and significance p < 0.05. Data were analyzed on BioEstat 5.3 software.

RESULTS

The SLC6A3 40 bp-VNTR was associated with AD, allelic, and genotypic frequencies were significantly different, respectively (A9 vs. A10: OR = 1.88; p = 0.01; A9/A9 vs. A10/A10: OR = 6.25; p = 0.02; A9/A9 vs. A9/A10 + A10A10: OR = 5.44; p = 0.03). However, there was no statistically significant difference when the allelic (p = 0.10) and genotypic (p > 0.05) frequencies for DRD2/ANKK1 Taq1A were compared.

CONCLUSIONS

These findings suggest that A9 allele and A9/A9 genotype of the SLC6A3 40 bp-VNTR are involved in the vulnerability to AD in the population studied. However, for the DRD2/ANKK1 SNP does not present contributions to the development of AD.

Breast and ovarian cancer risks in a large series of clinically ascertained families with a high proportion of BRCA1 and BRCA2 Dutch founder mutations

BACKGROUND

BRCA1 or BRCA2 mutations confer increased risks of breast and ovarian cancer, but risks have been found to vary across studies and populations.

METHODS

We ascertained pedigree data of 582 BRCA1 and 176 BRCA2 families and studied the variation in breast and ovarian cancer risks using a modified segregation analysis model.

RESULTS

The average cumulative breast cancer risk by age 70 years was estimated to be 45% (95% CI 36 to 52%) for BRCA1 and 27% (95% CI 14 to 38%) for BRCA2 mutation carriers. The corresponding cumulative risks for ovarian cancer were 31% (95% CI 17 to 43%) for BRCA1 and 6% (95% CI 2 to 11%) for BRCA2 mutation carriers. In BRCA1 families, breast cancer relative risk (RR) increased with more recent birth cohort (p heterogeneity = 0.0006) and stronger family histories of breast cancer (p heterogeneity < 0.001). For BRCA1, our data suggest a significant association between the location of the mutation and the ratio of breast to ovarian cancer (p<0.001). By contrast, in BRCA2 families, no evidence was found for risk heterogeneity by birth cohort, family history or mutation location.

CONCLUSIONS

BRCA1 mutation carriers conferred lower overall breast and ovarian cancer risks than reported so far, while the estimates of BRCA2 mutations were among the lowest. The low estimates for BRCA1 might be due to older birth cohorts, a moderate family history, or founder mutations located within specific regions of the gene. These results are important for a more accurate counselling of BRCA1/2 mutation carriers.

A genome-wide copy number association study of osteoporotic fractures points to the 6p25.1 locus

BACKGROUND

Osteoporosis is a systemic skeletal disease characterised by reduced bone mineral density and increased susceptibility to fracture; these traits are highly heritable. Both common and rare copy number variants (CNVs) potentially affect the function of genes and may influence disease risk.

AIM

To identify CNVs associated with osteoporotic bone fracture risk.

METHOD

We performed a genome-wide CNV association study in 5178 individuals from a prospective cohort in the Netherlands, including 809 osteoporotic fracture cases, and performed in silico lookups and de novo genotyping to replicate in several independent studies.

RESULTS

A rare (population prevalence 0.14%, 95% CI 0.03% to 0.24%) 210 kb deletion located on chromosome 6p25.1 was associated with the risk of fracture (OR 32.58, 95% CI 3.95 to 1488.89; p = 8.69 × 10(-5)). We performed an in silico meta-analysis in four studies with CNV microarray data and the association with fracture risk was replicated (OR 3.11, 95% CI 1.01 to 8.22; p = 0.02). The prevalence of this deletion showed geographic diversity, being absent in additional samples from Australia, Canada, Poland, Iceland, Denmark, and Sweden, but present in the Netherlands (0.34%), Spain (0.33%), USA (0.23%), England (0.15%), Scotland (0.10%), and Ireland (0.06%), with insufficient evidence for association with fracture risk.

CONCLUSIONS

These results suggest that deletions in the 6p25.1 locus may predispose to higher risk of fracture in a subset of populations of European origin; larger and geographically restricted studies will be needed to confirm this regional association. This is a first step towards the evaluation of the role of rare CNVs in osteoporosis.

Delineating the genetic heterogeneity of ALS using targeted high-throughput sequencing

BACKGROUND

Over 100 genes have been implicated in the aetiology of amyotrophic lateral sclerosis (ALS). A detailed understanding of their independent and cumulative contributions to disease burden may help guide various clinical and research efforts.

METHODS

Using targeted high-throughput sequencing, we characterised the variation of 10 Mendelian and 23 low penetrance/tentative ALS genes within a population-based cohort of 444 Irish ALS cases (50 fALS, 394 sALS) and 311 age-matched and geographically matched controls.

RESULTS

Known or potential high-penetrance ALS variants were identified within 17.1% of patients (38% of fALS, 14.5% of sALS). 12.8% carried variants of Mendelian disease genes (C9orf72 8.78%; SETX 2.48%; ALS2 1.58%; FUS 0.45%; TARDBP 0.45%; OPTN 0.23%; VCP 0.23%. ANG, SOD1, VAPB 0%), 4.7% carried variants of low penetrance/tentative ALS genes and 9.7% (30% of fALS, 7.1% of sALS) carried previously described ALS variants (C9orf72 8.78%; FUS 0.45%; TARDBP 0.45%). 1.6% of patients carried multiple known/potential disease variants, including all identified carriers of an established ALS variant (p<0.01); TARDBP:c.859G>A(p.[G287S]) (n=2/2 sALS). Comparison of our results with those from studies of other European populations revealed significant differences in the spectrum of disease variation (p=1.7×10(-4)).

CONCLUSIONS

Up to 17% of Irish ALS cases may carry high-penetrance variants within the investigated genes. However, the precise nature of genetic susceptibility differs significantly from that reported within other European populations. Certain variants may not cause disease in isolation and concomitant analysis of disease genes may prove highly important.