Several secondary metabolite gene clusters in the genomes of ten Penicillium spp. raise the risk of multiple mycotoxin occurrence in chestnuts.

Penicillium spp. produce a great variety of secondary metabolites, including several mycotoxins, on food substrates. Chestnuts represent a favorable substrate for Penicillium spp. development. In this study, the genomes of ten Penicillium species, virulent on chestnuts, were sequenced and annotated: P. bialowiezense. P. pancosmium, P. manginii, P. discolor, P. crustosum, P. palitans, P. viridicatum, P. glandicola, P. taurinense and P. terrarumae. Assembly size ranges from 27.5 to 36.8 Mb and the number of encoded genes ranges from 9,867 to 12,520. The total number of predicted biosynthetic gene clusters (BGCs) in the ten species is 551. The most represented families of BGCs are non ribosomal peptide synthase (191) and polyketide synthase (175), followed by terpene synthases (87). Genome-wide collections of gene phylogenies (phylomes) were reconstructed for each of the newly sequenced Penicillium species allowing for the prediction of orthologous relationships among our species, as well as other 20 annotated Penicillium species available in the public domain. We investigated in silico the presence of BGCs for 10 secondary metabolites, including 5 mycotoxins, whose production was validated in vivo through chemical analyses. Among the clusters present in this set of species we found andrastin A and its related cluster atlantinone A, mycophenolic acid, patulin, penitrem A and the cluster responsible for the synthesis of roquefortine C/glandicoline A/glandicoline B/meleagrin. We confirmed the presence of these clusters in several of the Penicillium species conforming our dataset and verified their capacity to synthesize them in a chestnut-based medium with chemical analysis. Interestingly, we identified mycotoxin clusters in some species for the first time, such as the andrastin A cluster in P. flavigenum and P. taurinense, and the roquefortine C cluster in P. nalgiovense and P. taurinense. Chestnuts proved to be an optimal substrate for species of Penicillium with different mycotoxigenic potential, opening the door to risks related to the occurrence of multiple mycotoxins in the same food matrix.

Genome-wide identification and characterization of DTX family genes highlighting their locations, functions, and regulatory factors in banana (Musa acuminata).

The detoxification efflux carriers (DTX) are a significant group of multidrug efflux transporter family members that play diverse functions in all kingdoms of living organisms. However, genome-wide identification and characterization of DTX family transporters have not yet been performed in banana, despite its importance as an economic fruit plant. Therefore, a detailed genome-wide analysis of DTX family transporters in banana (Musa acuminata) was conducted using integrated bioinformatics and systems biology approaches. In this study, a total of 37 DTX transporters were identified in the banana genome and divided into four groups (I, II, III, and IV) based on phylogenetic analysis. The gene structures, as well as their proteins' domains and motifs, were found to be significantly conserved. Gene ontology (GO) annotation revealed that the predicted DTX genes might play a vital role in protecting cells and membrane-bound organelles through detoxification mechanisms and the removal of drug molecules from banana cells. Gene regulatory analyses identified key transcription factors (TFs), cis-acting elements, and post-transcriptional regulators (miRNAs) of DTX genes, suggesting their potential roles in banana. Furthermore, the changes in gene expression levels due to pathogenic infections and non-living factor indicate that banana DTX genes play a role in responses to both biotic and abiotic stresses. The results of this study could serve as valuable tools to improve banana quality by protecting them from a range of environmental stresses.

Cis-regulatory evolution of the recently expanded Ly49 gene family.

Comparative genomics has revealed the rapid expansion of multiple gene families involved in immunity. Members within each gene family often evolved distinct roles in immunity. However, less is known about the evolution of their epigenome and cis-regulation. Here we systematically profile the epigenome of the recently expanded murine Ly49 gene family that mainly encode either inhibitory or activating surface receptors on natural killer cells. We identify a set of cis-regulatory elements (CREs) for activating Ly49 genes. In addition, we show that in mice, inhibitory and activating Ly49 genes are regulated by two separate sets of proximal CREs, likely resulting from lineage-specific losses of CRE activity. Furthermore, we find that some Ly49 genes are cross-regulated by the CREs of other Ly49 genes, suggesting that the Ly49 family has begun to evolve a concerted cis-regulatory mechanism. Collectively, we demonstrate the different modes of cis-regulatory evolution for a rapidly expanding gene family.

