Background

Estimating and accounting for hidden variables is widely practiced as an important step in molecular quantitative trait locus (molecular QTL, henceforth "QTL") analysis for improving the power of QTL identification. However, few benchmark studies have been performed to evaluate the efficacy of the various methods developed for this purpose.

Results

Here we benchmark popular hidden variable inference methods including surrogate variable analysis (SVA), probabilistic estimation of expression residuals (PEER), and hidden covariates with prior (HCP) against principal component analysis (PCA)-a well-established dimension reduction and factor discovery method-via 362 synthetic and 110 real data sets. We show that PCA not only underlies the statistical methodology behind the popular methods but is also orders of magnitude faster, better-performing, and much easier to interpret and use.

Conclusions

To help researchers use PCA in their QTL analysis, we provide an R package PCAForQTL along with a detailed guide, both of which are freely available at https://github.com/heatherjzhou/PCAForQTL . We believe that using PCA rather than SVA, PEER, or HCP will substantially improve and simplify hidden variable inference in QTL mapping as well as increase the transparency and reproducibility of QTL research.

When identifying differentially expressed genes between two conditions using human population RNA-seq samples, we found a phenomenon by permutation analysis: two popular bioinformatics methods, DESeq2 and edgeR, have unexpectedly high false discovery rates. Expanding the analysis to limma-voom, NOISeq, dearseq, and Wilcoxon rank-sum test, we found that FDR control is often failed except for the Wilcoxon rank-sum test. Particularly, the actual FDRs of DESeq2 and edgeR sometimes exceed 20% when the target FDR is 5%. Based on these results, for population-level RNA-seq studies with large sample sizes, we recommend the Wilcoxon rank-sum test.

In this response to the correspondence by Hejblum et al. [1], we clarify the reasons why we ran the Wilcoxon rank-sum test on the semi-synthetic RNA-seq samples without normalization, and why we could only run dearseq with its built-in normalization, in our published study [2]. We also argue that no normalization should be performed on the semi-synthetic samples. Hence, for a fairer method comparison and using the updated dearseq package by Hejblum et al., we re-run the six differential expression methods (DESeq2, edgeR, limma-voom, dearseq, NOISeq, and the Wilcoxon rank-sum test) without normalizing the semi-synthetic samples, i.e., under the “No normalization” scheme in [1]. Our updated results show that the Wilcoxon rank-sum test is still the best method in terms of false discovery rate (FDR) control and power performance under all settings investigated.

The efficacy and safety of combining open-configuration magnetic resonance imaging (MRI) and an optical surgical navigation system (OSNS) for guidance of percutaneous liver cryoablation were evaluated in a pig model. This approach was successfully executed in eight test pigs under general anesthesia. Another two pigs spared from ablation served as controls. MRI studies, blood testing, and pathologic examinations were performed at various time points. MRI-determined minimal and maximal cryolesion diameters were 19.4 ± 1.5 mm and 44.3 ± 1.8 mm, respectively (mean volume, 8.32 ± 2.41 cm³). The lesions were readily visible by MRI on postablation Days 2 and 7, all presenting as teardrop-shaped regions of low signal intensity in T1-weighted images and high signal intensity in T2-weighted images. Mean serum alanine aminotransferase concentration transiently increased postablation, and mean blood platelet count transiently declined (both, P < 0.05). The mean necrotic tissue volume generated by pathologic assessment was less than that determined by MRI (4.24 ± 2.3 cm³ vs. 8.32 ± 2.41 cm³, P < 0.05). Cell necrosis isotherms in MRI views were within 6.0 ± 0.8 mm of cryolesion edges. Inflammatory cell infiltrates and exuberant granulation tissue surrounded necrotic areas on postablation Day 7. Combined open-configuration MRI and OSNS guidance of percutaneous liver cryoablation proved effective and safe in our animal model. This concept seems ideal for precision image-guided hepatic ablation in a clinical setting.

Purpose

