Large Gene Cluster Research Articles

Cross-Talking of Pathway-Specific Regulators in Glycopeptide Antibiotics (Teicoplanin and A40926) Production.

Teicoplanin and A40926 (natural precursor of dalbavancin) are clinically relevant glycopeptide antibiotics (GPAs) produced by Actinoplanes teichomyceticus NRRL B-16726 and Nonomuraea gerenzanensis ATCC 39727. Their biosynthetic enzymes are coded within large biosynthetic gene clusters (BGCs), named tei for teicoplanin and dbv for A40926, whose expression is strictly regulated by pathway-specific transcriptional regulators (PSRs), coded by cluster-situated regulatory genes (CSRGs). Herein, we investigated the "cross-talk" between the CSRGs from tei and dbv, through the analysis of GPA production levels in A. teichomyceticus and N. gerenzanensis strains, with knockouts of CSRGs cross-complemented by the expression of heterologous CSRGs. We demonstrated that Tei15* and Dbv4 StrR-like PSRs, although orthologous, were not completely interchangeable: tei15* and dbv4 were only partially able or unable to cross-complement N. gerenzanensis knocked out in dbv4 and A. teichomyceticus knocked out in tei15*, implying that the DNA-binding properties of these PSRs are more different in vivo than it was believed before. At the same time, the unrelated LuxR-like PSRs Tei16* and Dbv3 were able to cross-complement corresponding N. gerenzanensis knocked out in dbv3 and A. teichomyceticus knocked out in tei16*. Moreover, the heterologous expression of dbv3 in A. teichomyceticus led to a significant increase in teicoplanin production. Although the molecular background of these events merits further investigations, our results contribute to a deeper understanding of GPA biosynthesis regulation and offer novel biotechnological tools to improve their production.

Five copper homeostasis gene clusters encode the Cu-efflux resistome of the highly copper-tolerant Methylorubrum extorquens AM1.

In the last decade, the use of copper has reemerged as a potential strategy to limit healthcare-associated infections and to control the spread of multidrug-resistant pathogens. Numerous environmental studies have proposed that most opportunistic pathogens have acquired antimicrobial resistance in their nonclinical primary habitat. Thus, it can be presumed that copper-resistant bacteria inhabiting a primary commensal niche might potentially colonize clinical environments and negatively affect the bactericidal efficacy of Cu-based treatments. The use of copper in agricultural fields is one of the most important sources of Cu pollution that may exert selection pressure for the increase of copper resistance in soil and plant-associated bacteria. To assess the emergence of copper-resistant bacteria in natural habitats, we surveyed a laboratory collection of bacterial strains belonging to the order Rhizobiales. This study proposes that Methylorubrum extorquens AM1 is an environmental isolate well adapted to thrive in copper-rich environments that could act as a reservoir of copper resistance genes. The minimal inhibitory concentrations (MICs) of CuCl2 were used to estimate the copper tolerance of eight plant-associated facultative diazotrophs (PAFD) and five pink-pigmented facultative methylotrophs (PPFM) belonging to the order Rhizobiales presumed to come from nonclinical and nonmetal-polluted natural habitats based on their reported source of isolation. Their sequenced genomes were used to infer the occurrence and diversity of Cu-ATPases and the copper efflux resistome of Mr. extorquens AM1. These bacteria exhibited minimal inhibitory concentrations (MICs) of CuCl2 ranging between 0.020 and 1.9 mM. The presence of multiple and quite divergent Cu-ATPases per genome was a prevalent characteristic. The highest copper tolerance exhibited by Mr. extorquens AM1 (highest MIC of 1.9 mM) was similar to that found in the multimetal-resistant model bacterium Cupriavidus metallidurans CH34 and in clinical isolates of Acinetobacter baumannii. The genome-predicted copper efflux resistome of Mr. extorquens AM1 consists of five large (6.7 to 25.7 kb) Cu homeostasis gene clusters, three clusters share genes encoding Cu-ATPases, CusAB transporters, numerous CopZ chaperones, and enzymes involved in DNA transfer and persistence. The high copper tolerance and the presence of a complex Cu efflux resistome suggest the presence of relatively high copper tolerance in environmental isolates of Mr. extorquens.

Chromosome-level Genome Assembly of Euphorbia peplus, a Model System for Plant Latex, Reveals that Relative Lack of Ty3 Transposons Contributed to Its Small Genome Size.

