Understanding the complex genetic architecture connecting rheumatoid arthritis, osteoporosis and inflammation: discovering causal pathways.

Rheumatoid arthritis (RA) and osteoporosis (OP) are two comorbid complex inflammatory conditions with evidence of shared genetic background and causal relationships. We aimed to clarify the genetic architecture underlying RA and various OP phenotypes while additionally considering an inflammatory component, C-reactive protein (CRP). Genome-wide association study summary statistics were acquired from the GEnetic Factors for OSteoporosis Consortium, Cohorts for Heart and Aging Research Consortium and UK Biobank. Mendelian randomization (MR) was used to detect the presence of causal relationships. Colocalization analysis was performed to determine shared genetic variants between CRP and OP phenotypes. Analysis of pleiotropy between traits owing to shared causal single nucleotide polymorphisms (SNPs) was performed using PL eiotropic A nalysis under CO mposite null hypothesis (PLACO). MR analysis was suggestive of horizontal pleiotropy between RA and OP traits. RA was a significant causal risk factor for CRP (ß = 0.027, 95% confidence interval = 0.016-0.038). There was no evidence of CRP→OP causal relationship, but horizontal pleiotropy was apparent. Colocalization established shared genomic regions between CRP and OP, including GCKR and SERPINA1 genes. Pleiotropy arising from shared causal SNPs revealed through the colocalization analysis was all confirmed by PLACO. These genes were found to be involved in the same molecular function 'protein binding' (GO:0005515) associated with RA, OP and CRP. We identified three major components explaining the epidemiological relationship among RA, OP and inflammation: (1) Pleiotropy explains a portion of the shared genetic relationship between RA and OP, albeit polygenically; (2) RA contributes to CRP elevation and (3) CRP, which is influenced by RA, demonstrated pleiotropy with OP.

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.

Relationships between Circulating Biomarkers and Body Composition Parameters in Patients with Metabolic Syndrome: A Community-Based Study.

Metabolic syndrome (MetS) is a complex disease involving multiple physiological, biochemical, and metabolic abnormalities. The search for reliable biomarkers may help to better elucidate its pathogenesis and develop new preventive and therapeutic strategies. In the present population-based study, we looked for biomarkers of MetS among obesity- and inflammation-related circulating factors and body composition parameters in 1079 individuals (with age range between 18 and 80) belonging to an ethnically homogeneous population. Plasma levels of soluble markers were measured by using ELISA. Body composition parameters were assessed using bioimpedance analysis (BIA). Statistical analysis, including mixed-effects regression, with MetS as a dependent variable, revealed that the most significant independent variables were mainly adipose tissue-related phenotypes, including fat mass/weight (FM/WT) [OR (95% CI)], 2.77 (2.01-3.81); leptin/adiponectin ratio (L/A ratio), 1.50 (1.23-1.83); growth and differentiation factor 15 (GDF-15) levels, 1.32 (1.08-1.62); inflammatory markers, specifically monocyte to high-density lipoprotein cholesterol ratio (MHR), 2.53 (2.00-3.15), and a few others. Additive Bayesian network modeling suggests that age, sex, MHR, and FM/WT are directly associated with MetS and probably affect its manifestation. Additionally, MetS may be causing the GDF-15 and L/A ratio. Our novel findings suggest the existence of complex, age-related, and possibly hierarchical relationships between MetS and factors associated with obesity.

Glycoprotein Acetyls Is a Novel Biomarker Predicting Cardiovascular Complications in Rheumatoid Arthritis.

The relationship between rheumatoid arthritis (RA) and early onset atherosclerosis is well depicted, each with an important inflammatory component. Glycoprotein acetyls (GlycA), a novel biomarker of inflammation, may play a role in the manifestation of these two inflammatory conditions. The present study examined a potential mediating role of GlycA within the RA-atherosclerosis relationship to determine whether it accounts for the excess risk of cardiovascular disease over that posed by lipid risk factors. The UK Biobank dataset was acquired to establish associations among RA, atherosclerosis, GlycA, and major lipid factors: total cholesterol (TC), high- and low-density lipoprotein (HDL, LDL) cholesterol, and triglycerides (TGs). Genome-wide association study summary statistics were collected from various resources to perform genetic analyses. Causality among variables was tested using Mendelian Randomization (MR) analysis. Genes of interest were identified using colocalization analysis and gene enrichment analysis. MR results appeared to indicate that the genetic relationship between GlycA and RA and also between RA and atherosclerosis was explained by horizontal pleiotropy (p-value = 0.001 and <0.001, respectively), while GlycA may causally predict atherosclerosis (p-value = 0.017). Colocalization analysis revealed several functionally relevant genes shared between GlycA and all the variables assessed. Two loci were apparent in all relationships tested and included the HLA region as well as SLC22A1. GlycA appears to mediate the RA-atherosclerosis relationship through several possible pathways. GlycA, although pleiotropically related to RA, appears to causally predict atherosclerosis. Thus, GlycA is suggested as a significant factor in the etiology of atherosclerosis development in RA.

