Somatic Evolution in Non-neoplastic IBD-Affected Colon.

Inflammatory bowel disease (IBD) is a chronic inflammatory disease associated with increased risk of gastrointestinal cancers. We whole-genome sequenced 446 colonic crypts from 46 IBD patients and compared these to 412 crypts from 41 non-IBD controls from our previous publication on the mutation landscape of the normal colon. The average mutation rate of affected colonic epithelial cells is 2.4-fold that of healthy colon, and this increase is mostly driven by acceleration of mutational processes ubiquitously observed in normal colon. In contrast to the normal colon, where clonal expansions outside the confines of the crypt are rare, we observed widespread millimeter-scale clonal expansions. We discovered non-synonymous mutations in ARID1A, FBXW7, PIGR, ZC3H12A, and genes in the interleukin 17 and Toll-like receptor pathways, under positive selection in IBD. These results suggest distinct selection mechanisms in the colitis-affected colon and that somatic mutations potentially play a causal role in IBD pathogenesis.

Human SNP links differential outcomes in inflammatory and infectious disease to a FOXO3-regulated pathway.

The clinical course and eventual outcome, or prognosis, of complex diseases varies enormously between affected individuals. This variability critically determines the impact a disease has on a patient's life but is very poorly understood. Here, we exploit existing genome-wide association study data to gain insight into the role of genetics in prognosis. We identify a noncoding polymorphism in FOXO3A (rs12212067: T > G) at which the minor (G) allele, despite not being associated with disease susceptibility, is associated with a milder course of Crohn's disease and rheumatoid arthritis and with increased risk of severe malaria. Minor allele carriage is shown to limit inflammatory responses in monocytes via a FOXO3-driven pathway, which through TGFß1 reduces production of proinflammatory cytokines, including TNFα, and increases production of anti-inflammatory cytokines, including IL-10. Thus, we uncover a shared genetic contribution to prognosis in distinct diseases that operates via a FOXO3-driven pathway modulating inflammatory responses.

shaPRS: Leveraging shared genetic effects across traits or ancestries improves accuracy of polygenic scores.

We present shaPRS, a method that leverages widespread pleiotropy between traits or shared genetic effects across ancestries, to improve the accuracy of polygenic scores. The method uses genome-wide summary statistics from two diseases or ancestries to improve the genetic effect estimate and standard error at SNPs where there is homogeneity of effect between the two datasets. When there is significant evidence of heterogeneity, the genetic effect from the disease or population closest to the target population is maintained. We show via simulation and a series of real-world examples that shaPRS substantially enhances the accuracy of polygenic risk scores (PRSs) for complex diseases and greatly improves PRS performance across ancestries. shaPRS is a PRS pre-processing method that is agnostic to the actual PRS generation method, and as a result, it can be integrated into existing PRS generation pipelines and continue to be applied as more performant PRS methods are developed over time.

Somatic mutations provide important and unique insights into the biology of complex diseases.

Somatic evolution of cells within the body is well known to lead to cancers. However, spread of somatic mutations within a tissue over time may also contribute to the pathogenesis of non-neoplastic diseases. Recent years have seen the publication of many studies aiming to characterize somatic evolution in healthy tissues. A logical next step is to extend such work to diseased conditions. As our understanding of the interplay between somatic mutations and non-neoplastic disease grows, opportunities for the joint study of germline and somatic variants will present themselves. Here, we present our thoughts on the utility of somatic mutations for understanding both the causes and consequences of common complex disease and the challenges that remain for the joint study of the soma and germline.

An Integrated Taxonomy for Monogenic Inflammatory Bowel Disease.

