hic分析之分步工作原理

分步工作原理：

先看test文件

/PATH/TO/HiC-Pro_2.11.1/test-op/test_data
.
|-- dixon_2M
|   |-- SRR400264_00_R1.fastq.gz
|   `-- SRR400264_00_R2.fastq.gz
`-- dixon_2M_2
    |-- SRR400264_01_R1.fastq.gz
    `-- SRR400264_01_R2.fastq.gz

再看看opt文件

tree
.
|-- bowtie_results
|   |-- bwt2
|   |   |-- dixon_2M
|   |   |   |-- SRR400264_00_R1_hg19.bwt2merged.bam
|   |   |   |-- SRR400264_00_R1_hg19.mapstat
|   |   |   |-- SRR400264_00_R2_hg19.bwt2merged.bam
|   |   |   |-- SRR400264_00_R2_hg19.mapstat
|   |   |   |-- SRR400264_00_hg19.bwt2pairs.bam
|   |   |   `-- SRR400264_00_hg19.bwt2pairs.pairstat
|   |   `-- dixon_2M_2
|   |       |-- SRR400264_01_R1_hg19.bwt2merged.bam
|   |       |-- SRR400264_01_R1_hg19.mapstat
|   |       |-- SRR400264_01_R2_hg19.bwt2merged.bam
|   |       |-- SRR400264_01_R2_hg19.mapstat
|   |       |-- SRR400264_01_hg19.bwt2pairs.bam
|   |       |-- SRR400264_01_hg19.bwt2pairs.pairstat
|   |-- bwt2_global
|   |   |-- dixon_2M
|   |   |   |-- SRR400264_00_R1_hg19.bwt2glob.bam
|   |   |   |-- SRR400264_00_R1_hg19.bwt2glob.unmap.fastq
|   |   |   |-- SRR400264_00_R2_hg19.bwt2glob.bam
|   |   |   |-- SRR400264_00_R2_hg19.bwt2glob.unmap.fastq
|   |   |-- dixon_2M_2
|   |       |-- SRR400264_01_R1_hg19.bwt2glob.bam
|   |       |-- SRR400264_01_R1_hg19.bwt2glob.unmap.fastq
|   |       |-- SRR400264_01_R2_hg19.bwt2glob.bam
|   |       |-- SRR400264_01_R2_hg19.bwt2glob.unmap.fastq
|   |-- bwt2_local
|       |-- dixon_2M
|       |   |-- SRR400264_00_R1_hg19.bwt2glob.unmap_bwt2loc.bam
|       |   |-- SRR400264_00_R1_hg19.bwt2glob.unmap_trimmed.fastq
|       |   |-- SRR400264_00_R2_hg19.bwt2glob.unmap_bwt2loc.bam
|       |   |-- SRR400264_00_R2_hg19.bwt2glob.unmap_trimmed.fastq
|       |-- dixon_2M_2
|           |-- SRR400264_01_R1_hg19.bwt2glob.unmap_bwt2loc.bam
|           |-- SRR400264_01_R1_hg19.bwt2glob.unmap_trimmed.fastq
|           |-- SRR400264_01_R2_hg19.bwt2glob.unmap_bwt2loc.bam
|           |-- SRR400264_01_R2_hg19.bwt2glob.unmap_trimmed.fastq
|-- config_test_latest.txt
|-- hic_results
|   |-- data
|   |   |-- dixon_2M
|   |   |   |-- SRR400264_00_hg19.bwt2pairs.DEPairs
|   |   |   |-- SRR400264_00_hg19.bwt2pairs.DumpPairs
|   |   |   |-- SRR400264_00_hg19.bwt2pairs.FiltPairs
|   |   |   |-- SRR400264_00_hg19.bwt2pairs.REPairs
|   |   |   |-- SRR400264_00_hg19.bwt2pairs.RSstat
|   |   |   |-- SRR400264_00_hg19.bwt2pairs.SCPairs
|   |   |   |-- SRR400264_00_hg19.bwt2pairs.SinglePairs
|   |   |   |-- SRR400264_00_hg19.bwt2pairs.validPairs
|   |   |   |-- SRR400264_00_hg19.bwt2pairs_interaction.bam
|   |   |   `-- dixon_2M.allValidPairs
|   |   `-- dixon_2M_2
|   |       |-- SRR400264_01_hg19.bwt2pairs.DEPairs
|   |       |-- SRR400264_01_hg19.bwt2pairs.DumpPairs
|   |       |-- SRR400264_01_hg19.bwt2pairs.FiltPairs
|   |       |-- SRR400264_01_hg19.bwt2pairs.REPairs
|   |       |-- SRR400264_01_hg19.bwt2pairs.RSstat
|   |       |-- SRR400264_01_hg19.bwt2pairs.SCPairs
|   |       |-- SRR400264_01_hg19.bwt2pairs.SinglePairs
|   |       |-- SRR400264_01_hg19.bwt2pairs.validPairs
|   |       |-- SRR400264_01_hg19.bwt2pairs_interaction.bam
|   |       \-- dixon_2M_2.allValidPairs
|   |-- matrix
|   |   |-- dixon_2M
|   |   |   |-- iced
|   |   |   |   |-- 500000
|   |   |   |-- raw
|   |   |       |-- 1000000
|   |   |       |   |-- dixon_2M_1000000.matrix
|   |   |       |   |-- dixon_2M_1000000_abs.bed
|   |   |       |-- 500000
|   |   |           |-- dixon_2M_500000.matrix
|   |   |           |-- dixon_2M_500000_abs.bed
|   |   |-- dixon_2M_2
|   |       |-- raw
|   |           |-- 1000000
|   |           |   |-- dixon_2M_2_1000000.matrix
|   |           |   |-- dixon_2M_2_1000000_abs.bed
|   |           |-- 500000
|   |               |-- dixon_2M_2_500000.matrix
|   |               |-- dixon_2M_2_500000_abs.bed
|   |-- pic
|   |   |-- dixon_2M
|   |   |   |-- plotHiCContactRanges_dixon_2M.pdf
|   |   |   |-- plotHiCFragmentSize_dixon_2M.pdf
|   |   |   |-- plotHiCFragment_dixon_2M.pdf
|   |   |   |-- plotMappingPairing_dixon_2M.pdf
|   |   |   |-- plotMapping_dixon_2M.pdf
|   |   |-- dixon_2M_2
|   |       |-- plotHiCContactRanges_dixon_2M_2.pdf
|   |       |-- plotHiCFragmentSize_dixon_2M_2.pdf
|   |       |-- plotHiCFragment_dixon_2M_2.pdf
|   |       |-- plotMappingPairing_dixon_2M_2.pdf
|   |       |-- plotMapping_dixon_2M_2.pdf
|   |-- stats
|       |-- dixon_2M
|       |   |-- dixon_2M.mRSstat
|       |   |-- dixon_2M.mpairstat
|       |   |-- dixon_2M_R1.mmapstat
|       |   |-- dixon_2M_R2.mmapstat
|       |   |-- dixon_2M_allValidPairs.mergestat
|       |-- dixon_2M_2
|           |-- dixon_2M_2.mRSstat
|           |-- dixon_2M_2.mpairstat
|           |-- dixon_2M_2_R1.mmapstat
|           |-- dixon_2M_2_R2.mmapstat
|           |-- dixon_2M_2_allValidPairs.mergestat
|-- logs
|   |-- dixon_2M
|   |   |-- SRR400264_00_R1_bowtie2.log
|   |   |-- SRR400264_00_R1_hg19.bwt2glob.unmap_bowtie2.log
|   |   |-- SRR400264_00_R1_hg19.bwt2glob.unmap_readsTrimming.log
|   |   |-- SRR400264_00_R2_bowtie2.log
|   |   |-- SRR400264_00_R2_hg19.bwt2glob.unmap_bowtie2.log
|   |   |-- SRR400264_00_R2_hg19.bwt2glob.unmap_readsTrimming.log
|   |   |-- build_raw_maps.log
|   |   |-- ice_500000.log
|   |   |-- make_Rplots.log
|   |   |-- mapped_2hic_fragments.log
|   |   |-- mapping_combine.log
|   |   |-- mapping_stats.log
|   |   |-- mapping_step1.log
|   |   |-- mapping_step2.log
|   |   |-- mergeSAM.log
|   |   |-- merge_stats.log
|   |   |-- merge_valid_interactions.log
|   |   |-- plot_hic_contacts.Rout
|   |   |-- plot_hic_fragment.Rout
|   |   |-- plot_mapping_portion.Rout
|   |   |-- plot_pairing_portion.Rout
|   |-- dixon_2M_2
|       |-- SRR400264_01_R1_bowtie2.log
|       |-- SRR400264_01_R1_hg19.bwt2glob.unmap_bowtie2.log
|       |-- SRR400264_01_R1_hg19.bwt2glob.unmap_readsTrimming.log
|       |-- SRR400264_01_R2_bowtie2.log
|       |-- SRR400264_01_R2_hg19.bwt2glob.unmap_bowtie2.log
|       |-- SRR400264_01_R2_hg19.bwt2glob.unmap_readsTrimming.log
|       |-- build_raw_maps.log
|       |-- make_Rplots.log
|       |-- mapped_2hic_fragments.log
|       |-- mapping_combine.log
|       |-- mapping_stats.log
|       |-- mapping_step1.log
|       |-- mapping_step2.log
|       |-- mergeSAM.log
|       |-- merge_stats.log
|       |-- merge_valid_interactions.log
|       |-- plot_hic_contacts.Rout
|       |-- plot_hic_fragment.Rout
|       |-- plot_mapping_portion.Rout
|       |-- plot_pairing_portion.Rout
|-- rawdata -> /PATH/TO/HiC-Pro_2.11.1/test-op/test_data
|-- tmp

工作原理介绍：

image.png

1.Reads Mapping

将fastq文件用bowtie2 mapping一次，看log：

cat mapping_step1.log

###mapping fastq to reference genome , 将unmapping的read 另存为 unmap.fastq
/PATH/TO/bowtie2/bowtie2-2.3.3.1/bowtie2 --very-sensitive -L 30 --score-min L,-0.6,-0.2 --end-to-end --reorder --un bowtie_results/bwt2_global/dixon_2M/SRR400264_00_R2_hg19.bwt2glob.unmap.fastq --rg-id BMG --rg SM:SRR400264_00_R2 -p 20 -x /PATH/TO/HiC-Pro_2.11.1/Bowtie2Index/hg19 -U rawdata/dixon_2M/SRR400264_00_R2.fastq.gz 2>> logs/dixon_2M/SRR400264_00_R2_bowtie2.log| /PATH/TO/samtools/samtools-1.6/bin/samtools view -F 4 -bS - > bowtie_results/bwt2_global/dixon_2M/SRR400264_00_R2_hg19.bwt2glob.bam
/PATH/TO/bowtie2/bowtie2-2.3.3.1/bowtie2 --very-sensitive -L 30 --score-min L,-0.6,-0.2 --end-to-end --reorder --un bowtie_results/bwt2_global/dixon_2M/SRR400264_00_R1_hg19.bwt2glob.unmap.fastq --rg-id BMG --rg SM:SRR400264_00_R1 -p 20 -x /PATH/TO/HiC-Pro_2.11.1/Bowtie2Index/hg19 -U rawdata/dixon_2M/SRR400264_00_R1.fastq.gz 2>> logs/dixon_2M/SRR400264_00_R1_bowtie2.log| /PATH/TO/samtools/samtools-1.6/bin/samtools view -F 4 -bS - > bowtie_results/bwt2_global/dixon_2M/SRR400264_00_R1_hg19.bwt2glob.bam

得到这样的文件：

|   |-- bwt2_global
|   |   |-- dixon_2M
|   |   |   |-- SRR400264_00_R1_hg19.bwt2glob.bam
|   |   |   |-- SRR400264_00_R1_hg19.bwt2glob.unmap.fastq
|   |   |   |-- SRR400264_00_R2_hg19.bwt2glob.bam
|   |   |   |-- SRR400264_00_R2_hg19.bwt2glob.unmap.fastq
|   |   `-- dixon_2M_2
|   |       |-- SRR400264_01_R1_hg19.bwt2glob.bam
|   |       |-- SRR400264_01_R1_hg19.bwt2glob.unmap.fastq
|   |       |-- SRR400264_01_R2_hg19.bwt2glob.bam
|   |       |-- SRR400264_01_R2_hg19.bwt2glob.unmap.fastq

