Zinc finger nuclease-mediated targeting of multiple transgenes to an endogenous soybean genomic locus via non-homologous end joining.

Emerging genome editing technologies hold great promise for the improvement of agricultural crops. Several related genome editing methods currently in development utilize engineered, sequence-specific endonucleases to generate DNA double strand breaks (DSBs) at user-specified genomic loci. These DSBs subsequently result in small insertions/deletions (indels), base substitutions or incorporation of exogenous donor sequences at the target site, depending on the application. Targeted mutagenesis in soybean (Glycine max) via non-homologous end joining (NHEJ)-mediated repair of such DSBs has been previously demonstrated with multiple nucleases, as has homology-directed repair (HDR)-mediated integration of a single transgene into target endogenous soybean loci using CRISPR/Cas9. Here we report targeted integration of multiple transgenes into a single soybean locus using a zinc finger nuclease (ZFN). First, we demonstrate targeted integration of biolistically delivered DNA via either HDR or NHEJ to the FATTY ACID DESATURASE 2-1a (FAD2-1a) locus of embryogenic cells in tissue culture. We then describe ZFN- and NHEJ-mediated, targeted integration of two different multigene donors to the FAD2-1a locus of immature embryos. The largest donor delivered was 16.2 kb, carried four transgenes, and was successfully transmitted to T1 progeny of mature targeted plants obtained via somatic embryogenesis. The insertions in most plants with a targeted, 7.1 kb, NHEJ-integrated donor were perfect or near-perfect, demonstrating that NHEJ is a viable alternative to HDR for gene targeting in soybean. Taken together, these results show that ZFNs can be used to generate fertile transgenic soybean plants with NHEJ-mediated targeted insertions of multigene donors at an endogenous genomic locus.

Expression characterization of the herbicide tolerance gene Aryloxyalkanoate Dioxygenase (aad-1) controlled by seven combinations of regulatory elements.

BACKGROUND: Availability of well characterized maize regulatory elements for gene expression in a variety of tissues and developmental stages provides effective alternatives for single and multigene transgenic concepts. We studied the expression of the herbicide tolerance gene aryloxyalkanoate dioxygenase (aad-1) driven by seven different regulatory element construct designs including the ubiquitin promoters of maize and rice, the actin promoters of melon and rice, three different versions of the Sugarcane Bacilliform Badnavirus promoters in association with other regulatory elements of gene expression. RESULTS: Gene expression of aad-1 was characterized at the transcript and protein levels in a collection of maize tissues and developmental stages. Protein activity against its target herbicide was characterized by herbicide dosage response. Although differences in transcript and protein accumulation were observed among the different constructs tested, all events were tolerant to commercially relevant rates of quizalafop-P-ethyl compared to non-traited maize under greenhouse conditions. DISCUSSION: The data reported demonstrate how different regulatory elements affect transcript and protein accumulation and how these molecular characteristics translate into the level of herbicide tolerance. The level of transcript detected did not reflect the amount of protein quantified in a particular tissue since protein accumulation may be influenced not only by levels of transcript produced but also by translation rate, post-translational regulation mechanisms and protein stability. The amount of AAD-1 enzyme produced with all constructs tested showed sufficient enzymatic activity to detoxify the herbicide and prevent most herbicidal damage at field-relevant levels without having a negative effect on plant health. CONCLUSIONS: Distinctive profiles of aad-1 transcript and protein accumulation were observed when different regulatory elements were utilized in the constructs under study. The ZmUbi and the SCBV constructs showed the most consistent robust tolerance, while the melon actin construct provided the lowest level of tolerance compared to the other regulatory elements used in this study. These data provide insights into the effects of differing levels of gene expression and how these molecular characteristics translate into the level of herbicide tolerance. Furthermore, these data provide valuable information to optimize future designs of single and multiple gene constructs for maize research and crop improvement.

Identification and use of the sugarcane bacilliform virus enhancer in transgenic maize.

BACKGROUND: Transcriptional enhancers are able to increase transcription from heterologous promoters when placed upstream, downstream and in either orientation, relative to the promoter. Transcriptional enhancers have been used to enhance expression of specific promoters in transgenic plants and in activation tagging studies to help elucidate gene function. RESULTS: A transcriptional enhancer from the Sugarcane Bacilliform Virus - Ireng Maleng isolate (SCBV-IM) that can cause increased transcription when integrated into the the genome near maize genes has been identified. In transgenic maize, the SCBV-IM promoter was shown to be comparable in strength to the maize ubiquitin 1 promoter in young leaf and root tissues. The promoter was dissected to identify sequences that confer high activity in transient assays. Enhancer sequences were identified and shown to increase the activity of a heterologous truncated promoter. These enhancer sequences were shown to be more active when arrayed in 4 copy arrays than in 1 or 2 copy arrays. When the enhancer array was transformed into maize plants it caused an increase in accumulation of transcripts of genes near the site of integration in the genome. CONCLUSIONS: The SCBV-IM enhancer can activate transcription upstream or downstream of genes and in either orientation. It may be a useful tool to activate enhance from specific promoters or in activation tagging.

