Clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 gene editing technology has been widely used to facilitate efficient genome editing. Current popular sgRNA design tools only consider the sgRNA perfectly matched to the target site and provide the results without any on-target mismatch. We suppose taking on-target gRNA-DNA mismatches into consideration might provide better sgRNA with similar binding activity and reduced off-target sites. Here, we trained a seq2seq-attention model with feedback-loop architecture, to automatically generate sgRNAs with on-target mismatches. Dual-luciferase reporter experiment showed that multiple sgRNAs with three mismatches could achieve the 80% of the relative activity of the perfect matched sgRNA. Meanwhile, it could reduce the number of off-target sites using sgRNAs with on-target mismatches. Finally, we provided a freely accessible web server sgRNA design tool named ExsgRNA. Users could submit their target sequence to this server and get optimal sgRNAs with less off-targets and similar on-target activity compared with the perfect-matched sgRNA.

The triangular correlation heatmap aiming to visualize the linkage disequilibrium (LD) pattern and haplotype block structure of SNPs is ubiquitous component of population-based genetic studies. However, current tools suffered from the problem of time and memory consuming. Here, we developed LDBlockShow, an open source software, for visualizing LD and haplotype blocks from variant call format files. It is time and memory saving. In a test dataset with 100 SNPs from 60 000 subjects, it was at least 10.60 times faster and used only 0.03-13.33% of physical memory as compared with other tools. In addition, it could generate figures that simultaneously display additional statistical context (e.g. association P-values) and genomic region annotations. It can also compress the SVG files with a large number of SNPs and support subgroup analysis. This fast and convenient tool will facilitate the visualization of LD and haplotype blocks for geneticists.

Bioinformatics (2020) doi: 10.1093/bioinformatics/btaa567 In the originally published version of this manuscript, author affiliation 1 was incorrect. The incorrect affiliation was as follows: “Key Laboratory of Biomedical Information Engineering of Ministry of Education, Biomedical Informatics & Genomics Center, School of Life Science and Technology, Zhongnan Hospital of Wuhan University, Wuhan 430071, China.” The correct affiliation should read: “Key Laboratory of Biomedical Information Engineering of Ministry of Education, Biomedical Informatics & Genomics Center, School of Life Science and Technology, Xi'an Jiaotong University, Xi'an 710049, P. R. China.” This error has been corrected online.

MOTIVATION: CircRNAs are an abundant class of non-coding RNAs with widespread, cell-/tissue-specific patterns. Previous work suggested that epigenetic features might be related to circRNA expression. However, the contribution of epigenetic changes to circRNA expression has not been investigated systematically. Here, we built a machine learning framework named CIRCScan, to predict circRNA expression in various cell lines based on the sequence and epigenetic features. RESULTS: The predicted accuracy of the expression status models was high with area under the curve of receiver operating characteristic (ROC) values of 0.89-0.92 and the false-positive rates of 0.17-0.25. Predicted expressed circRNAs were further validated by RNA-seq data. The performance of expression-level prediction models was also good with normalized root-mean-square errors of 0.28-0.30 and Pearson's correlation coefficient r over 0.4 in all cell lines, along with Spearman's correlation coefficient ρ of 0.33-0.46. Noteworthy, H3K79me2 was highly ranked in modeling both circRNA expression status and levels across different cells. Further analysis in additional nine cell lines demonstrated a significant enrichment of H3K79me2 in circRNA flanking intron regions, supporting the potential involvement of H3K79me2 in circRNA expression regulation. AVAILABILITY AND IMPLEMENTATION: The CIRCScan assembler is freely available online for academic use at https://github.com/johnlcd/CIRCScan. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Genome-wide association studies (GWASs) are an effective strategy to identify susceptibility loci for human complex diseases. However, missing heritability is still a big problem. Most GWASs single-nucleotide polymorphisms (SNPs) are located in noncoding regions, which has been considered to be the unexplored territory of the genome. Recently, data from the Encyclopedia of DNA Elements (ENCODE) and Roadmap Epigenomics projects have shown that many GWASs SNPs in the noncoding regions fall within regulatory elements. In this study, we developed a pipeline named functional disease-associated SNPs prediction (FDSP), to identify novel susceptibility loci for complex diseases based on the interpretation of the functional features for known disease-associated variants with machine learning. We applied our pipeline to predict novel susceptibility SNPs for type 2 diabetes (T2D) and hypertension. The predicted SNPs could explain heritability beyond that explained by GWAS-associated SNPs. Functional annotation by expression quantitative trait loci analyses showed that the target genes of the predicted SNPs were significantly enriched in T2D or hypertension-related pathways in multiple tissues. Our results suggest that combining GWASs and regulatory features data could identify additional functional susceptibility SNPs for complex diseases. We hope FDSP could help to identify novel susceptibility loci for complex diseases and solve the missing heritability problem.

MOTIVATION: Linkage disequilibrium (LD) decay is of great interest in population genetic studies. However, no tool is available now to do LD decay analysis from variant call format (VCF) files directly. In addition, generation of pair-wise LD measurements for whole genome SNPs usually resulting in large storage wasting files. RESULTS: We developed PopLDdecay, an open source software, for LD decay analysis from VCF files. It is fast and is able to handle large number of variants from sequencing data. It is also storage saving by avoiding exporting pair-wise results of LD measurements. Subgroup analyses are also supported. AVAILABILITY AND IMPLEMENTATION: PopLDdecay is freely available at https://github.com/BGI-shenzhen/PopLDdecay.