Euphorbia peplus (petty spurge) is a small, fast-growing plant that is native to Eurasia and has become a naturalized weed in North America and Australia. Euphorbia peplus is not only medicinally valuable, serving as a source for the skin cancer drug ingenol mebutate, but also has great potential as a model for latex production owing to its small size, ease of manipulation in the laboratory, and rapid reproductive cycle. To help establish E. peplus as a new model, we generated a 267.2-Mb Hi-C-anchored PacBio HiFi nuclear genome assembly with a BUSCO score of 98.5%, a genome annotation based on RNA-seq data from six organs, and publicly accessible tools including a genome browser and an interactive organ-specific expression atlas. Chromosome number is highly variable across Euphorbia species. Using a comparative analysis of our newly sequenced E. peplus genome with other Euphorbiaceae genomes, we show that variation in Euphorbia chromosome number between E. peplus and Euphorbia lathyris is likely due to fragmentation and rearrangement rather than chromosomal duplication followed by diploidization of the duplicated sequence. Moreover, we found that the E. peplus genome is relatively compact compared with related members of the genus in part due to restricted expansion of the Ty3 transposon family. Finally, we identify a large gene cluster that contains many previously identified enzymes in the putative ingenol mebutate biosynthesis pathway, along with additional gene candidates for this biosynthetic pathway. The genomic resources we have created for E. peplus will help advance research on latex production and ingenol mebutate biosynthesis in the commercially important Euphorbiaceae family.

Neuron-restrictive silencer factor/repressor element 1-silencing transcription factor (NRSF/REST) controls spatial K+ buffering in primary cortical astrocytes.

Neuron-restrictive silencer factor/repressor element 1 (RE1)-silencing transcription factor (NRSF/REST) is a transcriptional repressor of a large cluster of neural genes containing RE1 motifs in their promoter region. NRSF/REST is ubiquitously expressed in non-neuronal cells, including astrocytes, while it is down-regulated during neuronal differentiation. While neuronal NRSF/REST homeostatically regulates intrinsic excitability and synaptic transmission, the role of the high NRSF/REST expression levels in the homeostatic functions of astrocytes is poorly understood. Here, we investigated the functional consequences of NRSF/REST deletion in primary cortical astrocytes derived from NRSF/REST conditional knockout mice (KO). We found that NRSF/REST KO astrocyte displayed a markedly reduced activity of inward rectifying K+ channels subtype 4.1 (Kir4.1) underlying spatial K+ buffering that was associated with a decreased expression and activity of the glutamate transporter-1 (GLT-1) responsible for glutamate uptake by astrocytes. The effects of the impaired astrocyte homeostatic functions on neuronal activity were investigated by co-culturing wild-type hippocampal neurons with NRSF/REST KO astrocytes. Interestingly, neurons experienced increased neuronal excitability at high firing rates associated with decrease after hyperpolarization and increased amplitude of excitatory postsynaptic currents. The data indicate that astrocytic NRSF/REST directly participates in neural circuit homeostasis by regulating intrinsic excitability and excitatory transmission and that dysfunctions of NRSF/REST expression in astrocytes may contribute to the pathogenesis of neurological disorders.

Mechanic Insight into the Distinct and Common Roles of Ovariectomy Versus Adrenalectomy on Adipose Tissue Remodeling in Female Mice.

Dysfunctions of the ovaries and adrenal glands are both evidenced to cause aberrant adipose tissue (AT) remodeling and resultant metabolic disorders, but their distinct and common roles are poorly understood. In this study, through biochemical, histological and RNA-seq analyses, we comprehensively explored the mechanisms underpinning subcutaneous (SAT) and visceral adipose tissue (VAT) remodeling, in response to ovariectomy (OVX) versus adrenalectomy (ADX) in female mice. OVX promoted adipocyte differentiation and fat accumulation in both SAT and VAT, by potentiating the Pparg signaling, while ADX universally prevented the cell proliferation and extracellular matrix organization in both SAT and VAT, likely by inactivating the Nr3c1 signaling, thus causing lipoatrophy in females. ADX, but not OVX, exerted great effects on the intrinsic difference between SAT and VAT. Specifically, ADX reversed a large cluster of genes differentially expressed between SAT and VAT, by activating 12 key transcription factors, and thereby caused senescent cell accumulation, massive B cell infiltration and the development of selective inflammatory response in SAT. Commonly, both OVX and ADX enhance circadian rhythmicity in VAT, and impair cell proliferation, neurogenesis, tissue morphogenesis, as well as extracellular matrix organization in SAT, thus causing dysfunction of adipose tissues and concomitant metabolic disorders.

