Langbahn Team – Weltmeisterschaft

Translation (biology)

Overview of eukaryotic messenger RNA (mRNA) translation
Translation of mRNA and ribosomal protein synthesis
Initiation and elongation stages of translation involving RNA nucleobases, the ribosome, transfer RNA, and amino acids
The three phases of translation: (1) in initiation, the small ribosomal subunit binds to the RNA strand and the initiator tRNA–amino acid complex binds to the start codon, culminating in attachment of the large subunit; (2) elongation occurs as a cycle, in which codons are sequentially recognized by charged tRNAs, followed by peptide bond formation with transfer of the polypeptide between tRNAs within the ribosome and finally translocation of the ribosome to the next codon; (3) termination, when a stop codon is reached, a release factor binds and the polypeptide is released (note that labels for translocation and peptide bond formation are reversed in this image)

In biology, translation is the process in living cells in which proteins are produced using RNA molecules as templates. The generated protein is a sequence of amino acids. This sequence is determined by the sequence of nucleotides in the RNA. The nucleotides are considered three at a time. Each such triple results in addition of one specific amino acid to the protein being generated. The matching from nucleotide triple to amino acid is called the genetic code. The translation is performed by a large complex of functional RNA and proteins called ribosomes. The entire process is called gene expression.

In translation, messenger RNA (mRNA) is decoded in a ribosome, outside the nucleus, to produce a specific amino acid chain, or polypeptide. The polypeptide later folds into an active protein and performs its functions in the cell. The ribosome facilitates decoding by inducing the binding of complementary transfer RNA (tRNA) anticodon sequences to mRNA codons. The tRNAs carry specific amino acids that are chained together into a polypeptide as the mRNA passes through and is "read" by the ribosome.

Translation proceeds in three phases:

  1. Initiation: The ribosome assembles around the target mRNA. The first tRNA is attached at the start codon.
  2. Elongation: The last tRNA validated by the small ribosomal subunit (accommodation) transfers the amino acid. It carries to the large ribosomal subunit which binds it to one of the preceding admitted tRNA (transpeptidation). The ribosome then moves to the next mRNA codon to continue the process (translocation), creating an amino acid chain.
  3. Termination: When a stop codon is reached, the ribosome releases the polypeptide. The ribosomal complex remains intact and moves on to the next mRNA to be translated.

In prokaryotes (bacteria and archaea), translation occurs in the cytosol, where the large and small subunits of the ribosome bind to the mRNA. In eukaryotes, translation occurs in the cytoplasm or across the membrane of the endoplasmic reticulum in a process called co-translational translocation. In co-translational translocation, the entire ribosome/mRNA complex binds to the outer membrane of the rough endoplasmic reticulum (ER), and the new protein is synthesized and released into the ER; the newly created polypeptide can be stored inside the ER for future vesicle transport and secretion outside the cell, or immediately secreted.

Many types of transcribed RNA, such as tRNA, ribosomal RNA, and small nuclear RNA, do not undergo a translation into proteins.

Several antibiotics act by inhibiting translation. These include anisomycin, cycloheximide, chloramphenicol, tetracycline, streptomycin, erythromycin, and puromycin. Prokaryotic ribosomes have a different structure from that of eukaryotic ribosomes, and thus antibiotics can specifically target bacterial infections without any harm to a eukaryotic host's cells.

Basic mechanisms

A ribosome translating a protein that is secreted into the endoplasmic reticulum (tRNAs colored dark blue)
Tertiary structure of tRNA (CCA tail in yellow, Acceptor stem in purple, Variable loop in orange, D arm in red, Anticodon arm in blue with Anticodon in black, T arm in green)

The basic process of protein production is the addition of one amino acid at a time to the end of a protein. This operation is performed by a ribosome.[1] A ribosome is made up of two subunits, a small subunit, and a large subunit. These subunits come together before the translation of mRNA into a protein to provide a location for translation to be carried out and a polypeptide to be produced.[2] The choice of amino acid type to add is determined by a messenger RNA (mRNA) molecule. Each amino acid added is matched to a three-nucleotide subsequence of the mRNA. For each such triplet possible, the corresponding amino acid is accepted. The successive amino acids added to the chain are matched to successive nucleotide triplets in the mRNA. In this way, the sequence of nucleotides in the template mRNA chain determines the sequence of amino acids in the generated amino acid chain.[3] The addition of an amino acid occurs at the C-terminus of the peptide; thus, translation is said to be amine-to-carboxyl directed.[4]

