Modeling associations between genetic markers using Bayesian networks

Motivation: Understanding the patterns of association between polymorphisms at different loci in a population ( linkage disequilibrium, LD) is of fundamental importance in various genetic studies. Many coefficients were proposed for measuring the degree of LD, but they provide only a static view of the current LD structure. Generative models (GMs) were proposed to go beyond these measures, giving not only a description of the actual LD structure but also a tool to help understanding the process that generated such structure. GMs based in coalescent theory have been the most appealing because they link LD to evolutionary factors. Nevertheless, the inference and parameter estimation of such models is still computationally challenging. Results: We present a more practical method to build GM that describe LD. The method is based on learning weighted Bayesian network structures from haplotype data, extracting equivalence structure classes and using them to model LD. The results obtained in public data from the HapMap database showed that the method is a promising tool for modeling LD. The associations represented by the learned models are correlated with the traditional measure of LD D`. The method was able to represent LD blocks found by standard tools. The granularity of the association blocks and the readability of the models can be controlled in the method. The results suggest that the causality information gained by our method can be useful to tell about the conservability of the genetic markers and to guide the selection of subset of representative markers.

A Study on Unconditionally Stable FDTD Methods for the Modeling of Metamaterials

We assess the performance of three unconditionally stable finite-difference time-domain (FDTD) methods for the modeling of doubly dispersive metamaterials: 1) locally one-dimensional FDTD; 2) locally one-dimensional FDTD with Strang splitting; and (3) alternating direction implicit FDTD. We use both double-negative media and zero-index media as benchmarks.

Modeling and analysis of laminate composite plates with embedded active-passive piezoelectric networks

The objective of this work is to present the finite element modeling of laminate composite plates with embedded piezoelectric patches or layers that are then connected to active-passive resonant shunt circuits, composed of resistance, inductance and voltage source. Applications to passive vibration control and active control authority enhancement are also presented and discussed. The finite element model is based on an equivalent single layer theory combined with a third-order shear deformation theory. A stress-voltage electromechanical model is considered for the piezoelectric materials fully coupled to the electrical circuits. To this end, the electrical circuit equations are also included in the variational formulation. Hence, conservation of charge and full electromechanical coupling are guaranteed. The formulation results in a coupled finite element model with mechanical (displacements) and electrical (charges at electrodes) degrees of freedom. For a Graphite-Epoxy (Carbon-Fibre Reinforced) laminate composite plate, a parametric analysis is performed to evaluate optimal locations along the plate plane (xy) and thickness (z) that maximize the effective modal electromechanical coupling coefficient. Then, the passive vibration control performance is evaluated for a network of optimally located shunted piezoelectric patches embedded in the plate, through the design of resistance and inductance values of each circuit, to reduce the vibration amplitude of the first four vibration modes. A vibration amplitude reduction of at least 10 dB for all vibration modes was observed. Then, an analysis of the control authority enhancement due to the resonant shunt circuit, when the piezoelectric patches are used as actuators, is performed. It is shown that the control authority can indeed be improved near a selected resonance even with multiple pairs of piezoelectric patches and active-passive circuits acting simultaneously. (C) 2010 Elsevier Ltd. All rights reserved.

Hybrid modeling of convective laminar flow in a permeable tube associated with the cross-flow process

The confined flows in tubes with permeable surfaces arc associated to tangential filtration processes (microfiltration or ultrafiltration). The complexity of the phenomena do not allow for the development of exact analytical solutions, however, approximate solutions are of great interest for the calculation of the transmembrane outflow and estimate of the concentration, polarization phenomenon. In the present work, the generalized integral transform technique (GITT) was employed in solving the laminar and permanent flow in permeable tubes of Newtonian and incompressible fluid. The mathematical formulation employed the parabolic differential equation of chemical species conservation (convective-diffusive equation). The velocity profiles for the entrance region flow, which are found in the connective terms of the equation, were assessed by solutions obtained from literature. The velocity at the permeable wall was considered uniform, with the concentration at the tube wall regarded as variable with an axial position. A computational methodology using global error control was applied to determine the concentration in the wall and concentration boundary layer thickness. The results obtained for the local transmembrane flux and the concentration boundary layer thickness were compared against others in literature. (C) 2007 Elsevier B.V. All rights reserved.

A state-space modeling approach for active structural acoustic control

The demands for improvement in sound quality and reduction of noise generated by vehicles are constantly increasing, as well as the penalties for space and weight of the control solutions. A promising approach to cope with this challenge is the use of active structural-acoustic control. Usually, the low frequency noise is transmitted into the vehicle`s cabin through structural paths, which raises the necessity of dealing with vibro-acoustic models. This kind of models should allow the inclusion of sensors and actuators models, if accurate performance indexes are to be accessed. The challenge thus resides in deriving reasonable sized models that integrate structural, acoustic, electrical components and the controller algorithm. The advantages of adequate active control simulation strategies relies on the cost and time reduction in the development phase. Therefore, the aim of this paper is to present a methodology for simulating vibro-acoustic systems including this coupled model in a closed loop control simulation framework that also takes into account the interaction between the system and the control sensors/actuators. It is shown that neglecting the sensor/actuator dynamics can lead to inaccurate performance predictions.

