The UCSC Genome Browser database: 2018 update

Serial single-cell RNA sequencing unveils drug resistance and metastatic traits in stage IV breast cancer

Metastasis is a complex process that remains poorly understood at the molecular levels. We profiled single-cell transcriptomic, genomic, and epigenomic changes associated with cancer cell progression, chemotherapy resistance, and metastasis from a Stage IV breast cancer patient. Pretreatment- and posttreatment-specimens from the primary tumor and distant metastases were collected for single-cell RNA sequencing and subsequent cell clustering, copy number variation (CNV) estimation, transcriptomic factor estimation, and pseudotime analyses. CNV analysis revealed that a small population of pretreatment cancer cells resisted chemotherapy and expanded. New clones including Metastatic Precursor Cells (MPCs), emerged in the posttreatment primary tumors in CNV similar to metastatic cells. MPCs exhibited expression profiles indicative of epithelial-mesenchymal transition. Comparison of MPCs with metastatic cancer cells also revealed dynamic changes in transcription factors and calcitonin pathway gene expression. These findings demonstrate the utility of single-patient clinical sample analysis for understanding tumor drug resistance, regrowth, and metastasis.

A liver immune rheostat regulates CD8 T cell immunity in chronic HBV infection

Chronic hepatitis B virus (HBV) infection affects 300 million patients worldwide1,2, in whom virus-specific CD8 T cells by still ill-defined mechanisms lose their function and cannot eliminate HBV-infected hepatocytes3-7. Here we demonstrate that a liver immune rheostat renders virus-specific CD8 T cells refractory to activation and leads to their loss of effector functions. In preclinical models of persistent infection with hepatotropic viruses such as HBV, dysfunctional virus-specific CXCR6+ CD8 T cells accumulated in the liver and, as a characteristic hallmark, showed enhanced transcriptional activity of cAMP-responsive element modulator (CREM) distinct from T cell exhaustion. In patients with chronic hepatitis B, circulating and intrahepatic HBV-specific CXCR6+ CD8 T cells with enhanced CREM expression and transcriptional activity were detected at a frequency of 12-22% of HBV-specific CD8 T cells. Knocking out the inhibitory CREM/ICER isoform in T cells, however, failed to rescue T cell immunity. This indicates that CREM activity was a consequence, rather than the cause, of loss in T cell function, further supported by the observation of enhanced phosphorylation of protein kinase A (PKA) which is upstream of CREM. Indeed, we found that enhanced cAMP-PKA-signalling from increased T cell adenylyl cyclase activity augmented CREM activity and curbed T cell activation and effector function in persistent hepatic infection. Mechanistically, CD8 T cells recognizing their antigen on hepatocytes established close and extensive contact with liver sinusoidal endothelial cells, thereby enhancing adenylyl cyclase-cAMP-PKA signalling in T cells. In these hepatic CD8 T cells, which recognize their antigen on hepatocytes, phosphorylation of key signalling kinases of the T cell receptor signalling pathway was impaired, which rendered them refractory to activation. Thus, close contact with liver sinusoidal endothelial cells curbs the activation and effector function of HBV-specific CD8 T cells that target hepatocytes expressing viral antigens by means of the adenylyl cyclase-cAMP-PKA axis in an immune rheostat-like fashion.

Uncoupling the TFIIH Core and Kinase Modules Leads To Misregulated RNA Polymerase II CTD Serine 5 Phosphorylation

TFIIH is an essential transcription initiation factor for RNA polymerase II (RNApII). This multi-subunit complex comprises two modules that are physically linked by the subunit Tfb3 (MAT1 in metazoans). The TFIIH Core Module, with two DNA-dependent ATPases and several additional subunits, promotes DNA unwinding. The TFIIH Kinase Module phosphorylates Serine 5 of the C-terminal domain (CTD) of RNApII subunit Rpb1, a modification that coordinates exchange of initiation and early elongation factors. While it is not obvious why these two disparate activities are bundled into one factor, the connection may provide temporal coordination during early initiation. Here we show that Tfb3 can be split into two parts to uncouple the TFIIH modules. The resulting cells grow slower than normal, but are viable. Chromatin immunoprecipitation of the split TFIIH shows that the Core Module, but not the Kinase, is properly recruited to promoters. Instead of the normal promoter-proximal peak, high CTD Serine 5 phosphorylation is seen throughout transcribed regions. Therefore, coupling the TFIIH modules is necessary to localize and limit CTD kinase activity to early stages of transcription. These results are consistent with the idea that the two TFIIH modules began as independent functional entities that became connected by Tfb3 during early eukaryotic evolution.

Validated Negative Regions (VNRs) in the VISTA Database might be Truncated Forms of Bona Fide Enhancers

The VISTA enhancer database is a valuable resource for evaluating predicted enhancers in humans and mice. In addition to thousands of validated positive regions (VPRs) in the human and mouse genomes, the database also contains similar numbers of validated negative regions (VNRs). It is previously shown that the VPRs are on average half as long as predicted overlapping enhancers that are highly conserved and hypothesize that the VPRs may be truncated forms of long bona fide enhancers. Here, it is shown that like the VPRs, the VNRs also are under strong evolutionary constraints and overlap predicted enhancers in the genomes. The VNRs are also on average half as long as predicted overlapping enhancers that are highly conserved. Moreover, the VNRs and the VPRs display similar cell/tissue-specific modification patterns of key epigenetic marks of active enhancers. Furthermore, the VNRs and the VPRs show similar impact score spectra of in silico mutagenesis. These highly similar properties between the VPRs and the VNRs suggest that like the VPRs, the VNRs may also be truncated forms of long bona fide enhancers.

MFPINC: prediction of plant ncRNAs based on multi-source feature fusion

