Comparative epigenomic annotation of regulatory DNA.

Despite the explosive growth of genomic data, functional annotation of regulatory sequences remains difficult. Here, we introduce "comparative epigenomics"-interspecies comparison of DNA and histone modifications-as an approach for annotation of the regulatory genome. We measured in human, mouse, and pig pluripotent stem cells the genomic distributions of cytosine methylation, H2A.Z, H3K4me1/2/3, H3K9me3, H3K27me3, H3K27ac, H3K36me3, transcribed RNAs, and P300, TAF1, OCT4, and NANOG binding. We observed that epigenomic conservation was strong in both rapidly evolving and slowly evolving DNA sequences, but not in neutrally evolving sequences. In contrast, evolutionary changes of the epigenome and the transcriptome exhibited a linear correlation. We suggest that the conserved colocalization of different epigenomic marks can be used to discover regulatory sequences. Indeed, seven pairs of epigenomic marks identified exhibited regulatory functions during differentiation of embryonic stem cells into mesendoderm cells. Thus, comparative epigenomics reveals regulatory features of the genome that cannot be discerned from sequence comparisons alone.

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.

Conservation of chromatin conformation in carnivores.

High throughput chromatin conformation capture (Hi-C) of leukocyte DNA was used to investigate the evolutionary stability of chromatin conformation at the chromosomal level in 11 species from three carnivore families: Felidae, Canidae, and Ursidae. Chromosome-scale scaffolds (C-scaffolds) of each species were initially used for whole-genome alignment to a reference genome within each family. This approach established putative orthologous relationships between C-scaffolds among the different species. Hi-C contact maps for all C-scaffolds were then visually compared and found to be distinct for a given reference chromosome or C-scaffold within a species and indistinguishable for orthologous C-scaffolds having a 1:1 relationship within a family. The visual patterns within families were strongly supported by eigenvectors from the Hi-C contact maps. Analysis of Hi-C contact maps and eigenvectors across the three carnivore families revealed that most cross-family orthologous subchromosomal fragments have a conserved three-dimensional (3D) chromatin structure and thus have been under strong evolutionary constraint for ∼54 My of carnivore evolution. The most pronounced differences in chromatin conformation were observed for the X chromosome and the red fox genome, whose chromosomes have undergone extensive rearrangements relative to other canids. We also demonstrate that Hi-C contact map pattern analysis can be used to accurately identify orthologous relationships between C-scaffolds and chromosomes, a method we termed "3D comparative scaffotyping." This method provides a powerful means for estimating karyotypes in de novo sequenced species that have unknown karyotype and no physical mapping information.

Ethical, legal, and social issues in the Earth BioGenome Project.

The Earth BioGenome Project (EBP) is an audacious endeavor to obtain whole-genome sequences of representatives from all eukaryotic species on Earth. In addition to the project's technical and organizational challenges, it also faces complicated ethical, legal, and social issues. This paper, from members of the EBP's Ethical, Legal, and Social Issues (ELSI) Committee, catalogs these ELSI concerns arising from EBP. These include legal issues, such as sample collection and permitting; the applicability of international treaties, such as the Convention on Biological Diversity and the Nagoya Protocol; intellectual property; sample accessioning; and biosecurity and ethical issues, such as sampling from the territories of Indigenous peoples and local communities, the protection of endangered species, and cross-border collections, among several others. We also comment on the intersection of digital sequence information and data rights. More broadly, this list of ethical, legal, and social issues for large-scale genomic sequencing projects may be useful in the consideration of ethical frameworks for future projects. While we do not-and cannot-provide simple, overarching solutions for all the issues raised here, we conclude our perspective by beginning to chart a path forward for EBP's work.

Why sequence all eukaryotes?

Life on Earth has evolved from initial simplicity to the astounding complexity we experience today. Bacteria and archaea have largely excelled in metabolic diversification, but eukaryotes additionally display abundant morphological innovation. How have these innovations come about and what constraints are there on the origins of novelty and the continuing maintenance of biodiversity on Earth? The history of life and the code for the working parts of cells and systems are written in the genome. The Earth BioGenome Project has proposed that the genomes of all extant, named eukaryotes-about 2 million species-should be sequenced to high quality to produce a digital library of life on Earth, beginning with strategic phylogenetic, ecological, and high-impact priorities. Here we discuss why we should sequence all eukaryotic species, not just a representative few scattered across the many branches of the tree of life. We suggest that many questions of evolutionary and ecological significance will only be addressable when whole-genome data representing divergences at all of the branchings in the tree of life or all species in natural ecosystems are available. We envisage that a genomic tree of life will foster understanding of the ongoing processes of speciation, adaptation, and organismal dependencies within entire ecosystems. These explorations will resolve long-standing problems in phylogenetics, evolution, ecology, conservation, agriculture, bioindustry, and medicine.

