NusG controls transcription pausing and RNA polymerase translocation throughout the Bacillus subtilis genome.

Transcription is punctuated by RNA polymerase (RNAP) pausing. These pauses provide time for diverse regulatory events that can modulate gene expression. Transcription elongation factors dramatically affect RNAP pausing in vitro, but the genome-wide role of such factors on pausing has not been examined. Using native elongating transcript sequencing followed by RNase digestion (RNET-seq), we analyzed RNAP pausing in Bacillus subtilis genome-wide and identified an extensive role of NusG in pausing. This universally conserved transcription elongation factor is known as Spt5 in archaeal and eukaryotic organisms. B. subtilis NusG shifts RNAP to the posttranslocation register and induces pausing at 1,600 sites containing a consensus TTNTTT motif in the nontemplate DNA strand within the paused transcription bubble. The TTNTTT motif is necessary but not sufficient for NusG-dependent pausing. Approximately one-fourth of these pause sites were localized to untranslated regions and could participate in posttranscription initiation control of gene expression as was previously shown for tlrB and the trpEDCFBA operon. Most of the remaining pause sites were identified in protein-coding sequences. NusG-dependent pausing was confirmed for all 10 pause sites that we tested in vitro. Putative pause hairpins were identified for 225 of the 342 strongest NusG-dependent pause sites, and some of these hairpins were shown to function in vitro. NusG-dependent pausing in the ribD riboswitch provides time for cotranscriptional binding of flavin mononucleotide, which decreases the concentration required for termination upstream of the ribD coding sequence. Our phylogenetic analysis implicates NusG-dependent pausing as a widespread mechanism in bacteria.

Genome-wide DNA methylation patterns in LSH mutant reveals de-repression of repeat elements and redundant epigenetic silencing pathways.

Cytosine methylation is critical in mammalian development and plays a role in diverse biologic processes such as genomic imprinting, X chromosome inactivation, and silencing of repeat elements. Several factors regulate DNA methylation in early embryogenesis, but their precise role in the establishment of DNA methylation at a given site remains unclear. We have generated a comprehensive methylation map in fibroblasts derived from the murine DNA methylation mutant Hells(-/-) (helicase, lymphoid specific, also known as LSH). It has been previously shown that HELLS can influence de novo methylation of retroviral sequences and endogenous genes. Here, we describe that HELLS controls cytosine methylation in a nuclear compartment that is in part defined by lamin B1 attachment regions. Despite widespread loss of cytosine methylation at regulatory sequences, including promoter regions of protein-coding genes and noncoding RNA genes, overall relative transcript abundance levels in the absence of HELLS are similar to those in wild-type cells. A subset of promoter regions shows increases of the histone modification H3K27me3, suggesting redundancy of epigenetic silencing mechanisms. Furthermore, HELLS modulates CG methylation at all classes of repeat elements and is critical for repression of a subset of repeat elements. Overall, we provide a detailed analysis of gene expression changes in relation to DNA methylation alterations, which contributes to our understanding of the biological role of cytosine methylation.

CG hypomethylation in Lsh-/- mouse embryonic fibroblasts is associated with de novo H3K4me1 formation and altered cellular plasticity.

DNA methylation patterns are established in early embryogenesis and are critical for cellular differentiation. To investigate the role of CG methylation in potential enhancer formation, we assessed H3K4me1 modification in murine embryonic fibroblasts (MEFs) derived from the DNA methylation mutant Lsh(-/-) mice. We report here de novo formation of putative enhancer elements at CG hypomethylated sites that can be dynamically altered. We found a subset of differentially enriched H3K4me1 regions clustered at neuronal lineage genes and overlapping with known cis-regulatory elements present in brain tissue. Reprogramming of Lsh(-/-) MEFs into induced pluripotent stem (iPS) cells leads to increased neuronal lineage gene expression of premarked genes and enhanced differentiation potential of Lsh(-/-) iPS cells toward the neuronal lineage pathway compared with WT iPS cells in vitro and in vivo. The state of CG hypomethylation and H3K4me1 enrichment is partially maintained in Lsh(-/-) iPS cells. The acquisition of H3K27ac and activity of subcloned fragments in an enhancer reporter assay indicate functional activity of several of de novo H3K4me1-marked sequences. Our results suggest a functional link of H3K4me1 enrichment at CG hypomethylated sites, enhancer formation, and cellular plasticity.

The principal genetic determinants for nasopharyngeal carcinoma in China involve the HLA class I antigen recognition groove.

