Genetic Variation and Sensory Perception of a Pediatric Formulation of Ibuprofen: Can a Medicine Taste Too Good for Some?

There is wide variation in how individuals perceive the chemosensory attributes of liquid formulations of ibuprofen, encompassing both adults and children. To understand personal variation in the taste and chemesthesis properties of this medicine, and how to measure it, our first scientific strategy centered on utilizing trained adult panelists, due to the complex and time-consuming psychophysical tasks needed at this initial stage. We conducted a double-blind cohort study in which panelists underwent whole-genome-wide genotyping and psychophysically evaluated an over-the-counter pediatric medicine containing ibuprofen. Associations between sensory phenotypes and genetic variation near/within irritant and taste receptor genes were determined. Panelists who experienced the urge to cough or throat sensations found the medicine less palatable and sweet, and more irritating. Perceptions varied with genetic ancestry; panelists of African genetic ancestry had fewer chemesthetic sensations, rating the medicine sweeter, less irritating, and more palatable than did those of European genetic ancestry. We discovered a novel association between TRPA1 rs11988795 and tingling sensations, independent of ancestry. We also determined for the first time that just tasting the medicine allowed predictions of perceptions after swallowing, simplifying future psychophysical studies on diverse populations of different age groups needed to understand genetic, cultural-dietary, and epigenetic factors that influence individual perceptions of palatability and, in turn, adherence and the risk of accidental ingestion.

Epinephrine evokes shortening of human airway smooth muscle cells following ß2 adrenergic receptor desensitization.

Epinephrine (EPI), an endogenous catecholamine involved in the body's fight-or-flight responses to stress, activates α1-adrenergic receptors (α1ARs) expressed on various organs to evoke a wide range of physiological functions, including vasoconstriction. In the smooth muscle of human bronchi, however, the functional role of EPI on α1ARs remains controversial. Classically, evidence suggests that EPI promotes bronchodilation by stimulating ß2-adrenergic receptors (ß2ARs). Conventionally, the selective ß2AR agonism of EPI was thought to be, in part, due to a predominance of ß2ARs and/or a sparse, or lack of α1AR activity in human airway smooth muscle (HASM) cells. Surprisingly, we find that HASM cells express a high abundance of ADRA1B (the α1AR subtype B) and identify a spontaneous "switch-like" activation of α1ARs that evokes intracellular calcium, myosin light chain phosphorylation, and HASM cell shortening. The switch-like responses, and related EPI-induced biochemical and mechanical signals, emerged upon pharmacological inhibition of ß2ARs and/or under experimental conditions that induce ß2AR tachyphylaxis. EPI-induced procontractile effects were abrogated by an α1AR antagonist, doxazosin mesylate (DM). These data collectively uncover a previously unrecognized feed-forward mechanism driving bronchospasm via two distinct classes of G protein-coupled receptors (GPCRs) and provide a basis for reexamining α1AR inhibition for the management of stress/exercise-induced asthma and/or ß2-agonist insensitivity in patients with difficult-to-control, disease subtypes.

Genome-wide association study in minority children with asthma implicates DNAH5 in bronchodilator responsiveness.

Variability in response to short-acting ß2-agonists (e.g., albuterol) among patients with asthma from diverse racial/ethnic groups may contribute to asthma disparities. We sought to identify genetic variants associated with bronchodilator response (BDR) to identify potential mechanisms of drug response and risk factors for worse asthma outcomes. Genome-wide association studies of bronchodilator response (BDR) were performed using TOPMed Whole Genome Sequencing data of the Asthma Translational Genomic Collaboration (ATGC), which corresponded to 1136 Puerto Rican, 656 Mexican and 4337 African American patients with asthma. With the population-specific GWAS results, a trans-ethnic meta-analysis was performed to identify BDR-associated variants shared across the three populations. Replication analysis was carried out in three pediatric asthma cohorts, including CAMP (Childhood Asthma Management Program; n = 560), GACRS (Genetics of Asthma in Costa Rica Study; n = 967) and HPR (Hartford-Puerto Rico; n = 417). A genome-wide significant locus (rs35661809; P = 3.61 × 10-8) in LINC02220, a non-coding RNA gene, was identified in Puerto Ricans. While this region was devoid of protein-coding genes, capture Hi-C data showed a distal interaction with the promoter of the DNAH5 gene in lung tissue. In replication analysis, the GACRS cohort yielded a nominal association (1-tailed P < 0.05). No genetic variant was associated with BDR at the genome-wide significant threshold in Mexicans and African Americans. Our findings help inform genetic underpinnings of BDR for understudied minority patients with asthma, but the limited availability of genetic data for racial/ethnic minority children with asthma remains a paramount challenge.

