Wheat syntenome unveils new evidences of contrasted evolutionary plasticity between paleo- and neoduplicated subgenomes

Evolutionary trajectory of organelle-derived nuclear DNAs in the Triticum/Aegilops complex species

Organelle-derived nuclear DNAs, nuclear plastid DNAs (NUPTs), and nuclear mitochondrial DNAs (NUMTs) have been identified in plants. Most, if not all, genes residing in NUPTs/NUMTs (NUPGs/NUMGs) are known to be inactivated and pseudogenized. However, the role of epigenetic control in silencing NUPGs/NUMGs and the dynamic evolution of NUPTs/NUMTs with respect to organismal phylogeny remain barely explored. Based on the available nuclear and organellar genomic resources of wheat (genus Triticum) and goat grass (genus Aegilops) within Triticum/Aegilops complex species, we investigated the evolutionary fates of NUPTs/NUMTs in terms of their epigenetic silencing and their dynamic occurrence rates in the nuclear diploid genomes and allopolyploid subgenomes. NUPTs and NUMTs possessed similar genomic atlas, including (i) predominantly located in intergenic regions and preferential integration to gene regulation regions and (ii) generating sequence variations in the nuclear genome. Unlike nuclear indigenous genes, the alien NUPGs/NUMGs were associated with repressive epigenetic signals, namely high levels of DNA methylation and low levels of active histone modifications. Phylogenomic analyses suggested that the species-specific and gradual accumulation of NUPTs/NUMTs accompanied the speciation processes. Moreover, based on further pan-genomic analyses, we found significant subgenomic asymmetry in the NUPT/NUMT occurrence, which accumulated during allopolyploid wheat evolution. Our findings provide insight into the dynamic evolutionary fates of organelle-derived nuclear DNA in plants.

Genome-wide analysis of transmembrane 9 superfamily genes in wheat (Triticum aestivum) and their expression in the roots under nitrogen limitation and Bacillus amyloliquefaciens PDR1 treatment conditions

Introduction

Transmembrane 9 superfamily (TM9SF) proteins play significant roles in plant physiology. However, these proteins are poorly characterized in wheat (Triticum aestivum). The present study aimed at the genome-wide analysis of putative wheat TM9SF (TraesTM9SF) proteins and their potential involvement in response to nitrogen limitation and Bacillus amyloliquefaciens PDR1 treatments.

Methods

TraesTM9SF genes were retrieved from the wheat genome, and their physiochemical properties, alignment, phylogenetic, motif structure, cis-regulatory element, synteny, protein-protein interaction (PPI), and transcription factor (TF) prediction analyses were performed. Transcriptome sequencing and quantitative real-time polymerase reaction (qRT-PCR) were performed to detect gene expression in roots under single or combined treatments with nitrogen limitation and B. amyloliquefaciens PDR1.

Results and discussion

Forty-seven TraesTM9SF genes were identified in the wheat genome, highlighting the significance of these genes in wheat. TraesTM9SF genes were absent on some wheat chromosomes and were unevenly distributed on the other chromosomes, indicating that potential regulatory functions and evolutionary events may have shaped the TraesTM9SF gene family. Fifty-four cis-regulatory elements, including light-response, hormone response, biotic/abiotic stress, and development cis-regulatory elements, were present in the TraesTM9SF promoter regions. No duplication of TraesTM9SF genes in the wheat genome was recorded, and 177 TFs were predicted to target the 47 TraesTM9SF genes in a complex regulatory network. These findings offer valued data for predicting the putative functions of uncharacterized TM9SF genes. Moreover, transcriptome analysis and validation by qRT-PCR indicated that the TraesTM9SF genes are expressed in the root system of wheat and are potentially involved in the response of this plant to single or combined treatments with nitrogen limitation and B. amyloliquefaciens PDR1, suggesting their functional roles in plant growth, development, and stress responses.

Conclusion

These findings may be vital in further investigation of the function and biological applications of TM9SF genes in wheat.

Phased gap-free genome assembly of octoploid cultivated strawberry illustrates the genetic and epigenetic divergence among subgenomes

The genetic and epigenetic mechanisms underlying the coexistence and coordination of the four diverged subgenomes (ABCD) in octoploid strawberries (Fragaria × ananassa) remains poorly understood. In this study, we have assembled a haplotype-phased gap-free octoploid genome for the strawberry, which allowed us to uncover the sequence, structure, and epigenetic divergences among the subgenomes. The diploid progenitors of the octoploid strawberry, apart from subgenome A (Fragaria vesca), have been a subject of public controversy. Phylogenomic analyses revealed a close relationship between diploid species Fragaria iinumae and subgenomes B, C, and D. Subgenome A, closely related to F. vesca, retains the highest number of genes, exhibits the lowest content of transposable elements (TEs), experiences the strongest purifying selection, shows the lowest DNA methylation levels, and displays the highest expression level compared to the other three subgenomes. Transcriptome and DNA methylome analyses revealed that subgenome A-biased genes were enriched in fruit development biological processes. In contrast, although subgenomes B, C, and D contain equivalent amounts of repetitive sequences, they exhibit diverged methylation levels, particularly for TEs located near genes. Taken together, our findings provide valuable insights into the evolutionary patterns of subgenome structure, divergence and epigenetic dynamics in octoploid strawberries, which could be utilized in strawberry genetics and breeding research.

Genome shock in a synthetic allotetraploid wheat invokes subgenome-partitioned gene regulation, meiotic instability, and karyotype variation

It is becoming increasingly evident that interspecific hybridization at the homoploid level or coupled with whole-genome duplication (i.e. allopolyploidization) has played a major role in biological evolution. However, the direct impacts of hybridization and allopolyploidization on genome structure and function, phenotype, and fitness remains to be fully understood. Synthetic hybrids and allopolyploids are trackable experimental systems that can be used to address this issue. In this study, we resynthesized a pair of reciprocal F1 hybrids and corresponding reciprocal allotetraploids using the two diploid progenitor species of bread wheat (Triticum aestivum, BBAADD), namely T. urartu (AA) and Aegilops tauschii (DD). By comparing phenotypes related to growth, development, and fitness, and by analysing genome expression in both hybrids and allotetraploids in relation to the parents, we found that the types and trends of karyotype variation in the immediately formed allotetraploids were correlated with both instability of meiosis and chromosome- and subgenome-biased expression. We determined clear advantages of allotetraploids over diploid F1 hybrids in several morphological traits including fitness that mirrored the tissue- and developmental stage-dependent subgenome-partitioning of the allotetraploids. The allotetraploids were meiotically unstable primarily due to homoeologous pairing that varied dramatically among the chromosomes. Nonetheless, the manifestation of organismal karyotype variation and the occurrence of meiotic irregularity were not concordant, suggesting a role of functional constraints probably imposed by subgenome- and chromosome-biased gene expression. Our results provide new insights into the direct impacts and consequences of hybridization and allopolyploidization that are relevant to evolution and likely to be informative for future crop improvement approaches using synthetic polyploids.

Molecular Evolution and Genetic Variation of G2-Like Transcription Factor Genes in Wheat (Triticum aestivum L.)

The GOLDEN2-LIKE (G2-like) gene family members provide significant contributions to the growth and development of plants. In this study, a total of 76 wheat G2-like gene family members (TaG1-TaG76) were detected in the wheat genome and were categorized into three groups (including six subgroups) based on the gene structure and protein motif analyses. These genes were unevenly distributed in 19 of 21 wheat chromosomes. A total of 63 segmental duplication pairs of TaG2-like genes were identified in the wheat genome. The expression levels of all the TaG2-like genes indicated that TaG2-like genes showed different expression patterns in various organs and tissues. Moreover, the transcriptions of TaG2-like genes were significantly affected under abiotic stress (cold, ABA, NaCl, and PEG). This study offered valuable insights into the functional characterization of TaG2-like genes in wheat.

Tracing 100 million years of grass genome evolutionary plasticity

