© 2014. Published by The Company of Biologists Ltd | Disease Models & Mechanisms (2014) 7, 1143-1152 doi:10.1242/dmm.016337

RESEARCH ARTICLE

Phenotyping structural abnormalities in mouse embryos using high-resolution episcopic microscopy Wolfgang J. Weninger1,*, Stefan H. Geyer1, Alexandrine Martineau2, Antonella Galli3, David J. Adams3, Robert Wilson2 and Timothy J. Mohun2

ABSTRACT genetic manipulation with the ability to study uterine embryo The arrival of simple and reliable methods for 3D imaging of mouse development that more closely resembles that of humans. embryos has opened the possibility of analysing normal and Until recently, the majority of advances in understanding genetic abnormal development in a far more systematic and comprehensive controls of murine cellular differentiation and organogenesis have manner than has hitherto been possible. This will not only help to come from studies targeting individual or individual genetic extend our understanding of normal tissue and organ development pathways for genetic perturbation. However, it has now become but, by applying the same approach to embryos from genetically feasible to adopt a more systematic approach to such studies. Since modified mouse lines, such imaging studies could also transform our 2007, the International Consortium (IKMC, knowledge of function in embryogenesis and the aetiology of www.mousephenotype.org) has been coordinating an international developmental disorders. The International Mouse Phenotyping effort to produce knockout lines in which every gene in the mouse Consortium is coordinating efforts to phenotype single gene genome is individually targeted. This work is now well advanced knockouts covering the entire mouse genome, including and attention is now focused on the goal of phenotyping adult characterising developmental defects for those knockout lines that mice from the knockout lines, an ambitious task overseen by prove to be embryonic lethal. Here, we present a pilot study of 34 the International Mouse Phenotyping Consortium (IMPC, such lines, utilising high-resolution episcopic microscopy (HREM) for www.mousephenotype.org) which is expected to be complete by comprehensive 2D and 3D imaging of homozygous null embryos and 2021 (Bradley et al., 2012; Brown and Moore, 2012). their wild-type littermates. We present a simple phenotyping protocol The IKMC knockout programme is also of enormous potential that has been developed to take advantage of the high-resolution value for understanding the origins of developmental abnormalities, images obtained by HREM and that can be used to score tissue and as it provides the opportunity to identify genes whose activity is organ abnormalities in a reliable manner. Using this approach with essential to support normal mouse embryo development in a embryos at embryonic day 14.5, we show the wide range of structural comprehensive and unbiased manner. Taking advantage of the abnormalities that are likely to be detected in such studies and the evolutionary conservation in gene regulatory functions and utilizing variability in phenotypes between sibling homozygous null embryos. the data increasingly available from clinical genetics and genome sequencing studies, it is realistic to expect major advances in our KEY WORDS: Phenotype screen, HREM, Imaging, 3D, Episcopic understanding of developmental disease mechanisms from such work (Brown et al., 2006; Cox and Brown, 2003). Furthermore, the INTRODUCTION identification of genetically modified mouse lines that model aspects The developmental abnormalities in cellular, tissue or organ of human developmental abnormalities offers the possibility of structure and function that underlie congenital defects have complex investigating the complex interactions of environmental factors that genetic, epigenetic and environmental causes. Furthermore, their can affect disease prevalence and severity. incidence and severity is likely to be the result of interactions However, a precondition for such advances is the ability to between multiple contributory factors. Unravelling the aetiology of accurately phenotype the impact of individual gene knockouts on such defects therefore requires approaches that can both identify the embryonic development. The scale of such phenotyping is genetic programmes that drive normal embryo development and considerable; all studies to date indicate that at least one third of the permit experimental investigation of the impact of non-genetic 20,000 gene knockouts will result in embryonic or perinatal lethality factors. Although evolutionary conservation of genetic programmes and, whilst a proportion will be secondary to placental defects has ensured that much can be learned from a wide range of (Adams et al., 2013; Mohun et al., 2013), the overwhelming convenient animal models, such as fruit flies, fish and frogs, a majority of recessive lethals will result from abnormalities in the widely used experimental model for developmental disease is now development of embryo tissues. Comprehensive screening therefore the mouse because it combines the ability to undertake sophisticated requires embryos from ~7000 recessive lethal lines to be studied, with multiple embryos studied to take account of variable 1Centre for Anatomy and Cell Biology & MIC, Medical University of Vienna, 1090 phenotypic penetrance. Wien, Austria. 2MRC National Institute for Medical Research, London NW7 1AA, Two observations from previous work help simplify this 3 UK. Wellcome Trust Sanger Institute, Cambridge CB10 1SA, UK. potentially daunting task. Firstly, whatever the nature of the gene *Author for correspondence ([email protected]) targeted, it is almost invariably the case that embryonic lethal knockout lines show tissue or organ abnormalities that can be This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted identified by morphological criteria, without recourse to molecular use, distribution and reproduction in any medium provided that the original work is properly or functional analysis. When combined with existing and well- attributed. documented data (e.g. Horsch et al., 2008;

Received 17 March 2014; Accepted 19 August 2014 Rodriguez-Zas et al., 2008; Seltmann et al., 2005; Zhu et al., 2007), Disease Models & Mechanisms

1143 RESEARCH ARTICLE Disease Models & Mechanisms (2014) doi:10.1242/dmm.016337

