A six-amino-acid motif is a major determinant in functional evolution of HOX1 proteins.

Gene duplication and divergence is a major driver in the emergence of evolutionary novelties. How variations in amino acid sequences lead to loss of ancestral activity and functional diversification of proteins is poorly understood. We used cross-species functional analysis of Drosophila Labial and its mouse HOX1 orthologs (HOXA1, HOXB1, and HOXD1) as a paradigm to address this issue. Mouse HOX1 proteins display low (30%) sequence similarity with Drosophila Labial. However, substituting endogenous Labial with the mouse proteins revealed that HOXA1 has retained essential ancestral functions of Labial, while HOXB1 and HOXD1 have diverged. Genome-wide analysis demonstrated similar DNA-binding patterns of HOXA1 and Labial in mouse cells, while HOXB1 binds to distinct targets. Compared with HOXB1, HOXA1 shows an enrichment in co-occupancy with PBX proteins on target sites and exists in the same complex with PBX on chromatin. Functional analysis of HOXA1-HOXB1 chimeric proteins uncovered a novel six-amino-acid C-terminal motif (CTM) flanking the homeodomain that serves as a major determinant of ancestral activity. In vitro DNA-binding experiments and structural prediction show that CTM provides an important domain for interaction of HOXA1 proteins with PBX. Our findings show that small changes outside of highly conserved DNA-binding regions can lead to profound changes in protein function.

Shared retinoic acid responsive enhancers coordinately regulate nascent transcription of Hoxb coding and non-coding RNAs in the developing mouse neural tube.

Signaling pathways regulate the patterns of Hox gene expression that underlie their functions in the specification of axial identity. Little is known about the properties of cis-regulatory elements and underlying transcriptional mechanisms that integrate graded signaling inputs to coordinately control Hox expression. Here, we optimized a single molecule fluorescent in situ hybridization (smFISH) technique with probes spanning introns to evaluate how three shared retinoic acid response element (RARE)-dependent enhancers in the Hoxb cluster regulate patterns of nascent transcription in vivo at the level of single cells in wild-type and mutant embryos. We predominately detect nascent transcription of only a single Hoxb gene in each cell, with no evidence for simultaneous co-transcriptional coupling of all or specific subsets of genes. Single and/or compound RARE mutations indicate that each enhancer differentially impacts global and local patterns of nascent transcription, suggesting that selectivity and competitive interactions between these enhancers is important to robustly maintain the proper levels and patterns of nascent Hoxb transcription. This implies that rapid and dynamic regulatory interactions potentiate transcription of genes through combined inputs from these enhancers in coordinating the retinoic acid response.

The Hox protein conundrum: The "specifics" of DNA binding for Hox proteins and their partners.

Homeotic genes (Hox genes) are homeodomain-transcription factors involved in conferring segmental identity along the anterior-posterior body axis. Molecular characterization of HOX protein function raises some interesting questions regarding the source of the binding specificity of the HOX proteins. How do HOX proteins regulate common and unique target specificity across space and time? This review attempts to summarize and interpret findings in this area, largely focused on results from in vitro and in vivo studies in Drosophila and mouse systems. Recent studies related to HOX protein binding specificity compel us to reconsider some of our current models for transcription factor-DNA interactions. It is crucial to study transcription factor binding by incorporating components of more complex, multi-protein interactions in concert with small changes in binding motifs that can significantly impact DNA binding specificity and subsequent alterations in gene expression. To incorporate the multiple elements that can determine HOX protein binding specificity, we propose a more integrative Cooperative Binding model.

Analysis of lamprey meis genes reveals that conserved inputs from Hox, Meis and Pbx proteins control their expression in the hindbrain and neural tube.

