Summary-data based Mendelian randomization identifies gene expression regulatory polymorphisms associated with bovine paratuberculosis by modulation of the nuclear factor Kappa ß (NF-κß)-mediated inflammatory response.

Genome-wide association studies (GWAS) have identified host genetic variants associated with paratuberculosis (PTB) susceptibility. Most of the GWAS-identified SNPs are in non-coding regions. Connecting these non-coding variants and downstream affected genes is a challenge and, up to date, only a few functional mutations or expression quantitative loci (cis-eQTLs) associated with PTB susceptibility have been identified. In the current study, the associations between imputed whole-genome sequence genotypes and whole RNA-Sequencing data from peripheral blood (PB) and ileocecal valve (ICV) samples of Spanish Holstein cows (N = 16) were analyzed with TensorQTL. This approach allowed the identification of 88 and 37 cis-eQTLs regulating the expression levels of 90 and 37 genes in PB and ICV samples, respectively (False discorey rate, FDR ≤ 0.05). Next, we applied summary-based data Mendelian randomization (SMR) to integrate the cis-eQTL dataset with GWAS data obtained from a cohort of 813 culled cattle that were classified according to the presence or absence of PTB-associated histopathological lesions in gut tissues. After multiple testing corrections (FDR ≤ 0.05), we identified two novel cis-eQTLs affecting the expression of the early growth response factor 4 (EGR4) and the bovine neuroblastoma breakpoint family member 6-like protein isoform 2 (MGC134040) that showed pleiotropic associations with the presence of multifocal and diffuse lesions in gut tissues; P = 0.002 and P = 0.017, respectively. While EGR4 acts as a brake on T-cell proliferation and cytokine production through interaction with the nuclear factor Kappa ß (NF-κß), MGC134040 is a target gene of NF-κß. Our findings provide a better understanding of the genetic factors influencing PTB outcomes, confirm that the multifocal lesions are localized/confined lesions that have different underlying host genetics than the diffuse lesions, and highlight regulatory SNPs and regulated-gene targets to design future functional studies.

A meta-analysis of pre-pregnancy maternal body mass index and placental DNA methylation identifies 27 CpG sites with implications for mother-child health.

Higher maternal pre-pregnancy body mass index (ppBMI) is associated with increased neonatal morbidity, as well as with pregnancy complications and metabolic outcomes in offspring later in life. The placenta is a key organ in fetal development and has been proposed to act as a mediator between the mother and different health outcomes in children. The overall aim of the present work is to investigate the association of ppBMI with epigenome-wide placental DNA methylation (DNAm) in 10 studies from the PACE consortium, amounting to 2631 mother-child pairs. We identify 27 CpG sites at which we observe placental DNAm variations of up to 2.0% per 10 ppBMI-unit. The CpGs that are differentially methylated in placenta do not overlap with CpGs identified in previous studies in cord blood DNAm related to ppBMI. Many of the identified CpGs are located in open sea regions, are often close to obesity-related genes such as GPX1 and LGR4 and altogether, are enriched in cancer and oxidative stress pathways. Our findings suggest that placental DNAm could be one of the mechanisms by which maternal obesity is associated with metabolic health outcomes in newborns and children, although further studies will be needed in order to corroborate these findings.

Gluten-induced RNA methylation changes regulate intestinal inflammation via allele-specific XPO1 translation in epithelial cells.

OBJECTIVES: Coeliac disease (CD) is a complex autoimmune disorder that develops in genetically susceptible individuals. Dietary gluten triggers an immune response for which the only available treatment so far is a strict, lifelong gluten free diet. Human leucocyte antigen (HLA) genes and several non-HLA regions have been associated with the genetic susceptibility to CD, but their role in the pathogenesis of the disease is still essentially unknown, making it complicated to develop much needed non-dietary treatments. Here, we describe the functional involvement of a CD-associated single-nucleotide polymorphism (SNP) located in the 5'UTR of XPO1 in the inflammatory environment characteristic of the coeliac intestinal epithelium. DESIGN: The function of the CD-associated SNP was investigated using an intestinal cell line heterozygous for the SNP, N6-methyladenosine (m6A)-related knock-out and HLA-DQ2 mice, and human samples from patients with CD. RESULTS: Individuals harbouring the risk allele had higher m6A methylation in the 5'UTR of XPO1 RNA, rendering greater XPO1 protein amounts that led to downstream nuclear factor kappa B (NFkB) activity and subsequent inflammation. Furthermore, gluten exposure increased overall m6A methylation in humans as well as in in vitro and in vivo models. CONCLUSION: We identify a novel m6A-XPO1-NFkB pathway that is activated in CD patients. The findings will prompt the development of new therapeutic approaches directed at m6A proteins and XPO1, a target under evaluation for the treatment of intestinal disorders.

