Multiplex genotype determination at a DNA sequence polymorphism cluster in the human immunoglobulin heavy-chain region

NKX2-3 variant rs11190140 is associated with IBD and alters binding of NFAT

NKX2-3 SNP rs11190140 is associated with inflammatory bowel disease (IBD). The T allele is over-transmitted in IBD and the C allele represents a potential CpG methylation site. We hypothesize that genetic variation and/or methylation of SNP rs11190140 may play a role in NKX2-3 gene expression by affecting transcription factor binding. We studied 233 IBD cases and 250 unrelated healthy individuals from an IBD population from central Pennsylvania and performed genotype analyses of the genetic variation and methylation status analysis using PCR-based RFLP. For transcription factor binding, nuclear extracts from human B cells were incubated with biotin-labeled oligonucleotide sequences of the NKX2-3 promoter region containing the genetic variation of T, non-methylated C or methylated C at rs11190140, followed by biotin pull-down and Western blot analysis for transcription factors SP1, NFAT1, NF-κB, and ETS-1. In case-control analysis, the genetic variation was significantly associated with IBD (OR=0.503, 95% CI=0.330-0.764, p<0.001). Methylation status analyses revealed that the C allele is subject to modification by DNA methylation. transcription factor binding assay indicated distinct differential binding of NFAT1 to the NKX2-3 promoter sequence, with higher binding to those with non-methylated and methylated C than to T. The binding of NFAT1 to the NKX2-3 promoter region with rs1190140 was confirmed by ChIP assay. We speculate that the rs11190140 may regulate NKX2-3 expression and have a role in IBD pathogenesis.

A simple PCR-based genotyping method for M105I mutation of alpha-SNAP enhances the study of early pathological changes in hyh phenotype

alpha-SNAP is an essential component of the protein machinery responsible for membrane fusion events in different cell types. The hyh (hydrocephalus with hop gait) mouse carries a missense mutation in Napa gene that results in a point mutation (M105I) in alpha-SNAP protein. Homozygous animals for the mutant allele have been identified by the clinical and/or neuropathological phenotype, or by direct sequencing of PCR products. The aims of the present study were (i) to develop a high-throughput technique to genotype hyh mice, (ii) to correlate genotype-phenotype, and (iii) to analyze the earliest pathological changes of hyh mutant mice. As no restriction sites are affected by the hyh mutation, we resolved this problem by creating a BspHI restriction site with a modified (mismatch) polymerase chain reaction (PCR) primer in wild-type allele. This artificially created restriction site (ACRS)-PCR technique is a simple, rapid and reliable method to genotype hyh mice in a day-work procedure. Biochemical and histological analysis of genotyped hyh embryos at different developmental stages allowed us to identify and characterize the earliest brain pathological changes of the hyh phenotype, including the first signs of neuroepithelial disruption and neuronal ectopia. In addition, genotype-phenotype analysis of 327 animals confirmed that (i) hyh is a single-gene autosomal recessive disorder, and (ii) the disorder has 100% penetrance (i.e., the mutation was only present in affected mice). The genotyping method described here enhances the potentiality of hyh mouse as a unique in vivo model to study the role of membrane trafficking in different developmental and physiological processes.

Determination of gene organization in the human IGHV region on single chromosomes

Organization of the IGHV genes (n=108) on single human chromosomes has been determined by detecting these sequences in single sperm using multiplex PCR amplification followed by microarray detection. A total of 374 single sperm samples from five Caucasian males were studied. Three deletion/insertion polymorphisms (Del I-Del III) with deletion allele frequencies ranging from 0.1 to 0.3 were identified. Del I is a previously reported polymorphism affecting three IGHV genes (IGHV1-8, IGHV3-9, and IGHV2-10). Del II affects a region 2-18 kb containing two pseudogenes IGHV(II)-28.1 and IGHV3-29, and Del III spans approximately 21-53 kb involving genes IGHV4-39, IGHV7-40, IGHV(II)-40-1, and IGHV3-41. Deletion alleles of both Dels II and III were found in a heterozygous state, and therefore, could not be easily detected if haploid samples were not used in the study. Results of the present study indicate that deletions/insertions together with other possible chromosomal rearrangements may play an important role in forming the genetic structure of the IGHV region, and may significantly contribute to antibody diversity. Since these three polymorphisms are located within or next to the 3' half of the IGHV region, they may have an important role in the expressed IGHV gene repertoire during immune response.

