Search NASA⌕ Search

SEARCH · Search NASA

Results for “VERTEBRATES”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Scavenging dynamics on Guam and implications for invasive species management

Deployment of mouse carcasses laced with acetaminophen has become a common management tool to control invasive brown tree snakes (Boiga irregularis; BTS) on Guam. Additionally, anticoagulant rodenticides may be used to control invasive rats (Rattus spp.) if their populations increase due to predator release in the wake of BTS eradication. However, there has been little research examining how scavengers on Guam could be incidentally exposed to toxicants by scavenging carcasses of animals that die from these population control strategies. Furthermore, there is a limited understanding of how the proliferation of invasive species on Guam has influenced the composition of the scavenger community. We investigated these topics by examining scavenger consumption of mouse, rat, and BTS carcasses on Guam in both a coastal and upland site during the wet (May–Aug 2016) and dry season (Jan–Apr 2017). We documented carcass consumption by 9 species, which scavenged 48% of carcasses. Interactions between season, habitat, and carcass type influenced probability of scavenging, and appeared to be driven by consumption by the two main scavenger species, BTS and cane toads (Rhinella marina), both of which are invasive on Guam. Baiting programs should consider the potential for toxin exposure to land crabs (Coenobita spp., Birgus latro), native species that scavenged at every combination of carcass type, habitat, and season. Altogether, 60% of scavenging events were attributed to species considered pests that are recent introductions to Guam. Invasive species on Guam are the primary scavengers of small vertebrate carrion, suggesting a substantial role in trophic dynamics that extends beyond predation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Metabolic skinflint or spendthrift? Insights into ground sloth integument and thermophysiology revealed by biophysical modeling and clumped isotope paleothermometry

Abstract Remains of megatheres have been known since the 18th -century and were among the first megafaunal vertebrates to be studied. While several examples of preserved integument show a thick coverage of fur for smaller ground sloths living in cold climates such as Mylodon and Nothrotheriops , comparatively very little is known about megathere skin. Assuming a typical placental mammal metabolism, it was previously hypothesized that megatheres would have had little-to-no fur as they achieved giant body sizes. Here the “hairless model of integument” is tested using geochemical analyses to estimate body temperature to generate novel models of ground sloth metabolism, fur coverage, and paleoclimate with Niche Mapper software. The simulations assuming metabolic activity akin to those of modern xenarthrans suggest that sparse fur coverage would have resulted in cold stress across most latitudinal ranges inhabited by extinct ground sloths. Specifically, Eremotherium predominantly required dense 10 mm fur with implications for seasonal changes of coat depth in northernmost latitudes and sparse fur in the tropics; Megatherium required dense 30 mm fur year-round in its exclusive range of cooler, drier climates; Mylodon and Nothrotheriops required dense 10–50 mm fur to avoid thermal stress, matching the integument remains of both genera, and further implying the use of behavioral thermoregulation. Moreover, clumped isotope paleothermometry data from the preserved teeth of four genera of ground sloth yielded reconstructed body temperatures lower than those previously reported for large terrestrial mammals (29 ± 2°–32 ± 3° C). This combination of low metabolisms and thick fur allowed ground sloths to inhabit various environments.

Deak, Michael D.↗
DOE OSTI2022

Contrasting modes of macro and microsynteny evolution in a eukaryotic subphylum

Examination of the changes in order and arrangement of homologous genes is key for understanding the mechanisms of genome evolution in eukaryotes. Previous comparisons between eukaryotic genomes have revealed considerable conservation across species that diverged hundreds of millions of years ago (e.g., vertebrates, bilaterian animals, and filamentous fungi). However, understanding how genome organization evolves within and between eukaryotic major lineages remains underexplored. We analyzed high-quality genomes of 120 representative budding yeast species (subphylum Saccharomycotina) spanning ~400 million years of eukaryotic evolution to examine how their genome organization evolved and to compare it with the evolution of animal and plant genome organization. We, in this study, found that the decay of both macrosynteny (the conservation of homologous chromosomes) and microsynteny (the conservation of local gene content and order) was strongly associated with evolutionary divergence across budding yeast major clades. However, although macrosynteny decayed very fast, within ~100 million years, the microsynteny of many genes—especially genes in metabolic clusters (e.g., in the GAL gene cluster)—was much more deeply conserved both within major clades and across the subphylum. We further found that when genomes with similar evolutionary divergence times were compared, budding yeasts had lower macrosynteny conservation than animals and filamentous fungi but higher conservation than angiosperms. In contrast, budding yeasts had levels of microsynteny conservation on par with mammals, whereas angiosperms exhibited very low conservation. Our results provide new insight into the tempo and mode of the evolution of gene and genome organization across an entire eukaryotic subphylum.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

