Search

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

A photoredox protocol that uses a heteroleptic Cu (I) complex, [Cu­(dq)­(BINAP)]­BF4, has been developed for the photodeprotection of benzenesulfonyl-protected N-heterocycles. A range of substrates was examined, including indazoles, indoles, pyrazoles, and benzimidazole, featuring both electron-rich and electron-deficient substituents, giving good yields of the N-heterocycle products with broad functional group tolerance. This transformation was also found to be amenable to flow reaction conditions.

Site-directed nucleases (SDNs) used for targeted genome editing are powerful new tools to introduce precise genetic changes into plants. Like traditional approaches, such as conventional crossing and induced mutagenesis, genome editing aims to improve crop yield and nutrition. Next-generation sequencing studies demonstrate that across their genomes, populations of crop species typically carry millions of single nucleotide polymorphisms and many copy number and structural variants. Spontaneous mutations occur at rates of ∼10(−8) to 10(−9) per site per generation, while variation induced by chemical treatment or ionizing radiation results in higher mutation rates. In the context of SDNs, an off-target change or edit is an unintended, nonspecific mutation occurring at a site with sequence similarity to the targeted edit region. SDN-mediated off-target changes can contribute to a small number of additional genetic variants compared to those that occur naturally in breeding populations or are introduced by induced-mutagenesis methods. Recent studies show that using computational algorithms to design genome editing reagents can mitigate off-target edits in plants. Finally, crops are subject to strong selection to eliminate off-type plants through well-established multigenerational breeding, selection, and commercial variety development practices. Within this context, off-target edits in crops present no new safety concerns compared to other breeding practices. The current generation of genome editing technologies is already proving useful to develop new plant varieties with consumer and farmer benefits. Genome editing will likely undergo improved editing specificity along with new developments in SDN delivery and increasing genomic characterization, further improving reagent design and application.

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 transcriptional activity. Chromatin loops, made of two “anchors” flanking a loop “interior”, overlap with up to 90% of high-resolution interaction domains from a previous public maize interactome dataset. A majority of all anchors are shared between multiple loops, 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 analysis indicate that all anchors studied here co-locate with at least 81.5% of expressed genes and 74% of open chromatin regions. Up to 63% of all unique variants derived from a prior public maize eQTL datasets overlap with Hi-C loop anchors. Anchor annotation suggests that Conclusions A majority of expressed genes and open chromatin regions co-locate with a conserved set of chromatin loop anchors. The results presented here will be a useful reference to further investigate the function of chromatin loop anchors and of the formation of interaction regions in the regulation of gene expression in maize.

A photoredox protocol that uses a heteroleptic Cu (I) complex, [Cu(dq)(BINAP)]BF

C-H functionalization of electron-deficient heteroarenes using commercial unactivated alkyl halides through reductive quenching photoredox catalysis was developed. Mainstream approaches rely on the use of an excess of strong acids that result in regioselectivities dictated by the innate effect of the protonated heteroarene, leaving the functionalization of other carbons unexplored. We report a mild method under basic conditions that allows access to previously underexplored regioselectivities by relying on a combination of conjugate and halogen ortho-directing effects. Overall, this methodology gives quick access to a variety of alkylated heteroarenes that will be of interest to medicinal chemistry programs.

The photoredox cross-coupling of aryl halides and potassium alkyl trifluoroborates is a very effective means to form Csp3–Csp2 bonds. However, this transformation is inefficient for the coupling of unactivated primary trifluoroborates. We have developed a generally useful, continuous flow Csp3–Csp2 coupling procedure for the synthesis of diverse product sets that is compatible with both trifluoroborates and silicate reagents. This universal protocol provides diversity sets from both primary and secondary coupling partners. This easily scalable procedure widens the substrate scope of the coupling reaction and is efficient for producing a greater range of analogues bearing a high sp3 fraction.

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.

Long-read sequencing technologies have greatly facilitated assemblies of large eukaryotic genomes. In this paper, Oxford Nanopore sequences generated on a MinION sequencer are combined with Bionano Genomics Direct Label and Stain (DLS) optical maps to generate a chromosome-scale de novo assembly of the repeat-rich Sorghum bicolor Tx430 genome. The final assembly consists of 29 scaffolds, encompassing in most cases entire chromosome arms. It has a scaffold N50 of 33.28 Mbps and covers 90% of the expected genome length. A sequence accuracy of 99.85% is obtained after aligning the assembly against Illumina Tx430 data and 99.6% of the 34,211 public gene models align to the assembly. Comparisons of Tx430 and BTx623 DLS maps against the public BTx623 v3.0.1 genome assembly suggest substantial discrepancies whose origin remains to be determined. In summary, this study demonstrates that informative assemblies of complex plant genomes can be generated by combining nanopore sequencing with DLS optical maps., Assembly of large, repeat-rich eukaryotic genomes remains challenging. Here, the authors use BioNano Genomics DLS optical mapping and single-molecule nanopore sequencing to generate a chromosome-scale assembly of a new Sorghum bicolor accession and identify variation compared to the publicly available S. bicolor genome.

