Geo-evolutionary feedbacks: integrating rapid evolution and landscape change.

We develop a conceptual framework for geo-evolutionary feedbacks which describes the mutual interplay between landscape change and the evolution of traits of organisms residing on the landscape, with an emphasis on contemporary timeframes. Geo-evolutionary feedbacks can be realized via the direct evolution of geomorphic engineering traits or can be mediated by the evolution of trait variation that affects the population size and distribution of the specific geomorphic engineering organisms involved. Organisms that modify their local environments provide the basis for patch-scale geo-evolutionary feedbacks, whereas spatial self-organization provides a mechanism for geo-evolutionary feedbacks at the landscape scale. Understanding these likely prevalent geo-evolutionary feedbacks, that occur at timescales similar to anthropogenic climate change, will be essential to better predict landscape adaptive capacity and change.

Lack of Dosage Balance and Incomplete Dosage Compensation in the ZZ/ZW Gila Monster (Heloderma suspectum) Revealed by De Novo Genome Assembly.

Reptiles exhibit a variety of modes of sex determination, including both temperature-dependent and genetic mechanisms. Among those species with genetic sex determination, sex chromosomes of varying heterogamety (XX/XY and ZZ/ZW) have been observed with different degrees of differentiation. Karyotype studies have demonstrated that Gila monsters (Heloderma suspectum) have ZZ/ZW sex determination and this system is likely homologous to the ZZ/ZW system in the Komodo dragon (Varanus komodoensis), but little else is known about their sex chromosomes. Here, we report the assembly and analysis of the Gila monster genome. We generated a de novo draft genome assembly for a male using 10X Genomics technology. We further generated and analyzed short-read whole genome sequencing and whole transcriptome sequencing data for three males and three females. By comparing female and male genomic data, we identified four putative Z chromosome scaffolds. These putative Z chromosome scaffolds are homologous to Z-linked scaffolds identified in the Komodo dragon. Further, by analyzing RNAseq data, we observed evidence of incomplete dosage compensation between the Gila monster Z chromosome and autosomes and a lack of balance in Z-linked expression between the sexes. In particular, we observe lower expression of the Z in females (ZW) than males (ZZ) on a global basis, though we find evidence suggesting local gene-by-gene compensation. This pattern has been observed in most other ZZ/ZW systems studied to date and may represent a general pattern for female heterogamety in vertebrates.

Chromosome-Level Genome Assembly of the Blacktail Brush Lizard, Urosaurus nigricaudus, Reveals Dosage Compensation in an Endemic Lizard.

Urosaurus nigricaudus is a phrynosomatid lizard endemic to the Baja California Peninsula in Mexico. This work presents a chromosome-level genome assembly and annotation from a male individual. We used PacBio long reads and HiRise scaffolding to generate a high-quality genomic assembly of 1.87 Gb distributed in 327 scaffolds, with an N50 of 279 Mb and an L50 of 3. Approximately 98.4% of the genome is contained in 14 scaffolds, with 6 large scaffolds (334-127 Mb) representing macrochromosomes and 8 small scaffolds (63-22 Mb) representing microchromosomes. Using standard gene modeling and transcriptomic data, we predicted 17,902 protein-coding genes on the genome. The repeat content is characterized by a large proportion of long interspersed nuclear elements that are relatively old. Synteny analysis revealed some microchromosomes with high repeat content are more prone to rearrangements but that both macro- and microchromosomes are well conserved across reptiles. We identified scaffold 14 as the X chromosome. This microchromosome presents perfect dosage compensation where the single X of males has the same expression levels as two X chromosomes in females. Finally, we estimated the effective population size for U. nigricaudus was extremely low, which may reflect a reduction in polymorphism related to it becoming a peninsular endemic.

Identification of a novel polyomavirus in wild Sonoran Desert rodents of the family Heteromyidae.

Rodents are the largest and most diverse group of mammals. Covering a wide range of structural and functional adaptations, rodents successfully occupy virtually every terrestrial habitat, and they are often found in close association with humans, domestic animals, and wildlife. Although a significant amount of research has focused on rodents' prominence as known reservoirs of zoonotic viruses, there has been less emphasis on the viral ecology of rodents in general. Here, we utilized a viral metagenomics approach to investigate polyomaviruses in wild rodents from the Baja California peninsula, Mexico, using fecal samples. We identified a novel polyomavirus in fecal samples from two rodent species, a spiny pocket mouse (Chaetodipus spinatus) and a Dulzura kangaroo rat (Dipodomys simulans). These two polyomaviruses represent a new species in the genus Betapolyomavirus. Sequences of this polyomavirus cluster phylogenetically with those of other rodent polyomaviruses and two other non-rodent polyomaviruses (WU and KI) that have been identified in the human respiratory tract. Through our continued work on seven species of rodents, we endeavor to explore the viral diversity associated with wild rodents on the Baja California peninsula and expand on current knowledge of rodent viral ecology and evolution.

