The selfish yeast plasmid exploits a SWI/SNF-type chromatin remodeling complex for hitchhiking on chromosomes and ensuring high-fidelity propagation.

Extra-chromosomal selfish DNA elements can evade the risk of being lost at every generation by behaving as chromosome appendages, thereby ensuring high fidelity segregation and stable persistence in host cell populations. The yeast 2-micron plasmid and episomes of the mammalian gammaherpes and papilloma viruses that tether to chromosomes and segregate by hitchhiking on them exemplify this strategy. We document for the first time the utilization of a SWI/SNF-type chromatin remodeling complex as a conduit for chromosome association by a selfish element. One principal mechanism for chromosome tethering by the 2-micron plasmid is the bridging interaction of the plasmid partitioning proteins (Rep1 and Rep2) with the yeast RSC2 complex and the plasmid partitioning locus STB. We substantiate this model by multiple lines of evidence derived from genomics, cell biology and interaction analyses. We describe a Rep-STB bypass system in which a plasmid engineered to non-covalently associate with the RSC complex mimics segregation by chromosome hitchhiking. Given the ubiquitous prevalence of SWI/SNF family chromatin remodeling complexes among eukaryotes, it is likely that the 2-micron plasmid paradigm or analogous ones will be encountered among other eukaryotic selfish elements.

Disruption of cellular proteostasis by H1N1 influenza A virus causes α-synuclein aggregation.

Neurodegenerative diseases feature specific misfolded or misassembled proteins associated with neurotoxicity. The precise mechanisms by which protein aggregates first arise in the majority of sporadic cases have remained unclear. Likely, a first critical mass of misfolded proteins starts a vicious cycle of a prion-like expansion. We hypothesize that viruses, having evolved to hijack the host cellular machinery for catalyzing their replication, lead to profound disturbances of cellular proteostasis, resulting in such a critical mass of protein aggregates. Here, we investigated the effect of influenza virus (H1N1) strains on proteostasis of proteins associated with neurodegenerative diseases in Lund human mesencephalic dopaminergic cells in vitro and infection of Rag knockout mice in vivo. We demonstrate that acute H1N1 infection leads to the formation of α-synuclein and Disrupted-in-Schizophrenia 1 (DISC1) aggregates, but not of tau or TDP-43 aggregates, indicating a selective effect on proteostasis. Oseltamivir phosphate, an antiinfluenza drug, prevented H1N1-induced α-synuclein aggregation. As a cell pathobiological mechanism, we identified H1N1-induced blocking of autophagosome formation and inhibition of autophagic flux. In addition, α-synuclein aggregates appeared in infected cell populations connected to the olfactory bulbs following intranasal instillation of H1N1 in Rag knockout mice. We propose that H1N1 virus replication in neuronal cells can induce seeds of aggregated α-synuclein or DISC1 that may be able to initiate further detrimental downstream events and should thus be considered a risk factor in the pathogenesis of synucleinopathies or a subset of mental disorders. More generally, aberrant proteostasis induced by viruses may be an underappreciated factor in initiating protein misfolding.

MicroRNAs reinforce repression of PRC2 transcriptional targets independently and through a feed-forward regulatory network.

Gene expression can be regulated at multiple levels, but it is not known if and how there is broad coordination between regulation at the transcriptional and post-transcriptional levels. Transcription factors and chromatin regulate gene expression transcriptionally, whereas microRNAs (miRNAs) are small regulatory RNAs that function post-transcriptionally. Systematically identifying the post-transcriptional targets of miRNAs and the mechanism of transcriptional regulation of the same targets can shed light on regulatory networks connecting transcriptional and post-transcriptional control. We used individual-nucleotide resolution UV crosslinking and immunoprecipitation (iCLIP) for the RNA-induced silencing complex (RISC) component AGO2 and global miRNA depletion to identify genes directly targeted by miRNAs. We found that Polycomb repressive complex 2 (PRC2) and its associated histone mark, H3K27me3, is enriched at hundreds of miRNA-repressed genes. We show that these genes are directly repressed by PRC2 and constitute a significant proportion of direct PRC2 targets. For just over half of the genes corepressed by PRC2 and miRNAs, PRC2 promotes their miRNA-mediated repression by increasing expression of the miRNAs that are likely to target them. miRNAs also repress the remainder of the PRC2 target genes, but independently of PRC2. Thus, miRNAs post-transcriptionally reinforce silencing of PRC2-repressed genes that are inefficiently repressed at the level of chromatin, by either forming a feed-forward regulatory network with PRC2 or repressing them independently of PRC2.

