吉林省天然阔叶混交林生态系统多功能性及驱动因素

量化天然林生态系统的多功能性,分析不同功能间的权衡-协同关系及驱动因子,对于天然林保护及修复具有重要的意义。基于吉林省第8次森林资源清查天然阔叶混交林固定样地调查数据、土壤及气候数据,选取土壤保持、涵养水源、碳储量、气候调节、土壤肥力维持、生物多样性、生产力和木材生产8个生态系统功能来表征生态系统多功能性。利用平均值法中的最大值转换法计算多功能性指数。结果表明：（1）8个功能间权衡、协同和中性关系均存在,但以协同关系为主。生物多样性除与气候调节为权衡关系外,与其他功能均为协同关系;碳储量-木材生产协同关系最强（r=0.960,P<0.01）,气候调节-涵养水源间权衡关系最强（r=-0.934,P<0.01）;（2）吉林省天然阔叶混交林的多功能性指数在0.31-0.89之间,且生物多样性和气候调节为主导功能;（3）多功能性与驱动因子的结构方程模型确定系数为R²=0.795,多功能性的驱动因子的总路径系数依次为：林分密度指数（0.752） > 平均年龄（0.375） > 年降雨量（0.365） > 树种丰富度（0.101） > 土壤pH（0.064） > 结构多样性（-0.037） > 年均温（-0.105）,林分密度是最重要的驱动因子。结果对理解天然阔叶混交林的多功能形成及经营调控有一定的意义。     

吉林蛟河针阔混交林林层结构对生产力的影响

探讨不同林层物种多样性和林分结构与森林生产力之间的关系,为森林可持续经营管理策略优化提供理论支持。以吉林蛟河30 hm²针阔混交林固定监测样地中的38432株木本植物的野外调查数据为研究对象,利用结构方程模型探讨了不同林层物种多样性和林分结构多样性对森林生产力的作用机制。结果表明：（1）林冠层林分结构多样性和生产力之间呈负相关关系,路径系数为-0.138;林下层林分结构多样性和生产力之间呈正相关关系,路径系数为0.083。（2）林冠层物种多样性和生产力之间呈正相关关系,路径系数为0.099;林下层物种多样性和生产力之间关系不显著,且物种多样性在林冠层和林下层呈正相关关系,路径系数为0.147。物种多样性和林分结构多样性与森林生产力之间不存在唯一的关系,这种关系的大小、方向和潜在机制与林层有关。为此,在森林经营过程中要充分考虑不同林层物种多样性和林分结构多样性与森林生产力之间的关系。     

A penalized structural equation modeling method accounting for secondary phenotypes for variable selection on genetically regulated expression from PrediXcan for Alzheimer's disease

As the global burden of mental illness is estimated to become a severe issue in the near future, it demands the development of more effective treatments. Most psychiatric diseases are moderately to highly heritable and believed to involve many genes. Development of new treatment options demands more knowledge on the molecular basis of psychiatric diseases. Toward this end, we propose to develop new statistical methods with improved sensitivity and accuracy to identify disease‐related genes specialized for psychiatric diseases. The qualitative psychiatric diagnoses such as case control often suffer from high rates of misdiagnosis and oversimplify the disease phenotypes. Our proposed method utilizes endophenotypes, the quantitative traits hypothesized to underlie disease syndromes, to better characterize the heterogeneous phenotypes of psychiatric diseases. We employ the structural equation modeling using the liability‐index model to link multiple genetically regulated expressions from PrediXcan and the manifest variables including endophenotypes and case‐control status. The proposed method can be considered as a general method for multivariate regression, which is particularly helpful for psychiatric diseases. We derive penalized retrospective likelihood estimators to deal with the typical small sample size issue. Simulation results demonstrate the advantages of the proposed method and the real data analysis of Alzheimer's disease illustrates the practical utility of the techniques. Data used in preparation of this article were obtained from the Alzheimer's Disease Neuroimaging Initiative database.     

A detailed model and Monte Carlo simulation for predicting DIP genome length distribution in baculovirus infection of insect cells

Baculoviruses have enormous potential for use as biopesticides to control insect pest populations without the adverse environmental effects posed by the widespread use of chemical pesticides. However, continuous baculovirus production is susceptible to DNA mutation and the subsequent production of defective interfering particles (DIPs). The amount of DIPs produced and their genome length distribution are of great interest not only for baculoviruses but for many other DNA and RNA viruses. In this study, we elucidate this aspect of virus replication using baculovirus as an example system and both experimental and modeling studies. The existing mathematical models for the virus replication process consider DIPs as a lumped quantity and do not consider the genome length distribution of the DIPs. In this study, a detailed population balance model for the cell‐virus culture is presented, which predicts the genome length distribution of the DIP population along with their relative proportion. The model is simulated using the kinetic Monte Carlo algorithm, and the results agree well with the experimental results. Using this model, a practical strategy to maintain the DIP fraction to near to its maximum and minimum limits has been demonstrated.     

黄丘区特色治理开发小流域土壤侵蚀变化对景观格局演变的响应