Genetic Contribution of Endometriosis to the Risk of Developing Hormone-Related Cancers.

Endometriosis is a common gynecological disorder that has been associated with endometrial, breast and epithelial ovarian cancers in epidemiological studies. Since complex diseases are a result of multiple environmental and genetic factors, we hypothesized that the biological mechanism underlying their comorbidity might be explained, at least in part, by shared genetics. To assess their potential genetic relationship, we performed a two-sample mendelian randomization (2SMR) analysis on results from public genome-wide association studies (GWAS). This analysis confirmed previously reported genetic pleiotropy between endometriosis and endometrial cancer. We present robust evidence supporting a causal genetic association between endometriosis and ovarian cancer, particularly with the clear cell and endometrioid subtypes. Our study also identified genetic variants that could explain those associations, opening the door to further functional experiments. Overall, this work demonstrates the value of genomic analyses to support epidemiological data, and to identify targets of relevance in multiple disorders.

Identification of loci associated with susceptibility to bovine paratuberculosis and with the dysregulation of the MECOM, eEF1A2, and U1 spliceosomal RNA expression.

Although genome-wide association studies have identified single nucleotide polymorphisms (SNPs) associated with the susceptibility to Mycobacterium avium subsp. paratuberculosis (MAP) infection, only a few functional mutations for bovine paratuberculosis (PTB) have been characterized. Expression quantitative trait loci (eQTLs) are genetic variants typically located in gene regulatory regions that alter gene expression in an allele-specific manner. eQTLs can be considered as functional links between genomic variants, gene expression, and ultimately phenotype. In the current study, peripheral blood (PB) and ileocecal valve (ICV) gene expression was quantified by RNA-Seq from fourteen Holstein cattle with no lesions and with PTB-associated histopathological lesions in gut tissues. Genotypes were generated from the Illumina LD EuroG10K BeadChip. The associations between gene expression levels (normalized read counts) and genetic variants were analyzed by a linear regression analysis using R Matrix eQTL 2.2. This approach allowed the identification of 192 and 48 cis-eQTLs associated with the expression of 145 and 43 genes in the PB and ICV samples, respectively. To investigate potential relationships between these cis-eQTLs and MAP infection, a case-control study was performed using the genotypes for all the identified cis-eQTLs and phenotypical data (histopathology, ELISA for MAP-antibodies detection, tissue PCR, and bacteriological culture) of 986 culled cows. Our results suggested that the heterozygous genotype in the cis-eQTL-rs43744169 (T/C) was associated with the up-regulation of the MDS1 and EVI1 complex (MECOM) expression, with positive ELISA, PCR, and bacteriological culture results, and with increased risk of progression to clinical PTB. As supporting evidence, the presence of the minor allele was associated with higher MECOM levels in plasma samples from infected cows and with increased MAP survival in an ex-vivo macrophage killing assay. Moreover, the presence of the two minor alleles in the cis-eQTL-rs110345285 (C/C) was associated with the dysregulation of the eukaryotic elongation factor 1-α2 (eEF1A2) expression and with increased ELISA (OD) values. Finally, the presence of the minor allele in the cis-eQTL rs109859270 (C/T) was associated with the up-regulation of the U1 spliceosomal RNA expression and with an increased risk of progression to clinical PTB. The introduction of these novel functional variants into marker-assisted breeding programs is expected to have a relevant effect on PTB control.

Highly Reproducible Hyperthermia Response in Water, Agar, and Cellular Environment by Discretely PEGylated Magnetite Nanoparticles.

