Weighted gene co-expression network analysis reveals immune evasion related genes in Echinococcus granulosus sensu stricto.

Cystic echinococcosis (CE) is a zoonotic disease caused by the tapeworm Echinococcus granulosus sensu lato (s.l). In the intermediate host, this disease is characterized by the growth of cysts in viscera such as liver and lungs, inside of which the parasite develops to the next infective stage known as protoscoleces. There are records that the infected viscera affect the development and morphology of E. granulosus s.l. protoscolex in hosts such as buffalo or humans. However, the molecular mechanisms that drive these differences remains unknown. Weighted gene co-expression network analysis (WGCNA) using a set of RNAseq data obtained from E. granulosus sensu stricto (s.s.) protoscoleces found in liver and lung cysts reveals 34 modules in protoscoleces of liver origin, of which 12 have differential co-expression from protoscoleces of lung origin. Three of these twelve modules contain hub genes related to immune evasion: tegument antigen, tegumental protein, ubiquitin hydrolase isozyme L3, COP9 signalosome complex subunit 3, tetraspanin CD9 antigen, and the methyl-CpG-binding protein Mbd2. Also, two of the twelve modules contain only hypothetical proteins with unknown orthology, which means that there are a group of unknown function proteins co-expressed inside the protoscolex of liver CE cyst origin. This is the first evidence of gene expression differences in protoscoleces from CE cysts found in different viscera, with co-expression networks that are exclusive to protoscoleces from liver CE cyst samples. This should be considered in the control strategies of CE, as intermediate hosts can harbor CE cysts in liver, lungs, or both organs simultaneously.

A novel pectin methylesterase inhibitor, PMEI3, in common bean suggests a key role of pectin methylesterification in Pseudomonas resistance.

The mechanisms underlying susceptibility to and defense against Pseudomonas syringae (Pph) of the common bean (Phaseolus vulgaris) have not yet been clarified. To investigate these, 15-day-old plants of the variety Riñón were infected with Pph and the transcriptomic changes at 2 h and 9 h post-infection were analysed. RNA-seq analysis showed an up-regulation of genes involved in defense/signaling at 2 h, most of them being down-regulated at 9 h, suggesting that Pph inhibits the transcriptomic reprogramming of the plant. This trend was also observed in the modulation of 101 cell wall-related genes. Cell wall composition changes at early stages of Pph infection were associated with homogalacturonan methylation and the formation of egg boxes. Among the cell wall genes modulated, a pectin methylesterase inhibitor 3 (PvPMEI3) gene, closely related to AtPMEI3, was detected. PvPMEI3 protein was located in the apoplast and its pectin methylesterase inhibitory activity was demonstrated. PvPMEI3 seems to be a good candidate to play a key role in Pph infection, which was supported by analysis of an Arabidopsis pmei3 mutant, which showed susceptibility to Pph, in contrast to resistant Arabidopsis Col-0 plants. These results indicate a key role of the degree of pectin methylesterification in host resistance to Pph during the first steps of the attack.

DNA sequencing in the classroom: complete genome sequence of two earwig (Dermaptera; Insecta) species.

BACKGROUND: Despite representing the largest fraction of animal life, the number of insect species whose genome has been sequenced is barely in the hundreds. The order Dermaptera (the earwigs) suffers from a lack of genomic information despite its unique position as one of the basally derived insect groups and its importance in agroecosystems. As part of a national educational and outreach program in genomics, a plan was formulated to engage the participation of high school students in a genome sequencing project. Students from twelve schools across Chile were instructed to capture earwig specimens in their geographical area, to identify them and to provide material for genome sequencing to be carried out by themselves in their schools. RESULTS: The school students collected specimens from two cosmopolitan earwig species: Euborellia annulipes (Fam. Anisolabididae) and Forficula auricularia (Fam. Forficulidae). Genomic DNA was extracted and, with the help of scientific teams that traveled to the schools, was sequenced using nanopore sequencers. The sequence data obtained for both species was assembled and annotated. We obtained genome sizes of 1.18 Gb (F. auricularia) and 0.94 Gb (E. annulipes) with the number of predicted protein coding genes being 31,800 and 40,000, respectively. Our analysis showed that we were able to capture a high percentage (≥ 93%) of conserved proteins indicating genomes that are useful for comparative and functional analysis. We were also able to characterize structural elements such as repetitive sequences and non-coding RNA genes. Finally, functional categories of genes that are overrepresented in each species suggest important differences in the process underlying the formation of germ cells, and modes of reproduction between them, features that are one of the distinguishing biological properties that characterize these two distant families of Dermaptera. CONCLUSIONS: This work represents an unprecedented instance where the scientific and lay community have come together to collaborate in a genome sequencing project. The versatility and accessibility of nanopore sequencers was key to the success of the initiative. We were able to obtain full genome sequences of two important and widely distributed species of insects which had not been analyzed at this level previously. The data made available by the project should illuminate future studies on the Dermaptera.

