Search

Quantitative determination of neutralizing antibodies against Severe Acute Respiratory Syndrome Corona Virus-2 (SARS-CoV-2) is paramount in immunodiagnostics, vaccine efficacy testing, and immune response profiling among the vaccinated population. Cost-effective, rapid, easy-to-perform assays are essential to support the vaccine development process and immunosurveillance studies. We describe a bead-based screening assay for S1-neutralization using recombinant fluorescent proteins of hACE2 and SARS-CoV2-S1, immobilized on solid beads employing nanobodies/metal-affinity tags. Nanobody-mediated capture of SARS-CoV-2-Spike (S1) on agarose beads served as the trap for soluble recombinant ACE2-GFPSpark, inhibited by neutralizing antibody. The first approach demonstrates single-color fluorescent imaging of ACE2-GFPSpark binding to His-tagged S1-Receptor Binding Domain (RBD-His) immobilized beads. The second approach is dual-color imaging of soluble ACE2-GFPSpark to S1-Orange Fluorescent Protein (S1-OFPSpark) beads. Both methods showed a good correlation with the gold standard pseudovirion assay and can be adapted to any fluorescent platforms for screening.

Protein misfolding and aggregation play a significant role in several neurodegenerative diseases. In the present work, the spontaneous aggregation of hen egg-white lysozyme (HEWL) in an alkaline pH 12.2 at an ambient temperature was studied to obtain molecular insights. The time-dependent changes in spectral peaks indicated the formation of β sheets and their effects on the backbone and amino acids during the aggregation process. Introducing iodoacetamide revealed the crucial role of intermolecular disulphide bonds amidst monomers in the aggregation process. These findings were corroborated by Molecular Dynamics (MD) simulations and protein-docking studies. MD simulations helped establish and visualize the unfolding of the proteins when exposed to an alkaline pH. Protein docking revealed a preferential dimer formation between the HEWL monomers at pH 12.2 compared with the neutral pH. The combination of Raman spectroscopy and MD simulations is a powerful tool to study protein aggregation mechanisms.

Since its advent in December 2019, SARS-CoV-2 has diverged into multiple variants with differing levels of virulence owing to the accumulation of mutations in its genome. The structural changes induced by non-synonymous mutations in major drug targets of the virus are known to alter the binding of potential antagonistic inhibitors. Here, we analyzed the effects of non-synonymous mutations in major targets of SARS-CoV-2 in response to potential peptide inhibitors. We screened 12 peptides reported to have anti-viral properties against RBD and 5 peptides against Mpro of SARS-CoV-2 variants using molecular docking and simulation approaches. The mutational landscape of RBD among SARS-CoV-2 variants had 21 non-synonymous mutations across 18 distinct sites. Among these, 14 mutations were present in the RBM region directly interacting with the hACE2 receptor. However, Only 3 non-synonymous mutations were observed in Mpro. We found that LCB1 - a de novo-synthesized peptide has the highest binding affinity to RBD despite non-synonymous mutations in variants and engages key residues of RBD-hACE2 interaction such as K417, E484, N487, and N501. Similarly, an antimicrobial peptide; 2JOS, was identified against Mpro with high binding affinity as it interacts with key residues in dimerization sites such as E166 and F140 crucial for viral replication. MD simulations affirm the stability of RBD-LCB1 and Mpro-2JOS complexes with an average RMSD of 1.902 and 2.476 respectively. We ascertain that LCB1 and 2JOS peptides are promising inhibitors to combat emerging variants of SARS-CoV-2 and thus warrant further investigations using in-vitro and in-vivo analysis. Communicated by Ramaswamy H. Sarma

Intra-retinal axon guidance involves a coordinated expression of transcription factors, axon guidance genes, and secretory molecules within the retina. Pax6, the master regulator gene, has a spatio-temporal expression typically restricted till neurogenesis and fate-specification. However, our observation of persistent expression of Pax6 in mature RGCs led us to hypothesize that Pax6 could play a major role in axon guidance after fate specification. Here, we found significant alteration in intra-retinal axon guidance and fasciculation upon knocking out of Pax6 in E15.5 retina. Through unbiased transcriptome profiling between Pax6fl/fl and Pax6−/− retinas, we revealed the mechanistic insight of its role in axon guidance. Our results showed a significant increase in the expression of extracellular matrix molecules and decreased expression of retinal fate specification and neuron projection guidance molecules. Additionally, we found that EphB1 and Sema5B are directly regulated by Pax6 owing to the guidance defects and improper fasciculation of axons. We conclude that Pax6 expression post fate specification of RGCs is necessary for regulating the expression of axon guidance genes and most importantly for maintaining a conducive ECM through which the nascent axons get guided and fasciculate to reach the optic disc.

