PCR amplification of murine immunoglobulin germline V genes: strategies for minimization of recombination artefacts

Somatic mutation patterns at Ig and Non-Ig Loci

The reverse transcriptase (RT) model of immunoglobulin (Ig) somatic hypermutation (SHM) has received insufficient scientific attention. This is understandable given that DNA deamination mediated by activation-induced deaminase (AID), the initiating step of Ig SHM, has dominated experiments since 2002. We summarise some key history of the RT Ig SHM model dating to 1987. For example, it is now established that DNA polymerase η, the sole DNA repair polymerase involved in post-replication short-patch repair, is an efficient cellular RT. This implies that it is potentially able to initiate target site reverse transcription by RNA-directed DNA repair at AID-induced lesions. Recently, DNA polymerase θ has also been shown to be an efficient cellular RT. Since DNA polymerase θ plays no significant role in Ig SHM, it could serve a similar RNA-dependent DNA polymerase role as DNA polymerase η at non-Ig loci in the putative RNA-templated nucleotide excision repair of bulky adducts and other mutagenic lesions on the transcribed strand. A major yet still poorly recognised consequence of the proposed RT process in Ig SHM is the generation of significant and characteristic strand-biased mutation signatures at both deoxyadenosine/deoxythymidine and deoxyguanosine/deoxycytidine base pairs. In this historical perspective, we highlight how diagnostic strand-biased mutation signatures are detected in vivo during SHM at both Ig loci in germinal centre B lymphocytes and non-Ig loci in cancer genomes. These strand-biased signatures have been significantly obscured by technical issues created by improper use of the polymerase chain reaction technique. A heightened awareness of this fact should contribute to better data interpretation and somatic mutation pattern recognition both at Ig and non-Ig loci.

A BALB/c IGHV Reference Set, Defined by Haplotype Analysis of Long-Read VDJ-C Sequences From F1 (BALB/c x C57BL/6) Mice

The immunoglobulin genes of inbred mouse strains that are commonly used in models of antibody-mediated human diseases are poorly characterized. This compromises data analysis. To infer the immunoglobulin genes of BALB/c mice, we used long-read SMRT sequencing to amplify VDJ-C sequences from F1 (BALB/c x C57BL/6) hybrid animals. Strain variations were identified in the Ighm and Ighg2b genes, and analysis of VDJ rearrangements led to the inference of 278 germline IGHV alleles. 169 alleles are not present in the C57BL/6 genome reference sequence. To establish a set of expressed BALB/c IGHV germline gene sequences, we computationally retrieved IGHV haplotypes from the IgM dataset. Haplotyping led to the confirmation of 162 BALB/c IGHV gene sequences. A musIGHV398 pseudogene variant also appears to be present in the BALB/cByJ substrain, while a functional musIGHV398 gene is highly expressed in the BALB/cJ substrain. Only four of the BALB/c alleles were also observed in the C57BL/6 haplotype. The full set of inferred BALB/c sequences has been used to establish a BALB/c IGHV reference set, hosted at https://ogrdb.airr-community.org. We assessed whether assemblies from the Mouse Genome Project (MGP) are suitable for the determination of the genes of the IGH loci. Only 37 (43.5%) of the 85 confirmed IMGT-named BALB/c IGHV and 33 (42.9%) of the 77 confirmed non-IMGT IGHV were found in a search of the MGP BALB/cJ genome assembly. This suggests that current MGP assemblies are unsuitable for the comprehensive documentation of germline IGHVs and more efforts will be needed to establish strain-specific reference sets.

Inferred Allelic Variants of Immunoglobulin Receptor Genes: A System for Their Evaluation, Documentation, and Naming

