Distribution and Genetic Diversity of Genes from Brazilian Bacillus thuringiensis Strains Toxic to Agricultural Insect Pests Revealed by Real-Time PCR

Bacillus thuringiensis is a Gram-positive aerobic bacterium and the most used biopesticide worldwide. Given the importance of B. thuringiensis strain characterization for the development of new bioinsecticides or transgenic events and the identification and classification of new B. thuringiensis genes and strains to understand its distribution and diversity, this work is aimed at creating a gene identification system based on qPCR reactions utilizing core B. thuringiensis genes cry1, cry2, cry3, cry4, cry5, app6, cry7, cry8, cry9, cry10, cry11, vpb1, vpa2, vip3, cyt1, and cyt2 for the characterization of 257 strains of B. thuringiensis. This system was based on the Invertebrate Bacteria Collection from Embrapa Genetic Resources and Biotechnology and analyzed (a) the degree of correlation between the distribution of these strains and the origin of the substrate from which the strain was isolated and (b) between its distribution and geoclimatic conditions. This study made it possible to observe that the cry1, cry2, and vip3A/B genes occur homogeneously in the Brazilian territory, and some genes are found in specific regions. The biggest reservoir of variability is within B. thuringiensis strains in each region, and it is suggested that both geoclimatic conditions and regional crops interfere with the genetic diversity of the B. thuringiensis strains present in the region, and B. thuringiensis strains can constantly exchange genetic information.

Biocontrol potential of bacteria belonging to the Bacillus subtilis group against pests and diseases of agricultural interest through genome exploration

The usage of microorganisms as biocontrol agents and biofertilizers has been recommended and recognized as an ecologically correct alternative to maintaining the productivity and safety of crops. Thus, the objectives of this work were to characterize twelve strains belonging to Invertebrate Bacteria Collection of Embrapa Genetic Resources and Biotechnology by molecular, morphological, and biochemical methods and to evaluate the pathogenicity of these strains against pests and diseases of agricultural interest. The morphological characteristic of the strains was performed according to the principles of Bergy's Manual of Systematic Bacteriology. The genomes of the 12 strains were sequenced in Macrogen, Inc. (Seoul, Korea) using the HiSeq2000 and GS-FLX Plus high-performance platforms. In the determination of antibiotic sensibility profiles, disc-diffusion methods (Cefar Diagnótica Ltda) were adopted©. Selective bioassays were carried out with insects of the Lepidoptera (Spodoptera frugiperda, Helicoverpa armigera, and Chrysodeixis includens), Coleoptera (Anthonomus grandis), Diptera (Aedes aegypti) and Hemiptera (Euschistus heros) orders, and with the nematode Caenorhabditis elegans. In addition, the antagonistic action of the phytopathogens Fusarium oxysporum f. sp. vasinfectum and Sclerotinia sclerotiorum against the strains under study, and in vitro assays of phosphate solubilization were also performed. Sequencing of the complete genome of the 12 strains determined that all of them belonged to the Bacillus subtilis sensu lato group. In the strains genome were detected genic clusters responsible for encoding secondary metabolites such as surfactin, iturin, fengycins/plipastatin, bacillomycin, bacillisin, and siderophores. Due to the production of these compounds, there was a survival reduction of the Lepidoptera order insects and a reduction in the phytopathogens mycelial growth. These results show that the species of group B. subtilis s.l. can become promising microbiological alternatives to pest and disease control.

Proteomic Insights of Cowpea Response to Combined Biotic and Abiotic Stresses

The co-occurrence of biotic and abiotic stresses in agricultural areas severely affects crop performance and productivity. Drought is one of the most adverse environmental stresses, and its association with root-knot nematodes further limits the development of several economically important crops, such as cowpea. Plant responses to combined stresses are complex and require novel adaptive mechanisms through the induction of specific biotic and abiotic signaling pathways. Therefore, the present work aimed to identify proteins involved in the resistance of cowpea to nematode and drought stresses individually and combined. We used the genotype CE 31, which is resistant to the root-knot nematode Meloidogyne spp. And tolerant to drought. Three biological replicates of roots and shoots were submitted to protein extraction, and the peptides were evaluated by LC-MS/MS. Shotgun proteomics revealed 2345 proteins, of which 1040 were differentially abundant. Proteins involved in essential biological processes, such as transcriptional regulation, cell signaling, oxidative processes, and photosynthesis, were identified. However, the main defense strategies in cowpea against cross-stress are focused on the regulation of hormonal signaling, the intense production of pathogenesis-related proteins, and the downregulation of photosynthetic activity. These are key processes that can culminate in the adaptation of cowpea challenged by multiple stresses. Furthermore, the candidate proteins identified in this study will strongly contribute to cowpea genetic improvement programs.

Characterization of microRNAs and Target Genes in Musa acuminata subsp. burmannicoides, var. Calcutta 4 during Interaction with Pseudocercospora musae

Endogenous microRNAs (miRNAs) are small non-coding RNAs that perform post-transcriptional regulatory roles across diverse cellular processes, including defence responses to biotic stresses. Pseudocercospora musae, the causal agent of Sigatoka leaf spot disease in banana (Musa spp.), is an important fungal pathogen of the plant. Illumina HiSeq 2500 sequencing of small RNA libraries derived from leaf material in Musa acuminata subsp. burmannicoides, var. Calcutta 4 (resistant) after inoculation with fungal conidiospores and equivalent non-inoculated controls revealed 202 conserved miRNAs from 30 miR-families together with 24 predicted novel miRNAs. Conserved members included those from families miRNA156, miRNA166, miRNA171, miRNA396, miRNA167, miRNA172, miRNA160, miRNA164, miRNA168, miRNA159, miRNA169, miRNA393, miRNA535, miRNA482, miRNA2118, and miRNA397, all known to be involved in plant immune responses. Gene ontology (GO) analysis of gene targets indicated molecular activity terms related to defence responses that included nucleotide binding, oxidoreductase activity, and protein kinase activity. Biological process terms associated with defence included response to hormone and response to oxidative stress. DNA binding and transcription factor activity also indicated the involvement of miRNA target genes in the regulation of gene expression during defence responses. sRNA-seq expression data for miRNAs and RNAseq data for target genes were validated using stem-loop quantitative real-time PCR (qRT-PCR). For the 11 conserved miRNAs selected based on family abundance and known involvement in plant defence responses, the data revealed a frequent negative correlation of expression between miRNAs and target host genes. This examination provides novel information on miRNA-mediated host defence responses, applicable in genetic engineering for the control of Sigatoka leaf spot disease.

