Fixed‐effect variance and the estimation of repeatabilities and heritabilities: issues and solutions

Linear mixed‐effects models are frequently used for estimating quantitative genetic parameters, including the heritability, as well as the repeatability, of traits. Heritability acts as a filter that determines how efficiently phenotypic selection translates into evolutionary change, whereas repeatability informs us about the individual consistency of phenotypic traits. As quantities of biological interest, it is important that the denominator, the phenotypic variance in both cases, reflects the amount of phenotypic variance in the relevant ecological setting. The current practice of quantifying heritabilities and repeatabilities from mixed‐effects models frequently deprives their denominator of variance explained by fixed effects (often leading to upward bias of heritabilities and repeatabilities), and it has been suggested to omit fixed effects when estimating heritabilities in particular. We advocate an alternative option of fitting models incorporating all relevant effects, while including the variance explained by fixed effects into the estimation of the phenotypic variance. The approach is easily implemented and allows optimizing the estimation of phenotypic variance, for example by the exclusion of variance arising from experimental design effects while still including all biologically relevant sources of variation. We address the estimation and interpretation of heritabilities in situations in which potential covariates are themselves heritable traits of the organism. Furthermore, we discuss complications that arise in generalized and nonlinear mixed models with fixed effects. In these cases, the variance parameters on the data scale depend on the location of the intercept and hence on the scaling of the fixed effects. Integration over the biologically relevant range of fixed effects offers a preferred solution in those situations.     

Contribution of physical fitness component to health status in middle-aged and elderly females

This study determined the physical fitness component that contributes to improving and maintaining health status for each age group as well as quantifying the degree of the relationship between health status and physical fitness in middle-aged and elderly females. The participants were 2,371 females aged 30 to 69 years. Ten physical fitness tests and medical checkups were performed. The participants were divided into a healthy group and an unhealthy group according to health status. Multiple discriminant analysis was applied to the multivariate data. Correct discriminant probabilities of the multiple discriminant function to discriminate the healthy and unhealthy groups for females ranged from 63.0% to 77.5%. These results suggest that there is a relatively high relationship between health status and physical fitness level for middle-aged and elderly females. With each individual's discriminant score calculated by the obtained multiple discriminant function as the index of the degree of health, the Pearson's correlation coefficient of the discriminant score and the performance in each physical fitness test were calculated. The aging change from 30 to 69 years old was classified into four patterns according to the contribution. The result of this study is considered to be useful as objective data to prepare an exercise program considering the contribution of the physical fitness component of health status.     

Effect of conformational states on protein dynamical transition

In order to examine the properties specific to the folded protein, the effect of the conformational states on protein dynamical transition was studied by incoherent elastic neutron scattering for both wild type and a deletion mutant of staphylococcal nuclease. The deletion mutant of SNase which lacks C-terminal 13 residues takes a compact denatured structure, and can be regarded as a model of intrinsic unstructured protein. Incoherent elastic neutron scattering experiments were carried out at various temperature between 10 K and 300 K on IN10 and IN13 installed at ILL. Temperature dependence of mean-square displacements was obtained by the q-dependence of elastic scattering intensity. The measurements were performed on dried and hydrated powder samples. No significant differences were observed between wild type and the mutant for the hydrated samples, while significant differences were observed for the dried samples. A dynamical transition at ∼ 140 K observed for both dried and hydrated samples. The slopes of the temperature dependence of MSD before transition and after transition are different between wild type and the mutant, indicating the folding induces hardening. The hydration water activates a further transition at ∼ 240 K. The behavior of the temperature dependence of MSD is indistinguishable for wild type and the mutant, indicating that hydration water dynamics dominate the dynamical properties.     

Na+/Mg2+ transporter acts as a Mg2+ buffering mechanism in PC12 cells

Mg(2+) buffering mechanisms in PC12 cells were demonstrated with particular focus on the role of the Na(+)/Mg(2+) transporter by using a newly developed Mg(2+) indicator, KMG-20, and also a Na(+) indicator, Sodium Green. Carbonyl cyanide p-(trifluoromethoxy) phenylhydrazone (FCCP), a protonophore, induced a transient increase in the intracellular Mg(2+) concentration ([Mg(2+)](i)). The rate of decrease of [Mg(2+)](i) was slower in a Na(+)-free extracellular medium, suggesting the coupling of Na(+) influx and Mg(2+) efflux. Na(+) influxes were different for normal and imipramine- (a putative inhibitor of the Na(+)/Mg(2+) transporter) containing solutions. FCCP induced a rapid increase in [Na(+)](i) in the normal solution, while the increase was gradual in the imipramine-containing solution. The rate of decrease of [Mg(2+)](i) in the imipramine-containing solution was also slower than that in the normal solution. From these results, we show that the main buffering mechanism for excess Mg(2+) depends on the Na(+)/Mg(2+) transporter in PC12 cells.     

Molecular taxonomy of dermatophytes and related fungi by chitin synthase 1 (CHS1) gene sequences

In the present study, the nucleotide sequences of the CHS1 gene from dermatophytes and related fungi in the genera Chrysosporium, Epidermophyton, Microsporum and Trichophyton were investigated using molecular methods. About 440-bp genomic DNA fragments of the CHS1 gene from 21 species were amplified by polymerase chain reaction (PCR) and sequenced. The CHS1 nucleotide sequences of these fungi showed more than 83% similarity. The molecular taxonomy of the CHS1 gene sequences revealed that Microsporum was genetically distinct from Chrysosporium and Trichophyton, as classified by morphological characteristics. This revised version was published online in June 2006 with corrections to the Cover Date.     

Development of nephrotic ICGN mice--the origin, reproductive ability, and incidence of glomerulonephritis

ICGN is a strain of spontaneous nephrotic mice with nonproliferative glomerular lesions. It was derived from an outbred Yok: ICR colony in our laboratory. The renal disease constantly occurred in animals of the first to the tenth generations (greater than 13.0%; 70 days of age). When affected males were mated with unaffected females, the incidence of the disease in their offspring was 38.8% (n = 49) at 70 days after birth. When both parents were affected, their offspring were all affected (n = 12). The disease evenly progressed in both sexes. It usually began 40 to 150 days after birth and death occurred within two months after onset. The animals usually showed sufficient reproductive ability as long as unaffected females were used for mating.     

Radioimmunoassay for non-enzymatically glycated serum proteins

We have developed a radioimmunoassay (RIA) for nonenzymatically glycated serum proteins. The polyclonal antibodies prepared against reduced glycated human albumin were specific for the glucitollysine residues of serum proteins. Serum proteins from diabetic patients (n = 25) contained 5.3 +/- 2.8 nmoles of glucitollysine/mg protein, compared to 2.0 +/- 0.2 in controls (n = 20). The intra- and inter-assay variables were 3.2-6.2% and 4.4-8.6%, respectively. Results from this assay procedure correlated well with those from the boronate affinity chromatography procedure (r = 0.94; P less than 0.001). The data suggested that diabetic serum proteins contained at least 2.5 times as much immunochemically detectable glucitollysine residures as normal serum proteins after reduction of the proteins with sodium borohydride.