Whole genome sequencing analysis of Plasmodium vivax using whole genome capture

Background: Malaria caused by Plasmodium vivax is an experimentally neglected severe disease with a substantial burden on human health. Because of technical limitations, little is known about the biology of this important human pathogen. Whole genome analysis methods on patient-derived material are thus likely to have a substantial impact on our understanding of P. vivax pathogenesis and epidemiology. For example, it will allow study of the evolution and population biology of the parasite, allow parasite transmission patterns to be characterized, and may facilitate the identification of new drug resistance genes. Because parasitemias are typically low and the parasite cannot be readily cultured, on-site leukocyte depletion of blood samples is typically needed to remove human DNA that may be 1000X more abundant than parasite DNA. These features have precluded the analysis of archived blood samples and require the presence of laboratories in close proximity to the collection of field samples for optimal pre-cryopreservation sample preparation. Results: Here we show that in-solution hybridization capture can be used to extract P. vivax DNA from human contaminating DNA in the laboratory without the need for on-site leukocyte filtration. Using a whole genome capture method, we were able to enrich P. vivax DNA from bulk genomic DNA from less than 0.5% to a median of 55% (range 20%-80%). This level of enrichment allows for efficient analysis of the samples by whole genome sequencing and does not introduce any gross biases into the data. With this method, we obtained greater than 5X coverage across 93% of the P. vivax genome for four P. vivax strains from Iquitos, Peru, which is similar to our results using leukocyte filtration (greater than 5X coverage across 96% of the genome). Conclusion: The whole genome capture technique will enable more efficient whole genome analysis of P. vivax from a larger geographic region and from valuable archived sample collections.

SIS: a program to generate draft genome sequence scaffolds for prokaryotes

Background: Decreasing costs of DNA sequencing have made prokaryotic draft genome sequences increasingly common. A contig scaffold is an ordering of contigs in the correct orientation. A scaffold can help genome comparisons and guide gap closure efforts. One popular technique for obtaining contig scaffolds is to map contigs onto a reference genome. However, rearrangements that may exist between the query and reference genomes may result in incorrect scaffolds, if these rearrangements are not taken into account. Large-scale inversions are common rearrangement events in prokaryotic genomes. Even in draft genomes it is possible to detect the presence of inversions given sufficient sequencing coverage and a sufficiently close reference genome. Results: We present a linear-time algorithm that can generate a set of contig scaffolds for a draft genome sequence represented in contigs given a reference genome. The algorithm is aimed at prokaryotic genomes and relies on the presence of matching sequence patterns between the query and reference genomes that can be interpreted as the result of large-scale inversions; we call these patterns inversion signatures. Our algorithm is capable of correctly generating a scaffold if at least one member of every inversion signature pair is present in contigs and no inversion signatures have been overwritten in evolution. The algorithm is also capable of generating scaffolds in the presence of any kind of inversion, even though in this general case there is no guarantee that all scaffolds in the scaffold set will be correct. We compare the performance of SIS, the program that implements the algorithm, to seven other scaffold-generating programs. The results of our tests show that SIS has overall better performance. Conclusions: SIS is a new easy-to-use tool to generate contig scaffolds, available both as stand-alone and as a web server. The good performance of SIS in our tests adds evidence that large-scale inversions are widespread in prokaryotic genomes.

Non-classical HLA-E gene variability in Brazilians: a nearly invariable locus surrounded by the most variable genes in the human genome