Genome-wide identification and expression analysis of the KNOX family and its diverse roles in response to growth and abiotic tolerance in sweet potato and its two diploid relatives.

KNOXs, a type of homeobox genes that encode atypical homeobox proteins, play an essential role in the regulation of growth and development, hormonal response, and abiotic stress in plants. However, the KNOX gene family has not been explored in sweet potato. In this study, through sequence alignment, genomic structure analysis, and phylogenetic characterization, 17, 12 and 11 KNOXs in sweet potato (I. batatas, 2n = 6x = 90) and its two diploid relatives I. trifida (2n = 2x = 30) and I. triloba (2n = 2x = 30) were identified. The protein physicochemical properties, chromosome localization, phylogenetic relationships, gene structure, protein interaction network, cis-elements of promoters, tissue-specific expression and expression patterns under hormone treatment and abiotic stresses of these 40 KNOX genes were systematically studied. IbKNOX4, -5, and - 6 were highly expressed in the leaves of the high-yield varieties Longshu9 and Xushu18. IbKNOX3 and IbKNOX8 in Class I were upregulated in initial storage roots compared to fibrous roots. IbKNOXs in Class M were specifically expressed in the stem tip and hardly expressed in other tissues. Moreover, IbKNOX2 and - 6, and their homologous genes were induced by PEG/mannitol and NaCl treatments. The results showed that KNOXs were involved in regulating growth and development, hormone crosstalk and abiotic stress responses between sweet potato and its two diploid relatives. This study provides a comparison of these KNOX genes in sweet potato and its two diploid relatives and a theoretical basis for functional studies.

Genome mining of Lactiplantibacillus plantarum PA21: insights into its antimicrobial potential.

BACKGROUND: The dramatic increase of antimicrobial resistance in the healthcare realm has become inexorably linked to the abuse of antibiotics over the years. Therefore, this study seeks to identify potential postbiotic metabolites derived from lactic acid bacteria such as Lactiplantibacillus plantarum that could exhibit antimicrobial properties against multi-drug resistant pathogens. RESULTS: In the present work, the genome sequence of Lactiplantibacillus plantarum PA21 consisting of three contigs was assembled to a size of 3,218,706 bp. Phylogenomic analysis and average nucleotide identity (ANI) revealed L. plantarum PA21 is closely related to genomes isolated from diverse niches such as dairy products, food, and animals. Genome mining through the BAGEL4 and antiSMASH database revealed four bacteriocins in a single cluster and four regions of biosynthetic gene clusters responsible for the production of bioactive compounds. The potential probiotic genes indirectly responsible for postbiotic metabolites production were also identified. Additionally, in vitro studies showed that the L. plantarum PA21 cell-free supernatant exhibited antimicrobial activity against all nine methicillin-resistant Staphylococcus aureus (MRSA) and three out of 13 Klebsiella pneumoniae clinical isolates tested. CONCLUSION: Results in this study demonstrates that L. plantarum PA21 postbiotic metabolites is a prolific source of antimicrobials against multi-drug resistant pathogens with potential antimicrobial properties.

dgfr: an R package to assess sequence diversity of gene families.

