A physical, genetic and functional sequence assembly of the barley genome

Comparative molecular cytogenetic analyses of a major tandemly repeated DNA family and retrotransposon sequences in cultivated jute Corchorus species (Malvaceae)

BACKGROUND AND AIMS

The cultivated jute species Corchorus olitorius and Corchorus capsularis are important fibre crops. The analysis of repetitive DNA sequences, comprising a major part of plant genomes, has not been carried out in jute but is useful to investigate the long-range organization of chromosomes. The aim of this study was the identification of repetitive DNA sequences to facilitate comparative molecular and cytogenetic studies of two jute cultivars and to develop a fluorescent in situ hybridization (FISH) karyotype for chromosome identification.

METHODS

A plasmid library was generated from C. olitorius and C. capsularis with genomic restriction fragments of 100-500 bp, which was complemented by targeted cloning of satellite DNA by PCR. The diversity of the repetitive DNA families was analysed comparatively. The genomic abundance and chromosomal localization of different repeat classes were investigated by Southern analysis and FISH, respectively. The cytosine methylation of satellite arrays was studied by immunolabelling.

KEY RESULTS

Major satellite repeats and retrotransposons have been identified from C. olitorius and C. capsularis. The satellite family CoSat I forms two undermethylated species-specific subfamilies, while the long terminal repeat (LTR) retrotransposons CoRetro I and CoRetro II show similarity to the Metaviridea of plant retroelements. FISH karyotypes were developed by multicolour FISH using these repetitive DNA sequences in combination with 5S and 18S-5·8S-25S rRNA genes which enable the unequivocal chromosome discrimination in both jute species.

CONCLUSIONS

The analysis of the structure and diversity of the repeated DNA is crucial for genome sequence annotation. The reference karyotypes will be useful for breeding of jute and provide the basis for karyotyping homeologous chromosomes of wild jute species to reveal the genetic and evolutionary relationship between cultivated and wild Corchorus species.

A high density physical map of chromosome 1BL supports evolutionary studies, map-based cloning and sequencing in wheat

BACKGROUND

As for other major crops, achieving a complete wheat genome sequence is essential for the application of genomics to breeding new and improved varieties. To overcome the complexities of the large, highly repetitive and hexaploid wheat genome, the International Wheat Genome Sequencing Consortium established a chromosome-based strategy that was validated by the construction of the physical map of chromosome 3B. Here, we present improved strategies for the construction of highly integrated and ordered wheat physical maps, using chromosome 1BL as a template, and illustrate their potential for evolutionary studies and map-based cloning.

RESULTS

Using a combination of novel high throughput marker assays and an assembly program, we developed a high quality physical map representing 93% of wheat chromosome 1BL, anchored and ordered with 5,489 markers including 1,161 genes. Analysis of the gene space organization and evolution revealed that gene distribution and conservation along the chromosome results from the superimposition of the ancestral grass and recent wheat evolutionary patterns, leading to a peak of synteny in the central part of the chromosome arm and an increased density of non-collinear genes towards the telomere. With a density of about 11 markers per Mb, the 1BL physical map provides 916 markers, including 193 genes, for fine mapping the 40 QTLs mapped on this chromosome.

CONCLUSIONS

Here, we demonstrate that high marker density physical maps can be developed in complex genomes such as wheat to accelerate map-based cloning, gain new insights into genome evolution, and provide a foundation for reference sequencing.

Towards positional isolation of three quantitative trait loci conferring resistance to powdery mildew in two Spanish barley landraces

Three quantitative trait loci (QTL) conferring broad spectrum resistance to powdery mildew, caused by the fungus Blumeria graminis f. sp. hordei, were previously identified on chromosomes 7HS, 7HL and 6HL in the Spanish barley landrace-derived lines SBCC097 and SBCC145. In the present work, a genome-wide putative linear gene index of barley (Genome Zipper) and the first draft of the physical, genetic and functional sequence of the barley genome were used to go one step further in the shortening and explicit demarcation on the barley genome of these regions conferring resistance to powdery mildew as well as in the identification of candidate genes. First, a comparative analysis of the target regions to the barley Genome Zippers of chromosomes 7H and 6H allowed the development of 25 new gene-based molecular markers, which slightly better delimit the QTL intervals. These new markers provided the framework for anchoring of genetic and physical maps, figuring out the outline of the barley genome at the target regions in SBCC097 and SBCC145. The outermost flanking markers of QTLs on 7HS, 7HL and 6HL defined a physical area of 4 Mb, 3.7 Mb and 3.2 Mb, respectively. In total, 21, 10 and 16 genes on 7HS, 7HL and 6HL, respectively, could be interpreted as potential candidates to explain the resistance to powdery mildew, as they encode proteins of related functions with respect to the known pathogen defense-related processes. The majority of these were annotated as belonging to the NBS-LRR class or protein kinase family.

Genomics reveals new landscapes for crop improvement

