Pre-existing cross-reactive neutralizing activity against SARS-CoV-2 and seasonal coronaviruses prior to the COVID-19 pandemic (2014-2019) with limited immunity against recent emerging SARS-CoV-2 variants, Vietnam.

OBJECTIVES: SARS-CoV-2 transmission and epidemic potential is related to the population's immunity levels. As such, assessing different regions' preexisting immune responses to SARS-CoV-2 is important to understand the transmission potential of emerging SARS-CoV-2 variants. DESIGN: In 975 serum samples from Vietnam (2014 to 2019), anti-SARS-CoV-2 Immunoglobulin G levels were determined by enzyme-linked immunosorbent assay. Plaque reduction neutralization test (PRNT) was performed using Wuhan strain and variants of concern (VOCs). Cross-reactivity was confirmed by analyzing B-cell receptor (BCR) repertoire sequences and identifying BCR repertoire sequences-derived T-cell epitopes. RESULTS: Overall, 20.9% (n = 76/364) and 9.2% (n = 7) demonstrated SARS-CoV-2 neutralizing activity (PRNT50) against the Wuhan and Alpha strain, respectively. Neutralizing activity against Beta, Gamma, and Delta strains was absent (PRNT50<5) in all samples. Cross-reactive epitopes against SARS-CoV-2 and other coronavirus spike proteins were detected in the N-terminal domain, S2, and receptor-binding domain regions. CONCLUSIONS: Following BCR and major histocompatibility complex analysis, T-cell receptor-recognized epitope motif (TREM) among pathogenic coronaviruses and coronaviruses spike proteins were the top TREM peptide, suggesting that pre-existing immunity against SARS-CoV-2 in Vietnam was due to exposure to common cold coronaviruses. With limited immunity against emerging VOCs, further monitoring, and control of the epidemic, along with COVID-19 vaccine programs against VOCs, are necessary.

Idiotope-Driven T-Cell/B-Cell Collaboration-Based T-Cell Epitope Prediction Using B-Cell Receptor Repertoire Sequences in Infectious Diseases.

T-cell recognition of antigen epitopes is a crucial step for the induction of adaptive immune responses, and the identification of such T-cell epitopes is, therefore, important for understanding diverse immune responses and controlling T-cell immunity. A number of bioinformatic tools exist that predict T-cell epitopes; however, many of these methods highly rely on evaluating conventional peptide presentation by major histocompatibility complex (MHC) molecules, but they ignore epitope sequences recognized by T-cell receptor (TCR). Immunogenic determinant idiotopes are present on the variable regions of immunoglobulin molecules expressed on and secreted by B-cells. In idiotope-driven T-cell/B-cell collaboration, B-cells present the idiotopes on MHC molecules for recognition by idiotope-specific T-cells. According to the idiotype network theory formulated by Niels Jerne, such idiotopes found on anti-idiotypic antibodies exhibit molecular mimicry of antigens. Here, by combining these concepts and defining the patterns of TCR-recognized epitope motifs (TREMs), we developed a T-cell epitope prediction method that identifies T-cell epitopes derived from antigen proteins by analyzing B-cell receptor (BCR) sequences. This method allowed us to identify T-cell epitopes that contain the same TREM patterns between BCR and viral antigen sequences in two different infectious diseases caused by dengue virus and SARS-CoV-2 infection. The identified epitopes were among the T-cell epitopes detected in previous studies, and T-cell stimulatory immunogenicity was confirmed. Thus, our data support this method as a powerful tool for the discovery of T-cell epitopes from BCR sequences.

Pharmacological MEK inhibition promotes polyclonal T-cell reconstitution and suppresses xenogeneic GVHD.

Rapid immune reconstitution without developing graft-versus-host disease (GVHD) is required for the success of allogeneic hematopoietic stem cell transplantation. Here, we analyzed the effects of pharmacological MEK inhibition on human polyclonal T-cell reconstitution in a humanized mouse GVHD model utilizing deep sequencing-based T-cell receptor (TCR) repertoire analysis. GVHD mice exhibited a skewed TCR repertoire with a common clone within target organs. The MEK inhibitor trametinib ameliorated GVHD and enabled engraftment of diverse T-cell clones. Furthermore, trametinib also ameliorated GVHD sparing diverse T cell repertoire, even when it was given from day 15 through 28. Although tacrolimus also reduced development of GVHD, it disturbed diverse T cell reconstitution and resulted in skewed TCR repertoire. Thus, trametinib not only suppresses GVHD-inducing T cells but also promotes human T cell reconstitution in vivo, providing a novel rationale for translational studies targeting human GVHD.

