Comparative genomics of macaques and integrated insights into genetic variation and population history.

The crab-eating macaques ( Macaca fascicularis ) and rhesus macaques ( M. mulatta ) are widely studied nonhuman primates in biomedical and evolutionary research. Despite their significance, the current understanding of the complex genomic structure in macaques and the differences between species requires substantial improvement. Here, we present a complete genome assembly of a crab-eating macaque and 20 haplotype-resolved macaque assemblies to investigate the complex regions and major genomic differences between species. Segmental duplication in macaques is ∼42% lower, while centromeres are ∼3.7 times longer than those in humans. The characterization of ∼2 Mbp fixed genetic variants and ∼240 Mbp complex loci highlights potential associations with metabolic differences between the two macaque species (e.g., CYP2C76 and EHBP1L1 ). Additionally, hundreds of alternative splicing differences show post-transcriptional regulation divergence between these two species (e.g., PNPO ). We also characterize 91 large-scale genomic differences between macaques and humans at a single-base-pair resolution and highlight their impact on gene regulation in primate evolution (e.g., FOLH1 and PIEZO2 ). Finally, population genetics recapitulates macaque speciation and selective sweeps, highlighting potential genetic basis of reproduction and tail phenotype differences (e.g., STAB1 , SEMA3F , and HOXD13 ). In summary, the integrated analysis of genetic variation and population genetics in macaques greatly enhances our comprehension of lineage-specific phenotypes, adaptation, and primate evolution, thereby improving their biomedical applications in human diseases.

The complete and fully-phased diploid genome of a male Han Chinese.

Since the release of the complete human genome, the priority of human genomic study has now been shifting towards closing gaps in ethnic diversity. Here, we present a fully phased and well-annotated diploid human genome from a Han Chinese male individual (CN1), in which the assemblies of both haploids achieve the telomere-to-telomere (T2T) level. Comparison of this diploid genome with the CHM13 haploid T2T genome revealed significant variations in the centromere. Outside the centromere, we discovered 11,413 structural variations, including numerous novel ones. We also detected thousands of CN1 alleles that have accumulated high substitution rates and a few that have been under positive selection in the East Asian population. Further, we found that CN1 outperforms CHM13 as a reference genome in mapping and variant calling for the East Asian population owing to the distinct structural variants of the two references. Comparison of SNP calling for a large cohort of 8869 Chinese genomes using CN1 and CHM13 as reference respectively showed that the reference bias profoundly impacts rare SNP calling, with nearly 2 million rare SNPs miss-called with different reference genomes. Finally, applying the CN1 as a reference, we discovered 5.80 Mb and 4.21 Mb putative introgression sequences from Neanderthal and Denisovan, respectively, including many East Asian specific ones undetected using CHM13 as the reference. Our analyses reveal the advances of using CN1 as a reference for population genomic studies and paleo-genomic studies. This complete genome will serve as an alternative reference for future genomic studies on the East Asian population.

Convergent molecular evolution of thermogenesis and circadian rhythm in Arctic ruminants.

The muskox and reindeer are the only ruminants that have evolved to survive in harsh Arctic environments. However, the genetic basis of this Arctic adaptation remains largely unclear. Here, we compared a de novo assembled muskox genome with reindeer and other ruminant genomes to identify convergent amino acid substitutions, rapidly evolving genes and positively selected genes among the two Arctic ruminants. We found these candidate genes were mainly involved in brown adipose tissue (BAT) thermogenesis and circadian rhythm. Furthermore, by integrating transcriptomic data from goat adipose tissues (white and brown), we demonstrated that muskox and reindeer may have evolved modulating mitochondrion, lipid metabolism and angiogenesis pathways to enhance BAT thermogenesis. In addition, results from co-immunoprecipitation experiments prove that convergent amino acid substitution of the angiogenesis-related gene hypoxia-inducible factor 2alpha (HIF2A), resulting in weakening of its interaction with prolyl hydroxylase domain-containing protein 2 (PHD2), may increase angiogenesis of BAT. Altogether, our work provides new insights into the molecular mechanisms involved in Arctic adaptation.

Evolutionary analysis of a complete chicken genome.

