Cis- and Trans-Regulatory Variations in the Domestication of the Chili Pepper Fruit

Intraspecific gene regulation in cis- and trans

Changes in gene expression underlie much of evolution and occur via either cis-acting mutations, which lie near the affected gene and act in a context-specific manner, or trans-acting mutations, which may be far from the affected gene and act through diffusible molecules such as transcription factors. A commonly held view is that most expression variation within species is controlled in trans- while expression differences between species are largely controlled in cis-. Here, we summarize recent intraspecific gene regulation studies and find, contrary to this widely held view, that many studies in diverse taxa have revealed a large role for cis-acting mutations underlying expression variation within species. A review of the existing literature also shows that preparations using whole organisms rather than individual tissues may be biased toward identifying trans-regulation. Moreover, we note several examples of predominantly cis-acting regulation in recently diverged populations adapted to different environments. We highlight the challenges of drawing general conclusions from comparisons among studies that use different methodologies and we offer suggestions for studies that will address outstanding questions concerning the evolution of gene regulation.

Subgenome-informed statistical modeling of transcriptomes in 25 common wheat accessions reveals cis- and trans-regulation architectures

Common wheat is an allohexaploid, where it is difficult to obtain homoeolog-distinguished transcriptome data. Lasy-Seq, a type of 3' RNA-seq, is technology efficient at obtaining homoeolog-distinguished transcriptomes. Here, we applied Lasy-Seq to obtain transcriptome data from the seedlings, second leaves, and root tips of 25 common wheat lines mainly from East Asia. Roots and seedlings were similar to each other in transcriptome profiles, but they were different from the leaves. We then asked how three homoeologous genes from different subgenomes (i.e. triads) show different levels of expression. Specifically, we examined the effects of subgenomes, lines, and their interaction on the expression levels of each homoeolog triad, separately in each tissue. Of the 19 805 homoeolog triads, 51-55% showed significant effect of subgenome, suggesting cis-regulation, whereas 24-30% showed significant effect line, suggesting trans-regulation. We also found that 7.7-9.0% triads showed significant effects of the interaction. Hierarchical clustering and co-trans regulation network analysis of homoeolog triads revealed that the patterns of expression polymorphisms among the lines were shared in different genes. Our results also implied that expression variation between lines is caused by changes in a smaller number of common trans-factors. We performed gene ontology (GO)-term enrichment analysis using newly annotated and substantially improved GO annotations, which revealed that GO terms related to each tissue-type function were enriched in genes expressed in the leaves and roots. Our information provides fundamental knowledge for the future breeding of plants possessing complex gene regulatory networks such as common wheat.

Characterization of Metabolites and Transcriptome of Pepper Accessions from Four Southern Provinces of China

BACKGROUND/OBJECTIVES

Pepper (Capsicum annuum L.) is a widely grown vegetable and spice crop worldwide. This study aims to reveal the differences of metabolites among pepper accessions from different regions and explore candidate genes related to metabolites of pepper.

METHODS

The metabolome and transcriptome of 36 pepper accessions were determined by widely targeted metabolomics analysis and RNA sequencing technology, and the differential metabolites and differential genes among C. annuum from four important pepper production and consumption provinces of China, Hunan, Guizhou, Yunnan and Sichuan, were analyzed.

RESULTS

Flavonoids are the main characteristic metabolites that distinguish pepper accessions from Yun_Gui_Chuan Group and Hunan Group. The aglycones of characteristic flavonoids in each group are different; in Yun_Gui_Chuan Group mainly are luteolin, quercetin, chrysoeriol and isorhamnetin; in Hunan Group mainly are apigenin. Transcriptome data showed that two genes related to flavonoid 3'-monooxygenase differed significantly between the two groups of chili peppers, and we speculated that they may be the core enzymes regulating their flavonoid profile. And an SNP mutation located in gene Cgla06g001871 showed a strong correlation with pepper accessions from Yun_Gui_Chuan Group, which can be used as a DNA marker to identify pepper accessions from Yun_Gui_Chuan Group, and provide strong support for regional specialty variety conservation. In addition, we also analyzed the metabolites related to the taste and nutrition of pepper accessions in the four provinces, and the results showed that the sugar content of pepper accessions from Guizhou was low and the capsaicinoids content of pepper accessions from Sichuan was low, while no significant difference was found in acid and vitamin contents among pepper accessions from the four provinces.

