A heuristic method to identify runs of homozygosity associated with reduced performance in livestock

Although, for the most part, genome-wide metrics are currently used in managing livestock inbreeding, genomic data offer, in principle, the ability to identify functional inbreeding. Here, we present a heuristic method to identify haplotypes contained within a run of homozygosity (ROH) associated with reduced performance. Results are presented for simulated and swine data. The algorithm comprises 3 steps. Step 1 scans the genome based on marker windows of decreasing size and identifies ROH genotypes associated with an unfavorable phenotype. Within this stage, multiple aggregation steps reduce the haplotype to the smallest possible length. In step 2, the resulting regions are formally tested for significance with the use of a linear mixed model. Lastly, step 3 removes nested windows. The effect of the unfavorable haplotypes identified and their associated haplotype probabilities for a progeny of a given mating pair or an individual can be used to generate an inbreeding load matrix (ILM). Diagonals of ILM characterize the functional individual inbreeding load (IIL). We estimated the accuracy of predicting the phenotype based on IIL. We further compared the significance of the regression coefficient for IIL on phenotypes with genome-wide inbreeding metrics. We tested the algorithm using simulated scenarios (12 scenarios), combining different levels of linkage disequilibrium (LD) and number of loci impacting a quantitative trait. Additionally, we investigated 9 traits from 2 maternal purebred swine lines. In simulated data, as the LD in the population increased, the algorithm identified a greater proportion of the true unfavorable ROH effects. For example, the proportion of highly unfavorable true ROH effects identified rose from 32 to 41% for the low- to the high-LD scenario. In both simulated and real data, the haplotypes identified were contained within a much larger ROH (9.12-12.1 Mb). The IIL prediction accuracy was greater than 0 across all scenarios for simulated data (mean of 0.49 [95% confidence interval 0.47-0.52] for the high-LD scenario) and for nearly all swine traits (mean of 0.17 [SD 0.10]). On average, across simulated and swine data sets, the IIL regression coefficient was more closely related to progeny performance than any genome-wide inbreeding metric. A heuristic method was developed that identified ROH genotypes with reduced performance and characterized the combined effects of ROH genotypes within and across individuals.

The relationship between different measures of feed efficiency and feeding behavior traits in Duroc pigs

Utilization of feed in livestock species consists of a wide range of biological processes, and therefore, its efficiency can be expressed in various ways, including direct measurement, such as daily feed intake, as well as indicator measures, such as feeding behavior. Measuring feed efficiency is important to the swine industry, and its accuracy can be enhanced by using automated feeding systems, which record feed intake and associated feeding behavior of individual animals. Each automated feeder space is often shared among several pigs and therefore raises concerns about social interactions among pen mates with regard to feeding behavior. The study herein used a data set of 14,901 Duroc boars with individual records on feed intake, feeding behavior, and other off-test traits. These traits were modeled with and without the random spatial effect of Pen_Room, a concatenation of room and pen, or random social interaction among pen mates. The nonheritable spatial effect of common Pen-Room was observed for traits directly measuring feed intake and accounted for up to 13% of the total phenotypic variance in the average daily feeding rate. The social interaction effect explained larger proportions of phenotypic variation in all the traits studied, with the highest being 59% for ADFI in the group of feeding behaviors, 73% for residual feed intake (RFI; RFI4 and RFI6) in the feed efficiency traits, and 69% for intramuscular fat percentage in the off-test traits. After accounting for the social interaction effect, residual BW gain and RFI and BW gain (RIG) were found to have the heritability of 0.38 and 0.18, respectively, and had strong genetic correlations with growth and off-test traits. Feeding behavior traits were found to be moderately heritable, ranging from 0.14 (ADFI) to 0.52 (average daily occupation time), and some of them were strongly correlated with feed efficiency measures; for example, there was a genetic correlation of 0.88 between ADFI and RFI6. Our work suggested that accounting for the social common pen effect was important for estimating genetic parameters of traits recorded by the automated feeding system. Residual BW gain and RIG appeared to be two robust measures of feed efficiency. Feeding behavior measures are worth further investigation as indicators of feed efficiency.

A comparison of accuracy validation methods for genomic and pedigree-based predictions of swine litter size traits using Large White and simulated data

The objective of this study was to compare and determine the optimal validation method when comparing accuracy from single-step GBLUP (ssGBLUP) to traditional pedigree-based BLUP. Field data included six litter size traits. Simulated data included ten replicates designed to mimic the field data in order to determine the method that was closest to the true accuracy. Data were split into training and validation sets. The methods used were as follows: (i) theoretical accuracy derived from the prediction error variance (PEV) of the direct inverse (iLHS), (ii) approximated accuracies from the accf90(GS) program in the BLUPF90 family of programs (Approx), (iii) correlation between predictions and the single-step GEBVs from the full data set (GEBV_Full ), (iv) correlation between predictions and the corrected phenotypes of females from the full data set (Y_c ), (v) correlation from method iv divided by the square root of the heritability (Y_ch ) and (vi) correlation between sire predictions and the average of their daughters' corrected phenotypes (Y_cs ). Accuracies from iLHS increased from 0.27 to 0.37 (37%) in the Large White. Approximation accuracies were very consistent and close in absolute value (0.41 to 0.43). Both iLHS and Approx were much less variable than the corrected phenotype methods (ranging from 0.04 to 0.27). On average, simulated data showed an increase in accuracy from 0.34 to 0.44 (29%) using ssGBLUP. Both iLHS and Y_ch approximated the increase well, 0.30 to 0.46 and 0.36 to 0.45, respectively. GEBV_Full performed poorly in both data sets and is not recommended. Results suggest that for within-breed selection, theoretical accuracy using PEV was consistent and accurate. When direct inversion is infeasible to get the PEV, correlating predictions to the corrected phenotypes divided by the square root of heritability is adequate given a large enough validation data set.

The use of multiple imputation for the accurate measurements of individual feed intake by electronic feeders