The non-classical human leukocyte antigen (HLA) class I genes present a very low rate of variation. So far, only 10 HLA-E alleles encoding three proteins have been described, but only two are frequently found in worldwide populations. Because of its historical background, Brazilians are very suitable for population genetic studies. Therefore, 104 bone marrow donors from Brazil were evaluated for HLA-E exons 14. Seven variation sites were found, including two known single nucleotide polymorphisms (SNPs) at positions +424 and +756 and five new SNPs at positions +170 (intron 1), +1294 (intron 3), +1625, +1645 and +1857 (exon 4). Haplotyping analysis did show eight haplotypes, three of them known as E*01:01:01, E*01:03:01 and E*01:03:02:01 and five HLA-E new alleles that carry the new variation sites. The HLA-E*01:01:01 allele was the predominant haplotype (62.50%), followed by E*01:03:02:01 (24.52%). Selective neutrality tests have disclosed an interesting pattern of selective pressures in which balancing selection is probably shaping allele frequency distributions at an SNP at exon 3 (codon 107), sequence diversity at exon 4 and the non-coding regions is facing significant purifying pressure. Even in an admixed population such as the Brazilian one, the HLA-E locus is very conserved, presenting few polymorphic SNPs in the coding region.

Nuclear astrophysics: nucleosynthesis in the Universe

Nuclear astrophysics is a relatively young science; it is about half a century old. It is a multidisciplinary subject, since it combines nuclear physics with astrophysics and observations in astronomy. It also addresses fundamental issues in astrobiology through the formation of elements, in particular those required for a carbon-based life. In this paper, a rapid overview of nucleosynthesis is given, mainly from the point of view of nuclear physics. A short historical introduction is followed by the definition of the relevant nuclear parameters, such as nuclear reaction cross sections, astrophysical S-factors, the energy range defined by the Gamow peak and reaction rates. The different astrophysical scenarios that are the sites of nucleosynthesis, and different processes, cycles and chains that are responsible for the building of complex nuclei from the elementary hydrogen nuclei are then briefly described. Received 28 February 2012, accepted 5 April 2012, first published online 9 May 2012

Pollination systems in Pogonieae (Orchidaceae: Vanilloideae): A hypothesis of evolution among reward and rewardless flowers

The tribe Pogonieae of Vanilloideae (Orchidaceae) consists of six genera, including Pogoniopsis, a mycoheterotrophic taxon with morphological characteristics distinct from the remaining of the tribe. A hypothesis about the phylogeny of the tribe was inferred, involving all currently recognized genera, based on isolated and combined sequence data of 5.8S, 18S and 26S (nrDNA) regions using parsimony and Bayesian analyses. Phylogenetic analyses show that inclusion of Pogoniopsis turns the tribe Pogonieae paraphyletic. All analyses reveal that Pogoniopsis is closely related to members of Epidendroideae. The pantropical Vanilla is monophyletic if Dictyophyllaria is assumed as synonym of Vanilla. Members of Pogonieae are pollinated by several groups of solitary and social bees, two pollination systems being recognized: reward-producing and deceptive. The molecular phylogeny suggests that ancestrals related to Pogonieae gave rise to two evolutionary lines: a tropical one with reward production of flowers, and a predominantly temperate regions invading line with deceptive flowers. Reward-producing flowers characterize the South and Central American clade (=Cleistes), while deceptive pollination is prominent in the clade that includes North American-Asiatic taxa plus the Amazonian genus Duckeella. (C) 2012 Elsevier GmbH. All rights reserved.

Anomalous centrality evolution of two-particle angular correlations from Au-Au collisions at root s(NN)=62 and 200 GeV

We present two-dimensional (2D) two-particle angular correlations measured with the STAR detector on relative pseudorapidity eta and azimuth phi for charged particles from Au-Au collisions at root s(NN) = 62 and 200 GeV with transverse momentum p(t) >= 0.15 GeV/c, vertical bar eta vertical bar <= 1, and 2 pi in azimuth. Observed correlations include a same-side (relative azimuth <pi/2) 2D peak, a closely related away-side azimuth dipole, and an azimuth quadrupole conventionally associated with elliptic flow. The same-side 2D peak and away-side dipole are explained by semihard parton scattering and fragmentation (minijets) in proton-proton and peripheral nucleus-nucleus collisions. Those structures follow N-N binary-collision scaling in Au-Au collisions until midcentrality, where a transition to a qualitatively different centrality trend occurs within one 10% centrality bin. Above the transition point the number of same-side and away-side correlated pairs increases rapidly relative to binary-collision scaling, the eta width of the same-side 2D peak also increases rapidly (eta elongation), and the phi width actually decreases significantly. Those centrality trends are in marked contrast with conventional expectations for jet quenching in a dense medium. The observed centrality trends are compared to perturbative QCD predictions computed in HIJING, which serve as a theoretical baseline, and to the expected trends for semihard parton scattering and fragmentation in a thermalized opaque medium predicted by theoretical calculations and phenomenological models. We are unable to reconcile a semihard parton scattering and fragmentation origin for the observed correlation structure and centrality trends with heavy-ion collision scenarios that invoke rapid parton thermalization. If the collision system turns out to be effectively opaque to few-GeV partons the present observations would be inconsistent with the minijet picture discussed here. DOI: 10.1103/PhysRevC.86.064902