Evolution of the ancestral mammalian karyotype and syntenic regions.

Decrypting the rearrangements that drive mammalian chromosome evolution is critical to understanding the molecular bases of speciation, adaptation, and disease susceptibility. Using 8 scaffolded and 26 chromosome-scale genome assemblies representing 23/26 mammal orders, we computationally reconstructed ancestral karyotypes and syntenic relationships at 16 nodes along the mammalian phylogeny. Three different reference genomes (human, sloth, and cattle) representing phylogenetically distinct mammalian superorders were used to assess reference bias in the reconstructed ancestral karyotypes and to expand the number of clades with reconstructed genomes. The mammalian ancestor likely had 19 pairs of autosomes, with nine of the smallest chromosomes shared with the common ancestor of all amniotes (three still conserved in extant mammals), demonstrating a striking conservation of synteny for ∼320 My of vertebrate evolution. The numbers and types of chromosome rearrangements were classified for transitions between the ancestral mammalian karyotype, descendent ancestors, and extant species. For example, 94 inversions, 16 fissions, and 14 fusions that occurred over 53 My differentiated the therian from the descendent eutherian ancestor. The highest breakpoint rate was observed between the mammalian and therian ancestors (3.9 breakpoints/My). Reconstructed mammalian ancestor chromosomes were found to have distinct evolutionary histories reflected in their rates and types of rearrangements. The distributions of genes, repetitive elements, topologically associating domains, and actively transcribed regions in multispecies homologous synteny blocks and evolutionary breakpoint regions indicate that purifying selection acted over millions of years of vertebrate evolution to maintain syntenic relationships of developmentally important genes and regulatory landscapes of gene-dense chromosomes.

Standards recommendations for the Earth BioGenome Project.

A global international initiative, such as the Earth BioGenome Project (EBP), requires both agreement and coordination on standards to ensure that the collective effort generates rapid progress toward its goals. To this end, the EBP initiated five technical standards committees comprising volunteer members from the global genomics scientific community: Sample Collection and Processing, Sequencing and Assembly, Annotation, Analysis, and IT and Informatics. The current versions of the resulting standards documents are available on the EBP website, with the recognition that opportunities, technologies, and challenges may improve or change in the future, requiring flexibility for the EBP to meet its goals. Here, we describe some highlights from the proposed standards, and areas where additional challenges will need to be met.

Broad host range of SARS-CoV-2 predicted by comparative and structural analysis of ACE2 in vertebrates.

The novel coronavirus severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the cause of COVID-19. The main receptor of SARS-CoV-2, angiotensin I converting enzyme 2 (ACE2), is now undergoing extensive scrutiny to understand the routes of transmission and sensitivity in different species. Here, we utilized a unique dataset of ACE2 sequences from 410 vertebrate species, including 252 mammals, to study the conservation of ACE2 and its potential to be used as a receptor by SARS-CoV-2. We designed a five-category binding score based on the conservation properties of 25 amino acids important for the binding between ACE2 and the SARS-CoV-2 spike protein. Only mammals fell into the medium to very high categories and only catarrhine primates into the very high category, suggesting that they are at high risk for SARS-CoV-2 infection. We employed a protein structural analysis to qualitatively assess whether amino acid changes at variable residues would be likely to disrupt ACE2/SARS-CoV-2 spike protein binding and found the number of predicted unfavorable changes significantly correlated with the binding score. Extending this analysis to human population data, we found only rare (frequency <0.001) variants in 10/25 binding sites. In addition, we found significant signals of selection and accelerated evolution in the ACE2 coding sequence across all mammals, and specific to the bat lineage. Our results, if confirmed by additional experimental data, may lead to the identification of intermediate host species for SARS-CoV-2, guide the selection of animal models of COVID-19, and assist the conservation of animals both in native habitats and in human care.

Evolution of gene regulation in ruminants differs between evolutionary breakpoint regions and homologous synteny blocks.

