Exploring genetic variation in insecticide resistance loci in Anopheles gambiae s.l. populations from Ghana

問題の背景
マラリアは依然としてサブサハラアフリカにおける重大な公衆衛生上の負担であり、ガーナは世界の症例数に大きく寄与している。媒介蚊対策は、殺虫剤処理済みの蚊帳（ITN）や室内残留噴霧（IRS）に大きく依存しているが、その有効性は、Anopheles gambiae sensu lato (s.l.) 集団における広範な殺虫剤抵抗性によって、ますます脅かされている。vgsc、rdl、ace1における標的部位変異や、gste2における代謝性抵抗性といった抵抗性メカニズムは既知であるが、ガーナの異なる生態学的ゾーンにおける系統、遺伝的多様性、および抵抗性アレル（対立遺伝子）の具体的な分布に関しては、依然として重大な地理的空白が存在する。さらに、An. gambiae s.l. 種の形態学的同定は不十分な場合が多く、精密なサーベイランスには分子生物学的なツールが必要となる。

研究手法
本研究では、Oxford Nanopore Technologies (ONT) を用いた、低コストかつ高スループットなマルチプレックス・アンプリコン・シーケンシング手法を採用した。

• サンプリング: ガーナの3つの主要な生態学的ゾーン（沿岸サバンナ、森林、ギニアサバンナ）の18地点から、330個体のAn. gambiae s.l. 蚊を採集した。一部の蚊については、WHOチューブテストを用いたピレスロイドおよび有機リン系殺虫剤に対する感受性の表現型解析が事前に行われていた。
• シーケンシング戦略: 個々の蚊からDNAを抽出した。4つのドメインを持つ電位依存性ナトリウムチャネル（vgsc）、ace1、rdl、gste2の部分配列、および種同定マーカー（ITS1、ITS2、SINE200、cox1、nd4、28S、IGS）の計7つのゲノム領域を増幅するようにマルチプレックスPCRアッセイを設計した。
• バーコーディングとシーケンシング: サンプルのプールを可能にするため、プライマーに固有の8bpバーコードを付加した。ライブラリはNative Barcoding Kit 24 V14を用いて調製し、MinION Mk1Bデバイスとr10.4.1フローセルを使用してシーケンシングを行った。
• バイオインフォマティクス: リードのベースコール、デマルチプレックス、およびAn. gambiae参照ゲノム（AgamP4）へのアライメントを行った。変異の呼び出しにはFreebayesおよびGATK HaplotypeCallerを用い、品質によるフィルタリングとアノテーションを行った。種同定はSINE200挿入の有無に基づいて行った。系統解析および多様性解析は、ITS2およびミトコンドリア配列を用いて実施した。

主な貢献と結果

• 種の分布: 分子生物学的同定により、沿岸域（77.6%）および森林域（86.5%）ではAn. coluzziiが優占種であることが明らかになった一方、ギニアサバンナ域ではAn. gambiae s.s. が最も多く見られた（50.0%）。ハイブリッドも全ゾーンで検出された。
• 遺伝的多様性と遺伝子流動: ITS2およびミトコンドリア配列の系統解析の結果、核塩基多様性は低く、地理的または種による明確なクラスタリングは見られなかった。これは、ガーナのAn. gambië s.l. 集団間において高い遺伝子流動があることを示している。
• 抵抗性アレルの検出: 対象領域において202個のSNPを特定し、そのうち31個の非同義置換を検出した。
  • 既知の変異: 既知の抵抗性関連変異を13個検出した：vgscにおいて7個（高頻度のL995Fが75%、L995Sが6%を含む）、rdlにおいて2個（A296Gが51%）、ace1において1個（G280Sが22%）、およびgste2において3個（I114T、F120L、L119V）。
  • 新規および新興のバリアント: 17個の新規SNPを特定した。特に、ギニアビサウで最近観察された変異であるgste2のT154Sが24%の頻度で検出された。その他の新規変異には、vgscのA1553TおよびL1680Fが含まれる。
  • 連鎖不平衡: vgscの変異V402LとI1527Tの間には強い連鎖（r² = 0.90）が観察され、これらが頻繁に同一のハプロタイプ上に存在することを示唆している。
• 種特異的な変異: アレル頻度において種間の有意な差が認められた。例えば、vgsc L995FはAn. gambiae s.s.（90%）においてAn. coluzzii（72%）よりも高かった。また、vgsc V402LおよびI1527TはAn. coluzziiにおいて有意に高かった。
• 表現型と遺伝型の関連: 表現型解析を行った限定的なサブセットにおいて、ほとんどのSNPは抵抗性表現型との有意な関連を示さなかったが、vgsc L995SおよびI1527Tのヘテロ接合型は、それぞれデルタメトリンおよびペルメトリン抵抗性と有意に関連していた。

意義
本論文は、ナノポア技術を用いた高スループットなマルチプレックス・アンプリコン・シーケンシングが、マラリア媒介蚊のゲノム監視のためのスケーラブルでコスト効率の高い手法であることを実証している。本研究は、複雑な既知および新規のバリアントを明らかにし、抵抗性の進化のダイナミックな性質を浮き彫りにすることで、ガーナにおける殺虫剤抵抗性の遺伝学的景観を特徴づけることに成功した。これらの知見は、抵抗性アレルの迅速な拡散を促進する可能性のある高い遺伝子流動を考慮すると、標的を絞った殺虫剤抵抗性管理および適応的な媒介蚊対策を情報提供するために、継続的なゲノムモニタリングが必要であることを強調している。