Genetic Background of Hirschsprung Disease: A Bridge Between Basic Science and Clinical Application

Hirschsprung's disease (HSCR) is a congenital disorder, defined by partial or complete loss of the neuronal ganglion cells in the intestinal tract, which is caused by the failure of neural crest cells to migrate completely during intestinal development during fetal life. HSCR has a multifactorial etiology, and genetic factors play a key role in its pathogenesis; these include mutations within several gene loci. These have been identified by screening candidate genes, or by conducting genome wide association (GWAS) studies. However, only a small portion of them have been proposed as major genetic risk factors for the HSCR. In this review, we focus on those genes that have been identified as either low penetrant or high penetrant variants that determine the risk of Hirschsprung's disease. J. Cell. Biochem. 119: 28–33, 2018. © 2017 Wiley Periodicals, Inc.     

RET and NRG1 interplay in Hirschsprung disease

Hirschsprung disease (HSCR, aganglionic megacolon) is a complex genetic disorder of the enteric nervous system (ENS) characterized by the absence of enteric neurons along a variable length of the intestine. While rare variants (RVs) in the coding sequence (CDS) of several genes involved in ENS development lead to disease, the association of common variants (CVs) with HSCR has only been reported for RET (the major HSCR gene) and NRG1. Importantly, RVs in the CDS of these two genes are also associated with the disorder. To assess independent and joint effects between the different types of RET and NRG1 variants identified in HSCR patients, we used 254 Chinese sporadic HSCR patients and 143 ethnically matched controls for whom the RET and/or NRG1 variants genotypes (rare and common) were available. Four genetic risk factors were defined and interaction effects were modeled using conditional logistic regression analyses and pair-wise Kendall correlations. Our analysis revealed a joint effect of RET CVs with RET RVs, NRG1 CVs or NRG1 RVs. To assess whether the genetic interaction translated into functional interaction, mouse neural crest cells (NCCs; enteric neuron precursors) isolated from embryonic guts were treated with NRG1 (ErbB2 ligand) or/and GDNF (Ret ligand) and monitored during the subsequent neural differentiation process. Nrg1 inhibited the Gdnf-induced neuronal differentiation and Gdnf negatively regulated Nrg1-signaling by down-regulating the expression of its receptor, ErbB2. This preliminary data suggest that the balance neurogenesis/gliogenesis is critical for ENS development.     

Mutational Spectrum of Semaphorin 3A and Semaphorin 3D Genes in Spanish Hirschsprung patients

Hirschsprung disease (HSCR, OMIM 142623) is a developmental disorder characterized by the absence of ganglion cells along variable lengths of the distal gastrointestinal tract, which results in tonic contraction of the aganglionic colon segment and functional intestinal obstruction. The RET proto-oncogene is the major gene associated to HSCR with differential contributions of its rare and common, coding and noncoding mutations to the multifactorial nature of this pathology. In addition, many other genes have been described to be associated with this pathology, including the semaphorins class III genes SEMA3A (7p12.1) and SEMA3D (7q21.11) through SNP array analyses and by next-generation sequencing technologies. Semaphorins are guidance cues for developing neurons implicated in the axonal projections and in the determination of the migratory pathway for neural-crest derived neural precursors during enteric nervous system development. In addition, it has been described that increased SEMA3A expression may be a risk factor for HSCR through the upregulation of the gene in the aganglionic smooth muscle layer of the colon in HSCR patients. Here we present the results of a comprehensive analysis of SEMA3A and SEMA3D in a series of 200 Spanish HSCR patients by the mutational screening of its coding sequence, which has led to find a number of potentially deleterious variants. RET mutations have been also detected in some of those patients carrying SEMAs variants. We have evaluated the A131T-SEMA3A, S598G-SEMA3A and E198K-SEMA3D mutations using colon tissue sections of these patients by immunohistochemistry. All mutants presented increased protein expression in smooth muscle layer of ganglionic segments. Moreover, A131T-SEMA3A also maintained higher protein levels in the aganglionic muscle layers. These findings strongly suggest that these mutants have a pathogenic effect on the disease. Furthermore, because of their coexistence with RET mutations, our data substantiate the additive genetic model proposed for this rare disorder and further support the association of SEMAs genes with HSCR.     

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.     

