Identification of a wide spectrum of ciliary gene mutations in nonsyndromic biliary atresia patients implicates ciliary dysfunction as a novel disease mechanism

Background

Biliary atresia (BA) is the most common obstructive cholangiopathy in neonates, often progressing to end-stage cirrhosis. BA pathogenesis is believed to be multifactorial, but the genetic contribution, especially for nonsyndromic BA (common form: > 85%) remains poorly defined.

Methods

We conducted whole exome sequencing on 89 nonsyndromic BA trios to identify rare variants contributing to BA etiology. Functional evaluation using patients’ liver biopsies, human cell and zebrafish models were performed. Clinical impact on respiratory system was assessed with clinical evaluation, nasal nitric oxide (nNO), high speed video analysis and transmission electron microscopy.

Findings

We detected rare, deleterious de novo or biallelic variants in liver-expressed ciliary genes in 31.5% (28/89) of the BA patients. Burden test revealed 2.6-fold (odds ratio (OR) [95% confidence intervals (CI)]= 2.58 [1.15–6.07], adjusted p = 0.034) over-representation of rare, deleterious mutations in liver-expressed ciliary gene set in patients compared to controls. Functional analyses further demonstrated absence of cilia in the BA livers with KIF3B and TTC17 mutations, and knockdown of PCNT, KIF3B and TTC17 in human control fibroblasts and cholangiocytes resulted in reduced number of cilia. Additionally, CRISPR/Cas9-engineered zebrafish knockouts of KIF3B, PCNT and TTC17 displayed reduced biliary flow. Abnormally low level of nNO was detected in 80% (8/10) of BA patients carrying deleterious ciliary mutations, implicating the intrinsic ciliary defects.

Interpretation

Our findings support strong genetic susceptibility for nonsyndromic BA. Ciliary gene mutations leading to cholangiocyte cilia malformation and dysfunction could be a key biological mechanism in BA pathogenesis.

Funding

The study is supported by General Research Fund, HMRF Commissioned Paediatric Research at HKCH and Li Ka Shing Faculty of Medicine Enhanced New Staff Start-up Fund.

Exome chip association study excluded the involvement of rare coding variants with large effect sizes in the etiology of anorectal malformations

INTRODUCTION

Anorectal malformations (ARM) are rare congenital malformations, resulting from disturbed hindgut development. A genetic etiology has been suggested, but evidence for the involvement of specific genes is scarce. We evaluated the contribution of rare and low-frequency coding variants in ARM etiology, assuming a multifactorial model.

METHODS

We analyzed 568 Caucasian ARM patients and 1,860 population-based controls using the Illumina HumanExome Beadchip array, which contains >240,000 rare and low-frequency coding variants. GenomeStudio clustering and calling was followed by re-calling of 'no-calls' using zCall for patients and controls simultaneously. Single variant and gene-based analyses were performed to identify statistically significant associations, applying Bonferroni correction. Following an extra quality control step, candidate variants were selected for validation using Sanger sequencing.

RESULTS

When we applied a MAF of ≥1.0%, no variants or genes showed statistically significant associations with ARM. Using a MAF cut-off at 0.4%, 13 variants initially reached statistical significance, but had to be discarded upon further inspection: ten variants represented calling errors of the software, while the minor alleles of the remaining three variants were not confirmed by Sanger sequencing.

CONCLUSION

Our results show that rare and low-frequency coding variants with large effect sizes, present on the exome chip do not contribute to ARM etiology.

Conditional deletion of platelet derived growth factor receptor alpha (Pdgfra) in urorectal mesenchyme causes mesenchyme apoptosis and urorectal developmental anomalies in mice

