Hubbry Logo
Nucleic acid sequenceNucleic acid sequenceMain
Open search
Nucleic acid sequence
Community hub
Nucleic acid sequence
logo
8 pages, 0 posts
0 subscribers
Be the first to start a discussion here.
Be the first to start a discussion here.
Nucleic acid sequence
Nucleic acid sequence
Nucleic acid sequence
View on Wikipedia
from Wikipedia
Nucleic acid primary structureNucleic acid secondary structureNucleic acid tertiary structureNucleic acid quaternary structure
The image above contains clickable links
The image above contains clickable links
Interactive image of nucleic acid structure (primary, secondary, tertiary, and quaternary) using DNA helices and examples from the VS ribozyme and telomerase and nucleosome. (PDB: ADNA, 1BNA, 4OCB, 4R4V, 1YMO, 1EQZ​)

A nucleic acid sequence is a succession of bases within the nucleotides forming alleles within a DNA (using GACT) or RNA (GACU) molecule. This succession is denoted by a series of a set of five different letters that indicate the order of the nucleotides. By convention, sequences are usually presented from the 5' end to the 3' end. For DNA, with its double helix, there are two possible directions for the notated sequence; of these two, the sense strand is used. Because nucleic acids are normally linear (unbranched) polymers, specifying the sequence is equivalent to defining the covalent structure of the entire molecule. For this reason, the nucleic acid sequence is also termed the primary structure.

The sequence represents genetic information. Biological deoxyribonucleic acid represents the information which directs the functions of an organism.

Nucleic acids also have a secondary structure and tertiary structure. Primary structure is sometimes mistakenly referred to as "primary sequence". However there is no parallel concept of secondary or tertiary sequence.

Nucleotides

[edit]
Chemical structure of RNA
A series of codons in part of a mRNA molecule. Each codon consists of three nucleotides, usually representing a single amino acid.
Main article: Nucleotide

Nucleic acids consist of a chain of linked units called nucleotides. Each nucleotide consists of three subunits: a phosphate group and a sugar (ribose in the case of RNA, deoxyribose in DNA) make up the backbone of the nucleic acid strand, and attached to the sugar is one of a set of nucleobases. The nucleobases are important in base pairing of strands to form higher-level secondary and tertiary structures such as the famed double helix.

The possible letters are A, C, G, and T, representing the four nucleotide bases of a DNA strand – adenine, cytosine, guanine, thymine – covalently linked to a phosphodiester backbone. In the typical case, the sequences are printed abutting one another without gaps, as in the sequence AAAGTCTGAC, read left to right in the 5' to 3' direction. With regards to transcription, a sequence is on the coding strand if it has the same order as the transcribed RNA.

One sequence can be complementary to another sequence, meaning that they have the base on each position in the complementary (i.e., A to T, C to G) and in the reverse order. For example, the complementary sequence to TTAC is GTAA. If one strand of the double-stranded DNA is considered the sense strand, then the other strand, considered the antisense strand, will have the complementary sequence to the sense strand.

Notation

[edit]
Main article: Nucleic acid notation

While A, T, C, and G represent a particular nucleotide at a position, there are also letters that represent ambiguity which are used when more than one kind of nucleotide could occur at that position. The rules of the International Union of Pure and Applied Chemistry (IUPAC) are as follows:[1]

For example, W means that either an adenine or a thymine could occur in that position without impairing the sequence's functionality.

List of symbols
Symbol[2] Meaning/derivation Possible bases Complement
A Adenine A 1 T (or U)
C Cytosine C G
G Guanine G C
T Thymine T A
U Uracil U A
W Weak A T 2 S
S Strong C G W
M aMino A C K
K Keto G T M
R puRine A G Y
Y pYrimidine C T R
B not A (B comes after A) C G T 3 V
D not C (D comes after C) A G T H
H not G (H comes after G) A C T D
V not T (V comes after T and U) A C G B
N any Nucleotide (not a gap) A C G T 4 N
Z Zero 0 Z

These symbols are also valid for RNA, except with U (uracil) replacing T (thymine).[1]

Apart from adenine (A), cytosine (C), guanine (G), thymine (T) and uracil (U), DNA and RNA also contain bases that have been modified after the nucleic acid chain has been formed. In DNA, the most common modified base is 5-methylcytidine (m5C). In RNA, there are many modified bases, including pseudouridine (Ψ), dihydrouridine (D), inosine (I), ribothymidine (rT) and 7-methylguanosine (m7G).[3][4] Hypoxanthine and xanthine are two of the many bases created through mutagen presence, both of them through deamination (replacement of the amine-group with a carbonyl-group). Hypoxanthine is produced from adenine, and xanthine is produced from guanine.[5] Similarly, deamination of cytosine results in uracil.

Example of comparing and determining the % difference between two nucleotide sequences
  • AATCCGCTAG
  • AAACCCTTAG

Given the two 10-nucleotide sequences, line them up and compare the differences between them. Calculate the percent difference by taking the number of differences between the DNA bases divided by the total number of nucleotides. In this case there are three differences in the 10 nucleotide sequence. Thus there is a 30% difference.

Biological significance

[edit]
A depiction of the genetic code, by which the information contained in nucleic acids are translated into amino acid sequences in proteins.
Further information: Genetic code and Central dogma of molecular biology

In biological systems, nucleic acids contain information which is used by a living cell to construct specific proteins. The sequence of nucleobases on a nucleic acid strand is translated by cell machinery into a sequence of amino acids making up a protein strand. Each group of three bases, called a codon, corresponds to a single amino acid, and there is a specific genetic code by which each possible combination of three bases corresponds to a specific amino acid.

The central dogma of molecular biology outlines the mechanism by which proteins are constructed using information contained in nucleic acids. DNA is transcribed into mRNA molecules, which travel to the ribosome where the mRNA is used as a template for the construction of the protein strand. Since nucleic acids can bind to molecules with complementary sequences, there is a distinction between "sense" sequences which code for proteins, and the complementary "antisense" sequence, which is by itself nonfunctional, but can bind to the sense strand.

Sequence determination

[edit]
Electropherogram printout from automated sequencer for determining part of a DNA sequence
Main article: DNA sequencing

