On the inference of complex phylogenetic networks by Markov Chain Monte-Carlo

For various species, high quality sequences and complete genomes are nowadays available for many individuals. This makes data analysis challenging, as methods need not only to be accurate, but also time efficient given the tremendous amount of data to process. In this article, we introduce an efficient method to infer the evolutionary history of individuals under the multispecies coalescent model in networks (MSNC). Phylogenetic networks are an extension of phylogenetic trees that can contain reticulate nodes, which allow to model complex biological events such as horizontal gene transfer, hybridization and introgression. We present a novel way to compute the likelihood of biallelic markers sampled along genomes whose evolution involved such events. This likelihood computation is at the heart of a Bayesian network inference method called SnappNet, as it extends the Snapp method inferring evolutionary trees under the multispecies coalescent model, to networks. SnappNet is available as a package of the well-known beast 2 software.Recently, the MCMC_BiMarkers method, implemented in PhyloNet, also extended Snapp to networks. Both methods take biallelic markers as input, rely on the same model of evolution and sample networks in a Bayesian framework, though using different methods for computing priors. However, SnappNet relies on algorithms that are exponentially more time-efficient on non-trivial networks. Using simulations, we compare performances of SnappNet and MCMC_BiMarkers. We show that both methods enjoy similar abilities to recover simple networks, but SnappNet is more accurate than MCMC_BiMarkers on more complex network scenarios. Also, on complex networks, SnappNet is found to be extremely faster than MCMC_BiMarkers in terms of time required for the likelihood computation. We finally illustrate SnappNet performances on a rice data set. SnappNet infers a scenario that is consistent with previous results and provides additional understanding of rice evolution.     

Computing folding pathways between RNA secondary structures

Given an RNA sequence and two designated secondary structures A, B, we describe a new algorithm that computes a nearly optimal folding pathway from A to B. The algorithm, RNAtabupath, employs a tabu semi-greedy heuristic, known to be an effective search strategy in combinatorial optimization. Folding pathways, sometimes called routes or trajectories, are computed by RNAtabupath in a fraction of the time required by the barriers program of Vienna RNA Package. We benchmark RNAtabupath with other algorithms to compute low energy folding pathways between experimentally known structures of several conformational switches. The RNApathfinder web server, source code for algorithms to compute and analyze pathways and supplementary data are available at http://bioinformatics.bc.edu/clotelab/RNApathfinder.     

Allele-specific copy number profiling by next-generation DNA sequencing

The progression and clonal development of tumors often involve amplifications and deletions of genomic DNA. Estimation of allele-specific copy number, which quantifies the number of copies of each allele at each variant loci rather than the total number of chromosome copies, is an important step in the characterization of tumor genomes and the inference of their clonal history. We describe a new method, falcon, for finding somatic allele-specific copy number changes by next generation sequencing of tumors with matched normals. falcon is based on a change-point model on a bivariate mixed Binomial process, which explicitly models the copy numbers of the two chromosome haplotypes and corrects for local allele-specific coverage biases. By using the Binomial distribution rather than a normal approximation, falcon more effectively pools evidence from sites with low coverage. A modified Bayesian information criterion is used to guide model selection for determining the number of copy number events. Falcon is evaluated on in silico spike-in data and applied to the analysis of a pre-malignant colon tumor sample and late-stage colorectal adenocarcinoma from the same individual. The allele-specific copy number estimates obtained by falcon allows us to draw detailed conclusions regarding the clonal history of the individual''s colon cancer.     

Mechanism of primary and secondary ion-radical pair formation in photosystem I complexes

