Somatic mutational landscape of hereditary hematopoietic malignancies caused by germline variants in RUNX1, GATA2, and DDX41

Individuals with germ line variants associated with hereditary hematopoietic malignancies (HHMs) have a highly variable risk for leukemogenesis. Gaps in our understanding of premalignant states in HHMs have hampered efforts to design effective clinical surveillance programs, provide personalized preemptive treatments, and inform appropriate counseling for patients. We used the largest known comparative international cohort of germline RUNX1, GATA2, or DDX41 variant carriers without and with hematopoietic malignancies (HMs) to identify patterns of genetic drivers that are unique to each HHM syndrome before and after leukemogenesis. These patterns included striking heterogeneity in rates of early-onset clonal hematopoiesis (CH), with a high prevalence of CH in RUNX1 and GATA2 variant carriers who did not have malignancies (carriers-without HM). We observed a paucity of CH in DDX41 carriers-without HM. In RUNX1 carriers-without HM with CH, we detected variants in TET2, PHF6, and, most frequently, BCOR. These genes were recurrently mutated in RUNX1-driven malignancies, suggesting CH is a direct precursor to malignancy in RUNX1-driven HHMs. Leukemogenesis in RUNX1 and DDX41 carriers was often driven by second hits in RUNX1 and DDX41, respectively. This study may inform the development of HHM-specific clinical trials and gene-specific approaches to clinical monitoring. For example, trials investigating the potential benefits of monitoring DDX41 carriers-without HM for low-frequency second hits in DDX41 may now be beneficial. Similarly, trials monitoring carriers-without HM with RUNX1 germ line variants for the acquisition of somatic variants in BCOR, PHF6, and TET2 and second hits in RUNX1 are warranted.

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.

Novel modes of MPL activation in triple-negative myeloproliferative neoplasms

The identification of a somatic mutation associated with myeloid malignancy is of diagnostic importance in myeloproliferative neoplasms (MPNs). Individuals with no mutation detected in common screening tests for variants in JAK2, CALR, and MPL are described as 'triple-negative' and pose a diagnostic challenge if there is no other evidence of a clonal disorder. To identify potential drivers that might explain the clinical phenotype, we used an extended sequencing panel to characterise a cohort of 44 previously diagnosed triple-negative MPN patients for canonical mutations in JAK2, MPL and CALR at low variant allele frequency (found in 4/44 patients), less common variants in the JAK-STAT signalling pathway (12 patients), or other variants in recurrently mutated genes from myeloid malignancies (18 patients), including hotspot variants of potential clinical relevance in eight patients. In one patient with thrombocytosis we identified biallelic germline MPL variants. Neither MPL variant was activating in cell proliferation assays, and one of the variants was not expressed on the cell surface, yet co-expression of both variants led to thrombopoietin hypersensitivity. Our results highlight the clinical value of extended sequencing including germline variant analysis and illustrate the need for detailed functional assays to determine whether rare variants in JAK2 or MPL are pathogenic.

Shariant platform: Enabling evidence sharing across Australian clinical genetic-testing laboratories to support variant interpretation

Sharing genomic variant interpretations across laboratories promotes consistency in variant assertions. A landscape analysis of Australian clinical genetic-testing laboratories in 2017 identified that, despite the national-accreditation-body recommendations encouraging laboratories to submit genotypic data to clinical databases, fewer than 300 variants had been shared to the ClinVar public database. Consultations with Australian laboratories identified resource constraints limiting routine application of manual processes, consent issues, and differences in interpretation systems as barriers to sharing. This information was used to define key needs and solutions required to enable national sharing of variant interpretations. The Shariant platform, using both the GRCh37 and GRCh38 genome builds, was developed to enable ongoing sharing of variant interpretations and associated evidence between Australian clinical genetic-testing laboratories. Where possible, two-way automated sharing was implemented so that disruption to laboratory workflows would be minimized. Terms of use were developed through consultation and currently restrict access to Australian clinical genetic-testing laboratories. Shariant was designed to store and compare structured evidence, to promote and record resolution of inter-laboratory classification discrepancies, and to streamline the submission of variant assertions to ClinVar. As of December 2021, more than 14,000 largely prospectively curated variant records from 11 participating laboratories have been shared. Discrepant classifications have been identified for 11% (28/260) of variants submitted by more than one laboratory. We have demonstrated that co-design with clinical laboratories is vital to developing and implementing a national variant-interpretation sharing effort. This approach has improved inter-laboratory concordance and enabled opportunities to standardize interpretation practices.

