A High-Quality Blue Whale Genome, Segmental Duplications, and Historical Demography.

The blue whale, Balaenoptera musculus, is the largest animal known to have ever existed, making it an important case study in longevity and resistance to cancer. To further this and other blue whale-related research, we report a reference-quality, long-read-based genome assembly of this fascinating species. We assembled the genome from PacBio long reads and utilized Illumina/10×, optical maps, and Hi-C data for scaffolding, polishing, and manual curation. We also provided long read RNA-seq data to facilitate the annotation of the assembly by NCBI and Ensembl. Additionally, we annotated both haplotypes using TOGA and measured the genome size by flow cytometry. We then compared the blue whale genome with other cetaceans and artiodactyls, including vaquita (Phocoena sinus), the world's smallest cetacean, to investigate blue whale's unique biological traits. We found a dramatic amplification of several genes in the blue whale genome resulting from a recent burst in segmental duplications, though the possible connection between this amplification and giant body size requires further study. We also discovered sites in the insulin-like growth factor-1 gene correlated with body size in cetaceans. Finally, using our assembly to examine the heterozygosity and historical demography of Pacific and Atlantic blue whale populations, we found that the genomes of both populations are highly heterozygous and that their genetic isolation dates to the last interglacial period. Taken together, these results indicate how a high-quality, annotated blue whale genome will serve as an important resource for biology, evolution, and conservation research.

Enhancing biological signals and detection rates in single-cell RNA-seq experiments with cDNA library equalization.

Considerable effort has been devoted to refining experimental protocols to reduce levels of technical variability and artifacts in single-cell RNA-sequencing data (scRNA-seq). We here present evidence that equalizing the concentration of cDNA libraries prior to pooling, a step not consistently performed in single-cell experiments, improves gene detection rates, enhances biological signals, and reduces technical artifacts in scRNA-seq data. To evaluate the effect of equalization on various protocols, we developed Scaffold, a simulation framework that models each step of an scRNA-seq experiment. Numerical experiments demonstrate that equalization reduces variation in sequencing depth and gene-specific expression variability. We then performed a set of experiments in vitro with and without the equalization step and found that equalization increases the number of genes that are detected in every cell by 17-31%, improves discovery of biologically relevant genes, and reduces nuisance signals associated with cell cycle. Further support is provided in an analysis of publicly available data.

Ronin is essential for embryogenesis and the pluripotency of mouse embryonic stem cells.

Pluripotency is a unique biological state that allows cells to differentiate into any tissue type. Here we describe a candidate pluripotency factor, Ronin, that possesses a THAP domain, which is associated with sequence-specific DNA binding and epigenetic silencing of gene expression. Ronin is expressed primarily during the earliest stages of murine embryonic development, and its deficiency in mice produces periimplantational lethality and defects in the inner cell mass. Conditional knockout of Ronin prevents the growth of ES cells while forced expression of Ronin allows ES cells to proliferate without differentiation under conditions that normally do not promote self-renewal. Ectopic expression also partly compensates for the effects of Oct4 knockdown. We demonstrate that Ronin binds directly to HCF-1, a key transcriptional regulator. Our findings identify Ronin as an essential factor underlying embryogenesis and ES cell pluripotency. Its association with HCF-1 suggests an epigenetic mechanism of gene repression in pluripotent cells.

Rapid changes in chromatin structure during dedifferentiation of primary hepatocytes in vitro.

Primary hepatocytes are widely used in the pharmaceutical industry to screen drug candidates for hepatotoxicity, but hepatocytes quickly dedifferentiate and lose their mature metabolic function in culture. Attempts have been made to better recapitulate the in vivo liver environment in culture, but the full spectrum of signals required to maintain hepatocyte function ex vivo remains elusive. To elucidate molecular changes that accompany, and may contribute to dedifferentiation of hepatocytes ex vivo, we performed lineage tracing and comprehensive profiling of alterations in their gene expression profiles and chromatin landscape during culture. First, using genetically tagged hepatocytes we demonstrate that expression of the fetal gene alpha-fetoprotein in cultured hepatocytes comes from cells that previously expressed the mature gene albumin, and not from a population of albumin-negative precursor cells, proving mature hepatocytes undergo true dedifferentiation in culture. Next we studied the dedifferentiation process in detail through bulk RNA-sequencing of hepatocytes cultured over an extended period. We identified three distinct phases of dedifferentiation: an early phase, where mature hepatocyte genes are rapidly downregulated in a matter of hours; a middle phase, where fetal genes are activated; and a late phase, where initially rare contaminating non-parenchymal cells proliferate, taking over the culture. Lastly, to better understand the signaling events that result in the rapid downregulation of mature genes in hepatocytes, we examined changes in chromatin accessibility in these cells during the first 24 h of culture using Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq). We find that drastic and rapid changes in chromatin accessibility occur immediately upon the start of culture. Using binding motif analysis of the areas of open chromatin sharing similar temporal profiles, we identify several candidate transcription factors potentially involved in the dedifferentiation of primary hepatocytes in culture.

