Sequence characterization and comparative analysis of the gastrotropin gene in buffalo (Bubalus bubalis)

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Pilot survey of expressed sequence tags (ESTs) from the asexual blood stages of Plasmodium vivax in human patients

Abstract Background Plasmodium vivax is the most widely distributed human malaria, responsible for 70–80 million clinical cases each year and large socio-economical burdens for countries such as Brazil where it is the most prevalent species. Unfortunately, due to the impossibility of growing this parasite in continuous in vitro culture, research on P. vivax remains largely neglected. Methods A pilot survey of expressed sequence tags (ESTs) from the asexual blood stages of P. vivax was performed. To do so, 1,184 clones from a cDNA library constructed with parasites obtained from 10 different human patients in the Brazilian Amazon were sequenced. Sequences were automatedly processed to remove contaminants and low quality reads. A total of 806 sequences with an average length of 586 bp met such criteria and their clustering revealed 666 distinct events. The consensus sequence of each cluster and the unique sequences of the singlets were used in similarity searches against different databases that included P. vivax, Plasmodium falciparum, Plasmodium yoelii, Plasmodium knowlesi, Apicomplexa and the GenBank non-redundant database. An E-value of <10-30 was used to define a significant database match. ESTs were manually assigned a gene ontology (GO) terminology Results A total of 769 ESTs could be assigned a putative identity based upon sequence similarity to known proteins in GenBank. Moreover, 292 ESTs were annotated and a GO terminology was assigned to 164 of them. Conclusion These are the first ESTs reported for P. vivax and, as such, they represent a valuable resource to assist in the annotation of the P. vivax genome currently being sequenced. Moreover, since the GC-content of the P. vivax genome is strikingly different from that of P. falciparum, these ESTs will help in the validation of gene predictions for P. vivax and to create a gene index of this malaria parasite.

Parallel Layout Construction Algorithm for Irregular Shape Packing Problems

Cutting and packing problems arise in a variety of industries, including garment, wood and shipbuilding. Irregular shape packing is a special case which admits irregular items and is much more complex due to the geometry of items. In order to ensure that items do not overlap and no item from the layout protrudes from the container, the collision free region concept was adopted. It represents all possible translations for a new item to be inserted into a container with already placed items. To construct a feasible layout, collision free region for each item is determined through a sequence of Boolean operations over polygons. In order to improve the speed of the algorithm, a parallel version of the layout construction was proposed and it was applied to a simulated annealing algorithm used to solve bin packing problems. Tests were performed in order to determine the speed improvement of the parallel version over the serial algorithm

Models and Algorihtm for the Optimization of Real-World Routing and Logistics Problems

Logistics involves planning, managing, and organizing the flows of goods from the point of origin to the point of destination in order to meet some requirements. Logistics and transportation aspects are very important and represent a relevant costs for producing and shipping companies, but also for public administration and private citizens. The optimization of resources and the improvement in the organization of operations is crucial for all branches of logistics, from the operation management to the transportation. As we will have the chance to see in this work, optimization techniques, models, and algorithms represent important methods to solve the always new and more complex problems arising in different segments of logistics. Many operation management and transportation problems are related to the optimization class of problems called Vehicle Routing Problems (VRPs). In this work, we consider several real-world deterministic and stochastic problems that are included in the wide class of the VRPs, and we solve them by means of exact and heuristic methods. We treat three classes of real-world routing and logistics problems. We deal with one of the most important tactical problems that arises in the managing of the bike sharing systems, that is the Bike sharing Rebalancing Problem (BRP). We propose models and algorithms for real-world earthwork optimization problems. We describe the 3DP process and we highlight several optimization issues in 3DP. Among those, we define the problem related to the tool path definition in the 3DP process, the 3D Routing Problem (3DRP), which is a generalization of the arc routing problem. We present an ILP model and several heuristic algorithms to solve the 3DRP.

