Genomics analysis of the light-harvesting chlorophyll a/b-binding (Lhc) superfamily in cassava (Manihot esculenta Crantz)

The light-harvesting chlorophyll a/b-binding proteins of photosystem II family members are responsible for temperature sensitivity and leaf color phenotype in albino tea plant

INTRODUCTION

Light-harvesting chlorophyll a/b-binding (LHCB) protein complexes of photosystem II are integral to the formation of thylakoid structure and the photosynthetic process. They play an important role in photoprotection, a crucial process in leaf development under low-temperature stress. Nonetheless, potential key genes directly related to low-temperature response and albino phenotype have not been precisely identified in tea plant. Moreover, there are no studies simultaneously investigating multiple albino tea cultivars with different temperature sensitivity.

OBJECTIVES

The study aimed to clarify the basic characteristics of CsLHCB gene family members, and identify critical CsLHCB genes potentially influential in leaf color phenotypic variation and low-temperature stress response by contrasting green and albino tea cultivars. Concurrently, exploring the differential expression of the CsLHCB gene family across diverse temperature-sensitive albino tea cultivars.

METHODS

We identified 20 putative CsLHCB genes according to phylogenetic analysis. Evolutionary relationships, gene duplication, chromosomal localization, and structures were analyzed by TBtools; the physiological and biochemical characteristics were analyzed by protein analysis websites; the differences in coding sequences and protein accumulation in green and albino tea cultivars, gene expression with maturity were tested by molecular biology technology; and protein interaction was analyzed in the STRING database.

RESULTS

All genes were categorized into seven groups, mapping onto 7 chromosomes, including three tandem and one segmental duplications. They all own a conserved chlorophyll A/B binding protein domain. The expression of CsLHCB genes was tissue-specific, predominantly in leaves. CsLHCB5 may play a key role in the process of leaf maturation and senescence. In contrast to CsLHCB5, CsLHCB1.1, CsLHCB2, and CsLHCB3.2 were highly conserved in amino acid sequence between green and albino tea cultivars. In albino tea cultivars, unlike in green cultivars, the expression of CsLHCB1.1, CsLHCB1.2, and CsLHCB2 was down-regulated under low-temperature stress. The accumulation of CsLHCB1 and CsLHCB5 proteins was lower in albino tea cultivars. Greater accumulation of CsLHCB2 protein was detected in RX1 and RX2 compared to other albino cultivars.

CONCLUSIONS

CsLHCB1.1, CsLHCB1.2, and CsLHCB2 played a role in the response to low-temperature stress. The amino acid sequence site mutation of CsLHCB5 would distinguish the green and albino tea cultivars. The less accumulation of CsLHCB1 and CsLHCB5 had a potential influence on albino leaves. Albino cultivars more sensitive to temperature exhibited lower CsLHCB gene expression. CsLHCB2 may serve as an indicator of temperature sensitivity differences in albino tea cultivars. This study could provide a reference for further studies of the functions of the CsLHCB family and contribute to research on the mechanism of the albino in tea plant.

Comparative genomics of light harvesting chlorophyll (LHC) gene family and impact of chlorophyll-A contents under drought stress in Helianthus annuus

Drought is one of the main environmental stressors that can alter the water status of plants; negatively affect growth, assimilation, and photosynthesis; and eventually reduce crop yield. We explored the dependence of drought tolerance traits on chlorophyll-A content. Local sunflower cultivars (FH-01, FH-628, FH-633, FH-572, and FH-653) were grown in pots and subjected to drought by withholding water for 10, 15, or 20 d. One month after germination, the leaves of the treated and non-treated plants were collected and subjected to biochemical analyses. Under different water stress levels, the levels of peroxidase (POD), superoxide dismutase (SOD), catalase (CAT), and proline increased, whereas those of chlorophyll-A decreased. Regression analysis clearly found that proline (-0.442), POD (-0.528), SOD (-0.532), and CAT (-0.814) have negative beta coefficient values. Phylogenetic analysis revealed that the LHC gene family is divided into six clades. Subcellular locations indicated that most LHC genes were located in the chloroplast; however, only few genes were present in the peroxisomes and endoplasmic reticulum. Our research found that Arabidopsis thaliana LHC genes were highly homologous to the LHC genes of Helianthus annuus. Furthermore, the LHC genes of both species are located in the chloroplasts; therefore, they play a role in photosynthesis and renewable energy production. This study opens a new horizon for discussing the role of chlorophyll-A in the drought-related traits of sunflowers.

