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.

Platypus and echidna genomes reveal mammalian biology and evolution.

Egg-laying mammals (monotremes) are the only extant mammalian outgroup to therians (marsupial and eutherian animals) and provide key insights into mammalian evolution1,2. Here we generate and analyse reference genomes of the platypus (Ornithorhynchus anatinus) and echidna (Tachyglossus aculeatus), which represent the only two extant monotreme lineages. The nearly complete platypus genome assembly has anchored almost the entire genome onto chromosomes, markedly improving the genome continuity and gene annotation. Together with our echidna sequence, the genomes of the two species allow us to detect the ancestral and lineage-specific genomic changes that shape both monotreme and mammalian evolution. We provide evidence that the monotreme sex chromosome complex originated from an ancestral chromosome ring configuration. The formation of such a unique chromosome complex may have been facilitated by the unusually extensive interactions between the multi-X and multi-Y chromosomes that are shared by the autosomal homologues in humans. Further comparative genomic analyses unravel marked differences between monotremes and therians in haptoglobin genes, lactation genes and chemosensory receptor genes for smell and taste that underlie the ecological adaptation of monotremes.

Towards complete and error-free genome assemblies of all vertebrate species.

High-quality and complete reference genome assemblies are fundamental for the application of genomics to biology, disease, and biodiversity conservation. However, such assemblies are available for only a few non-microbial species1-4. To address this issue, the international Genome 10K (G10K) consortium5,6 has worked over a five-year period to evaluate and develop cost-effective methods for assembling highly accurate and nearly complete reference genomes. Here we present lessons learned from generating assemblies for 16 species that represent six major vertebrate lineages. We confirm that long-read sequencing technologies are essential for maximizing genome quality, and that unresolved complex repeats and haplotype heterozygosity are major sources of assembly error when not handled correctly. Our assemblies correct substantial errors, add missing sequence in some of the best historical reference genomes, and reveal biological discoveries. These include the identification of many false gene duplications, increases in gene sizes, chromosome rearrangements that are specific to lineages, a repeated independent chromosome breakpoint in bat genomes, and a canonical GC-rich pattern in protein-coding genes and their regulatory regions. Adopting these lessons, we have embarked on the Vertebrate Genomes Project (VGP), an international effort to generate high-quality, complete reference genomes for all of the roughly 70,000 extant vertebrate species and to help to enable a new era of discovery across the life sciences.

Divergent sensory and immune gene evolution in sea turtles with contrasting demographic and life histories.

Sea turtles represent an ancient lineage of marine vertebrates that evolved from terrestrial ancestors over 100 Mya. The genomic basis of the unique physiological and ecological traits enabling these species to thrive in diverse marine habitats remains largely unknown. Additionally, many populations have drastically declined due to anthropogenic activities over the past two centuries, and their recovery is a high global conservation priority. We generated and analyzed high-quality reference genomes for the leatherback (Dermochelys coriacea) and green (Chelonia mydas) turtles, representing the two extant sea turtle families. These genomes are highly syntenic and homologous, but localized regions of noncollinearity were associated with higher copy numbers of immune, zinc-finger, and olfactory receptor (OR) genes in green turtles, with ORs related to waterborne odorants greatly expanded in green turtles. Our findings suggest that divergent evolution of these key gene families may underlie immunological and sensory adaptations assisting navigation, occupancy of neritic versus pelagic environments, and diet specialization. Reduced collinearity was especially prevalent in microchromosomes, with greater gene content, heterozygosity, and genetic distances between species, supporting their critical role in vertebrate evolutionary adaptation. Finally, diversity and demographic histories starkly contrasted between species, indicating that leatherback turtles have had a low yet stable effective population size, exhibit extremely low diversity compared with other reptiles, and harbor a higher genetic load compared with green turtles, reinforcing concern over their persistence under future climate scenarios. These genomes provide invaluable resources for advancing our understanding of evolution and conservation best practices in an imperiled vertebrate lineage.

Pronounced early differentiation underlies zebra finch gonadal germ cell development.

