scholarly journals Variant Phasing and Haplotypic Expression from Single-molecule Long-read Sequencing in Maize

2019 ◽  
Author(s):  
Bo Wang ◽  
Elizabeth Tseng ◽  
Primo Baybayan ◽  
Kevin Eng ◽  
Michael Regulski ◽  
...  

AbstractHaplotype phasing of genetic variants in maize is important for interpretation of the genome, population genetic analysis and functional genomic analysis of allelic activity. Accordingly, accurate methods for phasing the full-length isoforms are essential for functional genomics studies. We performed an isoform-level phasing study in maize, using two inbred lines and their reciprocal crosses, based on the single-molecule full-length cDNA sequencing. To phase and analyze the full-length transcripts between hybrids and parents, we developed a tool called IsoPhase. Using this tool, we validated the majority of SNPs called against matching short-read data and identified cases of allele-specific, gene-level and isoform-level expression. Our results revealed that maize parental lines and hybrid lines exhibit different splicing activities. After phasing 6,907 genes in two reciprocal hybrids using embryo, endosperm and root tissues, we annotated the SNPs and identified large-effect genes. In addition, based on single-molecule sequencing, we identified parent-of-origin isoforms in maize hybrids, distinct novel isoforms in maize parent and hybrid lines, and imprinted genes from different tissues. Finally, we characterized variation in cis- and trans-regulatory effects. Our study provides measures of haplotypic expression that could increase accuracy in studies of allelic expression.


Author(s):  
Brian M Forde ◽  
Andrew Henderson ◽  
Elliott G Playford ◽  
David Looke ◽  
Belinda C Henderson ◽  
...  

Abstract Background Diphtheria is a potentially fatal respiratory disease caused by toxigenic Corynebacterium diphtheriae. Although resistance to erythromycin has been recognized, β-lactam resistance in toxigenic diphtheria has not been described. Here, we report a case of fatal respiratory diphtheria caused by toxigenic C. diphtheriae resistant to penicillin and all other β-lactam antibiotics, and describe a novel mechanism of inducible carbapenem resistance associated with the acquisition of a mobile resistance element. Methods Long-read whole-genome sequencing was performed using Pacific Biosciences Single Molecule Real-Time sequencing to determine the genome sequence of C. diphtheriae BQ11 and the mechanism of β-lactam resistance. To investigate the phenotypic inducibility of meropenem resistance, short-read sequencing was performed using an Illumina NextSeq500 sequencer on the strain both with and without exposure to meropenem. Results BQ11 demonstrated high-level resistance to penicillin (benzylpenicillin minimum inhibitory concentration [MIC] ≥ 256 μg/ml), β-lactam/β-lactamase inhibitors and cephalosporins (amoxicillin/clavulanic acid MIC ≥ 256 μg/mL; ceftriaxone MIC ≥ 8 μg/L). Genomic analysis of BQ11 identified acquisition of a novel transposon carrying the penicillin-binding protein (PBP) Pbp2c, responsible for resistance to penicillin and cephalosporins. When strain BQ11 was exposed to meropenem, selective pressure drove amplification of the transposon in a tandem array and led to a corresponding change from a low-level to a high-level meropenem-resistant phenotype. Conclusions We have identified a novel mechanism of inducible antibiotic resistance whereby isolates that appear to be carbapenem susceptible on initial testing can develop in vivo resistance to carbapenems with repeated exposure. This phenomenon could have significant implications for the treatment of C. diphtheriae infection, and may lead to clinical failure.



2020 ◽  
Author(s):  
Nil Aygün ◽  
Angela L. Elwell ◽  
Dan Liang ◽  
Michael J. Lafferty ◽  
Kerry E. Cheek ◽  
...  

