Deviation from baseline mutation burden provides powerful and robust rare-variants association test for complex diseases.

Identifying rare variants that contribute to complex diseases is challenging because of the low statistical power in current tests comparing cases with controls. Here, we propose a novel and powerful rare variants association test based on the deviation of the observed mutation burden of a gene in cases from a baseline predicted by a weighted recursive truncated negative-binomial regression (RUNNER) on genomic features available from public data. Simulation studies show that RUNNER is substantially more powerful than state-of-the-art rare variant association tests and has reasonable type 1 error rates even for stratified populations or in small samples. Applied to real case-control data, RUNNER recapitulates known genes of Hirschsprung disease and Alzheimer's disease missed by current methods and detects promising new candidate genes for both disorders. In a case-only study, RUNNER successfully detected a known causal gene of amyotrophic lateral sclerosis. The present study provides a powerful and robust method to identify susceptibility genes with rare risk variants for complex diseases.

Actionable pharmacogenetic variants in Hong Kong Chinese exome sequencing data and projected prescription impact in the Hong Kong population.

Preemptive pharmacogenetic testing has the potential to improve drug dosing by providing point-of-care patient genotype information. Nonetheless, its implementation in the Chinese population is limited by the lack of population-wide data. In this study, secondary analysis of exome sequencing data was conducted to study pharmacogenomics in 1116 Hong Kong Chinese. We aimed to identify the spectrum of actionable pharmacogenetic variants and rare, predicted deleterious variants that are potentially actionable in Hong Kong Chinese, and to estimate the proportion of dispensed drugs that may potentially benefit from genotype-guided prescription. The projected preemptive pharmacogenetic testing prescription impact was evaluated based on the patient prescription data of the public healthcare system in 2019, serving 7.5 million people. Twenty-nine actionable pharmacogenetic variants/ alleles were identified in our cohort. Nearly all (99.6%) subjects carried at least one actionable pharmacogenetic variant, whereas 93.5% of subjects harbored at least one rare deleterious pharmacogenetic variant. Based on the prescription data in 2019, 13.4% of the Hong Kong population was prescribed with drugs with pharmacogenetic clinical practice guideline recommendations. The total expenditure on actionable drugs was 33,520,000 USD, and it was estimated that 8,219,000 USD (24.5%) worth of drugs were prescribed to patients with an implicated actionable phenotype. Secondary use of exome sequencing data for pharmacogenetic analysis is feasible, and preemptive pharmacogenetic testing has the potential to support prescription decisions in the Hong Kong Chinese population.

Whole-genome analysis of noncoding genetic variations identifies multiscale regulatory element perturbations associated with Hirschsprung disease.

It is widely recognized that noncoding genetic variants play important roles in many human diseases, but there are multiple challenges that hinder the identification of functional disease-associated noncoding variants. The number of noncoding variants can be many times that of coding variants; many of them are not functional but in linkage disequilibrium with the functional ones; different variants can have epistatic effects; different variants can affect the same genes or pathways in different individuals; and some variants are related to each other not by affecting the same gene but by affecting the binding of the same upstream regulator. To overcome these difficulties, we propose a novel analysis framework that considers convergent impacts of different genetic variants on protein binding, which provides multiscale information about disease-associated perturbations of regulatory elements, genes, and pathways. Applying it to our whole-genome sequencing data of 918 short-segment Hirschsprung disease patients and matched controls, we identify various novel genes not detected by standard single-variant and region-based tests, functionally centering on neural crest migration and development. Our framework also identifies upstream regulators whose binding is influenced by the noncoding variants. Using human neural crest cells, we confirm cell stage-specific regulatory roles of three top novel regulatory elements on our list, respectively in the RET, RASGEF1A, and PIK3C2B loci. In the PIK3C2B regulatory element, we further show that a noncoding variant found only in the patients affects the binding of the gliogenesis regulator NFIA, with a corresponding up-regulation of multiple genes in the same topologically associating domain.

Assessing the safety of lipid-modifying medications among Chinese adolescents: a drug-target Mendelian randomization study.

