초파리 안구 시스템은 다양한 생물학적 과정, 특히 인간 신경 퇴행성 질환을 연구하는 데 유용한 도구입니다. 그러나 거친 눈 표현형의 수동 정량화는 편향되어 신뢰할 수 없을 수 있습니다. 여기에서는 ilastik과 Flynotyper를 사용하여 편견 없는 방식으로 눈 표현형을 정량화하는 방법을 설명합니다.
방법 논문
* These authors contributed equally
초파리 안구 시스템은 다양한 생물학적 과정, 특히 인간 신경 퇴행성 질환을 연구하는 데 유용한 도구입니다. 그러나 거친 눈 표현형의 수동 정량화는 편향되어 신뢰할 수 없을 수 있습니다. 여기에서는 ilastik과 Flynotyper를 사용하여 편견 없는 방식으로 눈 표현형을 정량화하는 방법을 설명합니다.
Drosophila melanogaster compound eye는 대칭 및 육각형 패턴을 나타내는 약 800개의 ommatidia로 구성된 잘 구조화되고 포괄적인 배열입니다. 이러한 규칙성과 관찰의 용이성으로 인해 Drosophila eye system은 다양한 인간 신경 퇴행성 질환을 모델링하는 강력한 도구입니다. 그러나 눈 심각도 점수의 수동 순위와 같은 비정상적인 표현형을 정량화하는 방법은 특히 눈 형태학의 약한 변화를 순위를 매길 때 제한이 있습니다. 이러한 한계를 극복하기 위해 Flynotyper와 같은 계산 접근 방식이 개발되었습니다. 링 라이트를 사용하면 개별 옴마티디아의 온전함에 접근하는 더 나은 질적 이미지를 얻을 수 있습니다. 그러나 이러한 이미지는 링 라이트에 의해 도입된 옴마티디아의 그림자로 인해 Flynotyper에서 직접 분석할 수 없습니다. 여기에서는 ilastik과 Flynotyper라는 두 가지 소프트웨어를 결합하여 초파리 질병 모델에서 관찰된 거친 눈 표현형을 편향되지 않게 정량화하는 방법을 설명합니다. ilastik으로 이미지를 전처리하면 Flynotyper로 거친 눈 표현형을 성공적으로 정량화할 수 있습니다.
Drosophila melanogaster 게놈은 인간 질병 관련 유전자 orthologs의 ~ 75 %를 포함합니다. 또한 초파리 눈 발달 중에는 게놈 유전자의 약 2/3가 발현되어 초파리 눈은 다양한 분자 및 세포 기능, 발달 및 질병 모델을 조사하는 뛰어난 유전 시스템입니다 1,2. 따라서 초파리의 눈 시스템은 다양한 생물학적 과정을 연구하는 데 유용한 실험 도구입니다.
Drosophila 겹눈은 대칭적이고 육각형 패턴을 나타내는 ~ 800 개의 옴 마티 디아 (ommatidia)의 잘 구조화되고 포괄적 인 배열입니다3. 이 육각형 패턴의 규칙성은 눈 형태에서 돌연변이 및 유전자 발현 변화의 영향을 추정하는 데 사용될 수 있다4. 눈 형태 평가가 필요한 이전 연구는 육안으로 감지된 눈 표현형의 중증도에 대한 수동 순위에 크게 의존했습니다. 눈 표현형의 순위를 매기기 위해, 외부 눈 형태학 이미지는 실체현미경 5,6에 의해 촬영된다. 각 그룹의 눈 표현형은 외부 눈을 4개의 영역으로 분할하고 각 영역의 퇴행 비율을 계산하여 평가합니다 5,6. 그런 다음, 이 값들은 대조군 파리7로부터 얻어진 값과 비교되는 평균을 계산하는데 사용된다. 점수는 융합 정도, 옴마티디아 손실 및 강모 조직 7,8을 기반으로 합니다. 실체현미경으로 촬영한 파리 눈 사진은 한 명의 연구원이 획득하고, 눈 표현형 분석은 다른 연구자가 삼중 검증 세트 7,8을 사용하여 수행합니다.
육안으로 눈의 형태가 약한 변화를 평가하는 데는 한계가 있다4. 이러한 한계를 극복하기 위해 FLEYE 및 Flynotyper와 같은 계산 접근 방식이 개발되었습니다 1,9. Flynotyper는 초파리 안구 시스템(Drosophila eye system)의 형태학적 변화를 정량적으로 추정하기 위한 새로운 계산 방법입니다1. Drosophila eye 및 개별 ommatidium을 자동으로 감지하여눈의 불규칙성을 기반으로 표현형 점수(P-Score)를 계산합니다 1. P-점수가 높을수록 파리의 눈이 더 퇴화되었음을 나타냅니다. 이 소프트웨어는 초파리 눈10의 이상을 정량화하는 데 성공적으로 사용되었습니다. Flynotyper는 자동화된 프로세스를 보장하지만 다양한 광학 현미경 방법으로 촬영한 일부 눈 이미지에는 여전히 성공적으로 적용할 수 없습니다.
질적으로, 우리는 각 옴마티듐을 더 정확하게 표현하기 때문에 단일 지점 광원에 비해 링 광원을 선호합니다. 그러나 링 라이트를 사용하면 옴마티듐의 반구형 모양으로 인해 각 옴마티듐의 상단에 고리 모양의 그림자가 생성됩니다. 이 고리 모양의 그림자는 Flynotyper에 의한 정확한 ommatidial 감지를 방해하여 P-Score의 잘못된 계산을 유발합니다.
이러한 문제를 극복하기 위해 다양한 분석을 위한 머신러닝 기반 도구인 ilastik을 구현하여 플라이 아이 이미지(fly eye image)에서 옴마티디아(ommatidia)를 분류했다11. 그런 다음 ilastik에서 생성된 이미지를 Flynotyper에 공급하여 P-Score를 계산했습니다. 이를 통해 Drosophila eye morphological defects를 편견 없이 정량화할 수 있습니다1.
1. 정량화 준비
2. ilastik을 사용하여 플라이 아이 이미지에서 옴마티디아 감지
3. ImageJ를 사용하여 Flynotyper용 사진 준비
4. Flynotyper를 사용하여 표현형 점수 계산
이전 연구에서는 이 프로토콜을 사용하여 전두측두엽 치매(ALS-FTD)를 동반한 근위축성 측삭 경화증과 관련된 돌연변이 VCP 단백질의 유전적 변형을 확인했습니다12. 또한, 이 방법은 구형 실체현미경13을 사용하는 경우에도S59L 매개 ALS-FTD CHCHD10 독성을 평가하기 위해 다른 논문에서도 사용되었습니다. 이러한 결과를 더욱 검증하기 위해 Baek et al.13과 동일한 ALS-FTD 질환 돌연변이에 ilastik과 Flynotyper를 사용했습니다. 그러나 훈련을 받지 않은 실험실 직원이 실험을 수행하여 유사한 결과를 얻었습니다. 따라서 ilastik과 Flynotyper를 결합하는 것은 Drosophila ommatidia의 변성을 편견 없이 점수화하는 강력한 방법입니다.
우리는 더 나은 정성적 평가를 가능하게 하기 때문에 링 라이트의 사용을 선호합니다(그림 3). 그러나 링 라이트를 사용하면 Flynotyper가 직접 분석할 수 없는 그림자가 발생합니다. 이를 극복하기 위해 Flynotyper로 P-Score를 정량화하기 전에 ilastik에서 이미지를 전처리했습니다. 전반적으로, 비퇴화(GMR-GAL4) 파리 눈의 P-점수는 중등도 퇴행된 눈(2x_S81L) 및 심하게 퇴행된 눈(3x_S81L)에 비해 유의하게 낮았습니다(그림 4A, B 및 표 1). 이는 Flynotyper에서 생성된 P-점수가 높을수록 눈이 더 퇴화되었음을 나타내기 때문에 예상되는 결과입니다1. ilastik을 사용하지 않았을 때, 퇴행되지 않은 눈과 적당히 퇴행된 눈 사이의 P-Score에는 무시할 수 있는 정도의 차이가 있었습니다(그림 4C 및 표 1). 이는 Flynotyper가 ilastik을 사용하지 않고 링 라이트 조건에서 옴마티디아 패턴을 정확하게 식별하는 데 어려움을 겪고 있음을 시사합니다. 이는 Flynotyper가 옴마티얼 중심1을 식별하기 위해 추가 조정을 통해 회색조로 변환된 이중 점 조명으로 촬영한 이미지를 사용하도록 설계되었기 때문입니다. Ilastik은 각 옴마티듐의 중앙에 점을 통합하여 이 프로세스를 복제하여 Flynotyper 내에서 정량화를 가능하게 합니다.

