Publications by Year: 2021

Central obesity is a leading health concern with a great burden carried by ethnic minority populations, especially Hispanics/Latinos. Genetic factors contribute to the obesity burden overall and to inter-population differences. We aimed to identify the loci associated with central adiposity measured as waist-to-hip ratio (WHR), waist circumference (WC) and hip circumference (HIP) adjusted for body mass index (adjBMI) by using the Hispanic Community Health Study/Study of Latinos (HCHS/SOL); determine if differences in associations differ by background group within HCHS/SOL and determine whether previously reported associations generalize to HCHS/SOL. Our analyses included 7472 women and 5200 men of mainland (Mexican, Central and South American) and Caribbean (Puerto Rican, Cuban and Dominican) background residing in the USA. We performed genome-wide association analyses stratified and combined across sexes using linear mixed-model regression. We identified 16 variants for waist-to-hip ratio adjusted for body mass index (WHRadjBMI), 22 for waist circumference adjusted for body mass index (WCadjBMI) and 28 for hip circumference adjusted for body mass index (HIPadjBMI), which reached suggestive significance (P < 1 × 10-6). Many loci exhibited differences in strength of associations by ethnic background and sex. We brought a total of 66 variants forward for validation in cohorts (N = 34 161) with participants of Hispanic/Latino, African and European descent. We confirmed four novel loci (P < 0.05 and consistent direction of effect, and P < 5 × 10-8 after meta-analysis), including two for WHRadjBMI (rs13301996, rs79478137); one for WCadjBMI (rs3168072) and one for HIPadjBMI (rs28692724). Also, we generalized previously reported associations to HCHS/SOL, (8 for WHRadjBMI, 10 for WCadjBMI and 12 for HIPadjBMI). Our study highlights the importance of large-scale genomic studies in ancestrally diverse Hispanic/Latino populations for identifying and characterizing central obesity susceptibility that may be ancestry-specific.

Arsenic exposure has been linked to poor pulmonary function, and inefficient arsenic metabolizers may be at increased risk. Dietary rice has recently been identified as a possible substantial route of exposure to arsenic, and it remains unknown whether it can provide a sufficient level of exposure to affect pulmonary function in inefficient metabolizers. Within 12,609 participants of HCHS/SOL, asthma diagnoses and spirometry-based measures of pulmonary function were assessed, and rice consumption was inferred from grain intake via a food frequency questionnaire. After stratifying by smoking history, the relationship between arsenic metabolism efficiency [percentages of inorganic arsenic (%iAs), monomethylarsenate (%MMA), and dimethylarsinate (%DMA) species in urine] and the measures of pulmonary function were estimated in a two-sample Mendelian randomization approach (genotype information from an Illumina HumanOmni2.5-8v1-1 array), focusing on participants with high inferred rice consumption. Among never-smoking high inferred consumers of rice (n = 1395), inefficient metabolism was associated with past asthma diagnosis and forced vital capacity below the lower limit of normal (LLN) (OR 1.40, p = 0.0212 and OR 1.42, p = 0.0072, respectively, for each percentage-point increase in %iAs; OR 1.26, p = 0.0240 and OR 1.24, p = 0.0193 for %MMA; OR 0.87, p = 0.0209 and OR 0.87, p = 0.0123 for the marker of efficient metabolism, %DMA). Among ever-smoking high inferred consumers of rice (n = 1127), inefficient metabolism was associated with peak expiratory flow below LLN (OR 1.54, p = 0.0108/percentage-point increase in %iAs, OR 1.37, p = 0.0097 for %MMA, and OR 0.83, p = 0.0093 for %DMA). Less efficient arsenic metabolism was associated with indicators of pulmonary dysfunction among those with high inferred rice consumption, suggesting that reductions in dietary arsenic could improve respiratory health.

STUDY OBJECTIVES: In an older African-American sample (n = 231) we tested associations of the household environment and in-bed behaviors with sleep duration, efficiency, and wakefulness after sleep onset (WASO).

