Exploring the multifactorial nature of autism through computational systems biology: Calcium and the Rho GTPase RAC1 under the spotlight

Autism is a neurodevelopmental disorder characterized by impaired social interaction and communication accompanied with repetitive behavioral patterns and unusual stereotyped interests. Autism is considered a highly heterogeneous disorder with diverse putative causes and associated factors giving rise to variable ranges of symptomatology. Incidence seems to be increasing with time, while the underlying pathophysiological mechanisms remain virtually uncharacterized (or unknown). By systematic review of the literature and a systems biology approach, our aims were to examine the multifactorial nature of autism with its broad range of severity, to ascertain the predominant biological processes, cellular components, and molecular functions integral to the disorder, and finally, to elucidate the most central contributions (genetic and/or environmental) in silico. With this goal, we developed an integrative network model for gene-environment interactions (GENVI model) where calcium (Ca2+) was shown to be its most relevant node. Moreover, considering the present data from our systems biology approach together with the results from the differential gene expression analysis of cerebellar samples from autistic patients, we believe that RAC1, in particular, and the RHO family of GTPases, in general, could play a critical role in the neuropathological events associated with autism. © 2013 Springer Science+Business Media New York.

Estudo da interação genótipo-ambiente em rebanhos da raça Nelore na Amazônia Legal

Dados de pesos padronizados aos 120 (P120), 210 (P210), 450 (P450) dias de idade, perímetro escrotal aos 450 dias de idade (PE450) e idade ao primeiro parto (IPP), de 211.744 registros de animais Nelore, provenientes de fazenda localizadas na região da Amazônia Legal, foram utilizados na análise. O efeito da interação genótipo-ambiente foi estudado por meio de estimativas de herdabilidade e de correlações entre classificações, comparando os animais da Amazônia Legal com a base geral de animais do PMGRN – Nelore Brasil. As análises bi-característica consideraram o P120 como característica-âncora, com P210, P450 e PE450. O modelo de análise considerou grupo de contemporâneos de P120, P210, P450 e da classe de idade da vaca ao parto (CIVP) como efeitos fixos e os efeitos aleatórios genéticos aditivos diretos, maternos e residuais. Nas análises de P120 com P450 e com PE450 foi desconsiderado o efeito materno no modelo. A covariância aditivo-materna foi fixada em zero, conforme protocolo estabelecido nas análises do PMGRN – Nelore Brasil. A característica IPP foi analisada separadamente, em análise de característica única e considerando como efeito fixo o GCIPP e como aleatórios os efeitos aditivos genéticos e residual. Para realizar a comparação entre classificações, por meio do procedimento PROC CORR opção spearman do SAS. As estimativas de herdabilidade para P120, P210, P450, PE450 e IPP nos dados da Amazônia Legal foram: 0,20 a 0,49; 0,21; 0,48; 0,45 e 0,21, respectivamente, e nos dados gerais do PMGRN – Nelore Brasil foram: 0,23; 0,25; 0,34; 0,43 e 0,11, respectivamente. As correlações entre classificações de rank para P120, P210, P450, PE450 e IPP foram iguais a 0,77; 0,79; 0,82; 0,78 e 0,38, respectivamente. As análises da interação genótipo-ambiente evidenciaram maiores efeitos sobre os aspectos maternos, de peso ao sobreano e IPP, enquanto que as correlações entre classificações mostraram fortes evidências de interação genótipo-ambiente em quase todas as características estudadas.

Application of Reactions Norms in Study of Genotype Environmental Interaction for Milk Yield of Buffaloes

