Transcriptomic investigation of the molecular mechanisms underlying resistance to the neonicotinoid thiamethoxam and the pyrethroid lambda-cyhalothrin in Euschistus heros (Hemiptera: Pentatomidae).

BACKGROUND: Laboratory-selected resistant strains of Euschistus heros to thiamethoxam (NEO) and lambda-cyhalothrin (PYR) were recently reported in Brazil. However, the mechanisms conferring resistance to these insecticides in E. heros remain unresolved. We utilized comparative transcriptome profiling and single nucleotide polymorphism (SNP) calling of susceptible and resistant strains of E. heros to investigate the molecular mechanism(s) underlying resistance. RESULTS: The E. heros transcriptome was assembled, generating 91 673 transcripts with a mean length of 720 bp and N50 of 1795 bp. Comparative gene expression analysis between the susceptible (SUS) and NEO strains identified 215 significantly differentially expressed (DE) transcripts. DE transcripts associated with the xenobiotic metabolism were all up-regulated in the NEO strain. The comparative analysis of the SUS and PYR strains identified 204 DE transcripts, including an esterase (esterase FE4), a glutathione-S-transferase, an ABC transporter (ABCC1) and aquaporins that were up-regulated in the PYR strain. We identified 9588 and 15 043 nonsynonymous SNPs in the PYR and NEO strains. One of the SNPs (D70N) detected in the NEO strain occurs in a subunit (α5) of the nAChRs, the target site of neonicotinoid insecticides. Nevertheless, this residue position in α5 is not conserved among insects. CONCLUSIONS: Neonicotinoid and pyrethroid resistance in laboratory-selected E. heros is associated with a potential metabolic resistance mechanism by the overexpression of proteins commonly involved in the three phases of xenobiotic metabolism. Together these findings provide insight into the potential basis of resistance in E. heros and will inform the development and implementation of resistance management strategies against this important pest. © 2023 Society of Chemical Industry.

Cross-species transcriptomes reveal species-specific and shared molecular adaptations for plants development on iron-rich rocky outcrops soils.

BACKGROUND: Canga is the Brazilian term for the savanna-like vegetation harboring several endemic species on iron-rich rocky outcrops, usually considered for mining activities. Parkia platycephala Benth. and Stryphnodendron pulcherrimum (Willd.) Hochr. naturally occur in the cangas of Serra dos Carajás (eastern Amazonia, Brazil) and the surrounding forest, indicating high phenotypic plasticity. The morphological and physiological mechanisms of the plants' establishment in the canga environment are well studied, but the molecular adaptative responses are still unknown. To understand these adaptative responses, we aimed to identify molecular mechanisms that allow the establishment of these plants in the canga environment. RESULTS: Plants were grown in canga and forest substrates collected in the Carajás Mineral Province. RNA was extracted from pooled leaf tissue, and RNA-seq paired-end reads were assembled into representative transcriptomes for P. platycephala and S. pulcherrimum containing 31,728 and 31,311 primary transcripts, respectively. We identified both species-specific and core molecular responses in plants grown in the canga substrate using differential expression analyses. In the species-specific analysis, we identified 1,112 and 838 differentially expressed genes for P. platycephala and S. pulcherrimum, respectively. Enrichment analyses showed that unique biological processes and metabolic pathways were affected for each species. Comparative differential expression analysis was based on shared single-copy orthologs. The overall pattern of ortholog expression was species-specific. Even so, we identified almost 300 altered genes between plants in canga and forest substrates with conserved responses in the two species. The genes were functionally associated with the response to light stimulus and the circadian rhythm pathway. CONCLUSIONS: Plants possess species-specific adaptative responses to cope with the substrates. Our results also suggest that plants adapted to both canga and forest environments can adjust the circadian rhythm in a substrate-dependent manner. The circadian clock gene modulation might be a central mechanism regulating the plants' development in the canga substrate in the studied legume species. The mechanism may be shared as a common mechanism to abiotic stress compensation in other native species.

Main Molecular Pathways Associated with Copper Tolerance Response in Imperata cylindrica by de novo Transcriptome Assembly.

