Jagged 1 Rescues the Duchenne Muscular Dystrophy Phenotype.

Duchenne muscular dystrophy (DMD), caused by mutations at the dystrophin gene, is the most common form of muscular dystrophy. There is no cure for DMD and current therapeutic approaches to restore dystrophin expression are only partially effective. The absence of dystrophin in muscle results in dysregulation of signaling pathways, which could be targets for disease therapy and drug discovery. Previously, we identified two exceptional Golden Retriever muscular dystrophy (GRMD) dogs that are mildly affected, have functional muscle, and normal lifespan despite the complete absence of dystrophin. Now, our data on linkage, whole-genome sequencing, and transcriptome analyses of these dogs compared to severely affected GRMD and control animals reveals that increased expression of Jagged1 gene, a known regulator of the Notch signaling pathway, is a hallmark of the mild phenotype. Functional analyses demonstrate that Jagged1 overexpression ameliorates the dystrophic phenotype, suggesting that Jagged1 may represent a target for DMD therapy in a dystrophin-independent manner. PAPERCLIP.

Identifying recent adaptations in large-scale genomic data.

Although several hundred regions of the human genome harbor signals of positive natural selection, few of the relevant adaptive traits and variants have been elucidated. Using full-genome sequence variation from the 1000 Genomes (1000G) Project and the composite of multiple signals (CMS) test, we investigated 412 candidate signals and leveraged functional annotation, protein structure modeling, epigenetics, and association studies to identify and extensively annotate candidate causal variants. The resulting catalog provides a tractable list for experimental follow-up; it includes 35 high-scoring nonsynonymous variants, 59 variants associated with expression levels of a nearby coding gene or lincRNA, and numerous variants associated with susceptibility to infectious disease and other phenotypes. We experimentally characterized one candidate nonsynonymous variant in Toll-like receptor 5 (TLR5) and show that it leads to altered NF-κB signaling in response to bacterial flagellin. PAPERFLICK:

Progressive Cactus is a multiple-genome aligner for the thousand-genome era.

New genome assemblies have been arriving at a rapidly increasing pace, thanks to decreases in sequencing costs and improvements in third-generation sequencing technologies1-3. For example, the number of vertebrate genome assemblies currently in the NCBI (National Center for Biotechnology Information) database4 increased by more than 50% to 1,485 assemblies in the year from July 2018 to July 2019. In addition to this influx of assemblies from different species, new human de novo assemblies5 are being produced, which enable the analysis of not only small polymorphisms, but also complex, large-scale structural differences between human individuals and haplotypes. This coming era and its unprecedented amount of data offer the opportunity to uncover many insights into genome evolution but also present challenges in how to adapt current analysis methods to meet the increased scale. Cactus6, a reference-free multiple genome alignment program, has been shown to be highly accurate, but the existing implementation scales poorly with increasing numbers of genomes, and struggles in regions of highly duplicated sequences. Here we describe progressive extensions to Cactus to create Progressive Cactus, which enables the reference-free alignment of tens to thousands of large vertebrate genomes while maintaining high alignment quality. We describe results from an alignment of more than 600 amniote genomes, which is to our knowledge the largest multiple vertebrate genome alignment created so far.

Using evolutionary constraint to define novel candidate driver genes in medulloblastoma.

Current knowledge of cancer genomics remains biased against noncoding mutations. To systematically search for regulatory noncoding mutations, we assessed mutations in conserved positions in the genome under the assumption that these are more likely to be functional than mutations in positions with low conservation. To this end, we use whole-genome sequencing data from the International Cancer Genome Consortium and combined it with evolutionary constraint inferred from 240 mammals, to identify genes enriched in noncoding constraint mutations (NCCMs), mutations likely to be regulatory in nature. We compare medulloblastoma (MB), which is malignant, to pilocytic astrocytoma (PA), a primarily benign tumor, and find highly different NCCM frequencies between the two, in agreement with the fact that malignant cancers tend to have more mutations. In PA, a high NCCM frequency only affects the BRAF locus, which is the most commonly mutated gene in PA. In contrast, in MB, >500 genes have high levels of NCCMs. Intriguingly, several loci with NCCMs in MB are associated with different ages of onset, such as the HOXB cluster in young MB patients. In adult patients, NCCMs occurred in, e.g., the WASF-2/AHDC1/FGR locus. One of these NCCMs led to increased expression of the SRC kinase FGR and augmented responsiveness of MB cells to dasatinib, a SRC kinase inhibitor. Our analysis thus points to different molecular pathways in different patient groups. These newly identified putative candidate driver mutations may aid in patient stratification in MB and could be valuable for future selection of personalized treatment options.

