Telomerecat: A ploidy-agnostic method for estimating telomere length from whole genome sequencing data.

Telomere length is a risk factor in disease and the dynamics of telomere length are crucial to our understanding of cell replication and vitality. The proliferation of whole genome sequencing represents an unprecedented opportunity to glean new insights into telomere biology on a previously unimaginable scale. To this end, a number of approaches for estimating telomere length from whole-genome sequencing data have been proposed. Here we present Telomerecat, a novel approach to the estimation of telomere length. Previous methods have been dependent on the number of telomeres present in a cell being known, which may be problematic when analysing aneuploid cancer data and non-human samples. Telomerecat is designed to be agnostic to the number of telomeres present, making it suited for the purpose of estimating telomere length in cancer studies. Telomerecat also accounts for interstitial telomeric reads and presents a novel approach to dealing with sequencing errors. We show that Telomerecat performs well at telomere length estimation when compared to leading experimental and computational methods. Furthermore, we show that it detects expected patterns in longitudinal data, repeated measurements, and cross-species comparisons. We also apply the method to a cancer cell data, uncovering an interesting relationship with the underlying telomerase genotype

Genetic determinants of risk in pulmonary arterial hypertension: international genome-wide association studies and meta-analysis

Background Rare genetic variants cause pulmonary arterial hypertension, but the contribution of common genetic variation to disease risk and natural history is poorly characterised. We tested for genome-wide association for pulmonary arterial hypertension in large international cohorts and assessed the contribution of associated regions to outcomes. Methods We did two separate genome-wide association studies (GWAS) and a meta-analysis of pulmonary arterial hypertension. These GWAS used data from four international case-control studies across 11744 individuals with European ancestry (including 2085 patients). One GWAS used genotypes from 5895 whole-genome sequences and the other GWAS used genotyping array data from an additional 5849 individuals. Cross-validation of loci reaching genome-wide significance was sought by meta-analysis. Conditional analysis corrected for the most significant variants at each locus was used to resolve signals for multiple associations. We functionally annotated associated variants and tested associations with duration of survival. All-cause mortality was the primary endpoint in survival analyses. Findings A locus near SOX17 (rs10103692, odds ratio 1·80 [95% CI 1·55–2·08], p=5·13×10– ¹⁵) and a second locus in HLA-DPA1 and HLA-DPB1 (collectively referred to as HLA-DPA1/DPB1 here; rs2856830, 1·56 [1·42–1·71], p=7·65×10– ²⁰) within the class II MHC region were associated with pulmonary arterial hypertension. The SOX17 locus had two independent signals associated with pulmonary arterial hypertension (rs13266183, 1·36 [1·25–1·48], p=1·69×10– ¹²; and rs10103692). Functional and epigenomic data indicate that the risk variants near SOX17 alter gene regulation via an enhancer active in endothelial cells. Pulmonary arterial hypertension risk variants determined haplotype-specific enhancer activity, and CRISPR-mediated inhibition of the enhancer reduced SOX17 expression. The HLA-DPA1/DPB1 rs2856830 genotype was strongly associated with survival. Median survival from diagnosis in patients with pulmonary arterial hypertension with the C/C homozygous genotype was double (13·50 years [95% CI 12·07 to >13·50]) that of those with the T/T genotype (6·97 years [6·02–8·05]), despite similar baseline disease severity. Interpretation This is the first study to report that common genetic variation at loci in an enhancer near SOX17 and in HLA-DPA1/DPB1 is associated with pulmonary arterial hypertension. Impairment of SOX17 function might be more common in pulmonary arterial hypertension than suggested by rare mutations in SOX17. Further studies are needed to confirm the association between HLA typing or rs2856830 genotyping and survival, and to determine whether HLA typing or rs2856830 genotyping improves risk stratification in clinical practice or trials. Funding UK NIHR, BHF, UK MRC, Dinosaur Trust, NIH/NHLBI, ERS, EMBO, Wellcome Trust, EU, AHA, ACClinPharm, Netherlands CVRI, Dutch Heart Foundation, Dutch Federation of UMC, Netherlands OHRD and RNAS, German DFG, German BMBF, APH Paris, INSERM, Université Paris-Sud, and French ANR

