Maria-Mercè Garcia-Barceló

Genome-wide association study identifies NRG1 as a susceptibility locus for Hirschsprung's disease

Hirschsprungs disease (HSCR), or aganglionic megacolon, is a congenital disorder characterized by the absence of enteric ganglia in variable portions of the distal intestine. RET is a well-established susceptibility locus, although existing evidence strongly suggests additional loci contributing to sporadic HSCR. To identify these additional genetic loci, we carried out a genome-wide association study using the Affymetrix 500K marker set. We successfully genotyped 293,836 SNPs in 181 Chinese subjects with sporadic HSCR and 346 ethnically matched control subjects. The SNPs most associated with HSCR were genotyped in an independent set of 190 HSCR and 510 control subjects. Aside from SNPs in RET, the strongest overall associations in plausible candidate genes were found for 2 SNPs located in intron 1 of the neuregulin1 gene (NRG1) on 8p12, with rs16879552 and rs7835688 yielding odds ratios of 1.68 [CI95%:(1.40, 2.00), P = 1.80 × 10−8] and 1.98 [CI95%:(1.59, 2.47), P = 1.12 × 10−9], respectively, for the heterozygous risk genotypes under an additive model. There was also a significant interaction between RET and NRG1 (P = 0.0095), increasing the odds ratio 2.3-fold to 19.53 for the RET rs2435357 risk genotype (TT) in the presence of the NRG1 rs7835688 heterozygote, indicating that NRG1 is a modifier of HSRC penetrance. Our highly significant association findings are backed-up by the important role of NRG1 as regulator of the development of the enteric ganglia precursors. The identification of NRG1 as an additional HSCR susceptibility locus not only opens unique fields of investigation into the mechanisms underlying the HSCR pathology, but also the mechanisms by which a discrete number of loci interact with each other to cause disease.

Identification of a HOXD13 mutation in a VACTERL patient

VACTERL acronym is assigned to a non‐random association of malformations in humans with poorly known etiology. It is comprised of vertebral defects (V), anal atresia (A), cardiac anomaly (C), tracheoesophageal fistula with esophageal atresia (TE), renal dysplasia (R) and limb lesions (L). Here, we report on, for the first time, a female patient with VACTERL association with a 21 base‐pair deletion in the exon 1 triplet repeats of HOXD13, a sonic hedgehog (SHH) downstream target. Our data provide the first piece of clinical evidence of the implication of the SHH pathway in VACTERL. Moreover, HOXD13 may not only be implicated in limb malformations but also in the development of gut and genitourinary structures, as predicted from the mouse models.

Genome-wide association study identifies a susceptibility locus for biliary atresia on 10q24.2

Biliary atresia (BA) is characterized by the progressive fibrosclerosing obliteration of the extrahepatic biliary system during the first few weeks of life. Despite early diagnosis and prompt surgical intervention, the disease progresses to cirrhosis in many patients. The current theory for the pathogenesis of BA proposes that during the perinatal period, a still unknown exogenous factor meets the innate immune system of a genetically predisposed individual and induces an uncontrollable and potentially self-limiting immune response, which becomes manifest in liver fibrosis and atresia of the extrahepatic bile ducts. Genetic factors that could account for the disease, let alone for its high incidence in Chinese, are to be investigated. To identify BA susceptibility loci, we carried out a genome-wide association study (GWAS) using the Affymetrix 5.0 and 500 K marker sets. We genotyped nearly 500 000 single-nucleotide polymorphisms (SNPs) in 200 Chinese BA patients and 481 ethnically matched control subjects. The 10 most BA-associated SNPs from the GWAS were genotyped in an independent set of 124 BA and 90 control subjects. The strongest overall association was found for rs17095355 on 10q24, downstream XPNPEP1, a gene involved in the metabolism of inflammatory mediators. Allelic chi-square test P-value for the meta-analysis of the GWAS and replication results was 6.94 x 10(-9). The identification of putative BA susceptibility loci not only opens new fields of investigation into the mechanisms underlying BA but may also provide new clues for the development of preventive and curative strategies.

A Germline Mutation (A339V) in Thyroid Transcription Factor-1 (TITF-1/NKX2.1) in Patients With Multinodular Goiter and Papillary Thyroid Carcinoma

