On estimation and identifiability issues of sex-linked inheritance with a case study of pigmentation in Swiss barn owl (Tyto alba)

Genetic evaluation using animal models or pedigree-based models generally assume only autosomal inheritance. Bayesian animal models provide a flexible framework for genetic evaluation, and we show how the model readily can accommodate situations where the trait of interest is influenced by both autosomal and sex-linked inheritance. This allows for simultaneous calculation of autosomal and sex-chromosomal additive genetic effects. Inferences were performed using integrated nested Laplace approximations (INLA), a nonsampling-based Bayesian inference methodology. We provide a detailed description of how to calculate the inverse of the X- or Z-chromosomal additive genetic relationship matrix, needed for inference. The case study of eumelanic spot diameter in a Swiss barn owl (Tyto alba) population shows that this trait is substantially influenced by variation in genes on the Z-chromosome (sigma(2)(z) = 0.2719 and sigma(2)(a) = 0.4405). Further, a simulation study for this study system shows that the animal model accounting for both autosomal and sex-chromosome-linked inheritance is identifiable, that is, the two effects can be distinguished, and provides accurate inference on the variance components.

Total hemoglobin mass--a new parameter to detect blood doping?

PURPOSE: All kinds of blood manipulations aim to increase the total hemoglobin mass (tHb-mass). To establish tHb-mass as an effective screening parameter for detecting blood doping, the knowledge of its normal variation over time is necessary. The aim of the present study, therefore, was to determine the intraindividual variance of tHb-mass in elite athletes during a training year emphasizing off, training, and race seasons at sea level. METHODS: tHb-mass and hemoglobin concentration ([Hb]) were determined in 24 endurance athletes five times during a year and were compared with a control group (n = 6). An analysis of covariance was used to test the effects of training phases, age, gender, competition level, body mass, and training volume. Three error models, based on 1) a total percentage error of measurement, 2) the combination of a typical percentage error (TE) of analytical origin with an absolute SD of biological origin, and 3) between-subject and within-subject variance components as obtained by an analysis of variance, were tested. RESULTS: In addition to the expected influence of performance status, the main results were that the effects of training volume (P = 0.20) and training phases (P = 0.81) on tHb-mass were not significant. We found that within-subject variations mainly have an analytical origin (TE approximately 1.4%) and a very small SD (7.5 g) of biological origin. CONCLUSION: tHb-mass shows very low individual oscillations during a training year (<6%), and these oscillations are below the expected changes in tHb-mass due to Herythropoetin (EPO) application or blood infusion (approximately 10%). The high stability of tHb-mass over a period of 1 year suggests that it should be included in an athlete's biological passport and analyzed by recently developed probabilistic inference techniques that define subject-based reference ranges.

Maternal and paternal contributions to pathogen resistance dependent on development stage in a whitefish (Salmonidae)

It is often assumed that maternal and paternal contributions to offspring phenotype change over the lifetime of an individual. However, studies on parental effects typically suffer from the problems that heritabilities and maternal environmental effects are difficult to separate, and that both may depend on environmental factors and developmental stage. In order to experimentally disentangle maternal from paternal contributions and the likely effects of developmental stage from ecological effects, we sampled a natural population of the whitefish Coregonus palaea, used gametes for full-factorial in vitro fertilizations, raised over 10000 of the resulting offspring singly at controlled conditions, and exposed them at different points during embryonic development to one of two strains of Pseudomonas fluorescens that differed in their virulence characteristics (only one caused mortality, while both delayed hatching and reduced growth). Vulnerability to infection increased markedly over embryo development. This change coincided with a distinct shift in the importance of maternal to additive genetic effects on survival. Timing of exposure also affected the variance components for hatching time and larval length, but in a less consistent direction than the variance components for mortality. No significant genetic variation was found for any reaction norms across time points of exposure, indicating a uniformity among genotypes in how susceptibility changed over development. Phenotypes were also typically correlated across time points, which could constrain the evolution of the reaction norms. Our experiment demonstrates that the relative maternal and paternal contributions to susceptibility to an infection, and hence the evolutionary potential to respond to pathogen-induced selection, depends not only on the kind of pathogenic stress but also on the timing of the challenge.

