Genome-wide analyses supported by RNA-Seq reveal non-canonical splice sites in plant genomes

Difference Analysis Among Six Kinds of Acceptor Splicing Sequences by the Dispersion Features of 6-mer Subsets in Human Genes

Identifying the sequence composition of different splicing modes is a challenge in current research. This study explored the dispersion distributions of 6-mer subsets in human acceptor splicing regions. Without differentiating acceptor splicing modes, obvious differences were observed across the upstream, core, and downstream regions of splicing sites for 16 dispersion distributions. These findings indicate that the dispersion value of each subset can effectively characterize the compositional properties of splicing sequences. When acceptor splicing sequences were classified into common, constitutive, and alternative modes, the differences in dispersion distributions for most of the XY1 6-mer subsets were significant among the three splicing modes. Furthermore, the alternative splicing mode was classified into normal, exonic, and intronic sub-modes, the differences in dispersion distributions for most of the XY1 6-mer subsets were also significant among the three splicing sub-modes. Our results indicate that dispersion values of XY1 6-mer subsets not only revealed the sequence composition patterns of acceptor splicing regions but also effectively identified the differences in base correlation among various acceptor splicing modes. Our research provides new insights into revealing and predicting different splicing modes.

A novel human protein-coding locus identified using a targeted RNA enrichment technique

BACKGROUND

Accurate and comprehensive genomic annotation, including the full list of protein-coding genes, is vital for understanding the molecular mechanisms of human biology. We have previously shown that the genome contains a multitude of yet hidden functional exons and transcripts, some of which might represent novel mRNAs. These results resonate with those from other groups and strongly argue that two decades after the completion of the first draft of the human genome sequence, the current annotation of human genes and transcripts remains far from being complete.

RESULTS

Using a targeted RNA enrichment technique, we showed that one of the novel functional exons previously discovered by us and currently annotated as part of a long non-coding RNA, is actually a part of a novel protein-coding gene, InSETG-4, which encodes a novel human protein with no known homologs or motifs. We found that InSETG-4 is induced by various DNA-damaging agents across multiple cell types and therefore might represent a novel component of DNA damage response. Despite its low abundance in bulk cell populations, InSETG-4 exhibited expression restricted to a small fraction of cells, as demonstrated by the amplification-based single-molecule fluorescence in situ hybridization (asmFISH) analysis.

CONCLUSIONS

This study argues that yet undiscovered human protein-coding genes exist and provides an example of how targeted RNA enrichment techniques can help to fill this major gap in our knowledge of the information encoded in the human genome.

The golden genome annotation of Ganoderma lingzhi reveals a more complex scenario of eukaryotic gene structure and transcription activity

BACKGROUND

It is generally accepted that nuclear genes in eukaryotes are located independently on chromosomes and expressed in a monocistronic manner. However, accumulating evidence suggests a more complex landscape of gene structure and transcription. Ganoderma lingzhi, a model medicinal fungus, currently lacks high-quality genome annotation, hindering genetic studies.

RESULTS

Here, we reported a golden annotation of G. lingzhi, featuring 14,147 high-confidence genes derived from extensive manual corrections. Novel characteristics of gene structure and transcription were identified accordingly. Notably, non-canonical splicing sites accounted for 1.99% of the whole genome, with the predominant types being GC-AG (1.85%), GT-AC (0.05%), and GT-GG (0.04%). 1165 pairs of genes were found to have overlapped transcribed regions, and 92.19% of which showed opposite directions of gene transcription. A total of 5,412,158 genetic variations were identified among 13 G. lingzhi strains, and the manually corrected gene sets resulted in enhanced functional annotation of these variations. More than 60% of G. lingzhi genes were alternatively spliced. In addition, we found that two or more protein-coding genes (PCGs) can be transcribed into a single RNA molecule, referred to as polycistronic genes. In total, 1272 polycistronic genes associated with 2815 PCGs were identified.

CONCLUSIONS

The widespread presence of polycistronic genes in G. lingzhi strongly complements the theory that polycistron is also present in eukaryotic genomes. The extraordinary gene structure and transcriptional activity uncovered through this golden annotation provide implications for the study of genes, genomes, and related studies in G. lingzhi and other eukaryotes.

Tracing the evolutionary and genetic footprints of atmospheric tillandsioids transition from land to air

