The identification of immune genes in the milk transcriptome of the Tasmanian devil (Sarcophilus harrisii)

Bioactive components in the marsupial pouch and milk

Marsupials give birth to immunologically naïve young after a relatively short gestation period compared with eutherians. Consequently, the joey relies significantly on maternal protection, which is the focus of the present review. The milk and the pouch environment are essential contributors to maternal protection for the healthy development of joeys. In this review, we discuss bioactive components found in the marsupial pouch and milk that form cornerstones of maternal protection. These bioactive components include immune cells, immunoglobulins, the S100 family of calcium-binding proteins, lysozymes, whey proteins, antimicrobial peptides and other immune proteins. Furthermore, we investigated the possibility of the presence of plurifunctional components in milk and pouches that are potentially bioactive. These compounds include caseins, vitamins and minerals, oligosaccharides, lipids and microRNAs. Where applicable, this review addresses variability in bioactive components during different phases of lactation, designed to fulfil the immunological needs of the growing pouch young. Yet, there are numerous additional research opportunities to pursue, including uncovering novel bioactive components and investigating their modes of action, dynamics, stability and ability to penetrate the gut epithelium to facilitate systemic effects.

Marsupial immune protection is shaped by enhancer sharing and gene cluster duplication of cathelicidin antimicrobial peptides

Marsupial neonates are born with immature immune systems, making them vulnerable to pathogens. While neonates receive maternal protection, they can also independently combat pathogens, though the mechanisms remain unknown. Using the sugar glider (Petaurus breviceps) as a model, we investigated immunological defense strategies of marsupial neonates. Cathelicidins, a family of antimicrobial peptides expanded in the genomes of marsupials, are highly expressed in developing neutrophils. Sugar glider cathelicidins reside in two genomic clusters and their coordinated expression is achieved by enhancer sharing within clusters and long-range physical interactions between clusters. These cathelicidins modulate immune responses and have potent antimicrobial effects, sufficient to provide protection in a mouse model of sepsis. Lastly, cathelicidins have a complex evolutionary history, where marsupials and monotremes are the only tetrapods that retained two cathelicidin clusters. Thus, cathelicidins are critical mediators of marsupial immunity, and their evolution reflects the life history-specific immunological needs of these animals.

Characterisation of defensins across the marsupial family tree

Defensins are antimicrobial peptides involved in innate immunity, and gene number differs amongst eutherian mammals. Few studies have investigated defensins in marsupials, despite their potential involvement in immunological protection of altricial young. Here we use recently sequenced marsupial genomes and transcriptomes to annotate defensins in nine species across the marsupial family tree. We characterised 35 alpha and 286 beta defensins; gene number differed between species, although Dasyuromorphs had the largest repertoire. Defensins were encoded in three gene clusters within the genome, syntenic to eutherians, and were expressed in the pouch and mammary gland. Marsupial beta defensins were closely related to eutherians, however marsupial alpha defensins were more divergent. We identified marsupial orthologs of human DEFB3 and 6, and several marsupial-specific beta defensin lineages which may have novel functions. Marsupial predicted mature peptides were highly variable in length and sequence composition. We propose candidate peptides for future testing to elucidate the function of marsupial defensins.

The evolution of two transmissible cancers in Tasmanian devils

Tasmanian devils have spawned two transmissible cancer lineages, named devil facial tumor 1 (DFT1) and devil facial tumor 2 (DFT2). We investigated the genetic diversity and evolution of these clones by analyzing 78 DFT1 and 41 DFT2 genomes relative to a newly assembled, chromosome-level reference. Time-resolved phylogenetic trees reveal that DFT1 first emerged in 1986 (1982 to 1989) and DFT2 in 2011 (2009 to 2012). Subclone analysis documents transmission of heterogeneous cell populations. DFT2 has faster mutation rates than DFT1 across all variant classes, including substitutions, indels, rearrangements, transposable element insertions, and copy number alterations, and we identify a hypermutated DFT1 lineage with defective DNA mismatch repair. Several loci show plausible evidence of positive selection in DFT1 or DFT2, including loss of chromosome Y and inactivation of MGA, but none are common to both cancers. This study reveals the parallel long-term evolution of two transmissible cancers inhabiting a common niche in Tasmanian devils.

Best genome sequencing strategies for annotation of complex immune gene families in wildlife

BACKGROUND

The biodiversity crisis and increasing impact of wildlife disease on animal and human health provides impetus for studying immune genes in wildlife. Despite the recent boom in genomes for wildlife species, immune genes are poorly annotated in nonmodel species owing to their high level of polymorphism and complex genomic organisation. Our research over the past decade and a half on Tasmanian devils and koalas highlights the importance of genomics and accurate immune annotations to investigate disease in wildlife. Given this, we have increasingly been asked the minimum levels of genome quality required to effectively annotate immune genes in order to study immunogenetic diversity. Here we set out to answer this question by manually annotating immune genes in 5 marsupial genomes and 1 monotreme genome to determine the impact of sequencing data type, assembly quality, and automated annotation on accurate immune annotation.