Identification of an Antiretroviral Small Molecule That Appears To Be a Host-Targeting Inhibitor of HIV-1 Assembly.

Given the projected increase in multidrug-resistant HIV-1, there is an urgent need for development of antiretrovirals that act on virus life cycle stages not targeted by drugs currently in use. Host-targeting compounds are of particular interest because they can offer a high barrier to resistance. Here, we report identification of two related small molecules that inhibit HIV-1 late events, a part of the HIV-1 life cycle for which potent and specific inhibitors are lacking. This chemotype was discovered using cell-free protein synthesis and assembly systems that recapitulate intracellular host-catalyzed viral capsid assembly pathways. These compounds inhibit replication of HIV-1 in human T cell lines and peripheral blood mononuclear cells, and are effective against a primary isolate. They reduce virus production, likely by inhibiting a posttranslational step in HIV-1 Gag assembly. Notably, the compound colocalizes with HIV-1 Gag in situ; however, unexpectedly, selection experiments failed to identify compound-specific resistance mutations in gag or pol, even though known resistance mutations developed upon parallel nelfinavir selection. Thus, we hypothesized that instead of binding to Gag directly, these compounds localize to assembly intermediates, the intracellular multiprotein complexes containing Gag and host factors that form during immature HIV-1 capsid assembly. Indeed, imaging of infected cells shows compound colocalized with two host enzymes found in assembly intermediates, ABCE1 and DDX6, but not two host proteins found in other complexes. While the exact target and mechanism of action of this chemotype remain to be determined, our findings suggest that these compounds represent first-in-class, host-targeting inhibitors of intracellular events in HIV-1 assembly.IMPORTANCE The success of antiretroviral treatment for HIV-1 is at risk of being undermined by the growing problem of drug resistance. Thus, there is a need to identify antiretrovirals that act on viral life cycle stages not targeted by drugs in use, such as the events of HIV-1 Gag assembly. To address this gap, we developed a compound screen that recapitulates the intracellular events of HIV-1 assembly, including virus-host interactions that promote assembly. This effort led to the identification of a new chemotype that inhibits HIV-1 replication at nanomolar concentrations, likely by acting on assembly. This compound colocalized with Gag and two host enzymes that facilitate capsid assembly. However, resistance selection did not result in compound-specific mutations in gag, suggesting that the chemotype does not directly target Gag. We hypothesize that this chemotype represents a first-in-class inhibitor of virus production that acts by targeting a virus-host complex important for HIV-1 Gag assembly.

The histone variant H2A.Z in yeast is almost exclusively incorporated into the +1 nucleosome in the direction of transcription.

Transcription start sites (TSS) in eukaryotes are characterized by a nucleosome-depleted region (NDR), which appears to be flanked upstream and downstream by strongly positioned nucleosomes incorporating the histone variant H2A.Z. H2A.Z associates with both active and repressed TSS and is important for priming genes for rapid transcriptional activation. However, the determinants of H2A.Z occupancy at specific nucleosomes and its relationship to transcription initiation remain unclear. To further elucidate the specificity of H2A.Z, we determined its genomic localization at single nucleosome resolution, as well as the localization of its chromatin remodelers Swr1 and Ino80. By analyzing H2A.Z occupancy in conjunction with RNA expression data that captures promoter-derived antisense initiation, we find that H2A.Z's bimodal incorporation on either side of the NDR is not a general feature of TSS, but is specifically a marker for bidirectional transcription, such that the upstream flanking -1 H2A.Z-containing nucleosome is more appropriately considered as a +1 H2A.Z nucleosome for antisense transcription. The localization of H2A.Z almost exclusively at the +1 nucleosome suggests that a transcription-initiation dependent process could contribute to its specific incorporation.

