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

Distinctive genetic variation of long-segment Hirschsprung's disease in Taiwan

BACKGROUND

Hirschsprung's disease (HSCR) is a congenital disorder with the absence of myenteric and submucosal ganglion cells within distal gut. Due to multigenic inheritance and interactions, we employed next-generation sequencing (NGS) to investigate genetic backgrounds of long-segment HSCR (L-HSCR) in Taiwan.

METHODS

Genomic DNA extracted from peripheral blood of L-HSCR patients was subjected to capture-based NGS, based on a 31-gene panel. The variants with allele frequency <0.05 and predicted by computational methods as deleterious were further validated by Sanger sequencing in patients and their family as well to tell de novo from inherited variants.

RESULTS

Between 2015/04 and 2018/05, this study enrolled 23 L-HSCR patients, including 15 (65.2%) sporadic cases and 8 (34.8%) familial patients in 4 different families. Six sporadic and seven familial cases showed possible harmful variants across eight different genes, accounting for an overall detection rate of 56.5%. These variants mainly resided in SEMA3C, followed by RET, NRG1, and NTRK1. Three sporadic and 2 familial cases exhibited strong pathogenic variants as a deletional frameshift or stop codon in RET, L1CAM or NRG1. In a HSCR family, the father passed on a pathogenic RET frameshift to two daughters; however, only one developed HSCR.

CONCLUSION

Using NGS, we disclosed deleterious mutations such as a frameshift or stop codon in either familial or sporadic patients. Our cases with isolated L-HSCR or even total colonic aganglionosis appeared to exhibit complex patterns of inheritance and incomplete penetrance even in families with the same genetic variants, reflecting the possible effects of environmental factors and genetic modifiers.

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.

Beta-Binomial Model for the Detection of Rare Mutations in Pooled Next-Generation Sequencing Experiments

Against diminishing costs, next-generation sequencing (NGS) still remains expensive for studies with a large number of individuals. As cost saving, sequencing genome of pools containing multiple samples might be used. Currently, there are many software available for the detection of single-nucleotide polymorphisms (SNPs). Sensitivity and specificity depend on the model used and data analyzed, indicating that all software have space for improvement. We use beta-binomial model to detect rare mutations in untagged pooled NGS experiments. We propose a multireference framework for pooled data with ability being specific up to two patients affected by neuromuscular disorders (NMD). We assessed the results comparing with The Genome Analysis Toolkit (GATK), CRISP, SNVer, and FreeBayes. Our results show that the multireference approach applying beta-binomial model is accurate in predicting rare mutations at 0.01 fraction. Finally, we explored the concordance of mutations between the model and software, checking their involvement in any NMD-related gene. We detected seven novel SNPs, for which the functional analysis produced enriched terms related to locomotion and musculature.

Exome-Wide Association Study Identified New Risk Loci for Hirschsprung's Disease

Hirschsprung disease (HSCR) is a rare congenital disease caused by impaired proliferation and migration of neural crest cells. In this study, we aimed to investigate the genetic loci involved in the pathogenesis of HSCR. The exome-wide scan was performed to screen the genetic variants with minor allele frequency (MAF) < 0.05 in exonic regions. Candidate mutation type and the wild type were overexpressed to investigate the affection on cell proliferation and migration. We found that ten variants were associated with HSCR at P < 10^-4 in the single-variant analysis while ten genes were also associated with HSCR at P < 10^-4 in the optimized sequence kernel association test (SKAT-O) test analysis. Among these SNPs, the missense variants catechol-O-methyltransferase (COMT) (rs6267) and armadillo repeat gene deleted in velocardiofacial syndrome (ARVCF) (rs80068543) indicated an ectopic expression in colon tissues of HSCR patients. The Ala72Ser variant in COMT induced proliferation suppression through NOTCH signal pathway, while the ARVCF affected cell migration via the downregulating of RHOA and ROC. In conclusion, this exome array study identified the COMT and ARVCF missense coding variants as candidate loci for HSCR. The finding implies the abnormal variant of COMT and ARVCF may account for the pathogenesis of HSCR.

Next-generation-based targeted sequencing as an efficient tool for the study of the genetic background in Hirschsprung patients

BACKGROUND

The development of next-generation sequencing (NGS) technologies has a great impact in the human variation detection given their high-throughput. These techniques are particularly helpful for the evaluation of the genetic background in disorders of complex genetic etiology such as Hirschsprung disease (HSCR). The purpose of this study was the design of a panel of HSCR associated genes as a rapid and efficient tool to perform genetic screening in a series of patients.

METHODS

We have performed NGS-based targeted sequencing (454-GS Junior) using a panel containing 26 associated or candidate genes for HSCR in a group of 11 selected HSCR patients.

RESULTS

The average percentage of covered bases was of 97%, the 91.4% of the targeted bases were covered with depth above 20X and the mean coverage was 422X. In addition, we have found a total of 13 new coding variants and 11 new variants within regulatory regions among our patients. These outcomes allowed us to re-evaluate the genetic component associated to HSCR in these patients.

CONCLUSIONS

Our validated NGS panel constitutes an optimum method for the identification of new variants in our patients. This approach could be used for a fast, reliable and more thorough genetic screening in future series of patients.

A permutation test for oligoset DNA pooling studies

Case-control association studies often suffer from population stratification bias. A previous triple combination strategy of stratum matching, genomic controlling, and multiple DNA pooling can correct the bias and save genotyping cost. However the method requires researchers to prepare a multitude of DNA pools-more than 30 case-control pooling sets in total (polyset). In this paper, the authors propose a permutation test for oligoset DNA pooling studies. Monte-Carlo simulations show that the proposed test has a type I error rate under control and a power comparable to that of individual genotyping. For a researcher on a tight budget, oligoset DNA pooling is a viable option.