TERRA ONTseq: a long-read-based sequencing pipeline to study the human telomeric transcriptome.

The long noncoding RNA TERRA is transcribed from telomeres in virtually all eukaryotes with linear chromosomes. In humans, TERRA transcription is driven in part by promoters comprising CpG dinucleotide-rich repeats of 29 bp repeats, believed to be present in half of the subtelomeres. Thus far, TERRA expression has been analyzed mainly using molecular biology-based approaches that only generate partial and somehow biased results. Here, we present a novel experimental pipeline to study human TERRA based on long-read sequencing (TERRA ONTseq). By applying TERRA ONTseq to different cell lines, we show that the vast majority of human telomeres produce TERRA and that the cellular levels of TERRA transcripts vary according to their chromosomes of origin. Using TERRA ONTseq, we also identified regions containing TERRA transcription start sites (TSSs) in more than half of human subtelomeres. TERRA TSS regions are generally found immediately downstream from 29 bp repeat-related sequences, which appear to be more widespread than previously estimated. Finally, we isolated a novel TERRA promoter from the highly expressed subtelomere of the long arm of Chromosome 7. With the development of TERRA ONTseq, we provide a refined picture of human TERRA biogenesis and expression and we equip the scientific community with an invaluable tool for future studies.

Reduced levels of prostaglandin I2 synthase: a distinctive feature of the cancer-free trichothiodystrophy.

The cancer-free photosensitive trichothiodystrophy (PS-TTD) and the cancer-prone xeroderma pigmentosum (XP) are rare monogenic disorders that can arise from mutations in the same genes, namely ERCC2/XPD or ERCC3/XPB Both XPD and XPB proteins belong to the 10-subunit complex transcription factor IIH (TFIIH) that plays a key role in transcription and nucleotide excision repair, the DNA repair pathway devoted to the removal of ultraviolet-induced DNA lesions. Compelling evidence suggests that mutations affecting the DNA repair activity of TFIIH are responsible for the pathological features of XP, whereas those also impairing transcription give rise to TTD. By adopting a relatives-based whole transcriptome sequencing approach followed by specific gene expression profiling in primary fibroblasts from a large cohort of TTD or XP cases with mutations in ERCC2/XPD gene, we identify the expression alterations specific for TTD primary dermal fibroblasts. While most of these transcription deregulations do not impact on the protein level, very low amounts of prostaglandin I2 synthase (PTGIS) are found in TTD cells. PTGIS catalyzes the last step of prostaglandin I2 synthesis, a potent vasodilator and inhibitor of platelet aggregation. Its reduction characterizes all TTD cases so far investigated, both the PS-TTD with mutations in TFIIH coding genes as well as the nonphotosensitive (NPS)-TTD. A severe impairment of TFIIH and RNA polymerase II recruitment on the PTGIS promoter is found in TTD but not in XP cells. Thus, PTGIS represents a biomarker that combines all PS- and NPS-TTD cases and distinguishes them from XP.

Protein instability associated with AARS1 and MARS1 mutations causes trichothiodystrophy.

Trichothiodystrophy (TTD) is a rare hereditary neurodevelopmental disorder defined by sulfur-deficient brittle hair and nails and scaly skin, but with otherwise remarkably variable clinical features. The photosensitive TTD (PS-TTD) forms exhibits in addition to progressive neuropathy and other features of segmental accelerated aging and is associated with impaired genome maintenance and transcription. New factors involved in various steps of gene expression have been identified for the different non-photosensitive forms of TTD (NPS-TTD), which do not appear to show features of premature aging. Here, we identify alanyl-tRNA synthetase 1 and methionyl-tRNA synthetase 1 variants as new gene defects that cause NPS-TTD. These variants result in the instability of the respective gene products alanyl- and methionyl-tRNA synthetase. These findings extend our previous observations that TTD mutations affect the stability of the corresponding proteins and emphasize this phenomenon as a common feature of TTD. Functional studies in skin fibroblasts from affected individuals demonstrate that these new variants also impact on the rate of tRNA charging, which is the first step in protein translation. The extension of reduced abundance of TTD factors to translation as well as transcription redefines TTD as a syndrome in which proteins involved in gene expression are unstable.