The advent of long-read sequencing technologies has greatly facilitated assemblies of large eukaryotic genomes. In this paper, Oxford Nanopore sequences generated on a MinION sequencer were combined with BioNano Genomics Direct Label and Stain (DLS) optical maps to generate a chromosome-scale de novo assembly of the repeat-rich Sorghum bicolor Tx430 genome. The final hybrid assembly consists of 29 scaffolds, encompassing in most cases entire chromosome arms. It has a scaffold N50 value of 33.28Mbps and covers >90% of Sorghum bicolor expected genome length. A sequence accuracy of 99.67% was obtained in unique regions after aligning contigs against Illumina Tx430 data. Alignments showed that 99.4% of the 34,211 public gene models are present in the assembly, including 94.2% mapping end-to-end. Comparisons of the DLS optical maps against the public Sorghum Bicolor v3.0.1 BTx623 genome assembly suggest the presence of substantial genomic rearrangements whose origin remains to be determined.

The photoredox cross-coupling of aryl halides and potassium alkyl trifluoroborates is a very effective means to form Csp

Background Soybean (Glycine max) production is significantly hampered by frequent droughts in many regions of the world including the United States. Identifying microRNA (miRNA)-controlled posttranscriptional gene regulation under drought will enhance our understanding of molecular basis of drought tolerance in this important cash crop. Indeed, miRNA profiles in soybean exposed to drought were studied but not from the primary root tips, which is not only a main zone of water uptake but also critical for water stress sensing and signaling. Methods Here we report miRNA profiles specifically from well-watered and water-stressed primary root tips (0 to 8 mm from the root apex) of soybean. Small RNA sequencing confirmed the expression of vastly diverse miRNA (303 individual miRNAs) population, and, importantly several conserved miRNAs were abundantly expressed in primary root tips. Results Notably, 12 highly conserved miRNA families were differentially regulated in response to water-deficit; six were upregulated while six others were downregulated at least by one fold (log2) change. Differentially regulated soybean miRNAs are targeting genes include auxin response factors, Cu/Zn Superoxide dismutases, laccases and plantacyanin and several others. Conclusions These results highlighted the importance of miRNAs in primary root tips both under control and water-deficit conditions; under control conditions, miRNAs could be important for cell division, cell elongation and maintenance of the root apical meristem activity including quiescent centre whereas under water stress differentially regulated miRNAs could decrease auxin signaling and oxidative stress as well as other metabolic processes that save energy and water. Electronic supplementary material The online version of this article (doi:10.1186/s12918-016-0374-0) contains supplementary material, which is available to authorized users.

Genome-scale data offer the opportunity to clarify phylogenetic relationships that are difficult to resolve with few loci, but they can also identify genomic regions with evolutionary history distinct from that of the species history. We collected whole-genome sequence data from 29 taxa in the legume genus Medicago, then aligned these sequences to the Medicago truncatula reference genome to confidently identify 87 596 variable homologous sites. We used this data set to estimate phylogenetic relationships among Medicago species, to investigate the number of sites needed to provide robust phylogenetic estimates and to identify specific genomic regions supporting topologies in conflict with the genome-wide phylogeny. Our full genomic data set resolves relationships within the genus that were previously intractable. Subsampling the data reveals considerable variation in phylogenetic signal and power in smaller subsets of the data. Even when sampling 5000 sites, no random sample of the data supports a topology identical to that of the genome-wide phylogeny. Phylogenetic relationships estimated from 500-site sliding windows revealed genome regions supporting several alternative species relationships among recently diverged taxa, consistent with the expected effects of deep coalescence or introgression in the recent history of Medicago .( Medicago; phylogenomics; whole-genome resequencing.)

The MinION is a portable single-molecule DNA sequencing instrument that was released by Oxford Nanopore Technologies in 2014, producing long sequencing reads by measuring changes in ionic flow when single-stranded DNA molecules translocate through the pores. While MinION long reads have an error rate substantially higher than the ones produced by short-read sequencing technologies, they can generate de novo assemblies of microbial genomes, after an initial correction step that includes alignment of Illumina sequencing data or detection of overlaps between Oxford Nanopore reads to improve accuracy. In this study, MinION reads were generated from the multi-chromosome genome of Agrobacterium tumefaciens strain LBA4404. Errors in the consensus two-directional (sense and antisense) “2D” sequences were first characterized by way of comparison with an internal reference assembly. Both Illumina-based correction and self-correction were performed and the resulting corrected reads assembled into high-quality hybrid and non-hybrid assemblies. Corrected read datasets and assemblies were subsequently compared. The results shown here indicate that both hybrid and non-hybrid methods can be used to assemble Oxford Nanopore reads into informative multi-chromosome assemblies, each with slightly different outcomes in terms of contiguity and accuracy.

