A multi-omic Nicotiana benthamiana resource for fundamental research and biotechnology.

Nicotiana benthamiana is an invaluable model plant and biotechnology platform with a ~3 Gb allotetraploid genome. To further improve its usefulness and versatility, we have produced high-quality chromosome-level genome assemblies, coupled with transcriptome, epigenome, microRNA and transposable element datasets, for the ubiquitously used LAB strain and a related wild accession, QLD. In addition, single nucleotide polymorphism maps have been produced for a further two laboratory strains and four wild accessions. Despite the loss of five chromosomes from the ancestral tetraploid, expansion of intergenic regions, widespread segmental allopolyploidy, advanced diploidization and evidence of recent bursts of Copia pseudovirus (Copia) mobility not seen in other Nicotiana genomes, the two subgenomes of N. benthamiana show large regions of synteny across the Solanaceae. LAB and QLD have many genetic, metabolic and phenotypic differences, including disparate RNA interference responses, but are highly interfertile and amenable to genome editing and both transient and stable transformation. The LAB/QLD combination has the potential to be as useful as the Columbia-0/Landsberg errecta partnership, utilized from the early pioneering days of Arabidopsis genomics to today.

Spatiotemporal changes in gibberellin content are required for soybean nodulation.

The plant hormone gibberellin (GA) is required at different stages of legume nodule development, with its spatiotemporal distribution tightly regulated. Transcriptomic and bioinformatic analyses established that several key GA biosynthesis and catabolism enzyme encoding genes are critical to soybean (Glycine max) nodule formation. We examined the expression of several GA oxidase genes and used a Förster resonance energy transfer-based GA biosensor to determine the bioactive GA content of roots inoculated with DsRed-labelled Bradyrhizobium diazoefficiens. We manipulated the level of GA by genetically disrupting the expression of GA oxidase genes. Moreover, exogenous treatment of soybean roots with GA3 induced the expression of key nodulation genes and altered infection thread and nodule phenotypes. GmGA20ox1a, GmGA3ox1a, and GmGA2ox1a are upregulated in soybean roots inoculated with compatible B. diazoefficiens. GmGA20ox1a expression is predominately localized to the transient meristem of soybean nodules and coincides with the spatiotemporal distribution of bioactive GA occurring throughout nodule organogenesis. GmGA2ox1a exhibits a nodule vasculature-specific expression pattern, whereas GmGA3ox1a can be detected throughout the nodule and root. Disruptions in the level of GA resulted in aberrant rhizobia infection and reduced nodule numbers. Collectively, our results establish a central role for GAs in root hair infection by symbiotic rhizobia and in nodule organogenesis.

Homo sapiens: The Superspreader of Plant Viral Diseases.

Plant viruses are commonly vectored by flying or crawling animals, such as aphids and beetles, and cause serious losses in major agricultural and horticultural crops. Controlling virus spread is often achieved by minimizing a crop's exposure to the vector, or by reducing vector numbers with compounds such as insecticides. A major, but less obvious, factor not controlled by these measures is Homo sapiens. Here, we discuss the inconvenient truth of how humans have become superspreaders of plant viruses on both a local and a global scale.

Are the current gRNA ranking prediction algorithms useful for genome editing in plants?

Introducing a new trait into a crop through conventional breeding commonly takes decades, but recently developed genome sequence modification technology has the potential to accelerate this process. One of these new breeding technologies relies on an RNA-directed DNA nuclease (CRISPR/Cas9) to cut the genomic DNA, in vivo, to facilitate the deletion or insertion of sequences. This sequence specific targeting is determined by guide RNAs (gRNAs). However, choosing an optimum gRNA sequence has its challenges. Almost all current gRNA design tools for use in plants are based on data from experiments in animals, although many allow the use of plant genomes to identify potential off-target sites. Here, we examine the predictive uniformity and performance of eight different online gRNA-site tools. Unfortunately, there was little consensus among the rankings by the different algorithms, nor a statistically significant correlation between rankings and in vivo effectiveness. This suggests that important factors affecting gRNA performance and/or target site accessibility, in plants, are yet to be elucidated and incorporated into gRNA-site prediction tools.