Genetic and environmental correlational structure among metabolic syndrome endophenotypes.

Metabolic syndrome (MetS) is diagnosed by the presence of high scores on three or more metabolic traits, including systolic and diastolic blood pressure (SBP, DBP), glucose and insulin levels, cholesterol and triglyceride (TG) levels, and central obesity. A diagnosis of MetS is associated with increased risk of cardiovascular disease and type 2 diabetes. The components of MetS have long been demonstrated to have substantial genetic components, but their genetic overlap is less well understood. The present paper takes a multi-prong approach to examining the extent of this genetic overlap. This includes the quantitative genetic and additive Bayesian network modeling of the large TwinsUK project and examination of the results of genome-wide association study (GWAS) of UK Biobank data through use of LD score regression and examination of the number of genes and pathways identified in the GWASes which overlap across MetS traits. Results demonstrate a modest genetic overlap, and the genetic correlations obtained from TwinsUK and UK Biobank are nearly identical. However, these correlations imply more genetic dissimilarity than similarity. Furthermore, examination of the extent of overlap in significant GWAS hits, both at the gene and pathway level, again demonstrates only modest but significant genetic overlap. This lends support to the idea that in clinical treatment of MetS, treating each of the components individually may be an important way to address MetS.

Revealing antibiotic cross-resistance patterns in hospitalized patients through Bayesian network modelling.

OBJECTIVES: Microbial resistance exhibits dependency patterns between different antibiotics, termed cross-resistance and collateral sensitivity. These patterns differ between experimental and clinical settings. It is unclear whether the differences result from biological reasons or from confounding, biasing results found in clinical settings. We set out to elucidate the underlying dependency patterns between resistance to different antibiotics from clinical data, while accounting for patient characteristics and previous antibiotic usage. METHODS: Additive Bayesian network modelling was employed to simultaneously estimate relationships between variables in a dataset of bacterial cultures derived from hospitalized patients and tested for resistance to multiple antibiotics. Data contained resistance results, patient demographics and previous antibiotic usage, for five bacterial species: Escherichia coli (n = 1054), Klebsiella pneumoniae (n = 664), Pseudomonas aeruginosa (n = 571), CoNS (n = 495) and Proteus mirabilis (n = 415). RESULTS: All links between resistance to the various antibiotics were positive. Multiple direct links between resistance of antibiotics from different classes were observed across bacterial species. For example, resistance to gentamicin in E. coli was directly linked with resistance to ciprofloxacin (OR = 8.39, 95% credible interval 5.58-13.30) and sulfamethoxazole/trimethoprim (OR = 2.95, 95% credible interval 1.97-4.51). In addition, resistance to various antibiotics was directly linked with previous antibiotic usage. CONCLUSIONS: Robust relationships among resistance to antibiotics belonging to different classes, as well as resistance being linked to having taken antibiotics of a different class, exist even when taking into account multiple covariate dependencies. These relationships could help inform choices of antibiotic treatment in clinical settings.

De novo mutations in Caudal Type Homeo Box transcription Factor 2 (CDX2) in patients with persistent cloaca.

