A Fast EM Algorithm for BayesA-Like Prediction of Genomic Breeding Values

Prediction accuracies of estimated breeding values for economically important traits are expected to benefit from genomic information. Single nucleotide polymorphism (SNP) panels used in genomic prediction are increasing in density, but the Markov Chain Monte Carlo (MCMC) estimation of SNP effects can be quite time consuming or slow to converge when a large number of SNPs are fitted simultaneously in a linear mixed model. Here we present an EM algorithm (termed “fastBayesA”) without MCMC. This fastBayesA approach treats the variances of SNP effects as missing data and uses a joint posterior mode of effects compared to the commonly used BayesA which bases predictions on posterior means of effects. In each EM iteration, SNP effects are predicted as a linear combination of best linear unbiased predictions of breeding values from a mixed linear animal model that incorporates a weighted marker-based realized relationship matrix. Method fastBayesA converges after a few iterations to a joint posterior mode of SNP effects under the BayesA model. When applied to simulated quantitative traits with a range of genetic architectures, fastBayesA is shown to predict GEBV as accurately as BayesA but with less computing effort per SNP than BayesA. Method fastBayesA can be used as a computationally efficient substitute for BayesA, especially when an increasing number of markers bring unreasonable computational burden or slow convergence to MCMC approaches.     

Accuracy of direct genomic breeding values for nationally evaluated traits in US Limousin and Simmental beef cattle

Background

In national evaluations, direct genomic breeding values can be considered as correlated traits to those for which phenotypes are available for traditional estimation of breeding values. For this purpose, estimates of the accuracy of direct genomic breeding values expressed as genetic correlations between traits and their respective direct genomic breeding values are required.

Methods

We derived direct genomic breeding values for 2239 registered Limousin and 2703 registered Simmental beef cattle genotyped with either the Illumina BovineSNP50 BeadChip or the Illumina BovineHD BeadChip. For the 264 Simmental animals that were genotyped with the BovineHD BeadChip, genotypes for markers present on the BovineSNP50 BeadChip were extracted. Deregressed estimated breeding values were used as observations in weighted analyses that estimated marker effects to derive direct genomic breeding values for each breed. For each breed, genotyped individuals were clustered into five groups using K-means clustering, with the aim of increasing within-group and decreasing between-group pedigree relationships. Cross-validation was performed five times for each breed, using four groups for training and the fifth group for validation. For each trait, we then applied a weighted bivariate analysis of the direct genomic breeding values of genotyped animals from all five validation sets and their corresponding deregressed estimated breeding values to estimate variance and covariance components.

Results

After minimizing relationships between training and validation groups, estimated genetic correlations between each trait and its direct genomic breeding values ranged from 0.39 to 0.76 in Limousin and from 0.29 to 0.65 in Simmental. The efficiency of selection based on direct genomic breeding values relative to selection based on parent average information ranged from 0.68 to 1.28 in genotyped Limousin and from 0.51 to 1.44 in genotyped Simmental animals. The efficiencies were higher for 323 non-genotyped young Simmental animals, born after January 2012, and ranged from 0.60 to 2.04.

Conclusions

Direct genomic breeding values show promise for routine use by Limousin and Simmental breeders to improve the accuracy of predicted genetic merit of their animals at a young age and increase response to selection. Benefits from selecting on direct genomic breeding values are greater for breeders who use natural mating sires in their herds than for those who use artificial insemination sires. Producers with unregistered commercial Limousin and Simmental cattle could also benefit from being able to identify genetically superior animals in their herds, an opportunity that has in the past been limited to seed stock animals.     

Surface Binding of Aflatoxin B1 by Lactic Acid Bacteria

Specific lactic acid bacterial strains remove toxins from liquid media by physical binding. The stability of the aflatoxin B₁ complexes formed with 12 bacterial strains in both viable and nonviable (heat- or acid-treated) forms was assessed by repetitive aqueous extraction. By the fifth extraction, up to 71% of the total aflatoxin B₁ remained bound. Nonviable bacteria retained the highest amount of aflatoxin B₁. Lactobacillus rhamnosus strain GG (ATCC 53103) and L. rhamnosus strain LC-705 (DSM 7061) removed aflatoxin B₁ from solution most efficiently and were selected for further study. The accessibility of bound aflatoxin B₁ to an antibody in an indirect competitive inhibition enzyme-linked immunosorbent assay suggests that surface components of these bacteria are involved in binding. Further evidence is the recovery of around 90% of the bound aflatoxin from the bacteria by solvent extraction. Autoclaving and sonication did not release any detectable aflatoxin B₁. Variation in temperature (4 to 37°C) and pH (2 to 10) did not have any significant effect on the amount of aflatoxin B₁ released. Binding of aflatoxin B₁ appears to be predominantly extracellular for viable and heat-treated bacteria. Acid treatment may permit intracellular binding. In all cases, binding is of a reversible nature, but the stability of the complexes formed depends on strain, treatment, and environmental conditions.     