Grasses derive from a family of monocotyledonous plants that includes crops of major economic importance such as wheat, rice, sorghum and barley, sharing a common ancestor some 100 million years ago. The genomic attributes of plant adaptation remain obscure and the consequences of recurrent whole genome duplications (WGD) or polyploidization events, a major force in plant evolution, remain largely speculative. We conducted a comparative analysis of omics data from ten grass species to unveil structural (inversions, fusions, fissions, duplications, substitutions) and regulatory (expression and methylation) basis of genome plasticity, as possible attributes of plant long lasting evolution and adaptation. The present study demonstrates that diverged polyploid lineages sharing a common WGD event often present the same patterns of structural changes and evolutionary dynamics, but these patterns are difficult to generalize across independent WGD events as a result of non-WGD factors such as selection and domestication of crops. Polyploidy is unequivocally linked to the evolutionary success of grasses during the past 100 million years, although it remains difficult to attribute this success to particular genomic consequences of polyploidization, suggesting that polyploids harness the potential of genome duplication, at least partially, in lineage-specific ways. Overall, the present study clearly demonstrates that post-polyploidization reprogramming is more complex than traditionally reported in investigating single species and calls for a critical and comprehensive comparison across independently polyploidized lineages.

The genome of Orychophragmus violaceus provides genomic insights into the evolution of Brassicaceae polyploidization and its distinct traits

Orychophragmus violaceus, referred to as "eryuelan" (February orchid) in China, is an early-flowering ornamental plant. The high oil content and abundance of unsaturated fatty acids in O. violaceus seeds make it a potential high-quality oilseed crop. Here, we generated a whole-genome assembly for O. violaceus using Nanopore and Hi-C sequencing technologies. The assembled genome of O. violaceus was ∼1.3 Gb in size, with 12 pairs of chromosomes. Through investigation of ancestral genome evolution, we determined that the genome of O. violaceus experienced a tetraploidization event from a diploid progenitor with the translocated proto-Calepineae karyotype. Comparisons between the reconstructed subgenomes of O. violaceus identified indicators of subgenome dominance, indicating that subgenomes likely originated via allotetraploidy. O. violaceus was phylogenetically close to the Brassica genus, and tetraploidy in O. violaceus occurred approximately 8.57 million years ago, close in time to the whole-genome triplication of Brassica that likely arose via an intermediate tetraploid lineage. However, the tetraploidization in Orychophragmus was independent of the hexaploidization in Brassica, as evidenced by the results from detailed phylogenetic analyses and comparisons of the break and fusion points of ancestral genomic blocks. Moreover, identification of multi-copy genes regulating the production of high-quality oil highlighted the contributions of both tetraploidization and tandem duplication to functional innovation in O. violaceus. These findings provide novel insights into the polyploidization evolution of plant species and will promote both functional genomic studies and domestication/breeding efforts in O. violaceus.

The suppression of TdMRP3 genes reduces the phytic acid and increases the nutrient accumulation in durum wheat grain

Micronutrient malnutrition affects more than half of the world population. Reduced bioavailability of microelements in the raw materials is considered one of the main causes of mineral deficiency in populations whose diet is largely based on the consumption of staple crops. In this context, the production of low phytic acid (lpa) cereals is a main goal of the breeding programs, as phytic acid (PA) binds essential mineral cations such as iron (Fe), zinc (Zn), manganese (Mn), potassium (K), calcium (Ca) and magnesium (Mg) precipitating in the form of phytate salts poorly digested by monogastric animals, including humans, due to the lack of phytases in the digestive tract. Since PA limits the bioavailability of microelements, it is widely recognized as an anti-nutritional compound. A Targeting Induced Local Lesions IN Genomes (TILLING) approach has been undertaken to silence the genes encoding the TdABCC13 proteins, known as Multidrug-Resistance associated Proteins 3 (TdMRP3), transporters involved in the accumulation of PA inside the vacuole in durum wheat. The TdMRP3 complete null genotypes showed a significant reduction in the content of PA and were able to accumulate a higher amount of essential micronutrients (Fe, Zn, Mn) compared to the control. The number of spikelets and seeds per spike, traits associated with the agronomic performances, were reduced compared to the control, but the negative effect was in part balanced by the increased grain weight. The TdMRP3 mutant lines showed morphological differences in the root apparatus such as a significant decrease in the number of root tips, root length, volume and surface area and an increase in root average diameter compared to the control plants. These materials represent a promising basis for obtaining new commercial durum wheats with higher nutritional value.

The Heterogeneity in the Landscape of Gene Dominance in Maize is Accompanied by Unique Chromatin Environments

Subgenome dominance after whole-genome duplication (WGD) has been observed in many plant species. However, the degree to which the chromatin environment affects this bias has not been explored. Here, we compared the dominant subgenome (maize1) and the recessive subgenome (maize2) with respect to patterns of sequence substitutions, genes expression, transposable element accumulation, small interfering RNAs, DNA methylation, histone modifications, and accessible chromatin regions (ACRs). Our data show that the degree of bias between subgenomes for all the measured variables does not vary significantly when both of the WGD genes are located in pericentromeric regions. Our data further indicate that the location of maize1 genes in chromosomal arms is pivotal for maize1 to maintain its dominance, but location has a less effect on maize2 homoeologs. In addition to homoeologous genes, we compared ACRs, which often harbor cis-regulatory elements, between the two subgenomes and demonstrate that maize1 ACRs have a higher level of chromatin accessibility, a lower level of sequence substitution, and are enriched in chromosomal arms. Furthermore, we find that a loss of maize1 ACRs near their nearby genes is associated with a reduction in purifying selection and expression of maize1 genes relative to their maize2 homoeologs. Taken together, our data suggest that chromatin environment and cis-regulatory elements are important determinants shaping the divergence and evolution of duplicated genes.

Subgenome dominance and its evolutionary implications in crop domestication and breeding

Polyploidization or whole-genome duplication (WGD) is a well-known speciation and adaptation mechanism in angiosperms, while subgenome dominance is a crucial phenomenon in allopolyploids, established following polyploidization. The dominant subgenomes contribute more to genome evolution and homoeolog expression bias, both of which confer advantages for short-term phenotypic adaptation and long-term domestication. In this review, we firstly summarize the probable mechanistic basis for subgenome dominance, including the effects of genetic [transposon, genetic incompatibility, and homoeologous exchange (HE)], epigenetic (DNA methylation and histone modification), and developmental and environmental factors on this evolutionary process. We then move to Brassica rapa, a typical allopolyploid with subgenome dominance. Polyploidization provides the B. rapa genome not only with the genomic plasticity for adapting to changeable environments, but also an abundant genetic basis for morphological variation, making it a representative species for subgenome dominance studies. According to the 'two-step theory', B. rapa experienced genome fractionation twice during WGD, in which most of the genes responding to the environmental cues and phytohormones were over-retained, enhancing subgenome dominance and consequent adaption. More than this, the pangenome of 18 B. rapa accessions with different morphotypes recently constructed provides further evidence to reveal the impacts of polyploidization and subgenome dominance on intraspecific diversification in B. rapa. Above and beyond the fundamental understanding of WGD and subgenome dominance in B. rapa and other plants, however, it remains elusive why subgenome dominance has tissue- and spatiotemporal-specific features and could shuffle between homoeologous regions of different subgenomes by environments in allopolyploids. We lastly propose acceleration of the combined application of resynthesized allopolyploids, omics technology, and genome editing tools to deepen mechanistic investigations of subgenome dominance, both genetic and epigenetic, in a variety of species and environments. We believe that the implications of genomic and genetic basis of a variety of ecologically, evolutionarily, and agriculturally interesting traits coupled with subgenome dominance will be uncovered and aid in making new discoveries and crop breeding.

Spatial and Temporal Transcriptomic Heredity and Asymmetry in an Artificially Constructed Allotetraploid Wheat (AADD)

Polyploidy, or whole-genome duplication (WGD), often induces dramatic changes in gene expression due to "transcriptome shock. " However, questions remain about how allopolyploidy (the merging of multiple nuclear genomes in the same nucleus) affects gene expression within and across multiple tissues and developmental stages during the initial foundation of allopolyploid plants. Here, we systematically investigated the immediate effect of allopolyploidy on gene expression variation in an artificial allopolyploidy system consisting of a constructed allotetraploid wheat (AADD genome, accession AT2) and its diploid progenitors Triticum urartu and Aegilops tauschii. We performed comprehensive RNA sequencing of 81 samples from different genotypes, tissues, and developmental stages. First, we found that intrinsic interspecific differences between the diploid parents played a major role in establishing the expression architecture of the allopolyploid. Nonetheless, allopolyploidy per se also induced dramatic and asymmetric patterns of differential gene expression between the subgenomes, and genes from the D subgenome exhibited a more drastic response. Second, analysis of homoeolog expression bias (HEB) revealed that the D subgenome exhibited significant expression bias and that de novo-generated HEB was attributed mainly to asymmetrical differential gene expression. Homoeolog-specific expression (HSE) analyses showed that the cis-only regulatory pattern was predominant in AT2, reflecting significant divergence between the parents. Co-expression network analysis revealed that homoeolog expression connectivity (HEC) was significantly correlated with sequence divergence in cis elements between subgenomes. Interestingly, allopolyploidy-induced reconstruction of network modules was also associated with different HSE patterns. Finally, a transcriptome atlas of spike development demonstrated that the phenotypic similarity of AT2 to T. urartu may be attributed to the combination of relatively stable expression of A-subgenome genes and drastic downregulation of their D-subgenome homoeologs. These findings provide a broad, multidimensional characterization of allopolyploidy-induced transcriptomic responses and suggest that allopolyploidy can have immediate and complex regulatory effects on the expression of nuclear genes.

