BOL: Related items

Pollard Lab

Fri, 25 Sep 2020 20:20:50 -0500

We are a bioinformatics research lab focused on developing novel methods and using them to study genome evolution, organization, and regulation. Our mission is to decode biomedical knowledge that is missed without rigorous statistical approaches.

http://docpollard.org/

Tools

http://docpollard.org/resources/software/

Ancient whole genome duplication (WGD) detection tools !

Rahul Nayak — Sun, 07 Mar 2021 00:32:44 -0600

There are two methods for ancient WGD detection, one is collinearity analysis, and the other is based on the Ks distribution map. Among them, Ks is defined as the average number of synonymous substitutions at each synonymous site, and there is also a Ka corresponding to it, which refers to the average number of non-synonymous substitutions at each non-synonymous site.

At present, some people have posted articles about the analysis process of WGD. I searched for the keyword "wgd pipeline" and found the following:

GenoDup: https:// github.com/MaoYafei/GenoDup-Pipeline
https://peerj.com/articles/6303/
WGDdetector: https:// github.com/yongzhiyang2 012/WGDdetector
https://bmcbioinformatics.biomedcentral.com/articles/10.1186/s12859-019-2670-3
wgd: https:// github.com/arzwa/wgd
https://bmcbioinformatics.biomedcentral.com/articles/10.1186/s12859-016-1142-2#Sec1
https://bmcbiol.biomedcentral.com/articles/10.1186/s12915-017-0399-x
GeNoGAP https://bmcbioinformatics.biomedcentral.com/articles/10.1186/s12859-016-1142-2
https://bmcbiol.biomedcentral.com/articles/10.1186/s12915-017-0399-x
https://github.com/dfguan/purge_dups
https://www.biorxiv.org/content/10.1101/2020.01.24.917997v1

This article introduces the usage of wgd.

Wgd cannot be installed directly with bioconda at present, so it is a little troublesome to install, because it depends on a lot of software. wgd depends on the following software

BLAST
MCL
MUSCLE/MAFFT/PRANK
PAML
PhyML/FastTree
i-ADHoRe

But the good news is that most of the software it depends on can be installed with bioconda

conda create -n wgd python=3.5 blast mcl muscle mafft prank paml fasttree cmake libpng mpi=1.0=mpich
conda activate wgd

Here mpi=1.0=mpich is selected, because i-adhore depends on mpich. If openmpi is installed, an error will appear while loading shared libraries: libmpi_cxx.so.40: cannot open shared object file: No such file or directory

After that, the installation is much simpler

git clone https://github.com/arzwa/wgd.git
cd wgd
pip install .
pip install git+https://github.com/arzwa/wgd.git
For i-ADHoRe, you need to register at http:// bioinformatics.psb.ugent.be /webtools/i-adhore/licensing/Agree to the license to download i-ADHoRe-3.0

Since my miniconda3 installed ~/opt/, the installation path is so~/opt/miniconda3/envs/wgd/

tar -zxvf i-adhore-3.0.01.tar.gz
cd i-adhore-3.0.01
mkdir -p build && cd build
cmake .. -DCMAKE_INSTALL_PREFIX=~/opt/miniconda3/envs/wgd/
make -j 4
make insatall

Take the sugarcane genome Saccharum spontaneum L as an example. The genome is 8-ploid with 32 chromosomes (2n = 4x8 = 32)

Download the tutorial for CDS and GFF annotation files

mkdir -p wgd_tutorial && cd wgd_tutorial
wget http://www.life.illinois.edu/ming/downloads/Spontaneum_genome/Sspon.v20190103.cds.fasta.gz
wget http://www.life.illinois.edu/ming/downloads/Spontaneum_genome/Sspon.v20190103.gff3.gz
gunzip *.gz

First conda activate wgdstart our analysis environment, and then start the analysis

Step 1 : Use to wgd mclidentify homologous genes in the genome

wgd mcl -n 20 --cds --mcl -s Sspon.v20190103.cds.fasta -o Sspon_cds.out

Step 2 : Use to wgd ksdbuild Ks distribution

wgd ksd --n_threads 80 Sspon_cds.out/Sspon.v20190103.cds.fasta.blast.tsv.mcl Sspon.v20190103.cds.fasta

Step 3 : If the quality of the genome is good, then wgd syncollinearity analysis can be used . It can help us find the collinearity block in the genome and the corresponding anchor point

wgd syn --feature gene --gene_attribute ID \
-ks wgd_ksd/Sspon.v20190103.cds.fasta.ks.tsv \
Sspon.v20190103.gff3 Sspon_cds.out/Sspon.v20190103.cds.fasta.blast.tsv.mcl

For more reading - There are 9 sub-modules in WGD