The specificity of H2A.Z occupancy in the yeast genome and its relationship to transcription.

The incorporation of histone variants into nucleosomes has important functional consequences in all aspects of eukaryotic chromatin biology. H2A.Z is a conserved histone variant found in all eukaryotes from yeast to mammals. Recent studies in yeast have shed light on the questions of where and how nucleosomes containing this variant are situated at promoters and in relation to genes, and what its specificity implies with regard to transcription. In yeast, H2A.Z appears to be primarily incorporated into the first nucleosome in the direction of transcription initiation, either of an mRNA transcript or a divergently transcribed upstream antisense non-coding RNA. This specificity of H2A.Z is due in part to the localization at promoters of SWR1, the ATP-dependent chromatin remodeler that incorporates H2A.Z into nucleosomes. Replacement of H2A.Z with canonical H2A is dependent on the function of the transcription pre-initiation complex. The recent studies summarized in this review reveal that the directionality of H2A.Z occupancy in relation to transcription thus reflects a balance of incorporation and eviction activities, which likely have varying contributions at distinct sets of genes across the genome.

Identification and removal of sequencing artifacts produced by mispriming during reverse transcription in multiple RNA-seq technologies.

The quality of RNA sequencing data relies on specific priming by the primer used for reverse transcription (RT-primer). Nonspecific annealing of the RT-primer to the RNA template can generate reads with incorrect cDNA ends and can cause misinterpretation of data (RT mispriming). This kind of artifact in RNA-seq based technologies is underappreciated and currently no adequate tools exist to computationally remove them from published data sets. We show that mispriming can occur with as little as two bases of complementarity at the 3' end of the primer followed by intermittent regions of complementarity. We also provide a computational pipeline that identifies cDNA reads produced from RT mispriming, allowing users to filter them out from any aligned data set. Using this analysis pipeline, we identify thousands of mispriming events in a dozen published data sets from diverse technologies including short RNA-seq, total/mRNA-seq, HITS-CLIP, and GRO-seq. We further show how RT mispriming can lead to misinterpretation of data. In addition to providing a solution to computationally remove RT-misprimed reads, we also propose an experimental solution to completely avoid RT-mispriming by performing RNA-seq using thermostable group II intron derived reverse transcriptase (TGIRT-seq).

The Determinants of Directionality in Transcriptional Initiation.

A new paradigm has emerged in recent years characterizing transcription initiation as a bidirectional process encompassing a larger proportion of the genome than previously thought. Past concepts of coding genes thinly scattered among a vast background of transcriptionally inert noncoding DNA have been abandoned. A richer picture has taken shape, integrating transcription of coding genes, enhancer RNAs (eRNAs), and various other noncoding transcriptional events. In this review we give an overview of recent studies detailing the mechanisms of RNA polymerase II (RNA Pol II)-based transcriptional initiation and discuss the ways in which transcriptional direction is established as well as its functional implications.

Formation of RNA Granule-Derived Capsid Assembly Intermediates Appears To Be Conserved between Human Immunodeficiency Virus Type 1 and the Nonprimate Lentivirus Feline Immunodeficiency Virus.

