Early Response of Radish to Heat Stress by Strand-Specific Transcriptome and miRNA Analysis

Radish is a crucial vegetable crop of the Brassicaceae family with many varieties and large cultivated area in China. Radish is a cool season crop, and there are only a few heat tolerant radish varieties in practical production with little information concerning the related genes in response to heat stress. In this work, some physiological parameter changes of young leaves under short-term heat stress were detected. Furthermore, we acquired 1802 differentially expressed mRNAs (including encoding some heat shock proteins, heat shock factor and heat shock-related transcription factors), 169 differentially expressed lncRNAs and three differentially expressed circRNAs (novel_circ_0000265, novel_circ_0000325 and novel_circ_0000315) through strand-specific RNA sequencing technology. We also found 10 differentially expressed miRNAs (ath-miR159b-3p, athmiR159c, ath-miR398a-3p, athmiR398b-3p, ath-miR165a-5p, ath-miR169g-3p, novel_86, novel_107, novel_21 and ath-miR171b-3p) by small RNA sequencing technology. Through function prediction and enrichment analysis, our results suggested that the significantly possible pathways/complexes related to heat stress in radish leaves were circadian rhythm-plant, photosynthesis—antenna proteins, photosynthesis, carbon fixation in photosynthetic organisms, arginine and proline metabolism, oxidative phosphorylation, peroxisome and plant hormone signal transduction. Besides, we identified one lncRNA–miRNA–mRNAs combination responsive to heat stress. These results will be helpful for further illustration of molecular regulation networks of how radish responds to heat stress.


Introduction
Due to the significant increase in greenhouse gas emissions from human activities, especially the oxides of carbon dioxide, methane, chlorofluorocarbons and nitrogen, the global warming and the rising temperature have made plants in danger of high-temperature stress [1]. When the ambient temperature is higher than the suitable growth temperature of plants, some catalytic enzymes in the plant slowly lose their activity, resulting in structural damage of cells, which may lead to abnormal growth, flowering and seed yield reduction [2,3].
As sessile organisms, plants have evolved lots of complex regulatory mechanisms to cope with environmental stresses [4,5]. The excess production of singlet oxygen, superoxide anion radical, hydrogen peroxide and hydroxyl, also known as reactive oxygen species (ROS), are induced in plants under abiotic stress, which undermines chloroplasts and even results in cell death [6,7]. To scavenge superfluous ROS, plants have evolved antioxidative enzymes system, including superoxide dismutase (SOD) and peroxidases (POD). Overexpression of SOD in plants can protect the physiological processes of plants by clearing superoxide. The hydrogen peroxide could be eliminated by a variety of POD [8,9]. In addition, membrane lipid peroxidation occurs under heat stress (HS), and malondialdehyde (MDA)

Physiological Response to Heat Stress in Leaves Of The Cultivar "Huoche"
POD and SOD activity increased significantly after 6 h of heat stress, and still remained high after 3 d of heat stress. The heat-induced POD and SOD activity declined to control level after recovery (Figure 2A-B). Chlorophyll, soluble sugar and MDA contents were not affected by short-term heat stress, but decreased significantly after 3 d of heat stress, even during the recovery period ( Figure 2C-E). We observed that heat stress did not significantly change the content of free proline ( Figure 2F). These data suggest that heat stress reduces photosynthesis, but alleviates oxidative damage by increasing antioxidant enzymes in the cultivar "Huoche".

Physiological Response to Heat Stress in Leaves Of The Cultivar "Huoche"
POD and SOD activity increased significantly after 6 h of heat stress, and still remained high after 3 d of heat stress. The heat-induced POD and SOD activity declined to control level after recovery (Figure 2A,B). Chlorophyll, soluble sugar and MDA contents were not affected by short-term heat stress, but decreased significantly after 3 d of heat stress, even during the recovery period ( Figure 2C-E). We observed that heat stress did not significantly change the content of free proline ( Figure 2F). These data suggest that heat stress reduces photosynthesis, but alleviates oxidative damage by increasing antioxidant enzymes in the cultivar "Huoche".

