Search

Biosynthesis of glycosylphosphatidylinositol (GPI) is required for anchoring proteins to the plasma membrane, and is essential for the integrity of the fungal cell wall. Here, we use a reporter gene-based screen in Saccharomyces cerevisiae for the discovery of antifungal inhibitors of GPI-anchoring of proteins, and identify the oligocyclopropyl-containing natural product jawsamycin (FR-900848) as a potent hit. The compound targets the catalytic subunit Spt14 (also referred to as Gpi3) of the fungal UDP-glycosyltransferase, the first step in GPI biosynthesis, with good selectivity over the human functional homolog PIG-A. Jawsamycin displays antifungal activity in vitro against several pathogenic fungi including Mucorales, and in vivo in a mouse model of invasive pulmonary mucormycosis due to Rhyzopus delemar infection. Our results provide a starting point for the development of Spt14 inhibitors for treatment of invasive fungal infections.

BackgroundThe giant squid (Architeuthis dux; Steenstrup, 1857) is an enigmatic giant mollusc with a circumglobal distribution in the deep ocean, except in the high Arctic and Antarctic waters. The elusiveness of the species makes it difficult to study. Thus, having a genome assembled for this deep-sea-dwelling species will allow several pending evolutionary questions to be unlocked.FindingsWe present a draft genome assembly that includes 200 Gb of Illumina reads, 4 Gb of Moleculo synthetic long reads, and 108 Gb of Chicago libraries, with a final size matching the estimated genome size of 2.7 Gb, and a scaffold N50 of 4.8 Mb. We also present an alternative assembly including 27 Gb raw reads generated using the Pacific Biosciences platform. In addition, we sequenced the proteome of the same individual and RNA from 3 different tissue types from 3 other species of squid (Onychoteuthis banksii, Dosidicus gigas, and Sthenoteuthis oualaniensis) to assist genome annotation. We annotated 33,406 protein-coding genes supported by evidence, and the genome completeness estimated by BUSCO reached 92%. Repetitive regions cover 49.17% of the genome.ConclusionsThis annotated draft genome of A. dux provides a critical resource to investigate the unique traits of this species, including its gigantism and key adaptations to deep-sea environments.

Following photosynthesis, sucrose is translocated to sink organs, where it provides the primary source of carbon and energy to sustain plant growth and development. Sugar transporters from the SWEET (sugar will eventually be exported transporter) family are rate-limiting factors that mediate sucrose transport across concentration gradients, sustain yields, and participate in reproductive development, plant senescence, stress responses, as well as support plant-pathogen interaction, the focus of this study. We identified 25 SWEET genes in the walnut genome and distinguished each by its individual gene structure and pattern of expression in different walnut tissues. Their chromosomal locations, cis-acting motifs within their 5' regulatory elements, and phylogenetic relationship patterns provided the first comprehensive analysis of the SWEET gene family of sugar transporters in walnut. This family is divided into four clades, the analysis of which suggests duplication and expansion of the SWEET gene family in Juglans regia. In addition, tissue-specific gene expression signatures suggest diverse possible functions for JrSWEET genes. Although these are commonly used by pathogens to harness sugar products from their plant hosts, little was known about their role during Xanthomonasarboricola pv. juglandis (Xaj) infection. We monitored the expression profiles of the JrSWEET genes in different tissues of "Chandler" walnuts when challenged with pathogen Xaj417 and concluded that SWEET-mediated sugar translocation from the host is not a trigger for walnut blight disease development. This may be directly related to the absence of type III secretion system-dependent transcription activator-like effectors (TALEs) in Xaj417, which suggests different strategies are employed by this pathogen to promote susceptibility to this major aboveground disease of walnuts.

Ciliates, a eukaryotic clade that is over 1 billion years old, are defined by division of genome function between transcriptionally inactive germline micronuclei and functional somatic macronuclei. To date, most analyses of gene family evolution have been limited to cultivable model lineages (e.g., Tetrahymena, Paramecium, Oxytricha, and Stylonychia). Here, we focus on the uncultivable Karyorelictea and its understudied sister class Heterotrichea, which represent two extremes in genome architecture. Somatic macronuclei within the Karyorelictea are described as nearly diploid, while the Heterotrichea have hyperpolyploid somatic genomes. Previous analyses indicate that genome architecture impacts ciliate gene family evolution as the most diverse and largest gene families are found in lineages with extensively processed somatic genomes (i.e., possessing thousands of gene-sized chromosomes). To further assess ciliate gene family evolution, we analyzed 43 single-cell transcriptomes from 33 ciliate species representing 10 classes. Focusing on conserved eukaryotic genes, we use estimates of transcript diversity as a proxy for the number of paralogs in gene families among four focal clades: Karyorelictea, Heterotrichea, extensive fragmenters (with gene-size somatic chromosomes), and non-extensive fragmenters (with more traditional somatic chromosomes), the latter two within the subphylum Intramacronucleata. Our results show that (i) the Karyorelictea have the lowest average transcript diversity, while Heterotrichea are highest among the four groups; (ii) proteins in Karyorelictea are under the highest functional constraints, and the patterns of selection in ciliates may reflect genome architecture; and (iii) stop codon reassignments vary among members of the Heterotrichea and Spirotrichea but are conserved in other classes.IMPORTANCE To further our understanding of genome evolution in eukaryotes, we assess the relationship between patterns of molecular evolution within gene families and variable genome structures found among ciliates. We combine single-cell transcriptomics with bioinformatic tools, focusing on understudied and uncultivable lineages selected from across the ciliate tree of life. Our analyses show that genome architecture correlates with patterns of protein evolution as lineages with more canonical somatic genomes, such as the class Karyorelictea, have more conserved patterns of molecular evolution compared to other classes. This study showcases the power of single-cell transcriptomics for investigating genome architecture and evolution in uncultivable microbial lineages and provides transcriptomic resources for further research on genome evolution.

Equine chromosome 24 microRNA cluster (C24MC), the ortholog of human C14MC, is a pregnancy-related miRNA cluster. This cluster is believed to be implicated in embryonic, fetal, and placental development. The current study aimed to characterize the expression profile of this cluster in maternal circulation throughout equine gestation. The expression profile of miRNAs belonging to this cluster was analyzed in the serum of non-pregnant (diestrus), pregnant (25 d, 45 d, 4 mo, 6 mo, 10 mo), and postpartum mares. Among the miRNAs examined, 11 miRNAs were differentially expressed across the analyzed time-points. Four of these miRNAs (eca-miR-1247-3p, eca-miR-134-5p, eca-miR-382-5p, and eca-miR-433-3p) were found to be enriched in the serum of pregnant mares at Day 25 relative to non-pregnant mares. To further assess the accuracy of these miRNAs in differentiating pregnant (25 d) from non-pregnant mares, receiver operating characteristic (ROC) analysis was performed for each of these miRNAs, revealing that eca-miR-1247-3p and eca-miR-134-5p had the highest accuracy (AUCROC = 0.92 and 0.91, respectively; p < 0.05). Moreover, eca-miR-1247-3p, eca-miR-134-5p, eca-miR-409-3p, and eca-miR-379-5p were enriched in the serum of Day 45 pregnant mares. Among those miRNAs, eca-miR-1247-3p and eca-miR-409-3p retained the highest accuracy as shown by ROC analysis. GO analysis revealed that these miRNAs are mainly implicated in nervous system development as well as organ development. Using in situ hybridization, we localized eca-miR-409-3p in the developing embryo (25 d) and extra-embryonic membranes (25 and 45 d). In conclusion, the present study is the first to elucidate the circulating maternal profile of C24MC-associated miRNAs throughout pregnancy and to suggest that serum eca-miR-1247-3p, eca-miR-134-5p, and eca-miR-409-3p could be used as pregnancy-specific markers during early gestation (25 and 45 d). Overall, the high abundance of these embryo-derived miRNAs in the maternal circulation suggests an embryo-maternal communication during the equine early pregnancy.

