Genome Annotation and RNA Sequencing Guide
Genome Annotation and RNA Sequencing Guide
Genome Annotation
What is Genome Annotation? Genome annotation is the process of identifying and describing
the functional elements within a genome. This includes identifying genes, regulatory regions,
and other functional sequences, and then attaching biological information to these elements.
The goal is to understand what a genome "does" and how its various parts contribute to the
organism's biology.
Levels of Genome Annotation:
Genome annotation generally occurs at two main levels:
1. Structural Annotation: This level focuses on identifying the physical locations of
genomic features. It involves:
○ Gene prediction: Locating protein-coding genes, tRNA genes, rRNA genes, and
other non-coding RNA genes.
○ Identification of regulatory elements: Finding promoters, enhancers, terminators,
and other DNA sequences that control gene expression.
○ Repeat identification: Locating repetitive sequences, which can include
transposable elements.
○ Identification of open reading frames (ORFs): Determining potential
protein-coding regions.
2. Functional Annotation: Once structural elements are identified, functional annotation
assigns biological meaning to them. This involves:
○ Assigning gene names and symbols: Based on sequence similarity to known
genes or functional characterization.
○ Predicting protein function: Inferring the biochemical activity, biological process,
and cellular component of proteins encoded by predicted genes. This often involves
comparing sequences to databases of known proteins (e.g., UniProt, Pfam, GO).
○ Identifying protein domains and motifs: Regions within proteins that have
specific structures or functions.
○ Mapping to metabolic pathways: Integrating genes into known biochemical
pathways (e.g., KEGG, MetaCyc).
○ Identifying post-translational modifications: Predicting sites where proteins
might be modified after translation.
Which are the two phases of Genome Annotation? The two main phases of genome
annotation are:
1. Structural Annotation: As described above, this phase focuses on identifying the
physical locations of genes and other functional elements within the genomic sequence.
2. Functional Annotation: This phase involves assigning biological functions and roles to
the identified structural elements.
RNA Sequencing
Discuss the steps for RNA Sequencing
RNA sequencing (RNA-Seq) is a powerful technology that uses next-generation sequencing to
reveal the presence and quantity of RNA in a biological sample at a given moment. The general
steps are:
1. Experimental Design:
○ Define research question: What biological question are you trying to answer?
○ Sample selection: Choose appropriate biological samples (e.g., treated vs. control,
different tissues, time points).
○ Replication: Include biological replicates to ensure statistical power and account
for biological variability.
○ Sequencing depth: Determine the desired number of reads per sample based on
the experimental goals and complexity of the transcriptome.
2. RNA Extraction and Quality Control:
○ Isolate high-quality total RNA from the chosen samples.
○ Assess RNA quantity (spectrophotometer) and integrity (Bioanalyzer/gel
electrophoresis). High-quality RNA is crucial for successful RNA-Seq.
3. Library Preparation (as discussed above):
○ Enrich for mRNA (poly(A) selection or rRNA depletion).
○ Fragment RNA.
○ Synthesize first-strand cDNA.
○ Synthesize second-strand cDNA.
○ End repair and A-tailing.
○ Ligate sequencing adapters.
○ Amplify the library via PCR.
○ Perform library quality control.
4. Sequencing:
○ Load the prepared cDNA library onto a sequencing platform (e.g., Illumina, Oxford
Nanopore, PacBio).
○ The sequencer performs massively parallel sequencing, generating millions to
billions of short sequence reads (or longer reads depending on the technology).
5. Data Analysis (Bioinformatics): This is a multi-step computational process:
○ Demultiplexing: If multiple samples were sequenced together (multiplexed), the
reads are separated based on their unique barcode sequences.
○ Quality Control (QC): Assess the quality of raw sequencing reads (e.g., using
FastQC). Remove low-quality reads, adapter sequences, and trim poor-quality
bases.
○ Read Alignment/Mapping: Align the high-quality reads to a reference genome or
transcriptome. This step identifies the genomic location from which each read
originated. Tools like STAR, HISAT2, Bowtie2 are commonly used.
○ Transcript Assembly (for de novo RNA-Seq): If a reference genome is not
available, reads are assembled de novo into transcripts.
○ Quantification: Count the number of reads that map to each gene or transcript.
This provides an estimate of gene expression levels (e.g., using HTSeq, Salmon,
Kallisto, featureCounts).
○ Normalization: Adjust raw read counts to account for differences in library size,
gene length, and other technical variations. This allows for meaningful comparisons
between samples.
○ Differential Expression Analysis: Identify genes that are significantly up- or
down-regulated between different experimental conditions (e.g., using DESeq2,
edgeR, Limma-Voom).
○ Downstream Analysis:
■ Functional enrichment analysis: Identify over-represented biological
pathways, GO terms, or functions among differentially expressed genes.
■ Clustering: Group genes with similar expression patterns or samples with
similar transcriptional profiles.
■ Pathway analysis: Map differentially expressed genes to known biological
pathways.
■ Visualization: Create plots and figures to represent gene expression data
(e.g., heatmaps, volcano plots, MA plots).
Gene Prediction
Why is gene prediction important? Gene prediction is paramount in genomics for several
reasons:
1. Unlocking Biological Function: Identifying genes is the first step towards understanding
the functional content of a genome. Genes encode proteins and functional RNAs, which
are the workhorses of the cell. Without knowing where genes are, it's impossible to
understand an organism's biology.
2. Foundation for Functional Annotation: Accurate gene prediction is the prerequisite for
functional annotation. Once a gene is predicted, its sequence can be compared to
databases of known proteins to infer its function, pathways, and potential roles in disease.
3. Disease Research: Many diseases are caused by mutations in genes. Identifying genes
allows researchers to pinpoint potential disease-causing mutations and develop
diagnostics and therapies.
4. Drug Discovery: Genes encoding proteins that are potential drug targets need to be
identified and characterized. Gene prediction helps in this process.
5. Evolutionary Studies: By comparing gene content and organization across different
species, researchers can infer evolutionary relationships and understand how genomes
evolve.
6. Agricultural Improvement: Identifying genes related to traits like yield, disease
resistance, and nutritional value in crops and livestock can facilitate breeding programs.
7. Synthetic Biology: Accurate gene prediction is essential for designing and synthesizing
new biological systems.
Define gene prediction? Gene prediction is the computational process of identifying
protein-coding genes, as well as non-coding RNA genes, within a genomic DNA sequence. It
involves analyzing the sequence for specific patterns and signals that indicate the presence of a
gene, such as start codons, stop codons, splice sites (in eukaryotes), ribosomal binding sites (in
prokaryotes), and open reading frames (ORFs).
What is comparative gene prediction? Comparative gene prediction is a method that
leverages sequence conservation across multiple related species to improve the accuracy of
gene prediction. It is based on the principle that functional genomic elements, including genes,
are often more conserved during evolution than non-functional sequences.
How it works:
● Genomes from several related species are aligned.
● Regions of high sequence conservation are identified. These conserved regions are more
likely to contain functional elements, including genes.
● Information from well-annotated genomes (e.g., known genes in a closely related species)
is used to guide gene prediction in less-annotated genomes.
● Algorithms search for gene-like signals (ORFs, splice sites) within these conserved
regions, using the evolutionary constraint as additional evidence.
Benefits of comparative gene prediction include higher accuracy, especially for identifying genes
with weaker intrinsic signals, and the ability to identify orthologous genes across species.
RAST Tool
Define Fig Fams in RAST tool? In the context of the RAST (Rapid Annotations Subsystems
Technology) tool, "FigFams" (Figurative Families) refer to a curated collection of protein
families. Each FigFam represents a group of homologous proteins that are believed to perform
the same function, based on sequence similarity and genomic context.
RAST uses these FigFams as a key component of its functional annotation process. When
annotating a bacterial or archaeal genome, RAST compares the predicted proteins to the
FigFams database. If a protein matches a FigFam, it's assigned the associated functional
annotation, which includes information about its role in specific metabolic pathways or cellular
processes. This approach helps to quickly and accurately assign functions to genes, particularly
within the context of Subsystems (see RAST definition below).
Define RAST tool. RAST (Rapid Annotations using Subsystems Technology) is a popular
and widely used bioinformatics tool and web server designed for rapid and accurate annotation
of bacterial and archaeal genomes. It automates much of the genome annotation pipeline,
making it accessible to researchers without extensive bioinformatics expertise.
Key features and functionalities of RAST include:
● Automated Gene Prediction: It predicts protein-coding genes, rRNA genes, and tRNA
genes.
● Functional Annotation: It assigns functions to predicted genes by comparing them to a
comprehensive database of curated protein families (FigFams) and integrating them into
Subsystems.
● Subsystems-based Annotation: A core concept in RAST is the "Subsystem." A
subsystem is a collection of functionally related genes that collectively carry out a specific
metabolic pathway or cellular process. RAST organizes annotations into these
subsystems, providing a more holistic view of an organism's metabolism.
● User-friendly Interface: It provides a web-based interface for uploading genomes,
running annotations, and visualizing the results.
● Speed: As its name suggests, RAST is designed for rapid annotation, making it suitable
for high-throughput genomic projects.
● Metabolic Reconstruction: By organizing genes into subsystems, RAST helps in
reconstructing metabolic pathways and understanding the physiological capabilities of the
organism.
RAST has been instrumental in accelerating the annotation of thousands of microbial genomes,
contributing significantly to our understanding of microbial diversity and function.
Phred
Define Phred Phred is a computer program used to analyze and assign quality scores to the
bases in DNA sequencing reads, particularly from Sanger sequencing, but its underlying
principles have influenced quality scoring in Next-Generation Sequencing (NGS) as well.
Key aspects of Phred:
● Base Calling: Phred processes the raw chromatogram data (fluorescence traces)
generated by DNA sequencers to determine the most likely base at each position (A, C,
G, or T).
● Quality Scores: For each base call, Phred assigns a Phred quality score (Q-score),
which is a logarithmic measure of the probability that a base call is incorrect. The higher
the Q-score, the lower the probability of error, and thus the higher the confidence in the
base call.
The Phred quality score is defined as: Q = -10 \log_{10} P where P is the estimated probability
of the base call being incorrect.
● Interpretation of Q-scores:
○ Q10: 1 in 10 chance of error (90% accuracy)
○ Q20: 1 in 100 chance of error (99% accuracy)
○ Q30: 1 in 1000 chance of error (99.9% accuracy)
○ Q40: 1 in 10000 chance of error (99.99% accuracy)
● Importance: Phred quality scores are critical for:
○ Filtering low-quality data: Bioinformatic pipelines use these scores to trim or
discard unreliable base calls and reads.
○ Variant calling: High-quality scores are essential for accurate identification of
SNPs and other genetic variants.
○ Assembly: Helps in resolving ambiguities during genome assembly by prioritizing
higher quality bases.
○ Comparison of sequencing data: Provides a standardized metric for assessing
the quality of sequencing runs.
Although initially developed for Sanger sequencing, the concept of Phred quality scores has
been widely adopted and adapted in NGS platforms (e.g., Illumina often reports Phred-scaled
quality scores).
Microarrays
Why are microarrays important? Microarrays, despite the rise of next-generation sequencing,
have been and remain important for several reasons:
1. High-Throughput Gene Expression Profiling: Microarrays allow for the simultaneous
measurement of the expression levels of thousands to tens of thousands of genes in a
single experiment. This high-throughput capability made them revolutionary for
understanding global gene expression changes.
2. Cost-Effectiveness (for targeted studies): For studies focusing on a pre-defined set of
genes (e.g., known pathways, specific disease markers), microarrays can still be more
cost-effective than RNA sequencing, especially when large numbers of samples are
involved.
3. Historical Data and Legacy Studies: A vast amount of publicly available gene
expression data was generated using microarrays. This historical data is invaluable for
meta-analyses, validating new findings, and exploring hypotheses without generating new
data.
4. Clinical Diagnostics and Prognostics: Microarrays have been used to identify gene
expression signatures for various diseases, aiding in diagnosis, prognosis, and predicting
response to therapy (e.g., in some cancer diagnostics).
5. Biomarker Discovery: They were instrumental in the early days of identifying potential
biomarkers for various biological processes and diseases.
6. Simplicity and Established Workflows: The technology and analysis pipelines for
microarrays are well-established and relatively straightforward compared to some NGS
applications.
7. Copy Number Variation (CNV) Analysis: Beyond gene expression, comparative
genomic hybridization (CGH) arrays (a type of microarray) are used for detecting copy
number variations across the genome.
While RNA-Seq offers advantages like novel transcript discovery, higher dynamic range, and
detection of splice variants, microarrays still hold their place for specific applications, particularly
when focusing on established gene sets or leveraging existing datasets.
Define Microarrays. A microarray (also known as a gene chip or biochip) is a laboratory tool
used to detect the expression of thousands of genes simultaneously or to detect variations in
DNA. It consists of a solid surface (typically a glass slide) to which thousands of microscopic
spots of DNA probes are attached in a predefined grid pattern. Each spot contains a specific
DNA sequence (probe) that is complementary to a known gene or DNA segment.
The basic principle of a gene expression microarray experiment involves:
1. Probe immobilization: Short, single-stranded DNA sequences (probes) representing
specific genes are synthesized and spotted onto the array.
2. Sample preparation: mRNA is extracted from experimental samples (e.g., diseased vs.
healthy tissue), converted into cDNA, and labeled with fluorescent dyes.
3. Hybridization: The labeled cDNA from the samples is denatured (made single-stranded)
and then incubated with the microarray. The labeled cDNA molecules will bind (hybridize)
only to their complementary probes on the array.
4. Detection: After washing away unbound molecules, the array is scanned with a laser to
detect the fluorescence intensity at each spot. The intensity of the fluorescence signal at a
particular spot is proportional to the abundance of the corresponding mRNA in the original
sample.
5. Data analysis: The raw intensity data is then processed and analyzed to determine
differential gene expression between samples.
There are different types of microarrays, including gene expression microarrays, comparative
genomic hybridization (CGH) arrays, and SNP arrays, each designed for specific applications.
SAM
What is SAM? In the context of next-generation sequencing and bioinformatics, SAM stands
for Sequence Alignment Map. It is a tab-delimited text file format used to store information
about sequence alignments to a reference genome.
Key characteristics of SAM:
● Alignment Information: Each line in a SAM file represents a single read and its
alignment to the reference. It includes details such as:
○ Read ID
○ Flags (e.g., whether the read is paired, mapped, reverse complemented)
○ Reference sequence name (chromosome/contig)
○ Start position of the alignment
○ Mapping quality (Phred-scaled score indicating confidence of alignment)
○ CIGAR string (Compact Idiosyncratic Gapped Alignment Report) which describes
how the read aligns to the reference (e.g., matches, mismatches, insertions,
deletions)
○ Mate information (for paired-end reads)
○ Sequence of the read
○ Quality scores of the read
○ Optional tags for additional information
● Human-Readable: SAM files are plain text, making them relatively easy to read and
parse, though they can be very large.
● BAM (Binary Alignment Map): Because SAM files can be enormous, they are typically
converted into their compressed binary equivalent, BAM. BAM files are significantly
smaller, indexed, and much faster for computational processing and visualization (e.g., in
genome browsers). Most bioinformatics tools prefer to work with BAM files.
SAM/BAM is a foundational format in NGS data analysis pipelines, serving as the standard
output for alignment tools (e.g., Bowtie2, BWA, STAR) and the input for downstream analyses
like variant calling, peak calling (for ChIP-seq), and quantification.
Gene Expression
What is gene expression? Gene expression is the process by which information from a gene
is used in the synthesis of a functional gene product, such as a protein or a functional RNA
(e.g., tRNA, rRNA, or non-coding RNA). It is the fundamental process that allows cells to
produce the molecules they need to carry out their functions.
The central dogma of molecular biology describes the primary flow of genetic information in
most organisms: DNA → RNA → Protein
Key stages of gene expression in eukaryotes include:
1. Transcription: The DNA sequence of a gene is copied into an RNA molecule
(messenger RNA or mRNA) by RNA polymerase.
2. RNA Processing (in eukaryotes): The newly synthesized primary RNA transcript
(pre-mRNA) undergoes modifications, including splicing (removal of introns and joining of
exons), capping at the 5' end, and polyadenylation at the 3' end. These modifications are
necessary for the mRNA to be stable, transported out of the nucleus, and translated
efficiently.
3. Translation: The processed mRNA molecule is transported to the cytoplasm, where
ribosomes "read" the mRNA sequence and synthesize a polypeptide chain (protein)
based on the genetic code.
4. Post-translational Modifications (for proteins): The newly synthesized protein may
undergo further modifications (e.g., folding, cleavage, addition of chemical groups) to
become fully functional.
Regulation of Gene Expression: Gene expression is tightly regulated at multiple levels
(transcription, RNA processing, translation, post-translational modification) to ensure that the
right genes are expressed at the right time, in the right cell type, and in the right amount,
responding to environmental cues and developmental programs. Measuring gene expression
(e.g., through RNA-Seq or microarrays) provides insights into the molecular state and function
of a cell or tissue.
SAGE
What is SAGE? SAGE stands for Serial Analysis of Gene Expression. It is a molecular
biology technique developed in the mid-1990s (by Velculescu et al.) that was used to
quantitatively measure the expression levels of thousands of genes simultaneously in a given
cell or tissue sample. It was a pioneering technique for transcriptome profiling before the advent
of microarrays and later, RNA sequencing.
The core principle of SAGE involves:
1. Generating Short Sequence Tags: SAGE relies on the idea that a short sequence tag
(typically 9-14 base pairs) at a defined position within an mRNA molecule is sufficient to
uniquely identify that mRNA.
2. Concatenating Tags: These tags are then enzymatically extracted, ligated together to
form long concatemer molecules, and cloned into vectors.
3. Sequencing Concatemers: The concatemers are then sequenced using traditional
Sanger sequencing methods.
4. Counting Tags: The individual tags are computationally extracted from the sequenced
concatemers, and the frequency of each unique tag is counted.
5. Quantification: The abundance of a specific tag directly correlates with the expression
level of the corresponding mRNA in the original sample.
Advantages of SAGE (historically):
● Quantitative: Provided an absolute measure of gene expression levels (tag counts).
● Discovery-oriented: Could identify novel transcripts without prior knowledge of their
sequences.
● High-throughput for its time: Enabled the analysis of thousands of genes.
Limitations of SAGE:
● Technical complexity: Involved many enzymatic steps, making it labor-intensive and
prone to errors.
● Requires relatively large amounts of RNA.
● Short tag length: Sometimes, the 9-14 bp tag was not long enough to unambiguously
identify a gene, especially in larger genomes with many paralogs.
● Reliance on Sanger sequencing: Limited throughput and high cost compared to NGS.
While SAGE has largely been superseded by RNA sequencing due to the latter's superior
throughput, dynamic range, and ability to capture full transcript information (including splice
variants), it was a significant advancement in its time and laid conceptual groundwork for later
high-throughput gene expression profiling technologies.
Illumina Sequencing
Discuss Illumina sequencing Illumina sequencing (also known as Sequencing by Synthesis,
SBS) is the most widely adopted and dominant next-generation sequencing (NGS) technology
used today. It is renowned for its high throughput, accuracy, and relatively low cost per base,
making it suitable for a vast array of genomic applications.
Here's a discussion of the key steps and principles of Illumina sequencing:
1. Library Preparation:
○ DNA or RNA (converted to cDNA) is fragmented into appropriately sized pieces
(typically 200-500 bp).
○ Specific adapters are ligated to both ends of each fragment. These adapters
contain sequences for binding to the flow cell, priming sequencing reactions, and
unique barcodes (indexes) for multiplexing.
2. Cluster Generation (Bridge Amplification):
○ The adapter-ligated library is loaded onto a flow cell. A flow cell is a glass slide
with millions of oligonucleotides (primers) complementary to the library adapters,
covalently attached to its surface.
○ Each single-stranded library fragment binds to a complementary oligo on the flow
cell.
○ Bridge amplification then occurs:
■ The bound fragment bends over to form a "bridge" and anneals to an
adjacent oligo.
■ DNA polymerase extends from both oligos, creating a double-stranded
bridge.
■ The double-stranded bridges are denatured, leaving two single-stranded
copies of the original fragment, each bound to a different oligo.
■ This process is repeated many times, creating millions of identical copies (a
cluster) of the original single fragment at a discrete spot on the flow cell.
Each cluster contains thousands of copies of a single unique fragment.
3. Sequencing by Synthesis (SBS):
○ After cluster generation, the flow cell is placed into the Illumina sequencer.
○ Sequencing primers bind to the adapters on the clusters.
○ Reversible terminator dNTPs (deoxyribonucleotide triphosphates) are added.
These dNTPs are fluorescently labeled (each base with a different color) and have
a reversible terminator group that prevents further incorporation.
○ DNA polymerase incorporates a single labeled dNTP complementary to the
template strand.
○ The flow cell is then imaged by a laser, capturing the fluorescence signal from each
cluster. The color detected indicates the incorporated base.
○ After imaging, a chemical cleavage step removes the fluorescent dye and the
reversible terminator, allowing the next base to be incorporated.
○ This cycle (incorporation, imaging, cleavage) is repeated tens to hundreds of times,
building up the sequence of each cluster base by base.
4. Data Analysis:
○ The images are processed to identify the sequence of each cluster, generating raw
reads (FASTQ files).
○ These reads are then demultiplexed (if indexed samples were run together),
quality-controlled, aligned to a reference genome, and subjected to various
downstream bioinformatics analyses depending on the application (e.g., variant
calling, gene expression quantification, ChIP-seq peak calling).
Key Advantages of Illumina Sequencing:
● High Throughput: Generates billions of reads and terabytes of data per run.
● High Accuracy: Low error rates, especially for single nucleotide polymorphisms.
● Cost-Effective: Low cost per base compared to older sequencing methods.
● Scalability: Various instrument models (MiSeq, NextSeq, NovaSeq) cater to different
throughput needs.
● Versatility: Applicable to a wide range of applications (whole genome sequencing, exome
sequencing, RNA-Seq, ChIP-Seq, metagenomics, etc.).
Limitations:
● Short Read Lengths: Traditionally, Illumina produces relatively short reads (50-300 bp),
which can be challenging for resolving highly repetitive regions or complex structural
variants. However, newer chemistries are pushing towards longer reads.
● GC Bias: Can exhibit some bias in regions of very high or low GC content.
● Bridge Amplification Issues: Can lead to some amplification bias and potential
homopolymer issues.
Next-Generation Sequencing
Define next generation sequencing Next-generation sequencing (NGS), also known as
high-throughput sequencing, refers to a suite of technologies that have revolutionized genomics
by enabling the sequencing of millions to billions of DNA or RNA molecules simultaneously and
in parallel. This massively parallel approach dramatically increases sequencing throughput and
reduces the cost per base compared to traditional Sanger sequencing.
Key characteristics that define NGS include:
● Massive Parallelism: Instead of sequencing one DNA fragment at a time (like Sanger),
NGS platforms sequence countless fragments in parallel, generating vast amounts of data
in a single run.
● High Throughput: Produces orders of magnitude more sequence data than conventional
methods.
● Reduced Cost: Significantly lowered the cost of sequencing, making large-scale genomic
studies feasible.
● Digital Data: Provides quantitative data in the form of read counts, enabling precise
measurement of molecular abundance.
● Library Preparation: All NGS technologies require the preparation of a "library" of DNA
or cDNA fragments that are specifically modified to be compatible with the sequencing
platform.
● Diverse Applications: NGS has enabled a wide range of applications, including:
○ Whole-genome sequencing (WGS)
○ Whole-exome sequencing (WES)
○ RNA sequencing (RNA-Seq) for gene expression profiling
○ ChIP-sequencing (ChIP-Seq) for protein-DNA interactions
○ Metagenomics for microbial community analysis
○ Variant calling (SNPs, indels, structural variants)
○ Epigenomic studies (e.g., methylation sequencing)
How it differs from Sanger Sequencing: Sanger sequencing, the "first generation" method,
relies on chain termination and typically sequences one fragment at a time, producing long but
relatively few reads. NGS, on the other hand, fragments the genome, sequences millions of
these fragments in parallel, and then computationally reassembles the genome or maps the
reads to a reference.
Examples of NGS technologies include Illumina Sequencing, Ion Torrent, Oxford Nanopore
Technologies, and PacBio SMRT Sequencing.
Spliceman Tool
Define Spliceman tool As of my last update, a prominent and widely recognized bioinformatics
tool specifically named "Spliceman" is not a standard or commonly cited tool in the broad field of
genomics or RNA sequencing. It's possible that:
● It's a very specialized tool used within a particular research group or niche.
● It's a newer tool that hasn't gained widespread recognition yet.
● There might be a misunderstanding or a typo in the name.
● It could be a component or a specific module within a larger bioinformatics suite.
If you have more context or can provide details about its function, I might be able to offer a more
specific answer.
However, given the name "Spliceman," it strongly suggests a tool involved in splicing analysis,
particularly related to:
● Alternative Splicing Detection: Identifying different mRNA isoforms produced from a
single gene.
● Splice Junction Identification: Locating the precise boundaries between exons and
introns.
● Splicing Aberration Detection: Finding abnormal splicing events that might be linked to
disease.
● Splicing Factor Binding Site Prediction: Identifying regulatory elements that control
splicing.
Tools that are commonly used for splicing analysis in RNA-Seq data include:
● TopHat/TopHat2 (older, but foundational for alignment across junctions)
● STAR (Spliced Transcripts Alignment to a Reference - highly popular for RNA-Seq
alignment)
● HISAT2 (another popular splice-aware aligner)
● DEXSeq, rMATS, MISO, LeafCutter, SUPPA2 (for differential splicing analysis)
● JunctionSeq
● ASprofile
If "Spliceman" is indeed a specific tool, it would likely fit into one of these categories for
analyzing splicing events from RNA sequencing data.
SNP
Discuss SNP. SNP stands for Single Nucleotide Polymorphism. It is the most common type
of genetic variation among people. A SNP is a variation at a single nucleotide position in a DNA
sequence among individuals or groups.
Key characteristics of SNPs:
● Single Base Change: It involves a change in a single base pair (e.g., at a specific
position in the genome, an 'A' might be present in one individual, while a 'G' is present in
another).
● Commonality: For a DNA variation to be classified as an SNP, it must occur in at least
1% of the population. If it occurs in less than 1%, it is usually considered a rare variant.
● Biallelic: Most SNPs are biallelic, meaning there are only two possible alleles (two
different nucleotides) at that specific genomic location.
● Abundance: SNPs are extremely abundant in the human genome, occurring
approximately every 100 to 300 base pairs. There are millions of SNPs in the human
genome.
● Location: SNPs can occur in various regions of the genome:
○ Coding regions (exons): Can lead to changes in amino acid sequence (missense,
nonsense) or silent changes.
○ Non-coding regions (introns, intergenic regions): May affect gene regulation
(e.g., alter promoter or enhancer activity), mRNA splicing, or have no known
functional consequence.
Causes of SNPs: SNPs primarily arise from spontaneous mutations during DNA replication
or repair processes. These mutations are then passed down through generations.
Environmental factors and mutagens can also contribute to the rate of mutations.
Example of an SNP: Consider a small segment of DNA sequence from two individuals:
Individual 1: A G C T **A** G T C Individual 2: A G C T **G** G T C
At the fifth position (highlighted), Individual 1 has an 'A', while Individual 2 has a 'G'. If this
difference is common (occurs in >1% of the population), it is an SNP. This particular SNP might
be referred to as rsID_12345 (an identifier from the dbSNP database), and its alleles would be A
and G.
Importance and Applications of SNPs:
● Genetic Markers: SNPs serve as powerful genetic markers for:
○ Disease Association Studies (GWAS): Identifying SNPs that are statistically
associated with complex diseases (e.g., diabetes, heart disease, cancer) or traits.
○ Pharmacogenomics: Predicting an individual's response to drugs based on their
genetic makeup.
○ Ancestry and Population Genetics: Tracing human migration patterns and
understanding population relationships.
○ Forensics: DNA fingerprinting and identification.
○ Agricultural Breeding: Identifying traits in crops and livestock.
● Understanding Disease Mechanisms: SNPs in coding or regulatory regions can directly
influence gene function, protein structure, or gene expression, thereby contributing to
disease susceptibility or resistance.
● Personalized Medicine: Informing tailored treatments and preventative strategies based
on an individual's unique genetic profile.
RFLP
What is RFLP? RFLP stands for Restriction Fragment Length Polymorphism. It is a
molecular biology technique that exploits variations in homologous DNA sequences
(polymorphisms) that result in different patterns of DNA fragment lengths after digestion with
restriction enzymes.
How it works (RFLP Technique):
1. DNA Isolation: Genomic DNA is extracted from a sample (e.g., blood, tissue).
2. Restriction Enzyme Digestion: The isolated DNA is cut into fragments using one or
more restriction enzymes. Restriction enzymes are bacterial enzymes that recognize
and cleave DNA at specific, short nucleotide sequences (restriction sites).
3. Agarose Gel Electrophoresis: The resulting DNA fragments are separated by size using
agarose gel electrophoresis. Smaller fragments migrate faster and further down the gel
than larger fragments.
4. Southern Blotting (Optional but common): The DNA fragments are then transferred
from the agarose gel to a solid membrane (e.g., nylon) through a process called Southern
blotting. This makes the DNA accessible for hybridization.
5. Hybridization with Probes: The membrane is incubated with a labeled DNA probe. This
probe is a short, single-stranded DNA sequence that is complementary to a specific
region of interest in the fragmented DNA. The probe will only bind (hybridize) to its
complementary fragments on the membrane.
6. Visualization: The hybridized probe is then detected (e.g., via autoradiography if
radiolabeled, or chemiluminescence if fluorescently labeled), revealing a unique pattern of
bands corresponding to the sizes of the restriction fragments that hybridized to the probe.
Why is it used?
RFLP is used because genetic variations (polymorphisms) can:
● Create a new restriction site: A base change might create a sequence recognized by a
restriction enzyme, leading to shorter fragments at that locus.
● Destroy an existing restriction site: A base change might eliminate a restriction site,
leading to longer fragments at that locus.
These changes in restriction sites result in different-sized DNA fragments (polymorphisms in
fragment lengths) among individuals, which can be visualized.
Benefits of RFLP:
● Genetic Mapping: Historically, RFLP markers were crucial for constructing the first
genetic linkage maps of human and other genomes. By tracking the inheritance of RFLP
patterns, researchers could determine the relative positions of genes and markers on
chromosomes.
● Disease Gene Mapping: It was used to identify the chromosomal location of genes
responsible for genetic diseases by observing co-segregation of an RFLP marker with the
disease trait in families.
● Paternity Testing and Forensics: Due to the unique RFLP patterns among individuals, it
was used for DNA fingerprinting in paternity testing and criminal investigations.
● Genetic Diversity Studies: Assessing genetic variation within and between populations.
● Diagnosis of Genetic Diseases: In some cases, RFLP could directly detect mutations
that alter restriction sites and are associated with a disease (e.g., Sickle Cell Anemia,
where a mutation abolishes an MstII restriction site).
Challenges of RFLP:
● Labor-Intensive and Time-Consuming: The multi-step process (digestion,
electrophoresis, Southern blotting, hybridization, detection) is laborious and takes several
days to complete.
● Requires Large Amounts of High-Quality DNA: Insufficient or degraded DNA can lead
to poor results.
● Radioactive Probes (Historically): Often relied on radioactively labeled probes,
requiring specialized safety protocols.
● Low Throughput: Only one or a few markers could be analyzed at a time, making it
unsuitable for large-scale genome-wide studies.
● Reliance on Specific Probes: Requires prior knowledge and availability of suitable
probes for the region of interest.
● Less Sensitive: Compared to PCR-based methods, RFLP is less sensitive and cannot
amplify small amounts of DNA.
Due to these challenges, RFLP has largely been superseded by more efficient, high-throughput,
and sensitive PCR-based methods (like PCR-RFLP, microsatellites, and SNP genotyping
arrays, and eventually NGS) for most applications. However, it remains a foundational technique
in molecular biology that paved the way for subsequent genomic advancements.
Immunoprecipitation
What is immunoprecipitation? How does it work?
What is Immunoprecipitation? Immunoprecipitation (IP) is a widely used laboratory
technique in molecular biology and biochemistry to isolate and purify a specific protein (or a
complex of proteins) from a complex mixture, such as a cell lysate, using an antibody that
specifically binds to that protein.
How does it work?
The principle of immunoprecipitation relies on the highly specific binding affinity between an
antibody and its target antigen (the protein of interest). Here's a step-by-step breakdown:
1. Preparation of Cell/Tissue Lysate:
○ Cells or tissues are lysed (broken open) using a suitable buffer to release their
intracellular contents, including the target protein. The lysis buffer usually contains
detergents to solubilize membranes, and protease inhibitors to prevent protein
degradation.
2. Antibody Binding (Incubation):
○ A specific antibody that recognizes and binds to the protein of interest is added to
the cell lysate.
○ The mixture is incubated (typically on a rocker or rotator for several hours to
overnight, at 4°C) to allow the antibody to bind to its target protein. This forms an
antibody-antigen complex.
3. Capture of Antibody-Antigen Complex (using beads):
○ To separate the antibody-antigen complex from the rest of the cell lysate, affinity
beads are added. These beads are typically made of agarose, Sepharose, or
magnetic particles and are coated with a protein (like Protein A or Protein G) that
has a high affinity for the Fc region of antibodies.
○ The beads bind to the antibody (which is already bound to the target protein),
effectively "pulling down" the entire antibody-antigen complex from the solution.
○ The mixture is again incubated to allow sufficient binding.
4. Washing:
○ The beads with the bound complexes are then extensively washed multiple times
with a mild buffer. This step is crucial to remove non-specifically bound proteins and
other cellular debris, ensuring the purity of the isolated protein.
○ For non-magnetic beads, this involves repeated centrifugation and resuspension.
For magnetic beads, a magnet is used to pellet the beads, and the supernatant is
removed.
5. Elution:
○ After washing, the target protein (and any co-purified proteins) needs to be
separated from the beads and the antibody. This is achieved by adding an elution
buffer that disrupts the antibody-antigen interaction or the antibody-bead interaction.
Common elution methods include:
■ Low pH buffer (e.g., glycine buffer): Disrupts antibody-antigen binding.
■ High pH buffer: Similar to low pH.
■ SDS-PAGE loading buffer with heating: Denatures the proteins and
dissociates them from the beads and antibody, typically for subsequent
SDS-PAGE analysis.
■ Competitive elution: Adding a high concentration of the specific antigen or a
peptide that competes for antibody binding.
6. Downstream Analysis:
○ The eluted protein sample is then analyzed using various techniques to confirm the
presence of the target protein, identify interacting partners, or quantify its amount.
Common downstream analyses include:
■ Western Blotting: To confirm the presence and size of the target protein.
■ Mass Spectrometry (MS): To identify the target protein and any
co-immunoprecipitated (interacting) proteins.
■ Enzymatic Assays: To test the activity of the isolated protein.
■ DNA/RNA analysis: In techniques like ChIP (Chromatin
Immunoprecipitation) or RIP (RNA Immunoprecipitation) where DNA or RNA
associated with a protein is isolated.
Variants of Immunoprecipitation:
● Co-Immunoprecipitation (Co-IP): Used to identify protein-protein interactions. If a
protein interacts with the target protein, it will also be "pulled down" with the target by the
antibody.
● Chromatin Immunoprecipitation (ChIP): Used to study protein-DNA interactions (e.g.,
transcription factors binding to DNA, histone modifications).
● RNA Immunoprecipitation (RIP): Used to study protein-RNA interactions.
Immunoprecipitation is a fundamental technique for studying protein function, protein
interactions, and post-translational modifications in a cellular context.