Mapping and Quantitative Assessment of Illumina Sequence
To comprehensively evaluate how HS influenced transcript profile, six strand-specific RNA libraries (detecting mRNA, lncRNA and circRNA) and six small RNA libraries were constructed. For strand-specific RNA sequencing, the numbers of raw reads of these six samples range from 87.47 to 103.85 million, yielding 12.62 to 15.14 G clean bases for RNA-Seq with a Q20 percentage over 98.06%, Q30 percentage over 94.83%, and a GC percentage between 42.26% and 43.25% (Table S1). For small RNA sequencing, the numbers of raw reads of these six samples range from 13.58 to 14.90 million, yielding 0.679 to 0.745G clean bases for RNA-Seq with a Q20 percentage over 97.73%, Q30 percentage over 94.87%, and a GC percentage between 49.83% and 50.13% (Table S2).
Each library was aligned with the released reference genome as mentioned above. For strand-specific RNA sequencing, over 91.82% clean reads were successfully mapped to the reference genome, in which at least 84.51% and 4.39% clean reads were uniquely mapped and multiple mapped, respectively (Table S3). For small RNA sequencing, over 84.97% reads were successfully mapped to the reference genome, in which at least 61.07% and 22.42% reads were mapped to chains with the same direction of reference sequences and chains with the opposite direction of reference sequences, respectively (Table S4).

Analysis of DE mRNAs
In QHC06vsQHC00, 1802 differentially expressed (DE) mRNAs (1040 up-regulated and 762 down-regulated) were detected ( Figure 3A, Table S5). In the meanwhile, we constructed and analyzed the Venn diagram to compare overlapping relationships between two groups ( Figure 3B), and the results indicated that most mRNAs were shared by QHC00 and QHC06. The changes in the activities of peroxidases (POD) and superoxide dismutase (SOD), the contents of chlorophyll, soluble sugar, MDA and free proline of the samples in response to heat stress (HS), respectively. Bars with different letters above the columns of figures indicate significant differences at p less than 0.05 (lowercase letters) or less than 0.01 (capital letters).

Mapping and Quantitative Assessment of Illumina Sequence
To comprehensively evaluate how HS influenced transcript profile, six strand-specific RNA libraries (detecting mRNA, lncRNA and circRNA) and six small RNA libraries were constructed. For strand-specific RNA sequencing, the numbers of raw reads of these six samples range from 87.47 to 103.85 million, yielding 12.62 to 15.14 G clean bases for RNA-Seq with a Q20 percentage over 98.06%, Q30 percentage over 94.83%, and a GC percentage between 42.26% and 43.25% (Table S1). For small RNA sequencing, the numbers of raw reads of these six samples range from 13.58 to 14.90 million, yielding 0.679 to 0.745G clean bases for RNA-Seq with a Q20 percentage over 97.73%, Q30 percentage over 94.87%, and a GC percentage between 49.83% and 50.13% (Table S2).
Each library was aligned with the released reference genome as mentioned above. For strand-specific RNA sequencing, over 91.82% clean reads were successfully mapped to the reference genome, in which at least 84.51% and 4.39% clean reads were uniquely mapped and multiple mapped, respectively (Table S3). For small RNA sequencing, over 84.97% reads were successfully mapped to the reference genome, in which at least 61.07% and 22.42% reads were mapped to chains with the same direction of reference sequences and chains with the opposite direction of reference sequences, respectively (Table S4).

Analysis of DE mRNAs
In QHC06vsQHC00, 1802 differentially expressed (DE) mRNAs (1040 up-regulated and 762 down-regulated) were detected ( Figure 3A, Table S5). In the meanwhile, we constructed and analyzed the Venn diagram to compare overlapping relationships between two groups ( Figure 3B), and the results indicated that most mRNAs were shared by QHC00 and QHC06.

GO and KEGG Enrichment Analysis of DE mRNA Corresponding Genes under HS
Based on gene ontology (GO) analysis of DE mRNA corresponding genes, the significantly enriched biological processes (BPs) were metabolic process (GO: 0008152), cellular process (GO: 0009987), organic substance metabolic process (GO: 0071704) and cellular metabolic process (GO: 0044237), the most significantly enriched cellular components (CCs) were cell (GO: 0005623), cell part (GO: 0044464), intracellular (GO: 0005622) and intracellular part (GO: 0044424). Structural molecule activity (GO: 0005198) and cofactor binding (GO: 0048037) were significantly enriched within the molecular function (MF) categories ( Figure 4).