Identification of cry genes in Bacillus thuringiensis by multiplex real-time PCR

Bacillus thuringiensis is an important bacterium of the group Bacillus cereus sensu lato due to its insecticidal properties. This microorganism has high genetic variability and its strains produce different Cry toxins, known as δ-endotoxins, which are mainly responsible for its toxic effect on insects that are agricultural pests or vector human diseases. Each strain can express a variety of cry genes, out of a total of 789 cry genes described so far. The detection of these genes is very important to characterize strains, as they may indicate their toxic potential. Several methods have been used to characterize B. thuringiensis strains, but one of the most common techniques is Polymerase Chain Reaction (PCR) from primers that detect the presence of cry genes. This technique has been optimized to make real-time multiplex quantitative PCR (qPCR) assays faster, more efficient, and safer, because the presence of three genes can be detected in a single reaction. In this work, a multiplex assay was developed to identify the presence of genes from the cry1A, cry1C, and cry1F families whose respective toxins are present in both bioinsecticides, and commercial transgenic plants used to control caterpillars. Specific primers were designed to identify the families of the cited genes and the system was validated with samples that were sequenced by next-generation sequencing (NGS). The system was implemented and used to characterize 214 strains. Of these, eight were submitted to conventional PCR, and the results matched, again validating the system. Thus, the application of the proposed technique allows the reliable evaluation through this system to detect the presence of the genes of the families cry1A, cry1C, and cry1F in samples of B. thuringiensis.

Integrated Omic Approaches Reveal Molecular Mechanisms of Tolerance during Soybean and Meloidogyne incognita Interactions

The root-knot nematode (RKN), Meloidogyne incognita, is a devastating soybean pathogen worldwide. The use of resistant cultivars is the most effective method to prevent economic losses caused by RKNs. To elucidate the mechanisms involved in resistance to RKN, we determined the proteome and transcriptome profiles from roots of susceptible (BRS133) and highly tolerant (PI 595099) Glycine max genotypes 4, 12, and 30 days after RKN infestation. After in silico analysis, we described major defense molecules and mechanisms considered constitutive responses to nematode infestation, such as mTOR, PI3K-Akt, relaxin, and thermogenesis. The integrated data allowed us to identify protein families and metabolic pathways exclusively regulated in tolerant soybean genotypes. Among them, we highlighted the phenylpropanoid pathway as an early, robust, and systemic defense process capable of controlling M. incognita reproduction. Associated with this metabolic pathway, 29 differentially expressed genes encoding 11 different enzymes were identified, mainly from the flavonoid and derivative pathways. Based on differential expression in transcriptomic and proteomic data, as well as in the expression profile by RT-qPCR, and previous studies, we selected and overexpressed the GmPR10 gene in transgenic tobacco to assess its protective effect against M. incognita. Transgenic plants of the T₂ generation showed up to 58% reduction in the M. incognita reproduction factor. Finally, data suggest that GmPR10 overexpression can be effective against the plant parasitic nematode M. incognita, but its mechanism of action remains unclear. These findings will help develop new engineered soybean genotypes with higher performance in response to RKN infections.

The early response of oil palm (Elaeis guineensis Jacq.) plants to water deprivation: Expression analysis of miRNAs and their putative target genes, and similarities with the response to salinity stress

Oil palm (Elaeis guineensis Jacq.) is a oilseed crop of great economic importance drastically affected by abiotic stresses. MicroRNAs (miRNAs) play crucial roles in transcription and post-transcription regulation of gene expression, being essential molecules in the response of plants to abiotic stress. To better understand the molecular mechanisms behind the response of young oil palm plants to drought stress, this study reports on the prediction and characterization of miRNAs and their putative target genes in the apical leaf of plants subjected to 14 days of water deprivation. Then, the data from this study were compared to the data from a similar study that focused on salinity stress. Both, the drought-and salt-responsive miRNAs and their putative target genes underwent correlation analysis to identify similarities and dissimilarities among them. Among the 81 identified miRNAs, 29 are specific for oil palm, including two (egu-miR28ds and egu-miR29ds) new ones - described for the first time. As for the expression profile, 62 miRNAs were significantly differentially expressed under drought stress, being five up-regulated (miR396e, miR159b, miR529b, egu-miR19sds, and egu-miR29ds) and 57 down-regulated. Transcription factors, such as MYBs, HOXs, and NF-Ys, were predicted as putative miRNA-target genes in oil palm under water deprivation; making them the most predominant group of such genes. Finally, the correlation analysis study revealed a group of putative target genes with similar behavior under salt and drought stresses. Those genes that are upregulated by these two abiotic stresses encode lncRNAs and proteins linked to stress tolerance, stress memory, modulation of ROS signaling, and defense response regulation to abiotic and biotic stresses. In summary, this study provides molecular evidence for the possible involvement of miRNAs in the drought stress response in oil palm. Besides, it shows that, at the molecular level, there are many similarities in the response of young oil palm plants to these two abiotic stresses.

The Identification of Small RNAs Differentially Expressed in Apple Buds Reveals a Potential Role of the Mir159-MYB Regulatory Module during Dormancy

Winter dormancy is an adaptative mechanism that temperate and boreal trees have developed to protect their meristems against low temperatures. In apple trees (Malus domestica), cold temperatures induce bud dormancy at the end of summer/beginning of the fall. Apple buds stay dormant during winter until they are exposed to a period of cold, after which they can resume growth (budbreak) and initiate flowering in response to warmer temperatures in spring. It is well-known that small RNAs modulate temperature responses in many plant species, but however, how small RNAs are involved in genetic networks of temperature-mediated dormancy control in fruit tree species remains unclear. Here, we have made use of a recently developed ARGONAUTE (AGO)-purification technique to isolate small RNAs from apple buds. A small RNA-seq experiment resulted in the identification of 17 micro RNAs (miRNAs) that change their pattern of expression in apple buds during dormancy. Furthermore, the functional analysis of their predicted target genes suggests a main role of the 17 miRNAs in phenylpropanoid biosynthesis, gene regulation, plant development and growth, and response to stimulus. Finally, we studied the conservation of the Arabidopsis thaliana regulatory miR159-MYB module in apple in the context of the plant hormone abscisic acid homeostasis.

