ViroBLAST: a stand-alone BLAST web server for flexible queries of multiple databases and user's datasets

Analysis of HIV-1 subtypes and drug resistance mutations in people who inject drugs in Aizawl and nearby districts of Mizoram, India

Over the last decade, there has been a steady increase in HIV-1 prevalence in Mizoram, India. Importantly, this increase in HIV-1 prevalence is not only limited to the key population groups such as female sex workers (FSWs) or people who inject drugs (PWID), and has been witnessed in general population as well. Injecting drug use has long been one of the key drivers of HIV-1 epidemic across the north-eastern states of India. In this study, using HIV-1 pol gene region sequences from Aizawl and adjoining districts, we examined the HIV-1 subtypes, recombinant forms, drug resistance mutations and also the spatiotemporal dynamics of the potential unique recombinant forms. In our dataset, the dominant subtype was HIV-1 subtype C (94.91%). We could also identify the presence of CRF01_AE (1.69%) and BC recombinant forms (3.39%). Drug resistance mutation analysis revealed that resistance against non-nucleoside reverse transcriptase inhibitors was most common in the sequences having any resistance mutations. Evolutionary analysis of unique BC recombinants estimated the most recent common ancestor of these sequences around 2004-2005 and them having an ancestry of United States of America (USA) origin.

Genome-wide identification and analysis of the apple H+-ATPase gene family and its expression against iron deficiency stress

BACKGROUND

Iron (Fe) is indispensable to plants, playing a significant role in life activities such as respiration, chlorophyll biosynthesis, and photosynthetic electron transport. The decrease in pH caused by iron deficiency is related to the activation of H+-ATPase in the root plasma membrane, although the reaction mechanism of this enzyme is not well understood in apples. The H+-ATPase (HA) gene family has been extensively studied in Arabidopsis but is rarely reported in other species.

RESULTS

In this study, 14 HA genes were identified from the apple genome database through whole genome analysis. These apple H+-ATPase (MdHAs) genes were classified into four subsets (I, II, IV, V) based on phylogenetic analysis. Bioinformatics analysis revealed that these genes exhibited diversity in gene structure, chromosomal distribution, conserved motifs, and cis-acting elements. The qRT-PCR analysis revealed that iron deficiency stress significantly induced the upregulation of nine MdHA genes (MdHA5-MdHA14). Furthermore, in the roots of iron deficiency-resistant rootstock variety QZ1 compared to the non-resistant variety QZ2, the expression levels of nine genes (including MdHA1, MdHA6-MdHA13) in QZ2 were significantly lower than those in QZ1. This study lays a foundation for further study on the function of HA gene family.

CONCLUSIONS

These findings suggest that the HA gene family is involved in the apple's response to iron deficiency. This study not only lays the groundwork for future research but also highlights the role of MdHAs in iron stress tolerance.

VITALdb: to select the best viroinformatics tools for a desired virus or application

The recent pandemics of viral diseases, COVID-19/mpox (humans) and lumpy skin disease (cattle), have kept us glued to viral research. These pandemics along with the recent human metapneumovirus outbreak have exposed the urgency for early diagnosis of viral infections, vaccine development, and discovery of novel antiviral drugs and therapeutics. To support this, there is an armamentarium of virus-specific computational tools that are currently available. VITALdb (VIroinformatics Tools and ALgorithms database) is a resource of ~360 viroinformatics tools encompassing all major viruses (SARS-CoV-2, influenza virus, human immunodeficiency virus, papillomavirus, herpes simplex virus, hepatitis virus, dengue virus, Ebola virus, Zika virus, etc.) and several diverse applications [structural and functional annotation, antiviral peptides development, subspecies characterization, recognition of viral recombination, inhibitors identification, phylogenetic analysis, virus-host prediction, viral metagenomics, detection of mutation(s), primer designing, etc.]. Resources, tools, and other utilities mentioned in this article will not only facilitate further developments in the realm of viroinformatics but also provide tremendous fillip to translate fundamental knowledge into applied research. Most importantly, VITALdb is an inevitable tool for selecting the best tool(s) to carry out a desired task and hence will prove to be a vital database (VITALdb) for the scientific community. Database URL: https://compbio.iitr.ac.in/vitaldb.

Advances in bioinformatics and multi-omics integration: transforming viral infectious disease research in veterinary medicine

The world is changing due to factors like bioterrorism, massive environmental changes, globalization of trade and commerce, growing urbanization, changing climate, and pollution. Numerous diseases have emerged because of these factors, especially in companion and food-producing animals. Numerous pathogens have established themselves in naïve populations, harming reproduction, productivity, and health. Bioinformatics is considered a valuable tool in infectious disease research, as it provides a comprehensive overview of the identification of pathogens, their genetic makeup, and their evolutionary relationship. Therefore, there is an urgent need for a novel bioinformatics approach to help decipher and model viral epidemiology and informatics on domestic animals and livestock. With significant advancements in bioinformatics and NGS, researchers can now identify contigs, which are contiguous sequences of DNA that are assembled from overlapping fragments, assemble a complete genome, perform phylogenetic analysis to diagnose, investigate the risk of viral diseases in animals, handle and share large biological datasets across various species. Additionally, multi-omics data integration further deepens our understanding of homology, divergence, mutations, and evolutionary relationships, providing a comprehensive perspective on the molecular mechanisms driving animal pathogens infections. This review aims to reveal the importance of utilizing the multidisciplinary areas of bioinformatics, genomics, proteomics, transcriptomics, metabolomics, and metagenomics and their roles in studying viral infectious diseases in veterinary medicine that will eventually improve the health of animals.

scEccDNAdb: an integrated single-cell eccDNA resource for human and mouse

Extrachromosomal circular DNA (eccDNA), an extrachromosomal circular structured DNA, is extensively found in eukaryotes. Investigating eccDNA at the single-cell level is crucial for understanding cellular heterogeneity, evolution, development, and specific cellular functions. However, high-throughput identification methods for single-cell eccDNA are complex, and the lack of mature, widely applicable technologies has resulted in limited resources. To address this gap, we built scEccDNAdb, a database based on single-cell whole-genome sequencing data. It contains 3 195 464 single-cell eccDNA entries from human and mouse samples, with annotations including oncogenes, typical enhancers, super-enhancers, CCCTC-binding factor-binding sites, single nucleotide polymorphisms, chromatin accessibility, expression quantitative trait loci, transcription factor binding sites, motifs, and structural variants. Additionally, it provides nine online analysis and visualization tools, which enable the creation of publication-quality figures through user-uploaded files. Overall, scEccDNAdb is a comprehensive database for analyzing single-cell eccDNA data across diverse cell types, tissues, and species. Database URL: https://lcbb.swjtu.edu.cn/scEccDNAdb/.

Genomic reconstruction of unclassified microorganisms: Analysis of CRISPR arrays and genes involved in defense mechanisms

The constant battle between bacteria and viruses has led to the development of sophisticated antiviral defense strategies by bacteria to defend themselves against phages. This study analyzed a marshland metagenome to identify and characterize bacterial antiviral defense systems and phage interactions. We assembled 210 metagenome-assembled genomes (MAGs) from environmental DNA extracted from Pallikaranai marshland soil and 37 unclassified MAGs were filtered. MIMAG standards were followed, 2 high-quality and 15 medium-quality unclassified MAGs were picked. MINCED was used to identify 137 CRISPR arrays in the quality MAGs, and ViroBLAST was used to identify the phages that interact with the bacteria. About 242 spacer sequences were extracted from the CRISPR arrays, of which 54 had significant matches in the ViroBLAST database. 7 unverified bacteriophage species were also detected in the MAGs. The viral group of Caudoviricetes phage elements were identified as a frequent genome terminal repeat. The PADLOC identified 11 genes involved as a defense system in the MAGs. The PD-T4-6 defense system was found to be prevalent in 15 different unclassified MAGs. This study presents valuable insights intothe adaptations of unclassified bacteria to bacteriophages, as well as the genes used by these bacteria as a defense mechanism.

PungentDB: Bridging traditional Chinese medicine of medicine food homology and modern food flavor chemistry