The complete mitochondrial genome and novel gene arrangement in Nesodiprion zhejiangensis Zhou & Xiao (Hymenoptera: Diprionidae).

The complete mitochondrial genome (mitogenome) of the sawfly, Nesodiprion zhejiangensis Zhou & Xiao, was sequenced, assembled, and deposited in GenBank (Accession Number: OM501121). The 15,660bpN. zhejiangensis mitogenome encodes for 2 ribosomal RNAs (rrnL and rrnS), 22 transfer RNAs (tRNAs), 13 protein-coding genes (PCGs), and an AT-rich region of 450bp in length. The nucleotide composition is biased toward adenine and thymine (A + T = 81.8%). Each PCG is initiated by an ATN codon, except for cox2, which starts with a TTG. Of 13 PCGs, 9 have a TAA termination codon, while the remainder terminate with a TAG or a single T. All tRNAs have the classic cloverleaf structure, except for the dihydrouridine (DHU) arm of tRNAval, which forms a simple loop. There are 49 helices belonging to 6 domains in rrnL and 30 helices belonging to 4 domains in rrnS. In comparison to the ancestral architecture, N. zhejiangensis has the most rearranged mitogenome in Symphyta, in which rearrangement events of local inversion and transposition are identified in three gene clusters. Specifically, the main hotspot of gene rearrangement occurred between rrnS and trnY, and rearranged from rrnS-(AT-rich region)-I-Q-M-nd2-W-C-Y to rrnS-Q-W-C-nd2-I-M-(AT-rich region)-Y, involving a local inversion event of a large gene cluster and transposition events of some tRNAs. Transposition of trnA and trnR (rearranged from A-R to R-A) was observed at the nd3-nd5 gene junction while shuffling of trnP and trnT (rearranged from T-P to P-T) occurred at the nd4l-nd6 gene junction. While illegitimate inter-mtDNA recombination might explain the opposite orientations of transcription between rrnS and trnY, transposition events of tRNA in some gene blocks can be accounted for by the tandem duplication/random loss (TDRL) model. Our phylogenetic analysis suggests that N. zhejiangensis is closely related to congeneric species N. biremis and N. japonicus, which together form a sister lineage with the European pine sawfly, Neodiprion sertifer.

A New Face of the Old Gene: Deletion of the PssA, Encoding Monotopic Inner Membrane Phosphoglycosyl Transferase in Rhizobium leguminosarum, Leads to Diverse Phenotypes That Could Be Attributable to Downstream Effects of the Lack of Exopolysaccharide.

The biosynthesis of subunits of rhizobial exopolysaccharides is dependent on glycosyltransferases, which are usually encoded by large gene clusters. PssA is a member of a large family of phosphoglycosyl transferases catalyzing the transfer of a phosphosugar moiety to polyprenol phosphate; thus, it can be considered as priming glycosyltransferase commencing synthesis of the EPS repeating units in Rhizobium leguminosarum. The comprehensive analysis of PssA protein features performed in this work confirmed its specificity for UDP-glucose and provided evidence that PssA is a monotopic inner membrane protein with a reentrant membrane helix rather than a transmembrane segment. The bacterial two-hybrid system screening revealed interactions of PssA with some GTs involved in the EPS octasaccharide synthesis. The distribution of differentially expressed genes in the transcriptome of the ΔpssA mutant into various functional categories indicated complexity of cell response to the deletion, which can mostly be attributed to the lack of exopolysaccharide and downstream effects caused by such deficiency. The block in the EPS biosynthesis at the pssA step, potentially leading to an increased pool of UDP-glucose, is likely to be filtered through to other pathways, and thus the absence of EPS may indirectly affect the expression of proteins involved in these pathways.

The gap-free potato genome assembly reveals large tandem gene clusters of agronomical importance in highly repeated genomic regions

The gap-free potato genome assembly reveals large tandem gene clusters of agronomical importance in highly repeated genomic regions

SALL4: An Intriguing Therapeutic Target in Cancer Treatment.

