WEGO: a web tool for plotting GO annotations

Transcriptome analysis reveals key regulatory networks and genes involved in the acquisition of cold stress memory in pepper seedlings

Temperature is an important limiting factor in the counter-seasonal cultivation of pepper. Currently, there are no studies on transcriptomic analysis of 'cold stress memory' in pepper. In this study, in order to understand the mechanism of 'cold stress memory' in pepper (Capsicum annuum L.), seedlings were subjected to the following treatments: normal temperature treatment (P0), the first cold treatment for 3 days (P3), the recovery temperature treatment for 3 days (R3), and another cold treatment for 3 days (RP3). The results showed that P3 plants wilted the most, RP3 the second and R3 the least. Leaf reactive oxygen species (ROS) and electrolyte leakage were the most in P3, the second in RP3 and the least in R3. In addition, RP3 had the highest accumulation of zeaxanthin, violaxanthin and β-cryptoxanthin, followed by P3, and R3 had the least. These results suggest that pepper seedlings are characterized by 'cold stress memory'. Transcriptomics was used to analyze the key genes and transcription factors involved in the biosynthesis of zeaxanthin, violaxanthin and β-cryptoxanthin during the formation of 'cold stress memory'. This study provides candidate genes and transcription factors for an in-depth study of the cold tolerance mechanism in pepper.

3'UTR of tobacco vein mottling virus regulates downstream GFP expression and changes in host gene expression

Introduction

Tobacco vein mottling virus (TVMV) is a member of the family Potyviridae. The 3' untranslated region (3'UTR) of viral genomic RNA has been reported to significantly impact viral infection. Nevertheless, the role of the TVMV 3'UTR during viral infection remains unknown.

Methods

Here, a 3'UTR-GFP expression vector was transiently expressed in Nicotiana benthamiana, in which the 3'UTR of TVMV was introduced upstream of the green fluorescent protein (GFP) gene. Transcriptome sequencing was performed to analyze the genes associated with plant resistance. The effect of the TVMV 3'UTR on GFP expression was studied using an Agrobacterium-mediated transient expression assay, revealing that the TVMV 3'UTR significantly inhibited GFP expression. Transcriptome analysis of differentially expressed genes in 3'UTR-GFP in N. benthamiana was performed to elucidate the why the TVMV 3'UTR inhibited GFP expression.

Results

Eighty genes related to plant disease resistance were differentially expressed, including 29 upregulated and 51 downregulated genes. Significantly upregulated genes included those encoding the calcium-binding protein CML24, leucine-rich repeat receptor-like tyrosine-protein kinase, and respiratory burst oxidase homolog protein E. The significantly downregulated genes included calcium-binding protein 7, ethylene-responsive transcription factor 10, endoglucanase 5, and receptor-like protein kinase.

Discussion

These findings indicate that the 3'UTR of TVMV may inhibit the expression of GFP gene by inducing the expression of plant resistance genes. This study provides a theoretical basis for further research on the function and mechanism of the TVMV 3'UTR.

Unraveling potential gene biomarkers for dengue infection through RNA sequencing

Dengue virus hijacks host cell mechanisms and immune responses in order to replicate efficiently. The interaction between the host and the virus affects the host's gene expression, which remains largely unexplored. This pilot study aimed to profile the host transcriptome as a potential strategy for identifying specific biomarkers for dengue prediction and detection. High-throughput RNA sequencing (RNA-seq) was employed to generate host transcriptome profiles in 16 dengue patients and 10 healthy controls. Differentially expressed genes (DEGs) were identified in patients with severe dengue and those with dengue with warning signs compared to healthy individuals. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed to elucidate the functions of upregulated and downregulated genes. Compared to healthy controls, 6466 genes were significantly differentially expressed (p < 0.05) in the dengue with warning signs group and 3082 genes in the severe dengue group, with over half being upregulated. The major KEGG pathways implicated included transport and catabolism (14.4%-16.3%), signal transduction (6.6%-7.3%), global and overview maps (6.7%-7.1%), viral diseases (4.6%-4.8%), and the immune system (4.4%-4.6%). Several genes exhibited consistent and significant upregulation across all dengue patients, regardless of severity: Interferon alpha inducible protein 27 (IFI27), Potassium Channel Tetramerization Domain Containing 14 (KCTD14), Syndecan 1 (SDC1), DCC netrin 1 receptor (DCC), Ubiquitin C-terminal hydrolase L1 (UCHL1), Marginal zone B and B1 cell-specific protein (MZB1), Nestin (NES), C-C motif chemokine ligand 2 (CCL2), TNF receptor superfamily member 17 (TNFSF17), and TNF receptor superfamily member 13B (TNFRSF13B). Further analysis revealed potential biomarkers for severe dengue prediction, including TNF superfamily member 15 (TNFSF15), Plasminogen Activator Inhibitor-2 (SERPINB2), motif chemokine ligand 7 (CCL7), aconitate decarboxylase 1 (ACOD1), Metallothionein 1G (MT1G), and Myosin Light Chain Kinase (MYLK2), which were expressed 3.5 times, 2.9 times, 2.3 times, 2.1 times, 1.7 times, and 1.4 times greater, respectively, than dengue patients exhibiting warning signs. The identification of these host biomarkers through RNA-sequencing holds promising implications and potential to augment existing dengue detection algorithms, contributing significantly to improved diagnostic and prognostic capabilities.

Proteo-transcriptomic profiles reveal genetic mechanisms underlying primary hair follicle development in coarse sheep fetal skin

Long hair trait represents a valuable genetic asset in Qinghai Tibetan sheep, with its quality and yield being contingent upon the characteristics of hair follicles (HFs). This study aims to elucidate the genetic mechanism underlying primary hair follicles (PFs) formation through an integrated analysis of proteomics and transcriptomics. Samples were collected at key stages of fetal HF formation (E65 and E85) for histological observation, revealing significant alterations in the microstructure of PF (E65) during the developmental process. In this study, a comprehensive analysis revealed a total of 217 overlapping genes that exhibited concordant expression patterns at both the proteomic and transcriptomic levels. Furthermore, to ensure the reliability of our findings, we employed parallel response monitoring (PRM) to validate the obtained proteomic data. The protein-protein interaction (PPI) network diagram highlights five hub core proteins (TTN, IGTA2, F2, EGFR, and MYH14). These differentially expressed proteins (DEPs) play crucial roles in metabolic processes, cell adhesion, and diverse biological processes. The potential synergy between transcriptional regulation and post-translational modifications plays a pivotal role in governing the initiation PF development. The findings presented in this study offer innovative insights into the molecular mechanisms underlying HFs generation and establish a robust foundation for targeted breeding strategies aimed at augmenting wool traits in sheep. SIGNIFICANCE: The composition of coarse hair primarily consists of long, myelinated fibers originating from primary hair follicles. Sheep fetal skin initiates the formation of primary hair follicles around E65, followed by the development of secondary hair follicles around E85. Conducting differential proteomic and transcriptomic analyses during these developmental stages enhances our understanding of the molecular mechanisms underlying primary hair follicle development and offers valuable insights for sustainable utilization of high-quality germplasm resources.

Protein profile of extracellular vesicles derived from adult Parascaris spp

BACKGROUND

Parascaris spp. represent a significant threat to equine health worldwide, particularly in foals. The long-term survival of parasites in the host necessitates persistent modulation of the host immune response. Intercellular communication achieved through the exchange of molecules via extracellular vesicles (EVs) released from the parasite could be a crucial factor in this regard. This study aimed to isolate and characterize EVs released by adult male and female Parascaris worms and conduct a proteomic analysis to identify sex-specific proteins and potential immunomodulatory factors.

METHODS

Live adult Parascaris worms were collected, and EVs were isolated from spent culture media using differential ultracentrifugation. Nanoparticle tracking analysis and transmission electron microscopy confirmed the size, concentration, and morphology of the isolated EVs. Proteins within the isolated EVs were analyzed using mass spectrometry-based proteomics (LC-MS/MS).

RESULTS

Proteomic analysis revealed a total of 113 proteins in Parascaris EVs, with several proteins showing homology to known helminth exosome proteins and exhibiting immunomodulatory functions. Sex-specific differences in EV protein composition were observed, with a distinct abundance of C-type lectins in female EVs, suggesting potential sex-specific roles or regulation. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses revealed metabolic pathways shared between male and female Parascaris EVs, as well as differences in signal transduction, and cell growth and death pathways, indicating sex-specific variations.

CONCLUSIONS

These findings imply that Parascaris EVs and their protein cargo are complex. This data potentially opens avenues for discovering innovative approaches to managing and understanding helminth infection.

Analysis of Differential Gene Expression under Acute Lead or Mercury Exposure in Larval Zebrafish Using RNA-Seq

This study was first conducted to investigate the effects of acute lead exposure on developing zebrafish embryos or larvae from 24 to 120 h post-fertilization (hpf). Our data showed that treatment with 50-200 μM lead significantly affected larval survivability and morphology compared to the respective control. Second, we chose 120 hpf larvae treated with 12.5 μM lead for RNA sequencing due to its exposure level being sufficient to produce toxic effects with minimum death and lead bioaccumulation in developing zebrafish. A total of 137.45 million raw reads were obtained, and more than 86% of clean data were mapped to the zebrafish reference genome. Differential expression profiles generated 116 up- and 34 down-regulated genes upon lead exposure. The most enriched GO terms for representative DEGs were ion transport and lipid metabolism. Third, a comparison with the dataset of mercury-regulated gene expression identified 94 genes (64 up-regulated and 30 down-regulated) for exposure specific to lead, as well as 422 genes (338 up-regulated and 84 down-regulated) for exposure specific to mercury. In addition, 56 genes were co-regulated by micromolar mercury and lead treatment, and the expression of thirteen genes, including mt2, ctssb.1, prdx1, txn, sqrdl, tmprss13a, socs3a, trpv6, abcb6a, gsr, hbz, fads2, and zgc:92590 were validated by qRT-PCR. These genes were mainly associated with metal ion binding, proteolysis, antioxidant activity, signal transduction, calcium ion or oxygen transport, the fatty acid biosynthetic process, and protein metabolism. Taken together, these findings help better understand the genome-wide responses of developing zebrafish to lead or mercury and provide potential biomarkers for acute exposure to toxic metals.

