Mechanisms of gene duplication and amplification

Suppressors of dGTP Starvation in Escherichia coli

dGTP starvation, a newly discovered phenomenon in which Escherichia coli cells are starved specifically for the DNA precursor dGTP, leads to impaired growth and, ultimately, cell death. Phenomenologically, it represents an example of nutritionally induced unbalanced growth: cell mass amplifies normally as dictated by the nutritional status of the medium, but DNA content growth is specifically impaired. The other known example of such a condition, thymineless death (TLD), involves starvation for the DNA precursor dTTP, which has been found to have important chemotherapeutic applications. Experimentally, dGTP starvation is induced by depriving an E. coligpt optA1 strain of its required purine source, hypoxanthine. In our studies of this phenomenon, we noted the emergence of a relatively high frequency of suppressor mutants that proved resistant to the treatment. To study such suppressors, we used next-generation sequencing on a collection of independently obtained mutants. A significant fraction was found to carry a defect in the PurR transcriptional repressor, controlling de novo purine biosynthesis, or in its downstream purEK operon. Thus, upregulation of de novo purine biosynthesis appears to be a major mode of overcoming the lethal effects of dGTP starvation. In addition, another large fraction of the suppressors contained a large tandem duplication of a 250- to 300-kb genomic region that included the purEK operon as well as the acrAB-encoded multidrug efflux system. Thus, the suppressive effects of the duplications could potentially involve beneficial effects of a number of genes/operons within the amplified regions.IMPORTANCE Concentrations of the four precursors for DNA synthesis (2'-deoxynucleoside-5'-triphosphates [dNTPs]) are critical for both the speed of DNA replication and its accuracy. Previously, we investigated consequences of dGTP starvation, where the DNA precursor dGTP was specifically reduced to a low level. Under this condition, E. coli cells continued cell growth but eventually developed a DNA replication defect, leading to cell death due to formation of unresolvable DNA structures. Nevertheless, dGTP-starved cultures eventually resumed growth due to the appearance of resistant mutants. Here, we used whole-genome DNA sequencing to identify the responsible suppressor mutations. We show that the majority of suppressors can circumvent death by upregulating purine de novo biosynthesis, leading to restoration of dGTP to acceptable levels.

Plasmid and clonal interference during post horizontal gene transfer evolution

Plasmids are nucleic acid molecules that can drive their own replication in a living cell. They can be transmitted horizontally and can thrive in the host cell to high-copy numbers. Plasmid replication and gene expression consume cellular resources and cells carrying plasmids incur fitness costs. But many plasmids carry genes that can be beneficial under certain conditions, allowing the cell to endure in the presence of antibiotics, toxins, competitors or parasites. Horizontal transfer of plasmid-encoded genes can thus instantaneously confer differential adaptation to local or transient selection conditions. This conflict between cellular fitness and plasmid spread sets the scene for multilevel selection processes. We have engineered a system to study the short-term evolutionary impact of different synonymous versions of a plasmid-encoded antibiotic resistance gene. Applying experimental evolution under different selection conditions and deep sequencing allowed us to show rapid local adaptation to the presence of antibiotic and to the specific version of the resistance gene transferred. We describe the presence of clonal interference at two different levels: at the within-cell level, because a single cell can carry several plasmids, and at the between-cell level, because a bacterial population may contain several clones carrying different plasmids and displaying different fitness in the presence/absence of antibiotic. Understanding the within-cell and between-cell dynamics of plasmids after horizontal gene transfer is essential to unravel the dense network of mobile elements underlying the worldwide threat to public health of antibiotic resistance.

Fallacy of the Unique Genome: Sequence Diversity within Single Helicobacter pylori Strains

