Differential levels of circulating RNAs prior to endometrial cancer diagnosis.

Endometrial cancer (EC) is one of the most common female cancers and there is currently no routine screening strategy for early detection. An altered abundance of circulating microRNAs (miRNAs) and other RNA classes have the potential as early cancer biomarkers. We analyzed circulating RNA levels using small RNA sequencing, targeting RNAs in the size range of 17-47 nucleotides, in EC patients with samples collected prior to diagnosis compared to cancer-free controls. The analysis included 316 cases with samples collected 1-11 years prior to EC diagnosis, and 316 matched controls, both from the Janus Serum Bank cohort in Norway. We identified differentially abundant (DA) miRNAs, isomiRs, and small nuclear RNAs between EC cases and controls. The top EC DA miRNAs were miR-155-5p, miR-200b-3p, miR-589-5p, miR-151a-5p, miR-543, miR-485-5p, miR-625-p, and miR-671-3p. miR-200b-3p was previously reported to be among one of the top miRNAs with higher abundance in EC cases. We observed 47, 41, and 32 DA miRNAs for EC interacting with BMI, smoking status, and physical activity, respectively, including two miRNAs (miR-223-3p and miR-29b-3p) interacting with all three factors. The circulating RNAs are altered and show temporal dynamics prior to EC diagnosis. Notably, DA miRNAs for EC had the lowest q-value 4.39-6.66 years before diagnosis. Enrichment analysis of miRNAs showed that signaling pathways Fc epsilon RI, prolactin, toll-like receptor, and VEGF had the strongest associations.

The X-factor in ART: does the use of assisted reproductive technologies influence DNA methylation on the X chromosome?

BACKGROUND: Assisted reproductive technologies (ART) may perturb DNA methylation (DNAm) in early embryonic development. Although a handful of epigenome-wide association studies of ART have been published, none have investigated CpGs on the X chromosome. To bridge this knowledge gap, we leveraged one of the largest collections of mother-father-newborn trios of ART and non-ART (natural) conceptions to date to investigate sex-specific DNAm differences on the X chromosome. The discovery cohort consisted of 982 ART and 963 non-ART trios from the Norwegian Mother, Father, and Child Cohort Study (MoBa). To verify our results from the MoBa cohort, we used an external cohort of 149 ART and 58 non-ART neonates from the Australian 'Clinical review of the Health of adults conceived following Assisted Reproductive Technologies' (CHART) study. The Illumina EPIC array was used to measure DNAm in both datasets. In the MoBa cohort, we performed a set of X-chromosome-wide association studies ('XWASs' hereafter) to search for sex-specific DNAm differences between ART and non-ART newborns. We tested several models to investigate the influence of various confounders, including parental DNAm. We also searched for differentially methylated regions (DMRs) and regions of co-methylation flanking the most significant CpGs. Additionally, we ran an analogous model to our main model on the external CHART dataset. RESULTS: In the MoBa cohort, we found more differentially methylated CpGs and DMRs in girls than boys. Most of the associations persisted after controlling for parental DNAm and other confounders. Many of the significant CpGs and DMRs were in gene-promoter regions, and several of the genes linked to these CpGs are expressed in tissues relevant for both ART and sex (testis, placenta, and fallopian tube). We found no support for parental DNAm-dependent features as an explanation for the observed associations in the newborns. The most significant CpG in the boys-only analysis was in UBE2DNL, which is expressed in testes but with unknown function. The most significant CpGs in the girls-only analysis were in EIF2S3 and AMOT. These three loci also displayed differential DNAm in the CHART cohort. CONCLUSIONS: Genes that co-localized with the significant CpGs and DMRs associated with ART are implicated in several key biological processes (e.g., neurodevelopment) and disorders (e.g., intellectual disability and autism). These connections are particularly compelling in light of previous findings indicating that neurodevelopmental outcomes differ in ART-conceived children compared to those naturally conceived.

Modeling dependency structures in 450k DNA methylation data.