REALGAR: a web app of integrated respiratory omics data.

MOTIVATION: In the post genome-wide association study (GWAS) era, omics techniques have characterized information beyond genomic variants to include cell and tissue type-specific gene transcription, transcription factor binding sites, expression quantitative trait loci (eQTL) and many other biological layers. Analysis of omics data and its integration has in turn improved the functional interpretation of disease-associated genetic variants. Over 170 000 transcriptomic and epigenomic datasets corresponding to studies of various cell and tissue types under specific disease, treatment and exposure conditions are available in the Gene Expression Omnibus resource. Although these datasets are valuable to guide the design of experimental validation studies to understand the function of disease-associated genetic loci, in their raw form, they are not helpful to experimental researchers who lack adequate computational resources or experience analyzing omics data. We sought to create an integrated re-source of tissue-specific results from omics studies that is guided by disease-specific knowledge to facilitate the design of experiments that can provide biologically meaningful insights into genetic associations. RESULTS: We designed the Reducing Associations by Linking Genes and omics Results web app to provide multi-layered omics information based on results from GWAS, transcriptomic, epigenomic and eQTL studies for gene-centric analysis and visualization. With a focus on asthma datasets, the integrated omics results it contains facilitate the formulation of hypotheses related to airways disease-associated genes and can be addressed with experimental validation studies. AVAILABILITY AND IMPLEMENTATION: The REALGAR web app is available at: http://realgar.org/. The source code is available at: https://github.com/HimesGroup/realgar. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

CEBPD modulates the airway smooth muscle transcriptomic response to glucocorticoids.

BACKGROUND: CCAAT/Enhancer Binding Protein D (CEBPD), a pleiotropic glucocorticoid-responsive transcription factor, modulates inflammatory responses. Of relevance to asthma, expression of CEBPD in airway smooth muscle (ASM) increases with glucocorticoid exposure. We sought to characterize CEBPD-mediated transcriptomic responses to glucocorticoid exposure in ASM by measuring changes observed after knockdown of CEBPD and its impact on asthma-related ASM function. METHODS: Primary ASM cells derived from four donors were transfected with CEBPD or non-targeting (NT) siRNA and exposed to vehicle control, budesonide (100 nM, 18 h), TNFα (10 ng/ml, 18 h), or both budesonide and TNFα. Subsequently, RNA-Seq was used to measure gene expression levels, and pairwise differential expression results were obtained for exposures versus vehicle and knockdown versus control conditions. Weighted gene co-expression analysis was performed to identify groups of genes with similar expression patterns across the various experimental conditions (i.e., CEBPD knockdown status, exposures). RESULTS: CEBPD knockdown altered expression of 3037 genes under at least one exposure (q-value < 0.05). Co-expression analysis identified sets of 197, 152 and 290 genes that were correlated with CEBPD knockdown status, TNFα exposure status, and both, respectively. JAK-STAT signaling pathway genes, including IL6R and SOCS3, were among those influenced by both TNFα and CEBPD knockdown. Immunoblot assays revealed that budesonide-induced IL-6R protein expression and augmented IL-6-induced STAT3 phosphorylation levels were attenuated by CEBPD knockdown in ASM. CONCLUSIONS: CEBPD modulates glucocorticoid responses in ASM, in part via modulation of IL-6 receptor signaling.

