The RNA binding protein quaking regulates formation of circRNAs.

Circular RNAs (circRNAs), formed by non-sequential back-splicing of pre-mRNA transcripts, are a widespread form of non-coding RNA in animal cells. However, it is unclear whether the majority of circRNAs represent splicing by-products without function or are produced in a regulated manner to carry out specific cellular functions. We show that hundreds of circRNAs are regulated during human epithelial-mesenchymal transition (EMT) and find that the production of over one-third of abundant circRNAs is dynamically regulated by the alternative splicing factor, Quaking (QKI), which itself is regulated during EMT. Furthermore, by modulating QKI levels, we show the effect on circRNA abundance is dependent on intronic QKI binding motifs. Critically, the addition of QKI motifs is sufficient to induce de novo circRNA formation from transcripts that are normally linearly spliced. These findings demonstrate circRNAs are both purposefully synthesized and regulated by cell-type specific mechanisms, suggesting they play specific biological roles in EMT.

A Prox1 enhancer represses haematopoiesis in the lymphatic vasculature.

Transcriptional enhancer elements are responsible for orchestrating the temporal and spatial control over gene expression that is crucial for programming cell identity during development1-3. Here we describe a novel enhancer element that is important for regulating the expression of Prox1 in lymphatic endothelial cells. This evolutionarily conserved enhancer is bound by key lymphatic transcriptional regulators including GATA2, FOXC2, NFATC1 and PROX1. Genome editing of the enhancer to remove five nucleotides encompassing the GATA2-binding site resulted in perinatal death of homozygous mutant mice due to profound lymphatic vascular defects. Lymphatic endothelial cells in enhancer mutant mice exhibited reduced expression of genes characteristic of lymphatic endothelial cell identity and increased expression of genes characteristic of haemogenic endothelium, and acquired the capacity to generate haematopoietic cells. These data not only reveal a transcriptional enhancer element important for regulating Prox1 expression and lymphatic endothelial cell identity but also demonstrate that the lymphatic endothelium has haemogenic capacity, ordinarily repressed by Prox1.

Impact of additional genetic abnormalities at diagnosis of chronic myeloid leukemia for first-line imatinib-treated patients receiving proactive treatment intervention.

The BCR::ABL1 gene fusion initiates chronic myeloid leukemia (CML); however, evidence has accumulated from studies of highly selected cohorts that variants in other cancer-related genes are associated with treatment failure. Nevertheless, the true incidence and impact of additional genetic abnormalities (AGA) at diagnosis of chronic phase (CP)-CML is unknown. We sought to determine whether AGA at diagnosis in a consecutive imatinib-treated cohort of 210 patients enrolled in the TIDEL-II trial influenced outcome despite a highly proactive treatment intervention strategy. Survival outcomes including overall survival, progression-free survival, failure-free survival, and BCR::ABL1 kinase domain mutation acquisition were evaluated. Molecular outcomes were measured at a central laboratory and included major molecular response (MMR, BCR::ABL1 ≤0.1%IS), MR4 (BCR::ABL1 ≤0.01%IS), and MR4.5 (BCR::ABL1 ≤0.0032%IS). AGA included variants in known cancer genes and novel rearrangements involving the formation of the Philadelphia chromosome. Clinical outcomes and molecular response were assessed based on the patient's genetic profile and other baseline factors. AGA were identified in 31% of patients. Potentially pathogenic variants in cancer-related genes were detected in 16% of patients at diagnosis (including gene fusions and deletions) and structural rearrangements involving the Philadelphia chromosome (Ph-associated rearrangements) were detected in 18%. Multivariable analysis demonstrated that the combined genetic abnormalities plus the EUTOS long-term survival clinical risk score were independent predictors of lower molecular response rates and higher treatment failure. Despite a highly proactive treatment intervention strategy, first-line imatinib-treated patients with AGA had poorer response rates. These data provide evidence for the incorporation of genomically-based risk assessment for CML.

