Endothermic physiology of extinct megatooth sharks.

The evolution of the extinct megatooth shark, Otodus megalodon, and its close phylogenetic relatives remains enigmatic. A central question persists regarding the thermophysiological origins of these large predatory sharks through geologic time, including whether O. megalodon was ectothermic or endothermic (including regional endothermy), and whether its thermophysiology could help to explain the iconic shark's gigantism and eventual demise during the Pliocene. To address these uncertainties, we present unique geochemical evidence for thermoregulation in O. megalodon from both clumped isotope paleothermometry and phosphate oxygen isotopes. Our results show that O. megalodon had an overall warmer body temperature compared with its ambient environment and other coexisting shark species, providing quantitative and experimental support for recent biophysical modeling studies that suggest endothermy was one of the key drivers for gigantism in O. megalodon and other lamniform sharks. The gigantic body size with high metabolic costs of having high body temperatures may have contributed to the vulnerability of Otodus species to extinction when compared to other sympatric sharks that survived the Pliocene epoch.

Carbonate clumped isotope analysis (Δ47 ) of 21 carbonate standards determined via gas-source isotope-ratio mass spectrometry on four instrumental configurations using carbonate-based standardization and multiyear data sets.

RATIONALE: Clumped isotope geochemistry examines the pairing or clumping of heavy isotopes in molecules and provides information about the thermodynamic and kinetic controls on their formation. The first clumped isotope measurements of carbonate minerals were first published 15 years ago, and since then, interlaboratory offsets have been observed, and laboratory and community practices for measurement, data analysis, and instrumentation have evolved. Here we briefly review historical and recent developments for measurements, share Tripati Lab practices for four different instrument configurations, test a recently published proposal for carbonate-based standardization on multiple instruments using multi-year data sets, and report values for 21 different carbonate standards that allow for recalculations of previously published data sets. METHODS: We examine data from 4628 standard measurements on Thermo MAT 253 and Nu Perspective IS mass spectrometers, using a common acid bath (90°C) and small-sample (70°C) individual reaction vessels. Each configuration was investigated by treating some standards as anchors (working standards) and the remainder as unknowns (consistency standards). RESULTS: We show that different acid digestion systems and mass spectrometer models yield indistinguishable results when instrument drift is well characterized. For linearity correction, mixed gas-and-carbonate standardization or carbonate-only standardization yields similar results. No difference is observed in the use of three or eight working standards for the construction of transfer functions. CONCLUSIONS: We show that all configurations yield similar results if instrument drift is robustly characterized and validate a recent proposal for carbonate-based standardization using large multiyear data sets. Δ47 values are reported for 21 carbonate standards on both the absolute reference frame (ARF; also refered to as the Carbon Dioxide Equilibrated Scale or CDES) and the new InterCarb-Carbon Dioxide Equilibrium Scale (I-CDES) reference frame, facilitating intercomparison of data from a diversity of labs and instrument configurations and restandardization of a broad range of sample sets between 2006, when the first carbonate measurements were published, and the present.

Olfaction written in bone: cribriform plate size parallels olfactory receptor gene repertoires in Mammalia.

The evolution of mammalian olfaction is manifested in a remarkable diversity of gene repertoires, neuroanatomy and skull morphology across living species. Olfactory receptor genes (ORGs), which initiate the conversion of odorant molecules into odour perceptions and help an animal resolve the olfactory world, range in number from a mere handful to several thousand genes across species. Within the snout, each of these ORGs is exclusively expressed by a discrete population of olfactory sensory neurons (OSNs), suggesting that newly evolved ORGs may be coupled with new OSN populations in the nasal epithelium. Because OSN axon bundles leave high-fidelity perforations (foramina) in the bone as they traverse the cribriform plate (CP) to reach the brain, we predicted that taxa with larger ORG repertoires would have proportionately expanded footprints in the CP foramina. Previous work found a correlation between ORG number and absolute CP size that disappeared after accounting for body size. Using updated, digital measurement data from high-resolution CT scans and re-examining the relationship between CP and body size, we report a striking linear correlation between relative CP area and number of functional ORGs across species from all mammalian superorders. This correlation suggests strong developmental links in the olfactory pathway between genes, neurons and skull morphology. Furthermore, because ORG number is linked to olfactory discriminatory function, this correlation supports relative CP size as a viable metric for inferring olfactory capacity across modern and extinct species. By quantifying CP area from a fossil sabertooth cat (Smilodon fatalis), we predicted a likely ORG repertoire for this extinct felid.

