BOL: All site bookmarks

Bioinformatics in Thailand !

Shruti Paniwala — Wed, 28 Apr 2021 02:04:56 -0500

Our international PhD and master programs are designed for students who desire focused training in the elements of biology, computer science, and information technology needed for a successful career in the exciting new discipline of Bioinformatics & Systems Biology. Students in our program will receive comprehensive training in omics analysis, database design and management, software engineering and programming (including web-based development), simulation techniques and modeling, and data integration. Each student will apply their skills to a practical project, where they will design and implement a solution to a real-world problem under the guidance of an experienced mentor in industry or academia.

https://bioinformatics.kmutt.ac.th/about.html

Duangrudee Tanramluk (Ajarn Wi) uses computational biology and machine learning to tackle the key to drug design problems via MANORAA webserver.

https://mb.mahidol.ac.th/en/bioinformatics/

https://graduate.mahidol.ac.th/inter/

This international Doctorate programme is designed to further broaden students’ knowledge in Bioinformatics and Molecular Biology to their maximum capability.

http://www.mbb.psu.ac.th/programmes/phd

Ph.D. program in Bioinformatics and Computational Biology is a joint effort of the Faculty of Science and Faculty of Medicine, Chulalongkorn University. The program has study plans for both applicants who hold a bachelor’s degree and applicants who hold a master’s degree in any related fields of study.

http://www.bioinfo.sc.chula.ac.th/ph-d-program-specialization/

Additional detail

https://www.biotec.or.th/en/index.php/research/research-units/genome-technology-research-unit

https://tbrcnetwork.org/labtbrc/index.php/bioinformatics-and-chemoinformatics/

https://genomicsthailand.com/Genomic/home

Address of the bookmark: https://bioinformatics.kmutt.ac.th/

CoronaVIR: Computational Resources on Novel Coronavirus (SARS-CoV-2 or COVID-19)

LEGE — Tue, 27 Apr 2021 01:58:36 -0500

Aim of this web site is to facilitate the scientific community to fight against severe pandemic disease COVID-19 caused by SARS-CoV-2. Here, We have collected and organized information related to novel strain of coronavirus, i.e. SARS-CoV-2.and its resulting disease COVID-19 from the literature and other resources from the Internet. We are providing links to appropriate literature. Moreover, we are Bioinformatics Group, based on our knowledge and expertise, we are also proposing potential diagnostics primers, peptide and RNA based vaccine candidates and potential drug molecules. These are predicted candidates, need to be validated by experimental Researchers, who have appropriate infrastructure. It is an integrated multi-omics repository dedicated to current genomic, proteomic, diagnostic and therapeutic knowledge about coronaviruses particularly the recent strain, i.e. SARS-CoV-2 or 2019-nCoV. This web resource will be helpful for the researchers engaged in the development of therapies and drugs for the COVID-19. The information is collected from various available resources.
Cite: Patiyal, Sumeet, et al. “A Web-based Platform on COVID-19 to Maintain Predicted Diagnostic, Drug
and Vaccine Candidates.” Monoclon Antib Immunodiagn Immunother. doi.org/10.1089/mab.2020.0035

Address of the bookmark: https://webs.iiitd.edu.in/raghava/coronavir/

Modular, efficient and constant-memory single-cell RNA-seq preprocessing

Jit — Mon, 05 Apr 2021 11:19:43 -0500

With kallisto | bustools you can

Generate a cell x gene or cell x transcript equivalence class count matrix
Perform RNA velocity and single-nuclei RNA-seq analsis
Quantify data from numerous technologies such as 10x, inDrops, and Dropseq.
Customize workflows for new technologies and protocols.
Process feature barcoding data such as CITE-seq, REAP-seq, MULTI-seq, Clicktags, and Perturb-seq.
Obtain QC reports from single-cell RNA-seq data

The kallisto | bustools workflow is described in:

Páll Melsted*, A. Sina Booeshaghi*, Lauren Liu, Fan Gao, Lambda Lu, Kyung Hoi (Joseph) Min, Eduardo da Veiga Beltrame, Kristján Eldjárn Hjörleifsson, Jase Gehring & Lior Pachter† Modular and efficient pre-processing of single-cell RNA-seq, Nature Biotechnology (2021).

Documentation and tutorials for the kallisto bustools workflow are available at http://pachterlab.github.io/kallistobustools.

https://www.nature.com/articles/s41587-021-00870-2