The sequencing of large and complex genomes of crop species, facilitated by new sequencing technologies and bioinformatic approaches, has provided new opportunities for crop improvement. Current challenges include understanding how genetic variation translates into phenotypic performance in the field.

User guide for mapping-by-sequencing in Arabidopsis

Mapping-by-sequencing combines genetic mapping with whole-genome sequencing in order to accelerate mutant identification. However, application of mapping-by-sequencing requires decisions on various practical settings on the experimental design that are not intuitively answered. Following an experimentally determined recombination landscape of Arabidopsis and next generation sequencing-specific biases, we simulated more than 400,000 mapping-by-sequencing experiments. This allowed us to evaluate a broad range of different types of experiments and to develop general rules for mapping-by-sequencing in Arabidopsis. Most importantly, this informs about the properties of different crossing scenarios, the number of recombinants and sequencing depth needed for successful mapping experiments.

Distribution, functional impact, and origin mechanisms of copy number variation in the barley genome

Background

There is growing evidence for the prevalence of copy number variation (CNV) and its role in phenotypic variation in many eukaryotic species. Here we use array comparative genomic hybridization to explore the extent of this type of structural variation in domesticated barley cultivars and wild barleys.

Results

A collection of 14 barley genotypes including eight cultivars and six wild barleys were used for comparative genomic hybridization. CNV affects 14.9% of all the sequences that were assessed. Higher levels of CNV diversity are present in the wild accessions relative to cultivated barley. CNVs are enriched near the ends of all chromosomes except 4H, which exhibits the lowest frequency of CNVs. CNV affects 9.5% of the coding sequences represented on the array and the genes affected by CNV are enriched for sequences annotated as disease-resistance proteins and protein kinases. Sequence-based comparisons of CNV between cultivars Barke and Morex provided evidence that DNA repair mechanisms of double-strand breaks via single-stranded annealing and synthesis-dependent strand annealing play an important role in the origin of CNV in barley.

Conclusions

We present the first catalog of CNVs in a diploid Triticeae species, which opens the door for future genome diversity research in a tribe that comprises the economically important cereal species wheat, barley, and rye. Our findings constitute a valuable resource for the identification of CNV affecting genes of agronomic importance. We also identify potential mechanisms that can generate variation in copy number in plant genomes.

The inhibitor of wax 1 locus (Iw1) prevents formation of β- and OH-β-diketones in wheat cuticular waxes and maps to a sub-cM interval on chromosome arm 2BS

Glaucousness is described as the scattering effect of visible light from wax deposited on the cuticle of plant aerial organs. In wheat, two dominant genes lead to non-glaucous phenotypes: Inhibitor of wax 1 (Iw1) and Iw2. The molecular mechanisms and the exact extent (beyond visual assessment) by which these genes affect the composition and quantity of cuticular wax is unclear. To describe the Iw1 locus we used a genetic approach with detailed biochemical characterization of wax compounds. Using synteny and a large number of F2 gametes, Iw1 was fine-mapped to a sub-cM genetic interval on wheat chromosome arm 2BS, which includes a single collinear gene from the corresponding Brachypodium and rice physical maps. The major components of flag leaf and peduncle cuticular waxes included primary alcohols, β-diketones and n-alkanes. Small amounts of C19-C27 alkyl and methylalkylresorcinols that have not previously been described in wheat waxes were identified. Using six pairs of BC2 F3 near-isogenic lines, we show that Iw1 inhibits the formation of β- and hydroxy-β-diketones in the peduncle and flag leaf blade cuticles. This inhibitory effect is independent of genetic background or tissue, and is accompanied by minor but consistent increases in n-alkanes and C24 primary alcohols. No differences were found in cuticle thickness and carbon isotope discrimination in near-isogenic lines differing at Iw1.

Host-induced gene silencing in barley powdery mildew reveals a class of ribonuclease-like effectors

Obligate biotrophic pathogens of plants must circumvent or counteract defenses to guarantee accommodation inside the host. To do so, they secrete a variety of effectors that regulate host immunity and facilitate the establishment of pathogen feeding structures called haustoria. The barley powdery mildew fungus Blumeria graminis f. sp. hordei produces a large number of proteins predicted to be secreted from haustoria. Fifty of these Blumeria effector candidates (BEC) were screened by host-induced gene silencing (HIGS), and eight were identified that contribute to infection. One shows similarity to β-1,3 glucosyltransferases, one to metallo-proteases, and two to microbial secreted ribonucleases; the remainder have no similarity to proteins of known function. Transcript abundance of all eight BEC increases dramatically in the early stages of infection and establishment of haustoria, consistent with a role in that process. Complementation analysis using silencing-insensitive synthetic cDNAs demonstrated that the ribonuclease-like BEC 1011 and 1054 are bona fide effectors that function within the plant cell. BEC1011 specifically interferes with pathogen-induced host cell death. Both are part of a gene superfamily unique to the powdery mildew fungi. Structural modeling was consistent, with BEC1054 adopting a ribonuclease-like fold, a scaffold not previously associated with effector function.