Expression analysis of miRNAs and their putative target genes confirm a preponderant role of transcription factors in the early response of oil palm plants to salinity stress

BACKGROUND

Several mechanisms regulating gene expression contribute to restore and reestablish cellular homeostasis so that plants can adapt and survive in adverse situations. MicroRNAs (miRNAs) play roles important in the transcriptional and post-transcriptional regulation of gene expression, emerging as a regulatory molecule key in the responses to plant stress, such as cold, heat, drought, and salt. This work is a comprehensive and large-scale miRNA analysis performed to characterize the miRNA population present in oil palm (Elaeis guineensis Jacq.) exposed to a high level of salt stress, to identify miRNA-putative target genes in the oil palm genome, and to perform an in silico comparison of the expression profile of the miRNAs and their putative target genes.

RESULTS

A group of 79 miRNAs was found in oil palm, been 52 known miRNAs and 27 new ones. The known miRNAs found belonged to 28 families. Those miRNAs led to 229 distinct miRNA-putative target genes identified in the genome of oil palm. miRNAs and putative target genes differentially expressed under salinity stress were then selected for functional annotation analysis. The regulation of transcription, DNA-templated, and the oxidation-reduction process were the biological processes with the highest number of hits to the putative target genes, while protein binding and DNA binding were the molecular functions with the highest number of hits. Finally, the nucleus was the cellular component with the highest number of hits. The functional annotation of the putative target genes differentially expressed under salinity stress showed several ones coding for transcription factors which have already proven able to result in tolerance to salinity stress by overexpression or knockout in other plant species.

CONCLUSIONS

Our findings provide new insights into the early response of young oil palm plants to salinity stress and confirm an expected preponderant role of transcription factors - such as NF-YA3, HOX32, and GRF1 - in this response. Besides, it points out potential salt-responsive miRNAs and miRNA-putative target genes that one can utilize to develop oil palm plants tolerant to salinity stress.

Transcriptome Analysis and Knockdown of the Juvenile Hormone Esterase Gene Reveal Abnormal Feeding Behavior in the Sugarcane Giant Borer

The sugarcane giant borer (SGB), Telchin licus licus, is a pest that has strong economic relevance for sugarcane producers. Due to the endophytic behavior of the larva, current methods of management are inefficient. A promising biotechnological management option has been proposed based on RNA interference (RNAi), a process that uses molecules of double-stranded RNA (dsRNA) to specifically knock down essential genes and reduce insect survival. The selection of suitable target genes is often supported by omic sciences. Studies have shown that genes related to feeding adaptation processes are good candidates to be targeted by RNAi for pest management. Among those genes, esterases are highlighted because of their impact on insect development. In this study, the objective was to evaluate the transcriptome responses of the SGB’s gut in order to provide curated data of genes that could be used for pest management by RNAi in future studies. Further, we validated the function of an esterase-coding gene and its potential as a target for RNAi-based control. We sequenced the gut transcriptome of SGB larvae by Illumina HiSeq and evaluated its gene expression profiles in response to different diets (sugarcane stalk and artificial diet). We obtained differentially expressed genes (DEGs) involved in detoxification, digestion, and transport, which suggest a generalist mechanism of adaptation in SGB larvae. Among the DEGs, was identified and characterized a candidate juvenile hormone esterase gene (Tljhe). We knocked down the Tljhe gene by oral delivery of dsRNA molecules and evaluated gene expression in the gut. The survival and nutritional parameters of the larvae were measured along the developmental cycle of treated insects. We found that the gene Tljhe acts as a regulator of feeding behavior. The knockdown of Tljhe triggered a forced starvation state in late larval instars that significantly reduced the fitness of the larvae. However, the mechanism of action of this gene remains unclear, and the correlation between the expression of Tljhe and the levels of juvenile hormone (JH) metabolites in the hemolymph of the SGB must be assessed in future research.

Transcriptome Profiling-Based Analysis of Carbohydrate-Active Enzymes in Aspergillus terreus Involved in Plant Biomass Degradation

Given the global abundance of plant biomass residues, potential exists in biorefinery-based applications with lignocellulolytic fungi. Frequently isolated from agricultural cellulosic materials, Aspergillus terreus is a fungus efficient in secretion of commercial enzymes such as cellulases, xylanases and phytases. In the context of biomass saccharification, lignocellulolytic enzyme secretion was analyzed in a strain of A. terreus following liquid culture with sugarcane bagasse (SB) (1% w/v) and soybean hulls (SH) (1% w/v) as sole carbon source, in comparison to glucose (G) (1% w/v). Analysis of the fungal secretome revealed a maximum of 1.017 UI.mL^–1 xylanases after growth in minimal medium with SB, and 1.019 UI.mL^–1 after incubation with SH as carbon source. The fungal transcriptome was characterized on SB and SH, with gene expression examined in comparison to equivalent growth on G as carbon source. Over 8000 genes were identified, including numerous encoding enzymes and transcription factors involved in the degradation of the plant cell wall, with significant expression modulation according to carbon source. Eighty-nine carbohydrate-active enzyme (CAZyme)-encoding genes were identified following growth on SB, of which 77 were differentially expressed. These comprised 78% glycoside hydrolases, 8% carbohydrate esterases, 2.5% polysaccharide lyases, and 11.5% auxiliary activities. Analysis of the glycoside hydrolase family revealed significant up-regulation for genes encoding 25 different GH family proteins, with predominance for families GH3, 5, 7, 10, and 43. For SH, from a total of 91 CAZyme-encoding genes, 83 were also significantly up-regulated in comparison to G. These comprised 80% glycoside hydrolases, 7% carbohydrate esterases, 5% polysaccharide lyases, 7% auxiliary activities (AA), and 1% glycosyltransferases. Similarly, within the glycoside hydrolases, significant up-regulation was observed for genes encoding 26 different GH family proteins, with predominance again for families GH3, 5, 10, 31, and 43. A. terreus is a promising species for production of enzymes involved in the degradation of plant biomass. Given that this fungus is also able to produce thermophilic enzymes, this first global analysis of the transcriptome following cultivation on lignocellulosic carbon sources offers considerable potential for the application of candidate genes in biorefinery applications.

