Method Article

형태 분석 알고리즘은 2D에서 3D 세포 배양의 공간적 형태와 맥락을 정량화하여 표현형의 새로운 정량화를 수행합니다

DOI:

10.3791/71066

July 21st, 2026

In This Article

Summary

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

이 글에서는 세포와 오가노이드의 새로운 형태학적 특징을 선형화된 압축 극좌표(LCPC) 변환을 통해 어떻게 정량화할 수 있는지를 설명합니다. LCPC 변환은 면적, 부피, 표면적과 같은 전통적인 지표로는 포착할 수 없는 특성을 포착합니다. 형태 방향과 정렬에 대한 모범 사례가 논의됩니다.

Abstract

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

수많은 연구에서 세포나 오가노이드의 형태학적 표현형이 항암제에 대한 감수성과 상관관계가 있음을 보여주었습니다. 그러나 전통적인 표현형 측정 방법은 면적, 부피, 둘레, 신호 강도와 같은 공간적 지표에 의존하며, 이는 효과가 있지만 한계가 있습니다. 이러한 접근법은 좌손잡이와 우손잡이의 특성인 카이랄리티와 같은 공간적 맥락의 중요한 특징들을 측정하지 못한다. 부피는 키랄성을 기록하는 데 사용할 수 없는데, 왼쪽과 오른쪽 신발의 부피가 같기 때문입니다. 키랄성, 중력 방향, 극성 축과 같은 공간적 맥락은 인간에게 직관적이지만, 지금까지 세포 생물학자, 병리학자, 방사선학자, 머신러닝 실무자들이 사용하는 전통적인 지표로는 이러한 기본 개념을 포착하지 못합니다. 선형화된 압축 극좌표(LCPC) 변환은 다른 어떤 지표와도 비교할 수 없는 공간적 맥락을 포착할 수 있는 새로운 알고리즘입니다. LCPC 변환은 2차원(2D) 등고선을 이산 사인파 파형으로 변환하여, 등고선과 격자선 사이의 교차점을 추적하는 격자 시스템을 겹쳐 보입니다. 이 방법은 등고선을 독립 좌표 쌍의 연속으로 변환하며, 독립 좌표(x-좌표)는 2차원 공간 내 연속된 위치에 대응합니다. 각 종속 좌표(y-좌표)는 등고선과 격자선의 교차점에서 격자 시스템의 원점 또는 기준선까지의 거리로 구성됩니다. 이산 사인파 형태로 빠른 푸리에 변환이 데이터에 적용됩니다. 이렇게 2D 및 3D 세포 배양에서 세포의 형태가 체계적이고 다차원적으로 표현되어, 치료 저항성에 대한 통찰을 밝히는 견고한 정량적 층화가 가능해집니다.

Introduction

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

3차원(3D) 오가노이드 배양은 생 생물학을 모방하는 데 있어 2차원(2D) 세포 배양보다 우수하다는 것이 입증되었습니다 1,2. 오르가노이드는 효과적인 치료법 선별 검사와 질병 진행에 대한 통찰 확보를 위해 암 연구에 필수적인 존재가 되었습니다 3,4,5,6. 암 생물학에서는 형태학적 구분이 세포와 조직의 뚜렷한 생물학적 행동과 상관관계가 있다는 점은 잘 알려져 있지만, 이 분야는 범위가 제한된 전통적인 형태 지표에 의존하고 있습니다. 이 연구는 면적, 부피, 표면적과 같은 전통적인 접근법을 뛰어넘어, 미묘한 형태학적 특징을 전례 없는 정밀도로 객관적으로 정량화하는 새로운 공간 알고리즘을 설명합니다.

선형화된 압축 극좌표(LCPC) 변환은 조직병리학에서 관찰되는 복잡한 형태를 객관적이고 정량적으로 기술하기 위해 발명되었으며, 결장 폴립7의 등급 향상을 목표로 하였습니다. 이 접근법을 확장해 양극성 장애나 알츠하이머병과 같은 뇌 병리에서 거시적 MRI 기반 공간적 특징을 포착하려는 열망이 있었고, 이를 더욱 강화하는 증강 기술이 나왔습니다. 본 연구는 두 가지 버전의 LCPC 변환인 병렬 그리드 시스템과 방사형 그리드 시스템을 단계별 프로토콜과 모범 사례로 설명하며, 두 버전 모두 오픈 소스 스크립트와 튜토리얼 영상을 제공합니다.

전통적인 기하학적 방법 외에도, 세포와 장기의 복잡한 형태를 측정하는 데 효과적인 여러 추상화가 존재합니다. 여기에는 이심성, 왜곡성, 대비, 크로시스 측정과 함께 원형성, 비대칭성, 가장자리 불규칙성, 전역 등고선 구조 8,9,10을 포착하는 제르니케 모멘트가 포함됩니다. 프랙탈 분석은 복잡한 등고선을 스칼라 값으로 축소하는 데에도 유용했는데, 예를 들어 윤곽선의 불규칙성을 측정하는 프랙탈 차원과 등고선11,12에서 이질적으로 불규칙한 틈이 얼마나 나타나는지를 측정하는 결핍 등이 있다. 질감 측정을 위해 그레이 레벨 동현 매트릭스(GLCM) 방법이 대비, 에너지, 균일성, 상관관계, 엔트로피 측정에 널리 사용된다13,14. 이 모든 방법이 유익하지만, 중력 방향, 왼손잡이와 오른손잡이 비율, 형태 변화 방향에 영향을 미치는 중심 지지 구조의 위치 등 공간적 맥락을 포착하도록 설계된 것은 아닙니다. 더욱이 많은 이들이 공간 정보를 나타내기 위해 단일 스칼라 값 또는 소수의 스칼라 값을 생성하기 때문에, 인용된 연구들은 복잡한 형태를 평가할 때 종종 이들을 조합해 사용합니다.

LCPC 변환은 한 번의 측정에서 많은 특징을 생성하고, 중력 방향과 같은 공간적 맥락을 형태에 인코딩하는 공간 마커를 추가할 수 있도록 설계되었습니다. LCPC 변환7 과 가장 가까운 방법은15년 후에 발표되었으며, 출발점으로 같은 핵심 아이디어를 공유했습니다: 푸리에 변환을 셀 윤곽선에 적용하여 주파수 스펙트럼 형태의 형태로 형태를 측정하는 것입니다. 그러나 LCPC 변환은 측정 도형 외부의 공간적 맥락 지식과 함께 사용하도록 설계된 독립적인 격자 시스템을 통해 적용되도록 개발되었습니다. 더 나아가, LCPC 변환의 발명자는 이 원고에서 결과된 주파수 스펙트럼이 방대한 양의 숨겨진 공간 정보를 포함하고 있다고 설명합니다. 보조 파일 1 에는 인용된 방법을 사용할 때 종종 간과되는 생물학 공간적 맥락과 LCPC 변환을 어떻게 적용하여 이 맥락을 포착할 수 있는지에 대한 광범위한 논의가 포함되어 있습니다. 도형 데이터에서 대조군과 실험군 간에 차이가 없거나, 눈으로 도형 차이를 볼 수 있는지, 또는 도형 데이터에 큰 차이가 없지만 더 큰 차이가 있을 거라 의심하는 독자들은 LCPC 변환을 시도해보라고 권합니다.

LCPC 변환은 다차원 도형을 표현하여 공간 정보 측정에 전례 없는 정밀도를 제공합니다. 면적과 부피와 같은 전통적인 방법들이 각 도형당 하나의 스칼라 값(즉, 25 cm2)만 생성하는 것과 달리, LCPC 변환 결과는 여러 지수를 생성할 수 있으며, 각 지수는 도형의 서로 다른 형태학적 측면(예: 모서리의 둥근함 대 날카로움, 가장자리의 매끄러움 대 톱니 모양)과 상관관계가 있습니다. LCPC 변환은 2D 형태를 설명하기 위해 그 윤곽선과 교차하는 직선 격자를 겹쳐 표현합니다. 각 격자 시스템은 선형 거리를 측정할 수 있는 원점(그림 1) 또는 기준선(그림 2)을 가지고 있습니다. 알고리즘을 설명하는 단순화된 흐름도는 그림 1B2B 를 참조하십시오. 격자선과 도형의 윤곽선 사이의 모든 교차점을 감지합니다. 각 교차점의 거리는 기준선 또는 원점을 기준으로 계산됩니다. 이렇게 LCPC 변환은 2D 도형을 연속된 x-y 좌표 쌍의 연속으로 변환합니다. x-좌표는 0에서 무한대까지의 격자선 위치를 나타내며, y-좌표는 교차점에서 기준선 또는 원점까지의 거리를 나타냅니다. 이 형태는 이산 사인파이며, 빠른 푸리에 변환(FFT)을 적용하여 "위치 영역"에서 주파수 영역으로 데이터를 변환합니다. 만약 x 좌표가 시간을 나타낸다면, "위치 영역"은 FFT의 표준 응용에서 "시간 영역"과 동등하게 됩니다.