Immunoglobulins or antibodies are the main effector molecules of the B-cell lineage and are encoded by hundreds of variable (V), diversity (D), and joining (J) germline genes, which recombine to generate enormous IG diversity. Recently, high-throughput adaptive immune receptor repertoire sequencing (AIRR-seq) of recombined V-(D)-J genes has offered unprecedented insights into the dynamics of IG repertoires in health and disease. Faithful biological interpretation of AIRR-seq studies depends upon the annotation of raw AIRR-seq data, using reference germline gene databases to identify the germline genes within each rearrangement. Existing reference databases are incomplete, as shown by recent AIRR-seq studies that have inferred the existence of many previously unreported polymorphisms. Completing the documentation of genetic variation in germline gene databases is therefore of crucial importance. Lymphocyte receptor genes and alleles are currently assigned by the Immunoglobulins, T cell Receptors and Major Histocompatibility Nomenclature Subcommittee of the International Union of Immunological Societies (IUIS) and managed in IMGT®, the international ImMunoGeneTics information system® (IMGT). In 2017, the IMGT Group reached agreement with a group of AIRR-seq researchers on the principles of a streamlined process for identifying and naming inferred allelic sequences, for their incorporation into IMGT®. These researchers represented the AIRR Community, a network of over 300 researchers whose objective is to promote all aspects of immunoglobulin and T-cell receptor repertoire studies, including the standardization of experimental and computational aspects of AIRR-seq data generation and analysis. The Inferred Allele Review Committee (IARC) was established by the AIRR Community to devise policies, criteria, and procedures to perform this function. Formalized evaluations of novel inferred sequences have now begun and submissions are invited via a new dedicated portal (https://ogrdb.airr-community.org). Here, we summarize recommendations developed by the IARC-focusing, to begin with, on human IGHV genes-with the goal of facilitating the acceptance of inferred allelic variants of germline IGHV genes. We believe that this initiative will improve the quality of AIRR-seq studies by facilitating the description of human IG germline gene variation, and that in time, it will expand to the documentation of TR and IG genes in many vertebrate species.

Understanding and overcoming the pitfalls and biases of next-generation sequencing (NGS) methods for use in the routine clinical microbiological diagnostic laboratory

Recent advancements in next-generation sequencing (NGS) have provided the foundation for modern studies into the composition of microbial communities. The use of these NGS methods allows for the detection and identification of (‘difficult-to-culture’) microorganisms using a culture-independent strategy. In the field of routine clinical diagnostics however, the application of NGS is currently limited to microbial strain typing for epidemiological purposes only, even though the implementation of NGS for microbial community analysis may yield clinically important information. This lack of NGS implementation is due to many different factors, including issues relating to NGS method standardization and result reproducibility. In this review article, the authors provide a general introduction to the most widely used NGS methods currently available (i.e., targeted amplicon sequencing and shotgun metagenomics) and the strengths and weaknesses of each method is discussed. The focus of the publication then shifts toward 16S rRNA gene NGS methods, which are currently the most cost-effective and widely used NGS methods for research purposes, and are therefore more likely to be successfully implemented into routine clinical diagnostics in the short term. In this respect, the experimental pitfalls and biases created at each step of the 16S rRNA gene NGS workflow are explained, as well as their potential solutions. Finally, a novel diagnostic microbiota profiling platform (‘MYcrobiota’) is introduced, which was developed by the authors by taking into consideration the pitfalls, biases, and solutions explained in this article. The development of the MYcrobiota, and future NGS methodologies, will help pave the way toward the successful implementation of NGS methodologies into routine clinical diagnostics.

Uncoupling of sgRNAs from their associated barcodes during PCR amplification of combinatorial CRISPR screens

Many implementations of pooled screens in mammalian cells rely on linking an element of interest to a barcode, with the latter subsequently quantitated by next generation sequencing. However, substantial uncoupling between these paired elements during lentiviral production has been reported, especially as the distance between elements increases. We detail that PCR amplification is another major source of uncoupling, and becomes more pronounced with increased amounts of DNA template molecules and PCR cycles. To lessen uncoupling in systems that use paired elements for detection, we recommend minimizing the distance between elements, using low and equal template DNA inputs for plasmid and genomic DNA during PCR, and minimizing the number of PCR cycles. We also present a vector design for conducting combinatorial CRISPR screens that enables accurate barcode-based detection with a single short sequencing read and minimal uncoupling.

Reciprocal translocation of small numbers of inbred individuals rescues immunogenetic diversity

Genetic rescue can reduce inbreeding depression and increase fitness of small populations, even when the donor populations are highly inbred. In a recent experiment involving two inbred island populations of the New Zealand South Island robin, Petroica australis, reciprocal translocations improved microsatellite diversity and individual fitness. While microsatellite loci may reflect patterns of genome-wide diversity, they generally do not indicate the specific genetic regions responsible for increased fitness. We tested the effectiveness of this reciprocal translocation for rescuing diversity of two immunogenetic regions: Toll-like receptor (TLR) and major histocompatibility complex (MHC) genes. We found that the relatively small number of migrants (seven and ten per island) effectively brought the characteristic TLR gene diversity of each source population into the recipient population. However, when migrants transmitted TLR alleles that were already present at high frequency in the recipient population, it was possible for offspring of mixed heritage to have decreased gene diversity compared to recipient population diversity prior to translocation. In contrast to TLRs, we did not observe substantial changes in MHC allelic diversity following translocation, with limited evidence of a decrease in differentiation, perhaps because most MHC alleles were observed at both sites prior to the translocation. Overall, we conclude that small numbers of migrants may successfully restore the diversity of immunogenetic loci with few alleles, but that translocating larger numbers of animals would provide additional opportunity for the genetic rescue of highly polymorphic immunity regions, such as the MHC, even when the source population is inbred.

Long-Read Single Molecule Sequencing to Resolve Tandem Gene Copies: The Mst77Y Region on the Drosophila melanogaster Y Chromosome

The autosomal gene Mst77F of Drosophila melanogaster is essential for male fertility. In 2010, Krsticevic et al. (Genetics184: 295−307) found 18 Y-linked copies of Mst77F (“Mst77Y”), which collectively account for 20% of the functional Mst77F-like mRNA. The Mst77Y genes were severely misassembled in the then-available genome assembly and were identified by cloning and sequencing polymerase chain reaction products. The genomic structure of the Mst77Y region and the possible existence of additional copies remained unknown. The recent publication of two long-read assemblies of D. melanogaster prompted us to reinvestigate this challenging region of the Y chromosome. We found that the Illumina Synthetic Long Reads assembly failed in the Mst77Y region, most likely because of its tandem duplication structure. The PacBio MHAP assembly of the Mst77Y region seems to be very accurate, as revealed by comparisons with the previously found Mst77Y genes, a bacterial artificial chromosome sequence, and Illumina reads of the same strain. We found that the Mst77Y region spans 96 kb and originated from a 3.4-kb transposition from chromosome 3L to the Y chromosome, followed by tandem duplications inside the Y chromosome and invasion of transposable elements, which account for 48% of its length. Twelve of the 18 Mst77Y genes found in 2010 were confirmed in the PacBio assembly, the remaining six being polymerase chain reaction−induced artifacts. There are several identical copies of some Mst77Y genes, coincidentally bringing the total copy number to 18. Besides providing a detailed picture of the Mst77Y region, our results highlight the utility of PacBio technology in assembling difficult genomic regions such as tandemly repeated genes.

Soma-to-germline feedback is implied by the extreme polymorphism at IGHV relative to MHC: The manifest polymorphism of the MHC appears greatly exceeded at Immunoglobulin loci, suggesting antigen-selected somatic V mutants penetrate Weismann's Barrier

Soma-to-germline feedback is forbidden under the neo-Darwinian paradigm. Nevertheless, there is a growing realization it occurs frequently in immunoglobulin (Ig) variable (V) region genes. This is a surprising development. It arises from a most unlikely source in light of the exposure of co-author EJS to the haplotype data of RL Dawkins and others on the polymorphism of the Major Histocompatibility Complex, which is generally assumed to be the most polymorphic region in the genome (spanning ∼4 Mb). The comparison between the magnitude of MHC polymorphism with estimates for the human heavy chain immunoglobulin V locus (spanning ∼1 Mb), suggests IGHV could be many orders of magnitude more polymorphic than the MHC. This conclusion needs airing in the literature as it implies generational churn and soma-to-germline gene feedback. Pedigree-based experimental strategies to resolve the IGHV issue are outlined.

Effect of the enzyme and PCR conditions on the quality of high-throughput DNA sequencing results

Library preparation protocols for high-throughput DNA sequencing (HTS) include amplification steps in which errors can build up. In order to have confidence in the sequencing data, it is important to understand the effects of different Taq polymerases and PCR amplification protocols on the DNA molecules sequenced. We compared thirteen enzymes in three different marker systems: simple, single copy nuclear gene and complex multi-gene family. We also tested a modified PCR protocol, which has been suggested to reduce errors associated with amplification steps. We find that enzyme choice has a large impact on the proportion of correct sequences recovered. The most complex marker systems yielded fewer correct reads, and the proportion of correct reads was greatly affected by the enzyme used. Modified cycling conditions did reduce the number of incorrect sequences obtained in some cases, but enzyme had a much greater impact on the number of correct reads. Thus, the coverage required for the safe identification of genotypes using one of the low quality enzymes could be seven times larger than with more efficient enzymes in a biallelic system with equal amplification of the two alleles. Consequently, enzyme selection for downstream HTS has important consequences, especially in complex genetic systems.

Multi-template polymerase chain reaction

PCR is a formidable and potent technology that serves as an indispensable tool in a wide range of biological disciplines. However, due to the ease of use and often lack of rigorous standards many PCR applications can lead to highly variable, inaccurate, and ultimately meaningless results. Thus, rigorous method validation must precede its broad adoption to any new application. Multi-template samples possess particular features, which make their PCR analysis prone to artifacts and biases: multiple homologous templates present in copy numbers that vary within several orders of magnitude. Such conditions are a breeding ground for chimeras and heteroduplexes. Differences in template amplification efficiencies and template competition for reaction compounds undermine correct preservation of the original template ratio. In addition, the presence of inhibitors aggravates all of the above-mentioned problems. Inhibitors might also have ambivalent effects on the different templates within the same sample. Yet, no standard approaches exist for monitoring inhibitory effects in multitemplate PCR, which is crucial for establishing compatibility between samples.

External and semi-internal controls for PCR amplification of homologous sequences in mixed templates

In a mixed template, the presence of homologous target DNA sequences creates environments that almost inevitably give rise to artifacts and biases during PCR. Heteroduplexes, chimeras, and skewed template-to-product ratios are the exclusive attributes of mixed template PCR and never occur in a single template assay. Yet, multi-template PCR has been used without appropriate attention to quality control and assay validation, in spite of the fact that such practice diminishes the reliability of results. External and internal amplification controls became obligatory elements of good laboratory practice in different PCR assays. We propose the inclusion of an analogous approach as a quality control system for multi-template PCR applications. The amplification controls must take into account the characteristics of multi-template PCR and be able to effectively monitor particular assay performance. This study demonstrated the efficiency of a model mixed template as an adequate external amplification control for a particular PCR application. The conditions of multi-template PCR do not allow implementation of a classic internal control; therefore we developed a convenient semi-internal control as an acceptable alternative. In order to evaluate the effects of inhibitors, a model multi-template mix was amplified in a mixture with DNAse-treated sample. Semi-internal control allowed establishment of intervals for robust PCR performance for different samples, thus enabling correct comparison of the samples. The complexity of the external and semi-internal amplification controls must be comparable with the assumed complexity of the samples. We also emphasize that amplification controls should be applied in multi-template PCR regardless of the post-assay method used to analyze products.

A novel low temperature PCR assured high-fidelity DNA amplification

As previously reported, a novel low temperature (LoTemp) polymerase chain reaction (PCR) catalyzed by a moderately heat-resistant (MHR) DNA polymerase with a chemical-assisted denaturation temperature set at 85 °C instead of the conventional 94-96 °C can achieve high-fidelity DNA amplification of a target DNA, even after up to 120 PCR thermal cycles. Furthermore, such accurate amplification is not achievable with conventional PCR. Now, using a well-recognized L1 gene segment of the human papillomavirus (HPV) type 52 (HPV-52) as the template for experiments, we demonstrate that the LoTemp high-fidelity DNA amplification is attributed to an unusually high processivity and stability of the MHR DNA polymerase whose high fidelity in template-directed DNA synthesis is independent of non-existent 3'-5' exonuclease activity. Further studies and understanding of the characteristics of the LoTemp PCR technology may facilitate implementation of DNA sequencing-based diagnostics at the point of care in community hospital laboratories.

Molecular analysis of immunoglobulin genes reveals frequent clonal relatedness in double monoclonal gammopathies

Monoclonal gammopathies (MGs) are hematological diseases characterized by high levels of a monoclonal immunoglobulin (Ig) or M-protein. Within this group are patients with more than one M-protein, referred to as double MGs (DMGs). The M-proteins in DMG patients may have different heavy chain (HC) isotypes that are associated with different light chains (LCs), or different HCs that are LC matched. In this study, we examined the clonal relatedness of the M-proteins in the latter type in a cohort of 14 DMG patients. By using PCR, we identified 7/14 DMG patients that expressed two Ig HC isotypes with identical Ig HC variable (IGHV), diversity (IGHD), joining (IGHJ), and complementarity determining region (HCDR3) sequences. Two additional DMG patients had two Ig transcripts using the same IGHV, IGHD and IGHJ genes but with slight differences in variable region or HCDR3 mutations. LC analysis confirmed that a single LC was expressed in 3/7 DMG patients with identical HC transcripts and in the two DMGs with highly similar transcripts. The PCR findings were confirmed by immunofluorescence for HC and LC expression. Clonally related HC-dissimilar/LC-matched DMGs may occur often and defines a new subtype of MG that may serve as a tool for studies of disease pathogenesis.

Critical Analysis of Strand-Biased Somatic Mutation Signatures in TP53 versus Ig Genes, in Genome-Wide Data and the Etiology of Cancer

Previous analyses of rearranged immunoglobulin (Ig) variable genes (VDJs) concluded that the mechanism of Ig somatic hypermutation (SHM) involves the Ig pre-mRNA acting as a copying template resulting in characteristic strand biased somatic mutation patterns at A:T and G:C base pairs. We have since analysed cancer genome data and found the same mutation strand-biases, in toto or in part, in nonlymphoid cancers. Here we have analysed somatic mutations in a single well-characterised gene TP53. Our goal is to understand the genesis of the strand-biased mutation patterns in TP53—and in genome-wide data—that may arise by “endogenous” mechanisms as opposed to adduct-generated DNA-targeted strand-biased mutations caused by well-characterised “external” carcinogenic influences in cigarette smoke, UV-light, and certain dietary components. The underlying strand-biased mutation signatures in TP53, for many non-lymphoid cancers, bear a striking resemblance to the Ig SHM pattern. A similar pattern can be found in genome-wide somatic mutations in cancer genomes that have also mutated TP53. The analysis implies a role for base-modified RNA template intermediates coupled to reverse transcription in the genesis of many cancers. Thus Ig SHM may be inappropriately activated in many non-lymphoid tissues via hormonal and/or inflammation-related processes leading to cancer.

Effects of PCR cycle number and DNA polymerase type on the 16S rRNA gene pyrosequencing analysis of bacterial communities

The effects of PCR cycle number and DNA polymerase type on 16S rRNA gene pyrosequencing analysis were investigated using an artificially prepared bacterial community (mock community). The bacterial richness was overestimated at increased PCR cycle number mostly due to the occurence of chimeric sequences, and this was more serious with a DNA polymerase having proofreading activity than with Taq DNA polymerase. These results suggest that PCR cycle number must be kept as low as possible for accurate estimation of bacterial richness and that particular care must be taken when a DNA polymerase having proofreading activity is used.

Comprehensive assessment of potential multiple myeloma immunoglobulin heavy chain V-D-J intraclonal variation using massively parallel pyrosequencing

Multiple myeloma (MM) is characterized by the accumulation of malignant plasma cells (PCs) in the bone marrow (BM). MM is viewed as a clonal disorder due to lack of verified intraclonal sequence diversity in the immunoglobulin heavy chain variable region gene (IGHV). However, this conclusion is based on analysis of a very limited number of IGHV subclones and the methodology employed did not permit simultaneous analysis of the IGHV repertoire of non-malignant PCs in the same samples. Here we generated genomic DNA and cDNA libraries from purified MM BMPCs and performed massively parallel pyrosequencing to determine the frequency of cells expressing identical IGHV sequences. This method provided an unprecedented opportunity to interrogate the presence of clonally related MM cells and evaluate the IGHV repertoire of non-MM PCs. Within the MM sample, 37 IGHV genes were expressed, with 98.9% of all immunoglobulin sequences using the same IGHV gene as the MM clone and 83.0% exhibiting exact nucleotide sequence identity in the IGHV and heavy chain complementarity determining region 3 (HCDR3). Of interest, we observed in both genomic DNA and cDNA libraries 48 sets of identical sequences with single point mutations in the MM clonal IGHV or HCDR3 regions. These nucleotide changes were suggestive of putative subclones and therefore were subjected to detailed analysis to interpret: 1) their legitimacy as true subclones; and 2) their significance in the context of MM. Finally, we report for the first time the IGHV repertoire of normal human BMPCs and our data demonstrate the extent of IGHV repertoire diversity as well as the frequency of clonally-related normal BMPCs. This study demonstrates the power and potential weaknesses of in-depth sequencing as a tool to thoroughly investigate the phylogeny of malignant PCs in MM and the IGHV repertoire of normal BMPCs.