Obtaining accurate individual feed intake records is the key first step in achieving genetic progress toward more efficient nutrient utilization in pigs. Feed intake records collected by electronic feeding systems contain errors (erroneous and abnormal values exceeding certain cutoff criteria), which are due to feeder malfunction or animal-feeder interaction. In this study, we examined the use of a novel data-editing strategy involving multiple imputation to minimize the impact of errors and missing values on the quality of feed intake data collected by an electronic feeding system. Accuracy of feed intake data adjustment obtained from the conventional linear mixed model (LMM) approach was compared with 2 alternative implementations of multiple imputation by chained equation, denoted as MI (multiple imputation) and MICE (multiple imputation by chained equation). The 3 methods were compared under 3 scenarios, where 5, 10, and 20% feed intake error rates were simulated. Each of the scenarios was replicated 5 times. Accuracy of the alternative error adjustment was measured as the correlation between the true daily feed intake (DFI; daily feed intake in the testing period) or true ADFI (the mean DFI across testing period) and the adjusted DFI or adjusted ADFI. In the editing process, error cutoff criteria are used to define if a feed intake visit contains errors. To investigate the possibility that the error cutoff criteria may affect any of the 3 methods, the simulation was repeated with 2 alternative error cutoff values. Multiple imputation methods outperformed the LMM approach in all scenarios with mean accuracies of 96.7, 93.5, and 90.2% obtained with MI and 96.8, 94.4, and 90.1% obtained with MICE compared with 91.0, 82.6, and 68.7% using LMM for DFI. Similar results were obtained for ADFI. Furthermore, multiple imputation methods consistently performed better than LMM regardless of the cutoff criteria applied to define errors. In conclusion, multiple imputation is proposed as a more accurate and flexible method for error adjustments in feed intake data collected by electronic feeders.

Genetic relationship between purebred and crossbred sow longevity

Background

The overall breeding objective for a nucleus swine selection program is to improve crossbred commercial performance. Most genetic improvement programs are based on an assumed high degree of positive relationship between purebred performance in a nucleus herd and their relatives’ crossbred performance in a commercial herd. The objective of this study was to examine the relationship between purebred and crossbred sow longevity performance. Sow longevity was defined as a binary trait with a success occurring if a sow remained in the herd for a certain number of parities and including the cumulative number born alive as a measure of reproductive success. Heritabilities, genetic correlations, and phenotypic correlations were estimated using THRGIBBS1F90.

Results

Results indicated little to no genetic correlations between crossbred and purebred reproductive traits. This indicates that selection for longevity or lifetime performance at the nucleus level may not result in improved longevity and lifetime performance at the crossbred level. Early parity performance was highly correlated with lifetime performance indicating that an indicator trait at an early parity could be used to predict lifetime performance. This would allow a sow to have her own record for the selection trait before she has been removed from the herd.

Conclusions

Results from this study aid in quantifying the relationship between purebred and crossbred performance and provide information for genetic companies to consider when developing a selection program where the objective is to improve crossbred sow performance. Utilizing crossbred records in a selection program would be the best way to improve crossbred sow productivity.

Variance component estimates for alternative litter size traits in swine

Litter size at d 5 (LS5) has been shown to be an effective trait to increase total number born (TNB) while simultaneously decreasing preweaning mortality. The objective of this study was to determine the optimal litter size day for selection (i.e., other than d 5). Traits included TNB, number born alive (NBA), litter size at d 2, 5, 10, 30 (LS2, LS5, LS10, LS30, respectively), litter size at weaning (LSW), number weaned (NW), piglet mortality at d 30 (MortD30), and average piglet birth weight (BirthWt). Litter size traits were assigned to biological litters and treated as a trait of the sow. In contrast, NW was the number of piglets weaned by the nurse dam. Bivariate animal models included farm, year-season, and parity as fixed effects. Number born alive was fit as a covariate for BirthWt. Random effects included additive genetics and the permanent environment of the sow. Variance components were plotted for TNB, NBA, and LS2 to LS30 using univariate animal models to determine how variances changed over time. Additive genetic variance was minimized at d 7 in Large White and at d 14 in Landrace pigs. Total phenotypic variance for litter size traits decreased over the first 10 d and then stabilized. Heritability estimates increased between TNB and LS30. Genetic correlations between TNB, NBA, and LS2 to LS29 with LS30 plateaued within the first 10 d. A genetic correlation with LS30 of 0.95 was reached at d 4 for Large White and at d 8 for Landrace pigs. Heritability estimates ranged from 0.07 to 0.13 for litter size traits and MortD30. Birth weight had an h of 0.24 and 0.26 for Large White and Landrace pigs, respectively. Genetic correlations among LS30, LSW, and NW ranged from 0.97 to 1.00. In the Large White breed, genetic correlations between MortD30 with TNB and LS30 were 0.23 and -0.64, respectively. These correlations were 0.10 and -0.61 in the Landrace breed. A high genetic correlation of 0.98 and 0.97 was observed between LS10 and NW for Large White and Landrace breeds, respectively. This would indicate that NW could possibly be used as an effective maternal trait, given a low level of cross-fostering, to avoid back calculating litter size traits from piglet records. Litter size at d 10 would be a compromise between gain in litter size at weaning and minimizing the potentially negative effects of the nurse dam and direct additive genetics of the piglets, as they are expected to increase throughout lactation.

Feed intake, average daily gain, feed efficiency, and real-time ultrasound traits in Duroc pigs: II. Genomewide association

