Intronic non-CG DNA hydroxymethylation and alternative mRNA splicing in honey bees.

BACKGROUND: Previous whole-genome shotgun bisulfite sequencing experiments showed that DNA cytosine methylation in the honey bee (Apis mellifera) is almost exclusively at CG dinucleotides in exons. However, the most commonly used method, bisulfite sequencing, cannot distinguish 5-methylcytosine from 5-hydroxymethylcytosine, an oxidized form of 5-methylcytosine that is catalyzed by the TET family of dioxygenases. Furthermore, some analysis software programs under-represent non-CG DNA methylation and hydryoxymethylation for a variety of reasons. Therefore, we used an unbiased analysis of bisulfite sequencing data combined with molecular and bioinformatics approaches to distinguish 5-methylcytosine from 5-hydroxymethylcytosine. By doing this, we have performed the first whole genome analyses of DNA modifications at non-CG sites in honey bees and correlated the effects of these DNA modifications on gene expression and alternative mRNA splicing. RESULTS: We confirmed, using unbiased analyses of whole-genome shotgun bisulfite sequencing (BS-seq) data, with both new data and published data, the previous finding that CG DNA methylation is enriched in exons in honey bees. However, we also found evidence that cytosine methylation and hydroxymethylation at non-CG sites is enriched in introns. Using antibodies against 5-hydroxmethylcytosine, we confirmed that DNA hydroxymethylation at non-CG sites is enriched in introns. Additionally, using a new technique, Pvu-seq (which employs the enzyme PvuRts1l to digest DNA at 5-hydroxymethylcytosine sites followed by next-generation DNA sequencing), we further confirmed that hydroxymethylation is enriched in introns at non-CG sites. CONCLUSIONS: Cytosine hydroxymethylation at non-CG sites might have more functional significance than previously appreciated, and in honey bees these modifications might be related to the regulation of alternative mRNA splicing by defining the locations of the introns.

Spectrum of mutations in mut methylmalonic acidemia and identification of a common Hispanic mutation and haplotype.

Cobalamin nonresponsive methylmalonic acidemia (MMA, mut complementation class) results from mutations in the nuclear gene MUT, which codes for the mitochondrial enzyme methylmalonyl CoA mutase (MCM). To better elucidate the spectrum of mutations that cause MMA, the MUT gene was sequenced in 160 patients with mut MMA. Sequence analysis identified mutations in 96% of disease alleles. Mutations were found in all coding exons, but predominantly in exons 2, 3, 6, and 11. A total of 116 different mutations, 68 of which were novel, were identified. Of the 116 different mutations, 53% were missense mutations, 22% were deletions, duplications or insertions, 16% were nonsense mutations, and 9% were splice-site mutations. Sixty-one of the mutations have only been identified in one family. A novel mutation in exon 2, c.322C>T (p.R108C), was identified in 16 of 27 Hispanic patients. SNP genotyping data demonstrated that Hispanic patients with this mutation share a common haplotype. Three other mutations were seen exclusively in Hispanic patients: c.280G>A (p.G94R), c.1022dupA, and c.970G>A (p.A324T). Seven mutations were seen almost exclusively in black patients, including the previously reported c.2150G>T (p.G717V) mutation, which was identified in 12 of 29 black patients. Two mutations were seen only in Asian patients. Some frequently identified mutations were not population-specific and were identified in patients of various ethnic backgrounds. Some of these mutations were found in mutation clusters in exons 2, 3, 6, and 11, suggesting a recurrent mutation.

The c.419-420insA in the MTP gene is associated with abetalipoproteinemia among French-Canadians.

Abetalipoproteinemia (ABL) is a rare autosomal recessive disease characterised by the absence of apolipoprotein B (apoB) containing lipoproteins and, in consequence, very low triglyceride and cholesterol levels. Microsomal triglyceride transfer protein (MTP) has been associated with ABL. A search for sequence variants in the large subunit of MTP in a kindred of 10 individuals from Saguenay-Lac-St Jean area with a propositus exhibiting ABL as well as in four independent patients from the greater Quebec city area and exhibiting very low apoB and LDL-cholesterol levels identified 12 variations. Only one sequence variation, the c.419-420insA, was observed, in the homozygous form, in the abetalipoproteinemic patient. The -493G/-400A/-164T/282G/383T/419-420insA/453T/891C/969T/1151A/2884G haplotype carries the insertion and was found in all members of the family studied. In conclusion, the present study showed that the c.419-420insA alone, in the homozygous form, is a cause of classical recessive inherited ABL in the French-Canadian population.

A survey of genetic and epigenetic variation affecting human gene expression.

The identification of human sequence polymorphisms that regulate gene expression is key to understanding human genetic diseases. We report a survey of human genes that demonstrate allelic differences in gene expression, reflecting the presence of putative allele-specific cis-acting factors of either genetic or epigenetic nature. The expression of allelic transcripts in heterozygous samples is assessed directly by relative quantitation of intragenic marker alleles in messenger or heteronuclear RNA derived from cells or tissues. This survey used 193 single-nucleotide polymorphisms (SNPs) from 129 genes expressed in lymphoblastoid cell lines, to identify 23 genes (18%) with common allele-specific transcripts whose expression deviated from the expected equimolar ratio. A subset of these deviations, or "allelic imbalances," can be observed in multiple samples derived from reference CEPH ("Centre d'Etude du Polymorphisme Humain") pedigrees and demonstrate a spectrum of patterns of transmission, including cosegregation of allelic skewing across generations compatible with Mendelian inheritance as well as random monoallelic expression for three genes (IL1A, HTR2A, and FGB). Additional studies for BTN3A2 provide evidence of SNPs and haplotypes in complete linkage disequilibrium with high- and low-expressing transcripts. The pipeline described herein offers tools for efficient identification and characterization of allelic expression allowing identification of regulatory sequence variants as well as epigenetic variation affecting human gene expression.