Transcriptome analysis of the phosphate starvation response sheds light on strigolactone biosynthesis in rice.

Phosphorus (P) is a major element required for plant growth and development. To cope with P shortage, plants activate local and long-distance signaling pathways, such as an increase in the production and exudation of strigolactones (SLs). The role of the latter in mitigating P deficiency is, however, still largely unknown. To shed light on this, we studied the transcriptional response to P starvation and replenishment in wild-type rice and a SL mutant, dwarf10 (d10), and upon exogenous application of the synthetic SL GR24. P starvation resulted in major transcriptional alterations, such as the upregulation of P TRANSPORTER, SYG1/PHO81/XPR1 (SPX) and VACUOLAR PHOSPHATE EFFLUX TRANSPORTER. Gene Ontology (GO) analysis of the genes induced by P starvation showed enrichment in phospholipid catabolic process and phosphatase activity. In d10, P deficiency induced upregulation of genes enriched for sesquiterpenoid production, secondary shoot formation and metabolic processes, including lactone biosynthesis. Furthermore, several genes induced by GR24 treatment shared the same GO terms with P starvation-induced genes, such as oxidation reduction, heme binding and oxidoreductase activity, hinting at the role that SLs play in the transcriptional reprogramming upon P starvation. Gene co-expression network analysis uncovered a METHYL TRANSFERASE that displayed co-regulation with known rice SL biosynthetic genes. Functional characterization showed that this gene encodes an enzyme catalyzing the conversion of carlactonoic acid to methyl carlactonoate. Our work provides a valuable resource to further studies on the response of crops to P deficiency and reveals a tool for the discovery of SL biosynthetic genes.

Genome-wide expression analysis of vegetative organs during developmental and herbicide-induced whole plant senescence in Arabidopsis thaliana.

BACKGROUND: Whole plant senescence represents the final stage in the life cycle of annual plants, characterized by the decomposition of aging organs and transfer of nutrients to seeds, thereby ensuring the survival of next generation. However, the transcriptomic profile of vegetative organs during this death process remains to be fully elucidated, especially regarding the distinctions between natural programmed death and artificial sudden death induced by herbicide. RESULTS: Differential genes expression analysis using RNA-seq in leaves and roots of Arabidopsis thaliana revealed that natural senescence commenced in leaves at 45-52 days after planting, followed by roots initiated at 52-60 days. Additionally, both organs exhibited similarities with artificially induced senescence by glyphosate. Transcription factors Rap2.6L and WKRY75 appeared to serve as central mediators of regulatory changes during natural senescence, as indicated by co-expression networks. Furthermore, the upregulation of RRTF1, exclusively observed during natural death, suggested its role as a regulator of jasmonic acid and reactive oxygen species (ROS) responses, potentially triggering nitrogen recycling in leaves, such as the glutamate dehydrogenase (GDH) shunt. Root senescence was characterized by the activation of AMT2;1 and GLN1;3, facilitating ammonium availability for root-to-shoot translocation, likely under the regulation of PDF2.1. CONCLUSIONS: Our study offers valuable insights into the transcriptomic interplay between phytohormones and ROS during whole plant senescence. We observed distinct regulatory networks governing nitrogen utilization in leaf and root senescence processes. Furthermore, the efficient allocation of energy from vegetative organs to seeds emerges as a critical determinant of population sustainability of annual Arabidopsis.

Computational network analysis of host genetic risk variants of severe COVID-19.

BACKGROUND: Genome-wide association studies have identified numerous human host genetic risk variants that play a substantial role in the host immune response to SARS-CoV-2. Although these genetic risk variants significantly increase the severity of COVID-19, their influence on body systems is poorly understood. Therefore, we aim to interpret the biological mechanisms and pathways associated with the genetic risk factors and immune responses in severe COVID-19. We perform a deep analysis of previously identified risk variants and infer the hidden interactions between their molecular networks through disease mapping and the similarity of the molecular functions between constructed networks. RESULTS: We designed a four-stage computational workflow for systematic genetic analysis of the risk variants. We integrated the molecular profiles of the risk factors with associated diseases, then constructed protein-protein interaction networks. We identified 24 protein-protein interaction networks with 939 interactions derived from 109 filtered risk variants in 60 risk genes and 56 proteins. The majority of molecular functions, interactions and pathways are involved in immune responses; several interactions and pathways are related to the metabolic and cardiovascular systems, which could lead to multi-organ complications and dysfunction. CONCLUSIONS: This study highlights the importance of analyzing molecular interactions and pathways to understand the heterogeneous susceptibility of the host immune response to SARS-CoV-2. We propose new insights into pathogenicity analysis of infections by including genetic risk information as essential factors to predict future complications during and after infection. This approach may assist more precise clinical decisions and accurate treatment plans to reduce COVID-19 complications.

