Characterization of mitochondrial genotypes in the foundation herd of the Canchim beef cattle breed

The Canchim (5/8 Charolais + 3/8 Zebu) beef cattle breed was developed at Southeast-Embrapa Cattle to take advantage of hybrid vigor and to combine the higher growth rate and beef quality of Charolais with tropical adaptations of Zebu. The development of three lineages (old, new, and crossbred) has increased its genetic basis. The genotypic origin (Bos taurus or Bos indicus) of the mitochondrial DNA (mtDNA) of the Canchim breed was unknown. We characterized the mtDNA genotype of this founder herd by allele-specific polymerase chain reaction. The 173 founder Zebu females (62 Indubrasil, 3 Guzerat, and 108 Nellore) and their 6749 offspring were identified. The frequency of B. indicus mtDNA ranged from 1.15 to 2.05% among the descendants (N = 6404) of each maternal line with available DNA, and among animals that were alive (N = 689) in December 2007 among the three lineages. Though mtDNA characterization can be used to direct animal selection, the low frequency of B. indicus mtDNA impairs the evaluation of its effects on production traits in these animals. The high prevalence of B. taurus mtDNA in Canchim proves that the founder Zebu females from the Indubrasil, Guzerat and Nellore breeds were obtained from crosses of Zebu sires with local B. taurus dams.

Estimates of genetic trend for carcass traits in a commercial broiler line

Data from the slaughter of 24,001 chickens that were part of a selection program for the production of commercial broilers were used to estimate genetic trend for absolute carcass (CW), breast meat (BRW), and leg (LW) weights, and relative carcass (CY), breast meat (BRY), and leg (LY) weights. The components of (co) variance and breeding values of individuals were obtained by the restricted maximum likelihood method applied to animal models. The relationship matrix was composed of 132,442 birds. The models included as random effects, maternal additive genetic and permanent environmental for CW, BRW, LW, CY, and BRY, and only maternal permanent environmental for LY, besides the direct additive genetic and residual effects, and as fixed effects, hatch week, parents' mating group and sex. The estimates of genetic trend were obtained by average regression of breeding value on generation, and the average genetic trend was estimated by regression coefficients. The genetic trends for CW (+ 6.0336 g/generation), BRW (+ 3.6723 g/generation), LW (+ 1.5846 g/generation), CY (+ 0.1195%/generation), and BRY (+ 0.1388%/generation) were positive, and they were in accordance with the objectives of the selection program for these traits. The genetic trend for LY(-0.0019%/generation) was negative, possibly due to the strong emphasis on selection for BRY and the negative correlations between these two traits.

Genetic parameters for productive life traits and reproductive efficiency traits at 6 years in Nellore cattle

The objective of the present study was to estimate (co)variance components for length of productive life (LPL) and some alternative reproductive traits of 6-year-old Nellore cattle. The data set contained 57,410 records for age at first calving from Nellore females and was edited to remove animal records with uncertain paternity and cows with just one piece of calving information. Only animals with age at first calving ranging from 23 to 48 months and calving intervals between 11 and 24 months were kept for analysis. LPL and life production ( LP) were used to describe productive life. LPL was defined as the number of months a cow was kept in the herd until she was 6 years old, given that she was alive at first calving and LP was defined as total number of calves in that time. Four traits were used to describe reproductive traits: two breeding efficiencies on original scale were estimated using Wilcox and Tomar functions (BEW and BET, respectively), and two breeding efficiencies transformed (ASBEW and ASBET, respectively), using the function [arcsine (square root (BEi/100))]. Estimates of heritability for measures of LPL and LP were low and ranged from 0.04 to 0.05. Estimates of heritability for breeding efficiencies on original and transformed scales oscillated from 0.18 to 0.32. Estimates of genetic correlations ranged from -0.57 to 0.79 for LPL and other traits and from 0.28 to 0.63 for LP and other traits.

Genetic analysis of average annual productivity of Nellore breeding cows (COWPROD)