The aim of this study was to estimate genetic parameters for milk yield (MY) in buffaloes using reaction norms. Model included the additive direct effect as random and contemporary group (herd and year of birth) were included as fixed effects and cow age classes (linear) as covariables. The animal additive direct random effect was modeled through linear Legendre polynomials on environment gradient (EG) standardized means. Mean trends were taken into account by a linear regression on Legendre polynomials of environmental group means. Residual variance was modeled trough 6 heterogeneity classes (EG). These classes of residual variance was formed : EG1: mean = 866,93 kg (621,68 kg-1011,76 kg); EG2: mean = 1193,00 kg (1011,76 kg-1251,49 kg); EG3: mean = 1309,37 kg (1251,49 kg -1393,20 kg); EG4: mean = 1497,59 kg (1393,20 kg-1593,53 kg); EG5: mean = 1664,78 kg (1593,53 kg -1727,32kg) e EG6: mean = 1973,85 kg (1727,32 kg -2422,19 kg).(Co) variance functions were estimated by restricted maximum likelihood (REML) using the GIBBS3F90 package. The heritability estimates for MY raised as the environmental gradient increased, varying from 0.20 to 0.40. However, in intermediate to favorable environments, the heritability estimates obtained with Considerable genotype-environment interaction was found for MY using reaction norms. For genetic evaluation of MY is necessary to consider heterogeneity of variances to model the residual variance.

Genetic parameters and investigation of genotype x environment interactions in Nellore x Hereford crossbred for resistance to cattle ticks in different regions of Brazil

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Mapping of angular leaf spot resistance QTL in common bean (Phaseolus vulgaris L.) under different environments

Background: Common bean (Phaseolus vulgaris L.) is the most important grain legume for human diet worldwide and the angular leaf spot (ALS) is one of the most devastating diseases of this crop, leading to yield losses as high as 80%. In an attempt to breed resistant cultivars, it is important to first understand the inheritance mode of resistance and to develop tools that could be used in assisted breeding. Therefore, the aim of this study was to identify quantitative trait loci (QTL) controlling resistance to ALS under natural infection conditions in the field and under inoculated conditions in the greenhouse. Results: QTL analyses were made using phenotypic data from 346 recombinant inbreed lines from the IAC-UNA x CAL 143 cross, gathered in three experiments, two of which were conducted in the field in different seasons and one in the greenhouse. Joint composite interval mapping analysis of QTL x environment interaction was performed. In all, seven QTLs were mapped on five linkage groups. Most of them, with the exception of two, were significant in all experiments. Among these, ALS10.1(DG,UC) presented major effects (R-2 between 16% - 22%). This QTL was found linked to the GATS11b marker of linkage group B10, which was consistently amplified across a set of common bean lines and was associated with the resistance. Four new QTLs were identified. Between them the ALS5.2 showed an important effect (9.4%) under inoculated conditions in the greenhouse. ALS4.2 was another major QTL, under natural infection in the field, explaining 10.8% of the variability for resistance reaction. The other QTLs showed minor effects on resistance. Conclusions: The results indicated a quantitative inheritance pattern of ALS resistance in the common bean line CAL 143. QTL x environment interactions were observed. Moreover, the major QTL identified on linkage group B10 could be important for bean breeding, as it was stable in all the environments. Thereby, the GATS11b marker is a potential tool for marker assisted selection for ALS resistance.

Studying the genetic basis of drought tolerance in sorghum by managed stress trials and adjustments for phenological and plant height differences

Managed environments in the form of well watered and water stressed trials were performed to study the genetic basis of grain yield and stay green in sorghum with the objective of validating previously detected QTL. As variations in phenology and plant height may influence QTL detection for the target traits, QTL for flowering time and plant height were introduced as cofactors in QTL analyses for yield and stay green. All but one of the flowering time QTL were detected near yield and stay green QTL. Similar co-localization was observed for two plant height QTL. QTL analysis for yield, using flowering time/plant height cofactors, led to yield QTL on chromosomes 2, 3, 6, 8 and 10. For stay green, QTL on chromosomes 3, 4, 8 and 10 were not related to differences in flowering time/plant height. The physical positions for markers in QTL regions projected on the sorghum genome suggest that the previously detected plant height QTL, Sb-HT9-1, and Dw2, in addition to the maturity gene, Ma5, had a major confounding impact on the expression of yield and stay green QTL. Co-localization between an apparently novel stay green QTL and a yield QTL on chromosome 3 suggests there is potential for indirect selection based on stay green to improve drought tolerance in sorghum. Our QTL study was carried out with a moderately sized population and spanned a limited geographic range, but still the results strongly emphasize the necessity of corrections for phenology in QTL mapping for drought tolerance traits in sorghum.

Epistatic effects on grain yield of soybean [Glycine max (L.) Merrill]