The cloaca is an embryonic cavity that is divided into the urogenital sinus and rectum upon differentiation of the cloacal epithelium triggered by tissue-specific transcription factors including CDX2. Defective differentiation leads to persistent cloaca in humans (PC), a phenotype recapitulated in Cdx2 mutant mice. PC is linked to hypo/hyper-vitaminosis A. Although no gene has ever been identified, there is a strong evidence for a genetic contribution to PC. We applied whole-exome sequencing and copy-number-variants analyses to 21 PC patients and their unaffected parents. The damaging p.Cys132* and p.Arg237His de novo CDX2 variants were identified in two patients. These variants altered the expression of CYP26A1, a direct CDX2 target encoding the major retinoic acid (RA)-degrading enzyme. Other RA genes, including the RA-receptor alpha, were also mutated. Genes governing the development of cloaca-derived structures were recurrently mutated and over-represented in the basement-membrane components set (q-value < 1.65 × 10-6). Joint analysis of the patients' profile highlighted the extracellular matrix-receptor interaction pathway (MsigDBID: M7098, FDR: q-value < 7.16 × 10-9). This is the first evidence that PC is genetic, with genes involved in the RA metabolism at the lead. Given the CDX2 de novo variants and the role of RA, our observations could potentiate preventive measures. For the first time, a gene recapitulating PC in mouse models is found mutated in humans.

The Genes We Inherit and Those We Don't: Maternal Genetic Nurture and Child BMI Trajectories.

Recently, methods have been introduced using polygenic scores (PGS) to estimate the effects of genetic nurture, the environmentally-mediated effects of parental genotypes on the phenotype of their child above and beyond the effects of the alleles which are transmitted to the child. We introduce a simplified model for estimating genetic nurture effects and show, through simulation and analytical derivation, that our method provides unbiased estimates and offers an increase in power to detect genetic nurture of up to 1/3 greater than that of previous methods. Subsequently, we apply this method to data from the Avon Longitudinal Study of Parents and Children to estimate the effects of maternal genetic nurture on childhood body mass index (BMI) trajectories. Through mixed modeling, we observe a statistically significant age-dependent effect of maternal PGS on child BMI, such that the influence of maternal genetic nurture appears to increase throughout development.

Statistical Power and the Classical Twin Design.

Dr Nick Martin has made enormous contributions to the field of behavior genetics over the past 50 years. Of his many seminal papers that have had a profound impact, we focus on his early work on the power of twin studies. He was among the first to recognize the importance of sample size calculation before conducting a study to ensure sufficient power to detect the effects of interest. The elegant approach he developed, based on the noncentral chi-squared distribution, has been adopted by subsequent researchers for other genetic study designs, and today remains a standard tool for power calculations in structural equation modeling and other areas of statistical analysis. The present brief article discusses the main aspects of his seminal paper, and how it led to subsequent developments, by him and others, as the field of behavior genetics evolved into the present era.

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.

Genetic overlap between epilepsy and schizophrenia: Evidence from cross phenotype analysis in Hong Kong Chinese population.

Epilepsy and schizophrenia are common and typical neurological or mental illness respectively, and sometimes they comorbid in the same patients, however the underlying genetic relationship between the two brain diseases is still not fully understood. To investigate the possible genetic contribution to their comorbidity, we performed polygenic risk score (PRS) analyses and genetic correlation estimation so as to identify the overall genetic overlap between the two diseases. The global schizophrenia PRS is strongly associated with schizophrenia phenotype in Hong Kong population (odds ratio = 1.7, p = 2.26E-16), and focal epilepsy PRS is moderately associated with epilepsy phenotype in Hong Kong population (odds ratio = 1.14, p = 0.013). However the disease-specific PRS can only predict its own well-matched phenotype but not the other ones (p > 0.05). This pattern is further supported by non-significant pairwise genetic correlation and insufficient statistical power for PRS association from the cross-phenotype analyses. Our study reveals there's limited shared genetic aetiology between schizophrenia and epilepsy, and thus supports a model of shared environmental factors to explain the comorbidity between the two phenotypes.

Exome-chip association analysis reveals an Asian-specific missense variant in PAX4 associated with type 2 diabetes in Chinese individuals.