Multidrug-resistant Staphylococcus aureus is a leading concern worldwide. Coagulase-Negative Staphylococci are claimed to be the reservoir and source of important resistant elements in S. aureus. However, the origin and evolutionary route of resistant genes in S. aureus are still remaining unknown. Here, we performed a detailed phylogenomic analysis of 152 completely sequenced S. aureus strains in comparison with 7,529 non-Staphylococcus aureus reference bacterial genomes. Our results reveal that S. aureus has a large open pan-genome where 97 (55%) of its known resistant-related genes belonging to its accessory genome. Among these genes, 47 (27%) were located within the Staphylococcal Cassette Chromosome mec (SCCmec), a transposable element responsible for resistance against major classes of antibiotics including beta-lactams, macrolides, and aminoglycosides. However, the physically linked mec-box genes (MecA–MecR–MecI) that are responsible for the maintenance of SCCmec elements is not unique to S. aureus, instead it is widely distributed within Staphylococcaceae family. The phyletic patterns of SCCmec-encoded resistant genes in Staphylococcus species are significantly different from that of its core genes indicating frequent exchange of these genes between Staphylococcus species. Our in-depth analysis of SCCmec-resistant gene phylogenies reveals that genes such as blaZ, ble, kmA, and tetK that are responsible for beta-lactam, bleomycin, kanamycin, and tetracycline resistance in S. aureus were laterally transferred from non-Staphylococcus sources. In addition, at least 11 non-SCCmec-encoded resistant genes in S. aureus, were laterally acquired from distantly related species. Our study evidently shows that gene transfers played a crucial role in shaping the evolution of antibiotic resistance in S. aureus.

While the immunomodulatory pathways initiated in immune cells contribute to therapeutic response, their activation in cancer cells play a role in cancer progression. Also, many of the aberrantly expressed immunomodulators on cancer cells are considered as therapeutic targets. Here, we introduce host defense peptide (HDP), a known immuomodulator, as a therapeutic agent to target them. The cationic host defense peptides (HDPs), an integral part of the innate immune system, possess membranolytic activity, which imparts antimicrobial and antitumor efficacy to it. They act as immunomodulators by activating the immune cells. Though their antimicrobial function has been recently reassigned to immunoregulation, their antitumor activity is still attributed to its membranolytic activity. This membrane pore formation ability, which is proportional to the concentration of the peptide, also leads to side effects like hemolysis, limiting their therapeutic application. So, despite the identification of a variety of anticancer HDPs, their clinical utility is limited. Though HDPs are shown to exert the immunomodulatory activity through specific membrane targets on immune cells, their targets on cancer cells are unknown. We show that SSTP1, a novel HDP identified by shotgun cloning, binds to the active IL6/IL6Rα/gp130 complex on cancer cells, rearranging the active site residues. In contrast to the IL6 blockers inhibiting JAK/STAT activity, SSTP1 shifts the proliferative IL6/JAK/STAT signaling to the apoptotic IL6/JNK/AP1 pathway. In IL6Rα-overexpressing cancer cells, SSTP1 induces apoptosis at low concentration through JNK pathway, without causing significant membrane disruption. We highlight the importance of immunomodulatory pathways in cancer apoptosis, apart from its established role in immune cell regulation and cancer cell proliferation. Our study suggests that identification of the membrane targets for the promising anticancer HDPs might lead to the identification of new drugs for targeted therapy.

Quantitative determination of neutralizing antibodies against Severe Acute Respiratory Syndrome Corona Virus-2 (SARS-CoV-2) is paramount in immunodiagnostics, vaccine efficacy testing, and immune response profiling among the vaccinated population. Cost-effective, rapid, easy-to-perform assays are essential to support the vaccine development process and immunosurveillance studies. Here, we describe a bead-based screening assay for S1-neutralization using recombinant fluorescent proteins of hACE2 and SARS-CoV2-S1, immobilized on solid beads employing nanobodies /metal-affinity tags. Nanobody-mediated capture of SARS-CoV-2 - Spike (S1) on agarose beads served as the trap for soluble recombinant ACE2-GFPSpark, inhibited by neutralizing antibody. The first approach demonstrates single-color fluorescent imaging of ACE2–GFPspark binding to His-tagged S1-Receptor Binding Domain (RBD-His) immobilized beads. The second approach is dual-color imaging of soluble ACE2-GFPSpark to S1-Orange Fluorescent Protein (S1-OFPSpark) beads. Both methods showed a good correlation with the gold standard pseudovirion assay and can be adapted to any fluorescent platforms for screening. Life-time imaging of the ACE2-GFPSpark confirmed the interaction of ACE2 and S1-OFPSpark on beads. The self-renewable source of secreted recombinant proteins from stable cells and its direct use without necessitating purification renders the platform a cost-effective and rapid one than the popular pseudovirion assay and live virus-based assays. Any laboratory with minimum expertise can rapidly perform this bead assay for neutralizing antibody detection using stable engineered cells.