Studies addressing the estimation of genetic parameters in soybean have not emphasized the epistatic effects. The purpose of this study was to estimate the significance of these effects on soybean grain yield, based on the Modified Triple Test Cross design. Thirty-two inbred lines derived from a cross between two contrasting lines were used, which were crossed with two testers (L1 and L2). The experiments were carried out at two locations, in 10 x 10 triple lattice designs with 9 replications, containing 32 lines (Pi ), 64 crosses (32 Pi x L1 and 32 Pi x L2 ) and controls. The variation between ( ͞L1i + ͞L2i - ͞Pi ) revealed the presence of epistasis, as well as an interaction of epistasis x environment. Since the predominant component of epistasis in autogamous species is additive x additive (i type), we suggest postponing the selection for grain yield to later generations of inbreeding in order to exploit the beneficial effects of additive x additive epistasis.

STATISTICAL METHODS FOR RARE AND COMMON VARIANT ASSOCIATION STUDIES

Complex human diseases are a major challenge for biological research. The goal of my research is to develop effective methods for biostatistics in order to create more opportunities for the prevention and cure of human diseases. This dissertation proposes statistical technologies that have the ability of being adapted to sequencing data in family-based designs, and that account for joint effects as well as gene-gene and gene-environment interactions in the GWA studies. The framework includes statistical methods for rare and common variant association studies. Although next-generation DNA sequencing technologies have made rare variant association studies feasible, the development of powerful statistical methods for rare variant association studies is still underway. Chapter 2 demonstrates two adaptive weighting methods for rare variant association studies based on family data for quantitative traits. The results show that both proposed methods are robust to population stratification, robust to the direction and magnitude of the effects of causal variants, and more powerful than the methods using weights suggested by Madsen and Browning [2009]. In Chapter 3, I extended the previously proposed test for Testing the effect of an Optimally Weighted combination of variants (TOW) [Sha et al., 2012] for unrelated individuals to TOW &ndash F, TOW for Family &ndash based design. Simulation results show that TOW &ndash F can control for population stratification in wide range of population structures including spatially structured populations, is robust to the directions of effect of causal variants, and is relatively robust to percentage of neutral variants. In GWA studies, this dissertation consists of a two &ndash locus joint effect analysis and a two-stage approach accounting for gene &ndash gene and gene &ndash environment interaction. Chapter 4 proposes a novel two &ndash stage approach, which is promising to identify joint effects, especially for monotonic models. The proposed approach outperforms a single &ndash marker method and a regular two &ndash stage analysis based on the two &ndash locus genotypic test. In Chapter 5, I proposed a gene &ndash based two &ndash stage approach to identify gene &ndash gene and gene &ndash environment interactions in GWA studies which can include rare variants. The two &ndash stage approach is applied to the GAW 17 dataset to identify the interaction between KDR gene and smoking status.

Concept, design and implementation of a cardiovascular gene-centric 50 k SNP array for large-scale genomic association studies.

A wealth of genetic associations for cardiovascular and metabolic phenotypes in humans has been accumulating over the last decade, in particular a large number of loci derived from recent genome wide association studies (GWAS). True complex disease-associated loci often exert modest effects, so their delineation currently requires integration of diverse phenotypic data from large studies to ensure robust meta-analyses. We have designed a gene-centric 50 K single nucleotide polymorphism (SNP) array to assess potentially relevant loci across a range of cardiovascular, metabolic and inflammatory syndromes. The array utilizes a "cosmopolitan" tagging approach to capture the genetic diversity across approximately 2,000 loci in populations represented in the HapMap and SeattleSNPs projects. The array content is informed by GWAS of vascular and inflammatory disease, expression quantitative trait loci implicated in atherosclerosis, pathway based approaches and comprehensive literature searching. The custom flexibility of the array platform facilitated interrogation of loci at differing stringencies, according to a gene prioritization strategy that allows saturation of high priority loci with a greater density of markers than the existing GWAS tools, particularly in African HapMap samples. We also demonstrate that the IBC array can be used to complement GWAS, increasing coverage in high priority CVD-related loci across all major HapMap populations. DNA from over 200,000 extensively phenotyped individuals will be genotyped with this array with a significant portion of the generated data being released into the academic domain facilitating in silico replication attempts, analyses of rare variants and cross-cohort meta-analyses in diverse populations. These datasets will also facilitate more robust secondary analyses, such as explorations with alternative genetic models, epistasis and gene-environment interactions.

