Equivalence of multibreed animal models and hierarchical Bayes analysis for maternally influenced traits

Background

It has been argued that multibreed animal models should include a heterogeneous covariance structure. However, the estimation of the (co)variance components is not an easy task, because these parameters can not be factored out from the inverse of the additive genetic covariance matrix. An alternative model, based on the decomposition of the genetic covariance matrix by source of variability, provides a much simpler formulation. In this study, we formalize the equivalence between this alternative model and the one derived from the quantitative genetic theory. Further, we extend the model to include maternal effects and, in order to estimate the (co)variance components, we describe a hierarchical Bayes implementation. Finally, we implement the model to weaning weight data from an Angus × Hereford crossbred experiment.

Methods

Our argument is based on redefining the vectors of breeding values by breed origin such that they do not include individuals with null contributions. Next, we define matrices that retrieve the null-row and the null-column pattern and, by means of appropriate algebraic operations, we demonstrate the equivalence. The extension to include maternal effects and the estimation of the (co)variance components through the hierarchical Bayes analysis are then straightforward. A FORTRAN 90 Gibbs sampler was specifically programmed and executed to estimate the (co)variance components of the Angus × Hereford population.

Results

In general, genetic (co)variance components showed marginal posterior densities with a high degree of symmetry, except for the segregation components. Angus and Hereford breeds contributed with 50.26% and 41.73% of the total direct additive variance, and with 23.59% and 59.65% of the total maternal additive variance. In turn, the contribution of the segregation variance was not significant in either case, which suggests that the allelic frequencies in the two parental breeds were similar.

Conclusion

The multibreed maternal animal model introduced in this study simplifies the problem of estimating (co)variance components in the framework of a hierarchical Bayes analysis. Using this approach, we obtained for the first time estimates of the full set of genetic (co)variance components. It would be interesting to assess the performance of the procedure with field data, especially when interbreed information is limited.     

Managing and mining protein crystallization data

The crystallization of macromolecules remains a major bottleneck in structural biology. The routine screening of more than one thousand crystallization conditions and subsequent optimization by fine screening presents a challenge to conventional laboratory notebook keeping. In addition, the development of high-throughput robotic crystallization and imaging systems presents a pressing need for low-cost laboratory information management system (LIMS). Here we describe CLIMS2, a crystallization LIMS that features a simple, user-friendly graphical interface, allowing the storage, management, retrieval and mining of crystallization data. The CLIMS2 executable and documentation is freely available at http://clims.med.monash.edu.au.     

Phylogenetic analysis of the tenascin gene family: evidence of origin early in the chordate lineage

Background  

Tenascins are a family of glycoproteins found primarily in the extracellular matrix of embryos where they help to regulate cell proliferation, adhesion and migration. In order to learn more about their origins and relationships to each other, as well as to clarify the nomenclature used to describe them, the tenascin genes of the urochordate Ciona intestinalis, the pufferfish Tetraodon nigroviridis and Takifugu rubripes and the frog Xenopus tropicalis were identified and their gene organization and predicted protein products compared with the previously characterized tenascins of amniotes.     

Efficient large-scale protein sequence comparison and gene matching to identify orthologs and co-orthologs

Broadly, computational approaches for ortholog assignment is a three steps process: (i) identify all putative homologs between the genomes, (ii) identify gene anchors and (iii) link anchors to identify best gene matches given their order and context. In this article, we engineer two methods to improve two important aspects of this pipeline [specifically steps (ii) and (iii)]. First, computing sequence similarity data [step (i)] is a computationally intensive task for large sequence sets, creating a bottleneck in the ortholog assignment pipeline. We have designed a fast and highly scalable sort-join method (afree) based on k-mer counts to rapidly compare all pairs of sequences in a large protein sequence set to identify putative homologs. Second, availability of complex genomes containing large gene families with prevalence of complex evolutionary events, such as duplications, has made the task of assigning orthologs and co-orthologs difficult. Here, we have developed an iterative graph matching strategy where at each iteration the best gene assignments are identified resulting in a set of orthologs and co-orthologs. We find that the afree algorithm is faster than existing methods and maintains high accuracy in identifying similar genes. The iterative graph matching strategy also showed high accuracy in identifying complex gene relationships. Standalone afree available from http://vbc.med.monash.edu.au/～kmahmood/afree. EGM2, complete ortholog assignment pipeline (including afree and the iterative graph matching method) available from http://vbc.med.monash.edu.au/～kmahmood/EGM2.     

X-ray Crystal Structure and Specificity of the Plasmodium falciparum Malaria Aminopeptidase PfM18AAP

The malarial aminopeptidases have emerged as promising new drug targets for the development of novel antimalarial drugs. The M18AAP of Plasmodium falciparum malaria is a metallo-aminopeptidase that we show demonstrates a highly restricted specificity for peptides with an N-terminal Glu or Asp residue. Thus, the enzyme may function alongside other aminopeptidases in effecting the complete degradation or turnover of proteins, such as host hemoglobin, which provides a free amino acid pool for the growing parasite. Inhibition of PfM18AAP's function using antisense RNA is detrimental to the intra-erythrocytic malaria parasite and, hence, it has been proposed as a potential novel drug target. We report the X-ray crystal structure of the PfM18AAP aminopeptidase and reveal its complex dodecameric assembly arranged via dimer and trimer units that interact to form a large tetrahedron shape that completely encloses the 12 active sites within a central cavity. The four entry points to the catalytic lumen are each guarded by 12 large flexible loops that could control substrate entry into the catalytic sites. PfM18AAP thus resembles a proteasomal-like machine with multiple active sites able to degrade peptide substrates that enter the central lumen. The Plasmodium enzyme shows significant structural differences around the active site when compared to recently determined structures of its mammalian and human homologs, which provides a platform from which a rational approach to inhibitor design of new malaria-specific drugs can begin.     

The role of strand 1 of the C β-sheet in the structure and function of α1-antitrypsin

Serpins inhibit cognate serine proteases involved in a number of important processes including blood coagulation and inflammation. Consequently, loss of serpin function or stability results in a number of disease states. Many of the naturally occurring mutations leading to disease are located within strand 1 of the C beta-sheet of the serpin. To ascertain the structural and functional importance of each residue in this strand, which constitutes the so-called distal hinge of the reactive center loop of the serpin, an alanine scanning study was carried out on recombinant alpha(1)-antitrypsin Pittsburgh mutant (P1 = Arg). Mutation of the P10' position had no effect on its inhibitory properties towards thrombin. Mutations to residues P7' and P9' caused these serpins to have an increased tendency to act as substrates rather than inhibitors, while mutations at P6' and P8' positions caused the serpin to behave almost entirely as a substrate. Mutations at the P6' and P8' residues of the C beta-sheet, which are buried in the hydrophobic core in the native structure, caused the serpin to become highly unstable and polymerize much more readily. Thus, P6' and P8' mutants of alpha(1)-antitrypsin had melting temperatures 14 degrees lower than wild-type alpha(1)-antitrypsin. These results indicate the importance of maintaining the anchoring of the distal hinge to both the inhibitory mechanism and stability of serpins, the inhibitory mechanism being particularly sensitive to any perturbations in this region. The results of this study allow more informed analysis of the effects of mutations found at these positions in disease-associated serpin variants.