A single oscillating proto-hypothalamic neuron gates taxis behavior in the primitive chordate Ciona

Ciona larvae display a number of behaviors, including negative phototaxis. In negative phototaxis, the larvae first perform short spontaneous rhythmic casting swims. As larvae cast in a light field, their photoreceptors are directionally shaded by an associated pigment cell, providing a phototactic cue. This then evokes an extended negative taxis swim. We report here that the larval forebrain of Ciona has a previously uncharacterized single slow-oscillating inhibitory neuron (neuron cor-assBVIN78) that projects to the midbrain, where it targets key interneurons of the phototaxis circuit known as the photoreceptor relay neurons. The anatomical location, gene expression and oscillation of cor-assBVIN78 suggest homology to oscillating neurons of the vertebrate hypothalamus. Ablation of cor-assBVIN78 results in larvae showing extended phototaxis-like swims, but which occur in the absence of phototactic cues. These results indicate that cor-assBVIN78 has a gating activity on phototaxis by projecting temporally-oscillating inhibition to the photoreceptor relay neurons. However, in intact larvae the frequency of cor-assBVIN78 oscillation does not match that of the rhythmic spontaneous swims, indicating that the troughs in oscillations do not themselves initiate swims, but rather that cor-assBVIN78 may modulate the phototaxis circuit by filtering out low level inputs while restricting them temporally to the troughs in inhibition.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Discovery of additional ancient genome duplications in yeasts

Whole-genome duplication (WGD) has had profound macroevolutionary impacts on diverse lineages, preceding adaptive radiations in vertebrates, teleost fish, and angiosperms. In contrast to the many known ancient WGDs in animals, and especially plants, we are aware of evidence for only four WGDs in fungi. The oldest of these occurred ∼100 million years ago (mya) and is shared by ∼60 extant Saccharomycetales species, including the baker’s yeast Saccharomyces cerevisiae. Notably, this is the only known ancient WGD event in the yeast subphylum Saccharomycotina. The dearth of ancient WGD events in fungi remains a mystery. Some studies have suggested that fungal lineages that experience chromosome and genome duplication quickly go extinct, leaving no trace in the genomic record, while others contend that the lack of known WGDs is due to an absence of data. Under the second hypothesis, additional sampling and deeper sequencing of fungal genomes should lead to the discovery of more WGD events. Coupling hundreds of recently published genomes from nearly every described Saccharomycotina species, with three additional long-read assemblies, we discovered three novel WGD events. Although the functions of retained duplicate genes originating from these events are broad, they bear similarities to the well-known WGD that occurred in the Saccharomycetales. In conclusion, our results suggest that WGD may be a more common evolutionary force in fungi than previously believed.

convergent evolution↗
DOE OSTI2020

Uncovering DNA-PKcs ancient phylogeny, unique sequence motifs and insights for human disease