Background: The anticancer activities of host defense peptides (HDPs) are mainly attributed to their cell penetrating activity. Accordingly, several approaches based on machine learning, calculating the membrane pore formation ability, have been proposed for the anti-cancer peptide identification. Since the membranolytic activity can lead to nonspecific effects, like hemolysis, the therapeutic application of such molecules is limited. In this context, we considered the immunomodulatory activity of HDPs, which is regulated by specific membrane targets and immunomodulatory pathways in immune cells. As many immunomodulators are aberrantly expressed in cancer cells, the possibility of the activation of an immunomodulatory pathway was investigated for a novel anticancer HDP, SSTP1.Methods: The fourteen mature peptides identified by shotgun cloning from the frog-skin secretions were screened for cytotoxicity in oral cancer cells. The mechanism of action of the selected peptide, SSTP1 was investigated in comparison to its inactive mutant, SSTP2. An RNA-Seq coupled with pathway enrichment analysis was performed to identify the upstream signaling leading to the mitochondrial pathway of apoptosis. Since the activation of the IL6/IL6R pathway was suggested, we performed in silico docking studies to find the binding of SSTP1 to the IL6/IL6Rα/gp130 complex. The dynamic simulation predicted the conformational changes in the active site residues. The confocal co-localization studies, pull-down assay, FRET analysis, western blot and reporter assays were performed to elucidate the role of the IL6/IL6R pathway in SSTP1-induced apoptosis. Specific small molecule inhibitors and neutralizing antibodies were used to ascertain the role of the IL6/IL6Rα/gp130 complex in the activation of JNK/AP1 pathway-dependent cell death.Results: SSTP1, a novel temporin, modulates the IL6 pathway and induces apoptosis when it binds to IL6Rα on the active IL6/IL6Rα/gp130 complex, rearranging the active site residues. In contrast to the IL6 blockers inhibiting JAK/STAT activity, SSTP1 shifts the proliferative IL6/JAK/STAT signaling to the apoptotic IL6/JNK/AP1 pathway. In IL6Rα-overexpressing cancer cells, apoptosis is preferred over the membranolytic activity, upon SSTP1 treatment.Conclusions: Here, we provide the evidence of an HDP-induced signaling through immunomodulators, leading to apoptosis in cancer cells. Our study also implies the importance of identifying the targets of HDPs for their clinical application

Protein-protein interactions between virus and host are crucial for infection. SARS-CoV-2, the causative agent of COVID-19 pandemic is an RNA virus prone to mutations. Formation of a stable binding interface between the Spike (S) protein Receptor Binding Domain (RBD) of SARS-CoV-2 and Angiotensin-Converting Enzyme 2 (ACE2) of host actuates viral entry. Yet, how this binding interface evolves as virus acquires mutations during pandemic remains elusive. Here, using a high fidelity bioinformatics pipeline, we analysed 31,403 SARS-CoV-2 genomes across the globe, and identified 444 non-synonymous mutations that cause 49 distinct amino acid substitutions in the RBD. Molecular phylogenetic analysis suggested independent emergence of these RBD mutants during pandemic. In silico structure modelling of interfaces induced by mutations on residues which directly engage ACE2 or lie in the near vicinity revealed molecular rearrangements and binding energies unique to each RBD mutant. Comparative structure analysis using binding interface from mouse that prevents SARS-CoV-2 entry uncovered minimal molecular determinants in RBD necessary for the formation of stable interface. We identified that interfacial interaction involving amino acid residues N487 and G496 on either ends of the binding scaffold are indispensable to anchor RBD and are well conserved in all SARS-like corona viruses. All other interactions appear to be required to locally remodel binding interface with varying affinities and thus may decide extent of viral transmission and disease outcome. Together, our findings propose the modalities and variations in RBD-ACE2 interface formation exploited by SARS-CoV-2 for endurance.ImportanceCOVID-19, so far the worst hit pandemic to mankind, started in January 2020 and is still prevailing globally. Our study identified key molecular arrangements in RBD-ACE2 interface that help virus to tolerate mutations and prevail. In addition, RBD mutations identified in this study can serve as a molecular directory for experimental biologists to perform functional validation experiments. The minimal molecular requirements for the formation of RBD-ACE2 interface predicted using in silico structure models may help precisely design neutralizing antibodies, vaccines and therapeutics. Our study also proposes the significance of understanding evolution of protein interfaces during pandemic.

Metagenomic studies permit the exploration of microbial diversity in a defined habitat, and binning procedures enable phylogenomic analyses, taxon description, and even phenotypic characterizations in the absence of morphological evidence. Such lineages include asgard archaea, which were initially reported to represent archaea with eukaryotic cell complexity, although the first images of such an archaeon show simple cells with prokaryotic characteristics. However, these metagenome-assembled genomes (MAGs) might suffer from data quality problems not encountered in sequences from cultured organisms due to two common analytical procedures of bioinformatics: assembly of metagenomic sequences and binning of assembled sequences on the basis of innate sequence properties and abundance across samples. Consequently, genomic sequences of distantly related taxa, or domains, can in principle be assigned to the same MAG and result in chimeric sequences. The impacts of low-quality or chimeric MAGs on phylogenomic and metabolic prediction remain unknown. Debates that asgard archaeal data are contaminated with eukaryotic sequences are overshadowed by the lack of evidence indicating that individual asgard MAGs stem from the same chromosome. Here, we show that universal proteins including ribosomal proteins of asgard archaeal MAGs fail to meet the basic phylogenetic criterion fulfilled by genome sequences of cultured archaea investigated to date: These proteins do not share common evolutionary histories to the same extent as pure culture genomes do, pointing to a chimeric nature of asgard archaeal MAGs. Our analysis suggests that some asgard archaeal MAGs represent unnatural constructs, genome-like patchworks of genes resulting from assembly and/or the binning process.