The diversity of germ cell developmental strategies has been well documented across many vertebrate clades. However, much of our understanding of avian primordial germ cell (PGC) specification and differentiation has derived from only one species, the chicken (Gallus gallus). Of the three major classes of birds, chickens belong to Galloanserae, representing less than 4% of species, while nearly 95% of extant bird species belong to Neoaves. This represents a significant gap in our knowledge of germ cell development across avian species, hampering efforts to adapt genome editing and reproductive technologies developed in chicken to other birds. We therefore applied single-cell RNA sequencing to investigate inter-species differences in germ cell development between chicken and zebra finch (Taeniopygia castanotis), a Neoaves songbird species and a common model of vocal learning. Analysis of early embryonic male and female gonads revealed the presence of two distinct early germ cell types in zebra finch and only one in chicken. Both germ cell types expressed zebra finch Germline Restricted Chromosome (GRC) genes, present only in songbirds among birds. One of the zebra finch germ cell types expressed the canonical PGC markers, as did chicken, but with expression differences in several signaling pathways and biological processes. The second zebra finch germ cell cluster was marked by proliferation and fate determination markers, indicating beginning of differentiation. Notably, these two zebra finch germ cell populations were present in both male and female zebra finch gonads as early as HH25. Using additional chicken developmental stages, similar germ cell heterogeneity was identified in the more developed gonads of females, but not males. Overall, our study demonstrates a substantial heterochrony in zebra finch germ cell development compared to chicken, indicating a richer diversity of avian germ cell developmental strategies than previously known.

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.

Identification and characterization of primordial germ cells in a vocal learning Neoaves species, the zebra finch.

The zebra finch has been used as a valuable vocal learning animal model for human spoken language. It is representative of vocal learning songbirds specifically, which comprise half of all bird species, and of Neoaves broadly, which comprise 95% of all bird species. Although transgenesis in the zebra finch has been accomplished, it is with a very low efficiency of germ-line transmission and far from the efficiency with a more genetically tractable but vocal nonlearning species, the chicken (a Galloanseriformes). To improve germ-line transmission in the zebra finch, we identified and characterized its primordial germ cells (PGCs) and compared them with chicken. We found striking differences between the 2 species, including that zebra finch PGCs were more numerous, more widely distributed in early embryos before colonization into the gonads, had slower timing of colonization, and had a different developmental gene-expression program. We improved conditions for isolating and culturing zebra finch PGCs in vitro and were able to transfect them with gene-expression vectors and incorporate them into the gonads of host embryos. Our findings demonstrate important differences in the PGCs of the zebra finch and advance the first stage of creating PGC-mediated germ-line transgenics of a vocal learning species.-Jung, K. M., Kim, Y. M., Keyte, A. L., Biegler, M. T., Rengaraj, D., Lee, H. J., Mello, C. V., Velho, T. A. F., Fedrigo, O., Haase, B., Jarvis, E. D., Han, J. Y. Identification and characterization of primordial germ cells in a vocal learning Neoaves species, the zebra finch.

The Slicer Activity of ARGONAUTE1 Is Required Specifically for the Phasing, Not Production, of Trans-Acting Short Interfering RNAs in Arabidopsis.

ARGONAUTE1 (AGO1) mediates posttranscriptional silencing by microRNAs (miRNAs) and short interfering RNAS (siRNAs). AGO1-catalyzed RNA cleavage (slicing) represses miRNA targets, but current models also highlight the roles of slicing in formation of siRNAs and siRNA-AGO1 complexes. miRNA-guided slicing is required for biogenesis of phased, trans-acting siRNAs (tasiRNAs), whose cleaved precursor fragments are converted to double-stranded RNA by RNA-dependent RNA polymerase 6 (RDR6). In addition, unwinding of duplex siRNA bound to AGO1 requires passenger strand cleavage in vitro. In this study, we analyze how mutation of four metal ion-coordinating residues of Arabidopsis thaliana AGO1 affects slicer activity in vitro and siRNA function in vivo. We show that while all four residues are required for slicer activity, they do not contribute equally to catalysis. Moreover, passenger strand cleavage is required for assembly of active AGO1-siRNA complexes in vivo, and many AGO1-bound siRNAs are trimmed in the absence of slicer activity. Remarkably, seedlings defective in AGO1 slicer activity produce abundant siRNAs from tasiRNA loci in vivo. These siRNAs depend on RDR6 and SUPPRESSOR OF GENE SILENCING3, but unlike wild-type tasiRNAs, they are unphased. These results demonstrate that slicing is solely required for phase definition of tasiRNAs, and they strongly support recruitment of RDR6 by AGO1 rather than by cleavage fragments.

