A near complete genome assembly of chia assists in identification of key fatty acid desaturases in developing seeds.

Chia is an annual crop whose seeds have the highest content of α-linolenic acid (ALA) of any plant known to date. We generated a high-quality assembly of the chia genome using circular consensus sequencing (CCS) of PacBio. The assembled six chromosomes are composed of 21 contigs and have a total length of 361.7 Mb. Genome annotation revealed a 53.5% repeat content and 35,850 protein-coding genes. Chia shared a common ancestor with Salvia splendens ~6.1 million years ago. Utilizing the reference genome and two transcriptome datasets, we identified candidate fatty acid desaturases responsible for ALA biosynthesis during chia seed development. Because the seed of S. splendens contains significantly lower proportion of ALA but similar total contents of unsaturated fatty acids, we suggest that strong expression of two ShFAD3 genes are critical for the high ALA content of chia seeds. This genome assembly will serve as a valuable resource for breeding, comparative genomics, and functional genomics studies of chia.

Guard Cell Transcriptome Reveals Membrane Transport, Stomatal Development and Cell Wall Modifications as Key Traits Involved in Salinity Tolerance in Halophytic Chenopodium quinoa.

A comparative investigation was conducted to evaluate transcriptional changes in guard cells (GCs) of closely related halophytic (Chenopodium quinoa) and glycophytic (Spinacia oleracea) species. Plants were exposed to 3 weeks of 250 mM sodium chloride treatment, and GC-enriched epidermal fragments were mechanically prepared. In both species, salt-responsive genes were mainly related to categories of protein metabolism, secondary metabolites, signal transduction and transport systems. Genes related to abscisic acid (ABA) signaling and ABA biosynthesis were strongly induced in quinoa but not in spinach GCs. Also, expression of the genes encoding transporters of amino acids, proline, sugars, sucrose and potassium increased in quinoa GCs under salinity stress. Analysis of cell-wall-related genes suggests that genes involved in lignin synthesis (e.g. lignin biosynthesis LACCASE 4) were highly upregulated by salt in spinach GCs. In contrast, transcripts related to cell wall plasticity Pectin methylesterase3 (PME3) were highly induced in quinoa. Faster stomatal response to light and dark measured by observing kinetics of changes in stomatal conductance in quinoa might be associated with higher plasticity of the cell wall regulated by PME3 Furthermore, genes involved in the inhibition of stomatal development and differentiation were highly expressed by salt in quinoa, but not in spinach. These changes correlated with reduced stomatal density and index in quinoa, thus improving its water use efficiency. The fine modulation of transporters, cell wall modification and controlling stomatal development in GCs of quinoa may have resulted in high K+/Na+ ratio, lower stomatal conductance and higher stomatal speed for better adaptation to salinity stress in quinoa.

Recovery of chloroplast genomes from medieval millet grains excavated from the Areni-1 cave in southern Armenia.

Panicum miliaceum L. was domesticated in northern China at least 7000 years ago and was subsequentially adopted in many areas throughout Eurasia. One such locale is Areni-1 an archaeological cave site in Southern Armenia, where vast quantities archaeobotanical material were well preserved via desiccation. The rich botanical material found at Areni-1 includes P. miliaceum grains that were identified morphologically and14C dated to the medieval period (873 ± 36 CE and 1118 ± 35 CE). To investigate the demographic and evolutionary history of the Areni-1 millet, we used ancient DNA extraction, hybridization capture enrichment, and high throughput sequencing to assemble three chloroplast genomes from the medieval grains and then compared these sequences to 50 modern P. miliaceum chloroplast genomes. Overall, the chloroplast genomes contained a low amount of diversity with domesticated accessions separated by a maximum of 5 SNPs and little inference on demography could be made. However, in phylogenies the chloroplast genomes separated into two clades, similar to what has been reported for nuclear DNA from P. miliaceum. The chloroplast genomes of two wild (undomesticated) accessions of P. miliaceum contained a relatively large number of variants, 11 SNPs, not found in the domesticated accessions. These results demonstrate that P. miliaceum grains from archaeological sites can preserve DNA for at least 1000 years and serve as a genetic resource to study the domestication of this cereal crop.