SummaryInterpretation of the function of non-coding risk loci for neuropsychiatric disorders and brain-relevant traits via gene expression and alternative splicing is mainly performed in bulk post-mortem adult tissue. However, genetic risk loci are enriched in regulatory elements of cells present during neocortical differentiation, and regulatory effects of risk variants may be masked by heterogeneity in bulk tissue. Here, we map e/sQTLs and allele specific expression in primary human neural progenitors (n=85) and their sorted neuronal progeny (n=74). Using colocalization and TWAS, we uncover cell-type specific regulatory mechanisms underlying risk for these traits.



2021 ◽  
Vol 12 ◽  
Author(s):  
Fiza Liaquat ◽  
Muhammad Farooq Hussain Munis ◽  
Samiah Arif ◽  
Urooj Haroon ◽  
Jianxin Shi ◽  
...  

Schima superba (Theaceae) is a subtropical evergreen tree and is used widely for forest firebreaks and gardening. It is a plant that tolerates salt and typically accumulates elevated amounts of manganese in the leaves. With large ecological amplitude, this tree species grows quickly. Due to its substantial biomass, it has a great potential for soil remediation. To evaluate the thorough framework of the mRNA, we employed PacBio sequencing technology for the first time to generate S. Superba transcriptome. In this analysis, overall, 511,759 full length non-chimeric reads were acquired, and 163,834 high-quality full-length reads were obtained. Overall, 93,362 open reading frames were obtained, of which 78,255 were complete. In gene annotation analyses, the Kyoto Encyclopedia of Genes and Genomes (KEGG), Clusters of Orthologous Genes (COG), Gene Ontology (GO), and Non-Redundant (Nr) databases were allocated 91,082, 71,839, 38,914, and 38,376 transcripts, respectively. To identify long non-coding RNAs (lncRNAs), we utilized four computational methods associated with protein families (Pfam), Cooperative Data Classification (CPC), Coding Assessing Potential Tool (CPAT), and Coding Non-Coding Index (CNCI) databases and observed 8,551, 9,174, 20,720, and 18,669 lncRNAs, respectively. Moreover, nine genes were randomly selected for the expression analysis, which showed the highest expression of Gene 6 (Na_Ca_ex gene), and CAX (CAX-interacting protein 4) was higher in manganese (Mn)-treated group. This work provided significant number of full-length transcripts and refined the annotation of the reference genome, which will ease advanced genetic analyses of S. superba.



2019 ◽  
Vol 20 (24) ◽  
pp. 6350 ◽  
Author(s):  
Nan Deng ◽  
Chen Hou ◽  
Fengfeng Ma ◽  
Caixia Liu ◽  
Yuxin Tian

The limitations of RNA sequencing make it difficult to accurately predict alternative splicing (AS) and alternative polyadenylation (APA) events and long non-coding RNAs (lncRNAs), all of which reveal transcriptomic diversity and the complexity of gene regulation. Gnetum, a genus with ambiguous phylogenetic placement in seed plants, has a distinct stomatal structure and photosynthetic characteristics. In this study, a full-length transcriptome of Gnetum luofuense leaves at different developmental stages was sequenced with the latest PacBio Sequel platform. After correction by short reads generated by Illumina RNA-Seq, 80,496 full-length transcripts were obtained, of which 5269 reads were identified as isoforms of novel genes. Additionally, 1660 lncRNAs and 12,998 AS events were detected. In total, 5647 genes in the G. luofuense leaves had APA featured by at least one poly(A) site. Moreover, 67 and 30 genes from the bHLH gene family, which play an important role in stomatal development and photosynthesis, were identified from the G. luofuense genome and leaf transcripts, respectively. This leaf transcriptome supplements the reference genome of G. luofuense, and the AS events and lncRNAs detected provide valuable resources for future studies of investigating low photosynthetic capacity of Gnetum.



Forests ◽  
2020 ◽  
Vol 11 (8) ◽  
pp. 866
Author(s):  
Lei Kan ◽  
Qicong Liao ◽  
Zhiyao Su ◽  
Yushan Tan ◽  
Shuyu Wang ◽  
...  