Silicon era of carbon-based life: application of genomics and bioinformatics in crop stress research

Abiotic and biotic stresses lead to massive reprogramming of different life processes and are the major limiting factors hampering crop productivity. Omics-based research platforms allow for a holistic and comprehensive survey on crop stress responses and hence may bring forth better crop improvement strategies. Since high-throughput approaches generate considerable amounts of data, bioinformatics tools will play an essential role in storing, retrieving, sharing, processing, and analyzing them. Genomic and functional genomic studies in crops still lag far behind similar studies in humans and other animals. In this review, we summarize some useful genomics and bioinformatics resources available to crop scientists. In addition, we also discuss the major challenges and advancements in the "-omics" studies, with an emphasis on their possible impacts on crop stress research and crop improvement.

Divergent evolutionary mechanisms of co-located Tak/Lrk and Glu-D3 loci revealed by comparative analysis of grass genomes

Seed storage and disease resistance proteins are major traits of wheat. The study of their gene organization and evolution has some implications in breeding. In this study, we characterized the hexaploid wheat D-genome BAC clone TaBAC703A9 that contains a low molecular weight glutenin locus (Glu-D3) and a resistance gene analogue cluster. With a gene density of one gene per 4.8 kb, the cluster contains four resistance gene analogues, namely Tak703-1, Lrr703, Tak703, and Lrk703. This structural cluster unit was conserved across nine grass genomes, but divergent evolutionary mechanisms have been involved in shaping the Tak/Lrk loci in the different species. Gene duplication was the major force for the Tak/Lrk evolution in oats, maize, barley, wheat, sorghum, and Brachypodium, while tandem duplication drove the expansion of this locus in japonica rice. Despite the close proximity of the Glu-D3 and the Tak/Lrk loci in wheat, the evolutionary mechanisms that drove their amplification differ. The Glu-D3 region had a lower gene density, and its amplification was driven by retroelements.

Assembling the 20 Gb white spruce (Picea glauca) genome from whole-genome shotgun sequencing data

White spruce (Picea glauca) is a dominant conifer of the boreal forests of North America, and providing genomics resources for this commercially valuable tree will help improve forest management and conservation efforts. Sequencing and assembling the large and highly repetitive spruce genome though pushes the boundaries of the current technology. Here, we describe a whole-genome shotgun sequencing strategy using two Illumina sequencing platforms and an assembly approach using the ABySS software. We report a 20.8 giga base pairs draft genome in 4.9 million scaffolds, with a scaffold N50 of 20 356 bp. We demonstrate how recent improvements in the sequencing technology, especially increasing read lengths and paired end reads from longer fragments have a major impact on the assembly contiguity. We also note that scalable bioinformatics tools are instrumental in providing rapid draft assemblies.