Flux of transcript patterns during soybean seed development.

BACKGROUND: To understand gene expression networks leading to functional properties of the soybean seed, we have undertaken a detailed examination of soybean seed development during the stages of major accumulation of oils, proteins, and starches, as well as the desiccating and mature stages, using microarrays consisting of up to 27,000 soybean cDNAs. A subset of these genes on a highly-repetitive 70-mer oligonucleotide microarray was also used to support the results. RESULTS: It was discovered that genes related to cell growth and maintenance processes, as well as energy processes like photosynthesis, decreased in expression levels as the cotyledons approached the mature, dry stage. Genes involved with some storage proteins had their highest expression levels at the stage of highest fresh weight. However, genes encoding many transcription factors and DNA binding proteins showed higher expression levels in the desiccating and dry seeds than in most of the green stages. CONCLUSIONS: Data on 27,000 cDNAs have been obtained over five stages of soybean development, including the stages of major accumulation of agronomically-important products, using two different types of microarrays. Of particular interest are the genes found to peak in expression at the desiccating and dry seed stages, such as those annotated as transcription factors, which may indicate the preparation of pathways that will be needed later in the early stages of imbibition and germination.

Improved sensitivity of DNA microarrays using photonic crystal enhanced fluorescence.

DNA microarrays are used to profile changes in gene expression between samples in a high-throughput manner, but measurements of genes with low expression levels can be problematic with standard microarray substrates. In this work, we expand the detection capabilities of a standard microarray experiment using a photonic crystal (PC) surface that enhances fluorescence observed from microarray spots. This PC is inexpensively and uniformly fabricated using a nanoreplica molding technique, with very little variation in its optical properties within- and between-devices. By using standard protocols to process glass microarray substrates in parallel with PCs, we evaluated the impact of this substrate on a one-color microarray experiment comparing gene expression in two developmental stages of Glycine max. The PCs enhanced the signal-to-noise ratio observed from microarray spots by 1 order of magnitude, significantly increasing the number of genes detected above substrate fluorescence noise. PC substrates more than double the number of genes classified as differentially expressed, detecting changes in expression even for low expression genes. This approach increases the dynamic range of a surface-bound fluorescence-based assay to reliably quantify small quantities of DNA that would be impossible with standard substrates.

Targeted insertion of large DNA sequences by homology-directed repair or non-homologous end joining in engineered tobacco BY-2 cells using designed zinc finger nucleases.

Targeted integration of recombinant DNA fragments into plant genomes by DNA double-strand break (DSB) repair mechanisms has become a powerful tool for precision engineering of crops. However, many targeting platforms require the screening of many transgenic events to identify a low number of targeted events among many more random insertion events. We developed an engineered transgene integration platform (ETIP) that uses incomplete marker genes at the insertion site to enable rapid phenotypic screening and recovery of targeted events upon functional reconstitution of the marker genes. The two marker genes, encoding neomycin phosphotransferase II (nptII) and Discosoma sp. red fluorescent protein (DsRed) enable event selection on kanamycin-containing selective medium and subsequent screening for red fluorescent clones. The ETIP design allows targeted integration of donor DNA molecules either by homology-directed repair (HDR) or non-homologous end joining (NHEJ)-mediated mechanisms. Targeted donor DNA integration is facilitated by zinc finger nucleases (ZFN). The ETIP cassette was introduced into Nicotiana tabacum BY-2 suspension cells to generate target cell lines containing a single copy locus of the transgene construct. The utility of the ETIP platform has been demonstrated by targeting DNA constructs containing up to 25-kb payload. The success rate for clean targeted DNA integration was up to 21% for HDR and up to 41% for NHEJ based on the total number of calli analyzed by next-generation sequencing (NGS). The rapid generation of targeted events with large DNA constructs expands the utility of the nuclease-mediated gene addition platform both for academia and the commercial sector.

Specific elements of the glyoxylate pathway play a significant role in the functional transition of the soybean cotyledon during seedling development.

