Multiple-trait random regression models for the estimation of genetic parameters for milk, fat, and protein yield in buffaloes

In this study, genetic parameters for test-day milk, fat, and protein yield were estimated for the first lactation. The data analyzed consisted of 1,433 first lactations of Murrah buffaloes, daughters of 113 sires from 12 herds in the state of São Paulo, Brazil, with calvings from 1985 to 2007. Ten-month classes of lactation days were considered for the test-day yields. The (co)variance components for the 3 traits were estimated using the regression analyses by Bayesian inference applying an animal model by Gibbs sampling. The contemporary groups were defined as herd-year-month of the test day. In the model, the random effects were additive genetic, permanent environment, and residual. The fixed effects were contemporary group and number of milkings (1 or 2), the linear and quadratic effects of the covariable age of the buffalo at calving, as well as the mean lactation curve of the population, which was modeled by orthogonal Legendre polynomials of fourth order. The random effects for the traits studied were modeled by Legendre polynomials of third and fourth order for additive genetic and permanent environment, respectively, the residual variances were modeled considering 4 residual classes. The heritability estimates for the traits were moderate (from 0.21-0.38), with higher estimates in the intermediate lactation phase. The genetic correlation estimates within and among the traits varied from 0.05 to 0.99. The results indicate that the selection for any trait test day will result in an indirect genetic gain for milk, fat, and protein yield in all periods of the lactation curve. The accuracy associated with estimated breeding values obtained using multi-trait random regression was slightly higher (around 8%) compared with single-trait random regression. This difference may be because to the greater amount of information available per animal. © 2013 American Dairy Science Association.

Genetic parameter estimation for pre- and post-weaning traits in Brahman cattle in Brazil

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

A simple method based on ICP-MS for estimation of background levels of arsenic, cadmium, copper, manganese, nickel, lead, and selenium in blood of the Brazillian population

Throughout the world, biomonitoring has become the standard for assessing exposure of individuals to toxic elements as well as for responding to serious environmental public health problems. However, extensive biomonitoring surveys require rapid and simple analytical methods. Thus, a simple and high-throughput method is proposed for the determination of arsenic (As), cadmium (Cd), copper (Cu), manganese (Mn), nickel (Ni), lead (Pb), and selenium (Se) in blood samples by using inductively coupled plasma–mass spectrometry (ICPMS). Prior to analysis, 200 ml of blood samples was mixed with 500 ml of 10% v/v tetramethylammonium hydroxide (TMAH) solution, incubated for 10 min, and subsequently diluted to 10 ml with a solution containing 0.05% w/v ethylenediamine tetraacetic acid (EDTA) + 0.005% v/v Triton X-100. After that, samples were directly analyzed by ICP-MS (ELAN DRC II). Rhodium was selected as an internal standard with matrix-matching calibration. Method detection limits were 0.08, 0.04, 0.5, 0.09, 0.12, 0.04, and 0.1 mg//L for As, Cd, Cu, Mn, Ni, Pb, and Se, respectively. Validation data are provided based on the analysis of blood samples from the trace elements inter-\comparison program operated by the Institut National de Santé Publique du Quebec, Canada. Additional validation was provided by the analysis of human blood samples by the proposed method and by using electrothermal atomic absorption spectrometry (ETAAS). The method was subsequently applied for the estimation of background metal blood values in the Brazilian population. In general, the mean concentrations of As, Cd, Cu, Mn, Ni, Pb, and Se in blood were 1.1, 0.4, 890, 9.6, 2.1, 65.4, and 89.3 mg/L, respectively, and are in agreement with other global populations. Influences of age, gender, smoking habits, alcohol consumption, and geographical variation on the values were also considered. Smoking habits influenced the levels of Cd in blood. The levels of Cu, Mn, and Pb were significantly correlated with gender, whereas Cu and Pb were significantly correlated with age. There were also interesting differences in Mn and Se levels in the population living in the north of Brazil compared to the south.

