The Importance of Scale for Spatial-confounding Bias and Precision of Spatial Regression Estimators

Increasingly, regression models are used when residuals are spatially correlated. Prominent examples include studies in environmental epidemiology to understand the chronic health effects of pollutants. I consider the effects of residual spatial structure on the bias and precision of regression coefficients, developing a simple framework in which to understand the key issues and derive informative analytic results. When the spatial residual is induced by an unmeasured confounder, regression models with spatial random effects and closely-related models such as kriging and penalized splines are biased, even when the residual variance components are known. Analytic and simulation results show how the bias depends on the spatial scales of the covariate and the residual; bias is reduced only when there is variation in the covariate at a scale smaller than the scale of the unmeasured confounding. I also discuss how the scales of the residual and the covariate affect efficiency and uncertainty estimation when the residuals can be considered independent of the covariate. In an application on the association between black carbon particulate matter air pollution and birth weight, controlling for large-scale spatial variation appears to reduce bias from unmeasured confounders, while increasing uncertainty in the estimated pollution effect.

Variability between laboratories performing coagulation tests with identical platforms: a nationwide evaluation study.

BACKGROUND While the assessment of analytical precision within medical laboratories has received much attention in scientific enquiry, the degree of as well as the sources causing variation between them remains incompletely understood. In this study, we quantified the variance components when performing coagulation tests with identical analytical platforms in different laboratories and computed intraclass correlations coefficients (ICC) for each coagulation test. METHODS Data from eight laboratories measuring fibrinogen twice in twenty healthy subjects with one out of 3 different platforms and single measurements of prothrombin time (PT), and coagulation factors II, V, VII, VIII, IX, X, XI and XIII were analysed. By platform, the variance components of (i) the subjects, (ii) the laboratory and the technician and (iii) the total variance were obtained for fibrinogen as well as (i) and (iii) for the remaining factors using ANOVA. RESULTS The variability for fibrinogen measurements within a laboratory ranged from 0.02 to 0.04, the variability between laboratories ranged from 0.006 to 0.097. The ICC for fibrinogen ranged from 0.37 to 0.66 and from 0.19 to 0.80 for PT between the platforms. For the remaining factors the ICC's ranged from 0.04 (FII) to 0.93 (FVIII). CONCLUSIONS Variance components that could be attributed to technicians or laboratory procedures were substantial, led to disappointingly low intraclass correlation coefficients for several factors and were pronounced for some of the platforms. Our findings call for sustained efforts to raise the level of standardization of structures and procedures involved in the quantification of coagulation factors.

Reference Values and Factors Associated With Renal Resistive Index in a Family-Based Population Study

Increased renal resistive index (RRI) has been recently associated with target organ damage and cardiovascular or renal outcomes in patients with hypertension and diabetes mellitus. However, reference values in the general population and information on familial aggregation are largely lacking. We determined the distribution of RRI, associated factors, and heritability in a population-based study. Families of European ancestry were randomly selected in 3 Swiss cities. Anthropometric parameters and cardiovascular risk factors were assessed. A renal Doppler ultrasound was performed, and RRI was measured in 3 segmental arteries of both kidneys. We used multilevel linear regression analysis to explore the factors associated with RRI, adjusting for center and family relationships. Sex-specific reference values for RRI were generated according to age. Heritability was estimated by variance components using the ASSOC program (SAGE software). Four hundred women (mean age±SD, 44.9±16.7 years) and 326 men (42.1±16.8 years) with normal renal ultrasound had mean RRI of 0.64±0.05 and 0.62±0.05, respectively (P<0.001). In multivariable analyses, RRI was positively associated with female sex, age, systolic blood pressure, and body mass index. We observed an inverse correlation with diastolic blood pressure and heart rate. Age had a nonlinear association with RRI. We found no independent association of RRI with diabetes mellitus, hypertension treatment, smoking, cholesterol levels, or estimated glomerular filtration rate. The adjusted heritability estimate was 42±8% (P<0.001). In a population-based sample with normal renal ultrasound, RRI normal values depend on sex, age, blood pressure, heart rate, and body mass index. The significant heritability of RRI suggests that genes influence this phenotype.

