Mechanisms of gene duplication and amplification

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.

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.

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.

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.

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.

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.

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.

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.

MET Amplification as a Resistance Driver to TKI Therapies in Lung Cancer: Clinical Challenges and Opportunities

Targeted therapy has emerged as an important pillar for the standard of care in oncogene-driven non-small cell lung cancer (NSCLC), which significantly improved outcomes of patients whose tumors harbor oncogenic driver mutations. However, tumors eventually develop resistance to targeted drugs, and mechanisms of resistance can be diverse. MET amplification has been proven to be a driver of resistance to tyrosine kinase inhibitor (TKI)-treated advanced NSCLC with its activation of EGFR, ALK, RET, and ROS-1 alterations. The combined therapy of MET-TKIs and EGFR-TKIs has shown outstanding clinical efficacy in EGFR-mutated NSCLC with secondary MET amplification-mediated resistance in a series of clinical trials. In this review, we aimed to clarify the underlying mechanisms of MET amplification-mediated resistance to tyrosine kinase inhibitors, discuss the ways and challenges in the detection and diagnosis of MET amplifications in patients with metastatic NSCLC, and summarize the recently published clinical data as well as ongoing trials of new combination strategies to overcome MET amplification-mediated TKI resistance.

The road less travelled? Exploring the nuanced evolutionary consequences of duplicated genes

Duplicated genes have long been appreciated as both substrates and catalysts of evolutionary processes. From even the simplest cell to complex multicellular animals and plants, duplicated genes have made immeasurable contributions to the phenotypic evolution of all life on Earth. Not merely drivers of morphological innovation and speciation events, however, gene duplications sculpt the evolution of genetic architecture in ways we are only just coming to understand now we have the experimental tools to do so. As such, the present article revisits our understanding of the ways in which duplicated genes evolve, examining closely the various fates they can adopt in light of recent work that yields insights from studies of paralogues from across the tree of life that challenge the classical framework.

Architectural groups of a subtelomeric gene family evolve along distinct paths in Candida albicans

Subtelomeres are dynamic genomic regions shaped by elevated rates of recombination, mutation, and gene birth/death. These processes contribute to formation of lineage-specific gene family expansions that commonly occupy subtelomeres across eukaryotes. Investigating the evolution of subtelomeric gene families is complicated by the presence of repetitive DNA and high sequence similarity among gene family members that prevents accurate assembly from whole genome sequences. Here, we investigated the evolution of the telomere-associated (TLO) gene family in Candida albicans using 189 complete coding sequences retrieved from 23 genetically diverse strains across the species. Tlo genes conformed to the 3 major architectural groups (α/β/γ) previously defined in the genome reference strain but significantly differed in the degree of within-group diversity. One group, Tloβ, was always found at the same chromosome arm with strong sequence similarity among all strains. In contrast, diverse Tloα sequences have proliferated among chromosome arms. Tloγ genes formed 7 primary clades that included each of the previously identified Tloγ genes from the genome reference strain with 3 Tloγ genes always found on the same chromosome arm among strains. Architectural groups displayed regions of high conservation that resolved newly identified functional motifs, providing insight into potential regulatory mechanisms that distinguish groups. Thus, by resolving intraspecies subtelomeric gene variation, it is possible to identify previously unknown gene family complexity that may underpin adaptive functional variation.

Management and Treatment of Non-small Cell Lung Cancer with MET Alteration and Mechanisms of Resistance

OPINION STATEMENT

MET-driven tumors are a heterogenous group of non-small cell lung cancers (NSCLC) with activating mutations. Pathologic activation of MET can be achieved with increased number of gene copies overexpression, or decreased protein degradation through several mechanisms, including mutations, amplifications, or fusions. Besides its role as primary driver, MET activation might also mediate resistance to kinase inhibitors in NSCLC with various other actionable alterations. While checkpoint inhibitors have modest efficacy in MET-driven tumors, several approaches of targeted blockade are available. Among them the most promising are small tyrosine kinase inhibitors, antibody-drug conjugates, and bispecific antibodies. Unfortunately, resistance is virtually inevitable. Resistance to small kinase inhibitors might be mediated by kinase domain mutations or activation of shunting cascades. Various resistance mechanisms might be present in one patient, making it overcoming an unresolved problem.

Comparative genomics of Leuconostoc lactis strains isolated from human gastrointestinal system and fermented foods microbiomes

Background

