Exome Sequencing Identified NRG3 as a Novel Susceptible Gene of Hirschsprung’s Disease in a Chinese Population

Hirschsprung’s disease (HSCR) is a complex developmental defect characterized by the absence of enteric ganglia in the gastrointestinal tract. Although the genetic defect of enteric nervous system (ENS) was identified to play a critical role in the progress of HSCR, the systemic genetic dissection of HSCR still needs to be clarified. In this study, we firstly performed exome sequencing of two HSCR patients from a Han Chinese family, including the affected mother and son. After the initial quality filtering (coverage?≥?5X and SNP quality score?≥?40) of the raw data, we identified 13,948 and 13,856 single nucleotide variants (SNVs), respectively. We subsequently compared the SNVs against public databases (dbSNP130, HapMap, and 1000 Genome Project) and obtained a total of 15 novel nonsynonymous SNVs in 15 genes, which were shared between these two patients. Follow-up Sanger sequencing and bioinformatics analysis highlighted variant c.853G>A (p.E285K) in NRG3, a gene involved in the development of ENS. In the validation phase, we sequenced all nine exons of NRG3 in 96 additional sporadic HSCR cases and 110 healthy individuals and identified another nonsynonymous variant c.1329G>A (p.M443I) and two synonymous variants c.828G>A (p.T276T) and c.1365T>A (p.P455P) only in the cases. Our results indicated that NRG3 may be a susceptibility gene for HSCR in a Chinese population.     

Molecular evolution of the major capsid protein VP1 of enterovirus 70.

Nucleotide sequences of the genome RNA encoding capsid protein VP1 (918 nucleotides) of 18 enterovirus 70 (EV70) isolates collected from various parts of the world in 1971 to 1981 were determined, and nucleotide substitutions among them were studied. The genetic distances between isolates were calculated by the pairwise comparison of nucleotide difference. Regression analysis of the genetic distances against time of isolation of the strains showed that the synonymous substitution rate was very high at 21.53 x 10(-3) substitution per nucleotide per year, while the nonsynonymous rate was extremely low at 0.32 x 10(-3) substitution per nucleotide per year. The rate estimated by the average value of synonymous and nonsynonymous substitutions (W.-H. Li, C.-C. Wu, and C.-C. Luo, Mol. Biol. Evol. 2:150-174, 1985) was 5.00 x 10(-3) substitution per nucleotide per year. Taking the average value of synonymous and nonsynonymous substitutions as genetic distances between isolates, the phylogenetic tree was inferred by the unweighted pairwise grouping method of arithmetic average and by the neighbor-joining method. The tree indicated that the virus had evolved from one focal place, and the time of emergence was estimated to be August 1967 +/- 15 months, 2 years before first recognition of the pandemic of acute hemorrhagic conjunctivitis. By superimposing every nucleotide substitution on the branches of the phylogenetic tree, we analyzed nucleotide substitution patterns of EV70 genome RNA. In synonymous substitutions, the proportion of transitions, i.e., C<==>U and G<==>A, was found to be extremely frequent in comparison with that reported on other viruses or pseudogenes. In addition, parallel substitutions (independent substitutions at the same nucleotide position on different branches, i.e., different isolates, of the tree) were frequently found in both synonymous and nonsynonymous substitutions. These frequent parallel substitutions and the low nonsynonymous substitution rate despite the very high synonymous substitution rate described above imply a strong restriction on nonsynonymous substitution sites of VP1, probably due to the requirement for maintaining the rigid icosahedral conformation of the virus.     

Biallelic SUN5 Mutations Cause Autosomal-Recessive Acephalic Spermatozoa Syndrome