Population structure of indigenous inhabitants of Arabia.

Modern day Saudi Arabia occupies the majority of historical Arabia, which may have contributed to ancient waves of migration out of Africa. This ancient history has left a lasting imprint in the genetics of the region, including the diverse set of tribes that call Saudi Arabia their home. How these tribes relate to each other and to the world's major populations remains an unanswered question. In an attempt to improve our understanding of the population structure of Saudi Arabia, we conducted genomic profiling of 957 unrelated individuals who self-identify with 28 large tribes in Saudi Arabia. Consistent with the tradition of intra-tribal unions, the subjects showed strong clustering along tribal lines with the distance between clusters correlating with their geographical proximities in Arabia. However, these individuals form a unique cluster when compared to the world's major populations. The ancient origin of these tribal affiliations is supported by analyses that revealed little evidence of ancestral origin from within the 28 tribes. Our results disclose a granular map of population structure and have important implications for future genetic studies into Mendelian and common diseases in the region.

Exploring Antiviral Drugs on Monolayer Black Phosphorene: Atomistic Theory and Explainable Machine Learning-Assisted Platform.

Favipiravir (FP) and ebselen (EB) belong to a diverse class of antiviral drugs known for their significant efficacy in treating various viral infections. Utilizing molecular dynamics (MD) simulations, machine learning, and van der Waals density functional theory, we accurately elucidate the binding properties of these antiviral drugs on a phosphorene single-layer. To further investigate these characteristics, this study employs four distinct machine learning models-Random Forest, Gradient Boosting, XGBoost, and CatBoost. The Hamiltonian of antiviral molecules within a monolayer of phosphorene is appropriately trained. The key aspect of utilizing machine learning (ML) in drug design revolves around training models that are efficient and precise in approximating density functional theory (DFT). Furthermore, the study employs SHAP (SHapley Additive exPlanations) to elucidate model predictions, providing insights into the contribution of each feature. To explore the interaction characteristics and thermodynamic properties of the hybrid drug, we employ molecular dynamics and DFT calculations in a vacuum interface. Our findings suggest that this functionalized 2D complex exhibits robust thermostability, indicating its potential as an effective and enabled entity. The observed variations in free energy at different surface charges and temperatures suggest the adsorption potential of FP and EB molecules from the surrounding environment.

DeepSVP: integration of genotype and phenotype for structural variant prioritization using deep learning.

MOTIVATION: Structural genomic variants account for much of human variability and are involved in several diseases. Structural variants are complex and may affect coding regions of multiple genes, or affect the functions of genomic regions in different ways from single nucleotide variants. Interpreting the phenotypic consequences of structural variants relies on information about gene functions, haploinsufficiency or triplosensitivity and other genomic features. Phenotype-based methods to identifying variants that are involved in genetic diseases combine molecular features with prior knowledge about the phenotypic consequences of altering gene functions. While phenotype-based methods have been applied successfully to single nucleotide variants as well as short insertions and deletions, the complexity of structural variants makes it more challenging to link them to phenotypes. Furthermore, structural variants can affect a large number of coding regions, and phenotype information may not be available for all of them. RESULTS: We developed DeepSVP, a computational method to prioritize structural variants involved in genetic diseases by combining genomic and gene functions information. We incorporate phenotypes linked to genes, functions of gene products, gene expression in individual cell types and anatomical sites of expression, and systematically relate them to their phenotypic consequences through ontologies and machine learning. DeepSVP significantly improves the success rate of finding causative variants in several benchmarks and can identify novel pathogenic structural variants in consanguineous families. AVAILABILITY AND IMPLEMENTATION: https://github.com/bio-ontology-research-group/DeepSVP. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

From shallow to deep: some lessons learned from application of machine learning for recognition of functional genomic elements in human genome.