BACKGROUND & AIMS: Monogenic forms of inflammatory bowel disease (IBD) illustrate the essential roles of individual genes in pathways and networks safeguarding immune tolerance and gut homeostasis. METHODS: To build a taxonomy model, we assessed 165 disorders. Genes were prioritized based on penetrance of IBD and disease phenotypes were integrated with multi-omics datasets. Monogenic IBD genes were classified by (1) overlapping syndromic features, (2) response to hematopoietic stem cell transplantation, (3) bulk RNA-sequencing of 32 tissues, (4) single-cell RNA-sequencing of >50 cell subsets from the intestine of healthy individuals and patients with IBD (pediatric and adult), and (5) proteomes of 43 immune subsets. The model was validated by addition of newly identified monogenic IBD defects. As a proof-of-concept, we explore the intersection between immunometabolism and antimicrobial activity for a group of disorders (G6PC3/SLC37A4). RESULTS: Our quantitative integrated taxonomy defines the cellular landscape of monogenic IBD gene expression across 102 genes with high and moderate penetrance (81 in the model set and 21 genes in the validation set). We illustrate distinct cellular networks, highlight expression profiles across understudied cell types (e.g., CD8+ T cells, neutrophils, epithelial subsets, and endothelial cells) and define genotype-phenotype associations (perianal disease and defective antimicrobial activity). We illustrate processes and pathways shared across cellular compartments and phenotypic groups and highlight cellular immunometabolism with mammalian target of rapamycin activation as one of the converging pathways. There is an overlap of genes and enriched cell-specific expression between monogenic and polygenic IBD. CONCLUSION: Our taxonomy integrates genetic, clinical and multi-omic data; providing a basis for genomic diagnostics and testable hypotheses for disease functions and treatment responses.

Fine-mapping inflammatory bowel disease loci to single-variant resolution.

Inflammatory bowel diseases are chronic gastrointestinal inflammatory disorders that affect millions of people worldwide. Genome-wide association studies have identified 200 inflammatory bowel disease-associated loci, but few have been conclusively resolved to specific functional variants. Here we report fine-mapping of 94 inflammatory bowel disease loci using high-density genotyping in 67,852 individuals. We pinpoint 18 associations to a single causal variant with greater than 95% certainty, and an additional 27 associations to a single variant with greater than 50% certainty. These 45 variants are significantly enriched for protein-coding changes (n = 13), direct disruption of transcription-factor binding sites (n = 3), and tissue-specific epigenetic marks (n = 10), with the last category showing enrichment in specific immune cells among associations stronger in Crohn's disease and in gut mucosa among associations stronger in ulcerative colitis. The results of this study suggest that high-resolution fine-mapping in large samples can convert many discoveries from genome-wide association studies into statistically convincing causal variants, providing a powerful substrate for experimental elucidation of disease mechanisms.

Common and Rare Variant Prediction and Penetrance of IBD in a Large, Multi-ethnic, Health System-based Biobank Cohort.

BACKGROUND AND AIMS: Polygenic risk scores (PRS) may soon be used to predict inflammatory bowel disease (IBD) risk in prevention efforts. We leveraged exome-sequence and single nucleotide polymorphism (SNP) array data from 29,358 individuals in the multiethnic, randomly ascertained health system-based BioMe biobank to define effects of common and rare IBD variants on disease prediction and pathophysiology. METHODS: PRS were calculated from European, African American, and Ashkenazi Jewish (AJ) reference case-control studies, and a meta-GWAS run using all three association datasets. PRS were then combined using regression to assess which combination of scores best predicted IBD status in European, AJ, Hispanic, and African American cohorts in BioMe. Additionally, rare variants were assessed in genes associated with very early-onset IBD (VEO-IBD), by estimating genetic penetrance in each BioMe population. RESULTS: Combining risk scores based on association data from distinct ancestral populations improved IBD prediction for every population in BioMe and significantly improved prediction among European ancestry UK Biobank individuals. Lower predictive power for non-Europeans was observed, reflecting in part substantially lower African IBD case-control reference sizes. We replicated associations for two VEO-IBD genes, ADAM17 and LRBA, with high dominant model penetrance in BioMe. Autosomal recessive LRBA risk alleles are associated with severe, early-onset autoimmunity; we show that heterozygous carriage of an African-predominant LRBA protein-altering allele is associated with significantly decreased LRBA and CTLA-4 expression with T-cell activation. CONCLUSIONS: Greater genetic diversity in African populations improves prediction across populations, and generalizes some VEO-IBD genes. Increasing African American IBD case-collections should be prioritized to reduce health disparities and enhance pathophysiological insight.

