Dissecting quantitative trait nucleotides by saturation genome editing.

Genome editing technologies have the potential to transform our understanding of how genetic variation gives rise to complex traits through the systematic engineering and phenotypic characterization of genetic variants. However, there has yet to be a system with sufficient efficiency, fidelity, and throughput to comprehensively identify causal variants at the genome scale. Here we explored the ability of templated CRISPR editing systems to install natural variants genome-wide in budding yeast. We optimized several approaches to enhance homology-directed repair (HDR) with donor DNA templates, including donor recruitment to target sites, single-stranded donor production by bacterial retrons, and in vivo plasmid assembly. We uncovered unique advantages of each system that we integrated into a single superior system named MAGESTIC 3.0. We used MAGESTIC 3.0 to dissect causal variants residing in 112 quantitative trait loci across 32 environmental conditions, revealing an enrichment for missense variants and loci with multiple causal variants. MAGESTIC 3.0 will facilitate the functional analysis of the genome at single-nucleotide resolution and provides a roadmap for improving template-based genome editing systems in other organisms.

Striated muscle-specific base editing enables correction of mutations causing dilated cardiomyopathy.

Dilated cardiomyopathy is the second most common cause for heart failure with no cure except a high-risk heart transplantation. Approximately 30% of patients harbor heritable mutations which are amenable to CRISPR-based gene therapy. However, challenges related to delivery of the editing complex and off-target concerns hamper the broad applicability of CRISPR agents in the heart. We employ a combination of the viral vector AAVMYO with superior targeting specificity of heart muscle tissue and CRISPR base editors to repair patient mutations in the cardiac splice factor Rbm20, which cause aggressive dilated cardiomyopathy. Using optimized conditions, we repair >70% of cardiomyocytes in two Rbm20 knock-in mouse models that we have generated to serve as an in vivo platform of our editing strategy. Treatment of juvenile mice restores the localization defect of RBM20 in 75% of cells and splicing of RBM20 targets including TTN. Three months after injection, cardiac dilation and ejection fraction reach wild-type levels. Single-nuclei RNA sequencing uncovers restoration of the transcriptional profile across all major cardiac cell types and whole-genome sequencing reveals no evidence for aberrant off-target editing. Our study highlights the potential of base editors combined with AAVMYO to achieve gene repair for treatment of hereditary cardiac diseases.

A scalable, GMP-compatible, autologous organotypic cell therapy for Dystrophic Epidermolysis Bullosa.

Background: Gene editing in induced pluripotent stem (iPS) cells has been hailed to enable new cell therapies for various monogenetic diseases including dystrophic epidermolysis bullosa (DEB). However, manufacturing, efficacy and safety roadblocks have limited the development of genetically corrected, autologous iPS cell-based therapies. Methods: We developed Dystrophic Epidermolysis Bullosa Cell Therapy (DEBCT), a new generation GMP-compatible (cGMP), reproducible, and scalable platform to produce autologous clinical-grade iPS cell-derived organotypic induced skin composite (iSC) grafts to treat incurable wounds of patients lacking type VII collagen (C7). DEBCT uses a combined high-efficiency reprogramming and CRISPR-based genetic correction single step to generate genome scar-free, COL7A1 corrected clonal iPS cells from primary patient fibroblasts. Validated iPS cells are converted into epidermal, dermal and melanocyte progenitors with a novel 2D organoid differentiation protocol, followed by CD49f enrichment and expansion to minimize maturation heterogeneity. iSC product characterization by single cell transcriptomics was followed by mouse xenografting for disease correcting activity at 1 month and toxicology analysis at 1-6 months. Culture-acquired mutations, potential CRISPR-off targets, and cancer-driver variants were evaluated by targeted and whole genome sequencing. Findings: iPS cell-derived iSC grafts were reproducibly generated from four recessive DEB patients with different pathogenic mutations. Organotypic iSC grafts onto immune-compromised mice developed into stable stratified skin with functional C7 restoration. Single cell transcriptomic characterization of iSCs revealed prominent holoclone stem cell signatures in keratinocytes and the recently described Gibbin-dependent signature in dermal fibroblasts. The latter correlated with enhanced graftability. Multiple orthogonal sequencing and subsequent computational approaches identified random and non-oncogenic mutations introduced by the manufacturing process. Toxicology revealed no detectable tumors after 3-6 months in DEBCT-treated mice. Interpretation: DEBCT successfully overcomes previous roadblocks and represents a robust, scalable, and safe cGMP manufacturing platform for production of a CRISPR-corrected autologous organotypic skin graft to heal DEB patient wounds.

Fast and inexpensive whole-genome sequencing library preparation from intact yeast cells.