The role of chromosome rearrangements in driving evolution has been a long-standing question of evolutionary biology. Here we focused on ruminants as a model to assess how rearrangements may have contributed to the evolution of gene regulation. Using reconstructed ancestral karyotypes of Cetartiodactyls, Ruminants, Pecorans, and Bovids, we traced patterns of gross chromosome changes. We found that the lineage leading to the ruminant ancestor after the split from other cetartiodactyls was characterized by mostly intrachromosomal changes, whereas the lineage leading to the pecoran ancestor (including all livestock ruminants) included multiple interchromosomal changes. We observed that the liver cell putative enhancers in the ruminant evolutionary breakpoint regions are highly enriched for DNA sequences under selective constraint acting on lineage-specific transposable elements (TEs) and a set of 25 specific transcription factor (TF) binding motifs associated with recently active TEs. Coupled with gene expression data, we found that genes near ruminant breakpoint regions exhibit more divergent expression profiles among species, particularly in cattle, which is consistent with the phylogenetic origin of these breakpoint regions. This divergence was significantly greater in genes with enhancers that contain at least one of the 25 specific TF binding motifs and located near bovidae-to-cattle lineage breakpoint regions. Taken together, by combining ancestral karyotype reconstructions with analysis of cis regulatory element and gene expression evolution, our work demonstrated that lineage-specific regulatory elements colocalized with gross chromosome rearrangements may have provided valuable functional modifications that helped to shape ruminant evolution.

Fine-tuned adaptation of embryo-endometrium pairs at implantation revealed by transcriptome analyses in Bos taurus.

Interactions between embryo and endometrium at implantation are critical for the progression of pregnancy. These reciprocal actions involve exchange of paracrine signals that govern implantation and placentation. However, it remains unknown how these interactions between the conceptus and the endometrium are coordinated at the level of an individual pregnancy. Under the hypothesis that gene expression in endometrium is dependent on gene expression of extraembryonic tissues and genes expressed in extraembryonic tissues are dependent of genes expressed in the endometrium, we performed an integrative analysis of transcriptome profiles of paired extraembryonic tissue and endometria obtained from cattle (Bos taurus) pregnancies initiated by artificial insemination. We quantified strong dependence (|r| > 0.95, empirical false discovery rate [eFDR] < 0.01) in transcript abundance of genes expressed in the extraembryonic tissues and genes expressed in the endometrium. The profiles of connectivity revealed distinct coexpression patterns of extraembryonic tissues with caruncular and intercaruncular areas of the endometrium. Notably, a subset of highly coexpressed genes between extraembryonic tissue (n = 229) and caruncular areas of the endometrium (n = 218, r > 0.9999, eFDR < 0.001) revealed a blueprint of gene expression specific to each pregnancy. Gene ontology analyses of genes coexpressed between extraembryonic tissue and endometrium revealed significantly enriched modules with critical contribution for implantation and placentation, including "in utero embryonic development," "placenta development," and "regulation of transcription." Coexpressing modules were remarkably specific to caruncular or intercaruncular areas of the endometrium. The quantitative association between genes expressed in extraembryonic tissue and endometrium emphasize a coordinated communication between these two entities in mammals. We provide evidence that implantation in mammalian pregnancy relies on the ability of the extraembryonic tissue and the endometrium to develop a fine-tuned adaptive response characteristic of each pregnancy.

Earth BioGenome Project: Sequencing life for the future of life.

Increasing our understanding of Earth's biodiversity and responsibly stewarding its resources are among the most crucial scientific and social challenges of the new millennium. These challenges require fundamental new knowledge of the organization, evolution, functions, and interactions among millions of the planet's organisms. Herein, we present a perspective on the Earth BioGenome Project (EBP), a moonshot for biology that aims to sequence, catalog, and characterize the genomes of all of Earth's eukaryotic biodiversity over a period of 10 years. The outcomes of the EBP will inform a broad range of major issues facing humanity, such as the impact of climate change on biodiversity, the conservation of endangered species and ecosystems, and the preservation and enhancement of ecosystem services. We describe hurdles that the project faces, including data-sharing policies that ensure a permanent, freely available resource for future scientific discovery while respecting access and benefit sharing guidelines of the Nagoya Protocol. We also describe scientific and organizational challenges in executing such an ambitious project, and the structure proposed to achieve the project's goals. The far-reaching potential benefits of creating an open digital repository of genomic information for life on Earth can be realized only by a coordinated international effort.

Reconstruction and evolutionary history of eutherian chromosomes.

