Computing approximate standard errors for genetic parameters derived from random regression models fitted by average information REML

Approximate standard errors (ASE) of variance components for random regression coefficients are calculated from the average information matrix obtained in a residual maximum likelihood procedure. Linear combinations of those coefficients define variance components for the additive genetic variance at given points of the trajectory. Therefore, ASE of these components and heritabilities derived from them can be calculated. In our example, the ASE were larger near the ends of the trajectory.     

The mechanism of stabilization of the structure of nuclease-T by binding of ligands.

The rate of unfolding of Nuclease-T at pH 8,20 degrees was determined as a function of concentration of the ligands deoxythymidine 3',5'-diphosphate (pdTp) and Ca2+ on the basis of the rate of exchange between free fragment, Nuclease-T(50-149) and labeled fragment, Nuclease-T-(50-149) incorporated in the structure of nuclease-T (Taniuchi, H. (1973) J. Biol. Chem. 248, 5164-5174). The rate constant of unfolding of unliganded Nuclease-T' was 4.6 times 10-4s-1. Those of Nuclease-T' bound with pdTp, with Ca2+, and with both pdtp and Ca2+ were 9.0 times 10-5, 1.6 times 10-4, and 2.2 times 10-5s-1, respectively. The association constants of pdTp and Ca2+ with Nuclease-T' were found to be 1.0 times 10-4 and 2.0 times 10-2 m-1, respectively. Those of pdTp with Nuclease-T' plus Ca2+ and of Ca2+ with Nuclease-T' plus pdTp were 4 times 10-5 and 1.4 times 10-4M-1, respectively. The calculation of free energy change on the basis of the association constants shows that the magnitude of negative free energy change involved in the binding of either of the two ligands increases by approximately 2 kcal when the other ligand is already bound. There is a correlation between the free energy change and the specifically coupled with the cooperative interacions operating throught the three-dimensional structure resulting in strengthening of the interactions throughtout the structure, including those with the ligands, without a large change in conformation.     

ArrayIDer: automated structural re-annotation pipeline for DNA microarrays

Background  

Systems biology modeling from microarray data requires the most contemporary structural and functional array annotation. However, microarray annotations, especially for non-commercial, non-traditional biomedical model organisms, are often dated. In addition, most microarray analysis tools do not readily accept EST clone names, which are abundantly represented on arrays. Manual re-annotation of microarrays is impracticable and so we developed a computational re-annotation tool (ArrayIDer) to retrieve the most recent accession mapping files from public databases based on EST clone names or accessions and rapidly generate database accessions for entire microarrays.     

Effect of hypermethylation of CCWGG sequences in DNA of Mesembryanthemum crystallinum plants on their adaptation to salt stress

Under salt stress conditions, the level of CpNpG-methylation (N is any nucleoside) of the nuclear genome of the facultative halophyte Mesembryanthemum crystallinum in the CCWGG sequences (W = A or T) increases two-fold and is coupled with hypermethylation of satellite DNA on switching-over of C3-photosynthesis to the crassulacean acid metabolism (CAM) pathway of carbon dioxide assimilation. The methylation pattern of the CCWGG sequences is not changed in both the 5'-promoter region of the gene of phosphoenolpyruvate carboxylase, the key enzyme of C4-photosynthesis and CAM, and in the nuclear ribosomal DNA. Thus, a specific CpNpG-hypermethylation of satellite DNA has been found under conditions of expression of a new metabolic program. The functional role of the CpNpG-hypermethylation of satellite DNA is probably associated with formation of a specialized chromatin structure simultaneously regulating expression of a large number of genes in the cells of M. crystallinum plants on their adaptation to salt stress and switching-over to CAM metabolism.     

Accelerating the reconstruction of genome-scale metabolic networks

Background  

The genomic information of a species allows for the genome-scale reconstruction of its metabolic capacity. Such a metabolic reconstruction gives support to metabolic engineering, but also to integrative bioinformatics and visualization. Sequence-based automatic reconstructions require extensive manual curation, which can be very time-consuming. Therefore, we present a method to accelerate the time-consuming process of network reconstruction for a query species. The method exploits the availability of well-curated metabolic networks and uses high-resolution predictions of gene equivalency between species, allowing the transfer of gene-reaction associations from curated networks.     