Spalt-Like Transcription Factor 4 (SALL4) is a critical factor for self-renewal ability and pluripotency of stem cells. On the other hand, various reports show tight relation of SALL4 to cancer occurrence and metastasis. SALL4 exerts its effects not only by inducing gene expression but also repressing a large cluster of genes through interaction with various epigenetic modifiers. Due to high expression of SALL4 in cancer cells and its silence in almost all adult tissues, it is an ideal target for cancer therapy. However, targeting SALL4 meets various challenges. SALL4 is a transcription factor and designing appropriate drug to inhibit this intra-nucleus component is challenging. On the other hand, due to lack of our knowledge on structure of the protein and the suitable active sites, it becomes more difficult to reach the appropriate drugs against SALL4. In this review, we have focused on approaches applied yet to target this oncogene and discuss the potential of degrader systems as new therapeutics to target oncogenes.

Two large novel alpha-globin gene cluster deletions causing alpha(0)-thalassemia in two Chinese families

IntroductionMonosomy of terminal 16p13.3 is a relatively common subtelomeric abnormality, most affected individuals presented α-thalassemia, some also have mental retardation, developmental abnormalities and/or speech delay and facial dysmorphism, which is termed ATR-16 syndrome. Here, we reported two novel 16p13.3 deletions involving the α-globin gene cluster and multispecies conserved sequences (MCSs), causing only a phenotype of α-thalassemia. MethodsSamples were collected from members of the two families and were subjected to haematological and comprehensive genetic analysis. ResultsThe novel 108 Kb deletion in family A extends from the non-protein coding RNA gene (WASIR2) to the NPRL3 gene, removing MCS-R1 to R3. This deletion should arise de novo because it wasn’t detected in both parents. The novel 336 Kb deletion in family B should extend from telomere to ∼ chr16:336000, removing the entire α-globin gene cluster. Carriers of these two deletions presented with microcytosis and hypochromic red cells, in accordance with a phenotype of α0-thalassemia carrier. ConclusionOur study increases the mutation spectrum of α-thalassemia. MCSs deletion should be considered in clinical practice of thalassemia screening and diagnosis.

Breakthrough in efficient cloning and activation of large cryptic biosynthetic gene clusters from high GC actinobacteria

Breakthrough in efficient cloning and activation of large cryptic biosynthetic gene clusters from high GC actinobacteria

Characterization of gene cluster heterogeneity in single-cell transcriptomic data within and across cancer types.

ABSTRACTDespite the remarkable progress in probing tumor transcriptomic heterogeneity by single-cell RNA sequencing (sc-RNAseq) data, several gaps exist in prior studies. Tumor heterogeneity is frequently mentioned but not quantified. Clustering analyses typically target cells rather than genes, and differential levels of transcriptomic heterogeneity of gene clusters are not characterized. Relations between gene clusters inferred from multiple datasets remain less explored. We provided a series of quantitative methods to analyze cancer sc-RNAseq data. First, we proposed two quantitative measures to assess intra-tumoral heterogeneity/homogeneity. Second, we established a hierarchy of gene clusters from sc-RNAseq data, devised an algorithm to reduce the gene cluster hierarchy to a compact structure, and characterized the gene clusters with functional enrichment and heterogeneity. Third, we developed an algorithm to align the gene cluster hierarchies from multiple datasets to a small number of meta gene clusters. By applying these methods to nine cancer sc-RNAseq datasets, we discovered that cancer cell transcriptomes were more homogeneous within tumors than the accompanying normal cells. Furthermore, many gene clusters from the nine datasets were aligned to two large meta gene clusters, which had high and low heterogeneity and were enriched with distinct functions. Finally, we found the homogeneous meta gene cluster retained stronger expression coherence and associations with survival times in bulk level RNAseq data than the heterogeneous meta gene cluster, yet the combinatorial expression patterns of breast cancer subtypes in bulk level data were not preserved in single-cell data. The inference outcomes derived from nine cancer sc-RNAseq datasets provide insights about the contributing factors for transcriptomic heterogeneity of cancer cells and complex relations between bulk level and single-cell RNAseq data. They demonstrate the utility of our methods to enable a comprehensive characterization of co-expressed gene clusters in a wide range of sc-RNAseq data in cancers and beyond.

Genome mining of Burkholderia ambifaria strain T16, a rhizobacterium able to produce antimicrobial compounds and degrade the mycotoxin fusaric acid.