In mammals, urorectal development starts at early embryonic stage, defective urorectal development results in anorectal malformations, which are common congenital developmental defects of the anus and the urethra in newborns. The etiology and embryology of the defects are still largely unknown. Platelet-derived growth factor receptor alpha (Pdgfra) is a cell surface receptor tyrosine kinase, upon binding to its ligands (Pdgfa-d), mediates intracellular signaling and regulates embryonic development. The expression of Pdgfra is tightly regulated in the developing urorectal mesenchyme, and its dysregulation is associated with urorectal defects in animals with urorectal defects. Knockout of Pdgfra induces early embryo lethality which precludes investigation of Pdgfra in urorectal development. To address the temporal requirement of Pdgfra in urorectal development, we conditionally deleted Pdgfra in Pdgfra-expressing tissues using a tamoxifen inducible Cre-loxP approach in mice, examined the urorectal development in Pdgfra conditional knockout (Pdgfra-cKO) embryos. We showed that conditional deletion of Pdgfra in Pdgfra-expressing tissues at E10-E11 caused cloaca septation defect, anteriorly displaced anus, defective urogenital folds development and abnormal urethra tubularization in both male and female mice. Furthermore, we showed that Pdgfra was required for the survival of urorectal mesenchyme, deletion of Pdgfra caused apoptosis in the peri-cloacal, the peri-urethra and the urorectal septum mesenchyme of Pdgfra-cKO mutants, associated with an induction of p53, Ndrg1 and activation of caspase-3 in Pdgfra-cKO embryos. In conclusion, Pdgfra is required for the development and survival of the urorectal mesenchyme in embryo, dysregulated Pdgfra signaling induced urorectal defects in mice resembling human congenital diseases of anorectal malformations and hypospadias. Perturbation of PDGFRA signaling may contribute to anorectal malformations and hypospadias in human.

Trans-ethnic meta-analysis of genome-wide association studies for Hirschsprung disease

Hirschsprung disease (HSCR) is the most common cause of neonatal intestinal obstruction. It is characterized by the absence of ganglia in the nerve plexuses of the lower gastrointestinal tract. So far, three common disease-susceptibility variants at the RET, SEMA3 and NRG1 loci have been detected through genome-wide association studies (GWAS) in Europeans and Asians to understand its genetic etiologies. Here we present a trans-ethnic meta-analysis of 507 HSCR cases and 1191 controls, combining all published GWAS results on HSCR to fine-map these loci and narrow down the putatively causal variants to 99% credible sets. We also demonstrate that the effects of RET and NRG1 are universal across European and Asian ancestries. In contrast, we detected a European-specific association of a low-frequency variant, rs80227144, in SEMA3 [odds ratio (OR) = 5.2, P = 4.7 × 10-10]. Conditional analyses on the lead SNPs revealed a secondary association signal, corresponding to an Asian-specific, low-frequency missense variant encoding RET p.Asp489Asn (rs9282834, conditional OR = 20.3, conditional P = 4.1 × 10-14). When in trans with the RET intron 1 enhancer risk allele, rs9282834 increases the risk of HSCR from 1.1 to 26.7. Overall, our study provides further insights into the genetic architecture of HSCR and has profound implications for future study designs.

Patient complexity and genotype-phenotype correlations in biliary atresia: a cross-sectional analysis

Background

Biliary Atresia (BA) is rare and genetically complex, and the pathogenesis is elusive. The disease course is variable and can represent heterogeneity, which hinders effective disease management. Deciphering the BA phenotypic variance is a priority in clinics and can be achieved by the integrative analysis of genotype and phenotype. We aim to explore the BA phenotypic features and to delineate the source of its variance.

Methods

The study is a cross-sectional observational study collating with case/control association analysis. One-hundred-and-eighty-one type III non-syndromic BA patients and 431 controls were included for case–control association tests, including 89 patients (47.19% males, born June 15th, 1981 to September 17th, 2007) have detailed clinical records with follow-up of the disease course (median ~17.2 years). BA-association genes from the genome-wide gene-based association test on common genetic variants (CV) and rare copy-number-variants (CNVs) from the genome-wide survey, the later comprise only CNVs > 100 kb and found in the BA patients but not in the local population (N = 1,381) or the database (N = 11,943). Hereby comorbidity is defined as a chronic disease that affects the BA patients but has no known relationship with BA or with the BA treatment. We examined genotype-phenotype correlations of CNVs, connectivity of these novel variants with BA-associated CVs, and their role in the BA candidate gene network.

Results