With the aim of estimating the coefficient of heritability of average annual productivity of Nellore cows (COWPROD), a data set from 24,855 animals with known pedigree was analyzed. COWPROD is defined as the amount (in kilograms) of weaned calves produced yearly by one cow during her remaining time in herd ignoring a fixed period of 365 days. COWPROD was calculated regarding three standards: a) based on the post-weaning weight from the calves ignoring any kind of adjustment (COWPROD_NAJ), b) adjusted weight for the fixed effects (COWPROD_AJFIX) and c) adjusted weight for the fixed effects and for the genetic merit of the sire (COWPROD_AJFIN). The obtained heritabilities were 0.15, 0.15 and 0.16 for COWPROD_NAJ, COWPROD_AJFIX and COWPROD_AJFIN, respectively. A complete set composed of 105,158 COWPROD records on 130,740 animals in pedigree was also analyzed for predicting the genetic merit of all animals in the data set and for the calculation of the genetic, phenotypic and residual trends. Ranking correlation was high for the adjusted and non-adjusted data, yet, for some of the animals, the difference among the genetic values was large. This would be an indication that it would be better to work always with the adjusted weaning weights. The genetic trend was positive, but was of small magnitude (0.26% of the trait average) and the residual trend was negative as a consequence of the large intensification of the production system, which has been occurring in the last years in the farms studied. The phenotypic trend was also negative and intermediate between the genetic and the residual ones.

Effects of polymorphisms of LHR and FSHR genes on sexual precocity in a Bos taurus x Bos indicus beef composite population

The purpose of the present research was to investigate the effects of polymorphisms of luteinizing hormone receptor (LHR) and follicle-stimulating hormone receptor (FSHR) genes, evaluated by polymerase chain reaction-restriction fragment length polymorphism in European-Zebu composite beef heifers from six different breed compositions. The polymorphism site analysis from digestion with HhaI and AluI restriction endonucleases allowed the genotype identification for LHR (TT, CT and CC) and FSHR (GG, CG and CC) genes. A high frequency of heterozygous animals was recorded in all breed compositions for both genes, except in two compositions for LHR. The probability of pregnancy (PP) at first breeding was used to evaluate the polymorphism effect on sexual precocity. The PP was analyzed as a binary trait, with a value of 1 (success) assigned to heifers that were diagnosed pregnant by rectal palpation and a value of 0 (failure) assigned to those that were not pregnant at that time. Heterozygous heifers showed a higher pregnancy rate (67 and 66% for LHR and FSHR genes, respectively), but no significant effects were observed for the genes studied (P=0.9188 and 0.8831 for LHR and FSHR, respectively) on the PP. These results do not justify the inclusion of LHR and FSHR restriction fragment length polymorphism markers in selection programs for sexual precocity in beef heifers. Nevertheless, these markers make possible the genotype characterization and may be used in additional studies to evaluate the genetic structure in other bovine populations.

Additive genetic relationship of longevity with fertility and production traits in Nellore cattle based on bivariate models

Survival or longevity is an economically important trait in beef cattle. The main inconvenience for its inclusion in selection criteria is delayed recording of phenotypic data and the high computational demand for including survival in proportional hazard models. Thus, identification of a longevity-correlated trait that could be recorded early in life would be very useful for selection purposes. We estimated the genetic relationship of survival with productive and reproductive traits in Nellore cattle, including weaning weight (WW), post-weaning growth (PWG), muscularity (MUSC), scrotal circumference at 18 months (SC18), and heifer pregnancy (HP). Survival was measured in discrete time intervals and modeled through a sequential threshold model. Five independent bivariate Bayesian analyses were performed, accounting for cow survival and the five productive and reproductive traits. Posterior mean estimates for heritability (standard deviation in parentheses) were 0.55 (0.01) for WW, 0.25 (0.01) for PWG, 0.23 (0.01) for MUSC, and 0.48 (0.01) for SC18. The posterior mean estimates (95% confidence interval in parentheses) for the genetic correlation with survival were 0.16 (0.13-0.19), 0.30 (0.25-0.34), 0.31 (0.25-0.36), 0.07 (0.02-0.12), and 0.82 (0.78-0.86) for WW, PWG, MUSC, SC18, and HP, respectively. Based on the high genetic correlation and heritability (0.54) posterior mean estimates for HP, the expected progeny difference for HP can be used to select bulls for longevity, as well as for post-weaning gain and muscle score.

SMALL MAMMAL PREY SELECTION BY TWO OWL SPECIES IN SOUTHEASTERN BRAZIL

