Search

Frailty indexes (FIs) provide quantitative measurements of nonspecific health decline and are particularly useful as longitudinal monitors of morbidity in aging studies. For mouse studies, frailty assessments can be taken noninvasively, but they require handling and direct observation that is labor-intensive to the scientist and stress inducing to the animal. Here, we implement, evaluate, and provide a refined digital FI composed entirely of computational analyses of home-cage video and compare it to manually obtained frailty scores in both C57BL/6 and genetically heterogeneous Diversity Outbred mice. We show that the frailty scores assigned by our digital index correlate with both manually obtained frailty scores and chronological age. Thus, we provide an automated tool for frailty assessment that can be collected reproducibly, at scale, without substantial labor cost.

Diffuse idiopathic skeletal hyperostosis (DISH) is a condition where adjacent vertebrae become fused through formation of osteophytes. The genetic and epidemiological etiology of this condition is not well understood. Here, we implemented a machine learning algorithm to assess the prevalence and severity of the pathology in ~40,000 lateral DXA scans in the UK Biobank Imaging cohort. We find that DISH is highly prevalent, above the age of 45, ~20% of men and ~8% of women having multiple osteophytes. Surprisingly, we find strong phenotypic and genetic association of DISH with increased bone mineral density and content throughout the entire skeletal system. Genetic association analysis identified ten loci associated with DISH, including multiple genes involved in bone remodeling (RUNX2, IL11, GDF5, CCDC91, NOG, and ROR2). Overall, this study describes genetics of DISH and implicates the role of overactive osteogenesis as a key driver of the pathology.

The naked mole-rat (Heterocephalus glaber) is a mouse-sized rodent species, notable for its eusociality and long lifespan. Previously, we reported that demographic aging, i.e., the exponential increase of mortality hazard that accompanies advancing age in mammals and other organisms, does not occur in naked mole-rats (Ruby et al, 2018). The demographic data supporting that conclusion had taken over three decades to accumulate, starting with the original rearing ofH.glaberin captivity. In the five years following that study, we ∼doubled our quantity of demographic data. Here, we re-evaluated our prior conclusions in light of these new data and found them to be supported and indeed strengthened. We additionally provided insight into the social dynamics of captiveH.glaberwith data and analyses of body weight and colony size versus mortality. Finally, we provide a phylogenetically-proximal comparator in the form of lifespan data from our Damaraland mole-rat (Fukomys damarensis) colony and demographic meta-analysis of those data along with published data from Ansell’s mole-rat (Fukomys anselli). We foundFukomysmortality hazard to increase gradually with age, an observation with implications on the evolution of exceptional lifespan among mole-rats and the ecological factors that may have accompanied that evolution.

Calcification of large arteries is a high-risk factor in the development of cardiovascular diseases, however, due to the lack of routine monitoring, the pathology remains severely under-diagnosed and prevalence in the general population is not known. We have developed a set of machine learning methods to quantitate levels of abdominal aortic calcification (AAC) in the UK Biobank imaging cohort and carried out the largest to-date analysis of genetic, biochemical, and epidemiological risk factors associated with the pathology. In a genetic association study, we identified three novel loci associated with AAC (FGF9, NAV9, and APOE), and replicated a previously reported association at the TWIST1/HDAC9 locus. We find that AAC is a highly prevalent pathology, with ~ 1 in 10 adults above the age of 40 showing significant levels of hydroxyapatite build-up (Kauppila score > 3). Presentation of AAC was strongly predictive of future cardiovascular events including stenosis of precerebral arteries (HR~1.5), myocardial infarction (HR~1.3), ischemic heart disease (HR~1.3), as well as other diseases such as chronic obstructive pulmonary disease (HR~1.3). Significantly, we find that the risk for myocardial infarction from elevated AAC (HR ~1.4) was comparable to the risk of hypercholesterolemia (HR~1.4), yet most people who develop AAC are not hypercholesterolemic. Furthermore, the overwhelming majority (98%) of individuals who develop pathology do so in the absence of known pre-existing risk conditions such as chronic kidney disease and diabetes (0.6% and 2.7% respectively). Our findings indicate that despite the high cardiovascular risk, calcification of large arteries remains a largely under-diagnosed lethal condition, and there is a clear need for increased awareness and monitoring of the pathology in the general population.

Frailty indexes provide quantitative measurements of non-specific health decline and are particularly useful as longitudinal monitors of pre-mortal morbidity in aging studies. For mouse studies, frailty assessments can be taken non-invasively, but they require handling and direct observation that is labor-intensive to the scientist and stress-inducing to the animal. Here, we implement, evaluate, and provide a digital frailty index composed entirely of computational analyses of home-cage video and compare it to manually obtained frailty scores in genetically diverse mice. We show that the frailty scores assigned by our digital index correlate with both manually obtained frailty scores and chronological age. Thus, we provide a tool for frailty assessment that reduces stress to the animal and can be collected consistently, at scale, without substantial labor cost.