The metallophyte Imperata cylindrica inhabits copper (Cu) polluted soils in large areas from Central Chile. Here, we subjected clonal vegetative plantlets to 300 mg Cu kg-1 of substrate for 21 days to identify the main molecular pathways involved in the response to Cu stress. Transcriptomic analyses were performed for shoots and roots, with and without Cu supply. RNA-Seq and de novo transcriptome assembly were performed to identify the gene response associated with molecular mechanisms of Cu tolerance in I. cylindrica. De novo transcriptome revealed a total of 200,521 transcripts (1777 bp) comprising ~91% complete ultra-conserved genes in the eukaryote and Plantae database. The differentially expressed genes (DEGs) in roots were 7386, with 3558 of them being up-regulated and the other 3828 down-regulated. The transcriptome response in shoots was significantly less, showing only 13 up-regulated and 23 down-regulated genes. Interestingly, DEGs mainly related with actin and cytoskeleton formation, and to a minor degree, some DEGs associated with metal transporters and superoxide dismutase activity in root tissues were found. These transcriptomic results suggest that cytoskeleton could be acting as a mechanism of Cu-binding in the root, resulting in a high Cu tolerance response in this metallophyte, which deserve to be analyzed ultra-structurally. Our study contributes to reinforcing the potential of I. cylindrica as a candidate plant species to be used as a phytoremediation agent in Cu-contaminated environments.

Transcriptome of the Southern Muriqui Brachyteles arachnoides (Primates:Platyrrhini), a Critically Endangered New World Monkey: Evidence of Adaptive Evolution.

The southern muriqui (Brachyteles arachnoides) is the largest neotropical primate. This species is endemic to Brazil and is currently critically endangered due to its habitat destruction. The genetic basis underlying adaptive traits of New World monkeys has been a subject of interest to several investigators, with significant concern about genes related to the immune system. In the absence of a reference genome, RNA-seq and de novo transcriptome assembly have proved to be valuable genetic procedures for accessing gene sequences and testing evolutionary hypotheses. We present here a first report on the sequencing, assembly, annotation and adaptive selection analysis for thousands of transcripts of B. arachnoides from two different samples, corresponding to 13 different blood cells and fibroblasts. We assembled 284,283 transcripts with N50 of 2,940 bp, with a high rate of complete transcripts, with a median high scoring pair coverage of 88.2%, including low expressed transcripts, accounting for 72.3% of complete BUSCOs. We could predict and extract 81,400 coding sequences with 79.8% of significant BLAST hit against the Euarchontoglires SwissProt dataset. Of these 64,929 sequences, 34,084 were considered homologous to Supraprimate proteins, and of the remaining sequences (30,845), 94% were associated with a protein domain or a KEGG Orthology group, indicating potentially novel or specific protein-coding genes of B. arachnoides. We use the predicted protein sequences to perform a comparative analysis with 10 other primates. This analysis revealed, for the first time in an Atelid species, an expansion of APOBEC3G, extending this knowledge to all NWM families. Using a branch-site model, we searched for evidence of positive selection in 4,533 orthologous sets. This evolutionary analysis revealed 132 amino acid sites in 30 genes potentially evolving under positive selection, shedding light on primate genome evolution. These genes belonged to a wide variety of categories, including those encoding the innate immune system proteins (APOBEC3G, OAS2, and CEACAM1) among others related to the immune response. This work generated a set of thousands of complete sequences that can be used in other studies on molecular evolution and may help to unveil the evolution of primate genes. Still, further functional studies are required to provide an understanding of the underlying evolutionary forces modeling the primate genome.

A de novo transcriptome analysis revealed that photomorphogenic genes are required for carotenoid synthesis in the dark-grown carrot taproot.

Carotenoids are terpenoid pigments synthesized by all photosynthetic and some non-photosynthetic organisms. In plants, these lipophilic compounds are involved in photosynthesis, photoprotection, and phytohormone synthesis. In plants, carotenoid biosynthesis is induced by several environmental factors such as light including photoreceptors, such as phytochromes (PHYs) and negatively regulated by phytochrome interacting factors (PIFs). Daucus carota (carrot) is one of the few plant species that synthesize and accumulate carotenoids in the storage root that grows in darkness. Contrary to other plants, light inhibits secondary root growth and carotenoid accumulation suggesting the existence of new mechanisms repressed by light that regulate both processes. To identify genes induced by dark and repressed by light that regulate carotenoid synthesis and carrot root development, in this work an RNA-Seq analysis was performed from dark- and light-grown carrot roots. Using this high-throughput sequencing methodology, a de novo transcriptome model with 63,164 contigs was obtained, from which 18,488 were differentially expressed (DEG) between the two experimental conditions. Interestingly, light-regulated genes are preferably expressed in dark-grown roots. Enrichment analysis of GO terms with DEGs genes, validation of the transcriptome model and DEG analysis through qPCR allow us to hypothesize that genes involved in photomorphogenesis and light perception such as PHYA, PHYB, PIF3, PAR1, CRY2, FYH3, FAR1 and COP1 participate in the synthesis of carotenoids and carrot storage root development.

Revealing the Venomous Secrets of the Spider's Web.