Construction of a high-density bin-map and identification of fruit quality-related quantitative trait loci and functional genes in pear.

Pear (Pyrus spp.) is one of the most common fruit crops grown in temperate regions worldwide. Genetic enhancement of fruit quality is a fundamental goal of pear breeding programs. The genetic control of pear fruit quality traits is highly quantitative, and development of high-density genetic maps can facilitate fine-mapping of quantitative trait loci (QTLs) and gene identification. Bin-mapping is a powerful method of constructing high-resolution genetic maps from large-scale genotyping datasets. We performed whole-genome sequencing of pear cultivars 'Niitaka' and 'Hongxiangsu' and their 176 F 1 progeny to identify genome-wide single-nucleotide polymorphism (SNP) markers for constructing a high-density bin-map of pear. This analysis yielded a total of 1.93 million SNPs and a genetic bin-map of 3190 markers spanning 1358.5 cM, with an average adjacent interval of 0.43 cM. This bin-map, along with other high-density genetic maps in pear, improved the reference genome assembly from 75.5 to 83.7% by re-anchoring the scaffolds. A quantitative genetic analysis identified 148 QTLs for 18 fruit-related traits; among them, QTLs for stone cell content, several key monosaccharides, and fruit pulp acids were identified for the first time in pear. A gene expression analysis of six pear cultivars identified 399 candidates in the identified QTL regions, which showed expression specific to fruit developmental stages in pear. Finally, we confirmed the function of PbrtMT1, a tonoplast monosaccharide transporter-related gene responsible for the enhancement of fructose accumulation in pear fruit on linkage group 16, in a transient transformation experiment. This study provides genomic and genetic resources as well as potential candidate genes for fruit quality improvement in pear.

Genome-Wide Analysis of CqCrRLK1L and CqRALF Gene Families in Chenopodium quinoa and Their Roles in Salt Stress Response.

Chenopodium quinoa is a halophyte with exceptional nutritional qualities, and therefore it is potentially an ideal crop to grow in saline soils, not only addressing the problem of land salinization, but also providing nutrient food for the health of humans. Currently, the molecular mechanisms underlying salt tolerance in quinoa are still largely unknown. In Arabidopsis thaliana, Catharanthus roseus receptor-like kinase (CrRLK1Ls) FERONIA (FER) and its ligands rapid alkalinization factors (RALFs) have been reported that participate in the regulation of salt tolerance. Here, we performed a genome-wide analysis and identified 26 CqCrRLK1L and 18 CqRALF family genes in quinoa genome. Transcriptomic profiling of the leaf, root, stamen, and pistil tissues of quinoa reveals that different CqCrRLK1L and CqRALF genes exhibit tissue-specific expression patterns, which is consistent with that observed in other plant species. RNA-seq data show that three CqCrRLK1L genes are highly up-regulated after salt treatment, suggesting that some CqCrRLK1L family genes are transcriptionally responsive to salt stress in quinoa. Biochemical study indicates that CqRALF15, a paralog of Arabidopsis RALF22, is physically associated with CrRLK1L proteins CqFER and AtFER. CqRALF15 and AtRALF22 are functionally conserved in inducing the internalization of AtFER and in triggering root growth inhibition in both quinoa and Arabidopsis. Moreover, overexpression of CqRALF15 in Arabidopsis results in enhanced leaf bleaching under salt stress, indicating that CqRALF15 is involved in salt stress response. Together, our study characterizes CqCrRLK1L and CqRALF family genes in quinoa at genomic, transcriptional, and protein levels, and provides evidence to support their roles in salt stress response.

Genome-wide identification and comparative evolutionary analysis of sorbitol metabolism pathway genes in four Rosaceae species and three model plants.