RESULTS

Genome quality is directly linked to our ability to annotate complex immune gene families, with long reads and scaffolding technologies required to reassemble immune gene clusters and elucidate evolution, organisation, and true gene content of the immune repertoire. Draft-quality genomes generated from short reads with HiC or 10× Chromium linked reads were unable to achieve this. Despite mammalian BUSCOv5 scores of up to 94.1% amongst the 6 genomes, automated annotation pipelines incorrectly annotated up to 59% of manually annotated immune genes regardless of assembly quality or method of automated annotation.

CONCLUSIONS

Our results demonstrate that long reads and scaffolding technologies, alongside manual annotation, are required to accurately study the immune gene repertoire of wildlife species.

A reference genome for the critically endangered woylie, Bettongia penicillata ogilbyi

Biodiversity is declining globally, and Australia has one of the worst extinction records for mammals. The development of sequencing technologies means that genomic approaches are now available as important tools for wildlife conservation and management. Despite this, genome sequences are available for only 5% of threatened Australian species. Here we report the first reference genome for the woylie (Bettongia penicillata ogilbyi), a critically endangered marsupial from Western Australia, and the first genome within the Potoroidae family. The woylie reference genome was generated using Pacific Biosciences HiFi long-reads, resulting in a 3.39 Gbp assembly with a scaffold N50 of 6.49 Mbp and 86.5% complete mammalian BUSCOs. Assembly of a global transcriptome from pouch skin, tongue, heart and blood RNA-seq reads was used to guide annotation with Fgenesh++, resulting in the annotation of 24,655 genes. The woylie reference genome is a valuable resource for conservation, management and investigations into disease-induced decline of this critically endangered marsupial.

Opossum Mammary Maturation as It Relates to Immune Cell Infiltration and Nutritional Gene Transcription

The mammary gland has evolved to accommodate the developmental needs of offspring in species-specific ways. This is particularly true for marsupials. Marsupial milk content changes dramatically throughout lactation in ways appearing timed with neonatal ontogeny and behavior. Here we investigate morphological restructuring within the mammaries throughout lactation in the gray short-tailed opossum, Monodelphis domestica. Substantial remodeling of the mammaries occurs throughout the first half of active lactation. It is not until the latter half of lactation that opossum mammaries appear histologically similar to active eutherian mammaries. Noteworthy was the presence of eosinophils in early developing mammary tissue, which correlated with elevated abundance of transcripts encoding the chemokine IL-16. The presence and abundance of whey protein transcripts within the opossum mammaries were also quantified. Whey acidic protein (WAP) transcript abundance peaked in the latter half of lactation and remained elevated through weaning. Minimal transcripts for the marsupial-specific Early and Late Lactation Proteins (ELP/LLP) were detected during active lactation. Elevated abundance of LLP transcripts was only detected prior to parturition. Overall, the results support the role of eosinophils in mammary restructuring appearing early in mammalian evolution, and describe key similarities and differences in nutritional protein transcript abundance among marsupial species.

Didelphis albiventris: an overview of unprecedented transcriptome sequencing of the white-eared opossum

BACKGROUND

The white-eared opossum (Didelphis albiventris) is widely distributed throughout Brazil and South America. It has been used as an animal model for studying different scientific questions ranging from the restoration of degraded green areas to medical aspects of Chagas disease, leishmaniasis and resistance against snake venom. As a marsupial, D. albiventris can also contribute to the understanding of the molecular mechanisms that govern the different stages of organogenesis. Opossum joeys are born after only 13 days, and the final stages of organogenesis occur when the neonates are inside the pouch, depending on lactation. As neither the genome of this opossum species nor its transcriptome has been completely sequenced, the use of D. albiventris as an animal model is limited. In this work, we sequenced the D. albiventris transcriptome by RNA-seq to obtain the first catalogue of differentially expressed (DE) genes and gene ontology (GO) annotations during the neonatal stages of marsupial development.

RESULTS

The D. albiventris transcriptome was obtained from whole neonates harvested at birth (P0), at 5 days of age (P5) and at 10 days of age (P10). The de novo assembly of these transcripts generated 85,338 transcripts. Approximately 30% of these transcripts could be mapped against the amino acid sequences of M. domestica, the evolutionarily closest relative of D. albiventris to be sequenced thus far. Among the expressed transcripts, 2077 were found to be DE between P0 and P5, 13,780 between P0 and P10, and 1453 between P5 and P10. The enriched GO terms were mainly related to the immune system, blood tissue development and differentiation, vision, hearing, digestion, the CNS and limb development.

