Forty million years of independent evolution: a mitochondrial gene and its corresponding nuclear pseudogene in primates

Geographic Genetic Structure of Alectoris chukar in Türkiye: Post-LGM-Induced Hybridization and Human-Mediated Contaminations

Türkiye is considered an important evolutionary area for Chukar partridge (Alectoris chukar), since it is both a potential ancestral area and a diversification center for the species. Using 2 mitochondrial (Cty-b and D-loop) and 13 polymorphic microsatellite markers, we investigated the geographic genetic structure of A. chukar populations to determine how past climatic fluctuations and human activities have shaped the gene pool of this species in Türkiye. Our results indicate, firstly, that only A. chukar of the genus Alectoris is present in Türkiye (Anatolia and Thrace), with no natural or artificial gene flow from congenerics. Secondly, the geographic genetic structure of the species in Türkiye has been shaped by topographic heterogeneity, Pleistocene climatic fluctuations, and artificial transport by humans. Third, there appears to be three genetic clusters: Thracian, Eastern, and Western. Fourth, the post-LGM demographic expansion of the Eastern and Western populations has formed a hybrid zone in Central Anatolia (~8 kyBP). Fifth, the rate of China clade-B contamination in Türkiye is about 8% in mtDNA and about 12% in nuDNA, with the Southeastern Anatolian population having the highest contamination. Sixth, the Thracian population was the most genetically distinct, with the lowest genetic diversity and highest level of inbreeding and no China clad-B contamination. These results can contribute to the conservation regarding A. chukar populations, especially the Thracian population.

Re-recognition of pseudogenes: From molecular to clinical applications

Pseudogenes were initially regarded as "nonfunctional" genomic elements that did not have protein-coding abilities due to several endogenous inactivating mutations. Although pseudogenes are widely expressed in prokaryotes and eukaryotes, for decades, they have been largely ignored and classified as gene "junk" or "relics". With the widespread availability of high-throughput sequencing analysis, especially omics technologies, knowledge concerning pseudogenes has substantially increased. Pseudogenes are evolutionarily conserved and derive primarily from a mutation or retrotransposon, conferring the pseudogene with a "gene repository" role to store and expand genetic information. In contrast to previous notions, pseudogenes have a variety of functions at the DNA, RNA and protein levels for broadly participating in gene regulation to influence the development and progression of certain diseases, especially cancer. Indeed, some pseudogenes have been proven to encode proteins, strongly contradicting their "trash" identification, and have been confirmed to have tissue-specific and disease subtype-specific expression, indicating their own value in disease diagnosis. Moreover, pseudogenes have been correlated with the life expectancy of patients and exhibit great potential for future use in disease treatment, suggesting that they are promising biomarkers and therapeutic targets for clinical applications. In this review, we summarize the natural properties, functions, disease involvement and clinical value of pseudogenes. Although our knowledge of pseudogenes remains nascent, this field deserves more attention and deeper exploration.

Archaic mitochondrial DNA inserts in modern day nuclear genomes

Background

Traces of interbreeding of Neanderthals and Denisovans with modern humans in the form of archaic DNA have been detected in the genomes of present-day human populations outside sub-Saharan Africa. Up to now, only nuclear archaic DNA has been detected in modern humans; we therefore attempted to identify archaic mitochondrial DNA (mtDNA) residing in modern human nuclear genomes as nuclear inserts of mitochondrial DNA (NUMTs).

Results

We analysed 221 high-coverage genomes from Oceania and Indonesia using an approach which identifies reads that map both to the nuclear and mitochondrial DNA. We then classified reads according to the source of the mtDNA, and found one NUMT of Denisovan mtDNA origin, present in 15 analysed genomes; analysis of the flanking region suggests that this insertion is more likely to have happened in a Denisovan individual and introgressed into modern humans with the Denisovan nuclear DNA, rather than in a descendant of a Denisovan female and a modern human male.

Conclusions

Here we present our pipeline for detecting introgressed NUMTs in next generation sequencing data that can be used on genomes sequenced in the future. Further discovery of such archaic NUMTs in modern humans can be used to detect interbreeding between archaic and modern humans and can reveal new insights into the nature of such interbreeding events.

Possibility of numt co-amplification from gigantic genome of Orthoptera: testing efficiency of standard PCR protocol in producing orthologous COI sequences

Mitochondrial DNA has been the preferential genome biodiversity studies. However, several factors contribute to its inadequacy. Numts constitute one of the main complications that prevent obtaining orthologous mitochondrial sequences. Orthoptera have been a model group in numt studies because of their huge genome size. In this study we aimed to; (i) test efficiency of standard PCR protocol in producing orthologous sequences of cytochrome C oxidase, (ii) study presence/absence of numts in several unstudied Orthoptera species, (iii) test if there is a threshold between the length of mtDNA targeted for amplification and possibility of encountering numts, and (iv) estimate reliability of the sequences in databases in light of these findings. For these aims we studied 38 species of Orthoptera representing different sublineages and genome sizes. DNA extracted from each sample was used to amplify five different fragments of COI region by standard PCR protocol. Sequenced PCR amplicons were checked for numt possibility by several different numt criteria. No sequences without numt signs were obtained for the first fragment. The number of samples with numt signs for the other four fragments differed between the suborders Ensifera and Caelifera. The percentage of samples with numt signs was higher in Caelifera than Ensifera for all fragments. The numt percentage considerably decreased for the longest two fragments. Numts are more prevalent in families with larger genome size. We arrived at the following conclusions: (i) numts are common in all members of Orthoptera, but, their prevalence differs among intra-lineages, especially more prevalent in Caelifera, (ii) there seems a correlation between numt rate and genome size, (iii) there is no threshold to avoid numt co-amplification, but, a 1,000 bp length may be a threshold for Ensifera, (iv) Folmer region of COI doesn't seem an appropriate marker for animal barcoding. Additionally, a phylogenetic tree produced from the numt sequences of fragment four detected in genus Anterastes suggested a paleonumt gained in generic ancestor a 3.5-4 times slower divergence rate for numt sequences.

