Single-cell roadmap of human gonadal development.

Gonadal development is a complex process that involves sex determination followed by divergent maturation into either testes or ovaries1. Historically, limited tissue accessibility, a lack of reliable in vitro models and critical differences between humans and mice have hampered our knowledge of human gonadogenesis, despite its importance in gonadal conditions and infertility. Here, we generated a comprehensive map of first- and second-trimester human gonads using a combination of single-cell and spatial transcriptomics, chromatin accessibility assays and fluorescent microscopy. We extracted human-specific regulatory programmes that control the development of germline and somatic cell lineages by profiling equivalent developmental stages in mice. In both species, we define the somatic cell states present at the time of sex specification, including the bipotent early supporting population that, in males, upregulates the testis-determining factor SRY and sPAX8s, a gonadal lineage located at the gonadal-mesonephric interface. In females, we resolve the cellular and molecular events that give rise to the first and second waves of granulosa cells that compartmentalize the developing ovary to modulate germ cell differentiation. In males, we identify human SIGLEC15+ and TREM2+ fetal testicular macrophages, which signal to somatic cells outside and inside the developing testis cords, respectively. This study provides a comprehensive spatiotemporal map of human and mouse gonadal differentiation, which can guide in vitro gonadogenesis.

Cells of the human intestinal tract mapped across space and time.

The cellular landscape of the human intestinal tract is dynamic throughout life, developing in utero and changing in response to functional requirements and environmental exposures. Here, to comprehensively map cell lineages, we use single-cell RNA sequencing and antigen receptor analysis of almost half a million cells from up to 5 anatomical regions in the developing and up to 11 distinct anatomical regions in the healthy paediatric and adult human gut. This reveals the existence of transcriptionally distinct BEST4 epithelial cells throughout the human intestinal tract. Furthermore, we implicate IgG sensing as a function of intestinal tuft cells. We describe neural cell populations in the developing enteric nervous system, and predict cell-type-specific expression of genes associated with Hirschsprung's disease. Finally, using a systems approach, we identify key cell players that drive the formation of secondary lymphoid tissue in early human development. We show that these programs are adopted in inflammatory bowel disease to recruit and retain immune cells at the site of inflammation. This catalogue of intestinal cells will provide new insights into cellular programs in development, homeostasis and disease.

Flexible comparison of batch correction methods for single-cell RNA-seq using BatchBench.