DNA sequencing is the process of determining the nucleotide sequence of a given DNA fragment. The sequence of the DNA of a living thing encodes the necessary information for that living thing to survive and reproduce. Therefore, determining the sequence is useful in fundamental research into why and how organisms live, as well as in applied subjects. Because of the importance of DNA to living things, knowledge of a DNA sequence may be useful in practically any biological research. For example, in medicine it can be used to identify, diagnose and potentially develop treatments for genetic diseases. Similarly, research into pathogens may lead to treatments for contagious diseases. Biotechnology is a burgeoning discipline, with the potential for many useful products and services.

RNA is not sequenced directly. Instead, it is copied to a DNA by reverse transcriptase, and this DNA is then sequenced.

Current sequencing methods rely on the discriminatory ability of DNA polymerases, and therefore can only distinguish four bases. An inosine (created from adenosine during RNA editing) is read as a G, and 5-methyl-cytosine (created from cytosine by DNA methylation) is read as a C. With current technology, it is difficult to sequence small amounts of DNA, as the signal is too weak to measure. This is overcome by polymerase chain reaction (PCR) amplification.

Digital representation

[edit]
Genetic sequence in digital format.

Once a nucleic acid sequence has been obtained from an organism, it is stored in silico in digital format. Digital genetic sequences may be stored in sequence databases, be analyzed (see Sequence analysis below), be digitally altered and be used as templates for creating new actual DNA using artificial gene synthesis.

Sequence analysis

[edit]
Main article: Sequence analysis

Digital genetic sequences may be analyzed using the tools of bioinformatics to attempt to determine its function.

Genetic testing

[edit]
Main article: Genetic testing

The DNA in an organism's genome can be analyzed to diagnose vulnerabilities to inherited diseases, and can also be used to determine a child's paternity (genetic father) or a person's ancestry. Normally, every person carries two variations of every gene, one inherited from their mother, the other inherited from their father. The human genome is believed to contain around 20,000–25,000 genes. In addition to studying chromosomes to the level of individual genes, genetic testing in a broader sense includes biochemical tests for the possible presence of genetic diseases, or mutant forms of genes associated with increased risk of developing genetic disorders.

Genetic testing identifies changes in chromosomes, genes, or proteins.[6] Usually, testing is used to find changes that are associated with inherited disorders. The results of a genetic test can confirm or rule out a suspected genetic condition or help determine a person's chance of developing or passing on a genetic disorder. Several hundred genetic tests are currently in use, and more are being developed.[7][8]

Sequence alignment

[edit]
Main article: Sequence alignment

In bioinformatics, a sequence alignment is a way of arranging the sequences of DNA, RNA, or protein to identify regions of similarity that may be due to functional, structural, or evolutionary relationships between the sequences.[9] If two sequences in an alignment share a common ancestor, mismatches can be interpreted as point mutations and gaps as insertion or deletion mutations (indels) introduced in one or both lineages in the time since they diverged from one another. In sequence alignments of proteins, the degree of similarity between amino acids occupying a particular position in the sequence can be interpreted as a rough measure of how conserved a particular region or sequence motif is among lineages. The absence of substitutions, or the presence of only very conservative substitutions (that is, the substitution of amino acids whose side chains have similar biochemical properties) in a particular region of the sequence, suggest[10] that this region has structural or functional importance. Although DNA and RNA nucleotide bases are more similar to each other than are amino acids, the conservation of base pairs can indicate a similar functional or structural role.[11]

Computational phylogenetics makes extensive use of sequence alignments in the construction and interpretation of phylogenetic trees, which are used to classify the evolutionary relationships between homologous genes represented in the genomes of divergent species. The degree to which sequences in a query set differ is qualitatively related to the sequences' evolutionary distance from one another. Roughly speaking, high sequence identity suggests that the sequences in question have a comparatively young most recent common ancestor, while low identity suggests that the divergence is more ancient. This approximation, which reflects the "molecular clock" hypothesis that a roughly constant rate of evolutionary change can be used to extrapolate the elapsed time since two genes first diverged (that is, the coalescence time), assumes that the effects of mutation and selection are constant across sequence lineages. Therefore, it does not account for possible differences among organisms or species in the rates of DNA repair or the possible functional conservation of specific regions in a sequence. (In the case of nucleotide sequences, the molecular clock hypothesis in its most basic form also discounts the difference in acceptance rates between silent mutations that do not alter the meaning of a given codon and other mutations that result in a different amino acid being incorporated into the protein.) More statistically accurate methods allow the evolutionary rate on each branch of the phylogenetic tree to vary, thus producing better estimates of coalescence times for genes.

Sequence motifs

[edit]
Main article: Sequence motif

Frequently the primary structure encodes motifs that are of functional importance. Some examples of sequence motifs are: the C/D[12] and H/ACA boxes[13] of snoRNAs, Sm binding site found in spliceosomal RNAs such as U1, U2, U4, U5, U6, U12 and U3, the Shine-Dalgarno sequence,[14] the Kozak consensus sequence[15] and the RNA polymerase III terminator.[16]

Sequence entropy

[edit]

In bioinformatics, a sequence entropy, also known as sequence complexity or information profile,[17] is a numerical sequence providing a quantitative measure of the local complexity of a DNA sequence, independently of the direction of processing. The manipulations of the information profiles enable the analysis of the sequences using alignment-free techniques, such as for example in motif and rearrangements detection.[17][18][19]

See also

[edit]

References

[edit]
[edit]
Revisions and contributorsEdit on WikipediaRead on Wikipedia

Nucleic acid sequence

View on Grokipedia
from Grokipedia
A nucleic acid sequence is a polymer composed of nucleotides that forms the primary structure of deoxyribonucleic acid (DNA) or ribonucleic acid (RNA), serving as the fundamental carrier of genetic information in living organisms. Each nucleotide consists of a nitrogenous base, a five-carbon sugar (deoxyribose in DNA or ribose in RNA), and a phosphate group, with the bases—adenine (A), guanine (G), cytosine (C), thymine (T) in DNA, or uracil (U) replacing thymine in RNA—linked in a specific linear order that encodes instructions for biological processes. By convention, these sequences are written and read from the 5' end to the 3' end, reflecting the directionality of the phosphodiester bonds that connect the nucleotides. In DNA, the sequence typically forms a double-stranded helix where complementary bases pair (A with T, G with C), enabling stable storage of genetic data across generations, while RNA sequences are generally single-stranded and versatile, functioning in roles such as (mRNA) for protein synthesis or regulatory non-coding RNAs. The , for instance, comprises approximately 3 billion base pairs of DNA sequence, underscoring the immense informational capacity of these molecules. Variations in sequences, known as , can lead to , , and diseases, making central to fields like and . Discovered in the late by , were later recognized for their sequence-based coding of through landmark experiments in the mid-20th century.

