Diversity and Clonality of T Cell Receptor Repertoire and Antigen Specificities in Small Joints of Early Rheumatoid Arthritis.

OBJECTIVE: CD4+ T cells are implicated in rheumatoid arthritis (RA) pathology from the strong association between RA and certain HLA class II gene variants. This study was undertaken to examine the synovial T cell receptor (TCR) repertoire, T cell phenotypes, and T cell specificities in small joints of RA patients at time of diagnosis before therapeutic intervention. METHODS: Sixteen patients, of whom 11 patients were anti-citrullinated protein antibody (ACPA)-positive and 5 patients were ACPA-, underwent ultrasound-guided synovial biopsy of a small joint (n = 13) or arthroscopic synovial biopsy of a large joint (n = 3), followed by direct sorting of single T cells for paired sequencing of the αß TCR together with flow cytometry analysis. TCRs from expanded CD4+ T cell clones of 4 patients carrying an HLA-DRB1*04:01 allele were artificially reexpressed to study antigen specificity. RESULTS: T cell analysis demonstrated CD4+ dominance and the presence of peripheral helper T-like cells in both patient groups. We identified >4,000 unique TCR sequences, as well as 225 clonal expansions. Additionally, T cells with double α-chains were a recurring feature. We identified a biased gene usage of the Vß chain segment TRBV20-1 in CD4+ cells from ACPA+ patients. In vitro stimulation of T cell lines expressing selected TCRs with an extensive panel of citrullinated and viral peptides identified several different virus-specific TCRs (e.g., human cytomegalovirus and human herpesvirus 2). Still, the majority of clones remained orphans with unknown specificity. CONCLUSION: Minimally invasive biopsies of the RA synovium allow for single-cell TCR sequencing and phenotyping. Clonally expanded, viral-reactive T cells account for part of the diverse CD4+ T cell repertoire. TRBV20-1 bias in ACPA+ patients suggests recognition of common antigens.

Comparative analysis of antibody- and lipid-based multiplexing methods for single-cell RNA-seq.

BACKGROUND: Multiplexing of samples in single-cell RNA-seq studies allows a significant reduction of the experimental costs, straightforward identification of doublets, increased cell throughput, and reduction of sample-specific batch effects. Recently published multiplexing techniques using oligo-conjugated antibodies or -lipids allow barcoding sample-specific cells, a process called "hashing." RESULTS: Here, we compare the hashing performance of TotalSeq-A and -C antibodies, custom synthesized lipids and MULTI-seq lipid hashes in four cell lines, both for single-cell RNA-seq and single-nucleus RNA-seq. We also compare TotalSeq-B antibodies with CellPlex reagents (10x Genomics) on human PBMCs and TotalSeq-B with different lipids on primary mouse tissues. Hashing efficiency was evaluated using the intrinsic genetic variation of the cell lines and mouse strains. Antibody hashing was further evaluated on clinical samples using PBMCs from healthy and SARS-CoV-2 infected patients, where we demonstrate a more affordable approach for large single-cell sequencing clinical studies, while simultaneously reducing batch effects. CONCLUSIONS: Benchmarking of different hashing strategies and computational pipelines indicates that correct demultiplexing can be achieved with both lipid- and antibody-hashed human cells and nuclei, with MULTISeqDemux as the preferred demultiplexing function and antibody-based hashing as the most efficient protocol on cells. On nuclei datasets, lipid hashing delivers the best results. Lipid hashing also outperforms antibodies on cells isolated from mouse brain. However, antibodies demonstrate better results on tissues like spleen or lung.

Author Correction: Gene-associated markers provide tools for tackling illegal fishing and false eco-certification.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

DNA barcoding fishes from the Congo and the Lower Guinean provinces: Assembling a reference library for poorly inventoried fauna.

