Application of circular consensus sequencing and network analysis to characterize the bovine IgG repertoire

Analysis of immunoglobulin organization and complexity in mink (Neovison vison)

Mink are susceptible to viruses such as SARS-CoV-2, H1N1 and H9N2, so they are considered a potential animal model for studying human viral infections. Therefore, it is important to study the immune system of mink. Immunoglobulin (Ig) is an important component of humoral immunity and plays an important role in the body's immune defense. In this study, we described the gene loci structure of mink Ig germline by genome comparison, and analysed the mechanism of expression diversity of mink antibody library by 5'RACE and next-generation sequencing (NGS). The results were as follows: the IgH, Igκ and Igλ loci of mink were located on chromosome 13, chromosome 8 and chromosome 3, respectively, and they had 25, 36 and 7 V genes, 3, 5 and 7 J genes and 10 DH genes, respectively. Mink Ig heavy chain preferred the IGHV1, IGHD2 and IGHJ4 subgroups, κ chain mainly use the IGKV1, IGKJ1 and IGHL4 subgroups, and λ chain mainly use the IGLV3 and IGLJ3 subgroups. Linkage diversity analysis revealed that N nucleotide insertion was the main factor affecting the linkage diversity of mink Igs. On the mutation types of mink Ig Somatic Hypermutation (SHM), the high mutation types of heavy chain were mainly G > A, C > T, T > C, A > G, C > A, G > T, A > C, and T > G; the high mutation types of κ chain were G > A and T > C; and the high mutation types of λ chain were G > A and A > G. The objective of this study was to analyse the loci structure and expression diversity of Ig in mink. The results contribute to our comprehension of Ig expression patterns in mink and were valuable for advancing knowledge in mink immunogenetics, exploring the evolution of adaptive immune systems across different species, and conducting comparative genomics research.

Variations in antibody repertoires correlate with vaccine responses

An important challenge in vaccine development is to figure out why a vaccine succeeds in some individuals and fails in others. Although antibody repertoires hold the key to answering this question, there have been very few personalized immunogenomics studies so far aimed at revealing how variations in immunoglobulin genes affect a vaccine response. We conducted an immunosequencing study of 204 calves vaccinated against bovine respiratory disease (BRD) with the goal to reveal variations in immunoglobulin genes and somatic hypermutations that impact the efficacy of vaccine response. Our study represents the largest longitudinal personalized immunogenomics study reported to date across all species, including humans. To analyze the generated data set, we developed an algorithm for identifying variations of the immunoglobulin genes (as well as frequent somatic hypermutations) that affect various features of the antibody repertoire and titers of neutralizing antibodies. In contrast to relatively short human antibodies, cattle have a large fraction of ultralong antibodies that have opened new therapeutic opportunities. Our study reveals that ultralong antibodies are a key component of the immune response against the costliest disease of beef cattle in North America. The detected variants of the cattle immunoglobulin genes, which are implicated in the success/failure of the BRD vaccine, have the potential to direct the selection of individual cattle for ongoing breeding programs.

Screening and Validation of p38 MAPK Involved in Ovarian Development of Brachymystax lenok

Brachymystax lenok (lenok) is a rare cold-water fish native to China that is of high meat quality. Its wild population has declined sharply in recent years, and therefore, exploring the molecular mechanisms underlying the development and reproduction of lenoks for the purposes of artificial breeding and genetic improvement is necessary. The lenok comparative transcriptome was analyzed by combining single molecule, real-time, and next generation sequencing (NGS) technology. Differentially expressed genes (DEGs) were identified in five tissues (head kidney, spleen, liver, muscle, and gonad) between immature [300 days post-hatching (dph)] and mature [three years post-hatching (ph)] lenoks. In total, 234,124 and 229,008 full-length non-chimeric reads were obtained from the immature and mature sequencing data, respectively. After NGS correction, 61,405 and 59,372 non-redundant transcripts were obtained for the expression level and pathway enrichment analyses, respectively. Compared with the mature group, 719 genes with significantly increased expression and 1,727 genes with significantly decreased expression in all five tissues were found in the immature group. Furthermore, DEGs and pathways involved in the endocrine system and gonadal development were identified, and p38 mitogen-activated protein kinases (MAPKs) were identified as potentially regulating gonadal development in lenok. Inhibiting the activity of p38 MAPKs resulted in abnormal levels of gonadotropin-releasing hormone, follicle-stimulating hormone, and estradiol, and affected follicular development. The full-length transcriptome data obtained in this study may provide a valuable reference for the study of gene function, gene expression, and evolutionary relationships in B. lenok and may illustrate the basic regulatory mechanism of ovarian development in teleosts.