Comparative gut transcriptome analysis of Diatraea saccharalis in response to the dietary source

The sugarcane borer (Diatraea saccharalis, Fabricius, 1794) is a devastating pest that causes millions of dollars of losses each year to sugarcane producers by reducing sugar and ethanol yields. The control of this pest is difficult due to its endophytic behavior and rapid development. Pest management through biotechnological approaches has emerged in recent years as an alternative to currently applied methods. Genetic information about the target pests is often required to perform biotechnology-based management. The genomic and transcriptomic data for D. saccharalis are very limited. Herein, we report a tissue-specific transcriptome of D. saccharalis larvae and a differential expression analysis highlighting the physiological characteristics of this pest in response to two different diets: sugarcane and an artificial diet. Sequencing was performed on the Illumina HiSeq 2000 platform, and a de novo assembly was generated. A total of 27,626 protein-coding unigenes were identified, among which 1,934 sequences were differentially expressed between treatments. Processes such as defence, digestion, detoxification, signaling, and transport were highly represented among the differentially expressed genes (DEGs). Furthermore, seven aminopeptidase genes were identified as candidates to encode receptors of Cry proteins, which are toxins of Bacillus thuringiensis used to control lepidopteran pests. Since plant-insect interactions have produced a considerable number of adaptive responses in hosts and herbivorous insects, the success of phytophagous insects relies on their ability to overcome challenges such as the response to plant defences and the intake of nutrients. In this study, we identified metabolic pathways and specific genes involved in these processes. Thus, our data strongly contribute to the knowledge advancement of insect transcripts, which can be a source of target genes for pest management.

Evolutionarily conserved plant genes responsive to root-knot nematodes identified by comparative genomics

Root-knot nematodes (RKNs, genus Meloidogyne) affect a large number of crops causing severe yield losses worldwide, more specifically in tropical and sub-tropical regions. Several plant species display high resistance levels to Meloidogyne, but a general view of the plant immune molecular responses underlying resistance to RKNs is still lacking. Combining comparative genomics with differential gene expression analysis may allow the identification of widely conserved plant genes involved in RKN resistance. To identify genes that are evolutionary conserved across plant species, we used OrthoFinder to compared the predicted proteome of 22 plant species, including important crops, spanning 214 Myr of plant evolution. Overall, we identified 35,238 protein orthogroups, of which 6,132 were evolutionarily conserved and universal to all the 22 plant species (PLAnts Common Orthogroups-PLACO). To identify host genes responsive to RKN infection, we analyzed the RNA-seq transcriptome data from RKN-resistant genotypes of a peanut wild relative (Arachis stenosperma), coffee (Coffea arabica L.), soybean (Glycine max L.), and African rice (Oryza glaberrima Steud.) challenged by Meloidogyne spp. using EdgeR and DESeq tools, and we found 2,597 (O. glaberrima), 743 (C. arabica), 665 (A. stenosperma), and 653 (G. max) differentially expressed genes (DEGs) during the resistance response to the nematode. DEGs' classification into the previously characterized 35,238 protein orthogroups allowed identifying 17 orthogroups containing at least one DEG of each resistant Arachis, coffee, soybean, and rice genotype analyzed. Orthogroups contain 364 DEGs related to signaling, secondary metabolite production, cell wall-related functions, peptide transport, transcription regulation, and plant defense, thus revealing evolutionarily conserved RKN-responsive genes. Interestingly, the 17 DEGs-containing orthogroups (belonging to the PLACO) were also universal to the 22 plant species studied, suggesting that these core genes may be involved in ancestrally conserved immune responses triggered by RKN infection. The comparative genomic approach that we used here represents a promising predictive tool for the identification of other core plant defense-related genes of broad interest that are involved in different plant-pathogen interactions.

Transcriptome and gene expression analysis of three developmental stages of the coffee berry borer, Hypothenemus hampei

Coffee production is a global industry valued at approximately 173 billion US dollars. One of the main challenges facing coffee production is the management of the coffee berry borer (CBB), Hypothenemus hampei, which is considered the primary arthropod pest of coffee worldwide. Current control strategies are inefficient for CBB management. Although biotechnological alternatives, including RNA interference (RNAi), have been proposed in recent years to control insect pests, characterizing the genetics of the target pest is essential for the successful application of these emerging technologies. In this study, we employed RNA-seq to obtain the transcriptome of three developmental stages of the CBB (larva, female and male) to increase our understanding of the CBB life cycle in relation to molecular features. The CBB transcriptome was sequenced using Illumina Hiseq and assembled de novo. Differential gene expression analysis was performed across the developmental stages. The final assembly produced 29,434 unigenes, of which 4,664 transcripts were differentially expressed. Genes linked to crucial physiological functions, such as digestion and detoxification, were determined to be tightly regulated between the reproductive and nonreproductive stages of CBB. The data obtained in this study help to elucidate the critical roles that several genes play as regulatory elements in CBB development.

Spring Is Coming: Genetic Analyses of the Bud Break Date Locus Reveal Candidate Genes From the Cold Perception Pathway to Dormancy Release in Apple (Malus × domestica Borkh.)

