Enhancement of Carrimycin Production Via Traditional Mutagenesis with Metabolic Engineering in Streptomyces Spiramyceticus 54IA

Background: Carrimycin is a new approved class I antibiotic in China. The novel carrimycin producing strain, Streptomyces spiramyceticus 54IA, was constructed by CRISPR-Cas9 editing system without insertion of antibiotics resistant gene. The problem of low yield limits this strain in large scale fermentation. In this study, the carrimycin production was significantly improved by strain mutagenesis coupled metabolic engineering. Results: The sspD gene is responsible for degradation of triacylglycerol to provide precursors of the polyketide biosynthesis. The extra sspD gene controlled by the promoters of pks and bsm42 genes could moderately enhance carrimycin production. The Bsm42 was identified to play a pathway-specific positive regulator for carrimycin biosynthesis. Due to production of carrimycin significantly enhanced by bsm42 overexpression, the two different length promoters of bsm42 individually ligated with two reporter genes were used to monitor bsm42 expression for screening the higher carrimycin production mutants treated by plasma and ultraviolet. 47% of the 608 selected mutants had higher fermentation titer than the starting strain. The shorter promoter of bsm42 displayed more appropriate for selection of the carrimycin production improved mutants. The F2R-15 mutant had highest titer (1010±30 μg/mL), which was about 9 times higher than that of 54IA strain. Comparative analysis of transcriptome profiles of F2R-15 mutant and 54IA strains found 158 differential expression genes with more than 2 fold-changes. The up-regulated genes were associated with macrolide precursor biosynthesis, macrolide-inactivation, antibiotics transporter, oxidative phosphorylation; while the most down-regulated genes were referring to the primary metabolites synthetic genes and biosynthetic genes of other secondary metabolites. Conclusion: These results suggested that manipulation of the positive regulatory gene bsm42 and traditional mutagenesis coupled with reporter-guided mutant selection method facilitated selection of carrimycin high-yielding mutants.

New Insights into the Roles of Osmanthus Fragrans Heat-Shock Transcription Factors in Cold and Other Stress Responses

Sweet osmanthus (Osmanthus fragrans) is an evergreen woody plant that emits a floral aroma and is widely used in the landscape and fragrance industries. However, its application and cultivation regions are limited by cold stress. Heat-shock transcription factor (HSF) family members are widely present in plants and participate in, and regulate, the defense processes of plants under various abiotic stress conditions, but now, the role of this family in the responses of O. fragrans to cold stress is still not clear. Here, 46 OfHSF members were identified in the O. fragrans genome and divided into three subfamilies on the basis of a phylogenetic analysis. The promoter regions of most OfHSFs contained many cis-acting elements involved in multiple hormonal and abiotic stresses. RNA-seq data revealed that most of OfHSF genes were differentially expressed in various tissues, and some OfHSF members were induced by cold stress. The qRT-PCR analysis identified four OfHSFs that were induced by both cold and heat stresses, in which OfHSF11 and OfHSF43 had contrary expression trends under cold stress conditions and their expression patterns both showed recovery tendencies after the cold stress. OfHSF11 and OfHSF43 localized to the nuclei and their expression patterns were also induced under multiple abiotic stresses and hormonal treatments, indicating that they play critical roles in responses to multiple stresses. Furthermore, after a cold treatment, transient expression revealed that the malondialdehyde (MDA) content of OfHSF11-transformed tobacco significantly increased, and the expression levels of cold-response regulatory gene NbDREB3, cold response gene NbLEA5 and ROS detoxification gene NbCAT were significantly inhibited, implying that OfHSF11 is a negative regulator of cold responses in O. fragrans. Our study contributes to the further functional characterization of OfHSFs and will be useful in developing improved cold-tolerant cultivars of O. fragrans.

Identification and Evolutionary Characteristic Analysis of STARD Gene Family, and Overexpression VvSTARD5 Responses to Salt Stress in Tomato

