Next-generation sequencing of HIV-1 RNA genomes: determination of error rates and minimizing artificial recombination.

Next-generation sequencing (NGS) is a valuable tool for the detection and quantification of HIV-1 variants in vivo. However, these technologies require detailed characterization and control of artificially induced errors to be applicable for accurate haplotype reconstruction. To investigate the occurrence of substitutions, insertions, and deletions at the individual steps of RT-PCR and NGS, 454 pyrosequencing was performed on amplified and non-amplified HIV-1 genomes. Artificial recombination was explored by mixing five different HIV-1 clonal strains (5-virus-mix) and applying different RT-PCR conditions followed by 454 pyrosequencing. Error rates ranged from 0.04-0.66% and were similar in amplified and non-amplified samples. Discrepancies were observed between forward and reverse reads, indicating that most errors were introduced during the pyrosequencing step. Using the 5-virus-mix, non-optimized, standard RT-PCR conditions introduced artificial recombinants in a fraction of at least 30% of the reads that subsequently led to an underestimation of true haplotype frequencies. We minimized the fraction of recombinants down to 0.9-2.6% by optimized, artifact-reducing RT-PCR conditions. This approach enabled correct haplotype reconstruction and frequency estimations consistent with reference data obtained by single genome amplification. RT-PCR conditions are crucial for correct frequency estimation and analysis of haplotypes in heterogeneous virus populations. We developed an RT-PCR procedure to generate NGS data useful for reliable haplotype reconstruction and quantification.

Isolation and analysis of regulatory and structural mutations in the {\it recA\/} gene of {\it Escherichia coli}

The recA gene is essential for homologous recombination and for inducible DNA repair in Escherichia coli. The level of recA expression is important for these functions. The growth defect of a lambda phage carrying a recA-lacZ fusion was used to select mutations that reduced recA expression. Nine of these mutations were single base changes in the recA promoter; each reduced both induced and basal (repressed) levels of expression, indicating that only one promoter is used under both circumstances. Deletion analysis of the promoter region and S1 mapping of transcripts confirmed that there is only one promoter responsible for both basal and induced expression. Some of the mutants, however, displayed a ratio of induced to repressed expression that was much lower than wild-type. For one of these mutants (recA1270) LexA binding studies showed that this was not due to a change in the affinity of LexA repressor for the operator site. The extent of binding of RNA polymerase to this mutant promoter, however, was much reduced, and the complexes formed were qualitatively different. Further binding experiments provided some evidence that LexA does not block RNA polymerase binding to the recA promoter, but inhibits a later step in initiation. Behavior of the mutants with altered induction ratios could be explained if LexA binding to the operator actually increases RNA polymerase binding to the promoter in a closed complex compensating for defects in polymerase binding caused by the mutations.^ In a study of mutations in the recA structural gene, site-directed mutagenesis was used to replace cysteine codons at positions 90, 116, and 129 with a number of different codons. In vivo analysis of the replacements showed that none of the cysteines is absolutely essential and that they do not have a direct role as catalysts in ATP hydrolysis. Some amino acid substitutions abolished all RecA functions, while a few resulted in partial or altered function. Amino acids at positions 90 and 129 tended to affect all functions equally, while the amino acid at position 116 appeared to have a particular effect on the protease activity of the protein. ^

Phenotypic and genotypic analysis of mouse hepatitis virus (JHM) complementation group A temperature-sensitive mutants