Identification of genomic signals as indicators for functional genomic elements is one of the areas that received early and widespread application of machine learning methods. With time, the methods applied grew in variety and generally exhibited a tendency to improve their ability to identify some major genomic and transcriptomics signals. The evolution of machine learning in genomics followed a similar path to applications of machine learning in other fields. These were impacted in a major way by three dominant developments, namely an enormous increase in availability and quality of data, a significant increase in computational power available to machine learning applications, and finally, new machine learning paradigms, of which deep learning is the most well-known example. It is not easy in general to distinguish factors leading to improvements in results of applications of machine learning. This is even more so in the field of genomics, where the advent of next-generation sequencing and the increased ability to perform functional analysis of raw data have had a major effect on the applicability of machine learning in OMICS fields. In this paper, we survey the results from a subset of published work in application of machine learning in the recognition of genomic signals and regions in human genome and summarize some lessons learnt from this endeavor. There is no doubt that a significant progress has been made both in terms of accuracy and reliability of models. Questions remain however whether the progress has been sufficient and what these developments bring to the field of genomics in general and human genomics in particular. Improving usability, interpretability and accuracy of models remains an important open challenge for current and future research in application of machine learning and more generally of artificial intelligence methods in genomics.

The genome of Chenopodium quinoa.

Chenopodium quinoa (quinoa) is a highly nutritious grain identified as an important crop to improve world food security. Unfortunately, few resources are available to facilitate its genetic improvement. Here we report the assembly of a high-quality, chromosome-scale reference genome sequence for quinoa, which was produced using single-molecule real-time sequencing in combination with optical, chromosome-contact and genetic maps. We also report the sequencing of two diploids from the ancestral gene pools of quinoa, which enables the identification of sub-genomes in quinoa, and reduced-coverage genome sequences for 22 other samples of the allotetraploid goosefoot complex. The genome sequence facilitated the identification of the transcription factor likely to control the production of anti-nutritional triterpenoid saponins found in quinoa seeds, including a mutation that appears to cause alternative splicing and a premature stop codon in sweet quinoa strains. These genomic resources are an important first step towards the genetic improvement of quinoa.

Corrigendum: The genome of Chenopodium quinoa.

This corrects the article DOI: 10.1038/nature21370.

Glutamate-Sensing Genes Are Conserved among Populations Compared to Glutamate Metabolism Genes.

INTRODUCTION: Glutamate is a representative taste molecule with an umami flavor and is a major nutrient found abundantly in nature. Furthermore, it plays a significant role in the human body as a key metabolic intermediate and neurotransmitter. Therefore, the divergence of glutamate functions among populations during their evolution is of particular interest with a hypothesis that the genetic variation can lead to understanding divergence in taste perception. To elucidate variation in glutamate applications and to deepen our understanding of taste perception, we examined the nucleotide diversity of genes associated with glutamate sensing and metabolism among human populations. METHODS: We first established 67 genes related to glutamate sensing and metabolism based on the database and literature survey. Then, for those genes, we used a population genomics approach based on ten populations over 76,156 human genomes in the gnomAD database. RESULTS: Statistical tests of means and medians of the minor allele frequencies did not show any significant difference among populations. However, we observed substantial differences between two functional groups, glutamate sensing and glutamate metabolism, in populations of Latino/admixed American, Ashkenazi Jewish, and Others. Interestingly, we could find significant differences between the African population and the East Asian population at the single nucleotide polymorphism level of glutamate metabolism genes, but no clear differences were noted in glutamate-sensing genes. These suggest that glutamate-sensing genes are under the functional constraint compared to glutamate metabolism genes. CONCLUSION: Thus, glutamate-sensing genes and metabolism genes have a contrasting mode of the evolution, and glutamate-sensing genes are conservatively evolved, indicating its functional importance.

Behavioral and brain- transcriptomic synchronization between the two opponents of a fighting pair of the fish Betta splendens.

