Transcriptomic response to nitrogen availability reveals signatures of adaptive plasticity during tetraploid wheat domestication.

The domestication of crops, coupled with agroecosystem development, is associated with major environmental changes and provides an ideal model of phenotypic plasticity. Here, we examined 32 genotypes of three tetraploid wheat (Triticum turgidum L.) subspecies, wild emmer, emmer, and durum wheat, which are representative of the key stages in the domestication of tetraploid wheat. We developed a pipeline that integrates RNA-Seq data and population genomics to assess gene expression plasticity and identify selection signatures under diverse nitrogen availability conditions. Our analysis revealed differing gene expression responses to nitrogen availability across primary (wild emmer to emmer) and secondary (emmer to durum wheat) domestication. Notably, nitrogen triggered the expression of twice as many genes in durum wheat compared to that in emmer and wild emmer. Unique selection signatures were identified at each stage: primary domestication mainly influenced genes related to biotic interactions, whereas secondary domestication affected genes related to amino acid metabolism, in particular lysine. Selection signatures were found in differentially expressed genes (DEGs), notably those associated with nitrogen metabolism, such as the gene encoding glutamate dehydrogenase (GDH). Overall, our study highlights the pivotal role of nitrogen availability in the domestication and adaptive responses of a major food crop, with varying effects across different traits and growth conditions.

CRISPR-Cas9-based repeat depletion for high-throughput genotyping of complex plant genomes.

High-throughput genotyping enables the large-scale analysis of genetic diversity in population genomics and genome-wide association studies that combine the genotypic and phenotypic characterization of large collections of accessions. Sequencing-based approaches for genotyping are progressively replacing traditional genotyping methods because of the lower ascertainment bias. However, genome-wide genotyping based on sequencing becomes expensive in species with large genomes and a high proportion of repetitive DNA. Here we describe the use of CRISPR-Cas9 technology to deplete repetitive elements in the 3.76-Gb genome of lentil (Lens culinaris), 84% consisting of repeats, thus concentrating the sequencing data on coding and regulatory regions (single-copy regions). We designed a custom set of 566,766 gRNAs targeting 2.9 Gbp of repeats and excluding repetitive regions overlapping annotated genes and putative regulatory elements based on ATAC-seq data. The novel depletion method removed ∼40% of reads mapping to repeats, increasing those mapping to single-copy regions by ∼2.6-fold. When analyzing 25 million fragments, this repeat-to-single-copy shift in the sequencing data increased the number of genotyped bases of ∼10-fold compared to nondepleted libraries. In the same condition, we were also able to identify ∼12-fold more genetic variants in the single-copy regions and increased the genotyping accuracy by rescuing thousands of heterozygous variants that otherwise would be missed because of low coverage. The method performed similarly regardless of the multiplexing level, type of library or genotypes, including different cultivars and a closely related species (L. orientalis). Our results showed that CRISPR-Cas9-driven repeat depletion focuses sequencing data on single-copy regions, thus improving high-density and genome-wide genotyping in large and repetitive genomes.

Decreased metabolic diversity in common beans associated with domestication revealed by untargeted metabolomics, information theory, and molecular networking.

The process of crop domestication leads to a dramatic reduction in the gene expression associated with metabolic diversity. Genes involved in specialized metabolism appear to be particularly affected. Although there is ample evidence of these effects at the genetic level, a reduction in diversity at the metabolite level has been taken for granted despite having never been adequately accessed and quantified. Here we leveraged the high coverage of ultra high performance liquid chromatography-high-resolution mass spectrometry based metabolomics to investigate the metabolic diversity in the common bean (Phaseolus vulgaris). Information theory highlights a shift towards lower metabolic diversity and specialization when comparing wild and domesticated bean accessions. Moreover, molecular networking approaches facilitated a broader metabolite annotation than achieved to date, and its integration with gene expression data uncovers a metabolic shift from specialized metabolism towards central metabolism upon domestication of this crop.

A comprehensive metabolomics and lipidomics atlas for the legumes common bean, chickpea, lentil and lupin.