The topoisomerase I- and p53-binding protein topors is differentially expressed in normal and malignant human tissues and may function as a tumor suppressor

Topors was identified recently as a human protein that binds to topoisomerase I and p53. Topors contains a highly conserved RING domain and localizes in promyelocytic leukemia nuclear bodies. Relatively little is known regarding topors expression patterns or function. We now demonstrate that topors mRNA and protein are widely expressed in normal human tissues. By contrast, topors mRNA and protein levels are decreased or undetectable in colon adenocarcinomas relative to normal colon tissue, and expression of the topors protein is not detectable in several colon cancer cell lines. The human TOPORS gene is located on chromosome 9p21, with loss of heterozygosity in this region frequently observed in several different malignancies. While we were unable to detect loss of heterozygosity of the TOPORS gene in 16 sporadic colon cancer cases, increased methylation of a CpG island in the TOPORS promoter was evident in colon adenocarcinoma specimens relative to matched normal tissues. Additional studies indicate that forced expression of topors inhibits cellular proliferation and is associated with an accumulation of cells in the G(0)/G(1) phase of the cell cycle. This effect is independent of the topors RING domain and maps to a C-terminal region of the protein. These results suggest that topors functions as a negative regulator of cell growth, and possibly as a tumor suppressor.

Multiplex loss of heterozygosity analysis by using single or very few cells

Loss of heterozygosity (LOH) is an indication of tumor suppressor gene inactivation. However, loss of heterozygosity analysis has been limited to either a small scale or to very few genetic markers. To significantly increase the scale of study and to include a large number of markers in the analysis, experimental conditions were established for using single cells or single cell equivalent with 10 markers typed simultaneously. Under these conditions, the allele amplification failure rate was 3.7% when single tissue cultured human cells were used. When 30 cells from a 5- micro m paraffin-archived breast tumor tissue section were used, the failure rates were 0% for four of the five heterozygous loci and 10% for the fifth. Small amplification failure rates (6.1% and 6.7% on average) were observed when 5 or 10 cells from paraffin-archived breast tissue were used. These results indicate that with polymerase chain reaction (PCR) primers of high quality, it is possible to obtain reliable results by using single cells from fresh tissue or very few cells from paraffin-archived specimens. The results also show the importance of including replicates, using primers of high quality, and optimizing PCR conditions when a limited amount of material is used for the assay. The feasibility of LOH analysis with very few paraffin-embedded breast cancer tissues was demonstrated.

Diversifying selection of the tumor-growth promoter angiogenin in primate evolution

Diversifying selection drives the rapid differentiation of gene sequences and is one of the main forces behind adaptive evolution. Most genes known to be shaped by diversifying selection are those involved in host-pathogen or male-female interactions characterized as molecular "arms races." Here we report the unexpected detection of diversifying selection in the evolution of a tumor-growth promoter, angiogenin (ANG). A comparison among 11 primate species demonstrates that ANG has a significantly higher rate of nucleotide substitution at nonsynonymous sites than at synonymous sites, a hallmark of positive selection acting at the molecular level. Furthermore, we observed significant charge diversity at the molecular surface, suggesting the presence of selective pressures in the microenvironment of ANG, including its binding molecules. A population survey of ANG in chimpanzees, however, reveals no polymorphism, which may have resulted from a recent selective sweep of a charge-altering substitution in chimpanzee evolution. Functional assays of recombinant ANGs from the human and owl monkey indicate that primate ANGs retain angiogenic activity despite rapid evolution. Our study, together with findings of similar selection in the primate breast cancer suppressor gene, BRCA1, reveals an intriguing phenomenon of unusual selective pressures on, and adaptive evolution of, cancer-related genes in primate evolution.

Robust and accurate single nucleotide polymorphism genotyping by dynamic allele-specific hybridization (DASH): design criteria and assay validation