Of the 89 patients, 41.57% have comorbidities, including autoimmune-allergic disorders (22.47%). They carried 29 BA-private CNVs, including 3 CNVs underpinning the carriers’ immunity comorbidity and one JAG1 micro-deletion. The BA-CNV-intersected genes (N = 102) and the CV-tagged genes (N = 103) were both enriched with immune-inflammatory pathway genes (FDR q < 0.20), and the two gene sets were interconnected (permutation p = 0.039). The molecular network representing CVs and rare-CNV association genes fit into a core/periphery structure, the immune genes and their related modules are found at the coherence core of all connections, suggesting its dominant role in the BA pathogenesis pathway.

Conclusions

The study highlights a patient-complexity phenomenon as a novel BA phenotypic feature, which is underpinned by rare-CNVs that biologically converge with CVs into the immune-inflammatory pathway and drives the BA occurrence and the likely BA association with immune diseases in clinics.

Electronic supplementary material

The online version of this article (doi:10.1186/s12920-017-0259-0) contains supplementary material, which is available to authorized users.

Whole exome sequencing coupled with unbiased functional analysis reveals new Hirschsprung disease genes

Background

Hirschsprung disease (HSCR), which is congenital obstruction of the bowel, results from a failure of enteric nervous system (ENS) progenitors to migrate, proliferate, differentiate, or survive within the distal intestine. Previous studies that have searched for genes underlying HSCR have focused on ENS-related pathways and genes not fitting the current knowledge have thus often been ignored. We identify and validate novel HSCR genes using whole exome sequencing (WES), burden tests, in silico prediction, unbiased in vivo analyses of the mutated genes in zebrafish, and expression analyses in zebrafish, mouse, and human.

Results

We performed de novo mutation (DNM) screening on 24 HSCR trios. We identify 28 DNMs in 21 different genes. Eight of the DNMs we identified occur in RET, the main HSCR gene, and the remaining 20 DNMs reside in genes not reported in the ENS. Knockdown of all 12 genes with missense or loss-of-function DNMs showed that the orthologs of four genes (DENND3, NCLN, NUP98, and TBATA) are indispensable for ENS development in zebrafish, and these results were confirmed by CRISPR knockout. These genes are also expressed in human and mouse gut and/or ENS progenitors. Importantly, the encoded proteins are linked to neuronal processes shared by the central nervous system and the ENS.

Conclusions

Our data open new fields of investigation into HSCR pathology and provide novel insights into the development of the ENS. Moreover, the study demonstrates that functional analyses of genes carrying DNMs are warranted to delineate the full genetic architecture of rare complex diseases.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1174-6) contains supplementary material, which is available to authorized users.

HOXB5 binds to multi-species conserved sequence (MCS+9.7) of RET gene and regulates RET expression

RET gene is crucial for the development of enteric nervous system, and dys-regulation of RET expression causes Hirschsprung disease. HOXB5 regulates RET transcription, and perturbations in transcriptional regulation by HOXB5 caused reduced RET expression and defective enteric nervous system development in mice. The mechanisms by which HOXB5 regulate RET transcription are unclear. Thus, unraveling the regulatory mechanisms of HOXB5 on RET transcription could lead to a better understanding of the etiology of Hirschsprung disease. In this study, we identified and confirmed HOXB5 binding to the multi-species conserved sequence (MCS+9.7) in the first intron of the RET gene. We developed a RET mini-gene reporter system, and showed that MCS+9.7 enhanced HOXB5 trans-activation from RET promoter in human neuroblastoma SK-N-SH cells and in chick embryos. The deletion of HOXB5 binding site interfered with HOXB5 trans-activation. Furthermore, transfection of HOXB5 induced endogenous RET transcription, enhanced the co-precipitation of TATA-box binding protein with the transcription start site of RET, and induced histone H3K4 trimethylation in chromatin regions upstream and downstream of RET transcription start site. In conclusion, (i) HOXB5 physically interacted with MCS+9.7 and enhanced RET transcription, (ii) HOXB5 altered chromatin conformation and histone modification of RET locus, which could facilitate the formation of transcription complex, and enhance RET transcription, (iii) expression of RET was mediated by a complex regulatory network of transcription factors functioning in a synergistic, additive and/or independent manners. Hence, dys-regulation of RET expression by HOXB5 could result in insufficient RET expression and Hirschsprung disease.

Common genetic variants regulating ADD3 gene expression alter biliary atresia risk