Genetic association tests: a method for the joint analysis of family and case-control data.

With the trend in molecular epidemiology towards both genome-wide association studies and complex modelling, the need for large sample sizes to detect small effects and to allow for the estimation of many parameters within a model continues to increase. Unfortunately, most methods of association analysis have been restricted to either a family-based or a case-control design, resulting in the lack of synthesis of data from multiple studies. Transmission disequilibrium-type methods for detecting linkage disequilibrium from family data were developed as an effective way of preventing the detection of association due to population stratification. Because these methods condition on parental genotype, however, they have precluded the joint analysis of family and case-control data, although methods for case-control data may not protect against population stratification and do not allow for familial correlations. We present here an extension of a family-based association analysis method for continuous traits that will simultaneously test for, and if necessary control for, population stratification. We further extend this method to analyse binary traits (and therefore family and case-control data together) and accurately to estimate genetic effects in the population, even when using an ascertained family sample. Finally, we present the power of this binary extension for both family-only and joint family and case-control data, and demonstrate the accuracy of the association parameter and variance components in an ascertained family sample.

Genome-wide association analysis identifies variants associated with nonalcoholic fatty liver disease that have distinct effects on metabolic traits.

Nonalcoholic fatty liver disease (NAFLD) clusters in families, but the only known common genetic variants influencing risk are near PNPLA3. We sought to identify additional genetic variants influencing NAFLD using genome-wide association (GWA) analysis of computed tomography (CT) measured hepatic steatosis, a non-invasive measure of NAFLD, in large population based samples. Using variance components methods, we show that CT hepatic steatosis is heritable (∼26%-27%) in family-based Amish, Family Heart, and Framingham Heart Studies (n = 880 to 3,070). By carrying out a fixed-effects meta-analysis of genome-wide association (GWA) results between CT hepatic steatosis and ∼2.4 million imputed or genotyped SNPs in 7,176 individuals from the Old Order Amish, Age, Gene/Environment Susceptibility-Reykjavik study (AGES), Family Heart, and Framingham Heart Studies, we identify variants associated at genome-wide significant levels (p<5×10(-8)) in or near PNPLA3, NCAN, and PPP1R3B. We genotype these and 42 other top CT hepatic steatosis-associated SNPs in 592 subjects with biopsy-proven NAFLD from the NASH Clinical Research Network (NASH CRN). In comparisons with 1,405 healthy controls from the Myocardial Genetics Consortium (MIGen), we observe significant associations with histologic NAFLD at variants in or near NCAN, GCKR, LYPLAL1, and PNPLA3, but not PPP1R3B. Variants at these five loci exhibit distinct patterns of association with serum lipids, as well as glycemic and anthropometric traits. We identify common genetic variants influencing CT-assessed steatosis and risk of NAFLD. Hepatic steatosis associated variants are not uniformly associated with NASH/fibrosis or result in abnormalities in serum lipids or glycemic and anthropometric traits, suggesting genetic heterogeneity in the pathways influencing these traits.

A genomewide linkage study on suicidality in major depressive disorder confirms evidence for linkage to 2p12.