BACKGROUND: With increasing hypercholesterolemia prevalence in East Asian adolescents, pharmacologic interventions (e.g., HMGCR inhibitors (statins) and PCSK9 inhibitors) may have to be considered although their longer-term safety in the general adolescent population is unclear. This study aims to investigate the longer-term safety of HMGCR inhibitors and PCSK9 inhibitors among East Asian adolescents using genetics. METHODS: A drug-target Mendelian randomization study leveraging the Global Lipid Genetics Consortium (East Asian, n = 146,492) and individual-level data from Chinese participants in the Biobank clinical follow-up of Hong Kong's "Children of 1997" birth cohort (n = 3443, aged ~ 17.6 years). Safety outcomes (n = 100) included anthropometric and hematological traits, renal, liver, lung function, and other nuclear magnetic resonance metabolomics. Positive control outcomes were cholesterol markers from the "Children of 1997" birth cohort and coronary artery disease from Biobank Japan. RESULTS: Genetic inhibition of HMGCR and PCSK9 were associated with reduction in cholesterol-related NMR metabolomics, e.g., apolipoprotein B (HMGCR: beta [95% CI], - 1.06 [- 1.52 to - 0.60]; PCSK9: - 0.93 [- 1.56 to - 0.31]) and had the expected effect on the positive control outcomes. After correcting for multiple comparisons (p-value < 0.006), genetic inhibition of HMGCR was associated with lower linoleic acid - 0.79 [- 1.25 to - 0.35]. Genetic inhibition of PCSK9 was not associated with the safety outcomes assessed. CONCLUSIONS: Statins and PCSK9 inhibitors in East Asian adolescents appeared to be safe based on the outcomes concerned. Larger studies were warranted to verify these findings. This study serves as a proof of principle study to inform the medication safety among adolescents via genetics.

A polygenic risk score improves risk stratification of coronary artery disease: a large-scale prospective Chinese cohort study.

AIMS: To construct a polygenic risk score (PRS) for coronary artery disease (CAD) and comprehensively evaluate its potential in clinical utility for primary prevention in Chinese populations. METHODS AND RESULTS: Using meta-analytic approach and large genome-wide association results for CAD and CAD-related traits in East Asians, a PRS comprising 540 genetic variants was developed in a training set of 2800 patients with CAD and 2055 controls, and was further assessed for risk stratification for CAD integrating with the guideline-recommended clinical risk score in large prospective cohorts comprising 41 271 individuals. During a mean follow-up of 13.0 years, 1303 incident CAD cases were identified. Individuals with high PRS (the highest 20%) had about three-fold higher risk of CAD than the lowest 20% (hazard ratio 2.91, 95% confidence interval 2.43-3.49), with the lifetime risk of 15.9 and 5.8%, respectively. The addition of PRS to the clinical risk score yielded a modest yet significant improvement in C-statistic (1%) and net reclassification improvement (3.5%). We observed significant gradients in both 10-year and lifetime risk of CAD according to the PRS within each clinical risk strata. Particularly, when integrating high PRS, intermediate clinical risk individuals with uncertain clinical decision for intervention would reach the risk levels (10-year of 4.6 vs. 4.8%, lifetime of 17.9 vs. 16.6%) of high clinical risk individuals with intermediate (20-80%) PRS. CONCLUSION: The PRS could stratify individuals into different trajectories of CAD risk, and further refine risk stratification for CAD within each clinical risk strata, demonstrating a great potential to identify high-risk individuals for targeted intervention in clinical utility.

Genetics of Hirschsprung's disease.

Hirschsprung's disease (HSCR) is a classical model of enteric neuropathy, occurring in approximately 2-2.8 in 10,000 newborns. It is the commonest form of congenital bowel obstruction and is characterized by the absence of enteric ganglia in distal colon. Recent advances in genome-wide association analysis (GWAS) and next generation sequencing (NGS) studies have led to the discovery of a number of new HSCR candidate genes, thereby providing new insights into the genetic architecture and molecular mechanisms of the disease. Altogether, these findings indicated that genetic heterogeneity, variable penetrance and expressivity, and genetic interaction are the pervasive characteristics of HSCR genetics. In this review, we will provide an update on the genetic landscape of HSCR and discuss how the common and rare variants may act together to modulate the phenotypic manifestation. Translating the genetic findings to genetic risk prediction and to optimize clinical outcomes are undoubtedly the ultimate goals for genetic studies on HSCR. From this perspective, we will further discuss the major obstacles in the clinical translation of these latest genetic findings. Lastly, new measures to address these clinical challenges are suggested to advance precision medicine and to develop novel alternative therapies.