Madhuca pasquieri (Dubard) Lam. is a tree on the International Union for Conservation of Nature Red List and a national key protected wild plant (II) of China, known for its seed oil and timber. However, lacking of genomic and transcriptome data for this species hampers study of its reproduction, utilization, and conservation. Here, single-molecule long-read sequencing (PacBio) and next-generation sequencing (Illumina) were combined to obtain the transcriptome from five developmental stages of M. pasquieri. Overall, 25,339 transcript isoforms were detected by PacBio, including 24,492 coding sequences (CDSs), 9440 simple sequence repeats (SSRs), 149 long non-coding RNAs (lncRNAs), and 182 alternative splicing (AS) events, a majority was retained intron (RI). A further 1058 transcripts were identified as transcriptional factors (TFs) from 51 TF families. PacBio recovered more full-length transcript isoforms with a longer length, and a higher expression level, whereas larger number of transcripts (124,405) was captured in de novo from Illumina. Using Nr, Swissprot, KOG, and KEGG databases, 24,405 transcripts (96.31%) were annotated by PacBio. Functional annotation revealed a role for the auxin, abscisic acid, gibberellin, and cytokinine metabolic pathways in seed germination and post-germination. These findings support further studies on seed germination mechanism and genome of M. pasquieri, and better protection of this endangered species.



2021 ◽  
Author(s):  
Dafni A Glinos ◽  
Garrett Garborcauskas ◽  
Paul Hoffman ◽  
Nava Ehsan ◽  
Lihua Jiang ◽  
...  

SummaryRegulation of transcript structure generates transcript diversity and plays an important role in human disease. The advent of long-read sequencing technologies offers the opportunity to study the role of genetic variation in transcript structure. In this paper, we present a large human long-read RNA-seq dataset using the Oxford Nanopore Technologies platform from 88 samples from GTEx tissues and cell lines, complementing the GTEx resource. We identified just under 100,000 new transcripts for annotated genes, and validated the protein expression of a similar proportion of novel and annotated transcripts. We developed a new computational package, LORALS, to analyze genetic effects of rare and common variants on the transcriptome via allele-specific analysis of long reads. We called allele-specific expression and transcript structure events, providing novel insights into the specific transcript alterations caused by common and rare genetic variants and highlighting the resolution gained from long-read data. We were able to perturb transcript structure upon knockdown of PTBP1, an RNA binding protein that mediates splicing, thereby finding genetic regulatory effects that are modified by the cellular environment. Finally, we use this dataset to enhance variant interpretation and study rare variants leading to aberrant splicing patterns.



2017 ◽  
Author(s):  
Brian M. Forde ◽  
Lauren J. McAllister ◽  
James C. Paton ◽  
Adrienne W. Paton ◽  
Scott A. Beatson

AbstractShiga toxigenicEscherichia coli(STEC) are important food-borne pathogens and a major cause of haemorrhagic colitis and haemolytic-uremic syndrome (HUS) worldwide. In 1995 a severe HUS outbreak in Adelaide occurred. A recent genomic analysis of STEC O111:H-strains 95JB1 and 95NR1 from this outbreak found that the more virulent isolate, 95NR1, harboured two additional copies of the Shiga toxin 2 (Stx2) genes although the structure of the Stx2-converting prophages could not be fully resolved due to the fragmented assembly. In this study we have used Pacific Biosciences (PacBio) single molecule real-time (SMRT) long read sequencing to characterise the complete epigenome (genome and methylome) of 95JB1 and 95NR1. Using long reads we completely resolved the structure of two, tandemly inserted, stx2-converting phage in 95NR1. Our analysis of the methylome of 95NR1 and 95JB1 identified hemi-methylation of a novel motif (5’-CTGCm6AG-3’) in more than 4000 sites in the 95NR1 genome. These sites were entirely unmethalyted in the 95JB1, including at least 180 potential promoter regions that could explain regulatory differences between the strains. We identified a Type IIG methyltransferase encoded in both genomes in association with three additional genes in an operon-like arrangement. IS1203mediated disruption of this operon in 95JB1 is the likely cause of the observed differential patterns of methylation between 95NR1 and 95JB1. This study demonstrates the enormous potential of PacBio SMRT sequencing to resolve complex prophage regions and reveal the genetic and epigenetic heterogeneity within a clonal population of bacteria.



