The human mitochondrial transcriptome.

The human mitochondrial genome comprises a distinct genetic system transcribed as precursor polycistronic transcripts that are subsequently cleaved to generate individual mRNAs, tRNAs, and rRNAs. Here, we provide a comprehensive analysis of the human mitochondrial transcriptome across multiple cell lines and tissues. Using directional deep sequencing and parallel analysis of RNA ends, we demonstrate wide variation in mitochondrial transcript abundance and precisely resolve transcript processing and maturation events. We identify previously undescribed transcripts, including small RNAs, and observe the enrichment of several nuclear RNAs in mitochondria. Using high-throughput in vivo DNaseI footprinting, we establish the global profile of DNA-binding protein occupancy across the mitochondrial genome at single-nucleotide resolution, revealing regulatory features at mitochondrial transcription initiation sites and functional insights into disease-associated variants. This integrated analysis of the mitochondrial transcriptome reveals unexpected complexity in the regulation, expression, and processing of mitochondrial RNA and provides a resource for future studies of mitochondrial function (accessed at http://mitochondria.matticklab.com).

Heterozygous Variants in the Mechanosensitive Ion Channel TMEM63A Result in Transient Hypomyelination during Infancy.

Mechanically activated (MA) ion channels convert physical forces into electrical signals. Despite the importance of this function, the involvement of mechanosensitive ion channels in human disease is poorly understood. Here we report heterozygous missense mutations in the gene encoding the MA ion channel TMEM63A that result in an infantile disorder resembling a hypomyelinating leukodystrophy. Four unrelated individuals presented with congenital nystagmus, motor delay, and deficient myelination on serial scans in infancy, prompting the diagnosis of Pelizaeus-Merzbacher (like) disease. Genomic sequencing revealed that all four individuals carry heterozygous missense variants in the pore-forming domain of TMEM63A. These variants were confirmed to have arisen de novo in three of the four individuals. While the physiological role of TMEM63A is incompletely understood, it is highly expressed in oligodendrocytes and it has recently been shown to be a MA ion channel. Using patch clamp electrophysiology, we demonstrated that each of the modeled variants result in strongly attenuated stretch-activated currents when expressed in naive cells. Unexpectedly, the clinical evolution of all four individuals has been surprisingly favorable, with substantial improvements in neurological signs and developmental progression. In the three individuals with follow-up scans after 4 years of age, the myelin deficit had almost completely resolved. Our results suggest a previously unappreciated role for mechanosensitive ion channels in myelin development.

De Novo Variants in WDR37 Are Associated with Epilepsy, Colobomas, Dysmorphism, Developmental Delay, Intellectual Disability, and Cerebellar Hypoplasia.

WD40 repeat-containing proteins form a large family of proteins present in all eukaryotes. Here, we identified five pediatric probands with de novo variants in WDR37, which encodes a member of the WD40 repeat protein family. Two probands shared one variant and the others have variants in nearby amino acids outside the WD40 repeats. The probands exhibited shared phenotypes of epilepsy, colobomas, facial dysmorphology reminiscent of CHARGE syndrome, developmental delay and intellectual disability, and cerebellar hypoplasia. The WDR37 protein is highly conserved in vertebrate and invertebrate model organisms and is currently not associated with a human disease. We generated a null allele of the single Drosophila ortholog to gain functional insights and replaced the coding region of the fly gene CG12333/wdr37 with GAL4. These flies are homozygous viable but display severe bang sensitivity, a phenotype associated with seizures in flies. Additionally, the mutant flies fall when climbing the walls of the vials, suggesting a defect in grip strength, and repeat the cycle of climbing and falling. Similar to wall clinging defect, mutant males often lose grip of the female abdomen during copulation. These phenotypes are rescued by using the GAL4 in the CG12333/wdr37 locus to drive the UAS-human reference WDR37 cDNA. The two variants found in three human subjects failed to rescue these phenotypes, suggesting that these alleles severely affect the function of this protein. Taken together, our data suggest that variants in WDR37 underlie a novel syndromic neurological disorder.

Biallelic PI4KA variants cause neurological, intestinal and immunological disease.

Phosphatidylinositol 4-kinase IIIα (PI4KIIIα/PI4KA/OMIM:600286) is a lipid kinase generating phosphatidylinositol 4-phosphate (PI4P), a membrane phospholipid with critical roles in the physiology of multiple cell types. PI4KIIIα's role in PI4P generation requires its assembly into a heterotetrameric complex with EFR3, TTC7 and FAM126. Sequence alterations in two of these molecular partners, TTC7 (encoded by TTC7A or TCC7B) and FAM126, have been associated with a heterogeneous group of either neurological (FAM126A) or intestinal and immunological (TTC7A) conditions. Here we show that biallelic PI4KA sequence alterations in humans are associated with neurological disease, in particular hypomyelinating leukodystrophy. In addition, affected individuals may present with inflammatory bowel disease, multiple intestinal atresia and combined immunodeficiency. Our cellular, biochemical and structural modelling studies indicate that PI4KA-associated phenotypical outcomes probably stem from impairment of PI4KIIIα-TTC7-FAM126's organ-specific functions, due to defective catalytic activity or altered intra-complex functional interactions. Together, these data define PI4KA gene alteration as a cause of a variable phenotypical spectrum and provide fundamental new insight into the combinatorial biology of the PI4KIIIα-FAM126-TTC7-EFR3 molecular complex.

Patient-iPSC-Derived Kidney Organoids Show Functional Validation of a Ciliopathic Renal Phenotype and Reveal Underlying Pathogenetic Mechanisms.

Despite the increasing diagnostic rate of genomic sequencing, the genetic basis of more than 50% of heritable kidney disease remains unresolved. Kidney organoids differentiated from induced pluripotent stem cells (iPSCs) of individuals affected by inherited renal disease represent a potential, but unvalidated, platform for the functional validation of novel gene variants and investigation of underlying pathogenetic mechanisms. In this study, trio whole-exome sequencing of a prospectively identified nephronophthisis (NPHP) proband and her parents identified compound-heterozygous variants in IFT140, a gene previously associated with NPHP-related ciliopathies. IFT140 plays a key role in retrograde intraflagellar transport, but the precise downstream cellular mechanisms responsible for disease presentation remain unknown. A one-step reprogramming and gene-editing protocol was used to derive both uncorrected proband iPSCs and isogenic gene-corrected iPSCs, which were differentiated to kidney organoids. Proband organoid tubules demonstrated shortened, club-shaped primary cilia, whereas gene correction rescued this phenotype. Differential expression analysis of epithelial cells isolated from organoids suggested downregulation of genes associated with apicobasal polarity, cell-cell junctions, and dynein motor assembly in proband epithelial cells. Matrigel cyst cultures confirmed a polarization defect in proband versus gene-corrected renal epithelium. As such, this study represents a "proof of concept" for using proband-derived iPSCs to model renal disease and illustrates dysfunctional cellular pathways beyond the primary cilium in the setting of IFT140 mutations, which are established for other NPHP genotypes.

Workplace stress, common mental disorder and suicidal ideation in junior doctors.

BACKGROUND: Doctors-in-training report elevated rates of mental disorders and high levels of stress. Whilst a number of work-related sources of stress have been identified in the medical profession, it remains unclear as to the relative importance of workplace stressors for mental ill-health in junior doctors. AIMS: To examine workplace stressors reported by junior doctors and identify variables associated with adverse mental health outcomes. METHODS: Cross-sectional analysis of national 2013 survey of Australian doctors focussing on junior medical officers (N = 3053; 24.9% of total sample). Primary outcomes were caseness of common mental disorder (CMD) and suicidal ideation in the past year. RESULTS: Perceived level of conflict between study/career and family/personal responsibility (OR = 3.76, 95% CI: 2.61-5.43; P < 0.01) and sleep deprivation (OR = 2.19, 95% CI: 1.46-3.28; P < 0.01) were significantly associated with CMD, while perceived level of conflict between study/career and family/personal responsibility (OR = 3.13, 95% CI: 1.78-5.50; P < 0.01) and bullying (OR = 2.92, 95% CI: 1.42-6.03; P < 0.01) were most strongly associated with suicidal ideation in adjusted models. CONCLUSION: This study identifies modifiable workplace variables that are influential in junior doctors' mental health, and in doing so, provides meaningful evidence-informed targets for future interventions to prevent suicide and mental disorder in this population.

Help-seeking for depression among Australian doctors.

BACKGROUND: Depression is common among doctors. However, concerns remain that doctors are unlikely to ask for help when symptoms of depression arise. AIMS: To determine rates and patterns of help-seeking for depression among doctors and to identify predictors of and barriers to such behaviour. METHODS: A secondary analysis was conducted on a nation-wide survey of 12 252 Australian doctors. The study sample consisted of doctors who reported having ever felt seriously depressed (n = 4154; 33.9% of total sample). Rates of help-seeking, professional help-seeking behaviours and self-reported barriers were explored. Logistic regression was used to examine the association between professional help-seeking and predetermined predictive factors. RESULTS: Sixty percent (95% confidence interval (CI): 58.5-61.5) of doctors who have ever felt seriously depressed reported some form of professional help-seeking for depression. The most common barrier to help-seeking was 'privacy/confidentiality'. Females (odds ratio (OR) = 1.74; 95% CI: 1.50-2.01; P < 0.001), locally trained doctors (OR = 1.34; 95% CI: 1.12-1.59; P = 0.001) and senior doctors (OR = 1.35; 95% CI: 1.14-1.61; P = 0.001) were more likely to seek professional help than their counterparts. Compared with general practitioners, psychiatrists (OR = 1.565; 95% CI: 1.15-2.13; P = 0.004) were more likely to seek professional help while surgeons (OR = 0.518; 95% CI: 0.37-0.72; P < 0.001) and pathologists/radiologists (OR = 0.695; 95% CI: 0.49-0.99; P = 0.043) were less likely. CONCLUSION: While it is reassuring that the majority of depressed doctors were able to seek professional help, many were not. Major barriers to professional help-seeking, particularly concerns about confidentiality and impact on career, remain a problem. Male, overseas-trained, junior doctors, surgeons and pathologists/radiologists were less likely to seek help for depression. Targeted interventions are required to increase appropriate help-seeking for depression in doctors.

'The hardest job I've ever done': a qualitative exploration of the factors affecting junior doctors' mental health and well-being during medical training in Australia.

BACKGROUND: Medical practitioners can experience considerable stress and poor mental health during their careers, with doctors in training known to be particularly vulnerable. Previous research has documented work-related factors that may play a role in the mental health status of junior doctors. However, these and additional factors, need to be explored further by considering theory-driven, social, structural and contextual issues. This qualitative study aimed to explore the experiences of junior doctors working in Australian hospitals to identify factors that impact their mental health during medical training. METHOD: Semi-structured interviews were conducted with 12 junior medical officers (JMOs) employed across six hospitals in Australia. Transcribed de-identified interviews were analysed thematically using a data-driven inductive approach. RESULTS: Four interrelated main themes were identified: i) professional hierarchies; ii) occupational stress; iii) emotional labour, and iv) taking distress home; which detail the complex affective, relational and professional experiences of JMOs. The accounts demonstrate how the social, professional and organisational dimensions of these experiences impact upon trainee's well-being and mental health, both positively and negatively. Together, the findings document the dynamic, nuanced aspects of junior doctors' experiences of medical training and practice and highlights the importance of relational connections and the workplace environment in shaping JMOs' social and emotional well-being. CONCLUSION: The current study adds to the understanding of how junior doctors navigate medical training in Australian hospitals and highlights the complexities of this experience, particularly the ways in which mental health and well-being are shaped by different elements. These findings have important implications to inform new strategies to improve JMO mental health and to leverage work and non-work contexts to better support JMOs during medical training.