Meis genes are known to play important roles in the hindbrain and neural crest cells of jawed vertebrates. To explore the roles of Meis genes in head development during evolution of vertebrates, we have identified four meis genes in the sea lamprey genome and characterized their patterns of expression and regulation, with a focus on the hindbrain and pharynx. Each of the lamprey meis genes displays temporally and spatially dynamic patterns of expression, some of which are coupled to rhombomeric domains in the developing hindbrain and select pharyngeal arches. Studies of Meis loci in mouse and zebrafish have identified enhancers that are bound by Hox and TALE (Meis and Pbx) proteins, implicating these factors in the direct regulation of Meis expression. We examined the lamprey meis loci and identified a series of cis-elements conserved between lamprey and jawed vertebrate meis genes. In transgenic reporter assays we demonstrated that these elements act as neural enhancers in lamprey embryos, directing reporter expression in appropriate domains when compared to expression of their associated endogenous meis gene. Sequence alignments reveal that these conserved elements are in similar relative positions of the meis loci and contain a series of consensus binding motifs for Hox and TALE proteins. This suggests that ancient Hox and TALE-responsive enhancers regulated expression of ancestral vertebrate meis genes in segmental domains in the hindbrain and have been retained in the meis loci during vertebrate evolution. The presence of conserved Meis, Pbx and Hox binding sites in these lamprey enhancers links Hox and TALE factors to regulation of lamprey meis genes in the developing hindbrain, indicating a deep ancestry for these regulatory interactions prior to the divergence of jawed and jawless vertebrates.

Cancer-specific mutation of GATA3 disrupts the transcriptional regulatory network governed by Estrogen Receptor alpha, FOXA1 and GATA3.

Estrogen receptors (ER) are activated by the steroid hormone 17ß-estradiol. Estrogen receptor alpha (ER-α) forms a regulatory network in mammary epithelial cells and in breast cancer with the transcription factors FOXA1 and GATA3. GATA3 is one of the most frequently mutated genes in breast cancer and is capable of specifying chromatin localization of FOXA1 and ER-α. How GATA3 mutations found in breast cancer impact genomic localization of ER-α and the transcriptional network downstream of ER-α and FOXA1 remains unclear. Here, we investigate the function of a recurrent patient-derived GATA3 mutation (R330fs) on this regulatory network. Genomic analysis indicates that the R330fs mutant can disrupt localization of ER-α and FOXA1. Loci co-bound by all three factors are enriched for genes integral to mammary gland development as well as epithelial cell biology. This gene set is differentially regulated in GATA3 mutant cells in culture and in tumors bearing similar mutations in vivo. The altered distribution of ER-α and FOXA1 in GATA3-mutant cells is associated with altered chromatin architecture, which leads to differential gene expression. These results suggest an active role for GATA3 zinc finger 2 mutants in ER-α positive breast tumors.

HOXA1 and TALE proteins display cross-regulatory interactions and form a combinatorial binding code on HOXA1 targets.

Hoxa1 has diverse functional roles in differentiation and development. We identify and characterize properties of regions bound by HOXA1 on a genome-wide basis in differentiating mouse ES cells. HOXA1-bound regions are enriched for clusters of consensus binding motifs for HOX, PBX, and MEIS, and many display co-occupancy of PBX and MEIS. PBX and MEIS are members of the TALE family and genome-wide analysis of multiple TALE members (PBX, MEIS, TGIF, PREP1, and PREP2) shows that nearly all HOXA1 targets display occupancy of one or more TALE members. The combinatorial binding patterns of TALE proteins define distinct classes of HOXA1 targets, which may create functional diversity. Transgenic reporter assays in zebrafish confirm enhancer activities for many HOXA1-bound regions and the importance of HOX-PBX and TGIF motifs for their regulation. Proteomic analyses show that HOXA1 physically interacts on chromatin with PBX, MEIS, and PREP family members, but not with TGIF, suggesting that TGIF may have an independent input into HOXA1-bound regions. Therefore, TALE proteins appear to represent a wide repertoire of HOX cofactors, which may coregulate enhancers through distinct mechanisms. We also discover extensive auto- and cross-regulatory interactions among the Hoxa1 and TALE genes, indicating that the specificity of HOXA1 during development may be regulated though a complex cross-regulatory network of HOXA1 and TALE proteins. This study provides new insight into a regulatory network involving combinatorial interactions between HOXA1 and TALE proteins.