In vivo replacement of damaged bladder urothelium by Wolffian duct epithelial cells.

The bladder's remarkable regenerative capacity had been thought to derive exclusively from its own progenitors. While examining consequences of DNA methyltransferase 1 (Dnmt1) inactivation in mouse embryonic bladder epithelium, we made the surprising discovery that Wolffian duct epithelial cells can support bladder regeneration. Conditional Dnmt1 inactivation in mouse urethral and bladder epithelium triggers widespread apoptosis, depletes basal and intermediate bladder cells, and disrupts uroplakin protein expression. These events coincide with Wolffian duct epithelial cell recruitment into Dnmt1 mutant urethra and bladder where they are reprogrammed to express bladder markers, including FOXA1, keratin 5, P63, and uroplakin. This is evidence that Wolffian duct epithelial cells are summoned in vivo to replace damaged bladder epithelium and function as a reservoir of cells for bladder regeneration.

Epithelial DNA methyltransferase-1 regulates cell survival, growth and maturation in developing prostatic buds.

DNA methyltransferase 1 (DNMT1) is required for embryogenesis but roles in late forming organ systems including the prostate, which emerges from the urethral epithelium, have not been fully examined. We used a targeted genetic approach involving a Shhcre recombinase to demonstrate requirement of epithelial DNA methyltransferase-1 (Dnmt1) in mouse prostate morphogenesis. Dnmt1 mutant urethral cells exhibit DNA hypomethylation, DNA damage, p53 accumulation and undergo cell cycle arrest and apoptosis. Urethral epithelial cells are disorganized in Dnmt1 mutants, leading to impaired prostate growth and maturation and failed glandular development. We evaluated oriented cell division as a mechanism of bud elongation and widening by demonstrating that mitotic spindle axes typically form parallel or perpendicular to prostatic bud elongation axes. We then deployed a ShhcreERT allele to delete Dnmt1 from a subset of urethral epithelial cells, creating mosaic mutants with which to interrogate the requirement for cell division in specific prostatic bud epithelial populations. DNMT1- cell distribution within prostatic buds is not random as would be expected in a process where DNMT1 was not required. Instead, replication competent DNMT1 + cells primarily accumulate in prostatic bud margins and tips while replication impeded DNMT1- cells accumulate in prostatic bud cores. Together, these results highlight the role of DNMT1 in regulating epithelial bud formation by maintaining cell cycle progression and survival of rapidly dividing urethral epithelial cells, which can be extended to the study of other developing epithelial organs. In addition, our results show that prostatic buds consist of two epithelial cell populations with distinct molecular and functional characteristics that could potentially contribute to specialized lineages in the adult prostate.

SCnorm: robust normalization of single-cell RNA-seq data.

The normalization of RNA-seq data is essential for accurate downstream inference, but the assumptions upon which most normalization methods are based are not applicable in the single-cell setting. Consequently, applying existing normalization methods to single-cell RNA-seq data introduces artifacts that bias downstream analyses. To address this, we introduce SCnorm for accurate and efficient normalization of single-cell RNA-seq data.

Functional characterization of human pluripotent stem cell-derived arterial endothelial cells.

Here, we report the derivation of arterial endothelial cells from human pluripotent stem cells that exhibit arterial-specific functions in vitro and in vivo. We combine single-cell RNA sequencing of embryonic mouse endothelial cells with an EFNB2-tdTomato/EPHB4-EGFP dual reporter human embryonic stem cell line to identify factors that regulate arterial endothelial cell specification. The resulting xeno-free protocol produces cells with gene expression profiles, oxygen consumption rates, nitric oxide production levels, shear stress responses, and TNFα-induced leukocyte adhesion rates characteristic of arterial endothelial cells. Arterial endothelial cells were robustly generated from multiple human embryonic and induced pluripotent stem cell lines and have potential applications for both disease modeling and regenerative medicine.

Spatial patterns of gene expression are unveiled in the chick primitive streak by ordering single-cell transcriptomes.

