Pangenomic landscapes shape performances of a synthetic genetic circuit across Stutzerimonas species.

Engineering identical genetic circuits into different species typically results in large differences in performance due to the unique cellular environmental context of each host, a phenomenon known as the "chassis-effect" or "context-dependency". A better understanding of how genomic and physiological contexts underpin the chassis-effect will improve biodesign strategies across diverse microorganisms. Here, we combined a pangenomic-based gene expression analysis with quantitative measurements of performance from an engineered genetic inverter device to uncover how genome structure and function relate to the observed chassis-effect across six closely related Stutzerimonas hosts. Our results reveal that genome architecture underpins divergent responses between our chosen non-model bacterial hosts to the engineered device. Specifically, differential expression of the core genome, gene clusters shared between all hosts, was found to be the main source of significant concordance to the observed differential genetic device performance, whereas specialty genes from respective accessory genomes were not significant. A data-driven investigation revealed that genes involved in denitrification and components of trans-membrane transporter proteins were among the most differentially expressed gene clusters between hosts in response to the genetic device. Our results show that the chassis-effect can be traced along differences among the most conserved genome-encoded functions and that these differences create a unique biodesign space among closely related species.IMPORTANCEContemporary synthetic biology endeavors often default to a handful of model organisms to host their engineered systems. Model organisms such as Escherichia coli serve as attractive hosts due to their tractability but do not necessarily provide the ideal environment to optimize performance. As more novel microbes are domesticated for use as biotechnology platforms, synthetic biologists are urged to explore the chassis-design space to optimize their systems and deliver on the promises of synthetic biology. The consequences of the chassis-effect will therefore only become more relevant as the field of biodesign grows. In our work, we demonstrate that the performance of a genetic device is highly dependent on the host environment it operates within, promoting the notion that the chassis can be considered a design variable to tune circuit function. Importantly, our results unveil that the chassis-effect can be traced along similarities in genome architecture, specifically the shared core genome. Our study advocates for the exploration of the chassis-design space and is a step forward to empowering synthetic biologists with knowledge for more efficient exploration of the chassis-design space to enable the next generation of broad-host-range synthetic biology.

Genome wide identification of phenylalanine ammonia-lyase (PAL) gene family in Cucumis sativus (cucumber) against abiotic stress.

Phenylalanine ammonia lyase (PAL) is a widely studied enzyme in plant biology due to its role in connecting primary metabolism to secondary phenylpropanoid metabolism, significantly influencing plant growth, development, and stress response. Although PAL genes have been extensively studied in various plant species but their exploration in cucumber has been limited. This study successfully identified 11 CsPAL genes in Cucumis sativus (cucumber). These CsPAL genes were categorized based on their conserved sequences revealing patterns through MEME analysis and multiple sequence alignment. Interestingly, cis-elements related to stress were found in the promoter regions of CsPAL genes, indicating their involvement in responding to abiotic stress. Furthermore, these gene's promoters contained components associated with light, development and hormone responsiveness. This suggests that they may have roles in hormone developmental processes. MicroRNAs were identified as a key regulators for the CsPAL genes, playing a crucial role in modulating their expression. This discovery underscores the complex regulatory network involved in the plant's response to various stress conditions. The influence of these microRNAs further highlights the complicated mechanisms that plants use to manage stress. Gene expression patterns were analyzed using RNA-seq data. The significant upregulation of CsPAL9 during HT3h (heat stress for 3 h) and the heightened upregulation of both CsPAL9 and CsPAL7 under HT6h (heat stress for 6 h) in the transcriptome study suggest a potential role for these genes in cucumber's tolerance to heat stress. This comprehensive investigation aims to enhance our understanding of the PAL gene family's versatility, offering valuable insights for advancements in cucumber genetics.

Biological and clinical relevance of correlated expression levels of coding and long noncoding RNAs in HPV16 positive cervical cancers.