방사형 격자 시스템은 각 교차점의 원점까지의 거리를 측정하는 반면(그림 1A), 평행 격자 시스템은 항상 도형 왼쪽의 가상선을 기준으로 교차점 거리를 측정합니다(그림 2A). 오픈소스 스크립트는 윤곽선의 가장 왼쪽 픽셀을 찾아 그 위치에서 10픽셀 왼쪽으로 이동하여 이 가상의 선의 위치를 결정합니다. 여기서 이 위치의 x좌표가 모든 교차점을 계산하는 기준선이 됩니다. 이 10픽셀 규칙은 임의적이지만, 오픈소스 스크립트가 분석하는 모든 등고선은 네 변 모두에 최소 15픽셀의 여백이 있어야 하는 이유입니다.

접힘이나 여러 층이 있는 2D 도형의 경우, 격자선이 윤곽선을 한 번 이상 교차할 수 있습니다. 이 경우 격자선 내 모든 교차점의 거리를 하나의 값으로 합산합니다. 따라서 각 격자선은 단 하나의 x 좌표와 하나의 y 좌표만 가집니다. 이 합은 LCPC 변환이라는 명칭에 포함된 '압축된(compressed)'이라는 용어로 표현됩니다. 비선형 2D 형태를 이산 사인파로 렌더링하는 것은 알고리즘 이름에 "선형화(linearized)"라는 용어로 표현됩니다. 마지막으로, 알고리즘 이름에 '극좌표'라는 용어가 유래했는데, 이는 최초로 개념화된 격자 시스템이 극좌표7의 180도 방사형 격자였기 때문입니다. 극좌표와 데카르트 좌표가 서로 교환 가능하다는 것을 깨달은 후에도 알고리즘 이름은 그대로 유지되었습니다. 감사의 글 섹션에서는 LCPC 변환의 발명과 보강에 동기를 부여한 개인적인 이유를 설명합니다.

프로토콜 1은 그림 3B에 표시된 데이터를 얻는 순서였고, 프로토콜 2는 그림 3C–D에 표시된 데이터를 얻는 순서였습니다. 이 프로토콜들의 단계는 GitHub 저장소 "Pre-Processing-Tools-for-LCPC-Transform"에 있는 개별 Python 스크립트에 구현되어 있습니다. 프로토콜 1은 LCPC 변환을 준비하기 위해 "순수 형태"를 도출하기 위한 전처리 단계의 예입니다. 프로토콜 2는 도형을 "스케일"로 측정하기 위한 전처리 단계의 예로, 즉 서로 상대적으로 원래 스케일에서 측정하는 것입니다.

분할은 Mac의 Preview나 Microsoft Paint와 같은 이미지 처리 소프트웨어에서 수동으로 수행하거나 임계값 기반 방법을 사용할 수 있습니다. 수동으로 할 경우, 다음 네 가지 색상 중 하나를 선택해야 합니다: 파랑, 초록, 분홍/마젠타, 또는 빨강. GitHub 저장소 "Pre-Processing-Tools-for-LCPC-Transform"16 에는 "color extraction scripts"라는 폴더가 있습니다. 이 폴더에는 앞서 언급한 네 가지 색상을 추출해 흰색 배경에 파란색 마스크로 변환하는 Python 코드가 포함된 네 개의 Jupyter Notebook 파일이 들어 있습니다. 이 과정은 Fiji/ImageJ와 같은 이미지 처리 소프트웨어에서 윤곽선 색상을 임계값을 설정하는 방법으로도 수행할 수 있습니다. 경계가 닿거나 겹치는 경계선(예: 벤 다이어그램)이 포함된 윤곽선을 포함하는 이미지의 경우, 이미지 처리 소프트웨어를 사용해 여러 오르가노이드 윤곽선을 포함하는 합성 마스크 이미지 내에서 독립적인 객체가 되도록 분리해야 합니다. 이는 하나의 오가노이드만 포함된 이미지나 서로 접촉하지 않는 여러 오가노이드를 포함하는 이미지에는 적용되지 않습니다.

Protocol

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

본 연구에 사용된 1차 인간 세포는 윤리 승인과 사전 동의 절차에 따라 채취 및 배양되었으며, 이전에 보고된17명이 있습니다. 여기에 제시된 이미지들은 승인된 세포 배양에서 생성되었으며 처음으로 출판되었습니다. 조직 표본은 브리검 여성병원과 캘리포니아 대학교 샌프란시스코 캠퍼스에서 기관심의위원회(IRB) 승인 프로토콜에 따라 연구용 조직 사용에 대해 사전 동의를 제공한 환자들로부터 수집되었습니다. 이 연구를 위해 추가적인 인간 표본 채취는 수행되지 않았습니다.

1. LCPC 변환을 위한 2차원 오가노이드 윤곽선의 순수 형태 추출

  1. Microsoft Paint, MacOS Preview, Fiji/ImageJ 등 이미지 처리 소프트웨어를 사용해 각 오가노이드의 바깥 가장자리를 분할하세요.
  2. 오가노이드의 윤곽을 흰색 배경에 파란 선으로 그려진 열린 마스크로 추출합니다.
  3. 각 파란색 윤곽을 별도의 흰색 캔버스 위에 분리하세요.
    참고: 이 단계는 이미지에 여러 윤곽이 있을 때만 필요합니다.
  4. 각 윤곽을 네 면 모두에 1픽셀 간격으로 자르세요.
  5. 네 면에 100픽셀의 백백을 추가해 새로운 여백을 만듭니다.
  6. 가장 긴 내부 길이가 수평이 되도록 물체를 회전시킵니다.
    참고: 이 단계는 각 객체 내에서 가장 긴 내부 선을 객체의 너비로 변환합니다. 이 목적을 위해 두 가지 스크립트가 있습니다. 닫힌 도형의 경우, "Rotate CLOSED Object Outside by Longest Internal Length_v2.ipynb"라는 Jupyter 노트북을 사용하세요. 열린 도형이나 여러 구성 요소가 있는 도형의 경우, "Rotate Green Line_v3.ipynb"라는 Jupyter 노트북을 사용하세요. Rotate Green Line 도구에는 .ipynb 파일과 같은 폴더 안에 포함된 ReadMe 파일에 링크된 동영상 튜토리얼이 있습니다.
  7. 각 이미지를 파란 물체의 네 면 모두에 1픽셀 여백으로 잘라내세요.
    참고: 이 내용은 Step 1.4와 동일합니다. 하지만 이번에 트리밍하는 목적은 다음에 다룰 크기 조절 단계에서 파란색 물체뿐만 아니라 캔버스 전체의 크기를 조정하는 것입니다. 따라서 파란색 물체를 캔버스 자체와 거의 같은 너비와 높이로 만들면, 캔버스 너비를 400픽셀로 크기 조절하면 객체 크기도 거의 400픽셀 너비로 조정됩니다. "Trim margin to 1-pixel border.ipynb"라는 Jupyter 노트북이 있습니다.
  8. 각 이미지의 캔버스 너비를 400픽셀로 재조정하면서 화면 비율도 제한합니다.
    참고: 이 단계로 각 객체의 너비가 동일해집니다. 종횡비를 제한하면 크기 조정 시 객체가 왜곡되는 것을 방지합니다. "400픽셀로 폭을 재조정하지만 화면 비율을 제한하세요.ipynb"라는 Jupyter 노트북이 있습니다.
  9. 네 면에 100픽셀의 백백을 추가해 새로운 여백을 만듭니다.
  10. 각 이미지를 시계 방향 또는 반시계 방향으로 90도 회전시켜 모든 이미지에서 가장 긴 내부 길이를 수직으로 만듭니다.
    참고: 위 순서에서 순수 형태는 1.6, 1.7, 1.8, 1.9단계를 순차적으로 구현하여 얻습니다. 1.9단계에서 여백 공간을 추가하는 것은 다음에 수행할 LCPC 변환 스크립트의 오류를 방지하기 위한 것이지만, 순수 형태 추출에는 포함되지 않습니다. 도형을 "스케일"로 측정하고자 하는 사용자는 위 순서의 1.6, 1.7, 1.8단계를 수행할 필요가 없지만, LCPC 변환을 적용하기 전에 최적의 방향을 위해 마스크 이미지를 사전 처리해야 합니다.
  11. 각 등고선에 대해 LCPC 변환을 수행하세요. 방사형 그리드 메서드18 또는 병렬 그리드 메서드19에 제공된 Python 스크립트를 사용하세요.

2. LCPC 변환을 위한 2차원 오가노이드 윤곽선의 규모 형태 추출

