Transposon insertions regulate genome-wide allele-specific expression and underpin flower colour variations in apple (Malus spp.).

Allele-specific expression (ASE) can lead to phenotypic diversity and evolution. However, the mechanisms regulating ASE are not well understood, particularly in woody perennial plants. In this study, we investigated ASE genes in the apple cultivar 'Royal Gala' (RG). A high quality chromosome-level genome was assembled using a homozygous tetra-haploid RG plant, derived from anther cultures. Using RNA-sequencing (RNA-seq) data from RG flower and fruit tissues, we identified 2091 ASE genes. Compared with the haploid genome of 'Golden Delicious' (GD), a parent of RG, we distinguished the genomic sequences between the two alleles of 817 ASE genes, and further identified allele-specific presence of a transposable element (TE) in the upstream region of 354 ASE genes. These included MYB110a that encodes a transcription factor regulating anthocyanin biosynthesis. Interestingly, another ASE gene, MYB10 also showed an allele-specific TE insertion and was identified using genome data of other apple cultivars. The presence of the TE insertion in both MYB genes was positively associated with ASE and anthocyanin accumulation in apple petals through analysis of 231 apple accessions, and thus underpins apple flower colour evolution. Our study demonstrated the importance of TEs in regulating ASE on a genome-wide scale and presents a novel method for rapid identification of ASE genes and their regulatory elements in plants.

microRNA172 targets APETALA2 to regulate flavonoid biosynthesis in apple (Malus domestica).

MicroRNA172 (miR172) plays a role in regulating a diverse range of plant developmental processes, including flowering, fruit development and nodulation. However, its role in regulating flavonoid biosynthesis is unclear. In this study, we show that transgenic apple plants over-expressing miR172 show a reduction in red coloration and anthocyanin accumulation in various tissue types. This reduction was consistent with decreased expression of APETALA2 homolog MdAP2_1a (a miR172 target gene), MdMYB10, and targets of MdMYB10, as demonstrated by both RNA-seq and qRT-PCR analyses. The positive role of MdAP2_1a in regulating anthocyanin biosynthesis was supported by the enhanced petal anthocyanin accumulation in transgenic tobacco plants overexpressing MdAP2_1a, and by the reduction in anthocyanin accumulation in apple and cherry fruits transfected with an MdAP2_1a virus-induced-gene-silencing construct. We demonstrated that MdAP2_1a could bind directly to the promoter and protein sequences of MdMYB10 in yeast and tobacco, and enhance MdMYB10 promotor activity. In Arabidopsis, over-expression of miR172 reduced flavonoid (including anthocyanins and flavonols) concentration and RNA transcript abundance of flavonoid genes in plantlets cultured on medium containing 7% sucrose. The anthocyanin content and RNA abundance of anthocyanin genes could be partially restored by using a synonymous mutant of MdAP2_1a, which had lost the miR172 target sequences at mRNA level, but not restored by using a WT MdAP2_1a. These results indicate that miR172 inhibits flavonoid biosynthesis through suppressing the expression of an AP2 transcription factor that positively regulates MdMYB10.

Genomic analyses provide insights into peach local adaptation and responses to climate change.

The environment has constantly shaped plant genomes, but the genetic bases underlying how plants adapt to environmental influences remain largely unknown. We constructed a high-density genomic variation map of 263 geographically representative peach landraces and wild relatives. A combination of whole-genome selection scans and genome-wide environmental association studies (GWEAS) was performed to reveal the genomic bases of peach adaptation to diverse climates. A total of 2092 selective sweeps that underlie local adaptation to both mild and extreme climates were identified, including 339 sweeps conferring genomic pattern of adaptation to high altitudes. Using genome-wide environmental association studies (GWEAS), a total of 2755 genomic loci strongly associated with 51 specific environmental variables were detected. The molecular mechanism underlying adaptive evolution of high drought, strong UVB, cold hardiness, sugar content, flesh color, and bloom date were revealed. Finally, based on 30 yr of observation, a candidate gene associated with bloom date advance, representing peach responses to global warming, was identified. Collectively, our study provides insights into molecular bases of how environments have shaped peach genomes by natural selection and adds candidate genes for future studies on evolutionary genetics, adaptation to climate changes, and breeding.

