Pseudogenes can provide valuable insight into the evolutionary history of genomes which may be difficult to determine based on the examination of functional loci alone. This study presents the results of a comprehensive, two-stage genome-wide investigation into the identification and characterization of pseudogenes in Fusarium graminearum, the main causal agent of wheat head blight. By analyzing sequence homology between non-coding regions of the genome and predicted protein sequences, regions with homology to putative paralogous functional sequences were identified. These regions were characterized in terms of their matching sequence structure and position. Most of the pseudogenes identified were mapped within the fast-evolving genomic compartment and were derived from transposition events. The number of processed and putatively retroposed pseudogenes was found to be comparable. The low number of identified pseudogenes is consistent with the low number of gene duplicates in F. graminearum. No compelling evidence was found to suggest that pseudogene formation can be explained by evolutionary accidents during gene family expansion, nor that it is caused by RIP-associated mutagenic events. Notably, approximately one-third (144 of 436) of the pseudogenes were found to overlap with the untranslated or intron sequences of functional loci, suggesting that they could potentially be transcribed. Using comparative Fusarium genomics, we identified genomic regions with homology to genes lacking functional orthologues in F. graminearum. These were investigated as putative unitary pseudogenes, some of which had lost their original functions completely after the radiation of F. graminearum. Interestingly, the orthologous loci matching parents of the 18 loss-of-function pseudogenes showed homology to domains previously identified in proteins involved in pathogenesis.

Comparative analyses of several plant pathogens have revealed that genome plasticity could be associated with different genomic architectures. In certain species, genomic compartments are characterised by highly conserved regions that contain mainly housekeeping genes and rearranged regions that are enriched for genes related to virulence and adaptation. The compositional and structural characteristics of genomic regions have been significantly associated with compartment membership in single species, but little information is available on the covariation of these features between species. Here, the results of a comparative analysis of the genomic architectures of six agriculturally relevant Fusarium species, which differ for several biological and pathogenic characteristics, are presented. These include F. culmorum , F. fujikoroi , F. graminearum, F. oxysporum , F. solani, and F. verticillioides . The genome sequences of these species were partitioned into adjacent windows, with the average level of gene collinearity with the other species used as an index of compartment membership. High collinearity is typical of conserved regions, while low collinearity is typical of rearranged regions. Several genic and genomic variables were found to be consistently associated with compartment definition among all the Fusarium species that were investigated. The compartment that was characterised by lower collinearity (i.e., high genomic rearrangements) contained more relocated genes, species-specific genes and secreted protein-encoding genes than regions with low collinearity. Furthermore, several molecular evidence indicates that low-collinearity regions are more likely to be subjected to selective pressure than high-collinearity regions. Indeed, genes residing in the former regions exhibited higher rates of sequence evolution than in the latter, as indicated by the high non-synonymous-to-synonymous substitution rates. However, they exhibited signatures of selection to minimise the costs of transcription, as indicated by their high coding density. Our data suggests that although variable genomic compartments evolved mostly after species radiation, they share similar genomic features across related species and perhaps evolve with similar mechanisms. ### Competing Interest Statement The authors have declared no competing interest.

The knowledge of the organization of the domesticated gene pool of crop species is an essential requirement to understand crop evolution, to rationalize conservation programs, and to support practical decisions in plant breeding. Here, we integrate simple sequence repeat (SSR) analysis and phenotypic characterization to investigate a globe artichoke collection that comprises most of the varieties cultivated worldwide. We show that the cultivated gene pool of globe artichoke includes five distinct genetic groups associated with the major phenotypic typologies: Catanesi (which based on our analysis corresponds to Violetti di Provenza), Spinosi, Violetti di Toscana, Romaneschi, and Macau. We observed that 17 and 11% of the molecular and phenotypic variance, respectively, is between these groups, while within groups, strong linkage disequilibrium and heterozygote excess are evident. The divergence between groups for quantitative traits correlates with the average broad-sense heritability within the groups. The phenotypic divergence between groups for both qualitative and quantitative traits is strongly and positively correlated with SSR divergence (FST) between groups. All this implies a low population size and strong bottleneck effects, and indicates a long history of clonal propagation and selection during the evolution of the domesticated gene pool of globe artichoke. Moreover, the comparison between molecular and phenotypic population structures suggests that harvest time, plant architecture (i.e., plant height, stem length), leaf spininess, head morphology (i.e., head shape, bract shape, spininess) together with the number of heads per plant were the main targets of selection during the evolution of the cultivated germplasm. We emphasize our findings in light of the potential exploitation of this collection for association mapping studies.

