Finding the genes in genomic DNA

Genomic-transcriptomic analysis identifies the Syrian hamster as a superior animal model for human diseases

BACKGROUND

The Syrian hamster (Mesocricetus auratus) has shown promise as a human diseases model, recapitulating features of different human diseases including COVID-19. However, the landscape of its genome and transcriptome has not been systematically dissected, restricting its potential applications.

RESULTS

Here we provide a complete analysis of the genome and transcriptome of the Syrian hamster and found that its lineage diverged from that of the Chinese hamster (Cricetulus griseus) around 29.4 million years ago. 21,387 protein-coding genes were identified, with 90.03% of the 2.56G base pair sequence being anchored to 22 chromosomes. Further comparison of the transcriptomes from 15 tissues of the Syrian hamster revealed that the Syrian hamster shares a pattern of alternative splicing modes more similar to humans, compared to rats and mice. An integrated genomic-transcriptomic analysis revealed that the Syrian hamster also has genetic and biological advantages as a superior animal model for cardiovascular diseases. Strikingly, several genes involved in SARS-COV-2 infection, including ACE2, present a higher homology with humans compared to other rodents and show the same function as their human counterparts.

CONCLUSION

The detailed molecular characterisation of the Syrian hamster in the present study opens a wealth of fundamental resources from this small rodent for future research into human disease pathology and treatment.

Integrative multi-omics analysis reveals the contribution of neoVTX genes to venom diversity of Synanceia verrucosa

BACKGROUND

Animal venom systems are considered as valuable model for investigating the molecular mechanisms underlying phenotypic evolution. Stonefish are the most venomous and dangerous fish because of severe human envenomation and occasionally fatalities, whereas the genomic background of their venom has not been fully explored compared with that in other venomous animals.

RESULTS

In this study, we followed modern venomic pipelines to decode the Synanceia verrucosa venom components. A catalog of 478 toxin genes was annotated based on our assembled chromosome-level genome. Integrative analysis of the high-quality genome, the transcriptome of the venom gland, and the proteome of crude venom revealed mechanisms underlying the venom complexity in S. verrucosa. Six tandem-duplicated neoVTX subunit genes were identified as the major source for the neoVTX protein production. Further isoform sequencing revealed massive alternative splicing events with a total of 411 isoforms demonstrated by the six genes, which further contributed to the venom diversity. We then characterized 12 dominantly expressed toxin genes in the venom gland, and 11 of which were evidenced to produce the venom protein components, with the neoVTX proteins as the most abundant. Other major venom proteins included a presumed CRVP, Kuntiz-type serine protease inhibitor, calglandulin protein, and hyaluronidase. Besides, a few of highly abundant non-toxin proteins were also characterized and they were hypothesized to function in housekeeping or hemostasis maintaining roles in the venom gland. Notably, gastrotropin like non-toxin proteins were the second highest abundant proteins in the venom, which have not been reported in other venomous animals and contribute to the unique venom properties of S. verrucosa.

CONCLUSIONS

The results identified the major venom composition of S. verrucosa, and highlighted the contribution of neoVTX genes to the diversity of venom composition through tandem-duplication and alternative splicing. The diverse neoVTX proteins in the venom as lethal particles are important for understanding the adaptive evolution of S. verrucosa. Further functional studies are encouraged to exploit the venom components of S. verrucosa for pharmaceutical innovation.

Cancer detection with various classification models: A comprehensive feature analysis using HMM to extract a nucleotide pattern

This work presents a novel feature extraction method for identifying complex patterns in genomic sequences by employing the Hidden Markov Model (HMM). In this study, we use HMM to identify gene nucleotide patterns that are specific to malignant and non-malignant cells. Crucial genetic components DNA and RNA are involved in many biological processes that impact both healthy and malignant cells. Early patient identification is essential to successful cancer diagnosis and therapy. Varying nucleotide patterns indicate different cellular responses, which are important to understanding the molecular causes of cancer and associated disorders. We present a detailed study of nucleotide patterns in whole raw nucleotide sequences with variations in both protein sequence (CDS) and non-protein sequence (NCDS) in both malignant and non-malignant cells. Nucleotide prediction has been achieved while computational expenses are reduced by using the proposed HMM for feature extraction and selection. The classification models implemented in this work for cancer detection are Gradient-Boosted Decision Trees (GBDT), Random Forests (RF), Decision Trees (DT), and Support Vector Machines (SVM) with kernels. The suggested classification model's accuracy and 10-fold cross-validation have been validated via comprehensive case studies. The results reveal that DT and ensemble learning techniques significantly differentiate between malignant and non-malignant DNA sequences. SVM with suitable kernels improves cancer detection accuracy significantly. Combining feature reduction approaches with nucleotide pattern classifiers based on Hidden Markov models improves performance and ensures reliable cancer detection.

Microprotein-encoding RNA regulation in cells treated with pro-inflammatory and pro-fibrotic stimuli

BACKGROUND

Recent analysis of the human proteome via proteogenomics and ribosome profiling of the transcriptome revealed the existence of thousands of previously unannotated microprotein-coding small open reading frames (smORFs). Most functional microproteins were chosen for characterization because of their evolutionary conservation. However, one example of a non-conserved immunomodulatory microprotein in mice suggests that strict sequence conservation misses some intriguing microproteins.

RESULTS

We examine the ability of gene regulation to identify human microproteins with potential roles in inflammation or fibrosis of the intestine. To do this, we collected ribosome profiling data of intestinal cell lines and peripheral blood mononuclear cells and used gene expression of microprotein-encoding transcripts to identify strongly regulated microproteins, including several examples of microproteins that are only conserved with primates.

CONCLUSION

This approach reveals a number of new microproteins worthy of additional functional characterization and provides a dataset that can be queried in different ways to find additional gut microproteins of interest.

A Kernelized Classification Approach for Cancer Recognition Using Markovian Analysis of DNA Structure Patterns as Feature Mining

Nucleotide-based molecules called DNA and RNA are essential for several biological processes that affect both normal and cancerous cells. They contain the critical genetic material needed for normal cell growth and functioning. The DNA structure patterns that make up the genetic code affect cells' growth, behavior, and control. Different DNA structure patterns indicate different physiological effects in the cell. Knowledge of these patterns is necessary to identify the molecular origins of cancer and other disorders. Analyzing these patterns can help in the early detection of diseases, which is essential for the effectiveness of cancer research and therapy. The novelty of this study is to examine the patterns of dinucleotide structure in many genomic regions, including the non-coding region sequence (N-CDS), coding region sequence (CDS), and whole raw DNA sequence (W.R. sequence). It provides an in-depth discussion of dinucleotide patterns related to these diverse genetic environments and contains malignant and non-malignant DNA sequences. The Markovian modeling that predicts dinucleotide probabilities also reduces feature complexity and minimizes computational costs compared to the approaches of Kernelized Logistic Regression (KLR) and Support Vector Machine (SVM). This technique is effectively evaluated in essential case studies, as indicated by accuracy metrics and 10-fold cross-validation. The classifier and feature reduction, which are generated by Markovian probability, operate well together and can help predict cancer. Our findings successfully distinguish DNA sequences related to cancer from those diagnostics of non-cancerous diseases by analyzing the W.R. DNA sequence as well as its CDS and N-CDS regions.

Pangenomics of flax fungal parasite Fusarium oxysporum f. sp. lini

To assess the genomic diversity of Fusarium oxysporum f. sp. lini strains and compile a comprehensive gene repertoire, we constructed a pangenome using 13 isolates from four different clonal lineages, each exhibiting distinct levels of virulence. Syntenic analyses of two selected genomes revealed significant chromosomal rearrangements unique to each genome. A comprehensive examination of both core and accessory pangenome content and diversity points at an open genome state. Additionally, Gene Ontology (GO) enrichment analysis indicated that non-core pangenome genes are associated with pathogen recognition and immune signaling. Furthermore, the Folini pansecterome, encompassing secreted proteins critical for fungal pathogenicity, primarily consists of three functional classes: effector proteins, CAZYmes, and proteases. These three classes account for approximately 3.5% of the pangenome. Each functional class within the pansecterome was meticulously annotated and characterized with respect to pangenome category distribution, PFAM domain frequency, and strain virulence assessment. This analysis revealed that highly virulent isolates have specific types of PFAM domains that are exclusive to them. Upon examining the repertoire of SIX genes known for virulence in other formae speciales, it was found that all isolates had a similar gene content except for two, which lacked SIX genes entirely.