Genetic variability in mitochondrial and nuclear genes of Larus dominicanus (Charadriiformes, Laridae) from the Brazilian coast

Several phylogeographic studies of seabirds have documented low genetic diversity that has been attributed to bottleneck events or individual capacity for dispersal. Few studies have been done in seabirds on the Brazilian coast and all have shown low genetic differentiation on a wide geographic scale. The Kelp Gull is a common species with a wide distribution in the Southern Hemisphere. In this study, we used mitochondrial and nuclear markers to examine the genetic variability of Kelp Gull populations on the Brazilian coast and compared this variability with that of sub-Antarctic island populations of this species. Kelp Gulls showed extremely low genetic variability for nnitochondrial markers (cytb and ATPase) and high diversity for a nuclear locus (intron 7 of the beta-fibrinogen). The intraspecific evolutionary history of Kelp Gulls showed that the variability found in intron 7 of the beta-fibrinogen gene was compatible with the variability expected under neutral evolution but suggested an increase in population size during the last 10,000 years. However, none of the markers revealed evidence of a bottleneck population. These findings indicate that the recent origin of Kelp Gulls is the main explanation for their nuclear diversity, although selective pressure on the mtDNA of this species cannot be discarded.

Molecular characterization of a mariner-like element in the Atta sexdens rubropilosa genome

Mobile elements are widely present in eukaryotic genomes. They are repeated DNA segments that are able to move from one locus to another within the genome. They are divided into two main categories, depending on their mechanism of transposition, involving RNA (class I) or DNA (class II) molecules. The mariner-like elements are class II transposons. They encode their own transposase, which is necessary and sufficient for transposition in the absence of host factors. They are flanked by a short inverted terminal repeat and a TA dinucleotide target site, which is duplicated upon insertion. The transposase consists of two domains, an N-terminal inverted terminal repeat binding domain and a C-terminal catalytic domain. We identified a transposable element with molecular characteristics of a mariner-like element in Atta sexdens rubropilosa genome. Identification started from a PCR with degenerate primers and queen genomic DNA templates, with which it was possible to amplify a fragment with mariner transposable-element homology. Phylogenetic analysis demonstrated that this element belongs to the mauritiana subfamily of mariner-like elements and it was named Asmar1. We found that Asmar1 is homologous to a transposon described from another ant, Messor bouvieri. The predicted transposase sequence demonstrated that Asmar1 has a truncated transposase ORF. This study is part of a molecular characterization of mobile elements in the Atta spp genome. Our finding of mariner-like elements in all castes of this ant could be useful to help understand the dynamics of mariner-like element distribution in the Hymenoptera.

Evolution of a dense neutrino gas in matter and electromagnetic field

We describe the system of massive Weyl fields propagating in a background matter and interacting with an external electromagnetic field. The interaction with an electromagnetic field is due to the presence of anomalous magnetic moments. To canonically quantize this system first we develop the classical field theory treatment of Weyl spinors in frames of the Hamilton formalism which accounts for the external fields. Then, on the basis of the exact solution of the wave equation for a massive Weyl field in a background matter we obtain the effective Hamiltonian for the description of spin-flavor oscillations of Majorana neutrinos in matter and a magnetic field. Finally, we incorporate in our analysis the neutrino self-interaction which is essential when the neutrino density is sufficiently high. We also discuss the applicability of our results for the studies of collective effects in spin-flavor oscillations of supernova neutrinos in a dense matter and a strong magnetic field. (C) 2011 Elsevier B.V. All rights reserved.

