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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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

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.

MicroRNA Maturation and MicroRNA Target Gene Expression Regulation Are Severely Disrupted in Soybean dicer-like1 Double Mutants.

Small nonprotein-coding microRNAs (miRNAs) are present in most eukaryotes and are central effectors of RNA silencing-mediated mechanisms for gene expression regulation. In plants, DICER-LIKE1 (DCL1) is the founding member of a highly conserved family of RNase III-like endonucleases that function as core machinery proteins to process hairpin-like precursor transcripts into mature miRNAs, small regulatory RNAs, 21-22 nucleotides in length. Zinc finger nucleases (ZFNs) were used to generate single and double-mutants of putative soybean DCL1 homologs, DCL1a and DCL1b, to confirm their functional role(s) in the soybean miRNA pathway. Neither DCL1 single mutant, dcl1a or dcl1b plants, exhibited a pronounced morphological or molecular phenotype. However, the dcl1a/dcl1b double mutant expressed a strong morphological phenotype, characterized by reduced seed size and aborted seedling development, in addition to defective miRNA precursor transcript processing efficiency and deregulated miRNA target gene expression. Together, these findings indicate that the two soybean DCL1 paralogs, DCL1a and DCL1b, largely play functionally redundant roles in the miRNA pathway and are essential for normal plant development.

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.

Identical substitutions in magnesium chelatase paralogs result in chlorophyll-deficient soybean mutants.

The soybean [Glycine max (L.) Merr.] chlorophyll-deficient line MinnGold is a spontaneous mutant characterized by yellow foliage. Map-based cloning and transgenic complementation revealed that the mutant phenotype is caused by a nonsynonymous nucleotide substitution in the third exon of a Mg-chelatase subunit gene (ChlI1a) on chromosome 13. This gene was selected as a candidate for a different yellow foliage mutant, T219H (Y11y11), that had been previously mapped to chromosome 13. Although the phenotypes of MinnGold and T219H are clearly distinct, sequencing of ChlI1a in T219H identified a different nonsynonymous mutation in the third exon, only six base pairs from the MinnGold mutation. This information, along with previously published allelic tests, were used to identify and clone a third yellow foliage mutation, CD-5, which was previously mapped to chromosome 15. This mutation was identified in the ChlI1b gene, a paralog of ChlI1a. Sequencing of the ChlI1b allele in CD-5 identified a nonsynonymous substitution in the third exon that confers an identical amino acid change as the T219H substitution at ChlI1a. Protein sequence alignments of the two Mg-chelatase subunits indicated that the sites of amino acid modification in MinnGold, T219H, and CD-5 are highly conserved among photosynthetic species. These results suggest that amino acid alterations in this critical domain may create competitive inhibitory interactions between the mutant and wild-type ChlI1a and ChlI1b proteins.

Targeted mutagenesis for functional analysis of gene duplication in legumes.

Assessment of gene function oftentimes requires mutant populations that can be screened by forward or reverse genetic analysis. The situation becomes more complicated in polyploidy or paleopolyploid genomes that have two or more copies for most genes. Here we describe a method for engineering zinc-finger nucleases (ZFNs) for the purpose of creating targeted mutations in the paleopolyploid soybean genome. ZFNs are recombinant proteins composed of an engineered zinc-finger array fused to a nonspecific cleavage domain. When engineered to recognize a specific nucleotide sequence, the cleavage domain will generate highly mutagenic DNA double-strand breaks frequently resulting in insertions and deletions at the target locus. Depending on the number of target sites present within the genome, this method has the capacity to target either single- or multi-copy gene families. In this chapter, we describe an inexpensive, rapid, and user-friendly approach for ZFN assembly and application in soybean based on the previously described context-dependent assembly method.

Isolation and analysis of small RNAs from virus-infected plants.

In this chapter, we detail some of the methods available to the researcher for isolating and analyzing virus-derived small RNAs (vsRNAs). These methods have been successfully used for four plant viruses: Cucumber mosaic virus (CMV), including the CMV Y-Satellite, Turnip mosaic virus (TuMV), Potato leaf roll virus (PLRV), and Tomato spotted wilt virus (TSWV) from inoculated Arabidopsis thaliana plants (Fusaro et al. EMBO Rep 7:1168-1175, 2006; Curtin et al. FEBS Lett 582:2753-2760, 2008). The protocols presented here can also be employed for the isolation of non-virus related small RNAs such as microRNAs (miRNAs) and hairpin RNA (hpRNA).

DRB2 is required for microRNA biogenesis in Arabidopsis thaliana.

BACKGROUND: The Arabidopsis thaliana (Arabidopsis) DOUBLE-STRANDED RNA BINDING (DRB) protein family consists of five members, DRB1 to DRB5. The biogenesis of two developmentally important small RNA (sRNA) species, the microRNAs (miRNAs) and trans-acting small interfering RNAs (tasiRNAs) by DICER-LIKE (DCL) endonucleases requires the assistance of DRB1 and DRB4 respectively. The importance of miRNA-directed target gene expression in plant development is exemplified by the phenotypic consequence of loss of DRB1 activity (drb1 plants). PRINCIPAL FINDINGS: Here we report that the developmental phenotype of the drb235 triple mutant plant is the result of deregulated miRNA biogenesis in the shoot apical meristem (SAM) region. The expression of DRB2, DRB3 and DRB5 in wild-type seedlings is restricted to the SAM region. Small RNA sequencing of the corresponding tissue of drb235 plants revealed altered miRNA accumulation. Approximately half of the miRNAs detected remained at levels equivalent to those of wild-type plants. However, the accumulation of the remaining miRNAs was either elevated or reduced in the triple mutant. Examination of different single and multiple drb mutants revealed a clear association between the loss of DRB2 activity and altered accumulation for both the elevated and reduced miRNA classes. Furthermore, we show that the constitutive over-expression of DRB2 outside of its wild-type expression domain can compensate for the loss of DRB1 activity in drb1 plants. CONCLUSIONS/SIGNIFICANCE: Our results suggest that in the SAM region, DRB2 is both antagonistic and synergistic to the role of DRB1 in miRNA biogenesis, adding an additional layer of gene regulatory complexity in this developmentally important tissue.

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.