Human papillomavirus (HPV) drives cervical cancer (CaCx) pathogenesis and viral oncoproteins jeopardize global gene expression in such cancers. In this study, our aim was to identify differentially expressed coding (DEcGs) and long noncoding RNA genes (DElncGs) specifically sense intronic and Natural Antisense Transcripts as they are located in the genic regions and may have a direct influence on the expression pattern of their neighbouring coding genes. We compared HPV16-positive CaCx patients (N = 44) with HPV-negative normal individuals (N = 34) by employing strand-specific RNA-seq and determined the relationships between DEcGs and DElncGs and their clinical implications. By performing Gene set enrichment and protein-protein interaction (PPI) analyses of DEcGs, we identified enrichment of processes crucial for abortive virus life cycle and cancer progression. The DEcGs formed 16 gene clusters which we identified through Molecular Complex Detection (MCODE) plugin of Cytoscape. All the gene clusters portrayed cancer-related functions. We recorded significantly correlated expression levels of 79 DElncGs with DEcGs at proximal genomic loci based on Pearson's Correlation coefficients. Of these gene pairs, 24 pairs portrayed significantly altered correlation coefficients among patients, compared to normal individuals. Of these, 6 DEcGs of 6 such gene pairs, belonged to 5 of the identified gene clusters, one of which was survival-associated. Out of the 24 correlated DEcG: DElncG pairs, we identified 3 pairs, where expression of both members was significantly associated with patient overall survival. The findings justify the cooperative roles of these gene pairs, in patient prognostication, thereby bearing immense potential for translation. Thus, elucidation of correlative strengths between paired DElncGs and DEcGs in patient and normal samples, could serve as a foundation for identification of therapeutic and prognostic targets of HPV16-positive CaCx.

Arbuscular mycorrhizal fungi enhance drought resistance in Bombax ceiba by regulating SOD family genes.

The physiological activity facilitated by arbuscular mycorrhizal fungi (AMF) contributes to plants' ability to tolerate drought. Nevertheless, it is unclear if AMF colonization affects the expression of genes in the host plant that encode antioxidant enzymes in the superoxide dismutase (SOD) family, which help alleviate drought stress in plants. Here, we conducted a pot trial to determine whether colonization by the AMF Rhizophagus irregularis improves drought resistance in Bombax ceiba. We comprehensively analyzed the SOD gene family and evaluated genome-wide expression patterns of SODs and SOD activity in AMF-colonized and non-mycorrhizal plants under simulated drought. We identified a total of 13 SODs in the genome of B. ceiba, including three FeSODs (BcFSDs), three MnSODs (BcMSDs), and seven Cu/ZnSODs (BcCSDs). Phylogenetic analysis based on binding domain revealed that SOD genes from B. ceiba and various other plant species can be divided into three separate groups, showing significant bootstrap values. Our examination of gene composition and patterns suggests that most BcSOD genes in these three subgroups are significantly conserved. Additionally, it was noted that hormones and stress-responsive cis-regulatory elements were found in all BcSOD promoters. Expression profiling by qRT-PCR demonstrated that AMF increased relative expression levels of Cu/Zn-SODs in both roots and shoots under drought stress, except for BcCSD3 in roots. Furthermore, AMF colonization increased the relative expression of BcMSD1a and BcMSD1b in roots, augmenting SOD activities and increasing ROS scavenging during drought. In general, this work offers molecular evidence in support of the beneficial effect of AMF colonization on drought tolerance in B. ceiba. It also elucidates the expression patterns of SOD genes, which will support efforts to optimize mycorrhizal seedling cultivation under stressful conditions.

Genome-wide identification of gene families related to miRNA biogenesis in Mangifera indica L. and their possible role during heat stress.

Mango is a popular tropical fruit that requires quarantine hot water treatment (QHWT) for postharvest sanitation, which can cause abiotic stress. Plants have various defense mechanisms to cope with stress; miRNAs mainly regulate the expression of these defense responses. Proteins involved in the biogenesis of miRNAs include DICER-like (DCL), ARGONAUTE (AGO), HYPONASTIC LEAVES 1 (HYL1), SERRATE (SE), HUA ENHANCER1 (HEN1), HASTY (HST), and HEAT-SHOCK PROTEIN 90 (HSP90), among others. According to our analysis, the mango genome contains five DCL, thirteen AGO, six HYL, two SE, one HEN1, one HST, and five putative HSP90 genes. Gene structure prediction and domain identification indicate that sequences contain key domains for their respective gene families, including the RNase III domain in DCL and PAZ and PIWI domains for AGOs. In addition, phylogenetic analysis indicates the formation of clades that include the mango sequences and their respective orthologs in other flowering plant species, supporting the idea these are functional orthologs. The analysis of cis-regulatory elements of these genes allowed the identification of MYB, ABRE, GARE, MYC, and MeJA-responsive elements involved in stress responses. Gene expression analysis showed that most genes are induced between 3 to 6 h after QHWT, supporting the early role of miRNAs in stress response. Interestingly, our results suggest that mango rapidly induces the production of miRNAs after heat stress. This research will enable us to investigate further the regulation of gene expression and its effects on commercially cultivated fruits, such as mango, while maintaining sanitary standards.