Leuconostoc lactis forms a crucial member of the genus Leuconostoc and has been widely used in the fermentation industry to convert raw material into acidified and flavored products in dairy and plant-based food systems. Since the ecological niches that strains of Ln. lactis being isolated from were truly diverse such as the human gut, dairy, and plant environments, comparative genome analysis studies are needed to better understand the strain differences from a metabolic adaptation point of view across diverse sources of origin. We compared eight Ln. lactis strains of 1.2.28, aa_0143, BIOML-A1, CBA3625, LN19, LN24, WIKIM21, and WiKim40 using bioinformatics to elucidate genomic level characteristics of each strain for better utilization of this species in a broad range of applications in food industry.

Results

Phylogenomic analysis of twenty-nine Ln. lactis strains resulted in nine clades. Whole-genome sequence analysis was performed on eight Ln. lactis strains representing human gastrointestinal tract and fermented foods microbiomes. The findings of the present study are based on comparative genome analysis against the reference Ln. lactis CBA3625 genome. Overall, a ~ 41% of all CDS were conserved between all strains. When the coding sequences were assigned to a function, mobile genetic elements, mainly insertion sequences were carried by all eight strains. All strains except LN24 and WiKim40 harbor at least one intact putative prophage region, and two of the strains contained CRISPR-Cas system. All strains encoded Lactococcin 972 bacteriocin biosynthesis gene clusters except for CBA3625.

Conclusions

The findings in the present study put forth new perspectives on genomics of Ln. lactis via complete genome sequence based comparative analysis and further determination of genomic characteristics. The outcomes of this work could potentially pave the way for developing elements for future strain engineering applications.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-022-01074-6.

Global analysis of common bean multidrug and toxic compound extrusion transporters (PvMATEs): PvMATE8 and pinto bean seed coat darkening

In common bean (Phaseolus vulgaris L.), postharvest seed coat darkening is an undesirable trait that affects crop value. The increased accumulation of proanthocyanidins (PAs) in the seed coat results in darker seeds in many market classes of colored beans after harvest. The precursors of PAs are synthesized in the cytoplasm, and subsequently get glycosylated and then transported to the vacuoles where polymerization occurs. Thus, vacuolar transporters play an important role in the accumulation of PAs. Here, we report that common bean genome contains 59 multidrug and toxic compound extrusion genes (PvMATEs). Phylogenetic analysis of putative PvMATEs with functionally characterized MATEs from other plant species categorized them into substrate-specific clades. Our data demonstrate that a vacuolar transporter PvMATE8 is expressed at a higher level in the pinto bean cultivar CDC Pintium (regular darkening) compared to 1533-15 (slow darkening). PvMATE8 localizes in the vacuolar membrane and rescues the PA deficient (tt12) mutant phenotype in Arabidopsis thaliana. Analysis of PA monomers in transgenic seeds together with wild-type and mutants suggests a possible feedback regulation of PA biosynthesis and accumulation. Identification of PvMATE8 will help better understand the mechanism of PA accumulation in common bean.

B cell lymphoma 6A regulates immune development and function in zebrafish

BCL6A is a transcriptional repressor implicated in the development and survival of B and T lymphoctyes, which is also highly expressed in many non-Hodgkin’s lymphomas, such as diffuse large B cell lymphoma and follicular lymphoma. Roles in other cell types, including macrophages and non-hematopoietic cells, have also been suggested but require further investigation. This study sought to identify and characterize zebrafish BCL6A and investigate its role in immune cell development and function, with a focus on early macrophages. Bioinformatics analysis identified a homologue for BCL6A (bcl6aa), as well as an additional fish-specific duplicate (bcl6ab) and a homologue for the closely-related BCL6B (bcl6b). The human BCL6A and zebrafish Bcl6aa proteins were highly conserved across the constituent BTB/POZ, PEST and zinc finger domains. Expression of bcl6aa during early zebrafish embryogenesis was observed in the lateral plate mesoderm, a site of early myeloid cell development, with later expression seen in the brain, eye and thymus. Homozygous bcl6aa mutants developed normally until around 14 days post fertilization (dpf), after which their subsequent growth and maturation was severely impacted along with their relative survival, with heterozygous bcl6aa mutants showing an intermediate phenotype. Analysis of immune cell development revealed significantly decreased lymphoid and macrophage cells in both homozygous and heterozygous bcl6aa mutants, being exacerbated in homozygous mutants. In contrast, the number of neutrophils was unaffected. Only the homozygous bcl6aa mutants showed decreased macrophage mobility in response to wounding and reduced ability to contain bacterial infection. Collectively, this suggests strong conservation of BCL6A across evolution, including a role in macrophage biology.

Tandem duplication of a genomic region encoding glutathione S-transferase epsilon-2 and -4 genes in DDT-resistant Anopheles stephensi strain from India