Modeling Age-Specific Mortality For Marine Mammal Populations

A method is presented for estimating age-specific mortality based on minimal information: a model life table and an estimate of longevity. This approach uses expected patterns of mammalian survivorship to define a general model of age-specific mortality rates. One such model life table is based on data for northern fur seals (Callorhinus ursinus) using Siler’s (1979) 5-parameter competing risk model. Alternative model life tables are based on historical data for human females and on a published model for Old World monkeys. Survival rates for a marine mammal species are then calculated by scaling these models by the longevity of that species. By using a realistic model (instead of assuming constant mortality), one can see more easily the real biological limits to population growth. The mortality estimation procedure is illustrated with examples of spotted dolphins (Stenella attenuata) and harbor porpoise (Phocoena phocoena).

Estimation of pelvic tilt on anteroposterior X-rays--a comparison of six parameters

OBJECTIVE: To compare six different parameters described in literature for estimation of pelvic tilt on an anteroposterior pelvic radiograph and to create a simple nomogram for tilt correction of prosthetic cup version in total hip arthroplasty. DESIGN: Simultaneous anteroposterior and lateral pelvic radiographs are taken routinely in our institution and were analyzed prospectively. The different parameters (including three distances and three ratios) were measured and compared to the actual pelvic tilt on the lateral radiograph using simple linear regression analysis. PATIENTS: One hundred and four consecutive patients (41 men, 63 women with a mean age of 31.7 years, SD 9.2 years, range 15.7-59.1 years) were studied. RESULTS: The strongest correlation between pelvic tilt and one of the six parameters for both men and women was the distance between the upper border of the symphysis and the sacrococcygeal joint. The correlation coefficient was 0.68 for men (P<0.001) and 0.61 for women (P<0.001). Based on this linear correlation, a nomogram was created that enables fast, tilt-corrected cup version measurements in clinical routine use. CONCLUSION: This simple method for correcting variations in pelvic tilt on plain radiographs can potentially improve the radiologist's ability to diagnose and interpret malformations of the acetabulum (particularly acetabular retroversion and excessive acetabular overcoverage) and post-operative orientation of the prosthetic acetabulum.

MULTIPLE DISEASES IN CARRIER PROBABILITY ESTIMATION: ACCOUNTING FOR SURVIVING ALL CANCERS OTHER THAN BREAST AND OVARY IN BRCAPRO

Mendelian models can predict who carries an inherited deleterious mutation of known disease genes based on family history. For example, the BRCAPRO model is commonly used to identify families who carry mutations of BRCA1 and BRCA2, based on familial breast and ovarian cancers. These models incorporate the age of diagnosis of diseases in relatives and current age or age of death. We develop a rigorous foundation for handling multiple diseases with censoring. We prove that any disease unrelated to mutations can be excluded from the model, unless it is sufficiently common and dependent on a mutation-related disease time. Furthermore, if a family member has a disease with higher probability density among mutation carriers, but the model does not account for it, then the carrier probability is deflated. However, even if a family only has diseases the model accounts for, if the model excludes a mutation-related disease, then the carrier probability will be inflated. In light of these results, we extend BRCAPRO to account for surviving all non-breast/ovary cancers as a single outcome. The extension also enables BRCAPRO to extract more useful information from male relatives. Using 1500 familes from the Cancer Genetics Network, accounting for surviving other cancers improves BRCAPRO’s concordance index from 0.758 to 0.762 (p = 0.046), improves its positive predictive value from 35% to 39% (p < 10−6) without impacting its negative predictive value, and improves its overall calibration, although calibration slightly worsens for those with carrier probability < 10%. Copyright c 2000 John Wiley & Sons, Ltd.

Gauss-Seidel Estimation of Generalized Linear Mixed Models with Application to Poisson Modeling of Spatially Varying Disease Rates