CONCLUSIONS

The metabolome and transcriptome of 36 pepper accessions from four important pepper production and consumption provinces of China were determined, and the characteristic metabolites and expressed genes of pepper accessions from each province were analyzed.

Genetic Regulatory Perturbation of Gene Expression Impacted by Genomic Introgression in Fiber Development of Allotetraploid Cotton

Interspecific genomic introgression is an important evolutionary process with respect to the generation of novel phenotypic diversity and adaptation. A key question is how gene flow perturbs gene expression networks and regulatory interactions. Here, an introgression population of two species of allopolyploid cotton (Gossypium) to delineate the regulatory perturbations of gene expression regarding fiber development accompanying fiber quality change is utilized. De novo assembly of the recipient parent (G. hirsutum Emian22) genome allowed the identification of genomic variation and introgression segments (ISs) in 323 introgression lines (ILs) from the donor parent (G. barbadense 3-79). It documented gene expression dynamics by sequencing 1,284 transcriptomes of developing fibers and characterized genetic regulatory perturbations mediated by genomic introgression using a multi-locus model. Introgression of individual homoeologous genes exhibiting extreme low or high expression bias can lead to a parallel expression bias in their non-introgressed duplicates, implying a shared yet divergent regulatory fate of duplicated genes following allopolyploidy. Additionally, the IL N182 with improved fiber quality is characterized, and the candidate gene GhFLAP1 related to fiber length is validated. This study outlines a framework for understanding introgression-mediated regulatory perturbations in polyploids, and provides insights for targeted breeding of superior upland cotton fiber.

Anthocyanin biosynthesis in goji berry is inactivated by deletion in a bHLH transcription factor LrLAN1b promoter

Black goji berry (Lycium ruthenicum Murray) contains a rich source of health-promoting anthocyanins which are used in herbal medicine and nutraceutical foods in China. A natural variant producing white berries allowed us to identify two key genes involved in the regulation of anthocyanin biosynthesis in goji berries: one encoding a MYB transcription factor (LrAN2-like) and one encoding a basic helix-loop-helix (bHLH) transcription factor (LrAN1b). We previously found that LrAN1b expression was lost in the white berry variant, but the molecular basis for this phenotype was unknown. Here, we identified the molecular mechanism for loss of anthocyanins in white goji berries. In white goji, the LrAN1b promoter region has a 229 bp deletion that removes three MYB-binding elements and one bHLH-binding element, which are key to its expression. Complementation of the white goji berry LrAN1b allele with the LrAN1b promoter restored pigmentation. Virus-induced gene silencing of LrAN1b in black goji berry reduced fruit anthocyanin biosynthesis. Molecular analyses showed that LrAN2-like and another bHLH transcription factor LrJAF13 can activate LrAN1b by binding directly to the MYB-recognizing element and bHLH-recognizing element of its promoter-deletion region. LrAN1b expression is enhanced by the interaction of LrAN2-like with LrJAF13 and the WD40 protein LrAN11. LrAN2-like and LrAN11 interact with either LrJAF13 or LrAN1b to form two MYB-bHLH-WD40 complexes, which hierarchically regulate anthocyanin biosynthesis in black goji berry. This study on a natural variant builds a comprehensive anthocyanin regulatory network that may be manipulated to tailor goji berry traits.

Genetic analysis and QTL mapping of domestication-related traits in chili pepper (Capsicum annuum L.)

