Next generation sequencing for HLA loci in full heritage Pima Indians of Arizona, Part II: HLA-A, -B, and -C with selected non-classical loci at 4-field resolution from whole genome sequences.

While the samples and data from the Pima Indians of the Gila River Indian Community have been included in many international HLA workshops and conferences and have been the focus of numerous population reports and the source of novel alleles at the classical HLA loci, they have not been studied for the non-classical loci. In order to expand our HLA-disease association studies, we typed over 300 whole genome sequences from full Pima heritage members, controlled for first degree relationship, and employed recently developed computer algorithms to resolve HLA alleles. Both classical-HLA-A, -B, and -C- and non-classical- HLA-E, -F, -G, -J, -L, -W, -Y, -DPA2, -DPB2, -DMA, -DMB, -DOA, -DRB2, -DRB9, TAP1- loci were typed at the 4-field level of resolution. We present allele and selected haplotype frequencies, test the genotype distributions for population structure, discuss the issues that are created for tests of Hardy-Weinberg equilibrium over the four sample spaces of high resolution HLA typing, and address the implications for the evolution of non-classical pseudogenes that are no longer expressed in a phenotype subject to natural selection.

Mutation spectrum of NOD2 reveals recessive inheritance as a main driver of Early Onset Crohn's Disease.

Inflammatory bowel disease (IBD), clinically defined as Crohn's disease (CD), ulcerative colitis (UC), or IBD-unclassified, results in chronic inflammation of the gastrointestinal tract in genetically susceptible hosts. Pediatric onset IBD represents ≥ 25% of all IBD diagnoses and often presents with intestinal stricturing, perianal disease, and failed response to conventional treatments. NOD2 was the first and is the most replicated locus associated with adult IBD, to date. However, its role in pediatric onset IBD is not well understood. We performed whole-exome sequencing on a cohort of 1,183 patients with pediatric onset IBD (ages 0-18.5 years). We identified 92 probands with biallelic rare and low frequency NOD2 variants accounting for approximately 8% of our cohort, suggesting a Mendelian inheritance pattern of disease. Additionally, we investigated the contribution of recessive inheritance of NOD2 alleles in adult IBD patients from a large clinical population cohort. We found that recessive inheritance of NOD2 variants explains ~ 7% of cases in this adult IBD cohort, including ~ 10% of CD cases, confirming the observations from our pediatric IBD cohort. Exploration of EHR data showed that several of these adult IBD patients obtained their initial IBD diagnosis before 18 years of age, consistent with early onset disease. While it has been previously reported that carriers of more than one NOD2 risk alleles have increased susceptibility to Crohn's Disease (CD), our data formally demonstrate that recessive inheritance of NOD2 alleles is a mechanistic driver of early onset IBD, specifically CD, likely due to loss of NOD2 protein function. Collectively, our findings show that recessive inheritance of rare and low frequency deleterious NOD2 variants account for 7-10% of CD cases and implicate NOD2 as a Mendelian disease gene for early onset Crohn's Disease.

Exome Sequencing of 21 Bardet-Biedl Syndrome (BBS) Genes to Identify Obesity Variants in 6,851 American Indians.

OBJECTIVE: In an ongoing effort to identify the genetic variation that contributes to obesity in American Indians, known Bardet-Biedl syndrome (BBS) genes were analyzed for an effect on BMI and leptin signaling. METHODS: Potentially deleterious variants (Combined Annotation Dependent Depletion score > 20) in BBS genes were identified in whole-exome sequence data from 6,851 American Indians informative for BMI. Common variants (detected in ≥ 10 individuals) were analyzed for association with BMI; rare variants (detected in < 10 individuals) were analyzed for mean BMI of carriers. Functional assessment of variants' effect on signal transducer and activator of transcription 3 (STAT3) activity was performed in vitro. RESULTS: One common variant, rs59252892 (Thr549Ile) in BBS9, was associated with BMI (P = 0.0008, ß = 25% increase per risk allele). Among rare variants for which carriers had severe obesity (mean BMI > 40 kg/m2 ), four were in BBS9. In vitro analysis of BBS9 found the Ile allele at Thr549Ile had a 20% increase in STAT3 activity compared with the Thr allele (P = 0.01). Western blot analysis showed the Ile allele had a 15% increase in STAT3 phosphorylation (P = 0.006). Comparable functional results were observed with Ser545Gly and Val209Leu but not Leu665Phe and Lys810Glu. CONCLUSIONS: Potentially functional variants in BBS genes in American Indians are reported. However, functional evidence supporting a causal role for BBS9 in obesity is inconclusive.