The Congolese and Lower Guinean ichthyological provinces are understudied hotspots of the global fish diversity. Here, we barcoded 741 specimens from the Lower and Middle Congo River and from three major drainage basins of the Lower Guinean ichthyological province, Kouilou-Niari, Nyanga and Ogowe. We identified 194 morphospecies belonging to 82 genera and 25 families. Most morphospecies (92.8%) corresponded to distinct clusters of DNA barcodes. Of the four morphospecies present in both neighbouring ichthyological provinces, only one showed DNA barcode divergence <2.5%. A small fraction of the fishes barcoded here (12.9% of the morphospecies and 16.1% of the barcode clusters representing putative species) were also barcoded in a previous large-scale DNA analysis of freshwater fishes of the Lower Congo published in 2011 (191 specimens, 102 morphospecies). We compared species assignments before and after taxonomic updates and across studies performed by independent research teams and observed that most cases of inconsistent species assignments were due to unknown diversity (undescribed species and unknown intraspecific variation). Our results report more than 17 putative new species and show that DNA barcode data provide a measure of genetic variability that facilitates the inventory of underexplored ichthyofaunae. However, taxonomic scrutiny, associated with revisions and new species descriptions, is indispensable to delimit species and build a coherent reference library.

ACE-inhibition induces a cardioprotective transcriptional response in the metabolic syndrome heart.

Cardiovascular disease associated with metabolic syndrome has a high prevalence, but the mechanistic basis of metabolic cardiomyopathy remains poorly understood. We characterised the cardiac transcriptome in a murine metabolic syndrome (MetS) model (LDLR-/-; ob/ob, DKO) relative to the healthy, control heart (C57BL/6, WT) and the transcriptional changes induced by ACE-inhibition in those hearts. RNA-Seq, differential gene expression and transcription factor analysis identified 288 genes differentially expressed between DKO and WT hearts implicating 72 pathways. Hallmarks of metabolic cardiomyopathy were increased activity in integrin-linked kinase signalling, Rho signalling, dendritic cell maturation, production of nitric oxide and reactive oxygen species in macrophages, atherosclerosis, LXR-RXR signalling, cardiac hypertrophy, and acute phase response pathways. ACE-inhibition had a limited effect on gene expression in WT (55 genes, 23 pathways), and a prominent effect in DKO hearts (1143 genes, 104 pathways). In DKO hearts, ACE-I appears to counteract some of the MetS-specific pathways, while also activating cardioprotective mechanisms. We conclude that MetS and control murine hearts have unique transcriptional profiles and exhibit a partially specific transcriptional response to ACE-inhibition.

Describing the Reportable Range Is Important for Reliable Treatment Decisions: A Multiple Laboratory Study for Molecular Tumor Profiling Using Next-Generation Sequencing.

Because interpretation of next-generation sequencing (NGS) data remains challenging, optimization of the NGS process is needed to obtain correct sequencing results. Therefore, extensive validation and continuous monitoring of the quality is essential. NGS performance was compared with traditional detection methods and technical quality of nine NGS technologies was assessed. First, nine formalin-fixed, paraffin-embedded patient samples were analyzed by 114 laboratories by using different detection methods. No significant differences in performance were observed between analyses with NGS and traditional techniques. Second, two DNA control samples were analyzed for a selected number of variants by 26 participants with the use of nine different NGS technologies. Quality control metrics were analyzed from raw data files and a survey about routine procedures. Results showed large differences in coverages, but observed variant allele frequencies in raw data files were in line with predefined variant allele frequencies. Many false negative results were found because of low-quality regions, which were not reported as such. It is recommended to disclose the reportable range, the fraction of targeted genomic regions for which calls of acceptable quality can be generated, to avoid any errors in therapy decisions. NGS can be a reliable technique, only if essential quality control during analysis is applied and reported.

Neurogenomic Profiling Reveals Distinct Gene Expression Profiles Between Brain Parts That Are Consistent in Ophthalmotilapia Cichlids.

The detection of external and internal cues alters gene expression in the brain which in turn may affect neural networks that underly behavioral responses. Previous studies have shown that gene expression profiles differ between major brain regions within individuals and between species with different morphologies, cognitive abilities and/or behaviors. A detailed description of gene expression in all macroanatomical brain regions and in species with similar morphologies and behaviors is however lacking. Here, we dissected the brain of two cichlid species into six macroanatomical regions. Ophthalmotilapia nasuta and O. ventralis have similar morphology and behavior and occasionally hybridize in the wild. We use 3' mRNA sequencing and a stage-wise statistical testing procedure to identify differential gene expression between females that were kept in a social setting with other females. Our results show that gene expression differs substantially between all six brain parts within species: out of 11,577 assessed genes, 8,748 are differentially expressed (DE) in at least one brain part compared to the average expression of the other brain parts. At most 16% of these DE genes have |log2FC| significantly higher than two. Functional differences between brain parts were consistent between species. The majority (61-79%) of genes that are DE in a particular brain part were shared between both species. Only 32 genes show significant differences in fold change across brain parts between species. These genes are mainly linked to transport, transmembrane transport, transcription (and its regulation) and signal transduction. Moreover, statistical equivalence testing reveals that within each comparison, on average 89% of the genes show an equivalent fold change between both species. The pronounced differences in gene expression between brain parts and the conserved patterns between closely related species with similar morphologies and behavior suggest that unraveling the interactions between genes and behavior will benefit from neurogenomic profiling of distinct brain regions.

