Transmission from vaccinated individuals in a large SARS-CoV-2 Delta variant outbreak.

An outbreak of over 1,000 COVID-19 cases in Provincetown, Massachusetts (MA), in July 2021-the first large outbreak mostly in vaccinated individuals in the US-prompted a comprehensive public health response, motivating changes to national masking recommendations and raising questions about infection and transmission among vaccinated individuals. To address these questions, we combined viral genomic and epidemiological data from 467 individuals, including 40% of outbreak-associated cases. The Delta variant accounted for 99% of cases in this dataset; it was introduced from at least 40 sources, but 83% of cases derived from a single source, likely through transmission across multiple settings over a short time rather than a single event. Genomic and epidemiological data supported multiple transmissions of Delta from and between fully vaccinated individuals. However, despite its magnitude, the outbreak had limited onward impact in MA and the US overall, likely due to high vaccination rates and a robust public health response.

Recurrent and functional regulatory mutations in breast cancer.

Genomic analysis of tumours has led to the identification of hundreds of cancer genes on the basis of the presence of mutations in protein-coding regions. By contrast, much less is known about cancer-causing mutations in non-coding regions. Here we perform deep sequencing in 360 primary breast cancers and develop computational methods to identify significantly mutated promoters. Clear signals are found in the promoters of three genes. FOXA1, a known driver of hormone-receptor positive breast cancer, harbours a mutational hotspot in its promoter leading to overexpression through increased E2F binding. RMRP and NEAT1, two non-coding RNA genes, carry mutations that affect protein binding to their promoters and alter expression levels. Our study shows that promoter regions harbour recurrent mutations in cancer with functional consequences and that the mutations occur at similar frequencies as in coding regions. Power analyses indicate that more such regions remain to be discovered through deep sequencing of adequately sized cohorts of patients.

Detection of low-density Plasmodium falciparum infections using amplicon deep sequencing.

BACKGROUND: Deep sequencing of targeted genomic regions is becoming a common tool for understanding the dynamics and complexity of Plasmodium infections, but its lower limit of detection is currently unknown. Here, a new amplicon analysis tool, the Parallel Amplicon Sequencing Error Correction (PASEC) pipeline, is used to evaluate the performance of amplicon sequencing on low-density Plasmodium DNA samples. Illumina-based sequencing of two Plasmodium falciparum genomic regions (CSP and SERA2) was performed on two types of samples: in vitro DNA mixtures mimicking low-density infections (1-200 genomes/µl) and extracted blood spots from a combination of symptomatic and asymptomatic individuals (44-653,080 parasites/µl). Three additional analysis tools-DADA2, HaplotypR, and SeekDeep-were applied to both datasets and the precision and sensitivity of each tool were evaluated. RESULTS: Amplicon sequencing can contend with low-density samples, showing reasonable detection accuracy down to a concentration of 5 Plasmodium genomes/µl. Due to increased stochasticity and background noise, however, all four tools showed reduced sensitivity and precision on samples with very low parasitaemia (< 5 copies/µl) or low read count (< 100 reads per amplicon). PASEC could distinguish major from minor haplotypes with an accuracy of 90% in samples with at least 30 Plasmodium genomes/µl, but only 61% at low Plasmodium concentrations (< 5 genomes/µl) and 46% at very low read counts (< 25 reads per amplicon). The four tools were additionally used on a panel of extracted parasite-positive blood spots from natural malaria infections. While all four identified concordant patterns of complexity of infection (COI) across four sub-Saharan African countries, COI values obtained for individual samples differed in some cases. CONCLUSIONS: Amplicon deep sequencing can be used to determine the complexity and diversity of low-density Plasmodium infections. Despite differences in their approach, four state-of-the-art tools resolved known haplotype mixtures with similar sensitivity and precision. Researchers can therefore choose from multiple robust approaches for analysing amplicon data, however, error filtration approaches should not be uniformly applied across samples of varying parasitaemia. Samples with very low parasitaemia and very low read count have higher false positive rates and call for read count thresholds that are higher than current default recommendations.

Whole-genome and multisector exome sequencing of primary and post-treatment glioblastoma reveals patterns of tumor evolution.