Multiomics analysis identifies BIRC3 as a novel glucocorticoid response-associated gene.

BACKGROUND: Inhaled corticosteroid (ICS) response among patients with asthma is influenced by genetics, but biologically actionable insights based on associations have not been found. Various glucocorticoid response omics data sets are available to interrogate their biological effects. OBJECTIVE: We sought to identify functionally relevant ICS-response genetic associations by integrating complementary multiomics data sets. METHODS: Variants with P values less than 10-4 from a previous ICS-response genome-wide association study were reranked on the basis of integrative scores determined from (1) glucocorticoid receptor- and (2) RNA polymerase II-binding regions inferred from ChIP-Seq data for 3 airway cell types, (3) glucocorticoid response element motifs, (4) differentially expressed genes in response to glucocorticoid exposure according to 20 transcriptomic data sets, and (5) expression quantitative trait loci from GTEx. Candidate variants were tested for association with ICS response and asthma in 6 independent studies. RESULTS: Four variants had significant (q value < 0.05) multiomics integrative scores. These variants were in a locus consisting of 52 variants in high linkage disequilibrium (r2 ≥ 0.8) near glucocorticoid receptor-binding sites by the gene BIRC3. Variants were also BIRC3 expression quantitative trait loci in lung, and 2 were within/near putative glucocorticoid response element motifs. BIRC3 had increased RNA polymerase II occupancy and gene expression, with glucocorticoid exposure in 2 ChIP-Seq and 13 transcriptomic data sets. Some BIRC3 variants in the 52-variant locus were associated (P < .05) with ICS response in 3 independent studies and others with asthma in 1 study. CONCLUSIONS: BIRC3 should be prioritized for further functional studies of ICS response.

Inhibition of ABCC1 Decreases cAMP Egress and Promotes Human Airway Smooth Muscle Cell Relaxation.

In most living cells, the second-messenger roles for adenosine 3',5'-cyclic monophosphate (cAMP) are short-lived, confined to the intracellular space, and tightly controlled by the binary switch-like actions of Gαs (stimulatory G protein)-activated adenylyl cyclase (cAMP production) and cAMP-specific PDE (cAMP breakdown). Here, by using human airway smooth muscle (HASM) cells in culture as a model, we report that activation of the cell-surface ß2AR (ß2-adrenoceptor), a Gs-coupled GPCR (G protein-coupled receptor), evokes cAMP egress to the extracellular space. Increased extracellular cAMP levels ([cAMP]e) are long-lived in culture and are induced by receptor-dependent and receptor-independent mechanisms in such a way as to define a universal response class of increased intracellular cAMP levels ([cAMP]i). We find that HASM cells express multiple ATP-binding cassette (ABC) membrane transporters, with ABCC1 (ABC subfamily member C 1) being the most highly enriched transcript mapped to MRPs (multidrug resistance-associated proteins). We show that pharmacological inhibition or downregulation of ABCC1 with siRNA markedly reduces ß2AR-evoked cAMP release from HASM cells. Furthermore, inhibition of ABCC1 activity or expression decreases basal tone and increases ß-agonist-induced HASM cellular relaxation. These findings identify a previously unrecognized role for ABCC1 in the homeostatic regulation of [cAMP]i in HASM that may be conserved traits of the Gs-GPCRs (Gs-coupled family of GPCRs). Hence, the general features of this activation mechanism may uncover new disease-modifying targets in the treatment of airflow obstruction in asthma. Surprisingly, we find that serum cAMP levels are elevated in a small cohort of patients with asthma as compared with control subjects, which warrants further investigation.

Insights into glucocorticoid responses derived from omics studies.

Glucocorticoid drugs are commonly used in the treatment of several conditions, including autoimmune diseases, asthma and cancer. Despite their widespread use and knowledge of biological pathways via which they act, much remains to be learned about the cell type-specific mechanisms of glucocorticoid action and the reasons why patients respond differently to them. In recent years, human and in vitro studies have addressed these questions with genomics, transcriptomics and other omics approaches. Here, we summarize key insights derived from omics studies of glucocorticoid response, and we identify existing knowledge gaps related to mechanisms of glucocorticoid action that future studies can address.