Generalized linear mixed models (GLMMs) provide an elegant framework for the analysis of correlated data. Due to the non-closed form of the likelihood, GLMMs are often fit by computational procedures like penalized quasi-likelihood (PQL). Special cases of these models are generalized linear models (GLMs), which are often fit using algorithms like iterative weighted least squares (IWLS). High computational costs and memory space constraints often make it difficult to apply these iterative procedures to data sets with very large number of cases. This paper proposes a computationally efficient strategy based on the Gauss-Seidel algorithm that iteratively fits sub-models of the GLMM to subsetted versions of the data. Additional gains in efficiency are achieved for Poisson models, commonly used in disease mapping problems, because of their special collapsibility property which allows data reduction through summaries. Convergence of the proposed iterative procedure is guaranteed for canonical link functions. The strategy is applied to investigate the relationship between ischemic heart disease, socioeconomic status and age/gender category in New South Wales, Australia, based on outcome data consisting of approximately 33 million records. A simulation study demonstrates the algorithm's reliability in analyzing a data set with 12 million records for a (non-collapsible) logistic regression model.

Checking Assumptions in Latent Class Regression Models via a Markov Chain Monte Carlo Estimation Approach: An Application to Depression and Socio-Economic Status

Latent class regression models are useful tools for assessing associations between covariates and latent variables. However, evaluation of key model assumptions cannot be performed using methods from standard regression models due to the unobserved nature of latent outcome variables. This paper presents graphical diagnostic tools to evaluate whether or not latent class regression models adhere to standard assumptions of the model: conditional independence and non-differential measurement. An integral part of these methods is the use of a Markov Chain Monte Carlo estimation procedure. Unlike standard maximum likelihood implementations for latent class regression model estimation, the MCMC approach allows us to calculate posterior distributions and point estimates of any functions of parameters. It is this convenience that allows us to provide the diagnostic methods that we introduce. As a motivating example we present an analysis focusing on the association between depression and socioeconomic status, using data from the Epidemiologic Catchment Area study. We consider a latent class regression analysis investigating the association between depression and socioeconomic status measures, where the latent variable depression is regressed on education and income indicators, in addition to age, gender, and marital status variables. While the fitted latent class regression model yields interesting results, the model parameters are found to be invalid due to the violation of model assumptions. The violation of these assumptions is clearly identified by the presented diagnostic plots. These methods can be applied to standard latent class and latent class regression models, and the general principle can be extended to evaluate model assumptions in other types of models.

Non-invasive quantitative estimation of bone density in rats throughout the life cycle and in arthritic osteopenia: preliminary results

A longitudinal bone survey was conducted in 86 female Wistar rats in order to assess mineral density kinetics from young age (5 weeks: 115 g) till late adulthood (64 weeks: 586 g). In vivo quantitative radiographic scanning was performed on the caudal vertebrae, taking trabecular mass as the parameter. Measurements were expressed as Relative Optical Density (ROD) units by means of a high resolution densitometric device. Results showed a progressive increase in mineral density throughout the life cycle, with a tendency to level in the higher weight range, indicating that progressive mineral aposition occurs in rats in dependency of age. This phenomenon, however, should be always considered within the context of continuous skeletal growth and related changes typical of this species. Twelve different animals were also examined following induction of articular inflammation with Freund's adjuvant in six of them. Bone survey conducted 12 to 18 days after inoculation revealed a significant (P less than 0.01) reduction in trabecular bone mass of scanned vertebrae in comparison with the weight-matched untreated controls. It is concluded that the in vivo quantitative assessment of bone density illustrated in this report represents a sensitive and useful tool for the long-term survey of naturally occurring or experimentally induced bone changes. Scanning of the same part of the skeleton can be repeated, thereby avoiding sacrifice of the animal and time-consuming preparation of post-mortem material.

Climate spectrum estimation in the presence of timescale errors

