Systematic evaluation of B-cell clonal family inference approaches

The reconstruction of clonal families (CFs) in B-cell receptor (BCR) repertoire analysis is a crucial step to understand the adaptive immune system and how it responds to antigens. The BCR repertoire of an individual is formed throughout life and is diverse due to several factors such as gene recombination and somatic hypermutation. The use of Adaptive Immune Receptor Repertoire sequencing (AIRR-seq) using next generation sequencing enabled the generation of full BCR repertoires that also include rare CFs. The reconstruction of CFs from AIRR-seq data is challenging and several approaches have been developed to solve this problem. Currently, most methods use the heavy chain (HC) only, as it is more variable than the light chain (LC). CF reconstruction options include the definition of appropriate sequence similarity measures, the use of shared mutations among sequences, and the possibility of reconstruction without preliminary clustering based on V- and J-gene annotation. In this study, we aimed to systematically evaluate different approaches for CF reconstruction and to determine their impact on various outcome measures such as the number of CFs derived, the size of the CFs, and the accuracy of the reconstruction. The methods were compared to each other and to a method that groups sequences based on identical junction sequences and another method that only determines subclones. We found that after accounting for data set variability, in particular sequencing depth and mutation load, the reconstruction approach has an impact on part of the outcome measures, including the number of CFs. Simulations indicate that unique junctions and subclones should not be used as substitutes for CF and that more complex methods do not outperform simpler methods. Also, we conclude that different approaches differ in their ability to correctly reconstruct CFs when not considering the LC and to identify shared CFs. The results showed the effect of different approaches on the reconstruction of CFs and highlighted the importance of choosing an appropriate method.

TCRβ clones in muscle tissue share structural features in patients with idiopathic inflammatory myopathy and are associated with disease activity

Objectives

To characterize the T cell receptor (TCRβ) repertoire in peripheral blood and muscle tissues of treatment naïve patients with newly diagnosed idiopathic inflammatory myopathies (IIMs).

Methods

High throughput RNA sequencing of the TCRβ chain was performed in peripheral blood and muscle tissue in twenty newly-diagnosed treatment-naïve IIM patients (9 DM, 5 NM/OM, 5 IMNM and 1 ASyS) and healthy controls. Results thereof were correlated with markers of disease activity.

Results

Muscle tissue of IIM patients shows more expansion of TCRβ clones and decreased diversity when compared to peripheral blood of IIM as well as healthy controls (both p=0.0001). Several expanded TCRβ clones in muscle are tissue restricted and cannot be retrieved in peripheral blood. These clones have significantly longer CDR3 regions when compared to clones (also) found in circulation (p=0.0002), while their CDR3 region is more hydrophobic (p<0.01). Network analysis shows that clonal TCRβ signatures are shared between patients. Increased clonal expansion in muscle tissue is significantly correlated with increased CK levels (p=0.03), while it tends to correlate with decreased muscle strength (p=0.08).

Conclusion

Network analysis of clones in muscle of IIM patients shows shared clusters of sequences across patients. Muscle-restricted CDR3 TCRβ clones show specific structural features in their T cell receptor. Our results indicate that clonal TCRβ expansion in muscle tissue might be associated with disease activity. Collectively, these findings support a role for specific clonal T cell responses in muscle tissue in the pathogenesis of the IIM subtypes studied.

B-cell receptor profiling before and after IVIG monotherapy in newly diagnosed idiopathic inflammatory myopathies

OBJECTIVE

To unravel B-cell receptor (BcR) characteristics in muscle tissues and peripheral blood and gain more insight into BcR repertoire changes in peripheral blood in idiopathic inflammatory myopathies (IIMs), and study how this correlates to the clinical response to IVIG.

METHODS

Nineteen treatment-naive patients with newly diagnosed IIM were prospectively treated with IVIG monotherapy. RNA-based BcR repertoire sequencing was performed in muscle biopsies collected before, and in peripheral blood (PB) collected before and nine weeks after IVIG treatment. Results were correlated to patients' clinical improvement based on the total improvement score (TIS).

