Advances in Sequence Analysis: Theory, Methods, Applications

This book gives a general view of sequence analysis, the statistical study of successions of states or events. It includes innovative contributions on life course studies, transitions into and out of employment, contemporaneous and historical careers, and political trajectories. The approach presented in this book is now central to the life-course perspective and the study of social processes more generally. This volume promotes the dialogue between approaches to sequence analysis that developed separately, within traditions contrasted in space and disciplines. It includes the latest developments in sequential concepts, coding, atypical datasets and time patterns, optimal matching and alternative algorithms, survey optimization, and visualization. Field studies include original sequential material related to parenting in 19th-century Belgium, higher education and work in Finland and Italy, family formation before and after German reunification, French Jews persecuted in occupied France, long-term trends in electoral participation, and regime democratization. Overall the book reassesses the classical uses of sequences and it promotes new ways of collecting, formatting, representing and processing them. The introduction provides basic sequential concepts and tools, as well as a history of the method. Chapters are presented in a way that is both accessible to the beginner and informative to the expert.

Quantitative sequence analysis of FBN1 premature termination codons provides evidence for incomplete NMD in leukocytes.

We improved, evaluated, and used Sanger sequencing for quantification of single nucleotide polymorphism (SNP) variants in transcripts and gDNA samples. This improved assay resulted in highly reproducible relative allele frequencies (e.g., for a heterozygous gDNA 50.0+/-1.4%, and for a missense mutation-bearing transcript 46.9+/-3.7%) with a lower detection limit of 3-9%. It provided excellent accuracy and linear correlation between expected and observed relative allele frequencies. This sequencing assay, which can also be used for the quantification of copy number variations (CNVs), methylations, mosaicisms, and DNA pools, enabled us to analyze transcripts of the FBN1 gene in fibroblasts and blood samples of patients with suspected Marfan syndrome not only qualitatively but also quantitatively. We report a total of 18 novel and 19 known FBN1 sequence variants leading to a premature termination codon (PTC), 26 of which we analyzed by quantitative sequencing both at gDNA and cDNA levels. The relative amounts of PTC-containing FBN1 transcripts in fresh and PAXgene-stabilized blood samples were significantly higher (33.0+/-3.9% to 80.0+/-7.2%) than those detected in affected fibroblasts with inhibition of nonsense-mediated mRNA decay (NMD) (11.0+/-2.1% to 25.0+/-1.8%), whereas in fibroblasts without NMD inhibition no mutant alleles could be detected. These results provide evidence for incomplete NMD in leukocytes and have particular importance for RNA-based analyses not only in FBN1 but also in other genes.

Comprehensive analysis of Arabidopsis expression level polymorphisms with simple inheritance.

In Arabidopsis thaliana, gene expression level polymorphisms (ELPs) between natural accessions that exhibit simple, single locus inheritance are promising quantitative trait locus (QTL) candidates to explain phenotypic variability. It is assumed that such ELPs overwhelmingly represent regulatory element polymorphisms. However, comprehensive genome-wide analyses linking expression level, regulatory sequence and gene structure variation are missing, preventing definite verification of this assumption. Here, we analyzed ELPs observed between the Eil-0 and Lc-0 accessions. Compared with non-variable controls, 5' regulatory sequence variation in the corresponding genes is indeed increased. However, approximately 42% of all the ELP genes also carry major transcription unit deletions in one parent as revealed by genome tiling arrays, representing a >4-fold enrichment over controls. Within the subset of ELPs with simple inheritance, this proportion is even higher and deletions are generally more severe. Similar results were obtained from analyses of the Bay-0 and Sha accessions, using alternative technical approaches. Collectively, our results suggest that drastic structural changes are a major cause for ELPs with simple inheritance, corroborating experimentally observed indel preponderance in cloned Arabidopsis QTL.

Imatinib plasma concentrations variability in oncologic patients