참고: 객체의 스케일 형태를 단계별로 추출하는 방법은 위의 프로토콜 1에서 순수 형태를 추출하는 방법과 유사합니다. 두 프로토콜 모두에서 매칭을 용이하게 하기 위해 동일한 이름으로 명명되었습니다. 각 단계에 대해 프로토콜 1의 노트를 참고해 주세요. 순수 형태와 실제 크기 형태의 주요 차이점은 크기 형태는 크기 조절 단계인 1.6, 1.7, 1.8단계를 요구하지 않는다는 점입니다.

  1. 각 오가노이드의 바깥 가장자리를 분할합니다. 계산 방식을 사용하거나, 마이크로소프트의 페인트, 맥OS의 미리보기, 피지/이미지J 같은 기본 이미지 편집 소프트웨어를 수동으로 사용할 수 있습니다.
  2. 오가노이드의 윤곽을 흰색 캔버스 위에 파란 선으로 그려진 열린 마스크로 추출하세요. 마스크를 추출하려면 Fiji/ImageJ를 사용하거나 제공된 Python 스크립트16을 사용하세요.
  3. 각 파란색 윤곽을 별도의 흰색 캔버스 위에 분리하세요. 마스크를 추출하려면 Fiji/ImageJ를 사용하거나 제공된 Python 스크립트16을 사용하세요.
  4. 각 윤곽을 네 면 모두에 1픽셀 간격으로 자르세요. Fiji/ImageJ 또는 제공된 Python 스크립트16을 사용하세요.
  5. 네 면에 100픽셀의 백백을 추가해 새로운 여백을 만듭니다. Fiji/ImageJ 또는 제공된 Python 스크립트16을 사용하세요.
  6. 각 윤곽선을 회전시켜 가장 긴 내부 선이 수직이 되도록 합니다. Fiji/ImageJ 또는 제공된 Python 스크립트16을 사용하세요.
  7. 각 등고선에 대해 LCPC 변환을 수행하세요. 방사형 그리드 메서드18 또는 병렬 그리드 메서드19에 제공된 Python 스크립트를 사용하세요.

Results

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

방사형 격자 LCPC 변환을 유방암 오가노이드에 적용한 방법

3D 오가노이드와 종양은 다양한 형태를 취할 수 있으며, 그중 일부는 인간의 눈과는 명백히 다르지만, 전통적인 지표는 통계적으로 유의미하지 않은 차이를 보여줍니다. 반면, 오가노이드는 인간의 눈에 이질적으로 보이는 형태를 보일 수 있어, 반복적으로 나타나는 미묘한 형태를 가리며, 이는 서로 다른 하위 유형을 나타낸다. 인간 유방 일차 유방 오가노이드는 환자에서 추출하여 이전에 발표된 프로토콜17, 20, 21에 따라 3D 배양하였습니다. 그림 3A는 시각적 검사를 통해 접근 가능한 세 가지 형태학적 범주를 보여줍니다: 낭포성 표현형, 다낭성 표현형, 그리고 조밀 다낭성 표현형입니다. 그림 3B는 이 오가노이드들의 외부 윤곽선을 수작업으로 분할하고 순수한 형태를 추출하여 크기가 형태에 미치는 영향을 제거했음을 보여줍니다. 윤곽이 원형이었기 때문에, 순수 도형에 반경 격자 LCPC 변환을 적용한 후 PCA를 적용했습니다. 그림 3B에 겹쳐진 윤곽선들은 가장 긴 내부 선이 수직으로 회전되어 원본 이미지와 방향이 다르다는 점에 유의하세요. 그림 3B는 첫 두 주성분(PC)의 점플롯을 보여주며, K-평균 군집화는 0, 1, 2의 세 개의 군집을 식별했다. 실루엣 점수(SS)는 군집 0과 1(SS = 0.68), 군집 1과 2(SS = 0.69) 사이에 두 개의 별도 군집이 존재할 가능성을 확인시켜 주었으며, 군집 0과 2 사이의 분리가 약하다는 것을 나타냈다(SS = 0.27). 클러스터 0과 2 사이의 실루엣 점수가 약함에도 불구하고, 각 데이터 포인트의 순수한 형태로 데이터 포인트를 표시함으로써 클러스터 2가 가장 독특한 표현형을 나타낸다는 것을 정당화합니다: 매우 비대칭적이고 비원형적인 윤곽선(이 경우 윤곽의 두께는 중요하지 않음).

순수 형태는 크기가 형태에 미치는 영향을 제거하는데, 이는 면적과 부피를 측정할 때 불가능한 것이지만, 규모로 오가노이드 윤곽선을 측정하는 것도 흥미로운 통찰을 제공합니다. 그림 3C 는 크기 조정 없이 이미지에서 분할된 등고선에 대해 수행된 방사형 격자 LCPC 변환으로 생성된 첫 세 PC의 3D 도트 플롯을 보여줍니다. 이 방법은 LCPC 변환이 오가노이드의 형태와 크기를 동시에 측정한다는 것을 의미합니다. 그림 3C 는 다낭성 및 조밀 다낭성 오가노이드가 시낭성 오가노이드에 비해 덜 밀집되어 있음을 보여준다. 이 확장을 정량화하기 위해 그림 3D는 각 표현형이 공유하는 PCA 공간에서 산란된 모습을 보여준다. 각 표현형별로 평균 K-최근접 이웃 거리(KNND)가 계산되며, 이는 낭포성 오가노이드(KNND = 1.061)에서 다낭성 오가노이드(KNND = 1.575), 그리고 조밀한 다낭성 오가노이드(KNND = 2.247)로 점들의 밀도가 감소하는 것을 시각적으로 평가하는 것을 뒷받침합니다; 이 점들의 밀도 감소는 외부 윤곽의 형태학적 이질성 증가로도 해석될 수 있습니다. 마지막으로, LCPC 변환을 전통적인 면적 지표와 비교하기 위해, 그림 3E 는 낭포성 오가노이드와 조밀한 다낭성 오가노이드 면적(p = 0.016, Wilcoxon 랭크-합 검정) 사이에 통계적 유의성이 있음을 보여주지만, 면적에 기반한 세 가지 뚜렷한 형태학적 그룹이 순수 형태 LCPC 분석 후 K-평균 군집화로 식별됨을 나타냅니다(그림 3B)) 또는 대규모 LCPC 분석에서 나타나는 형태학적 이질성이 증가하고 있음을 알 수 있습니다(그림 3C–D). 그림 3E 는 세포성 오가노이드에서 다낭성, 그리고 조밀한 다낭성 오가노이드로 이동하는 데이터 점의 산란이 감소하는 모습을 보여주는데, 이는 대규모 LCPC 분석(그림 3C)에서 관찰되는 것과 반대로, LCPC 변환이 포착하는 공간 정보가 포획하는 영역과 다르다는 점을 강조합니다.

그림 3의 개별 2D 도형과 보충 그림 1 의 개별 형태에도 동일한 접근법이 적용되며, 2D 세포 배양 내 세포의 윤곽에도 적용됩니다: 단일 세포(예: 적혈구) 또는 세포 군집(예: 자궁경부 자궁경부 세포 검사)입니다. 각 사례에 대한 파일럿 연구는 영상으로 설명됩니다: 외부 윤곽선에서 다양한 정도의 접힘을 보이는 뇌 오가노이드의 형태 분석(22); 폐포의 형태 분석은 2D세포 배양 분석에도 동일한 원리를 적용할 수 있습니다. 자궁경부 세포진 검사에서 자궁경부 조직의 형태 분석으로, 핵과 세포막의 윤곽을 방사형 격자 LCPC 변환24를 통해 동시에 정량화합니다; 말라리아에 감염된 적혈구와 감염된 세포 인접하지만 감염되지 않은 세포의 형태 분석25; 겸상적혈구병에 영향을 받은 적혈구의 형태 분석26.

밀도 또는 질감을 공간 정보의 지수로 측정하는 병렬 격자 시스템

공간 정보는 일반적인 형태 개념, 예를 들어 정사각형, 둥근성, 곡률 이상의 의미를 포함하며, 질감과 밀도와 같은 개념을 포함합니다. 혈관과 같은 복잡한 관 모양 네트워크나 세포골격 필라멘트와 같은 가지형 네트워크는 원이나 다각형과 같은 기하학적 직관과 같은 형태 개념과 같은 의미에서 형태 개념에 맞지 않습니다. 하지만 밀도와 질감은 세포, 오가노이드, 종양의 구조와 기능을 이해하는 데 매우 유용한 공간 정보 개념입니다. 형광 영상에서 복잡한 네트워크를 정량화하는 방법을 다루기 위해 세포 생물학자들은 이미지 내 평균 신호 세기나 전체 면적에서 신호가 포함된 비율과 같은 특징을 측정합니다. 이러한 접근법은 공간 정보가 기하학적 모티프로 설명되는 것 이상의 의미가 있음을 인식합니다.