这里有mapping的，和unmapping的。因为read是嵌合基因。mapping上的是背景基因组。unmapping的是需要下一步处理的。

2.片段分配和过滤

先看log

cat mapping_step2.log
### 酶切
/PATH/TO/HiC-Pro/HiC-Pro_2.11.1/scripts/cutsite_trimming --fastq bowtie_results/bwt2_global/dixon_2M/SRR400264_00_R1_hg19.bwt2glob.unmap.fastq --cutsite AAGCTAGCTT --out bowtie_results/bwt2_local/dixon_2M/SRR400264_00_R1_hg19.bwt2glob.unmap_trimmed.fastq > logs/dixon_2M/SRR400264_00_R1_hg19.bwt2glob.unmap_readsTrimming.log 2>&1
/PATH/TO/HiC-Pro/HiC-Pro_2.11.1/scripts/cutsite_trimming --fastq bowtie_results/bwt2_global/dixon_2M/SRR400264_00_R2_hg19.bwt2glob.unmap.fastq --cutsite AAGCTAGCTT --out bowtie_results/bwt2_local/dixon_2M/SRR400264_00_R2_hg19.bwt2glob.unmap_trimmed.fastq > logs/dixon_2M/SRR400264_00_R2_hg19.bwt2glob.unmap_readsTrimming.log 2>&1
### 再用bowtie2 mapping一次
/PATH/TO/bowtie2/bowtie2-2.3.3.1/bowtie2 --very-sensitive -L 20 --score-min L,-0.6,-0.2 --end-to-end --reorder --rg-id BML --rg SM:SRR400264_00_R1_hg19.bwt2glob.unmap -p 20 -x /PATH/TO/HiC-Pro/HiC-Pro_2.11.1/Bowtie2Index/hg19 -U bowtie_results/bwt2_local/dixon_2M/SRR400264_00_R1_hg19.bwt2glob.unmap_trimmed.fastq 2>> logs/dixon_2M/SRR400264_00_R1_hg19.bwt2glob.unmap_bowtie2.log | /PATH/TO/samtools/samtools-1.6/bin/samtools view -bS - > bowtie_results/bwt2_local/dixon_2M/SRR400264_00_R1_hg19.bwt2glob.unmap_bwt2loc.bam
/PATH/TO/bowtie2/bowtie2-2.3.3.1/bowtie2 --very-sensitive -L 20 --score-min L,-0.6,-0.2 --end-to-end --reorder --rg-id BML --rg SM:SRR400264_00_R2_hg19.bwt2glob.unmap -p 20 -x /PATH/TO/HiC-Pro/HiC-Pro_2.11.1/Bowtie2Index/hg19 -U bowtie_results/bwt2_local/dixon_2M/SRR400264_00_R2_hg19.bwt2glob.unmap_trimmed.fastq 2>> logs/dixon_2M/SRR400264_00_R2_hg19.bwt2glob.unmap_bowtie2.log | /PATH/TO/samtools/samtools-1.6/bin/samtools view -bS - > bowtie_results/bwt2_local/dixon_2M/SRR400264_00_R2_hg19.bwt2glob.unmap_bwt2loc.bam