Multi-level biological effects of diverse alkyl chains phthalate esters on cotton seedlings (Gossypium hirsutum L.): Insights into individual, physiological-biochemical and molecular perspectives

Phthalate esters (PAEs) are organic contaminants that pose environmental threat and safety risks to soil health and crop production. However, the ecological toxicity of different PAEs to cotton and the underlying mechanisms are not clear. This study investigated the ecotoxic effects and potential mechanisms of different alkyl-chain PAEs, including dioctyl phthalate (DOP), dibutyl phthalate (DBP), and diethyl phthalate (DEP) on cotton seedlings at multiple levels. The results showed that PAEs significantly hindered the growth and development of cotton. The chlorophyll content decreased by 1.87-31.66 %, accompanied by non-stomatal photosynthetic inhibition. The antioxidant system was activated by the three PAEs in cotton seedlings, while the osmotic potential was boosted intracellularly. Additionally, PAEs significantly interfered with functional gene expression and exhibited genotoxicity. Risk assessment results indicated that the ecotoxicity was DOP >DBP >DEP, with a "dose-response" relationship. The affinity between the three PAEs and catalase increased as the alkyl chain length increased, further supporting the toxicity sequence. Surprisingly, the bioconcentration factors of short-chain DEP were 8.07 ± 5.89 times and 1837.49 ± 826.83 times higher than those of long-chain DBP and DOP, respectively. These results support the ecological risk assessment of PAEs in cotton and provide new insights into determining the toxicity levels of different PAEs.

Comparative genomics analysis provide insights into evolution and stress responses of Lhcb genes in Rosaceae fruit crops

BACKGROUND

Light-harvesting chlorophyll a/b b evelopment of higher plants and in response to abiotic stress. Previous works has demonstrated that that Lhcb genes were involved in the phytochrome regulation and responded to the different light and temperature conditions in Poaceae (such as maize). However, the evolution and functions of Lhcb genes remains poorly characterized in important Rosaceae species.

RESULTS

In this investigation, we conducted a genome-wide analysis and identified a total of 212 Lhcb genes across nine Rosaceae species. Specifically, we found 23 Lhcb genes in Fragaria vesca, 20 in Prunus armeniaca, 33 in Malus domestica 'Gala', 21 in Prunus persica, 33 in Rosa chinensis, 29 in Pyrus bretschneideri, 18 in Rubus occidentalis, 20 in Prunus mume, and 15 in Prunus salicina. Phylogenetic analysis revealed that the Lhcb gene family could be classified into seven major subfamilies, with members of each subfamily sharing similar conserved motifs. And, the functions of each subfamily was predicted based on the previous reports from other species. The Lhcb proteins were highly conserved within their respective subfamilies, suggesting similar functions. Interestingly, we observed similar peaks in Ks values (0.1-0.2) for Lhcb genes in apple and pear, indicating a recent whole genome duplication event (about 30 to 45 million years ago). Additionally, a few Lhcb genes underwent tandem duplication and were located across all chromosomes of nine species of Rosaceae. Furthermore, the analysis of the cis-acting elements in the 2000 bp promoter region upstream of the pear Lhcb gene revealed four main categories: light response correlation, stress response correlation, hormone response correlation, and plant growth. Quantitative expression analysis demonstrated that Lhcb genes exhibited tissue-specific expression patterns and responded differently to low-temperature stress in Rosaceae species.

CONCLUSIONS

These findings shed light on the evolution and phylogeny of Lhcb genes in Rosaceae and highlight the critical role of Lhcb in pear's response to low temperatures. The results obtained provide valuable insights for further investigations into the functions of Lhcb genes in Rosaceae, and these functional genes will be used for further fruit tree breeding and improvement to cope with the current climate changes.

Identification of the Light-Harvesting Chlorophyll a/b Binding Protein Gene Family in Peach (Prunus persica L.) and Their Expression under Drought Stress