The mechanisms of the ultrafast charge separation in reaction centers of photosystem I (PS I) complexes are discussed. A kinetic model of the primary reactions in PS I complexes is presented. The model takes into account previously calculated values of redox potentials of cofactors, reorganization energies of the primary P700⁺A ₀ ^- and secondary P700⁺A ₁ ^- ion-radical pairs formation, and the possibility of electron transfer via both symmetric branches A and B of redox-cofactors. The model assumes that the primary electron acceptor A₀ in PS I is represented by a dimer of chlorophyll molecules Chl2A/Chl3A and Chl2B/Chl3B in branches A and B of the cofactors. The characteristic times of formation of P700⁺A ₀ ^- and P700⁺A ₁ ^- calculated on the basis of the model are close to the experimental values obtained by pump-probe femtosecond absorption spectroscopy. It is demonstrated that a small difference in the values of redox potentials between the primary electron acceptors A_0A and A_0B in branches A and B leads to asymmetry of the electron transfer in a ratio of 70: 30 in favor of branch A. The secondary charge separation is thermodynamically irreversible in the submicrosecond range and is accompanied by additional increase in asymmetry between the branches of cofactors of PS I.     

Charge separation in Rhodobacter sphaeroides mutant reaction centers with increased midpoint potential of the primary electron donor

Primary charge separation dynamics in four mutant reaction centers (RCs) of the purple bacterium Rhodobacter sphaeroides with increased midpoint potential of the primary electron donor P (M160LH, L131LH, M197FH, and M160LH + L131LH + M197FH) have been studied by femtosecond transient absorption spectroscopy at room temperature. The decay of the excited singlet state in the wild-type and mutant RCs is complex and has two main exponential components, which indicates heterogeneity of electron transfer rates or the presence of reverse electron transfer reactions. The radical anion band of monomeric bacteriochlorophyll B_A at 1020 nm was first observed in transient absorbance difference spectra of single mutants. This band remains visible, although with somewhat reduced amplitude, even at delays up to tens of picoseconds when stimulated emission is absent and the reaction centers are in the P⁺H _A ^? state. The presence of this band in this time period indicates the existence of thermodynamic equilibrium between the P⁺B _A ^? H_A and P⁺B_AH _A ^? states. The data give grounds for assuming that the value of the energy difference between the states P*, P⁺B _A ^? H_A, and P⁺B_AH _A ^? at early times is of the same order of magnitude as the energy kT at room temperature. Besides, monomeric bacteriochlorophyll B_A is found to be an immediate electron acceptor in the single mutant RCs, where electron transfer is hampered due to increased energy of the P⁺B _A ^? state with respect to P*.     

Molecular and Biochemical Characterization of the Protein Template Controlling Biosynthesis of the Lipopeptide Lichenysin

Lichenysins are surface-active lipopeptides with antibiotic properties produced nonribosomally by several strains of Bacillus licheniformis. Here, we report the cloning and sequencing of an entire 26.6-kb lichenysin biosynthesis operon from B. licheniformis ATCC 10716. Three large open reading frames coding for peptide synthetases, designated licA, licB (three modules each), and licC (one module), could be detected, followed by a gene, licTE, coding for a thioesterase-like protein. The domain structure of the seven identified modules, which resembles that of the surfactin synthetases SrfA-A to -C, showed two epimerization domains attached to the third and sixth modules. The substrate specificity of the first, fifth, and seventh recombinant adenylation domains of LicA to -C (cloned and expressed in Escherichia coli) was determined to be Gln, Asp, and Ile (with minor Val and Leu substitutions), respectively. Therefore, we suppose that the identified biosynthesis operon is responsible for the production of a lichenysin variant with the primary amino acid sequence l-Gln–l-Leu–d-Leu–l-Val–l-Asp–d-Leu–l-Ile, with minor Leu and Val substitutions at the seventh position.Many strains of Bacillus are known to produce lipopeptides with remarkable surface-active properties (11). The most prominent of these powerful lipopeptides is surfactin from Bacillus subtilis (1). Surfactin is an acylated cyclic heptapeptide that reduces the surface tension of water from 72 to 27 mN m⁻¹ even in a concentration below 0.05% and shows some antibacterial and antifungal activities (1). Some B. subtilis strains are also known to produce other, structurally related lipoheptapeptides (Table ​(Table1),1), like iturin (32, 34) and bacillomycin (3, 27, 30), or the lipodecapeptides fengycin (50) and plipastatin (29).

TABLE 1

Lipoheptapeptide antibiotics of Bacillus spp.

Reverse engineering and analysis of large genome-scale gene networks