enormous amount of data potentially available for each embryo, it TRANSLATIONAL IMPACT will also require the adoption of analysis methods that are systematic Clinical issue and comprehensive, allowing the comparison of individual datasets This study is embedded in the Deciphering the Mechanisms of and reliable data mining of aggregated results. Crucially, such Developmental Disorders (DMDD) project. DMDD is part of a global methods must also be realistic for the scale of the undertaking that effort to identify mouse strains that have the potential to serve as models is required and tailored to the nature of the data that the imaging for challenging human congenital diseases. By using such models, the modality can yield. molecular and genetic pathways underlying the genesis and progression of developmental diseases can be elucidated. However, a precondition Here, we have investigated the feasibility of using high-resolution for these studies is the possibility to accurately screen the phenotypes episcopic microscopy (HREM), a simple block-face imaging of transgenic mouse lines in order to identify those defects that are procedure with resin-embedded embryos (Mohun and Weninger, associated with embryonic and pre- or perinatal lethality. In order to 2012b), to analyse the morphological phenotype of whole-mouse achieve this, high-resolution imaging of embryos should be combined embryos in a consistent and reproducible manner. We have with more advanced data analysis methods. previously shown the efficacy of this imaging method for both Results qualitative and quantitative analysis of organogenesis in mouse In this study, the high-resolution episcopic microscopy (HREM) method embryos, providing as it does, the highest resolution 3D analysis of was used for screening the phenotypes of embryos of 34 mouse strains embryo morphology currently available (Bamforth et al., 2013; that produce prenatally lethal homozygous offspring. The authors Dunlevy et al., 2010; Geyer and Weninger, 2012; Geyer et al., 2012; designed a protocol that allows ergonomic and comprehensive screening and scoring of all relevant phenotype features. They explored the range Geyer and Weninger, 2013; Mohun and Weninger, 2011; Pieles et of abnormalities that can be diagnosed with HREM and found that small al., 2007; Weninger et al., 2009). Using embryos from recessive to large defects, ranging from the histo-architectural to the macroscopic embryonic and perinatal lethal lines produced as part of the UK level, could be precisely identified with the aid of quick virtual contribution to the IKMC knockout programme, we have now resectioning techniques. A total of 58 phenotype abnormalities, which assessed the nature, range and probable severity of the structural were never described in mouse embryos, could be identified due to the phenotypes, identified using HREM, and outline a simple screening high resolution and quality of HREM data and the systematic phenotyping approach. protocol that enables HREM data to be reviewed in an efficient and timely manner. This protocol provides a way to harness the high Implications and future directions resolution of HREM imaging to support an effective and structured These results demonstrate that HREM combined with a newly designed method for undertaking systematic and comparative analysis of phenotyping protocol permits the identification of a wide range of previously unidentified structural malformations and tissue abnormalities embryo morphology. in prenatally lethal embryos of knockout mice. This indicates that a large number of defects that could be diagnosed with this approach have so RESULTS far escaped diagnosis using alternative imaging methods. Therefore, A standardised HREM screening protocol HREM should be considered to become the method that is routinely HREM image data resolves tissue and organ morphology at such a used for screening the phenotype of early mouse foetuses. This is likely high level of detail that establishing a focused, systematic and time- to advance our understanding of tissue and organ development and hopefully guide the design of new diagnostic and therapeutic regimes for efficient protocol for identifying abnormalities in structure is developmental defects. imperative if HREM is to be used in any screening procedure of comparative analysis. This precluded adopting the use of 3D modelling in our phenotyping procedure as a means of examining such malformations provide a clue to the developmental function of all organs and tissues because this would necessitate time- genes and an entry point for further research. In addition, about half consuming and technically challenging segmentation of datasets. of all affected embryos reach a stage of development (after 14 days The only use of 3D modelling was therefore in the first step of the of gestation) when normal organogenesis is largely complete and protocol, in which a volume-rendered model of the entire embryo much of the adult-like body plan is established. This therefore was used to define more precisely the embryo developmental stage, provides a useful single developmental window to screen for using the criteria of external features proposed by Theiler (Theiler, structural abnormalities, as recognised by the IMPC 1989). This was readily achieved using data subsampled to 25% recommendations for analysis of embryonic lethal gene knockouts with no appreciable loss in resolution. It also allowed screening for (Adams et al., 2013). surface abnormalities and major external malformations. Current in vivo imaging modalities lack the resolution to assess All subsequent analysis was performed using the cubic data set, embryo tissue and organ architecture that is necessary for visualising the captured 2D image set and its two digital, orthogonal phenotyping the range of defects associated with embryonic and plane counterparts (Box 1). Systematic analysis was best achieved perinatal lethality. Instead, analysis depends upon the use of ex vivo by analysing data in section sequence, using axial, sagittal and then examination of individual embryos. A variety of possible imaging coronal planes (supplementary material Figs S1-S3), rather than methods are available for such studies [such as optical projection organ by organ. Using this approach, a standard checklist of features tomography (OPT), micromagnetic resonance imaging (μMRI), readily discernible in HREM images was established (Box 1; microcomputed X-ray tomography (μCT), optical coherence supplementary material Fig. S4). This then formed the basis for tomography (OCT) and high-resolution episcopic microscopy systematic phenotyping of each embryo (including those showing (HREM)], differing in image resolution and technical complexity, multiple and complex abnormalities). Adopting this approach, but all sharing the key attribute of providing comprehensive 3D analysis and documentation of malformations of each embryo took datasets that encompass the entire embryo (Anderson et al., 2013; between 1 and 2.5 hours. Metscher, 2009; Schneider and Bhattacharya, 2004; Sharpe, 2003; Tschernig et al., 2013; Weninger et al., 2006; Wong et al., 2012; Range and penetrance of the abnormalities detected Zamyadi et al., 2010). However, success in phenotyping will depend Using the approach outlined above, 34 knockout lines were analysed

not only on the efficacy of the imaging modality. Given the (Table 1), with three homozygous null embryos compared with wild- Disease Models & Mechanisms

1144 RESEARCH ARTICLE Disease Models & Mechanisms (2014) doi:10.1242/dmm.016337

Box 1. Systematic phenotyping protocol for HREM data Step 1: Examine 3D model Step 2. Continued Size Sternum and ribs Surface abnormalities, jaws and closure defects Liver Pinna Great veins of lower body Limbs Superior mesenteric vein→vitelline vein→portal vein; umbilical vein→portal Theiler stage vein→ductus venosus (and valve)→hepatic veins→inferior vena cava→ iliac veins→hemiazygos and azygos veins Step 2: Examine axial sections Thoracal part of azygos vein and great veins of upper body Edema? Lymph sacs Brown adipose tissue Cardiac ventricles, interventricular septum and atrioventricular junction Subcutanous haemorrhage? Outflow tract and intrathoracic arteries Brain Atria and interatrial septum Ventricles→choroid plexus→tissue architecture→internal capsule Lung veins Spinal chord and spinal ganglia Pituitary gland Step 3: Examine sagittal resections Trigeminal ganglion and nerve Anus Eye, eye muscles and eye stalk and/or optic nerve Urethra Ear Proximal mesonephric duct and proximal paramesonephric duct External acoustic meatus→pharyngo tympanic tube→cochlea→semicircular Diaphragm and pericardium ducts→endolymphatic duct→auditory ossicles Possible reassessment: Nasal cavities, nasal septum and vomeronasal organ Proximal ureter, spleen, pancreas, inferior vena cava, superior Palatine mesenteric vein Tongue and hypoglossal nerve Tibia and fibula Oral cavity Ulna and radius Salivary glands Vertebral column and vertebrae (bodies→arcs) Mandible, styloid process and anlagen of incisor teeth Brain flexures Hyoid bone Innervation of vibrissae Larynx Fila olfactoria Trachea and lungs Thymus Olfactory bulb Thyroid gland Cranial nerves IX, X, XI Occipital bone and hypoglossal canal Anterior cerebral artery Topology of head arteries Pineal gland Vertebral artery→basilar artery→labyrinthine artery→posterior cerebral artery→ Commissural fibres middle cerebral artery→internal carotid artery Possible reassessment: Superior cervical ganglia Brain architecture, pituitary gland, eye muscles, larynx, edema Mammary glands Esophagus and intrathoracic topology of vagal nerve Step 4: Examine coronal resections Stomach Position of heart with respect to midline Intestine Semilunar valves Topology inside embryo→topology of umbilical hernia→position of caecum Cranial nerve VII Bile ducts and gall bladder Cranial nerve VIII Ventral and dorsal pancreas Number of ribs Spleen Possible reassessment: Kidney and proximal ureter Brain ventricles and tissue architecture, eyes, pituitary gland, palatine, Suprarenal gland lymph sacs, topology of organs, caudate nucleus Gonads Para-aortic bodies Urinary bladder and distal ureter After an initial assessment of developmental stage and external anatomy Distal mesonephric duct and distal paramesonephric duct based on 3D modelling of subsampled HREM data, the protocol outlines Umbilical and iliac arteries, abdominal aorta, superior mesenteric artery the order of organ analysis using the initial axial section plane, followed Osseous pelvis and femur by digitally derived sagittal and coronal resections. Examination of each Lower limb muscles organ proceeds from overall topology, through organ morphology to tissue Abdominal muscles architecture. For some organs, resection data provides a useful Scapula, humerus and clavicula opportunity to reassess any possible phenotypes scored earlier in the Upper limb muscles protocol and these steps are shown in italics. type littermates. Of the 102 homozygous null embryos examined, A striking result of this pilot study was the variation in phenotype 25 appeared morphologically normal, 27 showed only one or two penetrance, with some abnormalities evident in all three abnormalities, but approximately half (50) had multiple structural homozygotes, whereas the majority were found in only one or two defects (Fig. 1). In total, using the standardised phenotyping of the three (supplementary material Table S1). For some lines, this protocol, we were able to distinguish 166 different abnormalities, resulted from clear variations in severity with which individual 115 of which could be reasonably described using a set of 109 homozygotes were affected by the gene ablation, but for many, it mouse phenotype (MP) ontology terms (see supplementary material appears that variable penetrance of individual abnormalities Tables S1, S2). For a further 17 abnormalities, no appropriate MP produced a distinct but overlapping spectrum of lesions in each term was available, necessitating the use of the generic term ‘other affected embryo. At embryonic day 14.5 (E14.5), five lines showed defects’ (MP0005395). Thirty-four abnormalities were classified no detectable morphological abnormality, indicating either that the under the broad heading of ‘abnormal blood vessel morphology’ lesions responsible for embryonic and/or perinatal lethality are (MP0001614), these largely comprised the absence of segments undetectable by HREM imaging, or that they manifest themselves

