Systems biology of plants : from intraspecific genome variation to computational morphodynamics

Recently, the introduction of second generation sequencing and further advance-ments in confocal microscopy have enabled system-level studies for the functional characterization of genes. The degree of complexity intrinsic to these approaches needs the development of bioinformatics methodologies and computational models for extracting meaningful biological knowledge from the enormous amount of experi¬mental data which is continuously generated. This PhD thesis presents several novel bioinformatics methods and computational models to address specific biological questions in Plant Biology by using the plant Arabidopsis thaliana as a model system. First, a spatio-temporal qualitative analysis of quantitative transcript and protein profiles is applied to show the role of the BREVIS RADIX (BRX) protein in the auxin- cytokinin crosstalk for root meristem growth. Core of this PhD work is the functional characterization of the interplay between the BRX protein and the plant hormone auxin in the root meristem by using a computational model based on experimental evidence. Hyphotesis generated by the modelled to the discovery of a differential endocytosis pattern in the root meristem that splits the auxin transcriptional response via the plasma membrane to nucleus partitioning of BRX. This positional information system creates an auxin transcriptional pattern that deviates from the canonical auxin response and is necessary to sustain the expression of a subset of BRX-dependent auxin-responsive genes to drive root meristem growth. In the second part of this PhD thesis, we characterized the genome-wide impact of large scale deletions on four divergent Arabidopsis natural strains, through the integration of Ultra-High Throughput Sequencing data with data from genomic hybridizations on tiling arrays. Analysis of the identified deletions revealed a considerable portion of protein coding genes affected and supported a history of genomic rearrangements shaped by evolution. In the last part of the thesis, we showed that VIP3 gene in Arabidopsis has an evo-lutionary conserved role in the 3' to 5' mRNA degradation machinery, by applying a novel approach for the analysis of mRNA-Seq data from random-primed mRNA. Altogether, this PhD research contains major advancements in the study of natural genomic variation in plants and in the application of computational morphodynamics models for the functional characterization of biological pathways essential for the plant. - Récemment, l'introduction du séquençage de seconde génération et les avancées dans la microscopie confocale ont permis des études à l'échelle des différents systèmes cellulaires pour la caractérisation fonctionnelle de gènes. Le degrés de complexité intrinsèque à ces approches ont requis le développement de méthodologies bioinformatiques et de modèles mathématiques afin d'extraire de la masse de données expérimentale générée, des information biologiques significatives. Ce doctorat présente à la fois des méthodes bioinformatiques originales et des modèles mathématiques pour répondre à certaines questions spécifiques de Biologie Végétale en utilisant la plante Arabidopsis thaliana comme modèle. Premièrement, une analyse qualitative spatio-temporelle de profiles quantitatifs de transcripts et de protéines est utilisée pour montrer le rôle de la protéine BREVIS RADIX (BRX) dans le dialogue entre l'auxine et les cytokinines, des phytohormones, dans la croissance du méristème racinaire. Le noyau de ce travail de thèse est la caractérisation fonctionnelle de l'interaction entre la protéine BRX et la phytohormone auxine dans le méristème de la racine en utilisant des modèles informatiques basés sur des preuves expérimentales. Les hypothèses produites par le modèle ont mené à la découverte d'un schéma différentiel d'endocytose dans le méristème racinaire qui divise la réponse transcriptionnelle à l'auxine par le partitionnement de BRX de la membrane plasmique au noyau de la cellule. Cette information positionnelle crée une réponse transcriptionnelle à l'auxine qui dévie de la réponse canonique à l'auxine et est nécessaire pour soutenir l'expression d'un sous ensemble de gènes répondant à l'auxine et dépendant de BRX pour conduire la croissance du méristème. Dans la seconde partie de cette thèse de doctorat, nous avons caractérisé l'impact sur l'ensemble du génome des délétions à grande échelle sur quatre souches divergentes naturelles d'Arabidopsis, à travers l'intégration du séquençage à ultra-haut-débit avec l'hybridation génomique sur puces ADN. L'analyse des délétions identifiées a révélé qu'une proportion considérable de gènes codant était affectée, supportant l'idée d'un historique de réarrangement génomique modelé durant l'évolution. Dans la dernière partie de cette thèse, nous avons montré que le gène VÏP3 dans Arabidopsis a conservé un rôle évolutif dans la machinerie de dégradation des ARNm dans le sens 3' à 5', en appliquant une nouvelle approche pour l'analyse des données de séquençage d'ARNm issue de transcripts amplifiés aléatoirement. Dans son ensemble, cette recherche de doctorat contient des avancées majeures dans l'étude des variations génomiques naturelles des plantes et dans l'application de modèles morphodynamiques informatiques pour la caractérisation de réseaux biologiques essentiels à la plante. - Le développement des plantes est écrit dans leurs codes génétiques. Pour comprendre comment les plantes sont capables de s'adapter aux changements environnementaux, il est essentiel d'étudier comment leurs gènes gouvernent leur formation. Plus nous essayons de comprendre le fonctionnement d'une plante, plus nous réalisons la complexité des mécanismes biologiques, à tel point que l'utilisation d'outils et de modèles mathématiques devient indispensable. Dans ce travail, avec l'utilisation de la plante modèle Arabidopsis thalicinci nous avons résolu des problèmes biologiques spécifiques à travers le développement et l'application de méthodes informatiques concrètes. Dans un premier temps, nous avons investigué comment le gène BREVIS RADIX (BRX) régule le développement de la racine en contrôlant la réponse à deux hormones : l'auxine et la cytokinine. Nous avons employé une analyse statistique sur des mesures quantitatives de transcripts et de produits de gènes afin de démontrer que BRX joue un rôle antagonisant dans le dialogue entre ces deux hormones. Lorsque ce-dialogue moléculaire est perturbé, la racine primaire voit sa longueur dramatiquement réduite. Pour comprendre comment BRX répond à l'auxine, nous avons développé un modèle informatique basé sur des résultats expérimentaux. Les simulations successives ont mené à la découverte d'un signal positionnel qui contrôle la réponse de la racine à l'auxine par la régulation du mouvement intracellulaire de BRX. Dans la seconde partie de cette thèse, nous avons analysé le génome entier de quatre souches naturelles d'Arabidopsis et nous avons trouvé qu'une grande partie de leurs gènes étaient manquant par rapport à la souche de référence. Ce résultat indique que l'historique des modifications génomiques conduites par l'évolution détermine une disponibilité différentielle des gènes fonctionnels dans ces plantes. Dans la dernière partie de ce travail, nous avons analysé les données du transcriptome de la plante où le gène VIP3 était non fonctionnel. Ceci nous a permis de découvrir le rôle double de VIP3 dans la régulation de l'initiation de la transcription et dans la dégradation des transcripts. Ce rôle double n'avait jusqu'alors été démontrée que chez l'homme. Ce travail de doctorat supporte le développement et l'application de méthodologies informatiques comme outils inestimables pour résoudre la complexité des problèmes biologiques dans la recherche végétale. L'intégration de la biologie végétale et l'informatique est devenue de plus en plus importante pour l'avancée de nos connaissances sur le fonctionnement et le développement des plantes.