KCNQ1 and Long QT Syndrome in 1/45 Amish: The Road From Identification to Implementation of Culturally Appropriate Precision Medicine.

BACKGROUND: In population-based research exome sequencing, the path from variant discovery to return of results is not well established. Variants discovered by research exome sequencing have the potential to improve population health. METHODS: Population-based exome sequencing and agnostic ExWAS were performed 5521 Amish individuals. Additional phenotyping and in vitro studies enabled reclassification of a KCNQ1 variant from variant of unknown significance to pathogenic. Results were returned to participants in a community setting. RESULTS: A missense variant was identified in KCNQ1 (c.671C>T, p.T224M), a gene associated with long QT syndrome type 1, which can cause syncope and sudden cardiac death. The p.T224M variant, present in 1/45 Amish individuals is rare in the general population (1/248 566 in gnomAD) and was highly associated with QTc on electro-cardiogram (P=5.53E-24, ß=20.2 ms/allele). Because of the potential importance of this variant to the health of the population, additional phenotyping was performed in 88 p.T224M carriers and 54 noncarriers. There was stronger clinical evidence of long QT syndrome in carriers (38.6% versus 5.5%, P=0.0006), greater history of syncope (32% versus 17%, P=0.020), and higher rate of sudden cardiac death in first degree relatives

Exome Sequencing Identifies A Nonsense Variant in DAO Associated With Reduced Energy Expenditure in American Indians.

BACKGROUND: Obesity and energy expenditure (EE) are heritable and genetic variants influencing EE may contribute to the development of obesity. We sought to identify genetic variants that affect EE in American Indians, an ethnic group with high prevalence of obesity. METHODS: Whole-exome sequencing was performed in 373 healthy Pima Indians informative for 24-hour EE during energy balance. Genetic association analyses of all high-quality exonic variants (≥5 carriers) was performed, and those predicted to be damaging were prioritized. RESULTS: Rs752074397 introduces a premature stop codon (Cys264Ter) in DAO and demonstrated the strongest association for 24-hour EE, where the Ter allele associated with substantially lower 24-hour EE (mean lower by 268 kcal/d) and sleeping EE (by 135 kcal/d). The Ter allele has a frequency = 0.5% in Pima Indians, whereas is extremely rare in most other ethnic groups (frequency < 0.01%). In vitro functional analysis showed reduced protein levels for the truncated form of DAO consistent with increased protein degradation. DAO encodes D-amino acid oxidase, which is involved in dopamine synthesis which might explain its role in modulating EE. CONCLUSION: Our results indicate that a nonsense mutation in DAO may influence EE in American Indians. Identification of variants that influence energy metabolism may lead to new pathways to treat human obesity. CLINICAL TRIAL REGISTRATION NUMBER: NCT00340132.

Characterization of Exome Variants and Their Metabolic Impact in 6,716 American Indians from the Southwest US.