RESULTS

Prior to IVIG treatment, BcR clones found in muscle tissue could be retrieved in peripheral blood. Nine weeks after IVIG treatment, new patient-specific dominant BcR clones appeared in peripheral blood while pre-treatment dominant BcR clones disappeared. The cumulative frequency of all dominant BcR clones before treatment was significantly higher in individuals who responded to IVIG compared with those who did not respond to IVIG, and correlated with a higher CK. During follow-up, a decrease in the cumulative frequency of all dominant clones correlated with a higher TIS.

CONCLUSION

In treatment-naive patients with newly diagnosed IIM, muscle tissue and peripheral blood share expanded BcR clones. In our study a higher cumulative frequency of dominant BcR clones in blood before treatment was associated with a higher CK and better treatment response, suggesting that response to IVIG may depend on the composition of the pre-treatment BcR repertoire.

A public antibody class recognizes an S2 epitope exposed on open conformations of SARS-CoV-2 spike

Delineating the origins and properties of antibodies elicited by SARS-CoV-2 infection and vaccination is critical for understanding their benefits and potential shortcomings. Therefore, we investigate the SARS-CoV-2 spike (S)-reactive B cell repertoire in unexposed individuals by flow cytometry and single-cell sequencing. We show that ∼82% of SARS-CoV-2 S-reactive B cells harbor a naive phenotype, which represents an unusually high fraction of total human naive B cells (∼0.1%). Approximately 10% of these naive S-reactive B cells share an IGHV1-69/IGKV3-11 B cell receptor pairing, an enrichment of 18-fold compared to the complete naive repertoire. Following SARS-CoV-2 infection, we report an average 37-fold enrichment of IGHV1-69/IGKV3-11 B cell receptor pairing in the S-reactive memory B cells compared to the unselected memory repertoire. This class of B cells targets a previously undefined non-neutralizing epitope on the S2 subunit that becomes exposed on S proteins used in approved vaccines when they transition away from the native pre-fusion state because of instability. These findings can help guide the improvement of SARS-CoV-2 vaccines.

In rheumatoid arthritis inflamed joints share dominant patient-specific B-cell clones

Background

In patients with rheumatoid arthritis (RA) different joints were shown to share the same dominant T-cell clones, suggesting shared characteristics of the inflammatory process and indicating that strategies to selectively target the antigen receptor might be feasible. Since T- and B-lymphocytes closely interact in adaptive responses, we analysed to what extent different joints also share dominant B-cell clones.

Methods

In 11 RA patients, quantitative B-cell receptor (BCR) repertoire analysis was performed in simultaneously obtained samples from inflamed synovial tissue (ST) from distinct locations within one joint, from multiple joints, from synovial fluid (SF) and peripheral blood (PB).

Results

ST biopsies from different locations in the same joint showed clear overlap in the top-25 dominant BCR clones (16.7%, SD 12.5), in the same range as the overlap between ST and SF in the same joint (8.0%, SD 8.8) and the overlap between ST-ST between different joints (9.1%, SD 8.2), but clearly higher than the overlap between ST and PB (1.7%, SD 2.4; p<0.05) and SF and PB (2.7%, SD 4.1; p<0.05). Interestingly, these figures were substantially lower than the overlap observed in previous T-cell clonality studies.

Conclusions

We conclude that in RA BCR clonal responses may be more localized than TCR clonal responses, pointing to antigen-selective influx, proliferation and/or maturation of B-cells. B lineage cells in the SF may adequately represent the dominant BCR clones of the ST, which is in contrast to T-cells. Collectively, the presence of shared B- and especially T-cells in different joints from the same patient suggests that approaches might be feasible that aim to develop antigen-receptor specific targeting of lymphocyte clones in RA as an alternative to more generalized immunosuppressive strategies.

Longitudinal analysis of anti-drug antibody development in multiple sclerosis patients treated with interferon beta-1a (Rebif™) using B cell receptor repertoire analysis

