Search

A human betaretrovirus (HBRV) has been linked with the autoimmune liver disease, primary biliary cholangitis (PBC), and various cancers, including breast cancer and lymphoma. HBRV is closely related to the mouse mammary tumor virus, and represents the only exogenous betaretrovirus characterized in humans to date. Evidence of infection in patients with PBC has been demonstrated through the identification of proviral integration sites in lymphoid tissue, the major reservoir of infection, as well as biliary epithelium, which is the site of the disease process. Accordingly, we tested the hypothesis that patients with PBC harbor a transmissible betaretrovirus by co-cultivation of PBC patients' lymph node homogenates with the HS578T breast cancer line. Because of the low level of HBRV replication, betaretrovirus producing cells were subcloned to optimize viral isolation and production. Evidence of infection was provided by electron microscopy, RT-PCR, in situ hybridization, cloning of the HBRV proviral genome and demonstration of more than 3400 integration sites. Further evidence of viral transmissibility was demonstrated by infection of biliary epithelial cells. While HBRV did not show a preference for integration proximal to specific genomic features, analyses of common insertion sites revealed evidence of integration proximal to cancer associated genes. These studies demonstrate the isolation of HBRV with features similar to mouse mammary tumor virus and confirm that patients with PBC display evidence of a transmissible viral infection.

Cation and anion channelrhodopsins (CCRs and ACRs, respectively) primarily from two algal species, Chlamydomonas reinhardtii and Guillardia theta, have become widely used as optogenetic tools to control cell membrane potential with light. We mined algal and other protist polynucleotide sequencing projects and metagenomic samples to identify 75 channelrhodopsin homologs from four channelrhodopsin families, including one revealed in dinoflagellates in this study. We carried out electrophysiological analysis of 33 natural channelrhodopsin variants from different phylogenetic lineages and 10 metagenomic homologs in search of sequence determinants of ion selectivity, photocurrent desensitization, and spectral tuning in channelrhodopsins. Our results show that association of a reduced number of glutamates near the conductance path with anion selectivity depends on a wider protein context, because prasinophyte homologs with a glutamate pattern identical to that in cryptophyte ACRs are cation selective. Desensitization is also broadly context dependent, as in one branch of stramenopile ACRs and their metagenomic homologs, its extent roughly correlates with phylogenetic relationship of their sequences. Regarding spectral tuning, we identified two prasinophyte CCRs with red-shifted spectra to 585 nm. They exhibit a third residue pattern in their retinal-binding pockets distinctly different from those of the only two types of red-shifted channelrhodopsins known (i.e., the CCR Chrimson and RubyACRs). In cryptophyte ACRs we identified three specific residue positions in the retinal-binding pocket that define the wavelength of their spectral maxima. Lastly, we found that dinoflagellate rhodopsins with a TCP motif in the third transmembrane helix and a metagenomic homolog exhibit channel activity. IMPORTANCE Channelrhodopsins are widely used in neuroscience and cardiology as research tools and are considered prospective therapeutics, but their natural diversity and mechanisms remain poorly characterized. Genomic and metagenomic sequencing projects are producing an ever-increasing wealth of data, whereas biophysical characterization of the encoded proteins lags behind. In this study, we used manual and automated patch clamp recording of representative members of four channelrhodopsin families, including a family in dinoflagellates that we report in this study. Our results contribute to a better understanding of molecular determinants of ionic selectivity, photocurrent desensitization, and spectral tuning in channelrhodopsins.

Premise: Large disjunctions in species distributions provide excellent opportunities to study processes that shape biogeographic patterns. One such disjunction is the eastern Asia-eastern North America (EA-ENA) floristic disjunction. For many genera with this disjunction, species richness is greater in EA than in ENA; this pattern has been attributed, in part, to higher rates of molecular evolution and speciation in EA. Longer branch lengths have been found in some EA clades, relative to their ENA sister clades, suggesting that the EA lineages have evolved at a higher rate, possibly due to environmental heterogeneity, potentially contributing to the species richness anomaly., Methods: To evaluate whether rates of molecular evolution are elevated in EA relative to ENA, we used transcriptomes from species in 11 genera displaying this disjunction. Rates of molecular evolution were estimated for up to 385 orthologous nuclear loci per genus., Results: No statistically significant differences were identified in pairwise comparisons between EA and ENA sister species, suggesting equal rates of molecular evolution for both species; the data also suggest similar selection pressures in both regions. For larger genera, evidence likewise argues against more species-rich clades having higher molecular evolutionary rates, regardless of region. Our results suggest that genes across multiple gene ontology categories are evolving at similar rates under purifying selection in species in both regions., Conclusions: Our data support the hypothesis that greater species richness in EA than ENA is due to factors other than an overall increase in rates of molecular evolution in EA., (© 2020 Botanical Society of America.)

