在中短肛解剖，显微扫描和分析形态

摘要

许多生物结构缺乏容易定义的地标，从而难以应用现代形态的方法。在这里，我们举例说明的方法来研究的小鼠阴茎骨（在阴茎骨），包括夹层和显微扫描，随后计算方法来定义用于量化的尺寸和形状的变化的半标。

摘要

现代形态测量提供了强大的方法来量化的大小和形状的变化。一个基本要求是，确定地标坐标列表;然而，这样的坐标必须代表整个标本同源结构。虽然许多生物对象包括容易识别标志，以满足同源性的假设，许多缺乏这种结构。一个可能的解决方案是在数学上横跨试样代表相同形态区域中的对象在地方半标。在这里，我们说明了一个新近开发的管道数学定义从鼠标阴茎骨（阴茎骨）半地标。我们的方法应该适用于广泛范围的对象。

引言

形态测量的领域包括的方法来量化的大小和生物形式的形状，在科学调查^{1，2，3，4，5，6}的一个基本步骤的分集。传统上，大小和形状的统计分析开始于识别上的生物结构的地标，然后测量线性距离，角度和比例，这可以在一个多变量的框架进行分析。里程碑式的基于几何形态特征是保留的地标空间位置，通过分析和可视化⁵保存从数据收集几何信息的方法。广义普鲁克分析（GPA）可以应用到消除位置，规模，和地标旋转变化，以产生是m标本之间的对准inimizes他们的平方差-剩下的就是形状相异^7。

任何形态分析的一个重要概念是同源性，或一个可以可靠地识别代表该标本或结构之间对应生物学意义的和离散的特征的地标的想法。例如，人头骨有同源流程，椎间孔，缝合线，和管道，可以使形态学分析。不幸的是，对应的地标标识是在许多生物结构困难，特别是那些具有光滑表面或曲线^{8，9，10。}

我们下面的办法使用计算几何这个问题。的一般工作流程，以产生可以被表示为点的云的物体的三维扫描，然后旋转和转换，该转换点云，使得所有specimens被定向一个共同的坐标系上。然后我们数学定义从物体的特定区域的半标。放在这些地区离散半地标生物任意^11。开展GPA和随后的统计分析可以产生不良的文物^8，12，因为随意放置地标可能不是生物同源。因此，我们允许这些半地标数学"滑动"。此过程减少结构之间的电势差。如别处认为这里使用的滑动算法是适当的量化缺乏类似解剖区域容易识别相应的地标^{3，6，8，10，11，12。}这些方法各有其利mitations ^13，但应该适应不同尺寸和形状的物体。

这里，我们说明这个方法是如何在近一个鼠标阴茎骨¹⁴的研究应用，在已经获得和哺乳动物进化¹⁵期间丢失多个独立倍阴茎骨。我们讨论解剖和制备的特定骨的，所述阴茎骨（协议1），生成显微图像（协议2），并且这些图像以一种格式，使所有下游计算几何变换（协议3和4）。完成这些步骤后，由〜100K XYZ坐标每个试样表示。然后，我们通过一系列的变革，有效地调整所有标本到一个共同的方向（协议5）走，然后定义从对准标本半地标（协议6）。协议1-4应该是相似的，无论被分析的对象。协议5和6号议定书是SPEcifically专为阴茎骨，但我们希望通过详细介绍这些步骤，研究者可以想像的修改，这将是对有关他们感兴趣的对象。例如，这些方法改进了用于研究鲸骨盆和肋骨^16。

研究方案

所有的过程和人员由南加州学院的大学动物护理和使用委员会（IACUC），协议＃11394批准。

1.解剖短肛和准备

1. 通过二氧化碳过度曝光安乐死一个性成熟的雄性小鼠，根据相关机构的动物护理和使用委员会（IACUC）提出的协议。
2. 躺在仰卧位的动物，并通过施加压力拖延阴茎大拇指外侧包皮开幕。
3. 一旦阴茎旷日持久，通过包皮尽可能延长组织。
4. 用剪刀，剪断阴茎体近端，其中阴茎骨所在的龟头。
5. 转移解剖阴茎1.7毫升管中，加入200μl自来水。使某些阴茎在液体完全浸没。
6. 孵育在水中的组织在约50℃下3-5天。
7. 适当温育后，从阴茎骨除去周围组织，用解剖显微镜下用镊子。轻轻喷70％乙醇推下剩余的组织和清理骨头。
8. 放置解剖阴茎骨与盖打开一个新的离心管中。请假封顶开O / N为干骨。

