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.     

Mouse divalent metal transporter 1 is a copper transporter in HEK293 cells

Divalent Metal Transporter 1 (DMT1) is an apical Fe transporter in the duodenum and is involved in endosomal Fe export. Four protein isoforms have been described for DMT1, two from mRNA with an iron responsive element (IRE) and two from mRNA without it. The sets of two begin in exon 1A or 2. We have characterized copper transport using mouse 2/?IRE DMT1 during regulated ectopic expression. HEK293 cells carrying a TetR:Hyg element were stably transfected with pDEST31 containing a 2/?IRE construct. ⁶⁴Cu¹⁺ incorporation in doxycycline treated cells exhibited 18.6 and 30.0-fold increases in Cu content, respectively when were exposed to 10 and 100 μM of extracellular Cu. Cu content was ~4-fold above that of parent cells or cells carrying just the vector. ⁶⁴Cu uptake in transfected cells pre-incubated with 5 μM of Cu-His revealed a Vmax and Km of 11.98 ± 0.52 pmol mg protein^?1 min^?1 and 2.03 ± 0.03 μM, respectively. Doxycycline-stimulated Cu uptake was linear with time. The rates of apical Cu uptake decreased and transepithelial transport increased when intracellular Cu increased. The optimal pH for Cu transport was 6.5; uptake of Cu was temperature dependent. Silver does not inhibit Cu uptake in cells carrying the vector. In conclusion, Cu uptake in HEK293 cells that over-expressed the 2/?IRE isoform of DMT1 transporter supports our earlier contention that DMT1 transports Cu as Cu¹⁺.     

The Discovery of orally available thrombin inhibitors : Studies towards the optimisation of CGH1668

The chemical optimisation of CGH1668 1 is described employing an in vivo model of absorption to determine the influence on bioavailability of single point modifications to five key molecular templates. The discovery of an orally bioavailable and selective thrombin inhibitor, 24, highlights the utility of this approach.     

Surrogate transformation of perennial ryegrass, Lolium perenne, using genetically modified Acremonium endophyte.

Summary Conditions have been developed for transforming protoplasts of the perennial ryegrass endophyteAcremonium strain 187BB. Unlike most other ryegrass endophytes, this strain does not produce the lolitrem B neurotoxin and is therefore suitable as a host for surrogate introduction of foreign genes into grasses. Transformation frequencies of 700–800 transformants/g DNA were obtained for both linear and circular forms of pAN7-1, a hygromycin (hph) resistant plasmid. Up to 80% of the linear transformants were stable on further culturing but only 25% of the circular transformants retained hygromycin resistance. Integration of pAN7-1 into the genome was confirmed by Southern blotting and probing of genomic digests of transformant DNA. Both single and tandemly repeated copies of the plasmid were found in the genome and both the number and sites of integration varied among the transformants. At least 13 chromosomes were identified in 187BB using contour-clamped homogeneous electric field (CHEF) gel electrophoresis. Probing of Southern blots of these gels confirmed that pAN7-1 had integrated into different chromosomes. The -glucuronidase (GUS) gene,uidA, was also introduced into 187BB by co-transformation of pNOM-2 with pAN7-1. GUS activity was detected by growing the transformants on plates containing 5-bromo-4-chloro-3-indolyl -D-glucuronic acid and by enzyme assays of mycelial extracts. Severalhph- anduidA-containing transformants were reintroduced into ryegrass seedlings and expression of GUS visualized in vivo, demonstrating that 187BB can be used as a surrogate host to introduce foreign genes into perennial ryegrass. Molecular analysis of fungal isolates from the leaf sheath confirmed that the pattern of pAN7-1 and pNOM-2 hybridizing fragments was identical to that observed in the fungus used as inoculum.     

Response and inbreeding from a genomic selection experiment in layer chickens

Background

Genomic selection (GS) using estimated breeding values (GS-EBV) based on dense marker data is a promising approach for genetic improvement. A simulation study was undertaken to illustrate the opportunities offered by GS for designing breeding programs. It consisted of a selection program for a sex-limited trait in layer chickens, which was developed by deterministic predictions under different scenarios. Later, one of the possible schemes was implemented in a real population of layer chicken.

Methods

In the simulation, the aim was to double the response to selection per year by reducing the generation interval by 50 %, while maintaining the same rate of inbreeding per year. We found that GS with retraining could achieve the set objectives while requiring 75 % fewer reared birds and 82 % fewer phenotyped birds per year. A multi-trait GS scenario was subsequently implemented in a real population of brown egg laying hens. The population was split into two sub-lines, one was submitted to conventional phenotypic selection, and one was selected based on genomic prediction. At the end of the 3-year experiment, the two sub-lines were compared for multiple performance traits that are relevant for commercial egg production.

Results

Birds that were selected based on genomic prediction outperformed those that were submitted to conventional selection for most of the 16 traits that were included in the index used for selection. However, although the two programs were designed to achieve the same rate of inbreeding per year, the realized inbreeding per year assessed from pedigree was higher in the genomic selected line than in the conventionally selected line.

Conclusions

The results demonstrate that GS is a promising alternative to conventional breeding for genetic improvement of layer chickens.     

