Copy number variants and fixed duplications among 198 rhesus macaques (Macaca mulatta).

The rhesus macaque is an abundant species of Old World monkeys and a valuable model organism for biomedical research due to its close phylogenetic relationship to humans. Copy number variation is one of the main sources of genomic diversity within and between species and a widely recognized cause of inter-individual differences in disease risk. However, copy number differences among rhesus macaques and between the human and macaque genomes, as well as the relevance of this diversity to research involving this nonhuman primate, remain understudied. Here we present a high-resolution map of sequence copy number for the rhesus macaque genome constructed from a dataset of 198 individuals. Our results show that about one-eighth of the rhesus macaque reference genome is composed of recently duplicated regions, either copy number variable regions or fixed duplications. Comparison with human genomic copy number maps based on previously published data shows that, despite overall similarities in the genome-wide distribution of these regions, there are specific differences at the chromosome level. Some of these create differences in the copy number profile between human disease genes and their rhesus macaque orthologs. Our results highlight the importance of addressing the number of copies of target genes in the design of experiments and cautions against human-centered assumptions in research conducted with model organisms. Overall, we present a genome-wide copy number map from a large sample of rhesus macaque individuals representing an important novel contribution concerning the evolution of copy number in primate genomes.

GSEA of mouse and human mitochondriomes reveals fatty acid oxidation in astrocytes.

The prevalent view in neuroenergetics is that glucose is the main brain fuel, with neurons being mostly oxidative and astrocytes glycolytic. Evidence supporting that astrocyte mitochondria are functional has been overlooked. Here we sought to determine what is unique about astrocyte mitochondria by performing unbiased statistical comparisons of the mitochondriome in astrocytes and neurons. Using MitoCarta, a compendium of mitochondrial proteins, together with transcriptomes of mouse neurons and astrocytes, we generated cell-specific databases of nuclear genes encoding for mitochondrion proteins, ranked according to relative expression. Standard and in-house Gene Set Enrichment Analyses (GSEA) of five mouse transcriptomes revealed that genes encoding for enzymes involved in fatty acid oxidation (FAO) and amino acid catabolism are consistently more expressed in astrocytes than in neurons. FAO and oxidative-metabolism-related genes are also up-regulated in human cortical astrocytes versus the whole cortex, and in adult astrocytes versus fetal astrocytes. We thus present the first evidence of FAO in human astrocytes. Further, as shown in vitro, FAO coexists with glycolysis in astrocytes and is inhibited by glutamate. Altogether, these analyses provide arguments against the glucose-centered view of energy metabolism in astrocytes and reveal mitochondria as specialized organelles in these cells.

SeDuS: segmental duplication simulator.

SUMMARY: SeDuS is the first flexible and user-friendly forward-in-time simulator of patterns of molecular evolution within segmental duplications undergoing interlocus gene conversion and crossover. SeDuS introduces known features of interlocus gene conversion such as biased directionality and dependence on local sequence identity. Additionally, it includes aspects such as different selective pressures acting upon copy number and flexible crossover distributions. A graphical user interface allows fast fine-tuning of relevant parameters and straightforward real-time analysis of the evolution of duplicates. AVAILABILITY AND IMPLEMENTATION: SeDuS is implemented in C++ and can be run via command line or through a graphical user interface developed using Qt C++. Source code and binary executables for Linux, OS X and Windows are freely available at www.biologiaevolutiva.org/sedus/. A tutorial with a detailed description of implementation, parameters and output files is available online. CONTACT: arcadi.navarro@upf.edu.

Parallel evolution of amphioxus and vertebrate small-scale gene duplications.

BACKGROUND: Amphioxus are non-vertebrate chordates characterized by a slow morphological and molecular evolution. They share the basic chordate body-plan and genome organization with vertebrates but lack their 2R whole-genome duplications and their developmental complexity. For these reasons, amphioxus are frequently used as an outgroup to study vertebrate genome evolution and Evo-Devo. Aside from whole-genome duplications, genes continuously duplicate on a smaller scale. Small-scale duplicated genes can be found in both amphioxus and vertebrate genomes, while only the vertebrate genomes have duplicated genes product of their 2R whole-genome duplications. Here, we explore the history of small-scale gene duplications in the amphioxus lineage and compare it to small- and large-scale gene duplication history in vertebrates. RESULTS: We present a study of the European amphioxus (Branchiostoma lanceolatum) gene duplications thanks to a new, high-quality genome reference. We find that, despite its overall slow molecular evolution, the amphioxus lineage has had a history of small-scale duplications similar to the one observed in vertebrates. We find parallel gene duplication profiles between amphioxus and vertebrates and conserved functional constraints in gene duplication. Moreover, amphioxus gene duplicates show levels of expression and patterns of functional specialization similar to the ones observed in vertebrate duplicated genes. We also find strong conservation of gene synteny between two distant amphioxus species, B. lanceolatum and B. floridae, with two major chromosomal rearrangements. CONCLUSIONS: In contrast to their slower molecular and morphological evolution, amphioxus' small-scale gene duplication history resembles that of the vertebrate lineage both in quantitative and in functional terms.