Integrated Model of De Novo and Inherited Genetic Variants Yields Greater Power to Identify Risk Genes

De novo mutations affect risk for many diseases and disorders, especially those with early-onset. An example is autism spectrum disorders (ASD). Four recent whole-exome sequencing (WES) studies of ASD families revealed a handful of novel risk genes, based on independent de novo loss-of-function (LoF) mutations falling in the same gene, and found that de novo LoF mutations occurred at a twofold higher rate than expected by chance. However successful these studies were, they used only a small fraction of the data, excluding other types of de novo mutations and inherited rare variants. Moreover, such analyses cannot readily incorporate data from case-control studies. An important research challenge in gene discovery, therefore, is to develop statistical methods that accommodate a broader class of rare variation. We develop methods that can incorporate WES data regarding de novo mutations, inherited variants present, and variants identified within cases and controls. TADA, for Transmission And De novo Association, integrates these data by a gene-based likelihood model involving parameters for allele frequencies and gene-specific penetrances. Inference is based on a Hierarchical Bayes strategy that borrows information across all genes to infer parameters that would be difficult to estimate for individual genes. In addition to theoretical development we validated TADA using realistic simulations mimicking rare, large-effect mutations affecting risk for ASD and show it has dramatically better power than other common methods of analysis. Thus TADA''s integration of various kinds of WES data can be a highly effective means of identifying novel risk genes. Indeed, application of TADA to WES data from subjects with ASD and their families, as well as from a study of ASD subjects and controls, revealed several novel and promising ASD candidate genes with strong statistical support.     

Copy number variants in candidate genes are genetic modifiers of Hirschsprung disease

Hirschsprung disease (HSCR) is a neurocristopathy characterized by absence of intramural ganglion cells along variable lengths of the gastrointestinal tract. The HSCR phenotype is highly variable with respect to gender, length of aganglionosis, familiality and the presence of additional anomalies. By molecular genetic analysis, a minimum of 11 neuro-developmental genes (RET, GDNF, NRTN, SOX10, EDNRB, EDN3, ECE1, ZFHX1B, PHOX2B, KIAA1279, TCF4) are known to harbor rare, high-penetrance mutations that confer a large risk to the bearer. In addition, two other genes (RET, NRG1) harbor common, low-penetrance polymorphisms that contribute only partially to risk and can act as genetic modifiers. To broaden this search, we examined whether a set of 67 proven and candidate HSCR genes harbored additional modifier alleles. In this pilot study, we utilized a custom-designed array CGH with ～33,000 test probes at an average resolution of ～185 bp to detect gene-sized or smaller copy number variants (CNVs) within these 67 genes in 18 heterogeneous HSCR patients. Using stringent criteria, we identified CNVs at three loci (MAPK10, ZFHX1B, SOX2) that are novel, involve regulatory and coding sequences of neuro-developmental genes, and show association with HSCR in combination with other congenital anomalies. Additional CNVs are observed under relaxed criteria. Our research suggests a role for CNVs in HSCR and, importantly, emphasizes the role of variation in regulatory sequences. A much larger study will be necessary both for replication and for identifying the full spectrum of small CNV effects.     

Genetik des Prostatakarzinoms

Prostate cancer is the most frequent malignancy in males and its etiology is strongly influenced by genetic factors. Nevertheless, no mutated genes which could be used for diagnosis have been identified in a major proportion of familial cases. Three genes with germline mutations have been identified after linkage analysis (ELAC2, RNASEL, MSR1), but these mutations are very rare and their penetrance is not well defined. The association of most genes with genetic variants is weak, and only BRCA2/familial breast cancer is of clinical relevance. As a consequence of the extreme genetic heterogeneity, diagnostic tools are not available and genetic counseling has to rely on risk estimates from pedigree data in which a single affected first degree relative indicates a relevant risk.     

Genetics and aetiology of Pagetic disorders of bone