Temperature sensitive (ts) mutant viruses have helped elucidate replication processes in many viral systems. Several panels of replication-defective ts mutants in which viral RNA synthesis is abolished at the nonpermissive temperature (RNA$\sp{-})$ have been isolated for Mouse Hepatitis Virus, MHV (Robb et al., 1979; Koolen et al., 1983; Martin et al., 1988; Schaad et al., 1990). However, no one had investigated genetic or phenotypic relationships between these different mutant panels. In order to determine how the panel of MHV-JHM RNA$\sp{-}$ ts mutants (Robb et al., 1979) were genetically related to other described MHV RNA$\sp{-}$ ts mutants, the MHV-JHM mutants were tested for complementation with representatives from two different sets of MHV-A59 ts mutants (Koolen et al., 1983; Schaad et al., 1990). The three ts mutant panels together were found to comprise eight genetically distinct complementation groups. Of these eight complementation groups, three complementation classes are unique to their particular mutant panel; genetically equivalent mutants were not observed within the other two mutant panels. Two complementation groups were common to all three mutant panels. The three remaining complementation groups overlapped two of the three mutant sets. Mutants MHV-JHM tsA204 and MHV-A59 ts261 were shown to be within one of these overlapping complementation groups. The phenotype of the MHV-JHM mutants within this complementation class has been previously characterized (Leibowitz et al., 1982; Leibowitz et al, 1990). When these mutants were grown at the permissive temperature, then shifted up to the nonpermissive temperature at the start of RNA synthesis, genome-length RNA and leader RNA fragments accumulated, but no subgenomic mRNA was synthesized. MHV-A59 ts261 produced leader RNA fragments identical to those observed with MHV-JHM tsA204. Thus, these two MHV RNA$\sp{-}$ ts mutants that were genetically equivalent by complementation testing were phenotypically similar as well. Recombination frequencies obtained from crosses of MHV-A59 ts261 with several of the gene 1 MHV-A59 mutants indicated that the causal mutation(s) of MHV-A59 ts261 was located near the overlapping junction of ORF1a and ORF1b, in the 3$\sp\prime$ end of ORF1a, or the 5$\sp\prime$ end of ORF1b. Sequence analysis of this junction and 1400 nucleotides into the 5$\sp\prime$ end of ORF1b of MHV-A59 ts261 revealed one nucleotide change from the wildtype MHV-A59. This substitution at nucleotide 13,598 (A to G) was a silent mutation in the ORF1a reading frame, but resulted in an amino acid change in ORF1b gene product (I to V). This amino acid change would be expressed only in the readthrough translation product produced upon successful ribosome frameshifting. A revertant of MHV-A59 ts261 (R2) also retained this guanidine residue, but had a second substitution at nucleotide 14,475 in ORF1b. This mutation results in the substitution of valine for an isoleucine.^ The data presented here suggest that the mutation in MHV-A59 ts261 (nucleotide 13,598) would be responsible for the MHV-JHM complementation group A phenotype. A second-site reversion at nucleotide 14,475 may correct this defect in the revertant. Sequencing of gene 1 immediately upstream of nucleotide 13,296 and downstream of nucleotide 15,010 must be conducted to test this hypothesis. ^

Traditional linkage analysis in admixed families

Currently there is no general method to study the impact of population admixture within families on the assumptions of random mating and consequently, Hardy-Weinberg equilibrium (HWE) and linkage equilibrium (LE) and on the inference obtained from traditional linkage analysis. ^ First, through simulation, the effect of admixture of two populations on the log of the odds (LOD) score was assessed, using Prostate Cancer as the typical disease model. Comparisons between simulated mixed and homogeneous families were performed. LOD scores under both models of admixture (within families and within a data set of homogeneous families) were closest to the homogeneous family scores of the population having the highest mixing proportion. Random sampling of families or ascertainment of families with disease affection status did not affect this observation, nor did the mode of inheritance (dominant/recessive) or sample size. ^ Second, after establishing the effect of admixture on the LOD score and inference for linkage, the presence of induced disequilibria by population admixture within families was studied and an adjustment procedure was developed. The adjustment did not force all disequilibria to disappear but because the families were adjusted for the population admixture, those replicates where the disequilibria exist are no longer affected by the disequilibria in terms of maximization for linkage. Furthermore, the adjustment was able to exclude uninformative families or families that had such a high departure from HWE and/or LE that their LOD scores were not reliable. ^ Together these observations imply that the presence of families of mixed population ancestry impacts linkage analysis in terms of the LOD score and the estimate of the recombination fraction. ^

Whole genome analysis of Fusobacterium nucleatum, Rickettsia typhi and Francisella tularensis

The genomes of Fusobacterium nucleatum subspecies polymorphum strain ATCC 10953, Rickettsia typhi strain Wilmington, and Francisella tularensis subspecies holarctica strain OSU18 were sequenced, annotated, and analyzed. Each genome was then compared to the sequenced genomes of closely related bacteria. The genome of F. nucleatum ATCC 10953 was compared to two additional F. nucleatum subspecies, subspecies nucleatum and subspecies vincentii. This analysis revealed substantial evidence of horizontal gene transfer along with considerable genetic diversity within the species of F. nucleatum. R. typhi was compared to R. prowazekii and R. conorii. This analysis uncovered a hotspot for chromosomal rearrangements in the Spotted Fever Group but not the Typhus Group Rickettsia and revealed the close genetic relationship between the Typhus Group rickettsial species. F. tularensis OSU18 was compared to two additional F. tularensis strains. These comparisons uncovered significant chromosomal rearrangements between F. tularensis subspecies due to recombination between insertion sequence elements. ^