Evolution of pi(0) Suppression in Au plus Au Collisions from root s(NN)=39 to 200 GeV

Neutral-pion pi(0) spectra were measured at midrapidity (vertical bar y vertical bar < 0.35) in Au + Au collisions at root s(NN) = 39 and 62.4 GeV and compared with earlier measurements at 200 GeV in a transverse-momentum range of 1 < p(T) < 10 GeV/c. The high-p(T) tail is well described by a power law in all cases, and the powers decrease significantly with decreasing center-of-mass energy. The change of powers is very similar to that observed in the corresponding spectra for p + p collisions. The nuclear modification factors (RAA) show significant suppression, with a distinct energy, centrality, and p(T) dependence. Above p(T) = 7 GeV/c, R-AA is similar for root sNN = 62.4 and 200 GeV at all centralities. Perturbative-quantum-chromodynamics calculations that describe R-AA well at 200 GeV fail to describe the 39 GeV data, raising the possibility that, for the same p(T) region, the relative importance of initial-state effects and soft processes increases at lower energies. The p(T) range where pi(0) spectra in central Au + Au collisions have the same power as in p + p collisions is approximate to 5 and 7 GeV/c for root sNN = 200 and 62.4 GeV, respectively. For the root sNN = 39 GeV data, it is not clear whether such a region is reached, and the x(T) dependence of the x(T)-scaling power-law exponent is very different from that observed in the root sNN = 62 and 200 GeV data, providing further evidence that initial-state effects and soft processes mask the in-medium suppression of hardscattered partons to higher p(T) as the collision energy decreases.

Novel organization of the mitochondrial genome in the deep-sea coral, Madrepora oculata (Hexacorallia, Scleractinia, Oculinidae) and its taxonomic implications

Madrepora is one of the most ecologically important genera of reef-building scleractinians in the deep sea, occurring from tropical to high-latitude regions. Despite this, the taxonomic affinities and relationships within the genus Madrepora remain unclear. To clarify these issues, we sequenced the mitochondrial (mt) genome of the most widespread Madrepora species, M. oculata, and compared this with data for other scleractinians. The architecture of the M. oculara mt genome was very similar to that of other scleractinians, except for a novel gene rearrangement affecting only cox2 and cox3. This pattern of gene organization was common to four geographically distinct M. oculata individuals as well as the congeneric species M. minutiseptum, but was not shared by other genera that are closely related on the basis of cox1 sequence analysis nor other oculinids, suggesting that it might be unique to Madrepora. (C) 2012 Elsevier Inc. All rights reserved.

Profiling Three-Dimensional Nuclear Telomeric Architecture of Myelodysplastic Syndromes and Acute Myeloid Leukemia Defines Patient Subgroups

Purpose: Myelodysplastic syndromes (MDS) are a group of disorders characterized by cytopenias, with a propensity for evolution into acute myeloid leukemias (AML). This transformation is driven by genomic instability, but mechanisms remain unknown. Telomere dysfunction might generate genomic instability leading to cytopenias and disease progression. Experimental Design: We undertook a pilot study of 94 patients with MDS (56 patients) and AML (38 patients). The MDS cohort consisted of refractory cytopenia with multilineage dysplasia (32 cases), refractory anemia (12 cases), refractory anemia with excess of blasts (RAEB) 1 (8 cases), RAEB2 (1 case), refractory anemia with ring sideroblasts (2 cases), and MDS with isolated del(5q) (1 case). The AML cohort was composed of AML-M4 (12 cases), AML-M2 (10 cases), AML-M5 (5 cases), AML-M0 (5 cases), AML-M1 (2 cases), AML-M4eo (1 case), and AML with multidysplasia-related changes (1 case). Three-dimensional quantitative FISH of telomeres was carried out on nuclei from bone marrow samples and analyzed using TeloView. Results: We defined three-dimensional nuclear telomeric profiles on the basis of telomere numbers, telomeric aggregates, telomere signal intensities, nuclear volumes, and nuclear telomere distribution. Using these parameters, we blindly subdivided the MDS patients into nine subgroups and the AML patients into six subgroups. Each of the parameters showed significant differences between MDS and AML. Combining all parameters revealed significant differences between all subgroups. Three-dimensional telomeric profiles are linked to the evolution of telomere dysfunction, defining a model of progression from MDS to AML. Conclusions: Our results show distinct three-dimensional telomeric profiles specific to patients with MDS and AML that help subgroup patients based on the severity of telomere dysfunction highlighted in the profiles. Clin Cancer Res; 18(12); 3293-304. (C) 2012 AACR.