DNA sequencing in the classroom: complete genome sequence of two earwig (Dermaptera; Insecta) species

BACKGROUND: Despite representing the largest fraction of animal life, the number of insect species whose genome has been sequenced is barely in the hundreds. The order Dermaptera (the earwigs) suffers from a lack of genomic information despite its unique position as one of the basally derived insect groups and its importance in agroecosystems. As part of a national educational and outreach program in genomics, a plan was formulated to engage the participation of high school students in a genome sequencing project. Students from twelve schools across Chile were instructed to capture earwig specimens in their geographical area, to identify them and to provide material for genome sequencing to be carried out by themselves in their schools. RESULTS: The school students collected specimens from two cosmopolitan earwig species: Euborellia annulipes (Fam. Anisolabididae) and Forficula auricularia (Fam. Forficulidae). Genomic DNA was extracted and, with the help of scientific teams that traveled to the schools, was sequenced using nanopore sequencers. The sequence data obtained for both species was assembled and annotated. We obtained genome sizes of 1.18 Gb (F. auricularia) and 0.94 Gb (E. annulipes) with the number of predicted protein coding genes being 31,800 and 40,000, respectively. Our analysis showed that we were able to capture a high percentage (≥ 93%) of conserved proteins indicating genomes that are useful for comparative and functional analysis. We were also able to characterize structural elements such as repetitive sequences and non-coding RNA genes. Finally, functional categories of genes that are overrepresented in each species suggest important differences in the process underlying the formation of germ cells, and modes of reproduction between them, features that are one of the distinguishing biological properties that characterize these two distant families of Dermaptera. CONCLUSIONS: This work represents an unprecedented instance where the scientific and lay community have come together to collaborate in a genome sequencing project. The versatility and accessibility of nanopore sequencers was key to the success of the initiative. We were able to obtain full genome sequences of two important and widely distributed species of insects which had not been analyzed at this level previously. The data made available by the project should illuminate future studies on the Dermaptera.

Transcriptome and Gene Regulatory Network Analyses Reveal New Transcription Factors in Mature Fruit Associated with Harvest Date in Prunus persica.

Harvest date is a critical parameter for producers and consumers regarding agro-industrial performance. It involves a pleiotropic effect controlling the development of other fruit quality traits through finely controlling regulatory mechanisms. Fruit ripening is a process in which various signals and biological events co-occur and are regulated by hormone signaling that produces the accumulation/degradation of multiple compounds. However, the regulatory mechanisms that control the hormone signaling involved in fruit development and ripening are still unclear. To investigate the issue, we used individuals with early, middle and late harvest dates from a peach segregating population to identify regulatory candidate genes controlling fruit quality traits at the harvest stage and validate them in contrasting peach varieties for this trait. We identified 467 and 654 differentially expressed genes for early and late harvest through a transcriptomic approach. In addition, using the Arabidopsis DAP-seq database and network analysis, six transcription factors were selected. Our results suggest significant hormonal balance and cell wall composition/structure differences between early and late harvest samples. Thus, we propose that higher expression levels of the transcription factors HB7, ERF017 and WRKY70 in early harvest individuals would induce the expression of genes associated with the jasmonic acid pathway, photosynthesis and gibberellins inhibition. While on the other hand, the high expression levels of LHY, CDF3 and NAC083 in late harvest individuals would promote the induction of genes associated with abscisic acid biosynthesis, auxins and cell wall remodeling.

Response Mechanisms of "Hass" Avocado to Sequential 1-methylcyclopropene Applications at Different Maturity Stages during Cold Storage.

