Mechanisms of gene duplication and amplification

Streptococcus suis manganese transporter mutant as a live attenuated vaccine: Safety, efficacy, and virulence reversion mechanisms

Streptococcus suis is the leading cause of mortality in piglets and is responsible for severe economic losses in the global pork industry. Severe invasive diseases caused by S. suis include sepsis, meningitis, arthritis, and endocarditis. S. suis disease prevention is hampered by the lack of safe and efficacious vaccines. In this study, we constructed an S. suis live attenuated vaccine candidate lacking the major streptococcal manganese transporter, a known virulence determinant of this organism. The safety and efficacy of this live vaccine were evaluated in swine. Our clinical study results showed that when administered at a dose of 1010 CFU, the vaccine strain was safe and efficacious. However, a lower dose of 109 CFU failed to generate significant immune protection. To investigate if an adjuvant could enhance the efficacy of the vaccine at a lower dose, we spiked the vaccine with a polymeric adjuvant and evaluated its performance. Surprisingly, four pigs receiving the adjuvanted vaccine died during the vaccination phase. Pathology, microbiology, and genetic analyses suggested that the vaccine strain reverted to virulence in these animals. Functional genetic analysis found that the vaccine strain acquired compensatory mutations that upregulated the expression of a secondary manganese transporter, which in turn restored the virulence of the vaccine strain. Our results provide a new understanding of S. suis host adaptation mechanisms and useful information for the design of future live-attenuated vaccines.

CONNECTING THE DOTS: EPIGENETIC REGULATION OF EXTRACHROMOSOMAL AND INHERITED DNA AMPLIFICATIONS

DNA amplification has intrigued scientists for decades. Since its discovery, significant progress has been made in understanding the mechanisms promoting DNA amplification and their associated function(s). While DNA copy gains were once thought to be regulated purely by stochastic processes, recent findings have revealed the important role of epigenetic modifications in driving these amplifications and their integration into the genome. Furthermore, advances in genomic technology have enabled detailed characterization of these genomic events in terms of size, structure, formation, and regulation. This review highlights how our understanding of DNA amplifications has evolved over time, tracing its trajectory from initial discovery to the contemporary landscape. We describe how recent discoveries have started to uncover how these genomic events occur by controlled biological processes rather than stochastic mechanisms, presenting opportunities for therapeutic modulation.

Ecology, global diversity and evolutionary mechanisms in the Mycobacterium tuberculosis complex

With the COVID-19 pandemic receding, tuberculosis (TB) is again the number one cause of human death to a single infectious agent. TB is caused by bacteria that belong to the Mycobacterium tuberculosis complex (MTBC). Recent advances in genome sequencing have provided new insights into the ecology and evolution of the MTBC. This includes the discovery of new phylogenetic lineages within the MTBC, a deeper understanding of the host tropism among the various animal-adapted lineages, enhanced knowledge on the evolutionary dynamics of antimicrobial resistance and transmission, as well as a better grasp of the within-host MTBC diversity. Moreover, advances in long-read sequencing are increasingly highlighting the relevance of structural genomic variation in the MTBC. These findings not only shed new light on the biology and epidemiology of TB, but also give rise to new questions and research avenues. The purpose of this Review is to summarize these new insights and discuss their implications for global TB control.

To cleave or not and how? The DNA exonucleases and endonucleases in immunity

DNA exonucleases and endonucleases are key executors of the genome during many physiological processes. They generate double-stranded DNA by cleaving damaged endogenous or exogenous DNA, triggering the activation of the innate immune pathways such as cGAS-STING-IFN, and enabling the body to produce anti-viral or anti-tumor immune responses. This is of great significance for maintaining the stability of the genome and improving the therapeutic efficacy of tumors. In addition, genomic instability caused by exonuclease mutations contributes to the development of various autoimmune diseases. This review summarizes the DNA exonucleases and endonucleases which have critical functions in immunity and associated diseases.

Rational Design Assisted by Evolutionary Engineering Allows (De)Construction and Optimization of Complex Phenotypes in Pseudomonas putida KT2440

Beyond the rational construction of genetic determinants to encode target functions, complex phenotype engineering requires the contextualisation of their expression within the metabolic and genetic background of the host strain. Furthermore, wherever metabolic complexity is involved, phenotype engineering demands standard, reliable, plug-and-play tools. We introduce GENIO (GENome Integration and fitness Optimization platform for Pseudomonas putida), a framework to optimise genetic circuit performance by means of (i) chromosome-location-based differential gene expression and (ii) subsequent fitness improvement through evolutionary engineering if needed. Using gene expression strength and cell-to-cell variation, we characterised 10 P. putida chromosomal loci (ppLPS) to show that genome context rather than distance to ORI is the main factor driving differential expression performance. We further contextualised ppLPS gene expression against well-known chromosomal integration sites and plasmids displaying different copy numbers. GENIO supports comprehensive exploration of the gene expression space across P. putida's genome while unlocking performance optimization of complex heterologous metabolic pathways through evolutionary engineering. To demonstrate the usability of GENIO, we restored P. putida's aromatic hydrocarbon metabolism by (de)constructing the toluene/m-xylene catabolic pathway coded in the pWW0 plasmid. We also showed that engineering complex phenotypes requires accurate contextualisation of the synthetic pathways involved, a process that benefits from biological robustness.

Machine learning detection of heteroresistance in Escherichia coli

BACKGROUND

Heteroresistance (HR) is a significant type of antibiotic resistance observed for several bacterial species and antibiotic classes where a susceptible main population contains small subpopulations of resistant cells. Mathematical models, animal experiments and clinical studies associate HR with treatment failure. Currently used susceptibility tests do not detect heteroresistance reliably, which can result in misclassification of heteroresistant isolates as susceptible which might lead to treatment failure. Here we examined if whole genome sequence (WGS) data and machine learning (ML) can be used to detect bacterial HR.

METHODS

We classified 467 Escherichia coli clinical isolates as HR or non-HR to the often used β-lactam/inhibitor combination piperacillin-tazobactam using pre-screening and Population Analysis Profiling tests. We sequenced the isolates, assembled the whole genomes and created a set of predictors based on current knowledge of HR mechanisms. Then we trained several machine learning models on 80% of this data set aiming to detect HR isolates. We compared performance of the best ML models on the remaining 20% of the data set with a baseline model based solely on the presence of β-lactamase genes. Furthermore, we sequenced the resistant sub-populations in order to analyse the genetic mechanisms underlying HR.

FINDINGS

The best ML model achieved 100% sensitivity and 84.6% specificity, outperforming the baseline model. The strongest predictors of HR were the total number of β-lactamase genes, β-lactamase gene variants and presence of IS elements flanking them. Genetic analysis of HR strains confirmed that HR is caused by an increased copy number of resistance genes via gene amplification or plasmid copy number increase. This aligns with the ML model's findings, reinforcing the hypothesis that this mechanism underlies HR in Gram-negative bacteria.

INTERPRETATION

We demonstrate that a combination of WGS and ML can identify HR in bacteria with perfect sensitivity and high specificity. This improved detection would allow for better-informed treatment decisions and potentially reduce the occurrence of treatment failures associated with HR.

FUNDING

Funding provided to DIA from the Swedish Research Council (2021-02091) and NIH (1U19AI158080-01).

Heteroresistance in Enterobacter cloacae complex caused by variation in transient gene amplification events

Heteroresistance (HR) in bacteria describes a subpopulational phenomenon of antibiotic resistant cells of a generally susceptible population. Here, we investigated the molecular mechanisms and phenotypic characteristics underlying HR to ceftazidime (CAZ) in a clinical Enterobacter cloacae complex strain (ECC). We identified a plasmid-borne gene duplication-amplification (GDA) event of a region harbouring an ampC gene encoding a β-lactamase blaDHA-1 as the key determinant of HR. Individual colonies exhibited variations in the copy number of the genes resulting in resistance level variation which correlated with growth onset (lag times) and growth rates in the presence of CAZ. GDA copy number heterogeneity occurred within single resistant colonies, demonstrating heterogeneity of GDA on the single-cell level. The interdependence between GDA, lag time and antibiotic treatment and the strong plasticity underlying HR underlines the high risk for misdetection of antimicrobial HR and subsequent treatment failure.

In silico identification and analysis of paralogs encoding enzymes of carbohydrate metabolism in Drosophila melanogaster