Orb-weaving spiders use a highly strong, sticky and elastic web to catch their prey. These web properties alone would be enough for the entrapment of prey; however, these spiders may be hiding venomous secrets in the web, which current research is revealing. Here, we provide strong proteotranscriptomic evidence for the presence of toxin/neurotoxin-like proteins, defensins, and proteolytic enzymes on the web silk from Nephila clavipes spider. The results from quantitative-based transcriptomic and proteomic approaches showed that silk-producing glands produce an extensive repertoire of toxin/neurotoxin-like proteins, similar to those already reported in spider venoms. Meanwhile, the insect toxicity results demonstrated that these toxic components can be lethal and/or paralytic chemical weapons used for prey capture on the web, and the presence of fatty acids in the web may be a responsible mechanism opening the way to the web toxins for accessing the interior of prey's body, as shown here. Comparative phylogenomic-level evolutionary analyses revealed orthologous genes among two spider groups, Araneomorphae and Mygalomorphae, and the findings showed protein sequences similar to toxins found in the taxa Scorpiones and Hymenoptera in addition to Araneae. Overall, these data represent a valuable resource to further investigate other spider web toxin systems and also suggest that N. clavipes web is not a passive mechanical trap for prey capture, but it exerts an active role in prey paralysis/killing using a series of neurotoxins.

Transcript and metabolic adjustments triggered by drought in Ilex paraguariensis leaves.

MAIN CONCLUSION: Abscisic acid is involved in the drought response of Ilex paraguariensis. Acclimation includes root growth stimulation, stomatal closure, osmotic adjustment, photoprotection, and regulation of nonstructural carbohydrates and amino acid metabolisms. Ilex paraguariensis (yerba mate) is cultivated in the subtropical region of South America, where the occurrence of drought episodes limit yield. To explore the mechanisms that allow I. paraguariensis to overcome dehydration, we investigated (1) how gene expression varied between water-stressed and non-stressed plants and (2) in what way the modulation of gene expression was linked to physiological status and metabolite composition. A total of 4920 differentially expressed transcripts were obtained through RNA-Seq after water deprivation. Drought induced the expression of several transcripts involved in the ABA-signalling pathway. Stomatal closure and leaf osmotic adjustments were promoted to minimize water loss, and these responses were accompanied by a high transcriptional remodeling of stress perception, signalling and transcriptional regulation, the photoprotective and antioxidant systems, and other stress-responsive genes. Simultaneously, significant changes in metabolite contents were detected. Glutamine, phenylalanine, isomaltose, fucose, and malate levels were shown to be positively correlated with dehydration. Principal component analysis showed differences in the metabolic profiles of control and stressed leaves. These results provide a comprehensive overview of how I. paraguariensis responds to dehydration at transcriptional and metabolomic levels and provide further characterization of the molecular mechanisms associated with drought response in perennial subtropical species.

De novo assembly of Persea americana cv. 'Hass' transcriptome during fruit development.

BACKGROUND: Avocado (Persea americana Mill.) is a basal angiosperm from the Lauraceae family. This species has a diploid genome with an approximated size of ~ 920 Mbp and produces a climacteric, fleshy and oily fruit. The flowering and fruit set are particularly prolonged processes, lasting between one to three months, generating important differences in physiological ages of the fruit within the same tree. So far there is no detailed genomic information regarding this species, being the cultivar 'Hass' especially important for avocado growers worldwide. With the aim to explore the fruit avocado transcriptome and to identify candidate biomarkers to monitore fruit development, we carried out an RNA-Seq approach during 4 stages of 'Hass' fruit development: 150 days after fruit set (DAFS), 240 DAFS, 300 DAFS (harvest) and 390 DAFS (late-harvest). RESULTS: The 'Hass' de novo transcriptome contains 62,203 contigs (x̅=988 bp, N50 = 1050 bp). We found approximately an 85 and 99% of complete ultra-conserved genes in eukaryote and plantae database using BUSCO (Benchmarking Universal Single-Copy Orthologs) and CEGMA (Core Eukaryotic Gene Mapping Approach), respectively. Annotation was performed with BLASTx, resulting in a 58% of annotated contigs (90% of differentially expressed genes were annotated). Differentially expressed genes analysis (DEG; with False Discovery Rate ≤ 0.01) found 8672 genes considering all developmental stages. From this analysis, genes were clustered according to their expression pattern and 1209 genes show correlation with the four developmental stages. CONCLUSIONS: Candidate genes are proposed as possible biomarkers for monitoring the development of the 'Hass' avocado fruit associated with lipid metabolism, ethylene signaling pathway, auxin signaling pathway, and components of the cell wall.

