David E. Symer

The bonobo genome compared with the chimpanzee and human genomes

Two African apes are the closest living relatives of humans: the chimpanzee (Pan troglodytes) and the bonobo (Pan paniscus). Although they are similar in many respects, bonobos and chimpanzees differ strikingly in key social and sexual behaviours, and for some of these traits they show more similarity with humans than with each other. Here we report the sequencing and assembly of the bonobo genome to study its evolutionary relationship with the chimpanzee and human genomes. We find that more than three per cent of the human genome is more closely related to either the bonobo or the chimpanzee genome than these are to each other. These regions allow various aspects of the ancestry of the two ape species to be reconstructed. In addition, many of the regions that overlap genes may eventually help us understand the genetic basis of phenotypes that humans share with one of the two apes to the exclusion of the other.

Mutations in U4atac snRNA, a Component of the Minor Spliceosome, in the Developmental Disorder MOPD I

Minor RNA splicing defects can cause a major human developmental disorder. Small nuclear RNAs (snRNAs) are essential factors in messenger RNA splicing. By means of homozygosity mapping and deep sequencing, we show that a gene encoding U4atac snRNA, a component of the minor U12-dependent spliceosome, is mutated in individuals with microcephalic osteodysplastic primordial dwarfism type I (MOPD I), a severe developmental disorder characterized by extreme intrauterine growth retardation and multiple organ abnormalities. Functional assays showed that mutations (30G>A, 51G>A, 55G>A, and 111G>A) associated with MOPD I cause defective U12-dependent splicing. Endogenous U12-dependent but not U2-dependent introns were found to be poorly spliced in MOPD I patient fibroblast cells. The introduction of wild-type U4atac snRNA into MOPD I cells enhanced U12-dependent splicing. These results illustrate the critical role of minor intron splicing in human development.

Lysine methylation of HIV-1 Tat regulates transcriptional activity of the viral LTR

BackgroundThe rate of transcription of the HIV-1 viral genome is mediated by the interaction of the viral protein Tat with the LTR and other transcriptional machinery. These specific interactions can be affected by the state of post-translational modifications on Tat. Previously, we have shown that Tat can be phosphorylated and acetylated in vivo resulting in an increase in the rate of transcription. In the present study, we investigated whether Tat could be methylated on lysine residues, specifically on lysine 50 and 51, and whether this modification resulted in a decrease of viral transcription from the LTR.ResultsWe analyzed the association of Tat with histone methyltransferases of the SUV39-family of SET domain containing proteins in vitro. Tat was found to associate with both SETDB1 and SETDB2, two enzymes which exhibit methyltransferase activity. siRNA against SETDB1 transfected into cell systems with both transient and integrated LTR reporter genes resulted in an increase in transcription of the HIV-LTR in the presence of suboptimal levels of Tat. In vitro methylation assays with Tat peptides containing point mutations at lysines 50 and 51 showed an increased incorporation of methyl groups on lysine 51, however, both residues indicated susceptibility for methylation.ConclusionThe association of Tat with histone methyltransferases and the ability for Tat to be methylated suggests an interesting mechanism of transcriptional regulation through the recruitment of chromatin remodeling proteins to the HIV-1 promoter.

The methyltransferase Setdb2 mediates virus-induced susceptibility to bacterial superinfection

Immune responses are tightly regulated to ensure efficient pathogen clearance while avoiding tissue damage. Here we report that Setdb2 was the only protein lysine methyltransferase induced during infection with influenza virus. Setdb2 expression depended on signaling via type I interferons, and Setdb2 repressed expression of the gene encoding the neutrophil attractant CXCL1 and other genes that are targets of the transcription factor NF-κB. This coincided with occupancy by Setdb2 at the Cxcl1 promoter, which in the absence of Setdb2 displayed diminished trimethylation of histone H3 Lys9 (H3K9me3). Mice with a hypomorphic gene-trap construct of Setdb2 exhibited increased infiltration of neutrophils during sterile lung inflammation and were less sensitive to bacterial superinfection after infection with influenza virus. This suggested that a Setdb2-mediated regulatory crosstalk between the type I interferons and NF-κB pathways represents an important mechanism for virus-induced susceptibility to bacterial superinfection.

Multiple functional variants in long-range enhancer elements contribute to the risk of SNP rs965513 in thyroid cancer

