Analysis of 454 sequencing error rate, error sources, and artifact recombination for detection of Low-frequency drug resistance mutations in HIV-1 DNA.

BACKGROUND: 454 sequencing technology is a promising approach for characterizing HIV-1 populations and for identifying low frequency mutations. The utility of 454 technology for determining allele frequencies and linkage associations in HIV infected individuals has not been extensively investigated. We evaluated the performance of 454 sequencing for characterizing HIV populations with defined allele frequencies. RESULTS: We constructed two HIV-1 RT clones. Clone A was a wild type sequence. Clone B was identical to clone A except it contained 13 introduced drug resistant mutations. The clones were mixed at ratios ranging from 1% to 50% and were amplified by standard PCR conditions and by PCR conditions aimed at reducing PCR-based recombination. The products were sequenced using 454 pyrosequencing. Sequence analysis from standard PCR amplification revealed that 14% of all sequencing reads from a sample with a 50:50 mixture of wild type and mutant DNA were recombinants. The majority of the recombinants were the result of a single crossover event which can happen during PCR when the DNA polymerase terminates synthesis prematurely. The incompletely extended template then competes for primer sites in subsequent rounds of PCR. Although less often, a spectrum of other distinct crossover patterns was also detected. In addition, we observed point mutation errors ranging from 0.01% to 1.0% per base as well as indel (insertion and deletion) errors ranging from 0.02% to nearly 50%. The point errors (single nucleotide substitution errors) were mainly introduced during PCR while indels were the result of pyrosequencing. We then used new PCR conditions designed to reduce PCR-based recombination. Using these new conditions, the frequency of recombination was reduced 27-fold. The new conditions had no effect on point mutation errors. We found that 454 pyrosequencing was capable of identifying minority HIV-1 mutations at frequencies down to 0.1% at some nucleotide positions. CONCLUSION: Standard PCR amplification results in a high frequency of PCR-introduced recombination precluding its use for linkage analysis of HIV populations using 454 pyrosequencing. We designed a new PCR protocol that resulted in a much lower recombination frequency and provided a powerful technique for linkage analysis and haplotype determination in HIV-1 populations. Our analyses of 454 sequencing results also demonstrated that at some specific HIV-1 drug resistant sites, mutations can reliably be detected at frequencies down to 0.1%.

Lens epithelium-derived growth factor fusion proteins redirect HIV-1 DNA integration.

Lens epithelium-derived growth factor (LEDGF) fusion proteins can direct HIV-1 DNA integration to novel sites in the host genome. The C terminus of LEDGF contains an integrase binding domain (IBD), and the N terminus binds chromatin. LEDGF normally directs integrations to the bodies of expressed genes. Replacing the N terminus of LEDGF with chromatin binding domains (CBDs) from other proteins changes the specificity of HIV-1 DNA integration. We chose two well-characterized CBDs: the plant homeodomain (PHD) finger from ING2 and the chromodomain from heterochromatin binding protein 1alpha (HP1alpha). The ING2 PHD finger binds H3K4me3, a histone mark that is associated with the transcriptional start sites of expressed genes. The HP1alpha chromodomain binds H3K9me2,3, histone marks that are widely distributed throughout the genome. A fusion protein in which the ING2 PHD finger was linked to the LEDGF IBD directed integrations near the start sites of expressed genes. A similar fusion protein in which the HP1alpha chromodomain was linked to the LEDGF IBD directed integrations to sites that differed from both the PHD finger fusion-directed and LEDGF-directed integration sites. The ability to redirect HIV-1 DNA integration may help solve the problems associated with the activation of oncogenes when retroviruses are used in gene therapy.

Expressed antibody repertoires in human cord blood cells: 454 sequencing and IMGT/HighV-QUEST analysis of germline gene usage, junctional diversity, and somatic mutations.

