Gene Transfer Among Viruses Substantially Contributes to Gene Gain of Giant Viruses.

The phylum Nucleocytoviricota comprises a diverse group of double-stranded DNA viruses that display a wide range of gene repertoires. Although these gene repertoires determine the characteristics of individual viruses, the evolutionary processes that have shaped the gene repertoires of extant viruses since their common ancestor are poorly characterized. In this study, we aimed to address this gap in knowledge by using amalgamated likelihood estimation, a probabilistic tree reconciliation method that infers evolutionary scenarios by distinguishing origination, gene duplications, virus-to-virus horizontal gene transfer (vHGT), and gene losses. We analyzed over 4,700 gene families from 195 genomes spanning all known viral orders. The evolutionary reconstruction suggests a history of extensive gene gains and losses during the evolution of these viruses, notably with vHGT contributing to gene gains at a comparable level to duplications and originations. The vHGT frequently occurred between phylogenetically closely related viruses, as well as between distantly related viruses with an overlapping host range. We observed a pattern of massive gene duplications that followed vHGTs for gene families that was potentially related to host range control and virus-host arms race. These results suggest that vHGT represents a previously overlooked, yet important, evolutionary force that integrates the evolutionary paths of multiple viruses and affects shaping of Nucleocytoviricota virus gene repertoires.

Duplication and sub-functionalization of flavonoid biosynthesis genes plays important role in Leguminosae root nodule symbiosis evolution.

Gene innovation plays an essential role in trait evolution. Rhizobial symbioses, the most important N2-fixing agent in agricultural systems that exists mainly in Leguminosae, is one of the most attractive evolution events. However, the gene innovations underlying Leguminosae root nodule symbiosis (RNS) remain largely unknown. Here, we investigated the gene gain event in Leguminosae RNS evolution through comprehensive phylogenomic analyses. We revealed that Leguminosae-gain genes were acquired by gene duplication and underwent a strong purifying selection. Kyoto Encyclopedia of Genes and Genomes analyses showed that the innovated genes were enriched in flavonoid biosynthesis pathways, particular downstream of chalcone synthase (CHS). Among them, Leguminosae-gain type Ⅱ chalcone isomerase (CHI) could be further divided into CHI1A and CHI1B clades, which resulted from the products of tandem duplication. Furthermore, the duplicated CHI genes exhibited exon-intron structural divergences evolved through exon/intron gain/loss and insertion/deletion. Knocking down CHI1B significantly reduced nodulation in Glycine max (soybean) and Medicago truncatula; whereas, knocking down its duplication gene CHI1A had no effect on nodulation. Therefore, Leguminosae-gain type Ⅱ CHI participated in RNS and the duplicated CHI1A and CHI1B genes exhibited RNS functional divergence. This study provides functional insights into Leguminosae-gain genetic innovation and sub-functionalization after gene duplication that contribute to the evolution and adaptation of RNS in Leguminosae.

ZNF217 Gene Copy Number as a Marker of Response to Standard Therapy Drugs According to ERα Status in Breast Cancer.

Purpose: The therapeutic decision for the management of breast cancer (BC) patients is based on the evaluation of prognostic factors alongside clinical and pathological parameters. Despite the use of standard biomarkers, response and resistance to therapy represent a challenge for clinicians. Among the new potential biomarkers for BC the ZNF217 gene have gained importance in recent years. However, while associations between ZNF217 gene copy number and clinicopathological characteristics have been established, its correlation with treatment response remains unclear. Patients and Methods: This study aimed to evaluate the ZNF217 gene copy number and establish its associations with treatment response in estrogen receptor positive (ERα+) and ERα negative (ERα-) BC cell lines. In addition, a validation of the relationship between ZNF217 gene copy number and its prognostic value was performed using datasets of BC patients retrieved from the cBioPortal public database. Results: Our data show that in ERα+ cells, ZNF217 gene copy number increase (amplification), while cell proliferation decreases in response to standard drug treatments. In contrast, both ZNF217 gene copy number (gain) and cell proliferation increases in response to standard drug treatments in ERα- cells. The results obtained align with findings from the cBioPortal database analysis, demonstrating that ERα+/HER2- low proliferation patients, exhibiting ZNF217 gene amplification or gain, have a significantly higher survival probability after treatment, compared to ERα-/HER2- and HER2+ patients. Conclusion: Our results suggest that in ERα+ BC cells, ZNF217 gene amplification could be indicative of a favorable response, while in ERα- BC cells, ZNF217 gene gain could be postulated as a potential predictor of treatment resistance. A broader understanding of the role of ZNF217 gene in treatment response, together with prospective studies in BC patients, could contribute to confirming our data, as well as optimizing existing treatments and exploring novel approaches to improve overall cancer treatment outcomes.

