Summary
מבנים ביולוגיים רבים חסרים ציוני דרך להגדרה בקלות, מה שהופך אותו קשה ליישם שיטות morphometric המודרנית. כאן אנו מדגימים שיטות ללמוד את baculum העכבר (עצם הפין), כולל לנתיחה וסריקת microCT, ואחריו שיטות חישוביות להגדיר-ציוני למחצה המשמשים לכמת גודל וצורת וריאציה.
Abstract
morphometrics המודרנית מספקת שיטות רבות עצמה כדי לכמת וריאצית גודל וצורה. דרישה בסיסית רשימת קואורדינטות מגדירות ציוניות; אולם קואורדינטות כזה חייב לייצג מבנים הומולוגיים פני דגימות. בעוד אובייקטים ביולוגיים רבים מורכבים של ציוני דרך זיהו בקלות כדי לספק את ההנחה של הומולוגיה, רבים חסרים מבנים כאלה. פתרון אפשרי אחד הוא מתמטי-ציונים חצי מקום על אובייקט המייצג באותו האזור מורפולוגיים ברחבי דגימות. כאן, אנו ממחישים צינור שפותח לאחרונה להגדיר למחצה ציוני דרך מתמטית מן baculum העכבר (עצם הפין). השיטות שלנו צריכות לחול על מגוון רחב של חפצים.
Introduction
תחום morphometrics כולל מגוון של שיטות לכמת את הגודל והצורה של הטופס ביולוגי, צעד מהותי בחקירה מדעית 1, 2, 3, 4, 5, 6. באופן מסורתי, הניתוח הסטטיסטי של גודל וצורה מתחיל על ידי זיהוי נקודות ציון על מבנה ביולוגי, ולאחר מכן מדידת מרחקים ליניארי, זוויות יחסיות, אשר יכול להיות מנותחת מסגרת רבה-משתנית. לנדמרק המבוסס הגיאומטרית Morphometrics הוא גישה אשר שומרת את המיקום המרחבי של ציוני דרך, שימור מידע גיאומטרי של איסוף נתונים באמצעות ניתוח להדמיה 5. ניתוח סדום כללי (GPA) יכול להיות מיושם כדי להסיר וריאצית מיקום, סולם, וסיבוב של ציונים דרך לייצר יישור בין דגימות מinimizes ההבדלים בריבוע שלהם - מה שנותר הוא שוני בצורה 7.
מושג חשוב של כל ניתוח morphometric הוא הומולוגיה, או את הרעיון כי אחד מהימן יכול לזהות את האתרים המייצגים תכונות משמעות ביולוגית ודיסקרטיות המתאימים בין דגימות או מבנים. לדוגמא, יש גולגלות אדם תהליכים הומולוגיים, foramina, תפרים, ותעלות שיכולים לאפשר morphometric מנתח. למרבה הצער, את זיהוי של ציוני דרך מקבילים קשה על פני רבי מבנים ביולוגיים, במיוחד אלה עם משטחים חלקים או עקומות 8, 9, 10.
אנחנו מתקרבים בעיה זו מתחת שימוש בגיאומטריה חישובית. זרימת העבודה הכללית היא ליצור סריקה תלת ממדים של האובייקט שניתן לייצג כענן של נקודות, ולאחר מכן לסובב ולהפוך כי ענן נקודה כך שכל זהpecimens מכוון על מערכת קואורדינטות משותפת. אז אנחנו מתמטיים להגדיר-ציונים חצי מאזורים ספציפיים של האובייקט. חצי ציוני דרך דיסקרטית דגש על אזורים כאלה הן ביולוגיים שרירותיות 11. ניצוח GPA ניתוחים סטטיסטיים רצף יכול לייצר חפצים לא רצויים 8, 12 משום ציונים דרך להציב באופן שרירותי לא יכול להיות ביולוגי הומולוגי. לכן, אנו נאפשר-ציונים חצי אלה מתמטיים "להחליק". הליך זה מצמצם את הפרש הפוטנציאלים בין מבנים. כפי שטען במקום אחר האלגוריתם זז משמש כאן ראוי לכמת אזורי אנטומיים דומים חסרי זיהה בקלות מקבילה ציוני דרך 3, 6, 8, 10, 11, 12. שיטות אלה יש li שלהםmitations 13, אבל צריך להיות בעל יכולת הסתגלות לאובייקטים של גודל וצורה שונים.
כאן, אנו מדגימים כיצד שיטה זו יושמה במחקר שנערך לאחרונה של baculum העכבר 14, עצם הפין כי כבר צבר ואבד מספר פעמים עצמאיות במהלך האבולוציה של יונקי 15. אנחנו דנים לנתיחה והכנת עצם ספציפי, baculum (פרוטוקול 1), הדור של תמונות microCT (פרוטוקול 2), ואת המרה של תמונות אלה לתבנית המאפשרת לכל בגיאומטריה חישובית במורד הזרם (פרוטוקולים 3 ו -4). אחרי השלבים הבאים, כל דגימה מיוצגת על ידי ~ 100K XYZ קואורדינטות. לאחר מכן, אנו עוברים דרך סדרה של טרנספורמציות כי ביעילות ליישר את כל הדגימות לתוך אורינטציה משותפת (הפרוטוקול 5), ואז להגדיר-ציוני חצי דגימות מיושרות (פרוטוקול 6). פרוטוקולים 1-4 צריכים להיות דומים ללא קשר של האובייקט להיות מנותח. 5 פרוטוקול ופרוטוקול 6 הם SPEcifically מיועד baculum, אך תקוותנו היא כי באמצעות פירוט הפעולות הבאים, חוקרים יכולים לדמיין שינויים זה יהיה רלוונטי עבור מושא העניין שלהם. לדוגמה, שינויים בשיטות אלה יושמו ללמוד עצמות לווייתנים האגן ועצמות הצלע 16.
Subscription Required. Please recommend JoVE to your librarian.
Protocol
כל הנהלים ואנשי אושרו על ידי אוניברסיטת מכון דרום קליפורניה טיפול בבעלי חיים ועדת שימוש (IACUC), פרוטוקול # 11,394.
1. Baculum Dissection והכנה
- להרדים עכבר זכר בוגר מינית באמצעות פחמן דו חמצני חשיפת יתר, על פי פרוטוקולים שנקבעו על ידי הוועדה לטיפול המוסדית בעלי החיים והשימוש הרלוונטי (IACUC).
- הנח את החיה במצב שכיבה, ו מתארך הפין על ידי הפעלת לחץ בעזרת אגודלים לרוחב פתח preputial.
- לאחר הפין הוא ממושך, להאריך את הרקמה דרך העורלה ככל האפשר.
- עם מספריים, לחתוך את הפרוקסימלי גוף הפין אל העטרה שבו baculum מתגורר.
- מעביר את הפין הגזור לצינור מיליליטר 1.7 ולהוסיף 200 μL מי ברז. ודא כי הפין הוא שקוע לחלוטין בתוך הנוזל.
- דגירת הרקמות במים ב ~ 50 מעלות צלזיוס במשך 3-5 ימים.
- מניח את baculum גזור לתוך צינור microcentrifuge חדש עם הכובע פתוח. השאר מכסה פתוח O / N להתייבש עצם.
2. סריקת microCT
- לחץ בעל סריקה גלילי microCT לתוך לבנים של קצף פרחים ליצור גליל קצף פרחים.
- חלץ את גליל קצף הפרחים לחתוך פרוסות ~ 2-5 ס"מ עובי.
- Push היבש bacula לתוך קצף הפרחים, סביב בפריפריה של פרוסה בודדת כדי למזער הפרעה במהלך סריקה. האורינטציה המדויקת של עצמות חובה לציין המאפשר זיהוי נכון של דגימות בודדות בפרוטוקול 4.
- בעדינות במקום פרוסה עם העצמות מוטבעות לתוך מחזיק microCT.
- רוכשת סריקות microCT. במקרה של העכבר bacula 14
3. עיבוד microCT: המרת .DCM סטאק לקובץ .xyz יחיד
הערה: כל סריקת microCT מייצרת ערימה של .DCM, או "DICOM", קבצים מייצגים פרוסות תמונה נלקח דרך האובייקט. כל בגיאומטריה חישובית במורד דורשת קבצי .xyz שטוחים, אשר היא פשוט קובץ טקסט המכיל ארבעה עמודים - x, y, ו Z קואורדינטות של כל פיקסל, ואת עוצמת פיקסל, החל -5000 (שחור) עד 5000 (לבן). סף פיקסל מעל 3,000 בדרך כלל עובד גם סף להגדרת עצמות.
- התקן פייתון (www.python.org) ואת המודולים "פייתון" פקודות, DICOM, PYLAB, SYS, ו numpy.
- 01_pr פתח "ocess_dicom.py "{Figshare} עם כל עורך טקסט. על פי סעיף משתנים, נתיב השינוי, ספי פיקסל, ושמות הספרייה לפי הצורך.
- הפעל "פיתון 01_process_dicom.py". התקדמות תודפס למסך. בתוך כל תיקייה בשם בשלב 3.2, קובץ חדש מובא בשם; למשל, directory_name.PT3000.xyz, שבו PT3000 מציין את סף פיקסל שצוין בשלב 3.2.
4. עיבוד microCT: פילוח-אאוט פרט דגימת .xyz קבצים
- התקן R (https://www.r-project.org/) עם RGL הספרייה.
- פתחו את הקובץ '02_segment_dicoms.r' {Figshare} עם כל עורך טקסט. על פי הסעיף משתנה, לשנות את שם הנתיב להצביע על קובץ .xyz נוצר בפרוטוקול 3 לעיל.
- מתוך R, הפעל את הפקודה "המקור (" 02_segment_dicoms.r ') "(ללא מרכאות).
- לאחר שלושת ממדי התמונה של הקובץ .xyz נוצר בפרוטוקול 3 מופיע, הזן את number של דגימות בקובץ .xyz הכולל. עליו מדבקה ובחר את נקודות מכל הדגימה באמצעות פונקציות גלילה וזום.
הערה: ברקע, קבצי .xyz נפרדים ייעשו עבור כל דגימה. אלו מופיעות ספרייה בשם, למשל, XYZ_FILES_PT3000, שבו PT3000 מציין את סף פיקסל בשימוש.
5. "יישור" הדגימה .xyz קבצים לקואורדינטות נפוצות.
- פתח את סקריפט פייתון "03_transform.py" {Figshare}, המחייב את mattdean_modules.py מודול נוסף {Figshare}, וכן שני יישומים עומד עצמית: "rotate_translate_cylindrical" (https://github.com/timydaley/dean_cylindrical_tranform) ו "qconvex" (www.qhull.org/html/qconvex.htm) המשמשים התסריט הזה.
- על פי הסעיף משתנה, לזהות את שמות נתיב המלאים mattdean_modules.py, rotate_path ו qconvex_dir. בנוסף, לזהות את הנתיב המלא אל הספרייה המכילה את .xy הפרטקבצי z שנוצרו בשלב 4.
- 03_transform.py לרוץ, אשר יוצר קובץ חדש לכל דגימה עם הסיומת .TRANSFORMED.xyz.
6. דגימה "חיתוך" מיושרת .xyz קבצים לזהות-ציונית חצי.
- להרחיב ולהפעיל את סקריפט פייתון "04_identify_landmarks.py" {Figshare}. במקטע משתנה, לזהות שמות מלאי נתיב לספרייה המכילה את קבצי .TRANSFORMED.xyz. סקריפט זה מזהה 802-ציוני למחצה, שניתן להשתמש בהם כדי לכמת הגודל והצורה של המבנה.
Subscription Required. Please recommend JoVE to your librarian.
Representative Results
קואורדינטות xyz של האתרים ההסטוריים חצי יוצר הפרוטוקול 6 ניתן לייבא ישירות לתוך כל ניתוח morphometrics הגיאומטרי מבוסס ציון 17. הצינור חישובית מעל יושם ללמוד עכבר bacula 14, כמו גם אגן לווייתן עצמות צלעות 16. פרטים נוספים על ההגדרה חישובית של אתרים היסטוריים חצי מוצגים כאן, בניסיון לסייע לחוקרים לחזות צעדים שהיו יכולים להיות שונה כדי להכיל אובייקט של מושא העניין שלהם. Baculum מכיל כמה תכונות ייחודיות כי נוצלו כדי להפוך טרנספורמציות מסוימות. לדוגמה, לאחר חיתוך העצם המחשוב לשני חצאים לאורך ציר הפרוקסימלי-דיסטלי, זיהינו במחצית הפרוקסימלי פשוט על ידי השוואת המספר הכולל של נקודות (הפרוקסימלי יש יותר). כל עוד קיימים מאפיינים ייחודיים כגון זה, השיטות שלנו צריכות להיות בעל יכולת הסתגלות לכללְהִתְנַגֵד. בנוסף, יש להדגיש כי אנו אמפיריים נקבענו לסף מסוים, כגון "10% הפרוקסימלי" שביצועים שלהן במחקרי baculum שלנו, אבל בהחלט צריך להיות משוערך עבור אובייקטים אחרים.
החל פרוטוקול 5, צעד חישובית הראשון הוא לחשב את קמור (את האוסף המזערי של נקודות המכילות את כל נקודות אחרות במדגם) כדי לזהות את שתי נקודות מרוחקות קליעה אחד מהשני. שתי נקודות אלה (בתחומים אדומים, תרשים 1C) מתחילות להגדיר Z- ציר חדש (קו אדום, תרשים 1C) שיוצא הפרוקסימלי-דיסטלי דרך העצם. במקרה של baculum, במחצית של ענן הצבע שמכיל יותר נקודות מוגדרות סוף הפרוקסימלי.
שנית, ענן הצבע כולו עובר טרנספורמציה כך שהנקודה הפרוקסימלי לוקח על x, y, z קואורדינטות של 0,0,0 ואת t נקודת דיסטליakes על x, y, z קואורדינטות של 0,0, + z, שם + z הוא קצת ערך חיובי תלוי בגודל של העצם. בסוף שלב זה, Z- ציר עובר דרך באורך של העצם. לקבלת הנחיות מפורטות להלן, האורך מן המינימום z-קואורדינטות מרביות יהיה המכונה Zlength.
שלישית, כדי לתקן את השונות הקשורות מיקום המדויק של הנקודות הפרוקסימלית ומ דיסטלי שזוהו לעיל, 10% הפרוקסימלי ביותר ו -10% נקודות דיסטלי ביותר נדגמים בנפרד (איור 1D), זיהה centroids שלהם (ספירות אדומות, איור 1E) וענן נקודת טרנספורמציה כך centroid הפרוקסימלי הוא 0,0,0, ואת centroid דיסטלי הוא 0,0, + z, עם Z- ציר חדש שעובר דרך מרכז הדגימה (קו אדום, איור 1E) .
רביעית, ענן הצבע הוא הסתובב סביב ציר z ידי לקיחת הראשונה slicדואר של נקודות הפרוקסימלי 15 כדי 15.25% Zlength של המבנה (נקודות כחולות, איור 1E). פרוסת נקודות זה הוא שטוח בממד z (כלומר z-הקואורדינטות הם התעלמו פשוט), לקח את קמור, ואת המלבן התוחם מינימום (המלבן הקטן ביותר המכיל את כל נקודות אחרות) מחושב. דמיין קו מחבר midpoints של שני הצדדים הקצרים של מלבן תוחם מינימום זה. אנחנו לסובב את ענן הצבע עד שתיים האמצע אלה הופכים -x, 0, + z ו + x, 0, + z, בהתאמה, וכך הקו הזה הופך בציר Y חדש. לאחר השינוי, את המרחק בין ערכי x מקסימלי ומינימלי מכונים Xlength. קובץ חדש נוצר מן specimen.xyz כדי specimen.TRANSFORMED.xyz.
החמישי, מצביע בתוך 1% Xlength של ציר z הם חתוך החוצה (נקודות כחולות, האיור 1F), ונקודת דיסטלי האחת הכי קרובה והרווקים המזוהים ביותר מ SL המרכזי זהדיסטלי ו הפרוקסימלי קרח ומסומן, בהתאמה. אלה הם שני הציונים דרך למחצה הראשון מזוהים.
שישית, 50 פרוסות במרווחים שווים זה מזה נקודות נדגמים לאורך ציר z (נקודות אדומות, איור 1G). כל פרוסה היא בעובי של 1% Zlength. פרוסה כל בברוטליות ואז במימד z, ומחולק באופן שווה על ידי 7 קווים אנכיים (קווים אדומים, איור 1H). נקודות בתוך 2% Xlength של כל שורה נשמרות (נקודות אדומות, 1H האיור), אז נקודות עם המינימום והמקסימום y לתאם נשמרים, מוקרן כל שורת בהתאמה, ומסומן גחון ועל גב, בהתאמה. בנוסף, תוויות להכיל את המספר הפרוס ואת מספר השורה, למשל P15_VENTRAL4 הוא נקודת הגחון שנדגמו הקו האנכי 4 th של פרוסת ה -15. חשוב לציין, כל נקודה שכותרתו, למשל, P15_VENTRAL4, מתרחש פעם אחת ורק פעם אחת בכל הדגימות, שימור corres pondence. בנוסף לנקודות הגחון ועל הגב של כל אחד 7 הקווים (14 למחצה ציוני סך הכל), נקודות עם מקסימלית x-הערך המינימלי נדגמים ומסומן על ימין ועל שמאל, בהתאמה. קואורדינטות Y ו- Z של שמאל וימין הם החליקו באמצעות פונקצית lowess ב ר עבור baculum, סך של 16 למחצה ציוני דרך מוגדרים לכל פרוסה (ספירות אדומות, 1H האיור); עם 50 פרוסות בתוספת הפרוקסימלי הציוני-למחצה דיסטלי שהוגדר לעיל, 802-ציוני חצי נדגמים לכל דגימה (ספירות ירוקות, איור 1 ט). כל שאר הנקודות מן סריקת microCT המקורית מבוטלות.
יצוין כי למרות קוטביויות גחון / הגבי ו הפרוקסימלי / דיסטלי נקבעו באופן מתמטי, כל מערכי דגימת אושרו חזותיים מותאמים באופן ידני כנדרש. במדגם שלנו של 369 bacula, כ 10 צריך להיות מותאם באופן ידני.
er.within-page = "1"> 
איור 1: ייצוג חזותי של זרימת עבודת חישובית (פרוטוקול 4-6). (א) צילום מסך מתסריט 02_segment_dicoms.r (פרוטוקול 4), מראה את המשימה של ענני צבע ייחודיים פריטים שונים. (ב) תצוגה מוגדלת של baculum אחד, מיוצג על ידי ענן של נקודות xyz ~ 100K. (ג) זיהוי של שתי הנקודות רחוקות זו מזו (ספירות אדומות), המשמש להגדרת Z- ציר חדש שיוצא הפרוקסימלי-דיסטלי דרך העצם (הקו האדום). (ד) דגימה הפרוקסימלי-רוב 10% ו דיסטלי ביותר 10% של נקודות (נקודות אדומות) מספקים אמצעי כדי להתאים את השונות קלות המיקום המדויק של ציר z. (ה) centroids של הפרוקסימלי-רוב 10% ו דיסטלי ביותר 10% (בתחומים האדומים) משמשת להגדיר Z- ציר חדש (קו אדום). לאחר מכן, פרוס נקודות נופלים בין 15.00-15.25% של Z- ציר חדש זה (נקודות כחולות) נלקחו לחשב מלבן תוחם המינימום. ענן הנקודה מסובב עד הצד הארוך של מלבן תוחם המינימום מקביל בציר Y חדש. F) פרוס נקודות רצו לאורך קו האמצע (נקודות הכחולות) נדגמים ואת הפרוקסימלי-ביותר דיסטלי ביותר הנקודה מוגדרת כ-ציון למחצה. G) 50 פרוסות במרווחים שווים זה מזה נקודות נדגמים (נקודות אדומות), באות H) מראות פרוסה אחת כזו. 16 נקודות (עיגולים אדומים) מוגדרות כדי ללכוד את קווי המתאר של כל פרוסה. I) כאשר חזר בכל הפרוסות, סך של 802 כדורים למחצה ציוני דרך (ירוקים) להגדיר את המבנה והם משמשים בכל יישומי morphometric במורד הזרם. אנא לחץ כאן כדי לצפות בגרסה גדולה יותר של דמות זו.
Subscription Required. Please recommend JoVE to your librarian.
Discussion
השלבים הקריטיים בפרוטוקול לעיל הם 1) לנתח את bacula, 2) איסוף תמונות microCT, 3) המרת פלט microCT לקובץ שטוח של XYZ קואורדינטות, 4) פילוח מתוך כל ענן הנקודה של דגימה, 5) הפיכה כל דגימה סטנדרטי מערכת קואורדינטות, ו -6) הגדרה-ציוני למחצה. כל הצעדים הללו היא לשנות בקלות כדי להתאים חפצים שונים.
שיטות אלה ניתן ליישם עשויים חפץ שהוא בעצם "מוט בצורה", או לפחות לא מעוקל מדי. אובייקטים כי עקום בחזרה על עצמם להיות "בצורת U" לא ניתן לנתח כרגע, בגלל חיתוך (1G איור) יחזור נקודות מחלקים שונים של האובייקט. חפצים כאלה יכולים לאכלס בעתיד על ידי מחשוב מיישר את האובייקט לפני החיתוך.
הצגנו שיטה כללית-ציוני חצי מגדיר באופן מתמטי מן הצורות חסרות sציוני דרך Olid. שיטות כלליות אלה שונו על מנת ללמוד את ההתפתחות של עצמות האגן צלע לוויתן 16, אשר יש צורות שונות מאוד. השיטות חישוביות שלנו ציוניות דרך הגדרה צריכות לחול על כל סדרה של קואורדינטות XYZ. העסקנו סריקת microCT כאן, בהתחשב בגודל הקטן של עכבר bacula 14. עבור עצמות גדולות, כגון עצמות האגן צלע לווייתן, העסקנו סורק לייזר שיחזרו את פני השטח של העצמות 16. חשוב חזותי לבדוק כל הסטים של ציונים דרך חצי לאימות האיכות של השיטה. היתרון העיקרי של השיטות חישוביות שלנו הוא שהם דווקא לכמת את הגודל וצורת וריאציה, ולשמר את התכתובת בין האזורים השונים של האובייקט.
Subscription Required. Please recommend JoVE to your librarian.
Acknowledgments
טים דיילי ואנדרו סמית נתן הרבה דיונים חישובית שימושי במהלך הימים הראשונים; טים דיילי כתב התוכנית rotate_translate_cylindrical הכרחי פרוטוקול 5. משאבים חישובית נמסרו על ידי אשכול מחשוב עתיר ביצועים באוניברסיטת דרום קליפורניה. עבודה זו נתמכה על ידי GM098536 # מענק NIH (MDD).
Materials
| Name | Company | Catalog Number | Comments |
| Dissecting scissors | VWR | 470106-338 | Most sizes should work |
| Dissecting Forceps, Fine Tip, Curved | VWR | 82027-406 | |
| 1.7 mL microcentrifuge tube | VWR | 87003-294 | |
| Absolute Ethanol | Fisher Scientific | CAS 64-17-5 | To be diluted to 70% for dissections |
| Floral Foam | Wholesale Floral | 6002-48-07 | |
| uCT50 scanner | Scanco Medical AG, Bruttisellen, Switzerland |
References
- Slice, D. E. Geometrics morphometrics. Annu. Rev. Anthropol. 36, 261-281 (2007).
- Slice, D. E. Modern morphometrics in physical anthropology. 6, Springer. (2005).
- Zelditch, M. L., Swiderski, D. L., Sheets, H. D. Geometric morphometrics for biologists: a primer. , 2nd, Elsevier. (2012).
- Bookstein, F. Morphometric tools for landmark data: geometry and biology. , Cambridge University Press. (1991).
- Rohlf, F. J., Marcus, L. F. A Revolution in Morphometrics. Trends. Ecol. Evol. 8 (4), 129-132 (1993).
- Zelditch, M. L., Swiderski, D. L., Sheets, H. D., Fink, W. L. Geometric morphometrics for biologists: a primer. , Elsevier. (2004).
- Rohlf, F. J., Slice, D. E. Extensions of the Procrustes method for the optimal superimposition of landmarks. Syst. Zool. 39 (1), 40-59 (1990).
- Gunz, P., Mitteroecker, P. Semilandmarks: a method for quantifying curves and surfaces. Hystrix. 24 (1), 103-109 (2013).
- Gunz, P., Ramsier, M., Kuhrig, M., Hublin, J. J., Spoor, F. The mammalian bony labyrinth reconsidered, introducing a comprehensive geometric morphometric approach. J. Anat. 220 (6), 529-543 (2012).
- Mitteroecker, P., Gunz, P. Advances in geometric morphometrics. Evol. Biol. 36 (2), 235-247 (2009).
- Bookstein, F. J. Landmark methods for forms without landmarks: morphometrics of group differences in outline shape. Med. Im. Anal. 1 (3), 225-243 (1997).
- Gunz, P., Mitteroecker, P., Bookstein, F. Modern morphometrics in physical anthropology. Slice, D. E. , Kluwer Press. 73-98 (2005).
- Oxnard, C., O'Higgins, P. Biology Clearly Needs Morphometrics. Does Morphometrics Need Biology? Biological Theory. 4 (1), 84-97 (2009).
- Schultz, N. G., et al. The genetic basis of baculum size and shape variation in mice. G3. 6 (5), 1141-1151 (2016).
- Schultz, N. G., Lough-Stevens, M., Abreu, E., Orr, T. J., Dean, M. D. The baculum was gained and lost multiple times during mammalian evolution. Integr Comp Biol. 56 (4), 644-656 (2016).
- Dines, J. P., et al. Sexual selection targets cetacean pelvic bones. Evolution. 68 (11), 3296-3306 (2014).
- Adams, D. C., Otárola-Castillo, E. geomorph: an R package for the collection and analysis of geometric morphometric shape data. Methods Ecol. Evol. 4 (4), 393-399 (2013).
