Pigmentation in sand pear (Pyrus pyrifolia) fruit: biochemical characterization, gene discovery and expression analysis with exocarp pigmentation mutant

A transcription factor, PbWRKY24, contributes to russet skin formation in pear fruits by modulating lignin accumulation

Skin color is one of the major traits of fruit appearance quality in pear (Pyrus) that affects the fruit commodity value. Russet skin protects pear fruits from environmental stresses and its formation process is closely linked to lignin accumulation. However, the molecular regulatory networks underlying russet skin formation in pear fruits involve complex secondary metabolic pathways and remain elusive. Here, we explored the regulatory mechanisms underlying lignin accumulation in pear skin based on transcriptome sequencing, co-expression network analysis, and gene expression profiling. We identified a WRKY transcription factor gene, PbWRKY24, that regulates russet skin formation in pear fruits. The relative expression of PbWRKY24 in russet pear skin was significantly correlated with lignin content. We then verified the function of PbWRKY24 in lignin accumulation via genetic transformation. DNA affinity purification sequencing revealed that PbWRKY24 directly binds to the promoter of a lignin biosynthesis gene, PbPRX4. This binding was confirmed by yeast one-hybrid, dual-luciferase, and electrophoretic mobility shift assays. Overexpression of PbPRX4 in pear skin stimulated lignin accumulation and consequently promoted russet skin formation. This study provides a glimpse into the intricate lignin biosynthesis mechanisms during russet skin formation in pear fruits, which is of practical significance to pear breeding for fruit quality.

Budding mutation reprogrammed flavonoid biosynthesis in jujube by deploying MYB41 and bHLH93

Budding mutations are known to cause metabolic changes in new jujube varieties; however, the mechanisms underlying these changes are still unclear. Here, we performed muti-omics analysis to decipher the detailed metabolic landscape of "Saimisu 1" (S1) and its budding mutation line "Saimisu 2" (S2) at all fruit stages. We found that the genes involved in the biosyntheses of flavonoids, phenylpropanoids, and amino acids were upregulated in S2 fruits at all stages, especially PAL and DFR, resulting in increased accumulation of related compounds in S2 mature fruits. Further co-expression regulatory network analysis showed that the transcription factors MYB41 and bHLH93 potentially regulated the expression of PAL and DFR, respectively, by directly binding to their promoters. Moreover, the overexpression of MYB41 or bHLH93 induced their expression levels to redirect the flux of the flavonoid biosynthetic pathway, eventually leading to high levels of related compounds in S2 fruits. Overall, this study revealed the metabolic variations between S1 and S2 and contributed to the understanding of the mechanisms underlying budding mutation-mediated metabolic variations in plants, eventually providing the basis for breeding excellent jujube varieties using budding mutation lines.

Identification and Characterization of ABCG15-A Gene Required for Exocarp Color Differentiation in Pear

Exocarp color is a commercially essential quality for pear which can be divided into two types: green and russet. The occurrence of russet color is associated with deficiencies and defects in the cuticular and epidermal layers, which affect the structure of the cell wall and the deposition of suberin. Until now, the genetic basics triggering this trait have not been well understood, and limited genes have been identified for the trait. To figure out the gene controlling the trait of exocarp color, we perform a comprehensive genome-wide association study, and we describe the candidate genes. One gene encoding the ABCG protein has been verified to be associated with the trait, using an integrative analysis of the metabolomic and transcriptomic data. This review covers a variety of omics resources, which provide a valuable resource for identifying gene-controlled traits of interest. The findings in this study help to elucidate the genetic components responsible for the trait of exocarp color in pear, and the implications of these findings for future pear breeding are evaluated.

QTL Analysis and CAPS Marker Development Linked with Russet in Pear (Pyrus spp.)