We recently introduced a generic single nucleotide polymorphism (SNP) genotyping method, termed DASH (dynamic allele-specific hybridization), which entails dynamic tracking of probe (oligonucleotide) to target (PCR product) hybridization as reaction temperature is steadily increased. The reliability of DASH and optimal design rules have not been previously reported. We have now evaluated crudely designed DASH assays (sequences unmodified from genomic DNA) for 89 randomly selected and confirmed SNPs. Accurate genotype assignment was achieved for 89% of these worst-case-scenario assays. Failures were determined to be caused by secondary structures in the target molecule, which could be reliably predicted from thermodynamic theory. Improved design rules were thereby established, and these were tested by redesigning six of the failed DASH assays. This involved reengineering PCR primers to eliminate amplified target sequence secondary structures. This sophisticated design strategy led to complete functional recovery of all six assays, implying that SNPs in most if not all sequence contexts can be effectively scored by DASH. Subsequent empirical support for this inference has been evidenced by approximately 30 failure-free DASH assay designs implemented across a range of ongoing genotyping programs. Structured follow-on studies employed standardized assay conditions, and revealed that assay reproducibility (733 duplicated genotypes, six different assays) was as high as 100%, with an assay accuracy (1200 genotypes, three different assays) that exceeded 99.9%. No post-PCR assay failures were encountered. These findings, along with intrinsic low cost and high flexibility, validate DASH as an effective procedure for SNP genotyping.

Sequence variation at two eosinophil-associated ribonuclease loci in humans

Host defense against invading pathogens is of great importance to the survival of higher organisms. We have been studying the evolution of mammalian eosinophil-associated ribonucleases (EARs), which are members of the ribonuclease A superfamily with known antipathogen activities. Earlier studies showed that positive selection promoted rapid diversification of paralogous EAR genes in both primates and rodents. Intraspecifically, however, it is unknown whether these genes also have divergent alleles. The recent discovery that the gene repertoire of the EAR family is much larger in rodents than in primates has led us to consider the possibility that primates maintain a large number of polymorphic alleles to compensate for a smaller gene repertoire. Here we present sequences of 2417 nucleotides at the two EAR loci, the eosinophil-derived neurotoxin (EDN, RNase 2) and eosinophil cationic protein (ECP, RNase 3), from >50 human individuals. Our data demonstrate that the nucleotide diversities (0.06-0.11%) at these loci are typical for human nuclear genes, thus permitting us to reject this polymorphism hypothesis. No significant departure from neutrality is noted and no signs of overdominant selection are observed. Similar patterns were observed in a preliminary study of chimpanzees. In summary, our results suggest that the antipathogen functions of the primate EARs are conserved after they are established and that these proteins are not currently undergoing rapid diversification in response to challenge from invading microorganisms.

Evaluation of single nucleotide polymorphism typing with invader on PCR amplicons and its automation

Large-scale pharmacogenetics and complex disease association studies will require typing of thousands of single-nucleotide polymorphisms (SNPs) in thousands of individuals. Such projects would benefit from a genotyping system with accuracy >99% and a failure rate <5% on a simple, reliable, and flexible platform. However, such a system is not yet available for routine laboratory use. We have evaluated a modification of the previously reported Invader SNP-typing chemistry for use in a genotyping laboratory and tested its automation. The Invader technology uses a Flap Endonuclease for allele discrimination and a universal fluorescence resonance energy transfer (FRET) reporter system. Three hundred and eighty-four individuals were genotyped across a panel of 36 SNPs and one insertion/deletion polymorphism with Invader assays using PCR product as template, a total of 14,208 genotypes. An average failure rate of 2.3% was recorded, mostly associated with PCR failure, and the typing was 99.2% accurate when compared with genotypes generated with established techniques. An average signal-to-noise ratio (9:1) was obtained. The high degree of discrimination for single base changes, coupled with homogeneous format, has allowed us to deploy liquid handling robots in a 384-well microtitre plate format and an automated end-point capture of fluorescent signal. Simple semiautomated data interpretation allows the generation of approximately 25,000 genotypes per person per week, which is 10-fold greater than gel-based SNP typing and microsatellite typing in our laboratory. Savings on labor costs are considerable. We conclude that Invader chemistry using PCR products as template represents a useful technology for typing large numbers of SNPs rapidly and efficiently.

Fifty-million-year-old polymorphism at an immunoglobulin variable region gene locus in the rabbit evolutionary lineage