Acute Drug Effects on the Human Placental Tissue: The Development of a Placental Murine Xenograft Model.

OBJECTIVE: A pilot study was conducted to establish a human placental xenograft, which could serve as a model to evaluate the effect of toxic exposures during pregnancy. STUDY DESIGN: The protocol consisted of engraftment of third-trimester human placental tissue in immunocompromised mice, after induction of a pseudo-pregnancy state by ovariectomy and progesterone supplementation. To validate the model, the placental tissue before and after engraftment was examined by immunohistochemistry, fluorescence-activated cell sorting (FACS), single-nucleotide polymorphism (SNP) genotyping, and whole transcriptome sequencing (WTSS). The human chorion gonadotropin (hCG) production in serum and urine was examined by enzyme-linked immunosorbent assay. RESULTS: Microscopic evaluation of the placental tissue before and after engraftment revealed a stable morphology and preserved histological structure of the human tissue. Viable trophoblast was present after engraftment and remained stable over time. Vascularization and hormonal secretion (hCG) were present till 3 weeks after engraftment. Thirty-one SNPs were equally present, and there was a stable expression level for 56 451 genes evaluated by whole transcriptome sequencing. CONCLUSION: Although this human placental xenograft model cannot copy the unique uterine environment in which the placenta develops and interacts between the mother and the fetus, it could be a suitable tool to evaluate the acute impact and adaptive processes of the placental tissue to environmental changes.

Spatial dynamics and mixing of bluefin tuna in the Atlantic Ocean and Mediterranean Sea revealed using next-generation sequencing.

The Atlantic bluefin tuna is a highly migratory species emblematic of the challenges associated with shared fisheries management. In an effort to resolve the species' stock dynamics, a genomewide search for spatially informative single nucleotide polymorphisms (SNPs) was undertaken, by way of sequencing reduced representation libraries. An allele frequency approach to SNP discovery was used, combining the data of 555 larvae and young-of-the-year (LYOY) into pools representing major geographical areas and mapping against a newly assembled genomic reference. From a set of 184,895 candidate loci, 384 were selected for validation using 167 LYOY. A highly discriminatory genotyping panel of 95 SNPs was ultimately developed by selecting loci with the most pronounced differences between western Atlantic and Mediterranean Sea LYOY. The panel was evaluated by genotyping a different set of LYOY (n = 326), and from these, 77.8% and 82.1% were correctly assigned to western Atlantic and Mediterranean Sea origins, respectively. The panel revealed temporally persistent differentiation among LYOY from the western Atlantic and Mediterranean Sea (FST  = 0.008, p = .034). The composition of six mixed feeding aggregations in the Atlantic Ocean and Mediterranean Sea was characterized using genotypes from medium (n = 184) and large (n = 48) adults, applying population assignment and mixture analyses. The results provide evidence of persistent population structuring across broad geographic areas and extensive mixing in the Atlantic Ocean, particularly in the mid-Atlantic Bight and Gulf of St. Lawrence. The genomic reference and genotyping tools presented here constitute novel resources useful for future research and conservation efforts.

Antibiotic Resistance Plasmids Cointegrated into a Megaplasmid Harboring the blaOXA-427 Carbapenemase Gene.

OXA-427 is a new class D carbapenemase encountered in different species of Enterobacteriaceae in a Belgian hospital. To study the dispersal of this gene, we performed a comparative analysis of two plasmids containing the blaOXA-427 gene, isolated from a Klebsiella pneumoniae strain and an Enterobacter cloacae complex strain. The two IncA/C2 plasmids containing blaOXA-427 share the same backbone; in the K. pneumoniae strain, however, this plasmid is cointegrated into an IncFIb plasmid, forming a 321-kb megaplasmid with multiple multiresistance regions.