Glioblastoma (GBM) is a prototypical heterogeneous brain tumor refractory to conventional therapy. A small residual population of cells escapes surgery and chemoradiation, resulting in a typically fatal tumor recurrence ∼ 7 mo after diagnosis. Understanding the molecular architecture of this residual population is critical for the development of successful therapies. We used whole-genome sequencing and whole-exome sequencing of multiple sectors from primary and paired recurrent GBM tumors to reconstruct the genomic profile of residual, therapy resistant tumor initiating cells. We found that genetic alteration of the p53 pathway is a primary molecular event predictive of a high number of subclonal mutations in glioblastoma. The genomic road leading to recurrence is highly idiosyncratic but can be broadly classified into linear recurrences that share extensive genetic similarity with the primary tumor and can be directly traced to one of its specific sectors, and divergent recurrences that share few genetic alterations with the primary tumor and originate from cells that branched off early during tumorigenesis. Our study provides mechanistic insights into how genetic alterations in primary tumors impact the ensuing evolution of tumor cells and the emergence of subclonal heterogeneity.

Deformation-induced transitional myofibroblasts contribute to compensatory lung growth.

In many mammals, including humans, removal of one lung (pneumonectomy) results in the compensatory growth of the remaining lung. Compensatory growth involves not only an increase in lung size, but also an increase in the number of alveoli in the peripheral lung; however, the process of compensatory neoalveolarization remains poorly understood. Here, we show that the expression of α-smooth muscle actin (SMA)-a cytoplasmic protein characteristic of myofibroblasts-is induced in the pleura following pneumonectomy. SMA induction appears to be dependent on pleural deformation (stretch) as induction is prevented by plombage or phrenic nerve transection (P < 0.001). Within 3 days of pneumonectomy, the frequency of SMA+ cells in subpleural alveolar ducts was significantly increased (P < 0.01). To determine the functional activity of these SMA+ cells, we isolated regenerating alveolar ducts by laser microdissection and analyzed individual cells using microfluidic single-cell quantitative PCR. Single cells expressing the SMA (Acta2) gene demonstrated significantly greater transcriptional activity than endothelial cells or other discrete cell populations in the alveolar duct (P < 0.05). The transcriptional activity of the Acta2+ cells, including expression of TGF signaling as well as repair-related genes, suggests that these myofibroblast-like cells contribute to compensatory lung growth.

Novel loci for adiponectin levels and their influence on type 2 diabetes and metabolic traits: a multi-ethnic meta-analysis of 45,891 individuals.

Circulating levels of adiponectin, a hormone produced predominantly by adipocytes, are highly heritable and are inversely associated with type 2 diabetes mellitus (T2D) and other metabolic traits. We conducted a meta-analysis of genome-wide association studies in 39,883 individuals of European ancestry to identify genes associated with metabolic disease. We identified 8 novel loci associated with adiponectin levels and confirmed 2 previously reported loci (P = 4.5×10(-8)-1.2×10(-43)). Using a novel method to combine data across ethnicities (N = 4,232 African Americans, N = 1,776 Asians, and N = 29,347 Europeans), we identified two additional novel loci. Expression analyses of 436 human adipocyte samples revealed that mRNA levels of 18 genes at candidate regions were associated with adiponectin concentrations after accounting for multiple testing (p<3×10(-4)). We next developed a multi-SNP genotypic risk score to test the association of adiponectin decreasing risk alleles on metabolic traits and diseases using consortia-level meta-analytic data. This risk score was associated with increased risk of T2D (p = 4.3×10(-3), n = 22,044), increased triglycerides (p = 2.6×10(-14), n = 93,440), increased waist-to-hip ratio (p = 1.8×10(-5), n = 77,167), increased glucose two hours post oral glucose tolerance testing (p = 4.4×10(-3), n = 15,234), increased fasting insulin (p = 0.015, n = 48,238), but with lower in HDL-cholesterol concentrations (p = 4.5×10(-13), n = 96,748) and decreased BMI (p = 1.4×10(-4), n = 121,335). These findings identify novel genetic determinants of adiponectin levels, which, taken together, influence risk of T2D and markers of insulin resistance.

