Comparing genomic prediction accuracy from purebred,crossbred and combined purebred and crossbred reference populations in sheep

Background

The accuracy of genomic prediction depends largely on the number of animals with phenotypes and genotypes. In some industries, such as sheep and beef cattle, data are often available from a mixture of breeds, multiple strains within a breed or from crossbred animals. The objective of this study was to compare the accuracy of genomic prediction for several economically important traits in sheep when using data from purebreds, crossbreds or a combination of those in a reference population.

Methods

The reference populations were purebred Merinos, crossbreds of Border Leicester (BL), Poll Dorset (PD) or White Suffolk (WS) with Merinos and combinations of purebred and crossbred animals. Genomic breeding values (GBV) were calculated based on genomic best linear unbiased prediction (GBLUP), using a genomic relationship matrix calculated based on 48 599 Ovine SNP (single nucleotide polymorphisms) genotypes. The accuracy of GBV was assessed in a group of purebred industry sires based on the correlation coefficient between GBV and accurate estimated breeding values based on progeny records.

Results

The accuracy of GBV for Merino sires increased with a larger purebred Merino reference population, but decreased when a large purebred Merino reference population was augmented with records from crossbred animals. The GBV accuracy for BL, PD and WS breeds based on crossbred data was the same or tended to decrease when more purebred Merinos were added to the crossbred reference population. The prediction accuracy for a particular breed was close to zero when the reference population did not contain any haplotypes of the target breed, except for some low accuracies that were obtained when predicting PD from WS and vice versa.

Conclusions

This study demonstrates that crossbred animals can be used for genomic prediction of purebred animals using 50 k SNP marker density and GBLUP, but crossbred data provided lower accuracy than purebred data. Including data from distant breeds in a reference population had a neutral to slightly negative effect on the accuracy of genomic prediction. Accounting for differences in marker allele frequencies between breeds had only a small effect on the accuracy of genomic prediction from crossbred or combined crossbred and purebred reference populations.     

Wearing a raincoat: exocrine secretions contain anti-wetting agents in the oribatid mite, <Emphasis Type="Italic">Liacarus subterraneus</Emphasis> (Acari: Oribatida)

Liacarus subterraneus is a large, soil-dwelling oribatid mite species that possesses a conspicuously shiny, clean and not wettable cuticular surface. The exocrine cuticular chemistry of this species was investigated by means of gas chromatography–mass spectrometry. Besides a fraction of hydrocarbons and a terpene, hexane extracts of whole mite bodies exhibited free carboxylic acids and their glycerides as main components. The compounds were arranged in three distinct extract profiles. Based on data from individual extracts, (1) the majority (more than 3/4) of specimens showed large amounts of 1,2-dioctanoyl-glycerol (and three other related esters) but no (or only traces of) free carboxylic acids. (2) In about 1/8 of extracts, free acids (mainly octanoic (caprylic) acid) and glycerides were detected. This second type of profile highly varied with respect to the relative abundance of acids and esters. (3) The third profile (in about 7% of specimens) exclusively exhibited free acids and no (or only traces of) glycerides. In addition, a few extracts exhibited no components at all. The extract compounds most likely originate from the lipid layer of the cerotegument of L. subterraneus. The cuticle of individuals that possessed extractable cerotegumental compounds (profile I, II, III) exhibited strong water repellent properties, while the cuticle of individuals that possessed no components in their extract did not. After hexane extraction, water repellent properties got lost. The distinct extract profiles detected most likely portray the stepwise generation of an anti-wetting, exocrine surface lipid layer of glycerides: If this layer is lost, fatty acids may be discharged again (profile III) and may subsequently esterify (profile II) to larger and more stable esters (diacyl-glycerols), eventually building up the “raincoat” (mainly profile I) of L. subterraneus.     

Circles with two tandem long terminal repeats are specifically cleaved by pol gene-associated endonuclease from avian sarcoma and leukosis viruses: nucleotide sequences required for site-specific cleavage.