The mRNA carries genetic information encoded as a ribonucleotide sequence from the chromosomes to the ribosomes. The ribonucleotides are "read" by translational machinery in a sequence of nucleotide triplets called codons. Each of those triplets codes for a specific amino acid.[citation needed]

The ribosome molecules translate this code to a specific sequence of amino acids. The ribosome is a multisubunit structure containing ribosomal RNA (rRNA) and proteins. It is the "factory" where amino acids are assembled into proteins.

Transfer RNAs (tRNAs) are small noncoding RNA chains (74–93 nucleotides) that transport amino acids to the ribosome. The repertoire of tRNA genes varies widely between species, with some bacteria having between 20 and 30 genes while complex eukaryotes could have thousands.[5] tRNAs have a site for amino acid attachment, and a site called an anticodon. The anticodon is an RNA triplet complementary to the mRNA triplet that codes for their cargo amino acid.

Aminoacyl tRNA synthetases (enzymes) catalyze the bonding between specific tRNAs and the amino acids that their anticodon sequences call for. The product of this reaction is an aminoacyl-tRNA. The amino acid is joined by its carboxyl group to the 3' OH of the tRNA by an ester bond. When the tRNA has an amino acid linked to it, the tRNA is termed "charged". In bacteria, this aminoacyl-tRNA is carried to the ribosome by EF-Tu, where mRNA codons are matched through complementary base pairing to specific tRNA anticodons. Aminoacyl-tRNA synthetases that mispair tRNAs with the wrong amino acids can produce mischarged aminoacyl-tRNAs, which can result in inappropriate amino acids at the respective position in the protein. This "mistranslation"[6] of the genetic code naturally occurs at low levels in most organisms, but certain cellular environments cause an increase in permissive mRNA decoding, sometimes to the benefit of the cell.

The ribosome has two binding sites for tRNA. They are the aminoacyl site (abbreviated A), and the peptidyl site/ exit site (abbreviated P/E). Concerning the mRNA, the three sites are oriented 5' to 3' E-P-A, because ribosomes move toward the 3' end of mRNA. The A-site binds the incoming tRNA with the complementary codon on the mRNA. The P/E-site holds the tRNA with the growing polypeptide chain. When an aminoacyl-tRNA initially binds to its corresponding codon on the mRNA, it is in the A site. Then, a peptide bond forms between the amino acid of the tRNA in the A site and the amino acid of the charged tRNA in the P/E site. The growing polypeptide chain is transferred to the tRNA in the A site. Translocation occurs, moving the tRNA to the P/E site, now without an amino acid; the tRNA that was in the A site, now charged with the polypeptide chain, is moved to the P/E site and the uncharged tRNA leaves, and another aminoacyl-tRNA enters the A site to repeat the process.[7]

After the new amino acid is added to the chain, and after the tRNA is released out of the ribosome and into the cytosol, the energy provided by the hydrolysis of a GTP bound to the translocase EF-G (in bacteria) and a/eEF-2 (in eukaryotes and archaea) moves the ribosome down one codon towards the 3' end. The energy required for translation of proteins is significant. For a protein containing n amino acids, the number of high-energy phosphate bonds required to translate it is 4n-1.[8] The rate of translation varies; it is significantly higher in prokaryotic cells (up to 17–21 amino acid residues per second) than in eukaryotic cells (up to 6–9 amino acid residues per second).[9]

Initiation and termination of translation