BACKGROUND: Gene families are groups of homologous genes that often have similar biological functions. These families are formed by gene duplication events throughout evolution, resulting in multiple copies of an ancestral gene. Over time, these copies can acquire mutations and structural variations, resulting in members that may vary in size, motif ordering and sequence. Multigene families have been described in a broad range of organisms, from single-celled bacteria to complex multicellular organisms, and have been linked to an array of phenomena, such as host-pathogen interactions, immune evasion and embryonic development. Despite the importance of gene families, few approaches have been developed for estimating and graphically visualizing their diversity patterns and expression profiles in genome-wide studies. RESULTS: Here, we introduce an R package named dgfr, which estimates and enables the visualization of sequence divergence within gene families, as well as the visualization of secondary data such as gene expression. The package takes as input a multi-fasta file containing the coding sequences (CDS) or amino acid sequences from a multigene family, performs a pairwise alignment among all sequences, and estimates their distance, which is subjected to dimension reduction, optimal cluster determination, and gene assignment to each cluster. The result is a dataset that allows for the visualization of sequence divergence and expression within the gene family, an approximation of the number of clusters present in the family. CONCLUSIONS: dgfr provides a way to estimate and study the diversity of gene families, as well as visualize the dispersion and secondary profile of the sequences. The dgfr package is available at https://github.com/lailaviana/dgfr under the GPL-3 license.

Profiling expression strategies for a type III polyketide synthase in a lysate-based, cell-free system.

Some of the most metabolically diverse species of bacteria (e.g., Actinobacteria) have higher GC content in their DNA, differ substantially in codon usage, and have distinct protein folding environments compared to tractable expression hosts like Escherichia coli. Consequentially, expressing biosynthetic gene clusters (BGCs) from these bacteria in E. coli often results in a myriad of unpredictable issues with regard to protein expression and folding, delaying the biochemical characterization of new natural products. Current strategies to achieve soluble, active expression of these enzymes in tractable hosts can be a lengthy trial-and-error process. Cell-free expression (CFE) has emerged as a valuable expression platform as a testbed for rapid prototyping expression parameters. Here, we use a type III polyketide synthase from Streptomyces griseus, RppA, which catalyzes the formation of the red pigment flaviolin, as a reporter to investigate BGC refactoring techniques. We applied a library of constructs with different combinations of promoters and rppA coding sequences to investigate the synergies between promoter and codon usage. Subsequently, we assess the utility of cell-free systems for prototyping these refactoring tactics prior to their implementation in cells. Overall, codon harmonization improves natural product synthesis more than traditional codon optimization across cell-free and cellular environments. More importantly, the choice of coding sequences and promoters impact protein expression synergistically, which should be considered for future efforts to use CFE for high-yield protein expression. The promoter strategy when applied to RppA was not completely correlated with that observed with GFP, indicating that different promoter strategies should be applied for different proteins. In vivo experiments suggest that there is correlation, but not complete alignment between expressing in cell free and in vivo. Refactoring promoters and/or coding sequences via CFE can be a valuable strategy to rapidly screen for catalytically functional production of enzymes from BCGs, which advances CFE as a tool for natural product research.

Genome-wide characterization, functional analysis, and expression profiling of the Aux/IAA gene family in spinach.

BACKGROUND: The auxin/indole-3-acetic acid (Aux/IAA) gene family is a crucial element of the auxin signaling pathway, significantly influencing plant growth and development. Hence, we conducted a comprehensive investigation of Aux/IAAs gene family using the Sp75 and Monoe-Viroflay genomes in spinach. RESULTS: A total of 24 definitive Aux/IAA genes were identified, exhibiting diverse attributes in terms of amino acid length, molecular weight, and isoelectric points. This diversity underscores potential specific roles within the family, such as growth regulation and stress response. Structural analysis revealed significant variations in gene length and molecular weight. These variations indicate distinct roles within the Aux/IAA gene family. Chromosomal distribution analysis exhibited a dispersed pattern, with chromosomes 4 and 1 hosting the highest and lowest numbers of Aux/IAA genes, respectively. Phylogenetic analysis grouped the identified genes into distinct clades, revealing potential evolutionary relationships. Notably, the phylogenetic tree highlighted specific gene clusters suggesting shared genetic ancestry and potential functional synergies within spinach. Expression analysis under NAA treatment unveiled gene-specific and time-dependent responses, with certain genes exhibiting distinct temporal expression patterns. Specifically, SpoIAA5 displayed a substantial increase at 2 h post-NAA treatment, while SpoIAA7 and SpoIAA9 demonstrated continuous rises, peaking at the 4-hour time point. CONCLUSIONS: These observations indicate a complex interplay of gene-specific and temporal regulation in response to auxin. Moreover, the comparison with other plant species emphasized both shared characteristics and unique features in Aux/IAA gene numbers, providing insights into the evolutionary dynamics of this gene family. This comprehensive characterization of Aux/IAA genes in spinach not only establishes the foundation for understanding their specific functions in spinach development but also provides a valuable resource for experimental validation and further exploration of their roles in the intricate network of auxin signaling pathways.