During immature capsid assembly in cells, human immunodeficiency virus type 1 (HIV-1) Gag co-opts a host RNA granule, forming a pathway of intracellular assembly intermediates containing host components, including two cellular facilitators of assembly, ABCE1 and DDX6. A similar assembly pathway has been observed for other primate lentiviruses. Here we asked whether feline immunodeficiency virus (FIV), a nonprimate lentivirus, also forms RNA granule-derived capsid assembly intermediates. First, we showed that the released FIV immature capsid and a large FIV Gag-containing intracellular complex are unstable during analysis, unlike for HIV-1. We identified harvest conditions, including in situ cross-linking, that overcame this problem, revealing a series of FIV Gag-containing complexes corresponding in size to HIV-1 assembly intermediates. Previously, we showed that assembly-defective HIV-1 Gag mutants are arrested at specific assembly intermediates; here we identified four assembly-defective FIV Gag mutants, including three not previously studied, and demonstrated that they appear to be arrested at the same intermediate as the cognate HIV-1 mutants. Further evidence that these FIV Gag-containing complexes correspond to assembly intermediates came from coimmunoprecipitations demonstrating that endogenous ABCE1 and the RNA granule protein DDX6 are associated with FIV Gag, as shown previously for HIV-1 Gag, but are not associated with a ribosomal protein, at steady state. Additionally, we showed that FIV Gag associates with another RNA granule protein, DCP2. Finally, we validated the FIV Gag-ABCE1 and FIV Gag-DCP2 interactions with proximity ligation assays demonstrating colocalization in situ Together, these data support a model in which primate and nonprimate lentiviruses form intracellular capsid assembly intermediates derived from nontranslating host RNA granules.IMPORTANCE Like HIV-1 Gag, FIV Gag assembles into immature capsids; however, it is not known whether FIV Gag progresses through a pathway of immature capsid assembly intermediates derived from host RNA granules, as shown for HIV-1 Gag. Here we showed that FIV Gag forms complexes that resemble HIV-1 capsid assembly intermediates in size and in their association with ABCE1 and DDX6, two host facilitators of HIV-1 immature capsid assembly that are found in HIV-1 assembly intermediates. Our studies also showed that known and novel assembly-defective FIV Gag mutants fail to progress past putative intermediates in a pattern resembling that observed for HIV-1 Gag mutants. Finally, we used imaging to demonstrate colocalization of FIV Gag with ABCE1 and with the RNA granule protein DCP2. Thus, we conclude that formation of assembly intermediates derived from host RNA granules is likely conserved between primate and nonprimate lentiviruses and could provide targets for future antiviral strategies.

The ATP-dependent chromatin remodeler Chd1 is recruited by transcription elongation factors and maintains H3K4me3/H3K36me3 domains at actively transcribed and spliced genes.

Chd1 (Chromodomain Helicase DNA Binding Protein 1) is a conserved ATP-dependent chromatin remodeler that maintains the nucleosomal structure of chromatin, but the determinants of its specificity and its impact on gene expression are not well defined. To identify the determinants of Chd1 binding specificity in the yeast genome, we investigated Chd1 occupancy in mutants of several candidate factors. We found that several components of the PAF1 transcription elongation complex contribute to Chd1 recruitment to highly transcribed genes and identified Spt4 as a factor that appears to negatively modulate Chd1 binding to chromatin. We discovered that CHD1 loss alters H3K4me3 and H3K36me3 patterns throughout the yeast genome. Interestingly, the aberrant histone H3 methylation patterns were predominantly observed within 1 kb from the transcription start site, where both histone H3 methylation marks co-occur. A reciprocal change between the two marks was obvious in the absence of Chd1, suggesting a role for CHD1 in establishing or maintaining the boundaries of these largely mutually exclusive histone marks. Strikingly, intron-containing genes were most susceptible to CHD1 loss and exhibited a high degree of histone H3 methylation changes. Intron retention was significantly lower in the absence of CHD1, suggesting that CHD1 function as a chromatin remodeler could indirectly affect RNA splicing.

Subtype-specific addiction of the activated B-cell subset of diffuse large B-cell lymphoma to FOXP1.