Non-coding RNAs (ncRNAs) are recognized as pivotal players in the regulation of essential physiological processes such as nutrient homeostasis, development, and stress responses in plants. Common methods for predicting ncRNAs are susceptible to significant effects of experimental conditions and computational methods, resulting in the need for significant investment of time and resources. Therefore, we constructed an ncRNA predictor(MFPINC), to predict potential ncRNA in plants which is based on the PINC tool proposed by our previous studies. Specifically, sequence features were carefully refined using variance thresholding and F-test methods, while deep features were extracted and feature fusion were performed by applying the GRU model. The comprehensive evaluation of multiple standard datasets shows that MFPINC not only achieves more comprehensive and accurate identification of gene sequences, but also significantly improves the expressive and generalization performance of the model, and MFPINC significantly outperforms the existing competing methods in ncRNA identification. In addition, it is worth mentioning that our tool can also be found on Github ( https://github.com/Zhenj-Nie/MFPINC ) the data and source code can also be downloaded for free.

The rapidly evolving X-linked MIR-506 family fine-tunes spermatogenesis to enhance sperm competition

Despite rapid evolution across eutherian mammals, the X-linked MIR-506 family miRNAs are located in a region flanked by two highly conserved protein-coding genes (SLITRK2 and FMR1) on the X chromosome. Intriguingly, these miRNAs are predominantly expressed in the testis, suggesting a potential role in spermatogenesis and male fertility. Here, we report that the X-linked MIR-506 family miRNAs were derived from the MER91C DNA transposons. Selective inactivation of individual miRNAs or clusters caused no discernible defects, but simultaneous ablation of five clusters containing 19 members of the MIR-506 family led to reduced male fertility in mice. Despite normal sperm counts, motility, and morphology, the KO sperm were less competitive than wild-type sperm when subjected to a polyandrous mating scheme. Transcriptomic and bioinformatic analyses revealed that these X-linked MIR-506 family miRNAs, in addition to targeting a set of conserved genes, have more targets that are critical for spermatogenesis and embryonic development during evolution. Our data suggest that the MIR-506 family miRNAs function to enhance sperm competitiveness and reproductive fitness of the male by finetuning gene expression during spermatogenesis.

Yeast zinc cluster transcription factors involved in the switch from fermentation to respiration show interdependency for DNA binding revealing a novel type of DNA recognition

In budding yeast, fermentation is the most important pathway for energy production. Under low-glucose conditions, ethanol is used for synthesis of this sugar requiring a shift to respiration. This process is controlled by the transcriptional regulators Cat8, Sip4, Rds2 and Ert1. We characterized Gsm1 (glucose starvation modulator 1), a paralog of Rds2 and Ert1. Genome-wide analysis showed that Gsm1 has a DNA binding profile highly similar to Rds2. Binding of Gsm1 and Rds2 is interdependent at the gluconeogenic gene FBP1. However, Rds2 is required for Gsm1 to bind at other promoters but not the reverse. Gsm1 and Rds2 also bind to DNA independently of each other. Western blot analysis revealed that Rds2 controls expression of Gsm1. In addition, we showed that the DNA binding domains of Gsm1 and Rds2 bind cooperatively in vitro to the FBP1 promoter. In contrast, at the HAP4 gene, Ert1 cooperates with Rds2 for DNA binding. Mutational analysis suggests that Gsm1/Rds2 and Ert1/Rds2 bind to short common DNA stretches, revealing a novel mode of binding for this class of factors. Two-point mutations in a HAP4 site convert it to a Gsm1 binding site. Thus, Rds2 controls binding of Gsm1 at many promoters by two different mechanisms: regulation of Gsm1 levels and increased DNA binding by formation of heterodimers.

Mitotic Spindle Positioning (MISP) Facilitates Colorectal Cancer Progression by Forming a Complex with Opa Interacting Protein 5 (OIP5) and Activating the JAK2-STAT3 Signaling Pathway

Patients with inflammatory bowel disease (IBD) who experience long-term chronic inflammation of the colon are at an increased risk of developing colorectal cancer (CRC). Mitotic spindle positioning (MISP), an actin-binding protein, plays a role in mitosis and spindle positioning. MISP is found on the apical membrane of the intestinal mucosa and helps stabilize and elongate microvilli, offering protection against colitis. This study explored the role of MISP in colorectal tumorigenesis using a database, human CRC cells, and a mouse model for colitis-induced colorectal tumors triggered by azoxymethane (AOM)/dextran sodium sulfate (DSS) treatment. We found that MISP was highly expressed in colon cancer patient tissues and that reduced MISP expression inhibited cell proliferation. Notably, MISP-deficient mice showed reduced colon tumor formation in the AOM/DSS-induced colitis model. Furthermore, MISP was found to form a complex with Opa interacting protein 5 (OIP5) in the cytoplasm, influencing the expression of OIP5 in a unidirectional manner. We also observed that MISP increased the levels of phosphorylated STAT3 in the JAK2-STAT3 signaling pathway, which is linked to tumorigenesis. These findings indicate that MISP could be a risk factor for CRC, and targeting MISP might provide insights into the mechanisms of colitis-induced colorectal tumorigenesis.

The Rapidly Evolving X-linked miR-506 Family Finetunes Spermatogenesis to Enhance Sperm Competition

Despite rapid evolution across eutherian mammals, the X-linked miR-506 family miRNAs are located in a region flanked by two highly conserved protein-coding genes (Slitrk2 and Fmr1) on the X chromosome. Intriguingly, these miRNAs are predominantly expressed in the testis, suggesting a potential role in spermatogenesis and male fertility. Here, we report that the X-linked miR-506 family miRNAs were derived from the MER91C DNA transposons. Selective inactivation of individual miRNAs or clusters caused no discernable defects, but simultaneous ablation of five clusters containing nineteen members of the miR-506 family led to reduced male fertility in mice. Despite normal sperm counts, motility and morphology, the KO sperm were less competitive than wild-type sperm when subjected to a polyandrous mating scheme. Transcriptomic and bioinformatic analyses revealed that these X-linked miR-506 family miRNAs, in addition to targeting a set of conserved genes, have more targets that are critical for spermatogenesis and embryonic development during evolution. Our data suggest that the miR-506 family miRNAs function to enhance sperm competitiveness and reproductive fitness of the male by finetuning gene expression during spermatogenesis.

Computational prediction and experimental validation identify functionally conserved lncRNAs from zebrafish to human

Functional studies of long noncoding RNAs (lncRNAs) have been hindered by the lack of methods to assess their evolution. Here we present lncRNA Homology Explorer (lncHOME), a computational pipeline that identifies a unique class of long noncoding RNAs (lncRNAs) with conserved genomic locations and patterns of RNA-binding protein (RBP) binding sites (coPARSE-lncRNAs). Remarkably, several hundred human coPARSE-lncRNAs can be evolutionarily traced to zebrafish. Using CRISPR-Cas12a knockout and rescue assays, we found that knocking out many human coPARSE-lncRNAs led to cell proliferation defects, which were subsequently rescued by predicted zebrafish homologs. Knocking down coPARSE-lncRNAs in zebrafish embryos caused severe developmental delays that were rescued by human homologs. Furthermore, we verified that human, mouse and zebrafish coPARSE-lncRNA homologs tend to bind similar RBPs with their conserved functions relying on specific RBP-binding sites. Overall, our study demonstrates a comprehensive approach for studying the functional conservation of lncRNAs and implicates numerous lncRNAs in regulating vertebrate physiology.

PTEN-regulated PI3K-p110 and AKT isoform plasticity controls metastatic prostate cancer progression

PTEN loss, one of the most frequent mutations in prostate cancer (PC), is presumed to drive disease progression through AKT activation. However, two transgenic PC models with Akt activation plus Rb loss exhibited different metastatic development: Pten/RbPE:-/- mice produced systemic metastatic adenocarcinomas with high AKT2 activation, whereas RbPE:-/- mice deficient for the Src-scaffolding protein, Akap12, induced high-grade prostatic intraepithelial neoplasias and indolent lymph node dissemination, correlating with upregulated phosphotyrosyl PI3K-p85α. Using PC cells isogenic for PTEN, we show that PTEN-deficiency correlated with dependence on both p110β and AKT2 for in vitro and in vivo parameters of metastatic growth or motility, and with downregulation of SMAD4, a known PC metastasis suppressor. In contrast, PTEN expression, which dampened these oncogenic behaviors, correlated with greater dependence on p110α plus AKT1. Our data suggest that metastatic PC aggressiveness is controlled by specific PI3K/AKT isoform combinations influenced by divergent Src activation or PTEN-loss pathways.

Molecular identification of missense variants in SLC3A1 gene; an approach leading to computer-aided drug design for cystinuria

Cystinuria is a rare congenital disorder characterized by the formation of cystine stones in urinary system, mainly kidneys and urinary tract. It follows the autosomal recessive inheritance pattern, where both of the parents contain the mutant allele. Cystine is an oxidized dimeric form of amino acid cysteine, shining crystal of greenish-yellow color sized greater than 5 mm. A minor genetic defect in SLC3A1 gene, downregulate the cystine transporter, rBAT protein, to absorb cystine and other dibasic amino acids in proximal tubule of nephron, causing Cystinuria. Computational and molecular analysis of SLC3A1 gene was performed to identify the deleterious missense single nucleotide variations (mSNVs) linked with Cystinuria in Pakistani population. In silico analysis of whole SLC3A1 gene nsSNPs has revealed that the exon 1, 6 and 10 are the hotspot areas, which potentially alter the protein structure and function. Three SNVs including one synonymous SNV A186C in exon 1, and two mSNVs including G314T in exon 1 and G44972A in exons 10 were identified. Both mSNVs were confirmed by ARMS PCR in all the 68 samples. The results have shown that 10% of the patients have G314T, 16% have G44972A and 74% of the patients have both of these mSNVs. Both of these mSNVs were involved in the structural and functional deterioration of rBAT protein. Additionally, computer aided drug designing tools were used to design diaminobenzylpyrimidine drug around the mutant residues which exhibit the lowest binding affinity with the target as compared to the previously reported cystine binding thiol drugs. In future, the present study could be extended to a large scale for mass screening of reported SNVs and mSNVs which, ultimately, lead to the development of knockouts for the functional studies and treatments.

Epigenetic Dysregulation and Its Correlation with the Steroidogenic Machinery Impacting Breast Pathogenesis: Data Mining and Molecular Insights into Therapeutics

Breast cancer (BC) is a heterogeneous condition and comprises molecularly distinct subtypes. An imbalance in the levels of epigenetic histone deacetylases (HDACs), modulating estrogen accumulation, especially 17β-estradiol (E2), promotes breast tumorigenesis. In the present study, analyses of The Cancer Genome Atlas (TCGA) pan-cancer normalized RNA-Seq datasets revealed the dysregulation of 16 epigenetic enzymes (among a total of 18 members) in luminal BC subtypes, in comparison to their non-cancerous counterparts. Explicitly, genomic profiling of these epigenetic enzymes displayed increases in HDAC1, 2, 8, 10, 11, and Sirtuins (SIRTs) 6 and 7, and decreases in HDAC4-7, -9, and SIRT1-4 levels, respectively, in TCGA breast tumors. Kaplan-Meier plot analyses showed that these HDACs, with the exception of HDAC2 and SIRT2, were not correlated with the overall survival of BC patients. Additionally, disruption of the epigenetic signaling in TCGA BC subtypes, as assessed using both heatmaps and boxplots, was associated with the genomic expression of factors that are instrumental for cholesterol trafficking/utilization for accelerating estrogen/E2 levels, in which steroidogenic acute regulatory protein (STAR) mediates the rate-limiting step in steroid biosynthesis. TCGA breast samples showed diverse expression patterns of a variety of key steroidogenic markers and hormone receptors, including LIPE, CYP27A1, STAR, STARD3, CYP11A1, CYP19A1, ER, PGR, and ERBB2. Moreover, regulation of STAR-governed steroidogenic machinery was found to be influenced by various transcription factors, i.e., CREB1, CREM, SF1, NR4A1, CEBPB, SREBF1, SREBF2, SP1, FOS, JUN, NR0B1, and YY1. Along these lines, ingenuity pathway analysis (IPA) recognized a number of new targets and downstream effectors influencing BCs. Of note, genomic, epigenomic, transcriptional, and hormonal anomalies observed in human primary breast tumors were qualitatively similar in pertinent BC cell lines. These findings identify the functional correlation between dysregulated epigenetic enzymes and estrogen/E2 accumulation in human breast tumors, providing the molecular insights into more targeted therapeutic approaches involving the inhibition of HDACs for combating this life-threatening disease.

Rare variants in long non-coding RNAs are associated with blood lipid levels in the TOPMed whole-genome sequencing study

Long non-coding RNAs (lncRNAs) are known to perform important regulatory functions in lipid metabolism. Large-scale whole-genome sequencing (WGS) studies and new statistical methods for variant set tests now provide an opportunity to assess more associations between rare variants in lncRNA genes and complex traits across the genome. In this study, we used high-coverage WGS from 66,329 participants of diverse ancestries with measurement of blood lipids and lipoproteins (LDL-C, HDL-C, TC, and TG) in the National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) program to investigate the role of lncRNAs in lipid variability. We aggregated rare variants for 165,375 lncRNA genes based on their genomic locations and conducted rare-variant aggregate association tests using the STAAR (variant-set test for association using annotation information) framework. We performed STAAR conditional analysis adjusting for common variants in known lipid GWAS loci and rare-coding variants in nearby protein-coding genes. Our analyses revealed 83 rare lncRNA variant sets significantly associated with blood lipid levels, all of which were located in known lipid GWAS loci (in a ±500-kb window of a Global Lipids Genetics Consortium index variant). Notably, 61 out of 83 signals (73%) were conditionally independent of common regulatory variation and rare protein-coding variation at the same loci. We replicated 34 out of 61 (56%) conditionally independent associations using the independent UK Biobank WGS data. Our results expand the genetic architecture of blood lipids to rare variants in lncRNAs.

