Sequencing and assembly of low copy and genic regions of isolated Triticum aestivum chromosome arm 7DS

The genome of bread wheat (Triticum aestivum) is predicted to be greater than 16 Gbp in size and consist predominantly of repetitive elements, making the sequencing and assembly of this genome a major challenge. We have reduced genome sequence complexity by isolating chromosome arm 7DS and applied second‐generation technology and appropriate algorithmic analysis to sequence and assemble low copy and genic regions of this chromosome arm. The assembly represents approximately 40% of the chromosome arm and all known 7DS genes. Comparison of the 7DS assembly with the sequenced genomes of rice (Oryza sativa) and Brachypodium distachyon identified large regions of conservation. The syntenic relationship between wheat, B. distachyon and O. sativa, along with available genetic mapping data, has been used to produce an annotated draft 7DS syntenic build, which is publicly available at http://www.wheatgenome.info . Our results suggest that the sequencing of isolated chromosome arms can provide valuable information of the gene content of wheat and is a step towards whole‐genome sequencing and variation discovery in this important crop.     

Lr34 multi-pathogen resistance ABC transporter: molecular analysis of homoeologous and orthologous genes in hexaploid wheat and other grass species

The Triticum aestivum (bread wheat) disease resistance gene Lr34 confers durable, race non-specific protection against three fungal pathogens, and has been a highly relevant gene for wheat breeding since the green revolution. Lr34, located on chromosome 7D, encodes an ATP-binding cassette (ABC) transporter. Both wheat cultivars with and without Lr34-based resistance encode a putatively functional protein that differ by only two amino acid polymorphisms. In this study, we focused on the identification and characterization of homoeologous and orthologous Lr34 genes in hexaploid wheat and other grasses. In hexaploid wheat we found an expressed and putatively functional Lr34 homoeolog located on chromosome 4A, designated Lr34-B. Another homoeologous Lr34 copy, located on chromosome 7A, was disrupted by the insertion of repetitive elements. Protein sequences of LR34-B and LR34 were 97% identical. Orthologous Lr34 genes were detected in the genomes of Oryza sativa (rice) and Sorghum bicolor (sorghum). Zea mays (maize), Brachypodium distachyon and Hordeum vulgare (barley) lacked Lr34 orthologs, indicating independent deletion of this particular ABC transporter. Lr34 was part of a gene-rich island on the wheat D genome. We found gene colinearity on the homoeologous A and B genomes of hexaploid wheat, but little microcolinearity in other grasses. The homoeologous LR34-B protein and the orthologs from rice and sorghum have the susceptible haplotype for the two critical polymorphisms distinguishing the LR34 proteins from susceptible and resistant wheat cultivars. We conclude that the particular Lr34-haplotype found in resistant wheat cultivars is unique. It probably resulted from functional gene diversification that occurred after the polyploidization event that was at the origin of cultivated bread wheat.     

Sequencing wheat chromosome arm 7BS delimits the 7BS/4AL translocation and reveals homoeologous gene conservation

Complex Triticeae genomes pose a challenge to genome sequencing efforts due to their size and repetitive nature. Genome sequencing can reveal details of conservation and rearrangements between related genomes. We have applied Illumina second generation sequencing technology to sequence and assemble the low copy and unique regions of Triticum aestivum chromosome arm 7BS, followed by the construction of a syntenic build based on gene order in Brachypodium. We have delimited the position of a previously reported translocation between 7BS and 4AL with a resolution of one or a few genes and report approximately 13% genes from 7BS having been translocated to 4AL. An additional 13 genes are found on 7BS which appear to have originated from 4AL. The gene content of the 7DS and 7BS syntenic builds indicate a total of ~77,000 genes in wheat. Within wheat syntenic regions, 7BS and 7DS share 740 genes and a common gene conservation rate of ~39% of the genes from the corresponding regions in Brachypodium, as well as a common rate of colinearity with Brachypodium of ~60%. Comparison of wheat homoeologues revealed ~84% of genes previously identified in 7DS have a homoeologue on 7BS or 4AL. The conservation rates we have identified among wheat homoeologues and with Brachypodium provide a benchmark of homoeologous gene conservation levels for future comparative genomic analysis. The syntenic build of 7BS is publicly available at http://www.wheatgenome.info.     

Identification of transposons,retroelements, and a gene family predominantly expressed in floral tissues in chromosome 3DS of the hexaploid wheat progenitor <Emphasis Type="Italic">Aegilops tauschii</Emphasis>