DNA-dependent protein kinase catalytic subunit (DNA-PKcs) is a key member of the phosphatidylinositol-3 kinase-like (PIKK) family of protein kinases with critical roles in DNA-double strand break repair, transcription, metastasis, mitosis, RNA processing, and innate and adaptive immunity. The absence of DNA-PKcs from many model organisms has led to the assumption that DNA-PKcs is a vertebrate-specific PIKK. Here, we find that DNA-PKcs is widely distributed in invertebrates, fungi, plants, and protists, and that threonines 2609, 2638, and 2647 of the ABCDE cluster of phosphorylation sites are highly conserved amongst most Eukaryotes. Furthermore, we identify highly conserved amino acid sequence motifs and domains that are characteristic of DNA-PKcs relative to other PIKKs. These include residues in the Forehead domain and a novel motif we have termed YRPD, located in an α helix C-terminal to the ABCDE phosphorylation site loop. Combining sequence with biochemistry plus structural data on human DNA-PKcs unveils conserved sequence and conformational features with functional insights and implications. The defined generally progressive DNA-PKcs sequence diversification uncovers conserved functionality supported by Evolutionary Trace analysis, suggesting that for many organisms both functional sites and evolutionary pressures remain identical due to fundamental cell biology. The mining of cancer genomic data and germline mutations causing human inherited disease reveal that robust DNA-PKcs activity in tumors is detrimental to patient survival, whereas germline mutations compromising function are linked to severe immunodeficiency and neuronal degeneration. We anticipate that these collective results will enable ongoing DNA-PKcs functional analyses with biological and medical implications.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Relationships between maternally-transferred mercury and hatchling development, behavior, and survival in the American alligator ( Alligator mississippiensis )

Mercury is a toxic and pervasive environmental contaminant that can be transferred from mother to offspring during development. Consequences of maternally-transferred mercury have been observed in vertebrate taxa, including reduced clutch viability, reduced offspring size, and behavioral alterations. These sublethal effects have been assumed to decrease survivorship, though this is seldom assessed. Here, we examined how maternally-transferred mercury interacts with incubation temperature to influence reproductive success, offspring behavior, and subsequent survival in the American alligator (Alligator mississippiensis). We collected nine clutches of eggs from a mercury contaminated reservoir on the Savannah River Site, South Carolina, and incubated eggs at either female- or male-promoting temperatures. Clutch-averaged mercury in egg yolk was high relative to other studies in crocodilians and ranged from 0.248 to 0.554 ppm compared to 0.018–0.052 ppm at a site with low levels of mercury contamination; mercury levels in hatchling blood ranged from 0.090 to 0.490 ppm ($\overline{x}$= 0.240 ppm, n = 158). We found few, mostly negligible correlations between life history traits and mercury but noted a positive relationship with egg mass, possibly mediated by correlated maternal effects such as resource provisioning. Incubation temperature exerted strong effects on hatchling phenotypes, with warmer, male-promoting temperatures producing larger and bolder hatchlings. Presumptive females, produced from cooler incubation temperatures, spent more time in warm areas during behavior trials. Hatchlings were released 10–15 days post-hatch and surveyed over eight months to assess survival. Survivorship was positively correlated with hatchling size and negatively correlated with proportional time spent in warm areas. Presumptive females had much lower survival, and overall survivorship for the eight-month period was 0.185–0.208, depending on the modelling approach. Finally, our study suggests that, within the range of concentrations we observed, incubation temperature has a stronger effect on offspring behavior and survival than maternally-transferred mercury pollution in American alligators.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Vessel Rupture Thresholds for Vessel–Bubble Interactions Using an Earthworm Vasculature Model

Intravenous microbubble oscillation in the presence of ultrasound has the potential to yield a wide range of therapeutic benefits. However, the likelihood of vessel damage caused by mechanical effects has not been quantified as a function of the numerous important parameters in therapeutic ultrasound procedures. In this study, we examined the effects of microbubbles injected into the vasculature of the earthworm. It was found that the elastic properties of earthworm blood vessels are similar to those of arteries in older humans, and that earthworms are well suited to the large number of experiments necessary to investigate safety of procedures involving microbubble oscillation in sonicated vessels. Microbubbles were infused into earthworm vessels, and the rupture time during sonication was recorded as a function of ultrasound frequency, pulse repetition frequency and acoustic pressure. A modified mechanical index (MMI) was defined that successfully captured the trends in rupture probability and rupture time for the different parameter values, creating a database of vessel rupture thresholds. In the absence of bubbles, the product of MMI squared and rupture time was approximately constant, indicating a possible radiation-force effect. Here, the MMI was an effective correlating parameter in the presence of bubbles, though the mathematical dependence is not yet apparent. The results of the study are expected to be valuable in designing more refined studies in vertebrate models, as well as informing computational models.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Implication of rare genetic variants of NODAL and ACVR1B in congenital heart disease patients from Indian population