miR-200/375 control epithelial plasticity-associated alternative splicing by repressing the RNA-binding protein Quaking.

Members of the miR-200 family are critical gatekeepers of the epithelial state, restraining expression of pro-mesenchymal genes that drive epithelial-mesenchymal transition (EMT) and contribute to metastatic cancer progression. Here, we show that miR-200c and another epithelial-enriched miRNA, miR-375, exert widespread control of alternative splicing in cancer cells by suppressing the RNA-binding protein Quaking (QKI). During EMT, QKI-5 directly binds to and regulates hundreds of alternative splicing targets and exerts pleiotropic effects, such as increasing cell migration and invasion and restraining tumour growth, without appreciably affecting mRNA levels. QKI-5 is both necessary and sufficient to direct EMT-associated alternative splicing changes, and this splicing signature is broadly conserved across many epithelial-derived cancer types. Importantly, several actin cytoskeleton-associated genes are directly targeted by both QKI and miR-200c, revealing coordinated control of alternative splicing and mRNA abundance during EMT These findings demonstrate the existence of a miR-200/miR-375/QKI axis that impacts cancer-associated epithelial cell plasticity through widespread control of alternative splicing.

Pseudodiastrophic dysplasia expands the known phenotypic spectrum of defects in proteoglycan biosynthesis.

BACKGROUND: Pseudodiastrophic dysplasia (PDD) is a severe skeletal dysplasia associated with prenatal manifestation and early lethality. Clinically, PDD is classified as a 'dysplasia with multiple joint dislocations'; however, the molecular aetiology of the disorder is currently unknown. METHODS: Whole exome sequencing (WES) was performed on three patients from two unrelated families, clinically diagnosed with PDD, in order to identify the underlying genetic cause. The functional effects of the identified variants were characterised using primary cells and human cell-based overexpression assays. RESULTS: WES resulted in the identification of biallelic variants in the established skeletal dysplasia genes, B3GAT3 (family 1) and CANT1 (family 2). Mutations in these genes have previously been reported to cause 'multiple joint dislocations, short stature, and craniofacial dysmorphism with or without congenital heart defects' ('JDSCD'; B3GAT3) and Desbuquois dysplasia 1 (CANT1), disorders in the same nosological group as PDD. Follow-up of the B3GAT3 variants demonstrated significantly reduced B3GAT3/GlcAT-I expression. Downstream in vitro functional analysis revealed abolished biosynthesis of glycosaminoglycan side chains on proteoglycans. Functional evaluation of the CANT1 variant showed impaired nucleotidase activity, which results in inhibition of glycosaminoglycan synthesis through accumulation of uridine diphosphate. CONCLUSION: For the families described in this study, the PDD phenotype was caused by mutations in the known skeletal dysplasia genes B3GAT3 and CANT1, demonstrating the advantage of genomic analyses in delineating the molecular diagnosis of skeletal dysplasias. This finding expands the phenotypic spectrum of B3GAT3-related and CANT1-related skeletal dysplasias to include PDD and highlights the significant phenotypic overlap of conditions within the proteoglycan biosynthesis pathway.

Two monogenic disorders masquerading as one: severe congenital neutropenia with monocytosis and non-syndromic sensorineural hearing loss.

BACKGROUND: We report a large family with four successive generations, presenting with a complex phenotype of severe congenital neutropenia (SCN), partially penetrant monocytosis, and hearing loss of varying severity. METHODS: We performed whole exome sequencing to identify the causative variants. Sanger sequencing was used to perform segregation analyses on remaining family members. RESULTS: We identified and classified a pathogenic GFI1 variant and a likely pathogenic variant in MYO6 which together explain the complex phenotypes seen in this family. CONCLUSIONS: We present a case illustrating the benefits of a broad screening approach that allows identification of oligogenic determinants of complex human phenotypes which may have been missed if the screening was limited to a targeted gene panel with the assumption of a syndromic disorder. This is important for correct genetic diagnosis of families and disentangling the range and severity of phenotypes associated with high impact variants.