这一步是将unmapping的结果，用酶切位点分割，然后将分割下来的结果，再mapping一次。得到以下结果

|   |-- bwt2_local
|       |-- dixon_2M
|       |   |-- SRR400264_00_R1_hg19.bwt2glob.unmap_bwt2loc.bam
|       |   |-- SRR400264_00_R1_hg19.bwt2glob.unmap_trimmed.fastq
|       |   |-- SRR400264_00_R2_hg19.bwt2glob.unmap_bwt2loc.bam
|       |   |-- SRR400264_00_R2_hg19.bwt2glob.unmap_trimmed.fastq
|       |-- dixon_2M_2
|           |-- SRR400264_01_R1_hg19.bwt2glob.unmap_bwt2loc.bam
|           |-- SRR400264_01_R1_hg19.bwt2glob.unmap_trimmed.fastq
|           |-- SRR400264_01_R2_hg19.bwt2glob.unmap_bwt2loc.bam
|           |-- SRR400264_01_R2_hg19.bwt2glob.unmap_trimmed.fastq

3、mapping_combine.log

cat mapping_combine.log

### 将global的bam与local的bam用samtools merge
/PATH/TO/samtools/samtools-1.6/bin/samtools merge -@ 20 -n -f bowtie_results/bwt2/dixon_2M/SRR400264_00_R1_hg19.bwt2merged.bam bowtie_results/bwt2_global/dixon_2M/SRR400264_00_R1_hg19.bwt2glob.bam bowtie_results/bwt2_local/dixon_2M/SRR400264_00_R1_hg19.bwt2glob.unmap_bwt2loc.bam
/PATH/TO/samtools/samtools-1.6/bin/samtools merge -@ 20 -n -f bowtie_results/bwt2/dixon_2M/SRR400264_00_R2_hg19.bwt2merged.bam bowtie_results/bwt2_global/dixon_2M/SRR400264_00_R2_hg19.bwt2glob.bam bowtie_results/bwt2_local/dixon_2M/SRR400264_00_R2_hg19.bwt2glob.unmap_bwt2loc.bam

