Conserved Gene Cluster Research Articles

Gene clusters for biosynthesis of mycosporine-like amino acids in dinoflagellate nuclear genomes: Possible recent horizontal gene transfer between species of Symbiodiniaceae (Dinophyceae).

Global warming increases the temperature of the ocean surface, which can disrupt dinoflagellate‐coral symbioses and result in coral bleaching. Photosynthetic dinoflagellates of the family Symbiodiniaceae include bleaching‐tolerant and bleaching‐sensitive coral symbionts. Therefore, understanding the molecular mechanisms for changing symbiont diversity is potentially useful to assist recovery of coral holobionts (corals and their associated microbes, including multiple species of Symbiodiniaceae), although sexual reproduction has not been observed in the Symbiodiniaceae. Recent molecular phylogenetic analyses estimate that the Symbiodiniaceae appeared 160 million years ago and diversified into 15 groups, five genera of which now have available draft genomes (i.e., Symbiodinium, Durusdinium, Breviolum, Fugacium, and Cladocopium). Comparative genomic analyses have suggested that crown groups have fewer gene families than early‐diverging groups, although many genes that were probably acquired via gene duplications and horizontal gene transfers (HGTs) have been found in each decoded genome. Because UV stress is likely a contributor to coral bleaching, and because the highly conserved gene cluster for mycosporine‐like amino acid (MAA) biosynthesis has been found in thermal‐tolerant symbiont genomes, I reviewed genomic features of the Symbiodiniaceae, focusing on possible acquisition of a biosynthetic gene cluster for MAAs, which absorb UV radiation. On the basis of highly conserved noncoding sequences, I hypothesized that HGTs have occurred among members of the Symbiodiniaceae and have contributed to the diversification of Symbiodiniaceae–host relationships. Finally, I proposed that bleaching tolerance may be strengthened by multiple MAAs from both symbiotic dinoflagellates and corals.

Fucosylated Human Milk Oligosaccharide Foraging within the Species Bifidobacterium pseudocatenulatum Is Driven by Glycosyl Hydrolase Content and Specificity.

Human milk enriches members of the genus Bifidobacterium in the infant gut. One species, Bifidobacterium pseudocatenulatum, is found in the gastrointestinal tracts of adults and breastfed infants. In this study, B. pseudocatenulatum strains were isolated and characterized to identify genetic adaptations to the breastfed infant gut. During growth on pooled human milk oligosaccharides (HMOs), we observed two distinct groups of B. pseudocatenulatum, isolates that readily consumed HMOs and those that did not, a difference driven by variable catabolism of fucosylated HMOs. A conserved gene cluster for fucosylated HMO utilization was identified in several sequenced B. pseudocatenulatum strains. One isolate, B. pseudocatenulatum MP80, which uniquely possessed GH95 and GH29 α-fucosidases, consumed the majority of fucosylated HMOs tested. Furthermore, B. pseudocatenulatum SC585, which possesses only a single GH95 α-fucosidase, lacked the ability to consume the complete repertoire of linkages within the fucosylated HMO pool. Analysis of the purified GH29 and GH95 fucosidase activities directly on HMOs revealed complementing enzyme specificities with the GH95 enzyme preferring 1-2 fucosyl linkages and the GH29 enzyme favoring 1-3 and 1-4 linkages. The HMO-binding specificities of the family 1 solute-binding protein component linked to the fucosylated HMO gene cluster in both SC585 and MP80 are similar, suggesting differential transport of fucosylated HMO is not a driving factor in each strain's distinct HMO consumption pattern. Taken together, these data indicate the presence or absence of specific α-fucosidases directs the strain-specific fucosylated HMO utilization pattern among bifidobacteria and likely influences competitive behavior for HMO foraging in situ. IMPORTANCE Often isolated from the human gut, microbes from the bacterial family Bifidobacteriaceae commonly possess genes enabling carbohydrate utilization. Isolates from breastfed infants often grow on and possess genes for the catabolism of human milk oligosaccharides (HMOs), glycans found in human breast milk. However, catabolism of structurally diverse HMOs differs between bifidobacterial strains. This study identifies key gene differences between Bifidobacterium pseudocatenulatum isolates that may impact whether a microbe successfully colonizes an infant gut. In this case, the presence of complementary α-fucosidases may provide an advantage to microbes seeking residence in the infant gut. Such knowledge furthers our understanding of how diet drives bacterial colonization of the infant gut.

