Genetic Mapping of Flavonoid Grain Pigments in Durum Wheat.

Pigmented cereal grains with high levels of flavonoid compounds have attracted the attention of nutritional science backing the development of functional foods with claimed health benefits. In this study, we report results on the genetic factors controlling grain pigmentation in durum wheat using a segregant population of recombinant inbred lines (RILs) derived from a cross between an Ethiopian purple grain accession and an Italian amber grain cultivar. The RIL population was genotyped by the wheat 25K SNP array and phenotyped for total anthocyanin content (TAC), grain color, and the L*, a*, and b* color index of wholemeal flour, based on four field trials. The mapping population showed a wide variation for the five traits in the different environments, a significant genotype x environment interaction, and high heritability. A total of 5942 SNP markers were used for constructing the genetic linkage map, with an SNP density ranging from 1.4 to 2.9 markers/cM. Two quantitative trait loci (QTL) were identified for TAC mapping on chromosome arms 2AL and 7BS in the same genomic regions of two detected QTL for purple grain. The interaction between the two QTL was indicative of an inheritance pattern of two loci having complementary effects. Moreover, two QTL for red grain color were detected on chromosome arms 3AL and 3BL. The projection of the four QTL genomic regions on the durum wheat Svevo reference genome disclosed the occurrence of the candidate genes Pp-A3, Pp-B1, R-A1, and R-B1 involved in flavonoid biosynthetic pathways and encoding of transcription factors bHLH (Myc-1) and MYB (Mpc1, Myb10), previously reported in common wheat. The present study provides a set of molecular markers associated with grain pigments useful for the selection of essential alleles for flavonoid synthesis in durum wheat breeding programs and enhancement of the health-promoting quality of derived foods.

A Survey of the Transcriptomic Resources in Durum Wheat: Stress Responses, Data Integration and Exploitation.

Durum wheat (Triticum turgidum subsp. durum (Desf.) Husn.) is an allotetraploid cereal crop of worldwide importance, given its use for making pasta, couscous, and bulgur. Under climate change scenarios, abiotic (e.g., high and low temperatures, salinity, drought) and biotic (mainly exemplified by fungal pathogens) stresses represent a significant limit for durum cultivation because they can severely affect yield and grain quality. The advent of next-generation sequencing technologies has brought a huge development in transcriptomic resources with many relevant datasets now available for durum wheat, at various anatomical levels, also focusing on phenological phases and environmental conditions. In this review, we cover all the transcriptomic resources generated on durum wheat to date and focus on the corresponding scientific insights gained into abiotic and biotic stress responses. We describe relevant databases, tools and approaches, including connections with other "omics" that could assist data integration for candidate gene discovery for bio-agronomical traits. The biological knowledge summarized here will ultimately help in accelerating durum wheat breeding.

Whole-exome sequencing of selected bread wheat recombinant inbred lines as a useful resource for allele mining and bulked segregant analysis.

Although wheat (Triticum aestivum L.) is the main staple crop in the world and a major source of carbohydrates and proteins, functional genomics and allele mining are still big challenges. Given the advances in next-generation sequencing (NGS) technologies, the identification of causal variants associated with a target phenotype has become feasible. For these reasons, here, by combining sequence capture and target-enrichment methods with high-throughput NGS re-sequencing, we were able to scan at exome-wide level 46 randomly selected bread wheat individuals from a recombinant inbred line population and to identify and classify a large number of single nucleotide polymorphisms (SNPs). For technical validation of results, eight randomly selected SNPs were converted into Kompetitive Allele-Specific PCR (KASP) markers. This resource was established as an accessible and reusable molecular toolkit for allele data mining. The dataset we are making available could be exploited for novel studies on bread wheat genetics and as a foundation for starting breeding programs aimed at improving different key agronomic traits.

Polyphenolic Compound Variation in Globe Artichoke Cultivars as Affected by Fertilization and Biostimulants Application.