New genes helped acorn barnacles adapt to a sessile lifestyle

Barnacles are the only sessile lineages among crustaceans, and their sessile life begins with the settlement of swimming larvae (cyprids) and the formation of protective shells. These processes are crucial for adaptation to a sessile lifestyle, but the underlying molecular mechanisms remain poorly understood. While investigating these mechanisms in the acorn barnacle, Amphibalanus amphitrite, we discovered a new gene, bcs-6, which is involved in the energy metabolism of cyprid settlement and originated from a transposon by acquiring the promoter and cis-regulatory element. Unlike mollusks, the barnacle shell comprises alternate layers of chitin and calcite and requires another new gene, bsf, which generates silk-like fibers that efficiently bind chitin and aggregate calcite in the aquatic environment. Our findings highlight the importance of exploring new genes in unique adaptative scenarios, and the results will provide important insights into gene origin and material development.

The Influence of a Genetic Variant in CCDC78 on LMNA-Associated Skeletal Muscle Disease

Mutations in the LMNA gene-encoding A-type lamins can cause Limb-Girdle muscular dystrophy Type 1B (LGMD1B). This disease presents with weakness and wasting of the proximal skeletal muscles and has a variable age of onset and disease severity. This variability has been attributed to genetic background differences among individuals; however, such variants have not been well characterized. To identify such variants, we investigated a multigeneration family in which affected individuals are diagnosed with LGMD1B. The primary genetic cause of LGMD1B in this family is a dominant mutation that activates a cryptic splice site, leading to a five-nucleotide deletion in the mature mRNA. This results in a frame shift and a premature stop in translation. Skeletal muscle biopsies from the family members showed dystrophic features of variable severity, with the muscle fibers of some family members possessing cores, regions of sarcomeric disruption, and a paucity of mitochondria, not commonly associated with LGMD1B. Using whole genome sequencing (WGS), we identified 21 DNA sequence variants that segregate with the family members possessing more profound dystrophic features and muscle cores. These include a relatively common variant in coiled-coil domain containing protein 78 (CCDC78). This variant was given priority because another mutation in CCDC78 causes autosomal dominant centronuclear myopathy-4, which causes cores in addition to centrally positioned nuclei. Therefore, we analyzed muscle biopsies from family members and discovered that those with both the LMNA mutation and the CCDC78 variant contain muscle cores that accumulated both CCDC78 and RyR1. Muscle cores containing mislocalized CCDC78 and RyR1 were absent in the less profoundly affected family members possessing only the LMNA mutation. Taken together, our findings suggest that a relatively common variant in CCDC78 can impart profound muscle pathology in combination with a LMNA mutation and accounts for variability in skeletal muscle disease phenotypes.

A chromosome-level genome reveals genome evolution and molecular basis of anthraquinone biosynthesis in Rheum palmatum

BACKGROUND

Rhubarb is one of common traditional Chinese medicine with a diverse array of therapeutic efficacies. Despite its widespread use, molecular research into rhubarb remains limited, constraining our comprehension of the geoherbalism.

RESULTS

We assembled the genome of Rheum palmatum L., one of the source plants of rhubarb, to elucidate its genome evolution and unpack the biosynthetic pathways of its bioactive compounds using a combination of PacBio HiFi, Oxford Nanopore, Illumina, and Hi-C scaffolding approaches. Around 2.8 Gb genome was obtained after assembly with more than 99.9% sequences anchored to 11 pseudochromosomes (scaffold N50 = 259.19 Mb). Transposable elements (TE) with a continuous expansion of long terminal repeat retrotransposons (LTRs) is predominant in genome size, contributing to the genome expansion of R. palmatum. Totally 30,480 genes were predicted to be protein-coding genes with 473 significantly expanded gene families enriched in diverse pathways associated with high-altitude adaptation for this species. Two successive rounds of whole genome duplication event (WGD) shared by Fagopyrum tataricum and R. palmatum were confirmed. We also identified 54 genes involved in anthraquinone biosynthesis and other 97 genes entangled in flavonoid biosynthesis. Notably, RpALS emerged as a compelling candidate gene for the octaketide biosynthesis after the key residual screening.

CONCLUSION

Overall, our findings offer not only an enhanced understanding of this remarkable medicinal plant but also pave the way for future innovations in its genetic breeding, molecular design, and functional genomic studies.

DNA Origami-Enabled Gene Localization of Repetitive Sequences

Repetitive sequences, which make up over 50% of human DNA, have diverse applications in disease diagnosis, forensic identification, paternity testing, and population genetic analysis due to their crucial functions for gene regulation. However, representative detection technologies such as sequencing and fluorescence imaging suffer from time-consuming protocols, high cost, and inaccuracy of the position and order of repetitive sequences. Here, we develop a precise and cost-effective strategy that combines the high resolution of atomic force microscopy with the shape customizability of DNA origami for repetitive sequence-specific gene localization. "Tri-block" DNA structures were specifically designed to connect repetitive sequences to DNA origami tags, thereby revealing precise genetic information in terms of position and sequence for high-resolution and high-precision visualization of repetitive sequences. More importantly, we achieved the results of simultaneous detection of different DNA repetitive sequences on the gene template with a resolution of ∼6.5 nm (19 nt). This strategy is characterized by high efficiency, high precision, low operational complexity, and low labor/time costs, providing a powerful complement to sequencing technologies for gene localization of repetitive sequences.

Genome-Wide Analysis of the Cytochrome P450 Monooxygenases in the Lichenized Fungi of the Class Lecanoromycetes

Lichens are unique organisms that exhibit a permanent symbiosis between fungi and algae or fungi and photosynthetic bacteria. Lichens have been found to produce biotechnologically valuable secondary metabolites. A handful of studies showed that tailoring enzymes such as cytochrome P450 monooxygenases (CYPs/P450s) play a key role in synthesizing these metabolites. Despite the critical role of P450s in the biosynthesis of secondary metabolites, the systematic analysis of P450s in lichens has yet to be reported. This study is aimed to address this research gap. A genome-wide analysis of P450s in five lichens from the fungal class Lecanoromycetes revealed the presence of 434 P450s that are grouped into 178 P450 families and 345 P450 subfamilies. The study indicated that none of the P450 families bloomed, and 15 P450 families were conserved in all five Lecanoromycetes. Lecanoromycetes have more P450s and higher P450 family diversity compared to Pezizomycetes. A total of 73 P450s were found to be part of secondary metabolite gene clusters, indicating their potential involvement in the biosynthesis of secondary metabolites. Annotation of P450s revealed that CYP682BG1 and CYP682BG2 from Cladonia grayi and Pseudevernia furfuracea (physodic acid chemotype) are involved in the synthesis of grayanic acid and physodic acid, CYP65FQ2 from Stereocaulon alpinum is involved in the synthesis of atranorin, and CYP6309A2 from Cladonia uncialis is involved in the synthesis of usnic acid. This study serves as a reference for future annotation of P450s in lichens.

Characterising the Physiological Responses of Chinook Salmon (Oncorhynchus tshawytscha) Subjected to Heat and Oxygen Stress

In New Zealand, during the hottest periods of the year, some salmon farms in the Marlborough Sounds reach water temperatures above the optimal range for Chinook salmon. High levels of mortality are recorded during these periods, emphasising the importance of understanding thermal stress in this species. In this study, the responses of Chinook salmon (Oncorhynchus tshawytscha) to chronic, long-term changes in temperature and dissolved oxygen were investigated. This is a unique investigation due to the duration of the stress events the fish were exposed to. Health and haematological parameters were analysed alongside gene expression results to determine the effects of thermal stress on Chinook salmon. Six copies of heat shock protein 90 (HSP90) were discovered and characterised: HSP90AA1.1a, HSP90AA1.2a, HSP90AA1.1b, HSP90AA1.2b, HSP90AB1a and HSP90AB1b, as well as two copies of SOD1, named SOD1a and SOD1b. The amino acid sequences contained features similar to those found in other vertebrate HSP90 and SOD1 sequences, and the phylogenetic tree and synteny analysis provided conclusive evidence of their relationship to other vertebrate HSP90 and SOD1 genes. Primers were designed for qPCR to enable the expression of all copies of HSP90 and SOD1 to be analysed. The expression studies showed that HSP90 and SOD1 were downregulated in the liver and spleen in response to longer term exposure to high temperatures and lower dissolved oxygen. HSP90 was also downregulated in the gill; however, the results for SOD1 expression in the gill were not conclusive. This study provides important insights into the physiological and genetic responses of Chinook salmon to temperature and oxygen stress, which are critical for developing sustainable fish aquaculture in an era of changing global climates.

