본 연구는 염색체 서열이 이미 알려져 있거나 관심 없는 플라스미드의 나노포어 시퀀싱 워크플로우를 설명합니다. 초기 플라스미드 DNA 풍부화 단계는 최적의 플라스미드 서열 깊이를 보장하고 플로우 셀당 시퀀싱 가능한 샘플 수를 극대화합니다.
방법 논문
* These authors contributed equally
본 연구는 염색체 서열이 이미 알려져 있거나 관심 없는 플라스미드의 나노포어 시퀀싱 워크플로우를 설명합니다. 초기 플라스미드 DNA 풍부화 단계는 최적의 플라스미드 서열 깊이를 보장하고 플로우 셀당 시퀀싱 가능한 샘플 수를 극대화합니다.
본 연구는 플라스미드 DNA 수율과 시퀀싱 정확도를 극대화하기 위해 설계된 나노포어 시퀀싱 워크플로우를 설명합니다. 연구자들은 플라스미드 DNA 추출, 빠른 바코딩 및 어댑터 연결을 통한 라이브러리 준비, 라이브러리를 플로우 셀에 적재하여 시퀀싱, 염기 호출(실시간 또는 사후 호출), Autocycler를 이용한 서열 조립 등 공정의 각 단계를 다룹니다. 저자들은 또한 컴퓨터 요구사항을 명시하는데, 이는 주로 베이스콜링이 실시간으로 이루어지는지에 따라 달라집니다. 저자들은 이 워크플로우를 사용해 10가지 임상 균주에서 정제한 플라스미드를 시퀀싱했습니다. 숙주-세포 균주 차이를 변수로 제거하기 위해, 저자들은 이 임상적 균주들을 공통 수용자 균주와 결합시키고 환합 균주를 시퀀싱했습니다. 이 실험들은 겔 이동성으로 추정된 플라스미드 크기와 일치하는 플라스미드 조립체를 생성했으며, 3분의 2의 경우에는 단독 연마 조립체와도 일치하거나 0.05% 쌍별 서열 동일성 내에 있는 경우도 있었습니다. 플라스미드 DNA 서열의 풍부화는 플라스미드 시퀀싱의 효율성을 극대화하여 최대 24개의 샘플을 최적의 커버리지로 병렬 시퀀싱할 수 있게 합니다. 이 프로토콜은 4에서 174 kb까지의 플라스미드에 적합하며, 정확도에 영향을 주지 않고 최대 72%의 염색체 오염을 견딜 수 있습니다. 이 프로토콜은 염색체 서열이 이미 알려져 있거나 관심 없는 상황, 예를 들어 플라스미드 서열 검증, WGS로 이미 얻은 플라스미드 서열의 정확도 향상, 또는 접합으로 포획된 플라스미드 시퀀싱 등에 이상적입니다. 이 연구를 설명하는 데 사용된 유기체는 대 장균으로, 재조합 유전자 발현과 접합을 통한 플라스미드 포획에 사용되는 숙주를 대표합니다.
플라스미드는 박테리아에서 흔히 발견되는 자가 복제 DNA 조각입니다. 크기는 매우 다양하여 <800 bp에서 2.5 Mbp 사이이며, 일반적으로 이중분포를 따르며, 다중복제 비접합 플라스미드에 해당하는 작은 봉크와 저접합 접합 플라스미드에 해당하는 큰 봉크가 있습니다. 플라스미드는 미생물 생태계에서 중요한 역할을 하며, 숙주가 환경 문제에 적응하거나 새로운 생태적 지위를 확보하는 데 도움을 주는 유전자를 가지고 있습니다. 플라스미드는 또한 가속화된 진화의 플랫폼 역할을 하며, 유전적 가소성을 높이고 균주, 종, 속, 심지어 가족 간에 걸쳐 수평 유전자 전달을 매개합니다 3,4. 네트워크 및 비교 유전체 분석을 통해 환경 저항체가 임상적으로 중요한 항생제 내성 유전자(ARG)의 지속적인 저장소로서 역할을 하여, 생태계 환경 전반에 걸쳐 유지되고 확산되는 데 기여함을 확인했습니다5. 플라스미드와 그에 포함된 항생제 내성 및 병원성 유전자를 추적하는 것은 전파 경로를 식별하는 데 도움이 됩니다. 6. 따라서 플라스미드의 식별과 특성 분석은 다양한 공중보건, 임상 미생물학, 생명공학 응용 분야에서 매우 중요한 관심사입니다.
미지의 DNA 주형을 고처리량으로 시퀀싱할 수 있는 차세대 시퀀싱(NGS) 기술의 도입은 미생물 분리주의 유전체 특성 분석의 문을 열었습니다. 전체 유전체 시퀀싱(WGS)은 새로운 병원체, 약물 내성 균주, 백신 회피 변이의 출현, 확산 및 전파를 추적할 수 있게 합니다. 또한 임상 발병 시 전파 경로를 식별하는 데 점점 더 활용되고 있으며, 이는 향후 발병 예방을 위한 개입에 도움이 될 가능성이 있습니다. 7. 마지막으로, WGS는 ARG의 프로파일링을 가능하게 하여 치료 결정에 도움을 줄 수 있으나, ARG를 항생제 내성 표현형과 정성적 및 양적 연계하는 작업은 아직 진행 중입니다.
하지만 WGS가 대규모로 생산하는 단독(SR) 시퀀싱 데이터는 일반적으로 400–500 bp 범위에 속합니다. 이 크기는 대부분의 이동성 유전 요소(MGE)보다 짧아10을 분리할 수 없습니다. WGS 데이터에서 플라스미드 서열의 실리코 내 식별은 불완전한 범위와 플라스미드 간 및 박테리아 염색체와 다른 플라스미드 간의 서열 공유 때문에 더욱 복잡합니다10. 그 결과, 조립체는 종종 광범위하게 분해되고 불완전하여 그 하위 유전체 위치(염색체 또는 플라스미드)를10으로 귀속시키는 것이 불가능합니다.
롱리드(LR) 시퀀싱 기술은 더 많은 반복 과정에 걸쳐 매우 연속적인 조립체를 생성하여 완전한 플라스미드 조립체의 생산을 크게 용이하게 합니다. 본질적으로 두 가지 유형의 장독 단일 분자 기술이 존재합니다11,12. 이러한 접근법 중 하나인 나노포어 시퀀싱은 저렴한 비용, 빠른 처리 시간, 자원 제한 환경에 적합하다는 점으로 인해 미생물 유전체학에서 점점 인기를 얻고 있습니다13,14.
나노포어 시퀀싱 장치는 수천 개의 나노포어 배열을 포함하는 플로우셀을 사용합니다. DNA는 헬리케이스 활성을 가진 모터 단백질에 의해 열린 플로우 셀에 적재되며, 전류 차동에 의해 구동되는 나노 구멍을 통과합니다. 특정 단일 가닥의 뉴클레오타이드가 개별 나노기공을 통과할 때, 전류가 교란되어 각 기공에 연결된 전자 센서로 측정되는 특성 신호가 생성됩니다. 이 신호는 베이스콜링 알고리즘을 사용해 디코딩되고, 그 결과 데이터는 다양한 생물정보학 도구를 통해 하위 단계에서 수집 및 분석됩니다.
SR 시퀀싱 접근법과 비교할 때, 나노포어 시퀀싱은 낮은 읽기 품질(낮은 Q-점수)을 생성하고, 전기 신호 데이터 분석 중 동질중합체 서열에 삽입/결실을 발생시키는 경향이 있습니다(베이스콜링)15. 또한 자연 DNA에서 근처의 메틸화가 체계적인 염기 호출 오류를 일으킬 수 있다는 점도 관찰되었습니다16. 일반적으로 채택된 해결책은 조립 후 오류 수정을 위한 단판 데이터를 포함하는 것이었습니다. 하지만 이 과정('폴리싱'이라 불림)은 비용과 복잡성오버헤드를 증가시킵니다. 시퀀싱 커버리지를 늘리는 것도 정확도를 높일 수 있지만, 기술적으로 항상 가능한 것은 아닙니다.
2024년에는 옥스퍼드 나노포어 시퀀싱(ONT) R10.4.1 플로우 셀과 V14 화학이 G-균에 대해 고품질의 연속적인 조립체를 생성하는 반면, 단일 염기체 수준에서의 더 정밀한 분석도 SR 시퀀싱 데이터의 이점을 얻는 것으로 밝혀졌습니다18. 하지만 LR과 SR 시퀀싱 정확도 간의 격차는 점차 줄어들고 있습니다. 최근 다양한 개선 사항으로는 더 정밀한 데이터 샘플링 과정과 향상된 시퀀싱 속도, 천연 박테리아 DNA와 메틸화에 맞춘 새로운 염기호출자 도입(염기호출 정확도 향상)19, 그리고 위상 측정, 변이 호출, 별-대립유전자 호출에 관한 새로운 알고리즘 등 분석 도구의 업데이트가 포함됩니다. 또한 최신 Dorado 알고리즘 V.0.9.0 ONT (2025)는 원시 스퀴글 데이터의 리베이스 호출에 사용할 수 있다는 점에 유의하세요. 실제로 일부 WGS 연구에서는 이미 나노포어 시퀀싱을 통해 미생물 게놈이 생성되었으며, 이 품질은 단독 연마가 합의 서열을 크게 개선할 것으로 기대되지 않는다고 보고하고 있습니다20,21.
플라스미드는 특히 장내세균과(Enterobacteriaceae) 분리주에서 획득된 내성 유전자의 대부분을 포함하는 경우가 많습니다. 따라서 플라스미드 서열을 분해하면 AR 유전자의 맥락화가 가능하며, 이는 AMR 전파 추적에 핵심적인 역할을 합니다23,24. 나노포어 시퀀싱은 이미 blaKPC가 플라스미드를 보유한 미세진화를 조사하는 데 사용되었습니다. MinION 장치는 또한 병원 환경에서 임상 분리주에서 플라스미드의 실시간 벤치 상단 시퀀싱과 내성 유전자 검출에 사용되었습니다26. 그렇긴 하지만, 플라스미드 심문은 아직 일상적인 감염 통제 관행의 일부가 아닙니다.
LR 시퀀싱이 유전체서열 27에서 플라스미드를 구출하는 능력에서 SR 시퀀싱에 비해 크게 향상되었음에도 불구하고, 여전히 몇 가지 문제가 남아 있습니다. Flye28, Canu29 , Raven30 과 같은 롱리드 어셈블러는 고품질 염색체 조립체를 생산하지만 플라스미드 서열을 조립하는 데는 어려움을 겪습니다. 이는 부분적으로 플라스미드 DNA가 WGS 시퀀싱 데이터에서 과소대표되어 있어 범위 깊이가 부족하기 때문입니다31. 모듈형 구조, 풍부한 2차 구조, 반복 반복 현상으로 인해 플라스미드 DNA 서열의 고품질 조립은 염색체 DNA보다 더 높은 커버리지가 필요하다는 점에 유의하라(50배 읽기 깊이32).
플라스미드 서열로 DNA 샘플을 풍부하게 하는 것은 덮개 깊이를 높이는 효과적인 방법이지만, DNA 시퀀싱 전에 플라스미드 DNA를 분리하거나 풍부하게 하는 것은 임상 진단용으로 너무 비용이 많이 들거나 번거로울 수 있으며, 선택된 DNA 서열을 풍부하게 하여 편향을 초래할 수 있습니다. 2017년에는 임상 샘플에서 플라스미드 DNA의 완전한 조립을 시도하는 데 플라스미드 풍부 DNA를 사용하여ARGs 24 검출을 시도했습니다. 6종의 서로 다른 장관박테리아과에서 분리된 8종을 시퀀싱하여 개체군과 플라스미드 구조 다양성을 극대화하였습니다. 여러 조립기를 시도한 결과, Canu/Canu+Pilon이 다른 어셈블러보다 확실히 우수한 성능을 보였으며; 이 연구팀은 또한 메타 어셈블리 접근법을 시도하여 추가 플라스미드 구조의 해상도를 달성하여, 시퀀싱된 샘플 내 알려진 모든 플라스미드의 78%를 완전하게 조립하는 결과를 얻었습니다. 그러나 대규모 반복 구조의 존재와 (다소 적지만) DNA 추출 문제는 모든 단일 조직 플라스미드 조립체 검출에 상당한 어려움을 남겼다24.
최근에는 저점유 플라스미드 서열의 풍부화 보완적 접근법으로서 적응 샘플링이 박테리아 분리 샘플에서 염색체 서열을 거부하는 방식으로 연구되었으나, 이 방법은 알려진 참조 서열이 필요하다33.
여기서 저자들은 플라스미드 DNA 농축 샘플에 맞춘 ONT 시퀀싱 프로토콜을 설명합니다. 플라스미드 DNA의 표현성 증가와 최신 라이브러리 준비 프로토콜, 플로우 셀 기술, 서열 조립 도구를 결합하여, 이 워크플로우는 겔 이동성에서 추정된 크기와 일치하는 조립체를 생성했으며, 3분의 2의 경우에는 SR 시퀀싱 연마와 일치하거나 0.05% 쌍별 서열 동일성 내에 도달했습니다. 플라스미드 DNA 서열의 풍부화는 플라스미드 시퀀싱의 효율성을 극대화하여 최대 24개의 샘플을 최적의 커버리지로 병렬 시퀀싱할 수 있게 합니다. 이 프로토콜들은 표준 나노포어 WGS 접근법을 대체하려는 것이 아니라, 염색체 서열이 이미 알려져 있거나 관심 없는 상황에서 플라스미드 시퀀싱의 정확성과 효율성을 극대화하기 위한 것입니다.
워크플로우에는 플라스미드 DNA 추출, 바코딩을 이용한 라이브러리 준비, 모터 단백질과의 어댑터 연결, 플로우셀 점검, 라이브러리를 플로우셀에 적재하여 시퀀싱, 전기 신호 데이터 분석(실시간 또는 사후 분석), Autocycler34를 이용한 서열 조립, 그리고 다양한 후속 분석(그림 1)이 포함됩니다.

