Genes into geometry: imaging for mouse development in 3D

Morphological simulation tests the limits on phenotype discovery in 3D image analysis

In the past few decades, advances in 3D imaging have created new opportunities for reverse genetic screens. Rapidly growing datasets of 3D images of genetic knockouts require high-throughput, automated computational approaches for identifying and characterizing new phenotypes. However, exploratory, discovery-oriented image analysis pipelines used to discover these phenotypes can be difficult to validate because, by their nature, the expected outcome is not known a priori . Introducing known morphological variation through simulation can help distinguish between real phenotypic differences and random variation; elucidate the effects of sample size; and test the sensitivity and reproducibility of morphometric analyses. Here we present a novel approach for 3D morphological simulation that uses open-source, open-access tools available in 3D Slicer, SlicerMorph, and Advanced Normalization Tools in R (ANTsR). While we focus on diffusible-iodine contrast-enhanced micro-CT (diceCT) images, this approach can be used on any volumetric image. We then use our simulated datasets to test whether tensor-based morphometry (TBM) can recover our introduced differences; to test how effect size and sample size affect detectability; and to determine the reproducibility of our results. In our approach to morphological simulation, we first generate a simulated deformation based on a reference image and then propagate this deformation to subjects using inverse transforms obtained from the registration of subjects to the reference. This produces a new dataset with a shifted population mean while retaining individual variability because each sample deforms more or less based on how different or similar it is from the reference. TBM is a widely-used technique that statistically compares local volume differences associated with local deformations. Our results showed that TBM recovered our introduced morphological differences, but that detectability was dependent on the effect size, the sample size, and the region of interest (ROI) included in the analysis. Detectability of subtle phenotypes can be improved both by increasing the sample size and by limiting analyses to specific body regions. However, it is not always feasible to increase sample sizes in screens of essential genes. Therefore, methodical use of ROIs is a promising way to increase the power of TBM to detect subtle phenotypes. Generating known morphological variation through simulation has broad applicability in developmental, evolutionary, and biomedical morphometrics and is a useful way to distinguish between a failure to detect morphological difference and a true lack of morphological difference. Morphological simulation can also be applied to AI-based supervised learning to augment datasets and overcome dataset limitations.

Optical imaging of the small intestine immune compartment across scales

The limitations of 2D microscopy constrain our ability to observe and understand tissue-wide networks that are, by nature, 3-dimensional. Optical projection tomography (OPT) enables the acquisition of large volumes (ranging from micrometres to centimetres) in various tissues. We present a multi-modal workflow for the characterization of both structural and quantitative parameters of the mouse small intestine. As proof of principle, we evidence its applicability for imaging the mouse intestinal immune compartment and surrounding mucosal structures. We quantify the volumetric size and spatial distribution of Isolated Lymphoid Follicles (ILFs) and quantify the density of villi throughout centimetre-long segments of intestine. Furthermore, we exhibit the age and microbiota dependence for ILF development, and leverage a technique that we call reverse-OPT for identifying and homing in on regions of interest. Several quantification capabilities are displayed, including villous density in the autofluorescent channel and the size and spatial distribution of the signal of interest at millimetre-scale volumes. The concatenation of 3D imaging with reverse-OPT and high-resolution 2D imaging allows accurate localisation of ROIs and adds value to interpretations made in 3D. Importantly, OPT may be used to identify sparsely-distributed regions of interest in large volumes whilst retaining compatibility with high-resolution microscopy modalities, including confocal microscopy. We believe this pipeline to be approachable for a wide-range of specialties, and to provide a new method for characterisation of the mouse intestinal immune compartment.

Spatiotemporal imaging and analysis of mouse and human liver bud morphogenesis

BACKGROUND

The process of liver organogenesis has served as a paradigm for organ formation. However, there remains a lack of understanding regarding early mouse and human liver bud morphogenesis and early liver volumetric growth. Elucidating dynamic changes in liver volumes is critical for understanding organ development, implementing toxicological studies, and for modeling hPSC-derived liver organoid growth. New visualization, analysis, and experimental techniques are desperately needed.

RESULTS