Whole-genome assemblies of 19 placental mammals and two outgroup species were used to reconstruct the order and orientation of syntenic fragments in chromosomes of the eutherian ancestor and six other descendant ancestors leading to human. For ancestral chromosome reconstructions, we developed an algorithm (DESCHRAMBLER) that probabilistically determines the adjacencies of syntenic fragments using chromosome-scale and fragmented genome assemblies. The reconstructed chromosomes of the eutherian, boreoeutherian, and euarchontoglires ancestor each included >80% of the entire length of the human genome, whereas reconstructed chromosomes of the most recent common ancestor of simians, catarrhini, great apes, and humans and chimpanzees included >90% of human genome sequence. These high-coverage reconstructions permitted reliable identification of chromosomal rearrangements over ∼105 My of eutherian evolution. Orangutan was found to have eight chromosomes that were completely conserved in homologous sequence order and orientation with the eutherian ancestor, the largest number for any species. Ruminant artiodactyls had the highest frequency of intrachromosomal rearrangements, and interchromosomal rearrangements dominated in murid rodents. A total of 162 chromosomal breakpoints in evolution of the eutherian ancestral genome to the human genome were identified; however, the rate of rearrangements was significantly lower (0.80/My) during the first ∼60 My of eutherian evolution, then increased to greater than 2.0/My along the five primate lineages studied. Our results significantly expand knowledge of eutherian genome evolution and will facilitate greater understanding of the role of chromosome rearrangements in adaptation, speciation, and the etiology of inherited and spontaneously occurring diseases.

Massive dysregulation of genes involved in cell signaling and placental development in cloned cattle conceptus and maternal endometrium.

A major unresolved issue in the cloning of mammals by somatic cell nuclear transfer (SCNT) is the mechanism by which the process fails after embryos are transferred to the uterus of recipients before or during the implantation window. We investigated this problem by using RNA sequencing (RNA-seq) to compare the transcriptomes in cattle conceptuses produced by SCNT and artificial insemination (AI) at day (d) 18 (preimplantation) and d 34 (postimplantation) of gestation. In addition, endometrium was profiled to identify the communication pathways that might be affected by the presence of a cloned conceptus, ultimately leading to mortality before or during the implantation window. At d 18, the effects on the transcriptome associated with SCNT were massive, involving more than 5,000 differentially expressed genes (DEGs). Among them are 121 genes that have embryonic lethal phenotypes in mice, cause defects in trophoblast and placental development, and/or affect conceptus survival in mice. In endometria at d 18, <0.4% of expressed genes were affected by the presence of a cloned conceptus, whereas at d 34, ∼36% and <0.7% of genes were differentially expressed in intercaruncular and caruncular tissues, respectively. Functional analysis of DEGs in placental and endometrial tissues suggests a major disruption of signaling between the cloned conceptus and the endometrium, particularly the intercaruncular tissue. Our results support a "bottleneck" model for cloned conceptus survival during the periimplantation period determined by gene expression levels in extraembryonic tissues and the endometrial response to altered signaling from clones.

The future of livestock breeding: genomic selection for efficiency, reduced emissions intensity, and adaptation.

As the global population and global wealth both continue to increase, so will the demand for livestock products, especially those that are highly nutritious. However, competition with other uses for land and water resources will also intensify, necessitating more efficient livestock production. In addition, as climate change escalates, reduced methane emissions from cattle and sheep will be a critical goal. Application of new technologies, including genomic selection and advanced reproductive technologies, will play an important role in meeting these challenges. Genomic selection, which enables prediction of the genetic merit of animals from genome-wide SNP markers, has already been adopted by dairy industries worldwide and is expected to double genetic gains for milk production and other traits. Here, we review these gains. We also discuss how the use of whole-genome sequence data should both accelerate the rate of gain and enable rapid discovery and elimination of genetic defects from livestock populations.

Reference-assisted chromosome assembly.

One of the most difficult problems in modern genomics is the assembly of full-length chromosomes using next generation sequencing (NGS) data. To address this problem, we developed "reference-assisted chromosome assembly" (RACA), an algorithm to reliably order and orient sequence scaffolds generated by NGS and assemblers into longer chromosomal fragments using comparative genome information and paired-end reads. Evaluation of results using simulated and real genome assemblies indicates that our approach can substantially improve genomes generated by a wide variety of de novo assemblers if a good reference assembly of a closely related species and outgroup genomes are available. We used RACA to reconstruct 60 Tibetan antelope (Pantholops hodgsonii) chromosome fragments from 1,434 SOAPdenovo sequence scaffolds, of which 16 chromosome fragments were homologous to complete cattle chromosomes. Experimental validation by PCR showed that predictions made by RACA are highly accurate. Our results indicate that RACA will significantly facilitate the study of chromosome evolution and genome rearrangements for the large number of genomes being sequenced by NGS that do not have a genetic or physical map.