Intentions: a confident-based interaction design for smart spaces

The paradigm of ubiquitous computing has become a reference for the design of Smart Spaces. Current trends in Ambient Intelligence are increasingly related to the scope of Internet of Things. This paradigm has the potential to support cost-effective solutions in the fields of telecare, e-health and Ambient Assisted Living. Nevertheless, ubiquitous computing does not provide end users with a role for proactive interactions with the environment. Thus, the deployment of smart health care services at a private space like the home is still unsolved. This PhD dissertation aims to define a person-environment interaction model to foster acceptability and users confidence in private spaces by applying the concept of user-centred security and the human performance model of seven stages of action.

Case-only study of interactions between genetic polymorphisms of GSTM1, P1, T1 and Z1 and smoking in Parkinson's disease

Current opinion contends that complex interactions between genetic and environmental factors play a role in the etiology of Parkinson's disease (PD). Cigarette smoking is thought to reduce risk of PD, and emerging evidence suggests that genetic factors may modulate smoking's effect. We used a case-only design, an approach not previously used to study gene-environment interactions in PD, specifically to study interactions between glutathione-S-transferase (GST) gene polymorphisms and smoking in relation to PD. Four-hundred PD cases (age at onset: 60.0 +/- 10.7 years) were genotyped for common polymorphisms in GSTM1, PI, T1 and Z1 using well-established methods. Smoking exposure data were collected in face-to-face interviews. The independence of the studied GST genotypes and smoking exposure was confirmed by studying 402 healthy, aged individuals. No differences were observed in the distributions of GSTM1, T1 or Z1 polymorphisms between ever-smoked and never-smoked PD cases using logistic regression (all P > 0.43). However, GSTP1 *C haplotypes were over-represented among PD cases who ever smoked (odds ratio for interaction (ORi) = 2.00 (95% Cl: 1.11-3.60, P = 0.03)). Analysis revealed that ORi between smoking and the GSTP1-114Val carrier status increased with increasing smoking dose (P = 0.02 for trend). These data suggest that one or more GSTP1 polymorphisms may interact with cigarette smoking to influence the risk for PD. (C) 2004 Elsevier Ireland Ltd. All rights reserved.

Effect of environment and genotype on bread-making quality of spring-sown spring wheat cultivars in China

Improvement of end-use quality in bread wheat depends on a thorough understanding of current wheat quality and the influences of genotype (G), environment (E), and genotype by environment interaction (G x E) on quality traits. Thirty-nine spring-sown spring wheat (SSSW) cultivars and advanced lines from China were grown in four agro-ecological zones comprising seven locations during the 1998 and 1999 cropping seasons. Data on 12 major bread-making quality traits were used to investigate the effect of G, E, and G x E on these traits. Wide range variability for protein quantity and quality, starch quality parameters and milling quality in Chinese SSSW was observed. Genotype and environment were found to significantly influence all quality parameters as major effects. Kernel hardness, flour yield, Zeleny sedimentation value and mixograph properties were mainly influenced by the genetic variance components, while thousand kernel weight, test weight, and falling number were mostly influenced by the environmental variance components. Genotype, environment, and their interaction had important effects on test weight, mixing development time and RVA parameters. Cultivars originating from Zone VI (northeast) generally expressed high kernel hardness, good starch quality, but poor milling and medium to weak mixograph performance; those from Zone VII (north) medium to good gluten and starch quality, but low milling quality; those from Zone VIII (central northwest) medium milling and starch quality, and medium to strong mixograph performance; those from Zone IX (western/southwestern Qinghai-Tibetan Plateau) medium milling quality, but poor gluten strength and starch parameters; and those from Zone X (northwest) high milling quality, strong mixograph properties, but low protein content. Samples from Harbin are characterized by good gluten and starch quality, but medium to poor milling quality; those from Hongxinglong by strong mixograph properties, medium to high milling quality, but medium to poor starch quality and medium to low protein content; those from Hohhot by good gluten but poor milling quality; those from Linhe by weak gluten quality, medium to poor milling quality; those from Lanzhou by poor bread-making and starch quality; those from Yongning by acceptable bread-making and starch quality and good milling quality; and those from Urumqi by good milling quality, medium gluten quality and good starch pasting parameters. Our findings suggest that Chinese SSSW quality could be greatly enhanced through genetic improvement for targeted well-characterized production environments.

