Thiruvarangan Ramaraj

Genome, epigenome and RNA sequences of monozygotic twins discordant for multiple sclerosis

Monozygotic or ‘identical’ twins have been widely studied to dissect the relative contributions of genetics and environment in human diseases. In multiple sclerosis (MS), an autoimmune demyelinating disease and common cause of neurodegeneration and disability in young adults, disease discordance in monozygotic twins has been interpreted to indicate environmental importance in its pathogenesis. However, genetic and epigenetic differences between monozygotic twins have been described, challenging the accepted experimental model in disambiguating the effects of nature and nurture. Here we report the genome sequences of one MS-discordant monozygotic twin pair, and messenger RNA transcriptome and epigenome sequences of CD4+ lymphocytes from three MS-discordant, monozygotic twin pairs. No reproducible differences were detected between co-twins among ∼3.6 million single nucleotide polymorphisms (SNPs) or ∼0.2 million insertion-deletion polymorphisms. Nor were any reproducible differences observed between siblings of the three twin pairs in HLA haplotypes, confirmed MS-susceptibility SNPs, copy number variations, mRNA and genomic SNP and insertion-deletion genotypes, or the expression of ∼19,000 genes in CD4+ T cells. Only 2 to 176 differences in the methylation of ∼2 million CpG dinucleotides were detected between siblings of the three twin pairs, in contrast to ∼800 methylation differences between T cells of unrelated individuals and several thousand differences between tissues or between normal and cancerous tissues. In the first systematic effort to estimate sequence variation among monozygotic co-twins, we did not find evidence for genetic, epigenetic or transcriptome differences that explained disease discordance. These are the first, to our knowledge, female, twin and autoimmune disease individual genome sequences reported.

Genome and secretome analysis of the hemibiotrophic fungal pathogen, Moniliophthora roreri, which causes frosty pod rot disease of cacao: mechanisms of the biotrophic and necrotrophic phases.

BackgroundThe basidiomycete Moniliophthora roreri is the causal agent of Frosty pod rot (FPR) disease of cacao (Theobroma cacao), the source of chocolate, and FPR is one of the most destructive diseases of this important perennial crop in the Americas. This hemibiotroph infects only cacao pods and has an extended biotrophic phase lasting up to sixty days, culminating in plant necrosis and sporulation of the fungus without the formation of a basidiocarp.ResultsWe sequenced and assembled 52.3 Mb into 3,298 contigs that represent the M. roreri genome. Of the 17,920 predicted open reading frames (OFRs), 13,760 were validated by RNA-Seq. Using read count data from RNA sequencing of cacao pods at 30 and 60 days post infection, differential gene expression was estimated for the biotrophic and necrotrophic phases of this plant-pathogen interaction. The sequencing data were used to develop a genome based secretome for the infected pods. Of the 1,535 genes encoding putative secreted proteins, 1,355 were expressed in the biotrophic and necrotrophic phases. Analysis of the data revealed secretome gene expression that correlated with infection and intercellular growth in the biotrophic phase and invasive growth and plant cellular death in the necrotrophic phase.ConclusionsGenome sequencing and RNA-Seq was used to determine and validate the Moniliophthora roreri genome and secretome. High sequence identity between Moniliophthora roreri genes and Moniliophthora perniciosa genes supports the taxonomic relationship with Moniliophthora perniciosa and the relatedness of this fungus to other basidiomycetes. Analysis of RNA-Seq data from infected plant tissues revealed differentially expressed genes in the biotrophic and necrotrophic phases. The secreted protein genes that were upregulated in the biotrophic phase are primarily associated with breakdown of the intercellular matrix and modification of the fungal mycelia, possibly to mask the fungus from plant defenses. Based on the transcriptome data, the upregulated secreted proteins in the necrotrophic phase are hypothesized to be actively attacking the plant cell walls and plant cellular components resulting in necrosis. These genes are being used to develop a new understanding of how this disease interaction progresses and to identify potential targets to reduce the impact of this devastating disease.

Novel Meiotic miRNAs and Indications for a Role of PhasiRNAs in Meiosis