Genome-wide identification and characterization of SPXdomain-containing genes family in eggplant.

Phosphorus is one of the lowest elements absorbed and utilized by plants in the soil. SPX domain-containing genes family play an important role in plant response to phosphate deficiency signaling pathway, and related to seed development, disease resistance, absorption and transport of other nutrients. However, there are no reports on the mechanism of SPX domain-containing genes in response to phosphorus deficiency in eggplant. In this study, the whole genome identification and functional analysis of SPX domain-containing genes family in eggplant were carried out. Sixteen eggplant SPX domain-containing genes were identified and divided into four categories. Subcellular localization showed that these proteins were located in different cell compartments, including nucleus and membrane system. The expression patterns of these genes in different tissues as well as under phosphate deficiency with auxin were explored. The results showed that SmSPX1, SmSPX5 and SmSPX12 were highest expressed in roots. SmSPX1, SmSPX4, SmSPX5 and SmSPX14 were significantly induced by phosphate deficiency and may be the key candidate genes in response to phosphate starvation in eggplant. Among them, SmSPX1 and SmSPX5 can be induced by auxin under phosphate deficiency. In conclusion, our study preliminary identified the SPX domain genes in eggplant, and the relationship between SPX domain-containing genes and auxin was first analyzed in response to phosphate deficiency, which will provide theoretical basis for improving the absorption of phosphorus in eggplants through molecular breeding technology.

Genomic insights into CKX genes: key players in cotton fibre development and abiotic stress responses.

Cytokinin oxidase/dehydrogenase (CKX), responsible for irreversible cytokinin degradation, also controls plant growth and development and response to abiotic stress. While the CKX gene has been studied in other plants extensively, its function in cotton is still unknown. Therefore, a genome-wide study to identify the CKX gene family in the four cotton species was conducted using transcriptomics, quantitative real-time PCR (qRT-PCR) and bioinformatics. As a result, in G. hirsutum and G. barbadense (the tetraploid cotton species), 87 and 96 CKX genes respectively and 62 genes each in G. arboreum and G. raimondii, were identified. Based on the evolutionary studies, the cotton CKX gene family has been divided into five distinct subfamilies. It was observed that CKX genes in cotton have conserved sequence logos and gene family expansion was due to segmental duplication or whole genome duplication (WGD). Collinearity and multiple synteny studies showed an expansion of gene families during evolution and purifying selection pressure has been exerted. G. hirsutum CKX genes displayed multiple exons/introns, uneven chromosomal distribution, conserved protein motifs, and cis-elements related to growth and stress in their promoter regions. Cis-elements related to resistance, physiological metabolism and hormonal regulation were identified within the promoter regions of the CKX genes. Expression analysis under different stress conditions (cold, heat, drought and salt) revealed different expression patterns in the different tissues. Through virus-induced gene silencing (VIGS), the GhCKX34A gene was found to improve cold resistance by modulating antioxidant-related activity. Since GhCKX29A is highly expressed during fibre development, we hypothesize that the increased expression of GhCKX29A in fibres has significant effects on fibre elongation. Consequently, these results contribute to our understanding of the involvement of GhCKXs in both fibre development and response to abiotic stress.

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.