First Report of a Mitochondrial Pseudogene in Agnathan Vertebrates (Cyclostomata: Petromyzontidae)

We report herein the characterization of a nuclear paralog of a fragment of the mitochondrial genome (a numt) in two closely related species of lampreys (Ichthyomyzon spp.). Although numts have been characterized in several vertebrate taxa, numts have yet to be reported for fishes in general. Given the phylogenetic position of lampreys relative to other vertebrates, the presence of numts within the lamprey genome is either evidence of an ancestral trait lost in other fishes but uniquely retained in agnathans and amniotes, or (more intriguingly) a product of the genome rearrangements these animals undergo during development.

Mitochondrial Epigenetics and Environmental Exposure

The rising toll of chronic and debilitating diseases brought about by the exposure to an ever expanding number of environmental pollutants and socio-economic factors is calling for action. The understanding of the molecular mechanisms behind the effects of environmental exposures can lead to the development of biomarkers that can support the public health fields of both early diagnosis and intervention to limit the burden of environmental diseases. The study of mitochondrial epigenetics carries high hopes to provide important biomarkers of exposure and disease. Mitochondria are in fact on the frontline of the cellular response to the environment. Modifications of the epigenetic factors regulating the mitochondrial activity are emerging as informative tools that can effectively report on the effects of the environment on the phenotype. Here, we will discuss the emerging field of mitochondrial epigenetics. This review describes the main epigenetic phenomena that modify the activity of the mitochondrial DNA including DNA methylation, long and short non-coding RNAs. We will discuss the unique pattern of mitochondrial DNA methylation, describe the challenges of correctly measuring it, and report on the existing studies that have analysed the correlation between environmental exposures and mitochondrial DNA methylation. Finally, we provide a brief account of the therapeutic approaches targeting mitochondria currently under consideration.

Sheep mitochondrial heteroplasmy arises from tandem motifs and unspecific PCR amplification

The mitochondrial DNA control region (D-loop) is a widely used molecular marker in evolutionary and phylogeographic research. However, the occurrence of heteroplasmy of the D-loop region within individuals has rarely been investigated. In this study, a total of 85 Chinese sheep were used to amplify a partial D-loop region, and 15 heteroplasmic animals (17.64%) were identified. A comparative analysis of the PCR amplification and cloning of the D-loop sequences from the heteroplasmic samples revealed most of the sequencing profile from the heteroplasmic regions started at the beginning of a 75-bp random repeat motif. In addition, a total of 22 nonsyngeneic sequences with a D-loop were found in 61 of the clones obtained from the 4 random heteroplasmic and 3 homozygote animals, and their genomic locations were compared for homology. In summary, the D-Loop sequencing profiles appear to be heteroplasmic and could arise from tandem repeat motifs and unspecific replication during PCR amplification; however, they are not likely due to the presence of multiple mitochondrial genomes within an individual.

Unbiased Mitoproteome Analyses Confirm Non-canonical RNA, Expanded Codon Translations

Proteomic MS/MS mass spectrometry detections are usually biased towards peptides cleaved by experimentally added digestion enzyme(s). Hence peptides resulting from spontaneous degradation and natural proteolysis usually remain undetected. Previous analyses of tryptic human proteome data (cleavage after K, R) detected non-canonical tryptic peptides translated according to tetra- and pentacodons (codons expanded by silent mono- and dinucleotides), and from transcripts systematically (a) deleting mono-, dinucleotides after trinucleotides (delRNAs), (b) exchanging nucleotides according to 23 bijective transformations. Nine symmetric and fourteen asymmetric nucleotide exchanges (X ↔ Y, e.g. A ↔ C; and X → Y → Z → X, e.g. A → C → G → A) produce swinger RNAs. Here unbiased reanalyses of these proteomic data detect preferentially non-canonical tryptic peptides despite assuming random cleavage. Unbiased analyses couldn't reconstruct experimental tryptic digestion if most detected non-canonical peptides were false positives. Detected non-tryptic non-canonical peptides map preferentially on corresponding, previously described non-canonical transcripts, as for tryptic non-canonical peptides. Hence unbiased analyses independently confirm previous trypsin-biased analyses that showed translations of del- and swinger RNA and expanded codons. Accounting for natural proteolysis completes trypsin-biased mitopeptidome analyses, independently confirms non-canonical transcriptions and translations.

Chimeric mitochondrial peptides from contiguous regular and swinger RNA

Previous mass spectrometry analyses described human mitochondrial peptides entirely translated from swinger RNAs, RNAs where polymerization systematically exchanged nucleotides. Exchanges follow one among 23 bijective transformation rules, nine symmetric exchanges (X ↔ Y, e.g. A ↔ C) and fourteen asymmetric exchanges (X → Y → Z → X, e.g. A → C → G → A), multiplying by 24 DNA's protein coding potential. Abrupt switches from regular to swinger polymerization produce chimeric RNAs. Here, human mitochondrial proteomic analyses assuming abrupt switches between regular and swinger transcriptions, detect chimeric peptides, encoded by part regular, part swinger RNA. Contiguous regular- and swinger-encoded residues within single peptides are stronger evidence for translation of swinger RNA than previously detected, entirely swinger-encoded peptides: regular parts are positive controls matched with contiguous swinger parts, increasing confidence in results. Chimeric peptides are 200 × rarer than swinger peptides (3/100,000 versus 6/1000). Among 186 peptides with > 8 residues for each regular and swinger parts, regular parts of eleven chimeric peptides correspond to six among the thirteen recognized, mitochondrial protein-coding genes. Chimeric peptides matching partly regular proteins are rarer and less expressed than chimeric peptides matching non-coding sequences, suggesting targeted degradation of misfolded proteins. Present results strengthen hypotheses that the short mitogenome encodes far more proteins than hitherto assumed. Entirely swinger-encoded proteins could exist.