Dynamic regulation of Nanog and stem cell-signaling pathways by Hoxa1 during early neuro-ectodermal differentiation of ES cells.

Homeobox a1 (Hoxa1) is one of the most rapidly induced genes in ES cell differentiation and it is the earliest expressed Hox gene in the mouse embryo. In this study, we used genomic approaches to identify Hoxa1-bound regions during early stages of ES cell differentiation into the neuro-ectoderm. Within 2 h of retinoic acid treatment, Hoxa1 is rapidly recruited to target sites that are associated with genes involved in regulation of pluripotency, and these genes display early changes in expression. The pattern of occupancy of Hoxa1 is dynamic and changes over time. At 12 h of differentiation, many sites bound at 2 h are lost and a new cohort of bound regions appears. At both time points the genome-wide mapping reveals that there is significant co-occupancy of Nanog (Nanog homeobox) and Hoxa1 on many common target sites, and these are linked to genes in the pluripotential regulatory network. In addition to shared target genes, Hoxa1 binds to regulatory regions of Nanog, and conversely Nanog binds to a 3' enhancer of Hoxa1 This finding provides evidence for direct cross-regulatory feedback between Hoxa1 and Nanog through a mechanism of mutual repression. Hoxa1 also binds to regulatory regions of Sox2 (sex-determining region Y box 2), Esrrb (estrogen-related receptor beta), and Myc, which underscores its key input into core components of the pluripotential regulatory network. We propose a model whereby direct inputs of Nanog and Hoxa1 on shared targets and mutual repression between Hoxa1 and the core pluripotency network provides a molecular mechanism that modulates the fine balance between the alternate states of pluripotency and differentiation.

Hox genes: Downstream "effectors" of retinoic acid signaling in vertebrate embryogenesis.

One of the major regulatory challenges of animal development is to precisely coordinate in space and time the formation, specification, and patterning of cells that underlie elaboration of the basic body plan. How does the vertebrate plan for the nervous and hematopoietic systems, heart, limbs, digestive, and reproductive organs derive from seemingly similar population of cells? These systems are initially established and patterned along the anteroposterior axis (AP) by opposing signaling gradients that lead to the activation of gene regulatory networks involved in axial specification, including the Hox genes. The retinoid signaling pathway is one of the key signaling gradients coupled to the establishment of axial patterning. The nested domains of Hox gene expression, which provide a combinatorial code for axial patterning, arise in part through a differential response to retinoic acid (RA) diffusing from anabolic centers established within the embryo during development. Hence, Hox genes are important direct effectors of retinoid signaling in embryogenesis. This review focuses on describing current knowledge on the complex mechanisms and regulatory processes, which govern the response of Hox genes to RA in several tissue contexts including the nervous system during vertebrate development.

Dynamic transcriptional events in embryonic stem cells mediated by the super elongation complex (SEC).

Transcriptional regulation of developmentally controlled genes is at the heart of differentiation and organogenesis. In this study, we performed global genomic analyses in murine embryonic stem (ES) cells and in human cells in response to activation signals. We identified an essential role for the ELL (eleven-nineteen lysine-rich leukemia gene)/P-TEFb (positive transcription elongation factor)-containing super elongation complex (SEC) in the regulation of gene expression, including several genes bearing paused RNA polymerase II (Pol II). Paused Pol II has been proposed to be associated with loci that respond rapidly to environmental stimuli. However, our studies in ES cells also identified a requirement for SEC at genes without paused Pol II, which also respond dynamically to differentiation signals. Our findings suggest that SEC is a major class of active P-TEFb-containing complexes required for transcriptional activation in response to environmental cues such as differentiation signals.

Hoxa1 targets signaling pathways during neural differentiation of ES cells and mouse embryogenesis.