Many bacterial genomes are highly variable but nonetheless are typically published as a single assembled genome. Experiments tracking bacterial genome evolution have not looked at the variation present at a given point in time. Here, we analyzed the mouse-passaged Helicobacter pylori strain SS1 and its parent PMSS1 to assess intra- and intergenomic variability. Using high sequence coverage depth and experimental validation, we detected extensive genome plasticity within these H. pylori isolates, including movement of the transposable element IS607, large and small inversions, multiple single nucleotide polymorphisms, and variation in cagA copy number. The cagA gene was found as 1 to 4 tandem copies located off the cag island in both SS1 and PMSS1; this copy number variation correlated with protein expression. To gain insight into the changes that occurred during mouse adaptation, we also compared SS1 and PMSS1 and observed 46 differences that were distinct from the within-genome variation. The most substantial was an insertion in cagY, which encodes a protein required for a type IV secretion system function. We detected modifications in genes coding for two proteins known to affect mouse colonization, the HpaA neuraminyllactose-binding protein and the FutB α-1,3 lipopolysaccharide (LPS) fucosyltransferase, as well as genes predicted to modulate diverse properties. In sum, our work suggests that data from consensus genome assemblies from single colonies may be misleading by failing to represent the variability present. Furthermore, we show that high-depth genomic sequencing data of a population can be analyzed to gain insight into the normal variation within bacterial strains.IMPORTANCE Although it is well known that many bacterial genomes are highly variable, it is nonetheless traditional to refer to, analyze, and publish "the genome" of a bacterial strain. Variability is usually reduced ("only sequence from a single colony"), ignored ("just publish the consensus"), or placed in the "too-hard" basket ("analysis of raw read data is more robust"). Now that whole-genome sequences are regularly used to assess virulence and track outbreaks, a better understanding of the baseline genomic variation present within single strains is needed. Here, we describe the variability seen in typical working stocks and colonies of pathogen Helicobacter pylori model strains SS1 and PMSS1 as revealed by use of high-coverage mate pair next-generation sequencing (NGS) and confirmed by traditional laboratory techniques. This work demonstrates that reliance on a consensus assembly as "the genome" of a bacterial strain may be misleading.

Palindromic amplification of the ERBB2 oncogene in primary HER2-positive breast tumors

Oncogene amplification confers a growth advantage to tumor cells for clonal expansion. There are several, recurrently amplified oncogenes throughout the human genome. However, it remains unclear whether this recurrent amplification is solely a manifestation of increased fitness resulting from random amplification mechanisms, or if a genomic locus-specific amplification mechanism plays a role. Here we show that the ERBB2 oncogene at 17q12 is susceptible to palindromic gene amplification, a mechanism characterized by the inverted (palindromic) duplication of genomic segments, in HER2-positive breast tumors. We applied two genomic approaches to investigate amplification mechanisms: sequencing of DNA libraries enriched with tumor-derived palindromic DNA (Genome-wide Analysis of Palindrome Formation) and whole genome sequencing (WGS). We observed significant enrichment of palindromic DNA within amplified ERBB2 genomic segments. Palindromic DNA was particularly enriched at amplification peaks and at boundaries between amplified and normal copy-number regions. Thus, palindromic gene amplification shaped the amplified ERBB2 locus. The enrichment of palindromic DNA throughout the amplified segments leads us to propose that the ERBB2 locus is amplified through the mechanism that repeatedly generates palindromic DNA, such as Breakage-Fusion-Bridge cycles. The genomic architecture surrounding ERBB2 in the normal genome, such as segmental duplications, could promote the locus-specific mechanism.

Chromosomal Replication Complexity: A Novel DNA Metrics and Genome Instability Factor

As the ratio of the copy number of the most replicated to the unreplicated regions in the same chromosome, the definition of chromosomal replication complexity (CRC) appears to leave little room for variation, being either two during S-phase or one otherwise. However, bacteria dividing faster than they replicate their chromosome spike CRC to four and even eight. A recent experimental inquiry about the limits of CRC in Escherichia coli revealed two major reasons to avoid elevating it further: (i) increased chromosomal fragmentation and (ii) complications with subsequent double-strand break repair. Remarkably, examples of stable elevated CRC in eukaryotic chromosomes are well known under various terms like "differential replication," "underreplication," "DNA puffs," "onion-skin replication," or "re-replication" and highlight the phenomenon of static replication fork (sRF). To accurately describe the resulting "amplification by overinitiation," I propose a new term: "replification" (subchromosomal overreplication). In both prokaryotes and eukaryotes, replification, via sRF processing, causes double-strand DNA breaks and, with their repair elevating chromosomal rearrangements, represents a novel genome instability factor. I suggest how static replication bubbles could be stabilized and speculate that some tandem duplications represent such persistent static bubbles. Moreover, I propose how static replication bubbles could be transformed into tandem duplications, double minutes, or inverted triplications. Possible experimental tests of these models are discussed.