### sort
/PATH/TO/samtools/samtools-1.6/bin/samtools sort -@ 20 -n -T tmp/SRR400264_00_R2_hg19 -o bowtie_results/bwt2/dixon_2M/SRR400264_00_R2_hg19.bwt2merged.sorted.bam bowtie_results/bwt2/dixon_2M/SRR400264_00_R2_hg19.bwt2merged.bam
/PATH/TO/samtools/samtools-1.6/bin/samtools sort -@ 20 -n -T tmp/SRR400264_00_R1_hg19 -o bowtie_results/bwt2/dixon_2M/SRR400264_00_R1_hg19.bwt2merged.sorted.bam bowtie_results/bwt2/dixon_2M/SRR400264_00_R1_hg19.bwt2merged.bam

[bam_sort_core] merging from 0 files and 20 in-memory blocks...
[bam_sort_core] merging from 0 files and 20 in-memory blocks...

### 改名
mv bowtie_results/bwt2/dixon_2M/SRR400264_00_R1_hg19.bwt2merged.sorted.bam bowtie_results/bwt2/dixon_2M/SRR400264_00_R1_hg19.bwt2merged.bam
mv bowtie_results/bwt2/dixon_2M/SRR400264_00_R2_hg19.bwt2merged.sorted.bam bowtie_results/bwt2/dixon_2M/SRR400264_00_R2_hg19.bwt2merged.bam