Long Noncoding RNAs CUPID1 and CUPID2 Mediate Breast Cancer Risk at 11q13 by Modulating the Response to DNA Damage.

Breast cancer risk is strongly associated with an intergenic region on 11q13. We have previously shown that the strongest risk-associated SNPs fall within a distal enhancer that regulates CCND1. Here, we report that, in addition to regulating CCND1, this enhancer regulates two estrogen-regulated long noncoding RNAs, CUPID1 and CUPID2. We provide evidence that the risk-associated SNPs are associated with reduced chromatin looping between the enhancer and the CUPID1 and CUPID2 bidirectional promoter. We further show that CUPID1 and CUPID2 are predominantly expressed in hormone-receptor-positive breast tumors and play a role in modulating pathway choice for the repair of double-strand breaks. These data reveal a mechanism for the involvement of this region in breast cancer.

Identification of miR-29b targets using 3-cyanovinylcarbazole containing mimics.

MicroRNAs (miRNAs) are highly conserved ∼22 nt small noncoding RNAs that bind partially complementary sequences in target transcripts. MicroRNAs regulate both translation and transcript stability, and play important roles in development, cellular homeostasis, and disease. There are limited approaches available to agnostically identify microRNA targets transcriptome-wide, and methods using miRNA mimics, which in principle identify direct miRNA:transcript pairs, have low sensitivity and specificity. Here, we describe a novel method to identify microRNA targets using miR-29b mimics containing 3-cyanovinylcarbazole (CNVK), a photolabile nucleoside analog. We demonstrate that biotin-tagged, CNVK-containing miR-29b (CNVK-miR-29b) mimics are nontoxic in cell culture, associate with endogenous mammalian Argonaute2, are sensitive for known targets and recapitulate endogenous transcript destabilization. Partnering CNVK-miR-29b with ultra-low-input RNA sequencing, we recover ∼40% of known miR-29b targets and find conservation of the focal adhesion and apoptotic target pathways in mouse and human. We also identify hundreds of novel targets, including NRAS, HOXA10, and KLF11, with a validation rate of 71% for a subset of 73 novel target transcripts interrogated using a high-throughput luciferase assay. Consistent with previous reports, we show that both endogenous miR-29b and CNVK-miR-29b are trafficked to the nucleus, but find no evidence of nuclear-specific miR-29b transcript binding. This may indicate that miR-29b nuclear sequestration is a regulatory mechanism in itself. We suggest that CNVK-containing small RNA mimics may find applicability in other experimental models.

