Segmental and Tandem Duplications Driving the Recent NBS-LRR Gene Expansion in the Asparagus Genome

refseqR: an R package for common computational operations with records on RefSeq collection

Summary

We introduce refseqR, an R package that offers a user-friendly solution, enabling common computational operations on RefSeq entries (GenBank, NCBI). The package is specifically designed to interact with records curated from the RefSeq database. Most importantly, the interoperability and integration with several Bioconductor objects allow connections to be applied to other projects.

Availability and implementation

The package refseqR is implemented in R and published under the MIT open-source license. The source code, documentation, and usage instructions are available on CRAN (https://CRAN.R-project.org/package=refseqR).

Evolution patterns of NBS genes in the genus Dendrobium and NBS-LRR gene expression in D. officinale by salicylic acid treatment

BACKGROUND

Dendrobium officinale Kimura et Migo, which contains rich polysaccharides, flavonoids and alkaloids, is a Traditional Chinese Medicine (TCM) with important economic benefits, while various pathogens have brought huge losses to its industrialization. NBS gene family is the largest class of plant disease resistance (R) genes, proteins of which are widely distributed in the upstream and downstream of the plant immune systems and are responsible for receiving infection signals and regulating gene expression respectively. It is of great significance for the subsequent disease resistance breeding of D. officinale to identify NBS genes by using the newly published high-quality chromosome-level D. officinale genome.

RESULTS

In this study, a total of 655 NBS genes were uncovered from the genomes of D. officinale, D. nobile, D. chrysotoxum, V. planifolia, A. shenzhenica, P. equestris and A. thaliana. The phylogenetic results of CNL-type protein sequences showed that orchid NBS-LRR genes have significantly degenerated on branches a and b. The Dendrobium NBS gene homology analysis showed that the Dendrobium NBS genes have two obvious characteristics: type changing and NB-ARC domain degeneration. Because the NBS-LRR genes have both NB-ARC and LRR domains, 22 D. officinale NBS-LRR genes were used for subsequent analyses, such as gene structures, conserved motifs, cis-elements and functional annotation analyses. All these results suggested that D. officinale NBS-LRR genes take part in the ETI system, plant hormone signal transduction pathway and Ras signaling pathway. Finally, there were 1,677 DEGs identified from the salicylic acid (SA) treatment transcriptome data of D. officinale. Among them, six NBS-LRR genes (Dof013264, Dof020566, Dof019188, Dof019191, Dof020138 and Dof020707) were significantly up-regulated. However, only Dof020138 was closely related to other pathways from the results of WGCNA, such as pathogen identification pathways, MAPK signaling pathways, plant hormone signal transduction pathways, biosynthetic pathways and energy metabolism pathways.

CONCLUSION

Our results revealed that the NBS gene degenerations are common in the genus Dendrobium, which is the main reason for the diversity of NBS genes, and the NBS-LRR genes generally take part in D. officinale ETI system and signal transduction pathways. In addition, the D. officinale NBS-LRR gene Dof020138, which may have an important breeding value, is indirectly activated by SA in the ETI system.

Genome-Wide Identification of P450 Genes in Chironomid Propsilocerus akamusi Reveals Candidate Genes Involved in Gut Microbiota-Mediated Detoxification of Chlorpyrifos

Chironomids commonly dominate macroinvertebrate assemblages in aquatic habitats and these non-biting midges are known for their ability to tolerate contaminants. Studies regarding the interplay between gut microbiota and host detoxification ability is currently a point of interest. Cytochrome P450s (P450s) are critical metabolic enzymes in which a subset is involved in xenobiotic detoxification. In this study, we first conducted an integrated global investigation of P450s based on the whole genomic sequence of Propsilocerus akamusi and retrieved a series of 64 P450 genes which were further classified into 4 clans and 25 families on the basis of phylogenetic relationships. With assistance of RNA-Seq and RT-qPCR validation, the expression profile of screened PaP450s in guts was compared between chlorpyrifos-challenged larvae with deficient gut microbiota (GD) and those with a conventional gut community (CV). An increasing prevalence of chlorpyrifos from sublethal to lethal dosages induced a greater mortality rate of individuals coupled with remarkable downregulation of 14 P450s in GD larval guts when compared to CV ones. Moreover, it turned out that the decreased level of PaCYP3998B1 and PaCYP3987D1 might imply impaired host endogenous detoxification capability potentiated by gut dysbiosis, reflected by a remarkably severe mortality in GD larvae treated with lethal chlorpyrifos. Collectively, our study unveiled candidate P450 genes that might be mediated by gut symbionts in chlorpyrifos-challenged P. akamusi larvae, possibly facilitating further understanding of the detoxified mechanism that chironomids might employ to alleviate poisonousness.

Genome-Wide Characterization of the Aquaporin Gene Family in Radish and Functional Analysis of RsPIP2-6 Involved in Salt Stress

Aquaporins (AQPs) constitute a highly diverse family of channel proteins that transport water and neutral solutes. AQPs play crucial roles in plant development and stress responses. However, the characterization and biological functions of RsAQPs in radish (Raphanus sativus L.) remain elusive. In this study, 61 non-redundant members of AQP-encoding genes were identified from the radish genome database and located on nine chromosomes. Radish AQPs (RsAQPs) were divided into four subfamilies, including 21 plasma membrane intrinsic proteins (PIPs), 19 tonoplast intrinsic proteins (TIPs), 16 NOD-like intrinsic proteins (NIPs), and 5 small basic intrinsic proteins (SIPs), through phylogenetic analysis. All RsAQPs contained highly conserved motifs (motifs 1 and 4) and transmembrane regions, indicating the potential transmembrane transport function of RsAQPs. Tissue- and stage-specific expression patterns of AQP gene analysis based on RNA-seq data revealed that the expression levels of PIPs were generally higher than TIPs, NIPs, and SIPs in radish. In addition, quantitative real-time polymerase chain reaction (qRT-PCR) revealed that seven selected RsPIPs, according to our previous transcriptome data (e.g., RsPIP1-3, 1-6, 2-1, 2-6, 2-10, 2-13, and 2-14), exhibited significant upregulation in roots of salt-tolerant radish genotype. In particular, the transcriptional levels of RsPIP2-6 dramatically increased after 6 h of 150 mM NaCl treatment during the taproot thickening stage. Additionally, overexpression of RsPIP2-6 could enhance salt tolerance by Agrobacterium rhizogenes-mediated transgenic radish hairy roots, which exhibited the mitigatory effects of plant growth reduction, leaf relative water content (RWC) reduction and alleviation of O2- in cells, as shown by nitro blue tetrazolium (NBT) staining, under salt stress. These findings are helpful for deeply dissecting the biological function of RsAQPs on the salt stress response, facilitating practical application and genetic improvement of abiotic stress resistance in radish.

Genome-wide analysis of pathogenesis-related protein 1 (PR-1) gene family from Musa spp. and its role in defense response during stresses

Pathogenesis related protein-1 (PR-1) is the most abundantly produced protein during defense response against many biotic and abiotic stresses. However, knowledge on PR-1 gene family and its evolutionary relationship in banana is very limited. In order to study the potential role of PR-1 genes in banana, genome wide identification, structure analysis and expressions were performed. A total of 15 and 11 PR-1 genes were identified from A and B genomes of banana and the proteins encoded by this gene family are of varying lengths and harbor conserved domains and motifs. PR-1 genes are unevenly dispersed on 11 chromosomes with segmental duplication in both A and B genome, suggesting an important contribution of duplication in expansion of PR-1 gene family in banana. qRT-PCR analysis of PR-1 gene showed positive correlation with the RNAseq data under various stresses and examination of expression pattern of selected MaPR-1 genes in banana revealed its role in biotic and abiotic stresses in general and fusarium wilt in particular. This study provides significant insight into the functions of PR-1 genes which can be further exploited as a promising candidate for developing multiple stress tolerant banana varieties.

Identification and Characterization of NBS Resistance Genes in Akebia trifoliata

Akebia trifoliata is an important multiuse perennial plant that often suffers attacks from various pathogens due to its long growth cycle, seriously affecting its commercial value. The absence of research on the resistance (R) genes of A. trifoliata has greatly limited progress in the breeding of resistant varieties. Genes encoding proteins containing nucleotide binding sites (NBSs) and C-terminal leucine-rich repeats (LRRs), the largest family of plant resistance (R) genes, are vital for plant disease resistance. A comprehensive genome-wide analysis showed that there were only 73 NBS genes in the A. trifoliata genome, including three main subfamilies (50 coiled coil (CC)-NBS-LRR (CNL), 19 Toll/interleukin-1 receptor (TIR)-NBS-LRR (TNL) and four resistance to powdery mildew8 (RPW8)-NBS-LRR (RNL) genes). Additionally, 64 mapped NBS candidates were unevenly distributed on 14 chromosomes, most of which were assigned to the chromosome ends; 41 of these genes were located in clusters, and the remaining 23 genes were singletons. Both the CNLs and TNLs were further divided into four subgroups, and the CNLs had fewer exons than the TNLs. Structurally, all eight previously reported conserved motifs were identified in the NBS domains, and both their order and their amino acid sequences exhibited high conservation. Evolutionarily, tandem and dispersed duplications were shown to be the two main forces responsible for NBS expansion, producing 33 and 29 genes, respectively. A transcriptome analysis of three fruit tissues at four developmental stages showed that NBS genes were generally expressed at low levels, while a few of these genes showed relatively high expression during later development in rind tissues. Overall, this research is the first to identify and characterize A. trifoliata NBS genes and is valuable for both the development of new resistant cultivars and the study of molecular mechanisms of resistance.

Comparative Metabolome and Transcriptome Analyses of Susceptible Asparagus officinalis and Resistant Wild A. kiusianus Reveal Insights into Stem Blight Disease Resistance

Phomopsis asparagi is one of the most serious fungal pathogens, which causes stem blight disease in Asparagus officinalis (AO), adversely affecting its production worldwide. Recently, the development of novel asparagus varieties using wild Asparagus genetic resources with natural P. asparagi resistance has become a priority in Japan due to the lack of resistant commercial AO cultivars. In this study, comparative metabolome and transcriptome analyses of susceptible AO and resistant wild Asparagus kiusianus (AK) 24 and 48 h postinoculated (AOI_24 hpi, AOI_48 hpi, AKI_24 hpi and AKI_48 hpi, respectively) with P. asparagi were conducted to gain insights into metabolic and expression changes associated with AK species. Following infection, the resistant wild AK showed rapid metabolic changes with increased levels of flavonoids and steroidal saponins and decreased asparagusic acid glucose ester content, compared with the susceptible AO plants. Transcriptome data revealed a total of 21  differentially expressed genes (DEGs) as the core gene set that displayed upregulation in the resistant AK versus susceptible AO after infection with P. asparagi. Kyoto Encyclopedia of Genes and Genomes pathway analysis of these DEGs identified 11 significantly enriched pathways, including flavonoid biosynthesis and primary metabolite metabolism, in addition to plant signaling and defense-related pathways. In addition, comparative single-nucleotide polymorphism and Indel distributions in susceptible AO and resistant AK plants were evaluated using the latest AO reference genome Aspof.V1. The data generated in this study are important resources for advancing Asparagus breeding programs and for investigations of genetic linkage mapping, phylogenetic diversity and plant defense-related genes.

Genome-Wide Identification and Evolutionary Analysis of NBS-LRR Genes From Dioscorea rotundata

Dioscorea rotundata is an important food crop that is mainly cultivated in subtropical regions of the world. D. rotundata is frequently infected by various pathogens during its lifespan, which results in a substantial economic loss in terms of yield and quality. The disease resistance gene (R gene) profile of D. rotundata is largely unknown, which has greatly hampered molecular study of disease resistance in this species. Nucleotide-binding site–leucine-rich repeat (NBS-LRR) genes are the largest group of plant R genes, and they play important roles in plant defense responses to various pathogens. In this study, 167 NBS-LRR genes were identified from the D. rotundata genome. Subsequently, one gene was assigned to the resistance to powdery mildew8 (RPW8)-NBS-LRR (RNL) subclass and the other 166 genes to the coiled coil (CC)-NBS-LRR (CNL) subclass. None of the Toll/interleukin-1 receptor (TIR)-NBS-LRR (TNL) genes were detected in the genome. Among them, 124 genes are located in 25 multigene clusters and 43 genes are singletons. Tandem duplication serves as the major force for the cluster arrangement of NBS-LRR genes. Segmental duplication was detected for 18 NBS-LRR genes, although no whole-genome duplication has been documented for the species. Phylogenetic analysis revealed that D. rotundata NBS-LRR genes share 15 ancestral lineages with Arabidopsis thaliana genes. The NBS-LRR gene number increased by more than a factor of 10 during D. rotundata evolution. A conservatively evolved ancestral lineage was identified from D. rotundata, which is orthologs to the Arabidopsis RPM1 gene. Transcriptome analysis for four different tissues of D. rotundata revealed a low expression of most NBS-LRR genes, with the tuber and leaf displaying a relatively high NBS-LRR gene expression than the stem and flower. Overall, this study provides a complete set of NBS-LRR genes for D. rotundata, which may serve as a fundamental resource for mining functional NBS-LRR genes against various pathogens.

Chromosome-level genome assembly, annotation and evolutionary analysis of the ornamental plant Asparagus setaceus

Asparagus setaceus is a popular ornamental plant cultivated in tropical and subtropical regions globally. Here, we constructed a chromosome-scale reference genome of A. setaceus to facilitate the investigation of its genome characteristics and evolution. Using a combination of Nanopore long reads, Illumina short reads, 10× Genomics linked reads, and Hi-C data, we generated a high-quality genome assembly of A. setaceus covering 710.15 Mb, accounting for 98.63% of the estimated genome size. A total of 96.85% of the sequences were anchored to ten superscaffolds corresponding to the ten chromosomes. The genome of A. setaceus was predicted to contain 28,410 genes, 25,649 (90.28%) of which were functionally annotated. A total of 65.59% of the genome was occupied by repetitive sequences, among which long terminal repeats were predominant (42.51% of the whole genome). Evolutionary analysis revealed an estimated divergence time of A. setaceus from its close relative A. officinalis of ~9.66 million years ago, and A. setaceus underwent two rounds of whole-genome duplication. In addition, 762 specific gene families, 96 positively selected genes, and 76 resistance (R) genes were detected and functionally predicted in A. setaceus. These findings provide new knowledge about the characteristics and evolution of the A. setaceus genome, and will facilitate comparative genetic and genomic research on the genus Asparagus.

Population Genetics of the Highly Polymorphic RPP8 Gene Family

Plant nucleotide-binding domain and leucine-rich repeat containing (NLR) genes provide some of the most extreme examples of polymorphism in eukaryotic genomes, rivalling even the vertebrate major histocompatibility complex. Surprisingly, this is also true in Arabidopsis thaliana, a predominantly selfing species with low heterozygosity. Here, we investigate how gene duplication and intergenic exchange contribute to this extraordinary variation. RPP8 is a three-locus system that is configured chromosomally as either a direct-repeat tandem duplication or as a single copy locus, plus a locus 2 Mb distant. We sequenced 48 RPP8 alleles from 37 accessions of A. thaliana and 12 RPP8 alleles from Arabidopsis lyrata to investigate the patterns of interlocus shared variation. The tandem duplicates display fixed differences and share less variation with each other than either shares with the distant paralog. A high level of shared polymorphism among alleles at one of the tandem duplicates, the single-copy locus and the distal locus, must involve both classical crossing over and intergenic gene conversion. Despite these polymorphism-enhancing mechanisms, the observed nucleotide diversity could not be replicated under neutral forward-in-time simulations. Only by adding balancing selection to the simulations do they approach the level of polymorphism observed at RPP8. In this NLR gene triad, genetic architecture, gene function and selection all combine to generate diversity.