Identification of bi-allelic LFNG variants in three patients and further clinical and molecular refinement of spondylocostal dysostosis 3.

Spondylocostal dysostosis (SCD), a condition characterized by multiple segmentation defects of the vertebrae and rib malformations, is caused by bi-allelic variants in one of the genes involved in the Notch signaling pathway that tunes the "segmentation clock" of somitogenesis: DLL3, HES7, LFNG, MESP2, RIPPLY2, and TBX6. To date, seven individuals with LFNG variants have been reported in the literature. In this study we describe two newborns and one fetus with SCD, who were found by trio-based exome sequencing (trio-ES) to carry homozygous (c.822-5C>T) or compound heterozygous (c.[863dup];[1063G>A]) and (c.[521G>T];[890T>G]) variants in LFNG. Notably, the c.822-5C>T change, affecting the polypyrimidine tract of intron 5, is the first non-coding variant reported in LFNG. This study further refines the clinical and molecular features of spondylocostal dysostosis 3 and adds to the numerous investigations supporting the usefulness of trio-ES approach in prenatal and neonatal settings.

Alternative Splicing Changes Promoted by NOVA2 Upregulation in Endothelial Cells and Relevance for Gastric Cancer.

Angiogenesis is crucial for cancer progression. While several anti-angiogenic drugs are in use for cancer treatment, their clinical benefits are unsatisfactory. Thus, a deeper understanding of the mechanisms sustaining cancer vessel growth is fundamental to identify novel biomarkers and therapeutic targets. Alternative splicing (AS) is an essential modifier of human proteome diversity. Nevertheless, AS contribution to tumor vasculature development is poorly known. The Neuro-Oncological Ventral Antigen 2 (NOVA2) is a critical AS regulator of angiogenesis and vascular development. NOVA2 is upregulated in tumor endothelial cells (ECs) of different cancers, thus representing a potential driver of tumor blood vessel aberrancies. Here, we identified novel AS transcripts generated upon NOVA2 upregulation in ECs, suggesting a pervasive role of NOVA2 in vascular biology. In addition, we report that NOVA2 is also upregulated in ECs of gastric cancer (GC), and its expression correlates with poor overall survival of GC patients. Finally, we found that the AS of the Rap Guanine Nucleotide Exchange Factor 6 (RapGEF6), a newly identified NOVA2 target, is altered in GC patients and associated with NOVA2 expression, tumor angiogenesis, and poor patient outcome. Our findings provide a better understanding of GC biology and suggest that AS might be exploited to identify novel biomarkers and therapeutics for anti-angiogenic GC treatments.

Bi-allelic TARS Mutations Are Associated with Brittle Hair Phenotype.

Brittle and "tiger-tail" hair is the diagnostic hallmark of trichothiodystrophy (TTD), a rare recessive disease associated with a wide spectrum of clinical features including ichthyosis, intellectual disability, decreased fertility, and short stature. As a result of premature abrogation of terminal differentiation, the hair is brittle and fragile and contains reduced cysteine content. Hypersensitivity to UV light is found in about half of individuals with TTD; all of these individuals harbor bi-allelic mutations in components of the basal transcription factor TFIIH, and these mutations lead to impaired nucleotide excision repair and basal transcription. Different genes have been found to be associated with non-photosensitive TTD (NPS-TTD); these include MPLKIP (also called TTDN1), GTF2E2 (also called TFIIEß), and RNF113A. However, a relatively large group of these individuals with NPS-TTD have remained genetically uncharacterized. Here we present the identification of an NPS-TTD-associated gene, threonyl-tRNA synthetase (TARS), found by next-generation sequencing of a group of uncharacterized individuals with NPS-TTD. One individual has compound heterozygous TARS variants, c.826A>G (p.Lys276Glu) and c.1912C>T (p.Arg638∗), whereas a second individual is homozygous for the TARS variant: c.680T>C (p.Leu227Pro). We showed that these variants have a profound effect on TARS protein stability and enzymatic function. Our results expand the spectrum of genes involved in TTD to include genes implicated in amino acid charging of tRNA, which is required for the last step in gene expression, namely protein translation. We previously proposed that some of the TTD-specific features derive from subtle transcription defects as a consequence of unstable transcription factors. We now extend the definition of TTD from a transcription syndrome to a "gene-expression" syndrome.

