FGF20 and PGM2 variants are associated with childhood asthma in family-based whole-genome sequencing studies

BACKGROUND

Asthma is a heterogeneous common respiratory disease that remains poorly understood. The established genetic associations fail to explain the high estimated heritability, and the prevalence of asthma differs between populations and geographic regions. Robust association analyses incorporating different genetic ancestries and whole-genome sequencing data may identify novel genetic associations.

METHODS

We performed family-based genome-wide association analyses of childhood-onset asthma based on whole-genome sequencing (WGS) data for the 'The Genetic Epidemiology of Asthma in Costa Rica' study (GACRS) and the Childhood Asthma Management Program (CAMP). Based on parent-child trios with children diagnosed with asthma, we performed a single variant analysis using an additive and a recessive genetic model and a region-based association analysis of low-frequency and rare variants.

RESULTS

Based on 1180 asthmatic trios (894 GACRS trios and 286 CAMP trios, a total of 3540 samples with WGS data), we identified three novel genetic loci associated with childhood-onset asthma: rs4832738 on 4p14 ($P=1.72\ast{10}^{-9}$, recessive model), rs1581479 on 8p22 ($P=1.47\ast{10}^{-8}$, additive model) and rs73367537 on 10q26 ($P=1.21\ast{10}^{-8}$, additive model in GACRS only). Integrative analyses suggested potential novel candidate genes underlying these associations: PGM2 on 4p14 and FGF20 on 8p22.

CONCLUSION

Our family-based whole-genome sequencing analysis identified three novel genetic loci for childhood-onset asthma. Gene expression data and integrative analyses point to PGM2 on 4p14 and FGF20 on 8p22 as linked genes. Furthermore, region-based analyses suggest independent potential low-frequency/rare variant associations on 8p22. Follow-up analyses are needed to understand the functional mechanisms and generalizability of these associations.

Zinc finger protein 33B demonstrates sex interaction with atopy-related markers in childhood asthma

BACKGROUND

Sex differences related to immune responses can influence atopic manifestations in childhood asthma. While genome-wide association studies have investigated a sex-specific genetic architecture of the immune response, gene-by-sex interactions have not been extensively analysed for atopy-related markers including allergy skin tests, IgE and eosinophils in asthmatic children.

METHODS

We performed a genome-wide gene-by-sex interaction analysis for atopy-related markers using whole-genome sequencing data based on 889 trios from the Genetic Epidemiology of Asthma in Costa Rica Study (GACRS) and 284 trios from the Childhood Asthma Management Program (CAMP). We also tested the findings in UK Biobank participants with self-reported childhood asthma. Furthermore, downstream analyses in GACRS integrated gene expression to disentangle observed associations.

RESULTS

Single nucleotide polymorphism (SNP) rs1255383 at 10q11.21 demonstrated a genome-wide significant gene-by-sex interaction (pinteraction=9.08×10-10) for atopy (positive skin test) with opposite direction of effects between females and males. In the UK Biobank participants with a history of childhood asthma, the signal was consistently observed with the same sex-specific effect directions for high eosinophil count (pinteraction=0.0058). Gene expression of ZNF33B (zinc finger protein 33B), located at 10q11.21, was moderately associated with atopy in girls, but not in boys.

CONCLUSIONS

We report SNPs in/near a zinc finger gene as novel sex-differential loci for atopy-related markers with opposite effect directions in females and males. A potential role for ZNF33B should be studied further as an important driver of sex-divergent features of atopy in childhood asthma.

Novel recessive locus for body mass index in childhood asthma

Most genome-wide association studies of obesity and body mass index (BMI) have so far assumed an additive mode of inheritance in their analysis, although association testing supports a recessive effect for some of the established loci, for example, rs1421085 in FTO In two whole-genome sequencing (WGS) studies of children with asthma and their parents (892 Costa Rican trios and 286 North American trios), we discovered an association between a locus (rs9292139) in LOC102724122 and BMI that reaches genome-wide significance under a recessive model in the combined analysis. As the association does not achieve significance under an additive model, our finding illustrates the benefits of the recessive model in WGS analyses.

A unifying framework for rare variant association testing in family-based designs, including higher criticism approaches, SKATs, and burden tests

MOTIVATION