Blended Genome Exome (BGE) as a Cost Efficient Alternative to Deep Whole Genomes or Arrays.

Genomic scientists have long been promised cheaper DNA sequencing, but deep whole genomes are still costly, especially when considered for large cohorts in population-level studies. More affordable options include microarrays + imputation, whole exome sequencing (WES), or low-pass whole genome sequencing (WGS) + imputation. WES + array + imputation has recently been shown to yield 99% of association signals detected by WGS. However, a method free from ascertainment biases of arrays or the need for merging different data types that still benefits from deeper exome coverage to enhance novel coding variant detection does not exist. We developed a new, combined, "Blended Genome Exome" (BGE) in which a whole genome library is generated, an aliquot of that genome is amplified by PCR, the exome regions are selected and enriched, and the genome and exome libraries are combined back into a single tube for sequencing (33% exome, 67% genome). This creates a single CRAM with a low-coverage whole genome (2-3x) combined with a higher coverage exome (30-40x). This BGE can be used for imputing common variants throughout the genome as well as for calling rare coding variants. We tested this new method and observed >99% r 2 concordance between imputed BGE data and existing 30x WGS data for exome and genome variants. BGE can serve as a useful and cost-efficient alternative sequencing product for genomic researchers, requiring ten-fold less sequencing compared to 30x WGS without the need for complicated harmonization of array and sequencing data.

Genotypic analysis of RTS,S/AS01E malaria vaccine efficacy against parasite infection as a function of dosage regimen and baseline malaria infection status in children aged 5-17 months in Ghana and Kenya: a longitudinal phase 2b randomised controlled trial.

BACKGROUND: The first licensed malaria vaccine, RTS,S/AS01E, confers moderate protection against symptomatic disease. Because many malaria infections are asymptomatic, we conducted a large-scale longitudinal parasite genotyping study of samples from a clinical trial exploring how vaccine dosing regimen affects vaccine efficacy. METHODS: Between Sept 28, 2017, and Sept 25, 2018, 1500 children aged 5-17 months were randomly assigned (1:1:1:1:1) to receive four different RTS,S/AS01E regimens or a rabies control vaccine in a phase 2b open-label clinical trial in Ghana and Kenya. Participants in the four RTS,S groups received two full doses at month 0 and month 1 and either full doses at month 2 and month 20 (group R012-20); full doses at month 2, month 14, month 26, and month 38 (group R012-14); fractional doses at month 2, month 14, month 26, and month 38 (group Fx012-14; early fourth dose); or fractional doses at month 7, month 20, and month 32 (group Fx017-20; delayed third dose). We evaluated the time to the first new genotypically detected infection and the total number of new infections during two follow-up periods (12 months and 20 months) in more than 36 000 dried blood spot specimens from 1500 participants. To study vaccine effects on time to the first new infection, we defined vaccine efficacy as one minus the hazard ratio (HR; RTS,S vs control) of the first new infection. We performed a post-hoc analysis of vaccine efficacy based on malaria infection status at first vaccination and force of infection by month 2. This trial (MAL-095) is registered with ClinicalTrials.gov, NCT03281291. FINDINGS: We observed significant and similar vaccine efficacy (25-43%; 95% CI union 9-53) against first new infection for all four RTS,S/AS01E regimens across both follow-up periods (12 months and 20 months). Each RTS,S/AS01E regimen significantly reduced the mean number of new infections in the 20-month follow-up period by 1·1-1·6 infections (95% CI union 0·6-2·1). Vaccine efficacy against first new infection was significantly higher in participants who were infected with malaria (68%; 95% CI 50-80) than in those who were uninfected (37%; 23-48) at the first vaccination (p=0·0053). INTERPRETATION: All tested dosing regimens blocked some infections to a similar degree. Improved vaccine efficacy in participants infected during vaccination could suggest new strategies for highly efficacious malaria vaccine development and implementation. FUNDING: GlaxoSmithKline Biologicals SA, PATH, Bill & Melinda Gates Foundation, and the German Federal Ministry of Education and Research.

Baseline malaria infection status and RTS,S/AS01E malaria vaccine efficacy.