Darwinian genomics and diversity in the tree of life.

Genomics encompasses the entire tree of life, both extinct and extant, and the evolutionary processes that shape this diversity. To date, genomic research has focused on humans, a small number of agricultural species, and established laboratory models. Fewer than 18,000 of ∼2,000,000 eukaryotic species (<1%) have a representative genome sequence in GenBank, and only a fraction of these have ancillary information on genome structure, genetic variation, gene expression, epigenetic modifications, and population diversity. This imbalance reflects a perception that human studies are paramount in disease research. Yet understanding how genomes work, and how genetic variation shapes phenotypes, requires a broad view that embraces the vast diversity of life. We have the technology to collect massive and exquisitely detailed datasets about the world, but expertise is siloed into distinct fields. A new approach, integrating comparative genomics with cell and evolutionary biology, ecology, archaeology, anthropology, and conservation biology, is essential for understanding and protecting ourselves and our world. Here, we describe potential for scientific discovery when comparative genomics works in close collaboration with a broad range of fields as well as the technical, scientific, and social constraints that must be addressed.

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.

Chiropterans Are a Hotspot for Horizontal Transfer of DNA Transposons in Mammalia.

Horizontal transfer of transposable elements (TEs) is an important mechanism contributing to genetic diversity and innovation. Bats (order Chiroptera) have repeatedly been shown to experience horizontal transfer of TEs at what appears to be a high rate compared with other mammals. We investigated the occurrence of horizontally transferred (HT) DNA transposons involving bats. We found over 200 putative HT elements within bats; 16 transposons were shared across distantly related mammalian clades, and 2 other elements were shared with a fish and two lizard species. Our results indicate that bats are a hotspot for horizontal transfer of DNA transposons. These events broadly coincide with the diversification of several bat clades, supporting the hypothesis that DNA transposon invasions have contributed to genetic diversification of bats.

Comparative neurogenetics of dog behavior complements efforts towards human neuropsychiatric genetics.

Domestic dogs display a wide array of heritable behaviors that have intermediate genetic complexity thanks to a long history of human-influenced selection. Comparative genetics in dogs could address the scarcity of non-human neurogenetic systems relevant to human neuropsychiatric disorders, which are characterized by mental, emotional, and behavioral symptoms and involve vastly complex genetic and non-genetic risk factors. Our review describes the diverse behavioral "phenome" of domestic dogs, past and ongoing sources of behavioral selection, and the state of canine behavioral genetics. We highlight two naturally disordered behavioral domains that illustrate how dogs may prove useful as a comparative, forward neurogenetic system: canine age-related cognitive dysfunction, which can be examined more rapidly given the attenuated lifespan of dogs, and compulsive disorders, which may have genetic roots in purpose-bred behaviors. Growing community science initiatives aimed at the companion dog population will be well suited to investigating such complex behavioral phenotypes and offer a comparative resource that parallels human genomic initiatives in scale and dimensionality.

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.

How to Make a Rodent Giant: Genomic Basis and Tradeoffs of Gigantism in the Capybara, the World's Largest Rodent.

Gigantism results when one lineage within a clade evolves extremely large body size relative to its small-bodied ancestors, a common phenomenon in animals. Theory predicts that the evolution of giants should be constrained by two tradeoffs. First, because body size is negatively correlated with population size, purifying selection is expected to be less efficient in species of large body size, leading to increased mutational load. Second, gigantism is achieved through generating a higher number of cells along with higher rates of cell proliferation, thus increasing the likelihood of cancer. To explore the genetic basis of gigantism in rodents and uncover genomic signatures of gigantism-related tradeoffs, we assembled a draft genome of the capybara (Hydrochoerus hydrochaeris), the world's largest living rodent. We found that the genome-wide ratio of nonsynonymous to synonymous mutations (ω) is elevated in the capybara relative to other rodents, likely caused by a generation-time effect and consistent with a nearly neutral model of molecular evolution. A genome-wide scan for adaptive protein evolution in the capybara highlighted several genes controlling postnatal bone growth regulation and musculoskeletal development, which are relevant to anatomical and developmental modifications for an increase in overall body size. Capybara-specific gene-family expansions included a putative novel anticancer adaptation that involves T-cell-mediated tumor suppression, offering a potential resolution to the increased cancer risk in this lineage. Our comparative genomic results uncovered the signature of an intragenomic conflict where the evolution of gigantism in the capybara involved selection on genes and pathways that are directly linked to cancer.

Solenodon genome reveals convergent evolution of venom in eulipotyphlan mammals.

