高山松（Pinus densata）同倍性杂交物种形成的生态遗传学研究

种间杂交在自然界中普遍存在。杂交可以为群体带来多方面的遗传改变，具有重要的进化生物学意义。越来越多的证据表明，同倍性杂交物种形成是一种重要物种形成途径，在动、植物物种进化中扮演着重要的角色。同倍性杂交物种形成意味着一个稳定的、可育的、与亲本种间形成生殖隔离的杂种群体建立起来，但不改变染色体的数目。这种生殖隔离被认为是通过快速的染色体重组、生态分化、或空间分离形成的。对于现存杂交种的对比分析和对同倍性杂交物种形成的理论模拟都表明，适应性分化和生态转变在同倍性杂交物种形成中起着至关重要的作用。来自于生态成种的研究也表明生态选择自身就可以导致生殖隔离的形成和新物种的形成。尽管如此，我们对同倍性杂交物种形成过程的研究仍然较多地关注于染色体重组及相应的遗传效应；对于自然选择，特别是生态选择在同倍性杂交成种过程中的作用模式缺乏全面地了解。 本研究以同倍性杂交种－高山松为研究对象，利用地理信息学技术、基于天然群体取样的性状对比分析和人工杂交手段，全面地定量描述了高山松及其亲本种间的生态分化、生态转变下的天然群体生殖适应性、以及生态选择下的形态变异和种间杂交可配性。以此探讨了高山松的杂交成种过程，生态转变和自然选择的重要性以及它们与杂交种适应性、形态变异和生殖隔离之间的关系。利用地理信息学手段，对高山松及其亲本种间的生态分化对比分析表明：高山松与其亲本种间形成了明显的生态转变；在自然选择的作用下，高山松已经具备了对于青藏高原东南缘横断山区高山地带独特的生态适应性；相对于亲本种，它的适应性特征体现为混合型，在某些性状上它高出亲本种的适应性，在另一些性状上，它体现出居中、较低或者类似于亲本种适应性；高山松适生区由于地形地貌的复杂性，往往不连续，呈间断状。基于四个反映结实效率的性状、小孢子母细胞减数分裂异常现象及花粉萌发率在种间的对比分析表明，这些反映生殖状态的指标，在三个种及群体间并没有显著差异；相对于亲本种的天然群体，居于青藏高原的高山松群体，表现出了正常的授粉效率和生殖适应性。高山松天然群体正常的生殖效率表明，它经历了生态转变过程中的选择作用，已经具备了在新生境下，繁殖发展下去的生殖适应性。 在大规模群体取样的基础上，通过对反映植物个体发育和营养状态的球果和针叶形态性状的对比分析，我们发现，高山松及其亲本种间的形态变异存在以下特征：第一，绝大多数性状，在种间和种内群体间都存在着显著的差异，并且大部分性状的差异主要体现在种间；第二，高山松在大部分形态性状上介于两个亲本种之间，少数性状表现为超亲分离、类似于某个亲本或无种间差别；第三，在多数性状上，和亲本种相比，高山松群体都拥有较大的变异幅度。第四，在临近亲本种的群体中，高山松的形态特征更类似于所临近的亲本种。上述形态变异特征及其与气象、土壤、经纬度、海拔等因子间普遍的相关性，意味着生态选择在高山松成种过程中起了重要作用。总的来说，形态性状变异的总体模式揭示出：高山松在选择作用下的遗传和表型的精炼过程中已经形成了自身适应于高原环境的、稳定而有效的形态发育模式和资源利用方式，甚至是恰当的防御机制。 本研究通过云南松×油松的单交实验和油松×高山松部分双列杂交实验，检验了高山松及其亲本种间的杂交可配性。结果表明，就平均水平而言，高山松及其亲本种间的种间杂交结实可配性低于种内杂交，但个别种间杂交组合的结实可配性可能高于种间杂交；高山松及其亲本种间正在或已经形成了部分的合子后生殖障碍；现存的高山松群体及其亲本种间的生态地理隔离可能是很有效的生殖隔离形式，一方面做为合子前生殖隔离，它能避免种间相互传粉、阻碍种间杂交的发生，另一方面做为合子后生殖隔离，它能清除掉可能存在的劣势回交产物。 综合已有的证据，我们提出了高山松同倍性杂交物种形成模型：高山松的同倍性杂交成种是与青藏高原隆升过程以及隆升过程中所产生的全新生态空间紧密相关的；青藏高原隆升导致了生态环境的剧变，创造了全新的生态位，在选择的作用下，拥有综合了双亲遗传基础的杂种基因型开始占领这个全新的生态位；在分化的生态选择的持续作用下，杂种基因型在新的生境下稳定下来，形成了独立于亲本种的资源利用方式、形态特征和繁育系统，具备了完善的生殖和生态适应性；此时，一个与亲本种间以生态隔离存在的同倍性杂交种－高山松形成了；随着分化的生态选择的持续进行,高山松与亲本种间的内在生殖隔离不断加强。

