A Comparison of transgenic and wild type soybean seeds: analysis of transcriptome profiles using RNA-Seq

In Situ Proteomic Analysis of Herbicide-Resistant Soybean and Hybrid Seeds via Matrix-Assisted Laser Desorption/Ionization-Mass Spectrometry Imaging

Transgenic soybean is the commercial crop with the largest cultivation area worldwide. During transgenic soybean cultivation, exogenous genes may be transferred to wild relatives through gene flow, posing unpredictable ecological risks. Accordingly, an environmental risk assessment should focus on fitness changes and underlying mechanisms in hybrids between transgenic and wild soybeans (Glycine soja). Matrix-assisted laser desorption/ionization-mass spectrometry imaging (MALDI-MSI) was used for in situ detection and imaging of protein changes in the seeds of transgenic herbicide-resistant soybean harboring epsps and pat genes, non-transgenic soybean, wild soybean, and their F₂ hybrid. Protein data clearly distinguished wild soybeans, while the F₂ seeds had protein characteristics of both parents and were distinguished from wild soybean seeds. Using UPLC-Q-TOF-MS, 22 differentially expressed proteins (DEPs) were identified, including 13 specific to wild soybean. Sucrose synthase and stress response-related DEPs were differentially expressed in parental and offspring. Differences in these may underpin the greater adaptability of the latter. MSI revealed DEP distribution in transgenic, wild, and F₂ seeds. Identifying DEPs related to fitness may elucidate mechanisms underlying fitness differences among the studied varieties. Our study shows that MALDI-MSI has the potential to become a visual method for transgenic soybean analysis.

Phospholipases C and D and Their Role in Biotic and Abiotic Stresses

Plants, as sessile organisms, have adapted a fine sensing system to monitor environmental changes, therefore allowing the regulation of their responses. As the interaction between plants and environmental changes begins at the surface, these changes are detected by components in the plasma membrane, where a molecule receptor generates a lipid signaling cascade via enzymes, such as phospholipases (PLs). Phospholipids are the key structural components of plasma membranes and signaling cascades. They exist in a wide range of species and in different proportions, with conversion processes that involve hydrophilic enzymes, such as phospholipase-C (PLC), phospholipase-D (PLD), and phospholipase-A (PLA). Hence, it is suggested that PLC and PLD are highly conserved, compared to their homologous genes, and have formed clusters during their adaptive history. Additionally, they generate responses to different functions in accordance with their protein structure, which should be reflected in specific signal transduction responses to environmental stress conditions, including innate immune responses. This review summarizes the phospholipid systems associated with signaling pathways and the innate immune response.

Citrus Tristeza Virus Genotype Detection Using High-Throughput Sequencing

The application of high-throughput sequencing (HTS) has successfully been used for virus discovery to resolve disease etiology in many agricultural crops. The greatest advantage of HTS is that it can provide a complete viral status of a plant, including information on mixed infections of viral species or virus variants. This provides insight into the virus population structure, ecology, or evolution and can be used to differentiate among virus variants that may contribute differently toward disease etiology. In this study, the use of HTS for citrus tristeza virus (CTV) genotype detection was evaluated. A bioinformatic pipeline for CTV genotype detection was constructed and evaluated using simulated and real data sets to determine the parameters to discriminate between false positive read mappings and true genotype-specific genome coverage. A 50% genome coverage cut-off was identified for non-target read mappings. HTS with the associated bioinformatic pipeline was validated and proposed as a CTV genotyping assay.

Negligible transcriptome and metabolome alterations in RNAi insecticidal maize against Monolepta hieroglyphica