Modern molecular genetics has revolutionized gene discovery, genetic diagnoses, and precision medicine yet many patients remain unable to benefit from these advances as disease-causing variants remain elusive for up to half of Mendelian genetic disorders. Patient-derived induced pluripotent stem (iPS) cells and transcriptomics were used to identify the fate of unsolved ABCA4 alleles in patients with Stargardt disease.

Methods

Multiple independent iPS lines were generated from skin biopsies of three patients with Stargardt disease harboring a single identified pathogenic ABCA4 variant. Derived retinal pigment epithelial cells (dRPE) from a normal control and patient cells were subjected to RNA-Seq on the Novaseq6000 platform, analyzed using DESeq2 with calculation of allele specific imbalance from the pathogenic or a known linked variant. Protein analysis was performed using the automated Simple Western system.

Results

Nine dRPE samples were generated, with transcriptome analysis on eight. Allele-specific expression indicated normal transcripts expressed from splice variants albeit at low levels, and missense transcripts expressed at near-normal levels. Corresponding protein was not easily detected. Patient phenotype correlation indicated missense variants expressed at high levels have more deleterious outcomes. Transcriptome analysis suggests mitochondrial membrane biodynamics and the unfolded protein response pathway may be relevant in Stargardt disease.

Conclusions

Patient-specific iPS-derived RPE cells set the stage to assess non-expressing variants in difficult-to-detect genomic regions using easily biopsied tissue.

Translational relevance

This "Disease in a Dish" approach is likely to enhance the ability of patients to participate in and benefit from clinical trials while providing insights into perturbations in RPE biology.

DNA methylation regulates cell type-specific gene expression. Here, in a transgenic mouse model, we show that deletion of the gene encoding DNA methyltransferase Dnmt3a in hypothalamic AgRP neurons causes a sedentary phenotype characterized by reduced voluntary exercise and increased adiposity. Whole-genome bisulfite sequencing (WGBS) and transcriptional profiling in neuronal nuclei from the arcuate nucleus of the hypothalamus (ARH) reveal differentially methylated genomic regions and reduced expression of AgRP neuron-associated genes in knockout mice. We use read-level analysis of WGBS data to infer putative ARH neural cell types affected by the knockout, and to localize promoter hypomethylation and increased expression of the growth factor Bmp7 to AgRP neurons, suggesting a role for aberrant TGF-β signaling in the development of this phenotype. Together, these data demonstrate that DNA methylation in AgRP neurons is required for their normal epigenetic development and neuron-specific gene expression profiles, and regulates voluntary exercise behavior.

Despite increasing success in determining genetic diagnosis for patients with inherited retinal diseases (IRDs), mutations in about 30% of the IRD cases remain unclear or unsettled after targeted gene panel or whole exome sequencing. In this study, we aimed to investigate the contributions of structural variants (SVs) to settling the molecular diagnosis of IRD with whole-genome sequencing (WGS). A cohort of 755 IRD patients whose pathogenic mutations remain undefined were subjected to WGS. Four SV calling algorithms including include MANTA, DELLY, LUMPY and CNVnator were used to detect SVs throughout the genome. All SVs identified by any one of these four algorithms were included for further analysis. AnnotSV was used to annotate these SVs. SVs that overlap with known IRD-associated genes were examined with sequencing coverage, junction reads and discordant read pairs. Polymerase Chain Reaction (PCR) followed by Sanger sequencing was used to further confirm the SVs and identify the breakpoints. Segregation of the candidate pathogenic alleles with the disease was performed when possible. A total of 16 candidate pathogenic SVs were identified in 16 families, including deletions and inversions, representing 2.1% of patients with previously unsolved IRDs. Autosomal dominant, autosomal recessive and X-linked inheritance of disease-causing SVs were observed in 12 different genes. Among these, SVs in CLN3, EYS and PRPF31 were found in multiple families. Our study suggests that the contribution of SVs detected by short-read WGS is about 0.25% of our IRD patient cohort and is significantly lower than that of single nucleotide changes and small insertions and deletions.