Genome-wide identification and investigation of monosaccharide transporter gene family based on their evolution and expression analysis under abiotic stress and hormone treatments in maize (Zea mays L.).

BACKGROUND: Monosaccharide transporter (MST) family, as a carrier for monosaccharide transport, plays an important role in carbon partitioning and widely involves in plant growth and development, stress response, and signaling transduction. However, little information on the MST family genes is reported in maize (Zea mays), especially in response to abiotic stresses. In this study, the genome-wide identification of MST family genes was performed in maize. RESULT: A total of sixty-six putative members of MST gene family were identified and divided into seven subfamilies (including SPT, PMT, VGT, INT, pGlcT, TMT, and ERD) using bioinformatics approaches, and gene information, phylogenetic tree, chromosomal location, gene structure, motif composition, and cis-acting elements were investigated. Eight tandem and twelve segmental duplication events were identified, which played an important role in the expansion of the ZmMST family. Synteny analysis revealed the evolutionary features of MST genes in three gramineous crop species. The expression analysis indicated that most of the PMT, VGT, and ERD subfamilies members responded to osmotic and cadmium stresses, and some of them were regulated by ABA signaling, while only a few members of other subfamilies responded to stresses. In addition, only five genes were induced by NaCl stress in MST family. CONCLUSION: These results serve to understand the evolutionary relationships of the ZmMST family genes and supply some insight into the processes of monosaccharide transport and carbon partitioning on the balance between plant growth and development and stress response in maize.

In silico prediction of polyketide biosynthetic gene clusters in the genomes of Hypericum-borne endophytic fungi.

BACKGROUND: The search for new bioactive natural compounds with anticancer activity is still of great importance. Even though their potential for diagnostics and treatment of cancer has already been proved, the availability is still limited. Hypericin, a naphthodianthrone isolated essentially from plant source Hypericum perforatum L. along with other related anthraquinones and bisanthraquinones belongs to this group of compounds. Although it has been proven that hypericin is synthesized by the polyketide pathway in plants, none of the candidate genes coding for key enzymes has been experimentally validated yet. Despite the rare occurrence of anthraquinones in plants, their presence in microorganisms, including endophytic fungi, is quite common. Unlike plants, several biosynthetic genes grouped into clusters (BGCs) in fungal endophytes have already been characterized. RESULTS: The aim of this work was to predict, identify and characterize the anthraquinone BGCs in de novo assembled and functionally annotated genomes of selected endophytic fungal isolates (Fusarium oxysporum, Plectosphaerella cucumerina, Scedosporium apiospermum, Diaporthe eres, Canariomyces subthermophilus) obtained from different tissues of Hypericum spp. The number of predicted type I polyketide synthase (PKS) BGCs in the studied genomes varied. The non-reducing type I PKS lacking thioesterase domain and adjacent discrete gene encoding protein with product release function were identified only in the genomes of C. subthermophilus and D. eres. A candidate bisanthraquinone BGC was predicted in C. subthermophilus genome and comprised genes coding the enzymes that catalyze formation of the basic anthraquinone skeleton (PKS, metallo-beta-lactamase, decarboxylase, anthrone oxygenase), putative dimerization enzyme (cytochrome P450 monooxygenase), other tailoring enzymes (oxidoreductase, dehydrogenase/reductase), and non-catalytic proteins (fungal transcription factor, transporter protein). CONCLUSIONS: The results provide an insight into genetic background of anthraquinone biosynthesis in Hypericum-borne endophytes. The predicted bisanthraquinone gene cluster represents a basis for functional validation of the candidate biosynthetic genes in a simple eukaryotic system as a prospective biotechnological alternative for production of hypericin and related bioactive anthraquinones.

Control of a gene transfer agent cluster in Caulobacter crescentus by transcriptional activation and anti-termination.