Insights into early evolutionary adaptations of the Akkermansia genus to the vertebrate gut.

Akkermansia, a relevant mucin degrader from the vertebrate gut microbiota, is a member of the deeply branched Verrucomicrobiota, as well as the only known member of this phylum to be described as inhabitants of the gut. Only a few Akkermansia species have been officially described so far, although there is genomic evidence addressing the existence of more species-level variants for this genus. This niche specialization makes Akkermansia an interesting model for studying the evolution of microorganisms to their adaptation to the gastrointestinal tract environment, including which kind of functions were gained when the Akkermansia genus originated or how the evolutionary pressure functions over those genes. In order to gain more insight into Akkermansia adaptations to the gastrointestinal tract niche, we performed a phylogenomic analysis of 367 high-quality Akkermansia isolates and metagenome-assembled genomes, in addition to other members of Verrucomicrobiota. This work was focused on three aspects: the definition of Akkermansia genomic species clusters and the calculation and functional characterization of the pangenome for the most represented species; the evolutionary relationship between Akkermansia and their closest relatives from Verrucomicrobiota, defining the gene families which were gained or lost during the emergence of the last Akkermansia common ancestor (LAkkCA) and; the evaluation of the evolutionary pressure metrics for each relevant gene family of main Akkermansia species. This analysis found 25 Akkermansia genomic species clusters distributed in two main clades, divergent from their non-Akkermansia relatives. Pangenome analyses suggest that Akkermansia species have open pangenomes, and the gene gain/loss model indicates that genes associated with mucin degradation (both glycoside hydrolases and peptidases), (micro)aerobic metabolism, surface interaction, and adhesion were part of LAkkCA. Specifically, mucin degradation is a very ancestral innovation involved in the origin of Akkermansia. Horizontal gene transfer detection suggests that Akkermansia could receive genes mostly from unknown sources or from other Gram-negative gut bacteria. Evolutionary metrics suggest that Akkemansia species evolved differently, and even some conserved genes suffered different evolutionary pressures among clades. These results suggest a complex evolutionary landscape of the genus and indicate that mucin degradation could be an essential feature in Akkermansia evolution as a symbiotic species.

Chance Favors the Prepared Genomes: Horizontal Transfer Shapes the Emergence of Antibiotic Resistance Mutations in Core Genes.

Bacterial lineages acquire novel traits at diverse rates in part because the genetic background impacts the successful acquisition of novel genes by horizontal transfer. Yet, how horizontal transfer affects the subsequent evolution of core genes remains poorly understood. Here, we studied the evolution of resistance to quinolones in Escherichia coli accounting for population structure. We found 60 groups of genes whose gain or loss induced an increase in the probability of subsequently becoming resistant to quinolones by point mutations in the gyrase and topoisomerase genes. These groups include functions known to be associated with direct mitigation of the effect of quinolones, with metal uptake, cell growth inhibition, biofilm formation, and sugar metabolism. Many of them are encoded in phages or plasmids. Although some of the chronologies may reflect epidemiological trends, many of these groups encoded functions providing latent phenotypes of antibiotic low-level resistance, tolerance, or persistence under quinolone treatment. The mutations providing resistance were frequent and accumulated very quickly. Their emergence was found to increase the rate of acquisition of other antibiotic resistances setting the path for multidrug resistance. Hence, our findings show that horizontal gene transfer shapes the subsequent emergence of adaptive mutations in core genes. In turn, these mutations further affect the subsequent evolution of resistance by horizontal gene transfer. Given the substantial gene flow within bacterial genomes, interactions between horizontal transfer and point mutations in core genes may be a key to the success of adaptation processes.