Parametrical study of masonry walls subjected to in-plane loading through numerical modeling

This paper deals with the numerical assessment of the influence of parameters such as pre-compression level, aspect ratio, vertical and horizontal reinforcement ratios and boundary conditions on the lateral strength of masonry walls under in-plane loading. The numerical study is performed through the software DIANA (R) based on the Finite Element Method. The validation of the numerical model is carried out from a database of available experimental results on masonry walls tested under cyclic lateral loading. Numerical results revealed that boundary conditions play a central role on the lateral behavior of masonry walls under in-plane loading and determine the influence of level of pre-compression as well as the reinforcement ratio on the wall strength. The lateral capacity of walls decreases with the increase of aspect ratio and with the decrease of pre-compression. Vertical steel bars appear to have almost no influence in the shear strength of masonry walls and horizontal reinforcement only increases the lateral strength of masonry walls if the shear response of the walls is determinant for failure, which is directly related to the boundary conditions. (C) 2011 Elsevier Ltd. All rights reserved.

Performance and molecular evaluation of an anaerobic system with suspended biomass for treating wastewater with high fat content after enzymatic hydrolysis

The effect of a lipase-rich fungal enzymatic preparation, produced by a Penicillium sp. during solid-state fermentation, was evaluated in an anaerobic digester treating dairy wastewater with 1200 mg of oil and grease/L The oil and grease hydrolysis step was carried out with 0.1% (w/v) of solid enzymatic preparation at 30 degrees C for 24 h, and resulted in a final free acid concentration eight times higher than the initial value. The digester operated in sequential batches of 48 h at 30 degrees C for 245 days, and had high chemical oxygen demand (COD) removal efficiencies (around 90%) when fed with pre-hydrolyzed wastewater. However, when the pre-hydrolysis step was removed, the anaerobic digester performed poorly (with an average COD removal of 32%), as the oil and grease accumulated in the biomass and effluent oil and grease concentration increased throughout the operational period. PCR-DGGE analysis of the Bacteria and Archaea domains revealed remarkable differences in the microbial profiles in trials conducted with and without the pre-hydrolysis step, indicating that differences observed in overall parameters were intrinsically related to the microbial diversity of the anaerobic sludge. (C) 2009 Elsevier Ltd. All rights reserved.

Morphological and molecular characterization of cyanobacteria from a Brazilian facultative wastewater stabilization pond and evaluation of microcystin production

The cyanobacterial population in the Cajati waste stabilization pond system (WSP) from Sao Paulo State, Brazil was assessed by cell isolation and direct microscope counting techniques. Ten strains, belonging to five genera (Synechococcus, Merismopedia, Leptolyngbya, Limnothrix, and Nostoc), were isolated and identified by morphological and molecular analyses. Morphological identification of the isolated strains was congruent with their phylogenetic analyses based on 16S rDNA gene sequences. Six cyanobacterial genera (Synechocystis, Aphanocapsa, Merismopedia, Lyngbya, Phormidium, and Pseudanabaena) were identified by direct microscope inspection. Both techniques were complementary, since, of the six genera identified by direct microscopic inspection, only Merismopedia was isolated, and the four other isolated genera were not detected by direct inspection. Direct microscope counting of preserved cells showed that cyanobacteria were the dominant members (> 90%) of the phytoplankton community during both periods evaluated (summer and autumn). ELISA tests specific for hepatotoxicmicrocystins gave positive results for six strains (Synechococcus CENA108, Merismopedia CENA106, Leptolyngbya CENA103, Leptolyngbya CENA112, Limnothrix CENA109, and Limnothrix CENA110), and for wastewater samples collected from raw influent (3.70 mu g microcystins/l) and treated effluent (3.74 mu g microcystins/l) in summer. Our findings indicate that toxic cyanobacteria in WSP systems are of concern, since the treated effluent containing cyanotoxins will be discharged into rivers, irrigation channels, estuaries, or reservoirs, and can affect human and animal health.