Paget’s disease of bone (PDB) is a late-onset disorder characterised by focal areas of increased bone turnover containing enlarged hyperactive osteoclasts. The disease has a strong genetic predisposition and mutations in SQSTM1 have been associated with familial and sporadic disease in up to 40% of cases. Additional genetic loci have been associated in other cases, but genes are yet to be identified. Earlier-onset conditions with similar bone pathology (familial expansile osteolysis, expansile skeletal hyperphosphatasia and early-onset PDB) are caused by mutations in TNFRSF11A (RANK). The syndrome of inclusion body myositis, Paget’s disease and frontotemporal dementia is caused by mutations in VCP. Despite the increased knowledge about genes involved in PDB and related disorders, the etiology of the diseases remains puzzling. Presence of inclusion bodies appears to link Pagetic diseases mechanistically to diseases associated with presence of misfolded proteins or abnormalities in the ubiquitin-proteasomal, or autophagy pathways. Juvenile PDB, caused by osteoprotegerin deficiency, appears mechanistically distinct from the other Pagetic diseases. This review will discuss evidence from recent studies, including new animal models for Pagetic diseases.     

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.     

Deep Resequencing of GWAS Loci Identifies Rare Variants in CARD9, IL23R and RNF186 That Are Associated with Ulcerative Colitis

Genome-wide association studies and follow-up meta-analyses in Crohn''s disease (CD) and ulcerative colitis (UC) have recently identified 163 disease-associated loci that meet genome-wide significance for these two inflammatory bowel diseases (IBD). These discoveries have already had a tremendous impact on our understanding of the genetic architecture of these diseases and have directed functional studies that have revealed some of the biological functions that are important to IBD (e.g. autophagy). Nonetheless, these loci can only explain a small proportion of disease variance (∼14% in CD and 7.5% in UC), suggesting that not only are additional loci to be found but that the known loci may contain high effect rare risk variants that have gone undetected by GWAS. To test this, we have used a targeted sequencing approach in 200 UC cases and 150 healthy controls (HC), all of French Canadian descent, to study 55 genes in regions associated with UC. We performed follow-up genotyping of 42 rare non-synonymous variants in independent case-control cohorts (totaling 14,435 UC cases and 20,204 HC). Our results confirmed significant association to rare non-synonymous coding variants in both IL23R and CARD9, previously identified from sequencing of CD loci, as well as identified a novel association in RNF186. With the exception of CARD9 (OR = 0.39), the rare non-synonymous variants identified were of moderate effect (OR = 1.49 for RNF186 and OR = 0.79 for IL23R). RNF186 encodes a protein with a RING domain having predicted E3 ubiquitin-protein ligase activity and two transmembrane domains. Importantly, the disease-coding variant is located in the ubiquitin ligase domain. Finally, our results suggest that rare variants in genes identified by genome-wide association in UC are unlikely to contribute significantly to the overall variance for the disease. Rather, these are expected to help focus functional studies of the corresponding disease loci.     

Genome-wide searching of rare genetic variants in WTCCC data

Although they have demonstrated success in searching for common variants for complex diseases, genome-wide association (GWA) studies are less successful in detecting rare genetic variants because of the poor statistical power of most of current methods. We developed a two-stage method that can apply to GWA studies for detecting rare variants. Here we report the results of applying this two-stage method to the Wellcome Trust Case Control Consortium (WTCCC) dataset that include seven complex diseases: bipolar disorder, cardiovascular disease, hypertension (HT), rheumatoid arthritis, Crohn’s disease, type 1 diabetes and type 2 diabetes (T2D). We identified 24 genes or regions that reach genome wide significance. Eight of them are novel and were not reported in the WTCCC study. The cumulative risk (or protective) haplotype frequency for each of the 8 genes or regions is small, being at most 11%. For each of the novel genes, the risk (or protective) haplotype set cannot be tagged by the common SNPs available in chips (r ² < 0.32). The gene identified in HT was further replicated in the Framingham Heart Study, and is also significantly associated with T2D. Our analysis suggests that searching for rare genetic variants is feasible in current GWA studies and candidate gene studies, and the results can severe as guides to future resequencing studies to identify the underlying rare functional variants.     

Rare variants regulate expression of nearby individual genes in multiple tissues