Chili pepper (Capsicum annuum L.) is one of the oldest and most phenotypically diverse pre-Columbian crops of the Americas. Despite the abundance of genetic resources, the use of wild germplasm and landraces in chili pepper breeding is limited. A better understanding of the evolutionary history in chili peppers, particularly in the context of traits of agronomic interest, can contribute to future improvement and conservation of genetic resources. In this study, an F₂:₃ mapping population derived from a cross between a C. annuum wild accession (Chiltepin) and a cultivated variety (Puya) was used to identify genomic regions associated with 19 domestication and agronomic traits. A genetic map was constructed consisting of 1023 single nucleotide polymorphism (SNP) markers clustered into 12 linkage groups and spanning a total of 1,263.87 cM. A reciprocal translocation that differentiates the domesticated genome from its wild ancestor and other related species was identified between chromosomes 1 and 8. Quantitative trait locus (QTL) analysis detected 20 marker-trait associations for 13 phenotypes, from which 14 corresponded to previously identified loci, and six were novel genomic regions related to previously unexplored domestication-syndrome traits, including form of unripe fruit, seedlessness, deciduous fruit, and growth habit. Our results revealed that the genetic architecture of Capsicum domestication is similar to other domesticated species with few loci with large effects, the presence of QTLs clusters in different genomic regions, and the predominance of domesticated recessive alleles. Our analysis indicates the domestication process in chili pepper has also had an effect on traits not directly related to the domestication syndrome. The information obtained in this study provides a more complete understanding of the genetic basis of Capsicum domestication that can potentially guide strategies for the exploitation of wild alleles.

Dissecting the genetic basis of heterosis in elite super-hybrid rice

Y900 is one of the top hybrid rice (Oryza sativa) varieties, with its yield exceeding 15 t·hm-2. To dissect the mechanism of heterosis, we sequenced the male parent line R900 and female parent line Y58S using long-read and Hi-C technology. High-quality reference genomes of 396.41 Mb and 398.24 Mb were obtained for R900 and Y58S, respectively. Genome-wide variations between the parents were systematically identified, including 1,367,758 single-nucleotide polymorphisms, 299,149 insertions/deletions, and 4,757 structural variations. The level of variation between Y58S and R900 was the lowest among the comparisons of Y58S with other rice genomes. More than 75% of genes exhibited variation between the two parents. Compared with other two-line hybrids sharing the same female parent, the portion of Geng/japonica (GJ)-type genetic components from different male parents increased with yield increasing in their corresponding hybrids. Transcriptome analysis revealed that the partial dominance effect was the main genetic effect that constituted the heterosis of Y900. In the hybrid, both alleles from the two parents were expressed, and their expression patterns were dynamically regulated in different tissues. The cis-regulation was dominant for young panicle tissues, while trans-regulation was more common in leaf tissues. Overdominance was surprisingly prevalent in stems and more likely regulated by the trans-regulation mechanism. Additionally, R900 contained many excellent GJ haplotypes, such as NARROW LEAF1, Oryza sativa SQUAMOSA PROMOTER BINDING PROTEIN-LIKE13, and Grain number, plant height, and heading date8, making it a good complement to Y58S. The fine-tuned mechanism of heterosis involves genome-wide variation, GJ introgression, key functional genes, and dynamic gene/allele expression and regulation pattern changes in different tissues and growth stages.

Expression inheritance and constraints on cis- and trans-regulatory mutations underlying lotus color variation

