Intragenomic polymorphisms among high-copy loci: a genus-wide study of nuclear ribosomal DNA in Asclepias (Apocynaceae)

Despite knowledge that concerted evolution of high-copy loci is often imperfect, studies that investigate the extent of intragenomic polymorphisms and comparisons across a large number of species are rarely made. We present a bioinformatic pipeline for characterizing polymorphisms within an individual among copies of a high-copy locus. Results are presented for nuclear ribosomal DNA (nrDNA) across the milkweed genus, Asclepias. The 18S-26S portion of the nrDNA cistron of Asclepias syriaca served as a reference for assembly of the region from 124 samples representing 90 species of Asclepias. Reads were mapped back to each individual’s consensus and at each position reads differing from the consensus were tallied using a custom perl script. Low frequency polymorphisms existed in all individuals (mean = 5.8%). Most nrDNA positions (91%) were polymorphic in at least one individual, with polymorphic sites being less frequent in subunit regions and loops. Highly polymorphic sites existed in each individual, with highest abundance in the “noncoding” ITS regions. Phylogenetic signal was present in the distribution of intragenomic polymorphisms across the genus. Intragenomic polymorphisms in nrDNA are common in Asclepias, being found at higher frequency than any other study to date. The high and variable frequency of polymorphisms across species highlights concerns that phylogenetic applications of nrDNA may be error-prone. The new analytical approach provided here is applicable to other taxa and other high-copy regions characterized by low coverage genome sequencing (genome skimming).This is the publisher’s final pdf. The published article is copyrighted by the author(s) and published by PeerJ. The published article can be found at: https://peerj.com/.Keywords: Evolutionary studies, Bioinformatics, Intragenomic polymorphism, Plant science, Asclepias, Genomics, High-copy, Partial SNP (pSNP), ITS, Genetics, Concerted evolution, Nuclear ribosomal DNA (nrDNA), Intra-individual site polymorphism, Genome skimming, 2IS

Horizontal Transfer of DNA from the Mitochondrial to the Plastid Genome and Its Subsequent Evolution in Milkweeds (Apocynaceae)

Horizontal gene transfer (HGT) of DNA from the plastid to the nuclear and mitochondrial genomes of higher plants is a common phenomenon; however, plastid genomes (plastomes) are highly conserved and have generally been regarded as impervious to HGT. We sequenced the 158 kb plastome and the 690 kb mitochondrial genome of common milkweed (Asclepias syriaca [Apocynaceae]) and found evidence of intracellular HGT for a 2.4-kb segment of mitochondrial DNA to the rps2–rpoC2 intergenic spacer of the plastome. The transferred region contains an rpl2 pseudogene and is flanked by plastid sequence in the mitochondrial genome, including an rpoC2 pseudogene, which likely provided the mechanism for HGT back to the plastome through double-strand break repair involving homologous recombination. The plastome insertion is restricted to tribe Asclepiadeae of subfamily Asclepiadoideae, whereas the mitochondrial rpoC2 pseudogene is present throughout the subfamily, which confirms that the plastid to mitochondrial HGT event preceded the HGT to the plastome. Although the plastome insertion has been maintained in all lineages of Asclepiadoideae, it shows minimal evidence of transcription in A. syriaca and is likely nonfunctional. Furthermore, we found recent gene conversion of the mitochondrial rpoC2 pseudogene in Asclepias by the plastid gene, which reflects continued interaction of these genomes.Keywords: Mitochondrial genome, Phylogeny, Horizontal gene transfer, Gene conversion, Asclepias syriaca, PlastomeKeywords: Mitochondrial genome, Phylogeny, Horizontal gene transfer, Gene conversion, Asclepias syriaca, Plastom

Hyb-Seq: Combining target enrichment and genome skimming for plant phylogenomics

PREMISE OF THE STUDY: Hyb-Seq, the combination of target enrichment and genome skimming allows simultaneous data collection for low-copy nuclear genes and high-copy genomic targets for plant systematics and evolution studies. METHODS AND RESULTS: Genome and transcriptome assemblies for milkweed (Asclepias syriaca) were utilized to design enrichment probes for 3385 exons from 768 genes (>1.6 Mbp) followed by Illumina sequencing of enriched libraries. Hyb-Seq of twelve individuals (ten Asclepias species and two related genera) resulted in at least partial assembly of 92.6% of exons and 99.7% of genes and an average assembly length >2 Mbp. Importantly, complete plastomes and nrDNA cistrons were assembled using off-target reads. Phylogenomic analyses demonstrated signal conflict between genomes. CONCLUSIONS: The Hyb-Seq approach enables targeted sequencing of thousands of low-copy nuclear exons and flanking regions, as well as genome skimming of high-copy repeats and organellar genomes, to efficiently produce genome-scale datasets for phylogenomics.This is an author's peer-reviewed final manuscript, as accepted by the publisher. The article is in press and will be published in Applications in Plant Sciences, Vol. 2, no. 9, September 2014.Keywords: Genome skimming, Hyb-Seq, Target enrichment, Phylogenomics, Species tree, Nuclear lociKeywords: Genome skimming, Hyb-Seq, Target enrichment, Phylogenomics, Species tree, Nuclear loc