Current Understanding in the Clinical Characteristics and Molecular Mechanisms in Different Subtypes of Biliary Atresia.

Biliary atresia is a severe obliterative cholangiopathy in early infancy that is by far the most common cause of surgical jaundice and the most common indicator for liver transplantation in children. With the advanced knowledge gained from different clinical trials and the development of research models, a more precise clinical classification of BA (i.e., isolated BA (IBA), cystic BA (CBA), syndromic BA (SBA), and cytomegalovirus-associated BA (CMVBA)) is proposed. Different BA subtypes have similar yet distinguishable clinical manifestations. The clinical and etiological heterogeneity leads to dramatically different prognoses; hence, treatment needs to be specific. In this study, we reviewed the clinical characteristics of different BA subtypes and revealed the molecular mechanisms of their developmental contributors. We aimed to highlight the differences among these various subtypes of BA which ultimately contribute to the development of a specific management protocol for each subtype.

Isogenic Human-Induced Pluripotent Stem-Cell-Derived Cardiomyocytes Reveal Activation of Wnt Signaling Pathways Underlying Intrinsic Cardiac Abnormalities in Rett Syndrome.

Rett syndrome (RTT) is a severe neurodevelopmental disorder caused by MeCP2 mutations. Nonetheless, the pathophysiological roles of MeCP2 mutations in the etiology of intrinsic cardiac abnormality and sudden death remain unclear. In this study, we performed a detailed functional studies (calcium and electrophysiological analysis) and RNA-sequencing-based transcriptome analysis of a pair of isogenic RTT female patient-specific induced pluripotent stem-cell-derived cardiomyocytes (iPSC-CMs) that expressed either MeCP2wildtype or MeCP2mutant allele and iPSC-CMs from a non-affected female control. The observations were further confirmed by additional experiments, including Wnt signaling inhibitor treatment, siRNA-based gene silencing, and ion channel blockade. Compared with MeCP2wildtype and control iPSC-CMs, MeCP2mutant iPSC-CMs exhibited prolonged action potential and increased frequency of spontaneous early after polarization. RNA sequencing analysis revealed up-regulation of various Wnt family genes in MeCP2mutant iPSC-CMs. Treatment of MeCP2mutant iPSC-CMs with a Wnt inhibitor XAV939 significantly decreased the ß-catenin protein level and CACN1AC expression and ameliorated their abnormal electrophysiological properties. In summary, our data provide novel insight into the contribution of activation of the Wnt/ß-catenin signaling cascade to the cardiac abnormalities associated with MeCP2 mutations in RTT.

Actionable secondary findings in 1116 Hong Kong Chinese based on exome sequencing data.

The use of exome and genome sequencing has increased rapidly nowadays. After primary analysis, further analysis can be performed to identify secondary findings that offer medical benefit for patient care. Multiple studies have been performed to evaluate secondary findings in different ethnicities. However, relevant data are limited in Chinese. Here, with the use of a cohort consisted of 1116 Hong Kong Chinese exome sequencing data, we evaluated the secondary findings in the 59 genes recommended by the American College of Medical Genetics and Genomics. Fifteen unique pathogenic or likely pathogenic variants in 17 individuals were identified, representing a frequency of 1.52% in our cohort. Although 20 individuals harboured pathogenic or likely pathogenic variants in recessive conditions, none carried bi-allelic mutations in the same gene. Our finding was in accordance with the estimation from the American College of Medical Genetics and Genomics that about 1% individuals harbour secondary findings.

Roles of Enteric Neural Stem Cell Niche and Enteric Nervous System Development in Hirschsprung Disease.

The development of the enteric nervous system (ENS) is highly modulated by the synchronized interaction between the enteric neural crest cells (ENCCs) and the neural stem cell niche comprising the gut microenvironment. Genetic defects dysregulating the cellular behaviour(s) of the ENCCs result in incomplete innervation and hence ENS dysfunction. Hirschsprung disease (HSCR) is a rare complex neurocristopathy in which the enteric neural crest-derived cells fail to colonize the distal colon. In addition to ENS defects, increasing evidence suggests that HSCR patients may have intrinsic defects in the niche impairing the extracellular matrix (ECM)-cell interaction and/or dysregulating the cellular niche factors necessary for controlling stem cell behaviour. The niche defects in patients may compromise the regenerative capacity of the stem cell-based therapy and advocate for drug- and niche-based therapies as complementary therapeutic strategies to alleviate/enhance niche-cell interaction. Here, we provide a summary of the current understandings of the role of the enteric neural stem cell niche in modulating the development of the ENS and in the pathogenesis of HSCR. Deciphering the contribution of the niche to HSCR may provide important implications to the development of regenerative medicine for HSCR.

