Heterologous production of a new lanthipeptide boletupeptin using a cryptic biosynthetic gene cluster of the myxobacterium Melittangium boletus

Myxobacteria have comparatively large genomes that contain many biosynthetic genes with the potential to produce secondary metabolites. Based on genome mining, we discovered a new biosynthetic gene cluster of class III lanthipeptide in the genome of the myxobacterium Melittangium boletus. The biosynthetic gene cluster contained a precursor peptide-coding gene bolA, and a class III lanthipeptide synthetase-coding gene bolKC. The expression vector containing bolA and bolKC was constructed using synthetic DNA with codon-optimized sequences based on the commercially available vector pET29b. Co-expression of the two genes in the host Escherichia coli BL21(DE3) yielded a new class III lanthipeptide named boletupeptin. The structure of boletupeptin was proposed to have one unit of labionin, as determined by mass spectrometry experiments after reductive cleavage. This is the first report of a class III lanthipeptide from a myxobacterial origin.

Diversity and distribution of biosynthetic gene clusters in agricultural soil microbiomes

Bacterial secondary metabolites serve as an important source of molecules for drug discovery. They also play an important function in mediating the interactions of microbial producers with their living environment and surrounding organisms. However, little is known about the genetic novelty, distribution, and community-level impacts of soil bacterial biosynthetic potential on a large geographic scale. Here, we constructed the first catalog of 11,149 biosynthetic gene clusters (BGCs) from agricultural soils across China and unearthed hidden biosynthetic potential for new natural product discovery from the not-yet-cultivated soil bacteria. Notably, we revealed soil pH as the strongest environmental driver of BGC biogeography and predicted that soil acidification and global climate change could damage the biosynthetic potential of the soil microbiome. The co-occurrence network of bacterial genomes revealed two BGC-rich species, i.e., Nocardia niigatensis from Actinobacteriota and PSRF01 from Acidobacteriota, as the module hub and connector, respectively, indicating their keystone positions in the soil microbial communities. We also uncovered a dominant role of BGC-inferred biotic interactions over environmental drivers in structuring the soil microbiome. Overall, this study achieved novel insights into the BGC landscape in agricultural soils of China, substantially expanding our understanding of the diversity and novelty of bacterial secondary metabolism and the potential role of secondary metabolites in microbiota assembly.IMPORTANCEBacterial secondary metabolites not only serve as the foundation for numerous therapeutics (e.g., antibiotics and anticancer drugs), but they also play critical ecological roles in mediating microbial interactions (e.g., competition and communication). However, our knowledge of bacterial secondary metabolism is limited to only a small fraction of cultured strains, thus restricting our comprehensive understanding of their diversity, novelty, and potential ecological roles in soil ecosystems. Here, we used culture-independent metagenomics to explore biosynthetic potentials in agricultural soils of China. Our analyses revealed a high degree of genetic diversity and novelty within biosynthetic gene clusters in agricultural soil environments, offering valuable insights for biochemists seeking to synthesize novel bioactive products. Furthermore, we uncovered the pivotal role of BGC-rich species in microbial communities and the significant relationship between BGC richness and microbial phylogenetic turnover. This information emphasizes the importance of biosynthetic potential in the assembly of microbial communities.

Phylogenomics of the psychoactive mushroom genus Psilocybe and evolution of the psilocybin biosynthetic gene cluster

Psychoactive mushrooms in the genus Psilocybe have immense cultural value and have been used for centuries in Mesoamerica. Despite the recent surge of interest in these mushrooms due to the psychotherapeutic potential of their natural alkaloid psilocybin, their phylogeny and taxonomy remain substantially incomplete. Moreover, the recent elucidation of the psilocybin biosynthetic gene cluster is known for only five of ~165 species of Psilocybe, four of which belong to only one of two major clades. We set out to improve the phylogeny of Psilocybe using shotgun sequencing of fungarium specimens, from which we obtained 71 metagenomes including from 23 types, and conducting phylogenomic analysis of 2,983 single-copy gene families to generate a fully supported phylogeny. Molecular clock analysis suggests the stem lineage of Psilocybe arose ~67 mya and diversified ~56 mya. We also show that psilocybin biosynthesis first arose in Psilocybe, with 4 to 5 possible horizontal transfers to other mushrooms between 40 and 9 mya. Moreover, predicted orthologs of the psilocybin biosynthetic genes revealed two distinct gene orders within the biosynthetic gene cluster that corresponds to a deep split within the genus, possibly a signature of two independent acquisitions of the cluster within Psilocybe.

ABC-HuMi: the Atlas of Biosynthetic Gene Clusters in the Human Microbiome

The human microbiome has emerged as a rich source of diverse and bioactive natural products, harboring immense potential for therapeutic applications. To facilitate systematic exploration and analysis of its biosynthetic landscape, we present ABC-HuMi: the Atlas of Biosynthetic Gene Clusters (BGCs) in the Human Microbiome. ABC-HuMi integrates data from major human microbiome sequence databases and provides an expansive repository of BGCs compared to the limited coverage offered by existing resources. Employing state-of-the-art BGC prediction and analysis tools, our database ensures accurate annotation and enhanced prediction capabilities. ABC-HuMi empowers researchers with advanced browsing, filtering, and search functionality, enabling efficient exploration of the resource. At present, ABC-HuMi boasts a catalog of 19 218 representative BGCs derived from the human gut, oral, skin, respiratory and urogenital systems. By capturing the intricate biosynthetic potential across diverse human body sites, our database fosters profound insights into the molecular repertoire encoded within the human microbiome and offers a comprehensive resource for the discovery and characterization of novel bioactive compounds. The database is freely accessible at https://www.ccb.uni-saarland.de/abc_humi/.

Receptor tyrosine kinase gene expression profiling of orbital rhabdomyosarcoma unveils MET as a potential biomarker and therapeutic target

Receptor tyrosine kinases (RTKs) serve as molecular targets for the development of novel personalized therapies in many malignancies. In the present study, expression pattern of receptor tyrosine kinases and its clinical significance in orbital RMS has been explored. Eighteen patients with histopathologically confirmed orbital RMS formed part of this study. Comprehensive q-PCR gene expression profiles of 19 RTKs were generated in the cases and controls. The patients were followed up for 59.53 ± 20.93 years. Clustering and statistical analysis tools were applied to identify the significant combination of RTKs associated with orbital rhabdomyosarcoma patients. mRNA overexpression of RTKs which included MET, AXL, EGFR was seen in 60-80% of cases; EGFR3, IGFR2, FGFR1, RET, PDGFR1, VEGFR2, PDGFR2 in 30-60% of cases; and EGFR4, FGFR3,VEGFR3 and ROS,IGFR1, EGFR1, FGFR2, VEGFR1 in 10-30% of cases. Immunoexpression of MET was seen in 89% of cases. A significant association was seen between MET mRNA and its protein expression. In all the cases MET gene expression was associated with worst overall survival (P = 0.03).There was a significant correlation of MET mRNA expression with RET, ROS, AXL, FGFR1, FGFR3, PDGFR1, IGFR1, VEGFR2, and EGFR3 genes. Association between MET gene and collective expression of RTKs was further evaluated by semi-supervised gene cluster analysis and Principal component analysis, which showed well-separated tumor clusters. MET gene overexpression could be a useful biomarker for identifying high risk orbital rhabdomyosarcoma patients. Well-separated tumor clusters confirmed the association between MET gene and collective expression of RTK genes. Therefore, the therapeutic potential of multi-kinase inhibitors targeting MET and the 9 other significant RTKs needs to be explored.

The selection landscape and genetic legacy of ancient Eurasians

The Holocene (beginning around 12,000 years ago) encompassed some of the most significant changes in human evolution, with far-reaching consequences for the dietary, physical and mental health of present-day populations. Using a dataset of more than 1,600 imputed ancient genomes1, we modelled the selection landscape during the transition from hunting and gathering, to farming and pastoralism across West Eurasia. We identify key selection signals related to metabolism, including that selection at the FADS cluster began earlier than previously reported and that selection near the LCT locus predates the emergence of the lactase persistence allele by thousands of years. We also find strong selection in the HLA region, possibly due to increased exposure to pathogens during the Bronze Age. Using ancient individuals to infer local ancestry tracts in over 400,000 samples from the UK Biobank, we identify widespread differences in the distribution of Mesolithic, Neolithic and Bronze Age ancestries across Eurasia. By calculating ancestry-specific polygenic risk scores, we show that height differences between Northern and Southern Europe are associated with differential Steppe ancestry, rather than selection, and that risk alleles for mood-related phenotypes are enriched for Neolithic farmer ancestry, whereas risk alleles for diabetes and Alzheimer's disease are enriched for Western hunter-gatherer ancestry. Our results indicate that ancient selection and migration were large contributors to the distribution of phenotypic diversity in present-day Europeans.

Hox, homology, and parsimony: An organismal perspective