Here, we combine observational data with digital resources, new 3D imaging approaches, retrospective analysis of liver volume data, mathematical modeling, and experiments with hPSC-derived liver organoids. Mouse and human liver organogenesis, characterized by exponential growth, demonstrate distinct spatial features and growth curves over time, which we mathematically modeled using Gompertz models. Visualization of liver-epithelial and septum transversum mesenchyme (STM) interactions suggests extended interactions, which together with new spatial features may be responsible for extensive exponential growth. These STM interactions are modeled with a novel in vitro human pluripotent stem cell (hPSC)-derived hepatic organoid system that exhibits cell migration.

CONCLUSIONS

Our methods enhance our understanding of liver organogenesis, with new 3D visualization, analysis, mathematical modeling, and in vitro models with hPSCs. Our approach highlights mouse and human differences and provides potential hypothesis for further investigation in vitro and in vivo.

Midgut development in rat embryos using microcomputed tomography

The development of the mammalian gut was first described more than a century ago. Since then, it has been believed that a series of highly orchestrated developmental processes occur before the intestine achieves its final formation. The key steps include the formation of the umbilicus, the so-called "physiological herniation" of the midgut into the umbilical cord, an intestinal "rotation", and the "return of the gut" into the abdominal cavity. However, this sequence of events is predominantly based on histological sections of dissected embryos, a 2D technique with methodological limitations. For a better understanding of spatial relationships in the embryo, we utilized microcomputed tomography (µCT), a nondestructive 3D imaging method. Here, we show the detailed processes and mechanisms of intestinal development in rat embryos, including the development of the umbilicus, the formation of loops inside the umbilical coelom, and the subsequent shift of these loops into the abdominal cavity. Our 3D datasets of developing intestines will substantially advance the understanding of normal mammalian midgut embryology and offer new possibilities to reveal unknown mechanisms in the pathogenesis of congenital disorders.

In Utero MRI of Mouse Embryos

Genetically engineered mouse models are used extensively as models of human development and developmental diseases. Conventional histological approaches are static and two-dimensional, and do not provide a full understanding of the dynamic, spatiotemporal changes in developing mouse embryos. Magnetic resonance imaging (MRI) offers a noninvasive and longitudinal approach for three-dimensional in utero imaging of normal and mutant mouse embryos. In this chapter, we describe MRI approaches that have been developed for imaging the living embryonic mouse brain and vasculature. Details are provided on the animal preparation and setup, MRI equipment, acquisition and reconstruction methods that have been found to be most useful for in utero MRI, including examples of applications to fetal mouse neuroimaging.

Activated Hedgehog-GLI Signaling Causes Congenital Ureteropelvic Junction Obstruction

Intrinsic ureteropelvic junction obstruction is the most common cause of congenital hydronephrosis, yet the underlying pathogenesis is undefined. Hedgehog proteins control morphogenesis by promoting GLI-dependent transcriptional activation and inhibiting the formation of the GLI3 transcriptional repressor. Hedgehog regulates differentiation and proliferation of ureteric smooth muscle progenitor cells during murine kidney-ureter development. Histopathologic findings of smooth muscle cell hypertrophy and stroma-like cells, consistently observed in obstructing tissue at the time of surgical correction, suggest that Hedgehog signaling is abnormally regulated during the genesis of congenital intrinsic ureteropelvic junction obstruction. Here, we demonstrate that constitutively active Hedgehog signaling in murine intermediate mesoderm-derived renal progenitors results in hydronephrosis and failure to develop a patent pelvic-ureteric junction. Tissue obstructing the ureteropelvic junction was marked as early as E13.5 by an ectopic population of cells expressing Ptch2, a Hedgehog signaling target. Constitutive expression of GLI3 repressor in Ptch1-deficient mice rescued ectopic Ptch2 expression and obstructive hydronephrosis. Whole transcriptome analysis of isolated Ptch2⁺ cells revealed coexpression of genes characteristic of stromal progenitor cells. Genetic lineage tracing indicated that stromal cells blocking the ureteropelvic junction were derived from intermediate mesoderm-derived renal progenitors and were distinct from the smooth muscle or epithelial lineages. Analysis of obstructive ureteric tissue resected from children with congenital intrinsic ureteropelvic junction obstruction revealed a molecular signature similar to that observed in Ptch1-deficient mice. Together, these results demonstrate a Hedgehog-dependent mechanism underlying mammalian intrinsic ureteropelvic junction obstruction.