The avian retroviral pol gene-encoded DNA endonuclease (pol-endo) has been shown to selectively cleave the viral long terminal repeat sequences (LTRs) in single-stranded DNA substrates in a region known to be joined to host DNA during integration (G. Duyk, J. Leis, M. Longiaru, and A.M. Skalka, Proc. Natl. Acad. Sci. USA 80:6745-6749, 1983). The preferred sites of cleavage were mapped to the unique U5/U3 junctions found only in covalently closed circular DNA molecules containing two tandem LTRs. The cuts occurred three nucleotides 5' to the axis of symmetry of the 12-of-15-base-pair nearly perfect inverted repeat which marks the LTR junction. Experiments with double-stranded supercoiled DNA substrates revealed a similar specificity for nicking. Also, the endonuclease associated with the pol cleavage product, pp32, has the same specificity as the alpha beta form. The limits of sequence required for site-selective cleavage near the U5/U3 junction were established with single-stranded DNA substrates. A domain no larger than 44 base pairs allowed site-selective cleavage in each strand in vitro. Recognition of either strand appeared to be independent of the other, and in each case, the critical sequence was asymmetrically distributed with respect to the U5/U3 junction. The predominant contribution was from the U5 domain; this is consistent with its conservation in the LTR sequences of a number of avian sarcoma and leukosis viruses.     

Mutations that alter the activity of the Rous sarcoma virus protease.

Mutations designed by analysis of the Rous sarcoma virus (RSV) and human immunodeficiency virus (HIV)-1 protease (PR) crystal structures were introduced into 1) the substrate binding pocket, 2) the substrate enclosing "flaps," and 3) surface loops of RSV PR. Each mutant PR was expressed in Escherichia coli. Changes in activity were detected by following cleavage of a truncated (NC-PR) precursor polypeptide in E. coli and cleavage of synthetic peptide substrates representing RSV and HIV-1 PR cleavage sites in vitro. Mutations in the substrate binding pocket exchanged amino acid residues located close to the substrate in the HIV-1 PR for structurally equivalent residues in the RSV PR. Changing histidine 65 to glycine (H65G) gave an inactive enzyme, while a double mutant R105P,G106V, as well as the triple mutant, H65G,R105P,G106V, produced enzymes which showed significant activity toward a substrate that represented a HIV-1 cleavage site. Mutating the catalytic aspartate (D37S) or an adjacent conserved alanine to threonine (A40T), produced inactive enzymes. In contrast, the substitution A40S was active, but showed a reduced rate of catalysis. Mutations in the flaps of conserved glycines (G69L, G70L) produced inactive PRs. Two extended RSV PR surface loops were shortened to the size found in HIV-1 PR and resulted in drastically reduced activity. These results have confirmed some of the basic predictions made from structural models but have also revealed unexpected roles and interactions in the protein.     

The small-subunit ribosomal RNA gene sequences from the hypotrichous ciliates Oxytricha nova and Stylonychia pustulata

We have determined the complete nucleotide sequence of the small- subunit ribosomal RNA genes for the ciliate protozoans Stylonychia pustulata and Oxytricha nova. The sequences are homologous and sufficiently similar that these organisms must be closely related. In a phylogeny inferred from comparisons of several eukaryotic small-subunit ribosomal RNAs, the divergence of the ciliates from the eukaryotic line of descent is seen to coincide with the radiation of the plants, the animals, and the fungi. This radiation is preceded by the divergence of the slime mold, Dictyostelium discoideum.      

25-Hydroxyvitamin D levels of children are inversely related to adiposity assessed by body mass index

