Current status of community resources and priorities for weed genomics research.

Weeds are attractive models for basic and applied research due to their impacts on agricultural systems and capacity to swiftly adapt in response to anthropogenic selection pressures. Currently, a lack of genomic information precludes research to elucidate the genetic basis of rapid adaptation for important traits like herbicide resistance and stress tolerance and the effect of evolutionary mechanisms on wild populations. The International Weed Genomics Consortium is a collaborative group of scientists focused on developing genomic resources to impact research into sustainable, effective weed control methods and to provide insights about stress tolerance and adaptation to assist crop breeding.

Miniature type V-F CRISPR-Cas nucleases enable targeted DNA modification in cells.

Class 2 CRISPR systems are exceptionally diverse, nevertheless, all share a single effector protein that contains a conserved RuvC-like nuclease domain. Interestingly, the size of these CRISPR-associated (Cas) nucleases ranges from >1000 amino acids (aa) for Cas9/Cas12a to as small as 400-600 aa for Cas12f. For in vivo genome editing applications, compact RNA-guided nucleases are desirable and would streamline cellular delivery approaches. Although miniature Cas12f effectors have been shown to cleave double-stranded DNA, targeted DNA modification in eukaryotic cells has yet to be demonstrated. Here, we biochemically characterize two miniature type V-F Cas nucleases, SpCas12f1 (497 aa) and AsCas12f1 (422 aa), and show that SpCas12f1 functions in both plant and human cells to produce targeted modifications with outcomes in plants being enhanced with short heat pulses. Our findings pave the way for the development of miniature Cas12f1-based genome editing tools.

Chromatin loop anchors contain core structural components of the gene expression machinery in maize.

BACKGROUND: Three-dimensional chromatin loop structures connect regulatory elements to their target genes in regions known as anchors. In complex plant genomes, such as maize, it has been proposed that loops span heterochromatic regions marked by higher repeat content, but little is known on their spatial organization and genome-wide occurrence in relation to transcriptional activity. RESULTS: Here, ultra-deep Hi-C sequencing of maize B73 leaf tissue was combined with gene expression and open chromatin sequencing for chromatin loop discovery and correlation with hierarchical topologically-associating domains (TADs) and transcriptional activity. A majority of all anchors are shared between multiple loops from previous public maize high-resolution interactome datasets, suggesting a highly dynamic environment, with a conserved set of anchors involved in multiple interaction networks. Chromatin loop interiors are marked by higher repeat contents than the anchors flanking them. A small fraction of high-resolution interaction anchors, fully embedded in larger chromatin loops, co-locate with active genes and putative protein-binding sites. Combinatorial analyses indicate that all anchors studied here co-locate with at least 81.5% of expressed genes and 74% of open chromatin regions. Approximately 38% of all Hi-C chromatin loops are fully embedded within hierarchical TAD-like domains, while the remaining ones share anchors with domain boundaries or with distinct domains. Those various loop types exhibit specific patterns of overlap for open chromatin regions and expressed genes, but no apparent pattern of gene expression. In addition, up to 63% of all unique variants derived from a prior public maize eQTL dataset overlap with Hi-C loop anchors. Anchor annotation suggests that < 7% of all loops detected here are potentially devoid of any genes or regulatory elements. The overall organization of chromatin loop anchors in the maize genome suggest a loop modeling system hypothesized to resemble phase separation of repeat-rich regions. CONCLUSIONS: Sets of conserved chromatin loop anchors mapping to hierarchical domains contains core structural components of the gene expression machinery in maize. The data presented here will be a useful reference to further investigate their function in regard to the formation of transcriptional complexes and the regulation of transcriptional activity in the maize genome.

Mapping regulatory variants controlling gene expression in drought response and tolerance in maize.