Small RNAs (sRNA) add additional layers to the regulation of gene expression, with siRNAs directing gene silencing at the DNA level by RdDM (RNA-directed DNA methylation), and micro RNAs (miRNAs) directing post-transcriptional regulation of specific target genes, mostly by mRNA cleavage. We used manually isolated male meiocytes from maize (Zea mays) to investigate sRNA and DNA methylation landscapes during zygotene, an early stage of meiosis during which steps of meiotic recombination and synapsis of paired homologous chromosomes take place. We discovered two novel miRNAs from meiocytes, zma-MIR11969 and zma-MIR11970, and identified putative target genes. Furthermore, we detected abundant phasiRNAs of 21 and 24 nt length. PhasiRNAs are phased small RNAs which occur in 21 or 24 nt intervals, at a few hundred loci, specifically in male reproductive tissues in grasses. So far, the function of phasiRNAs remained elusive. Data from isolated meiocytes now revealed elevated DNA methylation at phasiRNA loci, especially in the CHH context, suggesting a role for phasiRNAs in cis DNA methylation. In addition, we consider a role of these phasiRNAs in chromatin remodeling/dynamics during meiosis. However, this is not well supported yet and will need more additional data. Here, we only lay out the idea due to other relevant literature and our additional observation of a peculiar GC content pattern at phasiRNA loci. Chromatin remodeling is also indicated by the discovery that histone genes were enriched for sRNA of 22 nt length. Taken together, we gained clues that lead us to hypothesize sRNA-driven DNA methylation and possibly chromatin remodeling during male meiosis in the monocot maize which is in line with and extends previous knowledge.

The transcriptome landscape of early maize meiosis

BackgroundA major step in the higher plant life cycle is the decision to leave the mitotic cell cycle and begin the progression through the meiotic cell cycle that leads to the formation of gametes. The molecular mechanisms that regulate this transition and early meiosis remain largely unknown. To gain insight into gene expression features during the initiation of meiotic recombination, we profiled early prophase I meiocytes from maize (Zea mays) using capillary collection to isolate meiocytes, followed by RNA-seq.ResultsWe detected ~2,000 genes as preferentially expressed during early meiotic prophase, most of them uncharacterized. Functional analysis uncovered the importance of several cellular processes in early meiosis. Processes significantly enriched in isolated meiocytes included proteolysis, protein targeting, chromatin modification and the regulation of redox homeostasis. The most significantly up-regulated processes in meiocytes were processes involved in carbohydrate metabolism. Consistent with this, many mitochondrial genes were up-regulated in meiocytes, including nuclear- and mitochondrial-encoded genes. The data were validated with real-time PCR and in situ hybridization and also used to generate a candidate maize homologue list of known meiotic genes from Arabidopsis.ConclusionsTaken together, we present a high-resolution analysis of the transcriptome landscape in early meiosis of an important crop plant, providing support for choosing genes for detailed characterization of recombination initiation and regulation of early meiosis. Our data also reveal an important connection between meiotic processes and altered/increased energy production.

Composition, Diversity and Abundance of Gut Microbiome in Prediabetes and Type 2 Diabetes.

Association between type 2 diabetes (T2DM) and compositional changes in the gut micro biota is established, however little is known about the dysbiosis in early stages of Prediabetes (preDM). The purpose of this investigation is to elucidate the characteristics of the gut micro biome in preDM and T2DM, compared to Non-Diabetic (nonDM) subjects. Forty nine subjects were recruited for this study, 15 nonDM, 20 preDM and 14 T2DM. Bacterial community composition and diversity were investigated in fecal DNA samples using Illumina sequencing of the V4 region within the 16S rRNA gene. The five most abundant phyla identified were: Bacteroidetes, Firmicutes, Proteobacteria, Verrucomicrobia, and Actinobacteria. Class Chloracido bacteria was increased in preDM compared to T2DM (p = 0.04). An unknown genus from family Pseudonocardiaceae was significantly present in preDM group compared to the others (p = 0.04). Genus Collinsella, and an unknown genus belonging to family Enterobacteriaceae were both found to be significantly increased in T2DM compared to the other groups (Collinsella, and p = 0.03, Enterobacteriaceae genus p = 0.02). PERMANOVA and Mantel tests performed did not reveal a relationship between overall composition and diagnosis group or HbA1C level. This study identified dysbiosis associated with both preDM and T2DM, specifically at the class and genus levels suggesting that earlier treatment in preDM could possibly have an impact on the intestinal micro flora transitioning to T2DM.

Characterization, correction and de novo assembly of an Oxford Nanopore genomic dataset from Agrobacterium tumefaciens

The MinION is a portable single-molecule DNA sequencing instrument that was released by Oxford Nanopore Technologies in 2014, producing long sequencing reads by measuring changes in ionic flow when single-stranded DNA molecules translocate through the pores. While MinION long reads have an error rate substantially higher than the ones produced by short-read sequencing technologies, they can generate de novo assemblies of microbial genomes, after an initial correction step that includes alignment of Illumina sequencing data or detection of overlaps between Oxford Nanopore reads to improve accuracy. In this study, MinION reads were generated from the multi-chromosome genome of Agrobacterium tumefaciens strain LBA4404. Errors in the consensus two-directional (sense and antisense) “2D” sequences were first characterized by way of comparison with an internal reference assembly. Both Illumina-based correction and self-correction were performed and the resulting corrected reads assembled into high-quality hybrid and non-hybrid assemblies. Corrected read datasets and assemblies were subsequently compared. The results shown here indicate that both hybrid and non-hybrid methods can be used to assemble Oxford Nanopore reads into informative multi-chromosome assemblies, each with slightly different outcomes in terms of contiguity and accuracy.