This study aimed to have a full understanding of the steroidogenic acute regulatory gene family member and evolutionary relationship in grape. 23 VvSTARD gene members were identified and divided into five groups in different species. Analyses of the gene codon preference, selective pressure, and tandem duplication of the VvSTARD, AtSTARD, and OsSTARD genes indicated that synteny relationship occurred in grapes, Arabidopsis thaliana, and rice genomes. The 8 lipid transporter proteins were found in the tertiary structure of the STARD gene family in grape. Expression profiles of the three species microarrays showed that the expression levels of the STARD genes in different organs and the response to abiotic stress in the same subgroup had similar characteristics. In addition, analysis of the VvSTARD genes expression levels was detected in response to different hormones and abiotic stresses by quantitative real-time polymerase chain reaction (qRT-PCR), and the results were the same as those predicted by the cis-elements and the expression profiles. Meanwhile, VvSTARD5 gene was screened in high concentration NaCl treatment by qRT-PCR. Furthermore, the VvSTARD5 was located at the nucleus by subcellular location. Through the function analysis of salt tolerance in transgenic tomato, overexpression VvSTARD5 obviously improved tolerance to salt stress. Taken together, our findings Preliminary identify the functions of VvSTARD gene family and vertify STARD5 that be likely involved in regulating salt tolerance, which may have potential application molecular breeding in grape.

Genetic Engineering of Streptomyces ghanaensis ATCC14672 for Improved Production of Moenomycins

Streptomycetes are soil-dwelling multicellular microorganisms famous for their unprecedented ability to synthesize numerous bioactive natural products (NPs). In addition to their rich arsenal of secondary metabolites, Streptomyces are characterized by complex morphological differentiation. Mostly, industrial production of NPs is done by submerged fermentation, where streptomycetes grow as a vegetative mycelium forming pellets. Often, suboptimal growth peculiarities are the major bottleneck for industrial exploitation. In this work, we employed genetic engineering approaches to improve the production of moenomycins (Mm) in Streptomyces ghanaensis, the only known natural direct inhibitors of bacterial peptidoglycan glycosyltransferses. We showed that in vivo elimination of binding sites for the pleiotropic regulator AdpA in the oriC region strongly influences growth and positively correlates with Mm accumulation. Additionally, a marker- and “scar”-less deletion of moeH5, encoding an amidotransferase from the Mm gene cluster, significantly narrows down the Mm production spectrum. Strikingly, antibiotic titers were strongly enhanced by the elimination of the pleiotropic regulatory gene wblA, involved in the late steps of morphogenesis. Altogether, we generated Mm overproducers with optimized growth parameters, which are useful for further genome engineering and chemoenzymatic generation of novel Mm derivatives. Analogously, such a scheme can be applied to other Streptomyces spp.

The MocR/GabR Ectoine and Hydroxyectoine Catabolism Regulator EnuR: Inducer and DNA Binding

The compatible solutes ectoine and 5-hydroxyectoine are widely synthesized by bacteria as osmostress protectants. These nitrogen-rich tetrahydropyrimidines can also be exploited as nutrients by microorganisms. Many ectoine/5-hydroxyectoine catabolic gene clusters are associated with a regulatory gene (enuR: ectoine nutrient utilization regulator) encoding a repressor protein belonging to the MocR/GabR sub-family of GntR-type transcription factors. Focusing on EnuR from the marine bacterium Ruegeria pomeroyi, we show that the dimerization of EnuR is mediated by its aminotransferase domain. This domain can fold independently from its amino-terminal DNA reading head and can incorporate pyridoxal-5′-phosphate (PLP) as cofactor. The covalent attachment of PLP to residue Lys302 of EnuR was proven by mass-spectrometry. PLP interacts with system-specific, ectoine and 5-hydroxyectoine-derived inducers: alpha-acetyldiaminobutyric acid (alpha-ADABA), and hydroxy-alpha-acetyldiaminobutyric acid (hydroxy-alpha-ADABA), respectively. These inducers are generated in cells actively growing with ectoines as sole carbon and nitrogen sources, by the EutD hydrolase and targeted metabolic analysis allowed their detection. EnuR binds these effector molecules with affinities in the low micro-molar range. Studies addressing the evolutionary conservation of EnuR, modelling of the EnuR structure, and docking experiments with the inducers provide an initial view into the cofactor and effector binding cavity. In this cavity, the two high-affinity inducers for EnuR, alpha-ADABA and hydroxy-alpha-ADABA, are positioned such that their respective primary nitrogen group can chemically interact with PLP. Purified EnuR bound with micro-molar affinity to a 48 base pair DNA fragment containing the sigma-70 type substrate-inducible promoter for the ectoine/5-hydroxyectoine importer and catabolic gene cluster. Consistent with the function of EnuR as a repressor, the core elements of the promoter overlap with two predicted EnuR operators. Our data lend themselves to a straightforward regulatory model for the initial encounter of EnuR-possessing ectoine/5-hydroxyectoine consumers with environmental ectoines and for the situation when the external supply of these compounds has been exhausted by catabolism.

