PmiREN2.0: from data annotation to functional exploration of plant microRNAs

Widespread changes to the translational landscape in a maize microRNA biogenesis mutant

MicroRNAs are short, non-coding RNAs that repress gene expression in both plants and animals and have diverse functions related to growth, development, and stress responses. The ribonuclease, DICER-LIKE1 (DCL1) is required for two steps in plant miRNA biogenesis: cleavage of the primary miRNAs (pri-miRNAs) to release a hairpin structure, called the precursor miRNA (pre-miRNA) and cleavage of the pre-miRNA to generate the miRNA/miRNA* duplex. The mature miRNA guides the RNA-induced silencing complex to target RNAs with complementary sequences, resulting in translational repression and/or RNA cleavage of target mRNAs. However, the relative contribution of translational repression versus mRNA degradation by miRNAs remains unknown at the genome-level in crops, especially in maize. The maize fuzzy tassel (fzt) mutant contains a hypomorphic mutation in DCL1 resulting in broad developmental defects. While most miRNAs are reduced in fzt, the levels of miRNA-targeted mRNAs are not dramatically increased, suggesting that translational regulation by miRNAs may be common. To gain insight into the repression mechanism of plant miRNAs, we combined ribosome profiling and RNA-sequencing to globally survey miRNA activities in maize. Our data indicate that translational repression contributes significantly to regulation of most miRNA targets and that approximately one-third of miRNA targets are regulated primarily at the translational level. Surprisingly, ribosomes appear altered in fzt mutants suggesting that DCL1 may also have a role in ribosome biogenesis. Thus, DICER-LIKE1 shapes the translational landscape in plants through both miRNA-dependent and miRNA-independent mechanisms.

The de novo, chromosome-level genome assembly of the sweet chestnut (Castanea sativa Mill.) Cv. Marrone Di Chiusa Pesio

OBJECTIVES

The sweet chestnut Castanea sativa Mill. is the only native Castanea species in Europe, and it is a tree of high economic value that provides appreciated fruits and valuable wood. In this study, we assembled a high-quality nuclear genome of the ancient Italian chestnut variety 'Marrone di Chiusa Pesio' using a combination of Oxford Nanopore Technologies long reads, whole-genome and Omni-C Illumina short reads.

DATA DESCRIPTION

The genome was assembled into 238 scaffolds with an N50 size of 21.8 Mb and an N80 size of 7.1 Mb for a total assembled sequence of 750 Mb. The BUSCO assessment revealed that 98.6% of the genome matched the embryophyte dataset, highlighting good completeness of the genetic space. After chromosome-level scaffolding, 12 chromosomes with a total length of 715.8 and 713.0 Mb were constructed for haplotype 1 and haplotype 2, respectively. The repetitive elements represented 37.3% and 37.4% of the total assembled genome in haplotype 1 and haplotype 2, respectively. A total of 57,653 and 58,146 genes were predicted in the two haplotypes, and approximately 73% of the genes were functionally annotated using the EggNOG-mapper. The assembled genome will be a valuable resource and reference for future chestnut breeding and genetic improvement.

Species-Specific miRNAs Contribute to the Divergence between Deciduous and Evergreen Species in Ilex

MicroRNAs (miRNAs) are pivotal regulators of gene expression, playing crucial roles in plant developmental processes and environmental responses. However, the function of miRNAs in influencing deciduous traits has been little explored. Here, we utilized sRNA-seq on two deciduous species, Ilex polyneura (Hand.-Mazz.) S. Y. Hu and Ilex asprella Champ. ex Benth., along with an evergreen species, Ilex latifolia Thunb., to identify and annotate miRNAs within these species. Our analysis revealed 162 species-specific miRNAs (termed SS-miRNAs) from 120 families, underscoring the fundamental roles and potential influence of SS-miRNAs on plant phenotypic diversity and adaptation. Notably, three SS-miRNAs in I. latifolia were found to target crucial genes within the abscission signaling pathway. Analysis of cis-regulatory elements suggested a novel regulatory relationship that may contribute to the evergreen phenotype of I. latifolia by modulating the abscission process in a light-independent manner. These findings propose a potential mechanism by which SS-miRNAs can influence the conserved abscission pathway, contributing to the phenotypic divergence between deciduous and evergreen species within the genus Ilex.

Improving plant miRNA-target prediction with self-supervised k-mer embedding and spectral graph convolutional neural network

Deciphering the targets of microRNAs (miRNAs) in plants is crucial for comprehending their function and the variation in phenotype that they cause. As the highly cell-specific nature of miRNA regulation, recent computational approaches usually utilize expression data to identify the most physiologically relevant targets. Although these methods are effective, they typically require a large sample size and high-depth sequencing to detect potential miRNA-target pairs, thereby limiting their applicability in improving plant breeding. In this study, we propose a novel miRNA-target prediction framework named kmerPMTF (k-mer-based prediction framework for plant miRNA-target). Our framework effectively extracts the latent semantic embeddings of sequences by utilizing k-mer splitting and a deep self-supervised neural network. We construct multiple similarity networks based on k-mer embeddings and employ graph convolutional networks to derive deep representations of miRNAs and targets and calculate the probabilities of potential associations. We evaluated the performance of kmerPMTF on four typical plant datasets: Arabidopsis thaliana, Oryza sativa, Solanum lycopersicum, and Prunus persica. The results demonstrate its ability to achieve AUPRC values of 84.9%, 91.0%, 80.1%, and 82.1% in 5-fold cross-validation, respectively. Compared with several state-of-the-art existing methods, our framework achieves better performance on threshold-independent evaluation metrics. Overall, our study provides an efficient and simplified methodology for identifying plant miRNA-target associations, which will contribute to a deeper comprehension of miRNA regulatory mechanisms in plants.

Identification of miRNA858 long-loop precursors in seed plants