HK3 stimulates immune cell infiltration to promote glioma deterioration

BACKGROUND

Glioma is the most common and lethal type of brain tumor, and it is characterized by unfavorable prognosis and high recurrence rates. The reprogramming of energy metabolism and an immunosuppressive tumor microenvironment (TME) are two hallmarks of tumors. Complex and dynamic interactions between neoplastic cells and the surrounding microenvironment can generate an immunosuppressive TME, which can accelerate the malignant progression of glioma. Therefore, it is crucial to explore associations between energy metabolism and the immunosuppressive TME and to identify new biomarkers for glioma prognosis.

METHODS

In our work, we analyzed the co-expression relationship between glycolytic genes and immune checkpoints based on the transcriptomic data from The Cancer Genome Atlas (TCGA) and Chinese Glioma Genome Atlas (CGGA) and found the correlation between HK3 expression and glioma tumor immune status. To investigate the biological role of HK3 in glioma, we performed bioinformatics analysis and established a mouse glioblastoma (GBM) xenograft model.

RESULTS

Our study showed that HK3 significantly stimulated immune cell infiltration into the glioma TME. Tissue samples with higher HK3 expressive level showed increasing levels of immune cells infiltration, including M2 macrophages, neutrophils, and various subtypes of activated memory CD4+ T cells. Furthermore, HK3 expression was significantly increasing along with the elevated tumor grade, had a higher level in the mesenchymal subtype compared with those in other subtypes of GBM and could independently predict poor outcomes of GBM patients.

CONCLUSION

The present work mainly concentrated on the biological role of HK3 in glioma and offered a novel insight of HK3 regulating the activation of immune cells in the glioma microenvironment. These findings could provide a new theoretical evidence for understanding the metabolic molecular within the glioma microenvironment and identifying new therapeutic targets.

CAG repeat expansion in the Huntington's disease gene shapes linear and circular RNAs biogenesis

Alternative splicing (AS) appears to be altered in Huntington's disease (HD), but its significance for early, pre-symptomatic disease stages has not been inspected. Here, taking advantage of Htt CAG knock-in mouse in vitro and in vivo models, we demonstrate a correlation between Htt CAG repeat length and increased aberrant linear AS, specifically affecting neural progenitors and, in vivo, the striatum prior to overt behavioral phenotypes stages. Remarkably, a significant proportion (36%) of the aberrantly spliced isoforms are not-functional and meant to non-sense mediated decay (NMD). The expanded Htt CAG repeats further reflect on a previously neglected, global impairment of back-splicing, leading to decreased circular RNAs production in neural progenitors. Integrative transcriptomic analyses unveil a network of transcriptionally altered micro-RNAs and RNA-binding proteins (Celf, hnRNPs, Ptbp, Srsf, Upf1, Ythd2) which might influence the AS machinery, primarily in neural cells. We suggest that this unbalanced expression of linear and circular RNAs might alter neural fitness, contributing to HD pathogenesis.

KDM6A epigenetically regulates subtype plasticity in small cell lung cancer

Small cell lung cancer (SCLC) exists broadly in four molecular subtypes: ASCL1, NEUROD1, POU2F3 and Inflammatory. Initially, SCLC subtypes were thought to be mutually exclusive, but recent evidence shows intra-tumoural subtype heterogeneity and plasticity between subtypes. Here, using a CRISPR-based autochthonous SCLC genetically engineered mouse model to study the consequences of KDM6A/UTX inactivation, we show that KDM6A inactivation induced plasticity from ASCL1 to NEUROD1 resulting in SCLC tumours that express both ASCL1 and NEUROD1. Mechanistically, KDM6A normally maintains an active chromatin state that favours the ASCL1 subtype with its loss decreasing H3K4me1 and increasing H3K27me3 at enhancers of neuroendocrine genes leading to a cell state that is primed for ASCL1-to-NEUROD1 subtype switching. This work identifies KDM6A as an epigenetic regulator that controls ASCL1 to NEUROD1 subtype plasticity and provides an autochthonous SCLC genetically engineered mouse model to model ASCL1 and NEUROD1 subtype heterogeneity and plasticity, which is found in 35-40% of human SCLCs.

Identification of Lutheran Blood Groups and Genetic Variants within KLF1 among Thai Blood Donors

Background

Lua and Lub are inherited as codominant allelic characters resulting from a single nucleotide variant (SNV) of the basal cell adhesion molecule (BCAM) gene. Red cells of the dominantly inherited suppressor of the Lutheran antigens In(Lu) phenotypically appear as Lu(a-b-) by the haemagglutination test. In(Lu) resulted from heterozygosity for mutations within the erythroid-specific Krüppel-like factor 1 (KLF1) gene. This study aimed to determine the frequency of the Lu(a) and Lu(b) phenotypes and genotypes and genetic variants of the distinct In(Lu) among Thai blood donors.

Material and Methods

Samples from 334 Thai donors were phenotyped with anti-Lua and anti-Lub. These DNA samples and an additional 1,370 donor DNA samples with unknown Lu(a)/Lu(b) phenotypes were genotyped using an in-house PCR-SSP. In the case of the three Lu(a-b-) donors, the BCAM and KLF1 genes were analysed by PCR and sequencing.

Results

A total of 331 of the 334 donors were Lu(a-b+), while the other observed phenotype, appearing as Lu(a-b-), was found among three donors. Of those three Lu(a-b-) donors with the LU*02/02 genotype, we identified KLF1 variant alleles, consisting of two variants: c.[304T>C, 1001C>G] and c.[304T>C, 519_525dupCGGCGCC], leading to the In(Lu) phenotype, and one homozygous variant (c.304T>C) mutation. Also, only one Thai donor was genotyped as LU*01/02, confirmed by serology test and DNA sequencing.

