WEGO 2.0: a web tool for analyzing and plotting GO annotations, 2018 update

Neural progenitor cells derived from fibroblasts induced by small molecule compounds under hypoxia for treatment of Parkinson's disease in rats

Neural progenitor cells (NPCs) capable of self-renewal and differentiation into neural cell lineages offer broad prospects for cell therapy for neurodegenerative diseases. However, cell therapy based on NPC transplantation is limited by the inability to acquire sufficient quantities of NPCs. Previous studies have found that a chemical cocktail of valproic acid, CHIR99021, and Repsox (VCR) promotes mouse fibroblasts to differentiate into NPCs under hypoxic conditions. Therefore, we used VCR (0.5 mM valproic acid, 3 μM CHIR99021, and 1 μM Repsox) to induce the reprogramming of rat embryonic fibroblasts into NPCs under a hypoxic condition (5%). These NPCs exhibited typical neurosphere-like structures that can express NPC markers, such as Nestin, SRY-box transcription factor 2, and paired box 6 (Pax6), and could also differentiate into multiple types of functional neurons and astrocytes in vitro. They had similar gene expression profiles to those of rat brain-derived neural stem cells. Subsequently, the chemically-induced NPCs (ciNPCs) were stereotactically transplanted into the substantia nigra of 6-hydroxydopamine-lesioned parkinsonian rats. We found that the ciNPCs exhibited long-term survival, migrated long distances, and differentiated into multiple types of functional neurons and glial cells in vivo. Moreover, the parkinsonian behavioral defects of the parkinsonian model rats grafted with ciNPCs showed remarkable functional recovery. These findings suggest that rat fibroblasts can be directly transformed into NPCs using a chemical cocktail of VCR without introducing exogenous factors, which may be an attractive donor material for transplantation therapy for Parkinson's disease.

Transcriptional memory of gene expression across generations participates in transgenerational plasticity of field pennycress in response to cadmium stress

Transgenerational plasticity (TGP) occurs when maternal environments influence the expression of traits in offspring, and in some cases may increase fitness of offspring and have evolutionary significance. However, little is known about the extent of maternal environment influence on gene expression of offspring, and its relationship with trait variations across generations. In this study, we examined TGP in the traits and gene expression of field pennycress (Thlaspi arvense) in response to cadmium (Cd) stress. In the first generation, along with the increase of soil Cd concentration, the total biomass, individual height, and number of seeds significantly decreased, whereas time to flowering, superoxide dismutase (SOD) activity, and content of reduced glutathione significantly increased. Among these traits, only SOD activity showed a significant effect of TGP; the offspring of Cd-treated individuals maintained high SOD activity in the absence of Cd stress. According to the results of RNA sequencing and bioinformatic analysis, 10,028 transcripts were identified as Cd-responsive genes. Among them, only 401 were identified as transcriptional memory genes (TMGs) that maintained the same expression pattern under normal conditions in the second generation as in Cd-treated parents in the first generation. These genes mainly participated in Cd tolerance-related processes such as response to oxidative stress, cell wall biogenesis, and the abscisic acid signaling pathways. The results of weighted correlation network analysis showed that modules correlated with SOD activity recruited more TMGs than modules correlated with other traits. The SOD-coding gene CSD2 was found in one of the modules correlated with SOD activity. Furthermore, several TMGs co-expressed with CSD2 were hub genes that were highly connected to other nodes and critical to the network's topology; therefore, recruitment of TMGs in offspring was potentially related to TGP. These findings indicated that, across generations, transcriptional memory of gene expression played an important role in TGP. Moreover, these results provided new insights into the trait evolution processes mediated by phenotypic plasticity.

Genome-wide identification and development of miniature inverted-repeat transposable elements and intron length polymorphic markers in tea plant (Camellia sinensis)

Marker-assisted breeding and tagging of important quantitative trait loci for beneficial traits are two important strategies for the genetic improvement of plants. However, the scarcity of diverse and informative genetic markers covering the entire tea genome limits our ability to achieve such goals. In the present study, we used a comparative genomic approach to mine the tea genomes of Camellia sinensis var. assamica (CSA) and C. sinensis var. sinensis (CSS) to identify the markers to differentiate tea genotypes. In our study, 43 and 60 Camellia sinensis miniature inverted-repeat transposable element (CsMITE) families were identified in these two sequenced tea genomes, with 23,170 and 37,958 putative CsMITE sequences, respectively. In addition, we identified 4912 non-redundant, Camellia sinensis intron length polymorphic (CsILP) markers, 85.8% of which were shared by both the CSS and CSA genomes. To validate, a subset of randomly chosen 10 CsMITE markers and 15 CsILP markers were tested and found to be polymorphic among the 36 highly diverse tea genotypes. These genome-wide markers, which were identified for the first time in tea plants, will be a valuable resource for genetic diversity analysis as well as marker-assisted breeding of tea genotypes for quality improvement.

Genome analysis of Phrixothrix hirtus (Phengodidae) railroad worm shows the expansion of odorant-binding gene families and positive selection on morphogenesis and sex determination genes

Among bioluminescent beetles of the Elateroidea superfamily, Phengodidae is the third largest family, with 244 bioluminescent species distributed only in the Americas, but is still the least studied from the phylogenetic and evolutionary points of view. The railroad worm Phrixothrix hirtus is an essential biological model and symbolic species due to its bicolor bioluminescence, being the only organism that produces true red light among bioluminescent terrestrial species. Here, we performed partial genome assembly of P. hirtus, combining short and long reads generated with Illumina sequencing, providing the first source of genomic information and a framework for comparative analyses of the bioluminescent system in Elateroidea. This is the largest genome described in the Elateroidea superfamily, with an estimated size of ∼3.4 Gb, displaying 32 % GC content, and 67 % transposable elements. Comparative genomic analyses showed a positive selection of genes and gene family expansion events of growths and morphogenesis gene products, which could be associated with the atypical anatomical development and morphogenesis found in paedomorphic females and underdeveloped males. We also observed gene family expansion among distinct odorant-binding receptors, which could be associated with the pheromone communication system typical of these beetles, and retrotransposable elements. Common genes putatively regulating bioluminescence production and control, including two luciferase genes corresponding to lateral lanterns green-emitting and head lanterns red-emitting luciferases with 7 exons and 6 introns, and genes potentially involved in luciferin biosynthesis were found, indicating that there are no clear differences about the presence or absence of gene families associated with bioluminescence in Elateroidea.

Transcriptome sequencing and DEG analysis in different developmental stages of floral buds induced by potassium chlorate in Dimocarpus longan

Potassium chlorate can promote off-season flowering in longan, but the molecular mechanisms are poorly understood. In this study, four-year-old 'Shixia' longan trees were injected in the trunk with potassium chlorate, and terminal buds were sampled and analyzed using transcriptomics and bioinformatics tools. To generate a reference longan transcriptome, we obtained 207,734 paired-end reads covering a total of 58,514,149 bp, which we assembled into 114,445 unigenes. Using this resource, we identified 3,265 differentially expressed genes (DEGs) that were regulated in longan terminal buds in response to potassium chlorate treatment for 2, 6 or 30 days, including 179 transcription factor genes. By reference to the Arabidopsis literature, we then defined 38 longan genes involved in flowering, from which we constructed the longan flowering pathway. According to RNA-seq data, at least 24 of these genes, which participate in multiple signaling pathways, are involved in potassium chlorate-stimulated floral induction, and the differential regulation in terminal buds of ten floral pathway genes (GI, CO, GID1, GA4, GA5, FLC, AP1, LFY, FT and SOC1) was confirmed by qRT-PCR. These data will contribute to an improved understanding of the functions of key genes involved in longan floral induction by potassium chlorate.

Transcriptome analysis of the response of Hypomyces chrysospermus to cadmium stress

Hypomyces chrysospermus is a fungal parasite that grows on Boletus species. One isolated strain of H. chrysospermus from B. griseus was obtained and proved of strong ability to tolerate and absorb cadmium (Cd) by previous research. However, the molecular mechanisms of underlying the resistance of H. chrysospermus to Cd stress have not been investigated. This study aimed to assess the effect of Cd stress on the global transcriptional regulation of H. chrysospermus. A total of 1,839 differentially expressed genes (DEGs) were identified under 120 mg/l Cd stress. Gene ontology (GO) enrichment analysis revealed that large amounts of DEGs were associated with cell membrane components, oxidoreductase activity, and transport activity. KEGG enrichment analysis revealed that these DEGs were mainly involved in the translation, amino acid metabolism, transport and catabolism, carbohydrate metabolism, and folding/sorting and degradation pathways under Cd stress. Moreover, the expression of DEGs encoding transporter proteins, antioxidant enzymes, nonenzymatic antioxidant proteins, detoxification enzymes, and transcription factors was associated with the Cd stress response. These results provide insights into the molecular mechanisms underlying Cd tolerance in H. chrysospermus and serve as a valuable reference for further studies on the detoxification mechanisms of heavy metal-tolerant fungi. Our findings may also facilitate the development of new and improved fungal bioremediation strategies.

An Integrative Transcriptional Network Revealed Spatial Molecular Interplay Underlying Alantolactone and Inulin Biosynthesis in Inula racemosa Hook f

