Crested wheatgrass (Agropyron cristatum L.) breeding programs aim to develop later maturing cultivars for extending early spring grazing in Western Canada. Plant maturity is a complex genetic trait, and little is known about genes associated with late maturity in this species. An attempt was made using RNA-Seq to profile the transcriptome of crested wheatgrass maturity and to analyze differentially expressed genes (DEGs) between early and late maturing lines. Three cDNA libraries for each line were generated by sampling leaves at the stem elongation stage, spikes at the boot and anthesis stages. A total of 75,218,230 and 74,015,092 clean sequence reads were obtained for early and late maturing lines, respectively. De novo assembly of all sequence reads generated 401,587 transcripts with a mean length of 546 bp and N50 length of 691 bp. Out of 13,133 DEGs detected, 22, 17, and eight flowering related DEGs were identified for the three stages, respectively. Twelve DEGs, including nine flowering related DEGs at the stem elongation stage were further confirmed by qRT-PCR. The analysis of homologous genes of the photoperiod pathway revealed their lower expression in the late maturing line at the stem elongation stage, suggesting that their differential expression contributed to late maturity in crested wheatgrass.

Small interfering RNA (siRNA) duplexes are short (usually 21 to 24 bp) double-stranded RNAs (dsRNAs) with several overhanging nucleotides at both 5'- and 3'-ends. It has been found that siRNA duplexes bind the RNA-induced silencing complex (RISC) and cleave the sense strands with endonucleases. In this study, for the first time, we detected siRNA duplexes induced by plant viruses on a large scale using next-generation sequencing (NGS) data. In addition, we used the detected 21 nucleotide (nt) siRNA duplexes with 2 nt overhangs to construct a dataset for future data mining. The analytical results of the features in the detected siRNA duplexes were consistent with those from previous studies. The investigation of siRNA duplexes is useful for a better understanding of the RNA interference (RNAi) mechanism. It can also help to improve the virus detection based on the small RNA sequencing (sRNA-seq) technologies and to rationally design siRNAs for RNAi experiments.

Aminoacyl-tRNA synthetases (aaRSs) play essential roles in protein translation. In addition, numerous aaRSs (mostly in vertebrates) have also been discovered to possess a range of non-canonical functions. Very few studies have been conducted to elucidate or characterize non-canonical functions of plant aaRSs. A genome-wide search for aaRS genes in Arabidopsis thaliana revealed a total of 59 aaRS genes. Among them, asparaginyl-tRNA synthetase (AsnRS) was found to possess a WHEP domain inserted into the catalytic domain in a plant-specific manner. This insertion was observed only in the cytosolic isoform. In addition, a long stretch of sequence that exhibited weak homology with histidine ammonia lyase (HAL) was found at the N-terminus of histidyl-tRNA synthetase (HisRS). This HAL-like domain has only been seen in plant HisRS, and only in cytosolic isoforms. Additionally, a number of genes lacking minor or major portions of the full-length aaRS sequence were found. These genes encode 14 aaRS fragments that lack key active site sequences and are likely catalytically null. These identified genes that encode plant-specific additional domains or aaRS fragment sequences are candidates for aaRSs possessing non-canonical functions.

Mgaloblishvili, a Vitis vinifera cultivar, exhibits unique resistance traits against Plasmopara viticola, the downy mildew agent. This offers the unique opportunity of exploring the molecular responses in compatible and incompatible plant-pathogen interaction. In this study, whole transcriptomes of Mgaloblishvili, Pinot noir (a V. vinifera susceptible cultivar), and Bianca (a resistant hybrid) leaves, inoculated and non-inoculated with the pathogen, were used to identify P. viticola effector-encoding genes and plant susceptibility/resistance genes. Multiple effector-encoding genes were identified in P. viticola transcriptome, with remarkable expression differences in relation to the inoculated grapevine cultivar. Intriguingly, five apoplastic effectors specifically associated with resistance in V. vinifera. Gene coexpression network analysis identified specific modules and metabolic changes occurring during infection in the three grapevine cultivars. Analysis of these data allowed, for the first time, the detection in V. vinifera of a putative P. viticola susceptibility gene, encoding a LOB domain-containing protein. Finally, the de novo assembly of Mgaloblishvili, Pinot noir, and Bianca transcriptomes and their comparison highlighted novel candidate genes that might be at the basis of the resistant phenotype. These results open the way to functional analysis studies and to new perspectives in molecular breeding of grapevine for resistance to P. viticola.