RNA-Based Targeted Gene Sequencing Improves the Diagnostic Yield of Mutant Detection in Chronic Myeloid Leukemia

Mutation detection is increasingly used for the management of hematological malignancies. Prior whole transcriptome and whole exome sequencing studies using total RNA and DNA identified diverse mutation types in cancer-related genes associated with treatment failure in patients with chronic myeloid leukemia. Variants included single-nucleotide variants and small insertions/deletions, plus fusion transcripts and partial or whole gene deletions. The hypothesis that all of these mutation types could be detected by a single cost-effective hybridization capture next-generation sequencing method using total RNA was assessed. A method was developed that targeted 130 genes relevant for myeloid and lymphoid leukemia. Retrospective samples with 121 precharacterized variants were tested using total RNA and/or DNA. Concordance of detection of precharacterized variants using RNA or DNA was 96%, whereas the enhanced sensitivity identified additional variants. Comparison between 24 matched DNA and RNA samples demonstrated 95.3% of 170 variants detectable using DNA were detected using RNA, including all but one variant predicted to activate nonsense-mediated decay. RNA identified an additional 10 variants, including fusion transcripts. Furthermore, the true effect of splice variants on RNA splicing was only evident using RNA. In conclusion, capture sequencing using total RNA alone is suitable for detecting a range of variants relevant in chronic myeloid leukemia and may be more broadly applied to other hematological malignancies where diverse variant types define risk groups.

Compound heterozygous variants in LAMC3 in association with posterior periventricular nodular heterotopia

Background

Periventricular nodular heterotopia (PNH) is a malformation of cortical development characterized by nodules of abnormally migrated neurons. The cause of posteriorly placed PNH is not well characterised and we present a case that provides insights into the cause of posterior PNH.

Case presentation

We report a fetus with extensive posterior PNH in association with biallelic variants in LAMC3. LAMC3 mutations have previously been shown to cause polymicrogyria and pachygyria in the occipital cortex, but not PNH. The occipital location of PNH in our case and the proposed function of LAMC3 in cortical development suggest that the identified LAMC3 variants may be causal of PNH in this fetus.

Conclusion

We hypothesise that this finding extends the cortical phenotype associated with LAMC3 and provides valuable insight into genetic cause of posterior PNH.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12920-021-00911-4.

Aberrant RAG-mediated recombination contributes to multiple structural rearrangements in lymphoid blast crisis of chronic myeloid leukemia

Blast crisis of chronic myeloid leukemia is associated with poor survival and the accumulation of genomic lesions. Using whole-exome and/or RNA sequencing of patients at chronic phase (CP, n = 49), myeloid blast crisis (MBC, n = 19), and lymphoid blast crisis (LBC, n = 20), we found 25 focal gene deletions and 14 fusions in 24 patients in BC. Deletions predominated in LBC (83% of structural variants). Transcriptional analysis identified the upregulation of genes involved in V(D)J recombination, including RAG1/2 and DNTT in LBC. RAG recombination is a reported mediator of IKZF1 deletion. We investigated the extent of RAG-mediated genomic lesions in BC. Molecular hallmarks of RAG activity; DNTT-mediated nucleotide insertions and a RAG-binding motif at structural variants were exclusively found in patients with high RAG expression. Structural variants in 65% of patients in LBC displayed these hallmarks compared with only 5% in MBC. RAG-mediated events included focal deletion and novel fusion of genes associated with hematologic cancer: IKZF1, RUNX1, CDKN2A/B, and RB1. Importantly, 8/8 patients with elevated DNTT at CP diagnosis progressed to LBC by 12 months, potentially enabling early prediction of LBC. This work confirms the central mutagenic role of RAG in LBC and describes potential clinical utility in CML management.

Aberrant Splicing of SDHC in Families With Unexplained Succinate Dehydrogenase-Deficient Paragangliomas

Context

Germline mutations in the succinate dehydrogenase genes (SDHA/B/C/D, SDHAF2-collectively, "SDHx") have been implicated in paraganglioma (PGL), renal cell carcinoma (RCC), gastrointestinal stromal tumor (GIST), and pituitary adenoma (PA). Negative SDHB tumor staining is indicative of SDH-deficient tumors, usually reflecting an underlying germline SDHx mutation. However, approximately 20% of individuals with SDH-deficient tumors lack an identifiable germline SDHx mutation.

