Search

The A-to-G point mutation at position 3243 in the human mitochondrial genome (m.3243A > G) is the most common pathogenic mtDNA variant responsible for disease in humans. It is widely accepted that m.3243A > G levels decrease in blood with age, and an age correction representing ~ 2% annual decline is often applied to account for this change in mutation level. Here we report that recent data indicate that the dynamics of m.3243A > G are more complex and depend on the mutation level in blood in a bi-phasic way. Consequently, the traditional 2% correction, which is adequate ‘on average’, creates opposite predictive biases at high and low mutation levels. Unbiased age correction is needed to circumvent these drawbacks of the standard model. We propose to eliminate both biases by using an approach where age correction depends on mutation level in a biphasic way to account for the dynamics of m.3243A > G in blood. The utility of this approach was further tested in estimating germline selection of m.3243A > G. The biphasic approach permitted us to uncover patterns consistent with the possibility of positive selection for m.3243A > G. Germline selection of m.3243A > G shows an ‘arching’ profile by which selection is positive at intermediate mutant fractions and declines at high and low mutant fractions. We conclude that use of this biphasic approach will greatly improve the accuracy of modelling changes in mtDNA mutation frequencies in the germline and in somatic cells during aging.

Every cell in our body contains a vibrant population of mitochondria, or, more precisely, of mitochondrial DNA molecules (mtDNAs). Just like members of any population mtDNAs multiply (by replication) and ‘die’ (i.e., are removed, either by degradation or by distribution into the sister cell in mitosis). An intriguing question is whether all mitochondria in this population are equal, especially whether some are responsible primarily for reproduction and some - for empowering the various jobs of the mitochondrion, oxidative phosphorylation in the first place. Importantly, because mtDNA is highly damaged such a separation of responsibilities could help greatly reduce the conversion of DNA damage into real inheritable mutations. An unexpected twist in the resolution of this problem has been brought about by a recent high-precision analysis of mtDNA mutations (Sanchez-Contreras et al. 2023). They discovered that certain transversion mutations, unlike more common transitions, are not accumulating with age in mice. We argue that this observation requires the existence of a permanent replicating subpopulation/lineage of mtDNA molecules, which are protected from DNA damage, a.k.a. the ‘stem’ mtDNA. This also implies the existence of its antipode i.e., the ‘worker’ mtDNA, which empowers OSPHOS, sustains damage and rarely replicates. The analysis of long HiFi reads of mtDNA performed by PacBio closed circular sequencing confirms this assertion.

The process of domestication is associated with decrease in effective population size, which in turn leads to accumulation of slightly-deleterious mutations due to genetic drift. To maintain genome quality at a high level, we propose to use a stress-induced strong purifying selection, which based on negative epistasis, can effectively eliminate organisms with an excess of deleterious variants. Here, to identify stress factors, which interact with the effect of deleterious mutations we performed a proof-of-principle experiment with several regimes of a heat shock. We observed that fitness of mutated versus wild-type carp lines drops stronger after heat shock, which is a signature of a negative epistasis. Although the observed trend is promising, the effect of the epistasis is weak and unstable from family to family. Thus, more deep tuning of heat shock regimes is needed to uncover the most efficient combination of factors (absolute temperature, duration, stage of the embryo development) aggravating the burden of deleterious mutations and thus exposing them to the selection.

The resilience of the mitochondrial genome (mtDNA) to a high mutational pressure depends, in part, on negative purifying selection in the germline. A paradigm in the field has been that such selection, at least in part, takes place in primordial germ cells (PGCs). Specifically, Floros et al. (Nature Cell Biology20: 144–51) reported an increase in the synonymity of mtDNA mutations (a sign of purifying selection) between pooled early-stage and late-stage PGCs. We re-analyzed Floros’ et al. pooled PGC data and noticed that their mutational dataset was significantly contaminated with single nucleotide variants (SNVs) derived from a nuclear sequence of mtDNA origin (NUMT) located on chromosome 5. Contamination was caused by co-amplification of the NUMT sequence by cross-specific PCR primers. Importantly,when we removed NUMT-derived SNVs, the evidence of purifying selection was abolished. In addition to pooled PGCs, Floros et al. reported the analysis ofsinglelate-stage PGCs, which were amplified with different sets of PCR primers that cannot amplify the NUMT sequence. Accordingly, we found no NUMT-derived SNVs among single PGCs mutations. Interestingly, single PGC mutations show adecreaseof synonymity with increased intracellular mutant fraction. This pattern is incompatible with predominantly negative selection. This suggests that germline selection of mtDNA mutations is a complex phenomenon and that the part of this process that takes place in PGCs may be predominantly positive. However counterintuitive, positive germline selection of detrimental mtDNA mutations has been reported previously and potentially may be evolutionarily advantageous.

Heat shock proteins in parallel with their main and originally discovered function – maintenance of folded proteins under stressful conditions, can play also background buffering role – by folding proteins with an excess of slightly-deleterious nonsynonymous variants (SDNV). Here we tested several scenarios of this buffering role. On the comparative species scale, we demonstrated that low-Ne species are characterized by a higher expression level of hsp90 which can be explained by the excess of SDNV. On the comparative tissue level, we showed that long-lived tissues have also a higher hsp90 expression level, which can be advantageous to maintain the functionality of proteins. On the comparative gene level, we demonstrated that purifying selection of hsp90 in low-Ne-species did not relax as strongly as it happens for control genes, similar to hsp90. Additionally, we demonstrated that hsp clients versus non-clients are characterised by decreased level of selective constraints; demonstrate stronger relaxation of purifying selection in low-Ne species; have an excess of slightly-deleterious variants associated with complex disease phenotypes in humans; have an excess of pathological variants associated with clinical phenotypes in humans, suggesting that clients, being buffered by hsp90 can degenerate a bit more as compared to non-clients. Altogether, our results show that the secondary role of hsp, buffering of SDNV, is widespread and universal affecting properties of species, tissues and genes. A deep understanding of the buffering role of hsp90 will help to predict the deleterious effect of each variant in the human genome more precisely as well as will extend the application of the effectively-neutral theory of molecular evolution.

