Viruses cause 10–15% of all human cancers. Massively parallel sequencing has recently proved effective for uncovering novel viruses and virus–tumour associations, but this approach has not yet been applied to comprehensive patient cohorts. Here we screen a diverse landscape of human cancer, encompassing 4,433 tumours and 19 cancer types, for known and novel expressed viruses based on >700 billion transcriptome sequencing reads from The Cancer Genome Atlas Research Network. The resulting map confirms and extends current knowledge. We observe recurrent fusion events, including human papillomavirus insertions in RAD51B and ERBB2. Patterns of coadaptation between host and viral gene expression give clues to papillomavirus oncogene function. Importantly, our analysis argues strongly against viral aetiology in several cancers where this has frequently been proposed. We provide a virus–tumour map of unprecedented scale that constitutes a reference for future studies of tumour-associated viruses using transcriptome sequencing data.

Mucins are proteins that cover and protect epithelial cells and are characterized by domains rich in proline, threonine, and serine that are heavily glycosylated (PTS or mucin domains). Because of their sequence polymorphism, these domains cannot be used for evolutionary analysis. Instead, we have made use of the von Willebrand D (VWD) and SEA domains, typical for mucins. A number of animal genomes were examined for these domains to identify mucin homologues, and domains of the resulting proteins were used in phylogenetic studies. The frog Xenopus tropicalis stands out because the number of gel-forming mucins has markedly increased to at least 25 as compared with 5 for higher animals. Furthermore, the frog Muc2 homologues contain unique PTS domains where cysteines are abundant. This animal also has a unique family of secreted mucin-like proteins with alternating PTS and SEA domains, a type of protein also identified in the fishes. The evolution of the Muc4 mucin seems to have occurred by recruitment of a PTS domain to AMOP, NIDO, and VWD domains from a sushi domain-containing family of proteins present in lower animals, and Xenopus is the most deeply branching animal where a protein similar to the mammalian Muc4 was identified. All transmembrane mucins seem to have appeared in the vertebrate lineage, and the MUC1 mucin is restricted to mammals. In contrast, proteins with properties of the gel-forming mucins were identified also in the starlet sea anemone Nematostella vectensis, demonstrating an early origin of this group of mucins. bioinformatics ͉ von Willebrand domain ͉ SEA domain ͉ protein evolution ͉ mucus

To advance our understanding of development, function and diseases in the kidney glomerulus, we have established and large-scale sequenced cDNA libraries from mouse glomeruli at different stages of development, resulting in a catalogue of 6053 different genes. The glomerular cDNA clones were arrayed and hybridized against a series of labeled targets from isolated glomeruli, non-glomerular kidney tissue, FACS-sorted podocytes and brain capillaries, which identified over 300 glomerular cell-enriched transcripts, some of which were further sublocalized to podocytes, mesangial cells and juxtaglomerular cells by in situ hybridization. For the earliest podocyte marker identified, Foxc2, knockout mice were used to analyze the role of this protein during glomerular development. We show that Foxc2 controls the expression of a distinct set of podocyte genes involved in podocyte differentiation and glomerular basement membrane maturation. The primary podocyte defects also cause abnormal differentiation and organization of the glomerular vascular cells. We surmise that studies on the other novel glomerulus-enriched transcripts identified in this study will provide new insight into glomerular development and pathomechanisms of disease.

An RNA hairpin structure referred to as the iron-responsive element (IRE) and iron regulatory proteins (IRPs) are key players in the control of iron metabolism in animal cells. They regulate translation initiation or mRNA stability, and the IRE is found in a variety of mRNAs, such as those encoding ferritin, transferrin receptor (Tfr), erythroid aminolevulinic acid synthase (eALAS), mitochondrial aconitase (mACO), ferroportin, and divalent metal transporter 1 (DMT1). We have studied the evolution of the IRE by considering all mRNAs previously known to be associated with this structure and by computationally examining its occurrence in a large variety of eukaryotic organisms. More than 100 novel sequences together with ;50 IREs that were previously reported resulted in a comprehensive view of the phylogenetic distribution of this element. A comparison of the different mRNAs shows that the IREs of eALAS and mACO are found in chordates, those of ferroportin and Tfr1 are found in vertebrates, and the IRE of DMT1 is confined to mammals. In contrast, the IRE of ferritin occurs in a majority of metazoa including lower metazoa such as sponges and Nematostella (sea anemone). These findings suggest that the ferritin IRE represents the ancestral version of this type of translational control and that during the evolution of higher animals the IRE structure was adopted by other genes. On the basis of primary sequence comparison between different organisms, we suggest that some of these IREs developed by ''convergent evolution'' through stepwise changes in sequence, rather than by recombination events.

RNases P and MRP are ribonucleoprotein complexes involved in tRNA and rRNA processing, respectively. The RNA subunits of these two enzymes are structurally related to each other and play an essential role in the enzymatic reaction. Both of the RNAs have a highly conserved helical region, P4, which is important in the catalytic reaction. We have used a bioinformatics approach based on conserved elements to computationally analyze available genomic sequences of eukaryotic organisms and have identified a large number of novel nuclear RNase P and MRP RNA genes. For MRP RNA for instance, this investigation increases the number of known sequences by a factor of three. We present secondary structure models of many of the predicted RNAs. Although all sequences are able to fold into the consensus secondary structure of P and MRP RNAs, a striking variation in size is observed, ranging from a Nosema locustae MRP RNA of 160 nt to much larger RNAs, e.g. a Plasmodium knowlesi P RNA of 696 nt. The P and MRP RNA genes appear in tandem in some protists, further emphasizing the close evolutionary relationship of these RNAs.