高山松的二倍体杂种起源-来自三个基因组的证据

利用松科植物特殊的遗传体系（叶绿体基因组一父系遗传、线粒体基因组—母系遗传、核基因组一双亲遗传），我们对高山松及其两个亲本种进行了广泛的群体取样，通过线粒体基因nadl、叶绿体基因rbcL和trnL-F基因间区以及低拷贝核基因4CL的序列分析或PCR-RFLP分析，为高山松同倍体杂种起源假说提供了翔实的遗传学证据，同时在个体水平上探讨了高山松不同群体的遗传组成、群体遗传结构、基因交流方向、群体建立过程以及杂种基因组的进化。具体结果如下： 1．细胞质基因组分析 1)线粒体基因nudl分析 本研究对油松、高山松和云南松的19个群体、295个个体的线粒体基因nadl的一个内含子进行了序列分析或PCR-RFLP分析，共检测到3种线粒体DNA单倍型-A、B和C。油松所有的取样群体仅含单倍型A;除BX群体外，所有的云南松群体仅含单倍型B; 10个高山松群体中，5个群体固定单倍型A，4个群体固定单倍型B，1个群体(ZD)分布有A和B两种单倍型。2)叶绿体rbcL基因分析 对同一组群体的rbcL基因进行序列分析或PCR-RFLP分析，共检测到两个变异位点和三种叶绿体单倍型（TT、TC和GC）。TT和GC分别是油松和云南松种特异性叶绿体单倍型，而在高山松群体里则三种单倍型均有分布，而且TC单倍型广泛地分布在7个杂种群体中，该单倍型很可能来源于点突变或第三个已灭绝的亲本。rbcL基因检测到的高山松群体分化系数很高(Gst=0.533)。 3)叶绿体trn L-F区序列分析 叶绿体trnL-F分子标记检测到的不同单倍型的差异主要是由引物“e”下游120碱基处一个多聚T结构的长度变异所致（叶绿体SSR位点）。10个高山松群体中共检测到5种叶绿体单倍型，其中两种主要的单倍型(9T和11T)分别为油松和云南松的种特异性单倍型，其他单倍型均为非典型单倍型。群体遗传结构分析表明：杂种群体表现最高的遗传多样性，而且trnL-F分析得到的高山松群体的分化系数也很高( Gst=0.443)。 总之，对高山松、油松和云南松的同一组群体取样进行的细胞质基因组分析表明：高山松群体分布有油松和云南松种特异性的线粒体和叶绿体单倍型，该细胞质DNA单倍型的地理分布为假说“高山松为油松和云南松的的二倍体杂种”提供了翔实的遗传学证据。油松和云南松在不同的杂种群体中分别做父本和母本，即两亲本在杂交过程中发生了双向基因交流。群体遗传结构分析发现高山松群体表现最高的遗传多样性，而且群体间的分化系数很高。不同的杂种群体在遗传组成上的差异表明他们经历过不同的建立和进化历史。从线粒体和叶绿体单倍型的地理分布可以看出杂种群体的建立曾经历强烈的奠基者效应和回交。青藏高原的隆升对高山松的起源、杂种群体的适应辐射以及保持产生了重要的影响。川西南和滇西北作为青藏高原的东边边界，很可能是当初云南松和油松分布的重叠区及杂交地带，即高山松的起源地。 2．核基因4CL分析 对高山松、油松和云南松的19个群体、32个个体的低拷贝核基因4CL进行了克隆及序列分析，获得的78条序列可分为两种类型（类型A和类型B）。这两种类型明显的差别是类型A相对于类型B在内含子区有- 20bp的缺失。以华山松的3条序列为外类群，对得到的78条序列进行基因谱系分析，发现所有的序列分成明显的两支，分别对应于类型A和类型B，而且每一支均包含三个种的部分序列，表明4CL基因在这三个种分化之前就已发生重复。另一个明显的特点是某个种的一条序列与另一个种的序列比其与同种的其他序列关系更近，可能因基因交流（杂交和渐渗）、非共祖、致同进化和重组等进化事件所致。三种松树中共检测到4CL基因序列的两种类型和六个亚类型，高山松群体中没有发现杂种独特的类型或亚类型。高山松和云南松共享三种序列亚类型以及最多的序列多态性，表明这两个种之间曾存在广泛的基因交流。