Necrotizing enterocolitis (NEC) is a devastating intestinal disease that occurs primarily in premature infants. We performed genome-resolved metagenomic analysis of 1163 fecal samples from premature infants to identify microbial features predictive of NEC. Features considered include genes, bacterial strain types, eukaryotes, bacteriophages, plasmids, and growth rates. A machine learning classifier found that samples collected before NEC diagnosis harbored significantly more Klebsiella, bacteria encoding fimbriae, and bacteria encoding secondary metabolite gene clusters related to quorum sensing and bacteriocin production. Notably, replication rates of all bacteria, especially Enterobacteriaceae, were significantly higher 2 days before NEC diagnosis. The findings uncover biomarkers that could lead to early detection of NEC and targets for microbiome-based therapeutics.

Microbial populations form intricate macroscopic colonies with diverse morphologies whose functions remain to be fully understood. Despite fungal colonies isolated from environmental and clinical samples revealing abundant intraspecies morphological diversity, it is unclear how this diversity affects fungal fitness and disease progression. Here we observe a notable effect of oxygen tension on the macroscopic and biofilm morphotypes of the human fungal pathogen Aspergillus fumigatus. A hypoxia-typic morphotype is generated through the expression of a subtelomeric gene cluster containing genes that alter the hyphal surface and perturb interhyphal interactions to disrupt in vivo biofilm and infection site morphologies. Consequently, this morphotype leads to increased host inflammation, rapid disease progression and mortality in a murine model of invasive aspergillosis. Taken together, these data suggest that filamentous fungal biofilm morphology affects fungal-host interactions and should be taken into consideration when assessing virulence and host disease progression of an isolated strain.

It is generally believed that exchange of secondary metabolite biosynthetic gene clusters (BGCs) among closely related bacteria is an important driver of BGC evolution and diversification. Applying this idea may help researchers efficiently connect many BGCs to their products and characterize the products' roles in various environments. However, existing genetic tools support only a small fraction of these efforts. Here, we present the development of chassis-independent recombinase-assisted genome engineering (CRAGE), which enables single-step integration of large, complex BGC constructs directly into the chromosomes of diverse bacteria with high accuracy and efficiency. To demonstrate the efficacy of CRAGE, we expressed three known and six previously identified but experimentally elusive non-ribosomal peptide synthetase (NRPS) and NRPS-polyketide synthase (PKS) hybrid BGCs from Photorhabdus luminescens in 25 diverse γ-Proteobacteria species. Successful activation of six BGCs identified 22 products for which diversity and yield were greater when the BGCs were expressed in strains closely related to the native strain than when they were expressed in either native or more distantly related strains. Activation of these BGCs demonstrates the feasibility of exploiting their underlying catalytic activity and plasticity, and provides evidence that systematic approaches based on CRAGE will be useful for discovering and identifying previously uncharacterized metabolites.

BackgroundSerratia plymuthica WS3236 was selected for whole genome sequencing based on preliminary genetic and chemical screening indicating the presence of multiple natural product pathways. This led to the identification of a putative sodorifen biosynthetic gene cluster (BGC). The natural product sodorifen is a volatile organic compound (VOC) with an unusual polymethylated hydrocarbon bicyclic structure (C16H26) produced by selected strains of S. plymuthica. The BGC encoding sodorifen consists of four genes, two of which (sodA, sodB) are homologs of genes encoding enzymes of the non-mevalonate pathway and are thought to enhance the amounts of available farnesyl pyrophosphate (FPP), the precursor of sodorifen. Proceeding from FPP, only two enzymes are necessary to produce sodorifen: an S-adenosyl methionine dependent methyltransferase (SodC) with additional cyclisation activity and a terpene-cyclase (SodD). Previous analysis of S. plymuthica found sodorifen production titers are generally low and vary significantly among different producer strains. This precludes studies on the still elusive biological function of this structurally and biosynthetically fascinating bacterial terpene.ResultsSequencing and mining of the S. plymuthica WS3236 genome revealed the presence of 38 BGCs according to antiSMASH analysis, including a putative sodorifen BGC. Further genome mining for sodorifen and sodorifen-like BGCs throughout bacteria was performed using SodC and SodD as queries and identified a total of 28 sod-like gene clusters. Using direct pathway cloning (DiPaC) we intercepted the 4.6 kb candidate sodorifen BGC from S. plymuthica WS3236 (sodA-D) and transformed it into Escherichia coli BL21. Heterologous expression under the control of the tetracycline inducible PtetO promoter firmly linked this BGC to sodorifen production. By utilizing this newly established expression system, we increased the production yields by approximately 26-fold when compared to the native producer. In addition, sodorifen was easily isolated in high purity by simple head-space sampling.ConclusionsGenome mining of all available genomes within the NCBI and JGI IMG databases led to the identification of a wealth of sod-like pathways which may be responsible for producing a range of structurally unknown sodorifen analogs. Introduction of the S. plymuthica WS3236 sodorifen BGC into the fast-growing heterologous expression host E. coli with a very low VOC background led to a significant increase in both sodorifen product yield and purity compared to the native producer. By providing a reliable, high-level production system, this study sets the stage for future investigations of the biological role and function of sodorifen and for functionally unlocking the bioinformatically identified putative sod-like pathways.

We report reference-quality genome assemblies and annotations for two accessions of soybean (Glycine max) and for one accession of Glycine soja, the closest wild relative of G. max. The G. max assemblies provided are for widely used US cultivars: the northern line Williams 82 (Wm82) and the southern line Lee. The Wm82 assembly improves the prior published assembly, and the Lee and G. soja assemblies are new for these accessions. Comparisons among the three accessions show generally high structural conservation, but nucleotide difference of 1.7 single-nucleotide polymorphisms (snps) per kb between Wm82 and Lee, and 4.7 snps per kb between these lines and G. soja. snp distributions and comparisons with genotypes of the Lee and Wm82 parents highlight patterns of introgression and haplotype structure. Comparisons against the US germplasm collection show placement of the sequenced accessions relative to global soybean diversity. Analysis of a pan-gene collection shows generally high conservation, with variation occurring primarily in genomically clustered gene families. We found approximately 40-42 inversions per chromosome between either Lee or Wm82v4 and G. soja, and approximately 32 inversions per chromosome between Wm82 and Lee. We also investigated five domestication loci. For each locus, we found two different alleles with functional differences between G. soja and the two domesticated accessions. The genome assemblies for multiple cultivated accessions and for the closest wild ancestor of soybean provides a valuable set of resources for identifying causal variants that underlie traits for the domestication and improvement of soybean, serving as a basis for future research and crop improvement efforts for this important crop species.

