A comprehensive catalogue of somatic mutations from a human cancer genome.

All cancers carry somatic mutations. A subset of these somatic alterations, termed driver mutations, confer selective growth advantage and are implicated in cancer development, whereas the remainder are passengers. Here we have sequenced the genomes of a malignant melanoma and a lymphoblastoid cell line from the same person, providing the first comprehensive catalogue of somatic mutations from an individual cancer. The catalogue provides remarkable insights into the forces that have shaped this cancer genome. The dominant mutational signature reflects DNA damage due to ultraviolet light exposure, a known risk factor for malignant melanoma, whereas the uneven distribution of mutations across the genome, with a lower prevalence in gene footprints, indicates that DNA repair has been preferentially deployed towards transcribed regions. The results illustrate the power of a cancer genome sequence to reveal traces of the DNA damage, repair, mutation and selection processes that were operative years before the cancer became symptomatic.

As old as the mountains: the radiations of the Ericaceae.

Mountains are often more species-rich than lowlands. This could be the result of migration from lowlands to mountains, of a greater survival rate in mountains, or of a higher diversification rate in mountains. We investigated this question in the globally distributed family Ericaceae, which includes c. 4426 species ranging from sea level to > 5000 m. We predict that the interaction of low specific leaf area (SLA) and montane habitats is correlated with increased diversification rates. A molecular phylogeny of Ericaceae based on rbcL and matK sequence data was built and dated with 18 fossil calibrations and divergence time estimates. We identified radiations using bamm and correlates of diversification rate changes using binary-state speciation and extinction (BiSSE) and multiple-state speciation and extinction (MuSSE) analyses. Analyses revealed six largely montane radiations. Lineages in mountains diversified faster than nonmountain lineages (higher speciation rate, but no difference in extinction rate), and lineages with low SLA diversified faster than high-SLA lineages. Further, habitat and trait had a positive interactive effect on diversification. Our results suggest that the species richness in mountains is the result of increased speciation rather than reduced extinction or increased immigration. Increased speciation in Ericaceae was facilitated by low SLA.

Accurate whole human genome sequencing using reversible terminator chemistry.

DNA sequence information underpins genetic research, enabling discoveries of important biological or medical benefit. Sequencing projects have traditionally used long (400-800 base pair) reads, but the existence of reference sequences for the human and many other genomes makes it possible to develop new, fast approaches to re-sequencing, whereby shorter reads are compared to a reference to identify intraspecies genetic variation. Here we report an approach that generates several billion bases of accurate nucleotide sequence per experiment at low cost. Single molecules of DNA are attached to a flat surface, amplified in situ and used as templates for synthetic sequencing with fluorescent reversible terminator deoxyribonucleotides. Images of the surface are analysed to generate high-quality sequence. We demonstrate application of this approach to human genome sequencing on flow-sorted X chromosomes and then scale the approach to determine the genome sequence of a male Yoruba from Ibadan, Nigeria. We build an accurate consensus sequence from >30x average depth of paired 35-base reads. We characterize four million single-nucleotide polymorphisms and four hundred thousand structural variants, many of which were previously unknown. Our approach is effective for accurate, rapid and economical whole-genome re-sequencing and many other biomedical applications.

Are pollen fossils useful for calibrating relaxed molecular clock dating of phylogenies? A comparative study using Myrtaceae.