Dicodon-based measures for modeling gene expression

MOTIVATION

Codon usage preference patterns have been associated with modulation of translation efficiency, protein folding, and mRNA decay. However, new studies support that codon pair usage has also a remarkable effect at the gene expression level. Here, we expand the concept of CAI to answer if codon pair usage patterns can be understood in terms of codon usage bias, or if they offer new information regarding coding translation efficiency.

RESULTS

Through the implementation of a weighting strategy to consider the dicodon contributions, we observe that the dicodon-based measure has greater correlations with gene expression level than CAI. Interestingly, we have noted that dicodons associated with a low value of adaptiveness are related to dicodons which mediate strong translational inhibition in yeast. We have also noticed that some codon-pairs have a smaller dicodon contribution than estimated by the product of the respective codon contributions.

AVAILABILITY AND IMPLEMENTATION

Scripts, implemented in Python, are freely available for download at https://zenodo.org/record/7738276#.ZBIDBtLMIdU.

Allele-aware chromosome-level genome assembly of the autohexaploid Diospyros kaki Thunb

Artificially improving persimmon (Diospyros kaki Thunb.), one of the most important fruit trees, remains challenging owing to the lack of reference genomes. In this study, we generated an allele-aware chromosome-level genome assembly for the autohexaploid persimmon 'Xiaoguotianshi' (Chinese-PCNA type) using PacBio CCS and Hi-C technology. The final assembly contained 4.52 Gb, with a contig N50 value of 5.28 Mb and scaffold N50 value of 44.01 Mb, of which 4.06 Gb (89.87%) of the assembly were anchored onto 90 chromosome-level pseudomolecules comprising 15 homologous groups with 6 allelic chromosomes in each. A total of 153,288 protein-coding genes were predicted, of which 98.60% were functionally annotated. Repetitive sequences accounted for 64.02% of the genome; and 110,480 rRNAs, 12,297 tRNAs, 1,483 miRNAs, and 3,510 snRNA genes were also identified. This genome assembly fills the knowledge gap in the autohexaploid persimmon genome, which is conducive in the study on the regulatory mechanisms underlying the major economically advantageous traits of persimmons and promoting breeding programs.

The first chromosome-level Fallopia multiflora genome assembly provides insights into stilbene biosynthesis

Fallopia multiflora (Thunb.) Harald, a vine belonging to the Polygonaceae family, is used in traditional medicine. The stilbenes contained in it have significant pharmacological activities in anti-oxidation and anti-aging. This study describes the assembly of the F. multiflora genome and presents its chromosome-level genome sequence containing 1.46 gigabases of data (with a contig N50 of 1.97 megabases), 1.44 gigabases of which was assigned to 11 pseudochromosomes. Comparative genomics confirmed that F. multiflora shared a whole-genome duplication event with Tartary buckwheat and then underwent different transposon evolution after separation. Combining genomics, transcriptomics, and metabolomics data to map a network of associated genes and metabolites, we identified two FmRS genes responsible for the catalysis of one molecule of p-coumaroyl-CoA and three molecules of malonyl-CoA to resveratrol in F. multiflora. These findings not only serve as the basis for revealing the stilbene biosynthetic pathway but will also contribute to the development of tools for increasing the production of bioactive stilbenes through molecular breeding in plants or metabolic engineering in microbes. Moreover, the reference genome of F. multiflora is a useful addition to the genomes of the Polygonaceae family.

Identification and analysis of proline-rich proteins and hybrid proline-rich proteins super family genes from Sorghum bicolor and their expression patterns to abiotic stress and zinc stimuli

Systematic genome-wide analysis of Sorghum bicolor revealed the identification of a total of 48 homologous genes comprising 21 proline-rich proteins (PRPs) and 27 hybrid proline-rich proteins (HyPRPs). Comprehensive scrutiny of these gene homologs was conducted for gene structure, phylogenetic investigations, chromosome mapping, and subcellular localization of proteins. Promoter analysis uncovered the regions rich with phosphorous- (BIHD), ammonium-, sulfur-responsive (SURE), and iron starvation-responsive (IRO2) along with biotic, abiotic, and development-specific cis-elements. Further, PRPs exhibit more methylation and acetylation sites in comparison with HyPRPs. miRNAs have been predicted which might play a role in cleavage and translation inhibition. Several of the SbPRP genes were stimulated in a tissue-specific manner under drought, salt, heat, and cold stresses. Additionally, exposure of plants to abscisic acid (ABA) and zinc (Zn) also triggered PRP genes in a tissue-dependent way. Among them, SbPRP17 has been found upregulated markedly in all tissues irrespective of the stress imposed. The expressions of SbHyPRPs, especially SbHyPRP2, SbHyPRP6, and SbHyPRP17 were activated under all stresses in all three tissues. On the other hand, SbHyPRP8 (root only) and SbHyPRP12 (all three tissues) were highly responsive to cold stress and ABA while SbHyPRP26 was induced by drought and Zn in the stem. Taken together, this study indicates the critical roles that SbPRPs and SbHyPRPs play during diverse abiotic stress conditions and notably the plausible roles that these genes play upon exposure to zinc, the crucial micronutrient in plants.

Evolution of a chordate-specific mechanism for myoblast fusion

Vertebrate myoblast fusion allows for multinucleated muscle fibers to compound the size and strength of mononucleated cells, but the evolution of this important process is unknown. We investigated the evolutionary origins and function of membrane-coalescing agents Myomaker and Myomixer in various groups of chordates. Here, we report that Myomaker likely arose through gene duplication in the last common ancestor of tunicates and vertebrates, while Myomixer appears to have evolved de novo in early vertebrates. Functional tests revealed a complex evolutionary history of myoblast fusion. A prevertebrate phase of muscle multinucleation driven by Myomaker was followed by the later emergence of Myomixer that enables the highly efficient fusion system of vertebrates. Evolutionary comparisons between vertebrate and nonvertebrate Myomaker revealed key structural and mechanistic insights into myoblast fusion. Thus, our findings suggest an evolutionary model of chordate fusogens and illustrate how new genes shape the emergence of novel morphogenetic traits and mechanisms.

False discovery rate: the Achilles' heel of proteogenomics

Proteogenomics refers to the integrated analysis of the genome and proteome that leverages mass-spectrometry (MS)-based proteomics data to improve genome annotations, understand gene expression control through proteoforms and find sequence variants to develop novel insights for disease classification and therapeutic strategies. However, proteogenomic studies often suffer from reduced sensitivity and specificity due to inflated database size. To control the error rates, proteogenomics depends on the target-decoy search strategy, the de-facto method for false discovery rate (FDR) estimation in proteomics. The proteogenomic databases constructed from three- or six-frame nucleotide database translation not only increase the search space and compute-time but also violate the equivalence of target and decoy databases. These searches result in poorer separation between target and decoy scores, leading to stringent FDR thresholds. Understanding these factors and applying modified strategies such as two-pass database search or peptide-class-specific FDR can result in a better interpretation of MS data without introducing additional statistical biases. Based on these considerations, a user can interpret the proteogenomics results appropriately and control false positives and negatives in a more informed manner. In this review, first, we briefly discuss the proteogenomic workflows and limitations in database construction, followed by various considerations that can influence potential novel discoveries in a proteogenomic study. We conclude with suggestions to counter these challenges for better proteogenomic data interpretation.

Pervasive translation in Mycobacterium tuberculosis