Partitioning of genetic (Rapd) variability among sexes and populations of the Barn Owl (Tyto alba) in Europe

The white Barn Owl subspecies (Tyto alba alba) is found in southern Europe and the reddish-brown subspecies (T a. guttata) in northern and eastern Europe. In central Europe, the two subspecies interbreed producing a large range of phenotypic variants. Because of the different ratios of the subspecies in different geographic regions, we predict that genetic variation should be greater in Switzerland than in Hungary. We tested this hypothesis by measuring genetic variation with the RAPD method. As predicted, the genetic differentiation within a Swiss population of Barn Owls was significantly greater than the variation within a Hungarian population. This suggests that gene flow is greater in central Europe than at the eastern limit of the Barn Owl distribution in Hungary. In both countries genetic variation was more pronounced in females than in males. As in other birds, this is probably because female Barn Owls are less philopatric than males. The number of migrants between Hungary and Switzerland is ca. 1 individual per generation; if calculated separately for the sexes, then 0.525 for males and ca. I for females (Nm values). The difference in the number of migrants between genders again is likely a consequence of higher male philopatry. The sexual differentiation is greater in the Swiss population than in the Hungarian and the genetic substructuring of the populations of the species is substantial. The reason for the considerable population substructuring could be the nonmigratory behavior and socially monogamous pairing of the species, as well as the geographical barriers (Alps) between the populations examined.

