Bread Wheat Genome Research Articles

Genome-wide association study for morphological, phenological, quality, and yield traits in einkorn (Triticum monococcum L. subsp. monococcum).

Einkorn (Triticum monococcum L. subsp. monococcum, 2n = 2× = 14, AmAm) is a diploid wheat whose cultivation was widespread in the Mediterranean and European area till the Bronze Age, before it was replaced by the more productive durum and bread wheats. Although scarcely cultivated nowadays, it has gained renewed interest due to its relevant nutritional properties and as source of genetic diversity for crop breeding. However, the molecular basis of many traits of interest in einkorn remain still unknown. A panel of 160 einkorn landraces, from different parts of the distribution area, was characterized for several phenotypic traits related to morphology, phenology, quality, and yield for 4 years in two locations. An approach based on co-linearity with the A genome of bread wheat, supported also by that with Triticum urartu genome, was exploited to perform association mapping, even without an einkorn anchored genome. The association mapping approach uncovered numerous marker-trait associations; for 37 of these, a physical position was inferred by homology with the bread wheat genome. Moreover, numerous associated regions were also assigned to the available T. monococcum contigs. Among the intervals detected in this work, three overlapped with regions previously described as involved in the same trait, while four other regions were localized in proximity of loci previously described and presumably refer to the same gene/QTL. The remaining associated regions identified in this work could represent a novel and useful starting point for breeding approaches to improve the investigated traits in this neglected species.

Crossover-active regions of the wheat genome are distinguished by DMC1, the chromosome axis, H3K27me3, and signatures of adaptation.

The hexaploid bread wheat genome comprises over 16 gigabases of sequence across 21 chromosomes. Meiotic crossovers are highly polarized along the chromosomes, with elevation in the gene-dense distal regions and suppression in the Gypsy retrotransposon-dense centromere-proximal regions. We profiled the genomic landscapes of the meiotic recombinase DMC1 and the chromosome axis protein ASY1 in wheat and investigated their relationships with crossovers, chromatin state, and genetic diversity. DMC1 and ASY1 chromatin immunoprecipitation followed by sequencing (ChIP-seq) revealed strong co-enrichment in the distal, crossover-active regions of the wheat chromosomes. Distal ChIP-seq enrichment is consistent with spatiotemporally biased cytological immunolocalization of DMC1 and ASY1 close to the telomeres during meiotic prophase I. DMC1 and ASY1 ChIP-seq peaks show significant overlap with genes and transposable elements in the Mariner and Mutator superfamilies. However, DMC1 and ASY1 ChIP-seq peaks were detected along the length of each chromosome, including in low-crossover regions. At the fine scale, crossover elevation at DMC1 and ASY1 peaks and genes correlates with enrichment of the Polycomb histone modification H3K27me3. This indicates a role for facultative heterochromatin, coincident with high DMC1 and ASY1, in promoting crossovers in wheat and is reflected in distalized H3K27me3 enrichment observed via ChIP-seq and immunocytology. Genes with elevated crossover rates and high DMC1 and ASY1 ChIP-seq signals are overrepresented for defense-response and immunity annotations, have higher sequence polymorphism, and exhibit signatures of selection. Our findings are consistent with meiotic recombination promoting genetic diversity, shaping host–pathogen co-evolution, and accelerating adaptation by increasing the efficiency of selection.

Integrating the physical and genetic map of bread wheat facilitates the detection of chromosomal rearrangements

The bread wheat genome harbors a high content of repetitive DNA, which is amenable to detection and characterization using fluorescence in situ hybridization (FISH) karyotyping. An integrated genetic map was derived from a recombinant inbred population bred from a cross between a synthetic hexaploid wheat and a commercial Chinese bread wheat cultivar, based on 28 variable FISH sites and >150 000 single nucleotide polymorphism (SNP) loci. The majority (20/28) of the variable FISH sites were physically located within a chromosomal region consistent with the genetic location inferred from that of their co-segregating SNP loci. The eight exceptions reflected the presence of either a translocation (1R/1B, 1A/7A) or a presumptive intra-chromosomal inversion (4A). For eight out of the nine FISH sites detected on the Chinese Spring (CS) karyotype, there was a good match with the reference genome sequence, indicating that the most recent assembly has dealt well with the problem of placing tandem repeats. The integrated genetic map produced for wheat is informative as to the location of blocks of tandemly repeated DNA and can aid in improving the quality of the genome sequence assembly in regions surrounding these blocks.

