Modeling and E-M estimation of haplotype-specific relative risks from genotype data for a case-control study of unrelated individuals

The US National Cancer Institute has recently sponsored the formation of a Cohort Consortium (http://2002.cancer.gov/scpgenes.htm) to facilitate the pooling of data on very large numbers of people, concerning the effects of genes and environment on cancer incidence. One likely goal of these efforts will be generate a large population-based case-control series for which a number of candidate genes will be investigated using SNP haplotype as well as genotype analysis. The goal of this paper is to outline the issues involved in choosing a method of estimating haplotype-specific risk estimates for such data that is technically appropriate and yet attractive to epidemiologists who are already comfortable with odds ratios and logistic regression. Our interest is to develop and evaluate extensions of methods, based on haplotype imputation, that have been recently described (Schaid et al., Am J Hum Genet, 2002, and Zaykin et al., Hum Hered, 2002) as providing score tests of the null hypothesis of no effect of SNP haplotypes upon risk, which may be used for more complex tasks, such as providing confidence intervals, and tests of equivalence of haplotype-specific risks in two or more separate populations. In order to do so we (1) develop a cohort approach towards odds ratio analysis by expanding the E-M algorithm to provide maximum likelihood estimates of haplotype-specific odds ratios as well as genotype frequencies; (2) show how to correct the cohort approach, to give essentially unbiased estimates for population-based or nested case-control studies by incorporating the probability of selection as a case or control into the likelihood, based on a simplified model of case and control selection, and (3) finally, in an example data set (CYP17 and breast cancer, from the Multiethnic Cohort Study) we compare likelihood-based confidence interval estimates from the two methods with each other, and with the use of the single-imputation approach of Zaykin et al. applied under both null and alternative hypotheses. We conclude that so long as haplotypes are well predicted by SNP genotypes (we use the Rh2 criteria of Stram et al. [1]) the differences between the three methods are very small and in particular that the single imputation method may be expected to work extremely well.     

ePAD,an oocyte and early embryo-abundant peptidylarginine deiminase-like protein that localizes to egg cytoplasmic sheets

Selected for its high relative abundance, a protein spot of MW approximately 75 kDa, pI 5.5 was cored from a Coomassie-stained two-dimensional gel of proteins from 2850 zona-free metaphase II mouse eggs and analyzed by tandem mass spectrometry (TMS), and novel microsequences were identified that indicated a previously uncharacterized egg protein. A 2.4-kb cDNA was then amplified from a mouse ovarian adapter-ligated cDNA library by RACE-PCR, and a unique 2043-bp open reading frame was defined encoding a 681-amino-acid protein. Comparison of the deduced amino acid sequence with the nonredundant database demonstrated that the protein was approximately 40% identical to the calcium-dependent peptidylarginine deiminase (PAD) enzyme family. Northern blotting, RT-PCR, and in situ hybridization analyses indicated that the protein was abundantly expressed in the ovary, weakly expressed in the testis, and absent from other tissues. Based on the homology with PADs and its oocyte-abundant expression pattern, the protein was designated ePAD, for egg and embryo-abundant peptidylarginine deiminase-like protein. Anti-recombinant ePAD monospecific antibodies localized the molecule to the cytoplasm of oocytes in primordial, primary, secondary, and Graafian follicles in ovarian sections, while no other ovarian cell type was stained. ePAD was also expressed in the immature oocyte, mature egg, and through the blastocyst stage of embryonic development, where expression levels began to decrease. Immunoelectron microscopy localized ePAD to egg cytoplasmic sheets, a unique keratin-containing intermediate filament structure found only in mammalian eggs and in early embryos, and known to undergo reorganization at critical stages of development. Previous reports that PAD-mediated deimination of epithelial cell keratin results in cytoskeletal remodeling suggest a possible role for ePAD in cytoskeletal reorganization in the egg and early embryo.     

CABYR, a novel calcium-binding tyrosine phosphorylation-regulated fibrous sheath protein involved in capacitation

To reach fertilization competence, sperm undergo an incompletely understood series of morphological and molecular maturational processes, termed capacitation, involving, among other processes, protein tyrosine phosphorylation and increased intracellular calcium. Hyperactivated motility and an ability to undergo the acrosome reaction serve as physiological end points to assess successful capacitation. We report here that acidic (pI 4.0) 86-kDa isoforms of a novel, polymorphic, testis-specific protein, designated calcium-binding tyrosine phosphorylation-regulated protein (CABYR), were tyrosine phosphorylated during in vitro capacitation and bound (45)Ca on 2D gels. Acidic 86-kDa calcium-binding forms of CABYR increased during in vitro capacitation, and calcium binding to these acidic forms was abolished by dephosphorylation with alkaline phosphatase. Six variants of CABYR containing two coding regions (CR-A and CR-B) were cloned from human testis cDNA libraries, including five variants with alternative splice deletions. A motif homologous to the RII dimerization domain of PK-A was present in the N-terminus of CR-A in four CABYR variants. A single putative EF handlike motif was noted in CR-A at aas 197-209, while seven potential tyrosine phosphorylation-like sites were noted in CR-A and four in CR-B. Pro-X-X-Pro (PXXP) modules were identified in the N- and C-termini of CR-A and CR-B. CABYR localizes to the principal piece of the human sperm flagellum in association with the fibrous sheath and is the first demonstration of a sperm protein that gains calcium-binding capacity when phosphorylated during capacitation.     

Isolation and identification of the major heparan sulfate proteoglycans in the developing bovine rib growth plate

Heparan sulfate proteoglycans are thought to mediate the action of growth factors. The heparan sulfate-containing proteoglycans in extracts of the bovine fetal rib growth plate were detected using the monoclonal antibody 3G10, which recognizes a neoepitope generated by heparitinase digestion (David, G., Bai, X. M., Van der Schueren, B., Cassiman, J. J., and Van den Berghe, H. (1992) J. Cell Biol. 119, 961-975). The heparan sulfate proteoglycans that react with this antibody were identified using antisera to known proteoglycans; purified using CsCl density gradient centrifugation, molecular sieve, and ion exchange chromatography; and then characterized. The major heparan sulfate proteoglycans in the growth plate had core proteins of 200 kDa and larger and were identified as perlecan and aggrecan. These two heparan sulfate proteoglycans could be effectively separated from each other by CsCl density gradient centrifugation alone. Perlecan contained 25% heparan sulfate and 75% chondroitin sulfate. The heparan sulfate chains on growth plate perlecan were considerably smaller than the chondroitin sulfate chains, and the heparan sulfate disaccharide content was different than that found for heparan sulfate from either kidney, tumor tissue, or growth plate aggrecan. Aggrecan contained only 0.1% heparan sulfate, which was localized to the CS-1 domain of aggrecan. These results indicate that perlecan and aggrecan would be the principal candidate proteoglycans involved in the action of heparan sulfate-binding proteins in the developing growth plate.     

Early embryo development in Fucus distichus is auxin sensitive

Auxin and polar auxin transport have been implicated in controlling embryo development in land plants. The goal of these studies was to determine if auxin and auxin transport are also important during the earliest stages of development in embryos of the brown alga Fucus distichus. Indole-3-acetic acid (IAA) was identified in F. distichus embryos and mature tissues by gas chromatography-mass spectroscopy. F. distichus embryos accumulate [(3)H]IAA and an inhibitor of IAA efflux, naphthylphthalamic acid (NPA), elevates IAA accumulation, suggesting the presence of an auxin efflux protein complex similar to that found in land plants. F. distichus embryos normally develop with a single unbranched rhizoid, but growth on IAA leads to formation of multiple rhizoids and growth on NPA leads to formation of embryos with branched rhizoids, at concentrations that are active in auxin accumulation assays. The effects of IAA and NPA are complete before 6 h after fertilization (AF), which is before rhizoid germination and cell division. The maximal effects of IAA and NPA are between 3.5 and 5 h AF and 4 and 5.5 h AF, respectively. Although, the location of the planes of cell division was significantly altered in NPA- and IAA-treated embryos, these abnormal divisions occurred after abnormal rhizoid initiation and branching was observed. The results of this study suggest that auxin acts in the formation of apical basal patterns in F. distichus embryo development.     

A lead-gill binding model to predict acute lead toxicity to rainbow trout (Oncorhynchus mykiss)

Rainbow trout (Oncorhynchus mykiss, approximately 2 g) were exposed to 0.6-1.0 microM Pb (125-200 microgl(-1)) for 3 h in ion-poor water. Complexing ligands (citrate, ethylenediamine, organic matter (OM)) or competing cations (Ca, Mg, Na) were added to the water. After exposure, trout gills were removed and analyzed for accumulated Pb. From these exposures, a conditional equilibrium binding constant (K) for Pb-gill binding was calculated (log K(Pb-gillPb)=6.0), plus conditional binding constants for cationic competition at the Pb binding sites and for Pb binding to OM in the water. These log K values were entered into the MINEQL+ aquatic chemistry equilibrium program, to calculate binding of Pb by trout gills. Two versions of the Pb-gill binding model were generated, one of which took into account OM quality as indicated by a simple measure of OM aromaticity, the specific absorption coefficient. The two model versions were tested against acute Pb toxicity (as the time to reach 50% fish mortality; LT50) during 1-week exposures of trout to 3.9 microM Pb in water collected from across southern Ontario. Both versions of the model generated highly significant correlations between the LT50 values and gill Pb concentrations calculated from measured exposure water chemistry, with the OM quality version correlating slightly better. Water pH also correlated well with the LT50 values, because the Pb exposures were in the pH range (7-8) where there is a nearly linear relationship between water pH and inorganic complexation of Pb. Advantages of the Pb-gill binding model include its completeness and the flexibility inherent in its conceptual framework, for example the inclusion of competition by Ca and H(+) for Pb binding sites on gills, and inclusion of complexation of Pb in the water column by natural OM and by carbonate.     

Proteomic analysis of rat soleus and tibialis anterior muscle following immobilization

A proteomic analysis was performed comparing normal slow twitch type fiber rat soleus muscle and normal fast twitch type fiber tibialis anterior muscle to immobilized soleus and tibialis anterior muscles at 0.5, 1, 2, 4, 6, 8 and 10 days post immobilization. Muscle mass measurements demonstrate mass changes throughout the period of immobilization. Proteomic analysis of normal and atrophied soleus muscle demonstrated statistically significant changes in the relative levels of 17 proteins. Proteomic analysis of normal and atrophied tibialis anterior muscle demonstrated statistically significant changes in the relative levels of 45 proteins. Protein identification using mass spectrometry was attempted for all differentially regulated proteins from both soleus and tibialis anterior muscles. Four differentially regulated soleus proteins and six differentially regulated tibialis anterior proteins were identified. The identified proteins can be grouped according to function as metabolic proteins, chaperone proteins, and contractile apparatus proteins. Together these data demonstrate that coordinated temporally regulated changes in the proteome occur during immobilization-induced atrophy in both slow twitch and fast twitch fiber type skeletal muscle.     

Proteomic analysis of rat soleus muscle undergoing hindlimb suspension-induced atrophy and reweighting hypertrophy

A proteomic analysis was performed comparing normal rat soleus muscle to soleus muscle that had undergone either 0.5, 1, 2, 4, 7, 10 and 14 days of hindlimb suspension-induced atrophy or hindlimb suspension-induced atrophied soleus muscle that had undergone 1 hour, 8 hour, 1 day, 2 day, 4 day and 7 days of reweighting-induced hypertrophy. Muscle mass measurements demonstrated continual loss of soleus mass occurred throughout the 21 days of hindlimb suspension; following reweighting, atrophied soleus muscle mass increased dramatically between 8 hours and 1 day post reweighting. Proteomic analysis of normal and atrophied soleus muscle demonstrated statistically significant changes in the relative levels of 29 soleus proteins. Reweighting following atrophy demonstrated statistically significant changes in the relative levels of 15 soleus proteins. Protein identification using mass spectrometry was attempted for all differentially regulated proteins from both atrophied and hypertrophied soleus muscle. Five differentially regulated proteins from the hindlimb suspended atrophied soleus muscle were identified while five proteins were identified in the reweighting-induced hypertrophied soleus muscles. The identified proteins could be generally grouped together as metabolic proteins, chaperone proteins and contractile apparatus proteins. Together these data demonstrate that coordinated temporally regulated changes in the skeletal muscle proteome occur during disuse-induced soleus muscle atrophy and reweighting hypertrophy.