Plant evolution is driven by key innovations of functional traits that enables their survivals in diverse ecological environments. However, plant adaptive evolution from land to atmospheric niches remains poorly understood. In this study, we use the epiphytic Tillandsioideae subfamily of Bromeliaceae as model plants to explore their origin, evolution and diversification. We provide a comprehensive phylogenetic tree based on nuclear transcriptomic sequences, indicating that core tillandsioids originated approximately 11.3 million years ago in the Andes. The geological uplift of the Andes drives the divergence of tillandsioids into tank-forming and atmospheric types. Our genomic and transcriptomic analyses reveal gene variations and losses associated with adaptive traits such as impounding tanks and absorptive trichomes. Furthermore, we uncover specific nitrogen-fixing bacterial communities in the phyllosphere of tillandsioids as potential source of nitrogen acquisition. Collectively, our study provides integrative multi-omics insights into the adaptive evolution of tillandsioids in response to elevated aerial habitats.

DRANetSplicer: A Splice Site Prediction Model Based on Deep Residual Attention Networks

The precise identification of splice sites is essential for unraveling the structure and function of genes, constituting a pivotal step in the gene annotation process. In this study, we developed a novel deep learning model, DRANetSplicer, that integrates residual learning and attention mechanisms for enhanced accuracy in capturing the intricate features of splice sites. We constructed multiple datasets using the most recent versions of genomic data from three different organisms, Oryza sativa japonica, Arabidopsis thaliana and Homo sapiens. This approach allows us to train models with a richer set of high-quality data. DRANetSplicer outperformed benchmark methods on donor and acceptor splice site datasets, achieving an average accuracy of (96.57%, 95.82%) across the three organisms. Comparative analyses with benchmark methods, including SpliceFinder, Splice2Deep, Deep Splicer, EnsembleSplice, and DNABERT, revealed DRANetSplicer's superior predictive performance, resulting in at least a (4.2%, 11.6%) relative reduction in average error rate. We utilized the DRANetSplicer model trained on O. sativa japonica data to predict splice sites in A. thaliana, achieving accuracies for donor and acceptor sites of (94.89%, 94.25%). These results indicate that DRANetSplicer possesses excellent cross-organism predictive capabilities, with its performance in cross-organism predictions even surpassing that of benchmark methods in non-cross-organism predictions. Cross-organism validation showcased DRANetSplicer's excellence in predicting splice sites across similar organisms, supporting its applicability in gene annotation for understudied organisms. We employed multiple methods to visualize the decision-making process of the model. The visualization results indicate that DRANetSplicer can learn and interpret well-known biological features, further validating its overall performance. Our study systematically examined and confirmed the predictive ability of DRANetSplicer from various levels and perspectives, indicating that its practical application in gene annotation is justified.

Where the minor things are: a pan-eukaryotic survey suggests neutral processes may explain much of minor intron evolution

Spliceosomal introns are gene segments removed from RNA transcripts by ribonucleoprotein machineries called spliceosomes. In some eukaryotes a second 'minor' spliceosome is responsible for processing a tiny minority of introns. Despite its seemingly modest role, minor splicing has persisted for roughly 1.5 billion years of eukaryotic evolution. Identifying minor introns in over 3000 eukaryotic genomes, we report diverse evolutionary histories including surprisingly high numbers in some fungi and green algae, repeated loss, as well as general biases in their positional and genic distributions. We estimate that ancestral minor intron densities were comparable to those of vertebrates, suggesting a trend of long-term stasis. Finally, three findings suggest a major role for neutral processes in minor intron evolution. First, highly similar patterns of minor and major intron evolution contrast with both functionalist and deleterious model predictions. Second, observed functional biases among minor intron-containing genes are largely explained by these genes' greater ages. Third, no association of intron splicing with cell proliferation in a minor intron-rich fungus suggests that regulatory roles are lineage-specific and thus cannot offer a general explanation for minor splicing's persistence. These data constitute the most comprehensive view of minor introns and their evolutionary history to date, and provide a foundation for future studies of these remarkable genetic elements.

Spliceosomal introns in the diplomonad parasite Giardia duodenalis revisited

