Comparative transcriptomics reveals patterns of selection in domesticated and wild tomato

Comprehensive analysis of the mitochondrial genome of Rehmannia glutinosa: insights into repeat-mediated recombinations and RNA editing-induced stop codon acquisition

Rehmannia glutinosa is an economically significant medicinal plant. Yet, the structure and sequence of its mitochondrial genome has not been published, which plays a crucial role in evolutionary analysis and regulating respiratory-related macromolecule synthesis. In this study, the R. glutinosa mitogenome was sequenced employing a combination of Illumina short reads and Nanopore long reads, with subsequent assembly using a hybrid strategy. We found that the predominant configuration of the R. glutinosa mitogenome comprises two circular chromosomes. The primary structure of the mitogenome encompasses two mitochondrial chromosomes corresponding to the two major configurations, Mac1-1 and Mac1-2. The R. glutinosa mitogenome encoded an angiosperm-typical set of 24 core genes, nine variable genes, three rRNA genes, and 15 tRNA genes. A phylogenetic analysis using the 16 shared protein-coding genes (PCG) yielded a tree consistent with the phylogeny of Lamiales species and two outgroup taxa. Mapping RNA-seq data to the coding sequences (CDS) of the PCGs revealed 507 C-to-U RNA editing sites across 31 PCGs of the R. glutinosa mitogenome. Furthermore, one start codon (nad4L) and two stop codons (rpl10 and atp6) were identified as products of RNA editing events in the R. glutinosa mitogenome.

Identification, Classification, and Transcriptional Analysis of Rab GTPase Genes from Tomato (Solanum lycopersicum) Reveals Salt Stress Response Genes

Salinity in plants generates an osmotic and ionic imbalance inside cells that compromises the viability of the plant. Rab GTPases, the largest family within the small GTPase superfamily, play pivotal roles as regulators of vesicular trafficking in plants, including the economically important and globally cultivated tomato (Solanum lycopersicum). Despite their significance, the specific involvement of these small GTPases in tomato vesicular trafficking and their role under saline stress remains poorly understood. In this work, we identified and classified 54 genes encoding Rab GTPases in cultivated tomato, elucidating their genomic distribution and structural characteristics. We conducted an analysis of duplication events within the S. lycopersicum genome, as well as an examination of gene structure and conserved motifs. In addition, we investigated the transcriptional profiles for these Rab GTPases in various tissues of cultivated and wild tomato species using microarray-based analysis. The results showed predominantly low expression in most of the genes in both leaves and vegetative meristem, contrasting with notably high expression levels observed in seedling roots. Also, a greater increase in gene expression in shoots from salt-tolerant wild tomato species was observed under normal conditions when comparing Solanum habrochaites, Solanum pennellii, and Solanum pimpinellifolium with S. lycopersicum. Furthermore, an expression analysis of Rab GTPases from Solanum chilense in leaves and roots under salt stress treatment were also carried out for their characterization. These findings revealed that specific Rab GTPases from the endocytic pathway and the trans-Golgi network (TGN) showed higher induction in plants exposed to saline stress conditions. Likewise, disparities in gene expression were observed both among members of the same Rab GTPase subfamily and between different subfamilies. Overall, this work emphasizes the high degree of conservation of Rab GTPases, their high functional diversification in higher plants, and the essential role in mediating salt stress tolerance and suggests their potential for further exploration of vesicular trafficking mechanisms in response to abiotic stress conditions.

Intergrative metabolomic and transcriptomic analyses reveal the potential regulatory mechanism of unique dihydroxy fatty acid biosynthesis in the seeds of an industrial oilseed crop Orychophragmus violaceus

BACKGROUND

Orychophragmus violaceus is a potentially important industrial oilseed crop due to the two 24-carbon dihydroxy fatty acids (diOH-FA) that was newly identified from its seed oil via a 'discontinuous elongation' process. Although many research efforts have focused on the diOH-FA biosynthesis mechanism and identified the potential co-expressed diacylglycerol acyltranferase (DGAT) gene associated with triacylglycerol (TAG)-polyestolides biosynthesis, the dynamics of metabolic changes during seed development of O. violaceus as well as its associated regulatory network changes are poorly understood.

RESULTS

In this study, by combining metabolome and transcriptome analysis, we identified that 1,003 metabolites and 22,479 genes were active across four stages of seed development, which were further divided into three main clusters based on the patterns of metabolite accumulation and/or gene expression. Among which, cluster2 was mostly related to diOH-FA biosynthesis pathway. We thus further constructed transcription factor (TF)-structural genes regulatory map for the genes associated with the flavonoids, fatty acids and diOH-FA biosynthesis pathway in this cluster. In particular, several TF families such as bHLH, B3, HD-ZIP, MYB were found to potentially regulate the metabolism associated with the diOH-FA pathway. Among which, multiple candidate TFs with promising potential for increasing the diOH-FA content were identified, and we further traced the evolutionary history of these key genes among species of Brassicaceae.

CONCLUSION

Taken together, our study provides new insight into the gene resources and potential relevant regulatory mechanisms of diOH-FA biosynthesis uniquely in seeds of O. violaceus, which will help to promote the downstream breeding efforts of this potential oilseed crop and advance the bio-lubricant industry.

Members of the tomato NRC4 h-NLR family augment each other in promoting basal immunity

Plants possess an efficient, two-tiered immune system to combat pathogens and pests. Several decades of research have characterized different features of these two well-known tiers, PTI and ETI (Pattern/ Effector-triggered Immunity). NLR (Nucleotide-binding domain Leucine-rich Repeat) receptors have been found to link PTI to ETI, and be required for full potentiation of plant immune responses in several systems. Intra-cellular helper-NLRs (h-NLRs) mediate ETI and have been focused on extensively in recent research. Previously, we investigated the roles of the h-NLR SlNRC4a in tomato immunity, finding that a specific mutation in this gene results in gain of function constitutive defense activation and broad disease resistance. Deletion of the entire NRC4 clade, which contains 3 genes, can compromise tomato immunity. Here, we decided to investigate the role of an additional clade member, SlNRC4b, in basal immunity. We generated a gain of function mutant in SlNRC4b using CRISPR-Cas9, as well as a double gain of function mutant in both genes. Similarly to the slnrc4a mutant, a slnrc4b mutant also possessed increased basal immunity and broad spectrum disease resistance. The double mutant displayed additive effects in some cases, with significant increases in resistance to fungal phytopathogens as compared with each of the single mutants. Our work confirms that the NRC4 family h-NLRs are important in the plant immune system, suggesting that this gene family has the potential to be promising in targeted agricultural adaptation in the Solanaceae family, promoting disease resistance and prevention of yield loss to pathogens.

A metabolic perspective of selection for fruit quality related to apple domestication and improvement

BACKGROUND

Apple is an economically important fruit crop. Changes in metabolism accompanying human-guided evolution can be revealed using a multiomics approach. We perform genome-wide metabolic analysis of apple fruits collected from 292 wild and cultivated accessions representing various consumption types.

RESULTS

We find decreased amounts of certain metabolites, including tannins, organic acids, phenolic acids, and flavonoids as the wild accessions transition to cultivated apples, while lysolipids increase in the "Golden Delicious" to "Ralls Janet" pedigree, suggesting better storage. We identify a total of 222,877 significant single-nucleotide polymorphisms that are associated with 2205 apple metabolites. Investigation of a region from 2.84 to 5.01 Mb on chromosome 16 containing co-mapping regions for tannins, organic acids, phenolic acids, and flavonoids indicates the importance of these metabolites for fruit quality and nutrition during breeding. The tannin and acidity-related genes Myb9-like and PH4 are mapped closely to fruit weight locus fw1 from 3.41 to 3.76 Mb on chromosome 15, a region under selection during domestication. Lysophosphatidylethanolamine (LPE) 18:1, which is suppressed by fatty acid desaturase-2 (FAD2), is positively correlated to fruit firmness. We find the fruit weight is negatively correlated with salicylic acid and abscisic acid levels. Further functional assays demonstrate regulation of these hormone levels by NAC-like activated by Apetala3/Pistillata (NAP) and ATP binding cassette G25 (ABCG25), respectively.