•

Chimeric peptides are translated from contiguous regular and swinger RNA

•

They are 200x rarer than mitochondrial swinger peptides

•

Chimeric peptides integrated in regular mitochondrial proteins are downregulated

•

Contiguous regular parts are matched positive controls for swinger parts

•

The last point validates results beyond other statistical tests for robustness

Rampant nuclear insertion of mtDNA across diverse lineages within Orthoptera (Insecta)

Nuclear mitochondrial pseudogenes (numts) are non-functional fragments of mtDNA inserted into the nuclear genome. Numts are prevalent across eukaryotes and a positive correlation is known to exist between the number of numts and the genome size. Most numt surveys have relied on model organisms with fully sequenced nuclear genomes, but such analyses have limited utilities for making a generalization about the patterns of numt accumulation for any given clade. Among insects, the order Orthoptera is known to have the largest nuclear genome and it is also reported to include several species with a large number of numts. In this study, we use Orthoptera as a case study to document the diversity and abundance of numts by generating numts of three mitochondrial loci across 28 orthopteran families, representing the phylogenetic diversity of the order. We discover that numts are rampant in all lineages, but there is no discernable and consistent pattern of numt accumulation among different lineages. Likewise, we do not find any evidence that a certain mitochondrial gene is more prone to nuclear insertion than others. We also find that numt insertion must have occurred continuously and frequently throughout the diversification of Orthoptera. Although most numts are the result of recent nuclear insertion, we find evidence of very ancient numt insertion shared by highly divergent families dating back to the Jurassic period. Finally, we discuss several factors contributing to the extreme prevalence of numts in Orthoptera and highlight the importance of exploring the utility of numts in evolutionary studies.

Feedback-based, system-level properties of vertebrate-microbial interactions

Background

Improved characterization of infectious disease dynamics is required. To that end, three-dimensional (3D) data analysis of feedback-like processes may be considered.

Methods

To detect infectious disease data patterns, a systems biology (SB) and evolutionary biology (EB) approach was evaluated, which utilizes leukocyte data structures designed to diminish data variability and enhance discrimination. Using data collected from one avian and two mammalian (human and bovine) species infected with viral, parasite, or bacterial agents (both sensitive and resistant to antimicrobials), four data structures were explored: (i) counts or percentages of a single leukocyte type, such as lymphocytes, neutrophils, or macrophages (the classic approach), and three levels of the SB/EB approach, which assessed (ii) 2D, (iii) 3D, and (iv) multi-dimensional (rotating 3D) host-microbial interactions.

Results

In all studies, no classic data structure discriminated disease-positive (D+, or observations in which a microbe was isolated) from disease-negative (D–, or microbial-negative) groups: D+ and D– data distributions overlapped. In contrast, multi-dimensional analysis of indicators designed to possess desirable features, such as a single line of observations, displayed a continuous, circular data structure, whose abrupt inflections facilitated partitioning into subsets statistically significantly different from one another. In all studies, the 3D, SB/EB approach distinguished three (steady, positive, and negative) feedback phases, in which D– data characterized the steady state phase, and D+ data were found in the positive and negative phases. In humans, spatial patterns revealed false-negative observations and three malaria-positive data classes. In both humans and bovines, methicillin-resistant Staphylococcus aureus (MRSA) infections were discriminated from non-MRSA infections.

Conclusions

More information can be extracted, from the same data, provided that data are structured, their 3D relationships are considered, and well-conserved (feedback-like) functions are estimated. Patterns emerging from such structures may distinguish well-conserved from recently developed host-microbial interactions. Applications include diagnosis, error detection, and modeling.

Factors affecting the relative abundance of nuclear copies of mitochondrial DNA (numts) in hominoids

Although nuclear copies of mitochondrial DNA (numts) can originate from any portion of the mitochondrial genome, evidence from humans suggests that more variable parts of the mitochondrial genome, such as the mitochondrial control region (MCR), are under-represented in the nucleus. This apparent deficit might arise from the erosion of sequence identity in numts originating from rapidly evolving mitochondrial sequences. However, the extent to which mitochondrial sequence properties impacts the number of numts detected in genomic surveys has not been evaluated. In order to address this question, we: (1) conducted exhaustive BLAST searches of MCR numts in three hominoid genomes; (2) assessed numt prevalence across the four MCR sub-domains (HV1, CCD, HV2, and MCR(F)); (3) estimated their insertion rates in great apes (Hominoidea); and (4) examined the relationship between mitochondrial DNA variability and numt prevalence in sequences originating from MCR and coding regions of the mitochondrial genome. Results indicate a marked deficit of numts from HV2 and MCR(F) MCR sub-domains in all three species. These MCR sub-domains exhibited the highest proportion of variable sites and the lowest number of detected numts per mitochondrial site. Variation in MCR insertion rate between lineages was also observed with a pronounced burst in recent integrations within chimpanzees and orangutans. A deficit of numts from HV2/MCR(F) was observed regardless of age, whereas HV1 is under-represented only in older numts (>25 million years). Finally, more variable mitochondrial genes also exhibit a lower identity with nuclear copies and because of this, appear to be under-represented in human numt databases.

Numts help to reconstruct the demographic history of the ocellated lizard (Lacerta lepida) in a secondary contact zone

