Complete Genome of Agrobacterium fabrum Strain 1D1416

This work reports the draft genome of Agrobacterium fabrum strain 1D1416. The assembled genome is composed of a 2,837,379-bp circular chromosome, a 2,043,296-bp linear chromosome, a 519,735-bp AT1 plasmid, a 188,396-bp AT2 plasmid, and a 196,706-bp Ti virulence plasmid. The nondisarmed strain produces gall-like structures in citrus tissue.

The CC-NB-LRR protein BSR1 from Brachypodium confers resistance to Barley stripe mosaic virus in gramineous plants by recognising TGB1 movement protein

Although some nucleotide binding, leucine-rich repeat immune receptor (NLR) proteins conferring resistance to specific viruses have been identified in dicot plants, NLR proteins involved in viral resistance have not been described in monocots. We have used map-based cloning to isolate the CC-NB-LRR (CNL) Barley stripe mosaic virus (BSMV) resistance gene barley stripe resistance 1 (BSR1) from Brachypodium distachyon Bd3-1 inbred line. Stable BSR1 transgenic Brachypodium line Bd21-3, barley (Golden Promise) and wheat (Kenong 199) plants developed resistance against BSMV ND18 strain. Allelic variation analyses indicated that BSR1 is present in several Brachypodium accessions collected from countries in the Middle East. Protein domain swaps revealed that the intact LRR domain and the C-terminus of BSR1 are required for resistance. BSR1 interacts with the BSMV ND18 TGB1 protein in planta and shows temperature-sensitive antiviral resistance. The R³⁹⁰ and T³⁹² residues of TGB1_ND (ND18 strain) and the G¹⁹⁶ and K¹⁹⁷ residues within the BSR1 P-loop motif are key amino acids required for immune activation. BSR1 is the first cloned virus resistance gene encoding a typical CNL protein in monocots, highlighting the utility of the Brachypodium model for isolation and analysis of agronomically important genes for crop improvement.

Aegilops tauschii genome assembly Aet v5.0 features greater sequence contiguity and improved annotation

Aegilops tauschii is the donor of the D subgenome of hexaploid wheat and an important genetic resource. The reference-quality genome sequence Aet v4.0 for Ae. tauschii acc. AL8/78 was therefore an important milestone for wheat biology and breeding. Further advances in sequencing acc. AL8/78 and release of the Aet v5.0 sequence assembly are reported here. Two new optical maps were constructed and used in the revision of pseudomolecules. Gaps were closed with Pacific Biosciences long-read contigs, decreasing the gap number by 38,899. Transposable elements and protein-coding genes were reannotated. The number of annotated high-confidence genes was reduced from 39,635 in Aet v4.0 to 32,885 in Aet v5.0. A total of 2245 biologically important genes, including those affecting plant phenology, grain quality, and tolerance of abiotic stresses in wheat, was manually annotated and disease-resistance genes were annotated by a dedicated pipeline. Disease-resistance genes encoding nucleotide-binding site domains, receptor-like protein kinases, and receptor-like proteins were preferentially located in distal chromosome regions, whereas those encoding transmembrane coiled-coil proteins were dispersed more evenly along the chromosomes. Discovery, annotation, and expression analyses of microRNA (miRNA) precursors, mature miRNAs, and phasiRNAs are reported, including miRNA target genes. Other small RNAs, such as hc-siRNAs and tRFs, were characterized. These advances enhance the utility of the Ae. tauschii genome sequence for wheat genetics, biotechnology, and breeding.

Characterization of Pm68, a new powdery mildew resistance gene on chromosome 2BS of Greek durum wheat TRI 1796

New powdery mildew resistance gene Pm68 was found in the terminal region of chromosome 2BS of Greek durum wheat TRI 1796. The co-segregated molecular markers could be used for MAS. Durum wheat (Triticum turgidum L. var. durum Desf.) is not only an important cereal crop for pasta making, but also a genetic resource for common wheat improvement. In the present study, a Greek durum wheat TRI 1796 was found to confer high resistance to all 22 tested isolates of Blumeria graminis f. sp. tritici (Bgt). Inheritance study on the F₁ plants and the F₂ population derived from the cross TRI 1796/PI 584832 revealed that the resistance in TRI 1796 was controlled by a single dominant gene, herein designated Pm68. Using the bulked segregant RNA-Seq (BSR-Seq) analysis combined with molecular analysis, Pm68 was mapped to the terminal part of the short arm of chromosome 2B and flanked by markers Xdw04 and Xdw12/Xdw13 with genetic distances of 0.22 cM each. According to the reference genome of durum wheat cv. Svevo, the corresponding physical region spanned the Pm68 locus was about 1.78-Mb, in which a number of disease resistance-related genes were annotated. This study reports the new powdery mildew resistance gene Pm68 that would be a valuable resource for improvement of both common wheat and durum wheat. The co-segregated markers (Xdw05-Xdw11) developed here would be useful tools for marker-assisted selection (MAS) in breeding.

A rare single nucleotide variant in Pm5e confers powdery mildew resistance in common wheat

Powdery mildew poses severe threats to wheat production. The most sustainable way to control this disease is through planting resistant cultivars. We report the map-based cloning of the powdery mildew resistance allele Pm5e from a Chinese wheat landrace. We applied a two-step bulked segregant RNA sequencing (BSR-Seq) approach in developing tightly linked or co-segregating markers to Pm5e. The first BSR-Seq used phenotypically contrasting bulks of recombinant inbred lines (RILs) to identify Pm5e-linked markers. The second BSR-Seq utilized bulks of genetic recombinants screened from a fine-mapping population to precisely quantify the associated genomic variation in the mapping interval, and identified the Pm5e candidate genes. The function of Pm5e was validated by transgenic assay, loss-of-function mutants and haplotype association analysis. Pm5e encodes a nucleotide-binding domain leucine-rich-repeat-containing (NLR) protein. A rare nonsynonymous single nucleotide variant (SNV) within the C-terminal leucine rich repeat (LRR) domain is responsible for the gain of powdery mildew resistance function of Pm5e, an allele endemic to wheat landraces of Shaanxi province of China. Results from this study demonstrate the value of landraces in discovering useful genes for modern wheat breeding. The key SNV associated with powdery mildew resistance will be useful for marker-assisted selection of Pm5e in wheat breeding programs.

Exploiting the reference genome sequence of hexaploid wheat: a proteomic study of flour proteins from the cultivar Chinese Spring

Although the economic value of wheat flour is determined by the complement of gluten proteins, these proteins have been challenging to study because of the complexity of the major protein groups and the tremendous sequence diversity among wheat cultivars. The completion of a high-quality wheat genome sequence from the reference wheat Chinese Spring recently facilitated the assembly and annotation of a complete set of gluten protein genes from a single cultivar, making it possible to link individual proteins in the flour to specific gene sequences. In a proteomic analysis of total wheat flour protein from Chinese Spring using quantitative two-dimensional gel electrophoresis combined with tandem mass spectrometry, gliadins or low-molecular-weight glutenin subunits were identified as the predominant proteins in 72 protein spots. Individual spots were associated with 40 of 56 Chinese Spring gene sequences, including 16 of 26 alpha gliadins, 10 of 11 gamma gliadins, six of seven omega gliadins, one of two delta gliadins, and nine of ten LMW-GS. Most genes that were not associated with protein spots were either expressed at low levels in endosperm or encoded proteins with high similarity to other proteins. A wide range of protein accumulation levels were observed and discrepancies between transcript levels and protein levels were noted. This work together with similar studies using other commercial cultivars should provide new insight into the molecular basis of wheat flour quality and allergenic potential.

Rapid evolution of α-gliadin gene family revealed by analyzing Gli-2 locus regions of wild emmer wheat

α-Gliadins are a major group of gluten proteins in wheat flour that contribute to the end-use properties for food processing and contain major immunogenic epitopes that can cause serious health-related issues including celiac disease (CD). α-Gliadins are also the youngest group of gluten proteins and are encoded by a large gene family. The majority of the gene family members evolved independently in the A, B, and D genomes of different wheat species after their separation from a common ancestral species. To gain insights into the origin and evolution of these complex genes, the genomic regions of the Gli-2 loci encoding α-gliadins were characterized from the tetraploid wild emmer, a progenitor of hexaploid bread wheat that contributed the AABB genomes. Genomic sequences of Gli-2 locus regions for the wild emmer A and B genomes were first reconstructed using the genome sequence scaffolds along with optical genome maps. A total of 24 and 16 α-gliadin genes were identified for the A and B genome regions, respectively. α-Gliadin pseudogene frequencies of 86% for the A genome and 69% for the B genome were primarily caused by C to T substitutions in the highly abundant glutamine codons, resulting in the generation of premature stop codons. Comparison with the homologous regions from the hexaploid wheat cv. Chinese Spring indicated considerable sequence divergence of the two A genomes at the genomic level. In comparison, conserved regions between the two B genomes were identified that included α-gliadin pseudogenes containing shared nested TE insertions. Analyses of the genomic organization and phylogenetic tree reconstruction indicate that although orthologous gene pairs derived from speciation were present, large portions of α-gliadin genes were likely derived from differential gene duplications or deletions after the separation of the homologous wheat genomes ~ 0.5 MYA. The higher number of full-length intact α-gliadin genes in hexaploid wheat than that in wild emmer suggests that human selection through domestication might have an impact on α-gliadin evolution. Our study provides insights into the rapid and dynamic evolution of genomic regions harboring the α-gliadin genes in wheat.

