Random regression models to estimate test-day milk yield genetic parameters Holstein cows in Southeastern Brazil

A total of 152,145 weekly test-day milk yield records from 7317 first lactations of Holstein cows distributed in 93 herds in southeastern Brazil were analyzed. Test-day milk yields were classified into 44 weekly classes of DIM. The contemporary groups were defined as herd-year-week of test-day. The model included direct additive genetic, permanent environmental and residual effects as random and fixed effects of contemporary group and age of cow at calving as covariable, linear and quadratic effects. Mean trends were modeled by a cubic regression on orthogonal polynomials of DIM. Additive genetic and permanent environmental random effects were estimated by random regression on orthogonal Legendre polynomials. Residual variances were modeled using third to seventh-order variance functions or a step function with 1, 6,13,17 and 44 variance classes. Results from Akaike`s and Schwarz`s Bayesian information criterion suggested that a model considering a 7th-order Legendre polynomial for additive effect, a 12th-order polynomial for permanent environment effect and a step function with 6 classes for residual variances, fitted best. However, a parsimonious model, with a 6th-order Legendre polynomial for additive effects and a 7th-order polynomial for permanent environmental effects, yielded very similar genetic parameter estimates. (C) 2008 Elsevier B.V. All rights reserved.

Intra- and interspecific divergence in the nuclear sequences of the clock gene period in species of the Drosophila buzzatii cluster

We have measured nucleotide variation in the CLOCK/CYCLE heterodimer inhibition domain (CCID) of the clock X-linked gene period in seven species belonging to the Drosophila buzzatii cluster, namely D. buzzatii, Drosophila koepferae, Drosophila antonietae, Drosophila serido, Drosophila gouveai, Drosophila seriema and Drosophila borborema. We detected that the purifying selection is the main force driving the sequence evolution in period, in agreement with the important role of CCID in clock machinery. Our survey revealed that period provides valuable phylogenetic information that allowed to resolve phylogenetic relationships among D. gouveai, D. borborema and D. seriema, which composed a polytomic clade in preliminary studies. The analysis of patterns of intraspecific variation revealed two different lineages of period in D. koepferae, probably reflecting introgressive hybridization from D. buzzatii, in concordance with previous molecular data.

Evolutionary History of Arabidopsis thaliana Aminoacyl-tRNA Synthetase Dual-Targeted Proteins

Aminoacyl-transfer RNA (tRNA) synthetases (aaRS) are key players in translation and act early in protein synthesis by mediating the attachment of amino acids to their cognate tRNA molecules. In plants, protein synthesis may occur in three subcellular compartments (cytosol, mitochondria, and chloroplasts), which requires multiple versions of the protein to be correctly delivered to its proper destination. The organellar aaRS are nuclear encoded and equipped with targeting information at the N-terminal sequence, which enables them to be specifically translocated to their final location. Most of the aaRS families present organellar proteins that are dual targeted to mitochondria and chloroplasts. Here, we examine the dual targeting behavior of aaRS from an evolutionary perspective. Our results show that Arabidopsis thaliana aaRS sequences are a result of a horizontal gene transfer event from bacteria. However, there is no evident bias indicating one single ancestor (Cyanobacteria or Proteobacteria). The dual-targeted aaRS phylogenetic relationship was characterized into two different categories (paralogs and homologs) depending on the state recovered for both dual-targeted and cytosolic proteins. Taken together, our results suggest that the dual-targeted condition is a gain-of-function derived from gene duplication. Selection may have maintained the original function in at least one of the copies as the additional copies diverged.

Models for genetic evaluation of Nelore cattle mature body weight