Conclusion

In this study, we identified KLF1 variants to be included in Lutheran typing analysis in Thai populations. Therefore, the application of genotyping and phenotyping methods has simultaneously been in use to screen and confirm the rare Lu(a+) and In(Lu) phenotypes.

Integrative mapping of the dog epigenome: Reference annotation for comparative intertissue and cross-species studies

Dogs have become a valuable model in exploring multifaceted diseases and biology relevant to human health. Despite large-scale dog genome projects producing high-quality draft references, a comprehensive annotation of functional elements is still lacking. We addressed this through integrative next-generation sequencing of transcriptomes paired with five histone marks and DNA methylome profiling across 11 tissue types, deciphering the dog's epigenetic code by defining distinct chromatin states, super-enhancer, and methylome landscapes, and thus showed that these regions are associated with a wide range of biological functions and cell/tissue identity. In addition, we confirmed that the phenotype-associated variants are enriched in tissue-specific regulatory regions and, therefore, the tissue of origin of the variants can be traced. Ultimately, we delineated conserved and dynamic epigenomic changes at the tissue- and species-specific resolutions. Our study provides an epigenomic blueprint of the dog that can be used for comparative biology and medical research.

Rare variants in long non-coding RNAs are associated with blood lipid levels in the TOPMed Whole Genome Sequencing Study

Long non-coding RNAs (lncRNAs) are known to perform important regulatory functions. Large-scale whole genome sequencing (WGS) studies and new statistical methods for variant set tests now provide an opportunity to assess the associations between rare variants in lncRNA genes and complex traits across the genome. In this study, we used high-coverage WGS from 66,329 participants of diverse ancestries with blood lipid levels (LDL-C, HDL-C, TC, and TG) in the National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) program to investigate the role of lncRNAs in lipid variability. We aggregated rare variants for 165,375 lncRNA genes based on their genomic locations and conducted rare variant aggregate association tests using the STAAR (variant-Set Test for Association using Annotation infoRmation) framework. We performed STAAR conditional analysis adjusting for common variants in known lipid GWAS loci and rare coding variants in nearby protein coding genes. Our analyses revealed 83 rare lncRNA variant sets significantly associated with blood lipid levels, all of which were located in known lipid GWAS loci (in a ±500 kb window of a Global Lipids Genetics Consortium index variant). Notably, 61 out of 83 signals (73%) were conditionally independent of common regulatory variations and rare protein coding variations at the same loci. We replicated 34 out of 61 (56%) conditionally independent associations using the independent UK Biobank WGS data. Our results expand the genetic architecture of blood lipids to rare variants in lncRNA, implicating new therapeutic opportunities.

PTEN regulated PI3K-p110 and AKT isoform plasticity controls metastatic prostate cancer progression

PTEN loss, one of the most frequent mutations in prostate cancer (PC), is presumed to drive disease progression through AKT activation. However, two transgenic PC models with Akt activation plus Rb loss exhibited different metastasis development: Pten/RbPE:-/- mice produced systemic metastatic adenocarcinomas with high AKT2 activation, whereas RbPE:-/- mice deficient for the Src-scaffolding protein, Akap12, induced high-grade prostatic intraepithelial neoplasias and indolent lymph node disseminations, correlating with upregulated phosphotyrosyl PI3K-p85α. Using PC cells isogenic for PTEN, we show that PTEN-deficiency correlated with dependence on both p110β and AKT2 for in vitro and in vivo parameters of metastatic growth or motility, and with downregulation of SMAD4, a known PC metastasis suppressor. In contrast, PTEN expression, which dampened these oncogenic behaviors, correlated with greater dependence on p110α plus AKT1. Our data suggest that metastatic PC aggressiveness is controlled by specific PI3K/AKT isoform combinations influenced by divergent Src activation or PTEN-loss pathways.

ChroKit: a Shiny-based framework for interactive analysis, visualization and integration of genomic data

We developed ChroKit (the Chromatin toolKit), an interactive web-based framework written in R that enables intuitive exploration, multidimensional analyses, and visualization of genomic data from ChIP-Seq, DNAse-Seq or any other NGS experiment that reports the enrichment of aligned reads over genomic regions. This program takes preprocessed NGS data and performs operations on genomic regions of interest, including resetting their boundaries, their annotation based on proximity to genomic features, the association to gene ontologies, and signal enrichment calculations. Genomic regions can be further refined or subsetted by user-defined logical operations and unsupervised classification algorithms. ChroKit generates a full range of plots that are easily manipulated by point and click operations, thus allowing 'on the fly' re-analysis and fast exploration of the data. Working sessions can be exported for reproducibility, accountability, and easy sharing within the bioinformatics community. ChroKit is multiplatform and can be deployed on a server to enhance computational speed and provide simultaneous access by multiple users. ChroKit is a fast and intuitive genomic analysis tool suited for a wide range of users due to its architecture and its user-friendly graphical interface. ChroKit source code is available at https://github.com/ocroci/ChroKit and the Docker image at https://hub.docker.com/r/ocroci/chrokit.

Leveraging base-pair mammalian constraint to understand genetic variation and human disease

Thousands of genomic regions have been associated with heritable human diseases, but attempts to elucidate biological mechanisms are impeded by an inability to discern which genomic positions are functionally important. Evolutionary constraint is a powerful predictor of function, agnostic to cell type or disease mechanism. Single-base phyloP scores from 240 mammals identified 3.3% of the human genome as significantly constrained and likely functional. We compared phyloP scores to genome annotation, association studies, copy-number variation, clinical genetics findings, and cancer data. Constrained positions are enriched for variants that explain common disease heritability more than other functional annotations. Our results improve variant annotation but also highlight that the regulatory landscape of the human genome still needs to be further explored and linked to disease.

Mammalian evolution of human cis-regulatory elements and transcription factor binding sites

Understanding the regulatory landscape of the human genome is a long-standing objective of modern biology. Using the reference-free alignment across 241 mammalian genomes produced by the Zoonomia Consortium, we charted evolutionary trajectories for 0.92 million human candidate cis-regulatory elements (cCREs) and 15.6 million human transcription factor binding sites (TFBSs). We identified 439,461 cCREs and 2,024,062 TFBSs under evolutionary constraint. Genes near constrained elements perform fundamental cellular processes, whereas genes near primate-specific elements are involved in environmental interaction, including odor perception and immune response. About 20% of TFBSs are transposable element-derived and exhibit intricate patterns of gains and losses during primate evolution whereas sequence variants associated with complex traits are enriched in constrained TFBSs. Our annotations illuminate the regulatory functions of the human genome.

Elevated enhancer-oncogene contacts and higher oncogene expression levels by recurrent CTCF inactivating mutations in acute T cell leukemia

Monoallelic inactivation of CCCTC-binding factor (CTCF) in human cancer drives altered methylated genomic states, altered CTCF occupancy at promoter and enhancer regions, and deregulated global gene expression. In patients with T cell acute lymphoblastic leukemia (T-ALL), we find that acquired monoallelic CTCF-inactivating events drive subtle and local genomic effects in nearly half of t(5; 14) (q35; q32.2) rearranged patients, especially when CTCF-binding sites are preserved in between the BCL11B enhancer and the TLX3 oncogene. These solitary intervening sites insulate TLX3 from the enhancer by inducing competitive looping to multiple binding sites near the TLX3 promoter. Reduced CTCF levels or deletion of the intervening CTCF site abrogates enhancer insulation by weakening competitive looping while favoring TLX3 promoter to BCL11B enhancer looping, which elevates oncogene expression levels and leukemia burden.

Phospholipase A and acyltransferase 4/retinoic acid receptor responder 3 at the intersection of tumor suppression and pathogen restriction

Phospholipase A and acyltransferase (PLAAT) 4 is a class II tumor suppressor with phospholipid metabolizing abilities. It was characterized in late 2000s, and has since been referred to as 'tazarotene-induced gene 3' (TIG3) or 'retinoic acid receptor responder 3' (RARRES3) as a key downstream effector of retinoic acid signaling. Two decades of research have revealed the complexity of its function and regulatory roles in suppressing tumorigenesis. However, more recent findings have also identified PLAAT4 as a key anti-microbial effector enzyme acting downstream of interferon regulatory factor 1 (IRF1) and interferons (IFNs), favoring protection from virus and parasite infections. Unveiling the molecular mechanisms underlying its action may thus open new therapeutic avenues for the treatment of both cancer and infectious diseases. Herein, we aim to summarize a brief history of PLAAT4 discovery, its transcriptional regulation, and the potential mechanisms in tumor prevention and anti-pathogen defense, and discuss potential future directions of PLAAT4 research toward the development of therapeutic approaches targeting this enzyme with pleiotropic functions.