Background: The only licensed malaria vaccine, RTS,S/AS01 E , confers moderate protection against symptomatic disease. Because many malaria infections are asymptomatic, we conducted a large-scale longitudinal parasite genotyping study of samples from a clinical trial exploring how vaccine dosing regimen affects vaccine efficacy (VE). Methods: 1,500 children aged 5-17 months were randomized to receive four different RTS,S/AS01 E regimens or a rabies control vaccine in a phase 2b clinical trial in Ghana and Kenya. We evaluated the time to the first new genotypically detected infection and the total number of new infections during two follow-up periods in over 36K participant specimens. We performed a post hoc analysis of VE based on malaria infection status at first vaccination and force of infection. Results: We observed significant and comparable VE (25-43%, 95% CI union 9-53%) against first new infection for all four RTS,S/AS01 E regimens across both follow-up periods (12 and 20 months). Each RTS,S/AS01 E regimen significantly reduced the number of new infections in the 20-month follow-up period (control mean 4.1 vs. RTS,S/AS01 E mean 2.6-3.0). VE against first new infection was significantly higher in participants who were malaria-infected (68%; 95% CI, 50 to 80%) versus uninfected (37%; 95% CI, 23 to 48%) at the first vaccination (P=0.0053) and in participants experiencing greater force of infection between dose 1 and 3 (P=0.059). Conclusions: All tested dosing regimens blocked some infections to a similar degree. Improved VE in participants infected during vaccination could suggest new strategies for highly efficacious malaria vaccine development and implementation. ( ClinicalTrials.gov number, NCT03276962 ).

Race-ethnic differences in the association of genetic loci with HbA1c levels and mortality in U.S. adults: the third National Health and Nutrition Examination Survey (NHANES III).

BACKGROUND: Hemoglobin A1c (HbA1c) levels diagnose diabetes, predict mortality and are associated with ten single nucleotide polymorphisms (SNPs) in white individuals. Genetic associations in other race groups are not known. We tested the hypotheses that there is race-ethnic variation in 1) HbA1c-associated risk allele frequencies (RAFs) for SNPs near SPTA1, HFE, ANK1, HK1, ATP11A, FN3K, TMPRSS6, G6PC2, GCK, MTNR1B; 2) association of SNPs with HbA1c and 3) association of SNPs with mortality. METHODS: We studied 3,041 non-diabetic individuals in the NHANES (National Health and Nutrition Examination Survey) III. We stratified the analysis by race/ethnicity (NHW: non-Hispanic white; NHB: non-Hispanic black; MA: Mexican American) to calculate RAF, calculated a genotype score by adding risk SNPs, and tested associations with SNPs and the genotype score using an additive genetic model, with type 1 error = 0.05. RESULTS: RAFs varied widely and at six loci race-ethnic differences in RAF were significant (p < 0.0002), with NHB usually the most divergent. For instance, at ATP11A, the SNP RAF was 54% in NHB, 18% in MA and 14% in NHW (p < .0001). The mean genotype score differed by race-ethnicity (NHW: 10.4, NHB: 11.0, MA: 10.7, p < .0001), and was associated with increase in HbA1c in NHW (ß = 0.012 HbA1c increase per risk allele, p = 0.04) and MA (ß = 0.021, p = 0.005) but not NHB (ß = 0.007, p = 0.39). The genotype score was not associated with mortality in any group (NHW: OR (per risk allele increase in mortality) = 1.07, p = 0.09; NHB: OR = 1.04, p = 0.39; MA: OR = 1.03, p = 0.71). CONCLUSION: At many HbA1c loci in NHANES III there is substantial RAF race-ethnic heterogeneity. The combined impact of common HbA1c-associated variants on HbA1c levels varied by race-ethnicity, but did not influence mortality.

Evidence of transmission from fully vaccinated individuals in a large outbreak of the SARS-CoV-2 Delta variant in Provincetown, Massachusetts.