4、merge

### 将R1与R2 merge为一个文件。
cat mergeSAM.log
/PATH/TO/python/python-2.7.10/bin/python /PATH/TO/HiC-Pro/HiC-Pro_2.11.1/scripts/mergeSAM.py -q 0 -t -v -f bowtie_results/bwt2/dixon_2M/SRR400264_00_R1_hg19.bwt2merged.bam -r bowtie_results/bwt2/dixon_2M/SRR400264_00_R2_hg19.bwt2merged.bam -o bowtie_results/bwt2/dixon_2M/SRR400264_00_hg19.bwt2pairs.bam
## mergeBAM.py
## forward= bowtie_results/bwt2/dixon_2M/SRR400264_00_R1_hg19.bwt2merged.bam
## reverse= bowtie_results/bwt2/dixon_2M/SRR400264_00_R2_hg19.bwt2merged.bam
## output= bowtie_results/bwt2/dixon_2M/SRR400264_00_hg19.bwt2pairs.bam
## min mapq= 0
## report_single= False
## report_multi= False
## verbose= True
## Merging forward and reverse tags ...

5、利用HiCPro的mapped_2hic_fragments.py程序将比对结果转化为Hi-C片段信息

cat mapped_2hic_fragments.log
/PATH/TO/python /PATH/TO/HiC-Pro/HiC-Pro_2.11.1/scripts/mapped_2hic_fragments.py -v -S -t 100 -m 100000 -s 100 -l 600 -a -f /PATH/TO/HiC-Pro_2.11.1/annotation/HindIII_resfrag_hg19.bed -r bowtie_results/bwt2/dixon_2M/SRR400264_00_hg19.bwt2pairs.bam -o hic_results/data/dixon_2M
## overlapMapped2HiCFragments.py
## mappedReadsFile= bowtie_results/bwt2/dixon_2M/SRR400264_00_hg19.bwt2pairs.bam
## fragmentFile= /PATH/TO/HiC-Pro/HiC-Pro_2.11.1/annotation/HindIII_resfrag_hg19.bed
## minInsertSize= 100
## maxInsertSize= 600
## minFragSize= 100
## maxFragSize= 100000
## allOuput= True
## SAM ouput= True
## verbose= True