Identification of a nonsense mutation in APAF1 that is likely causal for a decrease in reproductive efficiency in Holstein dairy cattle.

The HH1 haplotype on chromosome 5 is associated with a reduced conception rate and a deficit of homozygotes at the population level in Holstein cattle. The source HH1 haplotype was traced to the bull Pawnee Farm Arlinda Chief (Chief), who was born in 1962 and has sired more than 16,000 daughters. We identified a nonsense mutation in APAF1 (apoptotic protease activating factor 1;APAF1 p.Q579X) within HH1 using whole-genome resequencing of Chief and 3 of his sons. This mutation is predicted to truncate 670 AA (53.7%) of the encoded APAF1 protein that contains a WD40 domain critical to protein-protein interactions. Initial screening revealed no homozygous individuals for the mutation in 758 animals previously genotyped, whereas all 497 HH1 carriers possessed 1 copy of the mutant allele. Subsequent commercial genotyping of 246,773 Holsteins revealed 5,299 APAF1 heterozygotes and zero homozygotes for the mutation. The causative role of this mutation is also supported by functional data in mice that have demonstrated Apaf1 to be an essential molecule in the cytochrome-c-mediated apoptotic cascade and directly implicated in developmental and neurodegenerative disorders. In addition, most Apaf1 homozygous knockouts die by day 16.5 of development. We thus propose that the APAF1 p.Q579X nonsense mutation is the functional equivalent of the Apaf1 knockout. This mutation has caused an estimated 525,000 spontaneous abortions worldwide over the past 35 years, accounting for approximately $420 million in losses. With the mutation identified, selection against the deleterious allele in breeding schemes has aided in eliminating this defect from the population, reducing carrier frequency from 8% in past decades to 2% in 2015.

Nuclear reprogramming of cloned embryos and its implications for therapeutic cloning.

Therapeutic cloning, whereby somatic cell nuclear transfer (SCNT) is used to generate patient-specific embryonic stem cells (ESCs) from blastocysts cloned by nuclear transfer (ntESCs), holds great promise for the treatment of many human diseases. ntESCs have been derived in mice and cattle, but thus far there are no credible reports of human ntESCs. Here we review the recent literature on nuclear reprogramming by SCNT, including studies of gene expression, DNA methylation, chromatin remodeling, genomic imprinting and X chromosome inactivation. Reprogramming of genes expressed in the inner cell mass, from which ntESCs are derived, seems to be highly efficient. Defects in the extraembryonic lineage are probably the major cause of the low success rate of reproductive cloning but are not expected to affect the derivation of ntESCs. We remain optimistic that human therapeutic cloning is achievable and that the derivation of patient-specific ntESC lines will have great potential for regenerative medicine.

Whole-genome resequencing of two elite sires for the detection of haplotypes under selection in dairy cattle.

Using a combination of whole-genome resequencing and high-density genotyping arrays, genome-wide haplotypes were reconstructed for two of the most important bulls in the history of the dairy cattle industry, Pawnee Farm Arlinda Chief ("Chief") and his son Walkway Chief Mark ("Mark"), each accounting for ∼7% of all current genomes. We aligned 20.5 Gbp (∼7.3× coverage) and 37.9 Gbp (∼13.5× coverage) of the Chief and Mark genomic sequences, respectively. More than 1.3 million high-quality SNPs were detected in Chief and Mark sequences. The genome-wide haplotypes inherited by Mark from Chief were reconstructed using ∼1 million informative SNPs. Comparison of a set of 15,826 SNPs that overlapped in the sequence-based and BovineSNP50 SNPs showed the accuracy of the sequence-based haplotype reconstruction to be as high as 97%. By using the BovineSNP50 genotypes, the frequencies of Chief alleles on his two haplotypes then were determined in 1,149 of his descendants, and the distribution was compared with the frequencies that would be expected assuming no selection. We identified 49 chromosomal segments in which Chief alleles showed strong evidence of selection. Candidate polymorphisms for traits that have been under selection in the dairy cattle population then were identified by referencing Chief's DNA sequence within these selected chromosome blocks. Eleven candidate genes were identified with functions related to milk-production, fertility, and disease-resistance traits. These data demonstrate that haplotype reconstruction of an ancestral proband by whole-genome resequencing in combination with high-density SNP genotyping of descendants can be used for rapid, genome-wide identification of the ancestor's alleles that have been subjected to artificial selection.