A variety of factors influence prey selection by predators. Because Barn Owls (Tyto alba) and Burrowing Owls (Athene cunicularia) differ in size and foraging tactics, we expected differential predation on small mammal prey. We hypothesized that the Barn Owl, all active predator, would prey on smaller and younger individuals than the Burrowing Owl, a sit-and-wait predator. We used pellet analyses to evaluate selection of small mammals by the two owls in relation to prey), species, age, and size at the Ecological Station of Itirapina, state of Sao Paulo, in southeastern Brazil. Small mammals constituted most of the prey individuals and biomass in the diet of Barn Owls. Although Burrowing Owls consumed a wider range of taxa, small mammals represented one-third of all biomass consumed. With respect. to small mammals, Barn Owls foraged selectively relative to prey species, size, and age. Burrowing Owls foraged opportunistically relative to prey species, but selectively relative to prey size and age. Barn Owls selected smaller and younger (juvenile and subadult) individuals of the delicate vesper mouse (Calomys tener) and Burrowing Owls preyed more oil larger and older (subadult only) individuals. morphology and behavior of both prey and predators may explain this differential predation. Our data suggest that the active predator feeds oil smaller and younger prey, and the sit-and-wait predator took relatively larger and older prey.

The number of reproductive workers in highly eusocial Hymenoptera: monogyny and monandry

Haplodiploidy results in relatedness asymmetries between colony members of highly eusocial Hymenoptera. As a consequence, queen and reproductive workers are more related to their own sons than to each other's male offspring. Kin selection theory predicts multiple optima in male parentage: either the queen or the workers should produce all the males. Nevertheless, shared male parentage is common in highly eusocial hymenopterans. An inclusive fitness model was used to analyze the effect of the number of reproductive workers on male parentage shared by the queen and laying workers by isolating the male component from an inclusive fitness equation using the equal fitness through male condition for each pairwise combination of the three female classes comprised of the queen, laying workers and non-laying workers. The main result of the theoretical analyses showed that the fraction of males produced by workers increases asymptotically with the number of laying workers at an increasingly diminishing rate, tending to an asymptotic value of 0.67. In addition, as the number of laying workers increases, the share of male parentage converges to that of non-laying workers. The diminishing return effect on male parentage share depending on the number of reproductive workers leads us to expect the number of reproductive workers to be relatively small in a stingless bee colony, even in the absence of productivity costs. The available data confirms this hypothesis, as there is an unusually small number of reproductive workers in stingless bee colonies.

The taxonomic status of the endangered thin-spined porcupine, Chaetomys subspinosus (Olfers, 1818), based on molecular and karyologic data

Background: The thin-spined porcupine, also known as the bristle-spined rat, Chaetomys subspinosus (Olfers, 1818), the only member of its genus, figures among Brazilian endangered species. In addition to being threatened, it is poorly known, and even its taxonomic status at the family level has long been controversial. The genus Chaetomys was originally regarded as a porcupine in the family Erethizontidae, but some authors classified it as a spiny-rat in the family Echimyidae. Although the dispute seems to be settled in favor of the erethizontid advocates, further discussion of its affinities should be based on a phylogenetic framework. In the present study, we used nucleotide-sequence data from the complete mitochondrial cytochrome b gene and karyotypic information to address this issue. Our molecular analyses included one individual of Chaetomys subspinosus from the state of Bahia in northeastern Brazil, and other hystricognaths. Results: All topologies recovered in our molecular phylogenetic analyses strongly supported Chaetomys subspinosus as a sister clade of the erethizontids. Cytogenetically, Chaetomys subspinosus showed 2n = 52 and FN = 76. Although the sexual pair could not be identified, we assumed that the X chromosome is biarmed. The karyotype included 13 large to medium metacentric and submetacentric chromosome pairs, one small subtelocentric pair, and 12 small acrocentric pairs. The subtelocentric pair 14 had a terminal secondary constriction in the short arm, corresponding to the nucleolar organizer region (Ag-NOR), similar to the erethizontid Sphiggurus villosus, 2n = 42 and FN = 76, and different from the echimyids, in which the secondary constriction is interstitial. Conclusion: Both molecular phylogenies and karyotypical evidence indicated that Chaetomys is closely related to the Erethizontidae rather than to the Echimyidae, although in a basal position relative to the rest of the Erethizontidae. The high levels of molecular and morphological divergence suggest that Chaetomys belongs to an early radiation of the Erethizontidae that may have occurred in the Early Miocene, and should be assigned to its own subfamily, the Chaetomyinae.

Genetic Variation among Major Human Geographic Groups Supports a Peculiar Evolutionary Trend in PAX9

A total of 172 persons from nine South Amerindian, three African and one Eskimo populations were studied in relation to the Paired box gene 9 (PAX9) exon 3 (138 base pairs) as well as its 5' and 3' flanking intronic segments (232 bp and 220 bp, respectively) and integrated with the information available for the same genetic region from individuals of different geographical origins. Nine mutations were scored in exon 3 and six in its flanking regions; four of them are new South American tribe-specific singletons. Exon3 nucleotide diversity is several orders of magnitude higher than its intronic regions. Additionally, a set of variants in the PAX9 and 101 other genes related with dentition can define at least some dental morphological differences between Sub-Saharan Africans and non-Africans, probably associated with adaptations after the modern human exodus from Africa. Exon 3 of PAX9 could be a good molecular example of how evolvability works.