Distinct Effects of T-705 (Favipiravir) and Ribavirin on Influenza Virus Replication and Viral RNA Synthesis.

T-705 (favipiravir) is a new antiviral agent in advanced clinical development for influenza therapy. It is supposed to act as an alternative substrate for the viral polymerase, causing inhibition of viral RNA synthesis or virus mutagenesis. These mechanisms were also proposed for ribavirin, an established and broad antiviral drug that shares structural similarity with T-705. We here performed a comparative analysis of the effects of T-705 and ribavirin on influenza virus and host cell functions. Influenza virus-infected cell cultures were exposed to T-705 or ribavirin during single or serial virus passaging. The effects on viral RNA synthesis and infectious virus yield were determined and mutations appearing in the viral genome were detected by whole-genome virus sequencing. In addition, the cellular nucleotide pools as well as direct inhibition of the viral polymerase enzyme were quantified. We demonstrate that the anti-influenza virus effect of ribavirin is based on IMP dehydrogenase inhibition, which results in fast and profound GTP depletion and an imbalance in the nucleotide pools. In contrast, T-705 acts as a potent and GTP-competitive inhibitor of the viral polymerase. In infected cells, viral RNA synthesis is completely inhibited by T-705 or ribavirin at ≥50 µM, whereas exposure to lower drug concentrations induces formation of noninfectious particles and accumulation of random point mutations in the viral genome. This mutagenic effect is 2-fold higher for T-705 than for ribavirin. Hence, T-705 and ribavirin both act as purine pseudobases but profoundly differ with regard to the mechanism behind their antiviral and mutagenic effects on influenza virus.

A catalog of hemizygous variation in 127 22q11 deletion patients.

The 22q11.2 deletion syndrome is the most common microdeletion disorder, with wide phenotypic variability. To investigate variation within the non-deleted allele we performed targeted resequencing of the 22q11.2 region for 127 patients, identifying multiple deletion sizes, including two deletions with atypical breakpoints. We cataloged ~12,000 hemizygous variant positions, of which 84% were previously annotated. Within the coding regions 95 non-synonymous variants, three stop gains, and two frameshift insertions were identified, some of which we speculate could contribute to atypical phenotypes. We also catalog tolerability of 22q11 gene mutations based on related autosomal recessive disorders in man, embryonic lethality in mice, cross-species conservation and observations that some genes harbor more or less variants than expected. This extensive catalog of hemizygous variants will serve as a blueprint for future experiments to correlate 22q11DS variation with phenotype.

Polymerase specific error rates and profiles identified by single molecule sequencing.

DNA polymerases have an innate error rate which is polymerase and DNA context specific. Historically the mutational rate and profiles have been measured using a variety of methods, each with their own technical limitations. Here we used the unique properties of single molecule sequencing to evaluate the mutational rate and profiles of six DNA polymerases at the sequence level. In addition to accurately determining mutations in double strands, single molecule sequencing also captures direction specific transversions and transitions through the analysis of heteroduplexes. Not only did the error rates vary, but also the direction specific transitions differed among polymerases.

GBSX: a toolkit for experimental design and demultiplexing genotyping by sequencing experiments.

BACKGROUND: Massive parallel sequencing is a powerful tool for variant discovery and genotyping. To reduce costs, sequencing of restriction enzyme based reduced representation libraries can be utilized. This technology is generally referred to as Genotyping By Sequencing (GBS). To deal with GBS experimental design and initial processing specific bioinformatic tools are needed. RESULTS: GBSX is a package that assists in selecting the appropriate enzyme and the design of compatible in-line barcodes. Post sequencing, it performs optimized demultiplexing using these barcodes to create fastq files per barcode which can easily be plugged into existing variant analysis pipelines. Here we demonstrate the usability of the GBSX toolkit and demonstrate improved in-line barcode demultiplexing and trimming performance compared to existing tools. CONCLUSIONS: GBSX provides an easy to use suite of tools for designing and demultiplexing of GBS experiments.

Homozygous missense mutation in STYXL1 associated with moderate intellectual disability, epilepsy and behavioural complexities.

