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

Comprehensive Identification of Mitochondrial Pseudogenes (NUMTs) in the Human Telomere-to-Telomere Reference Genome

Practices related to mitochondrial research have long been hindered by the presence of mitochondrial pseudogenes within the nuclear genome (NUMTs). Even though partially assembled human reference genomes like hg38 have included NUMTs compilation, the exhaustive NUMTs within the only complete reference genome (T2T-CHR13) remain unknown. Here, we comprehensively identified the fixed NUMTs within the reference genome using human pan-mitogenome (HPMT) from GeneBank. The inclusion of HPMT serves the purpose of establishing an authentic mitochondrial DNA (mtDNA) mutational spectrum for the identification of NUMTs, distinguishing it from the polymorphic variations found in NUMTs. Using HPMT, we identified approximately 10% of additional NUMTs in three human reference genomes under stricter thresholds. And we also observed an approximate 6% increase in NUMTs in T2T-CHR13 compared to hg38, including NUMTs on the short arms of chromosomes 13, 14, and 15 that were not assembled previously. Furthermore, alignments based on 20-mer from mtDNA suggested the presence of more mtDNA-like short segments within the nuclear genome, which should be avoided for short amplicon or cell free mtDNA detection. Finally, through the assay of transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) on cell lines before and after mtDNA elimination, we concluded that NUMTs have a minimal impact on bulk ATAC-seq, even though 16% of sequencing data originated from mtDNA.

Evolution of eukaryotes as a story of survival and growth of mitochondrial DNA over two billion years

Mitochondria's significance in human diseases and in functioning, health and death of eukaryotic cell has been acknowledged widely. Yet our perspective in cell biology and evolution remains nucleocentric. Mitochondrial DNA, by virtue of its omnipresence and species-level conservation, is used as a barcode in animal taxonomy. This article analyses various levels of containment structures that enclose mitochondrial DNA and advocates a fresh perspective wherein evolution of organic structures of the eukarya domain seem to support and facilitate survival and proliferation of mitochondrial DNA by splitting containers as they age and by directing them along two distinct pathways: destruction of containers with more mutant mitochondrial DNA and rejuvenation of containers with less mutant mitochondrial DNA.

Characterization of nuclear mitochondrial insertions in the whole genomes of primates

The transfer and integration of whole and partial mitochondrial genomes into the nuclear genomes of eukaryotes is an ongoing process that has facilitated the transfer of genes and contributed to the evolution of various cellular pathways. Many previous studies have explored the impact of these insertions, referred to as NumtS, but have focused primarily on older events that have become fixed and are therefore present in all individual genomes for a given species. We previously developed an approach to identify novel Numt polymorphisms from next-generation sequence data and applied it to thousands of human genomes. Here, we extend this analysis to 79 individuals of other great ape species including chimpanzee, bonobo, gorilla, orang-utan and also an old world monkey, macaque. We show that recent Numt insertions are prevalent in each species though at different apparent rates, with chimpanzees exhibiting a significant increase in both polymorphic and fixed Numt sequences as compared to other great apes. We further assessed positional effects in each species in terms of evolutionary time and rate of insertion and identified putative hotspots on chromosome 5 for Numt integration, providing insight into both recent polymorphic and older fixed reference NumtS in great apes in comparison to human events.

Towards understanding the evolutionary dynamics of mtDNA

Historically, mtDNA was considered a selectively neutral marker that was useful for estimating the population genetic history of the maternal lineage. Over time there has been an increasing appreciation of mtDNA and mitochondria in maintaining cellular and organismal health. Beyond energy production, mtDNA and mitochondria have critical cellular roles in signalling. Here we briefly review the structure of mtDNA and the role of the mitochondrion in energy production. We then discuss the predictions that can be obtained from quaternary structure modelling and focus on mitochondrial complex I. Complex I is the primary entry point for electrons into the electron transport system is the largest respiratory complex of the chain and produces about 40% of the proton flux used to synthesize ATP. A focus of the review is Drosophila's utility as a model organism to study the selective advantage of specific mutations. However, we note that the incorporation of insights from a multitude of systems is necessary to fully understand the range of roles that mtDNA has in organismal fitness. We speculate that dietary changes can illicit stress responses that influence the selective advantage of specific mtDNA mutations and cause spatial and temporal fluctuations in the frequencies of mutations. We conclude that developing our understanding of the roles mtDNA has in determining organismal fitness will enable increased evolutionary insight and propose we can no longer assume it is evolving as a strictly neutral marker without testing this hypothesis.

Natural variation and evolutionary dynamics of transposable elements in Brassica oleracea based on next-generation sequencing data