Effect of Collapsed Duplications on Diversity Estimates: What to Expect.

The study of segmental duplications (SDs) and copy-number variants (CNVs) is of great importance in the fields of genomics and evolution. However, SDs and CNVs are usually excluded from genome-wide scans for natural selection. Because of high identity between copies, SDs and CNVs that are not included in reference genomes are prone to be collapsed-that is, mistakenly aligned to the same region-when aligning sequence data from single individuals to the reference. Such collapsed duplications are additionally challenging because concerted evolution between duplications alters their site frequency spectrum and linkage disequilibrium patterns. To investigate the potential effect of collapsed duplications upon natural selection scans we obtained expectations for four summary statistics from simulations of duplications evolving under a range of interlocus gene conversion and crossover rates. We confirm that summary statistics traditionally used to detect the action of natural selection on DNA sequences cannot be applied to SDs and CNVs since in some cases values for known duplications mimic selective signatures. As a proof of concept of the pervasiveness of collapsed duplications, we analyzed data from the 1,000 Genomes Project. We find that, within regions identified as variable in copy number, diversity between individuals with the duplication is consistently higher than between individuals without the duplication. Furthermore, the frequency of single nucleotide variants (SNVs) deviating from Hardy-Weinberg Equilibrium is higher in individuals with the duplication, which strongly suggests that higher diversity is a consequence of collapsed duplications and incorrect evaluation of SNVs within these CNV regions.

RNA editing independently occurs at three mir-376a-1 sites and may compromise the stability of the microRNA hairpin.

RNA editing is being recognized as an important post-transcriptional mechanism that may have crucial roles in introducing genetic variation and phenotypic diversity. Despite microRNA editing recurrence, defining its biological relevance is still under extended debate. To better understand microRNA editing function and regulation we performed an exhaustive characterization of the A-to-I site-specific patterns in mir-376a-1, a mammalian microRNA which RNA editing is involved in the regulation of development and in disease. Thorough an integrative approach based on high-throughput small RNA sequencing, Sanger sequencing and computer simulations we explored mir-376a-1 editing in samples from various individuals and primate species including human placenta and macaque, gorilla, chimpanzee and human brain cortex. We observed that mir-376a-1 editing is a common phenomenon in the mature and primary microRNA molecules and it is more frequently detected in brain than in placenta. Primary mir-376a-1 is edited at three positions, -1, +4 and +44. Editing frequency estimations and in silico simulations indicated that editing was not equally recurrent along the three mir-376a-1 sites, nevertheless no epistatic interactions among them were observed. Particularly, the +4 site, located in the seed region of the mature miR-376a-5p, reached the highest editing frequency in all samples. Secondary structure predictions revealed that the +4 position was the one that conferred the highest stability to the mir-376a-1 hairpin. We suggest that molecular stability might partially explain the editing recurrence observed in certain microRNAs and that editing events conferring new functional regulatory roles in particular tissues and species could have been conserved along evolution, as it might be the case of mir-376a-1 in primate brain cortex.

Interplay of interlocus gene conversion and crossover in segmental duplications under a neutral scenario.

Interlocus gene conversion is a major evolutionary force that drives the concerted evolution of duplicated genomic regions. Theoretical models successfully have addressed the effects of interlocus gene conversion and the importance of crossover in the evolutionary fate of gene families and duplications but have not considered complex recombination scenarios, such as the presence of hotspots. To study the interplay between interlocus gene conversion and crossover, we have developed a forward-time simulator that allows the exploration of a wide range of interlocus gene conversion rates under different crossover models. Using it, we have analyzed patterns of nucleotide variation and linkage disequilibrium within and between duplicate regions, focusing on a neutral scenario with constant population size and validating our results with the existing theoretical models. We show that the interaction of gene conversion and crossover is nontrivial and that the location of crossover junctions is a fundamental determinant of levels of variation and linkage disequilibrium in duplicated regions. We also show that if crossover activity between duplications is strong enough, recurrent interlocus gene conversion events can break linkage disequilibrium within duplicates. Given the complex nature of interlocus gene conversion and crossover, we provide a framework to explore their interplay to help increase knowledge on molecular evolution within segmental duplications under more complex scenarios, such as demographic changes or natural selection.