METHODS: Older adult participants completed a household-level sleep environment questionnaire, a sleep questionnaire, and underwent 7-day wrist actigraphy for objective measures of sleep. Perceived household environment (self-reported) was evaluated using questions regarding safety, physical comfort, temperature, noise, and light disturbances. In-bed behaviors included watching television, listening to radio/music, use of computer/tablet/phone, playing video games, reading books, and eating. To estimate the combined effect of the components in each domain (perceived household environment and in-bed behaviors), we calculated and standardized a weighted score per sleep outcome (e.g. duration, efficiency, WASO), with a higher score indicating worse conditions. The weights were derived from the coefficients of each component estimated from linear regression models predicting each sleep outcome while adjusting for covariates.

RESULTS: A standard deviation increase in an adverse household environment score was associated with lower self-reported sleep duration (β = -13.9 min, 95% confidence interval: -26.1, -1.7) and actigraphy-based sleep efficiency (β = -0.7%, -1.4, 0.0). A standard deviation increase in the in-bed behaviors score was associated with lower actigraphy-based sleep duration (β = -9.7 min, -18.0, -1.3), sleep efficiency (β = -1.2%, -1.9, -0.6), and higher WASO (5.3 min, 2.1, 8.6).

CONCLUSION: Intervening on the sleep environment, including healthy sleep practices, may improve sleep duration and continuity among African-Americans.

Whole-genome sequencing (WGS) and whole-exome sequencing studies have become increasingly available and are being used to identify rare genetic variants associated with health and disease outcomes. Investigators routinely use mixed models to account for genetic relatedness or other clustering variables (e.g., family or household) when testing genetic associations. However, no existing tests of the association of a rare variant with a binary outcome in the presence of correlated data control the type 1 error where there are (1) few individuals harboring the rare allele, (2) a small proportion of cases relative to controls, and (3) covariates to adjust for. Here, we address all three issues in developing a framework for testing rare variant association with a binary trait in individuals harboring at least one risk allele. In this framework, we estimate outcome probabilities under the null hypothesis and then use them, within the individuals with at least one risk allele, to test variant associations. We extend the BinomiRare test, which was previously proposed for independent observations, and develop the Conway-Maxwell-Poisson (CMP) test and study their properties in simulations. We show that the BinomiRare test always controls the type 1 error, while the CMP test sometimes does not. We then use the BinomiRare test to test the association of rare genetic variants in target genes with small-vessel disease (SVD) stroke, short sleep, and venous thromboembolism (VTE), in whole-genome sequence data from the Trans-Omics for Precision Medicine (TOPMed) program.

BACKGROUND: Sleep-disordered breathing is a common disorder associated with significant morbidity. The genetic architecture of sleep-disordered breathing remains poorly understood. Through the NHLBI Trans-Omics for Precision Medicine (TOPMed) program, we performed the first whole-genome sequence analysis of sleep-disordered breathing.

METHODS: The study sample was comprised of 7988 individuals of diverse ancestry. Common-variant and pathway analyses included an additional 13,257 individuals. We examined five complementary traits describing different aspects of sleep-disordered breathing: the apnea-hypopnea index, average oxyhemoglobin desaturation per event, average and minimum oxyhemoglobin saturation across the sleep episode, and the percentage of sleep with oxyhemoglobin saturation < 90%. We adjusted for age, sex, BMI, study, and family structure using MMSKAT and EMMAX mixed linear model approaches. Additional bioinformatics analyses were performed with MetaXcan, GIGSEA, and ReMap.

RESULTS: We identified a multi-ethnic set-based rare-variant association (p = 3.48 × 10-8) on chromosome X with ARMCX3. Additional rare-variant associations include ARMCX3-AS1, MRPS33, and C16orf90. Novel common-variant loci were identified in the NRG1 and SLC45A2 regions, and previously associated loci in the IL18RAP and ATP2B4 regions were associated with novel phenotypes. Transcription factor binding site enrichment identified associations with genes implicated with respiratory and craniofacial traits. Additional analyses identified significantly associated pathways.