In northwestern Iberia, two largely allopatric Lacerta lepida mitochondrial lineages occur, L5 occurring to the south of Douro River and L3 to the north, with a zone of putative secondary contact in the region of the Douro River valley. Cytochrome b sequence chromatograms with polymorphisms at nucleotide sites diagnostic for the two lineages were detected in individuals in the region of the Douro River and further north within the range of L3. We show that these polymorphisms are caused by the presence of four different numts (I-IV) co-occurring with the L3 genome, together with low levels of heteroplasmy. Two of the numts (I and II) are similar to the mitochondrial genome of L5 but are quite divergent from the mitochondrial genome of L3 where they occur. We show that these numts are derived from the mitochondrial genome of L5 and were incorporated in L3 through hybridization at the time of secondary contact between the lineages. The additional incidence of these numts to the north of the putative contact zone is consistent with an earlier postglacial northward range expansion of L5, preceding that of L3. We show that genetic exchange between the lineages responsible for the origin of these numts in L3 after secondary contact occurred prior to, or coincident with, the northward expansion of L3. This study shows that, in the context of phylogeographic analysis, numts can provide evidence for past demographic events and can be useful tools for the reconstruction of complex evolutionary histories.

Conservation of plastid sequences in the plant nuclear genome for millions of years facilitates endosymbiotic evolution

The nuclear genome of eukaryotes contains large amounts of cytoplasmic organelle DNA (nuclear integrants of organelle DNA [norgs]). The recent sequencing of many mitochondrial and chloroplast genomes has enabled investigation of the potential role of norgs in endosymbiotic evolution. In this article, we describe a new polymerase chain reaction-based method that allows the identification and evolutionary study of recent and older norgs in a range of eukaryotes. We tested this method in the genus Nicotiana and obtained sequences from seven nuclear integrants of plastid DNA (nupts) totaling 25 kb in length. These nupts were estimated to have been transferred 0.033 to 5.81 million years ago. The spectrum of mutations present in the potential protein-coding sequences compared with the noncoding sequences of each nupt revealed that nupts evolve in a nuclear-specific manner and are under neutral evolution. Indels were more frequent in noncoding regions than in potential coding sequences of former chloroplastic DNA, most probably due to the presence of a higher number of homopolymeric sequences. Unexpectedly, some potential protein-coding sequences within the nupts still contained intact open reading frames for up to 5.81 million years. These results suggest that chloroplast genes transferred to the nucleus have in some cases several millions of years to acquire nuclear regulatory elements and become functional. The different factors influencing this time frame and the potential role of nupts in endosymbiotic gene transfer are discussed.

Mitochondrial genome sequences effectively reveal the phylogeny of Hylobates gibbons

Background

Uniquely among hominoids, gibbons exist as multiple geographically contiguous taxa exhibiting distinctive behavioral, morphological, and karyotypic characteristics. However, our understanding of the evolutionary relationships of the various gibbons, especially among Hylobates species, is still limited because previous studies used limited taxon sampling or short mitochondrial DNA (mtDNA) sequences. Here we use mtDNA genome sequences to reconstruct gibbon phylogenetic relationships and reveal the pattern and timing of divergence events in gibbon evolutionary history.

Methodology/Principal Findings

We sequenced the mitochondrial genomes of 51 individuals representing 11 species belonging to three genera (Hylobates, Nomascus and Symphalangus) using the high-throughput 454 sequencing system with the parallel tagged sequencing approach. Three phylogenetic analyses (maximum likelihood, Bayesian analysis and neighbor-joining) depicted the gibbon phylogenetic relationships congruently and with strong support values. Most notably, we recover a well-supported phylogeny of the Hylobates gibbons. The estimation of divergence times using Bayesian analysis with relaxed clock model suggests a much more rapid speciation process in Hylobates than in Nomascus.

Conclusions/Significance

Use of more than 15 kb sequences of the mitochondrial genome provided more informative and robust data than previous studies of short mitochondrial segments (e.g., control region or cytochrome b) as shown by the reliable reconstruction of divergence patterns among Hylobates gibbons. Moreover, molecular dating of the mitogenomic divergence times implied that biogeographic change during the last five million years may be a factor promoting the speciation of Sundaland animals, including Hylobates species.

Haplowebs as a graphical tool for delimiting species: a revival of Doyle's "field for recombination" approach and its application to the coral genus Pocillopora in Clipperton

BACKGROUND

Usual methods for inferring species boundaries from molecular sequence data rely either on gene trees or on population genetic analyses. Another way of delimiting species, based on a view of species as "fields for recombination" (FFRs) characterized by mutual allelic exclusivity, was suggested in 1995 by Doyle. Here we propose to use haplowebs (haplotype networks with additional connections between haplotypes found co-occurring in heterozygous individuals) to visualize and delineate single-locus FFRs (sl-FFRs). Furthermore, we introduce a method to quantify the reliability of putative species boundaries according to the number of independent markers that support them, and illustrate this approach with a case study of taxonomically difficult corals of the genus Pocillopora collected around Clipperton Island (far eastern Pacific).

RESULTS

One haploweb built from intron sequences of the ATP synthase β subunit gene revealed the presence of two sl-FFRs among our 74 coral samples, whereas a second one built from ITS sequences turned out to be composed of four sl-FFRs. As a third independent marker, we performed a combined analysis of two regions of the mitochondrial genome: since haplowebs are not suited to analyze non-recombining markers, individuals were sorted into four haplogroups according to their mitochondrial sequences. Among all possible bipartitions of our set of samples, thirteen were supported by at least one molecular dataset, none by two and only one by all three datasets: this congruent pattern obtained from independent nuclear and mitochondrial markers indicates that two species of Pocillopora are present in Clipperton.

CONCLUSIONS

Our approach builds on Doyle's method and extends it by introducing an intuitive, user-friendly graphical representation and by proposing a conceptual framework to analyze and quantify the congruence between sl-FFRs obtained from several independent markers. Like delineation methods based on population-level statistical approaches, our method can distinguish closely-related species that have not yet reached reciprocal monophyly at most or all of their loci; like tree-based approaches, it can yield meaningful conclusions using a number of independent markers as low as three. Future efforts will aim to develop programs that speed up the construction of haplowebs from FASTA sequence alignments and help perform the congruence analysis outlined in this article.