Suicidal behavior is commonly associated with depression. Twin studies indicate that both suicidality and major depressive disorder (MDD) are heritable. However, epidemiological evidence suggests that the inheritance of suicidality is likely to be independent of the underlying psychiatric disorder, implying a distinct genetic contribution to suicidality. We conducted a genomewide linkage search aiming to detect genomic loci that may harbor susceptibility genes contributing to risk for suicidality in recurrent MDD. Affected sibling pair (ASP) variance components analysis was performed using the Depression Network cohort of 971 ASPs. The quantitative trait measuring suicidality as a broad phenotype, encompassing ideation and suicide attempts, was established from Schedules for Clinical Assessment in Neuropsychiatry interview items. We examined 1,060 genotyped microsatellite markers with an average spacing of 3.3 cM. Empirical thresholds for linkage evidence were set by whole-genome simulations (LOD = 2.71 for genomewide significance, 1.71 for suggestive linkage). No genomewide significant findings were found. Marker D3S1234 on 3p14 achieved suggestive linkage and yielded a maximum LOD of 1.853 (P = 0.0017), loci 9p24.3 and 18q22-q23 achieved LOD scores >1.5. We found some support for linkage to 2p12 (LOD = 1.2, P = 0.0087) which was previously implicated in linkage studies of suicidality. Our follow-up meta-analysis of five studies showed strong linkage to this region (P = 2 × 10(-6) ). In conclusion, this study analyzed suicidality as a continuous trait in MDD. We found modest evidence for linkage on 3p14. Our meta-analysis supports previous evidence of linkage to suicidality on 2p12. Some candidate genes in these regions may plausibly be implicated in suicidality.

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.

Multiple QTLs influencing triglyceride and HDL and total cholesterol levels identified in families with atherogenic dyslipidemia.

We conducted a genome-wide scan using variance components linkage analysis to localize quantitative-trait loci (QTLs) influencing triglyceride (TG), high density lipoprotein-cholesterol (HDL-C), low density lipoprotein-cholesterol, and total cholesterol (TC) levels in 3,071 subjects from 459 families with atherogenic dyslipidemia. The most significant evidence for linkage to TG levels was found in a subset of Turkish families at 11q22 [logarithm of the odds ratio (LOD)=3.34] and at 17q12 (LOD=3.44). We performed sequential oligogenic linkage analysis to examine whether multiple QTLs jointly influence TG levels in the Turkish families. These analyses revealed loci at 20q13 that showed strong epistatic effects with 11q22 (conditional LOD=3.15) and at 7q36 that showed strong epistatic effects with 17q12 (conditional LOD=3.21). We also found linkage on the 8p21 region for TG in the entire group of families (LOD=3.08). For HDL-C levels, evidence of linkage was identified on chromosome 15 in the Turkish families (LOD=3.05) and on chromosome 5 in the entire group of families (LOD=2.83). Linkage to QTLs for TC was found at 8p23 in the entire group of families (LOD=4.05) and at 5q13 in a subset of Turkish and Mediterranean families (LOD=3.72). These QTLs provide important clues for the further investigation of genes responsible for these complex lipid phenotypes. These data also indicate that a large proportion of the variance of TG levels in the Turkish population is explained by the interaction of multiple genetic loci.

The Learning Affect Monitor (LAM) - A computer-based system integrating quantitative and qualitative assessment of affective states in daily life.

The Learning Affect Monitor (LAM) is a new computer-based assessment system integrating basic dimensional evaluation and discrete description of affective states in daily life, based on an autonomous adapting system. Subjects evaluate their affective states according to a tridimensional space (valence and activation circumplex as well as global intensity) and then qualify it using up to 30 adjective descriptors chosen from a list. The system gradually adapts to the user, enabling the affect descriptors it presents to be increasingly relevant. An initial study with 51 subjects, using a 1 week time-sampling with 8 to 10 randomized signals per day, produced n = 2,813 records with good reliability measures (e.g., response rate of 88.8%, mean split-half reliability of .86), user acceptance, and usability. Multilevel analyses show circadian and hebdomadal patterns, and significant individual and situational variance components of the basic dimension evaluations. Validity analyses indicate sound assignment of qualitative affect descriptors in the bidimensional semantic space according to the circumplex model of basic affect dimensions. The LAM assessment module can be implemented on different platforms (palm, desk, mobile phone) and provides very rapid and meaningful data collection, preserving complex and interindividually comparable information in the domain of emotion and well-being.

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.

Bond strength tests of dental adhesive systems and their correlation with clinical results : a meta-analysis