Local heat generation from magnetic nanoparticles (MNPs) exposed to alternating magnetic fields can revolutionize cancer treatment. However, the application of MNPs as anticancer agents is limited by serious drawbacks. Foremost among these are the fast uptake and biodegradation of MNPs by cells and the unpredictable magnetic behavior of the MNPs when they accumulate within or around cells and tissues. In fact, several studies have reported that the heating power of MNPs is severely reduced in the cellular environment, probably due to a combination of increased viscosity and strong NP agglomeration. Herein, we present an optimized protocol to coat magnetite (Fe3O4) NPs larger than 20 nm (FM-NPs) with high molecular weight PEG molecules that avoid collective coatings, prevent the formation of large clusters of NPs and keep constant their high heating performance in environments with very different ionic strengths and viscosities (distilled water, physiological solutions, agar and cell culture media). The great reproducibility and reliability of the heating capacity of this FM-NP@PEG system in such different environments has been confirmed by AC magnetometry and by more conventional calorimetric measurements. The explanation of this behavior has been shown to lie in preserving as much as possible the magnetic single domain-type behavior of nearly isolated NPs. In vitro endocytosis experiments in a colon cancer-derived cell line indicate that FM-NP@PEG formulations with PEGs of higher molecular weight (20 kDa) are more resistant to endocytosis than formulations with smaller PEGs (5 kDa), showing quite large uptake mean-life (τ > 5 h) in comparison with other NP systems. The in vitro magnetic hyperthermia was performed at 21 mT and 650 kHz during 1 h in a pre-endocytosis stage and complete cell death was achieved 48 h posthyperthermia. These optimal FM-NP@PEG formulations with high resistance to endocytosis and predictable magnetic response will aid the progress and accuracy of the emerging era of theranostics.

Human mitochondrial DNA is extensively methylated in a non-CpG context.

Mitochondrial dysfunction plays critical roles in cancer development and related therapeutic response; however, exact molecular mechanisms remain unclear. Recently, alongside the discovery of mitochondrial-specific DNA methyltransferases, global and site-specific methylation of the mitochondrial genome has been described. Investigation of any functional consequences however remains unclear and debated due to insufficient evidence of the quantitative degree and frequency of mitochondrial DNA (mtDNA) methylation. This study uses WGBS to provide the first quantitative report of mtDNA methylation at single base pair resolution. The data show that mitochondrial genomes are extensively methylated predominantly at non-CpG sites. Importantly, these methylation patterns display notable differences between normal and cancer cells. Furthermore, knockdown of DNA methyltransferase enzymes resulted in a marked global reduction of mtDNA methylation levels, indicating these enzymes may be associated with the establishment and/or maintenance of mtDNA methylation. DNMT3B knockdown cells displayed a comparatively pronounced global reduction in mtDNA methylation with concomitant increases in gene expression, suggesting a potential functional link between methylation and gene expression. Together these results demonstrate reproducible, non-random methylation patterns of mtDNA and challenge the notion that mtDNA is lowly methylated. This study discusses key differences in methodology that suggest future investigations must allow for techniques that assess both CpG and non-CpG methylation.

A novel RT-QPCR-based assay for the relative quantification of residue specific m6A RNA methylation.

N6-methyladenosine (m6A) is the most common and abundant RNA modification. Recent studies have shown its importance in the regulation of several biological processes, including the immune response, and different approaches have been developed in order to map and quantify m6A marks. However, site specific detection of m6A methylation has been technically challenging, and existing protocols are long and tedious and often involve next-generation sequencing. Here, we describe a simple RT-QPCR based approach for the relative quantification of candidate m6A regions that takes advantage of the diminished capacity of BstI enzyme to retrotranscribe m6A residues. Using this technique, we have been able to confirm the recently described m6A methylation in the 3'UTR of SOCS1 and SOCS3 transcripts. Moreover, using the method presented here, we have also observed alterations in the relative levels of m6A in specific motifs of SOCS genes in celiac disease patients and in pancreatic ß-cells exposed to inflammatory stimuli.

The methylome of the celiac intestinal epithelium harbours genotype-independent alterations in the HLA region.

The Human Leucocyte Antigen (HLA) locus and other DNA sequence variants identified in Genome-Wide Association (GWA) studies explain around 50% of the heritability of celiac disease (CD). However, the pathogenesis of CD could be driven by other layers of genomic information independent from sequence variation, such as DNA methylation, and it is possible that allele-specific methylation explains part of the SNP associations. Since the DNA methylation landscape is expected to be different among cell types, we analyzed the methylome of the epithelial and immune cell populations of duodenal biopsies in CD patients and controls separately. We found a cell type-specific methylation signature that includes genes mapping to the HLA region, namely TAP1 and HLA-B. We also performed Immunochip SNP genotyping of the same samples and interrogated the expression of some of the affected genes. Our analysis revealed that the epithelial methylome is characterized by the loss of CpG island (CGI) boundaries, often associated to altered gene expression, and by the increased variability of the methylation across the samples. The overlap between differentially methylated positions (DMPs) and CD-associated SNPs or variants contributing to methylation quantitative trait loci (mQTLs) is minimal. In contrast, there is a notable enrichment of mQTLs among the most significant CD-associated SNPs. Our results support the notion that DNA methylation alterations constitute a genotype-independent event and confirm its role in the HLA region (apart from the well-known, DQ allele-specific effect). Finally, we find that a fraction of the CD-associated variants could exert its phenotypic effect through DNA methylation.