BACKGROUND: The soybean (Glycine max) cotyledon is a specialized tissue whose main function is to serve as a nutrient reserve that supplies the needs of the young plant throughout seedling development. During this process the cotyledons experience a functional transition to a mainly photosynthetic tissue. To identify at the genetic level the specific active elements that participate in the natural transition of the cotyledon from storage to photosynthetic activity, we studied the transcript abundance profile at different time points using a new soybean oligonucleotide chip containing 19,200 probes (70-mer long). RESULTS: After normalization and statistical analysis we determined that 3,594 genes presented a statistically significant altered expression in relation to the imbibed seed in at least one of the time points defined for the study. Detailed analysis of this data identified individual, specific elements of the glyoxylate pathway that play a fundamental role during the functional transition of the cotyledon from nutrient storage to photosynthesis. The dynamics between glyoxysomes and peroxisomes is evident during these series of events. We also identified several other genes whose products could participate co-ordinately throughout the functional transition and the associated mechanisms of control and regulation and we described multiple unknown genetic elements that by association have the potential to make a major contribution to this biological process. CONCLUSION: We demonstrate that the global transcript profile of the soybean cotyledon during seedling development is extremely active, highly regulated and dynamic. We defined the expression profiles of individual gene family members, enzymatic isoforms and protein subunits and classified them accordingly to their involvement in different functional activities relevant to seedling development and the cotyledonary functional transition in soybean, especially the ones associated with the glyoxylate cycle. Our data suggests that in the soybean cotyledon a very complex and synchronized system of control and regulation of several metabolic pathways is essential to carry out the necessary functions during this developmental process.

Microarray analysis of iron deficiency chlorosis in near-isogenic soybean lines.

BACKGROUND: Iron is one of fourteen mineral elements required for proper plant growth and development of soybean (Glycine max L. Merr.). Soybeans grown on calcareous soils, which are prevalent in the upper Midwest of the United States, often exhibit symptoms indicative of iron deficiency chlorosis (IDC). Yield loss has a positive linear correlation with increasing severity of chlorotic symptoms. As soybean is an important agronomic crop, it is essential to understand the genetics and physiology of traits affecting plant yield. Soybean cultivars vary greatly in their ability to respond successfully to iron deficiency stress. Microarray analyses permit the identification of genes and physiological processes involved in soybean's response to iron stress. RESULTS: RNA isolated from the roots of two near isogenic lines, which differ in iron efficiency, PI 548533 (Clark; iron efficient) and PI 547430 (IsoClark; iron inefficient), were compared on a spotted microarray slide containing 9,728 cDNAs from root specific EST libraries. A comparison of RNA transcripts isolated from plants grown under iron limiting hydroponic conditions for two weeks revealed 43 genes as differentially expressed. A single linkage clustering analysis of these 43 genes showed 57% of them possessed high sequence similarity to known stress induced genes. A control experiment comparing plants grown under adequate iron hydroponic conditions showed no differences in gene expression between the two near isogenic lines. Expression levels of a subset of the differentially expressed genes were also compared by real time reverse transcriptase PCR (RT-PCR). The RT-PCR experiments confirmed differential expression between the iron efficient and iron inefficient plants for 9 of 10 randomly chosen genes examined. To gain further insight into the iron physiological status of the plants, the root iron reductase activity was measured in both iron efficient and inefficient genotypes for plants grown under iron sufficient and iron limited conditions. Iron inefficient plants failed to respond to decreased iron availability with increased activity of Fe reductase. CONCLUSION: These experiments have identified genes involved in the soybean iron deficiency chlorosis response under iron deficient conditions. Single linkage cluster analysis suggests iron limited soybeans mount a general stress response as well as a specialized iron deficiency stress response. Root membrane bound reductase capacity is often correlated with iron efficiency. Under iron-limited conditions, the iron efficient plant had high root bound membrane reductase capacity while the iron inefficient plants reductase levels remained low, further limiting iron uptake through the root. Many of the genes up-regulated in the iron inefficient NIL are involved in known stress induced pathways. The most striking response of the iron inefficient genotype to iron deficiency stress was the induction of a profusion of signaling and regulatory genes, presumably in an attempt to establish and maintain cellular homeostasis. Genes were up-regulated that point toward an increased transport of molecules through membranes. Genes associated with reactive oxidative species and an ROS-defensive enzyme were also induced. The up-regulation of genes involved in DNA repair and RNA stability reflect the inhospitable cellular environment resulting from iron deficiency stress. Other genes were induced that are involved in protein and lipid catabolism; perhaps as an effort to maintain carbon flow and scavenge energy. The under-expression of a key glycolitic gene may result in the iron-inefficient genotype being energetically challenged to maintain a stable cellular environment. These experiments have identified candidate genes and processes for further experimentation to increase our understanding of soybeans' response to iron deficiency stress.