Functional analysis of BLAP75, a BLM-associated protein

Bloom syndrome (BS) is an autosomal recessive disorder characterized by dwarfism, immunodeficiency, impaired fertility, and most importantly, early development of a broad range of cancers. The hallmark of BS cells is hyper-recombination, characterized by a drastically elevated frequency of sister chromatid exchange (SCE). BLM, the gene mutated in BS, encodes a DNA helicase of the RecQ protein family. BLM is thought to participate in several DNA transactions and to interact with many proteins involved in DNA replication, recombination, and repair. However, the precise function of BLM and the BLM-dependent anti-tumor mechanism remain obscure. ^ A novel protein, BLAP75 (BLM-associated polypeptide, 75KD), was identified to form an evolutionarily conserved complex with BLM and DNA topoisomerase IIIα (Topo IIIα). Our work demonstrates that loss of BLAP75 destabilized BLM and Topo IIIα proteins. BLAP75 colocalized with BLM in subnuclear foci in response to DNA damage and the recruitment of BLM to these foci was BLAP75-dependent. Moreover, depletion of BLAP75 by siRNA resulted in an elevated SCE rate similar to cells depleted of BLM by siRNA. In addition, RNAi-mediated silencing of BLAP75 greatly diminished cell viability. This cellular deficiency was rescued by expression of wild type BLAP75 but not BLAP75 with mutated conserved domain III, which abrogated the interaction between BLAP75, BLM and Topo IIIα, suggesting that the integrity of BLM-Topo IIIα-BLAP75 complex might be critical for cell survival. Finally, I found that BLAP75 was phosphorylated during mitosis and upon various DNA-damaging agents, implying that BLAP75 might also function in mitosis and DNA damage response. ^ Taken together, this study has defined BLAP75 as an integral component of the BLM complex to maintain genome stability. Our findings provide insights into the molecular mechanisms of the BLM helicase pathway and tumorigenesis process associated with these mechanisms. ^

The role of ATM in the cell cycle control of V(D)J recombination

Lymphocyte development requires the assembly of diversified antigen receptor complexes generated by the genetically programmed V(D)J recombination event. Because germline DNA is cut, introducing potentially dangerous double-stranded breaks (DSBs) and rearranged prior to repair, its activity is limited to the non-cycling stages of the cell cycle, G0/G1. The potential involvement of a key mediator, Ataxia Telangiectasia Mutated or ATM, in the DNA damage response (DDR) and cell cycle checkpoints has been implicated in recombination, but its role is not fully understood. Thymic lymphomas from ATM deficient mice contain clonal chromosomal translocations involving the T-cell antigen receptor (TCR). A previous report found ATM and its downstream target p53 associated with V(D)J intermediates, suggesting the DDR senses recombination. In this study, we sought to understand the role of ATM in V(D)J recombination. Developing thymocytes from ATM deficient mice were analyzed according to the cell cycle to detect V(D)J intermediates. Examination of all TCR loci in the non-cycling (G0/G1) and cycling (S/G2/M) fractions revealed the persistence of intermediates in ATM deficient thymocytes, contrary to the wild-type in which intermediates are found only during G0/G1. Further analysis found no defect in end-joining of intermediates, nor were they detected in developed T-cells. Based upon the presence of persisting intermediates, the recombination initiating nuclease Rag-2 was examined; strict regulation limits it to G 0/G1. Rag-2 regulation was not affected by an ATM deficiency as Rag-2 expression remained contained within G0/G 1, indicating recombination is not continuous. To determine if an ATM deficiency affects recognition of V(D)J breaks, sites of recombination identified by a TCR locus or Rag expression were analyzed according to co-localization with a DDR factor phosphorylated immediately after DNA damage, phosphorylated H2AX (γH2AX). No differences in co-localization were found between the wild-type and ATM deficiency, demonstrating ATM deficient lymphocytes retain the ability to recognize DSBs. Together, these results suggest ATM is necessary in the cell cycle regulation of recombination but not essential for the identification of V(D)J breaks. ATM ensures the containment of intermediates within G0/G1 and maintains genomic stability of developing lymphocytes, emphasizing its fundamental role in preventing tumorigenesis.^