High diversity at PRDM9 in chimpanzees and bonobos

Background

The PRDM9 locus in mammals has increasingly attracted research attention due to its role in mediating chromosomal recombination and possible involvement in hybrid sterility and hence speciation processes. The aim of this study was to characterize sequence variation at the PRDM9 locus in a sample of our closest living relatives, the chimpanzees and bonobos.

Methodology/Principal Findings

PRDM9 contains a highly variable and repetitive zinc finger array. We amplified this domain using long-range PCR and determined the DNA sequences using conventional Sanger sequencing. From 17 chimpanzees representing three subspecies and five bonobos we obtained a total of 12 alleles differing at the nucleotide level. Based on a data set consisting of our data and recently published Pan PRDM9 sequences, we found that at the subspecies level, diversity levels did not differ among chimpanzee subspecies or between chimpanzee subspecies and bonobos. In contrast, the sample of chimpanzees harbors significantly more diversity at PRDM9 than samples of humans. Pan PRDM9 shows signs of rapid evolution including no alleles or ZnFs in common with humans as well as signals of positive selection in the residues responsible for DNA binding.

Conclusions and Significance

The high number of alleles specific to the genus Pan, signs of positive selection in the DNA binding residues, and reported lack of conservation of recombination hotspots between chimpanzees and humans suggest that PRDM9 could be active in hotspot recruitment in the genus Pan. Chimpanzees and bonobos are considered separate species and do not have overlapping ranges in the wild, making the presence of shared alleles at the amino acid level between the chimpanzee and bonobo species interesting in view of the hypothesis that PRDM9 plays a universal role in interspecific hybrid sterility.