Extreme flooding is one of the major risk factors for human health, and it can significantly influence the microbial communities and enhance the mobility of infectious disease agents within the affected areas. The flood crisis in 2018 was one of the severe natural calamities recorded in the southern state of India (Kerala) that significantly affected its economy and ecological habitat. We utilized a combination of shotgun metagenomics and bioinformatics approaches to understand the bacterial profile and the abundance of pathogenic and antibiotic-resistant bacteria in extremely flooded areas of Kuttanad, Kerala (4&ndash, 10 feet below sea level). Here we report the bacterial profiles of flooded sites that are abundant with virulent and resistant bacteria. The flooded sites were heavily contaminated with faecal contamination indicators such as Escherichia coli and Enterococcus faecalis and multidrug-resistant strains of Pseudomonas aeruginosa, Salmonella typhi/typhimurium, Klebsiella pneumoniae, Vibrio cholerae. The resistome of the flooded sites contains 103 known resistant genes, of which 38% are plasmid-encoded, where most of them are known to be associated with pathogenic bacteria. Our results reveal an overall picture of the bacterial profile and resistome of sites following a devastating flood event, which might increase the levels of pathogens and its associated risks.

The origin of mitochondria was a crucial event in eukaryote evolution. A recent report claimed to provide evidence, based on branch length variation in phylogenetic trees, that the mitochondrion came late in eukaryotic evolution. Here, we reinvestigate their claim with a reanalysis of the published data. We show that the analyses underpinning a late mitochondrial origin suffer from multiple fatal flaws founded in inappropriate statistical methods and analyses, in addition to erroneous interpretations.

Multi-drug resistant S. aureus is a leading concern worldwide. Coagulase-Negative Staphylococci (CoNS) are claimed to be the reservoir and source of important resistant elements in S. aureus. However, the origin and evolutionary route of resistant genes in S. aureus are still remaining unknown. Here, we performed a detailed phylogenomic analysis of 152 completely sequenced S. aureus strains in comparison with 7,529 non-S. aureus reference bacterial genomes. Our results reveals that S. aureus has a large open pan-genome where 97 (55%) of its known resistant related genes belonging to its accessory genome. Among these genes, 47 (27%) were located within the Staphylococcal Cassette Chromosome (SCCmec), a transposable element responsible for resistance against major classes of antibiotics including beta-lactams, macrolides and aminoglycosides. However, the physically linked mec-box genes (MecA-MecR-MecI) that are responsible for the maintenance of SCCmec elements is not unique to S. aureus, instead it is widely distributed within Staphylococcaceae family. The phyletic patterns of SCCmec encoded resistant genes in Staphylococcus species are significantly different from that of its core genes indicating frequent exchange of these genes between Staphylococcus species. Our in-depth analysis of SCCmec resistant gene phylogenies reveals that genes such as blaZ, ble, kmA and tetK that are responsible for beta-lactam, bleomycin, kanamycin and tetracycline resistance in S. aureus were laterally transferred from non-Staphylococcus sources. In addition, at least 11 non-SCCmec encoded resistant genes in S. aureus, mostly present in plasmid are laterally acquired from distantly related species. Our study evidently shows that gene transfers played a crucial role in shaping the evolution of antibiotic resistance in S. aureus.

SummaryMetagenomic studies have claimed the existence of novel lineages with unprecedented properties never before observed in prokaryotes. Such lineages include Asgard archaea1–3, which are purported to represent archaea with eukaryotic cell complexity, and the Candidate Phyla Radiation (CPR), a novel domain level taxon erected solely on the basis of metagenomic data4. However, it has escaped the attention of most biologists that these metagenomic sequences are not assembled into genomes by sequence overlap, as for cultured archaea and bacteria. Instead, short contigs are sorted into computer files by a process called binning in which they receive taxonomic assignment on the basis of sequence properties like GC content, dinucleotide frequencies, and stoichiometric co-occurrence across samples. Consequently, they are not genome sequences as we know them, reflecting the gene content of real organisms. Rather they are metagenome assembled genomes (MAGs). Debates that Asgard data are contaminated with individual eukaryotic sequences5–7are overshadowed by the more pressing issue that no evidence exists to indicate that any sequences in binned Asgard MAGs actually stem from the same chromosome, as opposed to simply stemming from the same environment. Here we show that Asgard and CPR MAGs fail spectacularly to meet the most basic phylogenetic criterion8fulfilled by genome sequences of all cultured prokaryotes investigated to date: the ribosomal proteins of Asgard and CPR MAGs do not share common evolutionary histories. Their phylogenetic behavior is anomalous to a degree never observed with genomes of real organisms. CPR and Asgard MAGs are binning artefacts, assembled from environments where up to 90% of the DNA is from dead cells9–12. Asgard and CPR MAGs are unnatural constructs, genome-like patchworks of genes that have been stitched together into computer files by binning.