MOTIVATION: DNA methylation has been shown to be spatially dependent across chromosomes. Previous studies have focused on the influence of genomic context on the dependency structure, while not considering differences in dependency structure between individuals. RESULTS: We modeled spatial dependency with a flexible framework to quantify the dependency structure, focusing on inter-individual differences by exploring the association between dependency parameters and technical and biological variables. The model was applied to a subset of the Finnish Twin Cohort study (N = 1611 individuals). The estimates of the dependency parameters varied considerably across individuals, but were generally consistent across chromosomes within individuals. The variation in dependency parameters was associated with bisulfite conversion plate, zygosity, sex and age. The age differences presumably reflect accumulated environmental exposures and/or accumulated small methylation differences caused by stochastic mitotic events, establishing recognizable, individual patterns more strongly seen in older individuals. AVAILABILITY AND IMPLEMENTATION: The twin dataset used in the current study are located in the Biobank of the National Institute for Health and Welfare, Finland. All the biobanked data are publicly available for use by qualified researchers following a standardized application procedure (https://thl.fi/en/web/thl-biobank/for-researchers). A R-script for fitting the dependency structure to publicly available DNA methylation data with the software used in this article is provided in supplementary data. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

The functional diversity of tissue-resident natural killer cells against infection.

For decades, studies of natural killer (NK) cells have focused on those found in peripheral blood (PBNK cells) as the prototype for NK cell biology. Only recently have researchers begun to explore the diversity of tissue-resident NK (tr-NK) cells. While tr-NK cells were initially identified from mice parabiosis and intravascular staining experiments, they can also be identified by tissue retention markers such as CD69, CD103 and others. More importantly, tr-NK cells have distinct functions compared to PBNK cells. Within the liver, there are diverse subsets of tr-NK cells expressing different combinations of tissue-retention markers and transcription factors, the clinical relevance of which are still unclear. Functionally, liver tr-NK are primed with immediate responsiveness to infection and equipped with regulatory mechanisms to prevent liver damage. When decidual NK (dNK) cells were first discovered, they were mainly characterized by their reduced cytotoxicity and functions related to placental development. Recent studies, however, revealed different mechanisms by which dNK cells prevent uterine infections. The lungs are one of the most highly exposed sites for infection due to their role in oxygen exchange. Upon influenza infection, lung tr-NK cells can degranulate and produce more inflammatory cytokines than PBNK cells. Less understood are gut tr-NK cells which were recently characterized in infants and adults for their functional differences. In this mini-review, we aim to provide a brief overview of the most recent discoveries on how several tr-NK cells are implicated in the immune response against infection.

The CRCbiome study: a large prospective cohort study examining the role of lifestyle and the gut microbiome in colorectal cancer screening participants.

BACKGROUND: Colorectal cancer (CRC) screening reduces CRC incidence and mortality. However, current screening methods are either hampered by invasiveness or suboptimal performance, limiting their effectiveness as primary screening methods. To aid in the development of a non-invasive screening test with improved sensitivity and specificity, we have initiated a prospective biomarker study (CRCbiome), nested within a large randomized CRC screening trial in Norway. We aim to develop a microbiome-based classification algorithm to identify advanced colorectal lesions in screening participants testing positive for an immunochemical fecal occult blood test (FIT). We will also examine interactions with host factors, diet, lifestyle and prescription drugs. The prospective nature of the study also enables the analysis of changes in the gut microbiome following the removal of precancerous lesions. METHODS: The CRCbiome study recruits participants enrolled in the Bowel Cancer Screening in Norway (BCSN) study, a randomized trial initiated in 2012 comparing once-only sigmoidoscopy to repeated biennial FIT, where women and men aged 50-74 years at study entry are invited to participate. Since 2017, participants randomized to FIT screening with a positive test result have been invited to join the CRCbiome study. Self-reported diet, lifestyle and demographic data are collected prior to colonoscopy after the positive FIT-test (baseline). Screening data, including colonoscopy findings are obtained from the BCSN database. Fecal samples for gut microbiome analyses are collected both before and 2 and 12 months after colonoscopy. Samples are analyzed using metagenome sequencing, with taxonomy profiles, and gene and pathway content as primary measures. CRCbiome data will also be linked to national registries to obtain information on prescription histories and cancer relevant outcomes occurring during the 10 year follow-up period. DISCUSSION: The CRCbiome study will increase our understanding of how the gut microbiome, in combination with lifestyle and environmental factors, influences the early stages of colorectal carcinogenesis. This knowledge will be crucial to develop microbiome-based screening tools for CRC. By evaluating biomarker performance in a screening setting, using samples from the target population, the generalizability of the findings to future screening cohorts is likely to be high. TRIAL REGISTRATION: ClinicalTrials.gov Identifier: NCT01538550 .

Blood-based epigenetic estimators of chronological age in human adults using DNA methylation data from the Illumina MethylationEPIC array.

BACKGROUND: Epigenetic clocks have been recognized for their precise prediction of chronological age, age-related diseases, and all-cause mortality. Existing epigenetic clocks are based on CpGs from the Illumina HumanMethylation450 BeadChip (450 K) which has now been replaced by the latest platform, Illumina MethylationEPIC BeadChip (EPIC). Thus, it remains unclear to what extent EPIC contributes to increased precision and accuracy in the prediction of chronological age. RESULTS: We developed three blood-based epigenetic clocks for human adults using EPIC-based DNA methylation (DNAm) data from the Norwegian Mother, Father and Child Cohort Study (MoBa) and the Gene Expression Omnibus (GEO) public repository: 1) an Adult Blood-based EPIC Clock (ABEC) trained on DNAm data from MoBa (n = 1592, age-span: 19 to 59 years), 2) an extended ABEC (eABEC) trained on DNAm data from MoBa and GEO (n = 2227, age-span: 18 to 88 years), and 3) a common ABEC (cABEC) trained on the same training set as eABEC but restricted to CpGs common to 450 K and EPIC. Our clocks showed high precision (Pearson correlation between chronological and epigenetic age (r) > 0.94) in independent cohorts, including GSE111165 (n = 15), GSE115278 (n = 108), GSE132203 (n = 795), and the Epigenetics in Pregnancy (EPIPREG) study of the STORK Groruddalen Cohort (n = 470). This high precision is unlikely due to the use of EPIC, but rather due to the large sample size of the training set. CONCLUSIONS: Our ABECs predicted adults' chronological age precisely in independent cohorts. As EPIC is now the dominant platform for measuring DNAm, these clocks will be useful in further predictions of chronological age, age-related diseases, and mortality.

