Comparative sequence analysis of the non-protein-coding mitochondrial DNA of inbred rat strains

Characterisation of the Complete Mitochondrial Genome of Critically Endangered Mustela lutreola (Carnivora: Mustelidae) and Its Phylogenetic and Conservation Implications

In this paper, a complete mitochondrial genome of the critically endangered European mink Mustela lutreola L., 1761 is reported. The mitogenome was 16,504 bp in length and encoded the typical 13 protein-coding genes, two ribosomal RNA genes and 22 transfer RNA genes, and harboured a putative control region. The A+T content of the entire genome was 60.06% (A > T > C > G), and the AT-skew and GC-skew were 0.093 and −0.308, respectively. The encoding-strand identity of genes and their order were consistent with a collinear gene order characteristic for vertebrate mitogenomes. The start codons of all protein-coding genes were the typical ATN. In eight cases, they were ended by complete stop codons, while five had incomplete termination codons (TA or T). All tRNAs had a typical cloverleaf secondary structure, except tRNA^Ser(AGC) and tRNA^Lys, which lacked the DHU stem and had reduced DHU loop, respectively. Both rRNAs were capable of folding into complex secondary structures, containing unmatched base pairs. Eighty-one single nucleotide variants (substitutions and indels) were identified. Comparative interspecies analyses confirmed the close phylogenetic relationship of the European mink to the so-called ferret group, clustering the European polecat, the steppe polecat and the black-footed ferret. The obtained results are expected to provide useful molecular data, informing and supporting effective conservation measures to save M. lutreola.

Single nanomolar doxorubicin exposure triggers compensatory mitochondrial responses in H9c2 cardiomyoblasts

Dose-dependent and cumulative cardiotoxicity associated with doxorubicin (DOX) is the main limitation of anticancer therapy. Pediatric cancer survivors are particularly vulnerable, and no effective prevention measures are available. The aim of the present study was to investigate the persistent effects of nanomolar DOX concentrations and determine whether a pretreatment would induce mitochondrial adaptations in H9c2 cardiomyoblasts. H9c2 cells were incubated with DOX (10 and 25 nM) for 24 h, followed by 9 days of recovery in drug-free medium. We found that the sub-therapeutic DOX treatment induced persistent hypertrophy and dose-dependent cell cycle arrest in G2/M. Glycolytic activity, indirectly based on extracellular acidification rate, and basal respiration were significantly decreased in DOX-treated cells compared to controls, although both groups showed similar maximal respiration. Additionally, nanomolar DOX pretreatment resulted in upregulation of mitochondrial DNA transcripts accompanied by a decrease in DNA methyltransferase 1 (DNMT1) and global methylation levels. Finally, the pretreatment with DOX ameliorated H9c2 cells resistance against a subsequent exposure to DOX. These results suggest that nanomolar DOX pretreatment induced a beneficial and possibly epigenetic-based mitochondrial adaptation, raising the possibility that an early sub-therapeutic DOX treatment can be used as a preconditioning and protective approach during anticancer therapies.

Droplet digital PCR shows the D-Loop to be an error prone locus for mitochondrial DNA copy number determination

Absolute quantification of mitochondrial DNA copy number (mCN) provides important insights in many fields of research including cancer, cardiovascular and reproductive health. Droplet digital PCR (ddPCR) natively reports absolute copy number, and we have developed a single-dye, multiplex assay to measure rat mCN that is accurate, precise and affordable. We demonstrate simple methods to optimize this assay and to determine nuclear reference pseudogene copy number to extend the range of mCN that can be measured with this assay. We evaluated two commonly used mitochondrial DNA reference loci to determine mCN, the ND1 gene and the D-Loop. Harnessing the absolute measures of ddPCR, we found that the D-Loop amplifies with a copy number of ~1.0–1.5 relative to other sites on the mitochondrial genome. This anomalous copy number varied significantly between rats and tissues (aorta, brain, heart, liver, soleus muscle). We advocate for avoiding the D-Loop as a mitochondrial reference in future studies of mCN. Further, we report a novel approach to quantifying immunolabelled mitochondrial DNA that provides single-cell estimates of mCN that closely agree with the population analyses by ddPCR. The combination of these assays represents a cost-effective and powerful suite of tools to study mCN.

Genomic analyses identify multiple Asian origins and deeply diverged mitochondrial clades in inbred brown rats (Rattus norvegicus)