A significant proportion of multiple sclerosis (MS) patients treated with interferon beta-1a (Rebif™) develop anti-drug antibodies (ADA) with a negative impact on treatment efficacy. We hypothesized that high-throughput B-cell receptor (BCR) repertoire analysis could be used to predict and monitor ADA development. To study this we analyzed 228 peripheral blood samples from 68 longitudinally followed patients starting on interferon beta-1a. Our results show that whole blood BCR analysis does not reflect, and does not predict ADA development in MS patients treated with interferon beta-1a. We propose that BCR analysis of phenotypically selected cell subsets or tissues might be more informative.

Proteogenomic analysis of the autoreactive B cell repertoire in blood and tissues of patients with Sjögren's syndrome

Objective

To comparatively analyse the aberrant affinity maturation of the antinuclear and rheumatoid factor (RF) B cell repertoires in blood and tissues of patients with Sjögren’s syndrome (SjS) using an integrated omics workflow.

Methods

Peptide sequencing of anti-Ro60, anti-Ro52, anti-La and RF was combined with B cell repertoire analysis at the DNA, RNA and single cell level in blood B cell subsets, affected salivary gland and extranodal marginal zone lymphomas of mucosa-associated lymphoid tissue (MALT) of patients with SjS.

Results

Affected tissues contained anti-Ro60, anti-Ro52, anti-La and RF clones as a small part of a polyclonal infiltrate. Anti-Ro60, anti-La and anti-Ro52 clones outnumbered RF clones. MALT lymphoma tissues contained monoclonal RF expansions. Autoreactive clones were not selected from a restricted repertoire in a circulating B cell subset. The antinuclear antibody (ANA) repertoires displayed similar antigen-dependent and immunoglobulin (Ig) G1-directed affinity maturation. RF clones displayed antigen-dependent, IgM-directed and more B cell receptor integrity-dependent affinity maturation. This coincided with extensive intra-clonal diversification in RF-derived lymphomas. Regeneration of clinical disease manifestations after rituximab coincided with large RF clones, which not necessarily belonged to the lymphoma clone, that displayed continuous affinity maturation and intra-clonal diversification.

Conclusion

The ANA and RF repertoires in patients with SjS display tissue-restricted, antigen-dependent and divergent affinity maturation. Affinity maturation of RF clones deviates further during RF clone derived lymphomagenesis and during regeneration of the autoreactive repertoire after temporary disruption by rituximab. These data give insight into the molecular mechanisms of autoreactive inflammation in SjS, assist MALT lymphoma diagnosis and allow tracking its response to rituximab.

Characterization and Monitoring of Antigen-Responsive T Cell Clones Using T Cell Receptor Gene Expression Analysis

High-throughput T-cell receptor repertoire sequencing constitutes a powerful tool to study T cell responses at the clonal level. However, it does not give information on the functional phenotype of the responding clones and lacks a statistical framework for quantitative evaluation. To overcome this, we combined datasets from different experiments, all starting from the same blood samples. We used a novel, sensitive, UMI-based protocol to perform repertoire analysis on experimental replicates. Applying established bioinformatic routines for transcriptomic expression analysis we explored the dynamics of antigen-induced clonal expansion after in vitro stimulation, identified antigen-responsive clones, and confirmed their activation status using the expression of activation markers upon antigen re-challenge. We demonstrate that the addition of IL-4 after antigen stimulation drives the expansion of T cell clones encoding unique receptor sequences. We show that our approach represents a scalable, high-throughput immunological tool, which can be used to identify and characterize antigen-responsive T cells at clonal level.

Novel loci for childhood body mass index and shared heritability with adult cardiometabolic traits