Identification of the RRM1 gene family in rice (Oryza sativa) and its response to rice blast.

To better understand RNA-binding proteins in rice, a comprehensive investigation was conducted on the RRM1 gene family of rice. It encompassed genome-wide identification and exploration of its role in rice blast resistance. The physicochemical properties of the rice OsRRM1 gene family were analyzed. There genes were also analyzed for their conserved domains, motifs, location information, gene structure, phylogenetic trees, collinearity, and cis-acting elements. Furthermore, alterations in the expression patterns of selected OsRRM1 genes were assessed using quantitative real-time PCR (qRT-PCR). A total of 212 members of the OsRRM1 gene family were identified, which were dispersed across 12 chromosomes. These genes all exhibit multiple exons and introns, all of which encompass the conserved RRM1 domain and share analogous motifs. This observation suggests a high degree of conservation within the encoded sequence domain of these genes. Phylogenetic analysis revealed the existence of five subfamilies within the OsRRM1 gene family. Furthermore, investigation of the promoter region identified cis-regulatory elements that are involved in nucleic acid binding and interaction with multiple transcription factors. By employing GO and KEGG analyses, four RRM1 genes were tentatively identified as crucial contributors to plant immunity, while the RRM1 gene family was also found to have a significant involvement in the complex of alternative splicing. The qRT-PCR results revealed distinct temporal changes in the expression patterns of OsRRM1 genes following rice blast infection. Additionally, gene expression analysis indicates that the majority of OsRRM1 genes exhibited constitutive expressions. These findings enrich our understanding of the OsRRM1 gene family. They also provide a foundation for further research on immune mechanisms rice and the management of rice blast.

A new synthetic biology system for investigating the biosynthesis of antibiotics and other secondary metabolites in streptomycetes.

We have created a novel synthetic biology expression system allowing easy refactoring of biosynthetic gene clusters (BGCs) as monocistronic transcriptional units. The system is based on a set of plasmids containing a strong kasOp* promoter, RBS and terminators. It allows the cloning of biosynthetic genes into transcriptional units kasOp*-gene(s)-terminator flanked by several rare restriction cloning sites that can be sequentially combined into the artificial BGC in three compatible Streptomyces integration vectors. They allow a simultaneous integration of these BGCs at three different attB sites in the Streptomyces chromosome. The system was validated with biosynthetic genes from two known BGCs for aromatic polyketides landomycin and mithramycin.

Metabologenomics reveals strain-level genetic and chemical diversity of Microcystis secondary metabolism.

Microcystis spp. are renowned for producing the hepatotoxin microcystin in freshwater cyanobacterial harmful algal blooms around the world, threatening drinking water supplies and public and environmental health. However, Microcystis genomes also harbor numerous biosynthetic gene clusters (BGCs) encoding the biosynthesis of other secondary metabolites, including many with toxic properties. Most of these BGCs are uncharacterized and currently lack links to biosynthesis products. However, recent field studies show that many of these BGCs are abundant and transcriptionally active in natural communities, suggesting potentially important yet unknown roles in bloom ecology and water quality. Here, we analyzed 21 xenic Microcystis cultures isolated from western Lake Erie to investigate the diversity of the biosynthetic potential of this genus. Through metabologenomic and in silico approaches, we show that these Microcystis strains contain variable BGCs, previously observed in natural populations, and encode distinct metabolomes across cultures. Additionally, we find that the majority of metabolites and gene clusters are uncharacterized, highlighting our limited understanding of the chemical repertoire of Microcystis spp. Due to the complex metabolomes observed in culture, which contain a wealth of diverse congeners as well as unknown metabolites, these results underscore the need to deeply explore and identify secondary metabolites produced by Microcystis beyond microcystins to assess their impacts on human and environmental health.IMPORTANCEThe genus Microcystis forms dense cyanobacterial harmful algal blooms (cyanoHABs) and can produce the toxin microcystin, which has been responsible for drinking water crises around the world. While microcystins are of great concern, Microcystis also produces an abundance of other secondary metabolites that may be of interest due to their potential for toxicity, ecological importance, or pharmaceutical applications. In this study, we combine genomic and metabolomic approaches to study the genes responsible for the biosynthesis of secondary metabolites as well as the chemical diversity of produced metabolites in Microcystis strains from the Western Lake Erie Culture Collection. This unique collection comprises Microcystis strains that were directly isolated from western Lake Erie, which experiences substantial cyanoHAB events annually and has had negative impacts on drinking water, tourism, and industry.