Conspecific male animals fight for resources such as food and mating opportunities but typically stop fighting after assessing their relative fighting abilities to avoid serious injuries. Physiologically, how the fighting behavior is controlled remains unknown. Using the fighting fish Betta splendens, we studied behavioral and brain-transcriptomic changes during the fight between the two opponents. At the behavioral level, surface-breathing, and biting/striking occurred only during intervals between mouth-locking. Eventually, the behaviors of the two opponents became synchronized, with each pair showing a unique behavioral pattern. At the physiological level, we examined the expression patterns of 23,306 brain transcripts using RNA-sequencing data from brains of fighting pairs after a 20-min (D20) and a 60-min (D60) fight. The two opponents in each D60 fighting pair showed a strong gene expression correlation, whereas those in D20 fighting pairs showed a weak correlation. Moreover, each fighting pair in the D60 group showed pair-specific gene expression patterns in a grade of membership analysis (GoM) and were grouped as a pair in the heatmap clustering. The observed pair-specific individualization in brain-transcriptomic synchronization (PIBS) suggested that this synchronization provides a physiological basis for the behavioral synchronization. An analysis using the synchronized genes in fighting pairs of the D60 group found genes enriched for ion transport, synaptic function, and learning and memory. Brain-transcriptomic synchronization could be a general phenomenon and may provide a new cornerstone with which to investigate coordinating and sustaining social interactions between two interacting partners of vertebrates.

Computational Studies of Auto-Active van der Waals Interaction Molecules on Ultra-Thin Black-Phosphorus Film.

Using the van der Waals density functional theory, we studied the binding peculiarities of favipiravir (FP) and ebselen (EB) molecules on a monolayer of black phosphorene (BP). We systematically examined the interaction characteristics and thermodynamic properties in a vacuum and a continuum, solvent interface for active drug therapy. These results illustrate that the hybrid molecules are enabled functionalized two-dimensional (2D) complex systems with a vigorous thermostability. We demonstrate in this study that these molecules remain flat on the monolayer BP system and phosphorus atoms are intact. It is inferred that the hybrid FP+EB molecules show larger adsorption energy due to the van der Waals forces and planar electrostatic interactions. The changes in Gibbs free energy at different surface charge fluctuations and temperatures imply that the FP and EB are allowed to adsorb from the gas phase onto the 2D film at high temperatures. Thereby, the results unveiled beneficial inhibitor molecules on two dimensional BP nanocarriers, potentially introducing a modern strategy to enhance the development of advanced materials, biotechnology, and nanomedicine.

Synergy of Small Antiviral Molecules on a Black-Phosphorus Nanocarrier: Machine Learning and Quantum Chemical Simulation Insights.

Favipiravir (FP) and Ebselen (EB) belong to a broad range of antiviral drugs that have shown active potential as medications against many viruses. Employing molecular dynamics simulations and machine learning (ML) combined with van der Waals density functional theory, we have uncovered the binding characteristics of these two antiviral drugs on a phosphorene nanocarrier. Herein, by using four different machine learning models (i.e., Bagged Trees, Gaussian Process Regression (GPR), Support Vector Regression (SVR), and Regression Trees (RT)), the Hamiltonian and the interaction energy of antiviral molecules in a phosphorene monolayer are trained in an appropriate way. However, training efficient and accurate models for approximating the density functional theory (DFT) is the final step in using ML to aid in the design of new drugs. To improve the prediction accuracy, the Bayesian optimization approach has been employed to optimize the GPR, SVR, RT, and BT models. Results revealed that the GPR model obtained superior prediction performance with an R2 of 0.9649, indicating that it can explain 96.49% of the data's variability. Then, by means of DFT calculations, we examine the interaction characteristics and thermodynamic properties in a vacuum and a continuum solvent interface. These results illustrate that the hybrid drug is an enabled, functionalized 2D complex with vigorous thermostability. The change in Gibbs free energy at different surface charges and temperatures implies that the FP and EB molecules are allowed to adsorb from the gas phase onto the 2D monolayer at different pH conditions and high temperatures. The results reveal a valuable antiviral drug therapy loaded by 2D biomaterials that may possibly open a new way of auto-treating different diseases, such as SARS-CoV, in primary terms.

A cautionary signal from the Red Sea on the impact of increased dust activity on marine microbiota.