RNAi-based genetically modified maize resistant to Monolepta hieroglyphica (Motschulsky) was demonstrated with negligible transcriptome and metabolome alterations compared to its unmodified equivalent. As one of the most prevalent insect pests afflicting various crops, Monolepta hieroglyphica (Motschulsky) causes severe loss of agricultural and economic productivity for many years in China. In an effort to reduce damages, in this study, an RNA interference (RNAi)-based genetically modified (GM) maize was developed. It was engineered to produce MhSnf7 double-stranded RNAs (dsRNAs), which can suppress the Snf7 gene expression and then lead M. hieroglyphica to death. Field trail analysis confirmed the robustly insecticidal ability of the MhSnf7 GM maize to resist damages by M. hieroglyphica. RNA sequencing analysis identified that only one gene was differentially expressed in the MhSnf7 GM maize compared to non-GM maize, indicating that the transcriptome in MhSnf7 GM maize is principally unaffected by the introduction of the MhSnf7 dsRNA expression vector. Likewise, metabolomics analysis identified that only 8 out of 5787 metabolites were significantly changed. Hence, the integration of transcriptomics and metabolomics demonstrates that there are negligible differences between MhSnf7 GM maize and its unmodified equivalent. This study not only presents a comprehensive assessment of cellular alteration in terms of gene transcription and metabolite abundance in RNAi-based GM maize, but also could be used as a reference for evaluating the unintended effect of GM crops.

Mapping Quantitative Trait Loci for Soybean Seedling Shoot and Root Architecture Traits in an Inter-Specific Genetic Population

Wild soybean species (Glycine soja Siebold & Zucc.) comprise a unique resource to widen the genetic base of cultivated soybean [Glycine max (L.) Merr.] for various agronomic traits. An inter-specific mapping population derived from a cross of cultivar Williams 82 and PI 483460B, a wild soybean accession, was utilized for genetic characterization of root architecture traits. The objectives of this study were to identify and characterize quantitative trait loci (QTL) for seedling shoot and root architecture traits, as well as to determine additive/epistatic interaction effects of identified QTLs. A total of 16,469 single nucleotide polymorphisms (SNPs) developed for the Illumina beadchip genotyping platform were used to construct a high resolution genetic linkage map. Among the 11 putative QTLs identified, two significant QTLs on chromosome 7 were determined to be associated with total root length (RL) and root surface area (RSA) with favorable alleles from the wild soybean parent. These seedling root traits, RL (BARC_020495_04641 ~ BARC_023101_03769) and RSA (SNP02285 ~ SNP18129_Magellan), could be potential targets for introgression into cultivated soybean background to improve both tap and lateral roots. The RL QTL region harbors four candidate genes with higher expression in root tissues: Phosphofructokinase (Glyma.07g126400), Snf7 protein (Glyma.07g127300), unknown functional gene (Glyma.07g127900), and Leucine Rich-Repeat protein (Glyma.07g127100). The novel alleles inherited from the wild soybean accession could be used as molecular markers to improve root system architecture and productivity in elite soybean lines.

Comprehensive transcriptomics and proteomics analyses of rice stripe virus-resistant transgenic rice

Stable transgenic rice line (named KRSV-1) with strong resistance against rice stripe virus was generated using the gene sequence of disease-specific protein by RNA interference. Comprehensive safety assessment of transgenic plants has turned into a significant field of genetic modification food safety. In this study, a safety assessment of KRSV-1 was carried out in a stepwise approach. The molecular analysis exhibited that KRSV-1 harbored one copy number of transgene, which was integrated into the intergenic non-coding region of chromosome 2 associated with inter-chromosomal translocations of 1.6-kb segments of chromosome 8. Then, transcriptomics and proteomics analyses were carried out to detect the unintended effects as a result of the integration of the transgene. Although 650 dramatically differentially expressed genes (DDEGs) and 357 differentially expressed proteins were detected between KRSV-1 and wild-type (WT) by transcriptomics and proteomics analyses, no harmful members in the form of toxic proteins and allergens were observed. Encouragingly, the nutritional compositions of seeds from KRSV-1 were comparable with WT seeds. The results of this entire study of molecular analysis, transcriptome and proteome profile of KRSV-1 revealed that no detrimental changes in the form of toxic proteins and allergens were detected in the transgenic rice line due to the integration of the transgene.