Transcription factors are important regulators of gene expression. They can orchestrate the activation or repression of hundreds or thousands of genes and control diverse processes in a coordinated way. This work explores the effect of a master regulator of plant development, BOLITA (BOL), in plant metabolism, with a special focus on specialized metabolism. For this, we used an Arabidopsis thaliana line in which the transcription factor activity can be induced. Fingerprinting metabolomic analyses of whole plantlets were performed at different times after induction. After 96 h, all induced replicas clustered as a single group, in contrast with all controls which did not cluster. Metabolomic analyses of shoot and root tissues enabled the putative identification of differentially accumulated metabolites in each tissue. Finally, the analysis of global gene expression in induced vs. non-induced root samples, together with enrichment analyses, allowed the identification of enriched metabolic pathways among the differentially expressed genes and accumulated metabolites after the induction. We concluded that the induction of BOL activity can modify the Arabidopsis metabolome. Future work should investigate whether its action is direct or indirect, and the implications of the metabolic changes for development regulation and bioprospection.

miR319 was the first plant miRNA discovered via forward genetic mutation screening. In this study, we found that miR319 family members had similar sequences but different expression patterns in Brassica campestris and Arabidopsis thaliana. RT-PCR analysis revealed that Bra-MIR319a and Bra-MIR319c had similar expression patterns and were widely expressed in plant development, whereas Bra-MIR319b could only be detected in stems. The overexpression of each Bra-MIR319 family member in Arabidopsis could inhibit cell division and function in leaf and petal morphogenesis. Bra-miR319a formed a new regulatory relationship after whole genome triplication, and Bra-MIR319a overexpressing in Arabidopsis led to the degradation of pollen content and affected the formation of intine, thereby causing pollen abortion. Our results suggest that Bra-MIR319 family members have functional similarity and difference in plant development.

Agave species are important crassulacean acid metabolism (CAM) plants and widely cultivated in tropical areas for producing tequila spirit and fiber. The hybrid H11648 of Agave ((A. amaniensis × A. angustifolia) × A. amaniensis) is the main cultivar for fiber production in Brazil, China, and African countries. Small Auxin Up-regulated RNA (SAUR) genes have broad effect on auxin signaling-regulated plant growth and development, while only few SAUR genes have been reported in Agave species. In this study, we identified 43, 60, 24, and 21 SAUR genes with full-length coding regions in A. deserti, A. tequilana, A. H11648, and A. americana, respectively. Although phylogenetic analysis revealed that rice contained a species-specific expansion pattern of SAUR gene, no similar phenomena were observed in Agave species. The in silico expression indicated that SAUR genes had a distinct expression pattern in A. H11648 compared with other Agave species; and four SAUR genes were differentially expressed during CAM diel cycle in A. americana. Additionally, an expression analysis was conducted to estimate SAUR gene expression during different leaf developmental stages, abiotic and biotic stresses in A. H11648. Together, we first characterized the SAUR genes of Agave based on previously published transcriptome datasets and emphasized the potential functions of SAUR genes in Agave's leaf development and stress responses. The identification of which further expands our understanding on auxin signaling-regulated plant growth and development in Agave species.