Components and Representation

Nucleotides

Nucleotides are the fundamental monomeric units that compose nucleic acid sequences, each consisting of a nitrogenous base, a five-carbon pentose sugar, and one or more phosphate groups attached to the sugar. The pentose sugar is ribose in ribonucleic acid (RNA) or 2'-deoxyribose in deoxyribonucleic acid (DNA), differing by the absence of a hydroxyl group at the 2' carbon position in deoxyribose. The phosphate group is typically linked to the 5' carbon of the sugar, forming a nucleotide monophosphate, though di- or triphosphate forms occur in metabolic contexts. The nitrogenous bases in nucleotides are heterocyclic aromatic compounds classified into two main types: purines and pyrimidines. Purines, and , feature a double-ring structure—a six-membered ring fused to a five-membered ring—with nitrogen atoms at positions 1, 3, 7, and 9. has an amino group at position 6, while has a at position 6 and an amino group at position 2. Pyrimidines, , , and uracil (U), possess a single six-membered ring with nitrogens at positions 1 and 3; has an amino group at position 4 and a at position 2, has a at position 5 and at positions 2 and 4, and uracil mirrors without the . In DNA, the canonical bases are , , , and , while in RNA, uracil substitutes for .
BaseTypeDNARNAKey Structural Features
Adenine (A)PurineYesYesFused pyrimidine-imidazole rings; amino group at C6
Guanine (G)PurineYesYesFused rings; carbonyl at C6, amino at C2
Cytosine (C)PyrimidineYesYesSingle ring; amino at C4, carbonyl at C2
Thymine (T)PyrimidineYesNoSingle ring; methyl at C5, carbonyls at C2 and C4
Uracil (U)PyrimidineNoYesSingle ring; carbonyls at C2 and C4
Nucleotides polymerize through phosphodiester bonds, where the 5' of one links to the 3' hydroxyl group of another via a , forming the sugar- backbone that provides structural integrity to the chain. This backbone alternates sugar and units, with the nitrogenous bases projecting inward or outward depending on the 's conformation. Beyond the canonical bases, nucleic acids contain non-canonical or modified bases, such as in (), which is an of with the base attached via a C-C rather than the standard N-glycosidic linkage. Other examples include dihydrouridine and , which arise from post-transcriptional modifications to standard bases.

Notation Systems

Nucleic acid sequences are symbolically represented using single-letter abbreviations for the four standard nucleotide bases. In deoxyribonucleic acid (DNA), these are A for adenine, C for cytosine, G for guanine, and T for thymine. For ribonucleic acid (RNA), uracil (U) replaces thymine, resulting in A, C, G, and U. These abbreviations, established as a compact notation for sequence description, facilitate clear communication in scientific literature and databases. By convention, nucleic acid sequences are written in the 5' to 3' direction, reflecting the polarity of the sugar-phosphate backbone where the 5' end terminates in a group attached to the 5' carbon of the , and the 3' end has a free hydroxyl group on the 3' carbon. This directionality aligns with the biochemical processes of replication and transcription, which proceed from 5' to 3'. For example, a short DNA sequence might be denoted as 5'-ATCG-3', indicating the order of bases from the 5' end to the 3' end. sequences follow the same convention, such as 5'-AUCG-3'. The prefixes 5' and 3' are often omitted when the direction is unambiguous, but explicit notation is used for clarity, especially in diagrams or when specifying strands. To handle uncertainty or variability in sequencing data, the International Union of Pure and Applied Chemistry (IUPAC) and the International Union of Biochemistry (IUB), now IUBMB, introduced ambiguity codes in their recommendations. These single-letter symbols represent groups of bases, allowing concise notation for degenerate or polymorphic sites. For instance, N denotes any base (A, C, G, or T/U), R specifies a (A or G), and Y indicates a (C or T/U). The full set of IUPAC codes is as follows:
SymbolBases RepresentedComplementary BasesOrigin of Designation
AAT
CCG
GGC
T (DNA)/U (RNA)T/UA/Uracil
RA or GY
YC or T/UR
MA or CKaMino
KG or T/UMKeto
SC or GSStrong (3 H-bonds)
WA or T/UWWeak (2 H-bonds)
HA or C or T/UDnot-G (H)
BC or G or T/UVnot-A (B)
VA or C or GBnot-T/U (V)
DA or G or T/UHnot-C (D)
NA or C or G or T/UNaNy
These codes ensure that complementary relationships are preserved, with each symbol mapping to its complement (e.g., complements ). The standardized notation evolved from early manual representations in the mid-20th century, which used full names or chemical formulas, to a unified system formalized by the IUPAC-IUB Commission on Biochemical Nomenclature in 1970. This addressed the growing need for consistency as techniques advanced, providing rules for abbreviations, sequence direction, and in publications. The 1970 recommendations were later refined in 1984 to incorporate expanded ambiguity symbols for incompletely specified sequences, reflecting improvements in sequencing accuracy. For double-stranded sequences, notation distinguishes the two antiparallel strands, which are connected via Watson-Crick base pairing: adenine pairs with thymine (A-T) in DNA or uracil (A-U) in RNA, and guanine pairs with cytosine (G-C). The sense strand (often the coding strand) is typically written 5' to 3' on the top line, with its complement below in the 3' to 5' direction to reflect the antiparallel orientation. For example: 5'-ATCG-3'
3'-TAGC-5' Ambiguity codes can be applied to either strand, with complementary symbols used for the opposite strand (e.g., an R on one strand corresponds to a Y on the complement). This format highlights base pairing and is essential for representing genomic regions or restriction sites.

Biological Roles

Genetic Information in DNA