Wild Origins of Macadamia Domestication Identified Through Intraspecific Chloroplast Genome Sequencing.

Identifying the geographic origins of crops is important for the conservation and utilization of novel genetic variation. Even so, the origins of many food crops remain elusive. The tree nut crop macadamia has a remarkable domestication history, from subtropical rain forests in Australia through Hawaii to global cultivation all within the last century. The industry is based primarily on Macadamia integrifolia and M. integrifolia-M. tetraphylla hybrid cultivars with Hawaiian cultivars the main contributors to world production. Sequence data from the chloroplast genome assembled using a genome skimming strategy was used to determine population structure among remnant populations of the main progenitor species, M. integrifolia. Phylogenetic analysis of a 506 bp chloroplast SNP alignment from 64 wild and cultivated accessions identified phylogeographic structure and deep divergences between clades providing evidence for historical barriers to seed dispersal. High levels of variation were detected among wild accessions. Most Hawaiian cultivars, however, shared a single chlorotype that was also present at two wild sites at Mooloo and Mt Bauple from the northernmost distribution of the species in south-east Queensland. Our results provide evidence for a maternal genetic bottleneck during early macadamia domestication, and pinpoint the likely source of seed used to develop the Hawaiian cultivars. The extensive variability and structuring of M. integrifolia chloroplast genomic variation detected in this study suggests much unexploited genetic diversity is available for improvement of this recently domesticated crop.

Integrated physical map of bread wheat chromosome arm 7DS to facilitate gene cloning and comparative studies.

Bread wheat (Triticum aestivum L.) is a staple food for a significant part of the world's population. The growing demand on its production can be satisfied by improving yield and resistance to biotic and abiotic stress. Knowledge of the genome sequence would aid in discovering genes and QTLs underlying these traits and provide a basis for genomics-assisted breeding. Physical maps and BAC clones associated with them have been valuable resources from which to generate a reference genome of bread wheat and to assist map-based gene cloning. As a part of a joint effort coordinated by the International Wheat Genome Sequencing Consortium, we have constructed a BAC-based physical map of bread wheat chromosome arm 7DS consisting of 895 contigs and covering 94% of its estimated length. By anchoring BAC contigs to one radiation hybrid map and three high resolution genetic maps, we assigned 73% of the assembly to a distinct genomic position. This map integration, interconnecting a total of 1713 markers with ordered and sequenced BAC clones from a minimal tiling path, provides a tool to speed up gene cloning in wheat. The process of physical map assembly included the integration of the 7DS physical map with a whole-genome physical map of Aegilops tauschii and a 7DS Bionano genome map, which together enabled efficient scaffolding of physical-map contigs, even in the non-recombining region of the genetic centromere. Moreover, this approach facilitated a comparison of bread wheat and its ancestor at BAC-contig level and revealed a reconstructed region in the 7DS pericentromere.

Na+/Ca2+ exchanger contributes to stool transport in mice with experimental diarrhea.

The Na+/Ca2+ exchanger (NCX) is a bidirectional transporter that is controlled by membrane potential and transmembrane gradients of Na+ and Ca2+. To reveal the functional role of NCX on gastrointestinal motility, we have previously used NCX1 and NCX2 heterozygote knockout mice (HET). We found that NCX1 and NCX2 play important roles in the motility of the gastric fundus, ileum and distal colon. Therefore, we believed that NCX1 and NCX2 play an important role in transport of intestinal contents. Here, we investigated the role of NCX in a mouse model of drug-induced diarrhea. The fecal consistencies in NCX1 HET and NCX2 HET were assessed using a diarrhea induced by magnesium sulfate, 5-hydroxytryptamine (5-HT) and prostaglandin E2 (PGE2). NCX2 HET, but not NCX1 HET, exacerbated magnesium sulfate-induced diarrhea by increasing watery fecals. Likewise, 5-HT-induced diarrheas were exacerbated in NCX2 HET, but not NCX1 HET. However, NCX1 HET and NCX2 HET demonstrated PGE2 induced diarrhea similar to those of wild-type mice (WT). As well as the result of the distal colon shown previously, in the proximal and transverse colons of WT, the myenteric plexus layers and the longitudinal and circular muscle layers were strongly immunoreactive to NCX1 and NCX2. In this study, we demonstrate that NCX2 has important roles in development of diarrhea.