Microchromosomes are prevalent in nonmammalian vertebrates [P. D. Waters et al., Proc. Natl. Acad. Sci. U.S.A. 118 (2021)], but a few of them are missing in bird genome assemblies. Here, we present a new chicken reference genome containing all autosomes, a Z and a W chromosome, with all gaps closed except for the W. We identified ten small microchromosomes (termed dot chromosomes) with distinct sequence and epigenetic features, among which six were newly assembled. Those dot chromosomes exhibit extremely high GC content and a high level of DNA methylation and are enriched for housekeeping genes. The pericentromeric heterochromatin of dot chromosomes is disproportionately large and continues to expand with the proliferation of satellite DNA and testis-expressed genes. Our analyses revealed that the 41-bp CNM repeat frequently forms higher-order repeats (HORs) at the centromeres of acrocentric chromosomes. The centromere core regions where the kinetochore attaches often encompass telomeric sequence (TTAGGG)n, and in a one of the dot chromosomes, the centromere core recruits an endogenous retrovirus (ERV). We further demonstrate that the W chromosome shares some common features with dot chromosomes, having large arrays of hypermethylated tandem repeats. Finally, using the complete chicken chromosome models, we reconstructed a fine picture of chordate karyotype evolution, revealing frequent chromosomal fusions before and after vertebrate whole-genome duplications. Our sequence and epigenetic characterization of chicken chromosomes shed insights into the understanding of vertebrate genome evolution and chromosome biology.

Evolution and diversification of Mountain voles (Rodentia: Cricetidae).

The systematics of the Cricetid genus Neodon have long been fraught with uncertainty due to sampling issues and a lack of comprehensive datasets. To gain better insights into the phylogeny and evolution of Neodon, we systematically sampled Neodon across the Hengduan and Himalayan Mountains, which cover most of its range in China. Analyses of skulls, teeth, and bacular structures revealed 15 distinct patterns corresponding to 15 species of Neodon. In addition to morphological analyses, we generated a high-quality reference genome for the mountain vole and generated whole-genome sequencing data for 47 samples. Phylogenomic analyses supported the recognition of six new species, revealing a long-term underestimation of Neodon diversity. We further identified positively selected genes potentially related to high-elevation adaptation. Together, our results illuminate how climate change caused the plateau to become the centre of Neodon origin and diversification and how mountain voles have adapted to the hypoxic high-altitude plateau environment.

A haplotype-resolved genome assembly of the Nile rat facilitates exploration of the genetic basis of diabetes.

BACKGROUND: The Nile rat (Avicanthis niloticus) is an important animal model because of its robust diurnal rhythm, a cone-rich retina, and a propensity to develop diet-induced diabetes without chemical or genetic modifications. A closer similarity to humans in these aspects, compared to the widely used Mus musculus and Rattus norvegicus models, holds the promise of better translation of research findings to the clinic. RESULTS: We report a 2.5 Gb, chromosome-level reference genome assembly with fully resolved parental haplotypes, generated with the Vertebrate Genomes Project (VGP). The assembly is highly contiguous, with contig N50 of 11.1 Mb, scaffold N50 of 83 Mb, and 95.2% of the sequence assigned to chromosomes. We used a novel workflow to identify 3613 segmental duplications and quantify duplicated genes. Comparative analyses revealed unique genomic features of the Nile rat, including some that affect genes associated with type 2 diabetes and metabolic dysfunctions. We discuss 14 genes that are heterozygous in the Nile rat or highly diverged from the house mouse. CONCLUSIONS: Our findings reflect the exceptional level of genomic resolution present in this assembly, which will greatly expand the potential of the Nile rat as a model organism.

Semi-automated assembly of high-quality diploid human reference genomes.