BACKGROUND: Global climate change together with growing desertification is leading to increased dust emissions to the atmosphere, drawing attention to possible impacts on marine ecosystems receiving dust deposition. Since microorganisms play important roles in maintaining marine homeostasis through nutrient cycling and carbon flow, detrimental changes in the composition of marine microbiota in response to increased dust input could negatively impact marine health, particularly so in seas located within the Global Dust Belt. Due to its strategic location between two deserts and unique characteristics, the Red Sea provides an attractive semi-enclosed "megacosm" to examine the impacts of large dust deposition on the vastly diverse microbiota in its exceptionally warm oligotrophic waters. RESULTS: We used culture-independent metagenomic approaches to assess temporal changes in the Red Sea microbiota in response to two severe sandstorms, one originated in the Nubian Desert in the summer 2016 and a second one originated in the Libyan Desert in the spring 2017. Despite differences in sandstorm origin and meteorological conditions, both sandstorms shifted bacterial and Archaeal groups in a similar mode. In particular, the relative abundance of autotrophic bacteria declined while those of heterotrophic bacteria, particularly Bacteroidetes, and Archaea increased. The changes peaked within six days from the start of sandstorms, and the community recovered the original assemblage within one month. CONCLUSION: Our results suggest that increased dust emission with expanding desertification could lead to undesirable impacts in ocean function, enhancing heterotrophic processes while reducing autotrophic ones, thereby affecting the marine food web in seas receiving dust deposition.

Role of C-Reactive Protein in Diabetic Inflammation.

Even though type 2 diabetes mellitus (T2DM) represents a worldwide chronic health issue that affects about 462 million people, specific underlying determinants of insulin resistance (IR) and impaired insulin secretion are still unknown. There is growing evidence that chronic subclinical inflammation is a triggering factor in the origin of T2DM. Increased C-reactive protein (CRP) levels have been linked to excess body weight since adipocytes produce tumor necrosis factor α (TNF-α) and interleukin 6 (IL-6), which are pivotal factors for CRP stimulation. Furthermore, it is known that hepatocytes produce relatively low rates of CRP in physiological conditions compared to T2DM patients, in which elevated levels of inflammatory markers are reported, including CRP. CRP also participates in endothelial dysfunction, the production of vasodilators, and vascular remodeling, and increased CRP level is closely associated with vascular system pathology and metabolic syndrome. In addition, insulin-based therapies may alter CRP levels in T2DM. Therefore, determining and clarifying the underlying CRP mechanism of T2DM is imperative for novel preventive and diagnostic procedures. Overall, CRP is one of the possible targets for T2DM progression and understanding the connection between insulin and inflammation may be helpful in clinical treatment and prevention approaches.

Genome sequencing and identification of cellulase genes in Bacillus paralicheniformis strains from the Red Sea.

BACKGROUND: Cellulolytic microorganisms are considered a key player in the degradation of plant biomass in various environments. These microorganisms can be isolated from various environments, such as soils, the insect gut, the mammalian rumen and oceans. The Red Sea exhibits a unique environment in terms of presenting a high seawater temperature, high salinity, low nutrient levels and high biodiversity. However, there is little information regarding cellulase genes in the Red Sea environment. This study aimed to examine whether the Red Sea can be a resource for the bioprospecting of microbial cellulases by isolating cellulase-producing microorganisms from the Red Sea environment and characterizing cellulase genes. RESULTS: Three bacterial strains were successfully isolated from the plankton fraction and the surface of seagrass. The isolated strains were identified as Bacillus paralicheniformis and showed strong cellulase activity. These results suggested that these three isolates secreted active cellulases. By whole genome sequencing, we found 10 cellulase genes from the three isolates. We compared the expression of these cellulase genes under cellulase-inducing and non-inducing conditions and found that most of the cellulase genes were generally upregulated during cellulolysis in the isolates. Our operon structure analysis also showed that cellulase genes form operons with genes involved in various kinds of cellular reactions, such as protein metabolism, which suggests the existence of crosstalk between cellulolysis and other metabolic pathways in the bacterial isolates. These results suggest that multiple cellulases are playing important roles in cellulolysis. CONCLUSIONS: Our study reports the isolation and characterization of cellulase-producing bacteria from the Red Sea. Our whole-genome sequencing classified our three isolates as Bacillus paralicheniformis, and we revealed the presence of ten cellulase orthologues in each of three isolates' genomes. Our comparative expression analysis also identified that most of the cellulase genes were upregulated under the inducing conditions in general. Although cellulases have been roughly classified into three enzyme groups of beta-glucosidase, endo-ß-1,4-glucanase and exoglucanase, these findings suggest the importance to consider microbial cellulolysis as a more complex reaction with various kinds of cellulase enzymes.

Eye-like ocelloids are built from different endosymbiotically acquired components.

