Clavulanic acid biosynthesis and genetic manipulation for its overproduction

Clavulanic acid, a β-lactamase inhibitor, is used together with β-lactam antibiotics to create drug mixtures possessing potent antimicrobial activity. In view of the clinical and industrial importance of clavulanic acid, identification of the clavulanic acid biosynthetic pathway and the associated gene cluster(s) in the main producer species, Streptomyces clavuligerus, has been an intriguing research question. Clavulanic acid biosynthesis was revealed to involve an interesting mechanism common to all of the clavam metabolites produced by the organism, but different from that of other β-lactam compounds. Gene clusters involved in clavulanic acid biosynthesis in S. clavuligerus occupy large regions of nucleotide sequence in three loci of its genome. In this review, clavulanic acid biosynthesis and the associated gene clusters are discussed, and clavulanic acid improvement through genetic manipulation is explained.     

Complete mitochondrial genome of the black‐tailed hornet,Vespa ducalis (Hymenoptera: Vespidae): Genomic comparisons in Vespoidea

We sequenced the complete mitochondrial genome (mitogenome) of the black‐tailed hornet, Vespa ducalis (Hymenoptera: Vespidae). The genome was 15,779‐bp long and contained typical sets of genes [13 protein‐coding genes (PCGs), 22 tRNAs, and 2 rRNAs]. The V. ducalis A + T‐rich region was 166‐bp long and was the shortest of all sequenced Vespoidea genomes, including Vespa. The genome was highly biased toward A/T nucleotides—80.1 % in the whole genome, 77.8 % in PCGs, 83.4–85.6 % in RNAs, and 92.8 % in the A + T‐rich region. These values are well within the typical range for genes and regions of Vespoidea mitogenomes. Start and stop codons in several Vespa species—including V. ducalis—were diversified, despite these species belonging to the same genus. In comparison with the ancestral mitogenomes, Vespa mitogenomes—including that of V. ducalis—showed substantial gene rearrangement; however, we detected no gene rearrangement among Vespa species. We conducted phylogenetic reconstruction based on concatenated sequences of 13 PCGs and two rRNAs (12,755 bp ) in available species of Vespoidea—21 species in six subfamilies in two families (Vespidae and Formicidae). The Bayesian inference and maximum likelihood (ML) methods revealed that each family formed strong monophyletic groups [Bayesian posterior probability (BPP) = 1; ML, 100 %]. Moreover, V. ducalis and V. mandarinia formed a strong sister group (BPP = 1; ML, 94 %).     

Glycerol Utilization Gene Cluster in Streptomyces clavuligerus

The Streptomyces clavuligerus ATCC 27064 glycerol cluster gylR-glpF1K1D1 is induced by glycerol but is not affected by glucose. S. clavuligerus growth and clavulanic acid production are stimulated by glycerol, but this does not occur in a glpK1-deleted mutant. Amplification of glpK1D1 results in transformants yielding larger amounts of clavulanic acid in the wild-type strain and in overproducer S. clavuligerus Gap15-7-30 or S. clavuligerus ΔrelA strains.Streptomyces clavuligerus ATCC 27064 produces the β-lactamase inhibitor clavulanic acid (CA). This compound is formed from arginine (17) and the three-carbon molecule glyceraldehyde-3-phosphate (6) which are condensed by the carboxyethylarginine synthase, the first enzyme of the pathway, encoded by ceaS2. Mutants disrupted in this gene do not produce CA in tryptic soy broth or starch-asparagine (SA) medium but form moderate amounts of CA in glycerol-supplemented media, probably due to glycerol utilization through the duplicated CeaS1 carboxyethylarginine synthase (10).The role of d-glyceraldehyde-3-phosphate as a CA precursor was further supported by the construction of a glyceraldehyde-3-phosphate dehydrogenase (gap1) mutant of S. clavuligerus, which produces 180 to 210% CA in comparison to the wild-type strain due to higher availability of the glyceraldehyde-3-phosphate precursor (9).Genes for glycerol utilization in Streptomyces coelicolor form an operon, gylCABX (15, 16), containing a gene for a putative glycerol transporter, a glycerol kinase, a glycerol-3-phosphate dehydrogenase, and a gene of unknown function. They are preceded by gylR (5), which encodes a glycerol-inducible repressor controlling both gylR and the gylCABX operon. Glycerol induction and glucose catabolite repression of the glp genes are thought to be directly related to the GylR protein in S. coelicolor (5).Due to the importance of the glycerol flow for CA production, we decided to analyze the glycerol-utilizing gene cluster of S. clavuligerus.     

