Tcbf

❤️ Tcbf = Topologically associating domain(TAD) Conservative Boundary Finder

中文文档

📣 Introduction

TADs are fundamental regulatory chromatin structures and are largely conserved across tissues and species. We developed a Python/R pipeline Tcbf to identify the conserved TAD boundaries between multiple genome.

✨ Pre-requisite:

MCL

The mcl clustering algorithm is available in the repositories of some Linux distributions and so can be installed in the same way as any other package. For example, on Ubuntu, Debian, Linux Mint:

sudo apt-get install mcl

Alternatively, it can be built from source which will likely require the 'build-essential' or equivalent package on the Linux distribution being used. Instructions are provided on the MCL webpage, http://micans.org/mcl/.

Minimap2

The default aligner is minimap2. In the future, we will consider supporting more aligner.

Mash

Mash is used to measure the genetics distance between different genome sequences and adjust the alignment parameters.

R language and the following packages

install.packages("BiocManager")
BiocManager::install("GenomicRanges")
BioManager::install("plyranges")
install.packages("tidyverse")

Python >= 3.6

This package is heavily dependent on f-string. As of Python 3.6, f-strings are a great new way to format strings. Not only are they more readable, more concise, and less prone to error than other ways of formatting, but they are also faster!

🔰 Installation

Install from source

git clone https://github.com/TcbfGroup/Tcbf
cd Tcbf
pip install -r requirements.txt
python setup.py install

//: # () //: # ()

//: # ()

//: # () //: # ()

//: # ()

📝 Usage

You can test the Tcbf pipeline with the example.

cd example;
bash download.sh;
Tcbf.py run -c config.txt  -o test

Advanced

For HPC users have multiple nodes, we also provide a step by step to accelerate work.

tcbf run -c config.txt  -o test --only_print_command

For more details, see [Here]()

Inputs

The config.txt is a table file with four columns. The first columns is path of genome sequence path in FASTA format.

The second columns is TAD annotaed file, which generated from HiTAD.

The third column is the individual name.

the fouth columns in gene annotation file in gff3 format

Of note, pay attention to the contigs or scaffolds. For organisms have well assembly, we usually recommand to remove all scaffolds.

genome1.fasta   genome1_tad.txt g1  genome1.gff3
genome2.fasta   genome2_tad.txt g2  genome2.gff3
genome3.fasta   genome3_tad.txt g3  genome3.gff3
genome4.fasta   genome4_tad.txt g4  genome4.gff3

we use the results from HiTAD as the TAD annotate file

#cat genome1_tad.txt
Chr01   0       2300000
Chr01   2300000 3950000
Chr01   3950000 4600000
Chr01   4600000 5900000
Chr01   5900000 6350000

Output

For conserved TAD boundaries among species, we provide a table file to show the the clustering results. This table has one TAD boundary group per line and one spcies per column and is ordered from largest orthogroup to smallest.

Visualization

The user can input a region of genome, and the software can automatically obtain the corresponding three-dimensional structural relationship. The red part in the figure represents the TAD range, and the yellow rectangle represents the TAD boundary. Light blue indicates collinear pairs.

tcbf plot-syn-pair -o test --reference human --chrom chr7 --start 127910000 --end 131410000 --plot example.pdf human macaque mouse

For heatmap visualization of multiple species, please refer to Here or online document.

Advanced features

The good performance of Tcbf depends on a few conditions. (1) It is challenging to accurately detect the conserved TAD boundaries in highly fragmented genomes with low-quality gene annotations. To minimize bias, we recommend using genomes with removal of unplaced contigs and scaffolds and high gene annotation quality. With the rapid advances in sequencing technologies for the assembly of high-quality reference genomes, this should not be an issue for most species. (2) In the TAD boundary identification process, we recommend using high-quality Hi-C or similar data with the same resolution in different species to reduce the input data bias of TAD coordinates. (3) For input data with variable genome sizes, users may select an appropriate TAD boundary range through the average TAD size to improve the quality of the results.

😉 Author

Tbcf are maintained by: * @HeXin

For more Help, Please leave a message in the issue, I will reply as soon as possible.

📃 License

MIT