Highlights: • Screening of 300 CHD patients revealed two rare missense variants (p.P51L & p.D327E) in NODAL. • A novel missense variation p.M345I was identified in ACVR1B. • Bioinformatic analyses suggested that these variants could have deleterious role. • Three-dimensional modelling suggested conformation changes in the mutated proteins compared to wild-type. • In vitro studies demonstrated that missense variants affect the function of NODAL and ACVR1B. NODAL signaling plays an essential role in vertebrate embryonic patterning and heart development. Accumulating evidences suggest that genetic mutations in TGF-β/NODAL signaling pathway can cause congenital heart disease in humans. To investigate the implication of NODAL signaling in isolated cardiovascular malformation, we have screened 300 non-syndromic CHD cases and 200 controls for NODAL and ACVR1B by Sanger sequencing and identified two rare missense (c.152C > T; p.P51L and c.981 T > A; p.D327E) variants in NODAL and a novel missense variant c.1035G > A; p.M345I in ACVR1B. All these variants are absent in 200 controls. Three-dimensional protein-modelling demonstrates that both p.P51L and p.D327E variations of NODAL and p.M345I mutation of ACVR1B, affect the tertiary structure of respective proteins. Variants of NODAL (p.P51L and p.D327E) and ACVR1B (p.M345I), significantly reduce the transactivation of AR3-Luc, (CAGA){sub 12}-Luc and (SBE){sub 4}-Luc promoters. Moreover, qRT-PCR results have also deciphered a reduction in the expression of cardiac-enriched transcription factors namely Gata4, Nkx2-5, and Tbx5 in both the mutants of NODAL. Decreased expression of, Gata4, Nkx2-5, Tbx5, and lefty is observed in p.M345I mutant of ACVR1B as well. Additionally, reduced phosphorylation of SMAD2/3 in response to these variants, suggests impaired NODAL signaling and possibly responsible for defective cell fate decision and differentiation of cardiomyocytes leading to CHD phenotype.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

TASOR is a pseudo-PARP that directs HUSH complex assembly and epigenetic transposon control

The HUSH complex represses retroviruses, transposons and genes to maintain the integrity of vertebrate genomes. HUSH regulates deposition of the epigenetic mark H3K9me3, but how its three core subunits — TASOR, MPP8 and Periphilin — contribute to assembly and targeting of the complex remains unknown. Here, we define the biochemical basis of HUSH assembly and find that its modular architecture resembles the yeast RNA-induced transcriptional silencing complex. TASOR, the central HUSH subunit, associates with RNA processing components. TASOR is required for H3K9me3 deposition over LINE-1 repeats and repetitive exons in transcribed genes. In the context of previous studies, this suggests that an RNA intermediate is important for HUSH activity. We dissect the TASOR and MPP8 domains necessary for transgene repression. Structure-function analyses reveal TASOR bears a catalytically-inactive PARP domain necessary for targeted H3K9me3 deposition. We conclude that TASOR is a multifunctional pseudo-PARP that directs HUSH assembly and epigenetic regulation of repetitive genomic targets.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

A spatio-temporally constrained gene regulatory network directed by PBX1/2 acquires limb patterning specificity via HAND2

A lingering question in developmental biology has centered on how transcription factors with widespread distribution in vertebrate embryos can perform tissue-specific functions. Here, using the murine hindlimb as a model, we investigate the elusive mechanisms whereby PBX TALE homeoproteins, viewed primarily as HOX cofactors, attain context-specific developmental roles despite ubiquitous presence in the embryo. We first demonstrate that mesenchymal-specific loss of PBX1/2 or the transcriptional regulator HAND2 generates similar limb phenotypes. By combining tissue-specific and temporally controlled mutagenesis with multi-omics approaches, we reconstruct a gene regulatory network (GRN) at organismal-level resolution that is collaboratively directed by PBX1/2 and HAND2 interactions in subsets of posterior hindlimb mesenchymal cells. Genome-wide profiling of PBX1 binding across multiple embryonic tissues further reveals that HAND2 interacts with subsets of PBX-bound regions to regulate limb-specific GRNs. Our research elucidates fundamental principles by which promiscuous transcription factors cooperate with cofactors that display domain-restricted localization to instruct tissue-specific developmental programs.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Structural basis for the allosteric modulation of rhodopsin by nanobody binding to its extracellular domain