High expression of the forkhead box P1 (FOXP1) transcription factor distinguishes the aggressive activated B cell (ABC) diffuse large B-cell lymphoma (DLBCL) subtype from the better prognosis germinal center B-cell (GCB)-DLBCL subtype and is highly correlated with poor outcomes. A genetic or functional role for FOXP1 in lymphomagenesis, however, remains unknown. Here, we report that sustained FOXP1 expression is vital for ABC-DLBCL cell-line survival. Genome-wide analyses revealed direct and indirect FOXP1 transcriptional enforcement of ABC-DLBCL hallmarks, including the classical NF-κB and MYD88 (myeloid differentiation primary response gene 88) pathways. FOXP1 promoted gene expression underlying transition of the GCB cell to the plasmablast--the transient B-cell stage targeted in ABC-DLBCL transformation--by antagonizing pathways distinctive of GCB-DLBCL, including that of the GCB "master regulator," BCL6 (B-cell lymphoma 6). Cell-line derived FOXP1 target genes that were highly correlated with FOXP1 expression in primary DLBCL accurately segregated the corresponding clinical subtypes of a large cohort of primary DLBCL isolates and identified conserved pathways associated with ABC-DLBCL pathology.

Correction: Quantitative genetics of CTCF binding reveal local sequence effects and different modes of X-chromosome association.

[This corrects the article DOI: 10.1371/journal.pgen.1004798.].

Ethylene induces combinatorial effects of histone H3 acetylation in gene expression in Arabidopsis.

BACKGROUND: Histone acetylation and deacetylation are essential for gene regulation and have been implicated in the regulation of plant hormone responses. Many studies have indicated the role of histone acetylation in ethylene signaling; however, few studies have investigated how ethylene signaling regulates the genomic landscape of chromatin states. Recently, we found that ethylene can specifically elevate histone H3K14 acetylation and the non-canonical histone H3K23 acetylation in etiolated seedlings and the gene activation is positively associated with the elevation of H3K14Ac and H3K23Ac in response to ethylene. To assess the role of H3K9, H3K14, and H3K23 histone modifications in the ethylene response, we examined how ethylene regulates histone acetylation and the transcriptome at global level and in ethylene regulated genes both in wild type (Col-0) and ein2-5 seedlings. RESULTS: Our results revealed that H3K9Ac, H3K14Ac, and H3K23Ac are preferentially enriched around the transcription start sites and are positively correlated with gene expression levels in Col-0 and ein2-5 seedlings both with and without ethylene treatment. In the absence of ethylene, no combinatorial effect of H3K9Ac, H3K14Ac, and H3K23Ac on gene expression was detected. In the presence of ethylene, however, combined enrichment of the three histone acetylation marks was associated with high gene expression levels, and this ethylene-induced change was EIN2 dependent. In addition, we found that ethylene-regulated genes are expressed at medium or high levels, and a group of ethylene regulated genes are marked by either one of H3K9Ac, H3K14Ac or H3K23Ac. In this group of genes, the levels of H3K9Ac were altered by ethylene, but in the absence of ethylene the levels of H3K9Ac and peak breadths are distinguished in up- and down- regulated genes. In the presence of ethylene, the changes in the peak breadths and levels of H3K14Ac and H3K23Ac are required for the alteration of gene expressions. CONCLUSIONS: Our study reveals that the plant hormone ethylene induces combinatorial effects of H3K9Ac, K14Ac and K23Ac histone acetylation in gene expression genome widely. Further, for a group of ethylene regulated genes, in the absence of ethylene the levels and the covered breadths of H3K9Ac are the preexist markers for distinguishing up- and down- regulated genes, the change in the peak breadths and levels of H3K14Ac and H3K23Ac are required for the alteration of gene expression in the presence of ethylene.

Divergent effects of Porcupine and Wntless on WNT1 trafficking, secretion, and signaling.