GO and KEGG Enrichment Analysis of DE mRNA Corresponding Genes under HS
Based on gene ontology (GO) analysis of DE mRNA corresponding genes, the significantly enriched biological processes (BPs) were metabolic process (GO: 0008152), cellular process (GO: 0009987), organic substance metabolic process (GO: 0071704) and cellular metabolic process (GO: 0044237), the most significantly enriched cellular components (CCs) were cell (GO: 0005623), cell part (GO: 0044464), intracellular (GO: 0005622) and intracellular part (GO: 0044424). Structural molecule activity (GO: 0005198) and cofactor binding (GO: 0048037) were significantly enriched within the molecular function (MF) categories ( Figure 4).

DE mRNAs Encoding Transcription Factor
In this study, a total of 165 DE mRNA were predicted to associate with 54 TF families (Table S7). The top four abundant types were the Orphans family with 13 DE mRNAs, the MYB family with 11 DE mRNAs, the AP2-EREBP family with 10 DE mRNAs, and the bZIP family with 10 DE mRNAs, respectively.

Analysis of DE lncRNAs, DE miRNAs and DE circRNAs
Compared with the control group, DE lncRNAs, miRNAs and circRNAs in the heat-treated group were displayed in the forms of a Volcano plot and Venn diagram ( Figure 9A-D and Figure 10A,B). The detailed data of the up-regulated and down-regulated lncRNAs between QHC06 and QHC00 are listed in Table S8. All DE miRNAs and circRNAs are listed in Tables 2 and 3, respectively. We identified 169 DE lncRNAs (117 up-regulated and 52 down-regulated), 10 DE miRNAs (two up-regulated and eight down-regulated) and three DE circRNAs (one up-regulated and two down-regulated) between QHC06 and QHC00, respectively.   We predicted the biological function of DE lncRNA through its co-location and co-expression with protein-coding genes. Figure 10C showed the Venn diagram of the intersection analysis between DE lncRNA targeted mRNA and DE mRNA. The numbers of up-regulated targeted mRNA of up-regulated lncRNAs, down-regulated targeted mRNA of down-regulated lncRNAs, up-regulated targeted mRNA of down-regulated lncRNAs and down-regulated targeted mRNA of up-regulated lncRNAs were 172, 103, 22 and 53, respectively (Table S9). Interestingly, a heat-related gene (HEAT repeat-containing protein 5B%2C) was also found in the list of up-regulated targeted mRNAs of up-regulated lncRNAs. We took another intersection analysis between DE miRNA targeted mRNA and DE mRNA based on psRNATarget online software and identified 18 up-regulated targeted mRNAs of down-regulated miRNAs and three down-regulated targeted mRNAs of up-regulated miRNAs under HS ( Figure 10D, Table 4). These results suggest that these mRNAs would be specially controlled by corresponding lncRNAs or miRNAs in the cultivar "Huoche" in response to heat stress.

Functional Prediction of DE lncRNA and DE miRNA in the Cultivar "Huoche"
On the basis of GO enrichment analysis of the targeted mRNA of DE lncRNAs, cellular process (GO: 0009987) and cellular metabolic process (GO: 0044237) were remarkably enriched among the BP categories. The most notably enriched CC was the cellular component (GO: 0005575). Among the CC categories, structural molecule activity (GO: 0005198) and structural constituent of ribosome (GO: 0003735) were significantly enriched ( Figure 11). On the basis of the GO enrichment analysis of the targeted mRNA of DE miRNAs, the most remarkably enriched BPs were negative regulation of apoptotic process (GO: 0043066), negative regulation of programmed cell death (GO: 0043069) and negative regulation of cell death (GO: 0060548). Anion binding (GO: 0043168), small molecule binding (GO: 0036094), nucleotide binding (GO: 0000166) and nucleoside phosphate binding (GO: 1901265) were significantly enriched among the MF categories ( Figure 12).
We utilized KEGG enrichment analysis to determine the most important biochemical metabolic pathways and signal transduction pathways involved in specific genes. The results showed that ribosome (ath03010) was the significantly enriched KEGG pathways of DE lncRNAs ( Figure 13A). The notably enriched KEGG pathways of DE miRNAs were propanoate metabolism (ath00640) and phagosome (ath04145; Figure 13B).
When comprehensively analyzing KEGG pathways among DE mRNA, DE lncRNA and DE miRNA, we found that their enriched pathways had high similarity. For instance, in the "circadian rhythm-plant" pathway, the gene in CCA1 node was enriched both in KEGG pathways of DE mRNA and DE miRNAs. The corresponding genes enriched in the "carbon fixation in photosynthetic organisms", "arginine and proline metabolism", "photosynthesis", "oxidative phosphorylation", "peroxisome" and "plant hormone signal transduction" of DE mRNA were almost enriched in the KEGG analysis of DE lncRNAs. Furthermore, the genes in EC:4.1.2.13 node and in EC:3.5.3.12 node were enriched among KEGG pathways of DE mRNA, DE lncRNA and DE miRNAs in the "carbon fixation in photosynthetic organisms" and "arginine and proline metabolism" pathways, respectively. All these results showed that these pathways might have something to do with thermotolerance of radish leaves.

