Relationship of porcine IGF2 imprinting status to DNA methylation at the H19 DMD and the IGF2 DMRs 1 and 2

BACKGROUND: Porcine IGF2 and the H19 genes are imprinted. The IGF2 is paternally expressed, while the H19 gene is maternally expressed. Extensive studies in mice established a boundary model indicating that the H19 differentially methylated domain (DMD) controls, upon binding with the CTCF protein, reciprocal imprinting of the IGF2 and the H19 genes. IGF2 transcription is tissue and development specific involving the use of 4 promoters. In the liver of adult Large White boars IGF2 is expressed from both parental alleles, whereas in skeletal muscle and kidney tissues we observed variable relaxation of IGF2 imprinting. We hypothesized that IGF2 expression from both paternal alleles and relaxation of IGF2 imprinting is reflected in differences in DNA methylation patterns at the H19 DMD and IGF2 differentially methylated regions 1 and 2 (DMR1 and DMR2). RESULTS: Bisulfite sequencing analysis did not show any differences in DNA methylation at the three porcine CTCF binding sites in the H19 DMD between liver, muscle and kidney tissues of adult pigs. A DNA methylation analysis using methyl-sensitive restriction endonuclease SacII and 'hot-stop' PCR gave consistent results with those from the bisulfite sequencing analysis. We found that porcine H19 DMD is distinctly differentially methylated, at least for the region formally confirmed by two SNPs, in liver, skeletal muscle and kidney of foetal, newborn and adult pigs, independent of the combined imprinting status of all IGF2 expressed transcripts. DNA methylation at CpG sites in DMR1 of foetal liver was significantly lower than in the adult liver due to the presence of hypomethylated molecules. An allele specific analysis was performed for IGF2 DMR2 using a SNP in the IGF2 3'-UTR. The maternal IGF2 DMR2 of foetal and newborn liver revealed a higher DNA methylation content compared to the respective paternal allele. CONCLUSIONS: Our results indicate that the IGF2 imprinting status is transcript-specific. Biallelic IGF2 expression in adult porcine liver and relaxation of IGF2 imprinting in porcine muscle were a common feature. These results were consistent with the IGF2 promoter P1 usage in adult liver and IGF2 promoter P2, P3 and P4 usages in muscle. The results showed further that bialellic IGF2 expression in liver and relaxation of imprinting in muscle and kidney were not associated with DNA methylation variation at and around at least one CTCF binding site in H19 DMD. The imprinting status in adult liver, muscle and kidney tissues were also not reflected in the methylation patterns of IGF2 DMRs 1 and 2.

Deep-sequencing identification of the genomic targets of the cytidine deaminase AID and its cofactor RPA in B lymphocytes

The cytidine deaminase AID hypermutates immunoglobulin genes but can also target oncogenes, leading to tumorigenesis. The extent of AID's promiscuity and its predilection for immunoglobulin genes are unknown. We report here that AID interacted broadly with promoter-proximal sequences associated with stalled polymerases and chromatin-activating marks. In contrast, genomic occupancy of replication protein A (RPA), an AID cofactor, was restricted to immunoglobulin genes. The recruitment of RPA to the immunoglobulin loci was facilitated by phosphorylation of AID at Ser38 and Thr140. We propose that stalled polymerases recruit AID, thereby resulting in low frequencies of hypermutation across the B cell genome. Efficient hypermutation and switch recombination required AID phosphorylation and correlated with recruitment of RPA. Our findings provide a rationale for the oncogenic role of AID in B cell malignancy.

Multilocus sequence analysis and rpoB sequencing of Mycobacterium abscessus (sensu lato) strains

Mycobacterium abscessus, Mycobacterium bolletii, and Mycobacterium massiliense (Mycobacterium abscessus sensu lato) are closely related species that currently are identified by the sequencing of the rpoB gene. However, recent studies show that rpoB sequencing alone is insufficient to discriminate between these species, and some authors have questioned their current taxonomic classification. We studied here a large collection of M. abscessus (sensu lato) strains by partial rpoB sequencing (752 bp) and multilocus sequence analysis (MLSA). The final MLSA scheme developed was based on the partial sequences of eight housekeeping genes: argH, cya, glpK, gnd, murC, pgm, pta, and purH. The strains studied included the three type strains (M. abscessus CIP 104536(T), M. massiliense CIP 108297(T), and M. bolletii CIP 108541(T)) and 120 isolates recovered between 1997 and 2007 in France, Germany, Switzerland, and Brazil. The rpoB phylogenetic tree confirmed the existence of three main clusters, each comprising the type strain of one species. However, divergence values between the M. massiliense and M. bolletii clusters all were below 3% and between the M. abscessus and M. massiliense clusters were from 2.66 to 3.59%. The tree produced using the concatenated MLSA gene sequences (4,071 bp) also showed three main clusters, each comprising the type strain of one species. The M. abscessus cluster had a bootstrap value of 100% and was mostly compact. Bootstrap values for the M. massiliense and M. bolletii branches were much lower (71 and 61%, respectively), with the M. massiliense cluster having a fuzzy aspect. Mean (range) divergence values were 2.17% (1.13 to 2.58%) between the M. abscessus and M. massiliense clusters, 2.37% (1.5 to 2.85%) between the M. abscessus and M. bolletii clusters, and 2.28% (0.86 to 2.68%) between the M. massiliense and M. bolletii clusters. Adding the rpoB sequence to the MLSA-concatenated sequence (total sequence, 4,823 bp) had little effect on the clustering of strains. We found 10/120 (8.3%) isolates for which the concatenated MLSA gene sequence and rpoB sequence were discordant (e.g., M. massiliense MLSA sequence and M. abscessus rpoB sequence), suggesting the intergroup lateral transfers of rpoB. In conclusion, our study strongly supports the recent proposal that M. abscessus, M. massiliense, and M. bolletii should constitute a single species. Our findings also indicate that there has been a horizontal transfer of rpoB sequences between these subgroups, precluding the use of rpoB sequencing alone for the accurate identification of the two proposed M. abscessus subspecies.