Loss-of-function studies have identified Porcupine (PORCN) and Wntless (WLS) as essential mediators of Wnt secretion and signaling. Whereas PORCN is thought to palmitoylate Wnt proteins, WLS is believed to transport palmitoylated Wnt proteins to the cell surface. However, little is known about how these two proteins cooperate to regulate Wnt palmitoylation, trafficking, secretion, and signaling. We first investigated possible interactions between PORCN, WLS, and WNT1, by carrying out co-immunoprecipitation studies. These studies demonstrate the existence of a complex containing PORCN and WLS. They further show that PORCN and WLS compete for binding to WNT1. Then, we used gain-of-function studies to investigate the cooperation between PORCN and WLS as well as possible biochemical interactions between PORCN, WLS, and WNT1. Consistent with the proposed roles for PORCN and WLS, we show that overexpression of PORCN promotes palmitoylation of WNT1 while overexpression of WLS does not. Overexpression of PORCN enhances the ability of WLS to promote WNT1 trafficking to the cell surface as well as secretion, but decreases the ability of WLS to activate WNT1 signaling in target cell. These observations suggest that the levels of WNT1 on the cell surface and in the media are not the sole determinants of the activation of Wnt signaling in target cells.

Quantitative genetics of CTCF binding reveal local sequence effects and different modes of X-chromosome association.

Associating genetic variation with quantitative measures of gene regulation offers a way to bridge the gap between genotype and complex phenotypes. In order to identify quantitative trait loci (QTLs) that influence the binding of a transcription factor in humans, we measured binding of the multifunctional transcription and chromatin factor CTCF in 51 HapMap cell lines. We identified thousands of QTLs in which genotype differences were associated with differences in CTCF binding strength, hundreds of them confirmed by directly observable allele-specific binding bias. The majority of QTLs were either within 1 kb of the CTCF binding motif, or in linkage disequilibrium with a variant within 1 kb of the motif. On the X chromosome we observed three classes of binding sites: a minority class bound only to the active copy of the X chromosome, the majority class bound to both the active and inactive X, and a small set of female-specific CTCF sites associated with two non-coding RNA genes. In sum, our data reveal extensive genetic effects on CTCF binding, both direct and indirect, and identify a diversity of patterns of CTCF binding on the X chromosome.

Dendritic cell fate is determined by BCL11A.

The plasmacytoid dendritic cell (pDC) is vital to the coordinated action of innate and adaptive immunity. pDC development has not been unequivocally traced, nor has its transcriptional regulatory network been fully clarified. Here we confirm an essential requirement for the BCL11A transcription factor in fetal pDC development, and demonstrate this lineage-specific requirement in the adult organism. Furthermore, we identify BCL11A gene targets and provide a molecular mechanism for its action in pDC commitment. Embryonic germ-line deletion of Bcl11a revealed an absolute cellular, molecular, and functional absence of pDCs in fetal mice. In adults, deletion of Bcl11a in hematopoietic stem cells resulted in perturbed yet continued generation of progenitors, loss of downstream pDC and B-cell lineages, and persisting myeloid, conventional dendritic, and T-cell lineages. Challenge with virus resulted in a marked reduction of antiviral response in conditionally deleted adults. Genome-wide analyses of BCL11A DNA binding and expression revealed that BCL11A regulates transcription of E2-2 and other pDC differentiation modulators, including ID2 and MTG16. Our results identify BCL11A as an essential, lineage-specific factor that regulates pDC development, supporting a model wherein differentiation into pDCs represents a primed "default" pathway for common dendritic cell progenitors.

Simultaneous mapping of transcript ends at single-nucleotide resolution and identification of widespread promoter-associated non-coding RNA governed by TATA elements.