Complete reference genomes, including correct feature annotations, are a fundamental aspect of genomic biology. In the case of protozoan species such as Giardia duodenalis, a major human and animal parasite worldwide, accurate genome annotation can deepen our understanding of the evolution of parasitism and pathogenicity by identifying genes underlying key traits and clinically relevant cellular mechanisms, and by extension, the development of improved prevention strategies and treatments. This study used bioinformatics analyses of Giardia mRNA libraries to characterize known introns and identify new intron candidates, working towards completion of the G. duodenalis assemblage A strain 'WB' genome and further elucidating Giardia's gene expression. By using a set of experimentally validated positive control loci to calibrate our intron detection pipeline, we were able to detect evidence of previously missed candidate splice junctions directly from expressed transcript data. These intron candidates were further studied in silico using NMDS (non-metric multidimensional scaling) clustering to determine shared characteristics and their relative importance such as secondary structure, splicing efficiency and motif conservation, and thus to refine intron models. Results from this study identified 34 new intron candidates, with several potential introns showing evidence that secondary structure of the mRNA molecule might play a more significant role in splicing than previously reported eukaryotic splicing activity mediated by a reduced spliceosome present in G. duodenalis.

Highly effective mlo-based powdery mildew resistance in hexaploid wheat without pleiotropic effects

Application of the mlo-based resistance in barley against powdery mildew attacks is a major success in crop breeding, since it confers durable disease resistance. Resistance caused by mutations in the Mlo gene seems to be ubiquitous across a range of species. This work addresses the introduction of mlo-based resistance into hexaploid wheat, which is complicated by the occurrence of three homoeologous genes: Mlo-A1, Mlo-B1 and Mlo-D1. EMS-generated mutant plants were screened for mutations in the three homoeologues. We selected and combined 6, 8, and 4 mutations, respectively, to obtain triple homozygous mlo mutant lines. Twenty-four mutant lines showed highly effective resistance towards attack by the powdery mildew pathogen under field conditions. All 18 mutations appeared to contribute to resistance; however, they had different effects on the occurrence of symptoms such as chlorotic and necrotic spots, which are pleiotropic to the mlo-based powdery mildew resistance. We conclude that to obtain highly effective powdery mildew resistance in wheat and to avoid detrimental pleiotropic effects, all three Mlo homoeologues should be mutated; however, at least one of the mutations should be of the weaker type in order to alleviate strong pleiotropic effects from the other mutations.

Sequence variations affect the 5' splice site selection of plant introns

Introns are noncoding sequences spliced out of pre-mRNAs by the spliceosome to produce mature mRNAs. The 5' ends of introns mostly begin with GU and have a conserved sequence motif of AG/GUAAGU that could base-pair with the core sequence of U1 snRNA of the spliceosome. Intriguingly, ∼ 1% of introns in various eukaryotic species begin with GC. This occurrence could cause misannotation of genes; however, the underlying splicing mechanism is unclear. We analyzed the sequences around the intron 5' splice site (ss) in Arabidopsis (Arabidopsis thaliana) and found sequences at the GC intron ss are much more stringent than those of GT introns. Mutational analysis at various positions of the intron 5' ss revealed that although mutations impair base pairing, different mutations at the same site can have different effects, suggesting that steric hindrance also affects splicing. Moreover, mutations of 5' ss often activate a hidden ss nearby. Our data suggest that the 5' ss is selected via a competition between the major ss and the nearby minor ss. This work not only provides insights into the splicing mechanism of intron 5' ss but also improves the accuracy of gene annotation and the study of the evolution of intron 5' ss.

Homoeologous non-reciprocal translocation explains a major QTL for seed lignin content in oilseed rape (Brassica napus L.)

A homoeologous non-reciprocal translocation was identified in the major QTL for seed lignin content in the low lignin line SGDH14. The lignin biosynthetic gene PAL4 was deleted. Oilseed rape is a major oil crop and a valuable protein source for animal and human nutrition. Lignin is a non-digestible, major component of the seed coat with negative effect on sensory quality, bioavailability and usage of oilseed rape's protein. Hence, seed lignin reduction is of economic and nutritional importance. In this study, the major QTL for reduced lignin content found on chromosome C05 in the DH population SGDH14 x Express 617 was further examined. SGDH14 had lower seed lignin content than Express 617. Harvested seeds from a F2 population of the same cross were additionally field tested and used for seed quality analysis. The F2 population showed a bimodal distribution for seed lignin content. F2 plants with low lignin content had thinner seed coats compared to high lignin lines. Both groups showed a dark seed colour with a slightly lighter colour in the low lignin group indicating that a low lignin content is not necessarily associated with yellow seed colour. Mapping of genomic long-reads from SGDH14 against the Express 617 genome assembly revealed a homoeologous non-reciprocal translocation (HNRT) in the confidence interval of the major QTL for lignin content. A homologous A05 region is duplicated and replaced the C05 region in SGDH14. As consequence several genes located in the C05 region were lost in SGDH14. Thus, a HNRT was identified in the major QTL region for reduced lignin content in the low lignin line SGDH14. The most promising candidate gene related to lignin biosynthesis on C05, PAL4, was deleted.