Testing Evolutionary and Dispersion Scenarios for the Settlement of the New World

Background: Discussion surrounding the settlement of the New World has recently gained momentum with advances in molecular biology, archaeology and bioanthropology. Recent evidence from these diverse fields is found to support different colonization scenarios. The currently available genetic evidence suggests a ""single migration'' model, in which both early and later Native American groups derive from one expansion event into the continent. In contrast, the pronounced anatomical differences between early and late Native American populations have led others to propose more complex scenarios, involving separate colonization events of the New World and a distinct origin for these groups. Methodology/Principal Findings: Using large samples of Early American crania, we: 1) calculated the rate of morphological differentiation between Early and Late American samples under three different time divergence assumptions, and compared our findings to the predicted morphological differentiation under neutral conditions in each case; and 2) further tested three dispersal scenarios for the colonization of the New World by comparing the morphological distances among early and late Amerindians, East Asians, Australo-Melanesians and early modern humans from Asia to geographical distances associated with each dispersion model. Results indicate that the assumption of a last shared common ancestor outside the continent better explains the observed morphological differences between early and late American groups. This result is corroborated by our finding that a model comprising two Asian waves of migration coming through Bering into the Americas fits the cranial anatomical evidence best, especially when the effects of diversifying selection to climate are taken into account. Conclusions: We conclude that the morphological diversity documented through time in the New World is best accounted for by a model postulating two waves of human expansion into the continent originating in East Asia and entering through Beringia.

Diversification and Species Boundaries of Rhinebothrium (Cestoda; Rhinebothriidea) in South American Freshwater Stingrays (Batoidea; Potamotrygonidae)

Background: Neotropical freshwater stingrays (Batoidea: Potamotrygonidae) host a diverse parasite fauna, including cestodes. Both cestodes and their stingray hosts are marine-derived, but the taxonomy of this host/parasite system is poorly understood. Methodology: Morphological and molecular (Cytochrome oxidase I) data were used to investigate diversity in freshwater lineages of the cestode genus Rhinebothrium Linton, 1890. Results were based on a phylogenetic hypothesis for 74 COI sequences and morphological analysis of over 400 specimens. Cestodes studied were obtained from 888 individual potamotrygonids, representing 14 recognized and 18 potentially undescribed species from most river systems of South America. Results: Morphological species boundaries were based mainly on microthrix characters observed with scanning electron microscopy, and were supported by COI data. Four species were recognized, including two redescribed (Rhinebothrium copianullum and R. paratrygoni), and two newly described (R. brooksi n. sp. and R. fulbrighti n. sp.). Rhinebothrium paranaensis Menoret & Ivanov, 2009 is considered a junior synonym of R. paratrygoni because the morphological features of the two species overlap substantially. The diagnosis of Rhinebothrium Linton, 1890 is emended to accommodate the presence of marginal longitudinal septa observed in R. copianullum and R. brooksi n. sp. Patterns of host specificity and distribution ranged from use of few host species in few river basins, to use of as many as eight host species in multiple river basins. Significance: The level of intra-specific morphological variation observed in features such as total length and number of proglottids is unparalleled among other elasmobranch cestodes. This is attributed to the large representation of host and biogeographical samples. It is unclear whether the intra-specific morphological variation observed is unique to this freshwater system. Nonetheless, caution is urged when using morphological discontinuities to delimit elasmobranch cestode species because the amount of variation encountered is highly dependent on sample size and/or biogeographical representation.

Detection of polymorphisms of the mtDNA control region of Caretta caretta (Testudines: Cheloniidae) by PCR-SSCP