BackgroundClustered protocadherins (PCDHs) map in tandem at human chromosome 5q31 and comprise three multi-genes clusters: α-, β- and γ-PCDH. The expression of this cluster consists of a complex mechanism involving DNA hub formation through DNA-CCTC binding factor (CTCF) interaction. Methylation alterations can affect this interaction, leading to transcriptional dysregulation. In cancer, clustered PCDHs undergo a mechanism of long-range epigenetic silencing by hypermethylation.ResultsIn this study, we detected frequent methylation alterations at CpG islands associated to these clustered PCDHs in all the solid tumours analysed (colorectal, gastric and biliary tract cancers, pilocytic astrocytoma), but not hematologic neoplasms such as chronic lymphocytic leukemia. Importantly, several altered CpG islands were associated with CTCF binding sites. Interestingly, our analysis revealed a hypomethylation event in pilocytic astrocytoma, suggesting that in neuronal tissue, where PCDHs are highly expressed, these genes become hypomethylated in this type of cancer. On the other hand, in tissues where PCDHs are lowly expressed, these CpG islands are targeted by DNA methylation. In fact, PCDH-associated CpG islands resulted hypermethylated in gastrointestinal tumours.ConclusionsOur study highlighted a strong alteration of the clustered PCDHs methylation pattern in the analysed solid cancers and suggested these methylation aberrations in the CpG islands associated with PCDH genes as powerful diagnostic biomarkers.

Objective: To analyse the transcriptional profiles of the pir multigene family of Plasmodium chabaudi chabaudi in male and female gametocytes isolated from the blood of infected mice. Results: Infected red blood cells containing female and male P. chabaudi gametocytes transcribe a distinct set of genes encoded by the multigene family pir. The overall patterns are similar to what has been observed in the close relative P. berghei, but here we show that gametocyte-associated pir genes are distinct from those involved in chronic blood-stage infection and highlight a male-associated pir gene which should be the focus of future studies. [ABSTRACT FROM AUTHOR]

α-Gliadins are a major group of gluten proteins in wheat flour that contribute to the end-use properties for food processing and contain major immunogenic epitopes that can cause serious health-related issues including celiac disease (CD). α-Gliadins are also the youngest group of gluten proteins and are encoded by a large gene family. The majority of the gene family members evolved independently in the A, B, and D genomes of different wheat species after their separation from a common ancestral species. To gain insights into the origin and evolution of these complex genes, the genomic regions of the Gli-2 loci encoding α-gliadins were characterized from the tetraploid wild emmer, a progenitor of hexaploid bread wheat that contributed the AABB genomes. Genomic sequences of Gli-2 locus regions for the wild emmer A and B genomes were first reconstructed using the genome sequence scaffolds along with optical genome maps. A total of 24 and 16 α-gliadin genes were identified for the A and B genome regions, respectively. α-Gliadin pseudogene frequencies of 86% for the A genome and 69% for the B genome were primarily caused by C to T substitutions in the highly abundant glutamine codons, resulting in the generation of premature stop codons. Comparison with the homologous regions from the hexaploid wheat cv. Chinese Spring indicated considerable sequence divergence of the two A genomes at the genomic level. In comparison, conserved regions between the two B genomes were identified that included α-gliadin pseudogenes containing shared nested TE insertions. Analyses of the genomic organization and phylogenetic tree reconstruction indicate that although orthologous gene pairs derived from speciation were present, large portions of α-gliadin genes were likely derived from differential gene duplications or deletions after the separation of the homologous wheat genomes ~ 0.5 MYA. The higher number of full-length intact α-gliadin genes in hexaploid wheat than that in wild emmer suggests that human selection through domestication might have an impact on α-gliadin evolution. Our study provides insights into the rapid and dynamic evolution of genomic regions harboring the α-gliadin genes in wheat.

Genome mining of biosynthetic pathways streamlines discovery of secondary metabolites but can leave ambiguities in the predicted structures, which must be rectified experimentally. Through coupling the reactivity predicted by biosynthetic gene clusters with verified structures, the origin of the β-hydroxyaspartic acid diastereomers in siderophores is reported herein. Two functional subtypes of nonheme Fe(II)/α-ketoglutarate-dependent aspartyl β-hydroxylases are identified in siderophore biosynthetic gene clusters, which differ in genomic organization-existing either as fused domains (IβHAsp) at the carboxyl terminus of a nonribosomal peptide synthetase (NRPS) or as stand-alone enzymes (TβHAsp)-and each directs opposite stereoselectivity of Asp β-hydroxylation. The predictive power of this subtype delineation is confirmed by the stereochemical characterization of β-OHAsp residues in pyoverdine GB-1, delftibactin, histicorrugatin, and cupriachelin. The l-threo (2S, 3S) β-OHAsp residues of alterobactin arise from hydroxylation by the β-hydroxylase domain integrated into NRPS AltH, while l-erythro (2S, 3R) β-OHAsp in delftibactin arises from the stand-alone β-hydroxylase DelD. Cupriachelin contains both l-threo and l-erythro β-OHAsp, consistent with the presence of both types of β-hydroxylases in the biosynthetic gene cluster. A third subtype of nonheme Fe(II)/α-ketoglutarate-dependent enzymes (IβHHis) hydroxylates histidyl residues with l-threo stereospecificity. A previously undescribed, noncanonical member of the NRPS condensation domain superfamily is identified, named the interface domain, which is proposed to position the β-hydroxylase and the NRPS-bound amino acid prior to hydroxylation. Through mapping characterized β-OHAsp diastereomers to the phylogenetic tree of siderophore β-hydroxylases, methods to predict β-OHAsp stereochemistry in silico are realized.

Flooding due to extreme weather threatens crops and ecosystems. To understand variation in gene regulatory networks activated by submergence, we conducted a high-resolution analysis of chromatin accessibility and gene expression at three scales of transcript control in four angiosperms, ranging from a dryland-adapted wild species to a wetland crop. The data define a cohort of conserved submergence-activated genes with signatures of overlapping cis regulation by four transcription factor families. Syntenic genes are more highly expressed than nonsyntenic genes, yet both can have the cis motifs and chromatin accessibility associated with submergence up-regulation. Whereas the flexible circuitry spans the eudicot-monocot divide, the frequency of specific cis motifs, extent of chromatin accessibility, and degree of submergence activation are more prevalent in the wetland crop and may have adaptive importance.

Oncodriver genes are usually identified when mutations recur in multiple tumours. Different drivers often converge in the activation or repression of key cancer-relevant pathways. However, as many pathways contain multiple members of the same gene family, individual mutations might be overlooked, as each family member would necessarily have a lower mutation frequency and thus not identified as significant in any one-gene-at-a-time analysis. Here, we looked for mutated, functional sequence positions in gene families that were mutually exclusive (in patients) with another gene in the same pathway, which identified both known and new candidate oncodrivers. For instance, many inactivating mutations in multiple G-protein (particularly Gi/o) coupled receptors, are mutually exclusive with Gαs oncogenic activating mutations, both of which ultimately enhance cAMP signalling. By integrating transcriptomics and interaction data, we show that the Gs pathway is upregulated in multiple cancer types, even those lacking known GNAS activating mutations. This suggests that cancer cells may develop alternative strategies to activate adenylate cyclase signalling in multiple cancer types. Our study provides a mechanistic interpretation for several rare somatic mutations in multi-gene oncodrivers, and offers possible explanations for known and potential off-label cancer treatments, suggesting new therapeutic opportunities.