BACKGROUND: Gene expression is a key determinant of cellular response. Natural variation in gene expression bridges genetic variation to phenotypic alteration. Identification of the regulatory variants controlling the gene expression in response to drought, a major environmental threat of crop production worldwide, is of great value for drought-tolerant gene identification. RESULTS: A total of 627 RNA-seq analyses are performed for 224 maize accessions which represent a wide genetic diversity under three water regimes; 73,573 eQTLs are detected for about 30,000 expressing genes with high-density genome-wide single nucleotide polymorphisms, reflecting a comprehensive and dynamic genetic architecture of gene expression in response to drought. The regulatory variants controlling the gene expression constitutively or drought-dynamically are unraveled. Focusing on dynamic regulatory variants resolved to genes encoding transcription factors, a drought-responsive network reflecting a hierarchy of transcription factors and their target genes is built. Moreover, 97 genes are prioritized to associate with drought tolerance due to their expression variations through the Mendelian randomization analysis. One of the candidate genes, Abscisic acid 8'-hydroxylase, is verified to play a negative role in plant drought tolerance. CONCLUSIONS: This study unravels the effects of genetic variants on gene expression dynamics in drought response which allows us to better understand the role of distal and proximal genetic effects on gene expression and phenotypic plasticity. The prioritized drought-associated genes may serve as direct targets for functional investigation or allelic mining.

Developmental and transcriptional responses of maize to drought stress under field conditions.

Drought is a common abiotic stress which significantly limits global crop productivity. Maize is an important staple crop and its yield is determined by successful development of the female inflorescence, the ear. We investigated drought stress responses across several developmental stages of the maize B73 inbred line under field conditions. Drought suppressed plant growth, but had little impact on progression through developmental stages. While ear growth was suppressed by drought, the process of spikelet initiation was not significantly affected. Tassel growth was reduced to a lesser extent compared to the observed reduction in ear growth under stress. Parallel RNA-seq profiling of leaves, ears, and tassels at several developmental stages revealed tissue-specific differences in response to drought stress. High temperature fluctuation was an additional environmental factor that also likely influenced gene expression patterns in the field. Drought induced significant transcriptional changes in leaves and ears but only minor changes in the tassel. Additionally, more genes were drought responsive in ears compared to leaves over the course of drought treatment. Genes that control DNA replication, cell cycle, and cell division were significantly down-regulated in stressed ears, which was consistent with inhibition of ear growth under drought. Inflorescence meristem genes were affected by drought to a lesser degree which was consistent with the minimal impact of drought on spikelet initiation. In contrast, genes that are involved in floret and ovule development were sensitive to stress, which is consistent with the detrimental effect of drought on gynoecium development and kernel set.

CRISPR-Cas9 Editing in Maize: Systematic Evaluation of Off-target Activity and Its Relevance in Crop Improvement.

CRISPR-Cas9 enabled genome engineering has great potential for improving agriculture productivity, but the possibility of unintended off-target edits has evoked some concerns. Here we employ a three-step strategy to investigate Cas9 nuclease specificity in a complex plant genome. Our approach pairs computational prediction with genome-wide biochemical off-target detection followed by validation in maize plants. Our results reveal high frequency (up to 90%) on-target editing with no evidence of off-target cleavage activity when guide RNAs were bioinformatically predicted to be specific. Predictable off-target edits were observed but only with a promiscuous guide RNA intentionally designed to validate our approach. Off-target editing can be minimized by designing guide RNAs that are different from other genomic locations by at least three mismatches in combination with at least one mismatch occurring in the PAM proximal region. With well-designed guides, genetic variation from Cas9 off-target cleavage in plants is negligible, and much less than inherent variation.

Characterization of Proteome Variation During Modern Maize Breeding.

The success of modern maize breeding has been demonstrated by remarkable increases in productivity with tremendous modification of agricultural phenotypes over the last century. Although the underlying genetic changes of the maize adaptation from tropical to temperate regions have been extensively studied, our knowledge is limited regarding the accordance of protein and mRNA expression levels accompanying such adaptation. Here we conducted an integrative analysis of proteomic and transcriptomic changes in a maize association panel. The minimum extent of correlation between protein and RNA levels suggests that variation in mRNA expression is often not indicative of protein expression at a population scale. This is corroborated by the observation that mRNA- and protein-based coexpression networks are relatively independent of each other, and many pQTLs arise without the presence of corresponding eQTLs. Importantly, compared with transcriptome, the subtypes categorized by the proteome show a markedly high accuracy to resemble the genomic subpopulation. These findings suggest that proteome evolved under a greater evolutionary constraint than transcriptome during maize adaptation from tropical to temperate regions. Overall, the integrated multi-omics analysis provides a functional context to interpret gene expression variation during modern maize breeding.