Identification of annotation artifacts concerning the chalcone synthase (CHS)

OBJECTIVE

Chalcone synthase (CHS) catalyzes the initial step of the flavonoid biosynthesis. The CHS encoding gene is well studied in numerous plant species. Rapidly growing sequence databases contain hundreds of CHS entries that are the result of automatic annotation. In this study, we evaluated apparent multiplication of CHS domains in CHS gene models of four plant species.

MAIN FINDINGS

CHS genes with an apparent triplication of the CHS domain encoding part were discovered through database searches. Such genes were found in Macadamia integrifolia, Musa balbisiana, Musa troglodytarum, and Nymphaea colorata. A manual inspection of the CHS gene models in these four species with massive RNA-seq data suggests that these gene models are the result of artificial fusions in the annotation process. While there are hundreds of seemingly correct CHS records in the databases, it is not clear why these annotation artifacts appeared.

Genetic mapping and genome-wide association study identify BhYAB4 as the candidate gene regulating seed shape in wax gourd (Benincasa hispida)

Wax gourd is an important vegetable crop of the Cucurbitaceae family. According to the shape and structure of the seed coat, the seeds of the wax gourd can be divided into bilateral and unilateral. Bilateral seeds usually germinate quickly and have a high germination rate than unilateral seeds. Thereby, wax gourd varieties with bilateral seeds are more welcomed by seed companies and growers. However, the genetic basis and molecular mechanism regulating seed shape remain unclear in the wax gourd. In this study, the genetic analysis demonstrated that the seed shape of wax gourd was controlled by a single gene, with bilateral dominant to unilateral. Combined with genetic mapping and genome-wide association study, Bhi04G000544 (BhYAB4), encoding a YABBY transcription factor, was identified as the candidate gene for seed shape determination in the wax gourd. A G/A single nucleotide polymorphism variation of BhYAB4 was detected among different germplasm resources, with BhYAB4^G specifically enriched in bilateral seeds and BhYAB4^A in unilateral seeds. The G to A mutation caused intron retention and premature stop codon of BhYAB4. Expression analysis showed that both BhYAB4^G and BhYAB4^A were highly expressed in seeds, while the nuclear localization of BhYAB4^A protein was disturbed compared with that of BhYAB4^G protein. Finally, a derived cleaved amplified polymorphic sequence marker that could efficiently distinguish between bilateral and unilateral seeds was developed, thereby facilitating the molecular marker-assisted breeding of wax gourd cultivars.

Mapping‑by‑Sequencing Reveals Genomic Regions Associated with Seed Quality Parameters in Brassica napus

Rapeseed (Brassica napus L.) is an important oil crop and has the potential to serve as a highly productive source of protein. This protein exhibits an excellent amino acid composition and has high nutritional value for humans. Seed protein content (SPC) and seed oil content (SOC) are two complex quantitative and polygenic traits which are negatively correlated and assumed to be controlled by additive and epistatic effects. A reduction in seed glucosinolate (GSL) content is desired as GSLs cause a stringent and bitter taste. The goal here was the identification of genomic intervals relevant for seed GSL content and SPC/SOC. Mapping by sequencing (MBS) revealed 30 and 15 new and known genomic intervals associated with seed GSL content and SPC/SOC, respectively. Within these intervals, we identified known but also so far unknown putatively causal genes and sequence variants. A 4 bp insertion in the MYB28 homolog on C09 shows a significant association with a reduction in seed GSL content. This study provides insights into the genetic architecture and potential mechanisms underlying seed quality traits, which will enhance future breeding approaches in B. napus.

ORFeome Phage Display Reveals a Major Immunogenic Epitope on the S2 Subdomain of SARS-CoV-2 Spike Protein

