Abordagem Bayesiana da curva de lactação de cabras Saanen de primeira e segunda ordem de parto

O objetivo deste trabalho foi utilizar o método Bayesiano no ajuste do modelo de Wood a dados de produção de leite de cabras da raça Saanen. Dois grupos de animais da primeira e segunda lactação foram considerados. Amostras das distribuições marginais a posteriori dos parâmetros do modelo de Wood e das funções de produção derivadas desses parâmetros - pico de produção, tempo do pico de produção, persistência e produção total de leite - foram obtidas pelo algoritmo Gibbs Sampler. As inferências foram feitas em cada população e os resultados mostraram diferenças na taxa de decréscimo da produção após o pico e na persistência, indicando maior produção nos animais de segunda lactação. Realizou-se um estudo de simulação de dados para avaliar o método Bayesiano sob diferentes estruturas de matrizes de covariâncias dos parâmetros. Os resultados desse estudo indicam que o método é eficiente no estudo das curvas de lactação quando a matriz de covariância apresenta alta correlação dos parâmetros.

Herdabilidade da habilidade de permanência no rebanho em fêmeas de bovinos da raça Canchim

O objetivo deste trabalho foi avaliar a possibilidade de a característica habilidade de permanência no rebanho de matrizes ser utilizada como critério de seleção na raça Canchim. A habilidade de permanência no rebanho foi definida como a probabilidade de uma matriz ter três partos ou mais até os 76 meses de idade, desde que ela tenha tido a oportunidade de se reproduzir pelo menos uma vez. Codificações binárias, com 0 para indicar fracasso e 1 para indicar sucesso, foram feitas para cada fêmea. Uma codificação alternativa, com quatro categorias, considerou os valores 0, 1, 2 e 3, respectivamente, para vacas com menos de três, com três, quatro ou cinco partos até os 76 meses de idade. As análises foram processadas pela aplicação de um modelo unicaracterístico de touro avô-materno. Utilizou-se a amostragem de Gibbs para estimar os componentes de variância e herdabilidades. As estimativas de herdabilidade foram de 0,07, para a característica definida pelo modelo binário, e de 0,08, para a característica definida pelo modelo com quatro categorias. Os resultados indicam pequenas diferenças na estimativa de herdabilidade, pelos modelos binário e com quatro categorias. Sugerem também que as estimativas de herdabilidade, para a característica habilidade de permanência, na raça Canchim, apresenta baixo potencial de resposta à seleção.

Zinc adsorption in highly weathered soils

The objective of this work was to assess the effects of pH and ionic strength upon zinc adsorption, in three highly weathered variable charge soils. Adsorption isotherms were elaborated from batch adsorption experiments, with increasing Zn concentrations (0-80 mg L-1), and adsorption envelopes were constructed through soil samples reactions with 0.01, 0.1 and 1 mol L-1 Ca(NO3)2 solutions containing 5 mg L-1 of Zn, with an increasing pH value from 3 to 8. Driving force of reaction was quantified by Gibbs free energy and separation factor. Isotherms were C-, H- and L-type and experimental results were fitted to nonlinear Langmuir model. Maximum adsorption ranged from 59-810 mg kg-1, and Zn affinity was greater in subsoil (0.13-0.81 L kg-1) than in the topsoil samples (0.01-0.34 L kg-1). Zinc adsorption was favorable and spontaneous, and showed sharply increase (20-90%) in the 4-6 pH range. No effect of ionic strength was observed at pH values below 5, because specific adsorption mechanisms predominated in the 3-5 pH range. Above pH 5, and in subsoil samples, Zn was adsorbed by electrostatic mechanisms, since ionic strength effect was observed. Despite depth and ionic strength effects, Zn adsorption depends mainly on the pH.

Modeling and control of electromechanical systems