Comparative human-horse sequence analysis of the CYP3A subfamily gene cluster

Cytochrome P450 enzymes (CYP450s) represent a superfamily of haem-thiolate proteins. CYP450s are most abundant in the liver, a major site of drug metabolism, and play key roles in the metabolism of a variety of substrates, including drugs and environmental contaminants. Interaction of two or more different drugs with the same enzyme can account for adverse effects and failure of therapy. Human CYP3A4 metabolizes about 50% of all known drugs, but little is known about the orthologous CYP450s in horses. We report here the genomic organization of the equine CYP3A gene cluster as well as a comparative analysis with the human CYP3A gene cluster. The equine CYP450 genes of the 3A family are located on ECA 13 between 6.97-7.53 Mb, in a region syntenic to HSA 7 99.05-99.35 Mb. Seven potential, closely linked equine CYP3A genes were found, in contrast to only four genes in the human genome. RNA was isolated from an equine liver sample, and the approximately 1.5-kb coding sequence of six CYP3A genes could be amplified by RT-PCR. Sequencing of the RT-PCR products revealed numerous hitherto unknown single nucleotide polymorphisms (SNPs) in these six CYP3A genes, and one 6-bp deletion compared to the reference sequence (EquCab2.0). The presence of the variants was confirmed in a sample of genomic DNA from the same horse. In conclusion, orthologous genes for the CYP3A family exist in horses, but their number differs from those of the human CYP3A gene family. CYP450 genes of the same family show high homology within and between mammalian species, but can be highly polymorphic.

Sequence analysis of the porcine IFNAR1 and IFNGR2 genes

A porcine BAC clone harboring the tightly linked IFNAR1 and IFNGR2 genes was identified by comparative analysis of the publicly available porcine BAC end sequences. The complete 168,835 bp insert sequence of this clone was determined. Sequence comparisons of the genomic sequence with EST sequences from public databases were performed and allowed a detailed annotation of the IFNAR1 and IFNGR2 genes. The analyzed genes showed a conserved genomic organization with their known mammalian orthologs, however the sequence conservation of these genes across species was relatively low. In addition to the IFNAR1 and IFNGR2 genes, which were completely sequenced, the analyzed BAC clone also contained parts of an orphan gene encoding a putative transmembrane protein (TMEM50B). In contrast to the IFNAR1 and IFNGR2 genes the sequence conservation of the TMEM50B gene across different mammalian species was extremely high.

Sequence analysis of a 212 kb defensin gene cluster on ECA 27q17

Defensins are a family of evolutionary ancient antimicrobial peptides consisting of three sub-families: alpha-, beta- and theta-defensins. This investigation was focused on the genomic characterization of equine beta-defensins and the investigation of the potential clustering of beta-defensin genes in the equine genome. Six genomic BAC clones were isolated from the CHORI-241 library and one of these was mapped by FISH to ECA 27q17. This location was confirmed by RH-mapping. The contiguous 212 kb sequence of this clone was determined. Sequence analysis revealed the identification of ten pseudogenes and nine genes, six of which were highly homologous to human beta-defensin DEFB4. Clustering of the beta-defensin genes was confirmed and the order of the genes on the analyzed BAC was related to the corresponding defensin cluster on HSA 8. The knowledge about the sequence and the genomic structure of the equine beta-defensin genes will improve the classification of different paralogous defensin genes and is a prerequisite for subsequent functional studies. Additionally, the first alpha-defensin-like sequence outside the groups of primates, lagomorphs and rodents (glires) was identified.

Genome sequence, comparative analysis, and population genetics of the domestic horse

We report a high-quality draft sequence of the genome of the horse (Equus caballus). The genome is relatively repetitive but has little segmental duplication. Chromosomes appear to have undergone few historical rearrangements: 53% of equine chromosomes show conserved synteny to a single human chromosome. Equine chromosome 11 is shown to have an evolutionary new centromere devoid of centromeric satellite DNA, suggesting that centromeric function may arise before satellite repeat accumulation. Linkage disequilibrium, showing the influences of early domestication of large herds of female horses, is intermediate in length between dog and human, and there is long-range haplotype sharing among breeds.