2.显微扫描

1. 按一个显微扫描圆柱支架插入花店泡沫砖创造花店泡沫的圆柱体。
2. 提取花店泡沫缸和切割片〜2-5厘米厚。
3. 干燥后推入的bacula花店泡沫，围绕单个片的外围扫描过程中，尽量减少干扰。骨骼的精确定位应该说允许在协议4个人标本的正确识别。
4. 轻轻将切片与嵌入式骨头进入显微持有人。
5. 获得显微扫描。在小鼠的bacula ¹⁴的情况下，我们根据以下设置在南加州大学分子影像中心使用的uCT50扫描器（SCANCO医疗AG，Bruttisellen，瑞士） 500毫秒，为15.5mm的体素尺寸。

3.显微加工:转换一个.DCM堆栈到单个名为.xyz文件

注意:每个显微扫描产生一个堆栈.DCM的，或"DICOM"，即表示通过对象拍摄图像片段文件。所有下游计算几何需要平面名为.xyz文件，这是一个简单的包含四列的文本文件 - 在x，y和每个像素的z坐标的，并且所述像素的强度，范围从-5000（黑色）到5000 （白色）。超过3000像素的阈值通常工作以及定义骨骼的阈值。

1. 安装Python（www.python.org）和Python模块命令，DICOM，PYLAB，SYS和numpy的。
2. 打开"01_process_dicom.py"{} Figshare用任何文本编辑器，在变量部分，改变路径，像素阈值和目录名是必要的。
3. 运行"蟒蛇01_process_dicom.py"。进展情况将被打印到屏幕上。在每个步骤3.2命名的目录，新文件由命名;例如，directory_name.PT3000.xyz，其中PT3000指示像素阈在步骤3.2表示。

4.显微加工:细分出单独的标本名为.xyz文件

1. 与库RGL安装R（https://www.r-project.org/）。
2. 用任何文本编辑器打开文件02_segment_dicoms.r'{} Figshare。根据变量部分，更改路径名称指向上述协议3中创建的名为.xyz文件。
3. 从内部R，运行命令"源（'02_segment_dicoms.r'）"（不包括双引号）。
4. 在协议3中创建的名为.xyz文件的三维图像出现后，输入numbe在整体名为.xyz文件标本河然后贴上标签并选择使用滚动和缩放功能，每个样本点。
   注:在后台，独立名为.xyz文件将每个样品上进行。这些出现在一个新的目录，例如，XYZ_FILES_PT3000，其中PT3000表示使用的像素阈值。

5."对齐"标本名为.xyz文件到通用坐标。

1. 打开Python脚本"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）受此脚本中使用。
2. 根据变量部分，确定全路径名mattdean_modules.py，rotate_path和qconvex_dir。此外，确定到包含个体.xy目录的完整路径步骤4中创建ž文件。
3. 运行03_transform.py，它创建每个试样的新文件的后缀.TRANSFORMED.xyz。

6."切片"不结盟标本名为.xyz文件识别半地标。

1. 打开并运行Python脚本"04_identify_landmarks.py"{} Figshare。在变量部分，确定全路径名包含.TRANSFORMED.xyz文件的目录。此脚本标识可以用来量化该结构的尺寸和形状802半标。

结果

在协议6中产生的半地标XYZ坐标可以直接导入到任何基于具有里程碑意义的几何形态测量分析^17。计算管道上面已经用于研究小鼠的bacula ^14，以及鲸骨盆和肋骨^16。上半的地标计算定义的更多细节，这里提出，以试图帮助研究人员可视可能被修改，以适应他们感兴趣的特定对象的步骤。该阴茎骨包含了被利用来自动化某?...

讨论

在上述协议中的关键步骤是:1）解剖的bacula，2）收集显微图像，3）显微输出转换到XYZ的平面文件坐标，4）分割出各试样的点云，5）将每个试样标准化的坐标系，和6）限定的半标。这些步骤很容易被修改，以适应不同的对象。

这些方法可能可以被应用到实质上是"棒状"，或者至少不太弯曲的任何对象。对象曲线回自己成为"U形"不能被目前分析，因为切片（ 图1G）...

材料

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