Vitamin D deficiency is associated with wide range of pathologies. Some evidences have shown that low vitamin D circulating levels in children and adolescent are related to fat mass and obesity. The objectives of the present study were to characterize vitamin D status in children and adolescents and to determine if serum 25-hydroxyvitamin D (25(OH)D) concentration is related to adiposity assessed by body mass index (BMI). Serum 25(OH)D levels were measured by LIAISON method in 471 children and adolescents (2 to 18 years age) and analyzed according to gender, pubertal period, age, and BMI. An overall prevalence of 25(OH)D insufficiency and deficiency was present in the 67.1%. Lower 25(OH)D levels were found in females (25.56 ± 14.03 vs 29.71 ± 17.10 ng ml^?1; P = 0.004) and pubertal children (25.52 ± 13.97 vs 29.21 ± 16.83 ng ml^?1; P = 0.011). In addition, an inverse relation of BMI and age on 25(OH)D concentrations was observed in children. In conclusion, low vitamin D status was highly prevalent among children and adolescents. Of note, a non-lineal regression model showed that 39.6% of vitamin D levels variability was explained by BMI. These results indicate that adiposity assessed by BMI impacts vitamin D status.     

Altered Rous sarcoma virus Gag polyprotein processing and its effects on particle formation.

Proteolytic processing of the Rous sarcoma virus (RSV) Gag precursor was altered in vivo through the introduction of amino acid substitutions into either the polyprotein cleavage junctions or the PR coding sequence. Single amino acid substitutions (V(P2)S and P(P4)G), which are predicted from in vitro peptide substrate cleavage data to decrease the rate of release of PR from the Gag polyprotein, were placed in the NC portion of the NC-PR junction. These substitutions do not affect the efficiency of release of virus-like particles from COS cells even though recovered particles contain significant amounts of uncleaved Pr76gag in addition to mature viral proteins. Single amino acid substitutions (A(P3)F and S(P1)Y), which increase the rate of PR release from Gag, also do not affect budding of virus-like particles from cells. Substitution of the inefficiently cleaved MA-p2 junction sequence in Gag by eight amino acids from the rapidly cleaved NC-PR sequence resulted in a significant increase in cleavage at the new MA-p2 junction, but again without an effect on budding. However, decreased budding was observed when the A(P3)F or S(P1)Y substitution was included in the NC-PR junction sequence between the MA and p2 proteins. A budding defect was also caused by substitution into Gag of a PR subunit containing three amino acid substitutions (R105P, G106V, and S107N) in the substrate binding pocket that increase the catalytic activity of PR. The defect appears to be the result of premature proteolytic processing that could be rescued by inactivating PR through substitution of a serine for the catalytic aspartic acid residue. This budding defect was also rescued by single amino acid substitutions in the NC-PR cleavage site which decrease the rate of release of PR from Gag. A similar budding defect was caused by replacing the Gag PR with two PR subunits covalently linked by four glycine residues. In contrast to the defect caused by the triply substituted PR, the budding defect observed with the linked PR dimer could not be rescued by NC-PR cleavage site mutations, suggesting that PR dimerization is a limiting step in the maturation process. Overall, these results are consistent with a model in which viral protein maturation occurs after PR subunits are released from the Gag polyprotein.     

The avian retroviral integration protein cleaves the terminal sequences of linear viral DNA at the in vivo sites of integration.

The purified integration protein (IN) of avian myeloblastosis virus is shown to nick double-stranded oligodeoxynucleotide substrates that mimic the ends of the linear form of viral DNA. In the presence of Mg2+, nicks are created 2 nucleotides from the 3' OH ends of both the U5 plus strand and the U3 minus strand. Similar cleavage is observed in the presence of Mn2+ but only when the extent of the reaction is limited. Neither the complementary strands nor sequences representing the termini of human immunodeficiency virus type 1 DNA were cleaved at analogous positions. Analysis of a series of substrates containing U5 base substitutions has defined the sequence requirements for site-selective nicking; nucleotides near the cleavage site are most critical for activity. The minimum substrate size required to demonstrate significant activity corresponds to the nearly perfect 15-base terminal inverted repeat. This in vitro activity of IN thus produces viral DNA ends that are joined to host DNA in vivo and corresponds to an expected early step in the integrative recombination reaction. These results provide the first enzymatic support using purified retroviral proteins for a linear DNA precursor to the integrated provirus.