Single-cell transcriptomes of the regenerating intestine reveal a revival stem cell.

The turnover of the intestinal epithelium is driven by multipotent LGR5+ crypt-base columnar cells (CBCs) located at the bottom of crypt zones1. However, CBCs are lost following injury, such as irradiation2, but the intestinal epithelium is nevertheless able to recover3. Thus, a second population of quiescent '+4' cells, or reserve stem cells (RSCs), has previously been proposed to regenerate the damaged intestine4-7. Although CBCs and RSCs were thought to be mutually exclusive4,8, subsequent studies have found that LGR5+ CBCs express RSC markers9 and that RSCs were dispensable-whereas LGR5+ cells were essential-for repair of the damaged intestine3. In addition, progenitors of absorptive enterocytes10, secretory cells11-15 and slow cycling LGR5+ cells16 have been shown to contribute to regeneration whereas the transcriptional regulator YAP1, which is important for intestinal regeneration, was suggested to induce a pro-survival phenotype in LGR5+ cells17. Thus, whether cellular plasticity or distinct cell populations are critical for intestinal regeneration remains unknown. Here we applied single-cell RNA sequencing to profile the regenerating mouse intestine and identified a distinct, damage-induced quiescent cell type that we term the revival stem cell (revSC). revSCs are marked by high clusterin expression and are extremely rare under homoeostatic conditions, yet give rise-in a temporal hierarchy-to all the major cell types of the intestine, including LGR5+ CBCs. After intestinal damage by irradiation, targeted ablation of LGR5+ CBCs, or treatment with dextran sodium sulfate, revSCs undergo a YAP1-dependent transient expansion, reconstitute the LGR5+ CBC compartment and are required to regenerate a functional intestine. These studies thus define a unique stem cell that is mobilized by damage to revive the homoeostatic stem cell compartment and regenerate the intestinal epithelium.

Genome-Wide Association Analysis of Young-Onset Stroke Identifies a Locus on Chromosome 10q25 Near HABP2.

BACKGROUND AND PURPOSE: Although a genetic contribution to ischemic stroke is well recognized, only a handful of stroke loci have been identified by large-scale genetic association studies to date. Hypothesizing that genetic effects might be stronger for early- versus late-onset stroke, we conducted a 2-stage meta-analysis of genome-wide association studies, focusing on stroke cases with an age of onset <60 years. METHODS: The discovery stage of our genome-wide association studies included 4505 cases and 21 968 controls of European, South-Asian, and African ancestry, drawn from 6 studies. In Stage 2, we selected the lead genetic variants at loci with association P<5×10(-6) and performed in silico association analyses in an independent sample of ≤1003 cases and 7745 controls. RESULTS: One stroke susceptibility locus at 10q25 reached genome-wide significance in the combined analysis of all samples from the discovery and follow-up stages (rs11196288; odds ratio =1.41; P=9.5×10(-9)). The associated locus is in an intergenic region between TCF7L2 and HABP2. In a further analysis in an independent sample, we found that 2 single nucleotide polymorphisms in high linkage disequilibrium with rs11196288 were significantly associated with total plasma factor VII-activating protease levels, a product of HABP2. CONCLUSIONS: HABP2, which encodes an extracellular serine protease involved in coagulation, fibrinolysis, and inflammatory pathways, may be a genetic susceptibility locus for early-onset stroke.

A high-performance computing toolset for relatedness and principal component analysis of SNP data.

Genome-wide association studies are widely used to investigate the genetic basis of diseases and traits, but they pose many computational challenges. We developed gdsfmt and SNPRelate (R packages for multi-core symmetric multiprocessing computer architectures) to accelerate two key computations on SNP data: principal component analysis (PCA) and relatedness analysis using identity-by-descent measures. The kernels of our algorithms are written in C/C++ and highly optimized. Benchmarks show the uniprocessor implementations of PCA and identity-by-descent are ∼8-50 times faster than the implementations provided in the popular EIGENSTRAT (v3.0) and PLINK (v1.07) programs, respectively, and can be sped up to 30-300-fold by using eight cores. SNPRelate can analyse tens of thousands of samples with millions of SNPs. For example, our package was used to perform PCA on 55 324 subjects from the 'Gene-Environment Association Studies' consortium studies.