Towards reducing the immunogenic potential of wheat flour: omega gliadins encoded by the D genome of hexaploid wheat may also harbor epitopes for the serious food allergy WDEIA

BACKGROUND

Omega-5 gliadins are a group of highly repetitive gluten proteins in wheat flour encoded on the 1B chromosome of hexaploid wheat. These proteins are the major sensitizing allergens in a severe form of food allergy called wheat-dependent exercise-induced anaphylaxis (WDEIA). The elimination of omega-5 gliadins from wheat flour through biotechnology or breeding approaches could reduce the immunogenic potential and adverse health effects of the flour.

RESULTS

A mutant line missing low-molecular weight glutenin subunits encoded at the Glu-B3 locus was selected previously from a doubled haploid population generated from two Korean wheat cultivars. Analysis of flour from the mutant line by 2-dimensional gel electrophoresis coupled with tandem mass spectrometry revealed that the omega-5 gliadins and several gamma gliadins encoded by the closely linked Gli-B1 locus were also missing as a result of a deletion of at least 5.8 Mb of chromosome 1B. Two-dimensional immunoblot analysis of flour proteins using sera from WDEIA patients showed reduced IgE reactivity in the mutant relative to the parental lines due to the absence of the major omega-5 gliadins. However, two minor proteins showed strong reactivity to patient sera in both the parental and the mutant lines and also reacted with a monoclonal antibody against omega-5 gliadin. Analysis of the two minor reactive proteins by mass spectrometry revealed that both proteins correspond to omega-5 gliadin genes encoded on chromosome 1D that were thought previously to be pseudogenes.

CONCLUSIONS

While breeding approaches can be used to reduce the levels of the highly immunogenic omega-5 gliadins in wheat flour, these approaches are complicated by the genetic linkage of different classes of gluten protein genes and the finding that omega-5 gliadins may be encoded on more than one chromosome. The work illustrates the importance of detailed knowledge about the genomic regions harboring the major gluten protein genes in individual wheat cultivars for future efforts aimed at reducing the immunogenic potential of wheat flour.

Analysis of the Gli-D2 locus identifies a genetic target for simultaneously improving the breadmaking and health-related traits of common wheat

Gliadins are a major component of wheat seed proteins. However, the complex homoeologous Gli-2 loci (Gli-A2, -B2 and -D2) that encode the α-gliadins in commercial wheat are still poorly understood. Here we analyzed the Gli-D2 locus of Xiaoyan 81 (Xy81), a winter wheat cultivar. A total of 421.091 kb of the Gli-D2 sequence was assembled from sequencing multiple bacterial artificial clones, and 10 α-gliadin genes were annotated. Comparative genomic analysis showed that Xy81 carried only eight of the α-gliadin genes of the D genome donor Aegilops tauschii, with two of them each experiencing a tandem duplication. A mutant line lacking Gli-D2 (DLGliD2) consistently exhibited better breadmaking quality and dough functionalities than its progenitor Xy81, but without penalties in other agronomic traits. It also had an elevated lysine content in the grains. Transcriptome analysis verified the lack of Gli-D2 α-gliadin gene expression in DLGliD2. Furthermore, the transcript and protein levels of protein disulfide isomerase were both upregulated in DLGliD2 grains. Consistent with this finding, DLGliD2 had increased disulfide content in the flour. Our work sheds light on the structure and function of Gli-D2 in commercial wheat, and suggests that the removal of Gli-D2 and the gliadins specified by it is likely to be useful for simultaneously enhancing the end-use and health-related traits of common wheat. Because gliadins and homologous proteins are widely present in grass species, the strategy and information reported here may be broadly useful for improving the quality traits of diverse cereal crops.

Modular assembly of transposable element arrays by microsatellite targeting in the guayule and rice genomes

Background

Guayule (Parthenium argentatum A. Gray) is a rubber-producing desert shrub native to Mexico and the United States. Guayule represents an alternative to Hevea brasiliensis as a source for commercial natural rubber. The efficient application of modern molecular/genetic tools to guayule improvement requires characterization of its genome.

Results

The 1.6 Gb guayule genome was sequenced, assembled and annotated. The final 1.5 Gb assembly, while fragmented (N₅₀ = 22 kb), maps > 95% of the shotgun reads and is essentially complete. Approximately 40,000 transcribed, protein encoding genes were annotated on the assembly. Further characterization of this genome revealed 15 families of small, microsatellite-associated, transposable elements (TEs) with unexpected chromosomal distribution profiles. These SaTar (Satellite Targeted) elements, which are non-autonomous Mu-like elements (MULEs), were frequently observed in multimeric linear arrays of unrelated individual elements within which no individual element is interrupted by another. This uniformly non-nested TE multimer architecture has not been previously described in either eukaryotic or prokaryotic genomes. Five families of similarly distributed non-autonomous MULEs (microsatellite associated, modularly assembled) were characterized in the rice genome. Families of TEs with similar structures and distribution profiles were identified in sorghum and citrus.

Conclusion

The sequencing and assembly of the guayule genome provides a foundation for application of current crop improvement technologies to this plant. In addition, characterization of this genome revealed SaTar elements with distribution profiles unique among TEs. Satar targeting appears based on an alternative MULE recombination mechanism with the potential to impact gene evolution.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4653-6) contains supplementary material, which is available to authorized users.

Dynamic Evolution of α-Gliadin Prolamin Gene Family in Homeologous Genomes of Hexaploid Wheat

Wheat Gli-2 loci encode complex groups of α-gliadin prolamins that are important for breadmaking, but also major triggers of celiac disease (CD). Elucidation of α-gliadin evolution provides knowledge to produce wheat with better end-use properties and reduced immunogenic potential. The Gli-2 loci contain a large number of tandemly duplicated genes and highly repetitive DNA, making sequence assembly of their genomic regions challenging. Here, we constructed high-quality sequences spanning the three wheat homeologous α-gliadin loci by aligning PacBio-based sequence contigs with BioNano genome maps. A total of 47 α-gliadin genes were identified with only 26 encoding intact full-length protein products. Analyses of α-gliadin loci and phylogenetic tree reconstruction indicate significant duplications of α-gliadin genes in the last ~2.5 million years after the divergence of the A, B and D genomes, supporting its rapid lineage-independent expansion in different Triticeae genomes. We showed that dramatic divergence in expression of α-gliadin genes could not be attributed to sequence variations in the promoter regions. The study also provided insights into the evolution of CD epitopes and identified a single indel event in the hexaploid wheat D genome that likely resulted in the generation of the highly toxic 33-mer CD epitope.

Gene Duplication and Evolution Dynamics in the Homeologous Regions Harboring Multiple Prolamin and Resistance Gene Families in Hexaploid Wheat

Improving end-use quality and disease resistance are important goals in wheat breeding. The genetic loci controlling these traits are highly complex, consisting of large families of prolamin and resistance genes with members present in all three homeologous A, B, and D genomes in hexaploid bread wheat. Here, orthologous regions harboring both prolamin and resistance gene loci were reconstructed and compared to understand gene duplication and evolution in different wheat genomes. Comparison of the two orthologous D regions from the hexaploid wheat Chinese Spring and the diploid progenitor Aegilops tauschii revealed their considerable difference due to the presence of five large structural variations with sizes ranging from 100 kb to 2 Mb. As a result, 44% of the Ae. tauschii and 71% of the Chinese Spring sequences in the analyzed regions, including 79 genes, are not shared. Gene rearrangement events, including differential gene duplication and deletion in the A, B, and D regions, have resulted in considerable erosion of gene collinearity in the analyzed regions, suggesting rapid evolution of prolamin and resistance gene families after the separation of the three wheat genomes. We hypothesize that this fast evolution is attributed to the co-evolution of the two gene families dispersed within a high recombination region. The identification of a full set of prolamin genes facilitated transcriptome profiling and revealed that the A genome contributes the least to prolamin expression because of its smaller number of expressed intact genes and their low expression levels, while the B and D genomes contribute similarly.

New insights into structural organization and gene duplication in a 1.75-Mb genomic region harboring the α-gliadin gene family in Aegilops tauschii, the source of wheat D genome