DNA methylation patterns vary in boar sperm cells with different levels of DNA fragmentation.

BACKGROUND: Sperm DNA integrity is considered essential for successful transmission of the paternal genome, fertilization and normal embryo development. DNA fragmentation index (DFI, %) has become a key parameter in the swine artificial insemination industry to assess sperm DNA integrity. Recently, in some elite Norwegian Landrace boars (boars with excellent field fertility records), a higher level of sperm DFI has been observed. In order to obtain a better understanding of this, and to study the complexity of sperm DNA integrity, liquid preserved semen samples from elite boars with contrasting DFI levels were examined for protamine deficiency, thiol profile and disulphide bonds. Additionally, the DNA methylation profiles of the samples were determined by reduced representation bisulphite sequencing (RRBS). RESULTS: In this study, different traits related to sperm DNA integrity were investigated (n = 18 ejaculates). Upon liquid storage, the levels of total thiols and disulphide bonds decreased significantly, while the DFI and protamine deficiency level increased significantly. The RRBS results revealed similar global patterns of low methylation from semen samples with different levels of DFI (low, medium and high). Differential methylation analyses indicated that the number of differentially methylated cytosines (DMCs) increased in the low-high compared to the low-medium and the medium-high DFI groups. Annotating the DMCs with gene and CpG features revealed clear differences between DFI groups. In addition, the number of annotated transcription starting sites (TSS) and associated pathways in the low-high comparison was greater than the other two groups. Pathway analysis showed that genes (based on the closest TSS to DMCs) corresponding to low-high DFI comparison were associated with important processes such as membrane function, metabolic cascade and antioxidant defence system. CONCLUSION: To our knowledge, this is the first study evaluating DNA methylation in boar sperm cells with different levels of DFI. The present study shows that sperm cells with varying levels of DNA fragmentation exhibit similar global methylation, but different site-specific DNA methylation signatures. Moreover, with increasing DNA fragmentation in spermatozoa, there is an increase in the number of potentially affected downstream genes and their respective regulatory pathways.

Long-Term Neurodevelopment of Low-Birthweight, Preterm Infants with Patent Ductus Arteriosus.

