Search

Extracellular vesicles (EVs) are cell-secreted nanovesicles that play an important role in long-range cell-cell communication. Although EVs pose a promising alternative to cell-based therapy, targeted in vivo delivery still falls short. Many studies have explored the surface modification of EVs to enhance their targeting capabilities. However, to our knowledge, there are no standardized practices to confirm the successful surface modification of EVs or calculate the degree of conjugation on EV surfaces (conjugation efficiency). These pieces of information are essential in the reproducibility of targeted EV therapeutics and the determination of optimized conjugation conditions for EVs to see significant therapeutic effects in vitro and in vivo. This review will discuss the vast array of techniques adopted, technologies developed, and efficiency definitions made by studies that have calculated EV/nanoparticle surface conjugation efficiency and how differences between studies may contribute to differently reported conjugation efficiencies.

Range expansions, whether they are biological invasions or climate change-mediated range shifts, may have profound ecological and evolutionary consequences for plant-soil interactions. Range-expanding plants encounter soil biota with which they have a limited coevolutionary history, especially when introduced to a new continent. Past studies have found mixed results on whether plants experience positive or negative soil feedback interactions in their novel range, and these effects often change over time. One important theoretical explanation is that plants locally adapt to the soil pathogens and mutualists in their novel range. We tested this hypothesis in Dittrichia graveolens, an annual plant that is both expanding its European native range, initially coinciding with climate warming, and rapidly invading in California after human introduction. In parallel greenhouse experiments on both continents, we used plant genotypes and soils from five locations at the core and edge of each range to compare plant growth in soil inhabited by D. graveolens and nearby control microsites as a measure of plant-soil feedback. Plant-soil interactions were highly idiosyncratic across each range. On average, plant-soil feedbacks were more positive in the native range than in the exotic range. In line with the strongly heterogeneous pattern of soil responses along our biogeographic gradients, we found no evidence for evolutionary differentiation between plant genotypes from the core to edge of either range. Our results suggest that the evolution of plant-soil interactions during range expansion may be more strongly driven by local evolutionary dynamics varying across the range than by large-scale biogeographic shifts.

Genomic evidence supports an important role for selection in shaping patterns of introgression along the genome, but frameworks for understanding the evolutionary dynamics within hybrid populations that underlie these patterns have been lacking. Due to the clock-like effect of recombination in hybrids breaking up parental haplotypes, drift and selection produce predictable patterns of ancestry variation at varying spatial genomic scales through time. Here, we develop methods based on the Discrete Wavelet Transform to study the genomic scale of local ancestry variation and its association with recombination rates and show that these methods capture temporal dynamics of drift and genome-wide selection after hybridization. We apply these methods to published datasets from hybrid populations of swordtail fish (Xiphophorus) and baboons (Papio) and to inferred Neanderthal introgression in modern humans. Across systems, upward of 20% of variation in local ancestry at the broadest genomic scales can be attributed to systematic selection against introgressed alleles, consistent with strong selection acting on early-generation hybrids. Signatures of selection at fine genomic scales suggest selection over longer time scales; however, we suggest that our ability to confidently infer selection at fine scales is likely limited by inherent biases in current methods for estimating local ancestry from contiguous segments of genomic similarity. Wavelet approaches will become widely applicable as genomic data from systems with introgression become increasingly available and can help shed light on generalities of the genomic consequences of interspecific hybridization.

The entomopathogenic fungus Conoideocrella luteorostrata has recently been implicated in natural epizootics among exotic elongate hemlock scale (EHS) insects in Christmas tree farms in the eastern United States. Since 1913, C. luteorostrata has been reported from various plant feeding Hemiptera in the southeastern United States, but comprehensive morphological and phylogenetic studies of U.S. populations are lacking. The recovery of multiple strains of C. luteorostrata from mycosed EHS in North Carolina provided an opportunity to conduct pathogenicity assays and morphological and phylogenetic studies to investigate genus- and species-level boundaries among the Clavicipitaceae. Pathogenicity assays confirmed that C. luteorostrata causes mortality of EHS crawlers, an essential first step in developing this fungus as a biocontrol. Morphological studies revealed that conidia aligned with previous measurements of the Paecilomyces-like asexual state of C. luteorostrata, with conidiophore morphology consistent with historical observations. Additionally, a Hirsutella-like synanamorph was observed in select C. luteorostrata strains. In both a four-locus, 54-taxon Clavicipitaceae-wide phylogenetic analysis including D1-D2 domains of the nuclear 28S rRNA region (28S), elongation factor 1 alpha (EF1-α), DNA-directed RNA polymerase II subunit 1 (RPB1), and DNA-directed RNA polymerase II subunit 2 (RPB2) and a two-locus, 38-taxon (28S and EF1-α) phylogenetic analysis, all three Conoideocrella species were resolved as strongly supported monophyletic lineages across all loci and both methods (maximum likelihood and Bayesian inference) of phylogenetic inference except for 28S for C. tenuis. Despite the strong support for individual Conoideocrella species, none of the analyses supported the monophyly of Conoideocrella with the inclusion of Dussiella. Due to the paucity of RPB1 and RPB2 sequence data, EF1-α provided superior delimitation of intraspecies groupings for Conoideocrella and should be used in future studies. Further development of C. luteorostrata as a biocontrol against EHS will require additional surveys across diverse Hemiptera and expanded pathogenicity testing to clarify host range and efficacy of this fungus.

Successful alien species may experience a period of quiescence, known as the lag phase, before becoming invasive and widespread. The existence of lags introduces severe uncertainty in risk analyses of aliens as the present state of species is a poor predictor of future distributions, invasion success and impact. Predicting a species' ability to invade and pose negative impacts requires a quantitative understanding of the commonality and magnitude of lags, environmental factors and mechanisms likely to terminate lag. Using herbarium and climate data, we analysed over 5,700 time series (species × regions) in 3,505 naturalized plant species from nine regions in temperate and tropical climates to quantify lags and test whether there have been shifts in the species' climatic space during the transition from the lag phase to the expansion phase. Lags were identified in 35% of the assessed invasion events. We detected phylogenetic signals for lag phases in temperate climate regions and that annual self-fertilizing species were less likely to experience lags. Where lags existed, they had an average length of 40 years and a maximum of 320 years. Lengthy lags (>100 years) were more likely to occur in perennial plants and less frequent in self-pollinating species. For 98% of the species with a lag phase, the climate spaces sampled during the lag period differed from those in the expansion phase based on the assessment of centroid shifts or degree of climate space overlap. Our results highlight the importance of functional traits for the onset of the expansion phase and suggest that climate discovery may play a role in terminating the lag phase. However, other possibilities, such as sampling issues and climate niche shifts, cannot be ruled out.