MicroRNAs (miRNAs) are a class of nonprotein-coding short transcripts that provide a layer of post-transcriptional regulation essential to many plant biological processes. MiR858, which targets the transcripts of MYB transcription factors, can affect a range of secondary metabolic processes. Although miR858 and its 187-nt precursor have been well studied in Arabidopsis (Arabidopsis thaliana), a systematic investigation of miR858 precursors and their functions across plant species is lacking due to a problem in identifying the transcripts that generate this subclass. By re-evaluating the transcript of miR858 and relaxing the length cut-off for identifying hairpins, we found in kiwifruit (Actinidia chinensis) that miR858 has long-loop hairpins (1,100 to 2,100 nt), whose intervening sequences between miRNA generating complementary sites were longer than all previously reported miRNA hairpins. Importantly, these precursors of miR858 containing long-loop hairpins (termed MIR858L) are widespread in seed plants including Arabidopsis, varying between 350 and 5,500 nt. Moreover, we showed that MIR858L has a greater impact on proanthocyanidin and flavonol levels in both Arabidopsis and kiwifruit. We suggest that an active MIR858L-MYB regulatory module appeared in the transition of early land plants to large upright flowering plants, making a key contribution to plant secondary metabolism.

Adaptive regulation of miRNAs/milRNAs in tissue-specific interaction between apple and Valsa mali

In plant-pathogen interactions, pathogens display tissue specificity, infecting and causing disease in particular tissues. However, the involvement of microRNAs/microRNA-like RNAs (miRNAs/milRNAs) in tissue-specific regulation during plant-pathogen interactions remains largely unexplored. This study investigates the differential expression of miRNAs/milRNAs, as well as their corresponding target genes, in interactions between Valsa mali (Vm) and different apple tissues. The results demonstrated that both apple miRNAs and Vm milRNAs exhibited distinct expression profiles when Vm infected bark and leaves, with functionally diverse corresponding target genes. Furthermore, one apple miRNA (Mdo-miR482a) and one Vm milRNA (Vm-milR57) were identified as exhibiting tissue-specific expression in interactions between Vm and apple bark or leaves. Mdo-miR482a was exclusively up-regulated in response to Vm infection in bark and target a nucleotide-binding leucine-rich repeat (NLR) gene of apple. When Mdo-miR482a was transiently over-expressed or silenced, the resistance was significantly reduced or improved. Similarly, transient expression of the NLR gene also showed an increase in resistance. Vm-milR57 could target two essential pathogenicity-related genes of Vm. During Vm infection in bark, the expression of Vm-milR57 was down-regulated to enhance the expression of the corresponding target gene to improve the pathogenicity. The study is the first to reveal tissue-specific characteristics of apple miRNAs and Vm milRNAs in interactions between Vm and different apple tissues, providing new insights into adaptive regulation in tissue-specific interactions between plants and fungi.

Genome-wide view and characterization of natural antisense transcripts in Cannabis Sativa L

Natural Antisense Transcripts (NATs) are a kind of complex regulatory RNAs that play crucial roles in gene expression and regulation. However, the NATs in Cannabis Sativa L., a widely economic and medicinal plant rich in cannabinoids remain unknown. In this study, we comprehensively predicted C. sativa NATs genome-wide using strand-specific RNA sequencing (ssRNA-Seq) data, and validated the expression profiles by strand-specific quantitative reverse transcription PCR (ssRT-qPCR). Consequently, a total of 307 NATs were predicted in C. sativa, including 104 cis- and 203 trans- NATs. Functional enrichment analysis demonstrated the potential involvement of the C. sativa NATs in DNA polymerase activity, RNA-DNA hybrid ribonuclease activity, and nucleic acid binding. Finally, 18 cis- and 376 trans- NAT-ST pairs were predicted to produce 621 cis- and 5,679 trans- small interfering RNA (nat-siRNAs), respectively. These nat-siRNAs were potentially involved in the biosynthesis of cannabinoids and cellulose. All these results will shed light on the regulation of NATs and nat-siRNAs in C. sativa.

The grapevine miR827a regulates the synthesis of stilbenes by targeting VqMYB14 and gives rise to susceptibility in plant immunity

Grapevine (Vitis vinifera L.) is an economically important fruit crop cultivated worldwide. In China, grapevine cultivation is very extensive, and a few Vitis grapes have excellent pathogen and stress resistance, but the molecular mechanisms underlying the grapevine response to stress remain unclear. In this study, a microRNA (miRNA; miR827a), which negatively regulates its target gene VqMYB14, a key regulatory role in the synthesis of stilbenes, was identified in Vitis quinquangularis (V. quinquangularis) using transcriptome sequencing. Using overexpression and silencing approaches, we found that miR827a regulates the synthesis of stilbenes by targeting VqMYB14. We used flagellin N-terminal 22-amino-acid peptide (flg22), the representative elicitor in plant basal immunity, as the elicitor to verify whether miR827a is involved in the basal immunity of V. quinquangularis. Furthermore, the promoter activity of miR827a was alleviated in transgenic grape protoplasts and Arabidopsis thaliana following treatment with flg22 and Pseudomonas syringae pv. Tomato DC3000 (Pst DC3000), respectively. In addition, yeast one-hybrid and dual luciferase reporter assay revealed that the ethylene transcription factor VqERF057 acted as a key regulator in the inhibition of miR827a transcription. These results will contribute to the understanding of the biological functions of miR827a in grapevine and clarify the molecular mechanism of the interaction between miR827a and VqMYB14.

sRNAminer: A multifunctional toolkit for next-generation sequencing small RNA data mining in plants

Small RNAs (sRNAs), found extensively in plants, play an essential role in plant growth and development. Although various sRNA analysis tools have been developed for plants, the use of most of them depends on programming and command-line environments, which is a challenge for many wet-lab biologists. Furthermore, current sRNA analysis tools mostly focus on the analysis of certain type of sRNAs and are resource-intensive, normally demanding an immense amount of time and effort to learn the use of numerous tools or scripts and assemble them into a workable pipeline to get the final results. Here, we present sRNAminer, a powerful stand-alone toolkit with a user-friendly interface that integrates all common functions for the analysis of three major types of plant sRNAs: microRNAs (miRNAs), phased small interfering RNAs (phasiRNAs), and heterochromatic siRNAs (hc-siRNAs). We constructed a curated or "golden" set of MIRNA and PHAS loci, which was used to assess the performance of sRNAminer in comparison to other existing tools. The results showed that sRNAminer outperformed these tools in multiple aspects, highlighting its functionality. In addition, to enable an efficient evaluation of sRNA annotation results, we developed Integrative Genomics Viewer (IGV)-sRNA, a modified genome browser optimized from IGV and we incorporated it as a functional module in sRNAminer. IGV-sRNA can display a wealth of sRNA-specific features, enabling a more comprehensive understanding of sRNA data. sRNAminer and IGV-sRNA are both platform-independent software that can be run under all operating systems. They are now freely available at https://github.com/kli28/sRNAminer and https://gitee.com/CJchen/IGV-sRNA.