Massively parallel haplotyping on microscopic beads for the high-throughput phase analysis of single molecules

In spite of the many advances in haplotyping methods, it is still very difficult to characterize rare haplotypes in tissues and different environmental samples or to accurately assess the haplotype diversity in large mixtures. This would require a haplotyping method capable of analyzing the phase of single molecules with an unprecedented throughput. Here we describe such a haplotyping method capable of analyzing in parallel hundreds of thousands single molecules in one experiment. In this method, multiple PCR reactions amplify different polymorphic regions of a single DNA molecule on a magnetic bead compartmentalized in an emulsion drop. The allelic states of the amplified polymorphisms are identified with fluorescently labeled probes that are then decoded from images taken of the arrayed beads by a microscope. This method can evaluate the phase of up to 3 polymorphisms separated by up to 5 kilobases in hundreds of thousands single molecules. We tested the sensitivity of the method by measuring the number of mutant haplotypes synthesized by four different commercially available enzymes: Phusion, Platinum Taq, Titanium Taq, and Phire. The digital nature of the method makes it highly sensitive to detecting haplotype ratios of less than 1:10,000. We also accurately quantified chimera formation during the exponential phase of PCR by different DNA polymerases.

Evaluation of two approaches to genotyping major histocompatibility complex class I in a passerine-CE-SSCP and 454 pyrosequencing

Genes of the highly dynamic major histocompatibility complex (MHC) are directly linked to individual fitness and are of high interest in evolutionary ecology and conservation genetics. Gene duplication and positive selection usually lead to high levels of polymorphism in the MHC region, making genotyping of MHC a challenging task. Here, we compare the performance of two methods for MHC class I genotyping in a passerine with highly duplicated MHC class I genes: capillary electrophoresis-single-strand conformation polymorphism (CE-SSCP) analysis and 454 GS FLX Titanium pyrosequencing. According to our findings, the number of MHC variants (called alleles for simplicity) detected by CE-SSCP is significantly lower than detected by 454. To resolve discrepancies between the two methods, we cloned and Sanger sequenced a MHC class I amplicon for an individual with high number of alleles. We found a perfect congruence between cloning/Sanger sequencing results and 454. Thus, in case of multi-locus amplification, CE-SSCP considerably underestimates individual MHC diversity. However, numbers of alleles detected by both methods are significantly correlated, although the correlation is weak (r = 0.32). Thus, in systems with highly duplicated MHC, 454 provides more reliable information on individual diversity than CE-SSCP.

Polymorphism and Balancing Selection of MHC Class II DAB Gene in 7 Selective Flounder (Paralichthys olivaceus) Families

In order to determine the genetic variation of the MHC class IIB exon2 allele in the offspring, 700 fry from seven families of Japanese flounder challenged with V. anguillarum were studied, and different mortality rates were found in those families. Five to ten surviving and dead fry from each of the seven families were selected to study the MHC class II B exon2 gene with PCR and a direct sequencing method. One hundred and sixteen different exon2 sequences were found and 116 different alleles were identified, while a minimum of four loci were revealed in the MHC class II B exon2 gene. The ratio (d_N/d_S) of nonsynonymous substitution (d_N) to synonymous substitutions (d_S) in the peptide-binding region (PBR) of the MHC class IIB gene was 6.234, which indicated that balancing selection is acting on the MHC class IIB genes. The MHC IIB alleles were thus being passed on to their progeny. Some alleles were significantly more frequent in surviving than dead individuals. All together our data suggested that the alleles Paol-DAB*4301, Paol-DAB*4601, Paol-DAB*4302, Paol-DAB*3803, and Paol-DAB*4101 were associated with resistance to V. anguillarum in flounder.

