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.

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.

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.

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.

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.

Heritability of cortisol response to confinement stress in European sea bass dicentrarchus labrax.

BACKGROUND: In fish, the most studied production traits in terms of heritability are body weight or growth, stress or disease resistance, while heritability of cortisol levels, widely used as a measure of response to stress, is less studied. In this study, we have estimated heritabilities of two growth traits (body weight and length) and of cortisol response to confinement stress in the European sea bass. FINDINGS: The F1 progeny analysed (n = 922) belonged to a small effective breeding population with contributions from an unbalanced family structure of just 10 males and 2 females. Heritability values ranged from 0.54 (± 0.21) for body weight to 0.65 (± 0.22) for standard body length and were low for cortisol response i.e. 0.08 (± 0.06). Genetic correlations were positive (0.94) between standard body length and body weight and negative between cortisol and body weight and between cortisol and standard body length (-0.60 and -0.55, respectively). CONCLUSION: This study confirms that in European sea bass, heritability of growth-related traits is high and that selection on such traits has potential. However, heritability of cortisol response to stress is low in European sea bass and since it is known to vary greatly among species, further studies are necessary to understand the reasons for these differences.

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.

Multilocus analyses of an Antarctic fish species flock (Teleostei, Notothenioidei, Trematominae): phylogenetic approach and test of the early-radiation event.

Clades that have undergone episodes of rapid cladogenesis are challenging from a phylogenetic point of view. They are generally characterised by short or missing internal branches in phylogenetic trees and by conflicting topologies among individual gene trees. This may be the case of the subfamily Trematominae, a group of marine teleosts of coastal Antarctic waters, which is considered to have passed through a period of rapid diversification. Despite much phylogenetic attention, the relationships among Trematominae species remain unclear. In contrast to previous studies that were mostly based on concatenated datasets of mitochondrial and/or single nuclear loci, we applied various single-locus and multilocus phylogenetic approaches to sequences from 11 loci (eight nuclear) and we also used several methods to assess the hypothesis of a radiation event in Trematominae evolution. Diversification rate analyses support the hypothesis of a period of rapid diversification during Trematominae history and only a few nodes in the hypothetical species tree were consistently resolved with various phylogenetic methods. We detected significant discrepancies among trees from individual genes of these species, most probably resulting from incomplete lineage sorting, suggesting that concatenation of loci is not the most appropriate way to investigate Trematominae species interrelationships. These data also provide information about the possible effects of historic climate changes on the diversification rate of this group of fish.

Telomere healing following DNA polymerase arrest-induced breakages is likely the main mechanism generating chromosome 4p terminal deletions.

Constitutional developmental disorders are frequently caused by terminal chromosomal deletions. The mechanisms and/or architectural features that might underlie those chromosome breakages remain largely unexplored. Because telomeres are the vital DNA protein complexes stabilizing linear chromosomes against chromosome degradation, fusion, and incomplete replication, those terminal-deleted chromosomes acquired new telomeres either by telomere healing or by telomere capture. To unravel the mechanisms leading to chromosomal breakage and healing, we sequenced nine chromosome 4p terminal deletion boundaries. A computational analysis of the breakpoint flanking region, including 12 previously published pure terminal breakage sites, was performed in order to identify architectural features that might be involved in this process. All terminal 4p truncations were likely stabilized by telomerase-mediated telomere healing. In the majority of breakpoints multiple genetic elements have a potential to induce secondary structures and an enrichment in replication stalling site motifs were identified. These findings suggest DNA replication stalling-induced chromosome breakage during early development is the first mechanistic step leading toward terminal deletion syndromes.

Differential modes of selection on the rhodopsin gene in coastal Baltic and North Sea populations of the sand goby, Pomatoschistus minutus.