Social determinants of health and disease in companion dogs: a cohort study from the Dog Aging Project.

Exposure to social environmental adversity is associated with health and survival across many social species, including humans. However, little is known about how these health and mortality effects vary across the lifespan and may be differentially impacted by various components of the environment. Here, we leveraged a relatively new and powerful model for human aging, the companion dog, to investigate which components of the social environment are associated with dog health and how these associations vary across the lifespan. We drew on comprehensive survey data collected on 21,410 dogs from the Dog Aging Project and identified five factors that together explained 33.7% of the variation in a dog's social environment. Factors capturing financial and household adversity were associated with poorer health and lower physical mobility in companion dogs, while factors that captured social support, such as living with other dogs, were associated with better health when controlling for dog age and weight. Notably, the effects of each environmental component were not equal: the effect of social support was 5× stronger than financial factors. The strength of these associations depended on the age of the dog, including a stronger relationship between the owner's age and the dog's health in younger as compared to older dogs. Taken together, these findings suggest the importance of income, stability and owner's age on owner-reported health outcomes in companion dogs and point to potential behavioral and/or environmental modifiers that can be used to promote healthy aging across species.

Incomplete dosage balance and dosage compensation in the ZZ/ZW Gila monster (Heloderma suspectum) revealed by de novo genome assembly.

Reptiles exhibit a variety of modes of sex determination, including both temperature-dependent and genetic mechanisms. Among those species with genetic sex determination, sex chromosomes of varying heterogamety (XX/XY and ZZ/ZW) have been observed with different degrees of differentiation. Karyotype studies have demonstrated that Gila monsters (Heloderma suspectum) have ZZ/ZW sex determination and this system is likely homologous to the ZZ/ZW system in the Komodo dragon (Varanus komodoensis), but little else is known about their sex chromosomes. Here, we report the assembly and analysis of the Gila monster genome. We generated a de novo draft genome assembly for a male using 10X Genomics technology. We further generated and analyzed short-read whole genome sequencing and whole transcriptome sequencing data for three males and three females. By comparing female and male genomic data, we identified four putative Z-chromosome scaffolds. These putative Z-chromosome scaffolds are homologous to Z-linked scaffolds identified in the Komodo dragon. Further, by analyzing RNAseq data, we observed evidence of incomplete dosage compensation between the Gila monster Z chromosome and autosomes and a lack of balance in Z-linked expression between the sexes. In particular, we observe lower expression of the Z in females (ZW) than males (ZZ) on a global basis, though we find evidence suggesting local gene-by-gene compensation. This pattern has been observed in most other ZZ/ZW systems studied to date and may represent a general pattern for female heterogamety in vertebrates.

Integrating Earth-life systems: a geogenomic approach.

For centuries, scientists have recognized and worked to understand how Earth's mutable landscape and climate shape the distribution and evolution of species. Here, we describe the emerging field of geogenomics, which uses the reciprocal and deep integration of geologic, climatic, and population genomic data to define and test cause-effect relationships between Earth and life at intermediate spatial and temporal scales (i.e., the mesoscale). Technological advances now power the detailed reconstruction of landscape and evolutionary histories, but transdisciplinary collaborations and new quantitative tools are needed to better integrate Earth-life data. Geogenomics can help build a more unified theory and characterize the boundary conditions under which geologic and climatic processes generate new biodiversity, how species' responses differ, and why.

Towards a unified framework to study causality in Earth-life systems.

There is considerable interest in better understanding how earth processes shape the generation and distribution of life on Earth. This question, at its heart, is one of causation. In this article I propose that at a regional level, earth processes can be thought of as behaving somewhat deterministically and may have an organized effect on the diversification and distribution of species. However, the study of how landscape features shape biology is challenged by pseudocongruent or collinear variables. I demonstrate that causal structures can be used to depict the cause-effect relationships between earth processes and biological patterns using recent examples from the literature about speciation and species richness in montane settings. This application shows that causal diagrams can be used to better decipher the details of causal relationships by motivating new hypotheses. Additionally, the abstraction of this knowledge into structural equation metamodels can be used to formulate theory about relationships within Earth-life systems more broadly. Causal structures are a natural point of collaboration between biologists and Earth scientists, and their use can mitigate against the risk of misassigning causality within studies. My goal is that by applying causal theory through application of causal structures, we can build a systems-level understanding of what landscape features or earth processes most shape the distribution and diversification of species, what types of organisms are most affected, and why.

