Genetic group x ambient temperature interaction effects on physiological responses and growth performance of rabbits

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Polymorphisms of DNA repair genes XRCC1 and XRCC3, interaction with environmental exposure and risk of chronic gastritis and gastric cancer

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Genotype x location x feeding interactions in Nellore cattle calculated through three-mode principal component analysis

The main objective of this study was to apply three-mode principal component analysis to assess the triple interaction (genotype x location x feeding) on direct genetic value for weight at 205 days of age. We used 60 sires with offspring in three regions of northeastern Brazil (Maranhao, Mata and Agreste, and Reconcavo Baiano) and raised on a pasture regime or with supplementation. There was no interaction between genotype and location, but there was a correlation between genotype and direct effect of feeding. The use of sires should be dictated according to the system of rearing of their offspring.

BAYESIAN STATISTICAL METHODS IN GENE-ENVIRONMENT AND GENE-GENE INTERACTION STUDIES

Complex diseases such as cancer result from multiple genetic changes and environmental exposures. Due to the rapid development of genotyping and sequencing technologies, we are now able to more accurately assess causal effects of many genetic and environmental factors. Genome-wide association studies have been able to localize many causal genetic variants predisposing to certain diseases. However, these studies only explain a small portion of variations in the heritability of diseases. More advanced statistical models are urgently needed to identify and characterize some additional genetic and environmental factors and their interactions, which will enable us to better understand the causes of complex diseases. In the past decade, thanks to the increasing computational capabilities and novel statistical developments, Bayesian methods have been widely applied in the genetics/genomics researches and demonstrating superiority over some regular approaches in certain research areas. Gene-environment and gene-gene interaction studies are among the areas where Bayesian methods may fully exert its functionalities and advantages. This dissertation focuses on developing new Bayesian statistical methods for data analysis with complex gene-environment and gene-gene interactions, as well as extending some existing methods for gene-environment interactions to other related areas. It includes three sections: (1) Deriving the Bayesian variable selection framework for the hierarchical gene-environment and gene-gene interactions; (2) Developing the Bayesian Natural and Orthogonal Interaction (NOIA) models for gene-environment interactions; and (3) extending the applications of two Bayesian statistical methods which were developed for gene-environment interaction studies, to other related types of studies such as adaptive borrowing historical data. We propose a Bayesian hierarchical mixture model framework that allows us to investigate the genetic and environmental effects, gene by gene interactions (epistasis) and gene by environment interactions in the same model. It is well known that, in many practical situations, there exists a natural hierarchical structure between the main effects and interactions in the linear model. Here we propose a model that incorporates this hierarchical structure into the Bayesian mixture model, such that the irrelevant interaction effects can be removed more efficiently, resulting in more robust, parsimonious and powerful models. We evaluate both of the 'strong hierarchical' and 'weak hierarchical' models, which specify that both or one of the main effects between interacting factors must be present for the interactions to be included in the model. The extensive simulation results show that the proposed strong and weak hierarchical mixture models control the proportion of false positive discoveries and yield a powerful approach to identify the predisposing main effects and interactions in the studies with complex gene-environment and gene-gene interactions. We also compare these two models with the 'independent' model that does not impose this hierarchical constraint and observe their superior performances in most of the considered situations. The proposed models are implemented in the real data analysis of gene and environment interactions in the cases of lung cancer and cutaneous melanoma case-control studies. The Bayesian statistical models enjoy the properties of being allowed to incorporate useful prior information in the modeling process. Moreover, the Bayesian mixture model outperforms the multivariate logistic model in terms of the performances on the parameter estimation and variable selection in most cases. Our proposed models hold the hierarchical constraints, that further improve the Bayesian mixture model by reducing the proportion of false positive findings among the identified interactions and successfully identifying the reported associations. This is practically appealing for the study of investigating the causal factors from a moderate number of candidate genetic and environmental factors along with a relatively large number of interactions. The natural and orthogonal interaction (NOIA) models of genetic effects have previously been developed to provide an analysis framework, by which the estimates of effects for a quantitative trait are statistically orthogonal regardless of the existence of Hardy-Weinberg Equilibrium (HWE) within loci. Ma et al. (2012) recently developed a NOIA model for the gene-environment interaction studies and have shown the advantages of using the model for detecting the true main effects and interactions, compared with the usual functional model. In this project, we propose a novel Bayesian statistical model that combines the Bayesian hierarchical mixture model with the NOIA statistical model and the usual functional model. The proposed Bayesian NOIA model demonstrates more power at detecting the non-null effects with higher marginal posterior probabilities. Also, we review two Bayesian statistical models (Bayesian empirical shrinkage-type estimator and Bayesian model averaging), which were developed for the gene-environment interaction studies. Inspired by these Bayesian models, we develop two novel statistical methods that are able to handle the related problems such as borrowing data from historical studies. The proposed methods are analogous to the methods for the gene-environment interactions on behalf of the success on balancing the statistical efficiency and bias in a unified model. By extensive simulation studies, we compare the operating characteristics of the proposed models with the existing models including the hierarchical meta-analysis model. The results show that the proposed approaches adaptively borrow the historical data in a data-driven way. These novel models may have a broad range of statistical applications in both of genetic/genomic and clinical studies.

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.

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.

