Search for: All records

Abstract Genomic resources across squamate reptiles (lizards and snakes) have lagged behind other vertebrate systems and high-quality reference genomes remain scarce. Of the 23 chromosome-scale reference genomes across the order, only 12 of the ~60 squamate families are represented. Within geckos (infraorder Gekkota), a species-rich clade of lizards, chromosome-level genomes are exceptionally sparse representing only two of the seven extant families. Using the latest advances in genome sequencing and assembly methods, we generated one of the highest-quality squamate genomes to date for the leopard gecko, Eublepharis macularius (Eublepharidae). We compared this assembly to the previous, short-read only, E. macularius reference genome published in 2016 and examined potential factors within the assembly influencing contiguity of genome assemblies using PacBio HiFi data. Briefly, the read N50 of the PacBio HiFi reads generated for this study was equal to the contig N50 of the previous E. macularius reference genome at 20.4 kilobases. The HiFi reads were assembled into a total of 132 contigs, which was further scaffolded using HiC data into 75 total sequences representing all 19 chromosomes. We identified 9 of the 19 chromosomal scaffolds were assembled as a near-single contig, whereas the other 10 chromosomes were each scaffolded together from multiple contigs. We qualitatively identified that the percent repeat content within a chromosome broadly affects its assembly contiguity prior to scaffolding. This genome assembly signifies a new age for squamate genomics where high-quality reference genomes rivaling some of the best vertebrate genome assemblies can be generated for a fraction of previous cost estimates. This new E. macularius reference assembly is available on NCBI at JAOPLA010000000. 

Synopsis When new land is created, initial microbial colonization lays the foundation for further ecological succession of plant and animal communities. Primary microbial succession of new aquatic habitats formed during volcanic activity has received little attention. The anchialine ecosystem, which includes coastal ponds in young lava flows, offers an opportunity to examine this process. Here, we characterized microbial communities of anchialine habitats in Hawaii that were created during volcanic eruptions in 2018. Benthic samples from three habitats were collected ∼2 years after their formation and at later time points spanning ∼1 year. Sequence profiling (16S and 18S) of prokaryotic and eukaryotic communities was used to test whether communities were similar to those from older, established anchialine habitats, and if community structure changed over time. Results show that microbial communities from the new habitats were unlike any from established anchialine microbial communities, having higher proportions of Planctomycetota and Chloroflexi but lower proportions of green algae. Each new habitat also harbored its own unique community relative to other habitats. While community composition in each habitat underwent statistically significant changes over time, they remained distinctive from established anchialine habitats. New habitats also had highly elevated temperatures compared to other habitats. These results suggest that idiosyncratic microbial consortia form during early succession of Hawaiian anchialine habitats. Future monitoring will reveal whether the early communities described here remain stable after temperatures decline and macro-organisms become more abundant, or if microbial communities will continue to change and eventually resemble those of established habitats. This work is a key first step in examining primary volcanic succession in aquatic habitats and suggests young anchialine habitats may warrant special conservation status. 

Abstract Rapid mutation rates are typical of mitochondrial genomes (mtDNAs) in animals, but it is not clear why. The difficulty of obtaining measurements of mtDNA mutation that are not biased by natural selection has stymied efforts to distinguish between competing hypotheses about the causes of high mtDNA mutation rates. Several studies which have measured mtDNA mutations in nematodes have yielded small datasets with conflicting conclusions about the relative abundance of different substitution classes (i.e., the mutation spectrum). We therefore leveraged Duplex Sequencing, a high-fidelity DNA sequencing technique, to characterize de novo mtDNA mutations in Caenorhabditis elegans. This approach detected nearly an order of magnitude more mtDNA mutations than documented in any previous nematode mutation study. Despite an existing extreme AT bias in the C. elegans mtDNA (75.6% AT), we found that a significant majority of mutations increase genomic AT content. Compared to some prior studies in nematodes and other animals, the mutation spectrum reported here contains an abundance of CG→AT transversions, supporting the hypothesis that oxidative damage may be a driver of mtDNA mutations in nematodes. Furthermore, we found an excess of G→T and C→T changes on the coding DNA strand relative to the template strand, consistent with increased exposure to oxidative damage. Analysis of the distribution of mutations across the mtDNA revealed significant variation among protein-coding genes and as well as among neighboring nucleotides. This high-resolution view of mitochondrial mutations in C. elegans highlights the value of this system for understanding relationships among oxidative damage, replication error, and mtDNA mutation. 