Burkholderia ambifaria T16 is a bacterium isolated from the rhizosphere of barley plants that showed a remarkable antifungal activity. This strain was also able to degrade fusaric acid (5-Butylpyridine-2-carboxylic acid) and detoxify this mycotoxin in inoculated barley seedlings. Genes and enzymes responsible for fusaric acid degradation have an important biotechnological potential in the control of fungal diseases caused by fusaric acid producers, or in the biodegradation/bio catalysis processes of pyridine derivatives. In this study, the complete genome of B. ambifaria T16 was sequenced and analyzed to identify genes involved in survival and competition in the rhizosphere, plant growth promotion, fungal growth inhibition, and degradation of aromatic compounds. The genomic analysis revealed the presence of several operons for the biosynthesis of antimicrobial compounds, such as pyrrolnitrin, ornibactin, occidiofungin and the membrane-associated AFC-BC11. These compounds were also detected in bacterial culture supernatants by mass spectrometry analysis. In addition, this strain has multiple genes contributing to its plant growth-promoting profile, including those for acetoin, 2,3-butanediol and indole-3-acetic acid production, siderophores biosynthesis, and solubilisation of organic and inorganic phosphate. A pan-genomic analysis demonstrated that the genome of strain T16 possesses large gene clusters that are absent in the genomes of B. ambifaria reference strains. According to predictions, most of these clusters would be involved in aromatic compounds degradation. One genomic region, encoding flavin-dependent monooxygenases of unknown function, is proposed as a candidate responsible for fusaric acid degradation.

Multiplexed Promoter Engineering for Improving Thaxtomin A Production in Heterologous Streptomyces Hosts.

Thaxtomin A is a potent bioherbicide in both organic and conventional agriculture; however, its low yield hinders its wide application. Here, we report the direct cloning and heterologous expression of the thaxtomin A gene cluster in three well-characterized Streptomyces hosts. Then, we present an efficient, markerless and multiplex large gene cluster editing method based on in vitro CRISPR/Cas9 digestion and yeast homologous recombination. With this method, we successfully engineered the thaxtomin A cluster by simultaneously replacing the native promoters of the txtED operon, txtABH operon and txtC gene with strong constitutive promoters, and the yield of thaxtomin A improved to 289.5 µg/mL in heterologous Streptomyces coelicolor M1154. To further optimize the biosynthetic pathway, we used constraint-based combinatorial design to build 27 refactored gene clusters by varying the promoter strength of every operon, and the highest titer of thaxtomin A production reached 504.6 μg/mL. Taken altogether, this work puts forward a multiplexed promoter engineering strategy to engineer secondary metabolism gene clusters for efficiently improving fermentation titers.

Genetics Behind the Glycosylation Patterns in the Biosynthesis of Dalbaheptides.

Glycopeptide antibiotics are valuable natural metabolites endowed with different pharmacological properties, among them are dalbaheptides used to treat different infections caused by multidrug-resistant Gram-positive pathogens. Dalbaheptides are produced by soil-dwelling high G-C Gram-positive actinobacteria. Their biosynthetic pathways are encoded within large biosynthetic gene clusters. A non-ribosomally synthesized heptapeptide aglycone is the common scaffold for all dalbaheptides. Different enzymatic tailoring steps, including glycosylation, are further involved in decorating it. Glycosylation of dalbaheptides is a crucial step, conferring them specific biological activities. It is achieved by a plethora of glycosyltransferases, encoded within the corresponding biosynthetic gene clusters, able to install different sugar residues. These sugars might originate from the primary metabolism, or, alternatively, their biosynthesis might be encoded within the biosynthetic gene clusters. Already installed monosaccharides might be further enzymatically modified or work as substrates for additional glycosylation. In the current minireview, we cover recent updates concerning the genetics and enzymology behind the glycosylation of dalbaheptides, building a detailed and consecutive picture of this process and of its biological evolution. A thorough understanding of how glycosyltransferases function in dalbaheptide biosynthesis might open new ways to use them in chemo-enzymatic synthesis and/or in combinatorial biosynthesis for building novel glycosylated antibiotics.

Stonewall prevents expression of ectopic genes in the ovary and accumulates at insulator elements in D. melanogaster.