高山松的群体历史及其适应性进化研究

物种形成一直以来是进化生物学的一个中心议题。因为杂交能够快速创造出新的遗传变异和提高基因组重组率，加速物种形成和适应性进化，所以杂交物种形成是其重要的组成部分。但是直到最近几年，人们才逐渐意识到同倍体杂种形成的重要性和普遍性。此外，人工杂交种群结果证实了杂种的适应性起源于亲本基因的超亲分离现象，这些基因之间通过加和效应或上位效应相互连锁，并在杂种后代中逐渐固定，因此，探讨维持这种适应机制的杂种的连锁不平衡模式就显得尤为重要。到目前为止，关于自然杂种的核苷酸多态性和连锁不平衡结构的报道非常稀见，对杂种的进化历史知之更少。高山松能生长于亲本种不能正常繁殖生长的青藏高原上，这种极端环境为杂种的适应性进化和与亲本的生态隔离造就了绝佳的条件，其生态和生殖的稳定性极大地方便了我们研究植物物种形成和适应性进化的遗传机理，对高山松开展的生态学、遗传学的系列分析已使高山松成为同倍体杂种形成的经典范例。 该论文通过随机挑选7个核基因位点对5个高山松群体、3个油松群体和3个云南松群体的164个大配子体单倍体基因组进行了核苷酸多态性分析。研究发现：所有基因座位的单倍型组成和基因谱系的拓扑结构都支持高山松的多次杂交起源。高山松的平均核苷酸多态性与油松相当，θW达到0.0107;比云南松高出一倍;高于已报道的裸子植物类群遗传变异水平。高山松如此高的遗传多态性与其杂交特性和有效群体大有关。我们运用分子钟和共祖模拟分析，推测出高山松的有效群体为7.32 x 105。等位基因共祖时间模拟分析发现，杂交过程早于青藏高原的隆升，也即是说，在杂种稳定成种以前，两个亲本种之间存在广泛的渐渗杂交。此外，高山松杂种的群体间分化严重，并且在不同的群体中，我们找到了多个偏离中性的基因座位。由此说明，高山松复杂的进化历史和适应的策略区域化。 在显著偏离中性的基因座位上，选取群体历史清晰的高山松群体，对其全基因序列多态性和连锁不平衡结构进行了深入调查，以期找到选择作用位点或区域。Ara-like和Dhn1两个基因在高山松与亲本之间都没有发现固定变异。在高山松中，Ara-like基因和Dhn1基因核苷酸多态性式样和连锁不平衡结构存在明显差异，这说明在两个基因上选择作用的程度和方向不尽相同。 我们对Dhn1基因的多态性分布、LD结构和中性偏离水平进行了分析，结果表明，该基因可能在PdNX和PdLZ群体中受到平衡选择的作用。从Dhn1基因谱系结构可以看出有两种来源的等位基因造成这种平衡多态性，一种属于祖先类型，另一种是从云南松继承衍生而来。这两种等位基因之间的分化很大，极有可能造成编码蛋白的亚功能化。事实上，Dhn1编码的脱水蛋白在植物对环境的抗逆过程中发挥着重要作用。由此我们推测高山松在高海拔极端环境中的适应性进化与Dhn1基因的平衡作用有关。 对Ara-like和Dhn1基因进行HKA检测，结果表明，与Dhn1基因相比，在Ara-like基因上，双维管束亚属与单维管束亚属之间存在显著分化，沉默突变位点的分化Ksil达到了0.1392，远远高于平均水平0.0508;在双维管束亚属共发现了43个固定突变位点，其中有6个能导致氨基酸突变，它们有可能导致了Ara-like基因功能的分化。结合裸子植物近缘物种间共享多态性的普遍性，我们推测Ara-like基因可能在单维管束亚属中的进化速率加快，暗示其在单维管束亚属中的适应性进化。 最后，基因内连锁不平衡分析结果显示，随机筛选的基因座位之间不存在连锁不平衡。高山松平均的基因内LD程度非常低，仅有18%的信息位点之间显著连锁。平均基因内LD在油松和高山松中的衰减速率很快，尤其在杂种中下降最快，在不到200bp以内就降到0.1以下。 LD的结果印证了高山松的有效群体大和多次起源特性。另外，我们也怀疑杂交物种形成过程中染色体组的重组和重排频繁发生，也是造成自然杂种现有群体的LD水平低的一个主要原因。 通过LD衰减曲线估计高山松的平均单位重组率,比亲本油松高17倍，比云南松高45倍。这个结果表明，杂种共适应的基因间要维持超亲分离，需要强烈的自然选择压力，才能保持它们之间的连锁不平衡。

