Publisher Correction: DNA sequencing at 40: past, present and future.

In this Review, the year of publication of reference 54 should be 2005, not 2015. In Box 2, "1982: GenBank ( https://www.ncbi.nlm.nih.gov/genbank/statistics/ )" should read "1982: Genbank/ENA/DDBJ" and "2007: NCBI Short Read Archive" should read "2007: NCBI and ENA Short Read Archives"; this is because the launches of these American, European and Japanese databases were coordinated. These errors have not been corrected.

DNA sequencing at 40: past, present and future.

This review commemorates the 40th anniversary of DNA sequencing, a period in which we have already witnessed multiple technological revolutions and a growth in scale from a few kilobases to the first human genome, and now to millions of human and a myriad of other genomes. DNA sequencing has been extensively and creatively repurposed, including as a 'counter' for a vast range of molecular phenomena. We predict that in the long view of history, the impact of DNA sequencing will be on a par with that of the microscope.

Effect of sodium phenylbutyrate/taurursodiol on tracheostomy/ventilation-free survival and hospitalisation in amyotrophic lateral sclerosis: long-term results from the CENTAUR trial.

BACKGROUND: Coformulated sodium phenylbutyrate/taurursodiol (PB/TURSO) was shown to prolong survival and slow functional decline in amyotrophic lateral sclerosis (ALS). OBJECTIVE: Determine whether PB/TURSO prolonged tracheostomy/ventilation-free survival and/or reduced first hospitalisation in participants with ALS in the CENTAUR trial. METHODS: Adults with El Escorial Definite ALS ≤18 months from symptom onset were randomised to PB/ TURSO or placebo for 6 months. Those completing randomised treatment could enrol in an open-label extension (OLE) phase and receive PB/TURSO for ≤30 months. Times to the following individual or combined key events were compared in the originally randomised treatment groups over a period spanning trial start through July 2020 (longest postrandomisation follow-up, 35 months): death, tracheostomy, permanent assisted ventilation (PAV) and first hospitalisation. RESULTS: Risk of any key event was 47% lower in those originally randomised to PB/TURSO (n=87) versus placebo (n=48, 71% of whom received delayed-start PB/TURSO in the OLE phase) (HR=0.53; 95% CI 0.35 to 0.81; p=0.003). Risks of death or tracheostomy/PAV (HR=0.51; 95% CI 0.32 to 0.84; p=0.007) and first hospitalisation (HR=0.56; 95% CI 0.34 to 0.95; p=0.03) were also decreased in those originally randomised to PB/TURSO. CONCLUSIONS: Early PB/TURSO prolonged tracheostomy/PAV-free survival and delayed first hospitalisation in ALS. TRIAL REGISTRATION NUMBER: NCT03127514; NCT03488524.

Regularities of context-dependent codon bias in eukaryotic genes.

Nucleotides surrounding a codon influence the choice of this particular codon from among the group of possible synonymous codons. The strongest influence on codon usage arises from the nucleotide immediately following the codon and is known as the N1 context. We studied the relative abundance of codons with N1 contexts in genes from four eukaryotes for which the entire genomes have been sequenced: Homo sapiens, Drosophila melanogaster, Caenorhabditis elegans and Arabidopsis thaliana. For all the studied organisms it was found that 90% of the codons have a statistically significant N1 context-dependent codon bias. The relative abundance of each codon with an N1 context was compared with the relative abundance of the same 4mer oligonucleotide in the whole genome. This comparison showed that in about half of all cases the context-dependent codon bias could not be explained by the sequence composition of the genome. Ranking statistics were applied to compare context-dependent codon biases for codons from different synonymous groups. We found regularities in N1 context-dependent codon bias with respect to the codon nucleotide composition. Codons with the same nucleotides in the second and third positions and the same N1 context have a statistically significant correlation of their relative abundances.

An exploratory randomised double-blind and placebo-controlled phase 2 study of a combination of baclofen, naltrexone and sorbitol (PXT3003) in patients with Charcot-Marie-Tooth disease type 1A.