miRNA92a targets KLF2 and the phosphatase PTEN signaling to promote human T follicular helper precursors in T1D islet autoimmunity.

Aberrant immune activation mediated by T effector cell populations is pivotal in the onset of autoimmunity in type 1 diabetes (T1D). T follicular helper (TFH) cells are essential in the induction of high-affinity antibodies, and their precursor memory compartment circulates in the blood. The role of TFH precursors in the onset of islet autoimmunity and signaling pathways regulating their differentiation is incompletely understood. Here, we provide direct evidence that during onset of islet autoimmunity, the insulin-specific target T-cell population is enriched with a C-X-C chemokine receptor type 5 (CXCR5)+CD4+ TFH precursor phenotype. During onset of islet autoimmunity, the frequency of TFH precursors was controlled by high expression of microRNA92a (miRNA92a). miRNA92a-mediated TFH precursor induction was regulated by phosphatase and tension homolog (PTEN) - phosphoinositol-3-kinase (PI3K) signaling involving PTEN and forkhead box protein O1 (Foxo1), supporting autoantibody generation and triggering the onset of islet autoimmunity. Moreover, we identify Krueppel-like factor 2 (KLF2) as a target of miRNA92a in regulating human TFH precursor induction. Importantly, a miRNA92a antagomir completely blocked induction of human TFH precursors in vitro. More importantly, in vivo application of a miRNA92a antagomir to nonobese diabetic (NOD) mice with ongoing islet autoimmunity resulted in a significant reduction of TFH precursors in peripheral blood and pancreatic lymph nodes. Moreover, miRNA92a antagomir application reduced immune infiltration and activation in pancreata of NOD mice as well as humanized NOD Scid IL2 receptor gamma chain knockout (NSG) human leucocyte antigen (HLA)-DQ8 transgenic animals. We therefore propose that miRNA92a and the PTEN-PI3K-KLF2 signaling network could function as targets for innovative precision medicines to reduce T1D islet autoimmunity.

The differential expression of alternatively polyadenylated transcripts is a common stress-induced response mechanism that modulates mammalian mRNA expression in a quantitative and qualitative fashion.

Stress adaptation plays a pivotal role in biological processes and requires tight regulation of gene expression. In this study, we explored the effect of cellular stress on mRNA polyadenylation and investigated the implications of regulated polyadenylation site usage on mammalian gene expression. High-confidence polyadenylation site mapping combined with global pre-mRNA and mRNA expression profiling revealed that stress induces an accumulation of genes with differentially expressed polyadenylated mRNA isoforms in human cells. Specifically, stress provokes a global trend in polyadenylation site usage toward decreased utilization of promoter-proximal poly(A) sites in introns or ORFs and increased utilization of promoter-distal polyadenylation sites in intergenic regions. This extensively affects gene expression beyond regulating mRNA abundance by changing mRNA length and by altering the configuration of open reading frames. Our study highlights the impact of post-transcriptional mechanisms on stress-dependent gene regulation and reveals the differential expression of alternatively polyadenylated transcripts as a common stress-induced mechanism in mammalian cells.

Distinct patterns of genetic variation at low-recombining genomic regions represent haplotype structure.

Genomic regions sometimes show patterns of genetic variation distinct from the genome-wide population structure. Such deviations have often been interpreted to represent effects of selection. However, systematic investigation of whether and how non-selective factors, such as recombination rates, can affect distinct patterns has been limited. Here, we associate distinct patterns of genetic variation with reduced recombination rates in a songbird, the Eurasian blackcap (Sylvia atricapilla), using a new reference genome assembly, whole-genome resequenc- ing data and recombination maps. We find that distinct patterns of genetic variation reflect haplotype structure at genomic regions with different prevalence of reduced recombination rate across populations. At low-recombining regions shared in most populations, distinct patterns reflect conspicuous haplotypes segregating in multiple populations. At low-recombining regions found only in a few populations, distinct patterns represent variance among cryptic haplotypes within the low-recombining populations. With simulations, we confirm that these distinct patterns evolve neutrally by reduced recombination rate, on which the effects of selection can be overlaid. Our results highlight that distinct patterns of genetic variation can emerge through evolutionary reduction of local recombination rate. The recombination landscape as an evolvable trait therefore plays an important role determining the heterogeneous distribution of genetic variation along the genome.