Efficient use of feed resources has become a clear challenge for the U.S. pork industry as feed costs continue to be the largest variable expense. The availability of the Illumina Porcine60K BeadChip has greatly facilitated whole-genome association studies to identify chromosomal regions harboring genes influencing those traits. The current study aimed at identifying genomic regions associated with variation in feed efficiency and several production traits in a Duroc terminal sire population, including ADFI, ADG, feed conversion ratio, residual feed intake (RFI), real-time ultrasound back fat thickness (BF), ultrasound muscle depth, intramuscular fat content (IMF), birth weight (BW at birth), and weaning weight (BW at weaning). Single trait association analyses were performed using Bayes B models with 35,140 SNP on 18 autosomes after quality control. Significance of nonoverlapping 1-Mb length windows (n = 2,380) were tested across 3 QTL inference methods: posterior distribution of windows variances from Monte Carlo Markov Chain, naive Bayes factor, and nonparametric bootstrapping. Genes within the informative QTL regions for the traits were annotated. A region ranging from166 to 140 Mb (4-Mb length) on SSC 1, approximately 8 Mb upstream of the MC4R gene, was significantly associated with ADFI, ADG, and BF, where SOCS6 and DOK6 are proposed as the most likely candidate genes. Another region affecting BW at weaning was identified on SSC 4 (84-85 Mb), harboring genes previously found to influence both human and cattle height: PLAG1, CHCHD7, RDHE2 (or SDR16C5), MOS, RPS20, LYN, and PENK. No QTL were identified for RFI, IMF, and BW at birth. In conclusion, we have identified several genomic regions associated with traits affecting nutrient utilization that could be considered for future genomic prediction to improve feed utilization.

Feed intake, average daily gain, feed efficiency, and real-time ultrasound traits in Duroc pigs: I. Genetic parameter estimation and accuracy of genomic prediction

The efficiency of producing salable products in the pork industry is largely determined by costs associated with feed and by the amount and quality of lean meat produced. The objectives of this paper were 1) to explore heritability and genetic correlations for growth, feed efficiency, and real-time ultrasound traits using both pedigree and marker information and 2) to assess accuracy of genomic prediction for those traits using Bayes A prediction models in a Duroc terminal sire population. Body weight at birth (BW at birth) and weaning (BW at weaning) and real-time ultrasound traits, including back fat thickness (BF), muscle depth (MD), and intramuscular fat content (IMF), were collected on the basis of farm protocol. Individual feed intake and serial BW records of 1,563 boars obtained from feed intake recording equipment (FIRE; Osborne Industries Inc., Osborne, KS) were edited to obtain growth, feed intake, and feed efficiency traits, including ADG, ADFI, feed conversion ratio (FCR), and residual feed intake (RFI). Correspondingly, 1,047 boars were genotyped using the Illumina PorcineSNP60 BeadChip. The remaining 516 boars, as an independent sample, were genotyped with a low-density GGP-Porcine BeadChip and imputed to 60K. Magnitudes of heritability from pedigree analysis were moderate for growth, feed intake, and ultrasound traits (ranging from 0.44 ± 0.11 for ADG to 0.58 ± 0.09 for BF); heritability estimates were 0.32 ± 0.09 for FCR but only 0.10 ± 0.05 for RFI. Comparatively, heritability estimates using marker information by Bayes A models were about half of those from pedigree analysis, suggesting "missing heritability." Moderate positive genetic correlations between growth and feed intake (0.32 ± 0.05) and back fat (0.22 ± 0.04), as well as negative genetic correlations between growth and feed efficiency traits (-0.21 ± 0.08, -0.05 ± 0.07), indicate selection solely on growth traits may lead to an undesirable increase in feed intake, back fat, and reduced feed efficiency. Genetic correlations among growth, feed intake, and FCR assessed by a multiple-trait Bayes A model resulted in increased genetic correlation between ADG and ADFI, a negative correlation between ADFI and FCR, and a positive correlation between ADG and FCR. Accuracies of genomic prediction for the traits investigated, ranging from 9.4% for RFI to 36.5% for BF, were reported that might provide new insight into pig breeding and future selection programs using genomic information.

Estimation of genetic parameters for birth weight, preweaning mortality, and hot carcass weight of crossbred pigs

Genetic parameters for birth weight (BWT), preweaning mortality (PWM), and HCW were estimated for a crossbred pig population to determine if BWT could be used as an early predictor for later performances. Sire genetic effects for those traits were estimated to determine if early selection of purebred sires used in crossbreeding could be improved. Data were recorded from 1 commercial farm between 2008 and 2010. Data were from 24,376 crossbred pigs from Duroc sires and crossbred Large White × Landrace dams and included 24,376 BWT and PWM records and 13,029 HCW records. For the analysis, PWM was considered as a binary trait (0 for live or 1 for dead piglet at weaning). A multitrait threshold-linear animal model was used, with animal effect divided into sire genetic and dam effects; the dam effects included both genetic and environmental variation due to the absence of pedigree information for crossbred dams. Fixed effects were sex and parity for all traits, contemporary groups for BWT and HCW, and age at slaughter as a linear covariable for HCW. Random effects were sire additive genetic, dam, litter, and residual effects for all traits and contemporary group for PWM. Heritability estimates were 0.04 for BWT, 0.02 for PWM, and 0.12 for HCW. The ratio between sire genetic and total estimated variances was 0.01 for BWT and PWM and 0.03 for HCW. Dam and litter variances explained, respectively, 14% and 15% of total variance for BWT, 2% and 10% for PWM, and 3% and 8% for HCW. Genetic correlations were -0.52 between BWT and PWM, 0.55 between BWT and HCW, and -0.13 between PWM and HCW. Selection of purebred sires for higher BWT of crossbreds may slightly improve survival until weaning and final market weight at the commercial level.

Q-Band resonance Raman investigation of turnip cytochrome f and Rhodobacter capsulatus cytochrome c1

The results of a comprehensive Q-band resonance Raman investigation of cytochrome c1 and cytochrome f subunits of bc1 and b6f complexes are presented. Q-band excitation provides a particularly effective probe of the local heme environments of these species. The effects of protein conformation (particularly axial ligation) on heme structure and function were further investigated by comparison of spectra obtained from native subunits to those of a site directed c1 mutant (M183L) and various pH-dependent species of horse heart cytochrome c. In general, all species examined displayed variability in their axial amino acid ligation that suggests a good deal of flexibility in their hemepocket conformations. Surprisingly, the large scale protein rearrangements that accompany axial ligand replacement have little or no effect on macrocycle geometry in these species. This indicates the identity and/or conformation of the peptide linkage between the two cysteines that are covalently linked to the heme periphery may determine heme geometry.