Genome-wide identification and expression analysis of the WRKY gene family in Rhododendron henanense subsp. lingbaoense.

Background: This work explored the characteristics of the WRKY transcription factor family in Rhododendron henanense subsp. lingbaoense (Rhl) and the expression patterns of these genes under abiotic stress by conducting bioinformatics and expression analyses. Methods: RhlWRKY genes were identified from a gene library of Rhl. Various aspects of these genes were analyzed, including genetic structures, conserved sequences, physicochemical properties, cis-acting elements, and chromosomal location. RNA-seq was employed to analyze gene expression in five different tissues of Rhl: roots, stems, leaves, flowers, and hypocotyls. Additionally, qRT-PCR was used to detect changes in the expression of five RhlWRKY genes under abiotic stress. Result: A total of 65 RhlWRKY genes were identified and categorized into three subfamilies based on their structural characteristics: Groups I, II, and III. Group II was further divided into five subtribes, with shared similar genetic structures and conserved motifs among members of the same subtribe. The physicochemical properties of these proteins varied, but the proteins are generally predicted to be hydrophilic. Most proteins are predicted to be in the cell nucleus, and distributed across 12 chromosomes. A total of 84 cis-acting elements were discovered, with many related to responses to biotic stress. Among the identified RhlWRKY genes, there were eight tandem duplicates and 97 segmental duplicates. The majority of duplicate gene pairs exhibited Ka/Ks values <1, indicating purification under environmental pressure. GO annotation analysis indicated that WRKY genes regulate biological processes and participate in a variety of molecular functions. Transcriptome data revealed varying expression levels of 66.15% of WRKY family genes in all five tissue types (roots, stems, leaves, flowers, and hypocotyls). Five RhlWRKY genes were selected for further characterization and there were changes in expression levels for these genes in response to various stresses. Conclusion: The analysis identified 65 RhlWRKY genes, among which the expression of WRKY_42 and WRKY_17 were mainly modulated by the drought and MeJA, and WRKY_19 was regulated by the low-temperature and high-salinity conditions. This insight into the potential functions of certain genes contributes to understanding the growth regulatory capabilities of Rhl.

Pangenome characterization and analysis of the NAC gene family reveals genes for Sclerotinia sclerotiorum resistance in sunflower (Helianthus annuus).

BACKGROUND: Sunflower (Helianthus annuus) is one of the most important economic crops in oilseed production worldwide. The different cultivars exhibit variability in their resistance genes. The NAC transcription factor (TF) family plays diverse roles in plant development and stress responses. With the completion of the H. annuus genome sequence, the entire complement of genes coding for NACs has been identified. However, the reference genome of a single individual cannot cover all the genetic information of the species. RESULTS: Considering only a single reference genome to study gene families will miss many meaningful genes. A pangenome-wide survey and characterization of the NAC genes in sunflower species were conducted. In total, 139 HaNAC genes are identified, of which 114 are core and 25 are variable. Phylogenetic analysis of sunflower NAC proteins categorizes these proteins into 16 subgroups. 138 HaNACs are randomly distributed on 17 chromosomes. SNP-based haplotype analysis shows haplotype diversity of the HaNAC genes in wild accessions is richer than in landraces and modern cultivars. Ten HaNAC genes in the basal stalk rot (BSR) resistance quantitative trait loci (QTL) are found. A total of 26 HaNAC genes are differentially expressed in response to Sclerotinia head rot (SHR). A total of 137 HaNAC genes are annotated in Gene Ontology (GO) and are classified into 24 functional groups. GO functional enrichment analysis reveals that HaNAC genes are involved in various functions of the biological process. CONCLUSIONS: We identified NAC genes in H. annuus (HaNAC) on a pangenome-wide scale and analyzed S. sclerotiorum resistance-related NACs. This study provided a theoretical basis for further genomic improvement targeting resistance-related NAC genes in sunflowers.

Genome­wide analysis of the MYB gene family in pumpkin.