Globe artichoke is an ancient herbaceous plant native to the Mediterranean Basin. The edible part of the plant is limited to the fleshy leaves (bracts) and receptacle of a large immature inflorescence (head) that has been shown to be a rich source of bioactive compounds. Nutritional and pharmacological properties of artichoke heads and leaves are attributed mainly to polyphenolic compounds and inulin present at high concentration. In this study, polyphenols were investigated in two artichoke cultivars (Opal and Madrigal) in response to four nitrogen rates and foliar applications of biostimulating products under drip irrigation. Field experiments were carried out over two growing seasons (2015-2016, 2016-2017) in Policoro (MT), Southern Italy, on a deep clay soil in sub-humid climate conditions. Phenolic compounds were isolated and characterized by means of high-performance liquid chromatography with photodiode array detection and electrospray ionization/mass spectrometry (HPLC-DAD-MS/MS) analysis. In both cultivars, caffeoylquinic acids were more abundant when a dose of 100 kg ha-1 of ammonium nitrate was provided, whereas apigenins were not affected by nitrogen fertilization. Luteolins increased in cv Opal and decreased in cv Madrigal following N fertilization. The application of biostimulants (3 L ha-1) favored the accumulation of polyphenols, in particular of caffeoylquinic acids and apigenin, in artichoke heads in both cultivars. The results obtained highlight some positive aspects related to the synergistic effect of nitrogen fertilization and biostimulant foliar application.

R2R3-MYBs in Durum Wheat: Genome-Wide Identification, Poaceae-Specific Clusters, Expression, and Regulatory Dynamics Under Abiotic Stresses.

MYB transcription factors (TFs) represent one of the biggest TF families in plants, being involved in various specific plant processes, such as responses to biotic and abiotic stresses. The implication of MYB TFs in the tolerance mechanisms to abiotic stress is particularly interesting for crop breeding, since environmental conditions can negatively affect growth and productivity. Wheat is a worldwide-cultivated cereal, and is a major source of plant-based proteins in human food. In particular, durum wheat plays an important role in global food security improvement, since its adaptation to hot and dry conditions constitutes the base for the success of wheat breeding programs in future. In the present study, a genome-wide identification of R2R3-MYB TFs in durum wheat was performed. MYB profile search and phylogenetic analyses based on homology with Arabidopsis and rice MYB TFs led to the identification of 233 R2R3-TdMYB (Triticum durum MYB). Three Poaceae-specific MYB clusters were detected, one of which had never been described before. The expression of eight selected genes under different abiotic stress conditions, revealed that most of them responded especially to salt and drought stress. Finally, gene regulatory network analyses led to the identification of 41 gene targets for three TdR2R3-MYBs that represent novel candidates for functional analyses. This study provides a detailed description of durum wheat R2R3-MYB genes and contributes to a deeper understanding of the molecular response of durum wheat to unfavorable climate conditions.

Sweet Cherry Diversity and Relationships in Modern and Local Varieties Based on SNP Markers.

The sweet cherry is an important fruit species that is widespread globally. In addition to the well-known traditional and modern varieties, a myriad of landraces is present in Europe, as well as in southern Italy. This study aims to evaluate the population structure, genetic relationships, and cases of duplicate samples in a collection of 143 accessions using GBS-derived SNP markers. The genetic material under investigation includes modern commercial varieties, ancient European and American varieties, landraces, and individuals retrieved from small orchards. Some of the known varieties were genetically analyzed here for the first time. In addition, several genotypes were collected from the Basilicata region (southern Italy), an area largely unexplored for sweet cherry genetic resources. The relationships among genotypes were assessed using four different methods: allele frequency and ancestry estimation, principal component analysis, Neighbor-Joining tree, and identity-by-state estimation. The analyses returned quite congruent results and highlighted the presence of four main genetic groups, namely: (i) American varieties, (ii) the 'Germersdorfer-Ferrovia' cluster, (iii) the 'Burlat' group, and (iv) the group of Italian landraces. The main drivers of clustering were ancestry, geographical distribution, and some important traits such as self-compatibility. The sweet cherries from Basilicata, herewith examined for the first time, were mostly distributed within the group of Italian landraces, being particularly linked to the autochthonous varieties of the Campania region. However, some genotypes were outside this group, thus suggesting the introduction of genetic material from other Italian regions or from European countries. The considerable amount of American and European modern varieties analyzed are genetically very closely related, suggesting a reduced genetic basis. In addition, we highlighted the discriminating ability of SNP markers to distinguish between an original variety and its mutant. Overall, our results may be useful in defining conservation strategies for sweet cherry germplasm and developing future breeding programs to enlarge the genetic basis of commercial varieties.

