Search

Previous genome-scale phylogenetic analyses of Fungi have under sampled taxa from Zoopagales; this order contains many predacious or parasitic genera, and most have never been grown in pure culture. We sequenced the genomes of 4 zoopagalean taxa that are predators of amoebae, nematodes, or rotifers and the genome of one taxon that is a parasite of amoebae using single cell sequencing methods with whole genome amplification. Each genome was a metagenome, which was assembled and binned using multiple techniques to identify the target genomes. We inferred phylogenies with both super matrix and coalescent approaches using 192 conserved proteins mined from the target genomes and performed ancestral state reconstructions to determine the ancestral trophic lifestyle of the clade. Our results indicate that Zoopagales is monophyletic. Ancestral state reconstructions provide moderate support for mycoparasitism being the ancestral state of the clade.

Verticillium dahliae is a broad host-range pathogen that causes vascular wilts in plants. Interactions between three hosts and specific V. dahliae genotypes result in severe defoliation. The underlying mechanisms of defoliation are unresolved. Genome resequencing, gene deletion and complementation, gene expression analysis, sequence divergence, defoliating phenotype identification, virulence analysis, and quantification of V. dahliae secondary metabolites were performed. Population genomics previously revealed that G-LSR2 was horizontally transferred from the fungus Fusarium oxysporum f. sp. vasinfectum to V. dahliae and is exclusively found in the genomes of defoliating (D) strains. Deletion of seven genes within G-LSR2, designated as VdDf genes, produced the nondefoliation phenotype on cotton, olive, and okra but complementation of two genes restored the defoliation phenotype. Genes VdDf5 and VdDf6 associated with defoliation shared homology with polyketide synthases involved in secondary metabolism, whereas VdDf7 shared homology with proteins involved in the biosynthesis of N-lauroylethanolamine (N-acylethanolamine (NAE) 12:0), a compound that induces defoliation. NAE overbiosynthesis by D strains also appears to disrupt NAE metabolism in cotton by inducing overexpression of fatty acid amide hydrolase. The VdDfs modulate the synthesis and overproduction of secondary metabolites, such as NAE 12:0, that cause defoliation either by altering abscisic acid sensitivity, hormone disruption, or sensitivity to the pathogen.

To prevent new infections with human immunodeficiency virus type 1 (HIV-1) in sub-Saharan Africa, UNAIDS recommends targeting interventions to populations that are at high risk of acquiring and passing on the virus. Yet it is often unclear who and where these 'source' populations are. Here we demonstrate how viral deep-sequencing can be used to reconstruct HIV-1 transmission networks and to infer the direction of transmission in these networks. We are able to deep-sequence virus from a large population-based sample of infected individuals in Rakai District, Uganda, reconstruct partial transmission networks, and infer the direction of transmission within them at an estimated error rate of 16.3% [8.8-28.3%]. With this error rate, deep-sequence phylogenetics cannot be used against individuals in legal contexts, but is sufficiently low for population-level inferences into the sources of epidemic spread. The technique presents new opportunities for characterizing source populations and for targeting of HIV-1 prevention interventions in Africa.

R-loops are abundant three-stranded nucleic-acid structures that form in cis during transcription. Experimental evidence suggests that R-loop formation is affected by DNA sequence and topology. However, the exact manner by which these factors interact to determine R-loop susceptibility is unclear. To investigate this, we developed a statistical mechanical equilibrium model of R-loop formation in superhelical DNA. In this model, the energy involved in forming an R-loop includes four terms-junctional and base-pairing energies and energies associated with superhelicity and with the torsional winding of the displaced DNA single strand around the RNA:DNA hybrid. This model shows that the significant energy barrier imposed by the formation of junctions can be overcome in two ways. First, base-pairing energy can favor RNA:DNA over DNA:DNA duplexes in favorable sequences. Second, R-loops, by absorbing negative superhelicity, partially or fully relax the rest of the DNA domain, thereby returning it to a lower energy state. In vitro transcription assays confirmed that R-loops cause plasmid relaxation and that negative superhelicity is required for R-loops to form, even in a favorable region. Single-molecule R-loop footprinting following in vitro transcription showed a strong agreement between theoretical predictions and experimental mapping of stable R-loop positions and further revealed the impact of DNA topology on the R-loop distribution landscape. Our results clarify the interplay between base sequence and DNA superhelicity in controlling R-loop stability. They also reveal R-loops as powerful and reversible topology sinks that cells may use to nonenzymatically relieve superhelical stress during transcription.

Liposarcomas (LPSs) are a group of malignant mesenchymal tumors showing adipocytic differentiation. Here, to gain insight into the enhancer dysregulation and transcriptional addiction in this disease, we chart super-enhancer structures in both LPS tissues and cell lines. We identify a bromodomain and extraterminal (BET) protein-cooperated FUS-DDIT3 function in myxoid LPS and a BET protein-dependent core transcriptional regulatory circuitry consisting of FOSL2, MYC, and RUNX1 in de-differentiated LPS. Additionally, SNAI2 is identified as a crucial downstream target that enforces both proliferative and metastatic potentials to de-differentiated LPS cells. Genetic depletion of BET genes, core transcriptional factors, or SNAI2 mitigates consistently LPS malignancy. We also reveal a compelling susceptibility of LPS cells to BET protein degrader ARV-825. BET protein depletion confers additional advantages to circumvent acquired resistance to Trabectedin, a chemotherapy drug for LPS. Moreover, this study provides a framework for discovering and targeting of core oncogenic transcriptional programs in human cancers.

Allele-specific protein-RNA binding is an essential aspect that may reveal functional genetic variants (GVs) mediating post-transcriptional regulation. Recently, genome-wide detection of in vivo binding of RNA-binding proteins is greatly facilitated by the enhanced crosslinking and immunoprecipitation (eCLIP) method. We developed a new computational approach, called BEAPR, to identify allele-specific binding (ASB) events in eCLIP-Seq data. BEAPR takes into account crosslinking-induced sequence propensity and variations between replicated experiments. Using simulated and actual data, we show that BEAPR largely outperforms often-used count analysis methods. Importantly, BEAPR overcomes the inherent overdispersion problem of these methods. Complemented by experimental validations, we demonstrate that the application of BEAPR to ENCODE eCLIP-Seq data of 154 proteins helps to predict functional GVs that alter splicing or mRNA abundance. Moreover, many GVs with ASB patterns have known disease relevance. Overall, BEAPR is an effective method that helps to address the outstanding challenge of functional interpretation of GVs.