Both cis- and trans-regulatory mutations drive changes in gene expression that underpin plant phenotypic evolution. However, how and why these two major types of regulatory mutations arise in different genes and how gene expression is inherited and associated with these regulatory changes are unclear. Here, by studying allele-specific expression in F1 hybrids of pink-flowered sacred lotus (Nelumbo nucifera) and yellow-flowered American lotus (N. lutea), we reveal the relative contributions of cis- and trans-regulatory changes to interspecific expression rewiring underlying petal color change and how the expression is inherited in hybrids. Although cis-only variants influenced slightly more genes, trans-only variants had a stronger impact on expression differences between species. In F1 hybrids, genes under cis-only and trans-only regulatory effects showed a propensity toward additive and dominant inheritance, respectively, whereas transgressive inheritance was observed in genes carrying both cis- and trans-variants acting in opposite directions. By investigating anthocyanin and carotenoid coexpression networks in petals, we found that the same category of regulatory mutations, particularly trans-variants, tend to rewire hub genes in coexpression modules underpinning flower color differentiation between species; we identified 45 known genes with cis- and trans-regulatory variants significantly correlated with flower coloration, such as ANTHOCYANIN 5-AROMATIC ACYLTRANSFERASE (ACT), GLUTATHIONE S-TRANSFERASE F11 (GSTF11), and LYCOPENE Ε-CYCLASE (LCYE). Notably, the relative abundance of genes in different categories of regulatory divergence was associated with the inferred magnitude of constraints like expression level and breadth. Overall, our study suggests distinct selective constraints and modes of gene expression inheritance among different regulatory mutations underlying lotus petal color divergence.

The role of non-additive gene action on gene expression variation in plant domestication

BACKGROUND

Plant domestication is a remarkable example of rapid phenotypic transformation of polygenic traits, such as organ size. Evidence from a handful of study cases suggests this transformation is due to gene regulatory changes that result in non-additive phenotypes. Employing data from published genetic crosses, we estimated the role of non-additive gene action in the modulation of transcriptional landscapes in three domesticated plants: maize, sunflower, and chili pepper. Using A. thaliana, we assessed the correlation between gene regulatory network (GRN) connectivity properties, transcript abundance variation, and gene action. Finally, we investigated the propagation of non-additive gene action in GRNs.

RESULTS

We compared crosses between domesticated plants and their wild relatives to a set of control crosses that included a pair of subspecies evolving under natural selection and a set of inbred lines evolving under domestication. We found abundance differences on a higher portion of transcripts in crosses between domesticated-wild plants relative to the control crosses. These transcripts showed non-additive gene action more often in crosses of domesticated-wild plants than in our control crosses. This pattern was strong for genes associated with cell cycle and cell fate determination, which control organ size. We found weak but significant negative correlations between the number of targets of trans-acting genes (Out-degree) and both the magnitude of transcript abundance difference a well as the absolute degree of dominance. Likewise, we found that the number of regulators that control a gene's expression (In-degree) is weakly but negatively correlated with the magnitude of transcript abundance differences. We observed that dominant-recessive gene action is highly propagable through GRNs. Finally, we found that transgressive gene action is driven by trans-acting regulators showing additive gene action.

CONCLUSIONS

Our study highlights the role of non-additive gene action on modulating domestication-related traits, such as organ size via regulatory divergence. We propose that GRNs are shaped by regulatory changes at genes with modest connectivity, which reduces the effects of antagonistic pleiotropy. Finally, we provide empirical evidence of the propagation of non-additive gene action in GRNs, which suggests a transcriptional epistatic model for the control of polygenic traits, such as organ size.

Combined population transcriptomic and genomic analysis reveals cis-regulatory differentiation of non-coding RNAs in maize

