Search

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has caused the pandemic of coronavirus disease 2019 (COVID-19). Great international efforts have been put into the development of prophylactic vaccines and neutralizing antibodies. However, the knowledge about the B cell immune response induced by the SARS-CoV-2 virus is still limited. Here, we report a comprehensive characterization of the dynamics of immunoglobin heavy chain (IGH) repertoire in COVID-19 patients. By using next-generation sequencing technology, we examined the temporal changes in the landscape of the patient's immunological status and found dramatic changes in the IGH within the patient's immune system after the onset of COVID-19 symptoms. Although different patients have distinct immune responses to SARS-CoV-2 infection, by employing clonotype overlap, lineage expansion, and clonotype network analyses, we observed a higher clonotype overlap and substantial lineage expansion of B cell clones 2-3 weeks after the onset of illness, which is of great importance to B-cell immune responses. Meanwhile, for preferences of V gene usage during SARS-CoV-2 infection, IGHV3-74 and IGHV4-34, and IGHV4-39 in COVID-19 patients were more abundant than those of healthy controls. Overall, we present an immunological resource for SARS-CoV-2 that could promote both therapeutic development as well as mechanistic research.

Currently circulating SARS-CoV-2 Omicron variants feature highly mutated spike proteins with extraordinary abilities in evading acute-infection-induced germline antibodies isolated earlier in the pandemic. We identified that memory B cells from Delta variant breakthrough-infection patients expressed antibodies with more extensive somatic hypermutations (SHMs) allowing isolation of a number of broadly neutralizing antibodies with activities against heterologous variants of concerns (VOCs) including Omicron variant. Structural studies identified that SHM introduced altered amino acids and highly unusual HCDR2 insertions respectively in two representative broadly neutralizing antibodies - YB9-258 and YB13-292. Previously, insertion/deletion were rarely reported for antiviral antibodies except for those induced by HIV-1 chronic infections. Identified SHMs involved heavily in epitope recognition, they broadened neutralization breadth by rendering antibodies resistant to VOC mutations highly detrimental to previously isolated antibodies targeting similar epitopes. These data provide molecular mechanisms for enhanced immunity to heterologous SARS-CoV-2 variants after repeated antigen exposures with implications for future vaccination strategy.

The field of antibody development is under pressure to meet rising demands for speed, cost-effectiveness, efficacy, reliability, and large-scale production. It is costly and time-consuming to immunize animals and build a single-domain antibody (sdAb) library for each target. Using the variable domain (VHH) of heavy-chain only antibodies (HcAbs) derived from blood samples of 75 non-immunized camelid animals (51 alpacas, 13 llamas, 11 Bactrian camels), and spleens from two Bactrian camels, a naïve sdAb library with extensive megadiversity and reusability was constructed. The library was evaluated using next-generation DNA sequencing (NGS) and was found to contain hundreds of billions of unique clones. To confirm the availability of target-specific VHHs, a naive library was screened for a variety of targets. At least two VHH candidates were extracted for each target using a 20-day selection pipeline. Some binders had ultrahigh potencies, with binding affinities in the nanomolar range. This naïve library, in particular, offers the possibility of acquiring unique antibodies targeting antigens of interest with low feasible dissociation constant (kD) without the time, effort, and price associated in producing antibodies in animals via antigen injection. Overall, the study shows that the megadiverse naïve library provides a rapid, adaptable, and easy platform for antibody creation, emphasizing its therapeutic and diagnostic implications.

Many approaches to identify therapeutically relevant neoantigens couple tumor sequencing with bioinformatic algorithms and inferred rules of tumor epitope immunogenicity. However, there are no reference data to compare these approaches, and the parameters governing tumor epitope immunogenicity remain unclear. Here, we assembled a global consortium wherein each participant predicted immunogenic epitopes from shared tumor sequencing data. 608 epitopes were subsequently assessed for T cell binding in patient-matched samples. By integrating peptide features associated with presentation and recognition, we developed a model of tumor epitope immunogenicity that filtered out 98% of non-immunogenic peptides with a precision above 0.70. Pipelines prioritizing model features had superior performance, and pipeline alterations leveraging them improved prediction performance. These findings were validated in an independent cohort of 310 epitopes prioritized from tumor sequencing data and assessed for T cell binding. This data resource enables identification of parameters underlying effective anti-tumor immunity and is available to the research community.