Diffusion-weighted MR imaging of upper abdominal organs: field strength and intervendor variability of apparent diffusion coefficients

PURPOSE To determine the variability of apparent diffusion coefficient (ADC) values in various anatomic regions in the upper abdomen measured with magnetic resonance (MR) systems from different vendors and with different field strengths. MATERIALS AND METHODS Ten healthy men (mean age, 36.6 years ± 7.7 [standard deviation]) gave written informed consent to participate in this prospective ethics committee-approved study. Diffusion-weighted (DW) MR imaging was performed in each subject with 1.5- and 3.0-T MR systems from each of three vendors at two institutions. Two readers independently measured ADC values in seven upper abdominal regions (left and right liver lobe, gallbladder, pancreas, spleen, and renal cortex and medulla). ADC values were tested for interobserver differences, as well as for differences related to field strength and vendor, with repeated-measures analysis of variance; coefficients of variation (CVs) and variance components were calculated. RESULTS Interreader agreement was excellent (intraclass coefficient, 0.876). ADC values were (77.5-88.8) ×10(-5) mm(2)/sec in the spleen and (250.6-278.5) ×10(-5) mm(2)/sec in the gallbladder. There were no significant differences between ADC values measured at 1.5 T and those measured at 3.0 T in any anatomic region (P >.10 for all). In two of seven regions at 1.5 T (left and right liver lobes, P < .023) and in four of seven regions at 3.0 T (left liver lobe, pancreas, and renal cortex and medulla, P < .008), intervendor differences were significant. CVs ranged from 7.0% to 27.1% depending on the anatomic location. CONCLUSION Despite significant intervendor differences in ADC values of various anatomic regions of the upper abdomen, ADC values of the gallbladder, pancreas, spleen, and kidney may be comparable between MR systems from different vendors and between different field strengths.

Heritabilities of health traits in Swiss Warmblood horses

REASONS FOR PERFORMING THE STUDY: There is a lack of evidence regarding genetic parameters of health traits in Swiss Warmblood horses. OBJECTIVES: To estimate heritabilities of equine sarcoid disease, horn quality of the hooves, prognathism and increased filling of talocrural joints as a possible indicator for osteochondrosis in Swiss Warmblood horses examined at the field tests for 3-year-olds between 2005 and 2013. STUDY DESIGN: Retrospective analysis of breed society database. METHODS: Swiss Warmblood horses were examined clinically by 13 veterinarians at field tests in Switzerland between 2005 and 2013. The presence of sarcoids, horn quality of the hooves, incisor occlusion and increased joint filling were assessed and recorded. Records of 3715 horses were integrated in a pedigree comprising 217,282 horses. Variance components and heritabilities were estimated on the liability scale using MTGSAM. RESULTS: The prevalences of the examined traits were rather low, ranging from 2.4 to 13.0%. Heritabilities estimated were 0.21 ± 0.07 for the occurrence of sarcoids, 0.04 ± 0.02 for hooves with markedly brittle and friable horn quality, 0.03 ± 0.01 for hooves with marked growth ring formation, 0.06 ± 0.03 for prognathism and 0.08 ± 0.04 for increased filling of the talocrural joint (an indicator of possible osteochondrosis). The influence of the examiner on the variance of these observations was considerable. CONCLUSIONS: With the exception of equine sarcoid disease, estimates for the heritabilities for the traits examined here were low. A standardised examination protocol may reduce the variance due to the examiner. This article is protected by copyright. All rights reserved. This article is protected by copyright. All rights reserved.

THE NATURAL AND ORTHOGONAL INTERACTION (NOIA) MODELS FOR QUANTITATIVE TRAITS (QTs) AND COMPLEX DISEASES