Imatinib (Glivec®) has transformed the treatment and short-term prognosis of chronic myeloid leukemia (CML) and gastrointestinal stromal tumor (GIST). However, the treatment must be taken indefinitely, it is not devoid of inconvenience and toxicity. Moreover, resistance or escape from disease control occurs in a significant number of patients. Imatinib is a substrate of the cytochromes P450 CYP3A4/5 and of the multidrug transporter P-glycoprotein (product of the MDR1 gene). Considering the large inter-individual differences in the expression and function of those systems, the disposition and clinical activity of imatinib can be expected to vary widely among patients, calling for dosage individualization. The aim of this exploratory study was to determine the average pharmacokinetic parameters characterizing the disposition of imatinib in the target population, to assess their inter-individual variability, and to identify influential factors affecting them. A total of 321 plasma concentrations, taken at various sampling times after the latest dose, were measured in 59 patients receiving Glivec at diverse regimens, using a validated HPLC-UV method developed for this study. The results were analyzed by non-linear mixed effect modeling (NONMEM). A one-compartment model with first-order absorption appeared appropriate to describe the data, with an average apparent clearance of 12.4 l/h, a distribution volume of 268 l and an absorption constant of 0.47 h-1. The clearance was affected by body weight, age and sex. No influences of interacting drugs were found. DNA samples were used for pharmacogenetic explorations. At present, only the MDR1 polymorphism has been assessed and seems to affect the pharmacokinetic parameters of imatinib. Large inter-individual variability remained unexplained by the demographic covariates considered, both on clearance (40 %) and distribution volume (71 %). Together with intra-patient variability (34 %), this translates into an 8-fold width of the 90 %-prediction interval of plasma concentrations expected under a fixed dosing regimen. This is a strong argument to further investigate the possible usefulness of a therapeutic drug monitoring program for imatinib. It may help to individualize the dosing regimen before overt disease progression or observation of treatment toxicity, thus improving both the long-term therapeutic effectiveness and tolerability of this drug.

Interindividual variability of the clinical pharmacokinetics of methadone: implications for the treatment of opioid dependence

Methadone is widely used for the treatment of opioid dependence. Although in most countries the drug is administered as a racemic mixture of (R)- and (S)- methadone, (R)-methadone accounts for most, if not all, of the opioid effects. Methadone can be detected in the blood 15-45 minutes after oral administration, with peak plasma concentration at 2.5-4 hours. Methadone has a mean bioavailability of around 75% (range 36-100%). Methadone is highly bound to plasma proteins, in particular to alpha(1)-acid glycoprotein. Its mean free fraction is around 13%, with a 4-fold interindividual variation. Its volume of distribution is about 4 L/kg (range 2-13 L/kg). The elimination of methadone is mediated by biotransformation, followed by renal and faecal excretion. Total body clearance is about 0.095 L/min, with wide interindividual variation (range 0.02-2 L/min). Plasma concentrations of methadone decrease in a biexponential manner, with a mean value of around 22 hours (range 5-130 hours) for elimination half-life. For the active (R)-enantiomer, mean values of around 40 hours have been determined. Cytochrome P450 (CYP) 3A4 and to a lesser extent 2D6 are probably the main isoforms involved in methadone metabolism. Rifampicin (rifampin), phenobarbital, phenytoin, carbamazepine, nevirapine, and efavirenz decrease methadone blood concentrations, probably by induction of CYP3A4 activity, which can result in severe withdrawal symptoms. Inhibitors of CYP3A4, such as fluconazole, and of CYP2D6, such as paroxetine, increase methadone blood concentrations. There is an up to 17-fold interindividual variation of methadone blood concentration for a given dosage, and interindividual variability of CYP enzymes accounts for a large part of this variation. Since methadone probably also displays large interindividual variability in its pharmacodynamics, methadone treatment must be individually adapted to each patient. Because of the high morbidity and mortality associated with opioid dependence, it is of major importance that methadone is used at an effective dosage in maintenance treatment: at least 60 mg/day, but typically 80-100 mg/day. Recent studies also show that a subset of patients might benefit from methadone dosages larger than 100 mg/day, many of them because of high clearance. In clinical management, medical evaluation of objective signs and subjective symptoms is sufficient for dosage titration in most patients. However, therapeutic drug monitoring can be useful in particular situations. In the case of non-response trough plasma concentrations of 400 microg/L for (R,S)-methadone or 250 microg/L for (R)-methadone might be used as target values.