Systemic lupus erythematous (SLE) is a heterogeneous autoimmune disease in which outcomes vary among different racial groups. Here, we aim to identify SLE subgroups within a multiethnic cohort using an unsupervised clustering approach based on the American College of Rheumatology (ACR) classification criteria. We identify three patient clusters that vary according to disease severity. Methylation association analysis identifies a set of 256 differentially methylated CpGs across clusters, including 101 CpGs in genes in the Type I Interferon pathway, and we validate these associations in an external cohort. A cis-methylation quantitative trait loci analysis identifies 744 significant CpG-SNP pairs. The methylation signature is enriched for ethnic-associated CpGs suggesting that genetic and non-genetic factors may drive outcomes and ethnic-associated methylation differences. Our computational approach highlights molecular differences associated with clusters rather than single outcome measures. This work demonstrates the utility of applying integrative methods to address clinical heterogeneity in multifactorial multi-ethnic disease settings.

Genome sequencing of environmental bacteria allows identification of biosynthetic gene clusters encoding unusual combinations of enzymes that produce unknown natural products. We identified a pathway in which a ribosomally synthesized small peptide serves as a scaffold for nonribosomal peptide extension and chemical modification. Amino acids are transferred to the carboxyl terminus of the peptide through adenosine triphosphate and amino acyl-tRNA-dependent chemistry that is independent of the ribosome. Oxidative rearrangement, carboxymethylation, and proteolysis of a terminal cysteine yields an amino acid-derived small molecule. Microcrystal electron diffraction demonstrates that the resulting product is isosteric to glutamate. We show that a similar peptide extension is used during the biosynthesis of the ammosamides, which are cytotoxic pyrroloquinoline alkaloids. These results suggest an alternative paradigm for biosynthesis of amino acid-derived natural products.

Among the many genes involved in the metabolism of therapeutic drugs, human arylamine N-acetyltransferases (NATs) genes have been extensively studied, due to their medical importance both in pharmacogenetics and disease epidemiology. One member of this small gene family, NAT2, is established as the locus of the classic human acetylation polymorphism in drug metabolism. Current hypotheses hold that selective processes favoring haplotypes conferring lower NAT2 activity have been operating in modern humans' recent history as an adaptation to local chemical and dietary environments. To shed new light on such hypotheses, we investigated the genetic diversity of the three members of the NAT gene family in seven hominid species, including modern humans, Neanderthals and Denisovans. Little polymorphism sharing was found among hominids, yet all species displayed high NAT diversity, but distributed in an opposite fashion in chimpanzees and bonobos (Pan genus) compared to modern humans, with higher diversity in Pan species at NAT1 and lower at NAT2, while the reverse is observed in humans. This pattern was also reflected in the results returned by selective neutrality tests, which suggest, in agreement with the predicted functional impact of mutations detected in non-human primates, stronger directional selection, presumably purifying selection, at NAT1 in modern humans, and at NAT2 in chimpanzees. Overall, the results point to the evolution of divergent functions of these highly homologous genes in the different primate species, possibly related to their specific chemical/dietary environment (exposome) and we hypothesize that this is likely linked to the emergence of controlled fire use in the human lineage.

E2F transcription factors are central regulators of cell division and cell fate decisions. E2F4 often represents the predominant E2F activity in cells. E2F4 is a transcriptional repressor implicated in cell cycle arrest and whose repressive activity depends on its interaction with members of the RB family. Here we show that E2F4 is important for the proliferation and the survival of mouse embryonic stem cells. In these cells, E2F4 acts in part as a transcriptional activator that promotes the expression of cell cycle genes. This role for E2F4 is independent of the RB family. Furthermore, E2F4 functionally interacts with chromatin regulators associated with gene activation and we observed decreased histone acetylation at the promoters of cell cycle genes and E2F targets upon loss of E2F4 in RB family-mutant cells. Taken together, our findings uncover a non-canonical role for E2F4 that provide insights into the biology of rapidly dividing cells.

In many animal species with a bilateral symmetry, Hox genes are clustered either at one or at several genomic loci. This organization has a functional relevance, as the transcriptional control applied to each gene depends upon its relative position within the gene cluster. It was previously noted that vertebrate Hox clusters display a much higher level of genomic organization than their invertebrate counterparts. The former are always more compact than the latter, they are generally devoid of repeats and of interspersed genes, and all genes are transcribed by the same DNA strand, suggesting that particular factors constrained these clusters toward a tighter structure during the evolution of the vertebrate lineage. Here, we investigate the importance of uniform transcriptional orientation by engineering several alleles within the HoxD cluster, such as to invert one or several transcription units, with or without a neighboring CTCF site. We observe that the association between the tight structure of mammalian Hox clusters and their regulation makes inversions likely detrimental to the proper implementation of this complex genetic system. We propose that the consolidation of Hox clusters in vertebrates, including transcriptional polarity, evolved in conjunction with the emergence of global gene regulation via the flanking regulatory landscapes, to optimize a coordinated response of selected subsets of target genes in cis.

The family Ampullariidae includes both aquatic and amphibious apple snails. They are an emerging model for evolutionary studies due to the high diversity, ancient history, and wide geographical distribution. Insight into drivers of ampullariid evolution is hampered, however, by the lack of genomic resources. Here, we report the genomes of four ampullariids spanning the Old World (Lanistes nyassanus) and New World (Pomacea canaliculata, P. maculata, and Marisa cornuarietis) clades. The ampullariid genomes have conserved ancient bilaterial karyotype features and a novel Hox gene cluster rearrangement, making them valuable in comparative genomic studies. They have expanded gene families related to environmental sensing and cellulose digestion, which may have facilitated some ampullarids to become notorious invasive pests. In the amphibious Pomacea, novel acquisition of an egg neurotoxin and a protein for making the calcareous eggshell may have been key adaptations enabling their transition from underwater to terrestrial egg deposition.

Embryonal tumors with multilayered rosettes (ETMRs) are highly lethal infant brain cancers with characteristic amplification of Chr19q13.41 miRNA cluster (C19MC) and enrichment of pluripotency factor LIN28A. Here we investigated C19MC oncogenic mechanisms and discovered a C19MC-LIN28A-MYCN circuit fueled by multiple complex regulatory loops including an MYCN core transcriptional network and super-enhancers resulting from long-range MYCN DNA interactions and C19MC gene fusions. Our data show that this powerful oncogenic circuit, which entraps an early neural lineage network, is potently abrogated by bromodomain inhibitor JQ1, leading to ETMR cell death.