Physiological Parameters and Transcriptomic Levels Reveal the Response Mechanism of Maize to Deep Sowing and the Mechanism of Exogenous MeJA to Alleviate Deep Sowing Stress

Deep sowing, as a method to mitigate drought and preserve soil moisture and seedlings, can effectively mitigate the adverse effects of drought stress on seedling growth. The elongation of the hypocotyl plays an important role in the emergence of maize seeds from deep-sowing stress. This study was designed to explore the function of exogenous methyl jasmonate (MeJA) in the growth of the maize mesocotyl and to examine its regulatory network. The results showed that the addition of a 1.5 μ mol L-1 MeJA treatment significantly increased the mesocotyl length (MES), mesocotyl and coleoptile length (MESCOL), and seedling length (SDL) of maize seedlings. Transcriptome analysis showed that exogenous MeJA can alleviate maize deep-sowing stress, and the differentially expressed genes (DEGs) mainly include ornithine decarboxylase, terpene synthase 7, ethylene responsive transcription factor 11, and so on. In addition, candidate genes that may regulate the length of maize hypocotyls were screened by Weighted Gene Co-expression Network Analysis (WGCNA). These genes may be involved in the growth of maize hypocotyls through transcriptional regulation, histones, ubiquitin protease, protein binding, and chlorophyll biosynthesis and play an important role in maize deep-sowing tolerance. Our research findings may provide a theoretical basis for determining the tolerance of maize to deep-sowing stress and the mechanism of exogenous hormone regulation of deep-sowing stress.

Hybrid de novo whole genome assembly of lipopeptide producing novel Bacillus thuringiensis strain NBAIR BtAr exhibiting antagonistic activity against Sclerotium rolfsii

Bacillus thuringiensis Berliner is recognized as a predominant bioinsecticide but its antifungal potential has been relatively underexplored. A novel B. thuringiensis strain NBAIR BtAr was isolated and morphologically characterized using light and scanning electron microscopy, revealing presence of bipyramidal, cuboidal, and spherical parasporal crystals. The crude form of lipopeptides was extracted from NBAIR BtAr and assessed for its antagonistic activity in vitro, and demonstrated 100 % inhibition of Sclerotium rolfsii Sacc. at a minimum inhibitory concentration of 50 μL of the crude lipopeptide extract per mL of potato dextrose agar. To identify the antagonistic genes responsible, we performed whole genome sequencing of NBAIR BtAr, revealing the presence of circular chromosome of 5,379,913 bp and 175,362 bp plasmid with 36.06 % guanine-cytosine content and 5814 protein-coding sequences. Average nucleotide identity and whole genome phylogenetic analysis delineated the NBAIR BtAr strain as konkukian serovar. Gene ontology analysis revealed associations of 1474, 1323, and 1833 genes with biological processes, molecular function, and cellular components, respectively. Antibiotics & secondary metabolite analysis shell analysis of the whole genome yielded secondary metabolites biosynthetic gene clusters with 100 %, 85 %, 40 %, and 35 % similarity for petrobactin, bacillibactin, fengycin, and paenilamicin, respectively. Also, novel biosynthetic gene clusters, along with antimicrobial genes, including zwittermicin A, chitinase, and phenazines, were identified. Moreover, the presence of eight bacteriophage sequences, 18 genomic islands, insertion sequences, and one CRISPR region indicated prior occurrences of genetic exchange and thus improved competitive fitness of the strain. Overall, the whole genome sequence of NBAIR BtAr is presented, with its taxonomic classification and critical genetic attributes that contribute to its strong antagonistic activity against S. rolfsii.

Transcriptome analysis reveals mechanisms of metabolic detoxification and immune responses following farnesyl acetate treatment in Metisa plana

Metisa plana is a widespread insect pest infesting oil palm plantations in Malaysia. Farnesyl acetate (FA), a juvenile hormone analogue, has been reported to exert in vitro and in vivo insecticidal activity against other insect pests. However, the insecticidal mechanism of FA on M. plana remains unclear. Therefore, this study aims to elucidate responsive genes in M. plana in response to FA treatment. The RNA-sequencing reads of FA-treated M. plana were de novo-assembled with existing raw reads from non-treated third instar larvae, and 55,807 transcripts were functionally annotated to multiple protein databases. Several insecticide detoxification-related genes were differentially regulated among the 321 differentially expressed transcripts. Cytochrome P450 monooxygenase, carboxylesterase, and ATP-binding cassette protein were upregulated, while peptidoglycan recognition protein was downregulated. Innate immune response genes, such as glutathione S-transferases, acetylcholinesterase, and heat shock protein, were also identified in the transcriptome. The findings signify that changes occurred in the insect's receptor and signaling, metabolic detoxification of insecticides, and immune responses upon FA treatment on M. plana. This valuable information on FA toxicity may be used to formulate more effective biorational insecticides for better M. plana pest management strategies in oil palm plantations.

Multi-Omics Revealed Regulatory Mechanisms Underlying the Flowering of Ferula sinkiangensis across Three Dimensions

Backgroud/Objectives: Ferula spp. is an essential crop in Central Asia with pronounced economic benefits governed by its flowering process. However, the mechanisms of the flowering phenotype remain unclear. Methods: In this study, using F. sinkiangensis as a model plant, we integrated transcriptome, proteome, and metabolome analyses to compare the multilayer differences in leaves and roots of plants with flowering and unflowering phenotypes. Results: We found that several variations in the transcriptome, proteome, and metabolome were closely associated with flowering. The Photosynthesis and Phenylpropanoid biosynthesis pathways in plants with the flowering phenotype were more active. Additionally, three flowering genes, named FL2-FL4, were upregulated in the leaves of flowering plants. Notably, six transcription factors were potentially responsible for regulating the expression of FL2-FL4 in the leaves to mediate flowering process of F. sinkiangensis. Moreover, genes relevant to Photosynthesis and Phenylpropanoid biosynthesis were also involved in regulating the expression of FL2-FL4 in flowering plants. Conclusions: The active regulation network together with Photosynthesis and Phenylpropanoid biosynthesis were essential for inducing the expression of flowering-related genes in leaves to promote the flowering process of F. sinkiangensis.

Unraveling the molecular landscape of breast muscle development in domestic Yuzhong pigeons and European meat pigeon: Insights from Iso-seq and RNA-seq analysis

The mechanisms governing gene regulation in domestic Yuzhong pigeon breast muscle development remain largely elusive. Here, we conducted a comparative analysis using Iso-seq and RNA-seq data from domestic Yuzhong pigeons and European meat pigeons to uncover signaling pathways and genes involved in breast muscle development. The Iso-seq data from domestic Yuzhong pigeons yielded 131,377,075 subreads, resulting in 16,587 non-redundant high-quality full-length transcripts post-correction. Furthermore, utilizing pfam, CPC, PLEK, and CPAT, we predicted 5575, 4973, 2333, and 4336 lncRNAs, respectively. Notably, several genes potentially implicated in breast muscle development were identified, including tropomyosin beta chain, myosin regulatory light chain 2, and myosin binding protein C. KEGG enrichment analysis revealed critical signaling pathways in breast muscle development, spanning carbon metabolism, biosynthesis of amino acids, glycolysis/gluconeogenesis, estrogen signaling, PI3K-AKT signaling, protein processing in the endoplasmic reticulum, oxidative phosphorylation, pentose phosphate pathway, fructose and mannose metabolism, and tight junctions. These findings offer insights into the biological processes driving breast muscle development in domestic Yuzhong pigeon, contributing to our understanding of this complex phenomenon.

Transcriptome analysis reveals genes associated with late blight resistance in potato

Late blight is a serious disease of potato worldwide. Our study aimed to unveil genes involved in late blight resistance in potato by RNA-seq analysis after artificial inoculation under controlled conditions. In this study, two potato somatic hybrids (P7 and Crd6) and three varieties such as Kufri Girdhari, Kufri Jyoti and Kufri Bahar (control) were used. Transcriptiome analysis revealed statistically significant (p < 0.05) differentially expressed genes (DEGs), which were analysed into up-regulated and down-regulated genes. Further, DEGs were functionally characterized by the Gene Ontology annotations and the Kyoto Encyclopedia of Genes and Genomes pathways. Overall, some of the up-regulated genes in resistant genotypes were disease resistance proteins such as CC-NBS-LRR resistance protein, ankyrin repeat family protein, cytochrome P450, leucine-rich repeat family protein/protein kinase family, and MYB transcription factor. Sequence diversity analysis based on 38 peptide sequences representing 18 genes showed distinct variation and the presence of three motifs in 15 amino acid sequences. Selected genes were also validated by real-time quantitative polymerase chain reaction analysis. Interestingly, gene expression markers were developed for late blight resistant genotypes. Our study elucidates genes involved in imparting late blight resistance in potato, which will be beneficial for its management strategies in the future.

Hepatic transcriptome, transcriptional effects and antioxidant responses in Poecilia vivipara exposed to sanitary sewage