The identification and application of reliable fossil calibrations represents a key component of many molecular studies of evolutionary timescales. In studies of plants, most paleontological calibrations are associated with macrofossils. However, the pollen record can also inform age calibrations if fossils matching extant pollen groups are found. Recent work has shown that pollen of the myrtle family, Myrtaceae, can be classified into a number of morphological groups that are synapomorphic with molecular groups. By assembling a data matrix of pollen morphological characters from extant and fossil Myrtaceae, we were able to measure the fit of 26 pollen fossils to a molecular phylogenetic tree using parsimony optimisation of characters. We identified eight Myrtaceidites fossils as appropriate for calibration based on the most parsimonious placements of these fossils on the tree. These fossils were used to inform age constraints in a Bayesian phylogenetic analysis of a sequence alignment comprising two sequences from the chloroplast genome (matK and ndhF) and one nuclear locus (ITS), sampled from 106 taxa representing 80 genera. Three additional analyses were calibrated by placing pollen fossils using geographic and morphological information (eight calibrations), macrofossils (five calibrations), and macrofossils and pollen fossils in combination (12 calibrations). The addition of new fossil pollen calibrations led to older crown ages than have previously been found for tribes such as Eucalypteae and Myrteae. Estimates of rate variation among lineages were affected by the choice of calibrations, suggesting that the use of multiple calibrations can improve estimates of rate heterogeneity among lineages. This study illustrates the potential of including pollen-based calibrations in molecular studies of divergence times.

Emergence and dynamics of self-producing information niches as a step towards pre-evolutionary organization.

As a step towards understanding pre-evolutionary organization in non-genetic systems, we develop a model to investigate the emergence and dynamics of proto-autopoietic networks in an interacting population of simple information processing entities (automata). Our simulations indicate that dynamically stable strongly connected networks of mutually producing communication channels emerge under specific environmental conditions. We refer to these distinct organizational steady states as information niches In each case, we measure the information content by the Shannon entropy, and determine the fitness landscape, robustness and transition pathways for information niches subjected to intermittent environmental perturbations under non-evolutionary conditions. By determining the information required to generate each niche, we show that niche transitions are only allowed if accompanied by an equal or increased level of information production that arises internally or via environmental perturbations that serve as an exogenous source of population diversification. Overall, our simulations show how proto-autopoietic networks of basic information processors form and compete, and under what conditions they persist over time or go extinct. These findings may be relevant to understanding how inanimate systems such as chemically communicating protocells can initiate the transition to living matter prior to the onset of contemporary evolutionary and genetic mechanisms.

Do Mediterranean-type ecosystems have a common history?--insights from the Buckthorn family (Rhamnaceae).

Mediterranean-type ecosystems (MTEs) are remarkable in their species richness and endemism, but the processes that have led to this diversity remain enigmatic. Here, we hypothesize that continent-dependent speciation and extinction rates have led to disparity in diversity between the five MTEs of the world: the Cape, California, Mediterranean Basin, Chile, and Western Australia. To test this hypothesis, we built a phylogenetic tree for 280 Rhamnaceae species, estimated divergence times using eight fossil calibrations, and used Bayesian methods and simulations to test for differences in diversification rates. Rhamnaceae lineages in MTEs generally show higher diversification rates than elsewhere, but speciation and extinction dynamics show a pattern of continent-dependence. We detected high speciation and extinction rates in California and significantly lower extinction rates in the Cape and Western Australia. The independent colonization of four of five MTEs may have occurred conterminously in the Oligocene/Early Miocene, but colonization of the Mediterranean Basin happened later, in the Late Miocene. This suggests that the in situ radiations of these clades were initiated before the onset of winter rainfall in these regions. These results indicate independent evolutionary histories of Rhamnaceae in MTEs, possibly related to the intensity of climate oscillations and the geological history of the regions.

Fossils and a large molecular phylogeny show that the evolution of species richness, generic diversity, and turnover rates are disconnected.

The magnitude and extent of global change during the Cenozoic is remarkable, yet the impacts of these global changes on the biodiversity and evolutionary dynamics of species diversification remain poorly understood. To investigate this question, we combine paleontological and neontological data for the angiosperm order Fagales, an ecologically important clade of about 1370 species of trees with an exceptional fossil record. We show differences in patterns of accumulation of generic diversity, species richness, and turnover rates for Fagales. Generic diversity evolved rapidly since the Late Cretaceous and peaked during the Eocene or Oligocene. Turnover rates were high during periods of extreme global climate change, but relatively low when the climate remained stable. Species richness accumulated gradually throughout the Cenozoic, possibly at an accelerated pace after the Middle Miocene. Species diversification occurred in new environments: Quercoids radiating in Oligocene subtropical seasonally arid habitats, Casuarinaceae in Australian pyrophytic biomes, and Betula in Late Neogene holarctic habitats. These radiations were counterbalanced by regional extinctions in Late Neogene mesic warm-temperate forests. Thus, the overall diversification at species level is linked to regional radiations of clades with appropriate ecologies exploiting newly available habitats.