Hox genes are important regulators in animal development. They often show a mosaic of conserved (e.g., longitudinal axis patterning) and lineage-specific novel functions (e.g., development of skeletal, sensory, or locomotory systems). Despite extensive research over the past decades, it remains controversial at which node in the animal tree of life the Hox cluster evolved. Its presence already in the last common metazoan ancestor has been proposed, although the genomes of both putative earliest extant metazoan offshoots, the ctenophores and the poriferans, are devoid of Hox sequences. The lack of Hox genes in the supposedly "simple"-built poriferans and their low number in cnidarians and the basally branching bilaterians, the xenacoelomorphs, seems to support the classical notion that the number of Hox genes is correlated with the degree of animal complexity. However, the 4-fold increase of the Hox cluster in xiphosurans, a basally branching chelicerate clade, as well as the situation in some teleost fishes that show a multitude of Hox genes compared to, e.g., human, demonstrates, that there is no per se direct correlation between organismal complexity and Hox number. Traditional approaches have tried to base homology on the morphological level on shared expression profiles of individual genes, but recent data have shown that, in particular with respect to Hox and other regulatory genes, complex gene-gene interactions rather than expression signatures of individual genes alone are responsible for shaping morphological traits during ontogeny. Accordingly, for sound homology assessments and reconstructions of character evolution on organ system level, additional independent datasets (e.g., morphological, developmental) need to be included in any such analyses. If supported by solid data, proposed structural homology should be regarded as valid and not be rejected solely on the grounds of non-parsimonious distribution of the character over a given phylogenetic topology.

The gene cluster associated with strong biofilm-formation capacity by histamine-producing Lentilactobacillus parabuchneri encodes a sortase-mediated pilus and is located on a plasmid

Histamine is a biogenic amine synthesized through the enzymatic decarboxylation of the amino acid histidine. It can accumulate at high concentrations in foods through the metabolism of certain bacteria, sometimes leading to adverse reactions in consumers. In cheese, histamine can accumulate at toxic levels; Lentilactobacillus parabuchneri has been identified the major cause of this problem. Previous studies have shown some L. parabuchneri strains to form biofilms on different surfaces, posing a contamination risk during cheese production, particularly for cheeses that are processed post-ripening (e.g., grating or slicing). The food contamination they cause can result in economic losses and even foodborne illness if histamine accumulates in the final product. The aim of the present work was to identify the genes of L. parabuchneri involved in biofilm formation, and to determine their function. The genomes of six strains with different biofilm-production capacities (strong, moderate and weak) were sequenced and analysed. A cluster of four genes, similar to those involved in sortase-mediated pilus formation, was identified in the strong biofilm-producers, suggesting it to have a role in surface adhesion. Cloning and heterologous expression in Lactococcus cremoris NZ9000 confirmed its functionality and involvement in adhesion and, therefore, in biofilm formation. PacBio sequencing showed this cluster to be located on a 33.4 kb plasmid, which might increase its chances of horizontal transmission. These findings provide insight into the genetic factors associated with biofilm formation in histamine-producing L. parabuchneri, and into the risks associated with this bacterium in cheese production.

Endophyte genomes support greater metabolic gene cluster diversity compared with non-endophytes in Trichoderma

Trichoderma is a cosmopolitan genus with diverse lifestyles and nutritional modes, including mycotrophy, saprophytism, and endophytism. Previous research has reported greater metabolic gene repertoires in endophytic fungal species compared to closely-related non-endophytes. However, the extent of this ecological trend and its underlying mechanisms are unclear. Some endophytic fungi may also be mycotrophs and have one or more mycoparasitism mechanisms. Mycotrophic endophytes are prominent in certain genera like Trichoderma, therefore, the mechanisms that enable these fungi to colonize both living plants and fungi may be the result of expanded metabolic gene repertoires. Our objective was to determine what, if any, genomic features are overrepresented in endophytic fungi genomes in order to undercover the genomic underpinning of the fungal endophytic lifestyle. Here we compared metabolic gene cluster and mycoparasitism gene diversity across a dataset of thirty-eight Trichoderma genomes representing the full breadth of environmental Trichoderma's diverse lifestyles and nutritional modes. We generated four new Trichoderma endophyticum genomes to improve the sampling of endophytic isolates from this genus. As predicted, endophytic Trichoderma genomes contained, on average, more total biosynthetic and degradative gene clusters than non-endophytic isolates, suggesting that the ability to create/modify a diversity of metabolites potential is beneficial or necessary to the endophytic fungi. Still, once the phylogenetic signal was taken in consideration, no particular class of metabolic gene cluster was independently associated with the Trichoderma endophytic lifestyle. Several mycoparasitism genes, but no chitinase genes, were associated with endophytic Trichoderma genomes. Most genomic differences between Trichoderma lifestyles and nutritional modes are difficult to disentangle from phylogenetic divergences among species, suggesting that Trichoderma genomes maybe particularly well-equipped for lifestyle plasticity. We also consider the role of endophytism in diversifying secondary metabolism after identifying the horizontal transfer of the ergot alkaloid gene cluster to Trichoderma.

Identification and characterization of a novel bacteriocin gene cluster in Lysinibacillus boronitolerans

Although the antibiotics inhibit or kill pathogens, the abuse leads to the resistance formation and even "Super Bacteria." Therefore, it is urgent to explore the natural and safe alternatives such as bacteriocin. In this study, an uncharacterized bacteriocin gene cluster for Lysinibacillus boronitolerans was first predicted by genome sequencing and bioinformatics analysis, of which including two biosynthetic genes, a regulatory gene, a transport-related gene, and six other genes. Subsequently, the 10.24-kb gene cluster was expressed in Escherichia coli BL21, and the lysate effectively inhibited the growths of pathogenic bacteria containing Bacillus pumilus, Bacillus velezensis, Pseudomonas syringae pv. tomato DC3000, and Xanthomonas axonopodis pv. manihotis. The antibacterial substance was purified by 70% ammonium sulfate precipitation and further identified by liquid chromatography-tandem mass spectrometry. The results showed that the antibacterial substance consisted of 44 amino acids and had 24.1% sequence identity with the cyanobacterin Piricyclamide 7005 E4 PirE4, a bacteriocin analogue. The minimal set of genes required for the biosynthesis of the antibacterial substance was determined by site-directed mutagenesis, suggesting both a transcriptional repressor and a phosphohydroxythreonine transaminase were essential. Subsequently, the evolution and conservation of the two proteins were analyzed among 22 Lysinibacillus species. Among them, the residues responsible for functions were identified. Collectively, our results set a solid foundation for investigation of the biosynthesis and application of bacteriocin.

The Hypervariable Tpr Multigene Family of Theileria Parasites, Defined by a Conserved, Membrane-Associated, C-Terminal Domain, Includes Several Copies with Defined Orthology Between Species

Multigene families often play an important role in host-parasite interactions. One of the largest multigene families in Theileria parva, the causative agent of East Coast fever, is the T. parva repeat (Tpr) gene family. The function of the putative Tpr proteins remains unknown. The initial publication of the T. parva reference genome identified 39 Tpr family open reading frames (ORFs) sharing a conserved C-terminal domain. Twenty-eight of these are clustered in a central region of chromosome 3, termed the "Tpr locus", while others are dispersed throughout all four nuclear chromosomes. The Tpr locus contains three of the four assembly gaps remaining in the genome, suggesting the presence of additional, as yet uncharacterized, Tpr gene copies. Here, we describe the use of long-read sequencing to attempt to close the gaps in the reference assembly of T. parva (located among multigene families clusters), characterize the full complement of Tpr family ORFs in the T. parva reference genome, and evaluate their evolutionary relationship with Tpr homologs in other Theileria species. We identify three new Tpr family genes in the T. parva reference genome and show that sequence similarity among paralogs in the Tpr locus is significantly higher than between genes outside the Tpr locus. We also identify sequences homologous to the conserved C-terminal domain in five additional Theileria species. Using these sequences, we show that the evolution of this gene family involves conservation of a few orthologs across species, combined with gene gains/losses, and species-specific expansions.

Efficient Protein Expression and Biosynthetic Gene Cluster Regulation in Bacillus subtilis Driven by a T7-BOOST System

Bacillus subtilis is a generally recognized as safe microorganism that is widely used for protein expression and chemical production, but has a limited number of genetic regulatory components compared with the Gram-negative model microorganism Escherichia coli. In this study, a two-module plug-and-play T7-based optimized output strategy for transcription (T7-BOOST) systems with low leakage expression and a wide dynamic range was constructed based on the inducible promoters Phy-spank and PxylA. The first T7 RNA polymerase-driven module was seamlessly integrated into the genome based on the CRISPR/Cpf1 system, while the second expression control module was introduced into low, medium, and high copy plasmids for characterization. As a proof of concept, the T7-BOOST systems were successfully employed for whole-cell catalysis production of γ-aminobutyric acid (109.8 g/L with a 98.0% conversion rate), expression of human αS1 casein and human lactoferrin, and regulation of exogenous lycopene biosynthetic gene cluster and endogenous riboflavin biosynthetic gene cluster. Overall, the T7-BOOST system serves as a stringent, controllable, and effective tool for regulating gene expression in B. subtilis.

