Genetic heterogeneity of residual variance : estimation of variance components using double hierarchical generalized linear models

Background: The sensitivity to microenvironmental changes varies among animals and may be under genetic control. It is essential to take this element into account when aiming at breeding robust farm animals. Here, linear mixed models with genetic effects in the residual variance part of the model can be used. Such models have previously been fitted using EM and MCMC algorithms. Results: We propose the use of double hierarchical generalized linear models (DHGLM), where the squared residuals are assumed to be gamma distributed and the residual variance is fitted using a generalized linear model. The algorithm iterates between two sets of mixed model equations, one on the level of observations and one on the level of variances. The method was validated using simulations and also by re-analyzing a data set on pig litter size that was previously analyzed using a Bayesian approach. The pig litter size data contained 10,060 records from 4,149 sows. The DHGLM was implemented using the ASReml software and the algorithm converged within three minutes on a Linux server. The estimates were similar to those previously obtained using Bayesian methodology, especially the variance components in the residual variance part of the model. Conclusions: We have shown that variance components in the residual variance part of a linear mixed model can be estimated using a DHGLM approach. The method enables analyses of animal models with large numbers of observations. An important future development of the DHGLM methodology is to include the genetic correlation between the random effects in the mean and residual variance parts of the model as a parameter of the DHGLM.

Genetic and environmental heterogeneity of residual variance of weight traits in Nellore beef cattle

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

Genetic Control of Residual Variance for Teat Number in Pigs

The genetic improvement in litter size in pigs has been substantial during the last 10-15 years. The number of teats on the sow must increase as well to meet the needs of the piglets, because each piglet needs access to its own teat. We applied a genetic heterogeneity model on teat numberin sows, and estimated medium-high heritability for teat number (0.5), but low heritability for residual variance (0.05), indicating that selection for reduced variance might have very limited effect. A numerically positive correlation (0.8) between additive genetic breeding values for mean and for variance was found, but because of the low heritability for residual variance, the variance will increase very slowly with the mean.

Genetic variability of residual variance of production traits in Nellore beef cattle

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

Genetic heterogeneity of within-family variance of body weight in Atlantic salmon (Salmo salar)

BACKGROUND: Canalization is defined as the stability of a genotype against minor variations in both environment and genetics. Genetic variation in degree of canalization causes heterogeneity of within-family variance. The aims of this study are twofold: (1) quantify genetic heterogeneity of (within-family) residual variance in Atlantic salmon and (2) test whether the observed heterogeneity of (within-family) residual variance can be explained by simple scaling effects. RESULTS: Analysis of body weight in Atlantic salmon using a double hierarchical generalized linear model (DHGLM) revealed substantial heterogeneity of within-family variance. The 95% prediction interval for within-family variance ranged from ~0.4 to 1.2 kg2, implying that the within-family variance of the most extreme high families is expected to be approximately three times larger than the extreme low families. For cross-sectional data, DHGLM with an animal mean sub-model resulted in severe bias, while a corresponding sire-dam model was appropriate. Heterogeneity of variance was not sensitive to Box-Cox transformations of phenotypes, which implies that heterogeneity of variance exists beyond what would be expected from simple scaling effects. CONCLUSIONS: Substantial heterogeneity of within-family variance was found for body weight in Atlantic salmon. A tendency towards higher variance with higher means (scaling effects) was observed, but heterogeneity of within-family variance existed beyond what could be explained by simple scaling effects. For cross-sectional data, using the animal mean sub-model in the DHGLM resulted in biased estimates of variance components, which differed substantially both from a standard linear mean animal model and a sire-dam DHGLM model. Although genetic differences in canalization were observed, selection for increased canalization is difficult, because there is limited individual information for the variance sub-model, especially when based on cross-sectional data. Furthermore, potential macro-environmental changes (diet, climatic region, etc.) may make genetic heterogeneity of variance a less stable trait over time and space.

Dissecting the genetic heterogeneity of depression through age at onset.