Over 500 strains of inbred brown rats (Rattus norvegicus) have been developed for use as a biomedical model organism. Most of these inbred lines were derived from the colony established at the Wistar Institute in 1906 or its descendants following worldwide distribution to research and breeding centers. The geographic source of the animals that founded the Wistar colony has been lost to history; thus, we compared 25 inbred rat strains to 326 wild rats from a global diversity dataset at 32 k SNPs, and 47 mitochondrial genomes to identify the source populations. We analyzed nuclear genomic data using principal component analyses and co-ancestry heat maps, and mitogenomes using phylogenetic trees and networks. In the nuclear genome, inbred rats clustered together indicating a single geographic origin for the strains studied and showed admixed ancestral variation with wild rats in eastern Asia and western North America. The Sprague Dawley derived, Wistar derived, and Brown Norway strains each had mitogenomes from different clades which diverged between 13 and 139 kya. Thus, we posit that rats originally collected for captive breeding had high mitochondrial diversity that became fixed through genetic drift and/or artificial selection. Our results show that these important medical models share common genomic ancestry from a few source populations, and opportunities exist to create new strains with diverse genomic backgrounds to provide novel insight into the genomic basis of disease phenotypes.

Polyamines in the life of Arabidopsis: profiling the expression of S-adenosylmethionine decarboxylase (SAMDC) gene family during its life cycle

BACKGROUND

Arabidopsis has 5 paralogs of the S-adenosylmethionine decarboxylase (SAMDC) gene. Neither their specific role in development nor the role of positive/purifying selection in genetic divergence of this gene family is known. While some data are available on organ-specific expression of AtSAMDC1, AtSAMDC2, AtSAMDC3 and AtSAMDC4, not much is known about their promoters including AtSAMDC5, which is believed to be non-functional.

RESULTS

(1) Phylogenetic analysis of the five AtSAMDC genes shows similar divergence pattern for promoters and coding sequences (CDSs), whereas, genetic divergence of 5'UTRs and 3'UTRs was independent of the promoters and CDSs; (2) while AtSAMDC1 and AtSAMDC4 promoters exhibit high activity (constitutive in the former), promoter activities of AtSAMDC2, AtSAMDC3 and AtSAMDC5 are moderate to low in seedlings (depending upon translational or transcriptional fusions), and are localized mainly in the vascular tissues and reproductive organs in mature plants; (3) based on promoter activity, it appears that AtSAMDC5 is both transcriptionally and translationally active, but based on it's coding sequence it seems to produce a non-functional protein; (4) though 5'-UTR based regulation of AtSAMDC expression through upstream open reading frames (uORFs) in the 5'UTR is well known, no such uORFs are present in AtSAMDC4 and AtSAMDC5; (5) the promoter regions of all five AtSAMDC genes contain common stress-responsive elements and hormone-responsive elements; (6) at the organ level, the activity of AtSAMDC enzyme does not correlate with the expression of specific AtSAMDC genes or with the contents of spermidine and spermine.

CONCLUSIONS

Differential roles of positive/purifying selection were observed in genetic divergence of the AtSAMDC gene family. All tissues express one or more AtSAMDC gene with significant redundancy, and concurrently, there is cell/tissue-specificity of gene expression, particularly in mature organs. This study provides valuable information about AtSAMDC promoters, which could be useful in future manipulation of crop plants for nutritive purposes, stress tolerance or bioenergy needs. The AtSAMDC1 core promoter might serve the need of a strong constitutive promoter, and its high expression in the gametophytic cells could be exploited, where strong male/female gametophyte-specific expression is desired; e.g. in transgenic modification of crop varieties.

Mitogenome Sequencing in the Genus Camelus Reveals Evidence for Purifying Selection and Long-term Divergence between Wild and Domestic Bactrian Camels

The genus Camelus is an interesting model to study adaptive evolution in the mitochondrial genome, as the three extant Old World camel species inhabit hot and low-altitude as well as cold and high-altitude deserts. We sequenced 24 camel mitogenomes and combined them with three previously published sequences to study the role of natural selection under different environmental pressure, and to advance our understanding of the evolutionary history of the genus Camelus. We confirmed the heterogeneity of divergence across different components of the electron transport system. Lineage-specific analysis of mitochondrial protein evolution revealed a significant effect of purifying selection in the concatenated protein-coding genes in domestic Bactrian camels. The estimated dN/dS < 1 in the concatenated protein-coding genes suggested purifying selection as driving force for shaping mitogenome diversity in camels. Additional analyses of the functional divergence in amino acid changes between species-specific lineages indicated fixed substitutions in various genes, with radical effects on the physicochemical properties of the protein products. The evolutionary time estimates revealed a divergence between domestic and wild Bactrian camels around 1.1 [0.58-1.8] million years ago (mya). This has major implications for the conservation and management of the critically endangered wild species, Camelus ferus.

Identification of candidate gene regions in the rat by co-localization of QTLs for bone density, size, structure and strength

