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).

Diversity of endogenous avian leukosis virus subgroup E (ALVE) insertions in indigenous chickens.

BACKGROUND: Avian leukosis virus subgroup E (ALVE) insertions are endogenous retroviruses (ERV) that are restricted to the domestic chicken and its wild progenitor. In commercial chickens, ALVE are known to have a detrimental effect on productivity and provide a source for recombination with exogenous retroviruses. The wider diversity of ALVE in non-commercial chickens and the role of these elements in ERV-derived immunity (EDI) are yet to be investigated. RESULTS: In total, 974 different ALVE were identified from 407 chickens sampled from village populations in Ethiopia, Iraq, and Nigeria, using the recently developed obsERVer bioinformatics identification pipeline. Eighty-eight percent of all identified ALVE were novel, bringing the known number of ALVE integrations to more than 1300 across all analysed chickens. ALVE content was highly lineage-specific and populations generally exhibited a large diversity of ALVE at low frequencies, which is typical for ERV involved in EDI. A significantly larger number of ALVE was found within or near coding regions than expected by chance, although a relative depletion of ALVE was observed within coding regions, which likely reflects selection against deleterious integrations. These effects were less pronounced than in previous analyses of chickens from commercial lines. CONCLUSIONS: Identification of more than 850 novel ALVE has trebled the known diversity of these retroviral elements. This work provides the basis for future studies to fully quantify the role of ALVE in immunity against exogenous ALV, and development of programmes to improve the productivity and welfare of chickens in developing economies.

Transcriptome profiling of granulosa and theca cells during dominant follicle development in the horse.

Several aspects of equine ovarian physiology are unique among domestic species. Moreover, follicular growth patterns are very similar between horses and humans. This study aimed to characterize, for the first time, global gene expression profiles associated with growth and preovulatory (PO) maturation of equine dominant follicles. Granulosa cells (GCs) and theca interna cells (TCs) were harvested from follicles (n = 5) at different stages of an ovulatory wave in mares corresponding to early dominance (ED; diameter ≥22 mm), late dominance (LD; ≥33 mm) and PO stage (34 h after administration of crude equine gonadotropins at LD stage), and separately analyzed on a horse gene expression microarray, followed by validation using quantitative PCR and immunoblotting/immunohistochemistry. Numbers of differentially expressed transcripts (DETs; ≥2-fold; P < 0.05) during the ED-LD and LD-PO transitions were 546 and 2419 in GCs and 5 and 582 in TCs. The most prominent change in GCs was the down-regulation of transcripts associated with cell division during both ED-LD and LD-PO. In addition, DET sets during LD-PO in GCs were enriched for genes involved in cell communication/adhesion, antioxidation/detoxification, immunity/inflammation, and cholesterol biosynthesis. In contrast, the largest change in TCs during the LD-PO transition was an up-regulation of genes involved in immune activation, with other DET sets mapping to GPCR/cAMP signaling, lipid/amino acid metabolism, and cell proliferation/survival and differentiation. In conclusion, distinct expression profiles were identified between growing and PO follicles and, particularly, between GCs and TCs within each stage. Several DETs were identified that have not been associated with follicle development in other species.

A chromosome-level genome assembly of a free-living white-crowned sparrow (Zonotrichia leucophrys gambelii).

The white-crowned sparrow, Zonotrichia leucophrys, is a passerine bird with a wide distribution and it is extensively adapted to environmental changes. It has historically acted as a model species in studies on avian ecology, physiology and behaviour. Here, we present a high-quality chromosome-level genome of Zonotrichia leucophrys using PacBio and OmniC sequencing data. Gene models were constructed by combining RNA-seq and Iso-seq data from liver, hypothalamus, and ovary. In total a 1,123,996,003 bp genome was generated, including 31 chromosomes assembled in complete scaffolds along with other, unplaced scaffolds. This high-quality genome assembly offers an important genomic resource for the research community using the white-crowned sparrow as a model for understanding avian genome biology and development, and provides a genomic basis for future studies, both fundamental and applied.

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.

Immune response of rats vaccinated orally with various plant-expressed recombinant cysteine proteinase constructs when challenged with Fasciola hepatica metacercariae.