Genotyping-by-Sequencing in Vigna unguiculata Landraces and Its Utility for Assessing Taxonomic Relationships.

Genotyping by sequencing (GBS) was used to analyze relationships among cowpea and asparagus bean landraces from southern Italy and to assess the utility of this technology to study taxonomy in a wider panel, including V. unguiculata cultigroups, subspecies, and other Vigna species. The analysis of SNPs derived from GBS highlighted that, among the cowpea landraces, the African samples were separated from the other material, while, for the Italian landraces, a certain clustering depending on seed color/pattern was observed in the dendrogram. When examining the V. unguiculata species complex, a clear separation between the two groups of wild subspecies, i.e., the allogamous wild perennials and the perennial out/inbreds, could be observed, the former representing the more ancestral wild progenitors of V. unguiculata. The species V. vexillata appeared more closely related to V. unguiculata than to the other Vigna species analyzed.

A draft genome of sweet cherry (Prunus avium L.) reveals genome-wide and local effects of domestication.

Sweet cherry (Prunus avium L.) trees are both economically important fruit crops but also important components of natural forest ecosystems in Europe, Asia and Africa. Wild and domesticated trees currently coexist in the same geographic areas with important questions arising on their historical relationships. Little is known about the effects of the domestication process on the evolution of the sweet cherry genome. We assembled and annotated the genome of the cultivated variety "Big Star*" and assessed the genetic diversity among 97 sweet cherry accessions representing three different stages in the domestication and breeding process (wild trees, landraces and modern varieties). The genetic diversity analysis revealed significant genome-wide losses of variation among the three stages and supports a clear distinction between wild and domesticated trees, with only limited gene flow being detected between wild trees and domesticated landraces. We identified 11 domestication sweeps and five breeding sweeps covering, respectively, 11.0 and 2.4 Mb of the P. avium genome. A considerable fraction of the domestication sweeps overlaps with those detected in the related species, Prunus persica (peach), indicating that artificial selection during domestication may have acted independently on the same regions and genes in the two species. We detected 104 candidate genes in sweep regions involved in different processes, such as the determination of fruit texture, the regulation of flowering and fruit ripening and the resistance to pathogens. The signatures of selection identified will enable future evolutionary studies and provide a valuable resource for genetic improvement and conservation programs in sweet cherry.

The Use of Nitrogen and Its Regulation in Cereals: Structural Genes, Transcription Factors, and the Role of miRNAs.

Cereals and, especially, rice, maize, and wheat, are essential commodities, on which human nutrition is based. Expanding population and food demand have required higher production which has been achieved by increasing fertilization, and especially nitrogen supply to cereal crops. In fact, nitrogen is a crucial nutrient for the plant, but excessive use poses serious environmental and health issues. Therefore, increasing nitrogen use efficiency in cereals is of pivotal importance for sustainable agriculture. The main steps in the use of nitrogen are uptake and transport, reduction and assimilation, and translocation and remobilization. Many studies have been carried out on the genes involved in these phases, and on transcription factors regulating these genes. Lately, increasing attention has been paid to miRNAs responding to abiotic stress, including nutrient deficiency. Many miRNAs have been found to regulate transcription factors acting on the expression of specific genes for nitrogen uptake or remobilization. Recent studies on gene regulatory networks have also demonstrated that miRNAs can interact with several nodes in the network, functioning as key regulators in nitrogen metabolism.

Overexpression of the MYB29 transcription factor affects aliphatic glucosinolate synthesis in Brassica oleracea.