Access to ultra-long IgG CDRH3 bovine antibody sequences using short read sequencing technology

The advances in high-throughput DNA sequencing and recombinant antibody technologies has presented new methods for characterizing antibody repertoires and significantly increased our understanding on the functional role of antibodies in immunity and their use in diagnostics, vaccine antigen design and as biological therapeutics. A subset of Bos taurus antibodies possesses unique ultra-long third complementary-determining region of the heavy chain (CDRH3) and are of special interest because they are thought to have unique functional abilities of broadly neutralizing properties - a functional role that has not been fully explored in vaccine development. Next generation sequencing technologies that are widely used to profile immunoglobulin (Ig) repertoires are based on short-read methods such as the Illumina technology. Although this technology has worked well in sequencing Ig V-D-J regions of most jawed vertebrates, it has faced serious technical challenges with sequencing regions in bovine Ig bearing ultra-long CDRH3 sequences, which are longer than 120 bp. To overcome this limitation, we have developed a sequencing strategy based on nested PCR products that allows sequence assembly of full-length bovine Ig heavy-chain (IgH) V-D-J regions. We have used this strategy to sequence IgH V-D-J regions of two Bos indicus breeds, Ankole and Boran. We confirm the presence of ultra-long CDRH3 sequences in IgG transcripts in both African cattle breeds, and provide preliminary evidence for differences and preferences in germline V_H, D_H and J_H allele gene usage as well as differences in the length of the VH region in the two bovine breeds. Our method provides tools that should allow more robust analyses of ultra-long CDRH3 sequences aiding antibody and epitope discovery in different cattle breeds and their role in mediating immunity.

Isoform Age - Splice Isoform Profiling Using Long-Read Technologies

Alternative splicing (AS) of RNA is a key mechanism that results in the expression of multiple transcript isoforms from single genes and leads to an increase in the complexity of both the transcriptome and proteome. Regulation of AS is critical for the correct functioning of many biological pathways, while disruption of AS can be directly pathogenic in diseases such as cancer or cause risk for complex disorders. Current short-read sequencing technologies achieve high read depth but are limited in their ability to resolve complex isoforms. In this review we examine how long-read sequencing (LRS) technologies can address this challenge by covering the entire RNA sequence in a single read and thereby distinguish isoform changes that could impact RNA regulation or protein function. Coupling LRS with technologies such as single cell sequencing, targeted sequencing and spatial transcriptomics is producing a rapidly expanding suite of technological approaches to profile alternative splicing at the isoform level with unprecedented detail. In addition, integrating LRS with genotype now allows the impact of genetic variation on isoform expression to be determined. Recent results demonstrate the potential of these techniques to elucidate the landscape of splicing, including in tissues such as the brain where AS is particularly prevalent. Finally, we also discuss how AS can impact protein function, potentially leading to novel therapeutic targets for a range of diseases.

Guidelines for Setting Up a mRNA Sequencing Experiment and Best Practices for Bioinformatic Data Analysis

RNA-sequencing, commonly referred to as RNA-seq, is the most recently developed method for the analysis of transcriptomes. It uses high-throughput next-generation sequencing technologies and has revolutionized our understanding of the complexity and dynamics of whole transcriptomes.In this chapter, we recall the key developments in transcriptome analysis and dissect the different steps of the general workflow that can be run by users to design and perform a mRNA-seq experiment as well as to process mRNA-seq data obtained by the Illumina technology. The chapter proposes guidelines for completing a mRNA-seq study properly and makes available recommendations for best practices based on recent literature and on the latest developments in technology and algorithms. We also remark the large number of choices available (especially for bioinformatic data analysis) in front of which the scientist may be in trouble.In the last part of the chapter we discuss the new frontiers of single-cell RNA-seq and isoform sequencing by long read technology.