The rapid decrease in sequencing cost has enabled genetic studies to discover rare variants associated with complex diseases and traits. Once this association is identified, the next step is to understand the genetic mechanism of rare variants on how the variants influence diseases. Similar to the hypothesis of common variants, rare variants may affect diseases by regulating gene expression, and recently, several studies have identified the effects of rare variants on gene expression using heritability and expression outlier analyses. However, identifying individual genes whose expression is regulated by rare variants has been challenging due to the relatively small sample size of expression quantitative trait loci studies and statistical approaches not optimized to detect the effects of rare variants. In this study, we analyze whole-genome sequencing and RNA-seq data of 681 European individuals collected for the Genotype-Tissue Expression (GTEx) project (v8) to identify individual genes in 49 human tissues whose expression is regulated by rare variants. To improve statistical power, we develop an approach based on a likelihood ratio test that combines effects of multiple rare variants in a nonlinear manner and has higher power than previous approaches. Using GTEx data, we identify many genes regulated by rare variants, and some of them are only regulated by rare variants and not by common variants. We also find that genes regulated by rare variants are enriched for expression outliers and disease-causing genes. These results suggest the regulatory effects of rare variants, which would be important in interpreting associations of rare variants with complex traits.     

Targeted Sequencing of the Mitochondrial Genome of Women at High Risk of Breast Cancer without Detectable Mutations in BRCA1/2

Breast Cancer is a complex multifactorial disease for which high-penetrance mutations have been identified. Approaches used to date have identified genomic features explaining about 50% of breast cancer heritability. A number of low- to medium penetrance alleles (per-allele odds ratio < 1.5 and 4.0, respectively) have been identified, suggesting that the remaining heritability is likely to be explained by the cumulative effect of such alleles and/or by rare high-penetrance alleles. Relatively few studies have specifically explored the mitochondrial genome for variants potentially implicated in breast cancer risk. For these reasons, we propose an exploration of the variability of the mitochondrial genome in individuals diagnosed with breast cancer, having a positive breast cancer family history but testing negative for BRCA1/2 pathogenic mutations. We sequenced the mitochondrial genome of 436 index breast cancer cases from the GENESIS study. As expected, no pathogenic genomic pattern common to the 436 women included in our study was observed. The mitochondrial genes MT-ATP6 and MT-CYB were observed to carry the highest number of variants in the study. The proteins encoded by these genes are involved in the structure of the mitochondrial respiration chain, and variants in these genes may impact reactive oxygen species production contributing to carcinogenesis. More functional and epidemiological studies are needed to further investigate to what extent variants identified may influence familial breast cancer risk.     

Most rare missense alleles are deleterious in humans: implications for complex disease and association studies

The accumulation of mildly deleterious missense mutations in individual human genomes has been proposed to be a genetic basis for complex diseases. The plausibility of this hypothesis depends on quantitative estimates of the prevalence of mildly deleterious de novo mutations and polymorphic variants in humans and on the intensity of selective pressure against them. We combined analysis of mutations causing human Mendelian diseases, of human-chimpanzee divergence, and of systematic data on human genetic variation and found that ~20% of new missense mutations in humans result in a loss of function, whereas ~27% are effectively neutral. Thus, the remaining 53% of new missense mutations have mildly deleterious effects. These mutations give rise to many low-frequency deleterious allelic variants in the human population, as is evident from a new data set of 37 genes sequenced in >1,500 individual human chromosomes. Surprisingly, up to 70% of low-frequency missense alleles are mildly deleterious and are associated with a heterozygous fitness loss in the range 0.001-0.003. Thus, the low allele frequency of an amino acid variant can, by itself, serve as a predictor of its functional significance. Several recent studies have reported a significant excess of rare missense variants in candidate genes or pathways in individuals with extreme values of quantitative phenotypes. These studies would be unlikely to yield results if most rare variants were neutral or if rare variants were not a significant contributor to the genetic component of phenotypic inheritance. Our results provide a justification for these types of candidate-gene (pathway) association studies and imply that mutation-selection balance may be a feasible evolutionary mechanism underlying some common diseases.     

Identification of Novel Candidate Genes for Early-Onset Colorectal Cancer Susceptibility