Regulatory Network in Response to HS
After the construction of lncRNA-miRNA-mRNA combinations, we built a competing endogenous RNA (ceRNA) regulation network with XR_001945247.1 as a decoy, ath-miR165a-5p as a center, XM_018579231.1 and XM_018585227.1 as the target. Since a small number of DE circRNAs were detected, no ceRNA molecular network of circRNA-miRNA-mRNA was constructed.

Discussion
Heat stress, as a dominating global concern, has a profound effect on the crop's growth and developments, particularly in agronomic yield [21]. In this study, we found that most of the leaves of the cultivar the cultivar "Huoche" withered and died under the heat stress of 40 °C, and some individual plants could survive. Furthermore, we detected the changes of SOD, POD and MDA under heat stress. To identify HS-related genes, we analyzed the transcriptome of young leave of the cultivar "Huoche" under HS, 1802 DE mRNAs were gained. Furthermore, some lncRNAs and

Regulatory Network in Response to HS
After the construction of lncRNA-miRNA-mRNA combinations, we built a competing endogenous RNA (ceRNA) regulation network with XR_001945247.1 as a decoy, ath-miR165a-5p as a center, XM_018579231.1 and XM_018585227.1 as the target. Since a small number of DE circRNAs were detected, no ceRNA molecular network of circRNA-miRNA-mRNA was constructed.

Discussion
Heat stress, as a dominating global concern, has a profound effect on the crop's growth and developments, particularly in agronomic yield [21]. In this study, we found that most of the leaves of the cultivar the cultivar "Huoche" withered and died under the heat stress of 40 • C, and some individual plants could survive. Furthermore, we detected the changes of SOD, POD and MDA under heat stress. To identify HS-related genes, we analyzed the transcriptome of young leave of the cultivar "Huoche" under HS, 1802 DE mRNAs were gained. Furthermore, some lncRNAs and miRNAs were found to be significantly induced under HS, and a regulatory combination of lncRNA-miRNA-mRNA was suggested.

TFs Response to HS
Transcription factors (TFs) are a group of DNA-binding proteins that regulate the expression of the relative genes, which play a pivotal part in plant tolerance to abiotic stress [22,23].
Although there are little known about their specific functions, a few Orphans transcription factors may regulate the sensitivity of heat shock [24]. MYB participates in response to heat by regulating downstream relative genes in plants [25]. AP2-EREBPs not only function as key regulators in many developmental processes but also are associated with plant responses to abiotic environmental stresses [26,27]. In plants, bZIP participates in many processes, such as stress signaling and plant stress tolerance [28]. In this study, the Orphans, MYB, AP2-EREBP and bZIP families ranked in the top four among all TFs, indicating that they may play key roles in response to HS in radish (Table S6).
Phytochrome interacting factors (PIFs) with bHLH domain have been involved in the signaling mechanism of heat response [29]. Jain et al. [30] identified that HB genes played a critical role in the development and abiotic stress response in rice. Some TCP genes could be significantly up-regulated during the early duration under heat shock in soybean [31]. The expression of SbTCP19 was restricted by heat in Sorghum [32]. Zhang et al. [33] concluded that heat stress could induce the expression of OsMSR15 in rice, and two C2H2-type zinc finger motifs were the component of OsMSR15. Furthermore, heat stress was able to reduce the expression and enzymatic activity of cinnamate 3-hydroxylase (C3H) [34]. In our experiment, we found that bHLH, HB, TCP, C2H2 and C3H transcription factors under HS were significantly differentially regulated, suggesting their potential role of heat tolerance in the radish cultivar "Huoche" (Table S7).
NAM, ATAF and CUC (NAC) transcription factors, a large protein family, could function by helping regulate plant abiotic stress responses, and they could be applied to enhance stress tolerance in plants by genetic engineering [35]. The overexpression of NAC TF JUNGBRUNNEN1 (JUB1; ANAC042) could promote the production of heat shock proteins and improved the thermotolerance of Arabidopsis thaliana [36,37]. Similarly, two NAC genes were differentially expressed were detected in our study (Table S7). Chen et al. [38] found that OsMADS87 might help to improve the thermal resilience of rice. Duan et al. [39] identified that heat stress upregulated five BrMADS genes in the Chinese cabbage. We also acquired one up-regulated MADS gene (agamous-like MADS-box protein AGL19%2C) in this work (Table S7).