Density-based binning of gene clusters to infer function or evolutionary history using GeneGrouper.

Identifying variant forms of gene clusters of interest in phylogenetically proximate and distant taxa can help to infer their evolutionary histories and functions. Conserved gene clusters may differ by only a few genes, but these small differences can in turn induce substantial phenotypes, such as by the formation of pseudogenes or insertions interrupting regulation. Particularly as microbial genomes and metagenomic assemblies become increasingly abundant, unsupervised grouping of similar, but not necessarily identical, gene clusters into consistent bins can provide a population-level understanding of their gene content variation and functional homology. We developed GeneGrouper, a command-line tool that uses a density-based clustering method to group gene clusters into bins. GeneGrouper demonstrated high recall and precision in benchmarks for the detection of the 23-gene Salmonella enterica LT2 Pdu gene cluster and four-gene Pseudomonas aeruginosa PAO1 Mex gene cluster among 435 genomes spanning mixed taxa. In a subsequent application investigating the diversity and impact of gene-complete and -incomplete LT2 Pdu gene clusters in 1130 S.enterica genomes, GeneGrouper identified a novel, frequently occurring pduN pseudogene. When investigated in vivo, introduction of the pduN pseudogene negatively impacted microcompartment formation. We next demonstrated the versatility of GeneGrouper by clustering distant homologous gene clusters and variable gene clusters found in integrative and conjugative elements. GeneGrouper software and code are publicly available at https://pypi.org/project/GeneGrouper/. Supplementary data are available at Bioinformatics online.

Approximate search for known gene clusters in new genomes using PQ-trees

Gene clusters are groups of genes that are co-locally conserved across various genomes, not necessarily in the same order. Their discovery and analysis is valuable in tasks such as gene annotation and prediction of gene interactions, and in the study of genome organization and evolution. The discovery of conserved gene clusters in a given set of genomes is a well studied problem, but with the rapid sequencing of prokaryotic genomes a new problem is inspired. Namely, given an already known gene cluster that was discovered and studied in one genomic dataset, to identify all the instances of the gene cluster in a given new genomic sequence. Thus, we define a new problem in comparative genomics, denoted PQ-Tree Search that takes as input a PQ-tree T representing the known gene orders of a gene cluster of interest, a gene-to-gene substitution scoring function h, integer arguments d_T and d_S, and a new sequence of genes S. The objective is to identify in S approximate new instances of the gene cluster; These instances could vary from the known gene orders by genome rearrangements that are constrained by T, by gene substitutions that are governed by h, and by gene deletions and insertions that are bounded from above by d_T and d_S, respectively. We prove that PQ-Tree Search is NP-hard and propose a parameterized algorithm that solves the optimization variant of PQ-Tree Search in O^*(2^{gamma }) time, where gamma is the maximum degree of a node in T and O^* is used to hide factors polynomial in the input size. The algorithm is implemented as a search tool, denoted PQFinder, and applied to search for instances of chromosomal gene clusters in plasmids, within a dataset of 1,487 prokaryotic genomes. We report on 29 chromosomal gene clusters that are rearranged in plasmids, where the rearrangements are guided by the corresponding PQ-trees. One of these results, coding for a heavy metal efflux pump, is further analysed to exemplify how PQFinder can be harnessed to reveal interesting new structural variants of known gene clusters.

H2-dependent formate production by hyperthermophilic Thermococcales: an alternative to sulfur reduction for reducing-equivalents disposal.