KEY MESSAGE: Overexpression of BoMYB29 gene up-regulates the aliphatic glucosinolate pathway in Brassica oleracea plants increasing the production of the anti-cancer metabolite glucoraphanin, and the toxic and pungent sinigrin. Isothiocyanates, the bio-active hydrolysis products of glucosinolates, naturally produced by several Brassicaceae species, play an important role in human health and agriculture. This study aims at correlating the content of aliphatic glucosinolates to the expression of genes involved in their synthesis in Brassica oleracea, and perform functional analysis of BoMYB29 gene. To this purpose, three genotypes were used: a sprouting broccoli, a cabbage, and a wild genotype (Winspit), a high glucosinolate containing accession. Winspit showed the highest transcript level of BoMYB28, BoMYB29 and BoAOP2 genes, and BoAOP2 expression was positively correlated with that of the two MYB genes. Further analyses of the aliphatic glucosinolates also showed a positive correlation between the expression of BoAOP2 and the production of sinigrin and gluconapin in Winspit. The Winspit BoMYB29 CDS was cloned and overexpressed in Winspit and in the DH AG1012 line. Overexpressing Winspit plants produced higher quantities of alkenyl glucosinolates, such as sinigrin. Conversely, the DH AG1012 transformants showed a higher production of methylsulphinylalkyl glucosinolates, including glucoraphanin, and, despite an up-regulation of the aliphatic glucosinolate genes, no increase in alkenyl glucosinolates. The latter may be explained by the absence of a functional AOP2 gene in DH AG1012. Nevertheless, an extract of DH AG1012 lines overexpressing BoMYB29 provided a chemoprotective effect on human colon cells. This work exemplifies how the genetic diversity of B. oleracea may be used by breeders to select for higher expression of transcription factors for glucosinolate biosynthesis to improve its natural, health-promoting properties.

Durum wheat genome highlights past domestication signatures and future improvement targets.

The domestication of wild emmer wheat led to the selection of modern durum wheat, grown mainly for pasta production. We describe the 10.45 gigabase (Gb) assembly of the genome of durum wheat cultivar Svevo. The assembly enabled genome-wide genetic diversity analyses revealing the changes imposed by thousands of years of empirical selection and breeding. Regions exhibiting strong signatures of genetic divergence associated with domestication and breeding were widespread in the genome with several major diversity losses in the pericentromeric regions. A locus on chromosome 5B carries a gene encoding a metal transporter (TdHMA3-B1) with a non-functional variant causing high accumulation of cadmium in grain. The high-cadmium allele, widespread among durum cultivars but undetected in wild emmer accessions, increased in frequency from domesticated emmer to modern durum wheat. The rapid cloning of TdHMA3-B1 rescues a wild beneficial allele and demonstrates the practical use of the Svevo genome for wheat improvement.

Genotyping-by-sequencing highlights patterns of genetic structure and domestication in artichoke and cardoon.

Exploiting the biodiversity of crops and their wild relatives is fundamental for maintaining and increasing food security. The species Cynara cardunculus includes three taxa: the globe artichoke, one of the most important Mediterranean vegetables, the leafy cardoon, and the wild cardoon. In this study, genotyping by sequencing (GBS) was successfully applied to reveal thousands of polymorphisms in a C. cardunculus germplasm collection, including 65 globe artichoke, 9 leafy cardoon, and 21 wild cardoon samples. The collection showed a strong population structure at K = 2, separating the globe artichoke from the leafy and wild cardoon. At higher K values, further substructures were observed, in which the wild cardoon was separated from the leafy cardoon, and the latter included the Spanish wild cardoons, while the wild sample from Portugal was admixed. Moreover, subpopulations within the globe artichoke set were highlighted. Structure analysis restricted to the globe artichoke dataset pointed out genetic differentiation between the ˝Catanesi˝ typology and all the other samples (K = 2). At higher values of K, the separation of the ˝Catanesi˝ group still held true, and green headed landraces from Apulia region, Italy (˝Green Apulian˝) formed a distinct subpopulation. ˝Romaneschi˝ artichoke types fell in a variable group with admixed samples, indicating that they should not be considered as a genetically uniform typology. The results of principal component analysis and Neighbor-Joining hierarchical clustering were consistent with structure results, and in addition provided a measure of genetic relationships among individual genotypes. Both analyses attributed the wild material from Spain and Portugal to the cultivated cardoon group, supporting the idea that this might be indeed a feral form of the leafy cardoon. Different reproductive habit and possibly selective pressure led to a slower LD decay in artichoke compared to cardoon. Genotyping by sequencing has proven a reliable methodology to obtain valuable SNPs and assess population genetics in C. cardunculus.

