Accurate sodium bisulfite sequencing in plants.

DNA cytosine methylation is a conserved epigenetic modification frequently correlating with transcriptional silencing in a wide variety of eukaryotic organisms. Sodium bisulfite treatment of DNA converts unmethylated cytosine to uracil, while 5-methylated cytosine is protected. We describe techniques that ensure reliable sequencing data following sodium bisulfite conversion and to avoid common pitfalls such as amplification of unconverted DNA and inclusion of sibling clones.

A whole-genome massively parallel sequencing analysis of BRCA1 mutant oestrogen receptor-negative and -positive breast cancers

BRCA1 encodes a tumour suppressor protein that plays pivotal roles in homologous recombination (HR) DNA repair, cell-cycle checkpoints, and transcriptional regulation. BRCA1 germline mutations confer a high risk of early-onset breast and ovarian cancer. In more than 80% of cases, tumours arising in BRCA1 germline mutation carriers are oestrogen receptor (ER)-negative; however, up to 15% are ER-positive. It has been suggested that BRCA1 ER-positive breast cancers constitute sporadic cancers arising in the context of a BRCA1 germline mutation rather than being causally related to BRCA1 loss-of-function. Whole-genome massively parallel sequencing of ER-positive and ER-negative BRCA1 breast cancers, and their respective germline DNAs, was used to characterize the genetic landscape of BRCA1 cancers at base-pair resolution. Only BRCA1 germline mutations, somatic loss of the wild-type allele, and TP53 somatic mutations were recurrently found in the index cases. BRCA1 breast cancers displayed a mutational signature consistent with that caused by lack of HR DNA repair in both ER-positive and ER-negative cases. Sequencing analysis of independent cohorts of hereditary BRCA1 and sporadic non-BRCA1 breast cancers for the presence of recurrent pathogenic mutations and/or homozygous deletions found in the index cases revealed that DAPK3, TMEM135, KIAA1797, PDE4D, and GATA4 are potential additional drivers of breast cancers. This study demonstrates that BRCA1 pathogenic germline mutations coupled with somatic loss of the wild-type allele are not sufficient for hereditary breast cancers to display an ER-negative phenotype, and has led to the identification of three potential novel breast cancer genes (ie DAPK3, TMEM135, and GATA4).

Construction of bovine adenovirus type 3 E1 and E3, substitution plasmids and the sequencing analysis of DNA pol and pTP /

The relative ease to concentrate and purify adenoviruses, their well characterized mid-sized genome, and the ability to delete non-essential regions from their genome to accommodate foreign gene, made adenoviruses a suitable candidate for the construction of vectors. The use of adenoviral vectors in gene therapy, vaccination, and as a general vector system for expressing foreign genes have been documented for some time. In this study, the objective was to rescue a BAV3 E1 or E3 recombinant vector carrying the kanamycin resistant gene, a dominant selectable marker with useful applications in studying vectored gene expression in mammalian cells. To accomplish the objective of this study, more information about BAV3 DNA sequences was required in order to make the manipulation of the virus genome accessible. Therefore, sequencing of the BAV3 genome from 1 1 .7% to 30.8% was carried out. Analysis of the determined sequences revealed the primary structure of important viral gene products coded by E2 including BAV3 DNA pol and precursor to terminal protein. Comparative analysis of these proteins with their counterparts from human and non human adenoviruses revealed important insights as to the evolutionary lineage of BAV3. In order to insert the kanamycin resistance gene in either E1 or E3, it was necessary to delete BAV3 sequences to accommodate the foreign gene so as not to exceed the limit of the packaging capacity of the virus. To construct a recombinant BAV3 in which a foreign gene was inserted in the deleted E1 region, an E1 shuttle vector was constructed. This involved the deletion from the viral sequences a region between 1.3% to 9% and inserting the kanamycin resistance gene to replace the deletion. The E1 shuttle vector contained the left (0%- 53.9%) segment of the genome and was expected to generate BAV3 recombinants that can be grown and propagated in cells that can complement the missing E1 functions. To construct a similar shuttle vector for E3 deletion, DNA sequences extending from 78.9% to 82.5% (1281 bp) were deleted from within the E3 region that had been cloned into a plasmid vector. The deleted region corresponds to those that have been shown to be non-essential for viral replication in cell culture. The resulting plasmid was used to construct another recombinant plasmid with BAV3 DNA sequences extending from 37.1% to 100% and with a deletion of E3 sequences that were replaced by kanamycin resistance gene. This shuttle plasmid was used in cotransfections with digested viral DNA in an attempt to rescue a recombinant BAV3 carrying the kanamycin resistance gene to replace the deleted E3. In spite of repeated attempts of transfection, El or E3 recombinant BAV3 were not isolated. It seems that other approaches should be applied to make a final conclusion on BAV3 infectivity.

