Small-scale gene duplications played a major role in the recent evolution of wheat chromosome 3B

A FLOWERING LOCUS T ortholog is associated with photoperiod-insensitive flowering in hemp (Cannabis sativa L.)

Hemp (Cannabis sativa L.) is an extraordinarily versatile crop, with applications ranging from medicinal compounds to seed oil and fibre products. Cannabis sativa is a short-day plant, and its flowering is highly controlled by photoperiod. However, substantial genetic variation exists for photoperiod sensitivity in C. sativa, and photoperiod-insensitive ("autoflower") cultivars are available. Using a bi-parental mapping population and bulked segregant analysis, we identified Autoflower2, a 0.5 Mbp locus significantly associated with photoperiod-insensitive flowering in hemp. Autoflower2 contains an ortholog of the central flowering time regulator FLOWERING LOCUS T (FT) from Arabidopsis thaliana which we termed CsFT1. We identified extensive sequence divergence between alleles of CsFT1 from photoperiod-sensitive and insensitive cultivars of C. sativa, including a duplication of CsFT1 and sequence differences, especially in introns. Furthermore, we observed higher expression of one of the CsFT1 copies found in the photoperiod-insensitive cultivar. Genotyping of several mapping populations and a diversity panel confirmed a correlation between CsFT1 alleles and photoperiod response, affirming that at least two independent loci involved in the photoperiodic control of flowering, Autoflower1 and Autoflower2, exist in the C. sativa gene pool. This study reveals the multiple independent origins of photoperiod insensitivity in C. sativa, supporting the likelihood of a complex domestication history in this species. By integrating the genetic relaxation of photoperiod sensitivity into novel C. sativa cultivars, expansion to higher latitudes will be permitted, thus allowing the full potential of this versatile crop to be reached.

Genome-Wide Analysis of MYB Gene Family in Chrysanthemum ×morifolium Provides Insights into Flower Color Regulation

MYBs constitute the second largest transcription factor (TF) superfamily in flowering plants with substantial structural and functional diversity, which have been brought into focus because they affect flower colors by regulating anthocyanin biosynthesis. Up to now, the genomic data of several Chrysanthemum species have been released, which provides us with abundant genomic resources for revealing the evolution of the MYB gene family in Chrysanthemum species. In the present study, comparative analyses of the MYB gene family in six representative species, including C. lavandulifolium, C. seticuspe, C. ×morifolium, Helianthus annuus, Lactuca sativa, and Arabidopsis thaliana, were performed. A total of 1104 MYBs, which were classified into four subfamilies and 35 lineages, were identified in the three Chrysanthemum species (C. lavandulifolium, C. seticuspe, and C. ×morifolium). We found that whole-genome duplication and tandem duplication are the main duplication mechanisms that drove the occurrence of duplicates in CmMYBs (particularly in the R2R3-MYB subfamily) during the evolution of the cultivated chrysanthemums. Sequence structure and selective pressure analyses of the MYB gene family revealed that some of R2R3-MYBs were subjected to positive selection, which are mostly located on the distal telomere segments of the chromosomes and contain motifs 7 and 8. In addition, the gene expression analysis of CmMYBs in different organs and at various capitulum developmental stages of C. ×morifolium indicated that CmMYBS2, CmMYB96, and CmMYB109 might be the negative regulators for anthocyanin biosynthesis. Our results provide the phylogenetic context for research on the genetic and functional evolution of the MYB gene family in Chrysanthemum species and deepen our understanding of the regulatory mechanism of MYB TFs on the flower color of C. ×morifolium.

Genomic insights into the chromosomal elongation in a family of Collembola

Collembola is a highly diverse and abundant group of soil arthropods with chromosome numbers ranging from 5 to 11. Previous karyotype studies indicated that the Tomoceridae family possesses an exceptionally long chromosome. To better understand chromosome size evolution in Collembola, we obtained a chromosome-level genome of Yoshiicerus persimilis with a size of 334.44 Mb and BUSCO completeness of 97.0% (n = 1013). Both genomes of Y. persimilis and Tomocerus qinae (recently published) have an exceptionally large chromosome (ElChr greater than 100 Mb), accounting for nearly one-third of the genome. Comparative genomic analyses suggest that chromosomal elongation occurred independently in the two species approximately 10 million years ago, rather than in the ancestor of the Tomoceridae family. The ElChr elongation was caused by large tandem and segmental duplications, as well as transposon proliferation, with genes in these regions experiencing weaker purifying selection (higher dN/dS) than conserved regions. Moreover, inter-genomic synteny analyses indicated that chromosomal fission/fusion events played a crucial role in the evolution of chromosome numbers (ranging from 5 to 7) within Entomobryomorpha. This study provides a valuable resource for investigating the chromosome evolution of Collembola.

Contribution of Duplicated Nucleotide-Binding Leucine-Rich Repeat (NLR) Genes to Wheat Disease Resistance

Wheat has a large and diverse repertoire of NLRs involved in disease resistance, with over 1500 NLRs detected in some studies. These NLR genes occur as singletons or clusters containing copies of NLRs from different phylogenetic clades. The number of NLRs and cluster size can differ drastically among ecotypes and cultivars. Primarily, duplication has led to the evolution and diversification of NLR genes. Among the various mechanisms, whole genome duplication (WGD) is the most intense and leading cause, contributing to the complex evolutionary history and abundant gene set of hexaploid wheat. Tandem duplication or recombination is another major mechanism of NLR gene expansion in wheat. The diversity and divergence of duplicate NLR genes are responsible for the broad-spectrum resistance of most plant species with limited R genes. Understanding the mechanisms underlying the rapid evolution and diversification of wheat NLR genes will help improve disease resistance in crops. The present review focuses on the diversity and divergence of duplicate NLR genes and their contribution to wheat disease resistance. Moreover, we provide an overview of disease resistance-associated gene duplication and the underlying strategies in wheat.

The evolutionary patterns, expression profiles, and genetic diversity of expanded genes in barley

Gene duplication resulting from whole-genome duplication (WGD), small-scale duplication (SSD), or unequal hybridization plays an important role in the expansion of gene families. Gene family expansion can also mediate species formation and adaptive evolution. Barley (Hordeum vulgare) is the world's fourth largest cereal crop, and it contains valuable genetic resources due to its ability to tolerate various types of environmental stress. In this study, 27,438 orthogroups in the genomes of seven Poaceae were identified, and 214 of them were significantly expanded in barley. The evolutionary rates, gene properties, expression profiles, and nucleotide diversity between expanded and non-expanded genes were compared. Expanded genes evolved more rapidly and experienced lower negative selection. Expanded genes, including their exons and introns, were shorter, they had fewer exons, their GC content was lower, and their first exons were longer compared with non-expanded genes. Codon usage bias was also lower for expanded genes than for non-expanded genes; the expression levels of expanded genes were lower than those of non-expanded genes, and the expression of expanded genes showed higher tissue specificity than that of non-expanded genes. Several stress-response-related genes/gene families were identified, and these genes could be used to breed barley plants with greater resistance to environmental stress. Overall, our analysis revealed evolutionary, structural, and functional differences between expanded and non-expanded genes in barley. Additional studies are needed to clarify the functions of the candidate genes identified in our study and evaluate their utility for breeding barley plants with greater stress resistance.

Location and Identification on Chromosome 3B of Bread Wheat of Genes Affecting Chiasma Number