Functional Diversification of the Four MARCKS Family Members in Zebrafish Neural Development

Myristoylated alanin-rich C-kinase substrate (MARCKS) and MARCKS-like 1, each encoded by a different gene, comprise a very small family of actin-modulating proteins with essential roles in mammalian neural development. We show here that four genes (two marcks and two marcksl1) are present in teleosts including zebrafish, while ancient actinopterigians, sarcopterigian fishes, and chondrichtyans only have two. No marcks genes were found in agnaths or invertebrates. All four zebrafish genes are expressed during development, and we show here how their early knockdown causes defects in neural development, with some phenotypical differences. Knockdown of marcksa generated embryos with smaller brain and eyes, while marcksb caused different morphogenetic defects, such as larger hindbrain ventricle and folded retina. marcksl1a and marcksl1b morpholinos also caused smaller eyes and brain, although marcksl1a alone generated larger brain ventricles. At 24 hpf, marcksb caused a wider angle of the hindbrain walls, while marcksl1a showed a "T-shaped" neural tube and alterations in neuroepithelium organization. The double knockdown surprisingly produced new features, which included an increased neuroepithelial disorganization and partial neural tube duplications evident at 48 hpf, suggesting defects in convergent extension. This disorganization was also evident in the retina, although retinal ganglion cells were still able to differentiate. marcksl1b morphants presented a unique retinal phenotype characterized by the occurrence of sporadic ectopic neuronal differentiation. Although only marcksl1a morphant had a clear "ciliary phenotype," all presented significantly shorter cilia. Altogether, our data show that all marcks genes have functions in zebrafish neural development, with some differences that suggest the onset of protein diversification.

Identification of Evolutionarily Conserved Md1 Splice Variants That Regulate Innate Immunity through Differential Induction of NF-кB

Although in mammals the TLR4/myeloid differentiation factor (MD)2/CD14 complex is responsible for the recognition of bacterial LPS, and it is known that the RP105/MD1 complex negatively regulates TLR4 signaling, the evolutionary history of LPS recognition remains enigmatic. Thus, zebrafish has orthologs of mammalian TLR4 (Tlr4a and Tlr4b), RP105, and MD1, but MD2 and CD14 seem to be absent from all fish genomes available to date. In addition, and to make the story more intriguing, zebrafish Tlr4a and Tlr4b do not recognize LPS, whereas the zebrafish Rp105/Md1 complex unexpectedly participates in the regulation of innate immunity and viral resistance. In this work, we report the identification of two novel splice variants of Md1, which are expressed at similar levels as full-length Md1 in the main immune-related organs of zebrafish and are highly induced upon viral infection. One of these splice variants, which is also expressed by mouse macrophages, lacks three conserved cysteine residues that have been shown to form disulfide bonds that are crucial for the three-dimensional structure of the MD-2–related lipid recognition domain of Md1. Functional studies in zebrafish demonstrate that this evolutionarily conserved splice variant shows higher antiviral activity than full-length Md1, but reduced proinflammatory activity, due to an impaired ability to activate the master regulator of inflammation, NF-κB. These results uncover a previously unappreciated evolutionarily conserved Md1 splice variant with important functions in the regulation of innate immunity and the antiviral response in zebrafish, and point to the need for additional functional studies in mammals on this little explored molecule.

Structural and functional analysis of the finished genome of the recently isolated toxic Anabaena sp. WA102

Background

Very few closed genomes of the cyanobacteria that commonly produce toxic blooms in lakes and reservoirs are available, limiting our understanding of the properties of these organisms. A new anatoxin-a-producing member of the Nostocaceae, Anabaena sp. WA102, was isolated from a freshwater lake in Washington State, USA, in 2013 and maintained in non-axenic culture.

Results