and/or abnormal topology of arteries (n=19) or veins (n=15). at later stages of embryo development. Disease Models & Mechanisms

1145 RESEARCH ARTICLE Disease Models & Mechanisms (2014) doi:10.1242/dmm.016337

Table 1. The gene knockout lines produced at the Sanger Institute that were analysed Gene symbol Genetic background

Agpat3 C57BL/6Brd-Tyr;C57BL/6Dnk;C57BL/6N Akap9 C57BL/6Brd-Tyr;C57BL/6Dnk;C57BL/6N Aldh18a1 C57BL/6N;C57BL/6NTac Amfr C57BL/6Brd-Tyr;C57BL/6N Arid2 C57BL/6Brd-Tyr;C57BL/6Dnk;C57BL/6N Arid4a C57BL/6Brd-Tyr;C57BL/6Dnk;C57BL/6N Asxl1 C57BL/6Brd-Tyr;C57BL/6Dnk;C57BL/6N Cbx1 C57BL/6Brd-Tyr;C57BL/6Dnk;C57BL/6N Cenpj C57BL/6Brd-Tyr;C57BL/6N Cntfr C57BL/6Brd-Tyr;C57BL/6Dnk;C57BL/6N Csrp2bp C57BL/6Brd-Tyr;C57BL/6N Dip2b C57BL/6Brd-Tyr;C57BL/6N Dot1l C57BL/6Dnk;C57BL/6N Dscc1 C57BL/6Brd-Tyr;C57BL/6N Ehd1 C57BL/6N;C57BL/6NTac Fam134c C57BL/6N;C57BL/6NTac Gap43 C57BL/6Brd-Tyr;C57BL/6Dnk;C57BL/6N Grin1 C57BL/6N;C57BL/6NTac Ikbkb C57BL/6Brd-Tyr;C57BL/6Dnk;C57BL/6N Jarid2 C57BL/6Dnk;C57BL/6N Lmnb2 C57BL/6Brd-Tyr;C57BL/6N Ltbp1 C57BL/6Brd-Tyr;C57BL/6N Mbtd1 C57BL/6Brd-Tyr;C57BL/6Dnk;C57BL/6N Mks1 C57BL/6Brd-Tyr;C57BL/6N Mysm1 C57BL/6Dnk;C57BL/6N Fig. 1. Multiple defects in a TS21+ mouse embryo (Suv420h1−/−). 3D Nf1 C57BL/6Brd-Tyr;C57BL/6Dnk;C57BL/6N volume rendered models. (A) External aspect of a wild-type (+/+) embryo. Pds5b C57BL/6Brd-Tyr;C57BL/6N Note the intestine (in) forming the umbilical hernia and the caecum Prkab1 C57BL/6N;C57BL/6NTac (arrowhead). (B) Coronally sectioned 3D model of this embryo. Note the Prmt3 C57BL/6Brd-Tyr;C57BL/6N mandibles (m) in the left and right lower jaw, the symmetric appearance of Psat1 C57BL/6Brd-Tyr;C57BL/6Dnk;C57BL/6N the tongue (to), the diaphragm (pr) cranial to the liver, and the position of the Rhot1 C57BL/6Brd-Tyr;C57BL/6Dnk;C57BL/6N organs. (C) External aspect of a Suv420h1−/− embryo. Note the Supv3l1 C57BL/6Brd-Tyr;C57BL/6Dnk;C57BL/6N thoracoschisis, the position of the heart (c), and the arrangement of the Suv420h1 C57BL/6Brd-Tyr;C57BL/6Dnk;C57BL/6N intestine. (D) Coronally sectioned 3D model of this embryo. Note the Zc3hc1 C57BL/6Brd-Tyr;C57BL/6Dnk;C57BL/6N following malformations: the left lower jaw is missing. The tongue is thin and These 43 lines showed some variations in the precise genetic background as fixed only on the right mandible. Cardiac ventricles (v) are outside the body a result of changes in the IKMC mouse knockout generation procedure over cavity. Stomach (sto) is on the right side. The liver (li) is shifted cranially into time. All HREM image data are available online (http://embryoimaging.org). the thoracal cavity. Scale bars: 1 mm.

Classifying the severity of abnormalities defects or were strong indicators that littermates of the affected In order to rank the abnormalities in terms of severity, we split them embryo were likely to show severe life-threatening defects. into five categories (Sc1-Sc5; supplementary material Table S2). Examples are right- or left-sided retro-esophageal subclavian artery, The first severity category (Sc1) encompassed 30 abnormalities that absence of lobes of the thyroid gland, subcutaneous edema and dual individually are highly likely to cause pre- or perinatal death. inferior vena cava (Fig. 5). Examples include severe heart defects, absence of the jaws, absence The last category (Sc5) included 54 abnormalities for which the of the corticospinal tracts, absence of the hypoglossal nerves, impact on embryo development or survival was difficult to assess. intramural bleeding in the aortic wall and absence of the tongue Examples include abdominal situs ambiguous, polydactyly, muscles (Figs 2, 3). herniated liver tissue, shortening of the gall ducts, rib abnormalities The second category (Sc2) contained a further 18 abnormalities and duplication of ductus venosus and bifid ureter (Fig. 6). that could well result in pre- or perinatal lethality, especially if combined with other defects. Examples include the presence of Comparison with previously reported phenotypes multiple fibroma in the brain, the draining of the mesenteric vein Several of the lines we have studied ablate genes that have into the umbilical vein, simple perimembranous ventricular septal previously been targeted in other mouse lines and for which some defects, blood filled cysts in the liver and the unilateral absence of phenotyping data are available from the the internal carotid artery (Fig. 4). (MGI) database. Direct comparisons are difficult as much of these The third category (Sc3) comprised 32 abnormalities that published data are derived from studies of neonates or adults and individually are unlikely to cause pre- or perinatal death, but that are encompasses behavioural, biochemical and neurological studies that likely to cause a diminished extra-uterine life span. Examples are have no counterpart in our screen. defects of the optic cup, the absence of eye muscles, simple In the case of Psd5b, previous analysis of homozygous null muscular ventricular septal defects, defects of the perirectal tissues embryos has been undertaken at multiple stages from E12.5 through and the absence of the neurohypophysis (Fig. 4). to term. Using a combination of histological analysis, Thirty of the abnormalities detected (category Sc4) were defects immunocytochemistry and in situ hybridisation, Zhang and colleagues

that are not lethal per se, but that were often combined with heart detected growth retardation, disproportional heads and limbs, facial Disease Models & Mechanisms

1146 RESEARCH ARTICLE Disease Models & Mechanisms (2014) doi:10.1242/dmm.016337

automated manner. The alternative is to use expert review of image data. Here, the challenge is to devise a procedure that is sufficiently structured to combine comprehensive phenotyping with standardised assessment, and which is sufficiently rapid in order to be useful in embryo screening protocols. In this study, we have described such an approach and present the results of a pilot study examining the morphological abnormalities in embryos at E14.5 from 34 recessive lethal knockout lines of the IKMC-coordinated gene-knockout screen. The pilot study not only provided a rigorous test of the phenotyping procedure but also offered an indication of the range, severity and penetrance of individual abnormalities that can be readily identified from HREM data. With modern 3D imaging techniques, such as HREM, accurate and reproducible phenotyping is not constrained by limitations in the available image data. Each embryo data set in our study comprised 3000-4000 serial images, with an isotropic resolution for the entire 3D data volume of only 3 μm. The range of tissues and structures that can be distinguished by using HREM enables Fig. 2. Abnormalities that cause pre- or perinatal death (Sc1). Axially abnormalities ranging from gross organ malformations to the loss of sectioned volume rendered 3D models. (A) Wild-type (+/+) embryo. Note the cranial nerves, white matter defects and the inappropriate branching whitish appearance of blood, e.g. in the right atrium (ra). (B) Zoom in on the or location of blood vessels to be detected. Nor, with modern heart shown in A. Note the blood-filled left ventricular outflow (ar), the cardiac computing resources, is there any technical difficulty in visualising cushion material (cm) that connects to the ventricle septum (vs), and the and reviewing such large data sets (~20-30 gigabytes), at least in the thickness of the ventricle walls (arrowheads). (C) Perimembranous form of 2D image stacks. Rather, the scale and detail of the image ventricular septal defect (vsd) in an Itbp1−/− embryo. (D) Double outlet right ventricle with transposition of ascending aorta (aa) and pulmonary trunk (pt) data poses the challenge of devising a simple procedure that in an Mks1−/− embryo. (E) Thin compact layer of the myocardium examines a broad range of organs and tissues in an easy and (arrowheads) of the left (lv) and right ventricles (rv) in a Mbtd1−/− embryo. br, sequential manner. bronchus; ca, caecum; da, descending aorta; e, esophagus; in, intestine; la, Having initially used 3D modelling to examine gross external left atrium; lol, lower limb; ra, right atrium; sc, spinal chord; ta, tail; tr, trachea. morphology, identify obvious malformations and assess Scale bar: 1 mm. developmental stage, our protocol only uses the captured dataset (axial or transverse) and its two computer-generated orthogonal dysmorphism, various skeleton anomalies, anomalies of the heart, counterparts, rather than the more demanding 3D modelling or abnormalities of the superior cervical ganglia, anomalies of the enteric creation of virtual, oblique 2D section views. Each orthogonal view nervous system of the large intestine and reduced numbers of germ proved to be useful for different aspects of morphology, and the cells in the gonads (Zhang et al., 2007). Likewise, in our study, we checklist of features to be scored (along with associated MP terms) have detected a similar spectrum of abnormalities, with developmental was designed to suit progressive viewing of each section series in delay and/or growth retardation, a variety heart defects, skeletal consecutive order. In choosing the features to score through the abnormalities and abnormalities of the superior cervical ganglia, protocol, we attempted to reconcile the wide-ranging nature of the which was identifiable in at least one of the three homozygotes defects that were detected by HREM imaging (ranging, for example, imaged (supplementary material Table S1). from gross organ malformations to loss of individual nerves or HREM imaging proved inappropriate for detecting abnormalities abnormalities in blood vessel branching) with the need for a in germ cells or the enteric nervous system; however, it was able to procedure that is suitably straightforward and rapid for use in the detect many additional defects, which were unreported in the earlier systematic screening of recessive lethal lines. Using the current study. These included severe abnormalities of the head arteries, the screening protocol, it is feasible for a single anatomist to analyse an suprarenal gland and derivatives of the vitelline veins, all of which entire knockout line in a day, including documenting scored constituted serious lesions that could impair embryo viability. phenotypes with formal descriptions [e.g. using MP or mouse Interestingly, these included the absence of the hypoglossal nerve in anatomical entity (MA) ontologies] (Eppig et al., 2012; Hayamizu one of the embryos. Other minor phenotypes evident in HREM et al., 2005; Smith and Eppig, 2009; Smith and Eppig, 2012) and included enlarged lymph sacs, eye lens defects, abnormalities of the simple image annotation. Perhaps counter-intuitively, unless detailed liver and an abnormal topology of the intestine. morphometry is envisaged, the bottleneck facing embryonic lethal screens is unlikely to be the phenotyping itself, but rather the rate DISCUSSION production of mutant lines and their embryos for analysis. By combining automated data acquisition with very high resolution, To describe the abnormalities in a standardised manner, we HREM provides a powerful way to examine embryo morphology in adopted the internationally acknowledged MP ontology system both 2D image stacks and 3D models. However, the extent of the (http://www.informatics.jax.org/searches/MP_form.shtml). For most structural detail it can yield raises important questions about how of the abnormalities, this ontology provided terms to describe the best to score embryo and organ structure in a manner that is suitable abnormalities with sufficient accuracy; however, significant for comparative analysis or screening studies. exceptions were the majority of blood vessel abnormalities, for One emerging possibility is the use of computer algorithms to which no appropriate MP terms were available. This was because compare embryos, an approach that can, in principle, be extended to many of the vessel abnormalities that we detected were topological HREM data (unpublished data), but which is currently rather limited rather than structural in nature, for example, vertebral arteries that

in the range of morphological structures that can be defined in an did not enter the spinal canal between the atlas and occipital bone. Disease Models & Mechanisms

1147 RESEARCH ARTICLE Disease Models & Mechanisms (2014) doi:10.1242/dmm.016337

Fig. 3. Abnormalities that cause pre- or perinatal death (Sc1). 2D orthogonal plane images. (A-C) Control sections (+/+). A corresponds to D, B to E, and C to F. Note the parasellar segment of the internal carotid artery in F. (D-F) Malformation of the adenohypophysis (ah) and absence of the neurohypophysis (nh). Note the visibility of multiple craniopharyngeal ducts (cd) in the inlays (Asxl1−/−). (D) Axial section (original section plane). Note the narrow basilar (ba) and vertebral (va) arteries. Scale bar: 1 mm. (E) Sagittal resection. Note the normally formed olfactory bulb (ob). (F) Coronal resection. Note the position of the palatine plate (pa) lateral to the tongue (to). (G) Missing tongue (to) innervation in a Dip2b−/− embryo. No hypoglossal nerve (arrowhead) exists lateral to the hypoglossus muscle (hgm). Axial section. (H) Double outlet right ventricle (rv), which is associated with transposition of the ascending aorta (aao) and pulmonary trunk (pt) and right-sided descending aorta (da) in an Mks1−/− embryo. Axial section. (I) Aortopulmonary window (apw) in a Psat1−/− embryo. Coronal resection. (J) Additional cusp of pulmonary valve (pva) in a Cbx1−/− embryo. Coronal resection. (K-N) Controls for G to J. (K) Normal position of the hypoglossal nerve (arrow in inlay). Note the vomeronasal organ (voo). (l) Normal position of the subvalvular aorta (sva). Note the left (lpb) and right (rpb) principal bronchii and the paravertebral sympathetic ganglion (spg). (M) Normal topology of ascending aorta and pulmonary trunk. (N) Normal number of cusps (*) of the pulmonary valve. ah, adenohypophysis; aov, aortic valve; apw, aortopulmonary window; ba, basilar arteries; cd, craniopharyngeal ducts; co, cochlea; da, descending aorta; hb, hyoid bone; III, third ventricle; la, larynx; lat, left atrial appendix; lpb, left principal bronchii; lv, left ventricle; ma, mandible; nh, neurohypophysis; ns, nasal septum; ob, olfactory bulb; pa, palatine plate; pt, pulmonary trunk; pva, pulmonary valve; rho, rhombencephalon; rpb, right principal bronchii; sag, salivary gland; sc, spinal chord; sgm, styloglossus muscle; sp, sphenoid bone; spg, paravertebral sympathetic ganglion; sva, subvalvular aorta; tg, trigeminal ganglion.

Similarly, for embryos in which no artery existed in the foramina of tissues. According to the gross appearance, we have classified these the transverse processus of the atlas or other cervical vertebrae, the defects as neurofibroma, but confirmation of this would require closest appropriate MP term was ‘abnormal vertebral artery immunohistological examination. morphology’ (MP 0011513). However, this does not in our opinion Our phenotyping protocol is based on qualitative judgments of adequately describe a defect resulting from agenesis of a normal morphology and cannot detect with any reliability any changes in vascular segment. Any phenotyping screen will therefore need to overall organ volume or shape that remain subtle in 2D images. Such generate additional agreed terms to meaningfully describe and quantitative measurements are of course possible from 3D data sets annotate the full range of abnormalities detected. and could add a further and truly objective aspect to embryo Furthermore, HREM data only permits the identification of phenotyping. Indeed, such an approach, initially pioneered for abnormal tissue arrangements and is not capable of providing comparison of MRI data has now been adapted to the higher information about abnormal cell morphology or abnormal resolution data sets that have been obtained using μCT. Although the composition of the non-fibrous components of the extracellular results of such studies are encouraging, they also demonstrate the matrix. Thus, it cannot reliably establish the true nature of tissues. challenges facing the development of automated phenotyping. For example, the abnormal structures inside the brain tissue shown Automated comparison of mutant 3D data sets with statistically

in Fig. 4 could in principle be caused by the proliferation of various averaged reference data requires manual identification and Disease Models & Mechanisms

1148 RESEARCH ARTICLE Disease Models & Mechanisms (2014) doi:10.1242/dmm.016337

Fig. 5. Abnormalities that are combined with or are indicators for life- threatening defects (Sc4). (A) Coronally sectioned, volume rendered, 3D model of a wild-type (+/+) embryo. Note the position of the stomach (sto), highlighted in C. Note the positions of the ventral and dorsal pancreas (pa). (B) Mks1−/− embryo with a right-sided stomach, highlighted in D. Note that there is only right-sided pancreas tissue. (E) Symmetric left- and right-sided (arrowheads) aortic arch, forming a vascular ring around the esophagus (Msk1−/−). Note the abnormal position of the atrium (at). (F) Hypoplasia of the left lobe (lty) of the thyroid gland in a Ltbp1−/− embryo. (G) Retro-esophageal left subclavian artery (al) in a Csrp2bp−/− embryo. Note the inverse position of Fig. 4. Severe and dangerous defects (Sc2, Sc3). Potential pre- or the descending aorta (da). (H) Retro-esophageal right subclavian artery in a perinatal life-threatening defects (Sc2 in A-F). Abnormalities that reduce life Mbtd1−/− embryo. e, esophagus; ki, kidney; li, liver; rty, right thyroid lobe; sa, span (Sc3 in G-L). (A) Axial section through a wild-type (+/+) embryo. Note subclavian artery; t, trachea; tg, trigeminal ganglion; th, thymus; ul, upper the intestine (in) forming the umbilical hernia and the liver (li). (B) Zoom into limb; v, ventricle. Scale bars: 1 mm. the junction of umbilical hernia and body. Note the umbilical vein (uv) right to the intestine and the mesenterial vessels (mv) between the intestinal segments leaving and entering the body. (C) A Pds5b−/− embryo in which the noted that because normal organs and tissues show differing degrees vitelline vein (mve) connects with the umbilical vein. (D) Hemangioma-like of variability in dimensions (in part, presumably, owing to a degree of lesion (bl) in the pericardium of an Mks1−/− embryo. Sagittal resection. heterochrony in development), the sensitivity and confidence with (E) Fibroma (arrowheads) in the telencephalon of an Mks1−/− embryo. Axial which size changes can be detected by automated procedures varies −/− section. (F) Muscular ventricular septal defect (vsd) in an Mks1 embryo. considerably, potentially necessitating the study of substantially larger Coronally sectioned 3D model. (G) Axial section through a wild-type embryo. mutant embryo cohorts and a consequent increase in screening costs. Note the form and position of the eyeball (eb), lens (le), optic nerve (on) and eye muscles (em). (H) Various eye defects are indicated with arrowheads. Furthermore, as previously noted (Anderson et al., 2013; Wong et al., Vacuolated lens at top (Csrp2bp−/−) left. Displaced lens (black arrowhead) 2012; Zamyadi et al., 2010), a number of organs are sufficiently and retro-lental bleeding (white arrowhead) are visible at the top right variable in morphology that they are unsuited for inclusion in (Ehd1−/−). A malformed eye ball (abnormal optic cup) is shown at the bottom automated detection. left (Zc3hc1−/−). An absent lens (vesicle) is marked at the bottom right A striking result of this study was the extent to which individual −/− −/− (Ehd1 ). (I) A missing eye ball (anophthalmia) in an Mks1 embryo. Note recessive lethal gene knockouts resulted in pleiotropic effects. Our that the eye muscles do exist. (J) Sagittal resection through a wild-type embryo. (K) Zoom in of the head region of embryo displayed in J. Note the finding that half of the mutant embryos studied had multiple pineal gland vesicle (p). (l) Missing pineal gland vesicle in an Mks1−/− morphological defects could have several explanations. One embryo. Note that also the rest of the diencephalon (di) is malformed. co, limitation inherent in screens of this nature is the adoption of a cochlea; gb, gall bladder; I, first brain ventricle; II, second brain ventricle; III, single time point (in this case E14.5) at which to compare mutant third brain ventricle; in, intestine; lu, lung; lv, left cardiac ventricle; me, and wild-type embryos, with the result that the malformations mesencephalon; rho, rhombencephalon; rv, right cardiac ventricle; sc, spinal detected could be both primary results of the gene knockout and chord; to, tongue. Scale bars: 1 mm. secondary abnormalities resulting from perturbed development. Multiple phenotypes could also indicate several distinct roles for the segmentation of individual organs, regions or tissues (Cleary et al., targeted gene during tissue differentiation or morphogenesis. 2011; Dorr et al., 2008; Kovačević et al., 2005; Wong et al., 2012), a Resolving these possibilities ultimately requires a detailed study of procedure which has only so far proved effective with a subset of the targeted gene expression and function; however, some clues embryo organs. The alternative is comparison that is independent of might be provided by the nature of the malformations themselves if organ identification, a procedure which can yield ‘heat maps’ of voxel they are consistent with a single developmental origin (e.g. neural pattern variation that can help identify regions of possible structural crest derivatives) or reminiscent of an established congenital variation (Wong et al., 2014). Such approaches could in principle be syndrome. Gene knockouts might also contribute to a lethal embryo applied to HREM data, taking advantage of the greatly increased phenotype through direct impairment of placental function or resolution of individual structures to expand the range of features that development and, for this reason, any comprehensive screen to

can be included in the automated comparisons. However, it should be analyse recessive lethals should also include systematic assessment Disease Models & Mechanisms

1149 RESEARCH ARTICLE Disease Models & Mechanisms (2014) doi:10.1242/dmm.016337

We used E14.5 embryos in our study and found a wide variation in the stages, as assessed by the more detailed criteria of Theiler, embryos ranging in fact from stages 21 to 23. Some apparent abnormalities (for example minor ventricular septal defects) could therefore result from developmental heterochrony or indeed from developmental delay, resulting from the gene deletion itself. Although distinguishing a genuine malformation from developmental delay or natural developmental variation will remain a subjective judgment, the key to its success will be to use high resolution imaging techniques to generate detailed maps of organ development over the transition from the embryonic to the fetal period, based on sufficient numbers to reveal the extent of normal variation (see for example http://embryoimaging.org) or the prevalence and penetrance of particular abnormalities that are associated with a particular inbred mouse background. Our pilot study was performed in the absence of such a comprehensive compendium of data and was further complicated by the variations in genetic background that have occurred as the IKMC knockout programme has evolved. A second challenge for maximising the effectiveness of phenotyping is to establish the most effective methods for achieving blood removal from embryos. In most imaging methods, including HREM, blood appears as dense material, obscuring adjacent structures, especially within the heart and great vessels. Because 25-30% of all knockout embryos are expected to show a cardiac phenotype at E14.5-15.5 (Adams et al., 2013), optimising exsanguination procedures will have a major impact on successful cardiac phenotyping. Fig. 6. Abnormalities that are not obvious indicators for severe defects A number of imaging modalities have previously been used to (Sc5). (A) Transparent volume rendered 3D model of a wild-type (+/+) screen mouse embryos for morphological defects, most notably embryo combined with surface rendered 3D models of the stomach (sto), μMRI (Bamforth et al., 2012; Berrios-Otero et al., 2009; Cleary et −/− pancreas (pa), cystic duct (bd) and intestine (in). (B) Agpat3 embryo of the al., 2011; Cleary et al., 2009; Schneider et al., 2004; Turnbull and same developmental stage and displayed in the same modality as A. Note the shortness of the cystic duct. (C) Zoom in of A. The cavity of the stomach Mori, 2007) and μCT (Degenhardt et al., 2010; Johnson et al., 2006; and the lumen of the intestine (in) are yellow, the wall of the intestine is dark Wong et al., 2012). In testing HREM as a screening procedure, we brown and indicates the umbilical hernia. (D) Zoom in of B. (E) Coronal attempted to identify the extent to which the enhanced imaging resection through the head of a wild-type embryo. Note the tongue (to) and resolution translated into improved phenotyping. As an the nasal septum (s) between the nasal cavities (nc). (F) Prmt3−/− embryo approximation, we have categorised the abnormalities that we with two nasal septa and a connective tissue space in between. (G) Dual detected into three groups, macroscopic (m), intermediate (i) and ductus venosus (dv) in an axial section of a Prkab1−/− embryo. (H) Axial section through a wild-type embryo. Note the position of the valve of the histologic (h), according to the level of resolution judged as being ductus venosus (va). (I) Axial section through a Psat1−/− embryo. Note the necessary for their detection (supplementary material Table S2). abnormal position of the valve of the ductus venosus at the connection with a Interestingly, approximately one third (58 of 160) of all the left-sided liver vein (liv). (J) Abnormal ureters in an Ehd1−/− embryo. Note two abnormalities identified are likely to be detected only with high- ureters on the right side and a bifid ureter on the left side (arrowheads). aa, resolution (HREM or histology) techniques. These include a large abdominal aorta; dv, ductus venosus; e, esophagus; I, first ventricle; II, number of life-threatening defects. An additional 45 abnormalities second ventricle; ki, kidney; le, lens; li, liver; lol, lower limb; lu, lung; me, might, in principle, be detectable in some form by using lower mesencephalon; pi, pinna; sc, spinal chord; sn, snout; tel, telencephalon; ul, upper limb; va, valve of ductus venosus; vci, vena cava inferior. Scale bars: resolution methods (μMRI, OPT or μCT), but certain identification 1 mm. would again require confirmation by a high-resolution technique. The adoption of high-resolution screening therefore directly of placental structure (Adams et al., 2013). As this was not included contributed to reproducible diagnosis of nearly two thirds of all the in our pilot study, we cannot assess the extent to which placental abnormalities that were identified in our screen. abnormalities contributed to the range of embryo defects observed. In principle, HREM could provide a way of imaging 3D placental MATERIALS AND METHODS structure; however, the prevalence of blood in this tissue limits the Embryos usefulness of this approach and 2D assessment using standard Recessive embryonic lethal lines were produced as part of the Mouse Genetics histology is likely to be more effective. Project at the Wellcome Trust Sanger Institute (https://www.sanger.ac.uk/ One of the challenges highlighted by our pilot study is the mouseportal/). Lines were identified as recessive embryonic or perinatal lethal comparison of mutant and wild-type data sets. Effective by the absence of any homozygote mutants amongst a minimum of 28 phenotyping needs to take into account the effects of normal progeny, genotyped 2 weeks after birth (14 days postpartum). For each of 34 knockout lines, three homozygous null embryos and a single wild-type heterochrony in embryo development such that, for example, littermate control were analysed. These came from a variety of genetic amongst wild-type embryos, one might show osseous structures that backgrounds, based upon C57BL/6 (see Table 1), and were identified by yolk are more developed than its littermate, whereas its nervous system sac genotyping. Although small, such numbers afford an 80% confidence of might be less developed. This complexity is compounded by detecting phenotypes that are only 50% penetrant and form the basis of the

variability in the precise developmental stage of harvested embryos. UK screening programme (Adams et al., 2013; Mohun et al., 2013). Embryos Disease Models & Mechanisms

1150 RESEARCH ARTICLE Disease Models & Mechanisms (2014) doi:10.1242/dmm.016337

were harvested at approximately E14.5 and immediately transferred to Bouin’s Bamforth, S. D., Schneider, J. E. and Bhattacharya, S. (2012). High-throughput fixative for a minimum of 24 hours. After extensive washing in phosphate- analysis of mouse embryos by magnetic resonance imaging. Cold Spring Harb. Protoc 2012, pdb.prot067538. buffered saline and dehydration through a methanol series (10-90% in 10% Bamforth, S. D., Chaudhry, B., Bennett, M., Wilson, R., Mohun, T. J., Van Mierop, increments followed by 95% and 100%; a minimum of 1 hour each step), all L. H., Henderson, D. J. and Anderson, R. H. (2013). Clarification of the identity of samples were embedded in JB4 methacrylate resin (Polysciences) containing the mammalian fifth pharyngeal arch artery. Clin. Anat. 26, 173-182. eosin and Acridine Orange, as previously described (Mohun and Weninger, Berrios-Otero, C. A., Wadghiri, Y. Z., Nieman, B. J., Joyner, A. L. and Turnbull, D. H. (2009). Three-dimensional micro-MRI analysis of cerebral artery development in 2012a). A 1-2 mm cushion of polymerised resin was used to ensure samples mouse embryos. Magn. Reson. Med. 62, 1431-1439. remained entirely embedded within the resin without contacting the block Bradley, A., Anastassiadis, K., Ayadi, A., Battey, J. F., Bell, C., Birling, M. C., surface. Embryos were oriented in the embedding process as consistently as Bottomley, J., Brown, S. D., Bürger, A., Bult, C. J. et al. (2012). The mammalian possible, such that transverse sections were obtained from crown to rump. gene function resource: the International Knockout Mouse Consortium. Mamm. Genome 23, 580-586. Blocks were left to polymerise overnight at room temperature and then baked Brown, S. D. and Moore, M. W. (2012). Towards an encyclopaedia of mammalian at 90°C for a minimum of 48 hours to ensure a uniform, hard texture prior to gene function: the International Mouse Phenotyping Consortium. Dis. Model. Mech. sectioning. Samples were subsequently cooled and stored at 4°C before 5, 289-292. Brown, S. D., Hancock, J. M. and Gates, H. (2006). Understanding mammalian imaging. genetic systems: the challenge of phenotyping in the mouse. PLoS Genet. 2, e118. Cleary, J. O., Price, A. N., Thomas, D. L., Scambler, P. J., Kyriakopoulou, V., Imaging McCue, K., Schneider, J. E., Ordidge, R. J. and Lythgoe, M. F. (2009). Cardiac phenotyping in ex vivo murine embryos using microMRI. NMR Biomed. 22, 857-866. Embryos were analysed by HREM, using 3-μm sections, green fluorescent Cleary, J. O., Modat, M., Norris, F. C., Price, A. N., Jayakody, S. A., Martinez- filters and a Hamamatsu Orca HR CCD camera to obtain the high- Barbera, J. P., Greene, N. D. E., Hawkes, D. J., Ordidge, R. J., Scambler, P. J. et resolution images, as previously described (Mohun and Weninger, 2012b). al. (2011). Magnetic resonance virtual histology for embryos: 3D atlases for Data sets comprised ~3000-4000 8-bit images (~35 GB total) and were automated high-throughput phenotyping. Neuroimage 54, 769-778. Cox, R. D. and Brown, S. D. (2003). Rodent models of genetic disease. Curr. Opin. normalised for any fluctuations in exposure. Pixel grey levels were adjusted Genet. Dev. 13, 278-283. using a Photoshop 6 (Adobe) macro to optimise tissue contrast, and each Degenhardt, K., Wright, A. C., Horng, D., Padmanabhan, A. and Epstein, J. A. dataset was subsampled to obtain cubic voxels. (2010). Rapid 3D phenotyping of cardiovascular development in mouse embryos by All HREM image data for the 34 mutant lines analysed is publically micro-CT with iodine staining. Circ. Cardiovasc. Imaging 3, 314-322. Dorr, A. E., Lerch, J. P., Spring, S., Kabani, N. and Henkelman, R. M. (2008). High available at full resolution (http://embryoimaging.org) along with summary resolution three-dimensional brain atlas using an average magnetic resonance movies of all three orthogonal section plane stacks. image of 40 adult C57Bl/6J mice. Neuroimage 42, 60-69. Dunlevy, L., Bennett, M., Slender, A., Lana-Elola, E., Tybulewicz, V. L., Fisher, E. M. and Mohun, T. (2010). Down’s syndrome-like cardiac developmental defects in Data analysis embryos of the transchromosomic Tc1 mouse. Cardiovasc. Res. 88, 287-295. Cubic (3×3×3 μm3) voxel datasets were analysed using Amira 5.4 software Eppig, J. T., Blake, J. A., Bult, C. J., Kadin, J. A., Richardson, J. E.; Mouse (Visage Imaging) using a 64-bit PC workstation (Microsoft Windows 7; Genome Database Group (2012). The Mouse Genome Database (MGD): comprehensive resource for genetics and genomics of the . Nucleic minimum of 72 GB RAM). For the routine morphology screen of the pilot Acids Res. 40, D881-D886. study, embryo morphology was assessed in a defined sequence using both Geyer, S. H. and Weninger, W. J. (2012). Some mice feature 5th pharyngeal arch 3D volume-rendered models and 2D images from each orthogonal section arteries and double-lumen aortic arch malformations. Cells Tissues Organs 196, 90- plane, using a standard protocol (Box 1 and below). 98. Geyer, S. H. and Weninger, W. J. (2013). Metric characterization of the aortic arch of For illustration, volume rendering was combined with arbitrary section early mouse fetuses and of a fetus featuring a double lumen aortic arch plane erosion in Amira to obtain the embryo models shown in the figures. malformation. Ann. Anat. 195, 175-182. In some cases, individual organs were also digitally segmented from image Geyer, S. H., Maurer, B., Pötz, L., Singh, J. and Weninger, W. J. (2012). High- stacks in Amira and used to generate surface-rendered pseudo-coloured 3D resolution episcopic microscopy data-based measurements of the arteries of mouse embryos: evaluation of significance and reproducibility under routine conditions. organ models, which were superimposed in an appropriate location upon Cells Tissues Organs 195, 524-534. semi-transparent volume rendering of the whole embryo. Hayamizu, T. F., Mangan, M., Corradi, J. P., Kadin, J. A. and Ringwald, M. (2005). The Adult Mouse Anatomical Dictionary: a tool for annotating and integrating data. Acknowledgements Genome Biol. 6, R29. Horsch, M., Schädler, S., Gailus-Durner, V., Fuchs, H., Meyer, H., de Angelis, M. H. We thank B. Maurer for technical assistance. and Beckers, J. (2008). Systematic gene expression profiling of mouse model series reveals coexpressed genes. Proteomics 8, 1248-1256. Competing interests Johnson, J. T., Hansen, M. S., Wu, I., Healy, L. J., Johnson, C. R., Jones, G. M., The authors declare no competing financial interests. Capecchi, M. R. and Keller, C. (2006). Virtual histology of transgenic mouse embryos for high-throughput phenotyping. PLoS Genet. 2, e61. Author contributions Kovačević, N., Henderson, J. T., Chan, E., Lifshitz, N., Bishop, J., Evans, A. C., Henkelman, R. M. and Chen, X. J. (2005). A three-dimensional MRI atlas of the T.J.M. and W.J.W. designed the pilot project. A.G. and D.J.A. generated the mouse mouse brain with estimates of the average and variability. Cereb. Cortex 15, 639- lines. T.J.M., A.M. and R.W. produced the HREM data. W.J.W. and S.H.G. 645. performed the phenoptyping. T.J.M. and W.J.W. wrote the paper. Metscher, B. D. (2009). MicroCT for developmental biology: a versatile tool for high- contrast 3D imaging at histological resolutions. Dev. Dyn. 238, 632-640. Funding Mohun, T. J. and Weninger, W. J. (2011). Imaging heart development using high- This work was supported by funding from the Wellcome Trust [WT100160] resolution episcopic microscopy. Curr. Opin. Genet. Dev. 21, 573-578. Mohun, T. J. and Weninger, W. J. (2012a). Embedding embryos for high-resolution supporting the strategic award programme ‘Deciphering the Mechanisms of episcopic microscopy (HREM). Cold Spring Harb. Protoc. 2012, pdb.prot069583. Developmental Disorders’, 2013-2018. T.M. was supported by the Medical Mohun, T. J. and Weninger, W. J. (2012b). Generation of volume data by episcopic Research Council [U117562103]. three-dimensional imaging of embryos. Cold Spring Harb. Protoc. 2012, 681-682. Mohun, T., Adams, D. J., Baldock, R., Bhattacharya, S., Copp, A. J., Hemberger, Supplementary material M., Houart, C., Hurles, M. E., Robertson, E., Smith, J. C. et al. (2013). Supplementary material available online at Deciphering the Mechanisms of Developmental Disorders (DMDD): a new programme for phenotyping embryonic lethal mice. Dis. Model. Mech. 6, 562-566. http://dmm.biologists.org/lookup/suppl/doi:10.1242/dmm.016337/-/DC1 Pieles, G., Geyer, S. H., Szumska, D., Schneider, J., Neubauer, S., Clarke, K., Dorfmeister, K., Franklyn, A., Brown, S. D., Bhattacharya, S. et al. (2007). References microMRI-HREM pipeline for high-throughput, high-resolution phenotyping of murine Adams, D., Baldock, R., Bhattacharya, S., Copp, A. J., Dickinson, M., Greene, N. embryos. J. Anat. 211, 132-137. D., Henkelman, M., Justice, M., Mohun, T., Murray, S. A. et al. (2013). Rodriguez-Zas, S. L., Ko, Y., Adams, H. A. and Southey, B. R. (2008). Advancing Bloomsbury report on mouse embryo phenotyping: recommendations from the IMPC the understanding of the embryo transcriptome co-regulation using meta-, functional, workshop on embryonic lethal screening. Dis. Model. Mech. 6, 571-579. and gene network analysis tools. Reproduction 135, 213-224. Anderson, G. A., Wong, M. D., Yang, J. and Henkelman, R. M. (2013). 3D imaging, Schneider, J. E. and Bhattacharya, S. (2004). Making the mouse embryo registration, and analysis of the early mouse embryonic vasculature. Dev. Dyn. 242, transparent: identifying developmental malformations using magnetic resonance 527-538. imaging. Birth Defects Res. C Embryo Today 72, 241-249. Disease Models & Mechanisms

1151 RESEARCH ARTICLE Disease Models & Mechanisms (2014) doi:10.1242/dmm.016337

Schneider, J. E., Böse, J., Bamforth, S. D., Gruber, A. D., Broadbent, C., Clarke, Weninger, W. J., Geyer, S. H., Mohun, T. J., Rasskin-Gutman, D., Matsui, T., K., Neubauer, S., Lengeling, A. and Bhattacharya, S. (2004). Identification of Ribeiro, I., Costa, L. F., Izpisúa-Belmonte, J. C. and Müller, G. B. (2006). High- cardiac malformations in mice lacking Ptdsr using a novel high-throughput magnetic resolution episcopic microscopy: a rapid technique for high detailed 3D analysis of resonance imaging technique. BMC Dev. Biol. 4, 16-27. gene activity in the context of tissue architecture and morphology. Anat. Embryol. Seltmann, M., Horsch, M., Drobyshev, A., Chen, Y., de Angelis, M. H. and Beckers, (Berl.) 211, 213-221. J. (2005). Assessment of a systematic expression profiling approach in ENU- Weninger, W. J., Maurer, B., Zendron, B., Dorfmeister, K. and Geyer, S. H. (2009). induced mouse mutant lines. Mamm. Genome 16, 1-10. Measurements of the diameters of the great arteries and semi-lunar valves of chick Sharpe, J. (2003). Optical projection tomography as a new tool for studying embryo and mouse embryos. J. Microsc. 234, 173-190. anatomy. J. Anat. 202, 175-181. Wong, M. D., Dorr, A. E., Walls, J. R., Lerch, J. P. and Henkelman, R. M. (2012). A Smith, C. L. and Eppig, J. T. (2009). The mammalian phenotype ontology: enabling novel 3D mouse embryo atlas based on micro-CT. Development 139, 3248-3256. Wong, M. D., Maezawa, Y., Lerch, J. P. and Henkelman, R. M. (2014). Automated robust annotation and comparative analysis. Wiley Interdiscip. Rev. Syst. Biol. Med. pipeline for anatomical phenotyping of mouse embryos using micro-CT. 1, 390-399. Development 141, 2533-2541. Smith, C. L. and Eppig, J. T. (2012). The Mammalian Phenotype Ontology as a Zamyadi, M., Baghdadi, L., Lerch, J. P., Bhattacharya, S., Schneider, J. E., unifying standard for experimental and high-throughput phenotyping data. Mamm. Henkelman, R. M. and Sled, J. G. (2010). Mouse embryonic phenotyping by Genome 23, 653-668. morphometric analysis of MR images. Physiol. Genomics 42A, 89-95. Theiler, K. (1989). The House Mouse: Atlas of Embryonic Development. New York, Zhang, B., Jain, S., Song, H., Fu, M., Heuckeroth, R. O., Erlich, J. M., Jay, P. Y. and NY: Springer-Verlag. Milbrandt, J. (2007). Mice lacking sister chromatid cohesion protein PDS5B exhibit Tschernig, T., Thrane, L., Jørgensen, T. M., Thommes, J., Pabst, R. and Yelbuz, T. developmental abnormalities reminiscent of Cornelia de Lange syndrome. M. (2013). An elegant technique for ex vivo imaging in experimental research – Development 134, 3191-3201. Optical Coherence Tomography (OCT). Ann. Anat. 195, 25-27. Zhu, H., Cabrera, R. M., Wlodarczyk, B. J., Bozinov, D., Wang, D., Schwartz, R. J. Turnbull, D. H. and Mori, S. (2007). MRI in mouse developmental biology. NMR and Finnell, R. H. (2007). Differentially expressed genes in embryonic cardiac Biomed. 20, 265-274. tissues of mice lacking Folr1 gene activity. BMC Dev. Biol. 7, 128. Disease Models & Mechanisms

1152