The genome sequence of Samia ricini,a new model species of lepidopteran insect

Samia ricini, a gigantic saturniid moth, has the potential to be a novel lepidopteran model species. Samia ricini is far more resistant to diseases than the current model species Bombyx mori, and therefore can be more easily reared. In addition, genetic resources available for S. ricini rival those for B. mori: at least 26 ecoraces of S. ricini are reported and S. ricini can hybridize with wild Samia species, which are distributed throughout Asian countries, and produce fertile progenies. Physiological traits such as food preference, integument colour and larval spot pattern differ among S. ricini strains and wild Samia species so that those traits can be targeted in forward genetic analyses. To facilitate genetic research in S. ricini, we determined its whole genome sequence. The assembled genome of S. ricini was 458 Mb with 155 scaffolds, and the scaffold N50 length of the assembly was ~ 21 Mb. In total, 16,702 protein coding genes were predicted. While the S. ricini genome was mostly collinear with that of B. mori with some rearrangements and few S. ricini‐specific genes were discovered, chorion genes and fibroin genes seemed to have expanded in the S. ricini lineage. As the first step of genetic analyses, causal genes for “Blue,” “Yellow,” “Spot,” and “Red cocoon” phenotypes were mapped to chromosomes.     

Chromosomal‐level genomes of three rice planthoppers provide new insights into sex chromosome evolution

The brown planthopper Nilaparvata lugens, white‐backed planthopper Sogatella furcifera, and small brown planthopper Laodelphax striatellus are three major insect pests of rice. They are genetically close; however, they differ in several ecological traits such as host range, migration capacity, and in their sex chromosomes. Though the draft genome of these three planthoppers have been previously released, the quality of genome assemblies need to be improved. The absence of chromosome‐level genome resources has hindered in‐depth research of these three species. Here, we performed a de novo genome assembly for N. lugens to increase its genome assembly quality with PacBio and Illumina platforms, increasing the contig N50 to 589.46 Kb. Then, with the new N. lugens genome and previously reported S. furcifera and L. striatellus genome assemblies, we generated chromosome‐level scaffold assemblies of these three planthopper species using HiC scaffolding technique. The scaffold N50s significantly increased to 77.63 Mb, 43.36 Mb and 29.24 Mb for N. lugens, S. furcifera and L. striatellus, respectively. To identify sex chromosomes of these three planthopper species, we carried out genome re‐sequencing of males and females and successfully determined the X and Y chromosomes for N. lugens, and X chromosome for S. furcifera and L. striatellus. The gene content of the sex chromosomes showed high diversity among these three planthoppers suggesting the rapid evolution of sex‐linked genes, and all chromosomes showed high synteny. The chromosome‐level genome assemblies of three planthoppers would provide a valuable resource for a broad range of future research in molecular ecology, and subsequently benefits development of modern pest control strategies.     

Chromosome‐level de novo genome assembly of Sarcophaga peregrina provides insights into the evolutionary adaptation of flesh flies