再对输出的valid pairs文件进行排序：

LANG=en; sort -T tmp -k2,2V -k3,3n -k5,5V -k6,6n -o hic_results/data/sample/sample_sample_genome_ref.bwt2pairs.validPairs hic_results/data/sample/sample_sample_genome_ref.bwt2pairs.validPairs

6 对所有的valid pairs进行合并

LANG=en; sort -T tmp -S 50% -k2,2V -k3,3n -k5,5V -k6,6n -m hic_results/data/dixon_2M_2/SRR400264_01_hg19.bwt2pairs.validPairs | awk -F"\t" 'BEGIN{c1=0;c2=0;s1=0;s2=0}(c1!=$2 || c2!=$5 || s1!=$3 || s2!=$6){print;c1=$2;c2=$5;s1=$3;s2=$6}' > hic_results/data/dixon_2M_2/dixon_2M_2.allValidPairs

这样就成了我们最终的文件，dixon_2M_2.allValidPairs。

这个文件还不能可视化，如果要可视化，一般在jucebox里面查看。我们需要用配套的hicpro2juicebox.sh来转化：

## e.g.1
HICPRO_PATH/bin/utils/hicpro2juicebox.sh -i hicpro_res/hic_results/data/dixon_2M/dixon_2M_allValidPairs -g hg19 -j /usr/local/juicebox/juicebox_clt_1.4.jar

很多人问，这个-j的 jar文件在哪里？
事实上用这个就可以。

/PATH/TO/juicer/juicer/UGER/scripts/juicer_tools.jar

下次，我们再讨论juicer与juicebox的用法。

最后编辑于：2019.08.26 11:02:29