BACKGROUND: Cysteine proteinases of Fasciola hepatica are important candidates for vaccine antigens because of their role in fluke biology and host-parasite relationships. In our previous experiments, we found that a recombinant cysteine proteinase cloned from adult F. hepatica (CPFhW) can protect rats against liver fluke infections when it is administered intramuscularly or intranasally in the form of cDNA. We also observed considerable protection upon challenge following mucosal vaccination with inclusion bodies containing recombinant CPFhW produced in Escherichia coli. In this study, we explore oral vaccination, which may be the desired method of delivery and is potentially capable of preventing infections at the site of helminth entry. To provide antigen encapsulation and to protect the vaccine antigen from degradation in the intestinal tract, transgenic plant-based systems are used. METHODOLOGY: In the present study, we aimed to evaluate the protective ability of mucosal vaccinations of 12-week-old rats with CPFhW produced in a transgenic-plant-based system. To avoid inducing tolerance and to maximise the immune response induced by oral immunisation, we used the hepatitis B virus (HBV) core protein (HBcAg) as a carrier. Animals were immunised with two doses of the antigen and challenged with 25 or 30 metacercariae of F. hepatica. CONCLUSIONS: We obtained substantial protection after oral administration of the plant-produced hybrids of CPFhW and HBcAg. The highest level of protection (65.4%) was observed in animals immunised with transgenic plants expressing the mature CPFhW enzyme flanked by Gly-rich linkers and inserted into c/e1 epitope of truncated HBcAg. The immunised rats showed clear IgG1 and IgM responses to CPFhW for 4 consecutive weeks after the challenge.

Binary Switching of Calendar Cells in the Pituitary Defines the Phase of the Circannual Cycle in Mammals.

Persistent free-running circannual (approximately year-long) rhythms have evolved in animals to regulate hormone cycles, drive metabolic rhythms (including hibernation), and time annual reproduction. Recent studies have defined the photoperiodic input to this rhythm, wherein melatonin acts on thyrotroph cells of the pituitary pars tuberalis (PT), leading to seasonal changes in the control of thyroid hormone metabolism in the hypothalamus. However, seasonal rhythms persist in constant conditions in many species in the absence of a changing photoperiod signal, leading to the generation of circannual cycles. It is not known which cells, tissues, and pathways generate these remarkable long-term rhythmic processes. We show that individual PT thyrotrophs can be in one of two binary states reflecting either a long (EYA3(+)) or short (CHGA(+)) photoperiod, with the relative proportion in each state defining the phase of the circannual cycle. We also show that a morphogenic cycle driven by the PT leads to extensive re-modeling of the PT and hypothalamus over the circannual cycle. We propose that the PT may employ a recapitulated developmental pathway to drive changes in morphology of tissues and cells. Our data are consistent with the hypothesis that the circannual timer may reside within the PT thyrotroph and is encoded by a binary switch timing mechanism, which may regulate the generation of circannual neuroendocrine rhythms, leading to dynamic re-modeling of the hypothalamic interface. In summary, the PT-ventral hypothalamus now appears to be a prime structure involved in long-term rhythm generation.

The variable part of the dnaK gene as an alternative marker for phylogenetic studies of rhizobia and related alpha Proteobacteria.

DnaK is the 70 kDa chaperone that prevents protein aggregation and supports the refolding of damaged proteins. Due to sequence conservation and its ubiquity this chaperone has been widely used in phylogenetic studies. In this study, we applied the less conserved part that encodes the so-called alpha-subdomain of the substrate-binding domain of DnaK for phylogenetic analysis of rhizobia and related non-symbiotic alpha-Proteobacteria. A single 330 bp DNA fragment was routinely amplified from DNA templates isolated from the species of the genera, Azorhizobium, Bradyrhizobium, Mesorhizobium, Rhizobium and Sinorhizobium, but also from some non-symbiotic alpha Proteobacteria such as Blastochloris, Chelatobacter and Chelatococcus. Phylogenetic analyses revealed high congruence between dnaK sequences and 16S rDNA trees, but they were not identical. In contrast, the partition homogeneity tests revealed that dnaK sequence data could be combined with other housekeeping genes such as recA, atpD or glnA. The dnaK trees exhibited good resolution in the cases of the genera Mesorhizobium, Sinorhizobium and Rhizobium, even better than usually shown by 16S rDNA phylogeny. The dnaK phylogeny supported the close phylogenetic relationship of Rhizobium galegae and Agrobacterium tumefaciens (R. radiobacter) C58, which together formed a separate branch within the fast-growing rhizobia, albeit closer to the genus Sinorhizobium. The Rhizobium and Sinorhizobium genera carried an insertion composed of two amino acids, which additionally supported the phylogenetic affinity of these two genera, as well as their distinctness from the Mesorhizobium genus. Consistently with the phylogeny shown by 16S-23S rDNA intergenic region sequences, the dnaK trees divided the genus Bradyrhizobium into three main lineages, corresponding to B. japonicum, B. elkanii, and photosynthetic Bradyrhizobium strains that infect Aeschynomene plants. Our results suggest that the 330 bp dnaK sequences could be used as an additional taxonomic marker for rhizobia and related species (alternatively to the 16S rRNA gene phylogeny).