그림 1: ilastik 워크플로우의 예. (A) ilastik에서 새 프로젝트 생성. (B) 표준 이미지 선택. (C) Feature Selection에서 시그마 값 선택. (D) 민감도를 위해 옴마티디아의 전경을 표시합니다. (E) 민감성을 위해 옴마티디아의 배경을 표시합니다. (F) 소프트웨어가 처리된 파일을 내보내는 방법을 선택합니다. (G) 분석할 실험용 플라이 아이 이미지 선택. (H) 어떤 실험용 플라이 아이 이미지를 분석할 것인지 보여줍니다. (I) 소프트웨어 실행 후 완성된 제품. 이 그림의 더 큰 버전을 보려면 여기를 클릭하십시오.

그림 2: Flynotyper 워크플로우의 예. (A) 추가 분석을 위해 파리의 눈만 선택합니다. (B) 새 창에서 자른 눈 이미지를 엽니다. (C) 분석할 절단 이미지와 분석을 위해 확인할 매개변수를 선택합니다. (D) 예시 결과. 표현형 점수는 맨 오른쪽에 표시됩니다. 이 그림의 더 큰 버전을 보려면 여기를 클릭하십시오.

그림 3: 정성적 초파리 눈 이미지. 링 라이트(왼쪽)로 촬영한 동일한 플라이 아이와 퇴화되지 않은(위) 및 퇴화된(아래) 조건에서 듀얼 포인트 라이트(오른쪽)로 촬영한 사진. 이 그림의 더 큰 버전을 보려면 여기를 클릭하십시오.