Susceptibility to osteoporotic fracture is influenced by genetic factors that can be dissected by whole-genome linkage analysis in experimental animal crosses. The aim of this study was to characterize quantitative trait loci (QTLs) for biomechanical and two-dimensional dual-energy X-ray absorptiometry (DXA) phenotypes in reciprocal F2 crosses between diabetic GK and normo-glycemic F344 rat strains and to identify possible co-localization with previously reported QTLs for bone size and structure. The biomechanical measurements of rat tibia included ultimate force, stiffness and work to failure while DXA was used to characterize tibial area, bone mineral content (BMC) and areal bone mineral density (aBMD). F2 progeny (108 males, 98 females) were genotyped with 192 genome-wide markers followed by sex- and reciprocal cross-separated whole-genome QTL analyses. Significant QTLs were identified on chromosome 8 (tibial area; logarithm of odds (LOD) = 4.7 and BMC; LOD = 4.1) in males and on chromosome 1 (stiffness; LOD = 5.5) in females. No QTLs showed significant sex-specific interactions. In contrast, significant cross-specific interactions were identified on chromosome 2 (aBMD; LOD = 4.7) and chromosome 6 (BMC; LOD = 4.8) for males carrying F344mtDNA, and on chromosome 15 (ultimate force; LOD = 3.9) for males carrying GKmtDNA, confirming the effect of reciprocal cross on osteoporosis-related phenotypes. By combining identified QTLs for biomechanical-, size- and qualitative phenotypes (pQCT and 3D CT) from the same population, overlapping regions were detected on chromosomes 1, 3, 4, 6, 8 and 10. These are strong candidate regions in the search for genetic risk factors for osteoporosis.

Genetic monitoring detects an overlooked cryptic species and reveals the diversity and distribution of three invasive Rattus congeners in South Africa

Background

South Africa's long and extensive trade activity has ensured ample opportunities for exotic species introduction. Whereas the rich biodiversity of endemic southern African fauna has been the focus of many studies, invasive vertebrates are generally overlooked despite potential impacts on biodiversity, health and agriculture. Genetic monitoring of commensal rodents in South Africa which uncovered the presence of Rattus tanezumi, a South-East Asian endemic not previously known to occur in Africa, provided the impetus for expanded studies on all invasive Rattus species present.

Results

To this end, intensified sampling at 28 South African localities and at one site in Swaziland, identified 149 Rattus specimens. Cytochrome b gene sequencing revealed the presence of two R. tanezumi, seven Rattus rattus and five Rattus norvegicus haplotypes in south Africa. Phylogenetic results were consistent with a single, recent R. tanezumi introduction and indicated that R. norvegicus and R. rattus probably became established following at least two and three independent introductions, respectively. Intra- and inter-specific diversity was highest in informal human settlements, with all three species occurring at a single metropolitan township site. Rattus norvegicus and R. rattus each occurred sympatrically with Rattus tanezumi at one and five sites, respectively. Karyotyping of selected R. rattus and R. tanezumi individuals identified diploid numbers consistent with those reported previously for these cryptic species. Ordination of bioclimatic variables and MaxEnt ecological niche modelling confirmed that the bioclimatic niche occupied by R. tanezumi in south Africa was distinct from that occupied in its naturalised range in south-east Asia suggesting that factors other than climate may influence the distribution of this species.

Conclusions

This study has highlighted the value of genetic typing for detecting cryptic invasive species, providing historical insights into introductions and for directing future sampling. The apparent ease with which a cryptic species can become established signals the need for broader implementation of genetic monitoring programmes. In addition to providing baseline data and potentially identifying high-risk introduction routes, the predictive power of ecological niche modelling is enhanced when species records are genetically verified.

Mitochondrial polymorphisms in rat genetic models of hypertension

Hypertension is a complex trait that has been studied extensively for genetic contributions of the nuclear genome. We examined mitochondrial genomes of the hypertensive strains: the Dahl Salt-Sensitive (S) rat, the Spontaneously Hypertensive Rat (SHR), and the Albino Surgery (AS) rat, and the relatively normotensive strains: the Dahl Salt-Resistant (R) rat, the Milan Normotensive Strain (MNS), and the Lewis rat (LEW). These strains were used previously for linkage analysis for blood pressure (BP) in our laboratory. The results provide evidence to suggest that variations in the mitochondrial genome do not account for observed differences in blood pressure between the S and R rats. However, variants were detected among the mitochondrial genomes of the various hypertensive strains, S, SHR, and AS, and also among the normotensive strains R, MNS, and LEW. A total of 115, 114, 106, 106, and 16 variations in mtDNA were observed between the comparisons S versus LEW, S versus MNS, S versus SHR, S versus AS, and SHR versus AS, respectively. Among the 13 genes coding for proteins of the electron transport chain, 8 genes had nonsynonymous variations between S, LEW, MNS, SHR, and AS. The lack of any sequence variants between the mitochondrial genomes of S and R rats provides conclusive evidence that divergence in blood pressure between these two inbred strains is exclusively programmed through their nuclear genomes. The variations detected among the various hypertensive strains provides the basis to construct conplastic strains and further evaluate the effects of these variants on hypertension and associated phenotypes.