Stage-specific dynamic reorganization of genome topology shapes transcriptional neighborhoods in developing human retinal organoids.

We have generated a high-resolution Hi-C map of developing human retinal organoids to elucidate spatiotemporal dynamics of genomic architecture and its relationship with gene expression patterns. We demonstrate progressive stage-specific alterations in DNA topology and correlate these changes with transcription of cell-type-restricted gene markers during retinal differentiation. Temporal Hi-C reveals a shift toward A compartment for protein-coding genes and B compartment for non-coding RNAs, displaying high and low expression, respectively. Notably, retina-enriched genes are clustered near lost boundaries of topologically associated domains (TADs), and higher-order assemblages (i.e., TAD cliques) localize in active chromatin regions with binding sites for eye-field transcription factors. These genes gain chromatin contacts at their transcription start site as organoid differentiation proceeds. Our study provides a global view of chromatin architecture dynamics associated with diversification of cell types during retinal development and serves as a foundational resource for in-depth functional investigations of retinal developmental traits.

PLAE Web App Enables Powerful Searching and Multiple Visualizations Across One Million Unified Single-Cell Ocular Transcriptomes.

Purpose: To create a high-performance reactive web application to query single-cell gene expression data across cell type, species, study, and other factors. Methods: We updated the content and structure of the underlying data (single cell Eye in a Disk [scEiaD]) and wrote the web application PLAE (https://plae.nei.nih.gov) to visualize and explore the data. Results: The new portal provides quick visualization of over a million individual cells from vertebrate eye and body transcriptomes encompassing four species, 60 cell types, six ocular tissues, and 23 body tissues across 35 publications. To demonstrate the value of this unified pan-eye dataset, we replicated known neurogenic and cone macula markers in addition to proposing six new cone human region markers. Conclusions: The PLAE web application offers the eye community a powerful and quick means to test hypotheses related to gene expression across a highly diverse, community-derived database. Translational Relevance: The PLAE resource enables any researcher or clinician to study and research gene expression patterning across a wide variety of curated ocular cell types with a responsive web app.

Protein kinase CK2 modulates the activity of Maf-family bZIP transcription factor NRL in rod photoreceptors of mammalian retina.

Maf-family basic motif leucine zipper protein NRL specifies rod photoreceptor cell fate during retinal development and, in concert with homeodomain protein CRX and other regulatory factors, controls the expression of most rod-expressed genes including the visual pigment gene Rhodopsin (Rho). Transcriptional regulatory activity of NRL is modulated by post-translational modifications, especially phosphorylation, and mutations at specific phosphosites can lead to retinal degeneration. During our studies to elucidate NRL-mediated transcriptional regulation, we identified protein kinase CK2 in NRL-enriched complexes bound to Rho promoter-enhancer regions and in NRL-enriched high molecular mass fractions from the bovine retina. The presence of CK2 in NRL complexes was confirmed by co-immunoprecipitation from developing and adult mouse retinal extracts. In vitro kinase assay and bioinformatic analysis indicated phosphorylation of NRL at Ser117 residue by CK2. Co-transfection of Csnk2a1 cDNA encoding murine CK2 with human NRL and CRX reduced the bovine Rho promoter-driven luciferase expression in HEK293 cells and mutagenesis of NRL-Ser117 residue to Ala restored the reporter gene activity. In concordance, overexpression of CK2 in the mouse retina in vivo by electroporation resulted in reduction of Rho promoter-driven DsRed reporter expression as well as the transcript level of many phototransduction genes. Thus, our studies demonstrate that CK2 can phosphorylate Ser117 of NRL. Modulation of NRL activity by CK2 suggests intricate interdependence of transcriptional and signaling pathways in maintaining rod homeostasis.

High-resolution genome topology of human retina uncovers super enhancer-promoter interactions at tissue-specific and multifactorial disease loci.

Chromatin organization and enhancer-promoter contacts establish unique spatiotemporal gene expression patterns in distinct cell types. Non-coding genetic variants can influence cellular phenotypes by modifying higher-order transcriptional hubs and consequently gene expression. To elucidate genomic regulation in human retina, we mapped chromatin contacts at high resolution and integrated with super-enhancers (SEs), histone marks, binding of CTCF and select transcription factors. We show that topologically associated domains (TADs) with central SEs exhibit stronger insulation and augmented contact with retinal genes relative to TADs with edge SEs. Merging genome-wide expression quantitative trait loci (eQTLs) with topology map reveals physical links between 100 eQTLs and corresponding eGenes associated with retinal neurodegeneration. Additionally, we uncover candidate genes for susceptibility variants linked to age-related macular degeneration and glaucoma. Our study of high-resolution genomic architecture of human retina provides insights into genetic control of tissue-specific functions, suggests paradigms for missing heritability, and enables the dissection of common blinding disease phenotypes.

Nicotinamide Promotes Formation of Retinal Organoids From Human Pluripotent Stem Cells via Enhanced Neural Cell Fate Commitment.

Retinal organoids (ROs) derived from human pluripotent stem cells (hPSCs) recapitulate key features of retinogenesis and provide a promising platform to study retinal development and disease in a human context. Although multiple protocols are currently in use, hPSCs exhibit tremendous variability in differentiation efficiency, with some cell lines consistently yielding few or even no ROs, limiting their utility in research. We report here that early nicotinamide (NAM) treatment significantly improves RO yield across 8 hPSC lines from different donors, including some that would otherwise fail to generate a meaningful number of ROs. NAM treatment promotes neural commitment of hPSCs at the expense of non-neural ectodermal cell fate, which in turn increases eye field progenitor generation. Further analysis suggests that this effect is partially mediated through inhibition of BMP signaling. Our data encourage a broader use of human ROs for disease modeling applications that require the use of multiple patient-specific cell lines.

A novel exome probe set captures phototransduction genes across birds (Aves) enabling efficient analysis of vision evolution.

The diversity of avian visual phenotypes provides a framework for studying mechanisms of trait diversification generally, and the evolution of vertebrate vision, specifically. Previous research has focused on opsins, but to fully understand visual adaptation, we must study the complete phototransduction cascade (PTC). Here, we developed a probe set that captures exonic regions of 46 genes representing the PTC and other light responses. For a subset of species, we directly compared gene capture between our probe set and low-coverage whole genome sequencing (WGS), and we discuss considerations for choosing between these methods. Finally, we developed a unique strategy to avoid chimeric assembly by using "decoy" reference sequences. We successfully captured an average of 64% of our targeted exome in 46 species across 14 orders using the probe set and had similar recovery using the WGS data. Compared to WGS or transcriptomes, our probe set: (1) reduces sequencing requirements by efficiently capturing vision genes, (2) employs a simpler bioinformatic pipeline by limiting required assembly and negating annotation, and (3) eliminates the need for fresh tissues, enabling researchers to leverage existing museum collections. We then utilized our vision exome data to identify positively selected genes in two evolutionary scenarios-evolution of night vision in nocturnal birds and evolution of high-speed vision specific to manakins (Pipridae). We found parallel positive selection of SLC24A1 in both scenarios, implicating the alteration of rod response kinetics, which could improve color discrimination in dim light conditions and/or facilitate higher temporal resolution.

Ecological mechanism of climate-mediated selection in a rapidly evolving invasive species.

Recurring seasonal changes can lead to the evolution of phenological cues. For example, many arthropods undergo photoperiodic diapause, a programmed developmental arrest induced by short autumnal day length. The selective mechanisms that determine the timing of autumnal diapause initiation have not been empirically identified. We quantified latitudinal clines in genetically determined diapause timing of an invasive mosquito, Aedes albopictus, on two continents. We show that variation in diapause timing within and between continents is explained by a novel application of a growing degree day (GDD) model that delineates a location-specific deadline after which it is not possible to complete an additional full life cycle. GDD models are widely used to predict spring phenology by modelling growth and development as physiological responses to ambient temperatures. Our results show that the energy accumulation dynamics represented by GDD models have also led to the evolution of an anticipatory life-history cue in autumn.

Gene Therapy of Dominant CRX-Leber Congenital Amaurosis using Patient Stem Cell-Derived Retinal Organoids.

Mutations in the photoreceptor transcription factor gene cone-rod homeobox (CRX) lead to distinct retinopathy phenotypes, including early-onset vision impairment in dominant Leber congenital amaurosis (LCA). Using induced pluripotent stem cells (iPSCs) from a patient with CRX-I138fs48 mutation, we established an in vitro model of CRX-LCA in retinal organoids that showed defective photoreceptor maturation by histology and gene profiling, with diminished expression of visual opsins. Adeno-associated virus (AAV)-mediated CRX gene augmentation therapy partially restored photoreceptor phenotype and expression of phototransduction-related genes as determined by single-cell RNA-sequencing. Retinal organoids derived from iPSCs of a second dominant CRX-LCA patient carrying K88N mutation revealed the loss of opsin expression as a common phenotype, which was alleviated by AAV-mediated augmentation of CRX. Our studies provide a proof-of-concept for developing gene therapy of dominant CRX-LCA and other CRX retinopathies.

Rapid adaptive evolution of the diapause program during range expansion of an invasive mosquito.

In temperate climates, the recurring seasonal exigencies of winter represent a fundamental physiological challenge for a wide range of organisms. In response, many temperate insects enter diapause, an alternative developmental program, including developmental arrest, that allows organisms to synchronize their life cycle with seasonal environmental variation. Geographic variation in diapause phenology contributing to local climatic adaptation is well documented. However, few studies have examined how the rapid evolution of a suite of traits expressed across the diapause program may contribute to climatic adaptation on a contemporary timescale. Here, we investigate the evolution of the diapause program over the past 35 years by leveraging a "natural experiment" presented by the recent invasion of the Asian tiger mosquito, Aedes albopictus, across the eastern United States. We sampled populations from two distinct climatic regions separated by 6° of latitude (∼700 km). Using common-garden experiments, we identified regional genetic divergence in diapause-associated cold tolerance, diapause duration, and postdiapause starvation tolerance. We also found regional divergence in nondiapause thermal performance. In contrast, we observed minimal regional divergence in nondiapause larval growth traits and at neutral molecular marker loci. Our results demonstrate rapid evolution of the diapause program and imply strong selection caused by differences in winter conditions.

Juvenile Hormone III but Not 20-Hydroxyecdysone Regulates the Embryonic Diapause of Aedes albopictus.

Diapause is an alternative developmental trajectory allowing insects to enter dormancy and persist through predictable periods of seasonally unfavorable conditions. This crucial ecological adaptation defines the geographic and seasonal abundance of many insect pollinators, pests, and vectors. Understanding the hormonal changes by which insects coordinate the perception of external, diapause-inducing cues with the physiological mechanisms that lead to developmental arrest is a long-standing goal in biology. The hormonal regulation of diapause tends to vary by the life stage at which development arrest occurs; for example, diapause is typically regulated by ecdysteroids in larvae and pupae, and by juvenile hormones in adults. However, little is known about the hormonal control of embryonic diapause, particularly in Diptera. To address this fundamental gap, we directly measured 20-hydroxyecdysone (20HE) (via LC-MS/MS) and juvenile hormone III (JH3) (via GC-MS) in diapause and non-diapause eggs of the Asian tiger mosquito, Aedes albopictus. While 20HE abundance did not differ, diapause eggs had lower JH3 abundance than non-diapause eggs. These results are corroborated by transcriptional and manipulative evidence suggesting that reduced JH3 regulates diapause in this medically important mosquito.

Diapause-associated changes in the lipid and metabolite profiles of the Asian tiger mosquito, Aedes albopictus.

Diapause is an alternative life-history strategy that allows organisms to enter developmental arrest in anticipation of unfavorable conditions. Diapause is widespread among insects and plays a key role in enhancing overwinter survival as well as defining the seasonal and geographic distributions of populations. Next-generation sequencing has greatly advanced our understanding of the transcriptional basis for this crucial adaptation but less is known about the regulation of embryonic diapause physiology at the metabolite level. Here, we characterized the lipid and metabolite profiles of embryonic diapause in the Asian tiger mosquito, Aedes albopictus We used an untargeted approach to capture the relative abundance of 250 lipids and 241 metabolites. We observed adjustments associated with increased energy storage, including an accumulation of lipids, the formation of larger lipid droplets and increased lipogenesis, as well as metabolite shifts suggesting reduced energy utilization. We also found changes in neuroregulatory- and insulin-associated metabolites with potential roles in diapause regulation. Finally, we detected a group of unidentified, diapause-specific metabolites which have physical properties similar to those of steroids/steroid derivatives and may be associated with the ecdysteroidal regulation of embryonic diapause in A.albopictus Together, these results deepen our understanding of the metabolic regulation of embryonic diapause and identify key targets for future investigations.

Guanine Limitation Results in CodY-Dependent and -Independent Alteration of Staphylococcus aureus Physiology and Gene Expression.

In Staphylococcus aureus, the global transcriptional regulator CodY modulates the expression of hundreds of genes in response to the availability of GTP and the branched-chain amino acids isoleucine, leucine, and valine (ILV). CodY DNA-binding activity is high when GTP and ILV are abundant. When GTP and ILV are limited, CodY's affinity for DNA drops, altering expression of CodY-regulated targets. In this work, we investigated the impact of guanine nucleotides (GNs) on S. aureus physiology and CodY activity by constructing a guaA null mutant (ΔguaA strain). De novo biosynthesis of guanine monophosphate is abolished due to the guaA mutation; thus, the mutant cells require exogenous guanosine for growth. We also found that CodY activity was reduced when we knocked out guaA, activating the Agr two-component system and increasing secreted protease activity. Notably, in a rich, complex medium, we detected an increase in alternative sigma factor B activity in the ΔguaA mutant, which results in a 5-fold increase in production of the antioxidant pigment staphyloxanthin. Under biologically relevant flow conditions, ΔguaA cells failed to form robust biofilms when limited for guanine or guanosine. Transcriptome sequencing (RNA-Seq) analysis of the S. aureus transcriptome during growth in guanosine-limited chemostats revealed substantial CodY-dependent and -independent alterations of gene expression profiles. Importantly, these changes increase production of proteases and δ-toxin, suggesting that S. aureus exhibits a more invasive lifestyle when limited for guanosine. Further, gene products upregulated under GN limitation, including those necessary for lipoic acid biosynthesis and sugar transport, may prove to be useful drug targets for treating Gram-positive infections.IMPORTANCEStaphylococcus aureus infections impose a serious economic burden on health care facilities and patients because of the emergence of strains resistant to last-line antibiotics. Understanding the physiological processes governing fitness and virulence of S. aureus in response to environmental cues is critical for developing efficient diagnostics and treatments. De novo purine biosynthesis is essential for both fitness and virulence in S. aureus since inhibiting production cripples S. aureus's ability to cause infection. Here, we corroborate these findings and show that blocking guanine nucleotide synthesis severely affects S. aureus fitness by altering metabolic and virulence gene expression. Characterizing pathways and gene products upregulated in response to guanine limitation can aid in the development of novel adjuvant strategies to combat S. aureus infections.

Helping struggling students in introductory biology: a peer-tutoring approach that improves performance, perception, and retention.

The high attrition rate among science, technology, engineering, and mathematics (STEM) majors has long been an area of concern for institutions and educational researchers. The transition from introductory to advanced courses has been identified as a particularly "leaky" point along the STEM pipeline, and students who struggle early in an introductory STEM course are predominantly at risk. Peer-tutoring programs offered to all students in a course have been widely found to help STEM students during this critical transition, but hiring a sufficient number of tutors may not be an option for some institutions. As an alternative, this study examines the viability of an optional peer-tutoring program offered to students who are struggling in a large-enrollment, introductory biology course. Struggling students who regularly attended peer tutoring increased exam performance, expert-like perceptions of biology, and course persistence relative to their struggling peers who were not attending the peer-tutoring sessions. The results of this study provide information to instructors who want to design targeted academic assistance for students who are struggling in introductory courses.