칩 내 조직, 종양, 장기 내 혈관 네트워크 패턴은 전통적인 기하학적 모티프로는 쉽게 정의되지 않지만, 질감과 밀도로 해석함으로써 이점을 얻을 수 있습니다. 종양 발생과 치료 내성에서의 혈관 재형성은 암의 잘 알려진 특징입니다(27번 검토). 458개의 원발성 신경모세포 종양의 조직병리학적 절편을분석하여 전체 혈관과 혈관 구간의 밀도, 크기, 형태를 특성화하였습니다. 그들은 예후가 좋지 않은 환자의 종양에서 혈관이 더 크고, 풍부하며, 더 불규칙한 형태를 보였다는 것을 발견했습니다. 흥미롭게도, '종양 군집을 포괄하는 혈관'(VETC)29라는 용어가 더 큰 간세포암 내에서 작은 종양 군집을 둘러싼 혈관 패턴을 설명하기 위해 만들어졌으며, 이는 전이와 재발률이 높음을 나타냅니다. 이 용어를 만든 동일한 그룹은 나중에 VET 양성 간세포암이 키나제 억제제 소라페닙30 치료에 대한 양성 반응과 상관관계가 있음을 발견했습니다. 따라서 혈관 패턴을 정량화하는 능력은 예후와 내성 기전 이해에 매우 중요합니다.

관, 섬, 가지의 밀집된 네트워크를 정량화하기 위해, LCPC 변환의 병렬 격자 시스템은 형태가 아닌 텍스처 밀도라고 부를 수 있는 것을 측정하는 데 유용합니다. 병렬 격자 LCPC 변환이 이를 구현할 수 있는 예로, 그림 4A 는 심장 혈관 연구에서 나온 심장 미세혈관 패턴을 보여줍니다.31: 육안으로는 경각된 부위와 손상 부위에서 멀리 떨어진 부위 간의 패턴을 객관적이고 정량적으로 구분하는 것이 매우 어렵습니다. 혈관의 마스크(그림 4A 1단계)를 추출하고 병렬 격자 LCPC 변환을 적용하면 정량적 차이를 관찰할 수 있습니다. 첫째, 데이터의 청결성을 높이기 위해 각 이미지는 균일한 직경의 원을 자르고, 이 원은 정사각형 이미지 중앙에 위치하게 하여 회전과 관계없이 각 이미지가 동일한 높이와 너비를 갖도록 했습니다(그림 4A 단계 2). 다음으로, 잘린 이미지를 회전시켜 혈관의 주축( 그림 4A 1단계의 빨간 이중 화살표)이 수직으로 정렬되도록 합니다(그림 4A 단계 3). 그 다음, 이미지 둘레에 파란색 선을 그어, 혈관 패턴이 드문드문한 이미지도 파란색 픽셀 직경이 같도록 합니다(그림 4A 4단계 4). 3단계와 4단계는 Fiji/ImageJ 소프트웨어(v2.14.0)를 사용하여 수행되었습니다. 마지막으로 병렬 그리드 LCPC 변환이 적용됩니다(그림 4A 단계 5). PCA를 통해 결과를 분석하고 처음 세 개의 PC를 그래프로 표시하면, 손상된 부위가 먼 지역과 별도로 집결되어 있음을 알 수 있습니다(그림 4B). 편면수 거리(CD) 계산은 그래프가 나타내는 시각적 패턴을 정량화합니다. 손상 후 1일째에 경결 부위는 기저 상태에서 더 멀리 떨어진 혈관 패턴을 보인다(기저 대 1Day_infarct: CD = 83.89), 비손상 부위(기저 대 1Day_remote: CD = 53.95). 부상 후 1일째 시점에 부상하지 않은 원격 부위는 기저 상태보다 부상 부위에 더 가깝습니다(1Day_remote 대 1Day_infarcxt: CD = 34.48). 손상 후 3일째에 손상 부위와 관련 비손상 원격부위의 혈관 패턴이 기저 상태와 서로 크게 다르게 나타납니다(기저 대 3Days_Infarct: CD = 240.63; 기저 3Days_Remote: CD = 67.28; 3Days_Remote 대 3Days_Infarct: CD = 168.35). 매우 흥미롭게도, 부상 후 7일째에 부상 부위와 부상하지 않은 원격 부위의 혈관 패턴은 기저 상태와 거의 직교인 3차원 공간을 차지하며, 부상 1일 후보다 더 가까운 거리에서 정렬되어 있습니다(7일차 CD = 24.76, 1일차 34.48). 데이터는 두 가지를 시사합니다: (1) 치유가 진행됨에 따라 손상 부위와 손상되지 않은 원격부위의 혈관 패턴이 기저 상태와는 구별되는 패턴을 채택한다; (2) 손상 부위에서 활성화된 재건 프로그램은 관련 비손상 원격지에서도 유사한 형태학적 변화를 유발한다.

LCPC 변환의 다재다능성을 평가하기 위해 혈관의 면적(5A)과 총 둘레(그림 5B)라는 전통적인 지표가 계산되었습니다. 면적과 주변은 시간 흐름 내 각 지점 간의 변화를 보여주지만, 선박의 뚜렷한 공간적 배열을 일관되게 정량화하지는 못합니다. 면적과 둘레는 다양한 혈관 배열을 나타내는 추상화로, 치료 조건과 시간 경과에 따라 혈관의 공간적 배열이 어떻게 변하는지에 대한 정보를 제공하지 못합니다. 각 그룹에 대해 n=3일 때, t-검정과 윌콕슨 검정은 신뢰할 수 없다. 클리프의 델타는 비모라미터 방식이며 효과 크기의 변화를 잘 보여줍니다. 혈관 면적을 측정함으로써 드러나는 주요 차이는 기저 상태와 7일 시점 사이이며; 혈관 전체 둘레의 차이는 기저 상태와 1일 상태, 기저 상태 대 7일 상태 간의 차이입니다. 그러나 경계와 경계는 감염된 부지와 부상당한 외딴 부지 간에 차이가 있다는 징후를 제공하지 않는데, 이는 감염된 부위가 고의로 손상된 점을 고려하면 매우 놀라운 일입니다. 여기서 LCPC 변환의 다양성이 전통적인 지표와 비교해 드러납니다. 각 샘플에 대해 LCPC 변환의 일부로서 FFT 결과는 주파수 집합과 연관된 크기 집합을 생성합니다. 각 빈도 빈의 크기는 평균화하여 처리 그룹 간 비교할 수 있습니다(보완 그림 1). 이를 통해 처리 그룹 간 가장 큰 차이를 보이는 빈도 빈을 간단히 평가할 수 있습니다. 그림 5C 는 다양한 치료 그룹 간, 특히 울폐된 부위와 관련 비손상 부위 간 효과 크기 차이가 분명한 6개의 빈을 보여주며, 해당 부위의 지역과 경계는 차이가 없음을 나타냅니다. 이 6개 빈의 독립적인 값들은 처리 조건 간 양적 차이를 보여주는 여러 인덱스를 제공합니다. 흥미롭게도, LCPC 변환 결과의 Bin 1에서 그룹 간 패턴은 전체 둘레 측정 결과와 유사하다. 반면 다른 지수는 면적과 주변 모두에 숨겨진 양적 차이를 드러냅니다. 기저 상태에 고유한 지수: Bin 49, Bin 56, Bin 21. 1일차 그룹 간 고유한 지수는 다음과 같습니다: Bin 1, Bin 6, Bin 49. 3일차 그룹에만 고유한 지수는 다음과 같습니다: Bin 6, Bin 112, Bin 49. 7일차 그룹에만 고유한 지수: Bin 6, Bin 49, Bin 112. 이 데이터는 LCPC 변환이 전통적인 면적과 경계 접근법에 비해 얼마나 다재다능하고 정밀한지를 보여줍니다.

figure-results-1
그림 1: LCPC 변환의 방사형 격자 시스템 설명. (A) 1단계: 관심 있는 형태를 분할하여 파란색 마스크로 추출합니다. 2단계: 마스크를 위한 경계 상자(즉, 물체를 네 면에서 감싸고 닿는 가장 작은 직사각형)를 결정합니다. 경계 상자 중심에서 시작된 방사형 격자 선 시스템이 마스크 위에 겹쳐집니다. 원점으로 물체의 중심이 아닌 바운딩 박스의 중심을 선택한 결정은 임의적이었습니다. 3단계: 격자선과 도형 사이의 각 교차점 거리를 원점에 대해 계산합니다. 교차점이 여러 개의 격자선은 그 모든 교차 거리를 하나의 값으로 합산합니다. 4단계: FFT를 적용하여 3단계에서 얻은 데이터를 주파수 그래프로 변환하며, 이는 1단계의 형태를 정적 다차원 표현(5단계)으로 나타냅니다. 6단계: 5단계부터 인덱스를 통해 스칼라 값을 생성할 수 있습니다. 계산할 지수 유형에 대한 논의는 그림 10 을 참조하세요. (B) 방사형 그리드 시스템에 대한 LCPC 변환 단계의 단순화된 설명을 보여주는 플로우차트. (C) 방사형 격자 시스템이 최적의 형태 테마의 예시. 도형들은 마이크로소프트 파워포인트로 그렸습니다. 이 이미지들은 저자가 제작했으며 출판 라이선스나 저작권 허가가 필요하지 않습니다. 이 그림의 더 큰 버전을 보시려면 여기를 클릭해 주세요.