In higher plants, light-harvesting chlorophyll a/b binding (Lhc) proteins play a vital role in photosynthetic processes and are widely involved in the regulation of plant growth, development, and response to abiotic stress. However, the Lhc gene family has not been well identified in peaches (Prunus persica L.). In this study, 19 PpLhc genes were identified in the peach genome database, which were unevenly distributed on all chromosomes. Phylogenetic analysis demonstrated that PpLhc proteins could be divided into three major subfamilies, each of whose members had different exon-intron structures but shared similar conserved motifs. A total of 17 different kinds of cis-regulatory elements were identified in the promoter regions of all PpLhc genes, which could be classified into three categories: plant growth and development, stress response, and phytohormone response. In addition, transcriptomic data analysis and RT-qPCR results revealed that the expression profiles of some PpLhc genes changed under drought treatment, suggesting the crucial roles of Lhc genes in the regulation of plant tolerance to drought stress. Taken together, these findings will provide valuable information for future functional studies of PpLhc genes, especially in response to drought stress.

Genome-wide analysis and identification of light-harvesting chlorophyll a/b binding (LHC) gene family and BSMV-VIGS silencing TaLHC86 reduced salt tolerance in wheat

The discovery and identification of gene families by using wide-genome and public databases is an effective way to gain initial insight into gene function, which also is one of the current hot spots of research. Chlorophyll ab-binding proteins (LHC) are important for photosynthesis and widely involved in plant adversity stress. However, the study in wheat has not been reported. In this study, we identified 127 TaLHC members from common wheat which were unevenly distributed on all chromosomes except 3B and 3D. All members divided into three subfamilies, LHC a, LHC b and the LHC t which was only discovered in wheat. All of them had maximum expression in leaves and contained multiple light-responsive cis-acting element, which were evidence of the extensive involvement of LHC families in photosynthesis. In addition, we also analyzed their collinear relationship, targeting relationship with miRNA and their responses under different stresses. Based on these analyses, it was found that TaLHC86 was an excellent candidate gene for stress resistance. The full-length ORF of TaLHC86 was 792 bp and was localized on the chloroplasts. The salt tolerance of wheat was reduced when BSMV-VIGS silenced TaLHC86, and the photosynthetic rate and electron transport were also seriously affected. This study made a comprehensive analysis of the TaLHC family and found that TaLHC86 was a good gene for salt tolerance.

De novo full-length transcriptome analysis of two ecotypes of Phragmites australis (swamp reed and dune reed) provides new insights into the transcriptomic complexity of dune reed and its long-term adaptation to desert environments

BACKGROUND

The extremely harsh environment of the desert is changing dramatically every moment, and the rapid adaptive stress response in the short term requires enormous energy expenditure to mobilize widespread regulatory networks, which is all the more detrimental to the survival of the desert plants themselves. The dune reed, which has adapted to desert environments with complex and variable ecological factors, is an ideal type of plant for studying the molecular mechanisms by which Gramineae plants respond to combinatorial stress of the desert in their natural state. But so far, the data on the genetic resources of reeds is still scarce, therefore most of their research has focused on ecological and physiological studies.

RESULTS

In this study, we obtained the first De novo non-redundant Full-Length Non-Chimeric (FLNC) transcriptome databases for swamp reeds (SR), dune reeds (DR) and the All of Phragmites australis (merged of iso-seq data from SR and DR), using PacBio Iso-Seq technology and combining tools such as Iso-Seq3 and Cogent. We then identified and described long non-coding RNAs (LncRNA), transcription factor (TF) and alternative splicing (AS) events in reeds based on a transcriptome database. Meanwhile, we have identified and developed for the first time a large number of candidates expressed sequence tag-SSR (EST-SSRs) markers in reeds based on UniTransModels. In addition, through differential gene expression analysis of wild-type and homogenous cultures, we found a large number of transcription factors that may be associated with desert stress tolerance in the dune reed, and revealed that members of the Lhc family have an important role in the long-term adaptation of dune reeds to desert environments.

CONCLUSIONS