Acephalic spermatozoa syndrome is a rare and severe form of teratozoospermia characterized by a predominance of headless spermatozoa in the ejaculate. Family clustering and consanguinity suggest a genetic origin; however, causative mutations have yet to be identified. We performed whole-exome sequencing in two unrelated infertile men and subsequent variant filtering identified one homozygous (c.824C>T [p.Thr275Met]) and one compound heterozygous (c.1006C>T [p.Arg356Cys] and c.485T>A [p.Met162Lys]) SUN5 (also named TSARG4) variants. Sanger sequencing of SUN5 in 15 additional unrelated infertile men revealed four compound heterozygous (c.381delA [p.Val128Serfs^∗7] and c.824C>T [p.Thr275Met]; c.381delA [p.Val128Serfs^∗7] and c.781G>A [p.Val261Met]; c.216G>A [p.Trp72^∗] and c.1043A>T [p.Asn348Ile]; c.425+1G>A/c.1043A>T [p.Asn348Ile]) and two homozygous (c.851C>G [p.Ser284^∗]; c.350G>A [p.Gly114Arg]) variants in six individuals. These 10 SUN5 variants were found in 8 of 17 unrelated men, explaining the genetic defect in 47.06% of the affected individuals in our cohort. These variants were absent in 100 fertile population-matched control individuals. SUN5 variants lead to absent, significantly reduced, or truncated SUN5, and certain variants altered SUN5 distribution in the head-tail junction of the sperm. In summary, these results demonstrate that biallelic SUN5 mutations cause male infertility due to autosomal-recessive acephalic spermatozoa syndrome.     

Mutation Screening of Retinal Dystrophy Patients by Targeted Capture from Tagged Pooled DNAs and Next Generation Sequencing

Purpose

Retinal dystrophies are genetically heterogeneous, resulting from mutations in over 200 genes. Prior to the development of massively parallel sequencing, comprehensive genetic screening was unobtainable for most patients. Identifying the causative genetic mutation facilitates genetic counselling, carrier testing and prenatal/pre-implantation diagnosis, and often leads to a clearer prognosis. In addition, in a proportion of cases, when the mutation is known treatment can be optimised and patients are eligible for enrolment into clinical trials for gene-specific therapies.

Methods

Patient genomic DNA was sheared, tagged and pooled in batches of four samples, prior to targeted capture and next generation sequencing. The enrichment reagent was designed against genes listed on the RetNet database (July 2010). Sequence data were aligned to the human genome and variants were filtered to identify potential pathogenic mutations. These were confirmed by Sanger sequencing.

Results

Molecular analysis of 20 DNAs from retinal dystrophy patients identified likely pathogenic mutations in 12 cases, many of them known and/or confirmed by segregation. These included previously described mutations in ABCA4 (c.6088C>T,p.R2030*; c.5882G>A,p.G1961E), BBS2 (c.1895G>C,p.R632P), GUCY2D (c.2512C>T,p.R838C), PROM1 (c.1117C>T,p.R373C), RDH12 (c.601T>C,p.C201R; c.506G>A,p.R169Q), RPGRIP1 (c.3565C>T,p.R1189*) and SPATA7 (c.253C>T,p.R85*) and new mutations in ABCA4 (c.3328+1G>C), CRB1 (c.2832_2842+23del), RP2 (c.884-1G>T) and USH2A (c.12874A>G,p.N4292D).

Conclusions

Tagging and pooling DNA prior to targeted capture of known retinal dystrophy genes identified mutations in 60% of cases. This relatively high success rate may reflect enrichment for consanguineous cases in the local Yorkshire population, and the use of multiplex families. Nevertheless this is a promising high throughput approach to retinal dystrophy diagnostics.     

Recurrent Muscle Weakness with Rhabdomyolysis,Metabolic Crises,and Cardiac Arrhythmia Due to Bi-allelic TANGO2 Mutations