Inula racemosa Hook. f. (Pushkarmula), a perennial Himalayan herb known for its aromatic and phytopharmaceutical attributes, is not yet explored at genomic/transcriptomic scale. In this study, efforts were made to unveil the global transcriptional atlas underlying organ-specific specialized metabolite biosynthesis by integrating RNA-Seq analysis of 433 million sequenced reads with the phytochemical analysis of leaf, stem, and root tissues. Overall, 7242 of 83,772 assembled nonredundant unigenes were identified exhibiting spatial expression in leaf (3761), root (2748), and stem (733). Subsequently, integration of the predicted transcriptional interactome network of 2541 unigenes (71,841 edges) with gene ontology and KEGG pathway enrichment analysis revealed isoprenoid, terpenoid, diterpenoid, and gibberellin biosynthesis with antimicrobial activities in root tissue. Interestingly, the root-specific expression of germacrene-mediated alantolactone biosynthesis (GAS, GAO, G8H, IPP, DMAP, and KAO) and antimicrobial activities (BZR1, DEFL, LTP) well-supported with both quantitative expression profiling and phytochemical accumulation of alantolactones (726.08 μg/10 mg) and isoalantolactones (988.59 μg/10 mg), which suggests “roots” as the site of alantolactone biosynthesis. A significant interaction of leaf-specific carbohydrate metabolism with root-specific inulin biosynthesis indicates source (leaf) to sink (root) regulation of inulin. Our findings comprehensively demonstrate the source-sink transcriptional regulation of alantolactone and inulin biosynthesis, which can be further extended for upscaling the targeted specialized metabolites. Nevertheless, the genomic resource created in this study can also be utilized for development of genome-wide functionally relevant molecular markers to expedite the breeding strategies for genetic improvement of I. racemosa.

Identifying Soybean Pod Borer (Leguminivora glycinivorella) Resistance QTLs and the Mechanism of Induced Defense Using Linkage Mapping and RNA-Seq Analysis

The soybean pod borer (Leguminivora glycinivorella) (SPB) is a major cause of soybean (Glycine max L.) yield losses in northeast Asia, thus it is desirable to elucidate the resistance mechanisms involved in soybean response to the SPB. However, few studies have mapped SPB-resistant quantitative trait loci (QTLs) and deciphered the response mechanism in soybean. Here, we selected two soybean varieties, JY93 (SPB-resistant) and K6 (SPB-sensitive), to construct F₂ and F_2:3 populations for QTL mapping and collected pod shells before and after SPB larvae chewed on the two parents to perform RNA-Seq, which can identify stable QTLs and explore the response mechanism of soybean to the SPB. The results show that four QTLs underlying SPB damage to seeds were detected on chromosomes 4, 9, 13, and 15. Among them, qESP-9-1 was scanned in all environments, hence it can be considered a stable QTL. All QTLs explained 0.79 to 6.09% of the phenotypic variation. Meanwhile, 2298 and 3509 DEGs were identified for JY93 and K6, respectively, after the SPB attack, and most of these genes were upregulated. Gene Ontology enrichment results indicated that the SPB-induced and differently expressed genes in both parents are involved in biological processes such as wound response, signal transduction, immune response, and phytohormone pathways. Interestingly, secondary metabolic processes such as flavonoid synthesis were only significantly enriched in the upregulated genes of JY93 after SPB chewing compared with K6. Finally, we identified 18 candidate genes related to soybean pod borer resistance through the integration of QTL mapping and RNA-Seq analysis. Seven of these genes had similar expression patterns to the mapping parents in four additional soybean germplasm after feeding by the SPB. These results provide additional knowledge of the early response and induced defense mechanisms against the SPB in soybean, which could help in breeding SPB-resistant soybean accessions.

Study of the migration of Fasciola hepatica juveniles across the intestinal barrier of the host by quantitative proteomics in an ex vivo model

Fasciola hepatica is a trematode parasite that infects animals and humans causing fasciolosis, a worldwide-distributed disease responsible for important economic losses and health problems. This disease is of growing public health concern since parasite isolates resistant to the current treatment (triclabendazole) have increasingly been described. F. hepatica infects its vertebrate host after ingestion of the encysted parasite (metacercariae), which are found in the water or attached to plants. Upon ingestion, newly excysted juveniles of F. hepatica (FhNEJ) emerge in the intestinal lumen and cross the intestinal barrier, reach the peritoneum and migrate to the biliary ducts, where adult worms fully develop. Despite the efforts made to develop new therapeutic and preventive tools, to date, protection against F. hepatica obtained in different animal models is far from optimal. Early events of host-FhNEJ interactions are of paramount importance for the infection progress in fasciolosis, especially those occurring at the host-parasite interface. Nevertheless, studies of FhNEJ responses to the changing host environment encountered during migration across host tissues are still scarce. Here, we set-up an ex vivo model coupled with quantitative SWATH-MS proteomics to study early host-parasite interaction events in fasciolosis. After comparing tegument and somatic fractions from control parasites and FhNEJ that managed to cross a mouse intestinal section ex vivo, a set of parasite proteins whose expression was statistically different were found. These included upregulation of cathepsins L3 and L4, proteolytic inhibitor Fh serpin 2, and a number of molecules linked with nutrient uptake and metabolism, including histone H4, H2A and H2B, low density lipoprotein receptor, tetraspanin, fatty acid binding protein a and glutathione-S-transferase. Downregulated proteins in FhNEJ after gut passage were more numerous than the upregulated ones, and included the heath shock proteins HSP90 and alpha crystallin, amongst others. This study brings new insights into early host-parasite interactions in fasciolosis and sheds light on the proteomic changes in FhNEJ triggered upon excystment and intestinal wall crossing, which could serve to define new targets for the prevention and treatment of this widespread parasitic disease.

Fasciolosis caused by the helminth parasite Fasciola hepatica is a serious health and economic problem worldwide. Treatment and prevention of this disease pose several drawbacks that have so far not been solved. The definition of suitable parasite molecular targets to overcome such drawbacks should be based on thoroughly deciphering host-parasite interactions, and in this regard most studies have focused on the adult stages of F. hepatica. Nevertheless, in this context, the study of the transient juvenile stages of this parasite could be of higher utility due to the importance of early interactions with the host for parasite migration and the successful establishment of infection. In this work, we set-up an ex vivo model and performed a quantitative proteomics approach to study the changes in F. hepatica juveniles upon gut passage. We found that the parasite tegument and somatic compartments experienced deep changes in their composition and showed that the host triggers the expression of specific molecules that are important for parasite migration and survival at this stage. The molecules described here could serve to better understand host-parasite interactions and to define new targets against fasciolosis.

Proteomic analysis response of rice (Oryza sativa) leaves to ultraviolet-B radiation stress

Rice (Oryza sativa) is a human staple food and serves as a model organism for genetic and molecular studies. Few studies have been conducted to determine the effects of ultraviolet-B (UV-B) stress on rice. UV-B stress triggers morphological and physiological changes in plants. However, the underlying mechanisms governing these integrated responses are unknown. In this study, we conducted a proteomic response of rice leaves to UV-B stress using two-dimensional gel electrophoresis and identified the selected proteins by mass spectrometry analysis. Four levels of daily biologically effective UV-B radiation intensities were imposed to determine changes in protein accumulation in response to UV-B stress: 0 (control), 5, 10, and 15 kJ m^-2 d^-1in two cultivars, i.e., IR6 and REX. To mimic the natural environment, we conducted this experiment in Sunlit Soil-Plant-Atmosphere-Research (SPAR) chambers. Among the identified proteins, 11% of differentially expressed proteins were found in both cultivars. In the Rex cultivar, only 45% of proteins are differentially expressed, while only 27.5% were expressed in IR6. The results indicate that REX is more affected by UV-B stress than IR6 cultivars. The identified protein TSJT1 (spot 16) in both cultivars plays a crucial role in plant growth and development during stress treatment. Additionally, we found that UV-B stress altered many antioxidant enzymes associated with redox homeostasis and cell defense response. Another enzyme, the glyceraldehyde-3-phosphate dehydrogenase (GAPDH), has been identified as spot 15, which plays an essential role in glycolysis and cellular energy production. Another vital protein identified is glycosyl hydrolase (GH) as spot 9, which catalyzes the hydrolysis of glycosidic bonds in cell wall polymers and significantly affects cell wall architecture. Some identified proteins are related to photosynthesis, protein biosynthesis, signal transduction, and stress response. The findings of our study provide new insights into understanding how rice plants are tailored to UV-B stress via modulating the expression of UV-B responsive proteins, which will help develop superior rice breeds in the future to combat UV-B stress. Data are available via ProteomeXchange with identifier PXD032163.

MSD2-mediated ROS metabolism fine-tunes the timing of floral organ abscission in Arabidopsis

The timely removal of end-of-purpose flowering organs is as essential for reproduction and plant survival as timely flowering. Despite much progress in understanding the molecular mechanisms of floral organ abscission, little is known about how various environmental factors are integrated into developmental programmes that determine the timing of abscission. Here, we investigated whether reactive oxygen species (ROS), mediators of various stress-related signalling pathways, are involved in determining the timing of abscission and, if so, how they are integrated with the developmental pathway in Arabidopsis thaliana. MSD2, encoding a secretory manganese superoxide dismutase, was preferentially expressed in the abscission zone of flowers, and floral organ abscission was accelerated by the accumulation of ROS in msd2 mutants. The expression of the genes encoding the receptor-like kinase HAESA (HAE) and its cognate peptide ligand INFLORESCENCE DEFICIENT IN ABSCISSION (IDA), the key signalling components of abscission, was accelerated in msd2 mutants, suggesting that MSD2 acts upstream of IDA-HAE. Further transcriptome and pharmacological analyses revealed that abscisic acid and nitric oxide facilitate abscission by regulating the expression of IDA and HAE during MSD2-mediated signalling. These results suggest that MSD2-dependent ROS metabolism is an important regulatory point integrating environmental stimuli into the developmental programme leading to abscission.