Methods

We performed whole-exome sequencing (WES) of germline and tumor DNA followed by Sanger sequencing validation, transcriptome analysis, metabolomic studies, and haplotype analysis in 2 Italian-Australian families with SDH-deficient PGLs and various neoplasms, including RCC, GIST, and PA.

Results

Germline WES revealed a novel SDHC intronic variant, which had been missed during previous routine testing, in 4 affected siblings of the index family. Transcriptome analysis demonstrated aberrant SDHC splicing, with the retained intronic segment introducing a premature stop codon. WES of available tumors in this family showed chromosome 1 deletion with loss of wild-type SDHC in a PGL and a somatic gain-of-function KIT mutation in a GIST. The SDHC intronic variant identified was subsequently detected in the second family, with haplotype analysis indicating a founder effect.

Conclusions

This is the deepest intronic variant to be reported among the SDHx genes. Intronic variants beyond the limits of standard gene sequencing analysis should be considered in patients with SDH-deficient tumors but negative genetic test results.

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.

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.

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.

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.

The Genomic Landscape of Sporadic Prolactinomas

Somatic GNAS and USP8 mutations have been implicated in sporadic somatotrophinomas and corticotrophinomas, respectively. However, no genes are known to be recurrently mutated in sporadic prolactinomas. The prevalence of copy number variants (CNV), which is emerging as a mechanism of tumorigenesis in sporadic pituitary adenomas in general, is also unclear in prolactinomas. To characterize the genetic events underpinning sporadic prolactinomas, we performed whole exome sequencing of paired tumor and germline DNA from 12 prolactinoma patients. We observed recurrent large-scale CNV, most commonly in the form of copy number gains. We also identified sequence variants of interest in 15 genes. This included the DRD2, PRL, TMEM67, and MLH3 genes with plausible links to prolactinoma formation. Of the 15 genes of interest, CNV was seen at the gene locus in the corresponding tumor in 10 cases, and pituitary expression of eight genes was in the top 10% of tissues. However, none of our shortlisted somatic variants appeared to be classical driver mutations as no variant was found in more than one tumor. Future directions of research include mechanistic studies to investigate how CNV may contribute to prolactinoma formation, larger studies of relevant prolactinoma subsets according to clinical characteristics, and additional genetic investigations for aberrations not captured by whole exome sequencing.

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.

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.

BAM-matcher: a tool for rapid NGS sample matching

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.

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.

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.

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.

ARMC5 mutations are common in familial bilateral macronodular adrenal hyperplasia

CONTEXT

Bilateral macronodular adrenal hyperplasia (BMAH) is a rare form of adrenal Cushing's syndrome. Familial cases have been reported, but at the time we conducted this study, the genetic basis of BMAH was unknown. Recently, germline variants of armadillo repeat containing 5 (ARMC5) in patients with isolated BMAH and somatic, second-hit mutations in tumor nodules, were identified.

OBJECTIVE

Our objective was to identify the genetic basis of familial BMAH.

DESIGN

We performed whole exome capture and sequencing of 2 affected individuals from each of 4 BMAH families (BMAH-01, BMAH-02, BMAH-03, and BMAH-05). Based on clinical evaluation, there were 7, 3, 3, and 4 affected individuals in these families, respectively. Sanger sequencing of ARMC5 was performed in 1 other BMAH kindred, BMAH-06.

RESULTS

Exome sequencing identified novel variants Chr16:g.31477540, c.2139delT, p.(Thr715Leufs*1) (BMAH-02) and Chr16:g.31473811, c.943C→T, p.(Arg315Trp) (BMAH-03) in ARMC5 (GRch37/hg19), validated by Sanger sequencing. BMAH-01 had a recently reported mutation Chr16:g.31476121, c.1777C→T, p.(Arg593Trp). Sanger sequencing of ARMC5 in BMAH-06 identified a previously reported mutation, Chr16:g. 31473688; c.799C→T, p.(Arg267*). The genetic basis of BMAH in BMAH-05 was not identified.

CONCLUSIONS