1-Methylcyclopropene (1-MCP) is used for extending the postharvest life of the avocado during storage. Evaluated the effect of 1-MCP application at different times after harvest, i.e., 0, 7, 14, and 21 d at 5 °C, to identify the threshold of the ethylene inhibition response in "Hass" avocado. Our results showed that fruits from two maturity stages at harvest: low dry matter (20-23%) and high dry matter (27%). Changes in ethylene production rates and transcript accumulation of genes involved in ethylene metabolism were measured at harvest and during storage. 1-MCP treated fruit up to 14 d of storage showed similar values of firmness and skin color as fruit treated at harvest time. In contrast, when the application was performed after 21 d, the fruit showed ripening attributes similar to those of the untreated ones. To further understand the molecular mechanisms responsible for the lack of response to 1-MCP at 21 d of storage, transcriptomic analysis was performed. Gene ontology analyses based on the DEG analysis showed enrichment of transcripts involved in the 'response to ethylene' for both maturity stages. All genes evaluated showed similar expression profiles induced by cold storage time, with a peak at 21 d of storage and an increased softening of the fruit and peel color. This was a two-year field study, and results were consistent across the two experimental years. Our results should help growers and markets in selecting the optimal timing of 1-MCP application in "Hass" avocados and should contribute to a deeper understanding of the molecular mechanisms of the avocado ripening process.

Transcriptome analysis of Echinococcus granulosus sensu stricto protoscoleces reveals differences in immune modulation gene expression between cysts found in cattle and sheep.

Cystic Echinococcosis (CE), a zoonotic parasitic disease, is caused by the cestode Echinococcus granulosus sensu lato. CE inflicts severe damage in cattle, sheep, and human hosts worldwide. Fertile CE cysts are characterized by the presence of viable protoscoleces. These parasite forms are studied with minimal contamination with host molecules. Hosts, cattle and sheep, show differences in their CE cyst fertility. The effect of the host in protoscolex transcriptome is not known. We genotyped and performed transcriptomic analysis on sheep protoscoleces obtained from liver and lung CE cysts. The transcriptomic data of Echinococcus granulosus sensu stricto protoscoleces from 6 lung CE cysts and 6 liver CE cysts were Collected. For host comparison analysis, 4 raw data files belonging to Echinococcus granulosus sensu stricto protoscoleces from cattle liver CE cysts were obtained from the NCBI SRA database. Principal component and differential expression analysis did not reveal any statistical differences between protoscoleces obtained from liver or lung cysts, either within the same sheep or different sheep hosts. Conversely, there are significant differences between cattle and sheep protoscolex samples. We found differential expression of immune-related genes. In cattle, 7 genes were upregulated in protoscoleces from liver cysts. In sheep, 3 genes were upregulated in protoscoleces from liver and lung CE cysts. Noteworthy, are the differential expression of antigen B, tegument antigen, and arginase-2 in samples obtained from sheep CE cysts, and basigin in samples from cattle CE cysts. These findings suggest that the host species is an important factor involved in the differential expression of immune related genes, which in turn is possibly related to the fertility of Echinococcus granulosus sensu stricto cysts.

Unravelling the Molecular Regulation Mechanisms of Slow Ripening Trait in Prunus persica.

Fruit development is a complex process that involves the interplay of cell division, expansion, and differentiation. As a model to study fruit development, nectarines incapable of ripening were described as slow ripening. Slow ripening fruits remained firm and exhibited no rise in CO2 or ethylene production rates for one month or more at 20 °C. Different studies suggest that this trait is controlled by a single gene (NAC072). Transcriptome analysis between normal and slow ripening fruits showed a total of 157, 269, 976, and 5.224 differentially expressed genes in each fruit developmental stage analyzed (T1, T2, T3, and T7, respectively), and no expression of NAC072 was found in the slow ripening individuals. Using this transcriptomic information, we identified a correlation of NAC072 with auxin-related genes and two genes associated with terpene biosynthesis. On the other hand, significant differences were observed in hormonal biosynthetic pathways during fruit development between the normal and slow ripening individuals (gibberellin, ethylene, jasmonic acid and abscisic acid). These results suggest that the absence of NAC072 by the direct or indirect expression control of auxins or terpene-related genes prevents normal peach fruit development.

Identification of DNA Methylation and Transcriptomic Profiles Associated With Fruit Mealiness in Prunus persica (L.) Batsch.