Characterization of proteins utilized in the desulfurization of petroleum products by matrix-assisted laser desorption ionization time-of-flight mass spectrometry

Matrix-assisted laser desorption ionization time-of-flight mass spectrometry (MALDI/TOF/MS) with delayed extraction is utilized in linear, reflected-ion and post-source decay (PSD) modes to directly characterize enzymes being developed for use in a petroleum desulfurization process. The DNA sequence for the genes isolated from Rhodococcus sp. strain IGTS8 that produce three of the four enzymes under study had been previously reported with a discrepancy in residue assignments for one of the enzymes, dsz-C. The use of proteolytic digests followed by MALDI/TOF/MS with delayed extraction in the reflected-ion mode provided sequence-specific information with mass accuracies exceeding 40 ppm over a range of masses and signal-to-noise values. Peptide mapping of >80% of the residues was accomplished for all four proteins. The use of PSD established the true sequence for dsz-C, resolving the discrepancy in the literature. A posttranslational loss of N-terminal methionine was observed for each of the four proteins in linear MALDI/MS and was reconfirmed by peptide mapping for three of the proteins.

Molecular mechanisms of biocatalytic desulfurization of fossil fuels

The development of biocatalytic desulfurization of petroleum fractions may allow its use in place of conventional hydrodesulfurization (HDS). Dibenzothiophene (DBT) is representative of a broad range of sulfur heterocycles found in petroleum that are recalcitrant to desulfurization via HDS. Rhodococcus sp. strain IGTS8 has the ability to convert DBT to 2-hydroxybiphenyl (HBP) with the release of inorganic sulfur. The conversion of DBT to HBP is catalyzed by a multienzyme pathway consisting of two monooxygenases and a desulfinase. The final reaction catalyzed by the desulfinase appears to be the rate limiting step in the pathway. Each of the enzymes has been purified to homogeneity and their kinetic and physical properties studied. Neither monooxygenase has a tightly bound cofactor and each requires an NADH-FMN oxidoreductase for activity. An NADH-FMN oxidoreductase has been purified from Rhodococcus and is a protein of approximately 25,000 molecular weight with no apparent sequence homology to any other protein in the databases. We describe a unique sulfur acquisition system that Rhodococcus uses to obtain sulfur from very stable heterocyclic molecules.

Axial heme ligation in the cytochrome bc1 complexes of mitochondrial and photosynthetic membranes. A near-infrared magnetic circular dichroism and electron paramagnetic resonance study

The combination of EPR and low-temperature near-IR magnetic circular dichroism spectroscopies have been used to investigate the axial ligation of the cytochromes in the cytochrome bc1 complexes from bovine heart mitochondria, Rhodobacter capsulatus, Rhodobacter sphaeroides, and Rhodospirillum rubrum, and the purified cytochromes c1 from bovine heart mitochondria, Rb. capsulatus and Rb. sphaeroides. The possibility of axial ligation of cytochrome c1 by the amino terminus of the polypeptide was also assessed by acetylating the N-terminus of Rb. capsulatus cytochrome c1 and comparing the properties of the acetylated and unmodified samples. The results are consistent with bis-histidine axial ligation for the high- and low-potential b-type cytochromes and histidine/methionine axial ligation for the c1-type cytochrome in the intact cytochrome bc1 complexes. Purified samples of cytochrome c1 are mixtures of two forms, one with histidine/methionine and the other with bis-histidine axial ligation. The form with bis-histidine axial ligation is also assembled in the M183L mutant of the Rb. capsulatus cyt bc1 complex in which the methionine residue coordinating cyt c1 is replaced by a leucine. The bis-histidine form appears to be an artifact of dissociation of cytochrome c1 from the cytochrome bc1 complex and is greatly enhanced particularly in the bacterial cytochromes c1 by sample handling and the addition of 50% (v/v) ethylene glycol or glycerol.

Effects of cold storage on the subsequent structure and function of microvenous autografts

A study was undertaken to evaluate the structure and patency of cold stored rabbit femoral veins following reinsertion as autografts for 3 weeks. The periods of cold ischaemic storage were 1, 2, 3 and 4 weeks (n = 10/gp) and 6 and 10 weeks (n = 6/gp). All rabbits were subject to 3 operations under general anaesthesia. In the first, a 4 cm segment of left femoral vein was harvested and stored at 4 degrees C for the specified ischaemic interval. Following storage the graft was microsurgically reinserted at a second operation into the right femoral artery of the donor rabbit. Three weeks later, graft patency was assessed by surgical exploration and the graft processed for light and electron microscopy. Patency rates remained over 80% after 3 weeks in all groups except the 10 week storage group where only 1 of 6 (17%) grafts was patent at 3 weeks. In all groups normal vein structure was absent, being replaced by thin walled necrotic areas or by neointimal ingrowth. The excellent patency rates achieved indicate it is possible to cold preserve extra lengths of vein grafts harvested at initial operation for reuse should regrafting be necessary.

A look at adolescent treatment in a time of change

Multiple changes in the authors' inpatient adolescent unit--including unit mergers, increased patient acuity, and shorter length of stay--created the need and opportunity to implement a new Behavioral Milieu Program with two significant components: A positive Behavioral Points System and 15 nurse-led groups. The authors describe this new program, which provides a safe and effective treatment modality and promotes each adolescent's sense of responsibility for his own behavior and life. Dramatic outcomes include decreased use of mechanical restraints by 56% in the first three months of implementation and by 82% in the second three months. This program has applications for inpatient units, community residences, and schools, and can even be adapted for home use by parents.

Rhodobacter capsulatus contains a novel cb-type cytochrome c oxidase without a CuA center