Aquatic environments are subject to threats from multiple human activities, particularly through the release of untreated sanitary sewage into the coastal environments. These effluents contain a large group of natural or synthetic compounds referred to as emerging contaminants. Monitoring the types and quantities of toxic substances in the environment, especially complex mixtures, is an exhausting and challenging task. Integrative effect-based tools, such as biomarkers, are recommended for environmental quality monitoring programs. In this study, fish Poecilia vivipara were exposed for 24 and 96 h to raw untreated sewage diluted 33 % (v/v) in order to identify hepatic genes to be used as molecular biomarkers. Through a de novo hepatic transcriptome assembly, using Illumina MiSeq, 54,285 sequences were assembled creating a reference transcriptome for this guppy species. Transcripts involved in biotransformation systems, antioxidant defenses, ABC transporters, nuclear and xenobiotic receptors were identified and evaluated by qPCR. Sanitary sewage induced transcriptional changes in AhR, PXR, CYP2K1, CYP3A30, NQO1, UGT1A1, GSTa3, GSTmu, ST1C1, SOD, ABCC1 and SOX9 genes from liver of fish, particularly after 96 h of exposure. Changes in hepatic enzyme activities were also observed. The enzymes showed differences in fish exposed to both periods, while in the gills there was a prevalence of significant results after 96 h. The observed differences were associated to gender and/or to sewage exposure. The obtained results support the use of P. vivipara as sentinel and model organism for ecotoxicological studies and evidence the importance of understanding the differential responses associated to gender.

N(alpha)-acetyltransferase 40-mediated histone acetylation plays an important role in ecdysone regulation of metamorphosis in the red flour beetle, Tribolium castaneum

Histone acetylation, a crucial epigenetic modification, is governed by histone acetyltransferases (HATs), that regulate many biological processes. Functions of HATs in insects are not well understood. We identified 27 HATs and determined their functions using RNA interference (RNAi) in the model insect, Tribolium castaneum. Among HATs studied, N-alpha-acetyltransferase 40 (NAA40) knockdown caused a severe phenotype of arrested larval development. The steroid hormone, ecdysone induced NAA40 expression through its receptor, EcR (ecdysone receptor). Interestingly, ecdysone-induced NAA40 regulates EcR expression. NAA40 acetylates histone H4 protein, associated with the promoters of ecdysone response genes: EcR, E74, E75, and HR3, and causes an increase in their expression. In the absence of ecdysone and NAA40, histone H4 methylation by arginine methyltransferase 1 (ART1) suppressed the above genes. However, elevated ecdysone levels at the end of the larval period induced NAA40, promoting histone H4 acetylation and increasing the expression of ecdysone response genes. NAA40 is also required for EcR, and steroid-receptor co-activator (SRC) mediated induction of E74, E75, and HR3. These findings highlight the key role of ecdysone-induced NAA40-mediated histone acetylation in the regulation of metamorphosis.

Transcriptomic and Metatranscriptomic Analyses Provide New Insights into the Response of the Pea Aphid Acyrthosiphon pisum (Hemiptera: Aphididae) to Acetamiprid

Acetamiprid is a broad-spectrum neonicotinoid insecticide used in agriculture to control aphids. While recent studies have documented resistance to acetamiprid in several aphid species, the underlying mechanisms are still not fully understood. In this study, we analyzed the transcriptome and metatranscriptome of a laboratory strain of the pea aphid, Acyrthosiphon pisum (Harris, 1776), with reduced susceptibility to acetamiprid after nine generations of exposure to identify candidate genes and the microbiome involved in the adaptation process. Sequencing of the transcriptome of both selected (RS) and non-selected (SS) strains allowed the identification of 14,858 genes and 4938 new transcripts. Most of the differentially expressed genes were associated with catalytic activities and metabolic pathways involving carbon and fatty acids. Specifically, alcohol-forming fatty acyl-CoA reductase (FAR) and acyl-CoA synthetase (ACSF2), both involved in the synthesis of epidermal wax layer components, were significantly upregulated in RS, suggesting that adaptation to acetamiprid involves the synthesis of a thicker protective layer. Metatranscriptomic analyses revealed subtle shifts in the microbiome of RS. These results contribute to a deeper understanding of acetamiprid adaptation by the pea aphid and provide new insights for aphid control strategies.

Transcriptome profiling and characterization of genes associated with tuberization under high temperature in aeroponics in potato cv. Kufri Anand

BACKGROUND

High temperature stress is an important abiotic factor, which affects tuberization and ultimately causes heavy yield reduction in potato.

OBJECTIVES

Identification and characterization of genes associated with tuberization under high temperature stress is essential for future management through biotechnology.

METHODOLOGY

Two contrasting potato varieties Kufri Anand (profuse tuber-bearing) versus Kufri Frysona (very less/scanty tuber-bearing, control) were cultivated in aeroponics under high temperature stress, and transcriptomes were analyzed.

RESULTS

Potato cv. Kufri Anand was found superior over control (Kufri Frysona) for tuber yield and its component traits along with root morphology under aeroponics. Transcriptomes of tuber and leaf tissues were analyzed. Statistically significant (p < 0.05) differentially expressed genes (DEGs) were categorised into up-regulated (> 2 log2 fold change, FC) and down-regulated (< -2 log2 FC) genes. DEGs were annotated by gene ontology and KEGG pathways. A few selected up-regulated genes of both tissues were identified, and phylogeny tree and motif analysis were analysed based on 36 peptide sequences representing 15 selected DEGs in this study. Further, gene expression markers were developed and validated by real time qPCR analysis for the identification of high temperature tolerant genotypes.

CONCLUSION

A few key genes associated in tuberization under high temperature conditions were heat shock proteins (e.g. 18.5 kDa class I heat shock protein), sugar metabolism (e.g. glucosyltransferase), transcription factor (e.g. WRKY), and phytohormones (e.g. auxin-induced beta-glucosidase). Our study provides an overview of key genes involved in tuberization under high temperature stress in potato cv. Kufri Anand under aeroponics.

Comparative Genomics Reveal Phylogenetic Relationship and Chromosomal Evolutionary Events of Eight Cervidae Species

Cervidae represents a family that is not only rich in species diversity but also exhibits a wide range of karyotypes. The controversies regarding the phylogeny and classification of Cervidae still persist. The flourishing development of the genomic era has made it possible to address these issues at the genomic level. Here, the genomes of nine species were used to explore the phylogeny and chromosomal evolutionary events of Cervidae. By conducting whole-genome comparisons, we identified single-copy orthologous genes across the nine species and constructed a phylogenetic tree based on the single-copy orthologous genes sequences, providing new insights into the phylogeny of Cervidae, particularly the phylogenetic relationship among sika deer, red deer, wapiti and Tarim red deer. Gene family analysis revealed contractions in the olfactory receptor gene family and expansions in the histone gene family across eight Cervidae species. Furthermore, synteny analysis was used to explore the chromosomal evolutionary events of Cervidae species, revealing six chromosomal fissions during the evolutionary process from Bovidae to Cervidae. Notably, specific chromosomal fusion events were found in four species of Cervus, and a unique chromosomal fusion event was identified in Muntiacus reevesi. Our study further completed the phylogenetic relationship within the Cervidae and demonstrated the feasibility of inferring species phylogeny at the whole-genome level. Additionally, our findings on gene family evolution and the chromosomal evolutionary events in eight Cervidae species lay a foundation for comprehensive research of the evolution of Cervidae.

De novo transcriptomic analysis of Doum Palm (Hyphaene compressa) revealed an insight into its potential drought tolerance

BACKGROUND

Doum palms (Hyphaene compressa) perform a crucial starring role in the lives of Kenya's arid and semi-arid people for empowerment and sustenance. Despite the crop's potential for economic gain, there is a lack of genetic resources and detailed information about its domestication at the molecular level. Given the doum palm's vast potential as a widely distributed plant in semi-arid and arid climates and a source of many applications, coupled with the current changing climate scenario, it is essential to understand the molecular processes that provide drought resistance to this plant.

RESULTS

Assembly of the first transcriptome of doum palms subjected to water stress generated about 39.97 Gb of RNA-Seq data. The assembled transcriptome revealed 193,167 unigenes with an average length of 1655 bp, with 128,708 (66.63%) successfully annotated in seven public databases. Unigenes exhibited significant differentially expressed genes (DEGs) in well-watered and stressed-treated plants, with 45071 and 42457 accounting for up-regulated and down-regulated DEGs, respectively. GO term, KEGG, and KOG analysis showed that DEGs were functionally enriched cellular processes, metabolic processes, cellular and catalytic activity, metabolism, genetic information processing, signal transduction mechanisms, and posttranslational modification pathways. Transcription factors (TF), such as the MYB, WRKY, NAC family, FAR1, B3, bHLH, and bZIP, were the prominent TF families identified as doum palm DEGs encoding drought stress tolerance.

CONCLUSIONS

This study provides a complete understanding of DEGs involved in drought stress at the transcriptome level in doum palms. This research is, therefore, the foundation for the characterization of potential genes, leading to a clear understanding of its drought stress responses and providing resources for improved genetic modification.

Transcriptome analysis of the bloom-forming dinoflagellate Prorocentrum donghaiense exposed to Ginkgo biloba leaf extract, with an emphasis on photosynthesis

Ginkgo biloba leaf extract (GBE) can effectively treat bloom-forming freshwater algae. However, there is limited information about the underlying suppression mechanism of the marine bloom-forming Prorocentrum donghaiense-the most dominant algal bloom species in the East China Sea. We investigated the effect of GBE on P. donghaiense in terms of its response to photosynthesis at the molecular/omic level. In total, 93,743 unigenes were annotated using six functional databases. Furthermore, 67,203 differentially expressed genes (DEGs) were identified in algae treated with 1.8 g∙L-1 GBE. Among these DEGs, we identified the genes involved in photosynthesis. PsbA, PsbB and PsbD in photosystem II, PsaA in photosystem I, and PetB and PetD in the cytochrome b6/f complex were downregulated. Other related genes, such as PsaC, PsaE, and PsaF in photosystem I; PetA in the cytochrome b6/f complex; and atpA, atpD, atpH, atpG, and atpE in the F-type H+-ATPase were upregulated. These results suggest that the structure and activity of the complexes were destroyed by GBE, thereby inhibiting the electron flow between the primary and secondary quinone electron acceptors, primary quinone electron acceptor, and oxygen-evolving complex in the PSII complex, and interrupting the electron flow between PSII and PSI, ultimately leading to a decline in algal cell photosynthesis. These findings provide a basis for understanding the molecular mechanisms underlying P. donghaiense exposure to GBE and a theoretical basis for the prevention and control of harmful algal blooms.

