AlliumDB: a central portal for comparative and functional genomics in Allium

The genus Allium belongs to the botanical family Amaryllidaceae and includes economically important crops such as onion, garlic, bunching onion, and leek, used as vegetables, spices, and traditional medicines. The large sizes of Allium genomes hamper the genetic dissection of agronomically important traits and molecular breeding. With the growing accumulation of genomic, resequencing, transcriptome, and phenotypic data, the demand for an integrative Allium database is increasing. Here we present a user-friendly database, AlliumDB (https://allium.qau.edu.cn), as a functional genomics hub integrating public and in-house data. The database contains all currently available nuclear and organelle genomes for Allium species, with genes comprehensively annotated based on Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses, orthology, gene families, protein families (Pfam), and non-coding RNA families (Rfam). Transcriptome and variation profiles are integrated into dynamic visualization tools. We took phenotypic photographs and generated trait records for hundreds of Allium germplasms collected worldwide, which are included in the database. We incorporated JBrowse for the visualization of gene structures, RNA sequencing data, and variation data. Analysis tools such as the basic local alignment search tool (BLAST), sequence fetch, enrichment, and motif analyses are available to explore potential gene functions. This database incorporates comprehensive Allium genotypic and phenotypic datasets. As the community assembles new genomes and generates resequencing data for Allium germplasms, the database will be improved and continuously updated with these multi-omics data and comparative genomic studies. We expect the AlliumDB database to become a key resource for the study of Allium crops.

ASAP: a platform for gene functional analysis in Angelica sinensis

BACKGROUND

Angelica sinensis (Danggui), a renowned medicinal orchid, has gained significant recognition for its therapeutic effects in treating a wide range of ailments. Genome information serves as a valuable resource, enabling researchers to gain a deeper understanding of gene function. In recent times, the availability of chromosome-level genomes for A. sinensis has opened up vast opportunities for exploring gene functionality. Integrating multiomics data can allow researchers to unravel the intricate mechanisms underlying gene function in A. sinensis and further enhance our knowledge of its medicinal properties.

RESULTS

In this study, we utilized genomic and transcriptomic data to construct a coexpression network for A. sinensis. To annotate genes, we aligned them with sequences from various databases, such as the NR, TAIR, trEMBL, UniProt, and SwissProt databases. For GO and KEGG annotations, we employed InterProScan and GhostKOALA software. Additionally, gene families were predicted using iTAK, HMMER, OrholoFinder, and KEGG annotation. To facilitate gene functional analysis in A. sinensis, we developed a comprehensive platform that integrates genomic and transcriptomic data with processed functional annotations. The platform includes several tools, such as BLAST, GSEA, Heatmap, JBrowse, and Sequence Extraction. This integrated resource and approach will enable researchers to explore the functional aspects of genes in A. sinensis more effectively.

CONCLUSION

We developed a platform, named ASAP, to facilitate gene functional analysis in A. sinensis. ASAP ( www.gzybioinformatics.cn/ASAP ) offers a comprehensive collection of genome data, transcriptome resources, and analysis tools. This platform serves as a valuable resource for researchers conducting gene functional research in their projects, providing them with the necessary data and tools to enhance their studies.

High-throughput sequencing yields a complete mitochondrial genome of Emberiza godlewskii (aves, emberidae)

Emberiza godlewskii (Taczanowski, 1874) is a passerine bird of eastern Asia which belongs to the genus Emberiza in the bunting family Emberizidae. Obtaining the complete mitochondrial genome sequence of E. godlewskii is helpful to understand the species delimitation for further study of E. cia/godlewskii complex. The circular genome (16,839 bp in length) contains 13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes, and 1 control region. The base composition shows that purine percentage (52.87%) is slightly higher than pyrimidine (47.13%). The phylogenetic analysis based on the published data of the mitochondrial genome showed that E. godlewskii is closely related to E. cioides. This new mitochondrial genome data will provide more essential molecular data for further study related to E. cia/godlewskii complex.

RoseAP: an analytical platform for gene function of Rosa rugosa

Rosa rugosa, a perennial shrub belonging to family Rosaceae, is a well-known ornamental plant. Its petals contain an abundance of essential oils and anthocyanins with enormous economic and health benefits when used as edible or cosmetic ingredients. The whole genome of R. rugosa was sequenced in 2021, which provided opportunities and challenges for gene regulation. However, many gene functions remain unknown. Therefore, an analytical platform named RoseAP (http://www.gzybioinformatics.cn/RoseAP/index.php) for the functional analysis of R. rugosa genes was constructed. It improved the gene annotation rate by integrating and analyzing genomic and transcriptomic datasets. First, 38,815 genes, covering 97.76% of the coding genes, were annotated functionally and structurally using a variety of algorithms and rules. Second, a total of 33 transcriptome samples were integrated, including 23 samples from our lab and 10 samples from the SRA database. A co-expression network containing approximately 29,657 positive or negative gene pairs, covering 74.7% of the coding genes, was constructed based on PCC and MR algorithms. Network analysis revealed that the DFR function was closely related to anthocyanin metabolism. It demonstrated the reliability of the network. Several SAUR genes of R. rugosa shared similar expression patterns. RoseAP was used to determine the sequence, structure, functional annotation, expression profile, regulatory network, and functional modules at the transcriptional and protein levels by inputting gene IDs. In addition, auxiliary analytical tools, including BLAST, gene set enrichment, orthologue conversion, gene sequence extraction, gene expression value extraction, and JBrowse, were utilized. Regular updates to RoseAP are expected to facilitate mining of gene function and promote genetic improvement in R. rugosa.

Comparative genomics analysis of WAK/WAKL family in Rosaceae identify candidate WAKs involved in the resistance to Botrytis cinerea

BACKGROUND

Wall associated kinase (WAK) and WAK-like (WAKL) are typical pattern recognition receptors act as the first sentry of plant defense. But little of WAK/WAKL family is known in Rosaceae.

RESULTS

In this study, 131 WAK/WAKL genes from apple, peach and strawberry were identified using a bioinformatics approach. Together with 68 RcWAK/RcWAKL in rose, we performed a comparative analysis of 199 WAK/WAKL in four Rosaceae crops. The phylogenetic analysis divided all the WAK/WAKL into five clades. Among them, the cis-elements of Clade II and Clade V promoters were enriched in jasmonic acid (JA) signaling and abiotic stress, respectively. And this can also be verified by the rose transcriptome responding to different hormone treatments. WAK/WAKL families have experienced a considerable proportion of purifying selection during evolution, but still 26 amino acid sites evolved under positive selection, which focused on extracellular conserved domains. WAK/WAKL genes presented collinearity relationship within and between crops, throughout four crops we mined four orthologous groups (OGs). The WAK/WAKL genes in OG1 and OG4 were speculated to involve in plant-Botrytis cinerea interaction, which were validated in rose via VIGS as well as strawberry by qRT-PCR.

CONCLUSIONS

These results not only provide genetic resources and valuable information for the evolutionary relationship of WAK/WAKL gene family, but also offer a reference for future in-depth studies of Rosaceae WAK/WAKL genes.

WheatCENet: A Database for Comparative Co-expression Networks Analysis of Allohexaploid Wheat and Its Progenitors

Genetic and epigenetic changes after polyploidization events could result in variable gene expression and modified regulatory networks. Here, using large-scale transcriptome data, we constructed co-expression networks for diploid, tetraploid, and hexaploid wheat species, and built a platform for comparing co-expression networks of allohexaploid wheat and its progenitors, named WheatCENet. WheatCENet is a platform for searching and comparing specific functional co-expression networks, as well as identifying the related functions of the genes clustered therein. Functional annotations like pathways, gene families, protein-protein interactions, microRNAs (miRNAs), and several lines of epigenome data are integrated into this platform, and Gene Ontology (GO) annotation, gene set enrichment analysis (GSEA), motif identification, and other useful tools are also included. Using WheatCENet, we found that the network of WHEAT ABERRANT PANICLE ORGANIZATION 1 (WAPO1) has more co-expressed genes related to spike development in hexaploid wheat than its progenitors. We also found a novel motif of CCWWWWWWGG (CArG) specifically in the promoter region of WAPO-A1, suggesting that neofunctionalization of the WAPO-A1 gene affects spikelet development in hexaploid wheat. WheatCENet is useful for investigating co-expression networks and conducting other analyses, and thus facilitates comparative and functional genomic studies in wheat. WheatCENet is freely available at http://bioinformatics.cpolar.cn/WheatCENet and http://bioinformatics.cau.edu.cn/WheatCENet.

Comparative Transcriptome Analysis Revealed the Freezing Tolerance Signaling Events in Winter Rapeseed (Brassica rapa L.)

Winter rapeseed (Brassica rapa L.) is an important oilseed crop in northwest China. Freezing stress severely limits its production and geographical distribution, and frequent extreme freezing events caused by climate change are increasing the chances of winter freeze-injury. However, the underlying mechanism of B. rapa response to freezing stress remains elusive. Here, B. rapa genome (v3.0) was used as a reference for the comparative transcriptomic analysis of Longyou 6 and Tianyou 2 (strong and weak cold tolerance, respectively) under different freezing stress. Before and after freezing stress, 5,982 and 11,630 unique differentially expressed genes (DEGs) between two cultivars were identified, respectively. After freezing stress, the GO terms in Tianyou 2 were mainly involved in "macromolecule biosynthetic process", and those in Longyou 6 were involved in "response to stimulus" and "oxidoreductase activity". Morphological and physiological results indicated that Longyou 6 retained a higher basal freezing resistance than Tinayou 2, and that cold acclimation could strengthen the basal freezing resistance. Freezing stress could activate the MAPK signal cascades, and the phosphorylation level of Longyou 6 showed a higher increase in response to freezing treatment than Tianyou 2. Based on our findings, it was speculated that the cell membrane of B. rapa perceives external signals under freezing stress, which are then transmitted to the nucleus through the cold-activated MAPK cascades and Ca²⁺-related protein kinase pathway, thus leading to activation of downstream target genes to enhance the freezing resistance of B. rapa.

Systems biology-based analysis indicates that PHO1;H10 positively modulates high light-induced anthocyanin biosynthesis in Arabidopsis leaves

Arabidopsis PHO1;H10 is a member of the PHO1 gene family with SPX and EXS domains, and its functions remain largely unknown. As shown in PCSD database, the upstream region of PHO1;H10 gene is in the active chromatin states, with high DHS accessibility and binding sites of multiple transcription factors, especially ABI5, SPCH and HY5. Co-expression network and data-mining analyses showed PHO1;H10 and co-expression genes were with activation under high light stress. We did wet-lab experiments, and found that the detached leaves of PHO1;H10 overexpression lines accumulated more anthocyanin than those of WT and mutant under high light treatment. RNA-seq results showed overexpression of PHO1;H10 up-regulated many anthocyanin biosynthetic genes. The GSEA analysis result showed that the functional module related to anthocyanin pathway was significantly enriched. In summary, we conducted systems biology approach, combining dry- and wet-lab analyses, and discovered that PHO1;H10 might play an essential role during modulating high light-induced anthocyanin accumulation in the Arabidopsis detached leaves.

GURFAP: A Platform for Gene Function Analysis in Glycyrrhiza Uralensis

Glycyrrhiza uralensis (Licorice), which belongs to Leguminosae, is famous for the function of pharmacologic action and natural sweetener with its dried roots and rhizomes. In recent years, the whole-genome sequence of G. uralensis has been completed, which will help to lay the foundation for the study of gene function. Here, we integrated the available genomic and transcriptomic data of G. uralensis and constructed the G. uralensis gene co-expression network. We then annotated gene functions of G. uralensis via aligning with public databases. Furthermore, gene families of G. uralensis were predicted by tools including iTAK (Plant Transcription factor and Protein kinase Identifier and Classifier), HMMER (hidden Markov models), InParanoid, and PfamScan. Finally, we constructed a platform for gene function analysis in G. uralensis (GURFAP, www.gzybioinfoormatics.cn/GURFAP). For analyzed and predicted gene function, we introduced various tools including BLAST (Basic local alignment search tool), GSEA (Gene set enrichment analysis), Motif, Heatmap, and JBrowse. Our analysis based on this platform indicated that the biosynthesis of glycyrrhizin might be regulated by MYB and bHLH. We also took CYP88D6, CYP72A154, and bAS gene in the synthesis pathway of glycyrrhizin as examples to demonstrate the reliability and availability of our platform. Our platform GURFAP will provide convenience for researchers to mine the gene function of G. uralensis and thus discover more key genes involved in the biosynthetic pathway of active ingredients.

GelFAP: Gene Functional Analysis Platform for Gastrodia elata

Gastrodia elata, also named Tianma, is a valuable traditional Chinese herbal medicine. It has numerous important pharmacological roles such as in sedation and lowering blood pressure and as anticonvulsant and anti-aging, and it also has effects on the immune and cardiovascular systems. The whole genome sequencing of G. elata has been completed in recent years, which provides a strong support for the construction of the G. elata gene functional analysis platform. Therefore, in our research, we collected and processed 39 transcriptome data of G. elata and constructed the G. elata gene co-expression networks, then we identified functional modules by the weighted correlation network analysis (WGCNA) package. Furthermore, gene families of G. elata were identified by tools including HMMER, iTAK, PfamScan, and InParanoid. Finally, we constructed a gene functional analysis platform for G. elata . In our platform, we introduced functional analysis tools such as BLAST, gene set enrichment analysis (GSEA), and cis-elements (motif) enrichment analysis tool. In addition, we analyzed the co-expression relationship of genes which might participate in the biosynthesis of gastrodin and predicted 19 mannose-binding lectin antifungal proteins of G. elata. We also introduced the usage of the G. elata gene function analysis platform (GelFAP) by analyzing CYP51G1 and GFAP4 genes. Our platform GelFAP may help researchers to explore the gene function of G. elata and make novel discoveries about key genes involved in the biological processes of gastrodin.

AppleMDO: A Multi-Dimensional Omics Database for Apple Co-Expression Networks and Chromatin States

As an economically important crop, apple is one of the most cultivated fruit trees in temperate regions worldwide. Recently, a large number of high-quality transcriptomic and epigenomic datasets for apple were made available to the public, which could be helpful in inferring gene regulatory relationships and thus predicting gene function at the genome level. Through integration of the available apple genomic, transcriptomic, and epigenomic datasets, we constructed co-expression networks, identified functional modules, and predicted chromatin states. A total of 112 RNA-seq datasets were integrated to construct a global network and a conditional network (tissue-preferential network). Furthermore, a total of 1,076 functional modules with closely related gene sets were identified to assess the modularity of biological networks and further subjected to functional enrichment analysis. The results showed that the function of many modules was related to development, secondary metabolism, hormone response, and transcriptional regulation. Transcriptional regulation is closely related to epigenetic marks on chromatin. A total of 20 epigenomic datasets, which included ChIP-seq, DNase-seq, and DNA methylation analysis datasets, were integrated and used to classify chromatin states. Based on the ChromHMM algorithm, the genome was divided into 620,122 fragments, which were classified into 24 states according to the combination of epigenetic marks and enriched-feature regions. Finally, through the collaborative analysis of different omics datasets, the online database AppleMDO (http://bioinformatics.cau.edu.cn/AppleMDO/) was established for cross-referencing and the exploration of possible novel functions of apple genes. In addition, gene annotation information and functional support toolkits were also provided. Our database might be convenient for researchers to develop insights into the function of genes related to important agronomic traits and might serve as a reference for other fruit trees.

[Development and evaluation of a system for pulse generation and signal processing control in a plasma atomic fluorescence spectrometer with hollow cathode lamp as excitation source]

A Z80 processor was replaced with a 386 microcomputer to control pulse generation for hollow cathode lamp and fluorescence signal processing in a Baird Plasma AFS-2000 atomic fluorescence spectrometer. An interface board was designed for control purpose. The system was evaluated with detection limits, relative standard deviations and linear regression coefficients of 15 elements, which were in good agreement with the corresponding values in the manufacturer's manual. With the higher-level microprocessor in the system, the performance of the instrument was improved.