Identification of Eya3 and TAC1 as long-day signals in the sheep pituitary.

Seasonally breeding mammals such as sheep use photoperiod, encoded by the nocturnal secretion of the pineal hormone melatonin, as a critical cue to drive hormone rhythms and synchronize reproduction to the most optimal time of year. Melatonin acts directly on the pars tuberalis (PT) of the pituitary, regulating expression of thyrotropin, which then relays messages back to the hypothalamus to control reproductive circuits. In addition, a second local intrapituitary circuit controls seasonal prolactin (PRL) release via one or more currently uncharacterized low-molecular-weight peptides, termed "tuberalins," of PT origin. Studies in birds have identified the transcription factor Eya3 as the first molecular response activated by long photoperiod (LP). Using arrays and in situ hybridization studies, we demonstrate here that Eya3 is the strongest LP-activated gene in sheep, revealing a common photoperiodic molecular response in birds and mammals. We also demonstrate TAC1 (encoding the tachykinins substance P and neurokinin A) to be strongly activated by LP within the sheep PT. We show that these PRL secretagogues act on primary pituitary cells and thus are candidates for the elusive PT-expressed tuberalin seasonal hormone regulator.

Immunoprotective properties of transgenic plants expressing E2 glycoprotein from CSFV and cysteine protease from Fasciola hepatica.

Immune responses were elicited in laboratory animals after oral vaccination by transgenic plants (lettuce and alfalfa) expressing the E2 glycoprotein of Classical Swine Fever Virus (CSFV) or cysteine protease from Fasciola hepatica. ELISA analyses demonstrated that the oral route is effective in inducing a specific antibody response against these antigens in mice.

Low sequence similarity and gene content of symbiotic clusters of Bradyrhizobium sp. WM9 (Lupinus) indicate early divergence of "lupin" lineage in the genus Bradyrhizobium.

Two sequenced nodulation regions of lupin Bradyrhizobium sp. WM9 carried the majority of genes involved in the Nod factor production. The nod region I harbored: nolA, nodD, nodA, nodB, nodC, nodS, nodI, nodJ, nolO, nodZ, fixR, nifA, fixA, nodM, nolK and noeL. This gene arrangement resembled that found in the nodulation region of Bradyrhizobium japonicum USDA110, however strain WM9 harbored only one nodD gene copy, while the nodM, nolK and noeL genes had no counterparts in the 410 kb symbiotic region of strain USDA110. Region II harbored nolL and nodW, but lacked an nodV gene. Both regions carried ORFs that lacked similarity to the published USDA110 sequences, though they had homologues in symbiotic regions of Rhizobium etli, Sinorhizobium sp. NGR234 and Mesorhizobium loti. These differences in gene content, as well as a low average sequence identity (70%) of symbiotic genes with respect to B. japonicum USDA110 were in contrast with the phylogenetic relationship of USDA110 and WM9 revealed by the analysis of 16S rDNA and dnaK sequences. This most likely reflected an early divergence of symbiotic loci, and possible co-speciation with distinct legumes. During this process the loss of a noel gene and the acquisition of a nolL gene could be regarded as an adaptation towards these legumes that responded to Nod factors carrying 4-O-acetylfucose rather than 2-O-methylfucose. This explained various responses of lupins and serradella plants to infection by mutants in nodZ and nolL genes, knowing that serradella is a stringent legume while lupins are more promiscuous legumes.