High regional climate sensitivity over continental China constrained by glacial-recent changes in temperature and the hydrological cycle.

The East Asian monsoon is one of Earth's most significant climatic phenomena, and numerous paleoclimate archives have revealed that it exhibits variations on orbital and suborbital time scales. Quantitative constraints on the climate changes associated with these past variations are limited, yet are needed to constrain sensitivity of the region to changes in greenhouse gas levels. Here, we show central China is a region that experienced a much larger temperature change since the Last Glacial Maximum than typically simulated by climate models. We applied clumped isotope thermometry to carbonates from the central Chinese Loess Plateau to reconstruct temperature and water isotope shifts from the Last Glacial Maximum to present. We find a summertime temperature change of 6-7 °C that is reproduced by climate model simulations presented here. Proxy data reveal evidence for a shift to lighter isotopic composition of meteoric waters in glacial times, which is also captured by our model. Analysis of model outputs suggests that glacial cooling over continental China is significantly amplified by the influence of stationary waves, which, in turn, are enhanced by continental ice sheets. These results not only support high regional climate sensitivity in Central China but highlight the fundamental role of planetary-scale atmospheric dynamics in the sensitivity of regional climates to continental glaciation, changing greenhouse gas levels, and insolation.

Body temperatures of modern and extinct vertebrates from (13)C-(18)O bond abundances in bioapatite.

The stable isotope compositions of biologically precipitated apatite in bone, teeth, and scales are widely used to obtain information on the diet, behavior, and physiology of extinct organisms and to reconstruct past climate. Here we report the application of a new type of geochemical measurement to bioapatite, a "clumped-isotope" paleothermometer, based on the thermodynamically driven preference for (13)C and (18)O to bond with each other within carbonate ions in the bioapatite crystal lattice. This effect is dependent on temperature but, unlike conventional stable isotope paleothermometers, is independent from the isotopic composition of water from which the mineral formed. We show that the abundance of (13)C-(18)O bonds in the carbonate component of tooth bioapatite from modern specimens decreases with increasing body temperature of the animal, following a relationship between isotope "clumping" and temperature that is statistically indistinguishable from inorganic calcite. This result is in agreement with a theoretical model of isotopic ordering in carbonate ion groups in apatite and calcite. This thermometer constrains body temperatures of bioapatite-producing organisms with an accuracy of 1-2 degrees C. Analyses of fossilized tooth enamel of both Pleistocene and Miocene age yielded temperatures within error of those derived from similar modern taxa. Clumped-isotope analysis of bioapatite represents a new approach in the study of the thermophysiology of extinct species, allowing the first direct measurement of their body temperatures. It will also open new avenues in the study of paleoclimate, as the measurement of clumped isotopes in phosphorites and fossils has the potential to reconstruct environmental temperatures.

Environmental DNA metabarcoding reveals distinct fish assemblages supported by seagrass (Zostera marina and Zostera pacifica) beds in different geographic settings in Southern California.

Seagrass beds are disappearing at a record pace despite their known value to our oceans and coastal communities. Simultaneously, our coastlines are under the constant pressure of climate change which is impacting their chemical, physical and biological characteristics. It is thus pertinent to evaluate and record habitat use so we can understand how these different environments contribute to local biodiversity. This study evaluates the assemblages of fish found at five Zostera beds in Southern California using environmental DNA (eDNA) metabarcoding. eDNA is a powerful biodiversity monitoring tool that offers key advantages to conventional monitoring. Results from our eDNA study found 78 species of fish that inhabit these five beds around Southern California representing embayment, open coastal mainland and open coastal island settings. While each bed had the same average number of species found throughout the year, the composition of these fish assemblages was strongly site dependent. There were 35 fish that were found at both open coast and embayment seagrass beds, while embayment seagrass sites had 20 unique fish and open coast sites had 23 unique fish. These results demonstrate that seagrass fish assemblages are heterogenous based on their geographic positioning and that marine managers must take this into account for holistic conservation and restoration efforts.

Equilibrated Gas and Carbonate Standard-Derived Dual (Δ47 and Δ48) Clumped Isotope Values.