A long noncoding RNA associated with susceptibility to celiac disease.

Recent studies have implicated long noncoding RNAs (lncRNAs) as regulators of many important biological processes. Here we report on the identification and characterization of a lncRNA, lnc13, that harbors a celiac disease-associated haplotype block and represses expression of certain inflammatory genes under homeostatic conditions. Lnc13 regulates gene expression by binding to hnRNPD, a member of a family of ubiquitously expressed heterogeneous nuclear ribonucleoproteins (hnRNPs). Upon stimulation, lnc13 levels are reduced, thereby allowing increased expression of the repressed genes. Lnc13 levels are significantly decreased in small intestinal biopsy samples from patients with celiac disease, which suggests that down-regulation of lnc13 may contribute to the inflammation seen in this disease. Furthermore, the lnc13 disease-associated variant binds hnRNPD less efficiently than its wild-type counterpart, thus helping to explain how these single-nucleotide polymorphisms contribute to celiac disease.

Expression analysis in intestinal mucosa reveals complex relations among genes under the association peaks in celiac disease.

Celiac disease is a chronic immune-mediated disorder with an important genetic component. To date, there are 57 independent association signals from 39 non-HLA loci, and a total of 66 candidate genes have been proposed. We aimed to scrutinize the functional implication of 45 of those genes by analyzing their expression in the disease tissue of celiac patients (at diagnosis/treatment) compared with non-celiac controls. Moreover, we investigated the SNP genotype effect in gene expression and performed coexpression analyses. Several genes showed differential expression among disease groups, most of them related to immune response. Multiple trans-eQTLs but only four cis-eQTLs were found, and surprisingly the genotype effect seems to be stimulus dependent as it differs among groups. Coexpression levels vary from higher to lower levels in active patients at diagnosis, treated patients and non-celiac controls respectively. A subset of 18 genes tightly correlated in both groups of patients but not in controls was identified. Interestingly, this subset of genes was influenced by the genotype of three SNPs. One of the SNPs, rs1018326 on chromosome two is on top of a known lincRNA whose function is not yet described, and whose expression seems to be upregulated in active disease when comparing biopsy pairs from the same individuals. Our results strongly suggest that the effects of disease-associated SNPs go far beyond the oversimplistic idea of transcriptional control at a nearby locus. Further investigations are needed to determine how each variant disrupts fine-tuning mechanisms in the genome that eventually lead to disease.

Celiac disease in the Mediterranean area.

BACKGROUND: The World Gastroenterology Organization recommends developing national guidelines for the diagnosis of Celiac Disease (CD): hence a profile of the diagnosis of CD in each country is required. We aim to describe a cross-sectional picture of the clinical features and diagnostic facilities in 16 countries of the Mediterranean basin. Since a new ESPGHAN diagnostic protocol was recently published, our secondary aim is to estimate how many cases in the same area could be identified without a small intestinal biopsy. METHODS: By a stratified cross-sectional retrospective study design, we examined clinical, histological and laboratory data from 749 consecutive unselected CD children diagnosed by national referral centers. RESULTS: The vast majority of cases were diagnosed before the age of 10 (median: 5 years), affected by diarrhea, weight loss and food refusal, as expected. Only 59 cases (7.8%) did not suffer of major complaints. Tissue transglutaminase (tTG) assay was available, but one-third of centers reported financial constraints in the regular purchase of the assay kits. 252 cases (33.6%) showed tTG values over 10 times the local normal limit. Endomysial antibodies and HLA typing were routinely available in only half of the centers. CD was mainly diagnosed from small intestinal biopsy, available in all centers. Based on these data, only 154/749 cases (20.5%) would have qualified for a diagnosis of CD without a small intestinal biopsy, according to the new ESPGHAN protocol. CONCLUSIONS: This cross-sectional study of CD in the Mediterranean referral centers offers a puzzling picture of the capacities to deal with the emerging epidemic of CD in the area, giving a substantive support to the World Gastroenterology Organization guidelines.

Detection of atomic scale changes in the free volume void size of three-dimensional colorectal cancer cell culture using positron annihilation lifetime spectroscopy.