Human cord blood cell-derived IgM antibodies are important for the neonate immune responses and construction of germline-based immunoglobulin libraries. Several previous studies of a relatively small number of sequences found that they exhibit restrictions in the usage of germline genes and in the diversity of the variable heavy chain complementarity determining region 3 compared to adults. To further characterize such restrictions on a larger scale and to compare the early B-cell diversity to adult IgM repertoires, we performed 454 sequencing and IMGT/HighV-QUEST analysis of cord blood IG libraries from two babies and determined germline gene usage, V-D-J rearrangement, VHCDR3 diversity, and somatic mutations to characterize human neonate repertoire. Most of the germline subgroups were identified with frequencies comparable to those present in the adult IgM repertoire except for the IGHV1-2 gene that was preferentially expressed in the cord blood cells. The gene usage diversity contributed to 1,430 unique IGH V-D-J rearrangement patterns while the exonuclease trimming and N region addition at the V-D-J junctions along with gene diversity created a wide range of VHCDR3 with different lengths and sequence variability. We observed a lower degree of somatic mutations in the CDR and framework regions of antibodies from cord blood cells compared to adults. These results provide insights into the characteristics of human cord blood antibody repertoires, which have gene usage diversity and VHCDR3 lengths similar to that of the adult IgM repertoire but differ significantly in some of the gene usages, V-D-J rearrangements, junctional diversity, and somatic mutations.

Detailed DNA methylation profiles of the E-cadherin promoter in the NCI-60 cancer cells.

E-cadherin (E-cad) is a transmembrane adhesion glycoprotein, the expression of which is often reduced in invasive or metastatic tumors. To assess E-cad's distribution among different types of cancer cells, we used bisulfite-sequencing for detailed, base-by-base measurement of CpG methylation in E-cad's promoter region in the NCI-60 cell lines. The mean methylation levels of the cell lines were distributed bimodally, with values pushed toward either the high or low end of the methylation scale. The 38 epithelial cell lines showed substantially lower (28%) mean methylation levels compared with the nonepithelial cell lines (58%). The CpG site at -143 with respect to the transcriptional start was commonly methylated at intermediate levels, even in cell lines with low overall DNA methylation. We also profiled the NCI-60 cell lines using Affymetrix U133 microarrays and found E-cad expression to be correlated with E-cad methylation at highly statistically significant levels. Above a threshold of approximately 20% to 30% mean methylation, the expression of E-cad was effectively silenced. Overall, this study provides a type of detailed analysis of methylation that can also be applied to other cancer-related genes. As has been shown in recent years, DNA methylation status can serve as a biomarker for use in choosing therapy.

Genetic polymorphisms and spontaneous preterm birth.

OBJECTIVE: To examine whether selected genetic polymorphisms in the infant are associated with spontaneous preterm birth (less than 37 weeks) among children with or without later-diagnosed cerebral palsy. METHODS: Exploratory case-control study investigating the relationship of gestational age at delivery to 31 single nucleotide polymorphisms measured in newborn screening bloodspots. Among all 443 children with later-diagnosed cerebral palsy born to white women in South Australia in 1986-1999, 234 were born after spontaneous onset of labor, and 108 of these were preterm (gestational age less than 37 weeks). The comparison group was 549 infants born after spontaneous onset of labor, of whom 147 were preterm. Distributions of genotypic frequencies were examined in preterm compared with term infants with and without cerebral palsy. Genotyping was performed using a Taqman assay. RESULTS: In children without cerebral palsy, preterm birth after spontaneous onset of labor was more frequent in association with a variant of the beta2 adrenergic receptor gene (ADRB2 Q27E, P=.003), inducible nitric oxide synthase (iNOS or NOS2A, P=.042), or thrombomodulin (G127A, P=.006). Among children with cerebral palsy, preterm birth was associated with polymorphisms in genes for endothelial nitric oxide synthase (eNOS -922, P=.012), plasminogen activator inhibitor-2 (P=.015 and .019), and alpha adducin (ADD1, P=.047). CONCLUSION: We confirm previous observations that variants of the beta adrenergic receptor and of nitric oxide synthase are associated with prematurity, and suggest that genetic variants of the placental antifibrinolytic plasminogen activator inhibitor-2, and thrombomodulin and alpha adducin may be contributors to risk of spontaneous preterm birth. LEVEL OF EVIDENCE: II.

Development of multiplex DNA electronic microarrays using a universal adaptor system for detection of single nucleotide polymorphisms.

The NanoChip electronic microarray is designed for the rapid detection of genetic variation in research and clinical diagnosis. We have developed a multiplex electronic microarray assay, specific for single nucleotide polymorphism (SNP) genotyping and mutation detection, using universal adaptor sequences tailed to the 5' end of PCR primers specific to each target. PCR products, amplified by primers directed to the universal adaptor sequence, are immobilized on the microarray either directly or via capture oligonucleotides complementary to the universal adaptor sequence. This simple modification results in a significant increase in fidelity with improved specificity and accuracy. In addition, the multiplexing of genetic variant detection allows increased throughput and significantly reduced cost per assay. This general schema can also be applied to other microarray and macroarray formats.