Approximately 25–30% of colorectal cancer (CRC) cases are expected to result from a genetic predisposition, but in only 5–10% of these cases highly penetrant germline mutations are found. The remaining CRC heritability is still unexplained, and may be caused by a hitherto-undefined set of rare variants with a moderately penetrant risk. Here we aimed to identify novel risk factors for early-onset CRC using whole-exome sequencing, which was performed on a cohort of CRC individuals (n = 55) with a disease onset before 45 years of age. We searched for genes that were recurrently affected by rare variants (minor allele frequency ≤0.001) with potentially damaging effects and, subsequently, re-sequenced the candidate genes in a replication cohort of 174 early-onset or familial CRC individuals. Two functionally relevant genes with low frequency variants with potentially damaging effects, PTPN12 and LRP6, were found in at least three individuals. The protein tyrosine phosphatase PTP-PEST, encoded by PTPN12, is a regulator of cell motility and LRP6 is a component of the WNT-FZD-LRP5-LRP6 complex that triggers WNT signaling. All variants in LRP6 were identified in individuals with an extremely early-onset of the disease (≤30 years of age), and two of the three variants showed increased WNT signaling activity in vitro. In conclusion, we present PTPN12 and LRP6 as novel candidates contributing to the heterogeneous susceptibility to CRC.     

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 genetics of Hirschsprung disease: a model of multigenic neurocristopathy

Hirschsprung's disease (HSCR, aganglionic megacolon) is a frequent congenital malformation regarded as a multigenic neurocristopathy. Three susceptibility genes have been recently identified in HSCR, namely the RET proto-oncogene, the endothelin B receptor (EDNRB) gene, and the endothelin 3 (EDN3) gene. RET gene mutations were found in significant proportions of familial (50%) and sporadic (15-20%) HSCR, while homozygosity for EDNRB or EDN3 mutations accounted for the rare HSCR-Waardenburg syndrome (WS) association. More recently, heterozygous EDNRB an EDN3 missense mutations have been reported in isolated HSCR patients. Some of these results were obtained after the identification of mouse genes whose natural or site-directed mutations resulted in megacolon and coat color spotting. There is also conclusive evidence for the involvement of other independent loci in HSCR. In particular, the recent identification of neurotrophic factors acting as RET ligands (GDNF and Neurturin) provide additional candidate genes for HSCR. The dissection of the genetic etiology of HSCR disease may then provide a unique opportunity to distinguish between a polygenic and a genetically heterogeneous disease, thereby helping to understand other complex disorders and congenital malformations hitherto considered as multifactorial in origin. Finally, the study of the molecular bases of HSCR is also a step towards the understanding of developmental genetics of the enteric nervous system giving support to the role of the tyrosine kinase and endothelin-signaling pathways in the development of neural crest-derived enteric neurons in human.     

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 HSCRSince 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.     

The genetic architecture of autism spectrum disorders(ASDs) and the potential importance of common regulatory genetic variants

Currently, there is great interest in identifying genetic variants that contribute to the risk of developing autism spectrum disorders(ASDs), due in part to recent increases in the frequency of diagnosis of these disorders worldwide. While there is nearly universal agreement that ASDs are complex diseases, with multiple genetic and environmental contributing factors, there is less agreement concerning the relative importance of common vs rare genetic variants in ASD liability. Recent observations that rare mutations and copy number variants(CNVs) are frequently associated with ASDs, combined with reduced fecundity of individuals with these disorders, has led to the hypothesis that ASDs are caused primarily by de novo or rare genetic mutations. Based on this model, large-scale whole-genome DNA sequencing has been proposed as the most appropriate method for discovering ASD liability genes. While this approach will undoubtedly identify many novel candidate genes and produce important new insights concerning the genetic causes of these disorders, a full accounting of the genetics of ASDs will be incomplete absent an understanding of the contributions of common regulatory variants, which are likely to influence ASD liability by modifying the effects of rare variants or, by assuming unfavorable combinations, directly produce these disorders. Because it is not yet possible to identify regulatory genetic variants by examination of DNA sequences alone, their identification will require experimentation. In this essay, I discuss these issues and describe the advantages of measurements of allelic expression imbalance(AEI) of m RNA expression for identifying cis-acting regulatory variants that contribute to ASDs.