Isolation and Characterization of the Flavonol Regulator CcMYB12 From the Globe Artichoke [Cynara cardunculus var. scolymus (L.) Fiori].

Flavonoids are a well-studied group of secondary metabolites, belonging to the phenylpropanoid pathway. Flavonoids are known to exhibit health promoting effects such as antioxidant capacities, anti-cancer and anti-inflammatory activity. Globe artichoke is an important source of bioactive phenolic compounds, including flavonoids. To study the regulation of their biosynthesis, a R2R3-MYB transcription factor, CcMYB12, was isolated from artichoke leaves. Phylogenetic analysis showed that this protein belongs to the MYB subgroup 7 (flavonol-specific MYB), which includes Arabidopsis AtMYB12, grapevine VvMYBF1, and tomato SlMYB12. CcMYB12 transcripts were detected specifically in artichoke immature inflorescence and young leaves and overlapped with the profiles of flavonol biosynthetic genes. Electrophoretic mobility shift assays (EMSAs) revealed that recombinant CcMYB12 protein is able to bind to ACII element, a DNA binding site ubiquitously present in the promoters of genes encoding flavonol biosynthetic enzymes. In transgenic Arabidopsis plants, the overexpression of CcMYB12 activated the expression of endogenous flavonol biosynthesis genes, leading to an increase of flavonol accumulation and a decrease of anthocyanins in leaves. Likewise, in transgenic tobacco petals and leaves, the overexpression of CcMYB12 decreased anthocyanin levels and increased flavonols.

MicroRNAs in durum wheat seedlings under chronic and short-term nitrogen stress.

Nitrogen is an essential macronutrient for plant growth and reproduction. In durum wheat, an appropriate nitrogen soil availability is essential for an optimal seed development. miRNAs contribute to the environmental change adaptation of plants through the regulation of important genes involved in stress processes. In this work, nitrogen stress response was evaluated in durum wheat seedlings of Ciccio and Svevo cultivars. Eight small RNA libraries from leaves and roots of chronically stressed plants were sequenced to detect conserved and novel miRNAs. A total of 294 miRNAs were identified, 7 of which were described here for the first time. The expression level of selected miRNAs and target genes was analyzed by qPCR in seedlings subjected to chronic (Ciccio and Svevo, leaves and roots) or short-term (Svevo roots) stress conditions. Some miRNAs showed an immediate stress response, and their level of expression was either maintained or returned to a basal level during a long-term stress. Other miRNAs showed a gradual up- or downregulation during the short-term stress. The newly identified miRNA ttu-novel-106 showed an immediate strongly downregulation after nitrogen stress, which was negatively correlated with the expression of MYB-A, its putative target gene. PHO2 gene was significantly upregulated after 24-48-h stress, corresponding to a downregulation of miR399b. Ttu-miR399b putative binding sites in the 5' UTR region of the Svevo PHO2 gene were identified in the A and B genomes. Both MYB-A and PHO2 genes were validated for their cleavage site using 5' RACE assay.

Durum wheat miRNAs in response to nitrogen starvation at the grain filling stage.

Durum wheat highly depends on nitrogen for seed development and yield, and the obtainment of varieties with a better nitrogen use efficiency is crucial to reduce production costs and environmental pollution. In this study, sequencing of two small RNA libraries obtained from tissues of Ciccio and Svevo cultivars grown under nitrogen starvation conditions produced 84 novel, and 161 conserved miRNAs. Of these, 7 novel and 13 known miRNAs were newly identified in this work. Quantitative PCR analysis of selected miRNAs highlighted that the expression levels of some of them depends on the tissue and on the cultivar, Svevo being the most responsive to nitrogen starvation. A number of target genes were predicted to be involved in nitrogen metabolism. An inverse correlation for the qPCR expression data of miRNA/target pairs miR399b/PHO2, miR393c/AFB2, ttu-novel-61/CCAAT-TF was observed in specific tissues or cultivar. Especially, ttu-novel-61 was down-regulated and its target CCAAT-TF up-regulated in almost all tissues both in Svevo and in Ciccio. Moreover, CCAAT-TF was confirmed to be cleaved by ttu-novel-61 at the expected site. The discovery of miRNAs involved in the response to nitrogen stress represents an important step towards functional analyses, with the final aim to design strategies for improving nitrogen use efficiency in durum wheat.