In recent years, improved sequencing technology and computational tools have made de novo genome assembly more accessible. Many approaches, however, generate either an unphased or only partially resolved representation of a diploid genome, in which polymorphisms are detected but not assigned to one or the other of the homologous chromosomes. Yet chromosomal phase information is invaluable for the understanding of phenotypic trait inheritance in the cases of compound heterozygosity, allele-specific expression or cis-acting variants. Here we use a combination of tools and sequencing technologies to generate a de novo diploid assembly of the human primary cell line WI-38. First, data from PacBio single molecule sequencing and Bionano Genomics optical mapping were combined to generate an unphased assembly. Next, 10x Genomics linked reads were combined with the hybrid assembly to generate a partially phased assembly. Lastly, we developed and optimized methods to use short-read (Illumina) sequencing of flow cytometry-sorted metaphase chromosomes to provide phase information. The final genome assembly was almost fully (94%) phased with the addition of approximately 2.5-fold coverage of Illumina data from the sequenced metaphase chromosomes. The diploid nature of the final de novo genome assembly improved the resolution of structural variants between the WI-38 genome and the human reference genome. The phased WI-38 sequence data are available for browsing and download at wi38.research.calicolabs.com. Our work shows that assembling a completely phased diploid genome de novo from the DNA of a single individual is now readily achievable.

Funder: Calico Life Sciences, LLC, The naked mole‐rat (Heterocephalus glaber) has fascinated zoologists for at least half a century. It has also generated considerable biomedical interest not only because of its extraordinary longevity, but also because of unusual protective features (e.g. its tolerance of variable oxygen availability), which may be pertinent to several human disease states, including ischemia/reperfusion injury and neurodegeneration. A recent article entitled ‘Surprisingly long survival of premature conclusions about naked mole‐rat biology’ described 28 ‘myths’ which, those authors claimed, are a ‘perpetuation of beautiful, but falsified, hypotheses’ and impede our understanding of this enigmatic mammal. Here, we re‐examine each of these ‘myths’ based on evidence published in the scientific literature. Following Braude et al., we argue that these ‘myths’ fall into four main categories: (i) ‘myths’ that would be better described as oversimplifications, some of which persist solely in the popular press; (ii) ‘myths’ that are based on incomplete understanding, where more evidence is clearly needed; (iii) ‘myths’ where the accumulation of evidence over the years has led to a revision in interpretation, but where there is no significant disagreement among scientists currently working in the field; ( iv ) ‘myths’ where there is a genuine difference in opinion among active researchers, based on alternative interpretations of the available evidence. The term ‘myth’ is particularly inappropriate when applied to competing, evidence‐based hypotheses, which form part of the normal evolution of scientific knowledge. Here, we provide a comprehensive critical review of naked mole‐rat biology and attempt to clarify some of these misconceptions.

In recent years, improved sequencing technology and computational tools have made de novo genome assembly more accessible. Many approaches, however, generate either an unphased or only partially resolved representation of a diploid genome, in which polymorphisms are detected but not assigned to one or the other of the homologous chromosomes. Yet chromosomal phase information is invaluable for the understanding of phenotypic trait inheritance in the cases of compound heterozygosity, allele-specific expression or cis-acting variants. Here we use a combination of tools and sequencing technologies to generate a de novo diploid assembly of the human primary cell line WI-38. First, data from PacBio single molecule sequencing and Bionano Genomics optical mapping were combined to generate an unphased assembly. Next, 10x Genomics linked reads were combined with the hybrid assembly to generate a partially phased assembly. Lastly, we developed and optimized methods to use short-read (Illumina) sequencing of flow cytometry-sorted metaphase chromosomes to provide phase information. The final genome assembly was almost fully (94%) phased with the addition of approximately 2.5-fold coverage of Illumina data from the sequenced metaphase chromosomes. The diploid nature of the final de novo genome assembly improved the resolution of structural variants between the WI-38 genome and the human reference genome. The phased WI-38 sequence data are available for browsing and download at wi38.research.calicolabs.com. Our work shows that assembling a completely phased diploid genome de novo from the DNA of a single individual is now readily achievable.

BackgroundCalcification of the abdominal artery is an important contributor to cardiovascular disease in diabetic and chronic kidney disease (CKD) populations. However, prevalence of the pathology, risk factors, and long term disease outcomes in a general population have not been systematically analyzed.MethodWe developed machine learning models to quantify levels of abdominal aortic calcification (AAC) in 29,957 whole body dual-energy X-ray absorptiometry (DEXA) scans from the UK Biobank cohort. Using regression techniques we associated severity of calcification across a wide range of physiological parameters, clinical biomarkers, and environmental risk factors (406 in total). We performed a common variant genetic association study spanning 9,572,557 single-nucleotide polymorphisms to identify genetic loci relevant to AAC. We evaluated the prognostic value of AAC across 151 disease classes using Cox proportional hazard models. We further examined an epidemiological model of calcification on cardiovascular morbidity with and without LDL interactions.FindingsWe find evidence for AAC in >10.4% of the cohort despite low prevalence of diabetes (2.5%) and CKD (0.5%). Increased level of AAC is a strong prognostic indicator of cardiovascular outcomes for stenosis of precerebral arteries (HR~1.5), Myocardial Infarction (HR~1.5), & Ischemic Heart Disease (HR~1.33). We find that AAC is genetically correlated with cardiovascular-related traits and that the genetic signals are enriched in vascular and adipose tissue. We report three loci associated with AAC, with the strongest association occuring at theTWIST1/HDAC9locus (beta=0.078, p-value=1.4e-11) in a region also associated with coronary artery disease. Surprisingly, we find that elevated but still within clinically normal levels of serum phosphate and glycated hemoglobin are linked to increased vascular calcification. Furthermore, we show AAC arises in the absence of hypercholesterolemia. By our estimate, AAC is an LDL-independent risk factor for cardiovascular outcomes, with risk similar to elevated LDL.DataThis research has been conducted using the UK Biobank Resource.