Genome-wide studies in major depression have identified few replicated associations, potentially due to heterogeneity within the disorder. Several studies have suggested that age at onset (AAO) can distinguish sub-types of depression with specific heritable components. This paper investigates the role of AAO in the genetic susceptibility for depression using genome-wide association data on 2,746 cases and 1,594 screened controls from the RADIANT studies, with replication performed in 1,471 cases and 1,403 controls from two Munich studies. Three methods were used to analyze AAO: First a time-to-event analysis with controls censored, secondly comparing controls to case-subsets defined using AAO cut-offs, and lastly analyzing AAO as a quantitative trait. In the time-to-event analysis three SNPs reached suggestive significance (P < 5E-06), overlapping with the original case-control analysis of this study. In a case-control analysis using AAO thresholds, SNPs in 10 genomic regions showed suggestive association though again none reached genome-wide significance. Lastly, case-only analysis of AAO as a quantitative trait resulted in 5 SNPs reaching suggestive significance. Sex specific analysis was performed as a secondary analysis, yielding one SNP reaching genome-wide significance in early-onset males. No SNPs achieved significance in the replication study after correction for multiple testing. Analysis of AAO as a quantitative trait did suggest that, across all SNPs, common genetic variants explained a large proportion of the variance (51%, P = 0.04). This study provides the first focussed analysis of the genetic contribution to AAO in depression, and establishes a statistical framework that can be applied to a quantitative trait underlying any disorder. © 2012 Wiley Periodicals, Inc.

Molecular screening of ADAMTSL2 gene in 33 patients reveals the genetic heterogeneity of geleophysic dysplasia.

Background Geleophysic dysplasia (GD, OMIM 231050) is an autosomal recessive disorder characterised by short stature, small hands and feet, stiff joints, and thick skin. Patients often present with a progressive cardiac valvular disease which can lead to an early death. In a previous study including six GD families, we have mapped the disease gene on chromosome 9q34.2 and identified mutations in the A Disintegrin And Metalloproteinase with Thrombospondin repeats-like 2 gene (ADAMTSL2). Methods Following this study, we have collected the samples of 30 additional GD families, including 33 patients and identified ADAMTSL2 mutations in 14/33 patients, comprising 13 novel mutations. The absence of mutation in 19 patients prompted us to compare the two groups of GD patients, namely group 1, patients with ADAMTSL2 mutations (n=20, also including the 6 patients from our previous study), and group 2, patients without ADAMTSL2 mutations (n=19). Results The main discriminating features were facial dysmorphism and tip-toe walking, which were almost constantly observed in group 1. No differences were found concerning heart involvement, skin thickness, recurrent respiratory and ear infections, bronchopulmonary insufficiency, laryngo-tracheal stenosis, deafness, and radiographic features. Conclusions It is concluded that GD is a genetically heterogeneous condition. Ongoing studies will hopefully lead to the identification of another disease gene.

Genetic heterogeneity of the Coppock-like cataract: a mutation in CRYBB2 on chromosome 22q11.2.

PURPOSE: To identify the genetic defect for the Coppock-like cataract (CCL) affecting a Swiss family, which defect was unlinked to the chromosome 2q33-35 CCL locus. METHODS: A large family was characterized for linkage analysis by slit lamp examination or by the review of drawings made before cataract extraction. The affection status was attributed before genotyping, and the genotyping was masked to the affection status. Two-point and multipoint linkage analyses were performed using the MLINK and the LINKMAP components of the LINKAGE program package (ver. 5.1), respectively. Mutational analysis of candidate genes was performed by a combination of direct cycle sequencing and an amplification refractory mutation system assay. RESULTS: Ten individuals were affected with the CCL phenotype. The disease was autosomal dominant and appeared to be fully penetrant. A new CCL locus was identified on chromosome 22q11.2 within a 11.67-cM interval (maximum lod score [Zmax] = 4.14; theta = 0). Mutational analysis of the CRYBB2 candidate gene identified a disease-causing mutation in exon 6. This sequence change was identical with that previously described to be associated with the cerulean cataract, a clinically distinct entity. CONCLUSIONS: The CCL phenotype is genetically heterogeneous with a second gene on chromosome 22q11.2, CRYBB2. The CCL and the cerulean cataract are two distinct clinical entities associated with the same genetic defect. This work provides evidence for a modifier factor that influences cataract formation and that remains to be identified.

Genetic Heteroscedasticity for Domestic Animal Traits

Animal traits differ not only in mean, but also in variation around the mean. For instance, one sire’s daughter group may be very homogeneous, while another sire’s daughters are much more heterogeneous in performance. The difference in residual variance can partially be explained by genetic differences. Models for such genetic heterogeneity of environmental variance include genetic effects for the mean and residual variance, and a correlation between the genetic effects for the mean and residual variance to measure how the residual variance might vary with the mean. The aim of this thesis was to develop a method based on double hierarchical generalized linear models for estimating genetic heteroscedasticity, and to apply it on four traits in two domestic animal species; teat count and litter size in pigs, and milk production and somatic cell count in dairy cows. The method developed is fast and has been implemented in software that is widely used in animal breeding, which makes it convenient to use. It is based on an approximation of double hierarchical generalized linear models by normal distributions. When having repeated observations on individuals or genetic groups, the estimates were found to be unbiased. For the traits studied, the estimated heritability values for the mean and the residual variance, and the genetic coefficients of variation, were found in the usual ranges reported. The genetic correlation between mean and residual variance was estimated for the pig traits only, and was found to be favorable for litter size, but unfavorable for teat count.

