There are three main dietary groups in mammals: carnivores, omnivores, and herbivores. Currently, there is limited comparative genomics insight into the evolution of dietary specializations in mammals. Due to recent advances in sequencing technologies, we were able to perform in-depth whole genome analyses of representatives of these three dietary groups. We investigated the evolution of carnivory by comparing 18 representative genomes from across Mammalia with carnivorous, omnivorous, and herbivorous dietary specializations, focusing on Felidae (domestic cat, tiger, lion, cheetah, and leopard), Hominidae, and Bovidae genomes. We generated a new high-quality leopard genome assembly, as well as two wild Amur leopard whole genomes. In addition to a clear contraction in gene families for starch and sucrose metabolism, the carnivore genomes showed evidence of shared evolutionary adaptations in genes associated with diet, muscle strength, agility, and other traits responsible for successful hunting and meat consumption. Additionally, an analysis of highly conserved regions at the family level revealed molecular signatures of dietary adaptation in each of Felidae, Hominidae, and Bovidae. However, unlike carnivores, omnivores and herbivores showed fewer shared adaptive signatures, indicating that carnivores are under strong selective pressure related to diet. Finally, felids showed recent reductions in genetic diversity associated with decreased population sizes, which may be due to the inflexible nature of their strict diet, highlighting their vulnerability and critical conservation status. Our study provides a large-scale family level comparative genomic analysis to address genomic changes associated with dietary specialization. Our genomic analyses also provide useful resources for diet-related genetic and health research.

The cinereous vulture, Aegypius monachus, is the largest bird of prey and plays a key role in the ecosystem by removing carcasses, thus preventing the spread of diseases. Its feeding habits force it to cope with constant exposure to pathogens, making this species an interesting target for discovering functionally selected genetic variants. Furthermore, the presence of two independently evolved vulture groups, Old World and New World vultures, provides a natural experiment in which to investigate convergent evolution due to obligate scavenging. We sequenced the genome of a cinereous vulture, and mapped it to the bald eagle reference genome, a close relative with a divergence time of 18 million years. By comparing the cinereous vulture to other avian genomes, we find positively selected genetic variations in this species associated with respiration, likely linked to their ability of immune defense responses and gastric acid secretion, consistent with their ability to digest carcasses. Comparisons between the Old World and New World vulture groups suggest convergent gene evolution. We assemble the cinereous vulture blood transcriptome from a second individual, and annotate genes. Finally, we infer the demographic history of the cinereous vulture which shows marked fluctuations in effective population size during the late Pleistocene. We present the first genome and transcriptome analyses of the cinereous vulture compared to other avian genomes and transcriptomes, revealing genetic signatures of dietary and environmental adaptations accompanied by possible convergent evolution between the Old World and New World vultures.

Pakistan covers a key geographic area in human history, being both part of the Indus River region that acted as one of the cradles of civilization and as a link between Western Eurasia and Eastern Asia. This region is inhabited by a number of distinct ethnic groups, the largest being the Punjabi, Pathan (Pakhtuns), Sindhi, and Baloch. We analyzed the first ethnic male Pathan genome by sequencing it to 29.7-fold coverage using the Illumina HiSeq2000 platform. A total of 3.8 million single nucleotide variations (SNVs) and 0.5 million small indels were identified by comparing with the human reference genome. Among the SNVs, 129,441 were novel, and 10,315 nonsynonymous SNVs were found in 5,344 genes. SNVs were annotated for health consequences and high risk diseases, as well as possible influences on drug efficacy. We confirmed that the Pathan genome presented here is representative of this ethnic group by comparing it to a panel of Central Asians from the HGDP-CEPH panels typed for ~650 k SNPs. The mtDNA (H2) and Y haplogroup (L1) of this individual were also typical of his geographic region of origin. Finally, we reconstruct the demographic history by PSMC, which highlights a recent increase in effective population size compatible with admixture between European and Asian lineages expected in this geographic region. We present a whole-genome sequence and analyses of an ethnic Pathan from the north-west province of Pakistan. It is a useful resource to understand genetic variation and human migration across the whole Asian continent.

Background: MGISEQ-T7 is a new whole-genome sequencer developed by Complete Genomics and MGI utilizing DNA nanoball and combinatorial probe anchor synthesis technologies for generating short reads at a very large scale – up to 60 human genomes per day. However, it has not been objectively and systematically compared against Illumina short-read sequencers. Findings: By using the same KOREF sample, the Korean Reference Genome, we have compared seven sequencing platforms including BGISEQ-500, MGISEQ-T7, HiSeq2000, HiSeq2500, HiSeq4000, HiSeqX10, and NovaSeq6000. We measured sequencing quality by comparing sequencing statistics (base quality, duplication rate, and random error rate), mapping statistics (mapping rate, depth distribution, and %GC coverage), and variant statistics (transition/transversion ratio, dbSNP annotation rate, and concordance rate with SNP genotyping chip) across the seven sequencing platforms. We found that MGI platforms showed a higher concordance rate of SNP genotyping than HiSeq2000 and HiSeq4000. The similarity matrix of variant calls confirmed that the two MGI platforms have the most similar characteristics to the HiSeq2500 platform. Conclusions: Overall, MGI and Illumina sequencing platforms showed comparable levels of sequencing quality, uniformity of coverage, %GC coverage, and variant accuracy, thus we conclude that the MGI platforms can be used for a wide range of genomics research fields at approximately half the cost of the Illumina platforms.

Background: Unique among cnidarians, jellyfish have remarkable morphological and biochemical innovations that allow them to actively hunt in the water column. One of the first animals to become free-swimming, jellyfish employ pulsed jet propulsion and venomous tentacles to capture prey. Results: To understand these key innovations, we sequenced the genome of the giant Nomura’s jellyfish ( Nemopilema nomurai ), the transcriptomes of its bell and tentacles, and transcriptomes across tissues and developmental stages of the Sanderia malayensis jellyfish. Analyses of Nemopilema and other cnidarian genomes revealed adaptations associated with swimming, marked by codon bias in muscle contraction and expansion of neurotransmitter genes, along with expanded Myosin type II family and venom domains; possibly contributing to jellyfish mobility and active predation. We also identified gene family expansions of Wnt and posterior Hox genes, and discovered the important role of retinoic acid signaling in this ancient lineage of metazoans, which together may be related to the unique jellyfish body plan (medusa formation). Conclusions: Taken together, the jellyfish genome and transcriptomes genetically confirm their unique morphological and physiological traits that have combined to make these animals one of the world’s earliest and most successful multi-cellular predators.

Native to the mountains of East Asia, the long-tailed goral ( Naemorhedus caudatus ) is a vulnerable wild ungulate in the tribe Caprini. To understand key conservation issues related to fragmentation and subsequent endangerment of this montane species, we sequenced and analyzed the genome of a long-tailed goral to explore historical demography, contemporary levels of genetic diversity, and potential immune response. When compared to ten additional mammalian reference genomes, we identified 357 positively selected genes (PSGs) in the long-tailed goral and 364 PSGs in the Caprini lineage. Gene Ontology analyses showed statistical enrichment in biological processes related to immune function and in genes and pathways related to blood coagulation. We also identified low levels of heterozygosity (0.00114) in the long-tailed goral along with decreases in population size relative to other species of Caprini. Finally, we provide evidence for positive selection on the muscle development gene myostatin ( MSTN ) in the Caprini lineage, which may have increased the muscle development and climbing ability in the caprine common ancestor. Low effective population size and decreased heterozygosity of the long-tailed goral raise conservation concern about the effects of habitat fragmentation, over harvest, and inbreeding on this vulnerable species.