Identification of gene co-expression networks and key genes regulating flavonoid accumulation in apple (Malus × domestica) fruit skin.

Anthocyanin provides a red color for apple and health benefit for human. To better understand the molecular mechanisms of regulating apple color formation, we analyzed 27 transcriptomes of fruit skin from three cultivars 'Huashuo' (red-skinned), 'Hongcuibao' (red-skinned), and 'Golden Delicious' (yellow-skinned) at 0, 2, and 6 days after bag removal. Using pairwise comparisons and weighted gene co-expression network analyses (WGCNA), we constructed 17 co-expression modules. Among them, a specific module was negatively correlated to anthocyanin accumulation. The genes in the module are enriched in flavonoid biosynthesis pathways. These pathway genes were used to construct gene co-expression network of anthocyanin accumulation. Finally, a R2R3-MYB repressor designated MdMYB28 was identified as a key hub gene in the anthocyanin metabolism network. During the anthocyanin accumulation of apple fruit skin reaching a peak, MdMYB28 expression level was negatively correlated with the anthocyanin content. MdMYB28 was shown to directly bind to the promoter of MdMYB10 in yeast one-hybrid analyses. Over-expression of MdMYB28 decreased the anthocyanin biosynthesis in tobacco flower petals, suggesting that MdMYB28 acts as a negatively regulator of anthocyanin biosynthesis.

Transcriptome analysis of transgenic apple fruit overexpressing microRNA172 reveals candidate transcription factors regulating apple fruit development at early stages.

BACKGROUND: MicroRNA172 (miR172) has been proven to be critical for fruit growth, since elevated miR172 activity blocks the growth of apple (Malus x domestica Borkh.) fruit. However, it is not clear how overexpression of miR172 affects apple fruit developmental processes. METHODS: To answer this question, the present study, analyzed global transcriptional changes in miR172-overexpressing (miR172OX) and nongenetically modified wild-type (WT) apple fruit at two developmental stages and in different fruit tissues via RNA-seq. In addition, two cultivars, 'Hanfu' and 'M9', which have naturally fruit size variation, were included to identify miR172-dependent DEGs. qRT-PCRwas used to verify the reliability of our RNA-seq data. RESULTS: Overexpression of miR172 altered the expression levels of many cell proliferation- and cell expansion-related genes. Twenty-four libraries were generated, and 10,338 differentially expressed genes (DEGs) were detected between miR172OX and WT fruit tissues. 'Hanfu' and 'M9' are two common cultivars that bear fruit of different sizes (250 g and 75 g, respectively). Six libraries were generated, and 3,627 DEGs were detected between 'Hanfu' and 'M9'. After merging the two datasets, 6,888 candidate miR172-specific DEGs were identified. The potential networks associated with fruit size triggered traits were defined among genes belonging to the families of hormone synthesis, signaling pathways, and transcription factors. Our comparative transcriptome analysis provides insights into transcriptome responses to miR172 overexpression in apple fruit and a valuable database for future studies to validate functional genes and elucidate the fruit developmental mechanisms in apple.

Genomic analyses of an extensive collection of wild and cultivated accessions provide new insights into peach breeding history.