Reverse engineering the whole-genome networks of complex multicellular organisms continues to remain a challenge. While simpler models easily scale to large number of genes and gene expression datasets, more accurate models are compute intensive limiting their scale of applicability. To enable fast and accurate reconstruction of large networks, we developed Tool for Inferring Network of Genes (TINGe), a parallel mutual information (MI)-based program. The novel features of our approach include: (i) B-spline-based formulation for linear-time computation of MI, (ii) a novel algorithm for direct permutation testing and (iii) development of parallel algorithms to reduce run-time and facilitate construction of large networks. We assess the quality of our method by comparison with ARACNe (Algorithm for the Reconstruction of Accurate Cellular Networks) and GeneNet and demonstrate its unique capability by reverse engineering the whole-genome network of Arabidopsis thaliana from 3137 Affymetrix ATH1 GeneChips in just 9 min on a 1024-core cluster. We further report on the development of a new software Gene Network Analyzer (GeNA) for extracting context-specific subnetworks from a given set of seed genes. Using TINGe and GeNA, we performed analysis of 241 Arabidopsis AraCyc 8.0 pathways, and the results are made available through the web.     

Chromosome morphology of the czechoslovak species of the genusScorzonera

Karyotypic formulae of theScorzonera L. species are as follows:S. purpurea L.: K (2n)=14=8 A^m+6 B^sm and K (2n)=14+1=9 A^m+6 B^sm;S. austriaca Willd.: K (2n)=14=6 A^m+6 B^sm+2 C^st;S. humilis L.: K (2n)=14=12 A^m+2 B^sm;S. parviflora Jacq.: K (2n)=14=10 A^m+4 B^sm;S. hispanica L.: K (2n)=12 A^m+2 B^sm. The results of the study of the karyotypes of the generaScorzonera L.,Tragopogon L. andPodospermum DC. are summarized.     

Accurate and efficient reconstruction of deep phylogenies from structured RNAs

Ribosomal RNA (rRNA) genes are probably the most frequently used data source in phylogenetic reconstruction. Individual columns of rRNA alignments are not independent as a consequence of their highly conserved secondary structures. Unless explicitly taken into account, these correlation can distort the phylogenetic signal and/or lead to gross overestimates of tree stability. Maximum likelihood and Bayesian approaches are of course amenable to using RNA-specific substitution models that treat conserved base pairs appropriately, but require accurate secondary structure models as input. So far, however, no accurate and easy-to-use tool has been available for computing structure-aware alignments and consensus structures that can deal with the large rRNAs. The RNAsalsa approach is designed to fill this gap. Capitalizing on the improved accuracy of pairwise consensus structures and informed by a priori knowledge of group-specific structural constraints, the tool provides both alignments and consensus structures that are of sufficient accuracy for routine phylogenetic analysis based on RNA-specific substitution models. The power of the approach is demonstrated using two rRNA data sets: a mitochondrial rRNA set of 26 Mammalia, and a collection of 28S nuclear rRNAs representative of the five major echinoderm groups.     

Studies of the morphology of chromosomes of some selected species

