The Norway spruce genome sequence and conifer genome evolution.

Conifers have dominated forests for more than 200 million years and are of huge ecological and economic importance. Here we present the draft assembly of the 20-gigabase genome of Norway spruce (Picea abies), the first available for any gymnosperm. The number of well-supported genes (28,354) is similar to the >100 times smaller genome of Arabidopsis thaliana, and there is no evidence of a recent whole-genome duplication in the gymnosperm lineage. Instead, the large genome size seems to result from the slow and steady accumulation of a diverse set of long-terminal repeat transposable elements, possibly owing to the lack of an efficient elimination mechanism. Comparative sequencing of Pinus sylvestris, Abies sibirica, Juniperus communis, Taxus baccata and Gnetum gnemon reveals that the transposable element diversity is shared among extant conifers. Expression of 24-nucleotide small RNAs, previously implicated in transposable element silencing, is tissue-specific and much lower than in other plants. We further identify numerous long (>10,000 base pairs) introns, gene-like fragments, uncharacterized long non-coding RNAs and short RNAs. This opens up new genomic avenues for conifer forestry and breeding.

Population-scale sequencing reveals genetic differentiation due to local adaptation in Atlantic herring.

The Atlantic herring (Clupea harengus), one of the most abundant marine fishes in the world, has historically been a critical food source in Northern Europe. It is one of the few marine species that can reproduce throughout the brackish salinity gradient of the Baltic Sea. Previous studies based on few genetic markers have revealed a conspicuous lack of genetic differentiation between geographic regions, consistent with huge population sizes and minute genetic drift. Here, we present a cost-effective genome-wide study in a species that lacks a genome sequence. We first assembled a muscle transcriptome and then aligned genomic reads to the transcripts, creating an "exome assembly," capturing both exons and flanking sequences. We then resequenced pools of fish from a wide geographic range, including the Northeast Atlantic, as well as different regions in the Baltic Sea, aligned the reads to the exome assembly, and identified 440,817 SNPs. The great majority of SNPs showed no appreciable differences in allele frequency among populations; however, several thousand SNPs showed striking differences, some approaching fixation for different alleles. The contrast between low genetic differentiation at most loci and striking differences at others implies that the latter category primarily reflects natural selection. A simulation study confirmed that the distribution of the fixation index F(ST) deviated significantly from expectation for selectively neutral loci. This study provides insights concerning the population structure of an important marine fish and establishes the Atlantic herring as a model for population genetic studies of adaptation and natural selection.

Characterization of the viral microbiome in patients with severe lower respiratory tract infections, using metagenomic sequencing.

The human respiratory tract is heavily exposed to microorganisms. Viral respiratory tract pathogens, like RSV, influenza and rhinoviruses cause major morbidity and mortality from respiratory tract disease. Furthermore, as viruses have limited means of transmission, viruses that cause pathogenicity in other tissues may be transmitted through the respiratory tract. It is therefore important to chart the human virome in this compartment. We have studied nasopharyngeal aspirate samples submitted to the Karolinska University Laboratory, Stockholm, Sweden from March 2004 to May 2005 for diagnosis of respiratory tract infections. We have used a metagenomic sequencing strategy to characterize viruses, as this provides the most unbiased view of the samples. Virus enrichment followed by 454 sequencing resulted in totally 703,790 reads and 110,931 of these were found to be of viral origin by using an automated classification pipeline. The snapshot of the respiratory tract virome of these 210 patients revealed 39 species and many more strains of viruses. Most of the viral sequences were classified into one of three major families; Paramyxoviridae, Picornaviridae or Orthomyxoviridae. The study also identified one novel type of Rhinovirus C, and identified a number of previously undescribed viral genetic fragments of unknown origin.

Total RNA sequencing reveals nascent transcription and widespread co-transcriptional splicing in the human brain.

Transcriptome sequencing allows for analysis of mature RNAs at base pair resolution. Here we show that RNA-seq can also be used for studying nascent RNAs undergoing transcription. We sequenced total RNA from human brain and liver and found a large fraction of reads (up to 40%) within introns. Intronic RNAs were abundant in brain tissue, particularly for genes involved in axonal growth and synaptic transmission. Moreover, we detected significant differences in intronic RNA levels between fetal and adult brains. We show that the pattern of intronic sequence read coverage is explained by nascent transcription in combination with co-transcriptional splicing. Further analysis of co-transcriptional splicing indicates a correlation between slowly removed introns and alternative splicing. Our data show that sequencing of total RNA provides unique insight into the transcriptional processes in the cell, with particular importance for normal brain development.

Identification of novel exons and transcribed regions by chimpanzee transcriptome sequencing.