A systematic analysis of INO80 in DNA replication

ATP-dependent chromatin remodeling has been shown to be critical for transcription and DNA repair. However, the involvement of ATP-dependent chromatin remodeling in DNA replication remains poorly defined. Interestingly, we found that the INO80 chromatin-remodeling complex is directly involved in the DNA damage tolerance pathways activated during DNA replication. DNA damage tolerance is important for genomic stability and is controlled by formation of either mono-ubiquitinated or multi-ubiquitinated PCNA, which respectively induce error prone or error-free replication bypass of the lesions. In addition, homologous recombination (HR) mediated by the Rad51 pathway is also involved in the DNA damage tolerance pathways. ^ We found that INO80 is specifically recruited to replication origins during S phase in a genome-wide fashion. In addition, DNA combing analysis shows INO80 is required for the resumption of replication at stalled forks induced by methyl methane-sulfonate (MMS). Mechanistically, we find that INO80 is required for PCNA ubiquitination as well as for Rad51 mediated processing of replication forks after MMS treatment. Furthermore, chromatin immunoprecipitation at specific ARSs indicates INO80 is necessary for Rad18 and Rad51 recruitment to replication forks after MMS treatment. Moreover, 2D gel analysis shows INO80 is necessary to process Rad51 mediated intermediates at impeded replication forks. ^ In conclusion, our findings establish a novel role of a chromatin-remodeling complex in DNA damage tolerance pathways and suggest that chromatin remodeling is fundamentally important to ensure faithful replication of DNA and genome stability in eukaryotes. ^

Sub-bandgap spectral photo-response analysis of Ti supersaturated Si

We have analyzed the increase of the sheet conductance (ΔG□) under spectral illumination in high dose Ti implanted Si samples subsequently processed by pulsed-laser melting. Samples with Ti concentration clearly above the insulator-metal transition limit show a remarkably high ΔG□, even higher than that measured in a silicon reference sample. This increase in the ΔG□ magnitude is contrary to the classic understanding of recombination centers action and supports the lifetime recovery predicted for concentrations of deep levels above the insulator-metal transition.

Direct isolation of human BRCA2 gene by transformation-associated recombination in yeast

Mutant forms of the BRCA2 gene contribute significantly to hereditary breast cancer. Isolation of the normal and mutant forms of the BRCA2 gene with its natural promoter would greatly facilitate analysis of the gene and its contribution to breast cancer. We have accomplished the direct isolation of the 90-kb gene from total human DNA by transformation-associated recombination in yeast using a small amount of 5′ and 3′ BRCA2 sequence information. Because the entire isolation procedure of a single chromosomal gene could be accomplished in approximately 2 weeks, the transformation-associated recombination cloning approach is readily applicable to studies of chromosome alterations and human genetic diseases.

Recombination-dependent deletion formation in mammalian cells deficient in the nucleotide excision repair gene ERCC1

Nucleotide excision repair proteins have been implicated in genetic recombination by experiments in Saccharomyces cerevisiae and Drosophila melanogaster, but their role, if any, in mammalian cells is undefined. To investigate the role of the nucleotide excision repair gene ERCC1, the hamster homologue to the S. cerevisiae RAD10 gene, we disabled the gene by targeted knockout. Partial tandem duplications of the adenine phosphoribosyltransferase (APRT) gene then were constructed at the endogenous APRT locus in ERCC1− and ERCC1+ cells. To detect the full spectrum of gene-altering events, we used a loss-of-function assay in which the parental APRT+ tandem duplication could give rise to APRT− cells by homologous recombination, gene rearrangement, or point mutation. Measurement of rates and analysis of individual APRT− products indicated that gene rearrangements (principally deletions) were increased at least 50-fold, whereas homologous recombination was affected little. The formation of deletions is not caused by a general effect of the ERCC1 deficiency on gene stability, because ERCC1− cell lines with a single wild-type copy of the APRT gene yielded no increase in deletions. Thus, deletion formation is dependent on the tandem duplication, and presumably the process of homologous recombination. Recombination-dependent deletion formation in ERCC1− cells is supported by a significant decrease in a particular class of crossover products that are thought to arise by repair of a heteroduplex intermediate in recombination. We suggest that the ERCC1 gene product in mammalian cells is involved in the processing of heteroduplex intermediates in recombination and that the misprocessed intermediates in ERCC1− cells are repaired by illegitimate recombination.