A Distinct Genetic Cluster in Cultivated Chickpea as Revealed by Genome-wide Marker Discovery and Genotyping.

The accurate description of plant biodiversity is of utmost importance to efficiently address efforts in conservation genetics and breeding. Herein, we report the successful application of a genotyping-by-sequencing (GBS) approach in chickpea ( L.), resulting in the characterization of a cultivated germplasm collection with 3187 high-quality single nucleotide polymorphism (SNP) markers. Genetic structure inference, principal component analysis, and hierarchical clustering all indicated the identification of a genetic cluster corresponding to black-seeded genotypes traditionally cultivated in Southern Italy. Remarkably, this cluster was clearly distinct at both genetic and phenotypic levels from germplasm groups reflecting commercial chickpea classification into and seed types. Fixation index estimates for individual polymorphisms pointed out loci and genomic regions that might be of significance for the diversification of agronomic and commercial traits. Overall, our findings provide information on genetic relationships within cultivated chickpea and highlight a gene pool of great interest for the scientific community and chickpea breeding, which is limited by the low genetic diversity available in the primary gene pool.

Genetic diversity and accession structure in European Cynara cardunculus collections.

Understanding the distribution of genetic variations and accession structures is an important factor for managing genetic resources, but also for using proper germplasm in association map analyses and breeding programs. The globe artichoke is the fourth most important horticultural crop in Europe. Here, we report the results of a molecular analysis of a collection including globe artichoke and leafy cardoon germplasm present in the Italian, French and Spanish gene banks. The aims of this study were to: (i) assess the diversity present in European collections, (ii) determine the population structure, (iii) measure the genetic distance between accessions; (iv) cluster the accessions; (v) properly distinguish accessions present in the different national collections carrying the same name; and (vi) understand the diversity distribution in relation to the gene bank and the geographic origin of the germplasm. A total of 556 individuals grouped into 174 accessions of distinct typologies were analyzed by different types of molecular markers, i.e. dominant (ISSR and AFLP) and co-dominant (SSR). The data of the two crops (globe artichoke and leafy cardoon) were analyzed jointly and separately to compute, among other aims, the gene diversity, heterozygosity (He, Ho), fixation indexes, AMOVA, genetic distance and structure. The findings underline the huge diversity present in the analyzed material, and the existence of alleles that are able to discriminate among accessions. The accessions were clustered not only on the basis of their typology, but also on the basis of the gene bank they come from. Probably, the environmental conditions of the different field gene banks affected germplasm conservation. These outcomes will be useful in plant breeding to select accessions and to fingerprint varieties. Moreover, the results highlight the particular attention that should be paid to the method used to conserve the Cynara cardunculus germplasm and suggest to the preference of using accessions from different gene banks to run an association map.

Differential representation of albumins and globulins during grain development in durum wheat and its possible functional consequences.

