Analysis of mitochondrial length heteroplasmy in monozygous and non-monozygous siblings

Development of a multiplex assay for detection of autosomal and Y-chromosomal STRs, assessment of the degradation state of mitochondrial DNA and presence of mitochondrial length heteroplasmies

The current focus in most routine forensic casework is detection of autosomal or gonosomal Short Tandem Repeats (STRs). With increasing degradation, STR analysis tends to be less successful up to complete failure. For challenging samples such as telogen hair roots and shafts, touch DNA samples or skeletal remains, mitochondrial DNA (mtDNA) analysis provides a powerful tool. Determination of DNA quantity is an important part in the casework workflow. Several ready-to-use kits are commercially available for nuclear DNA targets. However, quantification of mtDNA targets requires the establishment of an in-house method. Some assays even contain assessment of degradation, which alleviates the choice of target enrichment for sequencing through medium or small amplicons. As Sanger-type Sequencing (STS) still remains the golden standard in many laboratories, identification of heteroplasmies in C-tract regions prior to the sequencing reaction is advantageous. Firstly, primer selection can be expanded with primers binding near the C-tract and secondly, determination of the dominant variant is straightforward. All those quantity (nuclear and mtDNA) and quality (degradation and length heteroplasmies) evaluations usually require at least two separate reactions. Therefore, the aim of this project was the combination of all these targets in one multiplex assay using capillary electrophoresis to spare valuable sample extract. Amplification of representative autosomal and Y-chromosomal STRs allows estimate of success of (Y-)STR analysis. Simultaneously, five length heteroplasmies in the mitochondrial control region are targeted as well as three conservative regions of differing fragment lengths for assessment of the mitochondrial degradation state. Based on the outcome of this assay, forensic examiners can decide if STR analysis may be suitable. In case of absent STR peaks, appropriate proceeding of mtDNA sequencing can be determined.

Multiplexed Microsphere Suspension-Array Assay for Urine Mitochondrial DNA Typing by C-Stretch Length in Hypervariable Regions

Background

The standard method for personal identification and verification of urine samples in doping control is short tandem repeat (STR) analysis using nuclear DNA (nDNA). The DNA concentration of urine is very low and decreases under most conditions used for sample storage; therefore, the amount of DNA from cryopreserved urine samples may be insufficient for STR analysis. We aimed to establish a multiplexed assay for urine mitochondrial DNA typing containing only trace amounts of DNA, particularly for Japanese populations.

Methods

A multiplexed suspension-array assay using oligo-tagged microspheres (Luminex MagPlex-TAG) was developed to measure C-stretch length in hypervariable region 1 (HV1) and 2 (HV2), five single nucleotide polymorphisms (SNPs), and one polymorphic indel. Based on these SNPs and the indel, the Japanese population can be classified into five major haplogroups (D4, B, M7a, A, D5). The assay was applied to DNA samples from urine cryopreserved for 1 - 1.5 years (n = 63) and fresh blood (n = 150).

Results

The assay with blood DNA enabled Japanese subjects to be categorized into 62 types, exhibiting a discriminatory power of 0.960. The detection limit for cryopreserved urine was 0.005 ng of nDNA. Profiling of blood and urine pairs revealed that 5 of 63 pairs showed different C-stretch patterns in HV1 or HV2.

Conclusions

The assay described here yields valuable information in terms of the verification of urine sample sources employing only trace amounts of recovered DNA. However, blood cannot be used as a reference sample.

Length heterogeneity at conserved sequence block 2 in human mitochondrial DNA acts as a rheostat for RNA polymerase POLRMT activity

The guanine (G)-tract of conserved sequence block 2 (CSB 2) in human mitochondrial DNA can result in transcription termination due to formation of a hybrid G-quadruplex between the nascent RNA and the nontemplate DNA strand. This structure can then influence genome replication, stability and localization. Here we surveyed the frequency of variation in sequence identity and length at CSB 2 amongst human mitochondrial genomes and used in vitro transcription to assess the effects of this length heterogeneity on the activity of the mitochondrial RNA polymerase, POLRMT. In general, increased G-tract length correlated with increased termination levels. However, variation in the population favoured CSB 2 sequences which produced efficient termination while particularly weak or strong signals were avoided. For all variants examined, the 3′ end of the transcripts mapped to the same downstream sequences and were prevented from terminating by addition of the transcription factor TEFM. We propose that CSB 2 length heterogeneity allows variation in the efficiency of transcription termination without affecting the position of the products or the capacity for regulation by TEFM.

