Interaction of Plectin with Keratins 5 and 14: Dependence on Several Plectin Domains and Keratin Quaternary Structure.

Plectin, a cytolinker of the plakin family, anchors the intermediate filament (IF) network formed by keratins 5 and 14 (K5/K14) to hemidesmosomes, junctional adhesion complexes in basal keratinocytes. Genetic alterations of these proteins cause epidermolysis bullosa simplex (EBS) characterized by disturbed cytoarchitecture and cell fragility. The mechanisms through which mutations located after the documented plectin IF-binding site, composed of the plakin-repeat domain (PRD) B5 and the linker, as well as mutations in K5 or K14, lead to EBS remain unclear. We investigated the interaction of plectin C terminus, encompassing four domains, the PRD B5, the linker, the PRD C, and the C extremity, with K5/K14 using different approaches, including a rapid and sensitive fluorescent protein-binding assay, based on enhanced green fluorescent protein-tagged proteins (FluoBACE). Our results demonstrate that all four plectin C-terminal domains contribute to its association with K5/K14 and act synergistically to ensure efficient IF binding. The plectin C terminus predominantly interacted with the K5/K14 coil 1 domain and bound more extensively to K5/K14 filaments compared with monomeric keratins or IF assembly intermediates. These findings indicate a multimodular association of plectin with K5/K14 filaments and give insights into the molecular basis of EBS associated with pathogenic mutations in plectin, K5, or K14 genes.Journal of Investigative Dermatology advance online publication, 10 July 2014; doi:10.1038/jid.2014.255.

Fine-scale population structure analysis of seven local Swiss sheep breeds using genome-wide SNP data

As part of the global sheep Hapmap project, 24 individuals from each of seven indigenous Swiss sheep breeds (Bundner Oberländer sheep (BOS), Engadine Red sheep (ERS), Swiss Black-Brown Mountain sheep (SBS), Swiss Mirror sheep (SMS), Swiss White Alpine (SWA) sheep, Valais Blacknose sheep (VBS) and Valais Red sheep (VRS)), were genotyped using Illumina’s Ovine SNP50 BeadChip. In total, 167 animals were subjected to a detailed analysis for genetic diversity using 45 193 informative single nucleotide polymorphisms. The results of the phylogenetic analyses supported the known proximity between populations such as VBS and VRS or SMS and SWA. Average genomic relatedness within a breed was found to be 12 percent (BOS), 5 percent (ERS), 9 percent (SBS), 10 percent (SMS), 9 percent (SWA), 12 percent (VBS) and 20 percent (VRS). Furthermore, genomic relationships between breeds were found for single individuals from SWA and SMS, VRS and VBS as well as VRS and BOS. In addition, seven out of 40 indicated parent–offspring pairs could not be confirmed. These results were further supported by results from the genome-wide population cluster analysis. This study provides a better understanding of fine-scale population structures within and between Swiss sheep breeds. This relevant information will help to increase the conservation activities of the local Swiss sheep breeds.

Prehistoric genomes reveal the genetic foundation and cost of horse domestication

The domestication of the horse revolutionized warfare, trade, and the exchange of people and ideas. This at least 5,500-y-long process, which ultimately transformed wild horses into the hundreds of breeds living today, is difficult to reconstruct from archeological data and modern genetics alone. We therefore sequenced two complete horse genomes, predating domestication by thousands of years, to characterize the genetic footprint of domestication. These ancient genomes reveal predomestic population structure and a significant fraction of genetic variation shared with the domestic breeds but absent from Przewalski’s horses. We find positive selection on genes involved in various aspects of locomotion, physiology, and cognition. Finally, we show that modern horse genomes contain an excess of deleterious mutations, likely representing the genetic cost of domestication.

Elucidation of the gas-phase structure of a sugar-modified DNA analogue