Comparative transcriptome analysis reveals potential regulatory mechanisms of genes and immune pathways following Vibrio harveyi infection in red drum (Sciaenops ocellatus)

Red drum (Sciaenops ocellatus), as an important economical marine fish, has been affected by various bacterial diseases in recent years. Vibrio harveyi cause fatal vibriosis in S. ocellatus, leading to massive mortality and causing significant setbacks in aquaculture. However, the regulatory mechanisms of S. ocellatus response to V. harveyi infection are poorly understood. In this regard, we performed transcriptomic analysis with head kidney tissues of S. ocellatus after V. harveyi infection from 12 h to 48 h to reveal genes, gene expression profiles, and pathways involved in immune and inflammation responses. Specifically, a total of 9,599, 5,728, and 7144 differentially expressed genes (DEGs) were identified after V. harveyi infection at 12 h, 24 h, and 48 h, respectively, and 1,848 shared DEGs have been identified from the above three comparison groups. Subsequent pathway analysis revealed that the shared DEGs following V. harveyi were involved in complement and coagulation cascades (C1R, C1QC, C3, C4, C5, C7, C8A, C8B, C8G, C9, CFB, CFH, and CFI), MAPK signaling pathway, chemokine signaling pathway (CCL19, CXCL8, CXCL12, CXCL14, CCR4, CCR7, and CXCR2), PPAR signaling pathway (PPAR-α, PPAR-γ and PPAR-β), and TNF signaling pathway. Finally, the expression patterns of DEGs in head kidney tissues and S. ocellatus macrophages were validated by qRT-PCR, suggesting the reliability of RNA sequencing for gene expression analysis. This dynamic transcriptome analyses provided insights into gene expression regulation and immune related pathways involved in S. ocellatus after V. harveyi infection, and provides useful information for further study on the immune defense mechanisms in S. ocellatus as well as other teleost species.

Multi-omics revealed the mechanisms of Codonopsis pilosula aqueous extract in improving UC through blocking abnormal activation of PI3K/Akt signaling pathway

ETHNOPHARMACOLOGICALRELEVANCE

Codonopsis pilosula (DS), a traditional Chinese medicine, had been used to regulate the immune, digestive and circulatory systems of the human, as well as protect the gastrointestinal tract, improve lung function.

AIM OF THE STUDY

The aim of study was to explore the effects and mechanism of Codonopsis pilosula aqueous extract (DS) intervention in improving ulcerative colitis (UC).

MATERIALS AND METHODS

UC model rats were established using combination of TNBS and ethanol. Tissue samples were collected for transcriptome and metabolomics analysis. Network pharmacology was performed on DS to identify bioactive compounds. Western blot was used to detect the key proteins involved in UC pathogenesis and PI3K/AKT pathways.

RESULTS

DS exerted the preventive and therapeutic effects in improving UC via inhibiting abnormal inflammatory responses and promoting antioxidant capacity. Levels of intestinal barrier, oxidative stress and inflammatory mediators were improved to nearly normal level in vivo by DS. Metabolome profiles showed that DS could restore the metabolic disorders associated with the UC pathogenesis. Further transcriptome results showed that DS mainly alleviate UC through inhibiting PI3K/Akt signaling pathway, and various related genes that dramatically expressed in UC Model rats were downregulated by DS. Typically, network pharmacology analysis identified that Glycitein was the hub compounds that involved in the mechanism of DS in improving UC.

CONCLUSIONS

The results show that Codonopsis pilosula (DS) was an potential excellent material in treating of UC depending on its suitable concentration. Possible therapeutic mechanisms of the DS involved in mitigating colonal inflammation, restoring metabolic disorders, promoting antioxidant capacity, and especially blocking the activation of PI3K/Akt pathway.

Systematic assessment of the antifungal mechanism of soil fumigant methyl isothiocyanate against Fusarium oxysporum

Fusarium oxysporum is an important phytopathogenic fungus, it can be controlled by the soil fumigant methyl isothiocyanate (MITC). However, the antimicrobial mechanism of MITC against F. oxysporum, especially at the transcriptional level, is still unclear. In this experiment, the antimicrobial mechanism of MITC against F. oxysporum was investigated. Our results indicated that when F. oxysporum was exposed to 6 mg/L MITC for 12 h, the inhibitory rate of MITC on F. oxysporum was 80%. Transmission electron microscopes showed that the cell wall and membrane of F. oxysporum had shrunk and folded, vacuoles increased, and mitochondria swelled and deformed. In addition, the enzyme activity of F. oxysporum treated with MITC showed a decrease of 32.50%, 8.28% and 74.04% in catalase, peroxidase and superoxide dismutase, respectively. Transcriptome sequencing of F. oxysporum was performed and the results showed that 1478 differentially expressed genes (DEGs) were produced in response to MITC exposure. GO and KEGG analysis showed that the DEGs identified were involved in substance and energy metabolism, signal transduction, transport and catalysis. MITC disrupted cell homeostasis by influencing the expression of some key genes involved in chitin synthase and detoxification enzymes production, but F. oxysporum also protected itself by up-regulating genes involved in energy synthesis (such as upregulating acnA, CS and LSC2 in TCA). qRT-PCR data validated the reliability of transcriptome data. Our research used biochemical and genetic techniques to identify molecular lesions in the mycelia of F. oxysporum exposed to MITC, and provide valuable insights into the toxic mechanism of pathogenic fungi mediated by MITC. These techniques are also likely to be useful for rapidly screening and identifying new, environmentally-friendly soil fumigants that are efficacious against fungal pathogens.

The Features of Shared Genes among Transcriptomes Probed in Atopic Dermatitis, Psoriasis, and Inflammatory Acne: S100A9 Selection as the Target Gene

BACKGROUND

Atopic dermatitis (AD), psoriasis (PS), and inflammatory acne (IA) are well-known as inflammatory skin diseases. Studies of the transcriptome with altered expression levels have reported a large number of dysregulated genes and gene clusters, particularly those involved in inflammatory skin diseases.

OBJECTIVE

To identify genes commonly shared in AD, PS, and IA that are potential therapeutic targets, we have identified consistently dysregulated genes and disease modules that overlap with AD, PS, and IA.

METHODS

Microarray data from AD, PS, and IA patients were downloaded from Gene Expression Omnibus (GEO), and identification of differentially expressed genes from microarrays of AD, PS, and IA was conducted. Subsequently, gene ontology and gene set enrichment analysis, detection of disease modules with known disease-associated genes, construction of the protein-protein interaction (PPI) network, and PPI sub-mapping analysis of shared genes were performed. Finally, the computational docking simulations between the selected target gene and inhibitors were conducted.

RESULTS

We identified 50 shared genes (36 up-regulated and 14 down-regulated) and disease modules for each disease. Among the shared genes, 20 common genes in PPI network were detected such as LCK, DLGAP5, SELL, CEP55, CDC20, RRM2, S100A7, S100A9, MCM10, AURKA, CCNB1, CHEK1, BTC, IL1F7, AGTR1, HABP4, SERPINB13, RPS6KA4, GZMB, and TRIP13. Finally, S100A9 was selected as the target gene for therapeutics. Docking simulations between S100A9 and known inhibitors indicated several key binding residues, and based on this result, we suggested several cannabinoids such as WIN-55212-2, JZL184, GP1a, Nabilone, Ajulemic acid, and JWH-122 could be potential candidates for a clinical study for AD, PS, and IA via inhibition of S100A9-related pathway.

CONCLUSION

Overall, our approach may become an effective strategy for discovering new disease candidate genes for inflammatory skin diseases with a reevaluation of clinical data.

Genetically modified soybean lines exhibit less transcriptomic variation compared to natural varieties

Genetically modified (GM) soybeans provide a huge amount of food for human consumption and animal feed. However, the possibility of unexpected effects of transgenesis has increased food safety concerns. High-throughput sequencing profiling provides a potential approach to directly evaluate unintended effects caused by foreign genes. In this study, we performed transcriptomic analyses to evaluate differentially expressed genes (DEGs) in individual soybean tissues, including cotyledon (C), germ (G), hypocotyl (H), and radicle (R), instead of using the whole seed, from four GM and three non-GM soybean lines. A total of 3,351 DEGs were identified among the three non-GM soybean lines. When the GM lines were compared with their non-GM parents, 1,836 to 4,551 DEGs were identified. Furthermore, Gene Ontology (GO) analysis of the DEGs showed more abundant categories of GO items (199) among non-GM lines than between GM lines and the non-GM natural varieties (166). Results of Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that most KEGG pathways were the same for the two types of comparisons. The study successfully employed RNA sequencing to assess the differences in gene expression among four tissues of seven soybean varieties, and the results suggest that transgenes do not induce massive transcriptomic alterations in transgenic soybeans compared with those that exist among natural varieties. This work offers empirical evidence to investigate the genomic-level disparities induced by genetic modification in soybeans, specifically focusing on seed tissues.

Comparative Transcriptome Analysis of the Skin and the Peritoneal Wall Layer of Triplophysa stenura Distributed in High Elevations