The glutathione S-transferases (GST) genes are a multigene family of enzymes involved in the metabolism of endogenous and xenobiotic compounds by catalysing the conjugation of the reduced form of glutathione to the substrate. The epsilon class of GST (GSTe), unique to arthropods, is known to be involved in the detoxification process of several classes of insecticides, and GSTe2 in particular is known to have DDT dehydrochlorinase activity. This communication reports a tandem duplication of a genomic region encoding GSTe2 and GSTe4 genes in a laboratory-colonized DDT-resistant Anopheles stephensi. We identified duplication breakpoints and the organization of gene duplication through Sanger sequencing performed on long-PCR products. Manual annotation of sequences revealed a tandemly-arrayed duplication of a 3.62 kb segment of GST epsilon gene clusters comprised of five genes: a partial GSTe1, GSTe2, GSTe2-pseudogene, GSTe4 and partial GSTe5, interconnected by a conserved 2.42 kb DNA insert segment major part of which is homologous to a genomic region located on a different chromosome. The tandemly duplicated array contained a total of two GSTe2 and three GSTe4 functional paralog genes. Read-depth coverage and split-read analysis of Illumina-based whole-genome sequence reads confirmed the presence of duplication in the corresponding region of the genome. The increased gene dose in mosquitoes as a result of the GSTe gene-duplication may be an adaptive process to increase levels of detoxifying enzymes to counter insecticide pressure.

A positive feedback loop between two C-type lectins originated from gene duplication and relish promotes the expression of antimicrobial peptides in Procambarus clarkii

Gene duplication (GD) leads to the expansion of gene families that contributes organisms adapting to stress or environment and dealing with the infection of various pathogens. C-type lectins (CTLs) in crustaceans undergo gene expansion and participate in various immune responses. However, the functions of different CTL produced by GD are not fully characterized. In the present study, two CTL genes (designated as PcLec-EPS and PcLec-QPS, respectively) were identified from Procambarus clarkii. PcLec-EPS and PcLec-QPS originate from GD and the main difference between them is exon 3. PcLec-EPS and PcLec-QPS respectively contains EPS and QPS motif in their carbohydrate recognition domain. The mRNA levels of PcLec-EPS and PcLec-QPS in hemocytes, gills, intestine and lymph underwent time-dependent enhancement after D-Mannose and D-Galactose challenge. Recombinant PcLec-EPS and PcLec-QPS could bind to carbohydrates and microbes, and agglutinate bacteria. The results of experiments on recombinant protein injection and RNA interference indicate that PcLec-EPS and PcLec-QPS can respectively strong recognize and bind D-Mannose and D-Galactose, activate the Relish transcriptional factor, and further upregulate the expression of different antimicrobial peptides (AMPs). In addition, these two CTLs and Relish could positively regulate the expression of each other, suggesting that there is a positive feedback loop between two CTLs and Relish that regulates the expression of AMPs. It may contribute to the expansion of the immune response for host quickly and efficiently eliminating pathogenic microorganisms. This study provides new knowledge for clear understanding the significance and function of different CTL generated by GD in immune defenses in crustacean.

Low levels of tetracyclines select for a mutation that prevents the evolution of high-level resistance to tigecycline

In a collection of Escherichia coli isolates, we discovered a new mechanism leading to frequent and high-level tigecycline resistance involving tandem gene amplifications of an efflux pump encoded by the tet(A) determinant. Some isolates, despite carrying a functional tet(A), could not evolve high-level tigecycline resistance by amplification due to the presence of a deletion in the TetR(A) repressor. This mutation impaired induction of tetA(A) (encoding the TetA(A) efflux pump) in presence of tetracyclines, with the strongest effect observed for tigecycline, subsequently preventing the development of tet(A) amplification-dependent high-level tigecycline resistance. We found that this mutated tet(A) determinant was common among tet(A)-carrying E. coli isolates and analysed possible explanations for this high frequency. First, while the mutated tet(A) was found in several ST-groups, we found evidence of clonal spread among ST131 isolates, which increases its frequency within E. coli databases. Second, evolution and competition experiments revealed that the mutation in tetR(A) could be positively selected over the wild-type allele at sub-inhibitory concentrations of tetracyclines. Our work demonstrates how low concentrations of tetracyclines, such as those found in contaminated environments, can enrich and select for a mutation that generates an evolutionary dead-end that precludes the evolution towards high-level, clinically relevant tigecycline resistance.