Antisense oligonucleotides are medical agents for the treatment of genetic diseases that are designed to interact specifically with mRNA. This interaction either induces enzymatic degradation of the targeted RNA or modifies processing pathways, e.g. by inducing alternative splicing of the pre-mRNA. The latter mechanism applies to the treatment of Duchenne muscular dystrophy with a sugar-modified DNA analogue called tricyclo-DNA (tcDNA). In tcDNA the ribose sugar-moiety is extended to a three-membered ring system, which augments the binding affinity and the selectivity of the antisense oligonucleotide for its target. The advent of chemically modified nucleic acids for antisense therapy presents a challenge to diagnostic tools, which must be able to cope with a variety of structural analogues. Mass spectrometry meets this demand for non-enzyme based sequencing methods ideally, because the technique is largely unaffected by structural modifications of the analyte. Sequence coverage of a fully modified tcDNA 15mer can be obtained in a single tandem mass spectrometric experiment. Beyond sequencing experiments, tandem mass spectrometry was applied to elucidate the gas-phase structure and stability of tcDNA:DNA and tcDNA:RNA hybrid duplexes. Most remarkable is the formation of truncated duplexes upon collision-induced dissociation of these structures. Our data suggest that the cleavage site within the duplex is directed by the modified sugar-moiety. Moreover, the formation of truncated duplexes manifests the exceptional stability of the hybrid duplexes in the gas-phase. This stability arises from the modified sugar-moiety, which locks the tcDNA single strand into a conformation that is similar to RNA in A-form duplexes. The conformational particularity of tcDNA in the gas-phase was confirmed by ion mobility-mass spectrometry experiments on tcDNA, DNA, and RNA.

Genetic and functional analyses of cell division proteins FtsZ and FtsA in Escherichia coli

In the current model for bacterial cell division, the FtsZ protein forms a ring that marks the division plane, creating a cytoskeletal framework for the subsequent action of other essential division proteins such as FtsA and ZipA. The putative protein complex ultimately generates the division septum. The essential cell division protein FtsZ is a functional and structural homolog of eukaryotic tubulin, and like tubulin, FtsZ hydrolyzes GTP and self-assembles into protein filaments in a strictly GTP-dependent manner. FtsA shares sequence similarity with members of the ATPase superfamily that include actin, but its actual function remains unknown. To test the division model and elucidate functions of the division proteins, this dissertation primarily focuses on the analysis of FtsZ and FtsA in Escherichia coli. ^ By tagging with green fluorescent protein, we first demonstrated that FtsA also exhibits a ring-like structure at the potential division site. The localization of FtsA was dependent on functional FtsZ, suggesting that FtsA is recruited to the septum by the FtsZ ring. In support of this idea, we showed that FtsA and FtsZ directly interact. Using a novel E. coli in situ assay, we found that the FtsA-FtsZ interaction appears to be species-specific, although an interspecies interaction could occur between FtsA and FtsZ proteins from two closely related organisms. In addition, mutagenesis of FtsA revealed that no single domain is solely responsible for its septal localization or interaction with FtsZ. To explore the function of FtsA, we purified FtsA protein and demonstrated that it has ATPase activity. Furthermore, purified FtsA stimulates disassembly of FtsZ polymers in a sedimentation assay but does not affect GTP hydrolysis of FtsZ. This result suggests that in the cell, FtsA may function similarly in regulating dynamic instability of the FtsZ ring during the cell division process. ^ To elucidate the structure-function relationship of FtsZ, we carried out thorough genetic and functional analyses of the mutagenized FtsZ derivatives. Our results indicate that the conserved N-terminal domain of FtsZ is necessary and sufficient for FtsZ self-assembly and localization. Moreover, we discovered a critical role for an extreme C-terminal domain of FtsZ that consists of only 12 residues. Truncated FtsZ derivatives lacking this domain, though able to polymerize and localize, are defective in ring formation in vivo as well as interaction with FtsA and ZipA. Alanine scanning mutagenesis of this region pinpointed at least five residues necessary for the function of FtsZ. Studies of protein levels and protein-protein interactions suggested that these residues may be involved in regulating protein stability and/or FtsZ-FtsA interactions. Interestingly, two of the point mutants exhibited dominant-negative phenotypes. ^ In summary, results from this thesis work have provided additional support for the division machinery model and will contribute to a better understanding of the coordinate functions of FtsA and FtsZ in the cell division process. ^

