Can the occurrence of rare insertion/deletion polymorphisms in human mtDNA be verified from phylogeny?

Due to its specific characteristics, such as maternal inheritance and absence of recombination, each mtDNA belongs to certain monophyletic clade in the rooted mtDNA tree (haplogroup) according to the mutations it harbors. Rare mutation (excluding parallel mutation) occurring at multiple times in different haplogroups could thus be a potential reading error according to the mtDNA phylogeny. This experience has been widely used in double-checking the credibility of the rare mutations in human mtDNA sequences. However, no test has been performed so far for the feasibility of applying this strategy to the rare insertion/deletion (indel) events in mtDNA sequences. In this study, we attempted to relate the rare indels in mtDNAs to their haplogroup status in a total of 2352 individuals from 50 populations in China. Our results show that the insertion of A at position 16259 is restricted to a subclade of haplogroup C and can be verified. The other indel polymorphisms, which occur in the repeat of the deleted or inserted nucleotide(s), may not be distinguished from phantom mutations from a phylogenetic point of view. Independently and multiply sequencing the fragment with the indel is the best and the most reliable way for confirmation.

Evolution of the DRD2 gene haplotype and its association with alcoholism in Mexican Americans

The human D2 dopamine receptor gene (DRD2) plays a central role in the neuromodulation of appetitive behaviors and is implicated in having a possible role in susceptibility to alcoholism. We genotyped an SNP in DRD2 Exon 8 in 251 nonalcoholic, unrelated, healthy controls and 200 alcoholic Mexican Americans. The DRD2 haplotypes were analyzed using the Exon 8 genotype in combination with five other SNP genotypes, which were obtained from our previous study. The ancestral origins of the DRD2 polymorphisms have been determined by sequencing the homologous region in other higher primates. Twenty DRD2 haplotypes, defined as H1 to H20 based on their frequency from high to low, were obtained in this major minority population. The ancestral haplotype "I-132-G-C-G-A1" and two one-step mutation haplotypes were absent in our study population. The haplotype H1, "I-B1-T-C-A-A1", with the highest frequency in the population, is a three-step mutation from the ancestral form. The first five or eight major haplotypes make up 87% or 95% of the entire population, respectively. The prevalence of the haplotype H1+ (H1/H1 and H1/Hn genotypes) is significantly higher in alcoholics and alcoholic subgroups, including early onset drinkers and benders, than in their respective control groups. The Promoter -141C allele is in linkage disequilibrium (LD) with five other loci in the nonalcoholic group, but not in the alcoholic group. All of the other five loci are in LD in both the alcoholic and control groups. The DRD2 TaqI B allele is in complete LD with the allele located in intron 6. Five SNPs, Promoter -141C, TaqI B (or Intron 6), Exon 7, Exon 8, and TaqI A, are sufficient to define the DRD2 haplotypes in Mexican Americans. Our data indicate that the DRD2 haplotypes are associated with alcoholism in Mexican Americans. (c) 2005 Elsevier Inc. All rights reserved.

Episodic evolution of prolactin gene in primates

In the present study, we obtained exon 2-5 of prolactin (PRL) gene from four primate species by PCR and sequencing. Adding other genes available in GenBank, we calculate amino acid substitution rates for prolactin gene in primate. Comparison of nonsynonymous substitution rate to synonymous substitution rate ratios shows no evidence of positive selection for any lineage of primate prolactin gene. According to this and the facts that (i) no sites under positive selection are inferred by using maximum-likelihood method; (ii) among 32 amino acid replacement that occurred along the rapid evolutionary phase, only two are included in the 40 functionally important residues, indicating that amino acid replacement tends to occur in those functionally unimportant residues; (iii) partial of prolactin function is replaced by placental lactogen in primate at the rapid evolutionary phase of prolactin gene, we thus deem that it is relaxation of purifying selection to some extent rather than positive selection that enforces the rapid evolution of primate prolactin gene.