We introduce an algorithm (called REDFITmc2) for spectrum estimation in the presence of timescale errors. It is based on the Lomb-Scargle periodogram for unevenly spaced time series, in combination with the Welch's Overlapped Segment Averaging procedure, bootstrap bias correction and persistence estimation. The timescale errors are modelled parametrically and included in the simulations for determining (1) the upper levels of the spectrum of the red-noise AR(1) alternative and (2) the uncertainty of the frequency of a spectral peak. Application of REDFITmc2 to ice core and stalagmite records of palaeoclimate allowed a more realistic evaluation of spectral peaks than when ignoring this source of uncertainty. The results support qualitatively the intuition that stronger effects on the spectrum estimate (decreased detectability and increased frequency uncertainty) occur for higher frequencies. The surplus information brought by algorithm REDFITmc2 is that those effects are quantified. Regarding timescale construction, not only the fixpoints, dating errors and the functional form of the age-depth model play a role. Also the joint distribution of all time points (serial correlation, stratigraphic order) determines spectrum estimation.

On trend estimation under monotone Gaussian subordination with long-memory: application to fossil pollen series

Fossil pollen data from stratigraphic cores are irregularly spaced in time due to non-linear age-depth relations. Moreover, their marginal distributions may vary over time. We address these features in a nonparametric regression model with errors that are monotone transformations of a latent continuous-time Gaussian process Z(T). Although Z(T) is unobserved, due to monotonicity, under suitable regularity conditions, it can be recovered facilitating further computations such as estimation of the long-memory parameter and the Hermite coefficients. The estimation of Z(T) itself involves estimation of the marginal distribution function of the regression errors. These issues are considered in proposing a plug-in algorithm for optimal bandwidth selection and construction of confidence bands for the trend function. Some high-resolution time series of pollen records from Lago di Origlio in Switzerland, which go back ca. 20,000 years are used to illustrate the methods.

Individual risk factors for Mycoplasma hyopneumoniae infections in suckling pigs at the age of weaning.

BACKGROUND In recent years, the occurrence and the relevance of Mycoplasma hyopneumoniae infections in suckling pigs has been examined in several studies. Whereas most of these studies were focused on sole prevalence estimation within different age groups, follow-up of infected piglets or assessment of pathological findings, none of the studies included a detailed analysis of individual and environmental risk factors. Therefore, the aim of the present study was to investigate the frequency of M. hyopneumoniae infections in suckling pigs of endemically infected herds and to identify individual risk factors potentially influencing the infection status of suckling pigs at the age of weaning. RESULTS The animal level prevalence of M. hyopneumoniae infections in suckling pigs examined in three conventional pig breeding herds was 3.6% (41/1127) at the time of weaning. A prevalence of 1.2% was found in the same pigs at the end of their nursery period. In a multivariable Poisson regression model it was found that incidence rate ratios (IRR) for suckling pigs are significantly lower than 1 when teeth grinding was conducted (IRR: 0.10). Moreover, high temperatures in the piglet nest during the first two weeks of life (occasionally >40°C) were associated with a decrease of the probability of an infection (IRR: 0.23-0.40). Contrary, the application of PCV2 vaccines to piglets was associated with an increased infection risk (IRR: 9.72). CONCLUSIONS Since single infected piglets are supposed to act as initiators for the transmission of this pathogen in nursery and fattening pigs, the elimination of the risk factors described in this study should help to reduce the incidence rate of M. hyopneumoniae infections and thereby might contribute to a reduced probability of high prevalences in older pigs.

Factors influencing the accuracy of fetal weight estimation with a focus on preterm birth at the limit of viability: a systematic literature review