The MYB gene family exerts significant influence over various biological processes and stress responses in plants. Despite this, a comprehensive analysis of this gene family in pumpkin remains absent. In this study, the MYB genes of Cucurbita moschata were identified and clustered into 33 groups (C1-33), with members of each group being highly conserved in terms of their motif composition. Furthermore, the distribution of 175 CmoMYB genes across all 20 chromosomes was found to be non-uniform. Examination of the promoter regions of these genes revealed the presence of cis-acting elements associated with phytohormone responses and abiotic/biotic stress. Utilizing quantitative real-time polymerase chain reaction (qRT-PCR), the expression patterns of 13 selected CmoMYB genes were validated, particularly in response to exogenous phytohormone exposure and various abiotic stressors, including ABA, SA, MeJA, and drought treatments. Expression analysis in different tissues showed that CmoMYB genes are expressed at different levels in different tissues, suggesting that they are functionally divergent in regulating growth and abiotic stresses. These results provide a basis for future studies to characterize the function of the MYB gene family under abiotic stresses in pumpkins.

Genome-wide identification and analysis of ascorbate peroxidase (APX) gene family in hemp (Cannabis sativa L.) under various abiotic stresses.

Ascorbate peroxidase (APX) plays a critical role in molecular mechanisms such as plant development and defense against abiotic stresses. As an important economic crop, hemp (Cannabis sativa L.) is vulnerable to adverse environmental conditions, such as drought, cold, salt, and oxidative stress, which lead to a decline in yield and quality. Although APX genes have been characterized in a variety of plants, members of the APX gene family in hemp have not been completely identified. In this study, we (1) identified eight members of the CsAPX gene family in hemp and mapped their locations on the chromosomes using bioinformatics analysis; (2) examined the physicochemical characteristics of the proteins encoded by these CsAPX gene family members; (3) investigated their intraspecific collinearity, gene structure, conserved domains, conserved motifs, and cis-acting elements; (4) constructed a phylogenetic tree and analyzed interspecific collinearity; and (5) ascertained expression differences in leaf tissue subjected to cold, drought, salt, and oxidative stresses using quantitative real-time-PCR (qRT-PCR). Under all four stresses, CsAPX6, CsAPX7, and CsAPX8 consistently exhibited significant upregulation, whereas CsAPX2 displayed notably higher expression levels under drought stress than under the other stresses. Taken together, the results of this study provide basic genomic information on the expression of the APX gene family and pave the way for studying the role of APX genes in abiotic stress.

Tools to make Stachybotrys chartarum genetically amendable: Key to unlocking cryptic biosynthetic gene clusters.

The soil and indoor fungus Stachybotrys chartarum can induce respiratory disorders, collectively referred to as stachybotryotoxicosis, owing to its prolific production of diverse bioactive secondary metabolites (SMs) or mycotoxins. Although many of these toxins responsible for the harmful effects on animals and humans have been identified in the genus Stachybotrys, however a number of SMs remain elusive. Through in silico analyses, we have identified 37 polyketide synthase (PKS) genes, highlighting that the chemical profile potential of Stachybotrys is far from being fully explored. Additionally, by leveraging phylogenetic analysis of known SMs produced by non-reducing polyketide synthases (NR-PKS) in other filamentous fungi, we showed that Stachybotrys possesses a rich reservoir of untapped SMs. To unravel natural product biosynthesis in S. chartarum, genetic engineering methods are crucial. For this purpose, we have developed a reliable protocol for the genetic transformation of S. chartarum and applied it to the ScPKS14 biosynthetic gene cluster. This cluster is homologous to the already known Claviceps purpurea CpPKS8 BGC, responsible for the production of ergochromes. While no novel SMs were detected, we successfully applied genetic tools, such as the generation of deletionand overexpression strains of single cluster genes. This toolbox can now be readily employed to unravel not only this particular BGC but also other candidate BGCs present in S. chartarum, making this fungus accessible for genetic engineering.

Oxytetracycline hyper-production through targeted genome reduction of Streptomyces rimosus.