Different patterns of evolution for duplicated DNA repair genes in bacteria of the Xanthomonadales group

Abstract Background DNA repair genes encode proteins that protect organisms against genetic damage generated by environmental agents and by-products of cell metabolism. The importance of these genes in life maintenance is supported by their high conservation, and the presence of duplications of such genes may be easily traced, especially in prokaryotic genomes. Results The genome sequences of two Xanthomonas species were used as the basis for phylogenetic analyses of genes related to DNA repair that were found duplicated. Although 16S rRNA phylogenetic analyses confirm their classification at the basis of the gamma proteobacteria subdivision, differences were found in the origin of the various genes investigated. Except for lexA, detected as a recent duplication, most of the genes in more than one copy are represented by two highly divergent orthologs. Basically, one of such duplications is frequently positioned close to other gamma proteobacteria, but the second is often positioned close to unrelated bacteria. These orthologs may have occurred from old duplication events, followed by extensive gene loss, or were originated from lateral gene transfer (LGT), as is the case of the uvrD homolog. Conclusions Duplications of DNA repair related genes may result in redundancy and also improve the organisms' responses to environmental challenges. Most of such duplications, in Xanthomonas, seem to have arisen from old events and possibly enlarge both functional and evolutionary genome potentiality.

Theoretical studies of Membrane Proteins : Properties, Prediction Methods and Genome-wide analysis

Membrane proteins are a large and important class of proteins. They are responsible for several of the key functions in a living cell, e.g. transport of nutrients and ions, cell-cell signaling, and cell-cell adhesion. Despite their importance it has not been possible to study their structure and organization in much detail because of the difficulty to obtain 3D structures. In this thesis theoretical studies of membrane protein sequences and structures have been carried out by analyzing existing experimental data. The data comes from several sources including sequence databases, genome sequencing projects, and 3D structures. Prediction of the membrane spanning regions by hydrophobicity analysis is a key technique used in several of the studies. A novel method for this is also presented and compared to other methods. The primary questions addressed in the thesis are: What properties are common to all membrane proteins? What is the overall architecture of a membrane protein? What properties govern the integration into the membrane? How many membrane proteins are there and how are they distributed in different organisms? Several of the findings have now been backed up by experiments. An analysis of the large family of G-protein coupled receptors pinpoints differences in length and amino acid composition of loops between proteins with and without a signal peptide and also differences between extra- and intracellular loops. Known 3D structures of membrane proteins have been studied in terms of hydrophobicity, distribution of secondary structure and amino acid types, position specific residue variability, and differences between loops and membrane spanning regions. An analysis of several fully and partially sequenced genomes from eukaryotes, prokaryotes, and archaea has been carried out. Several differences in the membrane protein content between organisms were found, the most important being the total number of membrane proteins and the distribution of membrane proteins with a given number of transmembrane segments. Of the properties that were found to be similar in all organisms, the most obvious is the bias in the distribution of positive charges between the extra- and intracellular loops. Finally, an analysis of homologues to membrane proteins with known topology uncovered two related, multi-spanning proteins with opposite predicted orientations. The predicted topologies were verified experimentally, providing a first example of "divergent topology evolution".

A high energy perspective of the co-evolution of black holes and their host galaxies