Analysis of rare variants in family-based studies remains a challenge. Transmission-based approaches provide robustness against population stratification, but the evaluation of the significance of test statistics based on asymptotic theory can be imprecise. Also, power will depend heavily on the choice of the test statistic and on the underlying genetic architecture of the locus, which will be generally unknown.

RESULTS

In our proposed framework, we utilize the FBAT haplotype algorithm to obtain the conditional offspring genotype distribution under the null hypothesis given the sufficient statistic. Based on this conditional offspring genotype distribution, the significance of virtually any association test statistic can be evaluated based on simulations or exact computations, without the need for asymptotic approximations. Besides standard linear burden-type statistics, this enables our approach to also evaluate other test statistics such as variance components statistics, higher criticism approaches, and maximum-single-variant-statistics, where asymptotic theory might be involved or does not provide accurate approximations for rare variant data. Based on these P-values, combined test statistics such as the aggregated Cauchy association test (ACAT) can also be utilized. In simulation studies, we show that our framework outperforms existing approaches for family-based studies in several scenarios. We also applied our methodology to a TOPMed whole-genome sequencing dataset with 897 asthmatic trios from Costa Rica.

AVAILABILITY AND IMPLEMENTATION

FBAT software is available at https://sites.google.com/view/fbatwebpage. Simulation code is available at https://github.com/julianhecker/FBAT_rare_variant_test_simulations. Whole-genome sequencing data for 'NHLBI TOPMed: The Genetic Epidemiology of Asthma in Costa Rica' is available at https://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?study_id=phs000988.v4.p1.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A novel locus for exertional dyspnoea in childhood asthma

Most children diagnosed with asthma have respiratory symptoms such as cough, dyspnoea and wheezing, which are also important markers of overall respiratory function. A decade of genome-wide association studies (GWAS) have investigated genetic susceptibility to asthma itself, but few have focused on important respiratory symptoms that characterise childhood asthma.Using whole-genome sequencing (WGS) data for 894 asthmatic trios from a Costa Rican cohort, we performed family-based association tests (FBATs) to assess the association between genetic variants and multiple asthma-relevant respiratory phenotypes: cough, phlegm, wheezing, exertional dyspnoea and exertional chest tightness. We tested whether genome-wide significant associations were replicated in two additional studies: 1) 286 asthmatic trios from the Childhood Asthma Management Program (CAMP), and 2) 2691 African American current or former smokers from the COPDGene study.In the 894 Costa Rican trios, we identified a genome-wide significant association (p=2.16×10^-9) between exertional dyspnoea and the single nucleotide polymorphism (SNP) rs10165869, located on chromosome 2q37.3, that was replicated in the CAMP cohort (p=0.023) with the same direction of association (combined p=3.28×10^-10). This association was not found in the African American participants from COPDGene. We also found suggestive evidence for an association between SNP rs10165869 and the atypical chemokine receptor 3 (ACKR3).Our finding encourages the secondary association analysis of a wider range of phenotypes that characterise respiratory symptoms in other airway diseases/studies.

Premature Menopause, Clonal Hematopoiesis, and Coronary Artery Disease in Postmenopausal Women

BACKGROUND

Premature menopause is an independent risk factor for cardiovascular disease in women, but mechanisms underlying this association remain unclear. Clonal hematopoiesis of indeterminate potential (CHIP), the age-related expansion of hematopoietic cells with leukemogenic mutations without detectable malignancy, is associated with accelerated atherosclerosis. Whether premature menopause is associated with CHIP is unknown.

METHODS

We included postmenopausal women from the UK Biobank (n=11 495) aged 40 to 70 years with whole exome sequences and from the Women's Health Initiative (n=8111) aged 50 to 79 years with whole genome sequences. Premature menopause was defined as natural or surgical menopause occurring before age 40 years. Co-primary outcomes were the presence of any CHIP and CHIP with variant allele frequency >0.1. Logistic regression tested the association of premature menopause with CHIP, adjusted for age, race, the first 10 principal components of ancestry, smoking, diabetes, and hormone therapy use. Secondary analyses considered natural versus surgical premature menopause and gene-specific CHIP subtypes. Multivariable-adjusted Cox models tested the association between CHIP and incident coronary artery disease.

RESULTS