In contrast to most land plant species, sorbitol, instead of sucrose, is the major photosynthetic product in many Rosaceae species. It has been well illustrated that three key functional genes encoding sorbitol-6-phosphate dehydrogenase (S6PDH), sorbitol dehydrogenase (SDH), and sorbitol transporter (SOT), are mainly responsible for the synthesis, degradation and transportation of sorbitol. In this study, the genome-wide identification of S6PDH, SDH and SOT genes was conducted in four Rosaceae species, peach, mei, apple and pear, and showed the sorbitol bio-pathway to be dominant (named sorbitol present group, SPG); another three related species, including tomato, poplar and Arabidopsis, showed a non-sorbitol bio-pathway (named sorbitol absent group, SAG). To understand the evolutionary differences of the three important gene families between SAG and SPG, their corresponding gene duplication, evolutionary rate, codon bias and positive selection patterns have been analyzed and compared. The sorbitol pathway genes in SPG were found to be expanded through dispersed and tandem gene duplications. Branch-specific model analyses revealed SDH and S6PDH clade A were under stronger purifying selection in SPG. A higher frequency of optimal codons was found in S6PDH and SDH than that of SOT in SPG, confirming the purifying selection effect on them. In addition, branch-site model analyses revealed SOT genes were under positive selection in SPG. Expression analyses showed diverse expression patterns of sorbitol-related genes. Overall, these findings provide new insights in the evolutionary characteristics for the three key sorbitol metabolism-related gene families in Rosaceae and other non-sorbitol dominant pathway species.

Liquid-liquid phase separation of RBGD2/4 is required for heat stress resistance in Arabidopsis.

As sessile organisms, plants are highly sensitive to environmental stresses. In response to stresses, globally repressed translation initiation leads to stress granule (SG) formation. Protein liquid-liquid phase separation (LLPS) contributes to SG formation, but a direct link between protein LLPS and stress resistance has not yet been found in plants. Here, we report that two RNA-binding proteins, RBGD2 and RBGD4, function redundantly to improve heat resistance in Arabidopsis. RBGD2 and RBGD4 undergo LLPS in vitro and condense into heat-induced SGs in vivo via tyrosine residue array (TRA). Importantly, disrupting LLPS by mutating TRA abolishes RBGD2/4 condensation in SGs and impairs their protective function against heat stress (HS). Further study found that upon HS, the RBGD2/4 interaction network expands with additional SG proteins and heat-responsive mRNA. Our work shows a mechanistic basis that underlies protein LLPS in HS response in plants and suggests manipulation of protein LLPS as a general strategy to improve plant stress resistance.

Salinity Effects on Guard Cell Proteome in Chenopodium quinoa.

Epidermal fragments enriched in guard cells (GCs) were isolated from the halophyte quinoa (Chenopodium quinoa Wild.) species, and the response at the proteome level was studied after salinity treatment of 300 mM NaCl for 3 weeks. In total, 2147 proteins were identified, of which 36% were differentially expressed in response to salinity stress in GCs. Up and downregulated proteins included signaling molecules, enzyme modulators, transcription factors and oxidoreductases. The most abundant proteins induced by salt treatment were desiccation-responsive protein 29B (50-fold), osmotin-like protein OSML13 (13-fold), polycystin-1, lipoxygenase, alpha-toxin, and triacylglycerol lipase (PLAT) domain-containing protein 3-like (eight-fold), and dehydrin early responsive to dehydration (ERD14) (eight-fold). Ten proteins related to the gene ontology term "response to ABA" were upregulated in quinoa GC; this included aspartic protease, phospholipase D and plastid-lipid-associated protein. Additionally, seven proteins in the sucrose-starch pathway were upregulated in the GC in response to salinity stress, and accumulation of tryptophan synthase and L-methionine synthase (enzymes involved in the amino acid biosynthesis) was observed. Exogenous application of sucrose and tryptophan, L-methionine resulted in reduction in stomatal aperture and conductance, which could be advantageous for plants under salt stress. Eight aspartic proteinase proteins were highly upregulated in GCs of quinoa, and exogenous application of pepstatin A (an inhibitor of aspartic proteinase) was accompanied by higher oxidative stress and extremely low stomatal aperture and conductance, suggesting a possible role of aspartic proteinase in mitigating oxidative stress induced by saline conditions.

Sugar Beet (Beta vulgaris) Guard Cells Responses to Salinity Stress: A Proteomic Analysis.