Among the wheat prolamins important for its end-use traits, α-gliadins are the most abundant, and are also a major cause of food-related allergies and intolerances. Previous studies of various wheat species estimated that between 25 and 150 α-gliadin genes reside in the Gli-2 locus regions. To better understand the evolution of this complex gene family, the DNA sequence of a 1.75-Mb genomic region spanning the Gli-2 locus was analyzed in the diploid grass, Aegilops tauschii, the ancestral source of D genome in hexaploid bread wheat. Comparison with orthologous regions from rice, sorghum, and Brachypodium revealed rapid and dynamic changes only occurring to the Ae. tauschii Gli-2 region, including insertions of high numbers of non-syntenic genes and a high rate of tandem gene duplications, the latter of which have given rise to 12 copies of α-gliadin genes clustered within a 550-kb region. Among them, five copies have undergone pseudogenization by various mutation events. Insights into the evolutionary relationship of the duplicated α-gliadin genes were obtained from their genomic organization, transcription patterns, transposable element insertions and phylogenetic analyses. An ancestral glutamate-like receptor (GLR) gene encoding putative amino acid sensor in all four grass species has duplicated only in Ae. tauschii and generated three more copies that are interspersed with the α-gliadin genes. Phylogenetic inference and different gene expression patterns support functional divergence of the Ae. tauschii GLR copies after duplication. Our results suggest that the duplicates of α-gliadin and GLR genes have likely taken different evolutionary paths; conservation for the former and neofunctionalization for the latter.

Sequencing and comparative analyses of Aegilops tauschii chromosome arm 3DS reveal rapid evolution of Triticeae genomes

Bread wheat (Triticum aestivum, AABBDD) is an allohexaploid species derived from two rounds of interspecific hybridizations. A high-quality genome sequence assembly of diploid Aegilops tauschii, the donor of the wheat D genome, will provide a useful platform to study polyploid wheat evolution. A combined approach of BAC pooling and next-generation sequencing technology was employed to sequence the minimum tiling path (MTP) of 3176 BAC clones from the short arm of Ae. tauschii chromosome 3 (At3DS). The final assembly of 135 super-scaffolds with an N50 of 4.2 Mb was used to build a 247-Mb pseudomolecule with a total of 2222 predicted protein-coding genes. Compared with the orthologous regions of rice, Brachypodium, and sorghum, At3DS contains 38.67% more genes. In comparison to At3DS, the short arm sequence of wheat chromosome 3B (Ta3BS) is 95-Mb large in size, which is primarily due to the expansion of the non-centromeric region, suggesting that transposable element (TE) bursts in Ta3B likely occurred there. Also, the size increase is accompanied by a proportional increase in gene number in Ta3BS. We found that in the sequence of short arm of wheat chromosome 3D (Ta3DS), there was only less than 0.27% gene loss compared to At3DS. Our study reveals divergent evolution of grass genomes and provides new insights into sequence changes in the polyploid wheat genome.

Rapid evolutionary dynamics in a 2.8-Mb chromosomal region containing multiple prolamin and resistance gene families in Aegilops tauschii

Prolamin and resistance gene families are important in wheat food use and in defense against pathogen attacks, respectively. To better understand the evolution of these multi-gene families, the DNA sequence of a 2.8-Mb genomic region, representing an 8.8 cM genetic interval and harboring multiple prolamin and resistance-like gene families, was analyzed in the diploid grass Aegilops tauschii, the D-genome donor of bread wheat. Comparison with orthologous regions from rice, Brachypodium, and sorghum showed that the Ae. tauschii region has undergone dramatic changes; it has acquired more than 80 non-syntenic genes and only 13 ancestral genes are shared among these grass species. These non-syntenic genes, including prolamin and resistance-like genes, originated from various genomic regions and likely moved to their present locations via sequence evolution processes involving gene duplication and translocation. Local duplication of non-syntenic genes contributed significantly to the expansion of gene families. Our analysis indicates that the insertion of prolamin-related genes occurred prior to the separation of the Brachypodieae and Triticeae lineages. Unlike in Brachypodium, inserted prolamin genes have rapidly evolved and expanded to encode different classes of major seed storage proteins in Triticeae species. Phylogenetic analyses also showed that the multiple insertions of resistance-like genes and subsequent differential expansion of each R gene family. The high frequency of non-syntenic genes and rapid local gene evolution correlate with the high recombination rate in the 2.8-Mb region with nine-fold higher than the genome-wide average. Our results demonstrate complex evolutionary dynamics in this agronomically important region of Triticeae species.

GSP: a web-based platform for designing genome-specific primers in polyploids

MOTIVATION

The sequences among subgenomes in a polyploid species have high similarity, making it difficult to design genome-specific primers for sequence analysis.

RESULTS

We present GSP, a web-based platform to design genome-specific primers that distinguish subgenome sequences in a polyploid genome. GSP uses BLAST to extract homeologous sequences of the subgenomes in existing databases, performs a multiple sequence alignment, and design primers based on sequence variants in the alignment. An interactive primers diagram, a sequence alignment viewer and a virtual electrophoresis are displayed as parts of the primer design result. GSP also designs specific primers from multiple sequences uploaded by users.

AVAILABILITY AND IMPLEMENTATION

GSP is a user-friendly and efficient web platform freely accessible at http://probes.pw.usda.gov/GSP Source code and command-line application are available at https://github.com/bioinfogenome/GSP CONTACTS: yong.gu@ars.usda.gov or devin.coleman-derr@ars.usda.gov

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Development of a D genome specific marker resource for diploid and hexaploid wheat

BACKGROUND

Mapping and map-based cloning of genes that control agriculturally and economically important traits remain great challenges for plants with complex highly repetitive genomes such as those within the grass tribe, Triticeae. Mapping limitations in the Triticeae are primarily due to low frequencies of polymorphic gene markers and poor genetic recombination in certain genetic regions. Although the abundance of repetitive sequence may pose common problems in genome analysis and sequence assembly of large and complex genomes, they provide repeat junction markers with random and unbiased distribution throughout chromosomes. Hence, development of a high-throughput mapping technology that combine both gene-based and repeat junction-based markers is needed to generate maps that have better coverage of the entire genome.

RESULTS

In this study, the available genomics resource of the diploid Aegilop tauschii, the D genome donor of bread wheat, were used to develop genome specific markers that can be applied for mapping in modern hexaploid wheat. A NimbleGen array containing both gene-based and repeat junction probe sequences derived from Ae. tauschii was developed and used to map the Chinese Spring nullisomic-tetrasomic lines and deletion bin lines of the D genome chromosomes. Based on these mapping data, we have now anchored 5,171 repeat junction probes and 10,892 gene probes, corresponding to 5,070 gene markers, to the delineated deletion bins of the D genome. The order of the gene-based markers within the deletion bins of the Chinese Spring can be inferred based on their positions on the Ae. tauschii genetic map. Analysis of the probe sequences against the Chinese Spring chromosome sequence assembly database facilitated mapping of the NimbleGen probes to the sequence contigs and allowed assignment or ordering of these sequence contigs within the deletion bins. The accumulated length of anchored sequence contigs is about 155 Mb, representing ~ 3.2 % of the D genome. A specific database was developed to allow user to search or BLAST against the probe sequence information and to directly download PCR primers for mapping specific genetic loci.

CONCLUSIONS