Understanding the relationships between regulatory factor binding, chromatin structure, cis-regulatory elements and RNA-regulation mechanisms relies on accurate information about transcription start sites (TSS) and polyadenylation sites (PAS). Although several approaches have identified transcript ends in yeast, limitations of resolution and coverage have remained, and definitive identification of TSS and PAS with single-nucleotide resolution has not yet been achieved. We developed SMORE-seq (simultaneous mapping of RNA ends by sequencing) and used it to simultaneously identify the strongest TSS for 5207 (90%) genes and PAS for 5277 (91%) genes. The new transcript annotations identified by SMORE-seq showed improved distance relationships with TATA-like regulatory elements, nucleosome positions and active RNA polymerase. We found 150 genes whose TSS were downstream of the annotated start codon, and additional analysis of evolutionary conservation and ribosome footprinting suggests that these protein-coding sequences are likely to be mis-annotated. SMORE-seq detected short non-coding RNAs transcribed divergently from more than a thousand promoters in wild-type cells under normal conditions. These divergent non-coding RNAs were less evident at promoters containing canonical TATA boxes, suggesting a model where transcription initiation at promoters by RNAPII is bidirectional, with TATA elements serving to constrain the directionality of initiation.

Host-rabies virus protein-protein interactions as druggable antiviral targets.

We present an unconventional approach to antiviral drug discovery, which is used to identify potent small molecules against rabies virus. First, we conceptualized viral capsid assembly as occurring via a host-catalyzed biochemical pathway, in contrast to the classical view of capsid formation by self-assembly. This suggested opportunities for antiviral intervention by targeting previously unappreciated catalytic host proteins, which were pursued. Second, we hypothesized these host proteins to be components of heterogeneous, labile, and dynamic multi-subunit assembly machines, not easily isolated by specific target protein-focused methods. This suggested the need to identify active compounds before knowing the precise protein target. A cell-free translation-based small molecule screen was established to recreate the hypothesized interactions involving newly synthesized capsid proteins as host assembly machine substrates. Hits from the screen were validated by efficacy against infectious rabies virus in mammalian cell culture. Used as affinity ligands, advanced analogs were shown to bind a set of proteins that effectively reconstituted drug sensitivity in the cell-free screen and included a small but discrete subfraction of cellular ATP-binding cassette family E1 (ABCE1), a host protein previously found essential for HIV capsid formation. Taken together, these studies advance an alternate view of capsid formation (as a host-catalyzed biochemical pathway), a different paradigm for drug discovery (whole pathway screening without knowledge of the target), and suggest the existence of labile assembly machines that can be rendered accessible as next-generation drug targets by the means described.

Technical note: The role of circulating low-density lipoprotein levels as a phenotypic marker for Holstein cholesterol deficiency in dairy cattle.

With the recent discovery of a Holstein cholesterol deficiency (HCD) haplotype, the USDA has labeled many dairy animals as HCD carriers based on haplotype and pedigree analysis. We set out to investigate the effect of HCD status on various cholesterol transport molecules, namely low-density lipoprotein cholesterol (LDL-C), high-density lipoprotein cholesterol, and triglycerides in both males and females. A genome-wide association study was also conducted to narrow down the genomic region correlated with varying LDL-C levels. In the study, 34 HCD carrier animals showed significantly lower cholesterol and LDL-C levels compared with their 34 closely related, non-HCD controls. The genome-wide association study based on 73 animals using 56,198 SNP markers revealed an association with chromosome 11 in the region of 66,218,925 to 66,946,746 bp. We also tested the effect of HCD status on sperm quality traits using fresh ejaculates and frozen-thawed semen samples, but did not find any discriminating effects. Our study has demonstrated the use of LDL-C as a key phenotypic marker for determining HCD status in dairy cattle and this is the first study that clearly shows a cause-effect relationship of the HCD haplotype on circulating LDL-C.

Genetic reconstruction of a functional transcriptional regulatory network.

Although global analyses of transcription factor binding provide one view of potential transcriptional regulatory networks, regulation also occurs at levels distinct from transcription factor binding. Here, we use a genetic approach to identify targets of transcription factors in yeast and reconstruct a functional regulatory network. First, we profiled transcriptional responses in S. cerevisiae strains with individual deletions of 263 transcription factors. Then we used directed-weighted graph modeling and regulatory epistasis analysis to identify indirect regulatory relationships between these transcription factors, and from this we reconstructed a functional transcriptional regulatory network. The enrichment of promoter motifs and Gene Ontology annotations provide insight into the biological functions of the transcription factors.