My dissertation focuses on developing methods for gene-gene/environment interactions and imprinting effect detections for human complex diseases and quantitative traits. It includes three sections: (1) generalizing the Natural and Orthogonal interaction (NOIA) model for the coding technique originally developed for gene-gene (GxG) interaction and also to reduced models; (2) developing a novel statistical approach that allows for modeling gene-environment (GxE) interactions influencing disease risk, and (3) developing a statistical approach for modeling genetic variants displaying parent-of-origin effects (POEs), such as imprinting. In the past decade, genetic researchers have identified a large number of causal variants for human genetic diseases and traits by single-locus analysis, and interaction has now become a hot topic in the effort to search for the complex network between multiple genes or environmental exposures contributing to the outcome. Epistasis, also known as gene-gene interaction is the departure from additive genetic effects from several genes to a trait, which means that the same alleles of one gene could display different genetic effects under different genetic backgrounds. In this study, we propose to implement the NOIA model for association studies along with interaction for human complex traits and diseases. We compare the performance of the new statistical models we developed and the usual functional model by both simulation study and real data analysis. Both simulation and real data analysis revealed higher power of the NOIA GxG interaction model for detecting both main genetic effects and interaction effects. Through application on a melanoma dataset, we confirmed the previously identified significant regions for melanoma risk at 15q13.1, 16q24.3 and 9p21.3. We also identified potential interactions with these significant regions that contribute to melanoma risk. Based on the NOIA model, we developed a novel statistical approach that allows us to model effects from a genetic factor and binary environmental exposure that are jointly influencing disease risk. Both simulation and real data analyses revealed higher power of the NOIA model for detecting both main genetic effects and interaction effects for both quantitative and binary traits. We also found that estimates of the parameters from logistic regression for binary traits are no longer statistically uncorrelated under the alternative model when there is an association. Applying our novel approach to a lung cancer dataset, we confirmed four SNPs in 5p15 and 15q25 region to be significantly associated with lung cancer risk in Caucasians population: rs2736100, rs402710, rs16969968 and rs8034191. We also validated that rs16969968 and rs8034191 in 15q25 region are significantly interacting with smoking in Caucasian population. Our approach identified the potential interactions of SNP rs2256543 in 6p21 with smoking on contributing to lung cancer risk. Genetic imprinting is the most well-known cause for parent-of-origin effect (POE) whereby a gene is differentially expressed depending on the parental origin of the same alleles. Genetic imprinting affects several human disorders, including diabetes, breast cancer, alcoholism, and obesity. This phenomenon has been shown to be important for normal embryonic development in mammals. Traditional association approaches ignore this important genetic phenomenon. In this study, we propose a NOIA framework for a single locus association study that estimates both main allelic effects and POEs. We develop statistical (Stat-POE) and functional (Func-POE) models, and demonstrate conditions for orthogonality of the Stat-POE model. We conducted simulations for both quantitative and qualitative traits to evaluate the performance of the statistical and functional models with different levels of POEs. Our results showed that the newly proposed Stat-POE model, which ensures orthogonality of variance components if Hardy-Weinberg Equilibrium (HWE) or equal minor and major allele frequencies is satisfied, had greater power for detecting the main allelic additive effect than a Func-POE model, which codes according to allelic substitutions, for both quantitative and qualitative traits. The power for detecting the POE was the same for the Stat-POE and Func-POE models under HWE for quantitative traits.

Incorporação de informações genômicas no desenvolvimento de índices econômicos para a seleção de bovinos leiteiros