HS-Responsive HSPs and HSFs
Heat stress induces the expression of genes encoding heat shock proteins (HSPs) activated by heat shock factors (HSFs), which interact with heat shock elements in the promoter of HSP genes [40]. HSPs are a class of molecular chaperones involved in heat stress [41], including HSP100, HSP90, HSP70, HSP60 and small HSPs (sHSPs) [42]. HSP70 is involved in the feedback control of HS and relates to the activity of HSFA1a [43]. Wang et al. [16] identified four up-regulated and 13 down-regulated HSP70 transcripts under HS in radish taproots. In this study, after 6 h under HS, one up-regulated HSP70 component (hsp70 nucleotide exchange factor fes1) and one down-regulated HSP70 component (hsp70 nucleotide exchange factor fes1-like) were identified in radish leaves (Table 1).
HSFs play crucial roles in heat stress response by mediating the expression of HSPs [44]. On the basis of the structural features of the oligomerization domains, plant HSFs are classified into three categories, i.e., HSFA, HSFB and HSFC [43]. In tomato, HSFA1 serves as a master regulator of induced thermotolerance [45]. HSFA3 is controlled by the DREB2A transcription factor and is important for the establishment of heat tolerance in Arabidopsis [46]. TaHsfA6f is a transcriptional activator regulating some heat stress protection genes in wheat [47]. In the current study, HSFA-8-like was up-regulated after 6 h under HS in radish (Table 1). Taken together, HSPs and HSFs play a key role in the response to high temperature, and the further research of radish HSFA-8-like is valuable.

HS-Induced miRNAs
In recent years, more and more ncRNAs have emerged as key regulatory molecules in response to high temperature, whose regulatory mechanisms have been revealed in some plants [48]. miRNAs that were predominantly 20-24 nucleotides function by silencing of target mRNAs [49].
Both Xin et al. [50] and Hivrale et al. [51] found that miR159 was up-regulated in response to HS. The overexpression of miR159 could inhibit MYB transcripts to reduce plant thermotolerance [48,52]. In our study, ath-miR159b-3p and ath-miR159c were also significantly induced under HS ( Table 2). We also found that ath-miR169g-3p and ath-miR171b-3p were down-regulated in radish leaves under HS, similar to the results in Populus tomentosa [53,54], but miR169 and miR171 were up-regulated in response to HS in Arabidopsis thaliana [55,56], suggesting different regulation networks of these two small RNA under HS between Arabidopsis and radish. The miRNA miR398 is important for response to different abiotic stresses, especially heat stress. There are two essential target genes of miR398 in Arabidopsis, CSDs (closely related Cu/Zn-SODs) and CCS1 (the copper chaperone for SOD) [57,58]. The CSDs could scavenge the superoxide radicals, and its expression is fine-tuned by the cleavage of miR398-directed mRNA. Under HS, CSD and CCS relate to the synthesis of HSF and HSP. Some studies have indicated that down-regulation of miR398 could up-regulate CSDs expression, which promotes the accumulation of HSF and HSP, helping plants resist heat stress [56,59]. In the present study, ath-miR398b-3p and ath-miR398a-3p were significantly down-regulated in radish leaves under HS (Table 2), similar to that in Brassica rapa and Populus tomentosa [54,60]. In addition, we identified that SOD activity in the radish was increased immediately after heat treatment, these results indicated that the miR398-CSD/CCS pathway could play a key role in response to HS in radish.
Based on the above data, there may be such an assumption for the response of radish leaf to heat stress: The chlorophyll content of radish leaves decreased under heat stress (Figure 2), and the expression of some genes in photosynthesis pathway and light capture protein complex was affected (Figures 7 and 8). At the same time, circadian rhythm-plant pathway was significantly affected by heat stress, and this pathway had a direct impact on the light capture protein complex ( Figure 6). Antioxidant enzymes POD and SOD increased rapidly, and membrane lipid peroxidation marker MDA decreased (Figure 2), which alleviated oxidative stress. Correspondingly, some genes of peroxisome were affected ( Figure S4). In addition, plant hormone signal transduction was also significantly affected by high temperature, especially the up-regulation of ABF in ABA transduction pathway ( Figure S5), which can regulate stomatal closure and alleviate heat stress. These results provided basic data for further clarifying how radish seedlings respond to heat stress at the molecular level.