In bread wheat, aneuploid stocks have been extensively used to assign markers linked with genes/traits to chromosomes, chromosome arms, and their specific bins. Through this study, we added thousands of markers to the existing wheat chromosome bin map, representing a significant step forward in providing a resource to navigate the wheat genome. The database website ( http://probes.pw.usda.gov/ATRJM/ ) provides easy access and efficient utilization of the data. The resources developed herein can aid map-based cloning of traits of interest and the sequencing of the D genome of hexaploid wheat.

Dynamic evolution of resistance gene analogs in the orthologous genomic regions of powdery mildew resistance gene MlIW170 in Triticum dicoccoides and Aegilops tauschii

Rapid evolution of powdery mildew resistance gene MlIW170 orthologous genomic regions in wheat subgenomes. Wheat is one of the most important staple grain crops in the world and also an excellent model for plant ploidy evolution research with different ploidy levels from diploid to hexaploid. Powdery mildew disease caused by Blumeria graminis f.sp. tritici can result in significant loss in both grain yield and quality in wheat. In this study, the wheat powdery mildew resistance gene MlIW170 locus located at the Triticum dicoccoides chromosome 2B short arm was further characterized by constructing and sequencing a BAC-based physical map contig covering a 0.3 cM genetic distance region (880 kb) and developing additional markers to delineate the resistance gene within a 0.16 cM genetic interval (372 kb). Comparative analyses of the T. dicoccoides 2BS region with the orthologous Aegilops tauschii 2DS region showed great gene colinearity, including the structure organization of both types of RGA1/2-like and RPS2-like resistance genes. Comparative analyses with the orthologous regions from Brachypodium and rice genomes revealed considerable dynamic evolutionary changes that have re-shaped this MlIW170 region in the wheat genome, resulting in a high number of non-syntenic genes including resistance-related genes. This result might reflect the rapid evolution in R-gene regions. Phylogenetic analysis on these resistance-related gene sequences indicated the duplication of these genes in the MlIW170 region, occurred before the separation of the wheat B and D genomes.

Fine physical and genetic mapping of powdery mildew resistance gene MlIW172 originating from wild emmer (Triticum dicoccoides)

Powdery mildew, caused by Blumeria graminis f. sp. tritici, is one of the most important wheat diseases in the world. In this study, a single dominant powdery mildew resistance gene MlIW172 was identified in the IW172 wild emmer accession and mapped to the distal region of chromosome arm 7AL (bin7AL-16-0.86-0.90) via molecular marker analysis. MlIW172 was closely linked with the RFLP probe Xpsr680-derived STS marker Xmag2185 and the EST markers BE405531 and BE637476. This suggested that MlIW172 might be allelic to the Pm1 locus or a new locus closely linked to Pm1. By screening genomic BAC library of durum wheat cv. Langdon and 7AL-specific BAC library of hexaploid wheat cv. Chinese Spring, and after analyzing genome scaffolds of Triticum urartu containing the marker sequences, additional markers were developed to construct a fine genetic linkage map on the MlIW172 locus region and to delineate the resistance gene within a 0.48 cM interval. Comparative genetics analyses using ESTs and RFLP probe sequences flanking the MlIW172 region against other grass species revealed a general co-linearity in this region with the orthologous genomic regions of rice chromosome 6, Brachypodium chromosome 1, and sorghum chromosome 10. However, orthologous resistance gene-like RGA sequences were only present in wheat and Brachypodium. The BAC contigs and sequence scaffolds that we have developed provide a framework for the physical mapping and map-based cloning of MlIW172.

Accelerated evolution of the mitochondrial genome in an alloplasmic line of durum wheat

BACKGROUND

Wheat is an excellent plant species for nuclear mitochondrial interaction studies due to availability of large collection of alloplasmic lines. These lines exhibit different vegetative and physiological properties than their parents. To investigate the level of sequence changes introduced into the mitochondrial genome under the alloplasmic condition, three mitochondrial genomes of the Triticum-Aegilops species were sequenced: 1) durum alloplasmic line with the Ae. longissima cytoplasm that carries the T. turgidum nucleus designated as (lo) durum, 2) the cytoplasmic donor line, and 3) the nuclear donor line.

RESULTS

The mitochondrial genome of the T. turgidum was 451,678 bp in length with high structural and nucleotide identity to the previously characterized T. aestivum genome. The assembled mitochondrial genome of the (lo) durum and the Ae. longissima were 431,959 bp and 399,005 bp in size, respectively. The high sequence coverage for all three genomes allowed analysis of heteroplasmy within each genome. The mitochondrial genome structure in the alloplasmic line was genetically distant from both maternal and paternal genomes. The alloplasmic durum and the Ae. longissima carry the same versions of atp6, nad6, rps19-p, cob and cox2 exon 2 which are different from the T. turgidum parent. Evidence of paternal leakage was also observed by analyzing nad9 and orf359 among all three lines. Nucleotide search identified a number of open reading frames, of which 27 were specific to the (lo) durum line.

CONCLUSIONS

Several heteroplasmic regions were observed within genes and intergenic regions of the mitochondrial genomes of all three lines. The number of rearrangements and nucleotide changes in the mitochondrial genome of the alloplasmic line that have occurred in less than half a century was significant considering the high sequence conservation between the T. turgidum and the T. aestivum that diverged from each other 10,000 years ago. We showed that the changes in genes were not limited to paternal leakage but were sufficiently significant to suggest that other mechanisms, such as recombination and mutation, were responsible. The newly formed ORFs, differences in gene sequences and copy numbers, heteroplasmy, and substoichiometric changes show the potential of the alloplasmic condition to accelerate evolution towards forming new mitochondrial genomes.

Structural analysis of the wheat genes encoding NADH-dependent glutamine-2-oxoglutarate amidotransferases and correlation with grain protein content

BACKGROUND

Nitrogen uptake and the efficient absorption and metabolism of nitrogen are essential elements in attempts to breed improved cereal cultivars for grain or silage production. One of the enzymes related to nitrogen metabolism is glutamine-2-oxoglutarate amidotransferase (GOGAT). Together with glutamine synthetase (GS), GOGAT maintains the flow of nitrogen from NH4 (+) into glutamine and glutamate, which are then used for several aminotransferase reactions during amino acid synthesis.

RESULTS

The aim of the present work was to identify and analyse the structure of wheat NADH-GOGAT genomic sequences, and study the expression in two durum wheat cultivars characterized by low and high kernel protein content. The genomic sequences of the three homoeologous A, B and D NADH-GOGAT genes were obtained for hexaploid Triticum aestivum and the tetraploid A and B genes of Triticum turgidum ssp. durum. Analysis of the gene sequences indicates that all wheat NADH-GOGAT genes are composed of 22 exons and 21 introns. The three hexaploid wheat homoeologous genes have high conservation of sequence except intron 13 which shows differences in both length and sequence. A comparative analysis of sequences among di- and mono-cotyledonous plants shows both regions of high conservation and of divergence. qRT-PCR performed with the two durum wheat cvs Svevo and Ciccio (characterized by high and low protein content, respectively) indicates different expression levels of the two NADH-GOGAT-3A and NADH-GOGAT-3B genes.

CONCLUSION

The three hexaploid wheat homoeologous NADH-GOGAT gene sequences are highly conserved - consistent with the key metabolic role of this gene. However, the dicot and monocot amino acid sequences show distinctive patterns, particularly in the transit peptide, the exon 16-17 junction, and the C-terminus. The lack of conservation in the transit peptide may indicate subcellular differences between the two plant divisions - while the sequence conservation within enzyme functional domains remains high. Higher expression levels of NADH-GOGAT are associated with higher grain protein content in two durum wheats.

The gene space in wheat: the complete γ-gliadin gene family from the wheat cultivar Chinese Spring

The complete set of unique γ-gliadin genes is described for the wheat cultivar Chinese Spring using a combination of expressed sequence tag (EST) and Roche 454 DNA sequences. Assemblies of Chinese Spring ESTs yielded 11 different γ-gliadin gene sequences. Two of the sequences encode identical polypeptides and are assumed to be the result of a recent gene duplication. One gene has a 3' coding mutation that changes the reading frame in the final eight codons. A second assembly of Chinese Spring γ-gliadin sequences was generated using Roche 454 total genomic DNA sequences. The 454 assembly confirmed the same 11 active genes as the EST assembly plus two pseudogenes not represented by ESTs. These 13 γ-gliadin sequences represent the complete unique set of γ-gliadin genes for cv Chinese Spring, although not ruled out are additional genes that are exact duplications of these 13 genes. A comparison with the ESTs of two other hexaploid cultivars (Butte 86 and Recital) finds that the most active genes are present in all three cultivars, with exceptions likely due to too few ESTs for detection in Butte 86 and Recital. A comparison of the numbers of ESTs per gene indicates differential levels of expression within the γ-gliadin gene family. Genome assignments were made for 6 of the 13 Chinese Spring γ-gliadin genes, i.e., one assignment from a match to two γ-gliadin genes found within a tetraploid wheat A genome BAC and four genes that match four distinct γ-gliadin sequences assembled from Roche 454 sequences from Aegilops tauschii, the hexaploid wheat D-genome ancestor.

Association of suppressor of cytokine signalling 3 polymorphisms with insulin resistance in patients with chronic hepatitis C

Suppressor of cytokine signalling 3 is thought to be associated with insulin resistance in patients with chronic hepatitis C. We evaluated the role of suppressor of cytokine signalling 3 polymorphisms in determining insulin resistance in patients with chronic hepatitis C. Two hundred and ninety untreated hepatitis C virus-infected patients without diabetes and cirrhosis were genotyped for the SNPs rs4969168, rs4969170 and rs12952093 of suppressor of cytokine signalling 3 using the TaqMan Genotyping Assay. We found that the rs4969170 AA genotype and rs4969170 A allele frequency were significantly more common in the insulin-resistant group than the non-insulin-resistant group (89.5% vs 76.1%, OR = 2.693, 95% CI: 1.221-5.939, P = 0.012 and 94.8% vs 88.0%, OR = 2.463, 95% CI: 1.151-5.271, P = 0.017, respectively). Haplotype G-C was likely associated with non-insulin resistance (adjusted P = 0.011). Multiple logistic regression analysis indicates that the independent risk factors for insulin resistance are the SNP rs4969170 AA genotype (OR = 3.005, 95% CI: 1.194-7.560, P = 0.019), HCV genotype 1 (OR = 2.524, 95% CI: 1.099-5.794, P = 0.029) and BMI (OR = 0.514, 95% CI: 0.265-0.999, P = 0.05).