OBJECTIVE: To evaluate whether the presence of patent ductus arteriosus (PDA) in preterm infants worsens long-term neurodevelopmental outcomes. STUDY DESIGN: This was a secondary observational analysis of data from 1090 preterm low-birthweight infants in the Infant Health and Development Program (IHDP), a multicenter longitudinal cohort study of outcomes assessed from 3 to 18 years of age. Multivariable analysis was adjusted for IHDP treatment group (intervention or follow-up), birth weight, maternal race, maternal education, infant sex, maternal preconception weight, Home Observation Measurement of the Environment (HOME) total score at 12 months, neonatal health index, and gestational age. RESULTS: Of the 1090 patients (49% male) included in the analysis, 135 had a PDA. Mean birth weight (1322 g vs 1871 g; P < .0001) and gestational age (30.2 weeks vs 33.4 weeks, P < .0001) were lower and mean ventilator days (11.8 vs 1.3; P < .0001), vasopressor use (12.6% vs 1.2%; P < .0001), and congestive heart failure (8.9% vs 0.1%; P < .0001) were higher in the PDA group. There were no differences between the PDA and no-PDA groups in maternal education level and HOME total score at age 12 months. Multivariable analysis demonstrated no between-group differences in cognitive development or behavioral competence at age 3, 8, and 18 years. CONCLUSIONS: The presence of a PDA in moderately preterm, low-birthweight infants does not impact long-term neurodevelopmental outcomes.

A comprehensive profile of circulating RNAs in human serum.

Non-coding RNA (ncRNA) molecules have fundamental roles in cells and many are also stable in body fluids as extracellular RNAs. In this study, we used RNA sequencing (RNA-seq) to investigate the profile of small non-coding RNA (sncRNA) in human serum. We analyzed 10 billion Illumina reads from 477 serum samples, included in the Norwegian population-based Janus Serum Bank (JSB). We found that the core serum RNA repertoire includes 258 micro RNAs (miRNA), 441 piwi-interacting RNAs (piRNA), 411 transfer RNAs (tRNA), 24 small nucleolar RNAs (snoRNA), 125 small nuclear RNAs (snRNA) and 123 miscellaneous RNAs (misc-RNA). We also investigated biological and technical variation in expression, and the results suggest that many RNA molecules identified in serum contain signs of biological variation. They are therefore unlikely to be random degradation by-products. In addition, the presence of specific fragments of tRNA, snoRNA, Vault RNA and Y_RNA indicates protection from degradation. Our results suggest that many circulating RNAs in serum can be potential biomarkers.

cnvScan: a CNV screening and annotation tool to improve the clinical utility of computational CNV prediction from exome sequencing data.

BACKGROUND: With advances in next generation sequencing technology and analysis methods, single nucleotide variants (SNVs) and indels can be detected with high sensitivity and specificity in exome sequencing data. Recent studies have demonstrated the ability to detect disease-causing copy number variants (CNVs) in exome sequencing data. However, exonic CNV prediction programs have shown high false positive CNV counts, which is the major limiting factor for the applicability of these programs in clinical studies. RESULTS: We have developed a tool (cnvScan) to improve the clinical utility of computational CNV prediction in exome data. cnvScan can accept input from any CNV prediction program. cnvScan consists of two steps: CNV screening and CNV annotation. CNV screening evaluates CNV prediction using quality scores and refines this using an in-house CNV database, which greatly reduces the false positive rate. The annotation step provides functionally and clinically relevant information using multiple source datasets. We assessed the performance of cnvScan on CNV predictions from five different prediction programs using 64 exomes from Primary Immunodeficiency (PIDD) patients, and identified PIDD-causing CNVs in three individuals from two different families. CONCLUSIONS: In summary, cnvScan reduces the time and effort required to detect disease-causing CNVs by reducing the false positive count and providing annotation. This improves the clinical utility of CNV detection in exome data.

A comparative study of ChIP-seq sequencing library preparation methods.