Markers of genetic susceptibility in human environmental hygiene and toxicology : the role of selected CYP, NAT and GST genes

Inherited genetic traits co-determine the susceptibility of an individual to a toxic chemical. Special emphasis has been put on individual responses to environmental and industrial carcinogens, but other chronic diseases are of increasing interest. Polymorphisms of relevant xenobiotic metabolising enzymes may be used as toxicological susceptibility markers. A growing number of genes encoding enzymes involved in biotransformation of toxicants and in cellular defence against toxicant-induced damage to the cells has been identified and cloned, leading to increased knowledge of allelic variants of genes and genetic defects that may result in a differential susceptibility toward environmental toxicants. "Low penetrating" polymorphisms in metabolism genes tend to be much more common in the population than allelic variants of "high penetrating" cancer genes, and are therefore of considerable importance from a public health point of view. Positive associations between cancer and CYP1A1 alleles, in particular the *2C I462V allele, were found for tissues following the aerodigestive tract. Again, in most cases, the effect of the variant CYP1A1 allele becomes apparent or clearer in connection with the GSTM1 null allele. The CYP1B1 codon 432 polymorphism (CYP1B1*3) has been identified as a susceptibility factor in smoking-related head-and-neck squameous cell cancer. The impact of this polymorphic variant of CYP1B1 on cancer risk was also reflected by an association with the frequency of somatic mutations of the p53 gene. Combined genotype analysis of CYP1B1 and the glutathione transferases GSTM1 or GSTT1 has also pointed to interactive effects. Of particular interest for the industrial and environmental field is the isozyme CYP2E1. Several genotypes of this isozyme have been characterised which seem to be associated with different levels of expression of enzyme activity. The acetylator status for NAT2 can be determined by genotyping or by phenotyping. In the pathogenesis of human bladder cancer due to occupational exposure to "classical" aromatic amines (benzidine, 4-aminodiphenyl, 1-naphthylamine) acetylation by NAT2 is regarded as a detoxication step. Interestingly, the underlying European findings of a higher susceptibility of slow acetylators towards aromatic amines are in contrast to findings in Chinese workers occupationally exposed to aromatic amines which points to different mechanisms of susceptibility between European and Chinese populations. Regarding human bladder cancer, the hypothesis has been put forward that genetic polymorphism of GSTM1 might be linked with the occurrence of this tumour type. This supports the hypothesis that exposure to PAH might causally be involved in urothelial cancers. The human polymorphic GST catalysing conjugation of halomethanes, dihalomethanes, ethylene oxide and a number of other industrial compounds could be characterised as a class theta enzyme (GSTT1) by means of molecular biology. "Conjugator" and "non-conjugator" phenotypes are coincident with the presence and absence of the GSTT1 gene. There are wide variations in the frequencies of GSTT1 deletion (GSTT1 *0/0) among different ethnicities. Human phenotyping is facilitated by the GST activity towards methyl bromide or ethylene oxide in erythrocytes which is representative of the metabolic GSTT1 competence of the entire organism. Inter-individual variations in xenobiotic metabolism capacities may be due to polymorphisms of the genes coding for the enzymes themselves or of the genes coding for the receptors or transcription factors which regulate the expression of the enzymes. Also, polymorphisms in several regions of genes may cause altered ligand affinity, transactivation activity or expression levels of the receptor subsequently influencing the expression of the downstream target genes. Studies of individual susceptibility to toxicants and gene-environment interaction are now emerging as an important component of molecular epidemiology.

Genetic insights into common pathways and complex relationships among iImmune-mediated diseases

Shared aetiopathogenic factors among immune-mediated diseases have long been suggested by their co-familiality and co-occurrence, and molecular support has been provided by analysis of human leukocyte antigen (HLA) haplotypes and genome-wide association studies. The interrelationships can now be better appreciated following the genotyping of large immune disease sample sets on a shared SNP array: the 'Immunochip'. Here, we systematically analyse loci shared among major immune-mediated diseases. This reveals that several diseases share multiple susceptibility loci, but there are many nuances. The most associated variant at a given locus frequently differs and, even when shared, the same allele often has opposite associations. Interestingly, risk alleles conferring the largest effect sizes are usually disease-specific. These factors help to explain why early evidence of extensive 'sharing' is not always reflected in epidemiological overlap. © 2013 Macmillan Publishers Limited. All rights reserved.