CONCLUSIONS

This study provides a metabolic perspective for selection on fruit quality during domestication and improvement, which is a valuable resource for investigating mechanisms controlling apple metabolite content and quality.

Diversity of tomato leaf form provides novel insights into breeding

Tomato (Solanum lycopersicum L.) is cultivated widely globally. The crop exhibits tremendous morphological variations because of its long breeding history. Apart from the commercial tomato varieties, wild species and heirlooms are grown in certain regions of the world. Since the fruit constitutes the edible part, much of the agronomical research is focused on it. However, recent studies have indicated that leaf morphology influences fruit quality. As leaves are specialized photosynthetic organs and the vascular systems transport the photosynthetic products to sink organs, the architectural characteristics of the leaves have a strong influence on the final fruit quality. Therefore, comprehensive research focusing on both the fruit and leaf morphology is required for further tomato breeding. This review summarizes an overview of knowledge of the basic tomato leaf development, morphological diversification, and molecular mechanisms behind them and emphasizes its importance in breeding. Finally, we discuss how these findings and knowledge can be applied to future tomato breeding.

Deleterious Variation in Natural Populations and Implications for Conservation Genetics

Deleterious mutations decrease reproductive fitness and are ubiquitous in genomes. Given that many organisms face ongoing threats of extinction, there is interest in elucidating the impact of deleterious variation on extinction risk and optimizing management strategies accounting for such mutations. Quantifying deleterious variation and understanding the effects of population history on deleterious variation are complex endeavors because we do not know the strength of selection acting on each mutation. Further, the effect of demographic history on deleterious mutations depends on the strength of selection against the mutation and the degree of dominance. Here we clarify how deleterious variation can be quantified and studied in natural populations. We then discuss how different demographic factors, such as small population size, nonequilibrium population size changes, inbreeding, and gene flow, affect deleterious variation. Lastly, we provide guidance on studying deleterious variation in nonmodel populations of conservation concern.

Toward kingdom-wide analyses of gene expression

Gene expression data for Archaeplastida are accumulating exponentially, with more than 300 000 RNA-sequencing (RNA-seq) experiments available for hundreds of species. The gene expression data stem from thousands of experiments that capture gene expression in various organs, tissues, cell types, (a)biotic perturbations, and genotypes. Advances in software tools make it possible to process all these data in a matter of weeks on modern office computers, giving us the possibility to study gene expression in a kingdom-wide manner for the first time. We discuss how the expression data can be accessed and processed and outline analyses that take advantage of cross-species analyses, allowing us to generate powerful and robust hypotheses about gene function and evolution.

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.

Evolutionary analysis and expression profiling of the HSP70 gene family in response to abiotic stresses in tomato (Solanum lycopersicum)

Heat shock protein 70 (HSP70) genes play essential roles in guarding plants against abiotic stresses, including heat, drought, and salt. In this study, the SlHSP70 gene family in tomatoes has been characterized using bioinformatic tools. 25 putative SlHSP70 genes in the tomato genome were found and classified into five subfamilies, with multi-subcellular localizations. Twelve pairs of gene duplications were identified, and segmental events were determined as the main factor for the gene family expansion. Based on public RNA-seq data, gene expression analysis identified the majority of genes expressed in the examined organelles. Further RNA-seq analysis and then quantitative RT-PCR validation showed that many SlHSP70 members are responsible for cellular feedback to heat, drought, and salt treatments, in which, at least five genes might be potential key players in the stress response. Our results provided a thorough overview of the SlHSP70 gene family in the tomato, which may be useful for the evolutionary and functional analysis of SlHSP70 under abiotic stress conditions.

Gene expression and DNA methylation altering lead to the high oil content in wild allotetraploid peanut (A. monticola)

INTRODUCTION

The wild allotetraploid peanut Arachis monticola contains a higher oil content than the cultivated allotetraploid Arachis hypogaea. Besides the fact that increasing oil content is the most important peanut breeding objective, a proper understanding of its molecular mechanism controlling oil accumulation is still lacking.

METHODS

We investigated this aspect by performing comparative transcriptomics from developing seeds between three wild and five cultivated peanut varieties.

RESULTS

The analyses not only showed species-specific grouping transcriptional profiles but also detected two gene clusters with divergent expression patterns between two species enriched in lipid metabolism. Further analysis revealed that expression alteration of lipid metabolic genes with co-expressed transcription factors in wild peanut led to enhanced activity of oil biogenesis and retarded the rate of lipid degradation. In addition, bisulfite sequencing was conducted to characterize the variation of DNA methylation between wild allotetraploid (245, WH 10025) and cultivated allotetraploid (Z16, Zhh 7720) genotypes. CG and CHG context methylation was found to antagonistically correlate with gene expression during seed development. Differentially methylated region analysis and transgenic assay further illustrated that variations of DNA methylation between wild and cultivated peanuts could affect the oil content via altering the expression of peroxisomal acyl transporter protein (Araip.H6S1B).

DISCUSSION

From the results, we deduced that DNA methylation may negatively regulate lipid metabolic genes and transcription factors to subtly affect oil accumulation divergence between wild and cultivated peanuts. Our work provided the first glimpse on the regulatory mechanism of gene expression altering for oil accumulation in wild peanut and gene resources for future breeding applications.

Environmental response in gene expression and DNA methylation reveals factors influencing the adaptive potential of Arabidopsis lyrata

Understanding what factors influence plastic and genetic variation is valuable for predicting how organisms respond to changes in the selective environment. Here, using gene expression and DNA methylation as molecular phenotypes, we study environmentally induced variation among Arabidopsis lyrata plants grown at lowland and alpine field sites. Our results show that gene expression is highly plastic, as many more genes are differentially expressed between the field sites than between populations. These environmentally responsive genes evolve under strong selective constraint – the strength of purifying selection on the coding sequence is high, while the rate of adaptive evolution is low. We find, however, that positive selection on cis-regulatory variants has likely contributed to the maintenance of genetically variable environmental responses, but such variants segregate only between distantly related populations. In contrast to gene expression, DNA methylation at genic regions is largely insensitive to the environment, and plastic methylation changes are not associated with differential gene expression. Besides genes, we detect environmental effects at transposable elements (TEs): TEs at the high-altitude field site have higher expression and methylation levels, suggestive of a broad-scale TE activation. Compared to the lowland population, plants native to the alpine environment harbor an excess of recent TE insertions, and we observe that specific TE families are enriched within environmentally responsive genes. Our findings provide insight into selective forces shaping plastic and genetic variation. We also highlight how plastic responses at TEs can rapidly create novel heritable variation in stressful conditions.

Vignette of Vigna domestication: From archives to genomics