Address of the bookmark: https://pachterlab.github.io/kallistobustools/

A simple tutorial for a complex ComplexHeatmap

Neel — Fri, 02 Apr 2021 06:18:32 -0500

ComplexHeatmap (Gu, Eils, and Schlesner (2016)) is an R Programming Language (R Core Team (2020)) package that is currently listed in the Bioconductor package repository.

install and load required packages

  require(RColorBrewer)
  require(ComplexHeatmap)
  require(circlize)
  require(digest)
  require(cluster)

If all load successfully, proceed to Part 3. Otherwise, go through the following code chunks in order to ensure that each package is installed and loaded properly.

BiocManager (Morgan (2019))

Address of the bookmark: https://github.com/kevinblighe/E-MTAB-6141

DEG 5.0: a database of essential genes in both prokaryotes and eukaryotes

Rahul Nayak — Tue, 30 Mar 2021 11:47:29 -0500

Essential genes are those indispensable for the survival of an organism, and their functions are therefore considered a foundation of life. Determination of a minimal gene set needed to sustain a life form, a fundamental question in biology, plays a key role in the emerging field, synthetic biology.

DEG is freely available at the website http://tubic.tju.edu.cn/deg or http://www.essentialgene.org.

Address of the bookmark: http://www.essentialgene.org/

DEG 5.0: a database of essential genes in both prokaryotes and eukaryotes

Rahul Nayak — Tue, 30 Mar 2021 11:47:28 -0500

DEG is freely available at the website http://tubic.tju.edu.cn/deg or http://www.essentialgene.org.

Address of the bookmark: http://www.essentialgene.org/

Public Databases for Bioinformatics !

Jit — Tue, 23 Mar 2021 05:32:15 -0500

https://www.nature.com/articles/s41467-020-17155-y

Server Infrastructure:

File Server:

dhara: Synology 3614 Storage Appliance
4 Core Xeon
108TB disk storage
10Gb ethernet to SCG3
Access atx: dhara:5000
Has btsync server (try it - its much better than dropbox)

Compute Servers:

nandi: Kundaje and Phi Server
24 intel cores
256GB RAM
500GB of SSD storage 
36TB RAID6 local storage
4 Intel Phi's (space for 4 more GPU's)


durga: Montgomery and sensitive data
24 intel cores
256GB RAM
500GB of SSD RAID0 storage 
60TB RAID6 local storage

mitra: Bassik and Web/DB Server
24 core
256GB RAM 
500GB of SSD RAID0 storage 
36TB RAID6 local storage

vayu: Kundaje GPU server
4 core
64GB RAM 
200GB of SSD storage 
8TB RAID10 local storage
4 Nvidia GTX 970 4GB GPUs

amold: Bickel and SGE server
32 AMD core
128GB RAM 
200GB of SSD storage 
12TB RAID5 local storage

wotan: Bickel and SGE server
64 AMD core
256GB RAM 
200GB of SSD storage 
12TB RAID5 local storage

Filesystem:

/users/$USER
default home directory
full backups nightly 
nfs mount to dhara
should store code, papers, and other highly processed data here

/mnt/data/
globally accessible data
should store common data here
e.g. genomes and indexes, annotations, ENCODE data  
if you dont want this to count towards your quote you must chown

/mnt/lab_data/$LAB/
lab accessible data
should store lab project data here 
e.g. ATAC-seq prediction data, enhancer prediction, motif calls

/srv/scratch/$USER
fast local storage
not backed up, but on raid and data will never be deleted
most analysis should be performed here

/srv/persistent/$USER
fast local storage
synced nightly, but not backed up
       ie if the hard drives fail or you delete something and notice 
       within 24 hours we can recover. Otherwise not. (vs home which is 
       properly backed up )  
intermediate analysis products that would be hard to recover should be stored here 
       e.g. stochastic analysis results that need to be kept so that paper 
       results can be reproduced

/srv/www/$LABNAME/
web accessible from mitra.stanford.edu
*NOT BACKED UP*

Some parallel programming patterns:

# gzip a bunch of files
parallel gzip -- *.FILESTOGZIP