The identification and characterization of gene paralogs is crucial to understand the functional contribution of individual genes/proteins to biological pathways. Here, we have identified 51 genes belonging to fifteen paralogous groups encoding enzymes involved in carbohydrate metabolism in Drosophila melanogaster . Strikingly, most paralogous groups comprise a single 'canonical' enzyme that is expressed ubiquitously and one or more variants expressed predominantly in the testis. Most of these testis-specific forms are predicted to be catalytically inactive, suggesting they may have adopted regulatory roles. This work will aid the planning and interpretation of experimental studies of several Drosophila metabolic pathways, including glycolysis, gluconeogenesis and the pentose phosphate pathway.

A more significant role for insertion sequences in large-scale rearrangements in bacterial genomes

Insertion sequences (ISs) are mobile pieces of DNA that are widespread in bacterial genomes. IS movements typically involve (i) excision of the IS element, (ii) cutting of target site DNA, and (iii) IS element insertion. This process generates a new copy of the IS element and a short duplication at the target site. It has been noted that, for some extant IS copies, no target site duplications (TSDs) are readily identifiable. TSD absence has been attributed to degeneration of the TSD after the insertion event, recombination between identical ISs, or adjacent deletions. Indeed, the latter two-recombination between ISs and adjacent deletions-are frequent causes for the absence of TSDs, which we demonstrate here in an analysis of genome sequence data from the Lenski long-term evolution experiment. Furthermore, we propose that some IS movements-namely, those that occur in association with large-scale genomic rearrangements-do not generate TSDs, and occur without evidence for recombination between ISs or adjacent deletions. In support of this hypothesis, we provide two direct, empirical observations of such IS transposition events: an IS5 movement plus a large deletion in Escherichia coli C, and an IS481 movement occurring with a large duplication in Pseudomonas fluorescens SBW25. Although unlikely, it is possible that the observed deletion and associated IS movement occurred in two successive events in one overnight culture. However, an IS at the center of a large-scale duplication is not readily explained, suggesting that IS element activity may promote both large-scale deletions and duplications.

IMPORTANCE

Insertion sequences are the most common mobile genetic elements found in bacterial genomes, and hence they significantly impact bacterial evolution. We observe insertion sequence movement at the center of large-scale deletions and duplications that occurred during laboratory evolution experiments with Escherichia coli and Pseudomonas fluorescens, involving three distinct types of transposase. We raise the possibility that the transposase does not mediate DNA cleavage but instead inserts into existing DNA breaks. Our research highlights the importance of insertion sequences for the generation of large-scale genomic rearrangements and raises questions concerning the mechanistic basis of these mutations.

Genetically engineered gas vesicle proteins with proliferative potential for synergistic targeted tumor therapy

Nanomedicine enables precision-targeted therapies through a non-invasive approach, and nanoparticles may be biologically affected during their colonization in vivo. Ensuring the efficient expression of their ex vivo performance in vivo, while ensuring biosafety, is of great significance. Previous studies have employed genetically engineered E. coli following in vivo entry as a genetically engineered targeting synergist, to enhance the effect of focused ultrasound ablation by exploiting its targeted colonization of tumor tissue. However, the proliferation process of the actual potentiating nanomaterials, i.e., the aerosol proteins produced by genetically engineered E. coli, in vivo has not been precisely observed. The authors of this paper demonstrate this spatiotemporal change in the expression of gas vesicle proteins while genetically engineered E. coli reproduces following tumor colonization. Based on their targeting and proliferative properties, the authors chose to intervene in the treatment at the maximal gas vesicle protein count to enhance the monitoring and utilization of the potentiator. By examining the therapeutic potential of the novel combination of genetic engineering and focused ultrasound, we present a robust strategy that improves the efficiency of non-invasive treatments.

Evolutionary dynamics and functional characterization of proximal duplicated sorbitol-6-phosphate dehydrogenase genes in Rosaceae

Sorbitol is a critical photosynthate and storage substance in the Rosaceae family. Sorbitol 6-phosphate dehydrogenase (S6PDH) functions as the pivotal rate-limiting enzyme in sorbitol synthesis. The origin and functional diversification of S6PDH in Rosaceae remain unclear, largely due to the complicated interplay of gene duplications. Here, we investigated the synteny relationships among all identified S6PDH genes in representative genomes within the Rosaceae family. By integrating phylogenetic analyses, we elucidated the lineage-specific expansion and syntenic conservation of S6PDH across diverse Rosaceae plant lineages. We found that S6PDH can be traced back to a pair of proximal duplicated genes of the common ancestor of the Rosaceae, and the further amplification of S6PDH in the Maleae primarily relies on WGD events in their lineages. In Rosaceae species, multiple copies of the S6PDH gene are preliminarily divided into two main clades (Clade 1 and Clade 2) based on sequence similarity. These clades have evolved to acquire different functional directions. In Clade 1, lineage-specific transposition events in the Amygdaloideae have led to changes in gene expression patterns and promoted lineage evolution. This is mainly characterized by a decrease in enzymatic activity and transcriptional expression in the leaves, but also includes specific functional diversification, such as sustained post-harvest fruit expression and enhanced expression under biotic stress in certain tissues. In contrast, S6PDH in the Rosoideae and Dryadoideae has not undergone additional duplications beyond early proximal duplication. The loss of exons and variations in exon length might the key factor leading to reduced enzymatic activity in the Clade 2 proximal gene pairs. Collectively, our findings illuminate the dynamic nature of S6PDH evolution and reveal the intricate interplay between duplication, transposition, and functional diversification. This work not only contributes valuable insights into the genetic mechanisms underlying sorbitol metabolism but also establish a crucial foundation for future investigations aimed at comprehensively characterizing the variations of sorbitol metabolism across different subfamilies within the Rosaceae family.

Antibiotic resistance mediated by gene amplifications

Apart from horizontal gene transfer and sequence-altering mutational events, antibiotic resistance can emerge due to the formation of tandem repeats of genomic regions. This phenomenon, also known as gene amplification, has been implicated in antibiotic resistance in both laboratory and clinical scenarios, where the evolution of resistance via amplifications can affect treatment efficacy. Antibiotic resistance mediated by gene amplifications is unstable and consequently can be difficult to detect, due to amplification loss in the absence of the selective pressure of the antibiotic. Further, due to variable copy numbers in a population, amplifications result in heteroresistance, where only a subpopulation is resistant to an antibiotic. While gene amplifications typically lead to resistance by increasing the expression of resistance determinants due to the higher copy number, the underlying mechanisms of resistance are diverse. In this review article, we describe the various pathways by which gene amplifications cause antibiotic resistance, from efflux and modification of the antibiotic, to target modification and bypass. We also discuss how gene amplifications can engender resistance by alternate mutational outcomes such as altered regulation and protein structure, in addition to just an increase in copy number and expression. Understanding how amplifications contribute to bacterial survival following antibiotic exposure is critical to counter their role in the rise of antimicrobial resistance.

Single-Cell Microfluidics: A Primer for Microbiologists

Recent advances in microfluidic technology have made it possible to image live bacterial cells with a high degree of precision and control. In particular, single-cell microfluidic designs have created new opportunities to study phenotypic variation in bacterial populations. However, the development and use of microfluidic devices require specialized resources, and these can be practical barriers to entry for microbiologists. With this review, our intentions are to help demystify the design, construction, and application of microfluidics. Our approach is to present design elements as building blocks from which a multitude of microfluidics applications can be imagined by the microbiologist.

Phage-mediated resolution of genetic conflict alters the evolutionary trajectory of Pseudomonas aeruginosa lysogens

The opportunistic human pathogen Pseudomonas aeruginosa is naturally infected by a large class of temperate, transposable, Mu-like phages. We examined the genotypic and phenotypic diversity of P. aeruginosa PA14 lysogen populations as they resolve clustered regularly interspaced short palindromic repeat (CRISPR) autoimmunity, mediated by an imperfect CRISPR match to the Mu-like DMS3 prophage. After 12 days of evolution, we measured a decrease in spontaneous induction in both exponential and stationary phase growth. Co-existing variation in spontaneous induction rates in the exponential phase depended on the way the coexisting strains resolved genetic conflict. Multiple mutational modes to resolve genetic conflict between host and phage resulted in coexistence in evolved populations of single lysogens that maintained CRISPR immunity to other phages and polylysogens that lost immunity completely. This work highlights a new dimension of the role of lysogenic phages in the evolution of their hosts.IMPORTANCEThe chronic opportunistic multi-drug-resistant pathogen Pseudomonas aeruginosa is persistently infected by temperate phages. We assess the contribution of temperate phage infection to the evolution of the clinically relevant strain UCBPP-PA14. We found that a low level of clustered regularly interspaced short palindromic repeat (CRISPR)-mediated self-targeting resulted in polylysogeny evolution and large genome rearrangements in lysogens; we also found extensive diversification in CRISPR spacers and cas genes. These genomic modifications resulted in decreased spontaneous induction in both exponential and stationary phase growth, increasing lysogen fitness. This work shows the importance of considering latent phage infection in characterizing the evolution of bacterial populations.