Our results provide a positive and usable genetic resource for Phragmites australis with a widespread adaptability and resistance, and provide a genetic database for subsequent reeds genome annotation and functional genomic studies.

Analysis of Lhc family genes reveals development regulation and diurnal fluctuation expression patterns in Cyperus esculentus, a Cyperaceae plant

Sixteen Lhc genes representing 13 phylogenetic groups were identified from the full-length transcriptome of tigernut, exhibiting development regulation and diurnal fluctuation expression patterns in leaves. Nuclear encoded light-harvesting chlorophyll a/b-binding (Lhc) proteins play indispensable roles in oxygenic photosynthesis. In this study, we present the first transcriptome-based characterization of Lhc family genes in tigernut (Cyperus esculentus L.), a Cyperaceae C₄ plant producing oil in underground tubers. A number of 16 Lhc genes representing 13 phylogenetic groups identified from the full-length tigernut transcriptome are equal to that found in both Carex littledalei (another Cyperaceae plant) and papaya, slightly more than 15 members present in both rice and jatropha, but relatively less than 18, 20, and 21 members present in sorghum, cassava, and Arabidopsis, respectively. Nevertheless, nearly one-vs-one orthologous relationship was observed in most groups, though some of them are no longer located in syntenic blocks and species-specific expansion was frequently found in Lhcb1. Comparative genomics analysis revealed that the loss of two groups (i.e., Lhca2 and Lhca5) in C. littledalei is species-specific, sometime after the split with tigernut, and the expansion of Lhcb1 was mainly contributed by tandem duplication as observed in most species. Interestingly, a transposed duplication, which appears to be shared by monocots, was also identified in Lhcb1. Further transcriptome profiling revealed a predominant expression pattern of most CeLhc family genes in photosynthetic tissues and enhanced transcription during leaf maturation, reflecting their key roles in light absorption. Moreover, qRT-PCR analysis revealed an apparent diurnal fluctuation expression pattern of 11 dominant CeLhc genes. These findings not only highlight species-specific evolution of Lhc genes in the Cyperaceae family as well as the monocot lineage, but also provide valuable information for further functional analysis and genetic improvement in tigernut.

Phylogenetic, Structural and Functional Evolution of the LHC Gene Family in Plant Species

Light-harvesting chlorophyll a/b-binding (LHC) superfamily proteins play a vital role in photosynthesis. Although the physiological and biochemical functions of LHC genes have been well-characterized, the structural evolution and functional differentiation of the products need to be further studied. In this paper, we report the genome-wide identification and phylogenetic analysis of LHC genes in photosynthetic organisms. A total of 1222 non-redundant members of the LHC family were identified from 42 species. According to the phylogenetic clustering of their homologues with Arabidopsis thaliana, they can be divided into four subfamilies. In the subsequent evolution of land plants, a whole-genome replication (WGD) event was the driving force for the evolution and expansion of the LHC superfamily, with its copy numbers rapidly increasing in angiosperms. The selection pressure of photosystem II sub-unit S (PsbS) and ferrochelatase (FCII) families were higher than other subfamilies. In addition, the transcriptional expression profiles of LHC gene family members in different tissues and their expression patterns under exogenous abiotic stress conditions significantly differed, and the LHC genes are highly expressed in mature leaves, which is consistent with the conclusion that LHC is mainly involved in the capture and transmission of light energy in photosynthesis. According to the expression pattern and copy number of LHC genes in land plants, we propose different evolutionary trajectories in this gene family. This study provides a basis for understanding the molecular evolutionary characteristics and evolution patterns of plant LHCs.

Genome-Wide Identification of the LHC Gene Family in Kiwifruit and Regulatory Role of AcLhcb3.1/3.2 for Chlorophyll a Content