BACKGROUND The genetic factors that determine the risk of papillary thyroid carcinoma (PTC) among patients with multinodular goiter (MNG) remain undefined. Because thyroid transcription factor-1 (TTF-1) is important to thyroid development, we evaluated whether the gene that encodes it, TITF-1/NKX2.1, is a genetic determinant of MNG/PTC predisposition. METHODS Twenty unrelated PTC patients with a history of MNG (MNG/PTC), 284 PTC patients without a history of MNG (PTC), and 349 healthy control subjects were screened for germline mutation(s) in TITF-1/NKX2.1 by sequencing of amplified DNA from blood. The effects of the mutation on the growth and differentiation of thyroid cells were demonstrated by ectopic expression of wild-type (WT) and mutant proteins in PCCL3 normal rat thyroid cells, followed by tests of cell proliferation, activation of cell growth pathways, and transcription of TTF-1 target genes. All statistical tests were two-sided. RESULTS A missense mutation (1016C>T) was identified in TITF-1/NKX2.1 that led to a mutant TTF-1 protein (A339V) in four of the 20 MNG/PTC patients (20%). These patients developed substantially more advanced tumors than MNG/PTC or PTC patients without the mutation (P = .022, Fisher exact test). Notably, this germline mutation was dominantly inherited in two families, with some members bearing the mutation affected with MNG, associated with either PTC or colon cancer. The mutation encoding the A339V substitution was not found among the 349 healthy control subjects nor among the 284 PTC patients who had no history of MNG. Overexpression of A339V TTF-1 in PCCL3 cells, as compared with overexpression of WT TTF-1, was associated with increased cell proliferation including thyrotropin-independent growth (average A339V proliferation rate = 134.27%, WT rate = 104.43%, difference = 34.3%, 95% confidence interval = 12.0% to 47.7%, P = .010), enhanced STAT3 activation, and impaired transcription of the thyroid-specific genes Tg, TSH-R, and Pax-8. CONCLUSION This is the first germline mutation identified in MNG/PTC patients. It could contribute to predisposition for MNG and/or PTC and to the pathogenesis of PTC.

ELF1 is associated with systemic lupus erythematosus in Asian populations

Systemic lupus erythematosus (SLE) is an autoimmune disease with a strong genetic involvement. The susceptibility genes identified so far can only explain a small proportion of disease heritability. Through a genome-wide association in a Hong Kong Chinese cohort and subsequent replication in two other Asian populations, with a total of 3164 patients and 4482 matched controls, we identified association of ELF1 (E74-like factor 1) with SLE (rs7329174, OR = 1.26, joint P= 1.47 × 10(-8)). ELF1 belongs to the ETS family of transcription factors and is known to be involved in T cell development and function. Database analysis revealed transcripts making use of three alternative exon1s for this gene. Near equivalent expression levels of distinct transcripts initiated from alternative exon1s were detected in peripheral blood mononuclear cells from both SLE patients and healthy controls. Although a direct association of rs7329174 with the three forms of transcripts for this gene was not detected, these findings support an important role of ELF1 in SLE susceptibility and suggest a potentially tight regulation for the expression of this gene.

Hedgehog/Notch-induced premature gliogenesis represents a new disease mechanism for Hirschsprung disease in mice and humans

Hirschsprung (HSCR) disease is a complex genetic disorder attributed to a failure of the enteric neural crest cells (ENCCs) to form ganglia in the hindgut. Hedgehog and Notch are implicated in mediating proliferation and differentiation of ENCCs. Nevertheless, how these signaling molecules may interact to mediate gut colonization by ENCCs and contribute to a primary etiology for HSCR are not known. Here, we report our pathway-based epistasis analysis of data generated by a genome-wide association study on HSCR disease, which indicates that specific genotype constellations of Patched (PTCH1) (which encodes a receptor for Hedgehog) and delta-like 3 (DLL3) (which encodes a receptor for Notch) SNPs confer higher risk to HSCR. Importantly, deletion of Ptch1 in mouse ENCCs induced robust Dll1 expression and activation of the Notch pathway, leading to premature gliogenesis and reduction of ENCC progenitors in mutant bowels. Dll1 integrated Hedgehog and Notch pathways to coordinate neuronal and glial cell differentiation during enteric nervous system development. In addition, Hedgehog-mediated gliogenesis was found to be highly conserved, such that Hedgehog was consistently able to promote gliogenesis of human neural crest-related precursors. Collectively, we defined PTCH1 and DLL3 as HSCR susceptibility genes and suggest that Hedgehog/Notch-induced premature gliogenesis may represent a new disease mechanism for HSCR.

Replication study of SNP associations for colorectal cancer in Hong Kong Chinese

Background:Recent genome-wide association studies of colorectal cancer (CRC) have identified common single-nucleotide polymorphisms (SNPs) mapping to 10 independent loci that confer modest increased risk. These studies have been conducted in European populations and it is unclear whether these observations generalise to populations with different ethnicities and rates of CRC.Methods:An association study was performed on 892 CRC cases and 890 controls recruited from the Hong Kong Chinese population, genotyping 32 SNPs, which were either associated with CRC in previous studies or are in close proximity to previously reported risk SNPs.Results:Twelve of the SNPs showed evidence of an association. The strongest associations were provided by rs10795668 on 10p14, rs4779584 on 15q14 and rs12953717 on 18q21.2. There was significant linear association between CRC risk and the number of independent risk variants possessed by an individual (P=2.29 × 10−5).Conclusion:These results indicate that some previously reported SNP associations also impact on CRC risk in the Chinese population. Possible reasons for failure of replication for some loci include inadequate study power, differences in allele frequency, linkage disequilibrium structure or effect size between populations. Our results suggest that many associations for CRC are likely to generalise across populations.