A total of 81,868 All-Unigenes were sequenced and assembled by the transcriptome in the dorsal skin, the lateral skin, and the peritoneal wall layer of Triplophysa stenura with a total assembly length of 123,827,585 bp, and 68,750 unigenes were annotated to seven functional databases. A total of 588 DEGs were screened between the dorsal and lateral skin, 17,097 DEGs were screened between the dorsal skin and the peritoneal wall layer, and 16,598 DEGs were screened between the lateral skin and the peritoneal wall layer. Most of DEGs in three tissues were annotated to GO terms related to cellular structures, binding, cellular processes, and catalytic activity. They were also annotated to KEGG pathways such as the MAPK signaling pathway, PI3K-Akt signaling pathway, Wnt signaling pathway, melanogenesis, tyrosine metabolism, and cell cycle. A total of twenty-three DEGs were found to be enriched in the melanin synthesis pathway by a local Blast comparison, of which nine DEGs were significantly upregulated in the peritoneal wall layer and six DEGs were significantly upregulated in the dorsal and lateral skin. The results suggest that these genes may be associated with the molecular mechanism of melanin synthesis in T. stenura, and the differential regulation of genes may be related to the differences in UVR intensity and tissue sites of melanin synthesis. Further investigation is needed on how these genes specifically regulate melanin synthesis.

A 4D Proteome Investigation of the Potential Mechanisms of SA in Triggering Resistance in Kiwifruit to Pseudomonas syringae pv. actinidiae

Kiwifruit bacterial cankers caused by Pseudomonas syringae pv. actinidiae (Psa) are a serious threat to the kiwifruit industry. Salicylic acid (SA) regulates plant defense responses and was previously found to enhance kiwifruit's resistance to Psa. However, the underlying mechanisms of this process remain unclear. In this study, we used 4D proteomics to investigate how SA enhances kiwifruit's resistance to Psa and found that both SA treatment and Psa infection induced dramatic changes in the proteomic pattern of kiwifruit. Psa infection triggered the activation of numerous resistance events, including the MAPK cascade, phenylpropanoid biosynthesis, and hormone signaling transduction. In most cases, the differential expression of a number of genes involved in the SA signaling pathway played a significant role in kiwifruit's responses to Psa. Moreover, SA treatment upregulated numerous resistance-related proteins, which functioned in defense responses to Psa, including phenylpropanoid biosynthesis, the MAPK cascade, and the upregulation of pathogenesis-related proteins. We also found that SA treatment could facilitate timely defense responses to Psa infection and enhance the activation of defense responses that were downregulated in kiwifruit during infection with Psa. Thus, our research deciphered the potential mechanisms of SA in promoting Psa resistance in kiwifruit and can provide a basis for the use of SA to enhance kiwifruit resistance and effectively control the occurrence of kiwifruit bacterial cankers.

Dispersal from the Qinghai-Tibet plateau by a high-altitude butterfly is associated with rapid expansion and reorganization of its genome

Parnassius glacialis is a typical "Out of the QTP" alpine butterfly that originated on the Qinghai-Tibet Plateau (QTP) and dispersed into relatively low-altitude mountainous. Here we assemble a chromosome-level genome of P. glacialis and resequence 9 populations in order to explore the genome evolution and local adaptation of this species. These results indicated that the rapid accumulation and slow unequal recombination of transposable elements (TEs) contributed to the formation of its large genome. Several ribosomal gene families showed extensive expansion and selective evolution through transposon-mediated processed pseudogenes. Additionally, massive structural variations (SVs) of TEs affected the genetic differentiation of low-altitude populations. These low-altitude populations might have experienced a genetic bottleneck in the past and harbor genes with selective signatures which may be responsible for the potential adaptation to low-altitude environments. These results provide a foundation for understanding genome evolution and local adaptation for "Out of the QTP" of P. glacialis.

Mapping Chromatin Occupancy of Ppp1r1b-lncRNA Genome-Wide Using Chromatin Isolation by RNA Purification (ChIRP)-seq

Long non-coding RNA (lncRNA) mediated transcriptional regulation is increasingly recognized as an important gene regulatory mechanism during development and disease. LncRNAs are emerging as critical regulators of chromatin state; yet the nature and the extent of their interactions with chromatin remain to be fully revealed. We have previously identified Ppp1r1b-lncRNA as an essential epigenetic regulator of myogenic differentiation in cardiac and skeletal myocytes in mice and humans. We further demonstrated that Ppp1r1b-lncRNA function is mediated by the interaction with the chromatin-modifying complex polycomb repressive complex 2 (PRC2) at the promoter of myogenic differentiation transcription factors, TBX5 and MyoD1. Herein, we employed unbiased chromatin isolation by RNA purification (ChIRP) and high throughput sequencing to map the repertoire of Ppp1r1b-lncRNA chromatin occupancy genome-wide in the mouse muscle myoblast cell line. We uncovered a total of 99732 true peaks corresponding to Ppp1r1b-lncRNA binding sites at high confidence (p-value < 1E-5) and enrichment score ≥ 10). The Ppp1r1b-lncRNA-binding sites averaged 558 bp in length and were distributed widely within the coding and non-coding regions of the genome. Approximately 46% of these true peaks were mapped to gene elements, of which 1180 were mapped to experimentally validated promoter sequences. Importantly, the promoter-mapped binding sites were enriched in myogenic transcription factors and heart development while exhibiting focal interactions with known motifs of proximal promoters and transcription initiation by RNA Pol-II, including TATA-box, transcription initiator motif, CCAAT-box, and GC-box, supporting Ppp1r1b-lncRNA role in transcription initiation of myogenic regulators. Remarkably, nearly 40% of Ppp1r1b-lncRNA-binding sites mapped to gene introns were enriched with the Homeobox family of transcription factors and exhibited TA-rich motif sequences, suggesting potential motif-specific Ppp1r1b-lncRNA-bound introns. Lastly, more than 136521 enhancer sequences were detected in Ppp1r1b-lncRNA-occupancy sites at high confidence. Among these enhancers, 3390 (12%) exhibited cell type/tissue-specific enrichment in fetal heart and muscles. Together, our findings provide further insights into the genome-wide Ppp1r1b-lncRNA: Chromatin interactome that may dictate its function in myogenic differentiation and potentially other cellular and biological processes.

siRNA Silencing of FpVtg Induces Ovarian Cell Apoptosis in Redtail Prawn, Fenneropenaeus penicillatus

Inadequate gonadal maturation and poor spawning performance increasingly threaten the sustainability of shrimp aquaculture. Unraveling the mechanisms regulating ovarian development and maturation hence is critical to address industry challenges. Vitellogenin (Vtg), a precursor of yolk protein found in the hepatopancreas and ovary of shrimp, plays a key role in facilitating shrimp's oocyte maturation and embryonic development after oviposition. This study found that FpVtg was specifically expressed in F. penicillatus hepatopancreas and ovary. FpVtg was localized predominantly in the oocyte cytoplasm and distributed uniformly in the hepatopancreas tissue. Silencing FpVtg led to apoptosis in both hepatopancreas and ovary tissues. Furthermore, FpVtg depletion upregulated the expression of ovarian peritrophin 1, ovarian peritrophin 2, serine proteinase inhibitor 6, and juvenile hormone esterase-like carboxylesterase 1, while downregulated that of vitellogenin, delta-9 desaturase, and insulin-like receptor. KEGG pathway analysis implicated such as PI3K-AKT signaling, RNA transport, ECM-receptor interaction, hippo signaling, oocyte meiosis, and apoptosis were enriched and involved in ovarian development. These findings have provided insights into the FpVtg's reproductive role and the associated regulatory genes and pathways in F. penicillatus. This knowledge can contribute to establishing strategies to improve the breeding and aquaculture production of F. penicillatus by elucidating its vitellogenesis regulation in redtail prawn and other penaeid species. Further characterization of the implicated pathways and genes will clarify the intricacies underlying ovarian maturation.

Mapping Chromatin Occupancy of Ppp1r1b-lncRNA Genome-Wide Using Chromatin Isolation by RNA Purification (ChIRP)-seq

Long non-coding RNA (lncRNA) mediated transcriptional regulation is increasingly recognized as an important gene regulatory mechanism during development and disease. LncRNAs are emerging as critical regulators of chromatin state; yet the nature and the extent of their interactions with chromatin remain to be fully revealed. We have previously identified Ppp1r1b-lncRNA as an essential epigenetic regulator of myogenic differentiation in cardiac and skeletal myocytes in mice and humans. We further demonstrated that Ppp1r1b-lncRNA function is mediated by the interaction with the chromatin-modifying complex polycomb repressive complex 2 (PRC2) at the promoter of myogenic differentiation transcription factors, TBX5 and MyoD1. Herein, we employed an unbiased chromatin isolation by RNA purification (ChIRP) and high throughput sequencing to map the repertoire of Ppp1r1b-lncRNA chromatin occupancy genome-wide in the mouse muscle myoblast cell line. We uncovered a total of 99732 true peaks corresponding to Ppp1r1b-lncRNA binding sites at high confidence (P-value < 1e-5 and enrichment score ≥ 10). The Ppp1r1b-lncRNA-binding sites averaged 558 bp in length and were distributed widely within the coding and non-coding regions of the genome. Approximately 46% of these true peaks were mapped to gene elements, of which 1180 were mapped to experimentally validated promoter sequences. Importantly, the promoter-mapped binding sites were enriched in myogenic transcription factors and heart development while exhibiting focal interactions with known motifs of proximal promoters and transcription initiation by RNA polII, including TATA, transcription initiator, CCAAT-box, and GC-box, supporting Ppp1r1b-lncRNA role in transcription initiation of myogenic regulators. Remarkably, nearly 40% of Ppp1r1b-lncRNA-binding sites mapped to gene introns, were enriched with the Homeobox family of transcription factors, and exhibited TA-rich motif sequences, suggesting potential motif specific Ppp1r1b-lncRNA-bound introns. Lastly, more than 136521enhancer sequences were detected in Ppp1r1b-lncRNA-occupancy sites at high confidence. Among these enhancers,12% exhibited cell type/tissue-specific enrichment in fetal heart and muscles. Together, our findings provide further insights into the genome-wide Ppp1r1b-lncRNA: Chromatin interactome that may potentially dictate its function in myogenic differentiation and potentially other cellular and biological processes.