Transfers of mitochondrial DNA to the nuclear genome in the wasp Nasonia vitripennis

Many organisms carry nuclear sequences of mitochondrial origin (NUMTs). We have identified 76 NUMTs in 25 genomic locations in the jewel wasp Nasonia vitripennis. The total amount of NUMTs in Nasonia is 42 972 bp exceeding over four-fold that found in Tribolium castaneum and almost fifty-fold that found in Drosophila melanogaster, whereas Apis mellifera has an even larger number of NUMTs in its genome (over 230 kb). The Nasonia NUMTs were inserted by multiple independent events and frequently involved large fragments spanning multiple mitochondrial genes. Most of the NUMTs are recent transfers that occurred less than one million years ago after the speciation of N. vitripennis. Duplications and rearrangements in the nucleus have also occurred. Data suggest that NUMTs may be more common in hymenoptera than in other insect genomes.

Population polymorphism of nuclear mitochondrial DNA insertions reveals widespread diploidy associated with loss of heterozygosity in Debaryomyces hansenii

Debaryomyces hansenii, a yeast that participates in the elaboration of foodstuff, displays important genetic diversity. Our recent phylogenetic classification of this species led to the subdivision of the species into three distinct clades. D. hansenii harbors the highest number of nuclear mitochondrial DNA (NUMT) insertions known so far for hemiascomycetous yeasts. Here we assessed the intraspecific variability of the NUMTs in this species by testing their presence/absence first in 28 strains, with 21 loci previously detected in the completely sequenced strain CBS 767(T), and second in a larger panel of 77 strains, with 8 most informative loci. We were able for the first time to structure populations in D. hansenii, although we observed little NUMT insertion variability within the clades. We determined the chronology of the NUMT insertions, which turned out to correlate with the previously defined taxonomy and provided additional evidence that colonization of nuclear genomes by mitochondrial DNA is a dynamic process in yeast. In combination with flow cytometry experiments, the NUMT analysis revealed the existence of both haploid and diploid strains, the latter being heterozygous and resulting from at least four crosses among strains from the various clades. As in the diploid pathogen Candida albicans, to which D. hansenii is phylogenetically related, we observed a differential loss of heterozygosity in the diploid strains, which can explain some of the large genetic diversity found in D. hansenii over the years.

Nucleotide polymorphism in the wheat transcriptional activator Spa influences its pattern of expression and has pleiotropic effects on grain protein composition, dough viscoelasticity, and grain hardness

Storage protein activator (SPA) is a key regulator of the transcription of wheat (Triticum aestivum) grain storage protein genes and belongs to the Opaque2 transcription factor subfamily. We analyzed the sequence polymorphism of the three homoeologous Spa genes in hexaploid wheat. The level of polymorphism in these genes was high particularly in the promoter. The deduced protein sequences of each homoeolog and haplotype show greater than 93% identity. Two major haplotypes were studied for each Spa gene. The three Spa homoeologs have similar patterns of expression during grain development, with a peak in expression around 300 degree days after anthesis. On average, Spa-B is 10 and seven times more strongly expressed than Spa-A and Spa-D, respectively. The haplotypes are associated with significant quantitative differences in Spa expression, especially for Spa-A and Spa-D. Significant differences were found in the quantity of total grain nitrogen allocated to the gliadin protein fractions for the Spa-A haplotypes, whereas the synthesis of glutenins is not modified. Genetic association analysis between Spa and dough viscoelasticity revealed that Spa polymorphisms are associated with dough tenacity, extensibility, and strength. Except for Spa-A, these associations can be explained by differences in grain hardness. No association was found between Spa markers and the average single grain dry mass or grain protein concentration. These results demonstrate that in planta Spa is involved in the regulation of grain storage protein synthesis. The associations between Spa and dough viscoelasticity and grain hardness strongly suggest that Spa has complex pleiotropic functions during grain development.

Epistatic interactions between Opaque2 transcriptional activator and its target gene CyPPDK1 control kernel trait variation in maize

Association genetics is a powerful method to track gene polymorphisms responsible for phenotypic variation, since it takes advantage of existing collections and historical recombination to study the correlation between large genetic diversity and phenotypic variation. We used a collection of 375 maize (Zea mays ssp. mays) inbred lines representative of tropical, American, and European diversity, previously characterized for genome-wide neutral markers and population structure, to investigate the roles of two functionally related candidate genes, Opaque2 and CyPPDK1, on kernel quality traits. Opaque2 encodes a basic leucine zipper transcriptional activator specifically expressed during endosperm development that controls the transcription of many target genes, including CyPPDK1, which encodes a cytosolic pyruvate orthophosphate dikinase. Using statistical models that correct for population structure and individual kinship, Opaque2 polymorphism was found to be strongly associated with variation of the essential amino acid lysine. This effect could be due to the direct role of Opaque2 on either zein transcription, zeins being major storage proteins devoid of lysine, or lysine degradation through the activation of lysine ketoglutarate reductase. Moreover, we found that a polymorphism in the Opaque2 coding sequence and several polymorphisms in the CyPPDK1 promoter nonadditively interact to modify both lysine content and the protein-versus-starch balance, thus revealing the role in quantitative variation in plants of epistatic interactions between a transcriptional activator and one of its target genes.

Primate numts and reticulate evolution of capped and golden leaf monkeys (Primates: Colobinae)