CONCLUSIONS: We have identified the first gene-based rare-variant associations with objectively measured sleep-disordered breathing traits. Our results increase the understanding of the genetic architecture of sleep-disordered breathing and highlight associations in genes that modulate lung development, inflammation, respiratory rhythmogenesis, and HIF1A-mediated hypoxic response.

Fine-mapping to plausible causal variation may be more effective in multi-ancestry cohorts, particularly in the MHC, which has population-specific structure. To enable such studies, we constructed a large (n = 21,546) HLA reference panel spanning five global populations based on whole-genome sequences. Despite population-specific long-range haplotypes, we demonstrated accurate imputation at G-group resolution (94.2%, 93.7%, 97.8% and 93.7% in admixed African (AA), East Asian (EAS), European (EUR) and Latino (LAT) populations). Applying HLA imputation to genome-wide association study data for HIV-1 viral load in three populations (EUR, AA and LAT), we obviated effects of previously reported associations from population-specific HIV studies and discovered a novel association at position 156 in HLA-B. We pinpointed the MHC association to three amino acid positions (97, 67 and 156) marking three consecutive pockets (C, B and D) within the HLA-B peptide-binding groove, explaining 12.9% of trait variance.

Mitochondrial DNA (mtDNA) is present in multiple copies in human cells. We evaluated cross-sectional associations of whole blood mtDNA copy number (CN) with several cardiometabolic disease traits in 408,361 participants of multiple ancestries in TOPMed and UK Biobank. Age showed a threshold association with mtDNA CN: among younger participants (<65 years of age), each additional 10 years of age was associated with 0.03 standard deviation (s.d.) higher level of mtDNA CN (P = 0.0014) versus a 0.14 s.d. lower level of mtDNA CN (P = 1.82 × 10-13) among older participants (≥65 years). At lower mtDNA CN levels, we found age-independent associations with increased odds of obesity (P = 5.6 × 10-238), hypertension (P = 2.8 × 10-50), diabetes (P = 3.6 × 10-7), and hyperlipidemia (P = 6.3 × 10-5). The observed decline in mtDNA CN after 65 years of age may be a key to understanding age-related diseases.

Mitochondrial DNA (mtDNA) is present in multiple copies in human cells. We evaluated cross-sectional associations of whole blood mtDNA copy number (CN) with several cardiometabolic disease traits in 408,361 participants of multiple ancestries in TOPMed and UK Biobank. Age showed a threshold association with mtDNA CN: among younger participants (<65 years of age), each additional 10 years of age was associated with 0.03 standard deviation (s.d.) higher level of mtDNA CN (P = 0.0014) versus a 0.14 s.d. lower level of mtDNA CN (P = 1.82 × 10-13) among older participants (≥65 years). At lower mtDNA CN levels, we found age-independent associations with increased odds of obesity (P = 5.6 × 10-238), hypertension (P = 2.8 × 10-50), diabetes (P = 3.6 × 10-7), and hyperlipidemia (P = 6.3 × 10-5). The observed decline in mtDNA CN after 65 years of age may be a key to understanding age-related diseases.

Mitochondrial DNA (mtDNA) is present in multiple copies in human cells. We evaluated cross-sectional associations of whole blood mtDNA copy number (CN) with several cardiometabolic disease traits in 408,361 participants of multiple ancestries in TOPMed and UK Biobank. Age showed a threshold association with mtDNA CN: among younger participants (<65 years of age), each additional 10 years of age was associated with 0.03 standard deviation (s.d.) higher level of mtDNA CN (P = 0.0014) versus a 0.14 s.d. lower level of mtDNA CN (P = 1.82 × 10-13) among older participants (≥65 years). At lower mtDNA CN levels, we found age-independent associations with increased odds of obesity (P = 5.6 × 10-238), hypertension (P = 2.8 × 10-50), diabetes (P = 3.6 × 10-7), and hyperlipidemia (P = 6.3 × 10-5). The observed decline in mtDNA CN after 65 years of age may be a key to understanding age-related diseases.