Leveraging Base Pair Mammalian Constraint to Understand Genetic Variation and Human Disease

Although thousands of genomic regions have been associated with heritable human diseases, attempts to elucidate biological mechanisms are impeded by a general inability to discern which genomic positions are functionally important. Evolutionary constraint is a powerful predictor of function that is agnostic to cell type or disease mechanism. Here, single base phyloP scores from the whole genome alignment of 240 placental mammals identified 3.5% of the human genome as significantly constrained, and likely functional. We compared these scores to large-scale genome annotation, genome-wide association studies (GWAS), copy number variation, clinical genetics findings, and cancer data sets. Evolutionarily constrained positions are enriched for variants explaining common disease heritability (more than any other functional annotation). Our results improve variant annotation but also highlight that the regulatory landscape of the human genome still needs to be further explored and linked to disease.

A CpG island-encoded mechanism protects genes from premature transcription termination

Transcription must be tightly controlled to regulate gene expression and development. However, our understanding of the molecular mechanisms that influence transcription and how these are coordinated in cells to ensure normal gene expression remains rudimentary. Here, by dissecting the function of the SET1 chromatin-modifying complexes that bind to CpG island-associated gene promoters, we discover that they play a specific and essential role in enabling the expression of low to moderately transcribed genes. Counterintuitively, this effect can occur independently of SET1 complex histone-modifying activity and instead relies on an interaction with the RNA Polymerase II-binding protein WDR82. Unexpectedly, we discover that SET1 complexes enable gene expression by antagonising premature transcription termination by the ZC3H4/WDR82 complex at CpG island-associated genes. In contrast, at extragenic sites of transcription, which typically lack CpG islands and SET1 complex occupancy, we show that the activity of ZC3H4/WDR82 is unopposed. Therefore, we reveal a gene regulatory mechanism whereby CpG islands are bound by a protein complex that specifically protects genic transcripts from premature termination, effectively distinguishing genic from extragenic transcription and enabling normal gene expression.

Long noncoding TSI attenuates aortic valve calcification by suppressing TGF-β1-induced osteoblastic differentiation of valve interstitial cells

INTRODUCTION

Calcific aortic valve disease (CAVD) is an active and cellular-driven fibrocalcific process characterised by differentiation of valve interstitial cells (VICs) towards an osteogenic-like phenotype. A recently identified lncRNA, lncTSI, has been reported to inhibit fibrogenesis through transforming growth factor (TGF)-β/Smad3 pathway. Here, the present study aimed to investigate the role of lncTSI in CAVD.

METHODS

The effect of TGF-β1 on lncTSI of VICs was measured. TGF-β1, RUNX2 and collagen I expression between calcified aortic valve tissue and normal samples by immunohistochemistry and western blotting. Human VICs were cultured and treated with TGF-β1. SiRNA and pcDNA3.1-lncTSI plasmid transfection were used to silence and overexpress lncTSI in VICs for 48 h, Smads phosphorylation, RUNX2 and collagen I expression were then verified by western blotting. In ApoE^-/- mice fed with 0.25 % high-cholesterol diet, AAV2-lncTSI were injected intravenously to observe their effect on the formation of aortic valve calcification.

RESULTS

lncTSI was highly expressed in VICs treated with TGF-β1. lncTSI was transcriptionally regulated by Smad3 and reversely inhibited TGF-β1-induced Smad3 phosphorylation and downregulated profibrotic gene expression. Silencing lncTSI increased TGF-β1-induced Smad3 phosphorylation, and subsequently, upregulated RUNX2 and collagen I expressions in VICs. While overexpression of lncTSI reversed the production of RUNX2 and collagen I in VICs. In a mouse CAVD model of 24 week 0.25 % high-cholesterol diet feeding, overexpression of lncTSI significantly reduced calcium deposition, RUNX2, pSmad3, and collagen I expression in aortic valve leaflets, with less aortic valve stenosis.

CONCLUSIONS

The novel findings of present study suggested that lncTSI alleviated aortic valve calcification through negative regulation of the TGF-β/Smad3 pathway. The results may help elucidate new diagnostic and therapeutic targets to prevent CAVD progression.

Changes of tRNA-Derived Fragments by Alzheimer's Disease in Cerebrospinal Fluid and Blood Serum

BACKGROUND

Alzheimer's disease (AD) is the most common type of dementia, affecting individuals over 65. AD is also a multifactorial disease, with disease mechanisms incompletely characterized, and disease-modifying therapies are marginally effective. Biomarker signatures may shed light on the diagnosis, disease mechanisms, and the development of therapeutic targets. tRNA-derived RNA fragments (tRFs), a family of recently discovered small non-coding RNAs, have been found to be significantly enhanced in human AD hippocampus tissues. However, whether tRFs change in body fluids is unknown.

OBJECTIVE

To investigate whether tRFs in body fluids are impacted by AD.

METHODS

We first used T4 polynucleotide kinase-RNA-seq, a modified next-generation sequencing technique, to identify detectable tRFs in human cerebrospinal fluid and serum samples. The detectable tRFs were then compared in these fluids from control, AD, and mild cognitive impairment patients using tRF qRT-PCR. The stability of tRFs in serum was also investigated by checking the change in tRFs in response to protein digestion or exosome lysis.

RESULTS

Among various tRFs, tRF5-ProAGG seemed to be impacted by AD in both cerebrospinal fluid and serum. AD-impacted serum tRF5-ProAGG showed a correlation with the AD stage. Putative targets of tRF5-ProAGG in the hippocampus were also predicted by a computational algorithm, with some targets being validated experimentally and one of them being in a negative correlation with tRF5-ProAGG even using a small size of samples.

CONCLUSIONS

tRF5-ProAGG showed the potential as an AD biomarker and may play a role in disease progression.

DriverMP enables improved identification of cancer driver genes

BACKGROUND

Cancer is widely regarded as a complex disease primarily driven by genetic mutations. A critical concern and significant obstacle lies in discerning driver genes amid an extensive array of passenger genes.

FINDINGS

We present a new method termed DriverMP for effectively prioritizing altered genes on a cancer-type level by considering mutated gene pairs. It is designed to first apply nonsilent somatic mutation data, protein‒protein interaction network data, and differential gene expression data to prioritize mutated gene pairs, and then individual mutated genes are prioritized based on prioritized mutated gene pairs. Application of this method in 10 cancer datasets from The Cancer Genome Atlas demonstrated its great improvements over all the compared state-of-the-art methods in identifying known driver genes. Then, a comprehensive analysis demonstrated the reliability of the novel driver genes that are strongly supported by clinical experiments, disease enrichment, or biological pathway analysis.

CONCLUSIONS

The new method, DriverMP, which is able to identify driver genes by effectively integrating the advantages of multiple kinds of cancer data, is available at https://github.com/LiuYangyangSDU/DriverMP. In addition, we have developed a novel driver gene database for 10 cancer types and an online service that can be freely accessed without registration for users. The DriverMP method, the database of novel drivers, and the user-friendly online server are expected to contribute to new diagnostic and therapeutic opportunities for cancers.

Noninvasive prediction of axillary lymph node status in breast cancer using promoter profiling of circulating cell-free DNA

BACKGROUND

Lymph node metastasis (LNM) is one of the most important factors affecting the prognosis of breast cancer. The accurate evaluation of lymph node status is useful to predict the outcomes of patients and guide the choice of cancer treatment. However, there is still lack of a low-cost non-invasive method to assess the status of axillary lymph node (ALN). Gene expression signature has been used to assess lymph node metastasis status of breast cancer. In addition, nucleosome footprint of cell-free DNA (cfDNA) carries gene expression information of its original tissues, so it may be used to evaluate the axillary lymph node status in breast cancer.

METHODS

In this study, we found that the cfDNA nucleosome footprints between the ALN-positive patients and ALN-negative patients showed different patterns by implementing whole-genome sequencing (WGS) to detect 15 ALN-positive and 15 ALN-negative patients. In order to further evaluate its potential for assessing ALN status, we developed a classifier with multiple machine learning models by using 330 WGS data of cfDNA from 162 ALN-positive and 168 ALN-negative samples to distinguish these two types of patients.

RESULTS