The material presented in the these notes covers the sessions Modelling of electromechanical systems, Passive control theory I and Passive control theory II of the II EURON/GEOPLEX Summer School on Modelling and Control of Complex Dynamical Systems.We start with a general description of what an electromechanical system is from a network modelling point of view. Next, a general formulation in terms of PHDS is introduced, and some of the previous electromechanical systems are rewritten in this formalism. Power converters, which are variable structure systems (VSS), can also be given a PHDS form.We conclude the modelling part of these lectures with a rather complex example, showing the interconnection of subsystems from several domains, namely an arrangement to temporally store the surplus energy in a section of a metropolitan transportation system based on dc motor vehicles, using either arrays of supercapacitors or an electric poweredflywheel. The second part of the lectures addresses control of PHD systems. We first present the idea of control as power connection of a plant and a controller. Next we discuss how to circumvent this obstacle and present the basic ideas of Interconnection and Damping Assignment (IDA) passivity-based control of PHD systems.

Análise genética de escores de avaliação visual de bovinos com modelos bayesianos de limiar e linear

O objetivo deste trabalho foi comparar as estimativas de parâmetros genéticos obtidas em análises bayesianas uni-característica e bi-característica, em modelo animal linear e de limiar, considerando-se as características categóricas morfológicas de bovinos da raça Nelore. Os dados de musculosidade, estrutura física e conformação foram obtidos entre 2000 e 2005, em 3.864 animais de 13 fazendas participantes do Programa Nelore Brasil. Foram realizadas análises bayesianas uni e bi-características, em modelos de limiar e linear. De modo geral, os modelos de limiar e linear foram eficientes na estimação dos parâmetros genéticos para escores visuais em análises bayesianas uni-características. Nas análises bi-características, observou-se que: com utilização de dados contínuos e categóricos, o modelo de limiar proporcionou estimativas de correlação genética de maior magnitude do que aquelas do modelo linear; e com o uso de dados categóricos, as estimativas de herdabilidade foram semelhantes. A vantagem do modelo linear foi o menor tempo gasto no processamento das análises. Na avaliação genética de animais para escores visuais, o uso do modelo de limiar ou linear não influenciou a classificação dos animais, quanto aos valores genéticos preditos, o que indica que ambos os modelos podem ser utilizados em programas de melhoramento genético.

Bayesian genetic parameters for body weight and survival of Nile tilapia farmed in Brazil

The objective of this study was to estimate genetic parameters for survival and weight of Nile tilapia (Oreochromis niloticus), farmed in cages and ponds in Brazil, and to predict genetic gain under different scenarios. Survival was recorded as a binary response (dead or alive), during harvest time in the 2008 grow-out period. Genetic parameters were estimated using a Bayesian mixed linear-threshold animal model via Gibbs sampling. The breeding population consisted of 2,912 individual fish, which were analyzed together with the pedigree of 5,394 fish. The heritabilities estimates, with 95% posterior credible intervals, for tagging weight, harvest weight and survival were 0.17 (0.09-0.27), 0.21 (0.12-0.32) and 0.32 (0.22-0.44), respectively. Credible intervals show a 95% probability that the true genetic correlations were in a favourable direction. The selection for weight has a positive impact on survival. Estimated genetic gain was high when selecting for harvest weight (5.07%), and indirect gain for tagging weight (2.17%) and survival (2.03%) were also considerable.

Modeling morphogen gradient formation from arbitrary realistically shaped sources.

Much of the analytical modeling of morphogen profiles is based on simplistic scenarios, where the source is abstracted to be point-like and fixed in time, and where only the steady state solution of the morphogen gradient in one dimension is considered. Here we develop a general formalism allowing to model diffusive gradient formation from an arbitrary source. This mathematical framework, based on the Green's function method, applies to various diffusion problems. In this paper, we illustrate our theory with the explicit example of the Bicoid gradient establishment in Drosophila embryos. The gradient formation arises by protein translation from a mRNA distribution followed by morphogen diffusion with linear degradation. We investigate quantitatively the influence of spatial extension and time evolution of the source on the morphogen profile. For different biologically meaningful cases, we obtain explicit analytical expressions for both the steady state and time-dependent 1D problems. We show that extended sources, whether of finite size or normally distributed, give rise to more realistic gradients compared to a single point-source at the origin. Furthermore, the steady state solutions are fully compatible with a decreasing exponential behavior of the profile. We also consider the case of a dynamic source (e.g. bicoid mRNA diffusion) for which a protein profile similar to the ones obtained from static sources can be achieved.