A Role for Human DNA Polymerase λ in Alternative Lengthening of Telomeres.

Telomerase negative cancer cell types use the Alternative Lengthening of Telomeres (ALT) pathway to elongate telomeres ends. Here, we show that silencing human DNA polymerase (Pol λ) in ALT cells represses ALT activity and induces telomeric stress. In addition, replication stress in the absence of Pol λ, strongly affects the survival of ALT cells. In vitro, Pol λ can promote annealing of even a single G-rich telomeric repeat to its complementary strand and use it to prime DNA synthesis. The noncoding telomeric repeat containing RNA TERRA and replication protein A negatively regulate this activity, while the Protection of Telomeres protein 1 (POT1)/TPP1 heterodimer stimulates Pol λ. Pol λ associates with telomeres and colocalizes with TPP1 in cells. In summary, our data suggest a role of Pol λ in the maintenance of telomeres by the ALT mechanism.

Characterization of the biological processes shaping the genetic structure of the Italian population.

BACKGROUND: The genetic structure of human populations is the outcome of the combined action of different processes such as demographic dynamics and natural selection. Several efforts toward the characterization of population genetic architectures and the identification of adaptation signatures were recently made. In this study, we provide a genome-wide depiction of the Italian population structure and the analysis of the major determinants of the current existing genetic variation. RESULTS: We defined and characterized 210 genomic loci associated with the first Principal Component calculated on the Italian genotypic data and correlated to the North-south genetic gradient. Using a gene-enrichment approach we identified the immune function as primarily involved in the Italian population differentiation and we described a locus on chromosome 13 showing combined evidence of North-south diversification in allele frequencies and signs of recent positive selection. In this region our bioinformatics analysis pinpointed an uncharacterized long intergenic non-coding (lincRNA), whose expression appeared specific for immune-related tissues suggesting its relevance for the immune function. CONCLUSIONS: Our study, combining population genetic analyses with biological insights provides a description of the Italian genetic structure that in future could contribute to the evaluation of complex diseases risk in the population context.

CorrelaGenes: a new tool for the interpretation of the human transcriptome.

BACKGROUND: The amount of gene expression data available in public repositories has grown exponentially in the last years, now requiring new data mining tools to transform them in information easily accessible to biologists. RESULTS: By exploiting expression data publicly available in the Gene Expression Omnibus (GEO) database, we developed a new bioinformatics tool aimed at the identification of genes whose expression appeared simultaneously altered in different experimental conditions, thus suggesting co-regulation or coordinated action in the same biological process. To accomplish this task, we used the 978 human GEO Curated DataSets and we manually performed the selection of 2,109 pair-wise comparisons based on their biological rationale. The lists of differentially expressed genes, obtained from the selected comparisons, were stored in a PostgreSQL database and used as data source for the CorrelaGenes tool. Our application uses a customized Association Rule Mining (ARM) algorithm to identify sets of genes showing expression profiles correlated with a gene of interest. The significance of the correlation is measured coupling the Lift, a well-known standard ARM index, and the χ(2) p value. The manually curated selection of the comparisons and the developed algorithm constitute a new approach in the field of gene expression profiling studies. Simulation performed on 100 randomly selected target genes allowed us to evaluate the efficiency of the procedure and to obtain preliminary data demonstrating the consistency of the results. CONCLUSIONS: The preliminary results of the simulation showed how CorrelaGenes could contribute to the characterization of molecular pathways and biological processes integrating data obtained from other applications and available in public repositories.

Are myocardial infarction--associated single-nucleotide polymorphisms associated with ischemic stroke?

