Genome assembly of Medicago truncatula accession SA27063 provides insight into spring black stem and leaf spot disease resistance.

Medicago truncatula, model legume and alfalfa relative, has served as an essential resource for advancing our understanding of legume physiology, functional genetics, and crop improvement traits. Necrotrophic fungus, Ascochyta medicaginicola, the causal agent of spring black stem (SBS) and leaf spot is a devasting foliar disease of alfalfa affecting stand survival, yield, and forage quality. Host resistance to SBS disease is poorly understood, and control methods rely on cultural practices. Resistance has been observed in M. truncatula accession SA27063 (HM078) with two recessively inherited quantitative-trait loci (QTL), rnpm1 and rnpm2, previously reported. To shed light on host resistance, we carried out a de novo genome assembly of HM078. The genome, referred to as MtHM078 v1.0, is comprised of 23 contigs totaling 481.19 Mbp. Notably, this assembly contains a substantial amount of novel centromere-related repeat sequences due to deep long-read sequencing. Genome annotation resulted in 98.4% of BUSCO fabales proteins being complete. The assembly enabled sequence-level analysis of rnpm1 and rnpm2 for gene content, synteny, and structural variation between SBS-resistant accession SA27063 (HM078) and SBS-susceptible accession A17 (HM101). Fourteen candidate genes were identified, and some have been implicated in resistance to necrotrophic fungi. Especially interesting candidates include loss-of-function events in HM078 because they fit the inverse gene-for-gene model, where resistance is recessively inherited. In rnpm1, these include a loss-of-function in a disease resistance gene due to a premature stop codon, and a 10.85 kbp retrotransposon-like insertion disrupting a ubiquitin conjugating E2. In rnpm2, we identified a frameshift mutation causing a loss-of-function in a glycosidase, as well as a missense and frameshift mutation altering an F-box family protein. This study generated a high-quality genome of HM078 and has identified promising candidates, that once validated, could be further studied in alfalfa to enhance disease resistance.

SELF PRUNING 3C is a flowering repressor that modulates seed germination, root architecture, and drought responses.

Allelic variation in the CETS (CENTRORADIALIS, TERMINAL FLOWER 1, SELF PRUNING) gene family controls agronomically important traits in many crops. CETS genes encode phosphatidylethanolamine-binding proteins that have a central role in the timing of flowering as florigenic and anti-florigenic signals. The great expansion of CETS genes in many species suggests that the functions of this family go beyond flowering induction and repression. Here, we characterized the tomato SELF PRUNING 3C (SP3C) gene, and show that besides acting as a flowering repressor it also regulates seed germination and modulates root architecture. We show that loss of SP3C function in CRISPR/Cas9-generated mutant lines increases root length and reduces root side branching relative to the wild type. Higher SP3C expression in transgenic lines promotes the opposite effects in roots, represses seed germination, and also improves tolerance to water stress in seedlings. These discoveries provide new insights into the role of SP paralogs in agronomically relevant traits, and support future exploration of the involvement of CETS genes in abiotic stress responses.

A Multipurpose Toolkit to Enable Advanced Genome Engineering in Plants.

We report a comprehensive toolkit that enables targeted, specific modification of monocot and dicot genomes using a variety of genome engineering approaches. Our reagents, based on transcription activator-like effector nucleases (TALENs) and the clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 system, are systematized for fast, modular cloning and accommodate diverse regulatory sequences to drive reagent expression. Vectors are optimized to create either single or multiple gene knockouts and large chromosomal deletions. Moreover, integration of geminivirus-based vectors enables precise gene editing through homologous recombination. Regulation of transcription is also possible. A Web-based tool streamlines vector selection and construction. One advantage of our platform is the use of the Csy-type (CRISPR system yersinia) ribonuclease 4 (Csy4) and tRNA processing enzymes to simultaneously express multiple guide RNAs (gRNAs). For example, we demonstrate targeted deletions in up to six genes by expressing 12 gRNAs from a single transcript. Csy4 and tRNA expression systems are almost twice as effective in inducing mutations as gRNAs expressed from individual RNA polymerase III promoters. Mutagenesis can be further enhanced 2.5-fold by incorporating the Trex2 exonuclease. Finally, we demonstrate that Cas9 nickases induce gene targeting at frequencies comparable to native Cas9 when they are delivered on geminivirus replicons. The reagents have been successfully validated in tomato (Solanum lycopersicum), tobacco (Nicotiana tabacum), Medicago truncatula, wheat (Triticum aestivum), and barley (Hordeum vulgare).