In legumes, pod shattering occurs when mature pods dehisce along the sutures, and detachment of the valves promotes seed dispersal. In Phaseolus vulgaris (L)., the major locus qPD5.1-Pv for pod indehiscence was identified recently. We developed a BC4/F4 introgression line population and narrowed the major locus down to a 22.5 kb region. Here, gene expression and a parallel histological analysis of dehiscent and indehiscent pods identified an AtMYB26 orthologue as the best candidate for loss of pod shattering, on a genomic region ~11 kb downstream of the highest associated peak. Based on mapping and expression data, we propose early and fine up-regulation of PvMYB26 in dehiscent pods. Detailed histological analysis establishes that pod indehiscence is associated with the lack of a functional abscission layer in the ventral sheath, and that the key anatomical modifications associated with pod shattering in common bean occur early during pod development. We finally propose that loss of pod shattering in legumes resulted from histological convergent evolution and that it is the result of selection at orthologous loci.

By selecting for prostrate growth habit of the juvenile phase of the cycle, durum wheat cultivars could be developed with improved competitive ability against weeds, and better soil coverage to reduce the soil water lost by evaporation. A panel of 184 durum wheat (Triticum turgidum subsp. durum) genotypes, previously genotyped with DArT-seq markers, was used to perform association mapping analysis of prostrate/erect growth habit trait and to identify candidate genes. Phenotypic data of plant growth habit were recorded during three consecutive growing seasons (2014–2016), two different growth conditions (field trial and greenhouse) and two sowing periods (autumn and spring). Genome-wide association study revealed significant marker-trait associations, twelve of which were specific for a single environment/year, 4 consistent in two environments, and two MTAs for the LSmeans were identified across all environments, on chromosomes 2B and 5A. The co-localization of some MTAs identified in this study with known vernalization and photoperiod genes demonstrated that the sensitivity to vernalization and photoperiod response are actually not only key components of spring/winter growth habit, but they play also an important role in defining the magnitude of the tiller angle during the tillering stage. Many zinc-finger transcription factors, such as C2H2 or CCCH-domain zinc finger proteins, known to be involved in plant growth habit and in leaf angle regulation were found as among the most likely candidate genes. The highest numbers of candidate genes putatively related to the trait were found on chromosomes 3A, 4B, 5A and 6A. Moreover, a bioinformatic approach has been considered to search for functional ortholog genes in wheat by using the sequence of rice and barley tiller angle-related genes. The information generated could be used to improve the understanding of the mechanisms that regulate the prostrate/erect growth habit in wheat and the adaptive potential of durum wheat under resource-limited environmental conditions.

Tomato (Solanum lycopersicum L.) is a widely used model plant species for dissecting out the genomic bases of complex traits to thus provide an optimal platform for modern "-omics" studies and genome-guided breeding. Genome-wide association studies (GWAS) have become a preferred approach for screening large diverse populations and many traits. Here, we present GWAS analysis of a collection of 115 landraces and 11 vintage and modern cultivars. A total of 26 conventional descriptors, 40 traits obtained by digital phenotyping, the fruit content of six carotenoids recorded at the early ripening (breaker) and red-ripe stages and 21 climate-related variables were analyzed in the context of genetic diversity monitored in the 126 accessions. The data obtained from thorough phenotyping and the SNP diversity revealed by sequencing of ripe fruit transcripts of 120 of the tomato accessions were jointly analyzed to determine which genomic regions are implicated in the expressed phenotypic variation. This study reveals that the use of fruit RNA-Seq SNP diversity is effective not only for identification of genomic regions that underlie variation in fruit traits, but also of variation related to additional plant traits and adaptive responses to climate variation. These results allowed validation of our approach because different marker-trait associations mapped on chromosomal regions where other candidate genes for the same traits were previously reported. In addition, previously uncharacterized chromosomal regions were targeted as potentially involved in the expression of variable phenotypes, thus demonstrating that our tomato collection is a precious reservoir of diversity and an excellent tool for gene discovery.