Adult-Diagnosed Nonsyndromic Nephronophthisis in Australian Families Caused by Biallelic NPHP4 Variants.

There is increasing appreciation of nephronophthisis (NPHP) as an autosomal recessive cause of kidney failure and earlier stages of chronic kidney disease among adults. We identified 2 families with presumed adult-diagnosed nonsyndromic NPHP and negative diagnostic genetic testing results from our Renal Genetics Clinic. Both had 2 affected siblings without extrarenal phenotypes. After informed consent, research whole-genome sequencing was undertaken. Biallelic NPHP4 variants were identified in trans and clinically confirmed in all 4 affected individuals, confirming a genetic diagnosis. Participant 1 of the first family (F1P1) had kidney failure diagnosed at 19 years of age. An affected younger sibling (F1P2) reached kidney failure at age 15 years after kidney biopsy suggested NPHP. Pathogenic variants detected in NPHP4 in this family were NM_015102.4:c.3766C>T (p.Gln1256*) and a 31-kb deletion affecting exons 12 to 16. In the second family, F2P3 reached kidney failure at age 27 years having undergone kidney biopsy suggesting NPHP. An affected younger sibling (F2P4) has chronic kidney disease stage 4 at age 39 years. The NPHP4 variants detected were NM_015102.4:c.1998_1999del (p.Tyr667Phefs*23) and c.3646G>T (p.Asp1216Tyr). The latter variant was initially missed in diagnostic sequencing due to inadequate NPHP4 coverage (94.3% exonic coverage). With these reports, we identify NPHP4 as an appreciable genetic cause for adult-diagnosed nonsyndromic NPHP that should be considered by adult nephrologists.