Low “penetrance” of phylogenetic knowledge in mitochondrial disease studies

An up-to-date view of the worldwide mitochondrial DNA (mtDNA) phylogeny together with an evaluation of the conservation of each site is a reliable tool for detecting errors in mtDNA studies and assessing the functional importance of alleged pathogenic mutations. However, most of the published studies on mitochondrial diseases make very little use of the phylogenetic knowledge that is currently available. This drawback has two inadvertent consequences: first, there is no sufficient a posteriori quality assessment of complete mtDNA sequencing efforts; and second, no feedback is provided for the general mtDNA database when apparently new mtDNA lineages are discovered. We demonstrate, by way of example, these issues by reanalysing three mtDNA sequencing attempts, two from Europe and another one from East Asia. To further validate our phylogenetic deductions, we completely sequenced two mtDNAs from healthy subjects that nearly match the mtDNAs of two patients, whose sequences gave problematic results. (c) 2005 Elsevier Inc. All rights reserved.

Genome 10K: A Proposal to Obtain Whole-Genome Sequence for 10 000 Vertebrate Species

The human genome project has been recently complemented by whole-genome assessment sequence of 32 mammals and 24 nonmammalian vertebrate species suitable for comparative genomic analyses. Here we anticipate a precipitous drop in costs and increase in sequ

Polymerase chain reaction based C4AQ0 and C4BQ0 genotyping: association with systemic lupus erythematosus in southwest Han Chinese

Objective: To investigate the association of complement C4 null genes (C4QO, including C4AQO and C4BQO) and C2 gene with systemic lupus erythematosus (SLE) in southwest Han Chinese; 136 patients with SLE and 174 matched controls were genotyped. Methods: C4 null genes were determined by a polymerase chain reaction (PCR) procedure with sequence specific primers (PCR-SSP). The 2 bp insertion in exon 29, which was previously identified in non-Chinese populations and caused defective C4A genes, was directly typed by sequencing the whole exon 29 using exon specific primers. The exon 6 of complement C2 was also sequenced in both the patients and controls. Results: The frequency of homozygous C4AQO allele was 12.5% (17/136) in patients with SLE compared with 1.1% (2/174) in controls (p<0.001, odds ratio (OR)=12.286, 95% confidence interval (95% CI) 2.786 to 54.170). There was no significant difference for homozygous C4BQO allele between patients with SLE and controls (p=0.699). Patients with the C4AQO gene had an increased risk of acquiring renal disorder, serositis, and anti-dsDNA antibodies compared with those without C4AQO (for renal disorder, p=0.018, OR=8.951, 95% Cl 1.132 to 70.804; for serositis, p=0.011, OR 4.891, 95% CI 1.574 to 15.198; for anti-dsDNA, p=0.004, OR 7.630, 95%Cl 1.636 to 35.584). None of the patients or controls had the 2 bp insertion in exon 29 of the C4 gene. The type I C2 deficiency was not detected in the 3 10 samples. Conclusion: It is suggested that deficiency of C4A (not due to a 2 bp insertion in exon 29), but not C4B or C2, may be a risk factor for acquiring SLE in south west Han Chinese; this results in increased risk of renal disorder, serositis, and anti-dsDNA antibodies in patients with SLE. Racial differences seem to be relevant in susceptibility to SLE.

Molecular characterization of six Chinese families with m.3460G > A and Leber hereditary optic neuropathy