Hypercholesterolemia, the driving force of atherosclerosis, accelerates the expansion and mobilization of hematopoietic stem and progenitor cells (HSPCs). The molecular determinants connecting hypercholesterolemia with hematopoiesis are unclear. Here, we report that a somite-derived prohematopoietic cue, AIBP, orchestrates HSPC emergence from the hemogenic endothelium, a type of specialized endothelium manifesting hematopoietic potential. Mechanistically, AIBP-mediated cholesterol efflux activates endothelial Srebp2, the master transcription factor for cholesterol biosynthesis, which in turn transactivates Notch and promotes HSPC emergence. Srebp2 inhibition impairs hypercholesterolemia-induced HSPC expansion. Srebp2 activation and Notch up-regulation are associated with HSPC expansion in hypercholesterolemic human subjects. Genome-wide chromatin immunoprecipitation followed by sequencing (ChIP-seq), RNA sequencing (RNA-seq), and assay for transposase-accessible chromatin using sequencing (ATAC-seq) indicate that Srebp2 transregulates Notch pathway genes required for hematopoiesis. Our studies outline an AIBP-regulated Srebp2-dependent paradigm for HSPC emergence in development and HPSC expansion in atherosclerotic cardiovascular disease.

Large structural variants (SVs) in the human genome are difficult to detect and study by conventional sequencing technologies. With long-range genome analysis platforms, such as optical mapping, one can identify large SVs (>2 kb) across the genome in one experiment. Analyzing optical genome maps of 154 individuals from the 26 populations sequenced in the 1000 Genomes Project, we find that phylogenetic population patterns of large SVs are similar to those of single nucleotide variations in 86% of the human genome, while ~2% of the genome has high structural complexity. We are able to characterize SVs in many intractable regions of the genome, including segmental duplications and subtelomeric, pericentromeric, and acrocentric areas. In addition, we discover ~60 Mb of non-redundant genome content missing in the reference genome sequence assembly. Our results highlight the need for a comprehensive set of alternate haplotypes from different populations to represent SV patterns in the genome.

Caenorhabditis elegans provides an amenable system to explore whether newly composed ribosomes are required to progress through development. Despite the complex pattern of tissues that are formed during embryonic development, we found that null homozygotes lacking any of the five different ribosomal proteins (RPs) can produce fully functional first-stage larvae, with similar developmental competence seen upon complete deletion of the multi-copy ribosomal RNA locus. These animals, relying on maternal but not zygotic contribution of ribosomal components, are capable of completing embryogenesis. In the absence of new ribosomal components, the resulting animals are arrested before progression from the first larval stage and fail in two assays for postembryonic plasticity of neuronal structure. Mosaic analyses of larvae that are a mixture of ribosome-competent and non-competent cells suggest a global regulatory mechanism in which ribosomal insufficiency in a subset of cells triggers organism-wide growth arrest.

The fragile X premutation is a CGG trinucleotide repeat expansion between 55 and 200 repeats in the 5'-untranslated region of the fragile X mental retardation 1 (FMR1) gene. Human carriers of the premutation allele are at risk of developing the late-onset neurodegenerative disorder, fragile X-associated tremor/ataxia syndrome (FXTAS). Characteristic neuropathology associated with FXTAS includes intranuclear inclusions in neurons and astroglia. Previous studies recapitulated these histopathological features in neurons in a knock-in mouse model, but without significant astroglial pathology. To determine the role of astroglia in FXTAS, we generated a transgenic mouse line (Gfa2-CGG99-eGFP) that selectively expresses a 99-CGG repeat expansion linked to an enhanced green fluorescent protein (eGFP) reporter in astroglia throughout the brain, including cerebellar Bergmann glia. Behaviorally these mice displayed impaired motor performance on the ladder-rung test, but paradoxically better performance on the rotarod. Immunocytochemical analysis revealed that CGG99-eGFP co-localized with GFAP and S-100ß, but not with NeuN, Iba1, or MBP, indicating that CGG99-eGFP expression is specific to astroglia. Ubiquitin-positive intranuclear inclusions were found in eGFP-expressing glia throughout the brain. In addition, intracytoplasmic ubiquitin-positive inclusions were found outside the nucleus in distal astrocyte processes. Intriguingly, intranuclear inclusions, in the absence of eGFP mRNA and eGFP fluorescence, were present in neurons of the hypothalamus and neocortex. Furthermore, intranuclear inclusions in both neurons and astrocytes displayed immunofluorescent labeling for the polyglycine peptide FMRpolyG, implicating FMRpolyG in the pathology found in Gfa2-CGG99 mice. Considered together, these results show that Gfa2-CGG99 expression in mice is sufficient to induce key features of FXTAS pathology, including formation of intranuclear inclusions, translation of FMRpolyG, and deficits in motor function.

Cancer cells develop mechanisms to escape immunosurveillance, among which modulating the expression of immune suppressive messenger RNAs is most well-documented. However, how this is molecularly achieved remains largely unresolved. Here, we develop an in vivo mouse model of liver cancer to study oncogene cooperation in immunosurveillance. We show that MYC overexpression (MYCTg) synergizes with KRASG12D to induce an aggressive liver tumor leading to metastasis formation and reduced mouse survival compared with KRASG12D alone. Genome-wide ribosomal footprinting of MYCTg;KRASG12 tumors compared with KRASG12D revealed potential alterations in translation of mRNAs, including programmed-death-ligand 1 (PD-L1). Further analysis revealed that PD-L1 translation is repressed in KRASG12D tumors by functional, non-canonical upstream open reading frames in its 5' untranslated region, which is bypassed in MYCTg;KRASG12D tumors to evade immune attack. We show that this mechanism of PD-L1 translational upregulation was effectively targeted by a potent, clinical compound that inhibits eIF4E phosphorylation, eFT508, which reverses the aggressive and metastatic characteristics of MYCTg;KRASG12D tumors. Together, these studies reveal how immune-checkpoint proteins are manipulated by distinct oncogenes at the level of mRNA translation, which can be exploited for new immunotherapies.

Genome-wide association studies have identified 40 ovarian cancer risk loci. However, the mechanisms underlying these associations remain elusive. In this study, we conducted a two-pronged approach to identify candidate causal SNPs and assess underlying biological mechanisms at chromosome 9p22.2, the first and most statistically significant associated locus for ovarian cancer susceptibility. Three transcriptional regulatory elements with allele-specific effects and a scaffold/matrix attachment region were characterized and, through physical DNA interactions, BNC2 was established as the most likely target gene. We determined the consensus binding sequence for BNC2 in vitro, verified its enrichment in BNC2 ChIP-seq regions, and validated a set of its downstream target genes. Fine-mapping by dense regional genotyping in over 15,000 ovarian cancer cases and 30,000 controls identified SNPs in the scaffold/matrix attachment region as among the most likely causal variants. This study reveals a comprehensive regulatory landscape at 9p22.2 and proposes a likely mechanism of susceptibility to ovarian cancer. SIGNIFICANCE: Mapping the 9p22.2 ovarian cancer risk locus identifies BNC2 as an ovarian cancer risk gene.See related commentary by Choi and Brown, p. 439.