During vertebrate development, progenitor cells give rise to tissues and organs through a complex choreography that commences at gastrulation. A hallmark event of gastrulation is the formation of the primitive streak, a linear assembly of cells along the anterior-posterior (AP) axis of the developing organism. To examine the primitive streak at a single-cell resolution, we measured the transcriptomes of individual chick cells from the streak or the surrounding tissue (the rest of the area pellucida) in Hamburger-Hamilton stage 4 embryos. The single-cell transcriptomes were then ordered by the statistical method Wave-Crest to deduce both the relative position along the AP axis and the prospective lineage of single cells. The ordered transcriptomes reveal intricate patterns of gene expression along the primitive streak.

A folic acid-enriched diet attenuates prostate involution in response to androgen deprivation.

BACKGROUND: Serum folate concentrations in the United States have risen since dietary folic acid fortification was first mandated in 1998. Although maternal folic acid offers protection against neural tube defects in conceptuses, its impact on other organ systems and life stages have not been fully examined. Here, we used a mouse model to investigate the impact of a Folic acid (FA) enriched diet on prostate homeostasis and response to androgen deprivation. METHODS: Male mice were fed a control diet (4 mg FA/kg feed) or a folic acid supplemented diet (24 mg FA/kg feed) beginning at conception and continuing through early adulthood, when mice were castrated. RESULTS: We made the surprising observation that dietary FA supplementation confers partial resistance to castration-mediated prostate involution. At 3, 10, and 14 days post-castration, FA enriched diet fed mice had larger prostates as assessed by wet weight, taller prostatic luminal epithelial cells, and more abundant RNAs encoding prostate secretory proteins than castrated control diet fed mice. Diet did not significantly affect prostate weights of intact mice or serum testosterone concentrations of castrated mice. RNA-Seq analysis revealed that the FA enriched diet was associated with a unique prostate gene expression signature, affecting several signaling and metabolic pathways. CONCLUSIONS: Continuous exposure to a FA enriched diet slows prostate involution in response to androgen deprivation. Prostates from FA diet mice have increased secretory gene expression and increased luminal cell heights. The influence of dietary FA supplementation on the prostate response to androgen deprivation raises a future need to consider how dietary folic acid supplementation affects efficacy of androgen-reducing therapies for treating prostate disease.

A Novel Approach to Single Cell RNA-Sequence Analysis Facilitates In Silico Gene Reporting of Human Pluripotent Stem Cell-Derived Retinal Cell Types.

Cell type-specific investigations commonly use gene reporters or single-cell analytical techniques. However, reporter line development is arduous and generally limited to a single gene of interest, while single-cell RNA (scRNA)-sequencing (seq) frequently yields equivocal results that preclude definitive cell identification. To examine gene expression profiles of multiple retinal cell types derived from human pluripotent stem cells (hPSCs), we performed scRNA-seq on optic vesicle (OV)-like structures cultured under cGMP-compatible conditions. However, efforts to apply traditional scRNA-seq analytical methods based on unbiased algorithms were unrevealing. Therefore, we developed a simple, versatile, and universally applicable approach that generates gene expression data akin to those obtained from reporter lines. This method ranks single cells by expression level of a bait gene and searches the transcriptome for genes whose cell-to-cell rank order expression most closely matches that of the bait. Moreover, multiple bait genes can be combined to refine datasets. Using this approach, we provide further evidence for the authenticity of hPSC-derived retinal cell types. Stem Cells 2018;36:313-324.

Trendy: segmented regression analysis of expression dynamics in high-throughput ordered profiling experiments.

BACKGROUND: High-throughput expression profiling experiments with ordered conditions (e.g. time-course or spatial-course) are becoming more common for studying detailed differentiation processes or spatial patterns. Identifying dynamic changes at both the individual gene and whole transcriptome level can provide important insights about genes, pathways, and critical time points. RESULTS: We present an R package, Trendy, which utilizes segmented regression models to simultaneously characterize each gene's expression pattern and summarize overall dynamic activity in ordered condition experiments. For each gene, Trendy finds the optimal segmented regression model and provides the location and direction of dynamic changes in expression. We demonstrate the utility of Trendy to provide biologically relevant results on both microarray and RNA-sequencing (RNA-seq) datasets. CONCLUSIONS: Trendy is a flexible R package which characterizes gene-specific expression patterns and summarizes changes of global dynamics over ordered conditions. Trendy is freely available on Bioconductor with a full vignette at https://bioconductor.org/packages/release/bioc/html/Trendy.html .