Channelrhodopsins guide algal phototaxis and are widely used as optogenetic probes for control of membrane potential with light. "Bacteriorhodopsin-like" cation channelrhodopsins (BCCRs) from cryptophytes differ in primary structure from other CCRs, lacking usual residues important for their cation conductance. Instead, the sequences of BCCR match more closely those of rhodopsin proton pumps, containing residues responsible for critical proton transfer reactions. We report 19 new BCCRs which, together with the earlier 6 known members of this family, form three branches (subfamilies) of a phylogenetic tree. Here, we show that the conductance mechanisms in two subfamilies differ with respect to involvement of the homolog of the proton donor in rhodopsin pumps. Two BCCRs from the genus Rhodomonas generate photocurrents that rapidly desensitize under continuous illumination. Using a combination of patch clamp electrophysiology, absorption, Raman spectroscopy, and flash photolysis, we found that the desensitization is due to rapid accumulation of a long-lived nonconducting intermediate of the photocycle with unusually blue-shifted absorption with a maximum at 330 nm. These observations reveal diversity within the BCCR family and contribute to deeper understanding of their independently evolved cation channel function. IMPORTANCE Cation channelrhodopsins, light-gated channels from flagellate green algae, are extensively used as optogenetic photoactivators of neurons in research and recently have progressed to clinical trials for vision restoration. However, the molecular mechanisms of their photoactivation remain poorly understood. We recently identified cryptophyte cation channelrhodopsins, structurally different from those of green algae, which have separately evolved to converge on light-gated cation conductance. This study reveals diversity within this new protein family and describes a subclade with unusually rapid desensitization that results in short transient photocurrents in continuous light. Such transient currents have not been observed in the green algae channelrhodopsins and are potentially useful in optogenetic protocols. Kinetic UV-visible (UV-vis) spectroscopy and photoelectrophysiology reveal that the desensitization is caused by rapid accumulation of a nonconductive photointermediate in the photochemical reaction cycle. The absorption maximum of the intermediate is 330 nm, the shortest wavelength reported in any rhodopsin, indicating a novel chromophore structure., (Copyright © 2020 Sineshchekov et al.)

Hornworts comprise a bryophyte lineage that diverged from other extant land plants >400 million years ago and bears unique biological features, including a distinct sporophyte architecture, cyanobacterial symbiosis and a pyrenoid-based carbon-concentrating mechanism (CCM). Here, we provide three high-quality genomes of Anthoceros hornworts. Phylogenomic analyses place hornworts as a sister clade to liverworts plus mosses with high support. The Anthoceros genomes lack repeat-dense centromeres as well as whole-genome duplication, and contain a limited transcription factor repertoire. Several genes involved in angiosperm meristem and stomatal function are conserved in Anthoceros and upregulated during sporophyte development, suggesting possible homologies at the genetic level. We identified candidate genes involved in cyanobacterial symbiosis and found that LCIB, a Chlamydomonas CCM gene, is present in hornworts but absent in other plant lineages, implying a possible conserved role in CCM function. We anticipate that these hornwort genomes will serve as essential references for future hornwort research and comparative studies across land plants.

Plants are the foundation of terrestrial ecosystems, and their colonization of land was probably facilitated by mutualistic associations with arbuscular mycorrhizal fungi. Following this founding event, plant diversification has led to the emergence of a tremendous diversity of mutualistic symbioses with microorganisms, ranging from extracellular associations to the most intimate intracellular associations, where fungal or bacterial symbionts are hosted inside plant cells. Here, through analysis of 271 transcriptomes and 116 plant genomes spanning the entire land-plant diversity, we demonstrate that a common symbiosis signalling pathway co-evolved with intracellular endosymbioses, from the ancestral arbuscular mycorrhiza to the more recent ericoid and orchid mycorrhizae in angiosperms and ericoid-like associations of bryophytes. By contrast, species forming exclusively extracellular symbioses, such as ectomycorrhizae, and those forming associations with cyanobacteria, have lost this signalling pathway. This work unifies intracellular symbioses, revealing conservation in their evolution across 450 million yr of plant diversification.

Premise: Phylogenetic trees of bryophytes provide important evolutionary context for land plants. However, published inferences of overall embryophyte relationships vary considerably. We performed phylogenomic analyses of bryophytes and relatives using both mitochondrial and plastid gene sets, and investigated bryophyte plastome evolution., Methods: We employed diverse likelihood-based analyses to infer large-scale bryophyte phylogeny for mitochondrial and plastid data sets. We tested for changes in purifying selection in plastid genes of a mycoheterotrophic liverwort (Aneura mirabilis) and a putatively mycoheterotrophic moss (Buxbaumia), and compared 15 bryophyte plastomes for major structural rearrangements., Results: Overall land-plant relationships conflict across analyses, generally weakly. However, an underlying (unrooted) four-taxon tree is consistent across most analyses and published studies. Despite gene coverage patchiness, relationships within mosses, liverworts, and hornworts are largely congruent with previous studies, with plastid results generally better supported. Exclusion of RNA edit sites restores cases of unexpected non-monophyly to monophyly for Takakia and two hornwort genera. Relaxed purifying selection affects multiple plastid genes in mycoheterotrophic Aneura but not Buxbaumia. Plastid genome structure is nearly invariant across bryophytes, but the tufA locus, presumed lost in embryophytes, is unexpectedly retained in several mosses., Conclusions: A common unrooted tree underlies embryophyte phylogeny, [(liverworts, mosses), (hornworts, vascular plants)]; rooting inconsistency across studies likely reflects substantial distance to algal outgroups. Analyses combining genomic and transcriptomic data may be misled locally for heavily RNA-edited taxa. The Buxbaumia plastome lacks hallmarks of relaxed selection found in mycoheterotrophic Aneura. Autotrophic bryophyte plastomes, including Buxbaumia, hardly vary in overall structure., (© 2019 Botanical Society of America.)

