Detecting genomic indel variants with exact breakpoints in single- and paired-end sequencing data using SplazerS

Motivation: The reliable detection of genomic variation in resequencing data is still a major challenge, especially for variants larger than a few base pairs. Sequencing reads crossing boundaries of structural variation carry the potential for their identification, but are difficult to map. Results: Here we present a method for ‘split’ read mapping, where prefix and suffix match of a read may be interrupted by a longer gap in the read-to-reference alignment. We use this method to accurately detect medium-sized insertions and long deletions with precise breakpoints in genomic resequencing data. Compared with alternative split mapping methods, SplazerS significantly improves sensitivity for detecting large indel events, especially in variant-rich regions. Our method is robust in the presence of sequencing errors as well as alignment errors due to genomic mutations/divergence, and can be used on reads of variable lengths. Our analysis shows that SplazerS is a versatile tool applicable to unanchored or single-end as well as anchored paired-end reads. In addition, application of SplazerS to targeted resequencing data led to the interesting discovery of a complete, possibly functional gene retrocopy variant. Availability: SplazerS is available from http://www.seqan.de/projects/ splazers

Whole genome sequencing identifies a duplicated region encompassing Xq13.2q13.3 in a large Iranian family with intellectual disability

Background The X chromosome has historically been one of the most thoroughly investigated chromosomes regarding intellectual disability (ID), whose etiology is attributed to many factors including copy number variations (CNVs). Duplications of the long arm of the X chromosome have been reported in patients with ID, short stature, facial anomalies, and in many cases hypoplastic genitalia and/or behavioral abnormalities. Methods Here, we report on a large Iranian family with X‐linked ID caused by a duplication on the X chromosome identified by whole genome sequencing in combination with linkage analysis. Results Seven affected males in different branches of the family presented with ID, short stature, seizures, facial anomalies, behavioral abnormalities (aggressiveness, self‐injury, anxiety, impaired social interactions, and shyness), speech impairment, and micropenis. The duplication of the region Xq13.2q13.3, which is ~1.8 Mb in size, includes seven protein‐coding OMIM genes. Three of these genes, namely SLC16A2, RLIM, and NEXMIF, if impaired, can lead to syndromes presenting with ID. Of note, this duplicated region was located within a linkage interval with a LOD score >3. Conclusion Our report indicates that CNVs should be considered in multi‐affected families where no candidate gene defect has been identified in sequencing data analysis

Mutations in the intellectual disability gene Ube2a cause neuronal dysfunction and impair parkin-dependent mitophagy

The prevalence of intellectual disability is around 3%; however, the etiology of the disease remains unclear in most cases. We identified a series of patients with X-linked intellectual disability presenting mutations in the Rad6a (Ube2a) gene, which encodes for an E2 ubiquitin-conjugating enzyme. Drosophila deficient for dRad6 display defective synaptic function as a consequence of mitochondrial failure. Similarly, mouse mRad6a (Ube2a) knockout and patient-derived hRad6a (Ube2a) mutant cells show defective mitochondria. Using in vitro and in vivo ubiquitination assays, we show that RAD6A acts as an E2 ubiquitin-conjugating enzyme that, in combination with an E3 ubiquitin ligase such as Parkin, ubiquitinates mitochondrial proteins to facilitate the clearance of dysfunctional mitochondria in cells. Hence, we identify RAD6A as a regulator of Parkin-dependent mitophagy and establish a critical role for RAD6A in maintaining neuronal function

A noncoding, regulatory mutation implicates HCFC1 in nonsyndromic intellectual disability

The discovery of mutations causing human disease has so far been biased toward protein-coding regions. Having excluded all annotated coding regions, we performed targeted massively parallel resequencing of the nonrepetitive genomic linkage interval at Xq28 of family MRX3. We identified in the binding site of transcription factor YY1 a regulatory mutation that leads to overexpression of the chromatin-associated transcriptional regulator HCFC1. When tested on embryonic murine neural stem cells and embryonic hippocampal neurons, HCFC1 overexpression led to a significant increase of the production of astrocytes and a considerable reduction in neurite growth. Two other nonsynonymous, potentially deleterious changes have been identified by X-exome sequencing in individuals with intellectual disability, implicating HCFC1 in normal brain function

Heterologous expression of a thermophilic diacylglycerol acyltransferase triggers triglyceride accumulation in Escherichia coli