Hoxa1 has important functional roles in neural crest specification, hindbrain patterning and heart and ear development, yet the enhancers and genes that are targeted by Hoxa1 are largely unknown. In this study, we performed a comprehensive analysis of Hoxa1 target genes using genome-wide Hoxa1 binding data in mouse ES cells differentiated with retinoic acid (RA) into neural fates in combination with differential gene expression analysis in Hoxa1 gain- and loss-of-function mouse and zebrafish embryos. Our analyses reveal that Hoxa1-bound regions show epigenetic marks of enhancers, occupancy of Hox cofactors and differential expression of nearby genes, suggesting that these regions are enriched for enhancers. In support of this, 80 of them mapped to regions with known reporter activity in transgenic mouse embryos based on the Vista enhancer database. Two additional enhancers in Dok5 and Wls1 were shown to mediate neural expression in developing mouse and zebrafish. Overall, our analysis of the putative target genes indicate that Hoxa1 has input to components of major signaling pathways, including Wnt, TGF-ß, Hedgehog and Hippo, and frequently does so by targeting multiple components of a pathway such as secreted inhibitors, ligands, receptors and down-stream components. We also identified genes implicated in heart and ear development, neural crest migration and neuronal patterning and differentiation, which may underlie major Hoxa1 mutant phenotypes. Finally, we found evidence for a high degree of evolutionary conservation of many binding regions and downstream targets of Hoxa1 between mouse and zebrafish. Our genome-wide analyses in ES cells suggests that we have enriched for in vivo relevant target genes and pathways associated with functional roles of Hoxa1 in mouse development.

Analysis of dynamic changes in retinoid-induced transcription and epigenetic profiles of murine Hox clusters in ES cells.

The clustered Hox genes, which are highly conserved across metazoans, encode homeodomain-containing transcription factors that provide a blueprint for segmental identity along the body axis. Recent studies have underscored that in addition to encoding Hox genes, the homeotic clusters contain key noncoding RNA genes that play a central role in development. In this study, we have taken advantage of genome-wide approaches to provide a detailed analysis of retinoic acid (RA)-induced transcriptional and epigenetic changes within the homeotic clusters of mouse embryonic stem cells. Although there is a general colinear response, our analyses suggest a lack of strict colinearity for several genes in the HoxA and HoxB clusters. We have identified transcribed novel noncoding RNAs (ncRNAs) and their cis-regulatory elements that function in response to RA and demonstrated that the expression of these ncRNAs from both strands represent some of the most rapidly induced transcripts in ES cells. Finally, we have provided dynamic analyses of chromatin modifications for the coding and noncoding genes expressed upon activation and suggest that active transcription can occur in the presence of chromatin modifications and machineries associated with repressed transcription state over the clusters. Overall, our data provide a resource for a better understanding of the dynamic nature of the coding and noncoding transcripts and their associated chromatin marks in the regulation of homeotic gene transcription during development.

The chemokine receptor CXCR3 promotes CD8+ T cell-dependent lung pathology during influenza pathogenesis.

The dual role of CD8+ T cells in influenza control and lung pathology is increasingly appreciated. To explore whether protective and pathological functions can be linked to specific subsets, we dissected CD8+ T responses in influenza-infected murine lungs. Our single-cell RNA-sequencing (scRNA-seq) analysis revealed notable diversity in CD8+ T subpopulations during peak viral load and infection-resolved state. While enrichment of a Cxcr3hi CD8+ T effector subset was associated with a more robust cytotoxic response, both CD8+ T effector and central memory exhibited equally potent effector potential. The scRNA-seq analysis identified unique regulons regulating the cytotoxic response in CD8+ T cells. The late-stage CD8+ T blockade in influenza-cleared lungs or continuous CXCR3 blockade mitigated lung injury without affecting viral clearance. Furthermore, adoptive transfer of wild-type CD8+ T cells exacerbated influenza lung pathology in Cxcr3-/- mice. Collectively, our data imply that CXCR3 interception could have a therapeutic effect in preventing influenza-linked lung injury.

Transcriptional responses of cancer cells to heat shock-inducing stimuli involve amplification of robust HSF1 binding.

