Integrating genomic and epidemiologic data to accelerate progress toward schistosomiasis elimination.

The global community has adopted ambitious goals to eliminate schistosomiasis as a public health problem, and new tools are needed to achieve them. Mass drug administration programs, for example, have reduced the burden of schistosomiasis, but the identification of hotspots of persistent and reemergent transmission threaten progress toward elimination and underscore the need to couple treatment with interventions that reduce transmission. Recent advances in DNA sequencing technologies make whole-genome sequencing a valuable and increasingly feasible option for population-based studies of complex parasites such as schistosomes. Here, we focus on leveraging genomic data to tailor interventions to distinct social and ecological circumstances. We consider two priority questions that can be addressed by integrating epidemiological, ecological, and genomic information: (1) how often do non-human host species contribute to human schistosome infection? and (2) what is the importance of locally acquired versus imported infections in driving transmission at different stages of elimination? These questions address processes that can undermine control programs, especially those that rely heavily on treatment with praziquantel. Until recently, these questions were difficult to answer with sufficient precision to inform public health decision-making. We review the literature related to these questions and discuss how whole-genome approaches can identify the geographic and taxonomic sources of infection, and how such information can inform context-specific efforts that advance schistosomiasis control efforts and minimize the risk of reemergence.

AP-2α and AP-2ß cooperatively function in the craniofacial surface ectoderm to regulate chromatin and gene expression dynamics during facial development.

The facial surface ectoderm is essential for normal development of the underlying cranial neural crest cell populations, providing signals that direct appropriate growth, patterning, and morphogenesis. Despite the importance of the ectoderm as a signaling center, the molecular cues and genetic programs implemented within this tissue are understudied. Here, we show that removal of two members of the AP-2 transcription factor family, AP-2α and AP-2ß, within the early embryonic ectoderm of the mouse leads to major alterations in the craniofacial complex. Significantly, there are clefts in both the upper face and mandible, accompanied by fusion of the upper and lower jaws in the hinge region. Comparison of ATAC-seq and RNA-seq analyses between controls and mutants revealed significant changes in chromatin accessibility and gene expression centered on multiple AP-2 binding motifs associated with enhancer elements within these ectodermal lineages. In particular, loss of these AP-2 proteins affects both skin differentiation as well as multiple signaling pathways, most notably the WNT pathway. We also determined that the mutant clefting phenotypes that correlated with reduced WNT signaling could be rescued by Wnt1 ligand overexpression in the ectoderm. Collectively, these findings highlight a conserved ancestral function for AP-2 transcription factors in ectodermal development and signaling, and provide a framework from which to understand the gene regulatory network operating within this tissue that directs vertebrate craniofacial development.

ProxECAT: Proxy External Controls Association Test. A new case-control gene region association test using allele frequencies from public controls.

A primary goal of the recent investment in sequencing is to detect novel genetic associations in health and disease improving the development of treatments and playing a critical role in precision medicine. While this investment has resulted in an enormous total number of sequenced genomes, individual studies of complex traits and diseases are often smaller and underpowered to detect rare variant genetic associations. Existing genetic resources such as the Exome Aggregation Consortium (>60,000 exomes) and the Genome Aggregation Database (~140,000 sequenced samples) have the potential to be used as controls in these studies. Fully utilizing these and other existing sequencing resources may increase power and could be especially useful in studies where resources to sequence additional samples are limited. However, to date, these large, publicly available genetic resources remain underutilized, or even misused, in large part due to the lack of statistical methods that can appropriately use this summary level data. Here, we present a new method to incorporate external controls in case-control analysis called ProxECAT (Proxy External Controls Association Test). ProxECAT estimates enrichment of rare variants within a gene region using internally sequenced cases and external controls. We evaluated ProxECAT in simulations and empirical analyses of obesity cases using both low-depth of coverage (7x) whole-genome sequenced controls and ExAC as controls. We find that ProxECAT maintains the expected type I error rate with increased power as the number of external controls increases. With an accompanying R package, ProxECAT enables the use of publicly available allele frequencies as external controls in case-control analysis.

Biliverdin Reductase B Dynamics Are Coupled to Coenzyme Binding.

Biliverdin reductase B (BLVRB) is a newly identified cellular redox regulator that catalyzes the NADPH-dependent reduction of multiple substrates. Through mass spectrometry analysis, we identified high levels of BLVRB in mature red blood cells, highlighting the importance of BLVRB in redox regulation. The BLVRB conformational changes that occur during conezyme/substrate binding and the role of dynamics in BLVRB function, however, remain unknown. Through a combination of NMR, kinetics, and isothermal titration calorimetry studies, we determined that BLVRB binds its coenzyme 500-fold more tightly than its substrate. While the active site of apo BLVRB is highly dynamic on multiple timescales, active site dynamics are largely quenched within holo BLVRB, in which dynamics are redistributed to other regions of the enzyme. We show that a single point mutation of Arg78➔Ala leads to both an increase in active site micro-millisecond motions and an increase in the microscopic rate constants of coenzyme binding. This demonstrates that altering BLVRB active site dynamics can directly cause a change in functional characteristics. Our studies thus address the solution behavior of apo and holo BLVRB and identify a role of enzyme dynamics in coenzyme binding.