As the cost of single-cell RNA-seq experiments has decreased, an increasing number of datasets are now available. Combining newly generated and publicly accessible datasets is challenging due to non-biological signals, commonly known as batch effects. Although there are several computational methods available that can remove batch effects, evaluating which method performs best is not straightforward. Here, we present BatchBench (https://github.com/cellgeni/batchbench), a modular and flexible pipeline for comparing batch correction methods for single-cell RNA-seq data. We apply BatchBench to eight methods, highlighting their methodological differences and assess their performance and computational requirements through a compendium of well-studied datasets. This systematic comparison guides users in the choice of batch correction tool, and the pipeline makes it easy to evaluate other datasets.

Cationic Geminoid Peptide Amphiphiles Inhibit DENV2 Protease, Furin, and Viral Replication.

Dengue is an important arboviral infectious disease for which there is currently no specific cure. We report gemini-like (geminoid) alkylated amphiphilic peptides containing lysines in combination with glycines or alanines (C15H31C(O)-Lys-(Gly or Ala)nLys-NHC16H33, shorthand notation C16-KXnK-C16 with X = A or G, and n = 0-2). The representatives with 1 or 2 Ala inhibit dengue protease and human furin, two serine proteases involved in dengue virus infection that have peptides with cationic amino acids as their preferred substrates, with IC50 values in the lower µM range. The geminoid C16-KAK-C16 combined inhibition of DENV2 protease (IC50 2.3 µM) with efficacy against replication of wildtype DENV2 in LLC-MK2 cells (EC50 4.1 µM) and an absence of toxicity. We conclude that the lysine-based geminoids have activity against dengue virus infection, which is based on their inhibition of the proteases involved in viral replication and are therefore promising leads to further developing antiviral therapeutics, not limited to dengue.

Transposon-driven transcription is a conserved feature of vertebrate spermatogenesis and transcript evolution.

Spermatogenesis is associated with major and unique changes to chromosomes and chromatin. Here, we sought to understand the impact of these changes on spermatogenic transcriptomes. We show that long terminal repeats (LTRs) of specific mouse endogenous retroviruses (ERVs) drive the expression of many long non-coding transcripts (lncRNA). This process occurs post-mitotically predominantly in spermatocytes and round spermatids. We demonstrate that this transposon-driven lncRNA expression is a conserved feature of vertebrate spermatogenesis. We propose that transposon promoters are a mechanism by which the genome can explore novel transcriptional substrates, increasing evolutionary plasticity and allowing for the genesis of novel coding and non-coding genes. Accordingly, we show that a small fraction of these novel ERV-driven transcripts encode short open reading frames that produce detectable peptides. Finally, we find that distinct ERV elements from the same subfamilies act as differentially activated promoters in a tissue-specific context. In summary, we demonstrate that LTRs can act as tissue-specific promoters and contribute to post-mitotic spermatogenic transcriptome diversity.

Large-scale analysis of microRNA expression, epi-transcriptomic features and biogenesis.

MicroRNAs are important genetic regulators in both animals and plants. They have a range of functions spanning development, differentiation, growth, metabolism and disease. The advent of next-generation sequencing technologies has made it a relatively straightforward task to detect these molecules and their relative expression via sequencing. There are a large number of published studies with deposited datasets. However, there are currently few resources that capitalize on these data to better understand the features, distribution and biogenesis of miRNAs. Herein, we focus on Human and Mouse for which the majority of data are available. We reanalyse sequencing data from 461 samples into a coordinated catalog of microRNA expression. We use this to perform large-scale analyses of miRNA function and biogenesis. These analyses include global expression comparison, co-expression of miRNA clusters and the prediction of miRNA strand-specificity and underlying constraints. Additionally, we report for the first time a global analysis of miRNA epi-transcriptomic modifications and assess their prevalence across tissues, samples and families. Finally, we report a list of potentially mis-annotated miRNAs in miRBase based on their aggregated modification profiles. The results have been collated into a comprehensive online repository of miRNA expression and features such as modifications and RNA editing events, which is available at: http://wwwdev.ebi.ac.uk/enright-dev/miratlas. We believe these findings will further contribute to our understanding of miRNA function in animals and benefit the miRNA community in general.

Mirnovo: genome-free prediction of microRNAs from small RNA sequencing data and single-cells using decision forests.

The discovery of microRNAs (miRNAs) remains an important problem, particularly given the growth of high-throughput sequencing, cell sorting and single cell biology. While a large number of miRNAs have already been annotated, there may well be large numbers of miRNAs that are expressed in very particular cell types and remain elusive. Sequencing allows us to quickly and accurately identify the expression of known miRNAs from small RNA-Seq data. The biogenesis of miRNAs leads to very specific characteristics observed in their sequences. In brief, miRNAs usually have a well-defined 5' end and a more flexible 3' end with the possibility of 3' tailing events, such as uridylation. Previous approaches to the prediction of novel miRNAs usually involve the analysis of structural features of miRNA precursor hairpin sequences obtained from genome sequence. We surmised that it may be possible to identify miRNAs by using these biogenesis features observed directly from sequenced reads, solely or in addition to structural analysis from genome data. To this end, we have developed mirnovo, a machine learning based algorithm, which is able to identify known and novel miRNAs in animals and plants directly from small RNA-Seq data, with or without a reference genome. This method performs comparably to existing tools, however is simpler to use with reduced run time. Its performance and accuracy has been tested on multiple datasets, including species with poorly assembled genomes, RNaseIII (Drosha and/or Dicer) deficient samples and single cells (at both embryonic and adult stage).

Extent, causes, and consequences of small RNA expression variation in human adipose tissue.

Small RNAs are functional molecules that modulate mRNA transcripts and have been implicated in the aetiology of several common diseases. However, little is known about the extent of their variability within the human population. Here, we characterise the extent, causes, and effects of naturally occurring variation in expression and sequence of small RNAs from adipose tissue in relation to genotype, gene expression, and metabolic traits in the MuTHER reference cohort. We profiled the expression of 15 to 30 base pair RNA molecules in subcutaneous adipose tissue from 131 individuals using high-throughput sequencing, and quantified levels of 591 microRNAs and small nucleolar RNAs. We identified three genetic variants and three RNA editing events. Highly expressed small RNAs are more conserved within mammals than average, as are those with highly variable expression. We identified 14 genetic loci significantly associated with nearby small RNA expression levels, seven of which also regulate an mRNA transcript level in the same region. In addition, these loci are enriched for variants significant in genome-wide association studies for body mass index. Contrary to expectation, we found no evidence for negative correlation between expression level of a microRNA and its target mRNAs. Trunk fat mass, body mass index, and fasting insulin were associated with more than twenty small RNA expression levels each, while fasting glucose had no significant associations. This study highlights the similar genetic complexity and shared genetic control of small RNA and mRNA transcripts, and gives a quantitative picture of small RNA expression variation in the human population.

Functional interlocked systems.

With the advent of supramolecular chemistry and later nanotechnology a great deal of research has been focused on new types of molecular structures, which are not held together by covalent bonds but by non-covalent mechanical interactions. Examples include the catenane, rotaxane, and knot interlocked structures. The design and synthesis of these architectures is an art by itself and as such is worth being reviewed. In this tutorial review we will focus, however, on the functional aspects of interlocked molecules and discuss how these can find applications, e.g. as artificial muscles, as molecular valves, as components of electronic devices, and as catalysts.

Structure-delivery relationships of lysine-based gemini surfactants and their lipoplexes.

The synthesis and properties of gemini surfactants of the type (R(1)(CO)-Lys(H)-NH)2(CH2)n are reported. For a spacer length of n = 6, the hydrophobic acyl tail was varied in length (R(1) = C8, C10, C12, C14, C16, and C18) and, for R(1) = C18, the degree of unsaturation. For R(1)(CO) = oleoyl (C18:1 Z) the spacer length (n = 2-8) and the stereochemistry of the lysine building block were varied; a 'half-gemini' derivative with a single oleoyl tail and head group was also prepared. The potential of the gemini surfactants to transfer polynucleotides across a cell membrane was investigated by transfection of HeLa cells with beta-galactosidase, both in the presence and absence of the helper lipid DOPE. Oleoyl was found to be by far the best hydrophobic tail for this biological activity, whereas the effect of the lysine stereochemistry was less pronounced. The effect of an optimum spacer length (n = 6) was observed only in the absence of helper lipid. The most active surfactant, i.e. the one with oleoyl chains and n = 6, formed liposomes with sizes in the range of 60-350 nm, and its lipoplex underwent a transition from a lamellar to a hexagonal morphology upon lowering the pH from 7 to 3.

Depletion of HPV16 early genes induces autophagy and senescence in a cervical carcinogenesis model, regardless of viral physical state.

In cervical carcinomas, high-risk human papillomavirus (HR-HPV) may be integrated into host chromosomes or remain extra-chromosomal (episomal). We used the W12 cervical keratinocyte model to investigate the effects of HPV16 early gene depletion on in vitro cervical carcinogenesis pathways, particularly effects shared by cells with episomal versus integrated HPV16 DNA. Importantly, we were able to study the specific cellular consequences of viral gene depletion by using short interfering RNAs known not to cause phenotypic or transcriptional off-target effects in keratinocytes. We found that while cervical neoplastic progression in vitro was characterized by dynamic changes in HPV16 transcript levels, viral early gene expression was required for cell survival at all stages of carcinogenesis, regardless of viral physical state, levels of early gene expression or histology in organotypic tissue culture. Moreover, HPV16 early gene depletion induced changes in host gene expression that were common to both episome-containing and integrant-containing cells. In particular, we observed up-regulation of autophagy genes, associated with enrichment of senescence and innate immune-response pathways, including the senescence-associated secretory phenotype (SASP). In keeping with these observations, HPV16 early gene depletion induced autophagy in both episome-containing and integrant-containing W12 cells, as evidenced by the appearance of autophagosomes, punctate expression of the autophagy marker LC3, conversion of LC3B-I to LC3B-II, and reduced levels of the autophagy substrate p62. Consistent with the reported association between autophagy and senescence pathways, HPV16 early gene depletion induced expression of the senescence marker beta-galactosidase and increased secretion of the SASP-related protein IGFBP3. Together, these data indicate that depleting HR-HPV early genes would be of potential therapeutic benefit in all cervical carcinogenesis pathways, regardless of viral physical state. In addition, the senescence/SASP response associated with autophagy induction may promote beneficial immune effects in bystander cells.

Kraken: a set of tools for quality control and analysis of high-throughput sequence data.

New sequencing technologies pose significant challenges in terms of data complexity and magnitude. It is essential that efficient software is developed with performance that scales with this growth in sequence information. Here we present a comprehensive and integrated set of tools for the analysis of data from large scale sequencing experiments. It supports adapter detection and removal, demultiplexing of barcodes, paired-end data, a range of read architectures and the efficient removal of sequence redundancy. Sequences can be trimmed and filtered based on length, quality and complexity. Quality control plots track sequence length, composition and summary statistics with respect to genomic annotation. Several use cases have been integrated into a single streamlined pipeline, including both mRNA and small RNA sequencing experiments. This pipeline interfaces with existing tools for genomic mapping and differential expression analysis.

Processive catalysis.

Nature's enzymes are an ongoing source of inspiration for scientists. The complex processes behind their selectivity and efficiency is slowly being unraveled, and these findings have spawned many biomimetic catalysts. However, nearly all focus on the conversion of small molecular substrates. Nature itself is replete with inventive catalytic systems which modify, replicate, or decompose entire polymers, often in a processive fashion. Such processivity can, for example, enhance the rate of catalysis by clamping to the polymer substrate, which imparts a large effective molarity. Reviewed herein are the various strategies for processivity in nature's arsenal and their properties. An overview of what has been achieved by chemists aiming to mimic one of nature's greatest tricks is also included.

Transcriptional signals of transformation in human cancer.

BACKGROUND: As normal cells transform into cancers, their cell state changes, which may drive cancer cells into a stem-like or more primordial, foetal, or embryonic cell state. The transcriptomic profile of this final state may encode information about cancer's origin and how cancers relate to their normal cell counterparts. METHODS: Here, we used single-cell atlases to study cancer transformation in transcriptional terms. We utilised bulk transcriptomes across a wide spectrum of adult and childhood cancers, using a previously established method to interrogate their relationship to normal cell states. We extend and validate these findings using single-cell cancer transcriptomes and organ-specific atlases of colorectal and liver cancer. RESULTS: Our bulk transcriptomic data reveals that adult cancers rarely return to an embryonic state, but that a foetal state is a near-universal feature of childhood cancers. This finding was confirmed with single-cell cancer transcriptomes. CONCLUSIONS: Our findings provide a nuanced picture of transformation in human cancer, indicating cancer-specific rather than universal patterns of transformation pervade adult epithelial cancers.

Targeted deletion of microRNA-22 promotes stress-induced cardiac dilation and contractile dysfunction.

BACKGROUND: Delineating the role of microRNAs (miRNAs) in the posttranscriptional gene regulation offers new insights into how the heart adapts to pathological stress. We developed a knockout of miR-22 in mice and investigated its function in the heart. METHODS AND RESULTS: Here, we show that miR-22-deficient mice are impaired in inotropic and lusitropic response to acute stress by dobutamine. Furthermore, the absence of miR-22 sensitized mice to cardiac decompensation and left ventricular dilation after long-term stimulation by pressure overload. Calcium transient analysis revealed reduced sarcoplasmic reticulum Ca(2+) load in association with repressed sarcoplasmic reticulum Ca(2+) ATPase activity in mutant myocytes. Genetic ablation of miR-22 also led to a decrease in cardiac expression levels for Serca2a and muscle-restricted genes encoding proteins in the vicinity of the cardiac Z disk/titin cytoskeleton. These phenotypes were attributed in part to inappropriate repression of serum response factor activity in stressed hearts. Global analysis revealed increased expression of the transcriptional/translational repressor purine-rich element binding protein B, a highly conserved miR-22 target implicated in the negative control of muscle expression. CONCLUSION: These data indicate that miR-22 functions as an integrator of Ca(2+) homeostasis and myofibrillar protein content during stress in the heart and shed light on the mechanisms that enhance propensity toward heart failure.

Long- and short-term selective forces on malaria parasite genomes.

Plasmodium parasites, the causal agents of malaria, result in more than 1 million deaths annually. Plasmodium are unicellular eukaryotes with small ∼23 Mb genomes encoding ∼5200 protein-coding genes. The protein-coding genes comprise about half of these genomes. Although evolutionary processes have a significant impact on malaria control, the selective pressures within Plasmodium genomes are poorly understood, particularly in the non-protein-coding portion of the genome. We use evolutionary methods to describe selective processes in both the coding and non-coding regions of these genomes. Based on genome alignments of seven Plasmodium species, we show that protein-coding, intergenic and intronic regions are all subject to purifying selection and we identify 670 conserved non-genic elements. We then use genome-wide polymorphism data from P. falciparum to describe short-term selective processes in this species and identify some candidate genes for balancing (diversifying) selection. Our analyses suggest that there are many functional elements in the non-genic regions of these genomes and that adaptive evolution has occurred more frequently in the protein-coding regions of the genome.

Yolk sac cell atlas reveals multiorgan functions during human early development.

The extraembryonic yolk sac (YS) ensures delivery of nutritional support and oxygen to the developing embryo but remains ill-defined in humans. We therefore assembled a comprehensive multiomic reference of the human YS from 3 to 8 postconception weeks by integrating single-cell protein and gene expression data. Beyond its recognized role as a site of hematopoiesis, we highlight roles in metabolism, coagulation, vascular development, and hematopoietic regulation. We reconstructed the emergence and decline of YS hematopoietic stem and progenitor cells from hemogenic endothelium and revealed a YS-specific accelerated route to macrophage production that seeds developing organs. The multiorgan functions of the YS are superseded as intraembryonic organs develop, effecting a multifaceted relay of vital functions as pregnancy proceeds.

Single-cell Atlas of common variable immunodeficiency shows germinal center-associated epigenetic dysregulation in B-cell responses.

Common variable immunodeficiency (CVID), the most prevalent symptomatic primary immunodeficiency, displays impaired terminal B-cell differentiation and defective antibody responses. Incomplete genetic penetrance and ample phenotypic expressivity in CVID suggest the participation of additional pathogenic mechanisms. Monozygotic (MZ) twins discordant for CVID are uniquely valuable for studying the contribution of epigenetics to the disease. Here, we generate a single-cell epigenomics and transcriptomics census of naïve-to-memory B cell differentiation in a CVID-discordant MZ twin pair. Our analysis identifies DNA methylation, chromatin accessibility and transcriptional defects in memory B-cells mirroring defective cell-cell communication upon activation. These findings are validated in a cohort of CVID patients and healthy donors. Our findings provide a comprehensive multi-omics map of alterations in naïve-to-memory B-cell transition in CVID and indicate links between the epigenome and immune cell cross-talk. Our resource, publicly available at the Human Cell Atlas, gives insight into future diagnosis and treatments of CVID patients.

Myeloid cell interferon secretion restricts Zika flavivirus infection of developing and malignant human neural progenitor cells.

Zika virus (ZIKV) can infect human developing brain (HDB) progenitors resulting in epidemic microcephaly, whereas analogous cellular tropism offers treatment potential for the adult brain cancer, glioblastoma (GBM). We compared productive ZIKV infection in HDB and GBM primary tissue explants that both contain SOX2+ neural progenitors. Strikingly, although the HDB proved uniformly vulnerable to ZIKV infection, GBM was more refractory, and this correlated with an innate immune expression signature. Indeed, GBM-derived CD11b+ microglia/macrophages were necessary and sufficient to protect progenitors against ZIKV infection in a non-cell autonomous manner. Using SOX2+ GBM cell lines, we found that CD11b+-conditioned medium containing type 1 interferon beta (IFNß) promoted progenitor resistance to ZIKV, whereas inhibition of JAK1/2 signaling restored productive infection. Additionally, CD11b+ conditioned medium, and IFNß treatment rendered HDB progenitor lines and explants refractory to ZIKV. These findings provide insight into neuroprotection for HDB progenitors as well as enhanced GBM oncolytic therapies.

Detecting microRNA binding and siRNA off-target effects from expression data.

Sylamer is a method for detecting microRNA target and small interfering RNA off-target signals in 3' untranslated regions from a ranked gene list, sorted from upregulated to downregulated, after a microRNA perturbation or RNA interference experiment. The output is a landscape plot that tracks occurrence biases using hypergeometric P-values for all words across the gene ranking. We demonstrated the utility, speed and accuracy of this approach on several datasets.