Genome-wide identification and in silico characterization of major RNAi gene families in date palm (Phoenix dactylifera)

BACKGROUND

Dates contain various minerals that are essential for good health. The major RNA interference (RNAi) gene families play a vital role in plant growth and development by controlling the expression of protein-coding genes against different biotic and abiotic stresses. However, these gene families for date palm are not yet studied. Therefore, this study has explored major RNAi genes and their characteristics in date palm.

RESULTS

We have identified 4 PdDCLs, 7 PdAGOs, and 3 PdRDRs as RNAi proteins from the date palm genome by using AtRNAi genes as query sequences in BLASTp search. Domain analysis of predicted RNAi genes has revealed the Helicase_C, Dicer_dimer, PAZ, RNase III, and Piwi domains that are associated with the gene silencing mechanisms. Most PdRNAi proteins have been found in the nucleus and cytosol associated with the gene silencing actions. The gene ontology (GO) enrichment analysis has revealed some important GO terms including RNA interference, dsRNA fragmentation, and ribonuclease_III activity that are related to the protein-coding gene silencing mechanisms. Gene regulatory network (GRN) analysis has identified PAZ and SNF2 as the transcriptional regulators of PdRNAi genes. Top-ranked 10 microRNAs including Pda-miR156b, Pda-miR396a, Pda-miR166a, Pda-miR167d, and Pda-miR529a have been identified as the key post-transcriptional regulators of PdRNAi genes that are associated with different biotic/abiotic stresses. The cis-acting regulatory element analysis of PdRNAi genes has detected some vital cis-acting elements including ABRE, MBS, MYB, MYC, Box-4, G-box, I-box, and STRE that are linked with different abiotic stresses.

CONCLUSION

The results of this study might be valuable resources for the improvement of different characteristics in date palm by further studies in wet-lab.

Genome-wide analysis of plant specific YABBY transcription factor gene family in carrot (Dacus carota) and its comparison with Arabidopsis

YABBY gene family is a plant-specific transcription factor with DNA binding domain involved in various functions i.e. regulation of style, length of flowers, and polarity development of lateral organs in flowering plants. Computational methods were utilized to identify members of the YABBY gene family, with Carrot (Daucus carota) 's genome as a foundational reference. The structure of genes, location of the chromosomes, protein motifs and phylogenetic investigation, syntony and transcriptomic analysis, and miRNA targets were analyzed to unmask the hidden structural and functional characteristics YABBY gene family in Carrots. In the following research, it has been concluded that 11 specific YABBY genes irregularly dispersed on all 9 chromosomes and proteins assembled into five subgroups i.e. AtINO, AtCRC, AtYAB5, AtAFO, and AtYAB2, which were created on the well-known classification of Arabidopsis. The wide ranges of YABBY genes in carrots were dispersed due to segmental duplication, which was detected as prevalent when equated to tandem duplication. Transcriptomic analysis showed that one of the DcYABBY genes was highly expressed during anthocyanin pigmentation in carrot taproots. The cis-regulatory elements (CREs) analysis unveiled elements that particularly respond to light, cell cycle regulation, drought induce ability, ABA hormone, seed, and meristem expression. Furthermore, a relative study among Carrot and Arabidopsis genes of the YABBY family indicated 5 sub-families sharing common characteristics. The comprehensive evaluation of YABBY genes in the genome provides a direction for the cloning and understanding of their functional properties in carrots. Our investigations revealed genome-wide distribution and role of YABBY genes in the carrots with best-fit comparison to Arabidopsis thaliana.

Molecular mechanism of miRNA mediated biosynthesis of secondary metabolites in medicinal plants

Plant secondary metabolites are important raw materials for the pharmaceutical industry, and their biosynthetic processes are subject to diverse and precise regulation by miRNA. The identification of miRNA molecules in medicinal plants and exploration of their mechanisms not only contribute to a deeper understanding of the molecular genetic mechanisms of plant growth, development and resistance to stress, but also provide a theoretical basis for elucidating the pharmacological effects of authentic medicinal materials and constructing bioreactors for the synthesis of medicinal secondary metabolite components. This paper summarizes the research reports on the discovery of miRNA in medicinal plants and their regulatory mechanisms on the synthesis of secondary metabolites by searching the relevant literature in public databases. It summarizes the currently discovered miRNA and their functions in medicinal plants, and summarizes the molecular mechanisms regulating the synthesis and degradation of secondary metabolites. Furthermore, it provides a prospect for the research and development of medicinal plant miRNA. The compiled information contributes to a comprehensive understanding of the research progress on miRNA in medicinal plants and provides a reference for the industrial development of related secondary metabolite biosynthesis.

Plant non-coding RNAs: The new frontier for the regulation of plant development and adaptation to stress

Most plant transcriptomes constitute functional non-coding RNAs (ncRNAs) that lack the ability to encode proteins. In recent years, more research has demonstrated that ncRNAs play important regulatory roles in almost all plant biological processes by modulating gene expression. Thus, it is important to study the biogenesis and function of ncRNAs, particularly in plant growth and development and stress tolerance. In this review, we systematically explore the process of formation and regulatory mechanisms of ncRNAs, particularly those of microRNAs (miRNAs), small interfering RNAs (siRNAs), long non-coding RNAs (lncRNAs), and circular RNAs (circRNAs). Additionally, we provide a comprehensive overview of the recent advancements in ncRNAs research, including their regulation of plant growth and development (seed germination, root growth, leaf morphogenesis, floral development, and fruit and seed development) and responses to abiotic and biotic stress (drought, heat, cold, salinity, pathogens and insects). We also discuss research challenges and provide recommendations to advance the understanding of the roles of ncRNAs in agronomic applications.

HetFCM: functional co-module discovery by heterogeneous network co-clustering