Rhodopsin is a prototypical G protein-coupled receptor (GPCR) critical for vertebrate vision. Research on GPCR signaling states has been facilitated using llama-derived nanobodies (Nbs), some of which bind to the intracellular surface to allosterically modulate the receptor. Extracellularly binding allosteric nanobodies have also been investigated, but the structural basis for their activity has not been resolved to date. Here, we report a library of Nbs that bind to the extracellular surface of rhodopsin and allosterically modulate the thermodynamics of its activation process. Crystal structures of Nb2 in complex with native rhodopsin reveal a mechanism of allosteric modulation involving extracellular loop 2 and native glycans. Nb2 binding suppresses Schiff base deprotonation and hydrolysis and prevents intracellular outward movement of helices five and six – a universal activation event for GPCRs. Nb2 also mitigates protein misfolding in a disease-associated mutant rhodopsin. Our data show the power of nanobodies to modulate the photoactivation of rhodopsin and potentially serve as therapeutic agents for disease-associated rhodopsin misfolding.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Conserved chromatin and repetitive patterns reveal slow genome evolution in frogs

Frogs are an ecologically diverse and phylogenetically ancient group of anuran amphibians that include important vertebrate cell and developmental model systems, notably the genus Xenopus. Here we report a high-quality reference genome sequence for the western clawed frog, Xenopus tropicalis, along with draft chromosome-scale sequences of three distantly related emerging model frog species, Eleutherodactylus coqui, Engystomops pustulosus, and Hymenochirus boettgeri. Frog chromosomes have remained remarkably stable since the Mesozoic Era, with limited Robertsonian (i.e., arm-preserving) translocations and end-to-end fusions found among the smaller chromosomes. Conservation of synteny includes conservation of centromere locations, marked by centromeric tandem repeats associated with Cenp-a binding surrounded by pericentromeric LINE/L1 elements. This work explores the structure of chromosomes across frogs, using a dense meiotic linkage map for X. tropicalis and chromatin conformation capture (Hi-C) data for all species. Abundant satellite repeats occupy the unusually long (~20 megabase) terminal regions of each chromosome that coincide with high rates of recombination. Both embryonic and differentiated cells show reproducible associations of centromeric chromatin and of telomeres, reflecting a Rabl-like configuration. Our comparative analyses reveal 13 conserved ancestral anuran chromosomes from which contemporary frog genomes were constructed.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Functional anatomy of zinc finger antiviral protein complexes

Abstract ZAP is an antiviral protein that binds to and depletes viral RNA, which is often distinguished from vertebrate host RNA by its elevated CpG content. Two ZAP cofactors, TRIM25 and KHNYN, have activities that are poorly understood. Here, we show that functional interactions between ZAP, TRIM25 and KHNYN involve multiple domains of each protein, and that the ability of TRIM25 to multimerize via its RING domain augments ZAP activity and specificity. We show that KHNYN is an active nuclease that acts in a partly redundant manner with its homolog N4BP1. The ZAP N-terminal RNA binding domain constitutes a minimal core that is essential for antiviral complex activity, and we present a crystal structure of this domain that reveals contacts with the functionally required KHNYN C-terminal domain. These contacts are remote from the ZAP CpG binding site and would not interfere with RNA binding. Based on our dissection of ZAP, TRIM25 and KHNYN functional anatomy, we could design artificial chimeric antiviral proteins that reconstitute the antiviral function of the intact authentic proteins, but in the absence of protein domains that are otherwise required for activity. Together, these results suggest a model for the RNA recognition and action of ZAP-containing antiviral protein complexes.

Science & Technology - Other Topics↗
DOE OSTI2025

Feeding state-dependent neuropeptidergic modulation of reciprocally interconnected inhibitory neurons biases sensorimotor decisions in Drosophila