Cell disposition of raltegravir and newer antiretrovirals in HIV-infected patients: high inter-individual variability in raltegravir cellular penetration.

Objectives The site of pharmacological activity of raltegravir is intracellular. Our aim was to determine the extent of raltegravir cellular penetration and whether raltegravir total plasma concentration (C(tot)) predicts cellular concentration (C(cell)). Methods Open-label, prospective, pharmacokinetic study on HIV-infected patients on a stable raltegravir-containing regimen. Plasma and peripheral blood mononuclear cells were simultaneously collected during a 12 h dosing interval after drug intake. C(tot) and C(cell) of raltegravir, darunavir, etravirine, maraviroc and ritonavir were measured by liquid chromatography coupled to tandem mass spectrometry after protein precipitation. Longitudinal mixed effects analysis was applied to the C(cell)/C(tot) ratio. Results Ten HIV-infected patients were included. The geometric mean (GM) raltegravir total plasma maximum concentration (C(max)), minimum concentration (C(min)) and area under the time-concentration curve from 0-12 h (AUC(0-12)) were 1068 ng/mL, 51.1 ng/mL and 4171 ng·h/mL, respectively. GM raltegravir cellular C(max), C(min) and AUC(0-12) were 27.5 ng/mL, 2.9 ng/mL and 165 ng·h/mL, respectively. Raltegravir C(cell) corresponded to 5.3% of C(tot) measured simultaneously. Both concentrations fluctuate in parallel, with C(cell)/C(tot) ratios remaining fairly constant for each patient without a significant time-related trend over the dosing interval. The AUC(cell)/AUC(tot) GM ratios for raltegravir, darunavir and etravirine were 0.039, 0.14 and 1.55, respectively. Conclusions Raltegravir C(cell) correlated with C(tot) (r = 0.86). Raltegravir penetration into cells is low overall (∼5% of plasma levels), with distinct raltegravir cellular penetration varying by as much as 15-fold between patients. The importance of this finding in the context of development of resistance to integrase inhibitors needs to be further investigated.

The region of mouse mammary tumor virus DNA containing the long terminal repeat includes a long coding sequence and signals for hormonally regulated transcription.

Starting from a biologically active recombinant DNA clone of exogenous unintegrated GR mouse mammary tumor virus, we have generated three subclones of PstI fragments of 1.45, 1.1, and 2.0 kb in the plasmid vector PBR322. The nucleotide sequence has been determined for the clone of 1.45 kb which includes almost the complete region of the long terminal repeat (LTR) plus an adjacent stretch of unique sequence DNA. A short region of the 2.0 kb clone, containing the beginning of the LTR, has also been sequenced. Starting with the A of an initiation codon outside the LTR, we detected an open reading frame of 960 nucleotides, potentially coding for a protein of 320 amino acids (36K). Two hundred nucleotides downstream from the termination codon, and approximately 25 nucleotides upstream from the presumptive initiation site of viral RNA synthesis, we found a promoter-like sequence. The sequence AGTAAA was detected approximately 15-20 nucleotides upstream from the 3' end of virion RNA and probably serves as a polyadenylation signal. The 1.45 kb PstI fragment has been transfected into Ltk- cells together with a plasmid containing the thymidine kinase gene of herpes simplex virus. The virus-specific RNA synthesis detected in a Tk+ cell clone was strongly stimulated by the addition of dexamethasone.