A multigene family expressed during early floral development was identified on the short arm of wheat chromosome 3D in the region of the Ph2 locus, a locus controlling homoeologous chromosome pairing in allohexaploid wheat. Physical, genetic and molecular characterisation of the Wheat Meiosis 1 (WM1) gene family identified seven members that localised within a region of 173-kb. WM1 gene family members were sequenced and they encode mainly type Ia plasma membrane-anchored leucine rich repeat-like receptor proteins. In situ expression profiling suggests the gene family is predominantly expressed in floral tissue. In addition to the WM1 gene family, a number of other genes, gene fragments and pseudogenes were identified. It has been predicted that there is approximately one gene every 19-kb and that this region of the wheat genome contains 23 repetitive elements including BARE-1 and Wis2-1 like sequences. Nearly 50% of the repetitive elements identified were similar to known transposons from the CACTA superfamily. Ty1-copia, Ty3-gypsy and Athila LTR retroelements were also prevalent within the region. The WM1 gene cluster is present on 3DS and on barley 3HS but missing from the A and B genomes of hexaploid wheat. This suggests either recent generation of the cluster or specific deletion of the cluster during wheat polyploidisation. The evolutionary significance of the cluster, its possible roles in disease response or floral and early meiotic development and its location at or near the Ph2 locus are discussed.     

Characterizing the composition and evolution of homoeologous genomes in hexaploid wheat through BAC-end sequencing on chromosome 3B

Bread wheat (Triticum aestivum) is one of the most important crops worldwide. However, because of its large, hexaploid, highly repetitive genome it is a challenge to develop efficient means for molecular analysis and genetic improvement in wheat. To better understand the composition and molecular evolution of the hexaploid wheat homoeologous genomes and to evaluate the potential of BAC-end sequences (BES) for marker development, we have followed a chromosome-specific strategy and generated 11 Mb of random BES from chromosome 3B, the largest chromosome of bread wheat. The sequence consisted of about 86% of repetitive elements, 1.2% of coding regions, and 13% remained unknown. With 1.2% of the sequence length corresponding to coding sequences, 6000 genes were estimated for chromosome 3B. New repetitive sequences were identified, including a Triticineae-specific tandem repeat (Fat) that represents 0.6% of the B-genome and has been differentially amplified in the homoeologous genomes before polyploidization. About 10% of the BES contained junctions between nested transposable elements that were used to develop chromosome-specific markers for physical and genetic mapping. Finally, sequence comparison with 2.9 Mb of random sequences from the D-genome of Aegilops tauschii suggested that the larger size of the B-genome is due to a higher content in repetitive elements. It also indicated which families of transposable elements are mostly responsible for differential expansion of the homoeologous wheat genomes during evolution. Our data demonstrate that BAC-end sequencing from flow-sorted chromosomes is a powerful tool for analysing the structure and evolution of polyploid and highly repetitive genomes.     

Genome analysis at different ploidy levels allows cloning of the powdery mildew resistance gene Pm3b from hexaploid wheat

In wheat, race-specific resistance to the fungal pathogen powdery mildew (Blumeria graminis f. sp. tritici) is controlled by the Pm genes. There are 10 alleles conferring resistance at the Pm3 locus (Pm3a to Pm3j) on chromosome 1AS of hexaploid bread wheat (Triticum aestivum L.). The genome of hexaploid wheat has a size of 1.6 x 1010 bp and contains more than 80% of repetitive sequences, making positional cloning difficult. Here, we demonstrate that the combined analysis of genomes from wheat species with different ploidy levels can be exploited for positional cloning in bread wheat. We have mapped the Pm3b gene in hexaploid wheat to a genetic interval of 0.97 centimorgan (cM). The diploid T. monococcum and the tetraploid T. turgidum ssp. durum provided models for the A genome of hexaploid wheat and allowed to establish a physical contig spanning the Pm3 locus. Although the haplotypes at the Pm3 locus differed markedly between the three species, a large resistance gene-like family specific to wheat group 1 chromosomes was consistently found at the Pm3 locus. A candidate gene for Pm3b was identified using partial sequence conservation between resistant line Chul and T. monococcum cv. DV92. A susceptible Pm3b mutant, carrying a single-base pair deletion in the coding region of the candidate gene was isolated. When tested in a single cell transformation assay, the Pm3b candidate gene conferred race-specific resistance to powdery mildew. These results demonstrate that the candidate gene, a member of the coiled-coil nucleotide binding site leucine-rich repeat (NBS-LRR) type of disease resistance genes, is the Pm3b gene.     

单个拟穴青蟹核内多拷贝线粒体COI假基因序列的研究(简报)

