Research on mulberry photosynthesis transcriptome

doi:10.14088/j.cnki.issn0439-8114.2019.23.023

Abstract

Abstract: The photosynthesis characteristics and photosynthetic related gene expression of high-yield mulberry varieties Esang 1 （E1） and control mulberry variety Husang 32（H32） were studied to provide a theoretical basis for high-efficiency breeding of mulberry trees and screening of germplasm resources. Using LI-6400XT mulberry leaf gas exchange and chlorophyll fluorescence parameters, chlorophyll content and ribulose carboxylase（RUBP） activity were measured simultaneously, transcriptome sequencing was performed on the peaks and valleys of mulberry photosynthesis curves, and the differences in photosynthetic related gene expression among different mulberry varieties were compared. The results of photosynthetic variation of E1 and H32 showed that the Pn peak of E1 leaves was significantly larger than that of H32（P<0.05）, and the peak of Tr of H32 leaves was significantly larger than that of E1 （P<0.05）. In the Pn-PAR and Pn-Ci responses of different mulberry leaves, CE and CSP were E1>H32, AQY, LCP, LSP and CCP were H32>E1. The ability of E1 to achieve photosynthetic product accumulation under weaker light intensity is stronger than that of H32, and the utilization efficiency of E1 to low concentration CO₂ is better than that of H32; among the chlorophyll fluorescence parameters of different mulberry varieties, ΦPSII and qP of E1 leaves were significantly higher than H32 （P<0.05）, and NPQ of H32 was significantly larger than E1 （P<0.05）, indicating that the light energy absorbed by H32 leaves was mainly dissipated in the form of heat energy; the chlorophyll a, chlorophyll b and chlorophyll content of H32 were significantly lower than E1（P<0.05）, and the RUBP activity of E1 was higher than that of H32（P>0.05）, indicating that the photosynthetic capacity of E1 was better than that of H32. 3 356 DEGs were found in different varieties and different time periods. The metabolic pathway analysis showed that there were 1 136 DEGs annotated to the KEEG database. The photosynthesis-related metabolic pathways significantly enriched by DEGs were mainly photosynthetic carbon sequestration, carbon metabolism, porphyrin and chlorophyll metabolism. Photosynthetically related metabolic pathways. A total of 10 new genes related to photosynthesis were screened for highly differentially expressed genes with FC≥2 and FDR<0.01. The four DEGs, Mulberry_newGene_3646, Mulberry_newGene_1405, Mulberry_newGene_2419 and Mulberry_newGene_320, may be the key genes causing the difference in net photosynthetic rate and yield between H32 and E1, which can be studied intensively.

Key words: mulberry, photosynthetic characteristics, chlorophyll fluorescence characteristics, transcriptome

CLC Number:

S888

LI Yong, YU Cui, MO Rong-li, DENG Wen, XIONG Chao, ZHUANG Chu-xiong, HU Xing-ming. Research on mulberry photosynthesis transcriptome[J]. HUBEI AGRICULTURAL SCIENCES, 2019, 58(23): 101-111.

References 29

