Introduction
In the face of potentially stressful challenges, activations of neural, neuroendocrine and immune systems promote adaption (termed ‘allostasis’) to maintain homeostasis. However, these allostatic systems can cause maladaption or disease if they are overactive or underactive over long periods. It thus appears that the stress response has both positive and negative effects. Despite the obvious positive effects of the stress response for homeostasis and health, the negative effects for the body have recently received increasing attention.
For laying hens reared under commercial conditions, stress is a common problem disrupting their physiological homeostasis and laying performance (Shini, Reference SHINI2003; Pohle and Cheng, Reference POHLE and CHENG2009). Such birds are exposed to various stressors such as repeated restraint (i.e., weighing, vaccination, moving to the layer shed, etc.) and a variety of environmental changes. These challenges typically trigger a cascade of endocrine secretions involving the hypothalamic-pituitary-adrenal (HPA) axis, which regulates secretion of glucocorticoids to mediate the necessary allostatic adjustments (Blas, Reference BLAS2015). Corticosterone is the predominant avian glucocorticoid and a hallmark of the adrenocortical response to stress in birds (Blas, Reference BLAS2015). Glucocorticoids transmit information about environmental conditions to the hypothalamic-pituitary-gonadal (HPG) axis, which ultimately influences egg production. Energy intake and distribution between organs or tissues are influenced.
Glucocorticoids and energy intake
evoked HPA axis responses: a consequence of energy deficiencies
A stressful environment causes an energetic challenge. The critical role of glucocorticoids in the regulation of energy balance in response to the energetic challenge caused by a stressor is well documented. Plasma levels of corticosterone, the primary glucocorticoid in birds (Holmes and Phillips, Reference HOLMES and PHILLIPS1976), are increased during periods of reduced food supply associated with cold and poor weather (Wingfield, Reference WINGFIELD1985; Astheimer et al., Reference ASTHEIMER, BUTTEMER and WINGFIELD1995; Jenni-Eiermann et al., Reference JENNI-EIERMANN, GLAUS, GRUEBLER, SCHWABL and JENNI2008) and during periods of experimental acute fasting (Lynn et al., Reference LYNN, BREUNER and WINGFIELD2003). Although corticosterone responses to fasting can vary among individuals and depend on energy reserves and experience (Richards, Reference RICHARDS2003), the majority of studies demonstrated an elevated corticosterone concentration during fasting challenge in chickens. An increase in plasma corticosterone is observed after food deprivation in both immature (Geris et al., Reference GERIS, BERGHMAN, KUHN and DARRAS1999) and adult chickens (Scanes et al., Reference SCANES, MERRILL, FORD, MAUSER and HOROWITZ1980). This means that food shortage acts as a potent activator of the HPA axis in birds, supporting the suggestion that a link exists between food availability and HPA axis activation in birds (Lynn et al., Reference LYNN, BREUNER and WINGFIELD2003). Given the critical role of glucocorticoids in the regulation of energy balance, it is possible that glucocorticoid secretion in response to energy constraints such as fasting act through a series of mechanisms, thereby suppressing egg-laying when conditions are not appropriate to support breeding (Perfito et al., Reference PERFITO, ZANN, BENTLEY and HAU2007).
Glucocorticoids and appetite
A number of reports document the effect of glucocorticoids on the food intake of chickens. Corticosterone induces a dose-dependent increase in food intake in poultry (Covasa and Forbes, Reference COVASA and FORBES1995; El-Lethey et al., Reference EL-LETHEY, JUNGI and HUBER-EICHER2001; Liu et al., Reference LIU, SONG, JIAO and LIN2014) and is involved in the control of appetite (Yuan et al., Reference YUAN, LIN, JIANG, JIAO and SONG2008). Moreover, the effect of corticosterone on food intake in birds depends on their diets, that is, corticosterone-challenge, causes a high-energy diet preference, such as high-protein diets (Bartov, Reference BARTOV1985) and high-fat diets (Liu et al., Reference LIU, SONG, KONG, JIAO and LIN2012a). Several features of high-fat foods, such as palatable flavour and delayed satiating effects, may contribute to the overeating response (Blundell and MacDiarmid, Reference BLUNDELL and MACDIARMID1997).
Many neuropeptides, so-called appetite regulatory peptides, are expressed in the hypothalamus and gastrointestinal tract and are involved in the regulation of food intake in birds (Richards and Proszkowiec-Weglarz, Reference RICHARDS and PROSZKOWIEC-WEGLARZ2007). A dynamic balance between appetite-associated peptides in the hypothalamus and gastrointestinal tract controls the feeding status of chickens under glucocorticoid treatment (Liu et al., Reference LIU, SONG, SHEIKHAHMADI, JIAO and LIN2012b). Some gastrointestinal tract peptides, such as ghrelin and cholecystokinin (CCK), may play pivotal roles in the reduction of food intake in both immature chickens (Liu et al., Reference LIU, LI, LI, JIAO, LIN, SHEIKHAHMADI, WANG and SONG2013) and laying hens (Song et al., Reference SONG, LIU, SHEIKHAHMADI, JIAO and LIN2012a) exposed to heat stress, whereas some hypothalamic neurons, such as neuropeptide Y (NPY), are associated with the regulation of appetite by glucocorticoids (Liu et al., Reference LIU, SONG, JIAO and LIN2014; Reference LIU, WANG, JIAO, ZHAO and LIN2015).
The hypothalamus is important in interpreting external environmental cues such as stressors and generates the appropriate responses for food intake in birds (Richards et al., Reference RICHARDS, ROSEBROUGH, COON and MCMURTRY2010). It is well-recognised that hypothalamic AMP-activated protein kinase (AMPK) signalling is involved in the control of food intake (Minokoshi et al., Reference MINOKOSHI, ALQUIER, FURUKAWA, KIM, LEE, XUE, MU, FOUFELLE, FERRE, BIRNBAUM, STUCK and KAHN2004). AMPK acts as a sensor of cellular energy changes and is activated by increasing AMP levels coupled with decreasing ATP levels. After activation, AMPK phosphorylates acetyl-CoA carboxylase (ACC) and initiates energy-producing pathways at the expense of energy-depleting processes (Kahn et al., Reference KAHN, ALQUIER, CARLING and HARDIE2005). Fasting evokes hypothalamic AMPK activation in chicks (Song et al., Reference SONG, LIU, YUE, JIAO, LIN, SHEIKHAHMADI, EVERAERT, DECUYPERE and BUYSE2012b). Activated AMPK stimulates NPY gene expression and food intake, whereas the inactivated AMPK induces anorexia in chicks (Kim et al., Reference KIM, PARK, NAMKOONG, JANG, RYU, SONG, YUN, NAMGOONG, HA, PARK, LEE, VIOLLET, YOUN, LEE and LEE2004; Liu et al., Reference LIU, SONG, JIAO and LIN2014). Moreover, a recent study indicated that corticosterone challenge causes a high-fat diet preference by enhancing NPY expression via the AMPK pathway in the hypothalamus. Liu et al. (Reference LIU, WANG, JIAO, ZHAO and LIN2015) suggested a potential role for the target of rapamycin (TOR) pathway of hypothalamus in the regulation of appetite-related genes by glucocorticoids, which must be further verified.
involvement of NPY in regulation of reproduction
Appetite-associated neuropeptides also integrate the reproductive and energy balance axes (Evans and Anderson, Reference EVANS and ANDERSON2012). Of these neuropeptides, NPY is considered to function as a switch in the brain that turns on appetite and turns off reproduction. The reciprocal control of appetite and reproduction is inherent in the NPY molecule itself (Kalra and Kalra, Reference KALRA and KALRA2004). NPY is an important neuromodulator affecting gonadal function in mammals and birds (Hilal et al., Reference HILAL, CHEN and SILVERMAN1996). In mammals, NPY expression causes a breeding failure (Kalra et al., Reference KALRA, CLARK, SAHU, DUBE and KALRA1988) that is manifested by inhibiting reproductive function (Kalra et al., Reference KALRA, XU, DUBE, MOLDAWER, MARTIN and KALRA1998) and suppressing sexual behaviour (Marin-Bivens et al., Reference MARIN-BIVENS, KALRA and OLSTER1998) in a receptor-dependent manner (Toufexis et al., Reference TOUFEXIS, KYRIAZIS and WOODSIDE2002; Gonzales et al., Reference GONZALES, VOIROL, GIACOMINI, GAILLARD, PEDRAZZINI and PRALONG2004; Klenke et al., Reference KLENKE, CONSTANTIN and WRAY2010). In chickens, the close association of NPY and gonadotropin-releasing hormone (GnRH), perikarya and fibres in the brain suggests a functional interaction between the two peptides in regulating gonadal function (Kuenzel, Reference KUENZEL2000). NPY may affect egg production rate through its role in the control of ovulation (Dunn et al., Reference DUNN, MIAO, MORRIS, ROMANOV, WILSON and WADDINGTON2004; Wu et al., Reference WU, LI, YAN, TANG, CHEN, WANG, GAO, TU, YU and ZHU2007; Li et al., Reference LI, ZHU, CHEN, WU, TANG, GAO, SONG and XU2009) and the secretion of hormones such as GnRH (Wang et al., Reference WANG, DAY and VASILATOS-YOUNKEN2001). NPY is the basis of a neuromodulatory link between nutrition and reproduction at the central nervous system level (Wójcik-Gładysz and Polkowska, 2006; Dhillon et al., Reference DHILLON, GINGERICH and BELSHAM2009; Klenke et al., Reference KLENKE, CONSTANTIN and WRAY2010). From an evolutionary perspective, neuropeptides that stimulate eating enhance survival during energetic challenges by inhibiting the HPG axis.
Taken together, glucocorticoids stimulate chickens to ingest more energy to cope with environmental changes, and furthermore, reproductive processes are inhibited when food availability is low or when increased energy demands are not met by compensatory food intake.
Glucocorticoids and energy distribution
The physiological mechanisms that control energy balance are integrated with those that control reproduction (Wade and Schneider, Reference WADE and SCHNEIDER1992; Wade et al., Reference WADE, SCHNEIDER and LI1996; Schneider and Watts, Reference SCHNEIDER and WATTS2002). Normally, these mechanisms optimise reproductive success under fluctuating metabolic conditions. Under stressful conditions, however, individuals should favour their own survival over that of their offspring, particularly when breeding becomes too costly in terms of body maintenance (Stearns, Reference STEARNS1992). A stressful environment can disrupt the reproduction of laying hens by triggering physiological and behavioural changes that redistribute the bodies energy reserves.
cross talk between HPA and HPG axis
The onset of breeding involves activation of the HPG axis. The hypothalamic secretion of GnRH, together with a newly discovered hypothalamic neuropeptide gonadotropin-inhibitory hormone (GnIH) which acts in opposition to GnRH-I (Tsutsui et al., Reference TSUTSUI, SAIGOH, UKENA, TERANISHI, FUJISAWA, KIKUCHI, ISHII and SHARP2000), modulates the release of the pituitary gonadotropins, luteinising hormone (LH) and follicle-stimulating hormone (FSH) (Dawson et al., Reference DAWSON, KING, BENTLEY and BALL2001), which in turn activates the gonadal development and release of sex steroids, including oestradiol and testosterone, thereby mediating the adoption of breeding behaviours (Ball, Reference BALL1993). The changes to the HPG axis during stress are dependent on the type and duration of stimulus. In the initial phase of an acute stress response, the pituitary gonadotropins (Euker et al., Reference EUKER, MEITES and RIEGLE1975; Turpen et al., Reference TURPEN, JOHNSON and DUNN1976; Mann et al., Reference MANN, EVANS, JACOBS and COLLINS1986) and sex steroid levels (Lennartsson et al., Reference LENNARTSSON, KUSHNIR, BERGQUIST, BILLIG and JONSDOTTIR2012) may increase. During prolonged periods of stress, however, the HPG axis and the production of sex steroids could be inhibited for the sake of saving energy. In laying hens, stress such as heat (Rozenboim et al., Reference ROZENBOIM, TAKO, GAL-GARBER, PROUDMAN and UNI2007), moulting and immune challenges (Sundaresan et al., Reference SUNDARESAN, ANISH, SASTRY, SAXENA, MOHAN and SAXENA2007; Ariyadi et al., Reference ARIYADI, ISOBE and YOSHIMURA2012; Zhao et al., Reference ZHAO, ZHANG, JIAO, WANG, JIANG, LUO and LIN2016), dexamethasone injection (Rories et al., Reference RORIES, LAU, FINK and SPELSBERG1989) and dietary corticosterone (Kim and Choi, Reference KIM and CHOI2014) can affect reproductive performance, including egg production, egg quality, gene expression and protein abundance in the ovary and oviduct.
The HPG axis is vulnerable to the actions of hormones associated with the activation of the HPA axis. Hypothalamic and pituitary peptides such as corticotropin-releasing hormone (CRH) and adrenocorticotropic hormone (ACTH) may act primarily at the central level (Gambacciani et al., Reference GAMBACCIANI, YEN and RASMUSSEN1986; Petraglia et al., Reference PETRAGLIA, SUTTON, VALE and PLOTSKY1987; Dobson et al., Reference DOBSON, ESSAWY and ALAM1988), while glucocorticoids may act at multiple levels to suppress the HPG axis. Strong evidence from both in vivo and in vitro systems in mammals demonstrate that glucocorticoids decrease hypothalamic GnRH release (Dubey and Plant, Reference DUBEY and PLANT1985; Oakley et al., Reference OAKLEY, BREEN, CLARKE, KARSCH, WAGENMAKER and TILBROOK2009), inhibit the basal or GnRH-stimulated release of LH from the pituitary (Saketos et al., Reference SAKETOS, SHARMA and SANTORO1993; Turner et al., Reference TURNER, HEMSWORTH, CANNY and TILBROOK1999a, b), and reduce gonadal responsiveness to LH and concentrations of LH receptors (Bambino and Hsueh, Reference BAMBINO and HSUEH1981; Hayashi and Moberg, Reference HAYASHI and MOBERG1990), or act directly on the testis/ovary by reducing steroidogenesis and/or gametogenesis (Welsh et al., Reference WELSH, BAMBINO and HSUEH1982; Whirledge and Cidlowski, Reference WHIRLEDGE and CIDLOWSKI2010). Likewise, in birds, corticosterone disturbs reproduction at the hypothalamic level by modulating GnRH/GnIH release (Son et al., Reference SON, UBUKA, NARIHIRO, FUKUDA, HASUNUMA, YAMAMOTO, BELSHAM and TSUTSUI2014), at the pituitary level by inhibiting LH secretion (Wilson and Lacassagne, 1978; Etches et al., Reference ETCHES, WILLIAMS and RZASA1984), and at the gonadal level by inhibiting testosterone /oestradiol release (Henriksen et al., Reference HENRIKSEN, GROOTHUIS and RETTENBACHER2011; Wang et al., Reference WANG, LI, SONG, GUO, JIAO, SONG and LIN2013) and ovarian function (Wilson and Lacassagne, 1978; Etches et al., Reference ETCHES, WILLIAMS and RZASA1984). We can see that glucocorticoids link the HPA axis to the HPG axis, with the secretion of glucocorticoids via the HPA axis transmitting information about environmental conditions to the HPG axis and acting at multiple possible sites along the HPG axis, which ultimately influences egg production. Thus, a corticosterone model has been widely used to study the effects of chronic and repeated stress during the rearing phase on reproductive system development (Williams et al., Reference WILLIAMS, ETCHES and RZASA1985; Shini et al., Reference SHINI, SHINI and HUFF2009), the onset of lay (Shini et al., Reference SHINI, SHINI and HUFF2009; Goutte et al., Reference GOUTTE, CLEMENT-CHASTEL, MOE, BECH, GABRIELSEN and CHASTEL2011) and performance of laying hens in subsequent laying period (Williams et al., Reference WILLIAMS, ETCHES and RZASA1985; Shini et al., Reference SHINI, SHINI and HUFF2009).
energy distribution and reproduction
Glucocorticoids trigger physiological and behavioural adjustments that shift energy investment away from reproduction and redirect it towards survival. In immature chickens, glucocorticoids stimulate energy intake (Lin et al., Reference LIN, DECUYPERE and BUYSE2004) and enhance hepatic lipogenesis and fat deposition in adipose tissues, indicating the redistribution of energy stores (Cai et al., Reference CAI, SONG, ZHANG, WANG, JIAO and LIN2009; Reference CAI, SONG, WANG, JIAO and LIN2011; Wang et al., Reference WANG, LIN, SONG and JIAO2010). In laying hens, corticosterone results in the liver and abdominal fat pad weighing significantly more whereas the ovary and oviduct is significantly regressed (Pilo et al., Reference PILO, ETCHES and GEORGE1985). We can see that glucocorticoids increase energy deposition at the expense of reproductive function. The excess energy storage in liver and adipose tissue causes deficits in the pool of fuels available for the reproductive system (Schneider, Reference SCHNEIDER2004). Specifically, glucocorticoids impede the utilisation of yolk precursors.
Glucocorticoids alter the hepatic synthesis of yolk precursors. The development of ovarian follicles is accompanied by the deposition of a large amount of yolk. During a laying cycle, follicular development is matched with a supply of yolk precursors. Triglyceride (TG) is the primary component of yolk lipid (Kuksis, Reference KUKSIS1992), which is mainly synthesised in the liver and secreted in the form of yolk-targeted very-low density lipoprotein (VLDLy). An elevated circulating VLDLy concentration was coincident with rapid yolk development and matched the variation in total yolky follicle mass (Challenger et al., Reference CHALLENGER, WILLIAMS, CHRISTIANS and VEZINA2001). VLDLy has two associated apo-lipoproteins, apo-lipoprotein B and apo-lipoprotein VLDL-II. Avian apo-lipoprotein B is a component of specialised VLDL particles that are produced by the liver in response to oestrogen (Kirchgessner et al., Reference KIRCHGESSNER, HEINZMANN, SVENSON, GORDON, NICOSIA, LEBHERZ, LUSIS and WILLIAMS1987); Apo-lipoprotein VLDL-II mRNA is specifically expressed in the liver and up-regulated after laying (Yen et al., Reference YEN, JIANG, SHEN, WONG, CHEN, CHEN, CHANG, TSAO and DING2005). The suppressed follicular development of corticosterone-treated hens is associated with the reduced availability of yolk precursor (low levels of VLDLy and vitellogenin) and the decreased expression of apo-lipoprotein B and apo-lipoprotein VLDL-II in the liver (Wang et al., Reference WANG, LI, SONG, GUO, JIAO, SONG and LIN2013). However, lipogenic activity as well as TG content in the liver was increased (Pilo et al., Reference PILO, ETCHES and GEORGE1985; Cai et al., Reference CAI, SONG, ZHANG, WANG, JIAO and LIN2009, Reference CAI, SONG, WANG, JIAO and LIN2011; Wang et al., Reference WANG, LI, SONG, GUO, JIAO, SONG and LIN2013). On the one hand, decreased apo-lipoprotein synthesis and VLDLy secretion by the liver are responsible for the decreased availability of circulating yolk precursor to oocytes; On the other hand, increased lipogenic activity with decreased lipid secretion contributes to the heavy deposition of lipid in liver.
Glucocorticoids reduce the uptake of VLDLy by follicles. In laying hens, VLDLy and vitellogenin are taken up by growing oocytes via receptor-mediated endocytosis, and both macromolecules bind to the VLDL receptors which localize on the oocyte plasma membrane (Nimpf and Schneider, Reference NIMPF and SCHNEIDER1991). Moreover, the transport of yolk precursors from the circulation to oocytes consists of five components, including granulosa cells, in which occludin- and claudin-based tight junctions are crucial for barrier function (Furuse et al., Reference FURUSE, HATA, FURUSE, YOSHIDA, HARATAKE, SUGITANI, NODA, KUBO and TSUKITA2002; Schuster et al., Reference SCHUSTER, SCHMIERER, SHKUMATAVA and KUCHLER2004; Schneider, Reference SCHNEIDER2009). Wang et al. (Reference WANG, LI, SONG, GUO, JIAO, SONG and LIN2013) observed a drop in VLDL receptor content and an increase in the expression of the tight junction proteins occludin and claudin1 in hierarchical follicles, indicating that enhanced tight junctions further prevent yolk deposition into oocytes. Furthermore, exogenous glucocorticoids perturb the reproduction of laying hens in an energy dependent-manner, that is, corticosterone suppresses follicular development when hens are in an energy deficit but not when they are in an energy sufficient state (Wang et al., Reference WANG, LI, SONG, GUO, JIAO, SONG and LIN2013).
energy redistribution in moulting
Moulting is a special period that is naturally initiated at the end of the lay cycle in birds, and is implicated in a regression of the reproductive system followed by a rejuvenation of egg-laying potential. In wild avian species, the initiation of seasonal feather moult frequently coincides with the incubation of eggs and brooding of offspring, and egg brooding is accompanied by weight loss and spontaneous anorexia (Sherry et al., Reference SHERRY, MROSOVSKY and HOGAN1980), reflecting a change in the regulation of energy balance during moulting. In commercial laying hens, induced moulting has been used as a practical procedure allowing a second or third cycle of egg production. Numerous studies show that food and water withdrawal, which is used in commercial practice to induce moulting, leads to regression of the ovary and the cessation of egg laying (Tanabe et al., Reference TANABE, OGAWA and NAKAMURA1981; Etches et al., Reference ETCHES, WILLIAMS and RZASA1984; Proszkowiec and Rzasa, Reference PROSZKOWIEC and RZASA2001). The mechanism and causes of natural moulting are not identical with those of artificial moult (Decuypere and Verheyen, Reference DECUYPERE and VERHEYEN1986), for example, prolactin is implicated in the development of broodiness and the subsequent moult when this occurs naturally, but not in the commercial laying hen induced to moult by fasting (Scanes et al., Reference SCANES, SHARP, HARVEY, GODDEN, CHADWICK and NEWCOMER1979; Verheyen et al., Reference VERHEYEN, DECUYPERE, KUHN, FONTAINE and DE GROOTE1983). In any case, moulting is a process that requires an investment of substantial amounts of energy and nutrients and is usually temporally separated from reproduction, another energy-demanding process.
A general increase in reproductive performance following induced moulting is referred to as rejuvenation, which may be mediated by an increase in tissue sensitivity or efficiency (Brake and Thaxton, Reference BRAKE and THAXTON1979). Dramatic changes take place in the liver and reproductive systems before, during and after moulting. During one year of egg production, the lipid content in the liver increases with age and shows a dramatic decrease after moulting (Garlich et al., Reference GARLICH, BRAKE, PARKHURST, THAXTON and MORGAN1984), and enhanced hepatic yolk precursor synthesis and secretion contributes to increased post-moulting laying performance (Meng et al., Reference MENG, HOU, WANG, JIAO, ZHAO, SONG and LIN2013). The loss of support provided by ovarian steroids results in regression of the oviduct and moult (Etches et al., Reference ETCHES, WILLIAMS and RZASA1984; Berry, Reference BERRY2003). Moult enhances the sensitivity of sex hormones in F1 follicles and facilitates yolk-precursor uptake by follicles, ultimately resulting in rejuvenation of reproductive performance (Meng et al., Reference MENG, HOU, WANG, JIAO, ZHAO, SONG and LIN2013). Initiation of moulting and the subsequent dynamics of energy reserves and recovery of tissue function, may be viewed as a complex physiological process, involving endocrine systems, such as corticosterone, ovarian steroids and thyroid hormones, as well as tissue structure, such as liver, ovary and oviduct (Decuypere and Verheyen, Reference DECUYPERE and VERHEYEN1986; Berry, Reference BERRY2003). These physiological mechanisms force hens to survive without food during the moult period, not only relying on catabolism of fat reserves, but also accessing resources by involution of the reproductive organs which spares other fat and protein reserves (Berry, Reference BERRY2003).
Conclusions
Stress in laying hens leads to rapid changes of endocrine physiology consistent with HPA axis activation and HPG axis deactivation. Reproductive failure and overeating are consequences of disproportionate energy storage induced by elevated corticosterone in response to stress stimulus. To attain the highest performance and optimum feed efficiency, it is important to be aware of the following factors. (1) Genetic breeding. The higher production performance and feed conversion efficiency make today's chickens more susceptible to stress than ever before. The increasing intensity of rearing makes it necessary to reconsider the selection strategy of today's commercial breeding programme in the long run. It would be interesting to systematically investigate whether strain-specific differences in stress resistance are accompanied by different function of some genes. (2) Nutritional strategies. Given the published results, stress impairs the laying performance of hens in an energy-dependent manner. The nutritional strategies of modulating energy level to satisfy the birds’ special needs during stress seem to alleviate the disadvantages of stress in the short term. (3) Production system. Over 90% of eggs all over the world are currently produced from caged layers. The research and development of new environmental control system, especially in new equipment and operation, is an urgent need for this intensive production.
Acknowledgments
This work was supported by grants from the National Key Research Programme of China (2016YFD0500510), the China Agriculture Research System (CARS-41), the National Natural Science Foundation of China (31672441, 31301993, 31272467 and 31472114), the Natural Science Foundation of Shandong Province (ZR2014CP028), the Funds of Shandong ‘Double Tops’ Programme, and the Taishan Scholars Programme (201511023).