Ambiguous nucleotide calls from population-based sequencing of HIV-1 are a marker for viral diversity and the age of infection

The time passed since the infection of a human immunodeficiency virus (HIV)-infected individual (the age of infection) is an important but often only poorly known quantity. We assessed whether the fraction of ambiguous nucleotides obtained from bulk sequencing as done for genotypic resistance testing can serve as a proxy of this parameter.

Targeted next generation sequencing as a diagnostic tool in epileptic disorders

Epilepsies have a highly heterogeneous background with a strong genetic contribution. The variety of unspecific and overlapping syndromic and nonsyndromic phenotypes often hampers a clear clinical diagnosis and prevents straightforward genetic testing. Knowing the genetic basis of a patient's epilepsy can be valuable not only for diagnosis but also for guiding treatment and estimating recurrence risks.

Identification of the bovine Arachnomelia mutation by massively parallel sequencing implicates sulfite oxidase (SUOX) in bone development

Arachnomelia is a monogenic recessive defect of skeletal development in cattle. The causative mutation was previously mapped to a approximately 7 Mb interval on chromosome 5. Here we show that array-based sequence capture and massively parallel sequencing technology, combined with the typical family structure in livestock populations, facilitates the identification of the causative mutation. We re-sequenced the entire critical interval in a healthy partially inbred cow carrying one copy of the critical chromosome segment in its ancestral state and one copy of the same segment with the arachnomelia mutation, and we detected a single heterozygous position. The genetic makeup of several partially inbred cattle provides extremely strong support for the causality of this mutation. The mutation represents a single base insertion leading to a premature stop codon in the coding sequence of the SUOX gene and is perfectly associated with the arachnomelia phenotype. Our findings suggest an important role for sulfite oxidase in bone development.

Generation and analysis of a mouse intestinal metatranscriptome through Illumina based RNA-sequencing

With the advent of high through-put sequencing (HTS), the emerging science of metagenomics is transforming our understanding of the relationships of microbial communities with their environments. While metagenomics aims to catalogue the genes present in a sample through assessing which genes are actively expressed, metatranscriptomics can provide a mechanistic understanding of community inter-relationships. To achieve these goals, several challenges need to be addressed from sample preparation to sequence processing, statistical analysis and functional annotation. Here we use an inbred non-obese diabetic (NOD) mouse model in which germ-free animals were colonized with a defined mixture of eight commensal bacteria, to explore methods of RNA extraction and to develop a pipeline for the generation and analysis of metatranscriptomic data. Applying the Illumina HTS platform, we sequenced 12 NOD cecal samples prepared using multiple RNA-extraction protocols. The absence of a complete set of reference genomes necessitated a peptide-based search strategy. Up to 16% of sequence reads could be matched to a known bacterial gene. Phylogenetic analysis of the mapped ORFs revealed a distribution consistent with ribosomal RNA, the majority from Bacteroides or Clostridium species. To place these HTS data within a systems context, we mapped the relative abundance of corresponding Escherichia coli homologs onto metabolic and protein-protein interaction networks. These maps identified bacterial processes with components that were well-represented in the datasets. In summary this study highlights the potential of exploiting the economy of HTS platforms for metatranscriptomics.

Real-time quantitative broad-range PCR assay for detection of the 16S rRNA gene followed by sequencing for species identification