Multiple summer events, including large indoor gatherings, in Provincetown, Massachusetts (MA), in July 2021 contributed to an outbreak of over one thousand COVID-19 cases among residents and visitors. Most cases were fully vaccinated, many of whom were also symptomatic, prompting a comprehensive public health response, motivating changes to national masking recommendations, and raising questions about infection and transmission among vaccinated individuals. To characterize the outbreak and the viral population underlying it, we combined genomic and epidemiological data from 467 individuals, including 40% of known outbreak-associated cases. The Delta variant accounted for 99% of sequenced outbreak-associated cases. Phylogenetic analysis suggests over 40 sources of Delta in the dataset, with one responsible for a single cluster containing 83% of outbreak-associated genomes. This cluster was likely not the result of extensive spread at a single site, but rather transmission from a common source across multiple settings over a short time. Genomic and epidemiological data combined provide strong support for 25 transmission events from, including many between, fully vaccinated individuals; genomic data alone provides evidence for an additional 64. Together, genomic epidemiology provides a high-resolution picture of the Provincetown outbreak, revealing multiple cases of transmission of Delta from fully vaccinated individuals. However, despite its magnitude, the outbreak was restricted in its onward impact in MA and the US, likely due to high vaccination rates and a robust public health response.

Single-Cell Transcriptional Profiling of Cells Derived From Regenerating Alveolar Ducts.

Lung regeneration occurs in a variety of adult mammals after surgical removal of one lung (pneumonectomy). Previous studies of murine post-pneumonectomy lung growth have identified regenerative "hotspots" in subpleural alveolar ducts; however, the cell-types participating in this process remain unclear. To identify the single cells participating in post-pneumonectomy lung growth, we used laser microdissection, enzymatic digestion and microfluidic isolation. Single-cell transcriptional analysis of the murine alveolar duct cells was performed using the C1 integrated fluidic circuit (Fluidigm) and a custom PCR panel designed for lung growth and repair genes. The multi-dimensional data set was analyzed using visualization software based on the tSNE algorithm. The analysis identified 6 cell clusters; 1 cell cluster was present only after pneumonectomy. This post-pneumonectomy cluster was significantly less transcriptionally active than 3 other clusters and may represent a transitional cell population. A provisional cluster identity for 4 of the 6 cell clusters was obtained by embedding bulk transcriptional data into the tSNE analysis. The transcriptional pattern of the 6 clusters was further analyzed for genes associated with lung repair, matrix production, and angiogenesis. The data demonstrated that multiple cell-types (clusters) transcribed genes linked to these basic functions. We conclude that the coordinated gene expression across multiple cell clusters is likely a response to a shared regenerative microenvironment within the subpleural alveolar ducts.

RNA sequence analysis reveals macroscopic somatic clonal expansion across normal tissues.

How somatic mutations accumulate in normal cells is poorly understood. A comprehensive analysis of RNA sequencing data from ~6700 samples across 29 normal tissues revealed multiple somatic variants, demonstrating that macroscopic clones can be found in many normal tissues. We found that sun-exposed skin, esophagus, and lung have a higher mutation burden than other tested tissues, which suggests that environmental factors can promote somatic mosaicism. Mutation burden was associated with both age and tissue-specific cell proliferation rate, highlighting that mutations accumulate over both time and number of cell divisions. Finally, normal tissues were found to harbor mutations in known cancer genes and hotspots. This study provides a broad view of macroscopic clonal expansion in human tissues, thus serving as a foundation for associating clonal expansion with environmental factors, aging, and risk of disease.

Aligning Single-Cell Developmental and Reprogramming Trajectories Identifies Molecular Determinants of Myogenic Reprogramming Outcome.

Cellular reprogramming through manipulation of defined factors holds great promise for large-scale production of cell types needed for use in therapy and for revealing principles of gene regulation. However, most reprogramming systems are inefficient, converting only a fraction of cells to the desired state. Here, we analyze MYOD-mediated reprogramming of human fibroblasts to myotubes, a well-characterized model system for direct conversion by defined factors, at pseudotemporal resolution using single-cell RNA-seq. To expose barriers to efficient conversion, we introduce a novel analytic technique, trajectory alignment, which enables quantitative comparison of gene expression kinetics across two biological processes. Reprogrammed cells navigate a trajectory with branch points that correspond to two alternative decision points, with cells that select incorrect branches terminating at aberrant or incomplete reprogramming outcomes. Analysis of these branch points revealed insulin and BMP signaling as crucial molecular determinants of reprogramming. Single-cell trajectory alignment enables rigorous quantitative comparisons between biological trajectories found in diverse processes in development, reprogramming, and other contexts.