Carbonate clumped isotope geochemistry has primarily focused on mass spectrometric determination of m/z 47 CO2 for geothermometry, but theoretical calculations and recent experiments indicate paired analysis of the m/z 47 (13C18O16O) and m/z 48 (12C18O18O) isotopologues (referred to as Δ47 and Δ48) can be used to study non-equilibrium isotope fractionations and refine temperature estimates. We utilize 5,448 Δ47 and 3,400 Δ48 replicate measurements of carbonate samples and standards, and 183 Δ47 and 195 Δ48 replicate measurements of gas standards from 2015 to 2021 from a multi-year and multi-instrument data set to constrain Δ47 and Δ48 values for 27 samples and standards, including Devils Hole cave calcite, and study equilibrium Δ47-Δ48, Δ47-temperature, and Δ48-temperature relationships. We compare results to previously published findings and calculate equilibrium regressions based on data from multiple laboratories. We report acid digestion fractionation factors, Δ*63-47 and Δ*64-48, and account for their dependence on the initial clumped isotope values of the mineral.

NKG2D ligand tumor expression and association with clinical outcome in early breast cancer patients: an observational study.

BACKGROUND: Cell surface NKG2D ligands (NKG2DL) bind to the activating NKG2D receptor present on NK cells and subsets of T cells, thus playing a role in initiating an immune response. We examined tumor expression and prognostic effect of NKG2DL in breast cancer patients. METHODS: Our study population (n = 677) consisted of all breast cancer patients primarily treated with surgery in our center between 1985 and 1994. Formalin-fixed paraffin-embedded tumor tissue was immunohistochemically stained with antibodies directed against MIC-A/MIC-B (MIC-AB), ULBP-1, ULBP-2, ULBP-3, ULBP-4, and ULBP-5. RESULTS: NKG2DL were frequently expressed by tumors (MIC-AB, 50% of the cases; ULBP-1, 90%; ULBP-2, 99%; ULBP-3, 100%; ULBP-4, 26%; ULBP-5, 90%) and often showed co-expression: MIC-AB and ULBP-4 (p = 0.043), ULBP-1 and ULBP-5 (p = 0.006), ULBP-4 and ULBP-5 (p < 0.001). MIC-AB (p = 0.001) and ULBP-2 (p = 0.006) expression resulted in a statistically significant longer relapse free period (RFP). Combined expression of these ligands showed to be an independent prognostic parameter for RFP (p < 0.001, HR 0.41). Combined expression of all ligands showed no associations with clinical outcome. CONCLUSIONS: We demonstrated for the first time that NKG2DL are frequently expressed and often co-expressed in breast cancer. Expression of MIC-AB and ULBP-2 resulted in a statistically significant beneficial outcome concerning RFP with high discriminative power. Combination of all NKG2DL showed no additive or interactive effect of ligands on each other, suggesting that similar and co-operative functioning of all NKG2DL can not be assumed. Our observations suggest that among driving forces in breast cancer outcome are immune activation on one site and tumor immune escape on the other site.

Trophic position of Otodus megalodon and great white sharks through time revealed by zinc isotopes.

Diet is a crucial trait of an animal's lifestyle and ecology. The trophic level of an organism indicates its functional position within an ecosystem and holds significance for its ecology and evolution. Here, we demonstrate the use of zinc isotopes (δ66Zn) to geochemically assess the trophic level in diverse extant and extinct sharks, including the Neogene megatooth shark (Otodus megalodon) and the great white shark (Carcharodon carcharias). We reveal that dietary δ66Zn signatures are preserved in fossil shark tooth enameloid over deep geologic time and are robust recorders of each species' trophic level. We observe significant δ66Zn differences among the Otodus and Carcharodon populations implying dietary shifts throughout the Neogene in both genera. Notably, Early Pliocene sympatric C. carcharias and O. megalodon appear to have occupied a similar mean trophic level, a finding that may hold clues to the extinction of the gigantic Neogene megatooth shark.

Cenozoic megatooth sharks occupied extremely high trophic positions.