BACKGROUND: Charcot-Marie-Tooth type 1A disease (CMT1A) is a rare orphan inherited neuropathy caused by an autosomal dominant duplication of a gene encoding for the structural myelin protein PMP22, which induces abnormal Schwann cell differentiation and dysmyelination, eventually leading to axonal suffering then loss and muscle wasting. We favour the idea that diseases can be more efficiently treated when targeting multiple disease-relevant pathways. In CMT1A patients, we therefore tested the potential of PXT3003, a low-dose combination of three already approved compounds (baclofen, naltrexone and sorbitol). Our study conceptually builds on preclinical experiments highlighting a pleiotropic mechanism of action that includes downregulation of PMP22. The primary objective was to assess safety and tolerability of PXT3003. The secondary objective aimed at an exploratory analysis of efficacy of PXT3003 in CMT1A, to be used for designing next clinical development stages (Phase 2b/3). METHODS: 80 adult patients with mild-to-moderate CMT1A received in double-blind for 1 year Placebo or one of the three increasing doses of PXT3003 tested, in four equal groups. Safety and tolerability were assessed with the incidence of related adverse events. Efficacy was assessed using the Charcot-Marie-Tooth Neuropathy Score (CMTNS) and the Overall Neuropathy Limitations Scale (ONLS) as main endpoints, as well as various clinical and electrophysiological outcomes. RESULTS: This trial confirmed the safety and tolerability of PXT3003. The highest dose (HD) showed consistent evidence of improvement beyond stabilization. CMTNS and ONLS, with a significant improvement of respectively of 8% (0.4% - 16.2%) and 12.1% (2% - 23.2%) in the HD group versus the pool of all other groups, appear to be the most sensitive clinical endpoints to treatment despite their quasi-stability over one year under Placebo. Patients who did not deteriorate over one year were significantly more frequent in the HD group. CONCLUSIONS: These results confirm that PXT3003 deserves further investigation in adults and could greatly benefit CMT1A-diagnosed children, usually less affected than adults. TRIAL REGISTRATION: EudraCT Number: 2010-023097-40. ClinicalTrials.gov Identifier: NCT01401257. The Committee for Orphan Medicinal Products issued in February 2014 a positive opinion on the application for orphan designation for PXT3003 (EMA/OD/193/13).

The most conserved genome segments for life detection on Earth and other planets.

On Earth, very simple but powerful methods to detect and classify broad taxa of life by the polymerase chain reaction (PCR) are now standard practice. Using DNA primers corresponding to the 16S ribosomal RNA gene, one can survey a sample from any environment for its microbial inhabitants. Due to massive meteoritic exchange between Earth and Mars (as well as other planets), a reasonable case can be made for life on Mars or other planets to be related to life on Earth. In this case, the supremely sensitive technologies used to study life on Earth, including in extreme environments, can be applied to the search for life on other planets. Though the 16S gene has become the standard for life detection on Earth, no genome comparisons have established that the ribosomal genes are, in fact, the most conserved DNA segments across the kingdoms of life. We present here a computational comparison of full genomes from 13 diverse organisms from the Archaea, Bacteria, and Eucarya to identify genetic sequences conserved across the widest divisions of life. Our results identify the 16S and 23S ribosomal RNA genes as well as other universally conserved nucleotide sequences in genes encoding particular classes of transfer RNAs and within the nucleotide binding domains of ABC transporters as the most conserved DNA sequence segments across phylogeny. This set of sequences defines a core set of DNA regions that have changed the least over billions of years of evolution and provides a means to identify and classify divergent life, including ancestrally related life on other planets.

The evolution of spliceosomal introns: patterns, puzzles and progress.

The origins and importance of spliceosomal introns comprise one of the longest-abiding mysteries of molecular evolution. Considerable debate remains over several aspects of the evolution of spliceosomal introns, including the timing of intron origin and proliferation, the mechanisms by which introns are lost and gained, and the forces that have shaped intron evolution. Recent important progress has been made in each of these areas. Patterns of intron-position correspondence between widely diverged eukaryotic species have provided insights into the origins of the vast differences in intron number between eukaryotic species, and studies of specific cases of intron loss and gain have led to progress in understanding the underlying molecular mechanisms and the forces that control intron evolution.

Complex early genes.

We use the pattern of intron conservation in 684 groups of orthologs from seven fully sequenced eukaryotic genomes to provide maximum likelihood estimates of the number of introns present in the same orthologs in various eukaryotic ancestors. We find: (i) intron density in the plant-animal ancestor was high, perhaps two-thirds that of humans and three times that of Drosophila; and (ii) intron density in the ancestral bilateran was also high, equaling that of humans and four times that of Drosophila. We further find that modern introns are generally very old, with two-thirds of modern bilateran introns dating to the ancestral bilateran and two-fifths of modern plant, animal, and fungus introns dating to the plant-animal ancestor. Intron losses outnumber gains over a large range of eukaryotic lineages. These results show that early eukaryotic gene structures were very complex, and that simplification, not embellishment, has dominated subsequent evolution.