Mechanisms by which members of the AP-1 family of transcription factors play non-redundant biological roles despite recognizing the same DNA sequence remain poorly understood. To address this question, here we investigate the molecular functions and genome-wide DNA binding patterns of AP-1 family members in primary and immortalized mouse macrophages. ChIP-sequencing shows overlapping and distinct binding profiles for each factor that were remodeled following TLR4 ligation. Development of a machine learning approach that jointly weighs hundreds of DNA recognition elements yields dozens of motifs predicted to drive factor-specific binding profiles. Machine learning-based predictions are confirmed by analysis of the effects of mutations in genetically diverse mice and by loss of function experiments. These findings provide evidence that non-redundant genomic locations of different AP-1 family members in macrophages largely result from collaborative interactions with diverse, locus-specific ensembles of transcription factors and suggest a general mechanism for encoding functional specificities of their common recognition motif.

Despite 25 years' effort, serious questions remain regarding the mechanism(s) underlying electron transfer through (or from) electrode-bound double-stranded DNA. In part this is because a control experiment regarding the putatively critical role of guanine bases in the most widely proposed transport mechanism (hopping from guanine to guanine through the π-stack) appears to be lacking from the prior literature. In response, we have employed chronoamperometry, which allows for high-precision determination of electron transfer rates, to characterize transfer to a redox reporter appended onto electrode-bound DNA duplexes. Specifically, we have measured the effects of guanines and base mismatches on the electron transfer rate associated with such constructs. Upon doing so, we find that, counter to prior reports, the transfer rate is, to within relatively tight experimental confidence intervals, unaffected by either. Parallel studies of the dependence of the electron transfer rate on the length of the DNA suggest that transfer from this system obeys a "collision" mechanism in which the redox reporter physically contacts the electrode surface prior to the exchange of electrons.

Transcriptional pausing underlies regulation of cellular RNA biogenesis. A consensus pause sequence that acts on RNA polymerases (RNAPs) from bacteria to mammals halts RNAP in an elemental paused state from which longer-lived pauses can arise. Although the structural foundations of pauses prolonged by backtracking or nascent RNA hairpins are recognized, the fundamental mechanism of the elemental pause is less well-defined. Here we report a mechanistic dissection that establishes the elemental pause signal (i) is multipartite; (ii) causes a modest conformational shift that puts γ-proteobacterial RNAP in an off-pathway state in which template base loading but not RNA translocation is inhibited; and (iii) allows RNAP to enter pretranslocated and one-base-pair backtracked states easily even though the half-translocated state observed in paused cryo-EM structures rate-limits pause escape. Our findings provide a mechanistic basis for the elemental pause and a framework to understand how pausing is modulated by sequence, cellular conditions, and regulators.

Recent studies suggest that closely related species can accumulate substantial genetic and phenotypic differences despite ongoing gene flow, thus challenging traditional ideas regarding the genetics of speciation. Baboons (genus Papio) are Old World monkeys consisting of six readily distinguishable species. Baboon species hybridize in the wild, and prior data imply a complex history of differentiation and introgression. We produced a reference genome assembly for the olive baboon (Papio anubis) and whole-genome sequence data for all six extant species. We document multiple episodes of admixture and introgression during the radiation of Papio baboons, thus demonstrating their value as a model of complex evolutionary divergence, hybridization, and reticulation. These results help inform our understanding of similar cases, including modern humans, Neanderthals, Denisovans, and other ancient hominins.

Rationale: Lung cancer is the leading cause of cancer death worldwide, and treatment options are limited to mainly cytotoxic agents. Here we reveal a novel role of miR-150 in non-small cell lung cancer (NSCLC) development and seek potential new therapeutic targets. Methods: The miR-150-mediated autophagy dysfunction in NSCLC cells were examined using molecular methods in vitro and in vivo. The upstream regulatory element and downstream target of miR-150 were identified in vitro and validated in vivo. Potential therapeutic methods (anti-c-myc or anti-miR-150) were tested in vitro and in vivo. Clinical relevance of the c-myc/miR-150/EPG5 axis in NSCLC was validated in human clinical samples and large genomics database. Results: miR-150 blocked the fusion of autophagosomes and lysosomes through directly repressing EPG5. The miR-150-mediated autophagy defect induced ER stress and increased cellular ROS levels and DNA damage response, and promoted NSCLC cell proliferation and tumor growth. Knockdown of EPG5 promoted NSCLC cell proliferation, and attenuated the effects of miR-150. c-myc gene was identified as a miR-150 transcriptional factor which increased miR-150 accumulation, therefore pharmacologically or genetically inhibiting c-myc/miR-150 expression significantly inhibited NSCLC cell growth in vitro and in vivo. Both c-myc and miR-150 were significantly over-expressed in NSCLC, while EPG5 was down-regulated in NSCLC. Expression levels of these molecules were well correlated, and also well correlated with patient survival. Conclusions: Our findings suggest that c-myc/miR-150/EPG5 mediated dysfunction of autophagy contributes to NSCLC development, which may provide a potential new diagnostic and therapeutic target in NSCLC.

Presence of the integrated endogenous banana streak virus (eBSV) in the B genome of plantain (AAB) is a major challenge for breeding and dissemination of hybrids. As the eBSV activates into infectious viral particles under stress, the progenitor Musa balbisiana and its derivants, having at least one B genome, cannot be used as parents for crop improvement. Here, we report a strategy to inactivate the eBSV by editing the virus sequences. The regenerated genome-edited events of Gonja Manjaya showed mutations in the targeted sites with the potential to prevent proper transcription or/and translational into functional viral proteins. Seventy-five percent of the edited events remained asymptomatic in comparison to the non-edited control plants under water stress conditions, confirming inactivation of eBSV into infectious viral particles. This study paves the way for the improvement of B genome germplasm and its use in breeding programs to produce hybrids that can be globally disseminated.

Intron sequences are removed from newly synthesized RNA and usually rapidly degraded. However, it now seems that introns have a surprising role — helping yeast cells survive when nutrients are scarce.

Given limited resources for managing invasive species, traditional survey methods may not be feasible to implement at a regional scale. Environmental DNA (eDNA) sampling has proven to be an effective method for detecting some invasive species, but comparisons between the detection probability of eDNA and traditional survey methods using modern occupancy modeling methods are rare. We developed a qPCR assay to detect two species of watersnake (Nerodia fasciata and Nerodia sipedon) introduced to California, USA, and we compared the efficacy of eDNA and aquatic trapping. We tested 3-9 water samples each from 30 sites near the known range of N. fasciata, and 61 sites near the known range of N. sipedon. We also deployed aquatic funnel traps at a subset of sites for each species. We detected N. fasciata eDNA in three of nine water samples from just one site, but captured N. fasciata in traps at three of ten sites. We detected N. sipedon eDNA in five of six water samples from one site, which was also the only site of nine at which this species was captured in traps. Traditional trapping surveys had a higher probability of detecting watersnakes than eDNA surveys, and both survey methods had higher detection probability for N. sipedon than N. fasciata. Occupancy models that integrated both trapping and eDNA surveys estimated that 5 sites (95% Credible Interval: 4-10) of 91 were occupied by watersnakes (both species combined), although snakes were only detected at four sites (three for N. fasciata, one for N. sipedon). Our study shows that despite the many successes of eDNA surveys, traditional sampling methods can have higher detection probability for some species. We recommend those tasked with managing species invasions explicitly compare eDNA and traditional survey methods in an occupancy framework to inform their choice of the best method for detecting nascent populations.