Responses of cells to stimuli are increasingly discovered to involve the binding of sequence-specific transcription factors outside of known target genes. We wanted to determine to what extent the genome-wide binding and function of a transcription factor are shaped by the cell type versus the stimulus. To do so, we induced the Heat Shock Response pathway in two different cancer cell lines with two different stimuli and related the binding of its master regulator HSF1 to nascent RNA and chromatin accessibility. Here, we show that HSF1 binding patterns retain their identity between basal conditions and under different magnitudes of activation, so that common HSF1 binding is globally associated with distinct transcription outcomes. HSF1-induced increase in DNA accessibility was modest in scale, but occurred predominantly at remote genomic sites. Apart from regulating transcription at existing elements including promoters and enhancers, HSF1 binding amplified during responses to stimuli may engage inactive chromatin.

Multiscale genetic architecture of donor-recipient differences reveals intronic LIMS1 mismatches associated with kidney transplant survival.

Donor-recipient (D-R) mismatches outside of human leukocyte antigens (HLAs) contribute to kidney allograft loss, but the mechanisms remain unclear, specifically for intronic mismatches. We quantified non-HLA mismatches at variant-, gene-, and genome-wide scales from single nucleotide polymorphism (SNP) data of D-Rs from 2 well-phenotyped transplant cohorts: Genomics of Chronic Allograft Rejection (GoCAR; n = 385) and Clinical Trials in Organ Transplantation-01/17 (CTOT-01/17; n = 146). Unbiased gene-level screening in GoCAR uncovered the LIMS1 locus as the top-ranked gene where D-R mismatches associated with death-censored graft loss (DCGL). A previously unreported, intronic, LIMS1 haplotype of 30 SNPs independently associated with DCGL in both cohorts. Haplotype mismatches showed a dosage effect, and minor-allele introduction to major-allele-carrying recipients showed greater hazard of DCGL. The LIMS1 haplotype and the previously reported LIMS1 SNP rs893403 are expression quantitative trait loci (eQTL) in immune cells for GCC2 (not LIMS1), which encodes a protein involved in mannose-6-phosphase receptor (M6PR) recycling. Peripheral blood and T cell transcriptome analyses associated the GCC2 gene and LIMS1 SNPs with the TGF-ß1/SMAD pathway, suggesting a regulatory effect. In vitro GCC2 modulation impacted M6PR-dependent regulation of active TGF-ß1 and downstream signaling in T cells. Together, our data link LIMS1 locus D-R mismatches to DCGL via GCC2 eQTLs that modulate TGF-ß1-dependent effects on T cells.

Utility of promoter hypermethylation in malignant risk stratification of intraductal papillary mucinous neoplasms.

BACKGROUND: Intraductal papillary mucinous neoplasms (IPMNs), a type of cystic pancreatic cancer (PC) precursors, are increasingly identified on cross-sectional imaging and present a significant diagnostic challenge. While surgical resection of IPMN-related advanced neoplasia, i.e., IPMN-related high-grade dysplasia or PC, is an essential early PC detection strategy, resection is not recommended for IPMN-low-grade dysplasia (LGD) due to minimal risk of carcinogenesis, and significant procedural risks. Based on their promising results in prior validation studies targeting early detection of classical PC, DNA hypermethylation-based markers may serve as a biomarker for malignant risk stratification of IPMNs. This study investigates our DNA methylation-based PC biomarker panel (ADAMTS1, BNC1, and CACNA1G genes) in differentiating IPMN-advanced neoplasia from IPMN-LGDs. METHODS: Our previously described genome-wide pharmaco-epigenetic method identified multiple genes as potential targets for PC detection. The combination was further optimized and validated for early detection of classical PC in previous case-control studies. These promising genes were evaluated among micro-dissected IPMN tissue (IPMN-LGD: 35, IPMN-advanced neoplasia: 35) through Methylation-Specific PCR. The discriminant capacity of individual and combination of genes were delineated through Receiver Operating Characteristics curve analysis. RESULTS: As compared to IPMN-LGDs, IPMN-advanced neoplasia had higher hypermethylation frequency of candidate genes: ADAMTS1 (60% vs. 14%), BNC1 (66% vs. 3%), and CACGNA1G (25% vs. 0%). We observed Area Under Curve (AUC) values of 0.73 for ADAMTS1, 0.81 for BNC1, and 0.63 for CACNA1G genes. The combination of the BNC1/ CACNA1G genes resulted in an AUC of 0.84, sensitivity of 71%, and specificity of 97%. Combining the methylation status of the BNC1/CACNA1G genes, blood-based CA19-9, and IPMN lesion size enhanced the AUC to 0.92. CONCLUSION: DNA-methylation based biomarkers have shown a high diagnostic specificity and moderate sensitivity for differentiating IPMN-advanced neoplasia from LGDs. Addition of specific methylation targets can improve the accuracy of the methylation biomarker panel and enable the development of noninvasive IPMN stratification biomarkers.