Using low-coverage, whole-genome sequences of trophectoderm biopsies from 11,610 human blastocyst-stage embryos, we analyzed the relationship between chromosomal abnormalities and mitochondrial (mt) DNA dynamics. Comparing 6,208 aneuploid and 5,402 euploid embryos in cohort studies, we found that mtDNA content in aneuploid embryos was significantly higher than that in euploid embryos. This outcome was confirmed through intrafamilial analyses of embryos with matched parents andin vitrofertilization cycles, and it occurred independent of maternal age. Additional human population-based studies uncovered a higher abundance of ultra-rare mtDNA variants located in never-altered positions in the human population in aneuploid compared to euploid embryos in both cohort- and family-based analyses. This maternal age-independent association of increased mtDNA content and aneuploidy in human embryos may reflect a novel mechanism of purifying selection against potentially deleterious mtDNA variants, which arise from germline or early developmental mtDNA damaging events, that occurs in human embryos prior to implantation.

The mutational spectrum of the mitochondrial DNA (mtDNA) does not resemble any of the known mutational signatures of the nuclear genome and variation in mtDNA mutational spectra between different organisms is still incomprehensible. Since mitochondria are responsible for aerobic respiration, it is expected that mtDNA mutational spectrum is affected by oxidative damage. Assuming that oxidative damage increases with age, we analyse mtDNA mutagenesis of different species in regards to their generation length. Analysing, (i) dozens of thousands of somatic mtDNA mutations in samples of different ages (ii) 70053 polymorphic synonymous mtDNA substitutions reconstructed in 424 mammalian species with different generation lengths and (iii) synonymous nucleotide content of 650 complete mitochondrial genomes of mammalian species we observed that the frequency of AH > GH substitutions (H: heavy strand notation) is twice bigger in species with high versus low generation length making their mtDNA more AH poor and GH rich. Considering that AH > GH substitutions are also sensitive to the time spent single-stranded (TSSS) during asynchronous mtDNA replication we demonstrated that AH > GH substitution rate is a function of both species-specific generation length and position-specific TSSS. We propose that AH > GH is a mitochondria-specific signature of oxidative damage associated with both aging and TSSS.

Evolution is a function of mutagenesis and selection. To analyse the role of mutagenesis on the structure of the SARS-CoV-2 genome, we reconstructed the mutational spectrum, which was highly C>U and G>U biased. This bias forces the SARS-CoV-2 genome to become increasingly U-rich unless selection cancels it. We analysed the consequences of this bias on the composition of the most neutral (four-fold degenerate synonymous substitutions) and the least neutral positions (nonsynonymous substitutions). The neutral nucleotide composition is already highly saturated by U and, according to our model, it is at equilibrium, suggesting that in the future, we don’t expect any more increase in U. However, nonsynonymous changes continue slowly evolve towards equilibrium substituting CG-rich amino-acids (“losers”) with U-rich ones (“gainers”). This process is universal for all genes of SARS-CoV-2 as well as for other coronaviridae species. In line with the direction mutation pressure hypothesis, we show that viral-specific amino acid content is associated with the viral-specific mutational spectrum due to the accumulation of effectively neutral slightly deleterious variants (losers to gainers) during the molecular evolution. The tuning of a protein space by the mutational process is expected to be typical for species with relaxed purifying selection, suggesting that the purging of slightly-deleterious variants in the SARS-CoV-2 population is not very effective, probably due to the fast expansion of the viral population during the pandemic. Understanding the mutational process can help to design more robust vaccines, based on gainer-rich motifs, close to the mutation-selection equilibrium.

A wealth of specialized neuroendocrine command systems intercalated within the hypothalamus control the most fundamental physiological needs in vertebrates1,2. Nevertheless, we lack a developmental blueprint that integrates the molecular determinants of neuronal and glial diversity along temporal and spatial scales of hypothalamus development3. Here we combine single-cell RNA sequencing of 51,199 mouse cells of ectodermal origin, gene regulatory network (GRN) screens in conjunction with genome-wide association study-based disease phenotyping, and genetic lineage reconstruction to show that nine glial and thirty-three neuronal subtypes are generated by mid-gestation under the control of distinct GRNs. Combinatorial molecular codes that arise from neurotransmitters, neuropeptides and transcription factors are minimally required to decode the taxonomical hierarchy of hypothalamic neurons. The differentiation of γ-aminobutyric acid (GABA) and dopamine neurons, but not glutamate neurons, relies on quasi-stable intermediate states, with a pool of GABA progenitors giving rise to dopamine cells4. We found an unexpected abundance of chemotropic proliferation and guidance cues that are commonly implicated in dorsal (cortical) patterning5 in the hypothalamus. In particular, loss of SLIT-ROBO signalling impaired both the production and positioning of periventricular dopamine neurons. Overall, we identify molecular principles that shape the developmental architecture of the hypothalamus and show how neuronal heterogeneity is transformed into a multimodal neural unit to provide virtually infinite adaptive potential throughout life.

The hypothesis that the evolution of humans involves hybridization between diverged species has been actively debated in recent years. We present the following novel evidence in support of this hypothesis: the analysis of nuclear pseudogenes of mtDNA (“NUMTs”). NUMTs are considered “mtDNA fossils” as they preserve sequences of ancient mtDNA and thus carry unique information about ancestral populations. Our comparison of a NUMT sequence shared by humans, chimpanzees, and gorillas with their mtDNAs implies that, around the time of divergence between humans and chimpanzees, our evolutionary history involved the interbreeding of individuals whose mtDNA had diverged as much as ~4.5 Myr prior. This large divergence suggests a distant interspecies hybridization. Additionally, analysis of two other NUMTs suggests that such events occur repeatedly. Our findings suggest a complex pattern of speciation in primate/human ancestors and provide one potential explanation for the mosaic nature of fossil morphology found at the emergence of the hominin lineage. A preliminary version of this manuscript was uploaded to the preprint server BioRxiv in 2017 (10.1101/134502).

The resilience of the mitochondrial genome to a high mutational pressure depends, in part, on purifying selection against detrimental mutations in the germline. It is crucial to understand the mechanisms of this process. Recently, Floros et al. concluded that much of the purifying selection takes place during the proliferation of primordial germ cells (PGCs) because, according to their analysis, the synonymity of mutations in late PGCs was seemingly increased compared to those in early PGCs. We re-analyzed the Floros et al. mutational data and discovered a high proportion of sequence variants that are not true mutations, but originate from NUMTs, the latter of which are segments of mitochondrial DNA (mtDNA) inserted into nuclear DNA, up to millions of years ago. This is a well-known artifact in mtDNA mutational analysis. Removal of these artifacts from the Floros et al. dataset abolishes the reported effect of purifying selection in PGCs. We therefore conclude that the mechanism of germline selection of mtDNA mutations remains open for debate, and more research is needed to fully elucidate the timing and nature of this process.