GWASTools: an R/Bioconductor package for quality control and analysis of genome-wide association studies.

GWASTools is an R/Bioconductor package for quality control and analysis of genome-wide association studies (GWAS). GWASTools brings the interactive capability and extensive statistical libraries of R to GWAS. Data are stored in NetCDF format to accommodate extremely large datasets that cannot fit within R's memory limits. The documentation includes instructions for converting data from multiple formats, including variants called from sequencing. GWASTools provides a convenient interface for linking genotypes and intensity data with sample and single nucleotide polymorphism annotation.

Identification of Two Genetic Loci Associated with Leukopenia after Chemotherapy in Patients with Breast Cancer.

PURPOSE: To identify molecular predictors of grade 3/4 neutropenic or leukopenic events (NLE) after chemotherapy using a genome-wide association study (GWAS). EXPERIMENTAL DESIGN: A GWAS was performed on patients in the phase III chemotherapy study SUCCESS-A (n = 3,322). Genotyping was done using the Illumina HumanOmniExpress-12v1 array. Findings were functionally validated with cell culture models and the genotypes and gene expression of possible causative genes were correlated with clinical treatment response and prognostic outcomes. RESULTS: One locus on chromosome 16 (rs4784750; NLRC5; P = 1.56E-8) and another locus on chromosome 13 (rs16972207; TNFSF13B; P = 3.42E-8) were identified at a genome-wide significance level. Functional validation revealed that expression of these two genes is altered by genotype-dependent and chemotherapy-dependent activity of two transcription factors. Genotypes also showed an association with disease-free survival in patients with an NLE. CONCLUSIONS: Two loci in NLRC5 and TNFSF13B are associated with NLEs. The involvement of the MHC I regulator NLRC5 implies the possible involvement of immuno-oncological pathways.

A multiplexed, next generation sequencing platform for high-throughput detection of SARS-CoV-2.

Population scale sweeps of viral pathogens, such as SARS-CoV-2, require high intensity testing for effective management. Here, we describe "Systematic Parallel Analysis of RNA coupled to Sequencing for Covid-19 screening" (C19-SPAR-Seq), a multiplexed, scalable, readily automated platform for SARS-CoV-2 detection that is capable of analyzing tens of thousands of patient samples in a single run. To address strict requirements for control of assay parameters and output demanded by clinical diagnostics, we employ a control-based Precision-Recall and Receiver Operator Characteristics (coPR) analysis to assign run-specific quality control metrics. C19-SPAR-Seq coupled to coPR on a trial cohort of several hundred patients performs with a specificity of 100% and sensitivity of 91% on samples with low viral loads, and a sensitivity of >95% on high viral loads associated with disease onset and peak transmissibility. This study establishes the feasibility of employing C19-SPAR-Seq for the large-scale monitoring of SARS-CoV-2 and other pathogens.

Next-generation RNA sequencing of archival formalin-fixed paraffin-embedded urothelial bladder cancer.

UNLABELLED: Molecular profiling of individual cancers is key to personalised medicine. While sequencing technologies have required stringent sample collection and handling, recent technical advances offer sequencing from tissues collected in routine practice and tissues already stored in archives. In this paper, we establish methods for whole-transcriptome RNA sequencing (RNA-seq) from formalin-fixed paraffin-embedded tissues. We obtain average RNA-seq reads of >100 million per sample using the Illumina HiSeq2000 platform. We find high concordance with results from matching fresh frozen samples (>0.8 Spearman correlation). For validation, we compared low- and high-grade bladder cancer transcriptomes in 49 tumour samples after transurethral resection of bladder tumour. We found 947 differentially expressed protein-coding genes. While high-grade lesions exhibited distinct intertumour transcriptome heterogeneity, the transcriptome of low-grade tumours was homogeneous. PATIENT SUMMARY: In this report, we show that it is now possible to use universally available bladder cancer samples that have been fixed in formalin to perform high-quality transcriptome analysis. This ability will facilitate the development of transcriptome-wide tests based on gene expression correlated with clinical outcome.