BACKGROUND AND PURPOSE: Ischemic stroke (IS) shares many common risk factors with coronary artery disease (CAD). We hypothesized that genetic variants associated with myocardial infarction (MI) or CAD may be similarly involved in the etiology of IS. To test this hypothesis, we evaluated whether single-nucleotide polymorphisms (SNPs) at 11 different loci recently associated with MI or CAD through genome-wide association studies were associated with IS. METHODS: Meta-analyses of the associations between the 11 MI-associated SNPs and IS were performed using 6865 cases and 11 395 control subjects recruited from 9 studies. SNPs were either genotyped directly or imputed; in a few cases a surrogate SNP in high linkage disequilibrium was chosen. Logistic regression was performed within each study to obtain study-specific ßs and standard errors. Meta-analysis was conducted using an inverse variance weighted approach assuming a random effect model. RESULTS: Despite having power to detect odds ratio of 1.09-1.14 for overall IS and 1.20-1.32 for major stroke subtypes, none of the SNPs were significantly associated with overall IS and/or stroke subtypes after adjusting for multiple comparisons. CONCLUSIONS: Our results suggest that the major common loci associated with MI risk do not have effects of similar magnitude on overall IS but do not preclude moderate associations restricted to specific IS subtypes. Disparate mechanisms may be critical in the development of acute ischemic coronary and cerebrovascular events.

DEAD-Box RNA Helicases DDX3X and DDX5 as Oncogenes or Oncosuppressors: A Network Perspective.

RNA helicases of the DEAD-box family are involved in several metabolic pathways, from transcription and translation to cell proliferation, innate immunity and stress response. Given their multiple roles, it is not surprising that their deregulation or mutation is linked to different pathological conditions, including cancer. However, while in some cases the loss of function of a given DEAD-box helicase promotes tumor transformation, indicating an oncosuppressive role, in other contexts the overexpression of the same enzyme favors cancer progression, thus acting as a typical oncogene. The roles of two well-characterized members of this family, DDX3X and DDX5, as both oncogenes and oncosuppressors have been documented in several cancer types. Understanding the interplay of the different cellular contexts, as defined by the molecular interaction networks of DDX3X and DDX5 in different tumors, with the cancer-specific roles played by these proteins could help to explain their apparently conflicting roles as cancer drivers or suppressors.

TERRA transcription destabilizes telomere integrity to initiate break-induced replication in human ALT cells.

Alternative Lengthening of Telomeres (ALT) is a Break-Induced Replication (BIR)-based mechanism elongating telomeres in a subset of human cancer cells. While the notion that spontaneous DNA damage at telomeres is required to initiate ALT, the molecular triggers of this physiological telomere instability are largely unknown. We previously proposed that the telomeric long noncoding RNA TERRA may represent one such trigger; however, given the lack of tools to suppress TERRA transcription in cells, our hypothesis remained speculative. We have developed Transcription Activator-Like Effectors able to rapidly inhibit TERRA transcription from multiple chromosome ends in an ALT cell line. TERRA transcription inhibition decreases marks of DNA replication stress and DNA damage at telomeres and impairs ALT activity and telomere length maintenance. We conclude that TERRA transcription actively destabilizes telomere integrity in ALT cells, thereby triggering BIR and promoting telomere elongation. Our data point to TERRA transcription manipulation as a potentially useful target for therapy.

Alterations in the expression, structure and function of progesterone receptor membrane component-1 (PGRMC1) in premature ovarian failure.

Premature ovarian failure (POF) is characterized by hypergonadotropic hypogonadism and amenorrhea before the age of 40. The condition has a heterogeneous background but genetic factors are demonstrated by the occurrence of familial cases. We identified a mother and daughter with POF both of whom carry an X;autosome translocation [t(X;11)(q24;q13)]. RNA expression studies of genes flanking the X-chromosome breakpoint revealed that both patients have reduced expression levels of the gene Progesterone Receptor Membrane Component-1 (PGRMC1). Mutation screening of 67 females with idiopathic POF identified a third patient with a missense mutation (H165R) located in the cytochrome b5 domain of PGRMC1. PGRMC1 mediates the anti-apoptotic action of progesterone in ovarian cells and it acts as a positive regulator of several cytochrome P450 (CYP)-catalyzed reactions. The CYPs are critical for intracellular sterol metabolism, including biosynthesis of steroid hormones. We show that the H165R mutation associated with POF abolishes the binding of cytochrome P450 7A1 (CYP7A1) to PGRMC1. In addition, the missense mutation attenuates PGRMC1's ability to mediate the anti-apoptotic action of progesterone in ovarian cells. These findings suggest that mutant or reduced levels of PGMRC1 may cause POF through impaired activation of the microsomal cytochrome P450 and increased apoptosis of ovarian cells.