Understanding meiotic crossover (CO) variation in crops like bread wheat (Triticum aestivum L.) is necessary as COs are essential to create new, original and powerful combinations of genes for traits of agronomical interest. We cytogenetically characterized a set of wheat aneuploid lines missing part or all of chromosome 3B to identify the most influential regions for chiasma formation located on this chromosome. We showed that deletion of the short arm did not change the total number of chiasmata genome-wide, whereas this latter was reduced by ~35% while deleting the long arm. Contrary to what was hypothesized in a previous study, deletion of the long arm does not disturb the initiation of the synaptonemal complex (SC) in early meiotic stages. However, progression of the SC is abnormal, and we never observed its completion when the long arm is deleted. By studying six different deletion lines (missing different parts of the long arm), we revealed that at least two genes located in both the proximal (C-3BL2-0.22) and distal (3BL7-0.63-1.00) deletion bins are involved in the control of chiasmata, each deletion reducing the number of chiasmata by ~15%. We combined sequence analyses of deletion bins with RNA-Seq data derived from meiotic tissues and identified a set of genes for which at least the homoeologous copy on chromosome 3B is expressed and which are involved in DNA processing. Among these genes, eight (CAP-E1/E2, DUO1, MLH1, MPK4, MUS81, RTEL1, SYN4, ZIP4) are known to be involved in the recombination pathway.

Identification of ankyrin-transmembrane-type subfamily genes in Triticeae species reveals TaANKTM2A-5 regulates powdery mildew resistance in wheat

The ankyrin-transmembrane (ANKTM) subfamily is the most abundant subgroup of the ANK superfamily, with critical roles in pathogen defense. However, the function of ANKTM proteins in wheat immunity remains largely unexplored. Here, a total of 381 ANKTMs were identified from five Triticeae species and Arabidopsis, constituting five classes. Among them, class a only contains proteins from Triticeae species and the number of ANKTM in class a of wheat is significantly larger than expected, even after consideration of the ploidy level. Tandem duplication analysis of ANKTM indicates that Triticum urartu, Triticum dicoccoides and wheat all had experienced tandem duplication events which in wheat-produced ANKTM genes all clustered in class a. The above suggests that not only did the genome polyploidization result in the increase of ANKTM gene number, but that tandem duplication is also a mechanism for the expansion of this subfamily. Micro-collinearity analysis of Triticeae ANKTMs indicates that some ANKTM type genes evolved into other types of ANKs in the evolution process. Public RNA-seq data showed that most of the genes in class d and class e are expressed, and some of them show differential responses to biotic stresses. Furthermore, qRT-PCR results showed that some ANKTMs in class d and class e responded to powdery mildew. Silencing of TaANKTM2A-5 by barley stripe mosaic virus-induced gene silencing compromised powdery mildew resistance in common wheat Bainongaikang58. Findings in this study not only help to understand the evolutionary process of ANKTM genes, but also form the basis for exploring disease resistance genes in the ANKTM gene family.

Genome-Wide Diversity of MADS-Box Genes in Bread Wheat is Associated with its Rapid Global Adaptability

MADS-box gene family members play multifarious roles in regulating the growth and development of crop plants and hold enormous promise for bolstering grain yield potential under changing global environments. Bread wheat (Triticum aestivum L.) is a key stable food crop around the globe. Until now, the available information concerning MADS-box genes in the wheat genome has been insufficient. Here, a comprehensive genome-wide analysis identified 300 high confidence MADS-box genes from the publicly available reference genome of wheat. Comparative phylogenetic analyses with Arabidopsis and rice MADS-box genes classified the wheat genes into 16 distinct subfamilies. Gene duplications were mainly identified in subfamilies containing unbalanced homeologs, pointing towards a potential mechanism for gene family expansion. Moreover, a more rapid evolution was inferred for M-type genes, as compared with MIKC-type genes, indicating their significance in understanding the evolutionary history of the wheat genome. We speculate that subfamily-specific distal telomeric duplications in unbalanced homeologs facilitate the rapid adaptation of wheat to changing environments. Furthermore, our in-silico expression data strongly proposed MADS-box genes as active guardians of plants against pathogen insurgency and harsh environmental conditions. In conclusion, we provide an entire complement of MADS-box genes identified in the wheat genome that could accelerate functional genomics efforts and possibly facilitate bridging gaps between genotype-to-phenotype relationships through fine-tuning of agronomically important traits.

Genome-wide identification and phylogenetic relationships of the Hsp70 gene family of Aegilops tauschii, wild emmer wheat (Triticum dicoccoides) and bread wheat (Triticum aestivum)

Heat shock protein 70 (Hsp70) plays an important role in plant development. It is closely related to the physiological process of cell development and the response to abiotic and biological stress. However, the classification and evolution of Hsp70 genes in bread wheat, wild emmer wheat and Aegilops tauschii are still unclear. Therefore, this study conducted a comprehensive bioinformatics analysis of Hsp70 gene in three species. Among these three species, 113, 79 and 36 Hsp70 genes were identified. They are divided into six subfamilies. Group vi-1 is different from Arabidopsis thaliana. It may be the result of early evolutionary segregation. The number of exons in different subfamilies (from 1 to 13) was different, but the distribution patterns of exons / introns in the same subfamily were similar. The results of Hsp70 promoter region analysis showed that the cis-regulatory elements of A. tauschii and wild emmer wheat were different from those of wheat. In addition, CpG island proportion of wild emmer Hsp70 was higher than that of wheat, which may be the molecular basis of heat resistance of wild wheat relative to cultivated wheat. Further comprehensive analysis of chromosome location and repeat events of Hsp70 gene showed that whole-genome duplication and tandem duplication events contributed to the evolution and expansion of Hsp70 gene in wheat. The results of non-synonymous substitution and synonymous substitution analysis showed that Hsp70 genes of three species had undergone purification selection. The expression profile analysis showed that Hsp70 gene was highly expressed in the roots during the vegetative growth period. In addition, TaHsp70 gene was highly expressed under various stress. The identification, classification and evolution of Hsp70 in wheat and its relatives provided a basis for further research on its evolution and its molecular mechanism in response to stress.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02639-5.

Homoeolog Inference Methods Requiring Bidirectional Best Hits or Synteny Miss Many Pairs

Homoeologs are pairs of genes or chromosomes in the same species that originated by speciation and were brought back together in the same genome by allopolyploidization. Bioinformatic methods for accurate homoeology inference are crucial for studying the evolutionary consequences of polyploidization, and homoeology is typically inferred on the basis of bidirectional best hit (BBH) and/or positional conservation (synteny). However, these methods neglect the fact that genes can duplicate and move, both prior to and after the allopolyploidization event. These duplications and movements can result in many-to-many and/or nonsyntenic homoeologs-which thus remain undetected and unstudied. Here, using the allotetraploid upland cotton (Gossypium hirsutum) as a case study, we show that conventional approaches indeed miss a substantial proportion of homoeologs. Additionally, we found that many of the missed pairs of homoeologs are broadly and highly expressed. A gene ontology analysis revealed a high proportion of the nonsyntenic and non-BBH homoeologs to be involved in protein translation and are likely to contribute to the functional repertoire of cotton. Thus, from an evolutionary and functional genomics standpoint, choosing a homoeolog inference method which does not solely rely on 1:1 relationship cardinality or synteny is crucial for not missing these potentially important homoeolog pairs.

Extensive structural variation in the Bowman-Birk inhibitor family in common wheat (Triticum aestivum L.)