The genome sequence of taurine cattle: a window to ruminant biology and evolution.

To understand the biology and evolution of ruminants, the cattle genome was sequenced to about sevenfold coverage. The cattle genome contains a minimum of 22,000 genes, with a core set of 14,345 orthologs shared among seven mammalian species of which 1217 are absent or undetected in noneutherian (marsupial or monotreme) genomes. Cattle-specific evolutionary breakpoint regions in chromosomes have a higher density of segmental duplications, enrichment of repetitive elements, and species-specific variations in genes associated with lactation and immune responsiveness. Genes involved in metabolism are generally highly conserved, although five metabolic genes are deleted or extensively diverged from their human orthologs. The cattle genome sequence thus provides a resource for understanding mammalian evolution and accelerating livestock genetic improvement for milk and meat production.

Abordagem bayesiana multivariada para características de crescimento, fertilidade e escores visuais de rebanhos da raça Brangus

O objetivo deste trabalho foi estimar parâmetros genéticos e tendências genéticas e fenotípicas de uma população da raça Brangus. As características peso, circunferência escrotal e escores visuais de conformação, precocidade, musculatura e umbigo, padronizadas para 550 dias de idade, foram avaliadas a partir de 6.789 registros de animais nascidos de 288 touros e 5.949 vacas, entre 1991 e 2001, em 49 fazendas localizadas nas regiões Centro‑Oeste, Sudeste e Sul do Brasil. Para a estimação dos parâmetros, das correlações e das tendências genéticas, foi adotado o modelo animal linear‑limiar hexacaracterística. As estimativas de herdabilidade direta foram de 0,39 e 0,27, para peso e circunferência escrotal, respectivamente, e de 0,22, 0,20, 0,23 e 0,33 para conformação, precocidade, musculatura e umbigo, o que indica considerável variação genética aditiva e que é possível obter ganho genético por meio da seleção. As correlações genéticas entre peso e circunferência escrotal com os escores de conformação, precocidade e musculatura mostram a possibilidade de resposta correlacionada. As tendências genéticas estimadas indicam grande contribuição de fontes de variação não genéticas para todas as características no período estudado, e apontam a necessidade de melhoria das condições ambientais, para que os animais expressem todo seu potencial genético.

Associação genética de escores de conformação com características de crescimento em bovinos da raça Angus

O objetivo deste trabalho foi determinar a associação genética entre escores visuais de conformação e as características de ganho de peso médio diário e de velocidade de crescimento em bovinos da raça Angus à desmama e ao sobreano. Os componentes de covariância foram estimados por modelo animal de análise tetracaracterística, com uso do método de inferência bayesiana, tendo-se assumido o modelo linear para: ganho de peso médio diário do nascimento à desmama (GMD) e da desmama ao sobreano (GMS); e velocidade de ganho de peso do nascimento à desmama (VD) e da desmama ao sobreano (VS). Um modelo não linear (de limiar) foi utilizado para os escores de conformação à desmama (CD) e ao sobreano (CS). As médias a posteriori, para a herdabilidade direta, foram: 0,12±0,023 (CD), 0,15±0,020 (GMD), 0,15±0,024 (VD), 0,17±0,020 (CS), 0,17±0,023(GMS), e 0,17±0,023 (VS). A correlação genética variou de -0,09±0,11 a 0,60±0,06, entre os escores CD e CS e as características de ganho médio diário de peso e velocidade de ganho de peso. A correlação entre CD e CS foi 0,52±0,089. A seleção direta para escores visuais de conformação, ganho médio diário e velocidade de ganho responde de forma lenta à seleção, tanto à desmama como ao sobreano.

Implicit methods for qualitative modeling of gene regulatory networks.