Chaperone-mediated autophagy regulates the pluripotency of embryonic stem cells.

Embryonic stem cells can propagate indefinitely in a pluripotent state, able to differentiate into all types of specialized cells when restored to the embryo. What sustains their pluripotency during propagation remains unclear. Here, we show that core pluripotency factors OCT4 and SOX2 suppress chaperone-mediated autophagy (CMA), a selective form of autophagy, until the initiation of differentiation. Low CMA activity promotes embryonic stem cell self-renewal, whereas its up-regulation enhances differentiation. CMA degrades isocitrate dehydrogenases IDH1 and IDH2 and reduces levels of intracellular α-ketoglutarate, an obligatory cofactor for various histone and DNA demethylases involved in pluripotency. These findings suggest that CMA mediates the effect of core pluripotency factors on metabolism, shaping the epigenetic landscape of stem cells and governing the balance between self-renewal and differentiation.

Gene Expression Profiles in Children With Suspected Sepsis.

STUDY OBJECTIVE: Sepsis recognition is a clinical challenge in children. We aim to determine whether peripheral blood gene expression profiles are associated with pathogen type and sepsis severity in children with suspected sepsis. METHODS: This was a prospective pilot observational study in a tertiary pediatric emergency department with a convenience sample of children enrolled. Participants were older than 56 days and younger than 18 years, had suspected sepsis, and had not received broad-spectrum antibiotics in the previous 4 hours. Primary outcome was source pathogen, defined as confirmed bacterial source from sterile body fluid or confirmed viral source. Secondary outcome was sepsis severity, defined as maximum therapy required for shock reversal in the first 3 hospital days. We drew peripheral blood for ribonucleic acid isolation at the sepsis protocol activation, obtained gene expression measures with the GeneChip Human Gene 2.0 ST Array, and conducted differential expression analysis. RESULTS: We collected ribonucleic acid samples from a convenience sample of 122 children with suspected sepsis and 12 healthy controls. We compared the 66 children (54%) with confirmed bacterial or viral infection and found 558 differentially expressed genes, many related to interferon signaling or viral immunity. We did not find statistically significant gene expression differences in patients according to sepsis severity. CONCLUSION: The study demonstrates feasibility of evaluating gene expression profiling data in children evaluated for sepsis in the pediatric emergency department setting. Our results suggest that gene expression profiling may facilitate identification of source pathogen in children with suspected sepsis, which could ultimately lead to improved tailoring of sepsis treatment and antimicrobial stewardship.

Airway Smooth Muscle-Specific Transcriptomic Signatures of Glucocorticoid Exposure.

Glucocorticoids, commonly used asthma controller medications, decrease symptoms in most patients, but some remain symptomatic despite high-dose treatment. The physiological basis underlying the glucocorticoid response, especially in asthma patients with severe, refractory disease, is not fully understood. We sought to identify differences between the transcriptomic response of airway smooth muscle (ASM) cells derived from donors with fatal asthma and donors without asthma to glucocorticoid exposure and to compare ASM-specific changes with those observed in other cell types. In cells derived from nine donors with fatal asthma and eight donors without asthma, RNA sequencing was used to measure ASM transcriptome changes after exposure to budesonide (100 nM 24 h) or control vehicle (DMSO). Differential expression results were obtained for this dataset, as well as 13 publicly available glucocorticoid-response transcriptomic datasets corresponding to seven cell types. Specific genes were differentially expressed in response to glucocorticoid exposure (7,835 and 6,957 in ASM cells derived from donors with fatal asthma and donors without asthma, respectively; adjusted P value < 0.05). Transcriptomic changes in response to glucocorticoid exposure were similar in ASM derived from donors with fatal asthma and donors without asthma, with enriched ontological pathways that included cytokine- and chemokine-related categories. A comparison of glucocorticoid-induced changes in the nonasthma ASM transcriptome with those observed in six other cell types showed that ASM has a distinct glucocorticoid-response signature that is also present in ASM cells from donors with fatal asthma.