CRISPR/Cas9 and TALENs generate heritable mutations for genes involved in small RNA processing of Glycine max and Medicago truncatula.

Processing of double-stranded RNA precursors into small RNAs is an essential regulator of gene expression in plant development and stress response. Small RNA processing requires the combined activity of a functionally diverse group of molecular components. However, in most of the plant species, there are insufficient mutant resources to functionally characterize each encoding gene. Here, mutations in loci encoding protein machinery involved in small RNA processing in soya bean and Medicago truncatula were generated using the CRISPR/Cas9 and TAL-effector nuclease (TALEN) mutagenesis platforms. An efficient CRISPR/Cas9 reagent was used to create a bi-allelic double mutant for the two soya bean paralogous Double-stranded RNA-binding2 (GmDrb2a and GmDrb2b) genes. These mutations, along with a CRISPR/Cas9-generated mutation of the M. truncatula Hua enhancer1 (MtHen1) gene, were determined to be germ-line transmissible. Furthermore, TALENs were used to generate a mutation within the soya bean Dicer-like2 gene. CRISPR/Cas9 mutagenesis of the soya bean Dicer-like3 gene and the GmHen1a gene was observed in the T0 generation, but these mutations failed to transmit to the T1 generation. The irregular transmission of induced mutations and the corresponding transgenes was investigated by whole-genome sequencing to reveal a spectrum of non-germ-line-targeted mutations and multiple transgene insertion events. Finally, a suite of combinatorial mutant plants were generated by combining the previously reported Gmdcl1a, Gmdcl1b and Gmdcl4b mutants with the Gmdrb2ab double mutant. Altogether, this study demonstrates the synergistic use of different genome engineering platforms to generate a collection of useful mutant plant lines for future study of small RNA processing in legume crops.

Validating Genome-Wide Association Candidates Controlling Quantitative Variation in Nodulation.

Genome-wide association (GWA) studies offer the opportunity to identify genes that contribute to naturally occurring variation in quantitative traits. However, GWA relies exclusively on statistical association, so functional validation is necessary to make strong claims about gene function. We used a combination of gene-disruption platforms (Tnt1 retrotransposons, hairpin RNA-interference constructs, and CRISPR/Cas9 nucleases) together with randomized, well-replicated experiments to evaluate the function of genes that an earlier GWA study in Medicago truncatula had identified as candidates contributing to variation in the symbiosis between legumes and rhizobia. We evaluated ten candidate genes found in six clusters of strongly associated single nucleotide polymorphisms, selected on the basis of their strength of statistical association, proximity to annotated gene models, and root or nodule expression. We found statistically significant effects on nodule production for three candidate genes, each validated in two independent mutants. Annotated functions of these three genes suggest their contributions to quantitative variation in nodule production occur through processes not previously connected to nodulation, including phosphorous supply and salicylic acid-related defense response. These results demonstrate the utility of GWA combined with reverse mutagenesis technologies to discover and validate genes contributing to naturally occurring variation in quantitative traits. The results highlight the potential for GWA to complement forward genetics in identifying the genetic basis of ecologically and economically important traits.

Genomic variation and DNA repair associated with soybean transgenesis: a comparison to cultivars and mutagenized plants.