A genomic basis of vocal rhythm in birds.

Vocal rhythm plays a fundamental role in sexual selection and species recognition in birds, but little is known of its genetic basis due to the confounding effect of vocal learning in model systems. Uncovering its genetic basis could facilitate identifying genes potentially important in speciation. Here we investigate the genomic underpinnings of rhythm in vocal non-learning Pogoniulus tinkerbirds using 135 individual whole genomes distributed across a southern African hybrid zone. We find rhythm speed is associated with two genes that are also known to affect human speech, Neurexin-1 and Coenzyme Q8A. Models leveraging ancestry reveal these candidate loci also impact rhythmic stability, a trait linked with motor performance which is an indicator of quality. Character displacement in rhythmic stability suggests possible reinforcement against hybridization, supported by evidence of asymmetric assortative mating in the species producing faster, more stable rhythms. Because rhythm is omnipresent in animal communication, candidate genes identified here may shape vocal rhythm across birds and other vertebrates.

Low mutation rate in epaulette sharks is consistent with a slow rate of evolution in sharks.

Sharks occupy diverse ecological niches and play critical roles in marine ecosystems, often acting as apex predators. They are considered a slow-evolving lineage and have been suggested to exhibit exceptionally low cancer rates. These two features could be explained by a low nuclear mutation rate. Here, we provide a direct estimate of the nuclear mutation rate in the epaulette shark (Hemiscyllium ocellatum). We generate a high-quality reference genome, and resequence the whole genomes of parents and nine offspring to detect de novo mutations. Using stringent criteria, we estimate a mutation rate of 7×10-10 per base pair, per generation. This represents one of the lowest directly estimated mutation rates for any vertebrate clade, indicating that this basal vertebrate group is indeed a slowly evolving lineage whose ability to restore genetic diversity following a sustained population bottleneck may be hampered by a low mutation rate.

An improved germline genome assembly for the sea lamprey Petromyzon marinus illuminates the evolution of germline-specific chromosomes.

Programmed DNA loss is a gene silencing mechanism that is employed by several vertebrate and nonvertebrate lineages, including all living jawless vertebrates and songbirds. Reconstructing the evolution of somatically eliminated (germline-specific) sequences in these species has proven challenging due to a high content of repeats and gene duplications in eliminated sequences and a corresponding lack of highly accurate and contiguous assemblies for these regions. Here, we present an improved assembly of the sea lamprey (Petromyzon marinus) genome that was generated using recently standardized methods that increase the contiguity and accuracy of vertebrate genome assemblies. This assembly resolves highly contiguous, somatically retained chromosomes and at least one germline-specific chromosome, permitting new analyses that reconstruct the timing, mode, and repercussions of recruitment of genes to the germline-specific fraction. These analyses reveal major roles of interchromosomal segmental duplication, intrachromosomal duplication, and positive selection for germline functions in the long-term evolution of germline-specific chromosomes.

The admixed brushtail possum genome reveals invasion history in New Zealand and novel imprinted genes.

Combining genome assembly with population and functional genomics can provide valuable insights to development and evolution, as well as tools for species management. Here, we present a chromosome-level genome assembly of the common brushtail possum (Trichosurus vulpecula), a model marsupial threatened in parts of their native range in Australia, but also a major introduced pest in New Zealand. Functional genomics reveals post-natal activation of chemosensory and metabolic genes, reflecting unique adaptations to altricial birth and delayed weaning, a hallmark of marsupial development. Nuclear and mitochondrial analyses trace New Zealand possums to distinct Australian subspecies, which have subsequently hybridised. This admixture allowed phasing of parental alleles genome-wide, ultimately revealing at least four genes with imprinted, parent-specific expression not yet detected in other species (MLH1, EPM2AIP1, UBP1 and GPX7). We find that reprogramming of possum germline imprints, and the wider epigenome, is similar to eutherian mammals except onset occurs after birth. Together, this work is useful for genetic-based control and conservation of possums, and contributes to understanding of the evolution of novel mammalian epigenetic traits.

A chromosome-level reference genome and pangenome for barn swallow population genomics.

