Plant material
‘Whangkeumbae’ pear fruit were picked in orchards located in Wulong and Laiyang,
Shandong province, People’s Republic of China. Hard-end fruit were picked from ten-year-old
‘Whangkeumbae’ pear trees in one orchard, and normal pears were harvested from healthy
trees in another orchard. Normal and hard-end fruit were sampled at 60, 90 and 120
days after anthesis, and sampled again at 0, 60 and 120 days after harvest when storaged
under 0 °C. Three biological replicates comprised of ten fruits each were used for
each condition (normal vs. hard-end) and at each sampling timepoint. The fruit tissues
near the bottom third of the calyx end were taken. After removal of the peel and/or
seed, fleshy tissues were sliced into small pieces (approximately 1 cm3) and immediately frozen in liquid nitrogen. Sample were stored at −70 °C until further
analysis.
RNA-seq analysis
The calyx pulp of normal and hard-end fruit at 120 days after anthesis were RNA-seq
analysis. Total RNA was extracted using an RNA extraction kit (Omega, Georgia, USA)
according to the manufacturer’s instructions. The integrity and quality of the total
RNA was evaluated using a 2100 Bioanalyzer RNA Nano chip device (Agilent, Santa Clara,
CA, USA). The poly A - mRNA fraction was enriched by treatment of the extracted RNA
with oligo (dT) beads and was then reverse-transcribed into first strand cDNA for
use in preparation of the sequencing libraries.
The cDNA libraries were sequenced using an Illumina HiSeq 2500 system at the Biomarker
Technologies Corporation (Beijing, China). Raw reads were first filtered to remove
adaptors and low quality sequences, and then mapped to the pear reference genome (https://www.rosaceae.org/species/pyrus/pyrus_communis/genome_v1.0) using TopHat software. A false discovery rate (FDR) < 0.01 and a fold change of
⩾2 were used to identify differentially expressed genes (DEGs). The predicted product
of each unigene sequence was aligned to a set of proteins retrieved from the NCBI
Nr, Swiss-Prot, Kyoto Encyclopedia of Genes and Genomes (KEGG), and Cluster of Orthologous
Groups of proteins (COG) databases. The Reads Per Kb per Million Fragments (RPKM)
was used to determine the expression level of genes. The total number of reads for
each unigene and gene length were normalized by RPKM. The formula used to calculate
was as follows: RPKM = total exon reads / (mapped reads (millions)×exon length (KB)).
KEGG pathway enrichment analysis was performed using KOBAS software and utilized an
adjusted P-value of <0.05. Transcription factors were identified and classified into
different families by reference to the NCBI Nr, Swiss-Prot, and COG databases. The
raw sequences generated for ‘Whangkeumbae’ in this study were deposited in NCBI (NCBI
BioProject Accession: SRP063324, http://www.ncbi.nlm.nih.gov/bioproject/PRJNA294723).
Reverse transcription-quantitative polymerase chain reaction (RT-qPCR)
Total RNA was extracted from pear flesh tissue using RNAplant Reagent (TianGen, Shanghai,
China) according to the manufacturer’s instructions. Tobacco leaf RNA was extracted
using an EASYspin Plant RNA Kit (Yuanpinghao, China) and genomic DNA was removed by
treatments with DNase (Fermentas, Vilnius, Lithuania). The cDNA was synthesized by
reverse transcription using the Prime Script™ RT reagent Kit (Takara, Dalian, China)
according to the manufacturer's instructions and was subsequently used as template
in the RT-qPCR analyses. RT-qPCR was performed on a Light Cycler® 480 instrument (Roche,
Switzerland). The protocol included annealing at 94 °C for 5 min, followed by 40 cycles
of 94 °C for 15 s, and 60 °C for 1 min. Actin genes from pear and tobacco were used
for the normalization of transcript levels. Gene-specific primers used in the RT-qPCR
analyses were designed with Primer 3 (http://bioinfo.ut.ee/primer3-0.4.0/) software and are listed in Additional File 1: Table S1. Mean expression level was calculated using the 2-ΔΔCt method [39]. The expression level in normal fruit at 60 days after anthesis was set
as 1 in the RT-qPCR analyses conducted on samples collected during fruit development,
and the day of harvest was set as 1 in the post-harvest analyses. Three biological
and three technical replicates were used in the RT-qPCR analysis of each gene at each
timepoint.
Cloning of PpNAC187
Total RNA isolation and cDNA synthesis followed the same protocol used in the RT-qPCR
analyses. The PCR primers used to clone PpNAC187 are shown in Additional File 1: Table S2. The PCR program was: 94 ºC for 5 min, 35 cycles of 94 ºC for 30 s, 60 ºC for 1 min,
and 72 ºC for 1 min, followed by an extension cycle at 72°C for 10 min and a final cycle
at 4 ºC. PCR products were cloned into PMD19-T vectors (Takara, Dalian, China). The open reading
frame (ORF) of PpNAC187 was amplified and cloned using Phusion® High-Fidelity DNA Polymerase (Thermo scientific,
Lithuania, EU).
Sequence alignment and phylogenetic analysis
The amino acid sequence alignment analysis of NACs was conducted using DNAMAN software.
A phylogenetic tree was reconstructed with Figtree (http://tree.bio.ed.ac.uk/software/figtree/) online software. The amino acid sequence alignment analysis of pear NACs were refered to Ahmad [20].
Construction of the expression vector
The ORF of PpNAC187 was ligated into the expression vector, pCambia1300, under the control of a 35S promoter.
The ORF fragment isolated by digestion with KpnI and HindⅢ was inserted into the expression
vector, pSuper1300, under the control of a 35S promoter. The vectors, pCambia1300-PpNAC187
and pSuper1300-PpNAC187, were transferred into Agrobacterium tumefaciens EHA105 using the freeze-thaw method [21]. The sequences of primers used to construct the
expression vector are listed in Additional File 1: Table S2.
Subcellular localization of the PpNAC187 transcription factor
The subcellular localization of gene expression was determined using the method described
by Sun with some modifications [22]. After incubation for 24 h at 28 °C in the dark,
fresh onion scales (1.5×1 cm) were placed on a 9 cm plate with their inner surface
submerged in a 10 mL Agrobacterium solution (OD600 = 0.6-0.8) supplemented with 20 mg acetosyringone/L for 15-20 min. The onion scales
were then transferred to a 1/2 MS solid medium amended with 20 mg acetosyringone/L
and cultured for 16-24 days at 28 °C. The onion scales were subsequently rinsed with
water and epidermal cell layers were peeled and directly transferred to glass slides.
Agrobacteria harboring the pCambia1300-PpNAC187 or the empty pCambia1300 vector were
used in the analysis of subcellular localization. The GFP of onion scales inoculated
with these vectors were observed under a confocal laser scanning microscope (TCSSP5Ⅱ,
Leica, Germany).
Transient expression of PpNAC187 in ‘Whangkeumbae’ pear
The method of transient expression of PpNAC187 in pear ‘Whangkeumbae’ followed the method described by Spolaore with some modifications
[23]. Holes were punched on the calyx end of hard-end fruits on the harvest day using
a sterile syringe needle. One ml of Agrobacterium solution (OD600 = 0.6-0.8) was then injected into the fruit via the holes using a syringe without
a needle and the injected fruit was stored in the dark. Fruits inoculated with pSuper1300-PpNAC187
(treated) or the empty pSuper1300 vector (control) were photographed at the sampled timepoints.
Samples were taken at 1, 3, 5 and 10 days after the injection, and were immediately
frozen in liquid nitrogen and stored at −70 °C until further processing.
Agrobacterium-mediated transformation of tobacco with PpNAC187
The empty pSuper1300 vector and pSuper1300-PpNAC187 were independently transformed
into tobacco plants using the Agrobacterium-mediated transformation method as described
by Zheng with some modifications [24]. Portions of tobacco leaves without veins were
cut into discs (1×1 cm) and pre-cultured on MS solid medium for 2 days at 28 °C in
the dark. The leaf discs were then submerged in 15 mL of Agrobacterium solution (OD600 = 0.6-0.8) supplemented with 20 mg acetosyringone/L for 15-20 min. Transgenic tobacco
plants were generated on selection media after a 24 h light treatment following the
method of Wang [25].
PCR verification of transformed PpNAC187 tobacco
DNA was extracted from tobacco leaf tissue using DNAplant Reagent (TianGen, Shanghai,
China) according to the manufacturer’s instructions. PpNAC187 primers listed in Additional File 1: Table S2 were used to verify the presence of PpNAC187. The PCR program utilized was: 94 ºC for 5 min, 35 cycles of 94 ºC for 30 s, 60 ºC
for 1 min, and 72 ºC for 1 min, followed by a 10 min extension at 72°C and a final
cycle at 4 ºC.
Wiesner staining and microscopy
Wiesner reagent (phloroglucinol/HCl) staining of plant tissue for 5 min was used to
visualize lignification [5]. Two grams of phloroglucinol were dissolved in 100 ml
of 95% alcohol and then filtered into 40 ml of concentrated hydrochloric acid. A razor
blade was used to dissect leaf tissue prior to observation. Lignified structures appeared
pink or fuchsia in color. Auto-fluorescence within stem sections was observed with
the aid of an EVOS smart fluorescence microscope (Thermo Fisher, America).
Lignin assay
Lignin content was assayed using the method described by Dyckmans with some modifications [26]. Samples were washed three times in a 10 ml solution (100 mM K2HPO4/KH2PO4,
0.5%Triton X-100, 0.5% PVP, PH 7.8), followed by an additional three washes in 100%
methanol. The samples of fruit tissues were then dried overnight and tissue samples
were then transferred into 1ml of solution composed of 2 M HCl and 0.1ml thioglycolic
acid. Lignin was extracted in this solution by placing samples in a boiling water
bath for 4 h. Pellets obtained by centrifugation were resuspended 2 ml 1M NaOH followed
by agitation for 18 h. After the addition of 0.2 ml HCl, the mixture was incubated
for 4 h at 4 °C. The end product was dissolved in 1 ml 1M NaOH and absorbance at 280
nm was recorded to estimate lignin content. All measurements were performed in triplicate.
Statistical analyses
Two-tailed t-test were performed to determine the statistical significance of differences between samples.
Figures were drawn using Origin 6.0 (Microcal Software Inc. Northampton, MA, USA).