Here, Ruby et al. analyze an unprecedented amount of public family tree data from Ancestry and determine that the heritability of human longevity was well below 10%, lower than the widely-held belief that lifespan..., Human life span is a phenotype that integrates many aspects of health and environment into a single ultimate quantity: the elapsed time between birth and death. Though it is widely believed that long life runs in families for genetic reasons, estimates of life span “heritability” are consistently low (∼15–30%). Here, we used pedigree data from Ancestry public trees, including hundreds of millions of historical persons, to estimate the heritability of human longevity. Although “nominal heritability” estimates based on correlations among genetic relatives agreed with prior literature, the majority of that correlation was also captured by correlations among nongenetic (in-law) relatives, suggestive of highly assortative mating around life span-influencing factors (genetic and/or environmental). We used structural equation modeling to account for assortative mating, and concluded that the true heritability of human longevity for birth cohorts across the 1800s and early 1900s was well below 10%, and that it has been generally overestimated due to the effect of assortative mating.

Non-traditional model organisms can facilitate discovery when their natural properties provide insight into biological mechanisms that are invariant across standard-use lab species. Long-lived naked mole-rats provide such insights for healthy aging.

We present a massive investigation into the genetic basis of human lifespan. Beginning with a genome-wide association (GWA) study using a de-identified snapshot of the unique AncestryDNA database – more than 300,000 genotyped individuals linked to pedigrees of over 400,000,000 people – we mapped six genome-wide significant loci associated with parental lifespan. We compared these results to a GWA analysis of the traditional lifespan proxy trait, age, and found only one locus, APOE, to be associated with both age and lifespan. By combining the AncestryDNA results with those of an independent UK Biobank dataset, we conducted a meta-analysis of more than 650,000 individuals and identified fifteen parental lifespan-associated loci. Beyond just those significant loci, our genome-wide set of polymorphisms accounts for up to 8% of the variance in human lifespan; this value represents a large fraction of the heritability estimated from phenotypic correlations between relatives.

Acoel-regeneration regulatory landscapes Some animals, including some types of worms, can undergo whole-body regeneration and replace virtually any missing cell type. Gehrke et al. sequenced and assembled the genome of Hofstenia miamia , a regenerative acoel worm species (see the Perspective by Alonge and Schatz). They identified a variable motif corresponding to regulation of the early growth response ( egr ) gene that was involved in regeneration. RNA interference experiments and validation in a second species showed that the protein Egr is a pioneer factor that stimulates regeneration. Science , this issue p. eaau6173 ; see also p. 1152

The longest-lived rodent, the naked mole-rat (Heterocephalus glaber), has a reported maximum lifespan of >30 years and exhibits delayed and/or attenuated age-associated physiological declines. We questioned whether these mouse-sized, eusocial rodents conform to Gompertzian mortality laws by experiencing an exponentially increasing risk of death as they get older. We compiled and analyzed a large compendium of historical naked mole-rat lifespan data with >3000 data points. Kaplan-Meier analyses revealed a substantial portion of the population to have survived at 30 years of age. Moreover, unlike all other mammals studied to date, and regardless of sex or breeding-status, the age-specific hazard of mortality did not increase with age, even at ages 25-fold past their time to reproductive maturity. This absence of hazard increase with age, in defiance of Gompertz’s law, uniquely identifies the naked mole-rat as a non-aging mammal, confirming its status as an exceptional model for biogerontology.

Although endogenous siRNAs (endo-siRNAs) have been described in many species, still little is known about their endogenous utility. Here, we show that Drosophila hairpin RNAs (hpRNAs) generate an endo-siRNA class with predominant expression in testes. Although hpRNAs are universally recently evolved, we identify highly complementary protein-coding targets for all hpRNAs. Importantly, we find broad evidence for evolutionary divergences that preferentially maintain compensatory pairing between hpRNAs and targets, serving as first evidence for adaptive selection for siRNA-mediated target regulation in metazoans. We demonstrate organismal impact of hpRNA activity, since knockout of hpRNA1 derepresses its target ATP synthase-β in testes and compromises spermatogenesis and male fertility. Moreover, we reveal surprising male-specific impact of RNAi factors on germ cell development and fertility, consistent with testis-directed function of the hpRNA pathway. Finally, the collected hpRNA loci chronicle an evolutionary timeline that reflects their origins from prospective target genes, mirroring a strategy described for plant miRNAs.

Deep sequencing was used to discover a novel rhabdovirus (Bas-Congo virus, or BASV) associated with a 2009 outbreak of 3 human cases of acute hemorrhagic fever in Mangala village, Democratic Republic of Congo (DRC), Africa. The cases, presenting over a 3-week period, were characterized by abrupt disease onset, high fever, mucosal hemorrhage, and, in two patients, death within 3 days. BASV was detected in an acute serum sample from the lone survivor at a concentration of 1.09 × 10(6) RNA copies/mL, and 98.2% of the genome was subsequently de novo assembled from ≈ 140 million sequence reads. Phylogenetic analysis revealed that BASV is highly divergent and shares less than 34% amino acid identity with any other rhabdovirus. High convalescent neutralizing antibody titers of1:1000 were detected in the survivor and an asymptomatic nurse directly caring for him, both of whom were health care workers, suggesting the potential for human-to-human transmission of BASV. The natural animal reservoir host or arthropod vector and precise mode of transmission for the virus remain unclear. BASV is an emerging human pathogen associated with acute hemorrhagic fever in Africa.