Applying exome sequencing to populations with unique genetic architecture has the potential to reveal novel genes and variants associated with traits and diseases. We sequenced and analyzed the exomes of 6,716 individuals from a Southwestern American Indian (SWAI) population with well-characterized metabolic traits. We found that the SWAI population has distinct allelic architecture compared to populations of European and East Asian ancestry, and there were many predicted loss-of-function (pLOF) and nonsynonymous variants that were highly enriched or private in the SWAI population. We used pLOF and nonsynonymous variants in the SWAI population to evaluate gene-burden associations of candidate genes from European genome-wide association studies (GWASs) for type 2 diabetes, body mass index, and four major plasma lipids. We found 19 significant gene-burden associations for 11 genes, providing additional evidence for prioritizing candidate effector genes of GWAS signals. Interestingly, these associations were mainly driven by pLOF and nonsynonymous variants that are unique or highly enriched in the SWAI population. Particularly, we found four pLOF or nonsynonymous variants in APOB, APOE, PCSK9, and TM6SF2 that are private or enriched in the SWAI population and associated with low-density lipoprotein (LDL) cholesterol levels. Their large estimated effects on LDL cholesterol levels suggest strong impacts on protein function and potential clinical implications of these variants in cardiovascular health. In summary, our study illustrates the utility and potential of exome sequencing in genetically unique populations, such as the SWAI population, to prioritize candidate effector genes within GWAS loci and to find additional variants in known disease genes with potential clinical impact.

Next generation sequencing and the classical HLA loci in full heritage Pima Indians of Arizona: Defining the core HLA variation for North American Paleo-Indians.

The Pima Indians of the Gila River Indian Community in Arizona have participated in a long-range study of type 2 diabetes mellitus since 1965 and have been the subject of HLA typing and population studies since the early days of serological assays. These data have been in numerous HLA workshops and conferences and have been the source of at least five novel alleles at the classical HLA loci. In recent time nearly the entire study group was subject to next generation sequencing by whole genome or exome technologies, which has allowed us to HLA type over 3000 full heritage persons with recently developed computer algorithms. We present here the results for the classical HLA Loci: HLA-A, B, C, DRA, DRB1, DRB3, DRB4, DRB5, DPA1, DPB1, DQA1, and DQB1 to the third field of resolution for synonymous alleles and type the likely four field resolution alleles from the subset of whole genome sequences. Allele frequencies, and haplotype frequencies at up to five loci, are presented as well as measures of population structure and heterozygosity. We define a core set of HLA variation that approximates the distribution for the Paleo-Indians and impute nine-locus, 4-field haplotypes that are expected to be common in full heritage peoples.

Reactivation of a developmental Bmp2 signaling center is required for therapeutic control of the murine periosteal niche.

Two decades after signals controlling bone length were discovered, the endogenous ligands determining bone width remain unknown. We show that postnatal establishment of normal bone width in mice, as mediated by bone-forming activity of the periosteum, requires BMP signaling at the innermost layer of the periosteal niche. This developmental signaling center becomes quiescent during adult life. Its reactivation however, is necessary for periosteal growth, enhanced bone strength, and accelerated fracture repair in response to bone-anabolic therapies used in clinical orthopedic settings. Although many BMPs are expressed in bone, periosteal BMP signaling and bone formation require only Bmp2 in the Prx1-Cre lineage. Mechanistically, BMP2 functions downstream of Lrp5/6 pathway to activate a conserved regulatory element upstream of Sp7 via recruitment of Smad1 and Grhl3. Consistent with our findings, human variants of BMP2 and GRHL3 are associated with increased risk of fractures.

Profiling and Leveraging Relatedness in a Precision Medicine Cohort of 92,455 Exomes.

Large-scale human genetics studies are ascertaining increasing proportions of populations as they continue growing in both number and scale. As a result, the amount of cryptic relatedness within these study cohorts is growing rapidly and has significant implications on downstream analyses. We demonstrate this growth empirically among the first 92,455 exomes from the DiscovEHR cohort and, via a custom simulation framework we developed called SimProgeny, show that these measures are in line with expectations given the underlying population and ascertainment approach. For example, within DiscovEHR we identified ∼66,000 close (first- and second-degree) relationships, involving 55.6% of study participants. Our simulation results project that >70% of the cohort will be involved in these close relationships, given that DiscovEHR scales to 250,000 recruited individuals. We reconstructed 12,574 pedigrees by using these relationships (including 2,192 nuclear families) and leveraged them for multiple applications. The pedigrees substantially improved the phasing accuracy of 20,947 rare, deleterious compound heterozygous mutations. Reconstructed nuclear families were critical for identifying 3,415 de novo mutations in ∼1,783 genes. Finally, we demonstrate the segregation of known and suspected disease-causing mutations, including a tandem duplication that occurs in LDLR and causes familial hypercholesterolemia, through reconstructed pedigrees. In summary, this work highlights the prevalence of cryptic relatedness expected among large healthcare population-genomic studies and demonstrates several analyses that are uniquely enabled by large amounts of cryptic relatedness.