A study on evolution of antimicrobial resistance reveals how sub-inhibitory concentrations of antibiotics enrich and select for a mutated allele that prevents evolution towards clinically significant levels of antibiotic resistance.

[Research Progress of Acquired Resistance Mediated by MET Amplification in Advanced Non-small Cell Lung Cancer]

Mesenchymal-epithelial transition factor (MET) amplification is an important driver of resistance in epidermal growth factor receptor (EGFR)-mutant non-small cell lung cancer (NSCLC), and the combination of MET proto-oncogene (MET) and EGFR-tyrosine kinase inhibitors (TKIs) has shown promise in overcoming this molecularly defined acquired resistance. Emerging data also demonstrate MET amplification as a resistance driver to TKIs-treated anaplastic lymphoma kinase (ALK)-, RET-, and ROS1-fusion NSCLC. Here, we review the literature on recent research progress of MET amplification as a resistance driver to targeted therapy in oncogene-driven NSCLC and summarize the progress of clinical strategies to overcome the resistance mechanism.
.

MDR Pumps as Crossroads of Resistance: Antibiotics and Bacteriophages

At present, antibiotic resistance represents a global problem in modern medicine. In the near future, humanity may face a situation where medicine will be powerless against resistant bacteria and a post-antibiotic era will come. The development of new antibiotics is either very expensive or ineffective due to rapidly developing bacterial resistance. The need to develop alternative approaches to the treatment of bacterial infections, such as phage therapy, is beyond doubt. The cornerstone of bacterial defense against antibiotics are multidrug resistance (MDR) pumps, which are involved in antibiotic resistance, toxin export, biofilm, and persister cell formation. MDR pumps are the primary non-specific defense of bacteria against antibiotics, while drug target modification, drug inactivation, target switching, and target sequestration are the second, specific line of their defense. All bacteria have MDR pumps, and bacteriophages have evolved along with them and use the bacteria’s need for MDR pumps to bind and penetrate into bacterial cells. The study and understanding of the mechanisms of the pumps and their contribution to the overall resistance and to the sensitivity to bacteriophages will allow us to either seriously delay the onset of the post-antibiotic era or even prevent it altogether due to phage-antibiotic synergy.

Comparative Genomics Reveals Insights into the Divergent Evolution of Astigmatic Mites and Household Pest Adaptations

Highly diversified astigmatic mites comprise many medically important human household pests such as house dust mites causing ∼1–2% of all allergic diseases globally; however, their evolutionary origin and diverse lifestyles including reversible parasitism have not been illustrated at the genomic level, which hampers allergy prevention and our exploration of these household pests. Using six high-quality assembled and annotated genomes, this study not only refuted the monophyly of mites and ticks, but also thoroughly explored the divergence of Acariformes and the diversification of astigmatic mites. In monophyletic Acariformes, Prostigmata known as notorious plant pests first evolved, and then rapidly evolving Astigmata diverged from soil oribatid mites. Within astigmatic mites, a wide range of gene families rapidly expanded via tandem gene duplications, including ionotropic glutamate receptors, triacylglycerol lipases, serine proteases and UDP glucuronosyltransferases. Gene diversification after tandem duplications provides many genetic resources for adaptation to sensing environmental signals, digestion, and detoxification in rapidly changing household environments. Many gene decay events only occurred in the skin-burrowing parasitic mite Sarcoptes scabiei. Throughout the evolution of Acariformes, massive horizontal gene transfer events occurred in gene families such as UDP glucuronosyltransferases and several important fungal cell wall lytic enzymes, which enable detoxification and digestive functions and provide perfect drug targets for pest control. This comparative study sheds light on the divergent evolution and quick adaptation to human household environments of astigmatic mites and provides insights into the genetic adaptations and even control of human household pests.

Recurrent Duplication and Diversification of Acrosomal Fertilization Proteins in Abalone