Exact moment calculations for genetic models with migration, mutation, and drift

Using properties of moment stationarity we develop exact expressions for the mean and covariance of allele frequencies at a single locus for a set of populations subject to drift, mutation, and migration. Some general results can be obtained even for arbitrary mutation and migration matrices, for example: (1) Under quite general conditions, the mean vector depends only on mutation rates, not on migration rates or the number of populations. (2) Allele frequencies covary among all pairs of populations connected by migration. As a result, the drift, mutation, migration process is not ergodic when any finite number of populations is exchanging genes. in addition, we provide closed form expressions for the mean and covariance of allele frequencies in Wright's finite-island model of migration under several simple models of mutation, and we show that the correlation in allele frequencies among populations can be very large for realistic rates of mutation unless an enormous number of populations are exchanging genes. As a result, the traditional diffusion approximation provides a poor approximation of the stationary distribution of allele frequencies among populations. Finally, we discuss some implications of our results for measures of population structure based on Wright's F-statistics.

Genetic integration of a nuclear -encoded mitochondrial gene with cardiac function

Cardiovascular disease (CVD) is the leading cause of death in the United States. One manifestation of CVD known to increase mortality is an enlarged, or hypertrophic heart. Hypertrophic cardiomyocytes adapt to increased contractile demand at the genetic level with a re-emergence of the fetal gene program and a downregulation of fatty acid oxidation genes with concomitant increased reliance on glucose-based metabolism. To understand the transcriptional regulatory pathways that implement hypertrophic directives we analyzed the upstream promoter region of the muscle specific isoform of the nuclear-encoded mitochondrial gene, carnitine palmitoyltransferase-1β (CPT-1β) in cultured rat neonatal cardiac myocytes. This enzyme catalyzes the rate-limiting step of fatty acid entry into β-oxidation and is downregulated in cardiac hypertrophy and failure, making it an attractive model for the study of hypertrophic gene regulation and metabolic adaptations. We demonstrate that the muscle-enriched transcription factors GATA-4 and SRF synergistically activate CPT-1β; moreover, DNA binding to cognate sites and intact protein structure are required. This mechanism coordinates upregulation of energy generating processes with activation of the energy consuming contractile promoter for cardiac α-actin. We hypothesized that fatty acid or glucose responsive transcription factors may also regulate CPT-1β. Oleate weakly stimulates CPT-1β activity; in contrast, the glucose responsive Upstream Stimulatory Factors (USF) dramatically depresses the CPT-1β reporter. USF regulates CPT-1β through a novel physical interaction with the cofactor PGC-1 and abrogation of MEF2A/PGC-1 synergistic stimulation. In this way, USF can inversely regulate metabolic gene programs and may play a role in the shift of metabolic substrate preference seen in hypertrophy. Failing hearts have elevated expression of the nuclear hormone receptor COUP-TF. We report that COUP-TF significantly suppresses reporter transcription independent of DNA binding and specific interactions with GATA-4, Nkx2.5 or USF. In summary, CPT-1β transcriptional regulation integrates mitochondrial gene expression with two essential cardiac functions: contraction and metabolic substrate oxidation. ^

Making genetic biodiversity measurable : a review of statistical multivariate methods to study variability at gene level