Light-harvesting chlorophyll a/b-binding (LHC) protein is a superfamily that plays a vital role in photosynthesis. However, the reported knowledge of LHCs in kiwifruit is inadequate and poorly understood. In this study, we identified 42 and 45 LHC genes in Actinidia chinensis (Ac) and A. eriantha (Ae) genomes. Phylogenetic analysis showed that the kiwifruit LHCs of both species were grouped into four subfamilies (Lhc, Lil, PsbS, and FCII). Expression profiles and qRT-PCR results revealed expression levels of LHC genes closely related to the light, temperature fluctuations, color changes during fruit ripening, and kiwifruit responses to Pseudomonas syringae pv. actinidiae (Psa). Subcellular localization analysis showed that AcLhcb1.5/3.1/3.2 were localized in the chloroplast while transient overexpression of AcLhcb3.1/3.2 in tobacco leaves confirmed a significantly increased content of chlorophyll a. Our findings provide evidence of the characters and evolution patterns of kiwifruit LHCs genes in kiwifruit and verify the AcLhcb3.1/3.2 genes controlling the chlorophyll a content.

Identification and expression analysis of chlorophyll a/b binding protein gene family in grape (Vitis vinifera)

In higher plants, light capture of chlorophyll a/b binding protein (Lhc) plays a crucial role in the plant's response to adverse environment. So far, the family has not been systematically identified in grapes. In this study, 20 VvLhcs were identified in the grape genome, which were distributed in 13 of 19 grape chromosomes and divided into 7 developing branches. The results of gene duplication analysis showed that 6 VvLhcs formed fragment duplication events, while there was no tandem duplication in VvLhcs. Exon-intron structure analysis showed that they had a wide number of exons. Protein conserved motif analysis showed that VvLhcs contained more similar motif structures in the same phylogenetic branch. The cis-acting elements in the VvLhcs promoter region mainly respond to light, plant hormones and abiotic stresses. In addition, qRT-PCR results showed that different proportions of salt stress and red-blue light affected the expression of VvLhcs and the expression patterns of genes in different grape varieties were different. The results for further study on different grape varieties in different combinations of red and blue light of the Lhc provide a theoretical basis.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-022-01204-5.

Identification of Differentially Expressed Genes Reveal Conserved Mechanisms in the Rice-Magnaporthe oryzae Interaction

Magnaporthe oryzae causes rice blast disease and is responsible for major losses in rice production worldwide. Although numerous studies have focused on the interactions between Oryza sativa and M. oryzae, to date, the conserved mechanisms remain in part unclear. In this study, a comparative analysis of transcriptomes of O. sativa L. ssp. japonica cv. 'Nipponbare' interacting with three M. oryzae strains (248, 235, and 163) were performed to explore the conserved molecular mechanisms. Differentially expressed genes with similar expression patterns in the interactions between cultivar 'Nipponbare' and three M. oryzae strains were defined as Conserved Differentially Expressed Genes (CDEGs). These included 3,647 O. sativa CDEGs and 3,655 M. oryzae CDEGs. Four rice CDEGs (LOC_Os03g19270, LOC_Os07g36600, LOC_Os05g28740, and LOC_Os01g32780) encoding universal stress protein (USP) were induced within 24 h post-inoculation (hpi) by three M. oryzae strains. Meanwhile, overexpression of LOC_Os07g36600 resulted in enhanced rice resistance against M. oryzae. Furthermore, four rice genes coding light-harvesting chlorophyll a/b-binding (LHC) protein (LOC_Os02g52650, LOC_Os09g12540, LOC_Os11g13850, LOC_Os05g22730) were also identified as CDEGs and were induced at 48 hpi, which might contribute to blast resistance through reactive oxygen species (ROS) accumulation. MoCDIP4 is M. oryzae effector inducing rice cell death and were verified that include AA9 CAZy domain (namely GH61 domain). In this study, we found seven MoCDIP4-homologous genes coding proteins with signal peptides and AA9 CAZy domains, which were continuously up-regulated across all infection stages relative to uninoculated control. This study uncovered that genes are required for conserved mechanisms of rice-M. oryzae interaction, which includes rice genes encoding USP proteins and LHC proteins, as well as M. oryzae genes encoding AA9 proteins. This study will help us to understand how O. sativa responds to M. oryzae infections and the molecular mechanisms of M. oryzae pathogenicity.

Characterization of cassava ORANGE proteins and their capability to increase provitamin A carotenoids accumulation