Gene Transfer Agents (GTAs) are phage-like particles that cannot self-multiply and be infectious. Caulobacter crescentus, a bacterium best known as a model organism to study bacterial cell biology and cell cycle regulation, has recently been demonstrated to produce bona fide GTA particles (CcGTA). Since C. crescentus ultimately die to release GTA particles, the production of GTA particles must be tightly regulated and integrated with the host physiology to prevent a collapse in cell population. Two direct activators of the CcGTA biosynthetic gene cluster, GafY and GafZ, have been identified, however, it is unknown how GafYZ controls transcription or how they coordinate gene expression of the CcGTA gene cluster with other accessory genes elsewhere on the genome for complete CcGTA production. Here, we show that the CcGTA gene cluster is transcriptionally co-activated by GafY, integration host factor (IHF), and by GafZ-mediated transcription anti-termination. We present evidence that GafZ is a transcription anti-terminator that likely forms an anti-termination complex with RNA polymerase, NusA, NusG, and NusE to bypass transcription terminators within the 14 kb CcGTA cluster. Overall, we reveal a two-tier regulation that coordinates the synthesis of GTA particles in C. crescentus.

Genome-wide identification of the MED25 BINDING RING-H2 PROTEIN gene family in foxtail millet (Setaria italica L.) and the role of SiMBR2 in resistance to abiotic stress in Arabidopsis.

MAIN CONCLUSION: The SiMBR genes in foxtail millet were identified and studied. Heterologous expression of SiMBR2 in Arabidopsis can improve plant tolerance to drought stress by decreasing the level of reactive oxygen species. Foxtail millet (Setaria italica L.), a C4 crop recognized for its exceptional resistance to drought stress, presents an opportunity to improve the genetic resilience of other crops by examining its unique stress response genes and understanding the underlying molecular mechanisms of drought tolerance. In our previous study, we identified several genes linked to drought stress by transcriptome analysis, including SiMBR2 (Seita.7G226600), a member of the MED25 BINDING RING-H2 PROTEIN (MBR) gene family, which is related to protein ubiquitination. Here, we have identified ten SiMBR genes in foxtail millet and conducted analyses of their structural characteristics, chromosomal locations, cis-acting regulatory elements within their promoters, and predicted transcription patterns specific to various tissues or developmental stages using bioinformatic approaches. Further investigation of the stress response of SiMBR2 revealed that its transcription is induced by treatments with salicylic acid and gibberellic acid, as well as by salt and osmotic stresses, while exposure to high or low temperatures led to a decrease in its transcription levels. Heterologous expression of SiMBR2 in Arabidopsis thaliana enhanced the plant's tolerance to water deficit by reducing the accumulation of reactive oxygen species under drought stress. In summary, this study provides support for exploring the molecular mechanisms associated with drought resistance of SiMBR genes in foxtail millet and contributing to genetic improvement and molecular breeding in other crops.

Biosynthetic gene clusters with biotechnological applications in novel Antarctic isolates from Actinomycetota.

Actinomycetota have been widely described as valuable sources for the acquisition of secondary metabolites. Most microbial metabolites are produced via metabolic pathways encoded by biosynthetic gene clusters (BGCs). Although many secondary metabolites are not essential for the survival of bacteria, they play an important role in their adaptation and interactions within microbial communities. This is how bacteria isolated from extreme environments such as Antarctica could facilitate the discovery of new BGCs with biotechnological potential. This study aimed to isolate rare Actinomycetota strains from Antarctic soil and sediment samples and identify their metabolic potential based on genome mining and exploration of biosynthetic gene clusters. To this end, the strains were sequenced using Illumina and Oxford Nanopore Technologies platforms. The assemblies were annotated and subjected to phylogenetic analysis. Finally, the BGCs present in each genome were identified using the antiSMASH tool, and the biosynthetic diversity of the Micrococcaceae family was evaluated. Taxonomic annotation revealed that seven strains were new and two were previously reported in the NCBI database. Additionally, BGCs encoding type III polyketide synthases (T3PKS), beta-lactones, siderophores, and non-ribosomal peptide synthetases (NRPS) have been identified, among others. In addition, the sequence similarity network showed a predominant type of BGCs in the family Micrococcaceae, and some genera were distinctly grouped. The BGCs identified in the isolated strains could be associated with applications such as antimicrobials, anticancer agents, and plant growth promoters, among others, positioning them as excellent candidates for future biotechnological applications and innovations. KEY POINTS: • Novel Antarctic rare Actinomycetota strains were isolated from soil and sediments • Genome-based taxonomic affiliation revealed seven potentially novel species • Genome mining showed metabolic potential for novel natural products.