Genome wide identification, structural characterization and phylogenetic analysis of High-Affinity potassium (HAK) ion transporters in common bean (Phaseolus vulgaris L.)

BACKGROUND

High-Affinity Potassium ions represent one of the most important and large group of potassium transporters. Although HAK genes have been studied in a variety of plant species, yet, remain unexplored in common bean.

RESULTS

In the current study, 20 HAK genes were identified in common bean genome. Super-family "K_trans" domain was found in all PvHAK genes. Signals for localization of PvHAK proteins were detected in cell membrane. Fifty three HAKs genes, across diverse plant species, were divided into 5 groups based on sequential homology. Twelve pairs of orthologs genes were found in various plant species. PvHAKs genes were distributed unequally on 7 chromosomes with maximum number (7) mapped on chromosome 2 while only 1 PvHAK found on each chromosome 1, 4, and 6. Tandem gene duplication was witnessed in 2 paralog pairs while 1 pair exhibited segmental gene duplication. Five groups were made in PvHAK gene family based on Phylogeny. Maximum PvHAKs (10) were detected in Group-V while group-II composed of only 1 PvHAK gene. Variation was witnessed in number and size of motifs, and structure of PvHAKs associated with different groups. Light and hormone responsive elements contributed 57 and 24% share, respectively, to cis regulatory elements. qRT-PCR based results revealed significant increase in expression of all 4 PvHAK genes under low-potassium stress.

CONCLUSION

The current study provides valuable information for further functional characterization and uncovering the molecular mechanism associated with Potassium transportation in plants.

Yeast-based heterologous production of the Colletochlorin family of fungal secondary metabolites

Transcriptomic studies have revealed that fungal pathogens of plants activate the expression of numerous biosynthetic gene clusters (BGC) exclusively when in presence of a living host plant. The identification and structural elucidation of the corresponding secondary metabolites remain challenging. The aim was to develop a polycistronic system for heterologous expression of fungal BGCs in Saccharomyces cerevisiae. Here we adapted a polycistronic vector for efficient, seamless and cost-effective cloning of biosynthetic genes using in vivo assembly (also called transformation-assisted recombination) directly in Escherichia coli followed by heterologous expression in S. cerevisiae. Two vectors were generated with different auto-inducible yeast promoters and selection markers. The effectiveness of these vectors was validated with fluorescent proteins. As a proof-of-principle, we applied our approach to the Colletochlorin family of molecules. These polyketide secondary metabolites were known from the phytopathogenic fungus Colletotrichum higginsianum but had never been linked to their biosynthetic genes. Considering the requirement for a halogenase, and by applying comparative genomics, we identified a BGC putatively involved in the biosynthesis of Colletochlorins in C. higginsianum. Following the expression of those genes in S. cerevisiae, we could identify the presence of the precursor Orsellinic acid, Colletochlorins and their non-chlorinated counterparts, the Colletorins. In conclusion, the polycistronic vectors described herein were adapted for the host S. cerevisiae and allowed to link the Colletochlorin compound family to their corresponding biosynthetic genes. This system will now enable the production and purification of infection-specific secondary metabolites of fungal phytopathogens. More widely, this system could be applied to any fungal BGC of interest.

Genome-wide identification and expression analysis of GA20ox and GA3ox genes during pod development in peanut

Background

Gibberellins (GAs) play important roles in regulating peanut growth and development. GA20ox and GA3ox are key enzymes involved in GA biosynthesis. These enzymes encoded by a multigene family belong to the 2OG-Fe (II) oxygenase superfamily. To date, no genome-wide comparative analysis of peanut AhGA20ox and AhGA3ox-encoding genes has been performed, and the roles of these genes in peanut pod development are not clear.

Methods

A whole-genome analysis of AhGA20ox and AhGA3ox gene families in peanut was carried out using bioinformatic tools. The expression of these genes at different stage of pod development was analyzed using qRT-PCR.

Results

In this study, a total of 15 AhGA20ox and five AhGA3ox genes were identified in peanut genome, which were distributed on 14 chromosomes. Phylogenetic analysis divided the GA20oxs and GA3oxs into three groups, but AhGA20oxs and AhGA3oxs in two groups. The conserved pattern of gene structure, cis-elements, and protein motifs further confirmed their evolutionary relationship in peanut. AhGA20ox and AhGA3ox genes were differential expressed at different stages of pod development. The strong expression of AhGA20ox1/AhGA20ox4, AhGA20ox12/AhGA20ox15, AhGA3ox1 and AhGA3ox4/AhGA3ox5 in S1-stage indicated that these genes could have a key role in controlling peg elongation. Furthermore, AhGA20ox and AhGA3ox also showed diverse expression patterns in different peanut tissues including leaves, main stems, flowers and inflorescences. Noticeably, AhGA20ox9/AhGA20ox11 and AhGA3o4/AhGA3ox5 were highly expressed in the main stem, whereas the AhGA3ox1 and AhGA20ox10 were strongly expressed in the inflorescence. The expression levels of AhGA20ox2/AhGA20ox3, AhGA20ox5/AhGA20ox6, AhGA20ox7/AhGA20ox8, AhGA20ox13/AhGA20ox14 and AhGA3ox2/AhGA3ox3 were high in the flowers, suggesting their involvement in flower development. These results provide a basis for deciphering the roles of AhGA20ox and AhGA3ox in peanut growth and development, especially in pod development.

Genome-wide identification, comparative analysis and functional roles in flavonoid biosynthesis of cytochrome P450 superfamily in pear (Pyrus spp.)

BACKGROUND

The cytochrome P450 (CYP) superfamily is the largest enzyme metabolism family in plants identified to date, and it is involved in many biological processes, including secondary metabolite biosynthesis, hormone metabolism and stress resistance. However, the P450 gene superfamily has not been well studied in pear (Pyrus spp.).

RESULTS

Here, the comprehensive identification and a comparative analysis of P450 superfamily members were conducted in cultivated and wild pear genomes. In total, 338, 299 and 419 P450 genes were identified in Chinese white pear, European pear and the wild pear, respectively. Based on the phylogenetic analyses, pear P450 genes were divided into ten clans, comprising 48 families. The motif and gene structure analyses further supported this classification. The expansion of the pear P450 gene family was attributed to whole-genome and single-gene duplication events. Several P450 gene clusters were detected, which have resulted from tandem and proximal duplications. Purifying selection was the major force imposed on the long-term evolution of P450 genes. Gene dosage balance, subfunctionalization and neofunctionalization jointly drove the retention and functional diversification of P450 gene pairs. Based on the association analysis between transcriptome expression profiles and flavonoid content during fruit development, three candidate genes were identified as being closely associated with the flavonoid biosynthesis, and the expression of one gene was further verified using qRT-PCR and its function was validated through transient transformation in pear fruit.

CONCLUSIONS

The study results provide insights into the evolution and biological functions of P450 genes in pear.

Functional in vitro and in vivo analysis of biosynthetic genes by heterologous expression in E. coli

Biosynthetic gene clusters of natural products often harbor genes of unknown function, which are difficult to characterize. Here, we present a protocol for the functional analysis in vitro and in vivo of these biosynthetic genes by heterologous expression in E. coli. We describe steps for the expression of genes of interest in an established E. coli strain optimized to heterologously express natural products. We then detail the expression of a His-tagged gene to deduce the specific function of the protein. For complete details on the use and execution of this protocol, please refer to Böhringer et al.1.

The complete sequence of a human Y chromosome

The human Y chromosome has been notoriously difficult to sequence and assemble because of its complex repeat structure that includes long palindromes, tandem repeats and segmental duplications1-3. As a result, more than half of the Y chromosome is missing from the GRCh38 reference sequence and it remains the last human chromosome to be finished4,5. Here, the Telomere-to-Telomere (T2T) consortium presents the complete 62,460,029-base-pair sequence of a human Y chromosome from the HG002 genome (T2T-Y) that corrects multiple errors in GRCh38-Y and adds over 30 million base pairs of sequence to the reference, showing the complete ampliconic structures of gene families TSPY, DAZ and RBMY; 41 additional protein-coding genes, mostly from the TSPY family; and an alternating pattern of human satellite 1 and 3 blocks in the heterochromatic Yq12 region. We have combined T2T-Y with a previous assembly of the CHM13 genome4 and mapped available population variation, clinical variants and functional genomics data to produce a complete and comprehensive reference sequence for all 24 human chromosomes.

Global analyses of biosynthetic gene clusters in phytobiomes reveal strong phylogenetic conservation of terpenes and aryl polyenes