Effects of fertilizer type and rate on partitioning of soil organic matter into pools of different stability

Type and rate of fertilizers influence the level of soil organic carbon (Corg) and total nitrogen (Nt) markedly, but the effect on C and N partitioning into different pools is open to question. The objectives of the present work were to: (i) quantify the impact of fertilizer type and rate on labile, intermediate and passive C and N pools by using a combination of biological, chemical and mathematical methods; (ii) explain previously reported differences in the soil organic matter (SOM) levels between soils receiving farmyard manure with or without biodynamic preparations by using Corg time series and information on SOM partitioning; and (iii) quantify the long-term and short-term dynamics of SOM in density fractions and microbial biomass as affected by fertilizer type and rate and determine the incorporation of crop residues into labile SOM fractions. Samples were taken from a sandy Cambisol from the long-term fertilization trial in Darmstadt, Germany, founded in 1980. The nine treatments (four field replicates) were: straw incorporation plus application of mineral fertilizer (MSI) and application of rotted farmyard manure with (DYN) or without (FYM) addition of biodynamic preparations, each at high (140 – 150 kg N ha-1 year-1; MSIH, DYNH, FYMH), medium (100 kg N ha-1 year-1; MSIM, DYNM, FYMM) and low (50 – 60 kg N ha-1 year-1; MSIL, DYNL, FYML) rates. The main findings were: (i) The stocks of Corg (t ha-1) were affected by fertilizer type and rate and increased in the order MSIL (23.6), MSIM (23.7), MSIH (24.2) < FYML (25.3) < FYMM (28.1), FYMH (28.1). Stocks of Nt were affected in the same way (C/N ratio: 11). Storage of C and N in the modelled labile pools (turnover times: 462 and 153 days for C and N, respectively) were not influenced by the type of fertilizer (FYM and MSI) but depended significantly (p ≤ 0.05) on the application rate and ranged from 1.8 to 3.2 t C ha 1 (7 – 13% of Corg) and from 90 to 140 kg N ha-1 (4-5% of Nt). In the calculated intermediate pool (C/N ratio 7), stocks of C were markedly higher in FYM treatments (15-18 t ha-1) compared to MSI treatments (12-14 t ha-1). This showed that differences in SOM stocks in the sandy Cambisol induced by fertilizer rate may be short-lived in case of changing management, but differences induced by fertilizer type may persist for decades. (ii) Crop yields, estimated C inputs (1.5 t ha-1 year-1) with crop residue, microbial bio¬mass C (Cmic, 118 – 150 mg kg-1), microbial biomass N (17 – 20 mg kg-1) and labile C and N pools did not differ significantly between FYM and DYN treatments. However, labile C increased linearly with application rate (R2 = 0.53) from 7 to 11% of Corg. This also applied for labile N (3.5 to 4.9% of Nt). The higher contents of Corg in DYN treatments existed since 1982, when the first sampling was conducted for all individual treatments. Contents of Corg between DYN and FYM treatments con-verged slightly since then. Furthermore, at least 30% of the difference in Corg was located in the passive pool where a treatment effect could be excluded. Therefore, the reported differences in Corg contents existed most likely since the beginning of the experiment and, as a single factor of biodynamic agriculture, application of bio-dynamic preparations had no effect on SOM stocks. (iii) Stocks of SOM, light fraction organic C (LFOC, ρ ≤ 2.0 g cm-3), light fraction organic N and Cmic decreased in the order FYMH > FYML > MSIH, MSIL for all sampling dates in 2008 (March, May, September, December). However, statistical significance of treatment effects differed between the dates, probably due to dif-ferences in the spatial variation throughout the year. The high proportion of LFOC on total Corg stocks (45 – 55%) highlighted the importance of selective preservation of OM as a stabilization mechanism in this sandy Cambisol. The apparent turnover time of LFOC was between 21 and 32 years, which agreed very well with studies with substantially longer vegetation change compared to our study. Overall, both approaches; (I) the combination of incubation, chemical fractionation and simple modelling and (II) the density fractionation; provided complementary information on the partitioning of SOM into pools of different stability. The density fractionation showed that differences in Corg stocks between FYM and MSI treatments were mainly located in the light fraction, i.e. induced by higher recalcitrance of the organic input in the FYM treatments. Moreover, the use of the combination of biological, chemical and mathematical methods indicated that effects of fertilizer rate on total Corg and Nt stocks may be short-lived, but that the effect of fertilizer type may persist for longer time spans in the sandy Cambisol.