Identification of Genes Associated With Hirschsprung Disease, Based on Whole-Genome Sequence Analysis, and Potential Effects on Enteric Nervous System Development.

BACKGROUND & AIMS: Hirschsprung disease, or congenital aganglionosis, is believed to be oligogenic-that is, caused by multiple genetic factors. We performed whole-genome sequence analyses of patients with Hirschsprung disease to identify genetic factors that contribute to disease development and analyzed the functional effects of these variants. METHODS: We performed whole-genome sequence analyses of 443 patients with short-segment disease, recruited from hospitals in China and Vietnam, and 493 ethnically matched individuals without Hirschsprung disease (controls). We performed genome-wide association analyses and gene-based rare-variant burden tests to identify rare and common disease-associated variants and study their interactions. We obtained induced pluripotent stem cell (iPSC) lines from 4 patients with Hirschsprung disease and 2 control individuals, and we used these to generate enteric neural crest cells for transcriptomic analyses. We assessed the neuronal lineage differentiation capability of iPSC-derived enteric neural crest cells using an in vitro differentiation assay. RESULTS: We identified 4 susceptibility loci, including 1 in the phospholipase D1 gene (PLD1) (P = 7.4 × 10-7). The patients had a significant excess of rare protein-altering variants in genes previously associated with Hirschsprung disease and in the ß-secretase 2 gene (BACE2) (P = 2.9 × 10-6). The epistatic effects of common and rare variants across these loci provided a sensitized background that increased risk for the disease. In studies of the iPSCs, we observed common and distinct pathways associated with variants in RET that affect risk. In functional assays, we found variants in BACE2 to protect enteric neurons from apoptosis. We propose that alterations in BACE1 signaling via amyloid ß precursor protein and BACE2 contribute to pathogenesis of Hirschsprung disease. CONCLUSIONS: In whole-genome sequence analyses of patients with Hirschsprung disease, we identified rare and common variants associated with disease risk. Using iPSC cells, we discovered some functional effects of these variants.

Whole-exome sequencing identifies MST1R as a genetic susceptibility gene in nasopharyngeal carcinoma.

Multiple factors, including host genetics, environmental factors, and Epstein-Barr virus (EBV) infection, contribute to nasopharyngeal carcinoma (NPC) development. To identify genetic susceptibility genes for NPC, a whole-exome sequencing (WES) study was performed in 161 NPC cases and 895 controls of Southern Chinese descent. The gene-based burden test discovered an association between macrophage-stimulating 1 receptor (MST1R) and NPC. We identified 13 independent cases carrying the MST1R pathogenic heterozygous germ-line variants, and 53.8% of these cases were diagnosed with NPC aged at or even younger than 20 y, indicating that MST1R germline variants are relevant to disease early-age onset (EAO) (age of ≤20 y). In total, five MST1R missense variants were found in EAO cases but were rare in controls (EAO vs. control, 17.9% vs. 1.2%, P = 7.94 × 10(-12)). The validation study, including 2,160 cases and 2,433 controls, showed that the MST1R variant c.G917A:p.R306H is highly associated with NPC (odds ratio of 9.0). MST1R is predominantly expressed in the tissue-resident macrophages and is critical for innate immunity that protects organs from tissue damage and inflammation. Importantly, MST1R expression is detected in the ciliated epithelial cells in normal nasopharyngeal mucosa and plays a role in the cilia motility important for host defense. Although no somatic mutation of MST1R was identified in the sporadic NPC tumors, copy number alterations and promoter hypermethylation at MST1R were often observed. Our findings provide new insights into the pathogenesis of NPC by highlighting the involvement of the MST1R-mediated signaling pathways.

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.

Actionable secondary findings from whole-genome sequencing of 954 East Asians.

