Mitochondrial DNA, segments I and II, its applications in population genetics and forensic sciences

The investigations of human mitochondrial DNA (mtDNA) have considerably contributed to human evolution and migration. The Middle East is considered to be an essential geographic area for human migrations out of Africa since it is located at the crossroads of Africa, and the rest of the world. United Arab Emirates (UAE) population inhabits the eastern part of Arabian Peninsula and was investigated in this study. Published data of 18 populations were included in the statistical analysis. The diversity indices showed (1) high genetic distance among African populations and (2) high genetic distance between African populations and non-African populations. Asian populations clustered together in the NJ tree between the African and European populations. MtDNA haplotypes database of the UAE population was generated. By incorporating UAE mtDNA dataset into the existing worldwide mtDNA database, UAE Forensic Laboratories will be able to analyze future mtDNA evidence in a more significant and consistent manner. ^

The relationship among serum levels of manganese superoxide dismutase and mitochondrial DNA 8-hydroxy-2'- deoxyguanosine, and dietary antioxidants intake in type 2 diabetes

Oxidative stress plays a key role in the development of Type 2 Diabetes (T2D). This cross-sectional study examined the relationship among serum levels of manganese superoxide dismutase (MnSOD), 8-hydroxy-2’-deoxyguanosine (8OHdG), dietary antioxidant intakes and glycemic control in African Americans (n=209) and Haitian Americans (n=234) with and without T2D. ^ African Americans had higher BMI (32.8 vs. 29.3 kg/m2), higher energy intake (2148 vs. 1770 kcal), and were more educated as compared to Haitian Americans; all variables were significant at p < .001. Serum levels of 8OHdG and MnSOD for African Americans (1691.0 ± 225.1 pg/ml, 2538.0 ± 1091.8 pg/ml; respectively) were significantly higher than for Haitian Americans (1626.2 ± 222.9, 2015.8 ± 656.3 pg/ml; respectively). 8OHdG was negatively correlated with MnSOD ( r = -.167, p < .001) in T2D. Having T2D was negatively correlated with MnSOD (r = -.337; p < .01) and positively correlated with 8OHdG (r = .500; p < .01). African Americans and Haitian Americans with T2D had fasting plasma glucose (FPG) levels of 143.0 ± 61.0 mg/dl and 157.6 ± 65.5 mg/dl, and A1C of 7.5 ± 1.8 % and 8.4 ± 2.4 %, respectively. African Americans and Haitian Americans without T2D had FPG levels of 95.8 ± 13.2 mg/dl and 98.7 ± 16.9 mg/dl, and A1C of 5.9 ± 0.4% and 6.0 ± 0.5%, respectively. Dietary intakes of vitamin C and vitamin D were negatively correlated with FPG (r = -.21; r = -.19, p < .05) respectively. Carotenoids negatively correlated with A1C (r = -.19, p < .05). Lower levels of MnSOD were associated with lower levels of zinc, r = .10, p < .05, and higher levels of carotenoids r = -.10, p < .05. Higher levels of 8OHdG were associated with lower levels of Vitamin D, r = -.14, p < .01, and carotenoids, r = -.09, p < .05. ^ The results demonstrate greater oxidative mtDNA damage in persons with T2D compared to those without T2D and in African Americans compared with Haitian Americans. The inverse relationship between dietary intake of antioxidants and oxidative stress implies a potential to reduce oxidative stress with diet. ^

Mitochondrial DNA damage induced autophagy, cell death, and disease.

Mammalian mitochondria contain multiple small genomes. While these organelles have efficient base excision removal of oxidative DNA lesions and alkylation damage, many DNA repair systems that work on nuclear DNA damage are not active in mitochondria. What is the fate of DNA damage in the mitochondria that cannot be repaired or that overwhelms the repair system? Some forms of mitochondrial DNA damage can apparently trigger mitochondrial DNA destruction, either via direct degradation or through specific forms of autophagy, such as mitophagy. However, accumulation of certain types of mitochondrial damage, in the absence of DNA ligase III (Lig3) or exonuclease G (EXOG), can directly trigger cell death. This review examines the cellular effects of persistent damage to mitochondrial genomes and discusses the very different cell fates that occur in response to different kinds of damage.