figure-results-2
그림 2: LCPC 변환의 병렬 그리드 시스템의 설명. (A) 1단계: 관심 있는 형태를 수동으로 분할하여 이진 마스크로 추출합니다. 2단계: 평행 격자 시스템을 도형에 겹쳐 놓고, 도형은 회전하여 베이스가 수평이 됩니다. 3단계: 격자선과 도형 사이의 각 교차점 거리를 기준선에 대해 계산합니다. 교차점이 여러 개의 격자선은 그 모든 교차 거리를 하나의 값으로 합산합니다. 4단계: FFT를 적용하여 3단계에서 얻은 데이터를 주파수 그래프, 즉 1단계의 도형을 정적으로 다차원으로 표현하는 그래프로 변환합니다. 6단계: 5단계부터 인덱스를 통해 스칼라 값을 생성할 수 있습니다. 계산할 지수 유형에 대한 논의는 그림 10 을 참조하세요. (B) 병렬 그리드 시스템의 LCPC 변환 단계를 단순화한 설명을 제공하는 플로우차트. (C) 예시는 평행 격자 시스템이 최적의 주제를 형성한다. 도형들은 마이크로소프트 파워포인트로 그렸습니다. 이 이미지들은 저자가 제작했으며 출판 라이선스나 저작권 허가가 필요하지 않습니다. 이 그림의 더 큰 버전을 보시려면 여기를 클릭해 주세요.

figure-results-3
그림 3: LCPC 변환은 유방암 오가노이드의 형태학적 아형을 객관적으로 식별합니다. (A) 1차 인간 유방 오가노이드는 이전에 보고된17, 20, 21보충 파일 1에 설명된 대로 유도되었다. 시각적 검사를 통해 세 가지 표현형이 관찰되고 분류될 수 있습니다: 낭포성, 다낭성, 그리고 조밀 다낭성. 스케일 바 = 50 μm. (B) 오가노이드의 외부 윤곽선을 수동으로 분할하여 (보충 그림 2; 보충 파일 1에 설명된 기준), 이 윤곽선들의 순수한 형태가 얻어졌습니다. 도표 내 윤곽선들은 모두 회전되어 가장 긴 내부 선이 수직임을 유의하세요; 따라서 원래 방향과 일치하지 않습니다. 그 후 방사형 격자 LCPC 변환을 적용했고, 결과를 PCA를 통해 분석되었습니다. K-평균 군집화는 세 개의 군집(0, 1, 2)을 확인했으며, 실루엣 점수(SS)를 계산하여 그 존재를 확인했습니다. 클러스터 0 대 1, SS: 0.68; 클러스터 0 대 2, SS: 0.27; 클러스터 1 대 클러스터 2, SS: 0.69. (C) 오가노이드의 외곽 윤곽선도 크기 조절 없이 크기 도형, 즉 서로 상대적으로 나타나는 도형에 대해 방사형 격자 LCPC 변환을 사용하여 측정되었습니다. 삽입 부분은 낭포성 오가노이드와 조밀한 다낭성 오가노이드의 느슨한 패킹을 보여주기 위해 확대되었습니다. B와 C의 경우: 초록색 = 낭포성, 보라색 = 다낭성 성, 노란색 = 조밀 다낭성 혈통. (D) 처음 세 개의 PC 값은 각 표현형별로 별도로 그래프로 표시되며, 평균 KNND는 낭포성 표현형이 가장 밀집되어 있는 반면, 다낭성 및 조밀 다낭성 표현형은 밀도가 낮음(즉, 더 이질성)을 보인다는 것을 나타냅니다. 각각 따로 그래프를 그리지만, 세 표현형 모두 동일한 PCA 공간을 공유합니다. 낭포성 표현형의 경우, 평균 KNND는 삽입 지점에서 계산되었는데, 이 표현형에만 존재하는 이상치가 계산에 영향을 미치기 때문입니다. 평균 KNNDs: 낭포성 1.061; 다낭성 1.575; 조밀 다낭성 2.247. (E) 오르가노이드 면적의 박스 플롯은 시포성 오가노이드에서 다포성 오가노이드, 그리고 조밀한 다낭성 오가노이드로 데이터 점이 점차 분산되고 있음을 보여주는데, 이는 대규모 LCPC 분석에서 관찰되는 것과 반대입니다(그림 3C). 이는 LCPC 변환이 포착한 공간 정보가 해당 영역이 포착하는 것과 다르다는 점을 강조합니다. 이 그림의 더 큰 버전을 보시려면 여기를 클릭해 주세요.

figure-results-4
그림 4: 평행 격자 LCPC 변환을 이용한 텍스처-밀도 측정. (A) 심장 조직의 미세혈관 사진은 이전에 다른 연구에서 얻은31,32. 이 이미지들은 이 연구들에서 직접 복제된 것이 아닙니다; 오히려 Gkontra가 제공한31개의 원래 3D 영상 데이터셋에서 생성된 새로운 뷰를 나타냅니다. 비록 이 데이터셋은 이전에 공개된 바 있지만, 여기 표시된 특정 이미지들은 처음으로 공개된 것입니다. 각 시점에서 부상하지 않은 원격 현장의 이미지는 공개되지 않습니다. 1단계: 2D 마스크를 생성하고, 혈관의 일반적인 정렬 축(빨간색 이중 화살표)을 결정합니다. 2단계: 이미지 중앙을 중심으로 원형 영역을 얻어 다음 단계에서 이미지가 다양한 정도로 회전된 후 각 이미지가 같은 높이를 가지도록 강제합니다(회전 정도에 따라 이미지의 높이가 변합니다). 3단계: 각 이미지를 1단계에서 결정한 정렬 축이 수직으로 회전되도록 합니다. 4단계: 각 이미지 주위에 균일한 파란색 테두리를 적용하여 이미지의 직경이 같은지 확인합니다. 밀도가 희박한 이미지는 더 짧을 수 있기 때문입니다(예: 3단계의 "심근경색 후 7일" 참조). 5단계: 병렬 그리드 LCPC 변환을 적용합니다. (B) 평행 격자 LCPC 변환 출력의 첫 세 PC의 점플롯 (보라색: 기저 상태, 초록색: 울색 부위, 노란색: 경색 부위에서 멀리 떨어진 부위). 세 조건 사이의 각 시점에서의 챔퍼 거리 계산에 관한 결과 섹션을 참조하세요. 이 그림의 더 큰 버전을 보시려면 여기를 클릭해 주세요.

figure-results-5
그림 5: LCPC 변환은 영역과 경계에 숨겨진 특징을 측정합니다. 전통적인 형태 지표와 비교하기 위해 (A) 혈관 네트워크의 면적백분율과 (B) 전체 둘레를 측정하였다. (C)LCPC 변환 결과에서 얻을 수 있는 정밀도 를 보여주기 위해, FFT로 생성된 개별 주파수 빈의 값을 비교했습니다(보충 그림 1). 치료 그룹 간 차이를 객관적으로 정량화하기 위한 지수로 여섯 개의 빈도가 선택되었습니다. 모든 패널에서 각 그룹의 표본 크기가 제한적(n = 3) 때문에, 효과 크기에 대한 비모수 클리프스 델타가 윌콕슨 랭크-합 검정의 p값 대신 계산되었습니다. 이 그림의 더 큰 버전을 보시려면 여기를 클릭해 주세요.

figure-results-6
그림 6: 깨끗한 데이터를 추출하기 위해 도형을 체계적으로 배치하는 것의 중요성. (A) 서로 거울상인 동일한 도형의 예시와 평행 격자 시스템이 이들의 뚜렷한 방향을 어떻게 결정하여 서로 다른 결과를 내는지. (B) 도형 2를 수평으로 뒤집어 방향이 도형 1과 일치하도록 하여 도형 1과 더 공정하게 비교할 수 있게 하는 예시. 이 구성에서 병렬 격자 LCPC 변환은 형태 1과 형태 2 모두에서 동일한 결과를 제공합니다. 이 방향 변경은 Shape 2를 수정하지 말아야 한다는 공간적 맥락이 없을 때만 이루어져야 합니다. (C) 숫자 "3"과 비슷하게 생긴 네 가지 변형 물체의 예시로, 일부는 잘못된 방향을 향하고 있습니다. (D) 병렬 그리드 LCPC 변환으로 얻을 데이터의 청결성을 극대화하기 위해 각 객체를 어떻게 방향화할 수 있는지에 대한 예시입니다. (E) 패널 D와 동일한 방향 규칙에 따라 패널 C의 물건들의 넓은 돌출부를 측정한 예시. (F) 패널 C의 물건들의 큰 전구를 측정하되 그 정렬 맥락을 유지한 예시. 패널 D에 정렬 규칙을 유지하는 동기는 넓은 돌출부의 공간 정보 기여도와 좁은 돌출부의 기여도를 분석하기 위함이다. C–F와 관련된 FFT 스펙트럼은 보충 그림 3 을 참조하십시오. A–B의 이미지는 Microsoft PowerPoint에서 가져온 아이콘입니다. C-F의 도형들은 Microsoft PowerPoint에서 그렸습니다. 이 이미지들은 저자가 제작했으며 출판 라이선스나 저작권 허가가 필요하지 않습니다. 이 그림의 더 큰 버전을 보시려면 여기를 클릭해 주세요.