Germline stem cells (GSCs) are the progenitor cells of the germline for the lifetime of an animal. In Drosophila, these cells reside in a cellular niche that is required for both their maintenance (self-renewal) and differentiation (asymmetric division resulting in a daughter cell that differs from the GSC). The stem cell—daughter cell transition is tightly regulated by a number of processes, including an array of proteins required for genome stability. The germline stem-cell maintenance factor Stonewall (Stwl) associates with heterochromatin, but its molecular function is poorly understood. We performed RNA-Seq on stwl mutant ovaries and found significant derepression of many transposon families but not heterochromatic genes. We also discovered inappropriate expression of multiple classes of genes. Most prominent are testis-enriched genes, including the male germline sex-determination switch Phf7, the differentiation factor bgcn, and a large testis-specific gene cluster on chromosome 2, all of which are upregulated or ectopically expressed in stwl mutant ovaries. Surprisingly, we also found that RNAi knockdown of stwl in somatic S2 cells results in ectopic expression of these testis genes. Using parallel ChIP-Seq and RNA-Seq experiments in S2 cells, we discovered that Stwl localizes upstream of transcription start sites and at heterochromatic sequences including repetitive sequences associated with telomeres. Stwl is also enriched at bgcn, suggesting that it directly regulates this essential differentiation factor. Finally, we identify Stwl binding motifs that are shared with known insulator binding proteins. We propose that Stwl affects gene regulation, including repression of male transcripts in the female germline, by binding insulators and establishing chromatin boundaries.

Cre/lox-Mediated Chromosomal Integration of Biosynthetic Gene Clusters for Heterologous Expression in Aspergillus nidulans.

Building strains of filamentous fungi for stable long-term heterologous expression of large biosynthetic pathways is limited by the low transformation efficiency or genetic stability of current methods. Here, we developed a system for targeted chromosomal integration of large biosynthetic gene clusters in Aspergillus nidulans based on site-specific recombinase-mediated cassette exchange. We built A.nidulans strains harboring a chromosomal landing pad for Cre/lox-mediated recombination and demonstrated efficient targeted integration of a 21 kb DNA fragment in a single step. We further evaluated the integration at two loci by analyzing the expression of a fluorescent reporter and the production of a heterologous polyketide metabolite. We compared chromosomal expression at those landing loci to episomal AMA1-based expression, which also shed light on uncharacterized aspects of episomal expression in filamentous fungi. This is the first demonstration of site-specific recombinase-mediated integration in filamentous fungi, setting the foundations for the further development of this tool.

Accessory Genome Dynamics of Local and Global Staphylococcus pseudintermedius Populations.

Staphylococcus pseudintermedius is a major bacterial colonizer and opportunistic pathogen in dogs. Methicillin-resistant S. pseudintermedius (MRSP) continues to emerge as a significant challenge to maintaining canine health. We sought to determine the phylogenetic relationships of S. pseudintermedius across five states in the New England region of the United States and place them in a global context. The New England dataset consisted of 125 previously published S. pseudintermedius genomes supplemented with 45 newly sequenced isolates. The core genome phylogenetic tree revealed many deep branching lineages consisting of 142 multi-locus sequence types (STs). In silico detection of the mecA gene revealed 40 MRSP and 130 methicillin-susceptible S. pseudintermedius (MSSP) isolates. MRSP were derived from five structural types of SCCmec, the mobile genetic element that carries the mecA gene conferring methicillin resistance. Although many genomes were MSSP, they nevertheless harbored genes conferring resistance to many other antibiotic classes, including aminoglycosides, macrolides, tetracyclines and penams. We compared the New England genomes to 297 previously published genomes sampled from five other states in the United States and 13 other countries. Despite the prevalence of the clonally expanding ST71 found worldwide and in other parts of the United States, we did not detect it in New England. We next sought to interrogate the combined New England and global datasets for the presence of coincident gene pairs linked to antibiotic resistance. Analysis revealed a large co-circulating accessory gene cluster, which included mecA as well as eight other resistance genes [aac (6′)-Ie-aph (2″)-Ia, aad (6), aph (3′)-IIIa, sat4, ermB, cat, blaZ, and tetM]. Furthermore, MRSP isolates carried significantly more accessory genes than their MSSP counterparts. Our results provide important insights to the evolution and geographic spread of high-risk clones that can threaten the health of our canine companions.

A large transposable element mediates metal resistance in the fungus Paecilomyces variotii