A Protein-Truncating HSD17B13 Variant and Protection from Chronic Liver Disease.

BACKGROUND: Elucidation of the genetic factors underlying chronic liver disease may reveal new therapeutic targets. METHODS: We used exome sequence data and electronic health records from 46,544 participants in the DiscovEHR human genetics study to identify genetic variants associated with serum levels of alanine aminotransferase (ALT) and aspartate aminotransferase (AST). Variants that were replicated in three additional cohorts (12,527 persons) were evaluated for association with clinical diagnoses of chronic liver disease in DiscovEHR study participants and two independent cohorts (total of 37,173 persons) and with histopathological severity of liver disease in 2391 human liver samples. RESULTS: A splice variant (rs72613567:TA) in HSD17B13, encoding the hepatic lipid droplet protein hydroxysteroid 17-beta dehydrogenase 13, was associated with reduced levels of ALT (P=4.2×10-12) and AST (P=6.2×10-10). Among DiscovEHR study participants, this variant was associated with a reduced risk of alcoholic liver disease (by 42% [95% confidence interval {CI}, 20 to 58] among heterozygotes and by 53% [95% CI, 3 to 77] among homozygotes), nonalcoholic liver disease (by 17% [95% CI, 8 to 25] among heterozygotes and by 30% [95% CI, 13 to 43] among homozygotes), alcoholic cirrhosis (by 42% [95% CI, 14 to 61] among heterozygotes and by 73% [95% CI, 15 to 91] among homozygotes), and nonalcoholic cirrhosis (by 26% [95% CI, 7 to 40] among heterozygotes and by 49% [95% CI, 15 to 69] among homozygotes). Associations were confirmed in two independent cohorts. The rs72613567:TA variant was associated with a reduced risk of nonalcoholic steatohepatitis, but not steatosis, in human liver samples. The rs72613567:TA variant mitigated liver injury associated with the risk-increasing PNPLA3 p.I148M allele and resulted in an unstable and truncated protein with reduced enzymatic activity. CONCLUSIONS: A loss-of-function variant in HSD17B13 was associated with a reduced risk of chronic liver disease and of progression from steatosis to steatohepatitis. (Funded by Regeneron Pharmaceuticals and others.).

MAPPIN: a method for annotating, predicting pathogenicity and mode of inheritance for nonsynonymous variants.

Nonsynonymous single nucleotide variants (nsSNVs) constitute about 50% of known disease-causing mutations and understanding their functional impact is an area of active research. Existing algorithms predict pathogenicity of nsSNVs; however, they are unable to differentiate heterozygous, dominant disease-causing variants from heterozygous carrier variants that lead to disease only in the homozygous state. Here, we present MAPPIN (Method for Annotating, Predicting Pathogenicity, and mode of Inheritance for Nonsynonymous variants), a prediction method which utilizes a random forest algorithm to distinguish between nsSNVs with dominant, recessive, and benign effects. We apply MAPPIN to a set of Mendelian disease-causing mutations and accurately predict pathogenicity for all mutations. Furthermore, MAPPIN predicts mode of inheritance correctly for 70.3% of nsSNVs. MAPPIN also correctly predicts pathogenicity for 87.3% of mutations from the Deciphering Developmental Disorders Study with a 78.5% accuracy for mode of inheritance. When tested on a larger collection of mutations from the Human Gene Mutation Database, MAPPIN is able to significantly discriminate between mutations in known dominant and recessive genes. Finally, we demonstrate that MAPPIN outperforms CADD and Eigen in predicting disease inheritance modes for all validation datasets. To our knowledge, MAPPIN is the first nsSNV pathogenicity prediction algorithm that provides mode of inheritance predictions, adding another layer of information for variant prioritization.