Synopsis Temperature is one of the most important environmental factors driving the genome-to-phenome relationship. Metabolic rates and related biological processes are predicted to increase with temperature due to the biophysical laws of chemical reactions. However, selection can also act on these processes across scales of biological organization, from individual enzymes to whole organisms. Although some studies have examined thermal responses across multiple scales, there is no general consensus on how these responses vary depending on the level of organization, or whether rates actually follow predicted theoretical patterns such as Arrhenius-like exponential responses or thermal performance curves (TPCs) that show peak responses. Here, we performed a meta-analysis on studies of ectotherms where biological rates were measured across the same set of temperatures, but at multiple levels of biological organization: enzyme activities, mitochondrial respiration, and/or whole-animal metabolic rates. Our final dataset consisted of 235 pairwise comparisons between levels of organization from 13 publications. Thermal responses differed drastically across levels of biological organization, sometimes showing completely opposite patterns. We developed a new effect size metric, “organizational disagreement” (OD) to quantify the difference in responses among levels of biological organization. Overall, rates at higher levels of biological organization (e.g., whole animal metabolic rates) increased more quickly with temperature than rates at lower levels, contrary to our predictions. Responses may differ across levels due to differing consequences of biochemical laws with increasing organization or due to selection for different responses. However, taxa and tissues examined generally did not affect OD. Theoretical TPCs, where rates increase to a peak value and then drop, were only rarely observed (12%), possibly because a broad range of test temperatures was rarely investigated. Exponential increases following Arrhenius predictions were more common (29%). This result suggests a classic assumption about thermal responses in biological rates is rarely observed in empirical datasets, although our results should be interpreted cautiously due to the lack of complete thermal profiles. We advocate for authors to explicitly address OD in their interpretations and to measure thermal responses across a wider, more incremental range of temperatures. These results further emphasize the complexity of connecting the genome to the phenome when environmental plasticity is incorporated: the impact of the environment on the phenotype can depend on the scale of organization considered. 

Abstract Microbiomes represent the collective bacteria, archaea, protist, fungi, and virus communities living in or on individual organisms that are typically multicellular eukaryotes. Such consortia have become recognized as having significant impacts on the development, health, and disease status of their hosts. Since understanding the mechanistic connections between an individual’s genetic makeup and their complete set of traits (i.e., genome to phenome) requires consideration at different levels of biological organization, this should include interactions with, and the organization of, microbial consortia. To understand microbial consortia organization, we elucidated the genetic constituents among phenotypically similar (and hypothesized functionally-analogous) layers (i.e., top orange, second orange, pink, and green layers) in the unique laminated orange cyanobacterial–bacterial crusts endemic to Hawaii’s anchialine ecosystem. High-throughput amplicon sequencing of ribosomal RNA hypervariable regions (i.e., Bacteria-specific V6 and Eukarya-biased V9) revealed microbial richness increasing by crust layer depth, with samples of a given layer more similar to different layers from the same geographic site than to their phenotypically-analogous layer from different sites. Furthermore, samples from sites on the same island were more similar to each other, regardless of which layer they originated from, than to analogous layers from another island. However, cyanobacterial and algal taxa were abundant in all surface and bottom layers, with anaerobic and chemoautotrophic taxa concentrated in the middle two layers, suggesting crust oxygenation from both above and below. Thus, the arrangement of oxygenated vs. anoxygenated niches in these orange crusts is functionally distinct relative to other laminated cyanobacterial–bacterial communities examined to date, with convergent evolution due to similar environmental conditions a likely driver for these phenotypically comparable but genetically distinct microbial consortia. 

Mitochondrial (mt) respiration depends on proteins encoded both by the mitochondrial and nuclear genomes. Variation in mt-DNA mutation rates exists across eukaryotes, although the functional consequences of elevated mt mutation rates in some lineages remain underexplored. In the angiosperm genus Silene , closely related, ecologically similar species have either ‘fast' or ‘slow' mt-DNA mutation rates. Here, we investigated the functional consequences of elevated mt-DNA mutation rates on mt respiration profiles of Silene mitochondria. Overall levels of respiration were similar among Species. Fast species had lower respiration efficiency than slow species and relied up to 48% more on nuclear-encoded respiratory enzymes alternative oxidase (AOX) and accessory dehydrogenases (DHex), which participate in stress responses in plants. However, not all fast species showed these trends. Respiratory profiles of some enzymes were correlated, most notably AOX and DHex. We conclude that subtle differences in mt physiology among Silene lineages with dramatically different mt mutation rates may underly similar phenotypes at higher levels of biological organization, betraying the consequences of mt mutations. 

Abstract The function and evolution of eukaryotic cells depend upon direct molecular interactions between gene products encoded in nuclear and cytoplasmic genomes. Understanding how these cytonuclear interactions drive molecular evolution and generate genetic incompatibilities between isolated populations and species is of central importance to eukaryotic biology. Plants are an outstanding system to investigate such effects because of their two different genomic compartments present in the cytoplasm (mitochondria and plastids) and the extensive resources detailing subcellular targeting of nuclear-encoded proteins. However, the field lacks a consistent classification scheme for mitochondrial- and plastid-targeted proteins based on their molecular interactions with cytoplasmic genomes and gene products, which hinders efforts to standardize and compare results across studies. Here, we take advantage of detailed knowledge about the model angiosperm Arabidopsis thaliana to provide a curated database of plant cytonuclear interactions at the molecular level. CyMIRA (Cytonuclear Molecular Interactions Reference for Arabidopsis) is available at http://cymira.colostate.edu/ and https://github.com/dbsloan/cymira and will serve as a resource to aid researchers in partitioning evolutionary genomic data into functional gene classes based on organelle targeting and direct molecular interaction with cytoplasmic genomes and gene products. It includes 11 categories (and 27 subcategories) of different cytonuclear complexes and types of molecular interactions, and it reports residue-level information for cytonuclear contact sites. We hope that this framework will make it easier to standardize, interpret, and compare studies testing the functional and evolutionary consequences of cytonuclear interactions.