그림 1: 옥스퍼드 나노포어 테크놀로지를 이용한 플라스미드 시퀀싱 워크플로우. 1단계: DNA 추출; 2단계: 도서관 준비; 3단계: 시퀀싱; 4단계: 조립 및 주석 작성. BioRender로 제작되었습니다. 코르테스, G. (2026) https://BioRender.com/u3b15st. 이 그림의 더 큰 버전을 보시려면 여기를 클릭해 주세요.
1. 플라스미드 DNA 추출
참고: 이 프로토콜은 상용 플라스미드 DNA 추출 키트35를 사용하여 임상용 대장균과 해당 트랜스콘가간에서 플라스미드 DNA를 추출하는 데 사용되었습니다. 활용에 대한 자세한 프로토콜은36에서 확인할 수 있습니다. 이 프로토콜은 ~3 kb에서 ~200 kb까지 고복사 및 저 복사 플라스미드 모두에 적합합니다. 대형 플라스미드는 복사 수가 적어 세포당 1-2개 정도 DNA 생산량이 제한되므로 특별한 주의가 필요합니다. 이를 보상하기 위해 프로토콜은 세포 질량과 용해 및 중화 완충제의 부피를 증가시켜 컬럼 결합에 필요한 충분한 플라스미드 DNA를 확보할 것을 권장합니다. 이 프로토콜을 사용하면 최대 200 kb의 플라스미드도 여전히 복구할 수 있습니다; 하지만 메가 플라스미드(>200 kb)에는 적용 범위가 제한적입니다.