A simple, coherent framework for partitioning uncertainty in climate predictions

A simple and coherent framework for partitioning uncertainty in multi-model climate ensembles is presented. The analysis of variance (ANOVA) is used to decompose a measure of total variation additively into scenario uncertainty, model uncertainty and internal variability. This approach requires fewer assumptions than existing methods and can be easily used to quantify uncertainty related to model-scenario interaction - the contribution to model uncertainty arising from the variation across scenarios of model deviations from the ensemble mean. Uncertainty in global mean surface air temperature is quantified as a function of lead time for a subset of the Coupled Model Intercomparison Project phase 3 ensemble and results largely agree with those published by other authors: scenario uncertainty dominates beyond 2050 and internal variability remains approximately constant over the 21st century. Both elements of model uncertainty, due to scenario-independent and scenario-dependent deviations from the ensemble mean, are found to increase with time. Estimates of model deviations that arise as by-products of the framework reveal significant differences between models that could lead to a deeper understanding of the sources of uncertainty in multi-model ensembles. For example, three models are shown diverging pattern over the 21st century, while another model exhibits an unusually large variation among its scenario-dependent deviations.

Partitioning the relative fitness effects of diet and trophic morphology in the threespine stickleback

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Niche occupancy and the relative role of micro-habitat and diet in resource partitioning among pond dwelling tadpoles.

We determined microhabitat and diet niche for tadpoles from two ponds in an agricultural landscape. Additionally, we verified the intraspecific variation in resource use, and if diet and microhabitat use were correlated. Tadpoles found in the two ponds differed in microhabitat use, because in the larger pond they explored deeper places far from the margin. There were three groups with high microhabitat niche overlap. In both ponds, plant cover was the best descriptor to explain interspecific variation in microhabitat use. Tadpoles of all species ingested mainly Bacillariophyceae and Trachellomonas however the diet differed intraspecifically in the species from the two ponds. Ten items in the temporary pond and 15 items in the permanent one were ingested by all species; however, the relative abundance of each item differed. Diet similarity was not correlated to similarity in microhabitat use. In this study, diet was as important as microhabitat use to explain resource partitioning.

Oviposition of a leaf-miner on Erythroxylum tortuosum (Erythroxylaceae) leaves: hierarchical variation of physical leaf traits

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Variance partitioning of deconstructed periphyton communities: does the use of biological traits matter?

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Disentangling the Phylogenetic and Ecological Components of Spider Phenotypic Variation

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Partitioning the relative contribution of one-phase and two-phase seed dispersal when evaluating seed dispersal effectiveness

Seed dispersal effectiveness (SDE) is a conceptual framework that aims at quantifying the contribution of seed dispersal vectors to plant fitness. While it is well recognized that diplochorous dispersal systems, characterized by two successive dispersal steps performed by two different vectors (Phase I=primary seed dispersal and Phase II=secondary seed dispersal) which are common in temperate and tropical regions, little attention has been given to distinguishing the relative contribution of one-phase and two-phase dispersal to overall SDE. This conceptual gap probably results from the lack of a clear methodology to include Phase II dispersal into the calculation of SDE and to quantify its relative contribution. We propose a method to evaluate the relative contribution of one-phase and two-phase dispersal to SDE and determine whether two seed dispersers are better than one. To do so, we used the SDE landscape and an extension of the SDE landscape, the Phase II effect landscape, which measures the direction and magnitude of the Phase II dispersal effect on overall SDE. We used simulated and empirical data from a diplochorous dispersal system in the Peruvian Amazon to illustrate this new approach. Our approach provides the relative contribution of one-phase SDE (SDE1) and two-phase SDE (SDE2) to overall SDE and quantifies how much SDE changes with the addition of Phase II dispersal. Considering that the seed dispersal process is context dependent so that Phase II depends on Phase I, we predict the possible range of variation of SDE according to the variation of the probability of Phase II dispersal. In our specific study system composed of two primate species as primary dispersal vectors and different species of dung beetles as secondary dispersal vectors, the relative contribution of SDE1 and SDE2 to overall SDE varied between plant species. We discuss the context dependency of the Phase II dispersal and the potential applications of our approach. This extension to the conceptual framework of SDE enables quantitative evaluation of the effect of Phase II dispersal on plant fitness and can be easily adapted to other biotic and/or abiotic diplochorous dispersal systems.