BACKGROUND: The safety of mutagenized and genetically transformed plants remains a subject of scrutiny. Data gathered and communicated on the phenotypic and molecular variation induced by gene transfer technologies will provide a scientific-based means to rationally address such concerns. In this study, genomic structural variation (e.g. large deletions and duplications) and single nucleotide polymorphism rates were assessed among a sample of soybean cultivars, fast neutron-derived mutants, and five genetically transformed plants developed through Agrobacterium based transformation methods. RESULTS: On average, the number of genes affected by structural variations in transgenic plants was one order of magnitude less than that of fast neutron mutants and two orders of magnitude less than the rates observed between cultivars. Structural variants in transgenic plants, while rare, occurred adjacent to the transgenes, and at unlinked loci on different chromosomes. DNA repair junctions at both transgenic and unlinked sites were consistent with sequence microhomology across breakpoints. The single nucleotide substitution rates were modest in both fast neutron and transformed plants, exhibiting fewer than 100 substitutions genome-wide, while inter-cultivar comparisons identified over one-million single nucleotide polymorphisms. CONCLUSIONS: Overall, these patterns provide a fresh perspective on the genomic variation associated with high-energy induced mutagenesis and genetically transformed plants. The genetic transformation process infrequently results in novel genetic variation and these rare events are analogous to genetic variants occurring spontaneously, already present in the existing germplasm, or induced through other types of mutagenesis. It remains unclear how broadly these results can be applied to other crops or transformation methods.

Selection-free zinc-finger-nuclease engineering by context-dependent assembly (CoDA).

Engineered zinc-finger nucleases (ZFNs) enable targeted genome modification. Here we describe context-dependent assembly (CoDA), a platform for engineering ZFNs using only standard cloning techniques or custom DNA synthesis. Using CoDA-generated ZFNs, we rapidly altered 20 genes in Danio rerio, Arabidopsis thaliana and Glycine max. The simplicity and efficacy of CoDA will enable broad adoption of ZFN technology and make possible large-scale projects focused on multigene pathways or genome-wide alterations.

Ideotype breeding and genome engineering for legume crop improvement.

Ideotype breeding is a strategy whereby traits are modeled a priori and then introduced into a model or crop species to assess their impact on yield. Thus, knowledge about the connection between genotype and phenotype is required for ideotype breeding to be deployed successfully. The growing understanding of the genetic basis of yield-related traits, combined with increasingly efficient genome engineering tools, improved transformation efficiency, and high-throughput genotyping of regenerants paves the way for the widespread adoption of ideotype breeding as a complement to conventional breeding. We briefly discuss how ideotype breeding, coupled with such state-of-the-art biotechnological tools, could contribute to knowledge-based legume breeding and accelerate yield gains to ensure food security in the coming decades.

Optimization of in vitro and ex vitro Agrobacterium rhizogenes-mediated hairy root transformation of soybean for visual screening of transformants using RUBY.

In vitro and ex vitro Agrobacterium rhizogenes-mediated hairy root transformation (HRT) assays are key components of the plant biotechnology and functional genomics toolkit. In this report, both in vitro and ex vitro HRT were optimized in soybean using the RUBY reporter. Different parameters including A. rhizogenes strain, optical density of the bacterial cell culture (OD600), co-cultivation media, soybean genotype, explant age, and acetosyringone addition and concentration were evaluated. Overall, the in vitro assay was more efficient than the ex vitro assay in terms of the percentage of induction of hairy roots and transformed roots (expressing RUBY). Nonetheless, the ex vitro technique was deemed faster and a less complicated approach. The highest transformation of RUBY was observed on 7-d-old cotyledons of cv. Bert inoculated for 30 minutes with the R1000 resuspended in » B5 medium to OD600 (0.3) and 150 µM of acetosyringone. The parameters of this assay also led to the highest percentage of RUBY through two-step ex vitro hairy root transformation. Finally, using machine learning-based modeling, optimal protocols for both assays were further defined. This study establishes efficient and reliable hairy root transformation protocols applicable for functional studies in soybean.

Medicago truncatula PHO2 genes have distinct roles in phosphorus homeostasis and symbiotic nitrogen fixation.

Three PHO2-like genes encoding putative ubiquitin-conjugating E2 enzymes of Medicago truncatula were characterized for potential roles in phosphorous (P) homeostasis and symbiotic nitrogen fixation (SNF). All three genes, MtPHO2A, B and C, contain miR399-binding sites characteristic of PHO2 genes in other plant species. Distinct spatiotemporal expression patterns and responsiveness of gene expression to P- and N-deprivation in roots and shoots indicated potential roles, especially for MtPHO2B, in P and N homeostasis. Phenotypic analysis of pho2 mutants revealed that MtPHO2B is integral to Pi homeostasis, affecting Pi allocation during plant growth under nutrient-replete conditions, while MtPHO2C had a limited role in controlling Pi homeostasis. Genetic analysis also revealed a connection between Pi allocation, plant growth and SNF performance. Under N-limited, SNF conditions, Pi allocation to different organs was dependent on MtPHO2B and, to a lesser extent, MtPHO2C and MtPHO2A. MtPHO2A also affected Pi homeostasis associated with nodule formation. Thus, MtPHO2 genes play roles in systemic and localized, i.e., nodule, P homeostasis affecting SNF.