Trophic position is a fundamental characteristic of animals, yet it is unknown in many extinct species. In this study, we ground-truth the 15N/14N ratio of enameloid-bound organic matter (δ15NEB) as a trophic level proxy by comparison to dentin collagen δ15N and apply this method to the fossil record to reconstruct the trophic level of the megatooth sharks (genus Otodus). These sharks evolved in the Cenozoic, culminating in Otodus megalodon, a shark with a maximum body size of more than 15 m, which went extinct 3.5 million years ago. Very high δ15NEB values (22.9 ± 4.4‰) of O. megalodon from the Miocene and Pliocene show that it occupied a higher trophic level than is known for any marine species, extinct or extant. δ15NEB also indicates a dietary shift in sharks of the megatooth lineage as they evolved toward the gigantic O. megalodon, with the highest trophic level apparently reached earlier than peak size.

Post-translational modification of the NKG2D ligand RAET1G leads to cell surface expression of a glycosylphosphatidylinositol-linked isoform.

NKG2D is an important activating receptor on lymphocytes. In human, it interacts with two groups of ligands: the major histocompatibility complex class I chain-related A/B (MICA/B) family and the UL-16 binding protein (ULBP) family, also known as retinoic acid early transcript (RAET1). MIC proteins are membrane-anchored, but all of the ULBP/RAET1 proteins, except for RAET1E and RAET1G, are glycosylphosphatidylinositol (GPI)-anchored. To address the reason for these differences we studied the association of RAET1G with the membrane. Using epitope-tagged RAET1G protein in conjunction with antibodies to different parts of the molecule and in pulse-chase experiments, we showed that the C terminus of the protein was cleaved soon after protein synthesis. Endoglycosidase H and peptide N-glycosidase treatment and cell surface immunoprecipitation indicated that most of the protein stayed in the endoplasmic reticulum, but some of the cleaved form was modified in the Golgi and transported to the cell surface. We examined the possibility of GPI anchoring of the protein in three ways: (i) Phosphatidylinositol (PI)-specific phospholipase C released the PI-linked form of the protein. (ii) The surface expression pattern of RAET1G decreased in cells defective in GPI anchoring through mutant GPI-amidase. (iii) Site-directed mutagenesis, to disrupt residues predicted to facilitate GPI-anchoring, resulted in diminished surface expression of RAET1G. Thus, a form of RAET1G is GPI-anchored, in line with most other ULBP/RAET1 family proteins. The cytoplasmic tail and transmembrane domains appear to result from gene duplication and frameshift mutation. Together with our previous results, our data suggest that RAET1G is regulated post-translationally to produce a GPI-anchored isoform.

Thermal stress reduces pocilloporid coral resilience to ocean acidification by impairing control over calcifying fluid chemistry.

The combination of thermal stress and ocean acidification (OA) can more negatively affect coral calcification than an individual stressors, but the mechanism behind this interaction is unknown. We used two independent methods (microelectrode and boron geochemistry) to measure calcifying fluid pH (pHcf) and carbonate chemistry of the corals Pocillopora damicornis and Stylophora pistillata grown under various temperature and pCO2 conditions. Although these approaches demonstrate that they record pHcf over different time scales, they reveal that both species can cope with OA under optimal temperatures (28°C) by elevating pHcf and aragonite saturation state (Ωcf) in support of calcification. At 31°C, neither species elevated these parameters as they did at 28°C and, likewise, could not maintain substantially positive calcification rates under any pH treatment. These results reveal a previously uncharacterized influence of temperature on coral pHcf regulation-the apparent mechanism behind the negative interaction between thermal stress and OA on coral calcification.

Patterns of Element Incorporation in Calcium Carbonate Biominerals Recapitulate Phylogeny for a Diverse Range of Marine Calcifiers.