Sarcophaga peregrina is considered to be of great ecological, medical and forensic significance, and has unusual biological characteristics such as an ovoviviparous reproductive pattern and adaptation to feed on carrion. The availability of a high‐quality genome will help to further reveal the mechanisms underlying these charcateristics. Here we present a de novo‐assembled genome at chromosome scale for S. peregrina. The final assembled genome was 560.31 Mb with contig N50 of 3.84 Mb. Hi‐C scaffolding reliably anchored six pseudochromosomes, accounting for 97.76% of the assembled genome. Moreover, 45.70% of repeat elements were identified in the genome. A total of 14,476 protein‐coding genes were functionally annotated, accounting for 92.14% of all predicted genes. Phylogenetic analysis indicated that S. peregrina and S. bullata diverged ~ 7.14 million years ago. Comparative genomic analysis revealed expanded and positively selected genes related to biological features that aid in clarifying its ovoviviparous reproduction and carrion‐feeding adaptations, such as lipid metabolism, olfactory receptor activity, antioxidant enzymes, proteolysis and serine‐type endopeptidase activity. Protein‐coding genes associated with ovoviparity, such as yolk proteins, transferrin and acid sphingomyelinase, were identified. This study provides a valuable genomic resource for S. peregrina, and sheds insight into further revealing the underlying molecular mechanisms of adaptive evolution.     

The complete chloroplast genome of Gracilariopsis lemaneiformis (Rhodophyta) gives new insight into the evolution of family Gracilariaceae

The complete chloroplast genome of Gracilariopsis lemaneiformis was recovered from a Next Generation Sequencing data set. Without quadripartite structure, this chloroplast genome (183,013 bp, 27.40% GC content) contains 202 protein‐coding genes, 34 tRNA genes, 3 rRNA genes, and 1 tmRNA gene. Synteny analysis showed plasmid incorporation regions in chloroplast genomes of three species of family Gracilariaceae and in Grateloupia taiwanensis of family Halymeniaceae. Combined with reported red algal plasmid sequences in nuclear and mitochondrial genomes, we postulated that red algal plasmids may have played an important role in ancient horizontal gene transfer among nuclear, chloroplast, and mitochondrial genomes. Substitution rate analysis showed that purifying selective forces maintaining stability of protein‐coding genes of nine red algal chloroplast genomes over long periods must be strong and that the forces acting on gene groups and single genes of nine red algal chloroplast genomes were similar and consistent. The divergence of Gp. lemaneiformis occurred ~447.98 million years ago (Mya), close to the divergence time of genus Pyropia and Porphyra (443.62 Mya).     

Genome sequence of M6, a diploid inbred clone of the high‐glycoalkaloid‐producing tuber‐bearing potato species Solanum chacoense,reveals residual heterozygosity

Cultivated potato (Solanum tuberosum L.) is a highly heterozygous autotetraploid that presents challenges in genome analyses and breeding. Wild potato species serve as a resource for the introgression of important agronomic traits into cultivated potato. One key species is Solanum chacoense and the diploid, inbred clone M6, which is self‐compatible and has desirable tuber market quality and disease resistance traits. Sequencing and assembly of the genome of the M6 clone of S. chacoense generated an assembly of 825 767 562 bp in 8260 scaffolds with an N50 scaffold size of 713 602 bp. Pseudomolecule construction anchored 508 Mb of the genome assembly into 12 chromosomes. Genome annotation yielded 49 124 high‐confidence gene models representing 37 740 genes. Comparative analyses of the M6 genome with six other Solanaceae species revealed a core set of 158 367 Solanaceae genes and 1897 genes unique to three potato species. Analysis of single nucleotide polymorphisms across the M6 genome revealed enhanced residual heterozygosity on chromosomes 4, 8 and 9 relative to the other chromosomes. Access to the M6 genome provides a resource for identification of key genes for important agronomic traits and aids in genome‐enabled development of inbred diploid potatoes with the potential to accelerate potato breeding.     

Characterization of maltotriose transporters from the Saccharomyces eubayanus subgenome of the hybrid Saccharomyces pastorianus lager brewing yeast strain Weihenstephan 34/70