Background

Genetic variants can modulate phenotypic outcomes via epigenetic intermediates, for example at methylation quantitative trait loci (mQTL). We present the first large-scale assessment of mQTL at human genomic regions selected for interindividual variation in CpG methylation, which we call correlated regions of systemic interindividual variation (CoRSIVs). These can be assayed in blood DNA and do not reflect interindividual variation in cellular composition.

Results

We use target-capture bisulfite sequencing to assess DNA methylation at 4086 CoRSIVs in multiple tissues from each of 188 donors in the NIH Gene-Tissue Expression (GTEx) program. At CoRSIVs, DNA methylation in peripheral blood correlates with methylation and gene expression in internal organs. We also discover unprecedented mQTL at these regions. Genetic influences on CoRSIV methylation are extremely strong (median R²=0.76), cumulatively comprising over 70-fold more human mQTL than detected in the most powerful previous study. Moreover, mQTL beta coefficients at CoRSIVs are highly skewed (i.e., the major allele predicts higher methylation). Both surprising findings are independently validated in a cohort of 47 non-GTEx individuals. Genomic regions flanking CoRSIVs show long-range enrichments for LINE-1 and LTR transposable elements; the skewed beta coefficients may therefore reflect evolutionary selection of genetic variants that promote their methylation and silencing. Analyses of GWAS summary statistics show that mQTL polymorphisms at CoRSIVs are associated with metabolic and other classes of disease.

Conclusions

A focus on systemic interindividual epigenetic variants, clearly enhanced in mQTL content, should likewise benefit studies attempting to link human epigenetic variation to the risk of disease.

This report describes the clinical and histological findings, genetic study, and treatment in a 1.3-year-old rhesus macaque with bilateral cataracts and unilateral secondary glaucoma. Intravitreal injection of gentamicin decreased the intraocular pressure from 56 to <2 mm Hg. A putative genetic cause of the cataracts was not identified.

BACKGROUND: Methylation of the p16 promoter resulting in epigenetic gene silencing-known as p16 epimutation-is frequently found in human colorectal cancer and is also common in normal-appearing colonic mucosa of aging individuals. Thus, to improve clinical care of colorectal cancer (CRC) patients, we explored the role of age-related p16 epimutation in intestinal tumorigenesis. METHODS: We established a mouse model that replicates two common genetic and epigenetic events observed in human CRCs: Apc mutation and p16 epimutation. We conducted long-term survival and histological analysis of tumor development and progression. Colonic epithelial cells and tumors were collected from mice and analyzed by RNA sequencing (RNA-seq), quantitative PCR, and flow cytometry. We performed single-cell RNA sequencing (scRNA-seq) to characterize tumor-infiltrating immune cells throughout tumor progression. We tested whether anti-PD-L1 immunotherapy affects overall survival of tumor-bearing mice and whether inhibition of both epigenetic regulation and immune checkpoint is more efficacious. RESULTS: Mice carrying combined Apc mutation and p16 epimutation had significantly shortened survival and increased tumor growth compared to those with Apc mutation only. Intriguingly, colon tumors with p16 epimutation exhibited an activated interferon pathway, increased expression of programmed death-ligand 1 (Pdl1), and enhanced infiltration of immune cells. scRNA-seq further revealed the presence of Foxp3+ Tregs and γδT17 cells, which contribute to an immunosuppressive tumor microenvironment (TME). Furthermore, we showed that a combined therapy using an inhibitor of DNA methylation and a PD-L1 immune checkpoint inhibitor is more effective for improving survival in tumor-bearing mice than blockade of either pathway alone. CONCLUSIONS: Our study demonstrated that age-dependent p16 epimutation creates a permissive microenvironment for malignant transformation of polyps to colon cancer. Our findings provide a mechanistic rationale for future targeted therapy in patients with p16 epimutation.