IDENTIFICATION OF GENES CONTROLLING BIOLOGICAL PROCESSES AND PATHWAYS THROUGH STATISTICAL ANALYSIS AND NETWORK RECONSTRUCTION

The developmental processes and functions of an organism are controlled by the genes and the proteins that are derived from these genes. The identification of key genes and the reconstruction of gene networks can provide a model to help us understand the regulatory mechanisms for the initiation and progression of biological processes or functional abnormalities (e.g. diseases) in living organisms. In this dissertation, I have developed statistical methods to identify the genes and transcription factors (TFs) involved in biological processes, constructed their regulatory networks, and also evaluated some existing association methods to find robust methods for coexpression analyses. Two kinds of data sets were used for this work: genotype data and gene expression microarray data. On the basis of these data sets, this dissertation has two major parts, together forming six chapters. The first part deals with developing association methods for rare variants using genotype data (chapter 4 and 5). The second part deals with developing and/or evaluating statistical methods to identify genes and TFs involved in biological processes, and construction of their regulatory networks using gene expression data (chapter 2, 3, and 6). For the first part, I have developed two methods to find the groupwise association of rare variants with given diseases or traits. The first method is based on kernel machine learning and can be applied to both quantitative as well as qualitative traits. Simulation results showed that the proposed method has improved power over the existing weighted sum method (WS) in most settings. The second method uses multiple phenotypes to select a few top significant genes. It then finds the association of each gene with each phenotype while controlling the population stratification by adjusting the data for ancestry using principal components. This method was applied to GAW 17 data and was able to find several disease risk genes. For the second part, I have worked on three problems. First problem involved evaluation of eight gene association methods. A very comprehensive comparison of these methods with further analysis clearly demonstrates the distinct and common performance of these eight gene association methods. For the second problem, an algorithm named the bottom-up graphical Gaussian model was developed to identify the TFs that regulate pathway genes and reconstruct their hierarchical regulatory networks. This algorithm has produced very significant results and it is the first report to produce such hierarchical networks for these pathways. The third problem dealt with developing another algorithm called the top-down graphical Gaussian model that identifies the network governed by a specific TF. The network produced by the algorithm is proven to be of very high accuracy.

Cladistic association analysis of Y chromosome effects on alcohol dependence and related personality traits

Association between Y chromosome haplotype variation and alcohol dependence and related personality traits was investigated in a large sample of psychiatrically diagnosed Finnish males. Haplotypes were constructed for 359 individuals using alleles at eight loci (seven microsatellite loci and a nucleotide substitution in the DYZ3 alphoid satellite locus). A cladogram linking the 102 observed haplotype configurations was constructed by using parsimony with a single-step mutation model. Then, a series of contingency tables nested according to the cladogram hierarchy were used to test for association between Y haplotype and alcohol dependence. Finally, using only alcohol-dependent subjects, we tested for association between Y haplotype and personality variables postulated to define subtypes of alcoholism—antisocial personality disorder, novelty seeking, harm avoidance, and reward dependence. Significant association with alcohol dependence was observed at three Y haplotype clades, with significance levels of P = 0.002, P = 0.020, and P = 0.010. Within alcohol-dependent subjects, no relationship was revealed between Y haplotype and antisocial personality disorder, novelty seeking, harm avoidance, or reward dependence. These results demonstrate, by using a fully objective association design, that differences among Y chromosomes contribute to variation in vulnerability to alcohol dependence. However, they do not demonstrate an association between Y haplotype and the personality variables thought to underlie the subtypes of alcoholism.

Biological significance of dead biomass retention trait in Mediterranean Basin species: an analysis between different successional niches and regeneration strategies as functional groups

Standing dead biomass retention is considered one of the most relevant fuel structural traits to affect plant flammability. However, very little is known about the biological significance of this trait and its distribution between different functional groups. Our aim was to analyse how the proportion of dead biomass produced in Mediterranean species is related to the successional niche of species (early-, mid- and late-successional stages) and the regeneration strategy of species (seeders and resprouters). We evaluated biomass distribution by size classes and standing dead biomass retention in nine dominant species from the Mediterranean Basin in different development stages (5, 9, 14 and 26 years since the last fire). The results revealed significant differences in the standing dead biomass retention of species that presented a distinct successional niche or regeneration strategy. These differences were restricted to the oldest ages studied (>9 years). Tree and small tree resprouters, typical in late-successional stages, presented slight variations with age and a less marked trend to retain dead biomass, while seeder shrubs and dwarf shrubs, characteristic of early-successional stages, showed high dead biomass loads. Our results suggest that the species that tend to retain more dead branches are colonising species that may promote fire in early-successional stages.

Production and body composition traits of broilers in relation to breast weight evaluated by path analysis