In host-pathogen interactions RNA interference (RNAi) has emerged as a pivotal mechanism to modify both, the immune responses of the host as well as the pathogenicity and virulence of the pathogen. In addition, in some fungi RNAi is also known to affect chromosome biology via its effect on chromatin conformation. Previous studies reported no effect of the RNAi machinery on the virulence of the fungal plant pathogen Zymoseptoria tritici however the role of RNAi is still poorly understood in this species. Herein, we elucidate whether the RNAi machinery is conserved within the genus Zymoseptoria. Moreover, we conduct functional analyses of Argonaute and Dicer-like proteins and test if the RNAi machinery affects chromosome stability. We show that the RNAi machinery is conserved among closely related Zymoseptoria species while an exceptional pattern of allelic diversity was possibly caused by introgression. The deletion of Ago1 reduced the ability of the fungus to produce asexual propagules in planta in a quantitative matter. Chromosome stability of the accessory chromosome of Z. tritici was not prominently affected by the RNAi machinery. These results indicate, in contrast to previous finding, a role of the RNAi pathway during host infection, but not in the stability of accessory chromosomes in Z. tritici.

Since its inception, RNA sequencing (RNA-seq) has become the most effective way to study gene expression. After more than a decade of development, numerous RNA-seq datasets have been created, and the full utilization of these datasets has emerged as a major issue. In this study, we built a comprehensive database named Grape-RNA, which is focused on the collection, evaluation, treatment, and data sharing of grape RNA-seq datasets. This database contains 1529 RNA-seq samples, 112 microRNA samples from the public platform, and 485 RNA-seq in-house datasets sequenced by our lab. We classified these data into 25 conditions and provide the sample information, cleaned raw data, expression level, assembled unigenes, useful tools, and other relevant information to the users. Thus, this study provides data and tools that should be beneficial for researchers by allowing them to easily use the RNA-seq. The provided information can greatly contribute to grape breeding and genomic and biological research. This study may improve the usage of RNA-seq.

Selection is expected to work differently in autosomal and X-linked genes because of their ploidy difference and the exposure of recessive X-linked mutations to haploid selection in males. However, it is not clear whether these expectations apply to recently evolved sex chromosomes, where many genes retain functional X- and Y-linked gametologs. We took advantage of the recently evolved sex chromosomes in the plant Silene latifolia and its closely related species to compare the selective pressures between hemizygous and non-hemizygous X-linked genes as well as between X-linked genes and autosomal genes. Our analysis, based on over 1000 genes, demonstrated that, similar to animals, X-linked genes in Silene evolve significantly faster than autosomal genes—the so-called faster-X effect. Contrary to expectations, faster-X divergence was detectable only for non-hemizygous X-linked genes. Our phylogeny-based analyses of selection revealed no evidence for faster adaptation in X-linked genes compared to autosomal genes. On the other hand, partial relaxation of purifying selection was apparent on the X-chromosome compared to the autosomes, consistent with a smaller genetic diversity in S. latifolia X-linked genes (πx = 0.016; πaut = 0.023). Thus, the faster-X divergence in S. latifolia appears to be a consequence of the smaller effective population size rather than of a faster adaptive evolution on the X-chromosome. We argue that this may be a general feature of “young” sex chromosomes, where the majority of X-linked genes are not hemizygous, preventing haploid selection in heterogametic sex.

Programmed cell death (PCD) is a critical process in plant immunity, enabling the targeted elimination of infected cells to prevent the spread of pathogens. The tight regulation of PCD within plant cells is well-documented; however, specific mechanisms remain elusive or controversial. Heterotrimeric G proteins are multifunctional signaling elements consisting of three distinct subunits, Gα, Gβ, and Gγ. In Arabidopsis, the Gβγ dimer serves as a positive regulator of plant defense. Conversely, in species such as rice, maize, cotton, and tomato, mutants deficient in Gβ exhibit constitutively active defense responses, suggesting a contrasting negative role for Gβ in defense mechanisms within these plants. Using a transient overexpression approach in addition to knockout mutants, we observed that Gβγ enhanced cell death progression and elevated the accumulation of reactive oxygen species in a similar manner across Arabidopsis, tomato, and Nicotiana benthamiana, suggesting a conserved G protein role in PCD regulation among diverse plant species. The enhancement of PCD progression was cooperatively regulated by Gβγ and one Gα, XLG2. We hypothesize that G proteins participate in two distinct mechanisms regulating the initiation and progression of PCD in plants. We speculate that G proteins may act as guardees, the absence of which triggers PCD. However, in Arabidopsis, this G protein guarding mechanism appears to have been lost in the course of evolution.