Novel mutations in the NRL gene and associated clinical findings in patients with dominant retinitis pigmentosa.

OBJECTIVES: To search for mutations in the neural retina leucine zipper (NRL) gene in patients with dominant retinitis pigmentosa and to compare the severity of disease in these patients with that observed previously in patients with dominant rhodopsin mutations. METHODS: Single-strand conformation analysis was used to survey 189 unrelated patients for mutations. The available relatives of index patients with mutations were also evaluated. In our clinical examination of patients, we measured visual acuity, final dark-adaptation threshold equivalent visual field diameter, and electroretinogram amplitudes among other parameters of visual function. We compared the clinical findings with those obtained earlier from similar evaluations of a group of 39 patients with the dominant rhodopsin mutation Pro23His and a group of 25 patients with the dominant rhodopsin mutation Pro347Leu. RESULTS: We identified 3 novel missense mutations in a total of 4 unrelated patients with dominant retinitis pigmentosa: Ser50Pro, Ser50Leu (2 patients), and Pro51Thr. Each mutation cosegregated with dominant retinitis pigmentosa. None of these mutations were found among 91 unrelated control individuals. The visual acuities among the 4 index patients and 3 relatives with NRL mutations who were clinically evaluated ranged from 20/20 (in a 9-year-old patient) to 20/200 (in a 73-year-old patient). All patients had bone-spicule pigment deposits in their fundi. Average rod-plus-cone and cone-isolated electroretinogram amplitudes were both decreased by 99% or more compared with normal amplitudes. The dark-adaptation thresholds, equivalent visual field diameters, and electroretinogram amplitudes (all corrected for age and refractive error) indicated that the disease caused by the NRL mutations was more severe than that caused by the dominant rhodopsin mutation Pro23His and was similar in severity to that produced by the rhodopsin mutation Pro347Leu. CONCLUSION: The 3 novel NRL mutations we discovered bring the total number of reported mutations in this gene to 6. Five of the 6 mutations affect residues 50 or 51, suggesting that these residues are important in a structural or functional domain of the encoded protein. CLINICAL RELEVANCE: Rod and cone function is affected to a similar degree in patients with these mutations. The disease caused by NRL mutations found in this study appears to be more severe than that caused by the rhodopsin mutation Pro23His and is similar in severity to that caused by the rhodopsin mutation Pro347Leu, even after correcting for age.

Novel mutations in the cellular retinaldehyde-binding protein gene (RLBP1) associated with retinitis punctata albescens: evidence of interfamilial genetic heterogeneity and fundus changes in heterozygotes.

OBJECTIVE: To evaluate the molecular genetic defects associated with retinitis punctata albescens (RPA) in 5 patients from 3 families with this disease. METHODS: We examined 3 probands and 2 clinically affected relatives with RPA. Clinical examinations included best-corrected visual acuity, visual field testing, electroretinography, dilated fundus examination, and fundus photography. Leukocyte DNA was analyzed for mutations in the exons of the genes encoding cellular retinaldehyde-binding protein 1 (RLBP1), 11-cis-retinol dehydrogenase (RDH5), interphotoreceptor retinoid-binding protein (RBP3), and photoreceptor all-trans-retinol dehydrogenase (RDH8). Not all patients were evaluated for mutations in each gene. The exons were individually amplified and screened for mutations by single-stranded conformational polymorphism analysis or direct genomic sequencing. RESULTS: The 3 probands had similar clinical findings, including a history of poor night vision, the presence of punctate white deposits in the retina, and substantially reduced or absent rod responses on electroretinogram testing. One of the probands (patient 2:III:2) had 2 novel mutations in the RLBP1 gene (Arg151Trp and Gly31[2-base pair deletion], [GGA-->G-]). Segregation analysis showed that the 2 mutations were allelic and that the patient was a compound heterozygote. Both parents of the proband manifested round white deposits in the retina. The other 2 probands had no detected pathogenic mutations in RLBP1 or in the other 3 genes evaluated. CONCLUSIONS: The identification of novel RLBP1 mutations in 1 of our 3 probands, all with RPA, is further evidence of genetic (nonallelic) heterogeneity in this disease. The presence of round white deposits in the retina may be observed in those heterozygous for RLBP1. Clinical Relevance Patients with a clinical presentation of RPA can have genetically different mutations. Drusen-like lesions may be observed in heterozygotes in families with this disease and a mutation in RLBP1.