A eficiência econômica da bovinocultura leiteira está relacionada à utilização de animais que apresentem, concomitantemente, bom desempenho quanto à produção, reprodução, saúde e longevidade. Nisto, o índice de seleção configura-se como ferramenta importante ao aumento da lucratividade nesse sistema, visto que permite a seleção de reprodutores para várias características simultaneamente, considerando a relação entre elas bem como a relevância econômica das mesmas. Com a recente disponibilidade de dados genômicos tornou-se ainda possível expandir a abrangência e acurácia dos índices de seleção por meio do aumento do número e qualidade das informações consideradas. Nesse contexto, dois estudos foram desenvolvidos. No primeiro, o objetivo foi estimar parâmetros genéticos e valores genéticos (VG) para características relacionadas à produção e qualidade do leite incluindo-se a informação genômica na avaliação genética. Foram utilizadas medidas de idade ao primeiro parto (IPP), produção de leite (PROD), teor de gordura (GOR), proteína (PROT), lactose, caseína, escore de células somáticas (ECS) e perfil de ácidos graxos de 4.218 vacas bem como os genótipos de 755 vacas para 57.368 polimorfismos de nucleotídeo único (SNP). Os componentes de variância e VG foram obtidos por meio de um modelo misto animal, incluindo-se os efeitos de grupos de contemporâneas, ordem de lactação, dias em lactação e os efeitos aditivo genético, ambiente permanente e residual. Duas abordagens foram desenvolvidas: uma tradicional, na qual a matriz de relacionamentos é baseada no pedigree; e uma genômica, na qual esta matriz é construída combinando-se a informação de pedigree e dos SNP. As herdabilidades variaram de 0,07 a 0,39. As correlações genéticas entre PROD e os componentes do leite variaram entre -0,45 e -0,13 enquanto correlações altas e positivas foram estimadas entre GOR e os ácidos graxos. O uso da abordagem genômica não alterou as estimativas de parâmetros genéticos; contudo, houve aumento entre 1,5% e 6,8% na acurácia dos VG, à exceção de IPP, para a qual houve uma redução de 1,9%. No segundo estudo, o objetivo foi incorporar a informação genômica no desenvolvimento de índices econômicos de seleção. Neste, os VG para PROD, GOR, PROT, teor de ácidos graxos insaturados (INSAT), ECS e vida produtiva foram combinados em índices de seleção ponderados por valores econômicos estimados sob três cenários de pagamento: exclusivamente por volume de leite (PAG1); por volume e por componentes do leite (PAG2); por volume e componentes do leite incluindo INSAT (PAG3). Esses VG foram preditos a partir de fenótipos de 4.293 vacas e genótipos de 755 animais em um modelo multi-característica sob as abordagens tradicional e genômica. O uso da informação genômica influenciou os componentes de variância, VG e a resposta à seleção. Entretanto, as correlações de ranking entre as abordagens foram altas nos três cenários, com valores entre 0,91 e 0,99. Diferenças foram principalmente observadas entre PAG1 e os demais cenários, com correlações entre 0,67 e 0,88. A importância relativa das características e o perfil dos melhores animais foram sensíveis ao cenário de remuneração considerado. Assim, verificou-se como essencial a consideração dos valores econômicos das características na avaliação genética e decisões de seleção.

Gasterosteus aculeatus egg traits and fish lengths vs. experiment setup

1) Our study addresses the role of non-genetic and genetic inheritance in shaping the adaptive potential of populations under a warming ocean scenario. We used a combined experimental approach (transgenerational plasticity and quantitative genetics) to partition the relative contribution of maternal vs. paternal (additive genetic) effects to offspring body size (a key component of fitness), and investigated a potential physiological mechanism (mitochondrial respiration capacities) underlying whole organism growth/size responses. 2) In very early stages of growth (up to 30 days), offspring body size of marine sticklebacks benefited from maternal transgenerational plasticity (TGP): offspring of mothers acclimated to17°C were larger when reared at 17°C, and offspring of mothers acclimated to 21°C were larger when reared at 21°C. The benefits of maternal TGP on body size were stronger and persisted longer (up to 60 days) for offspring reared in the warmer (21°C) environment, suggesting that maternal effects will be highly relevant for climate change scenarios in this system. 3) Mitochondrial respiration capacities measured on mature offspring (F1 adults) matched the pattern of TGP for juvenile body size, providing an intuitive mechanistic basis for the maternal acclimation persisting into adulthood. Size differences between temperatures seen at early growth stages remained in the F1 adults, linking offspring body size to maternal inheritance of mitochondria. 4) Lower maternal variance components in the warmer environment were mostly driven by mothers acclimated to ambient (colder) conditions, further supporting our tenet that maternal effects were stronger at elevated temperature. Importantly, all parent-offspring temperature combination groups showed genotype x environment (GxE) interactions, suggesting that reaction norms have the potential to evolve. 5) To summarise, transgenerational plasticity and genotype x environment interactions work in concert to mediate impacts of ocean warming on metabolic capacity and early growth of marine sticklebacks. TGP can buffer short-term detrimental effects of climate warming and may buy time for genetic adaptation to catch up, therefore markedly contributing to the evolutionary potential and persistence of populations under climate change.