The giant, single-celled organism Stentor coeruleus has a long history as a model system for studying pattern formation and regeneration in single cells. Stentor (Figure 1A,B [1,2]) is a heterotrichous ciliate distantly related to familiar ciliate models such as Tetrahymena or Paramecium. The primary distinguishing feature of Stentor is its incredible size: a single cell is 1 millimeter long. Early developmental biologists, including T.H. Morgan[3], were attracted to the system because of its regenerative abilities -- if large portions of a cell are surgically removed, the remnant reorganizes into a normal-looking but smaller cell with correct proportionality [2,3]. These biologists were also drawn to Stentor because it exhibits a rich repertoire of behaviors, including light avoidance, mechanosensitive contraction, food selection, and even the ability to habituate to touch, a simple form of learning usually seen in higher organisms [4]. While early microsurgical approaches demonstrated a startling array of regenerative and morphogenetic processes in this single-celled organism, Stentor was never developed as a molecular model system. We report the sequencing of the Stentor coeruleus macronuclear genome and reveal key features of the genome: First, we find that Stentor uses the standard genetic code, suggesting that ciliate specific genetic codes arose after Stentor branched from other ciliates. We also discover that ploidy correlates with Stentor’s cell size. Finally, in the Stentor genome, we discover the smallest spliceosomal introns reported for any species. The sequenced genome opens the door to molecular analysis of single-cell regeneration in Stentor.

Low-cost DNA sequencing technologies have expanded the role for direct nucleic acid sequencing in the analysis of genomes, transcriptomes, and the metagenomes of whole ecosystems. Human and machine comprehension of such large datasets can be simplified via synthesis of sequence fragments into long, contiguous blocks of sequence (contigs), but most of the progress in the field of assembly has focused on genomes in isolation rather than metagenomes. Here, we present software for paired-read iterative contig extension (PRICE), a strategy for focused assembly of particular nucleic acid species using complex metagenomic data as input. We describe the assembly strategy implemented by PRICE and provide examples of its application to the sequence of particular genes, transcripts, and virus genomes from complex multicomponent datasets, including an assembly of the BCBL-1 strain of Kaposi’s sarcoma-associated herpesvirus. PRICE is open-source and available for free download (derisilab.ucsf.edu/software/price/ or sourceforge.net/projects/pricedenovo/).

Meroterpenoids are mixed polyketide-terpenoid natural products with a broad range of biological activities. Herein, we present the structures of four new meroterpenoid antibiotics, merochlorins A-D, produced by the marine bacterium Streptomyces sp. strain CNH-189, which possess novel chemical skeletons unrelated to known bacterial agents. Draft genome sequencing, mutagenesis, and heterologous biosynthesis in the genome-minimized model actinomycete Streptomyces coelicolor provided the 57.6 kb merochlorin gene cluster that contains two genes encoding rare bacterial vanadium-dependent haloperoxidase (VHPO) genes. Pathway expression of two different fosmid clones that differ largely by the presence or absence of the VHPO gene mcl40 resulted in the differential biosynthesis of merochlorin C, suggesting that Mcl40 catalyzes an unprecedented 15-membered chloronium-induced macrocyclization reaction converting merochlorin D to merochlorin C.

Post-transcriptional gene regulation frequently occurs through elements in mRNA 3′ untranslated regions (UTRs)1, 2. Although crucial roles for 3′UTR-mediated gene regulation have been found in Caenorhabditis elegans3, 4, 5, most C. elegans genes have lacked annotated 3′UTRs6, 7. Here we describe a high-throughput method for reliable identification of polyadenylated RNA termini, and we apply this method, called poly(A)-position profiling by sequencing (3P-Seq), to determine C. elegans 3′UTRs. Compared to standard methods also recently applied to C. elegans UTRs8, 3P-Seq identified 8,580 additional UTRs while excluding thousands of shorter UTR isoforms that do not seem to be authentic. Analysis of this expanded and corrected data set suggested that the high A/U content of C. elegans 3′UTRs facilitated genome compaction, because the elements specifying cleavage and polyadenylation, which are A/U rich, can more readily emerge in A/U-rich regions. Indeed, 30% of the protein-coding genes have mRNAs with alternative, partially overlapping end regions that generate another 10,480 cleavage and polyadenylation sites that had gone largely unnoticed and represent potential evolutionary intermediates of progressive UTR shortening. Moreover, a third of the convergently transcribed genes use palindromic arrangements of bidirectional elements to specify UTRs with convergent overlap, which also contributes to genome compaction by eliminating regions between genes. Although nematode 3′UTRs have median length only one-sixth that of mammalian 3′UTRs, they have twice the density of conserved microRNA sites, in part because additional types of seed-complementary sites are preferentially conserved. These findings reveal the influence of cleavage and polyadenylation on the evolution of genome architecture and provide resources for studying post-transcriptional gene regulation., National Institutes of Health (U.S.) (Grant number GM067031), National Science Foundation (U.S.). Predoctural Fellowship, United States. Dept. of Energy. Computational Science Graduate Fellowship (Krell Institute)