Long intergenic non-coding RNA (lincRNA), cis-acting expression quantitative trait locus (cis-eQTL), maize, regulatory evolution. The law of genetic variation during domestication explains the evolutionary mechanism and provides a theoretical basis for improving existing varieties of maize. Previous studies focused on exploiting regulatory variations controlling the expression of protein-coding genes rather than of non-protein-coding genes. Here, we examined the genetic and evolutionary features of long non-coding RNAs from intergenic regions (long intergenic non-coding RNAs, lincRNAs) using population-scale transcriptome data and identified 1168 lincRNAs with cis-acting expression quantitative trait loci (cis-eQTLs). We found that lincRNAs are more likely to be regulated by cis-eQTLs, which exert stronger effects than the protein-coding genes. During maize domestication and improvement, upregulated alleles of lincRNAs, which originated from both standing variation and new mutation, accumulate more frequently and show larger effect sizes than the coding genes. A stronger signature of genetic differentiation was observed in their regulatory regions compared to those of randomly sampled lincRNAs. In addition, we found that cis-regulatory differentiation of lincRNAs is related to the sequence conservation of lincRNA transcripts. Non-conserved lincRNAs more tend to gain upregulated alleles and show a stronger relationship with selected traits than conserved lincRNAs between maize and its wild relatives. Our findings in maize improve the understanding of cis-regulatory variation in lincRNA genes during domestication and improvement and provide an effective approach for prioritizing candidates for further investigation.

Spatiotemporal Regulation of a Single Adaptively Evolving Trans-Regulatory Element Contributes to Spermatogenetic Expression Divergence in Drosophila

Due to extensive pleiotropy, trans-acting elements are often thought to be evolutionarily constrained. While the impact of trans-acting elements on gene expression evolution has been extensively studied, relatively little is understood about the contribution of a single trans regulator to interspecific expression and phenotypic divergence. Here, we disentangle the effects of genomic context and miR-983, an adaptively evolving young microRNA, on expression divergence between Drosophila melanogaster and D. simulans. We show miR-983 effects promote interspecific expression divergence in testis despite its antagonism with the often-predominant context effects. Single-cyst RNA-seq reveals that distinct sets of genes gain and lose miR-983 influence under disruptive or diversifying selection at different stages of spermatogenesis, potentially helping minimize antagonistic pleiotropy. At the round spermatid stage, the effects of miR-983 are weak and distributed, coincident with the transcriptome undergoing drastic expression changes. Knocking out miR-983 causes reduced sperm length with increased within-individual variation in D. melanogaster but not in D. simulans, and the D. melanogaster knockout also exhibits compromised sperm defense ability. Our results provide empirical evidence for the resolution of antagonistic pleiotropy and also have broad implications for the function and evolution of new trans regulators.

Transcriptome shock in interspecific F1 allotriploid hybrids between Brassica species

Interspecific hybridization drives the evolution of angiosperms and can be used to introduce novel alleles for important traits or to activate heterosis in crop breeding. Hybridization brings together gene expression networks from two different species, potentially causing global alterations of gene expression in the F1 plants which is called 'transcriptome shock'. Here, we explored such a transcriptome shock in allotriploid Brassica hybrids. We generated interspecific F1 allotriploid hybrids between the allotetraploid species Brassica napus and three accessions of the diploid species Brassica rapa. RNA-seq of the F1 hybrids and the parental plants revealed that 26.34-30.89% of genes were differentially expressed between the parents. We also analyzed expression level dominance and homoeolog expression bias between the parents and the F1 hybrids. The expression-level dominance biases of the Ar, An, and Cn subgenomes was genotype and stage dependent, whereas significant homoeolog expression bias was observed among three subgenomes from different parents. Furthermore, more genes were involved in trans regulation than in cis regulation in allotriploid F1 hybrids. Our findings provide new insights into the transcriptomic responses of cross-species hybrids and hybrids showing heterosis, as well as a new method for promoting the breeding of desirable traits in polyploid Brassica species.

Transcriptomic Variations and Network Hubs Controlling Seed Size and Weight During Maize Seed Development