DNA methylation of the promoter region at the CREB1 binding site is a mechanism for the epigenetic regulation of brain-specific PKMζ.

Protein kinase M zeta, PKMζ, is a brain enriched kinase with a well characterized role in Long-Term Potentiation (LTP), the activity-dependent strengthening of synapses involved in long-term memory formation. However, little is known about the molecular mechanisms that maintain the tissue specificity of this kinase. Here, we characterized the epigenetic factors, mainly DNA methylation, regulating PKMζ expression in the human brain. The PRKCZ gene has an upstream promoter regulating Protein kinase C ζ (PKCζ), and an internal promoter driving PKMζ expression. A demethylated region, including a canonical CREB binding site, situated at the internal promoter was only observed in human CNS tissues. The induction of site-specific hypermethylation of this region resulted in decreased CREB1 binding and downregulation of PKMζ expression. Noteworthy, CREB binding sites were absent in the upstream promoter of PRKCZ locus, suggesting a specific mechanism for regulating PKMζ expression. These observations were validated using a system of human neuronal differentiation from induced pluripotent stem cells (iPSCs). CREB1 binding at the internal promoter was detected only in differentiated neurons, where PKMζ is expressed. The same epigenetic mechanism in the context of CREB binding site was identified in other genes involved in neuronal differentiation and LTP. Additionally, aberrant DNA hypermethylation at the internal promoter was observed in cases of Alzheimer's disease, correlating with decreased expression of PKMζ in patient brains. Altogether, we present a conserved epigenetic mechanism regulating PKMζ expression and other genes enhanced in the CNS with possible implications in neuronal differentiation and Alzheimer's disease.

Age-dependent impairment in antibody responses elicited by a homologous CoronaVac booster dose.

The emergence of the SARS-CoV-2 Omicron sublineages resulted in increased transmission rates and reduced protection from vaccines. To counteract these effects, multiple booster strategies were used in different countries, although data comparing their efficiency in improving protective immunity remain sparse, especially among vulnerable populations, including older adults. The inactivated CoronaVac vaccine was among the most widely distributed vaccine worldwide and was essential in the early control of SARS-CoV-2-related hospitalizations and deaths. However, it is not well understood whether homologous versus heterologous booster doses in those fully vaccinated with CoronaVac induce distinct humoral responses or whether these responses vary across age groups. We analyzed plasma antibody responses from CoronaVac-vaccinated younger or older individuals who received a homologous CoronaVac or heterologous BNT162b2 or ChAdOx1 booster vaccine. All three evaluated boosters resulted in increased virus-specific IgG titers 28 days after the booster dose. However, we found that both IgG titers against SARS-CoV-2 Spike or RBD and neutralization titers against Omicron sublineages were substantially reduced in participants who received homologous CoronaVac compared with the heterologous BNT162b2 or ChAdOx1 booster. This effect was specifically prominent in recipients >50 years of age. In this group, the CoronaVac booster induced low virus-specific IgG titers and failed to elevate neutralization titers against any Omicron sublineage. Our results point to the notable inefficiency of CoronaVac immunization and boosting in mounting protective antiviral humoral immunity, particularly among older adults, during the Omicron wave. These observations also point to benefits of heterologous regimens in high-risk populations fully vaccinated with CoronaVac.