We found that the promoter profiling between the ALN-positive patients and ALN-negative patients showed distinct patterns. In addition, we observed 1071 genes with differential promoter coverage and their functions were closely related to tumorigenesis. We found that the predictive classifier based on promoter profiling with a support vector machine model, named PPCNM, produced the largest area under the curve of 0.897 (95% confidence interval 0.86-0.93).

CONCLUSIONS

These results indicate that promoter profiling can be used to distinguish ALN-positive patients from ALN-negative patients, which may be helpful to guide the choice of cancer treatment.

Methylation-Mediated Silencing of RBP7 Promotes Breast Cancer Progression through PPAR and PI3K/AKT Pathway

Retinoid-binding protein7 (RBP7) is a member of the cellular retinol-binding protein (CRBP) family, which is involved in the pathogenesis of breast cancer. The study aims to illustrate the prognostic value and the potential regulatory mechanisms of RBP7 expression in breast cancer. Bioinformatics analysis with the TCGA and CPTAC databases revealed that the mRNA and protein expression levels of RBP7 in normal were higher compared to breast cancer tissues. Survival analysis displayed that the lower expression of RBP7, the worse the prognosis in ER-positive (ER+) breast cancer patients. Genomic analysis showed that low expression of RBP7 correlates with its promoter hypermethylation in breast cancer. Functional enrichment analysis demonstrated that downregulation of RBP7 expression may exert its biological influence on breast cancer through the PPAR pathway and the PI3K/AKT pathway. In summary, we identified RBP7 as a novel biomarker that is helpful for the prognosis of ER+ breast cancer patients. Promoter methylation of RBP7 is involved in its gene silencing in breast cancer, thus regulating the occurrence and development of ER+ breast cancer through the PPAR and PI3K/AKT pathways.

Dynamic Expression of Palmitoylation Regulators across Human Organ Development and Cancers Based on Bioinformatics

Protein palmitoylation is a reversible modification process that links palmitate to cysteine residues via a reversible thioester bond. Palmitoylation exerts an important role in human organ development and tumor progression. However, a comprehensive landscape regarding the dynamic expression of palmitoylation regulators in human organ development remains unclear. In this study, we analyzed the dynamic expression of palmitoylation regulators in seven organ development and eight cancer types based on bioinformatics. We found that the expression levels of most palmitoylation regulators were altered after birth. In particular, ZDHHC7/20/21 exhibited converse expression patterns in multiple cancer types. Survival analysis showed that the poor prognosis in patients with kidney renal clear carcinoma (KIRC) is related to low expression of ZDHHC7/20/21, and a high expression of ZDHHC7/20/21 is related to worse survival in patients with liver hepatocellular carcinoma (LIHC). Furthermore, we found that the expression of ZDHHC7 is associated with infiltration levels of some types of immune cells in the tumor microenvironment (TME), and we explored the relationship between ZDHHC7 expression and immune checkpoint (ICP) genes across 33 cancer types. In addition, gene set enrichment analysis (GSEA) results indicated that ZDHHC7 might regulate different genes to mediate the same pathway in different organs. In summary, the comprehensive analysis of palmitoylation regulators reveals their functions in human organ development and cancer, which may provide new insights for developing new tumor markers.

Integrative Analysis Between Genome-Wide Association Study and Expression Quantitative Trait Loci Reveals Bovine Muscle Gene Expression Regulatory Polymorphisms Associated With Intramuscular Fat and Backfat Thickness

Understanding the architecture of gene expression is fundamental to unravel the molecular mechanisms regulating complex traits in bovine, such as intramuscular fat content (IMF) and backfat thickness (BFT). These traits are economically important for the beef industry since they affect carcass and meat quality. Our main goal was to identify gene expression regulatory polymorphisms within genomic regions (QTL) associated with IMF and BFT in Nellore cattle. For that, we used RNA-Seq data from 193 Nellore steers to perform SNP calling analysis. Then, we combined the RNA-Seq SNP and a high-density SNP panel to obtain a new dataset for further genome-wide association analysis (GWAS), totaling 534,928 SNPs. GWAS was performed using the Bayes B model. Twenty-one relevant QTL were associated with our target traits. The expression quantitative trait loci (eQTL) analysis was performed using Matrix eQTL with the complete SNP dataset and 12,991 genes, revealing a total of 71,033 cis and 36,497 trans-eQTL (FDR < 0.05). Intersecting with QTL for IMF, we found 231 eQTL regulating the expression levels of 117 genes. Within those eQTL, three predicted deleterious SNPs were identified. We also identified 109 eQTL associated with BFT and affecting the expression of 54 genes. This study revealed genomic regions and regulatory SNPs associated with fat deposition in Nellore cattle. We highlight the transcription factors FOXP4, FOXO3, ZSCAN2, and EBF4, involved in lipid metabolism-related pathways. These results helped us to improve our knowledge about the genetic architecture behind important traits in cattle.

Insights From a Large-Scale Whole-Genome Sequencing Study of Systolic Blood Pressure, Diastolic Blood Pressure, and Hypertension

BACKGROUND

The availability of whole-genome sequencing data in large studies has enabled the assessment of coding and noncoding variants across the allele frequency spectrum for their associations with blood pressure.

METHODS

We conducted a multiancestry whole-genome sequencing analysis of blood pressure among 51 456 Trans-Omics for Precision Medicine and Centers for Common Disease Genomics program participants (stage-1). Stage-2 analyses leveraged array data from UK Biobank (N=383 145), Million Veteran Program (N=318 891), and Reasons for Geographic and Racial Differences in Stroke (N=10 643) participants, along with whole-exome sequencing data from UK Biobank (N=199 631) participants.

RESULTS

Two blood pressure signals achieved genome-wide significance in meta-analyses of stage-1 and stage-2 single variant findings (P<5×10-8). Among them, a rare intergenic variant at novel locus, LOC100506274, was associated with lower systolic blood pressure in stage-1 (beta [SE]=-32.6 [6.0]; P=4.99×10-8) but not stage-2 analysis (P=0.11). Furthermore, a novel common variant at the known INSR locus was suggestively associated with diastolic blood pressure in stage-1 (beta [SE]=-0.36 [0.07]; P=4.18×10-7) and attained genome-wide significance in stage-2 (beta [SE]=-0.29 [0.03]; P=7.28×10-23). Nineteen additional signals suggestively associated with blood pressure in meta-analysis of single and aggregate rare variant findings (P<1×10-6 and P<1×10-4, respectively).

DISCUSSION

We report one promising but unconfirmed rare variant for blood pressure and, more importantly, contribute insights for future blood pressure sequencing studies. Our findings suggest promise of aggregate analyses to complement single variant analysis strategies and the need for larger, diverse samples, and family studies to enable robust rare variant identification.

Selection signature analyses and genome-wide association reveal genomic hotspot regions that reflect differences between breeds of horse with contrasting risk of degenerative suspensory ligament desmitis

Degenerative suspensory ligament desmitis is a progressive idiopathic condition that leads to scarring and rupture of suspensory ligament fibers in multiple limbs in horses. The prevalence of degenerative suspensory ligament desmitis is breed related. Risk is high in the Peruvian Horse, whereas pony and draft breeds have low breed risk. Degenerative suspensory ligament desmitis occurs in families of Peruvian Horses, but its genetic architecture has not been definitively determined. We investigated contrasts between breeds with differing risk of degenerative suspensory ligament desmitis and identified associated risk variants and candidate genes. We analyzed 670k single nucleotide polymorphisms from 10 breeds, each of which was assigned one of the four breed degenerative suspensory ligament desmitis risk categories: control (Belgian, Icelandic Horse, Shetland Pony, and Welsh Pony), low risk (Lusitano, Arabian), medium risk (Standardbred, Thoroughbred, Quarter Horse), and high risk (Peruvian Horse). Single nucleotide polymorphisms were used for genome-wide association and selection signature analysis using breed-assigned risk levels. We found that the Peruvian Horse is a population with low effective population size and our breed contrasts suggest that degenerative suspensory ligament desmitis is a polygenic disease. Variant frequency exhibited signatures of positive selection across degenerative suspensory ligament desmitis breed risk groups on chromosomes 7, 18, and 23. Our results suggest degenerative suspensory ligament desmitis breed risk is associated with disturbances to suspensory ligament homeostasis where matrix responses to mechanical loading are perturbed through disturbances to aging in tendon (PIN1), mechanotransduction (KANK1, KANK2, JUNB, SEMA7A), collagen synthesis (COL4A1, COL5A2, COL5A3, COL6A5), matrix responses to hypoxia (PRDX2), lipid metabolism (LDLR, VLDLR), and BMP signaling (GREM2). Our results do not suggest that suspensory ligament proteoglycan turnover is a primary factor in disease pathogenesis.