Significance Papillary thyroid carcinoma (PTC) displays a strong hereditary component that is, in part, due to the additive effects of numerous low-penetrance genes or variants, but virtually no mechanistic information is available. Here, we studied a well-known low-penetrance variant (SNP rs965513) located in a region devoid of coding genes. We show that at least four variants located in the immediate vicinity of rs965513 reside in enhancer elements that bind to the promoter region shared by two adjacent thyroid-related genes, forkhead box E1 (FOXE1) and PTC susceptibility candidate 2 (PTCSC2), regulating their expression. The role of intergenic regulatory variants in cancer predisposition and carcinogenesis is growing. Further mechanistic understanding of how these variants work such as described here needs to be acquired. The [A] allele of SNP rs965513 in 9q22 has been consistently shown to be highly associated with increased papillary thyroid cancer (PTC) risk with an odds ratio of ∼1.8 as determined by genome-wide association studies, yet the molecular mechanisms remain poorly understood. Previously, we noted that the expression of two genes in the region, forkhead box E1 (FOXE1) and PTC susceptibility candidate 2 (PTCSC2), is regulated by rs965513 in unaffected thyroid tissue, but the underlying mechanisms were not elucidated. Here, we fine-mapped the 9q22 region in PTC and controls and detected an ∼33-kb linkage disequilibrium block (containing the lead SNP rs965513) that significantly associates with PTC risk. Chromatin characteristics and regulatory element signatures in this block disclosed at least three regulatory elements functioning as enhancers. These enhancers harbor at least four SNPs (rs7864322, rs12352658, rs7847449, and rs10759944) that serve as functional variants. The variant genotypes are associated with differential enhancer activities and/or transcription factor binding activities. Using the chromosome conformation capture methodology, long-range looping interactions of these elements with the promoter region shared by FOXE1 and PTCSC2 in a human papillary thyroid carcinoma cell line (KTC-1) and unaffected thyroid tissue were found. Our results suggest that multiple variants coinherited with the lead SNP and located in long-range enhancers are involved in the transcriptional regulation of FOXE1 and PTCSC2 expression. These results explain the mechanism by which the risk allele of rs965513 predisposes to thyroid cancer.

SRGAP1 Is a Candidate Gene for Papillary Thyroid Carcinoma Susceptibility

BACKGROUND Papillary thyroid carcinoma (PTC) shows high heritability, yet efforts to find predisposing genes have been largely negative. OBJECTIVES The objective of this study was to identify susceptibility genes for PTC. METHODS A genome-wide linkage analysis was performed in 38 families. Targeted association study and screening were performed in 2 large cohorts of PTC patients and controls. Candidate DNA variants were tested in functional studies. RESULTS Linkage analysis and association studies identified the Slit-Robo Rho GTPase activating protein 1 gene (SRGAP1) in the linkage peak as a candidate gene. Two missense variants, Q149H and A275T, localized in the Fes/CIP4 homology domain segregated with the disease in 1 family each. One missense variant, R617C, located in the RhoGAP domain occurred in 1 family. Biochemical assays demonstrated that the ability to inactivate CDC42, a key function of SRGAP1, was severely impaired by the Q149H and R617C variants. CONCLUSIONS Our findings suggest that SRGAP1 is a candidate gene in PTC susceptibility. SRGAP1 is likely a low-penetrant gene, possibly of a modifier type.

An antisense promoter in mouse L1 retrotransposon open reading frame-1 initiates expression of diverse fusion transcripts and limits retrotransposition

Between 6 and 30% of human and mouse transcripts are initiated from transposable elements. However, the promoters driving such transcriptional activity are mostly unknown. We experimentally characterized an antisense (AS) promoter in mouse L1 retrotransposons for the first time, oriented antiparallel to the coding strand of L1 open reading frame-1. We found that AS transcription is mediated by RNA polymerase II. Rapid amplification of cDNA ends cloning mapped transcription start sites adjacent to the AS promoter. We identified >100 novel fusion transcripts, of which many were conserved across divergent mouse lineages, suggesting conservation of potential functions. To evaluate whether AS L1 transcription could regulate L1 retrotransposition, we replaced portions of native open reading frame-1 in donor elements by synonymously recoded sequences. The resulting L1 elements lacked AS promoter activity and retrotransposed more frequently than endogenous L1s. Overexpression of AS L1 transcripts also reduced L1 retrotransposition. This suppression of retrotransposition was largely independent of Dicer. Our experiments shed new light on how AS fusion transcripts are initiated from endogenous L1 elements across the mouse genome. Such AS transcription can contribute substantially both to natural transcriptional variation and to endogenous regulation of L1 retrotransposition.