Transcriptome changes in the phenylpropanoid pathway of Glycine max in response to Pseudomonas syringae infection.

BACKGROUND: Reports of plant molecular responses to pathogenic infections have pinpointed increases in activity of several genes of the phenylpropanoid pathway leading to the synthesis of lignin and flavonoids. The majority of those findings were derived from single gene studies and more recently from several global gene expression analyses. We undertook a global transcriptional analysis focused on the response of genes of the multiple branches of the phenylpropanoid pathway to infection by the Pseudomonas syringae pv. glycinea with or without the avirulence gene avrB to characterize more broadly the contribution of the multiple branches of the pathway to the resistance response in soybean. Transcript abundance in leaves was determined from analysis of soybean cDNA microarray data and hybridizations to RNA blots with specific gene probes. RESULTS: The majority of the genes surveyed presented patterns of increased transcript accumulation. Some increased rapidly, 2 and 4 hours after inoculation, while others started to accumulate slowly by 8-12 hours. In contrast, transcripts of a few genes decreased in abundance 2 hours post inoculation. Most interestingly was the opposite temporal fluctuation in transcript abundance between early responsive genes in defense (CHS and IFS1) and F3H, the gene encoding a pivotal enzyme in the synthesis of anthocyanins, proanthocyanidins and flavonols. F3H transcripts decreased rapidly 2 hours post inoculation and increased during periods when CHS and IFS transcripts decreased. It was also determined that all but one (CHS4) family member genes (CHS1, CHS2, CHS3, CHS5, CHS6 and CHS7/8) accumulated higher transcript levels during the defense response provoked by the avirulent pathogen challenge. CONCLUSION: Based on the mRNA profiles, these results show the strong bias that soybean has towards increasing the synthesis of isoflavonoid phytoalexins concomitant with the down regulation of genes required for the synthesis of anthocyanins and proanthocyanins. Although proanthocyanins are known to be toxic compounds, the cells in the soybean leaves seem to be programmed to prioritize the synthesis and accumulation of isoflavonoid and pterocarpan phytoalexins during the resistance response. It was known that CHS transcripts accumulate in great abundance rapidly after inoculation of the soybean plants but our results have demonstrated that all but one (CHS4) member of the gene family member genes accumulated higher transcript levels during the defense response.

Expression profiling soybean response to Pseudomonas syringae reveals new defense-related genes and rapid HR-specific downregulation of photosynthesis.

Transcript profiling during susceptible (S) and hypersensitive response-associated resistance (R) interactions was determined in soybean (Glycine max). Pseudomonas syringae pv. glycinea carrying or lacking the avirulence gene avrB, was infiltrated into cultivar Williams 82. Leaf RNA was sampled at 2, 8, and 24 h postinoculation (hpi). Significant changes in transcript abundance were observed for 3,897 genes during the experiment at P < or = 0.000005. Many of the genes showed a similar direction of increase or decrease in abundance in both the S and R responses, but the R response generally showed a significantly greater degree of differential expression. More than 25% of these responsive genes had not been previously reported as being associated with pathogen interactions, as 704 had no functional annotation and 378 had no homolog in National Center for Biotechnology Information databases. The highest number of transcriptional changes was noted at 8 hpi, including the downregulation of 94 chloroplast-associated genes specific to the R response. Photosynthetic measurements were consistent with an R-specific reduction in photosystem II operating efficiency (phiPSII) that was apparent at 8 hpi for the R response with little effect in the S or control treatments. Imaging analyses suggest that the decreased phiPSII was a result of physical damage to PSII reaction centers.

Nuclear genes that encode mitochondrial proteins for DNA and RNA metabolism are clustered in the Arabidopsis genome.

The plant mitochondrial genome is complex in structure, owing to a high degree of recombination activity that subdivides the genome and increases genetic variation. The replication activity of various portions of the mitochondrial genome appears to be nonuniform, providing the plant with an ability to modulate its mitochondrial genotype during development. These and other interesting features of the plant mitochondrial genome suggest that adaptive changes have occurred in DNA maintenance and transmission that will provide insight into unique aspects of plant mitochondrial biology and mitochondrial-chloroplast coevolution. A search in the Arabidopsis genome for genes involved in the regulation of mitochondrial DNA metabolism revealed a region of chromosome III that is unusually rich in genes for mitochondrial DNA and RNA maintenance. An apparently similar genetic linkage was observed in the rice genome. Several of the genes identified within the chromosome III interval appear to target the plastid or to be targeted dually to the mitochondria and the plastid, suggesting that the process of endosymbiosis likely is accompanied by an intimate coevolution of these two organelles for their genome maintenance functions.