Suitability of DNA extracted from archival specimens of fruit-eating bats of the genus Artibeus (Chiroptera, Phyllostomidae) for polymerase chain reaction and sequencing analysis

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Bacterial diversity of soil under eucalyptus assessed by 16S rDNA sequencing analysis

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Whole genome sequencing analysis of Plasmodium vivax using whole genome capture

Background: Malaria caused by Plasmodium vivax is an experimentally neglected severe disease with a substantial burden on human health. Because of technical limitations, little is known about the biology of this important human pathogen. Whole genome analysis methods on patient-derived material are thus likely to have a substantial impact on our understanding of P. vivax pathogenesis and epidemiology. For example, it will allow study of the evolution and population biology of the parasite, allow parasite transmission patterns to be characterized, and may facilitate the identification of new drug resistance genes. Because parasitemias are typically low and the parasite cannot be readily cultured, on-site leukocyte depletion of blood samples is typically needed to remove human DNA that may be 1000X more abundant than parasite DNA. These features have precluded the analysis of archived blood samples and require the presence of laboratories in close proximity to the collection of field samples for optimal pre-cryopreservation sample preparation. Results: Here we show that in-solution hybridization capture can be used to extract P. vivax DNA from human contaminating DNA in the laboratory without the need for on-site leukocyte filtration. Using a whole genome capture method, we were able to enrich P. vivax DNA from bulk genomic DNA from less than 0.5% to a median of 55% (range 20%-80%). This level of enrichment allows for efficient analysis of the samples by whole genome sequencing and does not introduce any gross biases into the data. With this method, we obtained greater than 5X coverage across 93% of the P. vivax genome for four P. vivax strains from Iquitos, Peru, which is similar to our results using leukocyte filtration (greater than 5X coverage across 96% of the genome). Conclusion: The whole genome capture technique will enable more efficient whole genome analysis of P. vivax from a larger geographic region and from valuable archived sample collections.

Comparative genomic sequencing analysis of a region in human chromosome 11p15.3, mouse chromosome 7 and analysis of the novel STK33, Stk33 gene

The comparative genomic sequence analysis of a region in human chromosome 11p15.3 and its homologous segment in mouse chromosome 7 between ST5 and LMO1 genes has been performed. 158,201 bases were sequenced in the mouse and compared with the syntenic region in human, partially available in the public databases. The analysed region exhibits the typical eukaryotic genomic structure and compared with the close neighbouring regions, strikingly reflexes the mosaic pattern distribution of (G+C) and repeats content despites its relative short size. Within this region the novel gene STK33 was discovered (Stk33 in the mouse), that codes for a serine/threonine kinase. The finding of this gene constitutes an excellent example of the strength of the comparative sequencing approach. Poor gene-predictions in the mouse genomic sequence were corrected and improved by the comparison with the unordered data from the human genomic sequence publicly available. Phylogenetical analysis suggests that STK33 belongs to the calcium/calmodulin-dependent protein kinases group and seems to be a novelty in the chordate lineage. The gene, as a whole, seems to evolve under purifying selection whereas some regions appear to be under strong positive selection. Both human and mouse versions of serine/threonine kinase 33, consists of seventeen exons highly conserved in the coding regions, particularly in those coding for the core protein kinase domain. Also the exon/intron structure in the coding regions of the gene is conserved between human and mouse. The existence and functionality of the gene is supported by the presence of entries in the EST databases and was in vivo fully confirmed by isolating specific transcripts from human uterus total RNA and from several mouse tissues. Strong evidence for alternative splicing was found, which may result in tissue-specific starting points of transcription and in some extent, different protein N-termini. RT-PCR and hybridisation experiments suggest that STK33/Stk33 is differentially expressed in a few tissues and in relative low levels. STK33 has been shown to be reproducibly down-regulated in tumor tissues, particularly in ovarian tumors. RNA in-situ hybridisation experiments using mouse Stk33-specific probes showed expression in dividing cells from lung and germinal epithelium and possibly also in macrophages from kidney and lungs. Preliminary experimentation with antibodies designed in this work, performed in parallel to the preparation of this manuscript, seems to confirm this expression pattern. The fact that the chromosomal region 11p15 in which STK33 is located may be associated with several human diseases including tumor development, suggest further investigation is necessary to establish the role of STK33 in human health.