BACKGROUND: ChIP-seq is the primary technique used to investigate genome-wide protein-DNA interactions. As part of this procedure, immunoprecipitated DNA must undergo "library preparation" to enable subsequent high-throughput sequencing. To facilitate the analysis of biopsy samples and rare cell populations, there has been a recent proliferation of methods allowing sequencing library preparation from low-input DNA amounts. However, little information exists on the relative merits, performance, comparability and biases inherent to these procedures. Notably, recently developed single-cell ChIP procedures employing microfluidics must also employ library preparation reagents to allow downstream sequencing. RESULTS: In this study, seven methods designed for low-input DNA/ChIP-seq sample preparation (Accel-NGS® 2S, Bowman-method, HTML-PCR, SeqPlex™, DNA SMART™, TELP and ThruPLEX®) were performed on five replicates of 1 ng and 0.1 ng input H3K4me3 ChIP material, and compared to a "gold standard" reference PCR-free dataset. The performance of each method was examined for the prevalence of unmappable reads, amplification-derived duplicate reads, reproducibility, and for the sensitivity and specificity of peak calling. CONCLUSIONS: We identified consistent high performance in a subset of the tested reagents, which should aid researchers in choosing the most appropriate reagents for their studies. Furthermore, we expect this work to drive future advances by identifying and encouraging use of the most promising methods and reagents. The results may also aid judgements on how comparable are existing datasets that have been prepared with different sample library preparation reagents.

Age-related somatic structural changes in the nuclear genome of human blood cells.

Structural variations are among the most frequent interindividual genetic differences in the human genome. The frequency and distribution of de novo somatic structural variants in normal cells is, however, poorly explored. Using age-stratified cohorts of 318 monozygotic (MZ) twins and 296 single-born subjects, we describe age-related accumulation of copy-number variation in the nuclear genomes in vivo and frequency changes for both megabase- and kilobase-range variants. Megabase-range aberrations were found in 3.4% (9 of 264) of subjects ≥60 years old; these subjects included 78 MZ twin pairs and 108 single-born individuals. No such findings were observed in 81 MZ pairs or 180 single-born subjects who were ≤55 years old. Recurrent region- and gene-specific mutations, mostly deletions, were observed. Longitudinal analyses of 43 subjects whose data were collected 7-19 years apart suggest considerable variation in the rate of accumulation of clones carrying structural changes. Furthermore, the longitudinal analysis of individuals with structural aberrations suggests that there is a natural self-removal of aberrant cell clones from peripheral blood. In three healthy subjects, we detected somatic aberrations characteristic of patients with myelodysplastic syndrome. The recurrent rearrangements uncovered here are candidates for common age-related defects in human blood cells. We anticipate that extension of these results will allow determination of the genetic age of different somatic-cell lineages and estimation of possible individual differences between genetic and chronological age. Our work might also help to explain the cause of an age-related reduction in the number of cell clones in the blood; such a reduction is one of the hallmarks of immunosenescence.

The Concordance and Heritability of Type 2 Diabetes in 34,166 Twin Pairs From International Twin Registers: The Discordant Twin (DISCOTWIN) Consortium.

Twin pairs discordant for disease may help elucidate the epigenetic mechanisms and causal environmental factors in disease development and progression. To obtain the numbers of pairs, especially monozygotic (MZ) twin pairs, necessary for in-depth studies while also allowing for replication, twin studies worldwide need to pool their resources. The Discordant Twin (DISCOTWIN) consortium was established for this goal. Here, we describe the DISCOTWIN Consortium and present an analysis of type 2 diabetes (T2D) data in nearly 35,000 twin pairs. Seven twin cohorts from Europe (Denmark, Finland, Norway, the Netherlands, Spain, Sweden, and the United Kingdom) and one from Australia investigated the rate of discordance for T2D in same-sex twin pairs aged 45 years and older. Data were available for 34,166 same-sex twin pairs, of which 13,970 were MZ, with T2D diagnosis based on self-reported diagnosis and medication use, fasting glucose and insulin measures, or medical records. The prevalence of T2D ranged from 2.6% to 12.3% across the cohorts depending on age, body mass index (BMI), and national diabetes prevalence. T2D discordance rate was lower for MZ (5.1%, range 2.9-11.2%) than for same-sex dizygotic (DZ) (8.0%, range 4.9-13.5%) pairs. Across DISCOTWIN, 720 discordant MZ pairs were identified. Except for the oldest of the Danish cohorts (mean age 79), heritability estimates based on contingency tables were moderate to high (0.47-0.77). From a meta-analysis of all data, the heritability was estimated at 72% (95% confidence interval 61-78%). This study demonstrated high T2D prevalence and high heritability for T2D liability across twin cohorts. Therefore, the number of discordant MZ pairs for T2D is limited. By combining national resources, the DISCOTWIN Consortium maximizes the number of discordant MZ pairs needed for in-depth genotyping, multi-omics, and phenotyping studies, which may provide unique insights into the pathways linking genes to the development of many diseases.