Soil salinity is a major environmental constraint affecting crop growth and threatening global food security. Plants adapt to salinity by optimizing the performance of stomata. Stomata are formed by two guard cells (GCs) that are morphologically and functionally distinct from the other leaf cells. These microscopic sphincters inserted into the wax-covered epidermis of the shoot balance CO2 intake for photosynthetic carbon gain and concomitant water loss. In order to better understand the molecular mechanisms underlying stomatal function under saline conditions, we used proteomics approach to study isolated GCs from the salt-tolerant sugar beet species. Of the 2088 proteins identified in sugar beet GCs, 82 were differentially regulated by salt treatment. According to bioinformatics analysis (GO enrichment analysis and protein classification), these proteins were involved in lipid metabolism, cell wall modification, ATP biosynthesis, and signaling. Among the significant differentially abundant proteins, several proteins classified as "stress proteins" were upregulated, including non-specific lipid transfer protein, chaperone proteins, heat shock proteins, inorganic pyrophosphatase 2, responsible for energized vacuole membrane for ion transportation. Moreover, several antioxidant enzymes (peroxide, superoxidase dismutase) were highly upregulated. Furthermore, cell wall proteins detected in GCs provided some evidence that GC walls were more flexible in response to salt stress. Proteins such as L-ascorbate oxidase that were constitutively high under both control and high salinity conditions may contribute to the ability of sugar beet GCs to adapt to salinity by mitigating salinity-induced oxidative stress.

Unbiased subgenome evolution following a recent whole-genome duplication in pear (Pyrus bretschneideri Rehd.).

Genome fractionation (also known as diploidization) frequently occurs following paleopolyploidization events. Biased fractionation between subgenomes has been found in some paleo-allopolyploids, while this phenomenon is absent in paleo-autopolyploids. Pear (Pyrus bretschneideri Rehd.) experienced a recent whole-genome duplication (WGD, ~30 million years ago); however, the evolutionary fate of the two subgenomes derived from this WGD event is not clear. In this study, we identified the two paleo-subgenomes in pear using peach (Prunus persica) as an outgroup and investigated differences in the gene loss rate, evolutionary rate, gene expression level, and DNA methylation level between these two subgenomes. Fractionation bias was not found between the two pear subgenomes, which evolved at similar evolutionary rates. The DNA methylation level of the two subgenomes showed little bias, and we found no expression dominance between the subgenomes. However, we found that singleton genes and homeologous genes within each subgenome showed divergent evolutionary patterns of selective constraints, expression and epigenetic modification. These results provide insights into subgenome evolution following paleopolyploidization in pear.

Gene duplication and evolution in recurring polyploidization-diploidization cycles in plants.

BACKGROUND: The sharp increase of plant genome and transcriptome data provide valuable resources to investigate evolutionary consequences of gene duplication in a range of taxa, and unravel common principles underlying duplicate gene retention. RESULTS: We survey 141 sequenced plant genomes to elucidate consequences of gene and genome duplication, processes central to the evolution of biodiversity. We develop a pipeline named DupGen_finder to identify different modes of gene duplication in plants. Genes derived from whole-genome, tandem, proximal, transposed, or dispersed duplication differ in abundance, selection pressure, expression divergence, and gene conversion rate among genomes. The number of WGD-derived duplicate genes decreases exponentially with increasing age of duplication events-transposed duplication- and dispersed duplication-derived genes declined in parallel. In contrast, the frequency of tandem and proximal duplications showed no significant decrease over time, providing a continuous supply of variants available for adaptation to continuously changing environments. Moreover, tandem and proximal duplicates experienced stronger selective pressure than genes formed by other modes and evolved toward biased functional roles involved in plant self-defense. The rate of gene conversion among WGD-derived gene pairs declined over time, peaking shortly after polyploidization. To provide a platform for accessing duplicated gene pairs in different plants, we constructed the Plant Duplicate Gene Database. CONCLUSIONS: We identify a comprehensive landscape of different modes of gene duplication across the plant kingdom by comparing 141 genomes, which provides a solid foundation for further investigation of the dynamic evolution of duplicate genes.