The variation in the mutational spectrum of the mitochondrial genome (mtDNA) among species is not well understood. Recently, we demonstrated an increase in A>G substitutions on a heavy chain (hereafter AH>GH) of mtDNA in aged mammals, interpreting it as a hallmark of age-related oxidative damage. In this study, we hypothesized that the occurrence of AH>GHsubstitutions may depend on the level of aerobic metabolism, which can be inferred from an organism’s body temperature. To test this hypothesis, we used body temperature in endotherms and environmental temperature in ectotherms as proxies for metabolic rate and reconstructed mtDNA mutational spectra for 1350 vertebrate species. Our results showed that temperature was associated with increased rates of AH>GHand asymmetry of AH>GHin different species of ray-finned fishes and within geographically distinct clades of European anchovy. Analysis of nucleotide composition in the most neutral synonymous sites of fishes revealed that warm-water species were expectedly more A-poor and G-rich compared to cold-water species. Finally, we extended our analyses to all vertebrates and observed higher AH>GHand increased asymmetry of AH>GHin warm-blooded (mammals and birds) compared to cold-blooded (Actinopterygii, amphibia, reptilia) vertebrate classes. Overall, our findings suggest that temperature, through its influence on metabolism and oxidative damage, shapes the mutational properties and nucleotide content of the mtDNA in all vertebrates.

Mitochondria are the central players in cellular energy production. Mitochondrial DNA (mtDNA) encodes core subunits of oxidative phosphorylation complex (Anderson et al. 1981) and thus it is expect...

Mitochondria is a powerhouse of all eukaryotic cells that have its own circular DNA (mtDNA) encoding various RNAs and proteins. Somatic perturbations of mtDNA are accumulating with age thus it is of great importance to uncover the main sources of mtDNA instability. Recent analyses demonstrated that somatic mtDNA deletions depend on imperfect repeats of various nature between distant mtDNA segments. However, till now there are no comprehensive databases annotating all types of imperfect repeats in numerous species with sequenced complete mitochondrial genome as well as there are no algorithms capable to call all types of imperfect repeats in circular mtDNA. We implemented naive algorithm of pattern recognition by analogy to standard dot-plot construction procedures allowing us to find both perfect and imperfect repeats of four main types: direct, inverted, mirror and complementary. Our algorithm is adapted to specific characteristics of mtDNA such as circularity and an excess of short repeats - it calls imperfect repeats starting from the length of 10 b.p. We constructed interactive web available database ImtRDB depositing perfect and imperfect repeats positions in mtDNAs of more than 3500 Vertebrate species. Additional tools, such as visualization of repeats within a genome, comparison of repeat densities among different genomes and a possibility to download all results make this database useful for many biologists. Our first analyses of the database demonstrated that mtDNA imperfect repeats (i) are usually short; (ii) associated with unfolded DNA structures; (iii) four types of repeats positively correlate with each other forming two equivalent pairs: direct and mirror versus inverted and complementary, with identical nucleotide content and similar distribution between species; (iv) abundance of repeats is negatively associated with GC content; (v) dinucleotides GC versus CG are overrepresented on light chain of mtDNA covered by repeats. ImtRDB is available at http://bioinfodbs.kantiana.ru/ImtRDB/ . It is accompanied by the software calling all types of interspersed repeats with different level of degeneracy in circular DNA. This database and software can become a very useful tool in various areas of mitochondrial and chloroplast DNA research.

Mitochondrial DNA (mtDNA) encodes core subunits of oxidative phosphorylation complexes and, as a result of intricate regulatory crosstalk between nuclear and mitochondrial genomes, the total number of mtDNA copies fits the requirements of each cell type. Deviations from the physiological number of mtDNA copies are expected to be deleterious and might cause some inherited diseases and normal ageing. We studied 46 obese patients with type 2 diabetes (T2DM) one year after a laparoscopic sleeve gastrectomy (LSG) and Roux-en-Y gastric bypass (RYGB). The results were compared with normal-weight patients without T2DM (control group 1) (body mass index (BMI) = 22.5 ± 3.01 kg/m 2 ) and patients with obesity without T2DM (control group 2) (BMI = 36 ± 3.45 kg/m 2 ). We detected an increase of mtDNA copy number in the cells of the buffy coat obtained from peripheral blood, sampled one year after bariatric surgery. We also found that average mtDNA copy number as well as its dynamics (before and after the surgery) are gender-specific. To the best of our knowledge, this is the first evidence for the restoration of mtDNA copy number in obese patients after LSG and RYGB.

Aging in postmitotic tissues is associated with clonal expansion of somatic mitochondrial deletions, the origin of which is not well understood. Deletions in mitochondrial DNA (mtDNA) are often flanked by direct nucleotide repeats, but this alone does not fully explain their distribution. Here, we hypothesized that the close proximity of direct repeats on single-stranded DNA might play a role in the formation of deletions. By analyzing human mtDNA deletions in the major arc of mtDNA, which is single-stranded during replication and is characterized by a high number of deletions, we found a non-uniform distribution with a "hot spot" where one deletion breakpoint occurred within the region of 6-9kb and another within 13-16kb of the mtDNA. This distribution was not explained by the presence of direct repeats, suggesting that other factors, such as the spatial proximity of these two regions can be the cause. In silico analyses revealed that the single-stranded major arc may be organized as a large-scale hairpin-like loop with a center close to 11kb and contacting regions between 6-9 kb and 13-16 kb, which would explain the high deletion activity in this contact zone. The direct repeats located within the contact zone, such as the well-known common repeat with a first arm at 8470-8482 bp and a second arm at 13447-13459 bp, are three times more likely to cause deletions compared to direct repeats located outside of the contact zone. An analysis of age- and disease-associated deletions demonstrated that the contact zone plays a crucial role in explaining the age-associated deletions, emphasizing its importance in the rate of healthy aging. Overall, we provide topological insights into the mechanism of age-associated deletion formation in human mtDNA, which could be used to predict somatic deletion burden and maximum lifespan in different human haplogroups and mammalian species.