Integrative genomic analysis reveals cancer-associated mutations at diagnosis of CML in patients with high-risk disease.

Genomic events associated with poor outcome in chronic myeloid leukemia (CML) are poorly understood. We performed whole-exome sequencing, copy-number variation, and/or RNA sequencing for 65 patients to discover mutations at diagnosis and blast crisis (BC). Forty-six patients with chronic-phase disease with the extremes of outcome were studied at diagnosis. Cancer gene variants were detected in 15 (56%) of 27 patients with subsequent BC or poor outcome and in 3 (16%) of 19 optimal responders (P = .007). Frequently mutated genes at diagnosis were ASXL1, IKZF1, and RUNX1 The methyltransferase SETD1B was a novel recurrently mutated gene. A novel class of variant associated with the Philadelphia (Ph) translocation was detected at diagnosis in 11 (24%) of 46 patients comprising fusions and/or rearrangement of genes on the translocated chromosomes, with evidence of fragmentation, inversion, and imperfect sequence reassembly. These were more frequent at diagnosis in patients with poor outcome: 9 (33%) of 27 vs 2 (11%) of 19 optimal responders (P = .07). Thirty-nine patients were tested at BC, and all had cancer gene variants, including ABL1 kinase domain mutations in 58%. However, ABL1 mutations cooccurred with other mutated cancer genes in 89% of cases, and these predated ABL1 mutations in 62% of evaluable patients. Gene fusions not associated with the Ph translocation occurred in 42% of patients at BC and commonly involved fusion partners that were known cancer genes (78%). Genomic analysis revealed numerous relevant variants at diagnosis in patients with poor outcome and all patients at BC. Future refined biomarker testing of specific variants will likely provide prognostic information to facilitate a risk-adapted therapeutic approach.

Paternal mosaicism for a novel PBX1 mutation associated with recurrent perinatal death: Phenotypic expansion of the PBX1-related syndrome.

Autosomal dominant (de novo) mutations in PBX1 are known to cause congenital abnormalities of the kidney and urinary tract (CAKUT), with or without extra-renal abnormalities. Using trio exome sequencing, we identified a PBX1 p.(Arg107Trp) mutation in a deceased one-day-old neonate presenting with CAKUT, asplenia, and severe bilateral diaphragmatic thinning and eventration. Further investigation by droplet digital PCR revealed that the mutation had occurred post-zygotically in the father, with different variant allele frequencies of the mosaic PBX1 mutation in blood (10%) and sperm (20%). Interestingly, the father had subclinical hydronephrosis in childhood. With an expected recurrence risk of one in five, chorionic villus sampling and prenatal diagnosis for the PBX1 mutation identified recurrence in a subsequent pregnancy. The family opted to continue the pregnancy and the second affected sibling was stillborn at 35 weeks, presenting with similar severe bilateral diaphragmatic eventration, microsplenia, and complete sex reversal (46, XY female). This study highlights the importance of follow-up studies for presumed de novo and low-level mosaic variants and broadens the phenotypic spectrum of developmental abnormalities caused by PBX1 mutations.

A putative role for the aryl hydrocarbon receptor (AHR) gene in a patient with cyclical Cushing's disease.