The sample included 19 606 women, including 418 (2.1%) with natural premature menopause and 887 (4.5%) with surgical premature menopause. Across cohorts, CHIP prevalence in postmenopausal women with versus without a history of premature menopause was 8.8% versus 5.5% (P<0.001), respectively. After multivariable adjustment, premature menopause was independently associated with CHIP (all CHIP: odds ratio, 1.36 [95% 1.10-1.68]; P=0.004; CHIP with variant allele frequency >0.1: odds ratio, 1.40 [95% CI, 1.10-1.79]; P=0.007). Associations were larger for natural premature menopause (all CHIP: odds ratio, 1.73 [95% CI, 1.23-2.44]; P=0.001; CHIP with variant allele frequency >0.1: odds ratio, 1.91 [95% CI, 1.30-2.80]; P<0.001) but smaller and nonsignificant for surgical premature menopause. In gene-specific analyses, only DNMT3A CHIP was significantly associated with premature menopause. Among postmenopausal middle-aged women, CHIP was independently associated with incident coronary artery disease (hazard ratio associated with all CHIP: 1.36 [95% CI, 1.07-1.73]; P=0.012; hazard ratio associated with CHIP with variant allele frequency >0.1: 1.48 [95% CI, 1.13-1.94]; P=0.005).

CONCLUSIONS

Premature menopause, especially natural premature menopause, is independently associated with CHIP among postmenopausal women. Natural premature menopause may serve as a risk signal for predilection to develop CHIP and CHIP-associated cardiovascular disease.

Inherited causes of clonal haematopoiesis in 97,691 whole genomes

Age is the dominant risk factor for most chronic human diseases; yet the mechanisms by which aging confers this risk are largely unknown.¹ Recently, the age-related acquisition of somatic mutations in regenerating hematopoietic stem cell populations leading to clonal expansion was associated with both hematologic cancer ^2–4 and coronary heart disease⁵, a phenomenon termed ‘Clonal Hematopoiesis of Indeterminate Potential’ (CHIP).⁶ Simultaneous germline and somatic whole genome sequence analysis now provides the opportunity to identify root causes of CHIP. Here, we analyze high-coverage whole genome sequences from 97,691 participants of diverse ancestries in the NHLBI TOPMed program and identify 4,229 individuals with CHIP. We identify associations with blood cell, lipid, and inflammatory traits specific to different CHIP genes. Association of a genome-wide set of germline genetic variants identified three genetic loci associated with CHIP status, including one locus at TET2 that was African ancestry specific. In silico-informed in vitro evaluation of the TET2 germline locus identified a causal variant that disrupts a TET2 distal enhancer resulting in increased hematopoietic stem cell self-renewal. Overall, we observe that germline genetic variation shapes hematopoietic stem cell function leading to CHIP through mechanisms that are both specific to clonal hematopoiesis and shared mechanisms leading to somatic mutations across tissues.

Allelic Heterogeneity at the CRP Locus Identified by Whole-Genome Sequencing in Multi-ancestry Cohorts

Whole-genome sequencing (WGS) can improve assessment of low-frequency and rare variants, particularly in non-European populations that have been underrepresented in existing genomic studies. The genetic determinants of C-reactive protein (CRP), a biomarker of chronic inflammation, have been extensively studied, with existing genome-wide association studies (GWASs) conducted in >200,000 individuals of European ancestry. In order to discover novel loci associated with CRP levels, we examined a multi-ancestry population (n = 23,279) with WGS (∼38× coverage) from the Trans-Omics for Precision Medicine (TOPMed) program. We found evidence for eight distinct associations at the CRP locus, including two variants that have not been identified previously (rs11265259 and rs181704186), both of which are non-coding and more common in individuals of African ancestry (∼10% and ∼1% minor allele frequency, respectively, and rare or monomorphic in 1000 Genomes populations of East Asian, South Asian, and European ancestry). We show that the minor (G) allele of rs181704186 is associated with lower CRP levels and decreased transcriptional activity and protein binding in vitro, providing a plausible molecular mechanism for this African ancestry-specific signal. The individuals homozygous for rs181704186-G have a mean CRP level of 0.23 mg/L, in contrast to individuals heterozygous for rs181704186 with mean CRP of 2.97 mg/L and major allele homozygotes with mean CRP of 4.11 mg/L. This study demonstrates the utility of WGS in multi-ethnic populations to drive discovery of complex trait associations of large effect and to identify functional alleles in noncoding regulatory regions.