Most bacterial ORFs are identified by automated prediction algorithms. However, these algorithms often fail to identify ORFs lacking canonical features such as a length of >50 codons or the presence of an upstream Shine-Dalgarno sequence. Here, we use ribosome profiling approaches to identify actively translated ORFs in Mycobacterium tuberculosis. Most of the ORFs we identify have not been previously described, indicating that the M. tuberculosis transcriptome is pervasively translated. The newly described ORFs are predominantly short, with many encoding proteins of ≤50 amino acids. Codon usage of the newly discovered ORFs suggests that most have not been subject to purifying selection, and hence are unlikely to contribute to cell fitness. Nevertheless, we identify 90 new ORFs (median length of 52 codons) that bear the hallmarks of purifying selection. Thus, our data suggest that pervasive translation of short ORFs in Mycobacterium tuberculosis serves as a rich source for the evolution of new functional proteins.

A Chromosome-Level Genome Assembly of Yellowtail Kingfish (Seriola lalandi)

Yellowtail kingfish (Seriola lalandi) is a pelagic marine piscivore with a circumglobal distribution. It is particularly suitable for open ocean aquaculture owing to its large body size, fast swimming, rapid growth, and high economic value. A high-precision genome is of great significance for future genetic breeding research and large-scale aquaculture in the open ocean. PacBio, Illumina, and Hi-C data were combined to assemble chromosome-level reference genome with the size of 648.34 Mb (contig N50: 28.52 Mb). 175 contigs was anchored onto 24 chromosomes with lengths ranging from 12.28 to 34.59 Mb, and 99.79% of the whole genome sequence was covered. The BUSCOs of genome and gene were 94.20 and 95.70%, respectively. Gene families associated with adaptive behaviors, such as olfactory receptors and HSP70 gene families, expanded in the genome of S. lalandi. An analysis of selection pressure revealed 652 fast-evolving genes, among which mkxb, popdc2, dlx6, and ifitm5 may be related to rapid growth traits. The data generated in this study provide a valuable resource for understanding the genetic basis of S. lalandi traits.

No one tool to rule them all: prokaryotic gene prediction tool annotations are highly dependent on the organism of study

MOTIVATION

The biases in CoDing Sequence (CDS) prediction tools, which have been based on historic genomic annotations from model organisms, impact our understanding of novel genomes and metagenomes. This hinders the discovery of new genomic information as it results in predictions being biased towards existing knowledge. To date, users have lacked a systematic and replicable approach to identify the strengths and weaknesses of any CDS prediction tool and allow them to choose the right tool for their analysis.

RESULTS

We present an evaluation framework (ORForise) based on a comprehensive set of 12 primary and 60 secondary metrics that facilitate the assessment of the performance of CDS prediction tools. This makes it possible to identify which performs better for specific use-cases. We use this to assess 15 ab initio- and model-based tools representing those most widely used (historically and currently) to generate the knowledge in genomic databases. We find that the performance of any tool is dependent on the genome being analysed, and no individual tool ranked as the most accurate across all genomes or metrics analysed. Even the top-ranked tools produced conflicting gene collections, which could not be resolved by aggregation. The ORForise evaluation framework provides users with a replicable, data-led approach to make informed tool choices for novel genome annotations and for refining historical annotations.

AVAILABILITY AND IMPLEMENTATION

Code and datasets for reproduction and customisation are available at https://github.com/NickJD/ORForise.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

High-quality reference genome sequences of two coconut cultivars provide insights into evolution of monocot chromosomes and differentiation of fiber content and plant height

BACKGROUND

Coconut is an important tropical oil and fruit crop whose evolutionary position renders it a fantastic species for the investigation of the evolution of monocot chromosomes and the subsequent differentiation of ancient plants.

RESULTS

Here, we report the assembly and annotation of reference-grade genomes of Cn. tall and Cn. dwarf, whose genome sizes are 2.40 Gb and 2.39 Gb, respectively. The comparative analysis reveals that the two coconut subspecies diverge about 2-8 Mya while the conserved Arecaceae-specific whole-genome duplication (ω WGD) occurs approximately 47-53 Mya. It additionally allows us to reconstruct the ancestral karyotypes of the ten ancient monocot chromosomes and the evolutionary trajectories of the 16 modern coconut chromosomes. Fiber synthesis genes in Cn. tall, related to lignin and cellulose synthesis, are found at a higher copy number and expression level than dwarf coconuts. Integrated multi-omics analysis reveals that the difference in coconut plant height is the result of altered gibberellin metabolism, with both the GA20ox copy number and a single-nucleotide change in the promoter together leading to the difference in plant height between Cn. tall and Cn. dwarf.

CONCLUSION

We provide high-quality coconut genomes and reveal the genetic basis of trait differences between two coconuts through multi-omics analysis. We also reveal that the selection of plant height has been targeted for the same gene for millions of years, not only in natural selection of ancient plant as illustrated in coconut, but also for artificial selection in cultivated crops such as rice and maize.

Whole-Genome and Transcriptome Sequencing of Phlebopus portentosus Reveals Its Associated Ectomycorrhizal Niche and Conserved Pathways Involved in Fruiting Body Development

Phlebopus portentosus (Berk. and Broome) Boedijin, a widely consumed mushroom in China and Thailand, is the first species in the order Boletaceae to have been industrially cultivated on a large scale. However, to date, the lignocellulose degradation system and molecular basis of fruiting body development in P. portentosus have remained cryptic. In the present study, genome and transcriptome sequencing of P. portentosus was performed during the mycelium (S), primordium (P), and fruiting body (F) stages. A genome of 32.74 Mb with a 48.92% GC content across 62 scaffolds was obtained. A total of 9,464 putative genes were predicted from the genome, of which the number of genes related to plant cell wall-degrading enzymes was much lower than that of some saprophytic mushrooms with specific ectomycorrhizal niches. Principal component analysis of RNA-Seq data revealed that the gene expression profiles at all three stages were different. The low expression of plant cell wall-degrading genes also confirmed the limited ability to degrade lignocellulose. The expression profiles also revealed that some conserved and specific pathways were enriched in the different developmental stages of P. portentosus. Starch and sucrose metabolic pathways were enriched in the mycelium stage, while DNA replication, the proteasome and MAPK signaling pathways may be associated with maturation. These results provide a new perspective for understanding the key pathways and hub genes involved in P. portentosus development.

Characterization and expression pattern of ZB and PS transposons in zebrafish

Despite comprising much of the genome, transposons were once thought of as junk. However, transposons play many roles in the eukaryotic genome, such as providing new proteins as domesticated genes, expressing during germline-soma differentiation, function in DNA rearrangement in the offspring, and so on. We sought to describe the distribution and structural organization of the two autonomous transposons (ZB and PS) in the zebrafish genome and examine their expression patterns in embryos and adult tissues. The intact copy of ZB and PS was queried by BLAST on NCBI and ENSEMBL using default parameters. Of the copies with coverage and identity, more than 90 % were downloaded to do structural analysis. Spatial and temporal expression patterns were detected by qRT-PCR and Whole-mount in situ hybridization (WISH). There are 19 intact copies of ZB, encoding 341 amino acid residues with DD34E catalytic domain and flanked by 201bp TIRs, and seven intact PS copies, containing 425 amino acid residues with DD35D catalytic domain flanked by 28bp TIRs, were detected in the genome of zebrafish respectively. Analysis of genomic insertions indicated that both ZB and PS transposons are prone to be retained in the intron and intergenic regions of the zebrafish genome. The sense and antisense transcripts of ZB and PS were detected during embryonic development stages and exhibited similar expression patterns. The difference is that the sense strand transcript of ZB was explicitly expressed in midbrain-hindbrain boundary (MHB) and otic vesicle (OV), and pharyngeal arches and pharyngeal pouches (PA&PP) at 48 hpf. In adult zebrafish, the expressions of ZB and PS in muscle and brain are much higher than in other tissues. Our study results indicate that ZB and PS transposons may be involved in the embryonic development and regulation of somatic cells of certain adult tissues, such as the brain and muscle.

Identification and expression analysis of the β-defensin genes in the goat small intestine