In metazoans Piwi-related Argonaute proteins have been linked to germ-line maintenance, and to a class of germ-line enriched small RNAs termed piRNAs. Here we show that an abundant class of 21-nucleotide small RNAs (21U-RNAs) are expressed in the C. elegans germ line, interact with the C. elegans Piwi-family member, PRG-1, and depend on PRG-1 activity for their accumulation. The PRG-1 protein is expressed throughout development and localizes to nuage-like structures called P-granules. Although 21U-RNA loci share a conserved upstream sequence motif, the mature 21U-RNAs are not conserved and, with few exceptions, fail to exhibit complementarity or evidence for direct regulation of other expressed sequences. Our findings demonstrate that 21U-RNAs are the piRNAs of C. elegans and link this class of small RNAs and their associated Piwi Argonaute to the maintenance of temperature-dependent fertility.

In contrast to microRNAs and Piwi-associated RNAs, short interfering RNAs (siRNAs) are seemingly dispensable for host-directed gene regulation in Drosophila. This notion is based on the fact that mutants lacking the core siRNA-generating enzyme Dicer-2 or the predominant siRNA effector Argonaute 2 are viable, fertile and of relatively normal morphology. Moreover, endogenous Drosophila siRNAs have not yet been identified. Here we report that siRNAs derived from long hairpin RNA genes (hpRNAs) programme Slicer complexes that can repress endogenous target transcripts. The Drosophila hpRNA pathway is a hybrid mechanism that combines canonical RNA interference factors (Dicer-2, Hen1 (known as CG12367) and Argonaute 2) with a canonical microRNA factor (Loquacious) to generate approximately 21-nucleotide siRNAs. These novel regulatory RNAs reveal unexpected complexity in the sorting of small RNAs, and open a window onto the biological usage of endogenous RNA interference in Drosophila.

MicroRNA (miRNA) genes give rise to small regulatory RNAs in a wide variety of organisms. We used computational methods to predict miRNAs conserved among Drosophila species and large-scale sequencing of small RNAs from Drosophila melanogaster to experimentally confirm and complement these predictions. In addition to validating 20 of our top 45 predictions for novel miRNA loci, the large-scale sequencing identified many miRNAs that had not been predicted. In total, 59 novel genes were identified, increasing our tally of confirmed fly miRNAs to 148. The large-scale sequencing also refined the identities of previously known miRNAs and provided insights into their biogenesis and expression. Many miRNAs were expressed in particular developmental contexts, with a large cohort of miRNAs expressed primarily in imaginal discs. Conserved miRNAs typically were expressed more broadly and robustly than were nonconserved miRNAs, and those conserved miRNAs with more restricted expression tended to have fewer predicted targets than those expressed more broadly. Predicted targets for the expanded set of microRNAs substantially increased and revised the miRNA-target relationships that appear conserved among the fly species. Insights were also provided into miRNA gene evolution, including evidence for emergent regulatory function deriving from the opposite arm of the miRNA hairpin, exemplified by mir-10, and even the opposite strand of the DNA, exemplified by mir-iab-4.

We sequenced approximately 400,000 small RNAs from Caenorhabditis elegans. Another 18 microRNA (miRNA) genes were identified, thereby extending to 112 our tally of confidently identified miRNA genes in C. elegans. Also observed were thousands of endogenous siRNAs generated by RNA-directed RNA polymerases acting preferentially on transcripts associated with spermatogenesis and transposons. In addition, a third class of nematode small RNAs, called 21U-RNAs, was discovered. 21U-RNAs are precisely 21 nucleotides long, begin with a uridine 5'-monophosphate but are diverse in their remaining 20 nucleotides, and appear modified at their 3'-terminal ribose. 21U-RNAs originate from more than 5700 genomic loci dispersed in two broad regions of chromosome IV-primarily between protein-coding genes or within their introns. These loci share a large upstream motif that enables accurate prediction of additional 21U-RNAs. The motif is conserved in other nematodes, presumably because of its importance for producing these diverse, autonomously expressed, small RNAs (dasRNAs).

Morphogenesis and pattern formation are vital processes in any organism, whether unicellular or multicellular. But in contrast to the developmental biology of plants and animals, the principles of morphogenesis and pattern formation in single cells remain largely unknown. Although all cells develop patterns, they are most obvious in ciliates; hence, we have turned to a classical unicellular model system, the giant ciliate Stentor coeruleus. Here we show that the RNA interference (RNAi) machinery is conserved in Stentor. Using RNAi, we identify the kinase coactivator Mob1--with conserved functions in cell division and morphogenesis from plants to humans-as an asymmetrically localized patterning protein required for global patterning during development and regeneration in Stentor. Our studies reopen the door for Stentor as a model regeneration system.

Berninamycin is a member of the pyridine-containing thiopeptide class of antibiotics that undergoes massive posttranslational modifications from ribosomally generated preproteins. Berninamycin has a 2-oxazolyl-3-thiazolyl-pyridine core embedded in a 35-atom macrocycle rather than typical trithiazolylpyridine cores embedded in 26-atom and 29-atom peptide macrocycles. We describe the cloning of an 11-gene berninamycin cluster from Streptomyces bernensis UC 5144, its heterologous expression in Streptomyces lividans TK24 and Streptomyces venezuelae ATCC 10712, and detection of variant and incompletely processed scaffolds. Posttranslational maturation in S. lividans of both the wild-type berninamycin prepeptide (BerA) and also a T3A mutant generates macrocyclic compounds as well as linear variants, which have failed to form the pyridine and the macrocycle. Expression of the gene cluster in S. venezuelae generates a variant of the 35-atom skeleton of berninamycin, containing a methyloxazoline in the place of a methyloxazole within the macrocyclic framework.