Populations composed of individuals descended from multiple distinct genetic lineages often feature significant differences in phenotypic frequencies. We considered hatchery production of steelhead, the migratory anadromous form of the salmonid species Oncorhynchus mykiss, and investigated how differences among genetic lineages and environmental variation impacted life history traits. We genotyped 23,670 steelhead returning to the four California Central Valley hatcheries over 9 years from 2011 to 2019, confidently assigning parentage to 13,576 individuals to determine age and date of spawning and rates of iteroparity and repeat spawning within each year. We found steelhead from different genetic lineages showed significant differences in adult life history traits despite inhabiting similar environments. Differences between coastal and Central Valley steelhead lineages contributed to significant differences in age at return, timing of spawning, and rates of iteroparity among programs. In addition, adaptive genomic variation associated with life history development in this species varied among hatchery programs and was associated with the age of steelhead spawners only in the coastal lineage population. Environmental variation likely contributed to variations in phenotypic patterns observed over time, as our study period spanned both a marine heatwave and a serious drought in California. Our results highlight evidence of a strong genetic component underlying known phenotypic differences in life history traits between two steelhead lineages.

Our understanding on the interplay between gene functionality and gene arrangement at different chromosome scales relies on a few Diptera and the honeybee, species with quality-reference genome assemblies, accurate gene annotations, and abundant transcriptome data. Using recently generated 'omics resources in the monarch butterfly D. plexippus, a species with many more and smaller chromosomes relative to Drosophila species and the honeybee, we examined the organization of genes preferentially expressed at broadly defined developmental stages (larva, pupa, adult males, and adult females) at both fine and whole-chromosome scales. We found that developmental stage-regulated genes do not form more clusters, but do form larger clusters, than expected by chance, a pattern consistent across the gene categories examined. Notably, out of the thirty chromosomes in the monarch genome, twelve of them, plus the fraction of the chromosome Z that corresponds to the ancestral Z in other Lepidoptera, were found enriched for developmental stage-regulated genes. These two levels of nonrandom gene organization are not independent as enriched chromosomes for developmental stage-regulated genes tend to harbor disproportionately large clusters of these genes. Further, although paralogous genes were overrepresented in gene clusters, their presence is not enough to explain two-thirds of the documented cases of whole-chromosome enrichment. The composition of the largest clusters often included paralogs from more than one multigene family as well as unrelated single-copy genes. Our results reveal intriguing patterns at the whole-chromosome scale in D. plexippus while shedding light on the interplay between gene expression and chromosome organization beyond Diptera and Hymenoptera.

Many symbionts are sexually transmitted and impact their host's development, ecology, and evolution. While the significance of symbionts that cause sexually transmitted diseases (STDs) is relatively well understood, the prevalence and potential significance of the sexual transmission of mutualists remain elusive. Here, we study the effects of sexually transmitted mutualist nematodes on their dung beetle hosts. Symbiotic Diplogastrellus monhysteroides nematodes are present on the genitalia of male and female Onthophagus beetles and are horizontally transmitted during mating and vertically passed on to offspring during oviposition. A previous study indicates that the presence of nematodes benefits larval development and life history in a single host species, Onthophagus taurus. However, Diplogastrellus nematodes can be found in association with a variety of beetle species. Here, we replicate these previous experiments, assess whether the beneficial effects extend to other host species, and test whether nematode-mediated effects differ between male and female host beetles. Rearing three relatively distantly related dung beetle species with and without nematodes, we find that the presence of nematodes benefits body size, but not development time or survival across all three species. Likewise, we found no difference in the benefit of nematodes to male compared to female beetles. These findings highlight the role of sexually transmitted mutualists in the evolution and ecology of dung beetles.

Plant-systemic neonicotinoid (NN) insecticides can exert non-target impacts on organisms like beneficial insects and soil microbes. NNs can affect plant microbiomes, but we know little about their effects on microbial communities that mediate plant-insect interactions, including nectar-inhabiting microbes (NIMs). Here we employed two approaches to assess the impacts of NN exposure on several NIM taxa. First, we assayed the in vitro effects of six NN compounds on NIM growth using plate assays. Second, we inoculated a standardised NIM community into the nectar of NN-treated canola (Brassica napus) and assessed microbial survival and growth after 24 h. With few exceptions, in vitro NN exposure tended to decrease bacterial growth metrics. However, the magnitude of the decrease and the NN concentrations at which effects were observed varied substantially across bacteria. Yeasts showed no consistent in vitro response to NNs. In nectar, we saw no effects of NN treatment on NIM community metrics. Rather, NIM abundance and diversity responded to inherent plant qualities like nectar volume. In conclusion, we found no evidence that NIMs respond to field-relevant NN levels in nectar within 24 h, but our study suggests that context, specifically assay methods, time and plant traits, is important in assaying the effects of NNs on microbial communities.

Ecosystems that are coupled by reciprocal flows of energy and nutrient subsidies can be viewed as a single "meta-ecosystem." Despite these connections, the reciprocal flow of subsidies is greatly asymmetrical and seasonally pulsed. Here, we synthesize existing literature on stream-riparian meta-ecosystems to quantify global patterns of the amount of subsidy consumption by organisms, known as "allochthony." These resource flows are important since they can comprise a large portion of consumer diets, but can be disrupted by human modification of streams and riparian zones. Despite asymmetrical subsidy flows, we found stream and riparian consumer allochthony to be equivalent. Although both fish and stream invertebrates rely on seasonally pulsed allochthonous resources, we find allochthony varies seasonally only for fish, being nearly three times greater during the summer and fall than during the winter and spring. We also find that consumer allochthony varies with feeding traits for aquatic invertebrates, fish, and terrestrial arthropods, but not for terrestrial vertebrates. Finally, we find that allochthony varies by climate for aquatic invertebrates, being nearly twice as great in arid climates than in tropical climates, but not for fish. These findings are critical to understanding the consequences of global change, as ecosystem connections are being increasingly disrupted.

Phenology varies widely over space and time because of its sensitivity to climate. However, whether phenological variation is primarily generated by rapid organismal responses (plasticity) or local adaptation remains unresolved. Here we used 1,038,027 herbarium specimens representing 1,605 species from the continental United States to measure flowering-time sensitivity to temperature over time (Stime) and space (Sspace). By comparing these estimates, we inferred how adaptation and plasticity historically influenced phenology along temperature gradients and how their contributions vary among species with different phenology and native climates and among ecoregions differing in species composition. Parameters Sspace and Stime were positively correlated (r = 0.87), of similar magnitude and more frequently consistent with plasticity than adaptation. Apparent plasticity and adaptation generated earlier flowering in spring, limited responsiveness in late summer and delayed flowering in autumn in response to temperature increases. Nonetheless, ecoregions differed in the relative contributions of adaptation and plasticity, from consistently greater importance of plasticity (for example, southeastern United States plains) to their nearly equal importance throughout the season (for example, Western Sierra Madre Piedmont). Our results support the hypothesis that plasticity is the primary driver of flowering-time variation along temperature gradients, with local adaptation having a widespread but comparatively limited role.