Cassava (Manihot esculenta Crantz) biofortification with provitamin A carotenoids is an ongoing process that aims to alleviate vitamin A deficiency. The moderate content of provitamin A carotenoids achieved so far limits the contribution to providing adequate dietary vitamin A levels. Strategies to increase carotenoid content focused on genes from the carotenoids biosynthesis pathway. In recent years, special emphasis was given to ORANGE protein (OR), which promotes the accumulation of carotenoids and their stability in several plants. The aim of this work was to identify, characterize and investigate the role of OR in the biosynthesis and stabilization of carotenoids in cassava and its relationship with phytoene synthase (PSY), the rate-limiting enzyme of the carotenoids biosynthesis pathway. Gene and protein characterization of OR, expression levels, protein amounts and carotenoids levels were evaluated in roots of one white (60444) and two yellow cassava cultivars (GM5309-57 and GM3736-37). Four OR variants were found in yellow cassava roots. Although comparable expression was found for three variants, significantly higher OR protein amounts were observed in the yellow varieties. In contrast, cassava PSY1 expression was significantly higher in the yellow cultivars, but PSY protein amount did not vary. Furthermore, we evaluated whether expression of one of the variants, MeOR_X1, affected carotenoid accumulation in cassava Friable Embryogenic Callus (FEC). Overexpression of maize PSY1 alone resulted in carotenoids accumulation and induced crystal formation. Co-expression with MeOR_X1 led to greatly increase of carotenoids although PSY1 expression was high in the co-expressed FEC. Our data suggest that posttranslational mechanisms controlling OR and PSY protein stability contribute to higher carotenoid levels in yellow cassava. Moreover, we showed that cassava FEC can be used to study the efficiency of single and combinatorial gene expression in increasing the carotenoid content prior to its application for the generation of biofortified cassava with enhanced carotenoids levels.

Genome-wide identification of the light-harvesting chlorophyll a/b binding (Lhc) family in Gossypium hirsutum reveals the influence of GhLhcb2.3 on chlorophyll a synthesis

Light-harvesting chlorophyll a/b binding (Lhc) family proteins play a significant role in photosynthetic processes. Our objective was systematic identification and analysis of the Lhc family in cotton, as well as the relationship between Lhc family genes and chlorophyll synthesis during photosynthetic processes. We used genome-wide identification, phylogenetic analysis, chromosomal distribution and collinearity to examine potential functions of Lhc superfamily genes in upland cotton. Subcellular localization, qRT-PCR, a yeast two hybrid (Y2H) , and Virus-induced gene silencing (VIGS) experiment were used to explore function of GhLhcb2.3. Focusing on GhLhc family, gene structural analysis of G. hirsutum Lhc family genes (GhLhc) indicated the conservation of selected Lhc family members. The expression pattern of GhLhc proteins shows that Lhc family proteins are important for photosynthetic processes in leaves. Results of subcellular localization and qRT-PCR in different cotton varieties showed that GhLhcb2.3 is closely related to chloroplast chlorophyll. Y2H found extensive heteromeric interactions between the GhLhcb2.3 and GhLhcb1.4. Subcellular localization revealed that GhLhcb1.4 is located in chloroplasts. VIGS showed that GhLhcb2.3 influenced chlorophyll a synthesis. We comprehensively identified Lhc family genes in cotton, characterized these genes and reveal the influence of GhLhcb2.3 on chlorophyll a synthesis.

Comparative physiological and transcriptomic analysis of pear leaves under distinct training systems

Canopy architecture is critical in determining the light interception and distribution, and subsequently the photosynthetic efficiency and productivity. However, the physiological responses and molecular mechanisms by which pear canopy architectural traits impact on photosynthesis remain poorly understood. Here, physiological investigations coupled with comparative transcriptomic analyses were performed in pear leaves under distinct training systems. Compared with traditional freestanding system, flat-type trellis system (DP) showed higher net photosynthetic rate (P_N) levels at the most time points throughout the entire monitored period, especially for the interior of the canopy in sunny side. Gene ontology analysis revealed that photosynthesis, carbohydrate derivative catabolic process and fatty acid metabolic process were over-represented in leaves of DP system with open-canopy characteristics. Weighted gene co-expression network analysis uncovered a significant network module positive correlated with P_N value. The hub genes (PpFKF1 and PpPRR5) of the module were enriched in circadian rhythm pathway, suggesting a functional role for circadian clock genes in mediating photosynthetic performance under distinct training systems. These results draw a link between pear photosynthetic response and specific canopy architectural traits, and highlight light harvesting and circadian clock network as potential targets for the input signals from the fluctuating light availability under distinct training systems.