Lipoprotein-dependent synthesis of prostacyclin in smooth muscle cells

Low (LDL) and high density lipoproteins (HDL) stimulated prostacycline (PGI2) synthesis in rabbit and human aorta smooth muscle cells growing in culture. The lipoproteins were added to the cells in concentrations equal to that of cholesterol. It was shown that HDL exerted a stronger stimulating effect as compared to LDL. The maximal effect was observed with HDL3. HDL3 isolated from blood serum of healthy volunteers appeared to be more active in PGI2 synthesis promotion than those of CDH patients with documented coronary atherosclerosis. Purified Apo A-1 stimulated the transformation of [14C]arachidonic acid into the products of its metabolism with increased accumulation of 6-keto-PGF1 alpha among labeled metabolites. Estradiol (1.10(-7) M) showed a stimulating effect; norepinephrine (1.10(-6) M) and progesterone (1.10(-7) M) showed an inhibiting effect, whereas corticosterone (1.10(-6) M) and deoxycorticosterone (1.10(-6) M) did not influence the rate of LDL-dependent PGI2 synthesis.     

Genome Analysis Traces Regional Dispersal of Rice in Taiwan and Southeast Asia

The dispersal of rice (Oryza sativa) following domestication influenced massive social and cultural changes across South, East, and Southeast (SE) Asia. The history of dispersal across islands of SE Asia, and the role of Taiwan and the Austronesian expansion in this process remain largely unresolved. Here, we reconstructed the routes of dispersal of O. sativa ssp. japonica rice to Taiwan and the northern Philippines using whole-genome resequencing of indigenous rice landraces coupled with archaeological and paleoclimate data. Our results indicate that japonica rice found in the northern Philippines diverged from Indonesian landraces as early as 3,500 years before present (BP). In contrast, rice cultivated by the indigenous peoples of the Taiwanese mountains has complex origins. It comprises two distinct populations, each best explained as a result of admixture between temperate japonica that presumably came from northeast Asia, and tropical japonica from the northern Philippines and mainland SE Asia, respectively. We find that the temperate japonica component of these indigenous Taiwan populations diverged from northeast Asia subpopulations at about 2,600 BP, whereas gene flow from the northern Philippines had begun before ∼1,300  BP. This coincides with a period of intensified trade established across the South China Sea. Finally, we find evidence for positive selection acting on distinct genomic regions in different rice subpopulations, indicating local adaptation associated with the spread of japonica rice.     

Blubber cortisol qualitatively reflects circulating cortisol concentrations in bottlenose dolphins

Stress hormones, released into circulation as a consequence of disturbance, are classically assayed from blood samples but may also be detected in a variety of matrices. Blubber and fecal samples can be remotely collected from free‐ranging cetaceans without the confounding hormone elevations associated with chase, capture, and handling required to collect blood samples. The relationship between cortisol concentrations in circulation with that of blubber and feces, however, is unknown. To assess these associations, we elevated cortisol by orally administering hydrocortisone for five days in five bottlenose dolphins. Voluntary blood and fecal samples were collected daily; blubber biopsies were collected on day one, just prior to hydrocortisone administration, and days three and five of hydrocortisone administration. We evaluated subsequent changes in several circulating stress hormones as well as cortisol and glucocorticoid metabolites in blubber and feces, respectively. There was a significant association between cortisol levels in serum and in blubber (F_1,12.7 = 14.3, P < 0.01, mR² = 0.57) despite substantial variability in blubber cortisol levels. Counterintuitively, fecal cortisol metabolite levels were inversely related to serum cortisol. The relationship between serum and blubber cortisol levels suggests blubber samples from remote sampling may be useful to detect stress loads in this species.