Elemental ratios in biogenic marine calcium carbonates are widely used in geobiology, environmental science, and paleoenvironmental reconstructions. It is generally accepted that the elemental abundance of biogenic marine carbonates reflects a combination of the abundance of that ion in seawater, the physical properties of seawater, the mineralogy of the biomineral, and the pathways and mechanisms of biomineralization. Here we report measurements of a suite of nine elemental ratios (Li/Ca, B/Ca, Na/Ca, Mg/Ca, Zn/Ca, Sr/Ca, Cd/Ca, Ba/Ca, and U/Ca) in 18 species of benthic marine invertebrates spanning a range of biogenic carbonate polymorph mineralogies (low-Mg calcite, high-Mg calcite, aragonite, mixed mineralogy) and of phyla (including Mollusca, Echinodermata, Arthropoda, Annelida, Cnidaria, Chlorophyta, and Rhodophyta) cultured at a single temperature (25°C) and a range of pCO2 treatments (ca. 409, 606, 903, and 2856 ppm). This dataset was used to explore various controls over elemental partitioning in biogenic marine carbonates, including species-level and biomineralization-pathway-level controls, the influence of internal pH regulation compared to external pH changes, and biocalcification responses to changes in seawater carbonate chemistry. The dataset also enables exploration of broad scale phylogenetic patterns of elemental partitioning across calcifying species, exhibiting high phylogenetic signals estimated from both uni- and multivariate analyses of the elemental ratio data (univariate: λ = 0-0.889; multivariate: λ = 0.895-0.99). Comparing partial R 2 values returned from non-phylogenetic and phylogenetic regression analyses echo the importance of and show that phylogeny explains the elemental ratio data 1.4-59 times better than mineralogy in five out of nine of the elements analyzed. Therefore, the strong associations between biomineral elemental chemistry and species relatedness suggests mechanistic controls over element incorporation rooted in the evolution of biomineralization mechanisms.

ULBP2 and RAET1E NKG2D ligands are independent predictors of poor prognosis in ovarian cancer patients.

The human activating immune receptor, NKG2D, binds to a diverse array of cellular ligands of the MIC and unique long 16 (UL16)-binding protein (ULBP)/retinoic acid early transcript (RAET) family. NKG2D is thought to participate in anticancer immune responses. By using tissue microarrays representing over 300 patients with defined clinicopathological factors, we present the first comprehensive screen of the expression of all NKG2D ligands in primary ovarian cancers. NKG2D ligands were expressed by the majority of tumors; however, the level of expression varied considerably. By categorizing each tumor as having negative, low or high expression, it was shown that high expression of several NKG2D ligands is inversely correlated with disease survival. Patients whose tumors had high expression of RAET1E (p = 0.037), ULBP1 (p = 0.036) and ULBP3 (p = 0.004) surviving a median of 11, 14 and 11 months, respectively, compared with disease-specific survival of 29, 30 and 25 months in patients whose tumors showed no expression of these ligands. These results contrast with previous findings showing that high level NKG2D ligand expression is associated with good prognosis in colorectal cancer and suggest a fundamental difference in the involvement of NKG2D-mediated immunity in these two types of cancer. By using multivariate analysis, the factors retaining independent prognostic significance were International Federation of Gynecologists and Obstetricians stage (p < 0.001), presence of residual disease (p = 0.003), ULBP2 (p = 0.042) and RAET1E (p = 0.030).

ULBP6/RAET1L is an additional human NKG2D ligand.

To date five ULBP/RAET (UL16-binding protein, also known as retinoic acid early transcript) genes, encoded on human chromosome 6q24.2-q25.3, have been shown to encode ligands of the activating immunoreceptor NKG2D. Here, we show that a sixth gene, ULBP6/RAET1L, is a polymorphic locus that expresses a functional transcript. ULBP6 had a more restricted expression profile in cell lines and primary human tissues than other NKG2D ligands, but expression was detected in several human papillomavirus-positive cervical carcinoma cell lines and was inducible on infection with human CMV. ULBP6 bound to recombinant NKG2D as well as the human CMV immune evasion molecule UL16. By confocal microscopy we show that UL16 retains ULBP6 inside the cell, preventing it from reaching the cell surface. Expression of ULBP6 on target cells induced a significant increase in NK-cell killing. Comparison of ULBP6 with ULBP4 and ULBP5 indicated that differences in recombinant NKG2D binding correlated with differences in NK-cell activation.

Regulation of calcification site pH is a polyphyletic but not always governing response to ocean acidification.

The response of marine-calcifying organisms to ocean acidification (OA) is highly variable, although the mechanisms behind this variability are not well understood. Here, we use the boron isotopic composition (δ11B) of biogenic calcium carbonate to investigate the extent to which organisms' ability to regulate pH at their site of calcification (pHCF) determines their calcification responses to OA. We report comparative δ11B analyses of 10 species with divergent calcification responses (positive, parabolic, threshold, and negative) to OA. Although the pHCF is closely coupled to calcification responses only in 3 of the 10 species, all 10 species elevate pHCF above pHsw under elevated pCO2. This result suggests that these species may expend additional energy regulating pHCF under future OA. This strategy of elevating pHCF above pHsw appears to be a polyphyletic, if not universal, response to OA among marine calcifiers-although not always the principal factor governing a species' response to OA.