Triglycerides (TAGs), the major storage molecules of metabolic energy and source of fatty acids, are produced as single cell oil by some oleogenic microorganisms. However, these microorganisms require strict culture conditions, show low carbon source flexibilities, lack efficient genetic modification tools and in some cases pose safety concerns. TAGs have essential applications such as behaving as a source for added-value fatty acids or giving rise to the production of biodiesel. Hence, new alternative methods are urgently required for obtaining these oils. In this work we describe TAG accumulation in the industrially appropriate microorganism Escherichia coli expressing the heterologous enzyme tDGAT, a wax ester synthase/triacylglycerol:acylCoA acyltranferase (WS/DGAT). With this purpose, we introduce a codon-optimized gene from the thermophilic actinomycete Thermomonospora curvata coding for a WS/DGAT into different E. coli strains, describe the metabolic effects associated to the expression of this protein and evaluate neutral lipid accumulation. We observe a direct relation between the expression of this WS/DGAT and TAG production within a wide range of culture conditions. More than 30% TAGs were detected within the bacterial neutral lipids in 90 minutes after induction. TAGs were observed to be associated with the hydrophobic enzyme while forming round intracytoplasmic bodies, which could represent a bottleneck for lipid accumulation in E. coli. We detected an increase of almost 3- fold in the monounsaturated fatty acids (MUFA) occurring in the recombinant strains. These MUFA were predominant in the accumulated TAGs achieving 46% of the TAG fatty acids. These results set the basis for further research on the achievement of a suitable method towards the sustainable production of these neutral lipids

Rectal Transmission of Transmitted/Founder HIV-1 Is Efficiently Prevented by Topical 1% Tenofovir in BLT Humanized Mice

Rectal microbicides are being developed to prevent new HIV infections in both men and women. We focused our in vivo preclinical efficacy study on rectally-applied tenofovir. BLT humanized mice (n = 43) were rectally inoculated with either the primary isolate HIV-1(JRCSF) or the MSM-derived transmitted/founder (T/F) virus HIV-1(THRO) within 30 minutes following treatment with topical 1% tenofovir or vehicle. Under our experimental conditions, in the absence of drug treatment we observed 50% and 60% rectal transmission by HIV-1(JRCSF) and HIV-1(THRO), respectively. Topical tenofovir reduced rectal transmission to 8% (1/12; log rank p = 0.03) for HIV-1(JRCSF) and 0% (0/6; log rank p = 0.02) for HIV-1(THRO). This is the first demonstration that any human T/F HIV-1 rectally infects humanized mice and that transmission of the T/F virus can be efficiently blocked by rectally applied 1% tenofovir. These results obtained in BLT mice, along with recent ex vivo, Phase 1 trial and non-human primate reports, provide a critically important step forward in the development of tenofovir-based rectal microbicides

Integration of Hi-C with short and long-read genome sequencing reveals the structure of germline rearranged genomes

Here the authors characterize structural variations (SVs) in a cohort of individuals with complex genomic rearrangements, identifying breakpoints by employing short- and long-read genome sequencing and investigate their impact on gene expression and the three-dimensional chromatin architecture. They find breakpoints are enriched in inactive regions and can result in chromatin domain fusions.Structural variants are a common cause of disease and contribute to a large extent to inter-individual variability, but their detection and interpretation remain a challenge. Here, we investigate 11 individuals with complex genomic rearrangements including germline chromothripsis by combining short- and long-read genome sequencing (GS) with Hi-C. Large-scale genomic rearrangements are identified in Hi-C interaction maps, allowing for an independent assessment of breakpoint calls derived from the GS methods, resulting in >300 genomic junctions. Based on a comprehensive breakpoint detection and Hi-C, we achieve a reconstruction of whole rearranged chromosomes. Integrating information on the three-dimensional organization of chromatin, we observe that breakpoints occur more frequently than expected in lamina-associated domains (LADs) and that a majority reshuffle topologically associating domains (TADs). By applying phased RNA-seq, we observe an enrichment of genes showing allelic imbalanced expression (AIG) within 100 kb around the breakpoints. Interestingly, the AIGs hit by a breakpoint (19/22) display both up- and downregulation, thereby suggesting different mechanisms at play, such as gene disruption and rearrangements of regulatory information. However, the majority of interpretable genes located 200 kb around a breakpoint do not show significant expression changes. Thus, there is an overall robustness in the genome towards large-scale chromosome rearrangements

Rhodococcus opacus B4: a promising bacterium for production of biofuels and biobased chemicals