kde: KDE fitting to the Ks distribution
ksd: Ks distribution construction
mcl: BLASP comparison of All-vs-ALl + MCL classification analysis.
mix: Hybrid modeling of Ks distribution.
pre: preprocess the CDS file
syn: Call I-ADHoRe 3.0 to use GFF files for collinearity analysis
viz: draw histogram and density plot
wf1: Ks standard analysis procedure of the whole genome paranome (paranome), call mcl, ksd and syn
wf2: Ks standard analysis procedure of one-vs-one homologous gene (ortholog), call wcl and kSD

Calling variants in non-diploid systems

Neel — Sat, 26 Jun 2021 15:37:49 -0500

The main challenge associated with non-diploid variant calling is the difficulty in distinguishing between the sequencing noise (abundant in all NGS platforms) and true low frequency variants. Some of the early attempts to do this well have been accomplished on human mitochondrial DNA although the same approaches will work equally good on viral and bacterial genomes (Rebolledo-Jaramillo et al. 2014, Li et al. 2015).

Address of the bookmark: https://training.galaxyproject.org/training-material/topics/variant-analysis/tutorials/non-dip/tutorial.html

MitoZ: a toolkit for animal mitochondrial genome assembly, annotation and visualization

Neel — Tue, 30 Nov 2021 23:23:57 -0600

MitoZ, consisting of independent modules of de novo assembly, findMitoScaf (find Mitochondrial Scaffolds), annotation and visualization, that can generate mitogenome assembly together with annotation and visualization results from HTS raw reads.

https://academic.oup.com/nar/article/47/11/e63/5377471

Address of the bookmark: https://github.com/linzhi2013/MitoZ

chromeister: An ultra fast, heuristic approach to detect conserved signals in extremely large pairwise genome comparisons.

Jit — Thu, 03 Feb 2022 04:01:55 -0600

chromeister: An ultra fast, heuristic approach to detect conserved signals in extremely large pairwise genome comparisons.

USAGE:

-query: sequence A in fasta format
-db: sequence B in fasta format
-out: output matrix
-kmer Integer: k>1 (default 32) Use 32 for chromosomes and genomes and 16 for small bacteria
-diffuse Integer: z>0 (default 4) Use 4 for everything - if using large plant genomes you can try using 1
-dimension Size of the output matrix and plot. Integer: d>0 (default 1000) Use 1000 for everything that is not full genome size, where 2000 is recommended

Address of the bookmark: https://github.com/estebanpw/chromeister

JBrowse 2: a modular genome browser with views of synteny and structural variation

Abhi — Tue, 25 Apr 2023 20:58:52 -0500

igvjs - a create-react-app with igv package from npm installed. the igv.js is instrumented to output "DONE" to the console when finished, and to have an increased fetchSizeLimit (which is otherwise git in CRAM longread tests)
jb2-web - stock instance of jbrowse-web v1.7.5
jb1 - stock instance of jbrowse 1 v1.16.11
jb2 embedded - a create-react-app with @jbrowse/react-linear-genome-view

Address of the bookmark: https://github.com/GMOD/jb2profile

bacLIFE: an automated genome mining tool for identification of lifestyle associated genes

BioStar — Fri, 15 Mar 2024 04:59:14 -0500

bacLIFE is a streamlined computational workflow that annotates bacterial genomes and performs large-scale comparative genomics to predict bacterial lifestyles and to pinpoint candidate genes, denominated lifestyle-associated genes (LAGs), and biosynthetic gene clusters associated with each lifestyle detected. This whole process is divided into different modules:

Clustering module Predicts, clusters and annotates the genes of every input genome
Lifestyle prediction Employs a machine learning model to forecast bacterial lifestyle or other specified metadata
Analitical module (Shiny app) Results from the previous modules are embedded in a user-friendly interface for comprehensive and interactive comparative genomics.