The facultative phototrophic bacterium Rhodobacter capsulatus is capable of growth in a wide range of environmental conditions using a highly branched electron-transfer chain. During respiratory growth of this organism reducing equivalents are conveyed to oxygen via two terminal oxidases, previously called "cyt b410" (cytochrome c oxidase) and "cyt b260" (quinol oxidase). The cytochrome c oxidase was purified to homogeneity from a semiaerobically grown R. capsulatus strain. The purified enzyme consumes oxygen at a rate of 600 s-1, oxidizes reduced equine cyt c and R. capsulatus cyt c2, and has high sensitivity to cyanide. The complex is composed of three major polypeptides of apparent molecular masses 45, 32, and 28 kDa on SDS-PAGE. The 32- and 28-kDa proteins also stain with tetramethylbenzidine, indicating that they are c-type cytochromes. Partial amino acid sequences obtained from each of the subunits reveal significant homology to the fixN, fixO, and fixP gene products of Bradyrhizobium japonicum and Rhizobium meliloti. The reduced enzyme has an optical absorption spectrum with distinct features near 550 and 560 nm and an asymmetric Soret band centered at 418 nm, indicating the presence of both c- and b-type cytochromes. Two electrochemically distinct cyt c are apparent, with redox midpoint potentials (Em7) of 265 and 320 mV, while the low-spin cyt b has an Em7 value of 385 mV. The enzyme binds carbon monoxide, and the CO difference spectrum indicates that CO binds to a high-spin cyt b. Pyridine hemochrome and HPLC analyses suggest that the complex contains 1 mol of heme C to 1 mol of protoheme and that neither heme O nor heme A is present.(ABSTRACT TRUNCATED AT 250 WORDS)

Requirement of histidine 217 for ubiquinone reductase activity (Qi site) in the cytochrome bc1 complex

Folding models suggest that the highly conserved histidine 217 of the cytochrome b subunit from the cytochrome bc1 complex is close to the quinone reductase (Qi) site. This histidine (bH217) in the cytochrome b polypeptide of the photosynthetic bacterium Rhodobacter capsulatus has been replaced with three other residues, aspartate (D), arginine (R), and leucine (L). bH217D and bH217R are able to grow photoheterotrophically and contain active cytochrome bc1 complexes (60% of wild-type activity), whereas the bH217L mutant is photosynthetically incompetent and contains a cytochrome bc1 complex that has only 10% of the wild-type activity. Single-turnover flash-activated electron transfer experiments show that cytochrome bH is reduced via the Qo site with near native rates in the mutant strains but that electron transfer between cytochrome bH and quinone bound at the Qi site is greatly slowed. These results are consistent with redox midpoint potential (Em) measurements of the cytochrome b subunit hemes and the Qi site quinone. The Em values of cyt bL and bH are approximately the same in the mutants and wild type, although the mutant strains have a larger relative concentration of what may be the high-potential form of cytochrome bH, called cytochrome b150. However, the redox properties of the semiquinone at the Qi site are altered significantly. The Qi site semiquinone stability constant of bH217R is 10 times higher than in the wild type, while in the other two strains (bH217D and bH217L) the stability constant is much lower than in the wild type. Thus H217 appears to have major effects on the redox properties of the quinone bound at the Qi site. These data are incorporated into a suggestion that H217 forms part of the binding pocket of the Qi site in a manner reminiscent of the interaction between quinone bound at the Qb site and H190 of the L subunit of the bacterial photosynthetic reaction center.

The functional and structural effects of hypothermic storage on ischaemic arterial grafts

The effects of hypothermic ischaemia on blood vessels are unknown. This study aimed to determine the 3 week patency rate and the pathology of 9 experimental groups of hypothermically stored ischaemic arteries and one control group in a rabbit femoral artery model. Ischaemia times were 0 h, 24 h, 1, 2, 4, 6, 8 and 10 weeks (Groups 1-8). Patency was over 80% in all groups after 3 weeks reinsertion. Following reinsertion control grafts maintained normal arterial structure, but cellular degeneration had occurred in all ischaemic grafts and appeared complete after 4 weeks ischaemia. The graft connective tissue framework frequently remained intact. Repair was evident in central graft regions after 2 weeks ischaemia and 3 weeks reinsertion, but occurred only adjacent to the anastomosis in 4-10 week ischaemic arteries. Four week ischaemic arteries (Groups 9 and 10) reinserted for 6 and 12 weeks respectively exhibited near complete repair but patency dropped to 60% in the 12 week group.

The bc1 complexes of Rhodobacter sphaeroides and Rhodobacter capsulatus

Photosynthetic bacteria offer excellent experimental opportunities to explore both the structure and function of the ubiquinol-cytochrome c oxidoreductase (bc1 complex). In both Rhodobacter sphaeroides and Rhodobacter capsulatus, the bc1 complex functions in both the aerobic respiratory chain and as an essential component of the photosynthetic electron transport chain. Because the bc1 complex in these organisms can be functionally coupled to the photosynthetic reaction center, flash photolysis can be used to study electron flow through the enzyme and to examine the effects of various amino acid substitutions. During the past several years, numerous mutations have been generated in the cytochrome b subunit, in the Rieske iron-sulfur subunit, and in the cytochrome c1 subunit. Both site-directed and random mutagenesis procedures have been utilized. Studies of these mutations have identified amino acid residues that are metal ligands, as well as those residues that are at or near either the quinol oxidase (Qo) site or the quinol reductase (Qi) site. The postulate that these two Q-sites are located on opposite sides of the membrane is supported by these studies. Current research is directed at exploring the details of the catalytic mechanism, the nature of the subunit interactions, and the assembly of this enzyme.

Mutagenesis of methionine-183 drastically affects the physicochemical properties of cytochrome c1 of the bc1 complex of Rhodobacter capsulatus