Exploring structural variation and gene family architecture with De Novo assemblies of 15 Medicago genomes

BackgroundPrevious studies exploring sequence variation in the model legume, Medicago truncatula, relied on mapping short reads to a single reference. However, read-mapping approaches are inadequate to examine large, diverse gene families or to probe variation in repeat-rich or highly divergent genome regions. De novo sequencing and assembly of M. truncatula genomes enables near-comprehensive discovery of structural variants (SVs), analysis of rapidly evolving gene families, and ultimately, construction of a pan-genome.ResultsGenome-wide synteny based on 15 de novo M. truncatula assemblies effectively detected different types of SVs indicating that as much as 22% of the genome is involved in large structural changes, altogether affecting 28% of gene models. A total of 63 million base pairs (Mbp) of novel sequence was discovered, expanding the reference genome space for Medicago by 16%. Pan-genome analysis revealed that 42% (180 Mbp) of genomic sequences is missing in one or more accession, while examination of de novo annotated genes identified 67% (50,700) of all ortholog groups as dispensable – estimates comparable to recent studies in rice, maize and soybean. Rapidly evolving gene families typically associated with biotic interactions and stress response were found to be enriched in the accession-specific gene pool. The nucleotide-binding site leucine-rich repeat (NBS-LRR) family, in particular, harbors the highest level of nucleotide diversity, large effect single nucleotide change, protein diversity, and presence/absence variation. However, the leucine-rich repeat (LRR) and heat shock gene families are disproportionately affected by large effect single nucleotide changes and even higher levels of copy number variation.ConclusionsAnalysis of multiple M. truncatula genomes illustrates the value of de novo assemblies to discover and describe structural variation, something that is often under-estimated when using read-mapping approaches. Comparisons among the de novo assemblies also indicate that different large gene families differ in the architecture of their structural variation.

Sequencing-based large-scale genomics approaches with small numbers of isolated maize meiocytes.

High-throughput sequencing has become the large-scale approach of choice to study global gene expression and the distribution of specific chromatin marks and features. However, the limited availability of large amounts of purified cells made it very challenging to apply sequencing-based techniques in plant meiosis research in the past. In this paper, we describe a method to isolate meiocytes from maize anthers and detailed protocols to successfully perform RNA-seq, smRNA-seq, H3K4me3-ChIP-seq, and DNA bisulfite conversion sequencing with 5000–30,000 isolated maize male meiotic cells. These methods can be adjusted for other flowering plant species as well.

Repertoire of virus-derived small RNAs produced by mosquito and mammalian cells in response to dengue virus infection

RNA interference (RNAi) is the major defense of many arthropods against arthropod-borne RNA viruses (arboviruses), but the role of RNAi in vertebrate immunity to arboviruses is not clear. RNA viruses can trigger RNAi in vertebrate cells, but the vertebrate interferon response may obscure this interaction. We quantified virus-derived small RNAs (vRNAs) generated by mosquito (U4.4) cells and interferon-deficient (Vero) and interferon-competent (HuH-7) mammalian cells infected with a single isolate of mosquito-borne dengue virus. Mosquito cells produced significantly more vRNAs than mammalian cells, and mosquito cell vRNAs were derived from both the positive- and negative-sense dengue genomes whereas mammalian cell vRNAs were derived primarily from positive-sense genome. Mosquito cell vRNAs were predominantly 21 nucleotides in length whereas mammalian cell vRNAs were between 12 and 36 nucleotides with a modest peak at 24 nucleotides. Hot-spots, regions of the virus genome that generated a disproportionate number of vRNAs, overlapped among the cell lines.

The Complete Chloroplast Genome Sequences for Four Amaranthus Species (Amaranthaceae)

Premise of the study: The amaranth genus contains many important grain and weedy species. We further our understanding of the genus through the development of a complete reference chloroplast genome. Methods and Results: A high-quality Amaranthus hypochondriacus (Amaranthaceae) chloroplast genome assembly was developed using long-read technology. This reference genome was used to reconstruct the chloroplast genomes for two closely related grain species (A. cruentus and A. caudatus) and their putative progenitor (A. hybridus). The reference genome was 150,518 bp and possesses a circular structure of two inverted repeats (24,352 bp) separated by small (17,941 bp) and large (83,873 bp) single-copy regions; it encodes 111 genes, 72 for proteins. Relative to the reference chloroplast genome, an average of 210 single-nucleotide polymorphisms (SNPs) and 122 insertion/deletion polymorphisms (indels) were identified across the analyzed genomes. Conclusions: This reference chloroplast genome, along with the reported simple sequence repeats, SNPs, and indels, is an invaluable genetic resource for studying the phylogeny and genetic diversity within the amaranth genus.