Non-redundant Protein Database Research Articles

Evolution of sequence, structural and functional diversity of the ubiquitous DNA/RNA-binding Alba domain

The DNA/RNA-binding Alba domain is prevalent across all kingdoms of life. First discovered in archaea, this protein domain has evolved from RNA- to DNA-binding, with a concomitant expansion in the range of cellular processes that it regulates. Despite its widespread presence, the full extent of its sequence, structural, and functional diversity remains unexplored. In this study, we employed iterative searches in PSI-BLAST to identify 15,161 unique Alba domain-containing proteins from the NCBI non-redundant protein database. Sequence similarity network (SSN) analysis clustered them into 13 distinct subgroups, including the archaeal Alba and eukaryotic Rpp20/Pop7 and Rpp25/Pop6 groups, as well as novel fungal and Plasmodium-specific Albas. Sequence and structural conservation analysis of the subgroups indicated high preservation of the dimer interface, with Alba domains from unicellular eukaryotes notably exhibiting structural deviations towards their C-terminal end. Finally, phylogenetic analysis, while supporting SSN clustering, revealed the evolutionary branchpoint at which the eukaryotic Rpp20- and Rpp25-like clades emerged from archaeal Albas, and the subsequent taxonomic lineage-based divergence within each clade. Taken together, this comprehensive analysis enhances our understanding of the evolutionary history of Alba domain-containing proteins across diverse organisms.

Comparative Genomic Analysis of an Apiotrichum cacaoliposimilis Strain Isolated from a Patient with Urinary Tract Infection.

Opportunistic infections caused by fungi, particularly those occurring in immunocompromised patients, are considered challenging worldwide. Therefore, a comprehensive understanding of pathogenic fungi is necessary. The present study reports the isolation of a strain of Apiotrichum cacaoliposimilis, which is difficult to detect using conventional clinical assays, from the sterile urine samples of a patient with a urinary tract infection. Sanger sequencing of the internal transcribed spacer regions confirmed the genus of the microbe, while whole-genome sequencing yielded the initial genome assembly of A. cacaoliposimilis. A total of 7,161 predicted proteincoding genes were mapped using multiple databases, including Gene Ontology, Kyoto Encyclopedia of Genes and Genomes, non-redundant protein database, Pathogen-Host Interactions Database, and Comprehensive Antibiotic Resistance Database. The phenotypic data, biochemical reactions, and antimicrobial susceptibility analyses were conducted to reveal the metabolic properties, virulence, and drug resistance profile of the isolated A. cacaoliposimilis. The rank-sum test revealed the differences in the intergeneric distribution of the highly virulent genes UgeB and Pem2. In addition, other genes exhibited significant overlap in terms of virulence factors with the clinical isolate Apiotrichum mycotoxinivorans GMU1709. Fortunately, similar to most fungi belonging to the Apiotrichum genus, the isolate investigated in the present study was also sensitive to the drug voriconazole (MIC = 0.06 μg/ml). In summary, the phylogenetic placement, potential pathogenic genes, drug sensitivity patterns, and morphological characteristics of the isolated A. cacaoliposimilis were determined precisely in the present study.

SoIR: a comprehensive Solanaceae information resource for comparative and functional genomic study.