The horizontal transfer of large gene clusters by mobile elements is a key driver of prokaryotic adaptation in response to environmental stresses. Eukaryotic microbes face similar stresses; however, a parallel role for mobile elements has not been established. A stress faced by many microorganisms is toxic metal ions in their environment. In fungi, identified mechanisms for protection against metals generally rely on genes that are dispersed within an organism's genome. Here, we discover a large (∼85 kb) region that confers tolerance to five metal/metalloid ions (arsenate, cadmium, copper, lead, and zinc) in the genomes of some, but not all, strains of a fungus, Paecilomyces variotii. We name this region HEPHAESTUS (Hφ) and present evidence that it is mobile within the P.variotii genome with features characteristic of a transposable element. HEPHAESTUS contains the greatest complement of host-beneficial genes carried by a transposable element in eukaryotes, suggesting that eukaryotic transposable elements might play a role analogous to bacteria in the horizontal transfer of large regions of host-beneficial DNA. Genes within HEPHAESTUS responsible for individual metal tolerances include those encoding a P-type ATPase transporter-PcaA-required for cadmium and lead tolerance, a transporter-ZrcA-providing tolerance to zinc, and a multicopper oxidase-McoA-conferring tolerance to copper. In addition, a subregion of Hφ confers tolerance to arsenate. The genome sequences of other fungi in the Eurotiales contain further examples of HEPHAESTUS, suggesting that it is responsible for independently assembling tolerance to a diverse array of ions, including chromium, mercury, and sodium.

Discovery of an Antarctic Ascidian-Associated Uncultivated Verrucomicrobia with Antimelanoma Palmerolide Biosynthetic Potential.

ABSTRACTThe Antarctic marine ecosystem harbors a wealth of biological and chemical innovation that has risen in concert over millennia since the isolation of the continent and formation of the Antarctic circumpolar current. Scientific inquiry into the novelty of marine natural products produced by Antarctic benthic invertebrates led to the discovery of a bioactive macrolide, palmerolide A, that has specific activity against melanoma and holds considerable promise as an anticancer therapeutic. While this compound was isolated from the Antarctic ascidian Synoicum adareanum, its biosynthesis has since been hypothesized to be microbially mediated, given structural similarities to microbially produced hybrid nonribosomal peptide-polyketide macrolides. Here, we describe a metagenome-enabled investigation aimed at identifying the biosynthetic gene cluster (BGC) and palmerolide A-producing organism. A 74-kbp candidate BGC encoding the multimodular enzymatic machinery (hybrid type I-trans-AT polyketide synthase-nonribosomal peptide synthetase and tailoring functional domains) was identified and found to harbor key features predicted as necessary for palmerolide A biosynthesis. Surveys of ascidian microbiome samples targeting the candidate BGC revealed a high correlation between palmerolide gene targets and a single 16S rRNA gene variant (R = 0.83 to 0.99). Through repeated rounds of metagenome sequencing followed by binning contigs into metagenome-assembled genomes, we were able to retrieve a nearly complete genome (10 contigs) of the BGC-producing organism, a novel verrucomicrobium within the Opitutaceae family that we propose here as “Candidatus Synoicihabitans palmerolidicus.” The refined genome assembly harbors five highly similar BGC copies, along with structural and functional features that shed light on the host-associated nature of this unique bacterium.IMPORTANCE Palmerolide A has potential as a chemotherapeutic agent to target melanoma. We interrogated the microbiome of the Antarctic ascidian, Synoicum adareanum, using a cultivation-independent high-throughput sequencing and bioinformatic strategy. The metagenome-encoded biosynthetic machinery predicted to produce palmerolide A was found to be associated with the genome of a member of the S. adareanum core microbiome. Phylogenomic analysis suggests the organism represents a new deeply branching genus, “Candidatus Synoicihabitans palmerolidicus,” in the Opitutaceae family of the Verrucomicrobia phylum. The Ca. Synoicihabitans palmerolidicus 4.29-Mb genome encodes a repertoire of carbohydrate-utilizing and transport pathways, a chemotaxis system, flagellar biosynthetic capacity, and other regulatory elements enabling its ascidian-associated lifestyle. The palmerolide producer’s genome also contains five distinct copies of the large palmerolide biosynthetic gene cluster that may provide structural complexity of palmerolide variants.