Overexpression of Na+/Ca2+ exchanger 1 display enhanced relaxation in the gastric fundus.

In gastric smooth muscles, the released Ca2+ activates the contractile proteins and Ca2+ taken up from the cytosol cause relaxation. The Na+/Ca2+ exchanger (NCX) is an antiporter membrane protein that controls Ca2+ influx and efflux across the membrane. However, the possible relation of NCX in gastric fundus motility is largely unknown. Here, we investigated electric field stimulation (EFS)-induced relaxations in the circular muscles of the gastric fundus in smooth muscle-specific NCX1 transgenic mice (Tg). EFS caused a bi-phasic response, transient and sustained relaxation. The sustained relaxation prolonged for an extended period after the end of the stimulus. EFS-induced transient relaxation and sustained relaxation were greater in Tg than in wild-type mice (WT). Disruption of nitric oxide component by N-nitro-l-arginine, EFS-induced transient and sustained relaxations caused still marked in Tg compared to WT. Inhibition of PACAP by antagonist, EFS-induced sustained relaxation in Tg was not seen, similar to WT. Nevertheless, transient relaxation remained more pronounced in Tg than in WT. Next, we examined responses to NO and PACAP in smooth muscles. The magnitudes of NOR-1, which generates NO, and PACAP-induced relaxations were greater in Tg than in WT. In this study, we demonstrate that NCX1 regulates gastric fundus motility.

Roles of Na(+)/Ca(2+) exchanger isoforms NCX1 and NCX2 in motility in mouse ileum.

The Na(+)/Ca(2+) exchanger (NCX) is a plasma membrane transporter that is involved in regulating intracellular Ca(2+) concentrations in various tissues. The physiological roles by which NCX influences gastrointestinal motility are incompletely understood, although its role in the heart, brain, and kidney has been widely investigated. In this study, we focused on the functions of the NCX isoforms, NCX1 and NCX2, in the motility of the ileum in the gastrointestinal tract. We investigated the response to electric field stimulation (EFS) in the longitudinal smooth muscle of the ileum obtained from wild-type mice (WT), NCX1-heterozygote knockout mice (NCX1 HET), NCX2 HET and smooth muscle-specific NCX1.3 transgenic mice (NCX1.3 Tg). EFS induced a phasic contraction that persisted during EFS and a tonic contraction that occurred after the end of EFS. We found that the amplitudes of the phasic and tonic contractions were significantly smaller in NCX2 HET, but not in NCX1 HET, compared to WT. Moreover, the magnitudes of acetylcholine (ACh)- and substance P (SP)-induced contractions of NCX2 HET, but not of NCX1 HET, were smaller compared to WT. In contrast, the amplitudes of the phasic and tonic contractions were greater in NCX1.3 Tg compared to WT. Similar to EFS, the magnitude of ACh-induced contraction was greater in NCX1.3 Tg than in WT. Taken together, our findings indicated that NCX1 and NCX2 play important roles in ileal motility and suggest that NCX1 and NCX2 regulate the motility in the ileum by controlling the sensitivity of smooth muscles to ACh and SP.

An efficient approach to BAC based assembly of complex genomes.

BACKGROUND: There has been an exponential growth in the number of genome sequencing projects since the introduction of next generation DNA sequencing technologies. Genome projects have increasingly involved assembly of whole genome data which produces inferior assemblies compared to traditional Sanger sequencing of genomic fragments cloned into bacterial artificial chromosomes (BACs). While whole genome shotgun sequencing using next generation sequencing (NGS) is relatively fast and inexpensive, this method is extremely challenging for highly complex genomes, where polyploidy or high repeat content confounds accurate assembly, or where a highly accurate 'gold' reference is required. Several attempts have been made to improve genome sequencing approaches by incorporating NGS methods, to variable success. RESULTS: We present the application of a novel BAC sequencing approach which combines indexed pools of BACs, Illumina paired read sequencing, a sequence assembler specifically designed for complex BAC assembly, and a custom bioinformatics pipeline. We demonstrate this method by sequencing and assembling BAC cloned fragments from bread wheat and sugarcane genomes. CONCLUSIONS: We demonstrate that our assembly approach is accurate, robust, cost effective and scalable, with applications for complete genome sequencing in large and complex genomes.