Legumes represent an important component of human and livestock diets; they are rich in macro- and micronutrients such as proteins, dietary fibers and polyunsaturated fatty acids. Whilst several health-promoting and anti-nutritional properties have been associated with grain content, in-depth metabolomics characterization of major legume species remains elusive. In this article, we used both gas chromatography-mass spectrometry (GC-MS) and liquid chromatography-mass spectrometry (LC-MS) to assess the metabolic diversity in the five legume species commonly grown in Europe, including common bean (Phaseolus vulgaris), chickpea (Cicer arietinum), lentil (Lens culinaris), white lupin (Lupinus albus) and pearl lupin (Lupinus mutabilis), at the tissue level. We were able to detect and quantify over 3400 metabolites covering major nutritional and anti-nutritional compounds. Specifically, the metabolomics atlas includes 224 derivatized metabolites, 2283 specialized metabolites and 923 lipids. The data generated here will serve the community as a basis for future integration to metabolomics-assisted crop breeding and facilitate metabolite-based genome-wide association studies to dissect the genetic and biochemical bases of metabolism in legume species.

Communication in oncology between healthcare providers, patients, the scientific community, and the media: recommendations from the Italian Association of Medical Oncology (AIOM).

AIM: To identify barriers between health and communication in oncology in order to promote the best possible practice. The areas of communication to be focused on are communication directly with the patient, communication within the scientific community, and communication with the media. MATERIAL AND METHODS: A working group including eminent experts from the national mass media, healthcare system, and patients' advocacy has been established on behalf of the Italian Association of Medical Oncology (AIOM), with the aim of developing suitable recommendations for the best communication in oncology. A literature search has been conducted selecting primary studies related to the best practices applied to communication in oncology. Subsequent to having identified the most representative statements, through a consensus conference using the RAND/University of California Los Angeles modified Delphi method, the panel evaluated the relevance of each statement to support useful strategies to develop effective communication between oncologist physicians and patients, communication within the scientific community, and communication with media outlets, including social media. RESULTS: A total of 292 statements have been extracted from 100 articles. Following an evaluation of relevance, it was found that among the 142 statements achieving the highest scores, 30 of these have been considered of particular interest by the panel. CONCLUSIONS: This consensus and the arising document represent an attempt to strengthen the strategic alliance between key figures in communication, identifying high-impact recommendations for the management of communication in oncology with respect to patients, the wider scientific community, and the media.

Genotype Combinations Drive Variability in the Microbiome Configuration of the Rhizosphere of Maize/Bean Intercropping System.

In an intercropping system, the interplay between cereals and legumes, which is strongly driven by the complementarity of below-ground structures and their interactions with the soil microbiome, raises a fundamental query: Can different genotypes alter the configuration of the rhizosphere microbial communities? To address this issue, we conducted a field study, probing the effects of intercropping and diverse maize (Zea mays L.) and bean (Phaseolus vulgaris L., Phaseolus coccineus L.) genotype combinations. Through amplicon sequencing of bacterial 16S rRNA genes from rhizosphere samples, our results unveil that the intercropping condition alters the rhizosphere bacterial communities, but that the degree of this impact is substantially affected by specific genotype combinations. Overall, intercropping allows the recruitment of exclusive bacterial species and enhances community complexity. Nevertheless, combinations of maize and bean genotypes determine two distinct groups characterized by higher or lower bacterial community diversity and complexity, which are influenced by the specific bean line associated. Moreover, intercropped maize lines exhibit varying propensities in recruiting bacterial members with more responsive lines showing preferential interactions with specific microorganisms. Our study conclusively shows that genotype has an impact on the rhizosphere microbiome and that a careful selection of genotype combinations for both species involved is essential to achieve compatibility optimization in intercropping.

Genetic control of pod morphological traits and pod edibility in a common bean RIL population.

KEY MESSAGE: QTL mapping, association analysis, and colocation study with previously reported QTL revealed three main regions controlling pod morphological traits and two loci for edible pod characteristics on the common bean chromosomes Pv01 and Pv06. Bean pod phenotype is a complex characteristic defined by the combination of different traits that determine the potential use of a genotype as a snap bean. In this study, the TUM RIL population derived from a cross between 'TU' (dry) and 'Musica' (snap) was used to investigate the genetic control of pod phenotype. The character was dissected into pod morphological traits (PMTs) and edible pod characteristics (EPC). The results revealed 35 QTL for PMTs located on seven chromosomes, suggesting a strong QTL colocation on chromosomes Pv01 and Pv06. Some QTL were colocated with previously reported QTL, leading to the mapping of 15 consensus regions associated with bean PMTs. Analysis of EPC of cooked beans revealed that two major loci with epistatic effect, located on chromosomes Pv01 and Pv06, are involved in the genetic control of this trait. An association study using a subset of the Spanish Diversity Panel (snap vs. non-snap) detected 23 genomic regions, with three regions being mapped at a position similar to those of two loci identified in the TUM population. The results demonstrated the relevant roles of Pv01 and Pv06 in the modulation of bean pod phenotype. Gene ontology enrichment analysis revealed a significant overrepresentation of genes regulating the phenylpropanoid metabolic process and auxin response in regions associated with PMTs and EPC, respectively. Both biological functions converged in the lignin biosynthetic pathway, suggesting the key role of the pathway in the genetic control of bean pod phenotype.