To be or not to be the odd one out--allele-specific transcription in pentaploid dogroses (Rosa L. sect. Caninae (DC.) Ser)

BACKGROUND

Multiple hybridization events gave rise to pentaploid dogroses which can reproduce sexually despite their uneven ploidy level by the unique canina meiosis. Two homologous chromosome sets are involved in bivalent formation and are transmitted by the haploid pollen grains and the tetraploid egg cells. In addition the egg cells contain three sets of univalent chromosomes which are excluded from recombination. In this study we investigated whether differential behavior of chromosomes as bivalents or univalents is reflected by sequence divergence or transcription intensity between homeologous alleles of two single copy genes (LEAFY, cGAPDH) and one ribosomal DNA locus (nrITS).

RESULTS

We detected a maximum number of four different alleles of all investigated loci in pentaploid dogroses and identified the respective allele with two copies, which is presumably located on bivalent forming chromosomes. For the alleles of the ribosomal DNA locus and cGAPDH only slight, if any, differential transcription was determined, whereas the LEAFY alleles with one copy were found to be significantly stronger expressed than the LEAFY allele with two copies. Moreover, we found for the three marker genes that all alleles have been under similar regimes of purifying selection.

CONCLUSIONS

Analyses of both molecular sequence evolution and expression patterns did not support the hypothesis that unique alleles probably located on non-recombining chromosomes are less functional than duplicate alleles presumably located on recombining chromosomes.

Estimating genetic variability in non-model taxa: a general procedure for discriminating sequence errors from actual variation

Genetic variation is the driving force of evolution and as such is of central interest for biologists. However, inadequate discrimination of errors from true genetic variation could lead to incorrect estimates of gene copy number, population genetic parameters, phylogenetic relationships and the deposition of gene and protein sequences in databases that are not actually present in any organism. Misincorporation errors in multi-template PCR cloning methods, still commonly used for obtaining novel gene sequences in non-model species, are difficult to detect, as no previous information may be available about the number of expected copies of genes belonging to multi-gene families. However, studies employing these techniques rarely describe in any great detail how errors arising in the amplification process were detected and accounted for. Here, we estimated the rate of base misincorporation of a widely-used PCR-cloning method, using a single copy mitochondrial gene from a single individual to minimise variation in the template DNA, as 1.62×10(-3) errors per site, or 9.26×10(-5) per site per duplication. The distribution of errors among sequences closely matched that predicted by a binomial distribution function. The empirically estimated error rate was applied to data, obtained using the same methods, from the Phospholipase A(2) toxin family from the pitviper Ovophis monticola. The distribution of differences detected closely matched the expected distribution of errors and we conclude that, when undertaking gene discovery or assessment of genetic diversity using this error-prone method, it will be informative to empirically determine the rate of base misincorporation.

Reducing chimera formation during PCR amplification to ensure accurate genotyping

Measurements of population diversity are fundamental to the reconstruction of the evolutionary and epidemiological history of organisms. Commonly used protocols to measure population diversity using the polymerase chain reaction (PCR) are prone to the introduction of artificial chimeras. These are often difficult to detect and can confound the correct interpretation of results due to the false generation of recombinants when the underlying DNA sample contains multiple distinct templates. This study presents a standardised procedure to suppress the formation of artificial chimeras during PCR amplification. The solution is based on the accurate determination of the efficiency and end point of the log-linear phase of a PCR. This procedure will facilitate the generation of data sets that more accurately reflect the underlying population diversity rather than artifacts introduced by the process itself.

Intra-isolate genome variation in arbuscular mycorrhizal fungi persists in the transcriptome

Arbuscular mycorrhizal fungi (AMF) are heterokaryotes with an unusual genetic makeup. Substantial genetic variation occurs among nuclei within a single mycelium or isolate. AMF reproduce through spores that contain varying fractions of this heterogeneous population of nuclei. It is not clear whether this genetic variation on the genome level actually contributes to the AMF phenotype. To investigate the extent to which polymorphisms in nuclear genes are transcribed, we analysed the intra-isolate genomic and cDNA sequence variation of two genes, the large subunit ribosomal RNA (LSU rDNA) of Glomus sp. DAOM-197198 (previously known as G. intraradices) and the POL1-like sequence (PLS) of Glomus etunicatum. For both genes, we find high sequence variation at the genome and transcriptome level. Reconstruction of LSU rDNA secondary structure shows that all variants are functional. Patterns of PLS sequence polymorphism indicate that there is one functional gene copy, PLS2, which is preferentially transcribed, and one gene copy, PLS1, which is a pseudogene. This is the first study that investigates AMF intra-isolate variation at the transcriptome level. In conclusion, it is possible that, in AMF, multiple nuclear genomes contribute to a single phenotype.

Functional copies of the Mst77F gene on the Y chromosome of Drosophila melanogaster

The Y chromosome of Drosophila melanogaster has <20 protein-coding genes. These genes originated from the duplication of autosomal genes and have male-related functions. In 1993, Russell and Kaiser found three Y-linked pseudogenes of the Mst77F gene, which is a testis-expressed autosomal gene that is essential for male fertility. We did a thorough search using experimental and computational methods and found 18 Y-linked copies of this gene (named Mst77Y-1-Mst77Y-18). Ten Mst77Y genes encode defective proteins and the other eight are potentially functional. These eight genes produce approximately 20% of the functional Mst77F-like mRNA, and molecular evolutionary analysis shows that they evolved under purifying selection. Hence several Mst77Y genes have all the features of functional genes. Mst77Y genes are present only in D. melanogaster, and phylogenetic analysis confirmed that the duplication is a recent event. The identification of functional Mst77Y genes reinforces the previous finding that gene gains play a prominent role in the evolution of the Drosophila Y chromosome.

A fish eye out of water: ten visual opsins in the four-eyed fish, Anableps anableps