Chemosensory proteins as putative semiochemical carriers in the desert isopod Hemilepistus reaumurii

The order Isopoda contains both aquatic and terrestrial species, among which Hemilepistus reaumurii, which lives in arid environments and is the most adapted to terrestrial life. Olfaction has been deeply investigated in insects while it has received very limited attention in other arthropods, particularly in terrestrial crustaceans. In insects, soluble proteins belonging to two main families, Odorant Binding Proteins (OBPs) and Chemosensory Proteins (CSPs), are contained in the olfactory sensillar lymph and are suggested to act as carriers of hydrophobic semiochemicals to or from membrane-bound olfactory receptors. Other protein families, namely Nieman-Pick type 2 (NPC2) and Lipocalins (LCNs) have been also reported as putative odorant carriers in insects and other arthropod clades. In this study, we have sequenced and analysed the transcriptomes of antennae and of the first pair of legs of H. reaumurii focusing on soluble olfactory proteins. Interestingly, we have found 13 genes encoding CSPs, whose sequences differ from those of the other arthropod clades, including non-isopod crustaceans, for the presence of two additional cysteine residues, besides the four conserved ones. Binding assays on two of these proteins showed strong affinities for fatty acids and long-chain unsaturated esters and aldehydes, putative semiochemicals for this species.

Arbuscular mycorrhizal fungi heterokaryons have two nuclear populations with distinct roles in host-plant interactions

Arbuscular mycorrhizal fungi (AMF) are prominent root symbionts that can carry thousands of nuclei deriving from two parental strains in a large syncytium. These co-existing genomes can also vary in abundance with changing environmental conditions. Here we assemble the nuclear genomes of all four publicly available AMF heterokaryons using PacBio high-fidelity and Hi-C sequencing. We find that the two co-existing genomes of these strains are phylogenetically related but differ in structure, content and epigenetics. We confirm that AMF heterokaryon genomes vary in relative abundance across conditions and show this can lead to nucleus-specific differences in expression during interactions with plants. Population analyses also reveal signatures of genetic exchange indicative of past events of sexual reproduction in these strains. This work uncovers the origin and contribution of two nuclear genomes in AMF heterokaryons and opens avenues for the improvement and environmental application of these strains.

DNA methylation mediates overgrazing-induced clonal transgenerational plasticity

Overgrazing generally induces dwarfism in grassland plants, and these phenotypic traits could be transmitted to clonal offspring even when overgrazing is excluded. However, the dwarfism-transmitted mechanism remains largely unknown, despite generally thought to be enabled by epigenetic modification. To clarify the potential role of DNA methylation on clonal transgenerational effects, we conducted a greenhouse experiment with Leymus chinensis clonal offspring from different cattle/sheep overgrazing histories via the demethylating agent 5-azacytidine. The results showed that clonal offspring from overgrazed (by cattle or sheep) parents were dwarfed and the auxin content of leaves significantly decreased compared to offspring from no-grazed parents'. The 5-azaC application generally increased the auxin content and promoted the growth of overgrazed offspring while inhibited no-grazed offspring growth. Meanwhile, there were similar trends in the expression level of genes related to auxin-responsive target genes (ARF7, ARF19), and signal transduction gene (AZF2). These results suggest that DNA methylation leads to overgrazing-induced plant transgenerational dwarfism via inhibiting auxin signal pathway.

Pseudocercospora fijiensis Conidial Germination Is Dominated by Pathogenicity Factors and Effectors

Conidia play a vital role in the survival and rapid spread of fungi. Many biological processes of conidia, such as adhesion, signal transduction, the regulation of oxidative stress, and autophagy, have been well studied. In contrast, the contribution of pathogenicity factors during the development of conidia in fungal phytopathogens has been poorly investigated. To date, few reports have centered on the pathogenicity functions of fungal phytopathogen conidia. Pseudocercospora fijiensis is a hemibiotrophic fungus and the causal agent of the black Sigatoka disease in bananas and plantains. Here, a conidial transcriptome of P. fijiensis was characterized computationally. Carbohydrates, amino acids, and lipid metabolisms presented the highest number of annotations in Gene Ontology. Common conidial functions were found, but interestingly, pathogenicity factors and effectors were also identified. Upon analysis of the resulting proteins against the Pathogen-Host Interaction (PHI) database, 754 hits were identified. WideEffHunter and EffHunter effector predictors identified 618 effectors, 265 of them were shared with the PHI database. A total of 1107 conidial functions devoted to pathogenesis were found after our analysis. Regarding the conidial effectorome, it was found to comprise 40 canonical and 578 non-canonical effectors. Effectorome characterization revealed that RXLR, LysM, and Y/F/WxC are the largest effector families in the P. fijiensis conidial effectorome. Gene Ontology classification suggests that they are involved in many biological processes and metabolisms, expanding our current knowledge of fungal effectors.

Arbuscular mycorrhizal fungi improve selenium uptake by modulating root transcriptome of rice (Oryza sativa L.)

Although selenium (Se) is an essential trace element in humans, the intake of Se from food is still generally inadequate throughout the world. Inoculation with arbuscular mycorrhizal fungi (AMF) improves the uptake of Se in rice (Oryza sativa L.). However, the mechanism by which AMF improves the uptake of Se in rice at the transcriptome level is unknown. Only a few studies have evaluated the effects of uptake of other elements in rice under the combined effects of Se and AMF. In this study, Se combined with the AMF Funneliformis mosseae (Fm) increased the biomass and Se concentration of rice plants, altered the pattern of ionomics of the rice roots and shoots, and reduced the antagonistic uptake of Se with nickel, molybdenum, phosphorus, and copper compared with the treatment of Se alone, indicating that Fm can enhance the effect of fertilizers rich in Se. Furthermore, a weighted gene co-expression network analysis (WGCNA) showed that the hub genes in modules significantly associated with the genes that contained Se and were related to protein phosphorylation, protein serine/threonine kinase activity, membrane translocation, and metal ion binding, suggesting that the uptake of Se by the rice roots may be associated with these genes when Fm and Se act in concert. This study provides a reference for the further exploration of genes related to Se uptake in rice under Fm treatment.

Decoding the mechanism of andrographolide to combat hepatocellular carcinoma: a network pharmacology integrated molecular docking and dynamics approach

HepatoCellular Carcinoma, being one of the most mortally convoluted malignancy with mounting number of occurrences across the world and being classified as the third most prevalent cause of cancer-associated mortalities and sixth most prevalent neoplasia. The active phytoconstituent andrographolide, derived from Andrographis paniculata is conveyed to reconcile a number of human ailments including various oncologies. However, the molecular mechanism underlying the anti-oncogenic effects of Andrographolide on HCC remains skeptical and unclear, emerging as a budding challenge for researchers and oncologists. The present study intends to analyze the underlying pharmacological mechanism of Andrographolide over HCC, established via assimilated approach of network pharmacology. Herein, the Network pharmacology stratagem was instigated to investigate potential HCC targets. The Andrographolide targets along with HCC targets were extracted from multiple databases. A total of 162 potential overlapping targets among HCC and Andrographolide were obtained and further subjected to gene ontology and Pathway enrichment analysis by employing OmicsBox and DAVID database, respectively. Subsequently, Protein-protein interaction network construction by Cytoscape software identified the top 10 hub nodes which were validated by survival and expression analysis. Further, the results derived from molecular docking and dynamic simulations by CB-Dock2 server and Desmond module (Schrodinger software) indicate ALB, CCND1, HIF1A, TNF, and VEGFA as potential Andrographolide related targets with high binding affinity and promising complex stability. Our findings not only reveal the antioncogenic role of andrographolide but also provide novel insights illuminating the identified targets as scientific foundation for anti-oncogenic clinical application of andrographolide in HCC therapeutics.Communicated by Ramaswamy H. Sarma.

Exogenous IAA application affects the specific characteristics of fluoranthene distribution in Arabidopsis

Indole-3-acetic acid (IAA) is a crucial growth regulator involved in the accumulation of polycyclic aromatic hydrocarbons (PAHs). However, the precise physiological and molecular mechanisms underlying IAA-mediated plant growth and PAH accumulation are not yet fully understood. In this study, two distinct IAA-sensitive genotypes of Arabidopsis thaliana (wild type and Axr5 mutant) were chosen to investigate the mechanisms of fluoranthene (Flu) uptake and accumulation in plant tissues (roots and leaves) through physiological and molecular analyses. The results revealed that the Flu concentration in Axr5 leaves was significantly higher than that in wild-type (WT) leaves. In roots, the Flu content decreased significantly with increasing IAA treatment, while no significant changes were observed with lower IAA treatment. Principal component analysis demonstrated that Flu accumulation in Arabidopsis roots was associated with IAA concentrations, whereas Flu accumulation in leaves was dependent on the genotype. Moreover, Flu accumulation showed a positive correlation with the activity of glutathione S-transferase (GST) and root length and a positive correlation with catalase (CAT) and peroxidase (POD) activity in the leaves. Transcriptome analysis confirmed that the expression of the ethylene-related gene ATERF6 and GST-related genes ATGSTF14 and ATGSTU27 in roots, as well as the POD-related genes AtPRX9 and AtPRX25 and CAT-related gene AtCAT3 in leaves, played a role in Flu accumulation. Furthermore, WRKY transcription factors (TFs) in roots and NAC TFs in leaves were identified as important regulators of Flu accumulation. Understanding the mechanisms of Flu uptake and accumulation in A. thaliana provides valuable insights for regulating PAH accumulation in plants.

Transcriptomic and Physiological Analysis Reveals Genes Associated with Drought Stress Responses in Populus alba × Populus glandulosa