Chilling requirement (CR) for bud dormancy completion determines the time of bud break in apple (Malus × domestica Borkh.). The molecular control of bud dormancy is highly heritable, suggesting a strong genetic control of the trait. An available Infinium II SNP platform for genotyping containing 8,788 single nucleotide polymorphic markers was employed, and linkage maps were constructed in a F₁ cross from the low CR M13/91 and the moderate CR cv. Fred Hough. These maps were used to identify quantitative trait loci (QTL) for bud break date as a trait related to dormancy release. A major QTL for bud break was detected at the beginning of linkage group 9 (LG9). This QTL remained stable during seven seasons in two different growing sites. To increase mapping efficiency in detecting contributing genes underlying this QTL, 182 additional SNP markers located at the locus for bud break were used. Combining linkage mapping and structural characterization of the region, the high proportion of the phenotypic variance in the trait explained by the QTL is related to the coincident positioning of Arabidopsis orthologs for ICE1, FLC, and PRE1 protein-coding genes. The proximity of these genes from the most explanatory markers of this QTL for bud break suggests potential genetic additive effects, reinforcing the hypothesis of inter-dependent mechanisms controlling dormancy induction and release in apple trees.

Analysis of the Transcriptome in Aspergillus tamarii During Enzymatic Degradation of Sugarcane Bagasse

The production of bioethanol from non-food agricultural residues represents an alternative energy source to fossil fuels for incorporation into the world's economy. Within the context of bioconversion of plant biomass into renewable energy using improved enzymatic cocktails, Illumina RNA-seq transcriptome profiling was conducted on a strain of Aspergillus tamarii, efficient in biomass polysaccharide degradation, in order to identify genes encoding proteins involved in plant biomass saccharification. Enzyme production and gene expression was compared following growth in liquid and semi-solid culture with steam-exploded sugarcane bagasse (SB) (1% w/v) and glucose (1% w/v) employed as contrasting sole carbon sources. Enzyme production following growth in liquid minimum medium supplemented with SB resulted in 0.626 and 0.711 UI.mL^-1 xylanases after 24 and 48 h incubation, respectively. Transcriptome profiling revealed expression of over 7120 genes, with groups of genes modulated according to solid or semi-solid culture, as well as according to carbon source. Gene ontology analysis of genes expressed following SB hydrolysis revealed enrichment in xyloglucan metabolic process and xylan, pectin and glucan catabolic process, indicating up-regulation of genes involved in xylanase secretion. According to carbohydrate-active enzyme (CAZy) classification, 209 CAZyme-encoding genes were identified with significant differential expression on liquid or semi-solid SB, in comparison to equivalent growth on glucose as carbon source. Up-regulated CAZyme-encoding genes related to cellulases (CelA, CelB, CelC, CelD) and hemicellulases (XynG1, XynG2, XynF1, XylA, AxeA, arabinofuranosidase) showed up to a 10-fold log2FoldChange in expression levels. Five genes from the AA9 (GH61) family, related to lytic polysaccharide monooxygenase (LPMO), were also identified with significant expression up-regulation. The transcription factor gene XlnR, involved in induction of hemicellulases, showed up-regulation on liquid and semi-solid SB culture. Similarly, the gene ClrA, responsible for regulation of cellulases, showed increased expression on liquid SB culture. Over 150 potential transporter genes were also identified with increased expression on liquid and semi-solid SB culture. This first comprehensive analysis of the transcriptome of A. tamarii contributes to our understanding of genes and regulatory systems involved in cellulose and hemicellulose degradation in this fungus, offering potential for application in improved enzymatic cocktail development for plant biomass degradation in biorefinery applications.

Anti-band 3 and anti-spectrin antibodies are increased in Plasmodium vivax infection and are associated with anemia

Clearance of non-infected red blood cells (nRBCs) is one of the main components of anemia associated with Plasmodium vivax malaria. Recently, we have shown that anemic patients with P. vivax infection had elevated levels of anti-RBCs antibodies, which could enhance in vitro phagocytosis of nRBCs and decrease their deformability. Using immunoproteomics, here we characterized erythrocytic antigens that are differentially recognized by autoantibodies from anemic and non-anemic patients with acute vivax malaria. Protein spots exclusively recognized by anemic P. vivax-infected patients were identified by mass spectrometry revealing band 3 and spectrin as the main targets. To confirm this finding, antibody responses against these specific proteins were assessed by ELISA. In addition, an inverse association between hemoglobin and anti-band 3 or anti-spectrin antibodies levels was found. Anemic patients had higher levels of IgG against both band 3 and spectrin than the non-anemic ones. To determine if these autoantibodies were elicited because of molecular mimicry, we used in silico analysis and identified P. vivax proteins that share homology with human RBC proteins such as spectrin, suggesting that infection drives autoimmune responses. These findings suggest that band 3 and spectrin are potential targets of autoantibodies that may be relevant for P. vivax malaria-associated anemia.

Systemic and sex-biased regulation of OBP expression under semiochemical stimuli

Constitutive expression of Odorant-Binding Proteins (OBPs) in antennae and other body parts has been examined mainly to infer their involvement in insect olfaction, while their regulation in response to semiochemical stimuli has remained poorly known. Previous studies of semiochemical response were basically done using electrophysiology, which integrates the response of the set of OBPs present in an antenna or sensillum, without revealing the regulation of OBPs or which ones might be involved. In this study we used boll weevil as a model and mined its OBPs by RNA-Seq to study their simultaneous antennal expression by qPCR under controlled semiochemical stimuli with aggregation pheromone and plant volatiles. In the absence of a semiochemical stimulus, 23 of 24 OBPs were constitutively expressed in the antenna in both sexes. Semiochemicals changed systemically the expression of OBPs in both sexes. There were different patterns of up- and down-regulation in female antennae for each semiochemical stimulus, consistent with female chemical ecology. On the other hand, the only response in males was down-regulation of some OBPs. We suggest that these systemic changes in OBP expression might be related to enhancing detection of the semiochemical stimuli and/or priming the olfactory system to detect other environmental chemicals.

Gene expression analysis in Musa acuminata during compatible interactions with Meloidogyne incognita

BACKGROUND AND AIMS

Endoparasitic root-knot nematodes (RKNs) ( Meloidogyne spp.) cause considerable losses in banana ( Musa spp.), with Meloidogyne incognita a predominant species in Cavendish sub-group bananas. This study investigates the root transcriptome in Musa acuminata genotypes 4297-06 (AA) and Cavendish Grande Naine (CAV; AAA) during early compatible interactions with M. incognita .

METHODS