A coastal coccolithophore maintains pH homeostasis and switches carbon sources in response to ocean acidification.

Ocean acidification will potentially inhibit calcification by marine organisms; however, the response of the most prolific ocean calcifiers, coccolithophores, to this perturbation remains under characterized. Here we report novel chemical constraints on the response of the widespread coccolithophore species Ochrosphaera neapolitana (O. neapolitana) to changing-CO2 conditions. We cultured this algae under three pCO2-controlled seawater pH conditions (8.05, 8.22, and 8.33). Boron isotopes within the algae's extracellular calcite plates show that this species maintains a constant pH at the calcification site, regardless of CO2-induced changes in pH of the surrounding seawater. Carbon and oxygen isotopes in the algae's calcite plates and carbon isotopes in the algae's organic matter suggest that O. neapolitana utilize carbon from a single internal dissolved inorganic carbon (DIC) pool for both calcification and photosynthesis, and that a greater proportion of dissolved CO2 relative to HCO3- enters the internal DIC pool under acidified conditions. These two observations may explain how O. neapolitana continues calcifying and photosynthesizing at a constant rate under different atmospheric-pCO2 conditions.

Comparison of Mediterranean Pteropod Shell Biometrics and Ultrastructure from Historical (1910 and 1921) and Present Day (2012) Samples Provides Baseline for Monitoring Effects of Global Change.

Anthropogenic carbon perturbation has caused decreases in seawater pH and increases in global temperatures since the start of the 20th century. The subsequent lowering of the saturation state of CaCO3 may make the secretion of skeletons more problematic for marine calcifiers. As organisms that precipitate thin aragonite shells, thecosome pteropods have been identified as being particularly vulnerable to climate change effects. Coupled with their global distribution, this makes them ideal for use as sentinel organisms. Recent studies have highlighted shell dissolution as a potential indicator of ocean acidification; however, this metric is not applicable for monitoring pH changes in supersaturated basins. In this study, the novel approach of high resolution computed tomography (CT) scanning was used to produce quantitative 3-dimensional renderings pteropod shells to assess the potential of using this method to monitor small changes in shell biometrics that may be driven by climate change drivers. An ontogenetic analysis of the shells of Cavolinia inflexa and Styliola subula collected from the Mediterranean was used to identify suitable monitoring metrics. Modern samples were then compared to historical samples of the same species, collected during the Mediterranean leg of the Thor (1910) and Dana (1921) cruises to assess whether any empirical differences could be detected. Shell densities were calculated and scanning electron microscopy was used to compare the aragonite crystal morphology. pH for the collection years was hind-cast using temperature and salinity time series with atmospheric CO2 concentrations from ice core data. Historical samples of S. subula were thicker than S. subula shells of the same size from 2012 and C. inflexa shells collected in 1910 were significantly denser than those from 2012. These results provide a baseline for future work to develop monitoring techniques for climate change in the oceans using the novel approach of high-resolution CT scanning.

Regulation of NKG2D ligand gene expression.

The activating immunoreceptor NKG2D has seven known host ligands encoded by the MHC class I chain-related MIC and ULBP/RAET genes. Why there is such diversity of NKG2D ligands is not known but one hypothesis is that they are differentially expressed in different tissues in response to different stresses. To explore this, we compared expression patterns and promoters of NKG2D ligand genes. ULBP/RAET genes were transcribed independent of each other in a panel of cell lines. ULBP/RAET gene expression was upregulated on infection with human cytomegalovirus; however, a clinical strain, Toledo, induced expression more slowly than did a laboratory strain, AD169. ULBP4/RAET1E was not induced by infection with either strain. To investigate the mechanisms behind the similarities and differences in NKG2D ligand gene expression a comparative sequence analysis of NKG2D ligand gene putative promoter regions was conducted. Sequence alignments demonstrated that there was significant sequence diversity; however, one region of high similarity between most of the genes is evident. This region contains a number of potential transcription factor binding sites, including those involved in shock responses and sites for retinoic acid-induced factors. Promoters of some NKG2D ligand genes are polymorphic and several sequence alterations in these alleles abolished putative transcription factor binding.