A recent phylogenetic study of langurs and leaf monkeys of South Asia suggested a reticulate evolution of capped and golden leaf monkeys through ancient hybridization between Semnopithecus and Trachypithecus .To test this hybridization scenario, I analysed nuclear copies of the mitochondrial cytochrome b gene (numts) from capped,golden and Phayre's leaf monkeys. These numts were aligned with mitochondrial cytochrome b sequences of various species belonging to the genera Semnopithecus and Trachypithecus .In the phylogenetic tree derived from this alignment,the numts fell into three distinct clades (A,B and C) suggesting three independent integration events.Clade A was basal to Semnopithecus, and clades B and C were basal to Trachypithecus. Among the numts in clades A and C were sequences derived from species not represented in their respective sister mitochondrial groups.This unusual placement of certain numts is taken as additional support for the hybridization scenario. Based on the molecular dating of these integration events,hybridization is estimated to have occurred around 7.1 to 3.4 million years ago.Capped and golden leaf monkeys might have to be assigned to a new genus to reconcile their unique evolutionary history. Additionally, northeast India appears to be a 'hot spot' for lineages that might have evolved through reticulate evolution.

Determining genetic background in captive stocks of cynomolgus macaques (Macaca fascicularis)

BACKGROUND

As in other model organisms, genetic background in the non-human primates Macaca mulatta and Macaca fascicularis is an experimental variable that affects the response of other study variables. Genetic background in model organisms is manipulated by breeding schemes but is generally pre-determined by the source population used to found captive stocks. In M. fascicularis three such sources have been distinguished, however, these are not routinely taken into consideration when designing research.

METHODS

We exemplify a mitochondrial DNA (mtDNA)-based strategy to trace the maternal geographic origins of M. fascicularis animals of unspecified origins.

RESULTS

Macaca fascicularis of unspecified origins kept at primate research centers carry mtDNA haplotypes representing all three major genetic subdivisions.

CONCLUSIONS

We suggest that the genetic background of study animals could be better specified in the future using an mtDNA-based approach, which would enable informed selection of study animals and help reduce variation within and among studies.

Mitochondrial COII Introgression into the Nuclear Genome of Gorilla gorilla

Numts are nonfunctional mitochondrial sequences that have translocated into nuclear DNA, where they evolve independently from the original mitochondrial DNA (mtDNA) sequence. Numts can be unintentionally amplified in addition to authentic mtDNA, complicating both the analysis and interpretation of mtDNA-based studies. Amplification of numts creates particular issues for studies on the noncoding, hypervariable 1 mtDNA region of gorillas. We provide data on putative numt sequences of the coding mitochondrial gene cytochrome oxidase subunit II (COII). Via polymerase chain reaction (PCR) and cloning, we obtained COII sequences for gorilla, orangutan, and human high-quality DNA and also from a gorilla fecal DNA sample. Both gorilla and orangutan samples yielded putative numt sequences. Phylogenetically more anciently transferred numts were amplified with a greater incidence from the gorilla fecal DNA sample than from the high-quality gorilla sample. Data on phylogenetically more recently transferred numts are equivocal. We further demonstrate the need for additional investigations into the use of mtDNA markers for noninvasively collected samples from gorillas and other primates.

Organization and evolution of mitochondrial gene clusters in human

Currently, the spatial patterns of mitochondrial genes and how the genomic localization of (pseudo)genes originated from mitochondrial DNA remain largely unexplained. The aim of this study was to elucidate the organization of mitochondrial (pseudo)genes given their evolutionary origin. We used a keyword finding method and a bootstrapping method to estimate parameter values that represent the distribution pattern of mitochondrial genes in the nuclear genome. Almost half of mitochondrial genes showing physical clusters were located in the pericentromeric and subtelomeric regions of the chromosome. Most interestingly, the size of these clusters ranged from 0.085 to 3.2 Mb (average+/-SD 1.3+/-0.73 Mb), which coincides with the size of the evolutionary pocket, or the average size of evolutionary breakpoint regions. Our findings imply that the localization of mitochondrial genes in the human genome is determined independent of adaptation.

The RHNumtS compilation: features and bioinformatics approaches to locate and quantify Human NumtS

Background

To a greater or lesser extent, eukaryotic nuclear genomes contain fragments of their mitochondrial genome counterpart, deriving from the random insertion of damaged mtDNA fragments. NumtS (Nuclear mt Sequences) are not equally abundant in all species, and are redundant and polymorphic in terms of copy number. In population and clinical genetics, it is important to have a complete overview of NumtS quantity and location. Searching PubMed for NumtS or Mitochondrial pseudo-genes yields hundreds of papers reporting Human NumtS compilations produced by in silico or wet-lab approaches. A comparison of published compilations clearly shows significant discrepancies among data, due both to unwise application of Bioinformatics methods and to a not yet correctly assembled nuclear genome. To optimize quantification and location of NumtS, we produced a consensus compilation of Human NumtS by applying various bioinformatics approaches.

Results

Location and quantification of NumtS may be achieved by applying database similarity searching methods: we have applied various methods such as Blastn, MegaBlast and BLAT, changing both parameters and database; the results were compared, further analysed and checked against the already published compilations, thus producing the Reference Human Numt Sequences (RHNumtS) compilation. The resulting NumtS total 190.

Conclusion

The RHNumtS compilation represents a highly reliable reference basis, which may allow designing a lab protocol to test the actual existence of each NumtS. Here we report preliminary results based on PCR amplification and sequencing on 41 NumtS selected from RHNumtS among those with lower score. In parallel, we are currently designing the RHNumtS database structure for implementation in the HmtDB resource. In the future, the same database will host NumtS compilations from other organisms, but these will be generated only when the nuclear genome of a specific organism has reached a high-quality level of assembly.

Population bottlenecks as a potential major shaping force of human genome architecture