Genome-wide analysis of the light-harvesting chlorophyll a/b-binding gene family in apple (Malus domestica) and functional characterization of MdLhcb4.3, which confers tolerance to drought and osmotic stress

In higher plants, the light-harvesting chlorophyll a/b-binding (Lhc) proteins function in multiple processes that are critical to plant growth, development, and abiotic stress response. However, the Lhc gene family has not been well characterized in the important fruit crop, apple (Malus × domestica Borkh.). In this study, we identified 27 Lhc genes in the apple genome. Phylogenetic analysis showed that the Lhc gene family could be classified into three major subfamilies, each of whose members shared similar conserved motifs. Evolutionary analysis indicated that duplicated MdLhc genes were primarily under purifying selection. MdLhcs were expressed at varying levels in all tissues examined and showed different expression patterns under drought stress. The overexpression of MdLhcb4.3 in transgenic Arabidopsis and apple callus enhanced their tolerance to drought and osmotic stress. Taken together, these results demonstrate the important role of Lhc proteins in the regulation of plant resistance to drought and osmotic stress and provide valuable information for further study of Lhc functions in apple.

Genes encoding light-harvesting chlorophyll a/b-binding proteins in papaya (Carica papaya L.) and insight into lineage-specific evolution in Brassicaceae

Light-harvesting chlorophyll a/b-binding (Lhc) proteins comprise a plant-specific superfamily involved in photosynthesis and stress responses. Despite their importance, little is known in papaya (Carica papaya), an economically important tree fruit crop as well as a species close to the model plant arabidopsis (Arabidopsis thaliana). This study reports a first genome-wide analysis of Lhc superfamily genes in papaya, and a total of 28 members that represent four defined families or 26 orthologous groups were identified from the papaya genome. The superfamily number is comparable to 28 or 27 reported in castor (Ricinus communis) and jatropha (Jatropha curcas), respectively, two Euphorbiaceous plants also without any recent whole-genome duplication (WGD), but relatively less than 35, 34, 32, 32, 37, 30 or 32 present in cassava (Manihot esculenta), arabidopsis, A. lyrata, A. halleri, Capsella rubella, C. grandiflora, and Eutrema salsugineum, respectively, representative species having experienced one or two recent WGDs. Local duplication was shown to play a predominant role in gene expansion in papaya, castor, and jatropha, which is only confined to the Lhcb1 group. By contrast, WGD plays a relatively more important role in cassava, arabidopsis, and other Brassicaceous plants. Further comparison of Brassicaceous plants revealed that loss of the SEP6 group in arabidopsis is lineage-specific, occurring sometime after papaya-arabidopsis divergence but before the radiation of Brassicaceous plants. Transcriptional profiling revealed a leaf-preferential expression pattern of most CpLhc superfamily genes and their transcript levels were markedly regulated by three abiotic stresses, i.e., mimicking drought, cold, and high salt. These findings not only facilitate further functional studies in papaya, but also improve our knowledge on lineage-specific evolution of this special gene superfamily in Brassicaceae.

Light-harvesting chlorophyll a/b-binding protein-coding genes in jatropha and the comparison with castor, cassava and arabidopsis