Phylogenomic analysis of the Porphyromonas gingivalis - Porphyromonas gulae duo: approaches to the origin of periodontitis.

Porphyromonas gingivalis is an oral human pathogen associated with the onset and progression of periodontitis, a chronic immune-inflammatory disease characterized by the destruction of the teeth-supporting tissue. P. gingivalis belongs to the genus Porphyromonas, which is characterized by being composed of Gram-negative, asaccharolytic, non-spore-forming, non-motile, obligatory anaerobic species, inhabiting niches such as the oral cavity, urogenital tract, gastrointestinal tract and infected wound from different mammals including humans. Among the Porphyromonas genus, P. gingivalis stands out for its specificity in colonizing the human oral cavity and its keystone pathogen role in periodontitis pathogenesis. To understand the evolutionary process behind P. gingivalis in the context of the Pophyoromonas genus, in this study, we performed a comparative genomics study with publicly available Porphyromonas genomes, focused on four main objectives: (A) to confirm the phylogenetic position of P. gingivalis in the Porphyromonas genus by phylogenomic analysis; (B) the definition and comparison of the pangenomes of P. gingivalis and its relative P. gulae; and (C) the evaluation of the gene family gain/loss events during the divergence of P. gingivalis and P. gulae; (D) the evaluation of the evolutionary pressure (represented by the calculation of Tajima-D values and dN/dS ratios) comparing gene families of P. gingivalis and P. gulae. Our analysis found 84 high-quality assemblies representing P. gingivalis and 14 P. gulae strains (from a total of 233 Porphyromonas genomes). Phylogenomic analysis confirmed that P. gingivalis and P. gulae are highly related lineages, close to P. loveana. Both organisms harbored open pangenomes, with a strong core-to-accessory ratio for housekeeping genes and a negative ratio for unknown function genes. Our analyses also characterized the gene set differentiating P. gulae from P. gingivalis, mainly associated with unknown functions. Relevant virulence factors, such as the FimA, Mfa1, and the hemagglutinins, are conserved in P. gulae, P. gingivalis, and P. loveana, suggesting that the origin of those factors occurred previous to the P. gulae - P. gingivalis divergence. These results suggest an unexpected evolutionary relationship between the P. gulae - P. gingivalis duo and P. loveana, showing more clues about the origin of the role of those organisms in periodontitis.

Analysis of multipartite bacterial genomes using alignment free and alignment-based pipelines.

Since the discovery of second chromosome in Rhodobacter sphaeroides 2.4.1 in 1989, multipartite genomes have been reported in over three hundred bacterial species under nine different phyla. This has shattered the unipartite (single chromosome) genome dogma in bacteria. Since then, many questions on various aspects of multipartite genomes in bacteria have been addressed. However, our understanding of how multipartite genomes emerge and evolve is still lacking. Importantly, the knowledge of genetic factors underlying the differences in multipartite and single-chromosome genomes is lacking. In this work, we have performed comparative evolutionary and functional genomics analyses to identify molecular factors that discriminate multipartite from unipartite bacteria, with the goal to decipher taxon-specific factors, and those that are prevalent across the taxa, underlying these traits. We assessed the roles of evolutionary mechanisms, specifically gene gain, in driving the divergence of bacteria with single and multiple chromosomes. In addition, we performed functional genomic analysis to garner support for our findings from comparative evolutionary analysis. We found genes such as those encoding conserved hypothetical proteins in Deinococcus radiodurans R1, and putative phage phi-C31 gp36 major capsid like and hypothetical proteins in Rhodobacter sphaeroides 2.4.1, which are located on accessory chromosomes in these bacteria but were not found in the inferred ancestral sequences, and on the primary chromosomes, as well as were not found in their closest relatives with single chromosome within the same clade. Our study shines a new light on the potential roles of the secondary chromosomes in helping bacteria with multipartite genomes to adapt to specialized environments or growth conditions.