Osmotic permeability of Novikoff hepatoma cells

The osmotic permeability coefficient (P_f) for water movement across Novikoff hepatoma cells was found to be 82 ± 3 (S.E.) · 10^?5 cm · s^?1 at 20°C. The corresponding diffusional permeability coefficient for ³HHO (P_d) was 97 ± 10 (S.E.) · 10^?5 cm · s^?1, therefore the ratio $\text{P}{}_{\mathrm{<mtext>f</mtext>}}\text{P}{}_{\mathrm{<mtext>d</mtext>}}$ is close to unity. The apparent activation energy for water filtration was 10.4 ± 0.4 (S.E.) kcal · mol^?1. This value is significantly greater than the activation energy for the self diffusion of water. The product of the hydraulic permeability coefficient and the viscosity coefficient for water was temperature-dependent. However, the product of the hydraulic permeability coefficient and the viscosity coefficient for membrane lipid did not vary with temperature. These data are interpreted as evidence for water movement across a lipid membrane barrier rather than through aqueous channels.     

Defining the cause of skewed X-chromosome inactivation in X-linked mental retardation by use of a mouse model

Extreme skewing of X-chromosome inactivation (XCI) is rare in the normal female population but is observed frequently in carriers of some X-linked mutations. Recently, it has been shown that various forms of X-linked mental retardation (XLMR) have a strong association with skewed XCI in female carriers, but the mechanisms underlying this skewing are unknown. ATR-X syndrome, caused by mutations in a ubiquitously expressed, chromatin-associated protein, provides a clear example of XLMR in which phenotypically normal female carriers virtually all have highly skewed XCI biased against the X chromosome that harbors the mutant allele. Here, we have used a mouse model to understand the processes causing skewed XCI. In female mice heterozygous for a null Atrx allele, we found that XCI is balanced early in embryogenesis but becomes skewed over the course of development, because of selection favoring cells expressing the wild-type Atrx allele. Unexpectedly, selection does not appear to be the result of general cellular-viability defects in Atrx-deficient cells, since it is restricted to specific stages of development and is not ongoing throughout the life of the animal. Instead, there is evidence that selection results from independent tissue-specific effects. This illustrates an important mechanism by which skewed XCI may occur in carriers of XLMR and provides insight into the normal role of ATRX in regulating cell fate.     

Joint estimation of migration rate and effective population size using the island model

Using the island model of population demography, I report that the demographic parameters migration rate and effective population size can be jointly estimated with equilibrium probabilities of identity in state calculated using a sample of genotypes collected at a single point in time from a single generation. The method, which uses moment-type estimators, applies to dioecious populations in which females and males have identical demography and monoecious populations with no selfing and requires that offspring genotypes are sampled following reproduction and prior to migration. I illustrate the estimation procedure using the infinite-island model with no mutation and the finite-island model with three kinds of mutation models. In the infinite-island model with no mutation, the estimators can be expressed as simple functions of estimates of the F-statistic parameters F(IT) and F(ST). In the finite-island model with mutation among k alleles, mutation rate, migration rate, and effective population size can be simultaneously estimated. The estimates of migration rate and effective population size are somewhat robust to violations in assumptions that may arise in empirical applications such as different kinds of mutation models and deviations from temporal equilibrium.     

Decoding Diffusivity in Multiple Sclerosis: Analysis of Optic Radiation Lesional and Non-Lesional White Matter

Objectives

Diffusion tensor imaging (DTI) has been suggested as a new promising tool in MS that may provide greater pathological specificity than conventional MRI, helping, therefore, to elucidate disease pathogenesis and monitor therapeutic efficacy. However, the pathological substrates that underpin alterations in brain tissue diffusivity are not yet fully delineated. Tract-specific DTI analysis has previously been proposed in an attempt to alleviate this problem. Here, we extended this approach by segmenting a single tract into areas bound by seemingly similar pathological processes, which may better delineate the potential association between DTI metrics and underlying tissue damage.

Method

Several compartments were segmented in optic radiation (OR) of 50 relapsing-remitting MS patients including T2 lesions, proximal and distal parts of fibers transected by lesion and fibers with no discernable pathology throughout the entire length of the OR.

Results

Asymmetry analysis between lesional and non-lesional fibers demonstrated a marked increase in Radial Diffusivity (RD), which was topographically limited to focal T2 lesions and potentially relates to the lesional myelin loss. A relative elevation of Axial Diffusivity (AD) in the distal part of the lesional fibers was observed in a distribution consistent with Wallerian degeneration, while diffusivity in the proximal portion of transected axons remained normal. A moderate, but significant elevation of RD in OR non-lesional fibers was strongly associated with the global (but not local) T2 lesion burden and is probably related to microscopic demyelination undetected by conventional MRI.

Conclusion

This study highlights the utility of the compartmentalization approach in elucidating the pathological substrates of diffusivity and demonstrates the presence of tissue-specific patterns of altered diffusivity in MS, providing further evidence that DTI is a sensitive marker of tissue damage in both lesions and NAWM. Our results suggest that, at least within the OR, parallel and perpendicular diffusivities are affected by tissue restructuring related to distinct pathological processes.