Facilitating stable gene integration expression and copy number amplification in Bacillus subtilis through a reversible homologous recombination switch

Strengthening the expression level of integrated genes on the genome is crucial for consistently expressing key enzymes in microbial cell factories for efficient bioproduction in synthetic biology. In comparison to plasmid-based multi-copy expression, the utilization of chromosomal multi-copy genes offers increased stability of expression level, diminishes the metabolic burden on host cells, and enhances overall genetic stability. In this study, we developed the "BacAmp", a stabilized gene integration expression and copy number amplification system for high-level expression in Bacillus subtilis, which was achieved by employing a combination of repressor and non-natural amino acids (ncAA)-dependent expression system to create a reversible switch to control the key gene recA for homologous recombination. When the reversible switch is turned on, genome editing and gene amplification can be achieved. Subsequently, the reversible switch was turned off therefore stabilizing the gene copy number. The stabilized gene amplification system marked by green fluorescent protein, achieved a 3-fold increase in gene expression by gene amplification and maintained the average gene copy number at 10 after 110 generations. When we implemented the gene amplification system for the regulation of N-acetylneuraminic acid (NeuAc) synthesis, the copy number of the critical gene increased to an average of 7.7, which yielded a 1.3-fold NeuAc titer. Our research provides a new avenue for gene expression in synthetic biology and can be applied in metabolic engineering in B. subtilis.

Dynamic Expansions and Retinal Expression of Spectrally Distinct Short-Wavelength Opsin Genes in Sea Snakes

The photopigment-encoding visual opsin genes that mediate color perception show great variation in copy number and adaptive function across vertebrates. An open question is how this variation has been shaped by the interaction of lineage-specific structural genomic architecture and ecological selection pressures. We contribute to this issue by investigating the expansion dynamics and expression of the duplicated Short-Wavelength-Sensitive-1 opsin (SWS1) in sea snakes (Elapidae). We generated one new genome, 45 resequencing datasets, 10 retinal transcriptomes, and 81 SWS1 exon sequences for sea snakes, and analyzed these alongside 16 existing genomes for sea snakes and their terrestrial relatives. Our analyses revealed multiple independent transitions in SWS1 copy number in the marine Hydrophis clade, with at least three lineages having multiple intact SWS1 genes: the previously studied Hydrophis cyanocinctus and at least two close relatives of this species; Hydrophis atriceps and Hydrophis fasciatus; and an individual Hydrophis curtus. In each lineage, gene copy divergence at a key spectral tuning site resulted in distinct UV and Violet/Blue-sensitive SWS1 subtypes. Both spectral variants were simultaneously expressed in the retinae of H. cyanocinctus and H. atriceps, providing the first evidence that these SWS1 expansions confer novel phenotypes. Finally, chromosome annotation for nine species revealed shared structural features in proximity to SWS1 regardless of copy number. If these features are associated with SWS1 duplication, expanded opsin complements could be more common in snakes than is currently recognized. Alternatively, selection pressures specific to aquatic environments could favor improved chromatic distinction in just some lineages.

Transcriptional, post-transcriptional, and post-translational regulation of polyunsaturated fatty acid synthase genes in Aurantiochytrium limacinum strain BL10: Responses to nitrogen starvation

Under nitrogen deficient conditions, the Aurantiochytrium limacinum strain BL10 greatly increases the production of docosahexaenoic acid (DHA) and n-6 docosapentaenoic acid. Researchers have yet to elucidate the mechanism by which BL10 promotes the activity of polyunsaturated fatty acid synthase (Pfa), which plays a key role in the synthesis of polyunsaturated fatty acid (PUFA). Analysis in the current study revealed that in nitrogen-depleted environments, BL10 boosts the transcription and synthesis of proteins by facilitating the expression of pfa genes via transcriptional regulation. It was also determined that BL10 adjusts the lengths of the 5'- and 3'-untranslated regions (suggesting post-transcriptional regulation) and modifies the ratio of two Pfa1 isoforms to favor PUFA production via post-translational regulation (ubiquitination). These findings clarify the exceptional DHA production of BL10 and provide additional insights into the regulatory mechanisms of PUFA biosynthesis in Aurantiochytrium.

Amplification editing enables efficient and precise duplication of DNA from short sequence to megabase and chromosomal scale

Duplication is a foundation of molecular evolution and a driver of genomic and complex diseases. Here, we develop a genome editing tool named Amplification Editing (AE) that enables programmable DNA duplication with precision at chromosomal scale. AE can duplicate human genomes ranging from 20 bp to 100 Mb, a size comparable to human chromosomes. AE exhibits activity across various cell types, encompassing diploid, haploid, and primary cells. AE exhibited up to 73.0% efficiency for 1 Mb and 3.4% for 100 Mb duplications, respectively. Whole-genome sequencing and deep sequencing of the junctions of edited sequences confirm the precision of duplication. AE can create chromosomal microduplications within disease-relevant regions in embryonic stem cells, indicating its potential for generating cellular and animal models. AE is a precise and efficient tool for chromosomal engineering and DNA duplication, broadening the landscape of precision genome editing from an individual genetic locus to the chromosomal scale.

IS1-mediated chromosomal amplification of the arn operon leads to polymyxin B resistance in Escherichia coli B strains

Polymyxins [colistin and polymyxin B (PMB)] comprise an important class of natural product lipopeptide antibiotics used to treat multidrug-resistant Gram-negative bacterial infections. These positively charged lipopeptides interact with lipopolysaccharide (LPS) located in the outer membrane and disrupt the permeability barrier, leading to increased uptake and bacterial cell death. Many bacteria counter polymyxins by upregulating genes involved in the biosynthesis and transfer of amine-containing moieties to increase positively charged residues on LPS. Although 4-deoxy-l-aminoarabinose (Ara4N) and phosphoethanolamine (PEtN) are highly conserved LPS modifications in Escherichia coli, different lineages exhibit variable PMB susceptibilities and frequencies of resistance for reasons that are poorly understood. Herein, we describe a mechanism prevalent in E. coli B strains that depends on specific insertion sequence 1 (IS1) elements that flank genes involved in the biosynthesis and transfer of Ara4N to LPS. Spontaneous and transient chromosomal amplifications mediated by IS1 raise the frequency of PMB resistance by 10- to 100-fold in comparison to strains where a single IS1 element located 90 kb away from the end of the arn operon has been deleted. Amplification involving IS1 becomes the dominant resistance mechanism in the absence of PEtN modification. Isolates with amplified arn operons gradually lose their PMB-resistant phenotype with passaging, consistent with classical PMB heteroresistance behavior. Analysis of the whole genome transcriptome profile showed altered expression of genes residing both within and outside of the duplicated chromosomal segment, suggesting complex phenotypes including PMB resistance can result from tandem amplification events.IMPORTANCEPhenotypic variation in susceptibility and the emergence of resistant subpopulations are major challenges to the clinical use of polymyxins. While a large database of genes and alleles that can confer polymyxin resistance has been compiled, this report demonstrates that the chromosomal insertion sequence (IS) content and distribution warrant consideration as well. Amplification of large chromosomal segments containing the arn operon by IS1 increases the Ara4N content of the lipopolysaccharide layer in Escherichia coli B lineages using a mechanism that is orthogonal to transcriptional upregulation through two-component regulatory systems. Altogether, our work highlights the importance of IS elements in modulating gene expression and generating diverse subpopulations that can contribute to phenotypic polymyxin B heteroresistance.

Suppression of the Escherichia coli rnpA49 conditionally lethal phenotype by different compensatory mutations

RNase P is an essential enzyme found across all domains of life that is responsible for the 5'-end maturation of precursor tRNAs. For decades, numerous studies have sought to elucidate the mechanisms and biochemistry governing RNase P function. However, much remains unknown about the regulation of RNase P expression, the turnover and degradation of the enzyme, and the mechanisms underlying the phenotypes and complementation of specific RNase P mutations, especially in the model bacterium, Escherichia coli In E. coli, the temperature-sensitive (ts) rnpA49 mutation in the protein subunit of RNase P has arguably been one of the most well-studied mutations for examining the enzyme's activity in vivo. Here, we report for the first time naturally occurring temperature-resistant suppressor mutations of E. coli strains carrying the rnpA49 allele. We find that rnpA49 strains can partially compensate the ts defect via gene amplifications of either RNase P subunit (rnpA49 or rnpB) or by the acquisition of loss-of-function mutations in Lon protease or RNase R. Our results agree with previous plasmid overexpression and gene deletion complementation studies, and importantly suggest the involvement of Lon protease in the degradation and/or regulatory pathway(s) of the mutant protein subunit of RNase P. This work offers novel insights into the behavior and complementation of the rnpA49 allele in vivo and provides direction for follow-up studies regarding RNase P regulation and turnover in E. coli.