Natural Killer cells demonstrate distinct eQTL and transcriptome-wide disease associations, highlighting their role in autoimmunity

Natural Killer cells are innate lymphocytes with central roles in immunosurveillance and are implicated in autoimmune pathogenesis. The degree to which regulatory variants affect Natural Killer cell gene expression is poorly understood. Here we perform expression quantitative trait locus mapping of negatively selected Natural Killer cells from a population of healthy Europeans (n = 245). We find a significant subset of genes demonstrate expression quantitative trait loci specific to Natural Killer cells and these are highly informative of human disease, in particular autoimmunity. A Natural Killer cell transcriptome-wide association study across five common autoimmune diseases identifies further novel associations at 27 genes. In addition to these cis observations, we find novel master-regulatory regions impacting expression of trans gene networks at regions including 19q13.4, the Killer cell Immunoglobulin-like Receptor region, GNLY, MC1R and UVSSA. Our findings provide new insights into the unique biology of Natural Killer cells, demonstrating markedly different expression quantitative trait loci from other immune cells, with implications for disease mechanisms.

Natural Killer cells are key mediators of anti-tumour immunosurveillance and anti-viral immunity. Here, the authors map regulatory genetic variation in primary Natural Killer cells, providing new insights into their role in human health and disease.

Profile of TREM2-Derived circRNA and mRNA Variants in the Entorhinal Cortex of Alzheimer’s Disease Patients

Genetic variants in TREM2, a microglia-related gene, are well-known risk factors for Alzheimer’s disease (AD). Here, we report that TREM2 originates from circular RNAs (circRNAs), a novel class of non-coding RNAs characterized by a covalent and stable closed-loop structure. First, divergent primers were designed to amplify circRNAs by RT-PCR, which were further assessed by Sanger sequencing. Then, additional primer sets were used to confirm back-splicing junctions. In addition, HMC3 cells were used to assess the microglial expression of circTREM2s. Three candidate circTREM2s were identified in control and AD human entorhinal samples. One of the circRNAs, circTREM2_1, was consistently amplified by all divergent primer sets in control and AD entorhinal cortex samples as well as in HMC3 cells. In AD cases, a moderate negative correlation (r = −0.434) was found between the global average area of Aβ deposits in the entorhinal cortex and circTREM2_1 expression level. In addition, by bioinformatics tools, a total of 16 miRNAs were predicted to join with circTREM2s. Finally, TREM2 mRNA corresponding to four isoforms was profiled by RT-qPCR. TREM2 mRNA levels were found elevated in entorhinal samples of AD patients with low or intermediate ABC scores compared to controls. To sum up, a novel circRNA derived from the TREM2 gene, circTREM2_1, has been identified in the human entorhinal cortex and TREM2 mRNA expression has been detected to increase in AD compared to controls. Unraveling the molecular genetics of the TREM2 gene may help to better know the innate immune response in AD.

Retrospective Natural History Study of RPGR-Related Cone- and Cone-Rod Dystrophies While Expanding the Mutation Spectrum of the Disease

Variants in the X-linked retinitis pigmentosa GTPase regulator gene (RPGR) and, specifically, in its retinal opening reading frame-15 isoform (RPGRORF15) may cause rod-cone (RCD), cone, and cone-rod dystrophies (CDs and CRDs). While RPGR-related RCDs have been frequently evaluated, the characteristics and progression of RPGR-related CD/CRDs are largely unknown. Therefore, the goal of our work was to perform genotype-phenotype correlations specifically in RPGRORF15-related CD/CRDs. This retrospective longitudinal study included 34 index patients and two affected relatives with a molecular diagnosis of RPGR-related CD/CRDs. Patients were recruited at the "Quinze-Vingts" Hospital, Paris, France and screened for mutations in RPGRORF15 at the Institut de la Vision, Paris, France. We identified 29 distinct variants, of which 27 were truncating. All were located in the 3' half of the RPGRORF15 transcript. Twenty of them were novel. Fifteen subjects were affected by CD, the remaining had CRD. When analyzing the longitudinal data, a progressive decline in visual acuity (VA) was noted, with more than 60% of the patients reaching VA ≥ 1 LogMar in the best eye after the fifth decade of life. To our knowledge, this is the largest described study of a cohort of CD/CRD patients affected by RPGRORF15 variants. Longitudinal data showed a rapidly progressive disease, possibly locating an optimal window of intervention for future therapies in younger ages.

TP53-inducible putative long noncoding RNAs encode functional polypeptides that suppress cell proliferation

Polypeptides encoded by long non-coding RNAs (lncRNAs) are a novel class of functional molecules. However, whether these hidden polypeptides participate in the TP53 pathway and play a significant biological role is still unclear. Here, we discover that TP53-regulated lncRNAs encode peptides, two of which are functional in various human cell lines. Using ribosome profiling and RNA-seq approaches in HepG2 cells, we systematically identified more than 300 novel TP53-regulated lncRNAs and further confirmed that fifteen of these TP53-regulated lncRNAs encode peptides. Furthermore, several peptides were validated by multiple mass spectrometry measures. Ten of the novel translational lncRNAs were directly inducible by TP53 in response to DNA damage. Notably, we showed that the TP53-inducible peptides TP53LC02 and TP53LC04, but not their lncRNAs, could suppress cell proliferation. TP53LC04 peptide also had a function associated with cell proliferation by regulating the cell cycle in response to DNA damage. This study demonstrates that TP53-inducible lncRNAs encode new functional peptides, leading to the enlargement of the components of TP53 tumor suppressor network and providing novel potential targets for cancer therapy.

Circular RNA circ-MTHFD1L induces HR repair to promote gemcitabine resistance via the miR-615-3p/RPN6 axis in pancreatic ductal adenocarcinoma

BACKGROUND

Chemoresistance of pancreatic cancer is the main reason for the poor treatment effect of pancreatic cancer patients. Exploring chemotherapy resistance-related genes has been a difficult and hot topic of oncology. Numerous studies implicate the key roles of circular RNAs (circRNAs) in the development of pancreatic cancer. However, the regulation of circRNAs in the process of pancreatic ductal adenocarcinoma (PDAC) chemotherapy resistance is not yet fully clear.

METHODS

Based on the cross-analysis of the Gene Expression Omnibus (GEO) database and the data of our center, we explored a new molecule, hsa_circ_0078297 (circ-MTHFD1L), related to chemotherapy resistance. QRT-PCR was used to detect the expression of circRNAs, miRNAs, and mRNAs in human PDAC tissues and their matched normal tissues. The interaction between circ-MTHFD1L and miR-615-3p/RPN6 signal axis was confirmed by a series of experiments such as Dual-luciferase reporter assay, fluorescence in situ hybridization (FISH) RNA immunoprecipitation (RIP) assays.

RESULTS

Circ-MTHFD1L was significantly increased in PDAC tissues and cells. And in PDAC patients, the higher the expression level of circ-MTHFD1L, the worse the prognosis. Mechanism analysis showed that circ-MTHFD1L, as an endogenous miR-615-3p sponge, upregulates the expression of RPN6, thereby promoting DNA damage repair and exerting its effect on enhancing gemcitabine chemotherapy resistance. More importantly, we also found that Silencing circ-MTHFD1L combined with olaparib can increase the sensitivity of pancreatic cancer to gemcitabine.

CONCLUSION

Circ-MTHFD1L maintains PDAC gemcitabine resistance through the miR-615-3p/RPN6 signal axis. Circ-MTHFD1L may be a molecular marker for the effective treatment of PDAC.

Cumulative Evidence for Associations Between Genetic Variants in Interleukin 6 Receptor Gene and Human Diseases and Phenotypes

Background

Genetic studies have linked polymorphisms in the interleukin 6 receptor (IL6R) gene to the risk of multiple human diseases and phenotypes, yet have reported inconsistent results. We aimed to synthesize current knowledge of variants in the IL6R gene on the risk of diseases and phenotypes.

Methods

We searched the Medline and Embase databases to identify relevant publications. Meta-analysis was performed utilizing DerSimonian and Laird random-effects model. We also graded cumulative evidence for significant associations. Furthermore, phenome-wide analyses and functional annotations were performed for variants with strong evidence.

Results

