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.

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.

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.

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.

Transcriptome Profiling of the Resistance Response of Musa acuminata subsp. burmannicoides, var. Calcutta 4 to Pseudocercospora musae.

Banana (Musa spp.), which is one of the world's most popular and most traded fruits, is highly susceptible to pests and diseases. Pseudocercospora musae, responsible for Sigatoka leaf spot disease, is a principal fungal pathogen of Musa spp., resulting in serious economic damage to cultivars in the Cavendish subgroup. The aim of this study was to characterize genetic components of the early immune response to P. musae in Musa acuminata subsp. burmannicoides, var. Calcutta 4, a resistant wild diploid. Leaf RNA samples were extracted from Calcutta 4 three days after inoculation with fungal conidiospores, with paired-end sequencing conducted in inoculated and non-inoculated controls using lllumina HiSeq 4000 technology. Following mapping to the reference M. acuminata ssp. malaccensis var. Pahang genome, differentially expressed genes (DEGs) were identified and expression representation analyzed on the basis of gene ontology enrichment, Kyoto Encyclopedia of Genes and Genomes orthology and MapMan pathway analysis. Sequence data mapped to 29,757 gene transcript models in the reference Musa genome. A total of 1073 DEGs were identified in pathogen-inoculated cDNA libraries, in comparison to non-inoculated controls, with 32% overexpressed. GO enrichment analysis revealed common assignment to terms that included chitin binding, chitinase activity, pattern binding, oxidoreductase activity and transcription factor (TF) activity. Allocation to KEGG pathways revealed DEGs associated with environmental information processing, signaling, biosynthesis of secondary metabolites, and metabolism of terpenoids and polyketides. With 144 up-regulated DEGs potentially involved in biotic stress response pathways, including genes involved in cell wall reinforcement, PTI responses, TF regulation, phytohormone signaling and secondary metabolism, data demonstrated diverse early-stage defense responses to P. musae. With increased understanding of the defense responses occurring during the incompatible interaction in resistant Calcutta 4, these data are appropriate for the development of effective disease management approaches based on genetic improvement through introgression of candidate genes in superior cultivars.

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 T2 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.

NBS-LRR-WRKY genes and protease inhibitors (PIs) seem essential for cowpea resistance to root-knot nematode.

Cowpea (Vigna unguiculata L. Walp) is a legume of great economic importance, however it is highly affected by nematodes. The present work aimed to identify proteins and genes involved in nematode resistance by proteomic and transcriptomic analysis. Plants of a genotype resistant (CE31) to root-knot nematode (Meloidogyne spp.) were collected 12 days after inoculation with Meloidogyne incognita and the total proteins and RNA were extracted from the root samples. Shotgun proteomic analysis was performed using an Orbitrap Elite mass spectrometer and the construction and sequencing of cDNA libraries were carried out in a Hi-Seq 2000 sequencing system. The proteomic and transcriptomic analyses revealed key processes involved in cowpea defense and some interesting candidates were further analyzed by RT-qPCR. Proteins and genes involved in essential biological processes were differentially accumulated such as, regulation of transcription, cell wall stiffening and microtubule-based process. However, the main defense strategies of Vigna unguiculata seem to be focused on the interaction of NBS-LRR and WRKY genes for the activation of R genes, production of protease inhibitors and maintenance of actin cytoskeleton. These are key processes that can culminate in the suppression of giant cell formation and consequently in the development of Meloidogyne incognita. SIGNIFICANCE: In this study, we identified proteins and transcripts regulated in cowpea resistant to the nematode Meloidogyne spp. upon inoculation. The results revealed key candidate genes involved in the activation of R genes, the production of protease inhibitors and maintenance of the actin cytoskeleton. These processes might be essential for cowpea resistance, as they can impede nematode nutrition, giant cell formation and consequently the development of Meloidogyne incognita.

Low levels of sulfur and cobalt during the pre- and periconceptional periods affect the oocyte yield of donors and the DNA methylome of preimplantation bovine embryos.

Maternal nutrition is critical in mammalian development, influencing the epigenetic reprogramming of gametes, embryos, and fetal programming. We evaluated the effects of different levels of sulfur (S) and cobalt (Co) in the maternal diet throughout the pre- and periconceptional periods on the biochemical and reproductive parameters of the donors and the DNA methylome of the progeny in Bos indicus cattle. The low-S/Co group differed from the control with respect to homocysteine, folic acid, B12, insulin growth factor 1, and glucose. The oocyte yield was lower in heifers from the low S/Co group than that in the control heifers. Embryos from the low-S/Co group exhibited 2320 differentially methylated regions (DMRs) across the genome compared with the control embryos. We also characterized candidate DMRs linked to the DNMT1 and DNMT3B genes in the blood and sperm cells of the adult progeny. A DMR located in DNMT1 that was identified in embryos remained differentially methylated in the sperm of the progeny from the low-S/Co group. Therefore, we associated changes in specific compounds in the maternal diet with DNA methylation modifications in the progeny. Our results help to elucidate the impact of maternal nutrition on epigenetic reprogramming in livestock, opening new avenues of research to study the effect of disturbed epigenetic patterns in early life on health and fertility in adulthood. Considering that cattle are physiologically similar to humans with respect to gestational length, our study may serve as a model for studies related to the developmental origin of health and disease in humans.

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.

Genome Expression Dynamics Reveal the Parasitism Regulatory Landscape of the Root-Knot Nematode Meloidogyne incognita and a Promoter Motif Associated with Effector Genes.

Root-knot nematodes (genus Meloidogyne) are the major contributor to crop losses caused by nematodes. These nematodes secrete effector proteins into the plant, derived from two sets of pharyngeal gland cells, to manipulate host physiology and immunity. Successful completion of the life cycle, involving successive molts from egg to adult, covers morphologically and functionally distinct stages and will require precise control of gene expression, including effector genes. The details of how root-knot nematodes regulate transcription remain sparse. Here, we report a life stage-specific transcriptome of Meloidogyne incognita. Combined with an available annotated genome, we explore the spatio-temporal regulation of gene expression. We reveal gene expression clusters and predicted functions that accompany the major developmental transitions. Focusing on effectors, we identify a putative cis-regulatory motif associated with expression in the dorsal glands, providing an insight into effector regulation. We combine the presence of this motif with several other criteria to predict a novel set of putative dorsal gland effectors. Finally, we show this motif, and thereby its utility, is broadly conserved across the Meloidogyne genus, and we name it Mel-DOG. Taken together, we provide the first genome-wide analysis of spatio-temporal gene expression in a root-knot nematode and identify a new set of candidate effector genes that will guide future functional analyses.

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.

Comparative Genomics Reveals Novel Target Genes towards Specific Control of Plant-Parasitic Nematodes.

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

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 F1 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.