The Anabaena sp. WA102 5.7 Mbp genome assembly has been closed with long-read, single-molecule sequencing and separately a draft genome assembly has been produced with short-read sequencing technology. The closed and draft genome assemblies are compared, showing a correlation between long repeats in the genome and the many gaps in the short-read assembly. Anabaena sp. WA102 encodes anatoxin-a biosynthetic genes, as does its close relative Anabaena sp. AL93 (also introduced in this study). These strains are distinguished by differences in the genes for light-harvesting phycobilins, with Anabaena sp. AL93 possessing a phycoerythrocyanin operon. Biologically relevant structural variants in the Anabaena sp. WA102 genome were detected only by long-read sequencing: a tandem triplication of the anaBCD promoter region in the anatoxin-a synthase gene cluster (not triplicated in Anabaena sp. AL93) and a 5-kbp deletion variant present in two-thirds of the population. The genome has a large number of mobile elements (160). Strikingly, there was no synteny with the genome of its nearest fully assembled relative, Anabaena sp. 90.

Conclusion

Structural and functional genome analyses indicate that Anabaena sp. WA102 has a flexible genome. Genome closure, which can be readily achieved with long-read sequencing, reveals large scale (e.g., gene order) and local structural features that should be considered in understanding genome evolution and function.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2738-7) contains supplementary material, which is available to authorized users.

Adaptation by Deletogenic Replication Slippage in a Nascent Symbiont

As a consequence of population level constraints in the obligate, host-associated lifestyle, intracellular symbiotic bacteria typically exhibit high rates of molecular sequence evolution and extensive genome degeneration over the course of their host association. While the rationale for genome degeneration is well understood, little is known about the molecular mechanisms driving this change. To understand these mechanisms we compared the genome of Sodalis praecaptivus, a nonhost associated bacterium that is closely related to members of the Sodalis-allied clade of insect endosymbionts, with the very recently derived insect symbiont Candidatus Sodalis pierantonius. The characterization of indel mutations in the genome of Ca. Sodalis pierantonius shows that the replication system in this organism is highly prone to deletions resulting from polymerase slippage events in regions encoding G+C-rich repetitive sequences. This slippage-prone phenotype is mechanistically associated with the loss of certain components of the bacterial DNA recombination machinery at an early stage in symbiotic life and is expected to facilitate rapid adaptation to the novel host environment. This is analogous to the emergence of mutator strains in both natural and laboratory populations of bacteria, which tend to reach high frequencies in clonal populations due to linkage between the mutator allele and the resulting adaptive mutations.

Evolution of mammalian sound localization circuits: A developmental perspective

Localization of sound sources is a central aspect of auditory processing. A unique feature of mammals is the smooth, tonotopically organized extension of the hearing range to high frequencies (HF) above 10kHz, which likely induced positive selection for novel mechanisms of sound localization. How this change in the auditory periphery is accompanied by changes in the central auditory system is unresolved. I will argue that the major VGlut2(+) excitatory projection neurons of sound localization circuits (dorsal cochlear nucleus (DCN), lateral and medial superior olive (LSO and MSO)) represent serial homologs with modifications, thus being paramorphs. This assumption is based on common embryonic origin from an Atoh1(+)/Wnt1(+) cell lineage in the rhombic lip of r5, same cell birth, a fusiform cell morphology, shared genetic components such as Lhx2 and Lhx9 transcription factors, and similar projection patterns. Such a parsimonious evolutionary mechanism likely accelerated the emergence of neurons for sound localization in all three dimensions. Genetic analyses indicate that auditory nuclei in fish, birds, and mammals receive contributions from the same progenitor lineages. Anatomical and physiological differences and the independent evolution of tympanic ears in vertebrate groups, however, argue for convergent evolution of sound localization circuits in tetrapods (amphibians, reptiles, birds, and mammals). These disparate findings are discussed in the context of the genetic architecture of the developing hindbrain, which facilitates convergent evolution. Yet, it will be critical to decipher the gene regulatory networks underlying development of auditory neurons across vertebrates to explore the possibility of homologous neuronal populations.

Reinterpreting Long-Term Evolution Experiments: Is Delayed Adaptation an Example of Historical Contingency or a Consequence of Intermittent Selection?