The current human reference genome, GRCh38, represents over 20 years of effort to generate a high-quality assembly, which has benefitted society1,2. However, it still has many gaps and errors, and does not represent a biological genome as it is a blend of multiple individuals3,4. Recently, a high-quality telomere-to-telomere reference, CHM13, was generated with the latest long-read technologies, but it was derived from a hydatidiform mole cell line with a nearly homozygous genome5. To address these limitations, the Human Pangenome Reference Consortium formed with the goal of creating high-quality, cost-effective, diploid genome assemblies for a pangenome reference that represents human genetic diversity6. Here, in our first scientific report, we determined which combination of current genome sequencing and assembly approaches yield the most complete and accurate diploid genome assembly with minimal manual curation. Approaches that used highly accurate long reads and parent-child data with graph-based haplotype phasing during assembly outperformed those that did not. Developing a combination of the top-performing methods, we generated our first high-quality diploid reference assembly, containing only approximately four gaps per chromosome on average, with most chromosomes within ±1% of the length of CHM13. Nearly 48% of protein-coding genes have non-synonymous amino acid changes between haplotypes, and centromeric regions showed the highest diversity. Our findings serve as a foundation for assembling near-complete diploid human genomes at scale for a pangenome reference to capture global genetic variation from single nucleotides to structural rearrangements.

Comparable Number of Genes Having Experienced Positive Selection among Great Ape Species.

Alleles that cause advantageous phenotypes with positive selection contribute to adaptive evolution. Investigations of positive selection in protein-coding genes rely on the accuracy of orthology, models, the quality of assemblies, and alignment. Here, based on the latest genome assemblies and gene annotations, we present a comparative analysis on positive selection in four great ape species and identify 211 high-confidence positively selected genes (PSGs). Even the differences in population size among these closely related great apes have resulted in differences in their ability to remove deleterious alleles and to adapt to changing environments, we found that they experienced comparable numbers of positive selection. We also uncovered that more than half of multigene families exhibited signals of positive selection, suggesting that imbalanced positive selection resulted in the functional divergence of duplicates. Moreover, at the expression level, although positive selection led to a more non-uniform pattern across tissues, the correlation between positive selection and expression patterns is diverse. Overall, this updated list of PSGs is of great significance for the further study of the phenotypic evolution in great apes.

Ancient and modern genomes unravel the evolutionary history of the rhinoceros family.

Only five species of the once-diverse Rhinocerotidae remain, making the reconstruction of their evolutionary history a challenge to biologists since Darwin. We sequenced genomes from five rhinoceros species (three extinct and two living), which we compared to existing data from the remaining three living species and a range of outgroups. We identify an early divergence between extant African and Eurasian lineages, resolving a key debate regarding the phylogeny of extant rhinoceroses. This early Miocene (∼16 million years ago [mya]) split post-dates the land bridge formation between the Afro-Arabian and Eurasian landmasses. Our analyses also show that while rhinoceros genomes in general exhibit low levels of genome-wide diversity, heterozygosity is lowest and inbreeding is highest in the modern species. These results suggest that while low genetic diversity is a long-term feature of the family, it has been particularly exacerbated recently, likely reflecting recent anthropogenic-driven population declines.

Evolutionary and biomedical insights from a marmoset diploid genome assembly.

The accurate and complete assembly of both haplotype sequences of a diploid organism is essential to understanding the role of variation in genome functions, phenotypes and diseases1. Here, using a trio-binning approach, we present a high-quality, diploid reference genome, with both haplotypes assembled independently at the chromosome level, for the common marmoset (Callithrix jacchus), an primate model system that is widely used in biomedical research2,3. The full spectrum of heterozygosity between the two haplotypes involves 1.36% of the genome-much higher than the 0.13% indicated by the standard estimation based on single-nucleotide heterozygosity alone. The de novo mutation rate is 0.43 × 10-8 per site per generation, and the paternal inherited genome acquired twice as many mutations as the maternal. Our diploid assembly enabled us to discover a recent expansion of the sex-differentiation region and unique evolutionary changes in the marmoset Y chromosome. In addition, we identified many genes with signatures of positive selection that might have contributed to the evolution of Callithrix biological features. Brain-related genes were highly conserved between marmosets and humans, although several genes experienced lineage-specific copy number variations or diversifying selection, with implications for the use of marmosets as a model system.

Value of peripheral blood neutrophil-to-lymphocyte ratio for clinical diagnosis and prognosis of elderly patients with chronic heart failure and atrial fibrillation.