The genus Vigna comprises fast-growing, diploid legumes, cultivated in tropical and subtropical parts of the world. It comprises more than 200 species among which Vigna angularis, Vigna radiata, Vigna mungo, Vigna aconitifolia, Vigna umbellata, Vigna unguiculata, and Vigna vexillata are of enormous agronomic importance. Human selection along with natural variability within these species encompasses a vital source for developing new varieties. The present review convokes the early domestication history of Vigna species based on archeological pieces of evidence and domestication-related traits (DRTs) together with genetics of domestication. Traces of early domestication of Vigna have been evidenced to spread across several temperate and tropical regions of Africa, Eastern Asia, and few parts of Europe. Several DRTs of Vigna species, such as pod shattering, pod and seed size, dormancy, seed coat, seed color, maturity, and pod dehiscence, can clearly differentiate wild species from their domesticates. With the advancement in next-generation high-throughput sequencing techniques, exploration of genetic variability using recently released reference genomes along with de novo sequencing of Vigna species have provided a framework to perform genome-wide association and functional studies to figure out different genes related to DRTs. In this review, genes and quantitative trait loci (QTLs) related to DRTs of different Vigna species have also been summarized. Information provided in this review will enhance the in-depth understanding of the selective pressures that causes crop domestication along with nature of evolutionary selection made in unexplored Vigna species. Furthermore, correlated archeological and domestication-related genetic evidence will facilitate Vigna species to be considered as suitable model plants.

cis-Regulatory variation affecting gene expression contributes to the improvement of maize kernel size

cis-Regulatory variations contribute to trait evolution and adaptation during crop domestication and improvement. As the most important harvested organ in maize (Zea mays L.), kernel size has undergone intensive selection for size. However, the associations between maize kernel size and cis-regulatory variations remain unclear. We chose two independent association populations to dissect the genetic architecture of maize kernel size together with transcriptomic and genotypic data. The resulting phenotypes reflected a strong influence of population structure on kernel size. Compared with genome-wide association studies (GWASs), which accounted for population structure and relatedness, GWAS based on a naïve or simple linear model revealed additional associated single-nucleotide polymorphisms significantly involved in the conserved pathways controlling seed size in plants. Regulation analyses through expression quantitative trait locus mapping revealed that cis-regulatory variations likely control kernel size by fine-tuning the expression of proximal genes, among which ZmKL1 (GRMZM2G098305) was transgenically validated. We also proved that the pyramiding of the favorable cis-regulatory variations has contributed to the improvement of maize kernel size. Collectively, our results demonstrate that cis-regulatory variations, together with their regulatory genes, provide excellent targets for future maize improvement.

CLASS-II KNOX genes coordinate spatial and temporal ripening in tomato

Fruits can be divided into dry and fleshy types. Dry fruits mature through senescence and fleshy fruits through ripening. Previous studies have indicated that partially common molecular networks could govern fruit maturation in these different fruit types. However, the nature of such networks remains obscure. CLASS-II KNOX genes were shown to regulate the senescence of the Arabidopsis (Arabidopsis thaliana) dry fruits, the siliques, but their roles in fleshy-fruit development are unknown. Here, we investigated the roles of the tomato (Solanum lycopersicum) CLASS-II KNOX (TKN-II) genes in fleshy fruit ripening using knockout alleles of individual genes and an artificial microRNA line (35S:amiR-TKN-II) simultaneously targeting all genes. 35S:amiR-TKN-II plants, as well as a subset of tkn-II single and double mutants, have smaller fruits. Strikingly, the 35S:amiR-TKN-II and tknII3 tknII7/+ fruits showed early ripening of the locular domain while their pericarp ripening was stalled. Further examination of the ripening marker-gene RIPENING INHIBITOR (RIN) expression and 35S:amiR-TKN-II rin-1 mutant fruits suggested that TKN-II genes arrest RIN activity at the locular domain and promote it in the pericarp. These findings imply that CLASS-II KNOX genes redundantly coordinate maturation in both dry and fleshy fruits. In tomato, these genes also control spatial patterns of fruit ripening, utilizing differential regulation of RIN activity at different fruit domains.

Comparative transcriptome analysis on the mangrove Acanthus ilicifolius and its two terrestrial relatives provides insights into adaptation to intertidal habitats

Acanthus is a unique genus covering both mangroves and terrestrial species, and thus is an ideal system to comparatively analyze the mechanisms of mangrove adaptation to intertidal habitats. We performed RNA sequencing of the mangrove plant Acanthus ilicifolius and its two terrestrial relatives, Acanthus leucostachyus and Acanthus mollis. A total of 91,125, 118,290, and 141,640 unigenes were obtained. Simple sequence repeats (SSR) analysis showed that A. ilicifolius had more SSRs, the highest frequency of distribution, and higher in polymorphism potential compared to the two terrestrial relatives. Phylogenetic analyses suggested a relatively recent split between A. ilicifolius and A. leucostachyus, i.e., about 16.76 million years ago (Mya), after their ancestor divergence with A. mollis (32.11 Mya), indicating that speciation of three Acanthus species occurred in the Early to Middle Miocene. Gene Ontology (GO) enrichment revealed that the unique unigenes in A. ilicifolius are predominantly related to rhythmic process, reproductive process and response to stimuli. The accelerated evolution and positive selection analyses indicated that the genus Acanthus migrated from terrestrial to intertidal habitats, where 311 pairs may be under positive selection. Functional enrichment analysis revealed that these genes associated with essential metabolism and biosynthetic pathways such as oxidative phosphorylation, plant hormone signal transduction, photosynthetic carbon fixation and arginine and proline metabolism, are related to the adaptation of A. ilicifolius to intertidal habitats, which are characterized by high salinity and hypoxia. Our results indicate the evolutionary processes and the mechanisms underlying the adaptability of Acanthus to various harsh environments from the arid terrestrial to intertidal habitats.

Internode elongation in energy cane shows remarkable clues on lignocellulosic biomass biosynthesis in Saccharum hybrids

Energy cane is a dedicated crop to high biomass production and selected during Saccharum breeding programs to fit specific industrial needs for 2G bioethanol production. Internode elongation is one of the most important characteristics in Saccharum hybrids due to its relationship with crop yield. In this study, we selected the third internode elongation of the energy cane. To characterize this process, we divided the internode into five sections and performed a detailed transcriptome analysis (RNA-Seq) and cell wall characterization. The histological analyses revealed a remarkable gradient that spans from cell division and protoxylem lignification to the internode maturation and complete vascular bundle lignification. RNA-Seq analysis revealed more than 11,000 differentially expressed genes between the sections internal. Gene ontology analyzes showed enriched categories in each section, as well as the most expressed genes in each section, presented different biological processes. We found that the internode elongation and division zones have a large number of unique genes. Evaluated the specific profile of genes related to primary and secondary cell wall formation, cellulose synthesis, hemicellulose, lignin, and growth-related genes. For each section these genes presented different profiles along the internode in elongation in energy cane. The results of this study provide an overview of the regulation of gene expression of an internode elongation in energy cane. Gene expression analysis revealed promising candidates for transcriptional regulation of energy cane lignification and evidence key genes for the regulation of internode development, which can serve as a basis for understanding the molecular regulatory mechanisms that support the growth and development of plants in the Saccahrum complex.

Comparative Transcriptomics and Metabolites Analysis of Two Closely Related Euphorbia Species Reveal Environmental Adaptation Mechanism and Active Ingredients Difference