The underlying genetic etiology of rhabdomyolysis remains elusive in a significant fraction of individuals presenting with recurrent metabolic crises and muscle weakness. Using exome sequencing, we identified bi-allelic mutations in TANGO2 encoding transport and Golgi organization 2 homolog (Drosophila) in 12 subjects with episodic rhabdomyolysis, hypoglycemia, hyperammonemia, and susceptibility to life-threatening cardiac tachyarrhythmias. A recurrent homozygous c.460G>A (p.Gly154Arg) mutation was found in four unrelated individuals of Hispanic/Latino origin, and a homozygous ∼34 kb deletion affecting exons 3–9 was observed in two families of European ancestry. One individual of mixed Hispanic/European descent was found to be compound heterozygous for c.460G>A (p.Gly154Arg) and the deletion of exons 3–9. Additionally, a homozygous exons 4–6 deletion was identified in a consanguineous Middle Eastern Arab family. No homozygotes have been reported for these changes in control databases. Fibroblasts derived from a subject with the recurrent c.460G>A (p.Gly154Arg) mutation showed evidence of increased endoplasmic reticulum stress and a reduction in Golgi volume density in comparison to control. Our results show that the c.460G>A (p.Gly154Arg) mutation and the exons 3–9 heterozygous deletion in TANGO2 are recurrent pathogenic alleles present in the Latino/Hispanic and European populations, respectively, causing considerable morbidity in the homozygotes in these populations.     

Association between Rare Variants in AP4E1, a Component of Intracellular Trafficking,and Persistent Stuttering

Stuttering is a common, highly heritable neurodevelopmental disorder characterized by deficits in the volitional control of speech. Whole-exome sequencing identified two heterozygous AP4E1 coding variants, c.1549G>A (p.Val517Ile) and c.2401G>A (p.Glu801Lys), that co-segregate with persistent developmental stuttering in a large Cameroonian family, and we observed the same two variants in unrelated Cameroonians with persistent stuttering. We found 23 other rare variants, including predicted loss-of-function variants, in AP4E1 in unrelated stuttering individuals in Cameroon, Pakistan, and North America. The rate of rare variants in AP4E1 was significantly higher in unrelated Pakistani and Cameroonian stuttering individuals than in population-matched control individuals, and coding variants in this gene are exceptionally rare in the general sub-Saharan West African, South Asian, and North American populations. Clinical examination of the Cameroonian family members failed to identify any symptoms previously reported in rare individuals carrying homozygous loss-of-function mutations in this gene. AP4E1 encodes the ε subunit of the heterotetrameric (ε-β4-μ4-σ4) AP-4 complex, involved in protein sorting at the trans-Golgi network. We found that the μ4 subunit of AP-4 interacts with NAGPA, an enzyme involved in the synthesis of the mannose 6-phosphate signal that targets acid hydrolases to the lysosome and the product of a gene previously associated with stuttering. These findings implicate deficits in intracellular trafficking in persistent stuttering.     

Variants of the PPARD Gene and Their Clinicopathological Significance in Colorectal Cancer

Background

Peroxisome proliferator-activated receptor delta (PPARD) is nuclear hormone receptor involved in colorectal cancer (CRC) differentiation and progression. The purpose of this study was to determine prevalence and spectrum of variants in the PPARD gene in CRC, and their contribution to clinicopathological endpoints.

Methods and Findings

Direct sequencing of the PPARD gene was performed in 303 primary tumors, in blood samples from 50 patients with ≥3 affected first-degree relatives, 50 patients with 2 affected first-degree relatives, 50 sporadic patients, 360 healthy controls, and in 6 colon cancer cell lines. Mutation analysis revealed 22 different transversions, 7 of them were novel. Three of all variants were somatic (c.548A>G, p.Y183C, c.425-9C>T, and c.628-16G>A). Two missense mutations (p.Y183C and p.R258Q) were pathogenic using in silico predictive program. Five recurrent variants were detected in/adjacent to the exons 4 (c.1-87T>C, c.1-67G>A, c.130+3G>A, and c.1-101-8C>T) and exon 7 (c.489T>C). Variant c.489C/C detected in tumors was correlated to worse differentiation (P = 0.0397).