Distribution and clinical impact of functional variants in 50,726 whole-exome sequences from the DiscovEHR study.

The DiscovEHR collaboration between the Regeneron Genetics Center and Geisinger Health System couples high-throughput sequencing to an integrated health care system using longitudinal electronic health records (EHRs). We sequenced the exomes of 50,726 adult participants in the DiscovEHR study to identify ~4.2 million rare single-nucleotide variants and insertion/deletion events, of which ~176,000 are predicted to result in a loss of gene function. Linking these data to EHR-derived clinical phenotypes, we find clinical associations supporting therapeutic targets, including genes encoding drug targets for lipid lowering, and identify previously unidentified rare alleles associated with lipid levels and other blood level traits. About 3.5% of individuals harbor deleterious variants in 76 clinically actionable genes. The DiscovEHR data set provides a blueprint for large-scale precision medicine initiatives and genomics-guided therapeutic discovery.

Differential contribution of cis-regulatory elements to higher order chromatin structure and expression of the CFTR locus.

Higher order chromatin structure establishes domains that organize the genome and coordinate gene expression. However, the molecular mechanisms controlling transcription of individual loci within a topological domain (TAD) are not fully understood. The cystic fibrosis transmembrane conductance regulator (CFTR) gene provides a paradigm for investigating these mechanisms.CFTR occupies a TAD bordered by CTCF/cohesin binding sites within which are cell-type-selective cis-regulatory elements for the locus. We showed previously that intronic and extragenic enhancers, when occupied by specific transcription factors, are recruited to the CFTR promoter by a looping mechanism to drive gene expression. Here we use a combination of CRISPR/Cas9 editing of cis-regulatory elements and siRNA-mediated depletion of architectural proteins to determine the relative contribution of structural elements and enhancers to the higher order structure and expression of the CFTR locus. We found the boundaries of the CFTRTAD are conserved among diverse cell types and are dependent on CTCF and cohesin complex. Removal of an upstream CTCF-binding insulator alters the interaction profile, but has little effect on CFTR expression. Within the TAD, intronic enhancers recruit cell-type selective transcription factors and deletion of a pivotal enhancer element dramatically decreases CFTR expression, but has minor effect on its 3D structure.

Chromatin Dynamics in the Regulation of CFTR Expression.

The contribution of chromatin dynamics to the regulation of human disease-associated loci such as the cystic fibrosis transmembrane conductance regulator (CFTR) gene has been the focus of intensive experimentation for many years. Recent technological advances in the analysis of transcriptional mechanisms across the entire human genome have greatly facilitated these studies. In this review we describe the complex machinery of tissue-specific regulation of CFTR expression, and put earlier observations in context by incorporating them into datasets generated by the most recent genomics methods. Though the gene promoter is required for CFTR expression, cell-type specific regulatory elements are located elsewhere in the gene and in flanking intergenic regions. Probably within its own topological domain established by the architectural proteins CTCF and cohesin, the CFTR locus utilizes chromatin dynamics to remodel nucleosomes, recruit cell-selective transcription factors, and activate intronic enhancers. These cis-acting elements are then brought to the gene promoter by chromatin looping mechanisms, which establish long-range interactions across the locus. Despite its complexity, the CFTR locus provides a paradigm for elucidating the critical role of chromatin dynamics in the transcription of individual human genes.

FOXA2 regulates a network of genes involved in critical functions of human intestinal epithelial cells.

