PCR 与分子克隆原理（PCR & Molecular Cloning）¶

┌─────────────────────────────────────────────────────────────┐
│  第1步：变性（Denaturation）  94-98°C，20-30秒              │
│  白话：高温把 DNA 双链"拉开"成两条单链                      │
│  类比：把拉链拉开                                           │
├─────────────────────────────────────────────────────────────┤
│  第2步：退火（Annealing）     50-65°C，20-30秒              │
│  白话：降温让引物"贴"到目标位置                              │
│  类比：书签精准插到对应页码                                  │
├─────────────────────────────────────────────────────────────┤
│  第3步：延伸（Extension）     72°C，30秒-数分钟              │
│  白话：DNA 聚合酶从引物开始，一个碱基一个碱基地"抄写"新链   │
│  类比：复印机从书签位置开始复印                              │
└─────────────────────────────────────────────────────────────┘

            重复 25-35 个循环
         ↓
    2^30 ≈ 10 亿份拷贝（理论值）

# 简单公式（适用于 <20bp 短引物）：
Tm = 2(A+T) + 4(G+C)    # 单位：°C

# 更准确的最近邻法（Nearest-Neighbor Method）：
Tm = ΔH / (ΔS + R × ln(C/4)) - 273.15
# ΔH = 焓变之和, ΔS = 熵变之和, R = 气体常数, C = 引物浓度

# Primer3 命令行版（适合批量设计）
# 安装（conda）
conda install -c bioconda primer3-py  # Python接口

# Python 调用示例
import primer3

# 设计引物
result = primer3.bindings.design_primers(
    seq_args={
        'SEQUENCE_ID': 'my_gene',
        'SEQUENCE_TEMPLATE': 'ATCGATCGATCG...（你的目标序列）',
        'SEQUENCE_INCLUDED_REGION': [50, 200],  # 在50-250bp范围内找引物
    },
    global_args={
        'PRIMER_OPT_SIZE': 20,        # 最优引物长度
        'PRIMER_MIN_SIZE': 18,        # 最小长度
        'PRIMER_MAX_SIZE': 25,        # 最大长度
        'PRIMER_OPT_TM': 60.0,       # 最优Tm值
        'PRIMER_MIN_TM': 55.0,       # 最小Tm值
        'PRIMER_MAX_TM': 65.0,       # 最大Tm值
        'PRIMER_MIN_GC': 40.0,       # 最小GC含量
        'PRIMER_MAX_GC': 60.0,       # 最大GC含量
        'PRIMER_PRODUCT_SIZE_RANGE': [[100, 300]],  # 产物大小范围
    }
)

# 输出结果
print(f"左引物: {result['PRIMER_LEFT_0_SEQUENCE']}")  # 正向引物序列
print(f"右引物: {result['PRIMER_RIGHT_0_SEQUENCE']}") # 反向引物序列
print(f"产物大小: {result['PRIMER_PAIR_0_PRODUCT_SIZE']} bp")
print(f"左引物Tm: {result['PRIMER_LEFT_0_TM']:.1f}°C")
print(f"右引物Tm: {result['PRIMER_RIGHT_0_TM']:.1f}°C")

使用流程：
1. 输入目标序列的 Accession Number 或粘贴序列
2. 指定 PCR 产物范围（如 Forward primer from: 100, to: 200）
3. 选择目标物种（Organism）
4. 设置引物参数（Tm、长度等）
5. 点击 "Get Primers" → 系统自动设计 + 特异性验证
6. 结果中 "Products on target templates" = 特异性好