The introduction of massive parallel sequencing has led to the identification of multiple novel genes for intellectual disability (ID) as well as epilepsy. Whereas dominant de novo mutations have been proven to be a leading cause for these disorders, they do not apply to families suggestive of an autosomal recessive inheritance pattern. In this study, we combined the use of linkage analysis with exome sequencing to elucidate the cause of moderate non-syndromic ID, epilepsy and behavioural problems in a consanguineous Asian family. A founder missense mutation was identified in STYXL1. We propose this as a novel candidate gene involved in ID, accompanied by seizures and behavioural problems. Our findings further confirm the genetic heterogeneity of cognitive disorders and genetic epilepsy.

Noninvasive prenatal testing using a novel analysis pipeline to screen for all autosomal fetal aneuploidies improves pregnancy management.

Noninvasive prenatal testing by massive parallel sequencing of maternal plasma DNA has rapidly been adopted as a mainstream method for detection of fetal trisomy 21, 18 and 13. Despite the relative high accuracy of current NIPT testing, a substantial number of false-positive and false-negative test results remain. Here, we present an analysis pipeline, which addresses some of the technical as well as the biologically derived causes of error. Most importantly, it differentiates high z-scores due to fetal trisomies from those due to local maternal CNVs causing false positives. This pipeline was retrospectively validated for trisomy 18 and 21 detection on 296 samples demonstrating a sensitivity and specificity of 100%, and applied prospectively to 1350 pregnant women in the clinical diagnostic setting with a result reported in 99.9% of cases. In addition, values indicative for trisomy were observed two times for chromosome 7 and once each for chromosomes 15 and 16, and once for a segmental trisomy 18. Two of the trisomies were confirmed to be mosaic, one of which contained a uniparental disomy cell line. As placental trisomies pose a risk for low-grade fetal mosaicism as well as uniparental disomy, genome-wide noninvasive aneuploidy detection is improving prenatal management.

Pseudoautosomal region 1 length polymorphism in the human population.

The human sex chromosomes differ in sequence, except for the pseudoautosomal regions (PAR) at the terminus of the short and the long arms, denoted as PAR1 and PAR2. The boundary between PAR1 and the unique X and Y sequences was established during the divergence of the great apes. During a copy number variation screen, we noted a paternally inherited chromosome X duplication in 15 independent families. Subsequent genomic analysis demonstrated that an insertional translocation of X chromosomal sequence into the Y chromosome generates an extended PAR [corrected].The insertion is generated by non-allelic homologous recombination between a 548 bp LTR6B repeat within the Y chromosome PAR1 and a second LTR6B repeat located 105 kb from the PAR boundary on the X chromosome. The identification of the reciprocal deletion on the X chromosome in one family and the occurrence of the variant in different chromosome Y haplogroups demonstrate this is a recurrent genomic rearrangement in the human population. This finding represents a novel mechanism shaping sex chromosomal evolution.

NGS-Logistics: federated analysis of NGS sequence variants across multiple locations.

As many personal genomes are being sequenced, collaborative analysis of those genomes has become essential. However, analysis of personal genomic data raises important privacy and confidentiality issues. We propose a methodology for federated analysis of sequence variants from personal genomes. Specific base-pair positions and/or regions are queried for samples to which the user has access but also for the whole population. The statistics results do not breach data confidentiality but allow further exploration of the data; researchers can negotiate access to relevant samples through pseudonymous identifiers. This approach minimizes the impact on data confidentiality while enabling powerful data analysis by gaining access to important rare samples. Our methodology is implemented in an open source tool called NGS-Logistics, freely available at https://ngsl.esat.kuleuven.be.

Microsatellite primers for the gynodioecious grassland perennial Saxifraga granulata (Saxifragaceae).

PREMISE OF THE STUDY: Nine polymorphic and 12 monomorphic microsatellite loci (simple sequence repeats [SSRs]) were isolated and characterized for the gynodioecious grassland perennial Saxifraga granulata. • METHODS AND RESULTS: Based on genomic screening of leaf material of four individuals from four populations, a total of 21 microsatellite primer pairs were designed for S. granulata. Nine loci were polymorphic and were optimized into two PCR multiplex reactions and tested on 100 individuals from five riparian populations from central Belgium. The number of alleles of the polymorphic loci ranged from three to 18, and gametic heterozygosity ranged from 0.26 to 0.94. • CONCLUSIONS: The markers that are presented here are the first microsatellite markers reported for S. granulata and will be used to assess how river systems shape the spatial genetic structure and diversity of riparian populations of this species.