Loss of seed shattering is a key trait in crop domestication, particularly for grain crops. For wild plants, seed shattering is a crucial mechanism to achieve greater fitness, although in the agricultural context, this mechanism reduces harvesting efficiency, especially under dry conditions. Loss of seed shattering was acquired independently in different monocotyledon and dicotyledon crop species by ‘convergent phenotypic evolution’, leading to similar low dehiscent and indehiscent phenotypes. Here, the main aim is to review the current knowledge about seed shattering in crops, in order to highlight the tissue modifications that underlie the convergent phenotypic evolution of reduced shattering in different types of fruit, from the silique of Brassicaceae species, to the pods of legumes and spikes of cereals. Emphasis is given to legumes, with consideration of recent data obtained for the common bean. The current review also discusses to what extent convergent phenotypes arose from parallel changes at the histological and/or molecular levels. For this reason, an overview is included of the main findings relating to the genetic control of seed shattering in the model species Arabidopsis thaliana and in other important crops.

Bakanae disease, caused by the fungal phytopathogen Fusarium fujikuroi, can be detected in most rice (Oryza sativa L.) growing areas worldwide. In this study, we investigated the population structure of this fungus in southern Lao PDR, a country located near the geographic origin of rice domestication. Microsatellites (SSRs) and mating type (MAT) analyses, pathogenicity and fungicide sensitivity tests were integrated in the study. The first key finding is that the population genetic structure of F. fujikuroi in Lao PDR is consistent with high clonal reproduction. Indeed, (i) "true" clones were identified; (ii) within populations, MAT types were frequently skewed from 1:1 ratio, (iii) linkage disequilibrium (among SSRs as also among SSRs and MAT) was present, and (iv) gene-flow between opposite MAT types within the same population is restricted. The presence of genetic divergence among areas and populations and the occurrence of positive spatial autocorrelation of genetic variation, indicate that migration is restricted, and that genetic drift plays an important role in the evolution of this fungus. Two main well-defined groups of isolates were detected (FST = 0.213) that display a non-random spatial distribution. They differ in the ability to induce seedlings death but not seedlings elongation (the typical Bakanae symptom) suggesting that the pathogen's ability to induce the two symptoms is under different genetic control. Finally, we compared two agroecosystems with contrasting characteristics: low-input and traditional (Lao PDR) vs high-input and modern (Italy). We found differences in the level of population structuring and of spatial autocorrelation. This suggests that the evolutionary potential of the fungus not only depends on its intrinsic characteristics, but is strongly influenced by other external factors, most likely by the dynamics of infested seed exchange. Thus, quarantine and chemical treatments are a way to reduce population connectivity and hence the evolutionary potential of this pathogen.

Understanding the population genetic structure of crop gene pool is a fundamental prerequisite in plant evolutionary biology and plant breeding. Based on this consideration, in this chapter we pursued two main aims. The first was to review the status of the research pointed at deciphering the population structure of the gene pools of globe artichoke. The second was to show how the use of a well-characterized germplasm collection together with deep genotyping could serve the identification of QTLs for traits of agronomic and economic relevance in globe artichoke.