A major quantitative trait locus for CD4-CD8 ratio is located on chromosome 11

CD4-CD8 ratio is an important diagnostic measure of immune system functioning. In particular, CD4-CD8 ratio predicts the time taken for progression of HIV infection to acquired immune deficiency syndrome (AIDS) and the long-term survival of AIDS patients. To map genes that regulate differences between healthy individuals in CD4-CD8 ratio, we typed 757 highly polymorphic microsatellite markers at an average spacing of similar to5 cM across the genome in 405 pairs of dizygotic twins at ages 12, 14 and 16. We used multipoint variance components linkage analysis to test for linkage between marker loci and CD4-CD8 ratio at each age. We found suggestive evidence of linkage on chromosome 11p in 12-year-old twins (LOD=2.55, P=0.00031) and even stronger evidence of linkage in the same region at age 14 (LOD 3.51, P=0.00003). Possible candidate genes include CD5 and CD6, which encode cell membrane proteins involved in the positive selection of thymocytes. We also found suggestive evidence of linkage at other areas of the genome including regions on chromosomes 1, 3, 4, 5, 6, 12, 13, 15, 17 and 22.

Multivariate QTL linkage analysis suggests a QTL for platelet count on chromosome 19q

Platelet count is a highly heritable trait with genetic factors responsible for around 80% of the phenotypic variance. We measured platelet count longitudinally in 327 monozygotic and 418 dizygotic twin pairs at 12, 14 and 16 years of age. We also performed a genome-wide linkage scan of these twins and their families in an attempt to localize QTLs that influenced variation in platelet concentrations. Suggestive linkage was observed on chromosome 19q13.13-19q13.31 at 12 (LOD=2.12, P=0.0009), 14 (LOD=2.23, P=0.0007) and 16 (LOD=1.01, P=0.016) years of age and multivariate analysis of counts at all three ages increased the LOD to 2.59 (P=0.0003). A possible candidate in this region is the gene for glycoprotein VI, a receptor involved in platelet aggregation. Smaller linkage peaks were also seen at 2p, 5p, 5q, 10p and 15q. There was little evidence for linkage to the chromosomal regions containing the genes for thrombopoietin (3q27) and the thrombopoietin receptor (1q34), suggesting that polymorphisms in these genes do not contribute substantially to variation in platelet count between healthy individuals.

A genome scan for eye color in 502 twin families: Most variation is due to a QTL on chromosome 15q

We have rated eye color on a 3-point scale (1=blue/grey, 2=hazel/green, 3=brown) in 502 twin families and carried out a 5-10 cM genome scan (400-757 markers). We analyzed eye color as a threshold trait and performed multipoint sib pair linkage analysis using variance components analysis in Mx. A lod of 19.2 was found at the marker D15S1002, less than 1 cM from OCA2, which has been previously implicated in eye color variation. We estimate that 74% of variance in eye color liability is due to this QTL and a further 18% due to polygenic effects. However, a large shoulder on this peak suggests that other loci affecting eye color may be telomeric of OCA2 and inflating the QTL estimate. No other peaks reached genome-wide significance, although lods >2 were seen on 5p and 14q and lods >1 were additionally seen on chromosomes 2, 3, 6, 7, 8, 9, 17 and 18. Most of these secondary peaks were reduced or eliminated when we repeated the scan as a two locus analysis with the 15q linkage included, although this does not necessarily exclude them as false positives. We also estimated the interaction between the 15q QTL and the other marker locus but there was only minor evidence for additive x additive epistasis. Elaborating the analysis to the full two-locus model including non-additive main effects and interactions did not strengthen the evidence for epistasis. We conclude that most variation in eye color in Europeans is due to polymorphism in OCA2 but that there may be modifiers at several other loci.