Plant Materials and HS Treatment
The seeds of the radish cultivar "Huoche" were soaked in water at 25 • C for three days. Germinated seeds were individually sown in plastic pots (35 cm × 22 cm) containing nutrient enriched peat soil in climate chamber with 16 h light at 25 • C and 8 h dark at 15 • C, 70% humidity and 4000 l× light (PQX-330B-30HM, Ningbo Life Technology Co. Ltd., Ningbo, China). Twenty-day-old seedlings were treated by heat shock at 40 • C.
After heat treatment of 0 h (Control) and 6 h, young leaves (the second and third leaf from the top) were collected for strand-specific RNA sequencing and small RNA sequencing (labeling QHC00 and QHC06). For physiological parameter measurement, radish seedlings were heat treated (16 h light at 40 • C) and 8 h dark at 30 • C) for 3 d, and then they were treated by normal temperature (16 h light at 25 • C and 8 h dark at 15 • C) for 3 d. Three biological repeats for RNA sequencing and physiological analysis were performed, and the collected samples were frozen immediately in liquid nitrogen and stored at −80 • C for further use. QHC001, QHC002 and QHC003 represent the three biological replicates of QHC00. QHC061, QHC062 and QHC063 represent the three biological replicates of QHC06.

Morphological and Physiological Analysis
According to the manufacturer's protocols (Nanjing Jiancheng Bioengineering Institute, Nanjing, China), the activities of POD, SOD and the contents of free proline, chlorophyll, soluble sugar and MDA of the samples were detected and determined using peroxidase assay kit, copper-zinc superoxide dismutase (CuZn-SOD) assay kit, proline assay kit, chlorophyll assay kit, plant soluble sugar content test kit, malondialdehyde (MDA) assay kit (TBA method), respectively.

RNA Isolation and RNA-Seq
The frozen samples were sent to Novogene Bioinformatics Technology Co. Ltd. (Beijing, China). Total RNA of each sample was isolated with Trizol reagents under the manufacturer's instruction (Thermo Fisher Scientific, Shanghai, China). Afterward, RNA quantification and qualification were checked with the methods described by Zhang et al. [61].

Differential Expression Analysis
DE mRNAs and DE lncRNAs were determined with the aid of the Ballgown [65]. DESeq2 was used to identified DE miRNAs and DE circRNAs [66].

GO and KEGG Enrichment Analysis
Gene ontology (GO) annotations and functional enrichment analysis of DE mRNA corresponding genes and DE ncRNA targeted genes were performed on the GO seq R package. GO terms are distributed into biological processes (BP), cellular components (CC) and molecular functions (MF) [67].
KEGG (Kyoto Encyclopedia of Genes and Genome) is the main public database related to the pathway. It is a systematic analysis of gene function and genomic information database, which is helpful for studying genes and expressions as a whole network (https://www.genome.jp/kegg/; 20-04-2018) [68]. By means of KOBAS v2.0 software [69], we detected the statistic enrichment of DE mRNA corresponding genes and DE ncRNA targeted genes in KEGG pathways.

Analysis of Transcription Factors
Transcription factors (TFs) are a group of proteins that can specifically bind to specific sequences upstream of the 5 end of the gene, thereby ensuring that the target gene is expressed at a specific time and space at a specific intensity. Plant transcription factor prediction was implemented by iTALK Light-harvesting complex I chlorophyll a/b binding protein 1 Lhcb1 Light-harvesting complex II chlorophyll a/b binding protein 1 Lhcb2 Light-harvesting complex II chlorophyll a/b binding protein 2 Lhcb4 Light-harvesting complex II chlorophyll a/b binding protein 4 Lhcb7 Light