Mutational spectrum of the mitochondrial genome (mtDNA) does not resemble signatures of any known mutagens and variation in mtDNA mutational spectra between different tissues and organisms is still incomprehensible. Since mitochondria is tightly involved in aerobic energy production, it is expected that mtDNA mutational spectra may be affected by the oxidative damage which is increasing with cellular and organismal aging. However, the well-documented mutational signature of the oxidative damage, G>T substitutions, is typical only for the nuclear genome while it is extremely rare and age-independent in mtDNA. Thus it is still unclear if there is a mitochondria - specific mutational signature of the oxidative damage. Here, reconstructing mtDNA mutational spectra for human cancers originated from 21 tissues with various cell turnover rate, human oocytes fertilized at different ages, and 424 mammalian species with variable generation length which is a proxy for oocyte age, we observed that the frequency of AH>GHsubstitutions (H- heavy chain notation) is positively correlated with cellular and organismal longevity. Moreover, this mutational bias from AHto GHaffects nucleotide content at the fourfold degenerative synonymous positions leading to a deficit of AHand excess of GH, which is especially pronounced in long-lived mammals. Taking into account additionally, that AH>GHis sensitive to time being single stranded during mtDNA asynchronous replication and A>G is associated with oxidative damage of single-stranded DNA in recent bacterial experiments we propose that AH>GHis a mutational signature of oxidative damage in mtDNA.

International audience; Copy-number changes in 16p11.2 contribute significantly to neuropsychiatric traits. Besides the 600 kb BP4-BP5 CNV found in 0.5%-1% of individuals with autism spectrum disorders and schizophrenia and whose rearrangement causes reciprocal defects in head size and body weight, a second distal 220 kb BP2-BP3 CNV is likewise a potent driver of neuropsychiatric, anatomical, and metabolic pathologies. These two CNVs are engaged in complex reciprocal chromatin looping, intimating a functional relationship between genes in these regions that might be relevant to pathomechanism. We assessed the drivers of the distal 16p11.2 duplication by overexpressing each of the nine encompassed genes in zebrafish. Only overexpression of LAT induced a reduction of brain proliferating cells and concomitant microcephaly. Consistently, suppression of the zebrafish ortholog induced an increase of proliferation and macrocephaly. These phenotypes were not unique to zebrafish; Lat knockout mice show brain volumetric changes. Consistent with the hypothesis that LAT dosage is relevant to the CNV pathology, we observed similar effects upon overexpression of CD247 and ZAP70, encoding members of the LAT signalosome. We also evaluated whether LAT was interacting with KCTD13, MVP, and MAPK3, major driver and modifiers of the proximal 16p11.2 600 kb BP4-BP5 syndromes, respectively. Co-injected embryos exhibited an increased microcephaly, suggesting the presence of genetic interaction. Correspondingly, carriers of 1.7 Mb BP1-BP5 rearrangements that encompass both the BP2-BP3 and BP4-BP5 loci showed more severe phenotypes. Taken together, our results suggest that LAT, besides its well-recognized function in T cell development, is a major contributor of the 16p11.2 220 kb BP2-BP3 CNV-associated neurodevelopmental phenotypes.

The majority of aneuploid fetuses are spontaneously miscarried. Nevertheless, some aneuploid individuals survive despite the strong genetic insult. Here, we investigate if the survival probability of aneuploid fetuses is affected by the genome-wide burden of slightly deleterious variants. We analyzed two cohorts of live-born Down syndrome individuals (388 genotyped samples and 16 fibroblast transcriptomes) and observed a deficit of slightly deleterious variants on Chromosome 21 and decreased transcriptome-wide variation in the expression level of highly constrained genes. We interpret these results as signatures of embryonic selection, and propose a genetic handicap model whereby an individual bearing an extremely severe deleterious variant (such as aneuploidy) could escape embryonic lethality if the genome-wide burden of slightly deleterious variants is sufficiently low. This approach can be used to study the composition and effect of the numerous slightly deleterious variants in humans and model organisms.

Genetic factors may play an important role in species extinction but their actual effect remains poorly understood, particularly because of a strong and potentially masking effect expected from ecological traits. We investigated the role of genetics in mammal extinction taking both ecological and genetic factors into account. As a proxy for the role of genetics we used the ratio of the rates of nonsynonymous (amino acid changing) to synonymous (leaving the amino acid unchanged) nucleotide substitutions, Ka / Ks. Because most nonsynonymous substitutions are likely to be slightly deleterious and thus selected against, this ratio is a measure of the inefficiency of selection: if large (but less than 1), it implies a low efficiency of selection against nonsynonymous mutations. As a result, nonsynonymous mutations may accumulate and thus contribute to extinction. As a proxy for the role of ecology we used body mass W, with which most extinction-related ecological traits strongly correlate. As a measure of extinction risk we used species’ affiliation with the five levels of extinction threat according to the IUCN Red List of Threatened Species. We calculated Ka / Ks for mitochondrial protein-coding genes of 211 mammalian species, each of which was characterized by body mass and the level of threat. Using logistic regression analysis, we then constructed a set of logistic regression models of extinction risk on ln(Ka / Ks) and lnW. We found that Ka / Ks and body mass are responsible for a 38% and a 62% increase in extinction risk, respectively. Given that the standard error of these values is 13%, the contribution of genetic factors to extinction risk in mammals is estimated to be one-quarter to one-half of the total of ecological and genetic effects. We conclude that the effect of genetics on extinction is significant, though it is almost certainly smaller than the effect of ecological traits.

The study of gene expression in mammalian single cells via genomic technologies now provides the possibility to investigate the patterns of allelic gene expression. We used single-cell RNA sequencing to detect the allele-specific mRNA level in 203 single human primary fibroblasts over 133,633 unique heterozygous single-nucleotide variants (hetSNVs). We observed that at the snapshot of analyses, each cell contained mostly transcripts from one allele from the majority of genes; indeed, 76.4% of the hetSNVs displayed stochastic monoallelic expression in single cells. Remarkably, adjacent hetSNVs exhibited a haplotype-consistent allelic ratio; in contrast, distant sites located in two different genes were independent of the haplotype structure. Moreover, the allele-specific expression in single cells correlated with the abundance of the cellular transcript. We observed that genes expressing both alleles in the majority of the single cells at a given time point were rare and enriched with highly expressed genes. The relative abundance of each allele in a cell was controlled by some regulatory mechanisms given that we observed related single-cell allelic profiles according to genes. Overall, these results have direct implications in cellular phenotypic variability.