The genetic background of childhood body mass index (BMI), and the extent to which the well-known associations of childhood BMI with adult diseases are explained by shared genetic factors, are largely unknown. We performed a genome-wide association study meta-analysis of BMI in 61,111 children aged between 2 and 10 years. Twenty-five independent loci reached genome-wide significance in the combined discovery and replication analyses. Two of these, located near NEDD4L and SLC45A3, have not previously been reported in relation to either childhood or adult BMI. Positive genetic correlations of childhood BMI with birth weight and adult BMI, waist-to-hip ratio, diastolic blood pressure and type 2 diabetes were detected (Rg ranging from 0.11 to 0.76, P-values <0.002). A negative genetic correlation of childhood BMI with age at menarche was observed. Our results suggest that the biological processes underlying childhood BMI largely, but not completely, overlap with those underlying adult BMI. The well-known observational associations of BMI in childhood with cardio-metabolic diseases in adulthood may reflect partial genetic overlap, but in light of previous evidence, it is also likely that they are explained through phenotypic continuity of BMI from childhood into adulthood.

A trans-ancestral meta-analysis of genome-wide association studies reveals loci associated with childhood obesity

Although hundreds of genome-wide association studies-implicated loci have been reported for adult obesity-related traits, less is known about the genetics specific for early-onset obesity and with only a few studies conducted in non-European populations to date. Searching for additional genetic variants associated with childhood obesity, we performed a trans-ancestral meta-analysis of 30 studies consisting of up to 13 005 cases (≥95th percentile of body mass index (BMI) achieved 2-18 years old) and 15 599 controls (consistently <50th percentile of BMI) of European, African, North/South American and East Asian ancestry. Suggestive loci were taken forward for replication in a sample of 1888 cases and 4689 controls from seven cohorts of European and North/South American ancestry. In addition to observing 18 previously implicated BMI or obesity loci, for both early and late onset, we uncovered one completely novel locus in this trans-ancestral analysis (nearest gene, METTL15). The variant was nominally associated with only the European subgroup analysis but had a consistent direction of effect in other ethnicities. We then utilized trans-ancestral Bayesian analysis to narrow down the location of the probable causal variant at each genome-wide significant signal. Of all the fine-mapped loci, we were able to narrow down the causative variant at four known loci to fewer than 10 single nucleotide polymorphisms (SNPs) (FAIM2, GNPDA2, MC4R and SEC16B loci). In conclusion, an ethnically diverse setting has enabled us to both identify an additional pediatric obesity locus and further fine-map existing loci.

N-Glycosylation Site Analysis of Citrullinated Antigen-Specific B-Cell Receptors Indicates Alternative Selection Pathways During Autoreactive B-Cell Development

Many autoimmune diseases are hallmarked by autoreactive B and plasma cell responses that are directly or indirectly involved in disease pathogenesis. These B-cell responses show large variability between diseases, both in terms of the secreted autoantibody repertoire and the dynamics and characteristics of the underlying B-cell responses. Hence, different mechanisms have been proposed to explain the emergence of autoreactive B cells in an otherwise self-tolerant immune system. Notably, most mechanistic insights have been obtained from murine studies using models harboring genetic modifications of B and T cells. Given recent technological advances that have rendered autoreactive human B cells accessible for analysis, we here discuss the phenomenon of extensive N-glycosylation of the B-cell receptor (BCR) variable domain of a prototypic human autoreactive B-cell response and its potential role in the generation of autoimmunity. Anti-citrullinated protein antibodies (ACPA) hallmark the most disease-specific autoimmune response in Rheumatoid Arthritis (RA). Remarkably, ACPA-IgG are heavily N-glycosylated in the variable domain due to somatic mutations that generate abundant N-glycosylation consensus sequences. These sites, obtained from full-length BCR sequences of ACPA-expressing B cells from 12 ACPA-positive RA patients, were here analyzed in detail. Sites that required a single nucleotide mutation to be generated were defined as single somatic hypermutation (s-SHM) sites, whereas sites requiring multiple mutations were defined as m-SHM sites. IgG sequences of 12 healthy donors were used as control. Computational modeling of the germinal center reaction (CLONE algorithm) was used with the germline counterparts of ACPA-IgG heavy chain (HC) sequences to simulate the germinal center response. Our analyses revealed an abundance of N-glycosylation sites in ACPA-IgG HC that frequently required multiple mutations and predominated in specific positions. Based on these data, and taking into account recent insights into the dynamics of the ACPA-response during disease development, we here discuss the hypothesis that N-glycosylation sites in ACPA-IgG variable domains could lead to alternative, possibly antibody affinity-independent selection forces. Presumably, this occurs during germinal center responses allowing these B cells to escape from putative tolerance checkpoints, thereby driving autoreactive B cell development in the pathogenesis of RA.