Genetic variation in recalcitrant repetitive regions of the Drosophila melanogaster genome

Many essential functions of organisms are encoded in highly repetitive genomic regions, including histones involved in DNA packaging, centromeres that are core components of chromosome segregation, ribosomal RNA comprising the protein translation machinery, telomeres that ensure chromosome integrity, piRNA clusters encoding host defenses against selfish elements, and virtually the entire Y chromosome. These regions, formed by highly similar tandem arrays, pose significant challenges for experimental and informatic study, impeding sequence-level descriptions essential for understanding genetic variation. Here, we report the assembly and variation analysis of such repetitive regions in Drosophila melanogaster, offering significant improvements to the existing community reference assembly. Our work successfully recovers previously elusive segments, including complete reconstructions of the histone locus and the pericentric heterochromatin of the X chromosome, spanning the Stellate locus to the distal flank of the rDNA cluster. To infer structural changes in these regions where alignments are often not practicable, we introduce landmark anchors based on unique variants that are putatively orthologous. These regions display considerable structural variation between different D. melanogaster strains, exhibiting differences in copy number and organization of homologous repeat units between haplotypes. In the histone cluster, although we observe minimal genetic exchange indicative of crossing over, the variation patterns suggest mechanisms such as unequal sister chromatid exchange. We also examine the prevalence and scale of concerted evolution in the histone and Stellate clusters and discuss the mechanisms underlying these observed patterns.

Genome-wide analysis of NAC transcription factors and exploration of candidate genes regulating selenium metabolism in Broussonetia papyrifera

MAIN CONCLUSION

Genome-wide identification revealed 79 BpNAC genes belonging to 16 subfamilies, and their gene structures and evolutionary relationships were characterized. Expression analysis highlighted their importance in plant selenium stress responses. Paper mulberry (Broussonetia papyrifera), a deciduous arboreal plant of the Moraceae family, is distinguished by its leaves, which are abundant in proteins, polysaccharides, and flavonoids, positioning it as a novel feedstock. NAC transcription factors, exclusive to plant species, are crucial in regulating growth, development, and response to biotic and abiotic stress. However, extensive characterization of the NAC family within paper mulberry is lacking. In this study, 79 BpNAC genes were identified from the paper mulberry genome, with an uneven distribution across 13 chromosomes. A comprehensive, genome-wide analysis of BpNACs was performed, including investigating gene structures, promoter regions, and chromosomal locations. Phylogenetic tree analysis, alongside comparisons with Arabidopsis thaliana NACs, allowed for categorizing these genes into 16 subfamilies in alignment with gene structure and motif conservation. Collinearity analysis suggested a significant homologous relationship between the NAC genes of paper mulberry and those in Morus notabilis, Ficus hispida, Antiaris toxicaria, and Cannabis sativa. Integrating transcriptome data and Se content revealed that 12 BpNAC genes were associated with selenium biosynthesis. Subsequent RT-qPCR analysis corroborated the correlation between BpNAC59, BpNAC62 with sodium selenate, and BpNAC55 with sodium selenite. Subcellular localization experiments revealed the nuclear functions of BpNAC59 and BpNAC62. This study highlights the potential BpNAC transcription factors involved in selenium metabolism, providing a foundation for strategically breeding selenium-fortified paper mulberry.

Three concurrent mechanisms generate gene copy number variation and transient antibiotic heteroresistance

Heteroresistance is a medically relevant phenotype where small antibiotic-resistant subpopulations coexist within predominantly susceptible bacterial populations. Heteroresistance reduces treatment efficacy across diverse bacterial species and antibiotic classes, yet its genetic and physiological mechanisms remain poorly understood. Here, we investigated a multi-resistant Klebsiella pneumoniae isolate and identified three primary drivers of gene dosage-dependent heteroresistance for several antibiotic classes: tandem amplification, increased plasmid copy number, and transposition of resistance genes onto cryptic plasmids. All three mechanisms imposed fitness costs and were genetically unstable, leading to fast reversion to susceptibility in the absence of antibiotics. We used a mouse gut colonization model to show that heteroresistance due to elevated resistance-gene dosage can result in antibiotic treatment failures. Importantly, we observed that the three mechanisms are prevalent among Escherichia coli bloodstream isolates. Our findings underscore the necessity for treatment strategies that address the complex interplay between plasmids, resistance cassettes, and transposons in bacterial populations.

Comparative genomics of sirenians reveals evolution of filaggrin and caspase-14 upon adaptation of the epidermis to aquatic life

The mammalian epidermis has evolved to protect the body in a dry environment. Genes of the epidermal differentiation complex (EDC), such as FLG (filaggrin), are implicated in the barrier function of the epidermis. Here, we investigated the molecular evolution of the EDC in sirenians (manatees and dugong), which have adapted to fully aquatic life, in comparison to the EDC of terrestrial mammals and aquatic mammals of the clade Cetacea (whales and dolphins). We show that the main subtypes of EDC genes are conserved or even duplicated, like late cornified envelope (LCE) genes of the dugong, whereas specific EDC genes have undergone inactivating mutations in sirenians. FLG contains premature stop codons in the dugong, and the ortholog of human CASP14 (caspase-14), which proteolytically processes filaggrin, is pseudogenized in the same species. As FLG and CASP14 have also been lost in whales, these mutations represent convergent evolution of skin barrier genes in different lineages of aquatic mammals. In contrast to the dugong, the manatee has retained functional FLG and CASP14 genes. FLG2 (filaggrin 2) is truncated in both species of sirenians investigated. We conclude that the land-to-water transition of sirenians was associated with modifications of the epidermal barrier at the molecular level.

Ancient Duplication and Lineage-Specific Transposition Determine Evolutionary Trajectory of ERF Subfamily across Angiosperms

AP2/ERF transcription factor family plays an important role in plant development and stress responses. Previous studies have shed light on the evolutionary trajectory of the AP2 and DREB subfamilies. However, knowledge about the evolutionary history of the ERF subfamily in angiosperms still remains limited. In this study, we performed a comprehensive analysis of the ERF subfamily from 107 representative angiosperm species by combining phylogenomic and synteny network approaches. We observed that the expansion of the ERF subfamily was driven not only by whole-genome duplication (WGD) but also by tandem duplication (TD) and transposition duplication events. We also found multiple transposition events in Poaceae, Brassicaceae, Poales, Brassicales, and Commelinids. These events may have had notable impacts on copy number variation and subsequent functional divergence of the ERF subfamily. Moreover, we observed a number of ancient tandem duplications occurred in the ERF subfamily across angiosperms, e.g., in Subgroup IX, IXb originated from ancient tandem duplication events within IXa. These findings together provide novel insights into the evolution of this important transcription factor family.

Insight on ecDNA-mediated tumorigenesis and drug resistance

Extrachromosomal DNAs (ecDNAs) are a pervasive feature found in cancer and contain oncogenes and their corresponding regulatory elements. Their unique structural properties allow a rapid amplification of oncogenes and alter chromatin accessibility, leading to tumorigenesis and malignant development. The uneven segregation of ecDNA during cell division enhances intercellular genetic heterogeneity, which contributes to tumor evolution that might trigger drug resistance and chemotherapy tolerance. In addition, ecDNA has the ability to integrate into or detach from chromosomal DNA, such progress results into structural alterations and genomic rearrangements within cancer cells. Recent advances in multi-omics analysis revealing the genomic and epigenetic characteristics of ecDNA are anticipated to make valuable contributions to the development of precision cancer therapy. Herein, we conclud the mechanisms of ecDNA generation and the homeostasis of its dynamic structure. In addition to the latest techniques in ecDNA research including multi-omics analysis and biochemical validation methods, we also discuss the role of ecDNA in tumor development and treatment, especially in drug resistance, and future challenges of ecDNA in cancer therapy.

Plasmid-mediated phenotypic noise leads to transient antibiotic resistance in bacteria