Brassica oleracea comprises various economically important vegetables and presents extremely diverse morphological variations. They provide a rich source of nutrition for human health and have been used as a model system for studying polyploidization. Transposable elements (TEs) account for nearly 40% of the B. oleracea genome and contribute greatly to genetic diversity and genome evolution. Although the proliferation of TEs has led to a large expansion of the B. oleracea genome, little is known about the population dynamics and evolutionary activity of TEs. A comprehensive mobilome profile of 45,737 TE loci was obtained from resequencing data from 121 diverse accessions across nine B. oleracea morphotypes. Approximately 70% (32,195) of the loci showed insertion polymorphisms between or within morphotypes. In particular, up to 1221 loci were differentially fixed among morphotypes. Further analysis revealed that the distribution of the population frequency of TE loci was highly variable across different TE superfamilies and families, implying a diverse expansion history during host genome evolution. These findings provide better insight into the evolutionary dynamics and genetic diversity of B. oleracea genomes and will potentially serve as a valuable resource for molecular markers and association studies between TE-based genomic variations and morphotype-specific phenotypic differentiation.

Exploring the causes of small effective population sizes in cyst nematodes using artificial Globodera pallida populations

The effective size of a population is the size of an ideal population which would undergo genetic drift at the same rate as the real population. The balance between selection and genetic drift depends on the effective population size ( Ne), rather than the real numbers of individuals in the population ( N). The objectives of the present study were to estimate Ne in the potato cyst nematode Globodera pallida and to explore the causes of a low Ne/ N ratio in cyst nematodes using artificial populations. Using a temporal analysis of 24 independent populations, the median Ne was 58 individuals (min Ne = 25 and max Ne = 228). Ne is commonly lower than N but in the case of cyst nematodes, the Ne/ N ratio was extremely low. Using artificial populations showed that this low ratio did not result from migration, selection and overlapping generations, but could be explain by the fact that G. pallida populations deviate in structure from the assumptions of the ideal population by having unequal sex ratios, high levels of inbreeding and a high variance in family sizes. The consequences of a low Ne, resulting in a strong intensity of genetic drift, could be important for their control because G. pallida populations will have a low capacity to adapt to changing environments.

NumtS colonization in mammalian genomes

The colonization of the nuclear genome by mitochondrial DNA is an ongoing process in eukaryotes and plays an important role in genomic variability. Notwithstanding the DNA sequence availability of about 100 complete eukaryotic genomes, up to now NumtS distribution has been fully reported for a small number of sequenced eukaryotic species. With the aim to clarify the time and way of NumtS evolution, we explored the genomic distribution of NumtS in 23 eukaryotic species using an intra/interspecies in silico approach based on a cross-species similarity search and deeply investigate the evolution of NumtS in mammals. The intra- and interspecies analysis underlined how some mitochondrial regions that populated nuclear genomes can be considered as hotspots. Considering the large amount of NumtS we found in platypus and opossum genomes, we hypothesized the occurrence of an earlier colonization that happened prior to the Prototherian/Therian mammal divergence, approximately 160-210 million years ago. These events are still detectable due to the species-specific dynamics that have affected these genomes. Phylogenetic analyses of NumtS derived from two different mitochondrial DNA loci allowed us to recognize the unusual NumtS evolution that acted differently on primate and non-primate species' genomes.

Numtogenesis as a mechanism for development of cancer

Transfer of genetic material from cytoplasmic organelles to the nucleus, an ongoing process, has implications in evolution, aging, and human pathologies such as cancer. The transferred mitochondrial DNA (mtDNA) fragments in the nuclear genome are called nuclear mtDNA or NUMTs. We have named the process numtogenesis, defining the term as the transfer of mtDNA into the nuclear genome, or, less specifically, the transfer of mitochondria or mitochondrial components into the nucleus. There is increasing evidence of the involvement of NUMTs in human biology and pathology. Although information pertaining to NUMTs and numtogenesis is sparse, the role of this aspect of mitochondrial biology to human cancers is apparent. In this review, we present available knowledge about the origin and mechanisms of numtogenesis, with special emphasis on the role of NUMTs in human malignancies. We describe studies undertaken in our laboratory and in others and discuss the influence of NUMTs in tumor initiation and progression and in survival of cancer patients. We describe suppressors of numtogenesis and evolutionary conserved mechanisms underlying numtogenesis in cancer. An understanding the emerging field of numtogenesis should allow comprehension of this process in various malignancies and other diseases and, more generally, in human health.

Extreme genomic erosion after recurrent demographic bottlenecks in the highly endangered Iberian lynx

BACKGROUND

Genomic studies of endangered species provide insights into their evolution and demographic history, reveal patterns of genomic erosion that might limit their viability, and offer tools for their effective conservation. The Iberian lynx (Lynx pardinus) is the most endangered felid and a unique example of a species on the brink of extinction.

RESULTS

We generate the first annotated draft of the Iberian lynx genome and carry out genome-based analyses of lynx demography, evolution, and population genetics. We identify a series of severe population bottlenecks in the history of the Iberian lynx that predate its known demographic decline during the 20th century and have greatly impacted its genome evolution. We observe drastically reduced rates of weak-to-strong substitutions associated with GC-biased gene conversion and increased rates of fixation of transposable elements. We also find multiple signatures of genetic erosion in the two remnant Iberian lynx populations, including a high frequency of potentially deleterious variants and substitutions, as well as the lowest genome-wide genetic diversity reported so far in any species.

CONCLUSIONS

The genomic features observed in the Iberian lynx genome may hamper short- and long-term viability through reduced fitness and adaptive potential. The knowledge and resources developed in this study will boost the research on felid evolution and conservation genomics and will benefit the ongoing conservation and management of this emblematic species.