Hybrid type I PKS/NRPS biosynthetic pathways typically proceed in a collinear manner wherein one molecular building block is enzymatically incorporated in a sequence that corresponds to gene arrangement. In this work, genome mining combined with the use of a fluorogenic azide-based click probe led to the discovery and characterization of vatiamides A-F, three structurally diverse alkynylated lipopeptides, and their brominated analogues, from the cyanobacterium Moorea producens ASI16Jul14-2. These derive from a unique combinatorial non-collinear PKS/NRPS system encoded by a 90 kb gene cluster in which an upstream PKS cassette interacts with three separate cognate NRPS partners. This is facilitated by a series of promiscuous intermodule PKS-NRPS docking motifs possessing identical amino acid sequences. This interaction confers a new type of combinatorial capacity for creating molecular diversity in microbial systems.

Plant endosymbiosis with nitrogen-fixing cyanobacteria has independently evolved in diverse plant lineages, offering a unique window to study the evolution and genetics of plant-microbe interaction. However, very few complete genomes exist for plant cyanobionts, and therefore little is known about their genomic and functional diversity. Here, we present four complete genomes of cyanobacteria isolated from bryophytes. Nanopore long-read sequencing allowed us to obtain circular contigs for all the main chromosomes and most of the plasmids. We found that despite having a low 16S rRNA sequence divergence, the four isolates exhibit considerable genome reorganizations and variation in gene content. Furthermore, three of the four isolates possess genes encoding vanadium (V)-nitrogenase (vnf), which is uncommon among diazotrophs and has not been previously reported in plant cyanobionts. In two cases, the vnf genes were found on plasmids, implying possible plasmid-mediated horizontal gene transfers. Comparative genomic analysis of vnf-containing cyanobacteria further identified a conserved gene cluster. Many genes in this cluster have not been functionally characterized and would be promising candidates for future studies to elucidate V-nitrogenase function and regulation.

Marine sponges are recognized as valuable sources of bioactive metabolites and renowned as petri dishes of the sea, providing specialized niches for many symbiotic microorganisms. Sponges of the family Dysideidae are well documented to be chemically talented, often containing high levels of polyhalogenated compounds, terpenoids, peptides, and other classes of bioactive small molecules. This group of tropical sponges hosts a high abundance of an uncultured filamentous cyanobacterium, Hormoscilla spongeliae Here, we report the comparative genomic analyses of two phylogenetically distinct Hormoscilla populations, which reveal shared deficiencies in essential pathways, hinting at possible reasons for their uncultivable status, as well as differing biosynthetic machinery for the production of specialized metabolites. One symbiont population contains clustered genes for expanded polybrominated diphenylether (PBDE) biosynthesis, while the other instead harbors a unique gene cluster for the biosynthesis of the dysinosin nonribosomal peptides. The hybrid sequencing and assembly approach utilized here allows, for the first time, a comprehensive look into the genomes of these elusive sponge symbionts.IMPORTANCE Natural products provide the inspiration for most clinical drugs. With the rise in antibiotic resistance, it is imperative to discover new sources of chemical diversity. Bacteria living in symbiosis with marine invertebrates have emerged as an untapped source of natural chemistry. While symbiotic bacteria are often recalcitrant to growth in the lab, advances in metagenomic sequencing and assembly now make it possible to access their genetic blueprint. A cell enrichment procedure, combined with a hybrid sequencing and assembly approach, enabled detailed genomic analysis of uncultivated cyanobacterial symbiont populations in two chemically rich tropical marine sponges. These population genomes reveal a wealth of secondary metabolism potential as well as possible reasons for historical difficulties in their cultivation.

Genes encoding cytochrome P450 (CYP; P450) enzymes occur widely in the Archaea, Bacteria, and Eukarya, where they play important roles in metabolism of endogenous regulatory molecules and exogenous chemicals. We now report that genes for multiple and unique P450s occur commonly in giant viruses in the Mimiviridae, Pandoraviridae, and other families in the proposed order Megavirales. P450 genes were also identified in a herpesvirus (Ranid herpesvirus 3) and a phage (Mycobacterium phage Adler). The Adler phage P450 was classified as CYP102L1, and the crystal structure of the open form was solved at 2.5 Å. Genes encoding known redox partners for P450s (cytochrome P450 reductase, ferredoxin and ferredoxin reductase, and flavodoxin and flavodoxin reductase) were not found in any viral genome so far described, implying that host redox partners may drive viral P450 activities. Giant virus P450 proteins share no more than 25% identity with the P450 gene products we identified in Acanthamoeba castellanii, an amoeba host for many giant viruses. Thus, the origin of the unique P450 genes in giant viruses remains unknown. If giant virus P450 genes were acquired from a host, we suggest it could have been from an as yet unknown and possibly ancient host. These studies expand the horizon in the evolution and diversity of the enormously important P450 superfamily. Determining the origin and function of P450s in giant viruses may help to discern the origin of the giant viruses themselves.

Extracellular matrix (ECM) is considered central to the evolution of metazoan multicellularity; however, the repertoire of ECM proteins in nonbilaterians remains unclear. Thrombospondins (TSPs) are known to be well conserved from cnidarians to vertebrates, yet to date have been considered a unique family, principally studied for matricellular functions in vertebrates. Through searches utilizing the highly conserved C-terminal region of TSPs, we identify undisclosed new families of TSP-related proteins in metazoans, designated mega-TSP, sushi-TSP, and poriferan-TSP, each with a distinctive phylogenetic distribution. These proteins share the TSP C-terminal region domain architecture, as determined by domain composition and analysis of molecular models against known structures. Mega-TSPs, the only form identified in ctenophores, are typically >2,700 aa and are also characterized by N-terminal leucine-rich repeats and central cadherin/immunoglobulin domains. In cnidarians, which have a well-defined ECM, Mega-TSP was expressed throughout embryogenesis in Nematostella vectensis, with dynamic endodermal expression in larvae and primary polyps and widespread ectodermal expression in adult Nematostella vectensis and Hydra magnipapillata polyps. Hydra Mega-TSP was also expressed during regeneration and siRNA-silencing of Mega-TSP in Hydra caused specific blockade of head regeneration. Molecular phylogenetic analyses based on the conserved TSP C-terminal region identified each of the TSP-related groups to form clades distinct from the canonical TSPs. We discuss models for the evolution of the newly defined TSP superfamily by gene duplications, radiation, and gene losses from a debut in the last metazoan common ancestor. Together, the data provide new insight into the evolution of ECM and tissue organization in metazoans.

Triacsins are a family of natural products having in common an N-hydroxytriazene moiety not found in any other known secondary metabolites. Though many studies have examined the biological activity of triacsins in lipid metabolism, their biosynthesis has remained unknown. Here we report the identification of the triacsin biosynthetic gene cluster in Streptomyces aureofaciens ATCC 31442. Bioinformatic analysis of the gene cluster led to the discovery of the tacrolimus producer Streptomyces tsukubaensis NRRL 18488 as a new triacsin producer. In addition to targeted gene disruption to identify necessary genes for triacsin production, stable isotope feeding was performed in vivo to advance the understanding of N-hydroxytriazene biosynthesis.