Breast weight has great economic importance in poultry industry, and may be associated with other variables. This work aimed to estimate phenotypic correlations between performance (live body weight at 7 and 28 days, and at slaughter, and depth of the breast muscle measured by ultrasonography), carcass (eviscerated body weight and leg weight) and body composition (heart, liver and abdominal fat weight) traits in a broiler line, and quantify the direct and indirect influence of these traits on breast weight. Path analysis was used by expanding the matrix of partial correlation in coefficients which give the direct influence of one trait on another, regardless the effect of the other traits. The simultaneous maintenance of live body weight at slaughter and eviscerated body weight in the matrix of correlations might be harmful for statistical analysis involving systems of normal equations, like path analysis, due to the observed multicollinearity. The live body weight at slaughter and the depth of the breast muscle as measured by ultrasonography directly affected breast weight and were identified as the most responsible factors for the magnitude of the correlation coefficients obtained between the studied traits and breast weight. Individual pre-selection for these traits could favor an increased breast weight in the future reproducer candidates of this line if the broilers' environmental conditions and housing are maintained, since the live body weight at slaughter and the depth of breast muscle measured by ultrasonography were directly related to breast weight.

Analysis of Bacteria, Parasites, and Heavy Metals in Lettuce (Lactuca sativa) and Rocket Salad (Eruca sativa L.) Irrigated with Treated Effluent from a Biological Wastewater Treatment Plant

This study aimed to evaluate the viability of using treated residuary water from the Biological Wastewater Treatment Plant of Ribeiro Preto to grow vegetables, through the characterization and quantification of parasites, coliforms, and heavy metals. Three equal cultivation areas were prepared. The first was irrigated with treated/chlorinated (0.2 mg L(-1)) wastewater, the second one with treated wastewater without chlorination, and the third site with potable water, which was the control group. The presence of Hymenolepis nana, Enterobius vermicularis, nematode larvae, and Entamoeba coli was verified in lettuce (Lactuca sativa) samples. Although nematode larvae were observed in rocket salad (Eruca sativa L.), no significant differences were found between the number of parasites and type of irrigation water used. No significant differences were found between the number of fecal coliforms in vegetables and the different types of irrigation. However, the vegetables irrigated with treated effluent without chlorination showed higher levels of fecal coliforms. The risk of pathogens is reduced with bleach addition to the treated effluent at 0.2 mg/L. Concentration of heavy metals in vegetables does not mean significant risks to human health, according with the parameters recommended by the World Health Organization.

Does using stepwise variable selection to build sequential path analysis models make sense?

Causal inference methods - mainly path analysis and structural equation modeling - offer plant physiologists information about cause-and-effect relationships among plant traits. Recently, an unusual approach to causal inference through stepwise variable selection has been proposed and used in various works on plant physiology. The approach should not be considered correct from a biological point of view. Here, it is explained why stepwise variable selection should not be used for causal inference, and shown what strange conclusions can be drawn based upon the former analysis when one aims to interpret cause-and-effect relationships among plant traits.

Power of the joint segregation analysis method for testing mixed major-gene and polygene inheritance models of quantitative traits

Understanding the genetic architecture of quantitative traits can greatly assist the design of strategies for their manipulation in plant-breeding programs. For a number of traits, genetic variation can be the result of segregation of a few major genes and many polygenes (minor genes). The joint segregation analysis (JSA) is a maximum-likelihood approach for fitting segregation models through the simultaneous use of phenotypic information from multiple generations. Our objective in this paper was to use computer simulation to quantify the power of the JSA method for testing the mixed-inheritance model for quantitative traits when it was applied to the six basic generations: both parents (P-1 and P-2), F-1, F-2, and both backcross generations (B-1 and B-2) derived from crossing the F-1 to each parent. A total of 1968 genetic model-experiment scenarios were considered in the simulation study to quantify the power of the method. Factors that interacted to influence the power of the JSA method to correctly detect genetic models were: (1) whether there were one or two major genes in combination with polygenes, (2) the heritability of the major genes and polygenes, (3) the level of dispersion of the major genes and polygenes between the two parents, and (4) the number of individuals examined in each generation (population size). The greatest levels of power were observed for the genetic models defined with simple inheritance; e.g., the power was greater than 90% for the one major gene model, regardless of the population size and major-gene heritability. Lower levels of power were observed for the genetic models with complex inheritance (major genes and polygenes), low heritability, small population sizes and a large dispersion of favourable genes among the two parents; e.g., the power was less than 5% for the two major-gene model with a heritability value of 0.3 and population sizes of 100 individuals. The JSA methodology was then applied to a previously studied sorghum data-set to investigate the genetic control of the putative drought resistance-trait osmotic adjustment in three crosses. The previous study concluded that there were two major genes segregating for osmotic adjustment in the three crosses. Application of the JSA method resulted in a change in the proposed genetic model. The presence of the two major genes was confirmed with the addition of an unspecified number of polygenes.

