• A simple intro to statistical distributions
• hypothesis testing
• linear models.

reading: http://compgenomr.github.io/book/stats.html

slides: https://github.com/BIMSBbioinfo/compgen2021/tree/main/week1/compgen2021_stats.pdf

exercises+code: https://github.com/BIMSBbioinfo/compgen2021/tree/main/week1/

Module 2: Unsupervised learning for genomics (9-15 August 2021)

• Understanding basic intuition behind machine learning approaches.
• Using unsupervised learning to cluster and visualise data points
• Dimension reduction techniques for visualisation and as input to clustering methods

reading: http://compgenomr.github.io/book/unsupervisedLearning.html

slides: https://github.com/BIMSBbioinfo/compgen2021/tree/main/week2/compgen2021_unsupervisedLearning.pdf

exercises+code: https://github.com/BIMSBbioinfo/compgen2021/tree/main/week2/

Module 3: Supervised learning for genomics (16-22 August 2021)

• Understanding and using supervised learning methods for predictive purposes
• How to measure prediction performance
• Understand and use cross-validation and related concepts

reading: http://compgenomr.github.io/book/supervisedLearning.html

slides: https://github.com/BIMSBbioinfo/compgen2021/tree/main/week3/compgen2021_supervisedLearning.pdf

exercises+code: https://github.com/BIMSBbioinfo/compgen2021/tree/main/week3/

https://github.com/BIMSBbioinfo/compgen2021

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

The vital metaphase stage of cell division, when the sister chromatids migrated to the centre and lined up in a row, and pulled apart using attached microtubules in such a way that half the DNA ends up in each daughter cell. However, before the mitotic spindle‐mediated movement gets start and pulled DNA apart, the chromosomes are free to undergo recombination which involves the exchange of genetic material either between multiple chromosomes or between different regions of the same chromosome.