Non-response to rituximab therapy in rheumatoid arthritis is associated with incomplete disruption of the B cell receptor repertoire

OBJECTIVE

To gain more insight into the dynamics of lymphocyte depletion and develop new predictors of clinical response to rituximab in rheumatoid arthritis (RA).

METHODS

RNA-based next-generation sequencing was used to analyse the B cell receptor (BCR) repertoire in peripheral blood and synovial tissue samples collected from 24 seropositive patients with RA treated with rituximab. Clonal expansion, mutation load and clonal overlap were assessed in samples collected before, at week 4 and at week 16 or 24 after treatment and correlated to the patients' clinical response.

RESULTS

After 4 weeks of rituximab-induced B cell depletion, the peripheral blood BCR repertoire of treated patients consisted of fewer, more dominant and more mutated BCR clones. No significant changes in the synovial tissue BCR repertoire were detected until week 16 post-treatment, when a reduced clonal overlap with baseline and an increased mutation load were observed. In patients who were non-responders at month 3 (n=5) using the European League Against Rheumatism response criteria, peripheral blood samples taken at week 4 after rituximab treatment showed more dominant clones compared with moderate responders (n=9) (median (IQR): 36 (27-52) vs 18 (16-26); p<0.01) and more clonal overlap with the baseline (median (IQR): 5% (2%-20%) vs 0% (0%-0%); p≤0.01).

CONCLUSION

Significant changes in BCR clonality are observed in peripheral blood of patients 4 weeks after rituximab treatment, while changes in synovial tissue were observed at later time points. Incomplete depletion of the dominant baseline peripheral blood BCR repertoire in the first month of treatment might predict clinical non-response at 3 months.

Maternal and fetal genetic effects on birth weight and their relevance to cardio-metabolic risk factors

Birth weight variation is influenced by fetal and maternal genetic and non-genetic factors, and has been reproducibly associated with future cardio-metabolic health outcomes. In expanded genome-wide association analyses of own birth weight (n = 321,223) and offspring birth weight (n = 230,069 mothers), we identified 190 independent association signals (129 of which are novel). We used structural equation modeling to decompose the contributions of direct fetal and indirect maternal genetic effects, then applied Mendelian randomization to illuminate causal pathways. For example, both indirect maternal and direct fetal genetic effects drive the observational relationship between lower birth weight and higher later blood pressure: maternal blood pressure-raising alleles reduce offspring birth weight, but only direct fetal effects of these alleles, once inherited, increase later offspring blood pressure. Using maternal birth weight-lowering genotypes to proxy for an adverse intrauterine environment provided no evidence that it causally raises offspring blood pressure, indicating that the inverse birth weight-blood pressure association is attributable to genetic effects, and not to intrauterine programming.

Computational Model Reveals Limited Correlation between Germinal Center B-Cell Subclone Abundancy and Affinity: Implications for Repertoire Sequencing