Identification of causal relationships among traits related to drought resistance in Stylosanthes scabra using QTL analysis

Previous studies have shown that a negative relationship exists between transpiration efficiency (TE) and carbon isotope discrimination (Delta) and between TE and specific leaf area (SLA) in Stylosanthes scabra, A glasshouse experiment was conducted to confirm these relationships in an F-2 population and to study the causal nature of these relationships through quantitative trait loci (QTL) analysis, One hundred and twenty F-2 genotypes from a cross between two genotypes within S. scabra were used. Three replications for each genotype were maintained through vegetative propagation, Water stress was imposed by maintaining plants at 40% of field capacity for about 45 d. To facilitate QTL analysis, a genetic linkage map consisting of 151 RAPD markers was developed, Results from this study show that Delta was significantly and negatively correlated with TE and biomass production. Similarly, SLA showed significant negative correlation with TE and biomass production, Most of the QTL for TE and Delta were present on linkage groups 5 and 11. Similarly, QTL for SLA, transpiration and biomass productivity traits were clustered on linkage groups 13 and 24, One unlinked marker was also associated with these traits, There were several markers coincident between different traits, At all the coincident QTL, the direction of QTL effects was consistent with phenotypic data, At the coincident markers between TE and Delta, high alleles of TE were associated with low alleles of Delta. Similarly, low alleles of SLA were associated with high alleles of biomass productivity traits and transpiration. At the coincident markers between trans-4-hydroxy-N-methyl proline (MHP) and relative water content (RWC), low alleles of MHP were associated with high alleles of RWC, This study suggests the causal nature of the relationship between TE and Delta. Phenotypic data and QTL, data show that SLA was more closely associated with biomass production than with TE, This study also shows that a cause-effect relationship may exist between SLA and biomass production.

Analysis of the microbial community structure and function of a laboratory scale enhanced biological phosphorus removal reactor

A laboratory scale sequencing batch reactor (SBR) operating for enhanced biological phosphorus removal (EBPR) and fed with a mixture of volatile fatty acids (VFAs) showed stable and efficient EBPR capacity over a four-year-period. Phosphorus (P), poly-beta-hydroxyalkanoate (PHA) and glycogen cycling consistent with classical anaerobic/aerobic EBPR were demonstrated with the order of anaerobic VFA uptake being propionate, acetate then butyrate. The SBR was operated without pH control and 63.67+/-13.86 mg P l(-1) was released anaerobically. The P% of the sludge fluctuated between 6% and 10% over the operating period (average of 8.04+/-1.31%). Four main morphological types of floc-forming bacteria were observed in the sludge during one year of in-tensive microscopic observation. Two of them were mainly responsible for anaerobic/aerobic P and PHA transformations. Fluorescence in situ hybridization (FISH) and post-FISH chemical staining for intracellular polyphosphate and PHA were used to determine that 'Candidatus Accumulibacter phosphatis' was the most abundant polyphosphate accumulating organism (PAO), forming large clusters of coccobacilli (1.0-1.5 mum) and comprising 53% of the sludge bacteria. Also by these methods, large coccobacillus-shaped gammaproteobacteria (2.5-3.5 mum) from a recently described novel cluster were glycogen-accumulating organisms (GAOs) comprising 13% of the bacteria. Tetrad-forming organisms (TFOs) consistent with the 'G bacterium' morphotype were alphaproteobacteria , but not Amaricoccus spp., and comprised 25% of all bacteria. According to chemical staining, TFOs were occasionally able to store PHA anaerobically and utilize it aerobically.

Development of a novel titration and off-gas analysis (TOGA) sensor for study of biological processes in wastewater treatment systems

