Adaptive sequence divergence forged new neurodevelopmental enhancers in humans.

Searches for the genetic underpinnings of uniquely human traits have focused on human-specific divergence in conserved genomic regions, which reflects adaptive modifications of existing functional elements. However, the study of conserved regions excludes functional elements that descended from previously neutral regions. Here, we demonstrate that the fastest-evolved regions of the human genome, which we term "human ancestor quickly evolved regions" (HAQERs), rapidly diverged in an episodic burst of directional positive selection prior to the human-Neanderthal split, before transitioning to constraint within hominins. HAQERs are enriched for bivalent chromatin states, particularly in gastrointestinal and neurodevelopmental tissues, and genetic variants linked to neurodevelopmental disease. We developed a multiplex, single-cell in vivo enhancer assay to discover that rapid sequence divergence in HAQERs generated hominin-unique enhancers in the developing cerebral cortex. We propose that a lack of pleiotropic constraints and elevated mutation rates poised HAQERs for rapid adaptation and subsequent susceptibility to disease.

Correcting signal biases and detecting regulatory elements in STARR-seq data.

High-throughput reporter assays such as self-transcribing active regulatory region sequencing (STARR-seq) have made it possible to measure regulatory element activity across the entire human genome at once. The resulting data, however, present substantial analytical challenges. Here, we identify technical biases that explain most of the variance in STARR-seq data. We then develop a statistical model to correct those biases and to improve detection of regulatory elements. This approach substantially improves precision and recall over current methods, improves detection of both activating and repressive regulatory elements, and controls for false discoveries despite strong local correlations in signal.

Handheld electronic device use in patient care: the emergency department patient perspective-a cross-sectional survey.

BACKGROUND: Staff use of smartphones and tablets in the healthcare setting is increasingly prevalent, but little is known about whether this use is acceptable to patients. Staff are concerned that the use of handheld electronic devices (HEDs) may be negatively misconstrued by patients. The HED can be a valuable tool, offering the emergency clinician access to a wealth of resources; it is therefore vital that patient views are addressed during their widespread adoption into clinical practice. METHODS: Patients, or those accompanying them, within the ED of the Royal Derby Hospital between April and June 2017 were asked to complete a survey consisting of 22 questions. Data collection took place to include all times of day and every day of the week. Every eligible individual within the department during a data collection period was approached. RESULTS: A total of 438 respondents successfully completed the survey with a response rate of 92%. Only 2% of those who observed staff using HEDs during their ED visit thought that they were being used for non-clinical purposes. 339 (78%) agreed that staff should be allowed to use HEDs in the workplace. Concerns expressed by respondents included devices being used for non-clinical purposes and data security. The main suggestion by respondents was that the purpose of the HEDs should be explained to patients to avoid misinterpretation. CONCLUSION: Our survey shows that the majority of survey respondents felt that clinical staff should be allowed to use HEDs in the workplace and that many of the concerns raised could be addressed with adequate patient information and clear governance.

Bayesian estimation of genetic regulatory effects in high-throughput reporter assays.

MOTIVATION: High-throughput reporter assays dramatically improve our ability to assign function to noncoding genetic variants, by measuring allelic effects on gene expression in the controlled setting of a reporter gene. Unlike genetic association tests, such assays are not confounded by linkage disequilibrium when loci are independently assayed. These methods can thus improve the identification of causal disease mutations. While work continues on improving experimental aspects of these assays, less effort has gone into developing methods for assessing the statistical significance of assay results, particularly in the case of rare variants captured from patient DNA. RESULTS: We describe a Bayesian hierarchical model, called Bayesian Inference of Regulatory Differences, which integrates prior information and explicitly accounts for variability between experimental replicates. The model produces substantially more accurate predictions than existing methods when allele frequencies are low, which is of clear advantage in the search for disease-causing variants in DNA captured from patient cohorts. Using the model, we demonstrate a clear tradeoff between variant sequencing coverage and numbers of biological replicates, and we show that the use of additional biological replicates decreases variance in estimates of effect size, due to the properties of the Poisson-binomial distribution. We also provide a power and sample size calculator, which facilitates decision making in experimental design parameters. AVAILABILITY AND IMPLEMENTATION: The software is freely available from www.geneprediction.org/bird. The experimental design web tool can be accessed at http://67.159.92.22:8080. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Human genome-wide measurement of drug-responsive regulatory activity.