figure-results-7
그림 7: 전통적인 형태 지표가 명확한 맥락을 인식하지 못하는 문제. (A) 체크 표시의 일련은 같은 모양이지만 왼쪽에서 오른쪽으로 크기가 커집니다. 이 시리즈에서는 각 체크마크의 면적이 다르지만 모양 개념은 동일합니다: 즉, 체크마크 스타일은 동일합니다. (B) 정사각형은 중심을 중심으로 45도 회전하여 그 형태 개념이 인간이 말하는 다이아몬드가 됩니다. 정사각형의 형태 개념은 다이아몬드와 매우 다르지만, 그 사이의 회전은 주관적으로 '기울어진 사각형' 또는 '기울어진 다이아몬드'로 해석됩니다. LCPC 변환은 이 급수의 각 단계별 차이를 객관적이고 정확하게 측정할 수 있습니다. 도형들은 마이크로소프트 파워포인트로 그렸습니다. 이 이미지들은 저자가 제작했으며 출판 라이선스나 저작권 허가가 필요하지 않습니다. 이 그림의 더 큰 버전을 보시려면 여기를 클릭해 주세요.

figure-results-8
그림 8: 순수 형태와 실제 축척 형태의 설명. "스케일"이라는 것은 물체의 원래 크기에 따라 서로 상대적으로 형태를 측정하는 것을 의미합니다. 반면 '순수 형태'는 크기가 형태에 미치는 영향을 제거하여 서로 다른 물체의 형태를 비슷한 규모로 측정하는 것을 의미합니다. 비교하는 객체들을 객관적으로 크기 조정하는 방법을 찾는 것이 중요합니다. 가장 긴 내부 선이 종종 최적의 선택입니다(Step 1). 두 물체 모두 가장 긴 내부 선이 수평이 되도록 회전시킵니다(Step 2). 이 경우 두 물체도 회전하여 움푹 들어간 면이 아래를 향하게 됩니다. 다음으로, 각 객체의 폭은 이제 각 객체의 가장 긴 내부 선과 같은 크기가 되며, 크기를 400픽셀로 재조정하고, 객체가 왜곡되지 않도록 종횡비를 제한합니다(3단계). 마지막으로, LCPC 변환을 적용하기 전에 크기 변경된 객체를 90도(4단계) 회전시키는 것이 선택 사항입니다(이 경우 반시계 방향). 이 그림의 더 큰 버전을 보시려면 여기를 클릭해 주세요.

figure-results-9
그림 9: 암 연구를 위한 공간 마커 예시. (A) 2D 배양에서 세포가 이동 방향으로 라멜리포디아를 확장하여 이동하는 가상의 예시. 상피세포의 이동 형태학 맥락에서는, 각 세포의 윤곽선을 이동 축과 일치하는 가장 긴 내부 선이 수직으로 정렬되도록 배치하는 것이 공간적 맥락을 포착하는 데 가장 합리적입니다. (B) 형태가 변하지 않는 외피 내에서 자라는 종양 또는 전암성 병변의 가상적 예시. 병변은 중심 지지 구조인 줄기에서 자랍니다. 병변이 성장하거나 시들면서 윤곽은 인피 내부 벽 쪽으로 또는 멀어지며 이동합니다. 점선은 인케이스의 얼마나 많은 부분을 따라 나가야 하는지(주황색 선)를 합리적이고 객관적으로 정의하는 방식을 나타내며, 너무 적거나 너무 많은 추가 정보가 추가되지 않도록 합니다. (C) 옵션 1은 병변의 윤곽선에 적용된 평행 격자 LCPC 변환을 보여준다. 옵션 2도 동일하게 보여주지만, 인케이스먼트가 위치한 공간적 맥락을 포함합니다. 예를 들어, 빨간 점은 격자선(검은 화살표, 초록 선)과 윤곽선의 교차점을 나타냅니다. 옵션 2는 각 격자선이 옵션 1에 없는 네 번째 교차점을 만들어내며, 이는 옵션 2의 결과에 정량적으로 공간적 맥락을 더하지만, 옵션 1에는 전혀 없는 맥락이 됩니다. 도형들은 마이크로소프트 파워포인트로 그렸습니다. 이 이미지들은 저자가 제작했으며 출판 라이선스나 저작권 허가가 필요하지 않습니다. 이 그림의 더 큰 버전을 보시려면 여기를 클릭해 주세요.

figure-results-10
그림 10: 가상의 FFT 결과 및 이로부터 도출할 수 있는 일반적인 지수들. (A) 두 가상의 도형에 대한 FFT에서 도출된 빈도 그래프가 표시된다. Bin 3과 Bin 4의 크기는 반대 경향을 보이므로, 이 두 빈 간의 비율은 두 형태 간의 차이를 나타내는 좋은 지수가 됩니다. 빨간색과 초록색 점선은 좋은 지수를 도출하는 데 사용할 수 있는 다른 주파수 빈을 나타냅니다. (B) 두 가상의 형태 형태의 차이를 나타내는 FFT 프로필이 단계적으로 감소하는 추세를 보이는 히스토그램. (C) 주파수 빈 구간의 중앙 영역에서 FFT 프로필이 가장 크게 다르게 나타나는 두 가상의 도형의 히스토그램. (D) 그 존재와 상관관계가 있는 지수를 생성할 수 있는 가상의 형태학적 특징의 예시. 도형들은 마이크로소프트 파워포인트로 그렸습니다. 이 이미지들은 저자가 제작했으며 출판 라이선스나 저작권 허가가 필요하지 않습니다. 이 그림의 더 큰 버전을 보시려면 여기를 클릭해 주세요.

보조 그림 1: 혈관 네트워크의 LCPC 변환 결과 전반에 걸쳐 다양한 주파수 빈의 크기를 비교한 박스-위스커 플롯. 그림 45에 나타난 혈관 네트워크의 LCPC 변환은 각 이미지를 나타내는 주파수 빈과 그 크기들을 제공합니다. 각 빈의 크기의 박스-위스커 플롯을 사용하여 모든 처리 그룹을 비교하며, LCPC 변환을 통해 측정된 질감-밀도 층화의 지표로 가장 적합한 빈을 식별합니다. 각 박스-위스커 플롯은 n = 3개의 독립적인 이미지로 구성됩니다. 이 파일을 다운로드하려면 여기를 클릭해 주세요.

보충 그림 2: 분절형 유방 종양 오가노이드 예시. 오르가노이드는 전문가가 수동으로 분절하여 정기적으로 배양하고 현미경으로 관찰했습니다. 이 그림은 분절된 오가노이드 등고선의 예시를 보여준다. 모든 현미경 이미지의 스케일 바는 50 μm입니다. 이 연구에 사용된 1차 인간 세포는 윤리 승인과 사전 동의 절차에 따라 채취 및 배양되었으며,21명이 보고되었습니다. 여기에 제시된 이미지들은 승인된 세포 배양에서 가져온 것이며, 이번이 처음으로 공개되었습니다. 추가적인 인간 표본 채취는 이루어지지 않았습니다. 이 파일을 다운로드하려면 여기를 클릭해 주세요.

보조 그림 3: 그림 6의 도형에 대한 FFT 스펙트럼. 그림 6의 도형에 대한 FFT 결과는 독자들이 도형 방향 변화가 공간 정보에 미치는 영향을 이해할 수 있도록 제공한다. 그래프의 왼쪽 열은 변경되지 않은 y축을 가지며, 오른쪽 열은 축이 잘려 있어 작은 빈 값을 시각화할 수 있습니다. (A) 그림 6C에 표시된 네 가지 도형의 평행 격자 LCPC 변환 스펙트럼으로, 재방향 조정이 없음을 의미합니다. (B) 그림 6D의 형태를 평행 격자 LCPC 변환에서 얻은 스펙트럼으로, 텍스트에 설명된 대로 재배치되었습니다. (C) 그림 6E에 나오는 형태의 병렬 격자 LCPC 변환 스펙트럼. (민주당) 그림 6F의 도형과 비교하기 위해 평행 그리드 LCPC 변환의 스펙트럼을 사용합니다. 이 파일을 다운로드하려면 여기를 클릭해 주세요.

보조 파일 1: LCPC 변환의 모범 사례에 대한 광범위한 세부 정보. 이 보조 파일은 모범 사례와 추가 방법론적 노트에 대한 광범위한 논의를 제공합니다. 문서의 개요는 다음과 같습니다: 1절 - 면적, 부피, 둘레, 거칠기가 포착할 수 없는 공간적 맥락. 섹션 2 - LCPC 변환 사용 및 결과 해석에 대한 모범 사례. 3절 - 자연의 고유한 극성을 이해하여 LCPC 변환을 극대화합니다. 4장 - LCPC 변환으로 생성된 FFT 스펙트럼 해석. 5장 - LCPC 변환으로 측정하기 전에 도형의 방향을 맞추는 방법. 섹션 6 - "순수 형태"와 "규모 형태"는 무엇인가요? 섹션 7 - 통계 방법에 관한 방법론적 노트. 8장 - 세포 배양 및 오가노이드의 분절에 관한 방법론적 노트. 섹션 9 - 자신만의 LCPC 변환을 코딩하는 사람들을 위한 FFT 매개변수. 이 파일을 다운로드하려면 여기를 클릭해 주세요.