To elucidate the mechanisms underlying seed development in maize, comprehensive RNA-seq analyses were conducted on Zhengdan1002 (ZD1002), Zhengdan958 (ZD958), and their parental lines during seven seed developmental stages. We found that gene expression levels were largely nonadditive in hybrids and that cis-only or trans × cis pattern played a large role in hybrid gene regulation during seed developmental stage. Weighted gene co-expression network (WGCNA) analysis showed that 36 modules were highly correlated (r = -0.90-0.92, p < 0.05) with kernel weight, length, and width during seed development. Forty-five transcription factors and 38 ribosomal protein genes were identified as major hub genes determining seed size/weight. We also described a network hub, Auxin Response Factor 12 of maize (ZmARF12), a member of a family of transcription factor that mediate gene expression in response to auxin, potentially links auxin signal pathways, cell division, and the size of the seeds. The ZmARF12 mutant exhibited larger seed size and higher grain weight. ZmARF12 transcription was negatively associated with cell division during seed development, which was confirmed by evaluating the yield of protoplasts that isolated from the kernels of the mutant and other inbred lines. Transient knock-down of ZmARF12 in maize plants facilitated cell expansion and division, whereas transient silencing of its potential interactor ZmIAA8 impaired cell division. ZmIAA8 expression was repressed in the ZmARF12 over-expressed protoplasts. The mutant phenotype and the genetics studies presented here illustrated evidence that ZmARF12 is a cell division repressor, and potentially determines the final seed size.

Inheritance of gene expression throughout fruit development in chili pepper

Gene expression is the primary molecular phenotype and can be estimated in specific organs or tissues at particular times. Here we analyzed genome-wide inheritance of gene expression in fruits of chili pepper (Capsicum annuum L.) in reciprocal crosses between a domesticated and a wild accession, estimating this parameter during fruit development. We defined a general hierarchical schema to classify gene expression inheritance which can be employed for any quantitative trait. We found that inheritance of gene expression is affected by both, the time of fruit development as well as the direction of the cross, and propose that such variations could be common in many developmental processes. We conclude that classification of inheritance patterns is important to have a better understanding of the mechanisms underlying gene expression regulation, and demonstrate that sets of genes with specific inheritance pattern at particular times of fruit development are enriched in different biological processes, molecular functions and cell components. All curated data and functions for analysis and visualization are publicly available as an R package.

Extracellular DNA: Insight of a Signal Molecule in Crop Protection

Simple Summary

Agriculture systems use multiple chemical treatments to prevent pests and diseases, and to fertilize plants and eliminate weeds around the crop. These practices are less accepted by the consumers each day, mostly because of the associated environmental, health, and ecological impact; thus, new sustainable green technologies are being developed to replace the use of chemical products. Among green technologies for agriculture practices, the use of plant elicitors represents an alternative with great potential, and extracellular DNA has shown beneficial effects on important production traits such as defence mechanisms, plant growth and development, and secondary metabolites production that results in yield increment and better-quality food. In this review, we reunite experimental evidence of the natural effect that extracellular DNA has on plants. We also aim to contribute a step closer to the agricultural application of extracellular DNA. Additionally, we suggest that extracellular DNA can have a biostimulant effect on plants, and can be applied as a highly sustainable treatment contributing to the circular economy of primary production.

Abstract

Agricultural systems face several challenges in terms of meeting everyday-growing quantities and qualities of food requirements. However, the ecological and social trade-offs for increasing agricultural production are high, therefore, more sustainable agricultural practices are desired. Researchers are currently working on diverse sustainable techniques based mostly on natural mechanisms that plants have developed along with their evolution. Here, we discuss the potential agricultural application of extracellular DNA (eDNA), its multiple functioning mechanisms in plant metabolism, the importance of hormetic curves establishment, and as a challenge: the technical limitations of the industrial scale for this technology. We highlight the more viable natural mechanisms in which eDNA affects plant metabolism, acting as a damage/microbe-associated molecular pattern (DAMP, MAMP) or as a general plant biostimulant. Finally, we suggest a whole sustainable system, where DNA is extracted from organic sources by a simple methodology to fulfill the molecular characteristics needed to be applied in crop production systems, allowing the reduction in, or perhaps the total removal of, chemical pesticides, fertilizers, and insecticides application.

Domestication of Crop Metabolomes: Desired and Unintended Consequences