人面猴
序言：七十年代末，一起剥皮案震惊了整个滨河市，随后出现的几起案子，更是在滨河造成了极大的恐慌，老刑警刘岩，带你破解...
沈念sama阅读 194,088评论 5赞 459
死咒
序言：滨河连续发生了三起死亡事件，死亡现场离奇诡异，居然都是意外死亡，警方通过查阅死者的电脑和手机，发现死者居然都...
沈念sama阅读 81,715评论 2赞 371
救了他两次的神仙让他今天三更去死
文/潘晓璐我一进店门，熙熙楼的掌柜王于贵愁眉苦脸地迎上来，“玉大人，你说我怎么就摊上这事。” “怎么了？”我有些...
开封第一讲书人阅读 141,361评论 0赞 319
道士缉凶录：失踪的卖姜人
文/不坏的土叔我叫张陵，是天一观的道长。经常有香客问我，道长，这世上最难降的妖魔是什么？我笑而不...
开封第一讲书人阅读 52,099评论 1赞 263
港岛之恋（遗憾婚礼）
正文为了忘掉前任，我火速办了婚礼，结果婚礼上，老公的妹妹穿的比我还像新娘。我一直安慰自己，他们只是感情好，可当我...
茶点故事阅读 60,987评论 4赞 355
恶毒庶女顶嫁案：这布局不是一般人想出来的
文/花漫我一把揭开白布。她就那样静静地躺着，像睡着了一般。火红的嫁衣衬着肌肤如雪。梳的纹丝不乱的头发上，一...
开封第一讲书人阅读 46,063评论 1赞 272
城市分裂传说
那天，我揣着相机与录音，去河边找鬼。笑死，一个胖子当着我的面吹牛，可吹牛的内容都是我干的。我是一名探鬼主播，决...
沈念sama阅读 36,486评论 3赞 381
双鸳鸯连环套：你想象不到人心有多黑
文/苍兰香墨我猛地睁开眼，长吁一口气：“原来是场噩梦啊……” “哼！你这毒妇竟也来了？” 一声冷哼从身侧响起，我...
开封第一讲书人阅读 35,175评论 0赞 253
万荣杀人案实录
序言：老挝万荣一对情侣失踪，失踪者是张志新（化名）和其女友刘颖，没想到半个月后，有当地人在树林里发现了一具尸体，经...
沈念sama阅读 39,440评论 1赞 290
护林员之死
正文独居荒郊野岭守林人离奇死亡，尸身上长有42处带血的脓包…… 初始之章·张勋以下内容为张勋视角年9月15日...
茶点故事阅读 34,518评论 2赞 309
白月光启示录
正文我和宋清朗相恋三年，在试婚纱的时候发现自己被绿了。大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
茶点故事阅读 36,305评论 1赞 326
活死人
序言：一个原本活蹦乱跳的男人离奇死亡，死状恐怖，灵堂内的尸体忽然破棺而出，到底是诈尸还是另有隐情，我是刑警宁泽，带...
沈念sama阅读 32,190评论 3赞 312
日本核电站爆炸内幕
正文年R本政府宣布，位于F岛的核电站，受9级特大地震影响，放射性物质发生泄漏。R本人自食恶果不足惜，却给世界环境...
茶点故事阅读 37,550评论 3赞 298
男人毒药：我在死后第九天来索命
文/蒙蒙一、第九天我趴在偏房一处隐蔽的房顶上张望。院中可真热闹，春花似锦、人声如沸。这庄子的主人今日做“春日...
开封第一讲书人阅读 28,880评论 0赞 17
一桩弑父案，背后竟有这般阴谋
文/苍兰香墨我抬头看了看天上的太阳。三九已至，却和暖如春，着一层夹袄步出监牢的瞬间，已是汗流浃背。一阵脚步声响...
开封第一讲书人阅读 30,152评论 1赞 250
情欲美人皮
我被黑心中介骗来泰国打工，没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留，地道东北人。一个月前我还...
沈念sama阅读 41,451评论 2赞 341
代替公主和亲
正文我出身青楼，却偏偏与公主长得像，于是被迫代替她去往敌国和亲。传闻我的和亲对象是个残疾皇子，可洞房花烛夜当晚...
茶点故事阅读 40,637评论 2赞 335