The importance of natural gene expression variation for human behavior is undisputed, but its impact on circadian physiology remains mostly unexplored. Using umbilical cord fibroblasts, we have determined by genome-wide association how common genetic variation impacts upon cellular circadian function. Gene set enrichment points to differences in protein catabolism as one major source of clock variation in humans. The two most significant alleles regulated expression of COPS7B, a subunit of the COP9 signalosome. We further show that the signalosome complex is imported into the nucleus in timed fashion to stabilize the essential circadian protein BMAL1, a novel mechanism to oppose its proteasome-mediated degradation. Thus, circadian clock properties depend in part upon a genetically-encoded competition between stabilizing and destabilizing forces, and genetic alterations in these mechanisms provide one explanation for human chronotype.

SummaryGenome-wide studies on the genetic basis of gene expression and the structural properties of chromatin have considerably advanced our understanding of the function of the human genome. However, it remains unclear how structure relates to function and, in this work, we aim at bridging both by assembling a dataset that combines the activity of regulatory elements (e.g. enhancers and promoters), expression of genes and genetic variations of 317 individuals and across two cell types. We show that the regulatory activity is structured within 12,583 Cis Regulatory Domains (CRDs) that are cell type specific and highly reflective of the local (i.e. Topologically Associating Domains) and global (i.e. A/B nuclear compartments) nuclear organization of the chromatin. These CRDs essentially delimit the sets of active regulatory elements involved in the transcription of most genes, thereby capturing complex regulatory networks in which the effects of regulatory variants are propagated and combined to finally mediate expression Quantitative Trait Loci. Overall, our analysis reveals the complexity and specificity of cis and trans regulatory networks and their perturbation by genetic variation.

The question: human evolution- gradual process or a rapid discontinuous change? Whether human origin was a gradual process or a result of rapid change has been a focus of intense debate. Of particular interest is the climate change ~2.9-2.5 Ma, thought to have precipitated the separation of the genus Homo (~2.8Ma). The debate mostly concerned continuity/punctuality of the fossil record, but of course the rate of the underlying genetic change is of ultimate interest/importance. Did hominid lineage experience an increased mutation rate when a large number of hominins emerged and eventually gave rise to the split between Australopitecus/Paranthropus and Homo? The obstacle: vague timing of conventional mutations. The difficulty in answering the above question lies in the way past mutations are timed. Conventional point mutations are assigned to specific branches of the DNA-derived phylogenetic trees. The essence of the problem is that mutations can be located within branch segments from branching point to branching point, but the exact position within the segment is principally unknown. Because the hominid DNA-derived phylogenetic tree is rather sparsely populated with branches, the precision of mutation timing is low, e.g., human-specific mutations can be positioned within ~6 My from separation from chimpanzee. The solution: NUMTs – mutations with an internal clock. NUMTs are insertions of mtDNA sequences into the nuclear genome. Unlike point mutation, each NUMTs actually represents a branch on the mtDNA phylogenic tree and thus its time of insertion can be determined as precise as their branching point can be positioned on the tree. In a sense, NUMTs are “mutations with an internal clock”, which is synchronized with the well-established mtDNA mutation evolution clock. By determining the NUMTs’ insertion time points, one can ask whether whether NUMTs were inserted at a constant rate over time or at increased rate during critical periods of evolution, according with the “punctuated evolution” model. Results: Hundreds of pseudogenes have been inserted into the human genome over the last ~60 My of which we considered the last 6 My. Various quality filters resulted in the selection of 18 NUMTs most suitable for phylogenetic analysis. Insertion times of these 18 NUMTs cluster around 2.8Ma. While timing of insertion of NUMTs is imprecise, the observation such a cluster is highly statistically significant. Discussion: It is tempting to hypothesize that accelerated insertion of NUMTs is somehow linked to the speciation process. NUMTs could be either "riders i.e., the rate of insertion could be increased by the overall higher genome flexibility during the speciation period, or "drivers", i.e. they are fixed in the population at increased rate during speciation due to increased selective pressures. If correct, the hypothesis of accelerated pseudogenization would support the idea that evolution of our genus might have been discontinuous.

Introduction: Increasingly, the emergence and evolution of our species is being tied to genetic exchange between divergent lineages within ~1Ma (e.g., Neanderthals, Denisovans). However, little is known about genetic exchange during earlier (pre-1Ma) human evolution and between more divergent lineages. Results: We present evidence of hybridization within human lineage, show that it likely happened between highly divergent (~4.5My) lineages, more than once. We use analysis of nuclear pseudogenes of mtDNA (“NUMTs”). NUMTs are considered “mtDNA fossils”, as they preserve sequences of ancient mtDNA because mutational rate in the nucleus is much lower than in mtDNA. We demonstrate that a NUMT on human chromosome 5, which is shared by chimpanzee and gorilla, had descended from a mitochondrial genome that had been divergent from our ancestor’s mtDNA by ~4.5% at the time of pseudogene insertion. This implies that this pseudogene should have been inserted in a hominid that at that time had been diverged by about 4.5My of evolution from the hominid that at that time carried our mtDNA lineage. In order for this pseudogene and our mtDNA to end up in the same body, these two hominids should have mated with each other. The large divergence implies a distant interspecies (or even inter-generic) hybridization. Additionally, analysis of two other NUMTs (on Chr11 and Chr7) suggests that hybridization events occurred repeatedly. To exclude the large ancestral population size effect we show that mtDNA divergence in extant ape populations does not depend on population size. Discussion: It is thought that within mammals, it takes ~2-4My to establish reproductive isolation. However, fertile inter-generic hybrids have been documented among several primates, separated by ca. 4My. Very recently, hybridization between Colobine genera separated by ~5 My was reported to involve a NUMT scenario similar to what we had proposed human ancestors. Interestingly, phylogenic analysis consistently places the chr5 NUMT insertion around the time of the Homo/Pan split. Intriguingly, certain hominin fossils of that epoch have been interpreted alternately as more human-like or more ape-like. Such morphological mosaicisity could potentially be explained by hybridization. Fixation of NUMTs in question within population should have been rather efficient, since these pseudogenes appear to have been fixed in more than one population. Thus their spread across populations might have been driven by selection. Indeed, NUMTs on chr5 and chr11 are located in 3’ regions of functional genes. Most intriguingly, Ps11 is located 3’ to the RNF141/ZNF230 gene, essential for spermatogenesis. NUMT might have served as an expression modifier for RNF141, resulting in reproductive advantage. Indeed, RNF141 demonstrates selectively driven expression shift in testis of the ancestor of hominines.