Association Between Titin Loss-of-Function Variants and Early-Onset Atrial Fibrillation

IMPORTANCE

Atrial fibrillation (AF) is the most common arrhythmia affecting 1% of the population. Young individuals with AF have a strong genetic association with the disease, but the mechanisms remain incompletely understood.

OBJECTIVE

To perform large-scale whole-genome sequencing to identify genetic variants related to AF.

DESIGN, SETTING, AND PARTICIPANTS

The National Heart, Lung, and Blood Institute's Trans-Omics for Precision Medicine Program includes longitudinal and cohort studies that underwent high-depth whole-genome sequencing between 2014 and 2017 in 18 526 individuals from the United States, Mexico, Puerto Rico, Costa Rica, Barbados, and Samoa. This case-control study included 2781 patients with early-onset AF from 9 studies and identified 4959 controls of European ancestry from the remaining participants. Results were replicated in the UK Biobank (346 546 participants) and the MyCode Study (42 782 participants).

EXPOSURES

Loss-of-function (LOF) variants in genes at AF loci and common genetic variation across the whole genome.

MAIN OUTCOMES AND MEASURES

Early-onset AF (defined as AF onset in persons <66 years of age). Due to multiple testing, the significance threshold for the rare variant analysis was P = 4.55 × 10-3.

RESULTS

Among 2781 participants with early-onset AF (the case group), 72.1% were men, and the mean (SD) age of AF onset was 48.7 (10.2) years. Participants underwent whole-genome sequencing at a mean depth of 37.8 fold and mean genome coverage of 99.1%. At least 1 LOF variant in TTN, the gene encoding the sarcomeric protein titin, was present in 2.1% of case participants compared with 1.1% in control participants (odds ratio [OR], 1.76 [95% CI, 1.04-2.97]). The proportion of individuals with early-onset AF who carried a LOF variant in TTN increased with an earlier age of AF onset (P value for trend, 4.92 × 10-4), and 6.5% of individuals with AF onset prior to age 30 carried a TTN LOF variant (OR, 5.94 [95% CI, 2.64-13.35]; P = 1.65 × 10-5). The association between TTN LOF variants and AF was replicated in an independent study of 1582 patients with early-onset AF (cases) and 41 200 control participants (OR, 2.16 [95% CI, 1.19-3.92]; P = .01).

CONCLUSIONS AND RELEVANCE

In a case-control study, there was a statistically significant association between an LOF variant in the TTN gene and early-onset AF, with the variant present in a small percentage of participants with early-onset AF (the case group). Further research is necessary to understand whether this is a causal relationship.

D-Dimer in African Americans: Whole Genome Sequence Analysis and Relationship to Cardiovascular Disease Risk in the Jackson Heart Study

OBJECTIVE

Plasma levels of the fibrinogen degradation product D-dimer are higher among African Americans (AAs) compared with those of European ancestry and higher among women compared with men. Among AAs, little is known of the genetic architecture of D-dimer or the relationship of D-dimer to incident cardiovascular disease.

APPROACH AND RESULTS

We measured baseline D-dimer in 4163 AAs aged 21 to 93 years from the prospective JHS (Jackson Heart Study) cohort and assessed association with incident cardiovascular disease events. In participants with whole genome sequencing data (n=2980), we evaluated common and rare genetic variants for association with D-dimer. Each standard deviation higher baseline D-dimer was associated with a 20% to 30% increased hazard for incident coronary heart disease, stroke, and all-cause mortality. Genetic variation near F3 was associated with higher D-dimer (rs2022030, β=0.284, P=3.24×10^-11). The rs2022030 effect size was nearly 3× larger among women (β=0.373, P=9.06×10^-13) than among men (β=0.135, P=0.06; P interaction =0.009). The sex by rs2022030 interaction was replicated in an independent sample of 10 808 multiethnic men and women (P interaction =0.001). Finally, the African ancestral sickle cell variant (HBB rs334) was significantly associated with higher D-dimer in JHS (β=0.507, P=1.41×10^-14), and this association was successfully replicated in 1933 AAs (P=2.3×10^-5).