Heterogeneous modeling of an anaerobic sequencing batch biofilm reactor (ASBBR)

The objective of this study was to estimate the first-order intrinsic kinetic constant (k(1)) and the liquid-phase mass transfer coefficient (k(c)) in a bench-scale anaerobic sequencing batch biofilm reactor (ASBBR) fed with glucose. A dynamic heterogeneous mathematical model, considering two phases (liquid and solid), was developed through mass balances in the liquid and solid phases. The model was adjusted to experimental data obtained from the ASBBR applied for the treatment of glucose-based synthetic wastewater with approximately 500 mg L-1 of glucose, operating in 8 h batch cycles, at 30 degrees C and 300 rpm. The values of the parameters obtained were 0.8911 min(-1) for k(1) and 0.7644 cm min(-1) for kc. The model was validated utilizing the estimated parameters with data obtained from the ASBBR operating in 3 h batch cycles, with a good representation of the experimental behavior. The solid-phase mass transfer flux was found to be the limiting step of the overall glucose conversion rate.

Solution of a heterogeneous modeling of a horizontal-flow anaerobic immobilized biomass (HAIB) reactor by the sequencing method

A modeling study was completed to develop a methodology that combines the sequencing and finite difference methods for the simulation of a heterogeneous model of a tubular reactor applied in the treatment of wastewater. The system included a liquid phase (convection diffusion transport) and a solid phase (diffusion reaction) that was obtained by completing a mass balance in the reactor and in the particle, respectively. The model was solved using a pilot-scale horizontal-flow anaerobic immobilized biomass (HAIB) reactor to treat domestic sewage, with the concentration results compared with the experimental data. A comparison of the behavior of the liquid phase concentration profile and the experimental results indicated that both the numerical methods offer a good description of the behavior of the concentration along the reactor. The advantage of the sequencing method over the finite difference method is that it is easier to apply and requires less computational time to model the dynamic simulation of outlet response of HAIB.

Ozone oxidation of pulp and paper wastewater and its impact on molecular weight distribution of organic matter

The impact of ozone oxidation on removing high molecular weight (HMW) organics in order to improve the biodegradability of alkaline bleach plant effluent was investigated using a semi-batch reactor under different initial pH (12 and 7). After the ozonation process, the ratio of BOD5/COD increased from 0.07 to 0.16 and 0.22 for initial pH 12 and 7, respectively. Also, the effluent color decreased by 48% and 61% at initial pH 12 and pH 7, respectively. These changes were primarily driven by reductions of the HMW fractions of the effluent during ozonation.

Modeling and Analysis of Piezoelectric Energy Harvesting From Aeroelastic Vibrations Using the Doublet-Lattice Method

Multifunctional structures are pointed out as an important technology for the design of aircraft with volume, mass, and energy source limitations such as unmanned air vehicles (UAVs) and micro air vehicles (MAVs). In addition to its primary function of bearing aerodynamic loads, the wing/spar structure of an UAV or a MAV with embedded piezoceramics can provide an extra electrical energy source based on the concept of vibration energy harvesting to power small and wireless electronic components. Aeroelastic vibrations of a lifting surface can be converted into electricity using piezoelectric transduction. In this paper, frequency-domain piezoaeroelastic modeling and analysis of a canti-levered platelike wing with embedded piezoceramics is presented for energy harvesting. The electromechanical finite-element plate model is based on the thin-plate (Kirchhoff) assumptions while the unsteady aerodynamic model uses the doublet-lattice method. The electromechanical and aerodynamic models are combined to obtain the piezoaeroelastic equations, which are solved using a p-k scheme that accounts for the electromechanical coupling. The evolution of the aerodynamic damping and the frequency of each mode are obtained with changing airflow speed for a given electrical circuit. Expressions for piezoaeroelastically coupled frequency response functions (voltage, current, and electrical power as well the vibratory motion) are also defined by combining flow excitation with harmonic base excitation. Hence, piezoaeroelastic evolution can be investigated in frequency domain for different airflow speeds and electrical boundary conditions. [DOI:10.1115/1.4002785]