HDAC1 and HDAC2 Regulate Intermediate Progenitor Positioning to Safeguard Neocortical Development.

Neural progenitors with distinct potential to generate progeny are associated with a spatially distinct microenvironment. Neocortical intermediate progenitors (IPs) in the subventricular zone (SVZ) of the developing brain generate neurons for all cortical layers and are essential for cortical expansion. Here, we show that spatial control of IP positioning is essential for neocortical development. We demonstrate that HDAC1 and HDAC2 regulate the spatial positioning of IPs to form the SVZ. Developmental stage-specific depletion of both HDAC1 and HDAC2 in radial glial progenitors results in mispositioning of IPs at the ventricular surface, where they divide and differentiate into neurons, thereby leading to the cortical malformation. We further identified the proneural gene Neurogenin2 as a key target of HDAC1 and HDAC2 for regulating IP positioning. Our results demonstrate the importance of the spatial positioning of neural progenitors in cortical development and reveal a mechanism underlying the establishment of the SVZ microenvironment.

The genome of broomcorn millet.

Broomcorn millet (Panicum miliaceum L.) is the most water-efficient cereal and one of the earliest domesticated plants. Here we report its high-quality, chromosome-scale genome assembly using a combination of short-read sequencing, single-molecule real-time sequencing, Hi-C, and a high-density genetic map. Phylogenetic analyses reveal two sets of homologous chromosomes that may have merged ~5.6 million years ago, both of which exhibit strong synteny with other grass species. Broomcorn millet contains 55,930 protein-coding genes and 339 microRNA genes. We find Paniceae-specific expansion in several subfamilies of the BTB (broad complex/tramtrack/bric-a-brac) subunit of ubiquitin E3 ligases, suggesting enhanced regulation of protein dynamics may have contributed to the evolution of broomcorn millet. In addition, we identify the coexistence of all three C4 subtypes of carbon fixation candidate genes. The genome sequence is a valuable resource for breeders and will provide the foundation for studying the exceptional stress tolerance as well as C4 biology.

Recent polyploidization events in three Saccharum founding species.

The complexity of polyploid Saccharum genomes hindered progress of genome research and crop improvement in sugarcane. To understand their genome structure, transcriptomes of 59 F1 individuals derived from S. officinarumLA Purple and S. robustum Molokai 5829 (2n = 80, x = 10 for both) were sequenced, yielding 11 157 and 8998 SNPs and 83 and 105 linkage groups, respectively. Most markers in each linkage group aligned to single sorghum chromosome. However, 71 interchromosomal rearrangements were detected between sorghum and S. officinarum or S. robustum, and 24 (33.8%) of them were shared between S. officinarum and S. robustum, indicating their occurrence before the speciation event that separated these two species. More than 2000 gene pairs from S. spontaneum, S. officinarum and S. robustum were analysed to estimate their divergence time. Saccharum officinarum and S. robustum diverged about 385 thousand years ago, and the whole-genome duplication events occurred after the speciation event because of shared interchromosomal rearrangements. The ancestor of these two species diverged from S. spontaneum about 769 thousand years ago, and the reduction in basic chromosome number from 10 to 8 in S. spontaneum occurred after the speciation event but before the two rounds of whole-genome duplication. Our results proved that S. officinarum is a legitimate species in its own right and not a selection from S. robustum during the domestication process in the past 10 000 years. Our findings rejected a long-standing hypothesis and clarified the timing of speciation and whole-genome duplication events in Saccharum.

Diversification and independent domestication of Asian and European pears.