Improved definition of the mouse transcriptome via targeted RNA sequencing.

Targeted RNA sequencing (CaptureSeq) uses oligonucleotide probes to capture RNAs for sequencing, providing enriched read coverage, accurate measurement of gene expression, and quantitative expression data. We applied CaptureSeq to refine transcript annotations in the current murine GRCm38 assembly. More than 23,000 regions corresponding to putative or annotated long noncoding RNAs (lncRNAs) and 154,281 known splicing junction sites were selected for targeted sequencing across five mouse tissues and three brain subregions. The results illustrate that the mouse transcriptome is considerably more complex than previously thought. We assemble more complete transcript isoforms than GENCODE, expand transcript boundaries, and connect interspersed islands of mapped reads. We describe a novel filtering pipeline that identifies previously unannotated but high-quality transcript isoforms. In this set, 911 GENCODE neighboring genes are condensed into 400 expanded gene models. Additionally, 594 GENCODE lncRNAs acquire an open reading frame (ORF) when their structure is extended with CaptureSeq. Finally, we validate our observations using current FANTOM and Mouse ENCODE resources.

Loss-of-function alanyl-tRNA synthetase mutations cause an autosomal-recessive early-onset epileptic encephalopathy with persistent myelination defect.

Mutations in genes encoding aminoacyl-tRNA synthetases are known to cause leukodystrophies and genetic leukoencephalopathies-heritable disorders that result in white matter abnormalities in the central nervous system. Here we report three individuals (two siblings and an unrelated individual) with severe infantile epileptic encephalopathy, clubfoot, absent deep tendon reflexes, extrapyramidal symptoms, and persistently deficient myelination on MRI. Analysis by whole exome sequencing identified mutations in the nuclear-encoded alanyl-tRNA synthetase (AARS) in these two unrelated families: the two affected siblings are compound heterozygous for p.Lys81Thr and p.Arg751Gly AARS, and the single affected child is homozygous for p.Arg751Gly AARS. The two identified mutations were found to result in a significant reduction in function. Mutations in AARS were previously associated with an autosomal-dominant inherited form of axonal neuropathy, Charcot-Marie-Tooth disease type 2N (CMT2N). The autosomal-recessive AARS mutations identified in the individuals described here, however, cause a severe infantile epileptic encephalopathy with a central myelin defect and peripheral neuropathy, demonstrating that defects of alanyl-tRNA charging can result in a wide spectrum of disease manifestations.