Conservation genetics of North Atlantic loggerhead sea turtles: analysis of nuclear and mitochondrial DNA

[EN] Complex population structure has been described for the loggerhead sea turtle (Caretta caretta), revealing lower levels of population genetic structure in nuclear compared to mitochondrial DNA assays. This may result from mating during spatially overlapping breeding migrations, or male-biased dispersal as previously found for the green turtle (Chelonia mydas). To further investigate these multiple possibilities, we carried out a comparative analysis from twelve newly developed microsatellite loci and the mitochondrial DNA control region (~804 bp) in adult females of the Cape Verde Islands (n=158), and Georgia, USA (n=17).

Maternal Copy Number Variants (CNV) transmission to their Autism Spectrum Disorder (ASD) sons correlates with phenotypic traits

Background: Common neurodevelopmental disorder, global prevalence ~1 %; Persistent deficits in social communication and social interaction; restricted and repetitive behavior, interests, or activities; Highly heterogeneous clinical presentation; Male to female ratio ~4:1.

Assessing the mitochondrial DNA diversity of the Chagas disease vector Triatoma sordida (Hemiptera: Reduviidae)

Triatoma sordida is a species that transmits Trypanosoma cruzi to humans. In Brazil, T. sordida currently deserves special attention because of its wide distribution, tendency to invade domestic environments and vectorial competence. For the planning and execution of control protocols to be effective against Triatominae, they must consider its population structure. In this context, this study aimed to characterise the genetic variability of T. sordida populations collected in areas with persistent infestations from Minas Gerais, Brazil. Levels of genetic variation and population structure were determined in peridomestic T. sordida by sequencing a polymorphic region of the mitochondrial cytochrome b gene. Low nucleotide and haplotype diversity were observed for all 14 sampled areas; π values ranged from 0.002-0.006. Most obtained haplotypes occurred at low frequencies, and some were exclusive to only one of the studied populations. Interpopulation genetic diversity analysis revealed strong genetic structuring. Furthermore, the genetic variability of Brazilian populations is small compared to that of Argentinean and Bolivian specimens. The possible factors related to the reduced genetic variability and strong genetic structuring obtained for studied populations are discussed in this paper.

Genome-wide copy number variation (CNV) detection in Nelore cattle reveals highly frequent variants in genome regions harboring QTLs affecting production traits.

Background: Copy number variations (CNVs) have been shown to account for substantial portions of observed genomic variation and have been associated with qualitative and quantitative traits and the onset of disease in a number of species. Information from high-resolution studies to detect, characterize and estimate population-specific variant frequencies will facilitate the incorporation of CNVs in genomic studies to identify genes affecting traits of importance. Results: Genome-wide CNVs were detected in high-density single nucleotide polymorphism (SNP) genotyping data from 1,717 Nelore (Bos indicus) cattle, and in NGS data from eight key ancestral bulls. A total of 68,007 and 12,786 distinct CNVs were observed, respectively. Cross-comparisons of results obtained for the eight resequenced animals revealed that 92 % of the CNVs were observed in both datasets, while 62 % of all detected CNVs were observed to overlap with previously validated cattle copy number variant regions (CNVRs). Observed CNVs were used for obtaining breed-specific CNV frequencies and identification of CNVRs, which were subsequently used for gene annotation. A total of 688 of the detected CNVRs were observed to overlap with 286 non-redundant QTLs associated with important production traits in cattle. All of 34 CNVs previously reported to be associated with milk production traits in Holsteins were also observed in Nelore cattle. Comparisons of estimated frequencies of these CNVs in the two breeds revealed 14, 13, 6 and 14 regions in high (>20 %), low (<20 %) and divergent (NEL > HOL, NEL < HOL) frequencies, respectively. Conclusions: Obtained results significantly enriched the bovine CNV map and enabled the identification of variants that are potentially associated with traits under selection in Nelore cattle, particularly in genome regions harboring QTLs affecting production traits.

Small-plasmid-mediated antibiotic resistance is enhanced by increases in plasmid copy number and bacterial fitness