A persistent concern with CRISPR-Cas9 gene editing has been the potential to generate mutations at off-target genomic sites. While CRISPR-engineering mice to delete a ~360 bp intronic enhancer, here we discovered a founder line that had marked immune dysregulation caused by a 24 kb tandem duplication of the sequence adjacent to the on-target deletion. Our results suggest unintended repair of on-target genomic cuts can cause pathogenic "bystander" mutations that escape detection by routine targeted genotyping assays.

Alternative splicing diversifies mRNA transcripts in human cells. While the spliceosome pairs exons with a high degree of accuracy, the rates of rare aberrant and non-canonical pre-mRNA splicing have not been evaluated at the nucleotide level to determine the quantity and identity of these events across splice junctions. Using ultra-deep sequencing the frequency of aberrant and non-canonical splicing events for three splice junctions flanking exon 7 of SMN1 were determined at single nucleotide resolution. After correction for background noise introduced by PCR amplification and sequencing steps, pre-mRNA splicing was shown to maintain a low overall rate of aberrant and non-canonically spliced events. Several previously unannotated splicing events across 3 exon|intron junctions in SMN1 were identified. Mutations within SMN exon 7 were shown to affect splicing fidelity by modulating RNA secondary structures, by altering the binding site of regulatory proteins and by changing the 5' splice site strength. Mutations also create a truncated SMN1 exon 7 through the introduction of a de novo non-canonical 5' splice site. The results from the ultra-deep sequencing approach highlight the impressive fidelity of pre-mRNA splicing and demonstrate that the immediate sequence context around splice sites is the main driving force behind non-canonical splice site pairing.