Partitioning k-space for cylindrical three-dimensional rapid acquisition with relaxation enhancement imaging in the mouse brain

Three-dimensional rapid acquisition with relaxation enhancement (RARE) scans require the assignment of each phase encode step in two dimensions to an echo in the echo train. Although this assignment is frequently made across the entire Cartesian grid, collection of only the central cylinder of k-space by eliminating the corners in each phase encode dimension reduces the scan time by ~22% with negligible impact on image quality. The recipe for the assignment of echoes to grid points for such an acquisition is less straightforward than for the simple full Cartesian acquisition case, and has important implications for image quality. We explored several methods of partitioning k-space-exploiting angular symmetry in one extreme or emulating a cropped Cartesian acquisition in the other-and acquired three-dimensional RARE magnetic resonance imaging (MRI) scans of the ex vivo mouse brain. We evaluated each partitioning method for sensitivity to artifacts and then further considered strategies to minimize these through averaging or interleaving of echoes and by empirical phase correction. All scans were collected 16 at a time with multiple-mouse MRI. Although all schemes considered could be used to generate images, the results indicate that the emulation of a standard Cartesian echo assignment, by partitioning preferentially along one dimension within the cylinder, is more robust to artifacts. Samples at the periphery of the bore showed larger phase deviations and higher sensitivity to artifacts, but images of good quality could still be obtained with an optimized acquisition protocol. A protocol for high-resolution (40 μm) ex vivo images using this approach is presented, and has been used routinely with a success rate of 99% in over 1000 images.

Recent Progress in Magnetic Resonance Imaging of the Embryonic and Neonatal Mouse Brain

The laboratory mouse has been widely used as a model system to investigate the genetic control mechanisms of mammalian brain development. Magnetic resonance imaging (MRI) is an important tool to characterize changes in brain anatomy in mutant mouse strains and injury progression in mouse models of fetal and neonatal brain injury. Progress in the last decade has enabled us to acquire MRI data with increasing anatomical details from the embryonic and neonatal mouse brain. High-resolution ex vivo MRI, especially with advanced diffusion MRI methods, can visualize complex microstructural organizations in the developing mouse brain. In vivo MRI of the embryonic mouse brain, which is critical for tracking anatomical changes longitudinally, has become available. Applications of these techniques may lead to further insights into the complex and dynamic processes of brain development.

Morphometrics, 3D Imaging, and Craniofacial Development

Recent studies have shown how volumetric imaging and morphometrics can add significantly to our understanding of morphogenesis, the developmental basis for variation, and the etiology of structural birth defects. On the other hand, the complex questions and diverse imaging data in developmental biology present morphometrics with more complex challenges than applications in virtually any other field. Meeting these challenges is necessary in order to understand the mechanistic basis for variation in complex morphologies. This chapter reviews the methods and theory that enable the application of modern landmark-based morphometrics to developmental biology and craniofacial development, in particular. We discuss the theoretical foundations of morphometrics as applied to development and review the basic approaches to the quantification of morphology. Focusing on geometric morphometrics, we discuss the principal statistical methods for quantifying and comparing morphological variation and covariation structure within and among groups. Finally, we discuss the future directions for morphometrics in developmental biology that will be required for approaches that enable quantitative integration across the genotype-phenotype map.

4D MEMRI atlas of neonatal FVB/N mouse brain development

The widespread use of the mouse as a model system to study brain development has created the need for noninvasive neuroimaging methods that can be applied to early postnatal mice. The goal of this study was to optimize in vivo three- (3D) and four-dimensional (4D) manganese (Mn)-enhanced MRI (MEMRI) approaches for acquiring and analyzing data from the developing mouse brain. The combination of custom, stage-dependent holders and self-gated (motion-correcting) 3D MRI sequences enabled the acquisition of high-resolution (100-μm isotropic), motion artifact-free brain images with a high level of contrast due to Mn-enhancement of numerous brain regions and nuclei. We acquired high-quality longitudinal brain images from two groups of FVB/N strain mice, six mice per group, each mouse imaged on alternate odd or even days (6 3D MEMRI images at each day) covering the developmental stages between postnatal days 1 to 11. The effects of Mn-exposure, anesthesia and MRI were assessed, showing small but significant transient effects on body weight and brain volume, which recovered with time and did not result in significant morphological differences when compared to controls. Metrics derived from deformation-based morphometry (DBM) were used for quantitative analysis of changes in volume and position of a number of brain regions. The cerebellum, a brain region undergoing significant changes in size and patterning at early postnatal stages, was analyzed in detail to demonstrate the spatiotemporal characterization made possible by this new atlas of mouse brain development. These results show that MEMRI is a powerful tool for quantitative analysis of mouse brain development, with great potential for in vivo phenotype analysis in mouse models of neurodevelopmental diseases.