The fruit skin types of pear (Pyrus spp.) are divided into russet, smooth, and intermediate. One of the important traits in pear breeding programs is russet on pear fruit skin because it affects the commercial value. In the present study, a high-density genetic linkage map of ‘Whangkeumbae’ (smooth) × ‘Minibae’ (russet) was constructed. In addition, quantitative trait loci (QTL) analysis was performed to identify russet related QTL and develop a cleaved amplified polymorphism sequence (CAPS) marker. Together with SNPs derived from Axiom Pear 70K Genotyping Array and genotyping-by-sequencing derived SNPs and SSRs generated in previous study, an integrated genetic linkage map of ‘Whangkeumbae’ × ‘Minibae’ was constructed. A total of 1263 markers were anchored in 17 linkage groups (LGs) with a total genetic distance of 1894.02 cM and an average marker density of 1.48 cM. The chromosome coverage of ‘Whangkeumbae’ × ‘Minibae’ map was improved because the SNPs derived from Axiom Pear 70K Genotyping Array were anchored. QTL analysis was performed using previous russet phenotype data evaluated with russet coverage and Hunter a. As a result of QTL analysis, russet coverage- and Hunter a-related QTLs were identified in LG8 of the ‘Whangkeumbae’ × ‘Minibae’ map, and SNPs located in the QTL region were heterozygous in the ‘Minibae’. Although the russet coverage- and Hunter a-related QTLs were commonly detected in LG8, the logarithm of odds values of SNPs in the QTL region were higher in QTL related to russet coverage than to Hunter a. The CAPS marker (CBp08ca01) was developed using an array SNP located in the russet coverage related QTL, and the genotype of CBp08ca01 showed a 1:1 ratio in ‘Whangkeumbae’ × ‘Minibae’ (χ2 = 0.65, p > 0.05). ‘Whangkeumbae’ and ‘Minibae’ were thought to have rr and Rr genotypes, respectively, and the genetic factors controlling the russet formation might be located in chromosome 8. The CBp08ca01 was able to select F1 individuals with less than 30% russet coverage. Thus, it will be a useful tool for marker-assisted selection in pears.

The Making of Plant Armor: The Periderm

The periderm acts as armor protecting the plant's inner tissues from biotic and abiotic stress. It forms during the radial thickening of plant organs such as stems and roots and replaces the function of primary protective tissues such as the epidermis and the endodermis. A wound periderm also forms to heal and protect injured tissues. The periderm comprises a meristematic tissue called the phellogen, or cork cambium, and its derivatives: the lignosuberized phellem and the phelloderm. Research on the periderm has mainly focused on the chemical composition of the phellem due to its relevance as a raw material for industrial processes. Today, there is increasing interest in the regulatory network underlying periderm development as a novel breeding trait to improve plant resilience and to sequester CO₂. Here, we discuss our current understanding of periderm formation, focusing on aspects of periderm evolution, mechanisms of periderm ontogenesis, regulatory networks underlying phellogen initiation and cork differentiation, and future challenges of periderm research.

An ARF1-binding factor triggering programmed cell death and periderm development in pear russet fruit skin

Plants have a cuticular membrane (CM) and periderm membrane (PM), which act as barriers to terrestrial stresses. The CM covers primary organs with a continuous hydrophobic layer of waxes embedded in cutin, while the PM stacks with suberized cells outermost to the secondary tissues. The formation of native periderm is regulated by a postembryonic meristem phellogen that produces suberized phellem (cork) outwardly. However, the mechanism controlling phellogen differentiation to phellem remains to be clarified. Here, map-based cloning in a pear F1 population with segregation for periderm development in fruit skin facilitated the identification of an aspartic acid repeat deletion in Pyrus Periderm Programmed Cell Death 1.1 (PyPPCD1.1) that triggers phellogen activity for cork formation in pear russet fruit skin. PyPPCD1.1 showed preferential expression in pear fruit skin, and the encoded protein shares a structural similarity to that of the viral capsid proteins. Asp deletion in PyPPCD1.1 weakened its nuclear localization but increased its accumulation in the chloroplast. Both PyPPCD1.1 and its recessive allele directly interact with ADP-ribosylation factor 1 (ARF1). PyPPCD1.1 triggered PCD in an ARF1-dependent manner. Thus, this study identified the switch gene for PCD and periderm development and provided a new molecular regulatory mechanism underlying the development of this trait.

Dissecting the molecular mechanism of russeting in sand pear (Pyrus pyrifolia Nakai) by metabolomics, transcriptomics, and proteomics