Rapid at-line early cell death quantification using capacitance spectroscopy.

Cell death is one of the failure modes of mammalian cell culture. Apoptosis is a regulated cell death process mainly observed in cell culture. Timely detection of apoptosis onset allows opportunities for preventive controls that ensure high productivity and consistent product quality. Capacitance spectroscopy captures the apoptosis-related cellular properties changes and thus quantifies the percentage of dying cells. This study demonstrated a quantification model that measures the percentage of apoptotic cells using a capacitance spectrometer in an at-line setup. When predicting the independent test set collected from bench-scale bioreactors, the root-mean-squared error of prediction was 8.8% (equivalent to 9.9% of the prediction range). The predicted culture evolution trajectory aligned with measured values from the flow cytometer. Furthermore, this method alarms cell death onset earlier than the traditional viability test, that is, the trypan blue exclusion test. Compared to flow cytometry (the traditional early cell death detection method), this method is rapid, simple, and less labor-intensive. In addition, this at-line setup can be easily transferred between scales (e.g., lab-scale for development to manufacturing scale), which benefits process transfers between facilities, scale-up, and other process transitions.

Correction: Insights into the genetic epidemiology of Crohn's and rare diseases in the Ashkenazi Jewish population.

[This corrects the article DOI: 10.1371/journal.pgen.1007329.].

HLA-DQA1*05 Carriage Associated With Development of Anti-Drug Antibodies to Infliximab and Adalimumab in Patients With Crohn's Disease.

BACKGROUND & AIMS: Anti-tumor necrosis factor (anti-TNF) therapies are the most widely used biologic drugs for treating immune-mediated diseases, but repeated administration can induce the formation of anti-drug antibodies. The ability to identify patients at increased risk for development of anti-drug antibodies would facilitate selection of therapy and use of preventative strategies. METHODS: We performed a genome-wide association study to identify variants associated with time to development of anti-drug antibodies in a discovery cohort of 1240 biologic-naïve patients with Crohn's disease starting infliximab or adalimumab therapy. Immunogenicity was defined as an anti-drug antibody titer ≥10 AU/mL using a drug-tolerant enzyme-linked immunosorbent assay. Significant association signals were confirmed in a replication cohort of 178 patients with inflammatory bowel disease. RESULTS: The HLA-DQA1*05 allele, carried by approximately 40% of Europeans, significantly increased the rate of immunogenicity (hazard ratio [HR], 1.90; 95% confidence interval [CI], 1.60-2.25; P = 5.88 × 10-13). The highest rates of immunogenicity, 92% at 1 year, were observed in patients treated with infliximab monotherapy who carried HLA-DQA1*05; conversely the lowest rates of immunogenicity, 10% at 1 year, were observed in patients treated with adalimumab combination therapy who did not carry HLA-DQA1*05. We confirmed this finding in the replication cohort (HR, 2.00; 95% CI, 1.35-2.98; P = 6.60 × 10-4). This association was consistent for patients treated with adalimumab (HR, 1.89; 95% CI, 1.32-2.70) or infliximab (HR, 1.92; 95% CI, 1.57-2.33), and for patients treated with anti-TNF therapy alone (HR, 1.75; 95% CI, 1.37-2.22) or in combination with an immunomodulator (HR, 2.01; 95% CI, 1.57-2.58). CONCLUSIONS: In an observational study, we found a genome-wide significant association between HLA-DQA1*05 and the development of antibodies against anti-TNF agents. A randomized controlled biomarker trial is required to determine whether pretreatment testing for HLA-DQA1*05 improves patient outcomes by helping physicians select anti-TNF and combination therapies. ClinicalTrials.gov ID: NCT03088449.