Rapid genome mapping in nanochannel arrays for highly complete and accurate de novo sequence assembly of the complex Aegilops tauschii genome

Next-generation sequencing (NGS) technologies have enabled high-throughput and low-cost generation of sequence data; however, de novo genome assembly remains a great challenge, particularly for large genomes. NGS short reads are often insufficient to create large contigs that span repeat sequences and to facilitate unambiguous assembly. Plant genomes are notorious for containing high quantities of repetitive elements, which combined with huge genome sizes, makes accurate assembly of these large and complex genomes intractable thus far. Using two-color genome mapping of tiling bacterial artificial chromosomes (BAC) clones on nanochannel arrays, we completed high-confidence assembly of a 2.1-Mb, highly repetitive region in the large and complex genome of Aegilops tauschii, the D-genome donor of hexaploid wheat (Triticum aestivum). Genome mapping is based on direct visualization of sequence motifs on single DNA molecules hundreds of kilobases in length. With the genome map as a scaffold, we anchored unplaced sequence contigs, validated the initial draft assembly, and resolved instances of misassembly, some involving contigs <2 kb long, to dramatically improve the assembly from 75% to 95% complete.

A new class of wheat gliadin genes and proteins

The utility of mining DNA sequence data to understand the structure and expression of cereal prolamin genes is demonstrated by the identification of a new class of wheat prolamins. This previously unrecognized wheat prolamin class, given the name δ-gliadins, is the most direct ortholog of barley γ3-hordeins. Phylogenetic analysis shows that the orthologous δ-gliadins and γ3-hordeins form a distinct prolamin branch that existed separate from the γ-gliadins and γ-hordeins in an ancestral Triticeae prior to the branching of wheat and barley. The expressed δ-gliadins are encoded by a single gene in each of the hexaploid wheat genomes. This single δ-gliadin/γ3-hordein ortholog may be a general feature of the Triticeae tribe since examination of ESTs from three barley cultivars also confirms a single γ3-hordein gene. Analysis of ESTs and cDNAs shows that the genes are expressed in at least five hexaploid wheat cultivars in addition to diploids Triticum monococcum and Aegilops tauschii. The latter two sequences also allow assignment of the δ-gliadin genes to the A and D genomes, respectively, with the third sequence type assumed to be from the B genome. Two wheat cultivars for which there are sufficient ESTs show different patterns of expression, i.e., with cv Chinese Spring expressing the genes from the A and B genomes, while cv Recital has ESTs from the A and D genomes. Genomic sequences of Chinese Spring show that the D genome gene is inactivated by tandem premature stop codons. A fourth δ-gliadin sequence occurs in the D genome of both Chinese Spring and Ae. tauschii, but no ESTs match this sequence and limited genomic sequences indicates a pseudogene containing frame shifts and premature stop codons. Sequencing of BACs covering a 3 Mb region from Ae. tauschii locates the δ-gliadin gene to the complex Gli-1 plus Glu-3 region on chromosome 1.

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.

Fine mapping of the Bsr1 barley stripe mosaic virus resistance gene in the model grass Brachypodium distachyon

The ND18 strain of Barley stripe mosaic virus (BSMV) infects several lines of Brachypodium distachyon, a recently developed model system for genomics research in cereals. Among the inbred lines tested, Bd3-1 is highly resistant at 20 to 25°C, whereas Bd21 is susceptible and infection results in an intense mosaic phenotype accompanied by high levels of replicating virus. We generated an F_6∶7 recombinant inbred line (RIL) population from a cross between Bd3-1 and Bd21 and used the RILs, and an F₂ population of a second Bd21 × Bd3-1 cross to evaluate the inheritance of resistance. The results indicate that resistance segregates as expected for a single dominant gene, which we have designated Barley stripe mosaic virus resistance 1 (Bsr1). We constructed a genetic linkage map of the RIL population using SNP markers to map this gene to within 705 Kb of the distal end of the top of chromosome 3. Additional CAPS and Indel markers were used to fine map Bsr1 to a 23 Kb interval containing five putative genes. Our study demonstrates the power of using RILs to rapidly map the genetic determinants of BSMV resistance in Brachypodium. Moreover, the RILs and their associated genetic map, when combined with the complete genomic sequence of Brachypodium, provide new resources for genetic analyses of many other traits.

Characterizing the citrus cultivar Carrizo genome through 454 shotgun sequencing