Transcription profiling during the cell cycle shows that a subset of Polycomb-targeted genes is upregulated during DNA replication.

Genome-wide gene expression analyses of the human somatic cell cycle have indicated that the set of cycling genes differ between primary and cancer cells. By identifying genes that have cell cycle dependent expression in HaCaT human keratinocytes and comparing these with previously identified cell cycle genes, we have identified three distinct groups of cell cycle genes. First, housekeeping genes enriched for known cell cycle functions; second, cell type-specific genes enriched for HaCaT-specific functions; and third, Polycomb-regulated genes. These Polycomb-regulated genes are specifically upregulated during DNA replication, and consistent with being epigenetically silenced in other cell cycle phases, these genes have lower expression than other cell cycle genes. We also find similar patterns in foreskin fibroblasts, indicating that replication-dependent expression of Polycomb-silenced genes is a prevalent but unrecognized regulatory mechanism.

DNA methylation and gene expression changes in monozygotic twins discordant for psoriasis: identification of epigenetically dysregulated genes.

Monozygotic (MZ) twins do not show complete concordance for many complex diseases; for example, discordance rates for autoimmune diseases are 20%-80%. MZ discordance indicates a role for epigenetic or environmental factors in disease. We used MZ twins discordant for psoriasis to search for genome-wide differences in DNA methylation and gene expression in CD4(+) and CD8(+) cells using Illumina's HumanMethylation27 and HT-12 expression assays, respectively. Analysis of these data revealed no differentially methylated or expressed genes between co-twins when analyzed separately, although we observed a substantial amount of small differences. However, combined analysis of DNA methylation and gene expression identified genes where differences in DNA methylation between unaffected and affected twins were correlated with differences in gene expression. Several of the top-ranked genes according to significance of the correlation in CD4(+) cells are known to be associated with psoriasis. Further, gene ontology (GO) analysis revealed enrichment of biological processes associated with the immune response and clustering of genes in a biological pathway comprising cytokines and chemokines. These data suggest that DNA methylation is involved in an epigenetic dysregulation of biological pathways involved in the pathogenesis of psoriasis. This is the first study based on data from MZ twins discordant for psoriasis to detect epigenetic alterations that potentially contribute to development of the disease.

Alternative splicing of placental lactogen (CSH2) in somatotroph pituitary adenomas.

OBJECTIVES: Somatotroph adenomas secrete supraphysiological amounts of GH, causing acromegaly. We have previously shown epithelial splicing regulator 1 (ESRP1) to play a role in epithelial mesenchymal transition (EMT) progression in these adenomas and account for poor treatment response. We evaluated if the mRNA levels of the GH/CSH gene cluster in somatotroph adenomas are associated with an epithelial phenotype and response to SA treatment. METHODS: We investigated the associations between ESRP1 and the growth hormone/chorionic somatomammotropin (GH/CSH) gene cluster by RNA sequencing (RNAseq). CSH2 isoform 3 mRNA was further evaluated in 65 somatotroph adenomas and associations with disease severity and treatment response. RESULTS: mRNA for all genes in the GH/CSH cluster were expressed, however, only chorionic somatomammotropin 2/placental lactogen 2 (CSH2) displayed an alternative splicing pattern. CSH2 isoform 3 was associated with a dense granulation pattern and an epithelial phenotype with high levels of ESRP1 and E-cadherin expression. Further, CSH2 isoform 3 was associated with reduced serum GH and IGF-I levels after somatostatin analog treatment. CONCLUSIONS: Attenuated CSH2 isoform 3 was associated with mesenchymal phenotype and a blunted clinical response to somatostatin analog treatment in patients with acromegaly.

A dominant STIM1 mutation causes Stormorken syndrome.