Identifying and understanding behavioral resistance to insecticides is vital for maintaining global food security, public health, and ecological balance. Behavioral resistance has been documented to occur in a multitude of insect taxa dating back to the 1940s, but has not received significant research attention due primarily to the complexities of studying insect behavior and a lack of any clear definition of behavioral resistance. In recent years, a systematic effort to investigate the mechanism(s) of behavioral resistance in pest taxa (e.g. the German cockroach and the house fly) has been undertaken. Here, we practically define behavioral resistance, describe the efforts taken by research groups to elucidate resistance mechanisms, and provide insight on designing appropriate bioassays for investigating behavioral resistance mechanisms in the future.

Anthropogenic climate change has increased the frequency of drought, wildfire, and invasions of non-native species. Although high-severity fires linked to drought can inhibit recovery of native vegetation in forested ecosystems, it remains unclear how drought impacts the recovery of other plant communities following wildfire. We leveraged an existing rainfall manipulation experiment to test the hypothesis that reduced precipitation, fuel load, and fire severity convert plant community composition from native shrubs to invasive grasses in a Southern California coastal sage scrub system. We measured community composition before and after the 2020 Silverado wildfire in plots with three rainfall treatments. Drought reduced fuel load and vegetation cover, which reduced fire severity. Native shrubs had greater prefire cover in added water plots compared to reduced water plots. Native cover was lower and invasive cover was higher in postfire reduced water plots compared to postfire added and ambient water plots. Our results demonstrate the importance of fuel load on fire severity and plant community composition on an ecosystem scale. Management should focus on reducing fire frequency and removing invasive species to maintain the resilience of coastal sage scrub communities facing drought. In these communities, controlled burns are not recommended as they promote invasive plants.

Carpenter ants in the genus Camponotus are large, conspicuous ants that are abundant and ecologically influential in many terrestrial ecosystems. The bicolored carpenter ant, Camponotus vicinus Mayr, is distributed across a wide range of elevations and latitudes in western North America, where it is a prominent scavenger and predator. Here, we present a high-quality genome assembly of C. vicinus from a sample collected in Sonoma County, California, near the type locality of the species. This genome assembly consists of 38 scaffolds spanning 302.74 Mb, with contig N50 of 15.9 Mb, scaffold N50 of 19.9 Mb, and BUSCO completeness of 99.2%. This genome sequence will be a valuable resource for exploring the evolutionary ecology of C. vicinus and carpenter ants generally. It also provides an important tool for clarifying cryptic diversity within the C. vicinus species complex, a genetically diverse set of populations, some of which are quite localized and of conservation interest.

The little pocket mouse, Perognathus longimembris, and its nine congeners are small heteromyid rodents found in arid and seasonally arid regions of Western North America. The genus is characterized by behavioral and physiological adaptations to dry and often harsh environments, including nocturnality, seasonal torpor, food caching, enhanced osmoregulation, and a well-developed sense of hearing. Here we present a genome assembly of Perognathus longimembris longimembris generated from PacBio HiFi long read and Omni-C chromatin-proximity sequencing as part of the California Conservation Genomics Project. The assembly has a length of 2.35 Gb, contig N50 of 11.6 Mb, scaffold N50 of 73.2 Mb, and includes 93.8% of the BUSCO Glires genes. Interspersed repetitive elements constitute 41.2% of the genome. A comparison with the highly endangered Pacific pocket mouse, P. l. pacificus, reveals broad synteny. These new resources will enable studies of local adaptation, genetic diversity, and conservation of threatened taxa.

Smoky rubyspot damselflies (Hetaerina titia Drury, 1773) are one of the most commonly encountered odonates along streams and rivers on both slopes of Central America and the Atlantic drainages in the United States and southern Canada. Owing to their highly variable wing pigmentation, they have become a model system for studying sexual selection and interspecific behavioral interference. Here, we sequence and assemble the genome of a female smoky rubyspot. Of the primary assembly (i.e. the principle pseudohaplotype), 98.8% is made up of 12 chromosomal pseudomolecules (2N = 22A + X). There are 75 scaffolds in total, an N50 of 120 Mb, a contig-N50 of 0.64 Mb, and a high arthropod BUSCO score [C: 97.6% (S: 97.3%, D: 0.3%), F: 0.8%, M: 1.6%]. We then compare our assembly to that of the blue-tailed damselfly genome (Ischnura elegans), the most complete damselfly assembly to date, and a recently published assembly for an American rubyspot damselfly (Hetaerina americana). Collectively, these resources make Hetaerina a genome-enabled genus for further studies of the ecological and evolutionary forces shaping biological diversity.

The Yuma myotis bat (Myotis yumanensis) is a small vespertilionid bat and one of 52 species of new world Myotis bats in the subgenus Pizonyx. While M. yumanensis populations currently appear relatively stable, it is one of 12 bat species known or suspected to be susceptible to white-nose syndrome, the fungal disease causing declines in bat populations across North America. Only two of these 12 species have genome resources available, which limits the ability of resource managers to use genomic techniques to track the responses of bat populations to white-nose syndrome generally. Here we present the first de novo genome assembly for Yuma myotis, generated as a part of the California Conservation Genomics Project. The M. yumanensis genome was generated using a combination of PacBio HiFi long reads and Omni-C chromatin-proximity sequencing technology. This high-quality genome is one of the most complete bat assemblies available, with a contig N50 of 28.03 Mb, scaffold N50 of 99.14 Mb, and BUSCO completeness score of 93.7%. The Yuma myotis genome provides a high-quality resource that will aid in comparative genomic and evolutionary studies, as well as inform conservation management related to white-nose syndrome.

Across internally fertilising species, males transfer ejaculate proteins that trigger wide-ranging changes in female behaviour and physiology. Much theory has been developed to explore the drivers of ejaculate protein evolution. The accelerating availability of high-quality genomes now allows us to test how these proteins are evolving at fine taxonomic scales. Here, we use genomes from 264 species to chart the evolutionary history of Sex Peptide (SP), a potent regulator of female post-mating responses in Drosophila melanogaster. We infer that SP first evolved in the Drosophilinae subfamily and has since followed markedly different evolutionary trajectories in different lineages. Outside of the Sophophora-Lordiphosa, SP exists largely as a single-copy gene with independent losses in several lineages. Within the Sophophora-Lordiphosa, the SP gene family has repeatedly and independently expanded. Up to seven copies, collectively displaying extensive sequence variation, are present in some species. Despite these changes, SP expression remains restricted to the male reproductive tract. Alongside, we document considerable interspecific variation in the presence and morphology of seminal microcarriers that, despite the critical role SP plays in microcarrier assembly in D. melanogaster, appears to be independent of changes in the presence/absence or sequence of SP. We end by providing evidence that SP's evolution is decoupled from that of its receptor, Sex Peptide Receptor, in which we detect no evidence of correlated diversifying selection. Collectively, our work describes the divergent evolutionary trajectories that a novel gene has taken following its origin and finds a surprisingly weak coevolutionary signal between a supposedly sexually antagonistic protein and its receptor.

