BOL: Related items

Edanchin Lab

Thu, 19 Nov 2020 08:00:07 -0600

My main topics of interest are:

The impact of non tree-like evolution such as horizontal gene transfers and hybridization on species biology
Evolution and adaptation of animals in the absence of sexual reproduction and the underlying mechanisms
Genomic signatures of adaptation to a parasitic life-style

More at https://edanchin.org/

Professor/Senior Lecturer of Comparative Genomics @ University of Glasgow

Fri, 06 Dec 2024 05:16:09 -0600

University of Glasgow
College of Medical, Veterinary and Life Sciences
School of Biodiversity, One Health and Veterinary Medicine

Professor/Senior Lecturer of Comparative Genomics
Vacancy Ref: 153610
Salary: Professor, Grade 10 will be within the Professorial range and
subject to negotiation
Senior Lecturer, Grade 9, 57,696 - 64,914 per annum

The School of Biodiversity, One Health and Veterinary Medicine has an
exciting opportunity to appoint a Professor/Senior Lecturer in Comparative
Genomics. You will make a substantial and positive contribution to the
strategic direction of the School/College through leading and contributing
to research of international standard, high quality teaching at both
undergraduate and postgraduate level, securing research funding, and
providing academic leadership and management within the School/College.

Applications are invited from candidates of international standing with
an appropriate record of academic achievement in comparative genomics
and associated omics technologies. We are looking for a candidate who
will complement our existing strengths in clinical veterinary medicine,
evolutionary biology, and animal physiology, with a demonstrable interest
in using domestic mammals among their study systems. We are particularly
interested in applications from candidates with a track record of
studying health related traits and their underlying genomic basis in
companion animals. Traits of specific interest include those related
to metabolism, ageing, and disease (e.g. cancer, autoimmune diseases,
neuromuscular disorders).

The School of Biodiversity, One Health and Veterinary Medicine is home to
researchers studying organismal biology and animal health across a diverse
range of systems, approaches and disciplines with existing strengths
in infectious disease, physiology, ageing, veterinary epidemiology, and
evolution among others. You will be based on the University of Glasgow's
Garscube campus, where the majority of veterinary teaching and research
infrastructure is located. This includes the Small Animal Hospital (a
recent 15M investment) and our Veterinary Diagnostic Services, offering
excellent opportunities for collaborative research at the clinical and
translational interface, especially with respect to companion animals.

We welcome applications from candidates with a Scottish Credit and
Qualification Framework level 12 (PhD) in animal biology, genomics and
health or related discipline with an extensive and established reputation
in research and significant teaching experience within the subject area.

This post is full time and open ended.

Visit our website for further information on The University of
Glasgow's, School of Biodiversity, One Health & Veterinary Medicine,
https://www.gla.ac.uk/schools/bohvm/

Informal Enquiries should be directed to Professor Roman Biek,
Roman.Biek@glasgow.ac.uk

Apply online at:
https://my.corehr.com/pls/uogrecruit/erq_jobspec_version_4.jobspec?p_id=153610

SATSUMA : Highly sensitive whole-genome synteny alignments.

Jit — Fri, 13 May 2016 05:25:26 -0500

Satsuma is a whole-genome synteny alignment program. It takes two genomes, computes alignments, and then keeps only the parts that are orthologous, i.e. following the conserved order and orientation of features, such as protein coding genes, non-coding genes, or neutral sequences. Satsuma does not require any pre-processing, such as repeat masking, since it will automatically detect ambiguous mappings.

Satsuma has parallelization built-in and is designed to run on multi-core architectures. The run-time for aligning two bird-size genomes (~1.2 Gb) is around two days on 24 CPUs.

You can find the manual here.
Download the latest source code from here.
Stable versions can also be downloaded from the Broad Institute's web site.

An incomplete list of questions and answers (yes, these have really been asked by our users! Please feel free to add your own by e-mailing us) is here.