A Cost-Effective Approach to Sequence Hundreds of Complete Mitochondrial Genomes

We present a cost-effective approach to sequence whole mitochondrial genomes for hundreds of individuals. Our approach uses small reaction volumes and unmodified (non-phosphorylated) barcoded adaptors to minimize reagent costs. We demonstrate our approach by sequencing 383 Fundulus sp. mitochondrial genomes (192 F. heteroclitus and 191 F. majalis). Prior to sequencing, we amplified the mitochondrial genomes using 4–5 custom-made, overlapping primer pairs, and sequencing was performed on an Illumina HiSeq 2500 platform. After removing low quality and short sequences, 2.9 million and 2.8 million reads were generated for F. heteroclitus and F. majalis respectively. Individual genomes were assembled for each species by mapping barcoded reads to a reference genome. For F. majalis, the reference genome was built de novo. On average, individual consensus sequences had high coverage: 61-fold for F. heteroclitus and 57-fold for F. majalis. The approach discussed in this paper is optimized for sequencing mitochondrial genomes on an Illumina platform. However, with the proper modifications, this approach could be easily applied to other small genomes and sequencing platforms.

Frequent somatic transfer of mitochondrial DNA into the nuclear genome of human cancer cells

Mitochondrial genomes are separated from the nuclear genome for most of the cell cycle by the nuclear double membrane, intervening cytoplasm, and the mitochondrial double membrane. Despite these physical barriers, we show that somatically acquired mitochondrial-nuclear genome fusion sequences are present in cancer cells. Most occur in conjunction with intranuclear genomic rearrangements, and the features of the fusion fragments indicate that nonhomologous end joining and/or replication-dependent DNA double-strand break repair are the dominant mechanisms involved. Remarkably, mitochondrial-nuclear genome fusions occur at a similar rate per base pair of DNA as interchromosomal nuclear rearrangements, indicating the presence of a high frequency of contact between mitochondrial and nuclear DNA in some somatic cells. Transmission of mitochondrial DNA to the nuclear genome occurs in neoplastically transformed cells, but we do not exclude the possibility that some mitochondrial-nuclear DNA fusions observed in cancer occurred years earlier in normal somatic cells.

The advent of genome-wide association studies for bacteria

Significant advances in sequencing technologies and genome-wide association studies (GWAS) have revealed substantial insight into the genetic architecture of human phenotypes. In recent years, the application of this approach in bacteria has begun to reveal the genetic basis of bacterial host preference, antibiotic resistance, and virulence. Here, we consider relevant differences between bacterial and human genome dynamics, apply GWAS to a global sample of Mycobacterium tuberculosis genomes to highlight the impacts of linkage disequilibrium, population stratification, and natural selection, and finally compare the traditional GWAS against phyC, a contrasting method of mapping genotype to phenotype based upon evolutionary convergence. We discuss strengths and weaknesses of both methods, and make suggestions for factors to be considered in future bacterial GWAS.

The genomic landscape of polymorphic human nuclear mitochondrial insertions

The transfer of mitochondrial genetic material into the nuclear genomes of eukaryotes is a well-established phenomenon that has been previously limited to the study of static reference genomes. The recent advancement of high throughput sequencing has enabled an expanded exploration into the diversity of polymorphic nuclear mitochondrial insertions (NumtS) within human populations. We have developed an approach to discover and genotype novel Numt insertions using whole genome, paired-end sequencing data. We have applied this method to a thousand individuals in 20 populations from the 1000 Genomes Project and other datasets and identified 141 new sites of Numt insertions, extending our current knowledge of existing NumtS by almost 20%. We find that recent Numt insertions are derived from throughout the mitochondrial genome, including the D-loop, and have integration biases that differ in some respects from previous studies on older, fixed NumtS in the reference genome. We determined the complete inserted sequence for a subset of these events and have identified a number of nearly full-length mitochondrial genome insertions into nuclear chromosomes. We further define their age and origin of insertion and present an analysis of their potential impact to ongoing studies of mitochondrial heteroplasmy and disease.

Transposable element dynamics among asymbiotic and ectomycorrhizal Amanita fungi

Transposable elements (TEs) are ubiquitous inhabitants of eukaryotic genomes and their proliferation and dispersal shape genome architectures and diversity. Nevertheless, TE dynamics are often explored for one species at a time and are rarely considered in ecological contexts. Recent work with plant pathogens suggests a link between symbiosis and TE abundance. The genomes of pathogenic fungi appear to house an increased abundance of TEs, and TEs are frequently associated with the genes involved in symbiosis. To investigate whether this pattern is general, and relevant to mutualistic plant-fungal symbioses, we sequenced the genomes of related asymbiotic (AS) and ectomycorrhizal (ECM) Amanita fungi. Using methods developed to interrogate both assembled and unassembled sequences, we characterized and quantified TEs across three AS and three ECM species, including the AS outgroup Volvariella volvacea. The ECM genomes are characterized by abundant numbers of TEs, an especially prominent feature of unassembled sequencing libraries. Increased TE activity in ECM species is also supported by phylogenetic analysis of the three most abundant TE superfamilies; phylogenies revealed many radiations within contemporary ECM species. However, the AS species Amanita thiersii also houses extensive amplifications of elements, highlighting the influence of additional evolutionary parameters on TE abundance. Our analyses provide further evidence for a link between symbiotic associations among plants and fungi, and increased TE activity, while highlighting the importance individual species’ natural histories may have in shaping genome architecture.