BACKGROUND: Human selection has a long history of transforming crop genomes. Peach (Prunus persica) has undergone more than 5000 years of domestication that led to remarkable changes in a series of agronomically important traits, but genetic bases underlying these changes and the effects of artificial selection on genomic diversity are not well understood. RESULTS: Here, we report a comprehensive analysis of peach evolution based on genome sequences of 480 wild and cultivated accessions. By focusing on a set of quantitative trait loci (QTLs), we provide evidence supporting that distinct phases of domestication and improvement have led to an increase in fruit size and taste and extended its geographic distribution. Fruit size was predominantly selected during domestication, and selection for large fruits has led to the loss of genetic diversity in several fruit weight QTLs. In contrast, fruit taste-related QTLs were successively selected for by domestication and improvement, with more QTLs selected for during improvement. Genome-wide association studies of 11 agronomic traits suggest a set of candidate genes controlling these traits and potential markers for molecular breeding. Candidate loci for genes that contributed to the adaption to low-chill regions were identified. Furthermore, the genomic bases of divergent selection for fruit texture and local breeding for different flavors between Asian and European/North American cultivars were also determined. CONCLUSIONS: Our results elucidate the genetic basis of peach evolution and provide new resources for future genomics-guided peach breeding.

Transcriptome analysis reveals novel genes involved in anthocyanin biosynthesis in the flesh of peach.

Peach is an important deciduous fruit tree species. Anthocyanins play an important role in fruit color formation and, through linkage analysis, previous studies have identified and mapped the key genes regulating anthocyanins' accumulation to chromosomes 3 and 5 in two different germplasms. To understand the overall regulatory network of anthocyanins biosynthesis, genes co-expressed with these key genes were identified in the red-fleshed 'Tianjin Shui Mi' and white-fleshed 'Hakuho' germplasms. Analysis of their flesh anthocyanin contents revealed differences 15 days before maturation. Therefore, transcriptome analysis of the flesh of fruits belonging to these two germplasms was performed to search for genes that were up-regulated at the late stage of development of 'Tianjin Shui Mi' but not of 'Hakuho', and identified 183 genes. These genes were also analyzed in the flesh transcriptomes of peach fruits belonging to 30 peach varieties with different anthocyanin contents at maturation, and the Pearson's correlation coefficients between their expression levels and anthocyanin contents were determined. The results showed that 66 genes were significantly correlated to anthocyanin contents, most of which previously reported as regulatory, biosynthetic, and transporter genes involved in anthocyanins' regulatory network. The results of this study enrich the understanding of key genes involved in the biological pathway regulating anthocyanins biosynthesis. The genes mostly associated with anthocyanins biosynthesis presented in this study are of great importance for molecular marker-assisted breeding.

Comparative Transcriptome and Microscopy Analyses Provide Insights into Flat Shape Formation in Peach (Prunus persica).

Fruit shape is an important external characteristic that consumers use to select preferred fruit cultivars. In peach, the flat fruit cultivars have become more and more popular worldwide. Genetic markers closely linking to the flat fruit trait have been identified and are useful for marker-assisted breeding. However, the cellular and genetic mechanisms underpinning flat fruit formation are still poorly understood. In this study, we have revealed the differences in fruit cell number, cell size, and in gene expression pattern between the traditional round fruit and modern flat fruit cultivars. Flat peach cultivars possessed significantly lower number of cells in the vertical axis because cell division in the vertical direction stopped early in the flat fruit cultivars at 15 DAFB (day after full bloom) than in round fruit cultivars at 35 DAFB. This resulted in the reduction in vertical development in the flat fruit. Significant linear relationship was observed between fruit vertical diameter and cell number in vertical axis for the four examined peach cultivars (R2 = 0.9964) at maturation stage, and was also observed between fruit vertical diameter and fruit weight (R2 = 0.9605), which indicated that cell number in vertical direction contributed to the flat shape formation. Furthermore, in RNA-seq analysis, 4165 differentially expressed genes (DEGs) were detected by comparing RNA-seq data between flat and round peach cultivars at different fruit development stages. In contrast to previous studies, we discovered 28 candidate genes potentially responsible for the flat shape formation, including 19 located in the mapping site and 9 downstream genes. Our study indicates that flat and round fruit shape in peach is primarily determined by the regulation of cell production in the vertical direction during early fruit development.