Records of 18,770 Nelore animals, born from 1975 to 2002, in 8 herds participating in the Nelore Cattle Breeding Program, were analyzed to estimate genetic parameters for mature BW. The mature BW were analyzed as a single BW taken closest to 4.5 yr of age for each cow in the data file, considering BW starting from 2 (W2Y_S), 3 (W3Y_S), or 4 (W4Y_S) yr of age or as repeated records, including all BW starting from 2 (W2Y_R), 3 (W3Y_R), or 4 (W4Y_R) yr of age. The variance components were estimated by restricted maximum likelihood, fitting univariate and bivariate animal models, including weaning weight. The heritability estimates were 0.29, 0.34, 0.36, 0.41, 0.44, and 0.46 for W2Y_S, W3Y_S, W4Y_S, W2Y_R, W3Y_R, and W4Y_R, respectively. The repeatability estimates for W2Y_R, W3Y_R, and W4Y_R were 0.59, 0.64, and 0.72, respectively. Larger accuracy values associated with the EBV were obtained in the repeated records models. The results indicated the bivariate repeated records model as the most appropriate for analyzing mature BW.

Genetic associations between accumulated productivity, and reproductive and growth traits in Nelore cattle

The total meat yield in a beef cattle production cycle is economically very important and depends on the number of calves born per year or birth season, being directly related to reproductive potential. Accumulated Productivity (ACP) is an index that expresses a cow`s capacity to give birth regularly at a young age and to wean animals of greater body weight. Using data from cattle participating in the ""Program for Genetic Improvement of the Nelore Breed"" (PMGRN - Nelore Brasil), bi-trait analyses were performed using the Restricted Maximum Likelihood method based on an ACP animal model and the following traits: age at first calving (AFC), female body weight adjusted for 365 (BW365) and 450 (BW450) days of age, and male scrotal circumference adjusted for 365 (SC365), 450 (SC450), 550 (SC550) and 730 (SC730) days of age. Median estimated ACP heritability was 0.19 and the genetic correlations with AFC, BW365, BW450, SC365, SC450, SC550 and SC730 were 0.33, 0.70, 0.65, 0.08, 0.07, 0.12 and 0.16, respectively. ACP increased and AFC decreased over time, revealing that the selection criteria genetically improved these traits. Selection based on ACP appears to favor the heaviest females at 365 and 450 days of age who showed better reproductive performance as regards AFC. Scrotal circumference was not genetically associated with ACP. (C) 2007 Elsevier B.V. All rights reserved.

A general hazard model for lifetime data in the presence of cure rate

Historically, the cure rate model has been used for modeling time-to-event data within which a significant proportion of patients are assumed to be cured of illnesses, including breast cancer, non-Hodgkin lymphoma, leukemia, prostate cancer, melanoma, and head and neck cancer. Perhaps the most popular type of cure rate model is the mixture model introduced by Berkson and Gage [1]. In this model, it is assumed that a certain proportion of the patients are cured, in the sense that they do not present the event of interest during a long period of time and can found to be immune to the cause of failure under study. In this paper, we propose a general hazard model which accommodates comprehensive families of cure rate models as particular cases, including the model proposed by Berkson and Gage. The maximum-likelihood-estimation procedure is discussed. A simulation study analyzes the coverage probabilities of the asymptotic confidence intervals for the parameters. A real data set on children exposed to HIV by vertical transmission illustrates the methodology.

Genetic associations between age at first calving and heifer body weight and scrotal circumference in Nelore cattle

P>Age at first calving (AFC) measures the entry of heifers into the beef cattle production system. This trait can be used as a selection criterion for earlier reproductive performance. Using data from Nelore cattle participating in the `Program for Genetic Improvement of the Nelore Breed` (PMGRN-Nelore Brazil), bi-trait analyses were performed using the restricted maximum likelihood method, based on an AFC animal model and the following traits: female body weight adjusted to 365 (BW365) and 450 (BW450) days of age, and male scrotal circumference adjusted to 365 (SC365), 450 (SC450), 550 (SC550) and 730 (SC730) days of age. The heritability estimates for AFC ranged from 0.02 +/- 0.02 to 0.04 +/- 0.02. The estimates of additive direct heritabilities (with standard error) for BW365, BW450, SC365, SC450, SC550 and SC730 were 0.36 +/- 0.07, 0.38 +/- 0.07, 0.48 +/- 0.07, 0.65 +/- 0.07, 0.64 +/- 0.07 and 0.42 +/- 0.07, respectively, and the genetic correlations with AFC were -0.38, -0.33, 0.10, -0.13, -0.13 and 0.06, respectively. In the herds studied, selection for SC365, SC450, SC550 or SC730 should not cause genetic changes in AFC. Selection based on BW365 or BW450 would favor smaller AFC breeding values. However, the low magnitude of direct heritability estimates for AFC in these farms indicates that changes in phenotypical expression depend mostly on non-genetic factors.

Genetic and environmental factors affecting ultrasound measures of longissimus muscle area and backfat thickness in Nelore cattle

The aim of the present study was to evaluate the genetic and environmental factors affecting records of longissimus muscle area (LMA) and back fat thickness (BF) obtained between the 12th and 13th ribs, and rump fat thickness (RF) between the hook and pin bones, measured by real-time ultrasound in Nelore cattle. Also, weight records of 22,778 animals born from 1998 to 2003, in ten farms across six Brazilian states were used. Carcass traits as measured by ultrasound of the live animal were recorded from 2002 to 2004 in 2590 males and females with ages varying from 450 to 599 days. Fixed models including farm, year and season of birth, sex and type of feed effects, and the covariates age of dam (AOD) and age of animal at measurement were used to study the effect of environmental factors on these traits. The genetic parameters for LMA, BF and RF were estimated with two and three-trait animal models with 120-day weights using a restricted maximum likelihood method. All environmental effects significantly affected carcass traits, with the exception of year of birth for BF and RF and AOD for LMA. The heritability estimates for LMA, BF and RF were 0.35, 0.51 and 0.39, respectively. Standard errors obtained in one-trait analyses were from 0.07 to 0.09. Genetic correlation estimates between LMA and the two traits of subcutaneous fat were low (close to zero) and 0.74 between BF and RF, indicating that the selection for LMA should not cause antagonism in the genetic improvement of subcutaneous fat measured by real-time ultrasound. (C) 2007 Elsevier B.V. All fights reserved.

Curupira-1 and Curupira-2, two novel Mutator-like DNA transposons from the genomes of human parasites Schistosoma mansoni and Schistosoma japonicum

Transposons of the Mutator superfamily have been widely described in plants, but only recently have metazoan organisms been shown to harbour them. In this work we describe novel Mutator superfamily transposons from the genomes of the human parasites Schistosoma mansoni and S. japonicum, which we name Curupira-1 and Curupira-2. Curupira elements do not have Terminal Inverted Repeats (TIRs) at their extremities and generate Target Site Duplications (TSDs) of 9 base pairs. Curupira-2 transposons code for a conserved transposase and SWIM zinc finger domains, while Curupira-1 elements comprise these same domains plus a WRKY zinc finger. Alignment of transcript sequences from both elements back to the genomes indicates that they are subject to splicing to produce mature transcripts. Phylogenetic analyses indicate that these transposons represent a new lineage of metazoan Mutator-like elements with characteristics that are distinct from the recently described Phantom elements. Description of these novel schistosome transposons provides new insights in the evolution of transposable elements in schistosomes.

Evolutionary history of dog rabies in Brazil

Although dogs are considered to be the principal transmitter of rabies in Brazil, dog rabies had never been recorded in South America before European colonization. In order to investigate the evolutionary history of dog rabies virus (RABV) in Brazil, we performed a phylogenetic analysis of carnivore RABV isolates from around the world and estimated the divergence times for dog RABV in Brazil. Our estimate for the time of introduction of dog RABV into Brazil was the late-19th to early-20th century, which was later than the colonization period but corresponded to a period of increased immigration from Europe to Brazil. In addition, dog RABVs appeared to have spread to indigenous animals in Brazil during the latter half of the 20th century, when the development and urbanization of Brazil occurred. These results suggest that the movement of rabid dogs, along with human activities since the 19th century, promoted the introduction and expansion of dog RABV in Brazil.

Detection and molecular characterization of a canine piroplasm from Brazil

In the beginning of the 20th century, a new canine disease was reported in Brazil under the name ""nambiuvu"", whose etiological agent was called Rangelia vitalii, a distinct piroplasm that was shown to parasitize not only erythrocytes, but also leucocytes and endothelial cells. In this new century, more publications on R. vitalii were reported from Brazil, including an extensive study on its ultrastructural analysis, in addition to clinical, pathological, and epidemiological data on nambiuvu. However, a molecular analysis of R. vitalii has not been performed to date. In the present study, we performed molecular phylogenetic analyses of R. vitalii based on fragments of the genes 18S rRNA and the heat shock protein 70 (hsp70), amplified by PCR performed on blood samples derived from five clinical cases of dogs presumably infected with R. vitalii in southern Brazil. In addition, we examined Giemsa-stained thin blood smears from these same dogs. DNA sequences (604-bp) of the 18S rRNA gene obtained from the five dogs were identical to each other, and by Blast analysis, this sequence shared the highest degree of sequence identity (95%) with Babesia sp. China-BQ1. DNA sequences (1056-bp) of the hsp70 gene obtained from the five dogs were identical to each other, and by Blast analysis, this sequence shared the highest degree of sequence identity (87%) with Babesia bigemina. Phylogenetic analyses inferred from either of the two genes resulted in the newly genotype being placed in the Babesia spp.sensu stricto clade with very high bootstrap support (95-100%) in three analyses (Neighbor-Joining, Maximum parsimony, and Maximum likelihood). Giemsa-stained thin blood smears from the dogs were shown to contain piroplasm organisms within erythrocytes, monocytes and neutrophils (individual forms), and schizont-like forms within neutrophils, in accordance with literature reports of R. vitalii. Based on these results, we conclude that R. vitalii, the etiological agent of ""nambiuvu"" in southern Brazil, is a valid species of piroplasm. Further studies are required to evaluate the validity of the genus Rangelia. (C) 2011 Elsevier B.V. All rights reserved.

Forest-obligate Sabethes mosquitoes suggest palaeoecological perturbations

The origin of tropical forest diversity has been hotly debated for decades. Although specific mechanisms vary, many such explanations propose some vicariance in the distribution of species during glacial cycles and several have been supported by genetic evidence in Neotropical taxa. However, no consensus exists with regard to the extent or time frame of the vicariance events. Here, we analyse the cytochrome oxidase II mitochondrial gene of 250 Sabethes albiprivus B mosquitoes sampled from western Sao Paulo in Brazil. There was very low population structuring among collection sites (Phi(ST) = 0.03, P = 0.04). Historic demographic analyses and the contemporary geographic distribution of genetic diversity suggest that the populations sampled are not at demographic equilibrium. Three distinct mitochondrial clades were observed in the samples, one of which differed significantly in its geographic distribution relative to the other two within a small sampling area (similar to 70 x 35 km). This fact, supported by the inability of maximum likelihood analyses to achieve adequate fits to simple models for the population demography of the species, suggests a more complex history, possibly involving disjunct forest refugia. This hypothesis is supported by a genetic signal of recent population growth, which is expected if population sizes of this forest-obligate insect increased during the forest expansions that followed glacial periods. Although a time frame cannot be reliably inferred for the vicariance event leading to the three genetic clades, molecular clock estimates place this at similar to 1 Myr before present.

Non-stationary frequency analysis of extreme daily rainfall in Sao Paulo, Brazil

This work is an assessment of frequency of extreme values (EVs) of daily rainfall in the city of Sao Paulo. Brazil, over the period 1933-2005, based on the peaks-over-threshold (POT) and Generalized Pareto Distribution (GPD) approach. Usually. a GPD model is fitted to a sample of POT Values Selected With a constant threshold. However. in this work we use time-dependent thresholds, composed of relatively large p quantities (for example p of 0.97) of daily rainfall amounts computed from all available data. Samples of POT values were extracted with several Values of p. Four different GPD models (GPD-1, GPD-2, GPD-3. and GDP-4) were fitted to each one of these samples by the maximum likelihood (ML) method. The shape parameter was assumed constant for the four models, but time-varying covariates were incorporated into scale parameter of GPD-2. GPD-3, and GPD-4, describing annual cycle in GPD-2. linear trend in GPD-3, and both annual cycle and linear trend in GPD-4. The GPD-1 with constant scale and shape parameters is the simplest model. For identification of the best model among the four models WC used rescaled Akaike Information Criterion (AIC) with second-order bias correction. This criterion isolates GPD-3 as the best model, i.e. the one with positive linear trend in the scale parameter. The slope of this trend is significant compared to the null hypothesis of no trend, for about 98% confidence level. The non-parametric Mann-Kendall test also showed presence of positive trend in the annual frequency of excess over high thresholds. with p-value being virtually zero. Therefore. there is strong evidence that high quantiles of daily rainfall in the city of Sao Paulo have been increasing in magnitude and frequency over time. For example. 0.99 quantiles of daily rainfall amount have increased by about 40 mm between 1933 and 2005. Copyright (C) 2008 Royal Meteorological Society

Copia Retrotransposon in the Zaprionus Genus: Another Case of Transposable Element Sharing with the Drosophila melanogaster Subgroup

Copia is a retrotransposon that appears to be distributed widely among the Drosophilidae subfamily. Evolutionary analyses of regulatory regions have indicated that the Copia retrotransposon evolved through both positive and purifying selection, and that horizontal transfer (HT) could also explain its patchy distribution of the among the subfamilies of the melanogaster subgroup. Additionally, Copia elements could also have transferred between melanogaster subgroup and other species of Drosophilidae-D. willistoni and Z. tuberculatus. In this study, we surveyed seven species of the Zaprionus genus by sequencing the LTR-ULR and reverse transcriptase regions, and by using RT-PCR in order to understand the distribution and evolutionary history of Copia in the Zaprionus genus. The Copia element was detected, and was transcriptionally active, in all species investigated. Structural and selection analysis revealed Zaprionus elements to be closely related to the most ancient subfamily of the melanogaster subgroup, and they seem to be evolving mainly under relaxed purifying selection. Taken together, these results allowed us to classify the Zaprionus sequences as a new subfamily-ZapCopia, a member of the Copia retrotransposon family of the melanogaster subgroup. These findings indicate that the Copia retrotransposon is an ancient component of the genomes of the Zaprionus species and broaden our understanding of the diversity of retrotransposons in the Zaprionus genus.

Phylogenetic relationships and infrageneric classification of Epidendrum subgenus Amphiglottium (Laeliinae, Orchidaceae)

Epidendrum L. is the largest genus of Orchidaceae in the Neotropical region; it has an impressive morphological diversification, which imposes difficulties in delimitation of both infrageneric and interspecific boundaries. In this study, we review infrageneric boundaries within the subgenus Amphiglottium and try to contribute to the understanding of morphological diversification and taxa delimitation within this group. We tested the monophyly of the subgenus Amphiglottium sect. Amphiglottium, expanding previous phylogenetic investigations and reevaluated previous infrageneric classifications proposed. Sequence data from the trnL-trnF region were analyzed with both parsimony and maximum likelihood criteria. AFLP markers were also obtained and analyzed with phylogenetic and principal coordinate analyses. Additionally, we obtained chromosome numbers for representative species within the group. The results strengthen the monophyly of the subgenus Amphiglottium but do not support the current classification system proposed by previous authors. Only section Tuberculata comprises a well-supported monophyletic group, with sections Carinata and Integra not supported. Instead of morphology, biogeographical and ecological patterns are reflected in the phylogenetic signal in this group. This study also confirms the large variability of chromosome numbers for the subgenus Amphiglottium (numbers ranging from 2n = 24 to 2n = 240), suggesting that polyploidy and hybridization are probably important mechanisms of speciation within the group.