Search

Abstract Whole genome sequencing (WGS) at high-depth (30X) allows the accurate discovery of variants in the coding and non-coding DNA regions and helps elucidate the genetic underpinnings of human health and diseases. Yet, due to the prohibitive cost of high-depth WGS, most large-scale genetic association studies use genotyping arrays or high-depth whole exome sequencing (WES). Here we propose a cost-effective method which we call “Whole Exome Genome Sequencing” (WEGS), that combines low-depth WGS and high-depth WES with up to 8 samples pooled and sequenced simultaneously (multiplexed). We experimentally assess the performance of WEGS with four different depth of coverage and sample multiplexing configurations. We show that the optimal WEGS configurations are 1.7–2.0 times cheaper than standard WES (no-plexing), 1.8–2.1 times cheaper than high-depth WGS, reach similar recall and precision rates in detecting coding variants as WES, and capture more population-specific variants in the rest of the genome that are difficult to recover when using genotype imputation methods. We apply WEGS to 862 patients with peripheral artery disease and show that it directly assesses more known disease-associated variants than a typical genotyping array and thousands of non-imputable variants per disease-associated locus.

Abstract The persistence of SARS-CoV-2 despite the development of vaccines and a degree of herd immunity is partly due to viral evolution reducing vaccine and treatment efficacy. Serial infections of wild-type (WT) SARS-CoV-2 in Balb/c mice yield mouse-adapted strains with greater infectivity and mortality. We investigate if passaging unmodified B.1.351 (Beta) and B.1.617.2 (Delta) 20 times in K18-ACE2 mice, expressing the human ACE2 receptor, in a BSL-3 laboratory without selective pressures, drives human health-relevant evolution and if evolution is lineage-dependent. Late-passage virus causes more severe disease, at organism and lung tissue scales, with late-passage Delta demonstrating antibody resistance and interferon suppression. This resistance co-occurs with a de novo spike S371F mutation, linked with both traits. S371F, an Omicron-characteristic mutation, is co-inherited at times with spike E1182G per Nanopore sequencing, existing in different within-sample viral variants at others. Both S371F and E1182G are linked to mammalian GOLGA7 and ZDHHC5 interactions, which mediate viral-cell entry and antiviral response. This study demonstrates SARS-CoV-2’s tendency to evolve with phenotypic consequences, its evolution varying by lineage, and suggests non-dominant quasi-species contribution.