Statistical estimates of multiple transcription factors binding in the model plant genomes based on ChIP-seq data

The development of high-throughput genomic sequencing coupled with chromatin immunoprecipitation technologies allows studying the binding sites of the protein transcription factors (TF) in the genome scale. The growth of data volume on the experimentally determined binding sites raises qualitatively new problems for the analysis of gene expression regulation, prediction of transcription factors target genes, and regulatory gene networks reconstruction. Genome regulation remains an insufficiently studied though plants have complex molecular regulatory mechanisms of gene expression and response to environmental stresses. It is important to develop new software tools for the analysis of the TF binding sites location and their clustering in the plant genomes, visualization, and the following statistical estimates. This study presents application of the analysis of multiple TF binding profiles in three evolutionarily distant model plant organisms. The construction and analysis of non-random ChIP-seq binding clusters of the different TFs in mammalian embryonic stem cells were discussed earlier using similar bioinformatics approaches. Such clusters of TF binding sites may indicate the gene regulatory regions, enhancers and gene transcription regulatory hubs. It can be used for analysis of the gene promoters as well as a background for transcription networks reconstruction. We discuss the statistical estimates of the TF binding sites clusters in the model plant genomes. The distributions of the number of different TFs per binding cluster follow same power law distribution for all the genomes studied. The binding clusters in Arabidopsis thaliana genome were discussed here in detail.

Runx2 rs59983488 polymorphism in class II malocclusion in the Indonesian subpopulation

Background: Class II malocclusion is one of the main orthodontic issues for patients in seeking treatment. The prevalence of class II malocclusion varies in different populations. Variation in skeletal profile is mainly controlled internally by a regulatory gene. Runt-related transcription factor-2 (Runx2) plays a role in osteoblast differentiation and is highly expressed during development. Purpose: This study aimed to evaluate the relation of regulatory gene variation in the Runx2 promoter with class II malocclusion. Methods: DNA samples were acquired from 95 orthodontic patients in Jakarta, Indonesia, who were divided into two groups: class I skeletal malocclusion (control group) and class II malocclusion. A single nucleotide polymorphism was investigated using the polymerase chain reaction and restriction fragment length polymorphism techniques. The distribution of alleles was assessed using the Hardy-Weinberg test. The relationship between polymorphism and skeletal variation was assessed with the Chi-Square test and logistic regression. Results: The frequency distributions of genotypes and alleles were tested for Hardy-Weinberg equilibrium and found to be slightly deviated. There was an equal distribution of G and T alleles throughout class II and class I skeletal malocclusions and the Chi-Square test showed that this relationship was not significant (p=0.5). Conclusion: Runx2 rs59983488 polymorphism was found in the Indonesian subpopulation; however, an association between Runx2 rs59983488 polymorphism and class II skeletal malocclusion was not found.

Identification and characterization of BrxR as a regulatory gene in the BREX phage restriction system