Single-nucleotide polymorphisms (SNPs,2000) were discovered by using RNA-seq and allele-specific sequencing approaches in pigeonpea (Cajanus cajan). For making the SNP genotyping cost-effective, successful competitive allele-specific polymerase chain reaction (KASPar) assays were developed for 1616 SNPs and referred to as PKAMs (pigeonpea KASPar assay markers). Screening of PKAMs on 24 genotypes [23 from cultivated species and 1 wild species (Cajanus scarabaeoides)] defined a set of 1154 polymorphic markers (77.4%) with a polymorphism information content (PIC) value from 0.04 to 0.38. One thousand and ninety-four PKAMs showed polymorphisms between parental lines of the reference mapping population (C. cajan ICP 28 × C. scarabaeoides ICPW 94). By using high-quality marker genotyping data on 167 F(2) lines from the population, a comprehensive genetic map comprising 875 PKAMs with an average inter-marker distance of 1.11 cM was developed. Previously mapped 35 simple sequence repeat markers were integrated into the PKAM map and an integrated genetic map of 996.21 cM was constructed. Mapped PKAMs showed a higher degree of synteny with the genome of Glycine max followed by Medicago truncatula and Lotus japonicus and least with Vigna unguiculata. These PKAMs will be useful for genetics research and breeding applications in pigeonpea and for utilizing genome information from other legume species.

The advent of next-generation DNA sequencing (NGS) technologies has led to the development of rapid genome-wide Single Nucleotide Polymorphism (SNP) detection applications in various plant species. Recent improvements in sequencing throughput combined with an overall decrease in costs per gigabase of sequence is allowing NGS to be applied to not only the evaluation of small subsets of parental inbred lines, but also the mapping and characterization of traits of interest in much larger populations. Such an approach, where sequences are used simultaneously to detect and score SNPs, therefore bypassing the entire marker assay development stage, is known as genotyping-by-sequencing (GBS). This review will summarize the current state of GBS in plants and the promises it holds as a genome-wide genotyping application.

A comprehensive transcriptome assembly for pigeonpea has been developed by analyzing 128.9 million short Illumina GA IIx single end reads, 2.19 million single end FLX/454 reads, and 18 353 Sanger expressed sequenced tags from more than 16 genotypes. The resultant transcriptome assembly, referred to as CcTA v2, comprised 21 434 transcript assembly contigs (TACs) with an N50 of 1510 bp, the largest one being ∼8 kb. Of the 21 434 TACs, 16 622 (77.5%) could be mapped on to the soybean genome build 1.0.9 under fairly stringent alignment parameters. Based on knowledge of intron junctions, 10 009 primer pairs were designed from 5033 TACs for amplifying intron spanning regions (ISRs). By using in silico mapping of BAC-end-derived SSR loci of pigeonpea on the soybean genome as a reference, putative mapping positions at the chromosome level were predicted for 6284 ISR markers, covering all 11 pigeonpea chromosomes. A subset of 128 ISR markers were analyzed on a set of eight genotypes. While 116 markers were validated, 70 markers showed one to three alleles, with an average of 0.16 polymorphism information content (PIC) value. In summary, the CcTA v2 transcript assembly and ISR markers will serve as a useful resource to accelerate genetic research and breeding applications in pigeonpea.

• Premise of the study: Next-generation sequencing (NGS) technologies are frequently used for resequencing and mining of single nucleotide polymorphisms (SNPs) by comparison to a reference genome. In crop species such as chickpea (Cicer arietinum) that lack a reference genome sequence, NGS-based SNP discovery is a challenge. Therefore, unlike probability-based statistical approaches for consensus calling and by comparison with a reference sequence, a coverage-based consensus calling (CbCC) approach was applied and two genotypes were compared for SNP identification. • Methods: A CbCC approach is used in this study with four commonly used short read alignment tools (Maq, Bowtie, Novoalign, and SOAP2) and 15.7 and 22.1 million Illumina reads for chickpea genotypes ICC4958 and ICC1882, together with the chickpea trancriptome assembly (CaTA). • Key results: A nonredundant set of 4543 SNPs was identified between two chickpea genotypes. Experimental validation of 224 randomly selected SNPs showed superiority of Maq among individual tools, as 50.0% of SNPs predicted by Maq were true SNPs. For combinations of two tools, greatest accuracy (55.7%) was reported for Maq and Bowtie, with a combination of Bowtie, Maq, and Novoalign identifying 61.5% true SNPs. SNP prediction accuracy generally increased with increasing reads depth. • Conclusions: This study provides a benchmark comparison of tools as well as read depths for four commonly used tools for NGS SNP discovery in a crop species without a reference genome sequence. In addition, a large number of SNPs have been identified in chickpea that would be useful for molecular breeding.

Pigeonpea is an important legume food crop grown primarily by smallholder farmers in many semi-arid tropical regions of the world. We used the Illumina next-generation sequencing platform to generate 237.2 Gb of sequence, which along with Sanger-based bacterial artificial chromosome end sequences and a genetic map, we assembled into scaffolds representing 72.7% (605.78 Mb) of the 833.07 Mb pigeonpea genome. Genome analysis predicted 48,680 genes for pigeonpea and also showed the potential role that certain gene families, for example, drought tolerance-related genes, have played throughout the domestication of pigeonpea and the evolution of its ancestors. Although we found a few segmental duplication events, we did not observe the recent genome-wide duplication events observed in soybean. This reference genome sequence will facilitate the identification of the genetic basis of agronomically important traits, and accelerate the development of improved pigeonpea varieties that could improve food security in many developing countries.