Influenza A (H1N1) continues to be a major public health threat due to possible emergence of a more virulent H1N1 strain resulting from dynamic changes in virus adaptability consequent to functional mutations and antigenic drift in the hemagglutinin (HA) and neuraminidase (NA) surface proteins. In this study, we describe the genetic and evolutionary characteristics of H1N1 strains that circulated in India over a period of nine years from 2009 to 2017 in relation to global strains. The finding is important from a global perspective since previous phylogenetic studies have suggested that the tropics contributed substantially to the global circulation of influenza viruses. Bayesian phylogenic analysis of HA sequences along with global strains indicated that there is a temporal pattern of H1N1 evolution and clustering of Indian isolates with globally circulating strains. Interestingly, we observed four new amino acid substitutions (S179N, I233T, S181T and I312V) in the HA sequence of H1N1 strains isolated during 2017 and two (S181T and I312V) were found to be unique in Indian isolates. Structurally these two unique mutations could lead to altered glycan specificity of the HA gene. Similarly, sequence and structural analysis of NA domain revealed that the presence of K432E mutation in H1N1 strains isolated after 2015 from India and in global strains found to induce a major loop shift in the vicinity of the catalytic site. The findings presented here offer an insight as to how these acquired mutations could be associated to an improved adaptability of the virus for efficient human transmissibility.

Pyruvate : ferredoxin oxidoreductase (PFO) and iron only hydrogenase ([Fe]-HYD) are common enzymes among eukaryotic microbes that inhabit anaerobic niches. Their function is to maintain redox balance by donating electrons from food oxidation via ferredoxin (Fd) to protons, generating H 2 as a waste product. Operating in series, they constitute a soluble electron transport chain of one-electron transfers between FeS clusters. They fulfil the same function—redox balance—served by two electron-transfers in the NADH- and O 2 -dependent respiratory chains of mitochondria. Although they possess O 2 -sensitive FeS clusters, PFO, Fd and [Fe]-HYD are also present among numerous algae that produce O 2 . The evolutionary persistence of these enzymes among eukaryotic aerobes is traditionally explained as adaptation to facultative anaerobic growth. Here, we show that algae express enzymes of anaerobic energy metabolism at ambient O 2 levels (21% v/v), Chlamydomonas reinhardtii expresses them with diurnal regulation. High O 2 environments arose on Earth only approximately 450 million years ago. Gene presence/absence and gene expression data indicate that during the transition to high O 2 environments and terrestrialization, diverse algal lineages retained enzymes of Fd-dependent one-electron-based redox balance, while the land plant and land animal lineages underwent irreversible specialization to redox balance involving the O 2 -insensitive two-electron carrier NADH.

In prokaryotes, known mechanisms of lateral gene transfer (transformation, transduction, conjugation and gene transfer agents) generate new combinations of genes among chromosomes during evolution. In eukaryotes, whose host lineage is descended from archaea, lateral gene transfer from organelles to the nucleus occurs at endosymbiotic events. Recent genome analyses studying gene distributions have uncovered evidence for sporadic, discontinuous events of gene transfer from bacteria to archaea during evolution. Other studies have used traditional birth-and-death phylogenetic models to investigate prokaryote genome evolution to claim that gene transfer to archaea was continuous during evolution, rather than involving occasional periodic mass gene influx events. Here we test the ability of Count, a birth-and-death based program, to recover known events of mass acquisition and differential loss using plastid genomes and eukaryotic protein families that were acquired from plastids. Count showed a strong bias towards reconstructed histories having gene acquisitions distributed uniformly across the tree. Sometimes as many as nine different acquisitions by plastid DNA were inferred for the same protein family. That is, Count recovered gradual and continuous lateral gene transfer among lineages, even when massive gains followed by gradual differential loss is the true evolutionary process that generated the gene distribution data.

