1
|
Yu D, Zhou M, Chen W, Ding Z, Wang C, Qian Y, Liu Y, He S, Yang L. Characterization of transcriptome changes in saline stress adaptation on Leuciscus merzbacheri using PacBio Iso-Seq and RNA-Seq. DNA Res 2024; 31:dsae019. [PMID: 38807352 PMCID: PMC11161863 DOI: 10.1093/dnares/dsae019] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2024] [Revised: 05/22/2024] [Accepted: 05/28/2024] [Indexed: 05/30/2024] Open
Abstract
Leuciscus merzbacheri is a native fish species found exclusively in the Junggar Basin in Xinjiang. It exhibits remarkable adaptability, thriving in varying water conditions such as the saline waters, the semi-saline water, and the freshwater. Despite its significant economic and ecological value, the underlying mechanisms of its remarkable salinity tolerance remain elusive. Our study marks the first time the full-length transcriptome of L. merzbacheri has been reported, utilizing RNA-Seq and PacBio Iso-Seq technologies. We found that the average length of the full-length transcriptome is 1,780 bp, with an N50 length of 2,358 bp. We collected RNA-Seq data from gill, liver, and kidney tissues of L. merzbacheri from both saline water and freshwater environments and conducted comparative analyses across these tissues. Further analysis revealed significant enrichment in several key functional gene categories and signalling pathways related to stress response and environmental adaptation. The findings provide a valuable genetic resource for further investigation into saline-responsive candidate genes, which will deepen our understanding of teleost adaptation to extreme environmental stress. This knowledge is crucial for the future breeding and conservation of native fish species.
Collapse
Affiliation(s)
- Dan Yu
- School of Ecology and Environment, Tibet University, Lhasa, 850000, China
- State Key Laboratory of Freshwater Ecology and Biotechnology, Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan, 430072, China
| | - Min Zhou
- School of Life Sciences, Jianghan Universily, Wuhan 430056, China
| | - Wenjun Chen
- State Key Laboratory of Freshwater Ecology and Biotechnology, Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan, 430072, China
- College of Advanced Agricultural Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
| | - Zufa Ding
- State Key Laboratory of Freshwater Ecology and Biotechnology, Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan, 430072, China
- College of Fisheries, Huazhong Agricultural University, Wuhan, 430070, China
| | - Cheng Wang
- State Key Laboratory of Freshwater Ecology and Biotechnology, Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan, 430072, China
- College of Advanced Agricultural Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
| | - Yuting Qian
- State Key Laboratory of Freshwater Ecology and Biotechnology, Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan, 430072, China
- College of Advanced Agricultural Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
| | - Yang Liu
- State Key Laboratory of Freshwater Ecology and Biotechnology, Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan, 430072, China
- College of Advanced Agricultural Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
| | - Shunping He
- State Key Laboratory of Freshwater Ecology and Biotechnology, Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan, 430072, China
| | - Liandong Yang
- State Key Laboratory of Freshwater Ecology and Biotechnology, Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan, 430072, China
| |
Collapse
|
2
|
Liu X, Wu J, Li M, Zuo F, Zhang G. A Comparative Full-Length Transcriptome Analysis Using Oxford Nanopore Technologies (ONT) in Four Tissues of Bovine Origin. Animals (Basel) 2024; 14:1646. [PMID: 38891695 PMCID: PMC11170998 DOI: 10.3390/ani14111646] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/10/2024] [Revised: 05/28/2024] [Accepted: 05/28/2024] [Indexed: 06/21/2024] Open
Abstract
The transcriptome complexity and splicing patterns in male and female cattle are ambiguous, presenting a substantial obstacle to genomic selection programs that seek to improve productivity, disease resistance, and reproduction in cattle. A comparative transcriptomic analysis using Oxford Nanopore Technologies (ONT) was conducted in bovine testes (TESTs), ovaries (OVAs), muscles (MUSCs), and livers (LIVs). An average of 5,144,769 full-length reads were obtained from each sample. The TESTs were found to have the greatest number of alternative polyadenylation (APA) events involved in processes such as sperm flagellum development and fertilization in male reproduction. In total, 438 differentially expressed transcripts (DETs) were identified in the LIVs in a comparison of females vs. males, and 214 DETs were identified in the MUSCs between females and males. Additionally, 14,735, 36,347, and 33,885 DETs were detected in MUSC vs. LIV, MUSC vs. TEST, and OVA vs. TEST comparisons, respectively, revealing the complexity of the TEST. Gene Set Enrichment Analysis (GSEA) showed that these DETs were mainly involved in the "spermatogenesis", "flagellated sperm motility", "spermatid development", "reproduction", "reproductive process", and "microtubule-based movement" KEGG pathways. Additional studies are necessary to further characterize the transcriptome in different cell types, developmental stages, and physiological conditions in bovines and ascertain the functions of the novel transcripts.
Collapse
Affiliation(s)
- Xinyue Liu
- College of Animal Science and Technology, Southwest University, Rongchang, Chongqing 402460, China; (X.L.); (J.W.); (M.L.); (F.Z.)
| | - Jiaxin Wu
- College of Animal Science and Technology, Southwest University, Rongchang, Chongqing 402460, China; (X.L.); (J.W.); (M.L.); (F.Z.)
| | - Meichen Li
- College of Animal Science and Technology, Southwest University, Rongchang, Chongqing 402460, China; (X.L.); (J.W.); (M.L.); (F.Z.)
| | - Fuyuan Zuo
- College of Animal Science and Technology, Southwest University, Rongchang, Chongqing 402460, China; (X.L.); (J.W.); (M.L.); (F.Z.)
- Beef Cattle Engineering and Technology Research Center of Chongqing, Southwest University, Rongchang, Chongqing 402460, China
| | - Gongwei Zhang
- College of Animal Science and Technology, Southwest University, Rongchang, Chongqing 402460, China; (X.L.); (J.W.); (M.L.); (F.Z.)
- Beef Cattle Engineering and Technology Research Center of Chongqing, Southwest University, Rongchang, Chongqing 402460, China
| |
Collapse
|
3
|
Lou F, Ren Z, Tang Y, Han Z. Full-length transcriptome reveals the circularly polarized light response-related molecular genetic characteristics of Oratosquilla oratoria. COMPARATIVE BIOCHEMISTRY AND PHYSIOLOGY. PART D, GENOMICS & PROTEOMICS 2024; 49:101183. [PMID: 38141370 DOI: 10.1016/j.cbd.2023.101183] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/25/2023] [Revised: 12/16/2023] [Accepted: 12/16/2023] [Indexed: 12/25/2023]
Abstract
The mantis shrimp is the only animal that can recognize circularly polarized light (CPL), but its molecular genetic characteristics are unclear. Multi-tissue level full-length (FL) transcriptome sequencing of Oratosquilla oratoria, a representative widely distributed mantis shrimp, was performed in the present study. We used comparative transcriptomics to explore the critical genes of O. oratoria selected by CPL and the GNβ gene associated with CPL signal transduction was hypothesized to be positively selected. Furthermore, the FL transcriptomes of O. oratoria compound eyes under five light conditions were sequenced and used to detect alternative splicing (AS). The ASs associated with CPL recognition mainly occurred in the LWS, ARR and TRPC regions. The number of FL transcripts with AS events and annotation information also provided evidence that O. oratoria could recognize LCPL. Additionally, 51 sequences belonging to the LWS, UV and Peropsin gene families were identified based on conserved 7tm domains. The LWS, UV and Peropsin opsins have similar 3D structures with seven domains across the cell membrane and conserved KSLRTPSN, DRY, and QAKK motifs. In conclusion, these results are undoubtedly valuable for perfecting the vision theory of O. oratoria and other mantis shrimp.
Collapse
Affiliation(s)
- Fangrui Lou
- School of Ocean, Yantai University, Yantai 264003, Shandong, China.
| | - Zhongjie Ren
- School of Ocean, Yantai University, Yantai 264003, Shandong, China
| | - Yongzheng Tang
- School of Ocean, Yantai University, Yantai 264003, Shandong, China
| | - Zhiqiang Han
- Fishery College, Zhejiang Ocean University, Zhoushan 316022, Zhejiang, China.
| |
Collapse
|
4
|
Zhao P, Gu L, Gao Y, Pan Z, Liu L, Li X, Zhou H, Yu D, Han X, Qian L, Liu GE, Fang L, Wang Z. Young SINEs in pig genomes impact gene regulation, genetic diversity, and complex traits. Commun Biol 2023; 6:894. [PMID: 37652983 PMCID: PMC10471783 DOI: 10.1038/s42003-023-05234-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2022] [Accepted: 08/09/2023] [Indexed: 09/02/2023] Open
Abstract
Transposable elements (TEs) are a major source of genetic polymorphisms and play a role in chromatin architecture, gene regulatory networks, and genomic evolution. However, their functional role in pigs and contributions to complex traits are largely unknown. We created a catalog of TEs (n = 3,087,929) in pigs and found that young SINEs were predominantly silenced by histone modifications, DNA methylation, and decreased accessibility. However, some transcripts from active young SINEs showed high tissue-specificity, as confirmed by analyzing 3570 RNA-seq samples. We also detected 211,067 dimorphic SINEs in 374 individuals, including 340 population-specific ones associated with local adaptation. Mapping these dimorphic SINEs to genome-wide associations of 97 complex traits in pigs, we found 54 candidate genes (e.g., ANK2 and VRTN) that might be mediated by TEs. Our findings highlight the important roles of young SINEs and provide a supplement for genotype-to-phenotype associations and modern breeding in pigs.
Collapse
Affiliation(s)
- Pengju Zhao
- Hainan Institute, Zhejiang University, Yongyou Industry Park, Yazhou Bay Sci-Tech City, Sanya, 572000, China
- College of Animal Sciences, Zhejiang University, Hangzhou, Zhejiang, 310058, China
| | - Lihong Gu
- Institute of Animal Science & Veterinary Medicine, Hainan Academy of Agricultural Sciences, No. 14 Xingdan Road, Haikou, 571100, China
| | - Yahui Gao
- Animal Genomics and Improvement Laboratory, Beltsville Agricultural Research Center, Agricultural Research Service, USDA, Beltsville, MD, 20705, USA
| | - Zhangyuan Pan
- Department of Animal Science, University of California, Davis, CA, 95616, USA
| | - Lei Liu
- Shenzhen Branch, Guangdong Laboratory of Lingnan Modern Agriculture, Genome Analysis Laboratory of the Ministry of Agriculture and Rural Affairs, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, 518124, China
| | - Xingzheng Li
- Shenzhen Branch, Guangdong Laboratory of Lingnan Modern Agriculture, Genome Analysis Laboratory of the Ministry of Agriculture and Rural Affairs, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, 518124, China
| | - Huaijun Zhou
- Department of Animal Science, University of California, Davis, CA, 95616, USA
| | - Dongyou Yu
- Hainan Institute, Zhejiang University, Yongyou Industry Park, Yazhou Bay Sci-Tech City, Sanya, 572000, China
- College of Animal Sciences, Zhejiang University, Hangzhou, Zhejiang, 310058, China
| | - Xinyan Han
- Hainan Institute, Zhejiang University, Yongyou Industry Park, Yazhou Bay Sci-Tech City, Sanya, 572000, China
- College of Animal Sciences, Zhejiang University, Hangzhou, Zhejiang, 310058, China
| | - Lichun Qian
- Hainan Institute, Zhejiang University, Yongyou Industry Park, Yazhou Bay Sci-Tech City, Sanya, 572000, China
- College of Animal Sciences, Zhejiang University, Hangzhou, Zhejiang, 310058, China
| | - George E Liu
- Animal Genomics and Improvement Laboratory, Beltsville Agricultural Research Center, Agricultural Research Service, USDA, Beltsville, MD, 20705, USA.
| | - Lingzhao Fang
- Center for Quantitative Genetics and Genomics, Aarhus University, Aarhus, 8000, Denmark.
| | - Zhengguang Wang
- Hainan Institute, Zhejiang University, Yongyou Industry Park, Yazhou Bay Sci-Tech City, Sanya, 572000, China.
- College of Animal Sciences, Zhejiang University, Hangzhou, Zhejiang, 310058, China.
| |
Collapse
|
5
|
Wang P, Zhang X, Huo H, Li W, Liu Z, Wang L, Li L, Sun YH, Huo J. Transcriptomic analysis of testis and epididymis tissues from Banna mini-pig inbred line boars with single-molecule long-read sequencing†. Biol Reprod 2023; 108:465-478. [PMID: 36477198 DOI: 10.1093/biolre/ioac216] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2022] [Revised: 05/04/2022] [Accepted: 11/30/2022] [Indexed: 12/14/2022] Open
Abstract
In mammals, testis and epididymis are critical components of the male reproductive system for androgen production, spermatogenesis, sperm transportation, as well as sperm maturation. Here, we report single-molecule real-time sequencing data from the testis and epididymis of the Banna mini-pig inbred line (BMI), a promising laboratory animal for medical research. We obtained high-quality full-length transcriptomes and identified 9879 isoforms and 8761 isoforms in the BMI testis and epididymis, respectively. Most of the isoforms we identified have novel exon structures that will greatly improve the annotation of testis- and epididymis-expressed genes in pigs. We also found that 3055 genes (over 50%) were shared between BMI testis and epididymis, indicating widespread expression profiles of genes related to reproduction. We characterized extensive alternative splicing events in BMI testis and epididymis and showed that 96 testis-expressed genes and 79 epididymis-expressed genes have more than six isoforms, revealing the complexity of alternative splicing. We accurately defined the transcribed isoforms in BMI testis and epididymis by combining Pacific Biotechnology Isoform-sequencing (PacBio Iso-Seq) and Illumina RNA Sequencing (RNA-seq) techniques. The refined annotation of some key genes governing male reproduction will facilitate further understanding of the molecular mechanisms underlying BMI male sterility. In addition, the high-confident identification of 548 and 669 long noncoding RNAs (lncRNAs) in these two tissues has established a candidate gene set for future functional investigations. Overall, our study provides new insights into the role of the testis and epididymis during BMI reproduction, paving the path for further studies on BMI male infertility.
Collapse
Affiliation(s)
- Pei Wang
- College of Animal Science and Technology, Yunnan Agricultural University, Kunming, China
| | - Xia Zhang
- College of Life Science, Lyuliang University, Lvliang, China
| | - Hailong Huo
- Yunnan Vocational and Technical college of Agriculture, Kunming, China
| | - Weizhen Li
- College of Veterinary Medicine, Yunnan Agricultural University, Kunming, China
| | - Zhipeng Liu
- College of Animal Science and Technology, Yunnan Agricultural University, Kunming, China
| | - Lina Wang
- College of Animal Science and Technology, Yunnan Agricultural University, Kunming, China
| | - Luogang Li
- College of Animal Science and Technology, Yunnan Agricultural University, Kunming, China
| | - Yu H Sun
- Department of Biology, University of Rochester, Rochester, NY, USA
| | - Jinlong Huo
- College of Animal Science and Technology, Yunnan Agricultural University, Kunming, China
- Department of Biology, University of Rochester, Rochester, NY, USA
| |
Collapse
|
6
|
Li J, Guan D, Halstead MM, Islas-Trejo AD, Goszczynski DE, Ernst CW, Cheng H, Ross P, Zhou H. Transcriptome annotation of 17 porcine tissues using nanopore sequencing technology. Anim Genet 2023; 54:35-44. [PMID: 36385508 DOI: 10.1111/age.13274] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/28/2022] [Revised: 10/20/2022] [Accepted: 11/01/2022] [Indexed: 11/18/2022]
Abstract
The annotation of animal genomes plays an important role in elucidating molecular mechanisms behind the genetic control of economically important traits. Here, we employed long-read sequencing technology, Oxford Nanopore Technology, to annotate the pig transcriptome across 17 tissues from two Yorkshire littermate pigs. More than 9.8 million reads were obtained from a single flow cell, and 69 781 unique transcripts at 50 108 loci were identified. Of these transcripts, 16 255 were found to be novel isoforms, and 22 344 were found at loci that were novel and unannotated in the Ensembl (release 102) and NCBI (release 106) annotations. Novel transcripts were mostly expressed in cerebellum, followed by lung, liver, spleen, and hypothalamus. By comparing the unannotated transcripts to existing databases, there were 21 285 (95.3%) transcripts matched to the NT database (v5) and 13 676 (61.2%) matched to the NR database (v5). Moreover, there were 4324 (19.4%) transcripts matched to the SwissProt database (v5), corresponding to 11 356 proteins. Tissue-specific gene expression analyses showed that 9749 transcripts were highly tissue-specific, and cerebellum contained the most tissue-specific transcripts. As the same samples were used for the annotation of cis-regulatory elements in the pig genome, the transcriptome annotation generated by this study provides an additional and complementary annotation resource for the Functional Annotation of Animal Genomes effort to comprehensively annotate the pig genome.
Collapse
Affiliation(s)
- Jinghui Li
- Department of Animal Science, University of California Davis, Davis, California, USA
| | - Dailu Guan
- Department of Animal Science, University of California Davis, Davis, California, USA
| | - Michelle M Halstead
- Department of Animal Science, University of California Davis, Davis, California, USA
| | - Alma D Islas-Trejo
- Department of Animal Science, University of California Davis, Davis, California, USA
| | - Daniel E Goszczynski
- Department of Animal Science, University of California Davis, Davis, California, USA
| | - Catherine W Ernst
- Department of Animal Science, Michigan State University, East Lansing, Michigan, USA
| | - Hao Cheng
- Department of Animal Science, University of California Davis, Davis, California, USA
| | - Pablo Ross
- Department of Animal Science, University of California Davis, Davis, California, USA
| | - Huaijun Zhou
- Department of Animal Science, University of California Davis, Davis, California, USA
| |
Collapse
|
7
|
Chen S, Liu S, Shi S, Jiang Y, Cao M, Tang Y, Li W, Liu J, Fang L, Yu Y, Zhang S. Comparative epigenomics reveals the impact of ruminant-specific regulatory elements on complex traits. BMC Biol 2022; 20:273. [PMID: 36482458 PMCID: PMC9730597 DOI: 10.1186/s12915-022-01459-0] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/25/2022] [Accepted: 11/07/2022] [Indexed: 12/13/2022] Open
Abstract
BACKGROUND Insights into the genetic basis of complex traits and disease in both human and livestock species have been achieved over the past decade through detection of genetic variants in genome-wide association studies (GWAS). A majority of such variants were found located in noncoding genomic regions, and though the involvement of numerous regulatory elements (REs) has been predicted across multiple tissues in domesticated animals, their evolutionary conservation and effects on complex traits have not been fully elucidated, particularly in ruminants. Here, we systematically analyzed 137 epigenomic and transcriptomic datasets of six mammals, including cattle, sheep, goats, pigs, mice, and humans, and then integrated them with large-scale GWAS of complex traits. RESULTS Using 40 ChIP-seq datasets of H3K4me3 and H3K27ac, we detected 68,479, 58,562, 63,273, 97,244, 111,881, and 87,049 REs in the liver of cattle, sheep, goats, pigs, humans and mice, respectively. We then systematically characterized the dynamic functional landscapes of these REs by integrating multi-omics datasets, including gene expression, chromatin accessibility, and DNA methylation. We identified a core set (n = 6359) of ruminant-specific REs that are involved in liver development, metabolism, and immune processes. Genes with more complex cis-REs exhibited higher gene expression levels and stronger conservation across species. Furthermore, we integrated expression quantitative trait loci (eQTLs) and GWAS from 44 and 52 complex traits/diseases in cattle and humans, respectively. These results demonstrated that REs with different degrees of evolutionary conservation across species exhibited distinct enrichments for GWAS signals of complex traits. CONCLUSIONS We systematically annotated genome-wide functional REs in liver across six mammals and demonstrated the evolution of REs and their associations with transcriptional output and conservation. Detecting lineage-specific REs allows us to decipher the evolutionary and genetic basis of complex phenotypes in livestock and humans, which may benefit the discovery of potential biomedical models for functional variants and genes of specific human diseases.
Collapse
Affiliation(s)
- Siqian Chen
- grid.22935.3f0000 0004 0530 8290Key Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture and Rural Affairs & National Engineering Laboratory for Animal Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, China
| | - Shuli Liu
- grid.22935.3f0000 0004 0530 8290Key Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture and Rural Affairs & National Engineering Laboratory for Animal Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, China ,grid.494629.40000 0004 8008 9315 School of Life Sciences, Westlake University, Hangzhou, China
| | - Shaolei Shi
- grid.22935.3f0000 0004 0530 8290Key Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture and Rural Affairs & National Engineering Laboratory for Animal Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, China
| | - Yifan Jiang
- grid.22935.3f0000 0004 0530 8290Key Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture and Rural Affairs & National Engineering Laboratory for Animal Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, China
| | - Mingyue Cao
- grid.22935.3f0000 0004 0530 8290Key Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture and Rural Affairs & National Engineering Laboratory for Animal Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, China
| | - Yongjie Tang
- grid.22935.3f0000 0004 0530 8290Key Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture and Rural Affairs & National Engineering Laboratory for Animal Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, China
| | - Wenlong Li
- grid.22935.3f0000 0004 0530 8290Key Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture and Rural Affairs & National Engineering Laboratory for Animal Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, China
| | - Jianfeng Liu
- grid.22935.3f0000 0004 0530 8290Key Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture and Rural Affairs & National Engineering Laboratory for Animal Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, China
| | - Lingzhao Fang
- grid.4305.20000 0004 1936 7988MRC Human Genetics Unit at the Institute of Genetics and Cancer, University of Edinburgh, Edinburgh, UK ,grid.7048.b0000 0001 1956 2722Center for Quantitative Genetics and Genomics (QGG), Aarhus University, Aarhus, Denmark
| | - Ying Yu
- grid.22935.3f0000 0004 0530 8290Key Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture and Rural Affairs & National Engineering Laboratory for Animal Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, China
| | - Shengli Zhang
- grid.22935.3f0000 0004 0530 8290Key Laboratory of Animal Genetics, Breeding and Reproduction, Ministry of Agriculture and Rural Affairs & National Engineering Laboratory for Animal Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, China
| |
Collapse
|
8
|
Hou L, Wang M, Zhu L, Ning M, Bi J, Du J, Kong X, Gu W, Meng Q. Full-length transcriptome sequencing and comparative transcriptome analysis of Eriocheir sinensis in response to infection by the microsporidian Hepatospora eriocheir. Front Cell Infect Microbiol 2022; 12:997574. [PMID: 36530442 PMCID: PMC9754153 DOI: 10.3389/fcimb.2022.997574] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/19/2022] [Accepted: 11/18/2022] [Indexed: 12/02/2022] Open
Abstract
As a new generation of high-throughput sequencing technology, PacBio Iso-Seq technology (Iso-Seq) provides a better alternative sequencing method for the acquisition of full-length unigenes. In this study, a total of 22.27 gigabyte (Gb) subread bases and 128,614 non-redundant unigenes (mean length: 2,324 bp) were obtained from six main tissues of Eriocheir sinensis including the heart, nerve, intestine, muscle, gills and hepatopancreas. In addition, 74,732 unigenes were mapped to at least one of the following databases: Non-Redundant Protein Sequence Database (NR), Gene Ontology (GO), Kyoto Encyclopaedia of Genes and Genomes (KEGG), KEGG Orthology (KO) and Protein family (Pfam). In addition, 6696 transcription factors (TFs), 28,458 long non-coding RNAs (lncRNAs) and 94,230 mRNA-miRNA pairs were identified. Hepatospora eriocheir is the primary pathogen of E. sinensis and can cause hepatopancreatic necrosis disease (HPND); the intestine is the main target tissue. Here, we attempted to identify the key genes related to H. eriocheir infection in the intestines of E. sinensis. By combining Iso-Seq and Illumina RNA-seq analysis, we identified a total of 12,708 differentially expressed unigenes (DEUs; 6,696 upregulated and 6,012 downregulated) in the crab intestine following infection with H. eriocheir. Based on the biological analysis of these DEUs, several key processes were identified, including energy metabolism-related pathways, cell apoptosis and innate immune-related pathways. Twelve selected genes from these DEUs were subsequently verified by quantitative real-time PCR (qRT-PCR) analysis. Our findings enhance our understanding of the E. sinensis transcriptome and the specific association between E. sinensis and H. eriocheir infection.
Collapse
Affiliation(s)
- Libo Hou
- Engineering Lab of Henan Province for Aquatic Animal Disease Control, College of Fisheries, Henan Normal University, Xinxiang, China
| | - Mengdi Wang
- Engineering Lab of Henan Province for Aquatic Animal Disease Control, College of Fisheries, Henan Normal University, Xinxiang, China
| | - Lei Zhu
- Engineering Lab of Henan Province for Aquatic Animal Disease Control, College of Fisheries, Henan Normal University, Xinxiang, China
| | - Mingxiao Ning
- Institution of Quality Standard and Testing Technology for Agro-product, Shandong Academy of Agricultural Science, Jinan, Shandong, China
| | - Jingxiu Bi
- Institution of Quality Standard and Testing Technology for Agro-product, Shandong Academy of Agricultural Science, Jinan, Shandong, China
| | - Jie Du
- Animal Husbandry and Veterinary College, Jiangsu Vocational College of Agriculture and Forestry, Jurong, Jiangsu, China
| | - Xianghui Kong
- Engineering Lab of Henan Province for Aquatic Animal Disease Control, College of Fisheries, Henan Normal University, Xinxiang, China
| | - Wei Gu
- Jiangsu Key Laboratory for Aquatic Crustacean Diseases, College of Marine Science and Engineering, Nanjing Normal University, Nanjing, Jiangsu, China
| | - Qingguo Meng
- Jiangsu Key Laboratory for Aquatic Crustacean Diseases, College of Marine Science and Engineering, Nanjing Normal University, Nanjing, Jiangsu, China,*Correspondence: Qingguo Meng,
| |
Collapse
|
9
|
Zhou K, Chen Z, Du X, Huang Y, Qin J, Wen L, Pan X, Lin Y. SMRT Sequencing Reveals Candidate Genes and Pathways With Medicinal Value in Cipangopaludina chinensis. Front Genet 2022; 13:881952. [PMID: 35783279 PMCID: PMC9243326 DOI: 10.3389/fgene.2022.881952] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/23/2022] [Accepted: 05/26/2022] [Indexed: 12/03/2022] Open
Abstract
Cipangopaludina chinensis is an economically important aquatic snail with high medicinal value. However, molecular biology research on C. chinensis is limited by the lack of a reference genome, so the analysis of its transcripts is an important step to study the regulatory genes of various substances in C. chinensis. Herein, we conducted the first full-length transcriptome analysis of C. chinensis using PacBio single-molecule real-time (SMRT) sequencing technology. We identified a total of 26,312 unigenes with an average length of 2,572 bp, of which the largest number of zf-c2h2 transcription factor families (120,18.24%) were found, and also observed that the majority of the 8,058 SSRs contained 4-7 repeat units, which provided data for subsequent work on snail genetics Subsequently, 91.86% (24,169) of the genes were successfully annotated to the four major databases, while the highest homology was observed with Pomacea canaliculata. Functional annotation revealed that the majority of transcripts were enriched in metabolism, signal transduction and Immune-related pathways, and several candidate genes involved in drug metabolism and immune response were identified (e.g., CYP1A1, CYP2J, CYP2U1, GST, ,PIK3, PDE3A, PRKAG). This study lays a foundation for future molecular biology research and provides a reference for studying genes associated with the medicinal value of C. chinensis.
Collapse
Affiliation(s)
| | | | | | | | | | | | | | - Yong Lin
- *Correspondence: Xianhui Pan, ; Yong Lin,
| |
Collapse
|
10
|
Tombácz D, Kakuk B, Torma G, Csabai Z, Gulyás G, Tamás V, Zádori Z, Jefferson VA, Meyer F, Boldogkői Z. In-Depth Temporal Transcriptome Profiling of an Alphaherpesvirus Using Nanopore Sequencing. Viruses 2022; 14:v14061289. [PMID: 35746760 PMCID: PMC9229804 DOI: 10.3390/v14061289] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/06/2022] [Revised: 06/05/2022] [Accepted: 06/08/2022] [Indexed: 12/10/2022] Open
Abstract
In this work, a long-read sequencing (LRS) technique based on the Oxford Nanopore Technology MinION platform was used for quantifying and kinetic characterization of the poly(A) fraction of bovine alphaherpesvirus type 1 (BoHV-1) lytic transcriptome across a 12-h infection period. Amplification-based LRS techniques frequently generate artefactual transcription reads and are biased towards the production of shorter amplicons. To avoid these undesired effects, we applied direct cDNA sequencing, an amplification-free technique. Here, we show that a single promoter can produce multiple transcription start sites whose distribution patterns differ among the viral genes but are similar in the same gene at different timepoints. Our investigations revealed that the circ gene is expressed with immediate–early (IE) kinetics by utilizing a special mechanism based on the use of the promoter of another IE gene (bicp4) for the transcriptional control. Furthermore, we detected an overlap between the initiation of DNA replication and the transcription from the bicp22 gene, which suggests an interaction between the two molecular machineries. This study developed a generally applicable LRS-based method for the time-course characterization of transcriptomes of any organism.
Collapse
Affiliation(s)
- Dóra Tombácz
- Department of Medical Biology, Albert Szent-Györgyi Medical School, University of Szeged, Somogyi u. 4, 6720 Szeged, Hungary; (D.T.); (B.K.); (G.T.); (Z.C.); (G.G.)
| | - Balázs Kakuk
- Department of Medical Biology, Albert Szent-Györgyi Medical School, University of Szeged, Somogyi u. 4, 6720 Szeged, Hungary; (D.T.); (B.K.); (G.T.); (Z.C.); (G.G.)
| | - Gábor Torma
- Department of Medical Biology, Albert Szent-Györgyi Medical School, University of Szeged, Somogyi u. 4, 6720 Szeged, Hungary; (D.T.); (B.K.); (G.T.); (Z.C.); (G.G.)
| | - Zsolt Csabai
- Department of Medical Biology, Albert Szent-Györgyi Medical School, University of Szeged, Somogyi u. 4, 6720 Szeged, Hungary; (D.T.); (B.K.); (G.T.); (Z.C.); (G.G.)
| | - Gábor Gulyás
- Department of Medical Biology, Albert Szent-Györgyi Medical School, University of Szeged, Somogyi u. 4, 6720 Szeged, Hungary; (D.T.); (B.K.); (G.T.); (Z.C.); (G.G.)
| | - Vivien Tamás
- Institute for Veterinary Medical Research, Centre for Agricultural Research, Hungária krt. 21, 1143 Budapest, Hungary; (V.T.); (Z.Z.)
| | - Zoltán Zádori
- Institute for Veterinary Medical Research, Centre for Agricultural Research, Hungária krt. 21, 1143 Budapest, Hungary; (V.T.); (Z.Z.)
| | - Victoria A. Jefferson
- Department of Biochemistry & Molecular Biology, Entomology & Plant Pathology, Mississippi State University, 408 Dorman P.O. Box 9655, 32 Creelman St., Starkville, MS 39762, USA; (V.A.J.); (F.M.)
| | - Florencia Meyer
- Department of Biochemistry & Molecular Biology, Entomology & Plant Pathology, Mississippi State University, 408 Dorman P.O. Box 9655, 32 Creelman St., Starkville, MS 39762, USA; (V.A.J.); (F.M.)
| | - Zsolt Boldogkői
- Department of Medical Biology, Albert Szent-Györgyi Medical School, University of Szeged, Somogyi u. 4, 6720 Szeged, Hungary; (D.T.); (B.K.); (G.T.); (Z.C.); (G.G.)
- Correspondence:
| |
Collapse
|
11
|
Karlsson M, Sjöstedt E, Oksvold P, Sivertsson Å, Huang J, Álvez MB, Arif M, Li X, Lin L, Yu J, Ma T, Xu F, Han P, Jiang H, Mardinoglu A, Zhang C, von Feilitzen K, Xu X, Wang J, Yang H, Bolund L, Zhong W, Fagerberg L, Lindskog C, Pontén F, Mulder J, Luo Y, Uhlen M. Genome-wide annotation of protein-coding genes in pig. BMC Biol 2022; 20:25. [PMID: 35073880 PMCID: PMC8788080 DOI: 10.1186/s12915-022-01229-y] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/01/2021] [Accepted: 01/07/2022] [Indexed: 12/25/2022] Open
Abstract
BACKGROUND There is a need for functional genome-wide annotation of the protein-coding genes to get a deeper understanding of mammalian biology. Here, a new annotation strategy is introduced based on dimensionality reduction and density-based clustering of whole-body co-expression patterns. This strategy has been used to explore the gene expression landscape in pig, and we present a whole-body map of all protein-coding genes in all major pig tissues and organs. RESULTS An open-access pig expression map ( www.rnaatlas.org ) is presented based on the expression of 350 samples across 98 well-defined pig tissues divided into 44 tissue groups. A new UMAP-based classification scheme is introduced, in which all protein-coding genes are stratified into tissue expression clusters based on body-wide expression profiles. The distribution and tissue specificity of all 22,342 protein-coding pig genes are presented. CONCLUSIONS Here, we present a new genome-wide annotation strategy based on dimensionality reduction and density-based clustering. A genome-wide resource of the transcriptome map across all major tissues and organs in pig is presented, and the data is available as an open-access resource ( www.rnaatlas.org ), including a comparison to the expression of human orthologs.
Collapse
Affiliation(s)
- Max Karlsson
- Department of Protein Science, Science for Life Laboratory, KTH-Royal Institute of Technology, Stockholm, Sweden
| | - Evelina Sjöstedt
- Department of Neuroscience, Karolinska Institutet, Stockholm, Sweden
- Department of Immunology, Genetics and Pathology, Uppsala University, Uppsala, Sweden
| | - Per Oksvold
- Department of Protein Science, Science for Life Laboratory, KTH-Royal Institute of Technology, Stockholm, Sweden
| | - Åsa Sivertsson
- Department of Protein Science, Science for Life Laboratory, KTH-Royal Institute of Technology, Stockholm, Sweden
| | - Jinrong Huang
- BGI-Shenzhen, Shenzhen, China
- Lars Bolund Institute of Regenerative Medicine, Qingdao-Europe Advanced Institute for Life Sciences, BGI-Qingdao, Qingdao, China
- Department of Biomedicine, Aarhus University, Aarhus, Denmark
| | - María Bueno Álvez
- Department of Protein Science, Science for Life Laboratory, KTH-Royal Institute of Technology, Stockholm, Sweden
| | - Muhammad Arif
- Department of Protein Science, Science for Life Laboratory, KTH-Royal Institute of Technology, Stockholm, Sweden
| | - Xiangyu Li
- Department of Protein Science, Science for Life Laboratory, KTH-Royal Institute of Technology, Stockholm, Sweden
| | - Lin Lin
- Department of Biomedicine, Aarhus University, Aarhus, Denmark
- Steno Diabetes Center Aarhus, Aarhus University Hospital, Aarhus, Denmark
| | - Jiaying Yu
- BGI-Shenzhen, Shenzhen, China
- Lars Bolund Institute of Regenerative Medicine, Qingdao-Europe Advanced Institute for Life Sciences, BGI-Qingdao, Qingdao, China
| | - Tao Ma
- MGI, BGI-Shenzhen, Shenzhen, China
| | - Fengping Xu
- BGI-Shenzhen, Shenzhen, China
- Lars Bolund Institute of Regenerative Medicine, Qingdao-Europe Advanced Institute for Life Sciences, BGI-Qingdao, Qingdao, China
| | - Peng Han
- Lars Bolund Institute of Regenerative Medicine, Qingdao-Europe Advanced Institute for Life Sciences, BGI-Qingdao, Qingdao, China
| | | | - Adil Mardinoglu
- Department of Protein Science, Science for Life Laboratory, KTH-Royal Institute of Technology, Stockholm, Sweden
| | - Cheng Zhang
- Department of Protein Science, Science for Life Laboratory, KTH-Royal Institute of Technology, Stockholm, Sweden
| | - Kalle von Feilitzen
- Department of Protein Science, Science for Life Laboratory, KTH-Royal Institute of Technology, Stockholm, Sweden
| | - Xun Xu
- BGI-Shenzhen, Shenzhen, China
| | | | | | - Lars Bolund
- BGI-Shenzhen, Shenzhen, China
- Lars Bolund Institute of Regenerative Medicine, Qingdao-Europe Advanced Institute for Life Sciences, BGI-Qingdao, Qingdao, China
- Department of Biomedicine, Aarhus University, Aarhus, Denmark
| | - Wen Zhong
- Department of Protein Science, Science for Life Laboratory, KTH-Royal Institute of Technology, Stockholm, Sweden
| | - Linn Fagerberg
- Department of Protein Science, Science for Life Laboratory, KTH-Royal Institute of Technology, Stockholm, Sweden
| | - Cecilia Lindskog
- Department of Immunology, Genetics and Pathology, Uppsala University, Uppsala, Sweden
| | - Fredrik Pontén
- Department of Immunology, Genetics and Pathology, Uppsala University, Uppsala, Sweden
| | - Jan Mulder
- Department of Neuroscience, Karolinska Institutet, Stockholm, Sweden
| | - Yonglun Luo
- BGI-Shenzhen, Shenzhen, China
- Lars Bolund Institute of Regenerative Medicine, Qingdao-Europe Advanced Institute for Life Sciences, BGI-Qingdao, Qingdao, China
- Department of Biomedicine, Aarhus University, Aarhus, Denmark
- Steno Diabetes Center Aarhus, Aarhus University Hospital, Aarhus, Denmark
| | - Mathias Uhlen
- Department of Protein Science, Science for Life Laboratory, KTH-Royal Institute of Technology, Stockholm, Sweden.
- Department of Neuroscience, Karolinska Institutet, Stockholm, Sweden.
| |
Collapse
|
12
|
Hao W, Yang Z, Sun Y, Li J, Zhang D, Liu D, Yang X. Characterization of Alternative Splicing Events in Porcine Skeletal Muscles with Different Intramuscular Fat Contents. Biomolecules 2022; 12:biom12020154. [PMID: 35204660 PMCID: PMC8961525 DOI: 10.3390/biom12020154] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/03/2021] [Revised: 01/11/2022] [Accepted: 01/14/2022] [Indexed: 02/06/2023] Open
Abstract
Meat quality is one of the most important economic traits in pig breeding and production. Intramuscular fat (IMF) is a major factor that improves meat quality. To better understand the alternative splicing (AS) events underlying meat quality, long-read isoform sequencing (Iso-seq) was used to identify differential (D)AS events between the longissimus thoracis (LT) and semitendinosus (ST), which differ in IMF content, together with short-read RNA-seq. Through Iso-seq analysis, we identified a total of 56,789 novel transcripts covering protein-coding genes, lncRNA, and fusion transcripts that were not previously annotated in pigs. We also identified 456,965 AS events, among which 3930 were DAS events, corresponding to 2364 unique genes. Through integrative analysis of Iso-seq and RNA-seq, we identified 1174 differentially expressed genes (DEGs), among which 122 were DAS genes, i.e., DE-DAS genes. There are 12 overlapped pathways between the top 20 DEGs and DE-DAS genes, as revealed by KEGG (Kyoto Encyclopedia of Genes and Genomes) analysis, indicating that DE-DAS genes play important roles in the differential phenotype of LT and ST. Further analysis showed that upregulated DE-DAS genes are more important than downregulated ones in IMF deposition. Fatty acid degradation and the PPAR (peroxisome proliferator-activated receptor) signaling pathway were found to be the most important pathways regulating the differential fat deposition of the two muscles. The results update the existing porcine genome annotations and provide data for the in-depth exploration of the mechanisms underlying meat quality and IMF deposition.
Collapse
Affiliation(s)
- Wanjun Hao
- College of Animal Science and Technology, Northeast Agricultural University, Harbin 150030, China; (W.H.); (Z.Y.); (Y.S.); (J.L.)
| | - Zewei Yang
- College of Animal Science and Technology, Northeast Agricultural University, Harbin 150030, China; (W.H.); (Z.Y.); (Y.S.); (J.L.)
| | - Yuanlu Sun
- College of Animal Science and Technology, Northeast Agricultural University, Harbin 150030, China; (W.H.); (Z.Y.); (Y.S.); (J.L.)
| | - Jiaxin Li
- College of Animal Science and Technology, Northeast Agricultural University, Harbin 150030, China; (W.H.); (Z.Y.); (Y.S.); (J.L.)
| | - Dongjie Zhang
- Institute of Animal Husbandry, Heilongjiang Academy of Agricultural Sciences, Harbin 150086, China;
| | - Di Liu
- Institute of Animal Husbandry, Heilongjiang Academy of Agricultural Sciences, Harbin 150086, China;
- Correspondence: (D.L.); (X.Y.); Tel.: +86-451-8667-7458 (D.L.); +86-451-5519-1738 (X.Y.)
| | - Xiuqin Yang
- College of Animal Science and Technology, Northeast Agricultural University, Harbin 150030, China; (W.H.); (Z.Y.); (Y.S.); (J.L.)
- Correspondence: (D.L.); (X.Y.); Tel.: +86-451-8667-7458 (D.L.); +86-451-5519-1738 (X.Y.)
| |
Collapse
|
13
|
Treitli SC, Peña-Diaz P, Hałakuc P, Karnkowska A, Hampl V. High quality genome assembly of the amitochondriate eukaryote Monocercomonoides exilis. Microb Genom 2021; 7. [PMID: 34951395 PMCID: PMC8767320 DOI: 10.1099/mgen.0.000745] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Abstract
Monocercomonoides exilis is considered the first known eukaryote to completely lack mitochondria. This conclusion is based primarily on a genomic and transcriptomic study which failed to identify any mitochondrial hallmark proteins. However, the available genome assembly has limited contiguity and around 1.5 % of the genome sequence is represented by unknown bases. To improve the contiguity, we re-sequenced the genome and transcriptome of M. exilis using Oxford Nanopore Technology (ONT). The resulting draft genome is assembled in 101 contigs with an N50 value of 1.38 Mbp, almost 20 times higher than the previously published assembly. Using a newly generated ONT transcriptome, we further improve the gene prediction and add high quality untranslated region (UTR) annotations, in which we identify two putative polyadenylation signals present in the 3′UTR regions and characterise the Kozak sequence in the 5′UTR regions. All these improvements are reflected by higher BUSCO genome completeness values. Regardless of an overall more complete genome assembly without missing bases and a better gene prediction, we still failed to identify any mitochondrial hallmark genes, thus further supporting the hypothesis on the absence of mitochondrion.
Collapse
Affiliation(s)
- Sebastian Cristian Treitli
- Department of Parasitology, Faculty of Science, Charles University, BIOCEV, Průmyslová 595, 252 42 Vestec, Czech Republic
| | - Priscila Peña-Diaz
- Department of Parasitology, Faculty of Science, Charles University, BIOCEV, Průmyslová 595, 252 42 Vestec, Czech Republic
| | - Paweł Hałakuc
- Institute of Evolutionary Biology, Faculty of Biology, Biological and Chemical Research Centre, University of Warsaw, Warsaw, Poland
| | - Anna Karnkowska
- Institute of Evolutionary Biology, Faculty of Biology, Biological and Chemical Research Centre, University of Warsaw, Warsaw, Poland
| | - Vladimír Hampl
- Department of Parasitology, Faculty of Science, Charles University, BIOCEV, Průmyslová 595, 252 42 Vestec, Czech Republic
| |
Collapse
|
14
|
Li T, Zhang X, Guo L, Qi T, Tang H, Wang H, Qiao X, Zhang M, Zhang B, Feng J, Zuo Z, Zhang Y, Xing C, Wu J. Single-molecule real-time transcript sequencing of developing cotton anthers facilitates genome annotation and fertility restoration candidate gene discovery. Genomics 2021; 113:4245-4253. [PMID: 34793949 DOI: 10.1016/j.ygeno.2021.11.014] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/03/2020] [Revised: 07/04/2021] [Accepted: 11/10/2021] [Indexed: 01/23/2023]
Abstract
Heterosis refers to the superior phenotypes observed in hybrids. Cytoplasmic male sterility (CMS) system plays an important role in cotton heterosis utilization. However, the global gene expression patterns of CMS-D2 and its interaction with the restorer gene Rf1 remain unclear. Here, the full-length transcript sequencing was performed in anthers of the CMS-D2 restorer line using PacBio single-molecule real-time sequencing technology. Combining PacBio SMRT long-read isoforms and Illumina RNA-seq data, 107,066 isoforms from 44,338 loci were obtained, including 10,086 novel isoforms of novel genes and 66,419 new isoforms of known genes. Totally 56,572 alternative splicing (AS) events, 1146 lncRNAs, 61 fusion transcripts and 10,466 genes exhibited alternative polyadenylation (APA), and 60,995 novel isoforms with predicted open reading frames (ORFs) were further identified. Furthermore, the specifically expressed genes in restorer line were selected and confirmed by qRT-PCR. These findings provide a basis for upland cotton genome annotation and transcriptome research, and will help to reveal the molecular mechanism of interaction between Rf1 and CMS-D2 cytoplasm.
Collapse
Affiliation(s)
- Ting Li
- Zhengzhou Research Base, State Key Laboratory of Cotton Biology, Zhengzhou University, Zhengzhou 450000, Henan, China
| | - Xuexian Zhang
- State Key Laboratory of Cotton Biology, Institute of Cotton Research of Chinese Academy of Agricultural Sciences, Key Laboratory for Cotton Genetic Improvement, Ministry of Agriculture, 38 Huanghe Dadao, Anyang 455000, Henan, China.
| | - Liping Guo
- State Key Laboratory of Cotton Biology, Institute of Cotton Research of Chinese Academy of Agricultural Sciences, Key Laboratory for Cotton Genetic Improvement, Ministry of Agriculture, 38 Huanghe Dadao, Anyang 455000, Henan, China.
| | - Tingxiang Qi
- State Key Laboratory of Cotton Biology, Institute of Cotton Research of Chinese Academy of Agricultural Sciences, Key Laboratory for Cotton Genetic Improvement, Ministry of Agriculture, 38 Huanghe Dadao, Anyang 455000, Henan, China.
| | - Huini Tang
- State Key Laboratory of Cotton Biology, Institute of Cotton Research of Chinese Academy of Agricultural Sciences, Key Laboratory for Cotton Genetic Improvement, Ministry of Agriculture, 38 Huanghe Dadao, Anyang 455000, Henan, China.
| | - Hailin Wang
- State Key Laboratory of Cotton Biology, Institute of Cotton Research of Chinese Academy of Agricultural Sciences, Key Laboratory for Cotton Genetic Improvement, Ministry of Agriculture, 38 Huanghe Dadao, Anyang 455000, Henan, China
| | - Xiuqin Qiao
- State Key Laboratory of Cotton Biology, Institute of Cotton Research of Chinese Academy of Agricultural Sciences, Key Laboratory for Cotton Genetic Improvement, Ministry of Agriculture, 38 Huanghe Dadao, Anyang 455000, Henan, China.
| | - Meng Zhang
- State Key Laboratory of Cotton Biology, Institute of Cotton Research of Chinese Academy of Agricultural Sciences, Key Laboratory for Cotton Genetic Improvement, Ministry of Agriculture, 38 Huanghe Dadao, Anyang 455000, Henan, China
| | - Bingbing Zhang
- State Key Laboratory of Cotton Biology, Institute of Cotton Research of Chinese Academy of Agricultural Sciences, Key Laboratory for Cotton Genetic Improvement, Ministry of Agriculture, 38 Huanghe Dadao, Anyang 455000, Henan, China
| | - Juanjuan Feng
- State Key Laboratory of Cotton Biology, Institute of Cotton Research of Chinese Academy of Agricultural Sciences, Key Laboratory for Cotton Genetic Improvement, Ministry of Agriculture, 38 Huanghe Dadao, Anyang 455000, Henan, China
| | - Zhidan Zuo
- State Key Laboratory of Cotton Biology, Institute of Cotton Research of Chinese Academy of Agricultural Sciences, Key Laboratory for Cotton Genetic Improvement, Ministry of Agriculture, 38 Huanghe Dadao, Anyang 455000, Henan, China
| | - Yongjie Zhang
- Zhengzhou Research Base, State Key Laboratory of Cotton Biology, Zhengzhou University, Zhengzhou 450000, Henan, China
| | - Chaozhu Xing
- State Key Laboratory of Cotton Biology, Institute of Cotton Research of Chinese Academy of Agricultural Sciences, Key Laboratory for Cotton Genetic Improvement, Ministry of Agriculture, 38 Huanghe Dadao, Anyang 455000, Henan, China.
| | - Jianyong Wu
- Zhengzhou Research Base, State Key Laboratory of Cotton Biology, Zhengzhou University, Zhengzhou 450000, Henan, China; State Key Laboratory of Cotton Biology, Institute of Cotton Research of Chinese Academy of Agricultural Sciences, Key Laboratory for Cotton Genetic Improvement, Ministry of Agriculture, 38 Huanghe Dadao, Anyang 455000, Henan, China.
| |
Collapse
|
15
|
Yuan Z, Ge L, Sun J, Zhang W, Wang S, Cao X, Sun W. Integrative analysis of Iso-Seq and RNA-seq data reveals transcriptome complexity and differentially expressed transcripts in sheep tail fat. PeerJ 2021; 9:e12454. [PMID: 34760406 PMCID: PMC8571958 DOI: 10.7717/peerj.12454] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2021] [Accepted: 10/18/2021] [Indexed: 01/22/2023] Open
Abstract
Background Nowadays, both customers and producers prefer thin-tailed fat sheep. To effectively breed for this phenotype, it is important to identify candidate genes and uncover the genetic mechanism related to tail fat deposition in sheep. Accumulating evidence suggesting that post-transcriptional modification events of precursor-messenger RNA (pre-mRNA), including alternative splicing (AS) and alternative polyadenylation (APA), may regulate tail fat deposition in sheep. Differentially expressed transcripts (DETs) analysis is a way to identify candidate genes related to tail fat deposition. However, due to the technological limitation, post-transcriptional modification events in the tail fat of sheep and DETs between thin-tailed and fat-tailed sheep remains unclear. Methods In the present study, we applied pooled PacBio isoform sequencing (Iso-Seq) to generate transcriptomic data of tail fat tissue from six sheep (three thin-tailed sheep and three fat-tailed sheep). By comparing with reference genome, potential gene loci and novel transcripts were identified. Post-transcriptional modification events, including AS and APA, and lncRNA in sheep tail fat were uncovered using pooled Iso-Seq data. Combining Iso-Seq data with six RNA-sequencing (RNA-Seq) data, DETs between thin- and fat-tailed sheep were identified. Protein protein interaction (PPI) network, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were implemented to investigate the potential functions of DETs. Results In the present study, we revealed the transcriptomic complexity of the tail fat of sheep, result in 9,001 potential novel gene loci, 17,834 AS events, 5,791 APA events, and 3,764 lncRNAs. Combining Iso-Seq data with RNA-Seq data, we identified hundreds of DETs between thin- and fat-tailed sheep. Among them, 21 differentially expressed lncRNAs, such as ENSOART00020036299, ENSOART00020033641, ENSOART00020024562, ENSOART00020003848 and 9.53.1 may regulate tail fat deposition. Many novel transcripts were identified as DETs, including 15.527.13 (DGAT2), 13.624.23 (ACSS2), 11.689.28 (ACLY), 11.689.18 (ACLY), 11.689.14 (ACLY), 11.660.12 (ACLY), 22.289.6 (SCD), 22.289.3 (SCD) and 22.289.14 (SCD). Most of the identified DETs have been enriched in GO and KEGG pathways related to extracellular matrix (ECM). Our result revealed the transcriptome complexity and identified many candidate transcripts in tail fat, which could enhance the understanding of molecular mechanisms behind tail fat deposition.
Collapse
Affiliation(s)
- Zehu Yuan
- Joint International Research Laboratory of Agriculture and Agri-Product Safety of Ministry of Education, Yangzhou University, Yangzhou, China
| | - Ling Ge
- College of Animal Science and Technology, Yangzhou University, Yangzhou, China
| | - Jingyi Sun
- College of Veterinary Medicine, Yangzhou University, Yangzhou, China
| | - Weibo Zhang
- College of Animal Science and Technology, Yangzhou University, Yangzhou, China
| | - Shanhe Wang
- College of Animal Science and Technology, Yangzhou University, Yangzhou, China
| | - Xiukai Cao
- Joint International Research Laboratory of Agriculture and Agri-Product Safety of Ministry of Education, Yangzhou University, Yangzhou, China
| | - Wei Sun
- Joint International Research Laboratory of Agriculture and Agri-Product Safety of Ministry of Education, Yangzhou University, Yangzhou, China.,College of Animal Science and Technology, Yangzhou University, Yangzhou, China
| |
Collapse
|
16
|
Li Y, Wang J, Elzo MA, Fan H, Du K, Xia S, Shao J, Lai T, Hu S, Jia X, Lai S. Molecular Profiling of DNA Methylation and Alternative Splicing of Genes in Skeletal Muscle of Obese Rabbits. Curr Issues Mol Biol 2021; 43:1558-1575. [PMID: 34698087 PMCID: PMC8929151 DOI: 10.3390/cimb43030110] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/08/2021] [Revised: 10/05/2021] [Accepted: 10/07/2021] [Indexed: 12/11/2022] Open
Abstract
DNA methylation and the alternative splicing of precursor messenger RNAs (pre-mRNAs) are two important genetic modification mechanisms. However, both are currently uncharacterized in the muscle metabolism of rabbits. Thus, we constructed the Tianfu black rabbit obesity model (obese rabbits fed with a 10% high-fat diet and control rabbits from 35 days to 70 days) and collected the skeletal muscle samples from the two groups for Genome methylation sequencing and RNA sequencing. DNA methylation data showed that the promoter regions of 599 genes and gene body region of 2522 genes had significantly differential methylation rates between the two groups, of which 288 genes had differential methylation rates in promoter and gene body regions. Analysis of alternative splicing showed 555 genes involved in exon skipping (ES) patterns, and 15 genes existed in differential methylation regions. Network analysis showed that 20 hub genes were associated with ubiquitinated protein degradation, muscle development pathways, and skeletal muscle energy metabolism. Our findings suggest that the two types of genetic modification have potential regulatory effects on skeletal muscle development and provide a basis for further mechanistic studies in the rabbit.
Collapse
Affiliation(s)
- Yanhong Li
- College of Animal Science and Technology, Sichuan Agricultural University, Chengdu 611130, China; (Y.L.); (J.W.); (H.F.); (K.D.); (S.X.); (J.S.); (T.L.); (S.H.); (X.J.)
| | - Jie Wang
- College of Animal Science and Technology, Sichuan Agricultural University, Chengdu 611130, China; (Y.L.); (J.W.); (H.F.); (K.D.); (S.X.); (J.S.); (T.L.); (S.H.); (X.J.)
| | - Mauricio A. Elzo
- Department of Animal Sciences, University of Florida, Gainesville, FL 32611, USA;
| | - Huimei Fan
- College of Animal Science and Technology, Sichuan Agricultural University, Chengdu 611130, China; (Y.L.); (J.W.); (H.F.); (K.D.); (S.X.); (J.S.); (T.L.); (S.H.); (X.J.)
| | - Kun Du
- College of Animal Science and Technology, Sichuan Agricultural University, Chengdu 611130, China; (Y.L.); (J.W.); (H.F.); (K.D.); (S.X.); (J.S.); (T.L.); (S.H.); (X.J.)
| | - Siqi Xia
- College of Animal Science and Technology, Sichuan Agricultural University, Chengdu 611130, China; (Y.L.); (J.W.); (H.F.); (K.D.); (S.X.); (J.S.); (T.L.); (S.H.); (X.J.)
| | - Jiahao Shao
- College of Animal Science and Technology, Sichuan Agricultural University, Chengdu 611130, China; (Y.L.); (J.W.); (H.F.); (K.D.); (S.X.); (J.S.); (T.L.); (S.H.); (X.J.)
| | - Tianfu Lai
- College of Animal Science and Technology, Sichuan Agricultural University, Chengdu 611130, China; (Y.L.); (J.W.); (H.F.); (K.D.); (S.X.); (J.S.); (T.L.); (S.H.); (X.J.)
| | - Shenqiang Hu
- College of Animal Science and Technology, Sichuan Agricultural University, Chengdu 611130, China; (Y.L.); (J.W.); (H.F.); (K.D.); (S.X.); (J.S.); (T.L.); (S.H.); (X.J.)
| | - Xianbo Jia
- College of Animal Science and Technology, Sichuan Agricultural University, Chengdu 611130, China; (Y.L.); (J.W.); (H.F.); (K.D.); (S.X.); (J.S.); (T.L.); (S.H.); (X.J.)
| | - Songjia Lai
- College of Animal Science and Technology, Sichuan Agricultural University, Chengdu 611130, China; (Y.L.); (J.W.); (H.F.); (K.D.); (S.X.); (J.S.); (T.L.); (S.H.); (X.J.)
- Correspondence:
| |
Collapse
|
17
|
Liu J, Bao Y, Zhong Y, Wang Q, Liu H. Genome-wide association study and transcriptome of olecranon-type traits in peach (Prunus persica L.) germplasm. BMC Genomics 2021; 22:702. [PMID: 34583632 PMCID: PMC8480057 DOI: 10.1186/s12864-021-08017-y] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2021] [Accepted: 09/16/2021] [Indexed: 01/24/2023] Open
Abstract
Background The top of the olecranon honey peach (Prunus persica L.) fruit appears similar to an eagle’s beak. In this study, a single olecranon honey peach with a round-type fruit was observed in our fruit orchard. To explore the genetic mechanism of olecranon formation, we performed full-length transcriptome sequencing analysis of olecranon and round peaches as well as a genome-wide association study of the association of olecranon-type trait loci. Results The gene locus was 26,924,482 base pairs in NC_034014.1. Transcriptome sequencing showed that the clean sequencing data of each sample reached 7.10GB, with 14,360 genes and 23,167 transcripts expressed in both the olecranon honey peach and round peach. Among the 11 differentially expressed genes selected as candidate genes, six were highly expressed in olecranon peach and named as LOC18775282, LOC18772209, LOC18773929, LOC18772013, LOC18773401, and ONT.13798.5. Five genes were highly expressed in round peach and named as LOC18773079, LOC18773525, LOC18773067, LOC18775244, and LOC18772236. Notably, ONT.13798.5 was not previously identified. The genes were within 1 Mb up- or down-stream of the main genome-wide association study locus for olecranon-type traits. Conclusions This study revealed loci associated with olecranon and provides useful information for analysis and breeding of olecranon honey peach. Supplementary Information The online version contains supplementary material available at 10.1186/s12864-021-08017-y.
Collapse
Affiliation(s)
- Jianliang Liu
- College of Light Industry and Food, Zhongkai University of Agriculture and Engineering, 510225, Guangzhou, Guangdong, China.,Guangzhou Key Laboratory for Research and Development of Crop Germplasm Resources, Zhongkai University of Agriculture and Engineering, 510225, Guangzhou, China.,Modern Agriculture Research Center, Zhongkai University of Agriculture and Engineering, 510225, Guangzhou, Guangdong, China
| | - Yao Bao
- College of Light Industry and Food, Zhongkai University of Agriculture and Engineering, 510225, Guangzhou, Guangdong, China.,Guangdong Provincial Key Laboratory of Lingnan Specialty Food Science and Technology, Guangzhou, Guangdong, China
| | - Yuming Zhong
- College of Environmental Science and Engineering, Zhongkai University of Agriculture and Engineering, 510225, Guangzhou, Guangdong, China
| | - Qin Wang
- College of Light Industry and Food, Zhongkai University of Agriculture and Engineering, 510225, Guangzhou, Guangdong, China.,Guangdong Provincial Key Laboratory of Lingnan Specialty Food Science and Technology, Guangzhou, Guangdong, China
| | - Huifan Liu
- College of Light Industry and Food, Zhongkai University of Agriculture and Engineering, 510225, Guangzhou, Guangdong, China. .,Guangdong Provincial Key Laboratory of Lingnan Specialty Food Science and Technology, Guangzhou, Guangdong, China.
| |
Collapse
|
18
|
Chang T, An B, Liang M, Duan X, Du L, Cai W, Zhu B, Gao X, Chen Y, Xu L, Zhang L, Gao H, Li J. PacBio Single-Molecule Long-Read Sequencing Provides New Light on the Complexity of Full-Length Transcripts in Cattle. Front Genet 2021; 12:664974. [PMID: 34527015 PMCID: PMC8437344 DOI: 10.3389/fgene.2021.664974] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/06/2021] [Accepted: 08/06/2021] [Indexed: 12/02/2022] Open
Abstract
Cattle (Bos taurus) is one of the most widely distributed livestock species in the world, and provides us with high-quality milk and meat which have a huge impact on the quality of human life. Therefore, accurate and complete transcriptome and genome annotation are of great value to the research of cattle breeding. In this study, we used error-corrected PacBio single-molecule real-time (SMRT) data to perform whole-transcriptome profiling in cattle. Then, 22.5 Gb of subreads was generated, including 381,423 circular consensus sequences (CCSs), among which 276,295 full-length non-chimeric (FLNC) sequences were identified. After correction by Illumina short reads, we obtained 22,353 error-corrected isoforms. A total of 305 alternative splicing (AS) events and 3,795 alternative polyadenylation (APA) sites were detected by transcriptome structural analysis. Furthermore, we identified 457 novel genes, 120 putative transcription factors (TFs), and 569 novel long non-coding RNAs (lncRNAs). Taken together, this research improves our understanding and provides new insights into the complexity of full-length transcripts in cattle.
Collapse
Affiliation(s)
- Tianpeng Chang
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, China
| | - Bingxing An
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, China
| | - Mang Liang
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, China
| | - Xinghai Duan
- College of Animal Science and Technology, Southwest University, Chongqing, China
| | - Lili Du
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, China
| | - Wentao Cai
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, China
| | - Bo Zhu
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, China
| | - Xue Gao
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, China
| | - Yan Chen
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, China
| | - Lingyang Xu
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, China
| | - Lupei Zhang
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, China
| | - Huijiang Gao
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, China
| | - Junya Li
- Laboratory of Molecular Biology and Bovine Breeding, Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, China
| |
Collapse
|
19
|
Ran X, Hu F, Mao N, Ruan Y, Yi F, Niu X, Huang S, Li S, You L, Zhang F, Tang L, Wang J, Liu J. Differences in gene expression and variable splicing events of ovaries between large and small litter size in Chinese Xiang pigs. Porcine Health Manag 2021; 7:52. [PMID: 34470660 PMCID: PMC8411529 DOI: 10.1186/s40813-021-00226-x] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/29/2021] [Accepted: 07/20/2021] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND Although lots of quantitative trait loci (QTLs) and genes present roles in litter size of some breeds, the information might not make it clear for the huge diversity of reproductive capability in pig breeds. To elucidate the inherent mechanisms of heterogeneity of reproductive capability in litter size of Xiang pig, we performed transcriptome analysis for the expression profile in ovaries using RNA-seq method. RESULTS We identified 1,419 up-regulated and 1,376 down-regulated genes in Xiang pigs with large litter size. Among them, 1,010 differentially expressed genes (DEGs) were differently spliced between two groups with large or small litter sizes. Based on GO and KEGG analysis, numerous members of genes were gathered in ovarian steroidogenesis, steroid biosynthesis, oocyte maturation and reproduction processes. CONCLUSIONS Combined with gene biological function, twelve genes were found out that might be related with the reproductive capability of Xiang pig, of which, eleven genes were recognized as hub genes. These genes may play a role in promoting litter size by elevating steroid and peptide hormones supply through the ovary and facilitating the processes of ovulation and in vivo fertilization.
Collapse
Affiliation(s)
- Xueqin Ran
- College of Animal Science, Institute of Agro-Bioengineering and Key Laboratory of Plant Resource Conservative and Germplam Innovation in Mountainous Region (Ministry of Education), Guizhou University, 550025, Guiyang, China
| | - Fengbin Hu
- College of Animal Science, Institute of Agro-Bioengineering and Key Laboratory of Plant Resource Conservative and Germplam Innovation in Mountainous Region (Ministry of Education), Guizhou University, 550025, Guiyang, China
| | - Ning Mao
- College of Animal Science, Institute of Agro-Bioengineering and Key Laboratory of Plant Resource Conservative and Germplam Innovation in Mountainous Region (Ministry of Education), Guizhou University, 550025, Guiyang, China
| | - Yiqi Ruan
- College of Animal Science, Institute of Agro-Bioengineering and Key Laboratory of Plant Resource Conservative and Germplam Innovation in Mountainous Region (Ministry of Education), Guizhou University, 550025, Guiyang, China
| | - Fanli Yi
- College of Animal Science, Institute of Agro-Bioengineering and Key Laboratory of Plant Resource Conservative and Germplam Innovation in Mountainous Region (Ministry of Education), Guizhou University, 550025, Guiyang, China
| | - Xi Niu
- College of Animal Science, Institute of Agro-Bioengineering and Key Laboratory of Plant Resource Conservative and Germplam Innovation in Mountainous Region (Ministry of Education), Guizhou University, 550025, Guiyang, China
| | - Shihui Huang
- College of Animal Science, Institute of Agro-Bioengineering and Key Laboratory of Plant Resource Conservative and Germplam Innovation in Mountainous Region (Ministry of Education), Guizhou University, 550025, Guiyang, China
| | - Sheng Li
- College of Animal Science, Institute of Agro-Bioengineering and Key Laboratory of Plant Resource Conservative and Germplam Innovation in Mountainous Region (Ministry of Education), Guizhou University, 550025, Guiyang, China
| | - Longjiang You
- College of Animal Science, Institute of Agro-Bioengineering and Key Laboratory of Plant Resource Conservative and Germplam Innovation in Mountainous Region (Ministry of Education), Guizhou University, 550025, Guiyang, China
| | - Fuping Zhang
- College of Animal Science, Institute of Agro-Bioengineering and Key Laboratory of Plant Resource Conservative and Germplam Innovation in Mountainous Region (Ministry of Education), Guizhou University, 550025, Guiyang, China
| | - Liangting Tang
- College of Animal Science, Institute of Agro-Bioengineering and Key Laboratory of Plant Resource Conservative and Germplam Innovation in Mountainous Region (Ministry of Education), Guizhou University, 550025, Guiyang, China
| | - Jiafu Wang
- College of Animal Science, Institute of Agro-Bioengineering and Key Laboratory of Plant Resource Conservative and Germplam Innovation in Mountainous Region (Ministry of Education), Guizhou University, 550025, Guiyang, China.
| | - Jianfeng Liu
- National Engineering Laboratory for Animal Breeding, College of Animal Science and Technology, China Agricultural University, 100193, Beijing, China
| |
Collapse
|
20
|
Naftaly AS, Pau S, White MA. Long-read RNA sequencing reveals widespread sex-specific alternative splicing in threespine stickleback fish. Genome Res 2021; 31:1486-1497. [PMID: 34131005 PMCID: PMC8327910 DOI: 10.1101/gr.274282.120] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/11/2020] [Accepted: 06/15/2021] [Indexed: 01/07/2023]
Abstract
Alternate isoforms are important contributors to phenotypic diversity across eukaryotes. Although short-read RNA-sequencing has increased our understanding of isoform diversity, it is challenging to accurately detect full-length transcripts, preventing the identification of many alternate isoforms. Long-read sequencing technologies have made it possible to sequence full-length alternative transcripts, accurately characterizing alternative splicing events, alternate transcription start and end sites, and differences in UTR regions. Here, we use Pacific Biosciences (PacBio) long-read RNA-sequencing (Iso-Seq) to examine the transcriptomes of five organs in threespine stickleback fish (Gasterosteus aculeatus), a widely used genetic model species. The threespine stickleback fish has a refined genome assembly in which gene annotations are based on short-read RNA sequencing and predictions from coding sequence of other species. This suggests some of the existing annotations may be inaccurate or alternative transcripts may not be fully characterized. Using Iso-Seq we detected thousands of novel isoforms, indicating many isoforms are absent in the current Ensembl gene annotations. In addition, we refined many of the existing annotations within the genome. We noted many improperly positioned transcription start sites that were refined with long-read sequencing. The Iso-Seq-predicted transcription start sites were more accurate and verified through ATAC-seq. We also detected many alternative splicing events between sexes and across organs. We found a substantial number of genes in both somatic and gonadal samples that had sex-specific isoforms. Our study highlights the power of long-read sequencing to study the complexity of transcriptomes, greatly improving genomic resources for the threespine stickleback fish.
Collapse
Affiliation(s)
- Alice S Naftaly
- Department of Genetics, University of Georgia, Athens, Georgia 30602, USA
| | - Shana Pau
- Department of Genetics, University of Georgia, Athens, Georgia 30602, USA
- Department of Biology, University of Texas Arlington, Arlington, Texas 76019, USA
| | - Michael A White
- Department of Genetics, University of Georgia, Athens, Georgia 30602, USA
| |
Collapse
|
21
|
Miao BB, Niu SF, Wu RX, Liang ZB, Tang BG, Zhai Y, Xu XQ. Gene Expression Profile and Co-Expression Network of Pearl Gentian Grouper under Cold Stress by Integrating Illumina and PacBio Sequences. Animals (Basel) 2021; 11:ani11061745. [PMID: 34208015 PMCID: PMC8230743 DOI: 10.3390/ani11061745] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/06/2021] [Revised: 06/07/2021] [Accepted: 06/09/2021] [Indexed: 12/14/2022] Open
Abstract
Simple Summary In this study, we investigated the liver transcriptomic responses of pearl gentian grouper towards cold stress. Some cold-related key genes and biological pathways were screened, of which energy-related metabolic pathways and genes had higher expression levels under cold stress. This suggested that energy homeostasis plays a crucial role in the physiological adjustments of pearl gentian grouper when exposed to the cold stress environment. Our results will expedite the understanding of different fishes adaptive mechanisms to profound environmental temperature changes and provide insights into the molecular breeding of cold-tolerant pearl gentian grouper varieties. Abstract Pearl gentian grouper (Epinephelus fuscoguttatus ♀ × Epinephelus lanceolatus ♂) is a fish of high commercial value in the aquaculture industry in Asia. However, this hybrid fish is not cold-tolerant, and its molecular regulation mechanism underlying cold stress remains largely elusive. This study thus investigated the liver transcriptomic responses of pearl gentian grouper by comparing the gene expression of cold stress groups (20, 15, 12, and 12 °C for 6 h) with that of control group (25 °C) using PacBio SMRT-Seq and Illumina RNA-Seq technologies. In SMRT-Seq analysis, a total of 11,033 full-length transcripts were generated and used as reference sequences for further RNA-Seq analysis. In RNA-Seq analysis, 3271 differentially expressed genes (DEGs), two low-temperature specific modules (tan and blue modules), and two significantly expressed gene sets (profiles 0 and 19) were screened by differential expression analysis, weighted gene co-expression networks analysis (WGCNA), and short time-series expression miner (STEM), respectively. The intersection of the above analyses further revealed some key genes, such as PCK, ALDOB, FBP, G6pC, CPT1A, PPARα, SOCS3, PPP1CC, CYP2J, HMGCR, CDKN1B, and GADD45Bc. These genes were significantly enriched in carbohydrate metabolism, lipid metabolism, signal transduction, and endocrine system pathways. All these pathways were linked to biological functions relevant to cold adaptation, such as energy metabolism, stress-induced cell membrane changes, and transduction of stress signals. Taken together, our study explores an overall and complex regulation network of the functional genes in the liver of pearl gentian grouper, which could benefit the species in preventing damage caused by cold stress.
Collapse
Affiliation(s)
- Ben-Ben Miao
- College of Fisheries, Guangdong Ocean University, Zhanjiang 524088, China; (B.-B.M.); (S.-F.N.); (Z.-B.L.); (B.-G.T.); (Y.Z.); (X.-Q.X.)
| | - Su-Fang Niu
- College of Fisheries, Guangdong Ocean University, Zhanjiang 524088, China; (B.-B.M.); (S.-F.N.); (Z.-B.L.); (B.-G.T.); (Y.Z.); (X.-Q.X.)
- Southern Marine Science and Engineering Guangdong Laboratory, Zhanjiang 524025, China
| | - Ren-Xie Wu
- College of Fisheries, Guangdong Ocean University, Zhanjiang 524088, China; (B.-B.M.); (S.-F.N.); (Z.-B.L.); (B.-G.T.); (Y.Z.); (X.-Q.X.)
- Southern Marine Science and Engineering Guangdong Laboratory, Zhanjiang 524025, China
- Correspondence:
| | - Zhen-Bang Liang
- College of Fisheries, Guangdong Ocean University, Zhanjiang 524088, China; (B.-B.M.); (S.-F.N.); (Z.-B.L.); (B.-G.T.); (Y.Z.); (X.-Q.X.)
| | - Bao-Gui Tang
- College of Fisheries, Guangdong Ocean University, Zhanjiang 524088, China; (B.-B.M.); (S.-F.N.); (Z.-B.L.); (B.-G.T.); (Y.Z.); (X.-Q.X.)
- Southern Marine Science and Engineering Guangdong Laboratory, Zhanjiang 524025, China
| | - Yun Zhai
- College of Fisheries, Guangdong Ocean University, Zhanjiang 524088, China; (B.-B.M.); (S.-F.N.); (Z.-B.L.); (B.-G.T.); (Y.Z.); (X.-Q.X.)
| | - Xue-Qi Xu
- College of Fisheries, Guangdong Ocean University, Zhanjiang 524088, China; (B.-B.M.); (S.-F.N.); (Z.-B.L.); (B.-G.T.); (Y.Z.); (X.-Q.X.)
| |
Collapse
|
22
|
Halstead MM, Islas-Trejo A, Goszczynski DE, Medrano JF, Zhou H, Ross PJ. Large-Scale Multiplexing Permits Full-Length Transcriptome Annotation of 32 Bovine Tissues From a Single Nanopore Flow Cell. Front Genet 2021; 12:664260. [PMID: 34093657 PMCID: PMC8173071 DOI: 10.3389/fgene.2021.664260] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/04/2021] [Accepted: 04/06/2021] [Indexed: 12/18/2022] Open
Abstract
A comprehensive annotation of transcript isoforms in domesticated species is lacking. Especially considering that transcriptome complexity and splicing patterns are not well-conserved between species, this presents a substantial obstacle to genomic selection programs that seek to improve production, disease resistance, and reproduction. Recent advances in long-read sequencing technology have made it possible to directly extrapolate the structure of full-length transcripts without the need for transcript reconstruction. In this study, we demonstrate the power of long-read sequencing for transcriptome annotation by coupling Oxford Nanopore Technology (ONT) with large-scale multiplexing of 93 samples, comprising 32 tissues collected from adult male and female Hereford cattle. More than 30 million uniquely mapping full-length reads were obtained from a single ONT flow cell, and used to identify and characterize the expression dynamics of 99,044 transcript isoforms at 31,824 loci. Of these predicted transcripts, 21% exactly matched a reference transcript, and 61% were novel isoforms of reference genes, substantially increasing the ratio of transcript variants per gene, and suggesting that the complexity of the bovine transcriptome is comparable to that in humans. Over 7,000 transcript isoforms were extremely tissue-specific, and 61% of these were attributed to testis, which exhibited the most complex transcriptome of all interrogated tissues. Despite profiling over 30 tissues, transcription was only detected at about 60% of reference loci. Consequently, additional studies will be necessary to continue characterizing the bovine transcriptome in additional cell types, developmental stages, and physiological conditions. However, by here demonstrating the power of ONT sequencing coupled with large-scale multiplexing, the task of exhaustively annotating the bovine transcriptome - or any mammalian transcriptome - appears significantly more feasible.
Collapse
Affiliation(s)
| | | | | | | | | | - Pablo J. Ross
- Department of Animal Science, University of California, Davis, Davis, CA, United States
| |
Collapse
|
23
|
Tombácz D, Moldován N, Torma G, Nagy T, Hornyák Á, Csabai Z, Gulyás G, Boldogkői M, Jefferson VA, Zádori Z, Meyer F, Boldogkői Z. Dynamic Transcriptome Sequencing of Bovine Alphaherpesvirus Type 1 and Host Cells Carried Out by a Multi-Technique Approach. Front Genet 2021; 12:619056. [PMID: 33897757 PMCID: PMC8059770 DOI: 10.3389/fgene.2021.619056] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2020] [Accepted: 03/16/2021] [Indexed: 12/12/2022] Open
Affiliation(s)
- Dóra Tombácz
- Department of Medical Biology, Faculty of Medicine, University of Szeged, Szeged, Hungary
| | - Norbert Moldován
- Department of Medical Biology, Faculty of Medicine, University of Szeged, Szeged, Hungary
| | - Gábor Torma
- Department of Medical Biology, Faculty of Medicine, University of Szeged, Szeged, Hungary
| | - Tibor Nagy
- Department of Biochemistry and Molecular Biology, Faculty of Medicine, University of Debrecen, Debrecen, Hungary
| | - Ákos Hornyák
- Institute for Veterinary Medical Research, Centre for Agricultural Research, Budapest, Hungary
| | - Zsolt Csabai
- Department of Medical Biology, Faculty of Medicine, University of Szeged, Szeged, Hungary
| | - Gábor Gulyás
- Department of Medical Biology, Faculty of Medicine, University of Szeged, Szeged, Hungary
| | - Miklós Boldogkői
- Department of Medical Biology, Faculty of Medicine, University of Szeged, Szeged, Hungary
| | - Victoria A Jefferson
- Department of Biochemistry & Molecular Biology, Entomology & Plant Pathology, Mississippi State University, Starkville, MS, United States
| | - Zoltán Zádori
- Institute for Veterinary Medical Research, Centre for Agricultural Research, Budapest, Hungary
| | - Florencia Meyer
- Department of Biochemistry & Molecular Biology, Entomology & Plant Pathology, Mississippi State University, Starkville, MS, United States
| | - Zsolt Boldogkői
- Department of Medical Biology, Faculty of Medicine, University of Szeged, Szeged, Hungary
| |
Collapse
|
24
|
Broseus L, Thomas A, Oldfield AJ, Severac D, Dubois E, Ritchie W. TALC: Transcript-level Aware Long-read Correction. Bioinformatics 2021; 36:5000-5006. [PMID: 32910174 DOI: 10.1093/bioinformatics/btaa634] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/03/2020] [Revised: 05/08/2020] [Accepted: 07/09/2020] [Indexed: 02/06/2023] Open
Abstract
MOTIVATION Long-read sequencing technologies are invaluable for determining complex RNA transcript architectures but are error-prone. Numerous 'hybrid correction' algorithms have been developed for genomic data that correct long reads by exploiting the accuracy and depth of short reads sequenced from the same sample. These algorithms are not suited for correcting more complex transcriptome sequencing data. RESULTS We have created a novel reference-free algorithm called Transcript-level Aware Long-Read Correction (TALC) which models changes in RNA expression and isoform representation in a weighted De Bruijn graph to correct long reads from transcriptome studies. We show that transcript-level aware correction by TALC improves the accuracy of the whole spectrum of downstream RNA-seq applications and is thus necessary for transcriptome analyses that use long read technology. AVAILABILITY AND IMPLEMENTATION TALC is implemented in C++ and available at https://github.com/lbroseus/TALC. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Lucile Broseus
- Department of Genome Dynamics, Institut de Génétique Humaine, Centre National de la Recherche Scientifique (CNRS), Université de Montpellier, Montpellier 34396, France
| | - Aubin Thomas
- Department of Genome Dynamics, Institut de Génétique Humaine, Centre National de la Recherche Scientifique (CNRS), Université de Montpellier, Montpellier 34396, France
| | - Andrew J Oldfield
- Department of Genome Dynamics, Institut de Génétique Humaine, Centre National de la Recherche Scientifique (CNRS), Université de Montpellier, Montpellier 34396, France
| | - Dany Severac
- MGX-Montpellier GenomiX, c/o Institut de Génomique Fonctionnelle, Montpellier Cedex 5 34094, France
| | - Emeric Dubois
- MGX-Montpellier GenomiX, c/o Institut de Génomique Fonctionnelle, Montpellier Cedex 5 34094, France
| | - William Ritchie
- Department of Genome Dynamics, Institut de Génétique Humaine, Centre National de la Recherche Scientifique (CNRS), Université de Montpellier, Montpellier 34396, France
| |
Collapse
|
25
|
Tu Z, Shen Y, Wen S, Liu H, Wei L, Li H. A Tissue-Specific Landscape of Alternative Polyadenylation, lncRNAs, TFs, and Gene Co-expression Networks in Liriodendron chinense. FRONTIERS IN PLANT SCIENCE 2021; 12:705321. [PMID: 34367224 PMCID: PMC8343429 DOI: 10.3389/fpls.2021.705321] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/05/2021] [Accepted: 06/28/2021] [Indexed: 05/08/2023]
Abstract
Liriodendron chinense is an economically and ecologically important deciduous tree species. Although the reference genome has been revealed, alternative polyadenylation (APA), transcription factors (TFs), long non-coding RNAs (lncRNAs), and co-expression networks of tissue-specific genes remain incompletely annotated. In this study, we used the bracts, petals, sepals, stamens, pistils, leaves, and shoot apex of L. chinense as materials for hybrid sequencing. On the one hand, we improved the annotation of the genome. We detected 13,139 novel genes, 7,527 lncRNAs, 1,791 TFs, and 6,721 genes with APA sites. On the other hand, we found that tissue-specific genes play a significant role in maintaining tissue characteristics. In total, 2,040 tissue-specific genes were identified, among which 9.2% of tissue-specific genes were affected by APA, and 1,809 tissue-specific genes were represented in seven specific co-expression modules. We also found that bract-specific hub genes were associated plant defense, leaf-specific hub genes were involved in energy metabolism. Moreover, we also found that a stamen-specific hub TF Lchi25777 may be involved in the determination of stamen identity, and a shoot-apex-specific hub TF Lchi05072 may participate in maintaining meristem characteristic. Our study provides a landscape of APA, lncRNAs, TFs, and tissue-specific gene co-expression networks in L. chinense that will improve genome annotation, strengthen our understanding of transcriptome complexity, and drive further research into the regulatory mechanisms of tissue-specific genes.
Collapse
Affiliation(s)
- Zhonghua Tu
- Key Laboratory of Forest Genetics & Biotechnology of Ministry of Education, Nanjing Forestry University, Nanjing, China
- Co-Innovation Center for Sustainable Forestry in Southern China, Nanjing Forestry University, Nanjing, China
| | - Yufang Shen
- Key Laboratory of Forest Genetics & Biotechnology of Ministry of Education, Nanjing Forestry University, Nanjing, China
- Co-Innovation Center for Sustainable Forestry in Southern China, Nanjing Forestry University, Nanjing, China
| | - Shaoying Wen
- Key Laboratory of Forest Genetics & Biotechnology of Ministry of Education, Nanjing Forestry University, Nanjing, China
- Co-Innovation Center for Sustainable Forestry in Southern China, Nanjing Forestry University, Nanjing, China
| | - Huanhuan Liu
- Key Laboratory of Forest Genetics & Biotechnology of Ministry of Education, Nanjing Forestry University, Nanjing, China
- Co-Innovation Center for Sustainable Forestry in Southern China, Nanjing Forestry University, Nanjing, China
| | - Lingmin Wei
- Key Laboratory of Forest Genetics & Biotechnology of Ministry of Education, Nanjing Forestry University, Nanjing, China
- Co-Innovation Center for Sustainable Forestry in Southern China, Nanjing Forestry University, Nanjing, China
| | - Huogen Li
- Key Laboratory of Forest Genetics & Biotechnology of Ministry of Education, Nanjing Forestry University, Nanjing, China
- Co-Innovation Center for Sustainable Forestry in Southern China, Nanjing Forestry University, Nanjing, China
- *Correspondence: Huogen Li,
| |
Collapse
|
26
|
Xu D, Yang H, Zhuo Z, Lu B, Hu J, Yang F. Characterization and analysis of the transcriptome in Opisina arenosella from different developmental stages using single-molecule real-time transcript sequencing and RNA-seq. Int J Biol Macromol 2020; 169:216-227. [PMID: 33340629 DOI: 10.1016/j.ijbiomac.2020.12.098] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2020] [Revised: 09/10/2020] [Accepted: 12/12/2020] [Indexed: 02/06/2023]
Abstract
Opisina arenosella is one of the main pests harming coconut trees. To date, there have been few studies on the molecular genetics, biochemistry and physiology of O. arenosella at the transcriptional level, and there are no available reference genomes. Here, Illumina RNA sequencing combined with PacBio single-molecule real-time analysis was applied to study the transcriptome of this pest at different developmental stages, providing reference data for transcript expression analysis. Twelve samples of O. arenosella from different stages of development were sequenced using Illumina RNA sequencing, and the pooled RNA samples were sequenced with PacBio technology (Iso-Seq). A full-length transcriptome with 41,938 transcripts was captured, and the N50 and N90 lengths were 3543 bp and 1646 bp, respectively. A total of 36,925 transcripts were annotated in public databases, 6493 of which were long noncoding RNAs, while 2510 represented alternative splicing events. There were significant differences in the gene expression profiles at different developmental stages, with high levels of differential gene expression associated with growth, development, carbohydrate metabolism and immunity. This work provides resources and information for the study of the transcriptome and gene function of O. arenosella and provides a valuable foundation for understanding the changes in gene expression during development.
Collapse
Affiliation(s)
- Danping Xu
- College of Life Science, China West Normal University, Nanchong 637002, China
| | - Hongjun Yang
- College of Forestry, Hainan University, Haikou 570228, China
| | - Zhihang Zhuo
- College of Life Science, China West Normal University, Nanchong 637002, China; College of Forestry, Hainan University, Haikou 570228, China.
| | - Baoqian Lu
- Environment and Plant Protection Institute, Chinese Academy of Tropical Agricultural Sciences, Haikou 571101, China
| | - Jiameng Hu
- College of Forestry, Hainan University, Haikou 570228, China
| | - Fan Yang
- College of Forestry, Hainan University, Haikou 570228, China
| |
Collapse
|
27
|
Time-course profiling of bovine alphaherpesvirus 1.1 transcriptome using multiplatform sequencing. Sci Rep 2020; 10:20496. [PMID: 33235226 PMCID: PMC7686369 DOI: 10.1038/s41598-020-77520-1] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2020] [Accepted: 11/04/2020] [Indexed: 12/11/2022] Open
Abstract
Long-read sequencing (LRS) has become a standard approach for transcriptome analysis in recent years. Bovine alphaherpesvirus 1 (BoHV-1) is an important pathogen of cattle worldwide. This study reports the profiling of the dynamic lytic transcriptome of BoHV-1 using two long-read sequencing (LRS) techniques, the Oxford Nanopore Technologies MinION, and the LoopSeq synthetic LRS methods, using multiple library preparation protocols. In this work, we annotated viral mRNAs and non-coding transcripts, and a large number of transcript isoforms, including transcription start and end sites, as well as splice variants of BoHV-1. Our analysis demonstrated an extremely complex pattern of transcriptional overlaps.
Collapse
|
28
|
Yang H, Xu D, Zhuo Z, Hu J, Lu B. Transcriptome and gene expression analysis of Rhynchophorus ferrugineus (Coleoptera: Curculionidae) during developmental stages. PeerJ 2020; 8:e10223. [PMID: 33194414 PMCID: PMC7643551 DOI: 10.7717/peerj.10223] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/05/2020] [Accepted: 09/29/2020] [Indexed: 01/15/2023] Open
Abstract
Background Red palm weevil, Rhynchophorus ferrugineus Olivier, is one of the most destructive pests harming palm trees. However, genomic resources for R. ferrugineus are still lacking, limiting the ability to discover molecular and genetic means of pest control. Methods In this study, PacBio Iso-Seq and Illumina RNA-seq were used to generate transcriptome from three developmental stages of R. ferrugineus (pupa, 7th-instar larva, adult) to increase the understanding of the life cycle and molecular characteristics of the pest. Results Sequencing generated 625,983,256 clean reads, from which 63,801 full-length transcripts were assembled with N50 of 3,547 bp. Expression analyses revealed 8,583 differentially expressed genes (DEGs). Moreover, gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed that these DEGs were mainly related to the peroxisome pathway which associated with metabolic pathways, material transportation and organ tissue formation. In summary, this work provides a valuable basis for further research on the growth and development, gene expression and gene prediction, and pest control of R. ferrugineus.
Collapse
Affiliation(s)
- Hongjun Yang
- College of Life Science, China West Normal University, Nanchong, Sichuan, China.,Key Laboratory of Genetics and Germplasm Innovation of Tropical Special Forest Trees and Ornamental Plants, Ministry of Education, Key Laboratory of Germplasm Resources Biology of Tropical Special Ornamental Plants of Hainan Province, College of Forestry, Hainan University, Haikou, Hainan,China
| | - Danping Xu
- College of Life Science, China West Normal University, Nanchong, Sichuan, China
| | - Zhihang Zhuo
- College of Life Science, China West Normal University, Nanchong, Sichuan, China.,Key Laboratory of Genetics and Germplasm Innovation of Tropical Special Forest Trees and Ornamental Plants, Ministry of Education, Key Laboratory of Germplasm Resources Biology of Tropical Special Ornamental Plants of Hainan Province, College of Forestry, Hainan University, Haikou, Hainan,China.,Key Laboratory of Integrated Pest Management on Crops in South China, Ministry of Agriculture, South China Agricultural University, Guangzhou, Guangdong, China
| | - Jiameng Hu
- Key Laboratory of Genetics and Germplasm Innovation of Tropical Special Forest Trees and Ornamental Plants, Ministry of Education, Key Laboratory of Germplasm Resources Biology of Tropical Special Ornamental Plants of Hainan Province, College of Forestry, Hainan University, Haikou, Hainan,China
| | - Baoqian Lu
- Key Laboratory of Integrated Pest Management on Tropical Crops, Ministry of Agriculture China, Environment and Plant Protection Institute, Chinese Academy of Tropical Agricultural Sciences, Haikou, Hainan, China
| |
Collapse
|
29
|
Chen D, Du Y, Fan X, Zhu Z, Jiang H, Wang J, Fan Y, Chen H, Zhou D, Xiong C, Zheng Y, Xu X, Luo Q, Guo R. Reconstruction and functional annotation of Ascosphaera apis full-length transcriptome utilizing PacBio long reads combined with Illumina short reads. J Invertebr Pathol 2020; 176:107475. [PMID: 32976816 DOI: 10.1016/j.jip.2020.107475] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/25/2019] [Revised: 08/26/2020] [Accepted: 09/16/2020] [Indexed: 01/19/2023]
Abstract
Ascosphaera apis is a widespread fungal pathogen of honeybee larvae that results in chalkbrood disease, leading to heavy losses for the beekeeping industry in China and many other countries. This work was aimed at generating a full-length transcriptome of A. apis using PacBio single-molecule real-time (SMRT) sequencing. Here, more than 23.97 Gb of clean reads was generated from long-read sequencing of A. apis mycelia, including 464,043 circular consensus sequences (CCS) and 394,142 full-length non-chimeric (FLNC) reads. In total, we identified 174,095 high-confidence transcripts covering 5141 known genes with an average length of 2728 bp. We also discovered 2405 genic loci and 11,623 isoforms that have not been annotated yet within the current reference genome. Additionally, 16,049, 10,682, 4520 and 7253 of the discovered transcripts have annotations in the Non-redundant protein (Nr), Clusters of Eukaryotic Orthologous Groups (KOG), Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. Moreover, 1205 long non-coding RNAs (lncRNAs) were identified, which have less exons, shorter exon and intron lengths, shorter transcript lengths, lower GC percent, lower expression levels, and fewer alternative splicing (AS) evens, compared with protein-coding transcripts. A total of 253 members from 17 transcription factor (TF) families were identified from our transcript datasets. Finally, the expression of A. apis isoforms was validated using a molecular approach. Overall, this is the first report of a full-length transcriptome of entomogenous fungi including A. apis. Our data offer a comprehensive set of reference transcripts and hence contributes to improving the genome annotation and transcriptomic study of A. apis.
Collapse
Affiliation(s)
- Dafu Chen
- College of Animal Sciences (College of Bee Science), Fujian Agriculture and Forestry University, 350002 Fuzhou, Fujian, China
| | - Yu Du
- College of Animal Sciences (College of Bee Science), Fujian Agriculture and Forestry University, 350002 Fuzhou, Fujian, China
| | - Xiaoxue Fan
- College of Animal Sciences (College of Bee Science), Fujian Agriculture and Forestry University, 350002 Fuzhou, Fujian, China
| | - Zhiwei Zhu
- College of Animal Sciences (College of Bee Science), Fujian Agriculture and Forestry University, 350002 Fuzhou, Fujian, China
| | - Haibin Jiang
- College of Animal Sciences (College of Bee Science), Fujian Agriculture and Forestry University, 350002 Fuzhou, Fujian, China
| | - Jie Wang
- College of Animal Sciences (College of Bee Science), Fujian Agriculture and Forestry University, 350002 Fuzhou, Fujian, China
| | - Yuanchan Fan
- College of Animal Sciences (College of Bee Science), Fujian Agriculture and Forestry University, 350002 Fuzhou, Fujian, China
| | - Huazhi Chen
- College of Animal Sciences (College of Bee Science), Fujian Agriculture and Forestry University, 350002 Fuzhou, Fujian, China
| | - Dingding Zhou
- College of Animal Sciences (College of Bee Science), Fujian Agriculture and Forestry University, 350002 Fuzhou, Fujian, China
| | - Cuiling Xiong
- College of Animal Sciences (College of Bee Science), Fujian Agriculture and Forestry University, 350002 Fuzhou, Fujian, China
| | - Yanzhen Zheng
- College of Animal Sciences (College of Bee Science), Fujian Agriculture and Forestry University, 350002 Fuzhou, Fujian, China
| | - Xijian Xu
- Jiangxi Province Institute of Apiculture, 330201 Nanchang, Jiangxi, China
| | - Qun Luo
- Jiangxi Province Institute of Apiculture, 330201 Nanchang, Jiangxi, China
| | - Rui Guo
- College of Animal Sciences (College of Bee Science), Fujian Agriculture and Forestry University, 350002 Fuzhou, Fujian, China; Engineering Research Center of Processing and Application of Bee Products of Ministry of Education, Fuzhou 350002, Fujian Province, China.
| |
Collapse
|
30
|
Full-length transcriptome sequencing combined with RNA-seq analysis revealed the immune response of fat greenling (Hexagrammos otakii) to Vibrio harveyi in early infection. Microb Pathog 2020; 149:104527. [PMID: 32980468 DOI: 10.1016/j.micpath.2020.104527] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/24/2020] [Revised: 09/20/2020] [Accepted: 09/21/2020] [Indexed: 01/20/2023]
Abstract
Fat greenling (Hexagrammos otakii) is an important commercial marine fish species cultured in northeast Asia, but its available gene sequences are limited. Vibrio harveyi is a causative agent of vibriosis in fat greenling and also causes severe losses to the aquaculture industry in China. In order to obtain more high-quality transcript information and investigate the early immune response of fat greenling against V. harveyi, the fish were artificially infected with V. harveyi, and five sampling points were set within 48 h. Iso-Seq combined with RNA-Seq were applied in the comprehensive transcriptome analysis of V. harveyi-infected fat greenling. Total 42,225 consensus isoforms were successfully extracted from the result of Iso-Seq, and more than 19,000 ORFs were predicted. In addition, total three modules were identified by WGCNA which significantly positive correlated to the infection time, and the KEGG analysis showed that the immune-related genes in these modules mainly enriched in TLR signaling pathway, NF-κB signaling pathway and Endocytosis. The activation of inflammation and endocytosis was the most significant characteristics of fat greenling immune response during the early infection. Based on the WGCNA, a series of high-degree nodes in the networks were identified as hub genes. The protein structures of cold-inducible RNA-binding protein (CIRBP), poly [ADP-ribose] polymerase 1 (PARP1) and protein arginine N-methyl transferase 1 (PRMT1) were subsequently found to be highly conserved in vertebrate, and the gene expression pattern of CIRBP, PARP1, PRMT1 and a part of TLR/NF-κB pathway-related genes indicated that these proteins might have similar biological functions in regulation of inflammatory response in teleost fish. The results of this study provided the first systematical full-length transcriptome profile of fat greenling and characterized its immune responses in early infection of V. harvey, which will serve as the foundation for further exploring the molecular mechanism of immune defense against bacterial infection in fat greenling.
Collapse
|
31
|
Zhang H, Xu H, Liu H, Pan X, Xu M, Zhang G, He M. PacBio single molecule long-read sequencing provides insight into the complexity and diversity of the Pinctada fucata martensii transcriptome. BMC Genomics 2020; 21:481. [PMID: 32660426 PMCID: PMC7359550 DOI: 10.1186/s12864-020-06894-3] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/13/2019] [Accepted: 07/07/2020] [Indexed: 01/22/2023] Open
Abstract
BACKGROUND The pearl oyster Pinctada fucata martensii is an economically valuable shellfish for seawater pearl production, and production of pearls depends on its growth. To date, the molecular mechanisms of the growth of this species remain poorly understood. The transcriptome sequencing has been considered to understanding of the complexity of mechanisms of the growth of P. f. martensii. The recently released genome sequences of P. f. martensii, as well as emerging Pacific Bioscience (PacBio) single-molecular sequencing technologies, provide an opportunity to thoroughly investigate these molecular mechanisms. RESULTS Herein, the full-length transcriptome was analysed by combining PacBio single-molecule long-read sequencing (PacBio sequencing) and Illumina sequencing. A total of 20.65 Gb of clean data were generated, including 574,561 circular consensus reads, among which 443,944 full-length non-chimeric (FLNC) sequences were identified. Through transcript clustering analysis of FLNC reads, 32,755 consensus isoforms were identified, including 32,095 high-quality consensus sequences. After removing redundant reads, 16,388 transcripts were obtained, and 641 fusion transcripts were derived by performing fusion transcript prediction of consensus sequences. Alternative splicing analysis of the 16,388 transcripts was performed after accounting for redundancy, and 9097 gene loci were detected, including 1607 new gene loci and 14,946 newly discovered transcripts. The original boundary of 11,235 genes on the chromosomes was corrected, 12,025 complete open reading frame sequences and 635 long non-coding RNAs (LncRNAs) were predicted, and functional annotation of 13,482 new transcripts was achieved. Two thousand three hundred eighteen alternative splicing events were detected. A total of 228 differentially expressed transcripts (DETs) were identified between the largest (L) and smallest (S) pearl oysters. Compared with the S, the L showed 99 and 129 significantly up-and down-regulated DETs, respectively. Six of these DETs were further confirmed by quantitative real-time RT-PCR (RT-qPCR) in independent experiment. CONCLUSIONS Our results significantly improve existing gene models and genome annotations, optimise the genome structure, and in-depth understanding of the complexity and diversity of the differential growth patterns of P. f. martensii.
Collapse
Affiliation(s)
- Hua Zhang
- CAS Key Laboratory of Tropical Marine Bio-resources and Ecology, Guangdong Provincial Key Laboratory of Applied Marine Biology, South China Sea Institute of Oceanology, Chinese Academy of Sciences, Guangzhou, 510301, China
| | - Hanzhi Xu
- CAS Key Laboratory of Tropical Marine Bio-resources and Ecology, Guangdong Provincial Key Laboratory of Applied Marine Biology, South China Sea Institute of Oceanology, Chinese Academy of Sciences, Guangzhou, 510301, China.,University of Chinese Academy of Sciences, Beijing, 100049, China
| | - Huiru Liu
- CAS Key Laboratory of Tropical Marine Bio-resources and Ecology, Guangdong Provincial Key Laboratory of Applied Marine Biology, South China Sea Institute of Oceanology, Chinese Academy of Sciences, Guangzhou, 510301, China.,University of Chinese Academy of Sciences, Beijing, 100049, China
| | - Xiaolan Pan
- CAS Key Laboratory of Tropical Marine Bio-resources and Ecology, Guangdong Provincial Key Laboratory of Applied Marine Biology, South China Sea Institute of Oceanology, Chinese Academy of Sciences, Guangzhou, 510301, China.,University of Chinese Academy of Sciences, Beijing, 100049, China
| | - Meng Xu
- CAS Key Laboratory of Tropical Marine Bio-resources and Ecology, Guangdong Provincial Key Laboratory of Applied Marine Biology, South China Sea Institute of Oceanology, Chinese Academy of Sciences, Guangzhou, 510301, China.,University of Chinese Academy of Sciences, Beijing, 100049, China
| | - Gege Zhang
- CAS Key Laboratory of Tropical Marine Bio-resources and Ecology, Guangdong Provincial Key Laboratory of Applied Marine Biology, South China Sea Institute of Oceanology, Chinese Academy of Sciences, Guangzhou, 510301, China.,University of Chinese Academy of Sciences, Beijing, 100049, China
| | - Maoxian He
- CAS Key Laboratory of Tropical Marine Bio-resources and Ecology, Guangdong Provincial Key Laboratory of Applied Marine Biology, South China Sea Institute of Oceanology, Chinese Academy of Sciences, Guangzhou, 510301, China.
| |
Collapse
|
32
|
Yang H, Xu D, Zhuo Z, Hu J, Lu B. SMRT sequencing of the full-length transcriptome of the Rhynchophorus ferrugineus (Coleoptera: Curculionidae). PeerJ 2020; 8:e9133. [PMID: 32509454 PMCID: PMC7246026 DOI: 10.7717/peerj.9133] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2020] [Accepted: 04/14/2020] [Indexed: 12/23/2022] Open
Abstract
Background Red palm weevil Rhynchophorus ferrugineus (Coleoptera: Curculionidae) is one of the most destructive insects for palm trees in the world. However, its genome resources are still in the blank stage, which limits the study of molecular and growth development analysis. Methods In this study, we used PacBio Iso-Seq and Illumina RNA-seq to first generate transcriptome from three developmental stages of R. ferrugineus (pupa, 7th larva, female and male) to increase our understanding of the life cycle and molecular characteristics of R. ferrugineus. Results A total of 63,801 nonredundant full-length transcripts were generated with an average length of 2,964 bp from three developmental stages, including the 7th instar larva, pupa, female adult and male adult. These transcripts showed a high annotation rate in seven public databases, with 54,999 (86.20%) successfully annotated. Meanwhile, 2,184 alternative splicing (AS) events, 2,084 transcription factors (TFs), 66,230 simple sequence repeats (SSR) and 9,618 Long noncoding RNAs (lncRNAs) were identified. In summary, our results provide a new source of full-length transcriptional data and information for the further study of gene expression and genetics in R. ferrugineus.
Collapse
Affiliation(s)
- Hongjun Yang
- Key Laboratory of Genetics and Germplasm Innovation of Tropical Special Forest Trees and Ornamental Plants, Ministry of Education, Key Laboratory of Germplasm Resources Biology of Tropical Special Ornamental Plants of Hainan Province, College of Forestry, Hainan University, Haikou, Hainan, China
| | - Danping Xu
- Sichuan Provincial Key Laboratory of Agricultural Products Processing and Preservative, College of Food Science, Sichuan Agricultural University, Yaan, Sichuan, China
| | - Zhihang Zhuo
- Key Laboratory of Genetics and Germplasm Innovation of Tropical Special Forest Trees and Ornamental Plants, Ministry of Education, Key Laboratory of Germplasm Resources Biology of Tropical Special Ornamental Plants of Hainan Province, College of Forestry, Hainan University, Haikou, Hainan, China.,Key Laboratory of Integrated Pest Management on Crops in South China, Ministry of Agriculture, South China Agricultural University, Guangzhou, Guangdong, China
| | - Jiameng Hu
- Key Laboratory of Genetics and Germplasm Innovation of Tropical Special Forest Trees and Ornamental Plants, Ministry of Education, Key Laboratory of Germplasm Resources Biology of Tropical Special Ornamental Plants of Hainan Province, College of Forestry, Hainan University, Haikou, Hainan, China
| | - Baoqian Lu
- Key Laboratory of Integrated Pest Management on Tropical Crops, Ministry of Agriculture China, Environment and Plant Protection Institute, Chinese Academy of Tropical Agricultural Sciences, Haikou, Hainan, China
| |
Collapse
|
33
|
Characterization and Functional Analysis of Polyadenylation Sites in Fast and Slow Muscles. BIOMED RESEARCH INTERNATIONAL 2020; 2020:2626584. [PMID: 32258109 PMCID: PMC7102456 DOI: 10.1155/2020/2626584] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 09/30/2019] [Revised: 12/01/2019] [Accepted: 01/16/2020] [Indexed: 12/05/2022]
Abstract
Many increasing documents have proved that alternative polyadenylation (APA) events with different polyadenylation sites (PAS) contribute to posttranscriptional regulation. However, little is known about the detailed molecular features of PASs and its role in porcine fast and slow skeletal muscles through microRNAs (miRNAs) and RNA binding proteins (RBPs). In this study, we combined single-molecule real-time sequencing and Illumina RNA-seq datasets to comprehensively analyze polyadenylation in pigs. We identified a total of 10,334 PASs, of which 8734 were characterized by reference genome annotation. 32.86% of PAS-associated genes were determined to have more than one PAS. Further analysis demonstrated that tissue-specific PASs between fast and slow muscles were enriched in skeletal muscle development pathways. In addition, we obtained 1407 target genes regulated by APA events through potential binding 69 miRNAs and 28 RBPs in variable 3′ UTR regions and some are involved in myofiber transformation. Furthermore, the de novo motif search confirmed that the most common usage of canonical motif AAUAAA and three types of PASs may be related to the strength of motifs. In summary, our results provide a useful annotation of PASs for pig transcriptome and suggest that APA may serve as a role in fast and slow muscle development under the regulation of miRNAs and RBPs.
Collapse
|
34
|
Summers KM, Bush SJ, Wu C, Su AI, Muriuki C, Clark EL, Finlayson HA, Eory L, Waddell LA, Talbot R, Archibald AL, Hume DA. Functional Annotation of the Transcriptome of the Pig, Sus scrofa, Based Upon Network Analysis of an RNAseq Transcriptional Atlas. Front Genet 2020; 10:1355. [PMID: 32117413 PMCID: PMC7034361 DOI: 10.3389/fgene.2019.01355] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/27/2019] [Accepted: 12/11/2019] [Indexed: 12/15/2022] Open
Abstract
The domestic pig (Sus scrofa) is both an economically important livestock species and a model for biomedical research. Two highly contiguous pig reference genomes have recently been released. To support functional annotation of the pig genomes and comparative analysis with large human transcriptomic data sets, we aimed to create a pig gene expression atlas. To achieve this objective, we extended a previous approach developed for the chicken. We downloaded RNAseq data sets from public repositories, down-sampled to a common depth, and quantified expression against a reference transcriptome using the mRNA quantitation tool, Kallisto. We then used the network analysis tool Graphia to identify clusters of transcripts that were coexpressed across the merged data set. Consistent with the principle of guilt-by-association, we identified coexpression clusters that were highly tissue or cell-type restricted and contained transcription factors that have previously been implicated in lineage determination. Other clusters were enriched for transcripts associated with biological processes, such as the cell cycle and oxidative phosphorylation. The same approach was used to identify coexpression clusters within RNAseq data from multiple individual liver and brain samples, highlighting cell type, process, and region-specific gene expression. Evidence of conserved expression can add confidence to assignment of orthology between pig and human genes. Many transcripts currently identified as novel genes with ENSSSCG or LOC IDs were found to be coexpressed with annotated neighbouring transcripts in the same orientation, indicating they may be products of the same transcriptional unit. The meta-analytic approach to utilising public RNAseq data is extendable to include new data sets and new species and provides a framework to support the Functional Annotation of Animals Genomes (FAANG) initiative.
Collapse
Affiliation(s)
- Kim M. Summers
- Mater Research Institute-University of Queensland, Translational Research Institute, Woolloongabba, QLD, Australia
| | - Stephen J. Bush
- Nuffield Department of Clinical Medicine, University of Oxford, Oxford, United Kingdom
| | - Chunlei Wu
- Department of Integrative Structural and Computational Biology, The Scripps Research Institute, La Jolla, CA, United States
| | - Andrew I. Su
- Department of Integrative Structural and Computational Biology, The Scripps Research Institute, La Jolla, CA, United States
| | - Charity Muriuki
- The Roslin Institute, University of Edinburgh, Midlothian, United Kingdom
| | - Emily L. Clark
- The Roslin Institute, University of Edinburgh, Midlothian, United Kingdom
| | | | - Lel Eory
- The Roslin Institute, University of Edinburgh, Midlothian, United Kingdom
| | - Lindsey A. Waddell
- The Roslin Institute, University of Edinburgh, Midlothian, United Kingdom
| | - Richard Talbot
- The Roslin Institute, University of Edinburgh, Midlothian, United Kingdom
| | - Alan L. Archibald
- The Roslin Institute, University of Edinburgh, Midlothian, United Kingdom
| | - David A. Hume
- Mater Research Institute-University of Queensland, Translational Research Institute, Woolloongabba, QLD, Australia
| |
Collapse
|
35
|
Bick JT, Zeng S, Robinson MD, Ulbrich SE, Bauersachs S. Mammalian Annotation Database for improved annotation and functional classification of Omics datasets from less well-annotated organisms. DATABASE-THE JOURNAL OF BIOLOGICAL DATABASES AND CURATION 2020; 2019:5539597. [PMID: 31353404 PMCID: PMC6661403 DOI: 10.1093/database/baz086] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 08/03/2018] [Revised: 03/08/2019] [Accepted: 06/05/2019] [Indexed: 02/02/2023]
Abstract
Next-generation sequencing technologies and the availability of an increasing number of mammalian and other genomes allow gene expression studies, particularly RNA sequencing, in many non-model organisms. However, incomplete genome annotation and assignments of genes to functional annotation databases can lead to a substantial loss of information in downstream data analysis. To overcome this, we developed Mammalian Annotation Database tool (MAdb, https://madb.ethz.ch) to conveniently provide homologous gene information for selected mammalian species. The assignment between species is performed in three steps: (i) matching official gene symbols, (ii) using ortholog information contained in Ensembl Compara and (iii) pairwise BLAST comparisons of all transcripts. In addition, we developed a new tool (AnnOverlappeR) for the reliable assignment of the National Center for Biotechnology Information (NCBI) and Ensembl gene IDs. The gene lists translated to gene IDs of well-annotated species such as a human can be used for improved functional annotation with relevant tools based on Gene Ontology and molecular pathway information. We tested the MAdb on a published RNA-seq data set for the pig and showed clearly improved overrepresentation analysis results based on the assigned human homologous gene identifiers. Using the MAdb revealed a similar list of human homologous genes and functional annotation results regardless of whether starting with gene IDs from NCBI or Ensembl. The MAdb database is accessible via a web interface and a Galaxy application.
Collapse
Affiliation(s)
- Jochen T Bick
- Animal Physiology, Institute of Agricultural Sciences, ETH Zurich, Zurich, Switzerland
| | - Shuqin Zeng
- Animal Physiology, Institute of Agricultural Sciences, ETH Zurich, Zurich, Switzerland.,Genetics and Functional Genomics, Vetsuisse Faculty Zurich, University of Zurich, Zurich, Switzerland
| | - Mark D Robinson
- Institute of Molecular Life Sciences and SIB Swiss Institute of Bioinformatics, University of Zurich, Zurich, Switzerland
| | - Susanne E Ulbrich
- Animal Physiology, Institute of Agricultural Sciences, ETH Zurich, Zurich, Switzerland
| | - Stefan Bauersachs
- Animal Physiology, Institute of Agricultural Sciences, ETH Zurich, Zurich, Switzerland.,Genetics and Functional Genomics, Vetsuisse Faculty Zurich, University of Zurich, Zurich, Switzerland
| |
Collapse
|
36
|
Feng X, Jia Y, Zhu R, Chen K, Chen Y. Characterization and analysis of the transcriptome in Gymnocypris selincuoensis on the Qinghai-Tibetan Plateau using single-molecule long-read sequencing and RNA-seq. DNA Res 2020; 26:353-363. [PMID: 31274170 PMCID: PMC6704404 DOI: 10.1093/dnares/dsz014] [Citation(s) in RCA: 25] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2019] [Accepted: 06/11/2019] [Indexed: 12/26/2022] Open
Abstract
The lakes on the Qinghai-Tibet Plateau (QTP) are the largest and highest lake group in the world. Gymnocypris selincuoensis is the only cyprinid fish living in lake Selincuo, the largest lake on QTP. However, its genetic resource is still blank, limiting studies on molecular and genetic analysis. In this study, the transcriptome of G. selincuoensis was first generated by using PacBio Iso-Seq and Illumina RNA-seq. A full-length (FL) transcriptome with 75,435 transcripts was obtained by Iso-Seq with N50 length of 3,870 bp. Among all transcripts, 75,016 were annotated to public databases, 64,710 contain complete open reading frames and 2,811 were long non-coding RNAs. Based on all- vs.-all BLAST, 2,069 alternative splicing events were detected, and 80% of them were validated by reverse transcription polymerase chain reaction (RT-PCR). Tissue gene expression atlas showed that the number of detected expressed transcripts ranged from 37,397 in brain to 19,914 in muscle, with 10,488 transcripts detected in all seven tissues. Comparative genomic analysis with other cyprinid fishes identified 77 orthologous genes with potential positive selection (Ka/Ks > 0.3). A total of 56,696 perfect simple sequence repeats were identified from FL transcripts. Our results provide valuable genetic resources for further studies on adaptive evolution, gene expression and population genetics in G. selincuoensis and other congeneric fishes.
Collapse
Affiliation(s)
- Xiu Feng
- The Key Laboratory of Aquatic Biodiversity and Conservation, Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan, China
| | - Yintao Jia
- The Key Laboratory of Aquatic Biodiversity and Conservation, Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan, China
| | - Ren Zhu
- The Key Laboratory of Aquatic Biodiversity and Conservation, Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan, China
| | - Kang Chen
- The Key Laboratory of Aquatic Biodiversity and Conservation, Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan, China.,University of Chinese Academy of Sciences, Beijing, China
| | - Yifeng Chen
- The Key Laboratory of Aquatic Biodiversity and Conservation, Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan, China
| |
Collapse
|
37
|
Lou F, Song N, Han Z, Gao T. Single-molecule real-time (SMRT) sequencing facilitates Tachypleus tridentatus genome annotation. Int J Biol Macromol 2020; 147:89-97. [PMID: 31923512 DOI: 10.1016/j.ijbiomac.2020.01.029] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/12/2019] [Revised: 01/04/2020] [Accepted: 01/04/2020] [Indexed: 12/19/2022]
Abstract
Tachypleus tridentatus is a keystone species in marine ecosystems. Its hemolymph also provides the limulus amebocyte lysate (LAL) for detection of bacterial endotoxin in human medical service. Here we combined SMRT sequencing and Illumina RNA-seq to characterize the novel isoforms, novel genetic loci, fusion isoforms formation and transcriptome structure and further to unveil the transcriptome complexity of T. tridentatus. We identified 26,705 non-redundancy isoforms form 10,919 genetic loci, including 25,713 novel isoforms, 2403 novel genes and 170 fusion isoforms. In addition, 1578 novel genes and 23,172 novel isoforms were annotated in the NR, Pfam, KOG, COG, eggNOG, Swiss-Prot, KEGG and GO databases. Meanwhile, we have obtained 4671 gene family clustering based on genetic loci. Furthermore, there are 17,296, 4887, 1054, and 1435 APAs, AS events, lncRNAs, and TFs were identified in the T. tridentatus long-read transcriptome and the target genes of 1054 lncRNA sequences were also predicted. Overall, our work firstly provided the long-read transcriptome and these data are very necessary to improve the annotation information of T. tridentatus genome and optimize the boundaries of 12,342 original reference annotated genes. Furthermore, these information are a potential resource to study LAL secretion mechanisms in T. tridentatus.
Collapse
Affiliation(s)
- Fangrui Lou
- Fishery College, Ocean University of China, Qingdao, Shandong 266003, China; Fishery College, Zhejiang Ocean University, Zhoushan, Zhejiang 316022, China
| | - Na Song
- Fishery College, Ocean University of China, Qingdao, Shandong 266003, China
| | - Zhiqiang Han
- Fishery College, Zhejiang Ocean University, Zhoushan, Zhejiang 316022, China.
| | - Tianxiang Gao
- Fishery College, Zhejiang Ocean University, Zhoushan, Zhejiang 316022, China.
| |
Collapse
|
38
|
Tian Y, Wen H, Qi X, Zhang X, Liu S, Li B, Sun Y, Li J, He F, Yang W, Li Y. Characterization of Full-Length Transcriptome Sequences and Splice Variants of Lateolabrax maculatus by Single-Molecule Long-Read Sequencing and Their Involvement in Salinity Regulation. Front Genet 2019; 10:1126. [PMID: 31803231 PMCID: PMC6873903 DOI: 10.3389/fgene.2019.01126] [Citation(s) in RCA: 23] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/20/2019] [Accepted: 10/17/2019] [Indexed: 12/17/2022] Open
Abstract
Transcriptome complexity plays crucial roles in regulating the biological functions of eukaryotes. Except for functional genes, alternative splicing and fusion transcripts produce a vast expansion of transcriptome diversity. In this study, we applied PacBio single-molecule long-read sequencing technology to unveil the whole transcriptome landscape of Lateolabrax maculatus. We obtained 28,809 high-quality non-redundant transcripts, including 18,280 novel isoforms covering 8,961 annotated gene loci within the current reference genome and 3,172 novel isoforms. A total of 10,249 AS events were detected, and intron retention was the predominant AS event. In addition, 1,359 alternative polyadenylation events, 3,112 lncRNAs, 29,609 SSRs, 365 fusion transcripts, and 1,194 transcription factors were identified in this study. Furthermore, we performed RNA-Seq analysis combined with Iso-Seq results to investigate salinity regulation mechanism at the transcripts level. A total of 518 transcripts were differentially expressed, which were further divided into 8 functional groups. Notably, transcripts from the same genes exhibited similar or opposite expression patterns. Our study provides a comprehensive view of the transcriptome complexity in L. maculatus, which significantly improves current gene models. Moreover, the diversity of the expression patterns of transcripts may enhance the understanding of salinity regulatory mechanism in L. maculatus and other euryhaline teleosts.
Collapse
Affiliation(s)
- Yuan Tian
- Key Laboratory of Mariculture, Ministry of Education, Ocean University of China, Qingdao, China
| | - Haishen Wen
- Key Laboratory of Mariculture, Ministry of Education, Ocean University of China, Qingdao, China
| | - Xin Qi
- Key Laboratory of Mariculture, Ministry of Education, Ocean University of China, Qingdao, China
| | - Xiaoyan Zhang
- Key Laboratory of Mariculture, Ministry of Education, Ocean University of China, Qingdao, China
| | - Shikai Liu
- Key Laboratory of Mariculture, Ministry of Education, Ocean University of China, Qingdao, China
| | - Bingyu Li
- Key Laboratory of Mariculture, Ministry of Education, Ocean University of China, Qingdao, China
| | - Yalong Sun
- Key Laboratory of Mariculture, Ministry of Education, Ocean University of China, Qingdao, China
| | - Jifang Li
- Key Laboratory of Mariculture, Ministry of Education, Ocean University of China, Qingdao, China
| | - Feng He
- Key Laboratory of Mariculture, Ministry of Education, Ocean University of China, Qingdao, China
| | - Wenzhao Yang
- Key Laboratory of Mariculture, Ministry of Education, Ocean University of China, Qingdao, China
| | - Yun Li
- Key Laboratory of Mariculture, Ministry of Education, Ocean University of China, Qingdao, China
| |
Collapse
|
39
|
Tombácz D, Moldován N, Balázs Z, Gulyás G, Csabai Z, Boldogkői M, Snyder M, Boldogkői Z. Multiple Long-Read Sequencing Survey of Herpes Simplex Virus Dynamic Transcriptome. Front Genet 2019; 10:834. [PMID: 31608102 PMCID: PMC6769088 DOI: 10.3389/fgene.2019.00834] [Citation(s) in RCA: 28] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2019] [Accepted: 08/13/2019] [Indexed: 12/12/2022] Open
Abstract
Long-read sequencing (LRS) has become increasingly important in RNA research due to its strength in resolving complex transcriptomic architectures. In this regard, currently two LRS platforms have demonstrated adequate performance: the Single Molecule Real-Time Sequencing by Pacific Biosciences (PacBio) and the nanopore sequencing by Oxford Nanopore Technologies (ONT). Even though these techniques produce lower coverage and are more error prone than short-read sequencing, they continue to be more successful in identifying polycistronic RNAs, transcript isoforms including splice and transcript end variants, as well as transcript overlaps. Recent reports have successfully applied LRS for the investigation of the transcriptome of viruses belonging to various families. These studies have substantially increased the number of previously known viral RNA molecules. In this work, we used the Sequel and MinION technique from PacBio and ONT, respectively, to characterize the lytic transcriptome of the herpes simplex virus type 1 (HSV-1). In most samples, we analyzed the poly(A) fraction of the transcriptome, but we also performed random oligonucleotide-based sequencing. Besides cDNA sequencing, we also carried out native RNA sequencing. Our investigations identified more than 2,300 previously undetected transcripts, including coding, and non-coding RNAs, multi-splice transcripts, as well as polycistronic and complex transcripts. Furthermore, we found previously unsubstantiated transcriptional start sites, polyadenylation sites, and splice sites. A large number of novel transcriptional overlaps were also detected. Random-primed sequencing revealed that each convergent gene pair produces non-polyadenylated read-through RNAs overlapping the partner genes. Furthermore, we identified novel replication-associated transcripts overlapping the HSV-1 replication origins, and novel LAT variants with very long 5' regions, which are co-terminal with the LAT-0.7kb transcript. Overall, our results demonstrated that the HSV-1 transcripts form an extremely complex pattern of overlaps, and that entire viral genome is transcriptionally active. In most viral genes, if not in all, both DNA strands are expressed.
Collapse
Affiliation(s)
- Dóra Tombácz
- Department of Medical Biology, Faculty of Medicine, University of Szeged, Szeged, Hungary
| | - Norbert Moldován
- Department of Medical Biology, Faculty of Medicine, University of Szeged, Szeged, Hungary
| | - Zsolt Balázs
- Department of Medical Biology, Faculty of Medicine, University of Szeged, Szeged, Hungary
| | - Gábor Gulyás
- Department of Medical Biology, Faculty of Medicine, University of Szeged, Szeged, Hungary
| | - Zsolt Csabai
- Department of Medical Biology, Faculty of Medicine, University of Szeged, Szeged, Hungary
| | - Miklós Boldogkői
- Department of Medical Biology, Faculty of Medicine, University of Szeged, Szeged, Hungary
| | - Michael Snyder
- Department of Genetics, School of Medicine, Stanford University, Stanford, CA, United States
| | - Zsolt Boldogkői
- Department of Medical Biology, Faculty of Medicine, University of Szeged, Szeged, Hungary
| |
Collapse
|
40
|
Mattucci F, Galaverni M, Lyons LA, Alves PC, Randi E, Velli E, Pagani L, Caniglia R. Genomic approaches to identify hybrids and estimate admixture times in European wildcat populations. Sci Rep 2019; 9:11612. [PMID: 31406125 PMCID: PMC6691104 DOI: 10.1038/s41598-019-48002-w] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2018] [Accepted: 07/25/2019] [Indexed: 12/22/2022] Open
Abstract
The survival of indigenous European wildcat (Felis silvestris silvestris) populations can be locally threatened by introgressive hybridization with free-ranging domestic cats. Identifying pure wildcats and investigating the ancestry of admixed individuals becomes thus a conservation priority. We analyzed 63k cat Single Nucleotide Polymorphisms (SNPs) with multivariate, Bayesian and gene-search tools to better evaluate admixture levels between domestic and wild cats collected in Europe, timing and ancestry proportions of their hybrids and backcrosses, and track the origin (wild or domestic) of the genomic blocks carried by admixed cats, also looking for possible deviations from neutrality in their inheritance patterns. Small domestic ancestry blocks were detected in the genomes of most admixed cats, which likely originated from hybridization events occurring from 6 to 22 generations in the past. We identified about 1,900 outlier coding genes with excess of wild or domestic ancestry compared to random expectations in the admixed individuals. More than 600 outlier genes were significantly enriched for Gene Ontology (GO) categories mainly related to social behavior, functional and metabolic adaptive processes (wild-like genes), involved in cognition and neural crest development (domestic-like genes), or associated with immune system functions and lipid metabolism (parental-like genes). These kinds of genomic ancestry analyses could be reliably applied to unravel the admixture dynamics in European wildcats, as well as in other hybridizing populations, in order to design more efficient conservation plans.
Collapse
Affiliation(s)
- Federica Mattucci
- Area per la Genetica della Conservazione (BIO-CGE), Istituto Superiore per la Protezione e la Ricerca Ambientale (ISPRA), Ozzano dell'Emilia, Italy.
| | | | - Leslie A Lyons
- Department of Veterinary Medicine and Surgery, College of Veterinary Medicine, University of Missouri, Columbia, USA
| | - Paulo C Alves
- Centro de Investigação em Biodiversidade e Recursos Genéticos (CIBIO), InBio - Laboratório Associado, Campus Agrário de Vairão, Vairão, Portugal
- Departamento de Biologia, Faculdade de Ciências, Universidade do Porto, Porto, Portugal
- Wildlife Biology Program, Department of Ecosystem and Conservation Sciences, University of Montana, Missoula, USA
| | - Ettore Randi
- Department of Biological, Geological and Environmental Sciences, University of Bologna, Bologna, Italy
- Department of Chemistry and Bioscience, Faculty of Engineering and Science, University of Aalborg, Aalborg, Denmark
| | - Edoardo Velli
- Area per la Genetica della Conservazione (BIO-CGE), Istituto Superiore per la Protezione e la Ricerca Ambientale (ISPRA), Ozzano dell'Emilia, Italy
| | - Luca Pagani
- Dipartimento di Biologia, Università degli Studi di Padova, Padua, Italy
- Estonian Biocentre, Institute of Genomics, University of Tartu, Tartu, Estonia
| | - Romolo Caniglia
- Area per la Genetica della Conservazione (BIO-CGE), Istituto Superiore per la Protezione e la Ricerca Ambientale (ISPRA), Ozzano dell'Emilia, Italy
| |
Collapse
|
41
|
Wang L, Jiang X, Wang L, Wang W, Fu C, Yan X, Geng X. A survey of transcriptome complexity using PacBio single-molecule real-time analysis combined with Illumina RNA sequencing for a better understanding of ricinoleic acid biosynthesis in Ricinus communis. BMC Genomics 2019; 20:456. [PMID: 31170917 PMCID: PMC6555039 DOI: 10.1186/s12864-019-5832-9] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/11/2018] [Accepted: 05/23/2019] [Indexed: 02/07/2023] Open
Abstract
Background Ricinus communis is a highly economically valuable oil crop plant from the spurge family, Euphorbiaceae. However, the available reference genomes are incomplete and to date studies on ricinoleic acid biosynthesis at the transcriptional level are limited. Results In this study, we combined PacBio single-molecule long read isoform and Illumina RNA sequencing to identify the alternative splicing (AS) events, novel isoforms, fusion genes, long non-coding RNAs (lncRNAs) and alternative polyadenylation (APA) sites to unveil the transcriptomic complexity of castor beans and identify critical genes related to ricinoleic acid biosynthesis. Here, we identified 11,285 AS-variants distributed in 21,448 novel genes and detected 520 fusion genes, 320 lncRNAs and 9511 (APA-sites). Furthermore, a total of 6067, 5983 and 4058 differentially expressed genes between developing beans of the R. communis lines 349 and 1115 with extremely different oil content were identified at 7, 14 and 21 days after flowering, respectively. Specifically, 14, 18 and 11 DEGs were annotated encoding key enzymes related to ricinoleic acid biosynthesis reflecting the higher castor oil content of 1115 compared than 349. Quantitative real-time RT-PCR further validated fifteen of these DEGs at three-time points. Conclusion Our results significantly improved the existed gene models of R. communis, and a putative model of key genes was built to show the differences between strains 349 and 1115, illustrating the molecular mechanism of castor oil biosynthesis. A multi-transcriptome database and candidate genes were provided to further improve the level of ricinoleic acid in transgenic crops. Electronic supplementary material The online version of this article (10.1186/s12864-019-5832-9) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Lijun Wang
- Oil Crops Research Institute of the Chinese Academy of Agricultural Sciences/Key Laboratory of Biology and Genetic, Improvement of Oil Crops, Ministry of Agriculture, Wuhan, China
| | - Xiaoling Jiang
- College of Life Science and Technology, Henan Institute of Science and Technology/Collaborative Innovation Center of Modern Biological Breeding, Xinxiang, China
| | - Lei Wang
- Oil Crops Research Institute of the Chinese Academy of Agricultural Sciences/Key Laboratory of Biology and Genetic, Improvement of Oil Crops, Ministry of Agriculture, Wuhan, China
| | - Wei Wang
- Oil Crops Research Institute of the Chinese Academy of Agricultural Sciences/Key Laboratory of Biology and Genetic, Improvement of Oil Crops, Ministry of Agriculture, Wuhan, China
| | - Chunling Fu
- Oil Crops Research Institute of the Chinese Academy of Agricultural Sciences/Key Laboratory of Biology and Genetic, Improvement of Oil Crops, Ministry of Agriculture, Wuhan, China
| | - Xingchu Yan
- Oil Crops Research Institute of the Chinese Academy of Agricultural Sciences/Key Laboratory of Biology and Genetic, Improvement of Oil Crops, Ministry of Agriculture, Wuhan, China.
| | - Xinxin Geng
- Applied Biotechnology Center, Wuhan University of Bioengineering, Wuhan, China.
| |
Collapse
|
42
|
Beiki H, Liu H, Huang J, Manchanda N, Nonneman D, Smith TPL, Reecy JM, Tuggle CK. Improved annotation of the domestic pig genome through integration of Iso-Seq and RNA-seq data. BMC Genomics 2019; 20:344. [PMID: 31064321 PMCID: PMC6505119 DOI: 10.1186/s12864-019-5709-y] [Citation(s) in RCA: 39] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2018] [Accepted: 04/17/2019] [Indexed: 01/07/2023] Open
Abstract
BACKGROUND Our understanding of the pig transcriptome is limited. RNA transcript diversity among nine tissues was assessed using poly(A) selected single-molecule long-read isoform sequencing (Iso-seq) and Illumina RNA sequencing (RNA-seq) from a single White cross-bred pig. RESULTS Across tissues, a total of 67,746 unique transcripts were observed, including 60.5% predicted protein-coding, 36.2% long non-coding RNA and 3.3% nonsense-mediated decay transcripts. On average, 90% of the splice junctions were supported by RNA-seq within tissue. A large proportion (80%) represented novel transcripts, mostly produced by known protein-coding genes (70%), while 17% corresponded to novel genes. On average, four transcripts per known gene (tpg) were identified; an increase over current EBI (1.9 tpg) and NCBI (2.9 tpg) annotations and closer to the number reported in human genome (4.2 tpg). Our new pig genome annotation extended more than 6000 known gene borders (5' end extension, 3' end extension, or both) compared to EBI or NCBI annotations. We validated a large proportion of these extensions by independent pig poly(A) selected 3'-RNA-seq data, or human FANTOM5 Cap Analysis of Gene Expression data. Further, we detected 10,465 novel genes (81% non-coding) not reported in current pig genome annotations. More than 80% of these novel genes had transcripts detected in > 1 tissue. In addition, more than 80% of novel intergenic genes with at least one transcript detected in liver tissue had H3K4me3 or H3K36me3 peaks mapping to their promoter and gene body, respectively, in independent liver chromatin immunoprecipitation data. CONCLUSIONS These validated results show significant improvement over current pig genome annotations.
Collapse
Affiliation(s)
- H Beiki
- Department of Animal Science, Iowa State University, 2255 Kildee Hall, Ames, IA, 50011, USA
| | - H Liu
- Department of Animal Science, Iowa State University, 2255 Kildee Hall, Ames, IA, 50011, USA
| | - J Huang
- Department of Animal Science, Iowa State University, 2255 Kildee Hall, Ames, IA, 50011, USA.,College of Animal Science and Technology, Jiangxi Agricultural University, Nanchang, Jiangxi, People's Republic of China
| | - N Manchanda
- Department of Ecology, Evolution, and Organismal Biology, Iowa State University, 819 Wallace Road, Ames, IA, 50011, USA
| | - D Nonneman
- USDA, ARS, U.S. Meat Animal Research Center, Clay Center, NE, 68933, USA
| | - T P L Smith
- USDA, ARS, U.S. Meat Animal Research Center, Clay Center, NE, 68933, USA
| | - J M Reecy
- Department of Animal Science, Iowa State University, 2255 Kildee Hall, Ames, IA, 50011, USA
| | - C K Tuggle
- Department of Animal Science, Iowa State University, 2255 Kildee Hall, Ames, IA, 50011, USA.
| |
Collapse
|