Phylogenetic marker development for target enrichment from transcriptome and genome skim data: the pipeline and its application in southern African Oxalis (Oxalidaceae)

Phylogenetics benefits from using a large number of putatively independent nuclear loci and their combination with other sources of information, such as the plastid and mitochondrial genomes. To facilitate the selection of orthologous low‐copy nuclear (LCN) loci for phylogenetics in nonmodel organisms, we created an automated and interactive script to select hundreds of LCN loci by a comparison between transcriptome and genome skim data. We used our script to obtain LCN genes for southern African Oxalis (Oxalidaceae), a speciose plant lineage in the Greater Cape Floristic Region. This resulted in 1164 LCN genes greater than 600 bp. Using target enrichment combined with genome skimming (Hyb‐Seq), we obtained on average 1141 LCN loci, nearly the whole plastid genome and the nrDNA cistron from 23 southern African Oxalis species. Despite a wide range of gene trees, the phylogeny based on the LCN genes was very robust, as retrieved through various gene and species tree reconstruction methods as well as concatenation. Cytonuclear discordance was strong. This indicates that organellar phylogenies alone are unlikely to represent the species tree and stresses the utility of Hyb‐Seq in phylogenetics

STAT3 in the systemic inflammation of cancer cachexia

Weight loss is diagnostic of cachexia, a debilitating syndrome contributing mightily to morbidity and mortality in cancer. Most research has probed mechanisms leading to muscle atrophy and adipose wasting in cachexia; however cachexia is a truly systemic phenomenon. Presence of the tumor elicits an inflammatory response and profound metabolic derangements involving not only muscle and fat, but also the hypothalamus, liver, heart, blood, spleen and likely other organs. This global response is orchestrated in part through circulating cytokines that rise in conditions of cachexia. Exogenous Interleukin-6 (IL6) and related cytokines can induce most cachexia symptomatology, including muscle and fat wasting, the acute phase response and anemia, while IL-6 inhibition reduces muscle loss in cancer. Although mechanistic studies are ongoing, certain of these cachexia phenotypes have been causally linked to the cytokine-activated transcription factor, STAT3, including skeletal muscle wasting, cardiac dysfunction and hypothalamic inflammation. Correlative studies implicate STAT3 in fat wasting and the acute phase response in cancer cachexia. Parallel data in non-cancer models and disease states suggest both pathological and protective functions for STAT3 in other organs during cachexia. STAT3 also contributes to cancer cachexia through enhancing tumorigenesis, metastasis and immune suppression, particularly in tumors associated with high prevalence of cachexia. This review examines the evidence linking STAT3 to multi-organ manifestations of cachexia and the potential and perils for targeting STAT3 to reduce cachexia and prolong survival in cancer patients

Mimosoid legume plastome evolution: IR expansion, tandem repeat expansions, and accelerated rate of evolution in clpP

The Leguminosae has emerged as a model for studying angiosperm plastome evolution because of its striking diversity of structural rearrangements and sequence variation. However, most of what is known about legume plastomes comes from few genera representing a subset of lineages in subfamily Papilionoideae. We investigate plastome evolution in subfamily Mimosoideae based on two newly sequenced plastomes (Inga and Leucaena) and two recently published plastomes (Acacia and Prosopis), and discuss the results in the context of other legume and rosid plastid genomes. Mimosoid plastomes have a typical angiosperm gene content and general organization as well as a generally slow rate of protein coding gene evolution, but they are the largest known among legumes. The increased length results from tandem repeat expansions and an unusual 13 kb IR-SSC boundary shift in Acacia and Inga. Mimosoid plastomes harbor additional interesting features, including loss of clpP intron1 in Inga, accelerated rates of evolution in clpP for Acacia and Inga, and dN/dS ratios consistent with neutral and positive selection for several genes. These new plastomes and results provide important resources for legume comparative genomics, plant breeding, and plastid genetic engineering, while shedding further light on the complexity of plastome evolution in legumes and angiosperms.This is the publisher’s final pdf. The published article is copyrighted by the author(s) and published by the Nature Publishing Group. The published article can be found at: http://www.nature.com/articles/srep1695

Building a model: developing genomic resources for common milkweed (Asclepias syriaca) with low coverage genome sequencing