BACKGROUND

Bowman-Birk inhibitors (BBI) are a family of serine-type protease inhibitors that modulate endogenous plant proteolytic activities during different phases of development. They also inhibit exogenous proteases as a component of plant defense mechanisms, and their overexpression can confer resistance to phytophagous herbivores and multiple fungal and bacterial pathogens. Dicot BBIs are multifunctional, with a "double-headed" structure containing two separate inhibitory loops that can bind and inhibit trypsin and chymotrypsin proteases simultaneously. By contrast, monocot BBIs have a non-functional chymotrypsin inhibitory loop, although they have undergone internal duplication events giving rise to proteins with multiple BBI domains.

RESULTS

We used a Hidden Markov Model (HMM) profile-based search to identify 57 BBI genes in the common wheat (Triticum aestivum L.) genome. The BBI genes are unevenly distributed, with large gene clusters in the telomeric regions of homoeologous group 1 and 3 chromosomes that likely arose through a series of tandem gene duplication events. The genomes of wheat progenitors also contain contiguous clusters of BBI genes, suggesting this family underwent expansion before the domestication of common wheat. However, the BBI gene family varied in size among different cultivars, showing this family remains dynamic. Because of these expansions, the BBI gene family is larger in wheat than other monocots such as maize, rice and Brachypodium. We found BBI proteins in common wheat with intragenic homologous duplications of cysteine-rich functional domains, including one protein with four functional BBI domains. This diversification may expand the spectrum of target substrates. Expression profiling suggests that some wheat BBI proteins may be involved in regulating endogenous proteases during grain development, while others were induced in response to biotic and abiotic stresses, suggesting a role in plant defense.

CONCLUSIONS

Genome-wide characterization reveals that the BBI gene family in wheat is subject to a high rate of homologous tandem duplication and deletion events, giving rise to a diverse set of encoded proteins. This information will facilitate the functional characterization of individual wheat BBI genes to determine their role in wheat development and stress responses, and their potential application in breeding.

The Landscape of the Genomic Distribution and the Expression of the F-Box Genes Unveil Genome Plasticity in Hexaploid Wheat during Grain Development and in Response to Heat and Drought Stress

FBX proteins are subunits of the SCF complex (Skp1-cullin-FBX) belonging to the E3 ligase family, which is involved in the ubiquitin-proteasome 26S (UPS) pathway responsible for the post-translational protein turnover. By targeting, in a selective manner, key regulatory proteins for ubiquitination and 26S proteasome degradation, FBX proteins play a major role in plant responses to diverse developmental and stress conditions. Although studies on the genomic organization of the FBX gene family in various species have been reported, knowledge related to bread wheat (Triticum aestivum) is scarce and needs to be broadened. Using the latest assembly of the wheat genome, we identified 3670 TaFBX genes distributed non-homogeneously within the three subgenomes (A, B and D) and between the 21 chromosomes, establishing it as one of the richest gene families among plant species. Based on the presence of the five different chromosomal regions previously identified, the present study focused on the genomic distribution of the TaFBX family and the identification of differentially expressed genes during the embryogenesis stages and in response to heat and drought stress. Most of the time, when comparing the expected number of genes (taking into account the formal gene distribution on the entire wheat genome), the TaFBX family harbors a different pattern at the various stratum of observation (subgenome, chromosome, chromosomal regions). We report here that the local gene expansion of the TaFBX family must be the consequence of multiple and complex events, including tandem and small-scale duplications. Regarding the differentially expressed TaFBX genes, while the majority of the genes are localized in the distal chromosomal regions (R1 and R3), differentially expressed genes are more present in the interstitial regions (R2a and R2b) than expected, which could be an indication of the preservation of major genes in those specific chromosomal regions.

New insights into homoeologous copy number variations in the hexaploid wheat genome

Bread wheat is an allohexaploid species originating from two successive and recent rounds of hybridization between three diploid species that were very similar in terms of chromosome number, genome size, TE content, gene content and synteny. As a result, it has long been considered that most of the genes were in three pairs of homoeologous copies. However, these so-called triads represent only one half of wheat genes, while the remaining half belong to homoeologous groups with various number of copies across subgenomes. In this study, we examined and compared the distribution, conservation, function, expression and epigenetic profiles of triads with homoeologous groups having undergone a deletion (dyads) or a duplication (tetrads) in one subgenome. We show that dyads and tetrads are mostly located in distal regions and have lower expression level and breadth than triads. Moreover, they are enriched in functions related to adaptation and more associated with the repressive H3K27me3 modification. Altogether, these results suggest that triads mainly correspond to housekeeping genes and are part of the core genome, while dyads and tetrads belong to the Triticeae dispensable genome. In addition, by comparing the different categories of dyads and tetrads, we hypothesize that, unlike most of the allopolyploid species, subgenome dominance and biased fractionation are absent in hexaploid wheat. Differences observed between the three subgenomes are more likely related to two successive and ongoing waves of post-polyploid diploidization, that had impacted A and B more significantly than D, as a result of the evolutionary history of hexaploid wheat.

Comparative Analysis of the Glutathione S-Transferase Gene Family of Four Triticeae Species and Transcriptome Analysis of GST Genes in Common Wheat Responding to Salt Stress

Glutathione S-transferases (GSTs) are ancient proteins encoded by a large gene family in plants, which play multiple roles in plant growth and development. However, there has been little study on the GST genes of common wheat (Triticum aestivum) and its relatives (Triticum durum, Triticum urartu, and Aegilops tauschii), which are four important species of Triticeae. Here, a genome-wide comprehensive analysis of this gene family was performed on the genomes of common wheat and its relatives. A total of 346 GST genes in T. aestivum, 226 in T. durum, 104 in T. urartu, and 105 in Ae. tauschii were identified, and all members were divided into ten classes. Transcriptome analysis was used to identify GST genes that respond to salt stress in common wheat, which revealed that the reaction of GST genes is not sensitive to low and moderate salt concentrations but is sensitive to severe concentrations of the stressor, and the GST genes related to salt stress mainly come from the Tau and Phi classes. Six GST genes which respond to different salt concentrations were selected and validated by a qRT-PCR assay. These findings will not only provide helpful information about the function of GST genes in Triticeae species but also offer insights for the future application of salt stress resistance breeding in common wheat.

Evaluation of genetic structure in European wheat cultivars and advanced breeding lines using high-density genotyping-by-sequencing approach

BACKGROUND

The genetic diversity and gene pool characteristics must be clarified for efficient genome-wide association studies, genomic selection, and hybrid breeding. The aim of this study was to evaluate the genetic structure of 509 wheat accessions representing registered varieties and advanced breeding lines via the high-density genotyping-by-sequencing approach.

RESULTS

More than 30% of 13,499 SNP markers representing 2162 clusters were mapped to genes, whereas 22.50% of 26,369 silicoDArT markers overlapped with coding sequences and were linked in 3527 blocks. Regarding hexaploidy, perfect sequence matches following BLAST searches were not sufficient for the unequivocal mapping to unique loci. Moreover, allelic variations in homeologous loci interfered with heterozygosity calculations for some markers. Analyses of the major genetic changes over the last 27 years revealed the selection pressure on orthologs of the gibberellin biosynthesis-related GA2 gene and the senescence-associated SAG12 gene. A core collection representing the wheat population was generated for preserving germplasm and optimizing breeding programs.

CONCLUSIONS