The Lhc (light-harvesting chlorophyll a/b-binding protein) superfamily represents a class of antennae proteins that play indispensable roles in capture of solar energy as well as photoprotection under stress conditions. Despite their importance, little information has been available beyond model plants. In this study, we presents a first genome-wide analysis of Lhc superfamily genes in jatropha (Jatropha curcas L., Euphorbiaceae), an oil-bearing plant for biodiesel purpose. A total of 27 members were identified from the jatropha genome, which were shown to distribute over nine out of the 11 chromosomes. The superfamily number is comparable to 28 present in castor (Ricinus communis, Euphorbiaceae), but relatively less than 35 in cassava (Manihot esculenta, Euphorbiaceae) and 34 in arabidopsis (Arabidopsis thaliana) that experienced one or two recent whole-genome duplications (WGDs), respectively. In contrast to a high number of paralogs present in cassava and arabidopsis, few duplicates were found in jatropha as observed in castor, corresponding to no recent WGD occurred in these two species. Nevertheless, 26 orthologous groups representing four defined families were found in jatropha, and nearly one-to-one orthologous relationship was observed between jatropha and castor. By contrast, a novel group named SEP6 was shown to have been lost in arabidopsis. Global transcriptome profiling revealed a predominant expression pattern of most JcLhc superfamily genes in green tissues, reflecting their key roles in photosynthesis. Moreover, their expression profiles upon hormones, drought, and salt stresses were also investigated. These findings not only improve our knowledge on species-specific evolution of the Lhc supergene family, but also provide valuable information for further studies in jatropha.

Genomics analysis of genes encoding respiratory burst oxidase homologs (RBOHs) in jatropha and the comparison with castor bean

Respiratory burst oxidase homologs (RBOHs), which catalyze the production of superoxide from oxygen and NADPH, play key roles in plant growth and development, hormone signaling, and stress responses. Compared with extensive studies in model plants arabidopsis and rice, little is known about RBOHs in other species. This study presents a genome-wide analysis of Rboh family genes in jatropha (Jatropha curcas) as well as the comparison with castor bean (Ricinus communis), another economically important non-food oilseed crop of the Euphorbiaceae family. The family number of seven members identified from the jatropha genome is equal to that present in castor bean, and further phylogenetic analysis assigned these genes into seven groups named RBOHD, -C, -B, -E, -F, -N, and -H. In contrast to a high number of paralogs present in arabidopsis and rice that experienced several rounds of recent whole-genome duplications, no duplicate was identified in both jatropha and castor bean. Conserved synteny and one-to-one orthologous relationship were observed between jatropha and castor bean Rboh genes. Although exon-intron structures are usually highly conserved between orthologs, loss of certain introns was observed for JcRbohB, JcRbohD, and RcRbohN, supporting their divergence. Global gene expression profiling revealed diverse patterns of JcRbohs over various tissues. Moreover, expression patterns of JcRbohs during flower development as well as various stresses were also investigated. These findings will not only improve our knowledge on species-specific evolution of the Rboh gene family, but also provide valuable information for further functional analysis of Rboh genes in jatropha.

Genome-wide comparison reveals divergence of cassava and rubber aquaporin family genes after the recent whole-genome duplication

Background

Aquaporins (AQPs) are a class of integral membrane proteins that facilitate the passive transport of water and other small solutes across biological membranes. Despite their importance, little information is available in cassava (Manihot esculenta), a perennial shrub of the Euphorbiaceae family that serves the sixth major staple crop in the world.

Results

This study presents a genome-wide analysis of the AQP gene family in cassava. The family of 42 members in this species could be divided into five subfamilies based on phylogenetic analysis, i.e., 14 plasma membrane intrinsic proteins (PIPs), 13 tonoplast intrinsic proteins (TIPs), nine NOD26-like intrinsic proteins (NIPs), four X intrinsic proteins (XIPs), and two small basic intrinsic proteins (SIPs). Best-reciprocal-hit-based sequence comparison and synteny analysis revealed 34 orthologous groups (OGs) present in the Euphorbiaceae ancestor, and nearly one-to-one or two-to-one orthologous relationships were observed between cassava with rubber/physic nut, reflecting the occurrence of one so-called ρ recent whole-genome duplication (WGD) in the last common ancestor of cassava and rubber. In contrast to a predominant role of the ρ WGD on family expansion in rubber, cassava AQP duplicates were derived from the WGD as well as local duplication. Species-specific gene loss was also observed in cassava, which includes the entire NIP4 group and/or six OGs. Comparison of conserved motifs and gene expression profiles revealed divergence of paralogs in cassava as observed in rubber.

Conclusions

Our findings will not only improve our knowledge on family evolution in Euphorbiaceae, but also provide valuable information for further functional analysis of AQP genes in cassava and rubber.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5780-4) contains supplementary material, which is available to authorized users.