Roots were analysed by brightfield light microscopy over a 35 d period to examine nematode penetration and morphological cell transformation. RNA samples were extracted 3, 7 and 10 days after inoculation (DAI) with nematode J2 juveniles, and cDNA libraries were sequenced using lllumina HiSeq technology. Sequences were mapped to the M. acuminata ssp. malaccensis var. Pahang genome sequence, differentially expressed genes (DEGs) identified and transcript representation determined by gene set enrichment and pathway mapping.

KEY RESULTS

Microscopic analysis revealed a life cycle of M. incognita completing in 24 d in CAV and 27 d in 4279-06. Comparable numbers of DEGs were up- and downregulated in each genotype, with potential involvement of many in early host defence responses involving reactive oxygen species and jasmonate/ethylene signalling. DEGs revealed concomitant auxin metabolism and cell wall modification processes likely to be involved in giant cell formation. Notable transcripts related to host defence included those coding for leucine-rich repeat receptor-like serine/threonine-protein kinases, peroxidases, thaumatin-like pathogenesis-related proteins, and DREB, ERF, MYB, NAC and WRKY transcription factors. Transcripts related to giant cell development included indole acetic acid-amido synthetase GH3.8 genes, involved in auxin metabolism, as well as genes encoding expansins and hydrolases, involved in cell wall modification.

CONCLUSIONS

Expression analysis in M. acuminata during compatible interactions with RKNs provides insights into genes modulated during infection and giant cell formation. Increased understanding of both defence responses to limit parasitism during compatible interactions and effector-targeted host genes in this complex interaction will facilitate the development of genetic improvement measures for RKNs.

Identification and expression profile of odorant-binding proteins in Halyomorpha halys (Hemiptera: Pentatomidae)

The brown marmorated stink bug, Halyomorpha halys, is a devastating invasive species in the USA. Similar to other insects, olfaction plays an important role in its survival and reproduction. As odorant-binding proteins (OBPs) are involved in the initial semiochemical recognition steps, we used RNA-Sequencing (RNA-Seq) to identify OBPs in its antennae, and studied their expression pattern in different body parts under semiochemical stimulation by either aggregation or alarm pheromone or food odorants. Thirty full-length putative HhalOBPs were identified, corresponding to 22 'classic' OBPs and eight 'Plus-C' OBPs. The similarity amongst them ranged from 4.95-70.92%, and with another 325 hemipteran OBPs similarity ranged from 1.94-91.51%, the highest levels being with other stink bug OBPs. Phylogenetic analysis confirmed the monophyly of seven groups of stink bug and other hemipteran OBPs. All 30 HhalOBPs were expressed and about 2/3 were expressed primarily in antennae. The expression of 21 HhalOBPs was higher in the antennae under alarm pheromone stimulus, indicating that multiple OBPs may be responding to this pheromone. Two were highest in antennae under aggregation pheromone stimulus. These findings should provide a basis for understanding the physiological functions of HhalOBPs and the chemosensory perception of this pest, which may help to uncover new control targets for behavioural interference.

Transcriptome-Based Identification of Highly Similar Odorant-Binding Proteins among Neotropical Stink Bugs and Their Egg Parasitoid

Olfaction plays a fundamental role in insect survival through resource location and intra and interspecific communications. We used RNA-Seq to analyze transcriptomes for odorant-binding proteins (OBPs) from major stink bug pest species in Brazil, Euschistus heros, Chinavia ubica, and Dichelops melacanthus, and from their egg parasitoid, Telenomus podisi. We identified 23 OBPs in E. heros, 25 OBPs in C. ubica, 9 OBPs in D. melacanthus, and 7 OBPs in T. podisi. The deduced amino acid sequences of the full-length OBPs had low intraspecific similarity, but very high similarity between two pairs of OBPs from E. heros and C. ubica (76.4 and 84.0%) and between two pairs of OBPs from the parasitoid and its preferred host E. heros (82.4 and 88.5%), confirmed by a high similarity of their predicted tertiary structures. The similar pairs of OBPs from E. heros and C. ubica may suggest that they have derived from a common ancestor, and retain the same biological function to bind a ligand perceived or produced in both species. The T. podisi OBPs similar to E. heros were not orthologous to any known hymenopteran OBPs, and may have evolved independently and converged to the host OBPs, providing a possible basis for the host location of T. podisi using E. heros semiochemical cues.

The genome sequence of Pseudoplusia includens single nucleopolyhedrovirus and an analysis of p26 gene evolution in the baculoviruses

BACKGROUND

Pseudoplusia includens single nucleopolyhedrovirus (PsinSNPV-IE) is a baculovirus recently identified in our laboratory, with high pathogenicity to the soybean looper, Chrysodeixis includens (Lepidoptera: Noctuidae) (Walker, 1858). In Brazil, the C. includens caterpillar is an emerging pest and has caused significant losses in soybean and cotton crops. The PsinSNPV genome was determined and the phylogeny of the p26 gene within the family Baculoviridae was investigated.

RESULTS

The complete genome of PsinSNPV was sequenced (Roche 454 GS FLX - Titanium platform), annotated and compared with other Alphabaculoviruses, displaying a genome apparently different from other baculoviruses so far sequenced. The circular double-stranded DNA genome is 139,132 bp in length, with a GC content of 39.3 % and contains 141 open reading frames (ORFs). PsinSNPV possesses the 37 conserved baculovirus core genes, 102 genes found in other baculoviruses and 2 unique ORFs. Two baculovirus repeat ORFs (bro) homologs, bro-a (Psin33) and bro-b (Psin69), were identified and compared with Chrysodeixis chalcites nucleopolyhedrovirus (ChchNPV) and Trichoplusia ni single nucleopolyhedrovirus (TnSNPV) bro genes and showed high similarity, suggesting that these genes may be derived from an ancestor common to these viruses. The homologous repeats (hrs) are absent from the PsinSNPV genome, which is also the case in ChchNPV and TnSNPV. Two p26 gene homologs (p26a and p26b) were found in the PsinSNPV genome. P26 is thought to be required for optimal virion occlusion in the occlusion bodies (OBs), but its function is not well characterized. The P26 phylogenetic tree suggests that this gene was obtained from three independent acquisition events within the Baculoviridae family. The presence of a signal peptide only in the PsinSNPV p26a/ORF-20 homolog indicates distinct function between the two P26 proteins.