Removal of reducing equivalents is an essential catabolic process for all microorganisms to maintain their internal redox balance. The electron disposal by chemoorganotrophic Thermococcales generates H2 by proton reduction or H2S in presence of S0. Although in the absence of S0 growth of these (hyper)thermopiles was previously described to be H2-limited, it remains unclear how Thermococcales could be present in H2-rich S0-depleted habitats. Here, we report that 12 of the 47 strains tested, distributed among all three orders of Thermococcales, could grow without S0 at 0.8 mM dissolved H2 and that tolerance to H2 was always associated with formate production. Two conserved gene clusters coding for a formate hydrogenlyase (FHL) and a putative formate dehydrogenase-NAD(P)H-oxidoreductase were only present in H2-dependent formate producers, and were both systematically associated with a formate dehydrogenase and a formate transporter. As the reaction involved in this alternative pathway for disposal of reducing equivalents was close to thermodynamic equilibrium, it was strongly controlled by the substrates-products concentration ratio even in the presence of S0. Moreover, experimental data and thermodynamic modelling also demonstrated that H2-dependent CO2 reduction to formate could occur within a large temperature range in contrasted hydrothermal systems, suggesting it could also provide an adaptive advantage.

Cis-Repression of Foxq1 Expression Affects Foxf2-Mediated Gene Expression in Palate Development.

Disruption of FOXF2, encoding a member of the Forkhead family transcription factors, has been associated with cleft palate in humans and mice. FOXF2 is located in a conserved gene cluster containing FOXQ1, FOXF2, and FOXC1. We found that expression of Foxq1 is dramatically upregulated in the embryonic palatal mesenchyme in Foxf2–/– mouse embryos. We show here that the Foxf2 promoter-deletion mutation caused dramatically increased expression of the cis-linked Foxq1 allele but had little effect on the Foxq1 allele in trans. We analyzed effects of the Foxf2 mutation on the expression of other neighboring genes and compared those effects with the chromatin domain structure and recently identified enhancer-promoter associations as well as H3K27ac ChIP-seq data. We show that the Foxf2 mutation resulted in significantly increased expression of the Foxq1 and Exoc2 genes located in the same topologically associated domain with Foxf2 but not the expression of the Foxc1 and Gmds genes located in the adjacent chromatin domain. We inactivated the Foxq1 gene in mice homozygous for a Foxf2 conditional allele using CRISPR genome editing and generated (Foxf2/Foxq1)+/– mice with loss-of-function mutations in Foxf2 and Foxq1 in cis. Whereas the (Foxf2/Foxq1)–/– mice exhibited cleft palate at birth similar as in the Foxf2–/– mice, systematic expression analyses of a large number of Foxf2-dependent genes revealed that the (Foxf2/Foxq1)–/– embryos exhibited distinct effects on the domain-specific expression of several important genes, including Foxf1, Shox2, and Spon1, in the developing palatal shelves compared with Foxf2–/– embryos. These results identify a novel cis-regulatory effect of the Foxf2 mutation and demonstrate that cis-regulation of Foxq1 contributed to alterations in palatal gene expression in Foxf2–/– embryos. These results have important implications for interpretation of results and mechanisms from studies of promoter- or gene-deletion alleles. In addition, the unique mouse lines generated in this study provide a valuable resource for understanding the cross-regulation and combinatorial functions of the Foxf2 and Foxq1 genes in development and disease.

GCsnap: Interactive Snapshots for the Comparison of Protein-Coding Genomic Contexts

The biological function and evolutionary history of protein-coding genes are not only written in their nucleotide sequences. The comparison of genomic contexts throughout different lineages may highlight genomic mechanisms in the generation of new protein families, while the conservation of gene clusters may unravel, for instance, metabolic pathways. Various tools and databases exist that allow for the analysis and comparison of genomic contexts, but each has its own limitations. Online databases allow for quick comparisons, but only for those genomes for which data were pre-calculated. More advanced tools may allow for the comparison of any genome, but are often limited to a given phylogenetic kingdom or provide only a snapshot of the genomic contexts without further information about the genes involved. Here, we introduce GCsnap, a flexible Python-based tool that allows for the interactive comparison of the genomic contexts of protein-coding genes from any genome at any taxonomic level, integrating them with functional and structural information. By connecting the output to different protein databases, users can navigate through the different genomic contexts from a simple interactive platform, facilitating the further analysis of the contexts found. GCsnap is not limited to a single input format, can perform batch jobs and accepts protein classification maps. Results are stored in detailed, human and machine-readable files, and customizable, publication-ready figures. GCsnap is freely available from https://github.com/JoanaMPereira/GCsnap and can be set up easily on any computer.