There are gaps in our understandings on how did the evolutionary relationships among members of the phytobiomes shape their ability to produce tremendously complex specialized metabolites under the influence of plant host. To determine these relationships, we investigated the phylogenetic conservation of biosynthetic gene clusters (BGCs) on a global collection of 4,519 high-quality and nonredundant (out of 12,181) bacterial isolates and metagenome-assembled genomes from 47 different plant hosts and soil, by adopting three independent phylogenomic approaches (D-test, Pagel's λ, and consenTRAIT). We report that the BGCs are phylogenetically conserved to varying strengths and depths in their different classes. We show that the ability to produce specialized metabolites qualifies as a complex trait, and the depth of conservation is equivalent to ecologically relevant complex microbial traits. Interestingly, terpene and aryl polyene BGCs had the strongest phylogenetic conservation in the phytobiomes, but not in the soil microbiomes. Furthermore, we showed that terpenes are largely uncharacterized in phytobiomes and pinpointed specific clades that harbor potentially novel terpenes. Taken together, this study sheds light on the evolution of specialized metabolites' biosynthesis potential in phytobiomes under the influence of plant hosts and presents strategies to rationally guide the discovery of potentially novel classes of metabolites. IMPORTANCE This study expands our understandings of the biosynthetic potential of phytobiomes by using such worldwide and extensive collection of microbiomes from plants and soil. Apart from providing such vital resource for the plant microbiome researchers, this study provides fundamental insights into the evolution of biosynthetic gene clusters (BGCs) in phytobiomes under the influence of plant host. Specifically, we report that the strength of phylogenetic conservation in microbiomes varies for different classes of BGCs and is influenced as a result of plant host association. Furthermore, our results indicate that biosynthetic potential of specialized metabolites is deeply conserved equivalent to other complex and ecologically relevant microbial traits. Finally, for the most conserved class of specialized metabolites (terpenes), we identified clades harboring potentially novel class of molecules. Future studies could focus on plant-microbe coevolution and interactions through specialized metabolites building upon these findings.

Multi-algorithm cooperation research of WRKY genes under nitrogen stress in Panax notoginseng

WRKY transcription factors play an important role in the immune system and the innate defense response of plants. WRKY transcription factors have great feedback on nitrogen stress. In this study, bioinformatics was used to detect the WRKYs of Panax notoginseng (PnWRKYs). The response of PnWRKYs under nitrogen stress was also well studied. PnWRKYs were distributed on 11 chromosomes. According to PnWRKY and Arabidopsis thaliana WRKY (AtWRKY) domains, these PnWRKY proteins were divided into three groups by phylogenetic analysis. MEME analysis showed that almost every member contained motif 1 and motif 2. PlantCARE online predicted the cis-acting elements of the promoter. PnWRKY gene family members obtained 22 pairs of repeat fragments by collinearity analysis. The expression levels of PnWRKYs in different parts (roots, flowers, and leafs) were analyzed by the gene expression pattern. They reflected tissue-specific expressions. The qRT-PCR experiments were used to detect 74 PnWRKYs under nitrogen stress. The results showed that the expression levels of 8 PnWRKYs were significantly induced. The PnWRKY gene family may be involved in biotic/abiotic stresses and hormone induction. This study will not only lay the foundation to explore the functions of PnWRKYs but also provide candidate genes for the future improvement of P. notoginseng.

Bacillus spp. as Bio-factories for Antifungal Secondary Metabolites: Innovation Beyond Whole Organism Formulations

Several fungi act as parasites for crops causing huge annual crop losses at both pre- and post-harvest stages. For years, chemical fungicides were the solution; however, their wide use has caused environmental contamination and human health problems. For this reason, the use of biofungicides has been in practice as a green solution against fungal phytopathogens. In the context of a more sustainable agriculture, microbial biofungicides have the largest share among the commercial biocontrol products that are available in the market. Precisely, the genus Bacillus has been largely studied for the management of plant pathogenic fungi because they offer a chemically diverse arsenal of antifungal secondary metabolites, which have spawned a heightened industrial engrossment of it as a biopesticide. In this sense, it is indispensable to know the wide arsenal that Bacillus genus has to apply these products for sustainable agriculture. Having this idea in our minds, in this review, secondary metabolites from Bacillus having antifungal activity are chemically and structurally described giving details of their action against several phytopathogens. Knowing the current status of Bacillus secreted antifungals is the base for the goal to apply these in agriculture and it is addressed in depth in the second part of this review.

Genome Organization of Four Brazilian Xanthomonas albilineans Strains Does Not Correlate with Aggressiveness

An integrative approach combining genomics, transcriptomics, and cell biology is presented to address leaf scald disease, a major problem for the sugarcane industry. To gain insight into the biology of the causal agent, the complete genome sequences of four Brazilian Xanthomonas albilineans strains with differing virulence capabilities are presented and compared to the GPEPC73 reference strain and FJ1. Based on the aggressiveness index, different strains were compared: Xa04 and Xa11 are highly aggressive, Xa26 is intermediate, and Xa21 is the least, while, based on genome structure, Xa04 shares most of its genomic features with Xa26, and Xa11 share most of its genomic features with Xa21. In addition to presenting more clustered regularly interspaced short palindromic repeats (CRISPR) clusters, four more novel prophage insertions are present than the previously sequenced GPEPC73 and FJ1 strains. Incorporating the aggressiveness index and in vitro cell biology into these genome features indicates that disease establishment is not a result of a single determinant factor, as in most other Xanthomonas species. The Brazilian strains lack the previously described plasmids but present more prophage regions. In pairs, the most virulent and the least virulent share unique prophages. In vitro transcriptomics shed light on the 54 most highly expressed genes among the 4 strains compared to ribosomal proteins (RPs), of these, 3 outer membrane proteins. Finally, comparative albicidin inhibition rings and in vitro growth curves of the four strains also do not correlate with pathogenicity. In conclusion, the results disclose that leaf scald disease is not associated with a single shared characteristic between the most or the least pathogenic strains. IMPORTANCE An integrative approach is presented which combines genomics, transcriptomics, and cell biology to address leaf scald disease. The results presented here disclose that the disease is not associated with a single shared characteristic between the most pathogenic strains or a unique genomic pattern. Sequence data from four Brazilian strains are presented that differ in pathogenicity index: Xa04 and Xa11 are highly virulent, Xa26 is intermediate, and Xa21 is the least pathogenic strain, while, based on genome structure, Xa04 shares with Xa26, and Xa11 shares with X21 most of the genome features. Other than presenting more CRISPR clusters and prophages than the previously sequenced strains, the integration of aggressiveness and cell biology points out that disease establishment is not a result of a single determinant factor as in other xanthomonads.

Heterologous Production and Structure Determination of a New Lanthipeptide Sinosporapeptin Using a Cryptic Gene Cluster in an Actinobacterium Sinosporangium siamense

Lipolanthine is a subclass of lanthipeptide that has the modification of lipid moiety at the N-terminus. A cryptic biosynthetic gene cluster comprising four genes (sinA, sinKC, sinD, and sinE) involved in the biosynthesis of lipolanthine was identified in the genome of an actinobacterium Sinosporangium siamense. Heterologous coexpression of a precursor peptide coding gene sinA and lanthipeptide synthetase coding gene sinKC in the host Escherichia coli strain BL21(DE3) resulted in the synthesis of a new lanthipeptide, sinosporapeptin. It contained unusual amino acids, including one labionin and two dehydrobutyrine residues, as determined using NMR and MS analyses. Another coexpression experiment with two additional genes of decarboxylase (sinD) and N-acetyl transferase (sinE) resulted in the production of a lipolanthine-like modified sinosporapeptin.

Genomic and Transcriptional Profiles of Kelch-like (klhl) Gene Family in Polyploid Carassius Complex

Genome duplication supplies raw genetic materials and has been thought to be essential for evolutionary innovation and ecological adaptation. Here, we select Kelch-like (klhl) genes to study the evolution of the duplicated genes in the polyploid Carassius complex, including amphidiploid C. auratus and amphitriploid C. gibelio. Phylogenetic, chromosomal location and read coverage analyses indicate that most of Carassius klhl genes exhibit a 2:1 relationship with zebrafish orthologs and confirm two rounds of polyploidy, an allotetraploidy followed by an autotriploidy, occurred during Carassius evolution. The lineage-specific expansion and biased retention/loss of klhl genes are also found in Carassius. Transcriptome analyses across eight adult tissues and seven embryogenesis stages reveal varied expression dominance and divergence between the two species. The expression of klhls in response to Carassius herpesvirus 2 infection shows different expression changes corresponding to distinct herpesvirus resistances in three C. gibelio gynogenetic clones. Finally, we find that most C. gibelio klhl genes possess three alleles except eight genes that have lost one or two alleles due to genome rearrangement. The allele expression bias is prosperous for Cgklhl genes and varies during embryogenesis owning to the sequential expression manner of the alleles. The current study provides global insights into the genomic and transcriptional evolution of duplicated genes in a given superfamily resulting from multiple rounds of polyploidization.

Complete integration of carbene-transfer chemistry into biosynthesis

Biosynthesis is an environmentally benign and renewable approach that can be used to produce a broad range of natural and, in some cases, new-to-nature products. However, biology lacks many of the reactions that are available to synthetic chemists, resulting in a narrower scope of accessible products when using biosynthesis rather than synthetic chemistry. A prime example of such chemistry is carbene-transfer reactions1. Although it was recently shown that carbene-transfer reactions can be performed in a cell and used for biosynthesis2,3, carbene donors and unnatural cofactors needed to be added exogenously and transported into cells to effect the desired reactions, precluding cost-effective scale-up of the biosynthesis process with these reactions. Here we report the access to a diazo ester carbene precursor by cellular metabolism and a microbial platform for introducing unnatural carbene-transfer reactions into biosynthesis. The α-diazoester azaserine was produced by expressing a biosynthetic gene cluster in Streptomyces albus. The intracellularly produced azaserine was used as a carbene donor to cyclopropanate another intracellularly produced molecule-styrene. The reaction was catalysed by engineered P450 mutants containing a native cofactor with excellent diastereoselectivity and a moderate yield. Our study establishes a scalable, microbial platform for conducting intracellular abiological carbene-transfer reactions to functionalize a range of natural and new-to-nature products and expands the scope of organic products that can be produced by cellular metabolism.

Phylogeny, transcriptional profile, and auxin-induced phosphorylation modification characteristics of conserved PIN proteins in Moso bamboo (Phyllostachys edulis)

Auxin polar transport is an important way for auxin to exercise its function, and auxin plays an irreplaceable role in the rapid growth of Moso bamboo. We identified and performed the structural analysis of PIN-FORMED auxin efflux carriers in Moso bamboo and obtained a total of 23 PhePIN genes from five gene subfamilies. We also performed chromosome localization and intra- and inter-species synthesis analysis. Phylogenetic analyses of 216 PIN genes showed that PIN genes are relatively conserved in the evolution of the Bambusoideae and have undergone intra-family segment replication in Moso bamboo. The PIN genes' transcriptional patterns showed that the PIN1 subfamily plays a major regulatory role. PIN genes and auxin biosynthesis maintain a high degree of consistency in spatial and temporal distribution. Phosphoproteomics analysis identified many phosphorylated protein kinases that respond to auxin regulation through autophosphorylation and phosphorylation of PIN proteins. The protein interaction network showed that there is a plant hormone interaction regulatory network with PIN protein as the core. We provide a comprehensive PIN protein analysis that complements the auxin regulatory pathway in Moso bamboo and paves the way for further auxin regulatory studies in bamboo.

Dissection of 3D chromosome organization in Streptomyces coelicolor A3(2) leads to biosynthetic gene cluster overexpression

The soil-dwelling filamentous bacteria, Streptomyces, is widely known for its ability to produce numerous bioactive natural products. Despite many efforts toward their overproduction and reconstitution, our limited understanding of the relationship between the host's chromosome three dimension (3D) structure and the yield of the natural products escaped notice. Here, we report the 3D chromosome organization and its dynamics of the model strain, Streptomyces coelicolor, during the different growth phases. The chromosome undergoes a dramatic global structural change from primary to secondary metabolism, while some biosynthetic gene clusters (BGCs) form special local structures when highly expressed. Strikingly, transcription levels of endogenous genes are found to be highly correlated to the local chromosomal interaction frequency as defined by the value of the frequently interacting regions (FIREs). Following the criterion, an exogenous single reporter gene and even complex BGC can achieve a higher expression after being integrated into the chosen loci, which may represent a unique strategy to activate or enhance the production of natural products based on the local chromosomal 3D organization.

Spruce giga-genomes: structurally similar yet distinctive with differentially expanding gene families and rapidly evolving genes

Spruces (Picea spp.) are coniferous trees widespread in boreal and mountainous forests of the northern hemisphere, with large economic significance and enormous contributions to global carbon sequestration. Spruces harbor very large genomes with high repetitiveness, hampering their comparative analysis. Here, we present and compare the genomes of four different North American spruces: the genome assemblies for Engelmann spruce (Picea engelmannii) and Sitka spruce (Picea sitchensis) together with improved and more contiguous genome assemblies for white spruce (Picea glauca) and for a naturally occurring introgress of these three species known as interior spruce (P. engelmannii × glauca × sitchensis). The genomes were structurally similar, and a large part of scaffolds could be anchored to a genetic map. The composition of the interior spruce genome indicated asymmetric contributions from the three ancestral genomes. Phylogenetic analysis of the nuclear and organelle genomes revealed a topology indicative of ancient reticulation. Different patterns of expansion of gene families among genomes were observed and related with presumed diversifying ecological adaptations. We identified rapidly evolving genes that harbored high rates of non-synonymous polymorphisms relative to synonymous ones, indicative of positive selection and its hitchhiking effects. These gene sets were mostly distinct between the genomes of ecologically contrasted species, and signatures of convergent balancing selection were detected. Stress and stimulus response was identified as the most frequent function assigned to expanding gene families and rapidly evolving genes. These two aspects of genomic evolution were complementary in their contribution to divergent evolution of presumed adaptive nature. These more contiguous spruce giga-genome sequences should strengthen our understanding of conifer genome structure and evolution, as their comparison offers clues into the genetic basis of adaptation and ecology of conifers at the genomic level. They will also provide tools to better monitor natural genetic diversity and improve the management of conifer forests. The genomes of four closely related North American spruces indicate that their high similarity at the morphological level is paralleled by the high conservation of their physical genome structure. Yet, the evidence of divergent evolution is apparent in their rapidly evolving genomes, supported by differential expansion of key gene families and large sets of genes under positive selection, largely in relation to stimulus and environmental stress response.

OCA-T1 and OCA-T2 are coactivators of POU2F3 in the tuft cell lineage

Tuft cells are a rare chemosensory lineage that coordinates immune and neural responses to foreign pathogens in mucosal tissues1. Recent studies have also revealed tuft-cell-like human tumours2,3, particularly as a variant of small-cell lung cancer. Both normal and neoplastic tuft cells share a genetic requirement for the transcription factor POU2F3 (refs. 2,4), although the transcriptional mechanisms that generate this cell type are poorly understood. Here we show that binding of POU2F3 to the uncharacterized proteins C11orf53 and COLCA2 (renamed here OCA-T1/POU2AF2 and OCA-T2/POU2AF3, respectively) is critical in the tuft cell lineage. OCA-T1 and OCA-T2 are paralogues of the B-cell-specific coactivator OCA-B; all three proteins are encoded in a gene cluster and contain a conserved peptide that binds to class II POU transcription factors and a DNA octamer motif in a bivalent manner. We demonstrate that binding between POU2F3 and OCA-T1 or OCA-T2 is essential in tuft-cell-like small-cell lung cancer. Moreover, we generated OCA-T1-deficient mice, which are viable but lack tuft cells in several mucosal tissues. These findings reveal that the POU2F3-OCA-T complex is the master regulator of tuft cell identity and a molecular vulnerability of tuft-cell-like small-cell lung cancer.

Evolutionarily stable gene clusters shed light on the common grounds of pathogenicity in the Acinetobacter calcoaceticus-baumannii complex

Nosocomial pathogens of the Acinetobacter calcoaceticus-baumannii (ACB) complex are a cautionary example for the world-wide spread of multi- and pan-drug resistant bacteria. Aiding the urgent demand for novel therapeutic targets, comparative genomics studies between pathogens and their apathogenic relatives shed light on the genetic basis of human-pathogen interaction. Yet, existing studies are limited in taxonomic scope, sensing of the phylogenetic signal, and resolution by largely analyzing genes independent of their organization in functional gene clusters. Here, we explored more than 3,000 Acinetobacter genomes in a phylogenomic framework integrating orthology-based phylogenetic profiling and microsynteny conservation analyses. We delineate gene clusters in the type strain A. baumannii ATCC 19606 whose evolutionary conservation indicates a functional integration of the subsumed genes. These evolutionarily stable gene clusters (ESGCs) reveal metabolic pathways, transcriptional regulators residing next to their targets but also tie together sub-clusters with distinct functions to form higher-order functional modules. We shortlisted 150 ESGCs that either co-emerged with the pathogenic ACB clade or are preferentially found therein. They provide a high-resolution picture of genetic and functional changes that coincide with the manifestation of the pathogenic phenotype in the ACB clade. Key innovations are the remodeling of the regulatory-effector cascade connecting LuxR/LuxI quorum sensing via an intermediate messenger to biofilm formation, the extension of micronutrient scavenging systems, and the increase of metabolic flexibility by exploiting carbon sources that are provided by the human host. We could show experimentally that only members of the ACB clade use kynurenine as a sole carbon and energy source, a substance produced by humans to fine-tune the antimicrobial innate immune response. In summary, this study provides a rich and unbiased set of novel testable hypotheses on how pathogenic Acinetobacter interact with and ultimately infect their human host. It is a comprehensive resource for future research into novel therapeutic strategies.

The spread of multi- and pan-drug resistant bacterial pathogens is a worldwide threat to human health. Understanding the genetics of host colonization and infection can substantially help in devising novel ways of treatment. Acinetobacter baumannii, a nosocomial pathogen ranked top by the World Health Organization in the list of bacteria for which novel therapeutic approaches are needed, is a prime example. Here, we have carved out the genetic make-up that distinguishes A. baumannii and its pathogenic next relatives from other and mostly apathogenic Acinetobacter species. We found a rich spectrum of pathways and regulatory modules that reveal how the pathogens have modified biofilm formation, iron scavenging, and their carbohydrate metabolism to adapt to their human host. Among these, the capability to metabolize kynurenine is particularly intriguing. Humans produce this substance to contain bacterial invaders and to fine-tune the innate immune response. But A. baumannii and closely related pathogens found a way to feed on kynurenine. This suggests that the pathogens might be able to dysregulate the human immune response. In summary, our study substantially deepens the understanding of how a highly critical pathogen interacts with its host, which substantially eases the identification of novel targets for innovative therapeutic strategies.

Exploring Newer Biosynthetic Gene Clusters in Marine Microbial Prospecting

Marine microbes genetically evolved to survive varying salinity, temperature, pH, and other stress factors by producing different bioactive metabolites. These microbial secondary metabolites (SMs) are novel, have high potential, and could be used as lead molecule. Genome sequencing of microbes revealed that they have the capability to produce numerous novel bioactive metabolites than observed under standard in vitro culture conditions. Microbial genome has specific regions responsible for SM assembly, termed biosynthetic gene clusters (BGCs), possessing all the necessary genes to encode different enzymes required to generate SM. In order to augment the microbial chemo diversity and to activate these gene clusters, various tools and techniques are developed. Metagenomics with functional gene expression studies aids in classifying novel peptides and enzymes and also in understanding the biosynthetic pathways. Genome shuffling is a high-throughput screening approach to improve the development of SMs by incorporating genomic recombination. Transcriptionally silent or lower level BGCs can be triggered by artificially knocking promoter of target BGC. Additionally, bioinformatic tools like antiSMASH, ClustScan, NAPDOS, and ClusterFinder are effective in identifying BGCs of existing class for annotation in genomes. This review summarizes the significance of BGCs and the different approaches for detecting and elucidating BGCs from marine microbes.

Horizontal transfer and evolution of the biosynthetic gene cluster for benzoxazinoids in plants

Benzoxazinoids are a class of protective and allelopathic plant secondary metabolites that have been identified in multiple grass species and are encoded by the Bx biosynthetic gene cluster (BGC) in maize. Data mining of 41 high-quality grass genomes identified complete Bx clusters (containing genes Bx1-Bx5 and Bx8) in three genera (Zea, Echinochloa, and Dichanthelium) of Panicoideae and partial clusters in Triticeae. The Bx cluster probably originated from gene duplication and chromosomal translocation of native homologs of Bx genes. An ancient Bx cluster that included additional Bx genes (e.g., Bx6) is presumed to have been present in ancestral Panicoideae. The ancient Bx cluster was putatively gained by the Triticeae ancestor via horizontal transfer (HT) from the ancestral Panicoideae and later separated into multiple segments on different chromosomes. Bx6 appears to have been under less constrained selection compared with the Bx cluster during the evolution of Panicoideae, as evidenced by the fact that it was translocated away from the Bx cluster in Zea mays, moved to other chromosomes in Echinochloa, and even lost in Dichanthelium. Further investigations indicate that purifying selection and polyploidization have shaped the evolutionary trajectory of Bx clusters in the grass family. This study provides the first candidate case of HT of a BGC between plants and sheds new light on the evolution of BGCs.

Pathogen-induced biosynthetic pathways encode defense-related molecules in bread wheat

Wheat is a widely grown food crop that suffers major yield losses due to attack by pests and pathogens. A better understanding of biotic stress responses in wheat is thus of major importance. The recently assembled bread wheat genome coupled with extensive transcriptomic resources provides unprecedented new opportunities to investigate responses to pathogen challenge. Here, we analyze gene coexpression networks to identify modules showing consistent induction in response to pathogen exposure. Within the top pathogen-induced modules, we identify multiple clusters of physically adjacent genes that correspond to six pathogen-induced biosynthetic pathways that share a common regulatory network. Functional analysis reveals that these pathways, all of which are encoded by biosynthetic gene clusters, produce various different classes of compounds—namely, flavonoids, diterpenes, and triterpenes, including the defense-related compound ellarinacin. Through comparative genomics, we also identify associations with the known rice phytoalexins momilactones, as well as with a defense-related gene cluster in the grass model plant Brachypodium distachyon. Our results significantly advance the understanding of chemical defenses in wheat and open up avenues for enhancing disease resistance in this agriculturally important crop. They also exemplify the power of transcriptional networks to discover the biosynthesis of chemical defenses in plants with large, complex genomes.

New Report: Genome Mining Untaps the Antibiotics Biosynthetic Gene Cluster of Pseudoalteromonas xiamenensis STKMTI.2 from a Mangrove Soil Sediment

The marine bacterium Pseudoalteromonas xiamenensis STKMTI.2 was isolated from a mangrove soil sediment on Setokok Island, Batam, Indonesia. The genome of this bacterium consisted of 4,563,326 bp (GC content: 43.2%) with 1 chromosome, 2 circular plasmids, 2 linear plasmids, 4,824 protein-coding sequences, 25 rRNAs, 104 tRNAs, 4 ncRNAs, and 1 clustered, regularly interspaced, short palindromic repeated (CRISPR). This strain possessed cluster genes which are responsible for the production of brominated marine pyrroles/phenols (bmp), namely, bmp8 and bmp9. Other gene clusters responsible for the synthesis of secondary metabolites were identified using antiSMASH and BAGEL4, which yielded five results, namely, non-ribosomal peptides, polyketide-like butyrolactone, Lant class I, and RiPP-like, detected in chromosome 1, while prodigiosin was detected in the unnamed plasmid 5. This suggests that these whole genome data will be of remarkable importance for the improved understanding of the biosynthesis of industrially important bioactive and antibacterial compounds produced by P. xiamenensis STKMTI.2.

Genome-Wide Analysis of LBD Transcription Factor Genes in Dendrobiumcatenatum

The LATERAL ORGAN BOUNDARIES DOMAIN (LBD) gene family comprises plant-specific transcription factors that control cell proliferation and differentiation during growth and development in many plant species. However, to date, no studies of the LBD gene family in Dendrobium catenatum have been reported. In this study, a genome-wide analysis of LBD genes was performed in D. catenatum and 24 LBD genes were identified. The genes were classified into two classes (I and II) based on phylogenetic relationships and motif structure. Subcellular localization analysis for DcaLBD6 and DcaLBD18 from class I and DcaLBD37 and DcaLBD41 from class II revealed that the proteins were localized in the nucleus. Transient expression analysis of DcaLBD6, DcaLBD18, DcaLBD37, and DcaLBD41 indicated that class I and class II members have opposite roles in regulating VASCULAR-RELATED NAC-DOMAIN 7 (VND7) expression. DcaLBD genes showed diverse expression patterns in response to different phytohormone treatments. Heat maps revealed diverse patterns of DcaLBD gene expression in different organs. These results lay the foundation for further detailed studies of the LBD gene family in D. catenatum.

Genome-Wide Identification and Characterization of SET Domain Family Genes in Brassica napus L

SET domain group encoding proteins function as histone lysine methyltransferases. These proteins are involved in various biological processes, including plant development and adaption to the environment by modifying the chromatin structures. So far, the SET domain genes (SDGs) have not been systematically investigated in Brassica napus (B. napus). In the current study, through genome-wide analysis, a total of 122 SDGs were identified in the B. napus genome. These BnSDGs were subdivided into seven (I-VII) classes based on phylogeny analysis, domain configurations, and motif distribution. Segmental duplication was involved in the evolution of this family, and the duplicated genes were under strong purifying selection. The promoter sequence of BnSDGs consisted of various growth, hormones, and stress-related cis-acting elements along with transcription factor binding sites (TFBSs) for 20 TF families in 59 of the 122 BnSDGs. The gene ontology (GO) analysis revealed that BnSDGs were closely associated with histone and non-histone methylation and metal binding capacity localized mostly in the nucleus. The in silico expression analysis at four developmental stages in leaf, stem root, floral organ, silique, and seed tissues showed a broad range of tissue and stage-specific expression pattern. The expression analysis under four abiotic stresses (dehydration, cold, ABA, and salinity) also provided evidence for the importance of BnSDGs in stress environments. Based on expression analysis, we performed reverse transcription-quantitative PCR for 15 target BnSDGs in eight tissues (young leaf, mature leaf, root, stem, carpel, stamen, sepal, and petals). Our results were in accordance with the in silico expression data, suggesting the importance of these genes in plant development. In conclusion, this study lays a foundation for future functional studies on SDGs in B. napus.

Protein with negative surface charge distribution, Bnr1, shows characteristics of a DNA-mimic protein and may be involved in the adaptation of Burkholderia cenocepacia

Adaptation of opportunistic pathogens to their host environment requires reprogramming of a vast array of genes to facilitate survival in the host. Burkholderia cenocepacia, a Gram-negative bacterium with a large genome of ∼8 Mb that colonizes environmental niches, is exquisitely adaptable to the hypoxic environment of the cystic fibrosis lung and survives in macrophages. We previously identified an immunoreactive acidic protein encoded on replicon 3, BCAS0292. Deletion of the BCAS0292 gene significantly altered the abundance of 979 proteins by 1.5-fold or more; 19 proteins became undetectable while 545 proteins showed ≥1.5-fold reduced abundance, suggesting the BCAS0292 protein is a global regulator. Moreover, the ∆BCAS0292 mutant showed a range of pleiotropic effects: virulence and host-cell attachment were reduced, antibiotic susceptibility was altered, and biofilm formation enhanced. Its growth and survival were impaired in 6% oxygen. In silico prediction of its three-dimensional structure revealed BCAS0292 presents a dimeric β-structure with a negative surface charge. The ΔBCAS0292 mutant displayed altered DNA supercoiling, implicated in global regulation of gene expression. Three proteins were identified in pull-downs with FLAG-tagged BCAS0292, including the Histone H1-like protein, HctB, which is recognized as a global transcriptional regulator. We propose that BCAS0292 protein, which we have named Burkholderia negatively surface-charged regulatory protein 1 (Bnr1), acts as a DNA-mimic and binds to DNA-binding proteins, altering DNA topology and regulating the expression of multiple genes, thereby enabling the adaptation of B. cenocepacia to highly diverse environments.

Cholesterol transport in the late endocytic pathway: Roles of ORP family proteins

Oxysterol-binding protein (OSBP) homologues, designated ORP or OSBPL proteins, constitute one of the largest families of intracellular lipid-binding/transfer proteins (LTP). This review summarizes the mounting evidence that several members of this family participate in the machinery facilitating cholesterol trafficking in the late endocytic pathway. There are indications that OSBP, besides acting as a cholesterol/phosphatidylinositol 4-phosphate (PI4P) exchanger at the endoplasmic reticulum (ER)-trans-Golgi network (TGN) membrane contact sites (MCS), also exchanges these lipids at ER-lysosome (Lys) contacts, increasing Lys cholesterol content. The long isoform of ORP1 (ORP1L), which also targets ER-late endosome (LE)/Lys MCS, has the capacity to mediate cholesterol transport either from ER to LE or in the opposite direction. Moreover, it regulates the motility, positioning and fusion of LE as well as autophagic flux. ORP2, the closest relative of ORP1, is mainly cytosolic, but also targets PI(4,5)P₂-rich endosomal compartments. Our latest data suggest that ORP2 transfers cholesterol from LE to recycling endosomes (RE) in exchange for PI(4,5)P₂, thus stimulating the recruitment of focal adhesion kinase (FAK) on the RE and cell adhesion. FAK activates phosphoinositide kinase on the RE to enhance PI(4,5)P₂ synthesis. ORP2 in turn transfers PI(4,5)P₂ from RE to LE, thus regulating LE tubule formation and transport activity.

Isolation and characterization of Streptomyces bacteriophages and Streptomyces strains encoding biosynthetic arsenals

The threat to public health posed by drug-resistant bacteria is rapidly increasing, as some of healthcare’s most potent antibiotics are becoming obsolete. Approximately two-thirds of the world’s antibiotics are derived from natural products produced by Streptomyces encoded biosynthetic gene clusters. Thus, to identify novel gene clusters, we sequenced the genomes of four bioactive Streptomyces strains isolated from the soil in San Diego County and used Bacterial Cytological Profiling adapted for agar plate culturing in order to examine the mechanisms of bacterial inhibition exhibited by these strains. In the four strains, we identified 104 biosynthetic gene clusters. Some of these clusters were predicted to produce previously studied antibiotics; however, the known mechanisms of these molecules could not fully account for the antibacterial activity exhibited by the strains, suggesting that novel clusters might encode antibiotics. When assessed for their ability to inhibit the growth of clinically isolated pathogens, three Streptomyces strains demonstrated activity against methicillin-resistant Staphylococcus aureus. Additionally, due to the utility of bacteriophages for genetically manipulating bacterial strains via transduction, we also isolated four new phages (BartholomewSD, IceWarrior, Shawty, and TrvxScott) against S. platensis. A genomic analysis of our phages revealed nearly 200 uncharacterized proteins, including a new site-specific serine integrase that could prove to be a useful genetic tool. Sequence analysis of the Streptomyces strains identified CRISPR-Cas systems and specific spacer sequences that allowed us to predict phage host ranges. Ultimately, this study identified Streptomyces strains with the potential to produce novel chemical matter as well as integrase-encoding phages that could potentially be used to manipulate these strains.

Genome-wide analysis of the NAAT, DMAS, TOM, and ENA gene families in maize suggests their roles in mediating iron homeostasis

BACKGROUND

Nicotianamine (NA), 2'-deoxymugineic acid (DMA), and mugineic acid (MA) are chelators required for iron uptake and transport in plants. Nicotianamine aminotransferase (NAAT), 2'-deoxymugineic acid synthase (DMAS), transporter of MAs (TOM), and efflux transporter of NA (ENA) are involved in iron uptake and transport in rice (Oryza sativa), wheat (Triticum aestivum), and barley (Hordeum vulgare); however, these families have not been fully identified and comprehensively analyzed in maize (Zea mays L.).

RESULTS

Here, we identified 5 ZmNAAT, 9 ZmDMAS, 11 ZmTOM, and 2 ZmENA genes by genome mining. RNA-sequencing and quantitative real-time PCR analysis revealed that these genes are expressed in various tissues and respond differently to high and low iron conditions. In particular, iron deficiency stimulated the expression of ZmDMAS1, ZmTOM1, ZmTOM3, and ZmENA1. Furthermore, we determined protein subcellular localization by transient expression of green fluorescent protein fusions in maize mesophyll protoplasts. ZmNAAT1, ZmNAAT-L4, ZmDMAS1, and ZmDMAS-L1 localized in the cytoplasm, whereas ZmTOMs and ZmENAs targeted to plasma and tonoplast membranes, endomembranes, and vesicles.

CONCLUSIONS

Our results suggest that the different gene expression profiles and subcellular localizations of ZmNAAT, ZmDMAS, ZmTOM, and ZmENA family members may enable specific regulation of phytosiderophore metabolism in different tissues and under different external conditions, shedding light on iron homeostasis in maize and providing candidate genes for breeding iron-rich maize varieties.

Transcriptional regulation by CRISPR/dCas9 in common wheat

The application of CRISPR/Cas9 system for gene editing, as a technical coup for biotechnology, is worldwide and encompasses multiple of species. The inactivation of catalytical site in Cas9 (dCas9) has been reprogrammed as an effective approach to regulate the transcriptional level of target genes, especially for the functionally essential genes and redundant genes. Here, we exploited the CRISPR/dCas9 system to manipulate the transcriptional level of target genes in common wheat. To improve target gene's expression, we generated transcriptional activator by fusing 6×TAL-VP128 activation domain to the C-terminus of dCas9 in frame. For target gene's repressing expression transcriptionally, 3×SRDX repression domain was conjugated to the C-terminus of dCas9 in frame. Our results showed that dCas9 fused activation or repression domain could increase or decrease the transcriptional level of target gene effectively in stable transgenic lines of wheat. The study on the tRNA-processing system in CRISPR/dCas9 based transcriptional regulation system demonstrated that this robust multiplex targeted tool can be incorporated to the CRISPR/dCas9 system to facilitate the target regulation of several genes' transcriptional level. Our data broaden the application of CRISPR/dCas9 based transcriptional regulation and provide great opportunities to investigate the function of essential genes in common wheat.

Genome-Wide Identification, Characterization and Expression Analysis of the CIPK Gene Family in Potato (Solanum tuberosum L.) and the Role of StCIPK10 in Response to Drought and Osmotic Stress

The potato (Solanum tuberosum L.), one of the most important food crops worldwide, is sensitive to environmental stresses. Sensor–responder complexes comprising calcineurin B-like (CBL) proteins and CBL-interacting protein kinases (CIPKs) not only modulate plant growth and development but also mediate numerous stress responses. Here, using a Hidden Markov Model and BLAST searches, 27 CIPK genes were identified in potato and divided into five groups by phylogenetic analysis and into two clades (intron-poor and intron-rich) by gene structure analysis. Quantitative reverse-transcription PCR (qRT-PCR) assays revealed that StCIPK genes play important roles in plant growth, development and abiotic stress tolerance. Up-regulated expression of StCIPK10 was significantly induced by drought, PEG6000 and ABA. StCIPK10 enhances both the ability of potato to scavenge reactive oxygen species and the content of corresponding osmoregulation substances, thereby strengthening tolerance to drought and osmotic stress. StCIPK10 is located at the intersection between the abscisic acid and abiotic stress signaling pathways, which control both root growth and stomatal closure in potato. In addition, StCIPK10 interacts with StCBL1, StCBL4, StCBL6, StCBL7, StCBL8, StCBL11 and StCBL12, and is specifically recruited to the plasma membrane by StCBL11.

An Oryza-specific hydroxycinnamoyl tyramine gene cluster contributes to enhanced disease resistance

Genomic clustering of non-homologous genes for the biosynthesis of plant defensive compounds is an emerging theme, but insights into their formation and physiological function remain limited. Here we report the identification of a newly discovered hydroxycinnamoyl tyramine (HT) gene cluster in rice. This cluster contains a pyridoxamine 5'-phosphate oxidase (OsPDX3) producing the cofactor pyridoxal 5'-phosphate (PLP), a PLP-dependent tyrosine decarboxylase (OsTyDC1), and two duplicated hydroxycinnamoyl transferases (OsTHT1 and OsTHT2). These members were combined to represent an enzymological innovation gene cluster. Natural variation analysis showed that the abundance of the toxic tyramine intermediate of the gene cluster among different rice accessions is mainly determined by the coordinated transcription of OsTyDC1 and OsTHT1. Further pathogen incubation assays demonstrated that the end products of the HT gene cluster displayed enhanced resistance to the bacterial pathogen Xanthomonas oryzae pv. Oryzae (Xoo) and fungal pathogen Magnaporthe oryzae (M. oryzae), and the enhanced resistance is associated with the boost of phytoalexins and the activation of defense response. The unique presence of the HT gene cluster in Oryza AA genome, together with the enrichment of transposon elements within this gene cluster region, provides an evolutionary background to accelerate cluster member combinations. Our study not only discovered a gene cluster involved in the phenylpropanoid metabolism but also addressed the key aspects of gene cluster formation. In addition, our results provide a new metabolic pool for plant defense against pathogens.

Comprehensive Analysis of the SUV Gene Family in Allopolyploid Brassica napus and Its Diploid Ancestors

SUV (the Suppressor of variegation [Su(var)] homologs and related) gene family is a subgroup of the SET gene family. According to the SRA domain and WIYLD domain distributions, it can be divided into two categories, namely SUVH (the Suppressor of variegation [Su(var)] homologs) and SUVR (the Suppressor of variegation [Su(var)] related). In this study, 139 SUV genes were identified in allopolyploid Brassica napus and its diploid ancestors, and their evolutionary relationships, protein properties, gene structures, motif distributions, transposable elements, cis-acting elements and gene expression patterns were analyzed. Our results showed that the SUV gene family of B. napus was amplified during allopolyploidization, in which the segmental duplication and TRD played critical roles. After the separation of Brassica and Arabidopsis lineages, orthologous gene analysis showed that many SUV genes were lost during the evolutionary process in B. rapa, B. oleracea and B. napus. The analysis of the gene and protein structures and expression patterns of 30 orthologous gene pairs which may have evolutionary relationships showed that most of them were conserved in gene structures and protein motifs, but only four gene pairs had the same expression patterns.

Genome-Wide Identification and Characterization of Soybean GmLOR Gene Family and Expression Analysis in Response to Abiotic Stresses

The LOR (LURP-one related) family genes encode proteins containing a conserved LOR domain. Several members of the LOR family genes are required for defense against Hyaloperonospora parasitica (Hpa) in Arabidopsis. However, there are few reports of LOR genes in response to abiotic stresses in plants. In this study, a genome-wide survey and expression levels in response to abiotic stresses of 36 LOR genes from Glycine max were conducted. The results indicated that the GmLOR gene family was divided into eight subgroups, distributed on 14 chromosomes. A majority of members contained three extremely conservative motifs. There were four pairs of tandem duplicated GmLORs and nineteen pairs of segmental duplicated genes identified, which led to the expansion of the number of GmLOR genes. The expansion patterns of the GmLOR family were mainly segmental duplication. A heatmap of soybean LOR family genes showed that 36 GmLOR genes exhibited various expression patterns in different tissues. The cis-acting elements in promoter regions of GmLORs include abiotic stress-responsive elements, such as dehydration-responsive elements and drought-inducible elements. Real-time quantitative PCR was used to detect the expression level of GmLOR genes, and most of them were expressed in the leaf or root except that GmLOR6 was induced by osmotic and salt stresses. Moreover, GmLOR4/10/14/19 were significantly upregulated after PEG and salt treatments, indicating important roles in the improvement of plant tolerance to abiotic stress. Overall, our study provides a foundation for future investigations of GmLOR gene functions in soybean.

Epromoters function as a hub to recruit key transcription factors required for the inflammatory response

Gene expression is controlled by the involvement of gene-proximal (promoters) and distal (enhancers) regulatory elements. Our previous results demonstrated that a subset of gene promoters, termed Epromoters, work as bona fide enhancers and regulate distal gene expression. Here, we hypothesized that Epromoters play a key role in the coordination of rapid gene induction during the inflammatory response. Using a high-throughput reporter assay we explored the function of Epromoters in response to type I interferon. We find that clusters of IFNa-induced genes are frequently associated with Epromoters and that these regulatory elements preferentially recruit the STAT1/2 and IRF transcription factors and distally regulate the activation of interferon-response genes. Consistently, we identified and validated the involvement of Epromoter-containing clusters in the regulation of LPS-stimulated macrophages. Our findings suggest that Epromoters function as a local hub recruiting the key TFs required for coordinated regulation of gene clusters during the inflammatory response.

In-silico genome wide analysis of Mitogen activated protein kinase kinase kinase gene family in C. sinensis

Mitogen activated protein kinase kinase kinase (MAPKKK) form the upstream component of MAPK cascade. It is well characterized in several plants such as Arabidopsis and rice however the knowledge about MAPKKKs in tea plant is largely unknown. In the present study, MAPKKK genes of tea were obtained through a genome wide search using Arabidopsis thaliana as the reference genome. Among 59 candidate MAPKKK genes in tea, 17 genes were MEKK-like, 31 genes were Raf-like and 11 genes were ZIK- like. Additionally, phylogenetic relationships were established along with structural analysis, which includes gene structure, its location as well as conserved motifs, cis-acting regulatory elements and functional domain signatures that were systematically examined. Also, on the basis of one orthologous gene found between tea and Arabidopsis, functional interaction was carried out in C. sinensis based on an Arabidopsis association model. The expressional profiles indicated major involvement of MAPKKK genes from tea in response to various abiotic stress factors. Taken together, this study provides the targets for additional inclusive identification, functional study, and provides comprehensive knowledge for a better understanding of the MAPKKK cascade regulatory network in C. sinensis.

Cell-type specialization is encoded by specific chromatin topologies

The three-dimensional (3D) structure of chromatin is intrinsically associated with gene regulation and cell function1-3. Methods based on chromatin conformation capture have mapped chromatin structures in neuronal systems such as in vitro differentiated neurons, neurons isolated through fluorescence-activated cell sorting from cortical tissues pooled from different animals and from dissociated whole hippocampi4-6. However, changes in chromatin organization captured by imaging, such as the relocation of Bdnf away from the nuclear periphery after activation7, are invisible with such approaches8. Here we developed immunoGAM, an extension of genome architecture mapping (GAM)2,9, to map 3D chromatin topology genome-wide in specific brain cell types, without tissue disruption, from single animals. GAM is a ligation-free technology that maps genome topology by sequencing the DNA content from thin (about 220 nm) nuclear cryosections. Chromatin interactions are identified from the increased probability of co-segregation of contacting loci across a collection of nuclear slices. ImmunoGAM expands the scope of GAM to enable the selection of specific cell types using low cell numbers (approximately 1,000 cells) within a complex tissue and avoids tissue dissociation2,10. We report cell-type specialized 3D chromatin structures at multiple genomic scales that relate to patterns of gene expression. We discover extensive 'melting' of long genes when they are highly expressed and/or have high chromatin accessibility. The contacts most specific of neuron subtypes contain genes associated with specialized processes, such as addiction and synaptic plasticity, which harbour putative binding sites for neuronal transcription factors within accessible chromatin regions. Moreover, sensory receptor genes are preferentially found in heterochromatic compartments in brain cells, which establish strong contacts across tens of megabases. Our results demonstrate that highly specific chromatin conformations in brain cells are tightly related to gene regulation mechanisms and specialized functions.