During recombination, the precise breakage of each strand, exchange between the strands, and sealing of the resulting recombined molecules happens. The “chromosomal breakpoints” refers to these places where they break. Mostly, this process occurs with a high degree of accuracy at high frequency in both eukaryotic and prokaryotic cells. But occasionally this “break and sealing/ break and reattach” process goes wrong and the reattachment happens in the wrong place which usually create disaster (with few exceptions).These chromosome disaster or abnormalities involve the gain, loss or rearrangement of visible amounts of genetic material during cell division. These abnormalities are of two type, the first one is numerical abnormalities  where severe disorders are caused by the loss or gain of whole chromosomes, which affect the copy number of hundreds or even thousands of genes. The second are structural abnormalities which can be unbalanced or balanced. The former are similar to numerical abnormalities in that genetic material is either gained or lost. The natural defects in chromosome segregation are linked to cancer and several genetic diseases (http://en.wikipedia.org/wiki/List_of_genetic_disorders). Therefore, the enzymes involved in regulating cell division are still the attractive drug targets for many diseases.

Apart from certain chromosome abnormalities, these “crossing over” of segments of maternal and paternal chromosomes to form hybrid chromosomes have some evolutionary importance and considered as a driver of genetic variation. Moreover, the chromosome breakage in evolution is considered to be non-random in nature(http://www.ploscompbiol.org/article/info%3Adoi%2F10.1371%2Fjournal.pcbi.0020014). In addition the study of breakpoint regions and non-breakpoint (stable) regions of chromosomes indicates both the regions evolved in distinctly different ways ( http://www.ncbi.nlm.nih.gov/pmc/articles/PMC2675965/). These breakage may lead to genetic diseases or participate to chromosomal rearranmgnets and contributed in development of new species.

I will try to explain the genome hotspots/Evolutionary Breakpoint Regions(EBRs)/fragile regions/weak fragments/  in my next blog.

Software for recombination detection:

RAT http://cbr.jic.ac.uk/dicks/software/RAT/

Breakpointer https://github.com/ruping/Breakpointer

DRP http://web.cbio.uct.ac.za/~darren/rdp.html

RB-finder http://www.ncbi.nlm.nih.gov/pubmed/18707535

LDhat2.0 http://ldhat.sourceforge.net/LDhat2.0/instructions.shtml

Reference:

http://www.nature.com/scitable/topicpage/genetic-recombination-514#

Image: Wikipedia , sciencelearn.org.nz

Recommended Articles:

http://www.friendshipcircle.org/blog/2012/05/22/13-chromosomal-disorders-youve-never-heard-of/

http://web.udl.es/usuaris/e4650869/docencia/segoncicle/genclin98/recursos_classe_%28pdf%29/revisionsPDF/chromosyndromes.pdf

http://www.ncbi.nlm.nih.gov/pmc/articles/PMC2775595/table/T2/

http://learn.genetics.utah.edu/content/disorders/chromosomal/

http://www.ncert.nic.in/html/learning_basket/biology/cc&cd.pdf

More at https://www.igib.res.in/bdmg/ScientistRecruitmentAdvt2023.pdf

Why Normalize RNA-Seq Data?

Normalization is a crucial step in RNA-Seq analysis for the following reasons:

• Sequencing depth: Different RNA-Seq experiments produce varying numbers of reads, making direct comparisons between samples misleading.

• Gene length: Longer genes inherently generate more reads, irrespective of their actual expression level.

• Bias reduction: Normalization mitigates technical biases, enabling meaningful biological interpretation.

TPM (Transcripts Per Million)

TPM measures the proportion of reads mapped to a transcript, normalized by transcript length and sequencing depth. It is calculated as:

Key Features:

1. Proportionality: TPM values sum to 1,000,000 across all transcripts in a sample, making it easier to compare between samples.

2. Intuitive interpretation: TPM values directly represent the abundance of transcripts in a sample.

3. Preferred for comparisons: TPM facilitates between-sample comparisons better than FPKM.

FPKM (Fragments Per Kilobase Million)

FPKM normalizes read counts by transcript length and sequencing depth, but without enforcing proportionality like TPM. It is defined as:

Key Features:

1. Historical significance: FPKM was one of the first normalization methods used for RNA-Seq.

2. Single-end vs. paired-end: In paired-end sequencing, FPKM becomes RPKM (Reads Per Kilobase Million).

3. Limited utility: FPKM values are not as robust as TPM for cross-sample comparisons due to lack of proportionality.

CPM (Counts Per Million)

CPM normalizes raw read counts by sequencing depth, without considering gene length. It is expressed as:

Key Features:

1. Simplicity: CPM is straightforward and computationally less intensive.

2. Application: Suitable for non-length-dependent analyses, such as comparing total expression levels or differential expression analysis.

3. Gene length agnostic: CPM does not correct for gene length, making it less ideal for measuring expression levels.

When to Use Each Method

• TPM: Best for comparing expression levels between samples, especially when transcript length and sequencing depth vary.

• FPKM: Useful for historical consistency but generally replaced by TPM.

• CPM: Ideal for differential expression analysis when gene length normalization is unnecessary.

Conclusion

Choosing the right normalization method depends on the specific objectives of your RNA-Seq analysis. TPM’s proportionality and robustness make it the preferred choice for most applications, while CPM serves well for differential expression studies. Although FPKM paved the way for RNA-Seq normalization, it has largely been supplanted by TPM in modern workflows. Understanding these methods and their nuances ensures accurate and meaningful interpretations of RNA-Seq data.

References:

1. Li, B., & Dewey, C. N. (2011). RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics.

2. Trapnell, C., et al. (2010). Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nature Biotechnology.

3. Law, C. W., et al. (2014). voom: precision weights unlock linear model analysis tools for RNA-seq read counts. Genome Biology.

Webinar on 'An integrated RNA and DNA approach to unravel genetic regulation in cancer'

Abstract

Whole exome DNA sequencing (WES) or whole genome DNA sequencing (WGS) allows detection of mutations and polymorphisms in all exonic and genomic regions, respectively, while messenger RNA sequencing (RNA-Seq) enables quantitative analysis of gene expression. Mutations in the genome result in diverse transcriptional aberrations that can be missed in a stand-alone WES/WGS analysis. An integration of DNA variant analysis and RNA-Seq analysis enables one to investigate the consequences of genomic changes in the RNA transcripts including germline and somatic changes, imprinting, RNA editing and allele specific expression (ASE). In this webinar, we will demonstrate this integrated approach using Strand NGS to identify high confidence mutations, RNA editing events and ASE in cancer.

Webinar Details

Register here: http://www.strand-ngs.com/webinar_registration

About Speaker:

Dr. Veena Hedatale, has a PhD in Plant Genetics from The Radboud University, Netherlands focused on meiosis and recombination. Her prior academic experience at Cornell University was on genetic mapping and gene transformation in Rice. She has worked with Monsanto, and contributed to data mining, database development as well as gene/promoter/pathway discovery for traits related to yield and stress in crop species. At Strand, Veena has worked on Pharmacogenomic analysis of targets and Gene family analysis projects. Currently, she is part of the Strand NGS Application Science team and is involved in the analysis of next generation sequencing data.

Please feel free to contact us 24/5, for availing free online training or if you have any questions.

Address of the bookmark: https://github.com/nedaz/kraken

Address of the bookmark: http://chmille4.github.io/Scribl/

http://www.phyloxml.org/archaeopteryx-js/adh.html

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

https://github.com/tanghaibao/jcvi

More at https://github.com/tanghaibao/jcvi/wiki

Address of the bookmark: https://github.com/tanghaibao/jcvi

Address of the bookmark: https://github.com/thackl/gggenomes