Key auxin response factor (ARF) genes constraining wheat tillering of mutant dmc

Tillering ability is a key agronomy trait for wheat (Triticum aestivum L.) production. Studies on a dwarf monoculm wheat mutant (dmc) showed that ARF11 played an important role in tillering of wheat. In this study, a total of 67 ARF family members were identified and clustered to two main classes with four subgroups based on their protein structures. The promoter regions of T. aestivum ARF (TaARF) genes contain a large number of cis-acting elements closely related to plant growth and development, and hormone response. The segmental duplication events occurred commonly and played a major role in the expansion of TaARFs. The gene collinearity degrees of the ARFs between wheat and other grasses, rice and maize, were significantly high. The evolution distances among TaARFs determine their expression profiles, such as homoeologous genes have similar expression profiles, like TaARF4-3A-1, TaARF4-3A-2 and their homoeologous genes. The expression profiles of TaARFs in various tissues or organs indicated TaARF3, TaARF4, TaARF9 and TaARF22 and their homoeologous genes played basic roles during wheat development. TaARF4, TaARF9, TaARF12, TaARF15, TaARF17, TaARF21, TaARF25 and their homoeologous genes probably played basic roles in tiller development. qRT-PCR analyses of 20 representative TaARF genes revealed that the abnormal expressions of TaARF11 and TaARF14 were major causes constraining the tillering of dmc. Indole-3-acetic acid (IAA) contents in dmc were significantly less than that in Guomai 301 at key tillering stages. Exogenous IAA application significantly promoted wheat tillering, and affected the transcriptions of TaARFs. These data suggested that TaARFs as well as IAA signaling were involved in controlling wheat tillering. This study provided valuable clues for functional characterization of ARFs in wheat.

Chromosome karyotype and stability of new synthetic hexaploid wheat

Synthetic hexaploid wheat offers breeders ready access to potentially novel genetic variation in wild ancestral species. In this study, we crossed MY3478 (2n = 4x = 28, AABB) as the maternal parent with the stripe rust-resistant SY41 (2n = 2x = 14, DD) as the paternal parent to construct the new hexaploid wheat line NA0928 through natural allopolyploidization. Agronomic traits and the cytology of the S8-S9 generations of NA0928 were analyzed. Abundant variation in agronomic traits was observed among each strain of NA0928 in the S8 generation. Agronomic traits were superior in strains resistant to stripe rust compared with those of highly susceptible strains. The rank order of the coefficients of variation were tiller number (55.3%) > spike length (15.3%) > number of spikelets (13.9%) > plant height (8.7). Number of tillers and spike length are important traits in wheat breeding to improve yield. Cytological observation and fluorescence in situ hybridization showed that the chromosome number and configuration showed rich variation among NA0928 strains in the S9 generation. Chromosome number ranged from 36 to 44. Variation in chromosome karyotype was detected in the A and B subgenomes. Meiotic chromosome behavior in pollen mother cells and multicolor genomic in situ hybridization revealed that two new synthetic hexaploid wheat strains showed genetic stability; one strain was resistant to stripe rust and developed multiple tillers, and the other strain was susceptible to stripe rust, but both showed improved thousand-kernel weight (TKW) weight and produced multiple tillers. The two strains will be valuable germplasm resources for use in wheat breeding.

Karyotype Reorganization in Wheat-Rye Hybrids Obtained via Unreduced Gametes: Is There a Limit to the Chromosome Number in Triticale?

To date, few data have been accumulated on the contribution of meiotic restitution to the formation of Triticum aestivum hybrid karyotypes. In this study, based on FISH and C-banding, karyotype reorganization was observed in three groups of F₅ wheat–rye hybrids 1R(1A) × R. Aberrations, including aneuploidy, telocentrics, and Robertsonian translocations, were detected in all groups. Some of the Group 1 plants and all of the Group 2 plants only had a 4R4R pair (in addition to 1R1R), which was either added or substituted for its homeolog in ABD subgenomes. In about 82% of meiocytes, 4R4R formed bivalents, which indicates its competitiveness. The rest of the Group 1 plants had 2R and 7R chromosomes in addition to 1R1R. Group 3 retained all their rye chromosomes, with a small aneuploidy on the wheat chromosomes. A feature of the meiosis in the Group 3 plants was asynchronous cell division and omission of the second division. Diploid gametes did not form because of the significant disturbances during gametogenesis. As a result, the frequency of occurrence of the formed dyads was negatively correlated (r = −0.73) with the seed sets. Thus, meiotic restitution in the 8n triticale does not contribute to fertility or increased ploidy in subsequent generations.

Unraveling the Complex Hybrid Ancestry and Domestication History of Cultivated Strawberry

Cultivated strawberry (Fragaria × ananassa) is one of our youngest domesticates, originating in early eighteenth-century Europe from spontaneous hybrids between wild allo-octoploid species (Fragaria chiloensis and Fragaria virginiana). The improvement of horticultural traits by 300 years of breeding has enabled the global expansion of strawberry production. Here, we describe the genomic history of strawberry domestication from the earliest hybrids to modern cultivars. We observed a significant increase in heterozygosity among interspecific hybrids and a decrease in heterozygosity among domesticated descendants of those hybrids. Selective sweeps were found across the genome in early and modern phases of domestication—59–76% of the selectively swept genes originated in the three less dominant ancestral subgenomes. Contrary to the tenet that genetic diversity is limited in cultivated strawberry, we found that the octoploid species harbor massive allelic diversity and that F. × ananassa harbors as much allelic diversity as either wild founder. We identified 41.8 M subgenome-specific DNA variants among resequenced wild and domesticated individuals. Strikingly, 98% of common alleles and 73% of total alleles were shared between wild and domesticated populations. Moreover, genome-wide estimates of nucleotide diversity were virtually identical in F. chiloensis,F. virginiana, and F. × ananassa (π = 0.0059–0.0060). We found, however, that nucleotide diversity and heterozygosity were significantly lower in modern F. × ananassa populations that have experienced significant genetic gains and have produced numerous agriculturally important cultivars.

Differences between diploid donors are the main contributing factor for subgenome asymmetry measured in either gene ratio or relative diversity in allopolyploids

Subgenome asymmetry (SA) has routinely been attributed to different responses between the subgenomes of a polyploid to various stimuli during evolution. Here, we compared subgenome differences in gene ratio and relative diversity between artificial and natural genotypes of several allopolyploid species. Surprisingly, consistent differences were not detected between these two types of polyploid genotypes, although they differ in times exposed to evolutionary selection. The estimated ratio of shared genes between a subgenome and its diploid donor was invariably higher for the artificial allopolyploid genotypes than those for the natural genotypes, which is expected as it is now well-known that many genes in a species are not shared among all individuals. As the exact diploid parent for a given subgenome is unknown, the estimated ratios of shared genes for the natural genotypes would also include difference among individual genotypes of the diploid donor species. Further, we detected the presence of SA in genotypes before the completion of the polyploidization events as well as in those which were not formed via polyploidization. These results indicate that SA may, to a large degree, reflect differences between its diploid donors or that changes occurred during polyploid evolution are defined by their donor genomes.

From Genetic Maps to QTL Cloning: An Overview for Durum Wheat