中国特有针叶树种对气候变化的敏感性研究

About 214 trees in 9 sampling sites, representing 5 endemic conifer species, were collected from the western Sichuan Province and eastern Qinghai Province, China. In this study, structure we try to investigate tree-ring sensitivity to climate in order to obtain primary information of reconstructing past climate from the trees in this region. All the 5 species present distinct ring boundaries ^ few ABS(absent rings) and available for cross-dating,which are all past the test by program COFECHA. Statistics for all the 8 tree-ring width residual chronologies present significant inter-correlation between series and high values of mean sensitivity. As well as the maximum latewood density of Picea crassifolia Kom and Pinus densata Mast. These results indicate usefulness of these chronologies for dendrochronological studies. Pearson correlation analyses were applied to provide a basic estimate of the causal relationships between tree-ring width and climate factors. We found some significant relationships between tree-ring width> maximum density and temperature as well as precipitation. Especially, there is high correlation between the maximum density of the Picea crassifolia Kom and the index of moisture, the ratio of precipitation and temperature, which can indicate well the climate; however the higher correlation can be see between the maximum density of Pinus densata Mast and the total temperature from June to September. Regardless of tree species, chronologies in our study region presented accordant variations of which may reveal strong common climate signal. Thus these chronologies are shown to be dependable for building tree-ring network in the nearly future. However, there are limitations in this study, only monthly mean of temperature and precipitation were available. Also, for this typical subtropical mountain system, meteorological stations are usually located in valley and biased to represent moisture conditions on the slopes. Thus the estimation of precipitation both in temporal and spatial domain was rather restricted. Further study, such as wood anatomy, physiology and densitometry, are needed for better understanding the environmental and climatic history in this area.

Long-tailed bats' use of a Pinus radiata stand in Kinleith Forest : Recommendations for monitoring

Targeted monitoring of threatened species within plantations is becoming more important due to forest certification programmes’ requirement to consider protection of threatened species, and to increase knowledge of the distribution of species. To determine patterns of long-tailed bat (Chalinolobus tuberculatus) activity in different habitat structures, with the aim of improving the likelihood of detection by targeting monitoring, we monitored one stand of 26 year-old Pinus radiata over seven months between December 2007 and June 2008 in Kinleith Forest, an exotic plantation forest centred around Tokoroa, South Waikato, New Zealand. Activity was determined by acoustic recording equipment, which is able to detect and record bats’ echolocation calls. We monitored activity from sunset to sunrise along a road through the stand, along stand edges, and in the interior of the stand. Bats were recorded on 80% of the 35 nights monitored. All activity throughout the monitoring period was detected on the edge of the stand or along the road. No bats were detected within the interior of the stand. Bat activity was highest along the road through the stand (40.4% of all passes), followed by an edge with stream running alongside (35.2%), along the road within a skidsite (19.8%), and along an edge without a stream (4.6%). There was a significant positive relationship between bat pass rate (bat passes h-1) and the feeding buzz rate (feeding buzzes h-1) indicating that bat activity was associated with feeding and not just commuting. Bat feeding activity was also highest along the road through the stand (59.2% of feeding buzzes), followed by the road within the skidsite (30.6%), and along the stream-side edge (10.2%). No feeding buzzes were recorded in either the interior or along the edge without the stream. Differences in overall feeding activity were significant only between the road and edge and between edges with and without a stream. Bat activity was detected each month and always by the second night of monitoring, and in this stand was highest during April. We recommend targeted monitoring for long-tailed bats be focused on road-side and stand edge habitat, and along streams, and that monitoring take place for at least three nights to maximise probability of detection.