Genetic heterogeneity of Escherichia coli isolated from pasteurized milk in State of Paraná, Brazil

Food contamination caused by enteric pathogens is a major cause of diarrheal disease worldwide, resulting in high morbidity and mortality and significant economic losses. Bacteria are important agents of foodborne diseases, particularly diarrheagenic Escherichia coli. The present study assessed the genetic diversity and antimicrobial resistance of E. coli isolates from pasteurized milk processed in 21 dairies in northwestern State of Parana, Brazil. The 95 E. coli isolates were subjected to antimicrobial susceptibility testing according to the recommendations of the Clinical and Laboratory Standards Institute and assessed genotypically by Enterobacterial Repetitive Intergenic Consensus-Polymerase Chain Reaction (ERIC-PCR). The highest rate of resistance was observed for cephalothin (55.78%). ERIC-PCR revealed high genetic diversity, clustering the 95 bacterial isolates into 90 different genotypic patterns. These results showed a heterogeneous population of E. coli in milk samples produced in the northwestern region of Parana and the need for good manufacturing practices throughout the processing of pasteurized milk to reduce the risk of foodborne illnesses.

Efeitos da Heterogeneidade de Variância Residual entre Grupos de Contemporâneos na Avaliação Genética de Bovinos de Corte

O objetivo deste estudo foi investigar, por meio de dados simulados, o efeito da heterogeneidade de variância residual entre grupos de contemporâneos (GC) sobre as avaliações genéticas de bovinos de corte, e comparar o uso de uma avaliação genética ponderada (R¹Isigmae²) em relação à avaliação que pressupõe homogeneidade de variância (R=Isigmae²). A característica estudada foi ganho de peso pós-desmame corrigido para 345 dias, sendo esta simulada com variância fenotípica de 300 kg² e herdabilidade igual a 0,4. A estrutura de um conjunto real de dados foi utilizada para fornecer os GC e os pais referentes às observações de cada animal. Cinco níveis de heterogeneidade de variância residual foram considerados de forma que os componentes de variância fossem, na média, iguais aos da situação de homogeneidade de variância. Na medida em que níveis mais acentuados de heterogeneidade de variância residual foram considerados, os animais foram selecionados dos GC com maior variabilidade, especialmente com pressão de seleção intensa. em relação à consistência de predição, os produtos e as vacas tiveram seus valores genéticos preditos mais afetados pela heterogeneidade de variância residual do que os touros. O fator de ponderação utilizado reduziu, mas não eliminou o efeito da heterogeneidade de variância. As avaliações genéticas ponderadas apresentaram resultados iguais ou superiores àqueles obtidos pelas avaliações que assumiram homogeneidade de variância. Mesmo quando não necessário, o uso de avaliações ponderadas produziu resultados não inferiores às avaliações que assumiram homogeneidade de variância.

Effect of pregnancy on the genetic evaluation of dairy cattle

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

Candidate gene linkage analysis indicates genetic heterogeneity in Marfan syndrome

Marfan syndrome (MFS) is an autosomal dominant disease of the connective tissue that affects the ocular, skeletal and cardiovascular systems, with a wide clinical variability. Although mutations in the FBN1 gene have been recognized as the cause of the disease, more recently other loci have been associated with MFS, indicating the genetic heterogeneity of this disease. We addressed the issue of genetic heterogeneity in MFS by performing linkage analysis of the FBN1 and TGFBR2 genes in 34 families (345 subjects) who met the clinical diagnostic criteria for the disease according to Ghent. Using a total of six microsatellite markers, we found that linkage with the FBN1 gene was observed or not excluded in 70.6% (24/34) of the families, and in 1 family the MFS phenotype segregated with the TGFBR2 gene. Moreover, in 4 families linkage with the FBN1 and TGFBR2 genes was excluded, and no mutations were identified in the coding region of TGFBR1, indicating the existence of other genes involved in MFS. Our results suggest that the genetic heterogeneity of MFS may be greater that previously reported.

Efeito da interação genótipo ambiente, sob a forma de heterogeneidade de variâncias entre rebanhos