CONCLUSIONS

These results highlight D-dimer as an important predictor of cardiovascular disease risk in AAs and suggest that sex-specific and African ancestral genetic effects of the F3 and HBB loci contribute to the higher levels of D-dimer among women and AAs.

Weight loss after laparoscopic adjustable gastric band and resolution of the metabolic syndrome and its components

BACKGROUND

Substantial weight loss in the setting of obesity has considerable metabolic benefits. Yet some studies have shown improvements in obesity-related metabolic comorbidities with more modest weight loss. By closely monitoring patients undergoing bariatric surgery, we aimed to determine the effects of weight loss on the metabolic syndrome and its components and determine the weight loss required for their resolution.

METHODS

We performed a prospective observational study of obese participants with metabolic syndrome (Adult Treatment Panel III criteria) who underwent laparoscopic adjustable gastric banding. Participants were assessed for all criteria of the metabolic syndrome monthly for the first 9 months, then 3-monthly until 24 months.

RESULTS

There were 89 participants with adequate longitudinal data. Baseline body mass index was 42.4±6.2 kg m^-2 with an average age was 48.2±10.7 years. There were 56 (63%) women. Resolution of the metabolic syndrome occurred in 60 of the 89 participants (67%) at 12 months and 60 of the 75 participants (80%) at 24 months. The mean weight loss when metabolic syndrome resolved was 10.9±7.7% total body weight loss (TBWL). The median weight loss at which prevalence of disease halved was 7.0% TBWL (17.5% excess weight loss (EWL)) for hypertriglyceridaemia; 11% TBWL (26.1-28% EWL) for high-density lipoprotein cholesterol and hyperglycaemia; 20% TBWL (59.5% EWL) for hypertension and 29% TBWL (73.3% EWL) for waist circumference. The odds ratio for resolution of the metabolic syndrome with 10-12.5% TBWL was 2.09 (P=0.025), with increasing probability of resolution with more substantial weight loss.

CONCLUSIONS

In obese participants with metabolic syndrome, a weight loss target of 10-12.5% TBWL (25-30% EWL) is a reasonable initial goal associated with significant odds of having metabolic benefits. If minimal improvements are seen with this initial target, additional weight loss substantially increases the probability of resolution.

Genome-wide association of familial prostate cancer cases identifies evidence for a rare segregating haplotype at 8q24.21

Previous genome-wide association studies (GWAS) of prostate cancer risk focused on cases unselected for family history and have reported over 100 significant associations. The International Consortium for Prostate Cancer Genetics (ICPCG) has now performed a GWAS of 2511 (unrelated) familial prostate cancer cases and 1382 unaffected controls from 12 member sites. All samples were genotyped on the Illumina 5M+exome single nucleotide polymorphism (SNP) platform. The GWAS identified a significant evidence for association for SNPs in six regions previously associated with prostate cancer in population-based cohorts, including 3q26.2, 6q25.3, 8q24.21, 10q11.23, 11q13.3, and 17q12. Of note, SNP rs138042437 (p = 1.7e(-8)) at 8q24.21 achieved a large estimated effect size in this cohort (odds ratio = 13.3). 116 previously sampled affected relatives of 62 risk-allele carriers from the GWAS cohort were genotyped for this SNP, identifying 78 additional affected carriers in 62 pedigrees. A test for an excess number of affected carriers among relatives exhibited strong evidence for co-segregation of the variant with disease (p = 8.5e(-11)). The majority (92 %) of risk-allele carriers at rs138042437 had a consistent estimated haplotype spanning approximately 100 kb of 8q24.21 that contained the minor alleles of three rare SNPs (dosage minor allele frequencies <1.7 %), rs183373024 (PRNCR1), previously associated SNP rs188140481, and rs138042437 (CASC19). Strong evidence for co-segregation of a SNP on the haplotype further characterizes the haplotype as a prostate cancer predisposition locus.

CONSERVATION. Genetic assignment of large seizures of elephant ivory reveals Africa's major poaching hotspots

