Fragmentation of contaminant and endogenous DNA in ancient samples determined by shotgun sequencing; prospects for human palaeogenomics

BACKGROUND: Despite the successful retrieval of genomes from past remains, the prospects for human palaeogenomics remain unclear because of the difficulty of distinguishing contaminant from endogenous DNA sequences. Previous sequence data generated on high-throughput sequencing platforms indicate that fragmentation of ancient DNA sequences is a characteristic trait primarily arising due to depurination processes that create abasic sites leading to DNA breaks. METHODOLOGY/PRINCIPALS FINDINGS: To investigate whether this pattern is present in ancient remains from a temperate environment, we have 454-FLX pyrosequenced different samples dated between 5,500 and 49,000 years ago: a bone from an extinct goat (Myotragus balearicus) that was treated with a depurinating agent (bleach), an Iberian lynx bone not subjected to any treatment, a human Neolithic sample from Barcelona (Spain), and a Neandertal sample from the El Sidrón site (Asturias, Spain). The efficiency of retrieval of endogenous sequences is below 1% in all cases. We have used the non-human samples to identify human sequences (0.35 and 1.4%, respectively), that we positively know are contaminants. CONCLUSIONS: We observed that bleach treatment appears to create a depurination-associated fragmentation pattern in resulting contaminant sequences that is indistinguishable from previously described endogenous sequences. Furthermore, the nucleotide composition pattern observed in 5' and 3' ends of contaminant sequences is much more complex than the flat pattern previously described in some Neandertal contaminants. Although much research on samples with known contaminant histories is needed, our results suggest that endogenous and contaminant sequences cannot be distinguished by the fragmentation pattern alone.     

Study of a major intermediate in the oxidative folding of leech carboxypeptidase inhibitor: contribution of the fourth disulfide bond

The oxidative folding pathway of leech carboxypeptidase inhibitor (LCI; four disulfide bonds) proceeds through the formation of two major intermediates (III-A and III-B) that contain three native disulfide bonds and act as strong kinetic traps in the folding process. The III-B intermediate lacks the Cys19-Cys43 disulfide bond that links the beta-sheet core with the alpha-helix in wild-type LCI. Here, an analog of this intermediate was constructed by replacing Cys19 and Cys43 with alanine residues. Its oxidative folding follows a rapid sequential flow through one, two, and three disulfide species to reach the native form; the low accumulation of two disulfide intermediates and three disulfide (scrambled) isomers accounts for a highly efficient reaction. The three-dimensional structure of this analog, alone and in complex with carboxypeptidase A (CPA), was determined by X-ray crystallography at 2.2A resolution. Its overall structure is very similar to that of wild-type LCI, although the residues in the region adjacent to the mutation sites show an increased flexibility, which is strongly reduced upon binding to CPA. The structure of the complex also demonstrates that the analog and the wild-type LCI bind to the enzyme in the same manner, as expected by their inhibitory capabilities, which were similar for all enzymes tested. Equilibrium unfolding experiments showed that this mutant is destabilized by approximately 1.5 kcal mol(-1) (40%) relative to the wild-type protein. Together, the data indicate that the fourth disulfide bond provides LCI with both high stability and structural specificity.     

The activation of 5-HT receptors in prefrontal cortex enhances dopaminergic activity

Atypical antipsychotics show preferential 5-HT 2A versus dopamine (DA) D2 receptor affinity. At clinical doses, they fully occupy cortical 5-HT2 receptors, which suggests a strong relationship with their therapeutic action. Half of the pyramidal neurones in the medial prefrontal cortex (mPFC) express 5-HT 2A receptors. Also, neurones excited through 5-HT 2A receptors project to the ventral tegmental area (VTA). We therefore hypothesized that prefrontal 5-HT 2A receptors can modulate DA transmission through excitatory mPFC-VTA inputs. In this study we used single unit recordings to examine the responses of DA neurones to local (in the mPFC) and systemic administration of the 5-HT 2A/2C agonist 1-[2,5-dimethoxy-4-iodophenyl-2-aminopropane] (DOI). Likewise, using microdialysis, we examined DA release in the mPFC and VTA (single/dual probe) in response to prefrontal and systemic drug administration. The local (in the mPFC) and systemic administration of DOI increased the firing rate and burst firing of DA neurones and DA release in the VTA and mPFC. The increase in VTA DA release was mimicked by the electrical stimulation of the mPFC. The effects of DOI were reversed by M100907 and ritanserin. These results indicate that the activity of VTA DA neurones is under the excitatory control of 5-HT 2A receptors in the mPFC. These observations may help in the understanding of the therapeutic action of atypical antipsychotics.     