We included 155 studies for evaluating the associations between 80 polymorphisms in the IL6R gene and the risk of 102 human diseases and 98 phenotypes. We conducted 58 main meta-analyses, and 41 significant associations were identified. Strong evidence was assigned to 29 associations that investigated ten variants (rs2228145, rs4129267, rs7529229, rs4537545, rs7518199, rs4845625, rs4553185, rs4845618, rs4845371, and rs6667434) related to the risk of four cardiovascular diseases (coronary heart disease, coronary artery disease, atherosclerosis, and abdominal aortic aneurysms), four inflammatory diseases (rheumatoid arthritis, Crohn’s disease, dermatitis, and asthma), and concentration of four phenotypes (C-reactive protein, fibrinogen, IL-6, and sIL-6R). Furthermore, phenome-wide analysis verified that rs2228145 associated with asthma and dermatitis risk. Functional analyses indicated that these polymorphisms fall within exon, enhancer regions.

Conclusions

Our study comprehensively summarizes current data on the genetic architecture of the IL6R gene and highlights the pharmacological targeting potential of IL-6R on cardiovascular and inflammatory diseases.

An Amish founder population reveals rare-population genetic determinants of the human lipidome

Identifying the genetic determinants of inter-individual variation in lipid species (lipidome) may provide deeper understanding and additional insight into the mechanistic effect of complex lipidomic pathways in CVD risk and progression beyond simple traditional lipids. Previous studies have been largely population based and thus only powered to discover associations with common genetic variants. Founder populations represent a powerful resource to accelerate discovery of previously unknown biology associated with rare population alleles that have risen to higher frequency due to genetic drift. We performed a genome-wide association scan of 355 lipid species in 650 individuals from the Amish founder population including 127 lipid species not previously tested. To the best of our knowledge, we report for the first time the lipid species associated with two rare-population but Amish-enriched lipid variants: APOB_rs5742904 and APOC3_rs76353203. We also identified novel associations for 3 rare-population Amish-enriched loci with several sphingolipids and with proposed potential functional/causal variant in each locus including GLTPD2_rs536055318, CERS5_rs771033566, and AKNA_rs531892793. We replicated 7 previously known common loci including novel associations with two sterols: androstenediol with UGT locus and estriol with SLC22A8/A24 locus. Our results show the double power of founder populations and detailed lipidome to discover novel trait-associated variants.

A GWAS of 355 lipid species in the Old Order Amish founder population reveals associations between Amish-enriched loci and several sphingolipids.

Starvation causes changes in the intestinal transcriptome and microbiome that are reversed upon refeeding

BACKGROUND

The ability of animals and their microbiomes to adapt to starvation and then restore homeostasis after refeeding is fundamental to their continued survival and symbiosis. The intestine is the primary site of nutrient absorption and microbiome interaction, however our understanding of intestinal adaptations to starvation and refeeding remains limited. Here we used RNA sequencing and 16S rRNA gene sequencing to uncover changes in the intestinal transcriptome and microbiome of zebrafish subjected to long-term starvation and refeeding compared to continuously fed controls.

RESULTS

Starvation over 21 days led to increased diversity and altered composition in the intestinal microbiome compared to fed controls, including relative increases in Vibrio and reductions in Plesiomonas bacteria. Starvation also led to significant alterations in host gene expression in the intestine, with distinct pathways affected at early and late stages of starvation. This included increases in the expression of ribosome biogenesis genes early in starvation, followed by decreased expression of genes involved in antiviral immunity and lipid transport at later stages. These effects of starvation on the host transcriptome and microbiome were almost completely restored within 3 days after refeeding. Comparison with published datasets identified host genes responsive to starvation as well as high-fat feeding or microbiome colonization, and predicted host transcription factors that may be involved in starvation response.

CONCLUSIONS

Long-term starvation induces progressive changes in microbiome composition and host gene expression in the zebrafish intestine, and these changes are rapidly reversed after refeeding. Our identification of bacterial taxa, host genes and host pathways involved in this response provides a framework for future investigation of the physiological and ecological mechanisms underlying intestinal adaptations to food restriction.

Gain of gene regulatory network interconnectivity at the origin of vertebrates

SignificanceIn this manuscript, we address an essential question in developmental and evolutionary biology: How have changes in gene regulatory networks contributed to the invertebrate-to-vertebrate transition? To address this issue, we perturbed four signaling pathways critical for body plan formation in the cephalochordate amphioxus and in zebrafish and compared the effects of such perturbations on gene expression and gene regulation in both species. Our data reveal that many developmental genes have gained response to these signaling pathways in the vertebrate lineage. Moreover, we show that the interconnectivity between these pathways is much higher in zebrafish than in amphioxus. We conclude that this increased signaling pathway complexity likely contributed to vertebrate morphological novelties during evolution.

Transposable element activation promotes neurodegeneration in a Drosophila model of Huntington's disease

Huntington's disease (HD) is an autosomal dominant disorder with progressive motor dysfunction and cognitive decline. The disease is caused by a CAG repeat expansion in the IT15 gene, which elongates a polyglutamine stretch of the HD protein, Huntingtin. No therapeutic treatments are available, and new pharmacological targets are needed. Retrotransposons are transposable elements (TEs) that represent 40% and 30% of the human and Drosophila genomes and replicate through an RNA intermediate. Mounting evidence suggests that mammalian TEs are active during neurogenesis and may be involved in diseases of the nervous system. Here we show that TE expression and mobilization are increased in a Drosophila melanogaster HD model. By inhibiting TE mobilization with Reverse Transcriptase inhibitors, polyQ-dependent eye neurodegeneration and genome instability in larval brains are rescued and fly lifespan is increased. These results suggest that TE activation may be involved in polyQ-induced neurotoxicity and a potential pharmacological target.

Decreased Expression of Estrogen Receptors Is Associated with Tumorigenesis in Papillary Thyroid Carcinoma

Papillary thyroid carcinomas (PTC), which is derived from thyroid follicular cells, is the most commonly differentiated thyroid cancer with sex disparity. However, the role of estrogen receptors (ERs) in the pathogenesis of PTC remains unclear. The present study aimed to determine the association of ER mRNA expression levels with clinicopathologic features in PTC. To that aim, the mRNA levels of ESR1 (ERα66), ESR1 (ERα36), ESR2, and G-protein-coupled estrogen receptor 1 (GPER1) in snap-frozen tissue samples from PTCs and adjacent normal thyroid tissues were determined using quantitative reverse transcription polymerase chain reaction (RT-qPCR), and the correlation between ER mRNA expression levels and clinicopathologic features was analyzed. The expression of ERα66, ERα36, ERβ, and GPER1 was lower in PTC specimens than in adjacent normal thyroid tissues. Moreover, low GPER1 expression was associated with extrathyroidal extension. There was no obvious difference in expression of ERs between PTC specimens from male and female patients. In conclusion, our findings highlight the importance of ERs in PTC tumorigenesis.

Identification of Germline Non-coding Deletions in XIAP Gene Causing XIAP Deficiency Reveals a Key Promoter Sequence

PURPOSE

X-linked inhibitor of apoptosis protein (XIAP) deficiency, also known as the X-linked lymphoproliferative syndrome of type 2 (XLP-2), is a rare immunodeficiency characterized by recurrent hemophagocytic lymphohistiocytosis, splenomegaly, and inflammatory bowel disease. Variants in XIAP including missense, non-sense, frameshift, and deletions of coding exons have been reported to cause XIAP deficiency. We studied three young boys with immunodeficiency displaying XLP-2-like clinical features. No genetic variation in the coding exons of XIAP was identified by whole-exome sequencing (WES), although the patients exhibited a complete loss of XIAP expression.

METHODS

Targeted next-generation sequencing (NGS) of the entire locus of XIAP was performed on DNA samples from the three patients. Molecular investigations were assessed by gene reporter expression assays in HEK cells and CRISPR-Cas9 genome editing in primary T cells.

RESULTS

NGS of XIAP identified three distinct non-coding deletions in the patients that were predicted to be driven by repetitive DNA sequences. These deletions share a common region of 839 bp that encompassed the first non-coding exon of XIAP and contained regulatory elements and marks specific of an active promoter. Moreover, we showed that among the 839 bp, the exon was transcriptionally active. Finally, deletion of the exon by CRISPR-Cas9 in primary cells reduced XIAP protein expression.

CONCLUSIONS

These results identify a key promoter sequence contained in the first non-coding exon of XIAP. Importantly, this study highlights that sequencing of the non-coding exons that are not currently captured by WES should be considered in the genetic diagnosis when no variation is found in coding exons.