Plasmids play a key role in the horizontal spread of antibiotic resistance determinants among bacterial pathogens. When an antibiotic resistance plasmid arrives in a new bacterial host, it produces a fitness cost, causing a competitive disadvantage for the plasmid-bearing bacterium in the absence of antibiotics. On the other hand, in the presence of antibiotics, the plasmid promotes the survival of the clone. The adaptations experienced by plasmid and bacterium in the presence of antibiotics during the first generations of coexistence will be crucial for the progress of the infection and the maintenance of plasmid-mediated resistance once the treatment is over. Here we developed a model system using the human pathogen Haemophilus influenzae carrying the small plasmid pB1000 conferring resistance to β-lactam antibiotics to investigate host and plasmid adaptations in the course of a simulated ampicillin therapy. Our results proved that plasmid-bearing clones compensated for the fitness disadvantage during the first 100 generations of plasmid-host adaptation. In addition, ampicillin treatment was associated with an increase in pB1000 copy number. The augmentation in both bacterial fitness and plasmid copy number gave rise to H. influenzae populations with higher ampicillin resistance levels. In conclusion, we show here that the modulations in bacterial fitness and plasmid copy number help a plasmid-bearing bacterium to adapt during antibiotic therapy, promoting both the survival of the host and the spread of the plasmid.

Species delimitation in sharpnose sharks (genus Rhizoprionodon) in the western Atlantic Ocean using mitochondrial DNA

Despite Springer’s (1964) revision of the sharpnose sharks (genus Rhizoprionodon), the taxonomic definition and ranges of Rhizoprionodon in the western Atlantic Ocean remains problematic. In particular, the distinction between Rhizoprionodon terraenovae and R. porosus, and the occurrence of R. terraenovae in South American waters are unresolved issues involving common and ecologically important species in need of fishery management in Caribbean and southwest Atlantic waters. In recent years, molecular markers have been used as efficient tools for the detection of cryptic species and to address controversial taxonomic issues. In this study 415 samples of the genus Rhizoprionodon captured in the western Atlantic Ocean from Florida to southern Brazil were examined for sequences of the COI gene and the D-loop and evaluated for nucleotide differences. The results on nucleotide composition, AMOVA tests, and relationship distances using Bayesian-likelihood method and haplotypes network, corroborates Springer’s (1964) morphometric and meristic finding and provide strong evidence that supports consideration of R. terraenovae and R. porosus as distinct species.

TiMMing: developing an innovative suite of bioinformatic tools to harmonize and track the origin of copy number alterations in the evolutive history of multiple myeloma

Multiple Myeloma (MM) is a hematologic cancer with heterogeneous and complex genomic landscape, where Copy Number Alterations (CNAs) play a key role in the disease's pathogenesis and prognosis. It is of biological and clinical interest to study the temporal occurrence of early alterations, as they play a disease "driver" function by deregulating key tumor pathways. This study presents an innovative bioinformatic tools suite created for harmonizing and tracing the origin of CNAs throughout the evolutionary history of MM. To this aim, large cohorts of newly-diagnosed MM (NDMM, N=1582) and Smoldering-MM (SMM, N=282) were aggregated. The tools developed in this study enable the harmonization of CNAs as obtained from different genomic platforms in such a way that a high statistical power can be obtained. By doing so, the high numerosity of those cohorts was harnessed for the identification of novel genes characterized as "driver" (NFKB2, NOTCH2, MAX, EVI5 and MYC-ME2-enhancer), and the generation of an innovative timing model, implemented with a statistical method to introduce confidence intervals in the CNAs-calls. By applying this model on both NDMM and SMM cohorts, it was possible to identify specific CNAs (1q(CKS1B)amp, 13q(RB1)del, 11q(CCND1)amp and 14q(MAX)del) and categorize them as "early"/ "driver" events. A high level of precision was guaranteed by the narrow confidence intervals in the timing estimates. These CNAs were proposed as critical MM alterations, which play a foundational role in the evolutionary history of both SMM and NDMM. Finally, a multivariate survival model was able to identify the independent genomic alterations with the greatest effect on patients’ survival, including RB1-del, CKS1B-amp, MYC-amp, NOTCH2-amp and TRAF3-del/mut. In conclusion, the alterations that were identified as both "early-drivers” and correlated with patients’ survival were proposed as biomarkers that, if included in wider survival models, could provide a better disease stratification and an improved prognosis definition.