BACKGROUND: Pear (Pyrus) is a globally grown fruit, with thousands of cultivars in five domesticated species and dozens of wild species. However, little is known about the evolutionary history of these pear species and what has contributed to the distinct phenotypic traits between Asian pears and European pears. RESULTS: We report the genome resequencing of 113 pear accessions from worldwide collections, representing both cultivated and wild pear species. Based on 18,302,883 identified SNPs, we conduct phylogenetics, population structure, gene flow, and selective sweep analyses. Furthermore, we propose a model for the divergence, dissemination, and independent domestication of Asian and European pears in which pear, after originating in southwest China and then being disseminated throughout central Asia, has eventually spread to western Asia, and then on to Europe. We find evidence for rapid evolution and balancing selection for S-RNase genes that have contributed to the maintenance of self-incompatibility, thus promoting outcrossing and accounting for pear genome diversity across the Eurasian continent. In addition, separate selective sweep signatures between Asian pears and European pears, combined with co-localized QTLs and differentially expressed genes, underline distinct phenotypic fruit traits, including flesh texture, sugar, acidity, aroma, and stone cells. CONCLUSIONS: This study provides further clarification of the evolutionary history of pear along with independent domestication of Asian and European pears. Furthermore, it provides substantive and valuable genomic resources that will significantly advance pear improvement and molecular breeding efforts.

Different Modes of Gene Duplication Show Divergent Evolutionary Patterns and Contribute Differently to the Expansion of Gene Families Involved in Important Fruit Traits in Pear (Pyrus bretschneideri).

Pear is an important fruit crop of the Rosaceae family and has experienced two rounds of ancient whole-genome duplications (WGDs). However, whether different types of gene duplications evolved differently after duplication remains unclear in the pear genome. In this study, we identified the different modes of gene duplication in pear. Duplicate genes derived from WGD, tandem, proximal, retrotransposed, DNA-based transposed or dispersed duplications differ in genomic distribution, gene features, selection pressure, expression divergence, regulatory divergence and biological roles. Widespread sequence, expression and regulatory divergence have occurred between duplicate genes over the 30-45 million years of evolution after the recent genome duplication in pear. The retrotransposed genes show relatively higher expression and regulatory divergence than other gene duplication modes. In contrast, WGD genes underwent a slower sequence divergence and may be influenced by abundant gene conversion events. Moreover, the different classes of duplicate genes exhibited biased functional roles. We also investigated the evolution and expansion patterns of the gene families involved in sugar and organic acid metabolism pathways, which are closely related to the fruit quality and taste in pear. Single-gene duplications largely account for the extensive expansion of gene families involved in the sorbitol metabolism pathway in pear. Gene family expansion was also detected in the sucrose metabolism pathway and tricarboxylic acid cycle pathways. Thus, this study provides insights into the evolutionary fates of duplicated genes.

Genome-wide characterization, evolution, and expression analysis of the leucine-rich repeat receptor-like protein kinase (LRR-RLK) gene family in Rosaceae genomes.

BACKGROUND: Leucine-rich repeat receptor-like protein kinase (LRR-RLK) is the largest gene family of receptor-like protein kinases (RLKs) and actively participates in regulating the growth, development, signal transduction, immunity, and stress responses of plants. However, the patterns of LRR-RLK gene family evolution in the five main Rosaceae species for which genome sequences are available have not yet been reported. In this study, we performed a comprehensive analysis of LRR-RLK genes for five Rosaceae species: Fragaria vesca (strawberry), Malus domestica (apple), Pyrus bretschneideri (Chinese white pear), Prunus mume (mei), and Prunus persica (peach), which contained 201, 244, 427, 267, and 258 LRR-RLK genes, respectively. RESULTS: All LRR-RLK genes were further grouped into 23 subfamilies based on the hidden Markov models approach. RLK-Pelle_LRR-XII-1, RLK-Pelle_LRR-XI-1, and RLK-Pelle_LRR-III were the three largest subfamilies. Synteny analysis indicated that there were 236 tandem duplicated genes in the five Rosaceae species, among which subfamilies XII-1 (82 genes) and XI-1 (80 genes) comprised 68.6%. CONCLUSIONS: Our results indicate that tandem duplication made a large contribution to the expansion of the subfamilies. The gene expression, tissue-specific expression, and subcellular localization data revealed that LRR-RLK genes were differentially expressed in various organs and tissues, and the largest subfamily XI-1 was highly expressed in all five Rosaceae species, suggesting that LRR-RLKs play important roles in each stage of plant growth and development. Taken together, our results provide an overview of the LRR-RLK family in Rosaceae genomes and the basis for further functional studies.

Genome-wide identification of the MADS-box transcription factor family in pear (Pyrus bretschneideri) reveals evolution and functional divergence.