[Genetic diversity of microorganisms]

Microorganisms are important components of the biosphere in maintaining the ecological balance. With the development of molecular biology techniques, researches on the microbial genetic diversity have been developed from morphological and/or protein levels to molecular level. The development of high-throughout sequencing and metagenomics technology not only provide more abundant information and powerful evidence for understanding microbial diversities, but also have great significance for rational utilization and protection of biological resources. The advances in research on genetic diversity of microorganisms, such as separation and identification, population genetic structure, speciation, phylogeny, and evolution of microorganisms, were discussed in this paper.

The evolution of ultraconserved elements with different phylogenetic origins

Background

Ultraconserved elements of DNA have been identified in vertebrate and invertebrate genomes. These elements have been found to have diverse functions, including enhancer activities in developmental processes. The evolutionary origins and functional roles of these elements in cellular systems, however, have not yet been determined.

Results

Here, we identified a wide range of ultraconserved elements common to distant species, from primitive aquatic organisms to terrestrial species with complicated body systems, including some novel elements conserved in fruit fly and human. In addition to a well-known association with developmental genes, these DNA elements have a strong association with genes implicated in essential cell functions, such as epigenetic regulation, apoptosis, detoxification, innate immunity, and sensory reception. Interestingly, we observed that ultraconserved elements clustered by sequence similarity. Furthermore, species composition and flanking genes of clusters showed lineage-specific patterns. Ultraconserved elements are highly enriched with binding sites to developmental transcription factors regardless of how they cluster.

Conclusion

We identified large numbers of ultraconserved elements across distant species. Specific classes of these conserved elements seem to have been generated before the divergence of taxa and fixed during the process of evolution. Our findings indicate that these ultraconserved elements are not the exclusive property of higher modern eukaryotes, but rather transmitted from their metazoan ancestors.

Next-generation sequencing of human mitochondrial reference genomes uncovers high heteroplasmy frequency

We describe methods for rapid sequencing of the entire human mitochondrial genome (mtgenome), which involve long-range PCR for specific amplification of the mtgenome, pyrosequencing, quantitative mapping of sequence reads to identify sequence variants and heteroplasmy, as well as de novo sequence assembly. These methods have been used to study 40 publicly available HapMap samples of European (CEU) and African (YRI) ancestry to demonstrate a sequencing error rate <5.63×10(-4), nucleotide diversity of 1.6×10(-3) for CEU and 3.7×10(-3) for YRI, patterns of sequence variation consistent with earlier studies, but a higher rate of heteroplasmy varying between 10% and 50%. These results demonstrate that next-generation sequencing technologies allow interrogation of the mitochondrial genome in greater depth than previously possible which may be of value in biology and medicine.

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.

Mammalian NUMT insertion is non-random

It is well known that remnants of partial or whole copies of mitochondrial DNA, known as Nuclear MiTochondrial sequences (NUMTs), are found in nuclear genomes. Since whole genome sequences have become available, many bioinformatics studies have identified putative NUMTs and from those attempted to infer the factors involved in NUMT creation. These studies conclude that NUMTs represent randomly chosen regions of the mitochondrial genome. There is less consensus regarding the nuclear insertion sites of NUMTs - previous studies have discussed the possible role of retrotransposons, but some recent ones have reported no correlation or even anti-correlation between NUMT sites and retrotransposons. These studies have generally defined NUMT sites using BLAST with default parameters. We analyze a redefined set of human NUMTs, computed with a carefully considered protocol. We discover that the inferred insertion points of NUMTs have a strong tendency to have high-predicted DNA curvature, occur in experimentally defined open chromatin regions and often occur immediately adjacent to A + T oligomers. We also show clear evidence that their flanking regions are indeed rich in retrotransposons. Finally we show that parts of the mitochondrial genome D-loop are under-represented as a source of NUMTs in primate evolution.

[Population genetics of plant pathogens]

Comparing to natural ecosystems, the evolution of plant pathogens in agricultural ecosystems is generally faster due to high-density monocultures, large-scale application of agrochemicals, and international trade in agricultural products. Knowledge of the population genetics and evolutionary biology of plant pathogens is necessary to understand disease epidemiology, effectively breed and use resistant cultivars, and control plant diseases. In this article, we outlined the aims of population genetic studies in plant pathogens, discuss contributions of five evolutionary forces (i.e., mutation, gene flow, recombination, random genetic drift, and natural selection) to origin, maintenance, and distribution of genetic variation in time and space, and gave an overview of current research status in this field.

The reference human nuclear mitochondrial sequences compilation validated and implemented on the UCSC genome browser

Background