CRISPR/Cas9-Mediated Enrichment Coupled to Nanopore Sequencing Provides a Valuable Tool for the Precise Reconstruction of Large Genomic Target Regions.

Complete and accurate identification of genetic variants associated with specific phenotypes can be challenging when there is a high level of genomic divergence between individuals in a study and the corresponding reference genome. We have applied the Cas9-mediated enrichment coupled to nanopore sequencing to perform a targeted de novo assembly and accurately reconstruct a genomic region of interest. This approach was used to reconstruct a 250-kbp target region on chromosome 5 of the common bean genome (Phaseolus vulgaris) associated with the shattering phenotype. Comparing a non-shattering cultivar (Midas) with the reference genome revealed many single-nucleotide variants and structural variants in this region. We cut five 50-kbp tiled sub-regions of Midas genomic DNA using Cas9, followed by sequencing on a MinION device and de novo assembly, generating a single contig spanning the whole 250-kbp region. This assembly increased the number of Illumina reads mapping to genes in the region, improving their genotypability for downstream analysis. The Cas9 tiling approach for target enrichment and sequencing is a valuable alternative to whole-genome sequencing for the assembly of ultra-long regions of interest, improving the accuracy of downstream genotype-phenotype association analysis.

Legume byproducts as ingredients for food applications: Preparation, nutrition, bioactivity, and techno-functional properties.

The demand for high-quality alternative food proteins has increased over the last few decades due to nutritional and environmental concerns, leading to the growing consumption of legumes such as common bean, chickpea, lentil, lupin, and pea. However, this has also increased the quantity of non-utilized byproducts (such as seed coats, pods, broken seeds, and wastewaters) that could be exploited as sources of ingredients and bioactive compounds in a circular economy. This review focuses on the incorporation of legume byproducts into foods when they are formulated as flours, protein/fiber or solid/liquid fractions, or biological extracts and uses an analytical approach to identify their nutritional, health-promoting, and techno-functional properties. Correlation-based network analysis of nutritional, technological, and sensory characteristics was used to explore the potential of legume byproducts in food products in a systematic manner. Flour is the most widely used legume-based food ingredient and is present at levels of 2%-30% in bakery products, but purified fractions and extracts should be investigated in more detail. Health beverages and vegan dressings with an extended shelf-life are promising applications thanks to the techno-functional features of legume byproducts (e.g., foaming and emulsifying behaviors) and the presence of polyphenols. A deeper exploration of eco-friendly processing techniques (e.g., fermentation and ohmic treatment) is necessary to improve the techno-functional properties of ingredients and the sensory characteristics of foods in a sustainable manner. The processing of legume byproducts combined with improved legume genetic resources could enhance the nutritional, functional, and technological properties of ingredients to ensure that legume-based foods achieve wider industrial and consumer acceptance.

The INCREASE project: Intelligent Collections of food-legume genetic resources for European agrofood systems.

Food legumes are crucial for all agriculture-related societal challenges, including climate change mitigation, agrobiodiversity conservation, sustainable agriculture, food security and human health. The transition to plant-based diets, largely based on food legumes, could present major opportunities for adaptation and mitigation, generating significant co-benefits for human health. The characterization, maintenance and exploitation of food-legume genetic resources, to date largely unexploited, form the core development of both sustainable agriculture and a healthy food system. INCREASE will implement, on chickpea (Cicer arietinum), common bean (Phaseolus vulgaris), lentil (Lens culinaris) and lupin (Lupinus albus and L. mutabilis), a new approach to conserve, manage and characterize genetic resources. Intelligent Collections, consisting of nested core collections composed of single-seed descent-purified accessions (i.e., inbred lines), will be developed, exploiting germplasm available both from genebanks and on-farm and subjected to different levels of genotypic and phenotypic characterization. Phenotyping and gene discovery activities will meet, via a participatory approach, the needs of various actors, including breeders, scientists, farmers and agri-food and non-food industries, exploiting also the power of massive metabolomics and transcriptomics and of artificial intelligence and smart tools. Moreover, INCREASE will test, with a citizen science experiment, an innovative system of conservation and use of genetic resources based on a decentralized approach for data management and dynamic conservation. By promoting the use of food legumes, improving their quality, adaptation and yield and boosting the competitiveness of the agriculture and food sector, the INCREASE strategy will have a major impact on economy and society and represents a case study of integrative and participatory approaches towards conservation and exploitation of crop genetic resources.