Defensins represent a family of cysteine-rich peptides that have broad-spectrum antimicrobial activities and serve as a typical kind of effector molecule in the immunity. Ruminant species have a large number of β-defensins in the absence of α- and θ-defensins. It is well-known that the genomes of sheep and cattle harbor at least 43 and 57 β-defensin genes, respectively. However, the repertoire of the goat β-defensin gene family has not been fully elucidated. In this study, we identified a total of 50 β-defensins from the goat genome, including 48 functional genes and 2 pseudogenes. Cross-species genomic and evolutionary analyses showed that all of the β-defensins in goat chromosomes 8, 13 and 23 present one-to-one orthologous relationships to their sheep and cattle counterparts, whereas some β-defensin genes in goat chromosome 27 are goat-specific. Moreover, we observed that some duplicated genes in goat chromosome 27 may be derived from gene copy number variation, and the annotation of sheep and cattle β-defensins appears to be incomplete in the genome. Importantly, real-time PCR analysis showed that 17 β-defensins are expressed in the small intestine with abundant cBD1s expression. These findings significant increased our knowledge of ruminant β-defensin and provided useful information for genetic studies, as well as providing a foundation for future research exploring the role of defensins in the immune response.

Multiple Cypridina Luciferase Genes in the Genome of Individual Ostracods, Vargula hilgendorfii (Cypridina hilgendorfii)

The genomic structure of the Cypridina luciferase gene in Vargula hilgendorfii (formerly Cypridina hilgendorfii) was determined with three λ phage clones (λ34, λ45, and λ61). The luciferase genes in clones λ34 and λ61 consisted of 13 exons and 12 introns, and clone λ45 only contained exons 1-5. The splicing sites of the luciferase genes in λ34 and λ61 were conserved completely with the consensus sequence. The translated luciferases had 555 amino acid residues, which were over 98.6% identical to those of cDNA clones as previously reported. In contrast, each intron in clones λ34, λ45, and λ61 varied significantly in length. To explain the variation of intron length among the three V. hilgendorfii luciferase genes, genomic DNA was isolated from a single V. hilgendorfii specimen and the regions from exon 1-3 of the luciferase gene were amplified by polymerase chain reaction (PCR). PCR products with various lengths were detected and were confirmed as the luciferase gene fragments by Southern blot analysis. Furthermore, DNA sequence analysis indicated that at least seven luciferase gene groups might be present in the genome of a single specimen. Thus, multiple Cypridina luciferase genes exist in the genome of a single V. hilgendorfii specimen.

Genome sequence and evolution of Betula platyphylla

Betula L. (birch) is a pioneer hardwood tree species with ecological, economic, and evolutionary importance in the Northern Hemisphere. We sequenced the Betula platyphylla genome and assembled the sequences into 14 chromosomes. The Betula genome lacks evidence of recent whole-genome duplication and has the same paleoploidy level as Vitis vinifera and Prunus mume. Phylogenetic analysis of lignin pathway genes coupled with tissue-specific expression patterns provided clues for understanding the formation of higher ratios of syringyl to guaiacyl lignin observed in Betula species. Our transcriptome analysis of leaf tissues under a time-series cold stress experiment revealed the presence of the MEKK1-MKK2-MPK4 cascade and six additional mitogen-activated protein kinases that can be linked to a gene regulatory network involving many transcription factors and cold tolerance genes. Our genomic and transcriptome analyses provide insight into the structures, features, and evolution of the B. platyphylla genome. The chromosome-level genome and gene resources of B. platyphylla obtained in this study will facilitate the identification of important and essential genes governing important traits of trees and genetic improvement of B. platyphylla.

Comparative genomic and transcriptomic analyses of chemosensory genes in the citrus fruit fly Bactrocera (Tetradacus) minax

The citrus fruit fly Bactrocera (Tetradacus) minax is a major and devastating agricultural pest in Asian subtropical countries. Previous studies have shown that B. minax interacts with plant hosts via the efficient chemosensory system. However, the molecular components of the B. minax chemosensory system have not been well characterized. Herein, we identified a total of 25 putative odorant-binding receptors (OBPs), 4 single-copy chemosensory proteins (CSPs) and 53 candidate odorant receptors (ORs) using a newly generated whole-genome dataset for B. minax. This study significantly extended the chemosensation-related gene profiles (particularly, OBPs and ORs) in six other tephritid species. Comparative transcriptome analysis of adult B. minax and Bactrocera dorsalis showed that there were 14 highly expressed OBPs (FPKM > 100) in B. dorsalis and 7 highly expressed ones in B. minax. The expression level of CSP3 gene and CSP4 gene was higher in B. dorsalis than that in B. minax. Comparative genomic and transcriptomic analyses of chemosensory genes in the citrus fruit fly B. minax provided new insights for preventive control of this agriculture important pest and closely related species.

Molecular identification and characterisation of a novel chicken leukocyte immunoglobulin-like receptor A5

1. Leukocyte immunoglobulin-like receptor A5 (LILRA5) is a key molecule that regulates the immune system. However, the LILRA5 gene has not been characterised in avian species, including chickens. The present study aimed to identify and functionally characterise LILRA5 identified from two genetically disparate chicken lines, viz., Marek's disease (MD)-resistant (R) line 6.3 and MD-susceptible (S) line 7.2. 2. Multiple sequence alignment and phylogenetic analyses confirmed that the identity and similarity homologies of amino acids of LILRA5 in chicken lines 6.3 and 7.2 ranged between 93% and 93.7%, whereas those between chicken and mammals ranged between 20.9% and 43.7% and 21.1% to 43.9%, respectively. The newly cloned LILRA5 from chicken lines 6.3 and 7.2 revealed high conservation and a close relationship with other known mammalian LILRA5 proteins. 3. The results indicated that LILRA5 from chicken lines 6.3 and 7.2 was associated with phosphorylation of Src kinases and protein tyrosine phosphatase non-receptor type 11 (SHP2), which play a regulatory role in immune functions. Moreover, the results demonstrated that LILRA5 in these lines was associated with the activation of major histocompatibility complex (MHC) class I and β2-microglobulin and induced the expression of the transporter associated with antigen processing. In addition, LILRA5 in both chicken lines activated and induced Janus kinase (JAK)-signal transducer and the activator of transcription (STAT), nuclear factor kappa B (NF-κB), phosphoinositide-3-kinase (PI3K)/protein kinase B (AKT) and the extracellular signal-regulated kinase (ERK)1/2 signalling pathways; toll-like receptors; and Th1-, Th2-, and Th17- cytokines. 4. The data suggested that LILRA5 has innate immune receptors essential for macrophage immune response and provide novel insights into the regulation of immunity and immunopathology.

Pathogenic Adaptations Revealed by Comparative Genome Analyses of Two Colletotrichum spp., the Causal Agent of Anthracnose in Rubber Tree

Colletotrichum siamense and Colletotrichum australisinense cause Colletotrichum leaf disease that differ in their symptoms in rubber tree (Hevea brasiliensis), and pathogenicity of these two fungal species is also not identical on different cultivars of rubber tree. This divergence is often attributed to pathogen virulence factors, namely carbohydrate-active enzymes (CAZymes), secondary metabolites (SM), and small-secreted protein (SSP) effectors. The draft genome assembly and functional annotation of potential pathogenicity genes of both species obtained here provide an important and timely genomic resource for better understanding the biology and lifestyle of Colletotrichum spp. This should pave the way for designing more efficient disease control strategies in plantations of rubber tree. In this study, the genes associated with these categories were manually annotated in the genomes of C. australisinense GX1655 and C. siamense HBCG01. Comparative genomic analyses were performed to address the evolutionary relationships among these gene families in the two species. First, the size of genome assembly, number of predicted genes, and some of the functional categories differed significantly between the two congeners. Second, from the comparative genomic analyses, we identified some specific genes, certain higher abundance of gene families associated with CAZymes, CYP450, and SM in the genome of C. siamense, and Nep1-like proteins (NLP) in the genome of C. australisinense.

Identification of Novel Alleles and Structural Haplotypes of Major Histocompatibility Complex Class I and DRB Genes in Domestic Cat (Felis catus) by a Newly Developed NGS-Based Genotyping Method