Insights into the evolution of non-model organisms are limited by the lack of reference genomes of high accuracy, completeness, and contiguity. Here, we present a chromosome-level, karyotype-validated reference genome and pangenome for the barn swallow (Hirundo rustica). We complement these resources with a reference-free multialignment of the reference genome with other bird genomes and with the most comprehensive catalog of genetic markers for the barn swallow. We identify potentially conserved and accelerated genes using the multialignment and estimate genome-wide linkage disequilibrium using the catalog. We use the pangenome to infer core and accessory genes and to detect variants using it as a reference. Overall, these resources will foster population genomics studies in the barn swallow, enable detection of candidate genes in comparative genomics studies, and help reduce bias toward a single reference genome.

Microarray-based multicycle-enrichment of genomic subsets for targeted next-generation sequencing.

The lack of efficient high-throughput methods for enrichment of specific sequences from genomic DNA represents a key bottleneck in exploiting the enormous potential of next-generation sequencers. Such methods would allow for a systematic and targeted analysis of relevant genomic regions. Recent studies reported sequence enrichment using a hybridization step to specific DNA capture probes as a possible solution to the problem. However, so far no method has provided sufficient depths of coverage for reliable base calling over the entire target regions. We report a strategy to multiply the enrichment performance and consequently improve depth and breadth of coverage for desired target sequences by applying two iterative cycles of hybridization with microfluidic Geniom biochips. Using this strategy, we enriched and then sequenced the cancer-related genes BRCA1 and TP53 and a set of 1000 individual dbSNP regions of 500 bp using Illumina technology. We achieved overall enrichment factors of up to 1062-fold and average coverage depths of 470-fold. Combined with high coverage uniformity, this resulted in nearly complete consensus coverages with >86% of target region covered at 20-fold or higher. Analysis of SNP calling accuracies after enrichment revealed excellent concordance, with the reference sequence closely mirroring the previously reported performance of Illumina sequencing conducted without sequence enrichment.

Induction of an immortalized songbird cell line allows for gene characterization and knockout by CRISPR-Cas9.

The zebra finch is one of the most commonly studied songbirds in biology, particularly in genomics, neuroscience and vocal communication. However, this species lacks a robust cell line for molecular biology research and reagent optimization. We generated a cell line, designated CFS414, from zebra finch embryonic fibroblasts using the SV40 large and small T antigens. This cell line demonstrates an improvement over previous songbird cell lines through continuous and density-independent growth, allowing for indefinite culture and monoclonal line derivation. Cytogenetic, genomic, and transcriptomic profiling established the provenance of this cell line and identified the expression of genes relevant to ongoing songbird research. Using this cell line, we disrupted endogenous gene sequences using S.aureus Cas9 and confirmed a stress-dependent localization response of a song system specialized gene, SAP30L. The utility of CFS414 cells enhances the comprehensive molecular potential of the zebra finch and validates cell immortalization strategies in a songbird species.

A high-quality, long-read genome assembly of the endangered ring-tailed lemur (Lemur catta).

BACKGROUND: The ring-tailed lemur (Lemur catta) is a charismatic strepsirrhine primate endemic to Madagascar. These lemurs are of particular interest, given their status as a flagship species and widespread publicity in the popular media. Unfortunately, a recent population decline has resulted in the census population decreasing to <2,500 individuals in the wild, and the species's classification as an endangered species by the IUCN. As is the case for most strepsirrhine primates, only a limited amount of genomic research has been conducted on L. catta, in part owing to the lack of genomic resources. RESULTS: We generated a new high-quality reference genome assembly for L. catta (mLemCat1) that conforms to the standards of the Vertebrate Genomes Project. This new long-read assembly is composed of Pacific Biosciences continuous long reads (CLR data), Optical Mapping Bionano reads, Arima HiC data, and 10X linked reads. The contiguity and completeness of the assembly are extremely high, with scaffold and contig N50 values of 90.982 and 10.570 Mb, respectively. Additionally, when compared to other high-quality primate assemblies, L. catta has the lowest reported number of Alu elements, which results predominantly from a lack of AluS and AluY elements. CONCLUSIONS: mLemCat1 is an excellent genomic resource not only for the ring-tailed lemur community, but also for other members of the Lemuridae family, and is the first very long read assembly for a strepsirrhine.