Chloroplast retrograde signaling networks are vital for chloroplast biogenesis, operation, and signaling, including excess light and drought stress signaling. To date, retrograde signaling has been considered in the context of land plant adaptation, but not regarding the origin and evolution of signaling cascades linking chloroplast function to stomatal regulation. We show that key elements of the chloroplast retrograde signaling process, the nucleotide phosphatase (SAL1) and 3'-phosphoadenosine-5'-phosphate (PAP) metabolism, evolved in streptophyte algae-the algal ancestors of land plants. We discover an early evolution of SAL1-PAP chloroplast retrograde signaling in stomatal regulation based on conserved gene and protein structure, function, and enzyme activity and transit peptides of SAL1s in species including flowering plants, the fern Ceratopteris richardii , and the moss Physcomitrella patens Moreover, we demonstrate that PAP regulates stomatal closure via secondary messengers and ion transport in guard cells of these diverse lineages. The origin of stomata facilitated gas exchange in the earliest land plants. Our findings suggest that the conquest of land by plants was enabled by rapid response to drought stress through the deployment of an ancestral SAL1-PAP signaling pathway, intersecting with the core abscisic acid signaling in stomatal guard cells., Competing Interests: The authors declare no conflict of interest., (Copyright © 2019 the Author(s). Published by PNAS.)

Premise of the Study: Studies of gene expression and polyploidy are typically restricted to characterizing differences in transcript concentration. Using diploid and autotetraploid Tolmiea, we present an integrated approach for cross-ploidy comparisons that account for differences in transcriptome size and cell density and make multiple comparisons of transcript abundance., Methods: We use RNA spike-in standards in concert with cell size and density to identify and correct for differences in transcriptome size and compare levels of gene expression across multiple scales: per transcriptome, per cell, and per biomass., Key Results: In total, ~17% of all loci were identified as differentially expressed (DEGs) between the diploid and autopolyploid species. The per-transcriptome normalization, the method researchers typically use, captured the fewest DEGs (58% of total DEGs) and failed to detect any DEGs not found by the alternative normalizations. When transcript abundance was normalized per biomass and per cell, ~66% and ~82% of the total DEGs were recovered, respectively. The discrepancy between per-transcriptome and per-cell recovery of DEGs occurs because per-transcriptome normalizations are concentration-based and therefore blind to differences in transcriptome size., Conclusions: While each normalization enables valid comparisons at biologically relevant scales, a holistic comparison of multiple normalizations provides additional explanatory power not available from any single approach. Notably, autotetraploid loci tend to conserve diploid-like transcript abundance per biomass through increased gene expression per cell, and these loci are enriched for photosynthesis-related functions., (© 2019 Botanical Society of America.)

Background: Inflammatory bowel diseases (IBD) are a group of complex and multifactorial disorders with unknown etiology. Chronic intestinal inflammation develops against resident intestinal bacteria in genetically susceptible hosts. We hypothesized that host intestinal immunoglobulin (Ig) G can be used to identify bacteria involved in IBD pathogenesis., Results: IgG-bound and -unbound microorganisms were collected from 32 pediatric terminal ileum aspirate washes during colonoscopy [non-IBD (n = 10), Crohn disease (n = 15), and ulcerative colitis (n = 7)], and composition was assessed using the Illumina MiSeq platform. In vitro analysis of invasive capacity was evaluated by fluorescence in situ hybridization and gentamicin invasion assay; immune activation was measured by qPCR. Despite considerable inter-individual variations, IgG binding favored specific and unique mucosa-associated species in pediatric IBD patients. Burkholderia cepacia, Flavonifractor plautii, and Rumminococcus sp. demonstrated increased IgG binding, while Pseudomonas ST29 demonstrated reduced IgG binding, in IBD. In vitro validation confirmed that B. cepacia, F. plautii, and Rumminococcus display invasive potential while Pseudomonas protogens did not., Conclusion: Using IgG as a marker of pathobionts in larger patient cohorts to identify microbes and elucidate their role in IBD pathogenesis will potentially underpin new strategies to facilitate development of novel, targeted diagnostic, and therapeutic approaches. Interestingly, this method can be used beyond the scope of this manuscript to evaluate altered gut pathobionts in a number of diseases associated with altered microbiota including arthritis, obesity, diabetes mellitus, alcoholic liver disease, cirrhosis, metabolic syndrome, and carcinomas.