AIM: We aimed to explore the value of peripheral blood neutrophil-to-lymphocyte ratio (NLR) for the clinical diagnosis and prognosis of elderly patients with chronic heart failure (CHF) and atrial fibrillation (AF). METHODS: A total of 248 eligible patients were followed up for five years, and divided into major adverse cardiovascular event (MACE) and non-MACE groups. The independent predictive factors for MACE were explored by multivariate logistic regression analysis. Based on quartile of NLR, they were divided into groups A to D. The duration of MACE was analysed using Kaplan-Meier survival curves. The diagnostic value of NLR for MACE was evaluated by receiver operating characteristic curves. RESULTS: Higher age, low-density lipoprotein cholesterol and NLR, lower left ventricular ejection fraction, diabetes and NYHA heart function class III and IV were independent predictive factors for MACE. The incidence of MACE rose with increasing NLR. Groups A to D had significantly different rates of acute myocardial infarction, severe arrhythmia and cardiac death (p < 0.05). The average duration of MACE in groups A to D were 49.31, 45.27, 43.63 and 40.34 months, respectively. CONCLUSIONS: The sensitivity and specificity of NLR for diagnosis of MACE were 72.39 and 86.18%, respectively. NLR was an independent predictive factor for MACE in these elderly patients with CHF and AF.

Efficient COI barcoding using high throughput single-end 400 bp sequencing.

BACKGROUND: Over the last decade, the rapid development of high-throughput sequencing platforms has accelerated species description and assisted morphological classification through DNA barcoding. However, the current high-throughput DNA barcoding methods cannot obtain full-length barcode sequences due to read length limitations (e.g. a maximum read length of 300 bp for the Illumina's MiSeq system), or are hindered by a relatively high cost or low sequencing output (e.g. a maximum number of eight million reads per cell for the PacBio's SEQUEL II system). RESULTS: Pooled cytochrome c oxidase subunit I (COI) barcodes from individual specimens were sequenced on the MGISEQ-2000 platform using the single-end 400 bp (SE400) module. We present a bioinformatic pipeline, HIFI-SE, that takes reads generated from the 5' and 3' ends of the COI barcode region and assembles them into full-length barcodes. HIFI-SE is written in Python and includes four function modules of filter, assign, assembly and taxonomy. We applied the HIFI-SE to a set of 845 samples (30 marine invertebrates, 815 insects) and delivered a total of 747 fully assembled COI barcodes as well as 70 Wolbachia and fungi symbionts. Compared to their corresponding Sanger sequences (72 sequences available), nearly all samples (71/72) were correctly and accurately assembled, including 46 samples that had a similarity score of 100% and 25 of ca. 99%. CONCLUSIONS: The HIFI-SE pipeline represents an efficient way to produce standard full-length barcodes, while the reasonable cost and high sensitivity of our method can contribute considerably more DNA barcodes under the same budget. Our method thereby advances DNA-based species identification from diverse ecosystems and increases the number of relevant applications.

Diverse coral reef invertebrates exhibit patterns of phylosymbiosis.

Microbiome assemblages of plants and animals often show a degree of correlation with host phylogeny; an eco-evolutionary pattern known as phylosymbiosis. Using 16S rRNA gene sequencing to profile the microbiome, paired with COI, 18S rRNA and ITS1 host phylogenies, phylosymbiosis was investigated in four groups of coral reef invertebrates (scleractinian corals, octocorals, sponges and ascidians). We tested three commonly used metrics to evaluate the extent of phylosymbiosis: (a) intraspecific versus interspecific microbiome variation, (b) topological comparisons between host phylogeny and hierarchical clustering (dendrogram) of host-associated microbial communities, and (c) correlation of host phylogenetic distance with microbial community dissimilarity. In all instances, intraspecific variation in microbiome composition was significantly lower than interspecific variation. Similarly, topological congruency between host phylogeny and the associated microbial dendrogram was more significant than would be expected by chance across all groups, except when using unweighted UniFrac distance (compared with weighted UniFrac and Bray-Curtis dissimilarity). Interestingly, all but the ascidians showed a significant positive correlation between host phylogenetic distance and associated microbial dissimilarity. Our findings provide new perspectives on the diverse nature of marine phylosymbioses and the complex roles of the microbiome in the evolution of marine invertebrates.