Discussion

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

모든 계산 도구와 마찬가지로, 출력의 품질은 입력이 도구가 설계된 예상 규칙을 준수하는지에 달려 있습니다. 제공된 Python 스크립트를 사용할 때 품질 관리 검사로 다음 규칙을 준수하고 사용해야 합니다. 먼저, 열린 마스크가 흰색 배경에 파란색 선으로 그려져 있는지 확인하세요. 파란색은 임의의 결정이며 의미가 없지만, LCPC 변환을 수행하는 스크립트들은 흰색 배경에서 파란색 픽셀을 찾도록 작성되었습니다. 둘째, 네 면에 10픽셀 미만의 백백 공간이 있는 이미지에는 LCPC 변환 스크립트를 실행하지 마세요. 그 이유에 대한 설명은 이 영상을 참고하세요: "여진 두께 매개변수 - 도형 개요 준비 방법"33. 셋째, 이 스크립트는 객체의 윤곽선이 흰색 배경 위의 파란색 픽셀이라고 가정합니다. 선의 두께는 5에서 10 픽셀 사이여야 합니다. 그 이유는 이 영상을 참고하세요: "선 두께 매개변수 - 도형 아웃라인 준비 방법"34. 넷째, 빈 이미지에 LCPC 변환 스크립트를 실행하지 마세요. 이로 인해 스크립트에 오류가 발생할 수 있습니다. 다섯째, 잔해 조각만 있는 이미지에 대해 LCPC 변환 스크립트를 실행하지 마세요. 이로 인해 격자 시스템이 점점을 감지하지 못해 빈 이미지를 분석하려 시도할 수 있어 오류가 발생합니다. 여섯째, 원치 않는 이물질이 포함된 이미지에 대해 LCPC 변환 스크립트를 실행하지 마세요. 이물질은 열린 마스크 객체의 높이와 너비를 변화시켜 부정확한 결과를 초래합니다. 일곱째, LCPC 변환으로 객체를 분석하기 전에 왜 그 크기를 이해해야 하는지 이해하기 위해 이 영상을 참고하세요: "스케일 매개변수 - 도형 개요 준비법"35.

LCPC 변환의 가장 놀라운 점 중 하나는 공간적 맥락을 포착하는 능력입니다. 생물학은 자연을 이해하는 데 필수적인 공간적 맥락으로 가득 차 있습니다. 예를 들어, 상피세포는 첨단 기저 극성이 36,37,38,39입니다. 더욱이, 미세중력 조건에서 변화하는 중력 방향은 우주비행사 생리학 40,41과 암세포 형태42,43,44 모두에 중요한 영향을 미칩니다. 마지막으로, 현장 내 전신포형성과 신생물증은 혈액 공급에 접근하는 중심 지지 구조 역할을 하는 줄기를 통해 내강으로 성장할 수 있습니다. 따라서 병변들이 후퇴할 때는 하내강에서 멀어지면서 줄기 쪽으로 모두 퇴행합니다45,46. LCPC 변환은 공간적 맥락에 매우 민감하여 사용자가 주의하지 않으면 추가적인 이질성을 초래할 수 있습니다. 그림 6은 LCPC 변환을 적용하기 전에 도형을 체계적으로 정렬하는 방법을 설명합니다. 그림 7은 전통적인 기하학이 인간에게 명백한 형태 개념을 측정할 수 없음을 설명합니다. 그림 8은 순수 형태 측정과 축척 형태 측정의 차이를 시각화한 도표를 제공합니다. 그림 9는 2D 배양, 조직학, 3D 영상의 2D 단면에서 형태를 정량화할 때 공간적 맥락이 종종 간과된다는 점을 설명합니다. 분량 제약으로 인해 보충 파일 1에는 5절의 그림 6그림 7, 그리고 6절의 그림 8그림 9에 대한 광범위한 논의가 포함되어 있습니다.

FFT는 다차원 스펙트럼 데이터를 생성하기 때문에, 사용자가 스펙트럼 내 반복되는 패턴에서 지수를 도출하는 데 익숙해지는 것이 중요합니다. 이 논문에서는 LCPC 변환을 FFT 적용 단계를 포함한다고 언급하지만, 실제 LCPC 변환은 FFT 단계 이전에 이루어집니다(그림 1A–B 4단계, 그림 2A–B 3단계 3). 그러나 FFT를 적용하면 LCPC 변환 출력의 주파수 영역에서 객관적이고 체계적인 통찰을 얻을 수 있습니다. 따라서 FFT를 최종 단계로 포함하지 않으면 LCPC 변환이 효과적이지 않을 것입니다. LCPC 변환에서 FFT는 일련의 빈과 그 크기들을 생성하는데, 이들은 모두 1단계의 윤곽선 입력을 정적으로 다차원으로 표현합니다. 정적인 다차원 정보 표현의 비유는 미국에서 개인의 전화 통신 채널을 나타내는 9자리 전화번호입니다. 도형의 스칼라 값 표현(예: 길이 85cm, 너비 56도 등)에 비해 번거롭지만, 다차원 데이터는 통찰을 추출하는 데 훨씬 더 많은 정보를 담고 있습니다. 그림 10 은 LCPC 변환을 적용할 때 관찰된 여러 스펙트럼 패턴과 이를 통해 지수를 도출하는 방법을 보여줍니다. 보조 파일 1 의 4절에서는 그림 10 을 더 자세히 다룹니다.

LCPC 변환의 주파수 영역은 육안으로는 쉽게 인지할 수 없는 2차원 형상의 객관적이고 정량적인 특징을 드러냅니다. 곡률, 톱니 모양, 비틀림의 매우 미묘한 차이도 견고하게 정량화할 수 있습니다. 그러나 이 방법으로 형태학적 하위 유형을 식별하면, 양적 차이와 상관관계가 있는 시각적으로 뚜렷한 주제를 관찰할 수 있습니다. 이 기능은 LCPC 변환을 머신러닝(ML) 모델의 공간 정보 추출 방법으로 해석 가능하게 하며, 결과 ML 모델 자체도 더 해석 가능하게 만듭니다. 신경망은 강력한 ML 방법이지만, 그 결과가 인간에게 해석되지 않는 경우가 많아 의료 인공지능(AI) 분야에 있어 큰 문제입니다. 의사들은 모델이 왜 잘못된 결정을 내리는지 이해하지 못하면 임상 의사결정에 ML 모델을 사용하는 것을 꺼립니다. 더불어, 어떤 인간도 왜 또는 어떻게 모델이 편향된 결정을 내렸는지 이해하지 못한다면, 이를 해결하는 것은 어려워집니다.

질감-밀도 측정에 관한 논의(그림 4 와 5)에 덧붙이자면, 인간의 눈은 착시나 무질서한 패턴에 쉽게 혼란스러워진다. 복잡한 패턴의 무질서한 배열은 종종 무작위적이거나 측정할 가치가 없다고 해석됩니다. LCPC 변환은 겉보기에는 무작위 패턴 속에서도 재현 가능하고 객관적이며 정량적인 차이를 측정할 수 있습니다. 하지만 사용자가 무작위로 위장한 혼돈되고 무질서한 자연의 패턴에 속지 않도록 훈련받지 않는 한, 예측 가능한 패턴을 밝혀내기 위해 무질서를 측정하려는 시도는 이루어지지 않습니다. 예측 가능한 패턴은 시스템이 완전히 무작위적이지 않다는 명확한 신호입니다. LCPC 변환은 육안으로는 미묘하게 다르거나 상상할 수 없을 정도로 다른 형태를 측정하기 위해 발명되었습니다. 자연의 많은 형태들이 예측 가능한 무질서의 하위 유형이라기보다는 의미 없는 혼돈 이종성에 가까운 미세한 차이를 보이기 때문에, 각 하위 유형은 서로 다른 생물학적 기능과 연관되어 있기 때문에, LCPC 변환과 함께 사용될 "4단계" 개념 틀이 개발되었습니다. 첫째, 복잡해 보이는 시스템이 무작위라는 생각에 저항하세요. 둘째, 혼돈적 이질성을 '무질서'라고 부르고, 합리적 기준 조건이 될 수 있는 질서 유형을 찾아보세요. 장애가 기준에서 벗어나고 있다고 가정하세요. 셋째, 참조 패턴을 객관적으로 측정할 방법을 찾으세요. 예를 들어 LCPC 변환입니다. 넷째, 이 동일한 객관적 방법을 시스템 내 장애로 간주되는 것을 측정하는 데 적용한다.

Disclosures

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