Co-expression of soybean Dicer-like genes in response to stress and development.

Regulation of gene transcription and post-transcriptional processes is critical for proper development, genome integrity, and stress responses in plants. Many genes involved in the key processes of transcriptional and post-transcriptional regulation have been well studied in model diploid organisms. However, gene and genome duplication may alter the function of the genes involved in these processes. To address this question, we assayed the stress-induced transcription patterns of duplicated gene pairs involved in RNAi and DNA methylation processes in the paleopolyploid soybean. Real-time quantitative PCR and Sequenom MassARRAY expression assays were used to profile the relative expression ratios of eight gene pairs across eight different biotic and abiotic stress conditions. The transcriptional responses to stress for genes involved in DNA methylation, RNAi processing, and miRNA processing were compared. The strongest evidence for pairwise co-expression in response to stresses was exhibited by non-paralogous Dicer-like (DCL) genes GmDCL2a-GmDCL3a and GmDCL1b-GmDCL2b, most profoundly in root tissues. Among homoeologous or paralogous DCL genes, the Dicer-like 2 (DCL2) gene pair exhibited the strongest response to stress and most conserved co-expression pattern. This was surprising because the DCL2 duplication event is more ancient than the other DCL duplications. Possible mechanisms that may be driving the DCL2 co-expression are discussed.

Targeted mutagenesis of duplicated genes in soybean with zinc-finger nucleases.

We performed targeted mutagenesis of a transgene and nine endogenous soybean (Glycine max) genes using zinc-finger nucleases (ZFNs). A suite of ZFNs were engineered by the recently described context-dependent assembly platform--a rapid, open-source method for generating zinc-finger arrays. Specific ZFNs targeting dicer-like (DCL) genes and other genes involved in RNA silencing were cloned into a vector under an estrogen-inducible promoter. A hairy-root transformation system was employed to investigate the efficiency of ZFN mutagenesis at each target locus. Transgenic roots exhibited somatic mutations localized at the ZFN target sites for seven out of nine targeted genes. We next introduced a ZFN into soybean via whole-plant transformation and generated independent mutations in the paralogous genes DCL4a and DCL4b. The dcl4b mutation showed efficient heritable transmission of the ZFN-induced mutation in the subsequent generation. These findings indicate that ZFN-based mutagenesis provides an efficient method for making mutations in duplicate genes that are otherwise difficult to study due to redundancy. We also developed a publicly accessible Web-based tool to identify sites suitable for engineering context-dependent assembly ZFNs in the soybean genome.

Alfalfa (Medicago sativa L.) pho2 mutant plants hyperaccumulate phosphate.

In this article, we describe a set of novel alfalfa (Medicago sativa L.) plants that hyper-accumulate Phosphate ion (Pi) at levels 3- to 6-fold higher than wild-type. This alfalfa germplasm will have practical applications reclaiming Pi from contaminated or enriched soil or be used in conservation buffer strips to protect waterways from Pi run-off. Hyper-accumulating alfalfa plants were generated by targeted mutagenesis of PHOSPHATE2 (PHO2) using newly created CRISPR/Cas9 reagents and an improved mutant screening strategy. PHO2 encodes a ubiquitin conjugating E2 enzyme (UBC24) previously characterized in Arabidopsis thaliana, Medicago truncatula, and Oryza sativa. Mutations of PHO2 disrupt Pi homeostasis resulting in Pi hyper-accumulation. Successful CRISPR/Cas9 editing of PHO2 demonstrates that this is an efficient mutagenesis tool in alfalfa despite its complex autotetraploid genome structure. Arabidopsis and M. truncatula ortholog genes were used to identify PHO2 haplotypes in outcrossing tetraploid M. sativa with the aim of generating heritable mutations in both PHO2-like genes (PHO2-B and PHO2-C). After delivery of the reagent and regeneration from transformed leaf explants, plants with mutations in all haplotypes of PHO2-B and PHO2-C were identified. These plants were evaluated for morphology, Pi accumulation, heritable transmission of targeted mutations, segregation of mutant haplotypes and removal of T-DNA(s). The Agrobacterium-mediated transformation assay and gene editing reagents reported here were also evaluated for further optimization for future alfalfa functional genomic studies.