GC-AG Introns Features in Long Non-coding and Protein-Coding Genes Suggest Their Role in Gene Expression Regulation.

Long non-coding RNAs (lncRNAs) are recognized as an important class of regulatory molecules involved in a variety of biological functions. However, the regulatory mechanisms of long non-coding genes expression are still poorly understood. The characterization of the genomic features of lncRNAs is crucial to get insight into their function. In this study, we exploited recent annotations by GENCODE to characterize the genomic and splicing features of long non-coding genes in comparison with protein-coding ones, both in human and mouse. Our analysis highlighted differences between the two classes of genes in terms of their gene architecture. Significant differences in the splice sites usage were observed between long non-coding and protein-coding genes (PCG). While the frequency of non-canonical GC-AG splice junctions represents about 0.8% of total splice sites in PCGs, we identified a significant enrichment of the GC-AG splice sites in long non-coding genes, both in human (3.0%) and mouse (1.9%). In addition, we found a positional bias of GC-AG splice sites being enriched in the first intron in both classes of genes. Moreover, a significant shorter length and weaker donor and acceptor sites were found comparing GC-AG introns to GT-AG introns. Genes containing at least one GC-AG intron were found conserved in many species, more prone to alternative splicing and a functional analysis pointed toward their enrichment in specific biological processes such as DNA repair. Our study shows for the first time that GC-AG introns are mainly associated with lncRNAs and are preferentially located in the first intron. Additionally, we discovered their regulatory potential indicating the existence of a new mechanism of non-coding and PCGs expression regulation.

BMP15 mutations associated with primary ovarian insufficiency cause a defective production of bioactive protein.

Bone morphogenetic protein-15 (BMP15) is selectively synthesized by oocytes as a pre-proprotein and is considered an ovarian follicle organizer whose adequate function is critical for female fertility. Missense mutations were reported in primary ovarian insufficiency (POI) but their biological impact remained unexplored. Here, screening of 300 unrelated idiopathic overt POI women with primary or secondary amenorrhea (SA) led to the identification of six heterozygous BMP15 variations in 29 of them. All alterations are nonconservative and include one insertion of three nucleotides (p.L262_L263insL) and five missense substitutions. Except for the p.S5R located in the signal sequence, the other variants (p.R68W, p.R138H, p.L148P, and p.A180T) localize in the proregion, which is essential for the processing and secretion of bioactive dimers. The mutations p.R68W, p.L148P, and the novel p.R138H lead to marked reductions of mature protein production. Their biological effects, evaluated by a novel luciferase-reporter assay in a human granulosa cell (GC) line, were significantly reduced. Cotransfection experiments of defective mutants with equal amounts of wild-type BMP15 cDNA, thus reproducing the heterozygous state seen in patients, did not generate a complete recovery of wild-type activity. No or minor deleterious effects were detected for the variants p.L262_L263insL, p.A180T, or p.S5R. In conclusion, heterozygous BMP15 mutations associated with the early onset of overt POI lead to defective secretion of bioactive dimers. These findings support the concept that an adequate amount of BMP15 secreted in the follicular fluid is critical for female fertility. We propose to consider the screening of BMP15 mutations among the analyses for the prediction of POI risk.

A large-scale association study to assess the impact of known variants of the human INHA gene on premature ovarian failure.