线粒体DNA(mtDNA)因它的高拷贝数、易扩增、高突变率、中性、低重组和母系遗传等特征,已广泛应用在系统进化和群体遗传研究方面[1,2].但是,由于mtDNA本身的特征,如异质性.     

Frequent gene movement and pseudogene evolution is common to the large and complex genomes of wheat, barley, and their relatives

All six arms of the group 1 chromosomes of hexaploid wheat (Triticum aestivum) were sequenced with Roche/454 to 1.3- to 2.2-fold coverage and compared with similar data sets from the homoeologous chromosome 1H of barley (Hordeum vulgare). Six to ten thousand gene sequences were sampled per chromosome. These were classified into genes that have their closest homologs in the Triticeae group 1 syntenic region in Brachypodium, rice (Oryza sativa), and/or sorghum (Sorghum bicolor) and genes that have their homologs elsewhere in these model grass genomes. Although the number of syntenic genes was similar between the homologous groups, the amount of nonsyntenic genes was found to be extremely diverse between wheat and barley and even between wheat subgenomes. Besides a small core group of genes that are nonsyntenic in other grasses but conserved among Triticeae, we found thousands of genic sequences that are specific to chromosomes of one single species or subgenome. By examining in detail 50 genes from chromosome 1H for which BAC sequences were available, we found that many represent pseudogenes that resulted from transposable element activity and double-strand break repair. Thus, Triticeae seem to accumulate nonsyntenic genes frequently. Since many of them are likely to be pseudogenes, total gene numbers in Triticeae are prone to pronounced overestimates.     

Genetic mapping of the wheat dimeric α-amylase inhibitor multi-gene family using allele-specific primers based on intergenomic SNPs

Single nucleotide polymorphisms (SNPs) identified in EST sequences can be used to map expressed genes. Though SNPs are useful markers for genetic mapping, SNP mapping of genes in common wheat is challenging because the genetic complement of wheat consists of three closely related genomes (designated A, B, and D), and most genes are present in triplicate sets. Mapping multi-gene family members is further complicated by the fact that it is difficult to distinguish SNP differences between the various paralogs from those between the different genomes. We have developed a PCR-based method for assigning wheat EST sequences to their proper genetic loci by first identifying and mapping SNPs that distinguish the three genomes. To develop this method, we focused on EST sequences encoding the dimeric α-amylase inhibitors (WDAI), The WDAI coding regions of hexaploid wheat were aligned. The sequences were classified into three groups based on nucleotide variations. Twenty-two SNPs were identified that distinguish the three groups. Group-specific primers based on these SNPs were designed to permit selective amplification of each group. The chromosomal location of each group was then determined using Group 3 ditelosomic lines of Chinese Spring. Groups 1 and 2 were assigned to chromosome locations 3DS and 3BS, respectively, whereas no sequence could be assigned to 3AS. A remarkable feature of this method is the ability to discriminate the location of homoeologous multigenes in the three genomes of wheat. This strategy can be useful for assigning unknown wheat EST sequences to specific chromosomes.     

Genetic and physical mapping of the patatin genes in potato and tomato

Summary Genes for the major storage protein of potato, patatin, have been mapped genetically and physically in both the potato and tomato genomes. In potato, all patatin genes detected by the cDNA clone pGM01 map to a single locus at the end of the long arm of chromosome 8. By means of pulsed field gel electrophoresis (PFGE) it was possible further to delimit this locus, containing 10–15 copies of the gene, to a maximum size of 1.4 million base pairs. Hybridizations with class-specific clones suggest that the locus is at least partially divided into domains containing the two major types of patatin genes, class I and II. In tomato, patatin-homologous sequences were found to reside at the orthologous locus at the end of chromosome 8. The approximately three copies in tomato were localized by PFGE to a single fragment of 300 kilobases. Whereas the class II-specific 5 promoter sequences reside in tomato at the same locus as the coding sequences, the single class I-specific copy of the 5 promoter sequences was localized on chromosome 3 with no coding sequence attached to it. A clone from this chromosome 3 locus of tomato was isolated and by restriction fragment length polymorphism mapping it could be further shown that a similar class I-specific sequence also exists on chromosome 3 of potato. As in tomato, this copy on chromosome 3 is not linked to a coding sequence for patatin. The results are discussed with respect to genome evolution and PFGE analysis of complex gene families.     

Cell organisation,sulphur metabolism and ion transport-related genes are differentially expressed in Paracoccidioides brasiliensis mycelium and yeast cells

Background

Bread wheat (Triticum aestivum) is an important staple food. However, wheat gluten proteins cause celiac disease (CD) in 0.5 to 1% of the general population. Among these proteins, the α-gliadins contain several peptides that are associated to the disease.