<p>Abstract</p> <p>Background</p> <p>Milkweeds (<it>Asclepias </it>L.) have been extensively investigated in diverse areas of evolutionary biology and ecology; however, there are few genetic resources available to facilitate and compliment these studies. This study explored how low coverage genome sequencing of the common milkweed (<it>Asclepias syriaca </it>L.) could be useful in characterizing the genome of a plant without prior genomic information and for development of genomic resources as a step toward further developing <it>A. syriaca </it>as a model in ecology and evolution.</p> <p>Results</p> <p>A 0.5× genome of <it>A. syriaca </it>was produced using Illumina sequencing. A virtually complete chloroplast genome of 158,598 bp was assembled, revealing few repeats and loss of three genes: <it>accD, clpP</it>, and <it>ycf1</it>. A nearly complete rDNA cistron (18S-5.8S-26S; 7,541 bp) and 5S rDNA (120 bp) sequence were obtained. Assessment of polymorphism revealed that the rDNA cistron and 5S rDNA had 0.3% and 26.7% polymorphic sites, respectively. A partial mitochondrial genome sequence (130,764 bp), with identical gene content to tobacco, was also assembled. An initial characterization of repeat content indicated that Ty1/<it>copia</it>-like retroelements are the most common repeat type in the milkweed genome. At least one <it>A. syriaca </it>microread hit 88% of <it>Catharanthus roseus </it>(Apocynaceae) unigenes (median coverage of 0.29×) and 66% of single copy orthologs (COSII) in asterids (median coverage of 0.14×). From this partial characterization of the <it>A. syriaca </it>genome, markers for population genetics (microsatellites) and phylogenetics (low-copy nuclear genes) studies were developed.</p> <p>Conclusions</p> <p>The results highlight the promise of next generation sequencing for development of genomic resources for any organism. Low coverage genome sequencing allows characterization of the high copy fraction of the genome and exploration of the low copy fraction of the genome, which facilitate the development of molecular tools for further study of a target species and its relatives. This study represents a first step in the development of a community resource for further study of plant-insect co-evolution, anti-herbivore defense, floral developmental genetics, reproductive biology, chemical evolution, population genetics, and comparative genomics using milkweeds, and <it>A. syriaca </it>in particular, as ecological and evolutionary models.</p

Mapping the Conformational Dynamics and Pathways of Spontaneous Steric Zipper Peptide Oligomerization

The process of protein misfolding and self-assembly into various, polymorphic aggregates is associated with a number of important neurodegenerative diseases. Only recently, crystal structures of several short peptides have provided detailed structural insights into -sheet rich aggregates, known as amyloid fibrils. Knowledge about early events of the formation and interconversion of small oligomeric states, an inevitable step in the cascade of peptide self-assembly, however, remains still limited

Prevalence and architecture of de novo mutations in developmental disorders.

The genomes of individuals with severe, undiagnosed developmental disorders are enriched in damaging de novo mutations (DNMs) in developmentally important genes. Here we have sequenced the exomes of 4,293 families containing individuals with developmental disorders, and meta-analysed these data with data from another 3,287 individuals with similar disorders. We show that the most important factors influencing the diagnostic yield of DNMs are the sex of the affected individual, the relatedness of their parents, whether close relatives are affected and the parental ages. We identified 94 genes enriched in damaging DNMs, including 14 that previously lacked compelling evidence of involvement in developmental disorders. We have also characterized the phenotypic diversity among these disorders. We estimate that 42% of our cohort carry pathogenic DNMs in coding sequences; approximately half of these DNMs disrupt gene function and the remainder result in altered protein function. We estimate that developmental disorders caused by DNMs have an average prevalence of 1 in 213 to 1 in 448 births, depending on parental age. Given current global demographics, this equates to almost 400,000 children born per year

Prevalence, phenotype and architecture of developmental disorders caused by de novo mutation: The Deciphering Developmental Disorders Study

Individuals with severe, undiagnosed developmental disorders (DDs) are enriched for damaging de novo mutations (DNMs) in developmentally important genes. We exome sequenced 4,293 families with individuals with DDs, and meta-analysed these data with published data on 3,287 individuals with similar disorders. We show that the most significant factors influencing the diagnostic yield of de novo mutations are the sex of the affected individual, the relatedness of their parents and the age of both father and mother. We identified 94 genes enriched for damaging de novo mutation at genome-wide significance (P < 7 × 10−7), including 14 genes for which compelling data for causation was previously lacking. We have characterised the phenotypic diversity among these genetic disorders. We demonstrate that, at current cost differentials, exome sequencing has much greater power than genome sequencing for novel gene discovery in genetically heterogeneous disorders. We estimate that 42% of our cohort carry pathogenic DNMs (single nucleotide variants and indels) in coding sequences, with approximately half operating by a loss-of-function mechanism, and the remainder resulting in altered-function (e.g. activating, dominant negative). We established that most haplo insufficient developmental disorders have already been identified, but that many altered-function disorders remain to be discovered. Extrapolating from the DDD cohort to the general population, we estimate that developmental disorders caused by DNMs have an average birth prevalence of 1 in 213 to 1 in 448 (0.22-0.47% of live births), depending on parental age