Structural Diversity of Ultralong CDRH3s in Seven Bovine Antibody Heavy Chains

Antigen recognition by mammalian antibodies represents the most diverse setting for protein-protein interactions, because antibody variable regions contain exceptionally diverse variable gene repertoires of DNA sequences containing combinatorial, non-templated junctional mutational diversity. Some animals use additional strategies to achieve structural complexity in the antibody combining site, and one of the most interesting of these is the formation of ultralong heavy chain complementarity determining region 3 loops in cattle. Repertoire sequencing studies of bovine antibody heavy chain variable sequences revealed that bovine antibodies can contain heavy chain complementarity determining region 3 (CDRH3) loops with 60 or more amino acids, with complex structures stabilized by multiple disulfide bonds. It is clear that bovine antibodies can achieve long, peculiarly structured CDR3s, but the range of diversity and complexity of those structures is poorly understood. We determined the atomic resolution structure of seven ultralong bovine CDRH3 loops. The studies, combined with five previous structures, reveal a large diversity of cysteine pairing variations, and highly diverse globular domains.

Current and Future Methods for mRNA Analysis: A Drive Toward Single Molecule Sequencing

The transcriptome encompasses a range of species including messenger RNA, and other noncoding RNA such as rRNA, tRNA, and short and long noncoding RNAs. Due to the huge role played by mRNA in development and disease, several methods have been developed to sequence and characterize mRNA, with RNA sequencing (RNA-Seq) emerging as the current method of choice particularly for large high-throughput studies. Short-read RNA-Seq which involves sequencing of short cDNA fragments and computationally assembling them to reconstruct the transcriptome, or aligning them to a reference is the most widely used approach. However, due to inherent limitations of this approach in de novo transcriptome assembly and isoform quantification, long-read RNA-Seq approaches, which also happen to be single molecule sequencing approaches, are increasingly becoming the standard for de novo transcriptome assembly and isoform quantification. In this chapter, we review the technical aspects of the current methods of RNA-Seq, both short and long-read approaches, and data analysis methods available. We discuss recent advances in single-cell RNA-Seq and direct RNA-Seq approaches, which perhaps will dominate the future of RNA-Seq.

How B-Cell Receptor Repertoire Sequencing Can Be Enriched with Structural Antibody Data

Next-generation sequencing of immunoglobulin gene repertoires (Ig-seq) allows the investigation of large-scale antibody dynamics at a sequence level. However, structural information, a crucial descriptor of antibody binding capability, is not collected in Ig-seq protocols. Developing systematic relationships between the antibody sequence information gathered from Ig-seq and low-throughput techniques such as X-ray crystallography could radically improve our understanding of antibodies. The mapping of Ig-seq datasets to known antibody structures can indicate structurally, and perhaps functionally, uncharted areas. Furthermore, contrasting naïve and antigenically challenged datasets using structural antibody descriptors should provide insights into antibody maturation. As the number of antibody structures steadily increases and more and more Ig-seq datasets become available, the opportunities that arise from combining the two types of information increase as well. Here, we review how these data types enrich one another and show potential for advancing our knowledge of the immune system and improving antibody engineering.

A generalized quantitative antibody homeostasis model: antigen saturation, natural antibodies and a quantitative antibody network

In a pair of articles, we present a generalized quantitative model for the homeostatic function of clonal humoral immune system. In this second paper, we describe how antibody production controls the saturation of antigens and the network of antibody interactions that emerges in the epitome space with the establishment of the immune system. Efficient control of antigens, be it self or foreign, requires the maintenance of antibody concentrations that saturate antigen to relevant levels. Simple calculations suggest that the observed diverse recognition of antigens by natural antibodies is only possible by cross-reactivity whereby particular clones of antibodies bind to diverse targets and shared recognition of particular antigens by multiple antibody clones contribute to the maintenance of antigen control. We also argue that natural antibodies are none else than the result of thymus-independent responses against immunological self. We interpret and explain antibody production and function in a virtual molecular interaction space and as a network of interactions. Indeed, the general quantitative (GQM) model we propose is in agreement with earlier models, confirms some assumptions and presumably provides the theoretical basis for the construction of a real antibody network using the sequence and interaction database data.