Through the increase in the capacity of sequencing machines massively parallel sequencing of thousands of samples in a single run is now possible. With the improved throughput and resulting drop in the price of sequencing, the cost and time for preparation of sequencing libraries have become the major bottleneck in large-scale experiments. Methods using a hyperactive variant of the Tn5 transposase efficiently generate libraries starting from cDNA or genomic DNA in a few hours and are highly scalable. For genome sequencing, however, the time and effort spent on genomic DNA isolation limit the practicability of sequencing large numbers of samples. Here, we describe a highly scalable method for preparing high-quality whole-genome sequencing libraries directly from Saccharomyces cerevisiae cultures in less than 3 h at 34 cents per sample. We skip the rate-limiting step of genomic DNA extraction by directly tagmenting lysed yeast spheroplasts and add a nucleosome release step prior to enrichment PCR to improve the evenness of genomic coverage. Resulting libraries do not show any GC bias and are comparable in quality to libraries processed from genomic DNA with a commercially available Tn5-based kit. We use our protocol to investigate CRISPR/Cas9 on- and off-target edits and reliably detect edited variants and shared polymorphisms between strains. Our protocol enables rapid preparation of unbiased and high-quality, sequencing-ready indexed libraries for hundreds of yeast strains in a single day at a low price. By adjusting individual steps of our workflow, we expect that our protocol can be adapted to other organisms.

Novel gene regulatory networks identified in response to nitro-conjugated linoleic acid in human endothelial cells.

Endothelial cell (EC) dysfunction is a crucial initiation event in the development of atherosclerosis and is associated with diabetes mellitus, hypertension, and heart failure. Both digestive and oxidative inflammatory conditions lead to the endogenous formation of nitrated derivatives of unsaturated fatty acids (FAs) upon generation of the proximal nitrating species nitrogen dioxide (·NO2) by nitric oxide (·NO) and nitrite-dependent reactions. Nitro-FAs (NO2-FAs) such as nitro-oleic acid (NO2-OA) and nitro-linoleic acid (NO2-LA) potently inhibit inflammation and oxidative stress, regulate cellular functions, and maintain cardiovascular homeostasis. Recently, conjugated linoleic acid (CLA) was identified as the preferential FA substrate of nitration in vivo. However, the functions of nitro-CLA (NO2-CLA) in ECs remain to be explored. In the present study, a distinct transcriptome regulated by NO2-CLA was revealed in primary human coronary artery endothelial cells (HCAECs) through RNA sequencing. Differential gene expression and pathway enrichment analysis identified numerous regulatory networks including those related to the modulation of inflammation, oxidative stress, cell cycle, and hypoxic responses by NO2-CLA, suggesting a diverse impact of NO2-CLA and other electrophilic nitrated FAs on cellular processes. These findings extend the understanding of the protective actions of NO2-CLA in cardiovascular diseases and provide new insight into the underlying mechanisms that mediate the pleiotropic cellular responses to NO2-CLA.

Rapid evolution of a retro-transposable hotspot of ovine genome underlies the alteration of BMP2 expression and development of fat tails.

BACKGROUND: Sheep have developed the ability to store fat in their tails, which is a unique way of reserving energy to survive a harsh environment. However, the mechanism underlying this adaptive trait remains largely unsolved. RESULTS: In the present study, we provide evidence for the genetic determinants of fat tails, based on whole genome sequences of 89 individual sheep. A genome-wide scan of selective sweep identified several candidate loci including a region at chromosome 13, a haplotype of which underwent rapid evolution and spread through fat-tailed populations in China and the Middle East. Sequence analysis revealed an inter-genic origin of this locus, which later became a hotspot of ruminant-specific retro-transposon named BovB. Additionally, the candidate locus was validated based on a fat- and thin-tailed cross population. The expression of an upstream gene BMP2 was differentially regulated between fat-tailed and thin-tailed individuals in tail adipose and several other tissue types. CONCLUSIONS: Our findings suggest the fixation of fat tails in domestic sheep is caused by a selective sweep near a retro-transposable hotspot at chromosome 13, the diversity of which specifically affects the expression of BMP2. The present study has shed light onto the understanding of fat metabolism.

Whole-genome sequences of 89 Chinese sheep suggest role of RXFP2 in the development of unique horn phenotype as response to semi-feralization.

Background: Animal domestication has been extensively studied, but the process of feralization remains poorly understood. Results: Here, we performed whole-genome sequencing of 99 sheep and identified a primary genetic divergence between 2 heterogeneous populations in the Tibetan Plateau, including 1 semi-feral lineage. Selective sweep and candidate gene analysis revealed local adaptations of these sheep associated with sensory perception, muscle strength, eating habit, mating process, and aggressive behavior. In particular, a horn-related gene, RXFP2, showed signs of rapid evolution specifically in the semi-feral breeds. A unique haplotype and repressed horn-related tissue expression of RXFP2 were correlated with higher horn length, as well as spiral and horizontally extended horn shape. Conclusions: Semi-feralization has an extensive impact on diverse phenotypic traits of sheep. By acquiring features like those of their wild ancestors, semi-feral sheep were able to regain fitness while in frequent contact with wild surroundings and rare human interventions. This study provides a new insight into the evolution of domestic animals when human interventions are no longer dominant.