Site-directed mutagenesis was used to investigate which of the highly conserved methionine residues (M183 and M205) provides the sixth axial ligand to the heme Fe in the cyt c1 subunit of the bc1 complex from the bacterium Rhodobacter capsulatus. These residues were changed to leucine (cM183L) and valine (cM205V). Two additional mutants were constructed, 1 in which a stop codon was inserted at M205 (cM205*) and the second in which 127 amino acids were deleted between the signal sequence and the putative C-terminal transmembrane alpha-helix (c delta SfuI). Only cM205V grew photosynthetically, and membranes isolated from this strain catalyzed quinol-dependent reduction of cyt c in amounts similar to that in a wild-type strain. Even though cM183L could not grow photosynthetically, it contained all the appropriate polypeptides and cofactors of the bc1 complex, as shown by SDS-PAGE and optical difference spectroscopy of intact membrane particles. Neither of the two deletion mutants contained a stable complex. Flash absorption spectroscopy using chromatophores showed no cytochrome c rereduction after oxidation by the reaction center in cM183L. The bc1 complex from each strain was isolated and characterized. Oxidation reduction midpoint potential titrations revealed that cyt c1 from cM183L had a dramatically shifted Em value (delta Em = -390 mV) compared with wild type and cM205V. While the optical absorption spectrum of cyt c1 from cM183L suggested that the c-type heme was low-spin, nonetheless it was able to react with the exogenous ligand carbon monoxide. The overall data support that M183, and not M205, is the sixth ligand to the heme Fe of cyt c1 of the bc1 complex.

Photochemical trapping of a bacteriopheophytin anion in site-specific reaction-center mutants from the photosynthetic bacterium Rhodobacter sphaeroides

The mutant YY in the reaction center of Rhodobacter sphaeroides, in which Phe181 on the L chain has been replaced by Tyr, and the double mutant FY, with Tyr210 on the M chain replaced by Phe and Phe181 on the L chain replaced by Tyr, have been constructed by site-directed mutagenesis. The studies described here were performed to complement a previous mutational analysis of mutant FF with Tyr210 replaced by Phe. Both new strains grow photoheterotrophically. The optical absorption spectra of reaction centers isolated from these mutants have band shifts attributable to the monomer bacteriochlorophylls in the vicinity of the substitutions. Photochemical trapping of the bacteriopheophytin anion (I-) indicates that the bacteriopheophytin on the B branch is reduced to a much greater extent in FF and FY as compared to YY and wild-type YF. Low temperature (77 K) absorption spectra clearly show that in the wild-type (YF) and YY reaction centers only the 545-nm-absorbing bacteriopheophytin is reduced while in the FF and FY reaction centers both the 535-nm and 545-nm-absorbing bacteriopheophytins are reduced. A simple kinetic analysis is used to explain these results. This analysis suggests that, in order for the observed trapping results to occur, a decrease in the 'cycling' time must take place, that is changes in the rate(s) of charge recombination must accompany the already known decrease in the forward electron transfer rate.

Secondary ischaemia in rabbit skin flaps: the roles played by thromboxane and free radicals

1. Biochemical mechanisms of ischaemia were investigated in rabbit skin flaps subjected to 2 h of primary ischaemia then, 24 h later, to 4 h of secondary ischaemia. During secondary ischaemia, flaps underwent either total ischaemia (arterial and venous blood supply occluded) or partial ischaemia (vein only occluded). Some of these flaps were treated at the time of reperfusion with the free-radical scavenger superoxide dismutase (EC 1.15.1.1) and/or the thromboxane synthetase inhibitor UK-38,485. 2. After 30 min of reperfusion, superoxide dismutase treatment significantly reduced blood thromboxane levels, elevated during ischaemia. Superoxide dismutase also reduced tissue levels of malonyldialdehyde and xanthine oxidase, indicators of free-radical damage, and restored the depleted tissue levels of superoxide dismutase. 3. UK-38,485 treatment failed to significantly alter any of these tissue free-radical parameters, although this agent significantly reduced blood thromboxane levels. 4. Combined superoxide dismutase plus UK-38,485 treatment was not significantly better than either treatment alone with respect to any parameter. 5. Partial ischaemia led to consistently higher levels of tissue free radicals and blood thromboxane than did total ischaemia. Thus partial ischaemia appears to result in greater free-radical damage than total ischaemia. 6. These results are consistent with the hypothesis that thromboxane acts as a mediator for free-radical damage in the ischaemic changes within the flap.

Resonance Raman characterization of Rhodobacter sphaeroides reaction centers bearing site-directed mutations at tyrosine M210

Resonance Raman (RR) spectroscopy and low-temperature absorption spectroscopy have been used to investigate the structural changes in the reaction centers (RCs) of Rhodobacter sphaeroides induced by site-specific mutations on the tyrosine (Y) M210 residue. RCs in which Y M210 has been genetically replaced with phenylalanine (F) or leucine (L) exhibit a 5-fold decrease in their primary electron-transfer kinetics (Finkele et al., 1990). The general similarity of RR spectra of the wild-type RCs as compared to those of the two mutant RCs indicates that no significant global structural changes occur upon mutation at the level of any of the six bacteriochlorin pigments. In the RR spectra of the two mutant RCs there is a conspicuous absence of contributions from the BPheM prosthetic group, which is interpreted in terms of a change in the resonance enhancement conditions of this chromophore. Low-temperature adsorption spectroscopy reveals marked shifts in the Qx absorption band of BPheM. This shift is interpreted as arising from a destabilization of the protein in the vicinity of BPheM and accounts for the change in resonance condition for this chromophore in its RR contributions. As well, there is a 3-nm red shift of the Qy absorption band of the BChls from 803 to 806 nm for the mutant RCs. Difference RR spectra yielding structural information concerning, selectively, the primary donor (P) indicate that the structure of the P binding pocket is conserved for these mutant RCs. The tyrosine M210 is not observed to be engaged in a hydrogen bond with either of the acetyl or keto carbonyls of P.

Interaction between thromboxane and free radical mechanisms in experimental ischaemic rabbit skin flaps