저자들은 공개할 이해 상충이 없습니다. LCPC 변환은 BrainScanology, Inc.에 의해 소프트웨어 제품 "Shape Genie"로 상업화되었으나, 특허 출원이 철회되었고 회사는 해산되었습니다. 이 연구는 BioRxiv에 사전 인쇄본으로 제공되며, 참고문헌47에 인용되어 있습니다. 그림 4 는 Gkontra가 제공한31 개의 3D 영상 데이터셋에서 생성된 이전에 공개되지 않은 뷰를 보여줍니다; 우리는 그들을 정당하게 인정했습니다. 그림 3보충 그림 1 은 공동 저자가 보고한21 건의 윤리 승인 및 사전 동의 절차에 따라 설정된 1차 인간 세포 배양에서 파생된 미공개 이미지를 보여줍니다. 이 연구를 위해 추가적인 인간 표본 채취는 수행되지 않았습니다.

Acknowledgements

Loading...
$$\rightleftharpoonup{xx}$$ $$\longleftharp{xx}$$, $$\longrightharp{xx}$$,

우리는 대장암으로 별세한 고등학교 과학 교사 듀안 니콜스를 기리고자 합니다. 이 연구는 결장 폴립의 형태를 특성화하기 위한 LCPC 변환의 발명에 영감을 주었습니다. 둘째로, 양극성 장애 2형을 앓았던 투안 트린에게 감사드리며, 공간적 맥락을 포착하는 마커 추가와 같은 LCPC 변환 개선에 영감을 주셨습니다. LCPC 변환은 비공식적으로 Nguyen-Nichols-Trinh(NNT) 변환이라고 불립니다. 셋째, LCPC 전환을 상업화하려 한 스타트업을 지원해 주신 수석 엔젤 투자자 폴 리얼에게 감사드립니다. 넷째, 이 연구를 위해 심장 혈관 사진의 사진사진을 공유해 주신 Polyxeni Gkontra 박사님께 감사드립니다. 다섯째, 조직 샘플 지원과 국립보건원(R01CA281361)의 조직 샘플 및 자금 지원을 제공해 주신 DF/HCC 유방 SPORE: 전문 연구 우수성 프로그램(SPORE)(NCI 1P50CA168504), UCSF 유방 치료 센터 외과 및 병리학 팀, 유방 치료 센터 인턴들에게 감사드립니다.

Materials

List of materials used in this article
NameCompanyCatalog NumberComments
Add 100 pixels to margins.ipynb조직 공간 기하학 연구실GitHub 저장소 이름: Pre-Processing-Tools-for-LCPC-Transform | GitHub 저장소 주소: https://github.com/dh2nguyen/Pre-Processing-Tools-for-LCPC-TransformJupyter Notebook의 Python 스크립트
Average K-means Nearest Neighbor via sklearn version 1.6.1SciKit Learnhttps://scikit-learn.org/stable/modules/generated/sklearn.neighbors.NearestNeighbors.htmlsklearn.neighbors 패키지 내의 NearestNeighbors 함수 사용
BLUE_Color Extraction.ipynb조직 공간 기하학 연구실GitHub 저장소 이름: Pre-Processing-Tools-for-LCPC-Transform | GitHub 저장소 주소: https://github.com/dh2nguyen/Pre-Processing-Tools-for-LCPC-TransformJupyter Notebook의 Python 스크립트
Chamfer Distance via scipy version 1.15.1SciPyhttps://docs.scipy.org/doc/scipy/reference/generated/scipy.spatial.cKDTree.htmlscipy.spatial 패키지 내의 cKDTree 함수 사용
Cliff's Delta조직 공간 기하학 연구실https://github.com/dh2nguyen/Cliffs-Delta파이썬 언어로 Cliff's Delta의 구현 작성
Enlarge Image Size by 3X.ipynb조직 공간 기하학 연구실GitHub 저장소 이름: Pre-Processing-Tools-for-LCPC-Transform | GitHub 저장소 주소: https://github.com/dh2nguyen/Pre-Processing-Tools-for-LCPC-TransformJupyter Notebook의 Python 스크립트
Fiji/ImageJ v2.14.0ImageJ.nethttps://imagej.net/
Flip images HORIZontally.ipynb조직 공간 기하학 연구실GitHub 저장소 이름: Pre-Processing-Tools-for-LCPC-Transform | GitHub 저장소 주소: https://github.com/dh2nguyen/Pre-Processing-Tools-for-LCPC-TransformJupyter Notebook의 Python 스크립트
Flip images VERTically.ipynb조직 공간 기하학 연구실GitHub 저장소 이름: Pre-Processing-Tools-for-LCPC-Transform | GitHub 저장소 주소: https://github.com/dh2nguyen/Pre-Processing-Tools-for-LCPC-TransformJupyter Notebook의 Python 스크립트
GREEN_Color Extraction.ipynb조직 공간 기하학 연구실GitHub 저장소 이름: Pre-Processing-Tools-for-LCPC-Transform | GitHub 저장소 주소: https://github.com/dh2nguyen/Pre-Processing-Tools-for-LCPC-TransformJupyter Notebook의 Python 스크립트
Isolate Objects Empty Blue.ipynb조직 공간 기하학 연구실GitHub 저장소 이름: Pre-Processing-Tools-for-LCPC-Transform | GitHub 저장소 주소: https://github.com/dh2nguyen/Pre-Processing-Tools-for-LCPC-TransformJupyter Notebook의 Python 스크립트
MacOS Preview Version 11.0 Apple Inchttps://support.apple.com/guide/preview/welcome/mac
Master Script_Parallel Grid LCPC_DN_071025.ipynb조직 공간 기하학 연구실GitHub 저장소 주소: https://github.com/dh2nguyen/Radial-Grid-LCPC-TransformJupyter Notebook의 Python 스크립트
Master Script_Radial Grid LCPC_bBoxCentered.ipynb조직 공간 기하학 연구실GitHub 저장소 주소: https://github.com/dh2nguyen/Parallel-Grid-LCPC-TransformJupyter Notebook의 Python 스크립트
Microsoft PaintMicrosofthttps://www.microsoft.com/en-us/windows/tips/paint수동 분할에 유용한 도구로 언급되었으나 이 연구에서는 사용되지 않음
Microsoft Powerpoint Version 16.109.3Microsofthttps://www.microsoft.com/en-us/microsoft-365/p/powerpoint/cfq7ttc0pbmc
PINK_Color Extraction.ipynb조직 공간 기하학 연구실GitHub 저장소 이름: Pre-Processing-Tools-for-LCPC-Transform | GitHub 저장소 주소: https://github.com/dh2nguyen/Pre-Processing-Tools-for-LCPC-TransformJupyter Notebook의 Python 스크립트
Principal Component Analysis (PCA) via sklearn version 1.6.1SciKit Learnhttps://scikit-learn.org/stable/api/sklearn.decomposition.htmlsklearn.decomposition 패키지의 PCA 함수 사용
RED_Color Extraction (Mac).ipynb조직 공간 기하학 연구실GitHub 저장소 이름: Pre-Processing-Tools-for-LCPC-Transform | GitHub 저장소 주소: https://github.com/dh2nguyen/Pre-Processing-Tools-for-LCPC-TransformJupyter Notebook의 Python 스크립트
RED_Color Extraction (Windows).ipynb조직 공간 기하학 연구실GitHub 저장소 이름: Pre-Processing-Tools-for-LCPC-Transform | GitHub 저장소 주소: https://github.com/dh2nguyen/Pre-Processing-Tools-for-LCPC-TransformJupyter Notebook의 Python 스크립트
Resize width to 400 pixels but constrain aspect ratio.ipynb조직 공간 기하학 연구실GitHub 저장소 이름: Pre-Processing-Tools-for-LCPC-Transform | GitHub 저장소 주소: https://github.com/dh2nguyen/Pre-Processing-Tools-for-LCPC-TransformJupyter Notebook의 Python 스크립트
Rotate CLOSED Object Horizontal by Longest Internal Length.ipynb조직 공간 기하학 연구실GitHub 저장소 이름: Pre-Processing-Tools-for-LCPC-Transform | GitHub 저장소 주소: https://github.com/dh2nguyen/Pre-Processing-Tools-for-LCPC-TransformJupyter Notebook의 Python 스크립트
Rotate Green Line_v3.ipynb조직 공간 기하학 연구실GitHub 저장소 이름: Pre-Processing-Tools-for-LCPC-Transform | GitHub 저장소 주소: https://github.com/dh2nguyen/Pre-Processing-Tools-for-LCPC-TransformJupyter Notebook의 Python 스크립트
Rotate images 90CCW.ipynb조직 공간 기하학 연구실GitHub 저장소 이름: Pre-Processing-Tools-for-LCPC-Transform | GitHub 저장소

Reprints and Permissions

Request permission to reuse the text or figures of this JoVE article

Request Permission

Tags

Bioengineering3D Organoidscell cultureShape AnalysisSpheroidsOrgans on a ChipLCPC Transform
Video Coming Soon

Related Articles