CONCLUSIONS

The elucidation of opossum transcriptomes provides an out-group for better understanding the distinct characteristics associated with the evolution of mammalian species. This study provides the first transcriptome sequences and catalogue of genes for a marsupial species at different neonatal stages, allowing the study of the mechanisms involved in organogenesis.

Transcriptomics of Tasmanian Devil (Sarcophilus Harrisii) Ear Tissue Reveals Homogeneous Gene Expression Patterns across a Heterogeneous Landscape

In an era of unprecedented global change, exploring patterns of gene expression among wild populations across their geographic range is crucial for characterizing adaptive potential. RNA-sequencing studies have successfully characterized gene expression differences among populations experiencing divergent environmental conditions in a wide variety of taxa. However, few of these studies have identified transcriptomic signatures to multivariate, environmental stimuli among populations in their natural environments. Herein, we aim to identify environmental and sex-driven patterns of gene expression in the Tasmanian devil (Sarcophilus harrisii), a critically endangered species that occupies a heterogeneous environment. We performed RNA-sequencing on ear tissue biopsies from adult male and female devils from three populations at the extremes of their geographic range. There were no transcriptome-wide patterns of differential gene expression that would be suggestive of significant, environmentally-driven transcriptomic responses. The general lack of transcriptome-wide variation in gene expression levels across the devil's geographic range is consistent with previous studies that documented low levels of genetic variation in the species. However, genes previously implicated in local adaptation to abiotic environment in devils were enriched for differentially expressed genes. Additionally, three modules of co-expressed genes were significantly associated with either population of origin or sex.

Transcriptome sequencing of the long-nosed bandicoot (Perameles nasuta) reveals conservation and innovation of immune genes in the marsupial order Peramelemorphia

Bandicoots are omnivorous marsupials of the order Peramelemorphia. Conservation concerns and their unique biological characteristics suggest peramelomorphs are worthy research subjects, but knowledge of their genetics and immunology has lagged behind that of other high-profile marsupials. Here, we characterise the transcriptome of the long-nose bandicoot (Perameles nasuta), the first high-throughput data set from any peramelomorph. We investigate the immune gene repertoire of the bandicoot, with a focus on key immune gene families, and compare to previously characterised marsupial and mammalian species. We find that the immune gene complement in bandicoot is often conserved with respect to other marsupials; however, the diversity of expressed transcripts in several key families, such as major histocompatibility complex, T cell receptor μ and natural killer cell receptors, appears greater in the bandicoot than other Australian marsupials, including devil and koala. This transcriptome is an important first step for future studies of bandicoots and the bilby, allowing for population level analysis and construction of bandicoot-specific immunological reagents and assays. Such studies will be critical to understanding the immunology and physiology of Peramelemorphia and to inform the conservation of these unique marsupials.

Antimicrobial Protection of Marsupial Pouch Young

Marsupials diverged from eutherian mammals about 148 million years ago and represent a unique lineage of mammals with distinctive morphological and reproductive characteristics. Marsupials have significantly shorter gestation periods than eutherians. Pregnancy typically ranges from 15 to 35 days, with young being born at a very early developmental stage and lacking differentiated lymphoid tissues and mature effector cells. Recent microbiome studies of the marsupial pouch revealed that marsupial young can face intense microbial challenges after birth, as the pouch contains a broad range of Gram-positive and Gram-negative bacteria. Antimicrobials are believed to play a significant role in the immune protection of marsupial newborns during their pouch life. The skin of the post-reproductive pouch secretes antimicrobial lysozyme and dermcidin, which may contribute to the decreased density of certain bacteria in the pouch. A range of antimicrobial agents, such as immunoglobulins, lysozyme, transferrin, and cathelicidins, have been identified in marsupial milk. Antimicrobial assays have revealed that marsupial cathelicidins have broad-spectrum activity against a variety of bacteria and fungi, including several multi-drug resistant strains. In this article, we will review the action mechanisms of these antimicrobial compounds and discuss how they protect marsupial newborns from potentially pathogenic bacteria inside the pouch. We will also discuss the potential of marsupial antimicrobial compounds as a source of novel antibiotics.

Characterization of the antimicrobial peptide family defensins in the Tasmanian devil (Sarcophilus harrisii), koala (Phascolarctos cinereus), and tammar wallaby (Macropus eugenii)