The Solanaceae family, which includes economically important crops such as tomatoes, potatoes and peppers, has experienced a rapid expansion in genomic data due to advancements in sequencing technologies. However, existing databases are limited by incomplete species representation, a lack of comprehensive comparative genomic tools and the absence of systematic pan-genomic analyses. To address these gaps, we developed the Solanaceae Information Resource (SoIR, https://soir.bio2db.com), a comprehensive genomics database for the Solanaceae family. SoIR integrates genomic data from 81 species and transcriptomic data from 41 species, encompassing a total of 3908408 gene annotations derived from Gene Ontology, nonredundant protein, Pfam, Swiss-Prot and TrEMBL databases. The resource also includes 3437115 CRISPR guide sequences, 212395 transcription factorsand 19086 genes associated with methylation modification. In addition to species-specific analyses, SoIR provides extensive bioinformatics tools for investigating gene family evolution, phylogenetic relationships and karyotype reconstruction across 25 fully sequenced genomes. With advanced tools such as Blast, Synteny and Sequence Alignment, the platform provides users with interactive and intuitive visualizations for conducting cross-species comparative genomics. As the first comprehensive pan-genomic resource for the entire Solanaceae family, SoIR facilitates in-depth cross-species analysis, supporting global research initiatives in plant evolution, functional genomics and crop improvement.

Transcriptome sequencing of medical herb Salvia Rosmarinus (Rosemary) revealed the phenylpropanoid biosynthesis pathway genes and their phylogenetic relationships.

The Salvia rosmarinus spenn. (rosemary) is considered an economically important ornamental and medicinal plant and is widely utilized in culinary and for treating several diseases. However, the procedure behind synthesizing secondary metabolites-based bioactive compounds at the molecular level in S. rosmarinus is not explored completely. We performed transcriptomic sequencing of the pooled sample from leaf and stem tissues on the Illumina HiSeqTM X10 platform. The transcriptomics analysis led to the generation of 29,523,608 raw reads, followed by data pre-processing which generated 23,208,592 clean reads, and de novo assembly of S. rosmarinus obtained 166,849 unigenes. Among them, nearly 75.1% of unigenes i.e., 28,757 were interpreted against a non-redundant protein database. The gene ontology-based annotation classified them into 3 main categories and 55 sub-categories, and clusters of orthologous genes annotation categorized them into 23 functional categories. The Kyoto Encyclopedia of Genes and Genomes database-based pathway analysis confirmed the involvement of 13,402 unigenes in 183 biochemical pathways, among these unigenes, 1,186 are involved in the 17 secondary metabolite production pathways. Several key enzymes involved in producing aromatic amino acids and phenylpropanoids were identified from the transcriptome database. Among the identified 48 families of transcription factors from coding unigenes, bHLH, MYB, WRKYs, NAC, C2H2, C3H, and ERF are involved in flavonoids and other secondary metabolites biosynthesis. The phylogenetic analysis revealed the evolutionary relationship between the phenylpropanoid pathway genes of rosemary with other members of Lamiaceae. Our work reveals a new molecular mechanism behind the biosynthesis of phenylpropanoids and their regulation in rosemary plants.

Transcriptomic analysis of nitrogen metabolism pathways in Klebsiella aerogenes under nitrogen-rich conditions.

The acceleration of the nitrogen cycle and the nitrogen excess observed in some coastal waters has increased interest into understanding the biochemical and molecular basis of nitrogen metabolism in various microorganisms. To investigate nitrogen metabolism of a novel heterotrophic nitrification and aerobic denitrification bacterium Klebsiella aerogenes strain (B23) under nitrogen-rich conditions, we conducted physiological and transcriptomic high-throughput sequencing analyses on strain B23 cultured on potassium nitrate-free or potassium nitrate-rich media. Overall, K. aerogenes B23 assimilated 82.47% of the nitrate present into cellular nitrogen. Further, 1,195 differentially expressed genes were observed between K. aerogenes B23 cultured on potassium nitrate-free media and those cultured on potassium nitrate-rich media. Gene annotation and metabolic pathway analysis of the transcriptome were performed using a series of bioinformatics tools, including Gene Ontology, Kyoto Encyclopedia of Genes and Genomes, and Non-Redundant Protein Database annotation. Accordingly, the nitrogen metabolism pathway of K. aerogenes B23 was analyzed; overall, 39 genes were determined to be involved in this pathway. Differential expression analysis of the genes involved in the nitrogen metabolism pathway demonstrated that, compared to the control, FNR, NarK/14945, fdx, gshA, proB, proA, gapA, argH, artQ, artJ, artM, ArgR, GAT1, prmB, pyrG, glnS, and Ca1 were significantly upregulated in the nitrogen-treated K. aerogenes B23; these genes have been established to be involved in the regulation of nitrate, arginine, glutamate, and ammonia assimilation. Further, norV, norR, and narI were also upregulated in nitrogen-treated K. aerogenes B23; these genes are involved in the regulation of NO metabolism. These differential expression results are important for understanding the regulation process of key nitrogen metabolism enzyme genes in K. aerogenes B23. Therefore, this study establishes a solid foundation for further research into the expression regulation patterns of nitrogen metabolism-associated genes in K. aerogenes B23 under nitrogen-rich conditions; moreover, this research provides essential insight into how K. aerogenes B23 utilizes nutritional elements.

Exploring a GtfB-Type 4,6-α-Glucanotransferase to Synthesize the (α1 → 6) Linkages in Linear Chain and Branching Points from Amylose and Enhance the Functional Property of Granular Corn Starches.

Starch-converting α-glucanotransferases of glycoside hydrolase family 70 (GH70) are promising enzymatic tools for the production of diverse α-glucans with (potential) commercial applications in food and health and as biomaterials. In this study, a novel GtfB enzyme from Weissella confusa MBF8-1 was screened in the National Center for Biotechnology Information (NCBI) nonredundant protein database. The enzyme (named WcMBF8-1 GtfB) displayed high conservation in motifs I-IV with other GtfB enzymes but possessed unique variations in several substrate-binding residues. Structural characterizations of its α-glucan products revealed that WcMBF8-1 GtfB exhibited an atypical 4,6-α-glucanotransferase activity and was capable of catalyzing, by cleaving off (α1 → 4)-linkages in starch-like substrates and the synthesis of linear (α1 → 6) linkages and (α1 → 4,6) branching points. The product specificity enlarges the diversity of α-glucans and facilitates recognition of the determinants of the linkage specificity in GtfB enzymes. Furthermore, the contents of slowly digestible starch and resistant starch of granular corn starches, modified by WcMBF8-1 GtfB, increased by 6.7%, which suggested the potential value for the utilization of WcMBF8-1 GtfB to prepare "clean-label" starch ingredients with improved functional attributes.

Whole-genome sequencing of a protease-producing strain isolated from fermented golden pompano and molecular docking of the protease to fish proteins

Whole-genome sequencing of a protease-producing strain isolated from fermented golden pompano and molecular docking of the protease to fish proteins

Transcriptome analysis and identification of differentially expressed genes between early and mature ovarian stages in the female mantis shrimp (Harpiosquilla raphidea) using RNA-Seq

Transcriptome analysis and identification of differentially expressed genes between early and mature ovarian stages in the female mantis shrimp (Harpiosquilla raphidea) using RNA-Seq

Microbial Diversity in Children with Gastroenteritis in the Amazon Region of Brazil: Development and Validation of a Molecular Method for Complete Sequencing of Viral Genomes.

Metagenomic sequencing is a powerful tool that is widely used in laboratories worldwide for taxonomic characterization of microorganisms in clinical and environmental samples. In this study, we utilized metagenomics to investigate comprehensively the microbial diversity in fecal samples of children over a four-year period. Our methods were carefully designed to ensure accurate and reliable results. Validated and analyzed were metagenomic data obtained from sequencing 27 fecal samples from children under 10 years old with gastroenteritis over a four-year period (2012-2016). The fecal specimens were collected from patients who received care at public health facilities in the northern region of Brazil. Sequencing libraries were prepared from cDNA and sequenced on the Illumina HiSeq. Kraken-2 was utilized to classify bacterial taxonomy based on the 16S rRNA gene, using the Silva rRNA database. Additionally, the Diamond program was used for mapping to the non-redundant protein database (NR database). Phylogenomic analyses were conducted using Geneious R10 and MEGA X software, and Bayesian estimation of phylogeny was performed using the MrBayes program. The results indicate significant heterogeneity among norovirus strains, with evidence of recombination and point mutations. This study presents the first complete genome of parechovirus 8 in the region. Additionally, it describes the bacterial populations and bacteriophages present in feces, with a high abundance of Firmicutes and Proteobacteria, including an increased proportion of the Enterobacteriaceae family. The presented data demonstrate the genetic diversity of microbial populations and provide a comprehensive report on viral molecular characterization. These findings are relevant for genomic studies in gastrointestinal infections. The metagenomic approach is a powerful tool for investigating microbial diversity in children with gastroenteritis. However, further studies are imperative to conduct genomic analysis of identified bacterial strains and thoroughly analyze antimicrobial resistance genes.

AI-assisted structural consensus-proteome prediction of human monkeypox viruses isolated within a year after the 2022 multi-country outbreak.

The 2022 outbreak of the monkeypox virus already involves, by April 2023, 110 countries with 86,956 confirmed cases and 119 deaths. Understanding an emerging disease on a molecular level is essential to study infection processes and eventually guide drug discovery at an early stage. To support this, we provide the so far most comprehensive structural proteome of the monkeypox virus, which includes 210 structural models, each computed with three state-of-the-art structure prediction methods. Instead of building on a single-genome sequence, we generated our models from a consensus of 3,713 high-quality genome sequences sampled from patients within 1 year of the outbreak. Therefore, we present an average structural proteome of the currently isolated viruses, including mutational analyses with a special focus on drug-binding sites. Continuing dynamic mutation monitoring within the structural proteome presented here is essential to timely predict possible physiological changes in the evolving virus.

The Draft Genome of the "Golden Tide" Seaweed, Sargassum horneri: Characterization and Comparative Analysis.

Sargassum horneri, a prevalent species of brown algae found along the coast of the northwest Pacific Ocean, holds significant importance as a valuable source of bioactive compounds. However, its rapid growth can lead to the formation of a destructive "golden tide", causing severe damage to the local economy and coastal ecosystems. In this study, we carried out de novo whole-genome sequencing of S. horneri using next-generation sequencing to unravel the genetic information of this alga. By utilizing a reference-guided de novo assembly pipeline with a closely related species, we successfully established a final assembled genome with a total length of 385 Mb. Repetitive sequences made up approximately 30.6% of this genome. Among the identified putative genes, around 87.03% showed homology with entries in the NCBI non-redundant protein database, with Ectocarpus siliculosus being the most closely related species for approximately one-third of these genes. One gene encoding an alkaline phosphatase family protein was found to exhibit positive selection, which could give a clue for the formation of S. horneri golden tides. Additionally, we characterized putative genes involved in fucoidan biosynthesis metabolism, a significant pathway in S. horneri. This study represents the first genome-wide characterization of a S. horneri species, providing crucial insights for future investigations, such as ecological genomic analyses.

NCRD: A non-redundant comprehensive database for detecting antibiotic resistance genes

Antibiotic resistance genes (ARGs) are emerging pollutants present in various environments. Identifying ARGs has become a growing concern in recent years. Several databases, including the Antibiotic Resistance Genes Database (ARDB), Comprehensive Antibiotic Resistance Database (CARD), and Structured Antibiotic Resistance Genes (SARG), have been applied to detect ARGs. However, these databases have limitations, which hinder the comprehensive profiling of ARGs in environmental samples. To address these issues, we constructed a non-redundant antibiotic resistance genes database (NRD) by consolidating sequences from ARDB, CARD, and SARG. We identified the homologous proteins of NRD from Non-redundant Protein Database (NR) and the Protein DataBank Database (PDB) and clustered them to establish a non-redundant comprehensive antibiotic resistance genes database (NCRD) with similarities of 100% (NCRD100) and 95% (NCRD95). To demonstrate the advantages of NCRD, we compared it with other databases by using metagenome datasets. Results revealed its strong ability in detecting potential ARGs.

First Report of Banana mild mosaic virus Infecting Banana in China.

Banana (Musa spp.) is an important fruit in tropical and subtropical regions and an essential food crop in some developing countries. China has a long history of banana cultivation and ranks second in global banana production, with a planting area exceeding 11 million hectares (FAOSTAT, 2023). Banana mild mosaic virus (BanMMV) is a flexuous filamentous virus infecting bananas and a banmivirus in the Betaflexiviridae family. Its infection often results in symptomless plants of Musa spp., and the virus has a worldwide distribution, which can explain its high prevalence (Kumar et al., 2015). BanMMV infection often causes transitory symptoms, such as mild chlorotic streaks and mosaics, on young leaves (Thomas, 2015). The mixed infection of BanMMV with other banana-infecting viruses such as banana streak viruses (BSV) and cucumber mosaic virus (CMV), can exacerbate the mosaic symptoms of BanMMV (Fidan et al., 2019). In October 2021, we collected twenty-six leaf samples of suspected viral disease of bananas from four cities (Huizhou, Qingyuan, Zhanjiang, and Yangjiang) in Guangdong province, two cities (Hekou and Jinghong) in Yunnan province, two cities (Yulin and Wuming) in Guangxi Zhuang autonomous region. After fully mixing these infected samples, we divided them into two pools and sent them to Shanghai Biotechnology Corporation (China) for metatranscriptome sequencing. Each sample contained about 5 g of leaves in total. Zymo-Seq RiboFree Total RNA Library Prep Kit (Zymo Research, USA) was used for ribosomal RNA depletion and library preparations. Illumina sequencing (Illumina NovaSeq 6000) was carried out by Shanghai Biotechnology Corporation (China). Paired-end (150 bp) sequencing of the RNA library was performed on an Illumina HiSeq 2000/2500 platform. Clean reads were assembled by a metagenomic de novo assembly using the CLC Genomics Workbench (version: 6.0.4). Then the non-redundant protein database in the National Center for Biotechnology Information (NCBI) was used for BLASTx annotation. A total of 79,528 contigs were generated from the clean reads (68,878,162) through de novo assembly. A contig of 7265 nucleotides (nts) showed the highest nucleotide sequence identity (90.08%) to the genome of BanMMV isolate EM4-2 (GenBank accession no. OL826745.1). We designed specific primers according to the BanMMV CP gene (Table S1), tested the twenty-six leaf samples collected from the above-mentioned eight cities, and found that only one sample of Fenjiao (Musa ABB Pisang Awak) in Guangzhou city was infected with this virus. The symptoms of banana leaves containing BanMMV were slight chlorosis and yellowing of leaf edges (Fig. S1). We failed to detect other banana viruses, such as BSV, CMV, and banana bunchy top virus (BBTV) in the BanMMV-infected banana leaves. RNA from the infected leaves was extracted, and the assembled contig was confirmed by overlapping PCR amplification across the whole sequence (Table S1). All ambiguous regions were amplified by PCR and RACE, and the products were subjected to Sanger sequencing. The complete genome of the virus candidate was 7310 nts in length, excluding the poly (A) tail. The sequence was deposited in GenBank under accession number ON227268 (isolate BanMMV-GZ from Guangzhou). A schematic representation of the genome organization of BanMMV-GZ is shown in Fig. S2. Its genome has five open reading frames (ORF) encoding RNA-dependent RNA polymerase (RdRp), three triple gene block proteins necessary for cell-to-cell movement (TGBp1 to TGBp3) and a coat protein (CP), similar to other BanMMV isolates (Kondo et al., 2021). Phylogenetic analyses of the complete nt sequence of the full genome and RdRp gene using the neighbor-joining (NJ) method also clearly placed the BanMMV-GZ firmly within all isolates of BanMMV (Fig. S3). To our knowledge, this is the first report of BanMMV infecting bananas in China, extending the geographical range of this viral disease around the world. Accordingly, larger-scale BanMMV investigations must be conducted to determine the distribution and prevalence of BanMMV in China.

More than mcr: canonical plasmid- and transposon-encoded mobilized colistin resistance genes represent a subset of phosphoethanolamine transferases.

Mobilized colistin resistance genes (mcr) may confer resistance to the last-resort antimicrobial colistin and can often be transmitted horizontally. mcr encode phosphoethanolamine transferases (PET), which are closely related to chromosomally encoded, intrinsic lipid modification PET (i-PET; e.g., EptA, EptB, CptA). To gain insight into the evolution of mcr within the context of i-PET, we identified 69,814 MCR-like proteins present across 256 bacterial genera (obtained by querying known MCR family representatives against the National Center for Biotechnology Information [NCBI] non-redundant protein database via protein BLAST). We subsequently identified 125 putative novel mcr-like genes, which were located on the same contig as (i) ≥1 plasmid replicon and (ii) ≥1 additional antimicrobial resistance gene (obtained by querying the PlasmidFinder database and NCBI's National Database of Antibiotic Resistant Organisms, respectively, via nucleotide BLAST). At 80% amino acid identity, these putative novel MCR-like proteins formed 13 clusters, five of which represented putative novel MCR families. Sequence similarity and a maximum likelihood phylogeny of mcr, putative novel mcr-like, and ipet genes indicated that sequence similarity was insufficient to discriminate mcr from ipet genes. A mixed-effect model of evolution (MEME) indicated that site- and branch-specific positive selection played a role in the evolution of alleles within the mcr-2 and mcr-9 families. MEME suggested that positive selection played a role in the diversification of several residues in structurally important regions, including (i) a bridging region that connects the membrane-bound and catalytic periplasmic domains, and (ii) a periplasmic loop juxtaposing the substrate entry tunnel. Moreover, eptA and mcr were localized within different genomic contexts. Canonical eptA genes were typically chromosomally encoded in an operon with a two-component regulatory system or adjacent to a TetR-type regulator. Conversely, mcr were represented by single-gene operons or adjacent to pap2 and dgkA, which encode a PAP2 family lipid A phosphatase and diacylglycerol kinase, respectively. Our data suggest that eptA can give rise to "colistin resistance genes" through various mechanisms, including mobilization, selection, and diversification of genomic context and regulatory pathways. These mechanisms likely altered gene expression levels and enzyme activity, allowing bona fide eptA to evolve to function in colistin resistance.

Genome sequencing, assembly, and characterization of Pichia fermentans Z9Y-3 as a non-Saccharomyces yeast with aroma enhancing potential

Genome sequencing, assembly, and characterization of Pichia fermentans Z9Y-3 as a non-Saccharomyces yeast with aroma enhancing potential

Phosphorylation of SiRAV1 at Ser31 regulates the SiCAT expression to enhance salt tolerance in Setaria italica

Phosphorylation of SiRAV1 at Ser31 regulates the SiCAT expression to enhance salt tolerance in Setaria italica

Genome Sequence Resource for Colletotrichum gloeosporioides, an Important Pathogenic Fungus Causing Anthracnose of Dioscorea alata.

Anthracose disease is one of the most important diseases of Dioscorea alata and many other food yams, which is caused by Colletotrichum gloeosporioides fugus from the Glomerellaceae family of the Sordariomycetes class. In the present study, a C. gloeosporioides starin named CgDa01 was isolated from D. alata, and its genome was sequenced based on Oxford Nanopore technology (ONT) and the Illumina sequencing platform. The high-quality genome of CgDa01 was assembled with a 62.78 Mb genome size and 15,845 predicted protein-coding genes. The proteins of predicted genes were annotated using multiple public databases, including the nonredundant protein database, the InterProScan databases, and Kyoto Encyclopedia of Genes and Genomes. Among the annotated protein-coding genes, 55 were predicted as potential virulence genes by the fungal virulence factor database. The C. gloeosporioides CgDa01 genome assembly described in this study can serve as a resource for better understanding the pathogenic mechanism of C. gloeosporioides on yam hosts.

Genome Sequence Resource of Ustilago crameri, a Fungal Pathogen Causing Millet Smut Disease of Foxtail Millet.

Genome Sequence Resource of Ustilago crameri, a Fungal Pathogen Causing Millet Smut Disease of Foxtail Millet.

Transcriptome Analysis of a Wild Eggplant Germplasm M239 in Response to Verticillium dahliae Infection

In this study, wild eggplant germplasm No. M239, which is highly susceptible to Verticillium wilt, was used as the experimental material. The physiological and biochemical indices (SOD, PAL, MDA and soluble protein) of M239 roots were measured at different times (0, 12, 24, 36, 48, 60 and 72 h) post inoculation with Verticillium dahliae, and the key time points for the M239 response to Verticillium wilt infection were screened. Then, RNA-Seq technology was used to screen the differentially expressed genes (DEGs) in M239 roots at 0, 12 and 48 h post-inoculation (hpi). The transcriptional results of M239 were also compared with those resistance genes from some reported wild relative Solanum species (S. sisymbriifolium and S. aculeatissimum). Then some DEGs were chosen for validation by qRT–PCR. The results showed that 12 and 48 hpi were the turning points in the changes in all physiological and biochemical indices. A total of 6,783 DEGs were identified by RNA-Seq, including 6,141 DEGs (3,046 upregulated and 3,095 downregulated) at the M_12 h vs. M_0 h, 1,903 DEGs (792 upregulated and 1,111 downregulated) at M_48 h vs. M_12 h, and 1,261 DEGs that appeared simultaneously in both stages. KEGG enrichment analysis showed that there were 5 metabolic pathways enriched from DEGs, which were mostly related to primary metabolism, such as glycolysis, amino acid and ribosome biogenesis. Compared with the NCBI non-redundant protein (NR) database, one Ve2 homologous gene and 8 PR protein-related genes were screened. Transcription factor analysis showed that there were a large number of DEGs, such as MYB, AP2-EREBP, bHLH, NAC and Orphans in the two stages. Compared with the reported Verticillium wilt-resistant wild eggplant species, it was found that there were fewer genes and enriched metabolic pathways in the M239 response to Verticillium wilt infection, and it also lacked the response of some known key resistance genes. These results proved that the above resistance genes and metabolic pathways played a key role in the wild eggplant response to V. dahliae infection.

Transcriptomic Identification and Expression Profile Analysis of Odorant-Degrading Enzymes from the Asian Corn Borer Moth, Ostrinia furnacalis.

The Asian corn borer moth Ostrinia furnacalis is an important lepidopteran pest of maize in Asia. Odorant-degrading enzymes (ODEs), including carboxylesterases (CCEs), glutathione S-transferases (GSTs), cytochrome P450s (CYPs), UDP-glycosyltransferases (UGTs), and aldehyde oxidases (AOXs), are responsible for rapid inactivation of odorant signals in the insect antennae. In this study, we performed a transcriptome assembly for the antennae of O. furnacalis to identify putative ODE genes. Transcriptome sequencing revealed 35,056 unigenes, and 21,012 (59.94%) of these were annotated by searching against the reference sequences in the NCBI non-redundant (NR) protein database. For functional classification, these unigenes were subjected to Gene Ontology (GO), Eukaryotic Orthologous Groups (KOG), and Kyoto Encyclopedia of Genes and Genomes (KEGG) annotations. We identified 79 genes encoding putative ODEs: 19 CCEs, 17 GSTs, 24 CYPs, 13 UGTs, and 6 AOXs. BLASTX best hit results indicated that these genes shared quite high amino acid identities with their respective orthologs from other lepidopteran species. Reverse transcription-quantitative PCR showed that OfurCCE2, OfurCCE5, and OfurCCE18 were enriched in male antennae, while OfurCCE7 and OfurCCE10 were enriched in female antennae. OfurCCE14 and OfurCCE15 were expressed at near-equal amounts in the antennae of both sexes. Our findings establish a solid foundation for future studies aimed at understanding the olfactory functions of these genes in O. furnacalis.