The forkhead box A (FOXA) family of pioneer transcription factors is critical for the development of many endoderm-derived tissues. Their importance in regulating biological processes in the lung and liver is extensively characterized, though much less is known about their role in intestine. Here we investigate the contribution of FOXA2 to coordinating intestinal epithelial cell function using postconfluent Caco2 cells, differentiated into an enterocyte-like model. FOXA2 binding sites genome-wide were determined by ChIP-seq and direct targets of the factor were validated by ChIP-qPCR and siRNA-mediated depletion of FOXA1/2 followed by RT-qPCR. Peaks of FOXA2 occupancy were frequent at loci contributing to gene ontology pathways of regulation of cell migration, cell motion, and plasma membrane function. Depletion of both FOXA1 and FOXA2 led to a significant reduction in the expression of multiple transmembrane proteins including ion channels and transporters, which form a network that is essential for maintaining normal ion and solute transport. One of the targets was the adenosine A2B receptor, and reduced receptor mRNA levels were associated with a functional decrease in intracellular cyclic AMP. We also observed that 30% of FOXA2 binding sites contained a GATA motif and that FOXA1/A2 depletion reduced GATA-4, but not GATA-6 protein levels. These data show that FOXA2 plays a pivotal role in regulating intestinal epithelial cell function. Moreover, that the FOXA and GATA families of transcription factors may work cooperatively to regulate gene expression genome-wide in the intestinal epithelium.

Architectural proteins CTCF and cohesin have distinct roles in modulating the higher order structure and expression of the CFTR locus.

Higher order chromatin structures across the genome are maintained in part by the architectural proteins CCCTC binding factor (CTCF) and the cohesin complex, which co-localize at many sites across the genome. Here, we examine the role of these proteins in mediating chromatin structure at the cystic fibrosis transmembrane conductance regulator (CFTR) gene. CFTR encompasses nearly 200 kb flanked by CTCF-binding enhancer-blocking insulator elements and is regulated by cell-type-specific intronic enhancers, which loop to the promoter in the active locus. SiRNA-mediated depletion of CTCF or the cohesin component, RAD21, showed that these two factors have distinct roles in regulating the higher order organization of CFTR. CTCF mediates the interactions between CTCF/cohesin binding sites, some of which have enhancer-blocking insulator activity. Cohesin shares this tethering role, but in addition stabilizes interactions between the promoter and cis-acting intronic elements including enhancers, which are also dependent on the forkhead box A1/A2 (FOXA1/A2) transcription factors (TFs). Disruption of the three-dimensional structure of the CFTR gene by depletion of CTCF or RAD21 increases gene expression, which is accompanied by alterations in histone modifications and TF occupancy across the locus, and causes internalization of the gene from the nuclear periphery.

Chromatin remodeling mediated by the FOXA1/A2 transcription factors activates CFTR expression in intestinal epithelial cells.

The forkhead box A transcription factors, FOXA1 and FOXA2, function as pioneer factors to open condensed chromatin and facilitate binding of other proteins. We showed previously that these factors are key components of a transcriptional network that drives enhancer function at the cystic fibrosis transmembrane conductance regulator (CFTR) locus in intestinal epithelial cells. The CFTR promoter apparently lacks tissue-specific regulatory elements and expression of the gene is controlled by multiple cis-acting elements, which coordinate gene expression in different cell types. Here we show that concurrent depletion of FOXA1 and FOXA2 represses CFTR expression and alters the three-dimensional architecture of the active locus by diminishing interactions between the promoter and intronic cis-acting elements. Reduction of FOXA1/A2 also modifies the enrichment profile of the active enhancer marks H3K27ac and H3K4me2 across the CFTR locus and alters chromatin accessibility at individual cis-elements. Moreover, loss of FOXA1/A2 suppresses the recruitment of other members of the transcriptional network including HNF1 and CDX2, to multiple cis-elements. These data reveal a complex molecular mechanism underlying the role of FOXA1/A2 in achieving high levels of CFTR expression in intestinal epithelial cells.

A genome-wide analysis of open chromatin in human epididymis epithelial cells reveals candidate regulatory elements for genes coordinating epididymal function.