Durum wheat is one of the most important cultivated cereal crops, providing nutrients to humans and domestic animals. Durum breeding programs prioritize the improvement of its main agronomic traits; however, the majority of these traits involve complex characteristics with a quantitative inheritance (quantitative trait loci, QTL). This can be solved with the use of genetic maps, new molecular markers, phenotyping data of segregating populations, and increased accessibility to sequences from next-generation sequencing (NGS) technologies. This allows for high-density genetic maps to be developed for localizing candidate loci within a few Kb in a complex genome, such as durum wheat. Here, we review the identified QTL, fine mapping, and cloning of QTL or candidate genes involved in the main traits regarding the quality and biotic and abiotic stresses of durum wheat. The current knowledge on the used molecular markers, sequence data, and how they changed the development of genetic maps and the characterization of QTL is summarized. A deeper understanding of the trait architecture useful in accelerating durum wheat breeding programs is envisioned.

DEFECTIVE ENDOSPERM-D1 (Dee-D1) is crucial for endosperm development in hexaploid wheat

Hexaploid wheat (Triticum aestivum L.) is a natural allopolyploid and provides a usable model system to better understand the genetic mechanisms that underlie allopolyploid speciation through the hybrid genome doubling. Here we aimed to identify the contribution of chromosome 1D in the development and evolution of hexaploid wheat. We identified and mapped a novel DEFECTIVE ENDOSPERM–D1 (Dee-D1) locus on 1DL that is involved in the genetic control of endosperm development. The absence of Dee-D1 leads to non-viable grains in distant crosses and alters grain shape, which negatively affects grain number and thousand-grain weight. Dee-D1 can be classified as speciation locus with a positive effect on the function of genes which are involved in endosperm development in hybrid genomes. The presence of Dee-D1 is necessary for the normal development of endosperm, and thus play an important role in the evolution and improvement of grain yield in hexaploid wheat.

Natalia Tikhenko et al. investigate the genetic contribution of the wheat chromosome 1D to its development and evolution. They find a novel locus, DEFECTIVE ENDOSPERM-D1, on the long arm of 1D that is required for normal endosperm development as its absence leads to non-viable grains and altered grain shape.

Resequencing of 145 Landmark Cultivars Reveals Asymmetric Sub-genome Selection and Strong Founder Genotype Effects on Wheat Breeding in China

Controlled pedigrees and the multi-decade timescale of national crop plant breeding programs offer a unique experimental context for examining how selection affects plant genomes. More than 3000 wheat cultivars have been registered, released, and documented since 1949 in China. In this study, a set of 145 elite cultivars selected from historical points of wheat breeding in China were re-sequenced. A total of 43.75 Tb of sequence data were generated with an average read depth of 17.94× for each cultivar, and more than 60.92 million SNPs and 2.54 million InDels were captured, based on the Chinese Spring RefSeq genome v1.0. Seventy years of breeder-driven selection led to dramatic changes in grain yield and related phenotypes, with distinct genomic regions and phenotypes targeted by different breeders across the decades. There are very clear instances illustrating how introduced Italian and other foreign germplasm was integrated into Chinese wheat programs and reshaped the genomic landscape of local modern cultivars. Importantly, the resequencing data also highlighted significant asymmetric breeding selection among the three sub-genomes: this was evident in both the collinear blocks for homeologous chromosomes and among sets of three homeologous genes. Accumulation of more newly assembled genes in newer cultivars implied the potential value of these genes in breeding. Conserved and extended sharing of linkage disequilibrium (LD) blocks was highlighted among pedigree-related cultivars, in which fewer haplotype differences were detected. Fixation or replacement of haplotypes from founder genotypes after generations of breeding was related to their breeding value. Based on the haplotype frequency changes in LD blocks of pedigree-related cultivars, we propose a strategy for evaluating the breeding value of any given line on the basis of the accumulation (pyramiding) of beneficial haplotypes. Collectively, our study demonstrates the influence of "founder genotypes" on the output of breeding efforts over many decades and also suggests that founder genotype perspectives are in fact more dynamic when applied in the context of modern genomics-informed breeding.

Characterization of the leaf rust responsive ARF genes in wheat (Triticum aestivum L.)

Genome-wide identification, classification, functional characterization and expression analysis of Auxin Responsive Factor (ARF) gene family in wheat reveal their attributes and role during leaf rust infection. Auxins are important plant growth regulators that also impact plant-pathogen interaction. Auxin responsive factors (ARF) are plant specific transcription factors that control responses to auxins. Whole genome investigation of ARF gene family is limited in allohexaploid wheat (Triticum aestivum L.). Comprehensive study of this gene family was carried out by employing the currently available reference genome sequence of wheat. In total, 27 ARF genes were identified and located on the wheat genome as well as were positioned on wheat chromosome arms. Additionally, examination of the predicted genes unveiled a decent degree of relatedness within and among the phylogenetic clades. Leaf rust, caused by the obligate biotrophic fungal pathogen Puccinia triticina, is responsible for drastic loss of wheat crop worldwide reducing grain yield by 10-90%. Expression profiling of ARF genes in retort to leaf rust infection indicated their differential regulation during this plant-pathogen interaction. Highest expression of ARF genes were observed at 12 hpi that was maintained up to 72 hpi during incompatible interaction, whereas the high expression levels receded at 48 hpi during compatible interactions. Few of the identified ARF genes were likely to be post-transcriptionally regulated by microRNAs. Many light and stress responsive elements were detected in the promoter regions of ARF genes. Microsynteny analysis showed the conservation of ARF genes within the members of the Poaceae family. This study provides fundamental details for understanding the different types of ARF genes in wheat and there putative roles during leaf rust-wheat interaction.

Preferential Subgenome Elimination and Chromosomal Structural Changes Occurring in Newly Formed Tetraploid Wheat-Aegilops ventricosa Amphiploid (AABBDvDvNvNv)

Artificial allopolyploids derived from the genera Triticum and Aegilops have been used as genetic resources for wheat improvement and are a classic example of evolution via allopolyploidization. In this study, we investigated chromosomes and subgenome transmission behavior in the newly formed allopolyploid of wheat group via multicolor Fluorescence in situ hybridization (mc-FISH), using pSc119.2, pTa535, and (GAA)₇ as probe combinations, to enabled us to precisely identify individual chromosomes in 381 S₃ and S₄ generations plants derived from reciprocal crosses between Ae. ventricosa (D^vD^vN^vN^v) and T. turgidum (AABB). A higher rate of aneuploidy, constituting 66.04–86.41% individuals, was observed in these two early generations. Of the four constituent subgenomes, D^v showed the highest frequency of elimination, followed by N^v and B, while A was the most stable. In addition, structural chromosomal changes occurred ubiquitously in the selfed progenies of allopolyploids. Among the constituent subgenomes, B showed the highest number of aberrations. In terms of chromosomal dynamics, there was no significant association between the chromosomal behavior model and the cytoplasm, with the exception of chromosomal loss in the D^v subgenome. The chromosome loss frequency in the D^v subgenome was significantly higher in the T. turgidum × Ae. ventricosa cross than in the Ae. ventricosa × T. turgidum cross. This result indicates that, although the D subgenome showed great instability, allopolyploids containing D subgenome could probably be maintained after a certain hybridization in which the D subgenome donor was used as the maternal parent at its onset stage. Our findings provide valuable information pertaining to the behavior patterns of subgenomes during allopolyploidization. Moreover, the allopolyploids developed here could be used as potential resources for the genetic improvement of wheat.

Suitability of Single-Nucleotide Polymorphism Arrays Versus Genotyping-By-Sequencing for Genebank Genomics in Wheat

Genebank genomics promises to unlock valuable diversity for plant breeding but first, one key question is which marker system is most suitable to fingerprint entire genebank collections. Using wheat as model species, we tested for the presence of an ascertainment bias and investigated its impact on estimates of genetic diversity and prediction ability obtained using three marker platforms: simple sequence repeat (SSR), genotyping-by-sequencing (GBS), and array-based SNP markers. We used a panel of 378 winter wheat genotypes including 190 elite lines and 188 plant genetic resources (PGR), which were phenotyped in multi-environmental trials for grain yield and plant height. We observed an ascertainment bias for the array-based SNP markers, which led to an underestimation of the molecular diversity within the population of PGR. In contrast, the marker system played only a minor role for the overall picture of the population structure and precision of genome-wide predictions. Interestingly, we found that rare markers contributed substantially to the prediction ability. This combined with the expectation that valuable novel diversity is most likely rare suggests that markers with minor allele frequency deserve careful consideration in the design of a pre-breeding program.

Genome-Wide Association Study and QTL Meta-Analysis Identified Novel Genomic Loci Controlling Potassium Use Efficiency and Agronomic Traits in Bread Wheat