The rise of antibiotic resistance is a critical public health concern, requiring an understanding of mechanisms that enable bacteria to tolerate antimicrobial agents. Bacteria use diverse strategies, including the amplification of drug-resistance genes. In this paper, we showed that multicopy plasmids, often carrying antibiotic resistance genes in clinical bacteria, can rapidly amplify genes, leading to plasmid-mediated phenotypic noise and transient antibiotic resistance. By combining stochastic simulations of a computational model with high-throughput single-cell measurements of blaTEM-1 expression in Escherichia coli MG1655, we showed that plasmid copy number variability stably maintains populations composed of cells with both low and high plasmid copy numbers. This diversity in plasmid copy number enhances the probability of bacterial survival in the presence of antibiotics, while also rapidly reducing the burden of carrying multiple plasmids in drug-free environments. Our results further support the tenet that multicopy plasmids not only act as vehicles for the horizontal transfer of genetic information between cells but also as drivers of bacterial adaptation, enabling rapid modulation of gene copy numbers. Understanding the role of multicopy plasmids in antibiotic resistance is critical, and our study provides insights into how bacteria can transiently survive lethal concentrations of antibiotics.

Large-scale duplication events underpin population-level flexibility in tRNA gene copy number in Pseudomonas fluorescens SBW25

The complement of tRNA genes within a genome is typically considered to be a (relatively) stable characteristic of an organism. Here, we demonstrate that bacterial tRNA gene set composition can be more flexible than previously appreciated, particularly regarding tRNA gene copy number. We report the high-rate occurrence of spontaneous, large-scale, tandem duplication events in laboratory populations of the bacterium Pseudomonas fluorescens SBW25. The identified duplications are up to ∼1 Mb in size (∼15% of the wildtype genome) and are predicted to change the copy number of up to 917 genes, including several tRNA genes. The observed duplications are inherently unstable: they occur, and are subsequently lost, at extremely high rates. We propose that this unusually plastic type of mutation provides a mechanism by which tRNA gene set diversity can be rapidly generated, while simultaneously preserving the underlying tRNA gene set in the absence of continued selection. That is, if a tRNA set variant provides no fitness advantage, then high-rate segregation of the duplication ensures the maintenance of the original tRNA gene set. However, if a tRNA gene set variant is beneficial, the underlying duplication fragment(s) may persist for longer and provide raw material for further, more stable, evolutionary change.

Uridine diphosphate glucosyltransferases are involved in spinosad resistance in western flower thrips Frankliniella occidentalis (Pergande)

Uridine diphosphate glucosyltransferases (UGTs) play crucial roles in the insect detoxification system and are associated with pesticide resistance. Our previous transcriptomic analysis of spinosad-susceptible (Ivf03) and resistant (NIL-R) Frankliniella occidentalis revealed numerous upregulated UGT genes in the NIL-R strain, suggesting their potential contribution to spinosad resistance. To investigate this hypothesis, here we conducted UGT activity assays and spinosad induction experiments, employing RNA interference (RNAi) techniques for gene function validation. We found significantly elevated UGT activity in the NIL-R strain compared to Ivf03, with 5-nitrouracil showing a substantial synergistic effect on the resistant strain. Eighteen UGT genes were identified in F. occidentalis, with gene expansion and duplication observed within families UGT466, 467, and 468. Ten out of the eighteen UGTs exhibited higher expression levels in NIL-R, specifically FoUGT466B1, FoUGT468A3, and FoUGT468A4 consistently being upregulated across nymphs, males, and females. RNAi-based functional validation targeting these three UGT genes led to increased susceptibility to spinosad in a life stage-, sex-, and dose-dependent manner. These results indicate that UGTs are indeed involved in spinosad resistance in F. occidentalis, and the effects are dependent on life stage, sex, and dose. Therefore, sustainable control for F. occidentalis resistance should always consider these differential responses.

Duplicated antibiotic resistance genes reveal ongoing selection and horizontal gene transfer in bacteria

Horizontal gene transfer (HGT) and gene duplication are often considered as separate mechanisms driving the evolution of new functions. However, the mobile genetic elements (MGEs) implicated in HGT can copy themselves, so positive selection on MGEs could drive gene duplications. Here, we use a combination of modeling and experimental evolution to examine this hypothesis and use long-read genome sequences of tens of thousands of bacterial isolates to examine its generality in nature. Modeling and experiments show that antibiotic selection can drive the evolution of duplicated antibiotic resistance genes (ARGs) through MGE transposition. A key implication is that duplicated ARGs should be enriched in environments associated with antibiotic use. To test this, we examined the distribution of duplicated ARGs in 18,938 complete bacterial genomes with ecological metadata. Duplicated ARGs are highly enriched in bacteria isolated from humans and livestock. Duplicated ARGs are further enriched in an independent set of 321 antibiotic-resistant clinical isolates. Our findings indicate that duplicated genes often encode functions undergoing positive selection and horizontal gene transfer in microbial communities.

Origin and early divergence of tandem duplicated sorbitol transporter genes in Rosaceae: insights from evolutionary analysis of the SOT gene family in angiosperms

Sorbitol is a critical photosynthate and storage substance in the Rosaceae family. Sorbitol transporters (SOTs) play a vital role in facilitating sorbitol allocation from source to sink organs and sugar accumulation in sink organs. While prior research has addressed gene duplications within the SOT gene family in Rosaceae, the precise origin and evolutionary dynamics of these duplications remain unclear, largely due to the complicated interplay of whole genome duplications and tandem duplications. Here, we investigated the synteny relationships among all identified Polyol/Monosaccharide Transporter (PLT) genes in 61 angiosperm genomes and SOT genes in representative genomes within the Rosaceae family. By integrating phylogenetic analyses, we elucidated the lineage-specific expansion and syntenic conservation of PLTs and SOTs across diverse plant lineages. We found that Rosaceae SOTs, as PLT family members, originated from a pair of tandemly duplicated PLT genes within Class III-A. Furthermore, our investigation highlights the role of lineage-specific and synergistic duplications in Amygdaloideae in contributing to the expansion of SOTs in Rosaceae plants. Collectively, our findings provide insights into the genomic origins, duplication events, and subsequent divergence of SOT gene family members. Such insights lay a crucial foundation for comprehensive functional characterizations in future studies.

Computational inference of eIF4F complex function and structure in human cancers

The canonical eukaryotic initiation factor 4F (eIF4F) complex, composed of eIF4G1, eIF4A1, and the cap-binding protein eIF4E, plays a crucial role in cap-dependent translation initiation in eukaryotic cells. An alternative cap-independent initiation can occur, involving only eIF4G1 and eIF4A1 through internal ribosome entry sites (IRESs). This mechanism is considered complementary to cap-dependent initiation, particularly in tumors under stress conditions. However, the selection and molecular mechanism of specific translation initiation remains poorly understood in human cancers. Thus, we analyzed gene copy number variations (CNVs) in TCGA tumor samples and found frequent amplification of genes involved in translation initiation. Copy number gains in EIF4G1 and EIF3E frequently co-occur across human cancers. Additionally, EIF4G1 expression strongly correlates with genes from cancer cell survival pathways including cell cycle and lipogenesis, in tumors with EIF4G1 amplification or duplication. Furthermore, we revealed that eIF4G1 and eIF4A1 protein levels strongly co-regulate with ribosomal subunits, eIF2, and eIF3 complexes, while eIF4E co-regulates with 4E-BP1, ubiquitination, and ESCRT proteins. Utilizing Alphafold predictions, we modeled the eIF4F structure with and without eIF4E binding. For cap-dependent initiation, our modeling reveals extensive interactions between the N-terminal eIF4E-binding domain of eIF4G1 and eIF4E. Furthermore, the eIF4G1 HEAT-2 domain positions eIF4E near the eIF4A1 N-terminal domain (NTD), resulting in the collaborative enclosure of the RNA binding cavity within eIF4A1. In contrast, during cap-independent initiation, the HEAT-2 domain directly binds the eIF4A1-NTD, leading to a stronger interaction between eIF4G1 and eIF4A1, thus closing the mRNA binding cavity without the involvement of eIF4E.

MET alterations detection platforms and clinical implications in solid tumors: a comprehensive review of literature

MET alterations, including MET exon 14 skipping variants, MET amplification, MET overexpression, and MET fusion, play pivotal roles in primary tumorigenesis and acquired resistance to targeted therapies, especially EGFR tyrosine kinase inhibitors. They represent important diagnostic, prognostic, and predictive biomarkers in many solid tumor types. However, the detection of MET alterations is challenging due to the complexity of MET alterations and the diversity of platform technologies. Therefore, techniques with high sensitivity, specificity, and reliable molecular detection accuracy are needed to overcome such hindrances and aid in biomarker-guided therapies. The current review emphasizes the role of MET alterations as oncogenic drivers in a variety of cancers and their involvement in the development of resistance to targeted therapies. Moreover, our review provides an overview of and recommendations on the selection of various cross-platform technologies for the detection of MET exon 14 skipping variants, MET amplification, MET overexpression, and MET fusion. Furthermore, challenges and hurdles underlying these common detection platforms are discussed.