Genome-wide discovery of human splicing branchpoints.

During the splicing reaction, the 5' intron end is joined to the branchpoint nucleotide, selecting the next exon to incorporate into the mature RNA and forming an intron lariat, which is excised. Despite a critical role in gene splicing, the locations and features of human splicing branchpoints are largely unknown. We use exoribonuclease digestion and targeted RNA-sequencing to enrich for sequences that traverse the lariat junction and, by split and inverted alignment, reveal the branchpoint. We identify 59,359 high-confidence human branchpoints in >10,000 genes, providing a first map of splicing branchpoints in the human genome. Branchpoints are predominantly adenosine, highly conserved, and closely distributed to the 3' splice site. Analysis of human branchpoints reveals numerous novel features, including distinct features of branchpoints for alternatively spliced exons and a family of conserved sequence motifs overlapping branchpoints we term B-boxes, which exhibit maximal nucleotide diversity while maintaining interactions with the keto-rich U2 snRNA. Different B-box motifs exhibit divergent usage in vertebrate lineages and associate with other splicing elements and distinct intron-exon architectures, suggesting integration within a broader regulatory splicing code. Lastly, although branchpoints are refractory to common mutational processes and genetic variation, mutations occurring at branchpoint nucleotides are enriched for disease associations.

Quantitative gene profiling of long noncoding RNAs with targeted RNA sequencing.

We compared quantitative RT-PCR (qRT-PCR), RNA-seq and capture sequencing (CaptureSeq) in terms of their ability to assemble and quantify long noncoding RNAs and novel coding exons across 20 human tissues. CaptureSeq was superior for the detection and quantification of genes with low expression, showed little technical variation and accurately measured differential expression. This approach expands and refines previous annotations and simultaneously generates an expression atlas.

A recurrent de novo mutation in TMEM106B causes hypomyelinating leukodystrophy.