Brown coloration and a rough appearance as russet and semi-russet (partial russet) are features unique to the popular Asian sand pear (Pyrus pyrifolia Nakai). The degree of russeting is different between different genotypes. Russeting is sensitive to water fluctuations, where excessive rainwater can trigger/stimulate its development. However, the molecular mechanism of russeting is currently unclear. Here, we employed multi-omics, i.e., metabolomics, transcriptomics, and proteomics, and analyzed the effect of different sand pear genotypes and artificial rainfall on russeting of pear fruits. This led to the identification of 79, 64, and 29 differentially produced/expressed metabolites, transcripts, and proteins that are involved in the biosynthesis of suberin, phenylpropane, cutin, and waxes. Further analysis of these differentially expressed genes and their encoded proteins revealed that four of them exhibited high expression at both transcript and protein levels. Transient expression of one such gene, PbHHT1 (accession number 103966555), which encodes ω-hydroxypalmitate-O-feruloyl transferase, in young green non-russet fruits triggered premature suberization in the russeting pear genotypes. This coincided with increased production of 16-feruloyloxypalmitic acid, a conjugated compound between phenols and esters during the polymerization for suberin formation. Collectively, our data from the combined three omics demonstrate that russeting in sand pear is a complex process involving the biosynthesis and transport of suberin and many other secondary metabolites.

Peridermal fruit skin formation in Actinidia sp. (kiwifruit) is associated with genetic loci controlling russeting and cuticle formation

BACKGROUND

The skin (exocarp) of fleshy fruit is hugely diverse across species. Most fruit types have a live epidermal skin covered by a layer of cuticle made up of cutin while a few create an outermost layer of dead cells (peridermal layer).

RESULTS

In this study we undertook crosses between epidermal and peridermal skinned kiwifruit, and showed that epidermal skin is a semi-dominant trait. Furthermore, backcrossing these epidermal skinned hybrids to a peridermal skinned fruit created a diverse range of phenotypes ranging from epidermal skinned fruit, through fruit with varying degrees of patches of periderm (russeting), to fruit with a complete periderm. Quantitative trait locus (QTL) analysis of this population suggested that periderm formation was associated with four loci. These QTLs were aligned either to ones associated with russet formation on chromosome 19 and 24, or cuticle integrity and coverage located on chromosomes 3, 11 and 24.

CONCLUSION

From the segregation of skin type and QTL analysis, it appears that skin development in kiwifruit is controlled by two competing factors, cuticle strength and propensity to russet. A strong cuticle will inhibit russeting while a strong propensity to russet can create a continuous dead skinned periderm.

An integrated metabolic and transcriptomic analysis reveals the mechanism through which fruit bagging alleviates exocarp semi-russeting in pear fruit

Fruit semi-russeting is an undesirable quality trait that occurs in fruit production. It is reported that preharvest fruit bagging could effectively alleviate fruit exocarp semi-russeting, but the physiological and molecular mechanisms remain unclear. In the present study, we performed an in-depth investigation into pear fruit semi-russeting from morphologic, metabolic and transcriptomic perspectives by comparing control (semi-russeted) and bagged (non-russeted) 'Cuiguan' pear fruits. The results showed that significant changes in cutin and suberin resulted in pear fruit semi-russeting. Compared with the skin of bagged fruits, the skin of the control fruits presented reduced cutin contents accompanied by an accumulation of suberin, which resulted in fruit semi-russeting; α, ω-dicarboxylic acids accounted for the largest proportion of typical suberin monomers. Moreover, combined transcriptomic and metabolic analysis revealed a series of genes involved in cutin and suberin biosynthesis, transport and polymerization differentially expressed between the two groups. Furthermore, the expression levels of genes involved in the stress response and in hormone biosynthesis and signaling were significantly altered in fruits with contrasting phenotypes. Finally, a number of transcription factors, including those of the MYB, NAC, bHLH and bZIP families, were differentially expressed. Taken together, the results suggest that the multilayered mechanism through which bagging alleviates pear fruit semi-russeting is complex, and the large number of candidate genes identified provides a good foundation for future functional studies.

Transcriptome and metabolomic analysis to reveal the browning spot formation of 'Huangguan' pear

BACKGROUND

Browning spot (BS) disorders seriously affect the appearance quality of 'Huangguan' pear and cause economic losses. Many studies on BS have mainly focused on physiological and biochemical aspects, and the molecular mechanism remains unclear.

RESULTS