Eukaryotic nuclear genomes contain fragments of mitochondrial DNA called NumtS (Nuclear mitochondrial Sequences), whose mode and time of insertion, as well as their functional/structural role within the genome are debated issues. Insertion sites match with chromosomal breaks, revealing that micro-deletions usually occurring at non-homologous end joining loci become reduced in presence of NumtS. Some NumtS are involved in recombination events leading to fragment duplication. Moreover, NumtS are polymorphic, a feature that renders them candidates as population markers. Finally, they are a cause of contamination during human mtDNA sequencing, leading to the generation of false heteroplasmies.

Results

Here we present RHNumtS.2, the most exhaustive human NumtSome catalogue annotating 585 NumtS, 97% of which were here validated in a European individual and in HapMap samples. The NumtS complete dataset and related features have been made available at the UCSC Genome Browser. The produced sequences have been submitted to INSDC databases. The implementation of the RHNumtS.2 tracks within the UCSC Genome Browser has been carried out with the aim to facilitate browsing of the NumtS tracks to be exploited in a wide range of research applications.

Conclusions

We aimed at providing the scientific community with the most exhaustive overview on the human NumtSome, a resource whose aim is to support several research applications, such as studies concerning human structural variation, diversity, and disease, as well as the detection of false heteroplasmic mtDNA variants. Upon implementation of the NumtS tracks, the application of the BLAT program on the UCSC Genome Browser has now become an additional tool to check for heteroplasmic artefacts, supported by data available through the NumtS tracks.

Resetting the evolution of marine reptiles at the Triassic-Jurassic boundary

Ichthyosaurs were important marine predators in the Early Jurassic, and an abundant and diverse component of Mesozoic marine ecosystems. Despite their ecological importance, however, the Early Jurassic species represent a reduced remnant of their former significance in the Triassic. Ichthyosaurs passed through an evolutionary bottleneck at, or close to, the Triassic-Jurassic boundary, which reduced their diversity to as few as three or four lineages. Diversity bounced back to some extent in the aftermath of the end-Triassic mass extinction, but disparity remained at less than one-tenth of pre-extinction levels, and never recovered. The group remained at low diversity and disparity for its final 100 Myr. The end-Triassic mass extinction had a previously unsuspected profound effect in resetting the evolution of apex marine predators of the Mesozoic.

Nonadaptive processes in primate and human evolution

Evolutionary biology has tended to focus on adaptive evolution by positive selection as the primum mobile of evolutionary trajectories in species while underestimating the importance of nonadaptive evolutionary processes. In this review, I describe evidence that suggests that primate and human evolution has been strongly influenced by nonadaptive processes, particularly random genetic drift and mutation. This is evidenced by three fundamental effects: a relative relaxation of selective constraints (i.e., purifying selection), a relative increase in the fixation of slightly deleterious mutations, and a general reduction in the efficacy of positive selection. These effects are observed in protein-coding, regulatory regions, and in gene expression data, as well as in an augmentation of fixation of large-scale mutations, including duplicated genes, mobile genetic elements, and nuclear mitochondrial DNA. The evidence suggests a general population-level explanation such as a reduction in effective population size (N(e)). This would have tipped the balance between the evolutionary forces of natural selection and random genetic drift toward genetic drift for variants having small selective effects. After describing these proximate effects, I describe the potential consequences of these effects for primate and human evolution. For example, an increase in the fixation of slightly deleterious mutations could potentially have led to an increase in the fixation rate of compensatory mutations that act to suppress the effects of slightly deleterious substitutions. The potential consequences of compensatory evolution for the evolution of novel gene functions and in potentially confounding the detection of positively selected genes are explored. The consequences of the passive accumulation of large-scale genomic mutations by genetic drift are unclear, though evidence suggests that new gene copies as well as insertions of transposable elements into genes can potentially lead to adaptive phenotypes. Finally, because a decrease in selective constraint at the genetic level is expected to have effects at the morphological level, I review studies that compare rates of morphological change in various mammalian and island populations where N(e) is reduced. Furthermore, I discuss evidence that suggests that craniofacial morphology in the Homo lineage has shifted from an evolutionary rate constrained by purifying selection toward a neutral evolutionary rate.

Evolutionary dynamics of transposable elements in a small RNA world

Transposable elements (TEs) are selfish elements that cause harmful mutations, contribute to the structure of regulatory networks and shape the architecture of genomes. Natural selection against their harmful effects has long been considered the dominant force limiting their spread. It is now clear that a genome defense system of RNA-mediated silencing also plays a crucial role in limiting TE proliferation. A full understanding of TE evolutionary dynamics must consider how these forces jointly determine their proliferation within genomes. Here I consider these forces from two perspectives - dynamics within populations and evolutionary games within the germline. The analysis of TE dynamics from these two perspectives promises to provide new insight into their role in evolution.

Human brain evolution: harnessing the genomics (r)evolution to link genes, cognition, and behavior

The evolution of the human brain has resulted in numerous specialized features including higher cognitive processes such as language. Knowledge of whole-genome sequence and structural variation via high-throughput sequencing technology provides an unprecedented opportunity to view human evolution at high resolution. However, phenotype discovery is a critical component of these endeavors and the use of nontraditional model organisms will also be critical for piecing together a complete picture. Ultimately, the union of developmental studies of the brain with studies of unique phenotypes in a myriad of species will result in a more thorough model of the groundwork the human brain was built upon. Furthermore, these integrative approaches should provide important insights into human diseases.