Evolutionary history of metazoan TMEM16 family.

Most of Transmembrane protein 16 (TMEM16) family members function as either a Ca2+-activated Cl- channel (CaCC) or phospholipid scramblase (CaPLSase) and play diverse physiological roles. It is well conserved in eukaryotes; however, the origin and evolution of different subfamilies in Metazoa are not yet understood. To uncover the evolutionary history of the TMEM16 family, we analyzed 398 proteins from 74 invertebrate species using evolutionary genomics. We found that the TMEM16C-F and J subfamilies are vertebrate-specific, but the TMEM16A/B, G, H, and K subfamilies are ancient and present in many, but not all metazoan species. The most ancient subfamilies in Metazoa, TMEM16L and M, are only maintained in limited species. TMEM16N and O are Cnidaria- and Ecdysozoa-specific subfamilies, respectively, and Ctenophora, Xenacoelomorpha, and Rotifera contain species-specific proteins. We also identified TMEM16 genes that are closely linked together in the genome, suggesting that they have been generated via recent gene duplication. The anoctamin domain structures of invertebrate-specific TMEM16 proteins predicted by AlphaFold2 contain conserved Ca2+-binding motifs and permeation pathways with either narrow or wide inner gates. The inner gate distance of TMEM16 protein may have frequently switched during metazoan evolution, and thus determined the function of the protein as either CaCC or CaPLSase. These results demonstrate that TMEM16 family has evolved by gene gain and loss in metazoans, and the genes have been generally under purifying selection to maintain protein structures and physiological functions.

Reconstruction of gene innovation associated with major evolutionary transitions in the kingdom Fungi.

BACKGROUND: Fungi exhibit astonishing diversity with multiple major phenotypic transitions over the kingdom's evolutionary history. As part of this process, fungi developed hyphae, adapted to land environments (terrestrialization), and innovated their sexual structures. These changes also helped fungi establish ecological relationships with other organisms (animals and plants), but the genomic basis of these changes remains largely unknown. RESULTS: By systematically analyzing 304 genomes from all major fungal groups, together with a broad range of eukaryotic outgroups, we have identified 188 novel orthogroups associated with major changes during the evolution of fungi. Functional annotations suggest that many of these orthogroups were involved in the formation of key trait innovations in extant fungi and are functionally connected. These innovations include components for cell wall formation, functioning of the spindle pole body, polarisome formation, hyphal growth, and mating group signaling. Innovation of mitochondria-localized proteins occurred widely during fungal transitions, indicating their previously unrecognized importance. We also find that prokaryote-derived horizontal gene transfer provided a small source of evolutionary novelty with such genes involved in key metabolic pathways. CONCLUSIONS: The overall picture is one of a relatively small number of novel genes appearing at major evolutionary transitions in the phylogeny of fungi, with most arising de novo and horizontal gene transfer providing only a small additional source of evolutionary novelty. Our findings contribute to an increasingly detailed portrait of the gene families that define fungal phyla and underpin core features of extant fungi.

Convergence Analysis of Rust Fungi and Anther Smuts Reveals Their Common Molecular Adaptation to a Phytoparasitic Lifestyle.

Convergent evolution between distantly related taxa often mirrors adaptation to similar environments. Rust fungi and anther smuts, which belong to different classes in Pucciniomycotina, have independently evolved a phytoparasitic lifestyle, representing an example of convergent evolution in the fungal kingdom. To investigate their adaptations and the genetic bases underlying their phytoparasitic lifestyles, we performed genome-wide convergence analysis of amino acid substitutions, evolutionary rates, and gene gains and losses. Convergent substitutions were detected in ATPeV0D and RP-S27Ae, two genes important for the generation of turgor pressure and ribosomal biosynthesis, respectively. A total of 51 positively selected genes were identified, including eight genes associated with translation and three genes related to the secretion pathway. In addition, rust fungi and anther smuts contained more proteins associated with oligopeptide transporters and vacuolar proteases than did other fungi. For rust fungi and anther smuts, these forms of convergence suggest four adaptive mechanisms for a phytoparasitic lifestyle: 1) reducing the metabolic demand for hyphal growth and penetration at the pre-penetration stage, 2) maintaining the efficiency of protein synthesis during colonization, 3) ensuring the normal secretion of rapidly evolving secreted proteins, and 4) improving the capacity for oligopeptide metabolism. Our results are the first to shed light on the genetic convergence mechanisms and molecular adaptation underlying phytoparasitic lifestyles in fungi.