Ectopic expression of ARGOS8 reveals a role for ethylene in root-lodging resistance in maize.

Ethylene plays a critical role in many diverse processes in plant development. Recent studies have demonstrated that overexpression of the maize ARGOS8 gene reduces the plant's response to ethylene by decreasing ethylene signaling and enhances grain yield in transgenic maize plants. The objective of this study was to determine the effects of ethylene on the development of nodal roots, which are primarily responsible for root-lodging resistance in maize. Exogenous application of the ethylene precursor 1-aminocyclopropane-1-carboxylic acid (ACC) was found to promote the emergence of nodal roots. Transcriptome analysis of nodal tissues revealed that the expression of genes involved in metabolic processes and cell wall biogenesis was upregulated in response to ACC treatment, supporting the notion that ethylene is a positive regulator for the outgrowth of young root primordia. In BSV::ARGOS8 transgenic plants with reduced ethylene sensitivity due to constitutive overexpression of ARGOS8, nodal root emergence was delayed and the promotional effect of ACC on nodal root emergence decreased. Field tests showed that the BSV::ARGOS8 plants had higher root lodging relative to non-transgenic controls. When ARGOS8 expression was controlled by the developmentally regulated promoter FTM1, which conferred ARGOS8 overexpression in adult plants but not in the nodal roots and nodes in juvenile plants, the FTM1::ARGOS8 plants had no significant difference in root lodging compared with the wild type but produced a higher grain yield. These results suggest that ethylene has a role in promoting nodal root emergence and that a delay in nodal root development has a negative effect on root-lodging resistance in maize.

Developing a flexible, high-efficiency Agrobacterium-mediated sorghum transformation system with broad application.

Sorghum is the fifth most widely planted cereal crop in the world and is commonly cultivated in arid and semi-arid regions such as Africa. Despite its importance as a food source, sorghum genetic improvement through transgenic approaches has been limited because of an inefficient transformation system. Here, we report a ternary vector (also known as cohabitating vector) system using a recently described pVIR accessory plasmid that facilitates efficient Agrobacterium-mediated transformation of sorghum. We report regeneration frequencies ranging from 6% to 29% in Tx430 using different selectable markers and single copy, backbone free 'quality events' ranging from 45% to 66% of the total events produced. Furthermore, we successfully applied this ternary system to develop transformation protocols for popular but recalcitrant African varieties including Macia, Malisor 84-7 and Tegemeo. In addition, we report the use of this technology to develop the first stable CRISPR/Cas9-mediated gene knockouts in Tx430.

Sequencing, de novo assembly and annotation of a pink bollworm larval midgut transcriptome.

BACKGROUND: The pink bollworm Pectinophora gossypiella (Saunders) (Lepidoptera: Gelechiidae) is one of the world's most important pests of cotton. Insecticide sprays and transgenic cotton producing toxins of the bacterium Bacillus thuringiensis (Bt) are currently used to manage this pest. Bt toxins kill susceptible insects by specifically binding to and destroying midgut cells, but they are not toxic to most other organisms. Pink bollworm is useful as a model for understanding insect responses to Bt toxins, yet advances in understanding at the molecular level have been limited because basic genomic information is lacking for this cosmopolitan pest. Here, we have sequenced, de novo assembled and annotated a comprehensive larval midgut transcriptome from a susceptible strain of pink bollworm. FINDINGS: A de novo transcriptome assembly for the midgut of P. gossypiella was generated containing 46,458 transcripts (average length of 770 bp) derived from 39,874 unigenes. The size of the transcriptome is similar to published midgut transcriptomes of other Lepidoptera and includes up to 91 % annotated contigs. The dataset is publicly available in NCBI and GigaDB as a resource for researchers. CONCLUSIONS: Foundational knowledge of protein-coding genes from the pink bollworm midgut is critical for understanding how this important insect pest functions. The transcriptome data presented here represent the first large-scale molecular resource for this species, and may be used for deciphering relevant midgut proteins critical for xenobiotic detoxification, nutrient digestion and allocation, as well as for the discovery of protein receptors important for Bt intoxication.