Measures of agro-ecosystems genetic variability are essential to sustain scientific-based actions and policies tending to protect the ecosystem services they provide. To build the genetic variability datum it is necessary to deal with a large number and different types of variables. Molecular marker data is highly dimensional by nature, and frequently additional types of information are obtained, as morphological and physiological traits. This way, genetic variability studies are usually associated with the measurement of several traits on each entity. Multivariate methods are aimed at finding proximities between entities characterized by multiple traits by summarizing information in few synthetic variables. In this work we discuss and illustrate several multivariate methods used for different purposes to build the datum of genetic variability. We include methods applied in studies for exploring the spatial structure of genetic variability and the association of genetic data to other sources of information. Multivariate techniques allow the pursuit of the genetic variability datum, as a unifying notion that merges concepts of type, abundance and distribution of variability at gene level.

The effects of cocaine self-administration on dendritic spine density in the rat hippocampus are dependent on genetic background

Chronic exposure to cocaine induces modifications to neurons in the brain regions involved in addiction. Hence, we evaluated cocaine-induced changes in the hippocampal CA1 field in Fischer 344 (F344) and Lewis (LEW) rats, 2 strains that have been widely used to study genetic predisposition to drug addiction, by combining intracellular Lucifer yellow injection with confocal microscopy reconstruction of labeled neurons. Specifically, we examined the effects of cocaine self-administration on the structure, size, and branching complexity of the apical dendrites of CA1 pyramidal neurons. In addition, we quantified spine density in the collaterals of the apical dendritic arbors of these neurons. We found differences between these strains in several morphological parameters. For example, CA1 apical dendrites were more branched and complex in LEW than in F344 rats, while the spine density in the collateral dendrites of the apical dendritic arbors was greater in F344 rats. Interestingly, cocaine self-administration in LEW rats augmented the spine density, an effect that was not observed in the F344 strain. These results reveal significant structural differences in CA1 pyramidal cells between these strains and indicate that cocaine self-administration has a distinct effect on neuron morphology in the hippocampus of rats with different genetic backgrounds.

Genetic characterization of chestnut (Castanea sativa Mill.) orchards and traditional nut varieties in El Bierzo, a glacial refuge and major cultivation site in Northwestern Spain

This paper presents a detailed genetic study of Castanea sativa in El Bierzo, a major nut production region with interesting features. It is located within a glacial refuge at one extreme of the distribution area (northwest Spain); it has a centenary tradition of chestnut management; and more importantly, it shows an unusual degree of genetic isolation. Seven nuclear microsatellite markers were selected to analyze the genetic variability and structure of 169 local trees grafted for nut production. We analyzed in the same manner 62 local nuts. The selected loci were highly discriminant for the genotypes studied, giving a combined probability of identity of 6.1 × 10−6. An unprecedented density of trees was sampled for this project over the entire region, and nuts were collected representing 18 cultivars marketed by local producers. Several instances of misclassification by local growers were detected. Fixation index estimates and analysis of molecular variance (AMOVA) data are supportive of an unexpectedly high level of genetic differentiation in El Bierzo, larger than that estimated in a previous study with broader geographical scope but based on limited local sampling (Pereira-Lorenzo et al., Tree Genet Genomes 6: 701–715, 2010a). Likewise, we have determined that clonality due to grafting had been previously overestimated. In line with these observations, no significant spatial structure was found using both a model-based Bayesian procedure and Mantel’s tests. Taken together, our results evidence the need for more fine-scale genetic studies if conservation strategies are to be efficiently improved.

Ecological and genetic determinants of Pepino mosaic virus emergence