Environmental stimuli commonly act via changes in gene regulation. Human-genome-scale assays to measure such responses are indirect or require knowledge of the transcription factors (TFs) involved. Here, we present the use of human genome-wide high-throughput reporter assays to measure environmentally-responsive regulatory element activity. We focus on responses to glucocorticoids (GCs), an important class of pharmaceuticals and a paradigmatic genomic response model. We assay GC-responsive regulatory activity across >108 unique DNA fragments, covering the human genome at >50×. Those assays directly detected thousands of GC-responsive regulatory elements genome-wide. We then validate those findings with measurements of transcription factor occupancy, histone modifications, chromatin accessibility, and gene expression. We also detect allele-specific environmental responses. Notably, the assays did not require knowledge of GC response mechanisms. Thus, this technology can be used to agnostically quantify genomic responses for which the underlying mechanism remains unknown.

Evaluation of chromatin accessibility in prefrontal cortex of individuals with schizophrenia.

Schizophrenia genome-wide association studies have identified >150 regions of the genome associated with disease risk, yet there is little evidence that coding mutations contribute to this disorder. To explore the mechanism of non-coding regulatory elements in schizophrenia, we performed ATAC-seq on adult prefrontal cortex brain samples from 135 individuals with schizophrenia and 137 controls, and identified 118,152 ATAC-seq peaks. These accessible chromatin regions in the brain are highly enriched for schizophrenia SNP heritability. Accessible chromatin regions that overlap evolutionarily conserved regions exhibit an even higher heritability enrichment, indicating that sequence conservation can further refine functional risk variants. We identify few differences in chromatin accessibility between cases and controls, in contrast to thousands of age-related differential accessible chromatin regions. Altogether, we characterize chromatin accessibility in the human prefrontal cortex, the effect of schizophrenia and age on chromatin accessibility, and provide evidence that our dataset will allow for fine mapping of risk variants.

Nuclease Footprints in Sperm Project Past and Future Chromatin Regulatory Events.

Nuclear remodeling to a condensed state is a hallmark of spermatogenesis. This is achieved by replacement of histones with protamines. Regions retaining nucleosomes may be of functional significance. To determine their potential roles, sperm from wild type and transgenic mice harboring a single copy insert of the human protamine cluster were subjected to Micrococcal Nuclease-seq. CENTIPEDE, a hierarchical Bayesian model, was used to identify multiple spatial patterns, "footprints", of MNase-seq reads along the sperm genome. Regions predicted by CENTIPEDE analysis to be bound by a regulatory factor in sperm were correlated with genomic landmarks and higher order chromatin structure datasets to identify potential roles for these factors in regulating either prior or post spermatogenic, i.e., early embryonic events. This approach linked robust endogenous protamine transcription and transgene suppression to its chromatin environment within topologically associated domains. Of the candidate enhancer-bound regulatory proteins, Ctcf, was associated with chromatin domain boundaries in testes and embryonic stem cells. The continuity of Ctcf binding through the murine germline may permit rapid reconstitution of chromatin organization following fertilization. This likely reflects its preparation for early zygotic genome activation and comparatively accelerated preimplantation embryonic development program observed in mouse as compared to human and bull.

Oxidative DNA damage in mouse sperm chromosomes: Size matters.