Receptor tyrosine kinases: biological functions and anticancer targeted therapy

Receptor tyrosine kinases (RTKs) are a class of protein kinases that play crucial roles in various cellular processes, including cell migration, morphological differentiation, cell growth, and angiogenesis. In humans, 58 RTKs have been identified and categorized into 20 distinct families based on the composition of their extracellular regions. RTKs are primarily activated by specific ligands that bind to their extracellular region. They not only regulate tumor transformation, proliferation, metastasis, drug resistance, and angiogenesis, but also initiate and maintain the self-renewal and cloning ability of cancer stem cells. Accurate diagnosis and grading of tumors with dysregulated RTKs are essential in clinical practice. There is a growing body of evidence supporting the benefits of RTKs-targeted therapies for cancer patients, and researchers are actively exploring new targets and developing targeted agents. However, further optimization of RTK inhibitors is necessary to effectively target the diverse RTK alterations observed in human cancers. This review provides insights into the classification, structure, activation mechanisms, and expression of RTKs in tumors. It also highlights the research advances in RTKs targeted anticancer therapy and emphasizes their significance in optimizing cancer diagnosis and treatment strategies.

Prevalent bee venom genes evolved before the aculeate stinger and eusociality

BACKGROUND

Venoms, which have evolved numerous times in animals, are ideal models of convergent trait evolution. However, detailed genomic studies of toxin-encoding genes exist for only a few animal groups. The hyper-diverse hymenopteran insects are the most speciose venomous clade, but investigation of the origin of their venom genes has been largely neglected.

RESULTS

Utilizing a combination of genomic and proteo-transcriptomic data, we investigated the origin of 11 toxin genes in 29 published and 3 new hymenopteran genomes and compiled an up-to-date list of prevalent bee venom proteins. Observed patterns indicate that bee venom genes predominantly originate through single gene co-option with gene duplication contributing to subsequent diversification.

CONCLUSIONS

Most Hymenoptera venom genes are shared by all members of the clade and only melittin and the new venom protein family anthophilin1 appear unique to the bee lineage. Most venom proteins thus predate the mega-radiation of hymenopterans and the evolution of the aculeate stinger.

Negative regulation of MurZ and MurA underlies the essentiality of GpsB- and StkP-mediated protein phosphorylation in Streptococcus pneumoniae D39

GpsB links peptidoglycan synthases to other proteins that determine the shape of the respiratory pathogen Streptococcus pneumoniae (pneumococcus; Spn) and other low-GC Gram-positive bacteria. GpsB is also required for phosphorylation of proteins by the essential StkP(Spn) Ser/Thr protein kinase. Here we report three classes of frequently arising chromosomal duplications (≈21-176 genes) containing murZ (MurZ-family homolog of MurA) or murA that suppress ΔgpsB or ΔstkP. These duplications arose from three different repeated sequences and demonstrate the facility of pneumococcus to modulate gene dosage of numerous genes. Overproduction of MurZ or MurA alone or overproduction of MurZ caused by ΔkhpAB mutations suppressed ΔgpsB or ΔstkP phenotypes to varying extents. ΔgpsB and ΔstkP were also suppressed by MurZ amino-acid changes distant from the active site, including one in commonly studied laboratory strains, and by truncation or deletion of the homolog of IreB(ReoM). Unlike in other Gram-positive bacteria, MurZ is predominant to MurA in pneumococcal cells. However, ΔgpsB and ΔstkP were not suppressed by ΔclpCP, which did not alter MurZ or MurA amounts. These results support a model in which regulation of MurZ and MurA activity, likely by IreB(Spn), is the only essential requirement for StkP-mediated protein phosphorylation in exponentially growing D39 pneumococcal cells.

Computational inference of eIF4F complex function and structure in human cancers

The canonical eukaryotic initiation factor 4F (eIF4F) complex, composed of eIF4G1, eIF4A1, and the cap-binding protein eIF4E, plays a crucial role in cap-dependent translation initiation in eukaryotic cells (1). However, cap-independent initiation can occur through internal ribosomal entry sites (IRESs), involving only eIF4G1 and eIF4A1 present, which is considered to be a complementary process to cap-dependent initiation in tumors under stress conditions (2). The selection and molecular mechanism of specific translation initiation in human cancers remains poorly understood. Thus, we analyzed gene copy number variations (CNVs) in TCGA tumor samples and found frequent amplification of genes involved in translation initiation. Copy number gains in EIF4G1 and EIF3E frequently co-occur across human cancers. Additionally, EIF4G1 expression strongly correlates with genes from cancer cell survival pathways including cell cycle and lipogenesis, in tumors with EIF4G1 amplification or duplication. Furthermore, we revealed that eIF4G1 and eIF4A1 protein levels strongly co-regulate with ribosomal subunits, eIF2, and eIF3 complexes, while eIF4E co-regulates with 4E-BP1, ubiquitination, and ESCRT proteins. Using Alphafold predictions, we modeled the eIF4F structure with and without eIF4G1-eIF4E binding. The modeling for cap-dependent initiation suggests that eIF4G1 interacts with eIF4E through its N-terminal eIF4E-binding domain, bringing eIF4E near the eIF4A1 mRNA binding cavity and closing the cavity with both eIF4G1 HEAT-2 domain and eIF4E. In the cap-independent mechanism, α-helix 5 of eIF4G1 HEAT-2 domain instead directly interacts with the eIF4A1 N-terminal domain to close the mRNA binding cavity without eIF4E involvement, resulting in a stronger interaction between eIF4G1 and eIF4A1.

Significance Statement

Translation initiation is primarily governed by eIF4F, employing a "cap-dependent" mechanism, but eIF4F dysregulation may lead to a "cap-independent" mechanism in stressed cancer cells. We found frequent amplification of translation initiation genes, and co-occurring copy number gains of EIF4G1 and EIF3E genes in human cancers. EIF4G1 amplification or duplication may be positively selected for its beneficial impact on the overexpression of cancer survival genes. The co-regulation of eIF4G1 and eIF4A1, distinctly from eIF4E, reveals eIF4F dysregulation favoring cap-independent initiation. Alphafold predicts changes in the eIF4F complex assembly to accommodate both initiation mechanisms. These findings have significant implications for evaluating cancer cell vulnerability to eIF4F inhibition and developing treatments that target cancer cells with dependency on the translation initiation mechanism.

Gene amplifications cause high-level resistance against albicidin in gram-negative bacteria

Antibiotic resistance is a continuously increasing concern for public healthcare. Understanding resistance mechanisms and their emergence is crucial for the development of new antibiotics and their effective use. The peptide antibiotic albicidin is such a promising candidate that, as a gyrase poison, shows bactericidal activity against a wide range of gram-positive and gram-negative bacteria. Here, we report the discovery of a gene amplification-based mechanism that imparts an up to 1000-fold increase in resistance levels against albicidin. RNA sequencing and proteomics data show that this novel mechanism protects Salmonella Typhimurium and Escherichia coli by increasing the copy number of STM3175 (YgiV), a transcription regulator with a GyrI-like small molecule binding domain that traps albicidin with high affinity. X-ray crystallography and molecular docking reveal a new conserved motif in the binding groove of the GyrI-like domain that can interact with aromatic building blocks of albicidin. Phylogenetic studies suggest that this resistance mechanism is ubiquitous in gram-negative bacteria, and our experiments confirm that STM3175 homologs can confer resistance in pathogens such as Vibrio vulnificus and Pseudomonas aeruginosa.

Oncogenic structural aberration landscape in gastric cancer genomes

Structural variants (SVs) are responsible for driver events in gastric cancer (GC); however, their patterns and processes remain poorly understood. Here, we examine 170 GC whole genomes to unravel the oncogenic structural aberration landscape in GC genomes and identify six rearrangement signatures (RSs). Non-random combinations of RSs elucidate distinctive GC subtypes comprising one or a few dominant RS that are associated with specific driver events (BRCA1/2 defects, mismatch repair deficiency, and TP53 mutation) and epidemiological backgrounds. Twenty-seven SV hotspots are identified as GC driver candidates. SV hotspots frequently constitute complexly clustered SVs involved in driver gene amplification, such as ERBB2, CCNE1, and FGFR2. Further deconstruction of the locally clustered SVs uncovers amplicon-generating profiles characterized by super-large SVs and intensive segmental amplifications, contributing to the extensive amplification of GC oncogenes. Comprehensive analyses using adjusted SV allele frequencies indicate the significant involvement of extra-chromosomal DNA in processes linked to specific RSs.