# fork example in python:
(for more detailed examples look at 
 https://github.com/nboley/grit/ grit/lib/multiprocessing_utils.py)

import os
import time
import random

import multiprocessing

class ProcessSafeOPStream( object ):
    def __init__( self, writeable_obj ):
        self.writeable_obj = writeable_obj
        self.lock = multiprocessing.Lock()
        self.name = self.writeable_obj.name
        return
    
    def write( self, data ):
        self.lock.acquire()
        self.writeable_obj.write( data )
        self.writeable_obj.flush()
        self.lock.release()
        return
    
    def close( self ):
        self.writeable_obj.close()

def worker(queue, ofp):
    # Try without this
    random.seed()
    while True:
        i = queue.get()
        if i == 'FINISHED': return
        # simulate an expensive function
        x = random.random()
        time.sleep(x/10)
        print i, x
        ofp.write("%i\t%s\n" % (i, x))

NSIMS = 10000
NPROC = 25

# populate queue
todo = multiprocessing.Queue()
for i in xrange(NSIMS): todo.put(i)
for i in xrange(NPROC): todo.put('FINISHED')

ofp = ProcessSafeOPStream( open("output.txt", "w") )

pids = []
for i in xrange(NPROC):
    pid = os.fork()
    if pid == 0:
       worker(todo, ofp)
       os._exit(0)
    else:
       pids.append(pid)  

for pid in pids:
    os.waitpid(pid, 0)

ofp.close()

print "FINISHED"

For use case 1 we obtained the following ENCODE and ROADMAP datasets https://www.encodeproject.org/files/ENCFF446WOD/@@download/ENCFF446WOD.bed.gz, https://www.encodeproject.org/files/ENCFF546PJU/@@download/ENCFF546PJU.bam, https://www.encodeproject.org/files/ENCFF059BEU/@@download/ENCFF059BEU.bam. Blacklisted regions were obtained from http://mitra.stanford.edu/kundaje/akundaje/release/blacklists/hg38-human/hg38.blacklist.bed.gz. The human genome version hg38 was obtained from http://hgdownload.cse.ucsc.edu/goldenPath/hg38/bigZips/hg38.fa.gz.

For use case 2 we used the set of narrowPeak files summarized in https://github.com/wkopp/janggu_usecases/tree/master/extra/urls.txt (archived version v1.0.1). The human genome version hg19 was obtained from http://hgdownload.cse.ucsc.edu/goldenPath/hg19/bigZips/hg19.fa.gz

For use case 3 we used the ENCODE datasets https://www.encodeproject.org/files/ENCFF591XCX/@@download/ENCFF591XCX.bam, https://www.encodeproject.org/files/ENCFF736LHE/@@download/ENCFF736LHE.bigWig, https://www.encodeproject.org/files/ENCFF177HHM/@@download/ENCFF177HHM.bam as we as the GENCODE annotation v29 from ftp://ftp.ebi.ac.uk/pub/databases/gencode/Gencode_human/release_29/gencode.v29.annotation.gtf.gz.

Address of the bookmark: http://mitra.stanford.edu/

Janggu - Deep learning for Genomics

Jit — Tue, 23 Mar 2021 05:14:43 -0500

Janggu is a python package that facilitates deep learning in the context of genomics. The package is freely available under a GPL-3.0 license.

Detail tutorial at https://janggu.readthedocs.io/en/latest/

USE cases

https://github.com/wkopp/janggu_usecases

Address of the bookmark: https://github.com/BIMSBbioinfo/janggu

Nucl2Vec: Local alignment of DNA sequences using Distributed Vector Representation

Jit — Tue, 16 Mar 2021 05:45:44 -0500

We demonstrate a novel approach forlocal alignment of DNA reads with respect to reference genome.For this process we have used Skip-gram model for creatingencoding(Nucl2Vec) and k-nearest neighbor for the alignment.With our new approach we have reduced computation cost forlocal alignment , while achieving accuracy comparable to existingdefacto standard BWA-MEM tool.

https://prakharg24.github.io/papers/401851.full.pdf

Address of the bookmark: https://prakharg24.github.io/papers/401851.full.pdf

DAVI: Deep learning-based tool for alignment and single nucleotide variant identification

Jit — Tue, 16 Mar 2021 05:41:33 -0500

DAVI consists of models for both global and local alignment and for variant calling. We have evaluated the performance of DAVI against existing state-of-the-art tool sets and found that its accuracy and performance is comparable to existing tools used for bench-marking. We further demonstrate that while existing tools are based on data generated from a specific sequencing technology, the models proposed in DAVI are generic and can be used across different NGS technologies as well as across different species

https://iopscience.iop.org/article/10.1088/2632-2153/ab7e19/pdf

Address of the bookmark: https://github.com/gguptaiitd/NEAT