Normal embryo and foetal development as well as the health of the progeny are mostly dependent on gamete nuclear integrity. In the present study, in order to characterize more precisely oxidative DNA damage in mouse sperm we used two mouse models that display high levels of sperm oxidative DNA damage, a common alteration encountered both in in vivo and in vitro reproduction. Immunoprecipitation of oxidized sperm DNA coupled to deep sequencing showed that mouse chromosomes may be largely affected by oxidative alterations. We show that the vulnerability of chromosomes to oxidative attack inversely correlated with their size and was not linked to their GC richness. It was neither correlated with the chromosome content in persisting nucleosomes nor associated with methylated sequences. A strong correlation was found between oxidized sequences and sequences rich in short interspersed repeat elements (SINEs). Chromosome position in the sperm nucleus as revealed by fluorescent in situ hybridization appears to be a confounder. These data map for the first time fragile mouse sperm chromosomal regions when facing oxidative damage that may challenge the repair mechanisms of the oocyte post-fertilization.

Chromatin and extracellular vesicle associated sperm RNAs.

A diverse pool of RNAs remain encapsulated within the transcriptionally silent spermatozoon despite the dramatic reduction in cellular and nuclear volume following cytoplasm/nucleoplasm expulsion. The impact of this pronounced restructuring on the distribution of transcripts inside the sperm essentially remains unknown. To define their compartmentalization, total RNA >100 nt was extracted from sonicated (SS) mouse spermatozoa and detergent demembranated sucrose gradient fractionated (Cs/Tx) sperm heads. Sperm RNAs predominately localized toward the periphery. The corresponding distribution of transcripts and thus localization and complexity were then inferred by RNA-seq. Interestingly, the number of annotated RNAs in the CsTx sperm heads exhibiting reduced peripheral enrichment was restricted. However this included Cabyr, the calcium-binding tyrosine phosphorylation-regulated protein encoded transcript. It is present in murine zygotes prior to the maternal to the zygotic transition yet absent in oocytes, consistent with the delivery of internally positioned sperm-borne RNAs to the embryo. In comparison, transcripts enriched in sonicated sperm contributed to the mitochondria and exosomes along with several nuclear transcripts including the metastasis associated lung adenocarcinoma transcript 1 (Malat1) and several small nucleolar RNAs. Their preferential peripheral localization suggests that chromatin remodeling during spermiogenesis is not limited to nucleoproteins as part of the nucleoprotein exchange.

Stability, delivery and functions of human sperm RNAs at fertilization.

Increasing attention has focused on the significance of RNA in sperm, in light of its contribution to the birth and long-term health of a child, role in sperm function and diagnostic potential. As the composition of sperm RNA is in flux, assigning specific roles to individual RNAs presents a significant challenge. For the first time RNA-seq was used to characterize the population of coding and non-coding transcripts in human sperm. Examining RNA representation as a function of multiple methods of library preparation revealed unique features indicative of very specific and stage-dependent maturation and regulation of sperm RNA, illuminating their various transitional roles. Correlation of sperm transcript abundance with epigenetic marks suggested roles for these elements in the pre- and post-fertilization genome. Several classes of non-coding RNAs including lncRNAs, CARs, pri-miRNAs, novel elements and mRNAs have been identified which, based on factors including relative abundance, integrity in sperm, available knockout data of embryonic effect and presence or absence in the unfertilized human oocyte, are likely to be essential male factors critical to early post-fertilization development. The diverse and unique attributes of sperm transcripts that were revealed provides the first detailed analysis of the biology and anticipated clinical significance of spermatozoal RNAs.

Local and global factors affecting RNA sequencing analysis.

High-throughput RNA sequencing (RNA-seq) continues to provide unparalleled insight into transcriptome complexity. Now the "gold standard" for assessing global transcript levels, RNA-seq is poised to revolutionize our understanding of transcription and posttranscriptional regulation of RNA. Despite significant advantages over prior experimental strategies, RNA-seq is not without pitfalls. We have identified a number of confounding factors that significantly affect sequencing coverage. These include regional GC content, preferential sites of fragmentation, and read "pile-up" due to primer affinity and transcript end effects. Independent of cell type and laboratory, when ignored, these factors can bias analyses. Understanding the underlying principles responsible for producing these artifacts is key to recognizing both their presence and how their effects may be controlled. Here we outline the causes of and strategies to avoid several previously unreported complicating factors common to RNA-seq experiments.