BioNano genome mapping of individual chromosomes supports physical mapping and sequence assembly in complex plant genomes.

The assembly of a reference genome sequence of bread wheat is challenging due to its specific features such as the genome size of 17 Gbp, polyploid nature and prevalence of repetitive sequences. BAC-by-BAC sequencing based on chromosomal physical maps, adopted by the International Wheat Genome Sequencing Consortium as the key strategy, reduces problems caused by the genome complexity and polyploidy, but the repeat content still hampers the sequence assembly. Availability of a high-resolution genomic map to guide sequence scaffolding and validate physical map and sequence assemblies would be highly beneficial to obtaining an accurate and complete genome sequence. Here, we chose the short arm of chromosome 7D (7DS) as a model to demonstrate for the first time that it is possible to couple chromosome flow sorting with genome mapping in nanochannel arrays and create a de novo genome map of a wheat chromosome. We constructed a high-resolution chromosome map composed of 371 contigs with an N50 of 1.3 Mb. Long DNA molecules achieved by our approach facilitated chromosome-scale analysis of repetitive sequences and revealed a ~800-kb array of tandem repeats intractable to current DNA sequencing technologies. Anchoring 7DS sequence assemblies obtained by clone-by-clone sequencing to the 7DS genome map provided a valuable tool to improve the BAC-contig physical map and validate sequence assembly on a chromosome-arm scale. Our results indicate that creating genome maps for the whole wheat genome in a chromosome-by-chromosome manner is feasible and that they will be an affordable tool to support the production of improved pseudomolecules.

Type B Heterotrimeric G Protein γ-Subunit Regulates Auxin and ABA Signaling in Tomato.

Heterotrimeric G proteins composed of α, ß, and γ subunits are central signal transducers mediating the cellular response to multiple stimuli in most eukaryotes. Gγ subunits provide proper cellular localization and functional specificity to the heterotrimer complex. Plant Gγ subunits, divided into three structurally distinct types, are more diverse than their animal counterparts. Type B Gγ subunits, lacking a carboxyl-terminal isoprenylation motif, are found only in flowering plants. We present the functional characterization of type B Gγ subunit (SlGGB1) in tomato (Solanum lycopersicum). We show that SlGGB1 is the most abundant Gγ subunit in tomato and strongly interacts with the Gß subunit. Importantly, the green fluorescent protein-SlGGB1 fusion protein as well as the carboxyl-terminal yellow fluorescent protein-SlGGB1/amino-terminal yellow fluorescent protein-Gß heterodimer were localized in the plasma membrane, nucleus, and cytoplasm. RNA interference-mediated silencing of SlGGB1 resulted in smaller seeds, higher number of lateral roots, and pointy fruits. The silenced lines were hypersensitive to exogenous auxin, while levels of endogenous auxins were lower or similar to those of the wild type. SlGGB1-silenced plants also showed strong hyposensitivity to abscisic acid (ABA) during seed germination but not in other related assays. Transcriptome analysis of the transgenic seeds revealed abnormal expression of genes involved in ABA sensing, signaling, and response. We conclude that the type B Gγ subunit SlGGB1 mediates auxin and ABA signaling in tomato.

Germin-like protein 2 gene promoter from rice is responsive to fungal pathogens in transgenic potato plants.