The hypothesis that the evolution of humans involved hybridization between diverged species has been actively debated in recent years. We present novel evidence in support of this hypothesis: the analysis of nuclear pseudogenes of mtDNA (“NUMTs”). NUMTs are considered “mtDNA fossils”, as they preserve sequences of ancient mtDNA and thus carry unique information about ancestral populations. Our comparison of a NUMT sequence shared by humans, chimpanzees, and gorillas with their mtDNAs implies that, around the time of divergence between humans and chimpanzees, our evolutionary history involved the interbreeding of individuals whose mtDNA had diverged as much as ~4.5 Myr prior. This large divergence suggests a distant interspecies hybridization. Additionally, analysis of two other NUMTs suggests that such events occurred repeatedly. Our findings suggest a complex pattern of speciation in primate human ancestors and provide a potential explanation for the mosaic nature of fossil morphology found at the emergence of the hominin lineage.

Gene expression levels can be subject to selection. We hypothesized that the age of gene origin is associated with expression constraints, given that it affects the level of gene integration into the functional cellular environment. By studying the genetic variation affecting gene expression levels (cis expression quantitative trait loci [cis-eQTLs]) and protein levels (cis protein QTLs [cis-pQTLs]), we determined that young, primate-specific genes are enriched in cis-eQTLs and cis-pQTLs. Compared to cis-eQTLs of old genes originating before the zebrafish divergence, cis-eQTLs of young genes have a higher effect size, are located closer to the transcription start site, are more significant, and tend to influence genes in multiple tissues and populations. These results suggest that the expression constraint of each gene increases throughout its lifespan. We also detected a positive correlation between expression constraints (approximated by cis-eQTL properties) and coding constraints (approximated by Ka/Ks) and observed that this correlation might be driven by gene age. To uncover factors associated with the increase in gene-age-related expression constraints, we demonstrated that gene connectivity, gene involvement in complex regulatory networks, gene haploinsufficiency, and the strength of posttranscriptional regulation increase with gene age. We also observed an increase in heritability of gene expression levels with age, implying a reduction of the environmental component. In summary, we show that gene age shapes key gene properties during evolution and is therefore an important component of genome function.

Somatic mtDNA mutations and deletions in particular are known to clonally expand within cells, eventually reaching detrimental intracellular concentrations. The possibility that clonal expansion is a slow process taking a lifetime had prompted an idea that founder mutations of mutant clones that cause mitochondrial dysfunction in the aged tissue might have originated early in life. If, conversely, expansion was fast, founder mutations should predominantly originate later in life. This distinction is important: indeed, from which mutations should we protect ourselves – those of early development/childhood or those happening at old age? Recently, high-resolution data describing the distribution of mtDNA deletions have been obtained using a novel, highly efficient method (Taylor et al., 2014). These data have been interpreted as supporting predominantly early origin of founder mutations. Re-analysis of the data implies that the data actually better fit mostly late origin of founders, although more research is clearly needed to resolve the controversy.

Background Mitochondria is a powerhouse of all eukaryotic cells that have its own circular DNA (mtDNA) encoding various RNAs and proteins. Somatic perturbations of mtDNA are accumulating with age thus it is of great importance to uncover the main sources of mtDNA instability. Recent analyses demonstrated that somatic mtDNA deletions depend on imperfect repeats of various nature between distant mtDNA segments. However, till now there are no comprehensive databases annotating all types of imperfect repeats in numerous species with sequenced complete mitochondrial genome as well as there are no algorithms capable to call all types of imperfect repeats in circular mtDNA. Results We implemented naïve algorithm of pattern recognition by analogy to standard dot-plot construction procedures allowing us to find both perfect and imperfect repeats of four main types: direct, inverted, mirror and complementary. Our algorithm is adapted to specific characteristics of mtDNA such as circularity and an excess of short repeats - it calls imperfect repeats starting from the length of 10 b.p. We constructed interactive web available database ImtRDB depositing perfect and imperfect repeats positions in mtDNAs of more than 3500 Vertebrate species. Additional tools, such as visualization of repeats within a genome, comparison of repeat densities among different genomes and a possibility to download all results make this database useful for many biologists. Our first analyses of the database demonstrated that mtDNA imperfect repeats (i) are usually short; (ii) associated with unfolded DNA structures; (iii) four types of repeats positively correlate with each other forming two equivalent pairs: direct and mirror versus inverted and complementary, with identical nucleotide content and similar distribution between species; (iv) abundance of repeats is negatively associated with GC content; (v) dinucleotides GC versus CG are overrepresented on light chain of mtDNA covered by repeats. Conclusions ImtRDB is available at http://bioinfodbs.kantiana.ru/ImtRDB/. It is accompanied by the software calling all types of interspersed repeats with different level of degeneracy in circular DNA. This database and software can become a very useful tool in various areas of mitochondrial and chloroplast DNA research.

Eukaryotic cells contain dynamic mitochondrial filaments: they fuse and divide. Here we summarize data on the protein machinery driving mitochondrial dynamics in yeast and also discuss the factors that affect the fusion-fission balance. Fission is a general stress response of cells, and in the case of yeast this response appears to be prosurvival. At the same time, even under normal conditions yeast mitochondria undergo continuous cycles of fusion and fission. This seems to be a futile cycle and also expensive from the energy point of view. Why does it exist? Benefits might be the same as in the case of sexual reproduction. Indeed, mixing and separating of mitochondrial content allows mitochondrial DNA to segregate and recombine randomly, leading to high variation in the numbers of mutations per individual mitochondrion. This opens a possibility for effective purifying selection-elimination of mitochondria highly contaminated by deleterious mutations. The beneficial action presumes a mechanism for removal of defective mitochondria. We argue that selective mitochondrial autophagy or asymmetrical distribution of mitochondria during cell division could be at the core of such mechanism.