In the present study, the structural characteristics of 'Huangguan' pear with BS were observed via scanning electron microscopy (SEM), the water loss and brown spots were evaluated, and transcriptomic and metabolomics analyses were conducted to reveal the molecular mechanism underlying 'Huangguan' pear skin browning disorder. The results showed that the occurrence of BS was accompanied by a decrease in the wax layer and an increase in lignified cells. Genes related to wax biosynthesis were downregulated in BS, resulting in a decrease in the wax layer in BS. Genes related to lignin were upregulated at the transcriptional level, resulting in upregulation of metabolites related to phenylpropanoid biosynthesis. Expression of calcium-related genes were upregulated in BS. Cold-induced genes may represent the key genes that induce the formation of BS. In addition, the results demonstrated that exogenous NaH₂PO₄·2H₂O and ABA treatment could inhibit the incidence of BS during harvest and storage time by increasing wax-related genes and calcium-related genes expression and increasing plant resistance, whereas the transcriptomics results indicated that GA₃ may accelerate the incidence and index of BS.

CONCLUSIONS

The results of this study indicate a molecular mechanism that could explain BS formation and elucidate the effects of different treatments on the incidence and molecular regulation of BS.

Yellow-fruited phenotype is caused by 573 bp insertion at 5' UTR of YFT1 allele in yft1 mutant tomato

The yft1 tomato mutant has a yellow-fruited phenotype controlled by a recessive gene of YFT1 allele, which has been shown by map-based cloning to be a homolog of ETHYLENE INSENSITIVE 2 (EIN2). Genetic lesion of YFT1 allele in yft1 is attributed to a 573 bp DNA fragment (IF₅₇₃) insertion at 1,200 bp downstream of the transcription start site. Transcriptomic analysis revealed that YFT1 lesion resulted in 5,053 differentially expressed genes (DEGs) in yft1 pericarp compared with the M82 wild type cultivar. These were annotated as being involved in ethylene synthesis, chromoplast development, and carotenoid synthesis. The YFT1 lesion caused a reduction in its own transcript levels in yft1 and impaired ethylene emission and signal transduction, delayed chromoplast development and decreased carotenoid accumulation. The molecular mechanism underlying the downregulated YFT1 allele in yft1 was examined at both RNA and DNA levels. The IF₅₇₃ event appeared to introduce two negative regulatory sequences located at -272 to -173 bp and -172 to -73 bp in the YFT1 allele promoter, causing alterative splicing due to introduction of aberrant splicing sites, and breaking upstream open reading frames (uORF) structure in the 5'-UTR. Those results a new provided insight into molecular regulation of color formation in tomato fruit.

Lignin and Quercetin Synthesis Underlies Berry Russeting in 'Sunshine Muscat' Grape

In order to further explore the mechanism of 'sunshine muscat' grape russet formation, transcriptomic and metabolomic analyses were performed on 'sunshine muscat' grape peels with and without russet. A total of 1491 differentially expressed genes (DEGs) were discovered based on these analyses. The phenylpropane synthesis pathway was the key metabolic pathway identified, and 28 DEGs related to phenylpropane synthesis pathway were screened, of which 16 were related to lignin synthesis. In addition, 60 differential metabolites were screened. There were 29 phenolic substances among the differential metabolites, which were all up-regulated and 10 were quercetin-related glycosides. Our results indicate that phenols likely play a dominant role in the formation of 'sunshine muscat' grape russet, and the synthesis of lignin and quercetin may be the key factors underlying russet formation.

Proteome and transcriptome profile analysis reveals regulatory and stress-responsive networks in the russet fruit skin of sand pear

The epidermal tissues of the cuticular membrane (CM) and periderm membrane (PM) confer first-line protection from environmental stresses in terrestrial plants. Although PM protection is essentially ubiquitous in plants, the protective mechanism, the function of many transcription factors and enzymes, and the genetic control of metabolic signaling pathways are poorly understood. Different microphenotypes and cellular components in russet (PM-covered) and green (CM-covered) fruit skins of pear were revealed by scanning and transmission electron microscopy. The two types of fruit skins showed distinct phytohormone accumulation, and different transcriptomic and proteomic profiles. The enriched pathways were detected by differentially expressed genes and proteins from the two omics analyses. A detailed analysis of the suberin biosynthesis pathways identified the regulatory signaling network, highlighting the general mechanisms required for periderm formation in russet fruit skin. The regulation of aquaporins at the protein level should play an important role in the specialized functions of russet fruit skin and PM-covered plant tissues.

Genome-wide association study and genomic prediction using parental and breeding populations of Japanese pear (Pyrus pyrifolia Nakai)