Nasopharyngeal carcinoma (NPC) is an epithelial malignancy facilitated by Epstein-Barr Virus infection. Here we resolve the major genetic influences for NPC incidence using a genome-wide association study (GWAS), independent cohort replication, and high-resolution molecular HLA class I gene typing including 4,055 study participants from the Guangxi Zhuang Autonomous Region and Guangdong province of southern China. We detect and replicate strong association signals involving SNPs, HLA alleles, and amino acid (aa) variants across the major histocompatibility complex-HLA-A, HLA -B, and HLA -C class I genes (P(HLA-A-aa-site-62) = 7.4 × 10(-29); P (HLA-B-aa-site-116) = 6.5 × 10(-19); P (HLA-C-aa-site-156) = 6.8 × 10(-8) respectively). Over 250 NPC-HLA associated variants within HLA were analyzed in concert to resolve separate and largely independent HLA-A, -B, and -C gene influences. Multivariate logistical regression analysis collapsed significant associations in adjacent genes spanning 500 kb (OR2H1, GABBR1, HLA-F, and HCG9) as proxies for peptide binding motifs carried by HLA- A*11:01. A similar analysis resolved an independent association signal driven by HLA-B*13:01, B*38:02, and B*55:02 alleles together. NPC resistance alleles carrying the strongly associated amino acid variants implicate specific class I peptide recognition motifs in HLA-A and -B peptide binding groove as conferring strong genetic influence on the development of NPC in China.

Genome-wide association study implicates PARD3B-based AIDS restriction.

BACKGROUND: Host genetic variation influences human immunodeficiency virus (HIV) infection and progression to AIDS. Here we used clinically well-characterized subjects from 5 pretreatment HIV/AIDS cohorts for a genome-wide association study to identify gene associations with rate of AIDS progression. METHODS: European American HIV seroconverters (n = 755) were interrogated for single-nucleotide polymorphisms (SNPs) (n = 700,022) associated with progression to AIDS 1987 (Cox proportional hazards regression analysis, co-dominant model). RESULTS: Association with slower progression was observed for SNPs in the gene PARD3B. One of these, rs11884476, reached genome-wide significance (relative hazard = 0.3; P =3. 370 × 10(-9)) after statistical correction for 700,022 SNPs and contributes 4.52% of the overall variance in AIDS progression in this study. Nine of the top-ranked SNPs define a PARD3B haplotype that also displays significant association with progression to AIDS (hazard ratio, 0.3; P = 3.220 × 10(-8)). One of these SNPs, rs10185378, is a predicted exonic splicing enhancer; significant alteration in the expression profile of PARD3B splicing transcripts was observed in B cell lines with alternate rs10185378 genotypes. This SNP was typed in European cohorts of rapid progressors and was found to be protective for AIDS 1993 definition (odds ratio, 0.43, P = .025). CONCLUSIONS: These observations suggest a potential unsuspected pathway of host genetic influence on the dynamics of AIDS progression.

Genetic associations of variants in genes encoding HIV-dependency factors required for HIV-1 infection.

BACKGROUND: High-throughput genome-wide techniques have facilitated the identification of previously unknown host proteins involved in cellular human immunodeficiency virus (HIV) infection. Recently, 3 independent studies have used small interfering RNA technology to silence each gene in the human genome to determine the importance of each in HIV infection. Genes conferring a significant effect were termed HIV-dependency factors (HDFs). METHODS: We assembled high-density panels of 6380 single-nucleotide polymorphisms (SNPs) in 278 HDF genes and tested for genotype associations with HIV infection and AIDS progression in 1633 individuals from clinical AIDS cohorts. RESULTS: After statistical correction for multiple tests, significant associations with HIV acquisition were found for SNPs in 2 genes, NCOR2 and IDH1. Weaker associations with AIDS progression were revealed for SNPs within the TM9SF2 and EGFR genes. CONCLUSIONS: This study independently verifies the influence of NCOR2 and IDH1 on HIV transmission, and its findings suggest that variation in these genes affects susceptibility to HIV infection in exposed individuals.

Examination of disease-based selection, demographic history and population structure in European Y-chromosome haplogroup I.

We attempted to refine the understanding of an association of Y-chromosomal haplogroup I (hg-I) with enhanced AIDS progression that had been previously reported. First, we compared the progression phenotype between hg-I and its phylogenetically closest haplogroup J. Then, we took a candidate gene approach resequencing DDX3Y, a crucial autoimmunity gene, in hg-I and other common European Y-chromosome haplogroups looking for functional variants. We extended the genetic analyses to CD24L4 and compared and contrasted the roles of disease-based selection, demographic history and population structure shaping the contemporary genetic landscape of hg-I chromosomes. Our results confirmed and refined the AIDS progression signal to hg-I, though no gene variant was identified that can explain the disease association. Molecular evolutionary and genetic analyses of the examined loci suggested a unique evolutionary history in hg-I, probably shaped by complex interactions of selection, demographic history and high geographical differentiation leading to the formation of distinct hg-I subhaplogroups that today are associated with HIV/AIDS onset. Clearly, further studies on Y-chromosome candidate loci sequencing to discover functional variants and discern the roles of evolutionary factors are warranted.