Peach (Prunus persica) fruits have a fast ripening process and a shelf-life of days, presenting a challenge for long-distance consuming markets. To prolong shelf-life, peach fruits are stored at low temperatures (0 to 7 °C) for at least two weeks, which can lead to the development of mealiness, a physiological disorder that reduces fruit quality and decreases consumer acceptance. Several studies have been made to understand this disorder, however, the molecular mechanisms underlying mealiness are not fully understood. Epigenetic factors, such as DNA methylation, modulate gene expression according to the genetic background and environmental conditions. In this sense, the aim of this work was to identify differentially methylated regions (DMRs) that could affect gene expression in contrasting individuals for mealiness. Peach flesh was studied at harvest time (E1 stage) and after cold storage (E3 stage) for 30 days. The distribution of DNA methylations within the eight chromosomes of P. persica showed higher methylation levels in pericentromeric regions and most differences between mealy and normal fruits were at Chr1, Chr4, and Chr8. Notably, differences in Chr4 co-localized with previous QTLs associated with mealiness. Additionally, the number of DMRs was higher in CHH cytosines of normal and mealy fruits at E3; however, most DMRs were attributed to mealy fruits from E1, increasing at E3. From RNA-Seq data, we observed that differentially expressed genes (DEGs) between normal and mealy fruits were associated with ethylene signaling, cell wall modification, lipid metabolism, oxidative stress and iron homeostasis. When integrating the annotation of DMRs and DEGs, we identified a CYP450 82A and an UDP-ARABINOSE 4 EPIMERASE 1 gene that were downregulated and hypermethylated in mealy fruits, coinciding with the co-localization of a transposable element (TE). Altogether, this study indicates that genetic differences between tolerant and susceptible individuals is predominantly affecting epigenetic regulation over gene expression, which could contribute to a metabolic alteration from earlier stages of development, resulting in mealiness at later stages. Finally, this epigenetic mark should be further studied for the development of new molecular tools in support of breeding programs.

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

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

Dormant but Active: Chilling Accumulation Modulates the Epigenome and Transcriptome of Prunus avium During Bud Dormancy.

Temperate deciduous fruit tree species like sweet cherry (Prunus avium) require long periods of low temperatures to trigger dormancy release and flowering. In addition to sequence-based genetic diversity, epigenetic variation may contribute to different chilling requirements among varieties. For the low chill variety 'Royal Dawn' and high chill variety 'Kordia', we studied the methylome of floral buds during chilling accumulation using MethylC-seq to identify differentially methylated regions (DMRs) during chilling hours (CH) accumulation, followed by transcriptome analysis to correlate changes in gene expression with DNA methylation. We found that during chilling accumulation, DNA methylation increased from 173 CH in 'Royal Dawn' and 443 CH in 'Kordia' and was mostly associated with the CHH context. In addition, transcriptional changes were observed from 443 CH in 'Kordia' with 1,210 differentially expressed genes, increasing to 4,292 genes at 1,295 CH. While 'Royal Dawn' showed approximately 5,000 genes differentially expressed at 348 CH and 516 CH, showing a reprogramming that was specific for each genotype. From conserved upregulated genes that overlapped with hypomethylated regions and downregulated genes that overlapped with hypermethylated regions in both varieties, we identified genes related to cold-sensing, cold-signaling, oxidation-reduction process, metabolism of phenylpropanoids and lipids, and a MADS-box SVP-like gene. As a complementary analysis, we used conserved and non-conserved DEGs that presented a negative correlation between DNA methylations and mRNA levels across all chilling conditions, obtaining Gene Ontology (GO) categories related to abiotic stress, metabolism, and oxidative stress. Altogether, this data indicates that changes in DNA methylation precedes transcript changes and may occur as an early response to low temperatures to increase the cold tolerance in the endodormancy period, contributing with the first methylome information about the effect of environmental cues over two different genotypes of sweet cherry.

Expression QTL (eQTLs) Analyses Reveal Candidate Genes Associated With Fruit Flesh Softening Rate in Peach [Prunus persica (L.) Batsch].

Significant differences in softening rate have been reported between melting flesh in peach and nectarine varieties. This trait seems to be controlled by several genes. We aimed to identify candidate genes involved in fruit softening rate by integrating quantitative trait loci (QTL) and expression QTL (eQTL) analyses, comparing siblings with contrasting softening rates. We used a segregating population derived from nectarine cv. 'Venus' selfing, which was phenotyped for softening rate during three seasons. Six siblings with high (HSR) and six with low softening rate (LSR) were sequenced using RNA-Seq. A group of 5,041 differentially expressed genes was identified. Also, we found a QTL with a LOD (logarithm of odds) score of 9.7 on LG4 in all analyzed seasons. Furthermore, we detected 1,062 eQTLs, of which 133 were found co-localizing with the identified QTL. Gene Ontology (GO) analysis showed 'Response to auxin' as one the main over-represented categories. Our findings suggest over-expression of auxin biosynthetic related genes in the HSR group, which implies a higher expression and/or accumulation of auxin, thereby triggering fast softening. Conversely, the LSR phenotype might be explained by an altered auxin-homeostasis associated with low auxin levels. This work will contribute to unraveling the genetic mechanisms responsible for the softening rate in peaches and nectarines and lead to the development of molecular markers.