DNA serves as the primary genetic material in most organisms, carrying the instructions necessary for development, functioning, growth, and . These sequences are organized into , which are compact structures consisting of DNA wrapped around proteins, enabling efficient storage and transmission during . In eukaryotic cells, the nuclear genome is divided among multiple linear s, while prokaryotes typically maintain a single circular . The central dogma of molecular biology posits that genetic information is stored in DNA and flows unidirectionally to RNA and then to proteins, with DNA acting as the stable repository for hereditary information. First articulated by Francis Crick in 1958, this framework emphasizes DNA's role in information storage and replication, ensuring the faithful transmission of genetic instructions across generations. Genes within the DNA sequence represent functional units that encode proteins or regulatory RNAs, structured with coding regions known as exons interspersed with non-coding introns, as well as upstream promoters that initiate transcription. For instance, the human β-globin gene consists of three exons separated by two introns, where exons contain the coding sequence (e.g., the first exon includes codons for the N-terminal amino acids of the β-globin protein), and the promoter features a TATA box consensus sequence like TATAAA approximately 25-30 base pairs upstream of the transcription start site to recruit RNA polymerase. DNA replication occurs via a semiconservative mechanism, in which each parental strand serves as a template for synthesizing a new complementary strand, resulting in two daughter molecules each containing one original and one newly synthesized strand. This process, experimentally demonstrated by Meselson and Stahl in 1958 using density-labeled DNA in E. coli, maintains sequence fidelity with an error rate of approximately 1 in 10^9 base pairs after proofreading and repair mechanisms. Mutations, such as point substitutions, insertions, or deletions, introduce variations in DNA sequences that drive evolutionary change by generating genetic diversity upon which natural selection acts. For example, a point mutation altering a single base can change an amino acid in a protein, while insertions or deletions may shift the reading frame, potentially leading to new traits or adaptations over time. In humans, the genome comprises about 3 billion base pairs across 23 pairs of chromosomes, encoding roughly 20,000 genes that collectively determine an individual's hereditary characteristics.

Functional Roles in RNA

RNA sequences play diverse functional roles in cellular processes, extending beyond their role as transcripts of DNA templates. These roles are often determined by specific sequence motifs that enable base-pairing, structural folding, and interactions with proteins or other nucleic acids. In eukaryotes and prokaryotes, RNA types such as (mRNA), (tRNA), (rRNA), and non-coding RNAs (ncRNAs) each exhibit sequence-dependent functions critical for and regulation. mRNA carries coding sequences from genes to for protein synthesis, with its untranslated regions (UTRs) containing regulatory sequences that influence stability, localization, and efficiency. For instance, in prokaryotes, the Shine-Dalgarno sequence—a purine-rich motif (typically AGGAGG) located 6-10 upstream of the —facilitates binding by base-pairing with the anti-Shine-Dalgarno sequence in 16S rRNA, enabling precise initiation. tRNA molecules feature anticodon sequences that base-pair with mRNA codons during , ensuring accurate incorporation; their cloverleaf secondary structure, formed by intramolecular base-pairing, is essential for recognition and function. rRNA forms the core of , where conserved sequence elements drive inter- and intramolecular base-pairing to create complex secondary and tertiary structures that position catalytic sites for formation. Many RNA functions rely on sequence-driven secondary structures, such as hairpins and loops, which arise from complementary base-pairing and modulate activity. Hairpins, consisting of a double-stranded stem and a single-stranded loop, are prevalent in ncRNAs like microRNAs (miRNAs), where the stem-loop structure is processed into mature miRNA for ; for example, the let-7 miRNA hairpin is recognized by LIN28A protein, inhibiting its maturation and thus regulating developmental timing. Loops, including apical or internal loops, often serve as binding sites for proteins or ligands, as seen in riboswitches—5' UTR sequences in bacterial mRNAs that fold into alternative conformations upon metabolite binding, thereby switching between terminator and antiterminator s to control transcription or translation. Post-transcriptional modifications like further diversify sequences and functions. Adenosine-to-inosine (A-to-I) editing, catalyzed by acting on () enzymes, targets double-stranded regions and is read as during , altering codons, splice sites, or miRNA targets to expand diversity and regulate innate immunity. In the , -mediated editing of subunits fine-tunes neuronal signaling, with editing levels varying by tissue and development stage. In viral RNA genomes, sequence features enable efficient replication and host interaction. Single-stranded RNA viruses like possess a ~30 kb positive-sense with open reading frames (ORFs) encoding structural and non-structural proteins; key features include a slippery sequence and in the ORF1ab region that induces -1 ribosomal frameshifting (15-60% efficiency), essential for producing the replicase polyprotein. The 's codon bias, favoring T/A-ending codons unlike the human host, optimizes viral while evading immune detection, and its 5' cap-like and 3' poly-A tail mimic host mRNAs for efficient expression. Regulatory roles of RNA sequences often involve ncRNAs acting as enhancers or silencers of . miRNAs, short ncRNAs (~22 nt), bind complementary sites in mRNA 3' UTRs via base-pairing, recruiting the (RISC) to repress or promote degradation, thereby fine-tuning gene networks in development and . In mRNA 3' UTRs, AU-rich elements ()—sequences like AUUUA repeats—act as silencers by accelerating deadenylation and decay, as in tumor necrosis factor-alpha (TNF-α) mRNA, where they limit inflammatory responses; binding proteins like tristetraprolin (TTP) mediate this instability. Conversely, stabilizing elements in UTRs, such as G-quadruplexes or stem-loops, can enhance expression by protecting against nucleases.

Sequencing Technologies

Historical Methods