Phyllosticta yuccae is an important plant pathogen causing leaf spot disease in Yucca gigantea Lem. It is imperative to note that the amount of information available about the mitogenome of this subject is severely limited. This must be addressed immediately, as it is crucial to our understanding and progress in this field. To better understand the mitogenomic characteristics of P. yuccae, we conducted its sequencing by MGISEQ. Afterwards, the mitogenome was assembled and annotated. The mitogenomic characteristics and phylogenetic placement of the P. yuccae strain KUMCC 6213 were analyzed. The study revealed that the mitogenome of P. yuccae is a circular DNA molecule, consisting of 178,540 base pairs. It contains a total of 64 genes, including 14 protein-coding genes (PCGs), 26 transfer RNA genes (tRNA), 2 ribosomal RNA genes (rRNA), and 22 open reading frame genes (ORF), accounting for 80.98% of the total size. Repetitive sequences accounted for 15.42% of the mitogenome. The analysis of codon usage indicated that the codon UUA was the most commonly utilized, whereas the amino acid Leu was the most frequently employed. A comparative analysis of mitogenomes between P. yuccae and Macrophomina phaseolina (Tassi) Goid. showed notable variations in the position and size of gene clusters, with cox1, nad4, and nad4L genes exhibiting relatively low conservation. Phylogenetic analysis based on the 14 PCGs revealed that P. yuccae has the closest genetic relationship with M. phaseolina (Botryosphaeriaceae, Botryosphaeriales). This study first reports the mitogenome of P. yuccae and validates its phylogenetic placement. The findings enhance the knowledge of mitogenomes in Botryosphaeriales, offering novel perspectives on the genetics and evolution of the plant pathogen P. yuccae. This is crucial for the accurate prevention and management of leaf spot disease in Y. gigantea.

Plant-parasitic nematodes cause extensive annual yield losses to worldwide agricultural production. Most cultivated plants have no known resistance against nematodes and the few bearing a resistance gene can be overcome by certain species. Chemical methods that have been deployed to control nematodes have largely been banned from use due to their poor specificity and high toxicity. Hence, there is an urgent need for the development of cleaner and more specific control methods. Recent advances in nematode genomics, including in phytoparasitic species, provide an unprecedented opportunity to identify genes and functions specific to these pests. Using phylogenomics, we compared 61 nematode genomes, including 16 for plant-parasitic species and identified more than 24,000 protein families specific to these parasites. In the genome of Meloidogyne incognita, one of the most devastating plant parasites, we found ca. 10,000 proteins with orthologs restricted only to phytoparasitic species and no further homology in protein databases. Among these phytoparasite-specific proteins, ca. 1000 shared the same properties as known secreted effectors involved in essential parasitic functions. Of these, 68 were novel and showed strong expression during the endophytic phase of the nematode life cycle, based on both RNA-seq and RT-qPCR analyses. Besides effector candidates, transcription-related and neuro-perception functions were enriched in phytoparasite-specific proteins, revealing interesting targets for nematode control methods. This phylogenomics analysis constitutes a unique resource for the further understanding of the genetic basis of nematode adaptation to phytoparasitism and for the development of more efficient control methods.