Van Hofwegen et al. demonstrated that Escherichia coli rapidly evolves the ability to use citrate when long selective periods are provided (D. J. Van Hofwegen, C. J. Hovde, and S. A. Minnich, J Bacteriol 198:1022-1034, 2016, http://dx.doi.org/10.1128/JB.00831-15). This contrasts with the extreme delay (15 years of daily transfers) seen in the long-term evolution experiments of Lenski and coworkers. Their idea of "historical contingency" may require reinterpretation. Rapid evolution seems to involve selection for duplications of the whole cit locus that are too unstable to contribute when selection is provided in short pulses.

Rapid Evolution of Citrate Utilization by Escherichia coli by Direct Selection Requires citT and dctA

The isolation of aerobic citrate-utilizing Escherichia coli (Cit(+)) in long-term evolution experiments (LTEE) has been termed a rare, innovative, presumptive speciation event. We hypothesized that direct selection would rapidly yield the same class of E. coli Cit(+) mutants and follow the same genetic trajectory: potentiation, actualization, and refinement. This hypothesis was tested with wild-type E. coli strain B and with K-12 and three K-12 derivatives: an E. coli ΔrpoS::kan mutant (impaired for stationary-phase survival), an E. coli ΔcitT::kan mutant (deleted for the anaerobic citrate/succinate antiporter), and an E. coli ΔdctA::kan mutant (deleted for the aerobic succinate transporter). E. coli underwent adaptation to aerobic citrate metabolism that was readily and repeatedly achieved using minimal medium supplemented with citrate (M9C), M9C with 0.005% glycerol, or M9C with 0.0025% glucose. Forty-six independent E. coli Cit(+) mutants were isolated from all E. coli derivatives except the E. coli ΔcitT::kan mutant. Potentiation/actualization mutations occurred within as few as 12 generations, and refinement mutations occurred within 100 generations. Citrate utilization was confirmed using Simmons, Christensen, and LeMaster Richards citrate media and quantified by mass spectrometry. E. coli Cit(+) mutants grew in clumps and in long incompletely divided chains, a phenotype that was reversible in rich media. Genomic DNA sequencing of four E. coli Cit(+) mutants revealed the required sequence of mutational events leading to a refined Cit(+) mutant. These events showed amplified citT and dctA loci followed by DNA rearrangements consistent with promoter capture events for citT. These mutations were equivalent to the amplification and promoter capture CitT-activating mutations identified in the LTEE.IMPORTANCE E. coli cannot use citrate aerobically. Long-term evolution experiments (LTEE) performed by Blount et al. (Z. D. Blount, J. E. Barrick, C. J. Davidson, and R. E. Lenski, Nature 489:513-518, 2012, http://dx.doi.org/10.1038/nature11514 ) found a single aerobic, citrate-utilizing E. coli strain after 33,000 generations (15 years). This was interpreted as a speciation event. Here we show why it probably was not a speciation event. Using similar media, 46 independent citrate-utilizing mutants were isolated in as few as 12 to 100 generations. Genomic DNA sequencing revealed an amplification of the citT and dctA loci and DNA rearrangements to capture a promoter to express CitT, aerobically. These are members of the same class of mutations identified by the LTEE. We conclude that the rarity of the LTEE mutant was an artifact of the experimental conditions and not a unique evolutionary event. No new genetic information (novel gene function) evolved.

Different Facets of Copy Number Changes: Permanent, Transient, and Adaptive

Chromosomal copy number changes are frequently associated with harmful consequences and are thought of as an underlying mechanism for the development of diseases. However, changes in copy number are observed during development and occur during normal biological processes. In this review, we highlight the causes and consequences of copy number changes in normal physiologic processes as well as cover their associations with cancer and acquired drug resistance. We discuss the permanent and transient nature of copy number gains and relate these observations to a new mechanism driving transient site-specific copy gains (TSSGs). Finally, we discuss implications of TSSGs in generating intratumoral heterogeneity and tumor evolution and how TSSGs can influence the therapeutic response in cancer.

Origin-Dependent Inverted-Repeat Amplification: Tests of a Model for Inverted DNA Amplification

DNA replication errors are a major driver of evolution—from single nucleotide polymorphisms to large-scale copy number variations (CNVs). Here we test a specific replication-based model to explain the generation of interstitial, inverted triplications. While no genetic information is lost, the novel inversion junctions and increased copy number of the included sequences create the potential for adaptive phenotypes. The model—Origin-Dependent Inverted-Repeat Amplification (ODIRA)—proposes that a replication error at pre-existing short, interrupted, inverted repeats in genomic sequences generates an extrachromosomal, inverted dimeric, autonomously replicating intermediate; subsequent genomic integration of the dimer yields this class of CNV without loss of distal chromosomal sequences. We used a combination of in vitro and in vivo approaches to test the feasibility of the proposed replication error and its downstream consequences on chromosome structure in the yeast Saccharomyces cerevisiae. We show that the proposed replication error—the ligation of leading and lagging nascent strands to create “closed” forks—can occur in vitro at short, interrupted inverted repeats. The removal of molecules with two closed forks results in a hairpin-capped linear duplex that we show replicates in vivo to create an inverted, dimeric plasmid that subsequently integrates into the genome by homologous recombination, creating an inverted triplication. While other models have been proposed to explain inverted triplications and their derivatives, our model can also explain the generation of human, de novo, inverted amplicons that have a 2:1 mixture of sequences from both homologues of a single parent—a feature readily explained by a plasmid intermediate that arises from one homologue and integrates into the other homologue prior to meiosis. Our tests of key features of ODIRA lend support to this mechanism and suggest further avenues of enquiry to unravel the origins of interstitial, inverted CNVs pivotal in human health and evolution.

Chromosomal aberration such as gene amplification is a common event in human diseases and is often selected during adaptation of microorganism to stress. We proposed a replication-based model to explain the formation of a particular type of genomic aberration: internal inverted DNA amplification with retention of the distal end of the chromosome. In this study, using yeast as a model, we test the feasibility of several of these steps for the formation of an inverted amplification: a specific DNA replication anomaly (1) leading to the formation of a palindromic extrachromosomal circular molecule (2) followed by the homologous reintegration of this molecule into the genome (3). A significant feature of this mode of amplification is that the amplified sequences contain one or more replication origins. The instability of the inverted junctions can lead, through homology driven processes, to more complex genomic structures that contain a partial triplication within a duplicated segment, a structure commonly found associated with human disease.

An Overview of the Molecular Mechanisms of Recombinational DNA Repair

Recombinational DNA repair is a universal aspect of DNA metabolism and is essential for genomic integrity. It is a template-directed process that uses a second chromosomal copy (sister, daughter, or homolog) to ensure proper repair of broken chromosomes. The key steps of recombination are conserved from phage through human, and an overview of those steps is provided in this review. The first step is resection by helicases and nucleases to produce single-stranded DNA (ssDNA) that defines the homologous locus. The ssDNA is a scaffold for assembly of the RecA/RAD51 filament, which promotes the homology search. On finding homology, the nucleoprotein filament catalyzes exchange of DNA strands to form a joint molecule. Recombination is controlled by regulating the fate of both RecA/RAD51 filaments and DNA pairing intermediates. Finally, intermediates that mature into Holliday structures are disjoined by either nucleolytic resolution or topological dissolution.

Regulation of recombination and genomic maintenance

Recombination is a central process to stably maintain and transmit a genome through somatic cell divisions and to new generations. Hence, recombination needs to be coordinated with other events occurring on the DNA template, such as DNA replication, transcription, and the specialized chromosomal functions at centromeres and telomeres. Moreover, regulation with respect to the cell-cycle stage is required as much as spatiotemporal coordination within the nuclear volume. These regulatory mechanisms impinge on the DNA substrate through modifications of the chromatin and directly on recombination proteins through a myriad of posttranslational modifications (PTMs) and additional mechanisms. Although recombination is primarily appreciated to maintain genomic stability, the process also contributes to gross chromosomal arrangements and copy-number changes. Hence, the recombination process itself requires quality control to ensure high fidelity and avoid genomic instability. Evidently, recombination and its regulatory processes have significant impact on human disease, specifically cancer and, possibly, neurodegenerative diseases.

Homologous recombination and human health: the roles of BRCA1, BRCA2, and associated proteins

Homologous recombination (HR) is a major pathway for the repair of DNA double-strand breaks in mammalian cells, the defining step of which is homologous strand exchange directed by the RAD51 protein. The physiological importance of HR is underscored by the observation of genomic instability in HR-deficient cells and, importantly, the association of cancer predisposition and developmental defects with mutations in HR genes. The tumor suppressors BRCA1 and BRCA2, key players at different stages of HR, are frequently mutated in familial breast and ovarian cancers. Other HR proteins, including PALB2 and RAD51 paralogs, have also been identified as tumor suppressors. This review summarizes recent findings on BRCA1, BRCA2, and associated proteins involved in human disease with an emphasis on their molecular roles and interactions.