The BLAST programs are widely used tools for searching protein and DNA databases for sequence similarities. For protein comparisons, a variety of definitional, algorithmic and statistical refinements described here permits the execution time of the BLAST programs to be decreased substantially while enhancing their sensitivity to weak similarities. A new criterion for triggering the extension of word hits, combined with a new heuristic for generating gapped alignments, yields a gapped BLAST program that runs at approximately three times the speed of the original. In addition, a method is introduced for automatically combining statistically significant alignments produced by BLAST into a position-specific score matrix, and searching the database using this matrix. The resulting Position-Specific Iterated BLAST (PSIBLAST) program runs at approximately the same speed per iteration as gapped BLAST, but in many cases is much more sensitive to weak but biologically relevant sequence similarities. PSI-BLAST is used to uncover several new and interesting members of the BRCT superfamily.

/pdf/gapped-blast-and-psi-blast-a-new-generation-of-protein-2ybl67x5yx.pdf

Gapped BLAST and PSI-BLAST: a new generation of protein database search programs.

The sensitivity of the commonly used progressive multiple sequence alignment method has been greatly improved for the alignment of divergent protein sequences. Firstly, individual weights are assigned to each sequence in a partial alignment in order to down-weight near-duplicate sequences and up-weight the most divergent ones. Secondly, amino acid substitution matrices are varied at different alignment stages according to the divergence of the sequences to be aligned. Thirdly, residue-specific gap penalties and locally reduced gap penalties in hydrophilic regions encourage new gaps in potential loop regions rather than regular secondary structure. Fourthly, positions in early alignments where gaps have been opened receive locally reduced gap penalties to encourage the opening up of new gaps at these positions. These modifications are incorporated into a new program, CLUSTAL W which is freely available.

/pdf/clustal-w-improving-the-sensitivity-of-progressive-multiple-2na2s9sq2h.pdf

Clustal w: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice

CLUSTAL X is a new windows interface for the widely-used progressive multiple sequence alignment program CLUSTAL W. The new system is easy to use, providing an integrated system for performing multiple sequence and profile alignments and analysing the results. CLUSTAL X displays the sequence alignment in a window on the screen. A versatile sequence colouring scheme allows the user to highlight conserved features in the alignment. Pull-down menus provide all the options required for traditional multiple sequence and profile alignment. New features include: the ability to cut-and-paste sequences to change the order of the alignment, selection of a subset of the sequences to be realigned, and selection of a sub-range of the alignment to be realigned and inserted back into the original alignment. Alignment quality analysis can be performed and low-scoring segments or exceptional residues can be highlighted. Quality analysis and realignment of selected residue ranges provide the user with a powerful tool to improve and refine difficult alignments and to trap errors in input sequences. CLUSTAL X has been compiled on SUN Solaris, IRIX5.3 on Silicon Graphics, Digital UNIX on DECstations, Microsoft Windows (32 bit) for PCs, Linux ELF for x86 PCs, and Macintosh PowerMac.

/pdf/the-clustal-x-windows-interface-flexible-strategies-for-3qez0eyvqt.pdf

The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools.

We report a major update of the MAFFT multiple sequence alignment program. This version has several new features, including options for adding unaligned sequences into an existing alignment, adjustment of direction in nucleotide alignment, constrained alignment and parallel processing, which were implemented after the previous major update. This report shows actual examples to explain how these features work, alone and in combination. Some examples incorrectly aligned by MAFFT are also shown to clarify its limitations. We discuss how to avoid misalignments, and our ongoing efforts to overcome such limitations.

https://www.science-open.com/document_file/00d79e6c-98e4-425c-9fd1-3a85d5fdefbc/PubMedCentral/00d79e6c-98e4-425c-9fd1-3a85d5fdefbc.pdf

MAFFT Multiple Sequence Alignment Software Version 7: Improvements in Performance and Usability

抗原变异可使得多种致病微生物易于逃避宿主免疫应答。表达在感染红细胞表面的恶性疟原虫红细胞表面蛋白1（PfPMP1）与感染红细胞、内皮细胞、树突状细胞以及胎盘的单个或多个受体作用，在黏附及免疫逃避中起关键的作用。每个单倍体基因组var基因家族编码约60种成员，通过启动转录不同的var基因变异体为抗原变异提供了分子基础。

疟原虫var基因转换速率变化导致抗原变异[英]／Paul H, Robert P, Christodoulou Z, et al//Proc Natl Acad Sci U S A

Chromatin accessibility mapping is a powerful approach to identify potential regulatory elements. In the popular ATAC-seq method, Tn5 transposase inserts sequencing adapters into accessible DNA ("tagmentation"). CUT&Tag is a tagmentation-based epigenomic profiling method in which antibody tethering of Tn5 to a chromatin epitope of interest profiles specific chromatin features in small samples and single cells. Here we show that by simply modifying the tagmentation conditions for histone H3K4me2/3 CUT&Tag, antibody-tethered tagmentation of accessible DNA sites is redirected to produce accessible DNA maps that are indistinguishable from the best ATAC-seq maps. Thus, DNA accessibility maps can be produced in parallel with CUT&Tag maps of other epitopes with all steps from nuclei to amplified sequencing-ready libraries performed in single PCR tubes in the laboratory or on a home workbench. As H3K4 methylation is produced by transcription at promoters and enhancers, our method identifies transcription-coupled accessible regulatory sites.