Genome-Wide Analysis of Alternative Splicing during Development and Drought Stress in Maize.

Alternative splicing plays a crucial role in plant development as well as stress responses. Although alternative splicing has been studied during development and in response to stress, the interplay between these two factors remains an open question. To assess the effects of drought stress on developmentally regulated splicing in maize (Zea mays), 94 RNA-seq libraries from ear, tassel, and leaf of the B73 public inbred line were constructed at four developmental stages under both well-watered and drought conditions. This analysis was supplemented with a publicly available series of 53 libraries from developing seed, embryo, and endosperm. More than 48,000 novel isoforms, often with stage- or condition-specific expression, were uncovered, suggesting that developmentally regulated alternative splicing occurs in thousands of genes. Drought induced large developmental splicing changes in leaf and ear but relatively few in tassel. Most developmental stage-specific splicing changes affected by drought were tissue dependent, whereas stage-independent changes frequently overlapped between leaf and ear. A linear relationship was found between gene expression changes in splicing factors and alternative spicing of other genes during development. Collectively, these results demonstrate that alternative splicing is strongly associated with tissue type, developmental stage, and stress condition.

Transcriptome analysis of near-isogenic lines provides molecular insights into starch biosynthesis in maize kernel.

Starch is the major component in maize kernels, providing a stable carbohydrate source for humans and livestock as well as raw material for the biofuel industry. Increasing maize kernel starch content will help meet industry demands and has the potential to increase overall yields. We developed a pair of maize near-isogenic lines (NILs) with different alleles for a starch quantitative trait locus on chromosome 3 (qHS3), resulting in different kernel starch content. To investigate the candidate genes for qHS3 and elucidate their effects on starch metabolism, RNA-Seq was performed for the developing kernels of the NILs at 14 and 21 d after pollination (DAP). Analysis of genomic and transcriptomic data identified 76 genes with nonsynonymous single nucleotide polymorphisms and 384 differentially expressed genes (DEGs) in the introgressed fragment, including a hexokinase gene, ZmHXK3a, which catalyzes the conversion of glucose to glucose-6-phosphate and may play a key role in starch metabolism. The expression pattern of all DEGs in starch metabolism shows that altered expression of the candidate genes for qHS3 promoted starch synthesis, with positive consequences for kernel starch content. These results expand the current understanding of starch biosynthesis and accumulation in maize kernels and provide potential candidate genes to increase starch content.

Southern-by-Sequencing: A Robust Screening Approach for Molecular Characterization of Genetically Modified Crops.

Molecular characterization of events is an integral part of the advancement process during genetically modified (GM) crop product development. Assessment of these events is traditionally accomplished by polymerase chain reaction (PCR) and Southern blot analyses. Southern blot analysis can be time-consuming and comparatively expensive and does not provide sequence-level detail. We have developed a sequence-based application, Southern-by-Sequencing (SbS), utilizing sequence capture coupled with next-generation sequencing (NGS) technology to replace Southern blot analysis for event selection in a high-throughput molecular characterization environment. SbS is accomplished by hybridizing indexed and pooled whole-genome DNA libraries from GM plants to biotinylated probes designed to target the sequence of transformation plasmids used to generate events within the pool. This sequence capture process enriches the sequence data obtained for targeted regions of interest (transformation plasmid DNA). Taking advantage of the DNA adjacent to the targeted bases (referred to as next-to-target sequence) that accompanies the targeted transformation plasmid sequence, the data analysis detects plasmid-to-genome and plasmid-to-plasmid junctions introduced during insertion into the plant genome. Analysis of these junction sequences provides sequence-level information as to the following: the number of insertion loci including detection of unlinked, independently segregating, small DNA fragments; copy number; rearrangements, truncations, or deletions of the intended insertion DNA; and the presence of transformation plasmid backbone sequences. This molecular evidence from SbS analysis is used to characterize and select GM plants meeting optimal molecular characterization criteria. SbS technology has proven to be a robust event screening tool for use in a high-throughput molecular characterization environment.

Hierarchal clustering yields insight into multidrug-resistant bacteria isolated from a cattle feedlot wastewater treatment system.