Virus emergence is a complex phenomenon, which generally involves spread to a new host from a wild host, followed by adaptation to the new host. Although viruses account for the largest fraction of emerging crop pathogens, knowledge about their emergence is incomplete. We address here the question of whether Pepino mosaic virus (PepMV) emergence as a major tomato pathogen worldwide could have involved spread from wild to cultivated plant species and host adaptation. For this, we surveyed natural populations of wild tomatoes in southern Peru for PepMV infection. PepMV incidence, genetic variation, population structure, and accumulation in various hosts were analyzed. PepMV incidence in wild tomatoes was high, and a strain not yet reported in domestic tomato was characterized. This strain had a wide host range within the Solanaceae, multiplying efficiently in most assayed Solanum species and being adapted to wild tomato hosts. Conversely, PepMV isolates from tomato crops showed evidence of adaptation to domestic tomato, possibly traded against adaptation to wild tomatoes. Phylogenetic reconstructions indicated that the most probable ancestral sequence came from a wild Solanum species. A high incidence of PepMV in wild tomato relatives would favor virus spread to crops and its efficient multiplication in different Solanum species, including tomato, allowing its establishment as an epidemic pathogen. Later, adaptation to tomato, traded off against adaptation to other Solanum species, would isolate tomato populations from those in other hosts.

A genetic algorithm for slope stability analyses with concave slip surfaces using custom operators

Heuristic methods are popular tools to find critical slip surfaces in slope stability analyses. A new genetic algorithm (GA) is proposed in this work that has a standard structure but a novel encoding and generation of individuals with custom-designed operators for mutation and crossover that produce kinematically feasible slip surfaces with a high probability. In addition, new indices to assess the efficiency of operators in their search for the minimum factor of safety (FS) are proposed. The proposed GA is applied to traditional benchmark examples from the literature, as well as to a new practical example. Results show that the proposed GA is reliable, flexible and robust: it provides good minimum FS estimates that are not very sensitive to the number of nodes and that are very similar for different replications

Unearthing the roots of degradation of Quercus pyrenaica coppices: an integrative perspective from clonal structure to carbon budgets