Foram simuladas estruturas de dados em modelos mistos representando o teste de 100 reprodutores, sendo cada reprodutor acasalado com 10 matrizes (total de 1000 matrizes), originando em cada acasalamento 2 proles, totalizando 2000 proles (vinte proles por reprodutor). De cada combinação reprodutor e matriz, dez proles tiveram seu fenótipo expresso no ambiente de baixa produção (Estrato 1) e, a outra metade, no ambiente de alta produção (Estrato 2). A simulação foi realizada de forma a representar diferentes situações de presença de heterogeneidade de variâncias, combinando-se as origens da heterogeneidade, de natureza genética e ambiental. Na presença de heterogeneidade residual, o valor estimado para o componente de variância residual, considerando homogeneidade de variâncias se aproximou do valor médio das variâncias entre os estratos. Houve superestimação, também, do componente de variância genético aditivo. Ao simular heterogeneidade de variância de origem genética, observou-se que a estimação desse componente situou-se em valor intermediário aos simulados. Nessa situação, o componente de variância residual estimado foi próximo do valor simulado, indicando que a heterogeneidade de variâncias quando proveniente de fatores genéticos, não interfere, substancialmente, sobre e estimação do componente de variância residual. Na simulação de dados com presença de heterogeneidade tanto de origem genética quanto ambiental (estrutura de dados 4), conduziu a estimação de componentes de variâncias intermediários aos valores simulados em cada estrato. Assim, observa-se que, mesmo quando os reprodutores apresentam proles bem distribuídas em ambos os estratos, a heterogeneidade de variância proveniente de fatores não genético provoca distorções sobre a estimação da variância genética aditiva. Mas por outro lado, quando a heterogeneidade de variância é decorrente de fatores genéticos, não há grande interferência sobre a estimativa da variância residual, tal comportamento pode ser explicado pela incorporação da matriz de parentesco na estimação do componente de variância genético aditivo, possibilitando discriminar melhor a origem da diferenças entre variâncias. Na estrutura onde a variância residual foi heterogênea a estimativa de herdabilidade foi menor em relação à estrutura de homogeneidade de variâncias. Por outro lado, quando somente a variância genética aditiva foi heterogênea, a estimativa de herdabilidade, considerando-se apenas o estrato de alta variabilidade genética, foi inflacionada pela superestimação da variância genética aditiva. No entanto, a estimativa de herdabilidade obtida, desconsiderando essa fonte de heterogeneidade de variância, foi próxima à situação de homogeneidade de variância, indicando que, quando os reprodutores possuem boa distribuição de proles em diferentes ambientes, as estimativas relacionadas ao efeito genético são ponderadas pelo desempenho dos animais em cada ambiente. As correlações de Spearman e de Pearson entre os valores genéticos preditos dos reprodutores, para todas as situações, foram maiores que 0,90. O resultado indica que, mesmo havendo presença de heterogeneidade de variância genética e/ou ambiental, se os reprodutores possuem proles bem distribuídas entre os ambientes (estratos heterogêneos) a classificação do mérito genético não se altera, o que era esperado, pois em análises unicarácter, quando ocorre uma fonte de viés na avaliação genética, ela é comum a todos os indivíduos. Na situação em que foi imposta a estrutura de dados à presença de heterogeneidade de variância residual com número de número desigual de proles por reprodutor nos estratos, provocou superestimação dos componentes de variância. Porém mesmo havendo alteração na magnitude dos valores genéticos preditos para os reprodutores, a heterogeneidade de variância não alterou a classificação entre os reprodutores todas as correlações de ordem foram próximas à unidade. O efeito da heterogeneidade de variância, oriunda de fatores ambientais, ocasiona em maiores distorções sobre a avaliação genética animal, em relação, quando a mesma é proveniente de causas genéticas. A conexidade genética entre diferentes ambientes, dilui o efeito da heterogeneidade de variância, tanto de origem genética, quanto ambiental, na predição de valores genéticos dos reprodutores.

Genetic history of hepatitis C virus in Pakistan

Hepatitis C virus (HCV) genotype 3a accounts for similar to 80% of HCV infections in Pakistan, where similar to 10 million people are HCV-infected. Here, we report analysis of the genetic heterogeneity of HCV NS3 and NS5b subgenomic regions from genotype 3a variants obtained from Pakistan. Phylogenetic analyses showed that Pakistani genotype 3a variants were as genetically diverse as global variants, with extensive intermixing. Bayesian estimates showed that the most recent ancestor for genotype 3a in Pakistan was last extant in similar to 1896-1914 C.E. (range: 1851-1932). This genotype experienced a population expansion starting from similar to 1905 to similar to 1970 after which the effective population leveled. Death/birth models suggest that HCV 3a has reached saturating diversity with decreasing turnover rate and positive extinction. Taken together, these observations are consistent with a long and complex history of HCV 3a infection in Pakistan. Published by Elsevier B.V.