Immunoglobulin repertoire sequencing has successfully been applied to identify expanded antigen-activated B-cell clones that play a role in the pathogenesis of immune disorders. One challenge is the selection of the Ag-specific B cells from the measured repertoire for downstream analyses. A general feature of an immune response is the expansion of specific clones resulting in a set of subclones with common ancestry varying in abundance and in the number of acquired somatic mutations. The expanded subclones are expected to have BCR affinities for the Ag higher than the affinities of the naive B cells in the background population. For these reasons, several groups successfully proceeded or suggested selecting highly abundant subclones from the repertoire to obtain the Ag-specific B cells. Given the nature of affinity maturation one would expect that abundant subclones are of high affinity but since repertoire sequencing only provides information about abundancies, this can only be verified with additional experiments, which are very labor intensive. Moreover, this would also require knowledge of the Ag, which is often not available for clinical samples. Consequently, in general we do not know if the selected highly abundant subclone(s) are also the high(est) affinity subclones. Such knowledge would likely improve the selection of relevant subclones for further characterization and Ag screening. Therefore, to gain insight in the relation between subclone abundancy and affinity, we developed a computational model that simulates affinity maturation in a single GC while tracking individual subclones in terms of abundancy and affinity. We show that the model correctly captures the overall GC dynamics, and that the amount of expansion is qualitatively comparable to expansion observed from B cells isolated from human lymph nodes. Analysis of the fraction of high- and low-affinity subclones among the unexpanded and expanded subclones reveals a limited correlation between abundancy and affinity and shows that the low abundant subclones are of highest affinity. Thus, our model suggests that selecting highly abundant subclones from repertoire sequencing experiments would not always lead to the high(est) affinity B cells. Consequently, additional or alternative selection approaches need to be applied.

Somatic Variation of T-Cell Receptor Genes Strongly Associate with HLA Class Restriction

Every person carries a vast repertoire of CD4⁺ T-helper cells and CD8⁺ cytotoxic T cells for a healthy immune system. Somatic VDJ recombination at genomic loci that encode the T-cell receptor (TCR) is a key step during T-cell development, but how a single T cell commits to become either CD4⁺ or CD8⁺ is poorly understood. To evaluate the influence of TCR sequence variation on CD4⁺/CD8⁺ lineage commitment, we sequenced rearranged TCRs for both α and β chains in naïve T cells isolated from healthy donors and investigated gene segment usage and recombination patterns in CD4⁺ and CD8⁺ T-cell subsets. Our data demonstrate that most V and J gene segments are strongly biased in the naïve CD4⁺ and CD8⁺ subsets with some segments increasing the odds of being CD4⁺ (or CD8⁺) up to five-fold. These V and J gene associations are highly reproducible across individuals and independent of classical HLA genotype, explaining ~11% of the observed variance in the CD4⁺ vs. CD8⁺ propensity. In addition, we identified a strong independent association of the electrostatic charge of the complementarity determining region 3 (CDR3) in both α and β chains, where a positively charged CDR3 is associated with CD4⁺ lineage and a negatively charged CDR3 with CD8⁺ lineage. Our findings suggest that somatic variation in different parts of the TCR influences T-cell lineage commitment in a predominantly additive fashion. This notion can help delineate how certain structural features of the TCR-peptide-HLA complex influence thymic selection.

Viral metagenomics in drug-naïve, first-onset schizophrenia patients with prominent negative symptoms

Although several studies suggest a virus or (endogenous) retrovirus involvement at the time of onset of schizophrenia, the unequivocal identification of one or more infectious agents, by means of an undirected catch-all technique, has never been conducted. In this study VIDISCA, a virus discovery method, was used in combination with Roche-454 high-throughput sequencing as a tool to determine the possible presence of viruses (known or unknown) in blood of first-onset drugs-naïve schizophrenic patients with prominent negative symptoms. Two viruses (the Anellovirus Torque Teno virus and GB virus C) were detected. Both viruses are commonly found in healthy individuals and no clear link with disease was ever established. Viruses from the family Anelloviridae were also identified in the control population (4.8%). Besides, one patient sample was positive for human endogenous retroviruses type K (HML-2) RNA but no specific predominant strain was detected, instead 119 different variants were found. In conclusion, these findings indicate no evidence for viral or endogenous retroviral involvement in sera at the time of onset of schizophrenia.

Streptococcus pneumoniae arginine synthesis genes promote growth and virulence in pneumococcal meningitis

Streptococcus pneumoniae (pneumococcus) is a major human pathogen causing pneumonia, sepsis and bacterial meningitis. Using a clinical phenotype based approach with bacterial whole-genome sequencing we identified pneumococcal arginine biosynthesis genes to be associated with outcome in patients with pneumococcal meningitis. Pneumococci harboring these genes show increased growth in human blood and cerebrospinal fluid (CSF). Mouse models of meningitis and pneumonia showed that pneumococcal strains without arginine biosynthesis genes were attenuated in growth or cleared, from lung, blood and CSF. Thus, S. pneumoniae arginine synthesis genes promote growth and virulence in invasive pneumococcal disease.