Breeding of fruit trees is hindered by their large size and long juvenile period. Genome-wide association study (GWAS) and genomic selection (GS) are promising methods for circumventing this hindrance, but preparing new large datasets for these methods may not always be practical. Here, we evaluated the potential of breeding populations evaluated routinely in breeding programs for GWAS and GS. We used a pear parental population of 86 varieties and breeding populations of 765 trees from 16 full-sib families, which were phenotyped for 18 traits and genotyped for 1,506 single nucleotide polymorphisms (SNPs). The power of GWAS and accuracy of genomic prediction were improved when we combined data from the breeding populations and the parental population. The accuracy of genomic prediction was improved further when full-sib data of the target family were available. The results suggest that phenotype data collected in breeding programs can be beneficial for GWAS and GS when they are combined with genome-wide marker data. The potential of GWAS and GS will be further extended if we can build a system for routine collection of the phenotype and marker genotype data for breeding populations.

A comparative transcriptomic approach to understanding the formation of cork

The transcriptome comparison of two oak species reveals possible candidates accounting for the exceptionally thick and pure cork oak phellem, such as those involved in secondary metabolism and phellogen activity. Cork oak, Quercus suber, differs from other Mediterranean oaks such as holm oak (Quercus ilex) by the thickness and organization of the external bark. While holm oak outer bark contains sequential periderms interspersed with dead secondary phloem (rhytidome), the cork oak outer bark only contains thick layers of phellem (cork rings) that accumulate until reaching a thickness that allows industrial uses. Here we compare the cork oak outer bark transcriptome with that of holm oak. Both transcriptomes present similitudes in their complexity, but whereas cork oak external bark is enriched with upregulated genes related to suberin, which is the main polymer responsible for the protective function of periderm, the upregulated categories of holm oak are enriched in abiotic stress and chromatin assembly. Concomitantly with the upregulation of suberin-related genes, there is also induction of regulatory and meristematic genes, whose predicted activities agree with the increased number of phellem layers found in the cork oak sample. Further transcript profiling among different cork oak tissues and conditions suggests that cork and wood share many regulatory mechanisms, probably reflecting similar ontogeny. Moreover, the analysis of transcripts accumulation during the cork growth season showed that most regulatory genes are upregulated early in the season when the cork cambium becomes active. Altogether our work provides the first transcriptome comparison between cork oak and holm oak outer bark, which unveils new regulatory candidate genes of phellem development.

Response of miR156-SPL Module during the Red Peel Coloration of Bagging-Treated Chinese Sand Pear (Pyrus pyrifolia Nakai)

MicroRNA156 is an evolutionarily highly conserved plant micro-RNA (miRNA) that controls an age-dependent flowering pathway. miR156 and its target SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) genes regulate anthocyanin accumulation in plants, but it is unknown whether this process is affected by light. Red Chinese sand pear (Pyrus pyrifolia) fruits exhibit a unique coloration pattern in response to bagging treatments, which makes them appropriate for studying the molecular mechanism underlying light-induced anthocyanin accumulation in fruit. Based on high-throughput miRNA and degradome sequencing data, we determined that miR156 was expressed in pear fruit peels, and targeted four SPL genes. Light-responsive elements were detected in the promoter regions of the miR156a and miR156ba precursors. We identified 19 SPL genes using the “Suli” pear (Pyrus pyrifolia Chinese White Pear Group) genome database, of which seven members were putative miR156 targets. The upregulated expression of anthocyanin biosynthetic and regulatory genes and downregulated expression of PpSPL2, PpSPL5, PpSPL7, PpSPL9, PpSPL10, PpSPL13, PpSPL16, PpSPL17, and PpSPL18 were observed in pear fruits after bags were removed from plants during the anthocyanin accumulation period. Additionally, miR156a/ba/g/s/sa abundance increased after bags were removed. Yeast two-hybrid results suggested that PpMYB10, PpbHLH, and PpWD40 could form a protein complex, probably involved in anthocyanin biosynthesis. Additionally, PpSPL10 and PpSPL13 interacted with PpMYB10. The results obtained in this study are helpful in understanding the possible role of miR156 and its target PpSPL genes in regulating light-induced red peel coloration and anthocyanin accumulation in pear.

The role of polyamines during exocarp formation in a russet mutant of 'Dangshansuli' pear (Pyrus bretschneideri Rehd.)