BACKGROUND: Apart from PRKAR1A mutations in a subset of cyclical Cushing's syndrome due to primary pigmented nodular adrenocortical disease, the molecular basis of cyclical Cushing's syndrome has not been investigated. We speculated that cyclical Cushing's syndrome may be due to mutations in the clock genes that govern circadian rhythms, including the hypothalamic-pituitary-adrenal axis. CASE PRESENTATION: A 47-year-old man presented with mass effects from a sellar lesion. He was ultimately diagnosed with cyclical Cushing's disease due to a giant corticotrophinoma. We performed whole exome sequencing of germline and tumour DNA, SNP array of tumour DNA and tumour immunohistochemistry in order to detect variants in candidate circadian/pituitary-associated genes. We identified a rare germline missense variant in the aryl hydrocarbon receptor (AHR) gene, which has previously been indirectly linked to pituitary tumorigenesis and clock system disruption. The AHR variant was found in a highly conserved site involved in phosphorylation. It was predicted to be damaging by multiple in silico tools and AHR tumour immunohistochemistry demonstrated loss of the normal nuclear staining pattern, suggestive of an inactivating mutation. We also found a novel, damaging germline missense variant in the retinoid X receptor gamma (RXRG) gene, multiple somatic chromosomal gains (including AHR), and a somatic mutational signature consistent with oncogenesis that may have acted synergistically with the AHR variant. CONCLUSIONS: This is the first report of an AHR variant with predicted pathogenicity in the pituitary adenoma setting. Our preliminary data suggest that the highly conserved AHR gene may represent a link between pituitary tumorigenesis, the hypothalamic-pituitary-adrenal axis and the clock system. Further research may indicate a role for the gene in the development of cyclical Cushing's disease.

Novel germ line DDX41 mutations define families with a lower age of MDS/AML onset and lymphoid malignancies.

Recently our group and others have identified DDX41 mutations both as germ line and acquired somatic mutations in families with multiple cases of late onset myelodysplastic syndrome (MDS) and/or acute myeloid leukemia (AML), suggesting that DDX41 acts as a tumor suppressor. To determine whether novel DDX41 mutations could be identified in families with additional types of hematologic malignancies, our group screened two cohorts of families with a diverse range of hematologic malignancy subtypes. Among 289 families, we identified nine (3%) with DDX41 mutations. As previously observed, MDS and AML were the most common malignancies, often of the erythroblastic subtype, and 1 family displayed early-onset follicular lymphoma. Five novel mutations were identified, including missense mutations within important functional domains and start-loss and splicing mutations predicted to result in truncated proteins. We also show that most asymptomatic mutation carriers have normal blood counts until malignancy develops. This study expands both the mutation and phenotypic spectra observed in families with germ line DDX41 mutations. With an increasing number of both inherited and acquired mutations in this gene being identified, further study of how DDX41 disruption leads to hematologic malignancies is critical.

BAM-matcher: a tool for rapid NGS sample matching.

UNLABELLED: The standard method used by high-throughput genome sequencing facilities for detecting mislabelled samples is to use independently generated high-density SNP data to determine sample identity. However, as it has now become commonplace to have multiple samples sequenced from the same source, such as for analysis of somatic variants using matched tumour and normal samples, we can directly use the genotype information inherent in the sequence data to match samples and thus bypass the need for additional laboratory testing. Here we present BAM-matcher, a tool that can rapidly determine whether two BAM files represent samples from the same biological source by comparing their genotypes. BAM-matcher is designed to be simple to use, provides easily interpretable results, and is suitable for deployment at early stages of data processing pipelines. AVAILABILITY AND IMPLEMENTATION: BAM-matcher is licensed under the Creative Commons by Attribution license, and is available from: https://bitbucket.org/sacgf/bam-matcher SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online. CONTACT: paul.wang@sa.gov.au.

T cell receptor assessment in autoimmune disease requires access to the most adjacent immunologically active organ.

Next generation sequencing of T and B cell receptors is emerging as a valuable and effective method to diagnose and monitor hematopoietic malignancies. So far, this approach has not been fully explored in regard to autoimmune diseases. T cells develop in the thymus where they undergo positive and negative selection, and the autoimmune regulator (Aire) is central in the establishment of immunological tolerance. Loss of Aire leads to severe multiorgan autoimmune disease with infiltration of autoreactive T cells in affected organs. Here, we have utilized next generation sequencing technology to investigate the T cell receptor repertoire in autoimmunity induced by immunization of mice with a self-antigen, myeloperoxidase. By investigating the T cell receptor repertoire in peripheral blood, spleen and lumbar lymph nodes from naïve and immunized Aire -/- mice and wild type littermates, changes in the usage of V and J genes were evident. Our results identify TCR clonotypes which could be potential targets for immune therapy. Also, Aire -/- autoimmunity is driven by a variety of autoantigens where the autoimmune response is highly polyclonal, and access to the most adjacent immunologically active tissue is required to identify T cell receptor sequences that are potentially unique to the antigen in Aire-/- immunized mice.

