Plasticity of repetitive sequences demonstrated by the complete mitochondrial genome of Eucalyptus camaldulensis.

The tree Eucalyptus camaldulensis is a ubiquitous member of the Eucalyptus genus, which includes several hundred species. Despite the extensive sequencing and assembly of nuclear genomes from various eucalypts, the genus has only one fully annotated and complete mitochondrial genome (mitogenome). Plant mitochondria are characterized by dynamic genomic rearrangements, facilitated by repeat content, a feature that has hindered the assembly of plant mitogenomes. This complexity is evident in the paucity of available mitogenomes. This study, to the best of our knowledge, presents the first E. camaldulensis mitogenome. Our findings suggest the presence of multiple isomeric forms of the E. camaldulensis mitogenome and provide novel insights into minor rearrangements triggered by nested repeat sequences. A comparative sequence analysis of the E. camaldulensis and E. grandis mitogenomes unveils evolutionary changes between the two genomes. A significant divergence is the evolution of a large repeat sequence, which may have contributed to the differences observed between the two genomes. The largest repeat sequences in the E. camaldulensis mitogenome align well with significant yet unexplained structural variations in the E. grandis mitogenome, highlighting the adaptability of repeat sequences in plant mitogenomes.

Beyond the exome: utility of long-read whole genome sequencing in exome-negative autosomal recessive diseases.

BACKGROUND: Long-read whole genome sequencing (lrWGS) has the potential to address the technical limitations of exome sequencing in ways not possible by short-read WGS. However, its utility in autosomal recessive Mendelian diseases is largely unknown. METHODS: In a cohort of 34 families in which the suspected autosomal recessive diseases remained undiagnosed by exome sequencing, lrWGS was performed on the Pacific Bioscience Sequel IIe platform. RESULTS: Likely causal variants were identified in 13 (38%) of the cohort. These include (1) a homozygous splicing SV in TYMS as a novel candidate gene for lethal neonatal lactic acidosis, (2) a homozygous non-coding SV that we propose impacts STK25 expression and causes a novel neurodevelopmental disorder, (3) a compound heterozygous SV in RP1L1 with complex inheritance pattern in a family with inherited retinal disease, (4) homozygous deep intronic variants in LEMD2 and SNAP91 as novel candidate genes for neurodevelopmental disorders in two families, and (5) a promoter SNV in SLC4A4 causing non-syndromic band keratopathy. Surprisingly, we also encountered causal variants that could have been identified by short-read exome sequencing in 7 families. The latter highlight scenarios that are especially challenging at the interpretation level. CONCLUSIONS: Our data highlight the continued need to address the interpretation challenges in parallel with efforts to improve the sequencing technology itself. We propose a path forward for the implementation of lrWGS sequencing in the setting of autosomal recessive diseases in a way that maximizes its utility.

LeafGo: Leaf to Genome, a quick workflow to produce high-quality de novo plant genomes using long-read sequencing technology.

Currently, different sequencing platforms are used to generate plant genomes and no workflow has been properly developed to optimize time, cost, and assembly quality. We present LeafGo, a complete de novo plant genome workflow, that starts from tissue and produces genomes with modest laboratory and bioinformatic resources in approximately 7 days and using one long-read sequencing technology. LeafGo is optimized with ten different plant species, three of which are used to generate high-quality chromosome-level assemblies without any scaffolding technologies. Finally, we report the diploid genomes of Eucalyptus rudis and E. camaldulensis and the allotetraploid genome of Arachis hypogaea.

Periodic Fever with Aphthous Stomatitis, Pharyngitis, and Cervical Adenitis Syndrome Is Associated with a CARD8 Variant Unable To Bind the NLRP3 Inflammasome.

Periodic fever with aphthous stomatitis, pharyngitis, and cervical adenitis (PFAPA) is a relatively common autoinflammatory condition that primarily affects children. Although tendencies were reported for this syndrome, genetic variations influencing risk and disease progression are poorly understood. In this study, we performed next-generation sequencing for 82 unrelated PFAPA patients and identified a frameshift variant in the CARD8 gene (CARD8-FS). Subsequently, we compared the frequency of CARD8-FS carriers in our PFAPA cohort (13.9%) with a healthy local population group (3.2%) and found a significant association between the CARD8-FS polymorphism and risk for PFAPA syndrome (p = 0.012; odds ratio: 4.96 [95% confidence interval, 1.33-18.47]). Moreover, CARD8-FS carriers display a distinct PFAPA phenotype that is characterized by a higher prevalence of symptoms out of flares and oral aphthosis (both p = 0.02 compared with PFAPA patients without the frameshift variant). CARD8 encodes a protein component of the NLRP3 inflammasome, which plays an important role in inflammation and contributes to the pathology of various autoinflammatory diseases. We found that the CARD8-FS variant led to a truncated CARD8 protein lacking the FIIND and CARD domains. As a result, the mutant CARD8 protein lost the ability to interact with the NOD domain of NLRP3. In summary, these results identify a new CARD8 variant associated with PFAPA and further suggest that disruption of the interaction between CARD8 and NLRP3 can regulate autoinflammation in patients.