The “four-eyed” fish Anableps anableps has numerous morphological adaptations that enable above and below-water vision. Here, as the first step in our efforts to identify molecular adaptations for aerial and aquatic vision in this species, we describe the A. anableps visual opsin repertoire. We used PCR, cloning, and sequencing to survey cDNA using unique primers designed to amplify eight sequences from five visual opsin gene subfamilies, SWS1, SWS2, RH1, RH2, and LWS. We also used Southern blotting to count opsin loci in genomic DNA digested with EcoR1 and BamH1. Phylogenetic analyses confirmed the identity of all opsin sequences and allowed us to map gene duplication and divergence events onto a tree of teleost fish. Each of the gene-specific primer sets produced an amplicon from cDNA, indicating that A. anableps possessed and expressed at least eight opsin genes. A second PCR-based survey of genomic and cDNA uncovered two additional LWS genes. Thus, A. anableps has at least ten visual opsins and all but one were expressed in the eyes of the single adult surveyed. Among these ten visual opsins, two have key site haplotypes not found in other fish. Of particular interest is the A. anableps-specific opsin in the LWS subfamily, S180γ, with a SHYAA five key site haplotype. Although A. anableps has a visual opsin gene repertoire similar to that found in other fishes in the suborder Cyprinodontoidei, the LWS opsin subfamily has two loci not found in close relatives, including one with a key site haplotype not found in any other fish species. A. anableps opsin sequence data will be used to design in situ probes allowing us to test the hypothesis that opsin gene expression differs in the distinct ventral and dorsal retinas found in this species.

Fungi from the roots of the common terrestrial orchid Gymnadenia conopsea

The fungal community associated with the terrestrial photosynthetic orchid Gymnadenia conopsea was characterized through PCR-amplification directly from root extracted DNA and cloning of the PCR products. Six populations in two geographically distinct regions in Germany were investigated. New ITS-primers amplifying a wide taxonomic range including Basidiomycetes and Ascomycetes revealed a high taxonomic and ecological diversity of fungal associates, including typical orchid mycorrhizas of the Tulasnellaceae and Ceratobasidiaceae as well as several ectomycorrhizal taxa of the Pezizales. The wide spectrum of potential mycorrhizal partners may contribute to this orchid's ability to colonize different habitat types with their characteristic microbial communities. The fungal community of G. conopsea showed a clear spatial structure. With 43 % shared taxa the species composition of the two regions showed only little overlap. Regardless of regions, populations were highly variable concerning taxon richness, varying between 5 and 14 taxa per population. The spatial structure and the continuous presence of mycorrhizal taxa on the one hand and the low specificity towards certain fungal taxa on the other hand suggest that the fungal community associated with G. conopsea is determined by multiple factors. In this context, germination as well as pronounced morphological and genetic differentiation within G. conopsea deserve attention as potential factors affecting the composition of the fungal community.

Molecular diversity of visual pigments in Stomatopoda (Crustacea)

Stomatopod crustaceans possess apposition compound eyes that contain more photoreceptor types than any other animal described. While the anatomy and physiology of this complexity have been studied for more than two decades, few studies have investigated the molecular aspects underlying the stomatopod visual complexity. Based on previous studies of the structure and function of the different types of photoreceptors, stomatopod retinas are hypothesized to contain up to 16 different visual pigments, with 6 of these having sensitivity to middle or long wavelengths of light. We investigated stomatopod middle- and long-wavelength-sensitive opsin genes from five species with the hypothesis that each species investigated would express up to six different opsin genes. In order to understand the evolution of this class of stomatopod opsins, we examined the complement of expressed transcripts in the retinas of species representing a broad taxonomic range (four families and three superfamilies). A total of 54 unique retinal opsins were isolated, resulting in 6–15 different expressed transcripts in each species. Phylogenetically, these transcripts form six distinct clades, grouping with other crustacean opsins and sister to insect long-wavelength visual pigments. Within these stomatopod opsin groups, intra- and interspecific clusters of highly similar transcripts suggest that there has been rampant recent gene duplication. Some of the observed molecular diversity is also due to ancient gene duplication events within the stem crustacean lineage. Using evolutionary trace analysis, 10 amino acid sites were identified as functionally divergent among the six stomatopod opsin clades. These sites form tight clusters in two regions of the opsin protein known to be functionally important: six in the chromophore-binding pocket and four at the cytoplasmic surface in loops II and III. These two clusters of sites indicate that stomatopod opsins have diverged with respect to both spectral tuning and signal transduction.

Major histocompatibility complex (MHC) heterozygote superiority to natural multi-parasite infections in the water vole (Arvicola terrestris)

The fundamental role of the major histocompatibility complex (MHC) in immune recognition has led to a general consensus that the characteristically high levels of functional polymorphism at MHC genes is maintained by balancing selection operating through host-parasite coevolution. However, the actual mechanism by which selection operates is unclear. Two hypotheses have been proposed: overdominance (or heterozygote superiority) and negative frequency-dependent selection. Evidence for these hypotheses was evaluated by examining MHC-parasite relationships in an island population of water voles (Arvicola terrestris). Generalized linear mixed models were used to examine whether individual variation at an MHC class II DRB locus explained variation in the individual burdens of five different parasites. MHC genotype explained a significant amount of variation in the burden of gamasid mites, fleas (Megabothris walkeri) and nymphs of sheep ticks (Ixodes ricinus). Additionally, MHC heterozygotes were simultaneously co-infected by fewer parasite types than homozygotes. In each case where an MHC-dependent effect on parasite burden was resolved, the heterozygote genotype was associated with fewer parasites, and the heterozygote outperformed each homozygote in two of three cases, suggesting an overall superiority against parasitism for MHC heterozygote genotypes. This is the first demonstration of MHC heterozygote superiority against multiple parasites in a natural population, a mechanism that could help maintain high levels of functional MHC genetic diversity in natural populations.

Spatio-temporal variation in the strength and mode of selection acting on major histocompatibility complex diversity in water vole (Arvicola terrestris) metapopulations

Patterns of spatio-temporal genetic variation at a class II major histocompatibility complex (MHC) locus and multiple microsatellite loci were analysed within and between three water vole metapopulations in Scotland, UK. Comparisons of MHC and microsatellite spatial genetic differentiation, based on standardised tests between two demographically asynchronous zones within a metapopulation, suggested that spatial MHC variation was affected by balancing selection, directional selection and random genetic drift, but that the relative effects of these microevolutionary forces vary temporally. At the metapopulation level, between-year differentiation for MHC loci was significantly correlated with that of microsatellites, signifying that neutral factors such as migration and drift were primarily responsible for overall temporal genetic change at the metapopulation scale. Between metapopulations, patterns of genetic differentiation implied that, at large spatial scales, MHC variation was primarily affected by directional selection and drift. Levels of MHC heterozygosity in excess of Hardy-Weinberg expectations were consistent with overdominant balancing selection operating on MHC variation within metapopulations. However, this effect was not constant among all samples, indicating temporal variation in the strength of selection relative to other factors. The results highlight the benefit of contrasting variation at MHC with neutral markers to separate the effects of stochastic and deterministic microevolutionary forces, and add to a growing body of evidence showing that the mode and relative strength of selection acting on MHC diversity varies both spatially and temporally.