BACKGROUND Fetal weight estimation (FWE) is an important factor for clinical management decisions, especially in imminent preterm birth at the limit of viability between 23(0/7) and 26(0/7) weeks of gestation. It is crucial to detect and eliminate factors that have a negative impact on the accuracy of FWE. DATA SOURCES In this systematic literature review, we investigated 14 factors that may influence the accuracy of FWE, in particular in preterm neonates born at the limit of viability. RESULTS We found that gestational age, maternal body mass index, amniotic fluid index and ruptured membranes, presentation of the fetus, location of the placenta and the presence of multiple fetuses do not seem to have an impact on FWE accuracy. The influence of the examiner's grade of experience and that of fetal gender were discussed controversially. Fetal weight, time interval between estimation and delivery and the use of different formulas seem to have an evident effect on FWE accuracy. No results were obtained on the impact of active labor. DISCUSSION This review reveals that only few studies investigated factors possibly influencing the accuracy of FWE in preterm neonates at the limit of viability. Further research in this specific age group on potential confounding factors is needed.

P-T-t estimation of syn-kinematic strain in low-grade quartz-feldspar bearing rocks using thermodynamic modeling and 40Ar/39Ar dating techniques: example of the Plan-de-Phasy shear zone unit (Briançonnais Zone, Western Alps)

Pressure–Temperature–time (P–T–t) estimates of the syn-kinematic strain at the peak-pressure conditions reached during shallow underthrusting of the Briançonnais Zone in the Alpine subduction zone was made by thermodynamic modelling and 40Ar/39Ar dating in the Plan-de-Phasy unit (SE of the Pelvoux Massif, Western Alps). The dated phengite minerals crystallized syn-kinematically in a shear zone indicating top-to-the-N motion. By combining X-ray mapping with multi-equilibrium calculations, we estimate the phengite crystallization conditions at 270 ± 50 °C and 8.1 ± 2 kbar at an age of 45.9 ± 1.1 Ma. Combining this P–T–t estimate with data from the literature allows us to constrain the timing and geometry of Alpine continental subduction. We propose that the Briançonnais units were scalped on top of the slab during ongoing continental subduction and exhumed continuously until collision.

Estimation of reference curves for the urinary steroid metabolome in the first year of life in healthy children: tracing the complexity of human postnatal steroidogenesis

CONTEXT Complex steroid disorders such as P450 oxidoreductase deficiency or apparent cortisone reductase deficiency may be recognized by steroid profiling using chromatographic mass spectrometric methods. These methods are highly specific and sensitive, and provide a complete spectrum of steroid metabolites in a single measurement of one sample which makes them superior to immunoassays. The steroid metabolome during the fetal-neonatal transition is characterized by a) the metabolites of the fetal-placental unit at birth, b) the fetal adrenal androgens until its involution 3-6 months postnatally, and c) the steroid metabolites produced by the developing endocrine organs. All these developmental events change the steroid metabolome in an age- and sex-dependent manner during the first year of life. OBJECTIVE The aim of this study was to provide normative values for the urinary steroid metabolome of healthy newborns at short time intervals in the first year of life. METHODS We conducted a prospective, longitudinal study to measure 67 urinary steroid metabolites in 21 male and 22 female term healthy newborn infants at 13 time-points from week 1 to week 49 of life. Urine samples were collected from newborn infants before discharge from hospital and from healthy infants at home. Steroid metabolites were measured by gas chromatography-mass spectrometry (GC-MS) and steroid concentrations corrected for urinary creatinine excretion were calculated. RESULTS 61 steroids showed age and 15 steroids sex specificity. Highest urinary steroid concentrations were found in both sexes for progesterone derivatives, in particular 20α-DH-5α-DH-progesterone, and for highly polar 6α-hydroxylated glucocorticoids. The steroids peaked at week 3 and decreased by ∼80% at week 25 in both sexes. The decline of progestins, androgens and estrogens was more pronounced than of glucocorticoids whereas the excretion of corticosterone and its metabolites and of mineralocorticoids remained constant during the first year of life. CONCLUSION The urinary steroid profile changes dramatically during the first year of life and correlates with the physiologic developmental changes during the fetal-neonatal transition. Thus detailed normative data during this time period permit the use of steroid profiling as a powerful diagnostic tool.