Lost and Found: The Family of NF-κB Inhibitors Is Larger than Assumed in Salmonid Fish

NF-κB signalling is largely controlled by the family of 'inhibitors of NF-κB' (IκB). The relevant databases indicate that the genome of rainbow trout contains multiple gene copies coding for iκbα (nfkbia), iκbε (nfkbie), iκbδ (nkfbid), iκbζ (nfkbiz), and bcl3, but it lacks iκbβ (nfkbib) and iκbη (ankrd42). Strikingly, three nfkbia paralogs are apparently present in salmonid fish, two of which share a high sequence identity, while the third putative nfkbia gene is significantly less like its two paralogs. This particular nfkbia gene product, iκbα, clusters with the human IκBβ in a phylogenetic analysis, while the other two iκbα proteins from trout associate with their human IκBα counterpart. The transcript concentrations were significantly higher for the structurally more closely related nfkbia paralogs than for the structurally less similar paralog, suggesting that iκbβ probably has not been lost from the salmonid genomes but has been incorrectly designated as iκbα. In the present study, two gene variants coding for iκbα (nfkbia) and iκbε (nfkbie) were prominently expressed in the immune tissues and, particularly, in a cell fraction enriched with granulocytes, monocytes/macrophages, and dendritic cells from the head kidney of rainbow trout. Stimulation of salmonid CHSE-214 cells with zymosan significantly upregulated the iκbα-encoding gene while elevating the copy numbers of the inflammatory markers interleukin-1-beta and interleukin-8. Overexpression of iκbα and iκbε in CHSE-214 cells dose-dependently quenched both the basal and stimulated activity of an NF-κB promoter suggesting their involvement in immune-regulatory processes. This study provides the first functional data on iκbε-versus the well-researched iκbα factor-in a non-mammalian model species.

Mimicking the last step of gene elongation: hints from the bacterial hisF gene

Gene elongation consists in an in-tandem duplication of a gene and divergence and fusion of the two copies, resulting in a gene constituted by two divergent paralogous modules. Many present-day proteins show internal repeats of amino acid sequences, generated by gene elongation events; however, gene elongation is still a poorly studied evolutionary molecular mechanism. The most documented case is that of the histidine biosynthetic genes hisA and hisF, which derive from the gene elongation of an ancestral gene half the size of the extant ones. The aim of this work was to experimentally simulate the possible last step of the gene elongation event occurred during hisF gene evolution under selective pressure conditions. Azospirillum brasilense hisF gene, carrying a single nucleotide mutation that generates a stop codon between the two halves of the gene, was used to transform the histidine-auxotrophic Escherichia coli strain FB182 (hisF892). The transformed strain was subjected to selective pressure (i.e., low concentration/absence of histidine in the growth medium) and the obtained mutants were characterized. The restoration of prototrophy was strongly dependent on the time of incubation and on the strength of the selective pressure. The mutations involved the introduced stop codon with a single base substitution and none of the mutants restored the wild-type codon. Possible correlations between the different mutations and i) E. coli codon usage, ii) three-dimensional structures of the mutated HisF proteins, and iii) growth ability of the mutants were investigated. On the contrary, when the experiment was repeated by mutating a more conserved codon, only a synonymous substitution was obtained. Thus, experiments performed in this study allowed to mimic a possible gene elongation event occurred during the evolution of hisF gene, evidencing the ability of bacterial cells to modify their genome in short times under selective conditions.

Transition from Transient DNA Rereplication to Inherited Gene Amplification Following Prolonged Environmental Stress

Cells require the ability to adapt to changing environmental conditions, however, it is unclear how these changes elicit stable permanent changes in genomes. We demonstrate that, in response to environmental metal exposure, the metallothionein (MT) locus undergoes DNA rereplication generating transient site-specific gene amplifications (TSSGs). Chronic metal exposure allows transition from MT TSSG to inherited MT gene amplification through homologous recombination within and outside of the MT locus. DNA rereplication of the MT locus is suppressed by H3K27me3 and EZH2. Long-term ablation of EZH2 activity eventually leads to integration and inheritance of MT gene amplifications without the selective pressure of metal exposure. The rereplication and inheritance of MT gene amplification is an evolutionarily conserved response to environmental metal from yeast to human. Our results describe a new paradigm for adaptation to environmental stress where targeted, transient DNA rereplication precedes stable inherited gene amplification.

Dosage balance acts as a time-dependent selective barrier to subfunctionalization

BACKGROUND

Gene duplication is an important process for genome expansion, sometimes allowing for new gene functions to develop. Duplicate genes can be retained through multiple processes, either for intermediate periods of time through processes such as dosage balance, or over extended periods of time through processes such as subfunctionalization and neofunctionalization.

RESULTS

Here, we built upon an existing subfunctionalization Markov model by incorporating dosage balance to describe the interplay between subfunctionalization and dosage balance to explore selective pressures on duplicate copies. Our model incorporates dosage balance using a biophysical framework that penalizes the fitness of genetic states with stoichiometrically imbalanced proteins. These imbalanced states cause increased concentrations of exposed hydrophobic surface areas, which cause deleterious mis-interactions. We draw comparison between our Subfunctionalization + Dosage-Balance Model (Sub + Dos) and the previous Subfunctionalization-Only (Sub-Only) Model. This comparison includes how the retention probabilities change over time, dependent upon the effective population size and the selective cost associated with spurious interaction of dosage-imbalanced partners. We show comparison between Sub-Only and Sub + Dos models for both whole-genome duplication and small-scale duplication events.

CONCLUSION

These comparisons show that following whole-genome duplication, dosage balance serves as a time-dependent selective barrier to the subfunctionalization process, by causing an overall delay but ultimately leading to a larger portion of the genome retained through subfunctionalization. This higher percentage of the genome that is ultimately retained is caused by the alternative competing process, nonfunctionalization, being selectively blocked to a greater extent. In small-scale duplication, the reverse pattern is seen, where dosage balance drives faster rates of subfunctionalization, but ultimately leads to a smaller portion of the genome retained as duplicates. This faster rate of subfunctionalization is because the dosage balance of interacting gene products is negatively affected immediately after duplication and the loss of a duplicate restores the stoichiometric balance. Our findings provide support that the subfunctionalization of genes that are susceptible to dosage balance effects, such as proteins involved in complexes, is not a purely neutral process. With stronger selection against stoichiometrically imbalanced gene partners, the rates of subfunctionalization and nonfunctionalization slow; however, this ultimately leads to a greater proportion of subfunctionalized gene pairs.

Dynamics of Gene Loss following Ancient Whole-Genome Duplication in the Cryptic Paramecium Complex

Whole-genome duplications (WGDs) have shaped the gene repertoire of many eukaryotic lineages. The redundancy created by WGDs typically results in a phase of massive gene loss. However, some WGD-derived paralogs are maintained over long evolutionary periods, and the relative contributions of different selective pressures to their maintenance are still debated. Previous studies have revealed a history of three successive WGDs in the lineage of the ciliate Paramecium tetraurelia and two of its sister species from the Paramecium aurelia complex. Here, we report the genome sequence and analysis of 10 additional P. aurelia species and 1 additional out group, revealing aspects of post-WGD evolution in 13 species sharing a common ancestral WGD. Contrary to the morphological radiation of vertebrates that putatively followed two WGD events, members of the cryptic P. aurelia complex have remained morphologically indistinguishable after hundreds of millions of years. Biases in gene retention compatible with dosage constraints appear to play a major role opposing post-WGD gene loss across all 13 species. In addition, post-WGD gene loss has been slower in Paramecium than in other species having experienced genome duplication, suggesting that the selective pressures against post-WGD gene loss are especially strong in Paramecium. A near complete lack of recent single-gene duplications in Paramecium provides additional evidence for strong selective pressures against gene dosage changes. This exceptional data set of 13 species sharing an ancestral WGD and 2 closely related out group species will be a useful resource for future studies on Paramecium as a major model organism in the evolutionary cell biology.

Chromosomal Duplications of MurZ (MurA2) or MurA (MurA1), Amino Acid Substitutions in MurZ (MurA2), and Absence of KhpAB Obviate the Requirement for Protein Phosphorylation in Streptococcus pneumoniae D39