Our studies have detected ARMC5 mutations in 4 of 5 BMAH families tested, confirming that these mutations are a frequent cause of BMAH. Two of the 4 families had novel mutations, indicating allelic heterogeneity. Preclinical evaluation did not predict mutation status. The ARMC5-negative family had unusual prominent hyperaldosteronism. Further studies are needed to determine the penetrance of BMAH in ARMC5 mutation-positive relatives of affected patients, the practical utility of genetic screening and genotype-phenotype correlations.

Sequencing error correction without a reference genome

BACKGROUND

Next (second) generation sequencing is an increasingly important tool for many areas of molecular biology, however, care must be taken when interpreting its output. Even a low error rate can cause a large number of errors due to the high number of nucleotides being sequenced. Identifying sequencing errors from true biological variants is a challenging task. For organisms without a reference genome this difficulty is even more challenging.

RESULTS

We have developed a method for the correction of sequencing errors in data from the Illumina Solexa sequencing platforms. It does not require a reference genome and is of relevance for microRNA studies, unsequenced genomes, variant detection in ultra-deep sequencing and even for RNA-Seq studies of organisms with sequenced genomes where RNA editing is being considered.

CONCLUSIONS

The derived error model is novel in that it allows different error probabilities for each position along the read, in conjunction with different error rates depending on the particular nucleotides involved in the substitution, and does not force these effects to behave in a multiplicative manner. The model provides error rates which capture the complex effects and interactions of the three main known causes of sequencing error associated with the Illumina platforms.

Nutlin-3a efficacy in sarcoma predicted by transcriptomic and epigenetic profiling

Nutlin-3a is a small-molecule antagonist of p53/MDM2 that is being explored as a treatment for sarcoma. In this study, we examined the molecular mechanisms underlying the sensitivity of sarcomas to Nutlin-3a. In an ex vivo tissue explant system, we found that TP53 pathway alterations (TP53 status, MDM2/MDM4 genomic amplification/mRNA overexpression, MDM2 SNP309, and TP53 SNP72) did not confer apoptotic or cytostatic responses in sarcoma tissue biopsies (n = 24). Unexpectedly, MDM2 status did not predict Nutlin-3a sensitivity. RNA sequencing revealed that the global transcriptomic profiles of these sarcomas provided a more robust prediction of apoptotic responses to Nutlin-3a. Expression profiling revealed a subset of TP53 target genes that were transactivated specifically in sarcomas that were highly sensitive to Nutlin-3a. Of these target genes, the GADD45A promoter region was shown to be hypermethylated in 82% of wild-type TP53 sarcomas that did not respond to Nutlin-3a, thereby providing mechanistic insight into the innate ability of sarcomas to resist apoptotic death following Nutlin-3a treatment. Collectively, our findings argue that the existing benchmark biomarker for MDM2 antagonist efficacy (MDM2 amplification) should not be used to predict outcome but rather global gene expression profiles and epigenetic status of sarcomas dictate their sensitivity to p53/MDM2 antagonists.

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.

Clusters of genes encoding fructan biosynthesizing enzymes in wheat and barley

Fructans are soluble carbohydrates with health benefits and possible roles in plant adaptation. Fructan biosynthetic genes were isolated using comparative genomics and physical mapping followed by BAC sequencing in barley. Genes encoding sucrose:sucrose 1-fructosyltransferase (1-SST), fructan:fructan 1-fructosyltransferase (1-FFT) and sucrose:fructan 6-fructosyltransferase (6-SFT) were clustered together with multiple copies of vacuolar invertase genes and a transposable element on two barley BAC. Intron-exon structures of the genes were similar. Phylogenetic analysis of the fructosyltransferases and invertases in the Poaceae showed that the fructan biosynthetic genes may have evolved from vacuolar invertases. Quantitative real-time PCR was performed using leaf RNA extracted from three wheat cultivars grown under different conditions. The 1-SST, 1-FFT and 6-SFT genes had correlated expression patterns in our wheat experiment and in existing barley transcriptome database. Single nucleotide polymorphism (SNP) markers were developed and successfully mapped to a major QTL region affecting wheat grain fructan accumulation in two independent wheat populations. The alleles controlling high- and low- fructan in parental lines were also found to be associated in fructan production in a diverse set of 128 wheat lines. To the authors' knowledge, this is the first report on the mapping and sequencing of a fructan biosynthetic gene cluster and in particular, the isolation of a novel 1-FFT gene from barley.