Prior to the development of direct sequencing methods in the 1970s, determining nucleic acid sequences relied on indirect approaches, such as hybridization probes, which allowed partial characterization by detecting complementary base pairing between known oligonucleotide probes and target DNA or RNA under controlled conditions. These techniques, often used in conjunction with restriction enzyme mapping, provided limited insights into sequence motifs or restriction sites but could not yield complete linear orders of nucleotides due to their reliance on inference rather than direct readout. For instance, early hybridization experiments in the 1960s and early 1970s helped infer short RNA sequences, like those in transfer RNAs, by comparing melting temperatures and specificity of probe binding. The first direct DNA sequencing methods emerged in 1977 with the independent publications of the Maxam-Gilbert chemical cleavage technique and Sanger's chain-termination method, marking the onset of practical nucleic acid sequencing. In the Maxam-Gilbert approach, DNA is labeled at one end and subjected to base-specific chemical treatments—such as dimethyl sulfate for guanine, hydrazine for pyrimidines, or formic acid for adenine and guanine—to induce strand breaks at particular nucleotides, generating a set of fragments whose sizes are resolved by polyacrylamide gel electrophoresis to infer the sequence from band patterns. This method enabled the sequencing of up to several hundred base pairs but required hazardous chemicals and radioactive labeling, limiting its scalability. Concurrently, and colleagues introduced the chain-termination method, also known as dideoxy sequencing, which enzymatically synthesizes strands using in the presence of normal deoxynucleotides (dNTPs) and chain-terminating dideoxynucleotides (ddNTPs) that lack a 3'-hydroxyl group, halting extension at random positions corresponding to each base. The resulting fragments are separated by , with Sanger's plus-minus variant initially using differential incorporation to generate overlapping reads, allowing . This enzymatic method proved more reproducible and safer than chemical cleavage, facilitating the first complete sequence of the φX174, a 5,386-nucleotide single-stranded , achieved by Sanger's team in 1977. Early applications of these techniques included sequencing biologically significant genes, such as the rat insulin gene reported in , which demonstrated their utility in elucidating eukaryotic regulatory elements and coding sequences. However, both methods were labor-intensive, requiring manual gel pouring, radioactive isotope handling, and film autoradiography for detection, often taking days per run. Read lengths were typically limited to 100-500 base pairs, with error rates around 1-2% due to compression artifacts and ambiguous band resolution, particularly in repetitive regions, restricting analyses to small genomes or targeted fragments. The transition to began in the with the replacement of radioactive labels by fluorescent dyes attached to ddNTPs, enabling four-color detection in a single lane and machine-based readout, which increased throughput and reduced manual effort while paving the way for larger-scale projects.

Contemporary Techniques

Contemporary techniques in nucleic acid sequencing have revolutionized through high-throughput, massively parallel methods that enable rapid and cost-effective analysis of DNA and RNA sequences. Next-generation sequencing (NGS), also known as second-generation sequencing, relies on amplifying and sequencing millions of DNA fragments simultaneously, achieving throughput in the gigabase range per run. These platforms have democratized sequencing, facilitating applications from personalized medicine to . Key NGS platforms include Illumina's sequencing by synthesis, which uses reversible terminator to detect incorporated bases via during stepwise , producing short reads (typically 50-300 ) with high accuracy (>99.9%). Ion Torrent employs technology to measure changes from released ions during incorporation, offering faster turnaround times but with read lengths around 200-400 . (PacBio) utilizes single-molecule real-time (SMRT) sequencing, where a incorporates fluorescently labeled in zero-mode waveguides, generating long reads up to 20 kb or more, ideal for resolving structural variants; high-fidelity (HiFi) reads achieve >99% accuracy through circular consensus, though raw reads have higher error rates (around 10-15%). Third-generation sequencing advances single-molecule analysis without amplification, reducing biases and enabling real-time data generation. Oxford Nanopore Technologies' platform sequences DNA or RNA by passing molecules through protein nanopores, detecting ionic current disruptions as bases translocate, which allows portable, real-time analysis with read lengths exceeding 1 Mb. Accuracy has reached >99% single-read and consensus levels as of 2025 through advanced basecalling algorithms and chemistry updates like R10.4.1, addressing early limitations in homopolymer resolution. These techniques support diverse applications, such as whole-genome sequencing, where the cost of sequencing fell below $1,000 by 2020 and further declined to approximately $200-$600 as of 2025, enabling population-scale studies. profiles microbial communities directly from environmental samples, while single-cell sequencing (scRNA-seq) captures transcriptomes from individual cells, revealing cellular heterogeneity in development and . Challenges persist in error correction, particularly for repetitive regions that confound short-read assembly; algorithms like those in the Canu assembler integrate long-read data to achieve near-complete . Recent advances as of 2025 include ultra-high-throughput platforms like Ultima Genomics UG 100, capable of sequencing at under $100 per , and novel spatial methods such as expansion genome sequencing for mapping DNA relative to cellular structures. CRISPR-based nucleic acid detection methods, such as SHERLOCK (using Cas13) and adaptations with CRISPR-Cas12a developed in the late and , amplify and detect specific sequences with high sensitivity for diagnostics, complementing sequencing in point-of-care applications. For RNA, direct sequencing without cDNA conversion preserves modifications like m6A, using platforms like Oxford Nanopore to sequence native strands and reveal epitranscriptomic features.

Digital Handling

Data Formats

Nucleic acid sequence data is stored and exchanged using standardized text-based and binary formats that encode sequences, metadata, and quality information for computational processing. These formats facilitate interoperability across sequencing platforms and analysis tools, building on basic notation systems for nucleotides. The is a simple, human-readable text format consisting of a header line beginning with a greater-than (>) followed by a sequence identifier and optional description, and subsequent lines containing the or sequence, typically limited to 60-80 characters per line for readability. It supports basic sequence representation without quality scores, though a variant called QUAL extends it by pairing a FASTA-like sequence file with a corresponding quality file. The builds on by incorporating per-base scores alongside , structured as four lines per record: a header starting with @, , a separator line with +, and a of equal to . scores are encoded using Phred values, where Q=−10log⁡10PQ = -10 \log_{10} P and PP is the estimated probability of an incorrect base call, allowing assessment of sequencing accuracy. This format originated at the Sanger Institute for sequencing and has variants for next-generation platforms like Illumina. For aligned sequence reads, the Sequence Alignment/Map (SAM) format provides a tab-delimited text structure detailing mappings to a , including fields for read name, flags (bitwise integers indicating properties like paired-end status), , position, mapping quality, and optional tags for additional metadata. Its binary counterpart, BAM, compresses SAM data losslessly for efficient storage and indexing, supporting flags that denote paired-end alignments and other read attributes. Specialized formats enrich sequences with annotations: uses a flat-file structure with sections for locus, definition, features (e.g., genes, exons), and the sequence itself, enabling detailed biological context like source organism and publication references. The General Feature Format (GFF), particularly GFF3, employs a nine-column tab-delimited layout per feature line to specify genomic elements such as genes or regulatory regions, with columns for sequence ID, source, type, coordinates, score, strand, phase, and attributes. To address growing data volumes, formats have evolved toward compression; for instance, CRAM (Compressed Reference-oriented Alignment Map) refines BAM by leveraging dependencies for lossy or lossless encoding, achieving typical file size reductions of 30-50% over BAM while maintaining compatibility with SAM tools. Broader compression techniques, including reference-based and , further reduce genomic dataset sizes by 50-90% in specialized implementations, balancing efficiency with accessibility for large-scale analyses.