Soil carbon balance in a tropical grassland: Estimation of soil respiration and its partitioning using a semi-empirical model

In savannah and tropical grasslands, which account for 60% of grasslands worldwide, a large share of ecosystem carbon is located below ground due to high root:shoot ratios. Temporal variations in soil CO2 efflux (R-S) were investigated in a grassland of coastal Congo over two years. The objectives were (1) to identify the main factors controlling seasonal variations in R-S and (2) to develop a semi-empirical model describing R-S and including a heterotrophic component (R-H) and an autotrophic component (R-A). Plant above-ground activity was found to exert strong control over soil respiration since 71% of seasonal R-S variability was explained by the quantity of photosynthetically active radiation absorbed (APAR) by the grass canopy. We tested an additive model including a parameter enabling R-S partitioning into R-A and R-H. Assumptions underlying this model were that R-A mainly depended on the amount of photosynthates allocated below ground and that microbial and root activity was mostly controlled by soil temperature and soil moisture. The model provided a reasonably good prediction of seasonal variations in R-S (R-2 = 0.85) which varied between 5.4 mu mol m(-2) s(-1) in the wet season and 0.9 mu mol m(-2) s(-1) at the end of the dry season. The model was subsequently used to obtain annual estimates of R-S, R-A and R-H. In accordance with results reported for other tropical grasslands, we estimated that R-H accounted for 44% of R-S, which represented a flux similar to the amount of carbon brought annually to the soil from below-ground litter production. Overall, this study opens up prospects for simulating the carbon budget of tropical grasslands on a large scale using remotely sensed data. (C) 2012 Elsevier B.V. All rights reserved.

Genotypic variation in Actinidia deliciosa fruit size and carbohydrate content