Transcriptomic sequencing reveals diverse adaptive gene expression responses of human vascular smooth muscle cells to nitro-conjugated linoleic acid.

Nitro-conjugated linoleic acid (NO2-CLA) is formed by metabolic and inflammatory reactions of nitric oxide and nitrite, and represents the most abundant nitro-fatty acid species in humans. These electrophilic fatty acid nitroalkene derivatives mediate pleiotropic cell signaling responses. Here, we report a systematic approach to investigate the effect of NO2-CLA on human coronary artery smooth muscle cells (hCASMC), based on the RNA-Seq and bioinformatics analysis. There were extensive differentially expressed genes in NO2-CLA vs. control (510) and NO2-CLA vs. CLA (272) treatment groups, respectively. Notably, only minimal alterations were observed in CLA vs. control conditions, indicating that the electrophilic character of NO2-CLA is requited to induce differential gene expression responses independently from native CLA. Functional enrichment analysis of differentially expressed genes reveals multiple cellular processes to be affected under NO2-CLA treatment, including cell proliferation, lipid metabolism, antioxidant and inflammatory-related gene expression responses. These findings reveal that nitro-fatty acid derivatives such as NO2-CLA regulate a broad array of adaptive gene expression responses by hCASMC.

Laminar Flow Attenuates Macrophage Migration Inhibitory Factor Expression in Endothelial Cells.

Macrophage migration inhibitory factor (MIF) is a non-canonical cytokine that is involved in multiple inflammatory diseases, including atherosclerosis. High MIF expression found in leukocytes which facilitates the initiation and progression of atherosclerosis. However, little is known about biomechanical forces in the induction of MIF in endothelial cells (ECs). Here, we show that laminar shear stress (LS) inhibits the expression of MIF in ECs. By profiling the whole transcriptome of human coronary artery ECs under different shear stress, we found that athero-protective LS attenuates the expression of MIF whereas pro-atherosclerotic oscillatory shear stress (OS) significantly increased the expression of MIF. En face staining of rabbit aorta revealed high MIF immunoreactivity in lesser curvature as well as arterial bifurcation areas where OS is predominant. Mechanistically, we found that Krüpple like factor 2 (KLF2) is required for inhibition of MIF expression in ECs in the context of shear stress. Knockdown of KLF2 abolishes LS-dependent MIF inhibition while overexpression of KLF2 significantly attenuated MIF expression. Overall, the present work showed that MIF is a shear stress-sensitive cytokine and is transcriptionally regulated by KLF2, suggesting that LS exerts its athero-protective effect in part by directly inhibiting pro-inflammatory MIF expression.

Adaptive evolution of proteins in hepatitis B virus during divergence of genotypes.

Hepatitis B virus (HBV) is classified into several genotypes, correlated with different geographic distributions, clinical outcomes and susceptible human populations. It is crucial to investigate the evolutionary significance behind the diversification of HBV genotypes, because it improves our understanding of their pathological differences and pathogen-host interactions. Here, we performed comprehensive analysis of HBV genome sequences collected from public database. With a stringent criteria, we generated a dataset of 2992 HBV genomes from eight major genotypes. In particular, we applied a specified classification of non-synonymous and synonymous variants in overlapping regions, to distinguish joint and independent gene evolutions. We confirmed the presence of selective constraints over non-synonymous variants in consideration of overlapping regions. We then performed the McDonald-Kreitman test and revealed adaptive evolutions of non-synonymous variants during genotypic differentiation. Remarkably, we identified strong positive selection that drove the differentiation of PreS1 domain, which is an essential regulator involved in viral transmission. Our study presents novel evidences for the adaptive evolution of HBV genotypes, which suggests that these viruses evolve directionally for maintenance or improvement of successful infections.

Chemical biology reveals CARF as a positive regulator of canonical Wnt signaling by promoting TCF/ß-catenin transcriptional activity.

Wnt/ß-catenin signaling regulates multiple biological processes and aberration of this pathway is frequently observed in human cancers. Previously, we uncovered NC043 as a small-molecule inhibitor of Wnt/ß-catenin signaling. Here, we identified CARF as the cellular target of NC043. We found that NC043 binds directly to CARF through forming a covalent bond with the Cys-516 residue of CARF. Further study revealed that CARF interacts with Dvl, which potentiates the Dvl-c-Jun-ß-catenin-TCF transcriptional complex and thus promotes Wnt signaling activation. NC043 could disrupt the interaction between CARF and Dvl, thereby impairing Wnt signal transduction. In line with this, knockdown of CARF in zebrafish leads to impairment of embryonic development, hematopoietic stem cell generation and caudal fin regeneration. Collectively, we identified CARF as the cellular target of NC043 and revealed CARF as a positive regulator of Wnt/ß-catenin signal transduction.