Availability: The Picea glauca genome sequencing and assembly data are available through NCBI (Accession#: ALWZ0100000000 PID: PRJNA83435). http://www.ncbi.nlm.nih.gov/bioproject/83435.

Contact:ibirol@bcgsc.ca

Supplementary information:Supplementary data are available at Bioinformatics online.

Comparative proteomic analysis of aluminum tolerance in tibetan wild and cultivated barleys

Aluminum (Al) toxicity is a major limiting factor for plant production in acid soils. Wild barley germplasm is rich in genetic diversity and may provide elite genes for crop Al tolerance improvement. The hydroponic-experiments were performed to compare proteomic and transcriptional characteristics of two contrasting Tibetan wild barley genotypes Al- resistant/tolerant XZ16 and Al-sensitive XZ61 as well as Al-resistant cv. Dayton. Results showed that XZ16 had less Al uptake and translocation than XZ61 and Dayton under Al stress. Thirty-five Al-tolerance/resistance-associated proteins were identified and categorized mainly in metabolism, energy, cell growth/division, protein biosynthesis, protein destination/storage, transporter, signal transduction, disease/defense, etc. Among them, 30 were mapped on barley genome, with 16 proteins being exclusively up-regulated by Al stress in XZ16, including 4 proteins (S-adenosylmethionine-synthase 3, ATP synthase beta subunit, triosephosphate isomerase, Bp2A) specifically expressed in XZ16 but not Dayton. The findings highlighted the significance of specific-proteins associated with Al tolerance, and verified Tibetan wild barley as a novel genetic resource for Al tolerance.

dbWFA: a web-based database for functional annotation of Triticum aestivum transcripts

The functional annotation of genes based on sequence homology with genes from model species genomes is time-consuming because it is necessary to mine several unrelated databases. The aim of the present work was to develop a functional annotation database for common wheat Triticum aestivum (L.). The database, named dbWFA, is based on the reference NCBI UniGene set, an expressed gene catalogue built by expressed sequence tag clustering, and on full-length coding sequences retrieved from the TriFLDB database. Information from good-quality heterogeneous sources, including annotations for model plant species Arabidopsis thaliana (L.) Heynh. and Oryza sativa L., was gathered and linked to T. aestivum sequences through BLAST-based homology searches. Even though the complexity of the transcriptome cannot yet be fully appreciated, we developed a tool to easily and promptly obtain information from multiple functional annotation systems (Gene Ontology, MapMan bin codes, MIPS Functional Categories, PlantCyc pathway reactions and TAIR gene families). The use of dbWFA is illustrated here with several query examples. We were able to assign a putative function to 45% of the UniGenes and 81% of the full-length coding sequences from TriFLDB. Moreover, comparison of the annotation of the whole T. aestivum UniGene set along with curated annotations of the two model species assessed the accuracy of the annotation provided by dbWFA. To further illustrate the use of dbWFA, genes specifically expressed during the early cell division or late storage polymer accumulation phases of T. aestivum grain development were identified using a clustering analysis and then annotated using dbWFA. The annotation of these two sets of genes was consistent with previous analyses of T. aestivum grain transcriptomes and proteomes.

Database URL:urgi.versailles.inra.fr/dbWFA/

Prospects for advancing defense to cereal rusts through genetical genomics

Rusts are one of the most severe threats to cereal crops because new pathogen races emerge regularly, resulting in infestations that lead to large yield losses. In 1999, a new race of stem rust, Puccinia graminis f. sp. tritici (Pgt TTKSK or Ug99), was discovered in Uganda. Most of the wheat and barley cultivars grown currently worldwide are susceptible to this new race. Pgt TTKSK has already spread northward into Iran and will likely spread eastward throughout the Indian subcontinent in the near future. This scenario is not unique to stem rust; new races of leaf rust (Puccinia triticina) and stripe rust (Puccinia striiformis) have also emerged recently. One strategy for countering the persistent adaptability of these pathogens is to stack complete- and partial-resistance genes, which requires significant breeding efforts in order to reduce deleterious effects of linkage drag. These varied resistance combinations are typically more difficult for the pathogen to defeat, since they would be predicted to apply lower selection pressure. Genetical genomics or expression Quantitative Trait Locus (eQTL) analysis enables the identification of regulatory loci that control the expression of many to hundreds of genes. Integrated deployment of these technologies coupled with efficient phenotyping offers significant potential to elucidate the regulatory nodes in genetic networks that orchestrate host defense responses. The focus of this review will be to present advances in genetical genomic experimental designs and analysis, particularly as they apply to the prospects for discovering partial disease resistance alleles in cereals.

Lunasin in cereal seeds: What is the origin?

Lunasin is a peptide from soybean seeds which has been demonstrated to have anticancer properties. It has also been reported in cereal seeds: wheat, rye, barley and Triticale. However, extensive searches of transcriptome and DNA sequence databases for wheat and other cereals have failed to identify sequences encoding either the lunasin peptide or a precursor protein. This raises the question of the origin of the lunasin reported in cereal grain.

Quantitative high resolution mapping of HvMLH3 foci in barley pachytene nuclei reveals a strong distal bias and weak interference

In barley (Hordeum vulgare L.), chiasmata (the physical sites of genetic crossovers) are skewed towards the distal ends of chromosomes, effectively consigning a large proportion of genes to recombination coldspots. This has the effect of limiting potential genetic variability, and of reducing the efficiency of map-based cloning and breeding approaches for this crop. Shifting the sites of recombination to more proximal chromosome regions by forward and reverse genetic means may be profitable in terms of realizing the genetic potential of the species, but is predicated upon a better understanding of the mechanisms governing the sites of these events, and upon the ability to recognize real changes in recombination patterns. The barley MutL Homologue (HvMLH3), a marker for class I interfering crossovers, has been isolated and a specific antibody has been raised. Immunolocalization of HvMLH3 along with the synaptonemal complex transverse filament protein ZYP1, used in conjunction with fluorescence in situ hybridization (FISH) tagging of specific barley chromosomes, has enabled access to the physical recombination landscape of the barley cultivars Morex and Bowman. Consistent distal localization of HvMLH3 foci throughout the genome, and similar patterns of HvMLH3 foci within bivalents 2H and 3H have been observed. A difference in total numbers of HvMLH3 foci between these two cultivars has been quantified, which is interpreted as representing genotypic variation in class I crossover frequency. Discrepancies between the frequencies of HvMLH3 foci and crossover frequencies derived from linkage analysis point to the existence of at least two crossover pathways in barley. It is also shown that interference of HvMLH3 foci is relatively weak compared with other plant species.

Capacity and plasticity of potassium channels and high-affinity transporters in roots of barley and Arabidopsis

The role of potassium (K(+)) transporters in high- and low-affinity K(+) uptake was examined in roots of intact barley (Hordeum vulgare) and Arabidopsis (Arabidopsis thaliana) plants by use of (42)K radiotracing, electrophysiology, pharmacology, and mutant analysis. Comparisons were made between results from barley and five genotypes of Arabidopsis, including single and double knockout mutants for the high-affinity transporter, AtHAK5, and the Shaker-type channel, AtAKT1. In Arabidopsis, steady-state K(+) influx at low external K(+) concentration ([K(+)]ext = 22.5 µm) was predominantly mediated by AtAKT1 when high-affinity transport was inhibited by ammonium, whereas in barley, by contrast, K(+) channels could not operate below 100 µm. Withdrawal of ammonium resulted in an immediate and dramatic stimulation of K(+) influx in barley, indicating a shift from active to passive K(+) uptake at low [K(+)]ext and yielding fluxes as high as 36 µmol g (root fresh weight)(-1) h(-1) at 5 mm [K(+)]ext, among the highest transporter-mediated K(+) fluxes hitherto reported. This ammonium-withdrawal effect was also established in all Arabidopsis lines (the wild types, atakt1, athak5, and athak5 atakt1) at low [K(+)]ext, revealing the concerted involvement of several transport systems. The ammonium-withdrawal effect coincided with a suppression of K(+) efflux and a significant hyperpolarization of the plasma membrane in all genotypes except athak5 atakt1, could be sustained over 24 h, and resulted in increased tissue K(+) accumulation. We discuss key differences and similarities in K(+) acquisition between two important model systems and reveal novel aspects of K(+) transport in planta.

A 4-gigabase physical map unlocks the structure and evolution of the complex genome of Aegilops tauschii, the wheat D-genome progenitor

The current limitations in genome sequencing technology require the construction of physical maps for high-quality draft sequences of large plant genomes, such as that of Aegilops tauschii, the wheat D-genome progenitor. To construct a physical map of the Ae. tauschii genome, we fingerprinted 461,706 bacterial artificial chromosome clones, assembled contigs, designed a 10K Ae. tauschii Infinium SNP array, constructed a 7,185-marker genetic map, and anchored on the map contigs totaling 4.03 Gb. Using whole genome shotgun reads, we extended the SNP marker sequences and found 17,093 genes and gene fragments. We showed that collinearity of the Ae. tauschii genes with Brachypodium distachyon, rice, and sorghum decreased with phylogenetic distance and that structural genome evolution rates have been high across all investigated lineages in subfamily Pooideae, including that of Brachypodieae. We obtained additional information about the evolution of the seven Triticeae chromosomes from 12 ancestral chromosomes and uncovered a pattern of centromere inactivation accompanying nested chromosome insertions in grasses. We showed that the density of noncollinear genes along the Ae. tauschii chromosomes positively correlates with recombination rates, suggested a cause, and showed that new genes, exemplified by disease resistance genes, are preferentially located in high-recombination chromosome regions.

Identification and expression analyses of cytosolic glutamine synthetase genes in barley (Hordeum vulgare L.)

Glutamine synthetase (GS) is a key enzyme in nitrogen (N) assimilation, particularly during seed development. Three cytosolic GS isoforms (HvGS1) were identified in barley (Hordeum vulgare L. cv Golden Promise). Quantitation of gene expression, localization and response to N supply revealed that each gene plays a non-redundant role in different tissues and during development. Localization of HvGS1_1 in vascular cells of different tissues, combined with its abundance in the stem and its response to changes in N supply, indicate that it is important in N transport and remobilization. HvGS1_1 is located on chromosome 6H at 72.54 cM, close to the marker HVM074 which is associated with a major quantitative trait locus (QTL) for grain protein content (GPC). HvGS1_1 may be a potential candidate gene to manipulate barley GPC. HvGS1_2 mRNA was localized to the leaf mesophyll cells, in the cortex and pericycle of roots, and was the dominant HvGS1 isoform in these tissues. HvGS1_2 expression increased in leaves with an increasing supply of N, suggesting its role in the primary assimilation of N. HvGS1_3 was specifically and predominantly localized in the grain, being highly expressed throughout grain development. HvGS1_3 expression increased specifically in the roots of plants grown on high NH(+)4, suggesting that it has a primary role in grain N assimilation and also in the protection against ammonium toxicity in roots. The expression of HvGS1 genes is directly correlated with protein and enzymatic activity, indicating that transcriptional regulation is of prime importance in the control of GS activity in barley.

Transcriptome-wide identification of R2R3-MYB transcription factors in barley with their boron responsive expression analysis

MYB family of transcription factors (TF) comprises one of the largest transcription factors in plants and is represented in all eukaryotes. They include highly conserved MYB repeats (1R, R2R3, 3R, and 4R) in the N-terminus. In addition to this, they have diverse C-terminal sequences which help the protein gain wide distinct functions, such as controlling development, secondary metabolism, hormonal regulation and response to biotic and abiotic stress. Stress-responsive roles of the MYB TFs were reported for drought, salt, wounding, cold, freezing, dehydration and osmotic stresses. This study describes the identification of barley R2R3-MYB TFs including their expression analysis in tissues under control and Boron (B) toxic conditions. Conserved motifs for MYB proteins were searched into barley full-transcriptome RNA-seq data and a total of 320 protein sequences were filtered as MYB TFs in which 51 of them corresponded to R2R3 MYB TFs. Using various bioinformatics tools, their conserved domain structures, chromosomal distributions, gene duplications, comparative functional analysis, as well as phylogenetic relations with Arabidopsis thaliana, were conducted. Beside the RNA-seq data-based expression pattern analysis of 51 R2R3 MYB TFs, quantitative analysis of selected R2R3 MYB TF genes was assessed in control and B-stressed root and leaf tissues. Critical B-induced R2R3 MYB TFs were identified. It was concluded that the results would be useful for functional characterizations of R2R3-type MYB transcription factors that are possibly involved in both B stress and divergent regulation mechanisms in plants.

Draft genome of the wheat A-genome progenitor Triticum urartu

Bread wheat (Triticum aestivum, AABBDD) is one of the most widely cultivated and consumed food crops in the world. However, the complex polyploid nature of its genome makes genetic and functional analyses extremely challenging. The A genome, as a basic genome of bread wheat and other polyploid wheats, for example, T. turgidum (AABB), T. timopheevii (AAGG) and T. zhukovskyi (AAGGA(m)A(m)), is central to wheat evolution, domestication and genetic improvement. The progenitor species of the A genome is the diploid wild einkorn wheat T. urartu, which resembles cultivated wheat more extensively than do Aegilops speltoides (the ancestor of the B genome) and Ae. tauschii (the donor of the D genome), especially in the morphology and development of spike and seed. Here we present the generation, assembly and analysis of a whole-genome shotgun draft sequence of the T. urartu genome. We identified protein-coding gene models, performed genome structure analyses and assessed its utility for analysing agronomically important genes and for developing molecular markers. Our T. urartu genome assembly provides a diploid reference for analysis of polyploid wheat genomes and is a valuable resource for the genetic improvement of wheat.

Advances in biotechnology and informatics to link variation in the genome to phenotypes in plants and animals

Advances in our understanding of genome structure provide consistent evidence for the existence of a core genome representing species classically defined by phenotype, as well as conditionally dispensable components of the genome that shows extensive variation between individuals of a given species. Generally, conservation of phenotypic features between species reflects conserved features of the genome; however, this is evidently not necessarily always the case as demonstrated by the analysis of the tunicate chordate Oikopleura dioica. In both plants and animals, the methylation activity of DNA and histones continues to present new variables for modifying (eventually) the phenotype of an organism and provides for structural variation that builds on the point mutations, rearrangements, indels, and amplification of retrotransposable elements traditionally considered. The translation of the advances in the structure/function analysis of the genome to industry is facilitated through the capture of research outputs in "toolboxes" that remain accessible in the public domain.

MORE SPIKELETS1 is required for spikelet fate in the inflorescence of Brachypodium

Grasses produce florets on a structure called a spikelet, and variation in the number and arrangement of both branches and spikelets contributes to the great diversity of grass inflorescence architecture. In Brachypodium (Brachypodium distachyon), the inflorescence is an unbranched spike with a terminal spikelet and a limited number of lateral spikelets. Spikelets are indeterminate and give rise to a variable number of florets. Here, we provide a detailed description of the stages of inflorescence development in Brachypodium. To gain insight into the genetic regulation of Brachypodium inflorescence development, we generated fast neutron mutant populations and screened for phenotypic mutants. Among the mutants identified, the more spikelets1 (mos1) mutant had an increased number of axillary meristems produced from inflorescence meristem compared with the wild type. These axillary meristems developed as branches with production of higher order spikelets. Using a candidate gene approach, mos1 was found to have a genomic rearrangement disrupting the expression of an ethylene response factor class of APETALA2 transcription factor related to the spikelet meristem identity genes branched silkless1 (bd1) in maize (Zea mays) and FRIZZY PANICLE (FZP) in rice (Oryza sativa). We propose MOS1 likely corresponds to the Brachypodium bd1 and FZP ortholog and that the function of this gene in determining spikelet meristem fate is conserved with distantly related grass species. However, MOS1 also appears to be involved in the timing of initiation of the terminal spikelet. As such, MOS1 may regulate the transition to terminal spikelet development in other closely related and agriculturally important species, particularly wheat (Triticum aestivum).

A graph-theoretical approach to the selection of the minimum tiling path from a physical map

The problem of computing the minimum tiling path (MTP) from a set of clones arranged in a physical map is a cornerstone of hierarchical (clone-by-clone) genome sequencing projects. We formulate this problem in a graph theoretical framework, and then solve by a combination of minimum hitting set and minimum spanning tree algorithms. The tool implementing this strategy, called FMTP, shows improved performance compared to the widely used software FPC. When we execute FMTP and FPC on the same physical map, the MTP produced by FMTP covers a higher portion of the genome, and uses a smaller number of clones. For instance, on the rice genome the MTP produced by our tool would reduce by about 11 percent the cost of a clone-by-clone sequencing project. Source code, benchmark data sets, and documentation of FMTP are freely available at >http://code.google.com/p/fingerprint-based-minimal-tiling-path/ under MIT license.

A root-specific wall-associated kinase gene, HvWAK1, regulates root growth and is highly divergent in barley and other cereals

Wall-associated receptor-like kinases (WAKs) are important candidates for directly linking the extracellular matrix with intracellular compartments and are involved in developmental processes and stress response. WAK gene family has been identified in plants such as Arabidopsis and rice. Here, we present a detailed analysis of the WAK1 gene from barley cv. Golden Promise, mapped to chromosome 5H. Three BAC clones corresponding to the WAK fragment were sequenced and the full-length WAK1 gene was characterized. The gene has three exons and two short introns with a coding region of 2,178 bp encoding a protein of 725 amino acids. A regulatory region was analyzed in -1,000 bp sequence upstream to start codon. Using conserved domains database and SMART, various conserved domains such as GUB WAK Bind, epidermal growth factor CA, and protein kinase C as well as other regions like signal peptides, active sites, and transmembrane domains were identified. The gene organization of HvWAK1 was compared with wheat (TaWAK1) and Arabidopsis (AtWAK1), suggesting that the WAK1 gene organization has remained highly conserved. Nonetheless, WAK1 was found to be highly divergent when compared with sequences available from barley cv. Haruna Nijo (50 %), rice (46 %), wheat (21 %), Arabidopsis (25 %), and maize (19 %). This divergence may have facilitated a better adaptation to surrounding environments due to its role in communication between the extracellular matrix, cell, and outer environment. Semiquantitative RT-PCR-based expression analysis indicates HvWAK1 expression is specific to roots. Significant differences in root growth between GP wild type and GP-Ds mutant seedlings were observed under control and salt stress conditions.

Draft genome sequence of the rubber tree Hevea brasiliensis

Background

Hevea brasiliensis, a member of the Euphorbiaceae family, is the major commercial source of natural rubber (NR). NR is a latex polymer with high elasticity, flexibility, and resilience that has played a critical role in the world economy since 1876.

Results

Here, we report the draft genome sequence of H. brasiliensis. The assembly spans ~1.1 Gb of the estimated 2.15 Gb haploid genome. Overall, ~78% of the genome was identified as repetitive DNA. Gene prediction shows 68,955 gene models, of which 12.7% are unique to Hevea. Most of the key genes associated with rubber biosynthesis, rubberwood formation, disease resistance, and allergenicity have been identified.

Conclusions

The knowledge gained from this genome sequence will aid in the future development of high-yielding clones to keep up with the ever increasing need for natural rubber.

An improved genotyping by sequencing (GBS) approach offering increased versatility and efficiency of SNP discovery and genotyping

Highly parallel SNP genotyping platforms have been developed for some important crop species, but these platforms typically carry a high cost per sample for first-time or small-scale users. In contrast, recently developed genotyping by sequencing (GBS) approaches offer a highly cost effective alternative for simultaneous SNP discovery and genotyping. In the present investigation, we have explored the use of GBS in soybean. In addition to developing a novel analysis pipeline to call SNPs and indels from the resulting sequence reads, we have devised a modified library preparation protocol to alter the degree of complexity reduction. We used a set of eight diverse soybean genotypes to conduct a pilot scale test of the protocol and pipeline. Using ApeKI for GBS library preparation and sequencing on an Illumina GAIIx machine, we obtained 5.5 M reads and these were processed using our pipeline. A total of 10,120 high quality SNPs were obtained and the distribution of these SNPs mirrored closely the distribution of gene-rich regions in the soybean genome. A total of 39.5% of the SNPs were present in genic regions and 52.5% of these were located in the coding sequence. Validation of over 400 genotypes at a set of randomly selected SNPs using Sanger sequencing showed a 98% success rate. We then explored the use of selective primers to achieve a greater complexity reduction during GBS library preparation. The number of SNP calls could be increased by almost 40% and their depth of coverage was more than doubled, thus opening the door to an increase in the throughput and a significant decrease in the per sample cost. The approach to obtain high quality SNPs developed here will be helpful for marker assisted genomics as well as assessment of available genetic resources for effective utilisation in a wide number of species.

Automated conserved non-coding sequence (CNS) discovery reveals differences in gene content and promoter evolution among grasses

Conserved non-coding sequences (CNS) are islands of non-coding sequence that, like protein coding exons, show less divergence in sequence between related species than functionless DNA. Several CNSs have been demonstrated experimentally to function as cis-regulatory regions. However, the specific functions of most CNSs remain unknown. Previous searches for CNS in plants have either anchored on exons and only identified nearby sequences or required years of painstaking manual annotation. Here we present an open source tool that can accurately identify CNSs between any two related species with sequenced genomes, including both those immediately adjacent to exons and distal sequences separated by >12 kb of non-coding sequence. We have used this tool to characterize new motifs, associate CNSs with additional functions, and identify previously undetected genes encoding RNA and protein in the genomes of five grass species. We provide a list of 15,363 orthologous CNSs conserved across all grasses tested. We were also able to identify regulatory sequences present in the common ancestor of grasses that have been lost in one or more extant grass lineages. Lists of orthologous gene pairs and associated CNSs are provided for reference inbred lines of arabidopsis, Japonica rice, foxtail millet, sorghum, brachypodium, and maize.

New insights into the dynamics of plant cell nuclei and chromosomes

The plant lamin-like protein NMCP/AtLINC and orthologues of the SUN-KASH complex across the nuclear envelope (NE) show the universality of nuclear structure in eukaryotes. However, depletion of components in the connection complex of the NE in plants does not induce severe defects, unlike in animals. Appearance of the Rabl configuration is not dependent on genome size in plant species. Topoisomerase II and condensin II are not essential for plant chromosome condensation. Plant endoreduplication shares several common characteristics with animals, including involvement of cyclin-dependent kinases and E2F transcription factors. Recent finding regarding endomitosis regulator GIG1 shed light on the suppression mechanism of endomitosis in plants. The robustness of plants, compared with animals, is reflected in their genome redundancy. Spatiotemporal functional analyses using chromophore-assisted light inactivation, super-resolution microscopy, and 4D (3D plus time) imaging will reveal new insights into plant nuclear and chromosomal dynamics.

Analysis of the bread wheat genome using whole-genome shotgun sequencing

Bread wheat (Triticum aestivum) is a globally important crop, accounting for 20 per cent of the calories consumed by humans. Major efforts are underway worldwide to increase wheat production by extending genetic diversity and analysing key traits, and genomic resources can accelerate progress. But so far the very large size and polyploid complexity of the bread wheat genome have been substantial barriers to genome analysis. Here we report the sequencing of its large, 17-gigabase-pair, hexaploid genome using 454 pyrosequencing, and comparison of this with the sequences of diploid ancestral and progenitor genomes. We identified between 94,000 and 96,000 genes, and assigned two-thirds to the three component genomes (A, B and D) of hexaploid wheat. High-resolution synteny maps identified many small disruptions to conserved gene order. We show that the hexaploid genome is highly dynamic, with significant loss of gene family members on polyploidization and domestication, and an abundance of gene fragments. Several classes of genes involved in energy harvesting, metabolism and growth are among expanded gene families that could be associated with crop productivity. Our analyses, coupled with the identification of extensive genetic variation, provide a resource for accelerating gene discovery and improving this major crop.

MIPS PlantsDB: a database framework for comparative plant genome research

The rapidly increasing amount of plant genome (sequence) data enables powerful comparative analyses and integrative approaches and also requires structured and comprehensive information resources. Databases are needed for both model and crop plant organisms and both intuitive search/browse views and comparative genomics tools should communicate the data to researchers and help them interpret it. MIPS PlantsDB (http://mips.helmholtz-muenchen.de/plant/genomes.jsp) was initially described in NAR in 2007 [Spannagl,M., Noubibou,O., Haase,D., Yang,L., Gundlach,H., Hindemitt, T., Klee,K., Haberer,G., Schoof,H. and Mayer,K.F. (2007) MIPSPlantsDB–plant database resource for integrative and comparative plant genome research. Nucleic Acids Res., 35, D834–D840] and was set up from the start to provide data and information resources for individual plant species as well as a framework for integrative and comparative plant genome research. PlantsDB comprises database instances for tomato, Medicago, Arabidopsis, Brachypodium, Sorghum, maize, rice, barley and wheat. Building up on that, state-of-the-art comparative genomics tools such as CrowsNest are integrated to visualize and investigate syntenic relationships between monocot genomes. Results from novel genome analysis strategies targeting the complex and repetitive genomes of triticeae species (wheat and barley) are provided and cross-linked with model species. The MIPS Repeat Element Database (mips-REdat) and Catalog (mips-REcat) as well as tight connections to other databases, e.g. via web services, are further important components of PlantsDB.

Integrating cereal genomics to support innovation in the Triticeae

The genomic resources of small grain cereals that include some of the most important crop species such as wheat, barley, and rye are attaining a level of completion that now is contributing to new structural and functional studies as well as refining molecular marker development and mapping strategies for increasing the efficiency of breeding processes. The integration of new efforts to obtain reference sequences in bread wheat and barley, in particular, is accelerating the acquisition and interpretation of genome-level analyses in both of these major crops.

Integrating cereal genomics to support innovation in the Triticeae

The genomic resources of small grain cereals that include some of the most important crop species such as wheat, barley, and rye are attaining a level of completion that now is contributing to new structural and functional studies as well as refining molecular marker development and mapping strategies for increasing the efficiency of breeding processes. The integration of new efforts to obtain reference sequences in bread wheat and barley, in particular, is accelerating the acquisition and interpretation of genome-level analyses in both of these major crops.