A common bean truncated CRINKLY4 kinase controls gene-for-gene resistance to the fungus Colletotrichum lindemuthianum.

Identifying the molecular basis of resistance to pathogens is critical to promote a chemical-free cropping system. In plants, nucleotide-binding leucine-rich repeat constitute the largest family of disease resistance (R) genes, but this resistance can be rapidly overcome by the pathogen, prompting research into alternative sources of resistance. Anthracnose, caused by the fungus Colletotrichum lindemuthianum, is one of the most important diseases of common bean. This study aimed to identify the molecular basis of Co-x, an anthracnose R gene conferring total resistance to the extremely virulent C. lindemuthianum strain 100. To that end, we sequenced the Co-x 58 kb target region in the resistant JaloEEP558 (Co-x) common bean and identified KTR2/3, an additional gene encoding a truncated and chimeric CRINKLY4 kinase, located within a CRINKLY4 kinase cluster. The presence of KTR2/3 is strictly correlated with resistance to strain 100 in a diversity panel of common beans. Furthermore, KTR2/3 expression is up-regulated 24 hours post-inoculation and its transient expression in a susceptible genotype increases resistance to strain 100. Our results provide evidence that Co-x encodes a truncated and chimeric CRINKLY4 kinase probably resulting from an unequal recombination event that occurred recently in the Andean domesticated gene pool. This atypical R gene may act as a decoy involved in indirect recognition of a fungal effector.

Pod indehiscence in common bean is associated with the fine regulation of PvMYB26.

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.

Comparative Analysis Based on Transcriptomics and Metabolomics Data Reveal Differences between Emmer and Durum Wheat in Response to Nitrogen Starvation.

Mounting evidence indicates the key role of nitrogen (N) on diverse processes in plant, including development and defense. Using a combined transcriptomics and metabolomics approach, we studied the response of seedlings to N starvation of two different tetraploid wheat genotypes from the two main domesticated subspecies: emmer and durum wheat. We found that durum wheat exhibits broader and stronger response in comparison to emmer as seen from the expression pattern of both genes and metabolites and gene enrichment analysis. They showed major differences in the responses to N starvation for transcription factor families, emmer showed differential reduction in the levels of primary metabolites while durum wheat exhibited increased levels of most of them to N starvation. The correlation-based networks, including the differentially expressed genes and metabolites, revealed tighter regulation of metabolism in durum wheat in comparison to emmer. We also found that glutamate and γ-aminobutyric acid (GABA) had highest values of centrality in the metabolic correlation network, suggesting their critical role in the genotype-specific response to N starvation of emmer and durum wheat, respectively. Moreover, this finding indicates that there might be contrasting strategies associated to GABA and glutamate signaling modulating shoot vs. root growth in the two different wheat subspecies.

Multi-tissue integration of transcriptomic and specialized metabolite profiling provides tools for assessing the common bean (Phaseolus vulgaris) metabolome.

Common bean (Phaseolus vulgaris L.) is an important legume species with a rich natural diversity of landraces that originated from the wild forms following multiple independent domestication events. After the publication of its genome, several resources for this relevant crop have been made available. A comprehensive characterization of specialized metabolism in P. vulgaris, however, is still lacking. In this study, we used a metabolomics approach based on liquid chromatography-mass spectrometry to dissect the chemical composition at a tissue-specific level in several accessions of common bean belonging to different gene pools. Using a combination of literature search, mass spectral interpretation, 13 C-labeling, and correlation analyses, we were able to assign chemical classes and/or putative structures for approximately 39% of all measured metabolites. Additionally, we integrated this information with transcriptomics data and phylogenetic inference from multiple legume species to reconstruct the possible metabolic pathways and identify sets of candidate genes involved in the biosynthesis of specialized metabolites. A particular focus was given to flavonoids, triterpenoid saponins and hydroxycinnamates, as they represent metabolites involved in important ecological interactions and they are also associated with several health-promoting benefits when integrated into the human diet. The data are presented here in the form of an accessible resource that we hope will set grounds for further studies on specialized metabolism in legumes.

Genomic dissection of pod shattering in common bean: mutations at non-orthologous loci at the basis of convergent phenotypic evolution under domestication of leguminous species.

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.