The epithelium lining the epididymis has a pivotal role in ensuring a luminal environment that can support normal sperm maturation. Many of the individual genes that encode proteins involved in establishing the epididymal luminal fluid are well characterized. They include ion channels, ion exchangers, transporters, and solute carriers. However, the molecular mechanisms that coordinate expression of these genes and modulate their activities in response to biological stimuli are less well understood. To identify cis-regulatory elements for genes expressed in human epididymis epithelial cells, we generated genome-wide maps of open chromatin by DNase-seq. This analysis identified 33,542 epididymis-selective DNase I hypersensitive sites (DHS), which were not evident in five cell types of different lineages. Identification of genes with epididymis-selective DHS at their promoters revealed gene pathways that are active in immature epididymis epithelial cells. These include processes correlating with epithelial function and also others with specific roles in the epididymis, including retinol metabolism and ascorbate and aldarate metabolism. Peaks of epididymis-selective chromatin were seen in the androgen receptor gene and the cystic fibrosis transmembrane conductance regulator (CFTR) gene, which has a critical role in regulating ion transport across the epididymis epithelium. In silico prediction of transcription factor binding sites that were overrepresented in epididymis-selective DHS identified epithelial transcription factors, including ELF5 and ELF3, the androgen receptor, Pax2, and Sox9, as components of epididymis transcriptional networks. Active genes, which are targets of each transcription factor, reveal important biological processes in the epididymis epithelium.

Coordinate regulation of the gel-forming mucin genes at chromosome 11p15.5.

Four of the genes that encode gel-forming mucins, which are major components of the mucus layer protecting many epithelial surfaces, are clustered at chromosome 11p15.5 and show both cell- and tissue-specific expression patterns. We aimed to determine whether the individual genes were coordinately regulated by mechanisms involving higher order chromatin structure. CCCTC-binding factor (CTCF) sites were predicted in silico and CTCF occupancy then evaluated by chromatin immunoprecipitation. CTCF was found at many sites across the gene cluster, and its binding was correlated with mucin gene expression. Next, siRNA-mediated depletion of CTCF was shown to increase MUC2 expression in A549 lung carcinoma cells and both MUC6 and MUC5AC expression in LS180 colon carcinoma cells. These changes correlated with loss of CTCF binding at multiple sites, although others retained occupancy. In cells actively expressing the mucins, the gene cluster was shown by chromosome conformation capture to form looped three-dimensional structures with direct interactions between the MUC2 promoter region, regions 30 kb 5' to it, close to the MUC6 promoter and others near the 3' end of MUC5AC, >170 kb away. Finally, to demonstrate the importance of CTCF binding to mucin gene expression, Calu-3 lung carcinoma cells were exposed to lipopolysaccharide (LPS). LPS increased the expression of MUC2 and MUC5AC and reduced MUC5B. CTCF occupancy was concurrently depleted at specific binding sites close to these genes. These data suggest that CTCF binding and cell type-specific long-range interactions across the 11p15.5 gene cluster are critical mechanisms for coordinating gel-forming mucin gene expression.

A genome-wide analysis of open chromatin in human tracheal epithelial cells reveals novel candidate regulatory elements for lung function.

BACKGROUND: Distal cell-type-specific regulatory elements may be located at very large distances from the genes that they control and are often hidden within intergenic regions or in introns of other genes. The development of methods that enable mapping of regions of open chromatin genome wide has greatly advanced the identification and characterisation of these elements. METHODS: Here we use DNase I hypersensitivity mapping followed by deep sequencing (DNase-seq) to generate a map of open chromatin in primary human tracheal epithelial (HTE) cells and use bioinformatic approaches to characterise the distribution of these sites within the genome and with respect to gene promoters, intronic and intergenic regions. RESULTS: Genes with HTE-selective open chromatin at their promoters were associated with multiple pathways of epithelial function and differentiation. The data predict novel cell-type-specific regulatory elements for genes involved in HTE cell function, such as structural proteins and ion channels, and the transcription factors that may interact with them to control gene expression. Moreover, the map of open chromatin can identify the location of potentially critical regulatory elements in genome-wide association studies (GWAS) in which the strongest association is with single nucleotide polymorphisms in non-coding regions of the genome. We demonstrate its relevance to a recent GWAS that identifies modifiers of cystic fibrosis lung disease severity. CONCLUSION: Since HTE cells have many functional similarities with bronchial epithelial cells and other differentiated cells in the respiratory epithelium, these data are of direct relevance to elucidating the molecular basis of normal lung function and lung disease.