Chantier qualité GA; International audience; Legumes (Fabaceae or Leguminosae) are unique among cultivated plants for their ability to carry out endosymbiotic nitrogen fixation with rhizobial bacteria, a process that takes place in a specialized structure known as the nodule. Legumes belong to one of the two main groups of eurosids, the Fabidae, which includes most species capable of endosymbiotic nitrogen fixation1. Legumes comprise several evolutionary lineages derived from a common ancestor 60 million years ago (Myr ago). Papilionoids are the largest clade, dating nearly to the origin of legumes and containing most cultivated species2. Medicago truncatula is a long-established model for the study of legume biology. Here we describe the draft sequence of the M. truncatula euchromatin based on a recently completed BAC assembly supplemented with Illumina shotgun sequence, together capturing ~94% of all M. truncatula genes. A whole-genome duplication (WGD) approximately 58 Myr ago had a major role in shaping the M. truncatula genome and thereby contributed to the evolution of endosymbiotic nitrogen fixation. Subsequent to the WGD, the M. truncatula genome experienced higher levels of rearrangement than two other sequenced legumes, Glycine max and Lotus japonicus. M. truncatula is a close relative of alfalfa (Medicago sativa), a widely cultivated crop with limited genomics tools and complex autotetraploid genetics. As such, the M. truncatula genome sequence provides significant opportunities to expand alfalfa’s genomic toolbox.

Studies have indicated that exon and intron size and intergenic distance are correlated with gene expression levels and expression breadth. Previous reports on these correlations in plants and animals have been conflicting. In this study, next-generation sequence data, which has been shown to be more sensitive than previous expression profiling technologies, were generated and analyzed from 14 tissues. Our results revealed a novel dichotomy. At the low expression level, an increase in expression breadth correlated with an increase in transcript size because of an increase in the number of exons and introns. No significant changes in intron or exon sizes were noted. Conversely, genes expressed at the intermediate to high expression levels displayed a decrease in transcript size as their expression breadth increased. This was due to smaller exons, with no significant change in the number of exons. Taking advantage of the known gene space of soybean, we evaluated the positioning of genes and found significant clustering of similarly expressed genes. Identifying the correlations between the physical parameters of individual genes could lead to uncovering the role of regulation owing to nucleotide composition, which might have potential impacts in discerning the role of the noncoding regions.

Monozygotic or 'identical' twins have been widely studied to dissect the relative contributions of genetics and environment in human diseases. In multiple sclerosis (MS), an autoimmune demyelinating disease and common cause of neurodegeneration and disability in young adults, disease discordance in monozygotic twins has been interpreted to indicate environmental importance in its pathogenesis. However, genetic and epigenetic differences between monozygotic twins have been described, challenging the accepted experimental model in disambiguating the effects of nature and nurture. Here we report the genome sequences of one MS-discordant monozygotic twin pair, and messenger RNA transcriptome and epigenome sequences of CD4(+) lymphocytes from three MS-discordant, monozygotic twin pairs. No reproducible differences were detected between co-twins among approximately 3.6 million single nucleotide polymorphisms (SNPs) or approximately 0.2 million insertion-deletion polymorphisms. Nor were any reproducible differences observed between siblings of the three twin pairs in HLA haplotypes, confirmed MS-susceptibility SNPs, copy number variations, mRNA and genomic SNP and insertion-deletion genotypes, or the expression of approximately 19,000 genes in CD4(+) T cells. Only 2 to 176 differences in the methylation of approximately 2 million CpG dinucleotides were detected between siblings of the three twin pairs, in contrast to approximately 800 methylation differences between T cells of unrelated individuals and several thousand differences between tissues or between normal and cancerous tissues. In the first systematic effort to estimate sequence variation among monozygotic co-twins, we did not find evidence for genetic, epigenetic or transcriptome differences that explained disease discordance. These are the first, to our knowledge, female, twin and autoimmune disease individual genome sequences reported.

Nodulation is the result of a mutualistic interaction between legumes and symbiotic soil bacteria (e.g. soybean [Glycine max] and Bradyrhizobium japonicum) initiated by the infection of plant root hair cells by the symbiont. Fewer than 20 plant genes involved in the nodulation process have been functionally characterized. Considering the complexity of the symbiosis, significantly more genes are likely involved. To identify genes involved in root hair cell infection, we performed a large-scale transcriptome analysis of B. japonicum-inoculated and mock-inoculated soybean root hairs using three different technologies: microarray hybridization, Illumina sequencing, and quantitative real-time reverse transcription-polymerase chain reaction. Together, a total of 1,973 soybean genes were differentially expressed with high significance during root hair infection, including orthologs of previously characterized root hair infection-related genes such as NFR5 and NIN. The regulation of 60 genes was confirmed by quantitative real-time reverse transcription-polymerase chain reaction. Our analysis also highlighted changes in the expression pattern of some homeologous and tandemly duplicated soybean genes, supporting their rapid specialization.