CONCLUSIONS

PsinSNPV has a genomic sequence apparently different from other baculoviruses sequenced so far. The complete genome sequence of PsinSNPV will provide a valuable resource, contributing to studies on its molecular biology and functional genomics, and will promote the development of this virus as an effective bioinsecticide.

A directed approach for the identification of transcripts harbouring the spliced leader sequence and the effect of trans-splicing knockdown in Schistosoma mansoni

Schistosomiasis is a major neglected tropical disease caused by trematodes from the genus Schistosoma. Because schistosomes exhibit a complex life cycle and numerous mechanisms for regulating gene expression, it is believed that spliced leader (SL) trans-splicing could play an important role in the biology of these parasites. The purpose of this study was to investigate the function of trans-splicing in Schistosoma mansoni through analysis of genes that may be regulated by this mechanism and via silencing SL-containing transcripts through RNA interference. Here, we report our analysis of SL transcript-enriched cDNA libraries from different S. mansoni life stages. Our results show that the trans-splicing mechanism is apparently not associated with specific genes, subcellular localisations or life stages. In cross-species comparisons, even though the sets of genes that are subject to SL trans-splicing regulation appear to differ between organisms, several commonly shared orthologues were observed. Knockdown of trans-spliced transcripts in sporocysts resulted in a systemic reduction of the expression levels of all tested trans-spliced transcripts; however, the only phenotypic effect observed was diminished larval size. Further studies involving the findings from this work will provide new insights into the role of trans-splicing in the biology of S. mansoni and other organisms. All Expressed Sequence Tags generated in this study were submitted to dbEST as five different libraries. The accessions for each library and for the individual sequences are as follows: (i) adult worms of mixed sexes (LIBEST_027999: JZ139310 - JZ139779), (ii) female adult worms (LIBEST_028000: JZ139780 - JZ140379), (iii) male adult worms (LIBEST_028001: JZ140380 - JZ141002), (iv) eggs (LIBEST_028002: JZ141003 - JZ141497) and (v) schistosomula (LIBEST_028003: JZ141498 - JZ141974).

Proteomic analysis of Trypanosoma cruzi response to ionizing radiation stress

Trypanosoma cruzi, the causative agent of Chagas disease, is extremely resistant to ionizing radiation, enduring up to 1.5 kGy of gamma rays. Ionizing radiation can damage the DNA molecule both directly, resulting in double-strand breaks, and indirectly, as a consequence of reactive oxygen species production. After a dose of 500 Gy of gamma rays, the parasite genome is fragmented, but the chromosomal bands are restored within 48 hours. Under such conditions, cell growth arrests for up to 120 hours and the parasites resume normal growth after this period. To better understand the parasite response to ionizing radiation, we analyzed the proteome of irradiated (4, 24, and 96 hours after irradiation) and non-irradiated T. cruzi using two-dimensional differential gel electrophoresis followed by mass spectrometry for protein identification. A total of 543 spots were found to be differentially expressed, from which 215 were identified. These identified protein spots represent different isoforms of only 53 proteins. We observed a tendency for overexpression of proteins with molecular weights below predicted, indicating that these may be processed, yielding shorter polypeptides. The presence of shorter protein isoforms after irradiation suggests the occurrence of post-translational modifications and/or processing in response to gamma radiation stress. Our results also indicate that active translation is essential for the recovery of parasites from ionizing radiation damage. This study therefore reveals the peculiar response of T. cruzi to ionizing radiation, raising questions about how this organism can change its protein expression to survive such a harmful stress.

Nucleotide excision repair in Trypanosoma brucei: specialization of transcription-coupled repair due to multigenic transcription

Nucleotide excision repair (NER) is a highly conserved genome repair pathway acting on helix distorting DNA lesions. NER is divided into two subpathways: global genome NER (GG-NER), which is responsible for repair throughout genomes, and transcription-coupled NER (TC-NER), which acts on lesions that impede transcription. The extent of the Trypanosoma brucei genome that is transcribed is highly unusual, since most genes are organized in multigene transcription units, each transcribed from a single promoter. Given this transcription organization, we have addressed the importance of NER to T. brucei genome maintenance by performing RNAi against all predicted contributing factors. Our results indicate that TC-NER is the main pathway of NER repair, but only CSB, XPBz and XPG contribute. Moreover, we show that UV lesions are inefficiently repaired in T. brucei, perhaps due to preferential use of RNA polymerase translesion synthesis. RNAi of XPC and DDB was found to be lethal, and we show that these factors act in inter-strand cross-link repair. XPD and XPB appear only to act in transcription, not repair. This work indicates that the predominance of multigenic transcription in T. brucei has resulted in pronounced adaptation of NER relative to the host and may be an attractive drug target.

Trypanosoma cruzi gene expression in response to gamma radiation

Trypanosoma cruzi is an organism highly resistant to ionizing radiation. Following a dose of 500 Gy of gamma radiation, the fragmented genomic DNA is gradually reconstructed and the pattern of chromosomal bands is restored in less than 48 hours. Cell growth arrests after irradiation but, while DNA is completely fragmented, RNA maintains its integrity. In this work we compared the transcriptional profiles of irradiated and non-irradiated epimastigotes at different time points after irradiation using microarray. In total, 273 genes were differentially expressed; from these, 160 were up-regulated and 113 down-regulated. We found that genes with predicted functions are the most prevalent in the down-regulated gene category. Translation and protein metabolic processes, as well as generation of precursor of metabolites and energy pathways were affected. In contrast, the up-regulated category was mainly composed of obsolete sequences (which included some genes of the kinetoplast DNA), genes coding for hypothetical proteins, and Retrotransposon Hot Spot genes. Finally, the tyrosyl-DNA phosphodiesterase 1, a gene involved in double-strand DNA break repair process, was up-regulated. Our study demonstrated the peculiar response to ionizing radiation, raising questions about how this organism changes its gene expression to manage such a harmful stress.

Evidence for reductive genome evolution and lateral acquisition of virulence functions in two Corynebacterium pseudotuberculosis strains

Background