A narrow group of monophyletic Tulasnella (Tulasnellaceae) symbiont lineages are associated with multiple species of Chiloglottis (Orchidaceae): Implications for orchid diversity.

PREMISE OF THE STUDY: The Orchidaceae is characterized by exceptional species diversity. Obligate orchid mycorrhizae are predicted to determine orchid distributions, and highly specific relationships between orchids and fungi may drive orchid diversification. In this study, mycorrhizal diversity was examined in the terrestrial, photosynthetic orchid genus Chiloglottis to test the hypothesis of mycorrhizal-mediated diversification in the genus Chiloglottis. This orchid genus secures pollination by sexual deception, an obligate and highly specific pollination strategy. Here we asked whether the obligate orchid-fungal interactions are also specific. • METHODS: Two sequenced loci, the internal transcribed spacer region (ITS) and mitochondrial large subunit (mtLSU), were used to identify fungal isolates and assess fungal species diversity. Symbiotic germination of two species Chiloglottis aff. jeanesii and C. valida were used to assess germination potential of isolates and confirm mycorrhizal association. • KEY RESULTS: Phylogenetic analyses revealed that six representative Chiloglottis species spanning a broad survey of the genus were all associated with a narrow group of monophyletic Tulasnella fungal lineages. • CONCLUSIONS: The Chiloglottis-Tulasnella interaction appears to be the first known case of such a narrow symbiont association across a broadly surveyed orchid genus. It appears that the specific pollination system of Chiloglottis, rather than specific orchid-fungal interactions has been the key driving force in the diversification of the genus. These findings also indicate that plant groups with highly specific mycorrhizal partners can have a widespread distribution.

Genomic diversity among drug sensitive and multidrug resistant isolates of Mycobacterium tuberculosis with identical DNA fingerprints.

BACKGROUND: Mycobacterium tuberculosis complex (MTBC), the causative agent of tuberculosis (TB), is characterized by low sequence diversity making this bacterium one of the classical examples of a genetically monomorphic pathogen. Because of this limited DNA sequence variation, routine genotyping of clinical MTBC isolates for epidemiological purposes relies on highly discriminatory DNA fingerprinting methods based on mobile and repetitive genetic elements. According to the standard view, isolates exhibiting the same fingerprinting pattern are considered direct progeny of the same bacterial clone, and most likely reflect ongoing transmission or disease relapse within individual patients. METHODOLOGY/PRINCIPAL FINDINGS: Here we further investigated this assumption and used massively parallel whole-genome sequencing to compare one drug-susceptible (K-1) and one multidrug resistant (MDR) isolate (K-2) of a rapidly spreading M. tuberculosis Beijing genotype clone from a high incidence region (Karakalpakstan, Uzbekistan). Both isolates shared the same IS6110 RFLP pattern and the same allele at 23 out of 24 MIRU-VNTR loci. We generated 23.9 million (K-1) and 33.0 million (K-2) paired 50 bp purity filtered reads corresponding to a mean coverage of 483.5 fold and 656.1 fold respectively. Compared with the laboratory strain H37Rv both Beijing isolates shared 1,209 SNPs. The two Beijing isolates differed by 130 SNPs and one large deletion. The susceptible isolate had 55 specific SNPs, while the MDR variant had 75 specific SNPs, including the five known resistance-conferring mutations. CONCLUSIONS: Our results suggest that M. tuberculosis isolates exhibiting identical DNA fingerprinting patterns can harbour substantial genomic diversity. Because this heterogeneity is not captured by traditional genotyping of MTBC, some aspects of the transmission dynamics of tuberculosis could be missed or misinterpreted. Furthermore, a valid differentiation between disease relapse and exogenous reinfection might be impossible using standard genotyping tools if the overall diversity of circulating clones is limited. These findings have important implications for clinical trials of new anti-tuberculosis drugs.