Mitochondrial DNA Variation and Heteroplasmy in Monozygotic Twins Clinically Discordant for Multiple Sclerosis

We examined the debated link between mitochondrial DNA (mtDNA) variation and multiple sclerosis (MS) using 49 monozygotic (MZ) twin pairs clinically discordant for MS, which enables to associate de novo mtDNA variants, skewed heteroplasmy, and mtDNA copy number with MS manifestation. Ultra-deep sequencing of blood-derived mtDNA revealed 25 heteroplasmic variants with potentially pathogenic features in 18 pairs. All variants were pair-specific and had low and/or similar heteroplasmy levels in both cotwins. In one pair, a confirmed pathogenic variant (m.11778G>A, heteroplasmy ∼50%) associated with Leber hereditary optic neuropathy was detected. Detailed diagnostic investigation revealed subclinical MS signs in the prior nondiseased cotwin. Moreover, neither mtDNA deletions nor copy-number variations were involved. Furthermore, the majority of heteroplasmic variants were shared among MZ twins and exhibited more similar heteroplasmy levels in the same tissue of MZ twins as compared with different tissues of the same individual. Heteroplasmy levels were also more similar within MZ twins compared with nonidentical siblings. Our analysis excludes mtDNA variation as a major driver of the discordant clinical manifestation of MS in MZ twins, and provides valuable insights into the occurrence and distribution of heteroplasmic variants within MZ twins and nonidentical siblings, and across different tissues.

Length heteroplasmy of the polyC-polyT-polyC stretch in the dog mtDNA control region

Previously, the mitochondrial control region of 214 Belgian dogs was sequenced. Analysis of this data indicated length heteroplasmy of the polyT stretch in the polyC-polyT-polyC stretch from positions 16661 to 16674. Nine polyC-polyT-polyC haplotype combinations were observed, consisting of seven major haplotypes (highest signal intensity) combined with minor haplotypes (lower signal intensity) one T shorter than the major haplotype in all but three dogs. The longer the polyT stretch, the smaller was the difference in signal intensity between the major and minor haplotype peaks. Additional sequencing, cloning, and PCR trap experiments were performed to further study the intra-individual variation of this mitochondrial DNA (mtDNA) region. Cloning experiments demonstrated that the proportion of clones displaying the minor haplotypes also increased with the length of the polyT stretch. Clone amplification showed that in vitro polymerase errors might contribute to the length heteroplasmy of polyT stretches with at least 10 Ts. Although major and minor polyC-polyT-polyC haplotypes did not differ intra-individually within and between tissues in this study, interpretation of polyT stretch variation should be handled with care in forensic casework.

Mitochondrial DNA heteroplasmy in diabetes and normal adults: role of acquired and inherited mutational patterns in twins

Heteroplasmy, the mixture of mitochondrial genomes (mtDNA), varies among individuals and cells. Heteroplasmy levels alter the penetrance of pathological mtDNA mutations, and the susceptibility to age-related diseases such as Parkinson's disease. Although mitochondrial dysfunction occurs in age-related type 2 diabetes mellitus (T2DM), the involvement of heteroplasmy in diabetes is unclear. We hypothesized that the heteroplasmic mutational (HM) pattern may change in T2DM. To test this, we used next-generation sequencing, i.e. massive parallel sequencing (MPS), along with PCR-cloning-Sanger sequencing to analyze HM in blood and skeletal muscle DNA samples from monozygotic (MZ) twins either concordant or discordant for T2DM. Great variability was identified in the repertoires and amounts of HMs among individuals, with a tendency towards more mutations in skeletal muscle than in blood. Whereas many HMs were unique, many were either shared among twin pairs or among tissues of the same individual, regardless of their prevalence. This suggested a heritable influence on even low abundance HMs. We found no clear differences between T2DM and controls. However, we found ~5-fold increase of HMs in non-coding sequences implying the influence of negative selection (P < 0.001). This negative selection was evident both in moderate to highly abundant heteroplasmy (>5% of the molecules per sample) and in low abundance heteroplasmy (<5% of the molecules). Although our study found no evidence supporting the involvement of HMs in the etiology of T2DM, the twin study found clear evidence of a heritable influence on the accumulation of HMs as well as the signatures of selection in heteroplasmic mutations.