Mixed model formulations for multi-environment trials

When studying genotype X environment interaction in multi-environment trials, plant breeders and geneticists often consider one of the effects, environments or genotypes, to be fixed and the other to be random. However, there are two main formulations for variance component estimation for the mixed model situation, referred to as the unconstrained-parameters (UP) and constrained-parameters (CP) formulations. These formulations give different estimates of genetic correlation and heritability as well as different tests of significance for the random effects factor. The definition of main effects and interactions and the consequences of such definitions should be clearly understood, and the selected formulation should be consistent for both fixed and random effects. A discussion of the practical outcomes of using the two formulations in the analysis of balanced data from multi-environment trials is presented. It is recommended that the CP formulation be used because of the meaning of its parameters and the corresponding variance components. When managed (fixed) environments are considered, users will have more confidence in prediction for them but will not be overconfident in prediction in the target (random) environments. Genetic gain (predicted response to selection in the target environments from the managed environments) is independent of formulation.

Trait physiology and crop modelling as a framework to link phenotypic complexity to underlying genetic systems

New tools derived from advances in molecular biology have not been widely adopted in plant breeding for complex traits because of the inability to connect information at gene level to the phenotype in a manner that is useful for selection. In this study, we explored whether physiological dissection and integrative modelling of complex traits could link phenotype complexity to underlying genetic systems in a way that enhanced the power of molecular breeding strategies. A crop and breeding system simulation study on sorghum, which involved variation in 4 key adaptive traits-phenology, osmotic adjustment, transpiration efficiency, stay-green-and a broad range of production environments in north-eastern Australia, was used. The full matrix of simulated phenotypes, which consisted of 547 location-season combinations and 4235 genotypic expression states, was analysed for genetic and environmental effects. The analysis was conducted in stages assuming gradually increased understanding of gene-to-phenotype relationships, which would arise from physiological dissection and modelling. It was found that environmental characterisation and physiological knowledge helped to explain and unravel gene and environment context dependencies in the data. Based on the analyses of gene effects, a range of marker-assisted selection breeding strategies was simulated. It was shown that the inclusion of knowledge resulting from trait physiology and modelling generated an enhanced rate of yield advance over cycles of selection. This occurred because the knowledge associated with component trait physiology and extrapolation to the target population of environments by modelling removed confounding effects associated with environment and gene context dependencies for the markers used. Developing and implementing this gene-to-phenotype capability in crop improvement requires enhanced attention to phenotyping, ecophysiological modelling, and validation studies to test the stability of candidate genetic regions.

Analysis of line x environment interactions for yield in navy beans. 2. Pattern analysis of lines and environment within years

Seven years of multi-environment yield trials of navy bean (Phaseolus vulgaris L.) grown in Queensland were examined. As is common with plant breeding evaluation trials, test entries and locations varied between years. Grain yield data were analysed for each year using cluster and ordination analyses (pattern analyses). These methods facilitate descriptions of genotype performance across environments and the discrimination among genotypes provided by the environments. The observed trends for genotypic yield performance across environments were partly consistent with agronomic and disease reactions at specific environments and also partly explainable by breeding and selection history. In some cases, similarities in discrimination among environments were related to geographic proximity, in others management practices, and in others similarities occurred between geographically widely separated environments which differed in management practices. One location was identified as having atypical line discrimination. The analysis indicated that the number of test locations was below requirements for adequate representation of line x environment interaction. The pattern analyses methods used were an effective aid in describing the patterns in data for each year and illustrated the variations in adaptive patterns from year to year. The study has implications for assessing the number and location of test sites for plant breeding multi-environment trials, and for the understanding of genetic traits contributing to line x environment interactions.