Despite rapid evolution across eutherian mammals, the X-linked MIR-506 family miRNAs are located in a region flanked by two highly conserved protein-coding genes (SLITRK2 and FMR1) on the X chromosome. Intriguingly, these miRNAs are predominantly expressed in the testis, suggesting a potential role in spermatogenesis and male fertility. Here, we report that the X-linked MIR-506 family miRNAs were derived from the MER91C DNA transposons. Selective inactivation of individual miRNAs or clusters caused no discernible defects, but simultaneous ablation of five clusters containing 19 members of the MIR-506 family led to reduced male fertility in mice. Despite normal sperm counts, motility, and morphology, the KO sperm were less competitive than wild-type sperm when subjected to a polyandrous mating scheme. Transcriptomic and bioinformatic analyses revealed that these X-linked MIR-506 family miRNAs, in addition to targeting a set of conserved genes, have more targets that are critical for spermatogenesis and embryonic development during evolution. Our data suggest that the MIR-506 family miRNAs function to enhance sperm competitiveness and reproductive fitness of the male by finetuning gene expression during spermatogenesis.

San Fernando Valley spineflower (Chorizanthe parryi var. fernandina [S. Watson] Jeps.) (Polygonaceae) is an herbaceous annual plant, endemic to California, and until rediscovered in 1999 had been thought to be extinct for almost seven decades. Historically documented at 10 locations, it currently persists at 2, separated by approximately 27 km. State listed as endangered, a description of its genetic diversity and structure is of conservation interest. After determining a lack of variation in ploidy, we examined genetic variation from samples within both populations: a common garden study for potentially adaptive genetic variation in selected growth and phenological traits and analysis of single nucleotide polymorphisms identified through restriction-site associated DNA sequencing. Both measures indicated that this highly restricted taxon nevertheless harbors substantial levels of genetic diversity and has significant between- and within-population genetic structure. Combining approaches from population genomics and common garden studies provided more insight into the patterns and basis of genetic diversity than is typical for studies of non-model species. Although local adaptation was not specifically studied (i.e., via reciprocal transplant studies), the differences determined from these two independent lines of evidence indicate that mixing gene pools between populations is not recommended at this time. Further, with significant differences revealed among subpopulations, we caution against mixing genotypes across subpopulations for the most part, and without much more evidence that this would not pose a risk of outbreeding depression. The importance of supporting pollinator health and diversity is highlighted. With genetic diversity—particularly with an annual species—being dynamic, fluctuating with the usual processes and with contributions from the soil seedbank, we recommend periodic resampling to monitor genetic diversity and structure. Climate change is anticipated to contribute to this variability.

Insectivorous birds can provide ecological and economic services by decreasing pest populations in agricultural systems, but the magnitude of effect is often debated. We modeled the capacity for birds to suppress pest population growth using a common tropical coffee pest, the coffee berry borer (CBB) (Hypothenemus hampei) as a study system. Previous field experiments show that birds play a role in suppressing CBB infestations through predation, but the degree to which birds can suppress population growth enough to control infestations is unknown since CBB are only vulnerable to predators when gravid females disperse. Using previously published data on CBB life-stage survivability, we constructed a female-only, daily time-step, deterministic Leslie matrix and projected CBB population growth for a single breeding season. Our goal was to assess the plausibility of CBB suppression by birds as a function of avian energy requirements, reported avian densities on coffee farms, prey composition of avian diets, estimated a caloric value of CBB, and the initial starting population size. Our model showed CBB population growth became exponential (λdaily = 1.042) and that at low, but not medium or high population sizes, were birds able to reduce population growth by 50%. In general, birds exert predation pressure on insect populations, but the ability to control infestations is complex, and is likely dependent on the initial CBB population size.

Understanding species distributions and predicting future range shifts requires considering all relevant abiotic factors and biotic interactions. Resource competition has received the most attention, but reproductive interference is another widespread biotic interaction that could influence species ranges. Rubyspot damselflies (Hetaerina spp.) exhibit a biogeographic pattern consistent with the hypothesis that reproductive interference has limited range expansion. Here, we use ecological niche models to evaluate whether this pattern could have instead been caused by niche differentiation. We found evidence for climatic niche differentiation, but the species that encounters the least reproductive interference has one of the narrowest and most peripheral niches. These findings strengthen the case that reproductive interference has limited range expansion and also provide a counterexample to the idea that release from negative species interactions triggers niche expansion. We propose that release from reproductive interference enables species to expand in range while specializing on the habitats most suitable for breeding.

BackgroundBiofilms in sulfide-rich springs present intricate microbial communities that play pivotal roles in biogeochemical cycling. We studied chemoautotrophically based biofilms that host diverse CPR bacteria and grow in sulfide-rich springs to investigate microbial controls on biogeochemical cycling.ResultsSulfide springs biofilms were investigated using bulk geochemical analysis, genome-resolved metagenomics, and scanning transmission X-ray microscopy (STXM) at room temperature and 87 K. Chemolithotrophic sulfur-oxidizing bacteria, including Thiothrix and Beggiatoa, dominate the biofilms, which also contain CPR Gracilibacteria, Absconditabacteria, Saccharibacteria, Peregrinibacteria, Berkelbacteria, Microgenomates, and Parcubacteria. STXM imaging revealed ultra-small cells near the surfaces of filamentous bacteria that may be CPR bacterial episymbionts. STXM and NEXAFS spectroscopy at carbon K and sulfur L2,3 edges show that filamentous bacteria contain protein-encapsulated spherical elemental sulfur granules, indicating that they are sulfur oxidizers, likely Thiothrix. Berkelbacteria and Moranbacteria in the same biofilm sample are predicted to have a novel electron bifurcating group 3b [NiFe]-hydrogenase, putatively a sulfhydrogenase, potentially linked to sulfur metabolism via redox cofactors. This complex could potentially contribute to symbioses, for example, with sulfur-oxidizing bacteria such as Thiothrix that is based on cryptic sulfur cycling. One Doudnabacteria genome encodes adjacent sulfur dioxygenase and rhodanese genes that may convert thiosulfate to sulfite. We find similar conserved genomic architecture associated with CPR bacteria from other sulfur-rich subsurface ecosystems.ConclusionsOur combined metagenomic, geochemical, spectromicroscopic, and structural bioinformatics analyses of biofilms growing in sulfide-rich springs revealed consortia that contain CPR bacteria and sulfur-oxidizing Proteobacteria, including Thiothrix, and bacteria from a new family within Beggiatoales. We infer roles for CPR bacteria in sulfur and hydrogen cycling. Video Abstract.