Our results confirmed considerable differences among wheat subgenomes A, B and D, with D characterized by the lowest diversity but the highest LD. They revealed genomic regions that have been targeted by breeding.

Large-Scale Cloning and Comparative Analysis of TaNAC Genes in Response to Stripe Rust and Powdery Mildew in Wheat (Triticum aestivum L.)

The NAM, ATAF1/2, and CUC2 (NAC) transcription factors (TFs) constitute the largest plant-specific TF superfamily, and play important roles in various physiological processes, including stress responses. Stripe rust and powdery mildew are the most damaging of the fungal diseases that afflict wheat (Triticum aestivum L.). However, studies on Triticum aestivum NAC (TaNAC)s’ role in resistance to the two diseases are still limited, especially in an overall comparative analysis of TaNACs responding or not to fungal stress. In the present study, 186 TaNAC transcripts were obtained from the resistant hexaploid wheat line N9134 under fungal stress, and 180 new transcripts were submitted to GenBank. Statistical results show that 35.1% (54/154) of TaNAC genes responded to stripe rust and powdery mildew in the seedling stage. “Abnormal” coding transcripts of differentially expressed (DE)-TaNAC genes in wheat responding to fungal stress were found in a significantly higher proportion (24/117 vs. 8/69, p = 0.0098) than in non-DE-NACs. This hinted that the alternative splicing of TaNAC genes was active in transcriptional or post-transcriptional regulation during plant-pathogen interactions. Full-length NAC proteins were classified into nine groups via phylogenetic analysis. Multiple-sequence alignment revealed diversity in the C-terminal structural organization, but the differentially expressed gene (DEG)-encoding proteins enriched in Subgroups VI and VII were conserved, with WV[L/V]CR amino acid residues in Motif 7 following the NAM domain. Our data that showed TaNAC TFs responded to fungal disease, which was affected by expression levels and by the regulation of multifarious transcript variants. These data for TaNAC responses to stripe rust and/or powdery mildew and their numerous structural variants provide a good resource for NAC function–mechanism analysis in the context of biotic-stress tolerance in wheat.

Identification and characterization of SET domain family genes in bread wheat (Triticum aestivum L.)

SET domain genes (SDGs) that are involved in histone methylation have been examined in many plant species, but have never been examined in bread wheat; the histone methylation caused due to SDGs is associated with regulation of gene expression at the transcription level. We identified a total of 166 bread wheat TaSDGs, which carry some interesting features including the occurrence of tandem/interspersed duplications, SSRs (simple sequence repeats), transposable elements, lncRNAs and targets for miRNAs along their lengths and transcription factor binding sites (TFBS) in the promoter regions. Only 130 TaSDGs encoded proteins with complete SET domain, the remaining 36 proteins had truncated SET domain. The TaSDG encoded proteins were classified into six classes (I–V and VII). In silico expression analysis indicated relatively higher expression (FPKM > 20) of eight of the 130 TaSDGs in different tissues, and downregulation of 30 TaSDGs under heat and drought at the seedling stage. qRT-PCR was also conducted to validate the expression of seven genes at the seedling stage in pairs of contrasting genotypes in response to abiotic stresses (water and heat) and biotic stress (leaf rust). These genes were generally downregulated in response to the three stresses examined.

Structural Variations Affecting Genes and Transposable Elements of Chromosome 3B in Wheats

Structural variations (SVs) such as copy number and presence–absence variations are polymorphisms that are known to impact genome composition at the species level and are associated with phenotypic variations. In the absence of a reference genome sequence, their study has long been hampered in wheat. The recent production of new wheat genomic resources has led to a paradigm shift, making possible to investigate the extent of SVs among cultivated and wild accessions. We assessed SVs affecting genes and transposable elements (TEs) in a Triticeae diversity panel of 45 accessions from seven tetraploid and hexaploid species using high-coverage shotgun sequencing of sorted chromosome 3B DNA and dedicated bioinformatics approaches. We showed that 23% of the genes are variable within this panel, and we also identified 330 genes absent from the reference accession Chinese Spring. In addition, 60% of the TE-derived reference markers were absent in at least one accession, revealing a high level of intraspecific and interspecific variability affecting the TE space. Chromosome extremities are the regions where we observed most of the variability, confirming previous hypotheses made when comparing wheat with the other grasses. This study provides deeper insights into the genomic variability affecting the complex Triticeae genomes at the intraspecific and interspecific levels and suggests a phylogeny with independent hybridization events leading to different hexaploid species.

Genome-wide identification and characterization of the fibrillin gene family in Triticum aestivum

Background

The fibrillin (FBN) gene family is highly conserved and widely distributed in the photosynthetic organs of plants. Members of this gene family are involved in the growth and development of plants and their response to biotic and abiotic stresses. Wheat (Triticum aestivum), an important food crop, has a complex genetic background and little progress has occurred in the understanding of its molecular mechanisms.

Methods

In this study, we identified 26 FBN genes in the whole genome of T. aestivum through bioinformatic tools and biotechnological means. These genes were divided into 11 subgroups and were distributed on 11 chromosomes of T. aestivum. Interestingly, most of the TaFBN genes were located on the chromosomes 2A, 2B and 2D. The gene structure of each subgroup of gene family members and the position and number of motifs were highly similar.

Results

The evolutionary analysis results indicated that the affinities of FBNs in monocots were closer together. The tissue-specific analysis revealed that TaFBN genes were expressed in different tissues and developmental stages. In addition, some TaFBNs were involved in one or more biotic and abiotic stress responses. These results provide a basis for further study of the biological function of FBNs.

Genome-wide characterization of the abscisic acid-, stress- and ripening-induced (ASR) gene family in wheat (Triticum aestivum L.)

BACKGROUND

Abscisic acid-, stress-, and ripening-induced (ASR) genes are a class of plant specific transcription factors (TFs), which play important roles in plant development, growth and abiotic stress responses. The wheat ASRs have not been described in genome-wide yet.

METHODS

We predicted the transmembrane regions and subcellular localization using the TMHMM server, and Plant-mPLoc server and CELLO v2.5, respectively. Then the phylogeny tree was built by MEGA7. The exon-intron structures, conserved motifs and TFs binding sites were analyzed by GSDS, MEME program and PlantRegMap, respectively.

RESULTS

In wheat, 33ASR genes were identified through a genome-wide survey and classified into six groups. Phylogenetic analyses revealed that the TaASR proteins in the same group tightly clustered together, compared with those from other species. Duplication analysis indicated that the TaASR gene family has expanded mainly through tandem and segmental duplication events. Similar gene structures and conserved protein motifs of TaASRs in wheat were identified in the same groups. ASR genes contained various TF binding cites associated with the stress responses in the promoter region. Gene expression was generally associated with the expected group-specific expression pattern in five tissues, including grain, leaf, root, spike and stem, indicating the broad conservation of ASR genes function during wheat evolution. The qRT-PCR analysis revealed that several ASRs were up-regulated in response to NaCl and PEG stress.

CONCLUSION

We identified ASR genes in wheat and found that gene duplication events are the main driving force for ASR gene evolution in wheat. The expression of wheat ASR genes was modulated in responses to multiple abiotic stresses, including drought/osmotic and salt stress. The results provided important information for further identifications of the functions of wheat ASR genes and candidate genes for high abiotic stress tolerant wheat breeding.

The default cyanobacterial linked genome: an interactive platform based on cyanobacterial linkage networks to assist functional genomics