The development of the new TOGA (titration and off-gas analysis) sensor for the detailed study of biological processes in wastewater treatment systems is outlined. The main innovation of the sensor is the amalgamation of titrimetric and off-gas measurement techniques. The resulting measured signals are: hydrogen ion production rate (HPR), oxygen transfer rate (OTR), nitrogen transfer rate (NTR), and carbon dioxide transfer rate (CTR). While OTR and NTR are applicable to aerobic and anoxic conditions, respectively, HPR and CTR are useful signals under all of the conditions found in biological wastewater treatment systems, namely, aerobic, anoxic and anaerobic. The sensor is therefore a powerful tool for studying the key biological processes under all these conditions. A major benefit from the integration of the titrimetric and off-gas analysis methods is that the acid/base buffering systems, in particular the bicarbonate system, are properly accounted for. Experimental data resulting from the TOGA sensor in aerobic, anoxic, and anaerobic conditions demonstrates the strength of the new sensor. In the aerobic environment, carbon oxidation (using acetate as an example carbon source) and nitrification are studied. Both the carbon and ammonia removal rates measured by the sensor compare very well with those obtained from off-line chemical analysis. Further, the aerobic acetate removal process is examined at a fundamental level using the metabolic pathway and stoichiometry established in the literature, whereby the rate of formation of storage products is identified. Under anoxic conditions, the denitrification process is monitored and, again, the measured rate of nitrogen gas transfer (NTR) matches well with the removal of the oxidised nitrogen compounds (measured chemically). In the anaerobic environment, the enhanced biological phosphorus process was investigated. In this case, the measured sensor signals (HPR and CTR) resulting from acetate uptake were used to determine the ratio of the rates of carbon dioxide production by competing groups of microorganisms, which consequently is a measure of the activity of these organisms. The sensor involves the use of expensive equipment such as a mass spectrometer and requires special gases to operate, thus incurring significant capital and operational costs. This makes the sensor more an advanced laboratory tool than an on-line sensor. (C) 2003 Wiley Periodicals, Inc.

Phenotypic correlations and path analysis for plant architecture traits and grain production in three generations of cowpea

ABSTRACT The objective of this study was to analyze the phenotypic correlation and path analysis of traits related to plant architecture, earliness and grain yield in F2, BC1 and BC2 generations, from crosses between cowpea cultivars BRS Carijó and BR14 Mulato. Most phenotypic correlations of the examined traits were concordant in statistical significance, with approximate values ​​among the examined generations. For the trait seed weight, significant and positive phenotypic correlations were observed in the three generations only for the trait number of secondary branches. The values ​​of the direct effects were in agreement with the values ​​of the phenotypic correlations, which indicate true association by the phenotypic correlation among the traits of grain yield examined. Path analysis indicated that the selection of productive plants will result in early plants and an increased number of secondary branches. In F2, plants with shorter length of the main branch and shorter length of secondary branches can be obtained. The causal model explained 15 to 30% of the total variation in grain weight in relation to the traits examined. The analyses indicated the possibility of selecting plants with a higher and early grain yield, shorter length of primary branches and lower number of nodes, which are important variables for mechanical or semi-mechanical harvesting.

Biological and genetic characteristics of uropathogenic Escherichia coli strains

The aim of the present study was to determine biological characteristics such as expression of fimbriae, Congo red binding, production of hemolysin and aerobactin, adhesion to HeLa and uroepithelial cells and invasion of HeLa cells by Escherichia coli isolates obtained from patients showing clinical signs of urinary tract infection (UTI). Also, the presence of genes (apa, afa, spa) for fimbria expression and cytotoxic necrotizing factors (CNF1, CNF2) was assayed using specific primers in PCR. The data obtained were compared with the clonal relationships obtained by analysis of multilocus enzyme electrophoresis (MLEE), restriction fragment length polymorphism (RFLP) of the rDNA (ribotyping) and enterobacterial repetitive intergenic consensus-PCR (ERIC-PCR). All isolates but one presented a combination of at least two of the characteristics studied, a fact suggesting the presence of pathogenicity islands (PAIs). Diffuse adherence type to HeLa cells was observed to occur in most of the strains, but adhesion to uroepithelial cells seems to be a more reliable test to verify pathogenicity. Although four strains seemed to be able to invade HeLa cells when assayed by light microscopy, electron microscopy studies demonstrated that these strains were not invasive. MLEE, RFLP and ERIC-PCR were able to group the isolates differently into main clusters that were not correlated with the presence of pathogenic traits.

Clonal structure of Trypanosoma cruzi Colombian strain (biodeme Type III): biological, isoenzymic and histopathological analysis of seven isolated clones

The clonal structure of the Colombian strain of Trypanosoma cruzi, biodeme Type III and zymodeme 1, was analyzed in order to characterize its populations and to establish its homogeneity or heterogeneity. Seven isolated clones presented the basic characteristics of Biodeme Type III, with the same patterns of parasitemic curves, tissue tropism to skeletal muscle and myocardium, high pathogenicity with extensive necrotic-inflammatory lesions from the 20th to 30th day of infection. The parental strain and its clones C1, C3, C4 and C6, determined the higher levels of parasitemia, 20 to 30 days of infection, with high mortality rate up to 30 days (79 to 100%); clones C2, C5 and C7 presented lower levels of parasitemia, with low mortality rates (7.6 to 23%). Isoenzymic patterns, characteristic of zymodeme 1, (Z1) were similar for the parental strain and its seven clones. Results point to a phenotypic homogeneity of the clones isolated from the Colombian strain and suggest the predominance of a principal clone, responsible for the biological behavior of the parental strain and clones.