研究土壤侵蚀与景观格局变化的关系对小流域的治理开发具有重要的指导意义。本研究以实施退耕还林草、生态农业、生态旅游及科技示范的黄土高原安塞南沟特色治理小流域为研究对象,基于GIS平台和通用土壤流失方程,分析小流域1981—2018年景观格局和土壤侵蚀量的时空演化特征,并利用主成分回归法,从斑块类型水平和景观水平两个尺度分析土壤侵蚀模数与3类9个景观格局指标的关系。结果表明: 研究期间,在5种景观类型中,耕地和林地面积的时空变化主导了南沟小流域景观格局的演化,并且影响整个小流域的聚集分散程度;南沟小流域的土壤侵蚀量逐年减少,1981—2018年土壤侵蚀面积减少29.7%,侵蚀模数减少61.2%,且有73.4%的区域土壤侵蚀强度减轻;耕地和林地面积的变化决定了整个小流域土壤侵蚀模数的变化,其景观格局指数的变化方向与该景观类型土壤侵蚀的变化方向一致;退耕还林草工程是流域景观格局变化、土壤侵蚀减轻的主要原因,特色开发治理可以减弱局部地区土壤侵蚀强度。景观类型的合理化配置能有效地防治小流域土壤侵蚀,将其与特色治理开发相结合有助于实现小流域可持续高质量发展。     

鼎湖山南亚热带常绿阔叶林植物功能性状变异与不同垂直层次个体生长的关联

建立植物功能性状与群落动态之间的关联是功能生态学的核心问题之一。本文基于鼎湖山1.44 ha塔吊样地的两次调查数据, 通过采集样地内所有4,142株个体的6种植物功能性状, 对比分析了个体水平植物功能性状和物种水平功能性状均值对不同垂直层次(灌木层、亚冠层和林冠层)个体生长的影响差异。首先, 分析了不同垂直层次下各植物功能性状的变化趋势; 其次, 计算了不同垂直层次下各植物功能性状的种内和种间变异水平; 最后, 运用结构方程模型探讨了植物功能性状、光竞争以及地下竞争对不同垂直层次树木生长的影响。结果表明: (1)不同垂直层次下的植物功能性状表现出明显的分异, 由灌木层至林冠层, 叶面积、比叶面积和能量供求关系指数显著降低, 而叶片厚度和叶片干物质含量显著升高; (2)不同垂直层次下植物功能性状的种间变异均大于种内变异, 且林冠层的种内功能性状变异均大于灌木层和亚冠层; (3)基于个体水平植物功能性状的结构方程模型较物种水平功能性状均值对生长具有更高的解释程度, 且个体水平植物功能性状的引入更有利于提高对灌木层个体生长的预测能力; (4)光竞争和地下竞争主要通过影响功能性状间接影响植物生长。由灌木层至林冠层, 同种间的相互作用逐渐减弱, 异种间的相互作用逐渐增强。综上, 将个体水平植物功能性状纳入分析有助于更好地理解群落的结构和动态。     

Validation of the Boussinesq equation for use in traction field determination

Cell traction force plays an important role in many biological processes. Several traction force microscopy methods have been developed to determine cell traction forces based on the Boussinesq solution. This approach, however, is rooted in a half-space assumption. The purpose of this study was to determine the error induced in the half-space assumption using a finite element method (FEM). It demonstrates that displacement error between the FEM and the Boussinesq equation can be used to measure the accuracy of the Boussinesq equation, although singularity exists in the loading point. For one concentrated force, significant difference between the FEM and the Boussinesq equation occurs in the whole field; this difference decreases with an increase in the plate thickness. However, in the case of the balanced forces, the offset of the balanced forces decreases the errors in the middle area. Overall, this study demonstrates that increasing the thickness of the polyacrylamide gel is important for reducing the error of the Boussinesq equation when determining the displacement field of the gel under loads.     

MODEL DEVELOPMENT AND PROCESS OPTIMIZATION FOR SOLVENT EXTRACTION OF POLYPHENOLS FROM RED GRAPES USING BOX–BEHNKEN DESIGN

The objective of the present study is to find out the optimum extraction conditions for extraction of polyphenols from red grapes using Box–Behnken design. Red grapes polyphenols were extracted using acid–ethanol solvent at various extraction temperature (40–60°C), extraction time (20–100 min) and different solid–liquid ratio (1:5–1:15 g:ml). The effect (main and interactive) of extraction conditions on total anthocyanin, phenolic and flavonoid content were studied using Box–Behnken design (three factors at three levels). The results showed that the contribution of the quadratic model was significant for all the responses. Second-order mathematical regression models were developed and were found to fit well with observed data. Derringer's desirability function methodology was performed to find out the optimal conditions based on both individual and combinations of all responses (extraction temperature: 57°C, time: 61 min, and solid–liquid ratio: 1:8.7 g:ml) were established. At this optimal condition, the anthocyanin yield, total phenolic and flavonoid content were 73.92 mg/100 g, 221.4 mg GAE/100 g, and 79.08 mg CE/100 g, respectively. A desirability value of 0.902 was achieved at this point.     

Body masses,functional responses and predator–prey stability

The stability of ecological communities depends strongly on quantitative characteristics of population interactions (type‐II vs. type‐III functional responses) and the distribution of body masses across species. Until now, these two aspects have almost exclusively been treated separately leaving a substantial gap in our general understanding of food webs. We analysed a large data set of arthropod feeding rates and found that all functional‐response parameters depend on the body masses of predator and prey. Thus, we propose generalised functional responses which predict gradual shifts from type‐II predation of small predators on equally sized prey to type‐III functional‐responses of large predators on small prey. Models including these generalised functional responses predict population dynamics and persistence only depending on predator and prey body masses, and we show that these predictions are strongly supported by empirical data on forest soil food webs. These results help unravelling systematic relationships between quantitative population interactions and large‐scale community patterns.