Storage and Databases

The primary repositories for nucleic acid sequences are maintained through the International Nucleotide Sequence Database Collaboration (INSDC), a longstanding partnership among (hosted by the , NCBI, in the United States, established in 1982), the European Nucleotide Archive (ENA, managed by the European Molecular Biology Laboratory's , EMBL-EBI), and the DNA Data Bank of Japan (DDBJ). These organizations synchronize their data daily to provide a unified, non-redundant view of global nucleotide sequence submissions, encompassing raw reads, assemblies, and annotations from diverse sources including genomic projects and research submissions. GenBank sequences are distributed in flat file formats that include detailed annotations such as locus identifiers, features, , and bibliographic references, enabling comprehensive metadata alongside the primary data. Comprehensive releases occur bimonthly, with daily incremental updates available via FTP to reflect ongoing submissions and ensure timely access. As of August 2025, alone holds over 47 trillion base pairs across nearly 6 billion records, reflecting driven by advances in sequencing technologies. Specialized repositories complement these core databases by focusing on niche aspects of nucleic acid data. RNAcentral serves as a centralized hub for sequences, aggregating data from 52 expert databases to provide unified access to ncRNA types such as miRNAs, lncRNAs, and snoRNAs across organisms. The , an international effort, maintains a detailed catalog of , including millions of single nucleotide polymorphisms (SNPs) and structural variants derived from sequencing over 2,500 individuals, supporting and disease association studies. Managing these repositories faces significant challenges due to the explosive growth of sequence data, projected to require up to 40 exabytes of storage capacity by 2025 for human genomics alone. concerns are paramount in human-related datasets, with regulations like the European Union's (GDPR) mandating strict controls on data sharing, consent, and re-identification risks to prevent misuse of sensitive genetic information. Versioning systems are essential to track updates and revisions in sequence records, ensuring reproducibility while handling the complexity of iterative assemblies and annotations. Access to these databases is facilitated by tools like NCBI's system, which integrates , protein, and genomic data for cross-database searching and retrieval. The Basic Local Alignment Search Tool (BLAST) enables rapid similarity searches against these repositories, supporting tasks from homology detection to functional . In the 2020s, integrations with and ontologies, such as the Sequence Alignment Ontology (SALON), have enhanced querying by enabling semantic searches and automated interpretation of alignments and metadata.

Analytical Approaches

Sequence Alignment

Sequence alignment is a fundamental computational technique in bioinformatics used to identify regions of similarity between sequences, which can indicate functional, structural, or evolutionary relationships. By comparing two or more sequences, alignments reveal conserved regions, insertions, deletions (indels), and substitutions, aiding in the inference of biological processes such as function prediction and phylogenetic . Pairwise sequence alignment compares two sequences to find the optimal arrangement that maximizes similarity. The Needleman-Wunsch algorithm, introduced in , performs global alignment by aligning entire sequences using dynamic programming, ensuring that the full length of both sequences is considered from end to end. This method constructs a scoring matrix where each cell represents the best alignment score up to that position, to recover the alignment path. For nucleotide sequences, it employs a to score matches and mismatches; a common simple scheme assigns +1 for identical bases and -1 for differences, though more sophisticated matrices like the NUC.4.4, derived from observed substitutions, can be used to account for transition/transversion biases. In contrast, the Smith-Waterman algorithm, developed in 1981, focuses on local alignment to detect high-similarity regions within longer sequences, which is particularly useful for identifying conserved domains in divergent nucleic acids. It modifies the Needleman-Wunsch approach by initializing the matrix with zeros and setting negative scores to zero, preventing penalties from propagating across unrelated regions. Both algorithms incorporate gap penalties to handle indels: linear penalties charge a constant cost (-d) per gap position, while affine penalties, introduced by Gotoh in 1982, distinguish gap opening (-a) from extension (-(g-1)d), better modeling biological insertion/deletion events by penalizing starts more heavily than continuations. The total alignment score can be expressed as: S=∑s(xi,yi)+∑g(k)S = \sum s(x_i, y_i) + \sum g(k) where s(xi,yi)s(x_i, y_i) is the substitution score for aligned positions, and g(k)g(k) is the gap penalty for each gap of length kk, typically negative. For comparing more than two sequences, multiple sequence alignment (MSA) extends pairwise methods to reveal patterns across a set. Progressive alignment strategies, a cornerstone of MSA, build alignments iteratively: first, a distance matrix is computed from pairwise scores, a guide tree is constructed via hierarchical clustering, and sequences are then aligned following the tree branches, starting with the most similar pairs. Clustal Omega, released in 2011, implements this approach with enhancements like mBed for large-scale alignments, enabling rapid processing of thousands of sequences while maintaining accuracy comparable to slower methods. Similarly, MAFFT, first described in 2002, uses fast Fourier transform to approximate distance calculations, accelerating progressive alignment and supporting iterative refinement for improved handling of divergent sequences. These alignments have key applications in detecting homology between sequences, where significant similarity suggests shared ancestry, and in constructing evolutionary trees by using alignment scores or distances to infer phylogenies. They are essential for managing indels in highly divergent sequences, as gap models allow flexible insertions without overly disrupting conserved regions. In the era of next-generation sequencing (NGS), particularly with long-read technologies like PacBio and Oxford Nanopore, alignment methods have evolved to accommodate error-prone, lengthy reads; tools such as Minimap2 employ seed-and-extend heuristics with affine gaps to map these efficiently against reference genomes, addressing challenges like structural variants that short-read aligners struggle with.

Motif Detection

Sequence motifs are short, recurring patterns in DNA or RNA sequences, typically 6-20 base pairs long, that often indicate functional elements due to their conservation across related sequences. For example, the in eukaryotic promoters, with the TATAAA, serves as a for the to initiate transcription. These motifs can vary slightly in sequence but maintain functional significance through evolutionary conservation. Motifs are classified by their biological roles, including regulatory motifs in enhancers that control gene expression, structural motifs such as ribosome binding sites (RBS) in mRNA, and protein-binding motifs like transcription factor binding sites (TFBS). Regulatory motifs, such as those in enhancers, recruit transcription factors to modulate gene activity in specific cellular contexts. Structural motifs, prominent in RNA, include the Shine-Dalgarno sequence (e.g., AGGAGG) upstream of start codons in prokaryotic mRNA, which facilitates ribosome assembly for translation initiation. Protein-binding motifs encompass TFBS in DNA, where sequence patterns enable specific interactions with regulatory proteins, and analogous sites in RNA for RNA-binding proteins. RNA motifs, often underemphasized, play critical roles in processes like splicing and RNA stability, with examples including internal ribosome entry sites (IRES) that direct cap-independent translation. Key tools for motif detection include (Multiple EM for Motif Elicitation), which uses expectation maximization to discover ungapped motifs in unaligned sequences by modeling them as position-specific scoring matrices. provides a database of documented patterns and profiles for identifying functional motifs in protein-coding nucleic acid sequences, aiding in the annotation of domains and sites. For motifs with positional variability, position weight matrices (PWMs) represent the likelihood of each nucleotide at every position, derived from aligned sequences. The score for a candidate sequence is calculated as the sum over positions jj of log⁡2(fj,bbb)\log_2 \left( \frac{f_{j,b}}{b_b} \right), where fj,bf_{j,b} is the observed frequency of base bb at position jj, and bbb_b is the background frequency; higher scores indicate better matches. Motif detection enables applications in predicting gene regulation, where identified patterns forecast enhancer activity or TFBS occupancy to model expression dynamics. It also supports functional annotation by linking motifs to biological roles, such as classifying non-coding regions as regulatory elements. Advances in machine learning, starting with DeepBind in 2015, employ convolutional neural networks to predict protein-DNA and protein-RNA binding specificities from sequence data, outperforming traditional PWM-based methods on large datasets. Post-2020 integrations of deep learning, including transformer-based models, have enhanced motif discovery in RNA contexts by incorporating structural features and improving accuracy in high-throughput data like CLIP-seq, addressing gaps in earlier approaches.

Complexity Measures

Complexity measures in nucleic acid sequences quantify the variability, randomness, and information content inherent in DNA or RNA strings, providing insights into their structural and functional properties independent of relational comparisons like alignments. These metrics assess how unpredictable or repetitive a sequence is, which correlates with biological constraints such as evolutionary pressures or mutational patterns. For instance, highly random sequences approach maximum entropy, indicating minimal redundancy, while repetitive or biased ones exhibit lower complexity, often reflecting functional adaptations. A primary measure is Shannon entropy, which evaluates the uncertainty or information content per base in a sequence. Defined as H=−∑pilog⁡2piH = -\sum p_i \log_2 p_i, where pip_i is the frequency of each base (A, C, G, T/U), it ranges from 0 bits per base for a fully repetitive sequence to 2 bits per base for a uniformly random one with equal base frequencies. This metric, adapted from , highlights ; for example, coding regions in genomes often show values around 1.8-1.9 bits due to codon usage biases, while non-coding repeats drop below 1 bit. Shannon entropy can be computed from base frequencies derived briefly from aligned sequences but applies to individual sequences as well. Other metrics complement by capturing different aspects of repetitiveness and . Lempel-Ziv approximates by counting distinct substrings in a during compression-like , yielding a normalized score between 0 (purely repetitive) and 1 (incompressible ); it is particularly useful for identifying low- regions in genomes, such as tandem repeats, where scores below 0.3 indicate high repetitiveness. , measured as the deviation from 50% guanine-cytosine proportion (e.g., via ∣GC%−50∣|\text{GC\%} - 50|), influences perceived by skewing base distributions, with extreme biases (e.g., >70% GC in vertebrate CpG islands) reducing and affecting evolutionary analyses. k-mer diversity assesses repetitiveness by counting unique substrings of length k (typically 3-6 bases), where lower diversity signals tandem repeats or segmental duplications, as observed in repetitive regions of eukaryotic genomes such as . These measures find applications in distinguishing functional genomic elements and analyzing . In distinguishing coding from non-coding regions, low and Lempel-Ziv scores (e.g., <0.4) mark non-coding areas with repeats, while higher values (~1.9 bits) typify protein-coding exons under selective pressure for diversity. For —clouds of mutant variants—Shannon quantifies intra-host diversity, with higher values indicating substantial mutation rates in variable regions. Compression efficiency, tied to Lempel-Ziv, optimizes storage of repetitive genomes. Biologically, low often arises in regulatory regions due to physicochemical constraints, such as secondary stability in promoters, where base pairing reduces variability to ~1.2 bits compared to 1.8 in intergenic spacers. Tools like the Entropy-One calculator from the Los Alamos HIV Database compute site-specific Shannon for aligned sequences, facilitating variability analysis in viral datasets. In , profiles reveal community diversity, with recent studies using vectors to encode microbial sequences for efficient assembly, capturing up to 95% of variability in gut microbiomes. Emerging AI models, such as those predicting long-range dependencies, infer from raw data, enhancing detection of cryptic regulatory motifs missed by traditional metrics.

References

  1. https://doi.org/10.1042/EBC20180038
  2. https://www.vanderbilt.edu/AnS/Chemistry/Rizzo/Chem220b/Ch28.pdf
  3. https://www.genome.gov/about-genomics/fact-sheets/DNA-Sequencing-Fact-Sheet
  4. https://www.genome.gov/genetics-glossary/Nucleotide
  5. https://pmc.ncbi.nlm.nih.gov/articles/PMC6822018/
  6. https://atdbio.com/nucleic-acids-book/Nucleic-acid-structure
  7. https://www.khanacademy.org/test-prep/mcat/biomolecules/dna/a/nucleic-acids
  8. https://www.genome.gov/genetics-glossary/Phosphate-Backbone
  9. https://pmc.ncbi.nlm.nih.gov/articles/PMC4615568/
  10. https://pmc.ncbi.nlm.nih.gov/articles/PMC7686126/
  11. https://iubmb.qmul.ac.uk/misc/naseq.html
  12. https://medlineplus.gov/genetics/understanding/basics/dna/
  13. https://www.genome.gov/about-genomics/fact-sheets/Deoxyribonucleic-Acid-Fact-Sheet
  14. https://www.ncbi.nlm.nih.gov/books/NBK9944/
  15. https://www.genome.gov/genetics-glossary/Central-Dogma
  16. https://www.dna.caltech.edu/courses/cs191/paperscs191/CrickCentralDogma1970.pdf
  17. https://www.ncbi.nlm.nih.gov/books/NBK10023/
  18. https://www.ncbi.nlm.nih.gov/books/NBK26887/
  19. https://pmc.ncbi.nlm.nih.gov/articles/PMC1557599/
  20. https://pmc.ncbi.nlm.nih.gov/articles/PMC528642/
  21. http://www.nature.com/scitable/topicpage/genetic-mutation-1127
  22. https://pmc.ncbi.nlm.nih.gov/articles/PMC4563715/
  23. https://pmc.ncbi.nlm.nih.gov/articles/PMC4829331/
  24. https://pmc.ncbi.nlm.nih.gov/articles/PMC7656661/
  25. https://www.nature.com/scitable/topicpage/riboswitches-a-common-rna-regulatory-element-14262702/
  26. https://pmc.ncbi.nlm.nih.gov/articles/PMC5731795/
  27. https://www.nature.com/articles/s41598-020-72533-2
  28. https://pmc.ncbi.nlm.nih.gov/articles/PMC1325018/
  29. https://pmc.ncbi.nlm.nih.gov/articles/PMC3474909/
  30. https://www.nature.com/articles/265687a0
  31. https://www.science.org/doi/10.1126/science.325648
  32. https://www.pacb.com/blog/long-read-sequencing-myths-debunked-part-1-hifi-sequencing/
  33. https://nanoporetech.com/platform/accuracy
  34. https://www.genome.gov/about-genomics/fact-sheets/DNA-Sequencing-Costs-Data
  35. https://frontlinegenomics.com/the-100-genome-wheres-the-limit/
  36. https://www.broadinstitute.org/news/new-technique-expands-cells-sequence-dna-and-capture-fine-structural-details
  37. https://www.nature.com/articles/s41596-019-0210-2
  38. https://nanoporetech.com/resource-centre/workflow-direct-rna-sequencing
  39. https://blast.ncbi.nlm.nih.gov/doc/blast-topics/
  40. https://pmc.ncbi.nlm.nih.gov/articles/PMC2847217/
  41. https://pmc.ncbi.nlm.nih.gov/articles/PMC2723002/
  42. https://samtools.github.io/hts-specs/SAMv1.pdf
  43. https://www.ncbi.nlm.nih.gov/genbank/samplerecord/
  44. http://www.ensembl.org/info/website/upload/gff3.html
  45. https://samtools.github.io/hts-specs/CRAMv3.pdf
  46. https://pmc.ncbi.nlm.nih.gov/articles/PMC8896640/
  47. https://www.insdc.org/
  48. https://www.ncbi.nlm.nih.gov/genbank/collab
  49. https://pubmed.ncbi.nlm.nih.gov/33166387/
  50. https://academic.oup.com/nar/article/53/D1/D62/7899524
  51. https://www.ncbi.nlm.nih.gov/genbank/release/current
  52. https://academic.oup.com/nar/article/53/D1/D56/7903376
  53. https://www.ncbi.nlm.nih.gov/genbank/
  54. https://ncbiinsights.ncbi.nlm.nih.gov/2025/08/26/genbank-release-268-0/
  55. https://rnacentral.org/
  56. https://www.biorxiv.org/content/10.1101/2025.09.19.677274v1
  57. http://www.internationalgenome.org/
  58. http://www.internationalgenome.org/data/
  59. https://medicalfuturist.com/the-genomic-data-challenges-of-the-future/
  60. https://law.illinois.edu/wp-content/uploads/2025/10/Sherkow-2025-Asociotechnicalapproachtogenomicdata.pdf
  61. https://pmc.ncbi.nlm.nih.gov/articles/PMC12573770/
  62. https://www.ncbi.nlm.nih.gov/Web/Search/entrezfs.html
  63. https://blast.ncbi.nlm.nih.gov/Blast.cgi
  64. https://bmcbioinformatics.biomedcentral.com/articles/10.1186/s12859-023-05190-7
  65. https://eprints.whiterose.ac.uk/id/eprint/229031/1/BIB-25-0517.R1_Proof_hi.pdf
  66. https://pubmed.ncbi.nlm.nih.gov/5420325/
  67. https://pubmed.ncbi.nlm.nih.gov/21988835/
  68. https://genomebiology.biomedcentral.com/articles/10.1186/s13059-023-02972-3
  69. https://open.oregonstate.education/appliedbioinformatics/chapter/chapter-2-sequence-motifs/
  70. https://pmc.ncbi.nlm.nih.gov/articles/PMC85905/
  71. https://www.ncbi.nlm.nih.gov/books/NBK26806/
  72. https://academic.oup.com/mbe/article/42/10/msaf237/8283777
  73. https://www.nature.com/articles/s41586-020-2077-3
  74. https://academic.oup.com/nar/article/34/suppl_2/W369/2505578
  75. https://europepmc.org/articles/pmc3531220
  76. https://pmc.ncbi.nlm.nih.gov/articles/PMC3820676/
  77. https://www.nature.com/articles/nbt.3300
  78. https://www.nature.com/articles/s12276-024-01243-w
  79. https://pmc.ncbi.nlm.nih.gov/articles/PMC3177615/
  80. https://pmc.ncbi.nlm.nih.gov/articles/PMC10124876/
  81. https://pmc.ncbi.nlm.nih.gov/articles/PMC441604/
  82. https://www.frontiersin.org/journals/genetics/articles/10.3389/fgene.2017.00016/full
  83. https://pmc.ncbi.nlm.nih.gov/articles/PMC11874746/
  84. https://pmc.ncbi.nlm.nih.gov/articles/PMC3358654/
  85. https://www.hiv.lanl.gov/content/sequence/ENTROPY/entropy_one.html
  86. https://pubmed.ncbi.nlm.nih.gov/40914969/
  87. https://deepmind.google/blog/alphagenome-ai-for-better-understanding-the-genome/
Add your contribution
Related Hubs
User Avatar
No comments yet.