The production of pentabromopseudilin and related brominated compounds by Pseudoalteromonas spp. has recently been linked to the bmp biosynthetic gene cluster. This study explored the distribution and evolutionary history of this gene cluster in the genus Pseudoalteromonas. A phylogeny of the genus revealed numerous clades that do not contain type strains, suggesting considerable species level diversity has yet to be described. Comparative genomics revealed four distinct versions of the gene cluster distributed among 19 of the 101 Pseudoalteromonas genomes examined. These were largely localized to the least inclusive clades containing the Pseudoalteromonas luteoviolacea and Pseudoalteromonas phenolica type strains and show clear evidence of gene and gene cluster loss in certain lineages. Bmp gene phylogeny is largely congruent with the Pseudoalteromonas species phylogeny, suggesting vertical inheritance within the genus. However, the gene cluster is found in three different genomic environments suggesting either chromosomal rearrangement or multiple acquisition events. Bmp conservation within certain lineages suggests the encoded products are highly relevant to the ecology of these bacteria.

The rice XA21-mediated immune response is activated on recognition of the RaxX peptide produced by the bacterium Xanthomonas oryzae pv. oryzae (Xoo). The 60-residue RaxX precursor is post-translationally modified to form a sulfated tyrosine peptide that shares sequence and functional similarity with the plant sulfated tyrosine (PSY) peptide hormones. The 5-kb raxX-raxSTAB gene cluster of Xoo encodes RaxX, the RaxST tyrosylprotein sulfotransferase, and the RaxA and RaxB components of a predicted type I secretion system. To assess raxX-raxSTAB gene cluster evolution and to determine its phylogenetic distribution, we first identified rax gene homologues in other genomes. We detected the complete raxX-raxSTAB gene cluster only in Xanthomonas spp., in five distinct lineages in addition to X. oryzae. The phylogenetic distribution of the raxX-raxSTAB gene cluster is consistent with the occurrence of multiple lateral (horizontal) gene transfer events during Xanthomonas speciation. RaxX natural variants contain a restricted set of missense substitutions, as expected if selection acts to maintain peptide hormone-like function. Indeed, eight RaxX variants tested all failed to activate the XA21-mediated immune response, yet retained peptide hormone activity. Together, these observations support the hypothesis that the XA21 receptor evolved specifically to recognize Xoo RaxX.

Glomeromycotina is a lineage of early diverging fungi that establish arbuscular mycorrhizal (AM) symbiosis with land plants. Despite their major ecological role, the genetic basis of their obligate mutualism remains largely unknown, hindering our understanding of their evolution and biology. We compared the genomes of Glomerales (Rhizophagus irregularis, Rhizophagus diaphanus, Rhizophagus cerebriforme) and Diversisporales (Gigaspora rosea) species, together with those of saprotrophic Mucoromycota, to identify gene families and processes associated with these lineages and to understand the molecular underpinning of their symbiotic lifestyle. Genomic features in Glomeromycotina appear to be very similar with a very high content in transposons and protein-coding genes, extensive duplications of protein kinase genes, and loss of genes coding for lignocellulose degradation, thiamin biosynthesis and cytosolic fatty acid synthase. Most symbiosis-related genes in R. irregularis and G. rosea are specific to Glomeromycotina. We also confirmed that the present species have a homokaryotic genome organisation. The high interspecific diversity of Glomeromycotina gene repertoires, affecting all known protein domains, as well as symbiosis-related orphan genes, may explain the known adaptation of Glomeromycotina to a wide range of environmental settings. Our findings contribute to an increasingly detailed portrait of genomic features defining the biology of AM fungi.

BackgroundGeographical allocation of interventions focusing on hotspots of human immunodeficiency virus (HIV) transmission has the potential to improve efficiency. We used phylogeographic analyses to identify hotspots of the HIV transmission in Cologne-Bonn, Germany.MethodsWe included 714 HIV-1 infected individuals, followed up at the University Hospitals Cologne and Bonn. Distance-based molecular network analyses were performed to infer putative relationships. Characteristics of genetically linked individuals and assortativity (shared characteristics) were analyzed. Geospatial diffusion (ie, viral gene flow) was evaluated using a Slatkin-Maddison approach. Geospatial dispersal was determined by calculating the average distance between the residences of linked individuals (centroids of 3-digit zip code).ResultsIn sum, 217/714 (30.4%) sequences had a putative genetic linkage, forming 77 clusters (size range: 2-8). Linked individuals were more likely to live in areas surrounding the city center (P = .043),

Thousands of natural products have been identified from cultured microorganisms, yet evidence of their production in the environment has proven elusive. Technological advances in mass spectrometry, combined with public databases, now make it possible to address this disparity by detecting compounds directly from environmental samples. Here, we used adsorbent resins, tandem mass spectrometry, and next-generation sequencing to assess the metabolome of marine sediments and its relationship to bacterial community structure. We identified natural products previously reported from cultured bacteria, providing evidence they are produced in situ, and compounds of anthropogenic origin, suggesting this approach can be used as an indicator of environmental impact. The bacterial metabolite staurosporine was quantified and shown to reach physiologically relevant concentrations, indicating that it may influence sediment community structure. Staurosporine concentrations were correlated with the relative abundance of the staurosporine-producing bacterial genus Salinispora and production confirmed in strains cultured from the same location, providing a link between compound and candidate producer. Metagenomic analyses revealed numerous biosynthetic gene clusters related to indolocarbazole biosynthesis, providing evidence for noncanonical sources of staurosporine and a path forward to assess the relationships between natural products and the organisms that produce them. Untargeted environmental metabolomics circumvents the need for laboratory cultivation and represents a promising approach to understanding the functional roles of natural products in shaping microbial community structure in marine sediments.IMPORTANCE Natural products are readily isolated from cultured bacteria and exploited for useful purposes, including drug discovery. However, these compounds are rarely detected in the environments from which the bacteria are obtained, thus limiting our understanding of their ecological significance. Here, we used environmental metabolomics to directly assess chemical diversity in marine sediments. We identified numerous metabolites and, in one case, isolated strains of bacteria capable of producing one of the compounds detected. Coupling environmental metabolomics with community and metagenomic analyses provides opportunities to link compounds and producers and begin to assess the complex interactions mediated by specialized metabolites in marine sediments.

Stochastic activation of clustered Protocadherin (Pcdh) α, β, and γ genes generates a cell-surface identity code in individual neurons that functions in neural circuit assembly. Here, we show that Pcdhα gene choice involves the activation of an antisense promoter located in the first exon of each Pcdhα alternate gene. Transcription of an antisense long noncoding RNA (lncRNA) from this antisense promoter extends through the sense promoter, leading to DNA demethylation of the CTCF binding sites proximal to each promoter. Demethylation-dependent CTCF binding to both promoters facilitates cohesin-mediated DNA looping with a distal enhancer (HS5-1), locking in the transcriptional state of the chosen Pcdhα gene. Uncoupling DNA demethylation from antisense transcription by Tet3 overexpression in mouse olfactory neurons promotes CTCF binding to all Pcdhα promoters, resulting in proximity-biased DNA looping of the HS5-1 enhancer. Thus, antisense transcription-mediated promoter demethylation functions as a mechanism for distance-independent enhancer/promoter DNA looping to ensure stochastic Pcdhα promoter choice.

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.