OBJECTIVE: To evaluate the variability of bond strength test results of adhesive systems (AS) and to correlate the results with clinical parameters of clinical studies investigating cervical restorations. MATERIALS AND METHODS: Regarding the clinical studies, the internal database which had previously been used for a meta-analysis on cervical restorations was updated with clinical studies published between 2008 and 2012 by searching the PubMed and SCOPUS databases. PubMed and the International Association for Dental Research abstracts online were searched for laboratory studies on microtensile, macrotensile and macroshear bond strength tests. The inclusion criteria were (1) dentin, (2) testing of at least four adhesive systems, (3) same diameter of composite and (4) 24h of water storage prior to testing. The clinical outcome variables were retention loss, marginal discoloration, detectable margins, and a clinical index comprising the three parameters by weighing them. Linear mixed models which included a random study effect were calculated for both, the laboratory and the clinical studies. The variability was assessed by calculating a ratio of variances, dividing the variance among the estimated bonding effects obtained in the linear mixed models by the sum of all variance components estimated in these models. RESULTS: Thirty-two laboratory studies fulfilled the inclusion criteria comprising 183 experiments. Of those, 86 used the microtensile test evaluating 22 adhesive systems (AS). Twenty-seven used the macrotensile test with 17 AS, and 70 used the macroshear test with 24 AS. For 28 AS the results from clinical studies were available. Microtensile and macrotensile (Spearman rho=0.66, p=0.007) were moderately correlated and also microtensile and macroshear (Spearman rho=0.51, p=0.03) but not macroshear and macrotensile (Spearman rho=0.34, p=0.22). The effect of the adhesive system was significant for microtensile and macroshear (p<0.001) but not for macrotensile. The effect of the adhesive system could explain 36% of the variability of the microtensile test, 27% of the macrotensile and 33% of the macroshear test. For the clinical trials, about 49% of the variability of retained restorations could be explained by the adhesive system. With respect to the correlation between bond strength tests and clinical parameters, only a moderate correlation between micro- and macrotensile test results and marginal discoloration was demonstrated. However, no correlation between these tests and a retention loss or marginal integrity was shown. The correlation improved when more studies were included compared to assessing only one study. SIGNIFICANCE: The high variability of bond strength test results highlights the need to establish individual acceptance levels for a given test institute. The weak correlation of bond-strength test results with clinical parameters leads to the conclusion that one should not rely solely on bond strength tests to predict the clinical performance of an adhesive system but one should conduct other laboratory tests like tests on the marginal adaptation of fillings in extracted teeth and the retention loss of restorations in non-retentive cavities after artificial aging.

Genetic and environmental components of phenotypic variation in immune response and body size of a colonial bird, Delichon urbica (the house martin).

Directional selection for parasite resistance is often intense in highly social host species. Using a partial cross-fostering experiment we studied environmental and genetic variation in immune response and morphology in a highly colonial bird species, the house martin (Delichon urbica). We manipulated intensity of infestation of house martin nests by the haematophagous parasitic house martin bug Oeciacus hirundinis either by spraying nests with a weak pesticide or by inoculating them with 50 bugs. Parasitism significantly affected tarsus length, T cell response, immunoglobulin and leucocyte concentrations. We found evidence of strong environmental effects on nestling body mass, body condition, wing length and tarsus length, and evidence of significant additive genetic variance for wing length and haematocrit. We found significant environmental variance, but no significant additive genetic variance in immune response parameters such as T cell response to the antigenic phytohemagglutinin, immunoglobulins, and relative and absolute numbers of leucocytes. Environmental variances were generally greater than additive genetic variances, and the low heritabilities of phenotypic traits were mainly a consequence of large environmental variances and small additive genetic variances. Hence, highly social bird species such as the house martin, which are subject to intense selection by parasites, have a limited scope for immediate microevolutionary response to selection because of low heritabilities, but also a limited scope for long-term response to selection because evolvability as indicated by small additive genetic coefficients of variation is weak.