Venom systems are key adaptations that have evolved throughout the tree of life and typically facilitate predation or defense. Despite venoms being model systems for studying a variety of evolutionary and physiological processes, many taxonomic groups remain understudied, including venomous mammals. Within the order Eulipotyphla, multiple shrew species and solenodons have oral venom systems. Despite morphological variation of their delivery systems, it remains unclear whether venom represents the ancestral state in this group or is the result of multiple independent origins. We investigated the origin and evolution of venom in eulipotyphlans by characterizing the venom system of the endangered Hispaniolan solenodon (Solenodon paradoxus). We constructed a genome to underpin proteomic identifications of solenodon venom toxins, before undertaking evolutionary analyses of those constituents, and functional assessments of the secreted venom. Our findings show that solenodon venom consists of multiple paralogous kallikrein 1 (KLK1) serine proteases, which cause hypotensive effects in vivo, and seem likely to have evolved to facilitate vertebrate prey capture. Comparative analyses provide convincing evidence that the oral venom systems of solenodons and shrews have evolved convergently, with the 4 independent origins of venom in eulipotyphlans outnumbering all other venom origins in mammals. We find that KLK1s have been independently coopted into the venom of shrews and solenodons following their divergence during the late Cretaceous, suggesting that evolutionary constraints may be acting on these genes. Consequently, our findings represent a striking example of convergent molecular evolution and demonstrate that distinct structural backgrounds can yield equivalent functions.

A comprehensive genomic history of extinct and living elephants.

Elephantids are the world's most iconic megafaunal family, yet there is no comprehensive genomic assessment of their relationships. We report a total of 14 genomes, including 2 from the American mastodon, which is an extinct elephantid relative, and 12 spanning all three extant and three extinct elephantid species including an ∼120,000-y-old straight-tusked elephant, a Columbian mammoth, and woolly mammoths. Earlier genetic studies modeled elephantid evolution via simple bifurcating trees, but here we show that interspecies hybridization has been a recurrent feature of elephantid evolution. We found that the genetic makeup of the straight-tusked elephant, previously placed as a sister group to African forest elephants based on lower coverage data, in fact comprises three major components. Most of the straight-tusked elephant's ancestry derives from a lineage related to the ancestor of African elephants while its remaining ancestry consists of a large contribution from a lineage related to forest elephants and another related to mammoths. Columbian and woolly mammoths also showed evidence of interbreeding, likely following a latitudinal cline across North America. While hybridization events have shaped elephantid history in profound ways, isolation also appears to have played an important role. Our data reveal nearly complete isolation between the ancestors of the African forest and savanna elephants for ∼500,000 y, providing compelling justification for the conservation of forest and savanna elephants as separate species.

Aquatic Adaptation and Depleted Diversity: A Deep Dive into the Genomes of the Sea Otter and Giant Otter.

Despite its recent invasion into the marine realm, the sea otter (Enhydra lutris) has evolved a suite of adaptations for life in cold coastal waters, including limb modifications and dense insulating fur. This uniquely dense coat led to the near-extinction of sea otters during the 18th-20th century fur trade and an extreme population bottleneck. We used the de novo genome of the southern sea otter (E. l. nereis) to reconstruct its evolutionary history, identify genes influencing aquatic adaptation, and detect signals of population bottlenecks. We compared the genome of the southern sea otter with the tropical freshwater-living giant otter (Pteronura brasiliensis) to assess common and divergent genomic trends between otter species, and with the closely related northern sea otter (E. l. kenyoni) to uncover population-level trends. We found signals of positive selection in genes related to aquatic adaptations, particularly limb development and polygenic selection on genes related to hair follicle development. We found extensive pseudogenization of olfactory receptor genes in both the sea otter and giant otter lineages, consistent with patterns of sensory gene loss in other aquatic mammals. At the population level, the southern sea otter and the northern sea otter showed extremely low genomic diversity, signals of recent inbreeding, and demographic histories marked by population declines. These declines may predate the fur trade and appear to have resulted in an increase in putatively deleterious variants that could impact the future recovery of the sea otter.

Natural selection and infectious disease in human populations.

The ancient biological 'arms race' between microbial pathogens and humans has shaped genetic variation in modern populations, and this has important implications for the growing field of medical genomics. As humans migrated throughout the world, populations encountered distinct pathogens, and natural selection increased the prevalence of alleles that are advantageous in the new ecosystems in both host and pathogens. This ancient history now influences human infectious disease susceptibility and microbiome homeostasis, and contributes to common diseases that show geographical disparities, such as autoimmune and metabolic disorders. Using new high-throughput technologies, analytical methods and expanding public data resources, the investigation of natural selection is leading to new insights into the function and dysfunction of human biology.