The hidden landscape: Evidence that sea-level change shaped the present population genomic patterns of marginal marine species.

Oftentimes, to understand the genetic relatedness and diversity of today's populations requires considering the ancient landscape on which those populations evolved. Nowhere is this clearer than along Earth's coastline, which has been in its present-day configuration for only about 6.5% of the past 800,000 years (Dolby et al., 2020; Miller et al., 2005). During ice ages when glaciers expanded in the Northern Hemisphere, they stored enough of the planet's water to drop global sea level by ~120 m below present levels ("lowstand", Figure 1a), and there have been at least eight of these 100,000-year cycles preceding today. When glaciers melted, ocean water reflooded shorelines, shifting and re-forming marginal marine habitats globally and shaping the relatedness of populations (Dolby et al., 2016). In a From the Cover article in this issue of Molecular Ecology, Stiller et al. (2020) integrate population genomic analysis of leafy seadragons in southern Australia with estimates of available seabed area to reveal that the expansion of habitat that accompanied this reflooding led to strong demographic expansions. With statistical models, they also show that western populations were eliminated and then recolonized because the continental shelf there is narrow, leaving little available habitat when sea level was low (Figure 1b). Their results document the dynamic and interrelated nature of a hidden, changing landscape and the evolution of species inhabiting it.

Immune and sex-biased gene expression in the threatened Mojave desert tortoise, Gopherus agassizii.

The immune system of ectotherms, particularly non-avian reptiles, remains poorly characterized regarding the genes involved in immune function, and their function in wild populations. We used RNA-Seq to explore the systemic response of Mojave desert tortoise (Gopherus agassizii) gene expression to three levels of Mycoplasma infection to better understand the host response to this bacterial pathogen. We found over an order of magnitude more genes differentially expressed between male and female tortoises (1,037 genes) than differentially expressed among immune groups (40 genes). There were 8 genes differentially expressed among both variables that can be considered sex-biased immune genes in this tortoise. Among experimental immune groups we find enriched GO biological processes for cysteine catabolism, regulation of type 1 interferon production, and regulation of cytokine production involved in immune response. Sex-biased transcription involves iron ion transport, iron ion homeostasis, and regulation of interferon-beta production to be enriched. More detailed work is needed to assess the seasonal response of the candidate genes found here. How seasonal fluctuation of testosterone and corticosterone modulate the immunosuppression of males and their susceptibility to Mycoplasma infection also warrants further investigation, as well as the importance of iron in the immune function and sex-biased differences of this species. Finally, future transcriptional studies should avoid drawing blood from tortoises via subcarapacial venipuncture as the variable aspiration of lymphatic fluid will confound the differential expression of genes.

Discovery of a New TLR Gene and Gene Expansion Event through Improved Desert Tortoise Genome Assembly with Chromosome-Scale Scaffolds.

Toll-like receptors (TLRs) are a complex family of innate immune genes that are well characterized in mammals and birds but less well understood in nonavian sauropsids (reptiles). The advent of highly contiguous draft genomes of nonmodel organisms enables study of such gene families through analysis of synteny and sequence identity. Here, we analyze TLR genes from the genomes of 22 tetrapod species. Findings reveal a TLR8 gene expansion in crocodilians and turtles (TLR8B), and a second duplication (TLR8C) specifically within turtles, followed by pseudogenization of that gene in the nonfreshwater species (desert tortoise and green sea turtle). Additionally, the Mojave desert tortoise (Gopherus agassizii) has a stop codon in TLR8B (TLR8-1) that is polymorphic among conspecifics. Revised orthology further reveals a new TLR homolog, TLR21-like, which is exclusive to lizards, snakes, turtles, and crocodilians. These analyses were made possible by a new draft genome assembly of the desert tortoise (gopAga2.0), which used chromatin-based assembly to yield draft chromosomal scaffolds (L50 = 26 scaffolds, N50 = 28.36 Mb, longest scaffold = 107 Mb) and an enhanced de novo genome annotation with 25,469 genes. Our three-step approach to orthology curation and comparative analysis of TLR genes shows what new insights are possible using genome assemblies with chromosome-scale scaffolds that permit integration of synteny conservation data.

Virus Discovery in Desert Tortoise Fecal Samples: Novel Circular Single-Stranded DNA Viruses.