BACKGROUND & AIMS

Biliary atresia (BA) is a rare and most severe cholestatic disease in neonates, but the pathogenic mechanisms are unknown. Through a previous genome wide association study (GWAS) on Han Chinese, we discovered association of the 10q24.2 region encompassing ADD3 and XPNPEP1 genes, which was replicated in Chinese and Thai populations. This study aims to fully characterize the genetic architecture at 10q24.2 and to reveal the link between the genetic variants and BA.

METHODS

We genotyped 107 single nucleotide polymorphisms (SNPs) in 10q24.2 in 339 Han Chinese patients and 401 matched controls using Sequenom. Exhaustive follow-up studies of the association signals were performed.

RESULTS

The combined BA-association p-value of the GWAS SNP (rs17095355) achieved 6.06×10(-10). Further, we revealed the common risk haplotype encompassing 5 tagging-SNPs, capturing the risk-predisposing alleles in 10q24.2 [p=5.32×10(-11); odds ratio, OR: 2.38; confidence interval, CI: (2.14-2.62)]. Through Sanger sequencing, no deleterious rare variants (RVs) residing in the risk haplotype were found, dismissing the theory of "synthetic" association. Moreover, in bioinformatics and in vivo genotype-expression investigations, the BA-associated potentially regulatory SNPs correlated with ADD3 gene expression (n=36; p=0.0030). Remarkably, the risk haplotype frequency coincides with BA incidences in the population, and, positive selection (favoring the derived alleles that arose from mutations) was evident at the ADD3 locus, suggesting a possible role for the BA-associated common variants in shaping the general population diversity.

CONCLUSIONS

Common genetic variants in 10q24.2 can alter BA risk by regulating ADD3 expression levels in the liver, and may exert an effect on disease epidemiology and on the general population.

Three SNPs in chromosome 11q23.3 are independently associated with systemic lupus erythematosus in Asians

Systemic lupus erythematosus (SLE) has a complex etiology and is affected by both genetic and environmental factors. Although more than 40 loci have shown robust association with SLE, the details of these loci, such as the independent contributors and the genes involved, are still unclear. In this study, we performed meta-analysis of two existing genome-wide association studies (GWASs) on Chinese Han populations from Hong Kong and Anhui, China, and followed the findings by further replication on three additional Chinese and Thailand cohorts with a total of 4254 cases and 6262 controls matched geographically and ethnically. We discovered multiple susceptibility variants for SLE in the 11q23.3 region, including variants in/near PHLDB1 (rs11603023, P(_combined) = 1.25E-08, OR = 1.20), DDX6 (rs638893, P(_combined) = 5.19E-07, OR = 1.22) and CXCR5 (rs10892301, P(_combined) = 2.51E-08, OR = 0.85). Genetic contributions from the newly identified variants were all independent of SNP rs4639966, whose association was reported from the previous GWAS. In addition, the three newly identified variants all showed independent association with the disease through modeling by both stepwise and conditional logistic regression. The presence of multiple independent variants in this region emphasizes its role in SLE susceptibility, and also hints the possibility that distinct biological mechanisms might be involved in the disease involving this genomic region.

Genetic Analyses of a Three Generation Family Segregating Hirschsprung Disease and Iris Heterochromia

We present the genetic analyses conducted on a three-generation family (14 individuals) with three members affected with isolated-Hirschsprung disease (HSCR) and one with HSCR and heterochromia iridum (syndromic-HSCR), a phenotype reminiscent of Waardenburg-Shah syndrome (WS4). WS4 is characterized by pigmentary abnormalities of the skin, eyes and/or hair, sensorineural deafness and HSCR. None of the members had sensorineural deafness. The family was screened for copy number variations (CNVs) using Illumina-HumanOmni2.5-Beadchip and for coding sequence mutations in WS4 genes (EDN3, EDNRB, or SOX10) and in the main HSCR gene (RET). Confocal microscopy and immunoblotting were used to assess the functional impact of the mutations. A heterozygous A/G transition in EDNRB was identified in 4 affected and 3 unaffected individuals. While in EDNRB isoforms 1 and 2 (cellular receptor) the transition results in the abolishment of translation initiation (M1V), in isoform 3 (only in the cytosol) the replacement occurs at Met91 (M91V) and is predicted benign. Another heterozygous transition (c.-248G/A; -predicted to affect translation efficiency-) in the 5'-untranslated region of EDN3 (EDNRB ligand) was detected in all affected individuals but not in healthy carriers of the EDNRB mutation. Also, a de novo CNVs encompassing DACH1 was identified in the patient with heterochromia iridum and HSCR Since the EDNRB and EDN3 variants only coexist in affected individuals, HSCR could be due to the joint effect of mutations in genes of the same pathway. Iris heterochromia could be due to an independent genetic event and would account for the additional phenotype within the family.