Facilitating Analysis of Publicly Available ChIP-Seq Data for Integrative Studies.

ChIP-Seq, a technique that allows for quantification of DNA sequences bound by transcription factors or histones, has been widely used to characterize genome-wide DNA-protein binding at baseline and induced by specific exposures. Integrating results of multiple ChIP-Seq datasets is a convenient approach to identify robust DNA- protein binding sites and determine their cell-type specificity. We developed brocade, a computational pipeline for reproducible analysis of publicly available ChIP-Seq data that creates R markdown reports containing information on datasets downloaded, quality control metrics, and differential binding results. Glucocorticoids are commonly used anti-inflammatory drugs with tissue-specific effects that are not fully understood. We demonstrate the utility of brocade via the analysis of five ChIP-Seq datasets involving glucocorticoid receptor (GR), a transcription factor that mediates glucocorticoid response, to identify cell type-specific and shared GR binding sites across the five cell types. Our results show that brocade facilitates analysis of individual ChIP-Seq datasets and comparative studies involving multiple datasets.

Genome-wide DNA methylation profiling in infants born to gestational diabetes mellitus.

BACKGROUND: Offspring exposed to gestational diabetes mellitus (GDM) are at a high risk for metabolic diseases. The mechanisms behind the association between offspring exposed to GDM in utero and an increased risk of health consequences later in life remain unclear. The aim of this study was to clarify the changes in methylation levels in the foetuses of women with GDM and to explore the possible mechanisms linking maternal GDM with a high risk of metabolic diseases in offspring later in life. METHODS: A genome-wide comparative methylome analysis on the umbilical cord blood of infants born to 30 women with GDM and 33 women with normal pregnancy was performed using Infinium HumanMethylation 450 BeadChip assays. A quantitative methylation analysis of 18 CpG dinucleotides was verified in the validation umbilical cord blood samples from 102 newborns exposed to GDM and 103 newborns who experienced normal pregnancy by MassARRAY EpiTYPER. RESULTS: A total of 4485 differentially methylated sites (DMSs), including 2150 hypermethylated sites and 2335 hypomethylated sites, with a mean ß-value difference of >0.05, were identified by the 450k array. Good agreement was observed between the massarray validation data and the 450k array data (R2 > 0.99; P < 0.0001). Thirty-seven CpGs (representing 20 genes) with a ß-value difference of > 0.15 between the GDM and healthy groups were identified and showed potential as clinical biomarkers for GDM. "hsa04940: Type I diabetes mellitus" was the most significant Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway, with a P-value = 3.20E-07 and 1.36E-02 in the hypermethylated and hypomethylated genepathway enrichment analyses, respectively. In the Gene Ontology (GO) pathway analyses, immune MHC (major histocompatibility complex)-related pathways and neuron development-related pathways were significantly enriched. CONCLUSIONS: Our results suggest that GDM has epigenetic effects on genes that are preferentially involved in the Type I diabetes mellitus pathway, immune MHC-related pathways and neuron development-related pathways, with consequences on fetal growth and development, and provide supportive evidence that DNA methylation is involved in fetal metabolic programming.

A functional splice variant associated with decreased asthma risk abolishes the ability of gasdermin B to induce epithelial cell pyroptosis.