A Large-Scale Genome-Based Survey of Acidophilic Bacteria Suggests That Genome Streamlining Is an Adaption for Life at Low pH.

The genome streamlining theory suggests that reduction of microbial genome size optimizes energy utilization in stressful environments. Although this hypothesis has been explored in several cases of low-nutrient (oligotrophic) and high-temperature environments, little work has been carried out on microorganisms from low-pH environments, and what has been reported is inconclusive. In this study, we performed a large-scale comparative genomics investigation of more than 260 bacterial high-quality genome sequences of acidophiles, together with genomes of their closest phylogenetic relatives that live at circum-neutral pH. A statistically supported correlation is reported between reduction of genome size and decreasing pH that we demonstrate is due to gene loss and reduced gene sizes. This trend is independent from other genome size constraints such as temperature and G + C content. Genome streamlining in the evolution of acidophilic bacteria is thus supported by our results. The analyses of predicted Clusters of Orthologous Genes (COG) categories and subcellular location predictions indicate that acidophiles have a lower representation of genes encoding extracellular proteins, signal transduction mechanisms, and proteins with unknown function but are enriched in inner membrane proteins, chaperones, basic metabolism, and core cellular functions. Contrary to other reports for genome streamlining, there was no significant change in paralog frequencies across pH. However, a detailed analysis of COG categories revealed a higher proportion of genes in acidophiles in the following categories: "replication and repair," "amino acid transport," and "intracellular trafficking". This study brings increasing clarity regarding the genomic adaptations of acidophiles to life at low pH while putting elements, such as the reduction of average gene size, under the spotlight of streamlining theory.

Novel patterns of expression and recruitment of new genes on the t-haplotype, a mouse selfish chromosome.

The t-haplotype of mice is a classical model for autosomal transmission distortion. A largely non-recombining variant of the proximal region of chromosome 17, it is transmitted to more than 90% of the progeny of heterozygous males through the disabling of sperm carrying a standard chromosome. While extensive genetic and functional work has shed light on individual genes involved in drive, much less is known about the evolution and function of the rest of its hundreds of genes. Here, we characterize the sequence and expression of dozens of t-specific transcripts and of their chromosome 17 homologues. Many genes showed reduced expression of the t-allele, but an equal number of genes showed increased expression of their t-copy, consistent with increased activity or a newly evolved function. Genes on the t-haplotype had a significantly higher non-synonymous substitution rate than their homologues on the standard chromosome, with several genes harbouring dN/dS ratios above 1. Finally, the t-haplotype has acquired at least two genes from other chromosomes, which show high and tissue-specific expression. These results provide a first overview of the gene content of this selfish element, and support a more dynamic evolutionary scenario than expected of a large genomic region with almost no recombination.

Comparative genome analysis of plant ascomycete fungal pathogens with different lifestyles reveals distinctive virulence strategies.

