Genome-wide analysis of trans-splicing in the nematode Pristionchus pacificus unravels conserved gene functions for germline and dauer development in divergent operons.

Discovery of trans-splicing in multiple metazoan lineages led to the identification of operon-like gene organization in diverse organisms, including trypanosomes, tunicates, and nematodes, but the functional significance of such operons is not completely understood. To see whether the content or organization of operons serves similar roles across species, we experimentally defined operons in the nematode model Pristionchus pacificus. We performed affinity capture experiments on mRNA pools to specifically enrich for transcripts that are trans-spliced to either the SL1- or SL2-spliced leader, using spliced leader-specific probes. We obtained distinct trans-splicing patterns from the analysis of three mRNA pools (total mRNA, SL1 and SL2 fraction) by RNA-seq. This information was combined with a genome-wide analysis of gene orientation and spacing. We could confirm 2219 operons by RNA-seq data out of 6709 candidate operons, which were predicted by sequence information alone. Our gene order comparison of the Caenorhabditis elegans and P. pacificus genomes shows major changes in operon organization in the two species. Notably, only 128 out of 1288 operons in C. elegans are conserved in P. pacificus. However, analysis of gene-expression profiles identified conserved functions such as an enrichment of germline-expressed genes and higher expression levels of operonic genes during recovery from dauer arrest in both species. These results provide support for the model that a necessity for increased transcriptional efficiency in the context of certain developmental processes could be a selective constraint for operon evolution in metazoans. Our method is generally applicable to other metazoans to see if similar functional constraints regulate gene organization into operons.

Exome sequencing helped the fine diagnosis of two siblings afflicted with atypical Timothy syndrome (TS2).

BACKGROUND: Long-QT syndrome (LQTS) causes a prolongation of the QT-interval in the ECG leading to life threatening tachyarrhythmia and ventricular fibrillation. One atypical form of LQTS, Timothy syndrome (TS), is associated with syndactyly, immune deficiency, cognitive and neurological abnormalities as well as distinct cranio-facial abnormalities. CASE PRESENTATION: On a family with both children diagnosed with clinical LQTS, we performed whole exome sequencing to comprehensively screen for causative mutations after a targeted candidate gene panel screen for Long-QT syndrome target genes failed to identify any underlying genetic defect. Using exome sequencing, we identified in both affected children, a p.402G > S mutation in exon 8 of the CACNA1C gene, a voltage-dependent Ca2+ channel. The mutation was inherited from their father, a mosaic mutation carrier. Based on this molecular finding and further more careful clinical examination, we refined the diagnosis to be Timothy syndrome (TS2) and thereby were able to present new therapeutic approaches. CONCLUSIONS: Our study highlights the difficulties in accurate diagnosis of patients with rare diseases, especially those with atypical clinical manifestation. Such challenge could be addressed with the help of comprehensive and unbiased mutation screening, such as exome sequencing.

Conserved miRNAs are candidate post-transcriptional regulators of developmental arrest in free-living and parasitic nematodes.

Animal development is complex yet surprisingly robust. Animals may develop alternative phenotypes conditional on environmental changes. Under unfavorable conditions, Caenorhabditis elegans larvae enter the dauer stage, a developmentally arrested, long-lived, and stress-resistant state. Dauer larvae of free-living nematodes and infective larvae of parasitic nematodes share many traits including a conserved endocrine signaling module (DA/DAF-12), which is essential for the formation of dauer and infective larvae. We speculated that conserved post-transcriptional regulatory mechanism might also be involved in executing the dauer and infective larvae fate. We used an unbiased sequencing strategy to characterize the microRNA (miRNA) gene complement in C. elegans, Pristionchus pacificus, and Strongyloides ratti. Our study raised the number of described miRNA genes to 257 for C. elegans, tripled the known gene set for P. pacificus to 362 miRNAs, and is the first to describe miRNAs in a Strongyloides parasite. Moreover, we found a limited core set of 24 conserved miRNA families in all three species. Interestingly, our estimated expression fold changes between dauer versus nondauer stages and infective larvae versus free-living stages reveal that despite the speed of miRNA gene set evolution in nematodes, homologous gene families with conserved "dauer-infective" expression signatures are present. These findings suggest that common post-transcriptional regulatory mechanisms are at work and that the same miRNA families play important roles in developmental arrest and long-term survival in free-living and parasitic nematodes.

Breakpoint analysis of balanced chromosome rearrangements by next-generation paired-end sequencing.

Characterisation of breakpoints in disease-associated balanced chromosome rearrangements (DBCRs), which disrupt or inactivate specific genes, has facilitated the molecular elucidation of a wide variety of genetic disorders. However, conventional methods for mapping chromosome breakpoints, such as in situ hybridisation with fluorescent dye-labelled bacterial artificial chromosome clones (BAC-FISH), are laborious, time consuming and often with insufficient resolution to unequivocally identify the disrupted gene. By combining DNA array hybridisation with chromosome sorting, the efficiency of breakpoint mapping has dramatically improved. However, this can only be applied when the physical properties of the derivative chromosomes allow them to be flow sorted. To characterise the breakpoints in all types of balanced chromosome rearrangements more efficiently and more accurately, we performed massively parallel sequencing using Illumina 1G analyser and ABI SOLiD systems to generate short sequencing reads from both ends of DNA fragments. We applied this method to four different DBCRs, including two reciprocal translocations and two inversions. By identifying read pairs spanning the breakpoints, we were able to map the breakpoints to a region of a few hundred base pairs that could be confirmed by subsequent PCR amplification and Sanger sequencing of the junction fragments. Our results show the feasibility of paired-end sequencing of systematic breakpoint mapping and gene finding in patients with disease-associated chromosome rearrangements.