Genes affected by mouse mammary tumor virus (MMTV) proviral insertions in mouse mammary tumors are deregulated or mutated in primary human mammary tumors.

The accumulation of mutations is a contributing factor in the initiation of premalignant mammary lesions and their progression to malignancy and metastasis. We have used a mouse model in which the carcinogen is the mouse mammary tumor virus (MMTV) which induces clonal premalignant mammary lesions and malignant mammary tumors by insertional mutagenesis. Identification of the genes and signaling pathways affected in MMTV-induced mouse mammary lesions provides a rationale for determining whether genetic alteration of the human orthologues of these genes/pathways may contribute to human breast carcinogenesis. A high-throughput platform for inverse PCR to identify MMTV-host junction fragments and their nucleotide sequences in a large panel of MMTV-induced lesions was developed. Validation of the genes affected by MMTV-insertion was carried out by microarray analysis. Common integration site (CIS) means that the gene was altered by an MMTV proviral insertion in at least two independent lesions arising in different hosts. Three of the new genes identified as CIS for MMTV were assayed for their capability to confer on HC11 mouse mammary epithelial cells the ability for invasion, anchorage independent growth and tumor development in nude mice. Analysis of MMTV induced mammary premalignant hyperplastic outgrowth (HOG) lines and mammary tumors led to the identification of CIS restricted to 35 loci. Within these loci members of the Wnt, Fgf and Rspo gene families plus two linked genes (Npm3 and Ddn) were frequently activated in tumors induced by MMTV. A second group of 15 CIS occur at a low frequency (2-5 observations) in mammary HOGs or tumors. In this latter group the expression of either Phf19 or Sdc2 was shown to increase HC11 cells invasion capability. Foxl1 expression conferred on HC11 cells the capability for anchorage-independent colony formation in soft agar and tumor development in nude mice. The published transcriptome and nucleotide sequence analysis of gene expression in primary human breast tumors was interrogated. Twenty of the human orthologues of MMTV CIS associated genes are deregulated and/or mutated in human breast tumors.

Candidate genes and cerebral palsy: a population-based study.

OBJECTIVE: The objective of this study was to examine whether selected genetic polymorphisms in the infant are associated with later-diagnosed cerebral palsy. METHODS: A population-based case-control study was conducted of 28 single-nucleotide polymorphisms measured in newborn screening blood spots. A total of 413 children with later-diagnosed cerebral palsy were born to white women in South Australia in 1986-1999, and there were 856 control children. Distributions of genotypic frequencies were examined in total cerebral palsy, in gestational age groups, and by types of cerebral palsy and gender. Genotyping was performed by using a TaqMan assay. RESULTS: For inducible nitric-oxide synthase, possession of the T allele was more common in all children with cerebral palsy and for heterozygotes who were born at term. For lymphotoxin alpha, homozygous variant status was associated with risk for cerebral palsy and with spastic hemiplegic or quadriplegic cerebral palsy. Among term infants, heterozygosity for the endothelial protein C receptor single-nucleotide polymorphism was more frequent in children with cerebral palsy. In preterm infants, the variant A allele of interleukin 8 and heterozygosity for the beta-2 adrenergic receptor were associated with cerebral palsy risk. Interleukin 8 heterozygote status was associated with spastic diplegia. Variants of several genes were associated with cerebral palsy in girls but not in boys. CONCLUSIONS: Two of the 28 single-nucleotide polymorphisms examined were associated with all types of spastic cerebral palsy in both gestational age groups and others with cerebral palsy in gestational age or cerebral palsy subgroups. Some of these associations support previous findings. There may be a genetic contribution to cerebral palsy risk, and additional investigation is warranted of genes and gene-environment interactions in cerebral palsy.

Human T-cell leukemia virus type 1 integration target sites in the human genome: comparison with those of other retroviruses.