Insect salivary glands play an important role for host feeding, specifically by secreting salivary proteins for digestion and potentially modulating host defenses. Compared to other hemipterans, the significance of salivary glands is less studied in the black-faced leafhopper, Graminella nigrifrons, a crop pest that vectors several agronomically important plant viruses. To identify functionally important genes in the salivary glands of the black-faced leafhopper, we compared transcriptomes between adult salivary glands (SG) and the remaining carcasses. We identified 14,297 salivary gland-enriched transcripts and 195 predicted secretory peptides (i.e., with a signal peptide and extracellular localization characteristics). Overall, the SG transcriptome included functions such as 'oxidoreduction', 'membrane transport', and 'ATP-binding', which might be important for the fundamental physiology of this tissue. We further evaluated transcripts with potential contributions in host feeding using RT-qPCR. Two SG-enriched transcripts (log2 fold change > 5), GnP19 and GnE63 (a putative calcium binding protein), were significantly upregulated in maize-fed adults relative to starved adults, validating their importance in feeding. The SG-enriched transcripts of the black-faced leafhopper could play a potential role for interacting with maize and could be targets of interest for further functional studies and improve pest control and disease transmission.

Vacuolar-type H+-ATPase (V-ATPase), a multisubunit proton pump located on the endomembrane, plays an important role in plant growth. The Arabidopsis thaliana V-ATPase d subunit (VHA-d) consists of two isoforms; AtVHA-d1 and AtVHA-d2. In this study, the function of AtVHA-d2 was investigated. Histochemical analysis revealed that the expression of AtVHA-d1 and AtVHA-d2 was generally highly overlapping in multiple tissues at different developmental stages of Arabidopsis. Subcellular localization revealed that AtVHA-d2 was mainly localized to the vacuole. AtVHA-d2 expression was significantly induced by oxidative stress. Analysis of phenotypic and H2O2 content showed that the atvha-d2 mutant was sensitive to oxidative stress. The noninvasive microtest monitoring demonstrated that the net H+ influx in the atvha-d2 roots was weaker than that in the wild-type under normal conditions. However, oxidative stress resulted in the H+ efflux in atvha-d2 roots, which was significantly different from that in the wild-type. RNA-seq combined with qPCR analysis showed that the expression of several members of the plasma membrane H+-ATPase gene (AtAHA) family in atvha-d2 was significantly different from that in the wild-type. Overall, our results indicate that AtVHA-d2 plays a role in Arabidopsis in response to oxidative stress by affecting H+ flux and AtAHA gene expression.

Drought stress requires plants to adjust their water balance to maintain tissue water levels. Isohydric plants ('water-savers') typically achieve this through stomatal closure, while anisohydric plants ('water-wasters') use osmotic adjustment and maintain stomatal conductance. Isohydry or anisohydry allows plant species to adapt to different environments. In this paper we show that both mechanisms occur in bread wheat (Triticum aestivum L.). Wheat lines with reproductive drought-tolerance delay stomatal closure and are temporarily anisohydric, before closing stomata and become isohydric at higher threshold levels of drought stress. Drought-sensitive wheat is isohydric from the start of the drought treatment. The capacity of the drought-tolerant line to maintain stomatal conductance correlates with repression of ABA synthesis in spikes and flag leaves. Gene expression profiling revealed major differences in the drought response in spikes and flag leaves of both wheat lines. While the isohydric drought-sensitive line enters a passive growth mode (arrest of photosynthesis, protein translation), the tolerant line mounts a stronger stress defence response (ROS protection, LEA proteins, cuticle synthesis). The drought response of the tolerant line is characterised by a strong response in the spike, displaying enrichment of genes involved in auxin, cytokinin and ethylene metabolism/signalling. While isohydry may offer advantages for longer term drought stress, anisohydry may be more beneficial when drought stress occurs during the critical stages of wheat spike development, ultimately improving grain yield.

Drechslera teres (D. teres) is an ascomycete, responsible for net blotch, the most serious barley disease causing an important economic impact. The cell wall is a crucial structure for the growth and development of fungi. Thus, understanding cell wall structure, composition and biosynthesis can help in designing new strategies for pest management. Despite the severity and economic impact of net blotch, this is the first study analyzing the cell wall-related genes in D. teres. We have identified key genes involved in the synthesis/remodeling of cell wall polysaccharides, namely chitin, β-(1,3)-glucan and mixed-linkage glucan synthases, as well as endo/exoglucanases and a mitogen-activated protein kinase. We have also analyzed the differential expression of these genes in D. teres spores and in the mycelium after cultivation on different media, as well as in the presence of Paraburkholderia phytofirmans strain PsJN, a plant growth-promoting bacterium (PGPB). The targeted gene expression analysis shows higher gene expression in the spores and in the mycelium with the application of PGPB. Besides analyzing key cell-wall-related genes, this study also identifies the most suitable reference genes to normalize qPCR results in D. teres, thus serving as a basis for future molecular studies on this ascomycete.