Forty-two percent of Escherichia coli and 58% of Enterococcus spp. isolated from cattle feedlot runoff and associated infiltration basin and constructed wetland treatment system were resistant to at least one antibiotic of clinical importance; a high level of multidrug resistance (22% of E. coli and 37% of Enterococcus spp.) was observed. Hierarchical clustering revealed a closely associated resistance cluster among drug-resistant E. coli isolates that included cephalosporins (ceftiofur, cefoxitin, and ceftriaxone), aminoglycosides (gentamycin, kanamycin, and amikacin), and quinolone nalidixic acid; antibiotics from these classes were used at the study site, and cross-resistance may be associated with transferrable multiple-resistance elements. For Enterococcus spp., co-resistance among vancomycin, linezolid, and daptomycin was common; these antibiotics are reserved for complicated clinical infections and have not been approved for animal use. Vancomycin resistance (n = 49) only occurred when isolates were resistant to linezolid, daptomycin, and all four of the MLSB (macrolide-lincosamide-streptogramin B) antibiotics tested (tylosin, erythromycin, lincomycin, and quinipristin/dalfopristin). This suggests that developing co-resistance to MLSB antibiotics along with cyclic lipopeptides and oxazolidinones may result in resistance to vancomycin as well. Effects of the treatment system on antibiotic resistance were pronounced during periods of no rainfall and low flow (long residence time). Increased hydraulic loading (short residence time) under the influence of rain caused antibiotic-resistant bacteria to be flushed through the treatment system. This presents concern for environmental discharge of multidrug-resistant organisms relevant to public health.

Genome-wide analysis of alternative splicing in Zea mays: landscape and genetic regulation.

Alternative splicing enhances transcriptome diversity in all eukaryotes and plays a role in plant tissue identity and stress adaptation. To catalog new maize (Zea mays) transcripts and identify genomic loci that regulate alternative splicing, we analyzed over 90 RNA-seq libraries from maize inbred lines B73 and Mo17, as well as Syn10 doubled haploid lines (progenies from B73 × Mo17). Transcript discovery was augmented with publicly available data from 14 maize tissues, expanding the maize transcriptome by more than 30,000 and increasing the percentage of intron-containing genes that undergo alternative splicing to 40%. These newly identified transcripts greatly increase the diversity of the maize proteome, sometimes coding for entirely different proteins compared with their most similar annotated isoform. In addition to increasing proteome diversity, many genes encoding novel transcripts gained an additional layer of regulation by microRNAs, often in a tissue-specific manner. We also demonstrate that the majority of genotype-specific alternative splicing can be genetically mapped, with cis-acting quantitative trait loci (QTLs) predominating. A large number of trans-acting QTLs were also apparent, with nearly half located in regions not shown to contain genes associated with splicing. Taken together, these results highlight the currently underappreciated role that alternative splicing plays in tissue identity and genotypic variation in maize.

tassel-less1 encodes a boron channel protein required for inflorescence development in maize.

tassel-less1 (tls1) is a classical maize (Zea mays) inflorescence mutant. Homozygous mutant plants have no tassels or very small tassels, and ear development is also impaired. Using a positional cloning approach, ZmNIP3;1 (a NOD26-like intrinsic protein) was identified as the candidate gene for tls1. The ZmNIP3;1 gene is completely deleted in the tls1 mutant genome. Two Mutator-insertional TUSC alleles of ZmNIP3;1 exhibited tls1-like phenotypes, and allelism tests confirmed that the tls1 gene encodes ZmNIP3;1. Transgenic plants with an RNA interference (RNAi) construct to down-regulate ZmNIP3;1 also showed tls1-like phenotypes, further demonstrating that TLS1 is ZmNIP3;1. Sequence analysis suggests that ZmNIP3;1 is a boron channel protein. Foliar application of boron could rescue the tls1 phenotypes and restore the normal tassel and ear development. Gene expression analysis indicated that in comparison with that of the wild type or tls1 plants treated with boron, the transition from the vegetative to reproductive phase or the development of the floral meristem is impaired in the shoot apical meristem of the tls1 mutant plants. It is concluded that the tls1 mutant phenotypes are caused by impaired boron transport, and boron is essential for inflorescence development in maize.