Evolution of microsatellite loci in the adaptive radiation of Hawaiian honeycreepers.

Previous studies have examined germ-line mutations to infer the processes that generate and maintain variability in microsatellite loci. Few studies, however, have examined patterns to infer processes that act on microsatellite loci over evolutionary time. Here, we examine changes in 8 dinucleotide loci across the adaptive radiation of Hawaiian honeycreepers. The loci were found to be highly variable across the radiation, and we did not detect ascertainment bias with respect to allelic diversity or allele size ranges. In examining patterns at the sequence level, we found that changes in flanking regions, repeat motifs, or repeat interruptions were often shared between closely related species and may be phylogenetically informative. Genetic distance measures based on microsatellites were strongly correlated with those based on mitochondrial DNA (mtDNA) sequences as well as with divergence time up to 3 My. Phylogenetic inferences based on microsatellite genetic distances consistently recovered 2 of the 4 honeycreeper clades observed in a tree based on mtDNA sequences but differed from the mtDNA tree in the relationships among clades. Our results confirm that microsatellite loci may be conserved over evolutionary time, making them useful in population-level studies of species that diverged from the species in which they were characterized as long as 5 Ma. Despite this, we found that their use in phylogenetic inference was limited to closely related honeycreeper species.

FIV cross-species transmission: an evolutionary prospective.

Feline and primate immunodeficiency viruses (FIVs, SIVs, and HIV) are transmitted via direct contact (e.g. fighting, sexual contact, and mother-offspring transmission). This dynamic likely poses a behavioral barrier to cross-species transmission in the wild. Recently, several host intracellular anti-viral proteins that contribute to species-specificity of primate lentiviruses have been identified revealing adaptive mechanisms that further limit spread of lentiviruses between species. Consistent with these inter-species transmission barriers, phylogenetic evidence supports the prediction that FIV transmission is an exceedingly rare event between free-ranging cat species, though it has occurred occasionally in captive settings. Recently we documented that puma and bobcats in Southern California share an FIV strain, providing an opportunity to evaluate evolution of both viral strains and host intracellular restriction proteins. These studies are facilitated by the availability of the 2x cat genome sequence annotation. In addition, concurrent viral and host genetic analyses have been used to track patterns of migration of the host species and barriers to transmission of the virus within the African lion. These studies illustrate the utility of FIV as a model to discover the variables necessary for establishment and control of lentiviral infections in new species.

Metagenomic Discovery of 83 New Human Papillomavirus Types in Patients with Immunodeficiency.

Several immunodeficiencies are associated with high susceptibility to persistent and progressive human papillomavirus (HPV) infection leading to a wide range of cutaneous and mucosal lesions. However, the HPV types most commonly associated with such clinical manifestations in these patients have not been systematically defined. Here, we used virion enrichment, rolling circle amplification, and deep sequencing to identify circular DNA viruses present in skin swabs and/or wart biopsy samples from 48 patients with rare genetic immunodeficiencies, including patients with warts, hypogammaglobulinemia, infections, myelokathexis (WHIM) syndrome, or epidermodysplasia verruciformis (EV). Their profiles were compared with the profiles of swabs from 14 healthy adults and warts from 6 immunologically normal children. Individual patients were typically infected with multiple HPV types; up to 26 different types were isolated from a single patient (multiple anatomical sites, one time point). Among these, we identified the complete genomes of 83 previously unknown HPV types and 35 incomplete genomes representing possible additional new types. HPV types in the genus Gammapapillomavirus were common in WHIM patients, whereas EV patients mainly shed HPVs from the genus Betapapillomavirus. Preliminary evidence based on three WHIM patients treated with plerixafor, a leukocyte mobilizing agent, suggest that longer-term therapy may correlate with decreased HPV diversity and increased predominance of HPV types associated with childhood skin warts.IMPORTANCE Although some members of the viral family Papillomaviridae cause benign skin warts (papillomas), many human papillomavirus (HPV) infections are not associated with visible symptoms. For example, most healthy adults chronically shed Gammapapillomavirus (Gamma) virions from apparently healthy skin surfaces. To further explore the diversity of papillomaviruses, we performed viromic surveys on immunodeficient individuals suffering from florid skin warts. Our results nearly double the number of known Gamma HPV types and suggest that WHIM syndrome patients are uniquely susceptible to Gamma HPV-associated skin warts. Preliminary results suggest that treatment with the drug plerixafor may promote resolution of the unusual Gamma HPV skin warts observed in WHIM patients.

Visualizing translocation dynamics and nascent transcript errors in paused RNA polymerases in vivo.

