Annotated Protein-coding Genes Research Articles

OysterDB: A Genome Database for Ostreidae.

The molluscan family Ostreidae, commonly known as oysters, is an important molluscan group due to its economic and ecological importance. In recent years, an abundance of genomic data of Ostreidae species has been generated and available in public domain. However, there is still a lack of a high-efficiency database platform to store and distribute these data with comprehensive tools. In this study, we developed an oyster genome database (OysterDB) to consolidate oyster genomic data. This database includes eight oyster genomes and 208,923 protein-coding gene annotations. Bioinformatic tools, such as BLAST and JBrowse, are integrated into the database to provide a user-friendly platform for homologous sequence searching, visualization of genomes, and screen for candidate gene information. Moreover, OysterDB will be continuously updated with ever-growing oyster genomic resources and facilitate future studies for comparative and functional genomic analysis of oysters ( http://oysterdb.com.cn/ ).

Senataxin deficiency disrupts proteostasis through nucleolar ncRNA-driven protein aggregation.

Senataxin is an evolutionarily conserved RNA-DNA helicase involved in DNA repair and transcription termination that is associated with human neurodegenerative disorders. Here, we investigated whether Senataxin loss affects protein homeostasis based on previous work showing R-loop-driven accumulation of DNA damage and protein aggregates in human cells. We find that Senataxin loss results in the accumulation of insoluble proteins, including many factors known to be prone to aggregation in neurodegenerative disorders. These aggregates are located primarily in the nucleolus and are promoted by upregulation of non-coding RNAs expressed from the intergenic spacer region of ribosomal DNA. We also map sites of R-loop accumulation in human cells lacking Senataxin and find higher RNA-DNA hybrids within the ribosomal DNA, peri-centromeric regions, and other intergenic sites but not at annotated protein-coding genes. These findings indicate that Senataxin loss affects the solubility of the proteome through the regulation of transcription-dependent lesions in the nucleus and the nucleolus.

High-quality genome resource of Lasiodiplodia pseudotheobromae associated with die-back on Eucalyptus trees.

Lasiodiplodia pseudotheobromae is an important fungal pathogen associated with die-back, canker and shoot blight in many plant hosts with a wide geographic distribution. The aim of our study was to provide high-quality genome assemblies and sequence annotation resources of L. pseudotheobromae, to facilitate future studies on the systematics, population genetics and genomics of the fungal pathogen L. pseudotheobromae. High-quality genomes of five L. pseudotheobromae isolates were sequenced based on Oxford Nanopore technology (ONT) and Illumina HiSeq sequencing platform. The total size of each assembly ranged from 43Mb to 43.86Mb and over 11,000 protein-coding genes were predicted from each genome. The proteins of predicted genes were annotated using multiple public databases, among the annotated protein-coding genes, more than 4,300 genes were predicted as potential virulence genes by the Pathogen Host Interactions (PHI) database. Moreover, the genome comparative analysis among L. pseudotheobromae and other closely related species revealed that 7,408 gene clusters were shared among them and 152 gene clusters unique to L. pseudotheobromae. This genome and associated datasets provided here will serve as a useful resource for further analyses of this fungal pathogen species.

Chromosomal-level assembly of Tokudaia osimensis, Tokudaia tokunoshimensis, and Tokudaia muenninki genomes

Herein, we present the first high-quality long-read-based chromosome-level genome assemblies and gene annotations of the genomes of three endangered Tokudaia species: Tokudaia osimensis, Tokudaia tokunoshimensis, and Tokudaia muenninki. These species, which are endemic to different islands of the Ryukyu Islands, Japan, exhibited unique karyotypes and sex chromosomal characteristics. The genome assemblies generated using PacBio, Illumina, and Hi-C sequence data consisted of 13 (corresponded to 12 autosomes and one X chromosome), 23 (corresponded to 22 autosomes and one X chromosome), and 23 (corresponded to 21 autosomes and the neo- and ancestral X regions) chromosome-level scaffolds that contained 2,445, 2,477, and 2,661 Mbp of sequence data, respectively. Annotations of protein-coding genes were performed using RNA-Seq-based, homology-based, and Ab initio methods. BUSCO completeness values for every species exceeded 96% for genomes and 98% for genes. These data can be an important resource for contributing to our understanding of species genomes resulting from allopatric speciation and provide insights into mammalian sex-determination mechanisms and sex chromosome evolution.

Reference genome of the nutrition-rich orphan crop chia (Salvia hispanica) and its implications for future breeding.

Chia (Salvia hispanica L.) is one of the most popular nutrition-rich foods and pseudocereal crops of the family Lamiaceae. Chia seeds are a rich source of proteins, polyunsaturated fatty acids (PUFAs), dietary fibers, and antioxidants. In this study, we present the assembly of the chia reference genome, which spans 303.6 Mb and encodes 48,090 annotated protein-coding genes. Our analysis revealed that ~42% of the chia genome harbors repetitive content, and identified ~3 million single nucleotide polymorphisms (SNPs) and 15,380 simple sequence repeat (SSR) marker sites. By investigating the chia transcriptome, we discovered that ~44% of the genes undergo alternative splicing with a higher frequency of intron retention events. Additionally, we identified chia genes associated with important nutrient content and quality traits, such as the biosynthesis of PUFAs and seed mucilage fiber (dietary fiber) polysaccharides. Notably, this is the first report of in-silico annotation of a plant genome for protein-derived small bioactive peptides (biopeptides) associated with improving human health. To facilitate further research and translational applications of this valuable orphan crop, we have developed the Salvia genomics database (SalviaGDB), accessible at https://salviagdb.org.

Identification and preliminary characterization of conserved uncharacterized proteins from Chlamydomonas reinhardtii, Arabidopsis thaliana, and Setaria viridis.

The rapid accumulation of sequenced plant genomes in the past decade has outpaced the still difficult problem of genome-wide protein-coding gene annotation. A substantial fraction of protein-coding genes in all plant genomes are poorly annotated or unannotated and remain functionally uncharacterized. We identified unannotated proteins in three model organisms representing distinct branches of the green lineage (Viridiplantae): Arabidopsis thaliana (eudicot), Setaria viridis (monocot), and Chlamydomonas reinhardtii (Chlorophyte alga). Using similarity searching, we identified a subset of unannotated proteins that were conserved between these species and defined them as Deep Green proteins. Bioinformatic, genomic, and structural predictions were performed to begin classifying Deep Green genes and proteins. Compared to whole proteomes for each species, the Deep Green set was enriched for proteins with predicted chloroplast targeting signals predictive of photosynthetic or plastid functions, a result that was consistent with enrichment for daylight phase diurnal expression patterning. Structural predictions using AlphaFold and comparisons to known structures showed that a significant proportion of Deep Green proteins may possess novel folds. Though only available for three organisms, the Deep Green genes and proteins provide a starting resource of high-value targets for further investigation of potentially new protein structures and functions conserved across the green lineage.

Partner Gene Associated Trinary Fusion and Orchestration of Multiple Molecular Events in Acute Promyelocytic Leukemia with Variant Rara Fusion Genes

Background All-trans retinoic acid (ATRA) has been used with great success in acute promyelocytic leukemia (APL) cases with PML::RARA fusion gene (FG). There are still PML::RARA negative cases manifest as APL. So far, 19 RARA fusion partner genes have been reported, with all reported as 5′ partners ( X::RARA). Intriguingly, a significant proportion of X::RARA cases were resistant to ATRA, and mechanisms are yet to be elucidated. Here we report the finding of a novel type of trinary fusion and the orchestrating of multiple molecular events that confer leukemogenesis and ATRA resistance in X::RARA-APL. Methods A total of 23 cases (Fig. 1A) with RNA-seq and WGS data from 10 centers were enrolled. We used Arriba for routine FG calling and manually investigated the nonconventional fusion sequence. RT-PCR, Sanger sequencing, and optical genome mapping (OGM) were used to validate the FGs and to determine whether the RARA 5′ and 3′ fusion events were located on the same cistron. The responsiveness of the fusion proteins to ATRA was evaluated using an optimized UAS/GAL4 reporter system. Results Among the cases enrolled, there were 17 males and 6 females, aged 3-66 (median 32) years. Morphological and immunophenotype analyses all showed features of APL. We first identified the STAT3::RARA fusion event in an ATRA-resistant case (the index case, #A1). Moreover, we noticed another fusion event of truncated RARA exon 9 and a 3′ partner sequence that did not map to any annotated protein-coding genes. Further analysis determined it was the RC sequence of the transposable element (TE) LTR40a with the annotated loci at 17p11.2. Both OGM (at the genome level) and RT-PCR (Fig. 1A) confirmed that the two RARA fusion events were in the same cistron. Thus, the fusion transcript being STAT3::RARA::LTR40a-RC spliced in tandem, but not separate transcripts of STAT3::RARA and RARA::LTR40a-RC. We named this novel form of tandem splicingtrinary fusion. We, therefore, recruit X::RARA-APL cases to investigate the recurrence of RARA 3′ fusion. We identified 9 RARA 5′ partner genes in 23 cases (Fig. 1B), with TTMV -ORF2, STAT5B, STAT3, IRF2BP1, THRAP3, and TNRC18 in multiple cases. Ten (43.5%) cases also harbor RARA 3′ fusion partners. The partner of case #A1_8 were TEs, and the partner of case #A9_10 was the same gene as their 5′ one. Thirteen cases did not harbor 3′ fusions. Notably, both TNRC18::RARA cases (#15_16) carry mutations on their RARA ligand binding domain (LBD) at the primary onset. Gene expression analysis confirmed the aberrant activation of the involved TEs. All TEs involved in the RARA 3' fusions confer a poly_A signal sequence, which was essential for a mRNA. In-depth analysis indicated that TEs participate in forming RARA-FGs through a transposition mechanism rather than a translocation mechanism. The molecular events addition to RARA 5′ splicing showed a 5′ partner associated pattern: all of the 4 STAT3::RARA and 4 STAT5B::RARA cases harbor TEs as their 3′ partner, both TNRC18::RARA cases carry LBD mutation, while all 7 TTMV -ORF2::RARA cases did not harbor 3′ splicing or LBD mutation. Remarkably, although the RARA 3' splicing sites varied, they all resulted in partial or complete exon 9 truncation. RARA exon 9 encodes helix 11_12 (H11_12) of its LBD, which is pivotal in allosteric response to ATRA. The RARA trinary fusion will result in LBD-H12 or H11_12 castration due to the 3' splicing, then lead to losing responsiveness to ATRA via an allosteric disability mechanism (Fig. 1C, 1D). Experiments confirmed that cells transfected with STAT3::RARA fusion with intact LBD responded well to ATRA, while both the fusion in case #A1 and artificial fusion with LBD-H12 castration showed no response to physiological and pharmacological concentrations of ATRA (Fig. 1E). Conclusions This is the first report of pathological trinary fusion, the TEs and transposition mechanism participating in composing pathological fusion genes, and the RARA LBD castration confers ATRA resistance in X::RARA-APL. Also, there might be another synergistic mechanism of RARA 5' fusion and LBD mutation. The formation of trinary fusion and gain of LBD mutation, in addition to X::RARA fusion, requires more genome splicing or transposition or mutation events, which explains the rarity of X::RARA-APL. This study highlights mechanisms of partner gene-associated trinary fusion and the orchestrating of multiple molecular events in the leukemogenesis in X::RARA-APL cases.

Cold stress induces differential gene expression of retained homeologs in Camelina sativa cv Suneson.

Camelina sativa (L.) Crantz, a member of the Brassicaceae, has potential as a biofuel feedstock which is attributable to the production of fatty acids in its seeds, its fast growth cycle, and low input requirements. While a genome assembly is available for camelina, it was generated from short sequence reads and is thus highly fragmented in nature. Using long read sequences, we generated a chromosome-scale, highly contiguous genome assembly (644,491,969 bp) for the spring biotype cultivar 'Suneson' with an N50 contig length of 12,031,512 bp and a scaffold N50 length of 32,184,682 bp. Annotation of protein-coding genes revealed 91,877 genes that encode 133,355 gene models. We identified a total of 4,467 genes that were significantly up-regulated under cold stress which were enriched in gene ontology terms associated with "response to cold" and "response to abiotic stress". Coexpression analyses revealed multiple coexpression modules that were enriched in genes differentially expressed following cold stress that had putative functions involved in stress adaptation, specifically within the plastid. With access to a highly contiguous genome assembly, comparative analyses with Arabidopsis thaliana revealed 23,625 A. thaliana genes syntenic with 45,453 Suneson genes. Of these, 24,960 Suneson genes were syntenic to 8,320 A. thaliana genes reflecting a 3 camelina homeolog to 1 Arabidopsis gene relationship and retention of all three homeologs. Some of the retained triplicated homeologs showed conserved gene expression patterns under control and cold-stressed conditions whereas other triplicated homeologs displayed diverged expression patterns revealing sub- and neo-functionalization of the homeologs at the transcription level. Access to the chromosome-scale assembly of Suneson will enable both basic and applied research efforts in the improvement of camelina as a sustainable biofuel feedstock.

The status of the human gene catalogue.

Scientists have been trying to identify every gene in the human genome since the initial draft was published in 2001. In the years since, much progress has been made in identifying protein-coding genes, currently estimated to number fewer than 20,000, with an ever-expanding number of distinct protein-coding isoforms. Here we review the status of the human gene catalogue and the efforts to complete it in recent years. Beside the ongoing annotation of protein-coding genes, their isoforms and pseudogenes, the invention of high-throughput RNA sequencing and other technological breakthroughs have led to a rapid growth in the number of reported non-coding RNA genes. For most of these non-coding RNAs, the functional relevance is currently unclear; we look at recent advances that offer paths forward to identifying their functions and towards eventually completing the human gene catalogue. Finally, we examine the need for a universal annotation standard that includes all medically significant genes and maintains their relationships with different reference genomes for the use of the human gene catalogue in clinical settings.

Galba: genome annotation with miniprot and AUGUSTUS

BackgroundThe Earth Biogenome Project has rapidly increased the number of available eukaryotic genomes, but most released genomes continue to lack annotation of protein-coding genes. In addition, no transcriptome data is available for some genomes.ResultsVarious gene annotation tools have been developed but each has its limitations. Here, we introduce GALBA, a fully automated pipeline that utilizes miniprot, a rapid protein-to-genome aligner, in combination with AUGUSTUS to predict genes with high accuracy. Accuracy results indicate that GALBA is particularly strong in the annotation of large vertebrate genomes. We also present use cases in insects, vertebrates, and a land plant. GALBA is fully open source and available as a docker image for easy execution with Singularity in high-performance computing environments.ConclusionsOur pipeline addresses the critical need for accurate gene annotation in newly sequenced genomes, and we believe that GALBA will greatly facilitate genome annotation for diverse organisms.

An insight into the functional genomics and species classification of Eudiplozoon nipponicum (Monogenea, Diplozoidae), a haematophagous parasite of the common carp Cyprinus carpio

BackgroundMonogenea (Platyhelminthes, Neodermata) are the most species-rich class within the Neodermata superclass of primarily fish parasites. Despite their economic and ecological importance, monogenean research tends to focus on their morphological, phylogenetic, and population characteristics, while comprehensive omics analyses aimed at describing functionally important molecules are few and far between. We present a molecular characterisation of monogenean representative Eudiplozoon nipponicum, an obligate haematophagous parasite infecting the gills of the common carp. We report its nuclear and mitochondrial genomes, present a functional annotation of protein molecules relevant to the molecular and biochemical aspect of physiological processes involved in interactions with the fish hosts, and re-examinate the taxonomic position of Eudiplozoon species within the Diplozoidae family.ResultsWe have generated 50.81 Gbp of raw sequencing data (Illumina and Oxford Nanopore reads), bioinformatically processed, and de novo assembled them into a genome draft 0.94 Gbp long, consisting of 21,044 contigs (N50 = 87 kbp). The final assembly represents 57% of the estimated total genome size (~ 1.64 Gbp), whereby repetitive and low-complexity regions account for ~ 64% of the assembled length. In total, 36,626 predicted genes encode 33,031 proteins and homology-based annotation of protein-coding genes (PCGs) and proteins characterises 14,785 (44.76%) molecules. We have detected significant representation of functional proteins and known molecular functions. The numbers of peptidases and inhibitors (579 proteins), characterised GO terms (16,016 unique assigned GO terms), and identified KEGG Orthology (4,315 proteins) acting in 378 KEGG pathways demonstrate the variety of mechanisms by which the parasite interacts with hosts on a macromolecular level (immunomodulation, feeding, and development). Comparison between the newly assembled E. nipponicum mitochondrial genome (length of 17,038 bp) and other diplozoid monogeneans confirms the existence of two distinct Eudiplozoon species infecting different fish hosts: Cyprinus carpio and Carassius spp.ConclusionsAlthough the amount of sequencing data and characterised molecules of monogenean parasites has recently increased, a better insight into their molecular biology is needed. The E. nipponicum nuclear genome presented here, currently the largest described genome of any monogenean parasite, represents a milestone in the study of monogeneans and their molecules but further omics research is needed to understand these parasites’ biological nature.

Genome Sequence Resource for Colletotrichum gloeosporioides, an Important Pathogenic Fungus Causing Anthracnose of Dioscorea alata.

Anthracose disease is one of the most important diseases of Dioscorea alata and many other food yams, which is caused by Colletotrichum gloeosporioides fugus from the Glomerellaceae family of the Sordariomycetes class. In the present study, a C. gloeosporioides starin named CgDa01 was isolated from D. alata, and its genome was sequenced based on Oxford Nanopore technology (ONT) and the Illumina sequencing platform. The high-quality genome of CgDa01 was assembled with a 62.78 Mb genome size and 15,845 predicted protein-coding genes. The proteins of predicted genes were annotated using multiple public databases, including the nonredundant protein database, the InterProScan databases, and Kyoto Encyclopedia of Genes and Genomes. Among the annotated protein-coding genes, 55 were predicted as potential virulence genes by the fungal virulence factor database. The C. gloeosporioides CgDa01 genome assembly described in this study can serve as a resource for better understanding the pathogenic mechanism of C. gloeosporioides on yam hosts.

Assembly of the Complete Mitochondrial Genome of Gelsemium elegans Revealed the Existence of Homologous Conformations Generated by a Repeat Mediated Recombination.

Gelsemium elegans (G. elegans) is a Chinese medicinal plant with substantial economic and feeding values. There is a lack of detailed studies on the mitochondrial genome of G. elegans. In this study, the mitochondrial genome of G. elegans was sequenced and assembled, and its substructure was investigated. The mitochondrial genome of G. elegans is represented by two circular chromosomes of 406,009 bp in length with 33 annotated protein-coding genes, 15 tRNA genes, and three rRNA genes. We detected 145 pairs of repeats and found that four pairs of repeats could mediate the homologous recombination into one major conformation and five minor conformations, and the presence of conformations was verified by PCR amplification and Sanger sequencing. A total of 124 SSRs were identified in the G. elegans mitochondrial genome. The homologous segments between the chloroplast and mitochondrial genomes accounted for 5.85% of the mitochondrial genome. We also predicted 477 RNA potential editing sites and found that the nad4 gene was edited 38 times, which was the most frequent occurrence. Taken together, the mitochondrial genome of G. elegans was assembled and annotated. We gained a more comprehensive understanding on the genome of this medicinal plant, which is vital for its effective utilization and genetic improvement, especially for cytoplasmic male sterility breeding and evolution analysis in G. elegans.

DeepTSS: multi-branch convolutional neural network for transcription start site identification from CAGE data

BackgroundThe widespread usage of Cap Analysis of Gene Expression (CAGE) has led to numerous breakthroughs in understanding the transcription mechanisms. Recent evidence in the literature, however, suggests that CAGE suffers from transcriptional and technical noise. Regardless of the sample quality, there is a significant number of CAGE peaks that are not associated with transcription initiation events. This type of signal is typically attributed to technical noise and more frequently to random five-prime capping or transcription bioproducts. Thus, the need for computational methods emerges, that can accurately increase the signal-to-noise ratio in CAGE data, resulting in error-free transcription start site (TSS) annotation and quantification of regulatory region usage. In this study, we present DeepTSS, a novel computational method for processing CAGE samples, that combines genomic signal processing (GSP), structural DNA features, evolutionary conservation evidence and raw DNA sequence with Deep Learning (DL) to provide single-nucleotide TSS predictions with unprecedented levels of performance.ResultsTo evaluate DeepTSS, we utilized experimental data, protein-coding gene annotations and computationally-derived genome segmentations by chromatin states. DeepTSS was found to outperform existing algorithms on all benchmarks, achieving 98% precision and 96% sensitivity (accuracy 95.4%) on the protein-coding gene strategy, with 96.66% of its positive predictions overlapping active chromatin, 98.27% and 92.04% co-localized with at least one transcription factor and H3K4me3 peak.ConclusionsCAGE is a key protocol in deciphering the language of transcription, however, as every experimental protocol, it suffers from biological and technical noise that can severely affect downstream analyses. DeepTSS is a novel DL-based method for effectively removing noisy CAGE signal. In contrast to existing software, DeepTSS does not require feature selection since the embedded convolutional layers can readily identify patterns and only utilize the important ones for the classification task.This study highlights the key role that DL can play in Molecular Biology, by removing the inherent flaws of experimental protocols, that form the backbone of contemporary research. Here, we show how DeepTSS can unleash the full potential of an already popular and mature method such as CAGE, and push the boundaries of coding and non-coding gene expression regulator research even further.

A chromosome-level genome of Brachymystax tsinlingensis provides resources and insights into salmonids evolution.

Brachymystax tsinlingensis Li, 1966 is an endangered freshwater fish with economic, ecological, and scientific values. Study of the genome of B. tsinlingensis might be particularly insightful given that this is the only Brachymystax species with genome. We present a high-quality chromosome-level genome assembly and protein-coding gene annotation for B. tsinlingensis with Illumina short reads, Nanopore long reads, Hi-C sequencing reads, and RNA-seq reads from 5 tissues/organs. The final chromosome-level genome size is 2,031,709,341 bp with 40 chromosomes. We found that the salmonids have a unique GC content and codon usage, have a slower evolutionary rate, and possess specific positively selected genes. We also confirmed the salmonids have undergone a whole-genome duplication event and a burst of transposon-mediated repeat expansion, and lost HoxAbβ Hox cluster, highly expressed genes in muscle may partially explain the migratory habits of B. tsinlingensis. The high-quality B. tsinlingensis assembled genome could provide a valuable reference for the study of other salmonids as well as aid the conservation of this endangered species.

Fine mapping and characterization of the Crinkled Dwarf gene in cotton

Cotton (Gossypium) is an important natural fiber crop worldwide. Plant height, which is an important commercial cotton breeding trait, is also useful for plant architecture research. Crinkled Dwarf (cr), a natural mutant of Sea Island cotton (Gossypium barbadense L.), simultaneously displays crinkled leaves and dwarfing. However, the genes controlling cr have not been characterized, and the molecular mechanism of the trait remains unknown. In this study, genetic analysis indicated that cr was likely controlled by a major recessive gene, although the cr phenotype was sensitive to environmental factors, including light intensity. Using bulked-segregant analysis and recombinant methods, the cr gene was mapped to a 35.8-kbp genomic interval on chromosome D03 that contained three annotated protein-coding genes. Comparative transcriptome analysis showed that cr enhanced stress/disease resistance and repressed plant development/growth. Endogenous phytohormone measurements showed that most auxin compounds were decreased in the cr mutant. Although the cr mutation was insensitive to exogenous auxin, gibberellin and cytokinin treatments promoted and repressed main stem elongation, respectively. However, no hormone treatment could fully rescue the crinkled leaf phenotype. These findings will contribute to the breeding of cotton lines with improved plant architecture or yield, and aid in the dissection of the genetic basis of cr.

Bioinformatic approach for the discovery of cis-eQTL signals during fruit ripening of a woody species as grape (Vitis vinifera L.).

Expression quantitative trait loci (eQTLs) are associations between genetic variants, such as Single Nucleotide Polymorphisms (SNPs), and gene expression. eQTLs are an important tool to understand the genetic variance of gene expression of complex phenotypes. eQTLs analyses are common in biomedical models but are scarce in woody crop species such as fruit trees or grapes. In this study, a comprehensive bioinformatic analysis was conducted leveraging with expression data from two different growth stages, around ripening onset, of 10 genotypes of grape (Vitis vinifera L.). A total of 2170 cis-eQTL were identified in 212 gene modulated at ripening onset. The 48% of these DEGs have a known function. Among the annotated protein-coding genes, terpene synthase, auxin-regulatory factors, GRFS, ANK_REP_REGION domain-containing protein, Kinesin motor domain-containing protein and flavonol synthase were noted. This new inventory of cis-eQTLs influencing gene expression during fruit ripening will be an important resource to examine variation for this trait and will help to elucidate the complex genetic architecture underlying this process in grape.

Annotation of a Newly Sequenced Genomic Region in Drosophila ananassae

Understanding evolutionary changes across Drosophila species can lead to a better understanding of how genomes have changed over time. Drosophila species have existed for millions of years, allowing the opportunity to study and identify the functions of genes orthologous to human genes. The etiology of numerous diseases can be better elucidated through the utilization of bioinformatic tools and genomic annotation. The Genome Education Partnership (GEP) is a nationwide community that provides undergraduate students the opportunity to gain experience in annotating eukaryotic genomes. Through the use of bioinformatics tools such as UCSC Genome Browser, NCBI Blast, and Flybase, students participate in the F‐element and Pathways projects. The F‐element project involves the annotation of protein‐coding genes located on chromosome four, commonly known as the Muller F element. Over time, the F element in some species of Drosophila, including Drosophila ananassae, has been expanding. Comparative genome annotation analyzing changes in the organization of this genomic region may lead to a better understanding of the complexity of these genomes. The Pathways project involves annotating specific genes involved in the Insulin Signaling Pathway. These genes play a critical role in the regulation of expression and metabolism in many eukaryotes. In this study, the annotation of a genomic region in Drosophila ananassae is discussed.A 25,000 base pair region of chromosome 3L in D. ananassae (Muller D element) was annotated for comparison to the F element. In this region, known as Contig 21, five protein‐coding genes were identified: CG10222, CG10713, CG32121, CG33263, and flr. These genes had high protein alignment (70% to 90% amino acid identity) compared to Drosophila melanogaster, however, the flr gene possessed some interesting characteristics. Upon the annotation of the flr gene in Drosophila ananassae, it was predicted that one of the isoforms was no longer expressed and/or functional. The gene flr encodes the Drosophila Actin Interacting Protein 1, which interacts with the gene product of trs (cofillin) to promote F‐actin disassembly. In Drosophila melanogaster, the flr gene encodes two isoforms, however in Drosophila ananassae, only the flr‐PA isoform appears to be expressed and/or functional. A deletion/insertion of a nucleotide within exon 1 in the flr gene leads to a frameshift resulting in a non‐functional protein due to a change in the original stop codon of isoform B in D. ananassae. The gene flr plays an important role in the development of the organism, therefore molecular analysis was utilized to confirm the annotation model above. RNA was extracted from Drosophila ananassae, and reverse transcription PCR (RT‐PCR) was utilized to amplify this region of the genome. Dideoxy sequencing was performed to analyze the expression of different isoforms of this gene. The sequencing data was compared to the Drosophila ananassae and Drosophila melanogaster genomes using BLAST and the UCSC Genome Browser. This data supported the expression of only one isoform of the flr gene in D. ananassae. This data also suggests that one of the flr gene isoforms is highly conserved in this species. Additional analysis of newly sequenced genomic data such as this will provide additional understanding on how genomes such as these have changed over time.

Dominant transcript expression profiles of human protein-coding genes interrogated with GTEx dataset

The discovery and quantification of mRNA transcripts using short-read next-generation sequencing (NGS) data is a complicated task. There are far more alternative mRNA transcripts expressed by human genes than can be identified from NGS transcriptome data and various bioinformatic pipelines, while the numbers of annotated human protein-coding genes has gradually declined in recent years. It is essential to learn more about the thorough tissue expression profiles of alternative transcripts in order to obtain their molecular modulations and actual functional significance. In this report, we present a bioinformatic database for interrogating the representative tissue of human protein-coding transcripts. The database allows researchers to visually explore the top-ranked transcript expression profiles in particular tissue types. Most transcripts of protein-coding genes were found to have certain tissue expression patterns. This observation demonstrated that many alternative transcripts were particularly modulated in different cell types. This user-friendly tool visually represents transcript expression profiles in a tissue-specific manner. Identification of tissue specific protein-coding genes and transcripts is a substantial advance towards interpreting their biological functions and further functional genomics studies.

A chromosome-level genome of the human blood fluke Schistosoma japonicum identifies the genomic basis of host-switching.

The evolution and adaptation of S.japonicum, a zoonotic parasite that causes human schistosomiasis, remain unclear because of the lack of whole-genome data. We construct a chromosome-level S.japonicum genome and analyze it together with 72 samples representing six populations of the entire endemic region. We observe a Taiwan zoophilic lineage splitting from zoonotic populations ∼45,000 years ago, consistent with the divergent history of their intermediate hosts. Interestingly, we detect a severe population bottleneck in S.japonicum, largely coinciding with human history in Asia during the last glacial maximum. We identify several genomic regions underlying natural selection, including GATAD2A and Lmln, both showing remarkable differentiation among different areas. RNAi knockdown suggests association of GATAD2A with parasite development and infection in definitive hosts, while Lmln relates to the specificity of the intermediate hosts. Our study provides insights into the evolution of S.japonicum and serves as a resource for further studies.