Defensins comprise a family of cysteine-rich antimicrobial peptides with important roles in innate and adaptive immune defense in vertebrates. We characterized alpha and beta defensin genes in three Australian marsupials: the Tasmanian devil (Sarcophilus harrisii), koala (Phascolarctos cinereus), and tammar wallaby (Macropus eugenii) and identified 48, 34, and 39 defensins, respectively. One hundred and twelve have the classical antimicrobial peptides characteristics required for pathogen membrane targeting, including cationic charge (between 1+ and 15+) and a high proportion of hydrophobic residues (>30%). Phylogenetic analysis shows that gene duplication has driven unique and species-specific expansions of devil, koala, and tammar wallaby beta defensins and devil alpha defensins. Defensin genes are arranged in three genomic clusters in marsupials, whereas further duplications and translocations have occurred in eutherians resulting in four and five gene clusters in mice and humans, respectively. Marsupial defensins are generally under purifying selection, particularly residues essential for defensin structural stability. Certain hydrophobic or positively charged sites, predominantly found in the defensin loop, are positively selected, which may have functional significance in defensin-target interaction and membrane insertion.

Cathelicidins in the Tasmanian devil (Sarcophilus harrisii)

Tasmanian devil joeys, like other marsupials, are born at a very early stage of development, prior to the development of their adaptive immune system, yet survive in a pathogen-laden pouch and burrow. Antimicrobial peptides, called cathelicidins, which provide innate immune protection during early life, are expressed in the pouch lining, skin and milk of devil dams. These peptides are active against pathogens identified in the pouch microbiome. Of the six characterised cathelicidins, Saha-CATH5 and 6 have broad-spectrum antibacterial activity and are capable of killing problematic human pathogens including methicillin-resistant S. aureus and vancomycin-resistant E. faecalis, while Saha-CATH3 is active against fungi. Saha-CATH5 and 6 were toxic to human A549 cells at 500 μg/mL, which is over seven times the concentration required to kill pathogens. The remaining devil cathelicidins were not active against tested bacterial or fungal strains, but are widely expressed throughout the body, such as in immune tissues, in digestive, respiratory and reproductive tracts, and in the milk and pouch, which indicates that they are likely also important components of the devil immune system. Our results suggest cathelicidins play a role in protecting naive young during pouch life by passive immune transfer in the milk and may modulate pouch microbe populations to reduce potential pathogens.

Characterisation of the immune compounds in koala milk using a combined transcriptomic and proteomic approach

Production of milk is a key characteristic of mammals, but the features of lactation vary greatly between monotreme, marsupial and eutherian mammals. Marsupials have a short gestation followed by a long lactation period, and milk constituents vary greatly across lactation. Marsupials are born immunologically naïve and rely on their mother’s milk for immunological protection. Koalas (Phascolarctos cinereus) are an iconic Australian species that are increasingly threatened by disease. Here we use a mammary transcriptome, two milk proteomes and the koala genome to comprehensively characterise the protein components of koala milk across lactation, with a focus on immune constituents. The most abundant proteins were well-characterised milk proteins, including β-lactoglobulin and lactotransferrin. In the mammary transcriptome, 851 immune transcripts were expressed, including immunoglobulins and complement components. We identified many abundant antimicrobial peptides, as well as novel proteins with potential antimicrobial roles. We discovered that marsupial VELP is an ortholog of eutherian Glycam1, and likely has an antimicrobial function in milk. We also identified highly-abundant koala endogenous-retrovirus sequences, identifying a potential transmission route from mother to young. Characterising the immune components of milk is key to understanding protection of marsupial young, and the novel immune compounds identified may have applications in clinical research.

Transcriptome Analysis of the Innate Immunity-Related Complement System in Spleen Tissue of Ctenopharyngodon idella Infected with Aeromonas hydrophila

The grass carp (Ctenopharyngodon idella) is an important commercial farmed herbivorous fish species in China, but is susceptible to Aeromonas hydrophila infections. In the present study, we performed de novo RNA-Seq sequencing of spleen tissue from specimens of a disease-resistant family, which were given intra-peritoneal injections containing PBS with or without a dose of A. hydrophila. The fish were sampled from the control group at 0 h, and from the experimental group at 4, 8, 12, 24, 48 and 72 h. 122.18 million clean reads were obtained from the normalized cDNA libraries; these were assembled into 425,260 contigs and then 191,795 transcripts. Of those, 52,668 transcripts were annotated with the NCBI Nr database, and 41,347 of the annotated transcripts were assigned into 90 functional groups. 20,569 unigenes were classified into six main categories, including 38 secondary KEGG pathways. 2,992 unigenes were used in the analysis of differentially expressed genes (DEGs). 89 of the putative DEGs were related to the immune system and 41 of them were involved in the complement and coagulation cascades pathway. This study provides insights into the complement and complement-related pathways involved in innate immunity, through expression profile analysis of the genomic resources in C. idella. We conclude that complement and complement-related genes play important roles during defense against A. hydrophila infection. The immune response is activated at 4 h after the bacterial injections, indicating that the complement pathways are activated at the early stage of bacterial infection. The study has improved our understanding of the immune response mechanisms in C. idella to bacterial pathogens.