This year marks the essential completion of the genome sequences of soybean ( Glycine max ), barrel medic ( Medicago truncatula ), and birdsfoot trefoil ( Lotus japonicus ). The impact of these assembled, annotated genomes will be enormous. Birdsfoot trefoil and barrel medic, both forage crops, are

Gram-negative soil bacteria (rhizobia) within the Rhizobiaceae phylogenetic family (alpha-proteobacteria) have the unique ability to infect and establish a nitrogen-fixing symbiosis on the roots of leguminous plants. This symbiosis is of agronomic importance, reducing the need for nitrogen fertilizer for agriculturally important plants (e.g. soybean and alfalfa). The establishment of the symbiosis involves a complex interplay between host and symbiont, resulting in the formation of a novel organ, the nodule, which the bacteria colonize as intracellular symbionts. This review focuses on the most recent discoveries relating to how this symbiosis is established. Two general developments have contributed to the recent explosion of research progress in this area: first, the adoption of two genetic model legumes, Medicago truncatula and Lotus japonicus, and second, the application of modern methods in functional genomics (e.g. transcriptomic, proteomic and metabolomic analyses).

An international consortium is sequencing the euchromatic genespace of Medicago truncatula. Extensive bioinformatic and database resources support the marker-anchored bacterial artificial chromosome (BAC) sequencing strategy. Existing physical and genetic maps and deep BAC-end sequencing help to guide the sequencing effort, while EST databases provide essential resources for genome annotation as well as transcriptome characterization and microarray design. Finished BAC sequences are joined into overlapping sequence assemblies and undergo an automated annotation process that integrates ab initio predictions with EST, protein, and other recognizable features. Because of the sequencing project's international and collaborative nature, data production, storage, and visualization tools are broadly distributed. This paper describes databases and Web resources for the project, which provide support for physical and genetic maps, genome sequence assembly, gene prediction, and integration of EST data. A central project Web site at medicago.org/genome provides access to genome viewers and other resources project-wide, including an Ensembl implementation at medicago.org, physical map and marker resources at mtgenome.ucdavis.edu, and genome viewers at the University of Oklahoma (www.genome.ou.edu), the Institute for Genomic Research (www.tigr.org), and Munich Information for Protein Sequences Center (mips.gsf.de).

Exposure of cell suspension cultures of Medicago truncatula Gaerth. to methyl jasmonate (MeJA) resulted in up to 50-fold induction of transcripts encoding the key triterpene biosynthetic enzyme beta-amyrin synthase (betaAS; EC 5.4.99.-). Transcripts reached maximum levels at 24 h post-elicitation with 0.5 mM MeJA. The entry point enzymes into the phenylpropanoid and flavonoid pathways, L: -phenylalanine ammonia-lyase (PAL; EC 4.3.1.5) and chalcone synthase (CHS; EC 2.3.1.74), respectively, were not induced by MeJA. In contrast, exposure of cells to yeast elicitor (YE) resulted in up to 45- and 14-fold induction of PAL and CHS transcripts, respectively, at only 2 h post-elicitation. betaAS transcripts were weakly induced at 12 h after exposure to YE. Over 30 different triterpene saponins were identified in the cultures, many of which were strongly induced by MeJA, but not by YE. In contrast, cinnamic acids, benzoic acids and isoflavone-derived compounds accumulated following exposure of cultures to YE, but few changes in phenylpropanoid levels were observed in response to MeJA. DNA microarray analysis confirmed the strong differential transcriptional re-programming of the cell cultures for multiple genes in the phenylpropanoid and triterpene pathways in response to MeJA and YE, and indicated different responses of individual members of gene families. This work establishes Medicago cell cultures as an excellent model for future genomics approaches to understand the regulation of legume secondary metabolism.

GC-MS-based metabolite profiling was used to analyse the response of Medicago truncatula cell cultures to elicitation with methyl jasmonate (MeJa), yeast elicitor (YE), or ultraviolet light (UV). Marked changes in the levels of primary metabolites, including several amino acids, organic acids, and carbohydrates, were observed following elicitation with MeJa. A similar, but attenuated response was observed following YE elicitation, whereas little response was observed following UV elicitation. MeJa induced the accumulation of the triterpene beta-amyrin, a precursor to the triterpene saponins, and LC-MS analysis confirmed the accumulation of triterpene saponins in MeJa-elicited samples. In addition, YE induced a slight, but significant accumulation of shikimic acid, an early precursor to the phenylpropanoid pathway, which was also demonstrated to be YE-inducible by LC-MS analyses. Correlation analyses of metabolite relationships revealed perturbation of the glycine, serine, and threonine biosynthetic pathway, and suggested the induction of threonine aldolase activity, an enzyme as yet uncharacterized from plants. Members of the branched chain amino acid pathway accumulated in a concerted fashion, with the strongest correlation being that between leucine and isoleucine (r2=0.941). While UV exposure itself had little effect on primary metabolites, the experimental procedure, as revealed by control treatments, induced changes in several metabolites which were similar to those following MeJa elicitation. Sucrose levels were lower in MJ- and YE-elicited samples compared with control samples, suggesting that a portion of the effects observed on the primary metabolic pool are a consequence of fundamental metabolic repartitioning of carbon resources rather than elicitor-specific induction. In addition, beta-alanine levels were elevated in all elicited samples, which, when viewed in the context of other elicitation responses, suggests the altered metabolism of coenzyme A and its esters, which are essential in secondary metabolism.