Conclusions

We found 7 novel variants among 22 inherited or acquired PPARD variants. Somatic and/or missense variants detected in CRC patients are rare but indicate the clinical importance of the PPARD gene.     

Biallelic Variants in TTLL5, Encoding a Tubulin Glutamylase,Cause Retinal Dystrophy

In a subset of inherited retinal degenerations (including cone, cone-rod, and macular dystrophies), cone photoreceptors are more severely affected than rods; ABCA4 mutations are the most common cause of this heterogeneous class of disorders. To identify retinal-disease-associated genes, we performed exome sequencing in 28 individuals with “cone-first” retinal disease and clinical features atypical for ABCA4 retinopathy. We then conducted a gene-based case-control association study with an internal exome data set as the control group. TTLL5, encoding a tubulin glutamylase, was highlighted as the most likely disease-associated gene; 2 of 28 affected subjects harbored presumed loss-of-function variants: c.[1586_1589delAGAG];[1586_1589delAGAG], p.[Glu529Valfs^∗2];[Glu529Valfs^∗2], and c.[401delT(;)3354G>A], p.[Leu134Argfs^∗45(;)Trp1118^∗]. We then inspected previously collected exome sequence data from individuals with related phenotypes and found two siblings with homozygous nonsense variant c.1627G>T (p.Glu543^∗) in TTLL5. Subsequently, we tested a panel of 55 probands with retinal dystrophy for TTLL5 mutations; one proband had a homozygous missense change (c.1627G>A [p.Glu543Lys]). The retinal phenotype was highly similar in three of four families; the sibling pair had a more severe, early-onset disease. In human and murine retinae, TTLL5 localized to the centrioles at the base of the connecting cilium. TTLL5 has been previously reported to be essential for the correct function of sperm flagella in mice and play a role in polyglutamylation of primary cilia in vitro. Notably, genes involved in the polyglutamylation and deglutamylation of tubulin have been associated with photoreceptor degeneration in mice. The electrophysiological and fundus autofluorescence imaging presented here should facilitate the molecular diagnosis in further families.     

Mutations in PIEZO2 Cause Gordon Syndrome,Marden-Walker Syndrome,and Distal Arthrogryposis Type 5

Gordon syndrome (GS), or distal arthrogryposis type 3, is a rare, autosomal-dominant disorder characterized by cleft palate and congenital contractures of the hands and feet. Exome sequencing of five GS-affected families identified mutations in piezo-type mechanosensitive ion channel component 2 (PIEZO2) in each family. Sanger sequencing revealed PIEZO2 mutations in five of seven additional families studied (for a total of 10/12 [83%] individuals), and nine families had an identical c.8057G>A (p.Arg2686His) mutation. The phenotype of GS overlaps with distal arthrogryposis type 5 (DA5) and Marden-Walker syndrome (MWS). Using molecular inversion probes for targeted sequencing to screen PIEZO2, we found mutations in 24/29 (82%) DA5-affected families and one of two MWS-affected families. The presence of cleft palate was significantly associated with c.8057G>A (Fisher’s exact test, adjusted p value < 0.0001). Collectively, although GS, DA5, and MWS have traditionally been considered separate disorders, our findings indicate that they are etiologically related and perhaps represent variable expressivity of the same condition.     

DEFB1 polymorphisms are involved in susceptibility to human papillomavirus infection in Brazilian gynaecological patients