An excellent model to elucidate the mechanisms and importance of evolution in the marine environment is the spectral tuning mechanism of the visual pigment in vertebrates. In the sand goby Pomatoschistus minutus (Teleostei; Gobiidae), a distribution-wide study showed that spatial variation at the rhodopsin gene (RH1) matches the characteristics of specific light environments. This match suggests that populations are locally adapted to selective light regimes targeting the RH1 gene. If so, then the direction of selection should depend on the regional spatial and temporal stability of the light conditions. We tested this prediction by comparing goby populations from two regions: the Baltic Sea, characterized by divergent, but temporally stable light conditions, and the North Sea, characterized by locally heterogeneous and temporally variable light conditions. RH1 sequences of 491 Pomatoschistus minutus individuals from 15 locations were analysed. We found that variation at the RH1 gene in the Baltic populations showed signatures of diversifying selection, whereas the RH1 gene in the North Sea showed signatures of stabilizing selection. These different modes of selection are consistent with the regional light conditions and hence support our predictions, but may also be influenced by migration between the open sea and more turbid estuarine environments. An interesting observation is that within one gene, synonymous and non-synonymous SNPs show a totally different pattern between populations. Population differentiation based on non-synonymous SNPs of the RH1 gene correlated with spectral variation of the local environment of the sand goby populations. In contrast, the differentiation based on synonymous SNPs of RH1 reflects more the neutral historical pattern of the species.

High molecular diversity in the rhodopsin gene in closely related goby fishes: A role for visual pigments in adaptive speciation?

The spectral tuning mechanism of visual pigments is an excellent model to elucidate the mechanisms of adaptive evolution and the importance of selection as an evolutionary force. Therefore, we use a phylogenetic approach to determine whether there is evidence for differential adaptive molecular evolution on the rhodopsin (RH1) gene among closely related 'sand goby' species (Teleostei, Gobiidae). Fragments of the RH1 gene (868 bp) were sequenced and analyzed for nine 'sand goby' species that inhabit different photic environments. A high level of interspecific polymorphism at the RH1 gene was observed, including non-synonymous mutations on amino acids known as spectral tuning sites. Clear indications for positive Darwinian selection were provided by three independent methods: (1) by linking functional variation on the RH1 gene to specific light environments of the different fish habitats; (2) by constructing and comparing phylogenies based on RH1 and the 'neutral' 12S and 16S mtDNA fragments; and (3) by performing statistical tests to detect signatures of directional selection on the RH1 gene. This study shows an unusual high variability in the gobiid visual RH1 pigment, and we therefore suggest a possible role for sensory genes in the adaptive radiation of 'sand goby' species.

To see in different seas: spatial variation in the rhodopsin gene of the sand goby (Pomatoschistus minutus).

Aquatic organisms living in a range of photic environments require specific mechanisms to tune their visual pigments. Maximum absorbance (lambda(max)) of retinal rods in populations of the marine demersal sand goby, (Pomatoschistus minutus; Gobiidae, Teleostei) correlates with the local optic environment. It has been shown that this is not regulated through a physiological response by exchanging the rhodopsin chromophore. To test for evolutionary adaptation, the sequence of the rhodopsin (RH1) gene was analysed in 165 Pomatoschistus minutus individuals from seven populations across its distribution range. Analysis showed a high level of intraspecific polymorphism at the RH1 gene, including nonsynonymous mutations on amino acids, known as spectral tuning sites. Population differentiation at these sites was in agreement with the observed differentiation in lambda(max) values. Analyses of d(N)/d(S) substitution rate ratios and likelihood ratio tests under site-specific models detected a significant signal of positive Darwinian selection on the RH1 gene. A strong discrepancy in differentiation was noticed between RH1 gene variation and the presumably neutral microsatellites and mitochondrial data. Samples did not cluster according to geographical or historical proximity with regards to RH1, but according to the general photic conditions of the habitat environment of the sand goby. This study highlights the usefulness of sensory genes, like rhodopsin, for studying the characteristics of local adaptation in marine nonmodel organisms.

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.

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.

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.

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.

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.

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.

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.

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.