The development of antibody therapies against SARS-CoV-2 remains a challenging task during the ongoing COVID-19 pandemic. All approved therapeutic antibodies are directed against the receptor binding domain (RBD) of the spike, and therefore lose neutralization efficacy against emerging SARS-CoV-2 variants, which frequently mutate in the RBD region. Previously, phage display has been used to identify epitopes of antibody responses against several diseases. Such epitopes have been applied to design vaccines or neutralize antibodies. Here, we constructed an ORFeome phage display library for the SARS-CoV-2 genome. Open reading frames (ORFs) representing the SARS-CoV-2 genome were displayed on the surface of phage particles in order to identify enriched immunogenic epitopes from COVID-19 patients. Library quality was assessed by both NGS and epitope mapping of a monoclonal antibody with a known binding site. The most prominent epitope captured represented parts of the fusion peptide (FP) of the spike. It is associated with the cell entry mechanism of SARS-CoV-2 into the host cell; the serine protease TMPRSS2 cleaves the spike within this sequence. Blocking this mechanism could be a potential target for non-RBD binding therapeutic anti-SARS-CoV-2 antibodies. As mutations within the FP amino acid sequence have been rather rare among SARS-CoV-2 variants so far, this may provide an advantage in the fight against future virus variants.

Case report: Altered pre-mRNA splicing caused by intronic variant c.1499 + 1G > A in the SLC4A4 gene

Proximal renal tubular acidosis (pRTA) with ocular abnormalities is an autosomal recessive disease caused by variants in the Solute Carrier Family 4 Member 4 (SLC4A4) gene. Patients present with metabolic acidosis and low plasma bicarbonate concentration (3∼17 mmol/L). In addition, they are often accompanied by ocular abnormalities, intellectual disability, and growth retardation. The patient underwent whole exome sequencing (WES) and bioinformatics analysis of variant pathogenicity in this study. Then, a minigene assay was conducted to analyze the splicing site variant further. Compound heterozygous variants in the SLC4A4 gene (NM_003759.3), c.145C > T (p.Arg49*) and c.1499 + 1G > A, were detected by WES. The minigene assay showed an mRNA splicing aberration caused by the c.1499 + 1G > A variant. Compared with the wild type, the mutant type caused 4-base insertion between exons 10 and 11 of SLC4A4 after expression in HEK293 cells. In conclusion, the c.1499 + 1G > A variant in the SLC4A4 gene may be one of the genetic causes in the patient. Moreover, our study provides the foundation for future gene therapy of such pathogenic variants.

Spliceator: multi-species splice site prediction using convolutional neural networks

Background

Ab initio prediction of splice sites is an essential step in eukaryotic genome annotation. Recent predictors have exploited Deep Learning algorithms and reliable gene structures from model organisms. However, Deep Learning methods for non-model organisms are lacking.

Results

We developed Spliceator to predict splice sites in a wide range of species, including model and non-model organisms. Spliceator uses a convolutional neural network and is trained on carefully validated data from over 100 organisms. We show that Spliceator achieves consistently high accuracy (89–92%) compared to existing methods on independent benchmarks from human, fish, fly, worm, plant and protist organisms.

Conclusions

Spliceator is a new Deep Learning method trained on high-quality data, which can be used to predict splice sites in diverse organisms, ranging from human to protists, with consistently high accuracy.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04471-3.

Expression of a human cDNA in moss results in spliced mRNAs and fragmentary protein isoforms

Production of biopharmaceuticals relies on the expression of mammalian cDNAs in host organisms. Here we show that the expression of a human cDNA in the moss Physcomitrium patens generates the expected full-length and four additional transcripts due to unexpected splicing. This mRNA splicing results in non-functional protein isoforms, cellular misallocation of the proteins and low product yields. We integrated these results together with the results of our analysis of all 32,926 protein-encoding Physcomitrella genes and their 87,533 annotated transcripts in a web application, physCO, for automatized optimization. A thus optimized cDNA results in about twelve times more protein, which correctly localizes to the ER. An analysis of codon preferences of different production hosts suggests that similar effects occur also in non-plant hosts. We anticipate that the use of our methodology will prevent so far undetected mRNA heterosplicing resulting in maximized functional protein amounts for basic biology and biotechnology.

Large scale genomic rearrangements in selected Arabidopsis thaliana T-DNA lines are caused by T-DNA insertion mutagenesis

BACKGROUND

Experimental proof of gene function assignments in plants is based on mutant analyses. T-DNA insertion lines provided an invaluable resource of mutants and enabled systematic reverse genetics-based investigation of the functions of Arabidopsis thaliana genes during the last decades.

RESULTS

