Method Article
여기, 우리 프로토콜 및 관련 만성 스트레스에 노출 된 후 변경 내용을 확인을 쥐 모델을 사용 하 여 메 틸 Seq, epigenomic 플랫폼의 구현을 설명 합니다. 결과 쥐 메 틸-Seq 플랫폼은 쥐에서 스트레스 노출에서 발생 하는 메 틸 화 차이 감지 수를 보여 줍니다.
다양 한 동물의 게놈을 사용할 수 있게,이 동물 모델에 동적 epigenetic 변화를 캡처할 수 있는 도구에 대 한 증가 필요가 있다. 쥐 epigenetic 도구 통찰력 기계 정보를 제공 하는 많은 약리학 및 행동 연구를 보완 수 있는 하나의 특정 모델 동물입니다. 이 위해, 우리는 쥐 쥐 게놈에서 DNA 메 틸 화 수준을 평가할 수 있는 SureSelect 캡처 시스템 (메 틸-Seq 라고 함)를 적응. 쥐 디자인 대상 발기인, CpG 섬, 섬 해안, 및 모든 RefSeq 유전자에서 GC 풍부한 지역.
쥐 실험에는 플랫폼을 구현 하기 위해 남성 Sprague Dawley 쥐 혈액 샘플 게놈 DNA 추출에 대 한 수집 된 후 3 주 동안 만성 변수 스트레스에 노출 되었다. 메 틸-Seq 라이브러리는 전단, 어댑터 결 찰, 대상 농축, 황산 변환 및 멀티플렉싱 하 여 쥐 DNA 견본에서 건설 되었다. 라이브러리는 다음-세대 시퀀싱 플랫폼에서 시퀀싱 했다 고 시퀀스 읽기 강조 하 고 강세가 없는 쥐의 DNA 사이 DMRs를 식별 하기 위해 분석 했다. 최고 후보자 DMRs 했다 독립적으로 유효성을 검사 하지 황산 pyrosequencing 플랫폼의 견고성을 확인 하 여.
결과 쥐 메 틸-Seq 플랫폼은 스트레스에 노출에 의해 유도 된 메 틸 화 변화를 캡처할 수 있는 유용한 epigenetic 도구를 보여 줍니다.
높 처리량 연속에 발전은 다양 한 모델 및 비 모형 유기 체 게놈 시퀀스에 이르렀다. 같은 시퀀스의 여부 유전학, 비교 유전체학, 및 transcriptomics 연구를 촉진 했다. 예를 들어, 사용 가능한 게놈 시퀀스는 매우 풍요롭게 히스톤 수정1또는 황산 시퀀싱와 연결에 따라 DNA 칩 Seq 실험에서 시퀀싱 데이터를 정렬 하는 데 유용는 DNA 메 틸 화에 의해 측정 unmethylated cytosines2의 황산 변환에서 형성 uracil 감지. 그러나, 유전자 기능에 영향을 미칠 수 있는 종의 규제 시퀀스의 주석된 데이터의 부족 때문에 그들의 디자인에서 사용 가능한 게놈 시퀀싱 데이터를 통합 하는 epigenomic 플랫폼의 구현에서 지연 되었습니다.
특히, DNA 메 틸 화 methylomic 플랫폼을 구축 하기 위한 사용 가능한 게놈 데이터를 활용할 수 있는 DNA에 가장 널리 공부 후 성적인 수정 중 하나입니다. 1 개의 그런 보기는 인간 methylome3, 정신과4,5종양학에서 다양 한 분야에서 널리 사용 되었습니다 배열 기반 플랫폼입니다. 불행 하 게도, 비인간적인 동물 모델에 대 한 비슷한 플랫폼은 거의 아무 널리 사용 플랫폼을 그들의 초기 디자인에서 게놈 시퀀스의 장점을 받아들으로 부족, 있습니다.
비 인간 동물 모델의 methylomic 프리를 평가 하는 일반적인 방법은 이다 감소 표현 황산 시퀀싱 (RRBS)6. 이 접근 방식을 제공 하는 포괄적인 methylomic 풍경, 낮은 읽기-심층 취재 비용 및 게놈2의 큰 유전자 가난한 지역에서 제한 된 기능 정보를 제공 하는 전체 게놈 황산 시퀀싱 비용 극복 . RRBS 제한 다이제스트 및 일반적으로 유전자 발기인 근처 발견 하 고 유전자 규칙7에 역할을 생각 하는 CpG 섬 등 매우 다양 한 GC 시퀀스에 대 한 풍부 하 게 genomic DNA의 크기 선택 포함 한다. 중요 한 연구의 숫자에 RRBS 메서드를 사용 하는 동안 제한 효소에 의존이 아니다 주목할 만한 도전과 한계 없이. 예를 들어, GC 풍부한 시퀀스 RRBS의 농축은 전기 이동 법에 의해 제한 효소와 후속 크기 선택에 의해 인식 하는 특정 순서의 존재에 전적으로 의존. 즉, 이러한 금지 사이트를 포함 하지 않는 모든 게놈 영역 크기 선택 시 제외 됩니다. 또한, 크로스-종 비교 동일한 제한 사이트는 다른 종 중 동일한 loci에 존재 하지 않는 도전 하.
RRBS의 한계를 극복 하는 한 가지 방법은 사용 하는 플랫폼의 설계에 게시 된 게놈 시퀀스를 활용 하는 농축 방법. 어레이 기반 인간 플랫폼 뇌관 프로브 대립 유전자 특정 (황산 변환 후 TG 대 CG) 대상 어 닐 링 및 뇌관 확장을 위한 특정 cpgs 내 겠에 대 한 설계를 사용 합니다. 그것의 디자인 뿐만 아니라 사용할 수 있는 인간 게놈 시퀀스, 하지만 실험적으로 확인 된 규제 영역 문의 인코딩 및 합8등의 여러 줄에서 인수를 반영 한다. 인간 methylomic 조사에 그것의 광범위 한 사용에도 불구 하 고 비슷한 플랫폼 모델 동물에 대 한 존재 하지 않습니다. 또한, 배열 기반 형식 프로브 배치에 사용할 수 있는 영역에 중요 한 제약 조건을 배치합니다. 지난 몇 년 동안, 노력 캡처 프로브 설계 및 차세대 시퀀싱의 높은 처리량 기능 제공 대상-특이성 결합 되었습니다. 이러한 노력 결과 마우스 게놈 (마우스 메 틸-Seq), 메 틸 화9,10에 두뇌 또는 glucocorticoid 유발 차이 식별 하는 데 사용 했다의 시퀀싱 기반 대상 농축 시스템입니다. 다른 모델과 아닌 모델 동물에 대 한 비슷한 플랫폼이이 동물에 epigenomic 연구를 촉진 하기 위하여 필요 합니다.
여기, 우리는 쥐에 methylomic 분석을 소설이 플랫폼의 구현 방법을 보여 줍니다. 쥐는 약리학, 신진 대사, 내분비학, 동작에 중요 한 동물 모델을 역임 했다. 예를 들어 약물 독성, 비만, 스트레스 반응, 또는 마약 중독을 기본 메커니즘을 이해 하는 증가 필요가 있다. 이러한 조건과 관련 된 methylomic 변경 캡처 수 있는 높은 처리량 플랫폼 메커니즘에 대 한 우리의 이해를 증가할 것입니다. 쥐 게놈은 아직 규제 지역에 대 한 주석 부족, 이후 비 중복 발기인, CpG 섬, 섬 해안11, 통합 하 고 이전 쥐 메 틸-Seq 플랫폼12GC 풍부한 시퀀스를 식별 합니다.
쥐 게놈에 대 한 SureSelect 대상 농축 (일반적으로 메 틸-Seq 라고도 함) 플랫폼의 성공적인 설계 및 구현 평가, 우리는 만성 변수 스트레스 (CVS)13 차동 methylated 식별 하의 쥐 모델 고용 강세와 스트레스 동물 사이의 영역입니다. 우리의 플랫폼, 프로토콜, 디자인과 구현 체 누구의 게놈 시퀀스는 이미 사용할 수 있지만 제대로 주석이 남아에 포괄적이 고 편견 후 조사를 실시 하 고 싶은 수 있습니다 수사에 대 한 유용할 수 있습니다.
모든 실험은 따라와 모든 관련 규제 및 제도 지침 준수, 기관 동물 관리 및 사용 위원회 의학의 존스 홉킨스 학교에서를 포함 하 여 완료 되었습니다.
1입니다. 동물
2. 만성 변수 스트레스
3. 내 분 비 분석 실험
4입니다. 동작
5. 쥐 메 틸-Seq의 디자인
6. 쥐 메 틸-Seq 라이브러리 Genomic DNA에서의 건설
참고: 일괄 처리 효과 제거 하려면, 동시에 여러 개의 샘플을 처리 하 고 따라 마스터 믹스 최대 규모. 상업적으로 이용 가능한 DNA 추출 키트를 사용 하 여 DNA를 추출 합니다. 열 또는 강 수 기반 방법의 얻을 높은 품질 genomic DNA (260/280 비율 ~ 1.8). 페 놀 기반 방법의 사용을 권장 하지 않습니다. Elute 또는 낮은 테 버퍼 (10 m m 테, 0.1 m m EDTA, pH 8.0)에서 DNA를 resuspend.
7. 다음-세대 Sequencer에서 시퀀싱
8. DMRs를 식별 하는 분석
9입니다. 황산 Pyrosequencing에 의해 유효성 검사
쥐 메 틸-Seq 플랫폼을 성공적으로 구현 하는 몇 가지 기준에 따라 달라 집니다. 그림 1 연구의 전체 워크플로 보여 줍니다 및 앞으로 이동 하기 전에 필요한 특정 품질 관리 (QC) 단계를 하이라이트. 고려해 야 할 첫 번째 요소 중 하나는 methylome에 걸쳐 발생 하는 epigenetic 변화의 크기를 결정 하는 동물 모델 및 스트레스 처방의 견고성입니다. 이후 우리의 동물 작품 corticosterone (CORT) 노출 변화 DNA 메 틸 화19,20에 이어질 수 있는 우리의 이전의 관찰에 입각 한, 우리의 만성 변수 스트레스 (CVS) 처방 필요한 생산 충분 한 사후 경직 높은 플라즈마 CORT 수준으로 쥐를 강조 했다. 전형적인 주간 CVS 처방은 표 1 에 나와 아침, 오후, 매일 스트레스의 구성 그리고 하룻밤 habituation를 방지 하기 위해 지속적으로 변경 되 고 스트레스 반응 감소. 3 주 처방을 통해 스트레스 동물 의미 플라즈마 CORT의 상당히 높은 수준의 전시 [일 4-21, 제어: 32.7 3.7 ng/mL, 스트레스: 103.0 11.9 ng/mL (말은 sem의), P = 10-4, 그림 2Ax 2.2]의 강세, 그 이상의 동물을 제어 합니다. 일관 되 게,이 동물 또한 보였다는 상승에 더 큰 불안 같은 행동 플러스 미로 (EPM), EPM 및 더 적은 시간 (그림 2B) 두 팔을 벌려에서 닫힌된 팔에는 훨씬 더 많은 시간으로 표시. 이러한 결과 CVS 노출 중요 일 분 비 및 행동 변화, 이러한 변화 했다 특정 DNA 메 틸 화 서명 관련 여부 조사 선도 이끌어 보여줍니다.
우리는 메 틸-Seq 라이브러리의 성공적인 건설에 대 한 중요 한 여러 검사점을 강조. 충분 한 양의 DNA로 시작 필요, 쥡니다, 여러 세척/정화, 대상 농축 이며 황산 변환 단계 연속적으로 완성 된 라이브러리에서 DNA의 양을 감소. 비록 몇 가지 PCR 증폭 단계 완화 DNA 템플렛의 손실, 과도 한 PCR 주기 숫자 높은 중복 읽기 발생할 수 있습니다. 현재 쥐 메 틸-Seq 연구를 위해 쥐 당 혈액 gDNA의 2 세대 사용 되었다. 우리 주의 메 틸-Seq 라이브러리 시작 DNA 금액 500으로 만들어질 수 있다 니. 작은 시작 물자 후속 시퀀싱에 대 한 라이브러리의 부족 한 금액을 생산의 위험이 증가 하지만 FACS (형광 활성화 된 세포 분류)에 의해 고립 된 DNA에서 라이브러리를 생성 하거나 바늘 펀치, 사용자 수 있습니다. QC는 DNA 분자 무게, 양, 및 몰 제공 하는 bioanalyzer에 샘플의 1 L의 전기 이동 법에 의해 수행 됩니다. bioanalyzer의 사용을 필요로 하는 3 개의 중요 한 단계는: 1) (~ 170 혈압, 레드, 그림 3); DNA의 충분 한 전단 되도록 쥡니다 단계를 다음 (~ 200 bp, 블루, 그림 3) PCR;에 의해 그들의 연속적인 증폭 되도록 전단된 DNA의 평균 크기에 변화에 의해 표시 2) 다음 어댑터 결 찰 단계 그리고 3) 수량 및 시퀀싱에 대 한 라이브러리의 크기를 되도록 최종 라이브러리 정화 단계를 다음과 같은.
R 패키지 BSSeq와 Bioconductor에 BSmooth 시퀀싱 데이터18황산을 분석 하는 데 사용 되었다. 그들은 도구와 정렬 시퀀스 읽습니다, 품질 관리를 수행 하기 위한 메서드를 포함 하 고 지역 (DMRs) 갠 차동 식별. BSmooth 소프트웨어는 황산 변환 게놈 시퀀스를 원시 입력된 읽기의 정렬에 의해 CpG 수준 측정 요약을 얻기 위해 내부 시퀀스 aligner로 Bowtie 2.016,17 를 호출 합니다. 정렬 된 읽기 다음 체계적인 시퀀싱 및 다운스트림 분석을 왜곡 수 있습니다 기본 호출 오류를 식별 하는 엄격한 품질 관리 절차를 통해 필터링 됩니다. 플롯의 시리즈는 시각적으로 필터링이 과정에 도움 생성 됩니다. 시퀀싱 통계와 같은 정렬된 읽기, % 대상의 포장 소비재 범위, 다른 (표 2) 중 당 문서 관련 정보에도 생성 됩니다. 일단 데이터 필터링, 스무 딩/정규화 알고리즘 수행, 모든 소비재 모든 QC에 따라 예상된 메 틸 화 값 할당 각 샘플에서 읽습니다 메 틸 화의 더 정확한 호출 되도록 cpgs 내 겠 이웃에서 추정 시퀀스 범위 낮은 경우에도 상태. 이 값에는 메 틸 화 CpG 사이트 각에서 확률의 매끄러운된 견적을 제공합니다. 두 치료 그룹 및 순위 이상에 가장 크게 다른 게놈 영역 사이 각 샘플의 부드러운된 메 틸 견적의 평균을 비교 하 여 DMRs의 목록이 생성 됩니다 (표 3).
스트레스와 강세 그룹 사이 DMR와 Rt1 m4, 쥐 중요 한 조직 적합성 유전자의 발기인에 있는 최고 강세 동물 (그림 4A) 보다 모든 cpgs 내 겠 전체 메 틸 화 상부를 전시 하는 동물을 강조 했다. 메 틸-Seq 플랫폼 및 데이터 분석의 성공적인 구현을 확인, 뇌관 DMR, 및 압박 하 고 강세가 없는 동물의 (메 틸-Seq에 의해 시퀀스 8 하지 시퀀스 8) 전체 코 호트에서 혈액 DNA 메 틸 화 수준에 대 한 설계 되었습니다. 황산 pyrosequencing에 의해 평가 됐다. 결과 분석 12 cpgs 내 겠 밖으로 10에 걸쳐 DNA 메 틸 화에 상당한 증가 입증 (5.1-10.4% 메 틸 화, P에에서 변화 < 0.037, 그림 4B). KEGG 통로 분석은 모든 스트레스와 관련 된 경로 식별 하기 위해 명목상 상당한 DMRs에서 수행 되었다. 일관 되 게, DMR 관련 경로 질병 당뇨병, 심혈 관 질환 및 암 (표 4)와 같은 만성 스트레스 노출과 관련 된 연루. 21 , 22 , 23 스트레스 후 데이터와 노출의 정도 사이 연결을 보여주는, 메 틸 화 CpG-10 레벨 각 동물에 대 한 평균 3 주 CORT 수준 비교 되었다. 결과 내 분 비 및 메 틸 화 데이터 사이 겸손 한 상관 관계를 보였다 (R2= 0.54, P = 0.001, 그림 5).
그림 1: 쥐 메 틸-Seq 플랫폼에 대 한 전반적인 회로도 워크플로. 1 g의 혈액에서 추출한 게놈 DNA의 스트레스 및 제어 쥐 먼저 시퀀싱, 분석 및 대상 식별에 대 한 메 틸-Seq 라이브러리 구성에 대 한 처리. 또 다른 100 DNA의 ng 황산 pyrosequencing에 의해 확인 된 후 대상의 독립적인 유효성 검사에 사용 됩니다. 이 그림의 더 큰 버전을 보려면 여기를 클릭 하십시오.
그림 2: 만성 변수 스트레스 (CVS)에 노출 쥐에 내 분 비 및 행동 변화를 리드. (A) 여러 samplings corticosterone (CORT)의 설명의 견고성 CVS 처방 3 주. 혈액 샘플은 매일 스트레스 처방 전에 아침에 수집 되었다. (B) 스트레스 동물 미로 (EPM) 플러스 시간에는 상승의 두 팔을 벌려 덜 닫힌된 무기에 더 많은 시간을 보냈다. Boxplots 각 동물에 대 한 데이터 요소와 함께 표시 됩니다. 스튜던트 T-검정은 통계적 의미에 대 한 수행 했습니다. * P < 0.05, * * P < 0.01, 그리고 * * * P < 0.001. 이 그림의 더 큰 버전을 보려면 여기를 클릭 하십시오.
그림 3:는 bioanalyzer에 전단 및 어댑터 출혈 쥐 DNA의 정량. 빨간색과 파란색 곡선 수량과 각각 등온 sonicator 및 어댑터 결 찰, 전단 뒤 genomic DNA (레드)의 크기를 표시 합니다. 각 줄 한 샘플 및 빨간색 나타내고 파란색 곡선 여러 단계 (최종 수리, 3'-adenylation, 및 샘플 정리) DNA의 두 손실 반영 및 bp 크기는 어댑터의 결 찰 때문에 증가. 날카로운 봉우리 25 bp와 1500 bp는 표준 마커를 로딩 버퍼에 추가 되었습니다. 이 그림의 더 큰 버전을 보려면 여기를 클릭 하십시오.
그림 4: CVS 유도 후 성적인 변화는 쥐 메 틸이 검색 쥐 메 틸-Seq 데이터의 (A) 분석 스트레스 (빨강) 및 컨트롤 (블루) 쥐 사이 차동 methylated 지역 (DMR)로 서 Rt1m4 유전자의 발기인 연루. Rt1m4 DMR (분홍색 음영된 지역)에 대 한 그래픽 출력 각 그룹 (빨간색 이나 파란색 라인), 그리고 각 동물 (빨간색 또는 파란색 점)에 대 한 % 메 틸 화 수준에 각 포장 소비재 (회색 세로선), 4 개의 샘플을 표시합니다. (B) 12 cpgs 내 겠은 DMR 내 황산 pyrosequencing에 의해 확인 했다. 막대 그래프는 SEM, 의미 표시 되 고 통계적 의미에 대 한 스튜던트 T-검정 수행 했다. * P < 0.05. 이 그림의 더 큰 버전을 보려면 여기를 클릭 하십시오.
그림 5: 선형 회귀 분석 보여 겸손 상관관계 % DNA CpG-10 Rt1m4 의 메 틸 화와 3 주 의미 플라즈마 CORT 수준을 모두 스트레스의 제어 동물 (N = 16). 데이터 스트레스 동물에서 빨간색 원으로 표시 됩니다. 이 그림의 더 큰 버전을 보려면 여기를 클릭 하십시오.
주 | 주 1 | 주 2 | 3 일 | 4 일 | 5 일째 | 주 6 | 7 일 |
오전 | 구속 | 수영 | 콜드 룸 | 수영 | 구속 | 쉐이 커 | 수영 |
오후 | 쉐이 커 | 케이지 기울기 | 구속 | 쉐이 커 | 콜드 룸 | 구속 | 콜드 룸 |
하룻밤 | 음식 제한 | 젖은 침구 | 절연 | 빛에 | 크롤 링 | 빛에 | 젖은 침구 |
표 1: 만성 변수 스트레스 처방 (CVS)의 전형적인 주간 일정.
시퀀싱 통계 | 스트레스1 | 제어1 |
(n = 4) | (n = 4) | |
끝 읽기 (당) 한 쌍 | 89,290,397 | 80,165,674 |
고유 하 게 매핑된 짝된 끝 (UMPER)를 읽으십시오 | 39,200,255 | 35,013,406 |
정렬 속도/매핑 효율 (UMPER / 1) | 44% | 44% |
중복 읽기 (UMPER %) | 73% | 65% |
Deduplicated UMPER | 10,481,031 | 12,306,018 |
평균 읽기 깊이 범위 (x) (ARDC) | 6 x | 6 x |
Cpgs 내 겠 (N) | 12,056,878 | 12,056,878 |
Cpgs 내 겠 ARDC (x) | 2 x | 2 x |
10 읽기 (N)으로 cpgs 내 겠 | 481,383 | 595,850 |
10 읽기로 cpgs 내 겠 ARDC (X) | 19 | 19 |
대상 cpgs 내 겠 (완전 한 조사 대상 지역 중복에) | 1,923,872 | 2,007,638 |
대상의 cpgs 내 겠 ARDC (x) | 7 x | 8 배속 |
10 읽기 (N)와 대상 cpgs 내 겠에 | 428,249 | 531,419 |
Cpgs 내 겠 적어도 10와 대상 ARDC (x)에 읽습니다. | 18 x | 18 x |
대상 (대상 지역 조사와 1 이상의 기본적인 쌍 중복 당) (UMPER) | 8,277,715 | 9,369,523 |
% (Deduplicated UMPER)의 대상에 | 78% | 77% |
대상 (총 기지 매핑)에 Mb | 125 mb | 128 mb |
에 대상 평균 읽기 깊이 범위 (x) (ARDC) | 9 x | 10 배 |
1 각 그룹에 과목 전체 평균에 따라 시퀀싱 통계 |
표 2: 시퀀싱 쥐 메 틸-Seq 플랫폼에서 얻은 통계.
chr | 시작 | 끝 | 유전자 | 거리 | areaStat | meanDiff | 스트레스 | 제어 | 방향 |
chr20 | 1,644,246 | 1,644,390 | RT1-M4 | in_gene | 93.03 | 0.22 | 0.33 | 0.11 | 이득 |
chr5 | 160,361,352 | 160,361,564 | LOC690911 | in_gene | -70.75 | -0.19 | 0.72 | 0.91 | 손실 |
chr3 | 61,138,281 | 61,138,330 | RGD1564319 | 265569 | 61.79 | 0.21 | 0.94 | 0.72 | 이득 |
chr2 | 143,064,811 | 143,065,010 | Ufm1 | 8569 | -59.48 | -0.11 | 0.13 | 0.24 | 손실 |
chr7 | 30,764,111 | 30,764,284 | Ntn4 | in_gene | 57.04 | 0.21 | 0.94 | 0.73 | 이득 |
chr17 | 12,469,112 | 12,469,218 | Idnk | 41996 | -50.91 | -0.13 | 0.74 | 0.88 | 손실 |
chr7 | 47,101,725 | 47,101,930 | Pawr | in_gene | -50.54 | -0.12 | 0.64 | 0.76 | 손실 |
chr5 | 76,111,248 | 76,111,822 | Txndc8 | 151703 | -50.38 | -0.11 | 0.85 | 0.96 | 손실 |
chr11 | 80,640,132 | 80,640,356 | Dgkg | in_gene | -50.07 | -0.16 | 0.73 | 0.89 | 손실 |
chr8 | 71,759,248 | 71,759,411 | Mir190 | 210226 | -47.84 | -0.17 | 0.58 | 0.75 | 손실 |
표 3: 최고 10 차동 갠 지역. 각 DMR에 대 한 출력 테이블 오른쪽 열에 왼쪽에서 보여줍니다: 의미, 염색체 위치 (chr), 좌표 (시작/종료), 유전자 이름, 거리는 전사에서 사이트를 시작, 차동 지역 통계 사이 강조 하 고 제어 그룹 (areaStat) 차동 메 틸 화 (meanDiff)에 대 한 각 DMR에서 말은 메 틸 화 수준의 스트레스와 제어 그룹 (스트레스/제어), 및 메 틸 화의 방향 컨트롤에서 변경.
KEGG 통로 용어 | 유전자 수 | % | P-값 | Benjamini |
당뇨병 | ||||
타입 II 당뇨병 mellitus | 12 | 0.1 | 3.6 x 10-4 | 9.8 x 10-3 |
심혈 관 질환 | ||||
혈관 평활 근 수축 | 18 | 0.1 | 1.6 x 10-3 | 3.6 x 10-2 |
Arrhythmogenic 오른쪽 심 실 심근 (ARVC) | 13 | 0.1 | 10-3 x 4.0 | 7.1 x 10-2 |
동 공이 확장 되어 심장 근육 병 증 | 14 | 0.1 | 7.6 x 10-3 | 1.2 x 10-1 |
신경 기능 | ||||
장기 potentiation | 11 | 0.1 | 1.5 x 10-2 | 1.4 x 10-1 |
신호 | ||||
MAPK 신호 통로 | 35 | 0.2 | 2.4 x 10-4 | 9.9 x 10-3 |
칼슘 신호 전달 경로 | 22 | 0.1 | 1.2 x 10-2 | 1.4 x 10-1 |
Chemokine 신호 통로 | 21 | 0.1 | 1.2 x 10-2 | 1.3 x 10-1 |
암 | ||||
암에 경로 | 42 | 0.3 | 4.1 x 10-5 | 3.4 x 10-3 |
Glioma | 15 | 0.1 | 4.4 x 10-5 | 2.4 x 10-3 |
비 작은 세포 폐 암 | 10 | 0.1 | 7.9 x 10-3 | 10-1 x 1.1 |
대 장 암 | 13 | 0.1 | 8.4 x 10-3 | 10-1 x 1.1 |
만성 골수성 백혈병 | 12 | 0.1 | 1.2 x 10-2 | 1.3 x 10-1 |
표 4: KEGG 통로 분석 DMRs 메 틸이 쥐에서 확인
이 연구에서 우리는 설계 및 쥐 게놈에 대 한 메 틸-Seq 플랫폼 구현. 스트레스의 쥐 모델의 유틸리티를 시연 함으로써 우리는 실험 및 분석 파이프라인 두 비교 그룹 사이 차동 갠 영역을 제공할 수 있습니다 설명 했다.
플랫폼의 성공적인 구현을 보장 하기 위해, 몇 가지 중요 한 단계를 관찰 해야 합니다. 첫 번째, 초기 DNA 품질 및 수량 최종 메 틸-Seq 라이브러리의 양과 품질에 큰 영향을 있다. 우리는 fluorometer 보다는 분 광 광도 계, 우리의 DNA 측정 현재 이중 가닥 DNA의 수량을 반영 되도록 사용. 분자 크기와 어댑터 결 찰 후 전단에 따라 DNA의 수량을 측정 하는 bioanalyzer 사용 되었다. 분자 크기 다음이 단계 사이 "변화"를 확인 하는 것은 후속 단계에서 PCR 어댑터 중재를 받게 될 것 이다 각 DNA 파편의 끝에 어댑터의 존재를 확인 중요 합니다. DNA는 어댑터 결 찰 단계의 끝에 남은 수량 또한 중요 하다, 적어도 100부터 ng 라이브러리 제품의 충분 한 수량 대상 농축 황산 변환 단계 후 사용할 수 있도록 하려면이 단계에서 필요 하다. 최종 높은 감도 측정 라이브러리는 제대로 다음-세대 시퀀서에 후속 클러스터링에 대 한 희석 수 있도록 생성 된 메 틸-Seq 라이브러리에서 수행 되었다. 마지막으로, 황산 pyrosequencing 높은 양적, 독립적인 방법으로 분석 파이프라인의 정확도 평가 하기 위해 고용 되었다. 원래 샘플 및 추가적인 동물을 사용 하 여 복제를 사용 하 여 최종 검증은 실험 DNA 메 틸 화에 생물학적으로 상당한 변화를 감지할 수 있도록 중요 한 단계입니다.
우리는 또한 프로토콜에서 편차 또는 경우 문제 발생 시 몇 가지 지침을 포함 합니다. 첫째, 그것은 최종 수리, 어댑터 결 찰, 또는 자석 구슬 정화 단계 동안 너무 많은 DNA를 잃을 수입니다. 시작 하는 양의 DNA 작은 수 또는 (< 200 ng) 제한 된 조직/DNA 가용성 또는 형광 활성화 된 세포 분류 등 다양 한 농축 방법의 구현. DNA의 과도 한 손실에 대 한 보상을 수 또는 도서관 건설 프로토콜에 걸쳐 DNA 양을 시작 낮은 사이클 수를 증가 하는 두 개의 라이브러리 증폭 단계 수 있습니다. 그러나, 더 이상 추가 2-3 주기는 권장, 과도 한 템플릿 증폭 시퀀싱 되 고 중복 읽기의 증가를 끌 것입니다. 이러한 중복 백분율 메 틸 화 계산에 바이어스를 방지 하기 위해 맞춤 단계에서 제외 됩니다. 둘째, 평균 DNA 크기 30 bps 이상으로 증가 하지 않습니다, 경우는 시 약은 새로운 T4 DNA 중 합 효소, Klenow, 또는 T4로 리가 나이가 있을 수 있습니다 있는지 확인 합니다. 상업적으로 이용 가능한 대체 시 약을 사용할 수 있습니다.
또한, 예측된 DMRs pyrosequencing, DNA 메 틸 화 차이 존재 하지 않는 또는 분석에 의해 예측 보다 훨씬 적은 있다 의해 확인 하지 수도 가능 하다. 후보 지역 가난한 정품 확인은 너무 때 pyrosequencing 결과 차동 메 틸 화를 확인 하지 않습니다 같은 많은 게놈 넓은 분석에 대 한 일반적인 문제 또는 효과 크기는 분석에 의해 예측 보다 훨씬 작습니다. BSmooth 여러 cpgs 내 겠 창 한 분석 패키지는 "부드럽게" 메 틸 화 수준입니다. 현재 실험에 대 한 BSmooth 연루는 DMR의 메 틸 화 수준의 황산 pyrosequencing에 의해 확인 되었다. 그러나, 거기 있을 것 이다 메 틸 화 수준 예측 BSmooth와 pyrosequencing에 의해 확인 사이의 불일치. 모든 연속 DNA 메 틸 화에 50% 이상 차이가 날 수 있습니다 cpgs 내 겠 또는 메 틸 화 값으로 인해 제외 된 cpgs 내 겠은 DMR 내 cpgs 내 겠 평균 메 틸 화 값을 추정 다듬기 기능에서 발생 하는 불일치 하위 임계값 깊이 읽기. R-패키지 MethylKit24 cpgs 내 겠 또는 심지어 단일 cpgs 내 겠 그 메 틸 화 수준 연관 강하게 pyrosequencing의 유효성을 검사 하는 그의 작은 윈도우를 식별 하기 위해 사용할 수 있습니다. 다른 패키지를 구현 하 고 pyrosequencing에 의해 그들의 예측된 영역 또는 차동 메 틸 화의 cpgs 내 겠 테스트 데이터의 안정성을 보장 합니다. 또는, 원래 메 틸-Seq 라이브러리 resequenced 하 고 읽기 깊이 증가 읽기 파일에 추가 될 수 있습니다. 때문에 메 틸 화 수준의 결정 반 양적 및 구술 읽기 수로 [(# of CpGs) / (TpGs + cpgs 내 겠 #)], 주어진된 포장 소비재의 백분율 메 틸 화 값의 정확도 높일 것에 대 한 읽기 깊이 증가. 이 연구에서 우리 메 틸 화 값 적어도 10 읽기에 의해 결정 되었다 고 각 포장 소비재에 대 한 19 x의 전체 읽기 범위 달성 cpgs 내 겠 고려.
쥐 메 틸-Seq 플랫폼의 제한 없이 아니다. 전체 게놈 황산 시퀀싱 보다 비용 효과적인 동안, 그것은 상당히 다른 방법 보다 더 비싸다입니다. 그럼에도 불구 하 고, 비용의 대부분을 캡처 시스템 그리고 시퀀서에 차선을 구매 했다. 크로스-조직 비교 덜 큰 (25-70%) 차이12 DNA 메 틸 화에 따른 요구와 필요한 읽기 깊이 따라 레인 당 더 많은 샘플을 멀티플렉싱 높은 용량 플랫폼을 사용 하 여 비용을 줄일 수 있습니다. 또한, 샘플 준비는 다른 방법 보다 더 많은 시간이 걸리는. 차세대 시퀀싱을 통합 하는 다른 풀 다운 방식와 마찬가지로, 하는 동안 추가 황산 변환 및 정화 단계 작업 부하를 추가 합니다. 전반적으로, 메 틸-Seq 플랫폼 전체 게놈 시퀀싱에 대 한 비용 효율적인 대안 이며 상당히 microarray 기반 플랫폼 분석 보다 더 이상의 2.3 백만 cpgs 내 겠에서 기본적인 쌍 해상도 제공. 날짜 하려면, 상용 인간과 마우스 메 틸-Seq 플랫폼 알코올 의존 변화 원숭이 뇌25,26, 마우스 뇌9, 및 혈액 neurodevelopmental 유전자를 사용 되었습니다. 스테로이드 제제10대상 또한, 제한 효소에 의해 시퀀스 인식에 특정 영역을 대상으로 기능 그것은 종 간 비교에 대 한 이상적인 플랫폼을 통해 있습니다. 이 연구를 위해 우리는 많은 약리학, 대사, 행동 실험 게놈 넓은 methylomic 도구의 이점 없이 수행 되는 쥐에 대 한 메 틸-Seq 플랫폼을 설계 되었습니다. 우리의 데이터는 스트레스의 쥐 모델에서 DMRs를 검색 하기 위해 사용할 수 있습니다을 전체 플라즈마 CORT 수준 같은 다른 생리 적인 매개 변수 상관을 보여준다.
메 틸-Seq 플랫폼 규제 영역을 문서화 하는 충분 한 실험적 증거를 없을 수 있습니다 시퀀스 유전자와 동물에 있는 후 성적인 실험에 이상적입니다. 이러한 영역은 사용할 수 있습니다 때 추가 지역 맞춤형 고 현재 버전에 연결 된 수 있습니다. 또한, 플랫폼 대상 농축 제한 효소 인식에 의해 제약 되지 않습니다 때문에, 비교 유전체학에 이상적입니다. 예를 들어, 그것은 특정 제한 사이트를 비호 하는 여부에 관심사의 어떤 유전자의 발기인 지구를 캡처할 수 있습니다. 마찬가지로, 마우스 또는 관심의 게놈에 보존은 인간에서 확인 된 그들과 같은 어떤 규제 영역을 캡처할 수 있다.
원고는 Agilent 기술에서 경연 대회 상금의 일부입니다.
이 연구는 NIH 그랜트 MH101392 (RSL) 및 다음 수상 및 재단에서 지원에 의해 투자 되었다: NARSAD 젊은 수 사관 상, 마가렛 앤 가격 조사 기금, 찰스 T. 바 우 어 재단, 베이커를 통해 제임스 Wah 기분 장애 학자 기금 재단, 그리고 프로젝트 일치 재단 (RSL).
Name | Company | Catalog Number | Comments |
Radioimmuno assay (RIA) | MP Biomedicals | 7120126 | Corticosterone, 125I labeled |
Master Pure DNA Purification Kit | Epicentre/Illumina | MC85200 | |
Thermal-LOK 2-Position Dry Heat Bath | USA Scientific | 2510-1102 | Used with 1.5 mL tubes |
Vortex Genie 2 | Fisher | 12-812 | Vortex Mixer |
Ethyl alcohol, Pure | Sigma-Aldrich | E7023 | 100% Ethanol, molecular grade |
Centrifuge 5424 R | Eppendorf | - | Must be capable of 20,000 x g |
Qubit 2.0 | ThermoFisher Scientific | Q32866 | Fluorometer |
Qubit dsDNA BR Assay Kit | ThermoFisher Scientific | Q32850 | |
Qubit dsDNA HS Assay Kit | ThermoFisher Scientific | Q32851 | High sensitivity DNA detection reagents |
Qubit Assay Tubes | ThermoFisher Scientific | Q32856 | |
SureSelectXT Rat Methyl-Seq Reagent Kit | Agilent Technologies | G9651A | Reagents for preparing the Methyl-Seq library |
SureSelect Rat Methyl-Seq Capture Library | Agilent Technologies | 931143 | RNA baits for enrichment of rat targets |
IDTE, pH 8.0 | IDT DNA | 11-05-01-09 | 10 mM TE, 0.1 mM EDTA |
DNA LoBind Tube 1.5 mL | Eppendorf | 22431021 | |
Covaris E-series or S-series | Covaris | - | Isothermal sonicator |
microTUBE AFA Fiber Pre-Slit Snap-Cap 6x16mm (25) | Covaris | 520045 | |
Water, Ultra Pure (Molecular Biology Grade) | Quality Biological | 351-029-721 | |
Veriti 96 Well-Thermal Cycler | Applied Biosystems | 4375786 | |
AMPure XP Beads | Beckman Coulter | A63880 | DNA-Binding magnetic beads |
96S Super Magnet | ALPAQUA | A001322 | Magnetic plate for purification steps |
2200 TapeStation | Agilent Technologies | G2965AA | Electrophresis-based bioanalyzer |
D1000 ScreenTape | Agilent Technologies | 5067-5582 | |
D1000 ScreenTape High Sensitivity | Agilent Technologies | 5067-5584 | |
D1000 Reagents | Agilent Technologies | 5067-5583 | |
D1000 Reagents High Sensitivity | Agilent Technologies | 5067-5585 | |
DNA110 SpeedVac | ThermoFisher Scientific | - | Vacuum Concentrator |
Dynabeads MyOne Streptavidin T1 magnetic beads | Invitrogen | 65601 | Streptavidin magnetic beads |
Labquake Tube Rotator | ThermoFisher Scientific | 415110Q | Nutator Mixer is also acceptable |
EZ DNA Methylation-Gold Kit | Zymo Research | D5006 | Bisulfite conversion kit. Contains Binding, Wash, Desulphonation, and Elution buffers |
Illumina Hi-Seq 2500 | Illumina | - | Next-generation sequencing machine |
PCR and Pyrosequencing Primers | IDT DNA | Variable | |
Taq DNA Polymerase with ThermoPol Buffer - 2,000 units | New England BioLabs | M0267L | |
Deoxynucleotide (dNTP) Solution Set | New England BioLabs | N0446S | |
Pyromark MD96 | QIAGEN | - | Pyrosequencing machine |
Ethyl Alcohol 200 Proof | Pharmco-Aaper | 111000200 | 70% Ethanol solution |
Sodium Hydroxide Pellets | Sigma-Aldrich | 221465 | 0.2 M NaOH denature buffer solution |
Tris (Base) from J.T. Baker | Fisher Scientific | 02-004-508 | 10 mM Tris Acetate Buffer wash buffer solution |
PyroMark Gold Q96 Reagents (50x96) | QIAGEN | 972807 | Reagents required for pyrosequencing |
PyroMark Annealing Buffer | QIAGEN | 979009 | |
PyroMark Binding Buffer (200 mL) | QIAGEN | 979006 | |
Streptavidin Sepharose High Performance Beads | GE Healthcare | 17-5113-01 | Streptavidin-coated sepharose beads |
PyroMark Q96 HS Plate | QIAGEN | 979101 | Pyrosequencing assay plate |
Eppendorf Thermomixer R | Fisher Scientific | 05-400-205 | Plate mixer. 96-well block sold separately (cat. No 05-400-207) |
SureDesign Website | Agilent Technologies | - | Target capture design software (https://earray.chem.agilent.com/suredesign/) |
UCSC Genome Browser | University of California Santa Cruz | - | rat Nov 2004 rn4 assembly |
Agilent Methyl-Seq Protocol | Agilent Technologies | - | https://www.agilent.com/cs/library/usermanuals/public/G7530-90002.pdf |
JoVE'article의 텍스트 или 그림을 다시 사용하시려면 허가 살펴보기
허가 살펴보기This article has been published
Video Coming Soon
Copyright © 2025 MyJoVE Corporation. 판권 소유