Bacterial lipids have relevant applications in the production of renewable fuels and biobased oleochemicals. The genus Rhodococcus is one of the most relevant lipid producers due to its capability to accumulate those compounds, mainly triacylglycerols (TAG), when cultivated on different defined substrates, namely sugars, organic acids and hydrocarbons but also on complex carbon sources present in industrial wastes. In this work, the production of storage lipids by Rhodococcus opacus B4 using glucose, acetate and hexadecane is reported for the first time and its productivity compared with Rhodococcus opacus PD630, the best TAG producer bacterium reported. Both strains accumulated mainly TAG from all carbon sources, being influenced by the carbon source itself and by the duration of the accumulation period. R. opacus B4 produced 0.09 and 0.14 g L1 at 24 and 72 h, with hexadecane as carbon source, which was 2 and 3.3 fold higher than the volumetric production obtained by R. opacus PD630. Both strains presented similar fatty acids (FA) profiles in intact cells while in TAG produced fraction, R. opacus B4 revealed a higher variability in fatty acid composition than R. opacus PD630, when both strains were cultivated on hexadecane. The obtained results open new perspectives for the use of R. opacus B4 to produce TAG, in particular using oily (alkane-contaminated) waste and wastewater as cheap raw-materials. Combining TAG production with hydrocarbons degradation is a promising strategy to achieve environmental remediation while producing added value compounds.This work was financially supported by the Portuguese Science Foundation (FCT) and European Social Fund (ESF, POPH-QREN) through the Grant given to A.R. Castro (SFRH/BD/64500/2009), the FCT Strategic Project of UID/BIO/04469/2013 unit and COMPETE 2020 (POCI-01-0145-FEDER-006684) and project RECI/BBB-EBI/0179/2012 (FCOMP-01-0124-FEDER-027462)

Nile red fluorescence screening facilitating neutral lipid phenotype determination in budding yeast, Saccharomyces cerevisiae, and the fission yeast Schizosaccharomyces pombe.

Investigation of yeast neutral lipid accumulation is important for biotechnology and also for modelling aberrant lipid metabolism in human disease. The Nile red (NR) method has been extensively utilised to determine lipid phenotypes of yeast cells via microscopic means. NR assays have been used to differentiate lipid accumulation and relative amounts of lipid in oleaginous species but have not been thoroughly validated for phenotype determination arising from genetic modification. A modified NR assay, first described by Sitepu et al. (J Microbiol Methods 91:321-328, 2012), was able to detect neutral lipid changes in Saccharomyces cerevisiae deletion mutants with sensitivity similar to more advanced methodology. We have also be able to, for the first time, successfully apply the NR assay to the well characterised fission yeast Schizosaccharomyces pombe, an increasingly important organism in biotechnology. The described NR fluorescence assay is suitable for increased throughput and rapid screening of genetically modified strains in both the biotechnology industry and for modelling ectopic lipid production for a variety of human diseases. This ultimately negates the need for labour intensive and time consuming lipid analyses of samples that may not yield a desirable lipid phenotype, whilst genetic modifications impacting significantly on the cellular lipid phenotype can be further promoted for more in depth analyses

Functional and clinical studies reveal pathophysiological complexity of CLCN4-related neurodevelopmental condition

Missense and truncating variants in the X-chromosome-linked CLCN4 gene, resulting in reduced or complete loss-of-function (LOF) of the encoded chloride/proton exchanger ClC-4, were recently demonstrated to cause a neurocognitive phenotype in both males and females. Through international clinical matchmaking and interrogation of public variant databases we assembled a database of 90 rare CLCN4 missense variants in 90 families: 41 unique and 18 recurrent variants in 49 families. For 43 families, including 22 males and 33 females, we collated detailed clinical and segregation data. To confirm causality of variants and to obtain insight into disease mechanisms, we investigated the effect on electrophysiological properties of 59 of the variants in Xenopus oocytes using extended voltage and pH ranges. Detailed analyses revealed new pathophysiological mechanisms: 25% (15/59) of variants demonstrated LOF, characterized by a “shift” of the voltage-dependent activation to more positive voltages, and nine variants resulted in a toxic gain-of-function, associated with a disrupted gate allowing inward transport at negative voltages. Functional results were not always in line with in silico pathogenicity scores, highlighting the complexity of pathogenicity assessment for accurate genetic counselling. The complex neurocognitive and psychiatric manifestations of this condition, and hitherto under-recognized impacts on growth, gastrointestinal function, and motor control are discussed. Including published cases, we summarize features in 122 individuals from 67 families with CLCN4-related neurodevelopmental condition and suggest future research directions with the aim of improving the integrated care for individuals with this diagnosis