The Arabidopsis thaliana double-stranded RNA binding protein DRB1 directs guide strand selection from microRNA duplexes.

In Arabidopsis thaliana (Arabidopsis), DICER-LIKE1 (DCL1) functions together with the double-stranded RNA binding protein (dsRBP), DRB1, to process microRNAs (miRNAs) from their precursor transcripts prior to their transfer to the RNA-induced silencing complex (RISC). miRNA-loaded RISC directs RNA silencing of cognate mRNAs via ARGONAUTE1 (AGO1)-catalyzed cleavage. Short interefering RNAs (siRNAs) are processed from viral-derived or transgene-encoded molecules of double-stranded RNA (dsRNA) by the DCL/dsRBP partnership, DCL4/DRB4, and are also loaded to AGO1-catalyzed RISC for cleavage of complementary mRNAs. Here, we use an artificial miRNA (amiRNA) technology, transiently expressed in Nicotiana benthamiana, to produce a series of amiRNA duplexes with differing intermolecular thermostabilities at the 5' end of duplex strands. Analyses of amiRNA duplex strand accumulation and target transcript expression revealed that strand selection (amiRNA and amiRNA*) is directed by asymmetric thermostability of the duplex termini. The duplex strand possessing a lower 5' thermostability was preferentially retained by RISC to guide mRNA cleavage of the corresponding target transgene. In addition, analysis of endogenous miRNA duplex strand accumulation in Arabidopsis drb1 and drb2345 mutant plants revealed that DRB1 dictates strand selection, presumably by directional loading of the miRNA duplex onto RISC for passenger strand degradation. Bioinformatic and Northern blot analyses of DCL4/DRB4-dependent small RNAs (miRNAs and siRNAs) revealed that small RNAs produced by this DCL/dsRBP combination do not conform to the same terminal thermostability rules as those governing DCL1/DRB1-processed miRNAs. This suggests that small RNA processing in the DCL1/DRB1-directed miRNA and DCL4/DRB4-directed sRNA biogenesis pathways operates via different mechanisms.

Potato improvement through genetic engineering.

Potato (Solanum tuberosum L.) is the third most important crop worldwide and a staple food for many people worldwide. Genetically, it poses many challenges for traditional breeding due to its autotetraploid nature and its tendency toward inbreeding depression. Breeding programs have focused on productivity, nutritional quality, and disease resistance. Some of these traits exist in wild potato relatives but their introgression into elite cultivars can take many years and, for traits such as pest resistance, their effect is often short-lasting. These problems can be addressed by genetic modification (GM) or gene editing (GE) and open a wide horizon for potato crop improvement. Current genetically modified and gene edited varieties include those with Colorado potato beetle and late blight resistance, reduction in acrylamide, and modified starch content. RNAi hairpin technology can be used to silence the haplo-alleles of multiple genes simultaneously, whereas optimization of newer gene editing technologies such as base and prime editing will facilitate the routine generation of advanced edits across the genome. These technologies will likely gain further relevance as increased target specificity and decreased off-target effects are demonstrated. In this Review, we discuss recent work related to these technologies in potato improvement.

The roles of plant dsRNA-binding proteins in RNAi-like pathways.