The major histocompatibility complex (MHC) is a highly polymorphic and duplicated genomic region that encodes transplantation and immune regulatory molecules. Although it is well-known that particular MHC allelic polymorphisms and haplotypes are genetically relate to immune-mediated diseases detailed information of the cat MHC (Feline Leukocyte Antigen; FLA) genetic and haplotypic structure and diversity is limited in comparison to humans and many other species. In this study, to better understand the degree and types of allele and allelic haplotype diversity of FLA-class I (FLA-I) and FLA-DRB loci in domestic cats, we identified six expressible FLA-I loci in peripheral white blood cells by in silico estimation of the coding exons and NGS-based amplicon sequencing using five unrelated cats. We then used a newly developed NGS-based genotyping method to genotype and annotate 32 FLA-I and 16 FLA-DRB sequences in two families of 20 domestic cats. A total of 14 FLA-I and seven FLA-DRB were identified as novel polymorphic sequences. Phylogenetic analyses grouped the sequences into six FLA-I (FLA-E/H/K, FLA-A, FLA-J, FLA-L, FLA-O and a tentatively named FLA-E/H/K_Rec) and four FLA-DRB (FLA-DRB1, FLA-DRB3, FLA-DRB4, and FLA-DRB5) lineages. Pedigree analysis of two cat families revealed eight distinct FLA structural haplotypes (Class I - DRB) with five to eight FLA-I and two to three FLA-DRB transcribed loci per haplotype. It is evident that the eight FLA haplotypes were generated by gene duplications and deletions, and rearrangements by genetic recombination with the accumulation and/or inheritance of novel polymorphisms. These findings are useful for further genetic diversity analysis and disease association studies among cat breeds and in veterinary medicine.

Genomic insights of body plan transitions from bilateral to pentameral symmetry in Echinoderms

Echinoderms are an exceptional group of bilaterians that develop pentameral adult symmetry from a bilaterally symmetric larva. However, the genetic basis in evolution and development of this unique transformation remains to be clarified. Here we report newly sequenced genomes, developmental transcriptomes, and proteomes of diverse echinoderms including the green sea urchin (L. variegatus), a sea cucumber (A. japonicus), and with particular emphasis on a sister group of the earliest-diverged echinoderms, the feather star (A. japonica). We learned that the last common ancestor of echinoderms retained a well-organized Hox cluster reminiscent of the hemichordate, and had gene sets involved in endoskeleton development. Further, unlike in other animal groups, the most conserved developmental stages were not at the body plan establishing phase, and genes normally involved in bilaterality appear to function in pentameric axis development. These results enhance our understanding of the divergence of protostomes and deuterostomes almost 500 Mya.

Li et al. investigate the evolution and genetic basis of the adult pentameral body plan in echinoderms using genomic, transcriptomic, and proteomic data. They determine that the last common ancestor of echinoderms contained an organized Hox cluster and endoskeleton genes, and suggest that cooption of bilateral development genes was involved in evolution of the pentameric body plan.

Echinoderms provide missing link in the evolution of PrRP/sNPF-type neuropeptide signalling

Neuropeptide signalling systems comprising peptide ligands and cognate receptors are evolutionarily ancient regulators of physiology and behaviour. However, there are challenges associated with determination of orthology between neuropeptides in different taxa. Orthologs of vertebrate neuropeptide-Y (NPY) known as neuropeptide-F (NPF) have been identified in protostome invertebrates, whilst prolactin-releasing peptide (PrRP) and short neuropeptide-F (sNPF) have been identified as paralogs of NPY/NPF in vertebrates and protostomes, respectively. Here we investigated the occurrence of NPY/NPF/PrRP/sNPF-related signalling systems in a deuterostome invertebrate phylum - the Echinodermata. Analysis of transcriptome/genome sequence data revealed loss of NPY/NPF-type signalling, but orthologs of PrRP-type neuropeptides and sNPF/PrRP-type receptors were identified in echinoderms. Furthermore, experimental studies revealed that the PrRP-type neuropeptide pQDRSKAMQAERTGQLRRLNPRF-NH₂ is a potent ligand for a sNPF/PrRP-type receptor in the starfish Asterias rubens. Our findings indicate that PrRP-type and sNPF-type signalling systems are orthologous and originated as a paralog of NPY/NPF-type signalling in Urbilateria.

Sex differences in testosterone reactivity and sensitivity in a non-model gerbil

Although testosterone (T) is a key regulator in vertebrate development, physiology, and behaviour in both sexes, studies suggest that its regulation may be sex-specific. We measured circulating T levels in Baluchistan gerbils (Gerbillus nanus) in the field and in the lab all year round and found no significant sex differences. However, we observed sex differences in circulating T levels following gonadotropin-releasing hormone (GnRH) challenge and T implants in this non-model species. Whereas only males elevated T following a GnRH challenge, females had higher serum T concentrations following T implant insertion. These differences may be a result of different points of regulation along the hypothalamic-pituitary-gonadal (HPG) axis. Consequently, we examined sex differences in the mRNA expression of the androgen receptor (AR) in multiple brain regions. We identified AR and β-actin sequences in assembled genomic sequences of members of the Gerbillinae, which were analogous to rat sequences, and designed primers for them. The distribution of the AR in G. nanus brain regions was similar to documented expression profiles in rodents. We found lower AR mRNA levels in females in the striatum. Additionally, G. nanus that experienced housing in mixed-sex pairs had higher adrenal AR expression than G. nanus that were housed alone. Regulation of the gerbil HPG axis may reflect evolutionary sex differences in life-history strategies, with males ready to reproduce when receptive females are available, while the possible reproductive costs associated with female T direct its regulation upstream.

A Markov chain-based feature extraction method for classification and identification of cancerous DNA sequences

Introduction: In recent decades, the growing rate of cancer incidence is a big concern for most societies. Due to the genetic origins of cancer disease, its internal structure is necessary for the study of this disease.

Methods: In this research, cancer data are analyzed based on DNA sequences. The transition probability of occurring two pairs of nucleotides in DNA sequences has Markovian property. This property inspires the idea of feature dimension reduction of DNA sequence for overcoming the high computational overhead of genes analysis. This idea is utilized in this research based on the Markovian property of DNA sequences. This mapping decreases feature dimensions and conserves basic properties for discrimination of cancerous and non-cancerous genes.

Results: The results showed that a non-linear support vector machine (SVM) classifier with RBF and polynomial kernel functions can discriminate selected cancerous samples from non-cancerous ones. Experimental results based on the 10-fold cross-validation and accuracy metrics verified that the proposed method has low computational overhead and high accuracy.

Conclusion: The proposed algorithm was successfully tested on related research case studies. In general, a combination of proposed Markovian-based feature reduction and non-linear SVM classifier can be considered as one of the best methods for discrimination of cancerous and non-cancerous genes.

A Novel pH-Regulated, Unusual 603 bp Overlapping Protein Coding Gene pop Is Encoded Antisense to ompA in Escherichia coli O157:H7 (EHEC)

Antisense transcription is well known in bacteria. However, translation of antisense RNAs is typically not considered, as the implied overlapping coding at a DNA locus is assumed to be highly improbable. Therefore, such overlapping genes are systematically excluded in prokaryotic genome annotation. Here we report an exceptional 603 bp long open reading frame completely embedded in antisense to the gene of the outer membrane protein ompA. An active σ70 promoter, transcription start site (TSS), Shine-Dalgarno motif and rho-independent terminator were experimentally validated, providing evidence that this open reading frame has all the structural features of a functional gene. Furthermore, ribosomal profiling revealed translation of the mRNA, the protein was detected in Western blots and a pH-dependent phenotype conferred by the protein was shown in competitive overexpression growth experiments of a translationally arrested mutant versus wild type. We designate this novel gene pop (pH-regulated overlapping protein-coding gene), thus adding another example to the growing list of overlapping, protein coding genes in bacteria.

The genome of jojoba (Simmondsia chinensis): A taxonomically isolated species that directs wax ester accumulation in its seeds

Seeds of the desert shrub, jojoba (Simmondsia chinensis), are an abundant, renewable source of liquid wax esters, which are valued additives in cosmetic products and industrial lubricants. Jojoba is relegated to its own taxonomic family, and there is little genetic information available to elucidate its phylogeny. Here, we report the high-quality, 887-Mb genome of jojoba assembled into 26 chromosomes with 23,490 protein-coding genes. The jojoba genome has only the whole-genome triplication (γ) shared among eudicots and no recent duplications. These genomic resources coupled with extensive transcriptome, proteome, and lipidome data helped to define heterogeneous pathways and machinery for lipid synthesis and storage, provided missing evolutionary history information for this taxonomically segregated dioecious plant species, and will support efforts to improve the agronomic properties of jojoba.