The possible relationship between increased blood levels of thromboxane (TXA2) and tissue levels of free radicals during ischaemia was investigated. Rabbit epigastric skin flaps were subjected to 4 h of body temperature ischaemia, then infused with either the TXA2 synthetase inhibitor UK-38,485, the free radical scavenger superoxide dismutase (SOD), or both immediately prior to reperfusion. After 30 min of reperfusion, increases in the tissue levels of xanthine oxidase (XO) and malonyldialdehyde (MDA), both of which are indices of free radical generation and decreases in the tissue levels of SOD were found. SOD treatment completely restored XO, MDA and SOD levels to normal, whereas UK-38,485 only partially improved all three parameters. None of these changes was statistically significant. Effluent blood thromboxane B2 (TXB2) levels from the flap increased significantly (P less than 0.01) after ischaemia and were reduced significantly by both UK-38,485 and SOD (P less than 0.05). Combined UK-38,485 and SOD treatment was no better than treatment with either agent alone. ATP levels and oedema, which decreased and increased respectively due to ischaemia, were not significantly altered by drug infusion. These results suggest that free radical damage may be related to TXA2-generated thrombosis in ischaemia/reperfusion injury.

Role of tyrosine M210 in the initial charge separation of reaction centers of Rhodobacter sphaeroides

Femtosecond spectroscopy was used in combination with site-directed mutagenesis to study the influence of tyrosine M210 (YM210) on the primary electron transfer in the reaction center of Rhodobacter sphaeroides. The exchange of YM210 to phenylalanine caused the time constant of primary electron transfer to increase from 3.5 +/- 0.4 ps to 16 +/- 6 ps while the exchange to leucine increased the time constant even more to 22 +/- 8 ps. The results suggest that tyrosine M210 is important for the fast rate of the primary electron transfer.

Initial characterization of site-directed mutants of tyrosine M210 in the reaction centre of Rhodobacter sphaeroides

A description of the properties of site-directed mutants of the reaction centre (RC) of Rhodobacter sphaeroides is presented. The residue tyrosine M210 (YM210) has been changed to phenylalanine (FM210) and leucine (LM210). Both mutants grew photoheterotrophically under conditions of high light but only the FM210 mutant grew under low light. Photobleaching spectra of chromatophores isolated from these mutants showed that the amount of functional RC was comparable to wild-type and that the spectral features were essentially unchanged. Shifts were observed in the absorption spectra in regions attributable to all the chromophores. An increase in intensity and a 3 nm red shift (from 803 to 806 nm) was observed in the Qy band of the monomer bacteriochlorophylls. A new extinction coefficient for the RC was determined at 806 nm (332 +/- 15 mM-1 cm-1). Linear dichroism (LD) spectra showed that there was no significant large scale change in the angles of the individual pigments relative to the C2 axis of symmetry. Cytochrome turnover assays were performed on isolated RC and light harvesting I complex (LH I)-RC (photosynthetic units, PSU) preparations. A turnover number of 120 cyt RC-1 s-1 was calculated for both the mutants while wild-type had a turnover number of 228 cyt RC-1 s-1. The cytochrome c2-mediated re-reduction kinetics of P+ were comparable to those observed in the wild-type. The half-time of charge recombination within the RC increased in the mutants to the wild-type (100 ms in the wild type, and 150 and 200 ms in FM210 and LM210, respectively).(ABSTRACT TRUNCATED AT 250 WORDS)

Spinach nitrite reductase: subunit composition and siroheme redox potential

The molecular weight of the spinach chloroplast, ferredoxin-dependent enzyme ferredoxin:nitrite oxidoreductase (EC 1.7.7.1) was shown to be 85,000 using gel filtration chromatography under nondenaturing conditions. In the presence of denaturants, either gel filtration chromatography or gel electrophoresis separated the enzyme into two subunits, with molecular weights of 61,000 and 24,000, respectively. Oxidation-reduction titrations of the siroheme prosthetic group of the putative native form (Mr 85,000) of the enzyme yielded a midpoint potential value of -305 mV, a value considerably more negative than the value of -30 mV obtained for siroheme in a modified, Mr 61,000 form of the enzyme.

Prosthetic group content and ligand binding properties of a spinach catalase

A recently discovered form of spinach catalase that contains both a novel heme and protoheme as prosthetic groups has been characterized using immunological and spectroscopic techniques. The enzyme appears to be a dimer of identical Mr 60,000 monomers. Extraction of the non-covalently bound prosthetic groups, followed by thin-layer chromatography of the extract, suggested that the novel heme contains four carboxylic acid side-chain groups. The resonance Raman spectrum of the resting enzyme indicates that the protoheme prosthetic group is five-coordinate and high-spin. The enzyme was shown to bind formate, azide and cyanide. Cyanide and azide binding to catalase are biphasic, suggesting the existence of two different binding sites for cyanide and azide in the enzyme. Results obtained from EPR and resonance Raman spectroscopies also support the hypothesis that two different ligand-binding sites are present in the enzyme. Western blots suggest that the Mr 60,000 peptide of the novel heme-containing catalase is similar or identical to that of a previously characterized, exclusively protoheme-containing, tetrameric catalase.

Complex formation between methylamine dehydrogenase and amicyanin from Paracoccus denitrificans

Two proteins isolated from Paracoccus denitrificans, the copper-containing electron carrier amicyanin and the pyrroloquinoline quinone-containing enzyme methylamine dehydrogenase, have been shown to form a complex. Complex formation between methylamine dehydrogenase and either oxidized or reduced amicyanin resulted in alterations in the absorbance spectrum of the pyrroloquinoline quinone prosthetic group of methylamine dehydrogenase. Binding of amicyanin to the enzyme exhibited positive cooperativity. Complex formation with methylamine dehydrogenase shifted the oxidation-reduction midpoint potential of amicyanin by 73 mV, from +294 to +221 mV, making electron transfer from amicyanin to cytochrome c551 (Em = +190 mV) thermodynamically possible.

Purification and characterization of two soluble cytochromes from the alkalophile Bacillus firmus RAB