Genome-wide identification, molecular evolution and gene expression of P450 gene family in Cyrtotrachelus buqueti.

BACKGROUND: Insect Cytochrome P450 monooxygenase (CYPs or P450s) plays an important role in detoxifying insecticides, causing insect populations to develop resistance. However, the molecular functions of P450 gene family in Cyrtotrachelus buqueti genome are still lacking. RESULTS: In this study, 71 CbuP450 genes have been identified. The amino acids length of CbuP450 proteins was between 183 aa ~ 1041 aa. They are proteins with transmembrane domains. The main component of their secondary structure is α-helix and random coils. Phylogenetic analysis showed that C. buqueti and Rhynchophorus ferrugineus were the most closely related. This gene family has 29 high-frequency codons, which tend to use A/T bases and A/T ending codons. Gene expression analysis showed that CbuP450_23 in the female adult may play an important role on high temperature resistance, and CbuP450_17 in the larval may play an important role on low temperature tolerance. CbuP450_10, CbuP450_17, CbuP450_23, CbuP450_10, CbuP450_16, CbuP450_20, CbuP450_23 and CbuP450_ 29 may be related to the regulation of bamboo fiber degradation genes in C. buqueti. Protein interaction analysis indicates that most CbuP450 proteins are mainly divided into three aspects: encoding the biosynthesis of ecdysteroids, participating in the decomposition of synthetic insecticides, metabolizing insect hormones, and participating in the detoxification of compounds. CONCLUSIONS: We systematically analyzed the gene and protein characteristics, gene expression, and protein interactions of CbuP450 gene family, revealing the key genes involved in the stress response of CbuP450 gene family in the resistance of C. buqueti to high or low temperature stress, and identified the key CbuP450 proteins involved in important life activity metabolism. These results provided a reference for further research on the function of P450 gene family in C. buqueti.

Comprehensive analysis of genomic variation, pan-genome and biosynthetic potential of Corynebacterium glutamicum strains.

Corynebacterium glutamicum is a non-pathogenic species of the Corynebacteriaceae family. It has been broadly used in industrial biotechnology for the production of valuable products. Though it is widely accepted at the industrial level, knowledge about the genomic diversity of the strains is limited. Here, we investigated the comparative genomic features of the strains and pan-genomic characteristics. We also observed phylogenetic relationships among the strains based on average nucleotide identity (ANI). We found diversity between strains at the genomic and pan-genomic levels. Less than one-third of the C. glutamicum pan-genome consists of core genes and soft-core genes. Whereas, a large number of strain-specific genes covered about half of the total pan-genome. Besides, C. glutamicum pan-genome is open and expanding, which indicates the possible addition of new gene families to the pan-genome. We also investigated the distribution of biosynthetic gene clusters (BGCs) among the strains. We discovered slight variations of BGCs at the strain level. Several BGCs with the potential to express novel bioactive secondary metabolites have been identified. Therefore, by utilizing the characteristic advantages of C. glutamicum, different strains can be potential applicants for natural drug discovery.

Identification and analysis of MAPK cascade gene families of Camellia oleifera and their roles in response to cold stress.