You can find the complete wiki here [https://github.com/Carrion-lab/bacLIFE/wiki/bacLIFE-wiki]

Address of the bookmark: https://github.com/Carrion-lab/bacLIFE

The Day We Began Reading the DNA of the World

Jitendra Narayan — Sat, 05 Sep 2026 01:43:21 -0500

Picture waking up and knowing that in one lab, a scientist is reading the DNA of a whale. In another, a team is sequencing a rare plant. Meanwhile, researchers in India are decoding the genomes of different human populations.

Many organisms. Many countries. All share one remarkable goal: to understand the genetic story of life.

DNA is nature's instruction manual, written with just four letters: A, T, C, and G. For years, scientists could only read small sections of this huge book. Now, new sequencing technology lets us read entire genomes like never before.

One of the most ambitious projects is the Earth BioGenome Project (EBP). Its goal is to create reference genomes for about 1.5 million known eukaryotic species. This project links genome research across continents, species, and scientific fields.

But this story goes beyond just animals and plants.

A New Chapter in Human Genomics

Large human genome projects are changing how we understand our own genetic diversity.

India's GenomeIndia project has sequenced thousands of people from different populations. This work is uncovering genetic variation that global databases have often missed.

Similar population-scale effSimilar large-scale projects are happening worldwide. The All of Us Research Program in the United States, Europe's 1+ Million Genomes initiative, South Korea's national genomic data program, and projects in Saudi Arabia, Australia, and Singapore are all building huge genomic resources.collect DNA.

The real goal is to link genomes with health, disease, and other biological information. This creates a base for more accurate research and, in time, more personalized medicine.

A Race Against Time

There is also an important twist.

We are sequencing Earth's biodiversity even as many species face growing environmental threats.

A genome alone cannot save a species, but it does keep important information about its biology, evolution, and genetic diversity. This knowledge helps scientists understand risks, guide conservation, and find traits that could help agriculture or medicine.

Projects like the Darwin Tree of Life, which studies species in Britain and Ireland, and the African BioGenome Project, which is growing genomics work across Africa, are key parts of this worldwide effort.

Projects to Watch

Earth BioGenome Project (EBP) - A global effort to sequence and annotate roughly 1.5 million known eukaryotic species and create a comprehensive reference library of Earth's biodiversity.
Vertebrate Genomes Project (VGP) - Focuses on producing high-quality reference genomes for vertebrate species, providing a major foundation for the wider Tree of Life.
Darwin Tree of Life - A UK and Ireland initiative aiming to sequence approximately 70,000 species, creating a detailed genomic map of regional biodiversity.
African BioGenome Project (AfricaBP) - A pan-African initiative focused on sequencing Africa's biodiversity while strengthening genomics and bioinformatics capacity across the continent.
European Reference Genome Atlas (ERGA) - A European effort to generate high-quality reference genomes for Europe's biodiversity and connect national sequencing efforts.
Canada BioGenome Project - Building reference genomes for Canadian biodiversity, supporting conservation, research and the study of ecosystems.
California Conservation Genomics Project - Uses genomics to understand and protect California's biodiversity and help inform conservation decisions.
Bat1K - An ambitious global effort to sequence the genomes of all known bat species, helping scientists study their evolution, longevity, immunity and unique biology.
Bird 10,000 (B10K) - A global project aiming to create genome resources covering the world's bird diversity and reconstruct avian evolutionary history.
Global Invertebrate Genomics Alliance (GIGA) - Builds genomic resources for the enormous and genetically diverse world of invertebrates.
GenomeIndia - India's national human-genome initiative, designed to capture the country's extraordinary population diversity and create a reference resource for Indian genomics.
All of Us Research Program - A US national research program combining genomic information with health and lifestyle data at population scale.
1+ Million Genomes / Genome of Europe - A European effort to enable secure, cross-border use of genomic and health data and develop genome-scale resources involving more than one million people.
Korea National Integrated Bio Big Data Project - South Korea's large-scale national effort combining genomic and health information from a very large population cohort.
Saudi Genome Program - A national genomics initiative focused on genetic diseases, population genomics and precision medicine in Saudi Arabia.
Australian Genomics / Genomics Health Futures Mission - Australia's growing investment in genomic medicine, rare disease research and population-scale health genomics.

The Library of Life

The fuThe future of genomics is not just about sequencing one genome at a time. Now, it is about creating a worldwide network of genomes—human and non-human, individual and species—linked by powerful databases and advanced computational tools.ine a library without shelves.

In this library, the books are genomes. The pages are made of DNA. Scientists and machines are the readers. And this library is still being written. Each genome sequenced today adds a new page to the story of life and helps us better understand our place in it.

33rd Annual Convention of Indian Association for Cancer Research from 13th to 15th February 2014

Tue, 27 Aug 2013 10:37:08 -0500

RGCB is organizing the 33rd Annual Convention of Indian Association for Cancer Research from 13th to 15th February 2014 with the theme "Discovery, Innovation and Translation in Cancer Research"

Kindly log on to conference website http://rgcb.res.in/IACR2014 for further details and timely updates and registration. We shall truly appreciate if the same be circulated among your friends, scholars and students encouraging them to participate in the meet.

http://210.212.237.38/iacrconference/

The Story of You: ENCODE and the human genome

Sat, 24 Aug 2013 18:49:03 -0500

Ever since a monk called Mendel started breeding pea plants we've been learning about our genomes. In 1953, Watson, Crick and Franklin described the structure of the molecule that makes up our genomes: the DNA double helix. Then, in 2001, scientists wrote down the entire 3-billion letter code contained in the average human genome. Now they're trying to interpret that code; to work out how it's used to make different types of cells and different people. The ENCODE project, as it's called, is the latest chapter in the story of you. To read the ENCODE research papers and more, visit http://www.nature.com/ENCODE