In prokaryotes, known mechanisms of lateral gene transfer (transformation, transduction, conjugation, and gene transfer agents) generate new combinations of genes among chromosomes during evolution. In eukaryotes, whose host lineage is descended from archaea, lateral gene transfer from organelles to the nucleus occurs at endosymbiotic events. Recent genome analyses studying gene distributions have uncovered evidence for sporadic, discontinuous events of gene transfer from bacteria to archaea during evolution. Other studies have used traditional models designed to investigate gene family size evolution (Count) to support claims that gene transfer to archaea was continuous during evolution, rather than involving occasional periodic mass gene influx events. Here, we show that the methodology used in analyses favoring continuous gene transfers to archaea was misapplied in other studies and does not recover known events of single simultaneous origin for many genes followed by differential loss in real data: plastid genomes. Using the same software and the same settings, we reanalyzed presence/absence pattern data for proteins encoded in plastid genomes and for eukaryotic protein families acquired from plastids. Contrary to expectations under a plastid origin model, we found that the methodology employed inferred that gene acquisitions occurred uniformly across the plant tree. Sometimes as many as nine different acquisitions by plastid DNA were inferred for the same protein family. That is, the methodology that recovered gradual and continuous lateral gene transfer among lineages for archaea obtains the same result for plastids, even though it is known that massive gains followed by gradual differential loss is the true evolutionary process that generated plastid gene distribution data. Our findings caution against the use of models designed to study gene family size evolution for investigating gene transfer processes, especially when transfers involving more than one gene per event are possible.

Chloroplasts arose from cyanobacteria, mitochondria arose from proteobacteria. Both organelles have conserved their prokaryotic biochemistry, but their genomes are reduced, and most organelle proteins are encoded in the nucleus. Endosymbiotic theory posits that bacterial genes in eukaryotic genomes entered the eukaryotic lineage via organelle ancestors. It predicts episodic influx of prokaryotic genes into the eukaryotic lineage, with acquisition corresponding to endosymbiotic events. Eukaryotic genome sequences, however, increasingly implicate lateral gene transfer, both from prokaryotes to eukaryotes and among eukaryotes, as a source of gene content variation in eukaryotic genomes, which predicts continuous, lineage-specific acquisition of prokaryotic genes in divergent eukaryotic groups. Here we discriminate between these two alternatives by clustering and phylogenetic analysis of eukaryotic gene families having prokaryotic homologues. Our results indicate (1) that gene transfer from bacteria to eukaryotes is episodic, as revealed by gene distributions, and coincides with major evolutionary transitions at the origin of chloroplasts and mitochondria; (2) that gene inheritance in eukaryotes is vertical, as revealed by extensive topological comparison, sparse gene distributions stemming from differential loss; and (3) that continuous, lineage-specific lateral gene transfer, although it sometimes occurs, does not contribute to long-term gene content evolution in eukaryotic genomes.

Endosymbiotic theory in eukaryotic-cell evolution rests upon a foundation of three cornerstone partners—the plastid (a cyanobacterium), the mitochondrion (a proteobacterium), and its host (an archaeon)—and carries a corollary that, over time, the majority of genes once present in the organelle genomes were relinquished to the chromosomes of the host (endosymbiotic gene transfer). However, notwithstanding eukaryote-specific gene inventions, single-gene phylogenies have never traced eukaryotic genes to three single prokaryotic sources, an issue that hinges crucially upon factors influencing phylogenetic inference. In the age of genomes, single-gene trees, once used to test the predictions of endosymbiotic theory, now spawn new theories that stand to eventually replace endosymbiotic theory with descriptive, gene tree-based variants featuring supernumerary symbionts: prokaryotic partners distinct from the cornerstone trio and whose existence is inferred solely from single-gene trees. We reason that the endosymbiotic ancestors of mitochondria and chloroplasts brought into the eukaryotic—and plant and algal—lineage a genome-sized sample of genes from the proteobacterial and cyanobacterial pangenomes of their respective day and that, even if molecular phylogeny were artifact-free, sampling prokaryotic pangenomes through endosymbiotic gene transfer would lead to inherited chimerism. Recombination in prokaryotes (transduction, conjugation, transformation) differs from recombination in eukaryotes (sex). Prokaryotic recombination leads to pangenomes, and eukaryotic recombination leads to vertical inheritance. Viewed from the perspective of endosymbiotic theory, the critical transition at the eukaryote origin that allowed escape from Muller’s ratchet—the origin of eukaryotic recombination, or sex—might have required surprisingly little evolutionary innovation.

The origin of plastids is best explained by endosymbiotic theory, which dates back to the early 1900s. Three lines of evidence based on protein import machineries and molecular phylogenies of eukaryote (host) and cyanobacterial (endosymbiont) genes point to a single origin of primary plastids, a unique and important event that successfully transferred two photosystems and oxygenic photosynthesis from prokaryotes to eukaryotes. The nature of the cyanobacterial lineage from which plastids originated has been a topic of investigation. Recent studies have focused on the branching position of the plastid lineage in the phylogeny based on cyanobacterial core genes, that is, genes shared by all cyanobacteria and plastids. These studies have delivered conflicting results, however. In addition, the core genes represent only a very small portion of cyanobacterial genomes and may not be a good proxy for the rest of the ancestral plastid genome. Information in plant nuclear genomes, where most genes that entered the eukaryotic lineage through acquisition from the plastid ancestor reside, suggests that heterocyst-forming cyanobacteria in Stanier’s sections IV and V are most similar to the plastid ancestor in terms of gene complement and sequence conservation, which is in agreement with models suggesting an important role of nitrogen fixation in symbioses involving cyanobacteria. Plastid origin is an ancient event that involved a prokaryotic symbiont and a eukaryotic host, organisms with different histories and genome evolutionary processes. The different modes of genome evolution in prokaryotes and eukaryotes bear upon our interpretations of plastid phylogeny.