The ancient reductive acetyl-CoA pathway is employed by acetogenic bacteria to form acetate from inorganic energy sources. Since the central pathway does not gain net ATP by substrate-level phosphorylation, chemolithoautotrophic growth relies on the additional formation of ATP via a chemiosmotic mechanism. Genome analyses indicated that some acetogens only have an energy-converting, ion-translocating hydrogenase (Ech) as a potential respiratory enzyme. Although the Ech-encoding genes are widely distributed in nature, the proposed function of Ech as an ion-translocating chemiosmotic coupling site has neither been demonstrated in bacteria nor has it been demonstrated that it can be the only energetic coupling sites in microorganisms that depend on a chemiosmotic mechanism for energy conservation. Here, we show that the Ech complex of the thermophilic acetogenic bacterium Thermoanaerobacter kivui is indeed a respiratory enzyme. Experiments with resting cells prepared from T. kivui cultures grown on carbon monoxide (CO) revealed CO oxidation coupled to H2 formation and the generation of a transmembrane electrochemical ion gradient ([Formula: see text]). Inverted membrane vesicles (IMVs) prepared from CO-grown cells also produced H2 and ATP from CO (via a loosely attached CO dehydrogenase) or a chemical reductant. Finally, we show that Ech activity led to the translocation of both H+ and Na+ across the membrane of the IMVs. The H+ gradient was then used by the ATP synthase for energy conservation. These data demonstrate that the energy-converting hydrogenase in concert with an ATP synthase may be the simplest form of respiration; it combines carbon dioxide fixation with the synthesis of ATP in an ancient pathway.

Fungi have the potential to produce a large repertoire of bioactive molecules, many of which can affect the growth and development of plants. Genomic survey of sequenced biofertilizer fungi showed many secondary metabolite gene clusters are anchored by iterative polyketide synthases (IPKSs), which are multidomain enzymes noted for generating diverse small molecules. Focusing on the biofertilizer Trichoderma harzianum t-22, we identified and characterized a cryptic IPKS-containing cluster that synthesizes tricholignan A, a redox-active ortho-hydroquinone. Tricholignan A is shown to reduce Fe(III) and may play a role in promoting plant growth under iron-deficient conditions. The construction of tricholignan by a pair of collaborating IPKSs was investigated using heterologous reconstitution and biochemical studies. A regioselective methylation step is shown to be a key step in formation of the ortho-hydroquinone. The responsible methyltransferase (MT) is fused with an N-terminal pseudo-acyl carrier protein (ψACP), in which the apo state of the ACP is essential for methylation of the growing polyketide chain. The ψACP is proposed to bind to the IPKS and enable the trans MT to access the growing polyketide. Our studies show that a genome-driven approach to discovering bioactive natural products from biofertilizer fungi can lead to unique compounds and biosynthetic knowledge.

Highbush blueberry (Vaccinium corymbosum) has long been consumed for its unique flavor and composition of health-promoting phytonutrients. However, breeding efforts to improve fruit quality in blueberry have been greatly hampered by the lack of adequate genomic resources and a limited understanding of the underlying genetics encoding key traits. The genome of highbush blueberry has been particularly challenging to assemble due, in large part, to its polyploid nature and genome size. Here, we present a chromosome-scale and haplotype-phased genome assembly of the cultivar "Draper," which has the highest antioxidant levels among a diversity panel of 71 cultivars and 13 wild Vaccinium species. We leveraged this genome, combined with gene expression and metabolite data measured across fruit development, to identify candidate genes involved in the biosynthesis of important phytonutrients among other metabolites associated with superior fruit quality. Genome-wide analyses revealed that both polyploidy and tandem gene duplications modified various pathways involved in the biosynthesis of key phytonutrients. Furthermore, gene expression analyses hint at the presence of a spatial-temporal specific dominantly expressed subgenome including during fruit development. These findings and the reference genome will serve as a valuable resource to guide future genome-enabled breeding of important agronomic traits in highbush blueberry.

Natural product discovery in the microbial world has historically been biased toward aerobes. Recent in silico analysis demonstrates that genomes of anaerobes encode unexpected biosynthetic potential for natural products, however, chemical data on natural products from the anaerobic world are extremely limited. Here, we review the current body of work on natural products isolated from strictly anaerobic microbes, including recent genome mining efforts to discover polyketides and non-ribosomal peptides from anaerobes. These known natural products of anaerobes have demonstrated interesting molecular scaffolds, biosynthetic logic, and/or biological activities, making anaerobes a promising reservoir for future natural product discovery.

Natural products are a large family of diverse and complex chemical molecules that have roles in both primary and secondary metabolism, and over 210,000 natural products have been described. Secondary metabolite natural products are of high commercial and societal value with therapeutic uses as antibiotics, antifungals, antitumor and antiparasitic products and in agriculture as products for crop protection and animal health. There is a resurgence of activity in exploring natural products for a wide range of applications, due to not only increasing antibiotic resistance, but the advent of next-generation genome sequencing and new technologies to interrogate and investigate natural product biosynthesis. Genome mining has revealed a previously undiscovered richness of biosynthetic potential in novel biosynthetic gene clusters for natural products. Complementing these computational processes are new experimental platforms that are being developed and deployed to access new natural products.

Disulfide-rich peptides (DRPs) play diverse physiological roles and have emerged as attractive sources of pharmacological tools and drug leads. Here we describe the 3D structure of a centipede venom peptide, U-SLPTX15-Sm2a, whose family defines a unique class of one of the most widespread DRP folds known, the cystine-stabilized α/β fold (CSαβ). This class, which we have named the two-disulfide CSαβ fold (2ds-CSαβ), contains only two internal disulfide bonds as opposed to at least three in all other confirmed CSαβ peptides, and constitutes one of the major neurotoxic peptide families in centipede venoms. We show the 2ds-CSαβ is widely distributed outside centipedes and is likely an ancient fold predating the split between prokaryotes and eukaryotes. Our results provide insights into the ancient evolutionary history of a widespread DRP fold and highlight the usefulness of 3D structures as evolutionary tools.

Testis-expressed X-linked genes typically evolve rapidly. Here, we report on a testis-expressed X-linked microRNA (miRNA) cluster that despite rapid alterations in sequence has retained its position in the Fragile-X region of the X chromosome in placental mammals. Surprisingly, the miRNAs encoded by this cluster (Fx-mir) have a predilection for targeting the immediately adjacent gene, Fmr1, an unexpected finding given that miRNAs usually act in trans, not in cis Robust repression of Fmr1 is conferred by combinations of Fx-mir miRNAs induced in Sertoli cells (SCs) during postnatal development when they terminate proliferation. Physiological significance is suggested by the finding that FMRP, the protein product of Fmr1, is downregulated when Fx-mir miRNAs are induced, and that FMRP loss causes SC hyperproliferation and spermatogenic defects. Fx-mir miRNAs not only regulate the expression of FMRP, but also regulate the expression of eIF4E and CYFIP1, which together with FMRP form a translational regulatory complex. Our results support a model in which Fx-mir family members act cooperatively to regulate the translation of batteries of mRNAs in a developmentally regulated manner in SCs.