About 15,000 angiosperms are dioecious, but the mechanisms of sex determination in plants remain poorly understood. In particular, how Y chromosomes evolve and degenerate, and whether dosage compensation evolves as a response, are matters of debate. Here, we focus on Coccinia grandis, a dioecious cucurbit with the highest level of X/Y heteromorphy recorded so far. We identified sex-linked genes using RNA sequences from a cross and a model-based method termed SEX-DETector. Parents and F1 individuals were genotyped, and the transmission patterns of SNPs were then analyzed. In the >1300 sex-linked genes studied, maximum X-Y divergence was 0.13-0.17, and substantial Y degeneration is implied by an average Y/X expression ratio of 0.63 and an inferred gene loss on the Y of ~40%. We also found reduced Y gene expression being compensated by elevated expression of corresponding genes on the X and an excess of sex-biased genes on the sex chromosomes. Molecular evolution of sex-linked genes in C. grandis is thus comparable to that in Silene latifolia, another dioecious plant with a strongly heteromorphic XY system, and cucurbits are the fourth plant family in which dosage compensation is described, suggesting it might be common in plants.

Quercus acutissima Carruth. is a Chinese important energy plant with high ecological and economic values. While the species chloroplast genome has been reported, its mitochondrial genome (mitogenome) is still unexplored. Here, we assembled and annotated the Q. acutissima mitogenome, and we compared its characteristic differences with several closely related species. The Q. acutissima mitogenome's main structure is branched with three distinguished contigs (linear molecule 1, circular molecule 2, and circular molecule 3) with 448,982 bp total length and 45.72% GC content. The mitogenome contained 51 genes, including 32 protein-coding, 16 tRNA and 3 rRNA genes. We examined codon usage, repeated sequences, genome recombination, chloroplast to mitochondrion DNA transformation, RNA editing, and synteny in the Q. acutissima mitogenome. Phylogenetic trees based on 29 species mitogenomes clarified the species classification. Our results provided comprehensive information of Q. acutissima mitogenome, and they are expected to provide valuable information for Fagaceae evolutionary biology and to promote the species germplasm utilization.

The thyroid gland is an important endocrine organ modulating development, growth, and metabolism, mainly by controlling the synthesis and secretion of thyroid hormones (THs). However, little is known about the pig thyroid transcriptome. Long non-coding RNAs (lncRNAs) regulate gene expression and play critical roles in many cellular processes. Yorkshire pigs have a higher growth rate but lower fat deposition than that of Jinhua pigs, and thus, these species are ideal models for studying growth and lipid metabolism. This study revealed higher levels of THs in the serum of Yorkshire pigs than in the serum of Jinhua pigs. By using Ribo-zero RNA sequencing-which can capture both polyA and non-polyA transcripts-the thyroid transcriptome of both breeds were analyzed and 22,435 known mRNAs were found to be expressed in the pig thyroid. In addition, 1189 novel mRNAs and 1018 candidate lncRNA transcripts were detected. Multiple TH-synthesis-related genes were identified among the 455 differentially-expressed known mRNAs, 37 novel mRNAs, and 52 lncRNA transcripts. Bioinformatics analysis revealed that differentially-expressed genes were enriched in the microtubule-based process, which contributes to THs secretion. Moreover, integrating analysis predicted 13 potential lncRNA-mRNA gene pairs. These data expanded the repertoire of porcine lncRNAs and mRNAs and contribute to understanding the possible molecular mechanisms involved in animal growth and lipid metabolism.