Molecular mechanisms involved in the adenosine A and A receptor-induced neuronal differentiation in neuroblastoma cells and striatal primary cultures

Adenosine A1 receptors (A1Rs) and adenosine A(2A) receptors (A(2A)Rs) are the major mediators of the neuromodulatory actions of adenosine in the brain. In the striatum A1Rs and A(2A)Rs are mainly co-localized in the GABAergic striatopallidal neurons. In this paper we show that agonist-induced stimulation of A1Rs and A(2A)Rs induces neurite outgrowth processes in the human neuroblastoma cell line SH-SY5Y and also in primary cultures of striatal neuronal precursor cells. The kinetics of adenosine-mediated neuritogenesis was faster than that triggered by retinoic acid. The triggering of the expression of TrkB neurotrophin receptor and the increase of cell number in the G1 phase by the activation of adenosine receptors suggest that adenosine may participate in early steps of neuronal differentiation. Furthermore, protein kinase C (PKC) and extracellular regulated kinase-1/2 (ERK-1/2) are involved in the A1R- and A(2A)R-mediated effects. Inhibition of protein kinase A (PKA) activity results in a total inhibition of neurite outgrowth induced by A(2A)R agonists but not by A1R agonists. PKA activation is therefore necessary for A(2A)R-mediated neuritogenesis. Co-stimulation does not lead to synergistic effects thus indicating that the neuritogenic effects of adenosine are mediated by either A1 or A(2A) receptors depending upon the concentration of the nucleoside. These results are relevant to understand the mechanisms by which adenosine receptors modulate neuronal differentiation and open new perspectives for considering the use of adenosine agonists as therapeutic agents in diseases requiring neuronal repair.     

NMR structural characterization and computational predictions of the major intermediate in oxidative folding of leech carboxypeptidase inhibitor

The III-A intermediate constitutes the major rate-determining step in the oxidative folding of leech carboxypeptidase inhibitor (LCI). In this work, III-A has been directly purified from the folding reaction and structurally characterized by NMR spectroscopy. This species, containing three native disulfides, displays a highly native-like structure; however, it lacks some secondary structure elements, making it more flexible than native LCI. III-A represents a structurally determined example of a disulfide-insecure intermediate; direct oxidation of this species to the fully native protein seems to be restricted by the burial of its two free cysteine residues inside a native-like structure. We also show that theoretical approaches based on topological constraints predict with good accuracy the presence of this folding intermediate. Overall, the derived results suggest that, as it occurs with non-disulfide bonded proteins, native-like interactions between segments of secondary structure rather than the crosslinking of disulfide bonds direct the folding of LCI.     

A common haplotype associated with the Basque 2362AG --> TCATCT mutation in the muscular calpain-3 gene

Limb-girdle muscular dystrophy type 2A (LGMD2A) is caused by any of over 150 mutations in the calpain-3 (CAPN3) gene. Of those, 2362AG --> TCATCT is particularly prevalent in Basque patients, and this mutation was hypothesized to have arisen in the Basque Country. To explore the natural history of this mutation, we genotyped 65 Basque and non-Basque patients with LGMD2A who carry the 2362AG --> TCATCT mutation for four microsatellites within or flanking the gene. A particular haplotype was found in three-fourths of the patients and was assumed to be ancestral. From the average number of recombinations and mutations accumulated from this ancestral haplotype, the age of the 2362AG ----> TCATCT mutation was estimated to be 50 generations (i.e., 1,250 years), which is more recent than the Paleolithic Basque heritage. The subsequent spread of the 2362AG --> TCATCT mutation can be related to gene flow out of the Basque Country, even across a cultural border.     