Dicers are associated with double-stranded RNA-binding proteins (dsRBPs) in animals. In the plant, Arabidopsis, there are four dicer-like (DCL) proteins and five potential dsRBPs. These DCLs act redundantly and hierarchically. However, we show there is little or no redundancy or hierarchy amongst the DRBs in their DCL interactions. DCL1 operates exclusively with DRB1 to produce micro (mi)RNAs, DCL4 operates exclusively with DRB4 to produce trans-acting (ta) siRNAs and 21nt siRNAs from viral RNA. DCL2 and DCL3 produce viral siRNAs without requiring assistance from any dsRBP. DRB2, DRB3 and DRB5 appear unnecessary for mi-, tasi-, viral si-, or heterochromatinising siRNA production but act redundantly in a developmental pathway.

Editing the Medicago truncatula Genome: Targeted Mutagenesis Using the CRISPR-Cas9 Reagent.

Medicago truncatula is an annual plant used for studying legume biology, in particular symbioses with nitrogen-fixing rhizobia and arbuscular mycorrhizal fungi. Efforts to decipher the genetic basis of these ecologically and economically important traits are a major goal of plant and crop biology. M. truncatula is an excellent model system for this purpose, as it has several publicly available sequenced genomes, has a rapid seed-to-seed generation time, and is highly transformable. Various mutagenesis platforms such as Tnt1 retrotransposons and RNAi knockdown have been used successfully in forward and reverse genetic studies to identify and functionally characterize candidate genes. The CRISPR/Cas9 reagent is the most recent mutagenesis platform and is highly effective at generating site-directed double-stranded breaks (DSB) in M. truncatula. This protocol will demonstrate the construction of reagents using two genome engineering platforms that have successfully generated mutant plants in M. truncatula, M. sativa, and soybean systems. The reagents are easy to assemble, can be quickly retrofitted to test novel regulatory sequences for improved efficiency, and can be used for more advanced genome engineering strategies such as gene insertion or gene replacement.

Design and Assembly of CRISPR/Cas9 Reagents for Gene Knockout, Targeted Insertion, and Replacement in Wheat.

Advances in cereal transformation along with the completion of the wheat genome sequence assembly have increased the demand for tools that perform targeted and specific modifications in this crop plant. This protocol demonstrates the construction of reagents using a comprehensive genome engineering kit to create single and multiple gene "knockouts," site-specific chromosome deletions and gene replacement or "knockins" including the use of geminivirus replicons (GVRs). The reagents allow for both easy construction of simple genome engineering vectors, and "mix and match" swapping of components such as the Cas9, guide RNA and donor template cassettes for gene targeting. In addition, a web-based tool greatly streamlines vector selection, primer design, and vector construction.

The evolution and diversification of Dicers in plants.

Most multicellular organisms regulate developmental transitions by microRNAs, which are generated by an enzyme, Dicer. Insects and fungi have two Dicer-like genes, and many animals have only one, yet the plant, Arabidopsis, has four. Examining the poplar and rice genomes revealed that they contain five and six Dicer-like genes, respectively. Analysis of these genes suggests that plants require a basic set of four Dicer types which were present before the divergence of mono- and dicotyledonous plants ( approximately 200 million years ago), but after the divergence of plants from green algae. A fifth type of Dicer seems to have evolved in monocots.

CRISPR/Cas mutagenesis of soybean and Medicago truncatula using a new web-tool and a modified Cas9 enzyme.

The CRISPR/Cas9 system is rapidly becoming the reagent of choice for targeted mutagenesis and gene editing in crop species. There are currently intense research efforts in the crop sciences to identify efficient CRISPR/Cas9 platforms to carry out targeted mutagenesis and gene editing projects. These efforts typically result in the incremental tweaking of various platform components including the identification of crop-specific promoters and terminators for optimal expression of the Cas9 enzyme and identification of promoters for expression of the CRISPR guide RNA. In this report, we demonstrate the development of an online web tool for fast identification of CRISPR/Cas9 target loci within soybean gene models, and generic DNA sequences. The web-tool described in this work can quickly identify a high number of potential CRISPR/Cas9 target sites, including restriction enzyme sites that can facilitate the detection of new mutations. In conjunction with the web tool, a soybean codon-optimized CRISPR/Cas9 platform was designed to direct double-stranded breaks to the targeted loci in hairy root transformed cells. The modified Cas9 enzyme was shown to successfully mutate target genes in somatic cells of 2 legume species, soybean and Medicago truncatula. These new tools may help facilitate targeted mutagenesis in legume and other plant species.