MADS-box transcription factors play significant roles in plant developmental processes such as floral organ conformation, flowering time, and fruit development. Pear (Pyrus), as the third-most crucial temperate fruit crop, has been fully sequenced. However, there is limited information about the MADS family and its functional divergence in pear. In this study, a total of 95 MADS-box genes were identified in the pear genome, and classified into two types by phylogenetic analysis. Type I MADS-box genes were divided into three subfamilies and type II genes into 14 subfamilies. Synteny analysis suggested that whole-genome duplications have played key roles in the expansion of the MADS family, followed by rearrangement events. Purifying selection was the primary force driving MADS-box gene evolution in pear, and one gene pairs presented three codon sites under positive selection. Full-scale expression information for PbrMADS genes in vegetative and reproductive organs was provided and proved by transcriptional and reverse transcription PCR analysis. Furthermore, the PbrMADS11(12) gene, together with partners PbMYB10 and PbbHLH3 was confirmed to activate the promoters of the structural genes in anthocyanin pathway of red pear through dual luciferase assay. In addition, the PbrMADS11 and PbrMADS12 were deduced involving in the regulation of anthocyanin synthesis response to light and temperature changes. These results provide a solid foundation for future functional analysis of PbrMADS genes in different biological processes, especially of pigmentation in pear.

Expansion and evolutionary patterns of cysteine-rich peptides in plants.

BACKGROUND: Cysteine-rich peptides (CRPs) are gaining recognition as regulators of cell-cell communication in plants. RESULTS: We identified 9556 CRPs in 12 plant species and analysed their evolutionary patterns. In most angiosperm plants, whole genome duplication and segmental duplication are the major factors driving the expansion of CRP family member genes, especially signal peptides. About 30% of the CRP genes were found clustered on the chromosomes, except in maize (Zea mays). Considerable collinearities between CRP genes between or within species reveal several syntenic regions on the chromosomes. Different subfamilies display diverse evolutionary rates, suggesting that these subfamilies are subjected to different selective pressures. CRPs in different duplication models also show contrasting evolutionary rates, although the underlying mechanism is unclear because of the complexity of gene evolution. The 1281 positively selected genes identified are probably generated within a certain period of time. While most of these belonged to maize and sorghum (Sorghum bicolor), new CRP functions would also be expected. Up-regulation of 10 CRPs was observed in self-pollinated pear pistils and pollen tubes under self S-RNase treatments in vitro. The expression divergence between different CRP gene duplication types suggests that different duplication mechanisms affected the fate of the duplicated CRPs. CONCLUSION: Our analyses of the evolution of the CRP gene family provides a unique view of the evolution of this large gene family.

Map-based cloning of the pear gene MYB114 identifies an interaction with other transcription factors to coordinately regulate fruit anthocyanin biosynthesis.

Red fruits are popular and widely accepted by consumers because of an enhanced appearance and enriched anthocyanins. The molecular mechanism of anthocyanin regulation in red-skinned pear (Pyrus) has been studied, and the genes encoding the biosynthetic steps and several transcription factors (TFs) have been characterized. In this study, a candidate R2R3 MYB TF, PyMYB114, was identified by linkage to the quantitative trait loci (QTL) for red skin color on linkage group 5 in a population of Chinese pear (Pyrus bretschneideri). The function of PyMYB114 was verified by transient transformation in tobacco (Nicotinana tabacum) leaves and strawberry (Fragaria) and pear fruits, resulting in the biosynthesis of anthocyanin. Suppression of PyMYB114 could inhibit anthocyanin biosynthesis in red-skinned pears. The ERF/AP2 TF PyERF3 was found to interact with PyMYB114 and its partner PybHLH3 to co-regulate anthocyanin biosynthesis, as shown by a dual luciferase reporter system and a yeast two-hybrid assay. In addition, the transcript abundance of PyMYB114 and PyMYB10 were correlated, and co-transformation of these two genes into tobacco and strawberry led to enhanced anthocyanin biosynthesis. This interaction network provides insight into the coloration of fruits and the interaction of different TFs to regulate anthocyanin biosynthesis.