Mutations in the latent TGF-beta binding protein 3 (LTBP3) gene cause brachyolmia with amelogenesis imperfecta

Inherited dental malformations constitute a clinically and genetically heterogeneous group of disorders. Here, we report on four families, three of them consanguineous, with an identical phenotype, characterized by significant short stature with brachyolmia and hypoplastic amelogenesis imperfecta (AI) with almost absent enamel. This phenotype was first described in 1996 by Verloes et al. as an autosomal recessive form of brachyolmia associated with AI. Whole-exome sequencing resulted in the identification of recessive hypomorphic mutations including deletion, nonsense and splice mutations, in the LTBP3 gene, which is involved in the TGF-beta signaling pathway. We further investigated gene expression during mouse development and tooth formation. Differentiated ameloblasts synthesizing enamel matrix proteins and odontoblasts expressed the gene. Study of an available knockout mouse model showed that the mutant mice displayed very thin to absent enamel in both incisors and molars, hereby recapitulating the AI phenotype in the human disorder.

Automated pipeline for anatomical phenotyping of mouse embryos using micro-CT

The International Mouse Phenotyping Consortium (IMPC) plans to phenotype 20,000 single-gene knockout mice to gain an insight into gene function. Approximately 30% of these knockout mouse lines will be embryonic or perinatal lethal. The IMPC has selected three-dimensional (3D) imaging to phenotype these mouse lines at relevant stages of embryonic development in an attempt to discover the cause of lethality using detailed anatomical information. Rate of throughput is paramount as IMPC production centers have been given the ambitious task of completing this phenotyping project by 2021. Sifting through the wealth of data within high-resolution 3D mouse embryo data sets by trained human experts is infeasible at this scale. Here, we present a phenotyping pipeline that identifies statistically significant anatomical differences in the knockout, in comparison with the wild type, through a computer-automated image registration algorithm. This phenotyping pipeline consists of three analyses (intensity, deformation, and atlas based) that can detect missing anatomical structures and differences in volume of whole organs as well as on the voxel level. This phenotyping pipeline was applied to micro-CT images of two perinatal lethal mouse lines: a hypomorphic mutation of the Tcf21 gene (Tcf21-hypo) and a knockout of the Satb2 gene. With the proposed pipeline we were able to identify the majority of morphological phenotypes previously published for both the Tcf21-hypo and Satb2 mutant mouse embryos in addition to novel phenotypes. This phenotyping pipeline is an unbiased, automated method that highlights only those structural abnormalities that survive statistical scrutiny and illustrates them in a straightforward fashion.

Neuroanatomical phenotypes in a mouse model of the 22q11.2 microdeletion

Recurrent deletions at the 22q11.2 locus have been established as a strong genetic risk factor for the development of schizophrenia and cognitive dysfunction. Individuals with 22q11.2 deletions have a range of well-defined volumetric abnormalities in a number of critical brain structures. A mouse model of the 22q11.2 deletion (Df(16)A(+/-)) has previously been utilized to characterize disease-associated abnormalities on synaptic, cellular, neurocircuitry, and behavioral levels. We performed a high-resolution MRI analysis of mutant mice compared with wild-type littermates. Our analysis revealed a striking similarity in the specific volumetric changes of Df(16)A(+/-) mice compared with human 22q11.2 deletion carriers, including in cortico-cerebellar, cortico-striatal and cortico-limbic circuits. In addition, higher resolution magnetic resonance imaging compared with neuroimaging in human subjects allowed the detection of previously unknown subtle local differences. The cerebellar findings in Df(16)A(+/-) mice are particularly instructive as they are localized to specific areas within both the deep cerebellar nuclei and the cerebellar cortex. Our study indicates that the Df(16)A(+/-)mouse model recapitulates most of the hallmark neuroanatomical changes observed in 22q11.2 deletion carriers. Our findings will help guide the design and interpretation of additional complementary studies and thereby advance our understanding of the abnormal brain development underlying the emergence of 22q11.2 deletion-associated psychiatric and cognitive symptoms.