Optogenetic inhibition of specific neuronal types in the brain enables analysis of neural circuitry and is promising for the treatment of a number of neurological disorders. Anion channelrhodopsins (ACRs) from the cryptophyte alga Guillardia theta generate larger photocurrents than other available inhibitory optogenetic tools, but more rapid channels are needed for temporally precise inhibition, such as single-spike suppression, of high-frequency firing neurons. Faster ACRs have been reported, but their potential advantages for time-resolved inhibitory optogenetics have not so far been verified in neurons. We report RapACR, nicknamed so for "rapid," an ACR from Rhodomonas salina , that exhibits channel half-closing times below 10 ms and achieves equivalent inhibition at 50-fold lower light intensity in lentivirally transduced cultured mouse hippocampal neurons as the second-generation engineered Cl - -conducting channelrhodopsin iC++. The upper limit of the time resolution of neuronal silencing with RapACR determined by measuring the dependence of spiking recovery after photoinhibition on the light intensity was calculated to be 100 Hz, whereas that with the faster of the two G. theta ACRs was 13 Hz. Further acceleration of RapACR channel kinetics was achieved by site-directed mutagenesis of a single residue in transmembrane helix 3 (Thr111 to Cys). We also show that mutation of another ACR (Cys to Ala at the same position) with a greatly extended lifetime of the channel open state acts as a bistable photochromic tool in mammalian neurons. These molecules extend the time domain of optogenetic neuronal silencing while retaining the high light sensitivity of Guillardia ACRs.

Ferns are the closest sister group to all seed plants, yet little is known about their genomes other than that they are generally colossal. Here, we report on the genomes of Azolla filiculoides and Salvinia cucullata (Salviniales) and present evidence for episodic whole-genome duplication in ferns-one at the base of 'core leptosporangiates' and one specific to Azolla. One fern-specific gene that we identified, recently shown to confer high insect resistance, seems to have been derived from bacteria through horizontal gene transfer. Azolla coexists in a unique symbiosis with N 2 -fixing cyanobacteria, and we demonstrate a clear pattern of cospeciation between the two partners. Furthermore, the Azolla genome lacks genes that are common to arbuscular mycorrhizal and root nodule symbioses, and we identify several putative transporter genes specific to Azolla-cyanobacterial symbiosis. These genomic resources will help in exploring the biotechnological potential of Azolla and address fundamental questions in the evolution of plant life.

Premise of the Study: For the past one billion years, green plants (Viridiplantae) have dominated global ecosystems, yet many key branches in their evolutionary history remain poorly resolved. Using the largest analysis of Viridiplantae based on plastid genome sequences to date, we examined the phylogeny and implications for morphological evolution at key nodes., Methods: We analyzed amino acid sequences from protein-coding genes from complete (or nearly complete) plastomes for 1879 taxa, including representatives across all major clades of Viridiplantae. Much of the data used was derived from transcriptomes from the One Thousand Plants Project (1KP); other data were taken from GenBank., Key Results: Our results largely agree with previous plastid-based analyses. Noteworthy results include (1) the position of Zygnematophyceae as sister to land plants (Embryophyta), (2) a bryophyte clade (hornworts, mosses + liverworts), (3) Equisetum + Psilotaceae as sister to Marattiales + leptosporangiate ferns, (4) cycads + Ginkgo as sister to the remaining extant gymnosperms, within which Gnetophyta are placed within conifers as sister to non-Pinaceae (Gne-Cup hypothesis), and (5) Amborella, followed by water lilies (Nymphaeales), as successive sisters to all other extant angiosperms. Within angiosperms, there is support for Mesangiospermae, a clade that comprises magnoliids, Chloranthales, monocots, Ceratophyllum, and eudicots. The placements of Ceratophyllum and Dilleniaceae remain problematic. Within Pentapetalae, two major clades (superasterids and superrosids) are recovered., Conclusions: This plastid data set provides an important resource for elucidating morphological evolution, dating divergence times in Viridiplantae, comparisons with emerging nuclear phylogenies, and analyses of molecular evolutionary patterns and dynamics of the plastid genome., (© 2018 Botanical Society of America.)

Tubulin is the target for many small-molecule natural compounds, which alter microtubules dynamics, and lead to cell cycle arrest and apoptosis. One of these compounds is colchicine, a plant alkaloid produced by Colchicum autumnale . While C. autumnale produces a potent cytotoxin, colchicine, and expresses its target protein, it is immune to colchicine's cytotoxic action and the mechanism of this resistance is hitherto unknown. In the present paper, the molecular mechanisms responsible for colchicine resistance in C. autumnale are investigated and compared to human tubulin. To this end, homology models for C. autumnale α-β tubulin heterodimer are created and molecular dynamics (MD) simulations together with molecular mechanics Poisson-Boltzmann calculations (MM/PBSA) are performed to determine colchicine's binding affinity for tubulin. Using our molecular approach, it is shown that the colchicine-binding site in C. autumnale tubulin contains a small number of amino acid substitutions compared to human tubulin. However, these substitutions induce significant reduction in the binding affinity for tubulin, and subsequently fewer conformational changes in its structure result. It is suggested that such small conformational changes are insufficient to profoundly disrupt microtubule dynamics, which explains the high resistance to colchicine by C. autumnale ., Competing Interests: The authors declare no conflict of interest.