Optimal sequence timing of CT angiography and perfusion CT in patients with stroke

OBJECTIVE Standard stroke CT protocols start with non-enhanced CT followed by perfusion-CT (PCT) and end with CTA. We aimed to evaluate the influence of the sequence of PCT and CTA on quantitative perfusion parameters, venous contrast enhancement and examination time to save critical time in the therapeutic window in stroke patients. METHODS AND MATERIALS Stroke CT data sets of 85 patients, 47 patients with CTA before PCT (group A) and 38 with CTA after PCT (group B) were retrospectively analyzed by two experienced neuroradiologists. Parameter maps of cerebral blood flow, cerebral blood volume, time to peak and mean transit time and contrast enhancements (arterial and venous) were compared. RESULTS Both readers rated contrast of brain-supplying arteries to be equal in both groups (p=0.55 (intracranial) and p=0.73 (extracranial)) although the extent of venous superimposition of the ICA was rated higher in group B (p=0.04). Quantitative perfusion parameters did not significantly differ between the groups (all p>0.18), while the extent of venous superimposition of the ICA was rated higher in group B (p=0.04). The time to complete the diagnostic CT examination was significantly shorter for group A (p<0.01). CONCLUSION Performing CTA directly after NECT has no significant effect on PCT parameters and avoids venous preloading in CTA, while examination times were significantly shorter.

Implicit task sequence learning in patients with Parkinson's disease, frontal lesions and amnesia: The critical role of fronto–striatal loops

The purpose of this study was to investigate the role of the fronto–striatal system for implicit task sequence learning. We tested performance of patients with compromised functioning of the fronto–striatal loops, that is, patients with Parkinson's disease and patients with lesions in the ventromedial or dorsolateral prefrontal cortex. We also tested amnesic patients with lesions either to the basal forebrain/orbitofrontal cortex or to thalamic/medio-temporal regions. We used a task sequence learning paradigm involving the presentation of a sequence of categorical binary-choice decision tasks. After several blocks of training, the sequence, hidden in the order of tasks, was replaced by a pseudo-random sequence. Learning (i.e., sensitivity to the ordering) was assessed by measuring whether this change disrupted performance. Although all the patients were able to perform the decision tasks quite easily, those with lesions to the fronto–striatal loops (i.e., patients with Parkinson's disease, with lesions in the ventromedial or dorsolateral prefrontal cortex and those amnesic patients with lesions to the basal forebrain/orbitofrontal cortex) did not show any evidence of implicit task sequence learning. In contrast, those amnesic patients with lesions to thalamic/medio-temporal regions showed intact sequence learning. Together, these results indicate that the integrity of the fronto–striatal system is a prerequisite for implicit task sequence learning.

OSL and TL dating of the Middle Stone Age sequence at Diepkloof Rock Shelter (South Africa): a clarification