/pdf/efficient-transcription-coupled-chromatin-accessibility-583jbvcozr.pdf

Efficient transcription-coupled chromatin accessibility mapping in situ

Chromatin accessibility mapping is a powerful approach to identify potential regulatory elements. A popular example is ATAC-seq, whereby Tn5 transposase inserts sequencing adapters into accessible DNA ("tagmentation"). CUT&Tag is a tagmentation-based epigenomic profiling method in which antibody tethering of Tn5 to a chromatin epitope of interest profiles specific chromatin features in small samples and single cells. Here we show that by simply modifying the tagmentation conditions for histone H3K4me2 or H3K4me3 CUT&Tag, antibody-tethered tagmentation of accessible DNA sites is redirected to produce chromatin accessibility maps that are indistinguishable from the best ATAC-seq maps. Thus, chromatin accessibility maps can be produced in parallel with CUT&Tag maps of other epitopes with all steps from nuclei to amplified sequencing-ready libraries performed in single PCR tubes in the laboratory or on a home workbench. As H3K4 methylation is produced by transcription at promoters and enhancers, our method identifies transcription-coupled accessible regulatory sites.

/pdf/efficient-chromatin-accessibility-mapping-in-situ-by-1gxq17c5gp.pdf

Efficient chromatin accessibility mapping in situ by nucleosome-tethered tagmentation

Lysine 27-to-methionine (K27M) mutations in the H3.1 or H3.3 histone genes are characteristic of pediatric diffuse midline gliomas (DMGs). These oncohistone mutations dominantly inhibit histone H3K27 trimethylation and silencing, but it is unknown if the oncohistone type affects gliomagenesis. We show that the genomic distributions of H3.1 and H3.3 oncohistones in human patient-derived DMG cells are consistent with the DNA replication-coupled deposition of histone H3.1 and the predominant replication-independent deposition of histone H3.3. Although H3K27 trimethylation is reduced for both oncohistone types, H3.3K27M-bearing cells retain some domains, and only H3.1K27M-bearing cells lack H3K27 trimethylation. Neither oncohistone interferes with PRC2 binding. Using Drosophila as a model, we demonstrate that inhibition of H3K27 trimethylation occurs only when H3K27M oncohistones are deposited into chromatin and only when expressed in cycling cells. We propose that oncohistones inhibit the H3K27 methyltransferase as chromatin patterns are being duplicated in proliferating cells, predisposing them to tumorigenesis.

https://www.biorxiv.org/content/biorxiv/early/2020/05/20/2020.05.18.101709.full.pdf

Histone deposition pathways determine the chromatin landscapes of H3.1 and H3.3 K27M oncohistones

We recently described CUT&Tag, a general strategy for epigenomic profiing in which antibody-tethered Tn5 transposase integrates DNA sequencing adapters at sites of specific chromatin protein binding or histone modification in intact cells or nuclei. Here we introduce a simplified CUT&Tag method that can be performed at home to help ameliorate the interruption of bench research caused by COVID-19 physical distancing requirements. All steps beginning with frozen nuclei are performed in single PCR tubes through to barcoded library amplication and clean-up, ready for pooling and DNA sequencing. Our CUT&Tag@home protocol has minimal equipment, reagent and supply needs and does not require handling of toxic or biologically active materials. We show that data quality and reproducibility for samples down to ~100 nuclei compare favorably to datasets produced using lab-based CUT&Tag and other chromatin profiling methods. We use CUT&Tag@home with antibodies to trimethylated histone H3 lysine-4, −36, −27 and-9 to comprehensively profile the epigenome of human K562 cells, consisting respectively of active gene regulatory elements, transcribed gene bodies, developmentally silenced domains and constitutively silenced parasitic elements.

/pdf/profiling-the-epigenome-at-home-w55o9qwdod.pdf

Profiling the epigenome at home

For more than a century, Formalin Fixed Paraffin Embedded (FFPE) sample preparation has been the preferred method for long-term preservation of biological material. However, the use of FFPE samples for epigenomic studies has been difficult because of chromatin damage from long exposure to high concentrations of formaldehyde. Previously, we introduced Cleavage Under Targeted Accessible Chromatin (CUTAC), an antibody-targeted chromatin accessibility mapping protocol based on CUT&Tag. Here we show that simple modifications of our single-tube CUTAC protocol are sufficient to produce high-resolution maps of paused RNA Polymerase II (RNAPII) at enhancers and promoters using FFPE samples. We find that transcriptional regulatory element differences produced by FFPE-CUTAC distinguish between mouse brain tumor specimens and identify regulatory element markers with high confidence and precision. Our simple work-flow is suitable for automation, making possible affordable epigenomic profiling of archived biological samples for biomarker identification, clinical applications and retrospective studies.

/pdf/epigenomic-analysis-of-formalin-fixed-paraffin-embedded-25c6qs3u.pdf

Jorja G. Henikoff

Papers

Efficient transcription-coupled chromatin accessibility mapping in situ

Efficient chromatin accessibility mapping in situ by nucleosome-tethered tagmentation

Histone deposition pathways determine the chromatin landscapes of H3.1 and H3.3 K27M oncohistones

Profiling the epigenome at home

Epigenomic analysis of Formalin-Fixed Paraffin-Embedded samples by CUT&Tag