The Genomic Scrapheap Challenge; Extracting Relevant Data from Unmapped Whole Genome Sequencing Reads, Including Strain Specific Genomic Segments, in Rats.

Unmapped next-generation sequencing reads are typically ignored while they contain biologically relevant information. We systematically analyzed unmapped reads from whole genome sequencing of 33 inbred rat strains. High quality reads were selected and enriched for biologically relevant sequences; similarity-based analysis revealed clustering similar to previously reported phylogenetic trees. Our results demonstrate that on average 20% of all unmapped reads harbor sequences that can be used to improve reference genomes and generate hypotheses on potential genotype-phenotype relationships. Analysis pipelines would benefit from incorporating the described methods and reference genomes would benefit from inclusion of the genomic segments obtained through these efforts.

Genomic landscape of rat strain and substrain variation.

BACKGROUND: Since the completion of the rat reference genome in 2003, whole-genome sequencing data from more than 40 rat strains have become available. These data represent the broad range of strains that are used in rat research including commonly used substrains. Currently, this wealth of information cannot be used to its full extent, because the variety of different variant calling algorithms employed by different groups impairs comparison between strains. In addition, all rat whole genome sequencing studies to date used an outdated reference genome for analysis (RGSC3.4 released in 2004). RESULTS: Here we present a comprehensive, multi-sample and uniformly called set of genetic variants in 40 rat strains, including 19 substrains. We reanalyzed all primary data using a recent version of the rat reference assembly (RGSC5.0 released in 2012) and identified over 12 million genomic variants (SNVs, indels and structural variants) among the 40 strains. 28,318 SNVs are specific to individual substrains, which may be explained by introgression from other unsequenced strains and ongoing evolution by genetic drift. Substrain SNVs may have a larger predicted functional impact compared to older shared SNVs. CONCLUSIONS: In summary we present a comprehensive catalog of uniformly analyzed genetic variants among 40 widely used rat inbred strains based on the RGSC5.0 assembly. This represents a valuable resource, which will facilitate rat functional genomic research. In line with previous observations, our genome-wide analyses do not show evidence for contribution of multiple ancestral founder rat subspecies to the currently used rat inbred strains, as is the case for mouse. In addition, we find that the degree of substrain variation is highly variable between strains, which is of importance for the correct interpretation of experimental data from different labs.

Lack of major genome instability in tumors of p53 null rats.

Tumorigenesis is often associated with loss of tumor suppressor genes (such as TP53), genomic instability and telomere lengthening. Previously, we generated and characterized a rat p53 knockout model in which the homozygous rats predominantly develop hemangiosarcomas whereas the heterozygous rats mainly develop osteosarcomas. Using genome-wide analyses, we find that the tumors that arise in the heterozygous and homozygous Tp53C273X mutant animals are also different in their genomic instability profiles. While p53 was fully inactivated in both heterozygous and homozygous knockout rats, tumors from homozygous animals show very limited aneuploidy and low degrees of somatic copy number variation as compared to the tumors from heterozygous animals. In addition, complex structural rearrangements such as chromothripsis and breakage-fusion-bridge cycles were never found in tumors from homozygous animals, while these were readily detectable in tumors from heterozygous animals. Finally, we measured telomere length and telomere lengthening pathway activity and found that tumors of homozygous animals have longer telomeres but do not show clear telomerase or alternative lengthening of telomeres (ALT) activity differences as compared to the tumors from heterozygous animals. Taken together, our results demonstrate that host p53 status in this rat p53 knockout model has a large effect on both tumor type and genomic instability characteristics, where full loss of functional p53 is not the main driver of large-scale structural variations. Our results also suggest that chromothripsis primarily occurs under p53 heterozygous rather than p53 null conditions.

Combined sequence-based and genetic mapping analysis of complex traits in outbred rats.

Genetic mapping on fully sequenced individuals is transforming understanding of the relationship between molecular variation and variation in complex traits. Here we report a combined sequence and genetic mapping analysis in outbred rats that maps 355 quantitative trait loci for 122 phenotypes. We identify 35 causal genes involved in 31 phenotypes, implicating new genes in models of anxiety, heart disease and multiple sclerosis. The relationship between sequence and genetic variation is unexpectedly complex: at approximately 40% of quantitative trait loci, a single sequence variant cannot account for the phenotypic effect. Using comparable sequence and mapping data from mice, we show that the extent and spatial pattern of variation in inbred rats differ substantially from those of inbred mice and that the genetic variants in orthologous genes rarely contribute to the same phenotype in both species.

Homozygous and heterozygous p53 knockout rats develop metastasizing sarcomas with high frequency.

The TP53 tumor suppressor gene is mutated in the majority of human cancers. Inactivation of p53 in a variety of animal models results in early-onset tumorigenesis, reflecting the importance of p53 as a gatekeeper tumor suppressor. We generated a mutant Tp53 allele in the rat using a target-selected mutagenesis approach. Here, we report that homozygosity for this allele results in complete loss of p53 function. Homozygous mutant rats predominantly develop sarcomas with an onset of 4 months of age with a high occurrence of pulmonary metastases. Heterozygous rats develop sarcomas starting at 8 months of age. Molecular analysis revealed that these tumors exhibit a loss-of-heterozygosity of the wild-type Tp53 allele. These unique features make this rat highly complementary to other rodent p53 knockout models and a versatile tool for investigating tumorigenesis processes as well as genotoxic studies.

Functional CD1d and/or NKT cell invariant chain transcript in horse, pig, African elephant and guinea pig, but not in ruminants.

CD1d-restricted invariant natural killer T cells (NKT cells) have been well characterized in humans and mice, but it is unknown whether they are present in other species. Here we describe the invariant TCR alpha chain and the full length CD1d transcript of pig and horse. Molecular modeling predicts that porcine (po) invariant TCR alpha chain/poCD1d/alpha-GalCer and equine (eq) invariant TCR alpha chain/eqCD1d/alpha-GalCer form complexes that are highly homologous to the human complex. Since a prerequisite for the presence of NKT cells is the expression of CD1d protein, we performed searches for CD1D genes and CD1d transcripts in multiple species. Previously, cattle and guinea pig have been suggested to lack CD1D genes. The CD1D genes of European taurine cattle (Bos taurus) are known to be pseudogenes because of disrupting mutations in the start codon and in the donor splice site of the first intron. Here we show that the same mutations are found in six other ruminants: African buffalo, sheep, bushbuck, bongo, N'Dama cattle, and roe deer. In contrast, intact CD1d transcripts were found in guinea pig, African elephant, horse, rabbit, and pig. Despite the discovery of a highly homologous NKT/CD1d system in pig and horse, our data suggest that functional CD1D and CD1d-restricted NKT cells are not universally present in mammals.