Potassium use efficiency, a complex trait, directly impacts the yield potential of crop plants. Low potassium efficiency leads to a high use of fertilizers, which is not only farmer unfriendly but also deteriorates the environment. Genome-wide association studies (GWAS) are widely used to dissect complex traits. However, most studies use single-locus one-dimensional GWAS models which do not provide true information about complex traits that are controlled by multiple loci. Here, both single-locus GWAS (MLM) and multi-locus GWAS (pLARmEB, FASTmrMLM, mrMLM, FASTmrEMMA) models were used with genotyping from 90 K Infinium SNP array and phenotype derived from four normal and potassium-stress environments, which identified 534 significant marker-trait associations (MTA) for agronomic and potassium related traits: pLARmEB = 279, FASTmrMLM = 213, mrMLM = 35, MLM = 6, FASTmrEMMA = 1. Further screening of these MTA led to the detection of eleven stable loci: q1A, q1D, q2B-1, q2B-2, q2D, q4D, q5B-1, q5B-2, q5B-3, q6D, and q7A. Moreover, Meta-QTL (MQTL) analysis of four independent QTL studies for potassium deficiency in bread wheat located 16 MQTL on 13 chromosomes. One locus identified in this study (q5B-1) colocalized with an MQTL (MQTL_{_11} ), while the other ten loci were novel associations. Gene ontology of these loci identified 20 putative candidate genes encoding functional proteins involved in key pathways related to stress tolerance, sugar metabolism, and nutrient transport. These findings provide potential targets for breeding potassium stress resistant wheat cultivars and advocate the advantages of multi-locus GWAS models for studying complex traits.

Introgression of a Novel Ug99-Effective Stem Rust Resistance Gene into Wheat and Development of Dasypyrum villosum Chromosome-Specific Markers via Genotyping-by-Sequencing (GBS)