Mutations in the mammalian target of rapamycin pathway regulators NPRL2 and NPRL3 cause focal epilepsy.

OBJECTIVE: Focal epilepsies are the most common form observed and have not generally been considered to be genetic in origin. Recently, we identified mutations in DEPDC5 as a cause of familial focal epilepsy. In this study, we investigated whether mutations in the mammalian target of rapamycin (mTOR) regulators, NPRL2 and NPRL3, also contribute to cases of focal epilepsy. METHODS: We used targeted capture and next-generation sequencing to analyze 404 unrelated probands with focal epilepsy. We performed exome sequencing on two families with multiple members affected with focal epilepsy and linkage analysis on one of these. RESULTS: In our cohort of 404 unrelated focal epilepsy patients, we identified five mutations in NPRL2 and five in NPRL3. Exome sequencing analysis of two families with focal epilepsy identified NPRL2 and NPRL3 as the top candidate-causative genes. Some patients had focal epilepsy associated with brain malformations. We also identified 18 new mutations in DEPDC5. INTERPRETATION: We have identified NPRL2 and NPRL3 as two new focal epilepsy genes that also play a role in the mTOR-signaling pathway. Our findings show that mutations in GATOR1 complex genes are the most significant cause of familial focal epilepsy identified to date, including cases with brain malformations. It is possible that deregulation of cellular growth control plays a more important role in epilepsy than is currently recognized.

Ectrodactyly and Lethal Pulmonary Acinar Dysplasia Associated with Homozygous FGFR2 Mutations Identified by Exome Sequencing.

Ectrodactyly/split hand-foot malformation is genetically heterogeneous with more than 100 syndromic associations. Acinar dysplasia is a rare congenital lung lesion of unknown etiology, which is frequently lethal postnatally. To date, there have been no reports of combinations of these two phenotypes. Here, we present an infant from a consanguineous union with both ectrodactyly and autopsy confirmed acinar dysplasia. SNP array and whole-exome sequencing analyses of the affected infant identified a novel homozygous Fibroblast Growth Factor Receptor 2 (FGFR2) missense mutation (p.R255Q) in the IgIII domain (D3). Expression studies of Fgfr2 in development show localization to the affected limbs and organs. Molecular modeling and genetic and functional assays support that this mutation is at least a partial loss-of-function mutation, and contributes to ectrodactyly and acinar dysplasia only in homozygosity, unlike previously reported heterozygous activating FGFR2 mutations that cause Crouzon, Apert, and Pfeiffer syndromes. This is the first report of mutations in a human disease with ectrodactyly with pulmonary acinar dysplasia and, as such, homozygous loss-of-function FGFR2 mutations represent a unique syndrome.

Maize maintains growth in response to decreased nitrate supply through a highly dynamic and developmental stage-specific transcriptional response.

Elucidation of the gene networks underlying the response to N supply and demand will facilitate the improvement of the N uptake efficiency of plants. We undertook a transcriptomic analysis of maize to identify genes responding to both a non-growth-limiting decrease in NO3- provision and to development-based N demand changes at seven representative points across the life cycle. Gene co-expression networks were derived by cluster analysis of the transcript profiles. The majority of NO3--responsive transcription occurred at 11 (D11), 18 (D18) and 29 (D29) days after emergence, with differential expression predominating in the root at D11 and D29 and in the leaf at D18. A cluster of 98 probe sets was identified, the expression pattern of which is similar to that of the high-affinity NO3- transporter (NRT2) genes across the life cycle. The cluster is enriched with genes encoding enzymes and proteins of lipid metabolism and transport, respectively. These are candidate genes for the response of maize to N supply and demand. Only a few patterns of differential gene expression were observed over the entire life cycle; however, the composition of the classes of the genes differentially regulated at individual time points was unique, suggesting tightly controlled regulation of NO3--responsive gene expression.

