- Research article
- Open Access
Formation of ultralong DH regions through genomic rearrangement
BMC Immunology volume 21, Article number: 30 (2020)
Cow antibodies are very unusual in having exceptionally long CDR H3 regions. The genetic basis for this length largely derives from long heavy chain diversity (DH) regions, with a single “ultralong” DH, IGHD8–2, encoding over 50 amino acids. Many bovine IGHD regions have sequence similarity but have several nucleotide repeating units that diversify their lengths. Genomically, most DH regions exist in three clusters that appear to have formed from DNA duplication events. However, the relationship between the genomic arrangement and long CDR lengths is unclear.
The DH cluster containing IGHD8–2 underwent a rearrangement and deletion event in relation to the other clusters in the region corresponding to IGHD8–2, with possible fusion of two DH regions and expansion of short repeats to form the ultralong IGHD8–2 gene.
Length heterogeneity within DH regions is a unique evolutionary genomic mechanism to create immune diversity, including formation of ultralong CDR H3 regions.
Adaptive immunity arose in vertebrates through the ability to somatically alter antigen receptor (antibody and T-cell receptor) genes to form diverse repertoires which are selected to bind and neutralize invading pathogens. A key component of this system is the ability to perform recombination of variable (V), diversity (D), and joining (J) gene segments through the process of V(D)J recombination [1,2,3]. A diversity of V, D, and J elements, along with imprecise joining at the V-D and D-J junctions enables different amino acids to be encoded in key paratopic regions which impact antigen binding.
The third complementary determining region of the heavy chain (CDR H3) is particularly important in antibody molecules as it contains the greatest diversity and also usually makes the most extensive contact with antigen. Long CDR H3 regions are often found in broadly neutralizing antibodies targeting human immunodeficiency (HIV), influenza, and polio viruses [4,5,6,7,8], and are also thought to be important in binding challenging antigens like G-protein coupled receptors and protease active sites [9, 10]. Thus, genetic mechanisms to form long CDR H3s may be very important in immune responses against key antigens.
In most organisms, the antibody CDR H3 forms a loop of 10–15 amino acids in length, and is encoded by the DH gene and associated recombinational junctions that form through VDJ recombination. Unusually long CDR H3s, such as those in rare broadly neutralizing anti-HIV antibodies from infected patients, are often over 20 amino acids in length [4, 11,12,13]. The major determinants of CDR H3 length are the length of the germline encoded DH region, as well as somatic insertion of nucleotides (e.g. N- or P- nucleotides) at the V-D and D-J junctions. In humans, the longest DH region, IGHD3–16, encodes 12 amino acids.
Bovines are remarkable in having very long CDR H3 regions [14,15,16,17,18,19,20,21,22,23,24], with an average length of 26 amino acids but with an exceptionally long subset of the repertoire (the “ultralong” CDR H3 antibodies) that can have CDR H3 lengths of up to 70 amino acids. These CDR H3 regions form their own independently folding mini domains comprised of a β-ribbon “stalk” that protrudes far from the typical paratope surface upon which sits a disulfide-bonded “knob” [21, 23, 25, 26]. Cows are the only species thus far investigated that can produce a broadly neutralizing antibody response against HIV when immunized with an engineered spike protein immunogen as an experimental vaccine. These antibodies are characterized by ultralong CDR H3 regions that penetrate the glycan shield of the spike protein to bind a conserved broadly neutralizing epitope in the CD4 binding region . Cows are therefore unusual in producing long CDR H3s, and this unique repertoire has major functional relevance in neutralizing an antigen that is extremely challenging for repertoires of other mammalian species. Therefore, understanding the natural genetic and evolutionary mechanisms behind ultralong CDR H3 generation would be important in vaccine generation as well as therapeutic antibody discovery and development.
At least two evolutionary genetic events occurred which enabled formation of ultralong CDR H3 antibodies in cows. First, a unique VH region evolved as a result of an 8-basepair duplication at the 3′ end of IGHV1–7 . This particular variable region is the only VH region used in ultralong CDR H3 antibodies, and the short duplication directly encodes the ascending strand of the stalk region characteristic of these antibodies. Second, a very long DH region is found in cattle, IGHD8–2, which encodes 49 amino acids [23, 28,29,30]. Antibodies with ultralong CDR H3 regions invariably use IGHV1–7 and IGHD8–2 [8, 16]. Here we examine the genetic features at work in the evolution of this unusually long DH region of cattle.
DH cluster 2 has a significant deletion
We analyzed the DH regions of the recent assembly of the Bos taurus immunoglobulin heavy chain locus for features associated with the ultralong IGHD8–2 region. Of particular note, the DH regions at the heavy chain locus are divided into “clusters” that arose from duplication events through evolution. The IMGT naming nomenclature for DH regions includes numerical designations for the family and cluster of each gene; for example, IGHD3–2 is in family 3 and located in cluster 2 [16, 31,32,33]. There are four clusters, with clusters 2–4 being homologous with nucleotide identities of 92% (cluster 2 vs cluster 3), 99.7% (cluster 3 vs cluster 4), and 92% (cluster 2 vs cluster 4). The sequences of the DH regions located within the clusters are also homologous, with DH regions occupying analogous locations being 96 to 100% identical at the nucleotide level (Supplemental Fig. 1). A major discrepancy in the cluster sequences, however, is that cluster 2 (3480 nucleotides) is 358 and 364 nucleotides shorter than clusters 3 (3838 nt) and 4 (3844 nt), respectfully. Additionally, cluster 2 is comprised of only five DH regions, with one of them being the ultralong IGHD8–2, whereas clusters 3 and 4 are comprised of six DH regions (Fig. 1). Thus, cluster 2 appears to have a significant genomic deletion in relation to the highly homologous clusters 3 and 4. We hypothesized that this deletion might be related to formation of the ultralong IGHD8–2 region located in cluster 2. In simplistic terms, one explanation for formation of an ultralong DH region would be by fusion of two DH regions through deletion of intragenic sequence, with the fusion maintaining recombination signal sequences of each DH at both the 5′ and 3′ ends.
Cluster 2 has a short chromosomal rearrangement
To evaluate the location of the deletion in cluster 2 relative to clusters 3 and 4, we performed a series of sequence alignments of the clusters, the DH regions, and the intergenic regions (between DH regions). Indeed, the deletion in cluster 2 in relation to clusters 3 and 4 occurred at IGHD8–2, however the deletion was also associated with a larger chromosomal rearrangement. In this regard, IGHD5–2 in cluster 2 appears to have replaced the paralog for IGHD3–3 (cluster 3) and IGHD3–4 (cluster 4)(Fig. 1, Supplemental Figs. 2–3). The IGHD5 homologs are immediately 5′ of the IGHD6 family members in clusters 3 and 4, however IGHD5–2 is situated immediately 3′ of IGHD2–2 and immediately 5′ of the ultralong IGHD8–2 region in cluster 2 (Fig. 1). There is no IGHD3 family member in cluster 2 (Supplemental Fig. 3), with the paralog of IGHD3–3 and IGHD3–4 either deleted or fused to the adjacent DH region, which would be a paralog of IGHD7–3 (cluster 3) or IGHD7–4 (cluster 4). Global alignments of the clusters show deleted nucleotides at IGHD8–2 as well as the position occupied by family 5 genes in clusters 3 and 4 (e.g. between IGHD7 and IGHD6). Alignments of the intergenic regions show that the intergenic region corresponding to the sequence between IGHD3–3 and IGHD7–3 in cluster 3 (or IGHD3–4 and IGHD7–4 in cluster 4) is deleted in cluster 2 (Supplemental Fig. 4). While IGHD5–2 has been transposed to a location 3′ to IGHD8–2, the actual genetic material deleted clearly includes IGHD3 and its 3′ intergenic region. Thus, one possibility is that the ultralong IGHD8–2 region resulted from a deletion and associated fusion of the cluster 2 paralogs of IGHD3–3 and IGHD7–3. However, local sequence alignment reveals that the 5′ end of IGHD6–3 is 91.2% identical to IGHD8–2 (89.4% for IGHD6–2) over the first 85 nucleotides, whereas IGHD3–3 (and IGHD3–4) is only 80% identical over the first 62 residues (Supplemental Fig. 7). Of note, IGHD6 family sequences share a cysteine in the same position as the conserved cysteine in IGHD8–2, which is highly conserved in deep sequenced ultralong CDR H3 antibodies, and participates in a conserved disulfide bond at the base of the ultralong CDR H3 stalk [23, 26]. Thus, donation of an IGHD6 to the 5′ end of an IGHD7 through a recombinational or gene conversion process is a likely mechanism to produce IGHD8–2. Given the high sequence similarity of many of the DH regions and intergenic regions, we cannot definitively identify exact chromosomal breakpoints and cannot rule out that other events could have occurred in conjunction with the deletion event of the intragenic region between IGHD3 and IGHD7. For example, gene conversion could alternatively have occurred between IGHD6 and IGHD7 paralogs, or a deletion event followed by insertions of repeats into an IGHD7 paralog could have occurred. However, RSS analysis indicates that the 5′ RSS of IGHD8–2 shares identity with either IGHD3 or IGHD6 families (Table 1), thus a fusion between IGHD6 and IGHD7 or gene conversion of IGHD6 into IGHD3 followed by fusion to IGHD7 are likely mechanisms to production of the IGHD8–2 gene through a fusion event. The 3′ RSS of IGHD8–2 is identical to IGHD7 genes, and local alignments show homology between IGHD7 and IGHD8–2, suggesting that a primordial IGHD7 paralog from cluster 2 now forms the 3′ region of IGHD8–2.
DH genes have expanded repeats
Bovine IGHD regions are comprised of multiple repeating short sequence motifs, with the major differences between several DH regions being length differences due to variable numbers of nucleotide repeats (Fig. 2). IGHD7–4 is the second longest DH region, and only differs from IGHD7–3 (its paralog in cluster 3) by one repeat of TGGTTA, which results in a two amino acid insertion. IGHD7–3, IGHD7–4 and IGHD8–2 (the ultralong DH region) are very similar in having several repeating units, but with IGHD8–2 being dramatically longer. The 3′ ends of IGHD7–3 and IGHD7–4 are 85.6 and 77.4% identical to IGHD8–2 over the last 96 nucleotides, respectively (Supplementary Fig. 5). The longer DH regions appear to be evolutionarily active in length evolution based on expanding or contracting repeats, as polymorphisms in Bos taurus IGHD8–2 differ in repeat lengths (Fig. 2, Supplemental Fig. 6). In this regard, two IGHD8–2 polymorphisms have been reported that differ in length and cysteine position, but share similar repeating nucleotide and amino acid sequences [29, 30, 34]. Related species like Bos grunniens (domestic Yak) and Bison bison (American buffalo) also have ultralong CDR H3 regions encoded by IGHD8–2 orthologs, but differ in their lengths due to apparent differences in hexanucleotide repeat expansion within the coding regions (Fig. 2, Supplemental Fig. 6). Thus, while two DH genes may have fused to form the long IGHD8–2 gene, nucleotide repeat expansion or contraction appears to also play a role in long DH region evolution in these species.
The antibody repertoire is a defining evolutionary feature of vertebrates. V(D)J recombination and its associated junctional diversity account for vast potential diversity in antibody receptors on naïve B-cells, with an ability to bind with low affinity to most antigens. Most species utilize many V, D, and J gene segments which produce a great combinatorial potential at the heavy and light chain loci. Some species, however, have fewer functional V, D, and J segments and may use additional mechanisms to add diversity to their repertoires. Cows, in particular, have only 10–12 functional VH regions in only one VH1 family [14, 17, 29], and 16 functional DH regions [29, 35], but have very long CDR H3 regions. The homologous DH regions are cysteine rich, and diversity can be generated through both germline and somatically generated cysteines, which can form a diverse array of potential disulfide bonded loops [23, 35, 36]. In the knob region of ultralong CDR H3s, a diversity of disulfide bond patterns has been observed in several crystal structures [23, 25, 26], and somatic mutations to and from cysteine have been confirmed through deep sequence analysis. Thus, novel cysteines encoded in the DH regions contributes to structural diversity in bovine antibodies.
The length of the DH regions in cows contributes to the overall increase in CDR H3 lengths in the antibody heavy chain repertoire. At the extreme, IGHD8–2 encodes 49 or 51 amino acids, depending on the polymorphic variant [29, 30, 34], and enables CDR H3 lengths of up to 70 amino acids. These CDR H3 regions form independently folding mini-domains comprised of a β-ribbon “stalk” and a disulfide-bonded “knob” that project far from the antigen surface. The sequence diversity of heavy chains with ultralong CDR H3 regions is enormous [16, 23], despite the fact that only one IGHD8–2 region (albeit with two polymorphic variants) is used in this entire class of antibodies. This vast diversity is explained by the fact that cattle utilize AID mediated somatic hypermutation in the pre-immune repertoire, as opposed to after antigen exposure as in other species, and this robust mutation induction substantially diversifies the repertoire through amino acid changes, cysteine mutations to alter disulfide loops, and a substantial proportion of deletional events which can also impact loop structures [16, 21]. All of these diversifying events use the germline IGHD regions as a template during repertoire formation. The IGHD templates feature multiple AID SH “hotspots” as well as nucleotide repeats that preferentially encode Ser, Gly, Tyr, and Cys, often in several repeating units like Gly-Tyr-Gly or Gly-Tyr-Ser. These amino acid repeats could be the result of nucleotide repeat expansion in the genome. Here we show that nucleotide repeating units differ between IGHD paralogs derived from different clusters, and that the unusual ultralong IGHD8–2 region likely formed from a DH-DH fusion in cluster 2 of primordial IGHD6 (or IGHD3) to IGHD7 family members. Clearly a substantial deletion event occurred in the region now encoding IGHD8–2, which can be explained by deletion of the 3′ end of primordial IGHD3, the intergenic region between IGHD3 and IGHD7, and the 5′ end of IGHD7. However, this event was also associated with a more substantial rearrangement that additionally replaced the IGHD3 paralog with IGHD5. Given that homologous variants of IGHD8–2 within Bos taurus, as well as in Bos grunniens and Bison bison differ in the length of IGHD8–2 through differences in the number of short repeats, it is likely that repeat expansion played a role in IGHD8–2 evolution either with a genetic fusion, or potentially with massive expansion in the absence of the fusion of two DH regions.
To summarize, our analysis indicates that the most likely origin of IGHD8–2 is through a fusion event comprising the 5′ end of IGHD6 with the 3′ end of IGHD7 based on sequence similarity analysis, as well as the preservation of the codon encoding the nearly completely conserved first cysteine of the knob domain of ultralong CDR H3 antibodies encoded by IGHD8–2. This event, however, was associated with a larger chromosomal rearrangement that replaced an IGHD3 paralog and its 3′ intergenic region in cluster 2 with IGHD5–2 (Fig. 3). The ability of the genome to diversify IGHD region lengths through genomic rearrangement and repeat expansion through evolutionary time provides a novel genetic mechanism for Darwinian diversification of the vertebrate immune system.
The DNA sequence encoding the bovine antibody heavy chain locus (accession no. KT723008)  was downloaded from the IMGT server (http://www.imgt.org/). Clusters 2, 3 and 4 were defined by the beginning of the Family 1 gene RSS to the end of the Family 6 gene RSS. Sequences of IGHD regions, intergenic regions, and clusters were derived using the Bioconductor program using the R statistical program language . Multiple sequence alignments were done using Clustal Omega, WebPrank, or Muscle (https://www.ebi.ac.uk/services). Local sequence alignments were done using Matcher (https://www.ebi.ac.uk/Tools/psa/). The Bison bison and Bos grunniens IGHD8–2 sequences were identified by BLAST search at the ensembl genome server (www.ensembl.org) using the Bos taurus IGHD8–2 gene as query. Bison bison and Bos grunniens IGHD8–2 genes were found within the genomic sequences with accession numbers XM_010833706.1 (Bison) and CM016710.1 (Yak).
Availability of data and materials
complementarity determining region
heavy chain diversity region
heavy chain variable region
recombination signal sequence
Fugmann SD, Lee AI, Shockett PE, Villey IJ, Schatz DG. The RAG proteins and V(D)J recombination: complexes, ends, and transposition. Annu Rev Immunol. 2000;18:495–527.
Smider V, Chu G. The end-joining reaction in V(D)J recombination. Semin Immunol. 1997;9(3):189–97.
Tonegawa S. Somatic generation of antibody diversity. Nature. 1983;302:575–81.
Burton DR, Hangartner L. Broadly neutralizing antibodies to HIV and their role in vaccine design. Annu Rev Immunol. 2016;34(1):635–59.
Burton DR, Poignard P, Stanfield RL, Wilson IA. Broadly neutralizing antibodies present new prospects to counter highly antigenically diverse viruses. Science. 2012;337(6091):183–6.
Kwong PD, Wilson IA. HIV-1 and influenza antibodies: seeing antigens in new ways. Nat Immunol. 2009;10(6):573–8.
Puligedda RD, Kouiavskaia D, Al-Saleem FH, Kattala CD, Nabi U, Yaqoob H, Bhagavathula VS, Sharma R, Chumakov K, Dessain SK. Characterization of human monoclonal antibodies that neutralize multiple poliovirus serotypes. Vaccine. 2017;35(41):5455–62.
Stanfield R, Wilson I. Antibody structure. Microbiolspec. 2014;2(2). https://doi.org/10.1128/microbiolspec.AID-0012-201310.1186/s12865-020-00359-8.
Douthwaite JA, Sridharan S, Huntington C, Hammersley J, Marwood R, Hakulinen JK, Ek M, Sjogren T, Rider D, Privezentzev C, et al. Affinity maturation of a novel antagonistic human monoclonal antibody with a long VH CDR3 targeting the class A GPCR formyl-peptide receptor 1. MAbs. 2015;7(1):152–66.
Nam DH, Rodriguez C, Remacle AG, Strongin AY, Ge X. Active-site MMP-selective antibody inhibitors discovered from convex paratope synthetic libraries. Proc Natl Acad Sci U S A. 2016;113(52):14970–5.
Bonsignori M, Hwang KK, Chen X, Tsao CY, Morris L, Gray E, Marshall DJ, Crump JA, Kapiga SH, Sam NE, et al. Analysis of a clonal lineage of HIV-1 envelope V2/V3 conformational epitope-specific broadly neutralizing antibodies and their inferred unmutated common ancestors. J Virol. 2011;85(19):9998–10009.
Doria-Rose NA, Schramm CA, Gorman J, Moore PL, Bhiman JN, DeKosky BJ, Ernandes MJ, Georgiev IS, Kim HJ, Pancera M, et al. Developmental pathway for potent V1V2-directed HIV-neutralizing antibodies. Nature. 2014;509(7498):55–62.
Walker LM, Huber M, Doores KJ, Falkowska E, Pejchal R, Julien JP, Wang SK, Ramos A, Chan-Hui PY, Moyle M, et al. Broad neutralization coverage of HIV by multiple highly potent antibodies. Nature. 2011;477(7365):466–70.
Berens SJ, Wylie DE, Lopez OJ. Use of a single VH family and long CDR3s in the variable region of cattle Ig heavy chains. Int Immunol. 1997;9(1):189–99.
de los Rios M, Criscitiello MF, Smider VV. Structural and genetic diversity in antibody repertoires from diverse species. Curr Opin Struct Biol. 2015;33:27–41.
Deiss TC, Vadnais M, Wang F, Chen PL, Torkamani A, Mwangi W, Lefranc M-P, Criscitiello MF, Smider VV. Immunogenetic factors driving formation of ultralong VH CDR3 in Bos taurus antibodies. Cell Mol Immunol. 2017;14:1–12.
Lopez O, Perez C, Wylie D. A single VH family and long CDR3s are the targets for hypermutation in bovine immunoglobulin heavy chains. Immunol Rev. 1998;162:55–66.
Saini SS, Allore B, Jacobs RM, Kaushik A. Exceptionally long CDR3H region with multiple cysteine residues in functional bovine IgM antibodies. Eur J Immunol. 1999;29(8):2420–6.
Saini SS, Farrugia W, Ramsland PA, Kaushik AK. Bovine IgM antibodies with exceptionally long complementarity-determining region 3 of the heavy chain share unique structural properties conferring restricted VH + Vlambda pairings. Int Immunol. 2003;15(7):845–53.
Saini SS, Kaushik A. Extensive CDR3H length heterogeneity exists in bovine foetal VDJ rearrangements. Scand J Immunol. 2002;55(2):140–8.
Stanfield RL, Haakenson J, Deiss TC, Criscitiello MF, Wilson IA, Smider VV. The unusual genetics and biochemistry of bovine Immunoglobulins. Adv Immunol. 2018;137:135–64.
Walther S, Czerny C-P, Diesterbeck US. Exceptionally long CDR3H are not isotype restricted in bovine immunoglobulins. PLoS One. 2013;8(5):e64234.
Wang F, Ekiert DC, Ahmad I, Yu W, Zhang Y, Bazirgan O, Torkamani A, Raudsepp T, Mwangi W, Criscitiello MF, et al. Reshaping antibody diversity. Cell. 2013;153(6):1379–93.
Zhao Y, Jackson SM, Aitken R. The bovine antibody repertoire. Dev Comp Immunol. 2006;30(1–2):175–86.
Dong J, Finn JA, Larsen PA, Smith TPL, Crowe JE Jr. Structural diversity of Ultralong CDRH3s in seven bovine antibody heavy chains. Front Immunol. 2019;10:558.
Stanfield RL, Wilson IA, Smider VV. Conservation and diversity in the ultralong third heavy-chain complementarity-determining region of bovine antibodies. Sci Immunol. 2016;1(1):aaf7962.
Sok D, Le KM, Vadnais M, Saye-Francisco KL, Jardine JG, Torres JL, Berndsen ZT, Kong L, Stanfield R, Ruiz J, et al. Rapid elicitation of broadly neutralizing antibodies to HIV by immunization in cows. Nature. 2017;548(7665):108–11.
Koti M, Kataeva G, Kaushik A. Organization of DH-gene locus is distinct in cattle. Dev Biol. 2008;132:307–13.
Ma L, Qin T, Chu D, Cheng X, Wang J, Wang X, Wang P, Han H, Ren L, Aitken R, et al. Internal duplications of DH, JH, and C region genes create an unusual IgH gene locus in cattle. J Immunol. 2016;196(10):4358–66.
Shojaei F, Saini SS, Kaushik AK. Unusually long germline DH genes contribute to large sized CDR3H in bovine antibodies. Mol Immunol. 2003;40(1):61–7.
Lefranc MP. Immunoglobulin and T cell receptor genes: IMGT((R)) and the birth and rise of Immunoinformatics. Front Immunol. 2014;5:22.
Lefranc MP, Pommie C, Ruiz M, Giudicelli V, Foulquier E, Truong L, Thouvenin-Contet V, Lefranc G. IMGT unique numbering for immunoglobulin and T cell receptor variable domains and Ig superfamily V-like domains. Dev Comp Immunol. 2003;27(1):55–77.
Lefranc M-P, Lefranc G. The immunoglobulin FactsBook. London: Academic Press; 2001.
Liljavirta J, Niku M, Pessa-Morikawa T, Ekman A, Iivanainen A. Expansion of the preimmune antibody repertoire by junctional diversity in Bos taurus. PLoS One. 2014;9(6):e99808.
Haakenson JK, Deiss TC, Warner GF, Mwangi W, Criscitiello MF, Smider VV. A broad role for Cysteines in bovine antibody diversity. Immunohorizons. 2019;3(10):478–87.
Haakenson JK, Huang R, Smider VV. Diversity in the cow Ultralong CDR H3 antibody repertoire. Front Immunol. 2018;9:1262.
Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, et al. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004;5(10):R80.
We are grateful for helpful conversations regarding this work with Ali Torkamani and Michael Criscitiello.
This work was funded by NIH grants R01 GM105826 and R01 HD088400 to V.V.S and funded salary support and supplies.
Consent for publication
The authors declare they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Alignment of Bos taurus DH regions.
Alignment of DH clusters 2, 3, and 4.
Phylogenetic of DH paralogs from clusters 2 and 3.
Phylogenetic analysis of intergenic regions from clusters 2 and 3.
Alignment of IGHD7 family members with IGHD8–2.
Identification ultralong IGHD region homologs in Bos taurus, Bos grunniens and Bison bison.
Alignment of IGHD3 and IGHD6 family members with IGHD8–2.
Fusion analysis of IGHD members to form IGHD8–2 homologs.
About this article
Cite this article
Smider, B.A., Smider, V.V. Formation of ultralong DH regions through genomic rearrangement. BMC Immunol 21, 30 (2020). https://doi.org/10.1186/s12865-020-00359-8
- Bovine immunoglobulin
- DH region
- Ultralong CDR3
- Cow antibody
- Antibody diversity
- Immunoglobulin repertoire
- CDR3 length