A New Dinoflagellate Genome Illuminates a Conserved Gene Cluster Involved in Sunscreen Biosynthesis.

Photosynthetic dinoflagellates of the Family Symbiodiniaceae live symbiotically with many organisms that inhabit coral reefs and are currently classified into fifteen groups, including seven genera. Draft genomes from four genera, Symbiodinium, Breviolum, Fugacium, and Cladocopium, which have been isolated from corals, have been reported. However, no genome is available from the genus Durusdinium, which occupies an intermediate phylogenetic position in the Family Symbiodiniaceae and is well known for thermal tolerance (resistance to bleaching). We sequenced, assembled, and annotated the genome of Durusdinium trenchii, isolated from the coral, Favia speciosa, in Okinawa, Japan. Assembled short reads amounted to 670 Mb with ∼47% GC content. This GC content was intermediate among taxa belonging to the Symbiodiniaceae. Approximately 30,000 protein-coding genes were predicted in the D. trenchii genome, fewer than in other genomes from the Symbiodiniaceae. However, annotations revealed that the D. trenchii genome encodes a cluster of genes for synthesis of mycosporine-like amino acids, which absorb UV radiation. Interestingly, a neighboring gene in the cluster encodes a glucose–methanol–choline oxidoreductase with a flavin adenine dinucleotide domain that is also found in Symbiodinium tridacnidorum. This conservation seems to partially clarify an ancestral genomic structure in the Symbiodiniaceae and its loss in late-branching lineages, including Breviolum and Cladocopium, after splitting from the Durusdinium lineage. Our analysis suggests that approximately half of the taxa in the Symbiodiniaceae may maintain the ability to synthesize mycosporine-like amino acids. Thus, this work provides a significant genomic resource for understanding the genomic diversity of Symbiodiniaceae in corals.

Characteristics of two myoviruses induced from the coastal photoheterotrophic bacterium Porphyrobacter sp. YT40.

In this study, we characterized two induced myoviruses from one marine photoheterotrophic bacterium Porphyrobacter sp. YT40 belonging to the Sphingomonadales family in Alphaproteobacteria. The genome sequence of prophage A is ∼36.9kb with an average GC content of 67.1%, and its core or functional genes are homologous to Mu or Mu-like phages. Furthermore, induced viral particles from prophage A show a knob-like neck structure, which is only found in bacteriophage Mu. The genome size of prophage B is ∼36.8kb with an average GC content of 65.3%. Prophage B contains a conserved gene cluster Q-P-O-N-M-L, which is unique in P2 phages. Induced viral particles from prophage B display an icosahedral head with a diameter of ∼55nm and a 130±5nm long contractile tail. To our knowledge, this is the first report that characterizes the induced P2-like phage in marine Alphaproteobacteria. Phylogeny analyses suggest that these two types of prophages are commonly found in sequenced bacteria of the Sphingomonadales family. This study sheds light on the ongoing interaction between marine bacteria and phages, and improves our understanding of bacterial genomic plasticity and evolution.

EvolClust: automated inference of evolutionary conserved gene clusters in eukaryotes.

MotivationThe evolution and role of gene clusters in eukaryotes is poorly understood. Currently, most studies and computational prediction programs limit their focus to specific types of clusters, such as those involved in secondary metabolism.ResultsWe present EvolClust, a python-based tool for the inference of evolutionary conserved gene clusters from genome comparisons, independently of the function or gene composition of the cluster. EvolClust predicts conserved gene clusters from pairwise genome comparisons and infers families of related clusters from multiple (all versus all) genome comparisons.Availability and implementation https://github.com/Gabaldonlab/EvolClust/.Supplementary information Supplementary data are available at Bioinformatics online.

Evolutionary and functional patterns of shared gene neighbourhood in fungi.