Consanguineous Marriage and Human Evolution

Mate choice among early human groups and in many historical populations was subject to both demographic and social constraints, ensuring that most unions were between couples who had coinherited substantial proportions of their genomes from common ancestors. Even in populations in which close consanguineous marriage was proscribed, community endogamy would have been sufficient to ensure high levels of homozygosity. Consanguineous marriage remains the choice of an estimated 10.4% of the global population, although there has been an overall decline in its popularity, especially in developed countries. Recent studies have indicated that the shift from consanguineous marriage to panmixia has been accompanied by a reduction in homozygosity. The concomitant predicted decrease in incidence of both recessive single-gene disorders and more common adult-onset diseases will have a significant impact on the health of future generations.

The migration of mitochondrial DNA fragments to the nucleus affects the chronological aging process of Saccharomyces cerevisiae

Migration of fragmented mitochondrial DNA (mtDNA) to the nucleus has been shown to occur in multiple species including yeast, plants, and mammals. Several human diseases, including Pallister-Hall syndrome and mucolipidosis, can be initiated by mtDNA insertion mutagenesis of nuclear DNA. In yeast, we demonstrated that the rate of mtDNA fragments translocating to the nucleus increases during chronological aging. The yeast chronological lifespan (CLS) is determined by the survival of nondividing cell populations. Whereas yeast strains with elevated migration rates of mtDNA fragments to the nucleus showed accelerated chronological aging, strains with decreased mtDNA transfer rates to the nucleus exhibited an extended CLS. Although one of the most popular theories of aging is the free radical theory, migration of mtDNA fragments to the nucleus may also contribute to the chronological aging process by possibly increasing nuclear genomic instability in cells with advanced age.

Origins, evolution, and phenotypic impact of new genes

Ever since the pre-molecular era, the birth of new genes with novel functions has been considered to be a major contributor to adaptive evolutionary innovation. Here, I review the origin and evolution of new genes and their functions in eukaryotes, an area of research that has made rapid progress in the past decade thanks to the genomics revolution. Indeed, recent work has provided initial whole-genome views of the different types of new genes for a large number of different organisms. The array of mechanisms underlying the origin of new genes is compelling, extending way beyond the traditionally well-studied source of gene duplication. Thus, it was shown that novel genes also regularly arose from messenger RNAs of ancestral genes, protein-coding genes metamorphosed into new RNA genes, genomic parasites were co-opted as new genes, and that both protein and RNA genes were composed from scratch (i.e., from previously nonfunctional sequences). These mechanisms then also contributed to the formation of numerous novel chimeric gene structures. Detailed functional investigations uncovered different evolutionary pathways that led to the emergence of novel functions from these newly minted sequences and, with respect to animals, attributed a potentially important role to one specific tissue--the testis--in the process of gene birth. Remarkably, these studies also demonstrated that novel genes of the various types significantly impacted the evolution of cellular, physiological, morphological, behavioral, and reproductive phenotypic traits. Consequently, it is now firmly established that new genes have indeed been major contributors to the origin of adaptive evolutionary novelties.

The origins and impact of primate segmental duplications

Duplicated sequences are substrates for the emergence of new genes and are an important source of genetic instability associated with rare and common diseases. Analyses of primate genomes have shown an increase in the proportion of interspersed segmental duplications (SDs) within the genomes of humans and great apes. This contrasts with other mammalian genomes that seem to have their recently duplicated sequences organized in a tandem configuration. In this review, we focus on the mechanistic origin and impact of this difference with respect to evolution, genetic diversity and primate phenotype. Although many genomes will be sequenced in the future, resolution of this aspect of genomic architecture still requires high quality sequences and detailed analyses.

A comparative approach shows differences in patterns of numt insertion during hominoid evolution

Nuclear integrations of mitochondrial DNA (numts) are widespread among eukaryotes, although their prevalence differs greatly among taxa. Most knowledge of numt evolution comes from analyses of whole-genome sequences of single species or, more recently, from genomic comparisons across vast phylogenetic distances. Here we employ a comparative approach using human and chimpanzee genome sequence data to infer differences in the patterns and processes underlying numt integrations. We identified 66 numts that have integrated into the chimpanzee nuclear genome since the human-chimp divergence, which is significantly greater than the 37 numts observed in humans. By comparing these closely related species, we accurately reconstructed the preintegration target site sequence and deduced nucleotide changes associated with numt integration. From >100 species-specific numts, we quantified the frequency of small insertions, deletions, duplications, and instances of microhomology. Most human and chimpanzee numt integrations were accompanied by microhomology and short indels of the kind typically observed in the nonhomologous end-joining pathway of DNA double-strand break repair. Human-specific numts have integrated into regions with a significant deficit of transposable elements; however, the same was not seen in chimpanzees. From a separate data set, we also found evidence for an apparent increase in the rate of numt insertions in the last common ancestor of humans and the great apes using a polymerase chain reaction-based screen. Last, phylogenetic analyses indicate that mitochondrial-numt alignments must be at least 500 bp, and preferably >1 kb in length, to accurately reconstruct hominoid phylogeny and recover the correct point of numt insertion.