Novel mutation in the CHST6 gene causes macular corneal dystrophy in a black South African family.

BACKGROUND: Macular corneal dystrophy (MCD) is a rare autosomal recessive disorder that is characterized by progressive corneal opacity that starts in early childhood and ultimately progresses to blindness in early adulthood. The aim of this study was to identify the cause of MCD in a black South African family with two affected sisters. METHODS: A multigenerational South African Sotho-speaking family with type I MCD was studied using whole exome sequencing. Variant filtering to identify the MCD-causal mutation included the disease inheritance pattern, variant minor allele frequency and potential functional impact. RESULTS: Ophthalmologic evaluation of the cases revealed a typical MCD phenotype and none of the other family members were affected. An average of 127 713 variants per individual was identified following exome sequencing and approximately 1.2 % were not present in any of the investigated public databases. Variant filtering identified a homozygous E71Q mutation in CHST6, a known MCD-causing gene encoding corneal N-acetyl glucosamine-6-O-sulfotransferase. This E71Q mutation results in a non-conservative amino acid change in a highly conserved functional domain of the human CHST6 that is essential for enzyme activity. CONCLUSION: We identified a novel E71Q mutation in CHST6 as the MCD-causal mutation in a black South African family with type I MCD. This is the first description of MCD in a black Sub-Saharan African family and therefore contributes valuable insights into the genetic aetiology of this disease, while improving genetic counselling for this and potentially other MCD families.

Heart structure-specific transcriptomic atlas reveals conserved microRNA-mRNA interactions.

MicroRNAs are short non-coding RNAs that regulate gene expression at the post-transcriptional level and play key roles in heart development and cardiovascular diseases. Here, we have characterized the expression and distribution of microRNAs across eight cardiac structures (left and right ventricles, apex, papillary muscle, septum, left and right atrium and valves) in rat, Beagle dog and cynomolgus monkey using microRNA sequencing. Conserved microRNA signatures enriched in specific heart structures across these species were identified for cardiac valve (miR-let-7c, miR-125b, miR-127, miR-199a-3p, miR-204, miR-320, miR-99b, miR-328 and miR-744) and myocardium (miR-1, miR-133b, miR-133a, miR-208b, miR-30e, miR-499-5p, miR-30e*). The relative abundance of myocardium-enriched (miR-1) and valve-enriched (miR-125b-5p and miR-204) microRNAs was confirmed using in situ hybridization. MicroRNA-mRNA interactions potentially relevant for cardiac functions were explored using anti-correlation expression analysis and microRNA target prediction algorithms. Interactions between miR-1/Timp3, miR-125b/Rbm24, miR-204/Tgfbr2 and miR-208b/Csnk2a2 were identified and experimentally investigated in human pulmonary smooth muscle cells and luciferase reporter assays. In conclusion, we have generated a high-resolution heart structure-specific mRNA/microRNA expression atlas for three mammalian species that provides a novel resource for investigating novel microRNA regulatory circuits involved in cardiac molecular physiopathology.

BarraCUDA - a fast short read sequence aligner using graphics processing units.

BACKGROUND: With the maturation of next-generation DNA sequencing (NGS) technologies, the throughput of DNA sequencing reads has soared to over 600 gigabases from a single instrument run. General purpose computing on graphics processing units (GPGPU), extracts the computing power from hundreds of parallel stream processors within graphics processing cores and provides a cost-effective and energy efficient alternative to traditional high-performance computing (HPC) clusters. In this article, we describe the implementation of BarraCUDA, a GPGPU sequence alignment software that is based on BWA, to accelerate the alignment of sequencing reads generated by these instruments to a reference DNA sequence. FINDINGS: Using the NVIDIA Compute Unified Device Architecture (CUDA) software development environment, we ported the most computational-intensive alignment component of BWA to GPU to take advantage of the massive parallelism. As a result, BarraCUDA offers a magnitude of performance boost in alignment throughput when compared to a CPU core while delivering the same level of alignment fidelity. The software is also capable of supporting multiple CUDA devices in parallel to further accelerate the alignment throughput. CONCLUSIONS: BarraCUDA is designed to take advantage of the parallelism of GPU to accelerate the alignment of millions of sequencing reads generated by NGS instruments. By doing this, we could, at least in part streamline the current bioinformatics pipeline such that the wider scientific community could benefit from the sequencing technology.BarraCUDA is currently available from http://seqbarracuda.sf.net.