Dasypyrum villosum is a wild relative of common wheat (Triticum aestivum L.) with resistance to Puccinia graminis f. tritici, the causal agent of stem rust, including the highly virulent race TTKSK (Ug99). In order to transfer resistance, T. durum-D. villosum amphiploids were initially developed and used as a bridge to create wheat-D. villosum introgression lines. Conserved ortholog set (COS) markers were used to identify D. villosum chromosome introgression lines, which were then subjected to seedling P. graminis f. tritici resistance screening with race TTKSK. A COS marker-verified line carrying chromosome 2V with TTKSK resistance was further characterized by combined genomic in situ and fluorescent in situ analyses to confirm a monosomic substitution line MS2V(2D) (20″ + 1' 2V[2D]). This is the first report of stem rust resistance on 2V, which was temporarily designated as SrTA10276-2V. To facilitate the use of this gene in wheat improvement, a complete set of previously developed wheat-D. villosum disomic addition lines was subjected to genotyping-by-sequencing analysis to develop D. villosum chromosome-specific markers. On average, 350 markers per chromosome were identified. These markers can be used to develop diagnostic markers for D. villosum-derived genes of interest in wheat improvement.

Linking the International Wheat Genome Sequencing Consortium bread wheat reference genome sequence to wheat genetic and phenomic data

The Wheat@URGI portal has been developed to provide the international community of researchers and breeders with access to the bread wheat reference genome sequence produced by the International Wheat Genome Sequencing Consortium. Genome browsers, BLAST, and InterMine tools have been established for in-depth exploration of the genome sequence together with additional linked datasets including physical maps, sequence variations, gene expression, and genetic and phenomic data from other international collaborative projects already stored in the GnpIS information system. The portal provides enhanced search and browser features that will facilitate the deployment of the latest genomics resources in wheat improvement.

Functional Divergence between Subgenomes and Gene Pairs after Whole Genome Duplications

Gene loss following whole genome duplication (WGD) is often biased, with one subgenome retaining more ancestral genes and the other sustaining more gene deletions. While bias toward the greater expression of gene copies on one subgenome can explain bias in gene loss, this raises the question to what drives differences in gene expression levels between subgenomes. Differences in chromatin modifications and epigenetic markers between subgenomes in several model species are now being identified, providing an explanation for bias in gene expression between subgenomes. WGDs can be classified into duplications with higher, biased gene loss and bias in gene expression between subgenomes versus those with lower, unbiased rates of gene loss and an absence of detectable bias between subgenomes; however, the originally proposed link between these two classes and whether WGD results from an allo- or autopolyploid event is inconsistent with recent data from the allopolyploid Capsella bursa-pastoris. The gene balance hypothesis can explain bias in the functional categories of genes retained following WGD, the difference in gene loss rates between unbiased and biased WGDs, and how plant genomes have avoided being overrun with genes encoding dose-sensitive subunits of multiprotein complexes. Comparisons of gene expression patterns between retained transcription factor pairs in maize suggest the high degree of retention for WGD-derived pairs of transcription factors may instead be explained by the older duplication-degeneration-complementation model.

Decomposing Additive Genetic Variance Revealed Novel Insights into Trait Evolution in Synthetic Hexaploid Wheat

Whole genome duplication (WGD) is an evolutionary phenomenon, which causes significant changes to genomic structure and trait architecture. In recent years, a number of studies decomposed the additive genetic variance explained by different sets of variants. However, they investigated diploid populations only and none of the studies examined any polyploid organism. In this research, we extended the application of this approach to polyploids, to differentiate the additive variance explained by the three subgenomes and seven sets of homoeologous chromosomes in synthetic allohexaploid wheat (SHW) to gain a better understanding of trait evolution after WGD. Our SHW population was generated by crossing improved durum parents (Triticum turgidum; 2n = 4x = 28, AABB subgenomes) with the progenitor species Aegilops tauschii (syn Ae. squarrosa, T. tauschii; 2n = 2x = 14, DD subgenome). The population was phenotyped for 10 fungal/nematode resistance traits as well as two abiotic stresses. We showed that the wild D subgenome dominated the additive effect and this dominance affected the A more than the B subgenome. We provide evidence that this dominance was not inflated by population structure, relatedness among individuals or by longer linkage disequilibrium blocks observed in the D subgenome within the population used for this study. The cumulative size of the three homoeologs of the seven chromosomal groups showed a weak but significant positive correlation with their cumulative explained additive variance. Furthermore, an average of 69% for each chromosomal group's cumulative additive variance came from one homoeolog that had the highest explained variance within the group across all 12 traits. We hypothesize that structural and functional changes during diploidization may explain chromosomal group relations as allopolyploids keep balanced dosage for many genes. Our results contribute to a better understanding of trait evolution mechanisms in polyploidy, which will facilitate the effective utilization of wheat wild relatives in breeding.

Characterization and Expression Patterns of Auxin Response Factors in Wheat

Auxin response factors (ARFs) are important transcription factors involved in both the auxin signaling pathway and the regulatory development of various plant organs. In this study, 23 TaARF members encoded by a total of 68 homeoalleles were isolated from 18 wheat chromosomes (excluding chromosome 4). The TaARFs, including their conserved domains, exon/intron structures, related microRNAs, and alternative splicing (AS) variants, were then characterized. Phylogenetic analysis revealed that members of the TaARF family share close homology with ARFs in other grass species. qRT-PCR analyses revealed that 20 TaARF members were expressed in different organs and tissues and that the expression of some members significantly differed in the roots, stems, and leaves of wheat seedlings in response to exogenous auxin treatment. Moreover, protein network analyses and co-expression results showed that TaTIR1-TaARF15/18/19-TaIAA13 may interact at both the protein and genetic levels. The results of subsequent evolutionary analyses showed that three transcripts of TaARF15 in the A subgenome of wheat exhibited high evolutionary rate and underwent positive selection. Transgenic analyses indicated that TaARF15-A.1 promoted the growth of roots and leaves of Arabidopsis thaliana and was upregulated in the overexpression plants after auxin treatment. Our results will provide reference information for subsequent research and utilization of the TaARF gene family.

Evolution of the Aux/IAA Gene Family in Hexaploid Wheat

The Aux/IAA (IAA) gene family, involved in the auxin signalling pathway, acts as an important regulator in plant growth and development. In this study, we explored the evolutionary trajectory of the IAA family in common wheat. The results showed ten pairs of paralogs among 34 TaIAA family members. Seven of the pairs might have undergone segmental duplication, and the other three pairs appear to have experienced tandem duplication. Except for TaIAA15-16, these duplication events occurred in the ancestral genomes before the divergence of Triticeae. After that point, two polyploidization events shaped the current TaIAA family consisting of three subgenomic copies. The structure or expression pattern of the TaIAA family begins to differentiate in the hexaploid genome, where TaIAAs in the D genome lost more genes (eight) and protein secondary structures (α1, α3 and β5) than did the other two genomes. Expression analysis showed that six members of the TaIAA family were not expressed, and members such as TaIAA8, 15, 16, 28 and 33 exhibited tissue-specific expression patterns. In addition, three of the ten pairs of paralogs (TaIAA5-12, TaIAA15-16 and TaIAA29-30) showed similar expression patterns, and another five paralog pairs displayed differential expression patterns. Phylogenetic analysis showed that paralog pairs with high rates of evolution (ω > ω ₀), particularly TaIAA15-16 and TaIAA29-30, experienced greater motif loss, with only zero to two interacting IAA proteins. In contrast, most paralogous genes with low ω, such as TaIAA5-12, had more complete motifs and higher degrees of interaction with other family members.

Genotyping-by-Sequencing (GBS) Revealed Molecular Genetic Diversity of Iranian Wheat Landraces and Cultivars

Background: Genetic diversity is an essential resource for breeders to improve new cultivars with desirable characteristics. Recently, genotyping-by-sequencing (GBS), a next-generation sequencing (NGS) technology that can simplify complex genomes, has now be used as a high-throughput and cost-effective molecular tool for routine breeding and screening in many crop species, including the species with a large genome. Results: We genotyped a diversity panel of 369 Iranian hexaploid wheat accessions including 270 landraces collected between 1931 and 1968 in different climate zones and 99 cultivars released between 1942 to 2014 using 16,506 GBS-based single nucleotide polymorphism (GBS-SNP) markers. The B genome had the highest number of mapped SNPs while the D genome had the lowest on both the Chinese Spring and W7984 references. Structure and cluster analyses divided the panel into three groups with two landrace groups and one cultivar group, suggesting a high differentiation between landraces and cultivars and between landraces. The cultivar group can be further divided into four subgroups with one subgroup was mostly derived from Iranian ancestor(s). Similarly, landrace groups can be further divided based on years of collection and climate zones where the accessions were collected. Molecular analysis of variance indicated that the genetic variation was larger between groups than within group. Conclusion: Obvious genetic diversity in Iranian wheat was revealed by analysis of GBS-SNPs and thus breeders can select genetically distant parents for crossing in breeding. The diverse Iranian landraces provide rich genetic sources of tolerance to biotic and abiotic stresses, and they can be useful resources for the improvement of wheat production in Iran and other countries.

Revisiting Pivotal-Differential Genome Evolution in Wheat

An interesting and possibly unique pattern of genome evolution following polyploidy can be observed among allopolyploids of the Triticum and Aegilops genera (wheat group). Most polyploids in this group are presumed to share a common unaltered (pivotal) subgenome (U, D, or A) together with one or two modified (differential) subgenomes, a status that has been referred to as 'pivotal-differential' genome evolution. In this review we discuss various mechanisms that could be responsible for this evolutionary pattern, as well as evidence for and against the putative evolutionary mechanisms involved. We suggest that, in light of recent advances in genome sequencing and related technologies in the wheat group, the time has come to reopen the investigation into pivotal-differential genome evolution.

Abiotic stress miRNomes in the Triticeae

The continued growth in world population necessitates increases in both the quantity and quality of agricultural production. Triticeae members, particularly wheat and barley, make an important contribution to world food reserves by providing rich sources of carbohydrate and protein. These crops are grown over diverse production environments that are characterized by a range of environmental or abiotic stresses. Abiotic stresses such as drought, heat, salinity, or nutrient deficiencies and toxicities cause large yield losses resulting in economic and environmental damage. The negative effects of abiotic stresses have increased at an alarming rate in recent years and are predicted to further deteriorate due to climate change, land degradation, and declining water supply. New technologies have provided an important tool with great potential for improving crop tolerance to the abiotic stresses: microRNAs (miRNAs). miRNAs are small regulators of gene expression that act on many different molecular and biochemical processes such as development, environmental adaptation, and stress tolerance. miRNAs can act at both the transcriptional and post-transcriptional levels, although post-transcriptional regulation is the most common in plants where miRNAs can inhibit the translation of their mRNA targets via complementary binding and cleavage. To date, expression of several miRNA families such as miR156, miR159, and miR398 has been detected as responsive to environmental conditions to regulate stress-associated molecular mechanisms individually and/or together with their various miRNA partners. Manipulation of these miRNAs and their targets may pave the way to improve crop performance under several abiotic stresses. Here, we summarize the current status of our knowledge on abiotic stress-associated miRNAs in members of the Triticeae tribe, specifically in wheat and barley, and the miRNA-based regulatory mechanisms triggered by stress conditions. Exploration of further miRNA families together with their functions under stress will improve our knowledge and provide opportunities to enhance plant performance to help us meet global food demand.

Wheat paleohistory created asymmetrical genomic evolution

Following the triplication reported in Brassiceae ∼10million years ago, and at the basis of rosids ∼100million years ago, bias in organization and regulation, known as subgenome dominance, has been reported between the three post-polyploidy compartments referenced to as less fractionated (LF), medium fractionated (MF1) and more fractionated (MF2), that have been proposed to derive from an hexaploidization event involving ancestors of 7-14-21 chromosomes. Modern bread wheat experienced similar paleohistory during the last half million year of evolution opening a new hypothesis where the wheat genome is at the earliest stages on the road of diploidization through subgenome dominance driving asymmetry in gene content, gene expression abundance, transposable element content as dynamics and epigenetic control between the A, B and D subgenomes.

Genome sequencing supports a multi-vertex model for Brassiceae species

The economically important Brassica genus is a good system for studying the evolution of polyploids. Brassica genomes have undergone whole genome triplication (WGT). Subgenome dominance phenomena such as biased gene fractionation and dominant gene expression were observed in tripled genomes of Brassica. The genome of radish (Raphanus sativus), another important crop of tribe Brassiceae, was derived from the same WGT event and shows similar subgenome dominance. These findings and molecular dating indicate that radish occupies a similar evolutionary origin as that of Brassica species. Here, we extended the Brassica "triangle of U" to a multi-vertex model. This model describes the relationships or the potential of using more Brassiceae mesohexaploids in the creation of new allotetraploid oil or vegetable crop species.

Reconciling the evolutionary origin of bread wheat (Triticum aestivum)

The origin of bread wheat (Triticum aestivum; AABBDD) has been a subject of controversy and of intense debate in the scientific community over the last few decades. In 2015, three articles published in New Phytologist discussed the origin of hexaploid bread wheat (AABBDD) from the diploid progenitors Triticum urartu (AA), a relative of Aegilops speltoides (BB) and Triticum tauschii (DD). Access to new genomic resources since 2013 has offered the opportunity to gain novel insights into the paleohistory of modern bread wheat, allowing characterization of its origin from its diploid progenitors at unprecedented resolution. We propose a reconciled evolutionary scenario for the modern bread wheat genome based on the complementary investigation of transposable element and mutation dynamics between diploid, tetraploid and hexaploid wheat. In this scenario, the structural asymmetry observed between the A, B and D subgenomes in hexaploid bread wheat derives from the cumulative effect of diploid progenitor divergence, the hybrid origin of the D subgenome, and subgenome partitioning following the polyploidization events.

Combined Genomic and Genetic Data Integration of Major Agronomical Traits in Bread Wheat (Triticum aestivum L.)

The high resolution integration of bread wheat genetic and genomic resources accumulated during the last decades offers the opportunity to unveil candidate genes driving major agronomical traits to an unprecedented scale. We combined 27 public quantitative genetic studies and four genetic maps to deliver an exhaustive consensus map consisting of 140,315 molecular markers hosting 221, 73, and 82 Quantitative Trait Loci (QTL) for respectively yield, baking quality, and grain protein content (GPC) related traits. Projection of the consensus genetic map and associated QTLs onto the wheat syntenome made of 99,386 genes ordered on the 21 chromosomes delivered a complete and non-redundant repertoire of 18, 8, 6 metaQTLs for respectively yield, baking quality and GPC, altogether associated to 15,772 genes (delivering 28,630 SNP-based makers) including 37 major candidates. Overall, this study illustrates a translational research approach in transferring information gained from grass relatives to dissect the genomic regions hosting major loci governing key agronomical traits in bread wheat, their flanking markers and associated candidate genes to be now considered as a key resource for breeding programs.

Epigenetic regulation of subgenome dominance following whole genome triplication in Brassica rapa

Subgenome dominance is an important phenomenon observed in allopolyploids after whole genome duplication, in which one subgenome retains more genes as well as contributes more to the higher expressing gene copy of paralogous genes. To dissect the mechanism of subgenome dominance, we systematically investigated the relationships of gene expression, transposable element (TE) distribution and small RNA targeting, relating to the multicopy paralogous genes generated from whole genome triplication in Brassica rapa. The subgenome dominance was found to be regulated by a relatively stable factor established previously, then inherited by and shared among B. rapa varieties. In addition, we found a biased distribution of TEs between flanking regions of paralogous genes. Furthermore, the 24-nt small RNAs target TEs and are negatively correlated to the dominant expression of individual paralogous gene pairs. The biased distribution of TEs among subgenomes and the targeting of 24-nt small RNAs together produce the dominant expression phenomenon at a subgenome scale. Based on these findings, we propose a bucket hypothesis to illustrate subgenome dominance and hybrid vigor. Our findings and hypothesis are valuable for the evolutionary study of polyploids, and may shed light on studies of hybrid vigor, which is common to most species.

Deciphering the evolutionary interplay between subgenomes following polyploidy: A paleogenomics approach in grasses

How did plant species emerge from their most recent common ancestors (MRCAs) 250 million years ago? Modern plant genomes help to address such key questions in unveiling precise species genealogies. The field of paleogenomics is undergoing a paradigm shift for investigating species evolution from the study of ancestral genomes from extinct species to deciphering the evolutionary forces (in terms of duplication, fusion, fission, deletion, and translocation) that drove present-day plant diversity (in terms of chromosome/gene number and genome size). In this review, inferred ancestral karyotype genomes are shown to be powerful tools to (1) unravel the past history of extant species by recovering the variations of ancestral genomic compartments and (2) accelerate translational research by facilitating the transfer of genomic information from model systems to species of agronomic interest.

Transcriptome shock invokes disruption of parental expression-conserved genes in tetraploid wheat

Allopolyploidy often triggers phenotypic novelty and gene expression remolding in the resulting polyploids. In this study, we employed multiple phenotypic and genetic approaches to investigate the nature and consequences of allotetraploidization between A- and S-subgenome of tetraploid wheat. Results showed that karyotype of the nascent allopolyploid plants (AT2) is stable but they showed clear novelty in multiple morphological traits which might have positively contributed to the initial establishment of the tetraploids. Further microarray-based transcriptome profiling and gene-specific cDNA-pyrosequencing have documented that transcriptome shock was exceptionally strong in AT2, but a substantial proportion of the induced expression changes was rapidly stabilized in early generations. Meanwhile, both additive and nonadditive expression genes showed extensive homeolog expression remodeling and which have led to the subgenome expression dominance in leaf and young inflorescence of AT2. Through comparing the homeolog-expressing patterns between synthetic and natural tetraploid wheats, it appears that the shock-induced expression changes at both the total expression level and subgenome homeolog partitioning are evolutionarily persistent. Together, our study shed new light on how gene expression changes have rapidly occurred at the initial stage following allotetraploidization, as well as their evolutionary relevance, which may have implications for wheat improvements.

DensityMap: a genome viewer for illustrating the densities of features

Background

Several tools are available for visualizing genomic data. Some, such as Gbrowse and Jbrowse, are very efficient for small genomic regions, but they are not suitable for entire genomes. Others, like Phenogram and CViT, can be used to visualise whole genomes, but are not designed to display very dense genomic features (eg: interspersed repeats). We have therefore developed DensityMap, a lightweight Perl program that can display the densities of several features (genes, ncRNA, cpg, etc.) along chromosomes on the scale of the whole genome. A critical advantage of DensityMap is that it uses GFF annotation files directly to compute the densities of features without needing additional information from the user. The resulting picture is readily configurable, and the colour scales used can be customized for a best fit to the data plotted.

Results

DensityMap runs on Linux architecture with few requirements so that users can easily and quickly visualize the distributions and densities of genomic features for an entire genome. The input is GFF3-formated data representing chromosomes (linkage groups or pseudomolecules) and sets of features which are used to calculate representations in density maps. In practise, DensityMap uses a tilling window to compute the density of one or more features and the number of bases covered by these features along chromosomes. The densities are represented by colour scales that can be customized to highlight critical points. DensityMap can compare the distributions of features; it calculates several chromosomal density maps in a single image, each of which describes a different genomic feature. It can also use the genome nucleotide sequence to compute and plot a density map of the GC content along chromosomes.

Conclusions

DensityMap is a compact, easily-used tool for displaying the distribution and density of all types of genomic features within a genome. It is flexible enough to visualize the densities of several types of features in a single representation. The images produced are readily configurable and their SVG format ensures that they can be edited.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-1055-0) contains supplementary material, which is available to authorized users.