The genome from the Saccharomyces pastorianus industrial lager brewing strain Weihenstephan 34/70, a natural Saccharomyces cerevisiae/Saccharomyces eubayanus hybrid, indicated the presence of two different maltotriose transporter genes: a new gene in the S. eubayanus subgenome with 81% of homology to the AGT1 permease from S. cerevisiae, and an amplification of the S. eubayanus MTY1 maltotriose permease previously identified in S. pastorianus yeasts. To characterize these S. eubayanus transporter genes, we used a S. cerevisiae strain deleted in the AGT1 permease and introduced the desired permease gene(s) into this locus through homologous recombination. Our results indicate that both the MTY1 and AGT1 genes from the S. eubayanus subgenome encode functional maltotriose transporters that allow fermentation of this sugar by yeast cells, despite their apparent differences in the kinetics of maltotriose‐H⁺ symport activity. The presence of two maltotriose transporters in the S. eubayanus subgenome not only highlights the importance of sugar transport for efficient maltotriose utilization by industrial yeasts, but these new genes can be used in breeding and/or selection programs aimed at increasing yeast fitness for the efficient fermentation of brewer's wort.     

Assembly and annotation of a draft genome sequence for Glycine latifolia,a perennial wild relative of soybean

Glycine latifolia (Benth.) Newell & Hymowitz (2n = 40), one of the 27 wild perennial relatives of soybean, possesses genetic diversity and agronomically favorable traits that are lacking in soybean. Here, we report the 939‐Mb draft genome assembly of G. latifolia (PI 559298) using exclusively linked‐reads sequenced from a single Chromium library. We organized scaffolds into 20 chromosome‐scale pseudomolecules utilizing two genetic maps and the Glycine max (L.) Merr. genome sequence. High copy numbers of putative 91‐bp centromere‐specific tandem repeats were observed in consecutive blocks within predicted pericentromeric regions on several pseudomolecules. No 92‐bp putative centromeric repeats, which are abundant in G. max, were detected in G. latifolia or Glycine tomentella. Annotation of the assembled genome and subsequent filtering yielded a high confidence gene set of 54 475 protein‐coding loci. In comparative analysis with five legume species, genes related to defense responses were significantly overrepresented in Glycine‐specific orthologous gene families. A total of 304 putative nucleotide‐binding site (NBS)‐leucine‐rich‐repeat (LRR) genes were identified in this genome assembly. Different from other legume species, we observed a scarcity of TIR‐NBS‐LRR genes in G. latifolia. The G. latifolia genome was also predicted to contain genes encoding 367 LRR‐receptor‐like kinases, a family of proteins involved in basal defense responses and responses to abiotic stress. The genome sequence and annotation of G. latifolia provides a valuable source of alternative alleles and novel genes to facilitate soybean improvement. This study also highlights the efficacy and cost‐effectiveness of the application of Chromium linked‐reads in diploid plant genome de novo assembly.     

Environmental systems biology of cold‐tolerant phenotype in Saccharomyces species adapted to grow at different temperatures

Temperature is one of the leading factors that drive adaptation of organisms and ecosystems. Remarkably, many closely related species share the same habitat because of their different temporal or micro‐spatial thermal adaptation. In this study, we seek to find the underlying molecular mechanisms of the cold‐tolerant phenotype of closely related yeast species adapted to grow at different temperatures, namely S. kudriavzevii CA111 (cryo‐tolerant) and S. cerevisiae 96.2 (thermo‐tolerant). Using two different systems approaches, i. thermodynamic‐based analysis of a genome‐scale metabolic model of S. cerevisiae and ii. large‐scale competition experiment of the yeast heterozygote mutant collection, genes and pathways important for the growth at low temperature were identified. In particular, defects in lipid metabolism, oxidoreductase and vitamin pathways affected yeast fitness at cold. Combining the data from both studies, a list of candidate genes was generated and mutants for two predicted cold‐favouring genes, GUT2 and ADH3, were created in two natural isolates. Compared with the parental strains, these mutants showed lower fitness at cold temperatures, with S. kudriavzevii displaying the strongest defect. Strikingly, in S. kudriavzevii, these mutations also significantly improve the growth at warm temperatures. In addition, overexpression of ADH3 in S. cerevisiae increased its fitness at cold. These results suggest that temperature‐induced redox imbalances could be compensated by increased glycerol accumulation or production of cytosolic acetaldehyde through the deletion of GUT2 or ADH3, respectively.     