Genome‐wide investigation and transcriptional profiling of the oxidosqualene cyclase (OSC) genes in wheat (Triticum aestivum)

AbstractTriterpene scaffolds are cyclized by differential functional genes of the oxidosqualene cyclase (OSC) gene family, which are widely present in plant species. Thirty‐nine OSCs were identified from the hexaploid bread wheat (Triticum aestivum L.) genome. The gene structure organization and the distribution of conserved motifs of OSCs gene family revealed that the TaOSC proteins are highly conserved. By phylogenetic analysis, all OSC proteins from wheat and other cereal crops were clustered into 11 subfamilies (I–XI). Based on the publicly available transcriptional profiling, the tissue‐specific expression patterns of TaOSC during several developmental stages revealed diverged expression patterns, indicating that the TaOSC family has undergone subfunctionalization. Furthermore, homeolog expression bias was commonly detected in the TaOSC subfamilies, such as the expression of Ta5A004900 and Ta5D011800 in the subfamily II were specifically expressed in the leaf sheath and flag leaf sheath. In addition, two genes were strongly induced by Pst, indicating that the two homolog genes are involved in the resistance to powdery mildew. Molecular evolution analysis indicated that the TaOSC family genes were undergone purifying selection. So, the powdery mildew response may be very conservative and important biological function of the Ta5A004900 and Ta5D011800 in wheat. In conclusion, these results can be used to explore the diversity of triterpenes and clarify the mechanism of triterpene cyclization, and also provide useful information for future studies on the function of the TaOSC gene family in bread wheat and other important crops.

Exploration of glutathione reductase for abiotic stress response in bread wheat (Triticum aestivum L.)

A total of seven glutathione reductase (GR) genes were identified in Triticum aestivum, which were used for comparative structural characterization, phylogenetic analysis and expression profiling with the GR genes of other cereal plants. The modulated gene expression and enzyme activity revealed the role of GRs in abiotic stress response in T. aestivum. Glutathione reductase (GR) is an enzymatic antioxidant that converts oxidized glutathione (GSSG) into reduced glutathione (GSH) through the ascorbate-glutathione cycle. In this study, a total of seven GR genes forming two homeologous groups were identified in the allohexaploid genome of bread wheat (Triticum aestivum). Besides, we identified three GR genes in each Aegilops tauschii, Brachypodium distachyon, Triticum urartu and Sorghum bicolor, which were used for comparative characterization. Phylogenetic analysis revealed the clustering of GR proteins into two groups; class I and class II, which were predicted to be localized in cytoplasm and chloroplast, respectively. The exon-intron and conserved motif patterns were almost conserved in each group, in which a maximum of 10 and 17 exons were present in chloroplastic and cytoplasmic GRs, respectively. The protein structure analysis confirmed the occurrence of conserved pyridine nucleotide disulfide oxidoreductase (Pyr_redox) and pyridine nucleotide disulfide oxidoreductase dimerization (Pyr_redox_dim) domains in each GR. The active site of GR proteins consisted of two conserved cysteine residues separated by four amino acid residues. Promoter analysis revealed the occurrence of growth and stress-related cis-active elements. Tissue-specific expression profiling suggested the involvement of GRs in both vegetative and reproductive tissue development in various plants. The differential expression of TaGR genes and enhanced GR enzyme activity suggested their roles under drought, heat, salt and arsenic stress. Interaction of GRs with other proteins and chemical compounds of the ascorbate-glutathione cycle revealed their coordinated functioning. The current study will provide a foundation for the validation of the precise role of each GR gene in future studies.

Annotation and Molecular Characterisation of the TaIRO3 and TaHRZ Iron Homeostasis Genes in Bread Wheat (Triticum aestivum L.).

