1
|
Thomas OO, Maga AM. Leveraging descriptor learning and functional map-based shape matching for automated anatomical Landmarking in mouse mandibles. J Anat 2025; 246:829-845. [PMID: 39814541 PMCID: PMC11996707 DOI: 10.1111/joa.14196] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2024] [Revised: 10/24/2024] [Accepted: 11/27/2024] [Indexed: 01/18/2025] Open
Abstract
Geometric morphometrics is used in the biological sciences to quantify morphological traits. However, the need for manual landmark placement hampers scalability, which is both time-consuming, labor-intensive, and open to human error. The selected landmarks embody a specific hypothesis regarding the critical geometry relevant to the biological question. Any adjustment to this hypothesis necessitates acquiring a new set of landmarks or revising them significantly, which can be impractical for large datasets. There is a pressing need for more efficient and flexible methods for landmark placement that can adapt to different hypotheses without requiring extensive human effort. This study investigates the precision and accuracy of landmarks derived from functional correspondences obtained through the functional map framework of geometry processing. We utilize a deep functional map network to learn shape descriptors, which enable us to achieve functional map-based and point-to-point correspondences between specimens in our dataset. Our methodology involves automating the landmarking process by interrogating these maps to identify corresponding landmarks, using manually placed landmarks from the entire dataset as a reference. We apply our method to a dataset of rodent mandibles and compare its performance to MALPACA's, a standard tool for automatic landmark placement. Our model demonstrates a speed improvement compared to MALPACA while maintaining a competitive level of accuracy. Although MALPACA typically shows the lowest RMSE, our models perform comparably well, particularly with smaller training datasets, indicating strong generalizability. Visual assessments confirm the precision of our automated landmark placements, with deviations consistently falling within an acceptable range for MALPACA estimates. Our results underscore the potential of unsupervised learning models in anatomical landmark placement, presenting a practical and efficient alternative to traditional methods. Our approach saves significant time and effort and provides the flexibility to adapt to different hypotheses about critical geometrical features without the need for manual re-acquisition of landmarks. This advancement can significantly enhance the scalability and applicability of geometric morphometrics, making it more feasible for large datasets and diverse biological studies.
Collapse
Affiliation(s)
- Oshane O Thomas
- Center for Development Biology and Regenerative Medicine, Seattle Children's Research Institute, Seattle, Washington, USA
| | - A Murat Maga
- Center for Development Biology and Regenerative Medicine, Seattle Children's Research Institute, Seattle, Washington, USA
- Division of Craniofacial Medicine, Department of Pediatrics, University of Washington, Seattle, Washington, USA
| |
Collapse
|
2
|
Rolfe SM, Mao D, Maga AM. Streamlining Asymmetry Quantification in Fetal Mouse Imaging: A Semi-Automated Pipeline Supported by Expert Guidance. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.10.31.621187. [PMID: 39554050 PMCID: PMC11565955 DOI: 10.1101/2024.10.31.621187] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Subscribe] [Scholar Register] [Indexed: 11/19/2024]
Abstract
Asymmetry is a key feature of numerous developmental disorders and in phenotypic screens is often used as a readout for environmental or genetic perturbations to normal development. A better understanding of the genetic basis of asymmetry and its relationship to disease susceptibility will help unravel the complex genetic and environmental factors and their interactions that increase risk in a range of developmental disorders. Large-scale imaging datasets offer opportunities to work with sample sizes needed to detect and quantify differences in morphology beyond severe deformities while also posing challenges to manual phenotyping protocols. In this work, we introduce a semi-automated open-source workflow to quantify abnormal asymmetry of craniofacial structures that integrates expert anatomical knowledge. We apply this workflow to explore the role of genes contributing to abnormal asymmetry by deep phenotyping 3D fetal microCT images from knockout strains acquired as part of the Knockout Mouse Phenotyping Program (KOMP2). Four knockout strains: Ccdc186, Acvr2a, Nhlh1, and Fam20c were identified with highly significant asymmetry in craniofacial regions, making them good candidates for further analysis into their potential roles in asymmetry and developmental disorders.
Collapse
Affiliation(s)
- S M Rolfe
- Center for Developmental Biology and Regenerative Medicine, Seattle Children's Research Institute, Seattle, WA, USA
| | - D Mao
- Department of Pediatrics, University of Washington, Seattle, WA, USA
| | - A M Maga
- Center for Developmental Biology and Regenerative Medicine, Seattle Children's Research Institute, Seattle, WA, USA
- Department of Pediatrics, University of Washington, Seattle, WA, USA
| |
Collapse
|
3
|
He Y, Mulqueeney JM, Watt EC, Salili-James A, Barber NS, Camaiti M, Hunt ESE, Kippax-Chui O, Knapp A, Lanzetti A, Rangel-de Lázaro G, McMinn JK, Minus J, Mohan AV, Roberts LE, Adhami D, Grisan E, Gu Q, Herridge V, Poon STS, West T, Goswami A. Opportunities and Challenges in Applying AI to Evolutionary Morphology. Integr Org Biol 2024; 6:obae036. [PMID: 40433986 PMCID: PMC12082097 DOI: 10.1093/iob/obae036] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2024] [Revised: 08/07/2024] [Accepted: 09/20/2024] [Indexed: 05/29/2025] Open
Abstract
Artificial intelligence (AI) is poised to revolutionize many aspects of science, including the study of evolutionary morphology. While classical AI methods such as principal component analysis and cluster analysis have been commonplace in the study of evolutionary morphology for decades, recent years have seen increasing application of deep learning to ecology and evolutionary biology. As digitized specimen databases become increasingly prevalent and openly available, AI is offering vast new potential to circumvent long-standing barriers to rapid, big data analysis of phenotypes. Here, we review the current state of AI methods available for the study of evolutionary morphology, which are most developed in the area of data acquisition and processing. We introduce the main available AI techniques, categorizing them into 3 stages based on their order of appearance: (1) machine learning, (2) deep learning, and (3) the most recent advancements in large-scale models and multimodal learning. Next, we present case studies of existing approaches using AI for evolutionary morphology, including image capture and segmentation, feature recognition, morphometrics, and phylogenetics. We then discuss the prospectus for near-term advances in specific areas of inquiry within this field, including the potential of new AI methods that have not yet been applied to the study of morphological evolution. In particular, we note key areas where AI remains underutilized and could be used to enhance studies of evolutionary morphology. This combination of current methods and potential developments has the capacity to transform the evolutionary analysis of the organismal phenotype into evolutionary phenomics, leading to an era of "big data" that aligns the study of phenotypes with genomics and other areas of bioinformatics.
Collapse
Affiliation(s)
- Y He
- Life Sciences, Natural History Museum, London, UK
| | - J M Mulqueeney
- Life Sciences, Natural History Museum, London, UK
- Department of Ocean & Earth Science, National Oceanography Centre Southampton, University of Southampton, Southampton, UK
| | - E C Watt
- Life Sciences, Natural History Museum, London, UK
- Division of Biosciences, University College London, London, UK
| | - A Salili-James
- AI and Innovation, Natural History Museum, London, UK
- Digital, Data and Informatics, Natural History Museum, London, UK
| | - N S Barber
- Life Sciences, Natural History Museum, London, UK
- Department of Anthropology, University College London, London, UK
| | - M Camaiti
- Life Sciences, Natural History Museum, London, UK
| | - E S E Hunt
- Life Sciences, Natural History Museum, London, UK
- Department of Life Sciences, Imperial College London, London, UK
- Grantham Institute, Imperial College London, London, UK
| | - O Kippax-Chui
- Life Sciences, Natural History Museum, London, UK
- Grantham Institute, Imperial College London, London, UK
- Department of Earth Science and Engineering, Imperial College London, London, UK
| | - A Knapp
- Life Sciences, Natural History Museum, London, UK
- Centre for Integrative Anatomy, University College London, London, UK
| | - A Lanzetti
- Life Sciences, Natural History Museum, London, UK
- School of Geography, Earth and Environmental Sciences, University of Birmingham, Birmingham, UK
| | - G Rangel-de Lázaro
- Life Sciences, Natural History Museum, London, UK
- School of Oriental and African Studies, London, UK
| | - J K McMinn
- Life Sciences, Natural History Museum, London, UK
- Department of Earth Sciences, University of Oxford, Oxford, UK
| | - J Minus
- Life Sciences, Natural History Museum, London, UK
- School of Biological and Behavioural Sciences, Queen Mary University of London, London, UK
| | - A V Mohan
- Life Sciences, Natural History Museum, London, UK
- Biodiversity Genomics Laboratory, Institute of Biology, University of Neuchâtel, Neuchâtel, Switzerland
| | - L E Roberts
- Life Sciences, Natural History Museum, London, UK
| | - D Adhami
- Life Sciences, Natural History Museum, London, UK
- Department of Life Sciences, Imperial College London, London, UK
- Imaging and Analysis Centre, Natural History Museum, London, UK
| | - E Grisan
- School of Engineering, London South Bank University, London, UK
| | - Q Gu
- AI and Innovation, Natural History Museum, London, UK
- Digital, Data and Informatics, Natural History Museum, London, UK
| | - V Herridge
- Life Sciences, Natural History Museum, London, UK
- School of Biosciences, University of Sheffield, Sheffield, UK
| | - S T S Poon
- AI and Innovation, Natural History Museum, London, UK
- Digital, Data and Informatics, Natural History Museum, London, UK
| | - T West
- Centre for Integrative Anatomy, University College London, London, UK
- Imaging and Analysis Centre, Natural History Museum, London, UK
| | - A Goswami
- Life Sciences, Natural History Museum, London, UK
| |
Collapse
|
4
|
Manuta N, Çakar B, Gündemir O, Spataru MC. Shape and Size Variations of Distal Phalanges in Cattle. Animals (Basel) 2024; 14:194. [PMID: 38254363 PMCID: PMC10812660 DOI: 10.3390/ani14020194] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2023] [Revised: 01/04/2024] [Accepted: 01/05/2024] [Indexed: 01/24/2024] Open
Abstract
Studies on the structure of the distal phalanx help explain the development of laminitis. Additionally, examining the structure of the distal phalanx from a taxonomic perspective also contributes to veterinary anatomy. In this study, we examined shape variation in the medial and lateral distal phalanx of both fore- and hindlimbs using the geometric morphometry method. We investigated whether the shape of the distal phalanx differed between phalanx positions and how much of the shape variation in this bone depends on size. For this purpose, distal phalanges from 20 Holstein cattle were used, and the bones were digitized in 3D. A draft containing 176 semi-landmarks was prepared for shape analysis, and this draft was applied to all samples using automated landmarking through point cloud alignment and correspondence analysis. A principal component analysis was performed to obtain general patterns of morphological variation. The centroid size (CS) was employed as an approximation of size. Although distal phalanx groups generally showed close variations, PC1 statistically separated the hindlimb lateral distal phalanx (HL) and the forelimb medial distal phalanx (FM) from each other in shape. While PC2 separated HL from other distal phalanx groups, PC3 separated fore- and hindlimb groups. The shape (Procrustes distance) of the hindlimb medial distal phalanx (HM) is markedly less variable than the other three phalanges. The smallest distal phalanx in size was HL. For both forelimb and hindlimb, the medial distal phalanges were larger than the lateral ones. Size (CS) was found to have an effect on PC1 and PC3. In this study, a reference model of the same breeds for distal phalanx was created. These results can provide useful information, especially in terms of veterinary anatomy, zooarchaeology, and paleontology.
Collapse
Affiliation(s)
- Nicoleta Manuta
- Institute of Graduate Studies, Istanbul University-Cerrahpasa, Istanbul 34320, Türkiye; (N.M.); (B.Ç.)
| | - Buket Çakar
- Institute of Graduate Studies, Istanbul University-Cerrahpasa, Istanbul 34320, Türkiye; (N.M.); (B.Ç.)
| | - Ozan Gündemir
- Department of Anatomy, Faculty of Veterinary Medicine, Istanbul University-Cerrahpasa, Istanbul 34320, Türkiye
| | - Mihaela-Claudia Spataru
- Department of Public Health, Faculty of Veterinary Medicine, Iasi University of Life Sciences, 700490 Iasi, Romania;
| |
Collapse
|
5
|
Korff A, Yang X, O’Donovan K, Gonzalez A, Teubner BJ, Nakamura H, Messing J, Yang F, Carisey AF, Wang YD, Patni T, Sheppard H, Zakharenko SS, Chook YM, Taylor JP, Kim HJ. A murine model of hnRNPH2-related neurodevelopmental disorder reveals a mechanism for genetic compensation by Hnrnph1. J Clin Invest 2023; 133:e160309. [PMID: 37463454 PMCID: PMC10348767 DOI: 10.1172/jci160309] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/23/2022] [Accepted: 05/24/2023] [Indexed: 07/20/2023] Open
Abstract
Mutations in HNRNPH2 cause an X-linked neurodevelopmental disorder with features that include developmental delay, motor function deficits, and seizures. More than 90% of patients with hnRNPH2 have a missense mutation within or adjacent to the nuclear localization signal (NLS) of hnRNPH2. Here, we report that hnRNPH2 NLS mutations caused reduced interaction with the nuclear transport receptor Kapβ2 and resulted in modest cytoplasmic accumulation of hnRNPH2. We generated 2 knockin mouse models with human-equivalent mutations in Hnrnph2 as well as Hnrnph2-KO mice. Knockin mice recapitulated clinical features of the human disorder, including reduced survival in male mice, impaired motor and cognitive functions, and increased susceptibility to audiogenic seizures. In contrast, 2 independent lines of Hnrnph2-KO mice showed no detectable phenotypes. Notably, KO mice had upregulated expression of Hnrnph1, a paralog of Hnrnph2, whereas knockin mice failed to upregulate Hnrnph1. Thus, genetic compensation by Hnrnph1 may counteract the loss of hnRNPH2. These findings suggest that HNRNPH2-related disorder may be driven by a toxic gain of function or a complex loss of HNRNPH2 function with impaired compensation by HNRNPH1. The knockin mice described here are an important resource for preclinical studies to assess the therapeutic benefit of gene replacement or knockdown of mutant hnRNPH2.
Collapse
Affiliation(s)
- Ane Korff
- Department of Cell and Molecular Biology, St. Jude Children’s Research Hospital, Memphis, Tennessee, USA
| | - Xiaojing Yang
- Department of Cell and Molecular Biology, St. Jude Children’s Research Hospital, Memphis, Tennessee, USA
| | - Kevin O’Donovan
- Department of Cell and Molecular Biology, St. Jude Children’s Research Hospital, Memphis, Tennessee, USA
| | - Abner Gonzalez
- Department of Pharmacology, University of Texas Southwestern Medical Center, Dallas, Texas, USA
| | | | - Haruko Nakamura
- Department of Cell and Molecular Biology, St. Jude Children’s Research Hospital, Memphis, Tennessee, USA
| | - James Messing
- Department of Cell and Molecular Biology, St. Jude Children’s Research Hospital, Memphis, Tennessee, USA
| | - Fen Yang
- Department of Cell and Molecular Biology, St. Jude Children’s Research Hospital, Memphis, Tennessee, USA
| | - Alexandre F. Carisey
- Department of Cell and Molecular Biology, St. Jude Children’s Research Hospital, Memphis, Tennessee, USA
| | - Yong-Dong Wang
- Department of Cell and Molecular Biology, St. Jude Children’s Research Hospital, Memphis, Tennessee, USA
| | | | - Heather Sheppard
- Department of Pathology, St. Jude Children’s Research Hospital, Memphis, Tennessee, USA
| | | | - Yuh Min Chook
- Department of Pharmacology, University of Texas Southwestern Medical Center, Dallas, Texas, USA
| | - J. Paul Taylor
- Department of Cell and Molecular Biology, St. Jude Children’s Research Hospital, Memphis, Tennessee, USA
- Howard Hughes Medical Institute, Chevy Chase, Maryland, USA
| | - Hong Joo Kim
- Department of Cell and Molecular Biology, St. Jude Children’s Research Hospital, Memphis, Tennessee, USA
| |
Collapse
|
6
|
Hegdé J, Tustison NJ, Parker WT, Branch F, Yanasak N, Stumpo LA. An Anatomical Template for the Normalization of Medical Images of Adult Human Hands. Diagnostics (Basel) 2023; 13:2010. [PMID: 37370905 DOI: 10.3390/diagnostics13122010] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2022] [Revised: 06/01/2023] [Accepted: 06/02/2023] [Indexed: 06/29/2023] Open
Abstract
During medical image analysis, it is often useful to align (or 'normalize') a given image of a given body part to a representative standard (or 'template') of that body part. The impact that brain templates have had on the analysis of brain images highlights the importance of templates in general. However, templates for human hands do not exist. Image normalization is especially important for hand images because hands, by design, readily change shape during various tasks. Here we report the construction of an anatomical template for healthy adult human hands. To do this, we used 27 anatomically representative T1-weighted magnetic resonance (MR) images of either hand from 21 demographically representative healthy adult subjects (13 females and 8 males). We used the open-source, cross-platform ANTs (Advanced Normalization Tools) medical image analysis software framework, to preprocess the MR images. The template was constructed using the ANTs standard multivariate template construction workflow. The resulting template image preserved all the essential anatomical features of the hand, including all the individual bones, muscles, tendons, ligaments, as well as the main branches of the median nerve and radial, ulnar, and palmar metacarpal arteries. Furthermore, the image quality of the template was significantly higher than that of the underlying individual hand images as measured by two independent canonical metrics of image quality.
Collapse
Affiliation(s)
- Jay Hegdé
- Department of Neuroscience and Regenerative Medicine, Medical College of Georgia, Augusta University, Augusta, GA 30912, USA
| | - Nicholas J Tustison
- Department of Radiology and Medical Imaging, School of Medicine, University of Virginia, Charlottesville, VA 22903, USA
| | - William T Parker
- Department of Radiology and Imaging, Medical College of Georgia, Augusta University, Augusta, GA 30912, USA
| | - Fallon Branch
- Department of Neuroscience and Regenerative Medicine, Medical College of Georgia, Augusta University, Augusta, GA 30912, USA
| | - Nathan Yanasak
- Department of Radiology and Imaging, Medical College of Georgia, Augusta University, Augusta, GA 30912, USA
| | - Lorie A Stumpo
- Department of Radiology and Imaging, Medical College of Georgia, Augusta University, Augusta, GA 30912, USA
| |
Collapse
|
7
|
Fournier G, Maret D, Telmon N, Savall F. An automated landmark method to describe geometric changes in the human mandible during growth. Arch Oral Biol 2023; 149:105663. [PMID: 36893681 DOI: 10.1016/j.archoralbio.2023.105663] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/23/2022] [Revised: 01/27/2023] [Accepted: 02/22/2023] [Indexed: 02/26/2023]
Abstract
OBJECTIVE The principal aim of this study was to assess an automatic landmarking approach to human mandibles based on the atlas method. The secondary aim was to identify the areas of greatest variation in the mandibles of middle-aged to older adults. DESIGN Our sample consisted of 160 mandibles from computed tomography scans of 80 men and 80 women aged between 40 and 79 years. Eleven anatomical landmarks were placed manually on mandibles. The automated landmarking through point cloud alignment and correspondence (ALPACA) method implemented in 3D Slicer was used to automatically place landmarks to all meshes. Euclidean distances, normalized centroid size, and Procrustes ANOVA were calculated for both methods. A pseudo-landmarks approach was followed using ALPACA to identify areas of changes among our sample. RESULTS The ALPACA method showed significant differences in Euclidean distances for all landmarks compared to the manual method. A mean Euclidean distance of 1.7 mm was found for the ALPACA method and 0.99 mm for the manual method. Both methods found that sex, age, and size had a significant effect on mandibular shape. The greatest variations were observed in the condyle, ramus, and symphysis regions. CONCLUSION The results obtained using the ALPACA method are acceptable and promising. This approach can automatically place landmarks with an average accuracy of less than 2 mm, which may be sufficient in most anthropometric analyses. In the light of our results, however, odontological application such as occlusal analysis is not recommended.
Collapse
Affiliation(s)
- G Fournier
- Faculté de Chirurgie Dentaire, Université Paul Sabatier, Centre Hospitalier Universitaire, Toulouse, France; Laboratory Centre for Anthropology and Genomics of Toulouse, Université Paul Sabatier, Toulouse, France.
| | - D Maret
- Faculté de Chirurgie Dentaire, Université Paul Sabatier, Centre Hospitalier Universitaire, Toulouse, France; Laboratory Centre for Anthropology and Genomics of Toulouse, Université Paul Sabatier, Toulouse, France
| | - N Telmon
- Laboratory Centre for Anthropology and Genomics of Toulouse, Université Paul Sabatier, Toulouse, France; Service de Médecine Légale, Hôpital de Rangueil, Toulouse, France
| | - F Savall
- Laboratory Centre for Anthropology and Genomics of Toulouse, Université Paul Sabatier, Toulouse, France; Service de Médecine Légale, Hôpital de Rangueil, Toulouse, France
| |
Collapse
|
8
|
Automated morphological phenotyping using learned shape descriptors and functional maps: A novel approach to geometric morphometrics. PLoS Comput Biol 2023; 19:e1009061. [PMID: 36656910 PMCID: PMC9970057 DOI: 10.1371/journal.pcbi.1009061] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2021] [Revised: 02/27/2023] [Accepted: 11/13/2022] [Indexed: 01/20/2023] Open
Abstract
The methods of geometric morphometrics are commonly used to quantify morphology in a broad range of biological sciences. The application of these methods to large datasets is constrained by manual landmark placement limiting the number of landmarks and introducing observer bias. To move the field forward, we need to automate morphological phenotyping in ways that capture comprehensive representations of morphological variation with minimal observer bias. Here, we present Morphological Variation Quantifier (morphVQ), a shape analysis pipeline for quantifying, analyzing, and exploring shape variation in the functional domain. morphVQ uses descriptor learning to estimate the functional correspondence between whole triangular meshes in lieu of landmark configurations. With functional maps between pairs of specimens in a dataset we can analyze and explore shape variation. morphVQ uses Consistent ZoomOut refinement to improve these functional maps and produce a new representation of shape variation, area-based and conformal (angular) latent shape space differences (LSSDs). We compare this new representation of shape variation to shape variables obtained via manual digitization and auto3DGM, an existing approach to automated morphological phenotyping. We find that LSSDs compare favorably to modern 3DGM and auto3DGM while being more computationally efficient. By characterizing whole surfaces, our method incorporates more morphological detail in shape analysis. We can classify known biological groupings, such as Genus affiliation with comparable accuracy. The shape spaces produced by our method are similar to those produced by modern 3DGM and to auto3DGM, and distinctiveness functions derived from LSSDs show us how shape variation differs between groups. morphVQ can capture shape in an automated fashion while avoiding the limitations of manually digitized landmarks, and thus represents a novel and computationally efficient addition to the geometric morphometrics toolkit.
Collapse
|
9
|
Zhang C, Porto A, Rolfe S, Kocatulum A, Maga AM. Automated landmarking via multiple templates. PLoS One 2022; 17:e0278035. [PMID: 36454982 PMCID: PMC9714854 DOI: 10.1371/journal.pone.0278035] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2022] [Accepted: 11/08/2022] [Indexed: 12/02/2022] Open
Abstract
Manually collecting landmarks for quantifying complex morphological phenotypes can be laborious and subject to intra and interobserver errors. However, most automated landmarking methods for efficiency and consistency fall short of landmarking highly variable samples due to the bias introduced by the use of a single template. We introduce a fast and open source automated landmarking pipeline (MALPACA) that utilizes multiple templates for accommodating large-scale variations. We also introduce a K-means method of choosing the templates that can be used in conjunction with MALPACA, when no prior information for selecting templates is available. Our results confirm that MALPACA significantly outperforms single-template methods in landmarking both single and multi-species samples. K-means based template selection can also avoid choosing the worst set of templates when compared to random template selection. We further offer an example of post-hoc quality check for each individual template for further refinement. In summary, MALPACA is an efficient and reproducible method that can accommodate large morphological variability, such as those commonly found in evolutionary studies. To support the research community, we have developed open-source and user-friendly software tools for performing K-means multi-templates selection and MALPACA.
Collapse
Affiliation(s)
- Chi Zhang
- Center for Development Biology and Regenerative Medicine, Seattle Children’s Research Institute, Seattle, Washington, United States of America
| | - Arthur Porto
- Department of Biological Sciences, Louisiana State University, Baton Rouge, Louisiana, United States of America
- Center for Computation and Technology, Louisiana State University, Baton Rouge, Louisiana, United States of America
| | - Sara Rolfe
- Center for Development Biology and Regenerative Medicine, Seattle Children’s Research Institute, Seattle, Washington, United States of America
- Friday Harbor Laboratories, University of Washington, San Juan Island, Washington, United States of America
| | - Altan Kocatulum
- Alfred University, Alfred, New York, United States of America
| | - A. Murat Maga
- Center for Development Biology and Regenerative Medicine, Seattle Children’s Research Institute, Seattle, Washington, United States of America
- Division of Craniofacial Medicine, Department of Pediatrics, University of Washington, Seattle, Washington, United States of America
- * E-mail:
| |
Collapse
|
10
|
Devine J, Vidal-García M, Liu W, Neves A, Lo Vercio LD, Green RM, Richbourg HA, Marchini M, Unger CM, Nickle AC, Radford B, Young NM, Gonzalez PN, Schuler RE, Bugacov A, Rolian C, Percival CJ, Williams T, Niswander L, Calof AL, Lander AD, Visel A, Jirik FR, Cheverud JM, Klein OD, Birnbaum RY, Merrill AE, Ackermann RR, Graf D, Hemberger M, Dean W, Forkert ND, Murray SA, Westerberg H, Marcucio RS, Hallgrímsson B. MusMorph, a database of standardized mouse morphology data for morphometric meta-analyses. Sci Data 2022; 9:230. [PMID: 35614082 PMCID: PMC9133120 DOI: 10.1038/s41597-022-01338-x] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/15/2021] [Accepted: 04/13/2022] [Indexed: 11/08/2022] Open
Abstract
Complex morphological traits are the product of many genes with transient or lasting developmental effects that interact in anatomical context. Mouse models are a key resource for disentangling such effects, because they offer myriad tools for manipulating the genome in a controlled environment. Unfortunately, phenotypic data are often obtained using laboratory-specific protocols, resulting in self-contained datasets that are difficult to relate to one another for larger scale analyses. To enable meta-analyses of morphological variation, particularly in the craniofacial complex and brain, we created MusMorph, a database of standardized mouse morphology data spanning numerous genotypes and developmental stages, including E10.5, E11.5, E14.5, E15.5, E18.5, and adulthood. To standardize data collection, we implemented an atlas-based phenotyping pipeline that combines techniques from image registration, deep learning, and morphometrics. Alongside stage-specific atlases, we provide aligned micro-computed tomography images, dense anatomical landmarks, and segmentations (if available) for each specimen (N = 10,056). Our workflow is open-source to encourage transparency and reproducible data collection. The MusMorph data and scripts are available on FaceBase ( www.facebase.org , https://doi.org/10.25550/3-HXMC ) and GitHub ( https://github.com/jaydevine/MusMorph ).
Collapse
Affiliation(s)
- Jay Devine
- Alberta Children's Hospital Research Institute, University of Calgary, 28 Oki Dr NW, Calgary, AB, T3B 6A8, Canada
- The McCaig Institute for Bone and Joint Health, University of Calgary, 3280 Hospital Dr NW, Calgary, AB, T2N 4Z6, Canada
- Department of Cell Biology and Anatomy, Cumming School of Medicine, University of Calgary, 3330 Hospital Dr NW, Calgary, AB, T2N 4N1, Canada
| | - Marta Vidal-García
- Alberta Children's Hospital Research Institute, University of Calgary, 28 Oki Dr NW, Calgary, AB, T3B 6A8, Canada
- The McCaig Institute for Bone and Joint Health, University of Calgary, 3280 Hospital Dr NW, Calgary, AB, T2N 4Z6, Canada
- Department of Cell Biology and Anatomy, Cumming School of Medicine, University of Calgary, 3330 Hospital Dr NW, Calgary, AB, T2N 4N1, Canada
| | - Wei Liu
- Alberta Children's Hospital Research Institute, University of Calgary, 28 Oki Dr NW, Calgary, AB, T3B 6A8, Canada
- The McCaig Institute for Bone and Joint Health, University of Calgary, 3280 Hospital Dr NW, Calgary, AB, T2N 4Z6, Canada
- Department of Cell Biology and Anatomy, Cumming School of Medicine, University of Calgary, 3330 Hospital Dr NW, Calgary, AB, T2N 4N1, Canada
| | - Amanda Neves
- Department of Biology, McMaster University, 1280 Main St W, Hamilton, ON, L8S 4L8, Canada
| | - Lucas D Lo Vercio
- Alberta Children's Hospital Research Institute, University of Calgary, 28 Oki Dr NW, Calgary, AB, T3B 6A8, Canada
- The McCaig Institute for Bone and Joint Health, University of Calgary, 3280 Hospital Dr NW, Calgary, AB, T2N 4Z6, Canada
- Department of Cell Biology and Anatomy, Cumming School of Medicine, University of Calgary, 3330 Hospital Dr NW, Calgary, AB, T2N 4N1, Canada
| | - Rebecca M Green
- School of Dental Medicine, University of Pittsburgh, 3501 Terrace St, Pittsburgh, PA, 15213, USA
| | - Heather A Richbourg
- Orthopaedic Trauma Institute, ZSFG, UCSF, 2550 23rd St, San Francisco, CA, 94110, USA
| | - Marta Marchini
- Alberta Children's Hospital Research Institute, University of Calgary, 28 Oki Dr NW, Calgary, AB, T3B 6A8, Canada
- The McCaig Institute for Bone and Joint Health, University of Calgary, 3280 Hospital Dr NW, Calgary, AB, T2N 4Z6, Canada
- Department of Cell Biology and Anatomy, Cumming School of Medicine, University of Calgary, 3330 Hospital Dr NW, Calgary, AB, T2N 4N1, Canada
| | - Colton M Unger
- The McCaig Institute for Bone and Joint Health, University of Calgary, 3280 Hospital Dr NW, Calgary, AB, T2N 4Z6, Canada
- Department of Biological Sciences, University of Calgary, 2500 University Dr NW, Calgary, AB, T2N 1N4, Canada
| | - Audrey C Nickle
- Center for Craniofacial Molecular Biology, Department of Biomedical Sciences, Herman Ostrow School of Dentistry, University of Southern California, Los Angeles, 2250 Alcazar St, Los Angeles, CA, 90033, USA
- Department of Biochemistry and Molecular Medicine, Keck School of Medicine, University of Southern California, Los Angeles, 1975 Zonal Ave, Los Angeles, CA, 90033, USA
| | - Bethany Radford
- Alberta Children's Hospital Research Institute, University of Calgary, 28 Oki Dr NW, Calgary, AB, T3B 6A8, Canada
- Department of Biochemistry and Molecular Biology, Cumming School of Medicine, University of Calgary, 3330 Hospital Dr NW, Calgary, AB, T2N 4N1, Canada
| | - Nathan M Young
- Orthopaedic Trauma Institute, ZSFG, UCSF, 2550 23rd St, San Francisco, CA, 94110, USA
| | - Paula N Gonzalez
- Institute for Studies in Neuroscience and Complex Systems (ENyS) CONICET, Av. Calchaquí, 5402, Florencio Varela, Buenos Aires, Argentina
| | - Robert E Schuler
- Information Sciences Institute, Viterbi School of Engineering, University of Southern California, 4676 Admiralty Way, Marina del Rey, CA, 90292, USA
| | - Alejandro Bugacov
- Information Sciences Institute, Viterbi School of Engineering, University of Southern California, 4676 Admiralty Way, Marina del Rey, CA, 90292, USA
| | - Campbell Rolian
- The McCaig Institute for Bone and Joint Health, University of Calgary, 3280 Hospital Dr NW, Calgary, AB, T2N 4Z6, Canada
- Department of Comparative Biology and Experimental Medicine, Faculty of Veterinary Medicine, University of Calgary, 3330 Hospital Dr NW, Calgary, AB, T2N 4N1, Canada
| | - Christopher J Percival
- Department of Anthropology, Stony Brook University, 100 Nicolls Rd, Stony Brook, NY, 11794, USA
| | - Trevor Williams
- Department of Craniofacial Biology, University of Colorado Anschutz Medical Campus, 12801 East 17th Ave, Aurora, CO, 80045, USA
| | - Lee Niswander
- Department of Molecular, Cellular and Developmental Biology, University of Colorado Boulder, Boulder, CO, 80309, USA
| | - Anne L Calof
- Department of Anatomy and Neurobiology, University of California, Irvine, Irvine, CA, 92697, USA
- Center for Complex Biological Systems, University of California, Irvine, Irvine, CA, 92697, USA
| | - Arthur D Lander
- Center for Complex Biological Systems, University of California, Irvine, Irvine, CA, 92697, USA
- Department of Developmental and Cell Biology, University of California, Irvine, Irvine, CA, 92697, USA
| | - Axel Visel
- Environmental Genomics and Systems Biology Division, Lawrence Berkeley National Laboratory, 1 Cyclotron Rd, Berkeley, CA, 94720, USA
- U.S. Department of Energy Joint Genome Institute, Lawrence Berkeley National Laboratory, 1 Cyclotron Rd, Berkeley, CA, 94720, USA
- School of Natural Sciences, University of California, Merced, 5200 Lake Rd, Merced, CA, 95343, USA
| | - Frank R Jirik
- Alberta Children's Hospital Research Institute, University of Calgary, 28 Oki Dr NW, Calgary, AB, T3B 6A8, Canada
- The McCaig Institute for Bone and Joint Health, University of Calgary, 3280 Hospital Dr NW, Calgary, AB, T2N 4Z6, Canada
- Department of Biochemistry and Molecular Biology, Cumming School of Medicine, University of Calgary, 3330 Hospital Dr NW, Calgary, AB, T2N 4N1, Canada
| | - James M Cheverud
- Department of Biology, Loyola University Chicago, 1032 W Sheridan Rd, Chicago, IL, 60660, USA
| | - Ophir D Klein
- Department of Orofacial Sciences and Program in Craniofacial Biology, University of California, San Francisco, 513 Parnassus Ave, San Francisco, CA, 94143, USA
- Department of Pediatrics and Institute for Human Genetics, University of California, San Francisco, 513 Parnassus Ave, San Francisco, CA, 94143, USA
- Department of Pediatrics, Cedars-Sinai Medical Center, 8700 Beverly Blvd, Los Angeles, CA, 90048, USA
| | - Ramon Y Birnbaum
- Department of Life Sciences, Faculty of Natural Sciences, The Ben-Gurion University of the Negev, David Ben Gurion Blvd 1, Be'er Sheva, Israel
| | - Amy E Merrill
- Center for Craniofacial Molecular Biology, Department of Biomedical Sciences, Herman Ostrow School of Dentistry, University of Southern California, Los Angeles, 2250 Alcazar St, Los Angeles, CA, 90033, USA
- Department of Biochemistry and Molecular Medicine, Keck School of Medicine, University of Southern California, Los Angeles, 1975 Zonal Ave, Los Angeles, CA, 90033, USA
| | - Rebecca R Ackermann
- Department of Archaeology, University of Cape Town, Rondebosch, Cape Town, 7700, South Africa
- Human Evolution Research Institute, University of Cape Town, Rondebosch, Cape Town, 7700, South Africa
| | - Daniel Graf
- School of Dentistry, Faculty of Medicine and Dentistry, University of Alberta, 116 St. and 85 Ave, Edmonton, AB, T6G 2R3, Canada
- Department of Medical Genetics, Faculty of Medicine and Dentistry, University of Alberta, 116 St. and 85 Ave, Edmonton, AB, T6G 2R3, Canada
| | - Myriam Hemberger
- Alberta Children's Hospital Research Institute, University of Calgary, 28 Oki Dr NW, Calgary, AB, T3B 6A8, Canada
- Department of Biochemistry and Molecular Biology, Cumming School of Medicine, University of Calgary, 3330 Hospital Dr NW, Calgary, AB, T2N 4N1, Canada
| | - Wendy Dean
- Alberta Children's Hospital Research Institute, University of Calgary, 28 Oki Dr NW, Calgary, AB, T3B 6A8, Canada
- Department of Cell Biology and Anatomy, Cumming School of Medicine, University of Calgary, 3330 Hospital Dr NW, Calgary, AB, T2N 4N1, Canada
| | - Nils D Forkert
- Alberta Children's Hospital Research Institute, University of Calgary, 28 Oki Dr NW, Calgary, AB, T3B 6A8, Canada
- Department of Radiology, Cumming School of Medicine, University of Calgary, 3330 Hospital Dr NW, Calgary, AB, T2N 4N1, Canada
| | | | - Henrik Westerberg
- Department of Bioimaging Informatics, MRC Harwell Institute, Oxfordshire, OX11 0RD, UK
| | - Ralph S Marcucio
- Orthopaedic Trauma Institute, ZSFG, UCSF, 2550 23rd St, San Francisco, CA, 94110, USA
| | - Benedikt Hallgrímsson
- Alberta Children's Hospital Research Institute, University of Calgary, 28 Oki Dr NW, Calgary, AB, T3B 6A8, Canada.
- The McCaig Institute for Bone and Joint Health, University of Calgary, 3280 Hospital Dr NW, Calgary, AB, T2N 4Z6, Canada.
- Department of Cell Biology and Anatomy, Cumming School of Medicine, University of Calgary, 3330 Hospital Dr NW, Calgary, AB, T2N 4N1, Canada.
| |
Collapse
|
11
|
Harper CM, Goldstein DM, Sylvester AD. Comparing and combining sliding semilandmarks and weighted spherical harmonics for shape analysis. J Anat 2022; 240:678-687. [PMID: 34747020 PMCID: PMC8930823 DOI: 10.1111/joa.13589] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/25/2021] [Revised: 10/27/2021] [Accepted: 10/29/2021] [Indexed: 11/28/2022] Open
Abstract
Quantifying morphological variation is critical for conducting anatomical research. Three-dimensional geometric morphometric (3D GM) landmark analyses quantify shape using homologous Cartesian coordinates (landmarks). Setting up a high-density landmark set and placing it on all specimens, however, can be a time-consuming task. Weighted spherical harmonics (SPHARM) provides an alternative method for analyzing the shape of such objects. Here we compare sliding semilandmark and SPHARM analyses of the calcaneus of Gorilla gorilla gorilla (n = 20), Pan troglodytes troglodytes (n = 20), and Homo sapiens (n = 20) to determine whether the SPHARM and sliding semilandmark analyses capture comparable levels of shape variation. We also compare both the sliding semilandmark and SPHARM analyses to a novel combination of the two methods, here termed SPHARM-sliding. In SPHARM-sliding, the vertices of the surface models produced from the SPHARM analysis (that are the same in number and relative location) are used as the starting landmark positions for a sliding semilandmark analysis. Calcaneal shape variation quantified by all three analyses was summarized using separate principal components analyses. Results were compared using the root mean square (RMS) and maximum distance between surface models of species averages scaled (up) to centroid size created from each analysis. The average RMS was 0.23 mm between sliding semilandmark and SPHARM average surface models, 0.19 mm between SPHARM and SPHARM sliding average surface models, and 0.22 mm between sliding semilandmark and SPHARM sliding average surface models. Although results indicate that all three analyses are comparable methods for 3D shape analysis, there are advantages and disadvantages to each. While the SPHARM analysis is less time-intensive, it is unable to capture the same level of detail around the sharp edges of articular facets on average surface models as the sliding semilandmark analysis. The SPHARM analysis also does not allow for individual articular facets to be analyzed in isolation. SPHARM-sliding, however, captures the same level of detail as the sliding semilandmark analysis, and (as in the sliding semilandmark analysis) allows for the evaluation of individual portions of bone. SPHARM is a comparable method to a 3D GM analysis for small, irregularly shaped bones, such as the calcaneus, and SPHARM-sliding allows for an expedited set up process for a sliding semilandmark analysis.
Collapse
Affiliation(s)
- Christine M. Harper
- Department of Biomedical SciencesCooper Medical School of Rowan UniversityCamdenNew JerseyUSA
- Center for Functional Anatomy and EvolutionThe Johns Hopkins University School of MedicineBaltimoreMarylandUSA
| | - Deanna M. Goldstein
- Center for Functional Anatomy and EvolutionThe Johns Hopkins University School of MedicineBaltimoreMarylandUSA
| | - Adam D. Sylvester
- Center for Functional Anatomy and EvolutionThe Johns Hopkins University School of MedicineBaltimoreMarylandUSA
| |
Collapse
|
12
|
Diamond KM, Rolfe SM, Kwon RY, Maga AM. Computational anatomy and geometric shape analysis enables analysis of complex craniofacial phenotypes in zebrafish. Biol Open 2022; 11:bio058948. [PMID: 35072203 PMCID: PMC8864294 DOI: 10.1242/bio.058948] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/22/2021] [Accepted: 01/11/2022] [Indexed: 11/20/2022] Open
Abstract
Due to the complexity of fish skulls, previous attempts to classify craniofacial phenotypes have relied on qualitative features or sparce 2D landmarks. In this work we aim to identify previously unknown 3D craniofacial phenotypes with a semiautomated pipeline in adult zebrafish mutants. We first estimate a synthetic 'normative' zebrafish template using MicroCT scans from a sample pool of wild-type animals using the Advanced Normalization Tools (ANTs). We apply a computational anatomy (CA) approach to quantify the phenotype of zebrafish with disruptions in bmp1a, a gene implicated in later skeletal development and whose human ortholog when disrupted is associated with Osteogenesis Imperfecta. Compared to controls, the bmp1a fish have larger otoliths, larger normalized centroid sizes, and exhibit shape differences concentrated around the operculum, anterior frontal, and posterior parietal bones. Moreover, bmp1a fish differ in the degree of asymmetry. Our CA approach offers a potential pipeline for high-throughput screening of complex fish craniofacial shape to discover novel phenotypes for which traditional landmarks are too sparce to detect. The current pipeline successfully identifies areas of variation in zebrafish mutants, which are an important model system for testing genome to phenome relationships in the study of development, evolution, and human diseases. This article has an associated First Person interview with the first author of the paper.
Collapse
Affiliation(s)
- Kelly M. Diamond
- Center for Developmental Biology and Regenerative Medicine, Seattle Children's Research Institute, Seattle, WA 98101, USA
| | - Sara M. Rolfe
- Center for Developmental Biology and Regenerative Medicine, Seattle Children's Research Institute, Seattle, WA 98101, USA
- Friday Harbor Marine Laboratories, University of Washington, San Juan, WA 98250, USA
| | - Ronald Y. Kwon
- Department of Orthopedics and Sports Medicine, University of Washington, Seattle, WA 98195, USA
- Institute for Stem Cell and Regenerative Medicine, University of Washington, Seattle, WA 98109, USA
| | - A. Murat Maga
- Center for Developmental Biology and Regenerative Medicine, Seattle Children's Research Institute, Seattle, WA 98101, USA
- Division of Craniofacial Medicine, Department of Pediatrics, University of Washington, Seattle, WA 98105, USA
| |
Collapse
|
13
|
Savriama Y, Tautz D. Testing the accuracy of 3D automatic landmarking via genome-wide association studies. G3 (BETHESDA, MD.) 2022; 12:jkab443. [PMID: 35100368 PMCID: PMC9210295 DOI: 10.1093/g3journal/jkab443] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 08/27/2021] [Accepted: 12/16/2021] [Indexed: 11/13/2022]
Abstract
Various advances in 3D automatic phenotyping and landmark-based geometric morphometric methods have been made. While it is generally accepted that automatic landmarking compromises the capture of the biological variation, no studies have directly tested the actual impact of such landmarking approaches in analyses requiring a large number of specimens and for which the precision of phenotyping is crucial to extract an actual biological signal adequately. Here, we use a recently developed 3D atlas-based automatic landmarking method to test its accuracy in detecting QTLs associated with craniofacial development of the house mouse skull and lower jaws for a large number of specimens (circa 700) that were previously phenotyped via a semiautomatic landmarking method complemented with manual adjustment. We compare both landmarking methods with univariate and multivariate mapping of the skull and the lower jaws. We find that most significant SNPs and QTLs are not recovered based on the data derived from the automatic landmarking method. Our results thus confirm the notion that information is lost in the automated landmarking procedure although somewhat dependent on the analyzed structure. The automatic method seems to capture certain types of structures slightly better, such as lower jaws whose shape is almost entirely summarized by its outline and could be assimilated as a 2D flat object. By contrast, the more apparent 3D features exhibited by a structure such as the skull are not adequately captured by the automatic method. We conclude that using 3D atlas-based automatic landmarking methods requires careful consideration of the experimental question.
Collapse
Affiliation(s)
- Yoland Savriama
- Department Evolutionary Genetics, Max-Planck Institute for Evolutionary Biology, 24306 Plön, Germany
| | - Diethard Tautz
- Department Evolutionary Genetics, Max-Planck Institute for Evolutionary Biology, 24306 Plön, Germany
| |
Collapse
|
14
|
Porto A, Rolfe S, Maga AM. ALPACA: A fast and accurate computer vision approach for automated landmarking of three-dimensional biological structures. Methods Ecol Evol 2021; 12:2129-2144. [PMID: 35874971 PMCID: PMC9291522 DOI: 10.1111/2041-210x.13689] [Citation(s) in RCA: 38] [Impact Index Per Article: 9.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 07/11/2021] [Indexed: 11/27/2022]
Abstract
Landmark-based geometric morphometrics has emerged as an essential discipline for the quantitative analysis of size and shape in ecology and evolution. With the ever-increasing density of digitized landmarks, the possible development of a fully automated method of landmark placement has attracted considerable attention. Despite the recent progress in image registration techniques, which could provide a pathway to automation, three-dimensional (3D) morphometric data are still mainly gathered by trained experts. For the most part, the large infrastructure requirements necessary to perform image-based registration, together with its system specificity and its overall speed, have prevented its wide dissemination.Here, we propose and implement a general and lightweight point cloud-based approach to automatically collect high-dimensional landmark data in 3D surfaces (Automated Landmarking through Point cloud Alignment and Correspondence Analysis). Our framework possesses several advantages compared with image-based approaches. First, it presents comparable landmarking accuracy, despite relying on a single, random reference specimen and much sparser sampling of the structure's surface. Second, it can be efficiently run on consumer-grade personal computers. Finally, it is general and can be applied at the intraspecific level to any biological structure of interest, regardless of whether anatomical atlases are available.Our validation procedures indicate that the method can recover intraspecific patterns of morphological variation that are largely comparable to those obtained by manual digitization, indicating that the use of an automated landmarking approach should not result in different conclusions regarding the nature of multivariate patterns of morphological variation.The proposed point cloud-based approach has the potential to increase the scale and reproducibility of morphometrics research. To allow ALPACA to be used out-of-the-box by users with no prior programming experience, we implemented it as a SlicerMorph module. SlicerMorph is an extension that enables geometric morphometrics data collection and 3D specimen analysis within the open-source 3D Slicer biomedical visualization ecosystem. We expect that convenient access to this platform will make ALPACA broadly applicable within ecology and evolution.
Collapse
Affiliation(s)
- Arthur Porto
- Department of Biological SciencesLouisiana State UniversityBaton RougeLAUSA
- Center for Computation and TechnologyLouisiana State UniversityBaton RougeLAUSA
| | - Sara Rolfe
- Friday Harbor LaboratoriesUniversity of WashingtonSan Juan IslandWAUSA
- Center for Development Biology and Regenerative MedicineSeattle Children's Research InstituteSeattleWAUSA
| | - A. Murat Maga
- Center for Development Biology and Regenerative MedicineSeattle Children's Research InstituteSeattleWAUSA
- Division of Craniofacial MedicineDepartment of PediatricsUniversity of WashingtonSeattleWAUSA
| |
Collapse
|
15
|
Rolfe S, Pieper S, Porto A, Diamond K, Winchester J, Shan S, Kirveslahti H, Boyer D, Summers A, Maga AM. SlicerMorph: An open and extensible platform to retrieve, visualize and analyze 3D morphology. Methods Ecol Evol 2021; 12:1816-1825. [PMID: 40401087 PMCID: PMC12094517 DOI: 10.1111/2041-210x.13669] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/15/2021] [Accepted: 06/22/2021] [Indexed: 11/30/2022]
Abstract
1. Large scale digitization projects such as #ScanAllFishes and oVert are generating high-resolution microCT scans of vertebrates by the thousands. Data from these projects are shared with the community using aggregate 3D specimen repositories like MorphoSource through various open licenses. We anticipate an explosion of quantitative research in organismal biology with the convergence of available data and the methodologies to analyze them. 2. Though the data are available, the road from a series of images to analysis is fraught with challenges for most biologists. It involves tedious tasks of data format conversions, preserving spatial scale of the data accurately, 3D visualization and segmentations, acquiring measurements and annotations. When scientists use commercial software with proprietary formats, a roadblock for data exchange, collaboration, and reproducibility is erected that hurts the efforts of the scientific community to broaden participation in research. 3. We developed SlicerMorph as an extension of 3D Slicer, a biomedical visualization and analysis ecosystem with extensive visualization and segmentation capabilities built on proven python-scriptable open-source libraries such as Visualization Toolkit and Insight Toolkit. In addition to the core functionalities of Slicer, SlicerMorph provides users with modules to conveniently retrieve open-access 3D models or import users own 3D volumes, to annotate 3D curve and patch-based landmarks, generate landmark templates, conduct geometric morphometric analyses of 3D organismal form using both landmark-driven and landmark-free approaches, and create 3D animations from their results. We highlight how these individual modules can be tied together to establish complete workflow(s) from image sequence to morphospace. Our software development efforts were supplemented with short courses and workshops that cover the fundamentals of 3D imaging and morphometric analyses as it applies to study of organismal form and shape in evolutionary biology. 4. Our goal is to establish a community of organismal biologists centered around Slicer and SlicerMorph to facilitate easy exchange of data and results and collaborations using 3D specimens. Our proposition to our colleagues is that using a common open platform supported by a large user and developer community ensures the longevity and sustainability of the tools beyond the initial development effort.
Collapse
Affiliation(s)
- Sara Rolfe
- University of Washington, Friday Harbor Marine Laboratories, San Juan, WA
- Seattle Children's Research Institute, Center for Developmental Biology and Regenerative Medicine, Seattle, WA
| | | | - Arthur Porto
- Department of Biological Sciences, Louisiana State University, Baton Rouge
- Center for Computation and Technology, Louisiana State University, Baton Rouge
| | - Kelly Diamond
- Seattle Children's Research Institute, Center for Developmental Biology and Regenerative Medicine, Seattle, WA
| | - Julie Winchester
- Duke University, Department of Evolutionary Anthropology, Durham, NC
| | - Shan Shan
- Mount Holyoke College, Department of Mathematics, South Hadley, MA
| | | | - Doug Boyer
- Department of Biological Sciences, Louisiana State University, Baton Rouge
| | - Adam Summers
- University of Washington, Friday Harbor Marine Laboratories, San Juan, WA
| | - A Murat Maga
- Seattle Children's Research Institute, Center for Developmental Biology and Regenerative Medicine, Seattle, WA
- University of Washington, Department of Pediatrics, Division of Craniofacial Medicine, Seattle WA
| |
Collapse
|
16
|
Rolfe S, Davis C, Maga AM. Comparing semi-landmarking approaches for analyzing three-dimensional cranial morphology. AMERICAN JOURNAL OF PHYSICAL ANTHROPOLOGY 2021; 175:227-237. [PMID: 33483951 DOI: 10.1002/ajpa.24214] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/10/2020] [Revised: 11/13/2020] [Accepted: 12/08/2020] [Indexed: 01/29/2023]
Abstract
OBJECTIVES Increased use of three-dimensional (3D) imaging data has led to a need for methods capable of capturing rich shape descriptions. Semi-landmarks have been demonstrated to increase shape information but placement in 3D can be time consuming, computationally expensive, or may introduce artifacts. This study implements and compares three strategies to more densely sample a 3D image surface. MATERIALS AND METHODS Three dense sampling strategies: patch, patch-thin-plate spline (TPS), and pseudo-landmark sampling, are implemented to analyze skulls from three species of great apes. To evaluate the shape information added by each strategy, the semi or pseudo-landmarks are used to estimate a transform between an individual and the population average template. The average mean root squared error between the transformed mesh and the template is used to quantify the success of the transform. RESULTS The landmark sets generated by each method result in estimates of the template that on average were comparable or exceeded the accuracy of using manual landmarks alone. The patch method demonstrates the most sensitivity to noise and missing data, resulting in outliers with large deviations in the mean shape estimates. Patch-TPS and pseudo-landmarking provide more robust performance in the presence of noise and variability in the dataset. CONCLUSIONS Each landmarking strategy was capable of producing shape estimations of the population average templates that were generally comparable to manual landmarks alone while greatly increasing the density of the shape information. This study highlights the potential trade-offs between correspondence of the semi-landmark points, consistent point spacing, sample coverage, repeatability, and computational time.
Collapse
Affiliation(s)
- Sara Rolfe
- Friday Harbor Labs, University of Washington, Friday Harbor, Washington, USA
| | | | - A Murat Maga
- Department of Pediatrics, University of Washington, Seattle, Washington, USA
- Center for Developmental Biology and Regenerative Medicine, Seattle Children's Research Institute, Seattle, Washington, USA
| |
Collapse
|
17
|
A Registration and Deep Learning Approach to Automated Landmark Detection for Geometric Morphometrics. Evol Biol 2020; 47:246-259. [PMID: 33583965 DOI: 10.1007/s11692-020-09508-8] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Abstract
Geometric morphometrics is the statistical analysis of landmark-based shape variation and its covariation with other variables. Over the past two decades, the gold standard of landmark data acquisition has been manual detection by a single observer. This approach has proven accurate and reliable in small-scale investigations. However, big data initiatives are increasingly common in biology and morphometrics. This requires fast, automated, and standardized data collection. We combine techniques from image registration, geometric morphometrics, and deep learning to automate and optimize anatomical landmark detection. We test our method on high-resolution, micro-computed tomography images of adult mouse skulls. To ensure generalizability, we use a morphologically diverse sample and implement fundamentally different deformable registration algorithms. Compared to landmarks derived from conventional image registration workflows, our optimized landmark data show up to a 39.1% reduction in average coordinate error and a 36.7% reduction in total distribution error. In addition, our landmark optimization produces estimates of the sample mean shape and variance-covariance structure that are statistically indistinguishable from expert manual estimates. For biological imaging datasets and morphometric research questions, our approach can eliminate the time and subjectivity of manual landmark detection whilst retaining the biological integrity of these expert annotations.
Collapse
|
18
|
The Occasional Perils of Reflection (Across the Midline; in Geometric Morphometrics). Evol Biol 2020. [DOI: 10.1007/s11692-020-09501-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
19
|
Khot R, McGettigan M, Patrie JT, Feuerlein S. Quantification of gas exchange-related upward motion of the liver during prolonged breathholding-potential reduction of motion artifacts in abdominal MRI. Br J Radiol 2019; 93:20190549. [PMID: 31778311 DOI: 10.1259/bjr.20190549] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023] Open
Abstract
OBJECTIVE To test the hypothesis that there is a measureable upward motion of the diaphragm during prolonged breath-holds that could have a detrimental effect on image quality in liver MRI and to identify factor that potentially influence the magnitude of this motion. METHODS 15 healthy volunteers underwent MRI examination using prolonged breath-holds in the maximum inspiratory position and a moderate inspiratory position. Coronal T1 weighted three-dimensional gradient echo sequences of the entire thorax were acquired every 6 s during breath-holding allowing the calculation of total lung volume and the measurement of the absolute position of the dome of the liver. The potential impact of subject's gender, body mass index, and total lung capacity on the change in lung volume/diaphragmatic motion was assessed using random coefficient regression. RESULTS All volunteers demonstrated a slow reduction of the total lung volume during prolonged breath-holding up to 123 ml. There was a measurable associated upward shift of the diaphragm, measuring up to 5.6 mm after 24 s. There was a positive correlation with female gender (p = 0.037) and total lung volume (p = 0.005) and a negative association with BMI (p = 0.012) for the maximum inspiratory position only. CONCLUSION There is a measureable reduction of lung volumes with consecutive upward shift of the diaphragm during prolonged breath-holding which likely contributes to motion artifacts in liver MRI. ADVANCES IN KNOWLEDGE There is a measureable gas exchange-related reduction of lung volumes with consecutive upward shift of the diaphragm during prolonged breath-holding which likely contributes to motion artifacts in liver MRI. Correcting for this predictable upward shift has potential to improve image quality.The magnitude of this effect does not seem to be related to gender, BMI or total lung capacity if a moderate inspiratory position is used.
Collapse
Affiliation(s)
- Rachita Khot
- Department of Radiology and Medical Imaging, University of Virginia, Charlottesville, VA, United States
| | - Melissa McGettigan
- H Lee Moffitt Cancer Center and Research Institute, Tampa, 33612, United States
| | - James T Patrie
- Department of Public Health Sciences, University of Virginia, Charlottesville, VA, United States
| | - Sebastian Feuerlein
- H Lee Moffitt Cancer Center and Research Institute, Tampa, 33612, United States
| |
Collapse
|
20
|
Percival CJ, Devine J, Darwin BC, Liu W, van Eede M, Henkelman RM, Hallgrimsson B. The effect of automated landmark identification on morphometric analyses. J Anat 2019; 234:917-935. [PMID: 30901082 DOI: 10.1111/joa.12973] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 02/07/2019] [Indexed: 01/20/2023] Open
Abstract
Morphometric analysis of anatomical landmarks allows researchers to identify specific morphological differences between natural populations or experimental groups, but manually identifying landmarks is time-consuming. We compare manually and automatically generated adult mouse skull landmarks and subsequent morphometric analyses to elucidate how switching from manual to automated landmarking will impact morphometric analysis results for large mouse (Mus musculus) samples (n = 1205) that represent a wide range of 'normal' phenotypic variation (62 genotypes). Other studies have suggested that the use of automated landmarking methods is feasible, but this study is the first to compare the utility of current automated approaches to manual landmarking for a large dataset that allows the quantification of intra- and inter-strain variation. With this unique sample, we investigated how switching to a non-linear image registration-based automated landmarking method impacts estimated differences in genotype mean shape and shape variance-covariance structure. In addition, we tested whether an initial registration of specimen images to genotype-specific averages improves automatic landmark identification accuracy. Our results indicated that automated landmark placement was significantly different than manual landmark placement but that estimated skull shape covariation was correlated across methods. The addition of a preliminary genotype-specific registration step as part of a two-level procedure did not substantially improve on the accuracy of one-level automatic landmark placement. The landmarks with the lowest automatic landmark accuracy are found in locations with poor image registration alignment. The most serious outliers within morphometric analysis of automated landmarks displayed instances of stochastic image registration error that are likely representative of errors common when applying image registration methods to micro-computed tomography datasets that were initially collected with manual landmarking in mind. Additional efforts during specimen preparation and image acquisition can help reduce the number of registration errors and improve registration results. A reduction in skull shape variance estimates were noted for automated landmarking methods compared with manual landmarking. This partially reflects an underestimation of more extreme genotype shapes and loss of biological signal, but largely represents the fact that automated methods do not suffer from intra-observer landmarking error. For appropriate samples and research questions, our image registration-based automated landmarking method can eliminate the time required for manual landmarking and have a similar power to identify shape differences between inbred mouse genotypes.
Collapse
Affiliation(s)
| | - Jay Devine
- Department of Cell Biology and Anatomy, University of Calgary, Calgary, AB, Canada
| | - Benjamin C Darwin
- Mouse Imaging Centre, The Hospital for Sick Children, Toronto, ON, Canada
| | - Wei Liu
- Department of Cell Biology and Anatomy, University of Calgary, Calgary, AB, Canada
| | - Matthijs van Eede
- Mouse Imaging Centre, The Hospital for Sick Children, Toronto, ON, Canada
| | - R Mark Henkelman
- Mouse Imaging Centre, The Hospital for Sick Children, Toronto, ON, Canada.,Department of Medical Biophysics, University of Toronto, Toronto, ON, Canada
| | - Benedikt Hallgrimsson
- Department of Cell Biology and Anatomy, University of Calgary, Calgary, AB, Canada.,Alberta Children's Hospital Research Institute for Child and Maternal Health, University of Calgary, Calgary, AB, Canada.,The McCaig Institute for Bone and Joint Health, University of Calgary, Calgary, AB, Canada
| |
Collapse
|
21
|
Roosenboom J, Lee MK, Hecht JT, Heike CL, Wehby GL, Christensen K, Feingold E, Marazita ML, Maga AM, Shaffer JR, Weinberg SM. Mapping genetic variants for cranial vault shape in humans. PLoS One 2018; 13:e0196148. [PMID: 29698431 PMCID: PMC5919379 DOI: 10.1371/journal.pone.0196148] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/16/2017] [Accepted: 04/07/2018] [Indexed: 01/17/2023] Open
Abstract
The shape of the cranial vault, a region comprising interlocking flat bones surrounding the cerebral cortex, varies considerably in humans. Strongly influenced by brain size and shape, cranial vault morphology has both clinical and evolutionary relevance. However, little is known about the genetic basis of normal vault shape in humans. We performed a genome-wide association study (GWAS) on three vault measures (maximum cranial width [MCW], maximum cranial length [MCL], and cephalic index [CI]) in a sample of 4419 healthy individuals of European ancestry. All measures were adjusted by sex, age, and body size, then tested for association with genetic variants spanning the genome. GWAS results for the two cohorts were combined via meta-analysis. Significant associations were observed at two loci: 15p11.2 (lead SNP rs2924767, p = 2.107 × 10−8) for MCW and 17q11.2 (lead SNP rs72841279, p = 5.29 × 10−9) for MCL. Additionally, 32 suggestive loci (p < 5x10-6) were observed. Several candidate genes were located in these loci, such as NLK, MEF2A, SOX9 and SOX11. Genome-wide linkage analysis of cranial vault shape in mice (N = 433) was performed to follow-up the associated candidate loci identified in the human GWAS. Two loci, 17q11.2 (c11.loc44 in mice) and 17q25.1 (c11.loc74 in mice), associated with cranial vault size in humans, were also linked with cranial vault size in mice (LOD scores: 3.37 and 3.79 respectively). These results provide further insight into genetic pathways and mechanisms underlying normal variation in human craniofacial morphology.
Collapse
Affiliation(s)
- Jasmien Roosenboom
- Center for Craniofacial and Dental Genetics, Department of Oral Biology, University of Pittsburgh, Pittsburgh, PA, United States of America
| | - Myoung Keun Lee
- Center for Craniofacial and Dental Genetics, Department of Oral Biology, University of Pittsburgh, Pittsburgh, PA, United States of America
| | - Jacqueline T. Hecht
- Department of Pediatrics, University of Texas McGovern Medical Center, Houston, TX, United States of America
| | - Carrie L. Heike
- Department of Pediatrics, Seattle Children’s Craniofacial Center, University of Washington, Seattle, WA, United States of America
| | - George L. Wehby
- Department of Health Management and Policy, University of Iowa, Iowa City, IA, United States of America
| | - Kaare Christensen
- Department of Epidemiology, Institute of Public Health, University of Southern Denmark, Odense, Denmark
| | - Eleanor Feingold
- Department of Human Genetics, University of Pittsburgh, Pittsburgh, PA, United States of America
| | - Mary L. Marazita
- Center for Craniofacial and Dental Genetics, Department of Oral Biology, University of Pittsburgh, Pittsburgh, PA, United States of America
- Department of Human Genetics, University of Pittsburgh, Pittsburgh, PA, United States of America
| | - A. Murat Maga
- Department of Pediatrics, Seattle Children’s Craniofacial Center, University of Washington, Seattle, WA, United States of America
- Center for Developmental Biology and Regenerative Medicine, Seattle Children’s Research Institute Seattle, WA, United States of America
| | - John R. Shaffer
- Center for Craniofacial and Dental Genetics, Department of Oral Biology, University of Pittsburgh, Pittsburgh, PA, United States of America
- Department of Human Genetics, University of Pittsburgh, Pittsburgh, PA, United States of America
| | - Seth M. Weinberg
- Center for Craniofacial and Dental Genetics, Department of Oral Biology, University of Pittsburgh, Pittsburgh, PA, United States of America
- Department of Human Genetics, University of Pittsburgh, Pittsburgh, PA, United States of America
- Department of Anthropology, University of Pittsburgh, Pittsburgh, PA, United States of America
- * E-mail:
| |
Collapse
|