Pittis and Gabaldón1 recently claimed that the mitochondrion came late in eukaryotic evolution, following an earlier phase of evolution in which the eukaryotic host lineage acquired genes from bacteria. Here we show that their paper has multiple fatal flaws founded in inappropriate statistical methods and analyses, in addition to erroneous interpretations.

Genomes record their own history. But if we want to look all the way back to life’s beginnings some 4 billion years ago, the record of microbial evolution that is preserved in prokaryotic genomes is not easy to read. Microbiology has a lot in common with geology in that regard. Geologists know that plate tectonics and erosion have erased much of the geological record, with ancient rocks being truly rare. The same is true of microbes. Lateral gene transfer (LGT) and sequence divergence have erased much of the evolutionary record that was once written in genomes, and it is not obvious which genes among sequenced genomes are genuinely ancient. Which genes trace to the last universal ancestor, LUCA? The classical approach has been to look for genes that are universally distributed. Another approach is to make all trees for all genes, and sift out the trees where signals have been overwritten by LGT. What is left ought to be ancient. If we do that, what do we find?

Life arose in a world without oxygen and the first organisms were anaerobes. Here we investigate the gene repertoire of the prokaryote common ancestor, estimating which genes it contained and to which lineages of modern prokaryotes it was most similar in terms of gene content. Using a phylogenetic approach we found that among trees for all 8779 protein families shared between 134 archaea and 1847 bacterial genomes, only 1045 have sequences from at least two bacterial and two archaeal groups and retain the ancestral archaeal–bacterial split. Among those, the genes shared by anaerobes were identified as candidate genes for the prokaryote common ancestor, which lived in anaerobic environments. We find that these anaerobic prokaryote common ancestor genes are today most frequently distributed among methanogens and clostridia, strict anaerobes that live from low free energy changes near the thermodynamic limit of life. The anaerobic families encompass genes for bifunctional acetyl-CoA-synthase/CO-dehydrogenase, heterodisulfide reductase subunits C and A, ferredoxins, and several subunits of the Mrp-antiporter/hydrogenase family, in addition to numerous S-adenosyl methionine (SAM) dependent methyltransferases. The data indicate a major role for methyl groups in the metabolism of the prokaryote common ancestor. The data furthermore indicate that the prokaryote ancestor possessed a rotor stator ATP synthase, but lacked cytochromes and quinones as well as identifiable redox-dependent ion pumping complexes. The prokaryote ancestor did possess, however, an Mrp-type H+/Na+ antiporter complex, capable of transducing geochemical pH gradients into biologically more stable Na+-gradients. The findings implicate a hydrothermal, autotrophic, and methyl-dependent origin of life. This article is part of a Special Issue entitled ‘EBEC 2016: 19th European Bioenergetics Conference, Riva del Garda, Italy, July 2–6, 2016’, edited by Prof. Paolo Bernardi., Highlights • Life arose without oxygen, the universal ancestor (Luca) was an anaerobe. • We used phylogenetic and physiological criteria to identify genes present in Luca. • An ancient core of 65 metabolic genes shed light on Luca's anaerobic lifestyle. • Ancient core genes are most widespread among modern methanogens and clostridia. • The data implicate a major role for methyl groups in Luca's anaerobic metabolism.

The mechanisms that underlie the origin of major prokaryotic groups are poorly understood. In principle, the origin of both species and higher taxa among prokaryotes should entail similar mechanisms--ecological interactions with the environment paired with natural genetic variation involving lineage-specific gene innovations and lineage-specific gene acquisitions. To investigate the origin of higher taxa in archaea, we have determined gene distributions and gene phylogenies for the 267,568 protein-coding genes of 134 sequenced archaeal genomes in the context of their homologues from 1,847 reference bacterial genomes. Archaeal-specific gene families define 13 traditionally recognized archaeal higher taxa in our sample. Here we report that the origins of these 13 groups unexpectedly correspond to 2,264 group-specific gene acquisitions from bacteria. Interdomain gene transfer is highly asymmetric, transfers from bacteria to archaea are more than fivefold more frequent than vice versa. Gene transfers identified at major evolutionary transitions among prokaryotes specifically implicate gene acquisitions for metabolic functions from bacteria as key innovations in the origin of higher archaeal taxa.

Enteric pathogens must overcome intestinal defenses to establish infection. In Drosophila, the ERK signaling pathway inhibits enteric virus infection. The intestinal microflora also impacts immunity but its role in enteric viral infection is unknown. Here we show that two signals are required to activate antiviral ERK signaling in the intestinal epithelium. One signal depends on recognition of peptidoglycan from the microbiota, particularly from the commensal Acetobacter pomorum, which primes the NF-kB-dependent induction of a secreted factor, Pvf2. However, the microbiota is not sufficient to induce this pathway; a second virus-initiated signaling involving release of transcriptional paused genes mediated by the kinase Cdk9 is also required for Pvf2 production. Pvf2 stimulates antiviral immunity by binding to the receptor tyrosine kinase PVR, which is necessary and sufficient for intestinal ERK responses. These findings demonstrate that sensing of specific commensals primes inflammatory signaling required for epithelial responses that restrict enteric viral infections.