Positron annihilation lifetime spectroscopy (PALS) provides a direct measurement of the free volume void sizes in polymers and biological systems. This free volume is critical in explaining and understanding physical and mechanical properties of polymers. Moreover, PALS has been recently proposed as a potential tool in detecting cancer at early stages, probing the differences in the subnanometer scale free volume voids between cancerous/healthy skin samples of the same patient. Despite several investigations on free volume in complex cancerous tissues, no positron annihilation studies of living cancer cell cultures have been reported. We demonstrate that PALS can be applied to the study in human living 3D cell cultures. The technique is also capable to detect atomic scale changes in the size of the free volume voids due to the biological responses to TGF-ß. PALS may be developed to characterize the effect of different culture conditions in the free volume voids of cells grown in vitro.

Coregulation and modulation of NFκB-related genes in celiac disease: uncovered aspects of gut mucosal inflammation.

It is known that the NFκB route is constitutively upregulated in celiac disease (CD), an immune-mediated disorder of the gut caused by intolerance to ingested gluten. Our aim was to scrutinize the expression patterns of several of the most biologically relevant components of the NFκB route in intestinal biopsies from active and treated patients and after in vitro gliadin challenge, and to assess normalization of the expression using an inhibitor of the MALT1 paracaspase. The expression of 93 NFκB genes was measured by RT-PCR in a set of uncultured active and treated CD and control biopsies, and in cultured biopsy series challenged with gliadin, the NFκB modulator, both compounds and none. Methylation of eight genes involved in NFκB signaling was analyzed by conventional pyrosequencing. Groups were compared and Pearson's correlation matrixes were constructed to check for coexpression and co-methylation. Our results confirm the upregulation of the NFκB pathway and show that constitutively altered genes usually belong to the core of the pathway and have central roles, whereas genes overexpressed only in active CD are more peripheral. Additionally, this is the first work to detect methylation level changes in celiac intestinal mucosa. Coexpression is very common in controls, whereas gliadin challenge and especially chronic inflammation present in untreated CD result in the disruption of the regulatory equilibrium. In contrast, co-methylation occurs more often in active CD. Importantly, NFκB modulation partially restores coregulation, opening the door to future therapeutic possibilities and targets.

Angiogenesis-related gene expression analysis in celiac disease.

Celiac disease (CD) involves disturbance of the small-bowel mucosal vascular network, and transglutaminase autoantibodies (TGA) have been related to angiogenesis disturbance, a complex phenomenon probably also influenced by common genetic variants in angiogenesis-related genes. A set of genes with "angiogenesis" GO term identified in a previous expression microarray experiment (SCG2, STAB1, TGFA, ANG, ERBB2, GNA13, PML, CASP8, ECGF1, JAG1, HIF1A, TNFSF13 and TGM2) was selected for genetic and functional studies. SNPs that showed a trend for association with CD in the first GWAS were genotyped in 555 patients and 541 controls. Gene expression of all genes was quantified in 15 pairs of intestinal biopsies (diagnosis vs. GFD) and in three-dimensional HUVEC and T84 cell cultures incubated with TGA-positive and negative serum. A regulatory SNP in TNFSF13 (rs11552708) is associated with CD (p = 0.01, OR = 0.7). Expression changes in biopsies pointed to TGM2 and PML as up-regulated antiangiogenic genes and to GNA13, TGFA, ERBB2 and SCG2 as down-regulated proangiogenic factors in CD. TGA seem to enhance TGM2 expression in both cell models, but PML expression was induced only in T84 enterocytes while GNA13 and ERBB2 were repressed in HUVEC endothelial cells, with several genes showing discordant effects in each model, highlighting the complexity of gene interactions in the pathogenesis of CD. Finally, cell culture models are useful tools to help dissect complex responses observed in human explants.

Upregulation of KIR3DL1 gene expression in intestinal mucosa in active celiac disease.

Killer cell immunoglobulin-like receptors (KIRs) modulate natural killer (NK) and T-cell function by human leukocyte antigen class I interaction and have been implicated in celiac disease (CD). Qualitative expression of 16 KIR genes was determined in biopsies from 22 CD patients at diagnosis and after >2 years on a gluten-free diet (GFD). Quantitative expression analysis of KIR2DL4, KIR3DL1, KIR3DL3, and KLRC2 (a marker of an NK-reprogrammed T-cell subpopulation augmented in CD) was performed in 35 additional CD biopsy pairs and 14 non-CD control biopsies. No specific KIR expression profile was observed in CD. KIR3DL1 was more frequently expressed in active CD compared with GFD (p = 0.0312) and controls (p = 0.0008), with slightly increased levels in active disease. KLRC2 was overexpressed in active (p = 0.0037) and GFD (p = 0.0469) patients compared with non-CD controls and coexpressed with KIR3DL1. Results suggest the participation of KIR3DL1 overexpression in the overall immune activation seen in CD mucosa, which could be partly explained by the NK-like T-cell subpopulation increase.