Controlled transgene expression via a promoter is particularly triggered in response to pathogen infiltration. This is significant for eliciting disease-resistant features in crops through genetic engineering. The germins and germin-like proteins (GLPs) are known to be associated with plant and developmental stages. The 1107-bp Oryza sativa root GLP2 (OsRGLP2) gene promoter fused to a ß-glucuronidase (GUS) reporter gene was transformed into potato plants through an Agrobacterium-mediated transformation. The OsRGLP2 promoter was activated in response to Fusarium solani (Mart.) Sacc. and Alternaria solani Sorauer. Quantitative real-time PCR results revealed 4-5-fold increase in promoter activity every 24 h following infection. There was a 15-fold increase in OsRGLP2 promoter activity after 72 h of F. solani (Mart.) Sacc. treatment and a 12-fold increase observed with A. solani Sorauer. Our results confirmed that the OsRGLP2 promoter activity was enhanced under fungal stress. Furthermore, a hyperaccumulation of H2O2 in transgenic plants is a clear signal for the involvement of OsRGLP2 promoter region in the activation of specific genes in the potato genome involved in H2O2-mediated defense response. The OsRGLP2 promoter evidently harbors copies of GT-I and Dof transcription factors (AAAG) that act in response to elicitors generated in the wake of pathogen infection.

Generation of the antioxidant yellow pigment derived from 4-methylthio-3-butenyl isothiocyanate in salted radish roots (takuan-zuke).

2-[3-(2-Thioxopyrrolidin-3-ylidene)methyl]-tryptophan (TPMT) is a yellow pigment of salted radish roots (takuan-zuke) derived from 4-methylthio-3-butenyl isothiocyanate (MTBITC), the pungent component of radish roots. Here, we prepared salted radish and analyzed the behavior of the yellow pigment and related substances in the dehydration process and long-term salting process. All salted radish samples turned yellow, and their b(*) values increased with time and temperature. The salted radish that was sun-dried and pickled at room temperature turned the brightest yellow, and the generation of TPMT was clearly confirmed. These results indicate that tissue shrinkage due to dehydration, salting temperature, and pH play important roles in the yellowing of takuan-zuke.

High-resolution skim genotyping by sequencing reveals the distribution of crossovers and gene conversions in Cicer arietinum and Brassica napus.

KEY MESSAGE: We characterise the distribution of crossover and non-crossover recombination in Brassica napus and Cicer arietinum using a low-coverage genotyping by sequencing pipeline SkimGBS. The growth of next-generation DNA sequencing technologies has led to a rapid increase in sequence-based genotyping for applications including diversity assessment, genome structure validation and gene-trait association. We have established a skim-based genotyping by sequencing method for crop plants and applied this approach to genotype-segregating populations of Brassica napus and Cicer arietinum. Comparison of progeny genotypes with those of the parental individuals allowed the identification of crossover and non-crossover (gene conversion) events. Our results identify the positions of recombination events with high resolution, permitting the mapping and frequency assessment of recombination in segregating populations.

Functional analysis of duplicated Symbiosis Receptor Kinase (SymRK) genes during nodulation and mycorrhizal infection in soybean (Glycine max).

Association between legumes and rhizobia results in the formation of root nodules, where symbiotic nitrogen fixation occurs. The early stages of this association involve a complex of signalling events between the host and microsymbiont. Several genes dealing with early signal transduction have been cloned, and one of them encodes the leucine-rich repeat (LRR) receptor kinase (SymRK; also termed NORK). The Symbiosis Receptor Kinase gene is required by legumes to establish a root endosymbiosis with Rhizobium bacteria as well as mycorrhizal fungi. Using degenerate primer and BAC sequencing, we cloned duplicated SymRK homeologues in soybean called GmSymRKα and GmSymRKß. These duplicated genes have high similarity of nucleotide (96%) and amino acid sequence (95%). Sequence analysis predicted a malectin-like domain within the extracellular domain of both genes. Several putative cis-acting elements were found in promoter regions of GmSymRKα and GmSymRKß, suggesting a participation in lateral root development, cell division and peribacteroid membrane formation. The mutant of SymRK genes is not available in soybean; therefore, to know the functions of these genes, RNA interference (RNAi) of these duplicated genes was performed. For this purpose, RNAi construct of each gene was generated and introduced into the soybean genome by Agrobacterium rhizogenes-mediated hairy root transformation. RNAi of GmSymRKß gene resulted in an increased reduction of nodulation and mycorrhizal infection than RNAi of GmSymRKα, suggesting it has the major activity of the duplicated gene pair. The results from the important crop legume soybean confirm the joint phenotypic action of GmSymRK genes in both mycorrhizal and rhizobial infection seen in model legumes.