Repair-mediated duplication by capture of proximal chromosomal DNA has shaped vertebrate genome evolution

DNA double-strand breaks (DSBs) are a common form of cellular damage that can lead to cell death if not repaired promptly. Experimental systems have shown that DSB repair in eukaryotic cells is often imperfect and may result in the insertion of extra chromosomal DNA or the duplication of existing DNA at the breakpoint. These events are thought to be a source of genomic instability and human diseases, but it is unclear whether they have contributed significantly to genome evolution. Here we developed an innovative computational pipeline that takes advantage of the repetitive structure of genomes to detect repair-mediated duplication events (RDs) that occurred in the germline and created insertions of at least 50 bp of genomic DNA. Using this pipeline we identified over 1,000 probable RDs in the human genome. Of these, 824 were intra-chromosomal, closely linked duplications of up to 619 bp bearing the hallmarks of the synthesis-dependent strand-annealing repair pathway. This mechanism has duplicated hundreds of sequences predicted to be functional in the human genome, including exons, UTRs, intron splice sites and transcription factor binding sites. Dating of the duplication events using comparative genomics and experimental validation revealed that the mechanism has operated continuously but with decreasing intensity throughout primate evolution. The mechanism has produced species-specific duplications in all primate species surveyed and is contributing to genomic variation among humans. Finally, we show that RDs have also occurred, albeit at a lower frequency, in non-primate mammals and other vertebrates, indicating that this mechanism has been an important force shaping vertebrate genome evolution.

The repair of DNA double-strand breaks (DSBs) is essential for the maintenance of genome integrity. The mechanisms by which DSBs are repaired have been the subject of intense experimental investigations. It has emerged that several imperfect repair pathways exist in eukaryotes that have the potential to result in chromosomal alterations, including genomic duplications. However, it remains unclear to what extent these imperfect repair events have contributed to shaping genomes throughout evolution. Here we introduce an innovative computational approach that takes advantage of the repetitive nature of eukaryotic genomes to identify repair-mediated duplications (RD) that occurred during evolution. We discovered over one thousand RDs in the human genome, with two-thirds resulting from the capture of a chromosomal DNA segment located in close proximity to the presumed site of the DSB, giving rise to local genomic duplications. Comparative genomic analyses reveal that the mechanism has operated continuously, but with decreasing intensity during primate evolution, generating species-specific duplications in all primates surveyed and generating genomic variation among humans. Finally, we show that RDs have also occurred in non-primate mammals and other vertebrates, indicating that this is a previously under-appreciated force shaping vertebrate genomes.

Independent mammalian genome contractions following the KT boundary

Although it is generally accepted that major changes in the earth's history are significant drivers of phylogenetic diversification and extinction, such episodes may also have long-lasting effects on genomic architecture. Here we show that widespread reductions in genome size have occurred in multiple lineages of mammals subsequent to the Cretaceous-Tertiary (KT) boundary, whereas there is no evidence for such changes in other vertebrate, invertebrate, or land plant lineages. Although the mechanisms remain unclear, such shifts in mammalian genome evolution may be a consequence of an increase in the efficiency of selection against excess DNA resulting from post-KT population size expansions. Independent historical changes in genome architecture in diverse lineages raise a significant challenge to the idea that genome size is finely tuned to achieve adaptive phenotypic modifications and suggest that attempts to use phylogenetic analysis to infer ancestral genome sizes may be problematical.

Large bottleneck size in Cauliflower Mosaic Virus populations during host plant colonization

The effective size of populations (Ne) determines whether selection or genetic drift is the predominant force shaping their genetic structure and evolution. Despite their high mutation rate and rapid evolution, this parameter is poorly documented experimentally in viruses, particularly plant viruses. All available studies, however, have demonstrated the existence of huge within-host demographic fluctuations, drastically reducing Ne upon systemic invasion of different organs and tissues. Notably, extreme bottlenecks have been detected at the stage of systemic leaf colonization in all plant viral species investigated so far, sustaining the general idea that some unknown obstacle(s) imposes a barrier on the development of all plant viruses. This idea has important implications, as it appoints genetic drift as a constant major force in plant virus evolution. By co-inoculating several genetic variants of Cauliflower mosaic virus into a large number of replicate host plants, and by monitoring their relative frequency within the viral population over the course of the host systemic infection, only minute stochastic variations were detected. This allowed the estimation of the CaMV Ne during colonization of successive leaves at several hundreds of viral genomes, a value about 100-fold higher than that reported for any other plant virus investigated so far, and indicated the very limited role played by genetic drift during plant systemic infection by this virus. These results suggest that the barriers that generate bottlenecks in some plant virus species might well not exist, or can be surmounted by other viruses, implying that severe bottlenecks during host colonization do not necessarily apply to all plant-infecting viruses.

Demography and weak selection drive patterns of transposable element diversity in natural populations of Arabidopsis lyrata