Usefulness of double locus sequence typing (DLST) for regional and international epidemiological surveillance of methicilin-resistant Staphylococcus aureus.

Methicillin-resistant Staphylococcus aureus (MRSA) is a major cause of nosocomial infections worldwide. To differentiate reliably among S. aureus isolates, we recently developed double locus sequence typing (DLST) based on the analysis of partial sequences of clfB and spa genes. In the present study, we evaluated the usefulness of DLST for epidemiological investigations of MRSA by routinely typing 1242 strains isolated in Western Switzerland. Additionally, particular local and international collections were typed by pulsed field gel electrophoresis (PFGE) and DLST to check the compatibility of DLST with the results obtained by PFGE, and for international comparisons. Using DLST, we identified the major MRSA clones of Western Switzerland, and demonstrated the close relationship between local and international clones. The congruence of 88% between the major PFGE and DLST clones indicated that our results obtained by DLST were compatible with earlier results obtained by PFGE. DLST could thus easily be incorporated in a routine surveillance procedure. In addition, the unambiguous definition of DLST types makes this method more suitable than PFGE for long-term epidemiological surveillance. Finally, the comparison of the results obtained by DLST, multilocus sequence typing, PFGE, Staphylococcal cassette chromosome mec typing and the detection of Panton-Valentine leukocidin genes indicated that no typing scheme should be used on its own. It is only the combination of data from different methods that gives the best chance of describing precisely the epidemiology and phylogeny of MRSA.

Imatinib plasma concentrations variability in oncologic patients

The anticancer drug imatinib has transformed the treatment and prognosis of chronic myeloid leukemia and gastrointestinal stromal tumor. However, the treatment must be taken indefinitely and is not devoid of inconveniences and toxicity. Moreover, resistance or escape from disease control are occurring. Considering the large interindividual differences in the function of the enzymatic and transport systems involved in imatinib disposition, exposure to this drug can be expected to vary widely among patients. This book describes an observational clinical trial aiming at exploring the influence of these covariates on imatinib pharmacokinetics and assessing the interindividual variability of the pharmacokinetic parameters of the drug. A large interindividual variability was observed, together with some preliminary concentration-effect relationships. These elements are arguments to further investigate the potential benefit of a therapeutic drug monitoring program to optimize the use of imatinib in patients. Such results should be especially useful to clinical oncologists or scientists involved in clinical oncology research.

Three-dimensional sequence stratigraphy and subtle stratigraphic traps associated with systems tracts: West Cameron region, offshore Louisiana, Gulf of Mexico

Three-dimensional sequence stratigraphy is a potent exploration and development tool for the discovery of subtle stratigraphic traps. Reservoir morphology, heterogeneity and subtle stratigraphic trapping mechanisms can be better understood through systematic horizontal identification of sedimentary facies of systems tracts provided by three-dimensional attribute maps used as an important complement to the sequential analysis on the two-dimensional seismic lines and the well log data. On new prospects as well as on already-producing fields, the additional input of sequential analysis on three-dimensional data enables the identification, location and precise delimitation of new potentially productive zones. The first part of this paper presents four typical horizontal seismic facies assigned to the successive systems tracts of a third- or fourth-order sequence deposited in inner to outer neritic conditions on a elastic shelf. The construction of this synthetic representative sequence is based on the observed reproducibility of the horizontal seismic facies response to cyclic eustatic events on more than 35 sequences registered in the Gulf coast Plio-Pleistocene and Late Miocene, offshore Louisiana in the West Cameron region of the Gulf of Mexico. The second part shows how three-dimensional sequence stratigraphy can contribute in localizing and understanding sedimentary facies associated with productive zones. A case study in the early Middle Miocene Cibicides opima sands shows multiple stacked gas accumulations in the top slope fan, prograding wedge and basal transgressive systems tract of the third-order sequence between SB15.5 and SB 13.8 Ma.

Variability of the relative afferent pupillary defect.