Poaching of elephants is now occurring at rates that threaten African populations with extinction. Identifying the number and location of Africa's major poaching hotspots may assist efforts to end poaching and facilitate recovery of elephant populations. We genetically assign origin to 28 large ivory seizures (≥0.5 metric tons) made between 1996 and 2014, also testing assignment accuracy. Results suggest that the major poaching hotspots in Africa may be currently concentrated in as few as two areas. Increasing law enforcement in these two hotspots could help curtail future elephant losses across Africa and disrupt this organized transnational crime.

Mapping epistatic quantitative trait loci

BACKGROUND

How to map quantitative trait loci (QTL) with epistasis efficiently and reliably has been a persistent problem for QTL mapping analysis. There are a number of difficulties for studying epistatic QTL. Linkage can impose a significant challenge for finding epistatic QTL reliably. If multiple QTL are in linkage and have interactions, searching for QTL can become a very delicate issue. A commonly used strategy that performs a two-dimensional genome scan to search for a pair of QTL with epistasis can suffer from low statistical power and also may lead to false identification due to complex linkage disequilibrium and interaction patterns.

RESULTS

To tackle the problem of complex interaction of multiple QTL with linkage, we developed a three-stage search strategy. In the first stage, main effect QTL are searched and mapped. In the second stage, epistatic QTL that interact significantly with other identified QTL are searched. In the third stage, new epistatic QTL are searched in pairs. This strategy is based on the consideration that most genetic variance is due to the main effects of QTL. Thus by first mapping those main-effect QTL, the statistical power for the second and third stages of analysis for mapping epistatic QTL can be maximized. The search for main effect QTL is robust and does not bias the search for epistatic QTL due to a genetic property associated with the orthogonal genetic model that the additive and additive by additive variances are independent despite of linkage. The model search criterion is empirically and dynamically evaluated by using a score-statistic based resampling procedure. We demonstrate through simulations that the method has good power and low false positive in the identification of QTL and epistasis.

CONCLUSION

This method provides an effective and powerful solution to map multiple QTL with complex epistatic pattern. The method has been implemented in the user-friendly computer software Windows QTL Cartographer. This will greatly facilitate the application of the method for QTL mapping data analysis.

Dependency effects in multi-locus match probabilities

For finite populations, differences in individual histories can cause between-locus allelic dependencies even for unlinked loci. The main motivation for this study is to quantify the effect of such dependencies on genotypic match probabilities. We compare the two-locus match probability, the probability that two individuals (four gametes) chosen at random will have the same genotype at both loci, with the probability computed as the product of the one-locus match probabilities. It is demonstrated that the product rule probability always underestimates the two-locus match probability. For highly mutable minisatellite loci, these probabilities can differ by an order of magnitude or more. A simplified three-locus problem is explored, providing evidence that the degree of under-estimation worsens for more loci.

Regulation of the "tetCD" genes of transposon Tn10

In addition to the genes involved in tetracycline resistance, the loop region of the composite transposon Tn10 contains two other known genes, tetC and tetD, whose functions are unclear. Using primarily a genetic approach, we examined tetCD gene expression and regulation. The tetC gene product, TetC, is a diffusible repressor of both tetC and tetD transcription. Despite an earlier claim by others, we do not detect induction of either tetC or tetD by tetracycline (Tc) or several of its analogs. Although the 5' ends of the tetC and tetD messages overlap due to transcription from convergent promoters, we find no evidence for anti-sense RNA control. The operator for the TetC repressor has been localized. We also demonstrate that transcription from the tetD promoter probably terminates within IS10-Right and does not apparently interfere with Tn10 or IS10-Right transposition or its regulation.

Stimulus context and infant orientation discrimination

In three experiments, the effect of additional "contextual" elements on the discrimination of the orientation of linear and curvilinear segments was investigated with 4-month-old infants. In Experiment 1, paired visual matrices (one which contained some irregularity in orientation of internal elements, vs one which contained no irregularities) were presented. Infants detected irregular matrices significantly better than chance, but such detection was not aided by contextual elements. Discrimination of orientation in Experiment 2 was assessed with a paired-comparison familiarization-novelty paradigm. It was found that the addition of elements here significantly aided discrimination of linear segment orientation, but not curvilinear segment orientation. Experiment 3 investigated why this effect was not found for curvilinear segments; after equating the curvilinear stimuli to linear ones used in Experiment 2 with respect to the closedness of figure, discrimination of curvilinear orientation was observed.