Genome-wide transcriptome profiling of transgenic hop (Humulus lupulus L.) constitutively overexpressing HlWRKY1 and HlWDR1 transcription factors

Background

The hop plant (Humulus lupulus L.) is a valuable source of several secondary metabolites, such as flavonoids, bitter acids, and essential oils. These compounds are widely implicated in the beer brewing industry and are having potential biomedical applications. Several independent breeding programs around the world have been initiated to develop new cultivars with enriched lupulin and secondary metabolite contents but met with limited success due to several constraints. In the present work, a pioneering attempt has been made to overexpress master regulator binary transcription factor complex formed by HlWRKY1 and HlWDR1 using a plant expression vector to enhance the level of prenylflavonoid and bitter acid content in the hop. Subsequently, we performed transcriptional profiling using high-throughput RNA-Seq technology in leaves of resultant transformants and wild-type hop to gain in-depth information about the genome-wide functional changes induced by HlWRKY1 and HlWDR1 overexpression.

Results

The transgenic WW-lines exhibited an elevated expression of structural and regulatory genes involved in prenylflavonoid and bitter acid biosynthesis pathways. In addition, the comparative transcriptome analysis revealed a total of 522 transcripts involved in 30 pathways, including lipids and amino acids biosynthesis, primary carbon metabolism, phytohormone signaling and stress responses were differentially expressed in WW-transformants. It was apparent from the whole transcriptome sequencing that modulation of primary carbon metabolism and other pathways by HlWRKY1 and HlWDR1 overexpression resulted in enhanced substrate flux towards secondary metabolites pathway. The detailed analyses suggested that none of the pathways or genes, which have a detrimental effect on physiology, growth and development processes, were induced on a genome-wide scale in WW-transgenic lines.

Conclusions

Taken together, our results suggest that HlWRKY1 and HlWDR1 simultaneous overexpression positively regulates the prenylflavonoid and bitter acid biosynthesis pathways in the hop and thus these transgenes are presented as prospective candidates for achieving enhanced secondary metabolite content in the hop.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5125-8) contains supplementary material, which is available to authorized users.

The importance of genotype identity, genetic heterogeneity, and bioinformatic handling for properly assessing genomic variation in transgenic plants

BACKGROUND

The advent of -omics technologies has enabled the resolution of fine molecular differences among individuals within a species. DNA sequence variations, such as single nucleotide polymorphisms or small deletions, can be tabulated for many kinds of genotype comparisons. However, experimental designs and analytical approaches are replete with ways to overestimate the level of variation present within a given sample. Analytical pipelines that do not apply proper thresholds nor assess reproducibility among samples are susceptible to calling false-positive variants. Furthermore, issues with sample genotype identity or failing to account for heterogeneity in reference genotypes may lead to misinterpretations of standing variants as polymorphisms derived de novo.

RESULTS

A recent publication that featured the analysis of RNA-sequencing data in three transgenic soybean event series appeared to overestimate the number of sequence variants identified in plants that were exposed to a tissue culture based transformation process. We reanalyzed these data with a stringent set of criteria and demonstrate three different factors that lead to variant overestimation, including issues related to the genetic identity of the background genotype, unaccounted genetic heterogeneity in the reference genome, and insufficient bioinformatics filtering.

CONCLUSIONS

This study serves as a cautionary tale to users of genomic and transcriptomic data that wish to assess the molecular variation attributable to tissue culture and transformation processes. Moreover, accounting for the factors that lead to sequence variant overestimation is equally applicable to samples derived from other germplasm sources, including chemical or irradiation mutagenesis and genome engineering (e.g., CRISPR) processes.

Comparative Transcriptome Analyses Reveal Potential Mechanisms of Enhanced Drought Tolerance in Transgenic Salvia Miltiorrhiza Plants Expressing AtDREB1A from Arabidopsis