Genetic and phenotypic architecture of metabolic syndrome-associated components in dyslipidemic and normolipidemic subjects: the GEMS Study.

Atherogenic dyslipidemia, manifest by low HDL-cholesterol and high TG levels, is an important component of ATP-III defined metabolic syndrome. Here, we dissected the phenotypic and genetic architecture of these traits by assessing their relationships with other metabolically relevant measures, including plasma adipo-cytokines, highly sensitive C-reactive protein (hsCRP) and LDL particle size, in a large family data set (n=2800) and in an independent set of dyslipidemic cases (n=716) and normolipidemic controls (n=1073). We explored the relationships among these phenotypes using variable clustering and then estimated their genetic heritabilities and cross-trait correlations. In families, four clusters explained 61% of the total variance, with one adiposity-related cluster (including hsCRP), one BP-related cluster, and two lipid-related clusters (HDL-C, TG, adiponectin and LDL particle size; apoB and non-HDL-C). A similar structure was observed in dyslipidemic cases and normolipidemic controls. The genetic correlations in the families largely paralleled the phenotype clustering results, suggesting that common genes having pleiotropic effects contributed to the correlations observed. In summary, our analyses support a model of metabolic syndrome with two major components, body fat and lipids, each with two subcomponents, and quantifies their degree of overlap with each other and with metabolic-syndrome related measures (adipokines, LDL particle size and hsCRP).

Age-specific fitness components and their temporal variation in the barn owl.

Theory predicts that temporal variability plays an important role in the evolution of life histories, but empirical studies evaluating this prediction are rare. In constant environments, fitness can be measured by the population growth rate lambda, and the sensitivity of lambda to changes in fitness components estimates selection on these traits. In variable environments, fitness is measured by the stochastic growth rate lambda(S), and stochastic sensitivities estimate selection pressure. Here we examine age-specific schedules for reproduction and survival in a barn owl population (Tyto alba). We estimated how temporal variability affected fitness and selection, accounting for sampling variance. Despite large sample sizes of old individuals, we found no strong evidence for senescence. The most variable fitness components were associated with reproduction. Survival was less variable. Stochastic simulations showed that the observed variation decreased fitness by about 30%, but the sensitivities of lambda and lambda(S) to changes in all fitness components were almost equal, suggesting that temporal variation had negligible effects on selection. We obtained these results despite high observed variability in the fitness components and relatively short generation time of the study organism, a situation in which temporal variability should be particularly important for natural selection and early senescence is expected.

Comparison of the polynomial model against explicit measurements of noise components for different mammography systems.

Given the adverse impact of image noise on the perception of important clinical details in digital mammography, routine quality control measurements should include an evaluation of noise. The European Guidelines, for example, employ a second-order polynomial fit of pixel variance as a function of detector air kerma (DAK) to decompose noise into quantum, electronic and fixed pattern (FP) components and assess the DAK range where quantum noise dominates. This work examines the robustness of the polynomial method against an explicit noise decomposition method. The two methods were applied to variance and noise power spectrum (NPS) data from six digital mammography units. Twenty homogeneously exposed images were acquired with PMMA blocks for target DAKs ranging from 6.25 to 1600 µGy. Both methods were explored for the effects of data weighting and squared fit coefficients during the curve fitting, the influence of the additional filter material (2 mm Al versus 40 mm PMMA) and noise de-trending. Finally, spatial stationarity of noise was assessed.Data weighting improved noise model fitting over large DAK ranges, especially at low detector exposures. The polynomial and explicit decompositions generally agreed for quantum and electronic noise but FP noise fraction was consistently underestimated by the polynomial method. Noise decomposition as a function of position in the image showed limited noise stationarity, especially for FP noise; thus the position of the region of interest (ROI) used for noise decomposition may influence fractional noise composition. The ROI area and position used in the Guidelines offer an acceptable estimation of noise components. While there are limitations to the polynomial model, when used with care and with appropriate data weighting, the method offers a simple and robust means of examining the detector noise components as a function of detector exposure.