Mechanism of somatic hypermutation: critical analysis of strand biased mutation signatures at A:T and G:C base pairs

The DNA sequence data of the somatic hypermutation (SHM) field published since 1984 has been critically reviewed. The analysis has revealed three strand biased mutation signatures. The first concerns the mutations generated at G:C base pairs in mice genetically deficient in uracil-DNA glycosylase and MSH2-MSH6-mediated mismatch repair. Such mice display the AID deaminase footprint and here C mutations exceed G mutations at least 1.5-fold. This supports earlier and more recent studies claiming that dC-to-dU deaminations occur preferentially in the single stranded DNA regions of the displaced nontranscribed strand (NTS) during transcription. The second concerns the signature generated in immunised mice where G mutations exceed C mutations by at least 1.7-fold. This is a newly identified strand bias which has previously gone undetected. It is consistent with the polynucleotide polymerisation signature of RNA polymerase II copying the template DNA strand carrying AID-mediated lesions generated at C bases, viz. uracils and abasic sites. A reverse transcription step would then need to intervene to fix the mutation pattern in DNA. The third concerns the long recognised strand biased signature generated in normal aged or actively immunised mice whereby A mutations exceed T mutations by two- to three-fold. It is argued that this pattern is best understood as a combination of adenosine-to-inosine (A-to-I) RNA editing followed by a reverse transcription step fixing the A-to-G, as well as A-to-T and A-to-C, as strand biased mutation signatures in DNA. The reasons why the AID-linked RNA polymerase II mutation signature had previously gone undetected are discussed with regard to limitations of standard PCR-based SHM assay techniques. It is concluded that the most economical SHM mechanism involves both DNA and RNA deaminations coupled to a reverse transcription process, most likely involving DNA polymerase eta acting in its reverse transcriptase mode. Experimental approaches to differentiate this RNA-based model from the standard DNA deamination model are discussed.

The molecular basis of color vision in colorful fish: four long wave-sensitive (LWS) opsins in guppies (Poecilia reticulata) are defined by amino acid substitutions at key functional sites

BACKGROUND

Comparisons of functionally important changes at the molecular level in model systems have identified key adaptations driving isolation and speciation. In cichlids, for example, long wavelength-sensitive (LWS) opsins appear to play a role in mate choice and male color variation within and among species. To test the hypothesis that the evolution of elaborate coloration in male guppies (Poecilia reticulata) is also associated with opsin gene diversity, we sequenced long wavelength-sensitive (LWS) opsin genes in six species of the family Poeciliidae.

RESULTS

Sequences of four LWS opsin genes were amplified from the guppy genome and from mRNA isolated from adult guppy eyes. Variation in expression was quantified using qPCR. Three of the four genes encode opsins predicted to be most sensitive to different wavelengths of light because they vary at key amino acid positions. This family of LWS opsin genes was produced by a diversity of duplication events. One, an intronless gene, was produced prior to the divergence of families Fundulidae and Poeciliidae. Between-gene PCR and DNA sequencing show that two of the guppy LWS opsins are linked in an inverted orientation. This inverted tandem duplication event occurred near the base of the poeciliid tree in the common ancestor of Poecilia and Xiphophorus. The fourth sequence has been uncovered only in the genus Poecilia. In the guppies surveyed here, this sequence is a hybrid, with the 5' end most similar to one of the tandem duplicates and the 3' end identical to the other.

CONCLUSION

Enhanced wavelength discrimination, a possible consequence of opsin gene duplication and divergence, might have been an evolutionary prerequisite for color-based sexual selection and have led to the extraordinary coloration now observed in male guppies and in many other poeciliids.

Wolbachia transmission dynamics in Formica wood ants

BACKGROUND

The role of Wolbachia endosymbionts in shaping the mitochondrial diversity of their arthropod host depends on the effects they have on host reproduction and on the mode of transmission of the bacteria. We have compared the sequence diversity of wsp (Wolbachia surface protein gene) and the host mtDNA in a group of Formica ant species that have diverged approximately 0.5 million years ago (MYA). The aim was to study the relationship of Wolbachia and its ant hosts in terms of vertical and horizontal transmission of the bacteria.

RESULTS

All studied ant species were doubly infected with two Wolbachia strains (wFex1 and wFex4) all over their geographical distribution area in Eurasia. The most common haplotypes of these strains were identical with strains previously described from a more distantly related Formica ant, with an estimated divergence time of 3.5 - 4 MYA. Some strain haplotypes were associated to the same or closely related mtDNA haplotypes as expected under vertical transmission. However, in several cases the wsp haplotypes coexisted with distant mtDNA haplotypes, a pattern which is more compatible with horizontal transmission of the bacteria.

CONCLUSION

Two lines of evidence suggest that the sharing of Wolbachia strains by all F. rufa species is rather due to horizontal than vertical transmission. First, the fact that endosymbiont strains identical to those of F. rufa ants have been found in another species that diverged 3.5-4 MYA strongly suggests that horizontal transfer can and does occur between Formica ants. Second, the frequent sharing of identical Wolbachia strains by distant mitochondrial lineages within the F. rufa group further shows that horizontal transmission has occurred repeatedly. Nevertheless, our dataset also provides some evidence for longer-term persistence of infection, indicating that Wolbachia infection within this host clade has been shaped by both horizontal and vertical transmission of symbionts. The fact that all the ants were infected irrespective of the family structure of their societies gives no support to the proposed hypotheses that the spreading of Wolbachia in ants might be associated to the types of their societies.

Potential inhibition of somatic hypermutation by nucleoside analogues

Somatic hypermutation, which occurs in antigen-activated germinal centre B lymphocytes, diversifies the genes that encode immunoglobulin variable regions and leads to the 'affinity maturation' of the humoral immune response. Hypermutation affects dC/dG and dA/dT pairs with approximately equal frequency in vivo. DNA polymerase-theta contributes to hypermutagenesis at dC/dG pairs and DNA polymerase-eta is substantially involved in the generation of hypermutations at dA/dT pairs. The biochemical properties of polymerases-theta and -eta indicate that their DNA synthetic activities are potentially susceptible to inhibition by nucleoside analogues, so it is feasible that nucleoside analogues reduce the accumulation of dC/dG- and dA/dT-targeted hypermutations in vivo. Nucleoside analogues could hence impair the humoral adaptive immune response of HIV-infected patients who are prescribed these chemotherapeutic agents.