Effective maintenance of plant iron (Fe) homoeostasis relies on a network of transcription factors (TFs) that respond to environmental conditions and regulate Fe uptake, translocation, and storage. The iron-related transcription factor 3 (IRO3), as well as haemerythrin motif-containing really interesting new gene (RING) protein and zinc finger protein (HRZ), are major regulators of Fe homeostasis in diploid species like Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa L.), but remain uncharacterised in hexaploid bread wheat (Triticum aestivum L.). In this study, we have identified, annotated, and characterised three TaIRO3 homoeologs and six TaHRZ1 and TaHRZ2 homoeologs in the bread wheat genome. Protein analysis revealed that TaIRO3 and TaHRZ proteins contain functionally conserved domains for DNA-binding, dimerisation, Fe binding, or polyubiquitination, and phylogenetic analysis revealed clustering of TaIRO3 and TaHRZ proteins with other monocot IRO3 and HRZ proteins, respectively. Quantitative reverse-transcription PCR analysis revealed that all TaIRO3 and TaHRZ homoeologs have unique tissue expression profiles and are upregulated in shoot tissues in response to Fe deficiency. After 24 h of Fe deficiency, the expression of TaHRZ homoeologs was upregulated, while the expression of TaIRO3 homoeologs was unchanged, suggesting that TaHRZ functions upstream of TaIRO3 in the wheat Fe homeostasis TF network.

Exploring the alpha-gliadin locus: the 33-mer peptide with six overlapping coeliac disease epitopes in Triticum aestivum is derived from a subgroup of Aegilops tauschii.

SummaryMost alpha‐gliadin genes of the Gli‐D2 locus on the D genome of hexaploid bread wheat (Triticum aestivum) encode for proteins with epitopes that can trigger coeliac disease (CD), and several contain a 33‐mer peptide with six partly overlapping copies of three epitopes, which is regarded as a remarkably potent T‐cell stimulator. To increase genetic diversity in the D genome, synthetic hexaploid wheat lines are being made by hybridising accessions of Triticum turgidum (AB genome) and Aegilops tauschii (the progenitor of the D genome). The diversity of alpha‐gliadins in A. tauschii has not been studied extensively. We analysed the alpha‐gliadin transcriptome of 51 A. tauschii accessions representative of the diversity in A. tauschii. We extracted RNA from developing seeds and performed 454 amplicon sequencing of the first part of the alpha‐gliadin genes. The expression profile of allelic variants of the alpha‐gliadins was different between accessions, and also between accessions of the Western and Eastern clades of A. tauschii. Generally, both clades expressed many allelic variants not found in bread wheat. In contrast to earlier studies, we detected the 33‐mer peptide in some A. tauschii accessions, indicating that it was introduced along with the D genome into bread wheat. In these accessions, transcripts with the 33‐mer peptide were present at lower frequencies than in bread wheat varieties. In most A. tauschii accessions, however, the alpha‐gliadins do not contain the epitope, and this may be exploited, through synthetic hexaploid wheats, to breed bread wheat varieties with fewer or no coeliac disease epitopes.

The role of gene flow and chromosomal instability in shaping the bread wheat genome.

Bread wheat (Triticum aestivum) is one of the world's most important crops; however, a low level of genetic diversity within commercial breeding accessions can significantly limit breeding potential. In contrast, wheat relatives exhibit considerable genetic variation and so potentially provide a valuable source of novel alleles for use in breeding new cultivars. Historically, gene flow between wheat and its relatives may have contributed novel alleles to the bread wheat pangenome. To assess the contribution made by wheat relatives to genetic diversity in bread wheat, we used markers based on single nucleotide polymorphisms to compare bread wheat accessions, created in the past 150 years, with 45 related species. We show that many bread wheat accessions share near-identical haplotype blocks with close relatives of wheat's diploid and tetraploid progenitors, while some show evidence of introgressions from more distant species and structural variation between accessions. Hence, introgressions and chromosomal rearrangements appear to have made a major contribution to genetic diversity in cultivar collections. As gene flow from relatives to bread wheat is an ongoing process, we assess the impact that introgressions might have on future breeding strategies.

Identification and development of novel salt-responsive candidate gene based SSRs (cg-SSRs) and MIR gene based SSRs (mir-SSRs) in bread wheat (Triticum aestivum)

