Structural and functional partitioning of bread wheat chromosome 3B

Molecular identification of the wheat male fertility gene Ms1 and its prospects for hybrid breeding

The current rate of yield gain in crops is insufficient to meet the predicted demands. Capturing the yield boost from heterosis is one of the few technologies that offers rapid gain. Hybrids are widely used for cereals, maize and rice, but it has been a challenge to develop a viable hybrid system for bread wheat due to the wheat genome complexity, which is both large and hexaploid. Wheat is our most widely grown crop providing 20% of the calories for humans. Here, we describe the identification of Ms1, a gene proposed for use in large-scale, low-cost production of male-sterile (ms) female lines necessary for hybrid wheat seed production. We show that Ms1 completely restores fertility to ms1d, and encodes a glycosylphosphatidylinositol-anchored lipid transfer protein, necessary for pollen exine development. This represents a key step towards developing a robust hybridization platform in wheat.Heterosis can rapidly boost yield in crop species but development of hybrid-breeding systems for bread wheat remains a challenge. Here, Tucker et al. describe the molecular identification of the wheat Ms1 gene and discuss its potential for large-scale hybrid seed production in wheat.

RMI1 and TOP3α limit meiotic CO formation through their C-terminal domains

At meiosis, hundreds of programmed DNA double-strand breaks (DSBs) form and are repaired by homologous recombination. From this large number of DSBs, only a subset yields crossovers (COs), with a minimum of one CO per chromosome pair. All DSBs must be repaired and every recombination intermediate must be resolved to avoid subsequent entanglement and chromosome breakage. The conserved BLM-TOP3α-RMI1 (BTR) complex acts on early and late meiotic recombination intermediates to both limit CO outcome and promote chromosome integrity. In Arabidopsis, the BLM homologues RECQ4A and RECQ4B act redundantly to prevent meiotic extra COs, but recombination intermediates are fully resolved in their absence. In contrast, TOP3α is needed for both processes. Here we show through the characterization of specific mutants that RMI1 is a major anti-CO factor, in addition to being essential to prevent chromosome breakage and entanglement. Further, our findings suggest a specific role of the C-terminal domains of RMI1 and TOP3α, that respectively contain an Oligo Binding domain (OB2) and ZINC finger motifs, in preventing extra-CO. We propose that these domains of TOP3α and RMI1 define a sub-domain of the BTR complex which is dispensable for the resolution of recombination intermediates but crucial to limit extra-COs.

Mapping stripe rust resistance gene YrZH22 in Chinese wheat cultivar Zhoumai 22 by bulked segregant RNA-Seq (BSR-Seq) and comparative genomics analyses

A stripe rust resistance gene YrZH22 was mapped by combined BSR-Seq and comparative genomics analyses to a 5.92 centimorgan (cM) genetic interval spanning a 4 Mb physical genomic region on wheat chromosome 4BL1. Stripe rust, caused by Puccinia striiformis f. sp. tritici (PST), is one of the most destructive diseases of wheat and severely threatens wheat production worldwide. The widely grown Chinese wheat cultivar Zhoumai 22 is highly resistant to the current prevailing PST race CYR34 (V26). Genetic analysis of F_5:6 and F_6:7 recombinant inbred line (RIL) populations indicated that adult-plant stripe rust resistance in Zhoumai 22 is controlled by a single gene, temporarily designated YrZH22. By applying bulked segregant RNA-Seq (BSR-Seq), 7 SNP markers were developed and SNP mapping showed that YrZH22 is located between markers WGGB105 and WGGB112 on chromosome arm 4BL. The corresponding genomic regions of the Chinese Spring 4BL genome assembly and physical map of Aegilops tauschii 4DL were selected for comparative genomics analyses to develop nine new polymorphic markers that were used to construct a high-resolution genetic linkage map of YrZH22. YrZH22 was delimited in a 5.92 cM genetic interval between markers WGGB133 and WGGB146, corresponding to 4.1 Mb genomic interval in Chinese Spring 4BL and a 2.2 Mb orthologous genomic region in Ae. tauschii 4DL. The genetic linkage map of YrZH22 will be valuable for fine mapping and positional cloning of YrZH22, and can be used for marker-assisted selection in wheat breeding.

Genome-wide identification, classification, evolutionary analysis and gene expression patterns of the protein kinase gene family in wheat and Aegilops tauschii

In this study we systematically identified and classified PKs in Triticum aestivum, Triticum urartu and Aegilops tauschii. Domain distribution and exon-intron structure analyses of PKs were performed, and we found conserved exon-intron structures within the exon phases in the kinase domain. Collinearity events were determined, and we identified various T. aestivum PKs from polyploidizations and tandem duplication events. Global expression pattern analysis of T. aestivum PKs revealed that some PKs might participate in the signaling pathways of stress response and developmental processes. QRT-PCR of 15 selected PKs were performed under drought treatment and with infection of Fusarium graminearum to validate the prediction of microarray. The protein kinase (PK) gene superfamily is one of the largest families in plants and participates in various plant processes, including growth, development, and stress response. To better understand wheat PKs, we conducted genome-wide identification, classification, evolutionary analysis and expression profiles of wheat and Ae. tauschii PKs. We identified 3269, 1213 and 1448 typical PK genes in T. aestivum, T. urartu and Ae. tauschii, respectively, and classified them into major groups and subfamilies. Domain distributions and gene structures were analyzed and visualized. Some conserved intron-exon structures within the conserved kinase domain were found in T. aestivum, T. urartu and Ae. tauschii, as well as the primitive land plants Selaginella moellendorffii and Physcomitrella patens, revealing the important roles and conserved evolutionary history of these PKs. We analyzed the collinearity events of T. aestivum PKs and identified PKs from polyploidizations and tandem duplication events. Global expression pattern analysis of T. aestivum PKs revealed tissue-specific and stress-specific expression profiles, hinting that some wheat PKs may regulate abiotic and biotic stress response signaling pathways. QRT-PCR of 15 selected PKs were performed under drought treatment and with infection of F. graminearum to validate the prediction of microarray. Our results will provide the foundational information for further studies on the molecular functions of wheat PKs.

Haplotype-resolved sweet potato genome traces back its hexaploidization history

Here we present the 15 pseudochromosomes of sweet potato, Ipomoea batatas, the seventh most important crop in the world and the fourth most significant in China. By using a novel haplotyping method based on genome assembly, we have produced a half haplotype-resolved genome from ~296 Gb of paired-end sequence reads amounting to roughly 67-fold coverage. By phylogenetic tree analysis of homologous chromosomes, it was possible to estimate the time of two recent whole-genome duplication events as occurring about 0.8 and 0.5 million years ago. This half haplotype-resolved hexaploid genome represents the first successful attempt to investigate the complexity of chromosome sequence composition directly in a polyploid genome, using sequencing of the polyploid organism itself rather than any of its simplified proxy relatives. Adaptation and application of our approach should provide higher resolution in future genomic structure investigations, especially for similarly complex genomes.

Cysteine proteases and wheat (Triticum aestivum L) under drought: A still greatly unexplored association

Bread wheat (Triticum aestivum L.) provides about 19% of global dietary energy. Environmental stress, such as drought, affects wheat growth causing premature plant senescence and ultimately plant death. A plant response to drought is an increase in protease-mediated proteolysis with rapid degradation of proteins required for metabolic processes. Among the plant proteases that are increased in their activity following stress, cysteine proteases are the best characterized. Very little is known about particular wheat cysteine protease sequences, their expression and also localization. The current knowledge on wheat cysteine proteases belonging to the five clans (CA, CD, CE, CF and CP) is outlined, in particular their expression and possible function under drought. The first successes in establishing an annotated wheat genome database are further highlighted which has allowed more detailed mining of cysteine proteases. We also share our thoughts on future research directions considering the growing availability of genomic resources of this very important food crop. Finally, we also outline future application of developed knowledge in transgenic wheat plants for environmental stress protection and also as senescence markers to monitor wheat growth under environmental stress conditions.

Transposable Elements Direct The Coevolution between Plants and Microbes