Functional molecular module (i.e., gene-miRNA co-modules and gene-miRNA-lncRNA triple-layer modules) analysis can dissect complex regulations underlying etiology or phenotypes. However, current module detection methods lack an appropriate usage and effective model of multi-omics data and cross-layer regulations of heterogeneous molecules, causing the loss of critical genetic information and corrupting the detection performance. In this study, we propose a heterogeneous network co-clustering framework (HetFCM) to detect functional co-modules. HetFCM introduces an attributed heterogeneous network to jointly model interplays and multi-type attributes of different molecules, and applies multiple variational graph autoencoders on the network to generate cross-layer association matrices, then it performs adaptive weighted co-clustering on association matrices and attribute data to identify co-modules of heterogeneous molecules. Empirical study on Human and Maize datasets reveals that HetFCM can find out co-modules characterized with denser topology and more significant functions, which are associated with human breast cancer (subtypes) and maize phenotypes (i.e., lipid storage, drought tolerance and oil content). HetFCM is a useful tool to detect co-modules and can be applied to multi-layer functional modules, yielding novel insights for analyzing molecular mechanisms. We also developed a user-friendly module detection and analysis tool and shared it at http://www.sdu-idea.cn/FMDTool.

Small RNA profiling reveals that an ovule-specific microRNA, cja-miR5179, targets a B-class MADS-box gene in Camellia japonica

BACKGROUND AND AIMS

The functional specialization of microRNA and its target genes is often an important factor in the establishment of spatiotemporal patterns of gene expression that are essential to plant development and growth. In different plant lineages, understanding the functional conservation and divergence of microRNAs remains to be explored.

METHODS

To identify small regulatory RNAs underlying floral patterning, we performed a tissue-specific profiling of small RNAs in various floral organs from single and double flower varieties (flowers characterized by multiple layers of petals) in Camellia japonica. We identified cja-miR5179, which belongs to a deeply conserved microRNA family that is conserved between angiosperms and basal plants but frequently lost in eudicots. We characterized the molecular function of cja-miR5179 and its target - a B-function MADS-box gene - through gene expression analysis and transient expression assays.

KEY RESULTS

We showed that cja-miR5179 is exclusively expressed in ovule tissues at the early stage of floral development. We found that cja-miR5179 targets the coding sequences of a DEFICIENS-like B-class gene (CjDEF) mRNA, which is located in the K motif of the MADS-box domain; and the target sites of miR5179/MADS-box were consistent in Camellia and orchids. Furthermore, through a petal transient-expression assay, we showed that the BASIC PENTACYSTEINE proteins bind to the GA-rich motifs in the cja-miR5179 promoter region and suppresses its expression.

CONCLUSIONS

We propose that the regulation between miR5179 and a B-class MADS-box gene in C. japonica has a deep evolutionary origin before the separation of monocots and dicots. During floral development of C. japonica, cja-miR5179 is specifically expressed in the ovule, which may be required for the inhibition of CjDEF function. This work highlights the evolutionary conservation as well as functional divergence of small RNAs in floral development.

Genetic Databases and Gene Editing Tools for Enhancing Crop Resistance against Abiotic Stress

Abiotic stresses extensively reduce agricultural crop production globally. Traditional breeding technology has been the fundamental approach used to cope with abiotic stresses. The development of gene editing technology for modifying genes responsible for the stresses and the related genetic networks has established the foundation for sustainable agriculture against environmental stress. Integrated approaches based on functional genomics and transcriptomics are now expanding the opportunities to elucidate the molecular mechanisms underlying abiotic stress responses. This review summarizes some of the features and weblinks of plant genome databases related to abiotic stress genes utilized for improving crops. The gene-editing tool based on clustered, regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9) has revolutionized stress tolerance research due to its simplicity, versatility, adaptability, flexibility, and broader applications. However, off-target and low cleavage efficiency hinder the successful application of CRISPR/Cas systems. Computational tools have been developed for designing highly competent gRNA with better cleavage efficiency. This powerful genome editing tool offers tremendous crop improvement opportunities, overcoming conventional breeding techniques' shortcomings. Furthermore, we also discuss the mechanistic insights of the CRISPR/Cas9-based genome editing technology. This review focused on the current advances in understanding plant species' abiotic stress response mechanism and applying the CRISPR/Cas system genome editing technology to develop crop resilience against drought, salinity, temperature, heavy metals, and herbicides.

Vv-circSIZ1 mediated by pre-mRNA processing machinery contributes to salt tolerance

CircRNAs exist widely in plants, but the regulatory mechanisms for the biogenesis and function of plant circRNAs remain largely unknown. Using extensive mutagenesis of expression plasmids and genetic transformation methods, we analyzed the biogenesis and anti-salt functions of a new grape circRNA Vv-circSIZ1. We identified Vv-circSIZ1 that is mainly expressed in the cytoplasm of xylem. CircSIZ1 is species-specific, and genomic circSIZ1-forming region of seven tested species could be backspliced in Nicotiana benthamiana, but not in Arabidopsis. The retention length of Vv-circSIZ1 flanking introns was significantly positively correlated with its generation efficiency. The precise splicing of Vv-circSIZ1 does not depend on its mature exon sequence or internal intron sequences, but on the AG/GT splicing signal sites and branch site of the flanking introns. The spliceosome activity was inversely proportional to the expression level of Vv-circSIZ1. Furthermore, RNA-binding proteins can regulate the expression of Vv-circSIZ1. The overexpression of Vv-circSIZ1 improved salt tolerance of grape and N. benthamiana. Additionally, Vv-circSIZ1 could relieve the repressive effect of VvmiR3631 on its target VvVHAc1. Vv-circSIZ1 also promoted transcription of its parental gene. Overall, these results broaden our understanding of circRNAs in plants.

Characterization of MicroRNAs and Gene Expression in ACC Oxidase RNA Interference-Based Transgenic Bananas

Banana (Musa acuminata, AAA group) is a typically respiratory climacteric fruit. Previously, genes encoding ACC oxidase, one of the key enzymes in ethylene biosynthesis, Mh-ACO1 and Mh-ACO2 in bananas were silenced individually using RNAi interference technology, and fruit ripening of transgenic bananas was postponed. Here, the differential expression of miRNAs and their targeted mRNAs were analyzed in the transcriptomes of fruits at the third ripening stage, peel color more green than yellow, from the untransformed and RNAi transgenic bananas. Five significantly differentially expressed miRNAs (mac-miR169a, mac-miR319c-3p, mac-miR171a, mac-miR156e-5p, and mac-miR164a-5p) were identified. The predicted miRNA target genes were mainly enriched in six KEGG pathways, including 'sulfur relay system', 'protein digestion and absorption', 'histidine metabolism', 'pathogenic E. coli infection', 'sulfur metabolism', and 'starch and sucrose metabolism'. After ethylene treatment, the expression of ACC oxidase silencing-associated miRNAs was down-regulated, and that of their target genes was up-regulated along with fruit ripening. The evolutionary clustering relationships of miRNA precursors among 12 gene families related to fruit ripening were analyzed. The corresponding expression patterns of mature bodies were mainly concentrated in flowers, fruits, and leaves. Our results indicated that ethylene biosynthesis is associated with miRNAs regulating the expression of sulfur metabolism-related genes in bananas.