Seeds are one of the most significant innovations in the land plant lineage, critical to the diversification and adaptation of plants to terrestrial environments. From perspective of seed evo-devo, the most crucial developmental stage in this innovation is seed maturation, which includes accumulation of storage reserves, acquisition of desiccation tolerance, and induction of dormancy. Based on previous studies of seed development in the model plant Arabidopsis thaliana , seed maturation is mainly controlled by the LAFL regulatory network, which includes LEAFY COTYLEDON1 ( LEC1 ) and LEC1-LIKE ( L1L ) of the NF-YB gene family, and ABSCISIC ACID INSENSITIVE3 ( ABI3 ), FUSCA3 ( FUS3 ), and LEC2 ( LEAFY COTYLEDON2 ) of the B3- AFL gene family. In the present study, molecular evolution of these LAFL genes was analyzed, using representative species from across the major plant lineages. Additionally, to elucidate the molecular mechanisms of the seed maturation program, co-expression pattern analyses of LAFL genes were conducted across vascular plants. The results show that the origin of AFL gene family dates back to a common ancestor of bryophytes and vascular plants, while LEC1 -type genes are only found in vascular plants. LAFL genes of vascular plants likely specify their co-expression in two different developmental phrases, spore and seed maturation, respectively, and expression patterns vary slightly across the major vascular plants lineages. All the information presented in this study will provide insights into the origin and diversification of seed plants.

Natural anion channelrhodopsins (ACRs) discovered in the cryptophyte alga Guillardia theta generate large hyperpolarizing currents at membrane potentials above the Nernst equilibrium potential for Cl - and thus can be used as efficient inhibitory tools for optogenetics. We have identified and characterized new ACR homologs in different cryptophyte species, showing that all of them are anion-selective, and thus expanded this protein family to 20 functionally confirmed members. Sequence comparison of natural ACRs and engineered Cl - -conducting mutants of cation channelrhodopsins (CCRs) showed radical differences in their anion selectivity filters. In particular, the Glu90 residue in channelrhodopsin 2, which needed to be mutated to a neutral or alkaline residue to confer anion selectivity to CCRs, is nevertheless conserved in all of the ACRs identified. The new ACRs showed a large variation of the amplitude, kinetics, and spectral sensitivity of their photocurrents. A notable variant, designated "ZipACR", is particularly promising for inhibitory optogenetics because of its combination of larger current amplitudes than those of previously reported ACRs and an unprecedentedly fast conductance cycle (current half-decay time 2-4 ms depending on voltage). ZipACR expressed in cultured mouse hippocampal neurons enabled precise photoinhibition of individual spikes in trains of up to 50 Hz frequency.

The anatomically simple plants that first colonized land must have acquired molecular and biochemical adaptations to drought stress. Abscisic acid (ABA) coordinates responses leading to desiccation tolerance in all land plants. We identified ABA nonresponsive mutants in the model bryophyte Physcomitrella patens and genotyped a segregating population to map and identify the ABA NON-RESPONSIVE (ANR) gene encoding a modular protein kinase comprising an N-terminal PAS domain, a central EDR domain, and a C-terminal MAPKKK-like domain. anr mutants fail to accumulate dehydration tolerance-associated gene products in response to drought, ABA, or osmotic stress and do not acquire ABA-dependent desiccation tolerance. The crystal structure of the PAS domain, determined to 1.7-Å resolution, shows a conserved PAS-fold that dimerizes through a weak dimerization interface. Targeted mutagenesis of a conserved tryptophan residue within the PAS domain generates plants with ABA nonresponsive growth and strongly attenuated ABA-responsive gene expression, whereas deleting this domain retains a fully ABA-responsive phenotype. ANR orthologs are found in early-diverging land plant lineages and aquatic algae but are absent from more recently diverged vascular plants. We propose that ANR genes represent an ancestral adaptation that enabled drought stress survival of the first terrestrial colonizers but were lost during land plant evolution., (© 2016 American Society of Plant Biologists. All rights reserved.)