BACKGROUND: Three variants of the human INHA gene have been reported to be associated with premature ovarian failure (POF) in case-control studies involving a small number of patients and controls. Since inhibin has a fundamental role in the control of ovarian function, it is important to establish the relevance of the reported variants for disease risk. METHODS: Three independent POF cohorts, recruited in Northern and Central Italy and in Germany consisting of a total of 611 patients and 1084 matched controls, were genotyped for the three variants: -16C > T, -124A > G and 769G > A. RESULTS: No significant difference was detected between allelic frequencies of the INHA promoter variants between POF patients and controls. The rare allele in the coding variant appeared to be more frequent among the control populations. CONCLUSIONS: The association between the INHA promoter variants and POF could not be replicated, and our results suggest that this discrepancy is likely to be due to the small sample size of previous studies. The rare allele of the coding variant seems to exert a protective effect against loss of ovarian function, which should be confirmed in additional large and ethnically diverse cohorts.

Identification of a novel antiviral micro-RNA targeting the NS1 protein of the H1N1 pandemic human influenza virus and a corresponding viral escape mutation.

The influenza A virus (IAV) NS1 protein is one of the major regulators of pathogenicity, being able to suppress innate immune response and host protein synthesis. In this study we identified the human micro RNA hsa-miR-1307-3p as a novel potent suppressor of NS1 expression and influenza virus replication. Transcriptomic analysis indicates that hsa-miR-1307-3p also negatively regulates apoptosis and promotes cell proliferation. In addition, we identified a novel mutation in the NS1 gene of A(H1N1)pdm09 strains circulating in Italy in the 2010-11 season, which enabled the virus to escape the hsa-miR-1307-3p inhibition, conferring replicative advantage to the virus in human cells. To the best of our knowledge, this is the first validation of suppression of IAV H1N1 NS1 by a human micro RNA and the first example of an escape mutation from micro RNA-mediated antiviral response for the A(H1N1)pdm09 virus.

Sequence variation at the human FOXO3 locus: a study of premature ovarian failure and primary amenorrhea.

BACKGROUND: The forkhead transcription factor Foxo3 is a master regulator and potent suppressor of primordial follicle activation. Loss of Foxo3 function in the mouse leads to premature ovarian failure (POF) due to global follicle activation. METHODS AND RESULTS: Here, we show that the mouse Foxo3 locus is haploinsufficient, and that Foxo3-/+ females undergo early reproductive senescence consistent with an increased rate of primordial follicle utilization. Then, to determine if heterozygous or homozygous polymorphisms or mutations of the human orthologue FOXO3 contribute to POF or idiopathic primary amenorrhea (PA), we sequenced the exons and flanking splice sequences of the gene in a large number of women with idiopathic POF (n = 273) or PA (n = 29). A total of eight single-nucleotide polymorphisms (SNPs) were identified, revealing a substantial amount of genetic variation at this locus. Allelic frequencies in control samples excluded several of these variants as causal. For the remaining variants, site-directed mutagenesis was performed to assess their functional impact. However, these rare sequence variants were not associated with significant decreases in FOXO3 activity. CONCLUSIONS: Taken together, our findings suggest that, despite the potential for FOXO3 haploinsufficiency to cause ovarian failure, FOXO3 mutations or common SNPs are not a common cause of either POF or PA.

Variation of hemoglobin levels in normal Italian populations from genetic isolates.

Normal hemoglobin levels vary greatly according to genetic and acquired factors. As a consequence there is no general agreement on the definition of anemia in terms of hemoglobin levels. Here we compare the hemoglobin levels of subjects recruited from normal genetically isolated Italian populations whose medical history, life style habits and results of laboratory tests are available. After the exclusion of pathological samples we analyzed the hemoglobin levels of 3,849 subjects (1,661 males and 2,188 females) and evaluated the hemoglobin heritability. Normal subjects of different age groups from a northern Italian isolate have significantly higher hemoglobin levels when compared to matched subjects of southern Italian isolates. The estimated heritability of hemoglobin levels ranges from 0.34 to 0.42 in the different isolates. Our study provides a dataset of hemoglobin levels for normal subjects of different geographical origin and indicate that hemoglobin levels are substantially influenced by heritable components.