Evolutionary trajectory of phytoalexin biosynthetic gene clusters in rice

Plants frequently possess operon‐like gene clusters for specialized metabolism. Cultivated rice, Oryza sativa, produces antimicrobial diterpene phytoalexins represented by phytocassanes and momilactones, and the majority of their biosynthetic genes are clustered on chromosomes 2 and 4, respectively. These labdane‐related diterpene phytoalexins are biosynthesized from geranylgeranyl diphosphate via ent‐copalyl diphosphate or syn‐copalyl diphosphate. The two gene clusters consist of genes encoding diterpene synthases and chemical‐modification enzymes including P450s. In contrast, genes for the biosynthesis of gibberellins, which are labdane‐related phytohormones, are scattered throughout the rice genome similar to other plant genomes. The mechanism of operon‐like gene cluster formation remains undefined despite previous studies in other plant species. Here we show an evolutionary insight into the rice gene clusters by a comparison with wild Oryza species. Comparative genomics and biochemical studies using wild rice species from the AA genome lineage, including Oryza barthii, Oryza glumaepatula, Oryza meridionalis and the progenitor of Asian cultivated rice Oryza rufipogon indicate that gene clustering for biosynthesis of momilactones and phytocassanes had already been accomplished before the domestication of rice. Similar studies using the species Oryza punctata from the BB genome lineage, the distant FF genome lineage species Oryza brachyantha and an outgroup species Leersia perrieri suggest that the phytocassane biosynthetic gene cluster was present in the common ancestor of the Oryza species despite the different locations, directions and numbers of their member genes. However, the momilactone biosynthetic gene cluster evolved within Oryza before the divergence of the BB genome via assembly of ancestral genes.     

Resistance gene enrichment sequencing (RenSeq) enables reannotation of the NB‐LRR gene family from sequenced plant genomes and rapid mapping of resistance loci in segregating populations

RenSeq is a NB‐LRR (nucleotide binding‐site leucine‐rich repeat) gene‐targeted, Resistance gene enrichment and sequencing method that enables discovery and annotation of pathogen resistance gene family members in plant genome sequences. We successfully applied RenSeq to the sequenced potato Solanum tuberosum clone DM, and increased the number of identified NB‐LRRs from 438 to 755. The majority of these identified R gene loci reside in poorly or previously unannotated regions of the genome. Sequence and positional details on the 12 chromosomes have been established for 704 NB‐LRRs and can be accessed through a genome browser that we provide. We compared these NB‐LRR genes and the corresponding oligonucleotide baits with the highest sequence similarity and demonstrated that ~80% sequence identity is sufficient for enrichment. Analysis of the sequenced tomato S. lycopersicum ‘Heinz 1706’ extended the NB‐LRR complement to 394 loci. We further describe a methodology that applies RenSeq to rapidly identify molecular markers that co‐segregate with a pathogen resistance trait of interest. In two independent segregating populations involving the wild Solanum species S. berthaultii (Rpi‐ber2) and S. ruiz‐ceballosii (Rpi‐rzc1), we were able to apply RenSeq successfully to identify markers that co‐segregate with resistance towards the late blight pathogen Phytophthora infestans. These SNP identification workflows were designed as easy‐to‐adapt Galaxy pipelines.     

A genomic island provides <Emphasis Type="Italic">Acidithiobacillus ferrooxidans</Emphasis> ATCC 53993 additional copper resistance: a possible competitive advantage