The citrus cultivar Carrizo is the single most important rootstock to the US citrus industry and has resistance or tolerance to a number of major citrus diseases, including citrus tristeza virus, foot rot, and Huanglongbing (HLB, citrus greening). A Carrizo genomic sequence database providing approximately 3.5× genome coverage (haploid genome size approximately 367 Mb) was populated through 454 GS FLX shotgun sequencing. Analysis of the repetitive DNA fraction indicated a total interspersed repeat fraction of 36.5%. Assembly and characterization of abundant citrus Ty3/gypsy elements revealed a novel type of element containing open reading frames encoding a viral RNA-silencing suppressor protein (RNA binding protein, rbp) and a plant cytokinin riboside 5′-monophosphate phosphoribohydrolase-related protein (LONELY GUY, log). Similar gypsy elements were identified in the Populus trichocarpa genome. Gene-coding region analysis indicated that 24.4% of the nonrepetitive reads contained genic regions. The depth of genome coverage was sufficient to allow accurate assembly of constituent genes, including a putative phloem-expressed gene. The development of the Carrizo database ( http://citrus.pw.usda.gov/ ) will contribute to characterization of agronomically significant loci and provide a publicly available genomic resource to the citrus research community.

Comparison of a high-density genetic linkage map to genome features in the model grass Brachypodium distachyon

The small annual grass Brachypodium distachyon (Brachypodium) is rapidly emerging as a powerful model system to study questions unique to the grasses. Many Brachypodium resources have been developed including a whole genome sequence, highly efficient transformation and a large germplasm collection. We developed a genetic linkage map of Brachypodium using single nucleotide polymorphism (SNP) markers and an F(2) mapping population of 476 individuals. SNPs were identified by targeted resequencing of single copy genomic sequences. Using the Illumina GoldenGate Genotyping platform we placed 558 markers into five linkage groups corresponding to the five chromosomes of Brachypodium. The unusually long total genetic map length, 1,598 centiMorgans (cM), indicates that the Brachypodium mapping population has a high recombination rate. By comparing the genetic map to genome features we found that the recombination rate was positively correlated with gene density and negatively correlated with repetitive regions and sites of ancestral chromosome fusions that retained centromeric repeat sequences. A comparison of adjacent genome regions with high versus low recombination rates revealed a positive correlation between interspecific synteny and recombination rate.

Annotation-based genome-wide SNP discovery in the large and complex Aegilops tauschii genome using next-generation sequencing without a reference genome sequence

BACKGROUND

Many plants have large and complex genomes with an abundance of repeated sequences. Many plants are also polyploid. Both of these attributes typify the genome architecture in the tribe Triticeae, whose members include economically important wheat, rye and barley. Large genome sizes, an abundance of repeated sequences, and polyploidy present challenges to genome-wide SNP discovery using next-generation sequencing (NGS) of total genomic DNA by making alignment and clustering of short reads generated by the NGS platforms difficult, particularly in the absence of a reference genome sequence.

RESULTS

An annotation-based, genome-wide SNP discovery pipeline is reported using NGS data for large and complex genomes without a reference genome sequence. Roche 454 shotgun reads with low genome coverage of one genotype are annotated in order to distinguish single-copy sequences and repeat junctions from repetitive sequences and sequences shared by paralogous genes. Multiple genome equivalents of shotgun reads of another genotype generated with SOLiD or Solexa are then mapped to the annotated Roche 454 reads to identify putative SNPs. A pipeline program package, AGSNP, was developed and used for genome-wide SNP discovery in Aegilops tauschii-the diploid source of the wheat D genome, and with a genome size of 4.02 Gb, of which 90% is repetitive sequences. Genomic DNA of Ae. tauschii accession AL8/78 was sequenced with the Roche 454 NGS platform. Genomic DNA and cDNA of Ae. tauschii accession AS75 was sequenced primarily with SOLiD, although some Solexa and Roche 454 genomic sequences were also generated. A total of 195,631 putative SNPs were discovered in gene sequences, 155,580 putative SNPs were discovered in uncharacterized single-copy regions, and another 145,907 putative SNPs were discovered in repeat junctions. These SNPs were dispersed across the entire Ae. tauschii genome. To assess the false positive SNP discovery rate, DNA containing putative SNPs was amplified by PCR from AL8/78 and AS75 and resequenced with the ABI 3730 xl. In a sample of 302 randomly selected putative SNPs, 84.0% in gene regions, 88.0% in repeat junctions, and 81.3% in uncharacterized regions were validated.

CONCLUSION

An annotation-based genome-wide SNP discovery pipeline for NGS platforms was developed. The pipeline is suitable for SNP discovery in genomic libraries of complex genomes and does not require a reference genome sequence. The pipeline is applicable to all current NGS platforms, provided that at least one such platform generates relatively long reads. The pipeline package, AGSNP, and the discovered 497,118 Ae. tauschii SNPs can be accessed at (http://avena.pw.usda.gov/wheatD/agsnp.shtml).

Nucleotide diversity maps reveal variation in diversity among wheat genomes and chromosomes

BACKGROUND

A genome-wide assessment of nucleotide diversity in a polyploid species must minimize the inclusion of homoeologous sequences into diversity estimates and reliably allocate individual haplotypes into their respective genomes. The same requirements complicate the development and deployment of single nucleotide polymorphism (SNP) markers in polyploid species. We report here a strategy that satisfies these requirements and deploy it in the sequencing of genes in cultivated hexaploid wheat (Triticum aestivum, genomes AABBDD) and wild tetraploid wheat (Triticum turgidum ssp. dicoccoides, genomes AABB) from the putative site of wheat domestication in Turkey. Data are used to assess the distribution of diversity among and within wheat genomes and to develop a panel of SNP markers for polyploid wheat.

RESULTS

Nucleotide diversity was estimated in 2114 wheat genes and was similar between the A and B genomes and reduced in the D genome. Within a genome, diversity was diminished on some chromosomes. Low diversity was always accompanied by an excess of rare alleles. A total of 5,471 SNPs was discovered in 1791 wheat genes. Totals of 1,271, 1,218, and 2,203 SNPs were discovered in 488, 463, and 641 genes of wheat putative diploid ancestors, T. urartu, Aegilops speltoides, and Ae. tauschii, respectively. A public database containing genome-specific primers, SNPs, and other information was constructed. A total of 987 genes with nucleotide diversity estimated in one or more of the wheat genomes was placed on an Ae. tauschii genetic map, and the map was superimposed on wheat deletion-bin maps. The agreement between the maps was assessed.

CONCLUSIONS

In a young polyploid, exemplified by T. aestivum, ancestral species are the primary source of genetic diversity. Low effective recombination due to self-pollination and a genetic mechanism precluding homoeologous chromosome pairing during polyploid meiosis can lead to the loss of diversity from large chromosomal regions. The net effect of these factors in T. aestivum is large variation in diversity among genomes and chromosomes, which impacts the development of SNP markers and their practical utility. Accumulation of new mutations in older polyploid species, such as wild emmer, results in increased diversity and its more uniform distribution across the genome.

Nucleotide diversity maps reveal variation in diversity among wheat genomes and chromosomes

BACKGROUND

A genome-wide assessment of nucleotide diversity in a polyploid species must minimize the inclusion of homoeologous sequences into diversity estimates and reliably allocate individual haplotypes into their respective genomes. The same requirements complicate the development and deployment of single nucleotide polymorphism (SNP) markers in polyploid species. We report here a strategy that satisfies these requirements and deploy it in the sequencing of genes in cultivated hexaploid wheat (Triticum aestivum, genomes AABBDD) and wild tetraploid wheat (Triticum turgidum ssp. dicoccoides, genomes AABB) from the putative site of wheat domestication in Turkey. Data are used to assess the distribution of diversity among and within wheat genomes and to develop a panel of SNP markers for polyploid wheat.

RESULTS

Nucleotide diversity was estimated in 2114 wheat genes and was similar between the A and B genomes and reduced in the D genome. Within a genome, diversity was diminished on some chromosomes. Low diversity was always accompanied by an excess of rare alleles. A total of 5,471 SNPs was discovered in 1791 wheat genes. Totals of 1,271, 1,218, and 2,203 SNPs were discovered in 488, 463, and 641 genes of wheat putative diploid ancestors, T. urartu, Aegilops speltoides, and Ae. tauschii, respectively. A public database containing genome-specific primers, SNPs, and other information was constructed. A total of 987 genes with nucleotide diversity estimated in one or more of the wheat genomes was placed on an Ae. tauschii genetic map, and the map was superimposed on wheat deletion-bin maps. The agreement between the maps was assessed.

CONCLUSIONS

In a young polyploid, exemplified by T. aestivum, ancestral species are the primary source of genetic diversity. Low effective recombination due to self-pollination and a genetic mechanism precluding homoeologous chromosome pairing during polyploid meiosis can lead to the loss of diversity from large chromosomal regions. The net effect of these factors in T. aestivum is large variation in diversity among genomes and chromosomes, which impacts the development of SNP markers and their practical utility. Accumulation of new mutations in older polyploid species, such as wild emmer, results in increased diversity and its more uniform distribution across the genome.

RJPrimers: unique transposable element insertion junction discovery and PCR primer design for marker development

Transposable elements (TE) exist in the genomes of nearly all eukaryotes. TE mobilization through ‘cut-and-paste’ or ‘copy-and-paste’ mechanisms causes their insertions into other repetitive sequences, gene loci and other DNA. An insertion of a TE commonly creates a unique TE junction in the genome. TE junctions are also randomly distributed along chromosomes and therefore useful for genome-wide marker development. Several TE-based marker systems have been developed and applied to genetic diversity assays, and to genetic and physical mapping. A software tool ‘RJPrimers’ reported here allows for accurate identification of unique repeat junctions using BLASTN against annotated repeat databases and a repeat junction finding algorithm, and then for fully automated high-throughput repeat junction-based primer design using Primer3 and BatchPrimer3. The software was tested using the rice genome and genomic sequences of Aegilops tauschii. Over 90% of repeat junction primers designed by RJPrimers were unique. At least one RJM marker per 10 Kb sequence of A. tauschii was expected with an estimate of over 0.45 million such markers in a genome of 4.02 Gb, providing an almost unlimited source of molecular markers for mapping large and complex genomes. A web-based server and a command line-based pipeline for RJPrimers are both available at http://wheat.pw.usda.gov/demos/RJPrimers/.

An SSR-based genetic linkage map of the model grass Brachypodium distachyon

The grass species Brachypodium distachyon (hereafter, Brachypodium) has been adopted as a model system for grasses. Here, we describe the development of a genetic linkage map of Brachypodium. The genetic linkage map was developed with an F2 population from a cross between the diploid Brachypodium lines Bd3-1 and Bd21. The map was populated with polymorphic simple sequence repeat (SSR) markers from Brachypodium expressed sequence tag (EST) and bacterial artificial chromosome (BAC) end sequences and conserved orthologous sequence (COS) markers from other grass species. The map is 1386 cM in length and consists of 139 marker loci distributed across 20 linkage groups. Five of the linkage groups exceed 100 cM in length, with the largest being 231 cM long. Assessment of colinearity between the Brachypodium linkage map and the rice genome sequence revealed significant regions of macrosynteny between the two genomes, as well as rearrangements similar to those reported in other grass comparative structural genomics studies. The Brachypodium genetic linkage map described here will serve as a new tool to pursue a range of molecular genetic analyses and other applications in this new model plant system.

A BAC-based physical map of Brachypodium distachyon and its comparative analysis with rice and wheat

Background

Brachypodium distachyon (Brachypodium) has been recognized as a new model species for comparative and functional genomics of cereal and bioenergy crops because it possesses many biological attributes desirable in a model, such as a small genome size, short stature, self-pollinating habit, and short generation cycle. To maximize the utility of Brachypodium as a model for basic and applied research it is necessary to develop genomic resources for it. A BAC-based physical map is one of them. A physical map will facilitate analysis of genome structure, comparative genomics, and assembly of the entire genome sequence.

Results

A total of 67,151 Brachypodium BAC clones were fingerprinted with the SNaPshot HICF fingerprinting method and a genome-wide physical map of the Brachypodium genome was constructed. The map consisted of 671 contigs and 2,161 clones remained as singletons. The contigs and singletons spanned 414 Mb. A total of 13,970 gene-related sequences were detected in the BAC end sequences (BES). These gene tags aligned 345 contigs with 336 Mb of rice genome sequence, showing that Brachypodium and rice genomes are generally highly colinear. Divergent regions were mainly in the rice centromeric regions. A dot-plot of Brachypodium contigs against the rice genome sequences revealed remnants of the whole-genome duplication caused by paleotetraploidy, which were previously found in rice and sorghum. Brachypodium contigs were anchored to the wheat deletion bin maps with the BES gene-tags, opening the door to Brachypodium-Triticeae comparative genomics.

Conclusion

The construction of the Brachypodium physical map, and its comparison with the rice genome sequence demonstrated the utility of the SNaPshot-HICF method in the construction of BAC-based physical maps. The map represents an important genomic resource for the completion of Brachypodium genome sequence and grass comparative genomics. A draft of the physical map and its comparisons with rice and wheat are available at .

ConservedPrimers 2.0: a high-throughput pipeline for comparative genome referenced intron-flanking PCR primer design and its application in wheat SNP discovery

Background

In some genomic applications it is necessary to design large numbers of PCR primers in exons flanking one or several introns on the basis of orthologous gene sequences in related species. The primer pairs designed by this target gene approach are called "intron-flanking primers" or because they are located in exonic sequences which are usually conserved between related species, "conserved primers". They are useful for large-scale single nucleotide polymorphism (SNP) discovery and marker development, especially in species, such as wheat, for which a large number of ESTs are available but for which genome sequences and intron/exon boundaries are not available. To date, no suitable high-throughput tool is available for this purpose.

Results

We have developed, the ConservedPrimers 2.0 pipeline, for designing intron-flanking primers for large-scale SNP discovery and marker development, and demonstrated its utility in wheat. This tool uses non-redundant wheat EST sequences, such as wheat contigs and singleton ESTs, and related genomic sequences, such as those of rice, as inputs. It aligns the ESTs to the genomic sequences to identify unique colinear exon blocks and predicts intron lengths. Intron-flanking primers are then designed based on the intron/exon information using the Primer3 core program or BatchPrimer3. Finally, a tab-delimited file containing intron-flanking primer pair sequences and their primer properties is generated for primer ordering and their PCR applications. Using this tool, 1,922 bin-mapped wheat ESTs (31.8% of the 6,045 in total) were found to have unique colinear exon blocks suitable for primer design and 1,821 primer pairs were designed from these single- or low-copy genes for PCR amplification and SNP discovery. With these primers and subsequently designed genome-specific primers, a total of 1,527 loci were found to contain one or more genome-specific SNPs.

Conclusion

The ConservedPrimers 2.0 pipeline for designing intron-flanking primers was developed and its utility demonstrated. The tool can be used for SNP discovery, genetic variation assays and marker development for any target genome that has abundant ESTs and a related reference genome that has been fully sequenced. The ConservedPrimers 2.0 pipeline has been implemented as a command-line tool as well as a web application. Both versions are freely available at .

Rapid development of PCR-based genome-specific repetitive DNA junction markers in wheat

In hexaploid wheat (Triticum aestivum L.) (AABBDD, C=17 000 Mb), repeat DNA accounts for approximately 90% of the genome, of which transposable elements (TEs) constitute 60%-80%. Despite the dynamic evolution of TEs, our previous study indicated that the majority of TEs are conserved and collinear between the homologous wheat genomes, based on identical insertion patterns. In this study, we exploited the unique and abundant TE insertion junction regions identified from diploid Aegilops tauschii to develop genome-specific repeat DNA junction markers (RJM) for use in hexaploid wheat. In this study, both BAC end and random shotgun sequences were used to search for RJM. Of the 300 RJM primer pairs tested, 269 (90%) amplified single bands from diploid Ae. tauschii. Of these 269 primer pairs, 260 (97%) amplified hexaploid wheat and 9 (3%) amplified Ae. tauschii only. Among the RJM primers that amplified hexaploid wheat, 88% were successfully assigned to individual chromosomes of the hexaploid D genome. Among the 38 RJM primers mapped on chromosome 6D, 31 (82%) were unambiguously mapped to delineated bins of the chromosome using various wheat deletion lines. Our results suggest that the unique RJM derived from the diploid D genome could facilitate genetic, physical, and radiation mapping of the hexaploid wheat D genome.

Development of SSR markers and analysis of diversity in Turkish populations of Brachypodium distachyon

BACKGROUND

Brachypodium distachyon (Brachypodium) is rapidly emerging as a powerful model system to facilitate research aimed at improving grass crops for grain, forage and energy production. To characterize the natural diversity of Brachypodium and provide a valuable new tool to the growing list of resources available to Brachypodium researchers, we created and characterized a large, diverse collection of inbred lines.

RESULTS

We developed 84 inbred lines from eight locations in Turkey. To enable genotypic characterization of this collection, we created 398 SSR markers from BAC end and EST sequences. An analysis of 187 diploid lines from 56 locations with 43 SSR markers showed considerable genotypic diversity. There was some correlation between SSR genotypes and broad geographic regions, but there was also a high level of genotypic diversity at individual locations. Phenotypic analysis of this new germplasm resource revealed considerable variation in flowering time, seed size, and plant architecture. The inbreeding nature of Brachypodium was confirmed by an extremely high level of homozygosity in wild plants and a lack of cross-pollination under laboratory conditions.

CONCLUSION

Taken together, the inbreeding nature and genotypic diversity observed at individual locations suggest a significant amount of long-distance seed dispersal. The resources developed in this study are freely available to the research community and will facilitate experimental applications based on natural diversity.

Novel microRNAs uncovered by deep sequencing of small RNA transcriptomes in bread wheat (Triticum aestivum L.) and Brachypodium distachyon (L.) Beauv

The small RNA transcriptomes of bread wheat and its emerging model Brachypodium distachyon were obtained by using deep sequencing technology. Small RNA compositions were analyzed in these two species. In addition to 70 conserved microRNAs (miRNAs) from 25 families, 23 novel wheat miRNAs were identified. For Brachypodium, 12 putative miRNAs were predicted from a limited number of expressed sequence tags, of which one was a potential novel miRNA. Also, 94 conserved miRNAs from 28 families were identified in this species. Expression validation was performed for several novel wheat miRNAs. RNA ligase-mediated 5' rapid amplification of complementary DNA ends experiments demonstrated their capability to cleave predicted target genes including three disease-resistant gene analogs. Differential expression of miRNAs was observed between Brachypodium vegetative and reproductive tissues, suggesting their different roles at the two growth stages. Our work significantly increases the novel miRNA numbers in wheat and provides the first set of small RNAs in B. distachyon.

Rapid development of PCR-based genome-specific repetitive DNA junction markers in wheat

In hexaploid wheat (Triticum aestivum L.) (AABBDD, C=17 000 Mb), repeat DNA accounts for approximately 90% of the genome, of which transposable elements (TEs) constitute 60%-80%. Despite the dynamic evolution of TEs, our previous study indicated that the majority of TEs are conserved and collinear between the homologous wheat genomes, based on identical insertion patterns. In this study, we exploited the unique and abundant TE insertion junction regions identified from diploid Aegilops tauschii to develop genome-specific repeat DNA junction markers (RJM) for use in hexaploid wheat. In this study, both BAC end and random shotgun sequences were used to search for RJM. Of the 300 RJM primer pairs tested, 269 (90%) amplified single bands from diploid Ae. tauschii. Of these 269 primer pairs, 260 (97%) amplified hexaploid wheat and 9 (3%) amplified Ae. tauschii only. Among the RJM primers that amplified hexaploid wheat, 88% were successfully assigned to individual chromosomes of the hexaploid D genome. Among the 38 RJM primers mapped on chromosome 6D, 31 (82%) were unambiguously mapped to delineated bins of the chromosome using various wheat deletion lines. Our results suggest that the unique RJM derived from the diploid D genome could facilitate genetic, physical, and radiation mapping of the hexaploid wheat D genome.

Structural characterization of Brachypodium genome and its syntenic relationship with rice and wheat

Brachypodium distachyon (Brachypodium) has been recently recognized as an emerging model system for both comparative and functional genomics in grass species. In this study, 55,221 repeat masked Brachypodium BAC end sequences (BES) were used for comparative analysis against the 12 rice pseudomolecules. The analysis revealed that approximately 26.4% of BES have significant matches with the rice genome and 82.4% of the matches were homologous to known genes. Further analysis of paired-end BES and approximately 1.0 Mb sequences from nine selected BACs proved to be useful in revealing conserved regions and regions that have undergone considerable genomic changes. Differential gene amplification, insertions/deletions and inversions appeared to be the common evolutionary events that caused variations of microcolinearity at different orthologous genomic regions. It was found that approximately 17% of genes in the two genomes are not colinear in the orthologous regions. Analysis of BAC sequences also revealed higher gene density (approximately 9 kb/gene) and lower repeat DNA content (approximately 13.1%) in Brachypodium when compared to the orthologous rice regions, consistent with the smaller size of the Brachypodium genome. The 119 annotated Brachypodium genes were BLASTN compared against the wheat EST database and deletion bin mapped wheat ESTs. About 77% of the genes retrieved significant matches in the EST database, while 9.2% matched to the bin mapped ESTs. In some cases, genes in single Brachypodium BACs matched to multiple ESTs that were mapped to the same deletion bins, suggesting that the Brachypodium genome will be useful for ordering wheat ESTs within the deletion bins and developing specific markers at targeted regions in the wheat genome.

BatchPrimer3: a high throughput web application for PCR and sequencing primer design

Background

Microsatellite (simple sequence repeat – SSR) and single nucleotide polymorphism (SNP) markers are two types of important genetic markers useful in genetic mapping and genotyping. Often, large-scale genomic research projects require high-throughput computer-assisted primer design. Numerous such web-based or standard-alone programs for PCR primer design are available but vary in quality and functionality. In particular, most programs lack batch primer design capability. Such a high-throughput software tool for designing SSR flanking primers and SNP genotyping primers is increasingly demanded.

Results

A new web primer design program, BatchPrimer3, is developed based on Primer3. BatchPrimer3 adopted the Primer3 core program as a major primer design engine to choose the best primer pairs. A new score-based primer picking module is incorporated into BatchPrimer3 and used to pick position-restricted primers. BatchPrimer3 v1.0 implements several types of primer designs including generic primers, SSR primers together with SSR detection, and SNP genotyping primers (including single-base extension primers, allele-specific primers, and tetra-primers for tetra-primer ARMS PCR), as well as DNA sequencing primers. DNA sequences in FASTA format can be batch read into the program. The basic information of input sequences, as a reference of parameter setting of primer design, can be obtained by pre-analysis of sequences. The input sequences can be pre-processed and masked to exclude and/or include specific regions, or set targets for different primer design purposes as in Primer3Web and primer3Plus. A tab-delimited or Excel-formatted primer output also greatly facilitates the subsequent primer-ordering process. Thousands of primers, including wheat conserved intron-flanking primers, wheat genome-specific SNP genotyping primers, and Brachypodium SSR flanking primers in several genome projects have been designed using the program and validated in several laboratories.

Conclusion

BatchPrimer3 is a comprehensive web primer design program to develop different types of primers in a high-throughput manner. Additional methods of primer design can be easily integrated into future versions of BatchPrimer3. The program with source code and thousands of PCR and sequencing primers designed for wheat and Brachypodium are accessible at .

BatchPrimer3: a high throughput web application for PCR and sequencing primer design

BACKGROUND

Microsatellite (simple sequence repeat - SSR) and single nucleotide polymorphism (SNP) markers are two types of important genetic markers useful in genetic mapping and genotyping. Often, large-scale genomic research projects require high-throughput computer-assisted primer design. Numerous such web-based or standard-alone programs for PCR primer design are available but vary in quality and functionality. In particular, most programs lack batch primer design capability. Such a high-throughput software tool for designing SSR flanking primers and SNP genotyping primers is increasingly demanded.

RESULTS

A new web primer design program, BatchPrimer3, is developed based on Primer3. BatchPrimer3 adopted the Primer3 core program as a major primer design engine to choose the best primer pairs. A new score-based primer picking module is incorporated into BatchPrimer3 and used to pick position-restricted primers. BatchPrimer3 v1.0 implements several types of primer designs including generic primers, SSR primers together with SSR detection, and SNP genotyping primers (including single-base extension primers, allele-specific primers, and tetra-primers for tetra-primer ARMS PCR), as well as DNA sequencing primers. DNA sequences in FASTA format can be batch read into the program. The basic information of input sequences, as a reference of parameter setting of primer design, can be obtained by pre-analysis of sequences. The input sequences can be pre-processed and masked to exclude and/or include specific regions, or set targets for different primer design purposes as in Primer3Web and primer3Plus. A tab-delimited or Excel-formatted primer output also greatly facilitates the subsequent primer-ordering process. Thousands of primers, including wheat conserved intron-flanking primers, wheat genome-specific SNP genotyping primers, and Brachypodium SSR flanking primers in several genome projects have been designed using the program and validated in several laboratories.

CONCLUSION

BatchPrimer3 is a comprehensive web primer design program to develop different types of primers in a high-throughput manner. Additional methods of primer design can be easily integrated into future versions of BatchPrimer3. The program with source code and thousands of PCR and sequencing primers designed for wheat and Brachypodium are accessible at http://wheat.pw.usda.gov/demos/BatchPrimer3/.

The nuclear genome of Brachypodium distachyon: analysis of BAC end sequences

Due in part to its small genome (approximately 350 Mb), Brachypodium distachyon is emerging as a model system for temperate grasses, including important crops like wheat and barley. We present the analysis of 10.9% of the Brachypodium genome based on 64,696 bacterial artificial chromosome (BAC) end sequences (BES). Analysis of repeat DNA content in BES revealed that approximately 11.0% of the genome consists of known repetitive DNA. The vast majority of the Brachypodium repetitive elements are LTR retrotransposons. While Bare-1 retrotransposons are common to wheat and barley, Brachypodium repetitive element sequence-1 (BRES-1), closely related to Bare-1, is also abundant in Brachypodium. Moreover, unique Brachypodium repetitive element sequences identified constitute approximately 7.4% of its genome. Simple sequence repeats from BES were analyzed, and flanking primer sequences for SSR detection potentially useful for genetic mapping are available at http://brachypodium.pw.usda.gov . Sequence analyses of BES indicated that approximately 21.2% of the Brachypodium genome represents coding sequence. Furthermore, Brachypodium BES have more significant matches to ESTs from wheat than rice or maize, although these species have similar sizes of EST collections. A phylogenetic analysis based on 335 sequences shared among seven grass species further revealed a closer relationship between Brachypodium and Triticeae than Brachypodium and rice or maize.

Rapid evolution and complex structural organization in genomic regions harboring multiple prolamin genes in the polyploid wheat genome

Genes encoding wheat prolamins belong to complicated multi-gene families in the wheat genome. To understand the structural complexity of storage protein loci, we sequenced and analyzed orthologous regions containing both gliadin and LMW-glutenin genes from the A and B genomes of a tetraploid wheat species, Triticum turgidum ssp. durum. Despite their physical proximity to one another, the gliadin genes and LMW-glutenin genes are organized quite differently. The gliadin genes are found to be more clustered than the LMW-glutenin genes which are separated from each other by much larger distances. The separation of the LMW-glutenin genes is the result of both the insertion of large blocks of repetitive DNA owing to the rapid amplification of retrotransposons and the presence of genetic loci interspersed between them. Sequence comparisons of the orthologous regions reveal that gene movement could be one of the major factors contributing to the violation of microcolinearity between the homoeologous A and B genomes in wheat. The rapid sequence rearrangements and differential insertion of repetitive DNA has caused the gene islands to be not conserved in compared regions. In addition, we demonstrated that the i-type LMW-glutenin originated from a deletion of 33-bps in the 5' coding region of the m-type gene. Our results show that multiple rounds of segmental duplication of prolamin genes have driven the amplification of the omega-gliadin genes in the region; such segmental duplication could greatly increase the repetitive DNA content in the genome depending on the amount of repetitive DNA present in the original duplicate region.

Construction and characterization of two BAC libraries from Brachypodium distachyon, a new model for grass genomics

Brachypodium is well suited as a model system for temperate grasses because of its compact genome and a range of biological features. In an effort to develop resources for genome research in this emerging model species, we constructed 2 bacterial artificial chromosome (BAC) libraries from an inbred diploid Brachypodium distachyon line, Bd21, using restriction enzymes HindIII and BamHI. A total of 73,728 clones (36,864 per BAC library) were picked and arrayed in 192,384-well plates. The average insert size for the BamHI and HindIII libraries is estimated to be 100 and 105 kb, respectively, and inserts of chloroplast origin account for 4.4% and 2.4%, respectively. The libraries individually represent 9.4- and 9.9-fold haploid genome equivalents with combined 19.3-fold genome coverage, based on a genome size of 355 Mb reported for the diploid Brachypodium, implying a 99.99% probability that any given specific sequence will be present in each library. Hybridization of the libraries with 8 starch biosynthesis genes was used to empirically evaluate this theoretical genome coverage; the frequency at which these genes were present in the library clones gave an estimated coverage of 11.6- and 19.6-fold genome equivalents. To obtain a first view of the sequence composition of the Brachypodium genome, 2185 BAC end sequences (BES) representing 1.3 Mb of random genomic sequence were compared with the NCBI GenBank database and the GIRI repeat database. Using a cutoff expectation value of E<10-10, only 3.3% of the BESs showed similarity to repetitive sequences in the existing database, whereas 40.0% had matches to the sequences in the EST database, suggesting that a considerable portion of the Brachypodium genome is likely transcribed. When the BESs were compared with individual EST databases, more matches hit wheat than maize, although their EST collections are of a similar size, further supporting the close relationship between Brachypodium and the Triticeae. Moreover, 122 BESs have significant matches to wheat ESTs mapped to individual chromosome bin positions. These BACs represent colinear regions containing the mapped wheat ESTs and would be useful in identifying additional markers for specific wheat chromosome regions.