No evidence for the use of DIR, D-D fusions, chromosome 15 open reading frames or VH replacement in the peripheral repertoire was found on application of an improved algorithm, JointML, to 6329 human immunoglobulin H rearrangements

Antibody diversity is created by imprecise joining of the variability (V), diversity (D) and joining (J) gene segments of the heavy and light chain loci. Analysis of rearrangements is complicated by somatic hypermutations and uncertainty concerning the sources of gene segments and the precise way in which they recombine. It has been suggested that D genes with irregular recombination signal sequences (DIR) and chromosome 15 open reading frames (OR15) can replace conventional D genes, that two D genes or inverted D genes may be used and that the repertoire can be further diversified by heavy chain V gene (VH) replacement. Safe conclusions require large, well-defined sequence samples and algorithms minimizing stochastic assignment of segments. Two computer programs were developed for analysis of heavy chain joints. JointHMM is a profile hidden Markow model, while JointML is a maximum-likelihood-based method taking the lengths of the joint and the mutational status of the VH gene into account. The programs were applied to a set of 6329 clonally unrelated rearrangements. A conventional D gene was found in 80% of unmutated sequences and 64% of mutated sequences, while D-gene assignment was kept below 5% in artificial (randomly permutated) rearrangements. No evidence for the use of DIR, OR15, multiple D genes or VH replacements was found, while inverted D genes were used in less than 1 per thousand of the sequences. JointML was shown to have a higher predictive performance for D-gene assignment in mutated and unmutated sequences than four other publicly available programs. An online version 1.0 of JointML is available at http://www.cbs.dtu.dk/services/VDJsolver.

Seven different genes encode a diverse mixture of isoforms of Bet v 1, the major birch pollen allergen

Background

Pollen of the European white birch (Betula pendula, syn. B. verrucosa) is an important cause of hay fever. The main allergen is Bet v 1, member of the pathogenesis-related class 10 (PR-10) multigene family. To establish the number of PR-10/Bet v 1 genes and the isoform diversity within a single tree, PCR amplification, cloning and sequencing of PR-10 genes was performed on two diploid B. pendula cultivars and one interspecific tetraploid Betula hybrid. Sequences were attributed to putative genes based on sequence identity and intron length. Information on transcription was derived by comparison with homologous cDNA sequences available in GenBank/EMBL/DDJB. PCR-cloning of multigene families is accompanied by a high risk for the occurrence of PCR recombination artifacts. We screened for and excluded these artifacts, and also detected putative artifact sequences among database sequences.

Results

Forty-four different PR-10 sequences were recovered from B. pendula and assigned to thirteen putative genes. Sequence homology suggests that three genes were transcribed in somatic tissue and seven genes in pollen. The transcription of three other genes remains unknown. In total, fourteen different Bet v 1-type isoforms were identified in the three cultivars, of which nine isoforms were entirely new. Isoforms with high and low IgE-reactivity are encoded by different genes and one birch pollen grain has the genetic background to produce a mixture of isoforms with varying IgE-reactivity. Allergen diversity is even higher in the interspecific tetraploid hybrid, consistent with the presence of two genomes.

Conclusion

Isoforms of the major birch allergen Bet v 1 are encoded by multiple genes, and we propose to name them accordingly. The present characterization of the Bet v 1 genes provides a framework for the screening of specific Bet v 1 genes among other B. pendula cultivars or Betula species, and for future breeding for trees with a reduced allergenicity. Investigations towards sensitization and immunotherapy should anticipate that patients are exposed to a mixture of Bet v 1 isoforms of different IgE-reactivity, even if pollen originates from a single birch tree.

Isolation and characterization of a MHC class II DRB locus in the European water vole (Arvicola terrestris)

In so-called model species, such as human and mouse, genes of the major histocompatibility complex (MHC) are characterized by extremely high levels of polymorphism, and it is considered that such diversity is maintained by balancing selection. ;There is now a recognized need to expand studies into nonmodel species to examine whether high MHC diversity is mirrored in natural populations, and to determine the ecological, ethological, and evolutionary processes that underpin balancing selection. To address such issues, a necessary prerequisite is the ability to characterize diversity at a single, expressed, polymorphic MHC locus on which selection may be acting. Here, we provide the first description of allelic diversity at exon 2 of an MHC class II DRB locus in the European water vole (Arvicola terrestris), characterize variation across four natural populations, and test whether the patterns of variation are consistent with the effects of balancing selection. Using single-strand conformation polymorphism analysis and subsequent DNA sequencing of gel excisions, five DRB alleles were resolved, each with a unique amino acid sequence, among 100 individuals from four geographically distinct populations. Reverse transcription polymerase chain reaction confirmed that the alleles were products from an expressed locus. Intra-allelic amino acid differences were high (10.5-33.3%), and the nonsynonymous substitution rate exceeded the synonymous substitution rate for the functional peptide-binding region (d (N):d (S)=3.91 and P<0.005). Phylogenetic comparison of resolved alleles with closely related homologues indicated that each allele represented a unique lineage preserved across speciation events. These results indicate that balancing selection has maintained diversity of DRB allelic lineages and amino acid function over evolutionary time scales, but may be less effective at preserving alleles in contemporary populations where stochastic microevolutionary processes may dominate.

Reconsidering the human immunoglobulin heavy-chain locus:

We have used a bioinformatics approach to evaluate the completeness and functionality of the reported human immunoglobulin heavy-chain IGHD gene repertoire. Using the hidden Markov-model-based iHMMune-align program, 1,080 relatively unmutated heavy-chain sequences were aligned against the reported repertoire. These alignments were compared with alignments to 1,639 more highly mutated sequences. Comparisons of the frequencies of gene utilization in the two databases, and analysis of features of aligned IGHD gene segments, including their length, the frequency with which they appear to mutate, and the frequency with which specific mutations were seen, were used to determine the reliability of alignments to the less commonly seen IGHD genes. Analysis demonstrates that IGHD4-23 and IGHD5-24, which have been reported to be open reading frames of uncertain functionality, are represented in the expressed gene repertoire; however, the functionality of IGHD6-25 must be questioned. Sequence similarities make the unequivocal identification of members of the IGHD1 gene family problematic, although all genes except IGHD1-14*01 appear to be functional. On the other hand, reported allelic variants of IGHD2-2 and of the IGHD3 gene family appear to be nonfunctional, very rare, or nonexistent. Analysis also suggests that the reported repertoire is relatively complete, although one new putative polymorphism (IGHD3-10*p03) was identified. This study therefore confirms a surprising lack of diversity in the available IGHD gene repertoire, and restriction of the germline sequence databases to the functional set described here will substantially improve the accuracy of IGHD gene alignments and therefore the accuracy of analysis of the V-D-J junction.