Motivation Accurate alignment of high-throughput RNA-seq data is a challenging and yet unsolved problem because of the non-contiguous transcript structure, relatively short read lengths and constantly increasing throughput of the sequencing technologies. Currently available RNA-seq aligners suffer from high mapping error rates, low mapping speed, read length limitation and mapping biases. Results To align our large (>80 billon reads) ENCODE Transcriptome RNA-seq dataset, we developed the Spliced Transcripts Alignment to a Reference (STAR) software based on a previously undescribed RNA-seq alignment algorithm that uses sequential maximum mappable seed search in uncompressed suffix arrays followed by seed clustering and stitching procedure. STAR outperforms other aligners by a factor of >50 in mapping speed, aligning to the human genome 550 million 2 × 76 bp paired-end reads per hour on a modest 12-core server, while at the same time improving alignment sensitivity and precision. In addition to unbiased de novo detection of canonical junctions, STAR can discover non-canonical splices and chimeric (fusion) transcripts, and is also capable of mapping full-length RNA sequences. Using Roche 454 sequencing of reverse transcription polymerase chain reaction amplicons, we experimentally validated 1960 novel intergenic splice junctions with an 80-90% success rate, corroborating the high precision of the STAR mapping strategy. Availability and implementation STAR is implemented as a standalone C++ code. STAR is free open source software distributed under GPLv3 license and can be downloaded from http://code.google.com/p/rna-star/.

/pdf/star-ultrafast-universal-rna-seq-aligner-r1vlt61g3c.pdf

STAR: ultrafast universal RNA-seq aligner

The human genome holds an extraordinary trove of information about human development, physiology, medicine and evolution. Here we report the results of an international collaboration to produce and make freely available a draft sequence of the human genome. We also present an initial analysis of the data, describing some of the insights that can be gleaned from the sequence.

/pdf/initial-sequencing-and-analysis-of-the-human-genome-5dapk1rsx1.pdf

Initial sequencing and analysis of the human genome.

limma is an R/Bioconductor software package that provides an integrated solution for analysing data from gene expression experiments. It contains rich features for handling complex experimental designs and for information borrowing to overcome the problem of small sample sizes. Over the past decade, limma has been a popular choice for gene discovery through differential expression analyses of microarray and high-throughput PCR data. The package contains particularly strong facilities for reading, normalizing and exploring such data. Recently, the capabilities of limma have been significantly expanded in two important directions. First, the package can now perform both differential expression and differential splicing analyses of RNA sequencing (RNA-seq) data. All the downstream analysis tools previously restricted to microarray data are now available for RNA-seq as well. These capabilities allow users to analyse both RNA-seq and microarray data with very similar pipelines. Second, the package is now able to go past the traditional gene-wise expression analyses in a variety of ways, analysing expression profiles in terms of co-regulated sets of genes or in terms of higher-order expression signatures. This provides enhanced possibilities for biological interpretation of gene expression differences. This article reviews the philosophy and design of the limma package, summarizing both new and historical features, with an emphasis on recent enhancements and features that have not been previously described.

/pdf/limma-powers-differential-expression-analyses-for-rna-4yl7ae3arh.pdf

limma powers differential expression analyses for RNA-sequencing and microarray studies

抗原变异可使得多种致病微生物易于逃避宿主免疫应答。表达在感染红细胞表面的恶性疟原虫红细胞表面蛋白1（PfPMP1）与感染红细胞、内皮细胞、树突状细胞以及胎盘的单个或多个受体作用，在黏附及免疫逃避中起关键的作用。每个单倍体基因组var基因家族编码约60种成员，通过启动转录不同的var基因变异体为抗原变异提供了分子基础。

疟原虫var基因转换速率变化导致抗原变异[英]／Paul H, Robert P, Christodoulou Z, et al//Proc Natl Acad Sci U S A

http://ggdpathway.wustl.edu/files/2014/08/ACMG-Guidelines-2015_Richards-et-al.pdf

Standards and guidelines for the interpretation of sequence variants: a joint consensus recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology.

Amino Acid Substitutions in the Factor VIII Coagulation Protein that Cause Hemophilia A in Man1

The study of transcription using genomic tiling arrays has lead to the identification of numerous additional exons. One example is the MECP2 gene on the X chromosome; using 5'RACE and RT-PCR in human tissues and cell lines, we have found more than 70 novel exons (RACEfrags) connecting to at least one annotated exon.. We sequenced all MECP2-connected exons and flanking sequences in 3 groups: 46 patients with the Rett syndrome and without mutations in the currently annotated exons of the MECP2 and CDKL5 genes; 32 patients with the Rett syndrome and identified mutations in the MECP2 gene; 100 control individuals from the same geoethnic group. Approximately 13 kb were sequenced per sample, (2.4 Mb of DNA resequencing). A total of 75 individuals had novel rare variants (mostly private variants) but no statistically significant difference was found among the 3 groups. These results suggest that variants in the newly discovered exons may not contribute to Rett syndrome. Interestingly however, there are about twice more variants in the novel exons than in the flanking sequences (44 vs. 21 for approximately 1.3 Mb sequenced for each class of sequences, p=0.0025). Thus the evolutionary forces that shape these novel exons may be different than those of neighboring sequences.

/pdf/variation-in-novel-exons-racefrags-of-the-mecp2-gene-in-rett-526j084d8p.pdf

Variation in novel exons (RACEfrags) of the MECP2 gene in Rett syndrome patients and controls.

Mutations in human genetic disease: nature and consequences

Perturbations of Heart Development and Function in Cardiomyocytes from hESC with Trisomy 21

Intellectual disability (ID) is a highly heterogeneous disorder with hundreds of associated genes. Despite progress in the identification of the genetic causes of ID following the introduction of high-throughput sequencing, about half of affected individuals still remain without a molecular diagnosis. Consanguineous families with affected individuals provide a unique opportunity to identify novel recessive causative genes.
In this report we describe a novel autosomal recessive neurodevelopmental disorder. We identified two consanguineous families with homozygous variants predicted to alter the splicing of 	ATP9A
 which encodes a transmembrane lipid flippase of the class II P4-ATPases. The three individuals homozygous for these putatively truncating variants presented with severe ID, motor and speech impairment, and behavioral anomalies. Consistent with a causative role of ATP9A in these patients, a previously described Atp9a-/- mouse model showed behavioral changes.

/pdf/biallelic-truncation-variants-in-atp9a-are-associated-with-a-1c9dw89wx2.pdf

Stylianos E. Antonarakis

Papers

Amino Acid Substitutions in the Factor VIII Coagulation Protein that Cause Hemophilia A in Man1

Variation in novel exons (RACEfrags) of the MECP2 gene in Rett syndrome patients and controls.

Mutations in human genetic disease: nature and consequences

Perturbations of Heart Development and Function in Cardiomyocytes from hESC with Trisomy 21

Biallelic truncation variants in ATP9A are associated with a novel autosomal recessive neurodevelopmental disorder