If you use Satsuma in your research, please cite:
Grabherr, M. G., Russell, P., Meyer, M., Mauceli, E., Alföldi, J., Di Palma, F., & Lindblad-Toh, K. (2010). Genome-wide synteny through highly sensitive sequence alignment: Satsuma. Bioinformatics, 26(9), 1145-51.

Tutorial at http://evomics.org/learning/genomics/satsuma/

Address of the bookmark: http://satsuma.sourceforge.net/

PRISM

Jit — Sat, 10 Dec 2016 15:19:40 -0600

PRISM is a software for split read (reads which span across a structrual variant -- SV ) mapping and SV calling from the mapping result. PRISM is able to detect small insertions and abitrary size deletions, inversions and tandom duplications with the direction of discordant read pairs. PRISM_CTX is a tool for detecting inter-chromosome trans-location events.

PRISM and PRISM_CTX were originally designed and written by Michael Brudno and Yue Jiang, The original PRISM publication can be found here.

The authors may be contacted via e-mail at: prism at cs.toronto.edu.

Additional information is available in the PRISM README file and PRISM_CTX README file.

http://compbio.cs.toronto.edu/prism/

Address of the bookmark: http://compbio.cs.toronto.edu/prism/

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.

dcGOR

Martin Jones — Sat, 08 Nov 2014 14:54:28 -0600