Recently, the American College of Medical Genetics (ACMG) recommended the return of actionable secondary findings detected from clinical sequencing. The reported frequency of secondary findings in Asian populations were highly variable and it is unclear whether the uniformity in coverage offered by whole-genome sequencing (WGS) may impact the estimate. In this analysis, we aimed to refine the rate of secondary findings on East Asians through a large-scale WGS study. We classified 1256 protein-altering or splicing variants of the 59 actionable genes detected from WGS of 954 East Asians in strict accordance with the ACMG and the Association for Molecular Pathology guidelines. A total of 21 pathogenic or likely pathogenic variants were detected in 24 of the 954 East Asian genomes with an estimate of 2.5% of East Asians carrying actionable variants. Although the overall estimate of secondary findings was consistent with those reported for non-East Asian ethnicities, genetic and allelic heterogeneity was observed. WGS offers a wider breadth of coverage over WES, which highlights the need to further investigate the variable sensitivity of WES and WGS in the detection of secondary findings. Identifying secondary findings in populations underrepresented in previous genetic literature might improve variant interpretation and has a profound impact on local decision-making with regard to the cost-effectiveness of returning the secondary findings from clinical sequencing.

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.

Targeted next-generation sequencing on Hirschsprung disease: a pilot study exploits DNA pooling.

To adopt an efficient approach of identifying rare variants possibly related to Hirschsprung disease (HSCR), a pilot study was set up to evaluate the performance of a newly designed protocol for next generation targeted resquencing. In total, 20 Chinese HSCR patients and 20 Chinese sex-matched individuals with no HSCR were included, for which coding sequences (CDS) of 62 genes known to be in signaling pathways relevant to enteric nervous system development were selected for capture and sequencing. Blood DNAs from eight pools of five cases or controls were enriched by PCR-based RainDance technology (RDT) and then sequenced on a 454 FLX platform. As technical validation, five patients from case Pool-3 were also independently enriched by RDT, indexed with barcode and sequenced with sufficient coverage. Assessment for CDS single nucleotide variants showed DNA pooling performed well (specificity/sensitivity at 98.4%/83.7%) at the common variant level; but relatively worse (specificity/sensitivity at 65.5%/61.3%) at the rare variant level. Further Sanger sequencing only validated five out of 12 rare damaging variants likely involved in HSCR. Hence more improvement at variant detection and sequencing technology is needed to realize the potential of DNA pooling for large-scale resequencing projects.

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.

Artificial intelligence-based approaches for the detection and prioritization of genomic mutations in congenital surgical diseases.

Genetic mutations are critical factors leading to congenital surgical diseases and can be identified through genomic analysis. Early and accurate identification of genetic mutations underlying these conditions is vital for clinical diagnosis and effective treatment. In recent years, artificial intelligence (AI) has been widely applied for analyzing genomic data in various clinical settings, including congenital surgical diseases. This review paper summarizes current state-of-the-art AI-based approaches used in genomic analysis and highlighted some successful applications that deepen our understanding of the etiology of several congenital surgical diseases. We focus on the AI methods designed for the detection of different variant types and the prioritization of deleterious variants located in different genomic regions, aiming to uncover susceptibility genomic mutations contributed to congenital surgical disorders.

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.

Mutations in the NRG1 gene are associated with Hirschsprung disease.

Hirschsprung disease (HSCR, congenital colon aganglionosis) is a relatively common complex genetic condition caused by abnormal development of the enteric nervous system (ENS). Through a recent genome-wide association study conducted on Chinese HSCR patients, we identified a new HSCR contributing locus, neuregulin 1 (NRG1; 8p12), a gene known to be involved in the development of the ENS. As genes in which disease-associated common variants are found are to be considered as candidates for the search of deleterious rare variants (RVs) in the coding sequences, we sequenced the NRG1 exons of 358 sporadic HSCR patients and 333 controls. We identified a total of 13 different heterozygous RVs including 8 non-synonymous (A28G, E134K, V266L, H347Y, P356L, V486M, A511T, P608A) and 3 synonymous amino acid substitutions (P24P, T169T, L483L), a frameshift (E239fsX10), and a c.503-4insT insertion. Functional analysis of the most conserved non-synonymous substitutions, H347Y and P356L, showed uneven intracellular distribution and aberrant expression of the mutant proteins. Except for T169T and V486M, all variants were exclusive to HSCR patients. Overall, there was a statistically significant over-representation of NRG1 RVs in HSCR patients (p = 0.008). We show here that not only common, but also rare variants of the NRG1 gene contribute to HSCR. This strengthens the role of NRG1.