In a global and increasingly competitive fresh produce market, more attention is being given to fruit quality traits and consumer satisfaction. Kiwifruit occupies a niche position in the worldwide market, when compared to apples, oranges or bananas. It is a fruit with extraordinarily good nutritional traits, and its benefits to human health have been widely described. Until recently, international trade in kiwifruit was restricted to a single cultivar, but different types of kiwifruit are now becoming available in the market. Effective programmes of kiwifruit improvement start by considering the requirements of consumers, and recent surveys indicate that sweeter fruit with better flavour are generally preferred. There is a strong correlation between at-harvest dry matter and starch content, and soluble solid concentration and flavour when fruit are eating ripe. This suggests that carbon accumulation strongly influences the development of kiwifruit taste. The overall aim of the present study was to determine what factors affect carbon accumulation during Actinidia deliciosa berry development. One way of doing this is by comparing kiwifruit genotypes that differ greatly in their ability to accumulate dry matter in their fruit. Starch is the major component of dry matter content. It was hypothesized that genotypes were different in sink strength. Sink strength, by definition, is the effect of sink size and sink activity. Chapter 1 reviews fruit growth, kiwifruit growth and development and carbon metabolism. Chapter 2 describes the materials and methods used. Chapter 3, 4, 5 and 6 describes different types of experimental work. Chapter 7 contains the final discussions and the conclusions Three Actinidia deliciosa breeding populations were analysed in detail to confirm that observed differences in dry matter content were genetically determined. Fruit of the different genotypes differed in dry matter content mainly because of differences in starch concentrations and dry weight accumulation rates, irrespective of fruit size. More detailed experiments were therefore carried out on genotypes which varied most in fruit starch concentrations to determine why sink strengths were so different. The kiwifruit berry comprises three tissues which differ in dry matter content. It was initially hypothesised that observed differences in starch content could be due to a larger proportion of one or other of these tissues, for example, of the central core which is highest in dry matter content. The study results showed that this was not the case. Sink size, intended as cell number or cell size, was then investigated. The outer pericarp makes up about 60% of berry weight in ‘Hayward’ kiwifruit. The outer pericarp contains two types of parenchyma cells: large cells with low starch concentration, and small cells with high starch concentration. Large cell, small cell and total cell densities in the outer pericarp were shown to be not correlated with either dry matter content or fruit size but further investigation of volume proportion among cell types seemed justified. It was then shown that genotypes with fruit having higher dry matter contents also had a higher proportion of small cells. However, the higher proportion of small cell volume could only explain half of the observed differences in starch content. So, sink activity, intended as sucrose to starch metabolism, was investigated. In transiently starch storing sinks, such as tomato fruit and potato tubers, a pivotal role in carbon metabolism has been attributed to sucrose cleaving enzymes (mainly sucrose synthase and cell wall invertase) and to ADP-glucose pyrophosphorylase (the committed step in starch synthesis). Studies on tomato and potato genotypes differing in starch content or in final fruit soluble solid concentrations have demonstrated a strong link with either sucrose synthase or ADP-glucose pyrophosphorylase, at both enzyme activity and gene expression levels, depending on the case. Little is known about sucrose cleaving enzyme and ADP-glucose pyrophosphorylase isoforms. The HortResearch Actinidia EST database was then screened to identify sequences putatively encoding for sucrose synthase, invertase and ADP-glucose pyrophosphorylase isoforms and specific primers were designed. Sucrose synthase, invertase and ADP-glucose pyrophosphorylase isoform transcript levels were anlayzed throughout fruit development of a selection of four genotypes (two high dry matter and two low dry matter). High dry matter genotypes showed higher amounts of sucrose synthase transcripts (SUS1, SUS2 or both) and higher ADP-glucose pyrophosphorylase (AGPL4, large subunit 4) gene expression, mainly early in fruit development. SUS1- like gene expression has been linked with starch biosynthesis in several crop (tomato, potato and maize). An enhancement of its transcript level early in fruit development of high dry matter genotypes means that more activated glucose (UDP-glucose) is available for starch synthesis. This can be then correlated to the higher starch observed since soon after the onset of net starch accumulation. The higher expression level of AGPL4 observed in high dry matter genotypes suggests an involvement of this subunit in drive carbon flux into starch. Changes in both enzymes (SUSY and AGPse) are then responsible of higher starch concentrations. Low dry matter genotypes showed generally higher vacuolar invertase gene expression (and also enzyme activity), early in fruit development. This alternative cleavage strategy can possibly contribute to energy loss, in that invertases’ products are not adenylated, and further reactions and transport are needed to convert carbon into starch. Although these elements match well with observed differences in starch contents, other factors could be involved in carbon metabolism control. From the microarray experiment, in fact, several kinases and transcription factors have been found to be differentially expressed. Sink strength is known to be modified by application of regulators. In ‘Hayward’ kiwifruit, the synthetic cytokinin CPPU (N-(2-Chloro-4-Pyridyl)-N-Phenylurea) promotes a dramatic increase in fruit size, whereas dry matter content decreases. The behaviour of CPPU-treated ‘Hayward’ kiwifruit was similar to that of fruit from low dry matter genotypes: dry matter and starch concentrations were lower. However, the CPPU effect was strongly source limited, whereas in genotype variation it was not. Moreover, CPPU-treated fruit gene expression (at sucrose cleavage and AGPase levels) was similar to that in high dry matter genotypes. It was therefore concluded that CPPU promotes both sink size and sink activity, but at different “speeds” and this ends in the observed decrease in dry matter content and starch concentration. The lower “speed” in sink activity is probably due to a differential partitioning of activated glucose between starch storage and cell wall synthesis to sustain cell expansion. Starch is the main carbohydrate accumulated in growing Actinidia deliciosa fruit. Results obtained in the present study suggest that sucrose synthase and AGPase enzymes contribute to sucrose to starch conversion, and differences in their gene expression levels, mainly early in fruit development, strongly affect the rate at which starch is therefore accumulated. This results are interesting in that starch and Actinidia deliciosa fruit quality are tightly connected.

Extension and Partitioning in an Oblique Subduction Zone, New Zealand: Constraints from Three-Dimensional Numerical Modeling