AIMS/HYPOTHESIS: Genome-wide association studies (GWASs) have identified many common type 2 diabetes-associated variants, mostly at the intronic or intergenic regions. Recent advancements of exome-array genotyping platforms have opened up a novel means for detecting the associations of low-frequency or rare coding variants with type 2 diabetes. We conducted an exomechip association analysis to identify additional type 2 diabetes susceptibility variants in the Chinese population. METHODS: An exome-chip association study was conducted by genotyping 5640 Chinese individuals from Hong Kong, using a custom designed exome array, the Asian Exomechip. Single variant association analysis was conducted on 77,468 single nucleotide polymorphisms (SNPs). Fifteen SNPs were subsequently genotyped for replication analysis in an independent Chinese cohort comprising 12,362 individuals from Guangzhou. A combined analysis involving 7189 cases and 10,813 controls was performed. RESULTS: In the discovery stage, an Asian-specific coding variant rs2233580 (p.Arg192His) in PAX4, and two variants at the known loci, CDKN2B-AS1 and KCNQ1, were significantly associated with type 2 diabetes with exome-wide significance (p discovery < 6.45 × 10-7). The risk allele (T) of PAX4 rs2233580 was associated with a younger age at diabetes diagnosis. This variant was replicated in an independent cohort and demonstrated a stronger association that reached genome-wide significance (p meta-analysis [p meta] = 3.74 × 10-15) in the combined analysis. CONCLUSIONS/INTERPRETATION: We identified the association of a PAX4 Asian-specific missense variant rs2233580 with type 2 diabetes in an exome-chip association analysis, supporting the involvement of PAX4 in the pathogenesis of type 2 diabetes. Our findings suggest PAX4 is a possible effector gene of the 7q32 locus, previously identified from GWAS in Asians.

Meta-analysis of GWAS on two Chinese populations followed by replication identifies novel genetic variants on the X chromosome associated with systemic lupus erythematosus.

Systemic lupus erythematosus (SLE) is a prototypic autoimmune disease that affects mainly females. What role the X chromosome plays in the disease has always been an intriguing question. In this study, we examined the genetic variants on the X chromosome through meta-analysis of two genome-wide association studies (GWAS) on SLE on Chinese Han populations. Prominent association signals from the meta-analysis were replicated in 4 additional Asian cohorts, with a total of 5373 cases and 9166 matched controls. We identified a novel variant in PRPS2 on Xp22.3 as associated with SLE with genome-wide significance (rs7062536, OR = 0.84, P = 1.00E-08). Association of the L1CAM-MECP2 region with SLE was reported previously. In this study, we identified independent contributors in this region in NAA10 (rs2071128, OR = 0.81, P = 2.19E-13) and TMEM187 (rs17422, OR = 0.75, P = 1.47E-15), in addition to replicating the association from IRAK1-MECP2 region (rs1059702, OR = 0.71, P = 2.40E-18) in Asian cohorts. The X-linked susceptibility variants showed higher effect size in males than that in females, similar to results from a genome-wide survey of associated SNPs on the autosomes. These results suggest that susceptibility genes identified on the X chromosome, while contributing to disease predisposition, might not contribute significantly to the female predominance of this prototype autoimmune disease.

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

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

Genome-wide search followed by replication reveals genetic interaction of CD80 and ALOX5AP associated with systemic lupus erythematosus in Asian populations.

OBJECTIVES: Genetic interaction has been considered as a hallmark of the genetic architecture of systemic lupus erythematosus (SLE). Based on two independent genome-wide association studies (GWAS) on Chinese populations, we performed a genome-wide search for genetic interactions contributing to SLE susceptibility. METHODS: The study involved a total of 1 659 cases and 3 398 controls in the discovery stage and 2 612 cases and 3 441 controls in three cohorts for replication. Logistic regression and multifactor dimensionality reduction were used to search for genetic interaction. RESULTS: Interaction of CD80 (rs2222631) and ALOX5AP (rs12876893) was found to be significantly associated with SLE (OR_int=1.16, P_int_all=7.7E-04 at false discovery rate<0.05). Single nuclear polymorphism rs2222631 was found associated with SLE with genome-wide significance (P_all=4.5E-08, OR=0.86) and is independent of rs6804441 in CD80, whose association was reported previously. Significant correlation was observed between expression of these two genes in healthy controls and SLE cases, together with differential expression of these genes between cases and controls, observed from individuals from the Hong Kong cohort. Genetic interactions between BLK (rs13277113) and DDX6 (rs4639966), and between TNFSF4 (rs844648) and PXK (rs6445975) were also observed in both GWAS data sets. CONCLUSIONS: Our study represents the first genome-wide evaluation of epistasis interactions on SLE and the findings suggest interactions and independent variants may help partially explain missing heritability for complex diseases.

Sacral agenesis: a pilot whole exome sequencing and copy number study.