Transposable elements (TEs) are the major component of most plant genomes, and characterizing their population dynamics is key to understanding plant genome complexity. Yet there have been few studies of TE population genetics in plant systems. To study the roles of selection, transposition, and demography in shaping TE population diversity, we generated a polymorphism dataset for six TE families in four populations of the flowering plant Arabidopsis lyrata. The TE data indicated significant differentiation among populations, and maximum likelihood procedures suggested weak selection. For strongly bottlenecked populations, the observed TE band-frequency spectra fit data simulated under neutral demographic models constructed from nucleotide polymorphism data. Overall, we propose that TEs are subjected to weak selection, the efficacy of which varies as a function of demographic factors. Thus, demographic effects could be a major factor driving distributions of TEs among plant lineages.

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.

Regulation of mammalian gene expression by retroelements and non-coding tandem repeats

Genomes of higher eukaryotes contain abundant non-coding repeated sequences whose overall biological impact is unclear. They comprise two categories. The first consists of retrotransposon-derived elements. These are three major families of retroelements (LINEs, SINEs and LTRs). SINEs are clustered in gene-rich regions and are found in promoters of genes while LINEs are concentrated in gene-poor regions and are depleted from promoters. The second class consists of non-coding tandem repeats (satellite DNAs and TTAGGG arrays), which are associated with mammalian centromeres, heterochromatin and telomeres. Terminal TTAGGG arrays are involved in telomere capping and satellite DNAs are located in heterochromatin, which is implicated in transcription silencing by gene repositioning (relocalization). It is unknown whether interstitial TTAGGG sequences, which are present in many vertebrates, have a function. Here, evidence will be presented that retroelements and TTAGGG arrays are involved in regulation of gene expression. Retroelements can provide binding sites for transcription factors and protect promoter CpG islands from repressive chromatin modifications, and may be also involved in nuclear compartmentalization of transcriptionally active and inactive domains. Interstitial telomere-like sequences can form dynamically maintained three-dimensional nuclear networks of transcriptionally inactive domains, which may be involved in transcription silencing like classic heterochromatin.

Human MHC architecture and evolution: implications for disease association studies

Major histocompatibility complex (MHC) variation is a key determinant of susceptibility and resistance to a large number of infectious, autoimmune and other diseases. Identification of the MHC variants conferring susceptibility to disease is problematic, due to high levels of variation and linkage disequilibrium. Recent cataloguing and analysis of variation over the complete MHC has facilitated localization of susceptibility loci for autoimmune diseases, and provided insight into the MHC's evolution. This review considers how the unusual genetic characteristics of the MHC impact on strategies to identify variants causing, or contributing to, disease phenotypes. It also considers the MHC in relation to novel mechanisms influencing gene function and regulation, such as epistasis, epigenetics and microRNAs. These developments, along with recent technological advances, shed light on genetic association in complex disease.

From human genetics and genomics to pharmacogenetics and pharmacogenomics: past lessons, future directions

A brief history of human genetics and genomics is provided, comparing recent progress in those fields with that in pharmacogenetics and pharmacogenomics, which are subsets of genetics and genomics, respectively. Sequencing of the entire human genome, the mapping of common haplotypes of single-nucleotide polymorphisms (SNPs), and cost-effective genotyping technologies leading to genome-wide association (GWA) studies - have combined convincingly in the past several years to demonstrate the requirements needed to separate true associations from the plethora of false positives. While research in human genetics has moved from monogenic to oligogenic to complex diseases, its pharmacogenetics branch has followed, usually a few years behind. The continuous discoveries, even today, of new surprises about our genome cause us to question reviews declaring that "personalized medicine is almost here" or that "individualized drug therapy will soon be a reality." As summarized herein, numerous reasons exist to show that an "unequivocal genotype" or even an "unequivocal phenotype" is virtually impossible to achieve in current limited-size studies of human populations. This problem (of insufficiently stringent criteria) leads to a decrease in statistical power and, consequently, equivocal interpretation of most genotype-phenotype association studies. It remains unclear whether personalized medicine or individualized drug therapy will ever be achievable by means of DNA testing alone.

Demographic histories of ERV-K in humans, chimpanzees and rhesus monkeys

We detected 19 complete endogenous retroviruses of the K family in the genome of rhesus monkey (Macaca mulatta; RhERV-K) and 12 full length elements in the genome of the common chimpanzee (Pan troglodytes; CERV-K). These sequences were compared with 55 human HERV-K and 20 CERV-K reported previously, producing a total data set of 106 full-length ERV-K genomes. Overall, 61% of the human elements compared to 21% of the chimpanzee and 47% of rhesus elements had estimated integration times less than 4.5 million years before present (MYBP), with an average integration times of 7.8 MYBP, 13.4 MYBP and 10.3 MYBP for HERV-K, CERV-K and RhERV-K, respectively. By excluding those ERV-K sequences generated by chromosomal duplication, we used 63 of the 106 elements to compare the population dynamics of ERV-K among species. This analysis indicated that both HERV-K and RhERV-K had similar demographic histories, including markedly smaller effective population sizes, compared to CERV-K. We propose that these differing ERV-K dynamics reflect underlying differences in the evolutionary ecology of the host species, such that host ecology and demography represent important determinants of ERV-K dynamics.