Genetic control of bone density and turnover: Role of the collagen 1α1, estrogen receptor, and vitamin D receptor genes

Genetic factors are known to influence both the peak bone mass and probably the rate of change in bone density. A range of regulatory and structural genes has been proposed to be involved including collagen 1α1 (COL1A1), the estrogen receptor (ER), and the vitamin D receptor (VDR), but the actual genes involved are uncertain. We therefore studied the role of the COL1A1 and VDR loci in control of bone density by linkage in 45 dizygotic twin pairs and 29 nuclear families comprising 120 individuals. The influences on bone density of polymorphisms of COL1A1, VDR, and ER were studied by association both cross-sectionally and longitudinally in 193 elderly postmenopausal women (average age, 69 years) over a mean follow-up time of 6.3 years. Weak linkage of the COL1A1 locus with bone density was observed in both twins and families (p = 0.02 in both data sets), confirming previous observations of linkage of this locus with bone density. Association between the MscI polymorphism of COL1A1 and rate of lumbar spine bone loss was observed with significant gene-environment interaction related to dietary calcium intake (p = 0.0006). In the lowest tertile of dietary calcium intake, carriers of "s" alleles lost more bone than "SS" homozygotes (p = 0.01), whereas the opposite was observed in the highest dietary calcium intake (p = 0.003). Association also was observed between rate of bone loss at both the femoral neck and the lumbar spine and the TaqI VDR polymorphism (p = 0.03). This association was strongest in those in the lowest tertile of calcium intake, also suggesting the presence of gene-environment interaction involving dietary calcium and VDR, influencing bone turnover. No significant association was observed between the PvuII ER polymorphism alone or in combination with VDR or COL1A1 genotypes, with either bone density or its rate of change. These data support the involvement of COL1A1 in determination of bone density and the interaction of both COL1A1 and VDR with calcium intake in regulation of change of bone density over time.

Estimation in a multiplicative mixed model involving a genetic relationship matrix.

Genetic models partitioning additive and non-additive genetic effects for populations tested in replicated multi-environment trials (METs) in a plant breeding program have recently been presented in the literature. For these data, the variance model involves the direct product of a large numerator relationship matrix A, and a complex structure for the genotype by environment interaction effects, generally of a factor analytic (FA) form. With MET data, we expect a high correlation in genotype rankings between environments, leading to non-positive definite covariance matrices. Estimation methods for reduced rank models have been derived for the FA formulation with independent genotypes, and we employ these estimation methods for the more complex case involving the numerator relationship matrix. We examine the performance of differing genetic models for MET data with an embedded pedigree structure, and consider the magnitude of the non-additive variance. The capacity of existing software packages to fit these complex models is largely due to the use of the sparse matrix methodology and the average information algorithm. Here, we present an extension to the standard formulation necessary for estimation with a factor analytic structure across multiple environments.

Physiological Basis for Genotypic Variation in Tolerance to and Recovery from Pre-flowering Drought in Peanut.

Drought during the pre-flowering stage can increase yield of peanut. There is limited information on genotypic variation for tolerance to and recovery from pre-flowering drought (PFD) and more importantly the physiological traits underlying genotypic variation. The objectives of this study were to determine the effects of moisture stress during the pre-flowering phase on pod yield and to understand some of the physiological responses underlying genotypic variation in response to and recovery from PFD. A glasshouse and field experiments were conducted at Khon Kaen University, Thailand. The glasshouse experiment was a randomized complete block design consisting of two watering regimes, i.e. fully-irrigated control and 1/3 available soil water from emergence to 40 days after emergence followed by adequate water supply, and 12 peanut genotypes. The field experiment was a split-plot design with two watering regimes as main-plots, and 12 peanut genotypes as sub-plots. Measurements of N-2 fixation, leaf area (LA) were made in both experiments. In addition, root growth was measured in the glasshouse experiment. Imposition of PFD followed by recovery resulted in an average increase in yield of 24 % (range from 10 % to 57 %) and 12 % (range from 2 % to 51 %) in the field and glasshouse experiments, respectively. Significant genotypic variation for N-2 fixation, LA and root growth was also observed after recovery. The study revealed that recovery growth following release of PFD had a stronger influence on final yield than tolerance to water deficits during the PFD. A combination of N-2 fixation, LA and root growth accounted for a major portion of the genotypic variation in yield (r = 0.68-0.93) suggesting that these traits could be used as selection criteria for identifying genotypes with rapid recovery from PFD. A combined analysis of glasshouse and field experiments showed that LA and N-2 fixation during the recovery had low genotype x environment interaction indicating potential for using these traits for selecting genotypes in peanut improvement programs.

