Fitting genetic models to twin data with binary and ordered categorical responses: A comparison of structural equation modelling and Bayesian hierarchical models

We compare Bayesian methodology utilizing free-ware BUGS (Bayesian Inference Using Gibbs Sampling) with the traditional structural equation modelling approach based on another free-ware package, Mx. Dichotomous and ordinal (three category) twin data were simulated according to different additive genetic and common environment models for phenotypic variation. Practical issues are discussed in using Gibbs sampling as implemented by BUGS to fit subject-specific Bayesian generalized linear models, where the components of variation may be estimated directly. The simulation study (based on 2000 twin pairs) indicated that there is a consistent advantage in using the Bayesian method to detect a correct model under certain specifications of additive genetics and common environmental effects. For binary data, both methods had difficulty in detecting the correct model when the additive genetic effect was low (between 10 and 20%) or of moderate range (between 20 and 40%). Furthermore, neither method could adequately detect a correct model that included a modest common environmental effect (20%) even when the additive genetic effect was large (50%). Power was significantly improved with ordinal data for most scenarios, except for the case of low heritability under a true ACE model. We illustrate and compare both methods using data from 1239 twin pairs over the age of 50 years, who were registered with the Australian National Health and Medical Research Council Twin Registry (ATR) and presented symptoms associated with osteoarthritis occurring in joints of the hand.

Genetic epidemiology of fecal egg excretion during Schistosoma mansoni infection in an endemic area in Minas Gerais, Brazil

There is considerable variation in the level of fecal egg excretion during Schistosoma mansoni infections. Within a single endemic area, the distribution of egg counts is typically overdispersed, with the majority of eggs excreted coming from a minority of residents. The purpose of this study was to quantify the influence of genetic factors on patterns of fecal egg excretion in a rural study sample in Brazil. Individual fecal egg excretions, expressed in eggs per gram of feces, were determined by the Kato-Katz method on stool samples collected on three different days. Detailed genealogic information was gathered at the time of sampling, which allowed assignment of 461 individuals to 14 pedigrees containing between 3 and 422 individuals. Using a maximum likelihood variance decomposition approach, we performed quantitative genetic analyses to determine if genetic factors could partially account for the observed pattern of fecal egg excretion. The quantitative genetic analysis indicated that between 21-37% of the variation in S. mansoni egg counts was attributable to additive genetic factors and that shared environment, as assessed by common household, accounted for a further 12-21% of the observed variation. A maximum likelihood heritability (h²) estimate of 0.44 ± 0.14 (mean ± SE) was found for the 9,604 second- and higher-degree pairwise relationships in the study sample, which is consistent with the upper limit (37%) of the genetic factor determined in the variance decomposition analysis. These analyses point to the significant influence of additive host genes on the pattern of S. mansoni fecal egg excretion in this endemic area.

Arterial properties in relation to genetic variations in the adducin subunits in a white population.

BACKGROUND: Adducin is a membrane skeleton protein, which consists of either alpha- and beta- or alpha- and gamma-subunits. We investigated whether arterial characteristics might be related to the genes encoding ADD1 (Gly460Trp-rs4961), ADD2 (C1797T-rs4984), and ADD3 (IVS11+386A>G-rs3731566). METHODS: We randomly recruited 1,126 Flemish subjects (mean age, 43.8 years; 50.3% women). Using a wall-tracking ultrasound system, we measured the properties of the carotid, femoral, and brachial arteries. We studied multivariate-adjusted phenotype-genotype associations, using a population- and family-based approach. RESULTS: In single-gene analyses, brachial diameter was 0.15 mm (P = 0.0022) larger, and brachial distensibility and cross-sectional compliance were 1.55 x 10(-3)/kPa (P = 0.013) and 0.017 mm(2)/kPa (P = 0.0029) lower in ADD3 AA than ADD3 GG homozygotes with an additive effect of the G allele. In multiple-gene analyses, the association of brachial diameter and distensibility with the ADD3 G allele occurred only in ADD1 GlyGly homozygotes. Otherwise, the associations between the arterial phenotypes in the three vascular beds and the ADD1 or ADD2 polymorphisms were not significant. In family-based analyses, the multivariate-adjusted heritability was 0.52, 0.38, and 0.30 for brachial diameter, distensibility, and cross-sectional compliance, respectively (P < 0.001). There was no evidence for population stratification (0.07 < or = P < or = 0.96). Transmission of the mutated ADD3 G allele was associated with smaller brachial diameter in 342 informative offspring (-0.12 +/- 0.04 mm; P = 0.0085) and in 209 offspring, who were ADD1 GlyGly homozygotes (-0.14 +/- 0.06 mm; P = 0.018). CONCLUSIONS: In ADD1 GlyGly homozygotes, the properties of the brachial artery are related to the ADD3 (A386G) polymorphism, but the underlying mechanism needs further clarification.