Abstract The American eel (Anguilla rostrata) has long been regarded as a panmictic fish and has been confirmed as such in the northern part of its range. In this paper, we tested for the first time whether panmixia extends to the tropical range of the species. To do so, we first assembled a reference genome (975 Mbp, 19 chromosomes) combining long (PacBio and Nanopore and short (Illumina paired‐end) reads technologies to support both this study and future research. To test for population structure, we estimated genotype likelihoods from low‐coverage whole‐genome sequencing of 460 American eels, collected at 21 sampling sites (in seven geographic regions) ranging from Canada to Trinidad and Tobago. We estimated genetic distance between regions, performed ADMIXTURE‐like clustering analysis and multivariate analysis, and found no evidence of population structure, thus confirming that panmixia extends to the tropical range of the species. In addition, two genomic regions with putative inversions were observed, both geographically widespread and present at similar frequencies in all regions. We discuss the implications of lack of genetic population structure for the species. Our results are key for the future genomic research in the American eel and the implementation of conservation measures throughout its geographic range. Additionally, our results can be applied to fisheries management and aquaculture of the species.

During the COVID-19 pandemic, the monitoring of SARS-CoV-2 RNA in wastewater was used to track the evolution and emergence of variant lineages and gauge infection levels in the community, informing appropriate public health responses without relying solely on clinical testing. As more sublineages were discovered, it increased the difficulty in identifying distinct variants in a mixed population sample, particularly those without a known lineage. Here, we compare the sequencing technology from Illumina and from Oxford Nanopore Technologies, in order to determine their efficacy at detecting variants of differing abundance, using 248 wastewater samples from various Quebec and Ontario cities. Our study used two analytical approaches to identify the main variants in the samples: the presence of signature and marker mutations and the co-occurrence of signature mutations within the same amplicon. We observed that each sequencing method detected certain variants at different frequencies as each method preferentially detects mutations of distinct variants. Illumina sequencing detected more mutations with a predominant lineage that is in low abundance across the population or unknown for that time period, while Nanopore sequencing had a higher detection rate of mutations that are predominantly found in the high abundance B.1.1.7 (Alpha) lineage as well as a higher sequencing rate of co-occurring mutations in the same amplicon. We present a workflow that integrates short-read and long-read sequencing to improve the detection of SARS-CoV-2 variant lineages in mixed population samples, such as wastewater.

Summary: The basal breast cancer subtype is enriched for triple-negative breast cancer (TNBC) and displays consistent large chromosomal deletions. Here, we characterize evolution and maintenance of chromosome 4p (chr4p) loss in basal breast cancer. Analysis of The Cancer Genome Atlas data shows recurrent deletion of chr4p in basal breast cancer. Phylogenetic analysis of a panel of 23 primary tumor/patient-derived xenograft basal breast cancers reveals early evolution of chr4p deletion. Mechanistically we show that chr4p loss is associated with enhanced proliferation. Gene function studies identify an unknown gene, C4orf19, within chr4p, which suppresses proliferation when overexpressed—a member of the PDCD10-GCKIII kinase module we name PGCKA1. Genome-wide pooled overexpression screens using a barcoded library of human open reading frames identify chromosomal regions, including chr4p, that suppress proliferation when overexpressed in a context-dependent manner, implicating network interactions. Together, these results shed light on the early emergence of complex aneuploid karyotypes involving chr4p and adaptive landscapes shaping breast cancer genomes.

Abstract Background Expression quantitative trait loci (eQTL) studies provide insights into regulatory mechanisms underlying disease risk. Expanding studies of gene regulation to underexplored populations and to medically relevant tissues offers potential to reveal yet unknown regulatory variants and to better understand disease mechanisms. Here, we performed eQTL mapping in subcutaneous (S) and visceral (V) adipose tissue from 106 Greek individuals (Greek Metabolic study, GM) and compared our findings to those from the Genotype-Tissue Expression (GTEx) resource. Results We identified 1,930 and 1,515 eGenes in S and V respectively, over 13% of which are not observed in GTEx adipose tissue, and that do not arise due to different ancestry. We report additional context-specific regulatory effects in genes of clinical interest (e.g. oncogene ST7) and in genes regulating responses to environmental stimuli (e.g. MIR21, SNX33). We suggest that a fraction of the reported differences across populations is due to environmental effects on gene expression, driving context-specific eQTLs, and suggest that environmental effects can determine the penetrance of disease variants thus shaping disease risk. We report that over half of GM eQTLs colocalize with GWAS SNPs and of these colocalizations 41% are not detected in GTEx. We also highlight the clinical relevance of S adipose tissue by revealing that inflammatory processes are upregulated in individuals with obesity, not only in V, but also in S tissue. Conclusions By focusing on an understudied population, our results provide further candidate genes for investigation regarding their role in adipose tissue biology and their contribution to disease risk and pathogenesis.

Aging leads to significant alteration in the gene expression of muscle stem cells. In vivo exposure of muscle stem cells from aged mice to a young niche environment restores the expression of a significant portion of age-altered genes in mice.

Abstract Background Juvenile Pilocytic Astrocytomas (JPAs) are one of the most common pediatric brain tumors, and they are driven by aberrant activation of the mitogen-activated protein kinase (MAPK) signaling pathway. RAF-fusions are the most common genetic alterations identified in JPAs, with the prototypical KIAA1549-BRAF fusion leading to loss of BRAF’s auto-inhibitory domain and subsequent constitutive kinase activation. JPAs are highly vascular and show pervasive immune infiltration, which can lead to low tumor cell purity in clinical samples. This can result in gene fusions that are difficult to detect with conventional omics approaches including RNA-Seq. Methods To this effect, we applied RNA-Seq as well as linked-read whole-genome sequencing and in situ Hi-C as new approaches to detect and characterize low-frequency gene fusions at the genomic, transcriptomic and spatial level. Results Integration of these datasets allowed the identification and detailed characterization of two novel BRAF fusion partners, PTPRZ1 and TOP2B, in addition to the canonical fusion with partner KIAA1549. Additionally, our Hi-C datasets enabled investigations of 3D genome architecture in JPAs which showed a high level of correlation in 3D compartment annotations between JPAs compared to other pediatric tumors, and high similarity to normal adult astrocytes. We detected interactions between BRAF and its fusion partners exclusively in tumor samples containing BRAF fusions. Conclusions We demonstrate the power of integrating multi-omic datasets to identify low frequency fusions and characterize the JPA genome at high resolution. We suggest that linked-reads and Hi-C could be used in clinic for the detection and characterization of JPAs.

Genomes from defective interfering (DI) particles following serial passaging of SARS-CoV-2 reveal a fusion protein that attenuates viral replication. Synthetic, recombinant DI genomes are designed to interfere with SARS-CoV-2 replication.

IntroductionMedulloblastoma is the most common type of malignant pediatric brain tumor with group 4 medulloblastomas (G4 MBs) accounting for 40% of cases. However, the molecular mechanisms that underlie this subgroup are still poorly understood. Point mutations are detected in a large number of genes at low incidence per gene while the detection of complex structural variants in recurrently affected genes typically requires the application of long-read technologies.MethodsHere, we applied linked-read sequencing, which combines the long-range genome information of long-read sequencing with the high base pair accuracy of short read sequencing and very low sample input requirements.ResultsWe demonstrate the detection of complex structural variants and point mutations in these tumors, and, for the first time, the detection of extrachromosomal DNA (ecDNA) with linked-reads. We provide further evidence for the high heterogeneity of somatic mutations in G4 MBs and add new complex events associated with it.DiscussionWe detected several enhancer-hijacking events, an ecDNA containing the MYCN gene, and rare structural rearrangements, such a chromothripsis in a G4 medulloblastoma, chromoplexy involving 8 different chromosomes, a TERT gene rearrangement, and a PRDM6 duplication.

AbstractGenetic sexing strains (GSS) are an important tool in support of sterile insect technique (SIT) applications against insect pests and disease vectors. The yet unknown temperature-sensitive lethaltslwhite pupaewpCeratitis capitatatsltslwpZwtsltsltslDrosophila melanogastertsl

Not all familial ovarian cancer (OC) cases are explained by pathogenic germline variants in known risk genes. A candidate gene approach involving DNA repair pathway genes was applied to identify rare recurring pathogenic variants in familial OC cases not associated with known OC risk genes from a population exhibiting genetic drift. Whole exome sequencing (WES) data of 15 OC cases from 13 families tested negative for pathogenic variants in known OC risk genes were investigated for candidate variants in 468 DNA repair pathway genes. Filtering and prioritization criteria were applied to WES data to select top candidates for further analyses. Candidates were genotyped in ancestry defined study groups of 214 familial and 998 sporadic OC or breast cancer (BC) cases and 1025 population-matched controls and screened for additional carriers in 605 population-matched OC cases. The candidate genes were also analyzed in WES data from 937 familial or sporadic OC cases of diverse ancestries. Top candidate variants in ERCC5, EXO1, FANCC, NEIL1 and NTHL1 were identified in 5/13 (39%) OC families. Collectively, candidate variants were identified in 7/435 (1.6%) sporadic OC cases and 1/566 (0.2%) sporadic BC cases versus 1/1025 (0.1%) controls. Additional carriers were identified in 6/605 (0.9%) OC cases. Tumour DNA from ERCC5, NEIL1 and NTHL1 variant carriers exhibited loss of the wild-type allele. Carriers of various candidate variants in these genes were identified in 31/937 (3.3%) OC cases of diverse ancestries versus 0-0.004% in cancer-free controls. The strategy of applying a candidate gene approach in a population exhibiting genetic drift identified new candidate OC predisposition variants in DNA repair pathway genes.

Microglia subpopulations may differentially contribute to pain. Here, the authors show that peripheral nerve injury induces time- and sex-specific transcriptional changes in mouse spinal cord microglia subpopulations and that ApoE is linked to neuropathic pain hypersensitivity in mice and humans.

Abstract Background Familial ovarian cancer (OC) cases not harbouring pathogenic variants in either of the BRCA1 and BRCA2 OC-predisposing genes, which function in homologous recombination (HR) of DNA, could involve pathogenic variants in other DNA repair pathway genes. Methods Whole exome sequencing was used to identify rare variants in HR genes in a BRCA1 and BRCA2 pathogenic variant negative OC family of French Canadian (FC) ancestry, a population exhibiting genetic drift. OC cases and cancer-free individuals from FC and non-FC populations were investigated for carrier frequency of FANCI c.1813C>T; p.L605F, the top-ranking candidate. Gene and protein expression were investigated in cancer cell lines and tissue microarrays, respectively. Results In FC subjects, c.1813C>T was more common in familial (7.1%, 3/42) than sporadic (1.6%, 7/439) OC cases (P = 0.048). Carriers were detected in 2.5% (74/2950) of cancer-free females though female/male carriers were more likely to have a first-degree relative with OC (121/5249, 2.3%; Spearman correlation = 0.037; P = 0.011), suggesting a role in risk. Many of the cancer-free females had host factors known to reduce risk to OC which could influence cancer risk in this population. There was an increased carrier frequency of FANCI c.1813C>T in BRCA1 and BRCA2 pathogenic variant negative OC families, when including the discovery family, compared to cancer-free females (3/23, 13%; OR = 5.8; 95%CI = 1.7–19; P = 0.005). In non-FC subjects, 10 candidate FANCI variants were identified in 4.1% (21/516) of Australian OC cases negative for pathogenic variants in BRCA1 and BRCA2, including 10 carriers of FANCI c.1813C>T. Candidate variants were significantly more common in familial OC than in sporadic OC (P = 0.04). Localization of FANCD2, part of the FANCI-FANCD2 (ID2) binding complex in the Fanconi anaemia (FA) pathway, to sites of induced DNA damage was severely impeded in cells expressing the p.L605F isoform. This isoform was expressed at a reduced level, destabilized by DNA damaging agent treatment in both HeLa and OC cell lines, and exhibited sensitivity to cisplatin but not to a poly (ADP-ribose) polymerase inhibitor. By tissue microarray analyses, FANCI protein was consistently expressed in fallopian tube epithelial cells and only expressed at low-to-moderate levels in 88% (83/94) of OC samples. Conclusions This is the first study to describe candidate OC variants in FANCI, a member of the ID2 complex of the FA DNA repair pathway. Our data suggest that pathogenic FANCI variants may modify OC risk in cancer families.

Abstract Background Québec was the Canadian province most impacted by COVID-19, with 401,462 cases as of September 24th, 2021, and 11,347 deaths due mostly to a very severe first pandemic wave. In April 2020, we assembled the Coronavirus Sequencing in Québec (CoVSeQ) consortium to sequence SARS-CoV-2 genomes in Québec to track viral introduction events and transmission within the province. Methods Using genomic epidemiology, we investigated the arrival of SARS-CoV-2 to Québec. We report 2921 high-quality SARS-CoV-2 genomes in the context of > 12,000 publicly available genomes sampled globally over the first pandemic wave (up to June 1st, 2020). By combining phylogenetic and phylodynamic analyses with epidemiological data, we quantify the number of introduction events into Québec, identify their origins, and characterize the spatiotemporal spread of the virus. Results Conservatively, we estimated approximately 600 independent introduction events, the majority of which happened from spring break until 2 weeks after the Canadian border closed for non-essential travel. Subsequent mass repatriations did not generate large transmission lineages (> 50 sequenced cases), likely due to mandatory quarantine measures in place at the time. Consistent with common spring break and “snowbird” destinations, most of the introductions were inferred to have originated from Europe via the Americas. Once introduced into Québec, viral lineage sizes were overdispersed, with a few lineages giving rise to most infections. Consistent with founder effects, the earliest lineages to arrive tended to spread most successfully. Fewer than 100 viral introductions arrived during spring break, of which 7–12 led to the largest transmission lineages of the first wave (accounting for 52–75% of all sequenced infections). These successful transmission lineages dispersed widely across the province. Transmission lineage size was greatly reduced after March 11th, when a quarantine order for returning travellers was enacted. While this suggests the effectiveness of early public health measures, the biggest transmission lineages had already been ignited prior to this order. Conclusions Combined, our results reinforce how, in the absence of tight travel restrictions or quarantine measures, fewer than 100 viral introductions in a week can ensure the establishment of extended transmission chains.

Full-length transcript sequencing remains a main goal of RNA sequencing. However, even the application of long-read sequencing technologies such as Oxford Nanopore Technologies still fail to yield full-length transcript sequencing for a significant portion of sequenced reads. Since these technologies can sequence reads that are far longer than the longest known processed transcripts, the lack of efficiency to obtain full-length transcripts from good quality RNAs stems from library preparation inefficiency rather than the presence of degraded RNA molecules. It has previously been shown that addition of inverted terminal repeats in cDNA during reverse transcription followed by single-primer PCR creates a PCR suppression effect that prevents amplification of short molecules thus enriching the library for longer transcripts. We adapted this method for Nanopore cDNA library preparation and show that not only is PCR efficiency increased but gene body coverage is dramatically improved. The results show that implementation of this simple strategy will result in better quality full-length RNA sequencing data and make full-length transcript sequencing possible for most of sequenced reads.

Abstract The Vero cell line is the most used continuous cell line for viral vaccine manufacturing with more than 40 years of accumulated experience in the vaccine industry. Additionally, the Vero cell line has shown a high affinity for infection by MERS-CoV, SARS-CoV, and recently SARS-CoV-2, emerging as an important discovery and screening tool to support the global research and development efforts in this COVID-19 pandemic. However, the lack of a reference genome for the Vero cell line has limited our understanding of host–virus interactions underlying such affinity of the Vero cell towards key emerging pathogens, and more importantly our ability to redesign high-yield vaccine production processes using Vero genome editing. In this paper, we present an annotated highly contiguous 2.9 Gb assembly of the Vero cell genome. In addition, several viral genome insertions, including Adeno-associated virus serotypes 3, 4, 7, and 8, have been identified, giving valuable insights into quality control considerations for cell-based vaccine production systems. Variant calling revealed that, in addition to interferon, chemokines, and caspases-related genes lost their functions. Surprisingly, the ACE2 gene, which was previously identified as the host cell entry receptor for SARS-CoV and SARS-CoV-2, also lost function in the Vero genome due to structural variations.

Nephric duct (ND)-derived ureteric buds (UB) form the kidney collecting duct system, while ureteric tips promote nephron formation. Here the authors use single-cell RNA-seq and introduce Cluster RNA-seq to identify four progenitor populations in developing ND/UB regulated by the transcription factors Tfap2a/b and Gata3.

Abstract The olive fruit fly, Bactrocera oleae, is the most important pest for the olive fruit but lacks adequate transcriptomic characterization that could aid in molecular control approaches. We apply nanopore long-read RNA-seq with internal RNA standards allowing absolute transcript quantification to analyze transcription dynamics during early embryo development for the first time in this organism. Sequencing on the MinION platform generated over 31 million reads. Over 50% of the expressed genes had at least one read covering its entire length validating our full-length approach. We generated a de novo transcriptome assembly and identified 1768 new genes and a total of 79,810 isoforms; a fourfold increase in transcriptome diversity compared to the current NCBI predicted transcriptome. Absolute transcript quantification per embryo allowed an insight into the dramatic re-organization of maternal transcripts. We further identified Zelda as a possible regulator of early zygotic genome activation in B. oleae and provide further insights into the maternal-to-zygotic transition. These data show the utility of long-read RNA in improving characterization of non-model organisms that lack a fully annotated genome, provide potential targets for sterile insect technic approaches, and provide the first insight into the transcriptome landscape of the developing olive fruit fly embryo.

Sonic Hedgehog medulloblastoma (Shh-MB) comprises four subtypes each with distinct clinical traits. Here the authors characterize the genome, transcriptome, and methylome of Shh-MB subtypes, revealing a complex fusion landscape and the molecular convergence of MYCN and cAMP signaling pathways.

The white pupae (wp) phenotype has been used for decades to selectively remove females of tephritid species in genetic sexing, but the determining gene is unknown. Here, the authors show that wp phenotype is produced by parallel mutations in a Major Facilitator Superfamily domain containing gene across multiple species.

Purpose The Canadian Longitudinal Study on Aging (CLSA) Comprehensive cohort was established to provide unique opportunities to study the genetic and environmental contributions to human disease as well as ageing process. The aim of this report was to describe the genomic data included in CLSA.Participants A total of 26 622 individuals from the CLSA Comprehensive cohort of men and women aged 45–85 recruited between 2010 and 2015 underwent genome-wide genotyping of DNA samples collected from blood. Comprehensive quality control metrics were measured for genetic markers and samples, respectively. The genotypes were imputed to the TOPMed reference panel. Sex chromosome abnormalities were identified by copy number profiling. Classical human leukocyte antigen gene haplotypes were imputed at two-field (four-digit).Findings to date Of the 26 622 genotyped participants, 24 655 (92.6%) were identified as having European ancestry. These genomic data were linked to physical, lifestyle, medical, economic, environmental and psychosocial factors collected longitudinally in CLSA. The combined analysis, including CLSA genomic data, uncovered over 100 novel loci associated with key parameters to define glaucoma. The CLSA genomic dataset validated the contribution of a polygenic risk score to screen individuals with high fracture risk. It is also a valuable resource to directly identify common genetic variations associated with conditions related to complex traits. Taking advantage of the comprehensive interview and physical information collected in CLSA, this genomic dataset has been linked to psychosocial factors to investigate both the independent and interactive effects on cardiovascular disease.Future plans The CLSA overall is ongoing. Follow-up data will continue to be collected from participants in the current genomic subcohort, including the DNA methylation and metabolomic data. Ongoing studies focus on elucidating the role of genetic factors in cognitive decline and cardiovascular diseases. This genomic data resource is available on request through the CLSA data access application process.

Noncoding transcripts contribute to the adaptation of cellular processes to oxygen levels. Here the authors characterize a hypoxia responsive lncRNA lincNORS and show that it has a role in cellular sterol homeostasis.

Glioblastoma is thought to arise from neural stem cells. Here, to investigate this, the authors use single-cell RNA-sequencing to compare glioblastoma to the fetal human brain, and find a similarity between glial progenitor cells and a subpopulation of glioblastoma cells.

Sample index hopping results in various artefacts in multiplexed scRNA-seq experiments. Here, the authors introduce a statistical model to estimate sample index hopping rate in droplet-based scRNA-seq data and show that artifacts can be corrected by purging phantom molecules from the data.

Abstract Background The olive fruit fly, Bactrocera oleae, is the most important pest in the olive fruit agribusiness industry. This is because female flies lay their eggs in the unripe fruits and upon hatching the larvae feed on the fruits thus destroying them. The lack of a high-quality genome and other genomic and transcriptomic data has hindered progress in understanding the fly’s biology and proposing alternative control methods to pesticide use. Results Genomic DNA was sequenced from male and female Demokritos strain flies, maintained in the laboratory for over 45 years. We used short-, mate-pair-, and long-read sequencing technologies to generate a combined male-female genome assembly (GenBank accession GCA_001188975.2). Genomic DNA sequencing from male insects using 10x Genomics linked-reads technology followed by mate-pair and long-read scaffolding and gap-closing generated a highly contiguous 489 Mb genome with a scaffold N50 of 4.69 Mb and L50 of 30 scaffolds (GenBank accession GCA_001188975.4). RNA-seq data generated from 12 tissues and/or developmental stages allowed for genome annotation. Short reads from both males and females and the chromosome quotient method enabled identification of Y-chromosome scaffolds which were extensively validated by PCR. Conclusions The high-quality genome generated represents a critical tool in olive fruit fly research. We provide an extensive RNA-seq data set, and genome annotation, critical towards gaining an insight into the biology of the olive fruit fly. In addition, elucidation of Y-chromosome sequences will advance our understanding of the Y-chromosome’s organization, function and evolution and is poised to provide avenues for sterile insect technique approaches.

The sex chromosome complement, XX or XY, determines sexual differentiation of the gonadal primordium into a testis or an ovary, which in turn directs differentiation of the germ cells into sperm and oocytes, respectively, in eutherian mammals. When the X monosomy or XY sex reversal occurs, XO and XY females exhibit subfertility and infertility in the mouse on the C57BL/6J genetic background, suggesting that functional germ cell differentiation requires the proper sex chromosome complement. Using these mouse models, we asked how the sex chromosome complement affects gene transcription in the oocytes during follicular growth. An oocyte accumulates cytoplasmic components such as mRNAs and proteins during follicular growth to support subsequent meiotic progression, fertilization, and early embryonic development without de novo transcription. However, how gene transcription is regulated during oocyte growth is not well understood. Our results revealed that XY oocytes became abnormal in chromatin configuration, mitochondria distribution, and de novo transcription compared to XX or XO oocytes near the end of growth phase. Therefore, we compared transcriptomes by RNA-sequencing among the XX, XO, and XY oocytes of 50–60 µm in diameter, which were still morphologically comparable. The results showed that the X chromosome dosage limited the X-linked and autosomal gene transcript levels in XO oocytes whereas many genes were transcribed from the Y chromosome and made the transcriptome in XY oocytes closer to that in XX oocytes. We then compared the transcript levels of 3 X-linked, 3 Y-linked and 2 autosomal genes in the XX, XO, and XY oocytes during the entire growth phase as well as at the end of growth phase using quantitative RT-PCR. The results indicated that the transcript levels of most genes increased with oocyte growth while largely maintaining the X chromosome dosage dependence. Near the end of growth phase, however, transcript levels of some X-linked genes did not increase in XY oocytes as much as XX or XO oocytes, rendering their levels much lower than those in XX oocytes. Thus, XY oocytes established a distinct transcriptome at the end of growth phase, which may be associated with abnormal chromatin configuration and mitochondria distribution.

AbstractDespite the commercial importance of Greenland Halibut (Reinhardtius hippoglossoides)

Population sequencing often requires collaboration across a distributed network of sequencing centers for the timely processing of thousands of samples. In such massive efforts, it is important that participating scientists can be confident that the accuracy of the sequence data produced is not affected by which center generates the data. A study was conducted across three established sequencing centers, located in Montreal, Toronto, and Vancouver, constituting Canada’s Genomics Enterprise (www.cgen.ca). Whole genome sequencing was performed at each center, on three genomic DNA replicates from three well-characterized cell lines. Secondary analysis pipelines employed by each site were applied to sequence data from each of the sites, resulting in three datasets for each of four variables (cell line, replicate, sequencing center, and analysis pipeline), for a total of 81 datasets. These datasets were each assessed according to multiple quality metrics including concordance with benchmark variant truth sets to assess consistent quality across all three conditions for each variable. Three-way concordance analysis of variants across conditions for each variable was performed. Our results showed that the variant concordance between datasets differing only by sequencing center was similar to the concordance for datasets differing only by replicate, using the same analysis pipeline. We also showed that the statistically significant differences between datasets result from the analysis pipeline used, which can be unified and updated as new approaches become available. We conclude that genome sequencing projects can rely on the quality and reproducibility of aggregate data generated across a network of distributed sites.

In this study, we used single-cell transcriptomic analysis to identify new specific biomarkers for nucleus pulposus (NP) and inner annulus fibrosis (iAF) cells, and to define cell populations within non-degenerating (nD) and degenerating (D) human intervertebral discs (IVD) of the same individual. Cluster analysis based on differential gene expression delineated 14 cell clusters. Gene expression profiles at single-cell resolution revealed the potential functional differences linked to degeneration, and among NP and iAF subpopulations. GO and KEGG analyses discovered molecular functions, biological processes, and transcription factors linked to cell type and degeneration state. We propose two lists of biomarkers, one as specific cell type, including C2orf40, MGP, MSMP, CD44, EIF1, LGALS1, RGCC, EPYC, HILPDA, ACAN, MT1F, CHI3L1, ID1, ID3 and TMED2. The second list proposes predictive IVD degeneration genes, including MT1G, SPP1, HMGA1, FN1, FBXO2, SPARC, VIM, CTGF, MGST1, TAF1D, CAPS, SPTSSB, S100A1, CHI3L2, PLA2G2A, TNRSF11B, FGFBP2, MGP, SLPI, DCN, MT-ND2, MTCYB, ADIRF, FRZB, CLEC3A, UPP1, S100A2, PRG4, COL2A1, SOD2 and MT2A. Protein and mRNA expression of MGST1, vimentin, SOD2 and SYF2 (p29) genes validated our scRNA-seq findings. Our data provide new insights into disc cells phenotypes and biomarkers of IVD degeneration that could improve diagnostic and therapeutic options.

RNA sequencing using next-generation sequencing technologies (NGS) is currently the standard approach for gene expression profiling, particularly for large-scale high-throughput studies. NGS technologies comprise high throughput, cost efficient short-read RNA-Seq, while emerging single molecule, long-read RNA-Seq technologies have enabled new approaches to study the transcriptome and its function. The emerging single molecule, long-read technologies are currently commercially available by Pacific Biosciences (PacBio) and Oxford Nanopore Technologies (ONT), while new methodologies based on short-read sequencing approaches are also being developed in order to provide long range single molecule level information—for example, the ones represented by the 10x Genomics linked read methodology. The shift toward long-read sequencing technologies for transcriptome characterization is based on current increases in throughput and decreases in cost, making these attractive for de novo transcriptome assembly, isoform expression quantification, and in-depth RNA species analysis. These types of analyses were challenging with standard short sequencing approaches, due to the complex nature of the transcriptome, which consists of variable lengths of transcripts and multiple alternatively spliced isoforms for most genes, as well as the high sequence similarity of highly abundant species of RNA, such as rRNAs. Here we aim to focus on single molecule level sequencing technologies and single-cell technologies that, combined with perturbation tools, allow the analysis of complete RNA species, whether short or long, at high resolution. In parallel, these tools have opened new ways in understanding gene functions at the tissue, network, and pathway levels, as well as their detailed functional characterization. Analysis of the epi-transcriptome, including RNA methylation and modification and the effects of such modifications on biological systems is now enabled through direct RNA sequencing instead of classical indirect approaches. However, many difficulties and challenges remain, such as methodologies to generate full-length RNA or cDNA libraries from all different species of RNAs, not only poly-A containing transcripts, and the identification of allele-specific transcripts due to current error rates of single molecule technologies, while the bioinformatics analysis on long-read data for accurate identification of 5′ and 3′ UTRs is still in development.

Information regarding transcriptome and metabolome has significantly contributed to identifying potential therapeutic targets for the management of a variety of cancers. Obesity has profound effects on both cancer cell transcriptome and metabolome that can affect the outcome of cancer therapy. The information regarding the potential effects of obesity on breast cancer (BC) transcriptome, metabolome, and its integration to identify novel pathways related to disease progression are still elusive. We assessed the whole blood transcriptome and serum metabolome, as circulating metabolites, of obese BC patients compared them with non-obese BC patients. In these patients' samples, 186 significant differentially expressed genes (DEGs) were identified, comprising 156 upregulated and 30 downregulated. The expressions of these gene were confirmed by qRT-PCR. Furthermore, 96 deregulated metabolites were identified as untargeted metabolomics in the same group of patients. These detected DEGs and deregulated metabolites enriched in many cellular pathways. Further investigation, by integration analysis between transcriptomics and metabolomics data at the pathway levels, revealed seven unique enriched pathways in obese BC patients when compared with non-obese BC patients, which may provide resistance for BC cells to dodge the circulating immune cells in the blood. In conclusion, this study provides information on the unique pathways altered at transcriptome and metabolome levels in obese BC patients that could provide an important tool for researchers and contribute further to knowledge on the molecular interaction between obesity and BC. Further studies are needed to confirm this and to elucidate the exact underlying mechanism for the effects of obesity on the BC initiation or/and progression.

The definitive characterization of common cancer stem cell (CSCs) subpopulations in breast cancer subtypes with distinct genotypic and phenotypic features remains an ongoing challenge. In this study, we have used a non-biased genome wide screening approach to identify transcriptional networks that may be specific to the CSC subpopulations in both luminal and basal breast cancer subtypes. In depth studies of three CSC-enriched breast cancer cell lines representing various subtypes of breast cancer revealed a striking hyperactivation of the mevalonate metabolic pathway in comparison to control cells. The upregulation of metabolic networks is a key feature of tumour cells securing growth and proliferative capabilities and dysregulated mevalonate metabolism has been associated with tumour malignancy and cellular transformation in breast cancer. Furthermore, accumulating evidence suggests that Simvastatin therapy, a mevalonate pathway inhibitor, could affect breast cancer progression and reduce breast cancer recurrence. When detailing the mevalonate pathway in breast cancer using a single-cell qPCR, we identified the mevalonate precursor enzyme, HMGCS1, as a specific marker of CSC-enriched subpopulations within both luminal and basal tumour subtypes. Down-regulation of HMGCS1 also decreased the CSC fraction and function in various model systems, suggesting that HMGCS1 is essential for CSC-activities in breast cancer in general. These data was supported by strong associations between HMGCS1 expression and aggressive features, such as high tumour grade, p53 mutations as well as ER-negativity in lymph node positive breast cancer. Importantly, loss of HMGCS1 also had a much more pronounced effect on CSC-activities compared to treatment with standard doses of Simvastatin. Taken together, this study highlights HMGCS1 as a potential gatekeeper for dysregulated mevalonate metabolism important for CSC-features in both luminal and basal breast cancer subtypes. Pharmacological inhibition of HMGCS1 could therefore be a superior novel treatment approach for breast cancer patients via additional CSC blocking functions.

Summary: Therapies targeting epidermal growth factor receptor (EGFR) have variable and unpredictable responses in breast cancer. Screening triple-negative breast cancer (TNBC) patient-derived xenografts (PDXs), we identify a subset responsive to EGFR inhibition by gefitinib, which displays heterogeneous expression of wild-type EGFR. Deep single-cell RNA sequencing of 3,500 cells from an exceptional responder identified subpopulations displaying distinct biological features, where elevated EGFR expression was significantly enriched in a mesenchymal/stem-like cellular cluster. Sorted EGFRhi subpopulations exhibited enhanced stem-like features, including ALDH activity, sphere-forming efficiency, and tumorigenic and metastatic potential. EGFRhi cells gave rise to EGFRhi and EGFRlo cells in primary and metastatic tumors, demonstrating an EGFR-dependent expansion and hierarchical state transition. Similar tumorigenic EGFRhi subpopulations were identified in independent PDXs, where heterogeneous EGFR expression correlated with gefitinib sensitivity. This provides new understanding for an EGFR-dependent hierarchy in TNBC and for patient stratification for therapeutic intervention. : Savage et al. demonstrate that sensitivity to EGFR inhibitor, gefitinib, in triple-negative breast cancer is paradoxically associated with EGFR heterogeneity. Using single-cell RNA sequencing in conjunction with functional assays, they identify TNBC tumors in which EGFR expression identifies cells with tumor-initiating capacity whose proliferative expansion is sensitive to EGFR inhibition. Keywords: breast cancer, tumor heterogeneity, patient-derived xenograft, single-cell RNA sequencing, EGFR inhibition, therapeutic response, tumor-initiating cell, cell hierarchy, BRCA1 mutation

Approximately 10% of breast cancer (BC) cases are hereditary BC (HBC), with HBC most commonly encountered in the context of hereditary breast and ovarian cancer (HBOC) syndrome. Although thousands of loss-of-function (LoF) alleles in over 20 genes have been associated with HBC susceptibility, the genetic etiology of approximately 50% of cases remains unexplained, even when polygenic risk models are considered. We focused on one of the least-studied European populations and applied whole-exome sequencing (WES) to 52 individuals from 17 Greek HBOC families, in which at least one patient was negative for known HBC risk variants. Initial screening revealed pathogenic variants in known cancer genes, including BARD1:p.Trp91* detected in a cancer-free individual, and MEN1:p.Glu260Lys detected in a BC patient. Gene- and variant-based approaches were applied to exome data to identify candidate risk variants outside of known risk genes. Findings were verified in a collection of Canadian HBOC patients of European ancestry (FBRCAX), in an independent group of Canadian BC patients (CHUM-BC) and controls (CARTaGENE), as well as in individuals from The Cancer Genome Atlas (TCGA) and the UK Biobank (UKB). Rare LoF variants were uncovered in MDM1 and NBEAL1 in Greek and Canadian HBOC patients. We also report prioritized missense variants SETBP1:c.4129G > C and C7orf34:c.248C > T. These variants comprise promising candidates whose role in cancer pathogenicity needs to be explored further.

Human hepatocellular carcinoma (HCC) is the most common and recurrent type of primary adult liver cancer without any effective therapy. Plant-derived compounds acting as anticancer agents can induce apoptosis by targeting several signaling pathways. Strigolactone (SL) is a novel class of phytohormone, whose analogues have been reported to possess anticancer properties on a panel of human cancer cell lines through inducing cell cycle arrest, destabilizing microtubular integrity, reducing damaged in the DNA repair machinery, and inducing apoptosis. In our previous study, we reported that a novel SL analogue, TIT3, reduces HepG2 cell proliferation, inhibits cell migration, and induces apoptosis. To decipher the mechanisms of TIT3-induced anticancer activity in HepG2, we performed RNA sequencing and the differential expression of genes was analyzed using different tools. RNA-Seq data showed that the genes responsible for microtubule organization such as TUBB, BUB1B, TUBG2, TUBGCP6, TPX2, and MAP7 were significantly downregulated. Several epigenetic modulators such as UHRF1, HDAC7, and DNMT1 were also considerably downregulated, and this effect was associated with significant upregulation of various proapoptotic genes including CASP3, TNF-α, CASP7, and CDKN1A (p21). Likewise, damaged DNA repair genes such as RAD51, RAD52, and DDB2 were also significantly downregulated. This study indicates that TIT3-induced antiproliferative and proapoptotic activities on HCC cells could involve several signaling pathways. Our results suggest that TIT3 might be a promising drug to treat HCC.

Cell labelling in a non-invasive and genetic engineering-free manner is crucial to cell biology applications. Here the authors develop cell labelling via photobleaching (CLaP), that uses laser illumination to label individual cells for genomics, cell-tracking, flow cytometry or ultra-microscopy.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Approximately 9% of cancer-related deaths are caused by colorectal cancer. Liver metastasis is a major factor for the high colorectal cancer mortality rate. However, the molecular mechanism underlying colorectal cancer liver metastasis remains unclear. Using a global and multidimensional integration approach, we studied sequencing data, protein-protein interactions, and regulation of transcription factor and non-coding RNAs in primary tumor samples and liver metastasis samples to unveil the potential bridging molecules and the regulators that functionally link different stages of colorectal cancer liver metastasis. Primary tumor samples and liver metastasis samples had modules with significant overlap and crosstalk from which we identified several bridging genes (e.g. KNG1 and COX5B), transcription factors (e.g. E2F4 and CDX2), microRNAs (e.g. miR-590-3p and miR-203) and lncRNAs (e.g. lincIRX5 and lincFOXF1) that may play an important role in the process of colorectal cancer liver metastasis. This study enhances our understanding of the genetic alterations and transcriptional regulation that drive the metastatic process, but also provides the methodology to guide the studies on other metastatic cancers.

Kaposi sarcoma-associated herpesvirus (KSHV) causes several tumors and hyperproliferative disorders. Hypoxia and hypoxia-inducible factors (HIFs) activate latent and lytic KSHV genes, and several KSHV proteins increase the cellular levels of HIF. Here, we used RNA sequencing, qRT-PCR, Taqman assays, and pathway analysis to explore the miRNA and mRNA response of uninfected and KSHV-infected cells to hypoxia, to compare this with the genetic changes seen in chronic latent KSHV infection, and to explore the degree to which hypoxia and KSHV infection interact in modulating mRNA and miRNA expression. We found that the gene expression signatures for KSHV infection and hypoxia have a 34% overlap. Moreover, there were considerable similarities between the genes up-regulated by hypoxia in uninfected (SLK) and in KSHV-infected (SLKK) cells. hsa-miR-210, a HIF-target known to have pro-angiogenic and anti-apoptotic properties, was significantly up-regulated by both KSHV infection and hypoxia using Taqman assays. Interestingly, expression of KSHV-encoded miRNAs was not affected by hypoxia. These results demonstrate that KSHV harnesses a part of the hypoxic cellular response and that a substantial portion of hypoxia-induced changes in cellular gene expression are induced by KSHV infection. Therefore, targeting hypoxic pathways may be a useful way to develop therapeutic strategies for KSHV-related diseases.

Interferon Regulatory Factor 5 (IRF5) plays a major role in setting up an inflammatory macrophage phenotype, but the molecular basis of its transcriptional activity is not fully understood. In this study, we conduct a comprehensive genome-wide analysis of IRF5 recruitment in macrophages stimulated with bacterial lipopolysaccharide and discover that IRF5 binds to regulatory elements of highly transcribed genes. Analysis of protein:DNA microarrays demonstrates that IRF5 recognizes the canonical IRF-binding (interferon-stimulated response element [ISRE]) motif in vitro. However, IRF5 binding in vivo appears to rely on its interactions with other proteins. IRF5 binds to a noncanonical composite PU.1:ISRE motif, and its recruitment is aided by RelA. Global gene expression analysis in macrophages deficient in IRF5 and RelA highlights the direct role of the RelA:IRF5 cistrome in regulation of a subset of key inflammatory genes. We map the RelA:IRF5 interaction domain and suggest that interfering with it would offer selective targeting of macrophage inflammatory activities.

Kaposi's sarcoma associated herpesvirus (KSHV) causes several tumors, including primary effusion lymphoma (PEL) and Kaposi's sarcoma (KS). Cellular and viral microRNAs (miRNAs) have been shown to play important roles in regulating gene expression. A better knowledge of the miRNA-mediated pathways affected by KSHV infection is therefore important for understanding viral infection and tumor pathogenesis. In this study, we used deep sequencing to analyze miRNA and cellular mRNA expression in a cell line with latent KSHV infection (SLKK) as compared to the uninfected SLK line. This approach revealed 153 differentially expressed human miRNAs, eight of which were independently confirmed by qRT-PCR. KSHV infection led to the dysregulation of ~15% of the human miRNA pool and most of these cellular miRNAs were down-regulated, including nearly all members of the 14q32 miRNA cluster, a genomic locus linked to cancer and that is deleted in a number of PEL cell lines. Furthermore, we identified 48 miRNAs that were associated with a total of 1,117 predicted or experimentally validated target mRNAs; of these mRNAs, a majority (73%) were inversely correlated to expression changes of their respective miRNAs, suggesting miRNA-mediated silencing mechanisms were involved in a number of these alterations. Several dysregulated miRNA-mRNA pairs may facilitate KSHV infection or tumor formation, such as up-regulated miR-708-5p, associated with a decrease in pro-apoptotic caspase-2 and leukemia inhibitory factor LIF, or down-regulated miR-409-5p, associated with an increase in the p53-inhibitor MDM2. Transfection of miRNA mimics provided further evidence that changes in miRNAs are driving some observed mRNA changes. Using filtered datasets, we also identified several canonical pathways that were significantly enriched in differentially expressed miRNA-mRNA pairs, such as the epithelial-to-mesenchymal transition and the interleukin-8 signaling pathways. Overall, our data provide a more detailed understanding of KSHV latency and guide further studies of the biological significance of these changes.

We performed array comparative genome hybridization (aCGH) analyses of five Arabidopsis thaliana mutants with genomic deletions ranging in size from 4 bp to >5 kb. We used the Roche NimbleGen Arabidopsis CGH 3 × 720 K whole genome custom tiling array to optimize deletion detection. Details of the microarray design and hybridization data have been deposited at the NCBI GEO repository with accession number GSE55327.

The Transforming Growth Factor-β (TGF-β) signaling pathway is one of the major pathways essential for normal embryonic development and tissue homeostasis, with anti-tumor but also pro-metastatic properties in cancer. This pathway directly regulates several target genes that mediate its downstream functions, however very few microRNAs (miRNAs) have been identified as targets. miRNAs are modulators of gene expression with essential roles in development and a clear association with diseases including cancer. Little is known about the transcriptional regulation of the primary transcripts (pri-miRNA, pri-miR) from which several mature miRNAs are often derived. Here we present the identification of miRNAs regulated by TGF-β signaling in mouse embryonic stem (ES) cells and early embryos. We used an inducible ES cell system to maintain high levels of the TGF-β activated/phosphorylated Smad2/3 effectors, which are the transcription factors of the pathway, and a specific inhibitor that blocks their activation. By performing short RNA deep-sequencing after 12 hours Smad2/3 activation and after 16 hours inhibition, we generated a database of responsive miRNAs. Promoter/enhancer analysis of a subset of these miRNAs revealed that the transcription of pri-miR-181c/d and the pri-miR-341∼3072 cluster were found to depend on activated Smad2/3. Several of these miRNAs are expressed in early mouse embryos, when the pathway is known to play an essential role. Treatment of embryos with TGF-β inhibitor caused a reduction of their levels confirming that they are targets of this pathway in vivo. Furthermore, we showed that pri-miR-341∼3072 transcription also depends on FoxH1, a known Smad2/3 transcription partner during early development. Together, our data show that miRNAs are regulated directly by the TGF-β/Smad2/3 pathway in ES cells and early embryos. As somatic abnormalities in functions known to be regulated by the TGF-β/Smad2/3 pathway underlie tumor suppression and metastasis, this research also provides a resource for miRNAs involved in cancer.

β-III spectrin is present in the brain and is known to be important in the function of the cerebellum. Heterozygous mutations in SPTBN2, the gene encoding β-III spectrin, cause Spinocerebellar Ataxia Type 5 (SCA5), an adult-onset, slowly progressive, autosomal-dominant pure cerebellar ataxia. SCA5 is sometimes known as "Lincoln ataxia," because the largest known family is descended from relatives of the United States President Abraham Lincoln. Using targeted capture and next-generation sequencing, we identified a homozygous stop codon in SPTBN2 in a consanguineous family in which childhood developmental ataxia co-segregates with cognitive impairment. The cognitive impairment could result from mutations in a second gene, but further analysis using whole-genome sequencing combined with SNP array analysis did not reveal any evidence of other mutations. We also examined a mouse knockout of β-III spectrin in which ataxia and progressive degeneration of cerebellar Purkinje cells has been previously reported and found morphological abnormalities in neurons from prefrontal cortex and deficits in object recognition tasks, consistent with the human cognitive phenotype. These data provide the first evidence that β-III spectrin plays an important role in cortical brain development and cognition, in addition to its function in the cerebellum; and we conclude that cognitive impairment is an integral part of this novel recessive ataxic syndrome, Spectrin-associated Autosomal Recessive Cerebellar Ataxia type 1 (SPARCA1). In addition, the identification of SPARCA1 and normal heterozygous carriers of the stop codon in SPTBN2 provides insights into the mechanism of molecular dominance in SCA5 and demonstrates that the cell-specific repertoire of spectrin subunits underlies a novel group of disorders, the neuronal spectrinopathies, which includes SCA5, SPARCA1, and a form of West syndrome.

Forkhead-box protein P2 is a transcription factor that has been associated with intriguing aspects of cognitive function in humans, non-human mammals, and song-learning birds. Heterozygous mutations of the human FOXP2 gene cause a monogenic speech and language disorder. Reduced functional dosage of the mouse version (Foxp2) causes deficient cortico-striatal synaptic plasticity and impairs motor-skill learning. Moreover, the songbird orthologue appears critically important for vocal learning. Across diverse vertebrate species, this well-conserved transcription factor is highly expressed in the developing and adult central nervous system. Very little is known about the mechanisms regulated by Foxp2 during brain development. We used an integrated functional genomics strategy to robustly define Foxp2-dependent pathways, both direct and indirect targets, in the embryonic brain. Specifically, we performed genome-wide in vivo ChIP-chip screens for Foxp2-binding and thereby identified a set of 264 high-confidence neural targets under strict, empirically derived significance thresholds. The findings, coupled to expression profiling and in situ hybridization of brain tissue from wild-type and mutant mouse embryos, strongly highlighted gene networks linked to neurite development. We followed up our genomics data with functional experiments, showing that Foxp2 impacts on neurite outgrowth in primary neurons and in neuronal cell models. Our data indicate that Foxp2 modulates neuronal network formation, by directly and indirectly regulating mRNAs involved in the development and plasticity of neuronal connections.

The basic unit of genome packaging is the nucleosome, and nucleosomes have long been proposed to restrict DNA accessibility both to damage and to transcription. Nucleosome number in cells was considered fixed, but recently aging yeast and mammalian cells were shown to contain fewer nucleosomes. We show here that mammalian cells lacking High Mobility Group Box 1 protein (HMGB1) contain a reduced amount of core, linker, and variant histones, and a correspondingly reduced number of nucleosomes, possibly because HMGB1 facilitates nucleosome assembly. Yeast nhp6 mutants lacking Nhp6a and -b proteins, which are related to HMGB1, also have a reduced amount of histones and fewer nucleosomes. Nucleosome limitation in both mammalian and yeast cells increases the sensitivity of DNA to damage, increases transcription globally, and affects the relative expression of about 10% of genes. In yeast nhp6 cells the loss of more than one nucleosome in four does not affect the location of nucleosomes and their spacing, but nucleosomal occupancy. The decrease in nucleosomal occupancy is non-uniform and can be modelled assuming that different nucleosomal sites compete for available histones. Sites with a high propensity to occupation are almost always packaged into nucleosomes both in wild type and nucleosome-depleted cells; nucleosomes on sites with low propensity to occupation are disproportionately lost in nucleosome-depleted cells. We suggest that variation in nucleosome number, by affecting nucleosomal occupancy both genomewide and gene-specifically, constitutes a novel layer of epigenetic regulation.

The integrated analysis of genotypic and expression data for association with complex traits could identify novel genetic pathways involved in complex traits. We profiled 19,573 expression probes in Epstein-Barr virus-transformed lymphoblastoid cell lines (LCLs) from 299 twins and correlated these with 44 quantitative traits (QTs). For 939 expressed probes correlating with more than one QT, we investigated the presence of eQTL associations in three datasets of 57 CEU HapMap founders and 86 unrelated twins. Genome-wide association analysis of these probes with 2.2 m SNPs revealed 131 potential eQTLs (1,989 eQTL SNPs) overlapping between the HapMap datasets, five of which were in cis (58 eQTL SNPs). We then tested 535 SNPs tagging the eQTL SNPs, for association with the relevant QT in 2,905 twins. We identified nine potential SNP-QT associations (P