Resolution of a deep animal divergence by the pattern of intron conservation.

The relationship between three biologically important groups, arthropods, nematodes, and deuterostomes, remains unresolved. It is unknown whether arthropods are more closely related to nematodes (consistent with the "ecdysozoa" hypothesis) or to deuterostomes (consistent with "coelomata"). We present a method in which we use the pattern of spliceosomal intron conservation to develop a series of inequalities that characterize each possible relationship. We find that only the ecdysozoa grouping satisfies these predictions, with P < 10(-6). Simulations show that our method, unlike some previous methods, is largely insensitive to rate variation between branches.

Rates of intron loss and gain: implications for early eukaryotic evolution.

We study the intron-exon structures of 684 groups of orthologs from seven diverse eukaryotic genomes and provide maximum likelihood estimates for rates and numbers of intron losses and gains in these same genes for a variety of lineages. Rates of intron loss vary from approximately 2 x 10(-9) to 2 x 10(-10) per year. Rates of gain vary from 6 x 10(-13) to 4 x 10(-12) per possible intron insertion site per year. There is an inverse correspondence between rates of intron loss and gain, leading to a 20-fold variation among lineages in the ratio of the rates of the two processes. The observed rates of intron gain are insufficient to explain the large number of introns estimated to have been present in the plant-animal ancestor, suggesting that introns present in early eukaryotes may have been created by a fundamentally different process than more recently gained introns.

The pattern of intron loss.

We studied intron loss in 684 groups of orthologous genes from seven fully sequenced eukaryotic genomes. We found that introns closer to the 3' ends of genes are preferentially lost, as predicted if introns are lost through gene conversion with a reverse transcriptase product of a spliced mRNA. Adjacent introns tend to be lost in concert, as expected if such events span multiple intron positions. Directly contrary to the expectations of some, introns that do not interrupt codons (phase zero) are more, not less, likely to be lost, an intriguing and previously unappreciated result. Adjacent introns with matching phases are not more likely to be retained, as would be expected if they enjoyed a relative selective advantage. The findings of 3' and phase zero intron loss biases are in direct contradiction to an extremely recent study of fungi intron evolution. All patterns are less pronounced in the lineage leading to Caenorhabditis elegans, suggesting that the process of intron loss may be qualitatively different in nematodes. Our results support a reverse transcriptase-mediated model of intron loss.

The universe of exons revisited.

We study the distribution of exons in eukaryotic genes to determine whether one can detect the reuse of exon sequences and to use the frequency of such reuse to estimate how many ancestral exon sequences there might have been. We use two databases of exons. One contained 56,276 internal exons from putatively unrelated genes (less than 20% sequence identity) and the second contained 8917 internal exons from regions of these genes that are homologous and colinear with prokaryotic genes; these are ancient conserved regions (ACRs). At the 95% significance level we find 3500 exon-sequence matches in the large database and 500 matches in the ACR database. These matches correspond to groups of similar sequences. The size-rank relationship for these groups follows a power law, the size falling off as the inverse square root of the rank. This form of the power law distribution leads us to make an estimate for the size of a possible universe of ancestral exons. Using the data corresponding to the ACR regions, that universe is estimated to be about 15,000-30,000 in size.

Phylogenetically older introns strongly correlate with module boundaries in ancient proteins.

The hypothesis that some (but not all) introns were used to construct ancient genes by exon shuffling of modules at the earliest stages of evolution is supported by the finding of an excess of phase-zero intron positions in the boundary regions of such modules in 276 ancient proteins (defined as common to eukaryotes and prokaryotes). Here we show further that as phase-zero intron positions are shared by distant taxa, and thus are truly phylogenetically ancient, their excess in the boundaries becomes greater, rising to an 80% excess if shared by four out of the five taxa: vertebrates, invertebrates, fungi, plants, and protists.

Large-scale comparison of intron positions in mammalian genes shows intron loss but no gain.

We compared intron-exon structures in 1,560 human-mouse orthologs and 360 mouse-rat orthologs. The origin of differences in intron positions between species was inferred by comparison with an outgroup, Fugu for human-mouse and human for mouse-rat. Among 10,020 intron positions in the human-mouse comparison, we found unequivocal evidence for five independent intron losses in the mouse lineage but no evidence for intron loss in humans or for intron gain in either lineage. Among 1,459 positions in rat-mouse comparisons, we found evidence for one loss in rat but neither loss in mouse nor gain in either lineage. In each case, the intron losses were exact, without change in the surrounding coding sequence, and involved introns that are extremely short, with an average of 200 bp, an order of magnitude shorter than the mammalian average. These results favor a model whereby introns are lost through gene conversion with intronless copies of the gene. In addition, the finding of widespread conservation of intron-exon structure, even over large evolutionary distances, suggests that comparative methods employing information about gene structures should be very successful in correctly predicting exon boundaries in genomic sequences.