Recent genome resequencing paraded COBRA-Like gene family roles in abiotic stress and wood formation in Poplar

A cell wall determines the mechanical properties of a cell, serves as a barrier against plant stresses, and allows cell division and growth processes. The COBRA-Like (COBL) gene family encodes a putative glycosylphosphatidylinositol (GPI)-anchored protein that controls cellulose deposition and cell progression in plants by contributing to the microfibril orientation of a cell wall. Despite being studied in different plant species, there is a dearth of the comprehensive global analysis of COBL genes in poplar. Poplar is employed as a model woody plant to study abiotic stresses and biomass production in tree research. Improved genome resequencing has enabled the comprehensive exploration of the evolution and functional capacities of PtrCOBLs (Poplar COBRA-Like genes) in poplar. Phylogeny analysis has discerned and classified PtrCOBLs into two groups resembling the Arabidopsis COBL family, and group I genes possess longer proteins but have fewer exons than group II. Analysis of gene structure and motifs revealed PtrCOBLs maintained a rather stable motif and exon-intron pattern across members of the same group. Synteny and collinearity analyses exhibited that the evolution of the COBL gene family was heavily influenced by gene duplication events. PtrCOBL genes have undergone both segmental duplication and tandem duplication, followed by purifying selection. Promotor analysis flaunted various phytohormone-, growth- and stress-related cis-elements (e.g., MYB, ABA, MeJA, SA, AuxR, and ATBP1). Likewise, 29 Ptr-miRNAs of 20 families were found targeting 11 PtrCOBL genes. PtrCOBLs were found localized at the plasma membrane and extracellular matrix, while gene ontology analysis showed their involvement in plant development, plant growth, stress response, cellulose biosynthesis, and cell wall biogenesis. RNA-seq datasets depicted the bulk of PtrCOBL genes expression being found in plant stem tissues and leaves, rendering mechanical strength and rejoinders to environmental cues. PtrCOBL2, 3, 10, and 11 manifested the highest expression in vasculature and abiotic stress, and resemblant expression trends were upheld by qRT-PCR. Co-expression network analysis identified PtrCOBL2 and PtrCOBL3 as hub genes across all abiotic stresses and wood developing tissues. The current study reports regulating roles of PtrCOBLs in xylem differentiating tissues, tension wood formation, and abiotic stress latency that lay the groundwork for future functional studies of the PtrCOBL genes in poplar breeding.

Characterization of HSP70 family in watermelon (Citrullus lanatus): identification, structure, evolution, and potential function in response to ABA, cold and drought stress

Watermelon (Citrullus lanatus) as a crop with important economic value, is widely cultivated around the world. The heat shock protein 70 (HSP70) family in plant is indispensable under stress conditions. However, no comprehensive analysis of watermelon HSP70 family is reported to date. In this study, 12 ClHSP70 genes were identified from watermelon, which were unevenly located in 7 out of 11 chromosomes and divided into three subfamilies. ClHSP70 proteins were predicted to be localized primarily in cytoplasm, chloroplast, and endoplasmic reticulum. Two pairs of segmental repeats and 1 pair of tandem repeats existed in ClHSP70 genes, and ClHSP70s underwent strong purification selection. There were many abscisic acid (ABA) and abiotic stress response elements in ClHSP70 promoters. Additionally, the transcriptional levels of ClHSP70s in roots, stems, true leaves, and cotyledons were also analyzed. Some of ClHSP70 genes were also strongly induced by ABA. Furthermore, ClHSP70s also had different degrees of response to drought and cold stress. The above data indicate that ClHSP70s may be participated in growth and development, signal transduction and abiotic stress response, laying a foundation for further analysis of the function of ClHSP70s in biological processes.

Delineating the tissue-mediated drought stress governed tuning of conserved miR408 and its targets in rice