Secondary binding site of the potato carboxypeptidase inhibitor. Contribution to its structure, folding, and biological properties

The contribution of each residue of the potato carboxypeptidase inhibitor (PCI) secondary binding site to the overall properties of this protein has been examined using alanine-scanning mutagenesis. Structural and enzymatic studies, performed on a series of PCI mutants, demonstrate that the proper positioning of the primary site for efficient binding and inhibition of carboxypeptidase A is significantly dependent on such a secondary contact region. The aromatic residues in this region play a key role in the stabilization of the PCI-enzyme complex, whereas polar residues contribute little to this task. A comparative study of the oxidative folding of these PCI mutants has been carried out using the disulfide quenching approach. The data, together with the structural characterization of some of these mutants, clearly indicate that noncovalent forces drive the refolding of this small disulfide-rich protein at the reshuffling stage, the rate-limiting step of the process. Moreover, it reveals that by introducing new noncovalent intramolecular contacts in PCI, we may create more stable variants, which also show improved folding efficiency. Taken together, the collected results clarify the folding determinants of the primary and secondary binding sites of PCI and their contribution to the inhibition of the carboxypeptidase, providing clues about PCI evolution and knowledge for its biotechnological redesign.     

Role of kinetic intermediates in the folding of leech carboxypeptidase inhibitor

The oxidative folding and reductive unfolding pathways of leech carboxypeptidase inhibitor (LCI; four disulfides) have been characterized in this work by structural and kinetic analysis of the acid-trapped folding intermediates. The oxidative folding of reduced and denatured LCI proceeds rapidly through a sequential flow of 1-, 2-, 3-, and 4-disulfide (scrambled) species to reach the native form. Folding intermediates of LCI comprise two predominant 3-disulfide species (designated as III-A and III-B) and a heterogeneous population of scrambled isomers that consecutively accumulate along the folding reaction. Our study reveals that forms III-A and III-B exclusively contain native disulfide bonds and correspond to stable and partially structured species that interconvert, reaching an equilibrium prior to the formation of the scrambled isomers. Given that these intermediates act as kinetic traps during the oxidative folding, their accumulation is prevented when they are destabilized, thus leading to a significant acceleration of the folding kinetics. III-A and III-B forms appear to have both native disulfides bonds and free thiols similarly protected from the solvent; major structural rearrangements through the formation of scrambled isomers are required to render native LCI. The reductive unfolding pathway of LCI undergoes an apparent all-or-none mechanism, although low amounts of intermediates III-A and III-B can be detected, suggesting differences in protection against reduction among the disulfide bonds. The characterization of III-A and III-B forms shows that the former intermediate structurally and functionally resembles native LCI, whereas the III-B form bears more resemblance to scrambled isomers.     

A Simple Geometrical Pattern for the Branching Distribution of the Bronchial Tree, Useful to Estimate Optimality Departures

The design of the bronchial tree has largely been proposed as a model of optimal design from a physical-functional perspective. However, the distributive function of the airway may be more related to a geometrical than a physical problem. The bronchial tree must distribute a three dimensional volume of inspired air on a two dimensional alveolar surface, included in a limited volume. It is thus valid to ask whether an optimal bronchial tree from a physical perspective is also optimum from a geometrical point of view. In this paper we generate a simple geometric model for the branching pattern of the bronchial tree, deducing relationships that permit estimation of the departures from the geometrical optimum of each bifurcation. We also, for comparative purposes, estimate the departures from the physical optimum. From the geometrical assumptions: i) a symmetrical dichotomic fractal design, ii) with minimum volume and iii) maximum dispersion of the terminal points; and several simulations we suggest that the optimality is characterized by a bifurcation angle theta approximately 60 degrees and a length reduction scale gamma = (1/2)(1/3) = 0.7937. We propose distances from the physical and geometrical optimality defined as Euclidean distances from the expected optima. We show how the advanced relationships and the distances can be used to estimate departures from the optimality in bronchographs of four species. We found lower physical and geometrical departures in the distal zone than those of the proximal zones, as well as lower physical than geometrical departures from optimality.     

Table of Contents of Volume 52 - 2004

Table of Contents

Table of Contents of Volume 52 - 2004