Diepkloof Rock Shelter offers an exceptional opportunity to study the onset and evolution of both Still Bay (SB) and Howiesons Poort (HP) techno-complexes. However, previous age estimates based on luminescence dating of burnt quartzites (Tribolo et al., 2009) and of sediments (Jacobs et al., 2008) were not in agreement. Here, we present new luminescence ages for 17 rock samples (equivalent dose estimated with a SAR-ITL protocol instead of classical MAAD-TL) as well as for 5 sediment samples (equivalent dose estimated with SAR-single grain OSL protocol) and an update of the 22 previous age estimates for burnt lithics (modified calibration and beta dose estimates). While a good agreement between the rock and sediment ages is obtained, these estimates are still significantly older than those reported by Jacobs et al. (2008). After our own analyses of the sediment from Diepkloof, it is suspected that these authors did not correctly chose the parameters for the equivalent dose determination, leading to an underestimate of the equivalent doses, and thus of the ages. From bottom to top, the mean ages are 100 ± 10 ka for stratigraphic unit (SU) Noël and 107 ± 11 ka for SU Mark (uncharacterized Lower MSA), 100 ± 10 ka for SU Lynn-Leo (Pre-SB type Lynn), 109 ± 10 ka for SUs Kim-Larry (SB), 105 ± 10 ka for SUs Kerry-Kate and 109 ± 10 ka for SU Jess (Early HP), 89 ± 8 ka for SU Jude (MSA type Jack), 77 ± 8 ka for SU John, 85 ± 9 ka for SU Fox, 83 ± 8 ka for SU Fred and 65 ± 8 ka for SU OB5 (Intermediate HP), 52 ± 5 ka for SUs OB2-4 (Late HP). This chronology, together with the technological analyses, greatly modifies the current chrono-cultural model regarding the SB and the HP and has important archaeological implications. Indeed, SB and HP no longer appear as short-lived techno-complexes with synchronous appearances for each and restricted to Oxygen Isotopic Stage (OIS) 4 across South Africa, as suggested by Jacobs et al. (2008, 2012). Rather, the sequence of Diepkloof supports a long chronology model with an early appearance of both SB and HP in the first half of OIS 5 and a long duration of the HP into OIS 3. These new dates imply that different technological traditions coexisted during OIS 5 and 4 in southern Africa and that SB and HP can no longer be considered as horizon markers.

Antigenic variation in Lyme disease spirochetes by segmental recombination of VMP-like sequence cassettes

Lyme disease is a multisystemic disorder caused by tick-borne infection of humans or other mammalian hosts with Borrelia burgdorferi. If untreated, the spirochetes can persist in the mammalian host for months or years. The mechanisms by which Lyme disease spirochetes evade the immune response have not been determined. In this study, we have identified and characterized an elaborate genetic system in the Lyme disease spirochete B. burgdorferi that promotes extensive antigenic variation of a 34-kDa surface-exposed lipoprotein, VlsE. A 28-kilobase linear plasmid of B. burgdorferi B31 (lp28-1) was found to contain a vmp-like sequence (vls) locus that closely resembles the variable major protein (vmp) system for antigenic variation of relapsing fever organisms. The presence of lp28-1 correlates with the high-infectivity phenotype in B. burgdorferi strains tested. Segments of the 15 non-expressed (silent) vls cassette sequences located upstream of vlsE are able to recombine into the centra vlsE cassette region during infection of C3H/HeN mice, resulting in antigenic variation of the expressed lipoprotein. When compared to parental VlsE, VlsE variants progressively accumulate sequence changes during the period of 4, 7, 14, 21, and 28 days post infection in C3H/HeN mice. However, no recombination was detected during the period of 28-day in vitro culture, suggesting in vivo induction of VlsE antigenic variation. Adaptive immune responses do not appear to play a significant role in this induction, since similar recombination events were also observed in immunodeficient SCID mice. The $5\sp\prime$ and $3\sp\prime$ noncassette regions of vlsE are apparently not subject to recombination and sequence variation. The structure and sequence of the silent vls cassette locus is preserved during the process of the VlsE antigenic variation, consistent with a nonreciprocal recombination mechanism. This combinatorial form of antigenic variation could potentially yield millions of VlsE variants in the mammalian host, and thereby contribute to immune evasion, long-term survival, and pathogenesis of B. burgdorferi. ^

Fast parallel kriging-based stepwise uncertainty reduction with application to the identification of an excursion set