The modern synthetic view of human evolution proposes that the fixation of novel mutations is driven by the balance among selective advantage, selective disadvantage, and genetic drift. When considering the global architecture of the human genome, the same model can be applied to understanding the rapid acquisition and proliferation of exogenous DNA. To explore the evolutionary forces that might have morphed human genome architecture, we investigated the origin, composition, and functional potential of numts (nuclear mitochondrial pseudogenes), partial copies of the mitochondrial genome found abundantly in chromosomal DNA. Our data indicate that these elements are unlikely to be advantageous, since they possess no gross positional, transcriptional, or translational features that might indicate beneficial functionality subsequent to integration. Using sequence analysis and fossil dating, we also show a probable burst of integration of numts in the primate lineage that centers on the prosimian–anthropoid split, mimics closely the temporal distribution of Alu and processed pseudogene acquisition, and coincides with the major climatic change at the Paleocene–Eocene boundary. We therefore propose a model according to which the gross architecture and repeat distribution of the human genome can be largely accounted for by a population bottleneck early in the anthropoid lineage and subsequent effectively neutral fixation of repetitive DNA, rather than positive selection or unusual insertion pressures.

Throughout evolutionary history, fragments of the mitochondrial genome, known as numts (for nuclear mitochondrial sequences), have been inserted into the nuclear genome. These fragments are distinct from all other classes of repetitive DNA found in nuclear genomes, not least because they are incapable of mediating their own proliferation. Taking advantage of their unique evolutionary properties, we have used numts to improve our understanding of the architecture of the human genome with special emphasis on the mechanism of acquisition and retention of repeat sequences, which comprise the bulk of nuclear DNA. We find that numts are unlikely to have any evolutionary benefit driving their retention. Moreover, numts are not acquired randomly during evolutionary time. Instead, their rate of acquisition spikes dramatically around pronounced population bottlenecks, in a manner reminiscent of other repeat classes. Therefore, we propose that the primary driving force of repeat acquisition in the genome is not selection, but random genetic drift, whose force becomes pronounced during profound reductions of population size. Our findings support the theory of neutral evolution, according to which random genetic drift exerts an influence on the acquisition of DNA changes that far outweighs the power of positive selection.

Cryptic ecological diversification of a planktonic estuarine copepod, Acartia tonsa

The recent discovery of cryptic species in marine holoplankton, organisms that 'drift' in oceanic currents throughout their life cycle, contrasts with their potential for long-distance passive dispersal and presumably high gene flow. These observations suggest that holoplankton species are adapting to surprisingly small-scale oceanographic features and imply either limited dispersal or strong selection gradients. Acartia tonsa is a widespread and numerically dominant estuarine copepod containing deep mitochondrial lineages within and among populations along the northwestern Atlantic coast. In this study, we intensively investigated A. tonsa populations in Chesapeake Bay with the goals of testing species status for the deep lineages and testing for their association with environmental features over space and time. Phylogenetic analyses of DNA sequences from mitochondrial cytochrome c oxidase I (mtCOI) and the nuclear ribosomal internal transcribed spacer (nITS) resolved two concordant monophyletic clades. Deep divergence between the two clades (13.7% uncorrected sequence divergence for mtCOI and 32.2% for nITS) and genealogical concordance within sympatric populations strongly suggest that the two clades represent reproductively isolated cryptic species. Based on restriction fragment length polymorphisms of mtCOI, representatives from the two clades were found consistently associated with contrasting salinity regimes (oligohaline vs. meso-polyhaline) with an overlap between 2 and 12 PSU in samples from 1995 to 2005. Finding these patterns in one of the best-known estuarine copepods reinforces the conclusion that marine biodiversity is underestimated, not only in terms of species numbers, but also with respect to niche partitioning and the potential importance of ecological divergence in marine holoplankton.

Extensive mitochondrial DNA transfer in a rapidly evolving rodent has been mediated by independent insertion events and by duplications

Mitochondrial DNA translocations to the nucleus (numt pseudogenes) are pervasive among eukaryotes, but copy number within the nuclear genome varies widely among taxa. As an increasing number of genomes are sequenced in their entirety, the origins, transfer mechanisms and insertion sites of numts are slowly being characterized. We investigated mitochondrial transfers within a genetically diverse rodent lineage and here report 15 numts totaling 21.8 kb that are harbored within the nuclear genome of the vole Microtus rossiaemeridionalis. The 15 numts total 21.8 kb and range from 0.39 to over 3.0 kb in length. Phylogenetic analyses revealed that these numts resulted from three independent insertions to the nucleus, two of which were followed by subsequent nuclear duplication events. The dates of the two translocations that led to subsequent duplications were estimated at 1.97 and 1.19 MYA, which coincide with the origin and radiation of the genus Microtus. Numt sequence data from five Microtus species were used to estimate an average rate of nucleotide substitution as 2.6x10(-8) subs/site/yr. This substitution rate is higher than in many other mammals, but is concordant with the elevated rate of mtDNA substitution in this lineage. Our data suggest that numt translocation in Microtus is more extensive than in either Mus or in Rattus, consistent with the elevated rate of speciation, karyotypic rearrangement, and mitochondrial DNA evolution in Microtus.

Distribution of nuclear mitochondrial DNA in cattle nuclear genome

The nuclear mitochondrial pseudogenes (numts), originated from mtDNA insertions into the nuclear genome, have been detected to exist in many species. However, the distribution of numts in cattle nuclear genome yet has not been fully reported. By referring to the whole cattle mtDNA sequence and to the recently released cattle nuclear genome by Human Genome Sequencing Center (HGSC), 303 numts were identified by BLAST with 55 numts unmapped to cattle nuclear genome. Further analysis found that the size of the numts ranges from 37 to 1932 bp, and the homologous identity between numts and their corresponding mtDNA fragments varies from 73 to 98%. Furthermore, the identified cattle numts cover nearly all the mitochondrial genes including mtDNA control region, distributing over all the chromosomes with the exception of the chromosome 23 and Y chromosome (in the latter the sequence data are not available). In the discovered numts, 29 relatively complete mitochondrial genes, which were distributed in 72 numts, were detected. Undoubtedly, this research would provide some valuable information for successive research related to mitochondrial genes and the evolution of cattle.