The Sonoran Desert tortoise Gopherus morafkai is adapted to the desert, and plays an important ecological role in this environment. There is limited information on the viral diversity associated with tortoises (family Testudinidae), and to date no DNA virus has been identified associated with these animals. This study aimed to assess the diversity of DNA viruses associated with the Sonoran Desert tortoise by sampling their fecal matter. A viral metagenomics approach was used to identify the DNA viruses in fecal samples from wild Sonoran Desert tortoises in Arizona, USA. In total, 156 novel single-stranded DNA viruses were identified from 40 fecal samples. Those belonged to two known viral families, the Genomoviridae (n = 27) and Microviridae (n = 119). In addition, 10 genomes were recovered that belong to the unclassified group of circular-replication associated protein encoding single-stranded (CRESS) DNA virus and five circular molecules encoding viral-like proteins.

How sea level change mediates genetic divergence in coastal species across regions with varying tectonic and sediment processes.

Plate tectonics and sediment processes control regional continental shelf topography. We examine the genetic consequences of how glacial-associated sea level change interacted with variable nearshore topography since the last glaciation. We reconstructed the size and distribution of areas suitable for tidal estuary formation from the last glacial maximum, ~20 thousand years ago, to present from San Francisco, California, USA (~38°N) to Reforma, Sinaloa, Mexico (~25°N). We assessed range-wide genetic structure and diversity of three codistributed tidal estuarine fishes (California Killifish, Shadow Goby, Longjaw Mudsucker) along ~4,600 km using mitochondrial control region and cytB sequence, and 16-20 microsatellite loci from a total of 524 individuals. Results show that glacial-associated sea level change limited estuarine habitat to few, widely separated refugia at glacial lowstand, and present-day genetic clades were sourced from specific refugia. Habitat increased during postglacial sea level rise and refugial populations admixed in newly formed habitats. Continental shelves with active tectonics and/or low sediment supply were steep and hosted fewer, smaller refugia with more genetically differentiated populations than on broader shelves. Approximate Bayesian computation favoured the refuge-recolonization scenarios from habitat models over isolation by distance and seaway alternatives, indicating isolation at lowstand is a major diversification mechanism among these estuarine (and perhaps other) coastal species. Because sea level change is a global phenomenon, we suggest this top-down physical control of extirpation-isolation-recolonization may be an important driver of genetic diversification in coastal taxa inhabiting other topographically complex coasts globally during the Mid- to Late Pleistocene and deeper timescales.

The Agassiz's desert tortoise genome provides a resource for the conservation of a threatened species.

Agassiz's desert tortoise (Gopherus agassizii) is a long-lived species native to the Mojave Desert and is listed as threatened under the US Endangered Species Act. To aid conservation efforts for preserving the genetic diversity of this species, we generated a whole genome reference sequence with an annotation based on deep transcriptome sequences of adult skeletal muscle, lung, brain, and blood. The draft genome assembly for G. agassizii has a scaffold N50 length of 252 kbp and a total length of 2.4 Gbp. Genome annotation reveals 20,172 protein-coding genes in the G. agassizii assembly, and that gene structure is more similar to chicken than other turtles. We provide a series of comparative analyses demonstrating (1) that turtles are among the slowest-evolving genome-enabled reptiles, (2) amino acid changes in genes controlling desert tortoise traits such as shell development, longevity and osmoregulation, and (3) fixed variants across the Gopherus species complex in genes related to desert adaptations, including circadian rhythm and innate immune response. This G. agassizii genome reference and annotation is the first such resource for any tortoise, and will serve as a foundation for future analysis of the genetic basis of adaptations to the desert environment, allow for investigation into genomic factors affecting tortoise health, disease and longevity, and serve as a valuable resource for additional studies in this species complex.

Sea-level driven glacial-age refugia and post-glacial mixing on subtropical coasts, a palaeohabitat and genetic study.

Using a novel combination of palaeohabitat modelling and genetic mixture analyses, we identify and assess a sea-level-driven recolonization process following the Last Glacial Maximum (LGM). Our palaeohabitat modelling reveals dramatic changes in estuarine habitat distribution along the coast of California (USA) and Baja California (Mexico). At the LGM (approx. 20 kya), when sea level was approximately 130 m lower, the palaeo-shoreline was too steep for tidal estuarine habitat formation, eliminating this habitat type from regions where it is currently most abundant, and limiting such estuaries to a northern and a southern refugium separated by 1000 km. We assess the recolonization of estuaries formed during post-LGM sea-level rise through examination of refugium-associated alleles and approximate Bayesian computation in three species of estuarine fishes. Results reveal sourcing of modern populations from both refugia, which admix in the newly formed habitat between the refuges. We infer a dramatic peak in habitat area between 15 and 10 kya with subsequent decline. Overall, this approach revealed a previously undocumented dynamic and integrated relationship between sea-level change, coastal processes and population genetics. These results extend glacial refugial dynamics to unglaciated subtropical coasts and have significant implications for biotic response to predicted sea-level rise.