Leaf transcriptome of two highly divergent genotypes of Urochloa humidicola (Poaceae), a tropical polyploid forage grass adapted to acidic soils and temporary flooding areas.

BACKGROUND: Urochloa humidicola (Koronivia grass) is a polyploid (6x to 9x) species that is used as forage in the tropics. Facultative apospory apomixis is present in most of the genotypes of this species, although one individual has been described as sexual. Molecular studies have been restricted to molecular marker approaches for genetic diversity estimations and linkage map construction. The objectives of the present study were to describe and compare the leaf transcriptome of two important genotypes that are highly divergent in terms of their phenotypes and reproduction modes: the sexual BH031 and the aposporous apomictic cultivar BRS Tupi. RESULTS: We sequenced the leaf transcriptome of Koronivia grass using an Illumina GAIIx system, which produced 13.09 Gb of data that consisted of 163,575,526 paired-end reads between the two libraries. We de novo-assembled 76,196 transcripts with an average length of 1,152 bp and filtered 35,093 non-redundant unigenes. A similarity search against the non-redundant National Center of Biotechnology Information (NCBI) protein database returned 65 % hits. We annotated 24,133 unigenes in the Phytozome database and 14,082 unigenes in the UniProtKB/Swiss-Prot database, assigned 108,334 gene ontology terms to 17,255 unigenes and identified 5,324 unigenes in 327 known metabolic pathways. Comparisons with other grasses via a reciprocal BLAST search revealed a larger number of orthologous genes for the Panicum species. The unigenes were involved in C4 photosynthesis, lignocellulose biosynthesis and flooding stress responses. A search for functional molecular markers revealed 4,489 microsatellites and 560,298 single nucleotide polymorphisms (SNPs). A quantitative real-time PCR analysis validated the RNA-seq expression analysis and allowed for the identification of transcriptomic differences between the two evaluated genotypes. Moreover, 192 unannotated sequences were classified as containing complete open reading frames, suggesting that the new, potentially exclusive genes should be further investigated. CONCLUSION: The present study represents the first whole-transcriptome sequencing of U. humidicola leaves, providing an important public information source of transcripts and functional molecular markers. The qPCR analysis indicated that the expression of certain transcripts confirmed the differential expression observed in silico, which demonstrated that RNA-seq is useful for identifying differentially expressed and unique genes. These results corroborate the findings from previous studies and suggest a hybrid origin for BH031.

De novo Assembly and Analysis of the Chilean Pencil Catfish Trichomycterus areolatus Transcriptome.

Trichomycterus areolatus is an endemic species of pencil catfish that inhabits the riffles and rapids of many freshwater ecosystems of Chile. Despite its unique adaptation to Chile's high gradient watersheds and therefore potential application in the investigation of ecosystem integrity and environmental contamination, relatively little is known regarding the molecular biology of this environmental sentinel. Here, we detail the assembly of the Trichomycterus areolatus transcriptome, a molecular resource for the study of this organism and its molecular response to the environment. RNA-Seq reads were obtained by next-generation sequencing with an Illumina® platform and processed using PRINSEQ. The transcriptome assembly was performed using TRINITY assembler. Transcriptome validation was performed by functional characterization with KOG, KEGG, and GO analyses. Additionally, differential expression analysis highlights sex-specific expression patterns, and a list of endocrine and oxidative stress related transcripts are included.

De novo transcriptome assembly of 'Angeleno' and 'Lamoon' Japanese plum cultivars (Prunus salicina).

Japanese plum (Prunus salicina L.) is a fruit tree of the Rosaceae family, which is an economically important stone fruit around the world. Currently, Japanese plum breeding programs combine traditional breeding and plant physiology strategies with genetic and genomic analysis. In order to understand the flavonoid pathway regulation and to develop molecular markers associated to the fuit skin color (EST-SSRs), we performed a next generation sequencing based on Illumina Hiseq2000 platform. A total of 22.4 GB and 21 GB raw data were obtained from 'Lamoon' and 'Angeleno' respectively, corresponding to 85,404,726 raw reads to 'Lamoon' and 79,781,666 to 'Angeleno'. A total of 139,775,975 reads were filtered after removing low-quality reads and trimming the adapter sequences. De novo transcriptome assembly was performed using CLC Genome Workbench software and a total of 54,584 unique contigs were generated, with an N50 of 1343 base pair (bp) and a mean length of 829 bp. This work contributed with a specific Japanese plum skin transcriptome, providing two libraries of contrasting fruit skin color phenotype (yellow and red) and increasing substantially the GB of raw data available until now for this specie.