Amino acid residues in five separate HLA genes can explain most of the known associations between the MHC and primary biliary cholangitis.

Primary Biliary Cholangitis (PBC) is a chronic autoimmune liver disease characterised by progressive destruction of intrahepatic bile ducts. The strongest genetic association is with HLA-DQA1*04:01, but at least three additional independent HLA haplotypes contribute to susceptibility. We used dense single nucleotide polymorphism (SNP) data in 2861 PBC cases and 8514 controls to impute classical HLA alleles and amino acid polymorphisms using state-of-the-art methodologies. We then demonstrated through stepwise regression that association in the HLA region can be largely explained by variation at five separate amino acid positions. Three-dimensional modelling of protein structures and calculation of electrostatic potentials for the implicated HLA alleles/amino acid substitutions demonstrated a correlation between the electrostatic potential of pocket P6 in HLA-DP molecules and the HLA-DPB1 alleles/amino acid substitutions conferring PBC susceptibility/protection, highlighting potential new avenues for future functional investigation.

Insights into the genetic epidemiology of Crohn's and rare diseases in the Ashkenazi Jewish population.

As part of a broader collaborative network of exome sequencing studies, we developed a jointly called data set of 5,685 Ashkenazi Jewish exomes. We make publicly available a resource of site and allele frequencies, which should serve as a reference for medical genetics in the Ashkenazim (hosted in part at https://ibd.broadinstitute.org, also available in gnomAD at http://gnomad.broadinstitute.org). We estimate that 34% of protein-coding alleles present in the Ashkenazi Jewish population at frequencies greater than 0.2% are significantly more frequent (mean 15-fold) than their maximum frequency observed in other reference populations. Arising via a well-described founder effect approximately 30 generations ago, this catalog of enriched alleles can contribute to differences in genetic risk and overall prevalence of diseases between populations. As validation we document 148 AJ enriched protein-altering alleles that overlap with "pathogenic" ClinVar alleles (table available at https://github.com/macarthur-lab/clinvar/blob/master/output/clinvar.tsv), including those that account for 10-100 fold differences in prevalence between AJ and non-AJ populations of some rare diseases, especially recessive conditions, including Gaucher disease (GBA, p.Asn409Ser, 8-fold enrichment); Canavan disease (ASPA, p.Glu285Ala, 12-fold enrichment); and Tay-Sachs disease (HEXA, c.1421+1G>C, 27-fold enrichment; p.Tyr427IlefsTer5, 12-fold enrichment). We next sought to use this catalog, of well-established relevance to Mendelian disease, to explore Crohn's disease, a common disease with an estimated two to four-fold excess prevalence in AJ. We specifically attempt to evaluate whether strong acting rare alleles, particularly protein-truncating or otherwise large effect-size alleles, enriched by the same founder-effect, contribute excess genetic risk to Crohn's disease in AJ, and find that ten rare genetic risk factors in NOD2 and LRRK2 are enriched in AJ (p < 0.005), including several novel contributing alleles, show evidence of association to CD. Independently, we find that genomewide common variant risk defined by GWAS shows a strong difference between AJ and non-AJ European control population samples (0.97 s.d. higher, p<10-16). Taken together, the results suggest coordinated selection in AJ population for higher CD risk alleles in general. The results and approach illustrate the value of exome sequencing data in case-control studies along with reference data sets like ExAC (sites VCF available via FTP at ftp.broadinstitute.org/pub/ExAC_release/release0.3/) to pinpoint genetic variation that contributes to variable disease predisposition across populations.

Terahertz Time of Flight Spectroscopy as a Coating Thickness Reference Method for Partial Least Squares Near Infrared Spectroscopy Models.

Near infrared spectroscopy (NIRS) is often used during the tablet coating process to assess coating thickness. As the coating process proceeds, the increase and decrease in NIRS signal from both the coating formulation and tablet core has been related to coating thickness. Partial least-squares models are often generated relating NIRS spectra to reference coating thickness measurements for in-line and/or at-line monitoring of the coating process. This study investigated the effect of the reference coating thickness measurements on the accuracy of the model. The two primary reference techniques used were weight gain-based coating thickness and terahertz-based coating thickness. Most NIRS coating thickness models currently use weight gain-based reference values; however, terahertz-time-of-flight spectroscopy (THz-TOF) offers a more direct reference coating thickness measurement. Results showed that the accuracy of the NIRS coating thickness model significantly improved when terahertz-based coating thickness measurements were used as reference when compared to weight gain-based coating thickness measurements. Therefore, the application of THz-TOF as a reference method is further demonstrated.

Factors Associated With Outcomes of Patients With Primary Sclerosing Cholangitis and Development and Validation of a Risk Scoring System.

We sought to identify factors that are predictive of liver transplantation or death in patients with primary sclerosing cholangitis (PSC), and to develop and validate a contemporaneous risk score for use in a real-world clinical setting. Analyzing data from 1,001 patients recruited to the UK-PSC research cohort, we evaluated clinical variables for their association with 2-year and 10-year outcome through Cox-proportional hazards and C-statistic analyses. We generated risk scores for short-term and long-term outcome prediction, validating their use in two independent cohorts totaling 451 patients. Thirty-six percent of the derivation cohort were transplanted or died over a cumulative follow-up of 7,904 years. Serum alkaline phosphatase of at least 2.4 × upper limit of normal at 1 year after diagnosis was predictive of 10-year outcome (hazard ratio [HR] = 3.05; C = 0.63; median transplant-free survival 63 versus 108 months; P < 0.0001), as was the presence of extrahepatic biliary disease (HR = 1.45; P = 0.01). We developed two risk scoring systems based on age, values of bilirubin, alkaline phosphatase, albumin, platelets, presence of extrahepatic biliary disease, and variceal hemorrhage, which predicted 2-year and 10-year outcomes with good discrimination (C statistic = 0.81 and 0.80, respectively). Both UK-PSC risk scores were well-validated in our external cohort and outperformed the Mayo Clinic and aspartate aminotransferase-to-platelet ratio index (APRI) scores (C statistic = 0.75 and 0.63, respectively). Although heterozygosity for the previously validated human leukocyte antigen (HLA)-DR*03:01 risk allele predicted increased risk of adverse outcome (HR = 1.33; P = 0.001), its addition did not improve the predictive accuracy of the UK-PSC risk scores. Conclusion: Our analyses, based on a detailed clinical evaluation of a large representative cohort of participants with PSC, furthers our understanding of clinical risk markers and reports the development and validation of a real-world scoring system to identify those patients most likely to die or require liver transplantation.

Association of Genetic Variants in NUDT15 With Thiopurine-Induced Myelosuppression in Patients With Inflammatory Bowel Disease.

Importance: Use of thiopurines may be limited by myelosuppression. TPMT pharmacogenetic testing identifies only 25% of at-risk patients of European ancestry. Among patients of East Asian ancestry, NUDT15 variants are associated with thiopurine-induced myelosuppression (TIM). Objective: To identify genetic variants associated with TIM among patients of European ancestry with inflammatory bowel disease (IBD). Design, Setting, and Participants: Case-control study of 491 patients affected by TIM and 679 thiopurine-tolerant unaffected patients who were recruited from 89 international sites between March 2012 and November 2015. Genome-wide association studies (GWAS) and exome-wide association studies (EWAS) were conducted in patients of European ancestry. The replication cohort comprised 73 patients affected by TIM and 840 thiopurine-tolerant unaffected patients. Exposures: Genetic variants associated with TIM. Main Outcomes and Measures: Thiopurine-induced myelosuppression, defined as a decline in absolute white blood cell count to 2.5 × 109/L or less or a decline in absolute neutrophil cell count to 1.0 × 109/L or less leading to a dose reduction or drug withdrawal. Results: Among 1077 patients (398 affected and 679 unaffected; median age at IBD diagnosis, 31.0 years [interquartile range, 21.2 to 44.1 years]; 540 [50%] women; 602 [56%] diagnosed as having Crohn disease), 919 (311 affected and 608 unaffected) were included in the GWAS analysis and 961 (328 affected and 633 unaffected) in the EWAS analysis. The GWAS analysis confirmed association of TPMT (chromosome 6, rs11969064) with TIM (30.5% [95/311] affected vs 16.4% [100/608] unaffected patients; odds ratio [OR], 2.3 [95% CI, 1.7 to 3.1], P = 5.2 × 10-9). The EWAS analysis demonstrated an association with an in-frame deletion in NUDT15 (chromosome 13, rs746071566) and TIM (5.8% [19/328] affected vs 0.2% [1/633] unaffected patients; OR, 38.2 [95% CI, 5.1 to 286.1], P = 1.3 × 10-8), which was replicated in a different cohort (2.7% [2/73] affected vs 0.2% [2/840] unaffected patients; OR, 11.8 [95% CI, 1.6 to 85.0], P = .03). Carriage of any of 3 coding NUDT15 variants was associated with an increased risk (OR, 27.3 [95% CI, 9.3 to 116.7], P = 1.1 × 10-7) of TIM, independent of TPMT genotype and thiopurine dose. Conclusions and Relevance: Among patients of European ancestry with IBD, variants in NUDT15 were associated with increased risk of TIM. These findings suggest that NUDT15 genotyping may be considered prior to initiation of thiopurine therapy; however, further study including additional validation in independent cohorts is required.

Feedforward and Feedback Control of a Pharmaceutical Coating Process.

This work demonstrates the use of a combination of feedforward and feedback loops to control the controlled release coating of theophylline granules. Feedforward models are based on the size distribution of incoming granules and are used to set values for the airflow in the fluid bed processor and the target coat weight to be applied to the granules. The target coat weight of the granules is controlled by a feedback loop using NIR spectroscopy to monitor the progress of the process. By combining feedforward and feedback loops, significant variation in the size distributions and ambient conditions were accommodated in the fluid bed coating of the granules and a desired dissolution profile was achieved. The feedforward component of the control system was specifically tested by comparing the performance of the control system with and without this element by Monte Carlo simulation.

Method Development and Validation of an Inline Process Analytical Technology Method for Blend Monitoring in the Tablet Feed Frame Using Raman Spectroscopy.

Inline process analytical technology sensors are the key elements to enable continuous manufacturing. They facilitate real-time monitoring of critical quality attributes of both intermediate materials and finished products. The aim of this study was to demonstrate method development and validation for inline and offline calibration strategies to determine the blend content during tablet compression via Raman spectroscopy. An inline principal component regression model was developed from Raman spectra collected in the feed frame. At the same time, an offline study was conducted over a small amount of the calibration blends using an in-house moving powder setup to simulate the environment of the feed frame. The model developed offline was able to predict the active ingredient content after a bias correction and used only a fraction of the material. The offline method can serve as a simple method to facilitate calibration development when the time and access to the press is limited. The study takes into consideration, the necessary components of method development and offers perspectives on the validation of an inline process analytics method. Method testing and validation was performed for the inline process analytical technology method. The established Raman method was demonstrated as suitable for the determination of bulk assay of the active ingredient in powders inside the feed frame for use during batch and continuous manufacturing processes.

Host-microbe interactions have shaped the genetic architecture of inflammatory bowel disease.

Crohn's disease and ulcerative colitis, the two common forms of inflammatory bowel disease (IBD), affect over 2.5 million people of European ancestry, with rising prevalence in other populations. Genome-wide association studies and subsequent meta-analyses of these two diseases as separate phenotypes have implicated previously unsuspected mechanisms, such as autophagy, in their pathogenesis and showed that some IBD loci are shared with other inflammatory diseases. Here we expand on the knowledge of relevant pathways by undertaking a meta-analysis of Crohn's disease and ulcerative colitis genome-wide association scans, followed by extensive validation of significant findings, with a combined total of more than 75,000 cases and controls. We identify 71 new associations, for a total of 163 IBD loci, that meet genome-wide significance thresholds. Most loci contribute to both phenotypes, and both directional (consistently favouring one allele over the course of human history) and balancing (favouring the retention of both alleles within populations) selection effects are evident. Many IBD loci are also implicated in other immune-mediated disorders, most notably with ankylosing spondylitis and psoriasis. We also observe considerable overlap between susceptibility loci for IBD and mycobacterial infection. Gene co-expression network analysis emphasizes this relationship, with pathways shared between host responses to mycobacteria and those predisposing to IBD.

Heterozygous loss-of-function mutations in YAP1 cause both isolated and syndromic optic fissure closure defects.

Exome sequence analysis of affected individuals from two families with autosomal-dominant inheritance of coloboma identified two different cosegregating heterozygous nonsense mutations (c.370C>T [p.Arg124*] and c. 1066G>T [p.Glu356*]) in YAP1. The phenotypes of the affected families differed in that one included no extraocular features and the other manifested with highly variable multisystem involvement, including hearing loss, intellectual disability, hematuria, and orofacial clefting. A combined LOD score of 4.2 was obtained for the association between YAP1 loss-of-function mutations and the phenotype in these families. YAP1 encodes an effector of the HIPPO-pathway-induced growth response, and whole-mount in situ hybridization in mouse embryos has shown that Yap1 is strongly expressed in the eye, brain, and fusing facial processes. RT-PCR showed that an alternative transcription start site (TSS) in intron 1 of YAP1 and Yap1 is widely used in human and mouse development, respectively. Transcripts from the alternative TSS are predicted to initiate at codon Met179 relative to the canonical transcript (RefSeq NM_001130145). In these alternative transcripts, the c.370C>T mutation in family 1305 is within the 5' UTR and cannot result in nonsense-mediated decay (NMD). The c. 1066G>T mutation in family 132 should result in NMD in transcripts from either TSS. Amelioration of the phenotype by the alternative transcripts provides a plausible explanation for the phenotypic differences between the families.

Monoallelic and biallelic mutations in MAB21L2 cause a spectrum of major eye malformations.

We identified four different missense mutations in the single-exon gene MAB21L2 in eight individuals with bilateral eye malformations from five unrelated families via three independent exome sequencing projects. Three mutational events altered the same amino acid (Arg51), and two were identical de novo mutations (c.151C>T [p.Arg51Cys]) in unrelated children with bilateral anophthalmia, intellectual disability, and rhizomelic skeletal dysplasia. c.152G>A (p.Arg51His) segregated with autosomal-dominant bilateral colobomatous microphthalmia in a large multiplex family. The fourth heterozygous mutation (c.145G>A [p.Glu49Lys]) affected an amino acid within two residues of Arg51 in an adult male with bilateral colobomata. In a fifth family, a homozygous mutation (c.740G>A [p.Arg247Gln]) altering a different region of the protein was identified in two male siblings with bilateral retinal colobomata. In mouse embryos, Mab21l2 showed strong expression in the developing eye, pharyngeal arches, and limb bud. As predicted by structural homology, wild-type MAB21L2 bound single-stranded RNA, whereas this activity was lost in all altered forms of the protein. MAB21L2 had no detectable nucleotidyltransferase activity in vitro, and its function remains unknown. Induced expression of wild-type MAB21L2 in human embryonic kidney 293 cells increased phospho-ERK (pERK1/2) signaling. Compared to the wild-type and p.Arg247Gln proteins, the proteins with the Glu49 and Arg51 variants had increased stability. Abnormal persistence of pERK1/2 signaling in MAB21L2-expressing cells during development is a plausible pathogenic mechanism for the heterozygous mutations. The phenotype associated with the homozygous mutation might be a consequence of complete loss of MAB21L2 RNA binding, although the cellular function of this interaction remains unknown.