Publisher Correction: Telomerecat: A ploidy-agnostic method for estimating telomere length from whole genome sequencing data.

A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has been fixed in the paper

miR-122 inhibition in a human liver organoid model leads to liver inflammation, necrosis, steatofibrosis and dysregulated insulin signaling

<div><p>To investigate the role of miR-122 in the development and regression of non-alcoholic fatty liver disease (NAFLD) <i>in vitro</i>, we used multicellular 3D human liver organoids developed in our laboratory. These organoids consist of primary human hepatocytes, Kupffer cells, quiescent stellate cells and liver sinusoidal endothelial cells. They remain viable and functional for 4 weeks expressing typical markers of liver function such as synthesis of albumin, urea, and alpha-1 p450 drug metabolism. Before mixing, hepatic cells were transduced with lentivirus to inhibit miR122 expression (ABM, CA). Immediately after the organoids were fully formed (day 4) or after 1 or 2 weeks of additional incubation (days 11 or 18), the organoids were analyzed using fluorescent live/dead staining and ATP production; total RNA was extracted for qPCR gene expression profiling. Our results show that miR-122 inhibition in liver organoids leads to inflammation, necrosis, steatosis and fibrosis. This was associated with increase in inflammatory cytokines (IL6, TNF), chemokines (CCL2, CCL3) and increase in a subset of Matrix Metaloproteinases (MMP8, MMP9). An altered expression of key genes in lipid metabolism (i.e LPL, LDLR) and insulin signaling (i.e GLUT4, IRS1) was also identified. <i>Conclusion</i>: Our results highlight the role of miR-122 inhibition in liver inflammation, steatofibrosis and dysregulation of insulin signaling. Patients with NAFLD are known to have altered levels of miR-122, therefore we suggest that miR-122 mimics could play a useful role in reversing liver steatofibrosis and insulin resistance seen in patients with NAFLD.</p></div

Sirius red staining of liver organoids.

<p>IHC sections of liver organoids stained with Sirius Red at 4 days post-infection (A-C), 11 days post-infection (D-F), and 18 days post-infection (G-I). Left column pictures (A, D, G) show control organoids, middle column show miR-122 inhibited organoids with MOI = 2 (B, E, H), right column show miR-122 inhibited organoids with MOI = 5 (C, F, I). Arrows indicate fibrosis.</p

Liver organoid live/dead cell assay.

<p>Liver organoids infected with either control lentivirus with MOI = 5 (A-C) or with lentivirus expressing miR-122 inhibitor at different weeks following infection. Organoids were fixed, with calcein Blue and Ethidium homodimer and images were taken using a Leica confocal microscope. Blue = cell nuclei, Red = dead cells, Green = lentivirus GFP expression. W denotes week.</p

H&E staining of liver organoids.

<p>IHC sections of liver organoids stained with H&E at 4 days post-infection (A-C), 11 days post-infection (D-F), and 18 days post-infection (G-I). Left column pictures (A, D, G) show control organoids, middle column show miR-122 inhibited organoids with MOI = 2 (B, E, H), right column show miR-122 inhibited organoids with MOI = 5 (C, F, I). Arrows indicate ballooned hepatocytes, arrowhead indicate Mallory’s hyaline.</p

Construction of liver organoids.

<p>(A) shows aggregation and compaction of cells at 24 hr. intervals. Organoids were fully formed after 4 days undisturbed culture. (B) Direct viability visualization of organoids over a 28-day period using calcein AM to stain live cells green and ethidium homodimer to stain dead cells red. Even by 28 days, the percentage of dead cells was <10%.</p