Fragments of mitochondrial DNA are known to get inserted into nuclear DNA to form NUMTs, i.e. nuclear pseudogenes of the mtDNA. The insertion of a NUMT is a rare event. Hundreds of pseudogenes have been cataloged in the human genome. NUMTs are, in essence, a special type of mutation with their own internal timer, which is synchronized with an established molecular clock, the mtDNA. Thus insertion of NUMTs can be timed with respect to evolution milestones such as the emergence of new species. We asked whether NUMTs were inserted uniformly over time or preferentially during certain periods of evolution, as implied by the “punctuated evolution” model. To our surprise, the NUMT insertion times do appear nonrandom with at least one cluster positioned at around 2.8 million years ago (Ma). Interestingly, 2.8 Ma closely corresponds to the time of emergence of the genus Homo, and to a well-documented period of major climate change ca. 2.9–2.5 Ma. It is tempting to hypothesize that the insertion of NUMTs is related to the speciation process. NUMTs could be either “riders”, i.e., their insertion could be facilitated by the overall higher genome rearrangement activity during speciation, or “drivers”, i.e. they may more readily get fixed in the population due to positive selection associated with speciation. If correct, the hypothesis would support the idea that evolution of our genus may have happened in a rapid, punctuated manner.

The mammalian mitochondrial genomes differ from the nuclear genomes by maternal inheritance, absence of recombination, and higher mutation rate. All these differences decrease the effective population size of mitochondrial genome and make it more susceptible to accumulation of slightly deleterious mutations. It was hypothesized that mitochondrial genes, especially in species with low effective population size, irreversibly degrade leading to decrease of organismal fitness and even to extinction of species through the mutational meltdown. To interrogate this hypothesis, we compared the purifying selections acting on the representative set of mitochondrial (potentially degrading) and nuclear (potentially not degrading) protein-coding genes in species with different effective population size. For 21 mammalian species, we calculated the ratios of accumulation of slightly deleterious mutations approximated by Kn/Ks separately for mitochondrial and nuclear genomes. The 75% of variation in Kn/Ks is explained by two independent variables: type of a genome (mitochondrial or nuclear) and effective population size of species approximated by generation time. First, we observed that purifying selection is more effective in mitochondria than in the nucleus that implies strong evolutionary constraints of mitochondrial genome. Mitochondrial de novo nonsynonymous mutations have at least 5-fold more harmful effect when compared with nuclear. Second, Kn/Ks of mitochondrial and nuclear genomes is positively correlated with generation time of species, indicating relaxation of purifying selection with decrease of species-specific effective population size. Most importantly, the linear regression lines of mitochondrial and nuclear Kn/Ks's from generation times of species are parallel, indicating congruent relaxation of purifying selection in both genomes. Thus, our results reveal that the distribution of selection coefficients of de novo nonsynonymous mitochondrial mutations has a similar shape with the distribution of de novo nonsynonymous nuclear mutations, but its mean is five times smaller. The harmful effect of mitochondrial de novo nonsynonymous mutations triggers highly effective purifying selection, which maintains the fitness of the mammalian mitochondrial genome.

Accepted for publication September 11, 2014. Supported by the Ellison Medical Foundation award (K.K.) and Novartis grant no. 14B065 (K.P.). Disclosures: K.P. receives funding from Novartis. Address correspondence to Konstantin Khrapko, Ph.D., 134 Mugar, Department of Biology, Northeastern University, 360 Huntington Ave, Boston, MA 02115 E-mail: k.khrapko@neu.edu The article by Muller-Hocker et al in this issue of The American Journal of Pathology explores mitochondrial DNA (mtDNA) mutations of benign hyperplastic nodules of the hyperfunctional parathyroid gland. The study compares mtDNA mutations in cells differing by the extent of mitochondrial proliferation (ie, clear chief cells, pre-oxyphil, and oxyphil cells). Mitochondrial proliferation in these cells is thought to be a compensatory effect caused by a defect in respiratory chain. The transition from clear chief to oxyphil cells is an example of the oncocytic transformation, which is also observed in other tissues. Muller-Hocker et al comparemtDNAmutations in normal cellswithmutations in cells defective in cytochrome c oxidase (COX). COX, a key enzyme of the respiratory chain, is partially encoded by the mtDNA. Of particular importance is their discovery of substantial load of COX-inactivating mtDNA mutations and their conclusion that oncocytic transformation in the parathyroid gland can be caused by COX mutations. This conclusion extends an existing view that oncocytic transformation is caused almost exclusively by mutations in another respiratory chain enzyme partially encoded by mtDNA [ie, NADH dehydrogenase (ND)]. Although appearing too specialized, the impact of the highlighted study is quite broad as oncocytic cells give rise to oncocytic tumors, which are a rare example of tumors where involvement of mtDNA mutations is certain. This study may thus provide one of themissing tiles in the complicated puzzle depicting the role of mtDNA mutations in carcinogenesis.

The mtDNA mutator mouse lacks the proofreading capacity of the sole mtDNA polymerase, leading to accumulation of somatic mtDNA mutations, and a profound premature aging phenotype including elevated oxidative stress and apoptosis, and reduced mitochondrial function. We have previously reported that endurance exercise alleviates the aging phenotype in the mutator mice, reduces oxidative stress, and enhances mitochondrial biogenesis. Here we summarize our findings, with the emphasis on the central role of p53 in these adaptations. We demonstrate that mtDNA in sedentary and exercised PolG mice carry similar amounts of mutations in muscle, but in addition to that sedentary mice have more non-mutational damage, which is mitigated by exercise. It follows therefore that the profound alleviation of the mtDNA mutator phenotype in muscle by exercise may not require a reduction in mtDNA mutational load, but rather a decrease of mtDNA damage and/or oxidative stress. We further hypothesize that the observed 'alleviation without a reduction of mutational load' implies that the oxidative stress in PolG muscle is maintained, at least in part, by the 'malicious cycle', a hypothetical positive feedback potentially driven by the 'transcriptional mutagenesis', that is the conversion of chemically modified nucleotides into mutant RNA bases by the mitochondrial RNA polymerase.

Basal cell carcinoma (BCC) of the skin is the most common malignant neoplasm in humans. BCC is primarily driven by the Sonic Hedgehog (Hh) pathway. However, its phenotypic variation remains unexplained. Our genetic profiling of 293 BCCs found the highest mutation rate in cancer (65 mutations/Mb). Eighty-five percent of the BCCs harbored mutations in Hh pathway genes (PTCH1, 73% or SMO, 20% (P = 6.6 × 10(-8)) and SUFU, 8%) and in TP53 (61%). However, 85% of the BCCs also harbored additional driver mutations in other cancer-related genes. We observed recurrent mutations in MYCN (30%), PPP6C (15%), STK19 (10%), LATS1 (8%), ERBB2 (4%), PIK3CA (2%), and NRAS, KRAS or HRAS (2%), and loss-of-function and deleterious missense mutations were present in PTPN14 (23%), RB1 (8%) and FBXW7 (5%). Consistent with the mutational profiles, N-Myc and Hippo-YAP pathway target genes were upregulated. Functional analysis of the mutations in MYCN, PTPN14 and LATS1 suggested their potential relevance in BCC tumorigenesis.