The human beta defensin 1 (hBD-1) antimicrobial peptide is a member of the innate immune system known to act in the first line of defence against microorganisms, including viruses such as human papillomavirus (HPV). In this study, five functional polymorphisms (namely g-52G>A, g-44C>G and g-20G>A in the 5’UTR and c.*5G>A and c.*87A>G in the 3’UTR) in the DEFB1 gene encoding for hBD-1 were analysed to investigate the possible involvement of these genetic variants in susceptibility to HPV infection and in the development of HPV-associated lesions in a population of Brazilian women. The DEFB1 g-52G>A and c.*5G>A single-nucleotide polymorphisms (SNPs) and the GCAAA haplotype showed associations with HPV-negative status; in particular, the c.*5G>A SNP was significantly associated after multiple test corrections. These findings suggest a possible role for the constitutively expressed beta defensin-1 peptide as a natural defence against HPV in the genital tract mucosa.     

Development of a mutation hotspot detection kit for the phenylalanine hydroxylase gene by ARMS-PCR combined with fluorescent probe technology

To develop a screening kit for detecting mutation hotspots of the phenylalanine hydroxylase (PAH) gene. Thirteen exons of the PAH gene were sequenced in 84 cases with phenylketonuria (PKU) diagnosed during neonatal genetic and metabolic disease screening in Shaanxi province, and their mutations were analyzed. We designed and developed a screening kit to detect nine mutation sites covering more than 50% of the PAH mutations found in Shaanxi province (c.728G>A, c.1197A>T, c.331C>T, c.1068C>A, c.611A>G, c.1238G>C, c.721C>T, c.442-1G>A, and c.158G>A) by using amplification refractory mutation system-polymerase chain reaction (ARMS-PCR) combined with fluorescent probe technology. Peripheral blood and dried blood samples from PKU families were used for clinical verification of the newly developed kit. PAH gene mutations were detected in 84 children diagnosed with PKU. A total of 159 mutant alleles were identified, consisting of 100 missense mutations, 28 shear mutations, 24 nonsense mutations, and 7 deletion mutations. Exon 7 had the highest mutation frequency (32.08%). Among them, the mutation frequency of p.R243Q was the highest, accounting for 20.13% of all mutations, followed by p.R111X, IVS4-1G>A, EX6-96A>G, and p.R413P; these five loci accounted for 47.17% (75/159) of all mutations. In addition, we identified three previously unreported PAH gene mutations (p.C334X, p.G46D, and p.G256D). Fifteen mutation sites were identified in the 47 PAH carriers identified by next-generation sequencing (NGS), which were verified by the newly developed kit, with an agreement rate of 100%. This newly developed kit based on ARMS-PCR combined with fluorescent probe technology can be used to detect common PAH gene mutations.     

Exome Sequencing Identifies GNB4 Mutations as a Cause of Dominant Intermediate Charcot-Marie-Tooth Disease

Charcot-Marie-Tooth disease (CMT) is a heterogeneous group of inherited neuropathies. Mutations in approximately 45 genes have been identified as being associated with CMT. Nevertheless, the genetic etiologies of at least 30% of CMTs have yet to be elucidated. Using a genome-wide linkage study, we previously mapped a dominant intermediate CMT to chromosomal region 3q28–q29. Subsequent exome sequencing of two affected first cousins revealed heterozygous mutation c.158G>A (p.Gly53Asp) in GNB4, encoding guanine-nucleotide-binding protein subunit beta-4 (Gβ₄), to cosegregate with the CMT phenotype in the family. Further analysis of GNB4 in an additional 88 unrelated CMT individuals uncovered another de novo mutation, c.265A>G (p.Lys89Glu), in this gene in one individual. Immunohistochemistry studies revealed that Gβ₄ was abundant in the axons and Schwann cells of peripheral nerves and that expression of Gβ₄ was significantly reduced in the sural nerve of the two individuals carrying the c.158G>A (p.Gly53Asp) mutation. In vitro studies demonstrated that both the p.Gly53Asp and p.Lys89Glu altered proteins impaired bradykinin-induced G-protein-coupled-receptor (GPCR) signaling, which was facilitated by the wild-type Gβ₄. This study identifies GNB4 mutations as a cause of CMT and highlights the importance of Gβ₄-related GPCR signaling in peripheral-nerve function in humans.     