Initiation involves the small subunit of the ribosome binding to the 5' end of mRNA with the help of initiation factors (IF). In bacteria and a minority of archaea, initiation of protein synthesis involves the recognition of a purine-rich initiation sequence on the mRNA called the Shine–Dalgarno sequence. The Shine–Dalgarno sequence binds to a complementary pyrimidine-rich sequence on the 3' end of the 16S rRNA part of the 30S ribosomal subunit. The binding of these complementary sequences ensures that the 30S ribosomal subunit is bound to the mRNA and is aligned such that the initiation codon is placed in the 30S portion of the P-site. Once the mRNA and 30S subunit are properly bound, an initiation factor brings the initiator tRNA–amino acid complex, f-Met-tRNA, to the 30S P site. The initiation phase is completed once a 50S subunit joins the 30S subunit, forming an active 70S ribosome.[10] Termination of the polypeptide occurs when the A site of the ribosome is occupied by a stop codon (UAA, UAG, or UGA) on the mRNA, creating the primary structure of a protein. tRNA usually cannot recognize or bind to stop codons. Instead, the stop codon induces the binding of a release factor protein[11] (RF1 & RF2) that prompts the disassembly of the entire ribosome/mRNA complex by the hydrolysis of the polypeptide chain from the peptidyl transferase center [1] of the ribosome.[12] Drugs or special sequence motifs on the mRNA can change the ribosomal structure so that near-cognate tRNAs are bound to the stop codon instead of the release factors. In such cases of 'translational readthrough', translation continues until the ribosome encounters the next stop codon.[13]

Errors in translation

Even though the ribosomes are usually considered accurate and processive machines, the translation process is subject to errors that can lead either to the synthesis of erroneous proteins or to the premature abandonment of translation, either because a tRNA couples to a wrong codon or because a tRNA is coupled to the wrong amino acid. [14] The rate of error in synthesizing proteins has been estimated to be between 1 in 105 and 1 in 103 misincorporated amino acids, depending on the experimental conditions.[15] The rate of premature translation abandonment, instead, has been estimated to be of the order of magnitude of 10−4 events per translated codon.[16][17]

Regulation

The process of translation is highly regulated in both eukaryotic and prokaryotic organisms. Regulation of translation can impact the global rate of protein synthesis which is closely coupled to the metabolic and proliferative state of a cell.

To delve deeper into this intricate process, scientists typically use a technique known as ribosome profiling.[18] This method enables researchers to take a snapshot of the translatome, showing which parts of the mRNA are being translated into proteins by ribosomes at a given time. Ribosome profiling provides valuable insights into translation dynamics, revealing the complex interplay between gene sequence, mRNA structure, and translation regulation. For example, research utilizing this method has revealed that genetic differences and their subsequent expression as mRNAs can also impact translation rate in an RNA-specific manner.[19]

Expanding on this concept, a more recent development is single-cell ribosome profiling, a technique that allows us to study the translation process at the resolution of individual cells.[20] This is particularly significant as cells, even those of the same type, can exhibit considerable variability in their protein synthesis. Single-cell ribosome profiling has the potential to shed light on the heterogeneous nature of cells, leading to a more nuanced understanding of how translation regulation can impact cell behavior, metabolic state, and responsiveness to various stimuli or conditions.

Clinical significance

Translational control is critical for the development and survival of cancer. Cancer cells must frequently regulate the translation phase of gene expression, though it is not fully understood why translation is targeted over steps like transcription. While cancer cells often have genetically altered translation factors, it is much more common for cancer cells to modify the levels of existing translation factors.[21] Several major oncogenic signaling pathways, including the RAS–MAPK, PI3K/AKT/mTOR, MYC, and WNT–β-catenin pathways, ultimately reprogram the genome via translation.[22] Cancer cells also control translation to adapt to cellular stress. During stress, the cell translates mRNAs that can mitigate the stress and promote survival. An example of this is the expression of AMPK in various cancers; its activation triggers a cascade that can ultimately allow the cancer to escape apoptosis (programmed cell death) triggered by nutrition deprivation. Future cancer therapies may involve disrupting the translation machinery of the cell to counter the downstream effects of cancer.[21]

Mathematical modeling of translation

Figure M0. Basic and the simplest model M0 of protein synthesis. Here, * M – amount of mRNA with translation initiation site not occupied by assembling ribosome, * F – amount of mRNA with translation initiation site occupied by assembling ribosome, * R – amount of ribosomes sitting on mRNA synthesizing proteins, * P – amount of synthesized proteins.[23]
Figure M1'. The extended model of protein synthesis M1 with explicit presentation of 40S, 60S and initiation factors (IF) binding.[23]

The transcription-translation process description, mentioning only the most basic "elementary" processes, consists of:

  1. production of mRNA molecules (including splicing),
  2. initiation of these molecules with help of initiation factors (e.g., the initiation can include the circularization step though it is not universally required),
  3. initiation of translation, recruiting the small ribosomal subunit,
  4. assembly of full ribosomes,
  5. elongation, (i.e. movement of ribosomes along mRNA with production of protein),
  6. termination of translation,
  7. degradation of mRNA molecules,
  8. degradation of proteins.