Most biosynthetic gene clusters (BGC) encoding the synthesis of important microbial secondary metabolites, such as antibiotics, are either silent or poorly expressed; therefore, to ensure a strong pipeline of novel antibiotics, there is a need to develop rapid and efficient strain development approaches. This study uses comparative genome analysis to instruct rational strain improvement, using Streptomyces rimosus, the producer of the important antibiotic oxytetracycline (OTC) as a model system. Sequencing of the genomes of two industrial strains M4018 and R6-500, developed independently from a common ancestor, identified large DNA rearrangements located at the chromosome end. We evaluated the effect of these genome deletions on the parental S. rimosus Type Strain (ATCC 10970) genome where introduction of a 145 kb deletion close to the OTC BGC in the Type Strain resulted in massive OTC overproduction, achieving titers that were equivalent to M4018 and R6-500. Transcriptome data supported the hypothesis that the reason for such an increase in OTC biosynthesis was due to enhanced transcription of the OTC BGC and not due to enhanced substrate supply. We also observed changes in the expression of other cryptic BGCs; some metabolites, undetectable in ATCC 10970, were now produced at high titers. This study demonstrated for the first time that the main force behind BGC overexpression is genome rearrangement. This new approach demonstrates great potential to activate cryptic gene clusters of yet unexplored natural products of medical and industrial value.IMPORTANCEThere is a critical need to develop novel antibiotics to combat antimicrobial resistance. Streptomyces species are very rich source of antibiotics, typically encoding 20-60 biosynthetic gene clusters (BGCs). However, under laboratory conditions, most are either silent or poorly expressed so that their products are only detectable at nanogram quantities, which hampers drug development efforts. To address this subject, we used comparative genome analysis of industrial Streptomyces rimosus strains producing high titers of a broad spectrum antibiotic oxytetracycline (OTC), developed during decades of industrial strain improvement. Interestingly, large-scale chromosomal deletions were observed. Based on this information, we carried out targeted genome deletions in the native strain S. rimosus ATCC 10970, and we show that a targeted deletion in the vicinity of the OTC BGC significantly induced expression of the OTC BGC, as well as some other silent BGCs, thus suggesting that this approach may be a useful way to identify new natural products.

Genome-wide identification of the Hsp70 gene family in Penaeus chinensis and their response to environmental stress.

The heat shock protein 70 (HSP70) gene family plays a crucial role in the response of organisms to environmental stress. However, it has not been systematically characterized in shrimp. In this study, we identified 25 PcHsp70 genes in the Penaeus chinensis genome. The encoded proteins were categorized into six subgroups based on phylogenetic relationships. Tandem duplication was the main driver of amplification in the PcHsp70 family, and the genes have experienced strong purifying selection during evolution. Transcriptome data analysis revealed that the 25 PcHsp70 members have different expression patterns in shrimp under conditions of low temperature, low salinity, and white spot syndrome virus infection. Among them, PcHsp70.11 was significantly induced under all three stress conditions, suggesting that this gene plays an important role in response to environmental stress in P. chinensis. To the best of our knowledge, this is the first study to systematically analyze the Hsp70 gene family in shrimp. The results provide important information on shrimp Hsp70s, contributing to a better understanding of the role of these genes in environmental stress and providing a basis for further functional studies.

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.

Decrypting biocontrol functions and application modes by genomes data of three Trichoderma Strains/Species.

Trichoderma is an excellent biocontrol agent, but most Trichoderma genomes remained at the scaffold level, which greatly limits the research of biocontrol mechanism. Here, we reported the chromosome-level genome of Trichoderma harzianum CGMCC20739 (Tha739), T. asperellum CGMCC11653 (Tas653) and T. atroviride CGMCC40488 (Tat488), they were assembled into 7 chromosomes, genome size were 40 Mb (10,611 genes), 37.3 Mb (10,102 genes) and 36.3 Mb (9,896 genes), respectively. The positive selected genes of three strains were associated to response to stimulus, signaling transduction, immune system and localization. Furthermore, the number of transcription factors in Tha739, Tas653 and Tat488 strains had significant difference, which may contribute to the differential biocontrol function and stress tolerance. The genes related to signal transduction and gene clusters related to antimicrobial compounds in Tha739 were more than those in Tas653 and Tat488, which showed Tha739 may keenly sense other fungi and quickly secret antimicrobial compounds to inhibit other fungi. Tha739 also contained more genes associated to detoxification, antioxidant and nutrition utilization, indicating it had higher stress-tolerance to hostile environments. And the substrate for synthesizing IAA in Tha739 was mainly 3-indole acetonitrile and indole acetaldehyde, but in Tat488, it was indole-3-acetamide, moreover, Tha739 secreted more phosphatase and phytase and was more related to soil phosphorus metabolism, Tat488 secreted more urease and was more related to soil nitrogen metabolism. These candidate genes related to biocontrol function and stress-tolerance laid foundations for construction of functional strains. All above proved the difference in biocontrol function of Tha739, Tas653 and Tat488 strains, however, the defects in individual strains could be compensated for through Trichoderma-biome during the commercial application process of biocontrol Trichoderma strains.

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.

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.