Retroviral integration into the host genome is not entirely random, and integration site preferences vary among different retroviruses. Human immunodeficiency virus (HIV) prefers to integrate within active genes, whereas murine leukemia virus (MLV) prefers to integrate near transcription start sites and CpG islands. On the other hand, integration of avian sarcoma-leukosis virus (ASLV) shows little preference either for genes, transcription start sites, or CpG islands. While host cellular factors play important roles in target site selection, the viral integrase is probably the major viral determinant. It is reasonable to hypothesize that retroviruses with similar integrases have similar preferences for target site selection. Although integration profiles are well defined for members of the lentivirus, spumaretrovirus, alpharetrovirus, and gammaretrovirus genera, no members of the deltaretroviruses, for example, human T-cell leukemia virus type 1 (HTLV-1), have been evaluated. We have mapped 541 HTLV-1 integration sites in human HeLa cells and show that HTLV-1, like ASLV, does not specifically target transcription units and transcription start sites. Comparing the integration sites of HTLV-1 with those of ASLV, HIV, simian immunodeficiency virus, MLV, and foamy virus, we show that global and local integration site preferences correlate with the sequence/structure of virus-encoded integrases, supporting the idea that integrase is the major determinant of retroviral integration site selection. Our results suggest that the global integration profiles of other retroviruses could be predicted from phylogenetic comparisons of the integrase proteins. Our results show that retroviruses that engender different insertional mutagenesis risks can have similar integration profiles.

Cytokine-independent growth and clonal expansion of a primary human CD8+ T-cell clone following retroviral transduction with the IL-15 gene.

Malignancies arising from retrovirally transduced hematopoietic stem cells have been reported in animal models and human gene therapy trials. Whether mature lymphocytes are susceptible to insertional mutagenesis is unknown. We have characterized a primary human CD8(+) T-cell clone, which exhibited logarithmic ex vivo growth in the absence of exogenous cytokine support for more than 1 year after transduction with a murine leukemia virus-based vector encoding the T-cell growth factor IL-15. Phenotypically, the clone was CD28(-), CD45RA(-), CD45RO(+), and CD62L(-), a profile consistent with effector memory T lymphocytes. After gene transfer with tumor-antigen-specific T-cell receptors, the clone secreted IFN-gamma upon encountering tumor targets, providing further evidence that they derived from mature lymphocytes. Gene-expression analyses revealed no evidence of insertional activation of genes flanking the retroviral insertion sites. The clone exhibited constitutive telomerase activity, and the presence of autocrine loop was suggested by impaired cell proliferation following knockdown of IL-15R alpha expression. The generation of this cell line suggests that nonphysiologic expression of IL-15 can result in the long-term in vitro growth of mature human T lymphocytes. The cytokine-independent growth of this line was a rare event that has not been observed in other IL-15 vector transduction experiments or with any other integrating vector system. It does not appear that the retroviral vector integration sites played a role in the continuous growth of this cell clone, but this remains under investigation.

A comprehensive SNP-based genetic analysis of inbred mouse strains.

Dense genetic maps of mammalian genomes facilitate a variety of biological studies including the mapping of polygenic traits, positional cloning of monogenic traits, mapping of quantitative or qualitative trait loci, marker association, allelic imbalance, speed congenic construction, and evolutionary or phylogenetic comparison. In particular, single nucleotide polymorphisms (SNPs) have proved useful because of their abundance and compatibility with multiple high-throughput technology platforms. SNP genotyping is especially suited for the genetic analysis of model organisms such as the mouse because biallelic markers remain fully informative when used to characterize crosses between inbred strains. Here we report the mapping and genotyping of 673 SNPs (including 519 novel SNPs) in 55 of the most commonly used mouse strains. These data have allowed us to construct a phylogenetic tree that correlates and expands known genealogical relationships and clarifies the origin of strains previously having an uncertain ancestry. All 55 inbred strains are distinguishable genetically using this SNP panel. Our data reveal an uneven SNP distribution consistent with a mosaic pattern of inheritance and provide some insight into the changing dynamics of the physical architecture of the genome. Furthermore, these data represent a valuable resource for the selection of markers and the design of experiments that require the genetic distinction of any pair of mouse inbred strains such as the generation of congenic mice, positional cloning, and the mapping of quantitative or qualitative trait loci.

Common sites of retroviral integration in mouse hematopoietic tumors identified by high-throughput, single nucleotide polymorphism-based mapping and bacterial artificial chromosome hybridization.

Retroviral insertional mutagenesis in mouse hematopoietic tumors provides a powerful cancer gene discovery tool. Here, we describe a high-throughput, single nucleotide polymorphism (SNP)-based method, for mapping retroviral integration sites cloned from mouse tumors, and a bacterial artificial chromosome (BAC) hybridization method, for localizing these retroviral integration sites to common sites of retroviral integration (CISs). Several new CISs were identified, including one CIS that mapped near Notch1, a gene that has been causally associated with human T-cell tumors. This mapping method is applicable to many different species, including ones where few genetic markers and little genomic sequence information are available. It can also be used to map endogenous proviruses.