Piezoaeroelastic Modeling and Analysis of a Generator Wing with Continuous and Segmented Electrodes

Unmanned air vehicles (UAVs) and micro air vehicles (MAVs) constitute unique application platforms for vibration-based energy harvesting. Generating usable electrical energy during their mission has the important practical value of providing an additional energy source to run small electronic components. Electrical energy can be harvested from aeroelastic vibrations of lifting surfaces of UAVs and MAVs as they tend to have relatively flexible wings compared to their larger counterparts. In this work, an electromechanically coupled finite element model is combined with an unsteady aerodynamic model to develop a piezoaeroelastic model for airflow excitation of cantilevered plates representing wing-like structures. The electrical power output and the displacement of the wing tip are investigated for several airflow speeds and two different electrode configurations (continuous and segmented). Cancelation of electrical output occurs for typical coupled bending-torsion aeroelastic modes of a cantilevered generator wing when continuous electrodes are used. Torsional motions of the coupled modes become relatively significant when segmented electrodes are used, improving the broadband performance and altering the flutter speed. Although the focus is placed on the electrical power that can be harvested for a given airflow speed, shunt damping effect of piezoelectric power generation is also investigated for both electrode configurations.

Molecular characterization of a miraculin-like gene differentially expressed during coffee development and coffee leaf miner infestation

The characterization of a coffee gene encoding a protein similar to miraculin-like proteins, which are members of the plant Kunitz serine trypsin inhibitor (STI) family of proteinase inhibitors (PIs), is described. PIs are important proteins in plant defence against insects and in the regulation of proteolysis during plant development. This gene has high identity with the Richadella dulcifica taste-modifying protein miraculin and with the tomato protein LeMir; and was named as CoMir (Coffea miraculin). Structural protein modelling indicated that CoMir had structural similarities with the Kunitz STI proteins, but suggested specific folding structures. CoMir was up-regulated after coffee leaf miner (Leucoptera coffella) oviposition in resistant plants of a progeny derived from crosses between C. racemosa (resistant) and C. arabica (susceptible). Interestingly, this gene was down-regulated during coffee leaf miner herbivory in susceptible plants. CoMir expression was up-regulated after abscisic acid application and wounding stress and was prominent during the early stages of flower and fruit development. In situ hybridization revealed that CoMir transcripts accumulated in the anther tissues that display programmed cell death (tapetum, endothecium and stomium) and in the metaxylem vessels of the petals, stigma and leaves. In addition, the recombinant protein CoMir shows inhibitory activity against trypsin. According to the present results CoMir may act in proteolytic regulation during coffee development and in the defence against L. coffeella. The similarity of CoMir with other Kunitz STI proteins and the role of CoMir in plant development and plant stress are discussed.

Application of molecular techniques to evaluate the methanogenic archaea and anaerobic bacteria in the presence of oxygen with different COD : sulfate ratios in a UASB reactor

In this paper, the microbial characteristics of the granular sludge in the presence of oxygen (3.0 +/- 0.7 mg O-2 1(-1)) were analyzed using molecular biology techniques. The granules were provided by an upflow anaerobic sludge blanket (UASB) operated over 469 days and fed with synthetic substrate. Ethanol and sulfate were added to obtain different COD/SO42- ratios (3.0, 2.0, and 1.6). The results of fluorescent in situ hybridization (FISH) analyses showed that archaeal cells, detected by the ARC915 probe, accounted for 77%, 84%, and 75% in the COD/SO42- ratios (3.0, 2.0, and 1.6, respectively). Methanosaeta sp. was the predominant acetoclastic archaea observed by optical microscopy and FISH analyses, and confirmed by sequencing of the excised bands of the DGGE gel with a similarity of 96%. The sulfate-reducing bacterium Desulfovibrio vulgaris subsp. vulgaris (similarity of 99%) was verified by sequencing of the DGGE band. Others identified microorganism were similar to Shewanella sp. and Desulfitobacterium hafniense, with similarities of 95% and 99%, respectively. These results confirmed that the presence of oxygen did not severely affect the metabolism of microorganisms that are commonly considered strictly anaerobic. We obtained mean efficiencies of organic matter conversion and sulfate reducing higher than 74%. (C) 2008 Elsevier Ltd. All rights reserved.