Roots of Euphorbia fischeriana and Euphorbia ebracteolata are recorded as the source plant of traditional Chinese medicine "Langdu," containing active ingredients with anticancer and anti-AIDS activity. However, the two species have specific patterns in the graphic distribution. Compared with E. ehracteolata, E. fischeriana distributes in higher latitude and lower temperature areas and might have experienced cold stress adaptation. To reveal the molecular mechanism of environmental adaptation, RNA-seq was performed toward the roots, stems, and leaves of E. fischeriana and E. ehracteolata. A total of 6,830 pairs of putative orthologs between the two species were identified. Estimations of non-synonymous or synonymous substitution rate ratios for these orthologs indicated that 533 of the pairs may be under positive selection (Ka/Ks > 0.5). Functional enrichment analysis revealed that significant proportions of the orthologs were in the TCA cycle, fructose and mannose metabolism, starch and sucrose metabolism, fatty acid biosynthesis, and terpenoid biosynthesis providing insights into how the two closely related Euphorbia species adapted differentially to extreme environments. Consistent with the transcriptome, a higher content of soluble sugars and proline was obtained in E. fischeriana, reflecting the adaptation of plants to different environments. Additionally, 5 primary or secondary metabolites were screened as the biomarkers to distinguish the two species. Determination of 4 diterpenoids was established and performed, showing jolkinolide B as a representative component in E. fischeriana, whereas ingenol endemic to E. ebracteolate. To better study population genetics, EST-SSR markers were generated and tested in 9 species of Euphorbia. A total of 33 of the 68 pairs were screened out for producing clear fragments in at least four species, which will furthermore facilitate the studies on the genetic improvement and phylogenetics of this rapidly adapting taxon. In this study, transcriptome and metabolome analyses revealed the evolution of genes related to cold stress tolerance, biosynthesis of TCA cycle, soluble sugars, fatty acids, and amino acids, consistent with the molecular strategy that genotypes adapting to environment. The key active ingredients of the two species were quantitatively analyzed to reveal the difference in pharmacodynamic substance basis and molecular mechanism, providing insights into rational crude drug use.

Variation of Photosynthetic Induction in Major Horticultural Crops Is Mostly Driven by Differences in Stomatal Traits

Under natural conditions, irradiance frequently fluctuates, causing net photosynthesis rate (A) to respond slowly and reducing the yields. We quantified the genotypic variation of photosynthetic induction in 19 genotypes among the following six horticultural crops: basil, chrysanthemum, cucumber, lettuce, tomato, and rose. Kinetics of photosynthetic induction and the stomatal opening were measured by exposing shade-adapted leaves (50 μmol m-2 s-1) to a high irradiance (1000 μmol m-2 s-1) until A reached a steady state. Rubisco activation rate was estimated by the kinetics of carboxylation capacity, which was quantified using dynamic A vs. [CO2] curves. Generally, variations in photosynthetic induction kinetics were larger between crops and smaller between cultivars of the same crop. Time until reaching 20-90% of full A induction varied by 40-60% across genotypes, and this was driven by a variation in the stomatal opening rather than Rubisco activation kinetics. Stomatal conductance kinetics were partly determined by differences in the stomatal size and density; species with densely packed, smaller stomata (e.g., cucumber) tended to open their stomata faster, adapting stomatal conductance more rapidly and efficiently than species with larger but fewer stomata (e.g., chrysanthemum). We conclude that manipulating stomatal traits may speed up photosynthetic induction and growth of horticultural crops under natural irradiance fluctuations.

Evolutionary Ecology of Plant-Arthropod Interactions in Light of the "Omics" Sciences: A Broad Guide

Aboveground plant-arthropod interactions are typically complex, involving herbivores, predators, pollinators, and various other guilds that can strongly affect plant fitness, directly or indirectly, and individually, synergistically, or antagonistically. However, little is known about how ongoing natural selection by these interacting guilds shapes the evolution of plants, i.e., how they affect the differential survival and reproduction of genotypes due to differences in phenotypes in an environment. Recent technological advances, including next-generation sequencing, metabolomics, and gene-editing technologies along with traditional experimental approaches (e.g., quantitative genetics experiments), have enabled far more comprehensive exploration of the genes and traits involved in complex ecological interactions. Connecting different levels of biological organization (genes to communities) will enhance the understanding of evolutionary interactions in complex communities, but this requires a multidisciplinary approach. Here, we review traditional and modern methods and concepts, then highlight future avenues for studying the evolution of plant-arthropod interactions (e.g., plant-herbivore-pollinator interactions). Besides promoting a fundamental understanding of plant-associated arthropod communities' genetic background and evolution, such knowledge can also help address many current global environmental challenges.

Transcriptome Analysis of Populus euphratica under Salt Treatment and PeERF1 Gene Enhances Salt Tolerance in Transgenic Populus alba × Populus glandulosa

Populus euphratica is mainly distributed in desert environments with dry and hot climate in summer and cold in winter. Compared with other poplars, P. euphratica is more resistant to salt stress. It is critical to investigate the transcriptome and molecular basis of salt tolerance in order to uncover stress-related genes. In this study, salt-tolerant treatment of P. euphratica resulted in an increase in osmo-regulatory substances and recovery of antioxidant enzymes. To improve the mining efficiency of candidate genes, the analysis combining both the transcriptome WGCNA and the former GWAS results was selected, and a range of key regulatory factors with salt resistance were found. The PeERF1 gene was highly connected in the turquoise modules with significant differences in salt stress traits, and the expression levels were significantly different in each treatment. For further functional verification of PeERF1, we obtained stable overexpression and dominant suppression transgenic lines by transforming into Populus alba × Populusglandulosa. The growth and physiological characteristics of the PeERF1 overexpressed plants were better than that of the wild type under salt stress. Transcriptome analysis of leaves of transgenic lines and WT revealed that highly enriched GO terms in DEGs were associated with stress responses, including abiotic stimuli responses, chemical responses, and oxidative stress responses. The result is helpful for in-depth analysis of the salt tolerance mechanism of poplar. This work provides important genes for poplar breeding with salt tolerance.

Comparative transcriptomics of wild and commercial Citrus during early ripening reveals how domestication shaped fruit gene expression

BACKGROUND

Interspecific hybridizations and admixtures were key in Citrus domestication, but very little is known about their impact at the transcriptomic level. To determine the effects of genome introgressions on gene expression, the transcriptomes of the pulp and flavedo of three pure species (citron, pure mandarin and pummelo) and four derived domesticated genetic admixtures (sour orange, sweet orange, lemon and domesticated mandarin) have been analyzed at color break.

RESULTS

Many genes involved in relevant physiological processes for domestication, such sugar/acid metabolism and carotenoid/flavonoid synthesis, were differentially expressed among samples. In the low-sugar, highly acidic species lemon and citron, many genes involved in sugar metabolism, the TCA cycle and GABA shunt displayed a reduced expression, while the P-type ATPase CitPH5 and most subunits of the vacuolar ATPase were overexpressed. The red-colored species and admixtures were generally characterized by the overexpression in the flavedo of specific pivotal genes involved in the carotenoid biosynthesis, including phytoene synthase, ζ-carotene desaturase, β-lycopene cyclase and CCD4b, a carotenoid cleavage dioxygenase. The expression patterns of many genes involved in flavonoid modifications, especially the flavonoid and phenylpropanoid O-methyltransferases showed extreme diversity. However, the most noticeable differential expression was shown by a chalcone synthase gene, which catalyzes a key step in the biosynthesis of flavonoids. This chalcone synthase was exclusively expressed in mandarins and their admixed species, which only expressed the mandarin allele. In addition, comparisons between wild and domesticated mandarins revealed that the major differences between their transcriptomes concentrate in the admixed regions.

CONCLUSION

In this work we present a first study providing broad evidence that the genome introgressions that took place during citrus domestication largely shaped gene expression in their fruits.

Identification and Characterization of Key Genes Responsible for Weedy and Cultivar Growth Types in Soybean