BACKGROUND: Pathogens have evolved diverse lifestyles and adopted pivotal new roles in both natural ecosystems and human environments. However, the molecular mechanisms underlying their adaptation to new lifestyles are obscure. Comparative genomics was adopted to determine distinct strategies of plant ascomycete fungal pathogens with different lifestyles and to elucidate their distinctive virulence strategies. RESULTS: We found that plant ascomycete biotrophs exhibited lower gene gain and loss events and loss of CAZyme-encoding genes involved in plant cell wall degradation and biosynthesis gene clusters for the production of secondary metabolites in the genome. Comparison with the candidate effectome detected distinctive variations between plant biotrophic pathogens and other groups (including human, necrotrophic and hemibiotrophic pathogens). The results revealed the biotroph-specific and lifestyle-conserved candidate effector families. These data have been configured in web-based genome browser applications for public display ( http://lab.malab.cn/soft/PFPG ). This resource allows researchers to profile the genome, proteome, secretome and effectome of plant fungal pathogens. CONCLUSIONS: Our findings demonstrated different genome evolution strategies of plant fungal pathogens with different lifestyles and explored their lifestyle-conserved and specific candidate effectors. It will provide a new basis for discovering the novel effectors and their pathogenic mechanisms.

Ancient Gene Capture and Recent Gene Loss Shape the Evolution of Orthopoxvirus-Host Interaction Genes.

The survival of viruses depends on their ability to resist host defenses and, of all animal virus families, the poxviruses have the most antidefense genes. Orthopoxviruses (ORPV), a genus within the subfamily Chordopoxvirinae, infect diverse mammals and include one of the most devastating human pathogens, the now eradicated smallpox virus. ORPV encode ∼200 genes, of which roughly half are directly involved in virus genome replication and expression as well as virion morphogenesis. The remaining ∼100 "accessory" genes are responsible for virus-host interactions, particularly counter-defense of innate immunity. Complete sequences are currently available for several hundred ORPV genomes isolated from a variety of mammalian hosts, providing a rich resource for comparative genomics and reconstruction of ORPV evolution. To identify the provenance and evolutionary trends of the ORPV accessory genes, we constructed clusters including the orthologs of these genes from all chordopoxviruses. Most of the accessory genes were captured in three major waves early in chordopoxvirus evolution, prior to the divergence of ORPV and the sister genus Centapoxvirus from their common ancestor. The capture of these genes from the host was followed by extensive gene duplication, yielding several paralogous gene families. In addition, nine genes were gained during the evolution of ORPV themselves. In contrast, nearly every accessory gene was lost, some on multiple, independent occasions in numerous lineages of ORPV, so that no ORPV retains them all. A variety of functional interactions could be inferred from examination of pairs of ORPV accessory genes that were either often or rarely lost concurrently. IMPORTANCE Orthopoxviruses (ORPV) include smallpox (variola) virus, one of the most devastating human pathogens, and vaccinia virus, comprising the vaccine used for smallpox eradication. Among roughly 200 ORPV genes, about half are essential for genome replication and expression as well as virion morphogenesis, whereas the remaining half consists of accessory genes counteracting the host immune response. We reannotated the accessory genes of ORPV, predicting the functions of uncharacterized genes, and reconstructed the history of their gain and loss during the evolution of ORPV. Most of the accessory genes were acquired in three major waves antedating the origin of ORPV from chordopoxviruses. The evolution of ORPV themselves was dominated by gene loss, with numerous genes lost at the base of each major group of ORPV. Examination of pairs of ORPV accessory genes that were either often or rarely lost concurrently during ORPV evolution allows prediction of different types of functional interactions.

Insights on the Evolutionary Genomics of the Blautia Genus: Potential New Species and Genetic Content Among Lineages.

Blautia, a genus established in 2008, is a relevantly abundant taxonomic group present in the microbiome of human and other mammalian gastrointestinal (GI) tracts. Several described (or proposed) Blautia species are available at this date. However, despite the increasing level of knowledge about Blautia, its diversity is still poorly understood. The increasing availability of Blautia genomic sequences in the public databases opens the possibility to study this genus from a genomic perspective. Here we report the pangenome analysis and the phylogenomic study of 225 Blautia genomes available in RefSeq. We found 33 different potential species at the genomic level, 17 of them previously undescribed; we also confirmed by genomic standards the status of 4 previously proposed new Blautia species. Comparative genomic analyses suggest that the Blautia pangenome is open, with a relatively small core genome (∼ 700-800 gene families). Utilizing a set of representative genomes, we performed a gene family gain/loss model for the genus, showing that despite terminal nodes suffered more massive gene gain events than internal nodes (i.e., predicted ancestors), some ancestors were predicted to have gained an important number of gene families, some of them associated with the possible acquisition of metabolic abilities. Gene loss events remained lower than gain events in most cases. General aspects regarding pangenome composition and gene gain/loss events are discussed, as well as the proposition of changes in the taxonomic assignment of B. coccoides TY and the proposition of a new species, "B. pseudococcoides.".

A Novel Auxiliary Agarolytic Pathway Expands Metabolic Versatility in the Agar-Degrading Marine Bacterium Colwellia echini A3T.

Marine microorganisms encode a complex repertoire of carbohydrate-active enzymes (CAZymes) for the catabolism of algal cell wall polysaccharides. While the core enzyme cascade for degrading agar is conserved across agarolytic marine bacteria, gain of novel metabolic functions can lead to the evolutionary expansion of the gene repertoire. Here, we describe how two less-abundant GH96 α-agarases harbored in the agar-specific polysaccharide utilization locus (PUL) of Colwellia echini strain A3T facilitate the versatility of the agarolytic pathway. The cellular and molecular functions of the α-agarases examined by genomic, transcriptomic, and biochemical analyses revealed that α-agarases of C. echini A3T create a novel auxiliary pathway. α-Agarases convert even-numbered neoagarooligosaccharides to odd-numbered agaro- and neoagarooligosaccharides, providing an alternative route for the depolymerization process in the agarolytic pathway. Comparative genomic analysis of agarolytic bacteria implied that the agarolytic gene repertoire in marine bacteria has been diversified during evolution, while the essential core agarolytic gene set has been conserved. The expansion of the agarolytic gene repertoire and novel hydrolytic functions, including the elucidated molecular functionality of α-agarase, promote metabolic versatility by channeling agar metabolism through different routes. IMPORTANCEColwellia echini A3T is an example of how the gain of gene(s) can lead to the evolutionary expansion of agar-specific polysaccharide utilization loci (PUL). C. echini A3T encodes two α-agarases in addition to the core ß-agarolytic enzymes in its agarolytic PUL. Among the agar-degrading CAZymes identified so far, only a few α-agarases have been biochemically characterized. The molecular and biological functions of two α-agarases revealed that their unique hydrolytic pattern leads to the emergence of auxiliary agarolytic pathways. Through the combination of transcriptomic, genomic, and biochemical evidence, we elucidate the complete α-agarolytic pathway in C. echini A3T. The addition of α-agarases to the agarolytic enzyme repertoire might allow marine agarolytic bacteria to increase competitive abilities through metabolic versatility.

Genome Analysis of Endotrypanum and Porcisia spp., Closest Phylogenetic Relatives of Leishmania, Highlights the Role of Amastins in Shaping Pathogenicity.

While numerous genomes of Leishmania spp. have been sequenced and analyzed, an understanding of the evolutionary history of these organisms remains limited due to the unavailability of the sequence data for their closest known relatives, Endotrypanum and Porcisia spp., infecting sloths and porcupines. We have sequenced and analyzed genomes of three members of this clade in order to fill this gap. Their comparative analyses revealed only minute differences from Leishmaniamajor genome in terms of metabolic capacities. We also documented that the number of genes under positive selection on the Endotrypanum/Porcisia branch is rather small, with the flagellum-related group of genes being over-represented. Most significantly, the analysis of gene family evolution revealed a substantially reduced repertoire of surface proteins, such as amastins and biopterin transporters BT1 in the Endotrypanum/Porcisia species when compared to amastigote-dwelling Leishmania. This reduction was especially pronounced for δ-amastins, a subfamily of cell surface proteins crucial in the propagation of Leishmania amastigotes inside vertebrate macrophages and, apparently, dispensable for Endotrypanum/Porcisia, which do not infect such cells.

Genome-wide identification and characterization of olfactory receptor genes in common carp (Cyprinus carpio).

The environment contains a large extent of chemical information, which could be detected as olfactory sense. Olfactory in vertebrates plays important roles on many aspects during life time, including localizing prey or food, avoiding predators, mating behavior and social communication. Considering the essential role of olfactory receptors in the specific recognition of diverse stimuli, understanding the evolutionary dynamics of olfactory receptors in teleost means a lot, especially in the allotetraploid common carp, who has undergone the fourth whole-genome duplication event. Here, we identified the whole set of olfactory receptor genes in representative teleosts and found a significant contraction in common carp when compared with other teleosts. Odorant receptor genes (OR) occupy the most among four groups of olfactory receptors, including 33 functional genes and 16 pseudogenes. Furthermore, 6 trace amine-associated receptor (TAAR) genes (including 1 pseudogene), 7 odorant-related-A receptor genes, and 10 olfactory C family receptor genes (including 3 pseudogenes) were identified in common carp. Phylogenetic and motif analysis were performed to illustrate the phylogenetic relationship and structural conservation of teleost olfactory receptors. Selection pressure analysis suggested that olfactory receptor groups in common carp were all under relaxed purifying-selection. Additionally, gene expression divergences for olfactory receptor genes were investigated during embryonic development stages of common carp. We aim to determine the abundance of common carp olfactory receptor genes, explore the evolutionary fate and expression dynamics, and provide some genomic clues for the evolution of polyploid olfactory after whole-genome duplication and for future studies of teleost olfactory.

Evolution of Microbial Genomics: Conceptual Shifts over a Quarter Century.

Prokaryote genomics started in earnest in 1995, with the complete sequences of two small bacterial genomes, those of Haemophilus influenzae and Mycoplasma genitalium. During the next quarter century, the prokaryote genome database has been growing exponentially, with no saturation in sight. For most of these 25 years, genome sequencing remained limited to cultivable microbes. Together with next-generation sequencing methods, advances in metagenomics and single-cell genomics have lifted this limitation, providing for an increasingly unbiased characterization of the global prokaryote diversity. Advances in computational genomics followed the progress of genome sequencing, even if occasionally lagging behind. Several major new branches of bacteria and archaea were discovered, including Asgard archaea, the apparent closest relatives of eukaryotes and expansive groups of bacteria and archaea with small genomes thought to be symbionts of other prokaryotes. Comparative analysis of numerous prokaryote genomes spanning a wide range of evolutionary distances changed the conceptual foundations of microbiology, supplanting the notion of species genomes with fixed gene sets with that of dynamic pangenomes and the notion of a single Tree of Life (ToL) with a statistical tree-like trend among individual gene trees. Strides were also made towards a theory and quantitative laws of prokaryote genome evolution.

Evolutionary Genomics of Structural Variation in Asian Rice (Oryza sativa) Domestication.

Structural variants (SVs) are a largely unstudied feature of plant genome evolution, despite the fact that SVs contribute substantially to phenotypes. In this study, we discovered SVs across a population sample of 347 high-coverage, resequenced genomes of Asian rice (Oryza sativa) and its wild ancestor (O. rufipogon). In addition to this short-read data set, we also inferred SVs from whole-genome assemblies and long-read data. Comparisons among data sets revealed different features of genome variability. For example, genome alignment identified a large (∼4.3 Mb) inversion in indica rice varieties relative to japonica varieties, and long-read analyses suggest that ∼9% of genes from the outgroup (O. longistaminata) are hemizygous. We focused, however, on the resequencing sample to investigate the population genomics of SVs. Clustering analyses with SVs recapitulated the rice cultivar groups that were also inferred from SNPs. However, the site-frequency spectrum of each SV type-which included inversions, duplications, deletions, translocations, and mobile element insertions-was skewed toward lower frequency variants than synonymous SNPs, suggesting that SVs may be predominantly deleterious. Among transposable elements, SINE and mariner insertions were found at especially low frequency. We also used SVs to study domestication by contrasting between rice and O. rufipogon. Cultivated genomes contained ∼25% more derived SVs and mobile element insertions than O. rufipogon, indicating that SVs contribute to the cost of domestication in rice. Peaks of SV divergence were enriched for known domestication genes, but we also detected hundreds of genes gained and lost during domestication, some of which were enriched for traits of agronomic interest.