Results

We obtained 230 distinct α-gliadin gene sequences from severaldiploid wheat species representing the ancestral A, B, and D genomes of the hexaploid bread wheat. The large majority of these sequences (87%) contained an internal stop codon. All α-gliadin sequences could be distinguished according to the genome of origin on the basis of sequence similarity, of the average length of the polyglutamine repeats, and of the differences in the presence of four peptides that have been identified as T cell stimulatory epitopes in CD patients through binding to HLA-DQ2/8. By sequence similarity, α-gliadins from the public database of hexaploid T. aestivum could be assigned directly to chromosome 6A, 6B, or 6D. T. monococcum (A genome) sequences, as well as those from chromosome 6A of bread wheat, almost invariably contained epitope glia-α9 and glia-α20, but never the intact epitopes glia-α and glia-α2. A number of sequences from T. speltoides, as well as a number of sequences fromchromosome 6B of bread wheat, did not contain any of the four T cell epitopes screened for. The sequences from T. tauschii (D genome), as well as those from chromosome 6D of bread wheat, were found to contain all of these T cell epitopes in variable combinations per gene. The differences in epitope composition resulted mainly from point mutations. These substitutions appeared to be genome specific.

Conclusion

Our analysis shows that α-gliadin sequences from the three genomes of bread wheat form distinct groups. The four known T cell stimulatory epitopes are distributed non-randomly across the sequences, indicating that the three genomes contribute differently to epitope content. A systematic analysis of all known epitopes in gliadins and glutenins will lead to better understanding of the differences in toxiCity among wheat varieties. On the basis of such insight, breeding strategies can be designed to generate less toxic varieties of wheat which may be tolerated by at least part of the CD patient population.     

Exploring the tertiary gene pool of bread wheat: sequence assembly and analysis of chromosome 5Mg of Aegilops geniculata

Next‐generation sequencing (NGS) provides a powerful tool for the discovery of important genes and alleles in crop plants and their wild relatives. Despite great advances in NGS technologies, whole‐genome shotgun sequencing is cost‐prohibitive for species with complex genomes. An attractive option is to reduce genome complexity to a single chromosome prior to sequencing. This work describes a strategy for studying the genomes of distant wild relatives of wheat by isolating single chromosomes from addition or substitution lines, followed by chromosome sorting using flow cytometry and sequencing of chromosomal DNA by NGS technology. We flow‐sorted chromosome 5M^g from a wheat/Aegilops geniculata disomic substitution line [DS5M^g (5D)] and sequenced it using an Illumina HiSeq 2000 system at approximately 50 × coverage. Paired‐end sequences were assembled and used for structural and functional annotation. A total of 4236 genes were annotated on 5M^g, in close agreement with the predicted number of genes on wheat chromosome 5D (4286). Single‐gene FISH indicated no major chromosomal rearrangements between chromosomes 5M^g and 5D. Comparing chromosome 5M^g with model grass genomes identified synteny blocks in Brachypodium distachyon, rice (Oryza sativa), sorghum (Sorghum bicolor) and barley (Hordeum vulgare). Chromosome 5M^g‐specific SNPs and cytogenetic probe‐based resources were developed and validated. Deletion bin‐mapped and ordered 5M^g SNP markers will be useful to track 5M‐specific introgressions and translocations. This study provides a detailed sequence‐based analysis of the composition of a chromosome from a distant wild relative of bread wheat, and opens up opportunities to develop genomic resources for wild germplasm to facilitate crop improvement.     

Chromosome mapping and phylogenetic analysis of the cytosolic acetyl-CoA carboxylase loci in wheat

The cytosolic isoform of plant acetyl-CoA carboxylase is a multidomain enzyme involved in the synthesis of very-long-chain fatty acids and in secondary metabolism. Chromosome mapping of wheat identified one locus containing cytosolic acetyl-CoA carboxylase genes (Acc-2) and a related partially processed pseudogene (Psi-Acc-2) in the distal region of the long arm of wheat homoeologous group 3 chromosomes. Multiple copies of the Acc-2 genes, whose presence was suggested by sequence analysis, are likely to be arranged in tandem repeats. At least three out of five genes cloned from hexaploid wheat map to this locus. Another locus containing Acc-2--related sequences is present in the distal region of the long arm of chromosome 5D. The identity of the hybridizing DNA present at this locus remains unknown. A system based on PCR-cloning and DNA sequence analysis of acetyl-CoA carboxylase genes was developed to address various phylogenetic and systematics questions in grasses. It was applied to reconstruct the phylogeny of the Acc-2 genes from D- and S-genome Aegilops and A-genome Triticum diploid species, AABB- and AAGG-genome tetraploid wheat, and AABBDD-genome hexaploid wheat, as well as from rye and barley. The combined cytogenetic and molecular evolution approach allowed assignment of gene sequences included in phylogenetic analysis to specific loci on homoeologous chromosomes. Recurring gene duplication followed by chromosome translocation and/or possible loss of some gene copies, as well as loss of introns, occurred in the gene family in different plant lineages. Two major Acc-2 clades appeared before the divergence of barley and rye. Nucleotide substitution rates in different parts of the Acc-2 gene were assessed. This analysis of the Acc-2 loci provides detailed information regarding evolutionary events at a low--copy-number locus containing important functional genes. These events are likely to be common and to play a significant role in shaping grass genomes.     