Merkel Cell Polyomavirus Exhibits Dominant Control of the Tumor Genome and Transcriptome in Virus-Associated Merkel Cell Carcinoma

ABSTRACT Merkel cell polyomavirus is the primary etiological agent of the aggressive skin cancer Merkel cell carcinoma (MCC). Recent studies have revealed that UV radiation is the primary mechanism for somatic mutagenesis in nonviral forms of MCC. Here, we analyze the whole transcriptomes and genomes of primary MCC tumors. Our study reveals that virus-associated tumors have minimally altered genomes compared to non-virus-associated tumors, which are dominated by UV-mediated mutations. Although virus-associated tumors contain relatively small mutation burdens, they exhibit a distinct mutation signature with observable transcriptionally biased kataegic events. In addition, viral integration sites overlap focal genome amplifications in virus-associated tumors, suggesting a potential mechanism for these events. Collectively, our studies indicate that Merkel cell polyomavirus is capable of hijacking cellular processes and driving tumorigenesis to the same severity as tens of thousands of somatic genome alterations. IMPORTANCE A variety of mutagenic processes that shape the evolution of tumors are critical determinants of disease outcome. Here, we sequenced the entire genome of virus-positive and virus-negative primary Merkel cell carcinomas (MCCs), revealing distinct mutation spectra and corresponding expression profiles. Our studies highlight the strong effect that Merkel cell polyomavirus has on the divergent development of viral MCC compared to the somatic alterations that typically drive nonviral tumorigenesis. A more comprehensive understanding of the distinct mutagenic processes operative in viral and nonviral MCCs has implications for the effective treatment of these tumors. A variety of mutagenic processes that shape the evolution of tumors are critical determinants of disease outcome. Here, we sequenced the entire genome of virus-positive and virus-negative primary Merkel cell carcinomas (MCCs), revealing distinct mutation spectra and corresponding expression profiles. Our studies highlight the strong effect that Merkel cell polyomavirus has on the divergent development of viral MCC compared to the somatic alterations that typically drive nonviral tumorigenesis. A more comprehensive understanding of the distinct mutagenic processes operative in viral and nonviral MCCs has implications for the effective treatment of these tumors.

Ultra-rare mutation in long-range enhancer predisposes to thyroid carcinoma with high penetrance.

Thyroid cancer shows high heritability but causative genes remain largely unknown. According to a common hypothesis the genetic predisposition to thyroid cancer is highly heterogeneous; being in part due to many different rare alleles. Here we used linkage analysis and targeted deep sequencing to detect a novel single-nucleotide mutation in chromosome 4q32 (4q32A>C) in a large pedigree displaying non-medullary thyroid carcinoma (NMTC). This mutation is generally ultra-rare; it was not found in 38 NMTC families, in 2676 sporadic NMTC cases or 2470 controls. The mutation is located in a long-range enhancer element whose ability to bind the transcription factors POU2F and YY1 is significantly impaired, with decreased activity in the presence of the C- allele compared with the wild type A-allele. An enhancer RNA (eRNA) is transcribed in thyroid tissue from this region and is greatly downregulated in NMTC tumors. We suggest that this is an example of an ultra-rare mutation predisposing to thyroid cancer with high penetrance.

A germline mutation in SRRM2, a splicing factor gene, is implicated in papillary thyroid carcinoma predisposition

Papillary thyroid carcinoma (PTC) displays strong but so far largely uncharacterized heritability. Here we studied genetic predisposition in a family with six affected individuals. We genotyped all available family members and conducted whole exome sequencing of blood DNA from two affected individuals. Haplotype analysis and other genetic criteria narrowed our list of candidates to a germline variant in the serine/arginine repetitive matrix 2 gene (SRRM2). This heterozygous variant, c.1037C > T (Ser346Phe or S346F; rs149019598) cosegregated with PTC in the family. It was not found in 138 other PTC families. It was found in 7/1,170 sporadic PTC cases and in 0/1,404 controls (p = 0.004). The encoded protein SRRM2 (also called SRm300) is part of the RNA splicing machinery. To evaluate the possibility that the S346F missense mutation affects alternative splicing, we compared RNA-Seq data in leukocytes from three mutation carriers and three controls. Significant differences in alternative splicing were identified for 1,642 exons, of which a subset of 7 exons was verified experimentally. The results confirmed a higher ratio of inclusion of exons in mutation carriers. These data suggest that the S346F mutation in SRRM2 predisposes to PTC by affecting alternative splicing of unidentified downstream target genes.