There is great interest in understanding how extremophilic biomining bacteria adapt to exceptionally high copper concentrations in their environment. Acidithiobacillus ferrooxidans ATCC 53993 genome possesses the same copper resistance determinants as strain ATCC 23270. However, the former strain contains in its genome a 160-kb genomic island (GI), which is absent in ATCC 23270. This GI contains, amongst other genes, several genes coding for an additional putative copper ATPase and a Cus system. A. ferrooxidans ATCC 53993 showed a much higher resistance to CuSO₄ (>100 mM) than that of strain ATCC 23270 (<25 mM). When a similar number of bacteria from each strain were mixed and allowed to grow in the absence of copper, their respective final numbers remained approximately equal. However, in the presence of copper, there was a clear overgrowth of strain ATCC 53993 compared to ATCC 23270. This behavior is most likely explained by the presence of the additional copper-resistance genes in the GI of strain ATCC 53993. As determined by qRT-PCR, it was demonstrated that these genes are upregulated when A. ferrooxidans ATCC 53993 is grown in the presence of copper and were shown to be functional when expressed in copper-sensitive Escherichia coli mutants. Thus, the reason for resistance to copper of two strains of the same acidophilic microorganism could be determined by slight differences in their genomes, which may not only lead to changes in their capacities to adapt to their environment, but may also help to select the more fit microorganisms for industrial biomining operations.     

Flip‐flop organization in the chloroplast genome of Capsosiphon fulvescens (Ulvophyceae,Chlorophyta)

To better understand organelle genome evolution of the ulvophycean green alga Capsosiphon fulvescens, we sequenced and characterized its complete chloroplast genome. The circular chloroplast genome was 111,561 bp in length with 31.3% GC content that contained 108 genes including 77 protein‐coding genes, two copies of rRNA operons, and 27 tRNAs. In this analysis, we found the two types of isoform, called heteroplasmy, were likely caused by a flip‐flop organization. The flip‐flop mechanism may have caused structural variation and gene conversion in the chloroplast genome of C. fulvescens. In a phylogenetic analysis based on all available ulvophycean chloroplast genome data, including a new C. fulvescens genome, we found three major conflicting signals for C. fulvescens and its sister taxon Pseudoneochloris marina within 70 individual genes: (i) monophyly with Ulotrichales, (ii) monophyly with Ulvales, and (iii) monophyly with the clade of Ulotrichales and Ulvales. Although the 70‐gene concatenated phylogeny supported monophyly with Ulvales for both species, these complex phylogenetic signals of individual genes need further investigations using a data‐rich approach (i.e., organelle genome data) from broader taxon sampling.     

Development of genomic resources for Nothofagus species using next‐generation sequencing data

Using next‐generation sequencing, we developed the first whole‐genome resources for two hybridizing Nothofagus species of the Patagonian forests that crucially lack genomic data, despite their ecological and industrial value. A de novo assembly strategy combining base quality control and optimization of the putative chloroplast gene map yielded ~32 000 contigs from 43% of the reads produced. With 12.5% of assembled reads, we covered ~96% of the chloroplast genome and ~70% of the mitochondrial gene content, providing functional and structural annotations for 112 and 52 genes, respectively. Functional annotation was possible on 15% of the contigs, with ~1750 potentially novel nuclear genes identified for Nothofagus species. We estimated that the new resources (13.41 Mb in total) included ~4000 gene regions representing ~6.5% of the expected genic partition of the genome, the remaining contigs potentially being nongenic DNA. A high‐quality single nucleotide polymorphisms resource was developed by comparing various filtering methods, and preliminary results indicate a strong conservation of cpDNA genomes in contrast to numerous exclusive nuclear polymorphisms in both species. Finally, we characterized 2274 potential simple sequence repeat (SSR) loci, designed primers for 769 of them and validated nine of 29 loci in 42 individuals per species. Nothofagus obliqua had more alleles (4.89) on average than N. nervosa (2.89), 8 SSRs were efficient to discriminate species, and three were successfully transferred in three other Nothofagus species. These resources will greatly help for future inferences of demographic, adaptive and hybridizing events in Nothofagus species, and for conserving and managing natural populations.