Gene clusters comprise genomically co-localized and potentially co-regulated genes that tend to be conserved across species. In eukaryotes, multiple examples of metabolic gene clusters are known, particularly among fungi and plants. However, little is known about how gene clustering patterns vary among taxa or with respect to functional roles. Furthermore, mechanisms of the formation, maintenance and evolution of gene clusters remain unknown. We surveyed 341 fungal genomes to discover gene clusters shared by different species, independently of their functions. We inferred 12,120 cluster families, which comprised roughly one third of the gene space and were enriched in genes associated with diverse cellular functions. Additionally, most clusters did not encode transcription factors, suggesting that they are regulated distally. We used phylogenomics to characterize the evolutionary history of these clusters. We found that most clusters originated once and were transmitted vertically, coupled to differential loss. However, convergent evolution-that is, independent appearance of the same cluster-was more prevalent than anticipated. Finally, horizontal gene transfer of entire clusters was somewhat restricted, with the exception of those associated with secondary metabolism. Altogether, our results provide insights on the evolution of gene clustering as well as a broad catalogue of evolutionarily conserved gene clusters whose function remains to be elucidated.

Abstract 261: L1-MET transcription silencing modulates MET and EGFR gene and their protein expression and induces apoptosis and cell-death in different types of cancer cells

Abstract The activation of the LINE-1 sequence located within the second intron of MET leads to the onset of L1-MET transcript. We recently characterized the full L1-MET structure in breast cancer and showed that high levels of this transcript recognize a subset of more aggressive breast carcinomas, mainly of triple negative phenotype. However, at present, the relationship between L1-MET and MET is still poorly understood. In order to elucidate this function, we silenced L1-MET transcription using cells expressing different levels of L1-MET/MET, including lung cancer (A549 and EBC1), gastric cancer (GTL16), and breast cancer (MDA-MB231), that were transiently transfected with Gapmers-LNA (Exiqon), specifically targeting L1-MET sequence. Cell viability and apoptosis were evaluated after 24 h by cell count, Cell Titer Glow (Promega) and propidium iodide/annexin based-cytofluorimeter assays. RNA was purified from sample and control cells to assess the L1-METsilencing by qRT-PCR and to evaluate the gene-expression of a subset of cancer-related genes using Nanostring Technology, whereas western blot analyses were carried out to measure the protein expressions. A significant decrease of cell viability was detected in A549, EBC1 and GTL16, but not in MDA-MB231 cells, characterized by the lowest level of L1-MET overall. In parallel, the highly expressing L1-MET cells showed an increased rate of early and late apoptosis together with a strong reduction of MET gene and its protein expression. On the contrary, in MDA-MB231 cells L1-MET silencing induced only a slight MET gene and protein impairment. Overall, L1-MET knock-down caused a decrease expression of a conserved gene cluster, including a marked reduction of EGFR protein expression. Moreover, L1-MET silenced cells showed lower MET and EGFR phosphorylation, with a downstream silencing effect on pERK and pAKT. Results of cell treatment with the inhibitors of lysosome and proteasome activity bafilomycin and MG-132 ruled out the interaction of L1-MET silencing with protein degradation pathway. This is the first study investigating the function of the L1-MET transcript in cancer models. Our results show that although L1-MET is unable to encode for a protein, its silencing exerts a strong phenotypic effect on different tumor cell types, suggesting potential regulations at the transcriptional level. Note: This abstract was not presented at the meeting. Citation Format: Enrico Berrino, Umberto Miglio, Valentina Miano, Letizia Lanzetti, Silvia Benvenuti, Carla Debernardi, MIchele De Bortoli, Caterina Marchio, Tiziana Venesio, Anna Sapino. L1-MET transcription silencing modulates MET and EGFR gene and their protein expression and induces apoptosis and cell-death in different types of cancer cells [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 261.

Complete Genomes of Symbiotic Cyanobacteria Clarify the Evolution of Vanadium-Nitrogenase.

Plant endosymbiosis with nitrogen-fixing cyanobacteria has independently evolved in diverse plant lineages, offering a unique window to study the evolution and genetics of plant–microbe interaction. However, very few complete genomes exist for plant cyanobionts, and therefore little is known about their genomic and functional diversity. Here, we present four complete genomes of cyanobacteria isolated from bryophytes. Nanopore long-read sequencing allowed us to obtain circular contigs for all the main chromosomes and most of the plasmids. We found that despite having a low 16S rRNA sequence divergence, the four isolates exhibit considerable genome reorganizations and variation in gene content. Furthermore, three of the four isolates possess genes encoding vanadium (V)-nitrogenase (vnf), which is uncommon among diazotrophs and has not been previously reported in plant cyanobionts. In two cases, the vnf genes were found on plasmids, implying possible plasmid-mediated horizontal gene transfers. Comparative genomic analysis of vnf-containing cyanobacteria further identified a conserved gene cluster. Many genes in this cluster have not been functionally characterized and would be promising candidates for future studies to elucidate V-nitrogenase function and regulation.

Characterization of l-Carnitine Metabolism in Sinorhizobium meliloti.

l-Carnitine is a trimethylammonium compound mostly known for its contribution to fatty acid transport into mitochondria. In bacteria, it is synthesized from γ-butyrobetaine (GBB) and can be used as a carbon source. l-Carnitine can be formed directly by GBB hydroxylation or synthesized via a biosynthetic route analogous to fatty acid degradation. However, this multistep pathway has not been experimentally characterized. In this work, we identified by gene context analysis a cluster of l-carnitine anabolic genes next to those involved in its catabolism and proceeded to the complete in vitro characterization of l-carnitine biosynthesis and degradation in Sinorhizobium meliloti The five enzymes catalyzing the seven steps that convert GBB to glycine betaine are described. Metabolomic analysis confirmed the multistage synthesis of l-carnitine in GBB-grown cells but also revealed that GBB is synthesized by S. meliloti To our knowledge, this is the first report of aerobic GBB synthesis in bacteria. The conservation of l-carnitine metabolism genes in different bacterial taxonomic classes underscores the role of l-carnitine as a ubiquitous nutrient.IMPORTANCE The experimental characterization of novel metabolic pathways is essential for realizing the value of genome sequences and improving our knowledge of the enzymatic capabilities of the bacterial world. However, 30% to 40% of genes of a typical genome remain unannotated or associated with a putative function. We used enzyme kinetics, liquid chromatography-mass spectroscopy (LC-MS)-based metabolomics, and mutant phenotyping for the characterization of the metabolism of l-carnitine in Sinorhizobium meliloti to provide an accurate annotation of the corresponding genes. The occurrence of conserved gene clusters for carnitine metabolism in soil, plant-associated, and marine bacteria underlines the environmental abundance of carnitine and suggests this molecule might make a significant contribution to ecosystem nitrogen and carbon cycling.

Functional reconstitution of the type IVa pilus assembly system from enterohaemorrhagic Escherichia coli.

Type 4a pili (T4aP) are long, thin and dynamic fibres displayed on the surface of diverse bacteria promoting adherence, motility and transport functions. Genomes of many Enterobacteriaceae contain conserved gene clusters encoding putative T4aP assembly systems. However, their expression has been observed only in few strains including Enterohaemorrhagic Escherichia coli (EHEC) and their inducers remain unknown. Here we used EHEC genomic DNA as a template to amplify and assemble an artificial operon composed of four gene clusters encoding 13 pilus assembly proteins. Controlled expressions of this operon in nonpathogenic E. coli strains led to efficient assembly of T4aP composed of the major pilin PpdD, as shown by shearing assays and immunofluorescence microscopy. When compared with PpdD pili assembled in a heterologous Klebsiella T2SS type 2 secretion system (T2SS) by using cryo-electron microscopy (cryoEM), these pili showed indistinguishable helical parameters, emphasizing that major pilins are the principal determinants of the fibre structure. Bacterial two-hybrid analysis identified several interactions of PpdD with T4aP assembly proteins, and with components of the T2SS that allow for heterologous fibre assembly. These studies lay ground for further characterization of the T4aP structure, function and biogenesis in enterobacteria.

Genomic discovery of the hypsin gene and biosynthetic pathways for terpenoids in Hypsizygus marmoreus

BackgroundHypsizygus marmoreus (Beech mushroom) is a popular ingredient in Asian cuisine. The medicinal effects of its bioactive compounds such as hypsin and hypsiziprenol have been reported, but the genetic basis or biosynthesis of these components is unknown.ResultsIn this study, we sequenced a reference strain of H. marmoreus (Haemi 51,987–8). We evaluated various assembly strategies, and as a result the Allpaths and PBJelly produced the best assembly. The resulting genome was 42.7 Mbp in length and annotated with 16,627 gene models. A putative gene (Hypma_04324) encoding the antifungal and antiproliferative hypsin protein with 75% sequence identity with the previously known N-terminal sequence was identified. Carbohydrate active enzyme analysis displayed the typical feature of white-rot fungi where auxiliary activity and carbohydrate-binding modules were enriched. The genome annotation revealed four terpene synthase genes responsible for terpenoid biosynthesis. From the gene tree analysis, we identified that terpene synthase genes can be classified into six clades. Four terpene synthase genes of H. marmoreus belonged to four different groups that implies they may be involved in the synthesis of different structures of terpenes. A terpene synthase gene cluster was well-conserved in Agaricomycetes genomes, which contained known biosynthesis and regulatory genes.ConclusionsGenome sequence analysis of this mushroom led to the discovery of the hypsin gene. Comparative genome analysis revealed the conserved gene cluster for terpenoid biosynthesis in the genome. These discoveries will further our understanding of the biosynthesis of medicinal bioactive molecules in this edible mushroom.

Comparative Genomics among Closely Related Streptomyces Strains Revealed Specialized Metabolite Biosynthetic Gene Cluster Diversity

Specialized metabolites are of great interest due to their possible industrial and clinical applications. The increasing number of antimicrobial resistant infectious agents is a major health threat and therefore, the discovery of chemical diversity and new antimicrobials is crucial. Extensive genomic data from Streptomyces spp. confirm their production potential and great importance. Genome sequencing of the same species strains indicates that specialized metabolite biosynthetic gene cluster (SMBGC) diversity is not exhausted, and instead, a pool of novel specialized metabolites still exists. Here, we analyze the genome sequence data from six phylogenetically close Streptomyces strains. The results reveal that the closer strains are phylogenetically, the number of shared gene clusters is higher. Eight specialized metabolites comprise the core metabolome, although some strains have only six core gene clusters. The number of conserved gene clusters common between the isolated strains and their closest phylogenetic counterparts varies from nine to 23 SMBGCs. However, the analysis of these phylogenetic relationships is not affected by the acquisition of gene clusters, probably by horizontal gene transfer events, as each strain also harbors strain-specific SMBGCs. Between one and 15 strain-specific gene clusters were identified, of which up to six gene clusters in a single strain are unknown and have no identifiable orthologs in other species, attesting to the existing SMBGC novelty at the strain level.

The complete mitochondrial genome of the early flowering plant Nymphaea colorata is highly repetitive with low recombination

BackgroundMitochondrial genomes of flowering plants (angiosperms) are highly dynamic in genome structure. The mitogenome of the earliest angiosperm Amborella is remarkable in carrying rampant foreign DNAs, in contrast to Liriodendron, the other only known early angiosperm mitogenome that is described as ‘fossilized’. The distinctive features observed in the two early flowering plant mitogenomes add to the current confusions of what early flowering plants look like. Expanded sampling would provide more details in understanding the mitogenomic evolution of early angiosperms. Here we report the complete mitochondrial genome of water lily Nymphaea colorata from Nymphaeales, one of the three orders of the earliest angiosperms.ResultsAssembly of data from Pac-Bio long-read sequencing yielded a circular mitochondria chromosome of 617,195 bp with an average depth of 601×. The genome encoded 41 protein coding genes, 20 tRNA and three rRNA genes with 25 group II introns disrupting 10 protein coding genes. Nearly half of the genome is composed of repeated sequences, which contributed substantially to the intron size expansion, making the gross intron length of the Nymphaea mitochondrial genome one of the longest among angiosperms, including an 11.4-Kb intron in cox2, which is the longest organellar intron reported to date in plants. Nevertheless, repeat mediated homologous recombination is unexpectedly low in Nymphaea evidenced by 74 recombined reads detected from ten recombinationally active repeat pairs among 886,982 repeat pairs examined. Extensive gene order changes were detected in the three early angiosperm mitogenomes, i.e. 38 or 44 events of inversions and translocations are needed to reconcile the mitogenome of Nymphaea with Amborella or Liriodendron, respectively. In contrast to Amborella with six genome equivalents of foreign mitochondrial DNA, not a single horizontal gene transfer event was observed in the Nymphaea mitogenome.ConclusionsThe Nymphaea mitogenome resembles the other available early angiosperm mitogenomes by a similarly rich 64-coding gene set, and many conserved gene clusters, whereas stands out by its highly repetitive nature and resultant remarkable intron expansions. The low recombination level in Nymphaea provides evidence for the predominant master conformation in vivo with a highly substoichiometric set of rearranged molecules.

Intrageneric chloroplast genome comparison in the genus Euglena (Phylum: Euglenophyta) with annotated chloroplast genomes of Euglena hiemalis and Euglena clara

The genus Euglena is composed of six subclades with a diversity of chloroplast morphologies, unlike that found in the other genera of Euglenaceae. This genus contains five published chloroplast genomes (cpGenome) and the colorless plastid genome (ptGenome) of the non-photosynthetic euglenoid Euglena longa within the same subclade as Euglena gracilis. Previous studies of cpGenomes revealed that Euglena viridis and E. gracilis, although in separate subclades, had few cpGenome differences, while Euglena mutabilis maintained the same gene order but was mirror-inverted except for the rRNA cluster. However, we expanded the number of cpGenomes available in Euglena by sequencing and annotating the cpGenomes of Euglena clara, the earliest diverging species in the E. gracilis and E. longa subclade and Euglena hiemalis, the putative sister species to E. longa. Analysis of these newly annotated cpGenomes showed them to be largely similar in gene content, conserved gene clusters (operons), and G + C/A + T percentage to previously published Euglena cpGenomes. The only extensive gene rearrangements observed were between cpGenomes and the ptGenome in subclade B. However, a unique feature of the subclade was multi-copies of the rRNA operon. Also, homologous twintrons in psbD and psbF to E. gracilis were observed in E. hiemalis. Overall, these results revealed a conserved intrageneric cpGenome for Euglena species.

Linkage of the Nit1C gene cluster to bacterial cyanide assimilation as a nitrogen source.

A genetic linkage between a conserved gene cluster (Nit1C) and the ability of bacteria to utilize cyanide as the sole nitrogen source was demonstrated for nine different bacterial species. These included three strains whose cyanide nutritional ability has formerly been documented (Pseudomonas fluorescens Pf11764, Pseudomonas putida BCN3 and Klebsiella pneumoniae BCN33), and six not previously known to have this ability [Burkholderia (Paraburkholderia) xenovorans LB400, Paraburkholderia phymatum STM815, Paraburkholderia phytofirmans PsJN, Cupriavidus (Ralstonia) eutropha H16, Gluconoacetobacter diazotrophicus PA1 5 and Methylobacterium extorquens AM1]. For all bacteria, growth on or exposure to cyanide led to the induction of the canonical nitrilase (NitC) linked to the gene cluster, and in the case of Pf11764 in particular, transcript levels of cluster genes (nitBCDEFGH) were raised, and a nitC knock-out mutant failed to grow. Further studies demonstrated that the highly conserved nitB gene product was also significantly elevated. Collectively, these findings provide strong evidence for a genetic linkage between Nit1C and bacterial growth on cyanide, supporting use of the term cyanotrophy in describing what may represent a new nutritional paradigm in microbiology. A broader search of Nit1C genes in presently available genomes revealed its presence in 270 different bacteria, all contained within the domain Bacteria, including Gram-positive Firmicutes and Actinobacteria, and Gram-negative Proteobacteria and Cyanobacteria. Absence of the cluster in the Archaea is congruent with events that may have led to the inception of Nit1C occurring coincidentally with the first appearance of cyanogenic species on Earth, dating back 400-500 million years.