The advent of next generation sequencing (NGS) has enabled investigations of the gut microbiome with unprecedented resolution and throughput. This has stimulated the development of sophisticated bioinformatics tools to analyze the massive amounts of data generated. Researchers therefore need a clear understanding of the key concepts required for the design, execution and interpretation of NGS experiments on microbiomes. We conducted a literature review and used our own data to determine which approaches work best. The two main approaches for analyzing the microbiome, 16S ribosomal RNA (rRNA) gene amplicons and shotgun metagenomics, are illustrated with analyses of libraries designed to highlight their strengths and weaknesses. Several methods for taxonomic classification of bacterial sequences are discussed. We present simulations to assess the number of sequences that are required to perform reliable appraisals of bacterial community structure. To the extent that fluctuations in the diversity of gut bacterial populations correlate with health and disease, we emphasize various techniques for the analysis of bacterial communities within samples (α-diversity) and between samples (β-diversity). Finally, we demonstrate techniques to infer the metabolic capabilities of a bacteria community from these 16S and shotgun data.

Retroviruses have been foundational in cancer research since early studies identified proto-oncogenes as targets for insertional mutagenesis. Integration of murine gamma-retroviruses into the host genome favours promoters and enhancers and entails interaction of viral integrase with host BET/bromodomain factors. We report that this integration pattern is conserved in feline leukaemia virus (FeLV), a gamma-retrovirus that infects many human cell types. Analysis of FeLV insertion sites in the MCF-7 mammary carcinoma cell line revealed strong bias towards active chromatin marks with no evidence of significant post-integration growth selection. The most prominent FeLV integration targets had little overlap with the most abundantly expressed transcripts, but were strongly enriched for annotated cancer genes. A meta-analysis based on several gamma-retrovirus integration profiling (GRIP) studies in human cells (CD34+, K562, HepG2) revealed a similar cancer gene bias but also remarkable cell-type specificity, with prominent exceptions including a universal integration hotspot at the long non-coding RNA MALAT1. Comparison of GRIP targets with databases of super-enhancers from the same cell lines showed that these have only limited overlap and that GRIP provides unique insights into the upstream drivers of cell growth. These observations elucidate the oncogenic potency of the gamma-retroviruses and support the wider application of GRIP to identify the genes and growth regulatory circuits that drive distinct cancer types.

Betalain pigments are unique to the Caryophyllales and structurally and biosynthetically distinct from anthocyanins. Two key enzymes within the betalain synthesis pathway have been identified: 4,5-dioxygenase (DODA) that catalyzes the formation of betalamic acid and CYP76AD1, a cytochrome P450 gene that catalyzes the formation of cyclo-DOPA. We performed phylogenetic analyses to reveal the evolutionary history of the DODA and CYP76AD1 lineages and in the context of an ancestral reconstruction of pigment states we explored the evolution of these genes in relation to the complex evolution of pigments in Caryophylalles. Duplications within the CYP76AD1 and DODA lineages arose just before the origin of betalain pigmentation in the core Caryophyllales. The duplications gave rise to DODA-α and CYP76AD1-α isoforms that appear specific to betalain synthesis. Both betalain-specific isoforms were then lost or downregulated in the anthocyanic Molluginaceae and Caryophyllaceae. Our findings suggest a single origin of the betalain synthesis pathway, with neofunctionalization following gene duplications in the CYP76AD1 and DODA lineages. Loss of DODA-α and CYP76AD1-α in anthocyanic taxa suggests that betalain pigmentation has been lost twice in Caryophyllales, and exclusion of betalain pigments from anthocyanic taxa is mediated through gene loss or downregulation. [Correction added after online publication 13 May 2015: in the last two paragraphs of the Summary the gene name CYP761A was changed to CYP76AD1.]., (© 2015 The Authors. New Phytologist © 2015 New Phytologist Trust.)

Plant phototropism, the ability to bend toward or away from light, is predominantly controlled by blue-light photoreceptors, the phototropins. Although phototropins have been well-characterized in Arabidopsis thaliana, their evolutionary history is largely unknown. In this study, we complete an in-depth survey of phototropin homologs across land plants and algae using newly available transcriptomic and genomic data. We show that phototropins originated in an ancestor of Viridiplantae (land plants + green algae). Phototropins repeatedly underwent independent duplications in most major land-plant lineages (mosses, lycophytes, ferns, and seed plants), but remained single-copy genes in liverworts and hornworts-an evolutionary pattern shared with another family of photoreceptors, the phytochromes. Following each major duplication event, the phototropins differentiated in parallel, resulting in two specialized, yet partially overlapping, functional forms that primarily mediate either low- or high-light responses. Our detailed phylogeny enables us to not only uncover new phototropin lineages, but also link our understanding of phototropin function in Arabidopsis with what is known in Adiantum and Physcomitrella (the major model organisms outside of flowering plants). We propose that the convergent functional divergences of phototropin paralogs likely contributed to the success of plants through time in adapting to habitats with diverse and heterogeneous light conditions.