We sequenced the genomes of 14 A. thaliana GABI-Kat T-DNA insertion lines, which eluded flanking sequence tag-based attempts to characterize their insertion loci, with Oxford Nanopore Technologies (ONT) long reads. Complex T-DNA insertions were resolved and 11 previously unknown T-DNA loci identified, resulting in about 2 T-DNA insertions per line and suggesting that this number was previously underestimated. T-DNA mutagenesis caused fusions of chromosomes along with compensating translocations to keep the gene set complete throughout meiosis. Also, an inverted duplication of 800 kbp was detected. About 10 % of GABI-Kat lines might be affected by chromosomal rearrangements, some of which do not involve T-DNA. Local assembly of selected reads was shown to be a computationally effective method to resolve the structure of T-DNA insertion loci. We developed an automated workflow to support investigation of long read data from T-DNA insertion lines. All steps from DNA extraction to assembly of T-DNA loci can be completed within days.

CONCLUSIONS

Long read sequencing was demonstrated to be an effective way to resolve complex T-DNA insertions and chromosome fusions. Many T-DNA insertions comprise not just a single T-DNA, but complex arrays of multiple T-DNAs. It is becoming obvious that T-DNA insertion alleles must be characterized by exact identification of both T-DNA::genome junctions to generate clear genotype-to-phenotype relations.

Gene Fusions Create Partner and Collateral Dependencies Essential to Cancer Cell Survival

Gene fusions frequently result from rearrangements in cancer genomes. In many instances, gene fusions play an important role in oncogenesis; in other instances, they are thought to be passenger events. Although regulatory element rearrangements and copy number alterations resulting from these structural variants are known to lead to transcriptional dysregulation across cancers, the extent to which these events result in functional dependencies with an impact on cancer cell survival is variable. Here we used CRISPR-Cas9 dependency screens to evaluate the fitness impact of 3,277 fusions across 645 cell lines from the Cancer Dependency Map. We found that 35% of cell lines harbored either a fusion partner dependency or a collateral dependency on a gene within the same topologically associating domain as a fusion partner. Fusion-associated dependencies revealed numerous novel oncogenic drivers and clinically translatable alterations. Broadly, fusions can result in partner and collateral dependencies that have biological and clinical relevance across cancer types. SIGNIFICANCE: This study provides insights into how fusions contribute to fitness in different cancer contexts beyond partner-gene activation events, identifying partner and collateral dependencies that may have direct implications for clinical care.

Two Cyc2CL transcripts (Cyc2CL-1 and Cyc2CL-2) may play key roles in the petal and stamen development of ray florets in chrysanthemum

BACKGROUND

Chrysanthemum morifolium is one of the most popular ornamental crops. The capitulum, which is the main ornamental part of chrysanthemum plants, consists of ligulate marginal ray florets, an attractive corolla (petals), and radially hermaphroditic disc florets, but no stamens. In Asteraceae species, the zygomorphic ray florets evolved from the actinomorphic disc florets. During this process, the zygomorphic ligulate corolla arose and the stamens were aborted. Although molecular genetic research has clarified ray floret development to some extent, the precise molecular mechanism underlying ray floret development in chrysanthemum remained unclear.

RESULTS

A CYC2-like gene, Cyc2CL, was cloned from C. morifolium 'Fenditan'. Subsequent analyses revealed that the alternative splicing of Cyc2CL, which occurred in the flower differentiation stage, resulted in the production of Cyc2CL-1 and Cyc2CL-2 in the apical buds. Prior to this stage, only Cyc2CL-1 was produced in the apical buds. A fluorescence in situ hybridization analysis of labeled Cyc2CL-1 and Cyc2CL-2 RNA indicated that Cyc2CL-2 was first expressed in the involucre tissue during the final involucre differentiation stage, but was subsequently expressed in the receptacle and floret primordia as the floral bud differentiation stage progressed. Moreover, Cyc2CL-2 was highly expressed in the inflorescence tissue during the corolla formation stage, and the expression remained high until the end of the floral bud differentiation stage. Furthermore, the overexpression of Cyc2CL-1 and Cyc2CL-2 in transgenic Arabidopsis inhibited stamen and petal development. Therefore, both Cyc2CL-1 and Cyc2CL-2 encode candidate regulators of petal development and stamen abortion and are important for the ray floret development in chrysanthemum.

CONCLUSION

In this study, we characterized the alternatively spliced transcripts of the CYC2-like gene that differ subtly regarding expression and function. The data presented herein will be useful for clarifying the regulatory mechanisms associated with the CYC2-like gene and may also be important for identifying the key genes and molecular mechanisms controlling the development of ray florets in chrysanthemum.

Spliceosomal Introns: Features, Functions, and Evolution