λ Rap protein is a structure-specific endonuclease involved in phage recombination

Bacteriophage λ encodes a number of genes involved in the recombinational repair of DNA double-strand breaks. The product of one of these genes, rap, has been purified. Truncated Rap proteins that copurify with the full-length form are derived, at least in part, from a ρ-dependent transcription terminator located within its coding sequence. Full-length and certain truncated Rap polypeptides bind preferentially to branched DNA substrates, including synthetic Holliday junctions and D-loops. In the presence of manganese ions, Rap acts as an endonuclease that cleaves at the branch point of Holliday and D-loop substrates. It shows no obvious sequence preference or symmetry of cleavage on a Holliday junction. The biochemical analysis of Rap gives an insight into how recombinants could be generated by the nicking of a D-loop without the formation of a classical Holliday junction.

Tid1/Rdh54 promotes colocalization of Rad51 and Dmc1 during meiotic recombination

Two RecA homologs, Rad51 and Dmc1, assemble as cytologically visible complexes (foci) at the same sites on meiotic chromosomes. Time course analysis confirms that co-foci appear and disappear as the single predominant form. A large fraction of co-foci are eliminated in a red1 mutant, which is expected as a characteristic of the interhomolog-specific recombination pathway. Previous studies suggested that normal Dmc1 loading depends on Rad51. We show here that a mutation in TID1/RDH54, encoding a RAD54 homolog, reduces Rad51-Dmc1 colocalization relative to WT. A rad54 mutation, in contrast, has relatively little effect on RecA homolog foci except when strains also contain a tid1/rdh54 mutation. The role of Tid1/Rdh54 in coordinating RecA homolog assembly may be very direct, because Tid1/Rdh54 is known to physically bind both Dmc1 and Rad51. Also, Dmc1 foci appear early in a tid1/rdh54 mutant. Thus, Tid1 may normally act with Rad51 to promote ordered RecA homolog assembly by blocking Dmc1 until Rad51 is present. Finally, whereas double-staining foci predominate in WT nuclei, a subset of nuclei with expanded chromatin exhibit individual Rad51 and Dmc1 foci side-by-side, suggesting that a Rad51 homo-oligomer and a Dmc1 homo-oligomer assemble next to one another at the site of a single double-strand break (DSB) recombination intermediate.

Does recombination improve selection on codon usage? Lessons from nematode and fly complete genomes

Understanding the factors responsible for variations in mutation patterns and selection efficacy along chromosomes is a prerequisite for deciphering genome sequences. Population genetics models predict a positive correlation between the efficacy of selection at a given locus and the local rate of recombination because of Hill–Robertson effects. Codon usage is considered one of the most striking examples that support this prediction at the molecular level. In a wide range of species including Caenorhabditis elegans and Drosophila melanogaster, codon usage is essentially shaped by selection acting for translational efficiency. Codon usage bias correlates positively with recombination rate in Drosophila, apparently supporting the hypothesis that selection on codon usage is improved by recombination. Here we present an exhaustive analysis of codon usage in C. elegans and D. melanogaster complete genomes. We show that in both genomes there is a positive correlation between recombination rate and the frequency of optimal codons. However, we demonstrate that in both species, this effect is due to a mutational bias toward G and C bases in regions of high recombination rate, possibly as a direct consequence of the recombination process. The correlation between codon usage bias and recombination rate in these species appears to be essentially determined by recombination-dependent mutational patterns, rather than selective effects. This result highlights that it is necessary to take into account the mutagenic effect of recombination to understand the evolutionary role and impact of recombination.

Ligand-activated site-specific recombination in mice.

Current mouse gene targeting technology is unable to introduce somatic mutations at a chosen time and/or in a given tissue. We report here that conditional site-specific recombination can be achieved in mice using a new version of the Cre/lox system. The Cre recombinase has been fused to a mutated ligand-binding domain of the human estrogen receptor (ER) resulting in a tamoxifen-dependent Cre recombinase, Cre-ERT, which is activated by tamoxifen, but not by estradiol. Transgenic mice were generated expressing Cre-ERT under the control of a cytomegalovirus promoter. We show that excision of a chromosomally integrated gene flanked by loxP sites can be induced by administration of tamoxifen to these transgenic mice, whereas no excision could be detected in untreated animals. This conditional site-specific recombination system should allow the analysis of knockout phenotypes that cannot be addressed by conventional gene targeting.