BACKGROUND: Caudal regression syndrome (CRS) or sacral agenesis is a rare congenital disorder characterized by a constellation of congenital caudal anomalies affecting the caudal spine and spinal cord, the hindgut, the urogenital system, and the lower limbs. CRS is a complex condition, attributed to an abnormal development of the caudal mesoderm, likely caused by the effect of interacting genetic and environmental factors. A well-known risk factor is maternal type 1 diabetes. METHOD: Whole exome sequencing and copy number variation (CNV) analyses were conducted on 4 Caucasian trios to identify de novo and inherited rare mutations. RESULTS: In this pilot study, exome sequencing and copy number variation (CNV) analyses implicate a number of candidate genes, including SPTBN5, MORN1, ZNF330, CLTCL1 and PDZD2. De novo mutations were found in SPTBN5, MORN1 and ZNF330 and inherited predicted damaging mutations in PDZD2 (homozygous) and CLTCL1 (compound heterozygous). Importantly, predicted damaging mutations in PTEN (heterozygous), in its direct regulator GLTSCR2 (compound heterozygous) and in VANGL1 (heterozygous) were identified. These genes had previously been linked with the CRS phenotype. Two CNV deletions, one de novo (chr3q13.13) and one homozygous (chr8p23.2), were detected in one of our CRS patients. These deletions overlapped with CNVs previously reported in patients with similar phenotype. CONCLUSION: Despite the genetic diversity and the complexity of the phenotype, this pilot study identified genetic features common across CRS patients.

Longitudinal heritability of childhood aggression.

The genetic and environmental contributions to the variation and longitudinal stability in childhood aggressive behavior were assessed in two large twin cohorts, the Netherlands Twin Register (NTR), and the Twins Early Development Study (TEDS; United Kingdom). In NTR, maternal ratings on aggression from the Child Behavior Checklist (CBCL) were available for 10,765 twin pairs at age 7, for 8,557 twin pairs at age 9/10, and for 7,176 twin pairs at age 12. In TEDS, parental ratings of conduct disorder from the Strength and Difficulty Questionnaire (SDQ) were available for 6,897 twin pairs at age 7, for 3,028 twin pairs at age 9 and for 5,716 twin pairs at age 12. In both studies, stability and heritability of aggressive behavioral problems was high. Heritability was on average somewhat, but significantly, lower in TEDS (around 60%) than in NTR (between 50% and 80%) and sex differences were slightly larger in the NTR sample. In both studies, the influence of shared environment was similar: in boys shared environment explained around 20% of the variation in aggression across all ages while in girls its influence was absent around age 7 and only came into play at later ages. Longitudinal genetic correlations were the main reason for stability of aggressive behavior. Individual differences in CBCL-Aggressive Behavior and SDQ-Conduct disorder throughout childhood are driven by a comparable but significantly different genetic architecture. © 2016 Wiley Periodicals, Inc.

Genome-wide copy number variation study in anorectal malformations.

Anorectal malformations (ARMs, congenital obstruction of the anal opening) are among the most common birth defects requiring surgical treatment (2-5/10 000 live-births) and carry significant chronic morbidity. ARMs present either as isolated or as part of the phenotypic spectrum of some chromosomal abnormalities or monogenic syndromes. The etiology is unknown. To assess the genetic contribution to ARMs, we investigated single-nucleotide polymorphisms and copy number variations (CNVs) at genome-wide scale. A total of 363 Han Chinese sporadic ARM patients and 4006 Han Chinese controls were included. Overall, we detected a 1.3-fold significant excess of rare CNVs in patients. Stratification of patients by presence/absence of other congenital anomalies showed that while syndromic ARM patients carried significantly longer rare duplications than controls (P = 0.049), non-syndromic patients were enriched with both rare deletions and duplications when compared with controls (P = 0.00031). Twelve chromosomal aberrations and 114 rare CNVs were observed in patients but not in 868 controls nor 11 943 healthy individuals from the Database of Genomic Variants. Importantly, these aberrations were observed in isolated ARM patients. Gene-based analysis revealed 79 genes interfered by CNVs in patients only. In particular, we identified a de novo DKK4 duplication. DKK4 is a member of the WNT signaling pathway which is involved in the development of the anorectal region. In mice, Wnt disruption results in ARMs. Our data suggest a role for rare CNVs not only in syndromic but also in isolated ARM patients and provide a list of plausible candidate genes for the disorder.