Oscope identifies oscillatory genes in unsynchronized single-cell RNA-seq experiments.

Oscillatory gene expression is fundamental to development, but technologies for monitoring expression oscillations are limited. We have developed a statistical approach called Oscope to identify and characterize the transcriptional dynamics of oscillating genes in single-cell RNA-seq data from an unsynchronized cell population. Applying Oscope to a number of data sets, we demonstrated its utility and also identified a potential artifact in the Fluidigm C1 platform.

OEFinder: a user interface to identify and visualize ordering effects in single-cell RNA-seq data.

UNLABELLED: A recent article identified an artifact in multiple single-cell RNA-seq (scRNA-seq) datasets generated by the Fluidigm C1 platform. Specifically, Leng et al. showed significantly increased gene expression in cells captured from sites with small or large plate output IDs. We refer to this artifact as an ordering effect (OE). Including OE genes in downstream analyses could lead to biased results. To address this problem, we developed a statistical method and software called OEFinder to identify a sorted list of OE genes. OEFinder is available as an R package along with user-friendly graphical interface implementations which allows users to check for potential artifacts in scRNA-seq data generated by the Fluidigm C1 platform. AVAILABILITY AND IMPLEMENTATION: OEFinder is freely available at https://github.com/lengning/OEFinder CONTACT: rstewart@morgridge.org or lengning1@gmail.com SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Efficient genome engineering in human pluripotent stem cells using Cas9 from Neisseria meningitidis.

Genome engineering in human pluripotent stem cells (hPSCs) holds great promise for biomedical research and regenerative medicine. Recently, an RNA-guided, DNA-cleaving interference pathway from bacteria [the type II clustered, regularly interspaced, short palindromic repeats (CRISPR)-CRISPR-associated (Cas) pathway] has been adapted for use in eukaryotic cells, greatly facilitating genome editing. Only two CRISPR-Cas systems (from Streptococcus pyogenes and Streptococcus thermophilus), each with their own distinct targeting requirements and limitations, have been developed for genome editing thus far. Furthermore, limited information exists about homology-directed repair (HDR)-mediated gene targeting using long donor DNA templates in hPSCs with these systems. Here, using a distinct CRISPR-Cas system from Neisseria meningitidis, we demonstrate efficient targeting of an endogenous gene in three hPSC lines using HDR. The Cas9 RNA-guided endonuclease from N. meningitidis (NmCas9) recognizes a 5'-NNNNGATT-3' protospacer adjacent motif (PAM) different from those recognized by Cas9 proteins from S. pyogenes and S. thermophilus (SpCas9 and StCas9, respectively). Similar to SpCas9, NmCas9 is able to use a single-guide RNA (sgRNA) to direct its activity. Because of its distinct protospacer adjacent motif, the N. meningitidis CRISPR-Cas machinery increases the sequence contexts amenable to RNA-directed genome editing.

Comparative RNA-seq analysis in the unsequenced axolotl: the oncogene burst highlights early gene expression in the blastema.

The salamander has the remarkable ability to regenerate its limb after amputation. Cells at the site of amputation form a blastema and then proliferate and differentiate to regrow the limb. To better understand this process, we performed deep RNA sequencing of the blastema over a time course in the axolotl, a species whose genome has not been sequenced. Using a novel comparative approach to analyzing RNA-seq data, we characterized the transcriptional dynamics of the regenerating axolotl limb with respect to the human gene set. This approach involved de novo assembly of axolotl transcripts, RNA-seq transcript quantification without a reference genome, and transformation of abundances from axolotl contigs to human genes. We found a prominent burst in oncogene expression during the first day and blastemal/limb bud genes peaking at 7 to 14 days. In addition, we found that limb patterning genes, SALL genes, and genes involved in angiogenesis, wound healing, defense/immunity, and bone development are enriched during blastema formation and development. Finally, we identified a category of genes with no prior literature support for limb regeneration that are candidates for further evaluation based on their expression pattern during the regenerative process.

Establishment of Transgene-Free Porcine Induced Pluripotent Stem Cells.