Interrogating the transgenic genome: development of an interspecies tiling array.

A single expressing copy of the human protamine domain was randomly inserted into an intron of Cyp2c38. The transgenic locus was shown to recapitulate the level of expression observed in normal human testis while not perturbing endogenous protamine expression. The development of an interspecies tiling array was pursued to enable direct comparison of the orthologous protamine domains in a single experiment. Probe design was adapted to generate species-specific high resolution probe sets that would tolerate repetitive elements. Results from competitive hybridizations demonstrate that interspecies tiling arrays are a valuable tool for parallel analysis of highly similar DNA sequences. This approach provides a rapid and reliable means of interrogating samples prior to deep sequencing analysis. These arrays should readily compliment most DNA isolation and analysis techniques such as ChIP, nuclease sensitivity and nuclear matrix association assays.

The sperm nucleus: chromatin, RNA, and the nuclear matrix.

Within the sperm nucleus, the paternal genome remains functionally inert and protected following protamination. This is marked by a structural morphogenesis that is heralded by a striking reduction in nuclear volume. Despite these changes, both human and mouse spermatozoa maintain low levels of nucleosomes that appear non-randomly distributed throughout the genome. These regions may be necessary for organizing higher order genomic structure through interactions with the nuclear matrix. The promoters of this transcriptionally quiescent genome are differentially marked by modified histones that may poise downstream epigenetic effects. This notion is supported by increasing evidence that the embryo inherits these differing levels of chromatin organization. In concert with the suite of RNAs retained in the mature sperm, they may synergistically interact to direct early embryonic gene expression. Irrespective, these features reflect the transcriptional history of spermatogenic differentiation. As such, they may soon be utilized as clinical markers of male fertility. In this review, we explore and discuss how this may be orchestrated.

Identification of human sperm transcripts as candidate markers of male fertility.

To date, there has been little progress towards identifying markers of normal male fertility. The need to supplement current subjective methods that rely on variable semen parameters to assess fertility status continues to be acknowledged. Several studies have shown that spermatozoal RNAs can describe characteristic failures of the spermatogenic pathway among infertile males. In spite of the inherent heterogeneity of semen that describe the "normal" fertile male, this holds the promise of developing markers that could help identify the ever elusive idiopathic infertile male. Through the analyses of the spermatozoal transcriptome from 24 donors of proven fertility, we identified a series of transcripts that were consistently present among all individuals. The heterogeneous nature of the samples, reflected by their semen parameters, was mirrored by the variability of the observed array signal. Nevertheless, clusters of invariable transcript pairs were identified. These were founded by a single central member that was linked in constant proportion even though the absolute level of each member of the transcript pair often varied among individuals. The presence of pairs of stable transcripts suggests that among the heterogeneity observed in the sperm transcriptome, a distinct set is strictly regulated.

Real-time PCR quantification using a variable reaction efficiency model.

The quantitative real-time polymerase chain reaction (PCR) remains a cornerstone technique in gene expression analysis and sequence characterization. Despite the importance of the approach to experimental biology, the confident assignment of reaction efficiency to the early cycles of real-time PCR reactions remains problematic. Considerable noise may be generated when few cycles in the amplification are available to estimate peak efficiency. An alternate approach that uses data from beyond the log-linear amplification phase is explored in this article with the aim of reducing noise and adding confidence to efficiency estimates. PCR reaction efficiency is regressed to estimate the per-cycle profile of an asymptotically departed peak efficiency even when this is not closely approximated in the measurable cycles. The process can be repeated over replicates to develop a robust estimate of peak reaction efficiency. This leads to an estimate of the maximum reaction efficiency that may be considered primer design specific. Using a series of biological scenarios, we demonstrate that this approach can provide an accurate estimate of initial template concentration.