Sequence diversity analyses of an improved rhesus macaque genome enhance its biomedical utility

The rhesus macaque (Macaca mulatta) is the most widely studied nonhuman primate (NHP) in biomedical research. We present an updated reference genome assembly (Mmul_10, contig N50 = 46 Mbp) that increases the sequence contiguity 120-fold and annotate it using 6.5 million full-length transcripts, thus improving our understanding of gene content, isoform diversity, and repeat organization. With the improved assembly of segmental duplications, we discovered new lineage-specific genes and expanded gene families that are potentially informative in studies of evolution and disease susceptibility. Whole-genome sequencing (WGS) data from 853 rhesus macaques identified 85.7 million single-nucleotide variants (SNVs) and 10.5 million indel variants, including potentially damaging variants in genes associated with human autism and developmental delay, providing a framework for developing noninvasive NHP models of human disease.

Analysis of conserved miRNAs in cynomolgus macaque genome using small RNA sequencing and homology searching

MicroRNAs (miRNAs) are important regulators that fine-tune diverse cellular activities. Cynomolgus macaques (Macaca fascicularis) are used extensively in biomedical and pharmaceutical research; however, substantially fewer miRNAs have been identified in this species than in humans. Consequently, we investigated conserved miRNA profiles in cynomolgus macaques by homology searching and small RNA sequencing. In total, 1,455 high-confidence miRNA gene loci were identified, 408 of which were also confirmed by RNA sequencing, including 73 new miRNA loci reported in cynomolgus macaques for the first time. Comparing miRNA expression with age, we found a positive correlation between sequence conservation and expression levels during miRNA evolution. Additionally, we found that the miRNA gene locations in cynomolgus macaque genome were very flexible. Most were embedded in intergenic spaces or introns and clustered together. Several miRNAs were found in certain gene locations, including 64 exon-resident miRNAs, six splice-site-overlapping miRNAs (SO-miRNAs), and two pairs of distinct mirror miRNAs. We also identified 78 miRNA clusters, 68 of which were conserved in the human genome, including 10 large miRNA clusters predicted to regulate diverse developmental and cellular processes in cynomolgus macaque. Thus, this study not only expands the number of identified miRNAs in cynomolgus macaques but also provides clues for future research on the differences in miRNA repertoire between macaques and humans.

Advancing the genetic utility of pre-clinical species through a high-quality assembly of the cynomolgus monkey (Macaca fascicularis) genome

The cynomolgus macaque is a non-human primate model, heavily used in biomedical research, but with outdated genomic resources. Here we have used the latest long-read sequencing technologies in order to assemble a fully phased, chromosome-level assembly for the cynomolgus macaque. We have built a hybrid assembly with PacBio, 10x Genomics, and HiC technologies, resulting in a diploid assembly that spans a length of 5.1 Gb with a total of 16,741 contigs (N50 of 0.86Mb) contained in 370 scaffolds (N50 of 138 Mb) positioned on 42 chromosomes (21 homologous pairs). This assembly is highly homologous to former assemblies and identifies novel inversions and provides higher confidence in the genetic architecture of the cynomolgus macaque genome. A demographic estimation is also able to capture the recent genetic bottleneck in the Mauritius population, from which the sequenced individual originates. We offer this resource as an enablement for genetic tools to be built around this important model for biomedical research.

Comparative Genomics Analysis Reveals High Levels of Differential Retrotransposition among Primates from the Hominidae and the Cercopithecidae Families

Mobile elements (MEs), making ∼50% of primate genomes, are known to be responsible for generating inter- and intra-species genomic variations and play important roles in genome evolution and gene function. Using a bioinformatics comparative genomics approach, we performed analyses of species-specific MEs (SS-MEs) in eight primate genomes from the families of Hominidae and Cercopithecidae, focusing on retrotransposons. We identified a total of 230,855 SS-MEs, with which we performed normalization based on evolutionary distances, and we also analyzed the most recent SS-MEs in these genomes. Comparative analysis of SS-MEs reveals striking differences in ME transposition among these primate genomes. Interesting highlights of our results include: 1) the baboon genome has the highest number of SS-MEs with a strong bias for SINEs, while the crab-eating macaque genome has a sustained extremely low transposition for all ME classes, suggesting the existence of a genome-wide mechanism suppressing ME transposition; 2) while SS-SINEs represent the dominant class in general, the orangutan genome stands out by having SS-LINEs as the dominant class; 3) the human genome stands out among the eight genomes by having the largest number of recent highly active ME subfamilies, suggesting a greater impact of ME transposition on its recent evolution; and 4) at least 33% of the SS-MEs locate to genic regions, including protein coding regions, presenting significant potentials for impacting gene function. Our study, as the first of its kind, demonstrates that mobile elements evolve quite differently among these primates, suggesting differential ME transposition as an important mechanism in primate evolution.