Here we determined the analytical sensitivities of broad-range real-time PCR-based assays employing one of three different genomic DNA extraction protocols in combination with one of three different primer pairs targeting the 16S rRNA gene to detect a panel of 22 bacterial species. DNA extraction protocol III, using lysozyme, lysostaphin, and proteinase K, followed by PCR with the primer pair Bak11W/Bak2, giving amplicons of 796 bp in length, showed the best overall sensitivity, detecting DNA of 82% of the strains investigated at concentrations of < or =10(2) CFU in water per reaction. DNA extraction protocols I and II, using less enzyme treatment, combined with other primer pairs giving shorter amplicons of 466 bp and 342 or 346 bp, respectively, were slightly more sensitive for the detection of gram-negative but less sensitive for the detection of gram-positive bacteria. The obstacle of detecting background DNA in blood samples spiked with bacteria was circumvented by introducing a broad-range hybridization probe, and this preserved the minimal detection limits observed in samples devoid of blood. Finally, sequencing of the amplicons generated using the primer pair Bak11W/Bak2 allowed species identification of the detected bacterial DNA. Thus, broad-spectrum PCR targeting the 16S rRNA gene in the quantitative real-time format can achieve an analytical sensitivity of 1 to 10 CFU per reaction in water, avoid detection of background DNA with the introduction of a broad-range probe, and generate amplicons that allow species identification of the detected bacterial DNA by sequencing. These prerequisites are important for its application to blood-containing patient samples.

Evaluation of the colorimetric VITEK 2 card for identification of gram-negative nonfermentative rods: comparison to 16S rRNA gene sequencing

Ninety strains of a collection of well-identified clinical isolates of gram-negative nonfermentative rods collected over a period of 5 years were evaluated using the new colorimetric VITEK 2 card. The VITEK 2 colorimetric system identified 53 (59%) of the isolates to the species level and 9 (10%) to the genus level; 28 (31%) isolates were misidentified. An algorithm combining the colorimetric VITEK 2 card and 16S rRNA gene sequencing for adequate identification of gram-negative nonfermentative rods was developed. According to this algorithm, any identification by the colorimetric VITEK 2 card other than Achromobacter xylosoxidans, Acinetobacter sp., Burkholderia cepacia complex, Pseudomonas aeruginosa, and Stenotrophomonas maltophilia should be subjected to 16S rRNA gene sequencing when accurate identification of nonfermentative rods is of concern.

The tumor suppressor gene hypermethylated in cancer 1 is transcriptionally regulated by E2F1

The Hypermethylated in Cancer 1 (HIC1) gene encodes a zinc finger transcriptional repressor that cooperates with p53 to suppress cancer development. We and others recently showed that HIC1 is a transcriptional target of p53. To identify additional transcriptional regulators of HIC1, we screened a set of transcription factors for regulation of a human HIC1 promoter reporter. We found that E2F1 strongly activates the full-length HIC1 promoter reporter. Promoter deletions and mutations identified two E2F responsive elements in the HIC1 core promoter region. Moreover, in vivo binding of E2F1 to the HIC1 promoter was shown by chromatin immunoprecipitation assays in human TIG3 fibroblasts expressing tamoxifen-activated E2F1. In agreement, activation of E2F1 in TIG3-E2F1 cells markedly increased HIC1 expression. Interestingly, expression of E2F1 in the p53(-/-) hepatocellular carcinoma cell line Hep3B led to an increase of endogenous HIC1 mRNA, although bisulfite genomic sequencing of the HIC1 promoter revealed that the region bearing the two E2F1 binding sites is hypermethylated. In addition, endogenous E2F1 induced by etoposide treatment bound to the HIC1 promoter. Moreover, inhibition of E2F1 strongly reduced the expression of etoposide-induced HIC1. In conclusion, we identified HIC1 as novel E2F1 transcriptional target in DNA damage responses.

Mutation screening of the medium-chain acyl-CoA dehydrogenase (MCAD) and the ornithine transcarbamylase (OTC) genes by multiplex PCR amplification and sequencing

BACKGROUND: Sequencing based mutation screening assays of genes encompassing large numbers of exons could be substantially optimized by multiplex PCR, which enables simultaneous amplification of many targets in one reaction. In the present study, a multiplex PCR protocol originally developed for fragment analysis was evaluated for sequencing based mutation screening of the ornithine transcarbamylase (OTC) and the medium-chain acyl-CoA dehydrogenase (MCAD) genes. METHODS: Single exon and multiplex PCR protocols were applied to generate PCR templates for subsequent DNA sequencing of all exons of the OTC and the MCAD genes. For each PCR protocol and using the same DNA samples, 66 OTC and 98 MCAD sequence reads were generated. The sequences derived from the two different PCR methods were compared at the level of individual signal-to-noise ratios of the four bases and the proportion of high-quality base-signals. RESULTS: The single exon and the multiplex PCR protocol gave qualitatively comparable results for the two genes. CONCLUSIONS: Many existing sequencing based mutation analysis protocols may be easily optimized with the proposed method, since the multiplex PCR protocol was successfully applied without any re-design of the PCR primers and other optimization steps for generating sequencing templates for the OTC and MCAD genes, respectively.