In cultivated plants, shoot morphology is an important factor that influences crop economic value. However, the effects of gene expression patterns on shoot morphology are not clearly understood. In this study, the molecular mechanism behind shoot morphology (including leaf, stem, and node) was analyzed using RNA sequencing to compare weedy (creeper) and cultivar (stand) growth types obtained in F7 derived from a cross of wild and cultivated soybeans. A total of 12,513 (in leaves), 14,255 (in stems), and 11,850 (in nodes) differentially expressed genes were identified among weedy and cultivar soybeans. Comparative transcriptome and expression analyses revealed 22 phytohormone-responsive genes. We found that GIBBERELLIN 2-OXIDASE 8 (GA2ox), SPINDLY (SPY), FERONIA (FER), AUXIN RESPONSE FACTOR 8 (ARF8), CYTOKININ DEHYDROGENASE-1 (CKX1), and ARABIDOPSIS HISTIDINE KINASE-3 (AHK3), which are crucial phytohormone response genes, were mainly regulated in the shoot of weedy and cultivar types. These results indicate that interactions between phytohormone signaling genes regulate shoot morphology in weedy and cultivar growth type plants. Our study provides insights that are useful for breeding and improving crops to generate high-yield soybean varieties.

Early Domestication History of Asian Rice Revealed by Mutations and Genome-Wide Analysis of Gene Genealogies

BACKGROUND

Asian rice (Oryza sativa L.) has been a model plant but its cultivation history is inadequately understood, and its origin still under debate. Several enigmas remain, including how this annual crop shifted its growth habit from its perennial ancestor, O. rufipogon, why genetic divergence between indica and japonica appears older than the history of human domestication, and why some domestication genes do not show signals of introgression between subgroups. Addressing these issues may benefit both basic research and rice breeding.

RESULTS

Gene genealogy-based mutation (GGM) analysis shows that history of Asian rice is divided into two phases (Phase I and II) of about equal lengths. Mutations occurred earlier than the partition of indica and japonica to Os genome mark Phase-I period. We diagnosed 91 such mutations among 101 genes sampled across 12 chromosomes of Asian rice and its wild relatives. Positive selection, detected more at 5' regions than at coding regions of some of the genes, involved 22 loci (e.g., An-1, SH4, Rc, Hd3a, GL3.2, OsMYB3, OsDFR, and OsMYB15), which affected traits from easy harvesting, grain color, flowering time, productivity, to likely taste and tolerance. Phase-I mutations of OsMYB3, OsHd3a and OsDFR were experimentally tested and all caused enhanced functions of the genes in vivo. Phase-II period features separate cultivations, lineage-specific selection, and expanded domestication to more genes. Further genomic analysis, along with phenotypic comparisons, indicates that O. sativa is hybrid progeny of O. rufipogon and O. nivara, inherited slightly more genes of O. rufipogon. Congruently, modern alleles of the sampled genes are approximately 6% ancient, 38% uni-specific, 40% bi-specific (mixed), and 15% new after accumulating significant mutations. Results of sequencing surveys across modern cultivars/landraces indicate locus-specific usages of various alleles while confirming the associated mutations.

CONCLUSIONS

Asian rice was initially domesticated as one crop and later separate selection mediated by human resulted in its major subgroups. This history and the hybrid origin well explain previous puzzles. Positive selection, particularly in 5' regions, was the major force underlying trait domestication. Locus-specific domestication can be characterized and the result may facilitate breeders in developing better rice varieties in future.

All-flesh fruit in tomato is controlled by reduced expression dosage of AFF through a structural variant mutation in the promoter

The formation of locule gel is an important process in tomato and is a typical characteristic of berry fruit. In this study, we examined a natural tomato mutant that produces all-flesh fruit (AFF) in which the locule tissue remains in a solid state during fruit development. We constructed different genetic populations to fine-map the causal gene for this trait and identified SlMBP3 as the locus conferring the locule gel formation, which we rename as AFF. We determined the causal mutation as a 416-bp deletion in the promoter region of AFF, which reduces its expression dosage. Generally, this sequence is highly conserved among Solanaceae, as well as within the tomato germplasm. Using BC6 near-isogenic lines, we determined that the reduced expression dosage of AFF did not affect the normal development of seeds, whilst producing unique, non-liquefied locule tissue that was distinct from that of normal tomatoes in terms of metabolic components. Combined analysis using mRNA-seq and metabolomics indicated the importance of AFF in locule tissue liquefaction. Our findings provide insights into fruit-type differentiation in Solanaceae crops and also present the basis for future applications of AFF in tomato breeding programs.

The landscape of transcriptional and translational changes over 22 years of bacterial adaptation

Organisms can adapt to an environment by taking multiple mutational paths. This redundancy at the genetic level, where many mutations have similar phenotypic and fitness effects, can make untangling the molecular mechanisms of complex adaptations difficult. Here, we use the Escherichia coli long-term evolution experiment (LTEE) as a model to address this challenge. To understand how different genomic changes could lead to parallel fitness gains, we characterize the landscape of transcriptional and translational changes across 12 replicate populations evolving in parallel for 50,000 generations. By quantifying absolute changes in mRNA abundances, we show that not only do all evolved lines have more mRNAs but that this increase in mRNA abundance scales with cell size. We also find that despite few shared mutations at the genetic level, clones from replicate populations in the LTEE are remarkably similar in their gene expression patterns at both the transcriptional and translational levels. Furthermore, we show that the majority of the expression changes are due to changes at the transcriptional level with very few translational changes. Finally, we show how mutations in transcriptional regulators lead to consistent and parallel changes in the expression levels of downstream genes. These results deepen our understanding of the molecular mechanisms underlying complex adaptations and provide insights into the repeatability of evolution.

Comparing Early Transcriptomic Responses of 18 Soybean (Glycine max) Genotypes to Iron Stress

Iron deficiency chlorosis (IDC) is an abiotic stress that negatively affects soybean (Glycine max [L.] Merr.) production. Much of our knowledge of IDC stress responses is derived from model plant species. Gene expression, quantitative trait loci (QTL) mapping, and genome-wide association studies (GWAS) performed in soybean suggest that stress response differences exist between model and crop species. Our current understanding of the molecular response to IDC in soybeans is largely derived from gene expression studies using near-isogenic lines differing in iron efficiency. To improve iron efficiency in soybeans and other crops, we need to expand gene expression studies to include the diversity present in germplasm collections. Therefore, we collected 216 purified RNA samples (18 genotypes, two tissue types [leaves and roots], two iron treatments [sufficient and deficient], three replicates) and used RNA sequencing to examine the expression differences of 18 diverse soybean genotypes in response to iron deficiency. We found a rapid response to iron deficiency across genotypes, most responding within 60 min of stress. There was little evidence of an overlap of specific differentially expressed genes, and comparisons of gene ontology terms and transcription factor families suggest the utilization of different pathways in the stress response. These initial findings suggest an untapped genetic potential within the soybean germplasm collection that could be used for the continued improvement of iron efficiency in soybean.

Deleterious protein-coding variants in diverse cattle breeds of the world

The domestication of wild animals has resulted in a reduction in effective population sizes, which can affect the deleterious mutation load of domesticated breeds. In addition, artificial selection contributes to the accumulation of deleterious mutations because of an increased rate of inbreeding among domesticated animals. Since founder population sizes and artificial selection differ between cattle breeds, their deleterious mutation load can vary. We investigated this question by using whole-genome data from 432 animals belonging to 54 worldwide cattle breeds. Our analysis revealed a negative correlation between genomic heterozygosity and nonsynonymous-to-silent diversity ratio, which suggests a higher proportion of single nucleotide variants (SNVs) affecting proteins in low-diversity breeds. Our results also showed that low-diversity breeds had a larger number of high-frequency (derived allele frequency (DAF) > 0.51) deleterious SNVs than high-diversity breeds. An opposite trend was observed for the low-frequency (DAF ≤ 0.51) deleterious SNVs. Overall, the number of high-frequency deleterious SNVs was larger in the genomes of taurine cattle breeds than of indicine breeds, whereas the number of low-frequency deleterious SNVs was larger in the genomes of indicine cattle than in those of taurine cattle. Furthermore, we observed significant variation in the counts of deleterious SNVs within taurine breeds. The variations in deleterious mutation load between taurine and indicine breeds could be attributed to the population sizes of the wild progenitors before domestication, whereas the variations observed within taurine breeds could be due to differences in inbreeding level, strength of artificial selection, and/or founding population size. Our findings imply that the incidence of genetic diseases can vary between cattle breeds.