Chromosome-levelde novoassembly of the pig-tailed macaque genome using linked-read sequencing and HiC proximity scaffolding

Old world monkey species share over 93% genome homology with humans and develop many disease phenotypes similar to those of humans, making them highly valuable animal models for the study of numerous human diseases. However, the quality of genome assembly and annotation for old world monkeys including macaque species lags behind the human genome effort. To close this gap and enhance functional genomics approaches, we employed a combination ofde novolinked-read assembly and scaffolding using proximity ligation assay (HiC) to assemble the pig-tailed macaque (Macaca nemestrina) genome. This combinatorial method yielded large scaffolds at chromosome-level with a scaffold N50 of 127.5 Mb; the 23 largest scaffolds covered 90% of the entire genome. This assembly revealed large-scale rearrangements between pig-tailed macaque chromosomes 7,12, and13 and human chromosomes 2, 14, and 15.

Comparative genomics analysis reveals high levels of differential DNA transposition among primates

ABSTRACTMobile elements generated via DNA transposition constitute ∼50% of the primate genomes. As a result of past and ongoing activity, DNA transposition is responsible for generating inter- and intra-species genomic variations, and it plays important roles in shaping genome evolution and impacting gene function. While limited analysis of mobile elements has been performed in many primate genomes, a large-scale comparative genomic analysis examining the impact of DNA transposition on primate evolution is still missing.Using a bioinformatics comparative genomics approach, we performed analysis of species-specific mobile elements (SS-MEs) in eight primate genomes, which include human, chimpanzee, gorilla, orangutan, green monkey, crab-eating macaque, rhesus monkey, and baboon. These species have good representations for the top two primate families, Hominidae (great apes) and the Cercopithecidae (old world monkeys), for which draft genome sequences are available.Our analysis identified a total of 230,855 SS-MEs from the eight primate genomes, which collectively contribute to ∼82 Mbp genome sequences, ranging from 14 to 25 Mbp for individual genomes. Several new interesting observations were made based on these SS-MEs. First, the DNA transposition activity level reflected by the numbers of SS-MEs was shown to be drastically different across species with the highest (baboon genome) being more than 30 times higher than the lowest (crab-eating macaque genome). Second, the compositions of SS-MEs, as well as the top active ME subfamilies, also differ significantly across genomes. By the copy numbers of SS-MEs divided into major ME classes, SINE represents the dominant class in all genomes, but more so in the Cercopithecidae genomes than in the Hominidae genomes in general with the orangutan genome being the outliner of this trend by having LINE as the dominant class. While AluY represents the major SINE groups in the Hominidae genomes, AluYRa1 is the dominant SINE in the Cercopithecidae genomes. For LINEs, each Hominidae genome seems to have a unique most active L1 subfamily, but all Cercopithecidae genomes have L1RS2 as the most active LINEs. While genomes with a high number of SS-MEs all have one or more very active ME subfamilies, the crab-eating macaque genome, being the one with an extremely low level of DNA transposition, has no single ME class being very active, suggesting the existence of a genome-wide mechanism suppressing DNA transposition. Third, DNA transposons, despite being considered dead in primate genomes, were in fact shown to have a certain level of activity in all genomes examined with a total of ∼2,400 entries as SS-MEs. Among these SS-MEs, at least 23% locate to genic regions, including exons and regulatory elements, presenting significant potentials for their impact on gene function. Very interestingly, our data demonstrate that, among the eight primates included in this study, the human genome is shown to be the most actively evolving genome via DNA transposition as having the highest most recent activity of many ME subfamilies, notably the AluYa5/Yb8/Yb9, L1HS, and SVA-D subfamilies.Representing the first of its kind, our large-scale comparative genomics study has shown that mobile elements evolved quite differently among different groups and species of primates, indicating that differential DNA transposition has served as an important mechanism in primate evolution.