Karyotypes of twelve species from twenty-four localities in southern Moravia and one locality in Slovakia were investigated. Their counts or karyotypic formulae are as follows:Chenopodium foliosum (Moench) Ascherson: K (2n)=18=16 A^m+2 B^sm;Astragalus austriacus Jacq.: K (2n)=16=8 A^m+8 B^sm;Astragalus exscapus L.: K (2n)=16=10 A^m+4 B^sm+2 C^st;Astragalus cicer L.: K (2n)=64;Astragalus onobrychis L.: K (2n=64 and K (2n)=64+1;Vicia dumetorum L.: K (2n=14=10 A^m+4 B^sm;Vicia sylvatica L.: K (2n)=14=2 A^m+10 B^sm+2 C^st;Vicia pisiformis L.: K (2n)=12=8 A^m+4 B^sm;Vicia cassubica L.: K (2n)=12=4 A^m+6 B^sm+2 C^st;Vicia cracca L. (from five localities in southern Moravia): K (2n)=28=4 A^m+12 B^sm+12 C^st and K (2n)=28+1=5 A^m+12 B^sm+12 C^st;Vicia cracca L. (from one locality in Slovakia): K (2n)=14=2 A^m+6 B^sm+6 C^st;Vicia tenuifolia Roth: K (2n)=24=4 A^m+16 B^sm+4 C^st;Serratula lycopifolia (Vill.) Kern.: K (2n)=60.     

Importing statistical measures into Artemis enhances gene identification in the <Emphasis Type="Italic">Leishmania</Emphasis> genome project

Background

Seattle Biomedical Research Institute (SBRI) as part of the Leishmania Genome Network (LGN) is sequencing chromosomes of the trypanosomatid protozoan species Leishmania major. At SBRI, chromosomal sequence is annotated using a combination of trained and untrained non-consensus gene-prediction algorithms with ARTEMIS, an annotation platform with rich and user-friendly interfaces.

Results

Here we describe a methodology used to import results from three different protein-coding gene-prediction algorithms (GLIMMER, TESTCODE and GENESCAN) into the ARTEMIS sequence viewer and annotation tool. Comparison of these methods, along with the CODON USAGE algorithm built into ARTEMIS, shows the importance of combining methods to more accurately annotate the L. major genomic sequence.

Conclusion

An improvised and powerful tool for gene prediction has been developed by importing data from widely-used algorithms into an existing annotation platform. This approach is especially fruitful in the Leishmania genome project where there is large proportion of novel genes requiring manual annotation.

     

Essentiality of c-di-AMP in Bacillus subtilis: Bypassing mutations converge in potassium and glutamate homeostasis

In order to adjust to changing environmental conditions, bacteria use nucleotide second messengers to transduce external signals and translate them into a specific cellular response. Cyclic di-adenosine monophosphate (c-di-AMP) is the only known essential nucleotide second messenger. In addition to the well-established role of this second messenger in the control of potassium homeostasis, we observed that glutamate is as toxic as potassium for a c-di-AMP-free strain of the Gram-positive model bacterium Bacillus subtilis. In this work, we isolated suppressor mutants that allow growth of a c-di-AMP-free strain under these toxic conditions. Characterization of glutamate resistant suppressors revealed that they contain pairs of mutations, in most cases affecting glutamate and potassium homeostasis. Among these mutations, several independent mutations affected a novel glutamate transporter, AimA (Amino acid importer A, formerly YbeC). This protein is the major transporter for glutamate and serine in B. subtilis. Unexpectedly, some of the isolated suppressor mutants could suppress glutamate toxicity by a combination of mutations that affect phospholipid biosynthesis and a specific gain-of-function mutation of a mechanosensitive channel of small conductance (YfkC) resulting in the acquisition of a device for glutamate export. Cultivation of the c-di-AMP-free strain on complex medium was an even greater challenge because the amounts of potassium, glutamate, and other osmolytes are substantially higher than in minimal medium. Suppressor mutants viable on complex medium could only be isolated under anaerobic conditions if one of the two c-di-AMP receptor proteins, DarA or DarB, was absent. Also on complex medium, potassium and osmolyte toxicity are the major bottlenecks for the growth of B. subtilis in the absence of c-di-AMP. Our results indicate that the essentiality of c-di-AMP in B. subtilis is caused by the global impact of the second messenger nucleotide on different aspects of cellular physiology.     

Engineering of Glarea lozoyensis for Exclusive Production of the Pneumocandin B0 Precursor of the Antifungal Drug Caspofungin Acetate

Pneumocandins produced by the fungus Glarea lozoyensis are acylated cyclic hexapeptides of the echinocandin family. Pneumocandin B₀ is the starting molecule for the first semisynthetic echinocandin antifungal drug, caspofungin acetate. In the wild-type strain, pneumocandin B₀ is a minor fermentation product, and its industrial production was achieved by a combination of extensive mutation and medium optimization. The pneumocandin biosynthetic gene cluster was previously elucidated by a whole-genome sequencing approach. Knowledge of the biosynthetic cluster suggested an alternative way to produce exclusively pneumocandin B₀. Disruption of GLOXY4, encoding a nonheme, α-ketoglutarate-dependent oxygenase, confirmed its involvement in l-leucine cyclization to form 4S-methyl-l-proline. The absence of 4S-methyl-l-proline abolishes pneumocandin A₀ production, and 3S-hydroxyl-l-proline occupies the hexapeptide core''s position 6, resulting in exclusive production of pneumocandin B₀. Retrospective analysis of the GLOXY4 gene in a previously isolated pneumocandin B₀-exclusive mutant (ATCC 74030) indicated that chemical mutagenesis disrupted the GLOXY4 gene function by introducing two amino acid mutations in GLOXY4. This one-step genetic manipulation can rationally engineer a high-yield production strain.     

Characterization of the Complete Mitochondrial Genome of Cerura menciana and Comparison with Other Lepidopteran Insects

The complete mitochondrial genome (mitogenome) of Cerura menciana (Lepidoptera: Notodontidae) was sequenced and analyzed in this study. The mitogenome is a circular molecule of 15,369 bp, containing 13 protein-coding genes (PCGs), two ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes and a A+T-rich region. The positive AT skew (0.031) indicated that more As than Ts were present. All PCGs were initiated by ATN codons, except for the cytochrome c oxidase subunit 1 (cox1) gene, which was initiated by CAG. Two of the 13 PCGs contained the incomplete termination codon T or TA, while the others were terminated with the stop codon TAA. The A+T-rich region was 372 bp in length and consisted of an ‘ATAGA’ motif followed by an 18 bp poly-T stretch, a microsatellite-like (AT)₈ and a poly-A element upstream of the trnM gene. Results examining codon usage indicated that Asn, Ile, Leu2, Lys, Tyr and Phe were the six most frequently occurring amino acids, while Cys was the rarest. Phylogenetic relationships, analyzed based on the nucleotide sequences of the 13 PCGs from other insect mitogenomes, confirmed that C. menciana belongs to the Notodontidae family.     

A Collection of Popcorn as a Reservoir of Genes for the Generation of Lineages

In the present study, we analyze the genetic structure and diversity among accessions of popcorn obtained from the CIMMYT International Maize and Wheat Improvement Center that represent the diversity available for current use by breeding programs. The main objectives were to identify SSR loci that were the best indicators of genetic diversity, to measure the genetic diversity within popcorn genotypes, and to analyze the genetic structure of the promising populations destined for use in breeding programs. The mean gene diversity of all SSR loci was 0.6352. An extremely high population differentiation level was detected (F _st  = 0.3152) with F _st for each locus ranging from 0.1125 (Umc1229) to 0.4870 (Umc1755). Analyzing the genetic structure of eight popcorn accessions was especially important for identifying both SSR loci with high levels of heterozygosity and genotypes showing high heterozygosity (Boya462 and Arzm13 050). This analysis should be the medium and long-term selection goal for the generation of inbred lines and the future production of new cultivars. Plant accessions Arzm05 083, Arzm13 050, and Urug298 may also be useful varieties that exhibit important agronomic characteristics that can be used through crosses to broaden the genetic basis of popcorn.     

Modeling Effect of a ��-Secretase Inhibitor on Amyloid-�� Dynamics Reveals Significant Role of an Amyloid Clearance Mechanism

Aggregation of the small peptide amyloid beta (A??) into oligomers and fibrils in the brain is believed to be a precursor to Alzheimer??s disease. A?? is produced via multiple proteolytic cleavages of amyloid precursor protein (APP), mediated by the enzymes ??- and ??-secretase. In this study, we examine the temporal dynamics of soluble (unaggregated) A?? in the plasma and cerebral-spinal fluid (CSF) of rhesus monkeys treated with different oral doses of a ??-secretase inhibitor. A dose-dependent reduction of A?? concentration was observed within hours of drug ingestion, for all doses tested. A?? concentration in the CSF returned to its predrug level over the monitoring period. In contrast, A?? concentration in the plasma exhibited an unexpected overshoot to as high as 200% of the predrug concentration, and this overshoot persisted as late as 72 hours post-drug ingestion. To account for these observations, we proposed and analyzed a minimal physiological model for A?? dynamics that could fit the data. Our analysis suggests that the overshoot arises from the attenuation of an A?? clearance mechanism, possibly due to the inhibitor. Our model predicts that the efficacy of A?? clearance recovers to its basal (pretreatment) value with a characteristic time of >48 hours, matching the time-scale of the overshoot. These results point to the need for a more detailed investigation of soluble A?? clearance mechanisms and their interaction with A??-reducing drugs.     

Colinearity and Similar Expression Pattern of Rice DREB1s Reveal Their Functional Conservation in the Cold-Responsive Pathway

The clustered genes C-repeat (CRT) binding factor (CBF)1/ dehydration-responsive element binding protein (DREB)1B, CBF2/DREB1C, and CBF3/DREB1A play a central role in cold acclimation and facilitate plant resistance to freezing in Arabidopsis thaliana. Rice (Oryza sativa L.) is very sensitive to low temperatures; enhancing the cold stress tolerance of rice is a key challenge to increasing its yield. In this study, we demonstrate chilling acclimation, a phenomenon similar to Arabidopsis cold acclimation, in rice. To determine whether rice CBF/DREB1 genes participate in this cold-responsive pathway, all putative homologs of Arabidopsis DREB1 genes were filtered from the complete rice genome through a BLASTP search, followed by phylogenetic, colinearity and expression analysis. We thereby identified 10 rice genes as putative DREB1 homologs: nine of these were located in rice genomic regions with some colinearity to the Arabidopsis CBF1–CBF4 region. Expression profiling revealed that six of these genes (Os01g73770, Os02g45450, Os04g48350, Os06g03670, Os09g35010, and Os09g35030) were similarly expressed in response to chilling acclimation and cold stress and were co-expressed with genes involved in cold signalling, suggesting that these DREB1 homologs may be involved in the cold response in rice. The results presented here serve as a prelude towards understanding the function of rice homologs of DREB1 genes in cold-sensitive crops.     

A Variable Sampling Interval Synthetic Xbar Chart for the Process Mean

The usual practice of using a control chart to monitor a process is to take samples from the process with fixed sampling interval (FSI). In this paper, a synthetic X¯ control chart with the variable sampling interval (VSI) feature is proposed for monitoring changes in the process mean. The VSI synthetic X¯ chart integrates the VSI X¯ chart and the VSI conforming run length (CRL) chart. The proposed VSI synthetic X¯ chart is evaluated using the average time to signal (ATS) criterion. The optimal charting parameters of the proposed chart are obtained by minimizing the out-of-control ATS for a desired shift. Comparisons between the VSI synthetic X¯ chart and the existing X¯, synthetic X¯, VSI X¯ and EWMA X¯ charts, in terms of ATS, are made. The ATS results show that the VSI synthetic X¯ chart outperforms the other X¯ type charts for detecting moderate and large shifts. An illustrative example is also presented to explain the application of the VSI synthetic X¯ chart.     

Hepatitis C Virus (HCV) Sequence Variation Induces an HCV-Specific T-Cell Phenotype Analogous to Spontaneous Resolution

Hepatitis C virus (HCV)-specific CD8⁺ T cells in persistent HCV infection are low in frequency and paradoxically show a phenotype associated with controlled infections, expressing the memory marker CD127. We addressed to what extent this phenotype is dependent on the presence of cognate antigen. We analyzed virus-specific responses in acute and chronic HCV infections and sequenced autologous virus. We show that CD127 expression is associated with decreased antigenic stimulation after either viral clearance or viral variation. Our data indicate that most CD8 T-cell responses in chronic HCV infection do not target the circulating virus and that the appearance of HCV-specific CD127⁺ T cells is driven by viral variation.Hepatitis C virus (HCV) persists in the majority of acutely infected individuals, potentially leading to chronic hepatitis, liver cirrhosis, and hepatocellular carcinoma. The cellular immune response has been shown to play a significant role in viral control and protection from liver disease. Phenotypic and functional studies of virus-specific T cells have attempted to define the determinants of a successful versus an unsuccessful T-cell response in viral infections (10). So far these studies have failed to identify consistent distinguishing features between a T-cell response that results in self-limiting versus chronic HCV infection; similarly, the impact of viral persistence on HCV-specific memory T-cell formation is poorly understood.Interleukin-7 (IL-7) receptor alpha chain (CD127) is a key molecule associated with the maintenance of memory T-cell populations. Expression of CD127 on CD8 T cells is typically only observed when the respective antigen is controlled and in the presence of significant CD4⁺ T-cell help (9). Accordingly, cells specific for persistent viruses (e.g., HIV, cytomegalovirus [CMV], and Epstein-Barr virus [EBV]) have been shown to express low levels of CD127 (6, 12, 14) and to be dependent on antigen restimulation for their maintenance. In contrast, T cells specific for acute resolving virus infections, such as influenza virus, respiratory syncytial virus (RSV), hepatitis B virus (HBV), and vaccinia virus typically acquire expression of CD127 rapidly with the control of viremia (5, 12, 14). Results for HCV have been inconclusive. The expected increase in CD127 levels in acute resolving but not acute persisting infection has been found, while a substantial proportion of cells with high CD127 expression have been observed in long-established chronic infection (2). We tried to reconcile these observations by studying both subjects with acute and chronic HCV infection and identified the presence of antigen as the determinant of CD127 expression.Using HLA-peptide multimers we analyzed CD8⁺ HCV-specific T-cell responses and CD127 expression levels in acute and chronic HCV infection. We assessed a cohort of 18 chronically infected subjects as well as 9 individuals with previously resolved infection. In addition, we longitudinally studied 9 acutely infected subjects (5 individuals who resolved infection spontaneously and 4 individuals who remain chronically infected) (Tables ​(Tables11 and ​and2).2). Informed consent in writing was obtained from each patient, and the study protocol conformed to the ethical guidelines of the 1975 Declaration of Helsinki, as reflected in a priori approval from the local institutional review boards. HLA-multimeric complexes were obtained commercially from Proimmune (Oxford, United Kingdom) and Beckman Coulter (CA). The staining and analysis procedure was as described previously (10). Peripheral blood mononuclear cells (PBMCs) were stained with the following antibodies: CD3 from Caltag; CD8, CD27, CCR7, CD127, and CD38 from BD Pharmingen; and PD-1 (kindly provided by Gordon Freeman). Primer sets were designed for different genotypes based on alignments of all available sequences from the public HCV database (http://hcvpub.ibcp.fr). Sequence analysis was performed as previously described (8).

TABLE 1.

Patient information and autologous sequence analysis for patients with chronic and resolved HCV infection

DIRAQ: scalable in situ data- and resource-aware indexing for optimized query performance

Scientific data analytics in high-performance computing environments has been evolving along with the advancement of computing capabilities. With the onset of exascale computing, the increasing gap between compute performance and I/O bandwidth has rendered the traditional post-simulation processing a tedious process. Despite the challenges due to increased data production, there exists an opportunity to benefit from “cheap” computing power to perform query-driven exploration and visualization during simulation time. To accelerate such analyses, applications traditionally augment, post-simulation, raw data with large indexes, which are then repeatedly utilized for data exploration. However, the generation of current state-of-the-art indexes involves a compute- and memory-intensive processing, thus rendering them inapplicable in an in situ context. In this paper we propose DIRAQ, a parallel in situ, in network data encoding and reorganization technique that enables the transformation of simulation output into a query-efficient form, with negligible runtime overhead to the simulation run. DIRAQ’s effective core-local, precision-based encoding approach incorporates an embedded compressed index that is 3–6 \(\times \) smaller than current state-of-the-art indexing schemes. Its data-aware index adjustmentation improves performance of group-level index layout creation by up to 35 % and reduces the size of the generated index by up to 27 %. Moreover, DIRAQ’s in network index merging strategy enables the creation of aggregated indexes that speed up spatial-context query responses by up to \(10\times \) versus alternative techniques. DIRAQ’s topology-, data-, and memory-aware aggregation strategy results in efficient I/O and yields overall end-to-end encoding and I/O time that is less than that required to write the raw data with MPI collective I/O.