Association of Mitochondrial DNA Polymerase γ Gene POLG1 Polymorphisms with Parkinsonism in Chinese Populations

Background

Mitochondrial DNA polymerase gamma (POLG1) mutations were associated with levodopa-responsive Parkinsonism. POLG1 gene contains a number of common nonsynonymous SNPs and intronic regulatory SNPs which may have functional consequences. It is of great interest to discover polymorphisms variants associated with Parkinson''s disease (PD), both in isolation and in combination with specific SNPs.

Materials and Methods

We conducted a case-control study and genotyped twenty SNPs and poly-Q polymorphisms of POLG1 gene including in 344 Chinese sporadic PD patients and 154 healthy controls. All the polymorphisms of POLG1 we found in this study were sequenced by PCR products with dye terminator methods using an ABI-3100 sequencer. Hardy-Weinberg equilibrium and linkage disequilibrium (LD) for association between twenty POLG1 SNPs and PD were calculated using the program Haploview.

Principal Results

We provided evidence for strong association of four intronic SNPs of the POLG1 gene (new report: c.2070-12T>A and rs2307439: c.2070-64G>A in intron 11, P = 0.00011, OR = 1.727; rs2302084: c.3105-11T>C and rs2246900: c.3105-36A>G in intron 19, P = 0.00031, OR = 1.648) with PD. However, we did not identify any significant association between ten exonic SNPs of POLG1 and PD. Linkage disequilibrium analysis indicated that c.2070-12T>A and c.2070-64G>A could be parsed into one block as Haplotype 1 as well as c.3105-11T>C and c.3105-36A>G in Haplotype 2. In addition, case and control study on association of POLG1 CAG repeat (poly-Q) alleles with PD showed a significant association (P = 0.03, OR = 2.16) of the non-10/11Q variants with PD. Although intronic SNPs associated with PD didn''t influence POLG1 mRNA alternative splicing, there was a strong association of c.2070-12T>A and c.2070-64G>A with decreased POLG1 mRNA level and protein levels.

Conclusions

Our findings indicate that POLG1 may play a role in the pathogenesis of PD in Chinese populations.     

Association of the FCN2 Gene Single Nucleotide Polymorphisms with Susceptibility to Pulmonary Tuberculosis

Ficolin-2 (FCN2) is an innate immune pattern recognition molecule that can activate the complement pathway, opsonophagocytosis, and elimination of the pathogens. The present study aimed to investigate the association of the FCN2 gene single nucleotide polymorphisms (SNPs) with susceptibility to pulmonary tuberculosis (TB). A total of seven SNPs in exon 8 (+6359 C>T and +6424 G>T) and in the promoter region (-986 G>A, -602 G>A, -557 A>G, -64 A>C and -4 A>G) of the FCN2 gene were genotyped using the PCR amplification and DNA sequencing methods in the healthy controls group (n = 254) and the pulmonary TB group (n = 282). The correlation between SNPs and pulmonary TB was analyzed using the logistic regression method. The results showed that there were no significant differences in the distribution of allelic frequencies of seven SNPs between the pulmonary TB group and the healthy controls group. However, the frequency of the variant homozygous genotype (P = 0.037, -557 A>G; P = 0.038, -64 A>C; P = 0.024, +6424 G>T) in the TB group was significantly lower than the control group. After adjustment for age and gender, these variant homozygous genotypes were found to be recessive models in association with pulmonary TB. In addition, -64 A>C (P = 0.047) and +6424 G>T (P = 0.03) were found to be codominant models in association with pulmonary TB. There was strong linkage disequilibrium (r² > 0.80, P < 0.0001) between 7 SNPs except the -602 G>A site. Therefore, -557 A>G, -64 A>C and +6424 G>T SNPs of the FCN2 gene were correlated with pulmonary TB, and may be protective factors for TB. This study provides a novel idea for the prevention and control of TB transmission from a genetics perspective.