Salt stress adversely affects the global wheat production and productivity. To improve salinity tolerance of crops, identification of robust molecular markers is highly imperative for development of salt-tolerant cultivars to mimic yield losses under saline conditions. In this study, we mined 171 salt-responsive genes (including 10 miRNAs) from bread wheat genome using the sequence information of functionally validated salt-responsive rice genes. Salt-stress, tissue and developmental stage-specific expression analysis of RNA-seq datasets revealed the constitutive as well as the inductive response of salt-responsive genes in different tissues of wheat. Fifty-four genotypes were phenotyped for salt stress tolerance. The stress tolerance index of the genotypes ranged from 0.30 to 3.18. In order to understand the genetic diversity, candidate gene based SSRs (cg-SSRs) and MIR gene based SSRs (miR-SSRs) were mined from 171 members of salt-responsive genes of wheat and validated among the contrasting panels of 54 tolerant as well as susceptible wheat genotypes. Among 53 SSR markers screened, 10 cg-SSRs and 8 miR-SSRs were found to be polymorphic. Polymorphic information content between the wheat genotypes ranged from 0.07 to 0.67, indicating the extant of wide genetic variation among the salt tolerant and susceptible genotypes at the DNA level. The genetic diversity analysis based on the allelic data grouped the wheat genotypes into three separate clusters of which single group encompassing most of the salt susceptible genotypes and two of them containing salt tolerance and moderately salt tolerance wheat genotypes were in congruence with penotypic data. Our study showed that both salt-responsive genes and miRNAs based SSRs were more diverse and can be effectively used for diversity analysis. This study reports the first extensive survey on genome-wide analysis, identification, development and validation of salt-responsive cg-SSRs and miR-SSRs in wheat. The information generated in the present study on genetic divergence among genotypes having a differential response to salt will help in the selection of suitable lines as parents for developing salt tolerant cultivars in wheat.

PaKman: A Scalable Algorithm for Generating Genomic Contigs on Distributed Memory Machines

De novo genome assembly is a fundamental problem in the field of bioinformatics, that aims to assemble the DNA sequence of an unknown genome from numerous short DNA fragments (aka reads) obtained from it. With the advent of high-throughput sequencing technologies, billions of reads can be generated in a matter of hours, necessitating efficient parallelization of the assembly process. While multiple parallel solutions have been proposed in the past, conducting a large-scale assembly at scale remains a challenging problem because of the inherent complexities associated with data movement, and irregular access footprints of memory and I/O operations. In this article, we present a novel algorithm, called PaKman , to address the problem of performing large-scale genome assemblies on a distributed memory parallel computer. Our approach focuses on improving performance through a combination of novel data structures and algorithmic strategies for reducing the communication and I/O footprint during the assembly process. PaKman presents a solution for the two most time-consuming phases in the full genome assembly pipeline, namely, k-mer counting and contig generation . A key aspect of our algorithm is its graph data structure (PaK-Graph), which comprises fat nodes (or what we call “macro-nodes”) that reduce the communication burden during contig generation. We present an extensive performance and qualitative evaluation of our algorithm across a wide range of genomes (varying in both size and species group), including comparisons to other state-of-the-art parallel assemblers. Our results demonstrate the ability to achieve near-linear speedups on up to 16K cores (tested) on the NERSC Cori supercomputer; perform better than or comparable to other state-of-the-art distributed memory and shared memory tools in terms of performance while delivering comparable (if not better) quality; and reduce time to solution significantly. For instance, PaKman is able to generate a high-quality set of assembled contigs for complex genomes such as the human and bread wheat genomes in under a minute on 16K cores. In addition, PaKman was able to successfully process a 3.1 TB simulated dataset of one of the largest known genomes (to date)- Ambystoma mexicanum (the axolotl), in just over 200 seconds on 16K cores.

Triticum population sequencing provides insights into wheat adaptation.

Bread wheat expanded its habitat from a core area of the Fertile Crescent to global environments within ~10,000 years. The genetic mechanisms of this remarkable evolutionary success are not well understood. By whole-genome sequencing of populations from 25 subspecies within the genera Triticum and Aegilops, we identified composite introgression from wild populations contributing to a substantial portion (4-32%) of the bread wheat genome, which increased the genetic diversity of bread wheat and allowed its divergent adaptation. Meanwhile, convergent adaptation to human selection showed 2- to 16-fold enrichment relative to random expectation-a certain set of genes were repeatedly selected in Triticum species despite their drastic differences in ploidy levels and growing zones, indicating the important role of evolutionary constraints in shaping the adaptive landscape of bread wheat. These results showed the genetic necessities of wheat as a global crop and provided new perspectives on transferring adaptive success across species for crop improvement.