Although protocols to generate authentic transgene-free mouse and human induced pluripotent stem cells (iPSCs) are now well established, standard methods for reprogramming porcine somatic cells still suffer from low efficiency and transgene retention. The Basic Protocol describes reprogramming procedures to establish transgene-free porcine iPSCs (PiPSCs) from porcine fibroblasts. This method uses episomal plasmids encoding POU5F1, SOX2, NANOG, KLF4, SV40LT, c-MYC, LIN28A, and microRNA-302/367, combined with an optimized medium, to establish PiPSC lines. Support protocols describe the establishment and characterization of clonal PiPSC lines, as well as the preparation of feeder cells and EBNA1 mRNA. This optimized, step-by-step approach tailored to this species enables the efficient derivation of PiPSCs in ∼4 weeks. The establishment of transgene-free PiPSCs provides a new and valuable model for studies of larger mammalian species' development, disease, and regenerative biology. © 2024 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol: Reprogramming of porcine fibroblasts with episomal plasmids Support Protocol 1: Preparation of mouse embryonic fibroblasts for feeder layer Support Protocol 2: Preparation of in vitro-transcribed EBNA1 mRNA Support Protocol 3: Establishment of clonal porcine induced pluripotent stem cell (PiPSC) lines Support Protocol 4: PiPSC characterization: Genomic DNA PCR and RT-PCR Support Protocol 5: PiPSC characterization: Immunostaining.

Chromosome level genome assembly of the Etruscan shrew Suncus etruscus.

Suncus etruscus is one of the world's smallest mammals, with an average body mass of about 2 grams. The Etruscan shrew's small body is accompanied by a very high energy demand and numerous metabolic adaptations. Here we report a chromosome-level genome assembly using PacBio long read sequencing, 10X Genomics linked short reads, optical mapping, and Hi-C linked reads. The assembly is partially phased, with the 2.472 Gbp primary pseudohaplotype and 1.515 Gbp alternate. We manually curated the primary assembly and identified 22 chromosomes, including X and Y sex chromosomes. The NCBI genome annotation pipeline identified 39,091 genes, 19,819 of them protein-coding. We also identified segmental duplications, inferred GO term annotations, and computed orthologs of human and mouse genes. This reference-quality genome will be an important resource for research on mammalian development, metabolism, and body size control.

Species-specific roles of the Notch ligands, receptors, and targets orchestrating the signaling landscape of the segmentation clock.

Somitogenesis is a hallmark feature of all vertebrates and some invertebrate species that involves the periodic formation of block-like structures called somites. Somites are transient embryonic segments that eventually establish the entire vertebral column. A highly conserved molecular oscillator called the segmentation clock underlies this periodic event and the pace of this clock regulates the pace of somite formation. Although conserved signaling pathways govern the clock in most vertebrates, the mechanisms underlying the species-specific divergence in various clock characteristics remain elusive. For example, the segmentation clock in classical model species such as zebrafish, chick, and mouse embryos tick with a periodicity of ∼30, ∼90, and ∼120 min respectively. This enables them to form the species-specific number of vertebrae during their overall timespan of somitogenesis. Here, we perform a systematic review of the species-specific features of the segmentation clock with a keen focus on mouse embryos. We perform this review using three different perspectives: Notch-responsive clock genes, ligand-receptor dynamics, and synchronization between neighboring oscillators. We further review reports that use non-classical model organisms and in vitro model systems that complement our current understanding of the segmentation clock. Our review highlights the importance of comparative developmental biology to further our understanding of this essential developmental process.

Efficient derivation of transgene-free porcine induced pluripotent stem cells enables in vitro modeling of species-specific developmental timing.

Sus scrofa domesticus (pig) has served as a superb large mammalian model for biomedical studies because of its comparable physiology and organ size to humans. The derivation of transgene-free porcine induced pluripotent stem cells (PiPSCs) will, therefore, benefit the development of porcine-specific models for regenerative biology and its medical applications. In the past, this effort has been hampered by a lack of understanding of the signaling milieu that stabilizes the porcine pluripotent state in vitro. Here, we report that transgene-free PiPSCs can be efficiently derived from porcine fibroblasts by episomal vectors along with microRNA-302/367 using optimized protocols tailored for this species. PiPSCs can be differentiated into derivatives representing the primary germ layers in vitro and can form teratomas in immunocompromised mice. Furthermore, the transgene-free PiPSCs preserve intrinsic species-specific developmental timing in culture, known as developmental allochrony. This is demonstrated by establishing a porcine in vitro segmentation clock model that, for the first time, displays a specific periodicity at ∼3.7 h, a timescale recapitulating in vivo porcine somitogenesis. We conclude that the transgene-free PiPSCs can serve as a powerful tool for modeling development and disease and developing transplantation strategies. We also anticipate that they will provide insights into conserved and unique features on the regulations of mammalian pluripotency and developmental timing mechanisms.