Merging traditional Chinese medicine's (TCM) principles of medicine-food homology with modern flavor chemistry, this research unveils PungentDB (http://www.pungentdb.org.cn/home), a database documenting 205 unique pungent flavor compounds from 231 TCMs. It provides detailed insights into their chemical attributes, biological targets (including IC50/EC50 values), and molecular structures (2D/3D), enriched with visualizations of target organ distribution and protein structures, exploring the pungent flavor space with the help of a feature-rich visual interface. This collection, derived from over 3249 sources and highlighting 9129 targets, delves into the compounds' unique pungent flavors-taste, aroma, and thermal sensations-and their interaction with taste and olfactory receptors. PungentDB bridges ancient wisdom and culinary innovation, offering a nuanced exploration of pungent flavors' role in enhancing food quality, safety, and sensory experiences. This initiative propels flavor chemistry forward, serving as a pivotal resource for food science advancement and the innovative application of pungent flavors.

The genomes of all lungfish inform on genome expansion and tetrapod evolution

The genomes of living lungfishes can inform on the molecular-developmental basis of the Devonian sarcopterygian fish-tetrapod transition. We de novo sequenced the genomes of the African (Protopterus annectens) and South American lungfishes (Lepidosiren paradoxa). The Lepidosiren genome (about 91 Gb, roughly 30 times the human genome) is the largest animal genome sequenced so far and more than twice the size of the Australian (Neoceratodus forsteri)1 and African2 lungfishes owing to enlarged intergenic regions and introns with high repeat content (about 90%). All lungfish genomes continue to expand as some transposable elements (TEs) are still active today. In particular, Lepidosiren's genome grew extremely fast during the past 100 million years (Myr), adding the equivalent of one human genome every 10 Myr. This massive genome expansion seems to be related to a reduction of PIWI-interacting RNAs and C2H2 zinc-finger and Krüppel-associated box (KRAB)-domain protein genes that suppress TE expansions. Although TE abundance facilitates chromosomal rearrangements, lungfish chromosomes still conservatively reflect the ur-tetrapod karyotype. Neoceratodus' limb-like fins still resemble those of their extinct relatives and remained phenotypically static for about 100 Myr. We show that the secondary loss of limb-like appendages in the Lepidosiren-Protopterus ancestor was probably due to loss of sonic hedgehog limb-specific enhancers.

C2CDB: an advanced platform integrating comprehensive information and analysis tools of cancer-related circRNAs

Motivation

Circular RNAs (circRNAs) play important roles in gene expression and their involvement in tumorigenesis is emerging. circRNA-related database is a powerful tool for researchers to investigate circRNAs. However, existing databases lack advanced platform integrating comprehensive information and analysis tools of cancer-related circRNAs.

Results

We developed a comprehensive platform called CircRNA to Cancer Database (C2CDB), encompassing 318 158 cancer-related circRNAs expressed in tumors and adjacent tissues across 30 types of cancers. C2CDB provides basic details such as sequence and expression levels of circRNAs, as well as crucial insights into biological mechanisms, including miRNA binding, RNA-binding protein interaction, coding potential, base modification, mutation, and secondary structure. Moreover, C2CDB collects an extensive compilation of published literature on cancer circRNAs, extracting and presenting pivotal content encompassing biological functions, underlying mechanisms, and molecular tools in these studies. Additionally, C2CDB offers integrated tools to analyse three potential mechanisms: circRNA-miRNA ceRNA interaction, circRNA encoding, and circRNA biogenesis, facilitating investigators with convenient access to highly reliable information. To enhance clarity and organization, C2CDB has meticulously curated and integrated the previously chaotic nomenclature of circRNAs, addressing the prevailing confusion and ambiguity surrounding their designations.

Availability and implementation

C2CDB is freely available at http://pengyonglab.com/c2cdb.

Expansion and Functional Diversification of Long-Wavelength-Sensitive Opsin in Anabantoid Fishes

Gene duplication is one of the most important sources of novel genotypic diversity and the subsequent evolution of phenotypic diversity. Determining the evolutionary history and functional changes of duplicated genes is crucial for a comprehensive understanding of adaptive evolution. The evolutionary history of visual opsin genes is very dynamic, with repeated duplication events followed by sub- or neofunctionalization. While duplication of the green-sensitive opsins rh2 is common in teleost fish, fewer cases of multiple duplication events of the red-sensitive opsin lws are known. In this study, we investigate the visual opsin gene repertoire of the anabantoid fishes, focusing on the five lws opsin genes found in the genus Betta. We determine the evolutionary history of the lws opsin gene by taking advantage of whole-genome sequences of nine anabantoid species, including the newly assembled genome of Betta imbellis. Our results show that at least two independent duplications of lws occurred in the Betta lineage. The analysis of amino acid sequences of the lws paralogs of Betta revealed high levels of diversification in four of the seven transmembrane regions of the lws protein. Amino acid substitutions at two key-tuning sites are predicted to lead to differentiation of absorption maxima (λmax) between the paralogs within Betta. Finally, eye transcriptomics of B. splendens at different developmental stages revealed expression shifts between paralogs for all cone opsin classes. The lws genes are expressed according to their relative position in the lws opsin cluster throughout ontogeny. We conclude that temporal collinearity of lws expression might have facilitated subfunctionalization of lws in Betta and teleost opsins in general.

Differential expression of "Candidatus Liberibacter solanacearum" genes and prophage loci in different life stages of potato psyllid

Psyllid species, including the potato psyllid (PoP) Bactericera cockerelli (Sulc) (Triozidae) serve as host and vector of "Candidatus Liberibacter spp." ("Ca. Liberibacter"), which also infects diverse plant hosts, including citrus and tomato. Psyllid transmission of "Ca. Liberibacter" is circulative and propagative. The time of "Ca. Liberibacter" acquisition and therefore vector life stage most competent for bacterial transmission varies by pathosystems. Here, the potato psyllid-"Ca. Liberibacter solanacearum" (CLso) pathosystem was investigated to dissect CLso-prophage interactions in the tomato plant and PoP-psyllid host by real-time quantitative reverse transcriptase amplification of CLso genes/loci with predicted involvement in host infection and psyllid-CLso transmission. Genes/loci analyzed were associated with (1) CLso-adhesion, -invasion, -pathogenicity, and -motility, (2) prophage-adhesion and pathogenicity, and (3) CLso-lysogenic cycle. Relative gene expression was quantified by qRT-PCR amplification from total RNA isolated from CLso-infected 1st-2nd and 4th-5th nymphs and teneral adults and CLso-infected tomato plants in which CLso infection is thought to occur without SC1-SC2 replication. Gene/loci expression was host-dependent and varied with the psyllid developmental stage. Loci previously associated with repressor-anti-repressor regulation in the "Ca Liberibacter asiaticus"-prophage pathosystem, which maintains the lysogenic cycle in Asian citrus psyllid Diaphorina citri, were expressed in CLso-infected psyllids but not in CLso-infected tomato plants.

Morphological and dietary changes encoded in the genome of Beroe ovata, a ctenophore-eating ctenophore