QTL mapping and validation of bread wheat flag leaf morphology across multiple environments in different genetic backgrounds.

Eight major and stably expressed QTL for flag leaf morphology across eleven environments were identified and validated using newly developed KASP markers in seven biparental populations with different genetic backgrounds. Flag leaf morphology is a determinant trait influencing plant architecture and yield potential in wheat (Triticum aestivum L.). A recombinant inbred line (RIL) population with a 55K SNP-based constructed genetic map was used to map quantitative trait loci (QTL) for flag leaf length (FLL), width (FLW), area (FLA), angle (FLANG), opening angle (FLOA), and bend angle (FLBA) in eleven environments. Eight major QTL were detected in 11 environments with 5.73-54.38% of explained phenotypic variation. These QTL were successfully verified using the newly developed Kompetitive Allele Specific PCR (KASP) markers in six biparental populations with different genetic backgrounds. Among these 8 major QTL, two co-located intervals were identified. Significant interactions for both FLL- and FLW-related QTL were detected. Comparison analysis showed that QFll.sau-SY-2B and QFla.sau-SY-2B are likely new loci. Significant relationships between flag leaf- and yield-related traits were observed and discussed. Several genes associated with leaf development including the ortholog of maize ZmRAVL1, a B3-domain transcription factor involved in regulation of leaf angle, were predicted in physical intervals harboring these major QTL on reference genomes of bread wheat 'Chinese spring', T. turgidum, and Aegilops tauschii. Taken together, these results broaden our understanding on genetic basis of flag leaf morphology and provide clues for fine mapping and marker-assisted breeding wheat with optimized plant architecture for promising loci.

Chromosome-Scale Assembly of the Bread Wheat Genome Reveals Thousands of Additional Gene Copies.

Bread wheat (Triticum aestivum) is a major food crop and an important plant system for agricultural genetics research. However, due to the complexity and size of its allohexaploid genome, genomic resources are limited compared to other major crops. The IWGSC recently published a reference genome and associated annotation (IWGSC CS v1.0, Chinese Spring) that has been widely adopted and utilized by the wheat community. Although this reference assembly represents all three wheat subgenomes at chromosome-scale, it was derived from short reads, and thus is missing a substantial portion of the expected 16Gbp of genomic sequence. We earlier published an independent wheat assembly (Triticum_aestivum_3.1, Chinese Spring) that came much closer in length to the expected genome size, although it was only a contig-level assembly lacking gene annotations. Here, we describe a reference-guided effort to scaffold those contigs into chromosome-length pseudomolecules, add in any missing sequence that was unique to the IWGSC CS v1.0 assembly, and annotate the resulting pseudomolecules with genes. Our updated assembly, Triticum_aestivum_4.0, contains 15.07 Gbp of nongap sequence anchored to chromosomes, which is 1.2 Gbps more than the previous reference assembly. It includes 108,639 genes unambiguously localized to chromosomes, including over 2000 genes that were previously unplaced. We also discovered >5700 additional gene copies, facilitating the accurate annotation of functional gene duplications including at the Ppd-B1 photoperiod response locus.

Molecular characterization revealed the role of catalases under abiotic and arsenic stress in bread wheat (Triticum aestivum L.)

Catalases are crucial antioxidant enzymes that reduce the excessive level of H2O2 caused by various environmental stresses and metal toxicity and hence protect the plant cells. In this study, a total of ten TaCAT genes, forming three homeologous groups, were identified in the genome of bread wheat (Triticum aestivum L.) and named as per the wheat gene symbolization guidelines. The identified catalases were characterized for various structural and physicochemical features. The proximal active-site (F(D/A)RERIPERVVHAKGASA) and heme-ligand (R(I/V)F(S/A)Y(A/S)DTQ) signature motifs, catalytic residues and peroxisomal targeting peptides were found conserved. Phylogenetic analysis clustered TaCATs into three classes, which showed conserved functional specialization based on their tissue specific expression. Modulated spatio-temporal expression of various TaCAT genes and alteration in total catalase enzyme activity during heat, drought, salt and arsenic (AsIII and AsV) treatment suggested their roles in abiotic stress response and arsenic tolerance. Molecular cloning and overexpression of TaCAT3-B gene in Escherichia coli showed tolerance against heat, drought, salt and varied concentrations of AsIII and AsV treatments. The results further confirmed their role in stress tolerance and recommended that these genes can be used in future stress management strategies for the development of abiotic and arsenic stress resistant transgenic crops.