Genetic variability in a population of arbuscular mycorrhizal fungi causes variation in plant growth.

Different species of arbuscular mycorrhizal fungi (AMF) alter plant growth and affect plant coexistence and diversity. Effects of within-AMF species or within-population variation on plant growth have received less attention. High genetic variation exists within AMF populations. However, it is unknown whether genetic variation contributes to differences in plant growth. In our study, a population of AMF was cultivated under identical conditions for several generations prior to the experiments thus avoiding environmental maternal effects. We show that genetically different Glomus intraradices isolates from one AMF population significantly alter plant growth in an axenic system and in greenhouse experiments. Isolates increased or reduced plant growth meaning that plants potentially receive benefits or are subject to costs by forming associations with different individuals in the AMF population. This shows that genetic variability in AMF populations could affect host-plant fitness and should be considered in future research to understand these important soil organisms.

Genetic polymorphisms of the main transcription factors for adiponectin gene promoter in regulation of adiponectin levels: association analysis in three European cohorts.

Adiponectin serum concentrations are an important biomarker in cardiovascular epidemiology with heritability etimates of 30-70%. However, known genetic variants in the adiponectin gene locus (ADIPOQ) account for only 2%-8% of its variance. As transcription factors are thought to play an under-acknowledged role in carrying functional variants, we hypothesized that genetic polymorphisms in genes coding for the main transcription factors for the ADIPOQ promoter influence adiponectin levels. Single nucleotide polymorphisms (SNPs) at these genes were selected based on the haplotype block structure and previously published evidence to be associated with adiponectin levels. We performed association analyses of the 24 selected SNPs at forkhead box O1 (FOXO1), sterol-regulatory-element-binding transcription factor 1 (SREBF1), sirtuin 1 (SIRT1), peroxisome-proliferator-activated receptor gamma (PPARG) and transcription factor activating enhancer binding protein 2 beta (TFAP2B) gene loci with adiponectin levels in three different European cohorts: SAPHIR (n = 1742), KORA F3 (n = 1636) and CoLaus (n = 5355). In each study population, the association of SNPs with adiponectin levels on log-scale was tested using linear regression adjusted for age, sex and body mass index, applying both an additive and a recessive genetic model. A pooled effect size was obtained by meta-analysis assuming a fixed effects model. We applied a significance threshold of 0.0033 accounting for the multiple testing situation. A significant association was only found for variants within SREBF1 applying an additive genetic model (smallest p-value for rs1889018 on log(adiponectin) = 0.002, β on original scale = -0.217 µg/ml), explaining ∼0.4% of variation of adiponectin levels. Recessive genetic models or haplotype analyses of the FOXO1, SREBF1, SIRT1, TFAPB2B genes or sex-stratified analyses did not reveal additional information on the regulation of adiponectin levels. The role of genetic variations at the SREBF1 gene in regulating adiponectin needs further investigation by functional studies.

Phenotypic plasticity in salmonid embryos : characterizing pathogen-induced stress effects on heritable variation in traits and reaction norms