Spliceosomal introns, which have been found in most eukaryotic genes, are non-coding sequences excised from pre-mRNAs by a special complex called spliceosome during mRNA splicing. Introns occur in both protein- and RNA-coding genes and can be found in coding and untranslated gene regions. Because intron sequences vary greatly due to a high rate of polymorphism, the functions of intron had been for a long time associated only with alternative splicing, while intron evolution had been viewed not as an evolution of an individual genomic element, but rather considered within a framework of the evolution of the gene intron-exon structure. Here, we review the theories of intron origin, evolutionary events in the exon-intron structure, such as intron gain, loss, and sliding, intron functions known to date, and mechanisms by which changes in the intron features (length and phase) can affect the regulation of gene-mediated processes.

The reuse of public datasets in the life sciences: potential risks and rewards

The 'big data' revolution has enabled novel types of analyses in the life sciences, facilitated by public sharing and reuse of datasets. Here, we review the prodigious potential of reusing publicly available datasets and the associated challenges, limitations and risks. Possible solutions to issues and research integrity considerations are also discussed. Due to the prominence, abundance and wide distribution of sequencing data, we focus on the reuse of publicly available sequence datasets. We define 'successful reuse' as the use of previously published data to enable novel scientific findings. By using selected examples of successful reuse from different disciplines, we illustrate the enormous potential of the practice, while acknowledging the respective limitations and risks. A checklist to determine the reuse value and potential of a particular dataset is also provided. The open discussion of data reuse and the establishment of this practice as a norm has the potential to benefit all stakeholders in the life sciences.

Automatic Identification of Players in the Flavonoid Biosynthesis with Application on the Biomedicinal Plant Croton tiglium

The flavonoid biosynthesis is a well-characterised model system for specialised metabolism and transcriptional regulation in plants. Flavonoids have numerous biological functions such as UV protection and pollinator attraction, but also biotechnological potential. Here, we present Knowledge-based Identification of Pathway Enzymes (KIPEs) as an automatic approach for the identification of players in the flavonoid biosynthesis. KIPEs combines comprehensive sequence similarity analyses with the inspection of functionally relevant amino acid residues and domains in subjected peptide sequences. Comprehensive sequence sets of flavonoid biosynthesis enzymes and knowledge about functionally relevant amino acids were collected. As a proof of concept, KIPEs was applied to investigate the flavonoid biosynthesis of the medicinal plant Croton tiglium on the basis of a transcriptome assembly. Enzyme candidates for all steps in the biosynthesis network were identified and matched to previous reports of corresponding metabolites in Croton species.

De novo assembly and annotation of transcriptomes from two cultivars of Cannabis sativa with different cannabinoid profiles

Cannabis has been cultivated for millennia for medicinal, industrial and recreational uses. Our long-term goal is to compare the transcriptomes of cultivars with different cannabinoid profiles for therapeutic purposes. Here we describe the de novo assembly, annotation and initial analysis of two cultivars of Cannabis, a high THC variety and a CBD plus THC variety. Cultivars were grown under different lighting conditions; flower buds were sampled over 71 days. Cannabinoid profiles were determined by ESI-LC/MS. RNA samples were sequenced using the HiSeq4000 platform. Transcriptomes were assembled using the DRAP pipeline and annotated using the BLAST2GO pipeline and other tools. Each transcriptome contained over twenty thousand protein encoding transcripts with ORFs and flanking sequence. Identification of transcripts for cannabinoid pathway and related enzymes showed full-length ORFs that align with the draft genomes of the Purple Kush and Finola cultivars. Two transcripts were found for olivetolic acid cyclase (OAC) that mapped to distinct locations on the Purple Kush genome suggesting multiple genes for OAC are expressed in some cultivars. The ability to make high quality annotated reference transcriptomes in Cannabis or other plants can promote rapid comparative analysis between cultivars and growth conditions in Cannabis and other organisms without annotated genome assemblies.

High Contiguity De Novo Genome Sequence Assembly of Trifoliate Yam (Dioscorea dumetorum) Using Long Read Sequencing

Trifoliate yam (Dioscorea dumetorum) is one example of an orphan crop, not traded internationally. Post-harvest hardening of the tubers of this species starts within 24 h after harvesting and renders the tubers inedible. Genomic resources are required for D. dumetorum to improve breeding for non-hardening varieties as well as for other traits. We sequenced the D. dumetorum genome and generated the corresponding annotation. The two haplophases of this highly heterozygous genome were separated to a large extent. The assembly represents 485 Mbp of the genome with an N50 of over 3.2 Mbp. A total of 35,269 protein-encoding gene models as well as 9941 non-coding RNA genes were predicted, and functional annotations were assigned.