Hypomyelinating leukodystrophies are a heterogeneous group of disorders with a clinical presentation that often includes early-onset nystagmus, ataxia and spasticity and a wide range of severity. Using next-generation sequencing techniques and GeneMatcher, we identified four unrelated patients with brain hypomyelination, all with the same recurrent dominant mutation, c.754G>A p.(Asp252Asn), in TMEM106B. The mutation was confirmed as de novo in three of the cases, and the mildly affected father of the fourth affected individual was confirmed as mosaic for this variant. The protein encoded by TMEM106B is poorly characterized but is reported to have a role in regulation of lysosomal trafficking. Polymorphisms in TMEM106B are thought to modify disease onset in frontotemporal dementia, but its relation to myelination is not understood. Clinical presentation in three of the four patients is remarkably benign compared to other hypomyelinating disorders, with congenital nystagmus and mild motor delay. These findings add TMEM106B to the growing list of genes causing hypomyelinating disorders and emphasize the essential role lysosomes play in myelination.

X-linked hypomyelination with spondylometaphyseal dysplasia (H-SMD) associated with mutations in AIFM1.

An X-linked condition characterized by the combination of hypomyelinating leukodystrophy and spondylometaphyseal dysplasia (H-SMD) has been observed in only four families, with linkage to Xq25-27, and recent genetic characterization in two families with a common AIFM1 mutation. In our study, 12 patients (6 families) with H-SMD were identified and underwent comprehensive assessment accompanied by whole-exome sequencing (WES). Pedigree analysis in all families was consistent with X-linked recessive inheritance. Presentation typically occurred between 12 and 36 months. In addition to the two disease-defining features of spondylometaphyseal dysplasia and hypomyelination on MRI, common clinical signs and symptoms included motor deterioration, spasticity, tremor, ataxia, dysarthria, cognitive defects, pulmonary hypertension, nystagmus, and vision loss due to retinopathy. The course of the disease was slowly progressive. All patients had maternally inherited or de novo mutations in or near exon 7 of AIFM1, within a region of 70 bp, including synonymous and intronic changes. AIFM1 mutations have previously been associated with neurologic presentations as varied as intellectual disability, hearing loss, neuropathy, and striatal necrosis, while AIFM1 mutations in this small region present with a distinct phenotype implicating bone. Analysis of cell lines derived from four patients identified significant reductions in AIFM1 mRNA and protein levels in osteoblasts. We hypothesize that AIFM1 functions in bone metabolism and myelination and is responsible for the unique phenotype in this condition.

Whole exome sequencing in patients with white matter abnormalities.

Here we report whole exome sequencing (WES) on a cohort of 71 patients with persistently unresolved white matter abnormalities with a suspected diagnosis of leukodystrophy or genetic leukoencephalopathy. WES analyses were performed on trio, or greater, family groups. Diagnostic pathogenic variants were identified in 35% (25 of 71) of patients. Potentially pathogenic variants were identified in clinically relevant genes in a further 7% (5 of 71) of cases, giving a total yield of clinical diagnoses in 42% of individuals. These findings provide evidence that WES can substantially decrease the number of unresolved white matter cases. Ann Neurol 2016;79:1031-1037.

A de novo mutation in the ß-tubulin gene TUBB4A results in the leukoencephalopathy hypomyelination with atrophy of the basal ganglia and cerebellum.

Hypomyelination with atrophy of the basal ganglia and cerebellum (H-ABC) is a rare hereditary leukoencephalopathy that was originally identified by MRI pattern analysis, and it has thus far defied all attempts at identifying the causal mutation. Only 22 cases are published in the literature to date. We performed exome sequencing on five family trios, two family quartets, and three single probands, which revealed that all eleven H-ABC-diagnosed individuals carry the same de novo single-nucleotide TUBB4A mutation resulting in nonsynonymous change p.Asp249Asn. Detailed investigation of one of the family quartets with the singular finding of an H-ABC-affected sibling pair revealed maternal mosaicism for the mutation, suggesting that rare de novo mutations that are initially phenotypically neutral in a mosaic individual can be disease causing in the subsequent generation. Modeling of TUBB4A shows that the mutation creates a nonsynonymous change at a highly conserved asparagine that sits at the intradimer interface of α-tubulin and ß-tubulin, and this change might affect tubulin dimerization, microtubule polymerization, or microtubule stability. Consistent with H-ABC's clinical presentation, TUBB4A is highly expressed in neurons, and a recent report has shown that an N-terminal alteration is associated with a heritable dystonia. Together, these data demonstrate that a single de novo mutation in TUBB4A results in H-ABC.