Drought stress affects plant productivity by altering plant responses at the morphological, physiological, and molecular levels. In this study, we identified physiological and genetic responses in Populus alba × Populus glandulosa hybrid clones 72-30 and 72-31 after 12 days of exposure to drought treatment. After 12 days of drought treatment, glucose, fructose, and sucrose levels were significantly increased in clone 72-30 under drought stress. The Fv/Fo and Fv/Fm values in both clones also decreased under drought stress. The changes in proline, malondialdehyde, and H2O2 levels were significant and more pronounced in clone 72-30 than in clone 72-31. The activities of antioxidant-related enzymes, such as catalase and ascorbate peroxidase, were significantly higher in the 72-31 clone. To identify drought-related genes, we conducted a transcriptomic analysis in P. alba × P. glandulosa leaves exposed to drought stress. We found 883 up-regulated and 305 down-regulated genes in the 72-30 clone and 279 and 303 in the 72-31 clone, respectively. These differentially expressed genes were mainly in synthetic pathways related to proline, abscisic acid, and antioxidants. Overall, clone 72-31 showed better drought tolerance than clone 72-30 under drought stress, and genetic changes also showed different patterns.

Development of genomic resources for Rhodes grass (Chloris gayana), draft genome and annotated variant discovery

Genomic resources for grasses, especially warm-season grasses are limited despite their commercial and environmental importance. Here, we report the first annotated draft whole genome sequence for diploid Rhodes grass (Chloris gayana), a tropical C4 species. Generated using long read nanopore sequencing and assembled using the Flye software package, the assembled genome is 603 Mbp in size and comprises 5,233 fragments that were annotated using the GenSas pipeline. The annotated genome has 46,087 predicted genes corresponding to 92.0% of the expected genomic content present via BUSCO analysis. Gene ontology terms and repetitive elements are identified and discussed. An additional 94 individual plant genotypes originating from three diploid and two tetraploid Rhodes grass cultivars were short-read whole genome resequenced (WGR) to generate a single nucleotide polymorphism (SNP) resource for the species that can be used to elucidate inter- and intra-cultivar relationships across both ploidy levels. A total of 75,777 high quality SNPs were used to generate a phylogenetic tree, highlighting the diversity present within the cultivars which agreed with the known breeding history. Differentiation was observed between diploid and tetraploid cultivars. The WGR data were also used to provide insights into the nature and evolution of the tetraploid status of the species, with results largely agreeing with the published literature that the tetraploids are autotetraploid.

Induced expression modes of genes related to Toll, Imd, and JAK/STAT signaling pathway-mediated immune response in Spodoptera frugiperda infected with Beauveria bassiana

Spodoptera frugiperda is one of the most harmful pests that attack maize and other major food crops and causes huge economic loss every year in China and other countries and regions. Beauveria bassiana, a kind of entomological fungus that is highly pathogenic to pests, is harmless to the environment and human beings. However, at present, S. frugiperda has gradually developed resistance to many pesticides and microbial insecticides. In this study, transcriptome sequencing was conducted to analyze the differences in gene expression between B. bassiana-infected and -uninfected S. frugiperda. More than 160 Gb of clean data were obtained as 150-bp paired-end reads using the Illumina HiSeq™ 4000 platform, and 2,767 and 2,892 DEGs were identified in LH36vsCK36 and LH144vsCK144, respectively. In order to explore the roles of JAK/STAT, Toll, and Imd signaling pathways in antifungal immune response in S. frugiperda against B. bassiana infection, the expression patterns of those signaling pathway-related genes in B. bassiana-infected S. frugiperda were analyzed by quantitative real-time PCR. In addition, antifungal activity experiments revealed that the suppression of JAK/STAT, Toll, and Imd signaling pathways by inhibitors could inhibit the antifungal activity to a large extent and lead to increased sensitivity of S. frugiperda to B. bassiana infection, indicating that JAK/STAT, Toll, and Imd signaling pathways and their associated genes might be involved in the synthesis and secretion of antifungal substances. This study implied that JAK/STAT, Toll, and Imd signaling pathways played crucial roles in the antifungal immune response of the S. frugiperda larvae, in which the related genes of these signaling pathways could play special regulatory roles in signal transduction. This study would improve our understanding of the molecular mechanisms underlying innate immunity and provide the basis for a wide spectrum of strategies against antifungal resistance of S. frugiperda.

Characterizing the circadian connectome of Ocimum tenuiflorum using an integrated network theoretic framework

Across the three domains of life, circadian clock is known to regulate vital physiological processes, like, growth, development, defence etc. by anticipating environmental cues. In this work, we report an integrated network theoretic methodology comprising of random walk with restart and graphlet degree vectors to characterize genome wide core circadian clock and clock associated raw candidate proteins in a plant for which protein interaction information is available. As a case study, we have implemented this framework in Ocimum tenuiflorum (Tulsi); one of the most valuable medicinal plants that has been utilized since ancient times in the management of a large number of diseases. For that, 24 core clock (CC) proteins were mined in 56 template plant genomes to build their hidden Markov models (HMMs). These HMMs were then used to identify 24 core clock proteins in O. tenuiflorum. The local topology of the interologous Tulsi protein interaction network was explored to predict the CC associated raw candidate proteins. Statistical and biological significance of the raw candidates was determined using permutation and enrichment tests. A total of 66 putative CC associated proteins were identified and their functional annotation was performed.

Transcriptome analysis provides insights into the stress response in cultivated peanut (Arachis hypogaea L.) subjected to drought-stress

BACKGROUND

Peanut (Arachis hypogaea L.) is one of the valuable oilseed crops grown in drought-prone areas worldwide. Drought severely limits peanut production and productivity significantly.

METHOD AND RESULTS

In order to decipher the drought tolerance mechanism in peanut under drought stress, RNA sequencing was performed in TAG - 24 (drought tolerant genotype) and JL-24 (drought susceptible genotype). Approximately 51 million raw reads were generated from four different libraries of two genotypes subjected to drought stress exerted by 20% PEG 6000 stress and control conditions, of which ~ 41 million (80.87%) filtered reads were mapped to the Arachis hypogaea L. reference genome. The transcriptome analysis detected 1,629 differentially expressed genes (DEGs), 186 genes encoding transcription factors (TFs) and 30,199 SSR among the identified DEGs. Among the differentially expressed TF encoding genes, the highest number of genes were WRKY followed by bZIP, C2H2, and MYB during drought stress. The comparative analysis between the two genotypes revealed that TAG-24 exhibits activation of certain key genes and transcriptional factors that are involved in essential biological processes. Specifically, TAG-24 showed activation of genes involved in the plant hormone signaling pathway such as PYL9, Auxin response receptor gene, and ABA. Additionally, genes related to water deprivation such as LEA protein and those involved in combating oxidative damage such as Glutathione reductase were also found to be activated in TAG-24.

CONCLUSION

This genome-wide transcription map, therefore, provides a valuable tool for future transcript profiling under drought stress and enriches the genetic resources available for this important oilseed crop.

Integrating Full-Length Transcriptome and RNA Sequencing of Siberian Wildrye (Elymus sibiricus) to Reveal Molecular Mechanisms in Response to Drought Stress

Drought is one of the most significant limiting factors affecting plant growth and development on the Qinghai-Tibet Plateau (QTP). Mining the drought-tolerant genes of the endemic perennial grass of the QTP, Siberian wildrye (Elymus sibiricus), is of great significance to creating new drought-resistant varieties which can be used in the development of grassland livestock and restoring natural grassland projects in the QTP. To investigate the transcriptomic responsiveness of E. sibiricus to drought stress, PEG-induced short- and long-term drought stress was applied to two Siberian wildrye genotypes (drought-tolerant and drought-sensitive accessions), followed by third- and second-generation transcriptome sequencing analysis. A total of 40,708 isoforms were detected, of which 10,659 differentially expressed genes (DEGs) were common to both genotypes. There were 2107 and 2498 unique DEGs in the drought-tolerant and drought-sensitive genotypes, respectively. Additionally, 2798 and 1850 DEGs were identified in the drought-tolerant genotype only under short- and long-term conditions, respectively. DEGs numbering 1641 and 1330 were identified in the drought-sensitive genotype only under short- and long-term conditions, respectively. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed that all the DEGs responding to drought stress in E. sibiricus were mainly associated with the mitogen-activated protein kinase (MAKP) signaling pathway, plant hormone signal transduction, the linoleic acid metabolism pathway, the ribosome pathway, and plant circadian rhythms. In addition, Nitrate transporter 1/Peptide transporter family protein 3.1 (NPF3.1) and Auxin/Indole-3-Acetic Acid (Aux/IAA) family protein 31(IAA31) also played an important role in helping E. sibiricus resist drought. This study used transcriptomics to investigate how E. sibiricus responds to drought stress, and may provide genetic resources and references for research into the molecular mechanisms of drought resistance in native perennial grasses and for breeding drought-tolerant varieties.

Comparative evolutionary analyses of eight whitefly Bemisia tabaci sensu lato genomes: cryptic species, agricultural pests and plant-virus vectors

BACKGROUND

The group of > 40 cryptic whitefly species called Bemisia tabaci sensu lato are amongst the world's worst agricultural pests and plant-virus vectors. Outbreaks of B. tabaci s.l. and the associated plant-virus diseases continue to contribute to global food insecurity and social instability, particularly in sub-Saharan Africa and Asia. Published B. tabaci s.l. genomes have limited use for studying African cassava B. tabaci SSA1 species, due to the high genetic divergences between them. Genomic annotations presented here were performed using the 'Ensembl gene annotation system', to ensure that comparative analyses and conclusions reflect biological differences, as opposed to arising from different methodologies underpinning transcript model identification.

RESULTS