┌──────────────────────────────────────────────────────────┐
│  Step 1: 酶切（Digestion）                               │
│  用限制性内切酶把载体和目标片段各切出匹配的"接口"         │
│  类比：把文件夹剪开一个口 + 给文件切出配套的形状          │
├──────────────────────────────────────────────────────────┤
│  Step 2: 连接（Ligation）                                │
│  用 DNA 连接酶（T4 Ligase）把片段和载体"粘"在一起        │
│  类比：用胶水把文件粘进文件夹                             │
├──────────────────────────────────────────────────────────┤
│  Step 3: 转化（Transformation）                          │
│  把重组质粒送进感受态细菌（通常是大肠杆菌 DH5α）         │
│  方法：热激法（42°C, 45秒）或电转化                      │
│  类比：把装好文件的文件夹交给复印店                       │
├──────────────────────────────────────────────────────────┤
│  Step 4: 筛选（Screening）                               │
│  用抗性平板筛选 + 菌落PCR/酶切验证/测序确认               │
│  类比：检查复印店有没有正确复印你的文件                   │
└──────────────────────────────────────────────────────────┘

        ┌──── 复制起点（ori）：决定质粒能在哪种宿主中复制
        │
        │     ┌── 抗性基因（如 AmpR）：筛选标记，只有含质粒的细菌能存活
        │     │
   ─────┼─────┼──────────────────────────
  │     │     │                          │
  │   [ori]  [AmpR]                      │
  │                                      │
  │         ┌── 多克隆位点（MCS）         │
  │         │   一排限制酶切割位点         │
  │        [MCS]                         │
  │                                      │
  │   [Promoter]──[Insert]──[Terminator] │
  │         │                            │
   ─────────┼────────────────────────────
             └── 目标基因插入区域

# 使用 UCSC In-Silico PCR（网页版）
# 网址：https://genome.ucsc.edu/cgi-bin/hgPcr
# 输入正反引物序列 → 选择基因组版本 → 查看预测扩增子

# 命令行工具：isPcr（Jim Kent 开发）
# 安装
conda install -c bioconda isPcr  # 或从 UCSC 下载二进制

# 运行 isPcr
isPcr genome.2bit primers.txt output.fa \
    -minPerfect=15 \    # 3'端完美匹配最少碱基数
    -minGood=15 \       # 总体最少匹配碱基数
    -maxSize=4000       # 最大产物大小

# primers.txt 格式（制表符分隔）：
# 引物名    正向引物序列    反向引物序列
# gene1    ATCGATCGATCG    GCTAGCTAGCTA

# 用 BLAST 验证引物特异性
blastn -query primer.fa -db nt \
    -task blastn-short \       # 短序列专用模式
    -evalue 1000 \             # 放宽阈值（引物短，需要宽松阈值）
    -word_size 7 \             # 短序列用小word_size
    -outfmt 6 \                # 表格输出
    -num_threads 4

# 或用 Bowtie2 快速比对引物到参考基因组
bowtie2 -x ref_genome -f primers.fa -S primers_aligned.sam \
    --very-sensitive \         # 高灵敏度模式
    -k 10                      # 报告最多10个比对位置（检查脱靶）

# QIIME2 扩增子分析核心流程
# 1. 导入数据
qiime tools import \
    --type 'SampleData[PairedEndSequencesWithQuality]' \
    --input-path manifest.tsv \
    --output-path demux.qza \
    --input-format PairedEndFastqManifestPhred33V2

# 2. 引物去除（cutadapt）
qiime cutadapt trim-paired \
    --i-demultiplexed-sequences demux.qza \
    --p-front-f CCTACGGGNGGCWGCAG \   # 341F 引物序列
    --p-front-r GACTACHVGGGTATCTAATCC \ # 805R 引物序列
    --o-trimmed-sequences trimmed.qza

# 3. DADA2 去噪（生成 ASV）
qiime dada2 denoise-paired \
    --i-demultiplexed-seqs trimmed.qza \
    --p-trunc-len-f 280 \     # 正向读长截断位置
    --p-trunc-len-r 200 \     # 反向读长截断位置
    --o-table table.qza \     # ASV 丰度表
    --o-representative-sequences rep-seqs.qza  # ASV 代表序列

PCR 扩增倍数 = 2^n（n = 循环数，理论值）
Tm ≈ 2(A+T) + 4(G+C)（简易公式，短引物）
退火温度 ≈ Tm - 5°C
Taq 延伸速度 ≈ 1 kb/min
Pfu/Q5 延伸速度 ≈ 0.5 kb/min（更慢但更准）
Taq 错误率 ≈ 1×10⁻⁴ / bp / cycle
Q5 错误率 ≈ 5×10⁻⁷ / bp / cycle