Differential genes of suberin, polyamine and transcription factors in transcriptome sequences and the contents of H 2 O 2 , spermidine, spermine, and putrescine changed significantly after treating with MGBG. Russeting is a commercially important process that restores the control of water loss through the skin via the formation of a waterproofing periderm just beneath the microcracked skin of pear primary fruit. A spontaneous russet skin mutant, the yellow-green 'Dangshansuli' pear, has been identified. To understand the role of polyamines in the formation of the russet skin of the mutant-type (MT) pear, it was treated with methylglyoxal-bis-(guanylhydrazone) (MGBG) for 4 weeks after full bloom. One week later, differentially expressed genes among the wild-type (WT), MT, and MGBG-treated MT pears were screened, hydrogen peroxide (H2O2) was localized using CeCl3, and the contents of H2O2 and polyamine were measured. A total of 57,086,772, 61,240,014, and 67,919,420 successful reads were generated from the transcriptomes of WT, MT, and MGBG-treated MT, with average unigene lengths of 701, 720, and 735 bp, respectively. Differentially expressed genes involved in polyamine metabolism and suberin synthesis were screened in 'Dangshansuli' and in the mutant libraries, and their relative expression was found to be significantly altered after treatment with MGBG, which was confirmed by real-time PCR. The expression patterns of differentially expressed transcription factors were identified and were found to be similar to those of the polyamine- and suberin-related genes. The results indicated that the H2O2 generated during polyamine metabolism might contribute to russet formation on the exocarp of the mutant pear. Furthermore, the contents of H2O2, spermidine, spermine, and putrescine and H2O2 localization provided a comprehensive transcriptomic view of russet formation in the mutant pear.

A major QTL controlling apple skin russeting maps on the linkage group 12 of 'Renetta Grigia di Torriana'

BACKGROUND

Russeting is a disorder developed by apple fruits that consists of cuticle cracking followed by the replacement of the epidermis by a corky layer that protects the fruit surface from water loss and pathogens. Although influenced by many environmental conditions and orchard management practices, russeting is under genetic control. The difficulty in classifying offspring and consequent variable segregation ratios have led several authors to conclude that more than one genetic determinant could be involved, although some evidence favours a major gene (Ru).

RESULTS

In this study we report the mapping of a major genetic russeting determinant on linkage group 12 of apple as inferred from the phenotypic observation in a segregating progeny derived from 'Renetta Grigia di Torriana', the construction of a 20 K Illumina SNP chip based genetic map, and QTL analysis. Recombination analysis in two mapping populations restricted the region of interest to approximately 400 Kb. Of the 58 genes predicted from the Golden Delicious sequence, a putative ABCG family transporter has been identified. Within a small set of russeted cultivars tested with markers of the region, only six showed the same haplotype of 'Renetta Grigia di Torriana'.

CONCLUSIONS

A major determinant (Ru_RGT) for russeting development putatively involved in cuticle organization is proposed as a candidate for controlling the trait. SNP and SSR markers tightly co-segregating with the Ru_RGT locus may assist the breeder selection. The observed segregations and the analysis of the 'Renetta Grigia di Torriana' haplotypic region in a panel of russeted and non-russeted cultivars may suggest the presence of other determinants for russeting in apple.

Apple russeting as seen through the RNA-seq lens: strong alterations in the exocarp cell wall

Russeting, a commercially important defect in the exocarp of apple (Malus × domestica), is mainly characterized by the accumulation of suberin on the inner part of the cell wall of the outer epidermal cell layers. However, knowledge on the underlying genetic components triggering this trait remains sketchy. Bulk transcriptomic profiling was performed on the exocarps of three russeted and three waxy apple varieties. This experimental design was chosen to lower the impact of genotype on the obtained results. Validation by qPCR was carried out on representative genes and additional varieties. Gene ontology enrichment revealed a repression of lignin and cuticle biosynthesis genes in russeted exocarps, concomitantly with an enhanced expression of suberin deposition, stress responsive, primary sensing, NAC and MYB-family transcription factors, and specific triterpene biosynthetic genes. Notably, a strong correlation (R(2) = 0.976) between the expression of a MYB93-like transcription factor and key suberin biosynthetic genes was found. Our results suggest that russeting is induced by a decreased expression of cuticle biosynthetic genes, leading to a stress response which not only affects suberin deposition, but also the entire structure of the cell wall. The large number of candidate genes identified in this study provides a solid foundation for further functional studies.