Novel crystal form of the ColE1 Rom protein.

The RNA I modulator protein (Rom) acts as a co-regulator of ColE1 plasmid copy number by binding to RNA kissing hairpins and stabilizing their interaction. The structure of Rom has been determined in a new crystal form from X-ray diffraction data to 2.5 A resolution. In this structure, a dimer of the 57-amino-acid protein is found in the asymmetric unit. Each subunit consists almost entirely of two antiparallel alpha-helices joined by a short hairpin bend. The dimer contains a non-crystallographic twofold axis and forms a highly regular four-alpha-helical bundle. The structural packing in this novel crystal form is different from previously known Rom structures. The asymmetric unit contains one dimer, giving a crystal volume per protein weight (V(M)) of 1.83 A(3) Da(-1) and a low solvent content of 30%. Strong packing interactions and low solvation are characteristic of the structure. The Rom protein was cocrystallized with the Tar-Tar* kissing hairpin RNA. Although the electron-density maps do not show bound RNA, altered conformations in the side chains of Rom that are known to be involved in RNA binding have been identified. These results provide additional information about Rom protein conformational flexibility and suggest that the presence of a highly charged polymer such as RNA can promote tight packing of an RNA-binding protein, even when the RNA itself is not observed in the crystal.

High-resolution structure of RNase P protein from Thermotoga maritima.

The structure of RNase P protein from the hyperthermophilic bacterium Thermotoga maritima was determined at 1.2-A resolution by using x-ray crystallography. This protein structure is from an ancestral-type RNase P and bears remarkable similarity to the recently determined structures of RNase P proteins from bacteria that have the distinct, Bacillus type of RNase P. These two types of protein span the extent of bacterial RNase P diversity, so the results generalize the structure of the bacterial RNase P protein. The broad phylogenetic conservation of structure and distribution of potential RNA-binding elements in the RNase P proteins indicate that all of these homologous proteins bind to their cognate RNAs primarily by interaction with the phylogenetically conserved core of the RNA. The protein is found to dimerize through an extensive, well-ordered interface. This dimerization may reflect a mechanism of thermal stability of the protein before assembly with the RNA moiety of the holoenzyme.

Crystal structures of unbound and aminooxyacetate-bound Escherichia coli gamma-aminobutyrate aminotransferase.

The X-ray crystal structures of Escherichia coli gamma-aminobutyrate aminotransferase unbound and bound to the inhibitor aminooxyacetate are reported. The enzyme crystallizes from ammonium sulfate solutions in the P3(2)21 space group with a tetramer in the asymmetric unit. Diffraction data were collected to 2.4 A resolution for the unliganded enzyme and 1.9 A resolution for the aminooxyacetate complex. The overall structure of the enzyme is similar to those of other aminotransferase subgroup II enzymes. The ability of gamma-aminobutyrate aminotransferase to act on primary amine substrates (gamma-aminobutyrate) in the first half-reaction and alpha-amino acids in the second is proposed to be enabled by the presence of Glu211, whose side chain carboxylate alternates between interactions with Arg398 in the primary amine half-reaction and an alternative binding site in the alpha-amino acid half-reaction, in which Arg398 binds the substrate alpha-carboxylate. The specificity for a carboxylate group on the substrate side chain is due primarily to the presence of Arg141, but also requires substantial local main chain rearrangements relative to the structurally homologous enzyme dialkylglycine decarboxylase, which is specific for small alkyl side chains. No iron-sulfur cluster is found in the bacterial enzyme as was found in the pig enzyme [Storici, P., De Biase, D., Bossa, F., Bruno, S., Mozzarelli, A., Peneff, C., Silverman, R. B., and Schirmer, T. (2004) J. Biol. Chem. 279, 363-73.]. The binding of aminooxyacetate causes remarkably small changes in the active site structure, and no large domain movements are observed. Active site structure comparisons with pig gamma-aminobutyrate aminotransferase and dialkylglycine decarboxylase are discussed.