Advancements in high-throughput technologies to measure increasingly complex biological phenomena at the genomic level are rapidly changing the face of biological research from the single-gene single-protein experimental approach to studying the behavior of a gene in the context of the entire genome (and proteome). This shift in research methodologies has resulted in a new field of network biology that deals with modeling cellular behavior in terms of network structures such as signaling pathways and gene regulatory networks. In these networks, different biological entities such as genes, proteins, and metabolites interact with each other, giving rise to a dynamical system. Even though there exists a mature field of dynamical systems theory to model such network structures, some technical challenges are unique to biology such as the inability to measure precise kinetic information on gene-gene or gene-protein interactions and the need to model increasingly large networks comprising thousands of nodes. These challenges have renewed interest in developing new computational techniques for modeling complex biological systems. This chapter presents a modeling framework based on Boolean algebra and finite-state machines that are reminiscent of the approach used for digital circuit synthesis and simulation in the field of very-large-scale integration (VLSI). The proposed formalism enables a common mathematical framework to develop computational techniques for modeling different aspects of the regulatory networks such as steady-state behavior, stochasticity, and gene perturbation experiments.

An Improved Hybrid Model for the Generic Hoist Scheduling Problem

Peer-reviewed

Multi-scale systems analysis of plant development : beyond cellular signaling and tissue morphodynamics