Pro-apoptotic protein Noxa regulates memory T cell population size and protects against lethal immunopathology

Memory T cells form a highly specific defense layer against reinfection with previously encountered pathogens. In addition, memory T cells provide protection against pathogens that are similar, but not identical to the original infectious agent. This is because each T cell response harbors multiple clones with slightly different affinities, thereby creating T cell memory with a certain degree of diversity. Currently, the mechanisms that control size, diversity, and cross-reactivity of the memory T cell pool are incompletely defined. Previously, we established a role for apoptosis, mediated by the BH3-only protein Noxa, in controlling diversity of the effector T cell population. This function might positively or negatively impact T cell memory in terms of function, pool size, and cross-reactivity during recall responses. Therefore, we investigated the role of Noxa in T cell memory during acute and chronic infections. Upon influenza infection, Noxa(-/-) mice generate a memory compartment of increased size and clonal diversity. Reinfection resulted in an increased recall response, whereas cross-reactive responses were impaired. Chronic infection of Noxa(-/-) mice with mouse CMV resulted in enhanced memory cell inflation, but no obvious pathology. In contrast, in a model of continuous, high-level T cell activation, reduced apoptosis of activated T cells rapidly led to severe organ pathology and premature death in Noxa-deficient mice. These results establish Noxa as an important regulator of the number of memory cells formed during infection. Chronic immune activation in the absence of Noxa leads to excessive accumulation of primed cells, which may result in severe pathology.

Provenance for distributed biomedical workflow execution

Scientific research has become very data and compute intensive because of the progress in data acquisition and measurement devices, which is particularly true in Life Sciences. To cope with this deluge of data, scientists use distributed computing and storage infrastructures. The use of such infrastructures introduces by itself new challenges to the scientists in terms of proper and efficient use. Scientific workflow management systems play an important role in facilitating the use of the infrastructure by hiding some of its complexity. Although most scientific workflow management systems are provenance-aware, not all of them come with provenance functionality out of the box. In this paper we describe the improvement and integration of a provenance system into an e-infrastructure for biomedical research based on the MOTEUR workflow management system. The main contributions of the paper are: presenting an OPM implementation using relational database backend for the provenance store, providing an e-infrastructure with a comprehensive provenance system, defining a generic approach to provenance implementation, potentially suitable for other workflow systems and application domains and demonstrating the value of this system based on use cases presenting the provenance data through a user-friendly web interface.

Looking ultra deep: short identical sequences and transcriptional slippage

Studying transcriptomes by ultra deep sequencing provides an in-depth picture of transcriptional regulation and it facilitates the detection of rare transcriptional events. Using ultra deep sequencing of amplicons we identified known isoforms and also various new low frequency variants. Most of these variants likely involve the splicing machinery except for two events that we named variations affecting multiple exons, which are mainly deletions affecting parts of adjacent exons and intra-exonic deletions. Both events involve short identical sequences of 1 to 8 nucleotides at the junction and canonical splice sites are missing. They were identified in different genes and species at very low frequencies. We excluded that they are an artifact of PCR, sequencing, or reverse transcription. We propose that these variants represent intramolecular slippage events that require short identical sequences for reannealing of dissociated transcripts.

Initial steps towards a production platform for DNA sequence analysis on the grid

Background

Bioinformatics is confronted with a new data explosion due to the availability of high throughput DNA sequencers. Data storage and analysis becomes a problem on local servers, and therefore it is needed to switch to other IT infrastructures. Grid and workflow technology can help to handle the data more efficiently, as well as facilitate collaborations. However, interfaces to grids are often unfriendly to novice users.

Results