Mystery of intron gain.

For nearly 15 years, it has been widely believed that many introns were recently acquired by the genes of multicellular organisms. However, the mechanism of acquisition has yet to be described for a single animal intron. Here, we report a large-scale computational analysis of the human, Drosophila melanogaster, Caenorhabditis elegans, and Arabidopsis thaliana genomes. We divided 147,796 human intron sequences into batches of similar lengths and aligned them with each other. Different types of homologies between introns were found, but none showed evidence of simple intron transposition. Also, 106,902 plant, 39,624 Drosophila, and 6021 C. elegans introns were examined. No single case of homologous introns in nonhomologous genes was detected. Thus, we found no example of transposition of introns in the last 50 million years in humans, in 3 million years in Drosophila and C. elegans, or in 5 million years in Arabidopsis. Either new introns do not arise via transposition of other introns or intron transposition must have occurred so early in evolution that all traces of homology have been lost.

Large-scale comparison of intron positions among animal, plant, and fungal genes.

We purge large databases of animal, plant, and fungal intron-containing genes to a 20% similarity level and then identify the most similar animal-plant, animal-fungal, and plant-fungal protein pairs. We identify the introns in each BLAST 2.0 alignment and score matched intron positions and slid (near-matched, within six nucleotides) intron positions automatically. Overall we find that 10% of the animal introns match plant positions, and a further 7% are "slides." Fifteen percent of fungal introns match animal positions, and 13% match plant positions. Furthermore, the number of alignments with high numbers of matches deviates greatly from the Poisson expectation. The 30 animal-plant alignments with the highest matches (for which 44% of animal introns match plant positions) when aligned with fungal genes are also highly enriched for triple matches: 39% of the fungal introns match both animal and plant positions. This is strong evidence for ancestral introns predating the animal-plant-fungal divergence, and in complete opposition to any expectations based on random insertion. In examining the slid introns, we show that at least half are caused by imperfections in the alignments, and are most likely to be actual matches at common positions. Thus, our final estimates are that approximately equal 14% of animal introns match plant positions, and that approximately equal 17-18% of fungal introns match animal or plant positions, all of these being likely to be ancestral in the eukaryotes.

The signal of ancient introns is obscured by intron density and homolog number.

In ancient genes whose products have known 3-dimensional structures, an excess of phase zero introns (those that lie between the codons) appear in the boundaries of modules, compact regions of the polypeptide chain. These excesses are highly significant and could support the hypothesis that ancient genes were assembled by exon shuffling involving compact modules. (Phase one and two introns, and many phase zero introns, appear to arise later.) However, as more genes, with larger numbers of homologs and intron positions, were examined, the effects became smaller, dropping from a 40% excess to an 8% excess as the number of intron positions increased from 570 to 3,328, even though the statistical significance remained strong. An interpretation of this behavior is that novel inserted positions appearing in homologs washed out the signal from a finite number of ancient positions. Here we show that this is likely to be the case. Analyses of intron positions restricted to those in genes for which relatively few intron positions from homologs are known, or to those in genes with a small number of known homologous gene structures, show a significant correlation of phase zero intron positions with the module structure, which weakens as the density of attributed intron positions or the number of homologs increases. These effects do not appear for phase one and phase two introns. This finding matches the expectation of the mixed model of intron origin, in which a fraction of phase zero introns are left from the assembly of the first genes, while other introns have been added in the course of evolution.

Do introns favor or avoid regions of amino acid conservation?

Are intron positions correlated with regions of high amino acid conservation? For a set of ancient conserved proteins, with intronless prokaryotic but intron-containing eukaryotic homologs, multiple sequence alignments identified residues invariant throughout evolution. Intron positions between codons show no preferences. However, introns lying after the first base of a codon prefer conserved regions, markedly in glycines. Because glycines are in excess in conserved regions, this behavior could reflect phase-one introns entering glycine residues randomly in the ancestral sequences. Examination of intron positions within codons of evolutionarily invariable amino acids showed that roughly 50% of these introns are bordered by guanines at both 5'- and 3'-ends, 25% have a G only before the intron, and 5% have a G only after the intron, whereas about 20% are bordered by nonguanine bases.