The characterization and comparative analysis of high-molecular-weight glutenin genes from genomes A and B of a hexaploid bread wheat

Summary Two high-molecular-weight subunit (HMWS) glutenin genes from the A and B genomes of the hexaploid bread wheat Triticum aestivum L. cv Cheyenne have been isolated and sequenced. Both of these genes are of the high M_r class (x-type) of HMW glutenins, and have not been previously reported. The entire set of six HMW genes from cultivar Cheyenne have now been isolated and characterized. An analysis of the Ax and Bx sequences shows that the Ax sequence is similar to the homoeologous gene from the D genome, while the Bx repeat structure is significantly different. The repetitive region of these proteins can be modelled as a series of interspersed copies of repeat modifs of 6, 9, and 15 amino acid residues. The evolution of these genes includes single-base substitutions over the entire coding region, plus insertion/deletions of single or blocks of repeats in the central repetitive domain.     

Homoeologous set of NBS-LRR genes located at leaf and stripe rust resistance loci on short arms of chromosome 1 of wheat

Homoeologous group 1 chromosomes of wheat contain important genes that confer resistance to leaf, stem and stripe rusts, powdery mildew and Russian wheat aphid. A disease resistance gene analog encoding nucleotide binding site-leucine rich repeat (NBS-LRR), designated RgaYr10, was previously identified at the stripe rust resistant locus, Yr10, located on chromosome 1BS distal to the storage protein, Gli-B1 locus. RgaYr10 identified gene members in the homoeologous region of chromosome 1DS cosegregating with the leaf rust resistance gene, Lr21, which originally was transferred from a diploid D genome progenitor. Four RgaYr10 gene members were isolated from chromosome 1DS and compared to two gene members previously isolated from the chromosome 1BS homeologue. NBS-LRR genes tightly linked to stripe rust resistance gene Yr10 on chromosome 1BS were closely related in sequence and structure to NBS-LRR genes tightly linked to leaf rust resistance gene Lr21 located within the homoeologous region on chromosome 1DS. The level of sequence homology was similar between NBS-LRR genes that were isolated from different genomes as compared to genes from the same genome. Electronic Publication     

Dynamics of the evolution of orthologous and paralogous portions of a complex locus region in two genomes of allopolyploid wheat

Two overlapping bacterial artificial chromosome (BAC) clones from the B genome of the tetraploid wheat Triticum turgidum were identified, each of which contains one of the two high-molecular-weight (HMW) glutenin genes, comprising the complex Glu-B1 locus. The complete sequence (285 506 bp of DNA) of this chromosomal region was determined. The two paralogous x-type ( Glu-1-1 ) and y-type ( Glu-1-2 ) HMW-glutenin genes of the complex Glu-B1 locus were found to be separated by ca. 168 000 bp instead of the 51 000 bp separation previously reported for the orthologous Glu-D1 locus of Aegilops tauschii, the D-genome donor of hexaploid wheat. This difference in intergene spacing is due almost entirely to be the insertion of clusters of nested retrotransposons. Otherwise, the orientation and order of the HMW glutenins and adjacent genes were identical in the two genomes. A comparison of these orthologous regions indicates modes and patterns of sequence divergence, with implications for the overall Triticeae genome structure and evolution. A duplicate globulin gene, found 5' of each HMW-glutenin gene, assists to tentatively define the original duplication event leading to the paralogous x- and y-type HMW-glutenin genes. The intergenic regions of the two loci are composed of different patterns and classes of retrotransposons, indicating that insertion times of these retroelements were after the divergence of the two wheat genomes. In addition, a putative receptor kinase gene near the y-type HMW-glutenin gene at the Glu-B1 locus is likely active as it matches recently reported ESTs from germinating barley endosperm. The presence of four genes represented only in the Triticeae endosperm ESTs suggests an endosperm-specific chromosome domain.