Systematic bias in high-throughput sequencing data and its correction by BEADS.

Genomic sequences obtained through high-throughput sequencing are not uniformly distributed across the genome. For example, sequencing data of total genomic DNA show significant, yet unexpected enrichments on promoters and exons. This systematic bias is a particular problem for techniques such as chromatin immunoprecipitation, where the signal for a target factor is plotted across genomic features. We have focused on data obtained from Illumina's Genome Analyser platform, where at least three factors contribute to sequence bias: GC content, mappability of sequencing reads, and regional biases that might be generated by local structure. We show that relying on input control as a normalizer is not generally appropriate due to sample to sample variation in bias. To correct sequence bias, we present BEADS (bias elimination algorithm for deep sequencing), a simple three-step normalization scheme that successfully unmasks real binding patterns in ChIP-seq data. We suggest that this procedure be done routinely prior to data interpretation and downstream analyses.

An assessment of histone-modification antibody quality.

We have tested the specificity and utility of more than 200 antibodies raised against 57 different histone modifications in Drosophila melanogaster, Caenorhabditis elegans and human cells. Although most antibodies performed well, more than 25% failed specificity tests by dot blot or western blot. Among specific antibodies, more than 20% failed in chromatin immunoprecipitation experiments. We advise rigorous testing of histone-modification antibodies before use, and we provide a website for posting new test results (http://compbio.med.harvard.edu/antibodies/).

Broad chromosomal domains of histone modification patterns in C. elegans.

Chromatin immunoprecipitation identifies specific interactions between genomic DNA and proteins, advancing our understanding of gene-level and chromosome-level regulation. Based on chromatin immunoprecipitation experiments using validated antibodies, we define the genome-wide distributions of 19 histone modifications, one histone variant, and eight chromatin-associated proteins in Caenorhabditis elegans embryos and L3 larvae. Cluster analysis identified five groups of chromatin marks with shared features: Two groups correlate with gene repression, two with gene activation, and one with the X chromosome. The X chromosome displays numerous unique properties, including enrichment of monomethylated H4K20 and H3K27, which correlate with the different repressive mechanisms that operate in somatic tissues and germ cells, respectively. The data also revealed striking differences in chromatin composition between the autosomes and between chromosome arms and centers. Chromosomes I and III are globally enriched for marks of active genes, consistent with containing more highly expressed genes, compared to chromosomes II, IV, and especially V. Consistent with the absence of cytological heterochromatin and the holocentric nature of C. elegans chromosomes, markers of heterochromatin such as H3K9 methylation are not concentrated at a single region on each chromosome. Instead, H3K9 methylation is enriched on chromosome arms, coincident with zones of elevated meiotic recombination. Active genes in chromosome arms and centers have very similar histone mark distributions, suggesting that active domains in the arms are interspersed with heterochromatin-like structure. These data, which confirm and extend previous studies, allow for in-depth analysis of the organization and deployment of the C. elegans genome during development.

Integrative analysis of the Caenorhabditis elegans genome by the modENCODE project.

We systematically generated large-scale data sets to improve genome annotation for the nematode Caenorhabditis elegans, a key model organism. These data sets include transcriptome profiling across a developmental time course, genome-wide identification of transcription factor-binding sites, and maps of chromatin organization. From this, we created more complete and accurate gene models, including alternative splice forms and candidate noncoding RNAs. We constructed hierarchical networks of transcription factor-binding and microRNA interactions and discovered chromosomal locations bound by an unusually large number of transcription factors. Different patterns of chromatin composition and histone modification were revealed between chromosome arms and centers, with similarly prominent differences between autosomes and the X chromosome. Integrating data types, we built statistical models relating chromatin, transcription factor binding, and gene expression. Overall, our analyses ascribed putative functions to most of the conserved genome.

DANGLE: A Bayesian inferential method for predicting protein backbone dihedral angles and secondary structure.

This paper introduces DANGLE, a new algorithm that employs Bayesian inference to estimate the likelihood of all possible values of the backbone dihedral angles phi and psi for each residue in a query protein, based on observed chemical shifts and the conformational preferences of each amino acid type. The method provides robust estimates of phi and psi within realistic boundary ranges, an indication of the degeneracy in the relationship between shift measurements and conformation at each site, and faithful secondary structure state assignments. When a simple degeneracy-based filtering procedure is applied, DANGLE offers an ideal compromise between accuracy and coverage when compared with other shift-based dihedral angle prediction methods. In addition, per residue analysis of shift/structure degeneracy has potential to be a useful new approach for studying the properties of unfolded proteins, with sufficient sensitivity to identify regions of residual structure in the acid denatured state of apomyoglobin.