Transcriptome-scale homoeolog-specific transcript assemblies of bread wheat

Background

Bread wheat is one of the world’s most important food crops and considerable efforts have been made to develop genomic resources for this species. This includes an on-going project by the International Wheat Genome Sequencing Consortium to assemble its large and complex genome, which is hexaploid and contains three closely related ‘homoeologous’ copies for each chromosome. This multi-national effort avoids the complications polyploidy entails for correct assembly of the genome by sequencing flow-sorted chromosome arms one at a time. Here we report on an alternate approach, a direct homoeolog-specific assembly of the expressed portion of the genome, the transcriptome.

Results

After assessment of the ability of various assemblers to generate homoeolog-specific assemblies, we employed a two-stage assembly process to produce a high-quality assembly of the transcriptome of hexaploid wheat from Roche-454 and Illumina GAII_x paired-end sequence reads. The assembly process made use of a rapid partitioning of expressed sequences into homoeologous clusters, followed by a parallel high-fidelity assembly of each cluster on a 1150-processor compute cloud. We assessed assembly quality through comparison to known wheat gene sequences and found that in ca. 98.5% of cases the assembly was sufficiently accurate for homoeologous triplets to be cleanly separated into either two or three separate contigs. Comparison to publicly available transcript collections suggests that the assembly covers ~75-80% of the complete transcriptome.

Conclusions

This work therefore describes the first homoeolog-specific sequence assembly of the wheat transcriptome and provides a reference transcriptome for future wheat research. Furthermore, our assembly methodology is transferable to other polyploid organisms.

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.

Comparative transcriptomics in the Triticeae

Background

Barley and particularly wheat are two grass species of immense agricultural importance. In spite of polyploidization events within the latter, studies have shown that genotypically and phenotypically these species are very closely related and, indeed, fertile hybrids can be created by interbreeding. The advent of two genome-scale Affymetrix GeneChips now allows studies of the comparison of their transcriptomes.

Results

We have used the Wheat GeneChip to create a "gene expression atlas" for the wheat transcriptome (cv. Chinese Spring). For this, we chose mRNA from a range of tissues and developmental stages closely mirroring a comparable study carried out for barley (cv. Morex) using the Barley1 GeneChip. This, together with large-scale clustering of the probesets from the two GeneChips into "homologous groups", has allowed us to perform a genomic-scale comparative study of expression patterns in these two species. We explore the influence of the polyploidy of wheat on the results obtained with the Wheat GeneChip and quantify the correlation between conservation in gene sequence and gene expression in wheat and barley. In addition, we show how the conservation of expression patterns can be used to elucidate, probeset by probeset, the reliability of the Wheat GeneChip.

Conclusion

While there are many differences in expression on the level of individual genes and tissues, we demonstrate that the wheat and barley transcriptomes appear highly correlated. This finding is significant not only because given small evolutionary distance between the two species it is widely expected, but also because it demonstrates that it is possible to use the two GeneChips for comparative studies. This is the case even though their probeset composition reflects rather different design principles as well as, of course, the present incomplete knowledge of the gene content of the two species. We also show that, in general, the Wheat GeneChip is not able to distinguish contributions from individual homoeologs. Furthermore, the comparison between the two species leads us to conclude that the conservation of both gene sequence as well as gene expression is positively correlated with absolute expression levels, presumably reflecting increased selection pressure on genes coding for proteins present at high levels. In addition, the results indicate the presence of a correlation between sequence and expression conservation within the Triticeae.

A framework for gene expression analysis

MOTIVATION

Global gene expression measurements as obtained, for example, in microarray experiments can provide important clues to the underlying transcriptional control mechanisms and network structure of a biological cell. In the absence of a detailed understanding of this gene regulation, current attempts at classification of expression data rely on clustering and pattern recognition techniques employing ad-hoc similarity criteria. To improve this situation, a better understanding of the expected relationships between expression profiles of genes associated by biological function is required.

RESULTS

It is shown that perturbation expansions familiar from biological systems theory make precise predictions for the types of relationships to be expected for expression profiles of biologically associated genes, even if the underlying biological factors responsible for this association are not known. Classification criteria are derived, most of which are not usually employed in clustering algorithms. The approach is illustrated by using the AtGenExpress Arabidopsis thaliana developmental expression map.