Multicellularity is often considered a prerequisite for morphological complexity, as seen in the camera-type eyes found in several groups of animals. A notable exception exists in single-celled eukaryotes called dinoflagellates, some of which have an eye-like 'ocelloid' consisting of subcellular analogues to a cornea, lens, iris, and retina. These planktonic cells are uncultivated and rarely encountered in environmental samples, obscuring the function and evolutionary origin of the ocelloid. Here we show, using a combination of electron microscopy, tomography, isolated-organelle genomics, and single-cell genomics, that ocelloids are built from pre-existing organelles, including a cornea-like layer made of mitochondria and a retinal body made of anastomosing plastids. We find that the retinal body forms the central core of a network of peridinin-type plastids, which in dinoflagellates and their relatives originated through an ancient endosymbiosis with a red alga. As such, the ocelloid is a chimaeric structure, incorporating organelles with different endosymbiotic histories. The anatomical complexity of single-celled organisms may be limited by the components available for differentiation, but the ocelloid shows that pre-existing organelles can be assembled into a structure so complex that it was initially mistaken for a multicellular eye. Although mitochondria and plastids are acknowledged chiefly for their metabolic roles, they can also be building blocks for greater structural complexity.

Sequencing effort dictates gene discovery in marine microbial metagenomes.

Massive metagenomic sequencing combined with gene prediction methods were previously used to compile the gene catalogue of the ocean and host-associated microbes. Global expeditions conducted over the past 15 years have sampled the ocean to build a catalogue of genes from pelagic microbes. Here we undertook a large sequencing effort of a perturbed Red Sea plankton community to uncover that the rate of gene discovery increases continuously with sequencing effort, with no indication that the retrieved 2.83 million non-redundant (complete) genes predicted from the experiment represented a nearly complete inventory of the genes present in the sampled community (i.e., no evidence of saturation). The underlying reason is the Pareto-like distribution of the abundance of genes in the plankton community, resulting in a very long tail of millions of genes present at remarkably low abundances, which can only be retrieved through massive sequencing. Microbial metagenomic projects retrieve a variable number of unique genes per Tera base-pair (Tbp), with a median value of 14.7 million unique genes per Tbp sequenced across projects. The increase in the rate of gene discovery in microbial metagenomes with sequencing effort implies that there is ample room for new gene discovery in further ocean and holobiont sequencing studies.

Mining biosynthetic gene clusters in Virgibacillus genomes.

BACKGROUND: Biosynthetic gene clusters produce a wide range of metabolites with activities that are of interest to the pharmaceutical industry. Specific interest is shown towards those metabolites that exhibit antimicrobial activities against multidrug-resistant bacteria that have become a global health threat. Genera of the phylum Firmicutes are frequently identified as sources of such metabolites, but the biosynthetic potential of its Virgibacillus genus is not known. Here, we used comparative genomic analysis to determine whether Virgibacillus strains isolated from the Red Sea mangrove mud in Rabigh Harbor Lagoon, Saudi Arabia, may be an attractive source of such novel antimicrobial agents. RESULTS: A comparative genomics analysis based on Virgibacillus dokdonensis Bac330, Virgibacillus sp. Bac332 and Virgibacillus halodenitrificans Bac324 (isolated from the Red Sea) and six other previously reported Virgibacillus strains was performed. Orthology analysis was used to determine the core genomes as well as the accessory genome of the nine Virgibacillus strains. The analysis shows that the Red Sea strain Virgibacillus sp. Bac332 has the highest number of unique genes and genomic islands compared to other genomes included in this study. Focusing on biosynthetic gene clusters, we show how marine isolates, including those from the Red Sea, are more enriched with nonribosomal peptides compared to the other Virgibacillus species. We also found that most nonribosomal peptide synthases identified in the Virgibacillus strains are part of genomic regions that are potentially horizontally transferred. CONCLUSIONS: The Red Sea Virgibacillus strains have a large number of biosynthetic genes in clusters that are not assigned to known products, indicating significant potential for the discovery of novel bioactive compounds. Also, having more modular synthetase units suggests that these strains are good candidates for experimental characterization of previously identified bioactive compounds as well. Future efforts will be directed towards establishing the properties of the potentially novel compounds encoded by the Red Sea specific trans-AT PKS/NRPS cluster and the type III PKS/NRPS cluster.