Transposable elements are powerful drivers of genome evolution in many eukaryotes. Although they are mostly considered as 'selfish' genetic elements, increasing evidence suggests that they contribute to genetic variability; particularly under stress conditions. Over the past few years, the role of transposable elements during host-microbe interactions has been recognised. It has been proposed that many pathogenic microbes have evolved a 'two-speed' genome with regions that show increased variability and that are enriched in transposable elements and pathogenicity-related genes. Plants similarly display structured genomes with transposable-element-rich regions that mediate accelerated evolution. Immune receptor genes typically reside in such regions. Various mechanisms have recently been identified through which transposable elements contribute to the coevolution between plants and their associated microbes.

In Silico Analysis of Fatty Acid Desaturase Genes and Proteins in Grasses

Fatty acid desaturases (FADs) catalyze the introduction of a double bond into acyl chains. Two FAD groups have been identified in plants: acyl-acyl carrier proteins (ACPs) and acyl-lipid or membrane-bound FAD. The former catalyze the conversion of 18:0 to 18:1 and to date have only been identified in plants. The latter are found in eukaryotes and bacteria and are responsible for multiple desaturations. In this study, we identified 82 desaturase gene and protein sequences from 10 grass species deposited in GenBank that were analyzed using bioinformatic approaches. Subcellular localization predictions of desaturase family revealed their localization in plasma membranes, chloroplasts, endoplasmic reticula, and mitochondria. The in silico mapping showed multiple chromosomal locations in most species. Furthermore, the presence of the characteristic histidine domains, the predicted motifs, and the finding of transmembrane regions strongly support the protein functionality. The identification of putative regulatory sites in the promotor and the expression profiles revealed the wide range of pathways in which fatty acid desaturases are involved. This study is an updated survey on desaturases of grasses that provides a comprehensive insight into diversity and evolution. This characterization is a necessary first step before considering these genes as candidates for new biotechnological approaches.

Genome-wide identification and analysis of the MADS-box gene family in bread wheat (Triticum aestivum L.)

The MADS-box genes encode transcription factors with key roles in plant growth and development. A comprehensive analysis of the MADS-box gene family in bread wheat (Triticum aestivum) has not yet been conducted, and our understanding of their roles in stress is rather limited. Here, we report the identification and characterization of the MADS-box gene family in wheat. A total of 180 MADS-box genes classified as 32 Mα, 5 Mγ, 5 Mδ, and 138 MIKC types were identified. Evolutionary analysis of the orthologs among T. urartu, Aegilops tauschii and wheat as well as homeologous sequences analysis among the three sub-genomes in wheat revealed that gene loss and chromosomal rearrangements occurred during and/or after the origin of bread wheat. Forty wheat MADS-box genes that were expressed throughout the investigated tissues and development stages were identified. The genes that were regulated in response to both abiotic stresses (i.e., phosphorus deficiency, drought, heat, and combined drought and heat) and biotic stresses (i.e., Fusarium graminearum, Septoria tritici, stripe rust and powdery mildew) were detected as well. A few notable MADS-box genes were specifically expressed in a single tissue and those showed relatively higher expression differences between the stress and control treatment. The expression patterns of considerable MADS-box genes differed from those of their orthologs in Brachypodium, rice, and Arabidopsis. Collectively, the present study provides new insights into the possible roles of MADS-box genes in response to stresses and will be valuable for further functional studies of important candidate MADS-box genes.

Spatiotemporal expression patterns of wheat amino acid transporters reveal their putative roles in nitrogen transport and responses to abiotic stress

Amino acid transporters have roles in amino acid uptake from soil, long-distance transport, remobilization from vegetative tissues and accumulation in grain. Critically, the majority of wheat grain nitrogen is derived from amino acids remobilized from vegetative organs. However, no systematic analysis of wheat AAT genes has been reported to date. Here, 283 full length wheat AAT genes representing 100 distinct groups of homeologs were identified and curated by selectively consolidating IWGSC CSSv2 and TGACv1 Triticum aestivum genome assemblies and reassembling or mapping of IWGSC CSS chromosome sorted reads to fill any gaps. Gene expression profiling was performed using public RNA-seq data from root, leaf, stem, spike, grain and grain cells (transfer cell (TC), aleurone cell (AL), and starchy endosperm (SE)). AATs highly expressed in roots are good candidates for amino acid uptake from soil whilst AATs highly expressed in senescing leaves and stems may be involved in translocation to grain. AATs in TC (TaAAP2 and TaAAP19) and SE (TaAAP13) may play important roles in determining grain protein content and grain yield. The expression levels of AAT homeologs showed unequal contributions in response to abiotic stresses and development, which may aid wheat adaptation to a wide range of environments.

Spatiotemporal expression patterns of wheat amino acid transporters reveal their putative roles in nitrogen transport and responses to abiotic stress

Amino acid transporters have roles in amino acid uptake from soil, long-distance transport, remobilization from vegetative tissues and accumulation in grain. Critically, the majority of wheat grain nitrogen is derived from amino acids remobilized from vegetative organs. However, no systematic analysis of wheat AAT genes has been reported to date. Here, 283 full length wheat AAT genes representing 100 distinct groups of homeologs were identified and curated by selectively consolidating IWGSC CSSv2 and TGACv1 Triticum aestivum genome assemblies and reassembling or mapping of IWGSC CSS chromosome sorted reads to fill any gaps. Gene expression profiling was performed using public RNA-seq data from root, leaf, stem, spike, grain and grain cells (transfer cell (TC), aleurone cell (AL), and starchy endosperm (SE)). AATs highly expressed in roots are good candidates for amino acid uptake from soil whilst AATs highly expressed in senescing leaves and stems may be involved in translocation to grain. AATs in TC (TaAAP2 and TaAAP19) and SE (TaAAP13) may play important roles in determining grain protein content and grain yield. The expression levels of AAT homeologs showed unequal contributions in response to abiotic stresses and development, which may aid wheat adaptation to a wide range of environments.

Exome sequence genotype imputation in globally diverse hexaploid wheat accessions

Imputing genotypes from the 90K SNP chip to exome sequence in wheat was moderately accurate. We investigated the factors that affect imputation and propose several strategies to improve accuracy. Imputing genetic marker genotypes from low to high density has been proposed as a cost-effective strategy to increase the power of downstream analyses (e.g. genome-wide association studies and genomic prediction) for a given budget. However, imputation is often imperfect and its accuracy depends on several factors. Here, we investigate the effects of reference population selection algorithms, marker density and imputation algorithms (Beagle4 and FImpute) on the accuracy of imputation from low SNP density (9K array) to the Infinium 90K single-nucleotide polymorphism (SNP) array for a collection of 837 hexaploid wheat Watkins landrace accessions. Based on these results, we then used the best performing reference selection and imputation algorithms to investigate imputation from 90K to exome sequence for a collection of 246 globally diverse wheat accessions. Accession-to-nearest-entry and genomic relationship-based methods were the best performing selection algorithms, and FImpute resulted in higher accuracy and was more efficient than Beagle4. The accuracy of imputing exome capture SNPs was comparable to imputing from 9 to 90K at approximately 0.71. This relatively low imputation accuracy is in part due to inconsistency between 90K and exome sequence formats. We also found the accuracy of imputation could be substantially improved to 0.82 when choosing an equivalent number of exome SNP, instead of 90K SNPs on the existing array, as the lower density set. We present a number of recommendations to increase the accuracy of exome imputation.

Physical Map of the Short Arm of Bread Wheat Chromosome 3D