Like biological species, languages change over time. As noted by Darwin, there are many parallels between language evolution and biological evolution. Insights into these parallels have also undergone change in the past 150 years. Just like genes, words change over time, and language evolution can be likened to genome evolution accordingly, but what kind of evolution? There are fundamental differences between eukaryotic and prokaryotic evolution. In the former, natural variation entails the gradual accumulation of minor mutations in alleles. In the latter, lateral gene transfer is an integral mechanism of natural variation. The study of language evolution using biological methods has attracted much interest of late, most approaches focusing on language tree construction. These approaches may underestimate the important role that borrowing plays in language evolution. Network approaches that were originally designed to study lateral gene transfer may provide more realistic insights into the complexities of language evolution.

Archaebacterial halophiles (Haloarchaea) are oxygen-respiring heterotrophs that derive from methanogens—strictly anaerobic, hydrogen-dependent autotrophs. Haloarchaeal genomes are known to have acquired, via lateral gene transfer (LGT), several genes from eubacteria, but it is yet unknown how many genes the Haloarchaea acquired in total and, more importantly, whether independent haloarchaeal lineages acquired their genes in parallel, or as a single acquisition at the origin of the group. Here we have studied 10 haloarchaeal and 1,143 reference genomes and have identified 1,089 haloarchaeal gene families that were acquired by a methanogenic recipient from eubacteria. The data suggest that these genes were acquired in the haloarchaeal common ancestor, not in parallel in independent haloarchaeal lineages, nor in the common ancestor of haloarchaeans and methanosarcinales. The 1,089 acquisitions include genes for catabolic carbon metabolism, membrane transporters, menaquinone biosynthesis, and complexes I–IV of the eubacterial respiratory chain that functions in the haloarchaeal membrane consisting of diphytanyl isoprene ether lipids. LGT on a massive scale transformed a strictly anaerobic, chemolithoautotrophic methanogen into the heterotrophic, oxygen-respiring, and bacteriorhodopsin-photosynthetic haloarchaeal common ancestor.

Language evolution is traditionally described in terms of family trees with ancestral languages splitting into descendent languages. However, it has long been recognized that language evolution also entails horizontal components, most commonly through lexical borrowing. For example, the English language was heavily influenced by Old Norse and Old French; eight per cent of its basic vocabulary is borrowed. Borrowing is a distinctly non-tree-like process—akin to horizontal gene transfer in genome evolution—that cannot be recovered by phylogenetic trees. Here, we infer the frequency of hidden borrowing among 2346 cognates (etymologically related words) of basic vocabulary distributed across 84 Indo-European languages. The dataset includes 124 (5%) known borrowings. Applying the uniformitarian principle to inventory dynamics in past and present basic vocabularies, we find that 1373 (61%) of the cognates have been affected by borrowing during their history. Our approach correctly identified 117 (94%) known borrowings. Reconstructed phylogenetic networks that capture both vertical and horizontal components of evolutionary history reveal that, on average, eight per cent of the words of basic vocabulary in each Indo-European language were involved in borrowing during evolution. Basic vocabulary is often assumed to be relatively resistant to borrowing. Our results indicate that the impact of borrowing is far more widespread than previously thought.

The concept of a last universal common ancestor of all cells (LUCA, or the progenote) is central to the study of early evolution and life's origin, yet information about how and where LUCA lived is lacking. We investigated all clusters and phylogenetic trees for 6.1 million protein coding genes from sequenced prokaryotic genomes in order to reconstruct the microbial ecology of LUCA. Among 286,514 protein clusters, we identified 355 protein families (∼0.1%) that trace to LUCA by phylogenetic criteria. Because these proteins are not universally distributed, they can shed light on LUCA's physiology. Their functions, properties and prosthetic groups depict LUCA as anaerobic, CO2-fixing, H2-dependent with a Wood-Ljungdahl pathway, N2-fixing and thermophilic. LUCA's biochemistry was replete with FeS clusters and radical reaction mechanisms. Its cofactors reveal dependence upon transition metals, flavins, S-adenosyl methionine, coenzyme A, ferredoxin, molybdopterin, corrins and selenium. Its genetic code required nucleoside modifications and S-adenosyl methionine-dependent methylations. The 355 phylogenies identify clostridia and methanogens, whose modern lifestyles resemble that of LUCA, as basal among their respective domains. LUCA inhabited a geochemically active environment rich in H2, CO2 and iron. The data support the theory of an autotrophic origin of life involving the Wood-Ljungdahl pathway in a hydrothermal setting.