Genome-Wide Association Mapping of Prostrate/Erect Growth Habit in Winter Durum Wheat.

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.

Adapting legume crops to climate change using genomic approaches.

Our agricultural system and hence food security is threatened by combination of events, such as increasing population, the impacts of climate change, and the need to a more sustainable development. Evolutionary adaptation may help some species to overcome environmental changes through new selection pressures driven by climate change. However, success of evolutionary adaptation is dependent on various factors, one of which is the extent of genetic variation available within species. Genomic approaches provide an exceptional opportunity to identify genetic variation that can be employed in crop improvement programs. In this review, we illustrate some of the routinely used genomics-based methods as well as recent breakthroughs, which facilitate assessment of genetic variation and discovery of adaptive genes in legumes. Although additional information is needed, the current utility of selection tools indicate a robust ability to utilize existing variation among legumes to address the challenges of climate uncertainty.

An Early Tailored Approach Is the Key to Effective Rehabilitation in the Intensive Care Unit.

OBJECTIVE: To investigate the effectiveness, feasibility, and safety of an evidence-based rehabilitation care pathway in the intensive care unit (ICU) in different patient populations. DESIGN: Observational prospective cohort study, with retrospective controls. SETTING: ICUs of a university hospital. PARTICIPANTS: Patients admitted between April 1, 2015, and June 30, 2015, were compared to a retrospective cohort admitted to the same ICUs during the same 3-month period in 2014. The number of patients studied (N=285) included 152 in the prospective group and 133 in the retrospective group. INTERVENTIONS: The prospective cohort benefited of a rehabilitation care pathway based on (1) interdisciplinary teamwork; (2) early customized and goal-oriented rehabilitation; (3) daily functional monitoring and treatment revision; (4) agreed discharge policy; and (5) continuity of care. The retrospective cohort underwent usual care. MAIN OUTCOME MEASURES: Included the following: (1) proportions of patients undergoing rehabilitation team evaluation; (2) latency between patient admission to ICUs and rehabilitation team assessment; (3) proportions of patients undergoing rehabilitation treatment during ICU stay; (4) latency between the patient admission to ICUs and rehabilitation start; (5) ICU stay and total acute hospital stay; and (5) proportion of ventilator-free days out of ICU stay. RESULTS: The novel rehabilitation care pathway led to (1) an increased proportion of patients receiving rehabilitative assessment (P<.0001); (2) a decreased latency from ICU admission to both rehabilitation team assessment and rehabilitation start (P<.0001); (3) an increased proportion of patients undergoing rehabilitation (P<.0001); (4) a shorter length of stay in ICUs (P<.0001) and in hospital (P=.047); and (5) a shorter mechanical ventilation duration (P<.02). A direct relationship between rehabilitation start latency and ICU length of stay was observed. CONCLUSIONS: An early, interdisciplinary team approach, providing a customized dynamic planning of physiotherapy programs, increases ventilator-free time and reduces total hospital stay, especially in patients admitted to the ICU after general surgery. This rehabilitation care pathway can be generalized to different geopolitical scenarios, being feasible, safe and cost effective.

Evolutionary Metabolomics Reveals Domestication-Associated Changes in Tetraploid Wheat Kernels.

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.

Decreased Nucleotide and Expression Diversity and Modified Coexpression Patterns Characterize Domestication in the Common Bean.

Using RNA sequencing technology and de novo transcriptome assembly, we compared representative sets of wild and domesticated accessions of common bean (Phaseolus vulgaris) from Mesoamerica. RNA was extracted at the first true-leaf stage, and de novo assembly was used to develop a reference transcriptome; the final data set consists of ∼190,000 single nucleotide polymorphisms from 27,243 contigs in expressed genomic regions. A drastic reduction in nucleotide diversity (∼60%) is evident for the domesticated form, compared with the wild form, and almost 50% of the contigs that are polymorphic were brought to fixation by domestication. In parallel, the effects of domestication decreased the diversity of gene expression (18%). While the coexpression networks for the wild and domesticated accessions demonstrate similar seminal network properties, they show distinct community structures that are enriched for different molecular functions. After simulating the demographic dynamics during domestication, we found that 9% of the genes were actively selected during domestication. We also show that selection induced a further reduction in the diversity of gene expression (26%) and was associated with 5-fold enrichment of differentially expressed genes. While there is substantial evidence of positive selection associated with domestication, in a few cases, this selection has increased the nucleotide diversity in the domesticated pool at target loci associated with abiotic stress responses, flowering time, and morphology.