BOL: Related items

Luo Lab | Symbiosis Genomics & Evolution

Wed, 09 Sep 2026 03:30:30 -0500

We study the evolutionary genomics of marine invertebrates to understand their origins and diversity. Our lab combines high-throughput sequencing and single-cell transcriptomics to explore a wide range of non-model systems. We are particularly interested in how evolutionary novelty arises, with a focus on animal development and photosymbiosis.

Research Directions

Stony corals: evolution of novelty and photosymbiosis

Symbiotic acoels: cell type evolution and photosymbiosis

Animal genomes: structural evolution and gene regulation

https://sgel.biodiv.tw/home

The Variant Everyone Ignored

BioStar — Mon, 05 Oct 2026 12:14:20 -0500

Consider a scenario in which a patient's genome has been sequenced. Among billions of DNA bases, a structural alteration may explain the patient's disease. Multiple advanced algorithms analyze the data, yet only one detects the variant, while the others do not. In standard bioinformatics workflows, such a solitary result is often regarded as unreliable and subsequently discarded. Although the solution exists within the data, prevailing computational protocols may overlook it.

A recent study published in Genome Biology (https://link.springer.com/article/10.1186/s13059-026-04280-y) addressed this challenge by introducing dicast (https://github.com/burgshrimps/dicast), a machine-learning approach for detecting structural variants in short-read sequencing data. Structural variants, such as large deletions, insertions, duplications, and inversions, can have significant biological and clinical implications, yet they are challenging to identify with short-read technologies. Because different detection methods frequently yield divergent results, researchers commonly employ consensus calling, considering a variant valid only if multiple tools detect it. While this approach reduces false positives, it relies on the potentially flawed assumption that the majority is always correct.

The researchers explored the impact of evaluating the supporting evidence for each variant, rather than simply tallying the number of algorithms that identified it. To establish a ground truth, they analyzed nine genomes using multiple sequencing technologies and 15 detection methods, initially identifying approximately 35 million potential variants. Through extensive filtering, evidence integration, and manual review of over 11,500 variants, they developed a robust benchmark comprising more than 236,000 structural variants. The findings underscored the complexity of the problem: short-read methods detected fewer than half of deletions and less than 10 percent of insertions, with performance declining markedly in repetitive genomic regions. In contrast, long-read technologies demonstrated superior detection capabilities. However, replacing the substantial volume of existing short-read data in clinical and research settings is not immediately feasible. Consequently, the researchers questioned whether short-read data might harbor more information than conventional analytical pipelines currently extract.

This line of inquiry led to the development of dicast. Rather than merely confirming agreement among multiple tools, dicast identifies patterns in sequencing data, including split and clipped reads, discordant read pairs, alignment characteristics, and the surrounding genomic context. An XGBoost machine-learning model evaluates which combinations of these signals are indicative of genuine structural variants. Thus, the approach shifts from tallying algorithmic consensus to interpreting the underlying evidence.

The researchers subsequently conducted a targeted evaluation by examining structural variants detected by only a single short-read tool, which are typically missed by consensus-based approaches. dicast successfully recovered approximately 81% of these single-caller deletions, insertions, and duplications. The signals for these variants were present in the data, but conventional filtering methods failed to integrate them effectively.

The utility of dicast was further demonstrated in cohorts with rare diseases, including congenital limb malformations, atrial fibrillation, and neuromuscular disorders. In one instance, dicast achieved a deletion recall rate of approximately 0.96, compared to 0.74 using consensus calling. The median number of variants requiring manual review was 29 per sample. Among 31 experimentally validated variants that standard filters would have missed, dicast identified 12, whereas consensus calling detected only one.

Overall, dicast identified approximately 20 percent more potential disease-causing deletions than consensus-based methods. While a 20 percent increase may appear modest, in clinical genomics such improvements can have significant practical implications. Missing a deletion may leave a case unresolved, whereas detecting a structural variant can provide critical diagnostic insights.

The study does not claim that machine learning has rendered short-read sequencing superior to long-read approaches. Instead, the results underscore the effectiveness of long-read sequencing for structural variant detection. However, dicast highlights a more nuanced perspective: substantial biological information may still be recoverable from the extensive short-read datasets already available.

The principal lesson extends beyond the detection of structural variants. For many years, bioinformatics pipelines have relied on threshold-based criteria, such as minimum coverage, quality scores, or support from multiple tools. While these rules are useful, biological phenomena do not always conform to rigid checklists; multiple weak signals, when considered collectively, can provide compelling evidence.

This perspective prompts consideration of the solitary variant: one algorithm identifies it, while several others do not. Traditional consensus techniques might have dismissed it, yet machine learning approaches evaluate the available evidence to determine whether the variant is plausible.

Occasionally, the most significant variant within a genome is the one that is almost universally overlooked.

Read more at https://link.springer.com/article/10.1186/s13059-026-04280-y

Genetic Test in India

Rahul Agarwal — Sun, 11 Aug 2013 10:54:35 -0500

1.Xcode Life Sciences Pvt. Ltd.
6B, Eldorado,
112, Nungambakkam High Road,
Nungambakkam, Chennai 600034
Tamil Nadu, India

2.Mapmygenome™
Royal Demeure,HUDA Techno Enclave,
Plot No. 12/2, Sector-1 500 081
Madhapur,Hyderabad
AP, India

3. DNA Labs India

http://www.dnalabsindia.com/lab.php

4.MedGenome Labs Pvt Ltd
(Division of SciGenom Labs Pvt Ltd.)
Plot no: 43A,SDF, 3rd floor
A Block,CSEZ, Kakanad, Cochin
Kerala - 682037
Phone: 0484 - 2413399
Fax: 0484 - 2413398
Email: info@medgenome.com

5.Narayana Nethralaya

Narayana Hrudayalaya Campus
Narayana Health City
# 258/A, Bommasandra, Hosur Road,
Bangalore - 560 099 - INDIA.
TEL: +91-80-66660655-0658
FAX: +91-80-66660650
Mobile: 9902 821128 (Emergency Only)
e-mail: info@narayananethralaya.com

6.BioAxis DNA Research Centre Private Limited
13-51,Sri Lakshmi Nagar colony,
Besides Big Bazar, Near Kamineni Hospitals
GSI Post BandalGuda (L B Nagar) Hydeabad-500068
Andhra Pradesh (India).
Phone : +91-40-24034503/+91-9246338983

7.Gene Guiide

8th Floor, Embassy Towers, 7 Bungalows Rd, Versova, Andheri West, Mumbai-61
09167 117799
info@geneguiide.com

See more at: http://www.geneguiide.com

8.INDIAN BIOSCIENCES
Regd. Office:
G-2 (Ground Floor Rear), Kailash Colony, New Delhi - 110048, India.
Phone: +91 (0)11 29236088, Email: info@inbdna.com.

9.SRL Limited

GP-26, MARUTI INDUSTRIAL ESTATE,

UDYOG VIHAR,SECTOR-18,

GURGAON - 122015

Tel: 0124-3001243 / 0124-3001209

SRL Limited
VASANT VIHAR, 8, PALAM MARG,
NEW DELHI - 110057
Tel: 011 - 4229 5333

Website: http://www.srlworld.com
National Customer care number:
Call Toll Free : 1800-222-660/1800-102-8282
E-mail id: customercare@srl.in

10.Tata Memorial Centre,

Advanced Centre for Treatment, Research and Education in Cancer

Kharghar, Navi Mumbai - 410 210, INDIA.

Tel: +91-22-2740 5000

Fax: +91-22-2740 5085

E-mail: mail@actrec.gov.in

NCBI Webinar

Jit — Sun, 08 Jun 2014 02:47:01 -0500

In less than two weeks, NCBI will offer a webinar entitled "Introducing 3 NCBI Resources to Navigate Testing for Disease Linked Variants: MedGen, GTR and ClinVar". This webinar will delve into the lifecycle of genetic testing and teach attendees how to navigate the NIH Genetic Testing Registry, ClinVar, and MedGen resources. These resources can be used to prepare for clinical cases, access detailed information about orderable genetic tests, interpret test results, and more.

More at https://attendee.gotowebinar.com/register/8452228815737989634

Sandelin group

Mon, 30 Sep 2013 19:12:58 -0500

Sandelin group have a deep interest in most biology, but focus on gene regulation and the many areas that are connected with this, including transcriptomics, epigenetics and technological and informatics aspects.

The group is both computational and experimental.

We ask biological questions to large datasets made using novel genomics techniques, with the help of computers. One of the strengths in the group are the many connections to high-profile experimental laboratories which supply data to be analyzed.

Lab webpage @ http://people.binf.ku.dk/albin/Sandelin_group_at_the_Bioinformatic_Centre/The_Sandelin_group.html

Molecular Genetics Lecture

Rahul Agarwal — Fri, 27 Sep 2013 04:24:45 -0500

"Robert Sapolsky makes interdisciplinary connections between behavioral biology and molecular genetic influences. He relates protein synthesis and point mutations to microevolutionary change, and discusses conflicting theories of gradualism and punctuated equilibrium and the influence of epigenetics on development theories."

"Robert Sapolsky is an American neuroendocrinologist, professor of biology, neuroscience, and neurosurgery at Stanford University, researcher and author" ----Wikipedia

Address of the bookmark: http://www.youtube.com/watch?v=_dRXA1_e30o

Prof. Chandrasekhar Kanduri Laboratory

Fri, 24 Jun 2016 16:59:43 -0500

Our lab has been interested in understanding how long noncoding RNAs control tumor initiation and progression, in addition to use them as potential biomarkers in diagnosis and therapy. We have been using neuroblastoma, a childhood cancer, as a model system to understand the functional role of long noncoding RNAs in cancer development and progression. By using new RNA sequencing technology on neuroblastoma tumors from a large group of Swedish children including both high-risk and low-risk neuroblastomas (108), we have identified several long noncoding RNAs that could have potential role in diagnosis and therapy. We are currently exploring the functional role of these differentially expressed long noncoding RNA in nuroblastoma progression and development.

Prof. Dr. rer. nat. Jeanette Erdmann's Lab

Sun, 17 Dec 2017 10:29:18 -0600

Her prestigious lab working on Integrative and Experimental genomics to understand the underlying mechanism thats able to define an inherited predisposition of coronary artery disease (CAD) and myocardial infarction (MI) via GWAS, different Sequencing technologies and Functional genomics. Her lab is located at Institute for Cardio-genetics Universität zu Lübeck, Maria-Goeppert-Str. 1, 23562 Lübeck.

Tel +49 451 3101 8301
cardiogenetics(at)uni-luebeck.de

http://cardiogenetics-luebeck.de/

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/

Ribbon: Visualizing complex genome alignments and structural variation:

Jit — Wed, 29 Nov 2017 07:40:22 -0600

Ribbon can be used for long reads, short reads, paired-end reads, and assembly/genome alignments. Instructions for each data format are available by clicking on "instructions" in each tab on the right.

Local installation:

You can install Ribbon locally from Github by following the instructions here: https://github.com/MariaNattestad/Ribbon

Address of the bookmark: http://genomeribbon.com/