Bacterial genes for molybdenum-containing and tungsten-containing enzymes are often differentially regulated depending on the metal availability in the environment. Here, we describe a new family of transcription factors with an unusual DNA-binding domain related to excisionases of bacteriophages. These transcription factors are associated with genes for various molybdate and tungstate-specific transporting systems as well as molybdo/tungsto-enzymes in a wide range of bacterial genomes. We used a combination of computational and experimental techniques to study a member of the TF family, named TaoR (for tungsten-containing aldehyde oxidoreductase regulator). In Desulfovibrio vulgaris Hildenborough, a model bacterium for sulfate reduction studies, TaoR activates expression of aldehyde oxidoreductase aor and represses tungsten-specific ABC-type transporter tupABC genes under tungsten-replete conditions. TaoR binding sites at aor promoter were identified by electrophoretic mobility shift assay and DNase I footprinting. We also reconstructed TaoR regulons in 45 Deltaproteobacteria by comparative genomics approach and predicted target genes for TaoR family members in other Proteobacteria and Firmicutes.

Pericyclases are an emerging family of enzymes catalyzing pericyclic reactions. A class of lipocalin-like enzymes recently characterized as Diels-Alderases (DAses) catalyze decalin formation through intramolecular Diels-Alder (IMDA) reactions between electron-rich dienes and electron-deficient dienophiles. Using this class of enzyme as a beacon for genome mining, we discovered a biosynthetic gene cluster from Penicillium variabile and identified that it encodes for the biosynthesis of varicidin A (1), a new antifungal natural product containing a cis-octahydrodecalin core. Biochemical analysis reveals a carboxylative deactivation strategy used in varicidin biosynthesis to suppress the nonenzymatic IMDA reaction of an early acyclic intermediate that favors trans-decalin formation. A P450 oxidizes the reactive intermediate to yield a relatively unreactive combination of an electron-deficient diene and an electron-deficient dienophile. The DAse PvhB catalyzes the final stage IMDA on the carboxylated intermediate to form the cis-decalin that is important for the antifungal activity.

Virulence mechanisms of the pathogenic fungus Aspergillus fumigatus are multifactorial and depend on the immune state of the host, but little is known about the fungal mechanism that develops during the process of lung invasion. In this study, microarray technology was combined with a histopathology evaluation of infected lungs so that the invasion strategy followed by the fungus could be described. To achieve this, an intranasal mice infection was performed to extract daily fungal samples from the infected lungs over four days post-infection. The pathological study revealed a heavy fungal progression throughout the lung, reaching the blood vessels on the third day after exposure and causing tissue necrosis. One percent of the fungal genome followed a differential expression pattern during this process. Strikingly, most of the genes of the intertwined fumagillin/pseurotin biosynthetic gene cluster were upregulated as were genes encoding lytic enzymes such as lipases, proteases (DppIV, DppV, Asp f 1 or Asp f 5) and chitinase (chiB1) as well as three genes related with pyomelanin biosynthesis process. Furthermore, we demonstrate that fumagillin is produced in an in vitro pneumocyte cell line infection model and that loss of fumagillin synthesis reduces epithelial cell damage. These results suggest that fumagillin contributes to tissue damage during invasive aspergillosis. Therefore, it is probable that A. fumigatus progresses through the lungs via the production of the mycotoxin fumagillin combined with the secretion of lytic enzymes that allow fungal growth, angioinvasion and the disruption of the lung parenchymal structure.

The increased interest in secondary metabolites (SMs) has driven a number of genome sequencing projects to elucidate their biosynthetic pathways. As a result, studies revealed that the number of secondary metabolite gene clusters (SMGCs) greatly outnumbers detected compounds, challenging current methods to dereplicate and categorize this amount of gene clusters on a larger scale. Here, we present an automated workflow for the genetic dereplication and analysis of secondary metabolism genes in fungi. Focusing on the secondary metabolite rich genus Aspergillus, we categorize SMGCs across genomes into SMGC families using network analysis. Our method elucidates the diversity and dynamics of secondary metabolism in section Nigri, showing that SMGC diversity within the section has the same magnitude as within the genus. Using our genome analysis we were able to predict the gene cluster responsible for biosynthesis of malformin, a potentiator of anti-cancer drugs, in 18 strains. To proof the general validity of our predictions, we developed genetic engineering tools in Aspergillus brasiliensis and subsequently verified the genes for biosynthesis of malformin.

BackgroundCopy Number Variations (CNVs) have becoming very significant variants, representing a major source of genomic variation. CNVs involvement in phenotypic expression and different diseases has been widely demonstrated in humans as well as in many domestic animals. However, genome wide investigation on these structural variations is still missing in Felis catus. The present work is the first CNV mapping from a large data set of Next Generation Sequencing (NGS) data in the domestic cat, performed within the 99 Lives Consortium.ResultsReads have been mapped on the reference assembly_6.2 by Maverix Biomics. CNV detection with cn.MOPS and CNVnator detected 592 CNVs. These CNVs were used to obtain 154 CNV Regions (CNVRs) with BedTools, including 62 singletons. CNVRs covered 0.26% of the total cat genome with 129 losses, 19 gains and 6 complexes. Cluster Analysis and Principal Component Analysis of the detected CNVRs showed that breeds tend to cluster together as well as cats sharing the same geographical origins. The 46 genes identified within the CNVRs were annotated.ConclusionThis study has improved the genomic characterization of 14 cat breeds and has provided CNVs information that can be used for studies of traits in cats. It can be considered a sound starting point for genomic CNVs identification in this species.

The goals of the KIR component of the 17th International HLA and Immunogenetics Workshop (IHIW) were to encourage and educate researchers to begin analyzing KIR at allelic resolution, and to survey the nature and extent of KIR allelic diversity across human populations. To represent worldwide diversity, we analyzed 1269 individuals from ten populations, focusing on the most polymorphic KIR genes, which express receptors having three immunoglobulin (Ig)-like domains (KIR3DL1/S1, KIR3DL2 and KIR3DL3). We identified 13 novel alleles of KIR3DL1/S1, 13 of KIR3DL2 and 18 of KIR3DL3. Previously identified alleles, corresponding to 33 alleles of KIR3DL1/S1, 38 of KIR3DL2, and 43 of KIR3DL3, represented over 90% of the observed allele frequencies for these genes. In total we observed 37 KIR3DL1/S1 allotypes, 40 for KIR3DL2 and 44 for KIR3DL3. As KIR allotype diversity can affect NK cell function, this demonstrates potential for high functional diversity worldwide. Allelic variation further diversifies KIR haplotypes. We determined KIR3DL3 ∼ KIR3DL1/S1 ∼ KIR3DL2 haplotypes from five of the studied populations, and observed multiple population-specific haplotypes in each. This included 234 distinct haplotypes in European Americans, 191 in Ugandans, 35 in Papuans, 95 in Egyptians and 86 in Spanish populations. For another 35 populations, encompassing 642,105 individuals we focused on KIR3DL2 and identified another 375 novel alleles, with approximately half of them observed in more than one individual. The KIR allelic level data gathered from this project represents the most comprehensive summary of global KIR allelic diversity to date, and continued analysis will improve understanding of KIR allelic polymorphism in global populations. Further, the wealth of new data gathered in the course of this workshop component highlights the value of collaborative, community-based efforts in immunogenetics research, exemplified by the IHIW.