BACKGROUND: Transcription elongation is frequently interrupted by pausing signals in DNA, with downstream effects on gene expression. Transcription errors also induce prolonged pausing, which can lead to a destabilized genome by interfering with DNA replication. Mechanisms of pausing associated with translocation blocks and misincorporation have been characterized in vitro, but not in vivo. RESULTS: We investigate the pausing pattern of RNA polymerase (RNAP) in Escherichia coli by a novel approach, combining native elongating transcript sequencing (NET-seq) with RNase footprinting of the transcripts (RNET-seq). We reveal that the G-dC base pair at the 5' end of the RNA-DNA hybrid interferes with RNAP translocation. The distance between the 5' G-dC base pair and the 3' end of RNA fluctuates over a three-nucleotide width. Thus, the G-dC base pair can induce pausing in post-translocated, pre-translocated, and backtracked states of RNAP. Additionally, a CpG sequence of the template DNA strand spanning the active site of RNAP inhibits elongation and induces G-to-A errors, which leads to backtracking of RNAP. Gre factors efficiently proofread the errors and rescue the backtracked complexes. We also find that pausing events are enriched in the 5' untranslated region and antisense transcription of mRNA genes and are reduced in rRNA genes. CONCLUSIONS: In E. coli, robust transcriptional pausing involves RNAP interaction with G-dC at the upstream end of the RNA-DNA hybrid, which interferes with translocation. CpG DNA sequences induce transcriptional pausing and G-to-A errors.

As the raven flies: using genetic data to infer the history of invasive common raven (Corvus corax) populations in the Mojave Desert.

Common raven (Corvus corax) populations in Mojave Desert regions of southern California and Nevada have increased dramatically over the past five decades. This growth has been attributed to increased human development in the region, as ravens have a commensal relationship with humans and feed extensively at landfills and on road-killed wildlife. Ravens, as a partially subsidized predator, also represent a problem for native desert wildlife, in particular threatened desert tortoises (Gopherus agassizii). However, it is unclear whether the more than 15-fold population increase is due to in situ population growth or to immigration from adjacent regions where ravens have been historically common. Ravens were sampled for genetic analysis at several local sites within five major areas: the West Mojave Desert (California), East Mojave Desert (southern Nevada), southern coastal California, northern coastal California (Bay Area), and northern Nevada (Great Basin). Analyses of mtDNA control region sequences reveal an increased frequency of raven 'Holarctic clade' haplotypes from south to north inland, with 'California clade' haplotypes nearly fixed in the California populations. There was significant structuring among regions for mtDNA, with high F(ST) values among sampling regions, especially between the Nevada and California samples. Analyses of eight microsatellite loci reveal a mostly similar pattern of regional population structure, with considerably smaller, but mostly significant, values. The greater mtDNA divergences may be due to lower female dispersal relative to males, lower N(e), or effects of high mutation rates on maximal values of F(ST). Analyses indicate recent population growth in the West Mojave Desert and a bottleneck in the northern California populations. While we cannot rule out in situ population growth as a factor, patterns of movement inferred from our data suggest that the increase in raven populations in the West Mojave Desert resulted from movements from southern California and the Central Valley. Ravens in the East Mojave Desert are more similar to ones from northern Nevada, indicating movement between those regions. If this interpretation of high gene flow into the Mojave Desert is correct, then efforts to manage raven numbers by local control may not be optimally effective.

Natural selection of the major histocompatibility complex (Mhc) in Hawaiian honeycreepers (Drepanidinae).

The native Hawaiian honeycreepers represent a classic example of adaptive radiation and speciation, but currently face one the highest extinction rates in the world. Although multiple factors have likely influenced the fate of Hawaiian birds, the relatively recent introduction of avian malaria is thought to be a major factor limiting honeycreeper distribution and abundance. We have initiated genetic analyses of class II beta chain Mhc genes in four species of honeycreepers using methods that eliminate the possibility of sequencing mosaic variants formed by cloning heteroduplexed polymerase chain reaction products. Phylogenetic analyses group the honeycreeper Mhc sequences into two distinct clusters. Variation within one cluster is high, with dN > dS and levels of diversity similar to other studies of Mhc (B system) genes in birds. The second cluster is nearly invariant and includes sequences from honeycreepers (Fringillidae), a sparrow (Emberizidae) and a blackbird (Emberizidae). This highly conserved cluster appears reminiscent of the independently segregating Rfp-Y system of genes defined in chickens. The notion that balancing selection operates at the Mhc in the honeycreepers is supported by transpecies polymorphism and strikingly high dN/dS ratios at codons putatively involved in peptide interaction. Mitochondrial DNA control region sequences were invariant in the i'iwi, but were highly variable in the 'amakihi. By contrast, levels of variability of class II beta chain Mhc sequence codons that are hypothesized to be directly involved in peptide interactions appear comparable between i'iwi and 'amakihi. In the i'iwi, natural selection may have maintained variation within the Mhc, even in the face of what appears to a genetic bottleneck.