Quercus pyrenaica es una especie rebrotadora de raíz intensa e históricamente aprovechada en monte bajo para la obtención de leñas, carbón y pastos. Debido al éxodo rural y a la aparición de nuevas fuentes energéticas, este aprovechamiento fue abandonado en la década de 1970. Desde entonces, las bajas producciones de madera y bellota y el puntisecado de los pies evidencian el generalizado estancamiento de estas masas. Uno de los mayores retos actuales de la selvicultura en el ámbito mediterráneo es encontrar usos alternativos para estos montes abandonados, siendo la conversión a monte alto una de las alternativas preferidas. Se han realizado resalveos de conversión, sin embrago, éstos se aplican sin un conocimiento integral de las causas de la degradación. En esta tesis doctoral, estudiamos un hipotético desequilibrio entre la parte radical y la parte aérea (R:S) de las cepas de rebollo como causa subyacente de su decaimiento. En una parcela experimental, aprovechada al menos desde el siglo XII, se realizaron análisis genéticos a priori para elucidar la estructura genética del rodal, y así estudiar la influencia del tamaño clonal en el funcionamiento de las cepas. Las cepas de mayor tamaño presentaron un menor crecimiento diametral de sus pies, así como mayores tasas de respiración radical, estimadas a partir de flujos internos de CO2 a través del xilema (FT) y de los flujos de CO2 del suelo. Estos resultados sugieren que el desequilibrio R:S aumenta con el tamaño clonal, dado que la eliminación periódica de órganos aéreos, al mismo tiempo que las raíces permanecen intactas, da lugar a un gran desarrollo del sistema radical que consume gran parte de los carbohidratos no estructurales (NSC) en respiración de mantenimiento, comprometiendo así el desarrollo de órganos aéreos. Se excavaron y pesaron dos cepas compuestas por cuatro y ocho pies, las cuales mostraron ratios R:S (0.5 y 1, respectivamente) superiores a los registrados en pies de origen sexual. Al igual que en otras especies rebrotadoras de raíz, se observaron altas concentraciones de NSC en las raíces (> 20% en primavera) y una gran proporción de albura en el sistema radical (52%) que alberga una notable reserva de NSC (87 kg en la cepa de mayor tamaño). En el sistema radical de dicha cepa, estimada mediante dataciones radiocarbónicas en 550 años de edad, se contaron 248 uniones radicales. La persistencia de sistemas radicales grandes, viejos, y altamente interconectados sugiere que la gran cantidad de recursos almacenados y consumidos en las raíces compensan un pobre desarrollo aéreo con una alta resiliencia vegetativa. Para un mejor entendimiento de los balances de carbono y del agotamiento de NSC en las cepas de rebollo, se midieron los flujos internos y externos de CO2 en troncos y los flujos de CO2 del suelo, y se estimó la respiración de órganos aéreos (RS) y subterráneos (RR). Estacionalmente, RS y RR reflejaron las dinámicas de flujo de savia y de crecimiento del tronco, y estuvieron determinadas principalmente por los flujos externos de CO2, dada la escasa contribución de FT a RS y RR (< 10% y < 2%, respectivamente). En una escala circadiana, la contribución de FT a RS aumentó hasta un 25% en momentos de alta transpiración. Las bajas concentraciones de CO2 en el xilema ([CO2] hasta un 0.11%) determinaron comparativamente unos bajos FT, probablemente causados por una limitada respiración del xilema y una baja resistencia a la difusión radial del CO2 impuestos por la sequía estival. Los pulsos de [CO2] observados tras las primeras lluvias de otoño apoyan esta idea. A lo largo del periodo vegetativo, el flujo medio de CO2 procedente del suelo (39 mol CO2 day-1) fue el mayor flujo respiratorio, tres y cuatro veces superior a RS (12 mol CO2 day-1) y RR (8-9 mol CO2 day-1), respectivamente. Ratios RR/RS menores que la unidad evidencian un importante peso de la respiración aérea como sumidero de carbono adicional. Finalmente, se ensayó el zanjado de raíces y el anillamiento de troncos como tratamientos selvícolas alternativos con el objetivo de aumentar las reservas de NSC en los troncos de las cepas. Los resultados preliminares desaconsejan el zanjado de raíces por el alto coste derivado posiblemente de la cicatrización de las heridas. El anillado de troncos imposibilitó el transporte de NSC a las raíces y aumentó la concentración de almidón por encima de la zona anillada, mientras que sistema radical se mantiene por los pies no anillados de la cepa. Son necesarias más mediciones y datos adicionales para comprobar el mantenimiento de esta respuesta positiva a largo plazo. Para concluir, destacamos la necesidad de estudios multidisciplinares que permitan una comprensión integral de la degradación de los rebollares ibéricos para poder aplicar a posteriori una gestión adecuada en estos montes bajos abandonados. ABSTRACT Quercus pyrenaica is a vigorous root-resprouting species intensively and historically coppiced for firewood, charcoal and woody pastures. Due to the rural exodus and the appearance of new energy sources, coppicing was abandoned towards 1970. Since then, tree overaging has resulted in stand stagnation displayed by slow stem growth, branch dieback, and scarce acorn production. The urgent need to find new alternative uses for abandoned coppices is recognized as one of the biggest challenges which currently faces Mediterranean silviculture; conversion into high forest by thinning is one of the preferred alternatives. For this aim, thinning has been broadly applied and seldom tested, although without a comprehensive understanding of the causes of stand stagnation. In this PhD study, we test the hypothesis of an imbalance between above- and below-ground organs, result of long term coppicing, as the underlying cause of Q. pyrenaica decay. In an experimental plot coppiced since at least the 12th century, genetic analyses were performed a priori to elucidate inconspicuous clonal structure of Q. pyrenaica to evaluate how clonal size affects the functioning of these multi-stemmed trees. Clonal size negatively affected diametric stem growth, whereas root respiration rates, measured by internal fluxes of CO2 through xylem (FT) and soil CO2 efflux, increased with clonal size. These results suggest root-to-shoot (R:S) imbalance intensifying with clonal size: periodic removal of aboveground organs whilst belowground organs remain undisturbed may have led to massive root systems which consume a great proportion of non-structural carbohydrates (NSC) for maintenance respiration, thus constraining aboveground performance. Furthermore, excavation of two multi-stemmed trees, composed by four and eight stems, revealed R:S ratios (0.5 and 1, respectively) greater than those reported for sexually regenerated trees. Moreover, as similarly observed in several root-resprouting species, NSC allocation to roots was favored ([NSC] > 20% in spring): a large proportion of sapwood maintained throughout the root system (52%) stored a remarkable NSC pool of 87 kg in the case of the largest clone. In this root system of the eight-stemmed tree, 248 root connections were counted and, by radiocarbon dating, its age was estimated to be 550-years-old. Persistence of massive, old and highly interconnected root systems suggests that enhanced belowground NSC storage and consumption reflects a trade-off between vegetative resilience and aboveground development. For a better understanding of tree carbon budget and the potential role of carbon starvation in Q. pyrenaica decay, internal and external stem CO2 fluxes and soil CO2 effluxes were monitored to evaluate respiratory costs above- and below-ground. On a seasonal scale, stem and root respiration (RS and RR) mirrored sap flow and stem growth dynamics. Respiration was determined to the greatest extent by external fluxes of CO2 to the atmosphere or soil, since FT accounted for a low proportion of RS and RR (< 10% and < 2%, respectively). On a diel scale, the contribution of FT to RS increased up to 25% at high transpiration rates. Comparatively low FT was determined by the low concentration of xylem CO2 registered ([CO2] as low as 0.11%), likely as a consequence of constrained xylem respiration and reduced resistance to CO2 radial diffusion imposed by summer drought. Xylem [CO2] pulses following first autumn rains support this idea. Averaged over the growing season, soil CO2 efflux was the greatest respiratory flux (39 mol CO2 day-1), three and four times greater than RS (12 mol CO2 day-1) and RR (8-9 mol CO2 day-1), respectively. Ratios of RR/RS below one evidence an additional and important weight of aboveground respiration as a tree carbon sink. Finally, root trenching and stem girdling were tested as complimentary treatments to thinning as a means to improve carbon reserves in stems of clonal trees. Preliminary results discouraged root trenching due to the high cost likely incurred for wound closure. Stem girdling successfully blocked NSC translocation downward, increasing starch concentrations above the girdled zone whilst the root system is fed by non-girdled stems within the clone. Further measurements and ancillary data are necessary to verify that this positive effect hold over time. To conclude, the need of multidisciplinary approaches for an integrative understanding on the functioning of abandoned Q pyrenaica coppices is highlighted for an appropriate management of these stands.