Unlabelled: •, Premise of the Study: Understanding fern (monilophyte) phylogeny and its evolutionary timescale is critical for broad investigations of the evolution of land plants, and for providing the point of comparison necessary for studying the evolution of the fern sister group, seed plants. Molecular phylogenetic investigations have revolutionized our understanding of fern phylogeny, however, to date, these studies have relied almost exclusively on plastid data.•, Methods: Here we take a curated phylogenomics approach to infer the first broad fern phylogeny from multiple nuclear loci, by combining broad taxon sampling (73 ferns and 12 outgroup species) with focused character sampling (25 loci comprising 35877 bp), along with rigorous alignment, orthology inference and model selection.•, Key Results: Our phylogeny corroborates some earlier inferences and provides novel insights; in particular, we find strong support for Equisetales as sister to the rest of ferns, Marattiales as sister to leptosporangiate ferns, and Dennstaedtiaceae as sister to the eupolypods. Our divergence-time analyses reveal that divergences among the extant fern orders all occurred prior to ∼200 MYA. Finally, our species-tree inferences are congruent with analyses of concatenated data, but generally with lower support. Those cases where species-tree support values are higher than expected involve relationships that have been supported by smaller plastid datasets, suggesting that deep coalescence may be reducing support from the concatenated nuclear data.•, Conclusions: Our study demonstrates the utility of a curated phylogenomics approach to inferring fern phylogeny, and highlights the need to consider underlying data characteristics, along with data quantity, in phylogenetic studies., (© 2015 Botanical Society of America, Inc.)