Marine turtles are increasingly being threatened worldwide by anthropogenic activities. Better understanding of their life cycle, behavior and population structure is imperative for the design of adequate conservation strategies. The mtDNA control region is a fast-evolving matrilineal marker that has been employed in the study of marine turtle populations. We developed and tested a simple molecular tracing system for Caretta caretta mtDNA haplotypes by polymerase chain reaction-single strand conformation polymorphism (PCR-SSCP). Using this technique, we were able to distinguish the SSCP patterns of 18 individuals of the haplotypes CC-A4, CC-A24 and CCxLO, which are commonly found in turtles sampled on the Brazilian coast. When we analyzed 15 turtles with previously unknown sequences, we detected two other haplotypes, in addition to the other four. Based on DNA sequencing, they were identified as the CC-A17 and CC-A1 haplotypes. Further analyses were made with the sea turtles, Chelonia mydas (N = 8), Lepidochelys olivacea (N = 3) and Eretmochelys imbricata (N = 1), demonstrating that the PCR-SSCP technique is able to distinguish intra-and interspecific variation in the family Cheloniidae. We found that this technique can be useful for identifying sea turtle mtDNA haplotypes, reducing the need for sequencing.

Genome Size, Karyotype Polymorphism and Chromosomal Evolution in Trypanosoma cruzi

Background: The Trypanosoma cruzi genome was sequenced from a hybrid strain (CL Brener). However, high allelic variation and the repetitive nature of the genome have prevented the complete linear sequence of chromosomes being determined. Determining the full complement of chromosomes and establishing syntenic groups will be important in defining the structure of T. cruzi chromosomes. A large amount of information is now available for T. cruzi and Trypanosoma brucei, providing the opportunity to compare and describe the overall patterns of chromosomal evolution in these parasites. Methodology/Principal Findings: The genome sizes, repetitive DNA contents, and the numbers and sizes of chromosomes of nine strains of T. cruzi from four lineages (TcI, TcII, TcV and TcVI) were determined. The genome of the TcI group was statistically smaller than other lineages, with the exception of the TcI isolate Tc1161 (Jose-IMT). Satellite DNA content was correlated with genome size for all isolates, but this was not accompanied by simultaneous amplification of retrotransposons. Regardless of chromosomal polymorphism, large syntenic groups are conserved among T. cruzi lineages. Duplicated chromosome-sized regions were identified and could be retained as paralogous loci, increasing the dosage of several genes. By comparing T. cruzi and T. brucei chromosomes, homologous chromosomal regions in T. brucei were identified. Chromosomes Tb9 and Tb11 of T. brucei share regions of syntenic homology with three and six T. cruzi chromosomal bands, respectively. Conclusions: Despite genome size variation and karyotype polymorphism, T. cruzi lineages exhibit conservation of chromosome structure. Several syntenic groups are conserved among all isolates analyzed in this study. The syntenic regions are larger than expected if rearrangements occur randomly, suggesting that they are conserved owing to positive selection. Mapping of the syntenic regions on T. cruzi chromosomal bands provides evidence for the occurrence of fusion and split events involving T. brucei and T. cruzi chromosomes.

Plasma Superoxide Dismutase-1 as a Surrogate Marker of Vivax Malaria Severity

Background: Severe outcomes have been described for both Plasmodium falciparum and P. vivax infections. The identification of sensitive and reliable markers of disease severity is fundamental to improving patient care. An intense pro-inflammatory response with oxidative stress and production of reactive oxygen species is present in malaria. Inflammatory cytokines such as tumor necrosis factor-alpha (TNF-alpha) and antioxidant agents such as superoxide dismutase-1 (SOD-1) are likely candidate biomarkers for disease severity. Here we tested whether plasma levels of SOD-1 could serve as a biomarker of severe vivax malaria. Methodology/Principal Findings: Plasma samples were obtained from residents of the Brazilian Amazon with a high risk for P. vivax transmission. Malaria diagnosis was made by both microscopy and nested PCR. A total of 219 individuals were enrolled: non-infected volunteers (n = 90) and individuals with vivax malaria: asymptomatic (n = 60), mild (n = 50) and severe infection (n = 19). SOD-1 was directly associated with parasitaemia, plasma creatinine and alanine amino-transaminase levels, while TNF-alpha correlated only with the later enzyme. The predictive power of SOD-1 and TNF-alpha levels was compared. SOD-1 protein levels were more effective at predicting vivax malaria severity than TNF-alpha. For discrimination of mild infection, elevated SOD-1 levels showed greater sensitivity than TNF-alpha (76% vs. 30% respectively; p < 0.0001), with higher specificity (100% vs. 97%; p < 0.0001). In predicting severe vivax malaria, SOD-1 levels exhibited higher sensitivity than TNF-alpha (80% vs. 56%, respectively; p < 0.0001; likelihood ratio: 7.45 vs. 3.14; p, 0.0001). Neither SOD-1 nor TNF-alpha could discriminate P. vivax infections from those caused by P. falciparum. Conclusion: SOD-1 is a powerful predictor of disease severity in individuals with different clinical presentations of vivax malaria.