SIZ1 regulation of phosphate starvation-induced root architecture remodeling involves the control of auxin accumulation

Phosphate (Pi) limitation causes plants to modulate the architecture of their root systems to facilitate the acquisition of Pi. Previously, we reported that the Arabidopsis (Arabidopsis thaliana) SUMO E3 ligase SIZ1 regulates root architecture remodeling in response to Pi limitation; namely, the siz1 mutations cause the inhibition of primary root (PR) elongation and the promotion of lateral root (LR) formation. Here, we present evidence that SIZ1 is involved in the negative regulation of auxin patterning to modulate root system architecture in response to Pi starvation. The siz1 mutations caused greater PR growth inhibition and LR development of seedlings in response to Pi limitation. Similar root phenotypes occurred if Pi-deficient wild-type seedlings were supplemented with auxin. N-1-Naphthylphthalamic acid, an inhibitor of auxin efflux activity, reduced the Pi starvation-induced LR root formation of siz1 seedlings to a level equivalent to that seen in the wild type. Monitoring of the auxin-responsive reporter DR5::uidA indicated that auxin accumulates in PR tips at early stages of the Pi starvation response. Subsequently, DR5::uidA expression was observed in the LR primordia, which was associated with LR elongation. The time-sequential patterning of DR5::uidA expression occurred earlier in the roots of siz1 as compared with the wild type. In addition, microarray analysis revealed that several other auxin-responsive genes, including genes involved in cell wall loosening and biosynthesis, were up-regulated in siz1 relative to wild-type seedlings in response to Pi starvation. Together, these results suggest that SIZ1 negatively regulates Pi starvation-induced root architecture remodeling through the control of auxin patterning.     

A combined long-range phasing and long haplotype imputation method to impute phase for SNP genotypes

Background

Knowing the phase of marker genotype data can be useful in genome-wide association studies, because it makes it possible to use analysis frameworks that account for identity by descent or parent of origin of alleles and it can lead to a large increase in data quantities via genotype or sequence imputation. Long-range phasing and haplotype library imputation constitute a fast and accurate method to impute phase for SNP data.

Methods

A long-range phasing and haplotype library imputation algorithm was developed. It combines information from surrogate parents and long haplotypes to resolve phase in a manner that is not dependent on the family structure of a dataset or on the presence of pedigree information.

Results

The algorithm performed well in both simulated and real livestock and human datasets in terms of both phasing accuracy and computation efficiency. The percentage of alleles that could be phased in both simulated and real datasets of varying size generally exceeded 98% while the percentage of alleles incorrectly phased in simulated data was generally less than 0.5%. The accuracy of phasing was affected by dataset size, with lower accuracy for dataset sizes less than 1000, but was not affected by effective population size, family data structure, presence or absence of pedigree information, and SNP density. The method was computationally fast. In comparison to a commonly used statistical method (fastPHASE), the current method made about 8% less phasing mistakes and ran about 26 times faster for a small dataset. For larger datasets, the differences in computational time are expected to be even greater. A computer program implementing these methods has been made available.

Conclusions

The algorithm and software developed in this study make feasible the routine phasing of high-density SNP chips in large datasets.     

Nine new polymorphic microsatellite loci for the amplification of archived otolith DNA of Atlantic cod,Gadus morhua L.

Nine out of 22 microsatellite primers tested were successfully amplified on three samples of cod Gadus morhua L. (two contemporary and one archived otolith samples). All loci were polymorphic (5–23 alleles/locus). The average observed heterozygosity across loci and samples was 0.625, ranging from 0.294 to 0.895 at each locus. All loci were under Hardy–Weinberg equilibrium, except PGmo56 that showed significant excess of heterozygotes in all studied samples. The isolated loci were suitable for degraded DNA and therefore useful for conducting a long‐term temporal study with DNA obtained from archived otoliths of cod.