In this study, a Gram-stain-positive, non-motile, oxidase- and catalase-negative, rod-shaped, bacterial strain (SG_E_30_P1T) that formed light yellow colonies was isolated from a groundwater sample of Sztaravoda spring, Hungary. Based on 16S rRNA phylogenetic and phylogenomic analyses, the strain was found to form a distinct linage within the family Microbacteriaceae. Its closest relatives in terms of near full-length 16S rRNA gene sequences are Salinibacterium hongtaonis MH299814 (97.72 % sequence similarity) and Leifsonia psychrotolerans GQ406810 (97.57 %). The novel strain grows optimally at 20-28 °C, at neutral pH and in the presence of NaCl (1-2 w/v%). Strain SG_E_30_P1T contains MK-7 and B-type peptidoglycan with diaminobutyrate as the diagnostic amino acid. The major cellular fatty acids are anteiso-C15 : 0, iso-C16 : 0 and iso-C14 : 0, and the polar lipid profile is composed of diphosphatidylglycerol and phosphatidylglycerol, as well as an unidentified aminoglycolipid, aminophospholipid and some unidentified phospholipids. The assembled draft genome is a contig with a total length of 2 897 968 bp and a DNA G+C content of 65.5 mol%. Amino acid identity values with it closest relatives with sequenced genomes of

Marine algae are central to global carbon fixation, and their productivity is dictated largely by resource availability. Reduced nutrient availability is predicted for vast oceanic regions as an outcome of climate change; however, there is much to learn regarding response mechanisms of the tiny picoplankton that thrive in these environments, especially eukaryotic phytoplankton. Here, we investigate responses of the picoeukaryote Micromonas commoda, a green alga found throughout subtropical and tropical oceans. Under shifting phosphate availability scenarios, transcriptomic analyses revealed altered expression of transfer RNA modification enzymes and biased codon usage of transcripts more abundant during phosphate-limiting versus phosphate-replete conditions, consistent with the role of transfer RNA modifications in regulating codon recognition. To associate the observed shift in the expression of the transfer RNA modification enzyme complement with the transfer RNAs encoded by M. commoda, we also determined the transfer RNA repertoire of this alga revealing potential targets of the modification enzymes. Codon usage bias was particularly pronounced in transcripts encoding proteins with direct roles in managing phosphate limitation and photosystem-associated proteins that have ill-characterized putative functions in "light stress." The observed codon usage bias corresponds to a proposed stress response mechanism in which the interplay between stress-induced changes in transfer RNA modifications and skewed codon usage in certain essential response genes drives preferential translation of the encoded proteins. Collectively, we expose a potential underlying mechanism for achieving growth under enhanced nutrient limitation that extends beyond the catalog of up- or downregulated protein-encoding genes to the cell biological controls that underpin acclimation to changing environmental conditions.