[1]	许大全. 光合作用效率[M].上海:上海科学技术出版社,2002.
[2]	CHEN F G,CHEN L G,ZHAO H X,et al.Sex-specific responses and tolerances of Populus cathayana to salinity[J].Physiol plant 2010,140:163-173.
[3]	PENUELAS J,FILELLA I,LLUSIA J,et al.Comparative field study of spring and summer leaf gas exchange and photobiology of the Mediterranean trees Quercusilex and Phillyrea latifolia[J].Journal of experiment botany,1998,49(319):229-238.
[4]	BERRY J A,DOWNTON W J S. Environmental regulation of photosynthesis[A].Govindjee photosynthesis[M].New York:Academic Press,1982.
[5]	LONG S P,BAKER N R,RAINS C A.Analyzing the responses of photosynthetic CO2 assimilation to long-term elevation of atmospheric CO2 concentration[J].Vegetatio,1993,104:33-45.
[6]	GENTY B,BRIANTIS J M,BAKER N R.The relationship between the quantum yield of photosynthetic electron transport and quenching of chlorophyll fluorescence[J].Biochim Biophys Acta,1989,894:183-192.
[7]	DEMMIG B,WINTER K,KRUGER A,et al.Photoinhibition and zeaxanthin formation in intact leaves:A possible role of the xanthophyllcycle in the dissipation of excess light[J].Plant physiology,1987,84:218-224.
[8]	PORRA RJ,THOMPSON WA,KRIEDEMANN PE.Determination of accurate extinction coefficients and simultaneous equations for assaying chlorophylls a and b extracted with four different solvents:verification of the concentration of chlorophyll standards by atomic absorption spectroscopy[J].Biochimica et Biophysica Acta Bioenergetics,1989,975:384-394.
[9]	BROUWER B,ZIOLKOWSKA A,BAGARD M, KEECH O,et al.The impact of light intensity on shade-induced leaf senescence[J].Plant cell and environment,2012,35:1084-1098.
[10]	ALTSCHUL SF,MADDEN TL,ZHANG J,et al.Gapped BLAST and PIBLASTS:A new gene ration of protein database search programs[J].Nucleic acids rsearche ,1997,25(17):3389-3402
[11]	DENG Y Y,LI J Q,WU S F,etal.Integrated nr database in protein annotation system and its localization[J].Computer engineering,2006,32(5):71-74.
[12]	APWEILER R,BAIROCH A,WU CH,etal.The universal protein knowledge base[J].Nucleic aids rsearchce,2004,32:D115-D119.
[13]	ASHBURNER M,BALL C A,BLAKE JA,et al.Gene ontology:tool for the unification of biology[J].Nature genetics,2000,25(1):25-29.
[14]	TATUSOV R L,GALPERIN M Y,NATALE D A.The COG database:A tool for genome scale analysis of protein functions and evolution[J].Nucleic aids research,2000,28(1):33-36.
[15]	KOONIN E V,FEDOROVA N D,JACKSON JD,et al.A comprehensive evolutionary classification of proteins encoded in complete eukaryotic genomes[J].Genomebiology,2004,5(2):R7.
[16]	FINN R D,BATEMAN A,CLEMENTS J,et al.Pfam:The protein families database[J].Nucleic aids research,2014,42:D222.
[17]	KANEHISA M,GOTO S,KAWASHIMa S,et al.The KEGG resource for deciphering the genome[J].Nucleic acids research,2004,32:D277-D280.
[18]	XIE C,MAO X,HUANG J,et al.KOBAS 2.0:A web server for annotation and identification of enriched pathways and diseases[J].Nucleic aids research,2011,39:W316-W322.
[19]	EDDY S R.Profile hidden Markov models[J].Bioinformatics,1998,14(9):755-763.
[20]	LiVAK KJ,SCHMITTGEN TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) method[J].Methods,2001,25:402-408.
[21]	CALNEL M G R,THOMLEY J H M.Temperature and CO2 responses of leaf and canopy photosynthesis:A clarification using the non-rectangularHyperbola model of photosynthesis.[J].Annals of botany,1998,82:883-892.
[22]	MORALES M,ABADÍA A,ABADÍA J.Photoinhibition and photoprotection under nutrient deficiencies,drought and salinity[A].DEMMIG A B,ADAMS W W,MATTOO A K. Photoprotection,gene regulation,and environment[M].Berlin:Springer-Verlag,2008.
[23]	LU C M,ZHANG J H.Heat-induced multiple effects on PS II in wheat plants[J].J Plant Physiol,2000,156:259-265.
[24]	EFEOLU B,EKMEKÇI Y,ÇIÇEK N. Physiological responses of three maize cultivars to drought stress and recovery[J].J S Afr Bot,2009,75:34-42.
[25]	SCHRADER S M,WISE R R,WACHOLTZ W F,et,al.Thylakoid membrane responses to moderately high leaf temperature in Pima cotton[J].Plant Cell Environ,2004,27:725-735.
[26]	SHARKEY T D.Effects of moderate heat stress on photosynthesis:Importance of thylakoid reactions, rubisco deactivation, reactive oxygen species,and thermotolerance provided by isoprene[J].Plant Cell Environ,2005,28:269-277.
[27]	FERREIRA S,HJERNØ K,LARSEN M,et al.Proteome profiling of Populus euphratica Oliv. upon heat stress[J].Ann Bot,2006, 98:361-77.
[28]	MASCLAUX-DAUBRESSE C,REISDORF-CREN M,ORSEL M.Leaf nitrogen remobilisation for plant development and grain filling[J].Plant biology,2008,10,23-36.
[29]	Parry M A.Rubisco activity and regulation as targets for crop improvement[J].J Exp Bot,2013,64(3):717-730.