Evaluating genetic risk for prostate cancer among Japanese and Latinos.

BACKGROUND: There have been few genome-wide association studies (GWAS) of prostate cancer among diverse populations. To search for novel prostate cancer risk variants, we conducted GWAS of prostate cancer in Japanese and Latinos. In addition, we tested prostate cancer risk variants and developed genetic risk models of prostate cancer for Japanese and Latinos. METHODS: Our first-stage GWAS of prostate cancer included Japanese (cases/controls = 1,033/1,042) and Latino (cases/controls = 1,043/1,057) from the Multiethnic Cohort (MEC). Significant associations from stage I (P < 1.0 × 10(-4)) were examined in silico in GWAS of prostate cancer (stage II) in Japanese (cases/controls = 1,583/3,386) and Europeans (cases/controls = 1,854/1,894). RESULTS: No novel stage I single-nucleotide polymorphism (SNP) outside of known risk regions reached genome-wide significance. For Japanese, in stage I, the most notable putative novel association was seen with 10 SNPs (P ≤ 8.0 × 10(-6)) at chromosome 2q33; however, this was not replicated in stage II. For Latinos, the most significant association was observed with rs17023900 at the known 3p12 risk locus (stage I: OR = 1.45; P = 7.01 × 10(-5) and stage II: OR = 1.58; P = 3.05 × 10(-7)). The majority of the established risk variants for prostate cancer, 79% and 88%, were positively associated with prostate cancer in Japanese and Latinos (stage I), respectively. The cumulative effects of these variants significantly influence prostate cancer risk (OR per allele = 1.10; P = 2.71 × 10(-25) and OR = 1.07; P = 1.02 × 10(-16) for Japanese and Latinos, respectively). CONCLUSION AND IMPACT: Our GWAS of prostate cancer did not identify novel genome-wide significant variants. However, our findings show that established risk variants for prostate cancer significantly contribute to risk among Japanese and Latinos.

Genome-wide association analysis of ischemic stroke in young adults.

Ischemic stroke (IS) is among the leading causes of death in Western countries. There is a significant genetic component to IS susceptibility, especially among young adults. To date, research to identify genetic loci predisposing to stroke has met only with limited success. We performed a genome-wide association (GWA) analysis of early-onset IS to identify potential stroke susceptibility loci. The GWA analysis was conducted by genotyping 1 million SNPs in a biracial population of 889 IS cases and 927 controls, ages 15-49 years. Genotypes were imputed using the HapMap3 reference panel to provide 1.4 million SNPs for analysis. Logistic regression models adjusting for age, recruitment stages, and population structure were used to determine the association of IS with individual SNPs. Although no single SNP reached genome-wide significance (P < 5 × 10(-8)), we identified two SNPs in chromosome 2q23.3, rs2304556 (in FMNL2; P = 1.2 × 10(-7)) and rs1986743 (in ARL6IP6; P = 2.7 × 10(-7)), strongly associated with early-onset stroke. These data suggest that a novel locus on human chromosome 2q23.3 may be associated with IS susceptibility among young adults.

Using cluster ensemble and validation to identify subtypes of pervasive developmental disorders.

Pervasive Developmental Disorders (PDD) are neurodevelopmental disorders characterized by impairments in social interaction, communication and behavior. Given the diversity and varying severity of PDD, diagnostic tools attempt to identify homogeneous subtypes within PDD. Identifying subtypes can lead to targeted etiology studies and to effective type-specific intervention. Cluster analysis can suggest coherent subsets in data; however, different methods and assumptions lead to different results. Several previous studies applied clustering to PDD data, varying in number and characteristics of the produced subtypes. Most studies used a relatively small dataset (fewer than 150 subjects), and all applied only a single clustering method. Here we study a relatively large dataset (358 PDD patients), using an ensemble of three clustering methods. The results are evaluated using several validation methods, and consolidated through an integration step. Four clusters are identified, analyzed and compared to subtypes previously defined by the widely used diagnostic tool DSM-IV.