Durum wheat (Triticum turgidum ssp. durum (Desf.) Husn.) is an economically important crop used for the production of semolina, which is the basis of pasta and other food products. Its grains provide proteins and starch for human consumption. Grain development is a key process in wheat physiology; it is highly affected by a number of enzymes that control the metabolic processes governing accumulation of starch and storage proteins and ultimately grain weight. Most of these enzymes are present in the albumin/globulin grain fraction, which represents about a quarter of total seed proteins. With the aim to describe the dynamic profile of the albumin/globulin fraction during durum wheat grain development, we performed a proteomic analysis of this subproteome using a two-dimensional differential gel electrophoresis (2D-DIGE)-based approach and compared six developmental stages. A total of 285 differentially (237 over- and 48 under-) represented spots was identified by nanoLC-ESI-LIT-MS/MS, which were associated with 217 non-redundant Triticum sequence entries. Quantitative protein dynamics demonstrated that carbon metabolism, energy, protein destination/storage, disease/defense and cell growth/division functional categories were highly affected during grain development, concomitantly with progressive grain size increase and starch/protein reserve accumulation. Bioinformatic interaction prediction revealed a complex network of differentially represented proteins mainly centered at enzymes involved in carbon and protein metabolism. A description of 18 proteins associated with wheat flour human allergies was also obtained; these components showed augmented levels at the last developmental stages. By providing a comprehensive understanding of the molecular basis of durum wheat grain development, yield and quality formation, this study provides the foundation and reveals potential biomarkers for further investigations of durum wheat breeding and semolina quality. BIOLOGICAL SIGNIFICANCE: A 2D-DIGE-based comparative analysis of the albumin/globulin fraction from durum wheat caryopses at six developmental stages was performed to describe the dynamic subproteomic changes associated with grain development. Quantitative variations of 217 differentially proteins demonstrated that highly affected are the functional categories of carbon metabolism, energy, protein destination/storage, disease/defense and cell growth/division, which displayed a general over-representation, consistently with concomitant occurrence of grain size increase and starch/protein reserve accumulation. Bioinformatics revealed a complex protein network centered mainly at enzymes involved in carbon and protein metabolism. Differentially represented proteins and corresponding functional categories highly resembled those previously identified as variable in developing bread wheat grain. This suggests that the main differences in kernel hardness between durum and bread wheat probably do not depend on proteomic changes in corresponding albumins/globulins, but on other specific factors affecting the interaction between the starch granules and the endosperm protein matrix in the kernel.

Transcriptomic response of durum wheat to nitrogen starvation.

Nitrogen (N) is a key macronutrient representing a limiting factor for plant growth and development and affects productivity in wheat. In this study, durum wheat response to N chronic starvation during grain filling was investigated through a transcriptomic approach in roots, leaves/stems, flag leaf and spikes of cv. Svevo. Nitrogen stress negatively influenced plant height, tillering, flag leaf area, spike and seed traits, and total N content. RNA-seq data revealed 4,626 differentially expressed genes (DEGs). Most transcriptomic changes were observed in roots, with 3,270 DEGs, while 963 were found in leaves/stems, 470 in flag leaf, and 355 in spike tissues. A total of 799 gene ontology (GO) terms were identified, 180 and 619 among the upregulated and downregulated genes, respectively. Among the most addressed GO categories, N compound metabolism, carbon metabolism, and photosynthesis were mostly represented. Interesting DEGs, such as N transporters, genes involved in N assimilation, along with transcription factors, protein kinases and other genes related to stress were highlighted. These results provide valuable information about the transcriptomic response to chronic N stress in durum wheat, which could be useful for future improvement of N use efficiency.

The miRNAome of durum wheat: isolation and characterisation of conserved and novel microRNAs and their target genes.

BACKGROUND: The allotetraploid durum wheat [Triticum turgidum subsp. durum (Desf.) Husn.] is a highly economically important species especially in the Mediterranean basin. However, its genomics, transcriptomics and in particular microRNAome are still largely unknown. RESULTS: In the present work, two small RNA libraries from durum wheat Ciccio and Svevo cultivars were generated from different tissues at the late milk (Z77) developmental stage. A total of 167 conserved and 98 potential novel miRNAs were identified in the two libraries and interestingly, three novel miRNAs were found to be derived from ribosomal RNA. Putative target genes were predicted for conserved and novel miRNAs, the majority of which interact with nucleic acids, according to GO terms relative to molecular function. Quantitative qPCR analysis showed that several miRNAs identified were differentially expressed in the mature (Z77) developmental stage compared to young (Z14) tissues. Moreover, target gene expression analysis suggested that in roots, the putative genes encoding for the SQUAMOSA SPL2 and TGA1 proteins are regulated by ttu-miR156n, while MYB3 transcription factor by ttu-miR319f. Additionally, the Photosystem II P680 chlorophyll A apoprotein gene showed an expression level negatively correlated to that of ttu-novel-48 in leaves. CONCLUSION: Our results suggest that, in durum wheat, these genes may play important roles in root/leaf development and are subjected to miRNA regulation. The prediction of novel miRNAs putatively derived from ribosomal RNA opens new perspectives on the study of plant miRNAs.