A coming of age: advanced imaging technologies for characterising the developing mouse

The immense challenge of annotating the entire mouse genome has stimulated the development of cutting-edge imaging technologies in a drive for novel information. These techniques promise to improve understanding of the genes involved in embryo development, at least one third of which have been shown to be essential. Aligning advanced imaging technologies with biological needs will be fundamental to maximising the number of phenotypes discovered in the coming years. International efforts are underway to meet this challenge through an integrated and sophisticated approach to embryo phenotyping. We review rapid advances made in the imaging field over the past decade and provide a comprehensive examination of the relative merits of current and emerging techniques. The aim of this review is to provide a guide to state-of-the-art embryo imaging that will enable informed decisions as to which technology to use and fuel conversations between expert imaging laboratories, researchers, and core mouse production facilities.

Mapping registration sensitivity in MR mouse brain images

Nonlinear registration algorithms provide a way to estimate structural (brain) differences based on magnetic resonance images. Their ability to align images of different individuals and across modalities has been well-researched, but the bounds of their sensitivity with respect to the recovery of salient morphological differences between groups are unclear. Here we develop a novel approach to simulate deformations on MR brain images to evaluate the ability of two registration algorithms to extract structural differences corresponding to biologically plausible atrophy and expansion. We show that at a neuroanatomical level registration accuracy is influenced by the size and compactness of structures, but do so differently depending on how much change is simulated. The size of structures has a small influence on the recovered accuracy. There is a trend for larger structures to be recovered more accurately, which becomes only significant as the amount of simulated change is large. More compact structures can be recovered more accurately regardless of the amount of simulated change. Both tested algorithms underestimate the full extent of the simulated atrophy and expansion. Finally we show that when multiple comparisons are corrected for at a voxelwise level, a very low rate of false positives is obtained. More interesting is that true positive rates average around 40%, indicating that the simulated changes are not fully recovered. Simulation experiments were run using two fundamentally different registration algorithms and we identified the same results, suggesting that our findings are generalizable across different classes of nonlinear registration algorithms.

MRI analysis of cerebellar and vestibular developmental phenotypes in Gbx2 conditional knockout mice

PURPOSE

Our aim in this study was to apply three-dimensional MRI methods to analyze early postnatal morphological phenotypes in a Gbx2 conditional knockout (Gbx2-CKO) mouse that has variable midline deletions in the central cerebellum, reminiscent of many human cerebellar hypoplasia syndromes.

METHODS

In vivo three-dimensional manganese-enhanced MRI at 100-µm isotropic resolution was used to visualize mouse brains between postnatal days 3 and 11, when cerebellum morphology undergoes dramatic changes. Deformation-based morphometry and volumetric analysis of manganese-enhanced MRI images were used to, respectively, detect and quantify morphological phenotypes in Gbx2-CKO mice. Ex vivo micro-MRI was performed after perfusion-fixation with supplemented gadolinium for higher resolution (50-µm) analysis.

RESULTS

In vivo manganese-enhanced MRI and deformation-based morphometry correctly identified known cerebellar defects in Gbx2-CKO mice, and novel phenotypes were discovered in the deep cerebellar nuclei and the vestibulo-cerebellum, both validated using histology. Ex vivo micro-MRI revealed subtle phenotypes in both the vestibulo-cerebellum and the vestibulo-cochlear organ, providing an interesting example of complementary phenotypes in a sensory organ and its associated brain region.

CONCLUSION

These results show the potential of three-dimensional MRI for detecting and analyzing developmental defects in mouse models of neurodevelopmental diseases.

High-resolution MRI of early-stage mouse embryos