Meta-analysis followed by replication identifies loci in or near CDKN1B, TET3, CD80, DRAM1, and ARID5B as associated with systemic lupus erythematosus in Asians

Systemic lupus erythematosus (SLE) is a prototype autoimmune disease with a strong genetic involvement and ethnic differences. Susceptibility genes identified so far only explain a small portion of the genetic heritability of SLE, suggesting that many more loci are yet to be uncovered for this disease. In this study, we performed a meta-analysis of genome-wide association studies on SLE in Chinese Han populations and followed up the findings by replication in four additional Asian cohorts with a total of 5,365 cases and 10,054 corresponding controls. We identified genetic variants in or near CDKN1B, TET3, CD80, DRAM1, and ARID5B as associated with the disease. These findings point to potential roles of cell-cycle regulation, autophagy, and DNA demethylation in SLE pathogenesis. For the region involving TET3 and that involving CDKN1B, multiple independent SNPs were identified, highlighting a phenomenon that might partially explain the missing heritability of complex diseases.

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⁻⁵), particularly for those encompassing genes (p = 5.00×10⁻⁶). Our study identified 246 rare-genic CNVs exclusive to patients. Among those, we detected a NRG3 deletion (p = 1.64×10⁻³). Subsequent follow-up (96 additional patients and 220 controls) on NRG3 revealed 9 deletions (combined p = 3.36×10⁻⁵) 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⁻⁵), non-syndromic patients were enriched in CNV number when compared to controls (p = 4.00×10⁻⁶) 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.

Copy number variations (CNVs) are significant genetic risk factors in disease pathogenesis and represent an important portion of missing heritability for some human diseases, making their discovery essential for the identification of genes and risk factors for a wide range of diseases, including Hirschsprung disease (HSCR, congenital colon aganglionosis). Since the discovery of the major HSCR gene, RET, a number of rare mutations have been reported in RET and other genes involved in the development of the enteric nervous system. However, these mutations contribute to only a small proportion of the disease susceptibility. Taking advantage of the recent technical and methodological advances, we have examined the contribution of CNVs to the disease. We have found that HSCR patients are enriched with CNVs encompassing genes. In particular, we found that deletions of NRG3, a paralog of the previously identified HSCR–susceptibility gene NRG1, were associated with the HSCR phenotype.

Functional analyses of RET mutations in Chinese Hirschsprung disease patients

BACKGROUND

Hirschsprung disease (HSCR) is a congenital disease characterized by the absence of ganglion cells in various length of distal digestive tract. The rearranged during transfection gene (RET) is considered the major gene in HSCR. Although an increasing number of HSCR-associated RET coding sequence (CDS) mutations have been identified in recent years, not many have been investigated for functional consequence on the RET protein.

METHODS AND RESULTS

We examined the functional implications of the de novo RET-CDS mutations V145G, Y483X, V636fsX1, and F961L that we first identified in sporadic Chinese patients with HSCR. The V145G disrupted RET glycosylation and F961L RET phosphorylation. Presumably, the truncation mutations would affect the translocation or the anchoring of the RET protein onto the cellular membrane.

CONCLUSION

The study of RET-CDS mutations that appear de novo is essential not only for understanding the mechanistic of the disease but also for penetrance and recurrence risk estimations, being the ultimate goal for the improvement in disease management and counseling.

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.

Fine mapping of the NRG1 Hirschsprung's disease locus

The primary pathology of Hirschsprung's 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.

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.

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.

Fine mapping of the 9q31 Hirschsprung's disease locus