Differential gene expression and identification of growth-related genes in the pituitary gland of South African goats

Growth and carcass quality are economically important traits in goat production. This study investigated differentially expressed genes from the caprine pituitary gland transcriptome of South African indigenous goat breeds of varying growth performances and carcass quality parameters. Tissues were harvested from the pituitary gland of three South African Boer goats and three village ecotype goats all raised under similar conditions simulating intensive commercial production systems. Three additional tissues were harvested from village ecotype goats that were raised extensively on village farms. Between breed differences were investigated by comparing differential gene expression among three South African Boer and three village goats that were both raised under intensive commercial production system at a research farm. Within-breed differences were investigated by comparing differential gene expression among three village goats raised under extensive conditions (on-farm in Pella, S.A. village farming community) and three village goats raised under intensive commercial production system (at ARC research farm in Pretoria, South Africa. Total RNA was isolated from the pituitary gland of 36-week-old animals (n = 9) and sequenced individually in triplicates. An average of 28,298,512 trimmed, and quality-controlled reads/animal were mapped to the goat genome (Capra_hircus.ARS1.94) using HiSat2 software. Transcript assembly and quantification yielded 104 differentially expressed genes for village goats raised under extensive system and 62 for village goats raised under the intensive production system at the false discovery rate (FRD) of ≤0.05 and a fold change of ≥2. Growth-related genes such as POU3F4 and TSHZ1 were highly expressed within breeds raised under both production systems. Conversely, growth-related genes such as FGFR2 and SMPX genes were highly expressed between breeds raised under similar production systems. Ballgown analysis revealed a high expression of GH1 and IGF1 in the intensively raised compared to extensively raised goats. Both genes were also highly expressed in the village goats when compared to the Boer. The differential gene expression data provided insights into genes and molecular mechanisms associated with growth and growth development in goats.

Proteomic analysis of Malaysian Horseshoe crab (Tachypleus gigas) hemocytes gives insights into its innate immunity host defence system and other biological processes

Horseshoe crabs are one of the most studied invertebrates due to their remarkable innate immunity mechanism and biological processes. In this work, the proteins of the lipopolysaccharides (LPS)-stimulated and non-stimulated hemocytes of Malaysian Tachypleus gigas were profiled using LC-MS/MS. A total of 154 proteins were identified in both types of samples. Additionally, seventy-seven proteins were commonly found in both conditions, while 52 and 25 proteins were uniquely found in the LPS-stimulated and non-stimulated hemocytes, respectively. ATP-dependent energy-generating proteins such as actins and BLTX actin-related proteins were detected in both stimulated and non-stimulated T. gigas hemocytes, but more of such proteins were found in the former type. Proteins such as tachylectin-2, coagulogen, c-reactive proteins, histones, hemocyanin, and DNA polymerase, which play key roles in the organism’s innate immunity, were differentially expressed in the hemocytes following LPS challenge. In conclusion, the proteins identified in the hemolymph of T. gigas are vital for the organism’s molecular functions, biological processes, and activation of innate immunity.

Comparative Transcriptome Analysis of Deep-Rooting and Shallow-Rooting Potato (Solanum tuberosum L.) Genotypes under Drought Stress

The selection and breeding of deep rooting and drought-tolerant varieties has become a promising approach for improving the yield and adaptability of potato (Solanum tuberosum L.) in arid and semiarid areas. Therefore, the discovery of root-development-related genes and drought tolerance signaling pathways in potato is important. In this study, we used deep-rooting (C119) and shallow-rooting (C16) potato genotypes, with different levels of drought tolerance, to achieve this objective. Both genotypes were treated with 150 mM mannitol for 0 h (T0), 2 h (T2), 6 h (T6), 12 h (T12), and 24 h (T24), and their root tissues were subjected to comparative transcriptome analysis. A total of 531, 1571, 1247, and 3540 differentially expressed genes (DEGs) in C16 and 1531, 1108, 674, and 4850 DEGs in C119 were identified in T2 vs. T0, T6 vs. T2, T12 vs. T6, and T24 vs. T12 comparisons, respectively. Gene expression analysis indicated that a delay in the onset of drought-induced transcriptional changes in C16 compared with C119. Functional enrichment analysis revealed genotype-specific biological processes involved in drought stress tolerance. The metabolic pathways of plant hormone transduction and MAPK signaling were heavily involved in the resistance of C16 and C119 to drought, while abscisic acid (ABA), ethylene, and salicylic acid signal transduction pathways likely played more important roles in C119 stress responses. Furthermore, genes involved in root cell elongation and division showed differential expression between the two genotypes under drought stress. Overall, this study provides important information for the marker-assisted selection and breeding of drought-tolerant potato genotypes.

Transcriptomic responses of the zearalenone (ZEN)-detoxifying yeast Apiotrichum mycotoxinivorans to ZEN exposure

Zearalenone (ZEN) is a potent oestrogenic mycotoxin that is mainly produced by Fusarium species and is a serious environmental pollutant in animal feeds. Apiotrichum mycotoxinivorans has been widely used as a feed additive to detoxify ZEN. However, the effects of ZEN on A. mycotoxinivorans and its detoxification mechanisms remain unclear. In this study, transcriptomic and bioinformatic analyses were used to investigate the molecular responses of A. mycotoxinivorans to ZEN exposure and the genetic basis of ZEN detoxification. We detected 1424 significantly differentially expressed genes (DEGs), of which 446 were upregulated and 978 were downregulated. Functional and enrichment analyses showed that ZEN-induced genes were significantly associated with xenobiotic metabolism, oxidative stress response, and active transport systems. However, ZEN-inhibited genes were mainly related to cell division, cell cycle, and fungal development. Subsequently, bioinformatic analysis identified candidate ZEN-detoxification enzymes. The Baeyer-Villiger monooxygenases and carboxylesterases, which are responsible for the formation and subsequent hydrolysis of a new ZEN lactone, respectively, were significantly upregulated. In addition, the expression levels of genes related to conjugation and transport involved in the xenobiotic detoxification pathway were significantly upregulated. Moreover, the expression levels of genes encoding enzymatic antioxidants and those related to growth and apoptosis were significantly upregulated and downregulated, respectively, which made it possible for A. mycotoxinivorans to survive in a highly toxic environment and efficiently detoxify ZEN. This is the first systematic report of ZEN tolerance and detoxification in A. mycotoxinivorans. We identified the metabolic enzymes that were potentially involved in detoxifying ZEN in the GMU1709 strain and found that ZEN-induced transcriptional regulation of genes is key to withstanding highly toxic environments. Hence, our results provide valuable information for developing enzymatic detoxification systems or engineering this detoxification pathway in other species.

Correlation of Moraxella catarrhalis macrolide susceptibility with the ability to adhere and invade human respiratory epithelial cells

Recently, the prevalence of macrolide-resistant Moraxella catarrhalis has been reported, especially among Chinese children. The fitness cost of resistance is reported to render the resistant bacteria less virulent. To investigate the correlation between macrolide susceptibility of M. catarrhalis and pathogenicity, the whole genome of 70 M. catarrhalis isolates belonging to four clonal complexes with different macrolide susceptibilities was sequenced. The gene products were annotated with the Gene Ontology terms. Based on 46 extracted essential virulence genes, 19 representative isolates were selected to infect type II alveolar cells (A549 cells). The ability of these isolates to adhere and invade human epithelial cells and to produce cytokines was comparatively analysed. Furthermore, mice were infected with a pair of M. catarrhalis isolates with different pathogenic behaviours and macrolide susceptibilities to examine pulmonary clearance, histological findings, and the production of cytokines. The percentages of annotations for binding, metabolic process, cellular process, and cell were non-significantly different between the macrolide-resistant and macrolide-susceptible groups. The presence of uspA2, uspA2H, pilO, lbpB, lex1, modM, mboIA, and mboIB significantly differed among the four clonal complexes and macrolide susceptibility groups. Furthermore, compared with those in macrolide-susceptible isolates, the adhesion ability was stronger (P = 0.0019) and the invasion ability was weaker (P < 0.0001) in the macrolide-resistant isolates. Mouse experiments revealed that pulmonary macrophages elicit immune responses against M. catarrhalis infection by significantly upregulating the Csf2, Il4, Il13, Il1b, Il6, Tnf, and Il18. Therefore, M. catarrhalis populations exhibited diverse pathogenicity in vitro and in vivo.

Bioinformatics Strategies to Identify Shared Molecular Biomarkers That Link Ischemic Stroke and Moyamoya Disease with Glioblastoma

Expanding data suggest that glioblastoma is accountable for the growing prevalence of various forms of stroke formation, such as ischemic stroke and moyamoya disease. However, the underlying deterministic details are still unspecified. Bioinformatics approaches are designed to investigate the relationships between two pathogens as well as fill this study void. Glioblastoma is a form of cancer that typically occurs in the brain or spinal cord and is highly destructive. A stroke occurs when a brain region starts to lose blood circulation and prevents functioning. Moyamoya disorder is a recurrent and recurring arterial disorder of the brain. To begin, adequate gene expression datasets on glioblastoma, ischemic stroke, and moyamoya disease were gathered from various repositories. Then, the association between glioblastoma, ischemic stroke, and moyamoya was established using the existing pipelines. The framework was developed as a generalized workflow to allow for the aggregation of transcriptomic gene expression across specific tissue; Gene Ontology (GO) and biological pathway, as well as the validation of such data, are carried out using enrichment studies such as protein–protein interaction and gold benchmark databases. The results contribute to a more profound knowledge of the disease mechanisms and unveil the projected correlations among the diseases.

Integrated metabolomic and transcriptomic analysis identifies benzo[a]pyrene-induced characteristic metabolic reprogramming during accumulation of lipids and reactive oxygen species in macrophages

Polycyclic aromatic hydrocarbon exposure is a major risk factor for cardiovascular diseases. Macrophage lipid accumulation is a characteristic molecular event in the pathophysiology of cardiovascular diseases. Metabolic reprogramming is an intervention target for diseases and toxic effects of environmental pollutants. However, comprehensive metabolic reprogramming related to BaP-induced macrophage lipid accumulation is currently unexplored. Therefore, metabolomics and transcriptomics were conducted to unveil relevant metabolic reprogramming in BaP-exposed macrophages, and to discover potential intervention targets. Metabolomics revealed that most amino acids, nucleotides, monosaccharides, and organic acids were significantly decreased, while most fatty acids and steroids accumulated in BaP-exposed macrophages. Transcriptomics showed that fatty acid synthesis and oxidation, and steroid synthesis and export were decreased, while import of fatty acids and steroids was increased, indicating potential roles of lipid transport in macrophage lipid accumulation following BaP exposure. Meanwhile, alanine, aspartate and glutamate metabolism, branched-chain amino acid degradation, nucleotide synthesis, monosaccharide import, pentose phosphate pathway, citrate synthesis, and glycolysis were decreased, while nucleotide degradation was increased, thus inducing decreases in most amino acids, nucleotides, monosaccharides, and organic acids in BaP-exposed macrophages. Additionally, increases in oxidative stress and the activation of antioxidant systems were observed in BaP-exposed macrophages, which was evinced by increases in reactive oxygen species, and the activation of Fenton reaction, Vdac2/3, Sod2, and Nrf2. Moreover, BaP-induced accumulation of reactive oxygen species and lipids in macrophages could be abolished by epigallocatechin-3-gallate. Quantitative PCR showed that BaP exposure activated aryl hydrocarbon receptor signaling and promoted the proinflammatory phenotype in macrophages, and these effects were inhibited or even abolished by the separate treatment with epigallocatechin-3-gallate or CH-223191, suggesting the regulatory role of aryl hydrocarbon receptor signaling in BaP-induced toxic effects. This study provides novel insights into the toxic effects of polycyclic aromatic hydrocarbons on macrophage metabolism and potential intervention targets.

Distribution Pattern of N6-Methyladenine DNA Modification in the Seashore Paspalum (Paspalum vaginatum) Genome

N6-methyladenine (6mA) DNA modification has been detected in several eukaryotic organisms, in some of them, it plays important role in the regulation process of stress-resistance response. However, the genome-wide distribution patterns and potential functions of 6mA DNA modification in halophyte Seashore paspalum (Paspalum vaginatum) remain largely unknown. Here, we examined the 6mA landscape in the P. vaginatum genome by adopting single molecule real-time sequencing technology and found that 6mA modification sites were broadly distributed across the P. vaginatum genome. We demonstrated distinct 6mA methylation levels and 6mA distribution patterns in different types of transcription genes, which hinted at different epigenetic rules. Furthermore, the moderate 6mA density genes in P. vaginatum functionally correlated with stress resistance, which also maintained a higher transcriptional level. On the other hand, a specific 6mA distribution pattern in the gene body and near TSS was observed in gene groups with higher RNA expression, which maybe implied some kind of regularity between 6mA site distribution and the protein coding genes transcription was possible. Our study provides new insights into the association between 6mA methylation and gene expression, which may also contribute to key agronomic traits in P. vaginatum.

Transcriptome analysis reveals putative pathogenesis genes in Alternaria panax during infecting Panax notoginseng leaves

BACKGROUND

Alternaria panax is the causative agent of black spot disease in Panax notoginseng, which causes significant yield loss. However, the molecular mechanisms underlying its pathogenicity remain mostly unknown.

OBJECTIVE

We sequenced the transcriptome of A. panax during infecting P. notoginseng leaves using next-generation RNA-seq to understand the molecular aspects of black spot disease.

METHODS

In this study, we sequenced the A. panax transcriptome during infecting P. notoginseng leaves through next-generation sequencing to explore the pathogenesis genes that may be responsible for black spot disease on P. notoginseng.

RESULT

The de novo transcriptome assembly of A. panax produced 23,036 unigenes, of which 18,096 genes were functionally annotated by at least one protein database. GO enrichment analysis and KEGG pathways of differentially up-regulated genes suggest that most genes are associated with metabolic processes, catalytic activity, starch, and sucrose metabolism during infection. Many pathogenesis-associated genes, including genes encoding secreted proteins, candidate secreted effectors, cell wall degrading enzymes, transcription factors, and transporters, were up-regulated in A. panax during infection. In addition, the secondary metabolite biosynthesis genes, including cytochrome P450, and nonribosomal peptide synthetases, were also identified in this study.

CONCLUSIONS

Differential gene expression analysis has confirmed that A. panax infection was mainly present in the middle and final stages. The findings show that these pathogenesis-associated genes in A. panax may be critical for the P. notoginseng black spots disease.

Comparative Genomic Analysis of SAUR Gene Family, Cloning and Functional Characterization of Two Genes (PbrSAUR13 and PbrSAUR52) in Pyrus bretschneideri

The SAUR (small auxin-up RNA) gene family is the biggest family of early auxin response genes in higher plants and has been associated with the control of a variety of biological processes. Although SAUR genes had been identified in several genomes, no systematic analysis of the SAUR gene family has been reported in Chinese white pear. In this study, comparative and systematic genomic analysis has been performed in the SAUR gene family and identified a total of 116 genes from the Chinese white pear. A phylogeny analysis revealed that the SAUR family could be classified into four groups. Further analysis of gene structure (introns/exons) and conserved motifs showed that they are diverse functions and SAUR-specific domains. The most frequent mechanisms are whole-genome duplication (WGD) and dispersed duplication (DSD), both of which may be important in the growth of the SAUR gene family in Chinese white pear. Moreover, cis-acting elements of the PbrSAUR genes were found in promoter regions associated with the auxin-responsive elements that existed in most of the upstream sequences. Remarkably, the qRT-PCR and transcriptomic data indicated that PbrSAUR13 and PbrSAUR52 were significantly expressed in fruit ripening. Subsequently, subcellular localization experiments revealed that PbrSAUR13 and PbrSAUR52 were localized in the nucleus. Moreover, PbrSAUR13 and PbrSAUR52 were screened for functional verification, and Dangshan pear and frandi strawberry were transiently transformed. Finally, the effects of these two genes on stone cells and lignin were analyzed by phloroglucinol staining, Fourier infrared spectroscopy, and qRT-PCR. It was found that PbrSAUR13 promoted the synthesis and accumulation of stone cells and lignin, PbrSAUR52 inhibited the synthesis and accumulation of stone cells and lignin. In conclusion, these results indicate that PbrSAUR13 and PbrSAUR52 are predominantly responsible for lignin inhibit synthesis, which provides a basic mechanism for further study of PbrSAUR gene functions.

Conversion between duplicated genes generated by polyploidization contributes to the divergence of poplar and willow

BACKGROUND

Gene conversion has an important effect on duplicate genes produced by polyploidization. Poplar (Populus trichocarpa) and willow (Salix brachista) are leading models and excellent green plants in the Salicaceae. Although much attention has been paid to the evolution of duplicated genes in poplar and willow, the role of conversion between duplicates generated from polyploidization remains poorly understood.

RESULTS

Here, through genomic synteny analyses, we identified duplicate genes generated by the Salicaceae common tetraploidization (SCT) in the poplar and willow genomes. We estimated that at least 0.58% and 0.25% of poplar and willow duplicates were affected by whole-gene conversion after the poplar-willow divergence, with more (5.73% and 2.66%) affected by partial-gene conversion. Moreover, we found that the converted duplicated genes were unevenly distributed on each chromosome in the two genomes, and the well-preserved homoeologous chromosome regions may facilitate the conversion of duplicates. Notably, we found that conversion maintained the similarity of duplicates, likely contributing to the conservation of certain sequences, but is essentially accelerated the rate of evolution and increased species divergence. In addition, we found that converted duplicates tended to have more similar expression patterns than nonconverted duplicates. We found that genes associated with multigene families were preferentially converted. We also found that the genes encoding conserved structural domains associated with specific traits exhibited a high frequency of conversion.

CONCLUSIONS

Extensive conversion between duplicate genes generated from the SCT contributes to the diversification of the family Salicaceae and has had long-lasting effects on those genes with important biological functions.

A Comprehensive Study of the WRKY Transcription Factor Family in Strawberry

WRKY transcription factors play critical roles in plant growth and development or stress responses. Using up-to-date genomic data, a total of 64 and 257 WRKY genes have been identified in the diploid woodland strawberry, Fragaria vesca, and the more complex allo-octoploid commercial strawberry, Fragaria × ananassa cv. Camarosa, respectively. The completeness of the new genomes and annotations has enabled us to perform a more detailed evolutionary and functional study of the strawberry WRKY family members, particularly in the case of the cultivated hybrid, in which homoeologous and paralogous FaWRKY genes have been characterized. Analysis of the available expression profiles has revealed that many strawberry WRKY genes show preferential or tissue-specific expression. Furthermore, significant differential expression of several FaWRKY genes has been clearly detected in fruit receptacles and achenes during the ripening process and pathogen challenged, supporting a precise functional role of these strawberry genes in such processes. Further, an extensive analysis of predicted development, stress and hormone-responsive cis-acting elements in the strawberry WRKY family is shown. Our results provide a deeper and more comprehensive knowledge of the WRKY gene family in strawberry.

Comparative transcriptome and tissue-specific expression analysis of genes reveal tissue-cultured plants as an alternative source for phenylethanoids and phenylpropanoids in Rhodiola imbricata (Edgew.)

Rhodiola imbricata (Crassulaceae) is a traditional trans-Himalayan endangered medicinal herb with immense therapeutic applications. Over the years, over-exploitation, un-managed harvesting, and lack of captive cultivation procedures persuaded threat to its wild habitat. Plant tissue culture and RNA-Seq-based molecular bioprospection of key regulatory genes aid the understanding of molecular dynamics involved in specialized metabolites (phenylethanoids and phenylpropanoids) biosynthesis and its sustainable production. Hence, comparative transcriptomic analysis was performed using leaf and root tissues from the wild and tissue-cultured plants, revealing tissue-specific production of salidroside and rosavin. The transcriptome profiling resulted in 345 million high-quality reads yielding 92,380 unique transcripts with an N₅₀ of 1260 bp. Tissue-specific gene expression analysis revealed that both phenylethanoids and phenylpropanoids biosynthesis are predominantly associated with the shikimate pathway. In addition to RNA-Seq data, the downstream biosynthesis pathways genes viz., phospho-2-dehydro-3-deoxyheptonate aldolase (DAHPS), 3-dehydroquinate synthase (DHQS), shikimate kinase (SK), chorismate mutase (CM), arogenate dehydrogenase (TYRAAT), aromatic-L-amino-acid decarboxylase (TDC), phenylalanine ammonia-lyase (PAL), 4-coumarate-CoA ligase (4-CL), cinnamoyl-CoA reductase (CCR), and cinnamyl alcohol dehydrogenase (CAD) showed higher expression pattern in wild plant tissues compared to tissue-cultured plants. The transcript fold expression determined by RT-qPCR results followed similar patterns as those observed in RNA-seq and targeted metabolite profiling data. Salidroside and rosavin content in wild plants exhibited 2.40 fold and 1.77 fold increase accumulation compared to the tissue-cultured plant. The present investigation explained the tissue and condition-specific significant differences between the expression of proposed biosynthetic pathway genes and salidroside and rosavin content. Additionally, NAC, bHLH, and ARF were the most abundant transcription factor families found in the transcriptomic analysis of R. imbricata. The generated transcriptome dataset provides a valuable gene(s)/transcription factors hub that can be used for the sustainable production of salidroside and rosavin in R. imbricata under tissue culture conditions.

A high-quality genome assembly and annotation of the dark-eyed junco Junco hyemalis, a recently diversified songbird

The dark-eyed junco (Junco hyemalis) is one of the most common passerines of North America, and has served as a model organism in studies related to ecophysiology, behavior, and evolutionary biology for over a century. It is composed of at least 6 distinct, geographically structured forms of recent evolutionary origin, presenting remarkable variation in phenotypic traits, migratory behavior, and habitat. Here, we report a high-quality genome assembly and annotation of the dark-eyed junco generated using a combination of shotgun libraries and proximity ligation Chicago and Dovetail Hi-C libraries. The final assembly is ∼1.03 Gb in size, with 98.3% of the sequence located in 30 full or nearly full chromosome scaffolds, and with a N50/L50 of 71.3 Mb/5 scaffolds. We identified 19,026 functional genes combining gene prediction and similarity approaches, of which 15,967 were associated to GO terms. The genome assembly and the set of annotated genes yielded 95.4% and 96.2% completeness scores, respectively when compared with the BUSCO avian dataset. This new assembly for J. hyemalis provides a valuable resource for genome evolution analysis, and for identifying functional genes involved in adaptive processes and speciation.

Identification and Characterization of Chemosensory Receptors in the Pheromone Gland-Ovipositor of Spodoptera frugiperda (J. E. Smith)

Simple Summary

Chemical cues are generally thought to be primarily detected by the cephalic organ antennae, maxillary palps, and proboscises in insects. Although several recent studies have reported the chemosensory roles of ovipositors in some moth species, the expression of chemosensory receptors and their functions in the ovipositor remain largely unknown. Here, we systematically analyzed the pheromone gland-ovipositor (PG-OV) transcriptome of the fall armyworm, Spodoptera frugiperda (Lepidoptera: Noctuidae). A total of 26 candidate chemosensory receptor genes were revealed, including 12 odorant receptors (ORs), 4 gustatory receptors (GRs), and 10 ionotropic receptors (IRs). Specific genes including pheromone receptors, ORco, CO₂ receptors, sugar receptors, and IR co-receptors were identified. Tissue expression profiling demonstrated that the annotated receptor genes were mainly expressed in the antennae (for ORs and IRs) or proboscis (for GRs), but two ORs, two GRs, and two IRs were also highly enriched in the PG-OV, with expression levels only slightly lower or even similar to those in the antennae/proboscis. This report provides the first large-scale description of chemosensory receptors in the PG-OV of S. frugiperda. It may inspire researchers to investigate how chemosensory receptors function in the ovipositor of S. frugiperda, as well as in the ovipositors of other moths.

Abstract

Chemoreception by moth ovipositors has long been suggested, but underlying molecular mechanisms are mostly unknown. To reveal such chemosensory systems in the current study, we sequenced and assembled the pheromone gland-ovipositor (PG-OV) transcriptome of females of the fall armyworm, Spodoptera frugiperda, a pest of many crops. We annotated a total of 26 candidate chemosensory receptor genes, including 12 odorant receptors (ORs), 4 gustatory receptors (GRs), and 10 ionotropic receptors (IRs). The relatedness of these chemosensory receptors with those from other insect species was predicted by phylogenetic analyses, and specific genes, including pheromone receptors, ORco, CO₂ receptors, sugar receptors, and IR co-receptors, were reported. Although real-time quantitative-PCR analyses of annotated genes revealed that OR and IR genes were mainly expressed in S. frugiperda antennae, two ORs and two IRs expressed in antennae were also highly expressed in the PG-OV. Similarly, GR genes were mainly expressed in the proboscis, but two were also highly expressed in the PG-OV. Our study provides the first large-scale description of chemosensory receptors in the PG-OV of S. frugiperda and provides a foundation for exploring the chemoreception mechanisms of PG-OV in S. frugiperda and in other moth species.

Identification of microRNAs in the Lyme Disease Vector Ixodes scapularis

MicroRNAs (miRNAs) are a class of small non-coding RNAs involved in many biological processes, including the immune pathways that control bacterial, parasitic, and viral infections. Pathogens probably modify host miRNAs to facilitate successful infection, so they might be useful targets for vaccination strategies. There are few data on differentially expressed miRNAs in the black-legged tick Ixodes scapularis after infection with Borrelia burgdorferi, the causative agent of Lyme disease in the United States. Small RNA sequencing and qRT-PCR analysis were used to identify and validate differentially expressed I. scapularis salivary miRNAs. Small RNA-seq yielded 133,465,828 (≥18 nucleotides) and 163,852,135 (≥18 nucleotides) small RNA reads from Borrelia-infected and uninfected salivary glands for downstream analysis using the miRDeep2 algorithm. As such, 254 miRNAs were identified across all datasets, 25 of which were high confidence and 51 low confidence known miRNAs. Further, 23 miRNAs were differentially expressed in uninfected and infected salivary glands: 11 were upregulated and 12 were downregulated upon pathogen infection. Gene ontology and network analysis of target genes of differentially expressed miRNAs predicted roles in metabolic, cellular, development, cellular component biogenesis, and biological regulation processes. Several Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, including sphingolipid metabolism; valine, leucine and isoleucine degradation; lipid transport and metabolism; exosome biogenesis and secretion; and phosphate-containing compound metabolic processes, were predicted as targets of differentially expressed miRNAs. A qRT-PCR assay was utilized to validate the differential expression of miRNAs. This study provides new insights into the miRNAs expressed in I. scapularis salivary glands and paves the way for their functional manipulation to prevent or treat B. burgdorferi infection.

Comprehensive Analysis Reveals the Genetic and Pathogenic Diversity of Ralstonia solanacearum Species Complex and Benefits Its Taxonomic Classification

Ralstonia solanacearum species complex (RSSC) is a diverse group of plant pathogens that attack a wide range of hosts and cause devastating losses worldwide. In this study, we conducted a comprehensive analysis of 131 RSSC strains to detect their genetic diversity, pathogenicity, and evolution dynamics. Average nucleotide identity analysis was performed to explore the genomic relatedness among these strains, and finally obtained an open pangenome with 32,961 gene families. To better understand the diverse evolution and pathogenicity, we also conducted a series of analyses of virulence factors (VFs) and horizontal gene transfer (HGT) in the pangenome and at the single genome level. The distribution of VFs and mobile genetic elements (MGEs) showed significant differences among different groups and strains, which were consistent with the new nomenclatures of the RSSC with three distinct species. Further functional analysis showed that most HGT events conferred from Burkholderiales and played a great role in shaping the genomic plasticity and genetic diversity of RSSC genomes. Our work provides insights into the genetic polymorphism, evolution dynamics, and pathogenetic variety of RSSC and provides strong supports for the new taxonomic classification, as well as abundant resources for studying host specificity and pathogen emergence.

Genome-Wide Analyses of MADS-Box Genes in Humulus lupulus L. Reveal Potential Participation in Plant Development, Floral Architecture, and Lupulin Gland Metabolism

MADS-box transcription factors (TFs) are involved in multiple plant development processes and are most known during the reproductive transition and floral organ development. Very few genes have been characterized in the genome of Humulus lupulus L. (Cannabaceae), an important crop for the pharmaceutical and beverage industries. The MADS-box family has not been studied in this species yet. We identified 65 MADS-box genes in the hop genome, of which 29 encode type-II TFs (27 of subgroup MIKC^C and 2 MIKC*) and 36 type-I proteins (26 α, 9 β, and 1 γ). Type-II MADS-box genes evolved more complex architectures than type-I genes. Interestingly, we did not find FLOWERING LOCUS C (FLC) homologs, a transcription factor that acts as a floral repressor and is negatively regulated by cold. This result provides a molecular explanation for a previous work showing that vernalization is not a requirement for hop flowering, which has implications for its cultivation in the tropics. Analysis of gene ontology and expression profiling revealed genes potentially involved in the development of male and female floral structures based on the differential expression of ABC homeotic genes in each whorl of the flower. We identified a gene exclusively expressed in lupulin glands, suggesting a role in specialized metabolism in these structures. In toto, this work contributes to understanding the evolutionary history of MADS-box genes in hop, and provides perspectives on functional genetic studies, biotechnology, and crop breeding.

Transcriptome changes in the developing sugarcane culm associated with high yield and early-season high sugar content

Sugarcane, with its exceptional carbon dioxide assimilation, biomass and sugar yield, has a high potential for the production of bio-energy, bio-plastics and high-value products in the food and pharmaceutical industries. A crucial challenge for long-term economic viability and environmental sustainability is also to optimize the production of biomass composition and carbon sequestration. Sugarcane varieties such as KQ228 and Q253 are highly utilized in the industry. These varieties are characterized by a high early-season sugar content associated with high yield. In order to investigate these correlations, 1,440 internodes were collected and combined to generate a set of 120 samples in triplicate across 24 sugarcane cultivars at five different development stages. Weighted gene co-expression network analysis (WGCNA) was used and revealed for the first time two sets of co-expressed genes with a distinct and opposite correlation between fibre and sugar content. Gene identification and metabolism pathways analysis was used to define these two sets of genes. Correlation analysis identified a large number of interconnected metabolic pathways linked to sugar content and fibre content. Unsupervised hierarchical clustering of gene expression revealed a stronger level of segregation associated with the genotypes than the stage of development, suggesting a dominant genetic influence on biomass composition and facilitating breeding selection. Characterization of these two groups of co-expressed key genes can help to improve breeding program for high fibre, high sugar species or plant synthetic biology.

Computational identification and characterization of vascular wilt pathogen (Fusarium oxysporum f. sp. lycopersici) CAZymes in tomato xylem sap

Fusarium oxysporum f. sp. lycopersici is a devastating plant pathogenic fungi known for wilt disease in the tomato plant and secrete cell wall degrading enzymes. These enzymes are collectively known as carbohydrate-active enzymes (CAZymes), crucial for growth, colonization and pathogenesis. Therefore, the present study was aimed to identify and annotate pathogen CAZymes in the xylem sap of a susceptible tomato variety using downstream proteomics and meta servers. Further, structural elucidation and conformational stability analysis of the selected CAZyme families were done through homology modeling and molecular dynamics simulation. Among all the fungal proteins identified, the carbohydrate metabolic process was found to be enriched. Most of the annotated CAZymes belonged to the hydrolase and oxidoreductase families, and 90% were soluble and extracellular. Moreover, using a publically available interactome database, interactions were observed between the families acting on chitin, hemicellulose and pectin. Subsequently, important catalytic residues were identified in the candidate CAZymes belonging to carbohydrate esterase (CE8) and glycosyl hydrolase (GH18 and GH28). Further, essential dynamics after molecular simulation of 100 ns revealed the overall behavior of these CAZymes with distinct global minima and transition states in CE8. Thus, our study identified some of the CAZyme families that assist in pathogenesis and growth through host cell wall deconstruction with further structural insight into the selected CAZyme families.Communicated by Ramaswamy H. Sarma.

Verticillium dahliae Secretes Small RNA to Target Host MIR157d and Retard Plant Floral Transition During Infection

Bidirectional trans-kingdom RNA silencing [or RNA interference (RNAi)] plays a key role in plant-pathogen interactions. It has been shown that plant hosts export specific endogenous miRNAs into pathogens to inhibit their virulence, whereas pathogens deliver small RNAs (sRNAs) into plant cells to disturb host immunity. Here, we report a trans-kingdom fungal sRNA retarding host plant floral transition by targeting a miRNA precursor. From Arabidopsis plants infected with Verticillium dahliae, a soil-borne hemibiotrophic pathogenic fungus that causes wilt diseases in a wide range of plant hosts, we obtained a number of possible trans-kingdom V. dahliae sRNAs (VdsRNAs) by sequencing AGO1-immunoprecipitated sRNAs. Among these, a 24-nt VdsRNA derived from V. dahliae rRNA, VdrsR-1, was shown to be an actual trans-kingdom VdsRNA that targets the miR157d precursor MIR157d, resulting in increased rather than reduced miR157d accumulation in V. dahliae-infected plants. Consistent with the miR157 family in the regulation of vegetative and floral transitions by targeting SPL genes in several plant species, we detected two SPL genes, SPL13A/B, that were notably reduced in V. dahliae-infected and VdrsR-1-expressing plants compared with control plants. Furthermore, V. dahliae-infected and VdrsR-1-expressing plants also displayed delayed vegetative phase change and floral transition compared to control plants. Taken together, we disclosed a novel mode of action for a trans-kingdom fungal sRNA, VdrsR-1, which was secreted into host cells to modulate plant floral transition by employing the miR157d/SPL13A/B regulatory module, leading to prolonged host vegetative growth that would undoubtedly benefit fungal propagation.

Characteristics of the MicroRNA Expression Profile of Exosomes Released by Vero Cells Infected with Porcine Epidemic Diarrhea Virus

Exosomes are nanoscale vesicles actively secreted by a variety of cells. They contain regulated microRNA (miRNA), allowing them to function in intercellular communication. In the present study, the role of exosomal miRNAs in porcine epidemic diarrhea virus (PEDV) infection was investigated using exosomes isolated from Vero cells infected with PEDV. The results of transmission electron microscopy observation showed that the exosomes are spherical in shape, uniform in size, and negatively stained in the membrane. Nanoparticle tracking analysis showed that the average exosome particle size is 130.5 nm. The results of miRNA sequencing showed that, compared with the control group, a total of 115 miRNAs are abnormally expressed in the exosomes of infected cells. Of these, 80 miRNAs are significantly upregulated and 35 miRNAs are significantly downregulated. Functional annotation analysis showed that the differentially expressed miRNAs are associated with PEDV infection through interaction with the cAMP, Hippo, TGF-beta, HIF-1, FoxO, MAPK, and Ras signaling pathways. Thus, our findings provide important information about the effects of PEDV infection on exosomal miRNA expression and will aid the search for potential anti-PEDV drug candidates.

Natural variation and improved genome annotation of the emerging biofuel crop field pennycress (Thlaspi arvense)

The Brassicaceae family comprises more than 3,700 species with a diversity of phenotypic characteristics, including seed oil content and composition. Recently, the global interest in Thlaspi arvense L. (pennycress) has grown as the seed oil composition makes it a suitable source for biodiesel and aviation fuel production. However, many wild traits of this species need to be domesticated to make pennycress ideal for cultivation. Molecular breeding and engineering efforts require the availability of an accurate genome sequence of the species. Here, we describe pennycress genome annotation improvements, using a combination of long- and short-read transcriptome data obtained from RNA derived from embryos of 22 accessions, in addition to public genome and gene expression information. Our analysis identified 27,213 protein-coding genes, as well as on average 6,188 biallelic SNPs. In addition, we used the identified SNPs to evaluate the population structure of our accessions. The data from this analysis support that the accession Ames 32872, originally from Armenia, is highly divergent from the other accessions, while the accessions originating from Canada and the United States cluster together. When we evaluated the likely signatures of natural selection from alternative SNPs, we found 7 candidate genes under likely recent positive selection. These genes are enriched with functions related to amino acid metabolism and lipid biosynthesis and highlight possible future targets for crop improvement efforts in pennycress.

Novel miRNA-SSRs for Improving Seed Hardness Trait of Pomegranate (Punica granatum L.)

Present research discovered novel miRNA-SSRs for seed type trait from 761 potential precursor miRNA sequences of pomegranate. SSR mining and BLASTx of the unique sequences identified 69 non-coding pre-miRNA sequences, which were then searched for BLASTn homology against Dabenzi genome. Sixty three true pri-miRNA contigs encoding 213 pre-miRNAs were predicted. Analysis of the resulting sequences enabled discovery of SSRs within pri-miRNA (227) and pre-miRNA sequences (79). A total of 132 miRNA-SSRs were developed for seed type trait from 63 true pri-miRNAs, of which 46 were specific to pre-miRNAs. Through ePCR, 123 primers were validated and mapped on eight Tunisia chromosomes. Further, 80 SSRs producing specific amplicons were ePCR-confirmed on multiple genomes i.e. Dabenzi, Taishanhong, AG2017 and Tunisia, yielding a set of 63 polymorphic SSRs (polymorphism information content ≥0.5). Of these, 32 miRNA-SSRs revealed higher polymorphism level (89.29%) when assayed on six pomegranate genotypes. Furthermore, target prediction and network analysis suggested a possible association of miRNA-SSRs i.e. miRNA_SH_SSR69, miRNA_SH_SSR36, miRNA_SH_SSR103, miRNA_SH_SSR35 and miRNA_SH_SSR53 with seed type trait. These miRNA-SSRs would serve as important genomic resource for rapid and targeted improvement of seed type trait of pomegranate.

Hierarchical multi-label classification based on LSTM network and Bayesian decision theory for LncRNA function prediction

Growing evidence shows that long noncoding RNAs (lncRNAs) play an important role in cellular biological processes at multiple levels, such as gene imprinting, immune response, and genetic regulation, and are closely related to diseases because of their complex and precise control. However, most functions of lncRNAs remain undiscovered. Current computational methods for exploring lncRNA functions can avoid high-throughput experiments, but they usually focus on the construction of similarity networks and ignore the certain directed acyclic graph (DAG) formed by gene ontology annotations. In this paper, we view the function annotation work as a hierarchical multilabel classification problem and design a method HLSTMBD for classification with DAG-structured labels. With the help of a mathematical model based on Bayesian decision theory, the HLSTMBD algorithm is implemented with the long-short term memory network and a hierarchical constraint method DAGLabel. Compared with other state-of-the-art algorithms, the results on GOA-lncRNA datasets show that the proposed method can efficiently and accurately complete the label prediction work.

Functional analysis of zona pellucida domain protein Dusky in Tribolium castaneum

The zona pellucida domain protein Dusky (Dy) plays a vital role in wing morphogenesis in insects, but little information on its function has been reported. In this study, we found that dy regulated wing cell size, larval and pupal duration, and the metabolism of amino acid and 20-hydroxyecdysone in Tribolium castaneum. Using RNA-seq, 413 differentially expressed genes were identified between physiological buffer-injected and dy-double-stranded RNA-treated larvae, including 88 downregulated genes and 325 upregulated genes. Among these genes, dy knockdown increased CYP18A1 expression to elevate the 26-hydroxylation of 20-hydroxyecdysone, which ultimately led to growth defects in wing cells. Silencing of dy upregulated the transcription of genes encoding tyrosine aminotransferase, 4-hydroxyphenylpyruvate dioxygenase, homogentisate 1, 2-dioxygenase, and Pale to promote the catabolism of tyrosine and phenylalanine, which eventually reduced amino acid content. Furthermore, dy knockdown upregulated 4E-BP expression, and 4E-BP silencing partially phenocopied dy RNA interference-mediated wing morphogenesis. These results suggest that Dy controls 20-hydroxyecdysone and amino acid metabolism to regulate wing morphogenesis in the insect.

Potential receptors in Fenneropenaeus merguiensis ovary and role of saxophone, the bone morphogenetic protein receptor, in ovarian development

Receptors, which play an initial role in signaling pathways in several physiological processes, including reproduction, are among the several molecular factors that control ovarian development in organisms. This study aimed to identify and study receptors potentially involved in controlling the reproductive process of female banana shrimp, Fenneropenaeus merguiensis. Ovarian transcriptomes derived from 4 developmental stages were generated by RNA sequencing. A total of 53,763 transcripts were obtained from the de novo assembled transcriptome, and 663 genes were identified as receptors. Among them, 185 receptors were differentially expressed during ovarian development. Fifteen of these differentially expressed receptors showed distinct expression patterns that were validated by RT-qPCR. Bone morphogenetic protein receptors (BMPR) and their signaling genes were investigated for their roles in shrimp vitellogenesis. The expressions of F. merguiensis saxophone (FmSax), a BMP type I receptor, and BMP type II receptor (FmBMPRII) as well as FmMad, FmMed, and FmSMAD3 were significantly altered during ovarian development. RNA interference was used to investigate the role of FmSax in vitellogenesis. The result indicated that the expression of vitellogenin (Vg) was significantly reduced in both ovary and hepatopancreas of FmSax-knockdown shrimp compared to control shrimp. Furthermore, in FmSax-silencing shrimp, FmBMPRII, FmMad, and FmMed expressions were decreased as well as Vg expression. These findings suggest that FmSax positively regulates Vg synthesis via the BMP signaling pathway.

Transcriptomic analysis of methyl jasmonate treatment reveals gene networks involved in drought tolerance in pearl millet

Water deficit stress at the early stage of development is one of the main factors limiting pearl millet production. One practice to counteract this limitation would be to resort to the application of hormones to stimulate plant growth and development at critical stages. Exogenous methyl jasmonate (MeJA) can improve drought tolerance by modulating signaling, metabolism, and photosynthesis pathways, therefore, we assumed that can occur in pearl millet during the early stage of development. To decipher the molecular mechanisms controlling these pathways, RNAseq was conducted in two pearl millet genotypes, drought-sensitive SosatC88 and drought-tolerant Souna3, in response to 200 μM of MeJA. Pairwise comparison between the MeJA-treated and non-treated plants revealed 3270 differentially expressed genes (DEGs) among 20,783 transcripts in SosatC88 and 127 DEGs out of 20,496 transcripts in Souna3. Gene ontology (GO) classification assigned most regulated DEGs in SosatC88 to heme binding, oxidation-reduction process, response to oxidative stress and membrane, and in Souna3 to terpene synthase activity, lyase activity, magnesium ion binding, and thylakoid. The Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis reveals that DEGs in SosatC88 are related to the oxidation-reduction process, the biosynthesis of other secondary metabolites, the signal transduction, and the metabolism of terpenoids, while in Souna3, DEGs are related to the metabolism of terpenoids and the energy metabolism. Two genes encoding a diterpenoid biosynthesis-related (Pgl_GLEAN_10009413) and a Glutathione S transferase T3 (Pgl_GLEAN_10034098) were contra-regulated between SosatC88 and Souna3. Additionally, five random genes differentially expressed by RNAseq were validated using qPCR, therefore, they are potential targets for the development of novel strategies breeding schemes for plant growth under water deficit stress. These insights into the molecular mechanisms of pearl millet genotype tolerance at the early stage of development contribute to the understanding of the role of hormones in adaptation to drought-prone environments.

An extremely streamlined macronuclear genome in the free-living protozoan Fabrea salina

Ciliated protists are among the oldest unicellular organisms with a heterotrophic lifestyle and share a common ancestor with Plantae. Unlike any other eukaryotes, there are two distinct nuclei in ciliates with separate germline and somatic cell functions. Here, we assembled a near-complete macronuclear genome of Fabrea salina, which belongs to one of the oldest clades of ciliates. Its extremely minimized genome (18.35 Mb) is the smallest among all free-living heterotrophic eukaryotes and exhibits typical streamlined genomic features, including high gene density, tiny introns, and shrinkage of gene paralogs. Gene families involved in hypersaline stress resistance, DNA replication proteins, and mitochondrial biogenesis are expanded, and the accumulation of phosphatidic acid may play an important role in resistance to high osmotic pressure. We further investigated the morphological and transcriptomic changes in the macronucleus during sexual reproduction and highlighted the potential contribution of macronuclear residuals to this process. We believe that the minimized genome generated in this study provides novel insights into the genome streamlining theory and will be an ideal model to study the evolution of eukaryotic heterotrophs.

Alternative splicing diversifies the transcriptome and proteome of the rice blast fungus during host infection

Alternative splicing (AS) contributes to diversifying and regulating cellular responses to environmental conditions and developmental cues by differentially producing multiple mRNA and protein isoforms from a single gene. Previous studies on AS in pathogenic fungi focused on profiling AS isoforms under a limited number of conditions. We analysed AS profiles in the rice blast fungus Magnaporthe oryzae, a global threat to rice production, using high-quality transcriptome data representing its vegetative growth (mycelia) and multiple host infection stages. We identified 4,270 AS isoforms derived from 2,413 genes, including 499 genes presumably regulated by infection-specific AS. AS appears to increase during infection, with 32.7% of the AS isoforms being produced during infection but absent in mycelia. Analysis of the isoforms observed at each infection stage showed that 636 AS isoforms were more abundant than corresponding annotated mRNAs, especially after initial hyphal penetration into host cell. Many such dominant isoforms were predicted to encode regulatory proteins such as transcription factors and phospho-transferases. We also identified the genes encoding distinct proteins via AS and confirmed the translation of some isoforms via a proteomic analysis, suggesting potential AS-mediated neo-functionalization of some genes during infection. Comprehensive profiling of the pattern of genome-wide AS during multiple stages of rice-M. oryzae interaction established a foundational resource that will help investigate the role and regulation of AS during rice infection.

An Insight Into the microRNA Profile of the Ectoparasitic Mite Varroa destructor (Acari: Varroidae), the Primary Vector of Honey Bee Deformed Wing Virus

The remarkably adaptive mite Varroa destructor is the most important honey bee ectoparasite. Varroa mites are competent vectors of deformed wing virus (DWV), and the Varroa-virus complex is a major determinant of annual honey bee colony mortality and collapse. MicroRNAs (miRNAs) are 22-24 nucleotide non-coding RNAs produced by all plants and animals and some viruses that influence biological processes through post-transcriptional regulation of gene expression. Knowledge of miRNAs and their function in mite biology remains limited. Here we constructed small RNA libraries from male and female V. destructor using Illumina's small RNA-Seq platform. A total of 101,913,208 and 91,904,732 small RNA reads (>18 nucleotides) from male and female mites were analyzed using the miRDeep2 algorithm. A conservative approach predicted 306 miRNAs, 18 of which were upregulated and 13 downregulated in female V. destructor compared with males. Quantitative real-time PCR validated the expression of selected differentially-expressed female Varroa miRNAs. This dataset provides a list of potential miRNA targets involved in regulating vital Varroa biological processes and paves the way for developing strategies to target Varroa and their viruses.

Linking genome wide RNA sequencing with physio-biochemical and cytological responses to catalogue key genes and metabolic pathways for alkalinity stress tolerance in lentil (Lens culinaris Medikus)

BACKGROUND

Alkaline soils cause low productivity in crop plants including lentil. Alkalinity adaptation strategies in lentil were revealed when morpho-anatomical and physio-biochemical observations were correlated with transcriptomics analysis in tolerant (PDL-1) and sensitive (L-4076) cultivars at seedling stage.

RESULTS

PDL-1 had lesser salt injury and performed better as compared to L-4076. Latter showed severe wilting symptoms and higher accumulation of Na⁺ and lower K⁺ in roots and shoots. PDL-1 performed better under high alkalinity stress which can be attributed to its higher mitotic index, more accumulation of K⁺ in roots and shoots and less aberrantly dividing cells. Also, antioxidant enzyme activities, osmolytes' accumulation, relative water content, membrane stability index and abscisic acid were higher in this cultivar. Differentially expressed genes (DEGs) related to these parameters were upregulated in tolerant genotypes compared to the sensitive one. Significantly up-regulated DEGs were found to be involved in abscisic acid (ABA) signalling and secondary metabolites synthesis. ABA responsive genes viz. dehydrin 1, 9-cis-epoxycarotenoid dioxygenase, ABA-responsive protein 18 and BEL1-like homeodomain protein 1 had log₂fold change above 4.0. A total of 12,836 simple sequence repeats and 4,438 single nucleotide polymorphisms were identified which can be utilized in molecular studies.

CONCLUSIONS

Phyto-hormones biosynthesis-predominantly through ABA signalling, and secondary metabolism are the most potent pathways for alkalinity stress tolerance in lentil. Cultivar PDL-1 exhibited high tolerance towards alkalinity stress and can be used in breeding programmes for improving lentil production under alkalinity stress conditions.

Comparative transcriptomes of three different skin sites for the Asiatic toad (Bufo gargarizans)

Toads release toxic dry secretions from glands in their skin. Toxin possesses a wide range of biological effects, but little is known about its specific gene expression pattern and regulatory mechanisms. The Asiatic toad (Bufo gargarizans) is widely used to produce toxin. Here, we explored the gene expression of 30 tissue samples from three different skin sites (parotoid gland, dorsal skin, and abdomen skin) of B. gargarizans. After de novo assembly, 783,130 unigenes with an average length of 489 bp (N50 = 556 bp) were obtained. A total of 9,248 significant differentially expressed genes (DEGs) were detected. There were 8,819 DEGs between the parotoid gland and abdomen skin and 1,299 DEGs between the dorsal skin and abdomen skin, while only 1,283 DEGs were obtained between the parotoid gland and dorsal skin. Through enrichment analysis, it was found that the detected differential gene expressions corresponded to the different functions of different skin sites. Our key findings were the genetic expression of toxin secretion, the protection function of skin, and the related genes such as HSD3B, Cyp2c, and CAT, LGALS9. In conclusion, we provide useful transcript resources to study the gene expression and gene function of B. gargarizans and other amphibians. The detected DEGs between different sites of the skin provided better insights into the genetic mechanisms of toxin secretion and the protection function of skin for amphibians.

iTRAQ Proteomic Analysis of Interactions Between 20E and Phospholipase C in Apolygus lucorum (Meyer-Dür)

Polyphagous Apolygus lucorum has become the dominant insect in Bacillus thuringiensis (Bt) cotton fields. Hormone 20-hydroxyecdysone (20E) regulates multiple insect development and physiology events. 20E responses are controlled by pathways triggered by phospholipase C (PLC)-associated proteins. However, 20E-modulated genes and related proteins that can be affected by PLC still remain unknown. Here, isobaric tag for relative and absolute quantitation (iTRAQ) and immunoblotting techniques were used to compare differentially expressed proteins (DEPs) in A. lucorum in response to the treatment of 20E and the PLC inhibitor U73122 as well as their combination. A total of 1,624 non-redundant proteins and 97, 248, 266 DEPs were identified in the 20E/control, U73122/control, and 20E + U73122/control groups, respectively. Only 8 DEPs, including pathogenesis-related protein 5-like, cuticle protein 19.8, trans-sialidase, larval cuticle protein A2B-like, cathepsin L1, hemolymph juvenile hormone-binding protein, ATP-dependent RNA helicase p62-like, and myosin-9 isoform X1, were detected in all three groups. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed that the DEPs were involved in diverse signaling pathways. The results were validated by immunoblotting, which highlighted the reliability of proteomics analysis. These findings provided novel insights into the function of PLC in 20E signaling pathway in A. lucorum.

PlantGF: an analysis and annotation platform for plant gene families

Gene families contain genes that come from the same ancestor and have similar sequences and structures. They perform certain specific functions within and among different species. Currently, there is no complete process or platform for the rapid analysis of plant gene families. In this study, a comprehensive query and analysis platform of plant gene families, the Plant Gene Family Platform (PlantGF), was constructed. The platform is composed of four main parts: Search, Tools, Statistics and Auxiliary. A total of 2 909 580 gene family members were identified from 138 plant species in PlantGF. The data can be queried in the Search section through a user-friendly interface. A general process for gene family analysis, having nine steps, is provided. The platform also includes four online tools (HMM-Search, BLAST, MAFFT and HMMER) in the Tools section for useful additional analyses. The statistical analysis of the relevant gene families is shown on the Statistics page. Auxiliary pages are provided for data downloading. The datasets for all 138 plant species' protein sequences and their gene families can be acquired on the Download page. A user's manual and some useful links are displayed on the Manual and Links pages, respectively. To the best of our knowledge, PlantGF is the first comprehensive platform for studying plant gene families, and it will make important contributions to plant gene family-related research. Database URL: http://biodb.sdau.edu.cn/PGF/index.html.

Insights from the analysis of draft genome sequence of Crocus sativus L

Saffron (Crocus sativus L.) is the low yielding plant of medicinal and economic importance. Therefore, it is of interest to report the draft genome sequence of C. sativus. The draft genome of C. sativus has been assembled using Illumina sequencing and is 3.01 Gb long covering 84.24% of genome. C. sativus genome annotation identified 53,546 functional genes (including 5726 transcription factors), 862,275 repeats and 964,231 SSR markers. The genes involved in the apocarotenoids biosynthesis pathway (crocin, crocetin, picrocrocin, and safranal) were found in the draft genome analysis.

MOET: a web-based gene set enrichment tool at the Rat Genome Database for multiontology and multispecies analyses

Biological interpretation of a large amount of gene or protein data is complex. Ontology analysis tools are imperative in finding functional similarities through overrepresentation or enrichment of terms associated with the input gene or protein lists. However, most tools are limited by their ability to do ontology-specific and species-limited analyses. Furthermore, some enrichment tools are not updated frequently with recent information from databases, thus giving users inaccurate, outdated or uninformative data. Here, we present MOET or the Multi-Ontology Enrichment Tool (v.1 released in April 2019 and v.2 released in May 2021), an ontology analysis tool leveraging data that the Rat Genome Database (RGD) integrated from in-house expert curation and external databases including the National Center for Biotechnology Information (NCBI), Mouse Genome Informatics (MGI), The Kyoto Encyclopedia of Genes and Genomes (KEGG), The Gene Ontology Resource, UniProt-GOA, and others. Given a gene or protein list, MOET analysis identifies significantly overrepresented ontology terms using a hypergeometric test and provides nominal and Bonferroni corrected P-values and odds ratios for the overrepresented terms. The results are shown as a downloadable list of terms with and without Bonferroni correction, and a graph of the P-values and number of annotated genes for each term in the list. MOET can be accessed freely from https://rgd.mcw.edu/rgdweb/enrichment/start.html.