The Ig heavy chain variable region (V(H)) genes encode the antigen-binding regions of antibodies. The rabbit genome contains more than 100 V(H) genes, but only one (V(H)1) is preferentially used in the VDJ gene rearrangement. Three highly divergent alleles occur at this V(H)1 locus in most rabbit populations. These three V(H) alleles are also present in snowshoe hare populations, indicating that the polymorphism of the V(H)1 alleles is trans-specific. Here we report the results of a phylogenetic analysis of rabbit Ig germ-line V(H) genes (alleles) together with V(H) genes from humans and mice. We have found that all rabbit V(H) genes belong to one mammalian V(H) group (group C), which also includes various human and mouse V(H) genes. Using the rate of nucleotide substitution obtained from human and mouse V(H) sequences, we have estimated that the V(H)1 polymorphism in the rabbit lineage has been maintained for about 50 million years. This extremely long persistence of V(H)1 polymorphism is apparently caused by overdominant selection, though the real mechanism is unclear.

An SP-B gene mutation responsible for SP-B deficiency in fatal congenital alveolar proteinosis: evidence for a mutation hotspot in exon 4

Mutations and polymorphisms within the human SP-B locus have been linked to fatal congenital alveolar proteinosis (CAP) and associated with respiratory distress syndrome (RDS), respectively. In the present study we used PCR and direct sequence analysis of the SP-B gene of three individuals from a family with CAP to search for additional SP-B mutations resulting in CAP and/or polymorphisms that could be used as markers in association studies of RDS and/or CAP. We found three novel mutations/polymorphisms in this family. One is a C/A substitution at nt 1013 at the splice junction of intron 2-exon 3. A second one is a single base T deletion at nt 1553 in exon 4. The single base (T) deletion at nucleotide 1553 (1553delT) shifts the reading frame at amino acid 122(122delT) and creates a premature termination codon at amino acid 214 in exon 6. The mutated gene produces no mature SP-B protein. Genotype analysis from the nuclear family carrying this mutation showed that both parents and three of the four living children are heterozygous for the mutation. One of the four living children is homozygous for the normal allele and a child that died in the perinatal period from CAP is homozygous for the mutation. A third change is a C/T substitution at nt 1580 in exon 4 that changes amino acid 131 from threonine to isoleucine (Thr131Ile). The location of a previously reported mutation, 121ins2 (1), is only 4 nt upstream of 122delT, and the missense mutation Thr131Ile (exon 4) is only 27 nt downstream of 122delT. These changes are within or in close proximity to a CCTG sequence and a poly(C) tract, both of which are shown in other systems to be mutation hotspots. The 122delT occurs within the CCTG and the poly(C) tract sequences, the Thr131Ile occurs 26 nt downstream from the CCTG sequence, and the 121ins2 occurs 2 nt upstream from CCTG sequence and within the poly(C) tract. The present observations suggest that the short SP-B sequence containing the CCTG motif and the poly(C) tract is a mutation hotspot.

Evolution by the birth-and-death process in multigene families of the vertebrate immune system

Concerted evolution is often invoked to explain the diversity and evolution of the multigene families of major histocompatibility complex (MHC) genes and immunoglobulin (Ig) genes. However, this hypothesis has been controversial because the member genes of these families from the same species are not necessarily more closely related to one another than to the genes from different species. To resolve this controversy, we conducted phylogenetic analyses of several multigene families of the MHC and Ig systems. The results show that the evolutionary pattern of these families is quite different from that of concerted evolution but is in agreement with the birth-and-death model of evolution in which new genes are created by repeated gene duplication and some duplicate genes are maintained in the genome for a long time but others are deleted or become nonfunctional by deleterious mutations. We found little evidence that interlocus gene conversion plays an important role in the evolution of MHC and Ig multigene families.

Multiplex genotype determination at a large number of gene loci

To facilitate large-scale genotype analysis, an efficient PCR-based multiplex approach has been developed. For simultaneously amplifying the target sequences at a large number of genetic loci, locus-specific primers containing 5' universal tails are used. Attaching the universal tails to the target sequences in the initial PCR steps allows replacement of all specific primers with a pair of primers identical to the universal tails and converts the multiplex amplification into "uniplex." Simultaneous amplification of 26 genetic loci with this approach is described. The multiplex amplification can be coupled with genotype determination. By incorporating a single-base mismatch between a primer and the template into the target sequences, a polymorphic site can be converted into a desirable restriction fragment length polymorphism when it is necessary. In this way, the allelic PCR products for the polymorphic loci can be discriminated by gel electrophoresis after restriction enzyme digestion. In this study, 32 loci were typed in such a multiplex way.