Allan-Herndon-Dudley syndrome with unusual profound sensorineural hearing loss.

The Allan-Herndon-Dudley syndrome is caused by mutations in the thyroid hormone transporter, Monocarboxylate transporter 8 (MCT8). It is characterized by profound intellectual disability and abnormal thyroid function. We report on a patient with Allan-Herndon-Dudley syndrome (AHDS) with profound sensorineural hearing loss which is not usually a feature of AHDS and which may have been due to a coexisting nonsense mutation in Microphthalmia-associated transcription factor (MITF).

A comparative analysis of algorithms for somatic SNV detection in cancer.

MOTIVATION: With the advent of relatively affordable high-throughput technologies, DNA sequencing of cancers is now common practice in cancer research projects and will be increasingly used in clinical practice to inform diagnosis and treatment. Somatic (cancer-only) single nucleotide variants (SNVs) are the simplest class of mutation, yet their identification in DNA sequencing data is confounded by germline polymorphisms, tumour heterogeneity and sequencing and analysis errors. Four recently published algorithms for the detection of somatic SNV sites in matched cancer-normal sequencing datasets are VarScan, SomaticSniper, JointSNVMix and Strelka. In this analysis, we apply these four SNV calling algorithms to cancer-normal Illumina exome sequencing of a chronic myeloid leukaemia (CML) patient. The candidate SNV sites returned by each algorithm are filtered to remove likely false positives, then characterized and compared to investigate the strengths and weaknesses of each SNV calling algorithm. RESULTS: Comparing the candidate SNV sets returned by VarScan, SomaticSniper, JointSNVMix2 and Strelka revealed substantial differences with respect to the number and character of sites returned; the somatic probability scores assigned to the same sites; their susceptibility to various sources of noise; and their sensitivities to low-allelic-fraction candidates. AVAILABILITY: Data accession number SRA081939, code at http://code.google.com/p/snv-caller-review/ CONTACT: david.adelson@adelaide.edu.au SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Cell-specific vacuolar calcium storage mediated by CAX1 regulates apoplastic calcium concentration, gas exchange, and plant productivity in Arabidopsis.

The physiological role and mechanism of nutrient storage within vacuoles of specific cell types is poorly understood. Transcript profiles from Arabidopsis thaliana leaf cells differing in calcium concentration ([Ca], epidermis <10 mM versus mesophyll >60 mM) were compared using a microarray screen and single-cell quantitative PCR. Three tonoplast-localized Ca(2+) transporters, CAX1 (Ca(2+)/H(+)-antiporter), ACA4, and ACA11 (Ca(2+)-ATPases), were identified as preferentially expressed in Ca-rich mesophyll. Analysis of respective loss-of-function mutants demonstrated that only a mutant that lacked expression of both CAX1 and CAX3, a gene ectopically expressed in leaves upon knockout of CAX1, had reduced mesophyll [Ca]. Reduced capacity for mesophyll Ca accumulation resulted in reduced cell wall extensibility, stomatal aperture, transpiration, CO(2) assimilation, and leaf growth rate; increased transcript abundance of other Ca(2+) transporter genes; altered expression of cell wall-modifying proteins, including members of the pectinmethylesterase, expansin, cellulose synthase, and polygalacturonase families; and higher pectin concentrations and thicker cell walls. We demonstrate that these phenotypes result from altered apoplastic free [Ca(2+)], which is threefold greater in cax1/cax3 than in wild-type plants. We establish CAX1 as a key regulator of apoplastic [Ca(2+)] through compartmentation into mesophyll vacuoles, a mechanism essential for optimal plant function and productivity.