We present here six new B. tabaci s.l. genomes from Africa and Asia, and two re-annotated previously published genomes, to provide evolutionary insights into these globally distributed pests. Genome sizes ranged between 616-658 Mb and exhibited some of the highest coverage of transposable elements reported within Arthropoda. Many fewer total protein coding genes (PCG) were recovered compared to the previously published B. tabaci s.l. genomes and structural annotations generated via the uniform methodology strongly supported a repertoire of between 12.8-13.2 × 10³ PCG. An integrative systematics approach incorporating phylogenomic analysis of nuclear and mitochondrial markers supported a monophyletic Aleyrodidae and the basal positioning of B. tabaci Uganda-1 to the sub-Saharan group of species. Reciprocal cross-mating data and the co-cladogenesis pattern of the primary obligate endosymbiont 'Candidatus Portiera aleyrodidarum' from 11 Bemisia genomes further supported the phylogenetic reconstruction to show that African cassava B. tabaci populations consist of just three biological species. We include comparative analyses of gene families related to detoxification, sugar metabolism, vector competency and evaluate the presence and function of horizontally transferred genes, essential for understanding the evolution and unique biology of constituent B. tabaci. s.l species.

CONCLUSIONS

These genomic resources have provided new and critical insights into the genetics underlying B. tabaci s.l. biology. They also provide a rich foundation for post-genomic research, including the selection of candidate gene-targets for innovative whitefly and virus-control strategies.

Screening of wild species and transcriptome profiling to identify differentially regulated genes in response to late blight resistance in potato

Late blight (Phytophthora infestans) is a serious disease of potatoes. The aim of this study was to screen wild potato species and identify differentially expressed genes (DEGs) associated with late blight resistance. Wild potato species such as PIN45 (Solanum pinnatisectum), CPH62 (Solanum cardiophyllum), JAM07 (Solanum jamesii), MCD24 (Solanum microdontum), PLD47 (Solanum polyadenium), and cv. Kufri Bahar (control) were tested by artificial inoculation of P. infestans under controlled conditions. Transcriptomes of the leaf tissues (96 h post-inoculation) were sequenced using the Illumina platform. Statistically significant (p < 0.05) DEGs were analyzed in wild species by comparison with the control, and upregulated (>2 log₂ fold change, FC) and downregulated (<-2 log₂ FC) genes were identified. DEGs were functionally characterized with Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Selected genes were validated by real-time PCR analysis to confirm RNA-seq results. We identified some upregulated genes associated with late blight resistance in wild species such as cytochrome P450, proline-rich protein, MYB transcription factor MYB139, ankyrin repeat-containing protein, and LRR receptor-like serine/threonine-protein kinase in PIN45; glucosyltransferase, fructose-bisphosphate aldolase, and phytophthora-inhibited protease 1 in CPH62; steroid binding protein and cysteine proteinase 3 in JAM07; glycine-rich cell wall structural protein 1 and RING finger protein in MCD24; and cysteine proteinase 3 and major latex protein in PLD47. On the other hand, downregulated genes in these species were snakin-2 and WRKY transcription factor 3 in PIN45; lichenase and phenylalanine ammonia-lyase 1 in CPH62; metallothionein and LRR receptor-like serine/threonine-protein kinase in JAM07; UDP-glucoronosyl/UDP-glucosyl transferase family protein and steroid binding protein in MCD24; and cytoplasmic small heat shock protein class I and phosphatase PLD47. Our study identified highly resistant wild potato species and underlying genes such as disease resistance, stress response, phytohormones, and transcription factors (e.g., MYB, WRKY, AP2/ERF, and AN1) associated with late blight resistance in wild potato species.

Molecular Mechanism Underlying the Sorghum sudanense (Piper) Stapf. Response to Osmotic Stress Determined via Single-Molecule Real-Time Sequencing and Next-Generation Sequencing

Drought, as a widespread environmental factor in nature, has become one of the most critical factors restricting the yield of forage grass. Sudangrass (Sorghum sudanense (Piper) Stapf.), as a tall and large grass, has a large biomass and is widely used as forage and biofuel. However, its growth and development are limited by drought stress. To obtain novel insight into the molecular mechanisms underlying the drought response and excavate drought tolerance genes in sudangrass, the first full-length transcriptome database of sudangrass under drought stress at different time points was constructed by combining single-molecule real-time sequencing (SMRT) and next-generation transcriptome sequencing (NGS). A total of 32.3 Gb of raw data was obtained, including 20,199 full-length transcripts with an average length of 1628 bp after assembly and correction. In total, 11,921 and 8559 up- and down-regulated differentially expressed genes were identified between the control group and plants subjected to drought stress. Additionally, 951 transcription factors belonging to 50 families and 358 alternative splicing events were found. A KEGG analysis of 158 core genes exhibiting continuous changes over time revealed that 'galactose metabolism' is a hub pathway and raffinose synthase 2 and β-fructofuranosidase are key genes in the response to drought stress. This study revealed the molecular mechanism underlying drought tolerance in sudangrass. Furthermore, the genes identified in this study provide valuable resources for further research into the response to drought stress.

Comparative transcriptome analysis to identify putative genes involved in carvacrol biosynthesis pathway in two species of Satureja, endemic medicinal herbs of Iran

Satureja is rich in phenolic monoterpenoids, mainly carvacrol, that is of interest due to diverse biological activities including antifungal and antibacterial. However, limited information is available regarding the molecular mechanisms underlying carvacrol biosynthesis and its regulation for this wonderful medicinal herb. To identify the putative genes involved in carvacrol and other monoterpene biosynthesis pathway, we generated a reference transcriptome in two endemic Satureja species of Iran, containing different yields (Satureja khuzistanica and Satureja rechingeri). Cross-species differential expression analysis was conducted between two species of Satureja. 210 and 186 transcripts related to terpenoid backbone biosynthesis were identified for S. khuzistanica and S. rechingeri, respectively. 29 differentially expressed genes (DEGs) involved in terpenoid biosynthesis were identified, and these DEGs were significantly enriched in monoterpenoid biosynthesis, diterpenoid biosynthesis, sesquiterpenoid and triterpenoid biosynthesis, carotenoid biosynthesis and ubiquinone and other terpenoid-quinone biosynthesis pathways. Expression patterns of S. khuzistanica and S. rechingeri transcripts involved in the terpenoid biosynthetic pathway were evaluated. In addition, we identified 19 differentially expressed transcription factors (such as MYC4, bHLH, and ARF18) that may control terpenoid biosynthesis. We confirmed the altered expression levels of DEGs that encode carvacrol biosynthetic enzymes using quantitative real-time PCR (qRT-PCR). This study is the first report on de novo assembly and transcriptome data analysis in Satureja which could be useful for an understanding of the main constituents of Satureja essential oil and future research in this genus.

Evolutionary history of an Alpine Archaeognath (Machilis pallida): Insights from different variant

Reconstruction of species histories is a central aspect of evolutionary biology. Patterns of genetic variation within and among populations can be leveraged to elucidate evolutionary processes and demographic histories. However, interpreting genetic signatures and unraveling the contributing processes can be challenging, in particular for non-model organisms with complex reproductive modes and genome organization. One way forward is the combined consideration of patterns revealed by different molecular markers (nuclear vs. mitochondrial) and types of variants (common vs. rare) that differ in their age, mode, and rate of evolution. Here, we applied this approach to RNAseq data generated for Machilis pallida (Archaeognatha), an Alpine jumping bristletail considered parthenogenetic and triploid. We generated de novo transcriptome and mitochondrial assemblies to obtain high-density data to investigate patterns of mitochondrial and common and rare nuclear variation in 17 M. pallida individuals sampled from all known populations. We find that the different variant types capture distinct aspects of the evolutionary history and discuss the observed patterns in the context of parthenogenesis, polyploidy, and survival during glaciation. This study highlights the potential of different variant types to gain insights into evolutionary scenarios even from challenging but often available data and the suitability of M. pallida and the genus Machilis as a study system for the evolution of sexual strategies and polyploidization during environmental change. We also emphasize the need for further research which will be stimulated and facilitated by these newly generated resources and insights.

Pest status, molecular evolution, and epigenetic factors derived from the genome assembly of Frankliniella fusca, a thysanopteran phytovirus vector

BACKGROUND

The tobacco thrips (Frankliniella fusca Hinds; family Thripidae; order Thysanoptera) is an important pest that can transmit viruses such as the tomato spotted wilt orthotospovirus to numerous economically important agricultural row crops and vegetables. The structural and functional genomics within the order Thysanoptera has only begun to be explored. Within the > 7000 known thysanopteran species, the melon thrips (Thrips palmi Karny) and the western flower thrips (Frankliniella occidentalis Pergrande) are the only two thysanopteran species with assembled genomes.

RESULTS

A genome of F. fusca was assembled by long-read sequencing of DNA from an inbred line. The final assembly size was 370 Mb with a single copy ortholog completeness of ~ 99% with respect to Insecta. The annotated genome of F. fusca was compared with the genome of its congener, F. occidentalis. Results revealed many instances of lineage-specific differences in gene content. Analyses of sequence divergence between the two Frankliniella species' genomes revealed substitution patterns consistent with positive selection in ~ 5% of the protein-coding genes with 1:1 orthologs. Further, gene content related to its pest status, such as xenobiotic detoxification and response to an ambisense-tripartite RNA virus (orthotospovirus) infection was compared with F. occidentalis. Several F. fusca genes related to virus infection possessed signatures of positive selection. Estimation of CpG depletion, a mutational consequence of DNA methylation, revealed that F. fusca genes that were downregulated and alternatively spliced in response to virus infection were preferentially targeted by DNA methylation. As in many other insects, DNA methylation was enriched in exons in Frankliniella, but gene copies with homology to DNA methyltransferase 3 were numerous and fragmented. This phenomenon seems to be relatively unique to thrips among other insect groups.

CONCLUSIONS

The F. fusca genome assembly provides an important resource for comparative genomic analyses of thysanopterans. This genomic foundation allows for insights into molecular evolution, gene regulation, and loci important to agricultural pest status.