In our previous study, drought-resistant transgenic plants of Salvia miltiorrhiza were produced via overexpression of the transcription factor AtDREB1A. To unravel the molecular mechanisms underpinning elevated drought tolerance in transgenic plants, in the present study we compared the global transcriptional profiles of wild-type (WT) and AtDREB1A-expressing transgenic plants using RNA-sequencing (RNA-seq). Using cluster analysis, we identified 3904 differentially expressed genes (DEGs). Compared with WT plants, 423 unigenes were up-regulated in pRD29A::AtDREB1A-31 before drought treatment, while 936 were down-regulated and 1580 and 1313 unigenes were up- and down-regulated after six days of drought. COG analysis revealed that the 'signal transduction mechanisms' category was highly enriched among these DEGs both before and after drought stress. Based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) annotation, DEGs associated with "ribosome", "plant hormone signal transduction", photosynthesis", "plant-pathogen interaction", "glycolysis/gluconeogenesis" and "carbon fixation" are hypothesized to perform major functions in drought resistance in AtDREB1A-expressing transgenic plants. Furthermore, the number of DEGs associated with different transcription factors increased significantly after drought stress, especially the AP2/ERF, bZIP and MYB protein families. Taken together, this study substantially expands the transcriptomic information for S. miltiorrhiza and provides valuable clues for elucidating the mechanism of AtDREB1A-mediated drought tolerance in transgenic plants.

Transcript Polymorphism Rates in Soybean Seed Tissue Are Increased in a Single Transformant of Glycine max

Transgenic crops have been utilized for decades to enhance agriculture and more recently have been applied as bioreactors for manufacturing pharmaceuticals. Recently, we investigated the gene expression profiles of several in-house transgenic soybean events, finding one transformant group to be consistently different from our controls. In the present study, we examined polymorphisms and sequence variations in the exomes of the same transgenic soybean events. We found that the previously dissimilar soybean line also exhibited markedly increased levels of polymorphisms within mRNA transcripts from seed tissue, many of which are classified as gene expression modifiers. The results from this work will direct future investigations to examine novel SNPs controlling traits of great interest for breeding and improving transgenic soybean crops. Further, this study marks the first work to investigate SNP rates in transgenic soybean seed tissues and demonstrates that while transgenesis may induce abundant unanticipated changes in gene expression and nucleotide variation, phenotypes and overall health of the plants examined remained unaltered.

De novo transcriptome analysis of Spodoptera exigua multiple nucleopolyhedrovirus (SeMNPV) genes in latently infected Se301 cells

Cells of the P8-Se301-C1 strain are Spodoptera exigua cell clones that each harbor a partial version of the S. exigua multiple nucleopolyhedrovirus (SeMNPV) genome and which are resistant to homologous SeMNPV infections. The cells produce no viral progeny, suggesting that the infection is a latent-like viral infection. To investigate the SeMNPV genes harbored in the P8-Se301-C1 cells, the de novo transcriptomes of P8-Se301-C1 cells and S. exigua Se301 cells were analyzed and compared. A total of 54,569,296 reads were obtained from the P8-Se301-C1 cells that yielded 112,565 final unigenes with a mean length of 1,093 nt. A total of 56,865,504 reads were obtained from the Se301 cells that yielded 102,996 final unigenes with a mean length of 1,082 nt. Ten SeMNPV gene transcripts (se5, se7, se8, se12, se43, se45, se89, se90, se124, and se126) were detected in the P8-Se301-C1 cells by RNA-Seq but not in the Se301 cells, which was verified by RTPCR. 5'/3' RACE analyses showed that the 3'- or 5'-end sequences of the viral transcripts are aligned to the host gene sequences in P8-Se301-C1 cells, suggesting that the SeMNPV genes may integrate into and be transcribed with the host genes in the P8-Se301-C1 cells. Furthermore, six additional viral gene transcripts, se11, se42, se44, se88, se91, and se127 (incorporated into chimeric fusion transcripts in the P8-Se301-C1 cells), were detected in the RACE analyses. Taken together, sixteen SeMNPV transcripts were identified in the P8-Se301-C1 cell strain. This study provides information to develop the understanding of baculovirus latent infections and superinfection exclusion.