Wheat miRNA ancestors: evident by transcriptome analysis of A, B, and D genome donors

MicroRNAs are critical players of post-transcriptional gene regulation with profound effects on the fundamental processes of cellular life. Their identification and characterization, together with their targets, hold great significance in exploring and exploiting their roles on a functional context, providing valuable clues into the regulation of important biological processes, such as stress tolerance or environmental adaptation. Wheat is a hardy crop, extensively harvested in temperate regions, and is a major component of the human diet. With the advent of the next generation sequencing technologies considerably decreasing sequencing costs per base-pair, genomic, and transcriptomic data from several wheat species, including the progenitors and wild relatives have become available. In this study, we performed in silico identification and comparative analysis of microRNA repertoires of bread wheat (Triticum aestivum L.) and its diploid progenitors and relatives, Aegilops sharonensis, Aegilops speltoides, Aegilops tauschii, Triticum monococcum, and Triticum urartu through the utilization of publicly available transcriptomic data. Over 200 miRNA families were identified, majority of which have not previously been reported. Ancestral relationships expanded our understanding of wheat miRNA evolution, while T. monococcum miRNAs delivered important clues on the effects of domestication on miRNA expression. Comparative analyses on wild Ae. sharonensis accessions highlighted candidate miRNAs that can be linked to stress tolerance. The miRNA repertoires of bread wheat and its diploid progenitors and relatives provide important insight into the diversification and distribution of miRNA genes, which should contribute to the elucidation of miRNA evolution of Poaceae family. A thorough understanding of the convergent and divergent expression profiles of miRNAs in different genetic backgrounds can provide unique opportunities to modulation of gene regulation for better crop performance.

Transcriptome asymmetry in synthetic and natural allotetraploid wheats, revealed by RNA-sequencing

Allopolyploidization has occurred frequently within the Triticum-Aegilops complex which provides a suitable system to investigate how allopolyploidization shapes the expression patterns of duplicated homeologs. We have conducted transcriptome-profiling of leaves and young inflorescences in wild and domesticated tetraploid wheats, Triticum turgidum ssp. dicoccoides (BBAA) and ssp. durum (BBAA), an extracted tetraploid (BBAA), and a synthetic tetraploid (S(l) S(l) AA) wheat together with its diploid parents, Aegilops longissima (S(l) S(l) ) and Triticum urartu (AA). The two diploid species showed tissue-specific differences in genome-wide ortholog expression, which plays an important role in transcriptome shock-mediated homeolog expression rewiring and hence transcriptome asymmetry in the synthetic tetraploid. Further changes of homeolog expression apparently occurred in natural tetraploid wheats, which led to novel transcriptome asymmetry between the two subgenomes. In particular, our results showed that extremely biased homeolog expression can occur rapidly upon the allotetraploidzation and this trend is further enhanced in the course of domestication and evolution of polyploid wheats. Our results suggest that allopolyploidization is accompanied by distinct phases of homeolog expression changes, with parental legacy playing major roles in the immediate rewiring of homeolog expression upon allopolyploidization, while evolution and domestication under allotetraploidy drive further homeolog-expression changes toward re-established subgenome expression asymmetry.

Heritable alteration of DNA methylation induced by whole-chromosome aneuploidy in wheat