In this study we reused a platform that was developed in the VL-e project for the analysis of medical images. Data transfer, workflow execution and job monitoring are operated from one graphical interface. We developed workflows for two sequence alignment tools (BLAST and BLAT) as a proof of concept. The analysis time was significantly reduced. All workflows and executables are available for the members of the Dutch Life Science Grid and the VL-e Medical virtual organizations All components are open source and can be transported to other grid infrastructures.

Conclusions

The availability of in-house expertise and tools facilitates the usage of grid resources by new users. Our first results indicate that this is a practical, powerful and scalable solution to address the capacity and collaboration issues raised by the deployment of next generation sequencers. We currently adopt this methodology on a daily basis for DNA sequencing and other applications. More information and source code is available via http://www.bioinformaticslaboratory.nl/

Human T-cell memory consists mainly of unexpanded clones

The immune system is able to respond to millions of antigens using adaptive receptors, including the alphabeta-T-cell receptor (TCR). Upon antigen encounter a T-cell may proliferate to produce a clone of TCR-identical cells, which develop a memory phenotype. Previous studies suggested that most memory clones are clearly expanded. In accordance, the beta-chain repertoire of T-cell memory subsets was reported to be 10 times less diverse than those of naive subsets, reflecting stringent selection. However, due to technological limitations detailed information was lacking regarding the size of clonal expansions and the diversity of the TCR-repertoire in naive and memory T-cell populations. Here, using high-throughput sequencing, we show that the memory repertoire in human peripheral blood contains only few expanded clones and consists mainly of low frequency clones. Additionally, the memory repertoire is much more diverse than expected. In two healthy persons we observed that only 2-7% of the CD4 and CD8 memory clones found were clearly expanded. In line with this observation we show that the beta-chains repertoire size of the CD4 memory compartment is only two times smaller, and that of the CD8 memory compartment is only 3-10 times smaller than the naive compartments. Our results show that the T-cell memory compartment has a very different distribution of clones than anticipated. This has important implications for the current dogma of immunological memory, and changes the interpretation of repertoire aberrations in (patho-)physiological situations such as ageing and auto-immunity. It raises new questions on the factors that steer maturation of memory phenotype and determine the size of memory clones.

The construction of genome-based transcriptional units

Gene-oriented sequence clusters (transcriptional units) have found many applications in genomics research including the construction of transcriptome maps and identification of splice variants. We developed a new method to construct transcriptional that uses the genomic sequence as a template. We present and discuss our method in detail together with an evaluation of the transcriptional units for human. We constructed 33,007 and 27,792 transcriptional units for human and mouse, respectively. The sensitivity (81%) and specificity (90%) of our method compares favorably to other established methods. We evaluated the representation of experimentally validated and predicted intergenic spliced transcripts in humans and show that we correctly represent a large fraction of these cases by single transcriptional units. Our method performs well, but the evaluation of the final set of transcriptional units show that improvements to the algorithm are still possible. However, because the precise number and types of errors are difficult to track, it is not obvious how to significantly improve the algorithm. We believe that ongoing research efforts are necessary to further improve current methods. This should include detailed documentation, comparison, and evaluation of current methods.

The human transcriptome map reveals extremes in gene density, intron length, GC content, and repeat pattern for domains of highly and weakly expressed genes

The chromosomal gene expression profiles established by the Human Transcriptome Map (HTM) revealed a clustering of highly expressed genes in about 30 domains, called ridges. To physically characterize ridges, we constructed a new HTM based on the draft human genome sequence (HTMseq). Expression of 25,003 genes can be analyzed online in a multitude of tissues (http://bioinfo.amc.uva.nl/HTMseq). Ridges are found to be very gene-dense domains with a high GC content, a high SINE repeat density, and a low LINE repeat density. Genes in ridges have significantly shorter introns than genes outside of ridges. The HTMseq also identifies a significant clustering of weakly expressed genes in domains with fully opposite characteristics (antiridges). Both types of domains are open to tissue-specific expression regulation, but the maximal expression levels in ridges are considerably higher than in antiridges. Ridges are therefore an integral part of a higher order structure in the genome related to transcriptional regulation.