Hirschsprung’s disease (HSCR) is a congenital disorder characterised by the absence of ganglia along variable lengths of the intestine. The RET gene is the major HSCR gene. Reduced penetrance of RET mutations and phenotypic variability suggest the involvement of additional modifying genes in the disease. A RET-dependent modifier locus was mapped to 9q31 in families bearing no coding sequence (CDS) RET mutations. Yet, the 9q31 causative locus is to be identified. To fine-map the 9q31 region, we genotyped 301 tag-SNPs spanning 7 Mb on 137 HSCR Dutch trios. This revealed two HSCR-associated regions that were further investigated in 173 Chinese HSCR patients and 436 controls using the genotype data obtained from a genome-wide association study recently conducted. Within one of the two identified regions SVEP1 SNPs were found associated with Dutch HSCR patients in the absence of RET mutations. This ratifies the reported linkage to the 9q31 region in HSCR families with no RET CDS mutations. However, this finding could not be replicated. In Chinese, HSCR was found associated with IKBKAP. In contrast, this association was stronger in patients carrying RET CDS mutations with p = 5.10 × 10⁻⁶ [OR = 3.32 (1.99, 5.59)] after replication. The HSCR-association found for IKBKAP in Chinese suggests population specificity and implies that RET mutation carriers may have an additional risk. Our finding is supported by the role of IKBKAP in the development of the nervous system.

Lack of association between nNOS -84G>A polymorphism and risk of infantile hypertrophic pyloric stenosis in a Chinese population

BACKGROUND

Infantile hypertrophic pyloric stenosis (IHPS) is one of the most common gastrointestinal obstructions in the infancy requiring surgery. Reduced expression of neuronal nitric oxide synthase (nNOS), which plays an important role in the regulation of the human pyloric muscle, is thought to underlie IHPS. The role of nNOS in IHPS has been supported by the genetic association of a functional regulatory nNOS polymorphism (-84G>A) with IHPS in whites. We reasoned that the corroboration of this association in a population of different ethnic origin would prompt follow-up studies and further investigation of the IHPS pathology at molecular level. Thus, we attempted to reproduce the original findings in a Chinese population of comparable size in what would be the first genetic study on IHPS conducted in Chinese.

METHODS

nNOS -84G>A genotypes were analyzed in 56 patients and 86 controls by polymerase chain reaction and DNA sequencing. Logistic regression was used to compute odds ratios.

RESULTS

Our study could not corroborate the association previously reported. Although the frequency of the IHPS-associated allele (-84A) in controls (0.205) was similar to that reported for white controls, there was a dramatic difference in -84A frequencies between white and Chinese patients (0.198). Similarly, there was no difference in the nNOS -84G>A genotype distribution between patients and controls, even when the GA and AA genotypes were combined to compare GG genotype (odds ratio, 1.01; 95% confidence interval, 0.47-2.19).

CONCLUSIONS

Failure to replicate the initial finding does not detract from its validity, because genetic effects may differ across populations. Differences across populations in linkage disequilibrium and/or allele frequencies may contribute to this lack of replication. The role nNOS in IHPS awaits further investigation.

MNX1 (HLXB9) mutations in Currarino patients

PURPOSE

The combination of partial absence of the sacrum, anorectal anomalies, and presacral mass constitutes Currarino syndrome (CS), which is associated with mutations in MNX1 motor neuron and pancreas homeobox 1 (previously HLXB9). Here, we report on the MNX1 mutations found in a family segregating CS and in 3 sporadic CS patients, as well as on the clinical characteristics of the affected individuals.

METHODS

MNX1 mutations were identified by direct sequencing the coding regions, intron/exon boundaries of MNX1 in 5 CS Japanese family members and 3 Chinese sporadic cases and their parents.

RESULTS

There were 2 novel (P18PfsX37, R243W) and 2 previously described (W288G and IVS2 + 1G > A) mutations. These mutations were not found in 198 control individuals and are predicted to impair the functioning of the MNX1 protein.

CONCLUSIONS

The variability of the CS phenotype among related or unrelated patients bearing the same mutation advocates for differences in the genetic background of each individual and invokes the implication of additional CS susceptibility genes.

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.