Development of a Bioinformatics Framework for the Detection of Gene Conversion and the Analysis of Combinatorial Diversity in Immunoglobulin Heavy Chains in Four Cattle Breeds

We have developed a new bioinformatics framework for the analysis of rearranged bovine heavy chain immunoglobulin (Ig) variable regions by combining and refining widely used alignment algorithms. This bioinformatics framework allowed us to investigate alignments of heavy chain framework regions (FRHs) and the separate alignments of FRHs and heavy chain complementarity determining regions (CDRHs) to determine their germline origin in the four cattle breeds Aubrac, German Black Pied, German Simmental, and Holstein Friesian. Now it is also possible to specifically analyze Ig heavy chains possessing exceptionally long CDR3Hs. In order to gain more insight into breed specific differences in Ig combinatorial diversity, somatic hypermutations and putative gene conversions of IgG, we compared the dominantly transcribed variable (IGHV), diversity (IGHD), and joining (IGHJ) segments and their recombination in the four cattle breeds. The analysis revealed the use of 15 different IGHV segments, 21 IGHD segments, and two IGHJ segments with significant different transcription levels within the breeds. Furthermore, there are preferred rearrangements within the three groups of CDR3H lengths. In the sequences of group 2 (CDR3H lengths (L) of 11–47 amino acid residues (aa)) a higher number of recombination was observed than in sequences of group 1 (L≤10 aa) and 3 (L≥48 aa). The combinatorial diversity of germline IGHV, IGHD, and IGHJ-segments revealed 162 rearrangements that were significantly different. The few preferably rearranged gene segments within group 3 CDR3H regions may indicate specialized antibodies because this length is unique in cattle. The most important finding of this study, which was enabled by using the bioinformatics framework, is the discovery of strong evidence for gene conversion as a rare event using pseudogenes fulfilling all definitions for this particular diversification mechanism.

Internal Duplications of DH, JH, and C Region Genes Create an Unusual IgH Gene Locus in Cattle

It has been suspected for many years that cattle possess two functional IgH gene loci, located on Bos taurus autosome (BTA) 21 and BTA11, respectively. In this study, based on fluorescence in situ hybridization and additional experiments, we showed that all functional bovine IgH genes were located on BTA21, and only a truncated μCH2 exon was present on BTA11. By sequencing of seven bacterial artificial chromosome clones screened from a Hostein cow bacterial artificial chromosome library, we generated a 678-kb continuous genomic sequence covering the bovine IGHV, IGHD, IGHJ, and IGHC genes, which are organized as IGHVn-IGHDn-IGHJn-IGHM1-(IGHDP-IGHV3-IGHDn)3-IGHJn-IGHM2-IGHD-IGHG3-IGHG1-IGHG2-IGHE-IGHA. Although both of two functional IGHM genes, IGHM1 and IGHM2, can be expressed via independent VDJ recombinations, the IGHM2 can also be expressed through class switch recombination. Likely because more IGHD segments can be involved in the expression of IGHM2, the IGHM2 gene was shown to be dominantly expressed in most tissues throughout different developmental stages. Based on the length and identity of the coding sequence, the 23 IGHD segments identified in the locus could be divided into nine subgroups (termed IGHD1 to IGHD9). Except two members of IGHD9 (14 nt in size), all other functional IGHD segments are longer than 30 nt, with the IGHD8 gene (149 bp) to be the longest. These remarkably long germline IGHD segments play a pivotal role in generating the exceptionally great H chain CDR 3 length variability in cattle.

High-accuracy de novo assembly and SNP detection of chloroplast genomes using a SMRT circular consensus sequencing strategy

