A whole-genome assembly of the domestic cow, Bos taurus

Background

The genome of the domestic cow, Bos taurus, was sequenced using a mixture of hierarchical and whole-genome shotgun sequencing methods.

Results

We have assembled the 35 million sequence reads and applied a variety of assembly improvement techniques, creating an assembly of 2.86 billion base pairs that has multiple improvements over previous assemblies: it is more complete, covering more of the genome; thousands of gaps have been closed; many erroneous inversions, deletions, and translocations have been corrected; and thousands of single-nucleotide errors have been corrected. Our evaluation using independent metrics demonstrates that the resulting assembly is substantially more accurate and complete than alternative versions.

Conclusions

By using independent mapping data and conserved synteny between the cow and human genomes, we were able to construct an assembly with excellent large-scale contiguity in which a large majority (approximately 91%) of the genome has been placed onto the 30 B. taurus chromosomes. We constructed a new cow-human synteny map that expands upon previous maps. We also identified for the first time a portion of the B. taurus Y chromosome.     

Partial inhibition of sarcoplasmic reticulum ca release evokes long-lasting ca release events in ventricular myocytes: role of luminal ca in termination of ca release

In cardiac myocytes, local sarcoplasmic reticulum (SR) Ca depletion during Ca sparks is believed to play an important role in the termination of SR Ca release. We tested whether decreasing the rate of SR Ca depletion by partially inhibiting SR Ca release channels (ryanodine receptors) delays Ca spark termination. In permeabilized cat ventricular myocytes, 0.7 mM tetracaine caused almost complete Ca spark inhibition followed by a recovery significantly below control level. The recovery was associated with increased SR Ca load and increased Ca spark duration. Additionally, SR Ca release events lasting several hundred milliseconds occurred consistently. These events had a significantly lower initial Ca release flux followed by a stable plateau, indicating delayed release termination and maintained SR Ca load. Increasing SR Ca load (without inhibiting SR Ca release rate) or decreasing SR Ca release rate (without increasing SR Ca load) both induced only a small increase in spark duration. These results show that the combination of decreased release flux and increased SR Ca load has synergistic effects and exerts major changes on the termination of Ca release events. Long-lasting Ca release events may originate from highly interconnected release junctions where Ca diffusion from neighboring sites partially compensates Ca depletion, thereby delaying SR Ca-dependent termination. Eventually, these events terminate by luminal Ca-independent mechanisms, such as inactivation, adaptation, or stochastic attrition.     

Structure, dynamics, and branch migration of a DNA Holliday junction: a single-molecule fluorescence and modeling study

The Holliday junction (HJ) is a central intermediate of various genetic processes, including homologous and site-specific DNA recombination and DNA replication. Elucidating the structure and dynamics of HJs provides the basis for understanding the molecular mechanisms of these genetic processes. Our previous single-molecule fluorescence studies led to a model according to which branch migration is a stepwise process consisting of consecutive migration and folding steps. These data led us to the conclusion that one hop can be more than 1 basepair (bp); moreover, we hypothesized that continuous runs over the entire sequence homology (5 bp) can occur. Direct measurements of the dependence of the fluorescence resonance energy transfer (FRET) value on the donor-acceptor (D-A) distance are required to justify this model and are the major goal of this article. To accomplish this goal, we performed single-molecule FRET experiments with a set of six immobile HJ molecules with varying numbers of bps between fluorescent dyes placed on opposite arms. The designs were made in such a way that the distances between the donor and acceptor were equal to the distances between the dyes formed upon 1-bp migration hops of a HJ having 10-bp homology. Using these designs, we confirmed our previous hypothesis that the migration of the junction can be measured with bp accuracy. Moreover, the FRET values determined for each acceptor-donor separation corresponded very well to the values for the steps on the FRET time trajectories, suggesting that each step corresponds to the migration of the branch at a defined depth. We used the dependence of the FRET value on the D-A distance to measure directly the size for each step on the FRET time trajectories. These data showed that one hop is not necessarily 1 bp. The junction is able to migrate over several bps, detected as one hop and confirming our model. The D-A distances extracted from the FRET properties of the immobile junctions formed the basis for modeling the HJ structures. The composite data fit a partially opened, side-by-side model with adjacent double-helical arms slightly kinked at the four-way junction and the junction as a whole adopting a global X-shaped form that mimics the coaxially stacked-X structure implicated in previous solution studies.     

Non-native glyceraldehyde-3-phosphate dehydrogenase can be an intrinsic component of amyloid structures

Interactions between different forms of glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and amyloid-beta peptide (1-42) were investigated by direct (surface plasmon resonance) and indirect (kinetics of spontaneous and GroEL/S-assisted reactivation of denatured GAPDH) methods. It was demonstrated that non-native forms of GAPDH obtained by different ways (cold denaturation, oxidation of the enzyme, and its unfolding in guanidine hydrochloride) efficiently bind to soluble amyloid-beta peptide (1-42) yielding a stable complex. Native tetrameric GAPDH does not interact with soluble amyloid-beta peptide (1-42), neither non-native forms of GAPDH interact with aggregated amyloid-beta peptide (1-42). The results suggest that non-native GAPDH species can be involved in the formation of amyloid structures during Alzheimer's disease, binding to soluble amyloid-beta peptide (1-42).     