Les plantes sont essentielles pour les sociétés humaines. Notre alimentation quotidienne, les matériaux de constructions et les sources énergétiques dérivent de la biomasse végétale. En revanche, la compréhension des multiples aspects développementaux des plantes est encore peu exploitée et représente un sujet de recherche majeur pour la science. L'émergence des technologies à haut débit pour le séquençage de génome à grande échelle ou l'imagerie de haute résolution permet à présent de produire des quantités énormes d'information. L'analyse informatique est une façon d'intégrer ces données et de réduire la complexité apparente vers une échelle d'abstraction appropriée, dont la finalité est de fournir des perspectives de recherches ciblées. Ceci représente la raison première de cette thèse. En d'autres termes, nous appliquons des méthodes descriptives et prédictives combinées à des simulations numériques afin d'apporter des solutions originales à des problèmes relatifs à la morphogénèse à l'échelle de la cellule et de l'organe. Nous nous sommes fixés parmi les objectifs principaux de cette thèse d'élucider de quelle manière l'interaction croisée des phytohormones auxine et brassinosteroïdes (BRs) détermine la croissance de la cellule dans la racine du méristème apical d'Arabidopsis thaliana, l'organisme modèle de référence pour les études moléculaires en plantes. Pour reconstruire le réseau de signalement cellulaire, nous avons extrait de la littérature les informations pertinentes concernant les relations entre les protéines impliquées dans la transduction des signaux hormonaux. Le réseau a ensuite été modélisé en utilisant un formalisme logique et qualitatif pour pallier l'absence de données quantitatives. Tout d'abord, Les résultats ont permis de confirmer que l'auxine et les BRs agissent en synergie pour contrôler la croissance de la cellule, puis, d'expliquer des observations phénotypiques paradoxales et au final, de mettre à jour une interaction clef entre deux protéines dans la maintenance du méristème de la racine. Une étude ultérieure chez la plante modèle Brachypodium dystachion (Brachypo- dium) a révélé l'ajustement du réseau d'interaction croisée entre auxine et éthylène par rapport à Arabidopsis. Chez ce dernier, interférer avec la biosynthèse de l'auxine mène à la formation d'une racine courte. Néanmoins, nous avons isolé chez Brachypodium un mutant hypomorphique dans la biosynthèse de l'auxine qui affiche une racine plus longue. Nous avons alors conduit une analyse morphométrique qui a confirmé que des cellules plus anisotropique (plus fines et longues) sont à l'origine de ce phénotype racinaire. Des analyses plus approfondies ont démontré que la différence phénotypique entre Brachypodium et Arabidopsis s'explique par une inversion de la fonction régulatrice dans la relation entre le réseau de signalisation par l'éthylène et la biosynthèse de l'auxine. L'analyse morphométrique utilisée dans l'étude précédente exploite le pipeline de traitement d'image de notre méthode d'histologie quantitative. Pendant la croissance secondaire, la symétrie bilatérale de l'hypocotyle est remplacée par une symétrie radiale et une organisation concentrique des tissus constitutifs. Ces tissus sont initialement composés d'une douzaine de cellules mais peuvent aisément atteindre des dizaines de milliers dans les derniers stades du développement. Cette échelle dépasse largement le seuil d'investigation par les moyens dits 'traditionnels' comme l'imagerie directe de tissus en profondeur. L'étude de ce système pendant cette phase de développement ne peut se faire qu'en réalisant des coupes fines de l'organe, ce qui empêche une compréhension des phénomènes cellulaires dynamiques sous-jacents. Nous y avons remédié en proposant une stratégie originale nommée, histologie quantitative. De fait, nous avons extrait l'information contenue dans des images de très haute résolution de sections transverses d'hypocotyles en utilisant un pipeline d'analyse et de segmentation d'image à grande échelle. Nous l'avons ensuite combiné avec un algorithme de reconnaissance automatique des cellules. Cet outil nous a permis de réaliser une description quantitative de la progression de la croissance secondaire révélant des schémas développementales non-apparents avec une inspection visuelle classique. La formation de pôle de phloèmes en structure répétée et espacée entre eux d'une longueur constante illustre les bénéfices de notre approche. Par ailleurs, l'exploitation approfondie de ces résultats a montré un changement de croissance anisotropique des cellules du cambium et du phloème qui semble en phase avec l'expansion du xylème. Combinant des outils génétiques et de la modélisation biomécanique, nous avons démontré que seule la croissance plus rapide des tissus internes peut produire une réorientation de l'axe de croissance anisotropique des tissus périphériques. Cette prédiction a été confirmée par le calcul du ratio des taux de croissance du xylème et du phloème au cours de développement secondaire ; des ratios élevés sont effectivement observés et concomitant à l'établissement progressif et tangentiel du cambium. Ces résultats suggèrent un mécanisme d'auto-organisation établi par un gradient de division méristématique qui génèrent une distribution de contraintes mécaniques. Ceci réoriente la croissance anisotropique des tissus périphériques pour supporter la croissance secondaire. - Plants are essential for human society, because our daily food, construction materials and sustainable energy are derived from plant biomass. Yet, despite this importance, the multiple developmental aspects of plants are still poorly understood and represent a major challenge for science. With the emergence of high throughput devices for genome sequencing and high-resolution imaging, data has never been so easy to collect, generating huge amounts of information. Computational analysis is one way to integrate those data and to decrease the apparent complexity towards an appropriate scale of abstraction with the aim to eventually provide new answers and direct further research perspectives. This is the motivation behind this thesis work, i.e. the application of descriptive and predictive analytics combined with computational modeling to answer problems that revolve around morphogenesis at the subcellular and organ scale. One of the goals of this thesis is to elucidate how the auxin-brassinosteroid phytohormone interaction determines the cell growth in the root apical meristem of Arabidopsis thaliana (Arabidopsis), the plant model of reference for molecular studies. The pertinent information about signaling protein relationships was obtained through the literature to reconstruct the entire hormonal crosstalk. Due to a lack of quantitative information, we employed a qualitative modeling formalism. This work permitted to confirm the synergistic effect of the hormonal crosstalk on cell elongation, to explain some of our paradoxical mutant phenotypes and to predict a novel interaction between the BREVIS RADIX (BRX) protein and the transcription factor MONOPTEROS (MP),which turned out to be critical for the maintenance of the root meristem. On the same subcellular scale, another study in the monocot model Brachypodium dystachion (Brachypodium) revealed an alternative wiring of auxin-ethylene crosstalk as compared to Arabidopsis. In the latter, increasing interference with auxin biosynthesis results in progressively shorter roots. By contrast, a hypomorphic Brachypodium mutant isolated in this study in an enzyme of the auxin biosynthesis pathway displayed a dramatically longer seminal root. Our morphometric analysis confirmed that more anisotropic cells (thinner and longer) are principally responsible for the mutant root phenotype. Further characterization pointed towards an inverted regulatory logic in the relation between ethylene signaling and auxin biosynthesis in Brachypodium as compared to Arabidopsis, which explains the phenotypic discrepancy. Finally, the morphometric analysis of hypocotyl secondary growth that we applied in this study was performed with the image-processing pipeline of our quantitative histology method. During its secondary growth, the hypocotyl reorganizes its primary bilateral symmetry to a radial symmetry of highly specialized tissues comprising several thousand cells, starting with a few dozens. However, such a scale only permits observations in thin cross-sections, severely hampering a comprehensive analysis of the morphodynamics involved. Our quantitative histology strategy overcomes this limitation. We acquired hypocotyl cross-sections from tiled high-resolution images and extracted their information content using custom high-throughput image processing and segmentation. Coupled with an automated cell type recognition algorithm, it allows precise quantitative characterization of vascular development and reveals developmental patterns that were not evident from visual inspection, for example the steady interspace distance of the phloem poles. Further analyses indicated a change in growth anisotropy of cambial and phloem cells, which appeared in phase with the expansion of xylem. Combining genetic tools and computational modeling, we showed that the reorientation of growth anisotropy axis of peripheral tissue layers only occurs when the growth rate of central tissue is higher than the peripheral one. This was confirmed by the calculation of the ratio of the growth rate xylem to phloem throughout secondary growth. High ratios are indeed observed and concomitant with the homogenization of cambium anisotropy. These results suggest a self-organization mechanism, promoted by a gradient of division in the cambium that generates a pattern of mechanical constraints. This, in turn, reorients the growth anisotropy of peripheral tissues to sustain the secondary growth.