The value of biodiversity in legume symbiotic nitrogen fixation and nodulation for biofuel and food production.

Much of modern agriculture is based on immense populations of genetically identical or near-identical varieties, called cultivars. However, advancement of knowledge, and thus experimental utility, is found through biodiversity, whether naturally-found or induced by the experimenter. Globally we are confronted by ever-growing food and energy challenges. Here we demonstrate how such biodiversity from the food legume crop soybean (Glycine max L. Merr) and the bioenergy legume tree Pongamia (Millettia) pinnata is a great value. Legume plants are diverse and are represented by over 18,000 species on this planet. Some, such as soybean, pea and medics are used as food and animal feed crops. Others serve as ornamental (e.g., wisteria), timber (e.g., acacia/wattle) or biofuel (e.g., Pongamia pinnata) resources. Most legumes develop root organs (nodules) after microsymbiont induction that serve as their habitat for biological nitrogen fixation. Through this, nitrogen fertiliser demand is reduced by the efficient symbiosis between soil Rhizobium-type bacteria and the appropriate legume partner. Mechanistic research into the genetics, biochemistry and physiology of legumes is thus strategically essential for future global agriculture. Here we demonstrate how molecular plant science analysis of the genetics of an established food crop (soybean) and an emerging biofuel P. pinnata feedstock contributes to their utility by sustainable production aided by symbiotic nitrogen fixation.

A chromosomal genomics approach to assess and validate the desi and kabuli draft chickpea genome assemblies.

With the expansion of next-generation sequencing technology and advanced bioinformatics, there has been a rapid growth of genome sequencing projects. However, while this technology enables the rapid and cost-effective assembly of draft genomes, the quality of these assemblies usually falls short of gold standard genome assemblies produced using the more traditional BAC by BAC and Sanger sequencing approaches. Assembly validation is often performed by the physical anchoring of genetically mapped markers, but this is prone to errors and the resolution is usually low, especially towards centromeric regions where recombination is limited. New approaches are required to validate reference genome assemblies. The ability to isolate individual chromosomes combined with next-generation sequencing permits the validation of genome assemblies at the chromosome level. We demonstrate this approach by the assessment of the recently published chickpea kabuli and desi genomes. While previous genetic analysis suggests that these genomes should be very similar, a comparison of their chromosome sizes and published assemblies highlights significant differences. Our chromosomal genomics analysis highlights short defined regions that appear to have been misassembled in the kabuli genome and identifies large-scale misassembly in the draft desi genome. The integration of chromosomal genomics tools within genome sequencing projects has the potential to significantly improve the construction and validation of genome assemblies. The approach could be applied both for new genome assemblies as well as published assemblies, and complements currently applied genome assembly strategies.

Mechanistic action of gibberellins in legume nodulation.

Legume plants are capable of entering into a symbiotic relationship with rhizobia bacteria. This results in the formation of novel organs on their roots, called nodules, in which the bacteria capture atmospheric nitrogen and provide it as ammonium to the host plant. Complex molecular and physiological changes are involved in the formation and establishment of such nodules. Several phytohormones are known to play key roles in this process. Gibberellins (gibberellic acids; GAs), a class of phytohormones known to be involved in a wide range of biological processes (i.e., cell elongation, germination) are reported to be involved in the formation and maturation of legume nodules, highlighted by recent transcriptional analyses of early soybean symbiotic steps. Here, we summarize what is currently known about GAs in legume nodulation and propose a model of GA action during nodule development. Results from a wide range of studies, including GA application, mutant phenotyping, and gene expression studies, indicate that GAs are required at different stages, with an optimum, tightly regulated level being key to achieve successful nodulation. Gibberellic acids appear to be required at two distinct stages of nodulation: (i) early stages of rhizobia infection and nodule primordium establishment; and (ii) later stages of nodule maturation.