Ambiguous splice sites distinguish circRNA and linear splicing in the human genome

MOTIVATION

Identification of splice sites is critical to gene annotation and to determine which sequences control circRNA biogenesis. Full-length RNA transcripts could in principle complete annotations of introns and exons in genomes without external ontologies, i.e., ab initio. However, whether it is possible to reconstruct genomic positions where splicing occurs from full-length transcripts, even if sampled in the absence of noise, depends on the genome sequence composition. If it is not, there exist provable limits on the use of RNA-Seq to define splice locations (linear or circular) in the genome.

RESULTS

We provide a formal definition of splice site ambiguity due to the genomic sequence by introducing equivalent junction, which is the set of local genomic positions resulting in the same RNA sequence when joined through RNA splicing. We show that equivalent junctions are prevalent in diverse eukaryotic genomes and occur in 88.64% and 78.64% of annotated human splice sites in linear and circRNA junctions, respectively. The observed fractions of equivalent junctions and the frequency of many individual motifs are statistically significant when compared against the null distribution computed via simulation or closed-form. The frequency of equivalent junctions establishes a fundamental limit on the possibility of ab initio reconstruction of RNA transcripts without appealing to the ontology of "GT-AG" boundaries defining introns. Said differently, completely ab initio is impossible in the vast majority of splice sites in annotated circRNAs and linear transcripts.

AVAILABILITY AND IMPLEMENTATION

Two python scripts generating an equivalent junction sequence per junction are available at: https://github.com/salzmanlab/Equivalent-Junctions.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Genome assembly provides insights into the genome evolution and flowering regulation of orchardgrass

Orchardgrass (Dactylis glomerata L.) is an important forage grass for cultivating livestock worldwide. Here, we report an ~1.84-Gb chromosome-scale diploid genome assembly of orchardgrass, with a contig N50 of 0.93 Mb, a scaffold N50 of 6.08 Mb and a super-scaffold N50 of 252.52 Mb, which is the first chromosome-scale assembled genome of a cool-season forage grass. The genome includes 40 088 protein-coding genes, and 69% of the assembled sequences are transposable elements, with long terminal repeats (LTRs) being the most abundant. The LTRretrotransposons may have been activated and expanded in the grass genome in response to environmental changes during the Pleistocene between 0 and 1 million years ago. Phylogenetic analysis reveals that orchardgrass diverged after rice but before three Triticeae species, and evolutionarily conserved chromosomes were detected by analysing ancient chromosome rearrangements in these grass species. We also resequenced the whole genome of 76 orchardgrass accessions and found that germplasm from Northern Europe and East Asia clustered together, likely due to the exchange of plants along the 'Silk Road' or other ancient trade routes connecting the East and West. Last, a combined transcriptome, quantitative genetic and bulk segregant analysis provided insights into the genetic network regulating flowering time in orchardgrass and revealed four main candidate genes controlling this trait. This chromosome-scale genome and the online database of orchardgrass developed here will facilitate the discovery of genes controlling agronomically important traits, stimulate genetic improvement of and functional genetic research on orchardgrass and provide comparative genetic resources for other forage grasses.

Isolation of a rhamnogalacturonan lyase expressed during ripening of the Chilean strawberry fruit and its biochemical characterization

Fragaria chiloensis (L.) Mill. fruit has exotic organoleptic properties however commercialization is a challenge due to its fast and intensive softening. Texture modifications associated to ripening are related to cell wall metabolism. Main cell wall polysaccharides metabolized in F. chiloensis fruit are pectins, being rhamnogalacturonan I (RG-I) an abundant pectin domain in strawberry. Several enzymes belonging to the fruit molecular machinery have been described to act on different cell wall polysaccharides in F. chiloensis, but none acting on the main chain of RG-I until now. A gene sequence coding for a rhamnogalacturonan endolyase (RG-lyase) (EC 4.2.2.23) was isolated from F. chiloensis. The FchRGL1 sequence belongs to Polysaccharide Lyase family 4 and contains the three functional domains of RG-lyases: RGL4 domain, fibronectin type III and the carbohydrate binding module. In addition, it contains key amino acid residues for activity and Ca²⁺ coordination. qRT-PCR analyses indicate that FchRGL1 transcripts increase in fruit throughout ripening. RG-lyase activity evidences a remarkable increase as the fruit ripens. The heterologous expression of FchRGL1 in Pichia pastoris provided an active protein that allows its biochemical characterization. RG-lyase activity is optimum at pH 5.0, 25-30 °C and 2 mM Ca²⁺. A K_M of 0.086 mg mL^-1 was determined for potato RG-I, and the enzyme undergoes inhibition at high substrate concentration. The enzyme is also able to degrade the mucilage of germinating A. thaliana's seeds. Finally, the properties of FchRGL1 and its expression pattern are congruent with a crucial role in cell wall re-organization during softening of F. chiloensis fruit.

Igh locus structure and evolution in Platyrrhines: new insights from a genomic perspective

Non-human primates have been used as animal models because of their phylogenetic closeness to humans. However, the genetic differences between humans and non-human primates must be considered to select the appropriate animal models. Recently, New World monkeys (Platyrrhines) have generated a higher interest in biomedical research, especially in assessing vaccine safety and immunogenicity. Given the continued and renewed interest in Platyrrhines as biomedical models, it is a necessary to have a better and more complete understanding of their immune system and its implications for research. Immunoglobulins (Ig) are the main proteins that mediate humoral immunity. These proteins have evolved as part of an adaptive immune response system derived from ancient vertebrates. There are at least four Ig classes in Prosimians, whereas five have been reported in Catarrhines. Information on the structure and evolution of the loci containing immunoglobulin heavy chain constant genes (Igh) in Platyrrhines, however, is limited. Here, Igh loci were characterized in 10 Platyrrhines using the available whole genome sequences. Human and Macaca Igh loci were also assessed to compare them with their Platyrrhines counterparts. Differences in Igh locus structure were observed between Platyrrhines and Catarrhines. Noteworthy changes occur in the γ gene, which encodes a key Ig involved in organism defense that would favor protection after vaccination. The remarkable differences between the immunoglobulin proteins of Platyrrhines and Catarrhines warrant a cautionary message to biomedical researchers.

GPRED-GC: a Gene PREDiction model accounting for 5 ′- 3′ GC gradient

Background

Gene is a key step in genome annotation. Ab initio gene prediction enables gene annotation of new genomes regardless of availability of homologous sequences. There exist a number of ab initio gene prediction tools and they have been widely used for gene annotation for various species. However, existing tools are not optimized for identifying genes with highly variable GC content. In addition, some genes in grass genomes exhibit a sharp 5 ^′- 3^′ decreasing GC content gradient, which is not carefully modeled by available gene prediction tools. Thus, there is still room to improve the sensitivity and accuracy for predicting genes with GC gradients.

Results

In this work, we designed and implemented a new hidden Markov model (HMM)-based ab initio gene prediction tool, which is optimized for finding genes with highly variable GC contents, such as the genes with negative GC gradients in grass genomes. We tested the tool on three datasets from Arabidopsis thaliana and Oryza sativa. The results showed that our tool can identify genes missed by existing tools due to the highly variable GC contents.

Conclusions

GPRED-GC can effectively predict genes with highly variable GC contents without manual intervention. It provides a useful complementary tool to existing ones such as Augustus for more sensitive gene discovery. The source code is freely available at https://sourceforge.net/projects/gpred-gc/.

SmPPT, a 4-hydroxybenzoate polyprenyl diphosphate transferase gene involved in ubiquinone biosynthesis, confers salt tolerance in Salvia miltiorrhiza