BACKGROUND: Genetic variants in the chromosomal region 17q21 are consistently associated with asthma. However, mechanistic studies have not yet linked any of the associated variants to a function that could influence asthma, and as a result, the identity of the asthma gene(s) remains elusive. OBJECTIVES: We sought to identify and characterize functional variants in the 17q21 locus. METHODS: We used the Exome Aggregation Consortium browser to identify coding (amino acid-changing) variants in the 17q21 locus. We obtained asthma association measures for these variants in both the Genetic Epidemiology Research in Adult Health and Aging (GERA) cohort (16,274 cases and 38,269 matched controls) and the EVE Consortium study (5,303 asthma cases and 12,560 individuals). Gene expression and protein localization were determined by quantitative RT-PCR and fluorescence immunostaining, respectively. Molecular and cellular studies were performed to determine the functional effects of coding variants. RESULTS: Two coding variants (rs2305480 and rs11078928) of the gasdermin B (GSDMB) gene in the 17q21 locus were associated with lower asthma risk in both GERA (odds ratio, 0.92; P = 1.01 × 10-6) and EVE (odds ratio, 0.85; joint PEVE = 1.31 × 10-13). In GERA, rs11078928 had a minor allele frequency (MAF) of 0.45 in unaffected (nonasthmatic) controls and 0.43 in asthma cases. For European Americans in EVE, the MAF of rs2305480 was 0.45 for controls and 0.39 for cases; for all EVE subjects, the MAF was 0.32 for controls and 0.27 for cases. GSDMB is highly expressed in differentiated airway epithelial cells, including the ciliated cells. We found that, when the GSDMB protein is cleaved by inflammatory caspase-1 to release its N-terminal fragment, potent pyroptotic cell death is induced. The splice variant rs11078928 deletes the entire exon 6, which encodes 13 amino acids in the critical N-terminus, and abolishes the pyroptotic activity of the GSDMB protein. CONCLUSIONS: Our study identified a functional asthma variant in the GSDMB gene of the 17q21 locus and implicates GSDMB-mediated epithelial cell pyroptosis in pathogenesis.

Integration of Transcriptomic Data Identifies Global and Cell-Specific Asthma-Related Gene Expression Signatures.

Over 140,000 transcriptomic studies performed in healthy and diseased cell and tissue types, at baseline and after exposure to various agents, are available in public repositories. Integrating results of transcriptomic datasets has been an attractive approach to identify gene expression signatures that are more robust than those obtained for individual datasets, especially datasets with small sample size. We developed Reproducible Analysis and Validation of Expression Data (RAVED), a pipeline that facilitates the creation of R Markdown reports detailing reproducible analysis of publicly available transcriptomic data, and used it to analyze asthma and glucocorticoid response microarray and RNA-Seq datasets. Subsequently, we used three approaches to integrate summary statistics of these studies and identify cell/tissue-specific and global asthma and glucocorticoid-induced gene expression changes. Transcriptomic integration methods were incorporated into an online app called REALGAR, where end-users can specify datasets to integrate and quickly obtain results that may facilitate design of experimental studies.

Using omics approaches to understand pulmonary diseases.

Omics approaches are high-throughput unbiased technologies that provide snapshots of various aspects of biological systems and include: 1) genomics, the measure of DNA variation; 2) transcriptomics, the measure of RNA expression; 3) epigenomics, the measure of DNA alterations not involving sequence variation that influence RNA expression; 4) proteomics, the measure of protein expression or its chemical modifications; and 5) metabolomics, the measure of metabolite levels. Our understanding of pulmonary diseases has increased as a result of applying these omics approaches to characterize patients, uncover mechanisms underlying drug responsiveness, and identify effects of environmental exposures and interventions. As more tissue- and cell-specific omics data is analyzed and integrated for diverse patients under various conditions, there will be increased identification of key mechanisms that underlie pulmonary biological processes, disease endotypes, and novel therapeutics that are efficacious in select individuals. We provide a synopsis of how omics approaches have advanced our understanding of asthma, chronic obstructive pulmonary disease (COPD), acute respiratory distress syndrome (ARDS), idiopathic pulmonary fibrosis (IPF), and pulmonary arterial hypertension (PAH), and we highlight ongoing work that will facilitate pulmonary disease precision medicine.

Multiregion sequencing reveals the intratumor heterogeneity of driver mutations in TP53-driven non-small cell lung cancer.