Bacteriophage exclusion (BREX) phage restriction systems are found in a wide range of bacteria. Various BREX systems encode unique combinations of proteins that usually include a site-specific methyltransferase; none appear to contain a nuclease. Here we describe the identification and characterization of a Type I BREX system from Acinetobacter and the effect of deleting each BREX ORF on growth, methylation and phage restriction. The analysis identified a previously uncharacterized gene at the 5-prime end of the BREX operon that is dispensable for methylation but involved in restriction. Biochemical and crystallographic analyses of this factor, which we term BrxR (BREX Regulator), demonstrate that it forms a homodimer and specifically binds a pseudo-palindromic DNA target site upstream of its transcription start site. Precise deletion of the BrxR gene causes cell toxicity, reduces phage restriction, and significantly increases the expression of BrxC. In contrast, the introduction of a premature stop codon into the BrxR gene has little effect, implying that the BrxR coding sequence and BrxR protein have independent functional roles in BREX regulation. We speculate that the BrxR coding sequence is involved in cis regulation of BREX activity and that the BrxR protein may play an additional regulatory role, perhaps during horizontal transfer of the system.

Lung Adenocarcinoma Transcriptomic Analysis Predicts Adenylate Kinase Signatures Contributing to Tumor Progression and Negative Patient Prognosis

The ability to detect and respond to hypoxia within a developing tumor appears to be a common feature amongst most cancers. This hypoxic response has many molecular drivers, but none as widely studied as Hypoxia-Inducible Factor 1 (HIF-1). Recent evidence suggests that HIF-1 biology within lung adenocarcinoma (LUAD) may be associated with expression levels of adenylate kinases (AKs). Using LUAD patient transcriptome data, we sought to characterize AK gene signatures related to lung cancer hallmarks, such as hypoxia and metabolic reprogramming, to identify conserved biological themes across LUAD tumor progression. Transcriptomic analysis revealed perturbation of HIF-1 targets to correlate with altered expression of most AKs, with AK4 having the strongest correlation. Enrichment analysis of LUAD tumor AK4 gene signatures predicts signatures involved in pyrimidine, and by extension, nucleotide metabolism across all LUAD tumor stages. To further discriminate potential drivers of LUAD tumor progression within AK4 gene signatures, partial least squares discriminant analysis was used at LUAD stage-stage interfaces, identifying candidate genes that may promote LUAD tumor growth or regression. Collectively, these results characterize regulatory gene networks associated with the expression of all nine human AKs that may contribute to underlying metabolic perturbations within LUAD and reveal potential mechanistic insight into the complementary role of AK4 in LUAD tumor development.

Prognostic Risk Model and Tumor Immune Environment Modulation of m5C-Related LncRNAs in Pancreatic Ductal Adenocarcinoma

RNA methylation modification is a key process in epigenetics that regulates posttranscriptional gene expression. With advances in next-generation sequencing technology, 5-methylcytosine (m5C) modification has also been found in multiple RNAs. Long non-coding RNAs (lncRNAs) were proved to have a key role in cancer progression and closely related to the tumor immune microenvironment. Thus, based on the PDAC patients’ clinical information and genetic transcriptome data from the TCGA database, we performed a detailed bioinformatic analysis to establish a m5C-related lncRNA prognostic risk model for PDAC patients and discovered the relationship between the risk model and PDAC immune microenvironment. Pearson correlation coefficient analysis was applied to conduct a m5C regulatory gene and m5C-related lncRNA co-expression network. Expression of m5C-related lncRNAs screened by univariate regression analysis with prognostic value showed a significant difference between pancreatic cancer and normal tissues. The least absolute shrinkage and selection operator (LASSO) Cox regression method was applied to determine an 8-m5C-related lncRNA prognostic risk model. We used principal component analysis to indicate that the risk model could distinguish all the samples clearly. The clinical nomogram also accurately predicted 1-, 1.5-, 2-, and 3-year survival time among PDAC patients. Additionally, this risk model was validated in the entire group and sub-test groups using KM analysis and ROC analysis. Combined with the clinical characteristics, the risk score was found to be an independent factor for predicting the survival of PDAC patients. Furthermore, the association between the risk model and tumor immune microenvironment was evaluated via the ESTIMATE R package and CIBERSORT method. Consequently, the results indicated that immune cells were associated with m5C-related lncRNA risk model scores and had different distribution in the high- and low-risk groups. Based on all these analyses, the m5C-related lncRNA risk model could be a reliable prognostic tool and therapeutic target for PDAC patients.