Identification of open reading frames unique to a select agent: Ralstonia solanacearum race 3 biovar 2

An 8x draft genome was obtained and annotated for Ralstonia solanacearum race 3 biovar 2 (R3B2) strain UW551, a United States Department of Agriculture Select Agent isolated from geranium. The draft UW551 genome consisted of 80,169 reads resulting in 582 contigs containing 5,925,491 base pairs, with an average 64.5% GC content. Annotation revealed a predicted 4,454 protein coding open reading frames (ORFs), 43 tRNAs, and 5 rRNAs; 2,793 (or 62%) of the ORFs had a functional assignment. The UW551 genome was compared with the published genome of R. solanacearum race 1 biovar 3 tropical tomato strain GMI1000. The two phylogenetically distinct strains were at least 71% syntenic in gene organization. Most genes encoding known pathogenicity determinants, including predicted type III secreted effectors, appeared to be common to both strains. A total of 402 unique UW551 ORFs were identified, none of which had a best hit or >45% amino acid sequence identity with any R. solanacearum predicted protein; 16 had strong (E < 10(-13)) best hits to ORFs found in other bacterial plant pathogens. Many of the 402 unique genes were clustered, including 5 found in the hrp region and 38 contiguous, potential prophage genes. Conservation of some UW551 unique genes among R3B2 strains was examined by polymerase chain reaction among a group of 58 strains from different races and biovars, resulting in the identification of genes that may be potentially useful for diagnostic detection and identification of R3B2 strains. One 22-kb region that appears to be present in GMI1000 as a result of horizontal gene transfer is absent from UW551 and encodes enzymes that likely are essential for utilization of the three sugar alcohols that distinguish biovars 3 and 4 from biovars 1 and 2.     

New binding site on common molecular scaffold provides HERG channel specificity of scorpion toxin BeKm-1

The scorpion toxin BeKm-1 is unique among a variety of known short scorpion toxins affecting potassium channels in its selective action on ether-a-go-go-related gene (ERG)-type channels. BeKm-1 shares the common molecular scaffold with other short scorpion toxins. The toxin spatial structure resolved by NMR consists of a short alpha-helix and a triple-stranded antiparallel beta-sheet. By toxin mutagenesis study we identified the residues that are important for the binding of BeKm-1 to the human ERG K+ (HERG) channel. The most critical residues (Tyr-11, Lys-18, Arg-20, Lys-23) are located in the alpha-helix and following loop whereas the "traditional" functional site of other short scorpion toxins is formed by residues from the beta-sheet. Thus the unique location of the binding site of BeKm-1 provides its specificity toward the HERG channel.     

Improving tolerance of Candida antarctica lipase B towards irreversible thermal inactivation through directed evolution

To expand the functionality of lipase B from Candida antarctica (CALB) we have used directed evolution to create CALB mutants with improved resistance towards irreversible thermal inactivation. Two mutants, 23G5 and 195F1, were generated with over a 20-fold increase in half-life at 70 degrees C compared with the wild-type CALB (WT-CALB). The increase in half-life was attributed to a lower propensity of the mutants to aggregate in the unfolded state and to an improved refolding. The first generation mutant, 23G5, obtained by error-prone PCR, had two amino acid mutations, V210I and A281E. The second generation mutant, 195F1, derived from 23G5 by error-prone PCR, had one additional mutation, V221D. Amino acid substitutions at positions 221 and 281 were determined to be critical for lipase stability, while the residue at position 210 had only a marginal effect. The catalytic efficiency of the mutants with p-nitrophenyl butyrate and 6,8-difluoro-4-methylumbelliferyl octanoate was also found to be superior to that of WT-CALB.     

The analysis of structure-anticancer and antiviral activity relationships for macrocyclic pyridinophanes and their analogues on the basis of 4D QSAR models (simplex representation of molecular structure)

A new 4D-QSAR approach has been considered. For all investigated molecules the 3D structural models have been created and the set of conformers (fourth dimension) have been used. Each conformer is represented as a system of different simplexes (tetratomic fragments of fixed structure, chirality and symmetry). The investigation of influence of molecular structure of macrocyclic pyridinophanes, their analogues and certain other compounds on anticancer and antiviral (anti-influenza, antiherpes and antiadenovirus) activity has been carried out by means of the 4D-QSAR. Statistic characteristics for QSAR of PLS (partial least squares) models are satisfactory (R = 0.92-0.97; CVR = 0.63-0.83). Molecular fragments increasing and decreasing biological activity were defined. This information may be useful for design, and direct synthesis of novel anticancer and antiviral agents.