Abstract An animal’s feeding state changes its behavioral priorities and thus influences even nonfeeding-related decisions. How the feeding state information is transmitted to nonfeeding-related circuits and what circuit mechanisms are involved in biasing nonfeeding-related decisions remain open questions. By combining calcium imaging, neuronal manipulations, behavioral analysis and computational modeling, we determined that the competition between different aversive responses to mechanical cues is biased by changes in the feeding state. We found that this effect is achieved by the differential modulation of two different types of reciprocally connected inhibitory neurons promoting opposing actions. This modulation results in a more frequent active type of response and, less frequently, a protective type of response if larvae are fed sugar than when they are fed a balanced diet. Information about the internal state is conveyed to inhibitory neurons through homologs of the vertebrate neuropeptide Y, which is known to be involved in regulating feeding behavior.

Science & Technology - Other Topics↗
DOE OSTI2025

PCM1 coordinates centrosome asymmetry with polarized endosome dynamics to regulate daughter cell fate

Vertebrate radial glia progenitors (RGPs) balance self-renewal and differentiation through asymmetric cell division (ACD), which involves unequal centrosome inheritance. How centrosome asymmetry directs cell fate remains poorly understood. Here, we identify Pericentriolar material 1 (Pcm1) as a key player in this process. In zebrafish embryonic RGPs, Pcm1 is asymmetrically associated with Cep83, a mother centrosome marker. Using in vivo time-lapse imaging and nanoscale-resolution expansion microscopy, we detect Pcm1 on Notch ligand-containing endosomes, where it interacts–either directly or indirectly–with Par-3 and dynein. Loss of pcm1 disrupts endosome dynamics, increasing neuronal differentiation at the expense of RGP self-renewal. Mechanistically, Pcm1 facilitates the transition from Rab5b to Rab11a and promotes the assembly of Par-3 and dynein macromolecular complexes on recycling endosomes. Furthermore, we find conserved PARD3-PCM1-CEP83-RAB11 associations in human cortical brain organoids. Our findings uncover that Pcm1 links centrosome asymmetry to polarized endosome trafficking, thereby regulating RGP fate decisions.

Cell fate and cell lineage↗
DOE OSTI2021

Metal Ion Fluxes Controlling Amphibian Fertilization

Mammalian oocytes undergo major changes in zinc content and localization to be fertilized, the most striking being the rapid exocytosis of over 10 billion zinc ions in what are known as zinc sparks. Here, we report that fertilization of amphibian Xenopus laevis eggs also initiates a zinc spark that progresses across the cell surface in coordination with dynamic calcium waves. This zinc exocytosis is accompanied by a newly recognized loss of intracellular manganese. Synchrotron-based X-ray fluorescence and analytical electron microscopy reveal that zinc and manganese are sequestered in a system of cortical granules that are abundant at the animal pole. Through electron-nuclear double-resonance studies, we rule out Mn2+ complexation with phosphate or nitrogenous ligands in intact eggs, but the data are consistent with a carboxylate coordination environment. Our observations suggest that zinc and manganese fluxes are a conserved feature of fertilization in vertebrates and that they function as part of a physiological block to polyspermy.

Seeler, John F.↗
DOE OSTI2023

Detecting macroevolutionary genotype–phenotype associations using error-corrected rates of protein convergence

On macroevolutionary timescales, extensive mutations and phylogenetic uncertainty mask the signals of genotype–phenotype associations underlying convergent evolution. To overcome this problem, we extended the widely used framework of non-synonymous to synonymous substitution rate ratios and developed the novel metric ω C , which measures the error-corrected convergence rate of protein evolution. While ω C distinguishes natural selection from genetic noise and phylogenetic errors in simulation and real examples, its accuracy allows an exploratory genome-wide search of adaptive molecular convergence without phenotypic hypothesis or candidate genes. Using gene expression data, we explored over 20 million branch combinations in vertebrate genes and identified the joint convergence of expression patterns and protein sequences with amino acid substitutions in functionally important sites, providing hypotheses on undiscovered phenotypes. We further extended our method with a heuristic algorithm to detect highly repetitive convergence among computationally non-trivial higher-order phylogenetic combinations. Our approach allows bidirectional searches for genotype–phenotype associations, even in lineages that diverged for hundreds of millions of years.

59 BASIC BIOLOGICAL SCIENCES↗