Aneuploidy causes changes in gene expression and phenotypes in all organisms studied. A previous study in the model plant Arabidopsis thaliana showed that aneuploidy-generated phenotypic changes can be inherited to euploid progenies and implicated an epigenetic underpinning of the heritable variations. Based on an analysis by amplified fragment length polymorphism and methylation-sensitive amplified fragment length polymorphism markers, we found that although genetic changes at the nucleotide sequence level were negligible, extensive changes in cytosine DNA methylation patterns occurred in all studied homeologous group 1 whole-chromosome aneuploid lines of common wheat (Triticum aestivum), with monosomic 1A showing the greatest amount of methylation changes. The changed methylation patterns were inherited by euploid progenies derived from the aneuploid parents. The aneuploidy-induced DNA methylation alterations and their heritability were verified at selected loci by bisulfite sequencing. Our data have provided empirical evidence supporting earlier suggestions that heritability of aneuploidy-generated, but aneuploidy-independent, phenotypic variations may have an epigenetic basis. That at least one type of aneuploidy - monosomic 1A - was able to cause significant epigenetic divergence of the aneuploid plants and their euploid progenies also lends support to recent suggestions that aneuploidy may have played an important and protracted role in polyploid genome evolution.

Genome wide identification of C1-2i zinc finger proteins and their response to abiotic stress in hexaploid wheat

The C1-2i wheat Q-type C2H2 zinc finger protein (ZFP) transcription factor subclass has been reported to play important roles in plant stress responses. This subclass of ZFPs has not been studied in hexaploid wheat (Triticum aestivum) and we aimed to identify all members of this subclass and evaluate their responses to different abiotic stresses causing oxidative stress. Exploiting the recently published wheat draft genome sequence, we identified 53 members (including homoeologs from A, B and D genomes) of the C1-2i wheat Q-type C2H2 ZFPs (TaZFPs) representing 21 genes. Evolution analysis revealed that 9 TaZFPs members are directly inherited from the parents Triticum urartu and Aegilops tauschii, while 15 diverged through neoploidization events. This TaZFP subclass is responsive to the oxidative stress generator H2O2 and to high light, drought stress and flooding. Most TaZFPs are responsive to H2O2 (37/53), high light (44/53), flooding (31/53) or drought (37/53); 32 TaZFPs were up-regulated by at least 3 stresses and 16 were responsive to all stresses tested. A large number of these TaZFPs were physically mapped on different wheat draft genome sequences with known markers useful for QTL mapping. Our results show that the C1-2i subclass of TaZFPs is associated with responses to different abiotic stresses and that most TaZFPs (30/53 or 57 %) are located on group 5 chromosomes known to be involved in environment adaptation. Detailed characterization of these novel wheat TaZFPs and their association to QTL or eQTL may help to design wheat cultivars with improved tolerance to abiotic stress.

Comparative genomic analysis of duplicated homoeologous regions involved in the resistance of Brassica napus to stem canker

All crop species are current or ancient polyploids. Following whole genome duplication, structural and functional modifications result in differential gene content or regulation in the duplicated regions, which can play a fundamental role in the diversification of genes underlying complex traits. We have investigated this issue in Brassica napus, a species with a highly duplicated genome, with the aim of studying the structural and functional organization of duplicated regions involved in quantitative resistance to stem canker, a disease caused by the fungal pathogen Leptosphaeria maculans. Genome-wide association analysis on two oilseed rape panels confirmed that duplicated regions of ancestral blocks E, J, R, U, and W were involved in resistance to stem canker. The structural analysis of the duplicated genomic regions showed a higher gene density on the A genome than on the C genome and a better collinearity between homoeologous regions than paralogous regions, as overall in the whole B. napus genome. The three ancestral sub-genomes were involved in the resistance to stem canker and the fractionation profile of the duplicated regions corresponded to what was expected from results on the B. napus progenitors. About 60% of the genes identified in these duplicated regions were single-copy genes while less than 5% were retained in all the duplicated copies of a given ancestral block. Genes retained in several copies were mainly involved in response to stress, signaling, or transcription regulation. Genes with resistance-associated markers were mainly retained in more than two copies. These results suggested that some genes underlying quantitative resistance to stem canker might be duplicated genes. Genes with a hydrolase activity that were retained in one copy or R-like genes might also account for resistance in some regions. Further analyses need to be conducted to indicate to what extent duplicated genes contribute to the expression of the resistance phenotype.

A high-density, SNP-based consensus map of tetraploid wheat as a bridge to integrate durum and bread wheat genomics and breeding

Consensus linkage maps are important tools in crop genomics. We have assembled a high-density tetraploid wheat consensus map by integrating 13 data sets from independent biparental populations involving durum wheat cultivars (Triticum turgidum ssp. durum), cultivated emmer (T. turgidum ssp. dicoccum) and their ancestor (wild emmer, T. turgidum ssp. dicoccoides). The consensus map harboured 30 144 markers (including 26 626 SNPs and 791 SSRs) half of which were present in at least two component maps. The final map spanned 2631 cM of all 14 durum wheat chromosomes and, differently from the individual component maps, all markers fell within the 14 linkage groups. Marker density per genetic distance unit peaked at centromeric regions, likely due to a combination of low recombination rate in the centromeric regions and even gene distribution along the chromosomes. Comparisons with bread wheat indicated fewer regions with recombination suppression, making this consensus map valuable for mapping in the A and B genomes of both durum and bread wheat. Sequence similarity analysis allowed us to relate mapped gene-derived SNPs to chromosome-specific transcripts. Dense patterns of homeologous relationships have been established between the A- and B-genome maps and between nonsyntenic homeologous chromosome regions as well, the latter tracing to ancient translocation events. The gene-based homeologous relationships are valuable to infer the map location of homeologs of target loci/QTLs. Because most SNP and SSR markers were previously mapped in bread wheat, this consensus map will facilitate a more effective integration and exploitation of genes and QTL for wheat breeding purposes.

Making the Bread: Insights from Newly Synthesized Allohexaploid Wheat

Bread wheat (or common wheat, Triticum aestivum) is an allohexaploid (AABBDD, 2n = 6x = 42) that arose by hybridization between a cultivated tetraploid wheat T. turgidum (AABB, 2n = 4x = 28) and the wild goatgrass Aegilops tauschii (DD, 2n = 2x = 14). Polyploidization provided niches for rigorous genome modification at cytogenetic, genetic, and epigenetic levels, rendering a broader spread than its progenitors. This review summarizes the latest advances in understanding gene regulation mechanisms in newly synthesized allohexaploid wheat and possible correlation with polyploid growth vigor and adaptation. Cytogenetic studies reveal persistent association of whole-chromosome aneuploidy with nascent allopolyploids, in contrast to the genetic stability in common wheat. Transcriptome analysis of the euploid wheat shows that small RNAs are driving forces for homoeo-allele expression regulation via genetic and epigenetic mechanisms. The ensuing non-additively expressed genes and those with expression level dominance to the respective progenitor may play distinct functions in growth vigor and adaptation in nascent allohexaploid wheat. Further genetic diploidization of allohexaploid wheat is not random. Regional asymmetrical gene distribution, rather than subgenome dominance, is observed in both synthetic and natural allohexaploid wheats. The combinatorial effects of diverged genomes, subsequent selection of specific gene categories, and subgenome-specific traits are essential for the successful establishment of common wheat.

Decoding the oak genome: public release of sequence data, assembly, annotation and publication strategies

The 1.5 Gbp/2C genome of pedunculate oak (Quercus robur) has been sequenced. A strategy was established for dealing with the challenges imposed by the sequencing of such a large, complex and highly heterozygous genome by a whole-genome shotgun (WGS) approach, without the use of costly and time-consuming methods, such as fosmid or BAC clone-based hierarchical sequencing methods. The sequencing strategy combined short and long reads. Over 49 million reads provided by Roche 454 GS-FLX technology were assembled into contigs and combined with shorter Illumina sequence reads from paired-end and mate-pair libraries of different insert sizes, to build scaffolds. Errors were corrected and gaps filled with Illumina paired-end reads and contaminants detected, resulting in a total of 17,910 scaffolds (>2 kb) corresponding to 1.34 Gb. Fifty per cent of the assembly was accounted for by 1468 scaffolds (N50 of 260 kb). Initial comparison with the phylogenetically related Prunus persica gene model indicated that genes for 84.6% of the proteins present in peach (mean protein coverage of 90.5%) were present in our assembly. The second and third steps in this project are genome annotation and the assignment of scaffolds to the oak genetic linkage map. In accordance with the Bermuda and Fort Lauderdale agreements and the more recent Toronto Statement, the oak genome data have been released into public sequence repositories in advance of publication. In this presubmission paper, the oak genome consortium describes its principal lines of work and future directions for analyses of the nature, function and evolution of the oak genome.