Thanks to the Chandra and XMM–Newton surveys, the hard X-ray sky is now probed down to a flux limit where the bulk of the X-ray background is almost completely resolved into discrete sources, at least in the 2–8 keV band. Extensive programs of multiwavelength follow-up observations showed that the large majority of hard X–ray selected sources are identified with Active Galactic Nuclei (AGN) spanning a broad range of redshifts, luminosities and optical properties. A sizable fraction of relatively luminous X-ray sources hosting an active, presumably obscured, nucleus would not have been easily recognized as such on the basis of optical observations because characterized by “peculiar” optical properties. In my PhD thesis, I will focus the attention on the nature of two classes of hard X-ray selected “elusive” sources: those characterized by high X-ray-to-optical flux ratios and red optical-to-near-infrared colors, a fraction of which associated with Type 2 quasars, and the X-ray bright optically normal galaxies, also known as XBONGs. In order to characterize the properties of these classes of elusive AGN, the datasets of several deep and large-area surveys have been fully exploited. The first class of “elusive” sources is characterized by X-ray-to-optical flux ratios (X/O) significantly higher than what is generally observed from unobscured quasars and Seyfert galaxies. The properties of well defined samples of high X/O sources detected at bright X–ray fluxes suggest that X/O selection is highly efficient in sampling high–redshift obscured quasars. At the limits of deep Chandra surveys (∼10−16 erg cm−2 s−1), high X/O sources are generally characterized by extremely faint optical magnitudes, hence their spectroscopic identification is hardly feasible even with the largest telescopes. In this framework, a detailed investigation of their X-ray properties may provide useful information on the nature of this important component of the X-ray source population. The X-ray data of the deepest X-ray observations ever performed, the Chandra deep fields, allows us to characterize the average X-ray properties of the high X/O population. The results of spectral analysis clearly indicate that the high X/O sources represent the most obscured component of the X–ray background. Their spectra are harder (G ∼ 1) than any other class of sources in the deep fields and also of the XRB spectrum (G ≈ 1.4). In order to better understand the AGN physics and evolution, a much better knowledge of the redshift, luminosity and spectral energy distributions (SEDs) of elusive AGN is of paramount importance. The recent COSMOS survey provides the necessary multiwavelength database to characterize the SEDs of a statistically robust sample of obscured sources. The combination of high X/O and red-colors offers a powerful tool to select obscured luminous objects at high redshift. A large sample of X-ray emitting extremely red objects (R−K >5) has been collected and their optical-infrared properties have been studied. In particular, using an appropriate SED fitting procedure, the nuclear and the host galaxy components have been deconvolved over a large range of wavelengths and ptical nuclear extinctions, black hole masses and Eddington ratios have been estimated. It is important to remark that the combination of hard X-ray selection and extreme red colors is highly efficient in picking up highly obscured, luminous sources at high redshift. Although the XBONGs do not present a new source population, the interest on the nature of these sources has gained a renewed attention after the discovery of several examples from recent Chandra and XMM–Newton surveys. Even though several possibilities were proposed in recent literature to explain why a relatively luminous (LX = 1042 − 1043erg s−1) hard X-ray source does not leave any significant signature of its presence in terms of optical emission lines, the very nature of XBONGs is still subject of debate. Good-quality photometric near-infrared data (ISAAC/VLT) of 4 low-redshift XBONGs from the HELLAS2XMMsurvey have been used to search for the presence of the putative nucleus, applying the surface-brightness decomposition technique. In two out of the four sources, the presence of a nuclear weak component hosted by a bright galaxy has been revealed. The results indicate that moderate amounts of gas and dust, covering a large solid angle (possibly 4p) at the nuclear source, may explain the lack of optical emission lines. A weak nucleus not able to produce suffcient UV photons may provide an alternative or additional explanation. On the basis of an admittedly small sample, we conclude that XBONGs constitute a mixed bag rather than a new source population. When the presence of a nucleus is revealed, it turns out to be mildly absorbed and hosted by a bright galaxy.