BACKGROUND: We profile the chimpanzee transcriptome by using deep sequencing of cDNA from brain and liver, aiming to quantify expression of known genes and to identify novel transcribed regions. RESULTS: Using stringent criteria for transcription, we identify 12,843 expressed genes, with a majority being found in both tissues. We further identify 9,826 novel transcribed regions that are not overlapping with annotated exons, mRNAs or ESTs. Over 80% of the novel transcribed regions map within or in the vicinity of known genes, and by combining sequencing data with de novo splice predictions we predict several of the novel transcribed regions to be new exons or 3' UTRs. For approximately 350 novel transcribed regions, the corresponding DNA sequence is absent in the human reference genome. The presence of novel transcribed regions in five genes and in one intergenic region is further validated with RT-PCR. Finally, we describe and experimentally validate a putative novel multi-exon gene that belongs to the ATP-cassette transporter gene family. This gene does not appear to be functional in human since one exon is absent from the human genome. In addition to novel exons and UTRs, novel transcribed regions may also stem from different types of noncoding transcripts. We note that expressed repeats and introns from unspliced mRNAs are especially common in our data. CONCLUSIONS: Our results extend the chimpanzee gene catalogue with a large number of novel exons and 3' UTRs and thus support the view that mammalian gene annotations are not yet complete.

Global and unbiased detection of splice junctions from RNA-seq data.

We have developed a new strategy for de novo prediction of splice junctions in short-read RNA-seq data, suitable for detection of novel splicing events and chimeric transcripts. When tested on mouse RNA-seq data, >31,000 splice events were predicted, of which 88% bridged between two regions separated by

Genome-wide analysis of chimpanzee genes with premature termination codons.

BACKGROUND: Premature termination codons (PTCs) cause mRNA degradation or a truncated protein and thereby contribute to the transcriptome and proteome divergence between species. Here we present the first genome-wide study of PTCs in the chimpanzee. By comparing the human and chimpanzee genome sequences we identify and characterize genes with PTCs, in order to understand the contribution of these mutations to the transcriptome diversity between the species. RESULTS: We have studied a total of 13,487 human-chimpanzee gene pairs and found that ~8% were affected by PTCs in the chimpanzee. A majority (764/1,109) of PTCs were caused by insertions or deletions and the remaining part was caused by substitutions. The distribution of PTC genes varied between chromosomes, with Y having the highest proportion. Furthermore, the density of PTC genes varied on a megabasepair scale within chromosomes and we found the density to be correlated both with indel divergence and proximity to the telomere. Within genes, PTCs were more common close to the 5' and 3' ends of the amino acid sequence. Gene Ontology classification revealed that olfactory receptor genes were over represented among the PTC genes. CONCLUSION: Our results showed that the density of PTC genes fluctuated across the genome depending on the local genomic context. PTCs were preferentially located in the terminal parts of the transcript, which generally have a lower frequency of functional domains, indicating that selection was operating against PTCs at sites central to protein function. The enrichment of GO terms associated with olfaction suggests that PTCs may have influenced the difference in the repertoire of olfactory genes between humans and chimpanzees. In summary, 8% of the chimpanzee genes were affected by PTCs and this type of variation is likely to have an important effect on the transcript and proteomic divergence between humans and chimpanzees.

Comparative genomic analysis of human and chimpanzee indicates a key role for indels in primate evolution.

Sequence comparison of humans and chimpanzees is of interest to understand the mechanisms behind primate evolution. Here we present an independent analysis of human chromosome 21 and the high-quality BAC clone sequences of the homologous chimpanzee chromosome 22. In contrast to previous studies, we have used global alignment methods and Ensembl predictions of protein coding genes (n = 224) for the analysis. Divergence due to insertions and deletions (indels) along with substitutions was examined separately for different genomic features (coding, noncoding genic, and intergenic sequence). The major part of the genomic divergence could be attributed to indels (5.07%), while the nucleotide divergence was estimated as 1.52%. Thus the total divergence was estimated as 6.58%. When excluding repeats and low-complexity DNA the total divergence decreased to 2.37%. The chromosomal distribution of nucleotide substitutions and indel events was significantly correlated. To further examine the role of indels in primate evolution we focused on coding sequences. Indels were found within the coding sequence of 13% of the genes and approximately half of the indels have not been reported previously. In 5% of the chimpanzee genes, indels or substitutions caused premature stop codons that rendered the affected transcripts nonfunctional. Taken together, our findings demonstrate that indels comprise the majority of the genomic divergence. Furthermore, indels occur frequently in coding sequences. Our results thereby support the hypothesis that indels may have a key role in primate evolution.