Identification of outcome-oriented cut-offs for copy number alterations in multiple myeloma: predictive biomarkers with improved prognostic accuracy

Aim of the present study was to develop a statistical approach to define the best cut-off Copy number alterations (CNAs) calling from genomic data provided by high throughput experiments, able to predict a specific clinical end-point (early relapse, 18 months) in the context of Multiple Myeloma (MM). 743 newly diagnosed MM patients with SNPs array-derived genomic and clinical data were included in the study. CNAs were called both by a conventional (classic, CL) and an outcome-oriented (OO) method, and Progression Free Survival (PFS) hazard ratios of CNAs called by the two approaches were compared. The OO approach successfully identified patients at higher risk of relapse and the univariate survival analysis showed stronger prognostic effects for OO-defined high-risk alterations, as compared to that defined by CL approach, statistically significant for 12 CNAs. Overall, 155/743 patients relapsed within 18 months from the therapy start. A small number of OO-defined CNAs were significantly recurrent in early-relapsed patients (ER-CNAs) - amp1q, amp2p, del2p, del12p, del17p, del19p -. Two groups of patients were identified either carrying or not ≥1 ER-CNAs (249 vs. 494, respectively), the first one with significantly shorter PFS and overall survivals (OS) (PFS HR 2.15, p<0001; OS HR 2.37, p<0.0001). The risk of relapse defined by the presence of ≥1 ER-CNAs was independent from those conferred both by R-IIS 3 (HR=1.51; p=0.01) and by low quality (< stable disease) clinical response (HR=2.59 p=0.004). Notably, the type of induction therapy was not descriptive, suggesting that ER is strongly related to patients’ baseline genomic architecture. In conclusion, the OO- approach employed allowed to define CNAs-specific dynamic clonality cut-offs, improving the CNAs calls’ accuracy to identify MM patients with the highest probability to ER. As being outcome-dependent, the OO-approach is dynamic and might be adjusted according to the selected outcome variable of interest.

Two distinct regions in 2q24.2-q24.3 associated with idiopathic epilepsy

P>Approximately 50% of all carriers of 2q21-q31 deletions present epileptic seizures. The band 2q24 constitutes the smallest commonly deleted segment in these patients, and contains the voltage-gated sodium channel genes SCN1A and SCN2A, associated with Dravet syndrome and benign familial neonatal-infantile seizures, respectively. A further putative locus involving epilepsy in the region was previously identified through disruption of the SLC4A10 gene by translocation. In the course of performing high-resolution DNA copy number analyses on syndromic mentally impaired individuals, we encountered three patients with overlapping deletions in chromosome region 2q24. Two of these patients exhibited epileptic seizures in addition to mental deficiency. The deletion in one of the epileptic patients did not include the SCN cluster, demonstrating that a less severe form of epilepsy maps to an adjacent genomic region. This second region comprises about 3 Mb and contains the candidate gene SLC4A10, providing further support for the potential role of this gene in epilepsy.

A study of adrenocortical tumors in children and adolescents by a comparative genomic hybridization technique

Adrenocortical tumors (ACT) are rare neoplasms of the adrenal glands accounting for 0.2% of all pediatric cancers. However, the incidence of ACT in South Brazilian children is 10 to 15 times greater than the worldwide incidence. Comparative genomic hybridization studies have revealed the presence of a high degree of chromosomal instability in ACT. We evaluated 16 ACT, 8 of them carcinomas and 8 adenomas. The presence of changes in DNA copy numbers was determined by comparative genomic hybridization, and the findings were validated by real-time polymerase chain reaction on the basis of IGF-II gene expression. The adenomas showed a mean of 19.7 imbalances per case, with the most frequent gain and loss being 4p15.1-p15.3 and 20p11.2-p13.2, respectively. The carcinomas presented with a mean of 35.5 imbalances per case, with the more frequent gain being 2q14.1-q24.3 and the more frequent losses being 3q21-q26.2, 20q12-qter, and 22q11.2-q13.3. The most frequent imbalances in both adenomas and carcinomas were gains of 1p21-p31.2, 2p12-p21 and loss of 20p11.2-p12. The expression of IGF-II mRNA (11p15.5) was higher in samples that presented with a gain of this region. It has been established that great genomic instability exists in pediatric ACT.