Intratumor heterogeneity (ITH) in non-small cell lung cancer (NSCLC) may account for resistance after a period of targeted therapies because drugs destroy only a portion of tumor cells. The recognition of ITH helps identify high-risk patients to make effective treatment decisions. However, ITH studies are confounded by interpatient heterogeneity in NSCLC and a large amount of passenger mutations. To address these issues, we recruited NSCLC patients carrying TP53 mutations and selected driver mutations within recurrently mutated genes in NSCLC. A total of 12-paired normal-tumor tissues were subjected to whole-genome/whole-exome sequencing. From these, 367 non-silent mutations were selected as driver mutations and deeply sequenced in 61 intratumoral microdissections. We identified a universal prevalence of heterogeneity in all 12 tumors, indicating branched evolution. Although TP53 mutations were observed in single biopsy of all 12 tumors, most tumors consist of both TP53 mutated and non-mutated cells in separate regions within the same tumor. This suggests the late molecular timing of the acquisition of TP53 mutations; therefore, the detection of TP53 mutations in a single biopsy may simply not reflect the early malignant potential. In addition, we identified regions of loss of heterozygosity surrounding TP53 and CDKN2A mutations in tumor 711, which also exhibited heterogeneity in different regional samples. Because the ITH of driver mutations likely has clinical consequences, further efforts are needed to limit the impact of ITH and to improve therapeutic efficiency, which will benefit NSCLC patients receiving targeted treatments.

Rare variant associations with waist-to-hip ratio in European-American and African-American women from the NHLBI-Exome Sequencing Project.

Waist-to-hip ratio (WHR), a relative comparison of waist and hip circumferences, is an easily accessible measurement of body fat distribution, in particular central abdominal fat. A high WHR indicates more intra-abdominal fat deposition and is an established risk factor for cardiovascular disease and type 2 diabetes. Recent genome-wide association studies have identified numerous common genetic loci influencing WHR, but the contributions of rare variants have not been previously reported. We investigated rare variant associations with WHR in 1510 European-American and 1186 African-American women from the National Heart, Lung, and Blood Institute-Exome Sequencing Project. Association analysis was performed on the gene level using several rare variant association methods. The strongest association was observed for rare variants in IKBKB (P=4.0 × 10(-8)) in European-Americans, where rare variants in this gene are predicted to decrease WHRs. The activation of the IKBKB gene is involved in inflammatory processes and insulin resistance, which may affect normal food intake and body weight and shape. Meanwhile, aggregation of rare variants in COBLL1, previously found to harbor common variants associated with WHR and fasting insulin, were nominally associated (P=2.23 × 10(-4)) with higher WHR in European-Americans. However, these significant results are not shared between African-Americans and European-Americans that may be due to differences in the allelic architecture of the two populations and the small sample sizes. Our study indicates that the combined effect of rare variants contribute to the inter-individual variation in fat distribution through the regulation of insulin response.

Leukocyte telomere length is associated with advanced age-related macular degeneration in the Han Chinese population.

Telomeres located at the ends of chromosomes are involved in genomic stability and play a key role in various cancers and age-related diseases. Age-related macular degeneration (AMD) is a late-onset, age-associated progressive neurodegenerative disease, which includes the geographic atrophy (GA) subtype and the choroidal neovascularization (CNV) subtype. To better understand how leukocyte telomere length (LTL) is related to AMD, we conducted an association study in 197 AMD patients and 259 healthy controls using the established quantitative PCR technique. Logistic regression was performed to evaluate the association of LTL and AMD with the age-adjusted ratio of the telomere length to the copy number of a single-copy gene (T/S). Notably, we found a significant association between AMD and LTL (OR=2.24; 95% CI=1.68-3.07; P=0.0001) after adjusting for age and sex. Furthermore, the results showed a strongly significant association between the GA subtype and the LTL (OR=4.81; 95% CI=3.15-7.82; P=0.0001) after adjusting for age and sex. Our findings provide evidence of the role that LTL plays in the pathological mechanisms of AMD, mainly in the GA subgroup but not the CNV subgroup.