Rapid emergence of SARS-CoV-2 Omicron variant is associated with an infection advantage over Delta in vaccinated persons.

Background: The SARS-CoV-2 Omicron variant became a global concern due to its rapid spread and displacement of the dominant Delta variant. We hypothesized that part of Omicron's rapid rise was based on its increased ability to cause infections in persons that are vaccinated compared to Delta. Methods: We analyzed nasal swab PCR tests for samples collected between December 12 and 16, 2021, in Connecticut when the proportion of Delta and Omicron variants was relatively equal. We used the spike gene target failure (SGTF) to classify probable Delta and Omicron infections. We fitted an exponential curve to the estimated infections to determine the doubling times for each variant. We compared the test positivity rates for each variant by vaccination status, number of doses, and vaccine manufacturer. Generalized linear models were used to assess factors associated with odds of infection with each variant among persons testing positive for SARS-CoV-2. Findings: For infections with high virus copies (Ct < 30) among vaccinated persons, we found higher odds that they were infected with Omicron compared to Delta, and that the odds increased with increased number of vaccine doses. Compared to unvaccinated persons, we found significant reduction in Delta positivity rates after two (43.4%-49.1%) and three vaccine doses (81.1%), while we only found a significant reduction in Omicron positivity rates after three doses (62.3%). Conclusion: The rapid rise in Omicron infections was likely driven by Omicron's escape from vaccine-induced immunity. Funding: This work was supported by the Centers for Disease Control and Prevention (CDC).

Comparative transmissibility of SARS-CoV-2 variants Delta and Alpha in New England, USA.

The SARS-CoV-2 Delta variant rose to dominance in mid-2021, likely propelled by an estimated 40%-80% increased transmissibility over Alpha. To investigate if this ostensible difference in transmissibility is uniform across populations, we partner with public health programs from all six states in New England in the United States. We compare logistic growth rates during each variant's respective emergence period, finding that Delta emerged 1.37-2.63 times faster than Alpha (range across states). We compute variant-specific effective reproductive numbers, estimating that Delta is 63%-167% more transmissible than Alpha (range across states). Finally, we estimate that Delta infections generate on average 6.2 (95% CI 3.1-10.9) times more viral RNA copies per milliliter than Alpha infections during their respective emergence. Overall, our evidence suggests that Delta's enhanced transmissibility can be attributed to its innate ability to increase infectiousness, but its epidemiological dynamics may vary depending on underlying population attributes and sequencing data availability.

Genome-Wide Binding Analyses of HOXB1 Revealed a Novel DNA Binding Motif Associated with Gene Repression.

Knowledge of the diverse DNA binding specificities of transcription factors is important for understanding their specific regulatory functions in animal development and evolution. We have examined the genome-wide binding properties of the mouse HOXB1 protein in embryonic stem cells differentiated into neural fates. Unexpectedly, only a small number of HOXB1 bound regions (7%) correlate with binding of the known HOX cofactors PBX and MEIS. In contrast, 22% of the HOXB1 binding peaks display co-occupancy with the transcriptional repressor REST. Analyses revealed that co-binding of HOXB1 with PBX correlates with active histone marks and high levels of expression, while co-occupancy with REST correlates with repressive histone marks and repression of the target genes. Analysis of HOXB1 bound regions uncovered enrichment of a novel 15 base pair HOXB1 binding motif HB1RE (HOXB1 response element). In vitro template binding assays showed that HOXB1, PBX1, and MEIS can bind to this motif. In vivo, this motif is sufficient for direct expression of a reporter gene and over-expression of HOXB1 selectively represses this activity. Our analyses suggest that HOXB1 has evolved an association with REST in gene regulation and the novel HB1RE motif contributes to HOXB1 function in part through a repressive role in gene expression.