An R package for analysing ontologies and protein domain annotations has been published in PLoS Computational Biology (http://dx.doi.org/10.1371/journal.pcbi.1003929). The package is distributed as part of CRAN (http://cran.r-project.org/package=dcGOR), and also at GitHub for version control.

The dedicated website is available in http://supfam.org/dcGOR, from which several demos are also provided:

1. Analysing SCOP domains: http://supfam.org/dcGOR/demo-Fang.html

2. Analysing Pfam domains: http://supfam.org/dcGOR/demo-Basu.html

3. Analysing InterPro domains: http://supfam.org/dcGOR/demo-Customisation.html

CANU: Assembling Large Genomes with Single-Molecule Sequencing and Locality Sensitive Hashing.

Jit — Tue, 26 Apr 2016 11:38:10 -0500

Canu is a fork of the Celera Assembler designed for high-noise single-molecule sequencing (such as the PacBio RSII or Oxford Nanopore MinION). The software is currently alpha level, feel free to use and report issues encountered.

Canu is a hierachical assembly pipeline which runs in four steps:

Detect overlaps in high-noise sequences using MHAP
Generate corrected sequence consensus
Trim corrected sequences
Assemble trimmed corrected sequences

Read the documentation

New release https://github.com/marbl/canu/releases

Address of the bookmark: https://github.com/marbl/canu

ORFfinder with smart BLAST

Jit — Tue, 17 May 2016 01:43:15 -0500

ORF Finder

ORFfinder is a graphical analysis tool for finding open reading frames (ORFs). We’ve been working on a few updates, and we’d like to find out what you think about them. Read on to find out what you can do with the new ORFfinder.

Smart BLAST (https://ncbiinsights.ncbi.nlm.nih.gov/2015/07/29/smartblast/)

Select one or a group of ORFs and BLAST several databases at once, and use the newly developed SmartBLAST to verify protein names. Looking for the traditional results from BLAST? They’re there too.

BBMap/BBTools package: Multipurpose tool designed for converting reads or other nucleotide data between different formats.

Jit — Mon, 13 Jun 2016 05:47:21 -0500

Reformatis a member of the BBMap/BBTools package. It is a multipurpose tool designed for converting reads or other nucleotide data between different formats. It supports, and can inter-convert:

fastq
fasta
fasta+qual
sam
scarf (an old Illumina format)
bam (if samtools is installed)
gzip
zip
ascii-33 (sanger)
ascii-64 (old Illumina)
paired files
interleaved files

It is multithreaded and can process data at over 500 megabytes per second, and can accept streams from standard in and write to standard out, allowing it to be easily dropped into the middle of a pipeline for format conversion. Reformat autodetects formats based on file extensions and content, making it very easy to use; and the autodetection can be overridden, allowing flexibility for people who don't like to follow naming conventions, or out-of-spec fastq files with qualities values like -17 or 120.

The program has been gradually expanded, and can now perform various other functions. None of these will break pairing, if the input is paired.

Quality trimming (either or both ends)
Quality filtering
Fixed-length trimming
Generation of histograms (base composition, quality, etc)
Subsampling (to a fraction of input reads, or an exact number of reads or bases)
Changing fasta line-wrapping length
Reverse-complementing (all reads or only read 2)
Adding /1 and /2 suffix to read names
GC-content filtering
Length-filtering
Testing for corrupted interleaved files

Reformat is compatible with any platform that supports Java 1.7 or higher. It also has a bash shellscript for simpler invocation. Typical usage examples:

Reformat fastq into fasta:
reformat.sh in=x.fq out=y.fa

Interleave paired reads:
reformat.sh in1=x1.fq in2=x2.fq out=y.fq

Note - you can actually use a shortcut if paired read files have the same name with a 1 and a 2. This is equivalent to the above command:
reformat.sh in=x#.fq out=y.fq

De-interleave reads:
reformat.sh in=x.fq out1=y1.fq out2=y2.fq

Verify that interleaving appears correct, assuming Illumina namimg conventions:
reformat.sh in=x.fq vint

Convert ASCII-33 to ASCII-64:
reformat.sh in=x.fq out=y.fq qin=33 qout=64

Quality-trim paired reads to Q10 on the left and right ends and discard reads shorter than 50bp after trimming:
reformat.sh in1=x1.fq in2=x2.fq out1=y1.fq out2=y2.fq outsingle=singletons.fq qtrim=rl trimq=10 minlength=50

Subsample 10% of the first 20000 pairs in an interleaved file:
reformat.sh in=x.fq out=y.fq reads=20000 samplerate=0.1 int=t
(in this case "int=t" overrides interleaving autodetection, to ensure reads are treated as pairs)

Pipe in a gzipped sam file and pipe out fasta:
reformat.sh in=stdin.sam.gz out=stdout.fa

Reverse-complement reads:
reformat.sh in=x.fq out=y.fq rcomp

For reformatting a file with very long sequences, Reformat will need more memory; just add the additional flag "-Xmx2g". For example, to change the line-wrapping length on the human genome (which has individual sequences over 200Mbp long) to 70 characters:
reformat.sh -Xmx2g in=HG19.fa.gz out=HG19_wrapped.fa.gz fastawrap=70

For additional functions, please run the shellscript with no arguments, or just read it with a text editor. If you have any questions, please post them in this thread.

For people using a non-bash terminal, you may need to type "bash reformat.sh" instead of just "reformat.sh".
For users of Windows or other platforms that do not support bash shellscripts, replace "reformat.sh" with "java -ea -Xmx200m /path/to/bbmap/current/ jgi.ReformatReads"
for example,
java -ea -Xmx200m C:\bbmap\current\ jgi.ReformatReads in=x.fq out=y.fa

Reformat can be downloaded with BBTools here:
https://sourceforge.net/projects/bbmap/

Mauve: a system for constructing multiple genome alignments in the presence of large-scale evolutionary events such as rearrangement and inversion

Jit — Sat, 24 Dec 2016 09:20:53 -0600

Mauve is a system for constructing multiple genome alignments in the presence of large-scale evolutionary events such as rearrangement and inversion. Multiple genome alignments provide a basis for research into comparative genomics and the study of genome-wide evolutionary dynamics.

Mauve has been developed with the idea that a multiple genome aligner should require only modest computational resources. It employs algorithmic techniques that scale well in the lengths of sequences being aligned. For example, a pair of Y. pestis genomes can be aligned in under a minute, while a group of 9 divergent Enterobacterial genomes can be aligned in a few hours. However, the current algorithm’s compute time (progressiveMauve) scales cubically in the number of genomes to align, making it unsuitable for datasets containing more than 50-100 bacterial genomes.

Address of the bookmark: http://darlinglab.org/mauve/mauve.html