Les pressions écologiques peuvent varier tant en nature qu'en intensité dans le temps et l'espace. C'est pourquoi, un phénotype unique ne peut pas forcément conférer la meilleure valeur sélective. La plasticité phénotypique peut être un moyen de s'accommoder de cette situation, en augmentant globalement la tolérance aux changements environnementaux. Comme pour tout trait de caractère, une variation génétique doit persister pour qu'évoluent les traits plastiques dans une population donnée. Cependant, les pressions extérieures peuvent affecter l'héritabilité, et la direction de ces changements peut dépendre du caractère en question, de l'espèce mais aussi du type de stress. Dans la présente thèse, nous avons cherché à élucider les effets des pressions pathogéniques sur les phénotypes et la génétique quantitative de plusieurs traits plastiques chez les embryons de deux salmonidés, la palée (Coregonus palaea), et la truite de rivière (Salmo trutta). Les salmonidés se prêtent à de telles études du fait de leur extraordinaire variabilité morphologique, comportementale et des traits d'histoire de vie. Par ailleurs, avec le déclin des salmonidés dans le monde, il est important de savoir combien la variabilité génétique persiste dans les normes de réaction afin d'aider à prédire leur capacité à répondre aux changements de leur milieu. Nous avons observé qu'une augmentation de la croissance des communautés microbiennes symbiotiques entraînait une mortalité accrue et une éclosion précoce chez la palée, et dévoilait la variance génétique additive pour ces deux caractères (Chapitres 1-2). Bien qu'aucune variation génétique n'ait été trouvée pour les normes de réaction, nous avons observé une variabilité de la plasticité d'éclosion. Néanmoins, on a trouvé que les temps d'éclosion étaient corrélés entre les environnements, ce qui pourrait limiter l'évolution de la norme de réaction. Le temps d'éclosion des embryons est lié à la taille des géniteurs mâles, ce qui indique des effets pléiotropiques. Dans le Chapitre 3, nous avons montré qu'une interaction triple entre la souche bactérienne {Pseudomonas fluorescens}, l'état de dévelopement de l'hôte ainsi que ses gènes ont une influence sur la mortalité, le temps d'éclosion et la taille des alevins de la palée. Nous avons démontré qu'une variation génétique subsistait généralement dans les normes de réaction des temps d'éclosion, mais rarement pour la taille des alevins, et jamais pour la mortalité. Dans le même temps, nous avons exhibé que des corrélations entre environnements dépendaient des caractères phénotypiques, mais contrairement au Chapitre 2, nous n'avons pas trouvé de preuve de corrélations transgénérationnelles. Le Chapitre 4 complète le chapitre précédent, en se plaçant du point de vue moléculaire, et décrit comment le traitement d'embryons avec P. fluorescens s'est traduit par une régulation négative d'expression du CMH-I indépendemment de la souche bactérienne. Nous avons non seulement trouvé une variation génétique des caractères phénotypiques moyens, mais aussi de la plasticité. Les deux derniers chapitres traitent de l'investigation, chez la truite de rivière, des différences spécifiques entre populations pour des normes de réaction induites par les pathogènes. Dans le Chapitre 5, nous avons illustré que le métissage entre des populations génétiquement distinctes n'affectait en rien la hauteur ou la forme des normes de réaction d'un trait précoce d'histoire de vie suite au traitement pathogénique. De surcroît, en dépit de l'éclosion tardive et de la réduction de la taille des alevins, le traitement n'a pas modifié la variation héritable des traits de caractère. D'autre part, dans le Chapitre 6, nous avons démontré que le traitement d'embryons avec des stimuli contenus dans l'eau de conspécifiques infectés a entraîné des réponses propre à chaque population en terme de temps d'éclosion ; néanmoins, nous avons observé peu de variabilité génétique des normes de réaction pour ce temps d'éclosion au sein des populations. - Ecological stressors can vary in type and intensity over space and time, and as such, a single phenotype may not confer the highest fitness. Phenotypic plasticity can act as a means to accommodate this situation, increasing overall tolerance to environmental change. As with any trait, for plastic traits to evolve in a population, genetic variation must persist. However, environmental stress can alter trait heritability, and the direction of this shift can be trait, species, and stressor-dependent. In this thesis, we sought to understand the effects of pathogen stressors on the phenotypes and genetic architecture of several plastic traits in the embryos of two salmonids, the whitefish (Coregonus palaea), and the brown trout (Salmo trutta). Salmonids lend themselves to such studies because their extraordinary variability in morphological, behavioral, and life-history traits. Also, with declines in salmonids worldwide, knowing how much genetic variability persists in reaction norms may help predict their ability to respond to environmental change. We found that increasing growth of symbiotic microbial communities increased mortality and induced hatching in whitefish, and released additive genetic variance for both traits (Chapters 1-2). While no genetic variation was found for survival reaction norms, we did find variability in hatching plasticity. Nevertheless, hatching time was correlated across environments, which could constrain evolution of the reaction norm. Hatching time in the induced environment was also correlated to sire size, indicating pleiotropic effects. In Chapter 3 we report that a three-way interaction between bacterial strain (Pseudomonas fluorescens), host developmental stage, and host genetics impacted mortality, hatching time, and hatchling size in whitefish. We also showed that genetic variation generally persisted in hatching age reaction norms, but rarely for hatchling length, and never for mortality. At the same time, we demonstrated that cross-environmental correlations were trait-dependent, and unlike Chapter 2, we found no evidence of cross-generational correlations. Chapter 4 expands on the previous chapter, moving to the molecular level, and describes how treatment of embryos with P. fluorescens resulted in strain-independent downregulation of MHC class I. Genetic variation was evident not only in trait means, but also in plasticity. In the last two chapters, we investigated population level differences in pathogen- induced reaction norms in brown trout. In Chapter 5, we found that interbreeding between genetically distinct populations did not affect the elevation or shapes of the reaction norms of early life-history traits after pathogen challenge. Moreover, despite delaying hatching and reducing larval length, treatment produced no discernable shifts in heritable variation in traits. On the other hand, in Chapter 6, we found that treatment of embryos with water-borne cues from infected conspecifics elicited population-specific responses in terms of hatching time; however, we found little evidence of genetic variability in hatching reaction norms within populations. We have made considerable progress in understanding how pathogen stressors affect various early life-history traits in salmonid embryos. We have demonstrated that the effect of a particular stressor on heritable variation in these traits can vary according to the trait and species under consideration, in addition to the developmental stage of the host. Moreover, we found evidence of genetic variability in some, but not all reaction norms in whitefish and brown trout.