A regulatory single nucleotide polymorphism in the ubiquitin D gene associated with celiac disease.

An aberrant immune response triggered by dietary gluten is the main driving force underlying celiac disease (CD), but other biologic pathways that are dysregulated also participate in disease development. Genetic variation within these pathways might influence expression, contributing to susceptibility to CD. We have investigated the implication of ubiquitin D (UBD), a member of the ubiquitin-proteasome system that is strongly upregulated in the intestinal mucosa of active CD. Reverse transcriptase-polymerase chain reaction analysis of intestinal biopsy sample pairs (at diagnosis vs treated) from 30 CD patients confirmed overexpression of UBD in active disease tissue (fold change = 8.3; p = 0.0022). In silico prediction tools identified rs11724 as a putative regulatory single nucleotide polymorphism and association analysis of 468 CD patients and 459 controls revealed that the minor rs11724*C allele was more frequent among patients (minor allele frequency = 0.44 vs 0.39; odds ratio [OR] = 1.23; p = 0.028) and suggested a dominant allele effect (OR = 1.49; p = 0.0045). Correlation of the rs11724 genotype and UBD mRNA levels (OR = 0.76; p = 0.0021) further supports its implication in disease development.

Combined functional and positional gene information for the identification of susceptibility variants in celiac disease.

BACKGROUND & AIMS: Celiac disease is a complex, immune-mediated disorder of the intestinal mucosa with a strong genetic component. HLA-DQ2 is the major determinant of risk, but other minor genes, still to be identified, also are involved. METHODS: We designed a strategy that combines gene expression profiling of intestinal biopsy specimens, linkage region information, and different bioinformatics tools for the selection of potentially regulatory single-nucleotide polymorphisms (SNPs) involved in the disease. We selected 361 SNPs from 71 genes that fulfilled stringent functional (changes in expression level) and positional criteria (located in regions that have been linked to the disease, other than HLA). These polymorphisms were genotyped in 262 celiac patients and 214 controls. RESULTS: We detected strong evidence of association with several SNPs (the most significant were rs6747096, P = 2.38 x 10(-5); rs7040561, P = 6.55 x 10(-5); and rs458046, P = 1.35 x 10(-4)) that pinpoint novel candidate determinants of predisposition to the disease in previously identified linkage regions (eg, SERPINE2 in 2q33, and PBX3 or PPP6C in 9q34). CONCLUSIONS: Our study shows that the combination of function and position is a valid strategy for the genetic dissection of complex traits.

Epigenetic defects of GNAS in patients with pseudohypoparathyroidism and mild features of Albright's hereditary osteodystrophy.

CONTEXT: Several endocrine disorders that share resistance to PTH are grouped under the term pseudohypoparathyroidism (PHP). PHP type I, associated with blunted PTH-induced nephrogenous cAMP formation and phosphate excretion, is subdivided according to the presence or absence of additional endocrine abnormalities, Albright's hereditary osteodystrophy (AHO), and reduced Gsalpha activity caused by GNAS mutations. OBJECTIVE: We sought to identify the molecular defect in four unrelated patients who were thought to have PHP-Ia because of PTH and TSH resistance and mild AHO features. METHODS: Gsalpha activity and mutation analysis, and assessment of GNAS haplotype, methylation, and gene expression were performed for probands and family members. RESULTS: Two patients showed modest decreases in erythrocyte Gsalpha activity. Instead of Gsalpha point mutations, however, all four patients showed methylation defects of the GNAS locus, a feature previously described only for PHP-Ib. Furthermore, one patient with an isolated loss of GNAS exon A/B methylation had the 3-kb STX16 deletion frequently identified in PHP-Ib patients. In all but one of the remaining patients, haplotype analysis excluded large deletions or uniparental disomy as the cause of the observed methylation changes. CONCLUSIONS: Our investigations indicate that an overlap may exist between molecular and clinical features of PHP-Ia and PHP-Ib. No current mechanisms can explain the AHO-like features of our patients, some of which may not be linked to GNAS. Therefore, patients with hormone resistance and AHO-like features in whom coding Gsalpha mutations have been excluded should be evaluated for epigenetic alterations within GNAS.