Glandular trichomes are a major site of plant natural product synthesis and accumulation for protection against insect predation. However, to date few studies have attempted to obtain a global view of trichome gene expression. Two contrasting approaches have been adopted to investigate genes expressed in glandular trichomes from alfalfa (Medicago sativa L.). In the first approach, 5,674 clones from an alfalfa glandular trichome cDNA library were sequenced. The most highly abundant expressed sequence tag (EST) corresponded to a lipid transfer protein. The presence of ESTs corresponding to enzymes for all steps in the biosynthesis of flavonoids suggests that these are important metabolites in alfalfa trichome biology, as confirmed by histochemistry and metabolite profiling. No ESTs corresponded to enzymes of cyclized terpenoid biosynthesis. In a second approach, microarray analysis was used to compare levels of alfalfa transcripts corresponding to 16,086 Medicago truncatula A17 genes in stems with and without trichomes. This revealed over 1,000 genes with strong preferential expression in the trichome fraction of the stem, 70% of which are of unknown function. These define a class of genes that are not trichome-specific, since M. truncatula A17 does not itself have glandular trichomes, but has potential importance for trichome function within the stem.

As a method for accurate quantitative and qualitative analysis of cell transcription, serial analysis of gene expression (SAGE) has been widely used in biological, medical, and pharmaceutical areas of research. The numerous steps required in the SAGE protocol, however, limit the efficiency of SAGE library construction. PCR amplification of ditags and concatemer separation are two such steps. Primer-2, which was used to amplify ditags in the original SAGE protocol, was prone to self-anneal and resulted in a decreased efficiency of the PCR amplification step. We have designed a new primer pair, which substantially enhances the efficiency of PCR amplification. Incorporating this modification reduces the number of reactions required for SAGE library construction. Sixty 50-μL reactions are now sufficient for SAGE library construction (up to an 80% reduction in comparison to the original protocol). Concatemer separation by column filtration and subsequent gel separation has been modified and has proven to be more efficient than simple heat treatment of the ligation mixture. Both of the above modifications may also be suitable for the newly developed LongSAGE protocol.

The genome sequence of Arabidopsis is complete and the genomes of plants representing legumes (Medicago truncatula) and grasses (rice) will soon follow. The rate at which new genes have been discovered has far outstripped the pace at which their function is determined. The greatest hurdle that plant biologists face in assigning gene function and in crop improvement is the lack of efficient and robust technologies to generate gene replacements or targeted gene knockouts. Many of the factors underlying these events remain to be elucidated. This review addresses the current status of plant gene targeting and what is known about the associated plant DNA repair mechanisms.

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.

Galanthamine is an Amaryllidaceae alkaloid used to treat the symptoms of Alzheimer’s disease. This compound is primarily isolated from daffodil (Narcissus spp.), snowdrop (Galanthus spp.), and summer snowflake (Leucojum aestivum). Despite its importance as a medicine, no genes involved in the biosynthetic pathway of galanthamine have been identified. This absence of genetic information on biosynthetic pathways is a limiting factor in the development of synthetic biology platforms for many important botanical medicines. The paucity of information is largely due to the limitations of traditional methods for finding biochemical pathway enzymes and genes in non-model organisms. A new bioinformatic approach using several recent technological improvements was applied to search for genes in the proposed galanthamine biosynthetic pathway, first targeting methyltransferases due to strong signature amino acid sequences in the proteins. Using Illumina sequencing, a de novo transcriptome assembly was constructed for daffodil. BLAST was used to identify sequences that contain signatures for plant O-methyltransferases in this transcriptome. The program HAYSTACK was then used to identify methyltransferases that fit a model for galanthamine biosynthesis in leaf, bulb and inflorescence tissues. One candidate gene for the methylation of norbelladine to 4′-O-methylnorbelladine in the proposed galanthamine biosynthetic pathway was identified. This methyltransferase cDNA was expressed in E. coli and the protein purified by affinity chromatography. The resulting protein was found to be a norbelladine 4′-O-methyltransferase (NpN4OMT) of the proposed galanthamine biosynthetic pathway.

Chimeric RNA/DNA and modified DNA oligonucleotides have been shown to direct gene-conversion events in vitro through a process involving proteins from several DNA-repair pathways. Recent experiments have extended the utility of these molecules to plants, and we previously demonstrated that plant cell-free extracts are competent to support oligonucleotide-directed genetic repair. Using this system, we are studying Arabidopsis DNA-repair mutants and the role of plant proteins in the DNA-repair process. Here we describe a method for investigating mechanisms of plastid DNA-repair pathways. Using a genetic readout system in bacteria and chimeric or modified DNA oligonucleotides designed to direct the conversion of mutations in antibiotic resistance genes, we have developed an assay for genetic repair of mutations in a spinach chloroplast lysate system. We report genetic repair of point and frameshift mutations directed by both types of modified oligonucleotides. This system enables the mechanistic study of plastid gene repair and facilitates the direct comparison between plant nuclear and organelle DNA-repair pathways.