The process of amino acid building to create protein in translation is a subject of various physic models for a long time starting from the first detailed kinetic models such as[24] or others taking into account stochastic aspects of translation and using computer simulations. Many chemical kinetics-based models of protein synthesis have been developed and analyzed in the last four decades.[25][26] Beyond chemical kinetics, various modeling formalisms such as Totally Asymmetric Simple Exclusion Process,[26] Probabilistic Boolean Networks, Petri Nets and max-plus algebra have been applied to model the detailed kinetics of protein synthesis or some of its stages. A basic model of protein synthesis that takes into account all eight 'elementary' processes has been developed,[23] following the paradigm that "useful models are simple and extendable".[27] The simplest model M0 is represented by the reaction kinetic mechanism (Figure M0). It was generalised to include 40S, 60S and initiation factors (IF) binding (Figure M1'). It was extended further to include effect of microRNA on protein synthesis.[28] Most of models in this hierarchy can be solved analytically. These solutions were used to extract 'kinetic signatures' of different specific mechanisms of synthesis regulation.

Genetic code

It is also possible to translate either by hand (for short sequences) or by computer (after first programming one appropriately, see section below); this allows biologists and chemists to draw out the chemical structure of the encoded protein on paper.

First, convert each template DNA base to its RNA complement (note that the complement of A is now U), as shown below. Note that the template strand of the DNA is the one the RNA is polymerized against; the other DNA strand would be the same as the RNA, but with thymine instead of uracil.

DNA -> RNA
 A  ->  U
 T  ->  A
 C  ->  G
 G  ->  C
 A=T-> A=U

Then split the RNA into triplets (groups of three bases). Note that there are 3 translation "windows", or reading frames, depending on where you start reading the code. Finally, use the table at Genetic code to translate the above into a structural formula as used in chemistry.

This will give the primary structure of the protein. However, proteins tend to fold, depending in part on hydrophilic and hydrophobic segments along the chain. Secondary structure can often still be guessed at, but the proper tertiary structure is often very hard to determine.

Whereas other aspects such as the 3D structure, called tertiary structure, of protein can only be predicted using sophisticated algorithms, the amino acid sequence, called primary structure, can be determined solely from the nucleic acid sequence with the aid of a translation table.

This approach may not give the correct amino acid composition of the protein, in particular if unconventional amino acids such as selenocysteine are incorporated into the protein, which is coded for by a conventional stop codon in combination with a downstream hairpin (SElenoCysteine Insertion Sequence, or SECIS).

There are many computer programs capable of translating a DNA/RNA sequence into a protein sequence. Normally this is performed using the Standard Genetic Code, however, few programs can handle all the "special" cases, such as the use of the alternative initiation codons which are biologically significant. For instance, the rare alternative start codon CTG codes for Methionine when used as a start codon, and for Leucine in all other positions.

Example: Condensed translation table for the Standard Genetic Code (from the NCBI Taxonomy webpage).[29]

 AAs    = FFLLSSSSYY**CC*WLLLLPPPPHHQQRRRRIIIMTTTTNNKKSSRRVVVVAAAADDEEGGGG
 Starts = ---M---------------M---------------M----------------------------
 Base1  = TTTTTTTTTTTTTTTTCCCCCCCCCCCCCCCCAAAAAAAAAAAAAAAAGGGGGGGGGGGGGGGG
 Base2  = TTTTCCCCAAAAGGGGTTTTCCCCAAAAGGGGTTTTCCCCAAAAGGGGTTTTCCCCAAAAGGGG
 Base3  = TCAGTCAGTCAGTCAGTCAGTCAGTCAGTCAGTCAGTCAGTCAGTCAGTCAGTCAGTCAGTCAG

The "Starts" row indicate three start codons, UUG, CUG, and the very common AUG. It also indicates the first amino acid residue when interpreted as a start: in this case it is all methionine.

Translation tables

Even when working with ordinary eukaryotic sequences such as the Yeast genome, it is often desired to be able to use alternative translation tables—namely for translation of the mitochondrial genes. Currently the following translation tables are defined by the NCBI Taxonomy Group for the translation of the sequences in GenBank:[29]

See also

References

  1. ^ a b Tirumalai MR, Rivas M, Tran Q, Fox GE (November 2021). "The Peptidyl Transferase Center: a Window to the Past". Microbiol Mol Biol Rev. 85 (4): e0010421. Bibcode:2021MMBR...85...21T. doi:10.1128/MMBR.00104-21. PMC 8579967. PMID 34756086.
  2. ^ Brooker RJ, Widmaier EP, Graham LE, Stiling PD (2014). Biology (Third international student ed.). New York, NY: McGraw Hill Education. p. 249. ISBN 978-981-4581-85-1.
  3. ^ Neill C (1996). Biology (Fourth ed.). The Benjamin/Cummings Publishing Company. pp. 309–310. ISBN 0-8053-1940-9.
  4. ^ Stryer L (2002). Biochemistry (Fifth ed.). W. H. Freeman and Company. p. 826. ISBN 0-7167-4684-0.
  5. ^ Santos, Fenícia Brito; Del-Bem, Luiz-Eduardo (2023). "The Evolution of tRNA Copy Number and Repertoire in Cellular Life". Genes. 14 (1): 27. doi:10.3390/genes14010027. ISSN 2073-4425. PMC 9858662. PMID 36672768.
  6. ^ Moghal A, Mohler K, Ibba M (November 2014). "Mistranslation of the genetic code". FEBS Letters. 588 (23): 4305–10. Bibcode:2014FEBSL.588.4305M. doi:10.1016/j.febslet.2014.08.035. PMC 4254111. PMID 25220850.
  7. ^ Griffiths A (2008). "9". Introduction to Genetic Analysis (9th ed.). New York: W.H. Freeman and Company. pp. 335–339. ISBN 978-0-7167-6887-6.
  8. ^ "Computational Analysis of Genomic Sequences utilizing Machine Learning". scholar.googleusercontent.com. Retrieved 2022-01-12.
  9. ^ Ross JF, Orlowski M (February 1982). "Growth-rate-dependent adjustment of ribosome function in chemostat-grown cells of the fungus Mucor racemosus". Journal of Bacteriology. 149 (2): 650–3. doi:10.1128/JB.149.2.650-653.1982. PMC 216554. PMID 6799491.
  10. ^ Nakamoto T (February 2011). "Mechanisms of the initiation of protein synthesis: in reading frame binding of ribosomes to mRNA". Molecular Biology Reports. 38 (2): 847–55. doi:10.1007/s11033-010-0176-1. PMID 20467902. S2CID 22038744.
  11. ^ Baggett NE, Zhang Y, Gross CA (March 2017). Ibba M (ed.). "Global analysis of translation termination in E. coli". PLOS Genetics. 13 (3): e1006676. doi:10.1371/journal.pgen.1006676. PMC 5373646. PMID 28301469.
  12. ^ Mora L, Zavialov A, Ehrenberg M, Buckingham RH (December 2003). "Stop codon recognition and interactions with peptide release factor RF3 of truncated and chimeric RF1 and RF2 from Escherichia coli". Molecular Microbiology. 50 (5): 1467–76. doi:10.1046/j.1365-2958.2003.03799.x. PMID 14651631.
  13. ^ Schueren F, Thoms S (August 2016). "Functional Translational Readthrough: A Systems Biology Perspective". PLOS Genetics. 12 (8): e1006196. doi:10.1371/JOURNAL.PGEN.1006196. PMC 4973966. PMID 27490485.
  14. ^ Ou X, Cao J, Cheng A, Peppelenbosch MP, Pan Q (March 2019). "Errors in translational decoding: tRNA wobbling or misincorporation?". PLOS Genetics. 15 (3): 2979–2986. doi:10.1371/journal.pgen.1008017. PMC 3158919. PMID 21930591.
  15. ^ Wohlgemuth I, Pohl C, Mittelstaet J, Konevega AL, Rodnina MV (October 2011). "Evolutionary optimization of speed and accuracy of decoding on the ribosome". Philosophical Transactions of the Royal Society of London. Series B, Biological Sciences. 366 (1580): 2979–86. doi:10.1098/rstb.2011.0138. PMC 6438450. PMID 30921315.
  16. ^ Sin C, Chiarugi D, Valleriani A (April 2016). "Quantitative assessment of ribosome drop-off in E. coli". Nucleic Acids Research. 44 (6): 2528–37. doi:10.1093/nar/gkw137. PMC 4824120. PMID 26935582.
  17. ^ Awad S, Valleriani A, Chiarugi D (April 2024). "A data-driven estimation of the ribosome drop-off rate in S. cerevisiae reveals a correlation with the genes length". NAR Genomics and Bioinformatics. 6 (2): lqae036. doi:10.1093/nargab/lqae036. PMC 11025885. PMID 38638702.
  18. ^ Ingolia NT, Ghaemmaghami S, Newman JR, Weissman JS (April 2009). "Genome-wide analysis in vivo of translation with nucleotide resolution using ribosome profiling". Science. 324 (5924): 218–23. Bibcode:2009Sci...324..218I. doi:10.1126/science.1168978. PMC 2746483. PMID 19213877.
  19. ^ Cenik C, Cenik ES, Byeon GW, Grubert F, Candille SI, Spacek D, et al. (November 2015). "Integrative analysis of RNA, translation, and protein levels reveals distinct regulatory variation across humans". Genome Research. 25 (11): 1610–21. doi:10.1101/gr.193342.115. PMC 4617958. PMID 26297486.
  20. ^ Ozadam H, Tonn T, Han CM, Segura A, Hoskins I, Rao S; et al. (2023). "Single-cell quantification of ribosome occupancy in early mouse development". Nature. 618 (7967): 1057–1064. Bibcode:2023Natur.618.1057O. doi:10.1038/s41586-023-06228-9. PMC 10307641. PMID 37344592.{{cite journal}}: CS1 maint: multiple names: authors list (link)
  21. ^ a b Xu Y, Ruggero D (March 2020). "The Role of Translation Control in Tumorigenesis and Its Therapeutic Implications". Annual Review of Cancer Biology. 4 (1): 437–457. doi:10.1146/annurev-cancerbio-030419-033420.
  22. ^ Truitt ML, Ruggero D (April 2016). "New frontiers in translational control of the cancer genome". Nature Reviews. Cancer. 16 (5): 288–304. doi:10.1038/nrc.2016.27. PMC 5491099. PMID 27112207.
  23. ^ a b c Gorban AN, Harel-Bellan A, Morozova N, Zinovyev A (July 2019). "Basic, simple and extendable kinetic model of protein synthesis". Mathematical Biosciences and Engineering. 16 (6): 6602–6622. arXiv:1204.5941. doi:10.3934/mbe.2019329. PMID 31698578.
  24. ^ MacDonald CT, Gibbs JH, Pipkin AC (1968). "Kinetics of biopolymerization on nucleic acid templates". Biopolymers. 6 (1): 1–5. doi:10.1002/bip.1968.360060102. PMID 5641411. S2CID 27559249.
  25. ^ Heinrich R, Rapoport TA (September 1980). "Mathematical modelling of translation of mRNA in eucaryotes; steady state, time-dependent processes and application to reticulocytes". Journal of Theoretical Biology. 86 (2): 279–313. Bibcode:1980JThBi..86..279H. doi:10.1016/0022-5193(80)90008-9. PMID 7442295.
  26. ^ a b Skjøndal-Bar N, Morris DR (January 2007). "Dynamic model of the process of protein synthesis in eukaryotic cells". Bulletin of Mathematical Biology. 69 (1): 361–93. doi:10.1007/s11538-006-9128-2. PMID 17031456. S2CID 83701439.
  27. ^ Coyte KZ, Tabuteau H, Gaffney EA, Foster KR, Durham WM (April 2017). "Reply to Baveye and Darnault: Useful models are simple and extendable". Proceedings of the National Academy of Sciences of the United States of America. 114 (14): E2804–E2805. Bibcode:2017PNAS..114E2804C. doi:10.1073/pnas.1702303114. PMC 5389313. PMID 28341710.
  28. ^ Morozova N, Zinovyev A, Nonne N, Pritchard LL, Gorban AN, Harel-Bellan A (September 2012). "Kinetic signatures of microRNA modes of action". RNA. 18 (9): 1635–55. doi:10.1261/rna.032284.112. PMC 3425779. PMID 22850425.
  29. ^ a b Elzanowski, Andrzej; Ostell, Jim (January 2019). "The Genetic Codes". National Center for Biotechnology Information (NCBI). Retrieved 31 May 2022.

Further reading