Transcriptome and Resequencing Analyses Provide Insight into Differences in Organic Acid Accumulation in Two Pear Varieties

Fruit acidity is one of the main determinants of fruit flavor and a target trait in fruit breeding. However, the genomic mechanisms governing acidity variation among different pear varieties remain poorly understood. In this study, two pear varieties with contrasting organic acid levels, ‘Dangshansuli’ (low-acidity) and ‘Amute’ (high-acidity), were selected, and a combination of transcriptome and population genomics analyses were applied to characterize their patterns of gene expression and genetic variation. Based on RNA-seq data analysis, differentially expressed genes (DEGs) involved in organic acid metabolism and accumulation were identified. Weighted correlation network analysis (WGCNA) revealed that nine candidate TCA (tricarboxylic acid)-related DEGs and three acid transporter-related DEGs were located in three key modules. The regulatory networks of the above candidate genes were also predicted. By integrating pear resequencing data, two domestication-related genes were found to be upregulated in ‘Amute’, and this trend was further validated for other pear varieties with high levels of organic acid, suggesting distinct selective sweeps during pear dissemination and domestication. Collectively, this study provides insight into organic acid differences related to expression divergence and domestication in two pear varieties, pinpointing several candidate genes for the genetic manipulation of acidity in pears.

Gene expression evolution in pattern-triggered immunity within Arabidopsis thaliana and across Brassicaceae species

Plants recognize surrounding microbes by sensing microbe-associated molecular patterns (MAMPs) to activate pattern-triggered immunity (PTI). Despite their significance for microbial control, the evolution of PTI responses remains largely uncharacterized. Here, by employing comparative transcriptomics of six Arabidopsis thaliana accessions and three additional Brassicaceae species to investigate PTI responses, we identified a set of genes that commonly respond to the MAMP flg22 and genes that exhibit species-specific expression signatures. Variation in flg22-triggered transcriptome responses across Brassicaceae species was incongruent with their phylogeny, while expression changes were strongly conserved within A. thaliana. We found the enrichment of WRKY transcription factor binding sites in the 5'-regulatory regions of conserved and species-specific responsive genes, linking the emergence of WRKY-binding sites with the evolution of gene expression patterns during PTI. Our findings advance our understanding of the evolution of the transcriptome during biotic stress.

Using precision phenotyping to inform de novo domestication

An update on the use of precision phenotyping to assess the potential of lesser cultivated species as candidates for de novo domestication or similar development for future agriculture.

Large-scale genomic analysis reveals the genetic cost of chicken domestication

BACKGROUND

Species domestication is generally characterized by the exploitation of high-impact mutations through processes that involve complex shifting demographics of domesticated species. These include not only inbreeding and artificial selection that may lead to the emergence of evolutionary bottlenecks, but also post-divergence gene flow and introgression. Although domestication potentially affects the occurrence of both desired and undesired mutations, the way wild relatives of domesticated species evolve and how expensive the genetic cost underlying domestication is remain poorly understood. Here, we investigated the demographic history and genetic load of chicken domestication.

RESULTS

We analyzed a dataset comprising over 800 whole genomes from both indigenous chickens and wild jungle fowls. We show that despite having a higher genetic diversity than their wild counterparts (average π, 0.00326 vs. 0.00316), the red jungle fowls, the present-day domestic chickens experienced a dramatic population size decline during their early domestication. Our analyses suggest that the concomitant bottleneck induced 2.95% more deleterious mutations across chicken genomes compared with red jungle fowls, supporting the "cost of domestication" hypothesis. Particularly, we find that 62.4% of deleterious SNPs in domestic chickens are maintained in heterozygous states and masked as recessive alleles, challenging the power of modern breeding programs to effectively eliminate these genetic loads. Finally, we suggest that positive selection decreases the incidence but increases the frequency of deleterious SNPs in domestic chicken genomes.

CONCLUSION

This study reveals a new landscape of demographic history and genomic changes associated with chicken domestication and provides insight into the evolutionary genomic profiles of domesticated animals managed under modern human selection.

Identifying genomic regions targeted during eggplant domestication using transcriptome data

Identifying genes and traits that have diverged during domestication provides key information of importance for maintaining and even increasing yield and nutrients in existing crops. A “bottom-up” population genetics approach was used to identify signatures of selection across the eggplant genome, to better understand the process of domestication. RNA-seq data were obtained for 4 wild eggplants (Solanum insanum L.) and 16 domesticated eggplants (S. melongena L.) and mapped to the eggplant genome. Single-nucleotide polymorphism (SNPs) exhibiting signatures of selection in domesticates were identified as those exhibiting high F_ST between the 2 populations (evidence of significant divergence) and low π for the domesticated population (indicative of a selective sweep). Some of these regions appear to overlap with previously identified quantitative trait loci for domestication traits. Genes in regions of linkage disequilibrium surrounding these SNPs were searched against the Arabidopsis thaliana and tomato genomes to find orthologs. Subsequent gene ontology (GO) enrichment analysis identified over-representation of GO terms related to photosynthesis and response to the environment. This work reveals genomic changes involved in eggplant domestication and improvement, and how this compares to observed changes in the tomato genome, revealing shared chromosomal regions involved in the domestication of both species.

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.

The Relevance of a Physiological-Stage Approach Study of the Molecular and Environmental Factors Regulating Seed Germination in Wild Plants

Germination represents the culmination of the seed developmental program and is affected by the conditions prevailing during seed maturation in the mother plant. During maturation, the dormancy condition and tolerance to dehydration are established. These characteristics are modulated by the environment to which they are subjected, having an important impact on wild species. In this work, a review was made of the molecular bases of the maturation, the processes of dormancy imposition and loss, as well as the germination process in different wild species with different life histories, and from diverse habitats. It is also specified which of these species present a certain type of management. The impact that the domestication process has had on certain characteristics of the seed is discussed, as well as the importance of determining physiological stages based on morphological characteristics, to face the complexities of the study of these species and preserve their genetic diversity and physiological responses.

Transcriptomic Analysis of Rice Plants Overexpressing PsGAPDH in Response to Salinity Stress

In plants, glyceraldehyde-3-phosphate dehydrogenase (GAPDH) is a main enzyme in the glycolytic pathway. It plays an essential role in glycerolipid metabolism and response to various stresses. To examine the function of PsGAPDH (Pleurotus sajor-caju GAPDH) in response to abiotic stress, we generated transgenic rice plants with single-copy/intergenic/homozygous overexpression PsGAPDH (PsGAPDH-OX) and investigated their responses to salinity stress. Seedling growth and germination rates of PsGAPDH-OX were significantly increased under salt stress conditions compared to those of the wild type. To elucidate the role of PsGAPDH-OX in salt stress tolerance of rice, an Illumina HiSeq 2000 platform was used to analyze transcriptome profiles of leaves under salt stress. Analysis results of sequencing data showed that 1124 transcripts were differentially expressed. Using the list of differentially expressed genes (DEGs), functional enrichment analyses of DEGs such as Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were performed. KEGG pathway enrichment analysis revealed that unigenes exhibiting differential expression were involved in starch and sucrose metabolism. Interestingly, trehalose-6-phosphate synthase (TPS) genes, of which expression was enhanced by abiotic stress, showed a significant difference in PsGAPDH-OX. Findings of this study suggest that PsGAPDH plays a role in the adaptation of rice plants to salt stress.

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