PURPOSE: Afferent asymmetry of visual function is detectable in both normal and pathologic conditions. With a computerized test, we assessed the variability in measuring afferent asymmetry of the pupillary light reflex, that is, the relative afferent pupillary defect. METHODS: In ten normal subjects, pupillary responses to an alternating light stimulus were recorded with computerized infrared pupillography. The relative afferent pupillary defect for each test was determined by using a new computer analysis. The 95% confidence interval of each determination of relative afferent pupillary defect was used to represent the short-term fluctuation in its measurement. To optimize the test for clinical use, we studied the influence of stimulus intensity, duration, and number on the variability of the relative afferent pupillary defect. RESULTS: When the relative afferent pupillary defect was based on only a few light alternations (stimulus pairs), there was excessive variability in its measurement (95% confidence interval > 0.5 log units). With approximately 200 stimulus pairs, the 95% confidence interval was reduced to less than 0.1 log unit (relative afferent pupillary defect +/- 0.05 log unit). Also, there was less variability when the dark interval between alternating light stimulation was less than one second. CONCLUSIONS: Computerized infrared pupillography can standardize the alternating light test and minimize the error in quantifying a relative afferent pupillary defect. A reproducible relative afferent pupillary defect measurement is desirable for defining afferent injury and following the course of disease.

Effect of the September 11, 2001 terrorist attack on a state highway patrol trooper's heart rate variability

BACKGROUND: On September 11, 2001, terrorists attacked the United States. By coincidence, a North Carolina highway patrol trooper was wearing an ambulatory ECG Holter monitor at this time as part of an air pollution study. METHODS: Heart rate variability parameters were analyzed: standard deviation of normal to normal beat intervals (SDNN) and percentage of interval differences >50 ms (PNN50). RESULTS: The trooper's heart rate variability changed immediately after learning about the terrorist attacks. Heart rate increased and PNN50 decreased, while SDNN increased strongly. CONCLUSIONS: These changes suggest strong emotional sympathetic stress associated with parasympathetic withdrawal in response to the news about the terrorist attack. [Authors]

Respiratory variability in mechanically ventilated patients. Critical Care 2011

Introduction: Increased respiratory pattern variability is associated with improved oxygenation. Pressure support (PS) is a widely used partial-assist mechanical ventilation (MV) mode, in which each breathing cycle is initiated by flow or pressure variation at the airway due to patient inspiratory effort. Neurally adjusted ventilatory assist (NAVA) is relatively new and uses the electrical activity of the diaphragm (Eadi) to deliver ventilatory support proportional to the patient's inspiratory demand. We hypothesize that respiratory variability should be greater with NAVA compared with PS.Methods: Twenty-two patients underwent 20 minutes of PS followed by 20 minutes of NAVA. Flow and Eadi curves were used to obtain tidal volume (Vt) and ∫Eadi for 300 to 400 breaths in each patient. Patient-specific cumulative distribution functions (CDF) show the percentage Vt and ∫Eadi within a clinically defined (±10%) variability band for each patient. Values are normalized to patient-specific medians for direct comparison. Variability in Vt (outcome) is thus expressed in terms of variability in ∫Eadi (demand) on the same plot.Results: Variability in Vt relative to variability in ∫Eadi is significantly greater for NAVA than PS (P = 0.00012). Hence, greater variability in outcome Vt is obtained for a given demand in ∫Eadi, under NAVA, as illustrated in Figure 1 for a typical patient. A Fisher 2 × 2 contingency analysis showed that 45% of patients under NAVA had a Vt variability in equal proportion to ∫Eadi variability, versus 0% for PS (P < 0.05).Conclusions: NAVA yields greater variability in tidal volume, relative to ∫Eadi demand, and a better match between Vt and ∫Eadi. These results indicate that NAVA could achieve improved oxygenation compared with PS when sufficient underlying variability in ∫Eadi is present, due to its ability to achieve higher tidal volume variability from a given variability in ∫Eadi.