Exploring the Etiology of the Association Between Birthweight and IQ in an Adolescent Twin Sample

The negative effects of very low birthweight on intellectual development have been well documented, and more recently this effect has been shown to generalise to birthweights within the normal range. In this study we investigate the etiology of this relationship by using a classical twin design to disentangle the contributions of genes and environment. A previous Dutch study (Boomsma et al., 2001) examining these effects indicated that genes were important in mediating the association of birthweight to full IQ measured at ages 7 and 10, but not at ages 5 and 12. Here the association between birthweight and IQ at age 16 is considered (N = 523 twin pairs). Using variance components modeling we found that the genetic variance in birthweight (4%) completely overlapped with that in verbal IQ but not performance or full IQ. Results further showed the importance of shared environmental effects on birthweight (similar to 60%) but not on IQ (with genes explaining up to 72% of IQ variance). Models incorporating a direction of causation parameter between birthweight and IQ provided adequate fit to the data in either causal direction for performance and full IQ, but the model with verbal 10 causing birthweight was preferred to one in which birthweight influenced verbal IQ. As the measurement of birthweight precedes the measurement of twins' IQ at age 16, the influence of verbal IQ might be better considered as a proxy for parents' 10 or education, and it is possible that brighter mothers provide better prenatal environments for their children.

Copulas in QTL mapping

The standard variance components method for mapping quantitative trait loci is derived on the assumption of normality. Unsurprisingly, statistical tests based on this method do not perform so well if this assumption is not satisfied. We use the statistical concept of copulas to relax the assumption of normality and derive a test that can perform well under any distribution of the continuous trait. In particular, we discuss bivariate normal copulas in the context of sib-pair studies. Our approach is illustrated by a linkage analysis of lipoprotein(a) levels, whose distribution is highly skewed. We demonstrate that the asymptotic critical levels of the test can still be calculated using the interval mapping approach. The new method can be extended to more general pedigrees and multivariate phenotypes in a similar way as the original variance components method.

Major quantitative trait locus for eosinophil count is located on chromosome 2q

Background: Eosinophils are granulocytic white blood cells implicated in asthma and atopic disease. The degree of eosinophilia in the blood of patients with asthma correlates with the severity of asthmatic symptoms. Quantitative trait loci (QTL) linkage analysis of eosinophil count may be a more powerful strategy of mapping genes involved in asthma than linkage analysis using affected relative pairs. 1 Objective: To identify QTLs responsible for variation in eosinophil count in adolescent twins. Methods: We measured eosinophil count longitudinally in 738 pairs of twins at 12, 14, and 16 years of age. We typed 757 highly polymorphic microsatellite markers at an average spacing of similar to5 centimorgans across the genome. We then used multipoint variance components linkage analysis to test for linkage between marker loci and eosinophil concentrations at each age across the genome. Results: We found highly significant linkage on chromosome 2q33 in 12-year-old twins (logarithm of the odds = 4.6; P = .000002) and suggestive evidence of linkage in the same region in 14-year-olds (logarithm of the odds = 1.0; P = .016). We also found suggestive evidence of linkage at other areas of the genome, including regions on chromosomes 2, 3, 4, 8, 9, 11, 12, 17, 20, and 22. Conclusion: A QTL for eosinophil count is present on chromosome 2q33. This QTL might represent a gene involved in asthma pathophysiology.

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.