The underlying pathway structure of biochemical reaction networks

Bioinformatics is yielding extensive, and in some cases complete, genetic and biochemical information about individual cell types and cellular processes, providing the composition of living cells and the molecular structure of its components. These components together perform integrated cellular functions that now need to be analyzed. In particular, the functional definition of biochemical pathways and their role in the context of the whole cell is lacking. In this study, we show how the mass balance constraints that govern the function of biochemical reaction networks lead to the translation of this problem into the realm of linear algebra. The functional capabilities of biochemical reaction networks, and thus the choices that cells can make, are reflected in the null space of their stoichiometric matrix. The null space is spanned by a finite number of basis vectors. We present an algorithm for the synthesis of a set of basis vectors for spanning the null space of the stoichiometric matrix, in which these basis vectors represent the underlying biochemical pathways that are fundamental to the corresponding biochemical reaction network. In other words, all possible flux distributions achievable by a defined set of biochemical reactions are represented by a linear combination of these basis pathways. These basis pathways thus represent the underlying pathway structure of the defined biochemical reaction network. This development is significant from a fundamental and conceptual standpoint because it yields a holistic definition of biochemical pathways in contrast to definitions that have arisen from the historical development of our knowledge about biochemical processes. Additionally, this new conceptual framework will be important in defining, characterizing, and studying biochemical pathways from the rapidly growing information on cellular function.