Energy-based modelling and simulation of the interconnection of a back-to-back converter and a doubly-fed induction machine

This paper describes the port interconnection of two subsystems: a power electronics subsystem (a back-to-back AC/CA converter (B2B), coupled to a phase of the power grid), and an electromechanical subsystem (a doubly-fed induction machine (DFIM). The B2B is a variable structure system (VSS), due to presence of control-actuated switches: however, from a modelling simulation, as well as a control-design, point of view, it is sensible to consider modulated transformers (MTF in the bond graph language) instead of the pairs of complementary switches. The port-Hamiltonian models of both subsystems are presented and, using a power-preserving interconnection, the Hamiltonian description of the whole system is obtained; detailed bond graphs of all subsystems and the complete system are also provided. Using passivity-based controllers computed in the Hamiltonian formalism for both subsystems, the whole model is simulated; simulations are run to rest the correctness and efficiency of the Hamiltonian network modelling approach used in this work.

Parâmetros genéticos do peso desde o nascimento até 730 dias de idade na raça Indubrasil

O objetivo deste trabalho foi determinar parâmetros genéticos do peso desde o nascimento até 730 dias de idade na raça Indubrasil. Registros de peso (82.481) de 20.890 animais do nascimento aos 730 dias de idade foram utilizados. Os parâmetros genéticos foram determinados por meio de regressão aleatória. Amostras da distribuição a posteriori dos parâmetros de interesse foram obtidas com o amostrador de Gibbs. Polinômios de Legendre ou segmentados foram utilizados para modelar a trajetória média de crescimento. Os efeitos genético aditivo direto e de ambiente permanente direto foram modelados por polinômios de Legendre. Um polinômio segmentado linear-linear, com nó aos 251 dias de idade, ajustou a trajetória média. Polinômios de Legendre quadráticos e quínticos ajustaram, respectivamente, os efeitos genético aditivo e de ambiente permanente. Três classes de idade foram suficientes para modelar a heterocedasticidade residual. Na maior parte do intervalo considerado, as médias a posteriori da herdabilidade foram crescentes, com valores entre 0,10 e 0,45; a proporção da variância de ambiente permanente em relação à variância fenotípica foi constante, em torno de 0,48; e as correlações genéticas dos pesos em diferentes idades foram altas, acima de 0,60. Há variabilidade genética quanto ao peso na raça Indubrasil, e a seleção pode aumentar a média desta característica.