The primary mutation m.3460G > A occurs with a very low frequency (similar to 1%) in Chinese patients with Leber hereditary optic neuropathy (LHON). Up to now, there is no comprehensive study of Chinese patients harboring this mutation. We characterized six unrelated probands with m.3460G > A in this study, which were identified from 1,626 patients with LHON or suspected with LHON. The overall penetrance of LHON (25.6% [10/39]) in four pedigrees with m.3460G > A was substantially lower than those families with m.11778G > A (33.3% [619/1859]) as reported in our previous study. Intriguingly, family Le688 with a heteroplasmic m.3460G > A presented a lower penetrance (12.5%) than the other three families with a homoplasmic mutation. There is an elevated gender bias (affected male to affected female = 4:1) in the four families with m.3460G > A compared to those LHON families with m.11778G > A (2.4:1). Complete mtDNA sequencing indicated that the six matrilines belonged to haplogroups B4d1, F2, A5b, M12a, D4b2b, and D4b2, respectively. We did not identify any potential secondary mutation(s) that will affect or be associated with the penetrance of LHON in the six probands by using an evolutionary analysis and protein secondary-structure prediction. Taken together, our results suggested that the m.3460G > A mutation occurred multiple times in Chinese LHON patients. The heteroplasmic status of mutation m.3460G > A might influence the penetrance of LHON in family Le688.

The Acquisition of an Inheritable 50-bp Deletion in the Human mtDNA Control Region Does Not Affect the mtDNA Copy Number in Peripheral Blood Cells

The mitochondrial DNA (mtDNA) control region is believed to play an important biological role in mtDNA replication. Large deletions in this region are rarely found, but when they do occur they might be expected to interfere with the replication of the molecule, thus leading to a reduction of mtDNA copy number. During a survey for mtDNA sequence variations in 5,559 individuals from the general Chinese population and 2,538 individuals with medical disorders, we identified a 50-bp deletion (m.298_347del50) in the mtDNA control region in a member of a healthy Han Chinese family belonging to haplogroup B4c1b2, as suggested by complete mtDNA genome sequencing. This deletion removes the conserved sequence block II (CSBII; region 299-315) and the replication primer location (region 317-321). However, quantification of the mtDNA copy number in this subject showed a value within a range that was observed in 20 healthy subjects without the deletion. The deletion was detected in the hair samples of the maternal relatives of the subject and exhibited variable heteroplasmy. Our current observation, together with a recent report for a benign 154-bp deletion in the mtDNA control region, suggests that the control of mtDNA replication may be more complex than we had thought. Hum Mutat 31:538-543, 2010. (C) 2010 Wiley-Liss, Inc.

Phantom mutation hotspots in human mitochondrial DNA

Phantom mutations are systematic artifacts generated in the course of the sequencing process. Contra common belief these artificial mutations are nearly ubiquitous in sequencing results, albeit at frequencies that may vary dramatically. The amount of artifacts depends not only on the sort of automated sequencer and sequencing chemistry employed, but also on other lab-specific factors. An experimental study executed on four samples under various combinations of sequencing conditions revealed a number of phantom mutations occurring at the same sites of mitochondrial DNA (mtDNA) repeatedly. To confirm these and identify further hotspots for artifacts, > 5000 mtDNA electropherograms were screened for artificial patterns. Further, > 30000 published hypervariable segment 1 sequences were compared at potential hotspots for phantom mutations, especially for variation at positions 16085 and 16197. Resequencing of several samples confirmed the artificial nature of these and other polymorphisms in the original publications. Single-strand sequencing, as typically executed in medical and anthropological studies, is thus highly vulnerable to this kind of artifacts. In particular, phantom mutation hotspots could easily lead to misidentification of somatic mutations and to misinterpretations in all kinds of clinical mtDNA studies.

Pseudomitochondrial genome haunts disease studies

The accidental amplification of nuclear mitochondrial pseudogenes (NUMTs) can pose a serious problem for mitochondrial disease studies. This report shows that the mutation spectrum left by spurious amplification of a NUMT can be detected because it usuall

Single base-resolution methylome of the silkworm reveals a sparse epigenomic map