Background - The process of generating raw genome sequence data continues to become cheaper, faster, and more accurate. However, assembly of such data into high-quality, finished genome sequences remains challenging. Many genome assembly tools are available, but they differ greatly in terms of their performance (speed, scalability, hardware requirements, acceptance of newer read technologies) and in their final output (composition of assembled sequence). More importantly, it remains largely unclear how to best assess the quality of assembled genome sequences. The Assemblathon competitions are intended to assess current state-of-the-art methods in genome assembly. Results - In Assemblathon 2, we provided a variety of sequence data to be assembled for three vertebrate species (a bird, a fish, and snake). This resulted in a total of 43 submitted assemblies from 21 participating teams. We evaluated these assemblies using a combination of optical map data, Fosmid sequences, and several statistical methods. From over 100 different metrics, we chose ten key measures by which to assess the overall quality of the assemblies. Conclusions - Many current genome assemblers produced useful assemblies, containing a significant representation of their genes, regulatory sequences, and overall genome structure. However, the high degree of variability between the entries suggests that there is still much room for improvement in the field of genome assembly and that approaches which work well in assembling the genome of one species may not necessarily work well for another., Comment: Additional files available at http://korflab.ucdavis.edu/Datasets/Assemblathon/Assemblathon2/Additional_files/ Major changes 1. Accessions for the 3 read data sets have now been included 2. New file: spreadsheet containing details of all Study, Sample, Run, & Experiment identifiers 3. Made miscellaneous changes to address reviewers comments. DOIs added to GigaDB datasets

Deep sequencing was used to discover a novel rhabdovirus (Bas-Congo virus, or BASV) associated with a 2009 outbreak of 3 human cases of acute hemorrhagic fever in Mangala village, Democratic Republic of Congo (DRC), Africa. The cases, presenting over a 3-week period, were characterized by abrupt disease onset, high fever, mucosal hemorrhage, and, in two patients, death within 3 days. BASV was detected in an acute serum sample from the lone survivor at a concentration of 1.09×106 RNA copies/mL, and 98.2% of the genome was subsequently de novo assembled from ∼140 million sequence reads. Phylogenetic analysis revealed that BASV is highly divergent and shares less than 34% amino acid identity with any other rhabdovirus. High convalescent neutralizing antibody titers of >1∶1000 were detected in the survivor and an asymptomatic nurse directly caring for him, both of whom were health care workers, suggesting the potential for human-to-human transmission of BASV. The natural animal reservoir host or arthropod vector and precise mode of transmission for the virus remain unclear. BASV is an emerging human pathogen associated with acute hemorrhagic fever in Africa., Author Summary We used deep sequencing, a method for generating millions of DNA sequence reads from clinical samples, to discover a novel rhabdovirus (Bas-Congo virus, or BASV) associated with a 2009 outbreak of 3 human cases of acute hemorrhagic fever in Mangala village, Democratic Republic of Congo (DRC), Africa. The cases, presenting over a 3-week period, were characterized by abrupt disease onset, high fever, bloody vomiting and diarrhea, and, in two patients, death within 3 days. BASV was present in the blood of the lone survivor at a concentration of over a million copies per milliliter. The genome of BASV, assembled from over 140 million sequence reads, reveals that it is very different from any other rhabdovirus. The lone survivor and a nurse caring for him (with no symptoms), both health care workers, were found to have high levels of antibodies to BASV, indicating that they both had been infected by the virus. Although the source of the virus remains unclear, our study findings suggest that BASV may be spread by human-to-human contact and is an emerging pathogen associated with acute hemorrhagic fever in Africa.

Inclusion body disease (IBD) is an infectious fatal disease of snakes typified by behavioral abnormalities, wasting, and secondary infections. At a histopathological level, the disease is identified by the presence of large eosinophilic cytoplasmic inclusions in multiple tissues. To date, no virus or other pathogen has been definitively characterized or associated with the disease. Using a metagenomic approach to search for candidate etiologic agents in snakes with confirmed IBD, we identified and de novo assembled the complete genomic sequences of two viruses related to arenaviruses, and a third arenavirus-like sequence was discovered by screening an additional set of samples. A continuous boa constrictor cell line was established and used to propagate and isolate one of the viruses in culture. Viral nucleoprotein was localized and concentrated within large cytoplasmic inclusions in infected cells in culture and tissues from diseased snakes. In total, viral RNA was detected in 6/8 confirmed IBD cases and 0/18 controls. These viruses have a typical arenavirus genome organization but are highly divergent, belonging to a lineage separate from that of the Old and New World arenaviruses. Furthermore, these viruses encode envelope glycoproteins that are more similar to those of filoviruses than to those of other arenaviruses. These findings implicate these viruses as candidate etiologic agents of IBD. The presence of arenaviruses outside mammals reveals that these viruses infect an unexpectedly broad range of species and represent a new reservoir of potential human pathogens., IMPORTANCE Inclusion body disease (IBD) is a common infectious disease of captive snakes. IBD is fatal and can cause the loss of entire animal collections. The cause of the disease has remained elusive, and no treatment exists. In addition to being important to pet owners, veterinarians, breeders, zoological parks, and aquariums, the study of animal disease is significant since animals are the source of virtually every emerging infectious human disease. We searched for candidate causative agents in snakes diagnosed with IBD and found a group of novel viruses distantly related mainly to arenaviruses but also to filoviruses, both of which can cause fatal hemorrhagic fevers when transmitted from animals to humans. In addition to providing evidence that strongly suggests that these viruses cause snake IBD, this discovery reveals a new and unanticipated domain of virus biology and evolution.