Both the availability of methods to manipulate genes and the completion of the mouse genome sequence have led to the generation of thousands of genetically modified mouse lines that provide a new platform for the study of mammalian development and developmental diseases. Phenotyping of mouse embryos has traditionally been performed on fixed embryos by the use of ex vivo histological, optical and high-resolution MRI techniques. Although potentially powerful, longitudinal imaging of individual animals is difficult or impossible with conventional optical methods because of the inaccessibility of mouse embryos inside the maternal uterus. To address this problem, we present a method of imaging the mouse embryo from stages as early as embryonic day (E)10.5, close to the onset of organogenesis in most physiological systems. This method uses a self-gated MRI protocol, combined with image registration, to obtain whole-embryo high-resolution (100 µm isotropic) three-dimensional images. Using this approach, we demonstrate high contrast in the cerebral vasculature, limbs, spine and central nervous system without the use of contrast agents. These results indicate the potential of MRI for the longitudinal imaging of developing mouse embryos in utero and for future applications in analyzing mutant mouse phenotypes.

Phenotype detection in morphological mutant mice using deformation features

Large-scale global efforts are underway to knockout each of the approximately 25,000 mouse genes and interpret their roles in shaping the mammalian embryo. Given the tremendous amount of data generated by imaging mutated prenatal mice, high-throughput image analysis systems are inevitable to characterize mammalian development and diseases. Current state-of-the-art computational systems offer only differential volumetric analysis of pre-defined anatomical structures between various gene-knockout mice strains. For subtle anatomical phenotypes, embryo phenotyping still relies on the laborious histological techniques that are clearly unsuitable in such big data environment. This paper presents a system that automatically detects known phenotypes and assists in discovering novel phenotypes in muCT images of mutant mice. Deformation features obtained from non-linear registration of mutant embryo to a normal consensus average image are extracted and analyzed to compute phenotypic and candidate phenotypic areas. The presented system is evaluated using C57BL/10 embryo images. All cases of ventricular septum defect and polydactyly, well-known to be present in this strain, are successfully detected. The system predicts potential phenotypic areas in the liver that are under active histological evaluation for possible phenotype of this mouse line.

A novel 3D mouse embryo atlas based on micro-CT

The goal of the International Mouse Phenotyping Consortium (IMPC) is to phenotype targeted knockout mouse strains throughout the whole mouse genome (23,000 genes) by 2021. A significant percentage of the generated mice will be embryonic lethal; therefore, phenotyping methods tuned to the mouse embryo are needed. Methods that are robust, quantitative, automated and high-throughput are attractive owing to the numbers of mice involved. Three-dimensional (3D) imaging is a useful method for characterizing morphological phenotypes. However, tools to automatically quantify morphological information of mouse embryos from 3D imaging have not been fully developed. We present a representative mouse embryo average 3D atlas comprising micro-CT images of 35 individual C57BL/6J mouse embryos at 15.5 days post-coitum. The 35 micro-CT images were registered into a consensus average image with our automated image registration software and 48 anatomical structures were segmented manually. We report the mean and variation in volumes for each of the 48 segmented structures. Mouse organ volumes vary by 2.6-4.2% on a linear scale when normalized to whole body volume. A power analysis of the volume data reports that a 9-14% volume difference can be detected between two classes of mice with sample sizes of eight. This resource will be crucial in establishing baseline anatomical phenotypic measurements for the assessment of mutant mouse phenotypes, as any future mutant embryo image can be registered to the atlas and subsequent organ volumes calculated automatically.

Structure-function relationship of cerebral networks in experimental neuroscience: contribution of magnetic resonance imaging

The analysis of neuronal networks, their interactions in resting condition as well as during brain activation have become of great interest for a better understanding of the signal processing of the brain during sensory stimulus or cognitive tasks. Parallel to the study of the functional networks and their dynamics, the underlying network structure is highly important as it provides the basis of the functional interaction. Moreover, under pathological conditions, some nodes in such a net may be impaired and the function of the whole network affected. Mechanisms such as functional deficit and improvement, and plastic reorganization are increasingly discussed in the context of existing structural and functional networks. While many of these aspects have been followed in human and clinical studies, the experimental range is limited for obvious reasons. Here, animal experimental studies are needed as they permit longer scan times and, moreover, comparison with invasive histology. Experimental non-invasive imaging modalities are now able to perform impressive contributions. In this review we try to highlight most recent new cutting-edge developments and applications in experimental neuroscience of functional and structural networks of the brain, relying on non-invasive imaging. We focus primarily on the potential of experimental Magnetic Resonance Imaging (MRI), but also touch upon micro positron emission tomography (μPET) and optical imaging developments where they are applicable to the topic of the present review.