A draft genome assembly of spotted hyena, Crocuta crocuta.

The spotted hyena (Crocuta crocuta), one of the largest terrestrial predators native to sub-Saharan Africa, is well known for its matriarchal social system and large-sized social group in which larger females dominate smaller males. Spotted hyenas are highly adaptable predators as they both actively hunt prey and scavenge kills by other predators, and possess an enhanced hypercarnivorous dentition that allows them to crack open bones and thereby feed on nearly all parts of a carcass. Here, we present a high-quality genome assembly of C. crocuta that was generated using a hybrid assembly strategy with Illumina multi-size libraries. A genome of about 2.3 Gb was generated with a scaffold N50 length of 7.2 Mb. More than 35.28% genome region was identified as repetitive elements, and 22,747 protein-coding genes were identified in the genome, with 97.45% of these annotated by databases. This high-quality genome will provide an opportunity to gain insight into the evolution of social behavior and social cognition in mammals, as well as for population genetics and metagenomics studies.

Genomic Analyses Provide Insights Into the Evolutionary History and Genetic Diversity of Auricularia Species.

Species in the genus Auricularia play important roles for people's food and nutrition especially Auricularia cornea and A. heimuer. To understand their evolutionary history, genome structure, and population-level genetic variation, we performed a high-quality genome sequencing of Auricularia cornea and the corresponding comparative genomic analysis. The genome size of A. cornea was similar to Auricularia subglabra, but 1.5 times larger than that of A. heimuer. Several factors were responsible for genome size variation including gene numbers, repetitive elements, and gene lengths. Phylogenomic analysis revealed that the estimated divergence time between A. heimuer and other Auricularia is ∼79.1 million years ago (Mya), while the divergence between A. cornea and A. subglabra occurred in ∼54.8 Mya. Population genomic analysis also provided insight into the demographic history of A. cornea and A. heimuer, indicating that their populations fluctuated over time with global climate change during Marine Isotope Stage 5-2. Moreover, despite the highly similar external morphologies of A. cornea and A. heimuer, their genomic properties were remarkably different. The A. cornea genome only shared 14% homologous syntenic blocks with A. heimuer and possessed more genes encoding carbohydrate-active enzymes and secondary metabolite biosynthesis proteins. The cross-taxa transferability rates of simple sequence repeat (SSR) and insertion or deletion (InDel) markers within the genus Auricularia were also lower than that previously observed for species within the same genus. Taken together, these results indicate a high level of genetic differentiation between these two Auricularia species. Consequently, our study provides new insights into the genomic evolution and genetic differentiation of Auricularia species that will facilitate future genetic breeding.

MitoZ: a toolkit for animal mitochondrial genome assembly, annotation and visualization.

Mitochondrial genome (mitogenome) plays important roles in evolutionary and ecological studies. It becomes routine to utilize multiple genes on mitogenome or the entire mitogenomes to investigate phylogeny and biodiversity of focal groups with the onset of High Throughput Sequencing (HTS) technologies. We developed a mitogenome toolkit MitoZ, consisting of independent modules of de novo assembly, findMitoScaf (find Mitochondrial Scaffolds), annotation and visualization, that can generate mitogenome assembly together with annotation and visualization results from HTS raw reads. We evaluated its performance using a total of 50 samples of which mitogenomes are publicly available. The results showed that MitoZ can recover more full-length mitogenomes with higher accuracy compared to the other available mitogenome assemblers. Overall, MitoZ provides a one-click solution to construct the annotated mitogenome from HTS raw data and will facilitate large scale ecological and evolutionary studies. MitoZ is free open source software distributed under GPLv3 license and available at https://github.com/linzhi2013/MitoZ.

Genome Sequencing of Cladobotryum protrusum Provides Insights into the Evolution and Pathogenic Mechanisms of the Cobweb Disease Pathogen on Cultivated Mushroom.