Comparative performance of Corymbia hybrids and parental species in subtropical Queensland and implications for breeding and deployment.

Eighty six full-sib Corymbia F1 hybrid families (crosses between C. torelliana and four spotted gum taxa: C. citriodora subsp. variegata, C. citriodora subsp. citriodora, C. henryi and C. maculata), were planted in six trials across six disparate sites in south-eastern Queensland to evaluate their productivity and determine their potential utility for plantation forestry. In each trial, the best-growing 20% of hybrid families grew significantly faster (P=0.05) than open-pollinated seedlots of the parent species Corymbia citriodora subsp. variegata, ranging from 107% to 181% and 127% to 287% of the height and diameter respectively. Relative performance of hybrid families growing on more than one site displayed consistency in ranking for growth across sites and analysis showed low genotype-by-environment interaction. Heritability estimates based on female and male parents across two sites at age six years for height and diameter at breast height, were high (0.62±0.28 to 0.64±0.35 and 0.31±0.21 to 0.69±0.37 respectively), and low to moderate (0.03±0.04 to 0.33±0.22) for stem straightness, branch size, incidence of ramicorns, and frost and disease resistance traits at ages one to three years. The proportion of dominance variance for height and diameter had reduced to zero by age six years. Based on these promising results, further breeding and pilot-scale family forestry and clonal forestry deployment is being undertaken. These results have also provided insights regarding the choice of a future hybrid breeding strategy.

Estimates of genetic parameters for growth and wood properties in Eucalyptus pellita F. Muell. to support tree breeding in Vietnam

Key message Eucalyptus pellita demonstrated good growth and wood quality traits in this study, with young plantation grown timber being suitable for both solid and pulp wood products. All traits examined were under moderate levels of genetic control with little genotype by environment interaction when grown on two contrasting sites in Vietnam. Context Eucalyptus pellita currently has a significant role in reforestation in the tropics. Research to support expanded of use of this species is needed: particularly, research to better understand the genetic control of key traits will facilitate the development of genetically improved planting stock. Aims This study aimed to provide estimates of the heritability of diameter at breast height over bark, wood basic density, Kraft pulp yield, modulus of elasticity and microfibril angle, and the genetic correlations among these traits, and understand the importance of genotype by environment interactions in Vietnam. Methods Data for diameter and wood properties were collected from two 10-year-old, open-pollinated progeny trials of E. pellita in Vietnam that evaluated 104 families from six native range and three orchard sources. Wood properties were estimated from wood samples using near-infrared (NIR) spectroscopy. Data were analysed using mixed linear models to estimate genetic parameters (heritability, proportion of variance between seed sources and genetic correlations). Results Variation among the nine sources was small compared to additive variance. Narrow-sense heritability and genetic correlation estimates indicated that simultaneous improvements in most traits could be achieved from selection among and within families as the genetic correlations among traits were either favourable or close to zero. Type B genetic correlations approached one for all traits suggesting that genotype by environment interactions were of little importance. These results support a breeding strategy utilizing a single breeding population advanced by selecting the best individuals across all seed sources. Conclusion Both growth and wood properties have been evaluated. Multi-trait selection for growth and wood property traits will lead to more productive populations of E. pellita both with improved productivity and improved timber and pulp properties.

Significant effects of Pgi genotype and body reserves on lifespan in the Glanville fritillary butterfly

Individuals with a particular variant of the gene phosphoglucose isomerase (Pgi ) have been shown to have superior dispersal capacity and fecundity in the Glanville fritillary butterfly (Melitaea cinxia), raising questions about the mechanisms that maintain polymorphism in this gene in the field. Here, we investigate how variation in the Pgi genotype affects female and male life history under controlled conditions. The most striking effect is the longer lifespan of genotypes with high dispersal capacity, especially in nonreproducing females. Butterflies use body reserves for somatic maintenance and reproduction, but
different resources (in thorax versus abdomen) are used under dissimilar conditions, with some interactions with the Pgi genotype. These results indicate life-history trade-offs that involve resource allocation and genotype!environment interactions, and these trade-offs are likely to contribute to the maintenance of Pgi polymorphism in the natural populations.

Complement Component 3: an assessment of association with AMD and analysis of gene-gene and gene-environment interactions in a Northern Irish cohort

Purpose: A non-synonymous single nucleotide polymorphism ( SNP) in complement component 3 has been shown to increase the risk of age-related macular degeneration (AMD). We assess its effect on AMD risk in a Northern Irish sample, test for gene-gene and gene-environment interaction, and review a risk prediction model.