Author(s):  
Chengcai Zhang ◽  
Huadong Ren ◽  
Xiaohua Yao ◽  
Kailiang Wang ◽  
Jun Chang

Abstract Pecan is rich in bioactive components such as fatty acids and flavonoids and is an important nut type worldwide. Therefore, the molecular mechanisms of phytochemical biosynthesis in pecan are a focus of research. Recently, a draft genome and several transcriptomes have been published. However, the full-length mRNA transcripts remain unclear, and the regulatory mechanisms behind the quality components biosynthesis and accumulation have not been fully investigated. In this study, single-molecule long read sequencing technology was used to obtain full-length transcripts of pecan kernels. In total, 37 504 isoforms of 16 702 genes were mapped to the reference genome. The numbers of known isoforms, new isoforms, and novel isoforms were 9013 (24.03%), 26 080 (69.54%), and 2411 (6.51%), respectively. Over 80% of the transcripts (30 751, 81.99%) had functional annotations. A total of 15 465 alternative splicing (AS) events and 65 761 alternative polyadenylation events were detected; wherein, the retained intron was the predominant type (5652, 36.55%) of AS. Furthermore, 1894 long non-coding RNAs and 1643 transcription factors were predicted using bioinformatics methods. Finally, the structural genes associated with fatty acid (FA) and flavonoid biosynthesis were characterized. A high frequency of AS accuracy (70.31%) was observed in FA synthesis-associated genes. The present study provides a full-length transcriptome dataset of pecan kernels, which will significantly enhance the understanding of the regulatory basis of phytochemical biosynthesis during pecan kernel maturation.



2021 ◽  
pp. jmedgenet-2021-108147
Author(s):  
Florentine Scharf ◽  
Rafaela Magalhaes Leal Silva ◽  
Monika Morak ◽  
Alex Hastie ◽  
Julia M A Pickl ◽  
...  

PurposeApproximately 20% of patients with clinical familial adenomatous polyposis (FAP) remain unsolved after molecular genetic analysis of the APC and other polyposis genes, suggesting additional pathomechanisms.MethodsWe applied multidimensional genomic analysis employing chromosomal microarray profiling, optical mapping, long-read genome and RNA sequencing combined with FISH and standard PCR of genomic and complementary DNA to decode a patient with an attenuated FAP that had remained unsolved by Sanger sequencing and multigene panel next-generation sequencing for years.ResultsWe identified a complex 3.9 Mb rearrangement involving 14 fragments from chromosome 5q22.1q22.3 of which three were lost, 1 reinserted into chromosome 5 and 10 inserted into chromosome 10q21.3 in a seemingly random order and orientation thus fulfilling the major criteria of chromothripsis. The rearrangement separates APC promoter 1B from the coding ORF (open reading frame) thus leading to allele-specific downregulation of APC mRNA. The rearrangement also involves three additional genes implicated in the APC–Axin–GSK3B–β-catenin signalling pathway.ConclusionsBased on comprehensive genomic analysis, we propose that constitutional chromothripsis dampening APC expression, possibly modified by additional APC–Axin–GSK3B–β-catenin pathway disruptions, underlies the patient’s clinical phenotype. The combinatorial approach we deployed provides a powerful tool set for deciphering unsolved familial polyposis and potentially other tumour syndromes and monogenic diseases.



Sign in / Sign up

Export Citation Format

Share Document