SmPPT, which encodes 4-hydroxybenzoate polyprenyl diphosphate transferase involved in ubiquinone biosynthesis, confers salt tolerance to S. miltiorrhiza through enhancing the activities of POD and CAT to scavenge ROS. Ubiquinone (UQ), also known as coenzyme Q (CoQ), is a key electron transporter in the mitochondrial respiratory system. UQ is composed of a benzene quinone ring and a polyisoprenoid side chain. Attachment of polyisoprenoid side chain to the benzene quinone ring is a rate-limiting step catalyzed by 4-hydroxybenzoate polyprenyl diphosphate transferase (PPT). So far, only a few plant PPT-encoding genes have been functionally analyzed. Through genome-wide analysis and subsequent molecular cloning, a PPT-encoding gene, termed SmPPT, was identified from an economically and academically important medicinal model plant, Salvia miltiorrhiza. SmPPT contained many putative cis-elements associated with abiotic stresses in the promoter region and were responsive to PEG-6000 and methyl jasmonate treatments. The deduced SmPPT protein contains the PT_UbiA conserved domain of polyprenyl diphosphate transferase and an N-terminal mitochondria transit peptide. Transient expression assay of SmPPT-GFP fusion protein showed that SmPPT was mainly localized in the mitochondria. SmPPT could functionally complement coq2 mutation and catalyzed UQ6 production in yeast cells. Overexpression of SmPPT increased UQ production and enhanced salt tolerance in S. miltiorrhiza. Under salinity stress conditions, transgenic plants accumulated less H₂O₂ and malondialdehyde and exhibited higher peroxidase (POD) and catalase (CAT) activities compared with wild-type plants. It indicates that SmPPT confers salt tolerance to S. miltiorrhiza at least partially through enhancing the activities of POD and CAT to scavenge ROS.

Towards the Complete Goat Pan-Genome by Recovering Missing Genomic Segments From the Reference Genome

It is broadly expected that next generation sequencing will ultimately generate a complete genome as is the latest goat reference genome (ARS1), which is considered to be one of the most continuous assemblies in livestock. However, the rich diversity of worldwide goat breeds indicates that a genome from one individual would be insufficient to represent the whole genomic contents of goats. By comparing nine de novo assemblies from seven sibling species of domestic goat with ARS1 and using resequencing and transcriptome data from goats for verification, we identified a total of 38.3 Mb sequences that were absent in ARS1. The pan-sequences contain genic fractions with considerable expression. Using the pan-genome (ARS1 together with the pan-sequences) as a reference genome, variation calling efficacy can be appreciably improved. A total of 56,657 spurious SNPs per individual were repressed and 24,414 novel SNPs per individual on average were recovered as a result of better reads mapping quality. The transcriptomic mapping rate was also increased by ∼1.15%. Our study demonstrated that comparing de novo assemblies from closely related species is an efficient and reliable strategy for finding missing sequences from the reference genome and could be applicable to other species. Pan-genome can serve as an improved reference genome in animals for a better exploration of the underlying genomic variations and could increase the probability of finding genotype-phenotype associations assessed by a comprehensive variation database containing much more differences between individuals. We have constructed a goat pan-genome web interface for data visualization (http://animal.nwsuaf.edu.cn/panGoat).

Repertoire-wide gene structure analyses: a case study comparing automatically predicted and manually annotated gene models

Background

The location and modular structure of eukaryotic protein-coding genes in genomic sequences can be automatically predicted by gene annotation algorithms. These predictions are often used for comparative studies on gene structure, gene repertoires, and genome evolution. However, automatic annotation algorithms do not yet correctly identify all genes within a genome, and manual annotation is often necessary to obtain accurate gene models and gene sets. As manual annotation is time-consuming, only a fraction of the gene models in a genome is typically manually annotated, and this fraction often differs between species. To assess the impact of manual annotation efforts on genome-wide analyses of gene structural properties, we compared the structural properties of protein-coding genes in seven diverse insect species sequenced by the i5k initiative.

Results

Our results show that the subset of genes chosen for manual annotation by a research community (3.5–7% of gene models) may have structural properties (e.g., lengths and exon counts) that are not necessarily representative for a species’ gene set as a whole. Nonetheless, the structural properties of automatically generated gene models are only altered marginally (if at all) through manual annotation. Major correlative trends, for example a negative correlation between genome size and exonic proportion, can be inferred from either the automatically predicted or manually annotated gene models alike. Vice versa, some previously reported trends did not appear in either the automatic or manually annotated gene sets, pointing towards insect-specific gene structural peculiarities.

Conclusions

In our analysis of gene structural properties, automatically predicted gene models proved to be sufficiently reliable to recover the same gene-repertoire-wide correlative trends that we found when focusing on manually annotated gene models only. We acknowledge that analyses on the individual gene level clearly benefit from manual curation. However, as genome sequencing and annotation projects often differ in the extent of their manual annotation and curation efforts, our results indicate that comparative studies analyzing gene structural properties in these genomes can nonetheless be justifiable and informative.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-6064-8) contains supplementary material, which is available to authorized users.

Genomic Identification and Expression Analysis of the Cathelicidin Gene Family of the Forest Musk Deer

Simple Summary

Cathelicidins are a group of host defense peptides in vertebrates with both antimicrobial and immunomodulatory activities. In the present study, we identified the entire repertoire of the cathelicidin gene family from the forest musk deer genome. Sequence comparison, phylogenetic topology, and gene and genomic organizations collectively suggest that all cathelicidin genes have already been fixed in the genome of forest musk deer before the split of moschidae and bovidae, while independent pseudogenization events have occurred after species divergence. In addition, real-time PCR analysis suggested that all functional cathelicidins play important roles in the immune system. The results of this study will be helpful for further evolutionary and functional studies.

Abstract

The forest musk deer (Moschus berezovskii) is a small-sized artiodactyl species famous for the musk secreted by adult males. In the captive population, this species is under the threat of infection diseases, which greatly limits the increase of individual numbers. In the present study, we computationally analyzed the repertoire of the cathelicidin (CATHL) family from the genome of forest musk deer and investigated their expression pattern by real-time PCR. Our results showed that the entire genome of forest musk deer encodes eight cathelicidins, including six functional genes and two pseudogenes. Phylogenetic analyses further revealed that all forest musk deer cathelicidin members have emerged before the split of the forest musk deer and cattle and that forest musk deer CATHL3L2 and CATHL9 are orthologous with two cattle pseudogenes. In addition, the gene expression results showed that the six functional genes are not only abundantly expressed in the spleen and lung, but are also differently expressed in response to abscesses, which suggests that forest musk deer cathelicidins may be involved in infections. Taken together, identification and characterization of the forest musk deer cathelicidins provide fundamental data for further investigating their evolutionary process and biological functions.

Functional Analysis of Polyprenyl Diphosphate Synthase Genes Involved in Plastoquinone and Ubiquinone Biosynthesis in Salvia miltiorrhiza

Polyprenyl diphosphate synthase (PPS) plays important roles in the biosynthesis of functionally important plastoquinone (PQ) and ubiquinone (UQ). However, only few plant PPS genes have been functionally characterized. Through genome-wide analysis, two PPS genes, termed SmPPS1 and SmPPS2, were identified from Salvia miltiorrhiza, an economically significant Traditional Chinese Medicine material and an emerging model medicinal plant. SmPPS1 and SmPPS2 belonged to different phylogenetic subgroups of plant trans-long-chain prenyltransferases and exhibited differential tissue expression and light-induced expression patterns. Computational prediction and transient expression assays showed that SmPPS1 was localized in the chloroplasts, whereas SmPPS2 was mainly localized in the mitochondria. SmPPS2, but not SmPPS1, could functionally complement the coq1 mutation in yeast cells and catalyzed the production of UQ-9 and UQ-10. Consistently, both UQ-9 and UQ-10 were detected in S. miltiorrhiza plants. Overexpression of SmPPS2 caused significant UQ accumulation in S. miltiorrhiza transgenics, whereas down-regulation resulted in decreased UQ content. Differently, SmPPS1 overexpression significantly elevated PQ-9 content in S. miltiorrhiza. Transgenic lines showing a down-regulation of SmPPS1 expression exhibited decreased PQ-9 level, abnormal chloroplast and trichome development, and varied leaf bleaching phenotypes. These results suggest that SmPPS1 is involved in PQ-9 biosynthesis, whereas SmPPS2 is involved in UQ-9 and UQ-10 biosynthesis.