Evaluation of the thyroid transcription factor-1 gene (TITF1) as a Hirschsprung's disease locus

Hirschsprung's disease (HSCR, colonic aganglionosis) is an oligogenic entity that usually requires mutations in RET and other interacting loci. Decreased levels of RET expression may lead to the manifestation of HSCR. We previously showed that RET transcription was decreased due to alteration of the TITF1 binding site by two HSCR-associated RET promoter single nucleotide polymorphisms (SNPs). This prompted us to investigate whether DNA alterations in TITF1 could play a role in HSCR by affecting the RET-regulatory properties of the TITF1 protein. Our initial study on 86 Chinese HSCR patients revealed a Gly322Ser amino acid substitution in the TITF1protein. In this study we have examined an additional 102 Chinese and 70 Caucasian patients, and 194 Chinese and 60 Caucasian unselected, unrelated, subjects as controls. The relevance of the DNA changes detected in TITF1 by direct sequencing were evaluated using bioinformatics, reporter and binding-assays, mouse neurosphere culture, immunohistochemistry and immunofluorescence techniques. Met3Leu and Pro48Pro were identified in 2 Caucasian patients and 1 Chinese patient, respectively. In vitro analysis showed that Met3Leu reduced the activity of the RET promoter by 100% in the presence of the wild-type or HSCR-associated RET promoter SNP alleles. The apparent binding affinity of the TITF1 mutated protein was not decreased. The Met3Leu mutation may affect the interaction of TITF1 with its protein partners. The absence of Titf1 expression in mouse gut but not in human gut suggests that the role of TITF1 in gut development differs between the two species. TITF1 mutations could contribute to HSCR by affecting RET expression through defective interactions with other transcription factors.

Correlation between genetic variations in Hox clusters and Hirschsprung's disease

Interactions between migrating neural crest cells and the environment of the gut are crucial for the development of the enteric nervous system (ENS). A key signalling mediator is the RET-receptor-tyrosine-kinase which, when defective, causes Hirschprung's disease (HSCR, colon aganglionosis). RET mutations alone cannot account for the variable HSCR phenotype, invoking interactions with as yet unknown, and probably inter-related, loci involved in ENS development. Homeobox (HOX) genes have a major role in gut development as depicted by the enteric Hox code. We investigated whether DNA alterations in HOX genes, either alone or in combination with RET, are implicated in HSCR. Genotyping effort was minimized by applying the HapMap data on Han Chinese from Beijing (CHB). 194 HSCR patients and 168 controls were genotyped using Sequenom technology for 72 tag, single nucleotide polymorphisms (SNPs) distributed along the HOX clusters. The HapMap frequencies were compared to those in our population and standard statistics were used for frequency comparisons. The multifactor-dimensionality-reduction method was used for multilocus analysis, in which RET promoter SNP genotypes were included. Genetic interactions were found between two HOX loci (5'-HOXA13 and 3'UTR-HOXB7) and the RET loci tested. Minor allele frequencies (MAF) of the SNPs tested in our sample were not significantly different from those reported by HapMap when the sample sizes of the populations compared were considered. This is the first evaluation of the HOX genes in HSCR and the first application of HapMap data in a Chinese population. The interacting HOX loci may affect the penetrance of the RET risk allele. HapMap data for the CHB population correlated well with the general Chinese population.

Dopamine D3 receptor gene and tardive dyskinesia in Chinese schizophrenic patients

Epidemiological studies have shown a lower prevalence of tardive dyskinesia (TD) among Chinese psychiatric patients compared to Caucasian and Black patient populations. It has been hypothesized that pharmacogenetic factors may underlie this cross-cultural difference. Due to the important implications of the dopamine D3 receptor gene (DRD3) in motor control, we investigated the frequency of polymorphic serine (ser) to glycine (gly) substitution of the gene DRD3 in Chinese schizophrenic patients. The sample size consisted of 65 patients with TD and 66 without TD. Patients were assessed for the severity of TD, the presence of akathisia and parkinsonian symptoms and were subsequently genotyped. We found no evidence that the dopamine D3 receptor gene is likely to confer susceptibility to the development of tardive dyskinesia in Chinese patients with schizophrenia.