그림 4: Drosophila의 거친 눈 표현형을 정확하게 정량화하기 위해 Flynotyper 전에 ilastik을 사용합니다. (A) 퇴화되지 않은 눈(GMR-GAL4), 중등도의 퇴행성 눈(2x_S81L) 및 심하게 퇴행된 눈(3x_S81L)의 대표적인 초파리 눈 이미지. (B) ilastik 사용 후 Flynotyper의 정량 결과. 양측 치환 t-검정의 P-값은 결과의 해석을 용이하게 하기 위해 계산되었습니다(GMR-GAL4 대 2x_S81L: p = 0, GMR-GAL4 대 3x_S81L: p = 0) (C) ilastik을 사용하지 않고 Flynotyper의 정량화 결과. 양측 치환 t-검정의 P-값은 결과 해석을 용이하게 하기 위해 계산되었습니다(GMR-GAL4 대 2x_S81L: p = 0.597, GMR-GAL4 대 3x_S81L: p = 0.687). 각 그룹에 대해 N = 3입니다. 데이터 포인트는 시각적 명확성을 위해 x축에 지터링되거나 누적됩니다. 수평 실선은 평균을 나타냅니다. 오차 막대는 표준 편차를 나타냅니다. 이 그림의 더 큰 버전을 보려면 여기를 클릭하십시오.
| GMR-갈4 | 2x_S81L | 3x_S81L | |||
| ilastik 함유 | ilastik 미함유 | ilastik 함유 | ilastik 미함유 | ilastik 함유 | ilastik 미함유 |
| 26.7 | 55.7 | 52.5 | 51.2 | 59.9 | 55.6 |
| 17.4 | 53.5 | 50 | 54.9 | 60.7 | 52 |
| 22.8 | 55.5 | 56.8 | 55.2 | 60.1 | 55.1 |
표 1: Flynotyper를 ilastik과 함께 사용할 때와 Flynotyper를 단독으로 사용할 때의 P-점수.
추가 그림 S1: 이미지 획득을 위해 설정된 플라이 헤드 마운트. 이것은 여러 장의 용지를 사용하여 이미지 획득을 설정하는 방법의 예입니다. 그러나 머리와 플라이 아이의 올바른 방향을 보장하는 모든 방법을 사용할 수 있습니다. 이 파일을 다운로드하려면 여기를 클릭하십시오.
보충 그림 S2: 기능 선택에서 시그마 값 2.5를 사용할 때 ilastik 및 Flynotyper의 정량적 결과. ilastik의 기능 선택에서 5.0 이외의 시그마 값을 사용하면 부정확한 정량화가 발생할 수 있습니다. 특히 시그마 값이 2.5이면 민감도가 손실됩니다. 양측 치환 t-검정의 P-값은 결과 해석을 용이하게 하기 위해 계산되었습니다(GMR-GAL4 대 2x_S81L: p = 0.192, GMR-GAL4 대 3x_S81L: p = 0). 각 그룹에 대해 N = 3입니다. 데이터 포인트는 시각적 명확성을 위해 x축에 지터링되거나 누적됩니다. 수평 실선은 평균을 나타냅니다. 오차 막대는 표준 편차를 나타냅니다. 이 파일을 다운로드하려면 여기를 클릭하십시오.
보충 그림 S3: 괴사성 패치로 인한 민감도 손실 및 부정확한 P-Score. 흰색 화살표는 괴사한 반점의 일부를 가리키고 있습니다. 이 파일을 다운로드하려면 여기를 클릭하십시오.
보충 표 S1: 기능 선택에서 시그마 값 2.5를 사용할 때의 P-점수. 이 파일을 다운로드하려면 여기를 클릭하십시오.
초파리의 옴마티디아는 다양한 생물학적 기능과 유전 질환을 연구하는 데 유용한 시스템으로 구성됩니다. 옴마티디아의 규칙성은 유전적 돌연변이의 영향을 조사하기 위한 좋은 척도이다4. 수동 순위와 같이 ommatidial 규칙성을 계산하는 몇 가지 방법이 존재하지만 이러한 방법은 심하게 편향될 수 있습니다. 이러한 편향된 접근 방식을 극복하기 위해 반자동 도구가 개발되었습니다 1,11.
정확한 눈 정량화를 위해서는 파리 준비 및 설정을 최적화해야 합니다. 각 실험실마다 설정이 있을 것이라는 것을 알고 있지만, 독자에게 이미지 획득을 위해 파리를 준비하는 방법에 대한 아이디어를 제공하기 위해 실험실에서 다음 설정을 제공하고 있습니다. 파리는 -80 °C에서 냉동 보관됩니다. 표본 크기 3을 사용한 이유는 우리의 경험에 따르면 GMR-GAL4가 S81L 돌연변이와 결합할 때 모집단 간에 변동이 거의 없기 때문입니다. 그러나 귀무가설(null-hypothesis) 유의성 검정을 수행하려면 표본 크기를 조정해야 합니다. 이미지 획득을 위해 파리를 장착하기 위해 파리 머리는 베개 역할을 하는 매우 작은 종이 위에 놓여 눈의 중심이 대략 초점면이 되도록 합니다(보충 그림 S1). 우리 연구실에서는 Z-stack 현미경과 카메라( 재료 표 참조)를 사용하여 전체적으로 더 선명한 이미지를 생성하며, 파리의 눈이 돔 모양이기 때문에 눈 전체를 가장 잘 시각화할 수 있습니다. 그러나 Z-stack 기능13이 없는 구형 실체현미경과 함께 이 방법을 성공적으로 사용했기 때문에 Z-stack 기능은 이 프로토콜에 대한 요구 사항이 아닙니다. 밝기를 위한 링 라이트 다이얼은 과다 노출을 피하기 위해 1/4 이상 위로 돌리지 않습니다. 샘플의 총 두께는 300에서 450 미크론 사이입니다.
이 연구에서 우리는 ilastik 및 Flynotyper를 사용한 파리 눈 분석이 편향되지 않고 신뢰할 수 있는 옴마티디아 정량화 결과를 생성한다는 것을 보여주었습니다. 머신러닝 기반 생물영상 분석 도구인 ilastik을 사용하여 Drosophila ommatidia 패턴을 분석하여 ommatidia를 표시한 다음, ommatidia 규칙성 1,11을 측정하는 데 사용되는 계산 방법인 Flynotyper를 사용하여 분석했습니다. 우리는 더 나은 정성적 평가를 위해 이미지를 얻기 위해 링 라이트를 사용하는 것을 선호합니다. 그러나 이 접근 방식으로는 Flynotyper가 이러한 이미지를 사용하여 정량적 결과를 생성할 수 없습니다. 따라서 ilastik은 이러한 이미지를 전처리하기 위해 구현되었으며, 이러한 처리된 이미지는 Flynotyper와 함께 성공적으로 사용되어 정량적 결과를 생성했습니다. Flynotyper 이전에 이미지를 ilastik으로 전처리하지 않았을 때, 전체 P-점수는 비퇴행성 눈과 퇴행성 눈 간에 유사했습니다. 이는 정량적 결과를 위해 링 라이트를 구현할 때 ilastik을 먼저 사용하는 것의 중요성을 강조합니다.
이 방법에는 몇 가지 제한 사항이 있습니다. ilastik의 설정은 크게 다를 수 있으며 유의하지 않은 데이터를 생성할 수 있습니다. 예를 들어, 동일한 이미지를 사용하여 다른 분석 세트를 수행했지만 ilastik의 전경 시그마 값을 5 대신 2.5로 조정했습니다. 이로 인해 중등도의 퇴행성 표현형에서 민감도가 손실되었습니다(보충 그림 S2 및 보충 표 S1). 또 다른 한계는 너무 많은 퇴행으로 인한 눈의 감도 손실입니다. 3x_S81L 눈과 비교했을 때, 괴사성 반점이 있는 눈은 질적으로 훨씬 더 퇴화되어 있습니다. 그러나 3x_S81L 눈의 평균 P-점수는 60.2인 반면, 괴사된 반점이 있는 눈의 평균 P-점수는 60.9였습니다. 따라서 괴사성 반점과 같이 퇴행이 심한 눈도 감도 상실을 초래할 수 있습니다(보충 그림 S3).
이 프로토콜의 중요한 단계는 ilastik을 훈련시키는 데 사용되는 적절한 표준 이미지를 선택하는 것입니다. 옴마티디아와 주변 색소 세포(예: 적색 또는 주황색 눈) 사이에 뚜렷한 대비가 있는 퇴화되지 않은 표준 이미지는 Flynotyper에서 보다 정확한 정량화를 제공합니다. 또 다른 중요한 단계는 ilastik의 Feature Selection으로, 시그마 값으로 표시되는 감도를 선택하여 기계 학습 소프트웨어를 미세 조정합니다(그림 1D). 너무 적은 기능을 선택하면 ilastik의 감도가 감소하여 ommatidia의 부정확한 마킹으로 이어질 수 있습니다. 너무 많은 피처를 선택하면 결과를 생성하는 데 필요한 시간과 처리 능력이 증가하여 절차의 효율성이 떨어집니다. 또는 눈에 띄는 정성적 퇴화가 있지만 분석 결과가 통계적으로 유의하지 않은 경우 선택한 특징의 시그마 값을 늘리는 것이 잠재적인 수정이 될 수 있습니다. 그러나 이 프로토콜에서 제공하는 기능 선택은 정기적으로 정확한 결과를 생성합니다. 설명된 대로 중요한 단계를 따르면 향후 가장 일반적인 문제와 문제 해결을 완화하는 데 도움이 됩니다.
요약하면, 이 방법은 더 나은 정성적 이미지를 제공하는 동시에 초파리 눈 시스템에 대한 편향되지 않고 신뢰할 수 있는 정량화를 제공합니다. 결과적으로, 이 방법은 미래에 퇴행성 질환을 모델링하기 위한 강력한 도구로서 눈 퇴행의 유용성을 향상시킬 수 있습니다.
저자는 공개할 이해 상충이 없습니다.
이 프로토콜에 사용된 현미경과 카메라를 사용해 주신 Pedro Fernandez-Funez에게 감사드립니다. 또한 프로토콜의 명확성에 대한 피드백을 제공해 주신 Ava Schapman에게도 감사드립니다. 재정 지원은 Wallin Neuroscience Discovery Fund에서 김남철에게 제공되었습니다.
| 이름 | 회사 | 카탈로그 번호 | 댓글 |
|---|---|---|---|
| 컴퓨터 사양 | Ryzen 5, 16GB RAM, Nvidia RTX 3070 Super, Windows 10 | ||
| Flynotyper | Iyer, J. et al. (2016) | 여기에서 소프트웨어 다운로드: https://flynotyper.sourceforge.net/imageJ.html | 오픈 소스 소프트웨어. Flynotyper 2.0을 사용하지 마십시오. 발표 당시 2.0은 상당히 새로운 것이었고 이 프로토콜은 Flynotyper의 원래 버전에 최적화되어 있습니다. |
| ilastik | Berg, S. et al. (2019) | 여기에서 소프트웨어 다운로드: https://www.ilastik.org/download.html | 오픈 소스 소프트웨어. 컴퓨터 운영 체제에 해당하는 Regular Builds에서 버전 1.4.0.post1을 다운로드하십시오. |
| ImageJ | 여기에서 소프트웨어를 다운로드하십시오 : https://imagej.net/ij/download.html | 오픈 소스 소프트웨어. 버전 1.53 및 1.54가 사용되었습니다. 1.54는 업데이트된 버전이며 기본 다운로드입니다. | |
| Leica Application Suite (LAS X) | Leica Microsystems | LASX Office 1.4.6 28433 | z-stack 획득에 사용되는 시스템 및 소프트웨어. |
| DMC2900 카메라가 있는 Leica Z16 APO 현미경 | Leica Microsystems | 10 447 173, 12 730 466 | 텍스트에서는 Z-stack 현미경 및 카메라라고 합니다. 이 제품은 현재 보관되어 있습니다. |
이 JoVE 논문의 텍스트 또는 그림 재사용 허가 요청
허가 요청