The complete or partial loss of shattering ability occurred independently during the domestication of several crops. Therefore, the study of this trait can provide an understanding of the link between phenotypic and molecular convergent evolution. The genetic dissection of 'pod shattering' in Phaseolus vulgaris is achieved here using a population of introgression lines and next-generation sequencing techniques. The 'occurrence' of the indehiscent phenotype (indehiscent versus dehiscent) depends on a major locus on chromosome 5. Furthermore, at least two additional genes are associated with the 'level' of shattering (number of shattering pods per plant: low versus high) and the 'mode' of shattering (non-twisting versus twisting pods), with all of these loci contributing to the phenotype by epistatic interactions. Comparative mapping indicates that the major gene identified on common bean chromosome 5 corresponds to one of the four quantitative trait loci for pod shattering in Vigna unguiculata. None of the loci identified comprised genes that are homologs of the known shattering genes in Glycine max. Therefore, although convergent domestication can be determined by mutations at orthologous loci, this was only partially true for P. vulgaris and V. unguiculata, which are two phylogenetically closely related crop species, and this was not the case for the more distant P. vulgaris and G. max. Conversely, comparative mapping suggests that the convergent evolution of the indehiscent phenotype arose through mutations in different genes from the same underlying gene networks that are involved in secondary cell-wall biosynthesis and lignin deposition patterning at the pod level.

Fusarium culmorum is the most commonly reported root rot pathogen in Tunisian durum wheat. Isolates of the pathogen from four durum wheat growing areas in the north of Tunisia were analyzed for their chemotypes. Two chemotypes were detected at unequal abundance (96% of 3-ADON and 4% of NIV). Distribution of a SNP mutation located at the position 34 bp after the first exon of the EF-1α partial sequence was analysed, to verify whether the haplotype was specifically associated to Fusarium root rot. A and T haplotypes were homogeneously distributed in three different Tunisian regions (Mateur, Beja and Bousalem) but not for the region of Bizerte, from which greatest number of A haplotype strains were detected. The isolates were tested for their virulence under glasshouse conditions, and a mean of 91% of crown and root infection was observed. Chemotype influenced virulence, but there was no significant influence of the geographical origin or haplotype on virulence. The distribution of three inter simple sequence repeats (ISSR) was examined, to better understand the structure of F. culmorum populations in Tunisia. A total of 27 fragments were obtained with eight polymorphic bands. Cluster analysis showed a high level of similarity between isolates. Analysis of molecular variance confirmed that there was little genetic differentiation among F. culmorum strains from different locations.

Seed shattering in crops is a key domestication trait due to its relevance for seed dispersal, yield, and fundamental questions in evolution (e.g., convergent evolution). Here, we focused on pod shattering in common bean (Phaseolus vulgaris L.), the most important legume crop for human consuption in the world. With this main aim, we developed a methodological pipeline that comprises a thorough characterization under field conditions, including also the chemical composition and histological analysis of the pod valves. The pipeline was developed based on the assumption that the shattering trait itself can be treated in principle as a “syndrome” (i.e., a set of correlated different traits) at the pod level. We characterized a population of 267 introgression lines that were developed ad-hoc to study shattering in common bean. Three main objectives were sought: (1) to dissect the shattering trait into its “components,” of level (percentage of shattering pods per plant) and mode (percentage of pods with twisting or non-twisting valves); (2) to test whether shattering is associated to the chemical composition and/or the histological characteristics of the pod valves; and (3) to test the associations between shattering and other plant traits. We can conclude the following: Very high shattering levels can be achieved in different modes; shattering resistance is mainly a qualitative trait; and high shattering levels is correlated with high carbon and lignin contents of the pod valves and with specific histological charaterstics of the ventral sheath and the inner fibrous layer of the pod wall. Our data also suggest that shattering comes with a “cost,” as it is associated with low pod size, low seed weight per pod, high pod weight, and low seed to pod-valves ratio; indeed, it can be more exaustively described as a syndrome at the pod level. Our work suggests that the valve chemical composition (i.e., carbon and lignin content) can be used for a high troughput phenotyping procedures for shattering phenotyping. Finally, we believe that the application of our pipeline will greatly facilitate comparative studies among legume crops, and gene tagging.