A circular consensus sequencing (CCS) strategy involving single molecule, real-time (SMRT) DNA sequencing technology was applied to de novo assembly and single nucleotide polymorphism (SNP) detection of chloroplast genomes. Chloroplast DNA was purified from enriched chloroplasts of pooled individuals to construct a shotgun library for each species. The sequencing reactions were performed on a PacBio RS platform. CCS sub-reads were generated from polymerase reads that passed the native dumbbell-shaped DNA templates multiple times. The complete chloroplast genome sequence was generated by mapping all reads to the draft sequence constructed in a step-by-step manner. The full-chain, PCR-free approach eliminates the possible context-specific biases in library construction and sequencing reaction. The chloroplast genome was easily and completely assembled using the data generated from one SMRT Cell without requiring a reference genome. Comparisons of the three assembled Fritillaria genomes to 34.1 kb of validation Sanger sequences revealed 100% concordance, and the detected intraspecies SNPs at a minimum variant frequency of 15% were all confirmed. This simple approach with potential for parallel sequencing yields high-quality chloroplast genomes for sensitive SNP detection and comparative analyses. We recommend this approach for its powerful applicability for evolutionary genetics and genomics studies in plants based on the sequences of chloroplast genomes.

The utility of PacBio circular consensus sequencing for characterizing complex gene families in non-model organisms

BACKGROUND

Molecular characterization of highly diverse gene families can be time consuming, expensive, and difficult, especially when considering the potential for relatively large numbers of paralogs and/or pseudogenes. Here we investigate the utility of Pacific Biosciences single molecule real-time (SMRT) circular consensus sequencing (CCS) as an alternative to traditional cloning and Sanger sequencing PCR amplicons for gene family characterization. We target vomeronasal gene receptors, one of the most diverse gene families in mammals, with the goal of better understanding intra-specific V1R diversity of the gray mouse lemur (Microcebus murinus). Our study compares intragenomic variation for two V1R subfamilies found in the mouse lemur. Specifically, we compare gene copy variation within and between two individuals of M. murinus as characterized by different methods for nucleotide sequencing. By including the same individual animal from which the M. murinus draft genome was derived, we are able to cross-validate gene copy estimates from Sanger sequencing versus CCS methods.

RESULTS

We generated 34,088 high quality circular consensus sequences of two diverse V1R subfamilies (here referred to as V1RI and V1RIX) from two individuals of Microcebus murinus. Using a minimum threshold of 7× coverage, we recovered approximately 90% of V1RI sequences previously identified in the draft M. murinus genome (59% being identical at all nucleotide positions). When low coverage sequences were considered (i.e. < 7× coverage) 100% of V1RI sequences identified in the draft genome were recovered. At least 13 putatively novel V1R loci were also identified using CCS technology.

CONCLUSIONS

Recent upgrades to the Pacific Biosciences RS instrument have improved the CCS technology and offer an alternative to traditional sequencing approaches. Our results suggest that the Microcebus murinus V1R repertoire has been underestimated in the draft genome. In addition to providing an improved understanding of V1R diversity in the mouse lemur, this study demonstrates the utility of CCS technology for characterizing complex regions of the genome. We anticipate that long-read sequencing technologies such as PacBio SMRT will allow for the assembly of multigene family clusters and serve to more accurately characterize patterns of gene copy variation in large gene families, thus revealing novel micro-evolutionary patterns within non-model organisms.

Systematic characterization and comparative analysis of the rabbit immunoglobulin repertoire

Rabbits have been used extensively as a model system for the elucidation of the mechanism of immunoglobulin diversification and for the production of antibodies. We employed Next Generation Sequencing to analyze Ig germline V and J gene usage, CDR3 length and amino acid composition, and gene conversion frequencies within the functional (transcribed) IgG repertoire of the New Zealand white rabbit (Oryctolagus cuniculus). Several previously unannotated rabbit heavy chain variable (VH) and light chain variable (VL) germline elements were deduced bioinformatically using multidimensional scaling and k-means clustering methods. We estimated the gene conversion frequency in the rabbit at 23% of IgG sequences with a mean gene conversion tract length of 59±36 bp. Sequencing and gene conversion analysis of the chicken, human, and mouse repertoires revealed that gene conversion occurs much more extensively in the chicken (frequency 70%, tract length 79±57 bp), was observed to a small, yet statistically significant extent in humans, but was virtually absent in mice.