Common Variants on Xq28 Conferring Risk of Schizophrenia in Han Chinese

Schizophrenia is a highly heritable, severe psychiatric disorder affecting approximately 1% of the world population. A substantial portion of heritability is still unexplained and the pathophysiology of schizophrenia remains to be elucidated. To identify more schizophrenia susceptibility loci, we performed a genome-wide association study (GWAS) on 498 patients with schizophrenia and 2025 controls from the Han Chinese population, and a follow-up study on 1027 cases and 1005 controls. In the follow-up study, we included 384 single nucleotide polymorphisms (SNPs) which were selected from the top hits in our GWAS (130 SNPs) and from previously implicated loci for schizophrenia based on the SZGene database, NHGRI GWAS Catalog, copy number variation studies, GWAS meta-analysis results from the international Psychiatric Genomics Consortium (PGC) and candidate genes from plausible biological pathways (254 SNPs). Within the chromosomal region Xq28, SNP rs2269372 in RENBP achieved genome-wide significance with a combined P value of 3.98 × 10(-8) (OR of allele A = 1.31). SNPs with suggestive P values were identified within 2 genes that have been previously implicated in schizophrenia, MECP2 (rs2734647, P combined = 8.78 × 10(-7), OR = 1.28; rs2239464, P combined = 6.71 × 10(-6), OR = 1.26) and ARHGAP4 (rs2269368, P combined = 4.74 × 10(-7), OR = 1.25). In addition, the patient sample in our follow-up study showed a significantly greater burden for pre-defined risk alleles based on the SNPs selected than the controls. This indicates the existence of schizophrenia susceptibility loci among the SNPs we selected. This also further supports multigenic inheritance in schizophrenia. Our findings identified a new schizophrenia susceptibility locus on Xq28, which harbor the genes RENBP, MECP2, and ARHGAP4.

Genome-wide copy number analysis uncovers a new HSCR gene: NRG3

Hirschsprung disease (HSCR) is a congenital disorder characterized by aganglionosis of the distal intestine. To assess the contribution of copy number variants (CNVs) to HSCR, we analysed the data generated from our previous genome-wide association study on HSCR patients, whereby we identified NRG1 as a new HSCR susceptibility locus. Analysis of 129 Chinese patients and 331 ethnically matched controls showed that HSCR patients have a greater burden of rare CNVs (p = 1.50×10−5), particularly for those encompassing genes (p = 5.00×10−6). Our study identified 246 rare-genic CNVs exclusive to patients. Among those, we detected a NRG3 deletion (p = 1.64×10−3). Subsequent follow-up (96 additional patients and 220 controls) on NRG3 revealed 9 deletions (combined p = 3.36×10−5) and 2 de novo duplications among patients and two deletions among controls. Importantly, NRG3 is a paralog of NRG1. Stratification of patients by presence/absence of HSCR–associated syndromes showed that while syndromic–HSCR patients carried significantly longer CNVs than the non-syndromic or controls (p = 1.50×10−5), non-syndromic patients were enriched in CNV number when compared to controls (p = 4.00×10−6) or the syndromic counterpart. Our results suggest a role for NRG3 in HSCR etiology and provide insights into the relative contribution of structural variants in both syndromic and non-syndromic HSCR. This would be the first genome-wide catalog of copy number variants identified in HSCR.

Fine Mapping of the NRG1 Hirschsprung's Disease Locus

The primary pathology of Hirschsprungs disease (HSCR, colon aganglionosis) is the absence of ganglia in variable lengths of the hindgut, resulting in functional obstruction. HSCR is attributed to a failure of migration of the enteric ganglion precursors along the developing gut. RET is a key regulator of the development of the enteric nervous system (ENS) and the major HSCR-causing gene. Yet the reduced penetrance of RET DNA HSCR-associated variants together with the phenotypic variability suggest the involvement of additional genes in the disease. Through a genome-wide association study, we uncovered a ∼350 kb HSCR-associated region encompassing part of the neuregulin-1 gene (NRG1). To identify the causal NRG1 variants contributing to HSCR, we genotyped 243 SNPs variants on 343 ethnic Chinese HSCR patients and 359 controls. Genotype analysis coupled with imputation narrowed down the HSCR-associated region to 21 kb, with four of the most associated SNPs (rs10088313, rs10094655, rs4624987, and rs3884552) mapping to the NRG1 promoter. We investigated whether there was correlation between the genotype at the rs10088313 locus and the amount of NRG1 expressed in human gut tissues (40 patients and 21 controls) and found differences in expression as a function of genotype. We also found significant differences in NRG1 expression levels between diseased and control individuals bearing the same rs10088313 risk genotype. This indicates that the effects of NRG1 common variants are likely to depend on other alleles or epigenetic factors present in the patients and would account for the variability in the genetic predisposition to HSCR.