Contraction, strike slip, and extension displacements along the Hikurangi margin northeast of the North Island of New Zealand coincide with large lateral gradients in material properties. We use a finite- difference code utilizing elastic and elastic-plastic rheologies to build large- scale, three-dimensional numerical models which investigate the influence of material properties on velocity partitioning within oblique subduction zones. Rheological variation in the oblique models is constrained by seismic velocity and attenuation information available for the Hikurangi margin. We compare the effect of weakly versus strongly coupled subduction interfaces on the development of extension and the partitioning of velocity components for orthogonal and oblique convergence and include the effect of ponded sediments beneath the Raukumara Peninsula. Extension and velocity partitioning occur if the subduction interface is weak, but neither develops if the subduction interface is strong. The simple mechanical model incorporating rheological variation based on seismic observations produces kinematics that closely match those published from the Hikurangi margin. These include extension within the Taupo Volcanic Zone, uplift over ponded sediments, and dextral contraction to the south.

Temporal variation in saproxylic beetle assemblages in a Mediterranean ecosystem

One of the main challenges in biological conservation has been to understand species distribution across space and time. Over the last decades, many diversity and conservation surveys have been conducted that have revealed that habitat heterogeneity acts as a major factor that determines saproxylic assemblages. However, temporal dynamics have been poorly studied, especially in Mediterranean forests. We analyzed saproxylic beetle distribution at inter and intra-annual scales in a “dehesa” ecosystem, which is a traditional Iberian agrosilvopastoral ecosystem that is characterized by the presence of old and scattered trees that dominate the landscape. Significant differences in effective numbers of families/species and species richness were found at the inter-annual scale, but this was not the case for composition. Temperature and relative humidity did not explain these changes which were mainly due to the presence of rare species. At the intra-annual scale, significant differences in the effective numbers of families/species, species richness and composition between seasons were found, and diversity partitioning revealed that season contributed significantly to gamma-diversity. Saproxylic beetle assemblages exhibited a marked seasonality in richness but not in abundance, with two peaks of activity, the highest between May and June, and the second between September and October. This pattern is mainly driven by the seasonality of the climate in the Mediterranean region, which influences ecosystem dynamics and imposes a marked seasonality on insect assemblages. An extended sampling period over different seasons allowed an overview of saproxylic dynamics, and revealed which families/species were restricted to particular seasons. Recognizing that seasons act as a driver in modelling saproxylic beetle assemblages might be a valuable tool in monitoring and for conservation strategies in Mediterranean forests.

The effects of habitat fragmentation due to forestry plantation establishment on the demography and genetic variation of a marsupial carnivore, Antechinus agilis

We conducted a demographic and genetic study to investigate the effects of fragmentation due to the establishment of an exotic softwood plantation on populations of a small marsupial carnivore, the agile antechinus (Antechinus agilis), and the factors influencing the persistence of those populations in the fragmented habitat. The first aspect of the study was a descriptive analysis of patch occupancy and population size, in which we found a patch occupancy rate of 70% among 23 sites in the fragmented habitat compared to 100% among 48 sites with the same habitat characteristics in unfragmented habitat. Mark-recapture analyses yielded most-likely population size estimates of between 3 and 85 among the 16 occupied patches in the fragmented habitat. Hierarchical partitioning and model selection were used to identify geographic and habitat-related characteristics that influence patch occupancy and population size. Patch occupancy was primarily influenced by geographic isolation and habitat quality (vegetation basal area). The variance in population size among occupied sites was influenced primarily by forest type (dominant Eucalyptus species) and, to a lesser extent, by patch area and topographic context (gully sites had larger populations). A comparison of the sex ratios between the samples from the two habitat contexts revealed a significant deficiency of males in the fragmented habitat. We hypothesise that this is due to male-biased dispersal in an environment with increased dispersal-associated mortality. The population size and sex ratio data were incorporated into a simulation study to estimate the proportion of genetic diversity that would have been lost over the known timescale since fragmentation if the patch populations had been totally isolated. The observed difference in genetic diversity (gene diversity and allelic richness at microsatellite and mitochondrial markers) between 16 fragmented and 12 unfragmented sites was extremely low and inconsistent with the isolation of the patch populations. Our results show that although the remnant habitat patches comprise approximately 2% of the study area, they can support non-isolated populations. However, the distribution of agile antechinus populations in the fragmented system is dependent on habitat quality and patch connectivity. (C) 2004 Elsevier Ltd. All rights reserved.