Low-cost short read sequencing technology has revolutionized genomics, though it is only just becoming practical for the high-quality de novo assembly of a novel large genome. We describe the Assemblathon 1 competition, which aimed to comprehensively assess the state of the art in de novo assembly methods when applied to current sequencing technologies. In a collaborative effort, teams were asked to assemble a simulated Illumina HiSeq data set of an unknown, simulated diploid genome. A total of 41 assemblies from 17 different groups were received. Novel haplotype aware assessments of coverage, contiguity, structure, base calling, and copy number were made. We establish that within this benchmark: (1) It is possible to assemble the genome to a high level of coverage and accuracy, and that (2) large differences exist between the assemblies, suggesting room for further improvements in current methods. The simulated benchmark, including the correct answer, the assemblies, and the code that was used to evaluate the assemblies is now public and freely available from http://www.assemblathon.org/.

MicroRNAs (miRNAs) are small regulatory RNAs that derive from distinctive hairpin transcripts. To learn more about the miRNAs of mammals, we sequenced 60 million small RNAs from mouse brain, ovary, testes, embryonic stem cells, three embryonic stages, and whole newborns. Analysis of these sequences confirmed 398 annotated miRNA genes and identified 108 novel miRNA genes. More than 150 previously annotated miRNAs and hundreds of candidates failed to yield sequenced RNAs with miRNA-like features. Ectopically expressing these previously proposed miRNA hairpins also did not yield small RNAs, whereas ectopically expressing the confirmed and newly identified hairpins usually did yield small RNAs with the classical miRNA features, including dependence on the Drosha endonuclease for processing. These experiments, which suggest that previous estimates of conserved mammalian miRNAs were inflated, provide a substantially revised list of confidently identified murine miRNAs from which to infer the general features of mammalian miRNAs. Our analyses also revealed new aspects of miRNA biogenesis and modification, including tissue-specific strand preferences, sequential Dicer cleavage of a metazoan precursor miRNA (pre-miRNA), consequential 5′ heterogeneity, newly identified instances of miRNA editing, and evidence for widespread pre-miRNA uridylation reminiscent of miRNA regulation by Lin28.

Canonical microRNAs (miRNAs) require two processing steps: the first by the Microprocessor, a complex of DGCR8 and Drosha, and the second by a complex of TRBP and Dicer. dgcr8Δ/Δ mouse embryonic stem cells (mESCs) have less severe phenotypes than dicer1Δ/Δ mESCs, suggesting a physiological role for Microprocessor-independent, Dicer-dependent small RNAs. To identify these small RNAs with unusual biogenesis, we performed high-throughput sequencing from wild-type, dgcr8Δ/Δ, and dicer1Δ/Δ mESCs. Several of the resulting DGCR8-independent, Dicer-dependent RNAs were noncanonical miRNAs. These derived from mirtrons and a newly identified subclass of miRNA precursors, which appears to be the endogenous counterpart of shRNAs. Our analyses also revealed endogenous siRNAs resulting from Dicer cleavage of long hairpins, the vast majority of which originated from one genomic locus with tandem, inverted short interspersed nuclear elements (SINEs). Our results extend the known diversity of mammalian small RNA-generating pathways and show that mammalian siRNAs exist in cell types other than oocytes.

High-throughput short-read technologies have revolutionized DNA sequencing by drastically reducing the cost per base of sequencing information. Despite producing gigabases of sequence per run, these technologies still present obstacles in resequencing and de novo assembly applications due to biased or insufficient target sequence coverage. We present here a simple sample preparation method termed the ‘‘long march’’ that increases both contig lengths and target sequence coverage using high-throughput short-read technologies. By incorporating a Type IIS restriction enzyme recognition motif into the sequencing primer adapter, successive rounds of restriction enzyme cleavage and adapter ligation produce a set of nested sub-libraries from the initial amplicon library. Sequence reads from these sub-libraries are offset from each other with enough overlap to aid assembly and contig extension. We demonstrate the utility of the long march in resequencing of the Plasmodium falciparum transcriptome, where the number of genomic bases covered was increased by 39%, as well as in metagenomic analysis of a serum sample from a patient with hepatitis B virus (HBV)-related acute liver failure, where the number of HBV bases covered was increased by 42%. We also offer a theoretical optimization of the long march for de novo sequence assembly.

Sequencing of multiple related species followed by comparative genomics analysis constitutes a powerful approach for the systematic understanding of any genome. Here, we use the genomes of 12 Drosophila species for the de novo discovery of functional elements in the fly. Each type of functional element shows characteristic patterns of change, or 'evolutionary signatures', dictated by its precise selective constraints. Such signatures enable recognition of new protein-coding genes and exons, spurious and incorrect gene annotations, and numerous unusual gene structures, including abundant stop-codon readthrough. Similarly, we predict non-protein-coding RNA genes and structures, and new microRNA (miRNA) genes. We provide evidence of miRNA processing and functionality from both hairpin arms and both DNA strands. We identify several classes of pre- and post-transcriptional regulatory motifs, and predict individual motif instances with high confidence. We also study how discovery power scales with the divergence and number of species compared, and we provide general guidelines for comparative studies.