The majority of the crops and vegetables of today were domesticated from their wild progenitors within the past 12 000 years. Considerable research effort has been expended on characterizing the genes undergoing positive and negative selection during the processes of crop domestication and improvement. Many studies have also documented how the contents of a handful of metabolites have been altered during human selection, but we are only beginning to unravel the true extent of the metabolic consequences of breeding. We highlight how crop metabolomes have been wittingly or unwittingly shaped by the processes of domestication, and highlight how we can identify new targets for metabolite engineering for the purpose of de novo domestication of crop wild relatives.

Biased allelic expression in tissues of F1 hybrids between tropical and temperate lotus (Nelumbo nuicfera)

The genome-wide allele-specific expression in F1 hybrids from the cross of tropical and temperate lotus unveils how cis-regulatory divergences affect genes in key pathways related to ecotypic divergence. Genetic variation, particularly cis-regulatory variation, plays a crucial role in phenotypic variation and adaptive evolution in plants. Temperate and tropical lotus, the two ecotypes of Nelumbo nucifera, show distinction in the degree of rhizome enlargement, which is associated with winter dormancy. To understand the roles of genome-wide cis-regulatory divergences on adaptive evolution of temperate and tropical lotus (Nelumbo nucifera), here we performed allele-specific expression (ASE) analyses on the tissues including flowers, leaves and rhizome from F1 hybrids of tropical and temperate lotus. For all investigated tissues in F1s, about 36% of genes showed ASE and about 3% of genes showed strong consistent ASE. Most of ASEs were biased towards the tropical parent in all surveyed samples, indicating that the tropical genome might be dominant over the temperate genome in gene expression of tissues from their F1 hybrids. We found that promoter sequences with similar allelic expression are more conserved than genes with significant or conditional ASE, suggesting the cis-regulatory sequence divergence underlie the allelic expression bias. We further uncovered biased genes being related to phenotypic differentiation between two lotus ecotypes, especially metabolic and phytohormone-related pathways in the rhizome. Overall, our study provides a global landscape of cis-regulatory variations between two lotus ecotypes and highlights their roles in rhizome growth variation for the climatic adaptation.

Transcriptome Analyses Throughout Chili Pepper Fruit Development Reveal Novel Insights into the Domestication Process

Chili pepper (Capsicum spp.) is an important crop, as well as a model for fruit development studies and domestication. Here, we performed a time-course experiment to estimate standardized gene expression profiles with respect to fruit development for six domesticated and four wild chili pepper ancestors. We sampled the transcriptomes every 10 days from flowering to fruit maturity, and found that the mean standardized expression profiles for domesticated and wild accessions significantly differed. The mean standardized expression was higher and peaked earlier for domesticated vs. wild genotypes, particularly for genes involved in the cell cycle that ultimately control fruit size. We postulate that these gene expression changes are driven by selection pressures during domestication and show a robust network of cell cycle genes with a time shift in expression, which explains some of the differences between domesticated and wild phenotypes.

Genomic, Transcriptomic and Epigenomic Tools to Study the Domestication of Plants and Animals: A Field Guide for Beginners

In the last decade, genomics and the related fields of transcriptomics and epigenomics have revolutionized the study of the domestication process in plants and animals, leading to new discoveries and new unresolved questions. Given that some domesticated taxa have been more studied than others, the extent of genomic data can range from vast to nonexistent, depending on the domesticated taxon of interest. This review is meant as a rough guide for students and academics that want to start a domestication research project using modern genomic tools, as well as for researchers already conducting domestication studies that are interested in following a genomic approach and looking for alternate strategies (cheaper or more efficient) and future directions. We summarize the theoretical and technical background needed to carry out domestication genomics, starting from the acquisition of a reference genome and genome assembly, to the sampling design for population genomics, paleogenomics, transcriptomics, epigenomics and experimental validation of domestication-related genes. We also describe some examples of the aforementioned approaches and the relevant discoveries they made to understand the domestication of the studied taxa.