Background The advent of proprotein convertase subtilisin/kexin type 9 (PCSK9)–inhibiting drugs have provided an effective, but extremely expensive treatment for the management of low density lipoprotein (LDL). Our aim was to explore a cost-effective application of camelid anti-PCSK9 single domain antibodies (sdAbs), which are high variable regions of the camelid heavy chain antibodies (VHHs), as a human PCSK9 (hPCSK9) inhibitor. One female llama was immunized with hPCSK9. Screening of high affinity anti-PCSK9 VHHs was carried out based on surface plasmon resonance (SPR) technology. We reported a lysate kinetic analysis method improving the screening efficiency. To increase the serum half-life and targeting properties, the constant region fragment of the human immunoglobulin gamma sub-type 4 (IgG4 Fc) was incorporated to form a novel llama-human chimeric molecule (VHH-hFc). Results The PCSK9 inhibiting effects of the VHH proteins were analyzed in two human liver hepatocellular cells (HepG2 and Huh7) and in the hPCSK9 transgenic Sprague–Dawley (SD) rat model. The hPCSK9 antagonistic potency of the bivalent VHH-hFc exceeded the monovalent VHH (P

Additional file 4: Table S1. The sequences of the sdAbs.

Sugar beet is widely grown throughout the world and represents the second largest crop used to produce sugar. Root rot in sugar beet, caused by Fusarium, significantly reduces yield, juice purity, and sugar concentration. Here, 307 Fusarium isolates were collected from sugar beet roots exhibiting typical root rot symptoms in eight provinces or autonomous regions of China from 2009 to 2012. Based on morphological characteristics and sequence data of the internal transcribed spacer (ITS) region of ribosomal DNA (rDNA) and the translation elongation factor 1α (EF-1α), Fusarium oxysporum (38.4%) was identified as the most prevalent species, followed by F. solani (20.9%), and F. equiseti (18.9%). These three species were widely distributed in all eight of the provinces and autonomous regions. F. tricinctum (5.9%), F. brachygibbosum (4.6%), F. redolens (3.3%), F. proliferatum (3.3%), F. graminearum (2.3%), F. verticillioides (1.6%), F. nygamai (0.7%), and F. culmorum (0.3%) were less frequently obtained. Of the 307 Fusarium isolates, 117 representing different species and geographic locations were demonstrated to cause tip rot and vascular discoloration in sugar beet roots, with disease incidence ranging from 84.2 to 100.0% and disease index ranging from 41.94 to 75.83. This is the first detailed report of Fusarium species, in particular F. tricinctum, F. brachygibbosum, F. redolens, F. proliferatum, F. nygamai, and F. culmorum, causing sugar beet root rot in China.

Sauropus spatulifolius Beille is one kind of Chinese herbal medicine with anti-inflammatory and analgesic activities. In this study, we reported the complete chloroplast genome of S. spatulifolius ...

Cytokeratin 18 (CK18), the main scaffold protein of keratinocyte, is distributed in epithelial cells. This structural protein maintains the integrity and continuity of epithelial tissue. Cytokeratin is also frequently used as an immunohistochemical marker of tumor growth. In recent years, immune repertoire (IR) evaluation using next-generation sequencing (NGS) have become increasingly efficient. Here we deep sequenced the mouse IR of the immunoglobulin heavy chain (IGH) after CK18 immunization. We comprehensively analyzed the IR based on complementarity determining region 3 (CDR3) abundance, germline gene usage polarization, clone diversity, and lineage. We found many convergence characteristics after CK18 immunization. Convergence represents a phenomenon that antigen stimulation or pathogen exposure induces the antigen specific clone expansion and enrichment. The convergence could be used for the immune evaluation and antibody screen. After immunization, the IGHV5 gene clusters became preponderant. The abundance and length of the most frequent CDR3 both increased, nevertheless the IR diversity level decreased. From the convergent IGH repertoires, we selected and expressed six antibodies with the most frequent CDR3s and IGH V-J combinations. The ELISA results suggested all screened six antibodies bound CK18 specifically. The most potential antibody had 9.424E-10M M affinity for the interaction with the CK18. Therefore, this is the NGS platform has been first used for anti-CK18 monoclonal antibodies (MAbs) discovery. These analyses methods could also be used for vaccine evaluation.

A total of 32,284 unigenes were obtained from the transcriptome of Alternaria tenuissima, a pathogenic fungus causing foliar disease in tomato, using next-generation sequencing (NGS) technology. In total, 24,670 unigenes were annotated using five databases, including NCBI non-redundant protein, Swiss-Prot, euKaryotic Orthologous Groups, Kyoto Encyclopedia of Genes and Genomes, and the Gene Ontology. A total of 1,140 simple sequence repeats were also identified for use as molecular markers. Sixteen of the simple sequence repeat loci were selected to study the population structure of A. tenuissima. A population genetic analysis of 191 A. tenuissima isolates, sampled from four geographic regions in China, indicated that A. tenuissima had a high level of genetic diversity, and that the selected simple sequence repeat markers could reliably capture the genetic variation. The null hypothesis of random mating was rejected for all four geographic regions in China. Isolation by distance was observed for the entire data set, but not within clusters, which is indicative of barriers to gene flow among geographic regions. The analyses of Bayesian and principal coordinates, however, did not separate four geographic regions into four separate genetic clusters. The different levels of historical migration rates suggest that isolation by distance did not represent a major biological obstacle to the spread of A. tenuissima. The potential epidemic spread of A. tenuissima in China may occur through the transport of plant products or other factors. The presented results provide a basis for a comprehensive understanding of the population genetics of A. tenuissima in China.

Summary Chondrichthyan (cartilaginous fish) occupies a key phylogenetic position and is important for investigating evolutionary processes of vertebrates. However, limited whole genomes impede our in-depth knowledge of important issues such as chromosome evolution and immunity. Here, we report the chromosome-level genome of white-spotted bamboo shark. Combing it with other shark genomes, we reconstructed 16 ancestral chromosomes of bamboo shark and illustrate a dynamic chromosome rearrangement process. We found that genes on 13 fast-evolving chromosomes can be enriched in immune-related pathways. And two chromosomes contain important genes that can be used to develop single-chain antibodies, which were shown to have high affinity to human disease markers by using enzyme-linked immunosorbent assay. We also found three bone formation-related genes were lost due to chromosome rearrangements. Our study highlights the importance of chromosome rearrangements, providing resources for understanding of cartilaginous fish diversification and potential application of single-chain antibodies., Graphical Abstract, Highlights • Inferred ancestral chromosome karyotypes of cartilaginous fish • Chromosome rearrangements resulted in fast-evolving chromosomes and immune genes • Chromosome rearrangements led to deletion of bone formation-related genes • Proved that single-domain antibodies in shark have great potential application, Biological Sciences; Genetics; Genomics; Phylogenetics; Evolutionary Biology

Objective To develop a set of methods of amplifying the natural paring heavy and light chain genes from one single B cell. Methods Yanhuang (YH) cells were the first whole genome sequenced human cells of Asian origin. With the immortalized cell lines as the raw material, single CD19

Predicting peptide binding affinity with human leukocyte antigen (HLA) is a crucial step in developing powerful antitumor vaccine for cancer immunotherapy. Currently available methods work quite well in predicting peptide binding affinity with HLA alleles such as HLA-A*0201, HLA-A*0101, and HLA-B*0702 in terms of sensitivity and specificity. However, quite a few types of HLA alleles that are present in the majority of human populations including HLA-A*0202, HLA-A*0203, HLA-A*6802, HLA-B*5101, HLA-B*5301, HLA-B*5401, and HLA-B*5701 still cannot be predicted with satisfactory accuracy using currently available methods. Furthermore, currently the most popularly used methods for predicting peptide binding affinity are inefficient in identifying neoantigens from a large quantity of whole genome and transcriptome sequencing data. Here we present a Position Specific Scoring Matrix (PSSM)-based software called PSSMHCpan to accurately and efficiently predict peptide binding affinity with a broad coverage of HLA class I alleles. We evaluated the performance of PSSMHCpan by analyzing 10-fold cross-validation on a training database containing 87 HLA alleles and obtained an average area under receiver operating characteristic curve (AUC) of 0.94 and accuracy (ACC) of 0.85. In an independent dataset (Peptide Database of Cancer Immunity) evaluation, PSSMHCpan is substantially better than the popularly used NetMHC-4.0, NetMHCpan-3.0, PickPocket, Nebula, and SMM with a sensitivity of 0.90, as compared to 0.74, 0.81, 0.77, 0.24, and 0.79. In addition, PSSMHCpan is more than 197 times faster than NetMHC-4.0, NetMHCpan-3.0, PickPocket, sNebula, and SMM when predicting neoantigens from 661 263 peptides from a breast tumor sample. Finally, we built a neoantigen prediction pipeline and identified 117 017 neoantigens from 467 cancer samples of various cancers from TCGA. PSSMHCpan is superior to the currently available methods in predicting peptide binding affinity with a broad coverage of HLA class I alleles.

Compared to classical antibodies, camel heavy chain antibodies (HCAbs) are smaller in size due to lack of the light chain and the first constant domain of the heavy chain (CH1 region). The variable regions of HCAbs (VHHs) are more soluble and stable than that of conventional antibodies (VHs). Even with such simple structure, they are still functional in antigen binding. Although HCAbs have been extensively investigated over the past two decades, most efforts have been based upon low throughput sequence analysis, and there are only limited reports trying to analyze and describe the complete immune repertoire (IR) of camel HCAbs. Here we leveraged the high-throughput data generated by Next Generation Sequencing (NGS) of the variable domains of the antibody heavy chains from three Bactrian camels to conduct in-depth comparative analyses of the immunoglobulin repertoire. These include analyses of the complementary determining region 3 (CDR3) length and distribution, mutation rate, antibody characteristic amino acids, the distribution of the cysteine (Cys) codons, and the non-classical VHHs. We found that there is higher diversity in the CDR2 than in the other sub-regions, and there is a higher mutation rate in the VHHs than in the VHs (P < 0.05). In addition to substitutions at amino acid (AA) residue positions NO.49/50/52 between VH and VHH clones, we also observed other substitutions at the positions NO.40/54/57/96/101 that could lead to additional structural alterations. We also found that VH-derived VHH clones, referred to as non-classical VHH clones in this study, accounted for about 8% of all clones. Further, only 5%-10% clones had the Trp > Arg AA substitution at the first position of framework 4 for all types of clones. We present, for the first time, a relatively complete picture of the Bactrian camel antibody immune repertoire, including conventional antibody (Ab) and HCAbs, using PCR and in silico analysis based on high-throughput NGS data.

Present whole-cell patch-clamp methodology has only moderate consistency and throughput, rendering impractical functional measurements on large numbers of ion channel ligands or on large numbers of unknown or mutant channel genes. In the population patch clamp (PPC) described herein, a single voltage-clamp amplifier sums the whole-cell currents from multiple cells at once, each sealed to a separate aperture in a planar substrate well. The resulting ensemble currents are more consistent from well to well, and the success rate for each recording attempt is >95%. The PPC was implemented by modifying the PatchPlate substrate and amplifiers in the IonWorks patch-clamp instrument. The increased data consistency and likelihood of a successful recording in each well, combined with 384-well measurements in parallel, allow the direct electrophysiological recording of thousands of ensemble ionic currents per day. Therapeutic groups in drug discovery programs require this order of throughput to screen directed compound libraries against ion channel targets. The potential for studying the function of large numbers of ion channel mutants may be realized with the technique. The procedure incorporates subtraction methods that correct for expected distortions and also reliably produces data that agree with previous patch-clamp studies.

Voltage-dependent movement of a sodium channel S4 segment was examined by cysteine scanning mutagenesis and testing accessibility of the residues to hydrophilic cysteine-modifying reagents. These experiments indicate that 2 charged S4 residues move completely from an internally accessible to an externally accessible location in response to depolarization by passage through a short “channel” in the protein. The energetic problems of S4 movement have thus been solved in the same way that many ion channels achieve highly selective and rapid ion permeation through an open pore, by restricting the contact region between the permion and its channel.

Mutations in the skeletal muscle voltage-gated Na+ channel alpha-subunit have been found in patients with two distinct hereditary disorders of sarcolemmal excitation: hyperkalemic periodic paralysis (HYPP) and paramyotonia congenita (PC). Six of these mutations have been functionally expressed in a heterologous cell line (tsA201 cells) using the recombinant human skeletal muscle Na+ channel alpha-subunit cDNA hSkM1. PC mutants from diverse locations in this subunit (T1313M, L1433R, R1448H, R1448C, A1156T) all exhibit a similar disturbance in channel inactivation characterized by reduced macroscopic rate, accelerated recovery, and altered voltage dependence. PC mutants had no significant abnormality in activation. In contrast, one HYPP mutation studied (T704M) has a normal inactivation rate but exhibits shifts in the midpoints of steady-state activation and inactivation along the voltage axis. These findings help to explain the phenotypic differences between HYPP and PC at the molecular and biophysical level and contribute to our understanding of Na+ channel structure and function.

T-type Ca(2+) channels encoded by voltage-gated Ca(2+) channel (Ca(v)) 3.1, 3.2, and 3.3 genes play important physiological roles and serve as therapeutic targets for neurological and cardiovascular disorders. Currently there is no selective T-channel blocker. To screen for such a blocker, we developed three stable cell lines expressing human recombinant Ca(v)3.1, 3.2, or 3.3 channels and then examined their usefulness in high throughput screens. All three cell lines displayed an increase in intracellular Ca(2+) in response to changes in extracellular Ca(2+) as detected with Ca(2+)-sensitive dyes using a fluorometric imaging plate reader (FLIPR [Molecular Devices, Sunnyvale, CA] or FlexStation [Molecular Devices]). The signal-to-noise ratio was 2-4. Co-expression of Ca(v)3.2 with a mouse leak K(+) channel, which by virtue of being open at rest hyperpolarizes the cell membrane, blocked the fluorescent signal. Co-addition of KCl to these cells induced a Ca(2+) signal that was similar to that observed in the cell line expressing Ca(v)3.2 alone. These results confirm that the detection of intracellular Ca(2+) increase in cells expressing Ca(v)3.2 alone results from Ca(2+) entry through channels that are open at the resting membrane potential of each cell line (i.e., window currents). Testing known drugs on Ca(v)3 channels showed that block could be reliably detected using the FlexStation assay, FLIPR assay, or voltage clamp recordings using the IonWorks HT system (Molecular Devices). These results support the use of the FLIPR window current assay for primary drug screening and high throughput patch recordings for secondary screening of novel T-channel blockers.

The mutation R1448C substitutes a cysteine for the outermost arginine in the fourth transmembrane segment (S4) of domain 4 in skeletal muscle sodium channels. We tested the accessibility of this cysteine residue to hydrophilic methanethiosulfonate reagents applied to the extracellular surface of cells expressing these mutant channels. The reagents irreversibly increase the rate of inactivation of R1 448C, but not wild-type, channels. Cysteine modification is voltage dependent, as if depolarization extends this residue into the extracellular space. The rate of cysteine modification increases with depolarization and has the voltage dependence and kinetics expected for the movement of a voltage sensor controlling channel gating.

The second and third basic residues of the S4 segment of domain 4 (D4:R2 and D4:R3) of the human skeletal muscle Na+ channel are known to be translocated from a cytoplasmic to an extracellular position during depolarization. Accessibilities of individual S4 residues were assayed by alteration of inactivation kinetics during modification of cysteine mutants by hydrophilic methanethiosulfonate reagents. The voltage dependences of the reaction rates are identical for extracellular application of cationic methanethiosulfonate-ethyltrimethylammonium (MTSET) and anionic methanethiosulfonate-ethylsulfonate (MTSES), suggesting that D4:R3C is situated outside the membrane electric field at depolarized voltages. The absolute rate of R3C modification is 281-fold greater for MTSET than for MTSES, however, suggesting that at depolarized voltages this S4 thiol resides in a negatively charged hydrophilic crevice. The two hydrophobic residues between D4:R2C and D4:R3C in the primary sequence (L1452 and A1453) are not externally exposed at any voltage. An alpha-helical representation of D4/S4 shows that the basic residues D4:R2 and D4:R3 are on the face opposite that of L1452 and A1453. We propose that in the depolarized conformation, the hydrophobic face of this portion of D4/S4 remains in contact with a hydrophobic region of the extracellular vestibule of the S4 channel.

Supplementary methods. (PDF 608 kb)

Supplementary figures. (PDF 1140 kb)

Supplementary figures. (PDF 1140 kb)

Supplementary methods. (PDF 608 kb)