BACKGROUND:Osteomyelitis is a bone infection, most often caused by Staphylococcus aureus, in which neutrophils play a key role. Cathepsin G (CTSG) is a bactericidal serine protease stored in the neutrophil azurophilic granules. CTSG regulates inflammation, activating matrix metalloproteinases (MMPs), and coagulation. Lactoferrin (LF), a neutrophil glycoprotein, increases CTSG catalytic activity and induces inflammation. The aim of this study was to analyze a potential association between a CTSG gene polymorphism (Asn125Ser or N125S, rs45567233), that modifies CTSG activity, and could affect susceptibility to, or outcome of, bacterial osteomyelitis. METHODS:CTSG N125S polymorphism was genotyped in 329 osteomyelitis patients and 415 controls), Blood coagulation parameters, serum CTSG activity, LF, MMP-1, MMP-13, and soluble receptor activator for nuclear factor κ B ligand (sRANKL) levels were assessed in carriers of the different CTSG genotypes. RESULTS:CTSG N125S (AG) genotype was significantly more frequent among osteomyelitis patients than controls (15.5% vs. 9.4%, p = 0.014). CTSG N125S variant G allele (AG +GG) was also more frequent among osteomyelitis patients (8.1% vs. 4.7%, p = 0.01). Serum CTSG activity and LF levels were significantly higher in osteomyelitis patients carrying the G allele compared to those with the AA genotype, (p

Transcriptomes consist of several classes of RNA that have wide-ranging but often poorly described functions and the deregulation of which leads to numerous diseases. Engineering of functionalized RNA-binding proteins (RBPs) could therefore have many applications. Our previous studies suggested that the RanBP2-type Zinc Finger (ZF) domain is a suitable scaffold to investigate the design of single-stranded RBPs. In the present work, we have analyzed the natural sequence specificity of various members of the RanBP2-type ZF family and characterized the interaction with their target RNA. Surprisingly, our data showed that natural RanBP2-type ZFs with different RNA-binding residues exhibit a similar sequence specificity and therefore no simple recognition code can be established. Despite this finding, different discriminative abilities were observed within the family. In addition, in order to target a long RNA sequence and therefore gain in specificity, we generated a 6-ZF array by combining ZFs from the RanBP2-type family but also from different families, in an effort to achieve a wider target sequence repertoire. We showed that this chimeric protein recognizes its target sequence (20 nucleotides), both in vitro and in living cells. Altogether, our results indicate that the use of ZFs in RBP design remains attractive even though engineering of specificity changes is challenging.

Clustering homologous sequences based on their similarity is a problem that appears in many bioinformatics applications. The fact that sequences cluster is ultimately the result of their phylogenetic relationships. Despite this observation and the natural ways in which a tree can define clusters, most applications of sequence clustering do not use a phylogenetic tree and instead operate on pairwise sequence distances. Due to advances in large-scale phylogenetic inference, we argue that tree-based clustering is under-utilized. We define a family of optimization problems that, given an arbitrary tree, return the minimum number of clusters such that all clusters adhere to constraints on their heterogeneity. We study three specific constraints, limiting (1) the diameter of each cluster, (2) the sum of its branch lengths, or (3) chains of pairwise distances. These three problems can be solved in time that increases linearly with the size of the tree, and for two of the three criteria, the algorithms have been known in the theoretical computer scientist literature. We implement these algorithms in a tool called TreeCluster, which we test on three applications: OTU clustering for microbiome data, HIV transmission clustering, and divide-and-conquer multiple sequence alignment. We show that, by using tree-based distances, TreeCluster generates more internally consistent clusters than alternatives and improves the effectiveness of downstream applications. TreeCluster is available at https://github.com/niemasd/TreeCluster.

Methyl-CpG-binding-Protein 2 (MeCP2) is an abundant nuclear protein highly enriched in neurons. Here we report live-cell single-molecule imaging studies of the kinetic features of mouse MeCP2 at high spatial-temporal resolution. MeCP2 displays dynamic features that are distinct from both highly mobile transcription factors and immobile histones. Stable binding of MeCP2 in living neurons requires its methyl-binding domain and is sensitive to DNA modification levels. Diffusion of unbound MeCP2 is strongly constrained by weak, transient interactions mediated primarily by its AT-hook domains, and varies with the level of chromatin compaction and cell type. These findings extend previous studies of the role of the MeCP2 MBD in high affinity DNA binding to living neurons, and identify a new role for its AT-hooks domains as critical determinants of its kinetic behavior. They suggest that limited nuclear diffusion of MeCP2 in live neurons contributes to its local impact on chromatin structure and gene expression.

Many strains of the spirochete Leptospira interrogans serovar Pomona express the osmotically inducible sphingomyelinase gene sph2 at much higher levels than strains from other serovars. We developed a new green fluorescent protein (GFP) reporter plasmid to examine sph2 gene expression determinants. The vector enables the fusion of the test promoter to the ribosome-binding site and coding region of gfp We fused the sph2 promoters from the L. interrogans serovar Lai strain 56601 and from the L. interrogans serovar Pomona strain LC82-25 to gfp to examine the molecular determinants of differential sph2 expression between the two strains. Similar to what was observed with the native sph2 genes, the introduction of the plasmids into the Lai 56601 strain resulted in near background levels of gfp expression from the Lai sph2 promoter, while the expression from the Pomona sph2 promoter was high. The expression of both fusions increased at physiologic levels of osmolarity achieved by adding sodium chloride to the culture medium. We examined the role of a 17-bp upstream element found in all L. interrogans strains expressing low basal levels of sph2 and missing from Pomona strains that express sph2 at high levels. When the 17-bp sequence present upstream of the Lai sph2 promoter was deleted or scrambled, the fusion expression increased substantially. Conversely, the insertion of the 17-bp sequence upstream of the Pomona sph2 promoter diminished fusion expression. In contrast, the removal of an insertion sequence-like element that is found only in the Pomona sph2 upstream sequence had no effect on the expression from the Pomona sph2 fusion in the Lai strain. These findings demonstrate the utility of the gfp reporter plasmid in analyzing gene expression in L. interrogansIMPORTANCE Genetic tools are needed to examine gene expression in the pathogen Leptospira interrogans We developed a reporter plasmid that replicates in L. interrogans with green fluorescent protein (GFP) as the readout of promoter activity. We demonstrated an application of the new reporter plasmid by identifying an upstream element responsible for the poor basal expression of the sph2 sphingomyelinase gene in an L. interrogans serovar Lai strain. This new tool is useful for the discovery of the molecular determinants of L. interrogans gene expression.

BackgroundA remarkable exception to the large genetic diversity often observed for bacteriophages infecting a specific bacterial host was found for the Cutibacterium acnes (formerly Propionibacterium acnes) phages, which are highly homogeneous. Phages infecting the related species, which is also a member of the Propionibacteriaceae family, Propionibacterium freudenreichii, a bacterium used in production of Swiss-type cheeses, have also been described and are common contaminants of the cheese manufacturing process. However, little is known about their genetic composition and diversity.ResultsWe obtained seven independently isolated bacteriophages that infect P. freudenreichii from Swiss-type cheese samples, and determined their complete genome sequences. These data revealed that all seven phage isolates are of similar genomic length and GC% content, but their genomes are highly diverse, including genes encoding the capsid, tape measure, and tail proteins. In contrast to C. acnes phages, all P. freudenreichii phage genomes encode a putative integrase protein, suggesting they are capable of lysogenic growth. This is supported by the finding of related prophages in some P. freudenreichii strains. The seven phages could further be distinguished as belonging to two distinct genomic types, or 'clusters', based on nucleotide sequences, and host range analyses conducted on a collection of P. freudenreichii strains show a higher degree of host specificity than is observed for the C. acnes phages.ConclusionsOverall, our data demonstrate P. freudenreichii bacteriophages are distinct from C. acnes phages, as evidenced by their higher genetic diversity, potential for lysogenic growth, and more restricted host ranges. This suggests substantial differences in the evolution of these related species from the Propionibacteriaceae family and their phages, which is potentially related to their distinct environmental niches.

OBJECTIVES:Thousands of cases of multidrug-resistant tuberculosis (TB) have been observed in the Philippines, but studies on the Mycobacterium tuberculosis (MTB) genotypes that underlie the observed drug resistance profiles are lacking. This study aimed to analyse the whole genomes of clinical MTB isolates representing various resistance profiles to identify single nucleotide polymorphisms (SNPs) in resistance-associated genes. METHODS:The genomes of ten MTB isolates cultured from banked sputum sources were sequenced. Bioinformatics analysis consisted of assembly, annotation and SNP identification in genes reported to be associated with resistance to isoniazid (INH), rifampicin (RIF), ethambutol (ETH), streptomycin, pyrazinamide (PZA) and fluoroquinolones (FQs). RESULTS:The draft assemblies covered an average of 97.08% of the expected genome size. Seven of the ten isolates belonged to the Indo-Oceanic lineage/EA12-Manila clade. Two isolates were classified into the Euro-American lineage, whilst the pre-XDR (pre-extensively drug-resistant) isolate was classified under the East Asian/Beijing clade. The SNPs katG Ser315Thr, rpoB Ser450Leu and embB Met306Val were found in INH- (4/7), RIF- (3/6) and ETH-resistant (2/6) isolates, respectively, but not in susceptible isolates. Mutations in the inhA promoter and in the pncA and gyrA genes known to be involved in resistance to INH, PZA and FQs, respectively, were also identified. CONCLUSIONS:This study represents the first effort to investigate the whole genomes of Philippine clinical strains of MTB exhibiting various multidrug resistance profiles. Whole-genome data can provide valuable insights to the mechanistic and epidemiological qualities of TB in a high-burden setting such as the Philippines.

Aspergillus section Nigri comprises filamentous fungi relevant to biomedicine, bioenergy, health, and biotechnology. To learn more about what genetically sets these species apart, as well as about potential applications in biotechnology and biomedicine, we sequenced 23 genomes de novo, forming a full genome compendium for the section (26 species), as well as 6 Aspergillus niger isolates. This allowed us to quantify both inter- and intraspecies genomic variation. We further predicted 17,903 carbohydrate-active enzymes and 2,717 secondary metabolite gene clusters, which we condensed into 455 distinct families corresponding to compound classes, 49% of which are only found in single species. We performed metabolomics and genetic engineering to correlate genotypes to phenotypes, as demonstrated for the metabolite aurasperone, and by heterologous transfer of citrate production to Aspergillus nidulans. Experimental and computational analyses showed that both secondary metabolism and regulation are key factors that are significant in the delineation of Aspergillus species.

BackgroundThe genomes of halophilic archaea (haloarchaea) often comprise multiple replicons. Genomic variation in haloarchaea has been linked to viral infection pressure and, in the case of Antarctic communities, can be caused by intergenera gene exchange. To expand understanding of genome variation and biogeography of Antarctic haloarchaea, here we assessed genomic variation between two strains of Halorubrum lacusprofundi that were isolated from Antarctic hypersaline lakes from different regions (Vestfold Hills and Rauer Islands). To assess variation in haloarchaeal populations, including the presence of genomic islands, metagenomes from six hypersaline Antarctic lakes were characterised.ResultsThe sequence of the largest replicon of each Hrr. lacusprofundi strain (primary replicon) was highly conserved, while each of the strains' two smaller replicons (secondary replicons) were highly variable. Intergenera gene exchange was identified, including the sharing of a type I-B CRISPR system. Evaluation of infectivity of an Antarctic halovirus provided experimental evidence for the differential susceptibility of the strains, bolstering inferences that strain variation is important for modulating interactions with viruses. A relationship was found between genomic structuring and the location of variation within replicons and genomic islands, demonstrating that the way in which haloarchaea accommodate genomic variability relates to replicon structuring. Metagenome read and contig mapping and clustering and scaling analyses demonstrated biogeographical patterning of variation consistent with environment and distance effects. The metagenome data also demonstrated that specific haloarchaeal species dominated the hypersaline systems indicating they are endemic to Antarctica.ConclusionThe study describes how genomic variation manifests in Antarctic-lake haloarchaeal communities and provides the basis for future assessments of Antarctic regional and global biogeography of haloarchaea.

A recent study in Genome Biology has characterized the evolution of candidate hominoid-specific liver enhancers by using massively parallel reporter assays (MPRAs).

BackgroundGlaciers cover ~ 10% of land but are among the least explored environments on Earth. The basal portion of glaciers often harbors unique aquatic microbial ecosystems in the absence of sunlight, and knowledge on the microbial community structures and their metabolic potential is very limited. Here, we provide insights into the microbial lifestyle present at the base of the Matanuska Glacier, Alaska.ResultsDNA and RNA were extracted from samples of the Matanuska Glacier basal ice. Using Illumina MiSeq and HiSeq sequencing, we investigated the microbial diversity with the metagenomic shotgun reads and 16S ribosomal RNA data. We further assembled 9 partial and draft bacterial genomes from the metagenomic assembly, and identified key metabolic pathways such as sulfur oxidation and nitrification. Collectively, our analyses suggest a prevalence of lithotrophic and heterotrophic metabolisms in the subglacial microbiome.ConclusionOur results present the first metagenomic assembly and bacterial draft genomes for a subglacial environment. These results extend our understanding of the chemical and biological processes in subglacial environments critically influenced by global climate change.

BackgroundBase editors are a class of genome editing tools with the ability to efficiently induce point mutations in genomic DNA, without inducing double-strand breaks or relying on homology-direct repair as in other such technologies. Recently, adenine base editors (ABEs) have been developed to mediate the conversion of A•T to G•C in genomic DNA of human cells, mice, and plants. Here, we investigated the activity and efficiency of several adenine base editors in zebrafish and showed that base editing can be used to create new models of pathogenic diseases caused by point mutations.ResultsThe original ABE7.10 exhibits almost no activity in zebrafish. After codon optimization, we found that a zABE7.10 variant could induce targeted conversion of adenine to guanine in zebrafish at multiple tested genomic loci, and all the target sites showed a high rate of germline targeting efficiency. Furthermore, using this system, we established a zebrafish model of 5q-Syndrome that contained a new point mutation in rps14. The further modification of zABE7.10 by a bipartite nuclear localization signals (bpNLS) resulted in 1.96-fold average improvement in ABE-mediated editing efficiency at four sites.ConclusionsCollectively, this system, designated as zABE7.10, provides a strategy to perform A•T to G•C base editing in zebrafish and enhances its capacity to model human diseases.

BackgroundTaxonomic classification of marker-gene sequences is an important step in microbiome analysis.ResultsWe present q2-feature-classifier ( https://github.com/qiime2/q2-feature-classifier ), a QIIME 2 plugin containing several novel machine-learning and alignment-based methods for taxonomy classification. We evaluated and optimized several commonly used classification methods implemented in QIIME 1 (RDP, BLAST, UCLUST, and SortMeRNA) and several new methods implemented in QIIME 2 (a scikit-learn naive Bayes machine-learning classifier, and alignment-based taxonomy consensus methods based on VSEARCH, and BLAST+) for classification of bacterial 16S rRNA and fungal ITS marker-gene amplicon sequence data. The naive-Bayes, BLAST+-based, and VSEARCH-based classifiers implemented in QIIME 2 meet or exceed the species-level accuracy of other commonly used methods designed for classification of marker gene sequences that were evaluated in this work. These evaluations, based on 19 mock communities and error-free sequence simulations, including classification of simulated "novel" marker-gene sequences, are available in our extensible benchmarking framework, tax-credit ( https://github.com/caporaso-lab/tax-credit-data ).ConclusionsOur results illustrate the importance of parameter tuning for optimizing classifier performance, and we make recommendations regarding parameter choices for these classifiers under a range of standard operating conditions. q2-feature-classifier and tax-credit are both free, open-source, BSD-licensed packages available on GitHub.

BackgroundThe neonatal intensive care unit (NICU) contains a unique cohort of patients with underdeveloped immune systems and nascent microbiome communities. Patients often spend several months in the same room, and it has been previously shown that the gut microbiomes of these infants often resemble the microbes found in the NICU. Little is known, however, about the identity, persistence, and absolute abundance of NICU room-associated bacteria over long stretches of time. Here, we couple droplet digital PCR (ddPCR), 16S rRNA gene surveys, and recently published metagenomics data from infant gut samples to infer the extent to which the NICU microbiome is shaped by its room occupants.ResultsOver 2832 swabs, wipes, and air samples were collected from 16 private-style NICU rooms housing very low birth weight (

BackgroundAdvances in genomic technologies have expanded our ability to accurately and exhaustively detect natural genomic variants that can be applied in crop improvement and to increase our knowledge of plant evolution and adaptation. Switchgrass (Panicum virgatum L.), an allotetraploid (2n = 4× = 36) perennial C4 grass (Poaceae family) native to North America and a feedstock crop for cellulosic biofuel production, has a large potential for genetic improvement due to its high genotypic and phenotypic variation. In this study, we analyzed single nucleotide polymorphism (SNP) variation in 372 switchgrass genotypes belonging to 36 accessions for 12 genes putatively involved in biomass production to investigate signatures of selection that could have led to ecotype differentiation and to population adaptation to geographic zones.ResultsA total of 11,682 SNPs were mined from ~ 15 Gb of sequence data, out of which 251 SNPs were retained after filtering. Population structure analysis largely grouped upland accessions into one subpopulation and lowland accessions into two additional subpopulations. The most frequent SNPs were in homozygous state within accessions. Sixty percent of the exonic SNPs were non-synonymous and, of these, 45% led to non-conservative amino acid changes. The non-conservative SNPs were largely in linkage disequilibrium with one haplotype being predominantly present in upland accessions while the other haplotype was commonly present in lowland accessions. Tajima's test of neutrality indicated that PHYB, a gene involved in photoperiod response, was under positive selection in the switchgrass population. PHYB carried a SNP leading to a non-conservative amino acid change in the PAS domain, a region that acts as a sensor for light and oxygen in signal transduction.ConclusionsSeveral non-conservative SNPs in genes potentially involved in plant architecture and adaptation have been identified and led to population structure and genetic differentiation of ecotypes in switchgrass. We suggest here that PHYB is a key gene involved in switchgrass natural selection. Further analyses are needed to determine whether any of the non-conservative SNPs identified play a role in the differential adaptation of upland and lowland switchgrass.

The spliceosome catalyzes the removal of introns from pre-messenger RNA (mRNA) and subsequent pairing of exons with remarkable fidelity. Some exons are known to be skipped or included in the mature mRNA in a cell type- or context-dependent manner (cassette exons), thereby contributing to the diversification of the human proteome. Interestingly, splicing is initiated (and sometimes completed) co-transcriptionally. Here, we develop a kinetic mathematical modeling framework to investigate alternative co-transcriptional splicing (CTS) and, specifically, the control of cassette exons' inclusion. We show that when splicing is co-transcriptional, default splice patterns of exon inclusion are more likely than when splicing is post-transcriptional, and that certain exons are more likely to be regulatable (i.e. cassette exons) than others, based on the exon-intron structure context. For such regulatable exons, transcriptional elongation rates may affect splicing outcomes. Within the CTS paradigm, we examine previously described hypotheses of co-operativity between splice sites of short introns (i.e. 'intron definition') or across short exons (i.e. 'exon definition'), and find that models encoding these faithfully recapitulate observations in the fly and human genomes, respectively.

The group II intron and the spliceosome share a common active site architecture and are thought to be evolutionarily related. Here we report the 3.7 Å crystal structure of a eukaryotic group II intron in the lariat-3' exon form, immediately preceding the second step of splicing, analogous to the spliceosomal P complex. This structure reveals the location of the intact 3' splice site within the catalytic core of the group II intron. The 3'-OH of the 5' exon is positioned in close proximity to the 3' splice site for nucleophilic attack and exon ligation. The active site undergoes conformational rearrangements with the catalytic triplex having different configurations before and after the second step of splicing. We describe a complete model for the second step of group II intron splicing that incorporates a dynamic catalytic triplex being responsible for creating the binding pocket for 3' splice site capture.

Cap-independent translation is believed to play an important role in eukaryotic protein synthesis, but the mechanisms of ribosomal recruitment and translation initiation remain largely unknown. Messenger RNA display was previously used to profile the human genome for RNA leader sequences that can enhance cap-independent translation. Surprisingly, many of the isolated sequences contain AUG triplets, suggesting a possible functional role for these motifs during translation initiation. Herein, we examine the sequence determinants of AUG triplets within a set of human translation enhancing elements (TEEs). Functional analyses performed in vitro and in cultured cells indicate that AUGs have the capacity to modulate mRNA translation either by serving as part of a larger ribosomal recruitment site or by directing the ribosome to defined initiation sites. These observations help constrain the functional role of AUG triplets in human TEEs and advance our understanding of this specific mechanism of cap-independent translation initiation.

Chloroviruses (family Phycodnaviridae) are dsDNA viruses found throughout the world's inland waters. The open reading frames in the genomes of 41 sequenced chloroviruses (330 ± 40 kbp each) representing three virus types were analyzed for evidence of evolutionarily conserved local genomic "contexts", the organization of biological information into units of a scale larger than a gene. Despite a general loss of synteny between virus types, we informatically detected a highly conserved genomic context defined by groups of three or more genes that we have termed "gene gangs". Unlike previously described local genomic contexts, the definition of gene gangs requires only that member genes be consistently co-localized and are not constrained by strand, regulatory sites, or intervening sequences (and therefore represent a new type of conserved structural genomic element). An analysis of functional annotations and transcriptomic data suggests that some of the gene gangs may organize genes involved in specific biochemical processes, but that this organization does not involve their coordinated expression.

The DNA base excision repair (BER) glycosylase MUTYH prevents DNA mutations by catalyzing adenine (A) excision from inappropriately formed 8-oxoguanine (8-oxoG):A mismatches. The importance of this mutation suppression activity in tumor suppressor genes is underscored by the association of inherited variants of MUTYH with colorectal polyposis in a hereditary colorectal cancer syndrome known as MUTYH-associated polyposis, or MAP. Many of the MAP variants encompass amino acid changes that occur at positions surrounding the two-metal cofactor-binding sites of MUTYH. One of these cofactors, found in nearly all MUTYH orthologs, is a [4Fe-4S]2+ cluster coordinated by four Cys residues located in the N-terminal catalytic domain. We recently uncovered a second functionally relevant metal cofactor site present only in higher eukaryotic MUTYH orthologs: a Zn2+ ion coordinated by three Cys residues located within the extended interdomain connector (IDC) region of MUTYH that connects the N-terminal adenine excision and C-terminal 8-oxoG recognition domains. In this work, we identified a candidate for the fourth Zn2+ coordinating ligand using a combination of bioinformatics and computational modeling. In addition, using in vitro enzyme activity assays, fluorescence polarization DNA binding assays, circular dichroism spectroscopy, and cell-based rifampicin resistance assays, the functional impact of reduced Zn2+ chelation was evaluated. Taken together, these results illustrate the critical role that the "Zn2+ linchpin motif" plays in MUTYH repair activity by providing for proper engagement of the functional domains on the 8-oxoG:A mismatch required for base excision catalysis. The functional importance of the Zn2+ linchpin also suggests that adjacent MAP variants or exposure to environmental chemicals may compromise Zn2+ coordination, and ability of MUTYH to prevent disease.

Approximately 1500 RNA-binding proteins (RBPs) profoundly impact mammalian cellular function by controlling distinct sets of transcripts, often using sequence-specific binding to 3' untranslated regions (UTRs) to regulate mRNA stability and translation. Aside from their individual effects, higher-order combinatorial interactions between RBPs on specific mRNAs have been proposed to underpin the regulatory network. To assess the extent of such co-regulatory control, we took a global experimental approach followed by targeted validation to examine interactions between two well-characterized and highly conserved RBPs, Argonaute2 (AGO2) and Pumilio (PUM1 and PUM2). Transcriptome-wide changes in AGO2-mRNA binding upon PUM knockdown were quantified by CLIP-seq, and the presence of PUM binding on the same 3'UTR corresponded with cooperative and antagonistic effects on AGO2 occupancy. In addition, PUM binding sites that overlap with AGO2 showed differential, weakened binding profiles upon abrogation of AGO2 association, indicative of cooperative interactions. In luciferase reporter validation of candidate 3'UTR sites where AGO2 and PUM colocalized, three sites were identified to host antagonistic interactions, where PUM counteracts miRNA-guided repression. Interestingly, the binding sites for the two proteins are too far for potential antagonism due to steric hindrance, suggesting an alternate mechanism. Our data experimentally confirms the combinatorial regulatory model and indicates that the mostly repressive PUM proteins can change their behavior in a context-dependent manner. Overall, the approach underscores the importance of further elucidation of complex interactions between RBPs and their transcriptome-wide extent.

Self-organization of discrete fates in human gastruloids is mediated by a hierarchy of signaling pathways. How these pathways are integrated in time, and whether cells maintain a memory of their signaling history remains obscure. Here, we dissect the temporal integration of two key pathways, WNT and ACTIVIN, which along with BMP control gastrulation. CRISPR/Cas9-engineered live reporters of SMAD1, 2 and 4 demonstrate that in contrast to the stable signaling by SMAD1, signaling and transcriptional response by SMAD2 is transient, and while necessary for pluripotency, it is insufficient for differentiation. Pre-exposure to WNT, however, endows cells with the competence to respond to graded levels of ACTIVIN, which induces differentiation without changing SMAD2 dynamics. This cellular memory of WNT signaling is necessary for ACTIVIN morphogen activity. A re-evaluation of the evidence gathered over decades in model systems, re-enforces our conclusions and points to an evolutionarily conserved mechanism.

Recruitment of poliovirus (PV) RNA to the human ribosome requires the coordinated interaction of the viral internal ribosome entry site (IRES) and several host cellular initiation factors and IRES trans-acting factors (ITAFs). Attenuated PV Sabin strains contain point mutations in the PV IRES domain V (dV) that inhibit viral translation. Remarkably, attenuation is most apparent in cells of the central nervous system, but the molecular basis to explain this is poorly understood. The dV contains binding sites for eukaryotic initiation factor 4G (eIF4G) and polypyrimidine tract-binding protein (PTB). Impaired binding of these proteins to the mutant IRESs has been observed, but these effects have not been quantitated. We used a fluorescence anisotropy assay to reveal that the Sabin mutants reduce the equilibrium dissociation constants of eIF4G and PTB to the PV IRES by up to 6-fold. Using the most inhibitory Sabin 3 mutant, we used a real-time fluorescence helicase assay to show that the apparent affinity of an active eIF4G/4A/4B helicase complex for the IRES is reduced by 2.5-fold. The Sabin 3 mutant did not alter the maximum rate of eIF4A-dependent helicase activity, suggesting that this mutant primarily reduces the affinity, rather than activity, of the unwinding complex. To confirm this affinity model of attenuation, we show that eIF4G overexpression in HeLa cells overcomes the attenuation of a Sabin 3 mutant PV-luciferase replicon. Our study provides a quantitative framework for understanding the mechanism of PV Sabin attenuation and provides an explanation for the previously observed cell type-specific translational attenuation.

Adipocytes arise from the commitment and differentiation of adipose precursors in white adipose tissue (WAT). In studying adipogenesis, precursor markers, including Pref-1 and PDGFRα, are used to isolate precursors from stromal vascular fractions of WAT, but the relation among the markers is not known. Here, we used the Pref-1 promoter-rtTA system in mice for labeling Pref-1+ cells and for inducible inactivation of the Pref-1 target Sox9. We show the requirement of Sox9 for the maintenance of Pref-1+ proliferative, early precursors. Upon Sox9 inactivation, these Pref-1+ cells become PDGFRα+ cells that express early adipogenic markers. Thus, we show that Pref-1+ cells precede PDGFRα+ cells in the adipogenic pathway and that Sox9 inactivation is required for WAT growth and expansion. Furthermore, we show that in maintaining early adipose precursors, Sox9 activates Meis1, which prevents adipogenic differentiation. Our study also demonstrates the Pref-1 promoter-rtTA system for inducible gene inactivation in early adipose precursor populations.

Hippo signaling has been recognized as a key tumor suppressor pathway. Here, we perform a comprehensive molecular characterization of 19 Hippo core genes in 9,125 tumor samples across 33 cancer types using multidimensional "omic" data from The Cancer Genome Atlas. We identify somatic drivers among Hippo genes and the related microRNA (miRNA) regulators, and using functional genomic approaches, we experimentally characterize YAP and TAZ mutation effects and miR-590 and miR-200a regulation for TAZ. Hippo pathway activity is best characterized by a YAP/TAZ transcriptional target signature of 22 genes, which shows robust prognostic power across cancer types. Our elastic-net integrated modeling further reveals cancer-type-specific pathway regulators and associated cancer drivers. Our results highlight the importance of Hippo signaling in squamous cell cancers, characterized by frequent amplification of YAP/TAZ, high expression heterogeneity, and significant prognostic patterns. This study represents a systems-biology approach to characterizing key cancer signaling pathways in the post-genomic era.

DNA methylation functions in gene silencing and the maintenance of genome integrity. In plants, non-CG DNA methylation is linked through a self-reinforcing loop with histone 3 lysine 9 dimethylation (H3K9me2). The plant-specific SUPPRESSOR OF VARIEGATION 3-9 HOMOLOG (SUVH) family H3K9 methyltransferases (MTases) bind to DNA methylation marks and catalyze H3K9 methylation. Here, we analyzed the structure and function of Arabidopsis thaliana SUVH6 to understand how this class of enzyme maintains methylation patterns in the genome. We reveal that SUVH6 has a distinct 5-methyl-dC (5mC) base-flipping mechanism involving a thumb loop element. Autoinhibition of H3 substrate entry is regulated by a SET domain loop, and a conformational transition in the post-SET domain upon cofactor binding may control catalysis. In vitro DNA binding and in vivo ChIP-seq data reveal that the different SUVH family H3K9 MTases have distinct DNA binding preferences, targeting H3K9 methylation to sites with different methylated DNA sequences, explaining the context biased non-CG DNA methylation in plants.

Telomerase adds telomeric repeats to chromosome ends by processive copying of a template within the telomerase RNA bound to telomerase reverse transcriptase. Telomerase RNAs have single-stranded regions that separate the template from a 5' stem and 3' pseudoknot, and mammals gained additional stem P2a.1 separating the template from the pseudoknot. Using human telomerase, we show that the length of template 3'-flanking single-stranded RNA is a determinant of repeat addition processivity whereas template 5'-flanking single-stranded RNA and P2a.1 are critical for activity but not processivity. In comparison, requirements for the template sequence itself are confounding: different substitutions of the same position have strikingly different consequences, from improved processivity and activity to complete inactivation. We discovered that some altered-template sequences stabilize an alternative RNA conformation that precludes the pseudoknot by base-pairing of one pseudoknot strand to the template 3' end. Using mutations to reduce over-stability of the alternative conformation, we restore high activity and processivity to otherwise inactive altered-template telomerase ribonucleoproteins. In cells, over-stabilization or destabilization of the alternative state severely inhibited biogenesis of active telomerase. Our findings delineate roles for human telomerase RNA template-flanking regions, establish a biologically relevant pseudoknot-alternative RNA conformation, and expand the repertoire of human telomerase repeat synthesis.

The growth of human immunodeficiency virus (HIV) sequence databases resulting from drug resistance testing has motivated efforts using phylogenetic methods to assess how HIV spreads1-4. Such inference is potentially both powerful and useful for tracking the epidemiology of HIV and the allocation of resources to prevention campaigns. We recently used simulation and a small number of illustrative cases to show that certain phylogenetic patterns are associated with different types of epidemiological linkage5. Our original approach was later generalized for large next-generation sequencing datasets and implemented as a free computational pipeline6. Previous work has claimed that direction and directness of transmission could not be established from phylogeny because one could not be sure that there were no intervening or missing links involved7-9. Here, we address this issue by investigating phylogenetic patterns from 272 previously identified HIV transmission chains with 955 transmission pairs representing diverse geography, risk groups, subtypes, and genomic regions. These HIV transmissions had known linkage based on epidemiological information such as partner studies, mother-to-child transmission, pairs identified by contact tracing, and criminal cases. We show that the resulting phylogeny inferred from real HIV genetic sequences indeed reveals distinct patterns associated with direct transmission contra transmissions from a common source. Thus, our results establish how to interpret phylogenetic trees based on HIV sequences when tracking who-infected-whom, when and how genetic information can be used for improved tracking of HIV spread. We also investigate limitations that stem from limited sampling and genetic time-trends in the donor and recipient HIV populations.