Nordic OCD & Related Disorders Consortium: Rationale, design, and methods.

Obsessive-compulsive disorder (OCD) is a debilitating psychiatric disorder, yet its etiology is unknown and treatment outcomes could be improved if biological targets could be identified. Unfortunately, genetic findings for OCD are lagging behind other psychiatric disorders. Thus, there is a pressing need to understand the causal mechanisms implicated in OCD in order to improve clinical outcomes and to reduce morbidity and societal costs. Specifically, there is a need for a large-scale, etiologically informative genetic study integrating genetic and environmental factors that presumably interact to cause the condition. The Nordic countries provide fertile ground for such a study, given their detailed population registers, national healthcare systems and active specialist clinics for OCD. We thus formed the Nordic OCD and Related Disorders Consortium (NORDiC, www.crowleylab.org/nordic), and with the support of NIMH and the Swedish Research Council, have begun to collect a large, richly phenotyped and genotyped sample of OCD cases. Our specific aims are geared toward answering a number of key questions regarding the biology, etiology, and treatment of OCD. This article describes and discusses the rationale, design, and methodology of NORDiC, including details on clinical measures and planned genomic analyses.

Low guanylyl cyclase activity in Weddell seals: implications for peripheral vasoconstriction and perfusion of the brain during diving.

Nitric oxide (NO) is a potent vasodilator, which improves perfusion and oxygen delivery during tissue hypoxia in terrestrial animals. The vertebrate dive response involves vasoconstriction in select tissues, which persists despite profound hypoxia. Using tissues collected from Weddell seals at necropsy, we investigated whether vasoconstriction is aided by downregulation of local hypoxia signaling mechanisms. We focused on NO-soluble guanylyl cyclase (GC)-cGMP signaling, a well-known vasodilatory transduction pathway. Seals have a lower GC protein abundance, activity, and capacity to respond to NO stimulation than do terrestrial mammals. In seal lung homogenates, GC produced less cGMP (20.1 ± 3.7 pmol·mg protein-1·min-1) than the lungs of dogs (-80 ± 144 pmol·mg protein-1·min-1 less than seals), sheep (-472 ± 96), rats (-664 ± 104) or mice (-1,160 ± 104, P < 0.0001). Amino acid sequences of the GC enzyme α-subunits differed between seals and terrestrial mammals, potentially affecting their structure and function. Vasoconstriction in diving Weddell seals is not consistent across tissues; perfusion is maintained in the brain and heart but decreased in other organs such as the kidney. A NO donor increased median GC activity 49.5-fold in the seal brain but only 27.4-fold in the kidney, consistent with the priority of cerebral perfusion during diving. Nos3 expression was high in the seal brain, which could improve NO production and vasodilatory potential. Conversely, Pde5a expression was high in the seal renal artery, which may increase cGMP breakdown and vasoconstriction in the kidney. Taken together, the results of this study suggest that alterations in the NO-cGMP pathway facilitate the diving response.

Analysis of the Human Mucosal Response to Cholera Reveals Sustained Activation of Innate Immune Signaling Pathways.

To better understand the innate immune response to Vibrio cholerae infection, we tracked gene expression in the duodenal mucosa of 11 Bangladeshi adults with cholera, using biopsy specimens obtained immediately after rehydration and 30 and 180 days later. We identified differentially expressed genes and performed an analysis to predict differentially regulated pathways and upstream regulators. During acute cholera, there was a broad increase in the expression of genes associated with innate immunity, including activation of the NF-κB, mitogen-activated protein kinase (MAPK), and Toll-like receptor (TLR)-mediated signaling pathways, which, unexpectedly, persisted even 30 days after infection. Focusing on early differences in gene expression, we identified 37 genes that were differentially expressed on days 2 and 30 across the 11 participants. These genes included the endosomal Toll-like receptor gene TLR8, which was expressed in lamina propria cells. Underscoring a potential role for endosomal TLR-mediated signaling in vivo, our pathway analysis found that interferon regulatory factor 7 and beta 1 and alpha 2 interferons were among the top upstream regulators activated during cholera. Among the innate immune effectors, we found that the gene for DUOX2, an NADPH oxidase involved in the maintenance of intestinal homeostasis, was upregulated in intestinal epithelial cells during cholera. Notably, the observed increases in DUOX2 and TLR8 expression were also modeled in vitro when Caco-2 or THP-1 cells, respectively, were stimulated with live V. cholerae but not with heat-killed organisms or cholera toxin alone. These previously unidentified features of the innate immune response to V. cholerae extend our understanding of the mucosal immune signaling pathways and effectors activated in vivo following cholera.