Construction of complete Tupaia belangeri transcriptome database by whole-genome and comprehensive RNA sequencing.

The northern tree shrew (Tupaia belangeri) possesses high potential as an animal model of human diseases and biology, given its genetic similarity to primates. Although genetic information on the tree shrew has already been published, some of the entire coding sequences (CDSs) of tree shrew genes remained incomplete, and the reliability of these CDSs remained difficult to determine. To improve the determination of tree shrew CDSs, we performed sequencing of the whole-genome, mRNA, and total RNA and integrated the resulting data. Additionally, we established criteria for the selection of reliable CDSs and annotated these sequences by comparison to the human transcriptome, resulting in the identification of complete CDSs for 12,612 tree shrew genes and yielding a more accurate tree shrew genome database (TupaiaBase: http://tupaiabase.org ). Transcriptome profiles in hepatitis B virus infected tree shrew livers were analyzed for validation. Gene ontology analysis showed enriched transcriptional regulation at 1 day post-infection, namely in the "type I interferon signaling pathway". Moreover, a negative regulator of type I interferon, SOCS3, was induced. This work, which provides a tree shrew CDS database based on genomic DNA and RNA sequencing, is expected to serve as a powerful tool for further development of the tree shrew model.

Epstein-Barr virus-related diffuse large B-cell lymphoma in mogamulizumab-treated adult T-cell leukemia with incomplete T-cell reconstitution.

Adult T-cell leukemia (ATL) is an aggressive mature T-cell malignancy with a poor prognosis. The anti-C-C motif chemokine receptor 4 (CCR4) antibody mogamulizumab (moga) reduces ATL cells and induces reconstitution of polyclonal T cells; however, ATL cases often remain resistant and moga sometimes causes fatal immunopathology. Epstein-Barr virus (EBV)-related B-cell lymphoma develops in severely immunocompromised subjects, and is particularly associated with impaired T-cell immunity. Here, we report an ATL patient who had received conventional chemotherapy plus moga, and subsequently developed EBV-related diffuse large B-cell lymphoma (DLBCL) of the central nervous system. Next-generation sequencing-based T-cell receptor repertoire analyses identified residual abnormal clones and revealed that reconstitution of polyclonal T cells was incomplete, even after moga treatment. Furthermore, a skin rash that developed after moga treatment was found to contain ATL clones. This case suggests that the limited therapeutic effects of moga and incomplete T-cell reconstitution are associated with severely impaired T-cell immunity and subsequent development of EBV-related DLBCL.

Next-Generation Immune Repertoire Sequencing as a Clue to Elucidate the Landscape of Immune Modulation by Host-Gut Microbiome Interactions.

The human immune system is a fine network consisted of the innumerable numbers of functional cells that balance the immunity and tolerance against various endogenous and environmental challenges. Although advances in modern immunology have revealed a role of many unique immune cell subsets, technologies that enable us to capture the whole landscape of immune responses against specific antigens have been not available to date. Acquired immunity against various microorganisms including host microbiome is principally founded on T cell and B cell populations, each of which expresses antigen-specific receptors that define a unique clonotype. Over the past several years, high-throughput next-generation sequencing has been developed as a powerful tool to profile T- and B-cell receptor repertoires in a given individual at the single-cell level. Sophisticated immuno-bioinformatic analyses by use of this innovative methodology have been already implemented in clinical development of antibody engineering, vaccine design, and cellular immunotherapy. In this article, we aim to discuss the possible application of high-throughput immune receptor sequencing in the field of nutritional and intestinal immunology. Although there are still unsolved caveats, this emerging technology combined with single-cell transcriptomics/proteomics provides a critical tool to unveil the previously unrecognized principle of host-microbiome immune homeostasis. Accumulation of such knowledge will lead to the development of effective ways for personalized immune modulation through deeper understanding of the mechanisms by which the intestinal environment affects our immune ecosystem.

Deep sequencing of the T cell receptor visualizes reconstitution of T cell immunity in mogamulizumab-treated adult T cell leukemia.

Although the anti-CCR4 antibody mogamulizumab (moga) shows striking antitumor activity against adult T cell leukemia (ATL), it can also cause fatal immunological pathology such as severe skin rash and graft-versus-host disease, which might be attributed to depletion of CCR4+ regulatory T cells. We previously showed that next generation sequencing enables precise analysis of the T cell receptor (TCR) repertoire, and we here used the technique to reveal the immunological dynamics in moga-treated ATL patients. Treatment with moga resulted in remarkable reduction or elimination of clonal cells, and enhanced reconstitution of non-tumor polyclonal CD4+ T cells and oligoclonal CD8+ T cells. Interestingly, cutaneous T cells infiltrating moga-related skin rashes did not share the same major clones in peripheral blood, which minimizes the possibility of cross-reaction. Thus, deep sequencing of the TCR can reveal the immune reconstitution of moga-treated ATL and provides powerful insights into its mode of action.

Highly functional T-cell receptor repertoires are abundant in stem memory T cells and highly shared among individuals.

To expand our knowledge of the ontogeny of the T-cell receptor (TCR) repertoire of antigen-specific T-cell subsets, we combined next-generation deep sequencing and single-cell multiplex clonotype analysis to evaluate the diversity and frequency of paired TCRs, their functions and whether clonotypic TCRs are shared among different individuals. Using an HLA-A*02-restricted cytomegalovirus (CMV) pp65-derived immunogenic peptide, we found that the more dominant pp65-specific TCR clonotypes in the blood of healthy donors have higher binding affinities for the CMV peptide and arise from clonotypes that are highly shared among individuals. Interestingly, these highly shared HLA-A*02-restricted CMV-specific TCRs were detected in a CMV-seronegative individual as well as in HLA-A*02-negative donors albeit at lower frequency. More intriguingly, these shared TCR clonotypes were abundant in the stem memory T-cell subset, and TCR diversity of the stem memory T-cell repertoire was significantly lower than in the central memory and effector memory T-cell repertoires. These results suggest that the stem memory T-cell subset may serve as a reservoir of highly shared and highly functional memory T-cells.

Genome sequence and analysis of the Japanese morning glory Ipomoea nil.

Ipomoea is the largest genus in the family Convolvulaceae. Ipomoea nil (Japanese morning glory) has been utilized as a model plant to study the genetic basis of floricultural traits, with over 1,500 mutant lines. In the present study, we have utilized second- and third-generation-sequencing platforms, and have reported a draft genome of I. nil with a scaffold N50 of 2.88 Mb (contig N50 of 1.87 Mb), covering 98% of the 750 Mb genome. Scaffolds covering 91.42% of the assembly are anchored to 15 pseudo-chromosomes. The draft genome has enabled the identification and cataloguing of the Tpn1 family transposons, known as the major mutagen of I. nil, and analysing the dwarf gene, CONTRACTED, located on the genetic map published in 1956. Comparative genomics has suggested that a whole genome duplication in Convolvulaceae, distinct from the recent Solanaceae event, has occurred after the divergence of the two sister families.

Extremotolerant tardigrade genome and improved radiotolerance of human cultured cells by tardigrade-unique protein.

Tardigrades, also known as water bears, are small aquatic animals. Some tardigrade species tolerate almost complete dehydration and exhibit extraordinary tolerance to various physical extremes in the dehydrated state. Here we determine a high-quality genome sequence of Ramazzottius varieornatus, one of the most stress-tolerant tardigrade species. Precise gene repertoire analyses reveal the presence of a small proportion (1.2% or less) of putative foreign genes, loss of gene pathways that promote stress damage, expansion of gene families related to ameliorating damage, and evolution and high expression of novel tardigrade-unique proteins. Minor changes in the gene expression profiles during dehydration and rehydration suggest constitutive expression of tolerance-related genes. Using human cultured cells, we demonstrate that a tardigrade-unique DNA-associating protein suppresses X-ray-induced DNA damage by ∼40% and improves radiotolerance. These findings indicate the relevance of tardigrade-unique proteins to tolerability and tardigrades could be a bountiful source of new protection genes and mechanisms.

Effector Regulatory T Cells Reflect the Equilibrium between Antitumor Immunity and Autoimmunity in Adult T-cell Leukemia.

The regulatory T cells (Treg) with the most potent immunosuppressive activity are the effector Tregs (eTreg) with a CD45RA(-)Foxp3(++)CCR4(+) phenotype. Adult T-cell leukemia (ATL) cells often share the Treg phenotype and also express CCR4. Although mogamulizumab, a monoclonal antibody to CCR4, shows marked antitumor effects against ATL and peripheral T-cell lymphoma, concerns have been raised that it may induce severe autoimmune immunopathology by depleting eTregs. Here, we present case reports for two patients with ATL who responded to mogamulizumab but developed a severe skin rash and autoimmune brainstem encephalitis. Deep sequencing of the T-cell receptor revealed that ATL cells and naturally occurring Tregs within the cell population with a Treg phenotype can be clearly distinguished according to CADM1 expression. The onset of skin rash and brainstem encephalitis was coincident with eTreg depletion from the peripheral blood, whereas ATL relapses were coincident with eTreg recovery. These results imply that eTreg numbers in the peripheral blood sensitively reflect the equilibrium between antitumor immunity and autoimmunity, and that mogamulizumab might suppress ATL until the eTreg population recovers. Close monitoring of eTreg numbers is crucial if we are to provide immunomodulatory treatments that target malignancy without severe adverse events. Cancer Immunol Res; 4(8); 644-9. ©2016 AACR.

Comparative genome and transcriptome analyses of the social amoeba Acytostelium subglobosum that accomplishes multicellular development without germ-soma differentiation.

BACKGROUND: Social amoebae are lower eukaryotes that inhabit the soil. They are characterized by the construction of a starvation-induced multicellular fruiting body with a spore ball and supportive stalk. In most species, the stalk is filled with motile stalk cells, as represented by the model organism Dictyostelium discoideum, whose developmental mechanisms have been well characterized. However, in the genus Acytostelium, the stalk is acellular and all aggregated cells become spores. Phylogenetic analyses have shown that it is not an ancestral genus but has lost the ability to undergo cell differentiation. RESULTS: We performed genome and transcriptome analyses of Acytostelium subglobosum and compared our findings to other available dictyostelid genome data. Although A. subglobosum adopts a qualitatively different developmental program from other dictyostelids, its gene repertoire was largely conserved. Yet, families of polyketide synthase and extracellular matrix proteins have not expanded and a serine protease and ABC transporter B family gene, tagA, and a few other developmental genes are missing in the A. subglobosum lineage. Temporal gene expression patterns are astonishingly dissimilar from those of D. discoideum, and only a limited fraction of the ortholog pairs shared the same expression patterns, so that some signaling cascades for development seem to be disabled in A. subglobosum. CONCLUSIONS: The absence of the ability to undergo cell differentiation in Acytostelium is accompanied by a small change in coding potential and extensive alterations in gene expression patterns.

A cDNA resource for gene expression studies of a hemichordate, Ptychodera flava.

Recent investigations into the evolution of deuterostomes and the origin of chordates have paid considerable attention to hemichordates (acorn worms), as hemichordates and echinoderms are the closest chordate relatives. The present study prepared cDNA libraries from Ptychodera flava, to study expression and function of genes involved in development of the hemichordate body plan. Expressed sequence tag (EST) analyses of nine cDNA libraries yielded 18,832 cloned genes expressed in eggs, 18,739 in blastulae, 18,539 in gastrulae, 18,811 in larvae, 18,978 in juveniles, 11,802 in adult proboscis, 17,259 in stomochord, 11,886 in gills, and 11,580 in liver, respectively. A set of 34,159 uni-gene clones of P. flava was obtained. This cDNA resource will be valuable for studying temporal and spatial expression of acorn worm genes during development.

Co-option of Sox3 as the male-determining factor on the Y chromosome in the fish Oryzias dancena.

Sex chromosomes harbour a primary sex-determining signal that triggers sexual development of the organism. However, diverse sex chromosome systems have been evolved in vertebrates. Here we use positional cloning to identify the sex-determining locus of a medaka-related fish, Oryzias dancena, and find that the locus on the Y chromosome contains a cis-regulatory element that upregulates neighbouring Sox3 expression in developing gonad. Sex-reversed phenotypes in Sox3(Y) transgenic fish, and Sox3(Y) loss-of-function mutants all point to its critical role in sex determination. Furthermore, we demonstrate that Sox3 initiates testicular differentiation by upregulating expression of downstream Gsdf, which is highly conserved in fish sex differentiation pathways. Our results not only provide strong evidence for the independent recruitment of Sox3 to male determination in distantly related vertebrates, but also provide direct evidence that a novel sex determination pathway has evolved through co-option of a transcriptional regulator potentially interacted with a conserved downstream component.

The ancestor of extant Japanese fancy mice contributed to the mosaic genomes of classical inbred strains.

Commonly used classical inbred mouse strains have mosaic genomes with sequences from different subspecific origins. Their genomes are derived predominantly from the Western European subspecies Mus musculus domesticus, with the remaining sequences derived mostly from the Japanese subspecies Mus musculus molossinus. However, it remains unknown how this intersubspecific genome introgression occurred during the establishment of classical inbred strains. In this study, we resequenced the genomes of two M. m. molossinus-derived inbred strains, MSM/Ms and JF1/Ms. MSM/Ms originated from Japanese wild mice, and the ancestry of JF1/Ms was originally found in Europe and then transferred to Japan. We compared the characteristics of these sequences to those of the C57BL/6J reference sequence and the recent data sets from the resequencing of 17 inbred strains in the Mouse Genome Project (MGP), and the results unequivocally show that genome introgression from M. m. molossinus into M. m. domesticus provided the primary framework for the mosaic genomes of classical inbred strains. Furthermore, the genomes of C57BL/6J and other classical inbred strains have long consecutive segments with extremely high similarity (>99.998%) to the JF1/Ms strain. In the early 20th century, Japanese waltzing mice with a morphological phenotype resembling that of JF1/Ms mice were often crossed with European fancy mice for early studies of "Mendelism," which suggests that the ancestor of the extant JF1/Ms strain provided the origin of the M. m. molossinus genome in classical inbred strains and largely contributed to its intersubspecific genome diversity.

A genomic overview of short genetic variations in a basal chordate, Ciona intestinalis.

BACKGROUND: Although the Ciona intestinalis genome contains many allelic polymorphisms, there is only limited data analyzed systematically. Establishing a dense map of genetic variations in C. intestinalis is necessary not only for linkage analysis, but also for other experimental biology including molecular developmental and evolutionary studies, because animals from natural populations are typically used for experiments. RESULTS: Here, we identified over three million candidate short genomic variations within a 110 Mb euchromatin region among five C. intestinalis individuals. The average nucleotide diversity was approximately 1.1%. Genetic variations were found at a similar density in intergenic and gene regions. Non-synonymous and nonsense nucleotide substitutions were found in 12,493 and 1,214 genes accounting for 81.9% and 8.0% of the entire gene set, respectively, and over 60% of genes in the single animal encode non-identical proteins between maternal and paternal alleles. CONCLUSIONS: Our results provide a framework for studying evolution of the animal genome, as well as a useful resource for a wide range of C. intestinalis researchers.

Centromere-targeted de novo integrations of an LTR retrotransposon of Arabidopsis lyrata.

The plant genome evolves with rapid proliferation of LTR-type retrotransposons, which is associated with their clustered accumulation in gene-poor regions, such as centromeres. Despite their major role for plant genome evolution, no mobile LTR element with targeted integration into gene-poor regions has been identified in plants. Here, we report such targeted integrations de novo. We and others have previously shown that an ATCOPIA93 family retrotransposon in Arabidopsis thaliana is mobilized when the DNA methylation machinery is compromised. Although ATCOPIA93 family elements are low copy number in the wild-type A. thaliana genome, high-copy-number related elements are found in the wild-type Arabidopsis lyrata genome, and they show centromere-specific localization. To understand the mechanisms for the clustered accumulation of the A. lyrata elements directly, we introduced one of them, named Tal1 (Transposon of Arabidopsis lyrata 1), into A. thaliana by transformation. The introduced Tal1 was retrotransposed in A. thaliana, and most of the retrotransposed copies were found in centromeric repeats of A. thaliana, suggesting targeted integration. The targeted integration is especially surprising because the centromeric repeat sequences differ considerably between A. lyrata and A. thaliana. Our results revealed unexpectedly dynamic controls for evolution of the transposon-rich heterochromatic regions.

Comprehensive functional analyses of expressed sequence tags in common wheat (Triticum aestivum).

About 1 million expressed sequence tag (EST) sequences comprising 125.3 Mb nucleotides were accreted from 51 cDNA libraries constructed from a variety of tissues and organs under a range of conditions, including abiotic stresses and pathogen challenges in common wheat (Triticum aestivum). Expressed sequence tags were assembled with stringent parameters after processing with inbuild scripts, resulting in 37,138 contigs and 215,199 singlets. In the assembled sequences, 10.6% presented no matches with existing sequences in public databases. Functional characterization of wheat unigenes by gene ontology annotation, mining transcription factors, full-length cDNA, and miRNA targeting sites were carried out. A bioinformatics strategy was developed to discover single-nucleotide polymorphisms (SNPs) within our large EST resource and reported the SNPs between and within (homoeologous) cultivars. Digital gene expression was performed to find the tissue-specific gene expression, and correspondence analysis was executed to identify common and specific gene expression by selecting four biotic stress-related libraries. The assembly and associated information cater a framework for future investigation in functional genomics.

Molecular cloning and characterization of the repetitive DNA sequences that comprise the constitutive heterochromatin of the W chromosomes of medaka fishes.

Among the medaka fishes of the genus Oryzias, most species have homomorphic sex chromosomes, while some species, such as Oryzias hubbsi and Oryzias javanicus, have heteromorphic ZW sex chromosomes. In this study, a novel family of repetitive sequence was molecularly cloned from O. hubbsi and characterized by chromosome in situ and filter hybridization, respectively. This repetitive element, which we designated as a BstNI family element, localized at heterochromatin regions on the W chromosome, as well as on two pairs of autosomes. Homologous sequences to this element were found only in O. javanicus, which is a sister species of O. hubbsi, suggesting that this repeated element originated in the common ancestor of these two species. However, the intensity of the hybridization signals was lower in O. javanicus than in O. hubbsi, and the chromosomal location of this element in O. javanicus was confined to heterochromatin regions on one pair of autosomes. Thus, we hypothesize that this repetitive element was extensively amplified in the O. hubbsi lineage, especially on its W chromosome, after the separation of the O. javanicus lineage. In addition, we also found the W chromosomal location of the 18S-28S ribosomal RNA genes in both O. hubbsi and O. javanicus. Our previous studies showed no linkage homology of the sex chromosomes in these species, indicating that the RNA genes were shared between W chromosomes of different origins. This situation may be explained by a translocation of the sex-determining region with the ribosomal RNA genes in either species or an independent accumulation of the RNA genes as a convergent process during W chromosome degeneration.

Easy-to-use phylogenetic analysis system for hepatitis B virus infection.

AIM: The molecular phylogenetic analysis has been broadly applied to clinical and virological study. However, the appropriate settings and application of calculation parameters are difficult for non-specialists of molecular genetics. In the present study, the phylogenetic analysis tool was developed for the easy determination of genotypes and transmission route. METHODS: A total of 23 patients of 10 families infected with hepatitis B virus (HBV) were enrolled and expected to undergo intrafamilial transmission. The extracted HBV DNA were amplified and sequenced in a region of the S gene. RESULTS: The software to automatically classify query sequence was constructed and installed on the Hepatitis Virus Database (HVDB). Reference sequences were retrieved from HVDB, which contained major genotypes from A to H. Multiple-alignments using CLUSTAL W were performed before the genetic distance matrix was calculated with the six-parameter method. The phylogenetic tree was output by the neighbor-joining method. User interface using WWW-browser was also developed for intuitive control. This system was named as the easy-to-use phylogenetic analysis system (E-PAS). Twenty-three sera of 10 families were analyzed to evaluate E-PAS. The queries obtained from nine families were genotype C and were located in one cluster per family. However, one patient of a family was classified into the cluster different from her family, suggesting that E-PAS detected the sample distinct from that of her family on the transmission route. CONCLUSIONS: The E-PAS to output phylogenetic tree was developed since requisite material was sequence data only. E-PAS could expand to determine HBV genotypes as well as transmission routes.