Engineering drought tolerance in rice needs to focus on regulators that enhance tolerance while boosting plant growth and vigor. The present study delineated the concealed function and tissue-mediated interplay of the miR408/target module in imparting drought stress tolerance in rice. The plant miR408 family comprises three dominant mature forms (21 nt), including a distinct monocot variant (F-7 with 5' C) and is divided into six groups. miR408 majorly cleaves genes belonging to the blue copper protein in addition to several other species-specific targets in plants. Comparative sequence analysis in 4726 rice accessions identified 22 sequence variants (SNP and InDELs) in its promoter (15) and pre-miR408 region. Haplotype analysis of the sequence variants indicated eight haplotypes (three: Japonica-specific and five: Indica-specific) of the miR408 promoter. In drought-tolerant Nagina 22, miR408 follows flag leaf preferential expression. Under drought conditions, its levels are upregulated in flag leaf and roots which seems to be regulated by a differential fraction of methylated cytosines (mCs) in the precursor region. The active pool of miR408 regulated targets under control and drought conditions is impacted by the tissue type. Comparative expression analysis of the miR408/target module under different sets of conditions features 83 targets exhibiting antagonistic expression in rice, out of which 12 genes, including four PLANTACYANINS (OsUCL6, 7, 9 and 30), PIRIN, OsLPR1, OsCHUP1, OsDOF12, OsBGLU1, glycine-rich cell wall gene, OsDUT, and OsERF7, are among the high confidence targets. Further, overexpression of MIR408 in drought-sensitive rice cultivar (PB1) leads to the massive enhancement of vegetative growth in rice with improved ETR and Y(II) and enhanced dehydration stress tolerance. The above results suggest that miR408 is likely to act as a positive regulator of growth and vigor, as well as dehydration stress, making it a potential candidate for engineering drought tolerance in rice.

Comparative Transcriptome Analysis Reveals the Effect of miR156a Overexpression on Mineral Nutrient Homeostasis in Nicotiana tabacum

Mineral nutrition plays an important role in crop growth, yield and quality. MiR156 is a regulatory hub for growth and development. To date, the understanding of miR156-mediated mineral homeostasis is limited. In this study, we overexpressed Nta-miR156a in the tobacco cultivar TN90 and analyzed the effects of miR156 on mineral element homeostasis in tobacco by comparative transcriptome analysis. The results showed that the overexpression of miR156a caused significant morphological changes in transgenic tobacco. Chlorophyll and three anti-resistance markers, proline, total phenolics, and total flavonoids, were altered due to increased miR156 expression levels. Interestingly, the distribution of Cu, Mn, Zn, and Fe in different tissues of transgenic tobacco was disordered compared with that of the wild type. Comparative transcriptome analysis showed that the overexpression of miR156 resulted in 2656 significantly differentially expressed genes. The expression levels of several metal-transport-related genes, such as NtABC, NtZIP, NtHMA, and NtCAX, were significantly increased or decreased in transgenic tobacco. These results suggest that miR156 plays an essential role in regulating mineral homeostasis. Our study provides a new perspective for the further study of mineral nutrient homeostasis in plants.

The NtSPL Gene Family in Nicotiana tabacum: Genome-Wide Investigation and Expression Analysis in Response to Cadmium Stress

The SQUAMOSA promoter binding protein-like (SPL)SPL family genes play an important role in regulating plant growth and development, synthesis of secondary metabolites, and resistance to stress. Understanding of the role of the SPL family in tobacco is still limited. In this study, 42 NtSPL genes were identified from the genome of the tobacco variety TN90. According to the results of the conserved motif and phylogenetic tree, the NtSPL genes were divided into eight subgroups, and the genes in the same subgroup showed similar gene structures and conserved domains. The cis-acting element analysis of the NtSPL promoters showed that the NtSPL genes were regulated by plant hormones and stresses. Twenty-eight of the 42 NtSPL genes can be targeted by miR156. Transcriptome data and qPCR results indicated that the expression pattern of miR156-targeted NtSPL genes was usually tissue specific. The expression level of miR156 in tobacco was induced by Cd stress, and the expression pattern of NtSPL4a showed a significant negative correlation with that of miR156. These results suggest that miR156-NtSPL4a may mediate the tobacco response to Cd stress. This study lays a foundation for further research on the function of the NtSPL gene and provides new insights into the involvement of NtSPL genes in the plant response to heavy metal stress.

LettuceGDB: The community database for lettuce genetics and omics

As a globally popular leafy vegetable and a representative plant of the Asteraceae family, lettuce has great economic and academic significance. In the last decade, high-throughput sequencing, phenotyping, and other multi-omics data in lettuce have accumulated on a large scale, thus increasing the demand for an integrative lettuce database. Here, we report the establishment of a comprehensive lettuce database, LettuceGDB (https://www.lettucegdb.com/). As an omics data hub, the current LettuceGDB includes two reference genomes with detailed annotations; re-sequencing data from over 1000 lettuce varieties; a collection of more than 1300 worldwide germplasms and millions of accompanying phenotypic records obtained with manual and cutting-edge phenomics technologies; re-analyses of 256 RNA sequencing datasets; a complete miRNAome; extensive metabolite information for representative varieties and wild relatives; epigenetic data on the genome-wide chromatin accessibility landscape; and various lettuce research papers published in the last decade. Five hierarchically accessible functions (Genome, Genotype, Germplasm, Phenotype, and O-Omics) have been developed with a user-friendly interface to enable convenient data access. Eight built-in tools (Assembly Converter, Search Gene, BLAST, JBrowse, Primer Design, Gene Annotation, Tissue Expression, Literature, and Data) are available for data downloading and browsing, functional gene exploration, and experimental practice. A community forum is also available for information sharing, and a summary of current research progress on different aspects of lettuce is included. We believe that LettuceGDB can be a comprehensive functional database amenable to data mining and database-driven exploration, useful for both scientific research and lettuce breeding.

Integrated analysis of lncRNAs, mRNAs, and TFs to identify network modules underlying diterpenoid biosynthesis in Salvia miltiorrhiza

Long non-coding RNAs (lncRNAs) are transcripts of more than 200 nucleotides (nt) in length, with minimal or no protein-coding capacity. Increasing evidence indicates that lncRNAs play important roles in the regulation of gene expression including in the biosynthesis of secondary metabolites. Salvia miltiorrhiza Bunge is an important medicinal plant in China. Diterpenoid tanshinones are one of the main active components of S. miltiorrhiza. To better understand the role of lncRNAs in regulating diterpenoid biosynthesis in S. miltiorrhiza, we integrated analysis of lncRNAs, mRNAs, and transcription factors (TFs) to identify network modules underlying diterpenoid biosynthesis based on transcriptomic data. In transcriptomic data, we obtained 6,651 candidate lncRNAs, 46 diterpenoid biosynthetic pathway genes, and 11 TFs involved in diterpenoid biosynthesis. Combining the co-expression and genomic location analysis, we obtained 23 candidate lncRNA-mRNA/TF pairs that were both co-expressed and co-located. To further observe the expression patterns of these 23 candidate gene pairs, we analyzed the time-series expression of S. miltiorrhiza induced by methyl jasmonate (MeJA). The results showed that 19 genes were differentially expressed at least a time-point, and four lncRNAs, two mRNAs, and two TFs formed three lncRNA-mRNA and/or TF network modules. This study revealed the relationship among lncRNAs, mRNAs, and TFs and provided new insight into the regulation of the biosynthetic pathway of S. miltiorrhiza diterpenoids.

Comprehensive analysis of the carboxylesterase gene reveals that NtCXE22 regulates axillary bud growth through strigolactone metabolism in tobacco

Carboxylesterases (CXE) are a class of hydrolytic enzymes with α/β-folding domains that play a vital role in plant growth, development, stress response, and activation of herbicide-active substances. In this study, 49 Nicotiana tabacum L. CXE genes (NtCXEs) were identified using a sequence homology search. The basic characteristics, phylogenetic evolution, gene structure, subcellular location, promoter cis-elements, and gene expression patterns of the CXE family were systematically analyzed. RNA-seq data and quantitative real-time PCR showed that the expression level of CXEs was associated with various stressors and hormones; gene expression levels were significantly different among the eight tissues examined and at different developmental periods. As a new class of hormones, strigolactones (SLs) are released from the roots of plants and can control the germination of axillary buds.NtCXE7, NtCXE9, NtCXE22, and NtCXE24 were homologous to Arabidopsis SLs hydrolase AtCXE15, and changes in their expression levels were induced by topping and by GR24 (a synthetic analogue of strigolactone). Further examination revealed that NtCXE22-mutant (ntcxe22) plants generated by CRISPR-Cas9 technology had shorter bud outgrowth with lower SLs content. Validation of NtCXE22 was also performed in NtCCD8-OE plants (with fewer axillary buds) and in ntccd8 mutant plants (with more axillary buds). The results suggest that NtCXE22 may act as an efficient SLs hydrolase and affects axillary bud development, thereby providing a feasible method for manipulating endogenous SLs in crops and ornamental plants.

Investigation and Expression Analysis of R2R3-MYBs and Anthocyanin Biosynthesis-Related Genes during Seed Color Development of Common Bean (Phaseolus vulgaris)

Anthocyanins are responsible for the coloration of common bean seeds, and their accumulation is positively correlated with the expression level of anthocyanin biosynthetic genes. The MBW (MYB-bHLH-WD40) complex is thought to regulate the expression of these genes, and MYB proteins, which are a key factor in activating anthocyanin pathway genes, have been identified in several plants. This study demonstrated gene structures, chromosomal placements, gene duplications of R2R3-MYBs, miRNAs associated with R2R3-MYBs, and the interaction of these genes with other flavonoid regulatory genes. qRT-PCR was used to investigate the role of specific R2R3-MYBs and flavonoid genes in common bean seed color development. As a result of a comprehensive analysis with the help of in silico tools, we identified 160 R2R3-MYB genes in the common bean genome. We divided these genes into 16 classes on the basis of their intron-exon and motif structures. Except for three, the rest of the common bean R2R3-MYB members were distributed to all chromosomes with different densities, primarily located on chromosomes 3 and 8. We identified a total of 44 duplicated gene pairs dispersed across 11 chromosomes and evolved under purifying selection (Ka/Ks  <  1), 19 of which were derived from a whole-genome duplication. Our research uncovered 25 putative repressor PvMYB proteins that contain the EAR motif. Additionally, fifty different cis-regulatory elements regulated by light, stress, and hormone were identified. Within the genome of the common bean, we discovered a total of 36 microRNAs that target a total of 72 R2R3-MYB transcripts. The effect of 16 R2R3-MYB genes and 16 phenylpropanoid pathway genes, selected on the basis of their interaction in the protein-protein interaction map, playing role in the regulation of seed coat color development was evaluated using qRT-PCR in 5 different tissues at different developmental stages. The results revealed that these specific genes have different expression levels during different developmental periods, with higher levels in the pod filling and early pod stages than in the rest of the developmental periods. Furthermore, it was shown that PvTT8 (bHLH), PvTT2 (PvMYB42), PvMYB113, PvTTG1, and PvWD68 genes have effects on the regulation of seed coat color. The findings of this study, which is the first to use whole-genome analysis to identify and characterize the R2R3-MYB genes in common bean, may serve as a reference for future functional research in the legume.

MIR390 Is Involved in Regulating Anthracnose Resistance in Apple

As an important cash crop in China, apple has a good flavor and is rich in nutrients. Fungal attacks have become a major obstacle in apple cultivation. Colletotrichum gloeosporioides is one of the most devastating fungal pathogens in apple. Thus, discovering resistance genes in response to C. gloeosporioides may aid in designing safer control strategies and facilitate the development of apple resistance breeding. A previous study reported that 'Hanfu' autotetraploid apple displayed higher C. gloeosporioides resistance than 'Hanfu' apple, and the expression level of mdm-MIR390b was significantly upregulated in autotetraploid plants compared to that in 'Hanfu' plants, as demonstrated by digital gene expression (DGE) analysis. It is still unclear, however, whether mdm-MIR390b regulates apple anthracnose resistance. Apple MIR390b was transformed into apple 'GL-3' plants to identify the functions of mdm-MIR390b in anthracnose resistance. C. gloeosporioides treatment analysis indicated that the overexpression of mdm-MIR390b reduced fungal damage to apple leaves and fruit. Physiology analysis showed that mdm-MIR390b increased C. gloeosporioides resistance by improving superoxide dismutase (SOD) and peroxidase (POD) activity to alleviate the damage caused by O₂^- and H₂O₂. Our results demonstrate that mdm-MIR390b can improve apple plants' anthracnose resistance.

What Do We Know about Barley miRNAs?

Plant miRNAs are powerful regulators of gene expression at the post-transcriptional level, which was repeatedly proved in several model plant species. miRNAs are considered to be key regulators of many developmental, homeostatic, and immune processes in plants. However, our understanding of plant miRNAs is still limited, despite the fact that an increasing number of studies have appeared. This systematic review aims to summarize our current knowledge about miRNAs in spring barley (Hordeum vulgare), which is an important agronomical crop worldwide and serves as a common monocot model for studying abiotic stress responses as well. This can help us to understand the connection between plant miRNAs and (not only) abiotic stresses in general. In the end, some future perspectives and open questions are summarized.

Identification of Species-Specific MicroRNAs Provides Insights into Dynamic Evolution of MicroRNAs in Plants

MicroRNAs (miRNAs) are an important class of regulatory small RNAs that program gene expression, mainly at the post-transcriptional level. Although sporadic examples of species-specific miRNAs (termed SS-miRNAs) have been reported, a genome-scale study across a variety of distant species has not been assessed. Here, by comprehensively analyzing miRNAs in 81 plant species phylogenetically ranging from chlorophytes to angiosperms, we identified 8048 species-specific miRNAs from 5499 families, representing over 61.2% of the miRNA families in the examined species. An analysis of the conservation from different taxonomic levels supported the high turnover rate of SS-miRNAs, even over short evolutionary distances. A comparison of the intrinsic features between SS-miRNAs and NSS-miRNAs (non-species-specific miRNAs) indicated that the AU content of mature miRNAs was the most striking difference. Our data further illustrated a significant bias of the genomic coordinates towards SS-miRNAs lying close to or within genes. By analyzing the 125,267 putative target genes for the 7966 miRNAs, we found the preferentially regulated functions of SS-miRNAs related to diverse metabolic processes. Collectively, these findings underscore the dynamic evolution of miRNAs in the species-specific lineages.

Identification and exploration of the GRF and GIF families in maize and foxtail millet

Plant growth and development are dependent on complex regulatory networks to adapt various environments. The growth regulatory factor (GRF) and GRF-interacting factor (GIF) families have been shown to control growth in various plant species. There are growing evidences that GRFs and GIFs can improve crop genetic transformation efficiency. In this study, we identified and classified 17 ZmGRFs, 10 SiGRFs, 4 ZmGIFs and 3 SiGIFs in maize (Zea mays L.) and foxtail millet (Setaria italica L.) using updated genome data. Many ABREs (Abscisic Acid-responsive elements) were present in the promoter regions of GRFs by analysis, and the expression levels of ZmGRF4, 9, 12, 14 and ZmGIF2 were associated with the Abscisic Acid (ABA) response. Furthermore, ZmGRF9 showed collinearity with AtGRF5 between Arabidopsis and maize. ZmGRF9 conservatively interacts with ZmGIF 2, 3, and 4. As a result, we systematically identified GRF and GIF family members, analyzed the regulatory network, and found that exogenous ABA inhibited the expression of GRFs, regulating responses to stress in the environment.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-022-01234-z.

Dynamically expressed small RNAs, substantially driven by genomic structural variants, contribute to transcriptomic changes during tomato domestication

Tomato has undergone extensive selections during domestication. Recent progress has shown that genomic structural variants (SVs) have contributed to gene expression dynamics during tomato domestication, resulting in changes of important traits. Here, we performed comprehensive analyses of small RNAs (sRNAs) from nine representative tomato accessions. We demonstrate that SVs substantially contribute to the dynamic expression of the three major classes of plant sRNAs: microRNAs (miRNAs), phased secondary short interfering RNAs (phasiRNAs), and 24-nucleotide heterochromatic siRNAs (hc-siRNAs). Changes in the abundance of phasiRNAs and 24-nucleotide hc-siRNAs likely contribute to the alteration of mRNA gene expression in cis during tomato domestication, particularly for genes associated with biotic and abiotic stress tolerance. We also observe that miRNA expression dynamics are associated with imprecise processing, alternative miRNA-miRNA* selections, and SVs. SVs mainly affect the expression of less-conserved miRNAs that do not have established regulatory functions or low abundant members in highly expressed miRNA families. Our data highlight different selection pressures on miRNAs compared to phasiRNAs and 24-nucleotide hc-siRNAs. Our findings provide insights into plant sRNA evolution as well as SV-based gene regulation during crop domestication. Furthermore, our dataset provides a rich resource for mining the sRNA regulatory network in tomato.

sRNAbench and sRNAtoolbox 2022 update: accurate miRNA and sncRNA profiling for model and non-model organisms

The NCBI Sequence Read Archive currently hosts microRNA sequencing data for over 800 different species, evidencing the existence of a broad taxonomic distribution in the field of small RNA research. Simultaneously, the number of samples per miRNA-seq study continues to increase resulting in a vast amount of data that requires accurate, fast and user-friendly analysis methods. Since the previous release of sRNAtoolbox in 2019, 55 000 sRNAbench jobs have been submitted which has motivated many improvements in its usability and the scope of the underlying annotation database. With this update, users can upload an unlimited number of samples or import them from Google Drive, Dropbox or URLs. Micro- and small RNA profiling can now be carried out using high-confidence Metazoan and plant specific databases, MirGeneDB and PmiREN respectively, together with genome assemblies and libraries from 441 Ensembl species. The new results page includes straightforward sample annotation to allow downstream differential expression analysis with sRNAde. Unassigned reads can also be explored by means of a new tool that performs mapping to microbial references, which can reveal contamination events or biologically meaningful findings as we describe in the example. sRNAtoolbox is available at: https://arn.ugr.es/srnatoolbox/.

The Construction and Exploration of a Comprehensive MicroRNA Centered Regulatory Network in Foxtail Millet (Setaria italica L.)

MicroRNA (miRNA) is an essential endogenous post-transcriptional regulatory factor, and foxtail millet (Setaria italica L.) is an ideal C4 model cereal that is a highly valuable crop in semiarid and arid areas. The Research on comprehensive and high confidence identification and annotation of foxtail millet miRNAs needs to be strengthened, and to our knowledge, there is no information on the regulatory network of foxtail millet miRNA. In this study, 136 high confidence miRNAs were identified through high-throughput sequencing of the small RNAs in seven tissues at the shooting and grain filling stages of foxtail millet. A total of 2,417 target genes were obtained by combining computational biology software and degradome sequencing methods. Furthermore, an analysis using transcriptome sequencing revealed the relationships between miRNAs and their target genes and simultaneously explored key regulatory modules in panicles during the grain filling stage. An miRNA regulatory network was constructed to explore the functions of miRNA in more detail. This network, centered on miRNAs and combining upstream transcriptional factors and downstream target genes, is primarily composed of feed forward loop motifs, which greatly enhances our knowledge of the potential functions of miRNAs and uncovers numerous previously unknown regulatory links. This study provides a solid foundation for research on the function and regulatory network of miRNAs in foxtail millet.

The 2022 Nucleic Acids Research database issue and the online molecular biology database collection

The 2022 Nucleic Acids Research Database Issue contains 185 papers, including 87 papers reporting on new databases and 85 updates from resources previously published in the Issue. Thirteen additional manuscripts provide updates on databases most recently published elsewhere. Seven new databases focus specifically on COVID-19 and SARS-CoV-2, including SCoV2-MD, the first of the Issue's Breakthrough Articles. Major nucleic acid databases reporting updates include MODOMICS, JASPAR and miRTarBase. The AlphaFold Protein Structure Database, described in the second Breakthrough Article, is the stand-out in the protein section, where the Human Proteoform Atlas and GproteinDb are other notable new arrivals. Updates from DisProt, FuzDB and ELM comprehensively cover disordered proteins. Under the metabolism and signalling section Reactome, ConsensusPathDB, HMDB and CAZy are major returning resources. In microbial and viral genomes taxonomy and systematics are well covered by LPSN, TYGS and GTDB. Genomics resources include Ensembl, Ensembl Genomes and UCSC Genome Browser. Major returning pharmacology resource names include the IUPHAR/BPS guide and the Therapeutic Target Database. New plant databases include PlantGSAD for gene lists and qPTMplants for post-translational modifications. The entire Database Issue is freely available online on the Nucleic Acids Research website (https://academic.oup.com/nar). Our latest update to the NAR online Molecular Biology Database Collection brings the total number of entries to 1645. Following last year's major cleanup, we have updated 317 entries, listing 89 new resources and trimming 80 discontinued URLs. The current release is available at http://www.oxfordjournals.org/nar/database/c/.