nf-core/isoseq: simple gene and isoform annotation with PacBio Iso-Seq long-read sequencing.

MOTIVATION: Iso-Seq RNA long-read sequencing enables the identification of full-length transcripts and isoforms, removing the need for complex analysis such as transcriptome assembly. However, the raw sequencing data need to be processed in a series of steps before annotation is complete. Here, we present nf-core/isoseq, a pipeline for automatic read processing and genome annotation. Following nf-core guidelines, the pipeline has few dependencies and can be run on any of platforms. AVAILABILITY AND IMPLEMENTATION: The pipeline is freely available online on the nf-core website (https://nf-co.re/isoseq) and on GitHub (https://github.com/nf-core/isoseq) under MIT License (DOI: 10.5281/zenodo.7116979).

Avian Immunome DB: an example of a user-friendly interface for extracting genetic information.

BACKGROUND: Genomic and genetic studies often require a target list of genes before conducting any hypothesis testing or experimental verification. With the ever-growing number of sequenced genomes and a variety of different annotation strategies, comes the potential for ambiguous gene symbols, making it cumbersome to capture the "correct" set of genes. In this article, we present and describe the Avian Immunome DB (AVIMM) for easy gene property extraction as exemplified by avian immune genes. The avian immune system is characterised by a cascade of complex biological processes underlaid by more than 1000 different genes. It is a vital trait to study particularly in birds considering that they are a significant driver in spreading zoonotic diseases. With the completion of phase II of the B10K ("Bird 10,000 Genomes") consortium's whole-genome sequencing effort, we have included 363 annotated bird genomes in addition to other publicly available bird genome data which serve as a valuable foundation for AVIMM. CONSTRUCTION AND CONTENT: A relational database with avian immune gene evidence from Gene Ontology, Ensembl, UniProt and the B10K consortium has been designed and set up. The foundation stone or the "seed" for the initial set of avian immune genes is based on the well-studied model organism chicken (Gallus gallus). Gene annotations, different transcript isoforms, nucleotide sequences and protein information, including amino acid sequences, are included. Ambiguous gene names (symbols) are resolved within the database and linked to their canonical gene symbol. AVIMM is supplemented by a command-line interface and a web front-end to query the database. UTILITY AND DISCUSSION: The internal mapping of unique gene symbol identifiers to canonical gene symbols allows for an ambiguous gene property search. The database is organised within core and feature tables, which makes it straightforward to extend for future purposes. The database design is ready to be applied to other taxa or biological processes. Currently, the database contains 1170 distinct avian immune genes with canonical gene symbols and 612 synonyms across 363 bird species. While the command-line interface readily integrates into bioinformatics pipelines, the intuitive web front-end with download functionality offers sophisticated search functionalities and tracks the origin for each record. AVIMM is publicly accessible at https://avimm.ab.mpg.de .

Illuminating the dark side of the human transcriptome with long read transcript sequencing.

BACKGROUND: The human transcriptome annotation is regarded as one of the most complete of any eukaryotic species. However, limitations in sequencing technologies have biased the annotation toward multi-exonic protein coding genes. Accurate high-throughput long read transcript sequencing can now provide additional evidence for rare transcripts and genes such as mono-exonic and non-coding genes that were previously either undetectable or impossible to differentiate from sequencing noise. RESULTS: We developed the Transcriptome Annotation by Modular Algorithms (TAMA) software to leverage the power of long read transcript sequencing and address the issues with current data processing pipelines. TAMA achieved high sensitivity and precision for gene and transcript model predictions in both reference guided and unguided approaches in our benchmark tests using simulated Pacific Biosciences (PacBio) and Nanopore sequencing data and real PacBio datasets. By analyzing PacBio Sequel II Iso-Seq sequencing data of the Universal Human Reference RNA (UHRR) using TAMA and other commonly used tools, we found that the convention of using alignment identity to measure error correction performance does not reflect actual gain in accuracy of predicted transcript models. In addition, inter-read error correction can cause major changes to read mapping, resulting in potentially over 6 K erroneous gene model predictions in the Iso-Seq based human genome annotation. Using TAMA's genome assembly based error correction and gene feature evidence, we predicted 2566 putative novel non-coding genes and 1557 putative novel protein coding gene models. CONCLUSIONS: Long read transcript sequencing data has the power to identify novel genes within the highly annotated human genome. The use of parameter tuning and extensive output information of the TAMA software package allows for in depth exploration of eukaryotic transcriptomes. We have found long read data based evidence for thousands of unannotated genes within the human genome. More development in sequencing library preparation and data processing are required for differentiating sequencing noise from real genes in long read RNA sequencing data.

Normalized long read RNA sequencing in chicken reveals transcriptome complexity similar to human.

BACKGROUND: Despite the significance of chicken as a model organism, our understanding of the chicken transcriptome is limited compared to human. This issue is common to all non-human vertebrate annotations due to the difficulty in transcript identification from short read RNAseq data. While previous studies have used single molecule long read sequencing for transcript discovery, they did not perform RNA normalization and 5'-cap selection which may have resulted in lower transcriptome coverage and truncated transcript sequences. RESULTS: We sequenced normalised chicken brain and embryo RNA libraries with Pacific Bioscience Iso-Seq. 5' cap selection was performed on the embryo library to provide methodological comparison. From these Iso-Seq sequencing projects, we have identified 60 k transcripts and 29 k genes within the chicken transcriptome. Of these, more than 20 k are novel lncRNA transcripts with ~3 k classified as sense exonic overlapping lncRNA, which is a class that is underrepresented in many vertebrate annotations. The relative proportion of alternative transcription events revealed striking similarities between the chicken and human transcriptomes while also providing explanations for previously observed genomic differences. CONCLUSIONS: Our results indicate that the chicken transcriptome is similar in complexity compared to human, and provide insights into other vertebrate biology. Our methodology demonstrates the potential of Iso-Seq sequencing to rapidly expand our knowledge of transcriptomics.

Predictors of Clergy's Ability to Fulfill a Suicide Prevention Gatekeeper Role.

Catholic, Jewish and Protestant clergy (n = 801) completed a survey to identify predictors of clergy's ability to fulfill a suicide gatekeeper role. Exploratory backward stepwise regression identified predictors of risk identification including suicide knowledge, religion, conducting suicide funerals, having an attitude that people have a right to die, age, and race. Predictors of ability to intervene include suicide knowledge, training, religion, right to die attitude, and ethnicity. Recommendations include more suicide training and clergy self-care.

Evolution of the avian ß-defensin and cathelicidin genes.

BACKGROUND: ß-defensins and cathelicidins are two families of cationic antimicrobial peptides (AMPs) with a broad range of antimicrobial activities that are key components of the innate immune system. Due to their important roles in host defense against rapidly evolving pathogens, the two gene families provide an ideal system for studying adaptive gene evolution. In this study we performed phylogenetic and selection analyses on ß-defensins and cathelicidins from 53 avian species representing 32 orders to examine the evolutionary dynamics of these peptides in birds. RESULTS AND CONCLUSIONS: Avian ß-defensins are found in a gene cluster consisting of 13 subfamiles. Nine of these are conserved as one to one orthologs in all birds, while the others (AvBD1, AvBD3, AvBD7 and AvBD14) are more subject to gene duplication or pseudogenisation events in specific avian lineages. Avian cathelicidins are found in a gene cluster consisting of three subfamilies with species-specific duplications and gene loss. Evidence suggested that the propiece and mature peptide domains of avian cathelicidins are possibly co-evolving in such a way that the cationicity of the mature peptide is partially neutralised by the negative charge of the propiece prior to peptide secretion (further evidence obtained by repeating the analyses on primate cathelicidins). Negative selection (overall mean dN < dS) was detected in most of the gene domains examined, conserving certain amino acid residues that may be functionally crucial for the avian ß-defensins and cathelicidins, while episodic positive selection was also involved in driving the diversification of specific codon sites of certain AMPs in avian evolutionary history. These findings have greatly improved our understanding of the molecular evolution of avian AMPs and will be useful to understand their role in the avian innate immune response. Additionally, the large dataset of ß-defensin and cathelicidin peptides may also provide a valuable resource for translational research and development of novel antimicrobial agents in the future.

Detection and characterization of small insertion and deletion genetic variants in modern layer chicken genomes.

BACKGROUND: Small insertions and deletions (InDels) constitute the second most abundant class of genetic variants and have been found to be associated with many traits and diseases. The present study reports on the detection and characterisation of about 883 K high quality InDels from the whole-genome analysis of several modern layer chicken lines from diverse breeds. RESULTS: To reduce the error rates seen in InDel detection, this study used the consensus set from two InDel-calling packages: SAMtools and Dindel, as well as stringent post-filtering criteria. By analysing sequence data from 163 chickens from 11 commercial and 5 experimental layer lines, this study detected about 883 K high quality consensus InDels with 93% validation rate and an average density of 0.78 InDels/kb over the genome. Certain chromosomes, viz, GGAZ, 16, 22 and 25 showed very low densities of InDels whereas the highest rate was observed on GGA6. In spite of the higher recombination rates on microchromosomes, the InDel density on these chromosomes was generally lower relative to macrochromosomes possibly due to their higher gene density. About 43-87% of the InDels were found to be fixed within each line. The majority of detected InDels (86%) were 1-5 bases and about 63% were non-repetitive in nature while the rest were tandem repeats of various motif types. Functional annotation identified 613 frameshift, 465 non-frameshift and 10 stop-gain/loss InDels. Apart from the frameshift and stopgain/loss InDels that are expected to affect the translation of protein sequences and their biological activity, 33% of the non-frameshift were predicted as evolutionary intolerant with potential impact on protein functions. Moreover, about 2.5% of the InDels coincided with the most-conserved elements previously mapped on the chicken genome and are likely to define functional elements. InDels potentially affecting protein function were found to be enriched for certain gene-classes e.g. those associated with cell proliferation, chromosome and Golgi organization, spermatogenesis, and muscle contraction. CONCLUSIONS: The large catalogue of InDels presented in this study along with their associated information such as functional annotation, estimated allele frequency, etc. are expected to serve as a rich resource for application in future research and breeding in the chicken.

Genome-wide analysis reveals the extent of EAV-HP integration in domestic chicken.

BACKGROUND: EAV-HP is an ancient retrovirus pre-dating Gallus speciation, which continues to circulate in modern chicken populations, and led to the emergence of avian leukosis virus subgroup J causing significant economic losses to the poultry industry. We mapped EAV-HP integration sites in Ethiopian village chickens, a Silkie, Taiwan Country chicken, red junglefowl Gallus gallus and several inbred experimental lines using whole-genome sequence data. RESULTS: An average of 75.22 ± 9.52 integration sites per bird were identified, which collectively group into 279 intervals of which 5 % are common to 90 % of the genomes analysed and are suggestive of pre-domestication integration events. More than a third of intervals are specific to individual genomes, supporting active circulation of EAV-HP in modern chickens. Interval density is correlated with chromosome length (P < 2.31(-6)), and 27 % of intervals are located within 5 kb of a transcript. Functional annotation clustering of genes reveals enrichment for immune-related functions (P < 0.05). CONCLUSIONS: Our results illustrate a non-random distribution of EAV-HP in the genome, emphasising the importance it may have played in the adaptation of the species, and provide a platform from which to extend investigations on the co-evolutionary significance of endogenous retroviral genera with their hosts.

Relapsing Polychondritis Presenting with Meningoencephalitis and Dementia: Correlation with Neuroimaging and Clinical Features.

PURPOSE: Relapsing polychondritis (RP) is a rare systemic autoimmune disease affecting cartilaginous and non-cartilaginous structures. Neurological involvement is rarer but results in profound disability. Early identification and treatment of underlying RP may promote neurological recovery. CASE REPORT: We illustrated a 53-year-old man diagnosed with dementia. Neuroimaging and cerebrospinal fluid studies disclosed meningoencephalitis. "Prominent ear sign" was evident on diffusion-weight magnetic resonance imaging. After glucocortisone administration, the improvement of clinical manifestations was closely correlated subsequent neuroimaging findings. CONCLUSION: The importance of better understanding of this disease in terms of the prevention of further tissue damage in patients with RP cannot be overemphasized.

The molecular basis of differential host responses to avian influenza viruses in avian species with differing susceptibility.

Introduction: Highly pathogenic avian influenza (HPAI) viruses, such as H5N1, continue to pose a serious threat to animal agriculture, wildlife and to public health. Controlling and mitigating this disease in domestic birds requires a better understanding of what makes some species highly susceptible (such as turkey and chicken) while others are highly resistant (such as pigeon and goose). Susceptibility to H5N1 varies both with species and strain; for example, species that are tolerant of most H5N1 strains, such as crows and ducks, have shown high mortality to emerging strains in recent years. Therefore, in this study we aimed to examine and compare the response of these six species, to low pathogenic avian influenza (H9N2) and two strains of H5N1 with differing virulence (clade 2.2 and clade 2.3.2.1) to determine how susceptible and tolerant species respond to HPAI challenge. Methods: Birds were challenged in infection trials and samples (brain, ileum and lung) were collected at three time points post infection. The transcriptomic response of birds was examined using a comparative approach, revealing several important discoveries. Results: We found that susceptible birds had high viral loads and strong neuro-inflammatory response in the brain, which may explain the neurological symptoms and high mortality rates exhibited following H5N1 infection. We discovered differential regulation of genes associated with nerve function in the lung and ileum, with stronger differential regulation in resistant species. This has intriguing implications for the transmission of the virus to the central nervous system (CNS) and may also indicate neuro-immune involvement at the mucosal surfaces. Additionally, we identified delayed timing of the immune response in ducks and crows following infection with the more deadly H5N1 strain, which may account for the higher mortality in these species caused by this strain. Lastly, we identified candidate genes with potential roles in susceptibility/resistance which provide excellent targets for future research. Discussion: This study has helped elucidate the responses underlying susceptibility to H5N1 influenza in avian species, which will be critical in developing sustainable strategies for future control of HPAI in domestic poultry.

Corrigendum: The molecular basis of differential host responses to avian influenza viruses in avian species with differing susceptibility.

[This corrects the article DOI: 10.3389/fcimb.2023.1067993.].

The swan genome and transcriptome, it is not all black and white.

BACKGROUND: The Australian black swan (Cygnus atratus) is an iconic species with contrasting plumage to that of the closely related northern hemisphere white swans. The relative geographic isolation of the black swan may have resulted in a limited immune repertoire and increased susceptibility to infectious diseases, notably infectious diseases from which Australia has been largely shielded. Unlike mallard ducks and the mute swan (Cygnus olor), the black swan is extremely sensitive to highly pathogenic avian influenza. Understanding this susceptibility has been impaired by the absence of any available swan genome and transcriptome information. RESULTS: Here, we generate the first chromosome-length black and mute swan genomes annotated with transcriptome data, all using long-read based pipelines generated for vertebrate species. We use these genomes and transcriptomes to show that unlike other wild waterfowl, black swans lack an expanded immune gene repertoire, lack a key viral pattern-recognition receptor in endothelial cells and mount a poorly controlled inflammatory response to highly pathogenic avian influenza. We also implicate genetic differences in SLC45A2 gene in the iconic plumage of the black swan. CONCLUSION: Together, these data suggest that the immune system of the black swan is such that should any avian viral infection become established in its native habitat, the black swan would be in a significant peril.

The Role of Dendritic Cells in the Host Response to Marek's Disease Virus (MDV) as Shown by Transcriptomic Analysis of Susceptible and Resistant Birds.

Despite the successful control of highly contagious tumorigenic Marek's disease (MD) by vaccination, a continuous increase in MD virus (MDV) virulence over recent decades has put emphasis on the development of more MD-resistant chickens. The cell types and genes involved in resistance therefore need to be recognized. The virus is primarily lymphotropic, but research should also focus on innate immunity, as innate immune cells are among the first to encounter MDV. Our previous study on MDV-macrophage interaction revealed significant differences between MHC-congenic lines 61 (MD-resistant) and 72 (MD-susceptible). To investigate the role of dendritic cells (DCs) in MD resistance, bone-marrow-derived DCs from these lines were infected with MDV in vitro. They were then characterized by cell sorting, and the respective transcriptomes analysed by RNA-seq. The differential expression (DE) of genes revealed a strong immune activation in DCs of the susceptible line, although an inherent immune supremacy was shown by the resistant line, including a significant expression of tumour-suppressor miRNA, gga-mir-124a, in line 61 control birds. Enrichment analysis of DE genes revealed high expression of an oncogenic transcription factor, AP-1, in the susceptible line following MDV challenge. This research highlights genes and pathways that may play a role in DCs in determining resistance or susceptibility to MDV infection.

A high-resolution single-molecule sequencing-based Arabidopsis transcriptome using novel methods of Iso-seq analysis.

BACKGROUND: Accurate and comprehensive annotation of transcript sequences is essential for transcript quantification and differential gene and transcript expression analysis. Single-molecule long-read sequencing technologies provide improved integrity of transcript structures including alternative splicing, and transcription start and polyadenylation sites. However, accuracy is significantly affected by sequencing errors, mRNA degradation, or incomplete cDNA synthesis. RESULTS: We present a new and comprehensive Arabidopsis thaliana Reference Transcript Dataset 3 (AtRTD3). AtRTD3 contains over 169,000 transcripts-twice that of the best current Arabidopsis transcriptome and including over 1500 novel genes. Seventy-eight percent of transcripts are from Iso-seq with accurately defined splice junctions and transcription start and end sites. We develop novel methods to determine splice junctions and transcription start and end sites accurately. Mismatch profiles around splice junctions provide a powerful feature to distinguish correct splice junctions and remove false splice junctions. Stratified approaches identify high-confidence transcription start and end sites and remove fragmentary transcripts due to degradation. AtRTD3 is a major improvement over existing transcriptomes as demonstrated by analysis of an Arabidopsis cold response RNA-seq time-series. AtRTD3 provides higher resolution of transcript expression profiling and identifies cold-induced differential transcription start and polyadenylation site usage. CONCLUSIONS: AtRTD3 is the most comprehensive Arabidopsis transcriptome currently. It improves the precision of differential gene and transcript expression, differential alternative splicing, and transcription start/end site usage analysis from RNA-seq data. The novel methods for identifying accurate splice junctions and transcription start/end sites are widely applicable and will improve single-molecule sequencing analysis from any species.

Clergy as Suicide Prevention Gatekeepers.

801 U.S. Catholic, Jewish and Protestant clergy reported on their suicide gatekeeping activities. Using vignettes, they identified suicide risk and selected interventions for three risk levels. Two-thirds of the sample who provide counseling reported at least one contact from a suicidal person per year. Clergy were significantly more concurrent with experts in identifying risk and selecting interventions with high risk but deviated more from the experts with low and medium risk. Most reported needing more training.

A high-quality genome and comparison of short- versus long-read transcriptome of the palaearctic duck Aythya fuligula (tufted duck).

BACKGROUND: The tufted duck is a non-model organism that experiences high mortality in highly pathogenic avian influenza outbreaks. It belongs to the same bird family (Anatidae) as the mallard, one of the best-studied natural hosts of low-pathogenic avian influenza viruses. Studies in non-model bird species are crucial to disentangle the role of the host response in avian influenza virus infection in the natural reservoir. Such endeavour requires a high-quality genome assembly and transcriptome. FINDINGS: This study presents the first high-quality, chromosome-level reference genome assembly of the tufted duck using the Vertebrate Genomes Project pipeline. We sequenced RNA (complementary DNA) from brain, ileum, lung, ovary, spleen, and testis using Illumina short-read and Pacific Biosciences long-read sequencing platforms, which were used for annotation. We found 34 autosomes plus Z and W sex chromosomes in the curated genome assembly, with 99.6% of the sequence assigned to chromosomes. Functional annotation revealed 14,099 protein-coding genes that generate 111,934 transcripts, which implies a mean of 7.9 isoforms per gene. We also identified 246 small RNA families. CONCLUSIONS: This annotated genome contributes to continuing research into the host response in avian influenza virus infections in a natural reservoir. Our findings from a comparison between short-read and long-read reference transcriptomics contribute to a deeper understanding of these competing options. In this study, both technologies complemented each other. We expect this annotation to be a foundation for further comparative and evolutionary genomic studies, including many waterfowl relatives with differing susceptibilities to avian influenza viruses.

Erratum: Chakraborty, P. et al. Macrophages from Susceptible and Resistant Chicken Lines have Different Transcriptomes following Marek's Disease Virus Infection. Genes 2019, 10, 74.

The authors wish to make the following correction to their paper published in Genes [...].

Macrophages from Susceptible and Resistant Chicken Lines have Different Transcriptomes following Marek's Disease Virus Infection.

Despite successful control by vaccination, Marek's disease (MD) has continued evolving to greater virulence over recent years. To control MD, selection and breeding of MD-resistant chickens might be a suitable option. MHC-congenic inbred chicken lines, 61 and 72, are highly resistant and susceptible to MD, respectively, but the cellular and genetic basis for these phenotypes is unknown. Marek's disease virus (MDV) infects macrophages, B-cells, and activated T-cells in vivo. This study investigates the cellular basis of resistance to MD in vitro with the hypothesis that resistance is determined by cells active during the innate immune response. Chicken bone marrow-derived macrophages from lines 61 and 72 were infected with MDV in vitro. Flow cytometry showed that a higher percentage of macrophages were infected in line 72 than in line 61. A transcriptomic study followed by in silico functional analysis of differentially expressed genes was then carried out between the two lines pre- and post-infection. Analysis supports the hypothesis that macrophages from susceptible and resistant chicken lines display a marked difference in their transcriptome following MDV infection. Resistance to infection, differential activation of biological pathways, and suppression of oncogenic potential are among host defense strategies identified in macrophages from resistant chickens.