Waxy Gene-Orthologs in Wheat × Thinopyrum Amphidiploids

Starch, as the main component of grain in cereals, serves as the major source of calories in staple food and as a raw material for industry. As the technological and digestive properties of starch depend on its content, the management of its components, amylose and amylopectin, is of great importance. The starch properties of wheat grain can be attuned using allelic variations of genes, including granule-bound starch synthase I (GBSS I), or Wx. The tertiary gene pool, including wheatgrass (Thinopyrum) species, provides a wide spectrum of genes-orthologs that can be used to increase the allelic diversity of wheat genes by wide hybridization. Octaploid partial wheat–wheatgrass hybrids (amphidiploids, WWGHs) combine the complete genome of bread wheat (BBAADD), and a mixed genome from the chromosomes of intermediate wheatgrass (Thinopyrum intermedium, genomic composition JrJrJvsJvsStSt) and tall wheatgrass (Th. ponticum, JJJJJJJsJsJsJs). Thus, WWGHs may carry Wx genes not only of wheat (Wx-B1, Wx-A1 and Wx-D1) but also of wheatgrass origin. We aimed to assess the level of amylose in starch and investigate the polymorphism of Wx genes in 12 accessions of WWGHs. Additionally, we characterized orthologous Wx genes in the genomes of wild wheat-related species involved in the development of the studied WWGHs, Th. intermedium and Th. ponticum, as well as in the putative donors of their subgenomes, bessarabian wheatgrass (Th. bessarabicum, JbJb) and bluebunch wheatgrass (Pseudoroegneria stipifolia, St1St1St2St2). Although no significant differences in amylose content were found between different WWGH accessions, SDS-PAGE demonstrated that at least two WWGHs have an additional band. We sequenced the Wx gene-orthologs in Th. bessarabicum, P. stipifolia, Th. intermedium and Th. ponticum, and developed a WXTH marker that can discriminate the Thinopyrum Wx gene in the wheat background, and localized it to the 7E chromosome in Th. elongatum. Using the WXTH marker we revealed the allelic polymorphism of the Thinopyrum Wx gene in the studied WWGHs. The applicability of Thinopyrum Wx genes in wheat breeding and their effect on starch quality are discussed.

Identification and characterization of large-scale genomic rearrangements during wheat evolution.

Following allopolyploidization, nascent polyploid wheat species react with massive genomic rearrangements, including deletion of transposable element-containing sequences. While such massive rearrangements are considered to be a prominent process in wheat genome evolution and speciation, their structure, extent, and underlying mechanisms remain poorly understood. In this study, we retrieved ~3500 insertions of a specific variant of Fatima, one of the most dynamic gypsy long-terminal repeat retrotransposons in wheat from the recently available high-quality genome drafts of Triticum aestivum (bread wheat) and Triticum turgidum ssp. dicoccoides or wild emmer, the allotetraploid mother of all modern wheats. The dynamic nature of Fatima facilitated the identification of large (i.e., up to ~ 1 million bases) Fatima-containing insertions/deletions (indels) upon comparison of bread wheat and wild emmer genomes. We characterized 11 such indels using computer-assisted analysis followed by PCR validation, and found that they might have occurred via unequal intra-strand recombination or double-strand break (DSB) events. Additionally, we observed one case of introgression of novel DNA fragments from an unknown source into the wheat genome. Our data thus indicate that massive large-scale DNA rearrangements might play a prominent role in wheat speciation.

Understanding DNA Methylation Patterns in Wheat.

The bread wheat genome is large (17Gb), allohexaploid, and highly repetitive (80-90% of the genome), which makes genomic and epigenomic analyses expensive to conduct and a challenge to analyze. Here we provide an overview of recent bioinformatic and experimental methods that have been developed to understand DNA methylation patterns in the complex polyploid genome of wheat.

Powdery mildew disease resistance and marker-assisted screening at the Pm60 locus in wild diploid wheat Triticum urartu