The process of domestication requires the rapid transformation of the wild morphology into the cultivated forms that humans select for. This process often takes place through changes in the regulation of genes, yet, there is no definite pattern on the role of cis- and trans-acting regulatory variations in the domestication of the fruit among crops. Using allele-specific expression and network analyses, we characterized the regulatory patterns and the inheritance of gene expression in wild and cultivated accessions of chili pepper, a crop with remarkable fruit morphological variation. We propose that gene expression differences associated to the cultivated form are best explained by cis-regulatory hubs acting through trans-regulatory cascades. We show that in cultivated chili, the expression of genes associated with fruit morphology is partially recessive with respect to those in the wild relative, consistent with the hybrid fruit phenotype. Decreased expression of fruit maturation and growth genes in cultivated chili suggest that selection for loss-of-function took place in its domestication. Trans-regulatory changes underlie the majority of the genes showing regulatory divergence and had larger effect sizes on gene expression than cis-regulatory variants. Network analysis of selected cis-regulated genes, including ARP9 and MED25, indicated their interaction with many transcription factors involved in organ growth and fruit ripening. Differentially expressed genes linked to cis-regulatory variants and their interactions with downstream trans-acting genes have the potential to drive the morphological differences observed between wild and cultivated fruits and provide an attractive mechanism of morphological transformation during the domestication of the chili pepper.

Comparative analysis of the transcriptomes of two rice subspecies during domestication

Two subspecies of rice, Oryza sativa ssp. indica and O. sativa ssp. japonica, with reproductive isolation and differences in morphology and phenotypic differences, were established during the process of rice domestication. To understand how domestication has changed the transcriptomes of the two rice subspecies and given rise to the phenotypic differences, we obtained approximately 700 Gb RNA-Seq data from 26 indica and 25 japonica accessions, and identified 97,005 transcribed fragments and 4579 novel transcriptionally active regions. The two rice subspecies had significantly different gene expression profiles, we identified 1,357 (3.3% in all genes) differentially expressed genes (DEGs) between indica and japonica rice. Combining existing gene function studies, it is found that some of these differential genes are related to the differentiation of the two subspecies, such as grain shape and cold tolerance, etc. Functional annotation of these DEGs indicates that they are involved in cell wall biosynthesis and reproductive processes. Furthermore, compared with the non-DEGs, the DEGs from both subspecies had more 5'flanking regions with low polymorphism to divergence ratios, indicating a stronger positive selection pressure on the regulation of the DEGs. This study improves our understanding of the rice genome by comparatively analyzing the transcriptomes of indica and japonica rice and identifies DEGs those may be responsible for the reproductive isolation and phenotypic differences between the two rice subspecies.

Integrating multi-omics data for crop improvement

Our agricultural systems are now in urgent need to secure food for a growing world population. To meet this challenge, we need a better characterization of plant genetic and phenotypic diversity. The combination of genomics, transcriptomics and metabolomics enables a deeper understanding of the mechanisms underlying the complex architecture of many phenotypic traits of agricultural relevance. We review the recent advances in plant genomics to see how these can be integrated with broad molecular profiling approaches to improve our understanding of plant phenotypic variation and inform crop breeding strategies.

Increased water use efficiency in miR396-downregulated tomato plants

MicroRNAs regulate plant development and responses to biotic and abiotic stresses but their impact on water use efficiency (WUE) is poorly known. Increasing WUE is a major task in crop improvement programs aimed to meet the challenges posed by the reduction in water availability associated with the ongoing climatic change. We have examined the physiological and molecular response to water stress of tomato (Solanum lycopersicum L.) plants downregulated for miR396 by target mimicry. In water stress conditions, miR396-downregulated plants displayed reduced transpiration and a less then proportional decrease in the photosynthetic rate that resulted in higher WUE. The increase in WUE was associated with faster foliar accumulation of abscisic acid (ABA), with the induction of several drought-protective genes and with the activation of the jasmonic acid (JA) and γ-aminobutyric acid (GABA) pathways. We propose a model in which the downregulation of miR396 leads to the activation of a complex molecular response to water stress. This response acts synergistically with a set of leaf morphological modifications to increase stomatal closure and preserve the efficiency of the photosynthetic activity, ultimately resulting in higher WUE.

Reproductive Proteins Evolve Faster Than Non-reproductive Proteins Among Solanum Species

Elevated rates of evolution in reproductive proteins are commonly observed in animal species, and are thought to be driven by the action of sexual selection and sexual conflict acting specifically on reproductive traits. Whether similar patterns are broadly observed in other biological groups is equivocal. Here, we examine patterns of protein divergence among wild tomato species (Solanum section Lycopersicon), to understand forces shaping the evolution of reproductive genes in this diverse, rapidly evolving plant clade. By comparing rates of molecular evolution among loci expressed in reproductive and non-reproductive tissues, our aims were to test if: (a) reproductive-specific loci evolve more rapidly, on average, than non-reproductive loci; (b) 'male'-specific loci evolve at different rates than 'female'-specific loci; (c) genes expressed exclusively in gametophytic (haploid) tissue evolve differently from genes expressed in sporophytic (diploid) tissue or in both tissue types; and (d) mating system variation (a potential proxy for the expected strength of sexual selection and/or sexual conflict) affects patterns of protein evolution. We observed elevated evolutionary rates in reproductive proteins. However, this pattern was most evident for female- rather than male-specific loci, both broadly and for individual loci inferred to be positively selected. These elevated rates might be facilitated by greater tissue-specificity of reproductive proteins, as faster rates were also associated with more narrow expression domains. In contrast, we found little evidence that evolutionary rates are consistently different in loci experiencing haploid selection (gametophytic-exclusive loci), or in lineages with quantitatively different mating systems. Overall while reproductive protein evolution is generally elevated in this diverse plant group, some specific patterns of evolution are more complex than those reported in other (largely animal) systems, and include a more prominent role for female-specific loci among adaptively evolving genes.

The Road to Sorghum Domestication: Evidence From Nucleotide Diversity and Gene Expression Patterns

Native African cereals (sorghum, millets) ensure food security to millions of low-income people from low fertility and drought-prone regions of Africa and Asia. In spite of their agronomic importance, the genetic bases of their phenotype and adaptations are still not well-understood. Here we focus on Sorghum bicolor, which is the fifth cereal worldwide for grain production and constitutes the staple food for around 500 million people. We leverage transcriptomic resources to address the adaptive consequences of the domestication process. Gene expression and nucleotide variability were analyzed in 11 domesticated and nine wild accessions. We documented a downregulation of expression and a reduction of diversity both in nucleotide polymorphism (30%) and gene expression levels (18%) in domesticated sorghum. These findings at the genome-wide level support the occurrence of a global reduction of diversity during the domestication process, although several genes also showed patterns consistent with the action of selection. Nine hundred and forty-nine genes were significantly differentially expressed between wild and domesticated gene pools. Their functional annotation points to metabolic pathways most likely contributing to the sorghum domestication syndrome, such as photosynthesis and auxin metabolism. Coexpression network analyzes revealed 21 clusters of genes sharing similar expression patterns. Four clusters (totaling 2,449 genes) were significantly enriched in differentially expressed genes between the wild and domesticated pools and two were also enriched in domestication and improvement genes previously identified in sorghum. These findings reinforce the evidence that the combined and intricated effects of the domestication and improvement processes do not only affect the behaviors of a few genes but led to a large rewiring of the transcriptome. Overall, these analyzes pave the way toward the identification of key domestication genes valuable for genetic resources characterization and breeding purposes.