Animal, Fungi, and Plant Genome Sequences Harbor Different Non-Canonical Splice Sites

Most protein-encoding genes in eukaryotes contain introns, which are interwoven with exons. Introns need to be removed from initial transcripts in order to generate the final messenger RNA (mRNA), which can be translated into an amino acid sequence. Precise excision of introns by the spliceosome requires conserved dinucleotides, which mark the splice sites. However, there are variations of the highly conserved combination of GT at the 5' end and AG at the 3' end of an intron in the genome. GC-AG and AT-AC are two major non-canonical splice site combinations, which have been known for years. Recently, various minor non-canonical splice site combinations were detected with numerous dinucleotide permutations. Here, we expand systematic investigations of non-canonical splice site combinations in plants across eukaryotes by analyzing fungal and animal genome sequences. Comparisons of splice site combinations between these three kingdoms revealed several differences, such as an apparently increased CT-AC frequency in fungal genome sequences. Canonical GT-AG splice site combinations in antisense transcripts are a likely explanation for this observation, thus indicating annotation errors. In addition, high numbers of GA-AG splice site combinations were observed in Eurytemoraaffinis and Oikopleuradioica. A variant in one U1 small nuclear RNA (snRNA) isoform might allow the recognition of GA as a 5' splice site. In depth investigation of splice site usage based on RNA-Seq read mappings indicates a generally higher flexibility of the 3' splice site compared to the 5' splice site across animals, fungi, and plants.

A Mapped Locus on LG A6 of Brassica juncea Line Tumida Conferring Resistance to White Rust Contains a CNL Type R Gene

White rust, causal agent oomycete Albugo candida, is a significant disease of the cultivated Brassica species. The Indian gene pool lines of oilseed mustard, Brassica juncea, are highly susceptible to the pathogen. Resistance to A. candida has been reported in the east European gene pool lines of mustard and mapped to LG A4 in line Heera and LG A5 in line Donskaja-IV. A new resistance-conferring locus to A. candida isolate AcB1 has been mapped to LG A6 of B. juncea line Tumida-a Chinese vegetable type mustard using an F1DH mapping population that has been developed from a Tumida × Varuna (susceptible Indian gene pool line) cross. A molecular map containing 8,303 genic and GBS markers was used to map the resistance trait to an interval of 63.0 cM-70.8 cM on LG A6. Genome assemblies of Tumida and Varuna were used to find the genes present within the flanking markers discerned by genetic mapping. The most likely candidate gene in the mapped interval is BjuA046215, a CC-NBS-LRR (CNL) type R gene that encodes a protein with all the specific subdomains of the proteins encoded by such genes. Alleles of BjuA046215 in Varuna and other lines of the Indian and the east European gene pools encode proteins that have truncated LRR domains. Analysis of the syntenic regions in some of the Brassicaceae genomes and phylogenetic analysis of CNL type R genes showed BjuA046215 to be closely related to a recently described white rust resistance-conferring R gene BjuWRR1 in B. juncea Donskaja-IV, both belonging to the CNL-D group of R genes. Related R genes in Arabidopsis thaliana confer resistance to another oomycete, Peronospora parasitica.

A chromosome-level sequence assembly reveals the structure of the Arabidopsis thaliana Nd-1 genome and its gene set

In addition to the BAC-based reference sequence of the accession Columbia-0 from the year 2000, several short read assemblies of THE plant model organism Arabidopsis thaliana were published during the last years. Also, a SMRT-based assembly of Landsberg erecta has been generated that identified translocation and inversion polymorphisms between two genotypes of the species. Here we provide a chromosome-arm level assembly of the A. thaliana accession Niederzenz-1 (AthNd-1_v2c) based on SMRT sequencing data. The best assembly comprises 69 nucleome sequences and displays a contig length of up to 16 Mbp. Compared to an earlier Illumina short read-based NGS assembly (AthNd-1_v1), a 75 fold increase in contiguity was observed for AthNd-1_v2c. To assign contig locations independent from the Col-0 gold standard reference sequence, we used genetic anchoring to generate a de novo assembly. In addition, we assembled the chondrome and plastome sequences. Detailed analyses of AthNd-1_v2c allowed reliable identification of large genomic rearrangements between A. thaliana accessions contributing to differences in the gene sets that distinguish the genotypes. One of the differences detected identified a gene that is lacking from the Col-0 gold standard sequence. This de novo assembly extends the known proportion of the A. thaliana pan-genome.