Runoff, sediment loss and water quality from forest roads in a southeast Queensland coastal plain Pinus plantation

Runoff and sediment loss from forest roads were monitored for a two-year period in a Pinus plantation in southeast Queensland. Two classes of road were investigated: a gravelled road, which is used as a primary daily haulage route for the logging area, and an ungravelled road, which provides the main access route for individual logging compartments and is intensively used as a haulage route only during the harvest of these areas (approximately every 30 years). Both roads were subjected to routine traffic loads and maintenance during the study. Surface runoff in response to natural rainfall was measured and samples taken for the determination of sediment and nutrient (total nitrogen, total phosphorus, dissolved organic carbon and total iron) loads from each road. Results revealed that the mean runoff coefficient (runoff depth/rainfall depth) was consistently higher from the gravelled road plot with 0.57, as compared to the ungravelled road with 0.38. Total sediment loss over the two-year period was greatest from the gravelled road plot at 5.7 t km−1 compared to the ungravelled road plot with 3.9 t km−1. Suspended solids contributed 86% of the total sediment loss from the gravelled road, and 72% from the ungravelled road over the two years. Nitrogen loads from the two roads were both relatively constant throughout the study, and averaged 5.2 and 2.9 kg km−1 from the gravelled and ungravelled road, respectively. Mean annual phosphorus loads were 0.6 kg km−1 from the gravelled road and 0.2 kg km−1 from the ungravelled road. Organic carbon and total iron loads increased in the second year of the study, which was a much wetter year, and are thought to reflect the breakdown of organic matter in roadside drains and increased sediment generation, respectively. When road and drain maintenance (grading) was performed runoff and sediment loss were increased from both road types. Additionally, the breakdown of the gravel road base due to high traffic intensity during wet conditions resulted in the formation of deep (10 cm) ruts which increased erosion. The Water Erosion Prediction Project (WEPP):Road model was used to compare predicted to observed runoff and sediment loss from the two road classes investigated. For individual rainfall events, WEPP:Road predicted output showed strong agreement with observed values of runoff and sediment loss. WEPP:Road predictions for annual sediment loss from the entire forestry road network in the study area also showed reasonable agreement with the extrapolated observed values.

Sediment and nutrient losses from exotic Pinus plantation management operations under simulated rainfall

Rainfall simulation experiments were carried out to measure runoff and soil water fluxes of suspended solids, total nitrogen, total phosphorus, dissolved organic carbon and total iron from sites in Pinus plantations on the coastal lowlands of south-eastern Queensland subjected to various operations (treatments). The operations investigated were cultivated and nil-cultivated site preparation, fertilised site preparation, clearfall harvesting and prescribed burning; these treatments were compared with an 8-y-old established plantation. Flow-weighted mean concentrations of total nitrogen and total phosphorus in surface runoff from the cultivated and nil-cultivated site-preparation, clearfall harvest, prescribed burning and 8-y-old established plantation treatments were very similar. However, both the soil water and the runoff from the fertilised site preparation treatment contained more nitrogen (N) and phosphorus (P) than the other treatments - with 3.10 mg N L-1 and 4.32 mg P L-1 (4 and 20 times more) in the runoff. Dissolved organic carbon concentrations in runoff from the nil-cultivated site-preparation and prescribed burn treatments were elevated. Iron concentrations were highest in runoff from the nil-cultivated site-preparation and 8-y-old established plantation treatments. Concentrations of suspended solids in runoff were higher from cultivated site preparation and prescribed burn treatments, and reflect the great disturbance of surface soil at these sites. The concentrations of all analytes were highest in initial runoff from plots, and generally decreased with time. Total nitrogen (mean 7.28, range 0.11-13.27 mg L-1) and total phosphorus (mean 11.60, range 0.06-83.99 mg L-1) concentrations in soil water were between 2 and 10 times greater than in surface runoff, which highlights the potential for nutrient fluxes in interflow (i.e. in the soil above the water table) through the general plantation area. Implications in regard to forest management are discussed, along with results of larger catchment-scale studies.