Made available in DSpace on 2014-12-03T13:10:41Z (GMT). No. of bitstreams: 0 Previous issue date: 2013-01-01Bitstream added on 2014-12-03T13:23:07Z : No. of bitstreams: 1 WOS000320155200010.pdf: 2040166 bytes, checksum: 8b98e8c7f6f77d18cab8c6df4dd9f7cb (MD5) Joint USDA/DOE Office of Science Feedstock genomics grant National Science Foundation NIH-NIGMS Background: Centromeres are essential for chromosome segregation, yet their DNA sequences evolve rapidly. In most animals and plants that have been studied, centromeres contain megabase-scale arrays of tandem repeats. Despite their importance, very little is known about the degree to which centromere tandem repeats share common properties between different species across different phyla. We used bioinformatic methods to identify high-copy tandem repeats from 282 species using publicly available genomic sequence and our own data.Results: Our methods are compatible with all current sequencing technologies. Long Pacific Biosciences sequence reads allowed us to find tandem repeat monomers up to 1,419 bp. We assumed that the most abundant tandem repeat is the centromere DNA, which was true for most species whose centromeres have been previously characterized, suggesting this is a general property of genomes. High-copy centromere tandem repeats were found in almost all animal and plant genomes, but repeat monomers were highly variable in sequence composition and length. Furthermore, phylogenetic analysis of sequence homology showed little evidence of sequence conservation beyond approximately 50 million years of divergence. We find that despite an overall lack of sequence conservation, centromere tandem repeats from diverse species showed similar modes of evolution.Conclusions: While centromere position in most eukaryotes is epigenetically determined, our results indicate that tandem repeats are highly prevalent at centromeres of both animal and plant genomes. This suggests a functional role for such repeats, perhaps in promoting concerted evolution of centromere DNA across chromosomes. Univ Calif Davis, Dept Mol & Cell Biol, Davis, CA 95616 USA Univ Calif Davis, Genome Ctr, Davis, CA 95616 USA Univ Calif Davis, Dept Plant Biol, Davis, CA 95616 USA USDA ARS, Western Reg Res Ctr, Albany, CA 94710 USA Univ Calif Davis, Dept Ecol & Evolut, Davis, CA 95616 USA Univ Calif San Francisco, Dept Biochem & Biophys, San Francisco, CA 94158 USA Pacific Biosci, Menlo Pk, CA 94025 USA Univ Estadual Paulista, Dept Anim Prod & Hlth, IAEA Collaborating Ctr Anim Genom & Bioinformat, BR-16050680 Aracatuba, SP, Brazil Howard Hughes Med Inst, Chevy Chase, MD 20815 USA USDA ARS, US Meat Anim Res Ctr, Clay Ctr, NE 68933 USA Univ Calif Davis, Dept Plant Sci, Ctr Populat Biol, Davis, CA 95616 USA Univ Estadual Paulista, Dept Anim Prod & Hlth, IAEA Collaborating Ctr Anim Genom & Bioinformat, BR-16050680 Aracatuba, SP, Brazil Joint USDA/DOE Office of Science Feedstock genomics grantDE-AI02-09ER64829 National Science FoundationIOS-0922703 National Science FoundationIOS-1026094 NIH-NIGMST32-GM008799

β-Human papillomavirus (β-HPV) DNA is present in some cutaneous squamous cell carcinomas (cuSCCs), but no mechanism of carcinogenesis has been determined. We used ultra-high-throughput sequencing of the cancer transcriptome to assess whether papillomavirus transcripts are present in these cancers. In all, 67 cuSCC samples were assayed for β-HPV DNA by PCR, and viral loads were measured with type-specific quantitative PCR. A total of 31 SCCs were selected for whole transcriptome sequencing. Transcriptome libraries were prepared in parallel from the HPV18-positive HeLa cervical cancer cell line and HPV16-positive primary cervical and periungual SCCs. Of the tumors, 30% (20/67) were positive for β-HPV DNA, but there was no difference in β-HPV viral load between tumor and normal tissue (P=0.310). Immunosuppression and age were significantly associated with higher viral load (P=0.016 for immunosuppression; P=0.0004 for age). Transcriptome sequencing failed to identify papillomavirus expression in any of the skin tumors. In contrast, HPV16 and HPV18 mRNA transcripts were readily identified in primary cervical and periungual cancers and HeLa cells. These data demonstrate that papillomavirus mRNA expression is not a factor in the maintenance of cuSCCs.

A colony of domestic rabbits in Tennessee, USA, experienced a high-mortality (~90%) outbreak of enterocolitis. The clinical characteristics were one to six days of lethargy, bloating, and diarrhea, followed by death. Heavy intestinal coccidial load was a consistent finding as was mucoid enteropathy with cecal impaction. Preliminary analysis by electron microscopy revealed the presence of virus-like particles in the stool of one of the affected rabbits. Analysis using the Virochip, a viral detection microarray, suggested the presence of an astrovirus, and follow-up PCR and sequence determination revealed a previously uncharacterized member of that family. Metagenomic sequencing enabled the recovery of the complete viral genome, which contains the characteristic attributes of astrovirus genomes. Attempts to propagate the virus in tissue culture have yet to succeed. Although astroviruses cause gastroenteric disease in other mammals, the pathogenicity of this virus and the relationship to this outbreak remains to be determined. This study therefore defines a viral species and a potential rabbit pathogen.