A genetic transformation system has been developed for three Mycosphaerella pathogens of banana and plantain (Musa spp.). Mycosphaerella fijiensis and Mycosphaerella musicola, the causal agents of black and yellow Sigatoka, respectively, and Mycosphaerella eumusae, which causes Septoria leaf spot of banana, were transformed with a construct carrying a synthetic gene encoding green fluorescent protein (GFP). Most single-spored transformants that expressed GFP constitutively were mitotically stable in the absence of selection for hygromycin B resistance. Transformants of all three species were pathogenic on the susceptible banana cultivar Grand Nain, and growth in planta was comparable to wild-type strains. GFP expression by transformants allowed us to observe extensive fungal growth within leaf tissue that eventually turned necrotic, at which point the fungi grew saprophytically on the dead tissue. Leaf chlorosis and necrosis were often observed in advance of saprophytic growth of the mycelium on necrotic tissue, which supports previous reports suggesting secretion of a phytotoxin.

The Medicago Genome Initiative (MGI) is a database of EST sequences of the model legume Medicago truncatula. The database is available to the public and has resulted from a collaborative research effort between the Samuel Roberts Noble Foundation and the National Center for Genome Resources to investigate the genome of M.truncatula. MGI is part of the greater integrated Medicago functional genomics program at the Noble Foundation (http://www.noble .org), which is taking a global approach in studying the genetic and biochemical events associated with the growth, development and environmental interactions of this model legume. Our approach will include: large-scale EST sequencing, gene expression profiling, the generation of M.truncatula activation-tagged and promoter trap insertion mutants, high-throughput metabolic profiling, and proteome studies. These multidisciplinary information pools will be interfaced with one another to provide scientists with an integrated, holistic set of tools to address fundamental questions pertaining to legume biology. The public interface to the MGI database can be accessed at http://www.ncgr.org/research/mgi.

The National Center for Genome Resources (NCGR) has developed a high-throughput DNA (deoxyribonucleic acid) sequence analysis pipeline, which allows researchers at remote sites to submit biological sequence information for rapid analysis, the results of which can be queried through a Web interface. Behind the browser interface is a relational database used to manage both the raw data and the results of the different analyses performed, and a server, which performs those analyses. The system allows multiple contributors to submit data and also allows the data to be marked as "private" or as available to the general public. The CPU-intensive part of the processing is done on a 40- processor domain of a Sun Enterprise 10000 computer, which is represented by a distributed system of software objects, implemented in CORBATM (Common Object Request Broker ArchitectureTM). In this paper we discuss the architecture of the pipeline, the database support, types of DNA sequence analysis used, the distributed analysis system, and the capabilities of the Web interface. As a case study, we present data from an ongoing collaborative project in which expressed sequence tags (ESTs) from Medicago truncatula are being processed. M. truncatula is a plant that is used as a research model for crops in the legume family, an economically important group of food and forage plants.

Nodulation is the result of a mutualistic interaction between legumes and symbiotic soil bacteria (e.g. soybean [Glycine max] and Bradyrhizobium japonicum) initiated by the infection of plant root hair cells by the symbiont. Fewer than 20 plant genes involved in the nodulation process have been functionally characterized. Considering the complexity of the symbiosis, significantly more genes are likely involved. To identify genes involved in root hair cell infection, we performed a large-scale transcriptome analysis of B. japonicum-inoculated and mock-inoculated soybean root hairs using three different technologies: microarray hybridization, Illumina sequencing, and quantitative real-time reverse transcription-polymerase chain reaction. Together, a total of 1,973 soybean genes were differentially expressed with high significance during root hair infection, including orthologs of previously characterized root hair infection-related genes such as NFR5 and NIN. The regulation of 60 genes was confirmed by quantitative real-time reverse transcription-polymerase chain reaction. Our analysis also highlighted changes in the expression pattern of some homeologous and tandemly duplicated soybean genes, supporting their rapid specialization.

Chimeric oligonucleotides are synthetic molecules comprised of RNA and DNA bases assembled in a double hairpin conformation. These molecules have been shown to direct gene conversion events in mammalian cells and animals through a process involving at least one protein from the DNA mismatch repair pathway. The mechanism of action for gene repair in mammalian cells has been partially elucidated through the use of a cell-free extract system. Recent experiments have expanded the utility of chimeric oligonucleotides to plants and have demonstrated genotypic and phenotypic conversion, as well as Mendelian transmission. Although these experiments showed correction of point and frameshift mutations, the biochemical and mechanistic aspects of the process were not addressed. In this paper, we describe the establishment of cell-free extract systems from maize (Zea mays), banana (Musa acuminata cv Rasthali), and tobacco (Nicotiana tabacum). Using a genetic readout system in bacteria and chimeric oligonucleotides designed to direct the conversion of mutations in antibiotic-resistant genes, we demonstrate gene repair of point and frameshift mutations. Whereas extracts from banana and maize catalyzed repair of mutations in a precise fashion, cell-free extracts prepared from tobacco exhibited either partial repair or non-targeted nucleotide conversion. In addition, an all-DNA hairpin molecule also mediated repair albeit in an imprecise fashion in all cell-free extracts tested. This system enables the mechanistic study of gene repair in plants and may facilitate the identification of DNA repair proteins operating in plant cells.

An abundant, catalytically active β-1,3-endoglucanase (EC 3.2.1.39) has been isolated from the pulp of ripe bananas. Biochemical analysis of the purified protein, molecular modelling, and molecular cloning of the corresponding gene indicate that this banana enzyme closely resembles previously characterized plant β-glucanases with respect to its amino-acid sequence, structure and biological activity. The results described in this paper demonstrate both the occurrence of an abundant active β-1,3-endoglucanases in fruits and also readdress the question of the possible involvement of these enzymes in the ripening and/or softening process.

Self-complementary chimeric oligonucleotides (COs) composed of DNA and modified RNA residues were evaluated as a means to ( i ) create stable, site-specific base substitutions in a nuclear gene and ( ii ) introduce a frameshift in a nuclear transgene in plant cells. To demonstrate the creation of allele-specific mutations in a member of a gene family, COs were designed to target the codon for Pro-196 of SuRA, a tobacco acetolactate synthase (ALS) gene. An amino acid substitution at Pro-196 of ALS confers a herbicide-resistance phenotype that can be used as a selectable marker in plant cells. COs were designed to contain a 25-nt homology domain comprised of a five-deoxyribonucleotide region (harboring a single base mismatch to the native ALS sequence) flanked by regions each composed of 10 ribonucleotides. After recovery of herbicide-resistant tobacco cells on selective medium, DNA sequence analyses identified base conversions in the ALS gene at the codon for Pro-196. To demonstrate a site-specific insertion of a single base into a targeted gene, COs were used to restore expression of an inactive green fluorescent protein transgene that had been designed to contain a single base deletion. Recovery of fluorescent cells confirmed the deletion correction. Our results demonstrate the application of a technology to modify individual genetic loci by catalyzing either a base substitution or a base addition to specific nuclear genes; this approach should have great utility in the area of plant functional genomics.

We have identified and characterized the abundant protein from the pulp of banana fruit (Musa acuminata cv. Grand Nain), and have isolated a cDNA clone encoding this protein. Comparison of the amino terminal sequence of the purified 31 kDa protein (P31) suggests that it is related to plant chitinases. Western analyses utilizing rabbit anti-P31 antiserum demonstrate that this protein is pulp-specific in banana. A full-length cDNA clone homologous to class III acidic chitinase genes has been isolated from a pulp cDNA library by differential screening. The identity of this clone as encoding P31 was verified by comparisons between the amino-terminal peptide sequence and the cDNA sequence and cross-hybridization of the translation product of the cDNA clone with P31 antiserum. Northern and western blot analyses of RNA and protein isolated from banana pulp at different stages of ripening indicate that the cDNA and protein are expressed at high levels in the pulp of unripe fruit, and that their abundance decreases as the fruit ripens. Based on its expression pattern and deduced amino acid sequence and composition, we hypothesize that the physiological role of P31 is not for plant protection, but as a storage protein in banana pulp.

Background Theobroma cacao L. cultivar Matina 1-6 belongs to the most cultivated cacao type. The availability of its genome sequence and methods for identifying genes responsible for important cacao traits will aid cacao researchers and breeders. Results We describe the sequencing and assembly of the genome of Theobroma cacao L. cultivar Matina 1-6. The genome of the Matina 1-6 cultivar is 445 Mbp, which is significantly larger than a sequenced Criollo cultivar, and more typical of other cultivars. The chromosome-scale assembly, version 1.1, contains 711 scaffolds covering 346.0 Mbp, with a contig N50 of 84.4 kbp, a scaffold N50 of 34.4 Mbp, and an evidence-based gene set of 29,408 loci. Version 1.1 has 10x the scaffold N50 and 4x the contig N50 as Criollo, and includes 111 Mb more anchored sequence. The version 1.1 assembly has 4.4% gap sequence, while Criollo has 10.9%. Through a combination of haplotype, association mapping and gene expression analyses, we leverage this robust reference genome to identify a promising candidate gene responsible for pod color variation. We demonstrate that green/red pod color in cacao is likely regulated by the R2R3 MYB transcription factor TcMYB113, homologs of which determine pigmentation in Rosaceae, Solanaceae, and Brassicaceae. One SNP within the target site for a highly conserved trans-acting siRNA in dicots, found within TcMYB113, seems to affect transcript levels of this gene and therefore pod color variation. Conclusions We report a high-quality sequence and annotation of Theobroma cacao L. and demonstrate its utility in identifying candidate genes regulating traits.