The extracellular matrix of scaly green flagellates consists of small organic scales consisting of polysaccharides and scale-associated proteins (SAPs). Molecular phylogenies have shown that these organisms represent the ancestral stock of flagellates from which all green plants (Viridiplantae) evolved. The molecular characterization of four different SAPs is presented. Three SAPs are type-2 membrane proteins with an arginine/alanine-rich short cytoplasmic tail and an extracellular domain that is most likely of bacterial origin. The fourth protein is a filamin-like protein. In addition, we report the presence of proteins similar to the integrin-associated proteins α-actinin (in transcriptomes of glaucophytes and some viridiplants), LIM-domain proteins, and integrin-associated kinase in transcriptomes of viridiplants, glaucophytes, and rhodophytes. We propose that the membrane proteins identified are the predicted linkers between scales and the cytoskeleton. These proteins are present in many green algae but are apparently absent from embryophytes. These proteins represent a new protein family we have termed gralins for green algal integrins. Gralins are absent from embryophytes. A model for the evolution of the cell surface proteins in Plantae is discussed., (© The Author(s) 2015. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Background and Aims: Cultivated flax (Linum usitatissimum) is known to have undergone a whole-genome duplication around 5-9 million years ago. The aim of this study was to investigate whether other whole-genome duplication events have occurred in the evolutionary history of cultivated flax. Knowledge of such whole-genome duplications will be important in understanding the biology and genomics of cultivated flax., Methods: Transcriptomes of 11 Linum species were sequenced using the Illumina platform. The short reads were assembled de novo and the DupPipe pipeline was used to look for signatures of polyploidy events from the age distribution of paralogues. In addition, phylogenies of all paralogues were assembled within an estimated age window of interest. These phylogenies were assessed for evidence of a paleopolyploidy event within the genus Linum., Key Results: A previously unknown paleopolyploidy event that occurred 20-40 million years ago was discovered and shown to be specific to a clade within Linum containing cultivated flax (L. usitatissimum) and other mainly blue-flowered species. The finding was supported by two lines of evidence. First, a significant change of slope (peak) was shown in the age distribution of paralogues that was phylogenetically restricted to, and ubiquitous in, this clade. Second, a large number of paralogue phylogenies were retrieved that are consistent with a polyploidy event occurring within that clade., Conclusions: The results show the utility of multi-species transcriptomics for detecting whole-genome duplication events and demonstrate that that multiple rounds of polyploidy have been important in shaping the evolutionary history of flax. Understanding and characterizing these whole-genome duplication events will be important for future Linum research.

Flax (Linum usitatissimum) is an ancient crop that is widely cultivated as a source of fiber, oil and medicinally relevant compounds. To accelerate crop improvement, we performed whole-genome shotgun sequencing of the nuclear genome of flax. Seven paired-end libraries ranging in size from 300 bp to 10 kb were sequenced using an Illumina genome analyzer. A de novo assembly, comprised exclusively of deep-coverage (approximately 94× raw, approximately 69× filtered) short-sequence reads (44-100 bp), produced a set of scaffolds with N(50) =694 kb, including contigs with N(50)=20.1 kb. The contig assembly contained 302 Mb of non-redundant sequence representing an estimated 81% genome coverage. Up to 96% of published flax ESTs aligned to the whole-genome shotgun scaffolds. However, comparisons with independently sequenced BACs and fosmids showed some mis-assembly of regions at the genome scale. A total of 43384 protein-coding genes were predicted in the whole-genome shotgun assembly, and up to 93% of published flax ESTs, and 86% of A. thaliana genes aligned to these predicted genes, indicating excellent coverage and accuracy at the gene level. Analysis of the synonymous substitution rates (K(s) ) observed within duplicate gene pairs was consistent with a recent (5-9 MYA) whole-genome duplication in flax. Within the predicted proteome, we observed enrichment of many conserved domains (Pfam-A) that may contribute to the unique properties of this crop, including agglutinin proteins. Together these results show that de novo assembly, based solely on whole-genome shotgun short-sequence reads, is an efficient means of obtaining nearly complete genome sequence information for some plant species., (© 2012 The Authors. The Plant Journal © 2012 Blackwell Publishing Ltd.)

Large-insert genome analysis (LIGAN) is a broadly applicable, high-throughput technology designed to characterize genome-scale structural variation. Fosmid paired-end sequences and DNA fingerprints from a query genome are compared to a reference sequence using the Genomic Variation Analysis (GenVal) suite of software tools to pinpoint locations of insertions, deletions, and rearrangements. Fosmids spanning regions that contain new structural variants can then be sequenced. Clonal pairs of Pseudomonas aeruginosa isolates from four cystic fibrosis patients were used to validate the LIGAN technology. Approximately 1.5 Mb of inserted sequences were identified, including 743 kb containing 615 ORFs that are absent from published P. aeruginosa genomes. Six rearrangement breakpoints and 220 kb of deleted sequences were also identified. Our study expands the "genome universe" of P. aeruginosa and validates a technology that complements emerging, short-read sequencing methods that are better suited to characterizing single-nucleotide polymorphisms than structural variation.

After the completion of a draft human genome sequence, the International Human Genome Sequencing Consortium has proceeded to finish and annotate each of the 24 chromosomes comprising the human genome. Here we describe the sequencing and analysis of human chromosome 3, one of the largest human chromosomes. Chromosome 3 comprises just four contigs, one of which currently represents the longest unbroken stretch of finished DNA sequence known so far. The chromosome is remarkable in having the lowest rate of segmental duplication in the genome. It also includes a chemokine receptor gene cluster as well as numerous loci involved in multiple human cancers such as the gene encoding FHIT, which contains the most common constitutive fragile site in the genome, FRA3B. Using genomic sequence from chimpanzee and rhesus macaque, we were able to characterize the breakpoints defining a large pericentric inversion that occurred some time after the split of Homininae from Ponginae, and propose an evolutionary history of the inversion.

Background: Comparative whole genome analysis of Mammalia can benefit from the addition of more species. The pig is an obvious choice due to its economic and medical importance as well as its evolutionary position in the artiodactyls., Results: We have generated approximately 3.84 million shotgun sequences (0.66X coverage) from the pig genome. The data are hereby released (NCBI Trace repository with center name "SDJVP", and project name "Sino-Danish Pig Genome Project") together with an initial evolutionary analysis. The non-repetitive fraction of the sequences was aligned to the UCSC human-mouse alignment and the resulting three-species alignments were annotated using the human genome annotation. Ultra-conserved elements and miRNAs were identified. The results show that for each of these types of orthologous data, pig is much closer to human than mouse is. Purifying selection has been more efficient in pig compared to human, but not as efficient as in mouse, and pig seems to have an isochore structure most similar to the structure in human., Conclusion: The addition of the pig to the set of species sequenced at low coverage adds to the understanding of selective pressures that have acted on the human genome by bisecting the evolutionary branch between human and mouse with the mouse branch being approximately 3 times as long as the human branch. Additionally, the joint alignment of the shot-gun sequences to the human-mouse alignment offers the investigator a rapid way to defining specific regions for analysis and resequencing.

The genome sequence of the Severe Acute Respiratory Syndrome (SARS)-associated virus provides essential information for the identification of pathogen(s), exploration of etiology and evolution, interpretation of transmission and pathogenesis, development of diagnostics, prevention by future vaccination, and treatment by developing new drugs. We report the complete genome sequence and comparative analysis of an isolate (BJ01) of the coronavirus that has been recognized as a pathogen for SARS. The genome is 29725 nt in size and has 11 ORFs (Open Reading Frames). It is composed of a stable region encoding an RNA-dependent RNA polymerase (composed of 2 ORFs) and a variable region representing 4 CDSs (coding sequences) for viral structural genes (the S, E, M, N proteins) and 5 PUPs (putative uncharacterized proteins). Its gene order is identical to that of other known coronaviruses. The sequence alignment with all known RNA viruses places this virus as a member in the family of Coronaviridae. Thirty putative substitutions have been identified by comparative analysis of the 5 SARS-associated virus genome sequences in GenBank. Fifteen of them lead to possible amino acid changes (non-synonymous mutations) in the proteins. Three amino acid changes, with predicted alteration of physical and chemical features, have been detected in the S protein that is postulated to be involved in the immunoreactions between the virus and its host. Two amino acid changes have been detected in the M protein, which could be related to viral envelope formation. Phylogenetic analysis suggests the possibility of non-human origin of the SARS-associated viruses but provides no evidence that they are man-made. Further efforts should focus on identifying the etiology of the SARS-associated virus and ruling out conclusively the existence of other possible SARS-related pathogen(s)., (© Science in China Press 2003.)