The promise and challenge of high-throughput sequencing of the antibody repertoire

Georgiou and colleagues discuss rapidly evolving methods for high-throughput sequencing of the antibody repertoire, and how the resulting data may be applied to answer basic and translational research questions.

Efforts to determine the antibody repertoire encoded by B cells in the blood or lymphoid organs using high-throughput DNA sequencing technologies have been advancing at an extremely rapid pace and are transforming our understanding of humoral immune responses. Information gained from high-throughput DNA sequencing of immunoglobulin genes (Ig-seq) can be applied to detect B-cell malignancies with high sensitivity, to discover antibodies specific for antigens of interest, to guide vaccine development and to understand autoimmunity. Rapid progress in the development of experimental protocols and informatics analysis tools is helping to reduce sequencing artifacts, to achieve more precise quantification of clonal diversity and to extract the most pertinent biological information. That said, broader application of Ig-seq, especially in clinical settings, will require the development of a standardized experimental design framework that will enable the sharing and meta-analysis of sequencing data generated by different laboratories.

Efficacy of a 3rd generation high-throughput sequencing platform for analyses of 16S rRNA genes from environmental samples

Longer sequences of the bacterial 16S rRNA gene could provide greater phylogenetic and taxonomic resolutions and advance knowledge of population dynamics within complex natural communities. We assessed the accuracy of a Pacific Biosciences (PacBio) single molecule, real time (SMRT) sequencing based on DNA polymerization, a promising 3rd generation high-throughput technique, and compared this to the 2nd generation Roche 454 pyrosequencing platform. Amplicons of the 16S rRNA gene from a known isolate, Shewanella oneidensis MR1, and environmental samples from two streambed habitats, rocks and sediments, and a riparian zone soil, were analyzed. On the PacBio we analyzed ~500 bp amplicons that covered the V1-V3 regions and the full 1500 bp amplicons of the V1-V9 regions. On the Roche 454 we analyzed the ~500 bp amplicons. Error rates associated with the isolate were lowest with the Roche 454 method (2%), increased by more than 2-fold for the 500 bp amplicons with the PacBio SMRT chip (4-5%), and by more than 8-fold for the full gene with the PacBio SMRT chip (17-18%). Higher error rates with the PacBio SMRT chip artificially inflated estimates of richness and lowered estimates of coverage for environmental samples. The 3rd generation sequencing technology we evaluated does not provide greater phylogenetic and taxonomic resolutions for studies of microbial ecology.

Exceptionally long CDR3H are not isotype restricted in bovine immunoglobulins

Exceptionally long third complementarity determining regions of the heavy chain (CDR3H) were previously described as a specificity of bovine IgG and IgM immunoglobulins. In addition, the genomic organization of the immunoglobulin heavy chain locus remains to be elucidated with a special focus on the number of variable segments (IGHV). By analyzing the variable regions according to the isotype-specific PCR using cDNA-PCR, we were able to prove the existence of exceptional long CDR3H in all bovine isotypes. The corresponding sequences of three distinct amplicons were grouped according to the length of the CDR3H. Sequences of CDR3H possessed 5 to 10, 12 to 31 or at least 48 amino acid residues. Long and mid-length CDR3H were composed of mainly hydrophilic amino acid residues, while short CDR3H also contained hydrophobic amino acid residues. All sequences with long CDR3H were related to the germline variable segment 10. Using the current genome assembly, Bos taurus NCBI build 6.1, the genomic organization of the bovine immunoglobulin heavy-chain locus was analyzed. A main locus was investigated on BTA21. Exons coding for variable, diversity, and joining segments, as well as for the constant regions of different isotypes, were also localized on BTA7, BTA8, and BTA20. Together with the information from unplaced contigs, 36 IGHV were detected of which 13 are putatively functional. Phylogenetic analysis revealed two bovine IGHV families (boVH1, boVH2). Thus, the existence of the two bovine families suggested was demonstrated, where boVH1 comprises all functional segments. This study substantially improves the understanding of the generation of immunoglobulin diversity in cattle.