Mitochondrial Introgressions into the Nuclear Genome of the Domestic Cat

Translocation of mtDNA into the nuclear genome, also referred to as numt, was first reported in the domestic cat (Felis catus) by Lopez et al. (1994). The Lopez-numt consisted of a translocation of 7.9 kbp of mtDNA that inserted into the domestic cat chromosome D2 around 1.8 million years ago. More than a decade later, the release of the domestic cat whole-genome shotgun sequences (1.9x coverage) provides the resource to obtain more comprehensive insight into the extent of mtDNA transfer over time in the domestic cat genome. MegaBLAST searches revealed that the cat genome harbors a wide variety of numts (298 320 bp), one-third of which likely correspond to the Lopez-numt tandem repeat, whereas the remaining numts are probably derived from multiple independent insertions, which in some cases were followed by segmental duplication after insertion in the nucleus. Numts were detected across most cat chromosomes, but the number of numts assigned to chromosomes is underestimated due to the relatively high number of numt sequences with insufficient flanking sequence to map. The catalog of cat numts provides a valuable resource for future studies in Felidae species, including its use as a tool to avoid numt contaminations that may confound population genetics and phylogenetic studies.

Molecular analyses of mitochondrial pseudogenes within the nuclear genome of arvicoline rodents

Nuclear sequences of mitochondrial origin (numts) are common among animals and plants. The mechanism(s) by which numts transfer from the mitochondrion to the nucleus is uncertain, but their insertions may be mediated in part by chromosomal repair mechanisms. If so, then lineages where chromosomal rearrangements are common should be good models for the study of numt evolution. Arvicoline rodents are known for their karyotypic plasticity and numt pseudogenes have been discovered in this group. Here, we characterize a 4 kb numt pseudogene in the arvicoline vole Microtus rossiaemeridionalis. This sequence is among the largest numts described for a mammal lacking a completely sequenced genome. It encompasses three protein-coding and six tRNA pseudogenes that span approximately 25% of the entire mammalian mitochondrial genome. It is bordered by a dinucleotide microsatellite repeat and contains four transposable elements within its sequence and flanking regions. To determine the phylogenetic distribution of this numt among the arvicolines, we characterized one of the mitochondrial pseudogenes (cytochrome b) in 21 additional arvicoline species. Average rates of nucleotide substitution in this arvicoline pseudogene are estimated as 2.3 x 10(-8) substitutions/per site/per year. Furthermore, we performed comparative analyses among all species to estimate the age of this mitochondrial transfer at nearly 4 MYA, predating the origin of most arvicolines.

Toothed whale monophyly reassessed by SINE insertion analysis: the absence of lineage sorting effects suggests a small population of a common ancestral species

Morphological data have indicated that toothed whales form a monophyletic group. However, research published in the last several years has made the issue of the monophyly or paraphyly of toothed whales a subject of debate. Our group previously characterized three independent loci in which SINE insertions were shared among dolphins and sperm whales, thus supporting the traditional, morphologically based hypothesis of toothed whale monophyly. Although in recent years a few additional molecular works proposed this topology, there is still skepticism over this monophyly from the view point of molecular systematics. When the phylogeny of rapidly radiated taxa is examined using the SINE method, it is important to consider the ascertainment bias that arises when choosing a particular taxon for SINE loci screening. To overcome this methodological problem specific to the SINE method, we examined all possible topologies among sperm whales, dolphins and baleen whales by extensively screening SINE loci from species of all three lineages. We characterized nine independent SINE loci from the genomes of sperm whales and dolphins, all of which cluster sperm whales and dolphins but exclude baleen whales. Furthermore, we characterized ten independent loci from baleen whales, all of which were amplified in a common ancestor of these whales. From these observations, we conclude that toothed whales form a monophyletic group and that no ancestral SINE polymorphisms hinder their phylogenetic assignment despite the short divergence times of the major lineages of extant whales during evolution. These results suggest that a small population of common ancestors of all toothed whales ultimately diverged into the lineages of sperm whales and dolphins.

Sauria SINEs: Novel short interspersed retroposable elements that are widespread in reptile genomes

SINEs are short interspersed retrotransposable elements that invade new genomic sites. Their retrotransposition depends on reverse transcriptase and endonuclease activities encoded by partner LINEs (long interspersed elements). Recent genomic research has demonstrated that retroposons account for at least 40% of the human genome. Hitherto, more than 30 families of SINEs have been characterized in mammalian genomes, comprising approximately 4600 extant species; the distribution and extent of SINEs in reptilian genomes, however, are poorly documented. With more than 7400 species of lizards and snakes, Squamata constitutes the largest and most diverse group of living reptiles. We have discovered and characterized a novel SINE family, Sauria SINEs, whose members are widely distributed among genomes of lizards, snakes, and tuataras. Sauria SINEs comprise a 5' tRNA-related region, a tRNA-unrelated region, and a 3' tail region (containing short tandem repeats) derived from LINEs. We distinguished eight Sauria SINE subfamilies in genomes of four major squamate lineages and investigated their evolutionary relationships. Our data illustrate the overall efficacy of Sauria SINEs as novel retrotransposable markers for elucidation of squamate evolutionary history. We show that all Sauria SINEs share an identical 3' sequence with Bov-B LINEs and propose that they utilize the enzymatic machinery of Bov-B LINEs for their own retrotransposition. This finding, along with the ubiquity of Bov-B LINEs previously demonstrated in squamate genomes, suggests that these LINEs have been an active partner of Sauria SINEs since this SINE family was generated more than 200 million years ago.