Evaluating sequence-derived mtDNA length heteroplasmy by amplicon size analysis

Length heteroplasmy (LH) in mitochondrial (mt)DNA is usually observed in homopolymeric tracts and manifest as mixture of various length variants. The generally used difference-coded annotation to report mtDNA haplotypes does not express the degree of LH variation present in a sample, even more so, it is sometimes difficult to establish which length variants are present and clearly distinguishable from background noise. It has therefore become routine practice for some researchers to call the dominant type, the "major molecule", which represents the LH variant that is most abundant in a DNA extract. In the majority of cases a clear single dominant variant can be identified. However, in some samples this interpretation is difficult, i.e. when (almost) equally quantitative LH variants are present or when multiple sequencing primers result in the presentation of different dominant types. To better understand those cases we designed amplicon sizing assays for the five most relevant LH regions in the mtDNA control region (around ntps 16,189, 310, 460, 573, and the AC-repeat between 514 and 524) to determine the ratio of the LH variants by fluorescence based amplicon sizing assays. For difficult LH constellations derived by Sanger sequencing (with Big Dye terminators) these assays mostly gave clear and unambiguous results. In the vast majority of cases we found agreement between the results of the sequence and amplicon analyses and propose this alternative method in difficult cases.

Behavioral genetic approach to the study of dyslexia

OBJECTIVE

Dyslexia is a prominent focus of practitioners, educators, and researchers because of the myriad consequences of failing to read proficiently. The aim of this study was to provide a brief overview of how twin studies can offer insight on the cause of many human behaviors and disorders including dyslexia, discuss common misconceptions regarding findings from behavioral genetic studies, briefly review the evidence on the relationship between genes, environment, and dyslexia, and finally present some findings from a large-scale twin study on reading and dyslexia.

METHOD

Participants were twins from a large ethnically and socioeconomically diverse twin sample in an ongoing longitudinal study of reading and dyslexia. Heritabilities of reading ability and dyslexia were calculated for 1,024 first grade twins on a standardized reading measure. Children were identified as dyslexic if they scored at the 15th percentile or below on a reading measure.

RESULTS

Relatively high heritabilities were observed for both reading ability and dyslexia indicating substantial genetic influences. Further, results indicated some overlap of genetic factors influencing reading ability and dyslexia.

CONCLUSIONS

Behavioral genetic studies offer a means of understanding the cause of dyslexia. This study extended research to a more diverse sample than extant studies and found lower heritability estimates of reading ability and dyslexia, but a similar pattern of results indicating possible genetic overlap. Twin studies provide perspective for discoveries of specific genes involved in dyslexia by quantifying the amount of variance waiting to be accounted for by genes while simultaneously providing an impetus to continue working on efforts for environmental intervention.

Sequence-length variation of mtDNA HVS-I C-stretch in Chinese ethnic groups

The purpose of this study was to investigate mitochondrial DNA (mtDNA) hypervariable segment-I (HVS-I) C-stretch variations and explore the significance of these variations in forensic and population genetics studies. The C-stretch sequence variation was studied in 919 unrelated individuals from 8 Chinese ethnic groups using both direct and clone sequencing approaches. Thirty eight C-stretch haplotypes were identified, and some novel and population specific haplotypes were also detected. The C-stretch genetic diversity (GD) values were relatively high, and probability (P) values were low. Additionally, C-stretch length heteroplasmy was observed in approximately 9% of individuals studied. There was a significant correlation (r=-0.961, P<0.01) between the expansion of the cytosine sequence length in the C-stretch of HVS-I and a reduction in the number of upstream adenines. These results indicate that the C-stretch could be a useful genetic maker in forensic identification of Chinese populations. The results from the Fst and dA genetic distance matrix, neighbor-joining tree, and principal component map also suggest that C-stretch could be used as a reliable genetic marker in population genetics.