A soluble cytochrome c and soluble cytochrome b were purified from the alkalophilic Bacillus firmus RAB. The cytochrome c, with an alpha band at 552 nm, had an apparent molecular weight of 16,500 and was acidic, with a pI of 3.4. At both pH 7.0 and 8.3, the midpoint potential of c-552 was +66 mV. Above pH 8.3, the cytochrome exhibited a pH-dependent decrease in midpoint potential. This property, among others, distinguished the cytochrome c-552 from other membrane-associated c-type cytochromes. The soluble cytochrome b, with an alpha band maximum at 558 nm, had a molecular weight of approx. 15,500 and was also an acidic protein, with a pI of 3.07. It exhibited a pH-independent midpoint potential of +28 mV.

Redox properties of the quinoprotein methylamine dehydrogenase from paracoccus denitrificans

Paracoccus denitrificans synthesizes a methylamine dehydrogenase that contains a covalently bound form of pyrroloquinoline quinone as a prosthetic group [Husain, M., & Davison, V.L. (1987) J. Bacteriol. 169, 1712-1717]. Anaerobic reductive titration of this enzyme with dithionite proceeded through a semiquinone intermediate with spectral properties quite distinct from those of the oxidized and reduced species. From these data the molar extinction coefficients were calculated at various wavelengths for the three redox states of this enzyme. The semiquinone was slowly reoxidized under aerobic conditions. The fully reduced enzyme was stable in the presence of oxygen and slowly reoxidized by ferricyanide. Reductive titration of methylamine dehydrogenase with methylamine proceeded directly to the fully reduced form of the enzyme without detectable formation of the semiquinone. Electrochemical titrations of the enzyme yielded an overall midpoint potential value for the two-electron couple (fully oxidized/fully reduced) of 100 +/- 4 mV and an n value of 2.15 +/- 0.15.

Measurement of the oxidation-reduction potentials of amicyanin and c-type cytochromes from Paracoccus denitrificans

The oxidation-reduction potentials of four periplasmic electron carrier proteins from Paracoccus denitrificans have been determined. Their midpoint potentials are: amicyanin, 294 +/- 6 mV; cytochrome c-550, 253 +/- 5 mV; cytochrome c-551i, 190 +/- 4 mV; and cytochrome c-553i, 148 +/- 5 mV. Although rapid amicyanin-mediated transfer of electrons from methylamine dehydrogenase to cytochrome c-551i was observed, reduced amicyanin did not reduce oxidized cytochrome c-551i in the absence of methylamine dehydrogenase.

An analysis of the reassembly of denatured creatine kinase from monkey brain

Creatine kinase isolated from monkey brain was characterized with respect to denaturation/inactivation and renaturation/reactivation/reassociation in order to determine the mechanism of reassembly. Enzyme unfolded in 8 M urea exhibits several characteristics of denatured protein: complete loss of enzymatic activity, decrease in intrinsic fluorescence with a red shift in the emission maximum and loss of circular dichroism at 220 nm. The renatured protein reassembles to its apparently native condition as judged by these criteria, but small differences of uncertain origin persist. Dependence of activity and fluorescence on denaturant concentration indicate that inactivation is more sensitive to urea than is unfolding; spectral changes at the intermediate urea concentrations suggest formation of aggregated protein. Upon dilution, enzyme previously exposed to 8 M urea for 40 min regains 70-80% native activity, independent of protein concentration over the range of 0.56-160 nM. Reactivation kinetics, measured using the assay mixture with and without trypsin, are independent of protein concentration, and are adequately described by a single rate constant, 3.2 X 10(-3) s-1 and 4.2 X 10(-3) s-1, respectively. Reactivation is completed 20-30 min after initiation of renaturation. Fluorescence changes during refolding are at least biphasic, exhibiting a rapid increase, then a slow decrease completed at approximately 15-20 min after initiating refolding. Reassociation is measured by competitive hybridization between dimerizing B subunits and M subunits to form MB heterodimer. The time dependent decay in heterodimer formation during competitive dimerization shows that reassociation is completed between 60 and 90 min after initiation of reassembly. These results indicate that the brain isozyme of creatine kinase, like the muscle form, is composed of subunits which do not require association for expression of catalytic activity. Furthermore, a comparison of spectral data and susceptibility to trypsin inactivation between the muscle and brain isozymes supports previous suggestions that in the native state, the brain isozyme is a conformationally looser, more open protein.

Purification and characterization of creatine kinase isozymes from the nurse shark Ginglymostoma cirratum

Creatine kinase from nurse shark brain and muscle has been purified to apparent homogeneity. In contrast to creatine kinases from most other vertebrate species, the muscle isozyme and the brain isozyme from nurse shark migrate closely in electrophoresis and, unusually, the muscle isozyme is anodal to the brain isozyme. The isoelectric points are 5.3 and 6.2 for the muscle and brain isozymes, respectively. The purified brain preparation also contains a second active protein with pI 6.0. The amino acid content of the muscle isozyme is compared with other isozymes of creatine kinase using the Metzger Difference Index as an estimation of compositional relatedness. All comparisons show a high degree of compositional similarity including arginine kinase from lobster muscle. The muscle isozyme is marginally more resistant to temperature inactivation than the brain isozyme; the muscle protein does not exhibit unusual stability towards high concentrations of urea. Kinetic analysis of the muscle isozyme reveals Michaelis constants of 1.6 mM MgATP, 12 mM creatine, 1.2 mM MgADP and 50 mM creatine phosphate. Dissociation constants for the same substrate from the binary and ternary enzyme-substrate complex do not differ significantly, indicating limited cooperatively in substrate binding. Enzyme activity is inhibited by small planar anions, most severely by nitrate. Shark muscle creatine kinase hybridizes in vitro with rabbit muscle or monkey brain creatine kinase; shark brain isozyme hybridizes with monkey brain or rabbit brain creatine kinase. Shark muscle and shark brain isozymes, under a wide range of conditions, failed to produce a detectable hybrid.