Crop species have been deeply affected by the domestication process, and there have been many efforts to identify selection signatures at the genome level. This knowledge will help geneticists to better understand the evolution of organisms, and at the same time, help breeders to implement successful breeding strategies. Here, we focused on domestication in the Mesoamerican gene pool of Phaseolus vulgaris by sequencing 49 gene fragments from a sample of 45 P. vulgaris wild and domesticated accessions, and as controls, two accessions each of the closely related species Phaseolus coccineus and Phaseolus dumosus. An excess of nonsynonymous mutations within the domesticated germplasm was found. Our data suggest that the cost of domestication alone cannot explain fully this finding. Indeed, the significantly higher frequency of polymorphisms in the coding regions observed only in the domesticated plants (compared to noncoding regions), the fact that these mutations were mostly nonsynonymous and appear to be recently derived mutations, and the investigations into the functions of their relative genes (responses to biotic and abiotic stresses), support a scenario that involves new functional mutations selected for adaptation during domestication. Moreover, consistent with this hypothesis, selection analysis and the possibility to compare data obtained for the same genes in different studies of varying sizes, data types, and methodologies allowed us to identify four genes that were strongly selected during domestication. Each selection candidate is involved in plant resistance/tolerance to abiotic stresses, such as heat, drought, and salinity. Overall, our study suggests that domestication acted to increase functional diversity at target loci, which probably controlled traits related to expansion and adaptation to new agro-ecological growing conditions.

A new era has begun for Phaseolus vulgaris, and other Phaseolus spp., with the release of the reference genomes of both the Mesoamerican and Andean genotypes. Exploiting the new genome sequences information and the derived tools, important insight in the common bean genomics can be achieved. A major breakthrough in common bean will be the identification of the molecular basis of the domestication syndrome, representing a main step towards our understanding of the evolutionary processes and a fundamental support for researchers and breeders involved in crop improvement. P. vulgaris, along with the other Phaseolus species, represents a unique model to study evolution, domestication and environmental adaptation, focusing on the major phenotypic changes occurring during its evolutionary histories and trajectories and unveiling the molecular mechanisms and genetic basis responsible for the observed changes. This chapter offers an overview of the current knowledge of the evolutionary history of common bean, and of the outcomes relating to the genetic bases of important domestication and adaptation traits. We provide an analysis of the process of domestication, with the focus on convergent phenotypic evolution, and a survey of the studies on common bean that have been specifically carried out on genes related to the domestication syndrome; with particular focus on studies that have compared wild and domesticated forms, highlighting findings on the genetic control of the domestication syndrome and on the genetic architecture of environmental and agronomic important traits.

Here, we aim to provide a comprehensive and up-to-date overview of the most significant outcomes in the literature regarding the origin of Phaseolus genus, the geographical distribution of the wild species, the domestication process, and the wide spread out of the centers of origin. Phaseolus can be considered as a unique model for the study of crop evolution, and in particular, for an understanding of the convergent phenotypic evolution that occurred under domestication. The almost unique situation that characterizes the Phaseolus genus is that five of its ∼70 species have been domesticated (i.e., Phaseolus vulgaris, P. coccineus, P. dumosus, P. acutifolius, and P. lunatus), and in addition, for P. vulgaris and P. lunatus, the wild forms are distributed in both Mesoamerica and South America, where at least two independent and isolated episodes of domestication occurred. Thus, at least seven independent domestication events occurred, which provides the possibility to unravel the genetic basis of the domestication process not only among species of the same genus, but also between gene pools within the same species. Along with this, other interesting features makes Phaseolus crops very useful in the study of evolution, including: (i) their recent divergence, and the high level of collinearity and synteny among their genomes; (ii) their different breeding systems and life history traits, from annual and autogamous, to perennial and allogamous; and (iii) their adaptation to different environments, not only in their centers of origin, but also out of the Americas, following their introduction and wide spread through different countries. In particular for P. vulgaris this resulted in the breaking of the spatial isolation of the Mesoamerican and Andean gene pools, which allowed spontaneous hybridization, thus increasing of the possibility of novel genotypes and phenotypes. This knowledge that is associated to the genetic resources that have been conserved ex situ and in situ represents a crucial tool in the hands of researchers, to preserve and evaluate this diversity, and at the same time, to identify the genetic basis of adaptation and to develop new improved varieties to tackle the challenges of climate change, and food security and sustainability.