Epigenetic inactivation of chalcone synthase-A transgene transcription in petunia leads to a reversion of the post-transcriptional gene silencing phenotype

Petunia plants that exhibit a white-flowering phenotype as a consequence of chalcone synthase transgene-induced silencing occasionally give rise to revertant branches that produce flowers with wild-type pigmentation. Transcription run-on assays confirmed that the production of white flowers is caused by post-transcriptional gene silencing (PTGS), and indicated that transgene transcription is repressed in the revertant plants, providing evidence that induction of PTGS depends on the transcription rate. Transcriptional repression of the transgene was associated with cytosine methylation at CpG, CpNpG and CpNpN sites, and the expression was restored by treatment with either 5-azacytidine or trichostatin A. These results demonstrate that epigenetic changes occurred in the PTGS line, and these changes interfere with the initiation of transgene transcription, leading to a reversion of the PTGS phenotype.

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.

The isolation and characterization of tescalcin, a novel gene differentially expressed in the mouse testis during the early stages of gonadal differentiation

Gonadal development is an ideal model to study organogenesis because a variety of developmental processes can be studied during the differentiation of the bipotential primordium into testis or ovary. To better understand this process, Representational Difference Analysis of cDNA was used to identify genes that are differentially expressed in mouse gonads at 13.5 days post-coitus. The analysis led to the identification of three testis specific genes and a sequence that was only expressed in the ovary. The male genes identified: renin, Col9a3, and a novel gene termed tescalcin had patterns of expression that suggested a role in testis determination. ^ Studies of the tescalcin gene revealed that it is organized into eight exons and seven introns. The gene was located at 64 cM in mouse chromosome 5, where it spans approximately 35 Kb. Three mRNA variants resulting from alternative splicing of intron 5 were identified in mouse tissues. Gel mobility shift assays demonstrated that Sp1 and Sp3 from Y-1, msc-1, and MIN-6 cells nuclear extracts bind the GC-boxes within the tescalcin proximal promoter. Bisulfite sequencing analysis of tescalcin CpG island revealed that it is differentially methylated in male and female mouse embryonic gonads, and that hypermethylation of this region represses expression of tescalcin in the β-TC3 cell line. ^ The major tescalcin mRNA encodes a protein with 214 amino acids that contains a consensus EF-hand Ca2+-binding domain and an N-myristoylation motif. The amino acid sequence of tescalcin is highly conserved among various species, and it showed the highest homology with calcineurin B homologous proteins 1 and 2, and calcineurin B. Western blot analysis using antibodies generated against the tescalcin protein confirmed its presence in specific mouse tissues and cell lines. Immunohistochemical analysis of mouse embryos confirmed the pattern of expression of tescalcin mRNA in fetal testis. Using pull-down assays, glyceraidehydes-3-phosphate dehydrogenase was identified as an interacting and potential functional partner of tescalcin. ^ The identification and characterization of tescalcin as a novel embryonic testicular marker will contribute to the elucidation of the genetic pathways involved in testis development and likely to the understanding of pathological conditions such as sex reversal and infertility. ^

Cloning, sequencing and expression analysis of cDNA encoding a constitutive heat shock protein 70 (HSC70) in Fenneropenaeus chinensis

The cDNA encoding hsc70 of Chinese shrimp Fenneropenaeus chinensis was cloned from hepatopancreas by RT-PCR based on its EST sequence. The full length cDNA of 2090 bp contained an open reading frame of 1956 nucleotides and partial 5'- and 3'-untranslated region(5'- and 3'-UTR). PCR amplification and sequencing analysis showed the existence of introns in the region of 1-547 bp, but they did not exist in the region of 548-2090 bp of hsc70 cDNA. When the deduced 652 amino acid sequence of HSC70 was compared with the members of HSP70 family from other organisms, the results showed 85.9% similarity with HSC71 from Oncorhynchus mykiss and HSC70 from Homo sapiens. It also exhibited 85.8% similarity with HSP70 from Mus musculu and 85.4% with HSC70 from Manduca sexta. Expression analysis showed that hsc70 mRNA was espressed constitutively in hepatopancreas, muscle, eyestalks, haemocytes, heart, ovary, intestine and gills in Fenneropenaeus chinensis. No difference could be detected on hsc70 mRNA level in muscle between heat-shocked and control animals.