Reproductive proteins mediating fertilization commonly exhibit rapid sequence diversification driven by positive selection. This pattern has been observed among nearly all taxonomic groups, including mammals, invertebrates, and plants, and is remarkable given the essential nature of the molecular interactions mediating fertilization. Gene duplication is another important mechanism that facilitates the generation of molecular novelty through functional divergence. Following duplication, paralogs may partition ancestral gene function (subfunctionalization) or acquire new roles (neofunctionalization). However, the contributions of duplication followed by sequence diversification to the molecular diversity of gamete recognition genes has been understudied in many models of fertilization. The marine gastropod mollusk abalone is a classic model for fertilization. Its two acrosomal proteins (lysin and sp18) are ancient gene duplicates with unique gamete recognition functions. Through detailed genomic and bioinformatic analyses we show how duplication events followed by sequence diversification has played an ongoing role in the evolution of abalone acrosomal proteins. The common ancestor of abalone had four members of its acrosomal protein family in a tandem gene array that repeatedly experienced positive selection. We find that both sp18 paralogs contain positively selected sites located in different regions of the paralogs, suggestive of functional divergence where selection acted upon distinct binding interfaces in each paralog. Further, a more recent species-specific duplication of both lysin and sp18 in the European abalone H. tuberculata is described. Despite clade-specific acrosomal protein paralogs, there are no concomitant duplications of egg coat proteins in H. tuberculata, indicating that duplication of egg proteins per se is not responsible for retention of duplicated acrosomal proteins. We hypothesize that, in a manner analogous to host/pathogen evolution, sperm proteins are selected for increased diversity through extensive sequence divergence and recurrent duplication driven by conflict mechanisms.

Evolution of rhizobial symbiosis islands through insertion sequence-mediated deletion and duplication

Symbiosis between organisms influences their evolution via adaptive changes in genome architectures. Immunity of soybean carrying the Rj2 allele is triggered by NopP (type III secretion system [T3SS]-dependent effector), encoded by symbiosis island A (SymA) in B. diazoefficiens USDA122. This immunity was overcome by many mutants with large SymA deletions that encompassed T3SS (rhc) and N2 fixation (nif) genes and were bounded by insertion sequence (IS) copies in direct orientation, indicating homologous recombination between ISs. Similar deletion events were observed in B. diazoefficiens USDA110 and B. japonicum J5. When we cultured a USDA122 strain with a marker gene sacB inserted into the rhc gene cluster, most sucrose-resistant mutants had deletions in nif/rhc gene clusters, similar to the mutants above. Some deletion mutants were unique to the sacB system and showed lower competitive nodulation capability, indicating that IS-mediated deletions occurred during free-living growth and the host plants selected the mutants. Among 63 natural bradyrhizobial isolates, 2 possessed long duplications (261-357 kb) harboring nif/rhc gene clusters between IS copies in direct orientation via homologous recombination. Therefore, the structures of symbiosis islands are in a state of flux via IS-mediated duplications and deletions during rhizobial saprophytic growth, and host plants select mutualistic variants from the resultant pools of rhizobial populations. Our results demonstrate that homologous recombination between direct IS copies provides a natural mechanism generating deletions and duplications on symbiosis islands.

Reconstructing Gene Gains and Losses with BadiRate

Estimating gene gain and losses is paramount to understand the molecular mechanisms underlying adaptive evolution. Despite the advent of high-throughput sequencing, such analyses have been so far hampered by the poor contiguity of genome assemblies. The increasing affordability of long-read sequencing technologies will however revolutionize our capacity to identify gene gains and losses at an unprecedented resolution, even in non-model organisms. To thoroughly exploit all such multigene family variation, the software BadiRate implements a collection of birth-and-death stochastic models, aiming at estimating by maximum likelihood the gene turnover rates along the internal and external branches of a given phylogenetic species tree. Its statistical framework also provides versatility for inferring the gene family content at the internal phylogenetic nodes (and to estimate the minimum number of gene gains and losses in each branch), for statistically contrasting competing hypotheses (e.g., accelerations of the gene turnover rates at pre-defined clades), and for pinpointing gene family expansions or contractions likely driven by natural selection. In this chapter we review the theoretical models implemented in BadiRate and illustrate their applicability by analyzing a hypothetical data set of 14 microbial species.

Multidrug Resistance Pumps as a Keystone of Bacterial Resistance

Antibiotic resistance is a global problem of modern medicine. A harbinger of the onset of the postantibiotic era is the complexity and high cost of developing new antibiotics as well as their inefficiency due to the rapidly developing resistance of bacteria. Multidrug resistance (MDR) pumps, involved in the formation of resistance to xenobiotics, the export of toxins, the maintenance of cellular homeostasis, and the formation of biofilms and persistent cells, are the keystone of bacterial protection against antibiotics. MDR pumps are the basis for the nonspecific protection of bacteria, while modification of the drug target, inactivation of the drug, and switching of the target or sequestration of the target is the second specific line of their protection. Thus, the nonspecific protection of bacteria formed by MDR pumps is a barrier that prevents the penetration of antibacterial substances into the cell, which is the main factor determining the resistance of bacteria. Understanding the mechanisms of MDR pumps and a balanced assessment of their contribution to total resistance, as well as to antibiotic sensitivity, will either seriously delay the onset of the postantibiotic era or prevent its onset in the foreseeable future.