GpsB links peptidoglycan synthases to other proteins that determine the shape of the respiratory pathogen Streptococcus pneumoniae (pneumococcus; Spn ) and other low-GC Gram-positive bacteria. GpsB is also required for phosphorylation of proteins by the essential StkP( Spn ) Ser/Thr protein kinase. Here we report three classes of frequently arising chromosomal duplications (≈21-176 genes) containing murZ (MurZ-family homolog of MurA) or murA that suppress Δ gpsB or Δ stkP . These duplications arose from three different repeated sequences and demonstrate the facility of pneumococcus to modulate gene dosage of numerous genes. Overproduction of MurZ or MurA alone or overexpression of MurZ caused by Δ khpAB mutations suppressed Δ gpsB or Δ stkP phenotypes to varying extents. Δ gpsB and Δ stkP were also suppressed by MurZ amino-acid changes distant from the active site, including one in commonly studied laboratory strains, and by truncation or deletion of the homolog of IreB(ReoM). Unlike in other Gram-positive bacteria, MurZ is predominant to MurA in pneumococcal cells. However, Δ gpsB and Δ stkP were not suppressed by Δ clpCP , which did not alter MurZ or MurA amounts. These results support a model in which regulation of MurZ and MurA activity, likely by IreB( Spn ), is the only essential requirement for protein phosphorylation in exponentially growing D39 pneumococcal cells.

Gene amplification mutations originate prior to selective stress in Acinetobacter baylyi

The controversial theory of adaptive amplification states gene amplification mutations are induced by selective environments where they are enriched due to the stress caused by growth restriction on unadapted cells. We tested this theory with three independent assays using an Acinetobacter baylyi model system that exclusively selects for cat gene amplification mutants. Our results demonstrate all cat gene amplification mutant colonies arise through a multistep process. While the late steps occur during selection exposure, these mutants derive from low-level amplification mutant cells that form before growth-inhibiting selection is imposed. During selection, these partial mutants undergo multiple secondary steps generating higher amplification over several days to multiple weeks to eventually form visible high-copy amplification colonies. Based on these findings, amplification in this Acinetobacter system can be explained by a natural selection process that does not require a stress response. These findings have fundamental implications to understanding the role of growth-limiting selective environments on cancer development. We suggest duplication mutations encompassing growth factor genes may serve as new genomic biomarkers to facilitate early cancer detection and treatment, before high-copy amplification is attained.

Genome-Wide Mining of the Tandem Duplicated Type III Polyketide Synthases and Their Expression, Structure Analysis of Senna tora

Senna tora is one of the homologous crops used as a medicinal food containing an abundance of anthraquinones. Type III polyketide synthases (PKSs) are key enzymes that catalyze polyketide formation; in particular, the chalcone synthase-like (CHS-L) genes are involved in anthraquinone production. Tandem duplication is a fundamental mechanism for gene family expansion. However, the analysis of the tandem duplicated genes (TDGs) and the identification and characterization of PKSs have not been reported for S. tora. Herein, we identified 3087 TDGs in the S. tora genome; the synonymous substitution rates (Ks) analysis indicated that the TDGs had recently undergone duplication. The Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed that the type III PKSs were the most enriched TDGs involved in the biosynthesis of the secondary metabolite pathways, as evidenced by 14 tandem duplicated CHS-L genes. Subsequently, we identified 30 type III PKSs with complete sequences in the S. tora genome. Based on the phylogenetic analysis, the type III PKSs were classified into three groups. The protein conserved motifs and key active residues showed similar patterns in the same group. The transcriptome analysis showed that the chalcone synthase (CHS) genes were more highly expressed in the leaves than in the seeds in S. tora. The transcriptome and qRT-PCR analysis showed that the CHS-L genes had a higher expression in the seeds than in other tissues, particularly seven tandem duplicated CHS-L2/3/5/6/9/10/13 genes. The key active-site residues and three-dimensional models of the CHS-L2/3/5/6/9/10/13 proteins showed slight variation. These results indicated that the rich anthraquinones in S. tora seeds might be ascribed to the PKSs' expansion from tandem duplication, and the seven key CHS-L2/3/5/6/9/10/13 genes provide candidate genes for further research. Our study provides an important basis for further research on the regulation of anthraquinones' biosynthesis in S. tora.

Ancient Darwinian replicators nested within eubacterial genomes

Integrative mobile genetic elements (MGEs), such as transposons and insertion sequences, propagate within bacterial genomes, but persistence times in individual lineages are short. For long-term survival, MGEs must continuously invade new hosts by horizontal transfer. Theoretically, MGEs that persist for millions of years in single lineages, and are thus subject to vertical inheritance, should not exist. Here we draw attention to an exception - a class of MGE termed REPIN. REPINs are non-autonomous MGEs whose duplication depends on non-jumping RAYT transposases. Comparisons of REPINs and typical MGEs show that replication rates of REPINs are orders of magnitude lower, REPIN population size fluctuations correlate with changes in available genome space, REPIN conservation depends on RAYT function, and REPIN diversity accumulates within host lineages. These data lead to the hypothesis that REPINs form enduring, beneficial associations with eubacterial chromosomes. Given replicative nesting, our hypothesis predicts conflicts arising from the diverging effects of selection acting simultaneously on REPINs and host genomes. Evidence in support comes from patterns of REPIN abundance and diversity in two distantly related bacterial species. Together this bolsters the conclusion that REPINs are the genetic counterpart of mutualistic endosymbiotic bacteria.

Very Low Rates of Spontaneous Gene Deletions and Gene Duplications in Dictyostelium discoideum

The study of spontaneous mutation rates has revealed a wide range of heritable point mutation rates across species, but there are comparatively few estimates for large-scale deletion and duplication rates. The handful of studies that have directly calculated spontaneous rates of deletion and duplication using mutation accumulation lines have estimated that genes are duplicated and deleted at orders of magnitude greater rates than the spontaneous point mutation rate. In our study, we tested whether spontaneous gene deletion and gene duplication rates are also high in Dictyostelium discoideum, a eukaryote with among the lowest point mutation rates (2.5 × 10-11 per site per generation) and an AT-rich genome (GC content of 22%). We calculated mutation rates of gene deletions and duplications using whole-genome sequencing data originating from a mutation accumulation experiment and determined the association between the copy number mutations and GC content. Overall, we estimated an average of 3.93 × 10-8 gene deletions and 1.18 × 10-8 gene duplications per gene per generation. While orders of magnitude greater than their point mutation rate, these rates are much lower compared to gene deletion and duplication rates estimated from mutation accumulation lines in other organisms (that are on the order of ~ 10-6 per gene/generation). The deletions and duplications were enriched in regions that were AT-rich even compared to the genomic background, in contrast to our expectations if low GC content was contributing to low mutation rates. The low deletion and duplication mutation rates in D. discoideum compared to other eukaryotes mirror their low point mutation rates, supporting previous work suggesting that this organism has high replication fidelity and effective molecular machinery to avoid the accumulation of mutations in their genome.

Differential regulation of the Tor gene homolog drives the red/green pigmentation phenotype in the aphid Myzuspersicae

In some aphid species, intraspecific variation in body colour is caused by differential carotenoid content: whilst green aphids contain only yellow carotenoids (β-, γ-, and β,γ-carotenes), red aphids additionally possess red carotenoids (torulene and 3,4-didehydrolycopene). Unusually, within animals who typically obtain carotenoids from their diet, ancestral horizontal gene transfer of carotenoid biosynthetic genes from fungi (followed by gene duplication), have imbued aphids with the intrinsic gene repertoire necessary to biosynthesise carotenoids. In the pea aphid, Acyrthosiphon pisum a lycopene (phytoene) desaturase gene (Tor) underpins the red/green phenotype, with this locus present in heterozygous form in red individuals but absent in green aphids, resulting in them being unable to convert lycopene into the red compounds 3,4-didehydrolycopene and torulene. The green peach aphid, Myzus persicae, separated from the pea aphid for ≈45MY also exists as distinct colour variable morphs, with both red and green individuals present. Here, we examined genomic data for both red and green morphs of M. persicae and identified an enlarged (compared to A. pisum) repertoire of 16 carotenoid biosynthetic genes (11 carotenoid desaturases and five carotenoid cyclase/synthase genes). From these, we identify the homolog of A. pisum Tor (here called carotene desaturase 2 or CDE-2) and show through 3D modelling that this homolog can accommodate the torulene precursor lycopene and, through RNA knockdown feeding experiments, demonstrate that disabling CDE-2 expression in red M. persicae clones results in green-coloured offspring. Unlike in A. pisum, we show that functional CDE-2 is present in the genomes of both red and green aphids. However, expression differences between the two colour morphs (350-700 fold CDE-2 overexpression in red clones), potentially driven by variants identified in upstream putative regulatory elements, underpin this phenotype. Thus, whilst aphids have a common origin of their carotenoid biosynthetic pathway, two aphid species separated for over 40MY have evolved very different drivers of intraspecific colour variation.