Domestication and breeding have influenced the genetic structure of plant populations due to selection for adaptation from natural habitats to agro-ecosystems. Here, we investigate the effects of selection on the contents of 51 primary kernel metabolites and their relationships in three Triticum turgidum L. subspecies (i.e., wild emmer, emmer, durum wheat) that represent the major steps of tetraploid wheat domestication. We present a methodological pipeline to identify the signature of selection for molecular phenotypic traits (e.g., metabolites and transcripts). Following the approach, we show that a reduction in unsaturated fatty acids was associated with selection during domestication of emmer (primary domestication). We also show that changes in the amino acid content due to selection mark the domestication of durum wheat (secondary domestication). These effects were found to be partially independent of the associations that unsaturated fatty acids and amino acids have with other domestication-related kernel traits. Changes in contents of metabolites were also highlighted by alterations in the metabolic correlation networks, indicating wide metabolic restructuring due to domestication. Finally, evidence is provided that wild and exotic germplasm can have a relevant role for improvement of wheat quality and nutritional traits.

Nuclear and chloroplast markers and phenotypic characters were integrated to analyse the population genetic structure of wild cardoon, Cynara cardunculus var. sylvestris, the ancestor of cultivated globe artichoke, Cynara cardunculus var. scolymus on the island of Sardinia, Italy. The spatial scale ranged from a few metres to ∼200km. Wild cardoon appears to be genetically fragmented, with significant genetic divergence at various scales, indicating that gene flow is insufficient to counterbalance the effects of genetic drift or founder effects. Divergence between populations was higher for chloroplast (40%) than for nuclear markers (15%), suggesting that gene flow via seed was lower than via pollen. Two main genetic groups were detected; these correlated with differences in flowering time, capitula size, glossiness, and anthocyanin pigmentation. A complex population structure of wild cardoon emerged over small spatial scales, likely resulting from the interplay between gene dispersal, colonisation history and selective forces. Indeed, Sardinia appears to be a 'hybrid zone' of different gene pools. The island has unique diverse germplasm that has originated from hybridisation among different gene pools. The sampling of seeds from a few plants but from many sites is suggested as the best strategy to harvest the genetic diversity of wild cardoon.

Here we present the first insights on a genome-wide association approach in globe artichoke, based on an association panel composed of 111 varietal types coming from the living collection maintained at AGRIS (Sardinia, Italy). The collection has been investigated through two complementary molecular approaches: (i) microsatellite (SSRs) genotyping, using 26 loci uniformly distributed throughout the genome to investigate the population structure of the core germplasm set; and (ii) genotyping by sequencing, using a two-enzymes RADseq (restriction-site associated DNA sequencing) technology, as a basis for an association mapping approach.

We have investigated the role of selection in the determination of the detected levels of introgression from modern maize hybrid varieties into maize landraces still cultivated in situ in Italy. We exploited the availability of a historical collection of landraces undertaken before the introduction and widespread use of modern maize, to analyse genomic changes that have occurred in these maize landraces over 50 years of co-existence with hybrid varieties. We have combined a previously published SSR dataset (n=21) with an AFLP loci dataset (n=168) to provide higher resolution power and to obtain a more detailed picture. We show that selection pressures for adaptation have favoured new alleles introduced by migration from hybrids. This shows the potential for analysis of historical introgression even over this short period of 50 years, for an understanding of the evolution of the genome and for the identification of its functionally important regions. Moreover, this demonstrates that landraces grown in situ represent almost unique populations for use for such studies when the focus is on the domesticated plant. This is due to their adaptation, which has arisen from their dynamic evolution under a continuously changing agro-ecological environment, and their capture of new alleles from hybridisation. We have also identified loci for which selection has inhibited introgression from modern germplasm and has enhanced the distinction between landraces and modern maize. These loci indicate that selection acted in the past, during the formation of the flint and dent gene pools. In particular, the locus showing the strongest signals of selection is a Misfit transposable element. Finally, molecular characterisation of the same samples with two different molecular markers has allowed us to compare their performances. Although the genetic-diversity and population-structure analyses provide the same global qualitative pattern, which thus provides the same inferences, there are differences related to their natures and characteristics.