Stepwise uncertainty reduction (SUR) strategies aim at constructing a sequence of points for evaluating a function  f in such a way that the residual uncertainty about a quantity of interest progressively decreases to zero. Using such strategies in the framework of Gaussian process modeling has been shown to be efficient for estimating the volume of excursion of f above a fixed threshold. However, SUR strategies remain cumbersome to use in practice because of their high computational complexity, and the fact that they deliver a single point at each iteration. In this article we introduce several multipoint sampling criteria, allowing the selection of batches of points at which f can be evaluated in parallel. Such criteria are of particular interest when f is costly to evaluate and several CPUs are simultaneously available. We also manage to drastically reduce the computational cost of these strategies through the use of closed form formulas. We illustrate their performances in various numerical experiments, including a nuclear safety test case. Basic notions about kriging, auxiliary problems, complexity calculations, R code, and data are available online as supplementary materials.

MATHEMATICAL STUDIES ON THE EVOLUTIONARY CHANGE OF DNA SEQUENCES

With the aim of understanding the mechanism of molecular evolution, mathematical problems on the evolutionary change of DNA sequences are studied. The problems studied and the results obtained are as follows: (1) Estimation of evolutionary distance between nucleotide sequences. Studying the pattern of nucleotide substitution for the case of unequal substitution rates, a new mathematical formula for estimating the average number of nucleotide substitutions per site between two homologous DNA sequences is developed. It is shown that this formula has a wider applicability than currently available formulae. A statistical method for estimating the number of nucleotide changes due to deletion and insertion is also developed. (2) Biases of the estimates of nucleotide substitutions obtained by the restriction enzyme method. The deviation of the estimate of nucleotide substitutions obtained by the restriction enzyme method from the true value is investigated theoretically. It is shown that the amount of the deviation depends on the nucleotides in the recognition sequence of the restriction enzyme used, unequal rates of substitution among different nucleotides, and nucleotide frequences, but the primary factor is the unequal rates of nucleotide substitution. When many different kinds of enzymes are used, however, the amount of average deviation is generally small. (3) Distribution of restriction fragment lengths. To see the effect of undetectable restriction fragments and fragment differences on the estimate of nucleotide differences, the theoretical distribution of fragment lengths is studied. This distribution depends on the type of restriction enzymes used as well as on the relative frequencies of four nucleotides. It is shown that undetectability of small fragments or fragment differences gives a serious underestimate of nucleotide substitutions when the length-difference method of estimation is used, but the extent of underestimation is small when the site-difference method is used. (4) Evolutionary relationships of DNA sequences in finite populations. A mathematical theory on the expected evolutionary relationships among DNA sequences (nucleons) randomly chosen from the same or different populations is developed under the assumption that the evolutionary change of nucleons is determined solely by mutation and random genetic drift. . . . (Author's abstract exceeds stipulated maximum length. Discontinued here with permission of author). UMI ^

Analyses of a Holocene soil/loess sequence in Southern Siberia

We conducted a high-resolution study of a unique Holocene sequence of wind-blown sediments and buried soils in Southern Siberia, far from marine environment influences. This was accomplished in order to assess the difference between North Atlantic marine and in-land climate variations. Relative wind strength was determined by grain size analyses of different stratigraphic units. Petromagnetic measurements were performed to provide a proxy for the relative extent of pedogenesis. An age model for the sections was built using the radiocarbon dating method. The windy periods are associated with the absence of soil formation and relatively low values of frequency dependence of magnetic susceptibility (FD), which appeared to be a valuable quantitative marker of pedogenic activity. These events correspond to colder intervals which registered reduced solar modulation and sun spot number. Events, where wind strength was lower, are characterized by soil formation with high FD values. Spectral analysis of our results demonstrates periodic changes of 1500, 1000 and 500 years of relatively warm and cold intervals during the Holocene of Siberia. We presume that the 1000 and 500 year climatic cycles are driven by increased solar insolation reaching the Earth surface and amplified by other still controversial mechanisms. The 1500 year cycle associated with the North Atlantic circulation appears only in the Late Holocene. Three time periods - 8400-9300 years BP, 3600-5100 years BP, and the last ~250 years BP - correspond to both the highest sun spot number and the most developed soil horizons in the studied sections