Hyperlipidemia-associated gene variations and expression patterns revealed by whole-genome and transcriptome sequencing of rabbit models.

The rabbit (Oryctolagus cuniculus) is an important experimental animal for studying human diseases, such as hypercholesterolemia and atherosclerosis. Despite this, genetic information and RNA expression profiling of laboratory rabbits are lacking. Here, we characterized the whole-genome variants of three breeds of the most popular experimental rabbits, New Zealand White (NZW), Japanese White (JW) and Watanabe heritable hyperlipidemic (WHHL) rabbits. Although the genetic diversity of WHHL rabbits was relatively low, they accumulated a large proportion of high-frequency deleterious mutations due to the small population size. Some of the deleterious mutations were associated with the pathophysiology of WHHL rabbits in addition to the LDLR deficiency. Furthermore, we conducted transcriptome sequencing of different organs of both WHHL and cholesterol-rich diet (Chol)-fed NZW rabbits. We found that gene expression profiles of the two rabbit models were essentially similar in the aorta, even though they exhibited different types of hypercholesterolemia. In contrast, Chol-fed rabbits, but not WHHL rabbits, exhibited pronounced inflammatory responses and abnormal lipid metabolism in the liver. These results provide valuable insights into identifying therapeutic targets of hypercholesterolemia and atherosclerosis with rabbit models.

Genome sequence of Anopheles sinensis provides insight into genetics basis of mosquito competence for malaria parasites.

BACKGROUND: Anopheles sinensis is an important mosquito vector of Plasmodium vivax, which is the most frequent and widely distributed cause of recurring malaria throughout Asia, and particularly in China, Korea, and Japan. RESULTS: We performed 454 next-generation sequencing and obtained a draft sequence of A. sinensis assembled into scaffolds spanning 220.8 million base pairs. Analysis of this genome sequence, we observed expansion and contraction of several immune-related gene families in anopheline relative to culicine mosquito species. These differences suggest that species-specific immune responses to Plasmodium invasion underpin the biological differences in susceptibility to Plasmodium infection that characterize these two mosquito subfamilies. CONCLUSIONS: The A. sinensis genome produced in this study, provides an important resource for analyzing the genetic basis of susceptibility and resistance of mosquitoes to Plasmodium parasites research which will ultimately facilitate the design of urgently needed interventions against this debilitating mosquito-borne disease.

The de novo sequence origin of two long non-coding genes from an inter-genic region.

BACKGROUND: The gene Polymorphic derived intron-containing, known as Pldi, is a long non-coding RNA (lncRNA) first discovered in mouse. Although parts of its sequence were reported to be conserved in rat and human, it can only be expressed in mouse testis with a mouse-specific transcription start site. The consensus sequence of Pldi is also part of an antisense transcript AK158810 expressed in a wide range of mouse tissues. RESULT: We focused on sequence origin of Pldi and Ak158810. We demonstrated that their sequence was originated from an inter-genic region and is only presented in mammalians. Transposable events and chromosome rearrangements were involved in the evolution of ancestral sequence. Moreover, we discovered high conservation in part of this region was correlated with chromosome rearrangements, CpG demethylation and transcriptional factor binding motif. These results demonstrated that multiple factors contributed to the sequence origin of Pldi. CONCLUSIONS: We comprehensively analyzed the sequence origin of Pldi-Ak158810 loci. We provided various factors, including rearrangement, transposable elements, contributed to the formation of the sequence.

Genome sequences of wild and domestic bactrian camels.

Bactrian camels serve as an important means of transportation in the cold desert regions of China and Mongolia. Here we present a 2.01 Gb draft genome sequence from both a wild and a domestic bactrian camel. We estimate the camel genome to be 2.38 Gb, containing 20,821 protein-coding genes. Our phylogenomics analysis reveals that camels shared common ancestors with other even-toed ungulates about 55-60 million years ago. Rapidly evolving genes in the camel lineage are significantly enriched in metabolic pathways, and these changes may underlie the insulin resistance typically observed in these animals. We estimate the genome-wide heterozygosity rates in both wild and domestic camels to be 1.0 × 10(-3). However, genomic regions with significantly lower heterozygosity are found in the domestic camel, and olfactory receptors are enriched in these regions. Our comparative genomics analyses may also shed light on the genetic basis of the camel's remarkable salt tolerance and unusual immune system.