Comparative Analysis Delineates the Transcriptional Resistance Mechanisms for Pod Borer Resistance in the Pigeonpea Wild Relative Cajanus scarabaeoides (L.) Thouars

Insect pests pose a serious threat to global food production. Pod borer (Helicoverpa armigera (Hübner)) is one of the most destructive pests of leguminous crops. The use of host resistance has been an effective, environmentally friendly and sustainable approach for controlling several agricultural pests. The exploitation of natural variations in crop wild relatives could yield pest-resistant crop varieties. In this study, we used a high-throughput transcriptome profiling approach to investigate the defense mechanisms of susceptible cultivated and tolerant wild pigeonpea genotypes against H. armigera infestation. The wild genotype displayed elevated pest-induced gene expression, including the enhanced induction of phytohormone and calcium/calmodulin signaling, transcription factors, plant volatiles and secondary metabolite genes compared to the cultivated control. The biosynthetic and regulatory processes associated with flavonoids, terpenes and glucosinolate secondary metabolites showed higher accumulations in the wild genotype, suggesting the existence of distinct tolerance mechanisms. This study provides insights into the molecular mechanisms underlying insect resistance in the wild pigeonpea genotype. This information highlights the indispensable role of crop wild relatives as a source of crucial genetic resources that could be important in devising strategies for crop improvement with enhanced pest resistance.

Plant TDP1 (Tyrosyl-DNA Phosphodiesterase 1): A Phylogenetic Perspective and Gene Expression Data Mining

The TDP1 (tyrosyl-DNA phosphodiesterase 1) enzyme removes the non-specific covalent intermediates between topoisomerase I and DNA, thus playing a crucial role in preventing DNA damage. While mammals possess only one TDP1 gene, in plants two genes (TDP1α and TDP1β) are present constituting a small gene subfamily. These display a different domain structure and appear to perform non-overlapping functions in the maintenance of genome integrity. Namely, the HIRAN domain identified in TDP1β is involved in the interaction with DNA during the recognition of stalled replication forks. The availability of transcriptomic databases in a growing variety of experimental systems provides new opportunities to fill the current gaps of knowledge concerning the evolutionary origin and the specialized roles of TDP1 genes in plants. Whereas a phylogenetic approach has been used to track the evolution of plant TDP1 protein, transcriptomic data from a selection of representative lycophyte, eudicots, and monocots have been implemented to explore the transcriptomic dynamics in different tissues and a variety of biotic and abiotic stress conditions. While the phylogenetic analysis indicates that TDP1α is of non-plant origin and TDP1β is plant-specific originating in ancient vascular plants, the gene expression data mining comparative analysis pinpoints for tissue- and stress-specific responses.

Identification and Characterization of PLATZ Transcription Factors in Wheat

The PLATZ (plant AT-rich protein and zinc-binding protein) transcription factor family is a class of plant-specific zinc-dependent DNA-binding proteins. PLATZ has essential roles in seed endosperm development, as well as promoting cell proliferation duration in the earlier stages of the crops. In the present study, 62 TaPLATZ genes were identified from the wheat genome, and they were unequally distributed on 15 chromosomes. According to the phylogenetic analysis, 62 TaPLATZ genes were classified into six groups, including two groups that were unique in wheat. Members in the same groups shared similar exon-intron structures. The polyploidization, together with genome duplication of wheat, plays a crucial role in the expansion of the TaPLATZs family. Transcriptome data indicated a distinct divergence expression pattern of TaPLATZ genes that could be clustered into four modules. The TaPLATZs in Module b possessed a seed-specific expression pattern and displayed obvious high expression in the earlier development stage of seeds. Subcellular localization data of TaPLATZs suggesting that they likely perform a function as a conventional transcription factor. This study provides insight into understanding the structure divergence, evolutionary features, expression profiles, and potential function of PLATZ in wheat.

Comparative transcriptome analysis reveals evolutionary divergence and shared network of cold and salt stress response in diploid D-genome cotton

BACKGROUND

Wild species of cotton are excellent resistance to abiotic stress. Diploid D-genome cotton showed abundant phenotypic diversity and was the putative donor species of allotetraploid cotton which produce the largest textile natural fiber.

RESULTS

A total of 41,053 genes were expressed in all samples by mapping RNA-seq Illumina reads of G. thurberi (D1), G. klotzschianum (D3-k), G. raimondii (D5) and G. trilobum (D8) to reference genome. The numbers of differently expressed genes (DEGs) were significantly higher under cold stress than salt stress. However, 34.1% DEGs under salt stress were overlapped with cold stress in four species. Notably, a potential shared network (cold and salt response, including 16 genes) was mined out by gene co-expression analysis. A total of 47,180-55,548 unique genes were identified in four diploid species by De novo assembly. Furthermore, 163, 344, 330, and 161 positively selected genes (PSGs) were detected in thurberi, G. klotzschianum, G. raimondii and G. trilobum by evolutionary analysis, respectively, and 9.5-17% PSGs of four species were DEGs in corresponding species under cold or salt stress. What's more, most of PSGs were enriched GO term related to response to stimulation. G. klotzschianum showed the best tolerance under both cold and salt stress. Interestingly, we found that a RALF-like protein coding gene not only is PSGs of G. klotzschianum, but also belongs to the potential shared network.

CONCLUSION

Our study provided new evidence that gene expression variations of evolution by natural selection were essential drivers of the morphological variations related to environmental adaptation during evolution. Additionally, there exist shared regulated networks under cold and salt stress, such as Ca2+ signal transduction and oxidation-reduction mechanisms. Our work establishes a transcriptomic selection mechanism for altering gene expression of the four diploid D-genome cotton and provides available gene resource underlying multi-abiotic resistant cotton breeding strategy.

Genomic inference of complex domestication histories in three Solanaceae species

Domestication is a human-induced selection process that imprints the genomes of domesticated populations over a short evolutionary time scale and that occurs in a given demographic context. Reconstructing historical gene flow, effective population size changes and their timing is therefore of fundamental interest to understand how plant demography and human selection jointly shape genomic divergence during domestication. Yet, the comparison under a single statistical framework of independent domestication histories across different crop species has been little evaluated so far. Thus, it is unclear whether domestication leads to convergent demographic changes that similarly affect crop genomes. To address this question, we used existing and new transcriptome data on three crop species of Solanaceae (eggplant, pepper and tomato), together with their close wild relatives. We fitted twelve demographic models of increasing complexity on the unfolded joint allele frequency spectrum for each wild/crop pair, and we found evidence for both shared and species-specific demographic processes between species. A convergent history of domestication with gene flow was inferred for all three species, along with evidence of strong reduction in the effective population size during the cultivation stage of tomato and pepper. The absence of any reduction in size of the crop in eggplant stands out from the classical view of the domestication process; as does the existence of a "protracted period" of management before cultivation. Our results also suggest divergent management strategies of modern cultivars among species as their current demography substantially differs. Finally, the timing of domestication is species-specific and supported by the few historical records available.

What does it take to be antiviral? An Argonaute-centered perspective on plant antiviral defense

RNA silencing is a major mechanism of constitutive antiviral defense in plants, mediated by a number of proteins, including the Dicer-like (DCL) and Argonaute (AGO) endoribonucleases. Both DCL and AGO protein families comprise multiple members. In particular, the AGO protein family has expanded considerably in different plant lineages, with different family members having specialized functions. Although the general mode of action of AGO proteins is well established, the properties that make different AGO proteins more or less efficient at targeting viruses are less well understood. In this report, we review methodologies used to study AGO antiviral activity and current knowledge about which AGO family members are involved in antiviral defense. In addition, we discuss what is known about the different properties of AGO proteins thought to be associated with this function.