그림 2: 원심분리 한 단계 후 흐린 용해상 상액. 여과 컬럼이 막히지 않도록 상원액이 맑아질 때까지 연속 원심분리가 필요합니다. 이 그림의 더 큰 버전을 보시려면 여기를 클릭해 주세요.
| 표준 분석 튜브 | 사용자 샘플 분석관 | |
| 작업 용액 (염료 및 버퍼, 1:200) | 190 μL | 180–199 μL |
| 스탠다드 (키트 기준) | 10 μL | – |
| 사용자 샘플 | – | 1–20 μL |
| 각 분석관 내 총 부피 | 200 μL | 200 μL |
표 1: 형광광량계용 시약 부피 및 양. 단일 샘플에 대한 양이 표시되며; 분석할 샘플 수에 따라 풀을 준비할 수 있습니다.
2. 방법 2: 도서관 준비
참고: 고려사항 추출한 플라스미드의 시퀀싱을 위한 라이브러리 준비는 바코딩 키트40에 제공된 지침에 따라 수행되었습니다. 이 프로토콜은 전이효소 기반 방식을 사용하여 DNA를 동시에 단편화하고 바코드를 처리하여, 빠르고 효율적인 워크플로우를 제공하여 직접 작업하는 시간을 줄이고 도서관 준비를 단순화합니다.
| 시약 | 상온에서 해동 | 스핀 다운 | 피펫팅으로 혼합 |
| 빠른 바코드 | 겨울왕국이 아니야 | ✓ | ✓ |
| 래피드 어댑터 | 겨울왕국이 아니야 | ✓ | ✓ |
| AMPure XP 비즈 | ✓ | ✓ | 사용 직전에 피펫팅 또는 소용돌이 처리로 혼합하세요 |
| 용출 버퍼 | ✓ | ✓ | ✓ |
| 어댑터 버퍼 | ✓ | ✓ | 소용돌이 혼합 |
표 2: 나노포어 라이브러리 준비를 위한 시약 부피 및 양. 각 시약에 대해 수행할 동작은 (✓)로 표시됩니다.
| 시약 | 볼륨/샘플 |
| 10 μL 템플릿 DNA (이전 단계에서 200 ng) | 10 μL |
| 1.5 μL 급속 바코드, 각 샘플마다 하나씩 | 1.5 μL |
| 총 부피 | 11.5 μL |
표 3: 빠른 바코딩을 위한 시약 및 부피. DNA 샘플당 부피가 표시되며; 풀링된 혼합물을 준비하고 DNA를 별도로 첨가할 수 있습니다.
| 1-24 바코드 | 48개의 바코드 | 72개의 바코드 | 96개의 바코드 | |
| 용출 버퍼 부피 | 15 μL | 30 μL | 45 μL | 60 μL |
표 4: 24개의 바코드당 사용되는 용출 버퍼 볼륨. 부피는 바코드 수에 따라 표시됩니다.
| 시약 | 볼륨 |
| 래피드 어댑터(RA) | 1.5 μL |
| 어댑터 버퍼 (ADB) | 3.5 μL |
| 총 부피 | 5μL |
표 5: 어댑터용 시약 및 양. 샘플 수에 따라 볼륨을 조절할 수 있습니다.
3. 플로우 셀 실행
참고: 중간 정도의 범위가 필요한 소수의 플라스미드를 시퀀싱할 때는 표준 출력 플로우 셀을 선택하세요. 나노포어는 활성 기공 수를 확인하기 위해 플로우 셀 검사를 수행할 것을 권장합니다; 나노포어 장치의 플로우셀은 보증 기준에 따라 최소 800개의 활성 기공을 가져야 합니다. 단일 MinION 플로우셀은 플라스미드 DNA로 풍부해진 여러 DNA 샘플의 시퀀싱에 충분한 출력을 제공하며(다중화에 관한 자세한 내용은 토론 섹션에서 참고), GRIDion이나 Promethean 플로우셀보다 더 낮은 비용과 더 큰 유연성을 제공합니다.
| 시약 | 플로우 셀당 부피 |
| 플로우셀 플러시(FCF) | 1,170 μL |
| 소 혈청 알부민(BSA) 50 mg/mL | 5 μL |
| 플로우 셀 테더(FCT) | 30 μL |
| 최종 튜브 내 총 부피 | 1,205 μL |
표 6: 플로우셀 프라이머 및 로딩을 위한 시약 및 부피. 각 플로우 셀마다 시약의 부피가 표시됩니다.
| 시약 | 플로우 셀당 부피 |
| 시퀀싱 버퍼(SB) | 37.5 μL |
| 라이브러리 비즈(LIB) 또는 라이브러리 솔루션(LIS), | 25.5 μL |
| DNA 라이브러리 | 12 μL |
| 총 | 75 μL |
표 7: 도서관 로딩을 위한 시약 및 부적. DNA 도서관의 권은 표시되어 있으며; 라이브러리 비즈(입찰)는 사용 직전에 반드시 섞어야 합니다.
4. 방법 4: 오토사이클러 조립
참고: 리드를 컨티그로 조립하는 방법은 여러 가지가 있습니다. 각 어셈블러는 장점과 편향이 있습니다. 이 프로토콜은 여러 어셈블리를 하나의 어셈블리로 결합하는 Autocycler 사용을 권장합니다. 이 방법은 더 많은 노동 집약적이고, 단일 조립기로 시간을 절약할 수 있습니다. 적절한 소프트웨어를 다운로드하려면 터미널에 Conda를 설치해 주세요. 보조 명령어 문서에는 시스템 설정을 위한 환경 파일 링크가 포함되어 있습니다.
5. 방법 5: 주석
참고: 박테리아에 대한 효과적인 주석 도구가 많이 있습니다. 자세한 주석을 위해 다음 프로그램들이 권장됩니다.
임상 대장균 균주
임상 균주 플라스미드 시퀀싱에 ONT를 사용하는 방법을 설명하기 위해, 저자들은 워싱턴 대학교의 스티븐 살리판테 박사가 후원하는 10개 균주를 시퀀싱했습니다. 이 균주는 워싱턴 대학교 의료센터45에서 일상 임상 진료 중 얻은 312개의 혈액 또는 소변 유래 대장균(ExPEC) 컬렉션의 일부였습니다.
플라스미드 접합 및 겔 이동성 분석법을 통한 특성 분석
연구진은 또한 대 장균 LMB100을 수혜자로 사용한 결합 검사도 수행했습니다. 접합은 한 세포(기증자)에 존재하는 플라스미드가 다른 세포(수혜자라고 불림)로 전달되는 과정입니다46. 이 이식은 두 세포 간의 물리적 접촉이 필요하며, 일반적으로 수용 세포 내에 접합성 플라스미드가 존재함을 나타내는 마커를 사용하여 검출됩니다. 이를 통해 저자들은 숙주 균주 차이를 변수로 제거하고, 플라스미드 포획(특정 생태적 지위에서 플라스미드를 프로파일링하는 방법)로 얻은 플라스미드 서열 분석 프로토콜의 사용을 설명할 수 있었습니다. 모든 결합 분석은 트랜스콘커건트를 생성하여 기증자 내에 접합 플라스미드가 존재함을 확인했습니다. 저자들이 따랐던 프로토콜은36쪽에 설명되어 있습니다. 연구팀은 기증자와 수혜자 균주에 존재하는 플라스미드를 확인하고, 펄스장 겔 전기영동(PFGE)을 통해 그 이동성을 추정했습니다. 표 8 에는 대표적인 10가지 공여 균주 목록이 나와 있으며, 이들 중 플라스미드의 추정 크기는 세 번째 열에 나와 있습니다. PFGE 겔의 해상도 범위를 벗어난 플라스미드는 이 표에서 제외되지만, 연구진은 약 ~4kb 정도의 작은 플라스미드 여러 개를 조립할 수 있었습니다(사진에는 표시되지 않음).
| 기증자 계통 | 플라스미드 | 추정 PFGE 크기 | 조립 길이 | 평균 깊이 | bp 제로 깊이 | 보장 % | 입장 변경 | SR 이전 정확도 % |
| blood_08 _0081 | pblood_08_ 0081_1 | 170000 | 171829 | 100 | 99.941803 | 54 | 0.0314 | 99.9686 |
| blood_08 _0081 | pBlood_08_ 0081_2 | 82000 | 85576 | 46.9 | 0 | 100 | 2 | 99.9977 |
| blood_d-08 _0094 | pBlood_d-08 _0094 | 179000 | 173334 | 382.6 | 3121 | 98.19943 | 2,347 | 98.646 |
| blood_08 _1447 | pBlood_08 _1447 | 97000 | 99624 | 36.2 | 0 | 100 | 0 | 100 |
| blood_10 _0913 | pBlood_10 _0913_1 | 70000 | 71788 | 49 | 4 | 99.994428 | 188 | 99.7381 |
| blood_10 _0913 | pBlood_10 _0913_2 | 112000 | 118727 | 41.2 | 126 | 99.893874 | 276 | 99.7675 |
| blood_1 1_184 | pBlood_11 _184_1 | 150000 | 146583 | 65.1 | 287 | 99.804206 | 134 | 99.9086 |
| blood_11_184 | pBlood_11 _184_2 | 75000 | 74263 | 82.1 | 0 | 100 | 32 | 99.9569 |
| blood_2011 _0238 | pBlood_ 2011_0238 | 170000 | 167286 | 110.9 | 0 | 100 | 0 | 100 |
| upec_108 | pUPEC_108 | 130500 | 131039 | 103.9 | 0 | 100 | 0 | 100 |
| upec_134 | pUPEC_134 | 160500 | 157829 | 169 | 0 | 100 | 0 | 100 |
| upec_271 | pUPEC_271 | 145500 | 144519 | 86 | 33311 | 76.950436 | 1,557 | 98.9226 |
| upec_90 | pUPEC_90 | 70000 | 74624 | 139.8 | 0 | 100 | 6 | 99.992 |
표 8: 10개 기증 균주 패널의 플라스미드 크기 추정 및 시퀀싱 정확도 지표. PFGE 또는 조립된 ONT 공여 플라스미드 서열을 기반으로 한 플라스미드 크기 추정치가 나열되어 있습니다. ONT 단단과 SR-서열 연마 플라스미드 조립체의 쌍별 비교도 정확도의 지표로 제시됩니다. 지표로는 어셈블리 길이, 평균 깊이, 제로 깊이 영역, 커버리지 비율, 위치 변경, SR 이전 정확도 비율 등이 포함됩니다.
ONT 기반 플라스미드 시퀀싱 출력 지표
시퀀싱 런은 24시간으로 설정되었습니다. 총 24개의 샘플이 다중화되었습니다. 각 샘플의 기술 데이터, 품질 점수(Q), 길이, 커버리지 등은 표 9에 나와 있습니다. 나노기공 시퀀싱은 분석된 모든 샘플에 대해 충분한 읽기 출력과 품질을 제공했습니다. 모든 라이브러리를 통틀어 시퀀싱은 샘플당 평균 134,115개의 리드를 생성했으며, 평균 수율은 576.2 Mb에 해당합니다. 리드의 평균 길이는 8,797 bp, 평균 N50은 3,693 bp로, 플라스미드 조립에 충분한 리드 연속성을 제공했습니다. 평균 읽기 품질 점수는 Q16.0, 중앙값은 Q16.6이었고, 개별 런은 Q13.7에서 Q18.0 사이였습니다.
| 서열 변형 | 읽기 (kb) | 베이스 (Mb) | 평균 길이 | 최대 길이 | 최소 길이 | N50 | 중간 길이 | 평균 Q | 중앙값 Q | |
| Blood_08_0081 기증자 | 109.8 | 261 | 5713 | 508418 | 78 | 2377 | 1351 | 14.5 | 15 | |
| Blood_08_0094 기증자 | 89.3 | 229 | 4177 | 114948 | 81 | 2565 | 1065 | 13.7 | 14 | |
| Blood_08_1447 기증자 | 87.5 | 246 | 5041 | 135147 | 75 | 2816 | 1609 | 15.1 | 15.7 | |
| Bblood_10_0913 기증자 | 247.7 | 1026 | 9260 | 174231 | 70 | 4143 | 1661 | 15.5 | 16.1 | |
| Blood_11_154 기증자 | 51.5 | 110 | 4172 | 91148 | 73 | 2127 | 1387 | 16.1 | 16.5 | |
| Blood_2011_0238 기증자 | 16.8 | 42 | 6707 | 104848 | 62 | 2534 | 845 | 15.3 | 15.2 | |
| UPEC_108 기증자 | 50.6 | 77 | 1915 | 90234 | 67 | 1525 | 1103 | 15.8 | 16.2 | |
| UPEC_134 기증자 | 50 | 143 | 7110 | 158393 | 72 | 2865 | 1133 | 16 | 16.2 | |
| UPEC_217 기증자 | 64 | 236 | 7512 | 144275 | 78 | 3678 | 1895 | 17.4 | 17.9 | |
| UPEC_90 기증자 | 52.9 | 157 | 5666 | 95944 | 72 | 2969 | 1593 | 16.8 | 17.2 | |
| Blood_08_0081 LMB100 환합 전합체 | 53.5 | 119 | 5162 | 225856 | 1 | 2232 | 1009 | 16 | 16.8 | |
| Blood_08_0094 LMB100 환합부 | 160 | 442 | 7257 | 252176 | 1 | 2761 | 1083 | 17.3 | 18.2 | |
| Blood_08_1447 LMB100 전공부액 | 238 | 1029 | 10249 | 209221 | 1 | 4320 | 1857 | 17.7 | 18.7 | |
| Blood_10_0913 LMB100 환결합 | 792.3 | 4341 | 11508 | 474325 | 1 | 5479 | 2731 | 17.7 | 18.8 | |
| Blood_11_184 LMB100 환합부 | 42.5 | 281 | 17931 | 146626 | 87 | 6615 | 2448 | 16 | 16.7 | |
| Blood_2011_0238 LMB100 환합부 | 147.7 | 585 | 9000 | 344203 | 1 | 3962 | 1781 | 18 | 18.9 | |
| UPEC_108 LMB100 환합 전공액 | 173.4 | 778 | 10127 | 168756 | 80 | 4485 | 2079 | 15.3 | 16.1 | |
| UPEC_134 LMB100 환결합 | 105.4 | 667 | 17568 | 349230 | 71 | 6323 | 2232 | 15.5 | 16.3 | |
| UPEC_271 LMB100 환합 | 74.2 | 336 | 14521 | 342195 | 83 | 4520 | 1306 | 14.6 | 14.8 | |
| UPEC_90 LMB100 환결합 | 75.2 | 418 | 15348 | 225661 | 83 | 5561 | 2016 | 15.7 | 16.4 | |
표 9: 10개 기증 균주 패널에 대한 NanoQ 시퀀싱 출력 지표 보고서. 시퀀싱 출력 지표가 나열되어 있습니다. 여기에는 리드, 염기, 길이(평균, 중앙값, 최대값, 최소값), 길이, N50, 그리고 기증자 및 환합자로부터 얻은 Q(평균 및 중앙값) 값이 포함됩니다.
플라스미드 조립 및 정확도
플라스미드 서열은 Autocycler34를 사용해 조립되었습니다. 10개 샘플의 조립된 플라스미드 서열의 크기는 표 8 네 번째 열에 나와 있으며, 모든 경우에 PFGE 이동도를 기반으로 한 추정치(세 번째 열에 나열됨)와 일치하였습니다. 정확도는 Medaka 연마된 롱리드 어셈블리와 Medaka(v2.2.0) 및 폴리폴리시(v0.6.1) 어셈블리를 비교하여 평가되었습니다. 서열 동일성은 MUMmer(v3.2.4)의 dnadiff를 사용해 계산하였다(표 8). 분석된 모든 플라스미드(n = 12)를 통틀어 조립체의 평균 길이는 120.4 kb였고, 평균 시퀀싱 깊이는 109배로 지원되었으며, 플라스미드 서열에 걸쳐 평균 97.9%의 적용률을 보였습니다. 연마 이전에는 조립체의 평균 서열 동일성이 SR 연마 조립체에 비해 99.74%였습니다.
염색체 오염 추정
플라스미드 범위와 염색체 오염을 추정하기 위해 레딩을 조립된 플라스미드 서열과 대장균 LMB100 염색체 서열로 매핑하였다. Basecalled Nanopore 리드는 minimap2(v2.2.6)를 사용하여 정렬되었습니다. 결과된 정렬은 SAMtools(v1.1.7)를 사용하여 정렬된 BAM 파일을 생성하고 매핑 통계를 계산했습니다. 플라스미드 커버리지와 평균 깊이는 SAMtools 커버리지를 사용하여 추정하였습니다. 플라스미드 시퀀싱 깊이는 샘플마다 달랐으나, 정확한 조립에 필요한 최소 커버리지(50–70 x)를 훨씬 상회했습니다. 리드 중 일부는 염색체에 매핑되어 플라스미드 추출 중 숙주 DNA의 운반을 나타냅니다(표 10).
| 플라스미드 | 총 읽기 횟수 | 염색체 판독 | % 염색체 판독 | 평균 염색체 커버리지 | 평균 플라스미드 적용 범위 | P:C 비율 |
| pBlood_08_0081_2 | 53491 | 20116 | 37.61% | 5 | 1068 | 236 |
| pBlood_d-08_0094 | 160024 | 74827 | 46.76% | 37 | 1618 | 44 |
| pBlood_08_1447 | 238031 | 127967 | 53.76% | 62 | 7842 | 127 |
| pBlood_10_0913_1 | 792263 | 221610 | 27.97% | 87 | 7801 | 90 |
| pBlood_10_0913_2 | 792263 | 221610 | 27.97% | 87 | 7772 | 89 |
| pBlood_11_184_1 | 42494 | 15919 | 37.46% | 7 | 1637 | 250 |
| pBlood_11_184_2 | 42494 | 15919 | 37.46% | 7 | 2188 | 335 |
| pBlood_2011_0238 | 147680 | 80526 | 54.53% | 56 | 2965 | 53 |
| pUPEC_108 | 173409 | 107884 | 62.21% | 61 | 4008 | 65 |
| pUPEC_134 | 105418 | 68645 | 65.12% | 29 | 3640 | 127 |
| pUPEC_271 | 74232 | 53554 | 72.14% | 34 | 1202 | 35 |
| pUPEC_90 | 75193 | 30052 | 39.97% | 16 | 4488 | 278 |
표 10: 염색체 DNA 오염 추정치. 수용자 균주에서 정제된 플라스미드의 염색체 오염은 플라스미드 리드에 비해 LMB100 염색체 서열에 매핑되는 리드의 비율로 정량화되었습니다. 각 플라스미드에 대해 플라스미드 리드, 총 리드, 염색체 리드, 평균 염색체 커버리지, 평균 플라스미드 커버리지, P:C 비율이 나열되어 있습니다. pBlood_08_0081_1 불필요한 이유는 수혜자에서 회수할 수 없었기 때문이며, 이 추정은 수혜자의 균일한 염색체 배경 때문에 수행되었다는 점에 유의하세요.
염색체 오염은 조립된 LMB100 염색체 서열에 매핑된 리드의 비율과 플라스미드 리드 간의 비율로 정량화되었습니다. 염색체 오염률은 72.1%에서 28.0% 사이였으며, 평균 46.9%였습니다(표 10). 이 방법은 플라스미드와 염색체 모두에 대한 매핑 리드가 염색체 전용 것으로 간주되어 염색체 오염을 과대평가할 가능성이 있지만, 상당한 염색체 운반을 나타냅니다. 염색체 운반에도 불구하고, 조립 작업 흐름은 모든 샘플에서 PFGE 이동성에서 예측된 크기와 일치하는 완전한 플라스미드 서열을 성공적으로 재구성했습니다(표 8, 3, 4열).
대표적인 기증자-수혜자 결합 플라스미드 쌍
업펙 균주는 2011년부터 2013년 사이에 요로 샘플에서 분리되어 계통그룹 B2와 서열 유형 ST131로 분류되었습니다. PFGE 분석 결과 약 145 kb 크기의 플라스미드에 해당하는 밴드가 존재함을 확인했습니다 [그림 3A, (기증자)]. 단독 서열 분석을 통해 IncF 복제체의 존재가 확인되었습니다.