A database of cyanobacterial linked genomes that can be accessed through an interactive platform (https://dfgm.ua.es/genetica/investigacion/cyanobacterial_genetics/Resources.html) was generated on the bases of conservation of gene neighborhood across 124 cyanobacterial species. It allows flexible generation of gene networks at different threshold values. The default cyanobacterial linked genome, whose global properties are analyzed here, connects most of the cyanobacterial core genes. The potential of the web tool is discussed in relation to other bioinformatics approaches based on guilty-by-association principles, with selected examples of networks illustrating its usefulness for genes found exclusively in cyanobacteria or in cyanobacteria and chloroplasts. We believe that this tool will provide useful predictions that are readily testable in Synechococcus elongatus PCC7942 and other model organisms performing oxygenic photosynthesis.

SNP markers for early identification of high molecular weight glutenin subunits (HMW-GSs) in bread wheat

A set of eight SNP markers was developed to facilitate the early selection of HMW-GS alleles in breeding programmes. In bread wheat (Triticum aestivum), the high molecular weight glutenin subunits (HMW-GSs) are the most important determinants of technological quality. Known to be very diverse, HMW-GSs are encoded by the tightly linked genes Glu-1-1 and Glu-1-2. Alleles that improve the quality of dough have been identified. Up to now, sodium dodecyl sulphate-polyacrylamide gel electrophoresis (SDS-PAGE) of grain proteins is the most widely used for their identification. To facilitate the early selection of HMW-GS alleles in breeding programmes, we developed DNA-based molecular markers. For each accession of a core collection (n = 364 lines) representative of worldwide bread wheat diversity, HMW-GSs were characterized by both genotyping and SDS-PAGE. Based on electrophoresis, we observed at least 8, 22 and 9 different alleles at the Glu-A1, Glu-B1 and Glu-D1 loci, respectively, including new variants. We designed a set of 17 single-nucleotide polymorphism (SNP) markers that were representative of the most frequent SDS-PAGE alleles at each locus. At Glu-A1 and Glu-D1, two and three marker-based haplotypes, respectively, captured the diversity of the SDS-PAGE alleles rather well. Discrepancies were found mainly for the Glu-B1 locus. However, statistical tests revealed that two markers at each Glu-B1 gene and their corresponding haplotypes were more significantly associated with the rheological properties of the dough than were the relevant SDS-PAGE alleles. To conclude, this study demonstrates that the SNP markers developed provide additional information on HMW-GS diversity. Two markers at Glu-A1, four at Glu-B1 and two at Glu-D1 constitute a useful toolbox for breeding wheat to improve end-use value.

Genotyping-by-Sequencing Derived Single Nucleotide Polymorphisms Provide the First Well-Resolved Phylogeny for the Genus Triticum (Poaceae)

Wheat (Triticum spp.) has been an important staple food crop for mankind since the beginning of agriculture. The genus Triticum L. is composed of diploid, tetraploid, and hexaploid species, majority of which have not yet been discriminated clearly, and hence their phylogeny and classification remain unresolved. Genotyping-by-sequencing (GBS) is an easy and affordable method that allows us to generate genome-wide single nucleotide polymorphism (SNP) markers. In this study, we used GBS to obtain SNPs covering all seven chromosomes from 283 accessions of Triticum-related genera. After filtering low-quality and redundant SNPs based on haplotype information, the GBS assay provided 14,188 high-quality SNPs that were distributed across the A (71%), B (26%), and D (2.4%) genomes. Cluster analysis and discriminant analysis of principal components (DAPC) allowed us to distinguish six distinct groups that matched well with Triticum species complexity. We constructed a Bayesian phylogenetic tree using 14,188 SNPs, in which 17 Triticum species and subspecies were discriminated. Dendrogram analysis revealed that the polyploid wheat species could be divided into groups according to the presence of A, B, D, and G genomes with strong nodal support and provided new insight into the evolution of spelt wheat. A total of 2,692 species-specific SNPs were identified to discriminate the common (T. aestivum) and durum (T. turgidum) wheat cultivar and landraces. In principal component analysis grouping, the two wheat species formed individual clusters and the SNPs were able to distinguish up to nine groups of 10 subspecies. This study demonstrated that GBS-derived SNPs could be used efficiently in genebank management to classify Triticum species and subspecies that are very difficult to distinguish by their morphological characters.

Genome-wide analysis of MIKC-type MADS-box genes in wheat: pervasive duplications, functional conservation and putative neofunctionalization

Wheat (Triticum aestivum) is one of the most important crops worldwide. Given a growing global population coupled with increasingly challenging cultivation conditions, facilitating wheat breeding by fine-tuning important traits is of great importance. MADS-box genes are prime candidates for this, as they are involved in virtually all aspects of plant development. Here, we present a detailed overview of phylogeny and expression of 201 wheat MIKC-type MADS-box genes. Homoeolog retention is significantly above the average genome-wide retention rate for wheat genes, indicating that many MIKC-type homoeologs are functionally important and not redundant. Gene expression is generally in agreement with the expected subfamily-specific expression pattern, indicating broad conservation of function of MIKC-type genes during wheat evolution. We also found extensive expansion of some MIKC-type subfamilies, especially those potentially involved in adaptation to different environmental conditions like flowering time genes. Duplications are especially prominent in distal telomeric regions. A number of MIKC-type genes show novel expression patterns and respond, for example, to biotic stress, pointing towards neofunctionalization. We speculate that conserved, duplicated and neofunctionalized MIKC-type genes may have played an important role in the adaptation of wheat to a diversity of conditions, hence contributing to the importance of wheat as a global staple food.

Identification and validation of a major and stably expressed QTL for spikelet number per spike in bread wheat

A major and stably expressed QTL for spikelet number per spike identified in a 2-cM interval on chromosome arm 2DS was validated using two populations with different genetic backgrounds. Spikelet number per spike (SNS) plays a key role in wheat yield improvement. Numerous genetic and environmental factors influencing SNS have been documented, but the number of major, stably expressed and validated loci underlying SNS is still limited. In this study, a recombinant inbred line (RIL) population derived from a normal spikelet cultivar and a multiple-spikelet wheat line (with a longer spike with more canonically oriented apical spikelets) was genotyped using a Wheat55K single-nucleotide polymorphism (SNP) array and simple sequence repeat (SSR) markers. SNS was measured for this RIL population in eight environments. Five QTL were each identified in two or more environments. One of them, QSns.sau-2D (LOD = 3.47-38.24, PVE = 10.16-45.68%), was detected in all the eight environments. The QTL was located in a 2-cM interval on chromosome arm 2DS flanked by the markers AX-109836946 and AX-111956072. This QTL, QSns.sau-2D, significantly increased SNS by up to 14.72%. Several genes associated with plant growth and development were identified in the physical interval of QSns.sau-2D. This QTL was further validated by the tightly linked Kompetitive Allele Specific PCR (KASP) marker, KASP-AX-94721936, in two other populations with different genetic backgrounds. The significant correlation between SNS and anthesis date, plant height, spike length, grain number per spike and thousand-grain weight were detected and discussed. These results lay the foundation for fine mapping and cloning gene(s) underlying QSns.sau-2D.

Serpins: Genome-Wide Characterisation and Expression Analysis of the Serine Protease Inhibitor Family in Triticum aestivum

The serine protease inhibitor (serpin) gene family is the largest family of protease inhibitors. Serine protease inhibitors have an active, but under-characterized, role in grain development and defense against pathogen attack in cereal crops. By exploiting publicly available genomic, transcriptomic and proteomic data for wheat (Triticum aestivum), we have identified and annotated the entire ’serpinome’ of wheat and constructed a high-quality and robust phylogenetic tree of the gene family, identifying paralogous and homeologous clades from the hexaploid wheat genome, including the Serpin-Z group that have been well characterized in barley. Using publicly available RNAseq data (http://www.wheat-expression.com/), expression profiles of the wheat serpins were explored across a variety of tissues from the developing grain, spikelet and spike. We show that the SERPIN-Z clade, among others, are highly expressed during grain development, and that there is homeologous and paralogous functional redundancy in this gene family. Further to their role in grain development, serpins play an important but under-explored role in response to fungal pathogens. Using 13 RNAseq datasets of wheat tissues infected by fungal pathogens, we identified 37 serpins with a significant disease response. The majority of the disease-responsive serpins were upregulated by Fusarium graminearum, a destructive fungal pathogen that attacks the spike and developing grain of wheat. As serpins are ubiquitous in wheat grain, the genes encoding serpins may be linked to grain development, with their disease response a result of pleiotropy.

Genome-Wide Association Study of Yield Component Traits in Intermediate Wheatgrass and Implications in Genomic Selection and Breeding

Intermediate wheatgrass (Thinopyrum intermedium, IWG) is a perennial grain crop with high biomass and grain yield, long seeds, and resistance to pests and diseases. It also reduces soil erosion, nitrate and mineral leaching into underground water tables, and sequesters carbon in its roots. The domestication timeline of IWG as a grain crop spans only 3 decades, hence it lags annual grain crops in yield and seed characteristics. One approach to improve its agronomic traits is by using molecular markers to uncover marker-trait associations. In this study, we performed association mapping on IWG breeding germplasm from the third recurrent selection cycle at the University of Minnesota. The IWG population was phenotyped in St Paul, MN in 2017 and 2018, and in Crookston, MN in 2018 for grain yield, seed length, width and weight, spike length and weight, and number of spikelets per spike. Strong positive correlations were observed among most trait pairs, with correlations as high as 0.76. Genotyping using high throughput sequencing identified 8,899 high-quality genome-wide SNPs which were combined with phenotypic data in association mapping to discover regions associated with the yield component traits. We detected 154 genetic loci associated with these traits of which 19 were shared between at least two traits. Prediction of breeding values using significant loci as fixed effects in genomic selection model improved predictive abilities by up to 14%. Genetic mapping of agronomic traits followed by using genomic selection to predict breeding values can assist breeders in selecting superior genotypes to accelerate IWG domestication.

Expansion and Functional Diversification of SKP1-Like Genes in Wheat (Triticum aestivum L.)

The ubiquitin proteasome 26S system (UPS), involving monomeric and multimeric E3 ligases is one of the most important signaling pathways in many organisms, including plants. The SCF (SKP1/Cullin/F-box) multimeric complex is particularly involved in response to development and stress signaling. The SKP1 protein (S-phase kinase-associated protein 1) is the core subunit of this complex. In this work, we firstly identified 92 and 87 non-redundant Triticum aestivumSKP1-like (TaSKP) genes that were retrieved from the latest release of the wheat genome database (International Wheat Genome Sequencing Consortium (IWGSC) RefSeq v1.0) and the genome annotation of the TGAC v1 respectively. We then investigated the structure, phylogeny, duplication events and expression patterns of the SKP1-like gene family in various tissues and environmental conditions using a wheat expression platform containing public data. TaSKP1-like genes were expressed differentially in response to stress conditions, displaying large genomic variations or short insertions/deletions which suggests functional specialization within TaSKP1-like genes. Finally, interactions between selected wheat FBX (F-box) proteins and putative ancestral TaSKP1-like proteins were tested using the yeast two-hybrid (Y2H) system to examine the molecular interactions. These observations suggested that six Ta-SKP1 genes are likely to be ancestral genes, having similar functions as ASK1 and ASK2 in Arabidopsis, OSK1 and OSK20 in rice and PpSKP1 and PpSKP2 in Physcomitrella patens.

Genome-wide analysis, expansion and expression of the NAC family under drought and heat stresses in bread wheat (T. aestivum L.)

The NAC family is one of the largest plant-specific transcription factor families, and some of its members are known to play major roles in plant development and response to biotic and abiotic stresses. Here, we inventoried 488 NAC members in bread wheat (Triticum aestivum). Using the recent release of the wheat genome (IWGS RefSeq v1.0), we studied duplication events focusing on genomic regions from 4B-4D-5A chromosomes as an example of the family expansion and neofunctionalization of TaNAC members. Differentially expressed TaNAC genes in organs and in response to abiotic stresses were identified using publicly available RNAseq data. Expression profiling of 23 selected candidate TaNAC genes was studied in leaf and grain from two bread wheat genotypes at two developmental stages in field drought conditions and revealed insights into their specific and/or overlapping expression patterns. This study showed that, of the 23 TaNAC genes, seven have a leaf-specific expression and five have a grain-specific expression. In addition, the grain-specific genes profiles in response to drought depend on the genotype. These genes may be considered as potential candidates for further functional validation and could present an interest for crop improvement programs in response to climate change. Globally, the present study provides new insights into evolution, divergence and functional analysis of NAC gene family in bread wheat.

Further studies on sugar transporter (SWEET) genes in wheat (Triticum aestivum L.)

SWEET proteins represent one of the largest sugar transporter family in the plant kingdom and play crucial roles in plant development and stress responses. In the present study, a total of 108 TaSWEET genes distributed on all the 21 wheat chromosomes were identified using the latest whole genome sequence (as against 59 genes reported in an earlier report). These 108 genes included 14 of the 17 types reported in Arabidopsis and also included three novel types. Tandem duplications (22) and segmental duplications (5) played a significant role in the expansion of TaSWEET family. A number of cis-elements were also identified in the promoter regions of TaSWEET genes, indicating response of TaSWEET genes during development and also during biotic/abiotic stresses. The TaSWEET proteins carried 4-7 trans-membrane helices (TMHs) showing diversity in structure. Phylogenetic analysis using SWEET proteins of wheat and 8 other species gave four well-known clusters. Expression analysis involving both in silico and in planta indicated relatively higher expression of TaSWEET genes in water/heat sensitive and leaf rust resistant genotypes. The results provided insights into the functional role of TaSWEETs in biotic and abiotic stresses, which may further help in planning strategies to develop high yielding wheat varieties tolerant to environmental stresses.

Identification of COS markers specific for Thinopyrum elongatum chromosomes preliminary revealed high level of macrosyntenic relationship between the wheat and Th. elongatum genomes

Thinopyrum elongatum (Host) D.R. Dewey has served as an important gene source for wheat breeding improvement for many years. The exact characterization of its chromosomes is important for the detailed analysis of prebreeding materials produced with this species. The major aim of this study was to identify and characterize new molecular markers to be used for the rapid analysis of E genome chromatin in wheat background. Sixty of the 169 conserved orthologous set (COS) markers tested on diverse wheat-Th. elongatum disomic/ditelosomic addition lines were assigned to various Th. elongatum chromosomes and will be used for marker-assisted selection. The macrosyntenic relationship between the wheat and Th. elongatum genomes was investigated using EST sequences. Several rearrangements were revealed in homoeologous chromosome groups 2, 5, 6 and 7, while chromosomes 1 and 4 were conserved. Molecular cytogenetic and marker analysis showed the presence of rearranged chromosome involved in 6ES and 2EL arms in the 6E disomic addition line. The selected chromosome arm-specific COS markers will make it possible to identify gene introgressions in breeding programmes and will also be useful in the development of new chromosome-specific markers, evolutionary analysis and gene mapping.

Genome-Wide Identification and Transcriptional Regulation of Aquaporin Genes in Bread Wheat (Triticum aestivum L.) under Water Stress

Aquaporins (AQPs) are transmembrane proteins essential for controlling the flow of water and other molecules required for development and stress tolerance in plants, including important crop species such as wheat (Triticum aestivum). In this study, we utilized a genomic approach for analyzing the information about AQPs available in public databases to characterize their structure and function. Furthermore, we validated the expression of a suite of AQP genes, at the transcriptional level, including accessions with contrasting responses to drought, different organs and water stress levels. We found 65 new AQP genes, from which 60% are copies expanded by polyploidization. Sequence analysis of the AQP genes showed that the purifying selection pressure acted on duplicate genes, which was related to a high conservation of the functions. This situation contrasted with the expression patterns observed for different organs, developmental stages or genotypes under water deficit conditions, which indicated functional divergence at transcription. Expression analyses on contrasting genotypes showed high gene transcription from Tonoplast Intrinsic Protein 1 (TIP1) and 2 (TIP2), and Plasma Membrane Intrinsic Protein 1 (PIP1) and 2 (PIP2) subfamilies in roots and from TIP1 and PIP1 subfamilies in leaves. Interestingly, during severe drought stress, 4 TIP genes analyzed in leaves of the tolerant accession reached up to 15-fold the level observed at the susceptible genotype, suggesting a positive relationship with drought tolerance. The obtained results extend our understanding of the structure and function of AQPs, particularly under water stress conditions.

Phylogenetic analyses and in-seedling expression of ammonium and nitrate transporters in wheat

Plants deploy several ammonium transporter (AMT) and nitrate transporter (NRT) genes to acquire NH₄⁺ and NO₃⁻ from the soil into the roots and then transport them to other plant organs. Coding sequences of wheat genes obtained from ENSEMBL were aligned to known AMT and NRT sequences of Arabidopsis, barley, maize, rice, and wheat to retrieve homologous genes. Bayesian phylogenetic relationships among these genes showed distinct classification of sequences with significant homology to NRT1, NRT2, and NRT3 (NAR2). Inter-species gene duplication analysis showed that eight AMT and 77 NRT genes were orthologous to the AMT and NRT genes of aforementioned plant species. Expression patterns of these genes were studied via whole transcriptome sequencing of 21-day old seedlings of five spring wheat lines. Eight AMT and 52 NRT genes were differentially expressed between root and shoot; and 131 genes did not express neither in root nor in shoot of 21-day old seedlings. Homeologous genes in the A, B, and D genomes, characterized by high sequence homology, revealed that their counterparts exhibited different expression patterns. This complement and evolutionary relationship of wheat AMT and NRT genes is expected to help in development of wheat germplasm with increased efficiency in nitrogen uptake and usage.

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

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

Identification and expression profiling of Pht1 phosphate transporters in wheat in controlled environments and in the field

Phosphorus (P) is an important macronutrient with critical functions in plants. Phosphate (Pi) transporters, which mediate Pi acquisition and Pi translocation within the plant, are key factors in Pi deficiency responses. However, their relevance for adaptation to long-term Pi limitation under agronomic conditions, particularly in wheat, remains unknown. Here, we describe the identification of the complete Pi transporter gene family (Pht1) in wheat (Triticum aestivum). Gene expression profiles were compared for hydroponic and field-grown plant tissues of wheat at multiple development stages. Cis-element analysis of selected Pht1 promoter regions was performed. A broad range of expression patterns of individual TaPht1 genes was observed in relation to tissue specificity and the nutrient supply in the soil or in liquid culture, as well as an influence of the experimental system. The expression patterns indicate the involvement of specific transporters in Pi uptake, and in Pi transport and remobilisation within the plant, at different growth developmental stages. Specifically, the influence of Pi nutrition indicates a complex regulatory pattern of TaPht1 gene transcript abundances as a response to low Pi availability in different culture systems, correlating with the existence of different cis-acting promoter elements.

High throughput SNP discovery and genotyping in hexaploid wheat

Because of their abundance and their amenability to high-throughput genotyping techniques, Single Nucleotide Polymorphisms (SNPs) are powerful tools for efficient genetics and genomics studies, including characterization of genetic resources, genome-wide association studies and genomic selection. In wheat, most of the previous SNP discovery initiatives targeted the coding fraction, leaving almost 98% of the wheat genome largely unexploited. Here we report on the use of whole-genome resequencing data from eight wheat lines to mine for SNPs in the genic, the repetitive and non-repetitive intergenic fractions of the wheat genome. Eventually, we identified 3.3 million SNPs, 49% being located on the B-genome, 41% on the A-genome and 10% on the D-genome. We also describe the development of the TaBW280K high-throughput genotyping array containing 280,226 SNPs. Performance of this chip was examined by genotyping a set of 96 wheat accessions representing the worldwide diversity. Sixty-nine percent of the SNPs can be efficiently scored, half of them showing a diploid-like clustering. The TaBW280K was proven to be a very efficient tool for diversity analyses, as well as for breeding as it can discriminate between closely related elite varieties. Finally, the TaBW280K array was used to genotype a population derived from a cross between Chinese Spring and Renan, leading to the construction a dense genetic map comprising 83,721 markers. The results described here will provide the wheat community with powerful tools for both basic and applied research.

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

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

Chromosomal assignment of centromere-specific histone CENH3 genes in rye (Secale cereale L.) and their phylogeny

Centromeres are essential for correct chromosome segregation during cell division and are determined by the presence of centromere-specific histone 3 (CENH3). Most of the diploid plant species, in which the structure and copy number of CENH3 genes have been determined, have this gene as a singleton; however, some cereal species in the tribe Triticeae have been found to have CENH3 in two variants. In this work, using the set of the wheat-rye addition lines we wanted to establish the chromosomal assignment of the CENH3 genes in the cultivated rye, Secale cereale (Linnaeus, 1753), in order to expand our knowledge about synteny conservation in the most important cereal species and about their chromosome evolution. To this end, we have also analyzed data in available genome sequencing databases. As a result, the αCENH3 and βCENH3 forms have been assigned to rye chromosomes 1R and 6R: specifically, the commonest variants αCENH3v1 and βCENH3v1 to chromosome 1R, and the rare variants, αCENH3v2 and probably βCENH3v2, to chromosome 6R. No other CENH3 variants have been found by analysis of the rye genome sequencing databases. Our chromosomal assignment of CENH3 in rye has been found to be the same as that in barley, suggesting that both main forms of CENH3 appeared in a Triticeae species before the barley and wheatrye lineages split.

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

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

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.

Genome-wide analysis of WRKY transcription factors in wheat (Triticum aestivum L.) and differential expression under water deficit condition

BACKGROUND

WRKY proteins, which comprise one of the largest transcription factor (TF) families in the plant kingdom, play crucial roles in plant development and stress responses. Despite several studies on WRKYs in wheat (Triticum aestivum L.), functional annotation information about wheat WRKYs is limited.

RESULTS

Here, 171 TaWRKY TFs were identified from the whole wheat genome and compared with proteins from 19 other species representing nine major plant lineages. A phylogenetic analysis, coupled with gene structure analysis and motif determination, divided these TaWRKYs into seven subgroups (Group I, IIa-e, and III). Chromosomal location showed that most TaWRKY genes were enriched on four chromosomes, especially on chromosome 3B. In addition, 85 (49.7%) genes were either tandem (5) or segmental duplication (80), which suggested that though tandem duplication has contributed to the expansion of TaWRKY family, segmental duplication probably played a more pivotal role. Analysis of cis-acting elements revealed putative functions of WRKYs in wheat during development as well as under numerous biotic and abiotic stresses. Finally, the expression of TaWRKY genes in flag leaves, glumes, and lemmas under water-deficit condition were analyzed. Results showed that different TaWRKY genes preferentially express in specific tissue during the grain-filling stage.

CONCLUSION

Our results provide a more extensive insight on WRKY gene family in wheat, and also contribute to the screening of more candidate genes for further investigation on function characterization of WRKYs under various stresses.

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.

VRN1 genes variability in tetraploid wheat species with a spring growth habit

BACKGROUND

Vernalization genes VRN1 play a major role in the transition from vegetative to reproductive growth in wheat. In di-, tetra- and hexaploid wheats the presence of a dominant allele of at least one VRN1 gene homologue (Vrn-A1, Vrn-B1, Vrn-G1 or Vrn-D1) determines the spring growth habit. Allelic variation between the Vrn-1 and vrn-1 alleles relies on mutations in the promoter region or the first intron. The origin and variability of the dominant VRN1 alleles, determining the spring growth habit in tetraploid wheat species have been poorly studied.

RESULTS

Here we analyzed the growth habit of 228 tetraploid wheat species accessions and 25 % of them were spring type. We analyzed the promoter and first intron regions of VRN1 genes in 57 spring accessions of tetraploid wheats. The spring growth habit of most studied spring accessions was determined by previously identified dominant alleles of VRN1 genes. Genetic experiments proof the dominant inheritance of Vrn-A1d allele which was widely distributed across the accessions of Triticum dicoccoides. Two novel alleles were discovered and designated as Vrn-A1b.7 and Vrn-B1dic. Vrn-A1b.7 had deletions of 20 bp located 137 bp upstream of the start codon and mutations within the VRN-box when compared to the recessive allele of vrn-A1. So far the Vrn-A1d allele was identified only in spring accessions of the T. dicoccoides and T. turgidum species. Vrn-B1dic was identified in T. dicoccoides IG46225 and had 11 % sequence dissimilarity in comparison to the promoter of vrn-B1. The presence of Vrn-A1b.7 and Vrn-B1dic alleles is a predicted cause of the spring growth habit of studied accessions of tetraploid species. Three spring accessions T. aethiopicum K-19059, T. turanicum K-31693 and T. turgidum cv. Blancal possess recessive alleles of both VRN-A1 and VRN-B1 genes. Further investigations are required to determine the source of spring growth habit of these accessions.

CONCLUSIONS

New allelic variants of the VRN-A1 and VRN-B1 genes were identified in spring accessions of tetraploid wheats. The origin and evolution of VRN-A1 alleles in di- and tetraploid wheat species was discussed.

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

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

Genomic Dissection and Expression Profiling Revealed Functional Divergence in Triticum aestivum Leucine Rich Repeat Receptor Like Kinases (TaLRRKs)

The leucine rich repeat receptor like kinases (LRRK) constitute the largest subfamily of receptor like kinases (RLK), which play critical roles in plant development and stress responses. Herein, we identified 531 TaLRRK genes in Triticum aestivum (bread wheat), which were distributed throughout the A, B, and D sub-genomes and chromosomes. These were clustered into 233 homologous groups, which were mostly located on either homeologous chromosomes from various sub-genomes or in proximity on the same chromosome. A total of 255 paralogous genes were predicted which depicted the role of duplication events in expansion of this gene family. Majority of TaLRRKs consisted of trans-membrane region and localized on plasma-membrane. The TaLRRKs were further categorized into eight phylogenetic groups with numerous subgroups on the basis of sequence homology. The gene and protein structure in terms of exon/intron ratio, domains, and motifs organization were found to be variably conserved across the different phylogenetic groups/subgroups, which indicated a potential divergence and neofunctionalization during evolution. High-throughput transcriptome data and quantitative real time PCR analyses in various developmental stages, and biotic and abiotic (heat, drought, and salt) stresses provided insight into modus operandi of TaLRRKs during these conditions. Distinct expression of majority of stress responsive TaLRRKs homologous genes suggested their specified role in a particular condition. These results provided a comprehensive analysis of various characteristic features including functional divergence, which may provide the way for future functional characterization of this important gene family in bread wheat.

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

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

Evolutionary Fates and Dynamic Functionalization of Young Duplicate Genes in Arabidopsis Genomes

Gene duplication is a primary means to generate genomic novelties, playing an essential role in speciation and adaptation. Particularly in plants, a high abundance of duplicate genes has been maintained for significantly long periods of evolutionary time. To address the manner in which young duplicate genes were derived primarily from small-scale gene duplication and preserved in plant genomes and to determine the underlying driving mechanisms, we generated transcriptomes to produce the expression profiles of five tissues in Arabidopsis thaliana and the closely related species Arabidopsis lyrata and Capsella rubella Based on the quantitative analysis metrics, we investigated the evolutionary processes of young duplicate genes in Arabidopsis. We determined that conservation, neofunctionalization, and specialization are three main evolutionary processes for Arabidopsis young duplicate genes. We explicitly demonstrated the dynamic functionalization of duplicate genes along the evolutionary time scale. Upon origination, duplicates tend to maintain their ancestral functions; but as they survive longer, they might be likely to develop distinct and novel functions. The temporal evolutionary processes and functionalization of plant duplicate genes are associated with their ancestral functions, dynamic DNA methylation levels, and histone modification abundances. Furthermore, duplicate genes tend to be initially expressed in pollen and then to gain more interaction partners over time. Altogether, our study provides novel insights into the dynamic retention processes of young duplicate genes in plant genomes.

Evolutionary Fates and Dynamic Functionalization of Young Duplicate Genes in Arabidopsis Genomes

Gene duplication is a primary means to generate genomic novelties, playing an essential role in speciation and adaptation. Particularly in plants, a high abundance of duplicate genes has been maintained for significantly long periods of evolutionary time. To address the manner in which young duplicate genes were derived primarily from small-scale gene duplication and preserved in plant genomes and to determine the underlying driving mechanisms, we generated transcriptomes to produce the expression profiles of five tissues in Arabidopsis thaliana and the closely related species Arabidopsis lyrata and Capsella rubella Based on the quantitative analysis metrics, we investigated the evolutionary processes of young duplicate genes in Arabidopsis. We determined that conservation, neofunctionalization, and specialization are three main evolutionary processes for Arabidopsis young duplicate genes. We explicitly demonstrated the dynamic functionalization of duplicate genes along the evolutionary time scale. Upon origination, duplicates tend to maintain their ancestral functions; but as they survive longer, they might be likely to develop distinct and novel functions. The temporal evolutionary processes and functionalization of plant duplicate genes are associated with their ancestral functions, dynamic DNA methylation levels, and histone modification abundances. Furthermore, duplicate genes tend to be initially expressed in pollen and then to gain more interaction partners over time. Altogether, our study provides novel insights into the dynamic retention processes of young duplicate genes in plant genomes.