IGFBP7 promoter methylation and gene expression analysis in prostate cancer

PURPOSE: IGFBP7 belongs to a family of insulin-like growth factor-1 regulatory binding proteins. IGFBP7 hypermethylation is associated with its down-regulation in various carcinomas. In prostate cancer IGFBP7 down-regulation has been widely reported but to our knowledge the mechanisms behind this event are unknown. We performed a denaturing high performance liquid chromatography screening and validation strategy to profile the methylation status of IGFBP7 in prostate cancer.

MATERIALS AND METHODS: We combined denaturing high performance liquid chromatography and bisulfite sequencing to examine IGFBP7 methylation in a panel of prostate cancer cell lines. Quantitative methylation specific polymerase chain reaction was used to determine methylation levels in prostate tissue specimens of primary prostate cancer, histologically benign prostate adjacent to tumor, high grade prostatic intraepithelial neoplasia and benign prostatic hyperplasia. IGFBP7 gene expression was measured by quantitative methylation specific polymerase chain reaction in cell lines and tissue specimens.

RESULTS: IGFBP7 was methylated in the 4 prostate cancer cell lines DU145, LNCaP, PC-3 and 22RV1. Quantitative methylation specific polymerase chain reaction analysis revealed that promoter methylation was associated with decreased IGFBP7 expression. Quantitative methylation specific polymerase chain reaction showed that IGFBP7 methylation was more frequently detected in prostate cancer (60% (31/52)) and high grade prostatic intraepithelial neoplasia (40% (6/15)) samples compared to histologically benign prostate adjacent to tumor (10%) and benign prostatic hyperplasia (0%) samples.

CONCLUSIONS: To our knowledge this is the first report of aberrant IGFBP7 promoter hypermethylation and concurrent IGFBP7 gene silencing in prostate cancer cell lines. Results demonstrate that CpG methylation of IGFBP7 may represent a novel biomarker of prostate cancer and pre-invasive neoplasms. Thus, future examination of IGFBP7 methylation and expression in a larger patient cohort, including bodily fluids, is justified to further evaluate its role in a diagnostic and prognostic setting.

Next-generation sequencing of cervical DNA detects human papillomavirus types not detected by commercial kits

Background Human papillomavirus (HPV) is the aetiological agent for cervical cancer and genital warts. Concurrent HPV and HIV infection in the South African population is high. HIV positive (+) women are often infected with multiple, rare and undetermined HPV types. Data on HPV incidence and genotype distribution are based on commercial HPV detection kits, but these kits may not detect all HPV types in HIV + women. The objectives of this study were to (i) identify the HPV types not detected by commercial genotyping kits present in a cervical specimen from an HIV positive South African woman using next generation sequencing, and (ii) determine if these types were prevalent in a cohort of HIV-infected South African women. Methods Total DNA was isolated from 109 cervical specimens from South African HIV + women. A specimen within this cohort representing a complex multiple HPV infection, with 12 HPV genotypes detected by the Roche Linear Array HPV genotyping (LA) kit, was selected for next generation sequencing analysis. All HPV types present in this cervical specimen were identified by Illumina sequencing of the extracted DNA following rolling circle amplification. The prevalence of the HPV types identified by sequencing, but not included in the Roche LA, was then determined in the 109 HIV positive South African women by type-specific PCR. Results Illumina sequencing identified a total of 16 HPV genotypes in the selected specimen, with four genotypes (HPV-30, 74, 86 and 90) not included in the commercial kit. The prevalence's of HPV-30, 74, 86 and 90 in 109 HIV positive South African women were found to be 14.6 %, 12.8 %, 4.6 % and 8.3 % respectively. Conclusions Our results indicate that there are HPV types, with substantial prevalence, in HIV positive women not being detected in molecular epidemiology studies using commercial kits. The significance of these types in relation to cervical disease remains to be investigated.