Stormorken syndrome is a rare autosomal-dominant disease with mild bleeding tendency, thrombocytopathy, thrombocytopenia, mild anemia, asplenia, tubular aggregate myopathy, miosis, headache, and ichthyosis. A heterozygous missense mutation in STIM1 exon 7 (c.910C>T; p.Arg304Trp) (NM_003156.3) was found to segregate with the disease in six Stormorken syndrome patients in four families. Upon sensing Ca(2+) depletion in the endoplasmic reticulum lumen, STIM1 undergoes a conformational change enabling it to interact with and open ORAI1, a Ca(2+) release-activated Ca(2+) channel located in the plasma membrane. The STIM1 mutation found in Stormorken syndrome patients is located in the coiled-coil 1 domain, which might play a role in keeping STIM1 inactive. In agreement with a possible gain-of-function mutation in STIM1, blood platelets from patients were in a preactivated state with high exposure of aminophospholipids on the outer surface of the plasma membrane. Resting Ca(2+) levels were elevated in platelets from the patients compared with controls, and store-operated Ca(2+) entry was markedly attenuated, further supporting constitutive activity of STIM1 and ORAI1. Thus, our data are compatible with a near-maximal activation of STIM1 in Stormorken syndrome patients. We conclude that the heterozygous mutation c.910C>T causes the complex phenotype that defines this syndrome.

Identification of copy number variants from exome sequence data.

BACKGROUND: With advances in next generation sequencing technologies and genomic capture techniques, exome sequencing has become a cost-effective approach for mutation detection in genetic diseases. However, computational prediction of copy number variants (CNVs) from exome sequence data is a challenging task. Whilst numerous programs are available, they have different sensitivities, and have low sensitivity to detect smaller CNVs (1-4 exons). Additionally, exonic CNV discovery using standard aCGH has limitations due to the low probe density over exonic regions. The goal of our study was to develop a protocol to detect exonic CNVs (including shorter CNVs that cover 1-4 exons), combining computational prediction algorithms and a high-resolution custom CGH array. RESULTS: We used six published CNV prediction programs (ExomeCNV, CONTRA, ExomeCopy, ExomeDepth, CoNIFER, XHMM) and an in-house modification to ExomeCopy and ExomeDepth (ExCopyDepth) for computational CNV prediction on 30 exomes from the 1000 genomes project and 9 exomes from primary immunodeficiency patients. CNV predictions were tested using a custom CGH array designed to capture all exons (exaCGH). After this validation, we next evaluated the computational prediction of shorter CNVs. ExomeCopy and the in-house modified algorithm, ExCopyDepth, showed the highest capability in detecting shorter CNVs. Finally, the performance of each computational program was assessed by calculating the sensitivity and false positive rate. CONCLUSIONS: In this paper, we assessed the ability of 6 computational programs to predict CNVs, focussing on short (1-4 exon) CNVs. We also tested these predictions using a custom array targeting exons. Based on these results, we propose a protocol to identify and confirm shorter exonic CNVs combining computational prediction algorithms and custom aCGH experiments.

Extensive variation and low heritability of DNA methylation identified in a twin study.

Disturbance of DNA methylation leading to aberrant gene expression has been implicated in the etiology of many diseases. Whereas variation at the genetic level has been studied extensively, less is known about the extent and function of epigenetic variation. To explore variation and heritability of DNA methylation, we performed bisulfite sequencing of 1760 CpG sites in 186 regions in the human major histocompatibility complex (MHC) in CD4+ lymphocytes from 49 monozygotic (MZ) and 40 dizygotic (DZ) twin pairs. Individuals show extensive variation in DNA methylation both between and within regions. In addition, many regions also have a complex pattern of variation. Globally, there appears to be a bimodal distribution of DNA methylation in the regions, but a significant fraction of the CpG sites are also heterogeneously methylated. Classification of regions into CpG islands (intragenic and intergenic), 5' end of genes not associated with a defined CpG island, conserved noncoding regions, and random CpG sites shows region-type differences in variation and heritability. Analyses revealed slightly lower intra-pair differences among MZ than among DZ pairs, suggesting some genetic influences on DNA methylation variation, with most of the variance attributed to nongenetic factors. Overall, heritability estimates of DNA methylation were low. Our heritability estimates are, however, somewhat deflated due to the presence of batch effects that artificially inflate the estimates of shared environment.