BACKGROUND: Low-temperature severely limits the growth and development of Camellia oleifera (C. oleifera). The mitogen-activated protein kinase (MAPK) cascade plays a key role in the response to cold stress. METHODS AND RESULTS: Our study aims to identify MAPK cascade genes in C. oleifera and reveal their roles in response to cold stress. In our study, we systematically identified and analyzed the MAPK cascade gene families of C. oleifera, including their physical and chemical properties, conserved motifs, and multiple sequence alignments. In addition, we characterized the interacting networks of MAPKK kinase (MAPKKK)-MAPK kinase (MAPKK)-MAPK in C. oleifera. The molecular mechanism of cold stress resistance of MAPK cascade genes in wild C. oleifera was analyzed by differential gene expression and real-time quantitative reverse transcription-PCR (qRT-PCR). CONCLUSION: In this study, 21 MAPKs, 4 MAPKKs and 55 MAPKKKs genes were identified in the leaf transcriptome of C. oleifera. According to the phylogenetic results, MAPKs were divided into 4 groups (A, B, C and D), MAPKKs were divided into 3 groups (A, B and D), and MAPKKKs were divided into 2 groups (MEKK and Raf). Motif analysis showed that the motifs in each subfamily were conserved, and most of the motifs in the same subfamily were basically the same. The protein interaction network based on Arabidopsis thaliana (A. thaliana) homologs revealed that MAPK, MAPKK, and MAPKKK genes were widely involved in C. oleifera growth and development and in responses to biotic and abiotic stresses. Gene expression analysis revealed that the CoMAPKKK5/CoMAPKKK43/CoMAPKKK49-CoMAPKK4-CoMAPK8 module may play a key role in the cold stress resistance of wild C. oleifera at a high-elevation site in Lu Mountain (LSG). This study can facilitate the mining and utilization of genetic resources of C. oleifera with low-temperature tolerance.

Aspergillus flavus pangenome (AflaPan) uncovers novel aflatoxin and secondary metabolite associated gene clusters.

BACKGROUND: Aspergillus flavus is an important agricultural and food safety threat due to its production of carcinogenic aflatoxins. It has high level of genetic diversity that is adapted to various environments. Recently, we reported two reference genomes of A. flavus isolates, AF13 (MAT1-2 and highly aflatoxigenic isolate) and NRRL3357 (MAT1-1 and moderate aflatoxin producer). Where, an insertion of 310 kb in AF13 included an aflatoxin producing gene bZIP transcription factor, named atfC. Observations of significant genomic variants between these isolates of contrasting phenotypes prompted an investigation into variation among other agricultural isolates of A. flavus with the goal of discovering novel genes potentially associated with aflatoxin production regulation. Present study was designed with three main objectives: (1) collection of large number of A. flavus isolates from diverse sources including maize plants and field soils; (2) whole genome sequencing of collected isolates and development of a pangenome; and (3) pangenome-wide association study (Pan-GWAS) to identify novel secondary metabolite cluster genes. RESULTS: Pangenome analysis of 346 A. flavus isolates identified a total of 17,855 unique orthologous gene clusters, with mere 41% (7,315) core genes and 59% (10,540) accessory genes indicating accumulation of high genomic diversity during domestication. 5,994 orthologous gene clusters in accessory genome not annotated in either the A. flavus AF13 or NRRL3357 reference genomes. Pan-genome wide association analysis of the genomic variations identified 391 significant associated pan-genes associated with aflatoxin production. Interestingly, most of the significantly associated pan-genes (94%; 369 associations) belonged to accessory genome indicating that genome expansion has resulted in the incorporation of new genes associated with aflatoxin and other secondary metabolites. CONCLUSION: In summary, this study provides complete pangenome framework for the species of Aspergillus flavus along with associated genes for pathogen survival and aflatoxin production. The large accessory genome indicated large genome diversity in the species A. flavus, however AflaPan is a closed pangenome represents optimum diversity of species A. flavus. Most importantly, the newly identified aflatoxin producing gene clusters will be a new source for seeking aflatoxin mitigation strategies and needs new attention in research.