The enzyme rubisco (ribulose-1,5-bisphosphate carboxylase/oxygenase) catalyzes the majority of biological carbon fixation on Earth. Although the vast majority of rubiscos across the tree of life assemble as homo-oligomers, the globally predominant form I enzyme-found in plants, algae, and cyanobacteria-forms a unique hetero-oligomeric complex. The recent discovery of a homo-oligomeric sister group to form I rubisco (named form I') has filled a key gap in our understanding of the enigmatic origins of the form I clade. However, to elucidate the series of molecular events leading to the evolution of form I rubisco, we must examine more distantly related sibling clades to contextualize the molecular features distinguishing form I and form I' rubiscos. Here, we present a comparative structural study retracing the evolutionary history of rubisco that reveals a complex structural trajectory leading to the ultimate hetero-oligomerization of the form I clade. We structurally characterize the oligomeric states of deep-branching form Iα and I'' rubiscos recently discovered from metagenomes, which represent key evolutionary intermediates preceding the form I clade. We further solve the structure of form I'' rubisco, revealing the molecular determinants that likely primed the enzyme core for the transition from a homo-oligomer to a hetero-oligomer. Our findings yield new insight into the evolutionary trajectory underpinning the adoption and entrenchment of the prevalent assembly of form I rubisco, providing additional context when viewing the enzyme family through the broader lens of protein evolution.

Life-history variation is the raw material of adaptation, and understanding its genetic and environmental underpinnings is key to designing effective conservation strategies. We used large-scale genetic pedigree reconstruction of anadromous steelhead trout (Oncorhynchus mykiss) from the Russian River, CA, USA, to elucidate sex-specific patterns of life-history traits and their heritability. SNP data from adults returning from sea over a 14-year period were used to identify 13,474 parent-offspring trios. These pedigrees were used to determine age structure, size distributions and family sizes for these fish, as well as to estimate the heritability of two key life-history traits, spawn date and age at maturity (first reproduction). Spawn date was highly heritable (h2  = 0.73) and had a cross-sex genetic correlation near unity. We provide the first estimate of heritability for age at maturity in ocean-going fish from this species and found it to be highly heritable (h2 from 0.29 to 0.62, depending on sex and method), with a much lower genetic correlation across sexes. We also evaluated genotypes at a migration-associated inversion polymorphism and found sex-specific correlations with age at maturity. The significant heritability of these two key reproductive traits in these imperiled fish, and their patterns of inheritance in the two sexes, is consistent with predictions of both natural and sexually antagonistic selection (sexes experience opposing selection pressures). This emphasizes the importance of anthropogenic factors, including hatchery practices and ecosystem modifications, in shaping the fitness of this species, thus providing important guidance for management and conservation efforts.

We explore how integrating behavioural ecology and macroecology can provide fundamental new insight into both fields, with particular relevance for understanding ecological responses to rapid environmental change. We outline the field of macrobehaviour, which aims to unite these disciplines explicitly, and highlight examples of research in this space. Macrobehaviour can be envisaged as a spectrum, where behavioural ecologists and macroecologists use new data and borrow tools and approaches from one another. At the heart of this spectrum, interdisciplinary research considers how selection in the context of large-scale factors can lead to systematic patterns in behavioural variation across space, time, and taxa, and in turn, influence macroecological patterns and processes. Macrobehaviour has the potential to enhance forecasts of future biodiversity change.

Migration is driven by a combination of environmental and genetic factors, but many questions remain about those drivers. Potential interactions between genetic and environmental variants associated with different migratory phenotypes are rarely the focus of study. We pair low coverage whole genome resequencing with a de novo genome assembly to examine population structure, inbreeding, and the environmental factors associated with genetic differentiation between migratory and resident breeding phenotypes in a species of conservation concern, the western burrowing owl (Athene cunicularia hypugaea). Our analyses reveal a dichotomy in gene flow depending on whether the population is resident or migratory, with the former being genetically structured and the latter exhibiting no signs of structure. Among resident populations, we observed significantly higher genetic differentiation, significant isolation-by-distance, and significantly elevated inbreeding. Among migratory breeding groups, on the other hand, we observed lower genetic differentiation, no isolation-by-distance, and substantially lower inbreeding. Using genotype-environment association analysis, we find significant evidence for relationships between migratory phenotypes (i.e., migrant versus resident) and environmental variation associated with cold temperatures during the winter and barren, open habitats. In the regions of the genome most differentiated between migrants and residents, we find significant enrichment for genes associated with the metabolism of fats. This may be linked to the increased pressure on migrants to process and store fats more efficiently in preparation for and during migration. Our results provide a significant contribution toward understanding the evolution of migratory behavior and vital insight into ongoing conservation and management efforts for the western burrowing owl.

Migratory animals follow seasonal cycles comprising linked phases often with different habitat requirements and demographic processes. Conservation of migratory species therefore must consider the full seasonal cycle to identify points limiting population viability. For western monarch butterflies, which have experienced significant declines, early spring is considered a critical period in the annual population cycle. However, records of western monarchs in early spring, when overall abundance is lowest, have historically been extremely limited. We used a community science initiative, the Western Monarch Mystery Challenge, to collect data on monarch distribution throughout the western United States between February 14th and April 22nd over 3 years. Using data from the Western Monarch Mystery Challenge and iNaturalist, we identified potential breeding habitat for western monarchs in early spring that spanned a large geographic area and several ecoregions. We observed monarchs in early spring that likely eclosed in the current year, suggesting that population expansion from overwintering sites reflects both movement and population growth. The number of records of western monarchs from early spring was higher during the Mystery Challenge (33.0/year) than earlier years (5.1/year). This study demonstrates the potential for and limitations of community science to increase our understanding of species at points in the life cycle when they are rare.

Large datasets along with sampling bias represent a challenge for phylodynamic reconstructions, particularly when the study data are obtained from various heterogeneous sources and/or through convenience sampling. In this study, we evaluate the presence of unbalanced sampled distribution by collection date, location, and risk group of human immunodeficiency virus Type 1 Subtype C using a comprehensive subsampling strategy and assess their impact on the reconstruction of the viral spatial and risk group dynamics using phylogenetic comparative methods. Our study shows that a most suitable dataset for ancestral trait reconstruction can be obtained through subsampling by all available traits, particularly using multigene datasets. We also demonstrate that sampling bias is inflated when considerable information for a given trait is unavailable or of poor quality, as we observed for the trait risk group. In conclusion, we suggest that, even if traits are not well recorded, including them deliberately optimizes the representativeness of the original dataset rather than completely excluding them. Therefore, we advise the inclusion of as many traits as possible with the aid of subsampling approaches in order to optimize the dataset for phylodynamic analysis while reducing the computational burden. This will benefit research communities investigating the evolutionary and spatio-temporal patterns of infectious diseases.

Background and scopeThe growth of experimental studies of crassulacean acid metabolism (CAM) in diverse plant clades, coupled with recent advances in molecular systematics, presents an opportunity to re-assess the phylogenetic distribution and diversity of species capable of CAM. It has been more than two decades since the last comprehensive lists of CAM taxa were published, and an updated survey of the occurrence and distribution of CAM taxa is needed to facilitate and guide future CAM research. We aimed to survey the phylogenetic distribution of these taxa, their diverse morphology, physiology and ecology, and the likely number of evolutionary origins of CAM based on currently known lineages.Results and conclusionsWe found direct evidence (in the form of experimental or field observations of gas exchange, day-night fluctuations in organic acids, carbon isotope ratios and enzymatic activity) for CAM in 370 genera of vascular plants, representing 38 families. Further assumptions about the frequency of CAM species in CAM clades and the distribution of CAM in the Cactaceae and Crassulaceae bring the currently estimated number of CAM-capable species to nearly 7 % of all vascular plants. The phylogenetic distribution of these taxa suggests a minimum of 66 independent origins of CAM in vascular plants, possibly with dozens more. To achieve further insight into CAM origins, there is a need for more extensive and systematic surveys of previously unstudied lineages, particularly in living material to identify low-level CAM activity, and for denser sampling to increase phylogenetic resolution in CAM-evolving clades. This should allow further progress in understanding the functional significance of this pathway by integration with studies on the evolution and genomics of CAM in its many forms.

We announce the assembly of the first de novo reference genome for the California Scrub-Jay (Aphelocoma californica). The genus Aphelocoma comprises four currently recognized species including many locally adapted populations across Mesoamerica and North America. Intensive study of Aphelocoma has revealed novel insights into the evolutionary mechanisms driving diversification in natural systems. Additional insights into the evolutionary history of this group will require continued development of high-quality, publicly available genomic resources. We extracted high molecular weight genomic DNA from a female California Scrub-Jay from northern California and generated PacBio HiFi long-read data and Omni-C chromatin conformation capture data. We used these data to generate a de novo partially phased diploid genome assembly, consisting of two pseudo-haplotypes, and scaffolded them using inferred physical proximity information from the Omni-C data. The more complete pseudo-haplotype assembly (arbitrarily designated "Haplotype 1") is 1.35 Gb in total length, highly contiguous (contig N50 = 11.53 Mb), and highly complete (BUSCO completeness score = 97%), with comparable scaffold sizes to chromosome-level avian reference genomes (scaffold N50 = 66.14 Mb). Our California Scrub-Jay assembly is highly syntenic with the New Caledonian Crow reference genome despite ~10 million years of divergence, highlighting the temporal stability of the avian genome. This high-quality reference genome represents a leap forward in publicly available genomic resources for Aphelocoma, and the family Corvidae more broadly. Future work using Aphelocoma as a model for understanding the evolutionary forces generating and maintaining biodiversity across phylogenetic scales can now benefit from a highly contiguous, in-group reference genome.

Rattlesnakes play important roles in their ecosystems by regulating prey populations, are involved in complex coevolutionary dynamics with their prey, and exhibit a variety of unusual adaptations, including maternal care, heat-sensing pit organs, hinged fangs, and medically-significant venoms. The western rattlesnake (Crotalus oreganus) is one of the widest ranging rattlesnake species, with a distribution from British Columbia, where it is listed as threatened, to Baja California and east across the Great Basin to western Wyoming, Colorado and New Mexico. Here, we report a new reference genome assembly for one of six currently recognized subspecies, C. oreganus helleri, as part of the California Conservation Genomics Project (CCGP). Consistent with the reference genomic sequencing strategy of the CCGP, we used Pacific Biosciences HiFi long reads and Hi-C chromatin-proximity sequencing technology to produce a de novo assembled genome. The assembly comprises a total of 698 scaffolds spanning 1,564,812,557 base pairs, has a contig N50 of 64.7 Mb, a scaffold N50 of 110.8 Mb, and BUSCO complete score of 90.5%. This reference genome will be valuable for studies on the genomic basis of venom evolution and variation within Crotalus, in resolving the taxonomy of C. oreganus and its relatives, and for the conservation and management of rattlesnakes in general.

Snakes in the family Colubridae include more than 2,000 currently recognized species, and comprise roughly 75% of the global snake species diversity on Earth. For such a spectacular radiation, colubrid snakes remain poorly understood ecologically and genetically. Two subfamilies, Colubrinae (788 species) and Dipsadinae (833 species), comprise the bulk of colubrid species richness. Dipsadines are a speciose and diverse group of snakes that largely inhabit Central and South America, with a handful of small-body-size genera that have invaded North America. Among them, the ring-necked snake, Diadophis punctatus, has an incredibly broad distribution with 14 subspecies. Given its continental distribution and high degree of variation in coloration, diet, feeding ecology, and behavior, the ring-necked snake is an excellent species for the study of genetic diversity and trait evolution. Within California, six subspecies form a continuously distributed "ring species" around the Central Valley, while a seventh, the regal ring-necked snake, Diadophis punctatus regalis is a disjunct outlier and Species of Special Concern in the state. Here, we report a new reference genome assembly for the San Diego ring-necked snake, D. p. similis, as part of the California Conservation Genomics Project. This assembly comprises a total of 444 scaffolds spanning 1,783 Mb and has a contig N50 of 8.0 Mb, scaffold N50 of 83 Mb, and BUSCO completeness score of 94.5%. This reference genome will be a valuable resource for studies of the taxonomy, conservation, and evolution of the ring-necked snake across its broad, continental distribution.

Acarospora socialis, the bright cobblestone lichen, is commonly found in southwestern North America. This charismatic yellow lichen is a species of key ecological significance as it is often a pioneer species in new environments. Despite their ecological importance virtually no research has been conducted on the genomics of A. socialis. To address this, we used long-read sequencing to generate the first high-quality draft genome of A. socialis. Lichen thallus tissue was collected from Pinkham Canyon in Joshua Tree National Park, California and deposited in the UC Riverside herbarium under accession #295874. The de novo assembly of the mycobiont partner of the lichen was generated from Pacific Biosciences HiFi long reads and Dovetail Omni-C chromatin capture data. After removing algal and bacterial contigs, the fungal genome was approximately 31.2 Mb consisting of 38 scaffolds with contig and scaffold N50 of 2.4 Mb. The BUSCO completeness score of the assembled genome was 97.5% using the Ascomycota gene set. Information on the genome of A. socialis is important for California conservation purposes given that this lichen is threatened in some places locally by wildfires due to climate change. This reference genome will be used for understanding the genetic diversity, population genomics, and comparative genomics of A. socialis species. Genomic resources for this species will support population and landscape genomics investigations, exploring the use of A. socialis as a bioindicator species for climate change, and in studies of adaptation by comparing populations that occur across aridity gradients in California.

The California Pipevine Swallowtail Butterfly, Battus philenor hirsuta, and its host plant, the California Pipevine or Dutchman's Pipe, Aristolochia californica Torr., are an important California endemic species pair. While this species pair is an ideal system to study co-evolution, genomic resources for both are lacking. Here, we report a new, chromosome-level assembly of B. philenor hirsuta as part of the California Conservation Genomics Project (CCGP). Following the sequencing and assembly strategy of the CCGP, we used Pacific Biosciences HiFi long reads and Hi-C chromatin proximity sequencing technology to produce a de novo assembled genome. Our genome assembly, the first for any species in the genus, contains 109 scaffolds spanning 443 mega base (Mb) pairs, with a contig N50 of 14.6 Mb, a scaffold N50 of 15.2 Mb, and BUSCO complete score of 98.9%. In combination with the forthcoming A. californica reference genome, the B. philenor hirsuta genome will be a powerful tool for documenting landscape genomic diversity and plant-insect co-evolution in a rapidly changing California landscape.

Diaphorina citri is the hemipteran pest and vector of a devastating bacterial pathogen of citrus worldwide. In addition to the two core bacterial endosymbionts of D. citri, Candidatus Carsonella ruddii and Candidatus Profftella armatura, the genome of a novel endosymbiont and as of yet undescribed microbe was discovered in a Hawaiian D. citri population through deep sequencing of multiple D. citri populations. Found to be closely related to the genus Asaia in the family Acetobacteraceae by 16S rRNA gene sequence analysis, it forms a sister clade along with other insect-associated 16S rRNA gene sequences from uncultured bacterium found associated with Aedes koreicus and Sogatella furcifera. Multilocus sequence analysis confirmed the phylogenetic placement sister to the Asaia clade. Despite the culturable Asaia clade being the closest phylogenetic neighbour, attempts to culture this newly identified bacterial endosymbiont were unsuccessful. On the basis of these distinct genetic differences, the novel endosymbiont is proposed to be classified into a candidate genus and species 'Candidatus Kirkpatrickella diaphorinae'. The full genome was deposited in GenBank (accession number CP107052; prokaryotic 16S rRNA OP600170).

Over 70% of soil organic carbon (SOC) is stored at a depth greater than 20 cm belowground. A portion of this deep SOC actively cycles on annual to decadal timescales and is sensitive to global change. However, deep SOC responses to global change likely differ from surface SOC responses because biotic controls on SOC cycling become weaker as mineral controls predominate with depth. Here, we synthesize the current information on deep SOC responses to the global change drivers of warming, shifting precipitation, elevated CO2, and land use and land cover change. Most deep SOC responses can only be hypothesized because few global change studies measure deep soils, and even fewer global change experiments manipulate deep soils. We call on scientists to incorporate deep soils into their manipulations, measurements, and models so that the response of deep SOC can be accounted for in projections of nature-based climate solutions and terrestrial feedbacks to climate change.

Many organisms actively manipulate the environment in ways that feed back on their own development, a process referred to as developmental niche construction. Yet, the role that constructed biotic and abiotic environments play in shaping phenotypic variation and its evolution is insufficiently understood. Here, we assess whether environmental modifications made by developing dung beetles impact the environment-sensitive expression of secondary sexual traits. Gazelle dung beetles both physically modify their ontogenetic environment and structure their biotic interactions through the vertical inheritance of microbial symbionts. By experimentally eliminating (i) physical environmental modifications and (ii) the vertical inheritance of microbes, we assess the degree to which (sym)biotic and physical environmental modifications shape the exaggeration of several traits varying in their degree and direction of sexual dimorphism. We expected the experimental reduction of a larva's ability to shape its environment to affect trait size and scaling, especially for traits that are sexually dimorphic and environmentally plastic. We find that compromised developmental niche construction indeed shapes sexual dimorphism in overall body size and the absolute sizes of male-limited exaggerated head horns, the strongly sexually dimorphic fore tibia length and width, as well as the weakly dimorphic elytron length and width. This suggests that environmental modifications affect sex-specific phenotypic variation in functional traits. However, most of these effects can be attributed to nutrition-dependent plasticity in size and non-isometric trait scaling rather than body-size-independent effects on the developmental regulation of trait size. Our findings suggest that the reciprocal relationship between developing organisms, their symbionts, and their environment can have considerable impacts on sexual dimorphism and functional morphology.

Phylogenetics has been foundational to SARS-CoV-2 research and public health policy, assisting in genomic surveillance, contact tracing, and assessing emergence and spread of new variants. However, phylogenetic analyses of SARS-CoV-2 have often relied on tools designed for de novo phylogenetic inference, in which all data are collected before any analysis is performed and the phylogeny is inferred once from scratch. SARS-CoV-2 data sets do not fit this mold. There are currently over 14 million sequenced SARS-CoV-2 genomes in online databases, with tens of thousands of new genomes added every day. Continuous data collection, combined with the public health relevance of SARS-CoV-2, invites an "online" approach to phylogenetics, in which new samples are added to existing phylogenetic trees every day. The extremely dense sampling of SARS-CoV-2 genomes also invites a comparison between likelihood and parsimony approaches to phylogenetic inference. Maximum likelihood (ML) and pseudo-ML methods may be more accurate when there are multiple changes at a single site on a single branch, but this accuracy comes at a large computational cost, and the dense sampling of SARS-CoV-2 genomes means that these instances will be extremely rare because each internal branch is expected to be extremely short. Therefore, it may be that approaches based on maximum parsimony (MP) are sufficiently accurate for reconstructing phylogenies of SARS-CoV-2, and their simplicity means that they can be applied to much larger data sets. Here, we evaluate the performance of de novo and online phylogenetic approaches, as well as ML, pseudo-ML, and MP frameworks for inferring large and dense SARS-CoV-2 phylogenies. Overall, we find that online phylogenetics produces similar phylogenetic trees to de novo analyses for SARS-CoV-2, and that MP optimization with UShER and matOptimize produces equivalent SARS-CoV-2 phylogenies to some of the most popular ML and pseudo-ML inference tools. MP optimization with UShER and matOptimize is thousands of times faster than presently available implementations of ML and online phylogenetics is faster than de novo inference. Our results therefore suggest that parsimony-based methods like UShER and matOptimize represent an accurate and more practical alternative to established ML implementations for large SARS-CoV-2 phylogenies and could be successfully applied to other similar data sets with particularly dense sampling and short branch lengths.

Nutrition-dependent growth of sexual traits is a major contributor to phenotypic diversity, and a large body of research documents insulin signalling as a major regulator of nutritional plasticity. However, findings across studies raise the possibility that the role of individual components within the insulin signalling pathway diverges in function among traits and taxa. Here, we use RNAi-mediated transcript depletion in the gazelle dung beetle to investigate the functions of forkhead box O (Foxo) and two paralogs of the insulin receptor (InR1 and InR2) in shaping nutritional plasticity in polyphenic male head horns, exaggerated fore legs, and weakly nutrition-responsive genitalia. Our functional genetic manipulations led to three main findings: FoxoRNAi reduced the length of exaggerated head horns in large males, while neither InR1 nor InR2 knock-downs resulted in measurable horn phenotypes. These results are similar to those documented previously for another dung beetle (Onthophagus taurus), but in stark contrast to findings in rhinoceros beetles. Secondly, knockdown of Foxo, InR1, and InR2 led to an increase in the intercept or slope of the scaling relationship of genitalia size. These findings are in contrast even to results documented previously for O. taurus. Lastly, while FoxoRNAi reduces male forelegs in D. gazella and O. taurus, the effects of InR1 and InR2 knockdowns diverged across dung beetle species. Our results add to the growing body of literature indicating that despite insulin signalling's conserved role as a regulator of nutritional plasticity, the functions of its components may diversify among traits and species, potentially fuelling the evolution of scaling relationships.

Human-driven environmental changes shape ecological communities from local to global scales. Within cities, landscape-scale patterns and processes and species characteristics generally drive local-scale wildlife diversity. However, cities differ in their structure, species pools, geographies and histories, calling into question the extent to which these drivers of wildlife diversity are predictive at continental scales. In partnership with the Urban Wildlife Information Network, we used occurrence data from 725 sites located across 20 North American cities and a multi-city, multi-species occupancy modelling approach to evaluate the effects of ecoregional characteristics and mammal species traits on the urbanization-diversity relationship. Among 37 native terrestrial mammal species, regional environmental characteristics and species traits influenced within-city effects of urbanization on species occupancy and community composition. Species occupancy and diversity were most negatively related to urbanization in the warmer, less vegetated cities. Additionally, larger-bodied species were most negatively impacted by urbanization across North America. Our results suggest that shifting climate conditions could worsen the effects of urbanization on native wildlife communities, such that conservation strategies should seek to mitigate the combined effects of a warming and urbanizing world.

Cancers rely on multiple, heterogeneous processes at different scales, pertaining to many biomedical fields. Therefore, understanding cancer is necessarily an interdisciplinary task that requires placing specialised experimental and clinical research into a broader conceptual, theoretical, and methodological framework. Without such a framework, oncology will collect piecemeal results, with scant dialogue between the different scientific communities studying cancer. We argue that one important way forward in service of a more successful dialogue is through greater integration of applied sciences (experimental and clinical) with conceptual and theoretical approaches, informed by philosophical methods. By way of illustration, we explore six central themes: (i) the role of mutations in cancer; (ii) the clonal evolution of cancer cells; (iii) the relationship between cancer and multicellularity; (iv) the tumour microenvironment; (v) the immune system; and (vi) stem cells. In each case, we examine open questions in the scientific literature through a philosophical methodology and show the benefit of such a synergy for the scientific and medical understanding of cancer.

Color signals which mediate behavioral interactions across taxa and contexts are often thought of as color 'patches' - parts of an animal that appear colorful compared to other parts of that animal. Color patches, however, cannot be considered in isolation because how a color is perceived depends on its visual background. This is of special relevance to the function and evolution of signals because backgrounds give rise to a fundamental tradeoff between color signal detectability and discriminability: as its contrast with the background increases, a color patch becomes more detectable, but discriminating variation in that color becomes more difficult. Thus, the signal function of color patches can only be fully understood by considering patch and background together as an integrated whole.

As climate changes, understanding the genetic basis of local adaptation in plants becomes an ever more pressing issue. Combining genotype-environment association (GEA) with genotype-phenotype association (GPA) analysis has an exciting potential to uncover the genetic basis of environmental responses. We use these approaches to identify genetic variants linked to local adaptation to drought in Pinus ponderosa. Over 4 million Single Nucleotide Polymorphisms (SNPs) were identified using 223 individuals from across the Sierra Nevada of California. 927,740 (22.3%) SNPs were retained after filtering for proximity to genes and used in our association analyses. We found 1374 associated with five major climate variables, with the largest number (1151) associated with April 1st snowpack. We also conducted a greenhouse study with various drought-tolerance traits measured in first-year seedlings of a subset of the genotyped trees grown in the greenhouse. 796 SNPs were associated with control-condition trait values, while 1149 were associated with responsiveness of these traits to drought. While no individual SNPs were associated with both the environmental variables and the measured traits, several annotated genes were associated with both, particularly those involved in cell wall formation, biotic and abiotic stress responses, and ubiquitination. However, the functions of many of the associated genes have not yet been determined due to the lack of gene annotation information for conifers. Future studies are needed to assess the developmental roles and ecological significance of these unknown genes.