그림 3: pUPEC271의 특성 분석. (A)기증자와 환결합액의 PFGE. 공여자(D)와 환결합부(T)의 PFGE 밴드가 표시되어 있습니다; 초록색 화살표는 pUPEC271에 해당하는 플라스미드 밴드를 나타냅니다. PFGE에 사용되는 프로토콜은 연구36에서 설명되어 있습니다. (B)pUPEC271 지도. 트랜스콘가겐 이 기증자에게 강렬하게 맞춰졌고, 100% 동일한 컨센서스가 초록색 링으로 표시됩니다. 이 표현은 Proksee48을 사용하여 생성되었습니다. 가장 바깥쪽 고리는 CDS 위치를 나타내고, 두 안쪽 고리는 GC 함량과 왜도를, 가장 안쪽 고리는 tRNA 유전자의 위치를 나타냅니다. 이 그림의 더 큰 버전을 보시려면 여기를 클릭해 주세요.
저자들이 pUPEC271이라 명명한 이 플라스미드의 전이는 앞서36 에서 설명한 프로토콜을 사용한 PFGE 분석을 통해 추가로 확인되었습니다 [그림 3A, T(transconjugant)]. 기증자 upec271 균주(pUPEC271_donor)에서 정제한 플라스미드 DNA는 나노포어를 사용해 시퀀싱되었습니다. PFGE 겔에서 밴드 크기와 일치하여, pUPEC271_donor 서열은 144,519 bp 크기의 플라스미드를 발견했습니다(그림 3B). 이 플라스미드는 IncF형 레플리콘을 가지고 있었으며, 여섯 가지 뚜렷한 항생제 군에 대한 내성 유전자를 보유하고 있었습니다( 보충표 1 참조). 저자들은 또한 트랜스콘이어 인사 중 하나(pUPEC271_TC)에서 DNA를 정제하여 기증자의 DNA와 동일한 서열을 얻었습니다. 최종 조립체는 pUPEC271_donor과 비교되어 그림 3B에 표시되어 있습니다. 이 결과, 기증자에 존재하는 접합 플라스미드가 트랜스인커겐 균주로 이전되었음을 확인하였습니다. 기술 데이터는 표 9에 나와 있습니다.
균주 Blood_10_0913는 여러 플라스미드를 가진 균주의 예입니다. 이 샘플은 2008년부터 2013년 사이에 혈액에서 채취되었으며, 단기 시퀀싱 분석 결과 IncF와 IncY 두 개의 복제 물질이 발견되었습니다. 이 균주에 대한 PFGE 분석에서 각각 약 70kb와 110kb의 두 플라스미드가 확인되었습니다(그림 4A). 이는 나노포어 시퀀싱을 통해 확인되었으며, 각각 71,182 bp와 111,226 bp 크기의 두 조립체가 생성되었습니다. 저자들은 이 두 플라스미드를 pBlood_10_0913_1와 pBlood_10_0913_2로 명명했다(그림 4B). 그러나 조립체들은 접합 기계를 가진 플라스미드 단면에서 광범위한 서열 상동성(광범위한 상동성의 부분은 그림 4B에 강조됨)으로 인해 처음에는 원형화되지 못했다. 붕대47 을 사용해 생성된 실패한 조립 도면은 그림 4C에 나와 있습니다. 두 플라스미드의 완전한 조립은 짧은 리드(< 3000 bp) 서열의 공격적인 리드 컬링 후에야 가능했습니다.