Bread wheat ( L.) is one of the most important crops worldwide. Although a reference genome sequence would represent a valuable resource for wheat improvement through genomics-assisted breeding and gene cloning, its generation has long been hampered by its allohexaploidy, high repeat content, and large size. As a part of a project coordinated by the International Wheat Genome Sequencing Consortium (IWGSC), a physical map of the short arm of wheat chromosome 3D (3DS) was prepared to facilitate reference genome assembly and positional gene cloning. It comprises 869 contigs with a cumulative length of 274.5 Mbp and represents 85.5% of the estimated chromosome arm size. Eighty-six Mbp of survey sequences from chromosome arm 3DS were assigned in silico to physical map contigs via next-generation sequencing of bacterial artificial chromosome pools, thus providing a high-density framework for physical map ordering along the chromosome arm. About 60% of the physical map was anchored in this single experiment. Finally, 1393 high-confidence genes were anchored to the physical map. Comparisons of gene space of the chromosome arm 3DS with genomes of closely related species [ (L.) P.Beauv., rice ( L.), and sorghum [ (L.) Moench] and homeologous wheat chromosomes provided information about gene movement on the chromosome arm.

Comparative Proteomic Analysis of Flag Leaves Reveals New Insight into Wheat Heat Adaptation

Hexaploid wheat (Triticum aestivum L.) is an important food crop but it is vulnerable to heat. The heat-responsive proteome of wheat remains to be fully elucidated because of previous technical and genomic limitations, and this has hindered our understanding of the mechanisms of wheat heat adaptation and advances in improving thermotolerance. Here, flag leaves of wheat during grain filling stage were subjected to high daytime temperature stress, and 258 heat-responsive proteins (HRPs) were identified with iTRAQ analysis. Enrichment analysis revealed that chlorophyll synthesis, carbon fixation, protein turnover, and redox regulation were the most remarkable heat-responsive processes. The HRPs involved in chlorophyll synthesis and carbon fixation were significantly decreased, together with severe membrane damage, demonstrating the specific effects of heat on photosynthesis of wheat leaves. In addition, the decrease in chlorophyll content may result from the decrease in HRPs involved in chlorophyll precursor synthesis. Further analysis showed that the accumulated effect of heat stress played a critical role in photosynthesis reduction, suggested that improvement in heat tolerance of photosynthesis, and extending heat tolerant period would be major research targets. The significantly accumulation of GSTs and Trxs in response to heat suggested their important roles in redox regulation, and they could be the promising candidates for improving wheat thermotolerance. In summary, our results provide new insight into wheat heat adaption and provide new perspectives on thermotolerance improvement.

Prediction and analysis of three gene families related to leaf rust (Puccinia triticina) resistance in wheat (Triticum aestivum L.)

BACKGROUND

The resistance to leaf rust (Lr) caused by Puccinia triticina in wheat (Triticum aestivum L.) has been well studied over the past decades with over 70 Lr genes being mapped on different chromosomes and numerous QTLs (quantitative trait loci) being detected or mapped using DNA markers. Such resistance is often divided into race-specific and race-nonspecific resistance. The race-nonspecific resistance can be further divided into resistance to most or all races of the same pathogen and resistance to multiple pathogens. At the molecular level, these three types of resistance may cover across the whole spectrum of pathogen specificities that are controlled by genes encoding different protein families in wheat. The objective of this study is to predict and analyze genes in three such families: NBS-LRR (nucleotide-binding sites and leucine-rich repeats or NLR), START (Steroidogenic Acute Regulatory protein [STaR] related lipid-transfer) and ABC (ATP-Binding Cassette) transporter. The focus of the analysis is on the patterns of relationships between these protein-coding genes within the gene families and QTLs detected for leaf rust resistance.

RESULTS

We predicted 526 ABC, 1117 NLR and 144 START genes in the hexaploid wheat genome through a domain analysis of wheat proteome. Of the 1809 SNPs from leaf rust resistance QTLs in seedling and adult stages of wheat, 126 SNPs were found within coding regions of these genes or their neighborhood (5 Kb upstream from transcription start site [TSS] or downstream from transcription termination site [TTS] of the genes). Forty-three of these SNPs for adult resistance and 18 SNPs for seedling resistance reside within coding or neighboring regions of the ABC genes whereas 14 SNPs for adult resistance and 29 SNPs for seedling resistance reside within coding or neighboring regions of the NLR gene. Moreover, we found 17 nonsynonymous SNPs for adult resistance and five SNPs for seedling resistance in the ABC genes, and five nonsynonymous SNPs for adult resistance and six SNPs for seedling resistance in the NLR genes. Most of these coding SNPs were predicted to alter encoded amino acids and such information may serve as a starting point towards more thorough molecular and functional characterization of the designated Lr genes. Using the primer sequences of 99 known non-SNP markers from leaf rust resistance QTLs, we found candidate genes closely linked to these markers, including Lr34 with distances to its two gene-specific markers being 1212 bases (to cssfr1) and 2189 bases (to cssfr2).

CONCLUSION

This study represents a comprehensive analysis of ABC, NLR and START genes in the hexaploid wheat genome and their physical relationships with QTLs for leaf rust resistance at seedling and adult stages. Our analysis suggests that the ABC (and START) genes are more likely to be co-located with QTLs for race-nonspecific, adult resistance whereas the NLR genes are more likely to be co-located with QTLs for race-specific resistance that would be often expressed at the seedling stage. Though our analysis was hampered by inaccurate or unknown physical positions of numerous QTLs due to the incomplete assembly of the complex hexaploid wheat genome that is currently available, the observed associations between (i) QTLs for race-specific resistance and NLR genes and (ii) QTLs for nonspecific resistance and ABC genes will help discover SNP variants for leaf rust resistance at seedling and adult stages. The genes containing nonsynonymous SNPs are promising candidates that can be investigated in future studies as potential new sources of leaf rust resistance in wheat breeding.

Characterization, Genomic Organization, Abundance, and Chromosomal Distribution of Ty1-copia Retrotransposons in Erianthus arundinaceus

Erianthus arundinaceus is an important wild species of the genus Saccharum with many valuable traits. However, the composition and structure of its genome are largely unknown, which have hindered its utilization in sugarcane breeding and evolutionary research. Retrotransposons constitute an appreciable fraction of plant genomes and may have played a significant role in the evolution and sequence organization of genomes. In the current study, we investigate the phylogenetic diversity and genomic abundance of Ty1-copia retrotransposons for the first time and inspect their chromosomal distribution patterns in E. arundinaceus. In total, 70 Ty1-copia reverse transcriptase (RT) sequences with significant levels of heterogeneity were obtained. The phylogenetic analysis revealed these Ty1-copia retrotransposons were classified into four distinct evolutionary lineages (Tork/TAR, Tork/Angela, Retrofit/Ale, and Sire/Maximus). Dot-blot analysis showed estimated the total copy number of Ty1-copia retrotransposons to be about 4.5 × 10³ in the E. arundinaceus genome, indicating they were a significant component. Fluorescence in situ hybridization revealed that Ty1-copia retrotransposons from the four lineages had strikingly similar patterns of chromosomal enrichment, being exclusively enriched in the subterminal heterochromatic regions of most E. arundinaceus chromosomes. This is the first clear evidence of the presence of Ty1-copia retrotransposons in the subterminal heterochromatin of E. arundinaceus. Altogether, these results promote the understanding of the diversification of Ty1-copia retrotransposons and shed light on their chromosomal distribution patterns in E. arundinaceus.

Massive expansion and differential evolution of small heat shock proteins with wheat (Triticum aestivum L.) polyploidization

Wheat (Triticum aestivum), one of the world’s most important crops, is facing unprecedented challenges due to global warming. To evaluate the gene resources for heat adaptation in hexaploid wheat, small heat shock proteins (sHSPs), the key plant heat protection genes, were comprehensively analysed in wheat and related species. We found that the sHSPs of hexaploid wheat were massively expanded in A and B subgenomes with intrachromosomal duplications during polyploidization. These expanded sHSPs were under similar purifying selection and kept the expressional patterns with the original copies. Generally, a strong purifying selection acted on the α-crystallin domain (ACD) and theoretically constrain conserved function. Meanwhile, weaker purifying selection and strong positive selection acted on the N-terminal region, which conferred sHSP flexibility, allowing adjustments to a wider range of substrates in response to genomic and environmental changes. Notably, in CI, CV, ER, MI and MII subfamilies, gene duplications, expression variations and functional divergence occurred before wheat polyploidization. Our results indicate the massive expansion of active sHSPs in hexaploid wheat may also provide more raw materials for evolving functional novelties and generating genetic diversity to face future global climate changes, and highlight the expansion of stress response genes with wheat polyploidization.

High-Resolution Mapping of Crossover Events in the Hexaploid Wheat Genome Suggests a Universal Recombination Mechanism

During meiosis, crossovers (COs) create new allele associations by reciprocal exchange of DNA. In bread wheat (Triticum aestivum L.), COs are mostly limited to subtelomeric regions of chromosomes, resulting in a substantial loss of breeding efficiency in the proximal regions, though these regions carry ∼60-70% of the genes. Identifying sequence and/or chromosome features affecting recombination occurrence is thus relevant to improve and drive recombination. Using the recent release of a reference sequence of chromosome 3B and of the draft assemblies of the 20 other wheat chromosomes, we performed fine-scale mapping of COs and revealed that 82% of COs located in the distal ends of chromosome 3B representing 19% of the chromosome length. We used 774 SNPs to genotype 180 varieties representative of the Asian and European genetic pools and a segregating population of 1270 F₆ lines. We observed a common location for ancestral COs (predicted through linkage disequilibrium) and the COs derived from the segregating population. We delineated 73 small intervals (<26 kb) on chromosome 3B that contained 252 COs. We observed a significant association of COs with genic features (73 and 54% in recombinant and nonrecombinant intervals, respectively) and with those expressed during meiosis (67% in recombinant intervals and 48% in nonrecombinant intervals). Moreover, while the recombinant intervals contained similar amounts of retrotransposons and DNA transposons (42 and 53%), nonrecombinant intervals had a higher level of retrotransposons (63%) and lower levels of DNA transposons (28%). Consistent with this, we observed a higher frequency of a DNA motif specific to the TIR-Mariner DNA transposon in recombinant intervals.

Rapid cloning of genes in hexaploid wheat using cultivar-specific long-range chromosome assembly

Cereal crops such as wheat and maize have large repeat-rich genomes that make cloning of individual genes challenging. Moreover, gene order and gene sequences often differ substantially between cultivars of the same crop species. A major bottleneck for gene cloning in cereals is the generation of high-quality sequence information from a cultivar of interest. In order to accelerate gene cloning from any cropping line, we report 'targeted chromosome-based cloning via long-range assembly' (TACCA). TACCA combines lossless genome-complexity reduction via chromosome flow sorting with Chicago long-range linkage to assemble complex genomes. We applied TACCA to produce a high-quality (N50 of 9.76 Mb) de novo chromosome assembly of the wheat line CH Campala Lr22a in only 4 months. Using this assembly we cloned the broad-spectrum Lr22a leaf-rust resistance gene, using molecular marker information and ethyl methanesulfonate (EMS) mutants, and found that Lr22a encodes an intracellular immune receptor homologous to the Arabidopsis thaliana RPM1 protein.

A TRIM insertion in the promoter of Ms2 causes male sterility in wheat

The male-sterile ms2 mutant has been known for 40 years and has become extremely important in the commercial production of wheat. However, the gene responsible for this phenotype has remained unknown. Here we report the map-based cloning of the Ms2 gene. The Ms2 locus is remarkable in several ways that have implications in basic biology. Beyond having no functional annotation, barely detectable transcription in fertile wild-type wheat plants, and accumulated destructive mutations in Ms2 orthologs, the Ms2 allele in the ms2 mutant has acquired a terminal-repeat retrotransposon in miniature (TRIM) element in its promoter. This TRIM element is responsible for the anther-specific Ms2 activation that confers male sterility. The identification of Ms2 not only unravels the genetic basis of a historically important breeding trait, but also shows an example of how a TRIM element insertion near a gene can contribute to genetic novelty and phenotypic plasticity.

Amplifying recombination genome-wide and reshaping crossover landscapes in Brassicas

Meiotic recombination by crossovers (COs) is tightly regulated, limiting its key role in producing genetic diversity. However, while COs are usually restricted in number and not homogenously distributed along chromosomes, we show here how to disrupt these rules in Brassica species by using allotriploid hybrids (AAC, 2n = 3x = 29), resulting from the cross between the allotetraploid rapeseed (B. napus, AACC, 2n = 4x = 38) and one of its diploid progenitors (B. rapa, AA, 2n = 2x = 20). We produced mapping populations from different genotypes of both diploid AA and triploid AAC hybrids, used as female and/or as male. Each population revealed nearly 3,000 COs that we studied with SNP markers well distributed along the A genome (on average 1 SNP per 1.25 Mbp). Compared to the case of diploids, allotriploid hybrids showed 1.7 to 3.4 times more overall COs depending on the sex of meiosis and the genetic background. Most surprisingly, we found that such a rise was always associated with (i) dramatic changes in the shape of recombination landscapes and (ii) a strong decrease of CO interference. Hybrids carrying an additional C genome exhibited COs all along the A chromosomes, even in the vicinity of centromeres that are deprived of COs in diploids as well as in most studied species. Moreover, in male allotriploid hybrids we found that Class I COs are mostly responsible for the changes of CO rates, landscapes and interference. These results offer the opportunity for geneticists and plant breeders to dramatically enhance the generation of diversity in Brassica species by disrupting the linkage drag coming from limits on number and distribution of COs.

An improved assembly and annotation of the allohexaploid wheat genome identifies complete families of agronomic genes and provides genomic evidence for chromosomal translocations

Advances in genome sequencing and assembly technologies are generating many high-quality genome sequences, but assemblies of large, repeat-rich polyploid genomes, such as that of bread wheat, remain fragmented and incomplete. We have generated a new wheat whole-genome shotgun sequence assembly using a combination of optimized data types and an assembly algorithm designed to deal with large and complex genomes. The new assembly represents >78% of the genome with a scaffold N50 of 88.8 kb that has a high fidelity to the input data. Our new annotation combines strand-specific Illumina RNA-seq and Pacific Biosciences (PacBio) full-length cDNAs to identify 104,091 high-confidence protein-coding genes and 10,156 noncoding RNA genes. We confirmed three known and identified one novel genome rearrangements. Our approach enables the rapid and scalable assembly of wheat genomes, the identification of structural variants, and the definition of complete gene models, all powerful resources for trait analysis and breeding of this key global crop.

Mapping of common bunt resistance gene Bt9 in wheat

KEY MESSAGE

The Bt9 resistance locus was mapped and shown to be distinct from the Bt10 locus. New markers linked to Bt9 have been identified and may be used to breed for resistance towards the seed-borne disease. Increasing organic wheat production in Denmark, and in other wheat-producing areas, in conjunction with legal requirements for organic seed production, may potentially lead to a rise in common bunt occurrence. As systemic pesticides are not used in organic farming, organic wheat production systems may benefit from genetic resistances. However, little is known about the underlying genetic mechanisms and locations of the resistance factors for common bunt resistance in wheat. A double haploid (DH) population segregating for common bunt resistance was used to identify the chromosomal location of common bunt resistance gene Bt9. DH lines were phenotyped in three environments and genotyped with DArTseq and SSR markers. The total length of the resulting linkage map was 2882 cM distributed across all 21 wheat chromosomes. Bt9 was mapped to the distal end of chromosome 6DL. Since wheat common bunt resistance gene Bt10 is also located on chromosome 6D, the possibility of their co-location was investigated. A comparison of marker sequences linked to Bt9 and Bt10 on physical maps of chromosome 6D confirmed that Bt9 and Bt10 are two distinct resistance factors located at the distal (6DL) and proximal (6DS) end, respectively, of chromosome 6D. Five new SSR markers Xgpw4005-1, Xgpw7433, Xwmc773, Xgpw7303 and Xgpw362 and many SNP and PAV markers flanking the Bt9 resistance locus were identified and they may be used in the future for marker-assisted selection.

microRNA172 plays a crucial role in wheat spike morphogenesis and grain threshability

Wheat domestication from wild species involved mutations in the Q gene. The q allele (wild wheats) is associated with elongated spikes and hulled grains, whereas the mutant Q allele (domesticated wheats) confers subcompact spikes and free-threshing grains. Previous studies showed that Q encodes an AP2-like transcription factor, but the causal polymorphism of the domestication traits remained unclear. Here, we show that the interaction between microRNA172 (miR172) and the Q allele is reduced by a single nucleotide polymorphism in the miRNA binding site. Inhibition of miR172 activity by a miRNA target mimic resulted in compact spikes and transition from glumes to florets in apical spikelets. By contrast, overexpression of miR172 was sufficient to induce elongated spikes and non-free-threshing grains, similar to those observed in three Q loss-of-function mutations. These lines showed transitions from florets to glumes in the basal spikelets. These localized homeotic changes were associated with opposing miR172/Q gradients along the spike. We propose that the selection of a nucleotide change at the miR172 binding site of Q contributed to subcompact spikes and free-threshing grains during wheat domestication.

Highlighted Article: A nucleotide change in the microRNA172 binding site of the AP2 gene Q played a critical role in wheat domestication and the origin of free-threshing modern wheats. See also Greenwood et al. in this issue.

Evolution of Gene Expression Balance Among Homeologs of Natural Polyploids

Polyploidy is a major evolutionary process in eukaryotes, yet the expression balance of homeologs in natural polyploids is largely unknown. To study this expression balance, the expression patterns of 2180 structurally well-characterized genes of wheat were studied, of which 813 had the expected three copies and 375 had less than three. Copy numbers of the remaining 992 ranged from 4 to 14, including homeologs, orthologs, and paralogs. Of the genes with three structural copies corresponding to homeologs, 55% expressed from all three, 38% from two, and the remaining 7% expressed from only one of the three copies. Homeologs of 76-87% of the genes showed differential expression patterns in different tissues, thus have evolved different gene expression controls, possibly resulting in novel functions. Homeologs of 55% of the genes showed tissue-specific expression, with the largest percentage (14%) in the anthers and the smallest (7%) in the pistils. The highest number (1.72/3) of homeologs/gene expression was in the roots and the lowest (1.03/3) in the anthers. As the expression of homeologs changed with changes in structural copy number, about 30% of the genes showed dosage dependence. Chromosomal location also impacted expression pattern as a significantly higher proportion of genes in the proximal regions showed expression from all three copies compared to that present in the distal regions.

iTRAQ-based quantitative proteomic analysis of wheat roots in response to salt stress

Salinity is a major abiotic stress that affects plant growth and development. Plant roots are the sites of salt uptake. Here, an isobaric tag for a relative and absolute quantitation based proteomic technique was employed to identify the differentially expressed proteins (DEPs) from seedling roots of the salt-tolerant genotype Han 12 and the salt-sensitive genotype Jimai 19 in response to salt treatment. A total of 121 NaCl-responsive DEPs were observed in Han 12 and Jimai 19. The main DEPs were ubiquitination-related proteins, transcription factors, pathogen-related proteins, membrane intrinsic protein transporters and antioxidant enzymes, which may work together to obtain cellular homeostasis in roots and to determine the overall salt tolerance of different wheat varieties in response to salt stress. Functional analysis of three salt-responsive proteins was performed in transgenic plants as a case study to confirm the salt-related functions of the detected proteins. Taken together, the results of this study may be helpful in further elucidating salt tolerance mechanisms in wheat.

A high-resolution reference genetic map positioning 8.8 K genes for the conifer white spruce: structural genomics implications and correspondence with physical distance

Over the last decade, extensive genetic and genomic resources have been developed for the conifer white spruce (Picea glauca, Pinaceae), which has one of the largest plant genomes (20 Gbp). Draft genome sequences of white spruce and other conifers have recently been produced, but dense genetic maps are needed to comprehend genome macrostructure, delineate regions involved in quantitative traits, complement functional genomic investigations, and assist the assembly of fragmented genomic sequences. A greatly expanded P. glauca composite linkage map was generated from a set of 1976 full-sib progeny, with the positioning of 8793 expressed genes. Regions with significant low or high gene density were identified. Gene family members tended to be mapped on the same chromosomes, with tandemly arrayed genes significantly biased towards specific functional classes. The map was integrated with transcriptome data surveyed across eight tissues. In total, 69 clusters of co-expressed and co-localising genes were identified. A high level of synteny was found with pine genetic maps, which should facilitate the transfer of structural information in the Pinaceae. Although the current white spruce genome sequence remains highly fragmented, dozens of scaffolds encompassing more than one mapped gene were identified. From these, the relationship between genetic and physical distances was examined and the genome-wide recombination rate was found to be much smaller than most estimates reported for angiosperm genomes. This gene linkage map shall assist the large-scale assembly of the next-generation white spruce genome sequence and provide a reference resource for the conifer genomics community.

Genetic and epigenetic variation of transposable elements in Arabidopsis

Transposable elements are mobile genetic elements that are prevalent in plant genomes and are silenced by epigenetic modification. Different epigenetic modification pathways play distinct roles in the control of transposable element transcription, replication and recombination. The Arabidopsis genome contains families of all of the major transposable element classes, which are differentially enriched in particular genomic regions. Whole genome sequencing and DNA methylation profiling of hundreds of natural Arabidopsis accessions has revealed that transposable elements exhibit significant intraspecific genetic and epigenetic variation, and that genetic variation often underlies epigenetic variation. Together, epigenetic modification and the forces of selection define the scope within which transposable elements can contribute to, and control, genome evolution.

mlo-based powdery mildew resistance in hexaploid bread wheat generated by a non-transgenic TILLING approach

Wheat is one of the most widely grown cereal crops in the world and is an important food grain source for humans. However, wheat yields can be reduced by many abiotic and biotic stress factors, including powdery mildew disease caused by Blumeria graminis f.sp. tritici (Bgt). Generating resistant varieties is thus a major effort in plant breeding. Here, we took advantage of the non-transgenic Targeting Induced Lesions IN Genomes (TILLING) technology to select partial loss-of-function alleles of TaMlo, the orthologue of the barley Mlo (Mildew resistance locus o) gene. Natural and induced loss-of-function alleles (mlo) of barley Mlo are known to confer durable broad-spectrum powdery mildew resistance, typically at the expense of pleiotropic phenotypes such as premature leaf senescence. We identified 16 missense mutations in the three wheat TaMlo homoeologues, TaMlo-A1, TaMlo-B1 and TaMlo-D1 that each lead to single amino acid exchanges. Using transient gene expression assays in barley single cells, we functionally analysed the different missense mutants and identified the most promising candidates affecting powdery mildew susceptibility. By stacking of selected mutant alleles we generated four independent lines with non-conservative mutations in each of the three TaMlo homoeologues. Homozygous triple mutant lines and surprisingly also some of the homozygous double mutant lines showed enhanced, yet incomplete, Bgt resistance without the occurrence of discernible pleiotropic phenotypes. These lines thus represent an important step towards the production of commercial non-transgenic, powdery mildew-resistant bread wheat varieties.

Identification and expression analysis of the NAC transcription factor family in durum wheat (Triticum turgidum L. ssp. durum)

The NAC (NAM, ATAF and CUC) proteins belong to one of the largest plant-specific transcription factor (TF) families and play important roles in plant development processes, response to biotic and abiotic cues and hormone signaling. Our analysis led to the identification of 168 NAC genes in durum wheat, including nine putative membrane-bound TFs and 48 homeologous genes pairs. Phylogenetic analyses of TtNACs along with their Arabidopsis, grape, barley and rice counterparts divided these proteins into 8 phylogenetic groups and allowed the identification of TtNAC-A7, TtNAC-B35, TtNAC-A68, TtNAC-B69 and TtNAC-A43 as homologs of OsNAC1, OsNAC8, OsNTL2, OsNTL5 and ANAC025/NTL14, respectively. In silico expression analysis, using RNA-seq data, revealed tissue-specific and stress responsive TtNAC genes. The expression of ten selected genes was analyzed under salt and drought stresses in two contrasting tolerance cultivars. This analysis is the first report of NAC gene family in durum wheat and will be useful for the identification and selection of candidate genes associated with stress tolerance.

Uncovering hidden variation in polyploid wheat

Comprehensive reverse genetic resources, which have been key to understanding gene function in diploid model organisms, are missing in many polyploid crops. Young polyploid species such as wheat, which was domesticated less than 10,000 y ago, have high levels of sequence identity among subgenomes that mask the effects of recessive alleles. Such redundancy reduces the probability of selection of favorable mutations during natural or human selection, but also allows wheat to tolerate high densities of induced mutations. Here we exploited this property to sequence and catalog more than 10 million mutations in the protein-coding regions of 2,735 mutant lines of tetraploid and hexaploid wheat. We detected, on average, 2,705 and 5,351 mutations per tetraploid and hexaploid line, respectively, which resulted in 35-40 mutations per kb in each population. With these mutation densities, we identified an average of 23-24 missense and truncation alleles per gene, with at least one truncation or deleterious missense mutation in more than 90% of the captured wheat genes per population. This public collection of mutant seed stocks and sequence data enables rapid identification of mutations in the different copies of the wheat genes, which can be combined to uncover previously hidden variation. Polyploidy is a central phenomenon in plant evolution, and many crop species have undergone recent genome duplication events. Therefore, the general strategy and methods developed herein can benefit other polyploid crops.

Genome-wide exploration of metal tolerance protein (MTP) genes in common wheat (Triticum aestivum): insights into metal homeostasis and biofortification

Metal transport process in plants is a determinant of quality and quantity of the harvest. Although it is among the most important of staple crops, knowledge about genes that encode for membrane-bound metal transporters is scarce in wheat. Metal tolerance proteins (MTPs) are involved in trace metal homeostasis at the sub-cellular level, usually by providing metal efflux out of the cytosol. Here, by using various bioinformatics approaches, genes that encode for MTPs in the hexaploid wheat genome (Triticum aestivum, abbreviated as Ta) were identified and characterized. Based on the comparison with known rice MTPs, the wheat genome contained 20 MTP sequences; named as TaMTP1-8A, B and D. All TaMTPs contained a cation diffusion facilitator (CDF) family domain and most members harbored a zinc transporter dimerization domain. Based on motif, phylogeny and alignment analysis, A, B and D genomes of TaMTP3-7 sequences demonstrated higher homology compared to TaMTP1, 2 and 8. With reference to their rice orthologs, TaMTP1s and TaMTP8s belonged to Zn-CDFs, TaMTP2s to Fe/Zn-CDFs and TaMTP3-7s to Mn-CDFs. Upstream regions of TaMTP genes included diverse cis-regulatory motifs, indicating regulation by developmental stage, tissue type and stresses. A scan of the coding sequences of 20 TaMTPs against published miRNAs predicted a total of 14 potential miRNAs, mainly targeting the members of most diverged groups. Expression analysis showed that several TaMTPs were temporally and spatially regulated during the developmental time-course. In grains, MTPs were preferentially expressed in the aleurone layer, which is known as a reservoir for high concentrations of iron and zinc. The work identified and characterized metal tolerance proteins in common wheat and revealed a potential involvement of MTPs in providing a sink for trace element storage in wheat grains.

The in silico identification and characterization of a bread wheat/Triticum militinae introgression line

The capacity of the bread wheat (Triticum aestivum) genome to tolerate introgression from related genomes can be exploited for wheat improvement. A resistance to powdery mildew expressed by a derivative of the cross-bread wheat cv. Tähti × T. militinae (Tm) is known to be due to the incorporation of a Tm segment into the long arm of chromosome 4A. Here, a newly developed in silico method termed rearrangement identification and characterization (RICh) has been applied to characterize the introgression. A virtual gene order, assembled using the GenomeZipper approach, was obtained for the native copy of chromosome 4A; it incorporated 570 4A DArTseq markers to produce a zipper comprising 2132 loci. A comparison between the native and introgressed forms of the 4AL chromosome arm showed that the introgressed region is located at the distal part of the arm. The Tm segment, derived from chromosome 7G, harbours 131 homoeologs of the 357 genes present on the corresponding region of Chinese Spring 4AL. The estimated number of Tm genes transferred along with the disease resistance gene was 169. Characterizing the introgression's position, gene content and internal gene order should not only facilitate gene isolation, but may also be informative with respect to chromatin structure and behaviour studies.

The iSelect 9 K SNP analysis revealed polyploidization induced revolutionary changes and intense human selection causing strong haplotype blocks in wheat

A Chinese wheat mini core collection was genotyped using the wheat 9 K iSelect SNP array. Total 2420 and 2396 polymorphic SNPs were detected on the A and the B genome chromosomes, which formed 878 haplotype blocks. There were more blocks in the B genome, but the average block size was significantly (P < 0.05) smaller than those in the A genome. Intense selection (domestication and breeding) had a stronger effect on the A than on the B genome chromosomes. Based on the genetic pedigrees, many blocks can be traced back to a well-known Strampelli cross, which was made one century ago. Furthermore, polyploidization of wheat (both tetraploidization and hexaploidization) induced revolutionary changes in both the A and the B genomes, with a greater increase of gene diversity compared to their diploid ancestors. Modern breeding has dramatically increased diversity in the gene coding regions, though obvious blocks were formed on most of the chromosomes in both tetraploid and hexaploid wheats. Tag-SNP markers identified in this study can be used for marker assisted selection using haplotype blocks as a wheat breeding strategy. This strategy can also be employed to facilitate genome selection in other self-pollinating crop species.

A High Resolution Radiation Hybrid Map of Wheat Chromosome 4A

Bread wheat has a large and complex allohexaploid genome with low recombination level at chromosome centromeric and peri-centromeric regions. This significantly hampers ordering of markers, contigs of physical maps and sequence scaffolds and impedes obtaining of high-quality reference genome sequence. Here we report on the construction of high-density and high-resolution radiation hybrid (RH) map of chromosome 4A supported by high-density chromosome deletion map. A total of 119 endosperm-based RH lines of two RH panels and 15 chromosome deletion bin lines were genotyped with 90K iSelect single nucleotide polymorphism (SNP) array. A total of 2316 and 2695 markers were successfully mapped to the 4A RH and deletion maps, respectively. The chromosome deletion map was ordered in 19 bins and allowed precise identification of centromeric region and verification of the RH panel reliability. The 4A-specific RH map comprises 1080 mapping bins and spans 6550.9 cR with a resolution of 0.13 Mb/cR. Significantly higher mapping resolution in the centromeric region was observed as compared to recombination maps. Relatively even distribution of deletion frequency along the chromosome in the RH panel was observed and putative functional centromere was delimited within a region characterized by two SNP markers.

Overexpression of TaWRKY146 Increases Drought Tolerance through Inducing Stomatal Closure in Arabidopsis thaliana

As a superfamily of transcription factors, the tryptophan-arginine-lysine-tyrosine (WRKY) transcription factors have been found to be essential for abiotic and biotic stress responses in plants. Currently, only 76 WRKY transcription factors in wheat could be identified in the NCBI database, among which only a few have been functionally analyzed. Herein, a total of 188 WRKY transcription factors were identified from the wheat genome database, which included 123 full-length coding sequences, and all of them were used for detailed evolution studies. By bioinformatics analysis, a WRKY transcription factor, named TaWRKY146, was found to be the homologous gene of AtWRKY46, overexpression of which leads to hypersensitivity to drought and salt stress in Arabidopsis. Consequently, the full length of TaWRKY146 was cloned, and the expression levels of TaWRKY146 were found significantly up-regulated in the leaves and roots of wheat seedlings, which were subjected to osmotic stress. Overexpression of TaWRKY146 in Arabidopsis was shown to enhance drought tolerance by the induction of stomatal closure that reduced the transpiration rate. All these results provide a firm foundation for further identification of WRKY transcription factors with important functions in wheat.

Enabling Molecular Technologies for Trait Improvement in Wheat

Wheat is the major staple food crop and a source of calories for humans worldwide. A steady increase in the wheat production is essential to meet the demands of an ever-increasing global population and to achieve food security. The large size and structurally intricate genome of polyploid wheat had hindered the genomic analysis. However, with the advent of new genomic technologies such as next generation sequencing has led to genome drafts for bread wheat and its progenitors and has paved the way to design new strategies for crop improvement. Here we provide an overview of the advancements made in wheat genomics together with the available "omics approaches" and bioinformatics resources developed for wheat research. Advances in genomic, transcriptomic, and metabolomic technologies are highlighted as options to circumvent existing bottlenecks in the phenotypic and genomic selection and gene transfer. The contemporary reverse genetics approaches, including the novel genome editing techniques to inform targeted manipulation of a single/multiple genes and strategies for generating marker-free transgenic wheat plants, emphasize potential to revolutionize wheat improvement shortly.

A Wheat R2R3-type MYB Transcription Factor TaODORANT1 Positively Regulates Drought and Salt Stress Responses in Transgenic Tobacco Plants

MYB transcription factors play important roles in plant responses to biotic and abiotic stress. In this study, TaODORANT1, a R2R3-MYB gene, was cloned from wheat (Triticum aestivum L.). TaODORANT1 was localized in the nucleus and functioned as a transcriptional activator. TaODORANT1 was up-regulated in wheat under PEG6000, NaCl, ABA, and H₂O₂ treatments. TaODORANT1-overexpressing transgenic tobacco plants exhibited higher relative water content and lower water loss rate under drought stress, as well as lower Na⁺ accumulation in leaves under salt stress. The transgenic plants showed higher CAT activity but lower ion leakage, H₂O₂ and malondialdehyde contents under drought and salt stresses. Besides, the transgenic plants also exhibited higher SOD activity under drought stress. Our results also revealed that TaODORANT1 overexpression up-regulated the expression of several ROS- and stress-related genes in response to both drought and salt stresses, thus enhancing transgenic tobacco plants tolerance. Our studies demonstrate that TaODORANT1 positively regulates plant tolerance to drought and salt stresses.

Molecular Characterization, Gene Evolution, and Expression Analysis of the Fructose-1, 6-bisphosphate Aldolase (FBA) Gene Family in Wheat (Triticum aestivum L.)

Fructose-1, 6-bisphosphate aldolase (FBA) is a key plant enzyme that is involved in glycolysis, gluconeogenesis, and the Calvin cycle. It plays significant roles in biotic and abiotic stress responses, as well as in regulating growth and development processes. In the present paper, 21 genes encoding TaFBA isoenzymes were identified, characterized, and categorized into three groups: class I chloroplast/plastid FBA (CpFBA), class I cytosol FBA (cFBA), and class II chloroplast/plastid FBA. By using a prediction online database and genomic PCR analysis of Chinese Spring nulli-tetrasomic lines, we have confirmed the chromosomal location of these genes in 12 chromosomes of four homologous groups. Sequence and genomic structure analysis revealed the high identity of the allelic TaFBA genes and the origin of different TaFBA genes. Numerous putative environment stimulus-responsive cis-elements have been identified in 1,500-bp regions of TaFBA gene promoters, of which the most abundant are the light-regulated elements (LREs). Phylogenetic reconstruction using the deduced protein sequence of 245 FBA genes indicated an independent evolutionary pathway for the class I and class II groups. Although, earlier studies have indicated that class II FBA only occurs in prokaryote and fungi, our results have demonstrated that a few class II CpFBAs exist in wheat and other closely related species. Class I TaFBA was predicted to be tetramers and class II to be dimers. Gene expression analysis based on microarray and transcriptome databases suggested the distinct role of TaFBAs in different tissues and developmental stages. The TaFBA 4-9 genes were highly expressed in leaves and might play important roles in wheat development. The differential expression patterns of the TaFBA genes in light/dark and a few abiotic stress conditions were also analyzed. The results suggested that LRE cis-elements of TaFBA gene promoters were not directly related to light responses. Most TaFBA genes had higher expression levels in the roots than in the shoots when under various stresses. Class I cytosol TaFBA genes, particularly TaFBA10/12/18 and TaFBA13/16, and three class II TaFBA genes are involved in responses to various abiotic stresses. Class I CpFBA genes in wheat are apparently sensitive to different stress conditions.

A High-Density Consensus Map of Common Wheat Integrating Four Mapping Populations Scanned by the 90K SNP Array

A high-density consensus map is a powerful tool for gene mapping, cloning and molecular marker-assisted selection in wheat breeding. The objective of this study was to construct a high-density, single nucleotide polymorphism (SNP)-based consensus map of common wheat (Triticum aestivum L.) by integrating genetic maps from four recombinant inbred line populations. The populations were each genotyped using the wheat 90K Infinium iSelect SNP assay. A total of 29,692 SNP markers were mapped on 21 linkage groups corresponding to 21 hexaploid wheat chromosomes, covering 2,906.86 cM, with an overall marker density of 10.21 markers/cM. Compared with the previous maps based on the wheat 90K SNP chip detected 22,736 (76.6%) of the SNPs with consistent chromosomal locations, whereas 1,974 (6.7%) showed different chromosomal locations, and 4,982 (16.8%) were newly mapped. Alignment of the present consensus map and the wheat expressed sequence tags (ESTs) Chromosome Bin Map enabled assignment of 1,221 SNP markers to specific chromosome bins and 819 ESTs were integrated into the consensus map. The marker orders of the consensus map were validated based on physical positions on the wheat genome with Spearman rank correlation coefficients ranging from 0.69 (4D) to 0.97 (1A, 4B, 5B, and 6A), and were also confirmed by comparison with genetic position on the previously 40K SNP consensus map with Spearman rank correlation coefficients ranging from 0.84 (6D) to 0.99 (6A). Chromosomal rearrangements reported previously were confirmed in the present consensus map and new putative rearrangements were identified. In addition, an integrated consensus map was developed through the combination of five published maps with ours, containing 52,607 molecular markers. The consensus map described here provided a high-density SNP marker map and a reliable order of SNPs, representing a step forward in mapping and validation of chromosomal locations of SNPs on the wheat 90K array. Moreover, it can be used as a reference for quantitative trait loci (QTL) mapping to facilitate exploitation of genes and QTL in wheat breeding.

Genome-metabolite associations revealed low heritability, high genetic complexity, and causal relations for leaf metabolites in winter wheat (Triticum aestivum)

We investigated associations between the metabolic phenotype, consisting of quantitative data of 76 metabolites from 135 contrasting winter wheat (Triticum aestivum) lines, and 17 372 single nucleotide polymorphism (SNP) markers. Metabolite profiles were generated from flag leaves of plants from three different environments, with average repeatabilities of 0.5-0.6. The average heritability of 0.25 was unaffected by the heading date. Correlations among metabolites reflected their functional grouping, highlighting the strict coordination of various routes of the citric acid cycle. Genome-wide association studies identified significant associations for six metabolic traits, namely oxalic acid, ornithine, L-arginine, pentose alcohol III, L-tyrosine, and a sugar oligomer (oligo II), with between one and 17 associated SNPs. Notable associations with genes regulating transcription or translation explained between 2.8% and 32.5% of the genotypic variance (pG). Further candidate genes comprised metabolite carriers (pG 32.5-38.1%), regulatory proteins (pG 0.3-11.1%), and metabolic enzymes (pG 2.5-32.5%). The combinatorial use of genomic and metabolic data to construct partially directed networks revealed causal inferences in the correlated metabolite traits and associated SNPs. The evaluated causal relationships will provide a basis for predicting the effects of genetic interferences on groups of correlated metabolic traits, and thus on specific metabolic phenotypes.

Overview of the Wheat Genetic Transformation and Breeding Status in China

In the past two decades, Chinese scientists have achieved significant progress on three aspects of wheat genetic transformation. First, the wheat transformation platform has been established and optimized to improve the transformation efficiency, shorten the time required from starting of transformation procedure to the fertile transgenic wheat plants obtained as well as to overcome the problem of genotype-dependent for wheat genetic transformation in wide range of wheat elite varieties. Second, with the help of many emerging techniques such as CRISPR/cas9 function of over 100 wheat genes has been investigated. Finally, modern technology has been combined with the traditional breeding technique such as crossing to accelerate the application of wheat transformation. Overall, the wheat end-use quality and the characteristics of wheat stress tolerance have been improved by wheat genetic engineering technique. So far, wheat transgenic lines integrated with quality-improved genes and stress tolerant genes have been on the way of Production Test stage in the field. The debates and the future studies on wheat transformation have been discussed, and the brief summary of Chinese wheat breeding research history has also been provided in this review.

Genome-wide nucleotide diversity of hatchery-reared Atlantic and Mediterranean strains of brown trout Salmo trutta compared to wild Mediterranean populations

A genome-wide assessment of diversity is provided for wild Mediterranean brown trout Salmo trutta populations from headwater tributaries of the Orb River and from Atlantic and Mediterranean hatchery-reared strains that have been used for stocking. Double-digest restriction-site-associated DNA sequencing (dd-RADseq) was performed and the efficiency of de novo and reference-mapping approaches to obtain individual genotypes was compared. Large numbers of single nucleotide polymorphism (SNP) markers with similar genome-wide distributions were discovered using both approaches (196 639 v. 121 016 SNPs, respectively), with c. 80% of the loci detected de novo being also found with reference mapping, using the Atlantic salmon Salmo salar genome as a reference. Lower mapping density but larger nucleotide diversity (π) was generally observed near extremities of linkage groups, consistent with regions of residual tetrasomic inheritance observed in salmonids. Genome-wide diversity estimates revealed reduced polymorphism in hatchery strains (π = 0·0040 and π = 0·0029 in Atlantic and Mediterranean strains, respectively) compared to wild populations (π = 0·0049), a pattern that was congruent with allelic richness estimated from microsatellite markers. Finally, pronounced heterozygote deficiency was found in hatchery strains (Atlantic F_IS = 0·18; Mediterranean F_IS = 0·42), indicating that stocking practices may affect the genetic diversity in wild populations. These new genomic resources will provide important tools to define better conservation strategies in S. trutta.

Pseudogenes and Their Genome-Wide Prediction in Plants

Pseudogenes are paralogs generated from ancestral functional genes (parents) during genome evolution, which contain critical defects in their sequences, such as lacking a promoter, having a premature stop codon or frameshift mutations. Generally, pseudogenes are functionless, but recent evidence demonstrates that some of them have potential roles in regulation. The majority of pseudogenes are generated from functional progenitor genes either by gene duplication (duplicated pseudogenes) or retro-transposition (processed pseudogenes). Pseudogenes are primarily identified by comparison to their parent genes. Bioinformatics tools for pseudogene prediction have been developed, among which PseudoPipe, PSF and Shiu's pipeline are publicly available. We compared these three tools using the well-annotated Arabidopsis thaliana genome and its known 924 pseudogenes as a test data set. PseudoPipe and Shiu's pipeline identified ~80% of A. thaliana pseudogenes, of which 94% were shared, while PSF failed to generate adequate results. A need for improvement of the bioinformatics tools for pseudogene prediction accuracy in plant genomes was thus identified, with the ultimate goal of improving the quality of genome annotation in plants.

Ca2+/Cation Antiporters (CaCA): Identification, Characterization and Expression Profiling in Bread Wheat (Triticum aestivum L.)

The Ca2+/cation antiporters (CaCA) superfamily proteins play vital function in Ca2+ ion homeostasis, which is an important event during development and defense response. Molecular characterization of these proteins has been performed in certain plants, but they are still not characterized in Triticum aestivum (bread wheat). Herein, we identified 34 TaCaCA superfamily proteins, which were classified into TaCAX, TaCCX, TaNCL, and TaMHX protein families based on their structural organization and evolutionary relation with earlier reported proteins. Since the T. aestivum comprises an allohexaploid genome, TaCaCA genes were derived from each A, B, and D subgenome and homeologous chromosome (HC), except chromosome-group 1. Majority of genes were derived from more than one HCs in each family that were considered as homeologous genes (HGs) due to their high similarity with each other. These HGs showed comparable gene and protein structures in terms of exon/intron organization and domain architecture. Majority of TaCaCA proteins comprised two Na_Ca_ex domains. However, TaNCLs consisted of an additional EF-hand domain with calcium binding motifs. Each TaCaCA protein family consisted of about 10 transmembrane and two α-repeat regions with specifically conserved signature motifs except TaNCL, which had single α-repeat. Variable expression of most of the TaCaCA genes during various developmental stages suggested their specified role in development. However, constitutively high expression of a few genes like TaCAX1-A and TaNCL1-B indicated their role throughout the plant growth and development. The modulated expression of certain genes during biotic (fungal infections) and abiotic stresses (heat, drought, salt) suggested their role in stress response. Majority of TaCCX and TaNCL family genes were found highly affected during various abiotic stresses. However, the role of individual gene needs to be established. The present study unfolded the opportunity for detail functional characterization of TaCaCA proteins and their utilization in future crop improvement programs.

Bioinformatic Analyses of Subgroup-A Members of the Wheat bZIP Transcription Factor Family and Functional Identification of TabZIP174 Involved in Drought Stress Response

Extensive studies in Arabidopsis and rice have demonstrated that Subgroup-A members of the bZIP transcription factor family play important roles in plant responses to multiple abiotic stresses. Although common wheat (Triticum aestivum) is one of the most widely cultivated and consumed food crops in the world, there are limited investigations into Subgroup A of the bZIP family in wheat. In this study, we performed bioinformatic analyses of the 41 Subgroup-A members of the wheat bZIP family. Phylogenetic and conserved motif analyses showed that most of the Subgroup-A bZIP proteins involved in abiotic stress responses of wheat, Arabidopsis, and rice clustered in Clade A1 of the phylogenetic tree, and shared a majority of conserved motifs, suggesting the potential importance of Clade-A1 members in abiotic stress responses. Gene structure analysis showed that TabZIP genes with close phylogenetic relationships tended to possess similar exon-intron compositions, and the positions of introns in the hinge regions of the bZIP domains were highly conserved, whereas introns in the leucine zipper regions were at variable positions. Additionally, eleven groups of homologs and two groups of tandem paralogs were also identified in Subgroup A of the wheat bZIP family. Expression profiling analysis indicated that most Subgroup-A TabZIP genes were responsive to abscisic acid and various abiotic stress treatments. TabZIP27, TabZIP74, TabZIP138, and TabZIP174 proteins were localized in the nucleus of wheat protoplasts, whereas TabZIP9-GFP fusion protein was simultaneously present in the nucleus, cytoplasm, and cell membrane. Transgenic Arabidopsis overexpressing TabZIP174 displayed increased seed germination rates and primary root lengths under drought treatments. Overexpression of TabZIP174 in transgenic Arabidopsis conferred enhanced drought tolerance, and transgenic plants exhibited lower water loss rates, higher survival rates, higher proline, soluble sugar, and leaf chlorophyll contents, as well as more stable osmotic potential under drought conditions. Additionally, overexpression of TabZIP174 increased the expression of stress-responsive genes (RD29A, RD29B, RAB18, DREB2A, COR15A, and COR47). The improved drought resistance might be attributed to the increased osmotic adjustment capacity. Our results indicate that TabZIP174 may participate in regulating plant response to drought stress and holds great potential for genetic improvement of abiotic stress tolerance in crops.

Chromosome-specific sequencing reveals an extensive dispensable genome component in wheat

The hexaploid wheat genotype Chinese Spring (CS) has been used worldwide as the reference base for wheat genetics and genomics, and significant resources have been used by the international community to generate a reference wheat genome based on this genotype. By sequencing flow-sorted 3B chromosome from a hexaploid wheat genotype CRNIL1A and comparing the obtained sequences with those available for CS, we detected that a large number of sequences in the former were missing in the latter. If the distribution of such sequences in the hexaploid wheat genome is random, CRNILA sequences missing in CS could be as much as 159.3 Mb even if only fragments of 50 bp or longer were considered. Analysing RNA sequences available in the public domains also revealed that dispensable genes are common in hexaploid wheat. Together with those extensive intra- and interchromosomal rearrangements in CS, the existence of such dispensable genes is another factor highlighting potential issues with the use of reference genomes in various studies. Strong deviation in distributions of these dispensable sequences among genotypes with different geographical origins provided the first evidence indicating that they could be associated with adaptation in wheat.

Rapid gene isolation in barley and wheat by mutant chromosome sequencing

Identification of causal mutations in barley and wheat is hampered by their large genomes and suppressed recombination. To overcome these obstacles, we have developed MutChromSeq, a complexity reduction approach based on flow sorting and sequencing of mutant chromosomes, to identify induced mutations by comparison to parental chromosomes. We apply MutChromSeq to six mutants each of the barley Eceriferum-q gene and the wheat Pm2 genes. This approach unambiguously identified single candidate genes that were verified by Sanger sequencing of additional mutants. MutChromSeq enables reference-free forward genetics in barley and wheat, thus opening up their pan-genomes to functional genomics.