The dynamics and regulators of cell fate decisions are revealed by pseudotemporal ordering of single cells.

Defining the transcriptional dynamics of a temporal process such as cell differentiation is challenging owing to the high variability in gene expression between individual cells. Time-series gene expression analyses of bulk cells have difficulty distinguishing early and late phases of a transcriptional cascade or identifying rare subpopulations of cells, and single-cell proteomic methods rely on a priori knowledge of key distinguishing markers. Here we describe Monocle, an unsupervised algorithm that increases the temporal resolution of transcriptome dynamics using single-cell RNA-Seq data collected at multiple time points. Applied to the differentiation of primary human myoblasts, Monocle revealed switch-like changes in expression of key regulatory factors, sequential waves of gene regulation, and expression of regulators that were not known to act in differentiation. We validated some of these predicted regulators in a loss-of function screen. Monocle can in principle be used to recover single-cell gene expression kinetics from a wide array of cellular processes, including differentiation, proliferation and oncogenic transformation.

Somatic mutation of CDKN1B in small intestine neuroendocrine tumors.

The diagnosed incidence of small intestine neuroendocrine tumors (SI-NETs) is increasing, and the underlying genomic mechanisms have not yet been defined. Using exome- and genome-sequence analysis of SI-NETs, we identified recurrent somatic mutations and deletions in CDKN1B, the cyclin-dependent kinase inhibitor gene, which encodes p27. We observed frameshift mutations of CDKN1B in 14 of 180 SI-NETs, and we detected hemizygous deletions encompassing CDKN1B in 7 out of 50 SI-NETs, nominating p27 as a tumor suppressor and implicating cell cycle dysregulation in the etiology of SI-NETs.

A genome-wide approach accounting for body mass index identifies genetic variants influencing fasting glycemic traits and insulin resistance.

Recent genome-wide association studies have described many loci implicated in type 2 diabetes (T2D) pathophysiology and ß-cell dysfunction but have contributed little to the understanding of the genetic basis of insulin resistance. We hypothesized that genes implicated in insulin resistance pathways might be uncovered by accounting for differences in body mass index (BMI) and potential interactions between BMI and genetic variants. We applied a joint meta-analysis approach to test associations with fasting insulin and glucose on a genome-wide scale. We present six previously unknown loci associated with fasting insulin at P < 5 × 10(-8) in combined discovery and follow-up analyses of 52 studies comprising up to 96,496 non-diabetic individuals. Risk variants were associated with higher triglyceride and lower high-density lipoprotein (HDL) cholesterol levels, suggesting a role for these loci in insulin resistance pathways. The discovery of these loci will aid further characterization of the role of insulin resistance in T2D pathophysiology.

Pathway analysis following association study.

Genome-wide association studies often emphasize single-nucleotide polymorphisms with the smallest p-values with less attention given to single-nucleotide polymorphisms not ranked near the top. We suggest that gene pathways contain valuable information that can enable identification of additional associations. We used gene set information to identify disease-related pathways using three methods: gene set enrichment analysis (GSEA), empirical enrichment p-values, and Ingenuity pathway analysis (IPA). Association tests were performed for common single-nucleotide polymorphisms and aggregated rare variants with traits Q1 and Q4. These pathway methods were evaluated by type I error, power, and the ranking of the VEGF pathway, the gene set used in the simulation model. GSEA and IPA had high power for detecting the VEGF pathway for trait Q1 (91.2% and 93%, respectively). These two methods were conservative with deflated type I errors (0.0083 and 0.0072, respectively). The VEGF pathway ranked 1 or 2 in 123 of 200 replicates using IPA and ranked among the top 5 in 114 of 200 replicates for GSEA. The empirical enrichment method had lower power and higher type I error. Thus pathway analysis approaches may be useful in identifying biological pathways that influence disease outcomes.