Additive effects of HLA alleles and innate immune genes determine viral outcome in HCV infection.

BACKGROUND: Chronic HCV infection is a leading cause of liver-related morbidity globally. The innate and adaptive immune responses are thought to be important in determining viral outcomes. Polymorphisms associated with the IFNL3 (IL28B) gene are strongly associated with spontaneous clearance and treatment outcomes. OBJECTIVE: This study investigates the importance of HLA genes in the context of genetic variation associated with the innate immune genes IFNL3 and KIR2DS3. DESIGN: We assess the collective influence of HLA and innate immune genes on viral outcomes in an Irish cohort of women (n=319) who had been infected from a single source as well as a more heterogeneous cohort (Swiss Cohort, n=461). In the Irish cohort, a number of HLA alleles are associated with different outcomes, and the impact of IFNL3-linked polymorphisms is profound. RESULTS: Logistic regression was performed on data from the Irish cohort, and indicates that the HLA-A*03 (OR 0.36 (0.15 to 0.89), p=0.027) -B*27 (OR 0.12 (0.03 to 0.45), p=<0.001), -DRB1*01:01 (OR 0.2 (0.07 to 0.61), p=0.005), -DRB1*04:01 (OR 0.31 (0.12 to 0.85, p=0.02) and the CC IFNL3 rs12979860 genotypes (OR 0.1 (0.04 to 0.23), p<0.001) are significantly associated with viral clearance. Furthermore, DQB1*02:01 (OR 4.2 (2.04 to 8.66), p=0.008), KIR2DS3 (OR 4.36 (1.62 to 11.74), p=0.004) and the rs12979860 IFNL3 'T' allele are associated with chronic infection. This study finds no interactive effect between IFNL3 and these Class I and II alleles in relation to viral clearance. There is a clear additive effect, however. Data from the Swiss cohort also confirms independent and additive effects of HLA Class I, II and IFNL3 genes in their prediction of viral outcome. CONCLUSIONS: This data supports a critical role for the adaptive immune response in the control of HCV in concert with the innate immune response.

Progressive erosion of genetic and epigenetic variation in callus-derived cocoa (Theobroma cacao) plants

Relatively little is known about the timing of genetic and epigenetic forms of somaclonal variation arising from callus growth. We surveyed for both types of change in cocoa (Theobroma cacao) plants regenerated from calli of various ages, and also between tissues from the source trees. For genetic change, we used 15 single sequence repeat (SSR) markers from four source trees and from 233 regenerated plants. For epigenetic change, we used 386 methylation-sensitive amplified polymorphism (MSAP) markers on leaf and explant (staminode) DNA from two source trees and on leaf DNA from 114 regenerants. Genetic variation within source trees was limited to one slippage mutation in one leaf. Regenerants were far more variable, with 35% exhibiting at least one mutation. Genetic variation initially accumulated with culture age but subsequently declined. MSAP (epigenetic) profiles diverged between leaf and staminode samples from source trees. Multivariate analysis revealed that leaves from regenerants occupied intermediate eigenspace between leaves and staminodes of source plants but became progressively more similar to source tree leaves with culture age. Statistical analysis confirmed this rather counterintuitive finding that leaves of ‘late regenerants’ exhibited significantly less genetic and epigenetic divergence from source leaves than those exposed to short periods of callus growth.

Genome-wide association mapping in a wild avian population identifies a link between genetic and phenotypic variation in a life-history trait

Understanding the genetic basis of traits involved in adaptation is a major challenge in evolutionary biology but remains poorly understood. Here, we use genome-wide association mapping using a custom 50 k single nucleotide polymorphism (SNP) array in a natural population of collared flycatchers to examine the genetic basis of clutch size, an important life-history trait in many animal species. We found evidence for an association on chromosome 18 where one SNP significant at the genome-wide level explained 3.9% of the phenotypic variance. We also detected two suggestive quantitative trait loci (QTLs) on chromosomes 9 and 26. Fitness differences among genotypes were generally weak and not significant, although there was some indication of a sex-by-genotype interaction for lifetime reproductive success at the suggestive QTL on chromosome 26. This implies that sexual antagonism may play a role in maintaining genetic variation at this QTL. Our findings provide candidate regions for a classic avian life-history trait that will be useful for future studies examining the molecular and cellular function of, as well as evolutionary mechanisms operating at, these loci.

Genetic parameter estimates for buffalo milk yield, milk quality and mozzarella production and Bayesian inference analysis of their relationships

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

Genetic parameters for buffalo milk yield and milk quality traits using Bayesian inference

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

Genetic association between body composition measured by ultrasound and visual scores in Brazilian Nelore cattle

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

Genetic parameter estimates for live weight and daily live weight gain obtained for Nellore bulls in a test station using different models

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

Quantitative variation in the tropical maize population, ESALQ-PB1

A boa produtividade e os valores intermediários para altura da planta e altura da espiga caracterizam a população de milho ESALQ-PB1 como agronomicamente promissora. São relatadas estimativas de parâmetros para 13 caracteres: altura da planta (PH), altura da espiga (EH), posição relativa da espiga (EP), comprimento do pendão (TL), peso do pendão (TW), número de ramificações do pendão (TB), peso de espigas (EW), peso de grãos (GW), comprimento da espiga (EL), diâmetro da espiga (ED), número de fileiras de grãos (RN), número de grãos por fileira (KR) e prolificidade (PR). Os resultados se referem a um único ambiente (um local e um ano). Foi detectada variação genética para todos os caracteres, e são apresentadas estimativas da variância genética aditiva. Os coeficientes de herdabilidade (indivíduos) variaram de 0,14 a 0,72 e foram considerados altos para PH, EH e TB; intermediários para EP, TL, TW, EL e ED, e baixos para EW, GW, KR e PR. O coeficiente de herdabilidade para médias de progênies mostrou aproximadamente a mesma tendência, variando de 0,40 a 0,75. O maior ganho esperado por seleção foi para TB (27% por ciclo) sob seleção massal e para TW (16,4%) por seleção entre progênies; o menor ganho esperado foi para ED, tanto por seleção massal (1,9%) como por seleção entre progênies (2,9%). Coeficientes de correlação aditiva (rA) 0.5

Estimation of (co)variance components and genetic parameters for weights of red-winged tinamou using random regression models

Com este trabalho objetivou-se determinar parâmetros genéticos para peso corporal de perdizes em cativeiro. Foram utilizados modelos de regressão aleatória na análise dos dados considerando os efeitos genéticos aditivos diretos (AD) e de ambiente permanente de animal (AP) como aleatórios. As variâncias residuais foram modeladas utilizando-se funções de variância de ordem 5. A curva média da população foi ajustada por polinômios ortogonais de Legendre de ordem 6. Os efeitos genéticos aditivos diretos e de ambiente permanente de animal foram modelados utilizando-se polinômios de Legendre de segunda a nona ordem. Os melhores resultados foram obtidos pelos modelos de ordem 6 de ajuste para os efeitos genéticos aditivos diretos e de ordem 3 para os de ambiente permanente pelo Critério de Informação de Akaike e ordem 3 para ambos os efeitos pelos Critério de Informação Bayesiano de Schwartz e Teste de Razão de Verossimilhança. As herdabilidades estimadas variaram de 0,02 a 0,57. O primeiro autovalor respondeu por 94 e 90% da variação decorrente de efeitos aditivos diretos e de ambiente permanente, respectivamente. A seleção de perdizes para peso corporal é mais efetiva a partir de 112 dias de idade.