그림 4: Blood_10_0913 균주 플라스미드의 식별 및 조립. (A) PFGE. 초록색 화살표는 플라스미드 밴드를 나타냅니다. PFGE에 사용되는 프로토콜 은36 (B) 서열 조립에 설명되어 있습니다. pBlood_10_0913 샘플에 존재하는 두 개의 작은 플라스미드인 pBlood_10_0913_2와 pBlood_10_0913_2의 조립 과정. (C) 조립 실패. 두 개의 서로 다른 어셈블리 간에 공유 시퀀스가 존재하여 어셈블리 실패를 표현하는 방식. 이 표현은 Bandage47을 사용해 생성되었습니다. 이 그림의 더 큰 버전을 보시려면 여기를 클릭해 주세요.
보충표 1: pUPEC271이 운반하는 내성 유전자. pUPEC271과 관련된 플라스미드 크기, 레플리콘, 내성 유전자가 나타났습니다. 이 파일을 다운로드하려면 여기를 클릭해 주세요.
보조 지휘부. 이 파일을 다운로드하려면 여기를 클릭해 주세요.
플라스미드의 식별 및 특성 분석은 공중보건, 임상 미생물학, 미생물 생태학에서 매우 중요한 과제로, 적응 유전자의 유지 및 확산을 생태 환경 전반에 촉진하고 가속화된 진화의 플랫폼 역할을 하기 때문입니다 1,5. NGS 기술의 도입으로 미생물 분리주와 집단 게놈의 고처리량 시퀀싱이 가능해졌으나, 완성된 플라스미드나 플라스미드 유래 컨티그의 조립 후 식별은 일반적으로 롱리드 서열이 필요한데, 이는 단독 시퀀싱 데이터가 리드 크기보다 큰 반복 영역을 해독하지 못하기 때문입니다. ONT는 전체 유전체 시퀀싱과 메타게노믹스를 넘어 플라스미드 서열 검증, 후생유전학적 변형의 검출 및 특성화, RNA 시퀀싱, 그리고 탠덤 유전자 배열의 내재적 불안정성 평가(12,49에서 검토됨)를 포함하는 매우 인기 있는 장기 독기 시퀀싱 접근법으로 자리 잡았습니다.
여기서 연구진은 염색체 서열이 이미 알려져 있거나 관심 없는 상황을 위해 설계된 플라스미드 ONT 시퀀싱 워크플로우를 제시합니다. 그 이유는 플라스미드 추출이 노동 집약적이고 비용이 많이 들며 시간이 많이 걸리기 때문입니다(숙주 집단을 배양 및 가공해야 하기 때문입니다). 그리고 나노포어 기술을 이용한 전체 유전체 시퀀싱만으로도 완전한 플라스미드 조립체를 얻으면서도 유전체적 맥락을 제공하는 데 충분하기 때문입니다. 프로토콜이 도움이 되는 예로는 플라스미드 유전공학맥락에서 플라스미드 서열 검증, WGS24에서 이미 얻은 플라스미드 서열의 정확도 향상, 임상 또는 환경 샘플에서 접합을 통해 포착된 플라스미드의 특성 규명 등이 있습니다.
임상 및 환경 샘플에서 DNA 추출은 일부 맞춤화가 필요할 수 있습니다. 여기서 저자들은 사용한 일부 임상 균주가 필터를 막는 바이오필름 매트릭스를 생성해 여러 번의 펠릿 부착 단계를 필요로 했고, 때로는 두 개의 필터를 추가해야 했다고 보고합니다. 또한, 접합 플라스미드를 시퀀싱하기 위해서는 다중 복제 플라스미드 시퀀싱에 비해 더 큰 배양 부피가 필요한데, 자연적으로 존재하는 중형 및 대형 플라스미드는 플라스미드 복제 번호2가 낮은 경향이 있기 때문입니다.
이 프로토콜은 4 kb에서 173 kb 크기의 플라스미드를 시퀀싱할 수 있게 했습니다. 이 범위는 제조사 사양(2에서 200kb 사이)에 정확히 부합합니다. 최근 NCBI에 침착된 23,000개의 플라스미드 서열을 조사한 연구에서는 출처에 따라 크기 차이가 크게 나타났으며, 인간기원의 플라스미드는 가장 작은 중간 크기(76 kb)를 가졌고, 토양과 식물 플라스미드는 각각 215 kb와 427 kb로 가장 컸습니다; 다른 모든 출처의 플라스미드는 중앙값이 79에서 147 kb 사이였으므로, 이 프로토콜은 대부분의 플라스미드를 포함할 수 있으며, 메가플라스미드52 와 토양 또는 식물 기원 플라스미드3을 제외하면 예외가 될 수 있습니다.
저자들은 플라스미드 DNA 농축 단계를 포함하여 플라스미드 DNA 조립체의 정확도를 높이고 있습니다. 대량의 염색체 DNA가 존재하면 시퀀싱 플로우 셀의 세포 구멍에 비표적 리드가 퍼져 플라스미드 DNA 깊이가 크게 감소합니다. 임상 샘플에서ARGs 24 검출을 위해 플라스미드 DNA 완전 조립을 시도하기 위해 플라스미드 DNA 농축 단계가 이미 시행되었습니다. 이 결과는 나노포어 단독 서열로 플라스미드 조립을 완료했을 때 효과가 크게 향상되었음을 보여주었습니다(벤치마크 대비 78% 완성). 이후50, 53, 54 등 다른 예시들도 있었습니다.
읽기 정확도를 높이기 위해 이 중 두 연구는 나노포어가 개발한 '의사쌍 리드에 대한 쌍 베이스콜링' 기법을 사용했습니다50,53. 이 방법은 동일한 DNA 분자의 센스 가닥과 안티센스 가닥에서 나오는 원시 전기 신호를 정렬하여 정확도가 약 100배 향상되는 것을 포함합니다. 의사쌍 리드에서의 쌍 염기호출은 같은 분자의 순방향 가닥과 역방향 가닥을 식별할 수 있는 전략이라면 다중화와 호환됩니다. 따라서 알려지지 않은 플라스미드 보체 샘플에서 정제되거나 플라스미드 혼합물이 포함된 플라스미드에는 사용할 수 없습니다. 유전체 내 플라스미드 수는 일반적으로 0에서 7 사이(특히 임상 분리주에서)로 범위를 가진다는 점에 유의하세요.
연구진은 샘플에서 높은 수준의 염색체 오염을 발견했습니다(표 10). 분석 대상 10개 샘플의 오염 수준은 리드 매핑에 따라 28%에서 72% 사이였습니다. 그러나 연구자들이 이 서열을 염색체 크기에 맞게 정규화했을 때, 플라스미드 서열이 염색체 서열보다 평균 144배 더 풍부하다는 사실이 밝혀졌습니다. 이 수치를 플라스미드 DNA 농축 수준으로 변환하려면 각 샘플 내 플라스미드 복사 수를 알아야 합니다.
MinION 플로우 셀은 일반적으로 ~30GB의 시퀀싱 출력을 생성하여, 이론적으로 한 번에 96개 이상의 플라스미드를 시퀀싱할 수 있습니다. 실제로 다중화는 96바코드 제한에 의해 제약받기 때문에, 플로우 셀당 실질적으로 96개의 플라스미드가 실질적인 상한선이 됩니다. 그러나 염색체 DNA 오염과 자연스러운 바코드 불균형은 샘플당 유효 읽기 깊이를 감소시킵니다. 각 샘플에 충분한 독독 깊이를 보장하기 위해 보다 보수적인 접근법은 플로우 셀당 약 24개의 플라스미드를 다중화하는 것입니다. 플라스미드 서열의 정확한 조립을 위한 권장 깊이는 50에서 70x32.56 사이입니다. 이 글에서 제시한 예시들에서, 저자들은 한 플로우 셀에서 21개의 다중화 샘플을 실행하여 평균 깊이(109x)를 얻었는데, 이는 플로우 셀별로 권장되는 샘플 최적값과 잘 일치합니다.
읽기 깊이가 과도할 경우, 통과/실패 읽기 파티션은 전체 읽기 품질 평가에 도움이 되지 않는다는 점에 유의하세요. 대신 길이 분포에서 도출된 지표는 데이터셋 품질을 더 명확하게 나타냅니다. 그 중 하나가 N50으로, 전체 염기의 50%를 차지하는 읽기 길이로 정의됩니다. 이 값은 품질 관리 단계에 매우 민감합니다. 예를 들어, Filtlong과 같은 도구로 저품질 또는 짧은 리드를 필터링하면 분포의 짧은 꼬리 부분의 상당 부분이 제거되어 N50이 부풀려집니다. 이 효과는 표 9에서 명확히 나타납니다.
저자들은 ONT의 Rapid Barcode 키트를 사용했는데, 이 키트는 트랜스포좀을 이용해 플라스미드 DNA를 분해하여 플라스미드의 제한 패턴에 대한 사전 지식이 필요하지 않습니다. 이전 연구에서는 소형 플라스미드(<20 kb)의 경우, 소플라스미드는 대부분의 DNA 추출 시 원형으로 남아 있어 어댑터 연결(adapter ligation 54)을 위한 자유 말단이 없기 때문에 효소 단절이 필요하다는 것을 발견했습니다. 그 연구는 PCR 증폭을 통해 추가 자유 말단을 생성하는 바코드를 사용하지 않았기 때문에, 이 경우에는 아마도 논쟁의 의미가 없을 것입니다. OnRamp는 일상적인 플라스미드 검증을 위한 고급 절차 중 하나로, 조립을 위해 전체 길이의 플라스미드 리드를 활용해 바코드 사용을 피하여 샘플 준비과정을 단순화합니다. 이 경우 저자들은 세 가지 이유로 바코드를 선호했습니다. 첫 번째는 플라스미드 길이와 상관없이 조립체를 만드는 것입니다. 두 번째는 염색체 DNA 오염을 견디는 것입니다. 세 번째는 앞서 언급한 대로 여러 샘플을 다중화할 수 있는 가능성을 열어준다는 점입니다.
이전 연구에서는 서로 다른 조립기의 집합이 개별 조립 플랫폼에 내재된 편향을 상쇄하여 최종 결과를 개선한다고 지적했다24. 이 메타분석은 저자들이 사용하는 어셈블리 프로그램인 Autocycler34에 통합된 기능입니다. 이전 버전인 트라이사이클러57에 비해 오토사이클러는 다소 더 자동화되어 있습니다; 읽기 중 여러 하위 집합을 사용하며, (앞서 언급했듯이) 여러 어셈블리를 사용해 더 정확한 합의 어셈블리를 만듭니다. 그러나 대부분의 조립자, 예를 들어 Flye28, Canu29, Raven30과 같은 최고의 조합기들조차도 광범위한 서열 상동성을 가진 플라스미드의 존재에 어려움을 겪거나, 잘못된 플라스미드 서열을 병합하거나, 단일 플라스미드의 여러 복제본을 생성하는 경우가 많으며, 작은 플라스미드58,59가 빠지는 경우가 많습니다(그림 4C 참조). 이럴 때 저품질 리드를 제거하고 짧은 리드를 제거하면 키메라 연결 문제를 해결할 수 있습니다.
시퀀싱 프로토콜의 전반적인 정확성을 확립하기 위해, 저자들은 ONT 전용 어셈블리와 롱리드와 숏리드를 결합해 다듬어진 어셈블리를 비교했는데, 이 어셈블리는 (PacBio 데이터를 제외하고) 현재 골드 스탠다드입니다. 저자들은 또한 플라스미드 수와 추정 크기를 보여주는 겔 데이터도 가지고 있었습니다(표 8, 세 번째 열). 연구진은 10개 샘플에서 발견된 13개의 플라스미드를 모두 성공적으로 조립했고, 조립된 서열과 일치하는 크기를 얻었습니다. 이는 매우 성공적인데, LR 서열 조립기는 구조적으로 복잡한 플라스미드를 조립하는 데 실패하는 경우가 많기 때문입니다. 더 세밀한 수준에서는 두 서열 집합 간의 쌍별 비교에서 정확도가 가변적임을 나타냈다(표 8). 플라스미드 조립체 중 다섯 개는 ONT 조립체에서 노출되지 않은 영역이 있었는데, 이는 ONT 전용 서열 조립 과정에서 추가된 서열 인공성 때문으로 추정됩니다. 나머지 7개 중 3개는 단일 염기체 다형성을 보였으며, 쌍별 차이는 0.0023%에서 0.05% 사이였고, 4개는 완벽한 일치를 보였다. 정확도 수준은 플라스미드 DNA 풍부도 수준이나 시퀀싱 깊이와 상관관계가 없으나, 한 샘플에 여러 플라스미드가 존재하면 상관관계가 있습니다. 이는 여러 플라스미드가 있을 경우 서열 복잡성이 증폭되기 때문일 수 있습니다(표 8). 요약하자면, 이 워크플로우는 플라스미드 서열 수율과 정확도를 최적화하도록 설계되었습니다. 플라스미드 서열을 조립할 수 있지만, 뉴클레오타이드 수준에서는 정확도가 변동적입니다. 이 프로토콜은 다양한 크기의 플라스미드를 시퀀싱할 수 있으며, 최소 72%의 염색체 DNA 오염을 견딜 수 있고, 플라스미드 크기와 플로우셀 상태에 따라 단일 플로우셀에 최소 24개의 플라스미드를 수용할 수 있습니다.
저자들은 공개할 것이 없습니다.
워싱턴 대학교의 스티브 살리판테가45년에 발표된 임상 대장균 균주를 아낌없이 공유해 주셔서 매우 감사드립니다. 이 연구는 UCSC에서 MC에 혁신 및 기업가 개발 이사직을 임명하기 위한 연구 지원 기금과, 캘리포니아 대학교 Alianza MX의 "멕시코 고부담 및 취약 인구 간 결핵 분리의 유전체 프로파일링: 전염 및 약물 내성에 대한 통찰"이라는 연구 지원 자금의 지원을 받았습니다.
| 이름 | 회사 | 카탈로그 번호 | 댓글 |
|---|---|---|---|
| 1.5 mL Eppendorf DNA LoBind tube | Sigma Aldrich | EP022431081 | 다른 브랜드로 교체 가능 |
| 1.5 mL tube racks | Thermo Fisher Scientific | 22-313630 | 다른 브랜드로 교체 가능 |
| 15 mL conical tubes | FALCON | 352096 | 다른 브랜드로 교체 가능 |
| 15 mL tube racks | Thermo Fisher Scientific | 8850 | 다른 브랜드로 교체 가능 |
| 2.0 mL microcentrifuge tubes | Fisherbrand | 05-408-138 | 다른 브랜드로 교체 가능 |
| 250 mL sterile flasks | PYREX | 5320 | 다른 브랜드로 교체 가능 |
| Carbenicillin | GOLD BIOTECHNOLOGY | C-103-50 | 다른 브랜드로 교체 가능 |
| Centrifuge | Eppendorf | 5424R-PR-R | 다른 브랜드로 교체 가능 |
| Centrifuge | Eppendorf | 5810R-P-R | 다른 브랜드로 교체 가능 |
| Centrifuge bottles | Sigma Aldrich | B1283-4EA | 다른 브랜드로 교체 가능 |
| Disposable Inoculating Loops: 1 μL | Fisherbrand | 22-363-595 | 다른 브랜드로 교체 가능 |
| Electronic Digital Scale | Denver Instrument | APX-2001 | 다른 브랜드로 교체 가능 |
| Ethanol | Fisher Scientific | 04-355-720 | 다른 브랜드로 교체 가능 |
| Ethanol-resistant markers | Sharpie | 37001; 37002; 37003 | 다른 브랜드로 교체 가능 |
| Flow cell (R10.4.1) | Oxford Nanopore Technologies | FLO-MIN114 | |
| Flow Cell Wash Kit | Oxford Nanopore Technologies | EXP-WSH004 | |
| Gloves | X-GEN | 44-100M | 모든 니트릴 장갑 브랜드 가능 |
| Ice bucket | Thermo Fisher Scientific | 432128 | 다른 브랜드로 교체 가능 |
| Isopropanol | Fisher Scientific | BP26184 | 다른 브랜드로 교체 가능 |
| MinION Sequencing Device | Oxford Nanopore Technologies | MIN-101B | |
| Native Barcoding Kit 96 V14 | Oxford Nanopore Technologies | SQK-NBD114.96 | |
| NucleoBond Xtra MIDI Plasmid DNA Extraction Kit | Takara | 740410.5 | |
| PCR tube racks | AXYGEN | R96PCRFSP | 다른 브랜드로 교체 가능 |
| PCR tubes: (0.2 mL) | AXYGEN | PCR-02-C | 다른 브랜드로 교체 가능 |
| Program autocycler v0.5.2 | Ryan Wick / University of Melbourne | ||
| Program bcftools v1.22 | Wellcome Sanger Institute (HTSlib / Samtools team) | ||
| Program bwa v0.7.19 | Heng Li, originally at the Broad Institute | ||
| Program canu v2.3 | University of California, Davis & Pacific Biosciences collaborators | ||
| Program flye v2.9.6 | University of California, San Diego (Pavel Pevzner lab) | ||
| Program medaka v2.1.1 | Oxford Nanopore Technologies | ||
| Program metamdbg v1.2 | INRIA | ||
| Program miniasm v0.3 | Heng Li | ||
| Program minimap2 v2.28 | Heng Li | ||
| Program minipolish v0.2.0 | Ryan Wick / University of Melbourne | ||
| Program porechop v0.2.4 | Ryan Wick / University of Melbourne | ||
| Program racon v1.5.0 | Genome Institute of Singapore | ||
| Program samtools v1.22.1 | Wellcome Sanger Institute (HTSlib / Samtools project) | ||
| Qubit® 4 Fluorometer | Invitrogen | Q33226 | |
| Qubit® Assay Tubes | Invitrogen | Q32856 | |
| Rapid Barcoding Kit 24 and 96 V14 | Oxford Nanopore Technologies | SQK-RBK114.96 | |
| Set of micropipettes that dispense between 1 – 10 μL (P10), 2 – 20 μL (P20), 20 – 200 μL (P200) and 200 – 1000 μL (P1000) | RAININ | 17008648; 17008650; 17008652; 17008653 | Micropipettes should be calibrated |
| Spectrophotometer | Thermo Scientific | 335905P | 다른 브랜드로 교체 가능 |
| Sterilized tips for 10 μL, 200 μL and 1000 μL | Eclipse | 1011-260-000-9; 1018-260-000; 1019-260-000-9 | |
| Thermal cycler C1000 Touch | Bio-Rad | 1851196 | 다른 브랜드로 교체 가능 |
| Water bath | Thermo Scientific | 51221052 | 다른 브랜드로 교체 가능 |