Corynebacterium pseudotuberculosis, a Gram-positive, facultative intracellular pathogen, is the etiologic agent of the disease known as caseous lymphadenitis (CL). CL mainly affects small ruminants, such as goats and sheep; it also causes infections in humans, though rarely. This species is distributed worldwide, but it has the most serious economic impact in Oceania, Africa and South America. Although C. pseudotuberculosis causes major health and productivity problems for livestock, little is known about the molecular basis of its pathogenicity.

Methodology and Findings

We characterized two C. pseudotuberculosis genomes (Cp1002, isolated from goats; and CpC231, isolated from sheep). Analysis of the predicted genomes showed high similarity in genomic architecture, gene content and genetic order. When C. pseudotuberculosis was compared with other Corynebacterium species, it became evident that this pathogenic species has lost numerous genes, resulting in one of the smallest genomes in the genus. Other differences that could be part of the adaptation to pathogenicity include a lower GC content, of about 52%, and a reduced gene repertoire. The C. pseudotuberculosis genome also includes seven putative pathogenicity islands, which contain several classical virulence factors, including genes for fimbrial subunits, adhesion factors, iron uptake and secreted toxins. Additionally, all of the virulence factors in the islands have characteristics that indicate horizontal transfer.

Conclusions

These particular genome characteristics of C. pseudotuberculosis, as well as its acquired virulence factors in pathogenicity islands, provide evidence of its lifestyle and of the pathogenicity pathways used by this pathogen in the infection process. All genomes cited in this study are available in the NCBI Genbank database (http://www.ncbi.nlm.nih.gov/genbank/) under accession numbers CP001809 and CP001829.

Ecology of malaria parasites infecting Southeast Asian macaques: evidence from cytochrome b sequences

Although malaria parasites infecting non-human primates are important models for human malaria, little is known of the ecology of infection by these parasites in the wild. We extensively sequenced cytochrome b (cytb) of malaria parasites (Apicomplexa: Haemosporida) from free-living southeast Asian monkeys Macaca nemestrina and Macaca fascicularis. The two most commonly observed taxa were Plasmodium inui and Hepatocystis sp., but certain other sequences did not cluster closely with any previously sequenced species. Most of the major clades of parasites were found in both Macaca species, and the two most commonly occurring parasite infected the two Macaca species at approximately equal levels. However, P. inui showed evidence of genetic differentiation between the populations infecting the two Macaca species, suggesting limited movement of this parasite among hosts. Moreover, coinfection with Plasmodium and Hepatocystis species occurred significantly less frequently than expected on the basis of the rates of infection with either taxon alone, suggesting the possibility of competitive exclusion. The results revealed unexpectedly complex communities of Plasmodium and Hepatocystis taxa infecting wild southeast Asian monkeys. Parasite taxa differed with respect to both the frequency of between-host movement and their frequency of coinfection.

Nucleotide sequence polymorphism at the apical membrane antigen-1 locus reveals population history of Plasmodium vivax in Thailand

Apical membrane antigen-1 is a candidate for inclusion in a vaccine for the human malaria parasite Plasmodium vivax. We collected 231 complete sequences of the gene encoding this antigen (pvama-1) from three regions of Thailand, the most extensive collection to date of sequences at this locus. The domain II loop (previously mentioned as a potential vaccine component) was almost completely conserved, with a single amino acid variant (I313R) observed in a single sequence. The 3' portion of the gene (domain II through the stop codon) showed significantly lower nucleotide diversity than the 5' portion (start codon through domain I); and a given domain I sequence might be found in a haplotype with more than one domain II sequence. These results imply a hotspot of recombination between domains I and II. We found significant geographic subdivision among the three regions of Thailand (NW, East, and South) in which collections were made in 2007. Numbers of P. vivax infections have experienced overall declines since 1990 in all three regions; but the decline has been most recent in the NW, and there has been a rebound in numbers of infections in the South since 2000. Consistent with population history, amino acid sequence diversity was greatest in the NW. The South, which had by far the lowest sequence diversity of the three regions, showed signs of a population that has expanded from a small number of founders after a bottleneck.

Polymorphism at the apical membrane antigen 1 locus reflects the world population history of Plasmodium vivax

Background

In malaria parasites (genus Plasmodium), ama-1 is a highly polymorphic locus encoding the Apical Membrane Protein-1, and there is evidence that the polymorphism at this locus is selectively maintained. We tested the hypothesis that polymorphism at the ama-1 locus reflects population history in Plasmodium vivax, which is believed to have originated in Southeast Asia and is widely geographically distributed. In particular, we tested for a signature of the introduction of P. vivax into the New World at the time of the European conquest and African slave trade and subsequent population expansion.

Results

One hundred and five ama-1 sequences were generated and analyzed from samples from six different Brazilian states and compared with database sequences from the Old World. Old World populations of P. vivax showed substantial evidence of population substructure, with high sequence divergence among localities at both synonymous and nonsynonymous sites, while Brazilian isolates showed reduced diversity and little population substructure.

Conclusion

These results show that genetic diversity in P. vivax AMA-1 reflects population history, with population substructure characterizing long-established Old World populations, whereas Brazilian populations show evidence of loss of diversity and recent population expansion.

Note

Nucleotide sequence data reported is this paper are available in the GenBank™ database under the accession numbers EF031154 – EF031216 and EF057446 – EF057487

Association between particular polymorphic residues on apical membrane antigen 1 (AMA-1) and platelet levels in patients with vivax malaria

Apical membrane antigen 1 (AMA-1) is an immunogenic type 1 integral membrane protein, present in all Plasmodium spp., that probably has a role in the initiation of the invasion process of the erythrocyte. The DNA sequence of variable domain I of the Plasmodium vivax ama1 gene was sequenced in Brazilian isolates obtained from thrombocytopenic patients (n = 32) and patients with normal platelet counts (n = 22). There was a significant negative correlation between parasite density and platelet counts. It was concluded that there is an additional effect of sequence on platelet counts. The presence of amino-acid residues Y(193) and S(210) was associated significantly with normal platelet counts in P. vivax malaria, independent of the level of parasitaemia (p <0.0001). These data have implications for AMA-1-based vaccine design and suggest the possible use of this molecule as a marker of morbidity.