Epigenetic regulation in insects may have effects on diverse biological processes. Here we survey the methylome of a model insect, the silkworm Bombyx mori, at single-base resolution using Illumina high-throughput bisulfite sequencing (MethylC-Seq). We conservatively estimate that 0.11% of genomic cytosines are methylcytosines, all of which probably occur in CG dinucleotides. CG methylation is substantially enriched in gene bodies and is positively correlated with gene expression levels, suggesting it has a positive role in gene transcription. We find that transposable elements, promoters and ribosomal DNAs are hypomethylated, but in contrast, genomic loci matching small RNAs in gene bodies are densely methylated. This work contributes to our understanding of epigenetics in insects, and in contrast to previous studies of the highly methylated genomes of Arabidopsis(1) and human(2), demonstrates a strategy for sequencing the epigenomes of organisms such as insects that have low levels of methylation.

The sequence and de novo assembly of the giant panda genome

Using next-generation sequencing technology alone, we have successfully generated and assembled a draft sequence of the giant panda genome. The assembled contigs (2.25 gigabases (Gb)) cover approximately 94% of the whole genome, and the remaining gaps (0.05 Gb) seem to contain carnivore-specific repeats and tandem repeats. Comparisons with the dog and human showed that the panda genome has a lower divergence rate. The assessment of panda genes potentially underlying some of its unique traits indicated that its bamboo diet might be more dependent on its gut microbiome than its own genetic composition. We also identified more than 2.7 million heterozygous single nucleotide polymorphisms in the diploid genome. Our data and analyses provide a foundation for promoting mammalian genetic research, and demonstrate the feasibility for using next-generation sequencing technologies for accurate, cost-effective and rapid de novo assembly of large eukaryotic genomes.

Large Gene Family Expansions and Adaptive Evolution for Odorant and Gustatory Receptors in the Pea Aphid, Acyrthosiphon pisum

Gaining insight into the mechanisms of chemoreception in aphids is of primary importance for both integrative studies on the evolution of host plant specialization and applied research in pest control management because aphids rely on their sense of smell

蓝氏贾第虫分离株tim 基因序列分析比较研究

　为了探讨来源于不同地理位置贾第虫分离株之间的遗传学关系, 采用tim 基因扩增、限制性酶切和 DNA 序列分析方法对来自我国(C1、C2、、CH2、CH3)、柬埔寨(CAM )、澳大利亚(A 1、A 2) 和美国(BP、CDC) 共9 株 蓝氏贾第虫的基因型进行了分析比较。结果表明,A 1、A 2 和CAM 属第1 型(WB) ; CH2 和CH3 属第2 型(JH) ; C1、 C2、BP 和CDC 属第3 型(GS)。本实验结果提示, 虫株间遗传学关系并非由地理位置的远近所决定。来源于同一地 理位置的虫株其遗传学特性可有很大的差异。相反, 来源于地理位置相隔很远的虫株其遗传学特性却可极为相近。 贾第虫分离株基因型的确定可为本虫分子系统进化和分子流行病学研究提供重要资料。

DNA sequence analysis of the triose phosphate isomerase gene from isolates of Giardia lamblia

Objective To confirm the genetic relation between Giardia lamblia (G. lamblia) isolates from different geographic regions of China and other countries. Methods Genomic DNA were extracted from the trophozoites or cysts of Giardia lamblia. The triose phosphate isomerase (tim) gene was amplified using polymerase chain reaction (PCR) technique. PCR products were digested with endonuclease and sequenced. The data of sequencing were analyzed with the DNAstar software and compared with that of the isolates acquired from GenBank. Results Of nine isolates of Giardia lamblia from China (C1, C2, CH2 and CH3), Cambodia (CAM), Australia (A1 and A2) and America (BP and CDC), respectively, 3 (A1, A2 and CAM) fit into Group 1 (WB), 2 (CH2 and CH3) into Group 2, and 4 (C1, C2, BP and CDC) into Group 3 (GS). The results confirmed the genetic relatedness of G. lamblia isolates from all over the world. Conclusion Genotyping isolates of G. Lamblia provides important information for establishing the phylogenetic relationship or for the epidemiological evaluation of the spreading of this organism.