Congruence in QTL for adventitious rooting in Pinus elliottii x Pinus caribaea hybrids resolves between and within-species effects

Targeting between-species effects for improvement in synthetic hybrid populations derived from outcrossing parental tree species may be one way to increase the efficacy and predictability of hybrid breeding. We present a comparative analysis of the quantitative trait loci (QTL) which resolved between from within-species effects for adventitious rooting in two populations of hybrids between Pinus elliottii and P. caribaea, an outbred F1 (n=287) and an inbred-like F2 family (n=357). Most small to moderate effect QTL (each explaining 2-5% of phenotypic variation, PV) were congruent (3 out of 4 QTL in each family) and therefore considered within-species effects as they segregated in both families. A single large effect QTL (40% PV) was detected uniquely in the F2 family and assumed to be due to a between-species effect, resulting from a genetic locus with contrasting alleles in each parental species. Oligogenic as opposed to polygenic architecture was supported in both families (60% and 20% PV explained by 4 QTL in the F 2 and F1 respectively). The importance of adventitious rooting for adaptation to survive water-logged environments was thought in part to explain oligogenic architecture of what is believed to be a complex trait controlled by many hundreds of genes.

Performance of seven hardwood species underplanted to Pinus elliottii in south-east Queensland

Seven hardwood species were tested as underplants under Pinus elliottii plantations on the coastal lowlands of south-east Queensland. The species tested were: Flindersia brayleyana (F. Muell) (Queensland maple), F. australis (R. Br.), (crow's ash), Swietenia macrophylla (King) (American mahogany), Grevillea robusta (A. cunn) (southern silky oak), Elaeocarpus grandis (F. Muell) (silver quandong), F. ifflaiana (F. Meull) (Cairns hickory) and Ceratopetalum apetalum (D. Don) (coachwood). Most species (except E. grandis) established successfully but slowly. Underplants suffered 9-16% mortality during thinning of the overstorey. By 2004 when aged c. 38 years, four underplanted species; F. brayleyana, S. macrophylla, F. ifflaiana and E. grandis, had attained predominant heights of 20 m and mean diameter at breast height of 25 cm or better. The presence of underplants increased total site productivity by up to 23% and did not have any detrimental effect on the development of the overwood.This experiment has demonstrated that some rainforest species will survive and grow healthily as underplants in exotic pine plantations plus produce small merchantable logs within a 38 year rotation. The results also indicated the importance of correct species selection if an underplanting option is to be pursued as some species have been a complete failure (notably G. robusta).

Spatial model of site index based on Y-ray spectrometry and a digital elevation model for two Pinus species in Tuan Toolara State Forest, Queensland, Australia

Site index prediction models are an important aid for forest management and planning activities. This paper introduces a multiple regression model for spatially mapping and comparing site indices for two Pinus species (Pinus elliottii Engelm. and Queensland hybrid, a P. elliottii x Pinus caribaea Morelet hybrid) based on independent variables derived from two major sources: g-ray spectrometry (potassium (K), thorium (Th), and uranium (U)) and a digital elevation model (elevation, slope, curvature, hillshade, flow accumulation, and distance to streams). In addition, interpolated rainfall was tested. Species were coded as a dichotomous dummy variable; interaction effects between species and the g-ray spectrometric and geomorphologic variables were considered. The model explained up to 60% of the variance of site index and the standard error of estimate was 1.9 m. Uranium, elevation, distance to streams, thorium, and flow accumulation significantly correlate to the spatial variation of the site index of both species, and hillshade, curvature, elevation and slope accounted for the extra variability of one species over the other. The predicted site indices varied between 20.0 and 27.3 m for P. elliottii, and between 23.1 and 33.1 m for Queensland hybrid; the advantage of Queensland hybrid over P. elliottii ranged from 1.8 to 6.8 m, with the mean at 4.0 m. This compartment-based prediction and comparison study provides not only an overview of forest productivity of the whole plantation area studied but also a management tool at compartment scale.