APOBEC3A and APOBEC3B, cytidine deaminases of the APOBEC family, are among the main factors causing mutations in human cancers. APOBEC deaminates cytosines in single-stranded DNA (ssDNA). A fraction of the APOBEC-induced mutations occur as clusters (“kataegis”) in single-stranded DNA produced during repair of double-stranded breaks (DSBs). However, the properties of the remaining 87% of nonclustered APOBEC-induced mutations, the source and the genomic distribution of the ssDNA where they occur, are largely unknown. By analyzing genomic and exomic cancer databases, we show that >33% of dispersed APOBEC-induced mutations occur on the lagging strand during DNA replication, thus unraveling the major source of ssDNA targeted by APOBEC in cancer. Although methylated cytosine is generally more mutation-prone than nonmethylated cytosine, we report that methylation reduces the rate of APOBEC-induced mutations by a factor of roughly two. Finally, we show that in cancers with extensive APOBEC-induced mutagenesis, there is almost no increase in mutation rates in late replicating regions (contrary to other cancers). Because late-replicating regions are depleted in exons, this results in a 1.3-fold higher fraction of mutations residing within exons in such cancers. This study provides novel insight into the APOBEC-induced mutagenesis and describes the peculiarity of the mutational processes in cancers with the signature of APOBEC-induced mutations.

(2002). Body size of Holopedium gibberum under the influence of fish predation in small subarctic lakes. SIL Proceedings, 1922-2010: Vol. 28, No. 1, pp. 204-209.

Trisomy 21 is the most frequent genetic cause of cognitive impairment. To assess the perturbations of gene expression in trisomy 21, and to eliminate the noise of genomic variability, we studied the transcriptome of fetal fibroblasts from a pair of monozygotic twins discordant for trisomy 21. Here we show that the differential expression between the twins is organized in domains along all chromosomes that are either upregulated or downregulated. These gene expression dysregulation domains (GEDDs) can be defined by the expression level of their gene content, and are well conserved in induced pluripotent stem cells derived from the twins' fibroblasts. Comparison of the transcriptome of the Ts65Dn mouse model of Down's syndrome and normal littermate mouse fibroblasts also showed GEDDs along the mouse chromosomes that were syntenic in human. The GEDDs correlate with the lamina-associated (LADs) and replication domains of mammalian cells. The overall position of LADs was not altered in trisomic cells; however, the H3K4me3 profile of the trisomic fibroblasts was modified and accurately followed the GEDD pattern. These results indicate that the nuclear compartments of trisomic cells undergo modifications of the chromatin environment influencing the overall transcriptome, and that GEDDs may therefore contribute to some trisomy 21 phenotypes.

Congenital heart defect (CHD) occurs in 40% of Down syndrome (DS) cases. While carrying three copies of chromosome 21 increases the risk for CHD, trisomy 21 itself is not sufficient to cause CHD. Thus, additional genetic variation and/or environmental factors could contribute to the CHD risk. Here we report genomic variations that in concert with trisomy 21, determine the risk for CHD in DS. This case-control GWAS includes 187 DS with CHD (AVSD = 69, ASD = 53, VSD = 65) as cases, and 151 DS without CHD as controls. Chromosome 21–specific association studies revealed rs2832616 and rs1943950 as CHD risk alleles (adjusted genotypic P-values

Perfect direct repeats and, in particular, the prominent 13 bp repeat, are thought to cause mitochondrial DNA (mtDNA) deletions, which have been associated with the aging process. Accordingly, individuals lacking the 13 bp repeat are highly prevalent among centenarians and overall number of perfect repeats in mammalian mitochondrial genomes negatively correlates with species' longevity. However, detailed examination of the distribution of mtDNA deletions challenges a special role of the 13 bp repeat in generating mtDNA deletions. Instead, deletions appear to depend on long and stable, albeit imperfect, duplexes between distant mtDNA segments. Furthermore, significant dissimilarities in breakpoint distributions suggest that multiple mechanisms are involved in creating mtDNA deletions.

Direct nucleotide repeats can facilitate deletions of segments of mitochondrial genome1, leading to a wide range of neuromuscular disorders1,2 as well as aging2,3 in humans. We hypothesized that the number of the direct perfect repeats in human mitochondrial genomes influences longevity through the formation of harmful mtDNA deletions in the somatic cells. The analysis of the complete mitochondrial genomes of 762 unrelated Japanese individuals4-6 reveals a negative correlation between the abundance of the direct perfect repeats and the expected longevity. This association is largely due to the disruption of the common repeat (8470,13447) by a point mutation 8473C which occurred at the origin of the D4a haplogroup characterized by extreme longevity in Japan7. Our results provide the first evidence for correlation between the number of nucleotide repeats and the lifespan on intraspecific level.

Background Mitochondrial tRNAs have been the subject of study for structural biologists interested in their secondary structure characteristics, evolutionary biologists have researched patterns of compensatory and structural evolution and medical studies have been directed towards understanding the basis of human disease. However, an up to date, manually curated database of mitochondrially encoded tRNAs from higher animals is currently not available. Description We obtained the complete mitochondrial sequence for 277 tetrapod species from GenBank and re-annotated all of the tRNAs based on a multiple alignment of each tRNA gene and secondary structure prediction made independently for each tRNA. The mitochondrial (mt) tRNA sequences and the secondary structure based multiple alignments are freely available as Supplemental Information online. Conclusion We compiled a manually curated database of mitochondrially encoded tRNAs from tetrapods with completely sequenced genomes. In the course of our work, we reannotated more than 10% of all tetrapod mt-tRNAs and subsequently predicted the secondary structures of 6060 mitochondrial tRNAs. This carefully constructed database can be utilized to enhance our knowledge in several different fields including the evolution of mt-tRNA secondary structure and prediction of pathogenic mt-tRNA mutations. In addition, researchers reporting novel mitochondrial genome sequences should check their tRNA gene annotations against our database to ensure a higher level of fidelity of their annotation.