As the sister group to all other animals, ctenophores (comb jellies) are important for understanding the emergence and diversification of numerous animal traits. Efforts to explore the evolutionary processes that promoted diversification within Ctenophora are hindered by undersampling genomic diversity within this clade. To address this gap, we present the sequence, assembly and initial annotation of the genome of Beroe ovata. Beroe possess unique morphology, behavior, ecology and development. Unlike their generalist carnivorous kin, beroid ctenophores feed exclusively on other ctenophores. Accordingly, our analyses revealed a loss of chitinase, an enzyme critical for the digestion of most non-ctenophore prey, but superfluous for ctenophorivores. Broadly, our genomic analysis revealed that extensive gene loss and changes in gene regulation have shaped the unique biology of B. ovata. Despite the gene losses in B. ovata, our phylogenetic analyses on photosensitive opsins and several early developmental regulatory genes show that these genes are conserved in B. ovata. This additional sampling contributes to a more complete reconstruction of the ctenophore ancestor and points to the need for extensive comparisons within this ancient and diverse clade of animals. To promote further exploration of these data, we present BovaDB (http://ryanlab.whitney.ufl.edu/bovadb/), a portal for the B. ovata genome.

Core proteome mediated subtractive approach for the identification of potential therapeutic drug target against the honeybee pathogen Paenibacillus larvae

Background & Objectives

American foulbrood (AFB), caused by the highly virulent, spore-forming bacterium Paenibacillus larvae, poses a significant threat to honey bee brood. The widespread use of antibiotics not only fails to effectively combat the disease but also raises concerns regarding honey safety. The current computational study was attempted to identify a novel therapeutic drug target against P. larvae, a causative agent of American foulbrood disease in honey bee.

Methods

We investigated effective novel drug targets through a comprehensive in silico pan-proteome and hierarchal subtractive sequence analysis. In total, 14 strains of P. larvae genomes were used to identify core genes. Subsequently, the core proteome was systematically narrowed down to a single protein predicted as the potential drug target. Alphafold software was then employed to predict the 3D structure of the potential drug target. Structural docking was carried out between a library of phytochemicals derived from traditional Chinese flora (n > 36,000) and the potential receptor using Autodock tool 1.5.6. Finally, molecular dynamics (MD) simulation study was conducted using GROMACS to assess the stability of the best-docked ligand.

Results

Proteome mining led to the identification of Ketoacyl-ACP synthase III as a highly promising therapeutic target, making it a prime candidate for inhibitor screening. The subsequent virtual screening and MD simulation analyses further affirmed the selection of ZINC95910054 as a potent inhibitor, with the lowest binding energy. This finding presents significant promise in the battle against P. larvae.

Conclusions

Computer aided drug design provides a novel approach for managing American foulbrood in honey bee populations, potentially mitigating its detrimental effects on both bee colonies and the honey industry.

TriticeaeSSRdb: a comprehensive database of simple sequence repeats in Triticeae

Microsatellites, known as simple sequence repeats (SSRs), are short tandem repeats of 1 to 6 nucleotide motifs found in all genomes, particularly eukaryotes. They are widely used as co-dominant markers in genetic analyses and molecular breeding. Triticeae, a tribe of grasses, includes major cereal crops such as bread wheat, barley, and rye, as well as abundant forage and lawn grasses, playing a crucial role in global food production and agriculture. To enhance genetic work and expedite the improvement of Triticeae crops, we have developed TriticeaeSSRdb, an integrated and user-friendly database. It contains 3,891,705 SSRs from 21 species and offers browsing options based on genomic regions, chromosomes, motif types, and repeat motif sequences. Advanced search functions allow personalized searches based on chromosome location and length of SSR. Users can also explore the genes associated with SSRs, design customized primer pairs for PCR validation, and utilize practical tools for whole-genome browsing, sequence alignment, and in silico SSR prediction from local sequences. We continually update TriticeaeSSRdb with additional species and practical utilities. We anticipate that this database will greatly facilitate trait genetic analyses and enhance molecular breeding strategies for Triticeae crops. Researchers can freely access the database at http://triticeaessrdb.com/.

The chromosome-level genome and functional database accelerate research about biosynthesis of secondary metabolites in Rosa roxburghii

Rosa roxburghii Tratt, a valuable plant in China with long history, is famous for its fruit. It possesses various secondary metabolites, such as L-ascorbic acid (vitamin C), alkaloids and poly saccharides, which make it a high nutritional and medicinal value. Here we characterized the chromosome-level genome sequence of R. roxburghii, comprising seven pseudo-chromosomes with a total size of 531 Mb and a heterozygosity of 0.25%. We also annotated 45,226 coding gene loci after masking repeat elements. Orthologs for 90.1% of the Complete Single-Copy BUSCOs were found in the R. roxburghii annotation. By aligning with protein sequences from public platform, we annotated 85.89% genes from R. roxburghii. Comparative genomic analysis revealed that R. roxburghii diverged from Rosa chinensis approximately 5.58 to 13.17 million years ago, and no whole-genome duplication event occurred after the divergence from eudicots. To fully utilize this genomic resource, we constructed a genomic database RroFGD with various analysis tools. Otherwise, 69 enzyme genes involved in L-ascorbate biosynthesis were identified and a key enzyme in the biosynthesis of vitamin C, GDH (L-Gal-1-dehydrogenase), is used as an example to introduce the functions of the database. This genome and database will facilitate the future investigations into gene function and molecular breeding in R. roxburghii.

DhuFAP: a platform for gene functional analysis in Dendrobium huoshanense

BACKGROUND

Dendrobium huoshanense, a traditional medicinal and food plant, has a rich history of use. Recently, its genome was decoded, offering valuable insights into gene function. However, there is no comprehensive gene functional analysis platform for D. huoshanense.

RESULT

To address this, we created a platform for gene function analysis and comparison in D. huoshanense (DhuFAP). Using 69 RNA-seq samples, we constructed a gene co-expression network and annotated D. huoshanense genes by aligning sequences with public protein databases. Our platform contained tools like Blast, gene set enrichment analysis, heatmap analysis, sequence extraction, and JBrowse. Analysis revealed co-expression of transcription factors (C2H2, GRAS, NAC) with genes encoding key enzymes in alkaloid biosynthesis. We also showcased the reliability and applicability of our platform using Chalcone synthases (CHS).

CONCLUSION

DhuFAP ( www.gzybioinformatics.cn/DhuFAP ) and its suite of tools represent an accessible and invaluable resource for researchers, enabling the exploration of functional information pertaining to D. huoshanense genes. This platform stands poised to facilitate significant biological discoveries in this domain.

Genome sequencing and protein modeling unraveled the 2AP biosynthesis in Bacillus cereus DB25

2-Acetyl-1-pyrroline (2AP) is an important and major flavor aroma compound responsible for the fragrance of basmati rice, cheese, wine, and several other food products. Biosynthesis of 2AP in aromatic rice and a few other plant species is associated with a recessive Betaine aldehyde dehydrogenase 2 (BADH2) gene. However, the literature is scant on the relationship between the functional BADH2 gene and 2AP biosynthesis in prokaryotic systems. Therefore, in the present study, we aimed to explore the functionality of the BADH2 gene for 2AP biosynthesis in 2AP synthesizing rice rhizobacterial isolate Bacillus cereus DB25 isolated from the rhizosphere of basmati rice (Oryza sativa L.). Full-length BcBADH2 sequence was obtained through whole genome sequencing (WGS) and further confirmed through traditional PCR and Sanger sequencing. Then the functionality of the BcBADH2 gene was evaluated in-silico through bioinformatics analysis and protein docking studies and further experimentally validated through enzyme assay. The sequencing and bioinformatics analysis results revealed a full-length 1485 bp BcBADH2 coding sequence without any deletion or premature stop codons. Full-length BcBADH2 was found to encode a fully functional protein of 54.08 kDa with pI of 5.22 and showed the presence of the conserved amino acids responsible for enzyme activity. The docking studies confirmed a good affinity between the protein and its substrate whereas the presence of BcBADH2 enzyme activity confirmed the functionality of BADH2 enzyme in B. cereus DB25. In conclusion, the findings of the present study suggest that B. cereus DB25 is able to synthesize 2AP despite a functional BADH2 gene and there may be a different molecular mechanism responsible for 2AP biosynthesis in bacterial systems, unlike that found in aromatic rice and other eukaryotic plant species.

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.

AgingReG: a curated database of aging regulatory relationships in humans

Aging and cellular senescence are characterized by a progressive loss of physiological integrity, which could be triggered by aging factors such as physiological, pathological and external factors. Numerous studies have shown that gene regulatory events play crucial roles in aging, increasing the need for a comprehensive repository of regulatory relationships during aging. Here, we established a manually curated database of aging factors (AgingReG, https://bio.liclab.net/Aging-ReG/), focusing on the regulatory relationships during aging with experimental evidence in humans. By curating thousands of published literature, 2157 aging factor entries (1345 aging gene entries, 804 external factor entries and eight aging-related pathway entries) and related regulatory information were manually curated. The regulatory relationships were classified into four types according to their functions: (i) upregulation, which indicates that aging factors upregulate the expression of target genes during aging; (ii) downregulation, which indicates that aging factors downregulate the expression of target genes during aging; (iii) activation, which indicates that aging factors influence the activity of target genes during aging and (iv) inhibition, which indicates that aging factors inhibit the activation of target molecule activity, leading to declined or lost target activity. AgingReG involves 651 upregulating pairs, 632 downregulating pairs, 330 activation-regulating pairs and 34 inhibition-regulating pairs, covering 195 disease types and more than 800 kinds of cells and tissues from 1784 published literature studies. AgingReG provides a user-friendly interface to query, browse and visualize detailed information about the regulatory relationships during aging. We believe that AgingReG will serve as a valuable resource database in the field of aging research. Database URL: https://bio.liclab.net/Aging-ReG/.

The pear genomics database (PGDB): a comprehensive multi-omics research platform for Pyrus spp

BACKGROUND

Pears are among the most important temperate fruit trees in the world, with significant research efforts increasing over the last years. However, available omics data for pear cannot be easily and quickly retrieved to enable further studies using these biological data.

DESCRIPTION

Here, we present a publicly accessible multi-omics pear resource platform, the Pear Genomics Database (PGDB). We collected and collated data on genomic sequences, genome structure, functional annotation, transcription factor predictions, comparative genomics, and transcriptomics. We provide user-friendly functional modules to facilitate querying, browsing and usage of these data. The platform also includes basic and useful tools, including JBrowse, BLAST, phylogenetic tree building, and additional resources providing the possibility for bulk data download and quick usage guide services.

CONCLUSIONS

The Pear Genomics Database (PGDB, http://pyrusgdb.sdau.edu.cn ) is an online data analysis and query resource that integrates comprehensive multi-omics data for pear. This database is equipped with user-friendly interactive functional modules and data visualization tools, and constitutes a convenient platform for integrated research on pear.

Phylogenomics of transcriptionally active AP2/ERF and bHLH transcription factors and study of their promoter regions in Nothapodytes nimmoniana (J.Graham) Mabb

Nothapodytes nimmoniana is a medicinally important plant producing anticancer monoterpene indole alkaloid (MIA), camptothecin (CPT). The CPT is synthesised through the strictosidine intermediate following the MIA pathway; however, transcriptional regulation of CPT pathway is still elusive in N. nimmoniana. Biosynthesis of MIA is regulated by various transcription factors (TFs) belonging to AP2/ERF, bHLH, MYB, and WRKY families. The present study identified transcriptionally active full-length 105 AP2/ERF and 68 bHLH family TFs from the N. nimmoniana. AP2/ERF TFs were divided into three subfamilies along with a soloist, while bHLH TFs were divided into 10 subfamilies according to their phylogenetic similarities. Three group IXa ERFs, Nn-ERF22, Nn-ERF29, and Nn-ERF41, one subfamily IVa TF Nn-bHLH7, and three subfamilies IIIe Nn-bHLH33, Nn-bHLH51, and Nn-bHLH52 clustered with the TFs regulating alkaloid biosynthesis in Catharanthus roseus, tomato, tobacco, and Artemisia annua. Expression of these TFs in N. nimmoniana was higher in roots, which is a primary CPT accumulating tissue. Moreover, genome skimming approach was used to reconstruct the promoter regions of candidate ERF genes to identify the cis-regulatory elements. The presence of G-boxes and other jasmonic acid-responsive elements in the promoter suggests the regulation of ERFs by bHLHs. The present study effectively generated and used genomics resource for characterisation of regulatory TFs from non-model medicinal plant.

Interplay between transcriptional regulators and VapBC toxin-antitoxin loci during thermal stress response in extremely thermoacidophilic archaea

Thermoacidophilic archaea lack sigma factors and the large inventory of heat shock proteins (HSPs) widespread in bacterial genomes, suggesting other strategies for handling thermal stress are involved. Heat shock transcriptomes for the thermoacidophilic archaeon Saccharolobus (f. Sulfolobus) solfataricus 98/2 revealed genes that were highly responsive to thermal stress, including transcriptional regulators YtrASs (Ssol_2420) and FadRSs (Ssol_0314), as well as type II toxin-antitoxin (TA) loci VapBC6 (Ssol_2337, Ssol_2338) and VapBC22 (Ssol_0819, Ssol_0818). The role, if any, of type II TA loci during stress response in microorganisms, such as Escherichia coli, is controversial. But, when genes encoding YtrASs , FadRSs , VapC22, VapB6, and VapC6 were systematically mutated in Sa. solfataricus 98/2, significant up-regulation of the other genes within this set was observed, implicating an interconnected regulatory network during thermal stress response. VapBC6 and VapBC22 have close homologues in other Sulfolobales, as well as in other archaea (e.g. Pyrococcus furiosus and Archaeoglobus fulgidus), and their corresponding genes were also heat shock responsive. The interplay between VapBC TA loci and heat shock regulators in Sa solfataricus 98/2 not only indicates a cellular mechanism for heat shock response that differs from bacteria but one that could have common features within the thermophilic archaea.

Complete Genome Sequence of Bacteriophage EO1, Which Infects Both Escherichia coli O157:H7 and Shigella

The lytic bacteriophage EO1 has been newly isolated. This phage infects Escherichia coli O157:H7 and has a broad antibacterial spectrum, including against Shigella. The complete genome sequence of phage EO1 was determined; its full length is 166,941 bp, and it has a G+C content of 35.46%.

MDSi: Multi-omics Database for Setaria italica

BACKGROUND

Foxtail millet (Setaria italica) harbors the small diploid genome (~ 450 Mb) and shows the high inbreeding rate and close relationship to several major foods, feed, fuel and bioenergy grasses. Previously, we created a mini foxtail millet, xiaomi, with an Arabidopsis-like life cycle. The de novo assembled genome data with high-quality and an efficient Agrobacterium-mediated genetic transformation system made xiaomi an ideal C₄ model system. The mini foxtail millet has been widely shared in the research community and as a result there is a growing need for a user-friendly portal and intuitive interface to perform exploratory analysis of the data.

RESULTS

Here, we built a Multi-omics Database for Setaria italica (MDSi, http://sky.sxau.edu.cn/MDSi.htm ), that contains xiaomi genome of 161,844 annotations, 34,436 protein-coding genes and their expression information in 29 different tissues of xiaomi (6) and JG21 (23) samples that can be showed as an Electronic Fluorescent Pictograph (xEFP) in-situ. Moreover, the whole-genome resequencing (WGS) data of 398 germplasms, including 360 foxtail millets and 38 green foxtails and the corresponding metabolic data were available in MDSi. The SNPs and Indels of these germplasms were called in advance and can be searched and compared in an interactive manner. Common tools including BLAST, GBrowse, JBrowse, map viewer, and data downloads were implemented in MDSi.

CONCLUSION

The MDSi constructed in this study integrated and visualized data from three levels of genomics, transcriptomics and metabolomics, and also provides information on the variation of hundreds of germplasm resources that can satisfies the mainstream requirements and supports the corresponding research community.

CRISPR detection in metagenome-assembled genomes (MAGs) of coal mine

Bacterial and archaeal CRISPR-Cas systems provide adaptive immune protection against foreign mobile genetic elements. When viruses infect bacteria, a small portion of the viral DNA is inserted into the bacterial DNA in a specific pattern to produce segments known as CRISPR arrays. Metagenome assembled genomes (MAGs) were used in our study to identify the CRISPR sequence for determining the interacted phage. Metagenomic data from a coal mine was used to perform a computational study. From raw reads, 206151 contigs were assembled. Then contigs were clustered into 150 Metagenome assembled genomes from which 78 non-redundant MAGs were selected. Using the CHECKM standard, seven MAGs were found to have >80 completeness and <20 contaminations. Those MAGs were analyzed for the presence of CRISPR elements. Out of seven MAGs, four MAGs have the CRISPR elements and are searched against the VIROblast database. CRISPR arrays have 4, 1, 3, and 7 spacer sequences in the MAGs of Burkholderia, Acinetobacter, Oxalobacteraceae, and Burkholderia multivorans respectively. The uncultured Caudovirales phage genomic regions were present in the genomes of Burkholderia, Oxalobacteriaceae, and Burkholderia multivorans. This study follows the unconventional metagenomics workflow to provide a better understanding of bacteria and phage interactions.

GelFAP v2.0: an improved platform for Gene functional analysis in Gastrodia elata

BACKGROUND

Gastrodia elata (tianma), a well-known medicinal orchid, is widely used to treat various kinds of diseases with its dried tuber. In recent years, new chromosome-level genomes of G.elata have been released in succession, which offer an enormous resource pool for understanding gene function. Previously we have constructed GelFAP for gene functional analysis of G.elata. As genomes are updated and transcriptome data is accumulated, collection data in GelFAP cannot meet the need of researchers.

RESULTS

Based on new chromosome-level genome and transcriptome data, we constructed co-expression network of G. elata, and then we annotated genes by aligning with sequences from NR, TAIR, Uniprot and Swissprot database. GO (Gene Ontology) and KEGG (Kyoto Encylopaedia of Genes and Genomes) annotations were predicted by InterProScan and GhostKOALA software. Gene families were further predicted by iTAK (Plant Transcription factor and Protein kinase Identifier and Classifier), HMMER (hidden Markov models), InParanoid. Finally, we developed an improved platform for gene functional analysis in G. elata (GelFAP v2.0) by integrating new genome, transcriptome data and processed functional annotation. Several tools were also introduced to platform including BLAST (Basic Local Alignment Search Tool), GSEA (Gene Set Enrichment Analysis), Heatmap, JBrowse, Motif analysis and Sequence extraction. Based on this platform, we found that the flavonoid biosynthesis might be regulated by transcription factors (TFs) such as MYB, HB and NAC. We also took C4H and GAFP4 as examples to show the usage of our platform.

CONCLUSION

An improved platform for gene functional analysis in G. elata (GelFAP v2.0, www.gzybioinformatics.cn/Gelv2 ) was constructed, which provides better genome data, more transcriptome resources and more analysis tools. The updated platform might be preferably benefit researchers to carry out gene functional research for their project.

BarleyExpDB: an integrative gene expression database for barley

BACKGROUND

RNA-sequencing (RNA-seq) has been widely used to study the dynamic expression patterns of transcribed genes, which can lead to new biological insights. However, processing and analyzing these huge amounts of histological data remains a great challenge for wet labs and field researchers who lack bioinformatics experience and computational resources.

RESULTS

We present BarleyExpDB, an easy-to-operate, free, and web-accessible database that integrates transcriptional profiles of barley at different growth and developmental stages, tissues, and stress conditions, as well as differential expression of mutants and populations to build a platform for barley expression and visualization. The expression of a gene of interest can be easily queried by searching by known gene ID or sequence similarity. Expression data can be displayed as a heat map, along with functional descriptions as well as Gene Ontology, Kyoto Encyclopedia of Genes and Genomes, Proteins Families Database, and Simple Modular Architecture Research Tool annotations.

CONCLUSIONS

BarleyExpDB will serve as a valuable resource for the barley research community to leverage the vast publicly available RNA-seq datasets for functional genomics research and crop molecular breeding.

eccDNA Atlas: a comprehensive resource of eccDNA catalog

Extrachromosomal circular DNA (eccDNA) represents a large category of non-mitochondrial and non-plasmid circular extrachromosomal DNA, playing an indispensable role in various aspects such as tumorigenesis, immune responses. However, the information of characteristics and functions about eccDNA is fragmented, hiding behind abundant literatures and massive whole-genome sequencing (WGS) data, which has not been sufficiently used for the identification of eccDNAs. Therefore, establishing an integrated repository portal is essential for identifying and analyzing eccDNAs. Here, we developed eccDNA Atlas (http://lcbb.swjtu.edu.cn/eccDNAatlas), a user-friendly database of eccDNAs that aims to provide a high-quality and integrated resource for browsing, searching and analyzing eccDNAs from multiple species. eccDNA Atlas currently contains 629 987 eccDNAs and 8221 ecDNAs manually curated from literatures and 1105 ecDNAs predicted by AmpliconArchitect based on WGS data involved in 66 diseases, 57 tissues and 319 cell lines. The content of each eccDNA entry includes multiple aspects such as sequence, disease, function, characteristic, validation strategies. Furthermore, abundant annotations and analyzing utilities were provided to explore existing eccDNAs in eccDNA Atlas or user-defined eccDNAs including oncogenes, typical enhancers, super enhancers, CTCF-binding sites, SNPs, chromatin accessibility, eQTLs, gene expression, survival and genome visualization. Overall, eccDNA Atlas provides an integrated eccDNA data warehouse and serves as an important tool for future research.

Novel configurations of type I-E CRISPR-Cas system in Corynebacterium striatum clinical isolates

Clustered regularly interspaced short palindromic repeats (CRISPR) are a prokaryotic adaptive immune system that, through Cas proteins, promote the degradation of foreign nucleic acids such as phages and plasmids. We analyzed 10 genomes of Corynebacterium striatum clinical isolates from a public hospital in Rio de Janeiro, Brazil, the most emergent multidrug-resistant Corynebacterium species. All isolates were submitted to antimicrobial susceptibility testing. The occurrence and diversity of the CRISPR system were investigated by bioinformatics tools. Our analysis revealed that the isolates exhibited type I-E gene arrangements, and 3 more multidrug-resistant isolates, alternative type I-E gene arrangements, showing a divergent gene arrangement within the cas operon. Phylogenetic analysis of the cas1 gene of this type I-E CRISPR-Cas system alternative arrangement, termed here type I-E', showed a cluster in a distinct clade of the type I-E CRISPR-Cas system. The systems' guanine-cytosine (GC) content is lower than the genomic DNA's GC content, and mobile genetic elements were found in some isolates near the CRISPR-Cas system. Most CRISPR spacers are unknown indicating that there is a reservoir of unexplored corynebacteriophages and plasmids. Some spacers showed perfect homologies with phage and plasmid sequences. Intact phage regions were found in 3 of our isolates, ranging from 9.1 to 43.8 kb, with regions showing similarity to Rhodococcus and Corynebacterium phages. Our results may contribute to research about the CRISPR-Cas system diversity in C. striatum, where there are no published data to date.

Artificial microRNA mediated silencing of cyclase and aldo-keto reductase genes reveal their involvement in the plumbagin biosynthetic pathway

Plumbagin and other naphthoquinone derivatives from the Plumbago zeylanica L. (Plumbaginaceae) are known for their anticancer and other medicinal properties. Previous reports suggest that 3-methyl-1,8-naphthalene-diol is an intermediate of the plumbagin biosynthetic pathway and is synthesized from hexaketide backbone; a reaction catalyzed by type III polyketide synthase (PKS) along with certain accessory enzymes. Our earlier transcriptomic and metabolomic studies suggest that along with PKS, putative cyclase and aldo-keto reductase might be involved in the formation of 3-methyl-1,8-naphthalene-diol. The present study probed young leaf transcriptome and identified cyclase and aldo-keto reductase like transcripts that might be involved in the intramolecular aldol condensation of hexaketide intermediate and decarboxylation, carbonyl reduction and hydroxyl elimination of keto or enol forms of hexaketide intermediates respectively. Moreover, sequence alignment of identified cyclase1 possesses signature β-α-β-β-α-α-β topology, which belongs to the dimeric α + β barrel (DABB) protein family and is involved in the C2-C11 and C4-C9 intramolecular aldol condensation of hexaketide intermediates. Along with cyclase1, we further identified and characterized P. zeylanica specific aldo-keto reductase1 (AKR1) which is a novel member of the aldo-keto reductase (AKR) multi-gene family that possesses the conserved Asp60, Tyr65, Lys91, and His132 residues and is proposed to be involved in the C1 decarboxylation, C3 carbonyl reduction and C7 hydroxyl elimination of keto or enol form of hexaketide intermediate to form 3-methyl-1,8-naphthalene-diol. Further, the functional characterization using the artificial microRNA mediated transient silencing approach confirmed the involvement of cyclase1 and AKR1 in the plumbagin biosynthetic pathway. This is the first study reporting the identification and functional characterization of cyclase1 and AKR1 genes involved in the plumbagin biosynthetic pathway and general plant polyketide biosynthesis.

Isolation and Characterization of Lytic Bacteriophages Active against Clinical Strains of E. coli and Development of a Phage Antimicrobial Cocktail

Pathogenic E. coli cause urinary tract, soft tissue and central nervous system infections, sepsis, etc. Lytic bacteriophages can be used to combat such infections. We investigated six lytic E. coli bacteriophages isolated from wastewater. Transmission electron microscopy and whole genome sequencing showed that the isolated bacteriophages are tailed phages of the Caudoviricetes class. One-step growth curves revealed that their latent period of reproduction is 20-30 min, and the average value of the burst size is 117-155. During co-cultivation with various E. coli strains, the phages completely suppressed bacterial host culture growth within the first 4 h at MOIs 10^-7 to 10^-3. The host range lysed by each bacteriophage varied from six to two bacterial strains out of nine used in the study. The cocktail formed from the isolated bacteriophages possessed the ability to completely suppress the growth of all the E. coli strains used in the study within 6 h and maintain its lytic activity for 8 months of storage. All the isolated bacteriophages may be useful in fighting pathogenic E. coli strains and in the development of phage cocktails with a long storage period and high efficiency in the treatment of bacterial infections.

Functional genomics-enabled characterization of CYP81B140 and CYP81B141 from Plumbago zeylanica L. substantiates their involvement in plumbagin biosynthesis

Novel cytochrome P450s, CYP81B140 and CYP81B141 from Plumbago zeylanica were functionally characterized to understand their involvement in polyketide plumbagin biosynthesis. Further, we propose 3-methyl-1-8-naphthalenediol and isoshinanolone as intermediates for plumbagin biosynthesis. Plumbago zeylanica L. (P. zeylanica) is a medicinally important plant belonging to the family Plumbaginaceae. It comprises the most abundant naphthoquinone plumbagin having anti-cancer activity. Only the polyketide synthase (PKS) enzyme has been identified from the biosynthetic pathway which catalyzes iterative condensation of acetyl-CoA and malonyl-CoA molecules. The plumbagin biosynthesis involves hydroxylation, oxidation, hydration and dehydration of intermediate compounds which are expected to be catalyzed by cytochrome P450s (CYPs). To identify the CYPs, co-expression analysis was carried out using PKS as a candidate gene. Out of the eight identified CYPs, CYP81B140 and CYP81B141 have similar expression with PKS and belong to the CYP81 family. Phylogenetic analysis suggested that CYP81B140 and CYP81B141 cluster with CYPs from CYP81B, CYP81D, CYP81E and CYP81AA subfamilies which are known to be involved in the hydroxylation and oxidation reactions. Moreover, artificial microRNA-mediated transient individual silencing and co-silencing of CYP81B140 and CYP81B141 significantly reduced plumbagin and increased the 3-methyl-1-8-naphthalenediol and isoshinanolone content. Based on metabolite analysis, we proposed that 3-methyl-1-8-naphthalenediol and isoshinanolone function as intermediates for plumbagin biosynthesis. Transient silencing, over-expression and docking analysis revealed that CYP81B140 is involved in C-1 oxidation, C-4 hydroxylation and [C2-C3] hydration of 3-methyl-1-8-naphthalenediol to form isoshinanolone, whereas CYP81B141 is catalyzing [C2-C3] dehydration and C-4 oxidation of isoshinanolone to form plumbagin. Our results indicated that both CYP81B140 and CYP81B141 are promiscuous and necessary for plumbagin biosynthesis. This is the first report of identification and functional characterization of P. zeylanica-specific CYPs involved in plumbagin biosynthetic pathway and in general hexaketide synthesis in plants.

Genome-wide identification and expression analysis of PIN gene family under phytohormone and abiotic stresses in Vitis Vinifera L

The auxin efflux transport proteins PIN-formed (PIN) has wide adaptability to hormone and abiotic stress, but the response mechanism of PINs in grape remains unclear. In this study, 12 members of VvPINs were identified and distributed on 8 chromosomes. The PIN genes of five species were divided into two subgroups, and the similarity of exons/introns and motifs of VvPIN genes were found in the same subgroup. Meanwhile, according to the examination of conserved motifs, the motif 3 included the conserved structure NPNTY. The promoter region of VvPIN gene family contained various cis-acting elements, which were related to light, abiotic stress, and hormones which are essential for growth and development. Additionally, VvPIN1, VvPIN9, and VvPIN11 proteins simultaneously interacted with the ARF, ABC, PINOID, GBF1, and VIT_08s0007g09010. The results of qRT-PCR revealed that the majority of the VvPINs were highly induced by NAA, GA3, ABA, MeJA, SA, NaCl, low-temperature (4 ℃), and PEG treatments, and the results were consistent with the prediction of the cis-acting elements. Moreover, the expression profile and quantitative real-time PCR (qRT-PCR) demonstrated that VvPIN genes were expressed in roots, stems, and leaves. The subcellular localization of VvPIN1 in Nicotiana benthamiana revealed that VvPIN1 was localized at the plasma membrane. Collectively, this study revealed that PIN genes could respond to various abiotic stresses, and provided a framework for regulating the expression of PIN genes to enhance the resistance of the grape.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-022-01239-8.

The TCP transcription factor HvTB2 heterodimerizes with VRS5 and controls spike architecture in barley

Understanding the molecular network, including protein-protein interactions, of VRS5 provide new routes towards the identification of other key regulators of plant architecture in barley. The TCP transcriptional regulator TEOSINTE BRANCHED 1 (TB1) is a key regulator of plant architecture. In barley, an important cereal crop, HvTB1 (also referred to as VULGARE SIX-ROWED spike (VRS) 5), inhibits the outgrowth of side shoots, or tillers, and grains. Despite its key role in barley development, there is limited knowledge on the molecular network that is utilized by VRS5. In this work, we performed protein-protein interaction studies of VRS5. Our analysis shows that VRS5 potentially interacts with a diverse set of proteins, including other class II TCP's, NF-Y TF, but also chromatin remodelers. Zooming in on the interaction capacity of VRS5 with other TCP TFs shows that VRS5 preferably interacts with other class II TCP TFs in the TB1 clade. Induced mutagenesis through CRISPR-Cas of one of the putative VRS5 interactors, HvTB2 (also referred to as COMPOSITUM 1 and BRANCHED AND INDETERMINATE SPIKELET 1), resulted in plants that have lost their characteristic unbranched spike architecture. More specifically, hvtb2 mutants exhibited branches arising at the main spike, suggesting that HvTB2 acts as inhibitor of branching. Our protein-protein interaction studies of VRS5 resulted in the identification of HvTB2 as putative interactor of VRS5, another key regulator of spike architecture in barley. The study presented here provides a first step to underpin the protein-protein interactome of VRS5 and to identify other, yet unknown, key regulators of barley plant architecture.

EfGD: the Erianthus fulvus genome database

Erianthus fulvus (TaxID: 154759) is a valuable germplasm resource in sugarcane breeding and research and has excellent agronomic traits, such as drought resistance, cold resistance, barren tolerance and high brix. With a stable chromosome number (2n = 20) and a small genome (0.9 Gb), it is an ideal candidate for research on sugarcane. Next-generation sequencing technology has enabled a growing number of studies to focus on genomics. Due to the large amount of omics data available, a centralized platform is necessary for ensuring the consistency, independence and maintainability of these large-scale datasets through storage, analysis and integration. Here, we present a comprehensive database for the E. fulvus genome, EfGD. By using the new high-quality reference genome and its annotations, the EfGD provides the largest whole-genome sequencing reference dataset for E. fulvus, which archives 27 165 protein-coding genes and 55 564 488 SNPs from 202 newly resequenced genomes. Furthermore, we created a user-friendly graphical interface for visualizing genomic diversity, population structure and evolution and provided other tools on an open platform.

Database URL: https://efgenome.ynau.edu.cn

Draft genome sequences data of two Rosemountvirus phages isolated from soil near poultry farm

Here we report the draft complete genome sequence of Escherichia phage vB_EcoM_IntR and Escherichia phage vB_EcoM_PiR isolated from a poultry farm soil sample. Bacteriophage isolation on the model of the pathogenic for birds Escherichia coli was carried out. Sequencing was carried out on the Illumina platform. Phage genomes were assembled using a software Geneious package. Phage sequences do not contain tRNA. Phages are a linear double-stranded DNA from Rosemountvirus genus of Myoviridae family with a genome 52,782 and 52,936 bp, respectively, containing 71 predicted open reading frames (ORFs). The phages have a standard icosahedral head and a contractive tail characteristic of the myovirus family.

High-Resolution Mapping of Barley mild mosaic virus Resistance Gene rym15

Barley yellow mosaic virus (BaYMV) and Barley mild mosaic virus (BaMMV), which are transmitted by the soil-borne plasmodiophorid Polymyxa graminis, cause high yield losses in barley. In previous studies, the recessive BaMMV resistance gene rym15, derived from the Japanese landrace Chikurin Ibaraki 1, was mapped on chromosome 6HS of Hordeum vulgare. In this study, 423 F4 segmental recombinant inbred lines (RILs) were developed from crosses of Chikurin Ibaraki 1 with two BaMMV-susceptible cultivars, Igri (139 RILs) and Uschi (284 RILs). A set of 32 competitive allele-specific PCR (KASP) assays, designed using single nucleotide polymorphisms (SNPs) from the barley 50 K Illumina Infinium iSelect SNP chip, genotyping by sequencing (GBS) and whole-genome sequencing (WGS), was used as a backbone for construction of two high-resolution maps. Using this approach, the target locus was narrowed down to 0.161 cM and 0.036 cM in the Igri × Chikurin Ibaraki 1 (I × C) and Chikurin Ibaraki 1 × Uschi (C × U) populations, respectively. Corresponding physical intervals of 11.3 Mbp and 0.281 Mbp were calculated for I × C and C × U, respectively, according to the Morex v3 genome sequence. In the 0.281 Mbp target region, six high confidence (HC) and two low confidence (LC) genes were identified. Genome assemblies of BaMMV-susceptible cultivars Igri and Golden Promise from the barley pan-genome, and a HiFi assembly of Chikurin Ibaraki 1 together with re-sequencing data for the six HC and two LC genes in susceptible parental cultivar Uschi revealed functional SNPs between resistant and susceptible genotypes only in two of the HC genes. These SNPs are the most promising candidates for the development of functional markers and the two genes represent promising candidates for functional analysis.

Comparative genomic analysis of Thermus provides insights into the evolutionary history of an incomplete denitrification pathway

Biological denitrification is a crucial process in the nitrogen biogeochemical cycle, and Thermus has been reported to be a significant heterotrophic denitrifier in terrestrial geothermal environments. However, neither the denitrification potential nor the evolutionary history of denitrification genes in the genus Thermus or phylum Deinococcota is well understood. Here, we performed a comparative analysis of 23 Thermus genomes and identified denitrification genes in 15 Thermus strains. We confirmed that Thermus harbors an incomplete denitrification pathway as none of the strains contain the nosZ gene. Ancestral character state reconstructions and phylogenetic analyses showed that narG, nirS, and norB genes were acquired by the last common ancestor of Thermales and were inherited vertically. In contrast, nirK of Thermales was acquired via two distinct horizontal gene transfers from Proteobacteria to the genus Caldithermus and from an unknown donor to the common ancestor of all known Thermus species except Thermus filiformis. This study expands our understanding of the genomic potential for incomplete denitrification in Thermus, revealing a largely vertical evolutionary history of the denitrification pathway in the Thermaceae, and supporting the important role for Thermus as an important heterotrophic denitrifier in geothermal environments.

SmGDB: genome database of Salvia miltiorrhiza, an important TCM Plant

BACKGROUND

Salvia miltiorrhiza is an important traditional Chinese medicinal (TCM) plant and a model plant in the genetic study of TCM. A series of omics related to Danshen have been published. Integrating, managing, storing, and sharing data has become an urgent problem to be solved in S. miltiorrhiza genetic studies.

OBJECTIVES

The genome database is the link for the exchange, acquisition, and use of different omics data between data producers and users, maximizing value and utilization of data.

METHODS

The genome database included DSS3 genome and five RNA-Seq data. The back-end performs data search and retrieval through the LAMP (Linux, Apache, MySQL, PHP) framework.

RESULTS

Here, we present SmGDB (S. miltiorrhiza genome database; http://8.140.162.85/ ), which houses the latest version of genome sequence and annotation data for S. miltiorrhiza, combining three unpublished RNA-Seq data from our group and two released RNA-Seq data. We also identified a novel gene cluster including seven CYP71D genes involved in the tanshinone synthesis pathway based on genome sequences and expression data. Besides, SmGDB provides user-friendly web interfaces for querying and browsing gene annotation, structure, location, and expression profiles for concerned genes. Popular bioinformatics tools such as 'BLAST', 'Search', 'Heatmap', 'JBrowse', etc., were also provided in SmGDB.

CONCLUSIONS

SmGDB will provide utility for characterizing the structure of the S. miltiorrhiza genome and better understanding gene functions and biological processes underlying complex secondary metabolism in Danshen.

Galbase: a comprehensive repository for integrating chicken multi-omics data

Background

Multi-omics data can provide a stereoscopic view to explore potential causal variations and genes, as well as underlying genetic mechanisms of complex traits. However, for many non-mammalian species, including chickens, these resources are poorly integrated and reused, greatly limiting genetic research and breeding processes of the species.

Results

Here, we constructed Galbase, an easily accessible repository that integrates public chicken multi-omics data from 928 re-sequenced genomes, 429 transcriptomes, 379 epigenomes, 15,275 QTL entries, and 7,526 associations. A total of 21.67 million SNPs, 2.71 million InDels, and 488,583 cis-regulatory elements were included. Galbase allows users to retrieve genomic variations in geographical maps, gene expression profiling in heatmaps, and epigenomic signals in peak patterns. It also provides modules for batch annotation of genes, regions, and loci based on multi-layered omics data. Additionally, a series of convenient tools, including the UCSC Genome Browser, WashU Epigenome Browser, BLAT, BLAST, and LiftOver, were also integrated to facilitate search, visualization, and analysis of sequence features.

Conclusion

Galbase grants new opportunities to research communities to undertake in-depth functional genomic studies on chicken. All features of Galbase make it a useful resource to identify genetic variations responsible for chicken complex traits. Galbase is publicly available at http://animal.nwsuaf.edu.cn/ChickenVar.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08598-2.

Quantitative Trait Loci Mapping of Adult Plant and Seedling Resistance to Stripe Rust (Puccinia striiformis Westend.) in a Multiparent Advanced Generation Intercross Wheat Population

Stripe rust caused by the biotrophic fungus Puccinia striiformis Westend. is one of the most important diseases of wheat worldwide, causing high yield and quality losses. Growing resistant cultivars is the most efficient way to control stripe rust, both economically and ecologically. Known resistance genes are already present in numerous cultivars worldwide. However, their effectiveness is limited to certain races within a rust population and the emergence of stripe rust races being virulent against common resistance genes forces the demand for new sources of resistance. Multiparent advanced generation intercross (MAGIC) populations have proven to be a powerful tool to carry out genetic studies on economically important traits. In this study, interval mapping was performed to map quantitative trait loci (QTL) for stripe rust resistance in the Bavarian MAGIC wheat population, comprising 394 F6 : 8 recombinant inbred lines (RILs). Phenotypic evaluation of the RILs was carried out for adult plant resistance in field trials at three locations across three years and for seedling resistance in a growth chamber. In total, 21 QTL for stripe rust resistance corresponding to 13 distinct chromosomal regions were detected, of which two may represent putatively new QTL located on wheat chromosomes 3D and 7D.

RGD v2.0: a major update of the ruminant functional and evolutionary genomics database

Ruminant Genome Database (RGD; http://animal.nwsuaf.edu.cn/RGD) provides visualization and analysis tools for ruminant comparative genomics and functional annotations. As more high-quality ruminant genome assemblies have become available, we have redesigned the user interface, integrated and expanded multi-omics data, and developed novel features to improve the database. The new version, RGD v2.0, houses 78 ruminant genomes; 110-species synteny alignments for major livestock (including cattle, sheep, goat) and wild ungulates; 21 012 orthologous gene clusters with Gene Ontology and pathway annotation; ∼8 600 000 conserved elements; and ∼1 000 000 cis-regulatory elements by utilizing 1053 epigenomic data sets. The transcriptome data in RGD v2.0 has nearly doubled, currently with 1936 RNA-seq data sets, and 155 174 phenotypic data sets have been newly added. New and updated features include: (i) The UCSC Genome Browser, BLAT, BLAST and Table Browser tools were updated for six available ruminant livestock species. (ii) The LiftOver tool was newly introduced into our browser to allow coordinate conversion between different ruminant assemblies. And (iii) tissue specificity index, tau, was calculated to facilitate batch screening of specifically expressed genes. The enhanced genome annotations and improved functionality in RGD v2.0 will be useful for study of genome evolution, environmental adaption, livestock breeding and biomedicine.

AprGPD: the apricot genomic and phenotypic database

BACKGROUND

Apricot is cultivated worldwide because of its high nutritive content and strong adaptability. Its flesh is delicious and has a unique and pleasant aroma. Apricot kernel is also consumed as nuts. The genome of apricot has been sequenced, and the transcriptome, resequencing, and phenotype data have been increasely generated. However, with the emergence of new information, the data are expected to integrate, and disseminate.

RESULTS

To better manage the continuous addition of new data and increase convenience, we constructed the apricot genomic and phenotypic database (AprGPD, http://apricotgpd.com ). At present, AprGPD contains three reference genomes, 1692 germplasms, 306 genome resequencing data, 90 RNA sequencing data. A set of user-friendly query, analysis, and visualization tools have been implemented in AprGPD. We have also performed a detailed analysis of 59 transcription factor families for the three genomes of apricot.

CONCLUSION

Six modules are displayed in AprGPD, including species, germplasm, genome, variation, product, tools. The data integrated by AprGPD will be helpful for the molecular breeding of apricot.

High-Throughput Sequencing of Small RNAs for the Sanitary Certification of Viruses in Grapevine

Biological indexing is the method generally recognized for the certification of propagative grapevines in many countries, and it is mandatory in the European Union. It consists of the evaluation of the plant material after grafting on indicators that are inspected for symptom development. This is a lengthy process that requires well-trained workers, testing field, etc. Alternative diagnostic methods such as serology and RT-qPCR have been discarded for certification because of their intrinsic drawbacks. In turn, high-throughput sequencing (HTS) of plant RNA has been proposed as a plausible alternative to bioassay, but before it is accepted, different aspects of this process must be evaluated. We have compared the HTS of small RNAs with bioassays and other diagnostic methods from a set of 40 grapevine plants submitted for certification. The results allowed the authors the identification of numerous grapevine viruses in the samples, as well as different variants. Besides, relationships between symptom expression and viromes were investigated, in particular leafroll-associated viruses. We compared HTS results using analytical and bioinformatics approaches in order to define minimum acceptable quality standards for certification schemes, resulting in a pipeline proposal. Finally, the comparison between HTS and bioassay resulted favorable for the former in terms of reliability, cost, and timing.

TCM-Blast for traditional Chinese medicine genome alignment with integrated resources

The traditional Chinese medicine (TCM) genome project aims to reveal the genetic information and regulatory network of herbal medicines, and to clarify their molecular mechanisms in the prevention and treatment of human diseases. Moreover, the TCM genome could provide the basis for the discovery of the functional genes of active ingredients in TCM, and for the breeding and improvement of TCM. The traditional Chinese Medicine Basic Local Alignment Search Tool (TCM-Blast) is a web interface for TCM protein and DNA sequence similarity searches. It contains approximately 40G of genome data on TCMs, including protein and DNA sequence for 36 TCMs with high medical value.The development of a publicly accessible TCM genome alignment database hosted on the TCM-Blast website ( http://viroblast.pungentdb.org.cn/TCM-Blast/viroblast.php ) has expanded to query multiple sequence databases to obtain TCM genome data, and provide user-friendly output for easy analysis and browsing of BLAST results. The genome sequencing of TCMs helps to elucidate the biosynthetic pathways of important secondary metabolites and provides an essential resource for gene discovery studies and molecular breeding. The TCMs genome provides a valuable resource for the investigation of novel bioactive compounds and drugs from these TCMs under the guidance of TCM clinical practice. Our database could be expanded to other TCMs after the determination of their genome data.

Identification and spatio-temporal expression analysis of barley genes that encode putative modular xylanolytic enzymes

Arabinoxylans are cell wall polysaccharides whose re-modelling and degradation during plant development are mediated by several classes of xylanolytic enzymes. Here, we present the identification and new annotation of twelve putative (1,4)-β-xylanase and six β-xylosidase genes, and their spatio-temporal expression patterns during vegetative and reproductive growth of barley (Hordeum vulgare cv. Navigator). The encoded xylanase proteins are all predicted to contain a conserved carbohydrate-binding module (CBM) and a catalytic glycoside hydrolase (GH) 10 domain. Additional domains in some xylanases define three discrete phylogenetic clades: one clade contains proteins with an additional N-terminal signal sequence, while another clade contains proteins with multiple CBMs. Homology modelling revealed that all fifteen xylanases likely contain a third domain, a β-sandwich folded from two non-contiguous sequence segments that bracket the catalytic GH domain, which may explain why the full length protein is required for correct folding of the active enzyme. Similarly, predicted xylosidase proteins share a highly conserved domain structure, each with an N-terminal signal peptide, a split GH 3 domain, and a C-terminal fibronectin-like domain. Several genes appear to be ubiquitously expressed during barley growth and development, while four newly annotated xylanase and xylosidase genes are expressed at extremely high levels, which may be of broader interest for industrial applications where cell wall degradation is necessary.

Fragrance in Pandanus amaryllifoliusRoxb. Despite the Presence of a Betaine Aldehyde Dehydrogenase 2

Pandanus amaryllifolius Roxb. accumulates the highest concentration of the major basmati aroma volatile 2-acetyl-1-pyrroline (2AP) in the plant kingdom. The expression of 2AP is correlated with the presence of a nonfunctional betaine aldehyde dehydrogenase 2(BADH2) in aromatic rice and other plant species. In the present study, a full-length BADH2 sequence was reconstructed from the transcriptome data of leaf tissue from P. amaryllifolius seedlings. Based on this sequence, a 1509 bp coding sequence was defined that encoded a 54 kD PaBADH2 protein. This revealed the presence of a full-length BADH2 protein in P. amaryllifolius. Moreover, quantitative real-time PCR analysis, combined with BADH2 enzyme activity, confirmed the expression and functionality of the PaBADH2 protein. To understand the apparent structural variation, docking analysis was carried out in which protein showed a good affinity with both betaine aldehyde (BAD) and γ-aminobutyraldehyde (GAB-ald) as substrates. Overall, the analysis showed the presence of a functional BADH2, along with substantial 2AP synthesis (4.38 ppm). Therefore, we conclude that unlike all other plants studied to date, 2AP biosynthesis in P. amaryllifolius is not due to the inactivation of BADH2.

Integrated Core Proteomics, Subtractive Proteomics, and Immunoinformatics Investigation to Unveil a Potential Multi-Epitope Vaccine against Schistosomiasis

Schistosomiasis is a parasitic infection that causes considerable morbidity and mortality in the world. Infections of parasitic blood flukes, known as schistosomes, cause the disease. No vaccine is available yet and thus there is a need to design an effective vaccine against schistosomiasis. Schistosoma japonicum, Schistosoma mansoni, and Schistosoma haematobium are the main pathogenic species that infect humans. In this research, core proteomics was combined with a subtractive proteomics pipeline to identify suitable antigenic proteins for the construction of a multi-epitope vaccine (MEV) against human-infecting Schistosoma species. The pipeline revealed two antigenic proteins-calcium binding and mycosubtilin synthase subunit C-as promising vaccine targets. T and B cell epitopes from the targeted proteins were predicted using multiple bioinformatics and immunoinformatics databases. Seven cytotoxic T cell lymphocytes (CTL), three helper T cell lymphocytes (HTL), and four linear B cell lymphocytes (LBL) epitopes were fused with a suitable adjuvant and linkers to design a 217 amino-acid-long MEV. The vaccine was coupled with a TLR-4 agonist (RS-09; Sequence: APPHALS) adjuvant to enhance the immune responses. The designed MEV was stable, highly antigenic, and non-allergenic to human use. Molecular docking, molecular dynamics (MD) simulations, and molecular mechanics/generalized Born surface area (MMGBSA) analysis were performed to study the binding affinity and molecular interactions of the MEV with human immune receptors (TLR2 and TLR4) and MHC molecules (MHC I and MHC II). The MEV expression capability was tested in an Escherichia coli (strain-K12) plasmid vector pET-28a(+). Findings of these computer assays proved the MEV as highly promising in establishing protective immunity against the pathogens; nevertheless, additional validation by in vivo and in vitro experiments is required to discuss its real immune-protective efficacy.

PID: An integrative and comprehensive platform of plant intron

Intron is a non-coding sequence of a broken gene and participates in important biological processes, such as transcription regulation, alternative splicing, and nuclear export. With the development of plant genomes, a comprehensive platform for intron analysis in plants must be established. Plant Intron Database (PID), a publicly available searchable database, was developed to efficiently store, query, analyze, and integrate intron resources in plants. The information of intron, exon, and gene can be searched by key words in PID. Users cannot only view intron length distribution pie chart and 5' and 3' splice site sequence feature maps in a statistical interface but can also browse intron information in a graphical visualization interface through JBrowse. ViroBlast for sequence homology searches, Intron detection and sequence interception tools were also provided. PID contains annotated genes from 118 sequenced plants, 24,782,048 introns, 30,843,049 exons, and 414 visual maps. This tool will greatly accelerate research on the distribution, length characteristics, and functions of introns in plants. PID is accessible at http://biodb.sdau.edu.cn/PID/index.php.