Triticum urartu, a diploid wild wheat and progenitor of the A genome of bread wheat, is an important resource for resistance to powdery mildew fungus caused by Blumeria graminis f. sp. tritici (Bgt). In the present study we systematically characterized the interaction between the Bgt fungus and T. urartu at the microscopic level. We also tested 227 T. urartu accessions for reaction to Bgt isolate E09 and discovered previously uncategorized powdery mildew resistance in this collection. Pm60 is a CC-NB-LRR type powdery mildew resistance gene that has at least three functional alleles, Pm60, Pm60a, and Pm60b. A marker-assisted screen targeting the Pm60 locus identified a non-functional allele of Pm60a, designated as Pm60a′. A sequence comparison of Pm60a′ and Pm60a revealed that they differed by 58 SNPs and one 3-nucleotide deletion. Based on the sequence variations two molecular markers were developed to differentiate the functional Pm60a allele from the non-functional Pm60a′. Our screen revealed the presence of a previously uncharacterized powdery mildew resistance in T. urartu and provides new insights into the Pm60 locus. We believe that the two molecular markers developed here and new T. urartu resistant accessions will facilitate further identification of novel powdery mildew resistance genes and benefit breeding for powdery mildew resistance.

LTR_FINDER_parallel: parallelization of LTR_FINDER enabling rapid identification of long terminal repeat retrotransposons

Annotation of plant genomes is still a challenging task due to the abundance of repetitive sequences, especially long terminal repeat (LTR) retrotransposons. LTR_FINDER is a widely used program for the identification of LTR retrotransposons but its application on large genomes is hindered by its single-threaded processes. Here we report an accessory program that allows parallel operation of LTR_FINDER, resulting in up to 8500X faster identification of LTR elements. It takes only 72 min to process the 14.5 Gb bread wheat (Triticum aestivum) genome in comparison to 1.16 years required by the original sequential version. LTR_FINDER_parallel is freely available at https://github.com/oushujun/LTR_FINDER_parallel.

Genome-wide identification and characterization of the GDP-L-galactose phosphorylase gene family in bread wheat

BackgroundAscorbate is a powerful antioxidant in plants and an essential micronutrient for humans. The GDP-L-galactose phosphorylase (GGP) gene encodes the rate-limiting enzyme of the L-galactose pathway—the dominant ascorbate biosynthetic pathway in plants—and is a promising gene candidate for increasing ascorbate in crops. In addition to transcriptional regulation, GGP production is regulated at the translational level through an upstream open reading frame (uORF) in the long 5′-untranslated region (5’UTR). The GGP genes have yet to be identified in bread wheat (Triticum aestivum L.), one of the most important food grain sources for humans.ResultsBread wheat chromosomal groups 4 and 5 were found to each contain three homoeologous TaGGP genes on the A, B, and D subgenomes (TaGGP2-A/B/D and TaGGP1-A/B/D, respectively) and a highly conserved uORF was present in the long 5’UTR of all six genes. Phylogenetic analyses demonstrated that the TaGGP genes separate into two distinct groups and identified a duplication event of the GGP gene in the ancestor of the Brachypodium/Triticeae lineage. A microsynteny analysis revealed that the TaGGP1 and TaGGP2 subchromosomal regions have no shared synteny suggesting that TaGGP2 may have been duplicated via a transposable element. The two groups of TaGGP genes have distinct expression patterns with the TaGGP1 homoeologs broadly expressed across different tissues and developmental stages and the TaGGP2 homoeologs highly expressed in anthers. Transient transformation of the TaGGP coding sequences in Nicotiana benthamiana leaf tissue increased ascorbate concentrations more than five-fold, confirming their functional role in ascorbate biosynthesis in planta.ConclusionsWe have identified six TaGGP genes in the bread wheat genome, each with a highly conserved uORF. Phylogenetic and microsynteny analyses highlight that a transposable element may have been responsible for the duplication and specialized expression of GGP2 in anthers in the Brachypodium/Triticeae lineage. Transient transformation of the TaGGP coding sequences in N. benthamiana demonstrated their activity in planta. The six TaGGP genes and uORFs identified in this study provide a valuable genetic resource for increasing ascorbate concentrations in bread wheat.