Cladobotryum protrusum is one of the mycoparasites that cause cobweb disease on cultivated edible mushrooms. However, the molecular mechanisms of evolution and pathogenesis of C. protrusum on mushrooms are largely unknown. Here, we report a high-quality genome sequence of C. protrusum using the single-molecule, real-time sequencing platform of PacBio and perform a comparative analysis with closely related fungi in the family Hypocreaceae. The C. protrusum genome, the first complete genome to be sequenced in the genus Cladobotryum, is 39.09 Mb long, with an N50 of 4.97 Mb, encoding 11,003 proteins. The phylogenomic analysis confirmed its inclusion in Hypocreaceae, with its evolutionary divergence time estimated to be ~170.1 million years ago. The genome encodes a large and diverse set of genes involved in secreted peptidases, carbohydrate-active enzymes, cytochrome P450 enzymes, pathogen⁻host interactions, mycotoxins, and pigments. Moreover, C. protrusum harbors arrays of genes with the potential to produce bioactive secondary metabolites and stress response-related proteins that are significant for adaptation to hostile environments. Knowledge of the genome will foster a better understanding of the biology of C. protrusum and mycoparasitism in general, as well as help with the development of effective disease control strategies to minimize economic losses from cobweb disease in cultivated edible mushrooms.

The Changing Pattern of Population Structure of Staphylococcus aureus from Bacteremia in China from 2013 to 2016: ST239-030-MRSA Replaced by ST59-t437.

To investigate the epidemiology and genetic structure of Staphylococcus aureus bacteremia in China, a total of 416 isolates from 22 teaching hospitals in 12 cities from 2013 and 2016 were characterized by antibiogram analysis, multilocus sequence typing (MLST), spa typing and staphylococcal cassette chromosome mec (SCCmec) typing. The predominant meticillin-susceptible (MSSA) genotypes in 2013 were ST188 (19.1%), ST7 (8.7%), and ST398 (7.8%), respectively, and they continued to be the main genotypes in 2016. The prevalence of meticillin-resistant S. aureus (MRSA) were 36.5% (66/181) and 36.6% (86/235) in 2013 and 2016, respectively. Interestingly, the susceptibility rates of MRSA to rifampicin and fluoroquinolones increased significantly from 2013 to 2016 (P < 0.01), and this was associated with changes in genetic structure. ST239-t030-MRSA, the predominant genotype among all MRSAs in 2013 (34.8%), was replaced by ST59-t437-MRSA (15.1%) in 2016. Further analysis revealed that the ST239-t030-MRSA were more resistant to rifampicin, tetracycline and fluoroquinolones than ST59-t437-MRSA (P < 0.01). To further gain insight into the mechanisms underlying the changes of genetic structure, in vitro competition and fitness measurements were performed. Importantly, ST239-t030-MRSA displayed lower growth rate and lower competitive advantage compared to ST59-t437-MRSA. Together, our findings reveal that fitness advantage of ST59-t437-MRSA over ST239-t030-MRSA may lead to changes in genetic structure and increased susceptibility of MRSA to rifampicin and fluoroquinolones in Chinese patients with S. aureus bacteremia. Our study supports temporal dynamics in MRSA clone diversities, further providing critical insights into the importance of continued monitoring of MRSA.

Characterization of viral RNA splicing using whole-transcriptome datasets from host species.

RNA alternative splicing (AS) is an important post-transcriptional mechanism enabling single genes to produce multiple proteins. It has been well demonstrated that viruses deploy host AS machinery for viral protein productions. However, knowledge on viral AS is limited to a few disease-causing viruses in model species. Here we report a novel approach to characterizing viral AS using whole transcriptome dataset from host species. Two insect transcriptomes (Acheta domesticus and Planococcus citri) generated in the 1,000 Insect Transcriptome Evolution (1KITE) project were used as a proof of concept using the new pipeline. Two closely related densoviruses (Acheta domesticus densovirus, AdDNV, and Planococcus citri densovirus, PcDNV, Ambidensovirus, Densovirinae, Parvoviridae) were detected and analyzed for AS patterns. The results suggested that although the two viruses shared major AS features, dramatic AS divergences were observed. Detailed analysis of the splicing junctions showed clusters of AS events occurred in two regions of the virus genome, demonstrating that transcriptome analysis could gain valuable insights into viral splicing. When applied to large-scale transcriptomics projects with diverse taxonomic sampling, our new method is expected to rapidly expand our knowledge on RNA splicing mechanisms for a wide range of viruses.