Characterization of major histocompatibility complex class I and class II genes from the Tasmanian devil (<Emphasis Type="Italic">Sarcophilus harrisii</Emphasis>)

The Tasmanian devil (Sarcophilus harrisii) is currently threatened by an emerging wildlife disease, devil facial tumour disease. The disease is decreasing devil numbers dramatically and may lead to the extinction of the species. At present, nothing is known about the immune genes or basic immunology of the devil. In this study, we report the construction of the first genetic library for the Tasmanian devil, a spleen cDNA library, and the isolation of full-length MHC Class I and Class II genes. We describe six unique Class II beta chain sequences from at least three loci, which belong to the marsupial Class II DA gene family. We have isolated 13 unique devil Class I sequences, representing at least seven Class I loci, two of which are most likely non-classical genes. The MHC Class I sequences from the devil have little heterogeneity, indicating recent divergence. The MHC genes described here are most likely involved in antigen presentation and are an important first step for studying MHC diversity and immune response in the devil.     

Genomics for conservation: a case study of behavioral genes in the Tasmanian devil

The increased availability of genomic resources for many species has expanded perspectives on problems in conservation by helping to design management strategies for threatened species. Tasmanian devils (Sarcophilus harrisii) are an iconic and endangered marsupial with an intensively managed breeding program aimed at preventing extinction in the wild caused by devil facial tumour disease. Between 2015 and 2017, 85 devils from this program were released to three sites in Tasmania to support wild populations. Of these, 26 were known to have been killed by vehicles shortly after release. A previous analysis indicated that increased generations in captivity was a positive predictor of vehicle strike, with possible behavioural change hypothesised. Here we use 39 resequenced devil genomes to characterise diversity at 35 behaviour-associated genes, which contained 826 single nucleotide polymorphisms (24 were non-synonymous). We tested for a predictor of survival by examining three genes (AVPR1B, OXT and SLC6A4) in 62 released devils with known fates (survived, N?=?39; died, N?=?23), and genome-wide associations via reduced-representation sequencing (1727 single nucleotide polymorphisms [SNPs]), in 55 devils with known fates (survived, N?=?38; died, N?=?17). Overall, there was little evidence of an association between genetic profile and probability of being struck by a vehicle. Despite previous evidence of low genetic diversity in devils, the 35 behaviour-associated genes contained variation that may influence their functions. Our dataset can be used for future research into devil behavioural ecology, and adds to the increasing body of research applying genomics to conservation problems.

     

Allorecognition in the Tasmanian devil (Sarcophilus harrisii), an endangered marsupial species with limited genetic diversity

Tasmanian devils (Sarcophilus harrisii) are on the verge of extinction due to a transmissible cancer, devil facial tumour disease (DFTD). This tumour is an allograft that is transmitted between individuals without immune recognition of the tumour cells. The mechanism to explain this lack of immune recognition and acceptance is not well understood. It has been hypothesized that lack of genetic diversity at the Major Histocompatibility Complex (MHC) allowed the tumour cells to grow in genetically similar hosts without evoking an immune response to alloantigens. We conducted mixed lymphocyte reactions and skin grafts to measure functional MHC diversity in the Tasmanian devil population. The limited MHC diversity was sufficient to produce measurable mixed lymphocyte reactions. There was a wide range of responses, from low or no reaction to relatively strong responses. The highest responses occurred when lymphocytes from devils from the east of Tasmania were mixed with lymphocytes from devils from the west of Tasmania. All of the five successful skin allografts were rejected within 14 days after surgery, even though little or no MHC I and II mismatches were found. Extensive T-cell infiltration characterised the immune rejection. We conclude that Tasmanian devils are capable of allogeneic rejection. Consequently, a lack of functional allorecognition mechanisms in the devil population does not explain the transmission of a contagious cancer.     

High MHC gene copy number maintains diversity despite homozygosity in a Critically Endangered single‐island endemic bird,but no evidence of MHC‐based mate choice

Small population sizes can, over time, put species at risk due to the loss of genetic variation and the deleterious effects of inbreeding. Losing diversity in the major histocompatibility complex (MHC) could be particularly harmful, given its key role in the immune system. Here, we assess MHC class I (MHC‐I) diversity and its effects on mate choice and survival in the Critically Endangered Raso lark Alauda razae, a species restricted to the 7 km² islet of Raso, Cape Verde, since ~1460, whose population size has dropped as low as 20 pairs. Exhaustively genotyping 122 individuals, we find no effect of MHC‐I genotype/diversity on mate choice or survival. However, we demonstrate that MHC‐I diversity has been maintained through extreme bottlenecks by retention of a high number of gene copies (at least 14), aided by cosegregation of multiple haplotypes comprising 2–8 linked MHC‐I loci. Within‐locus homozygosity is high, contributing to low population‐wide diversity. Conversely, each individual had comparably many alleles, 6–16 (average 11), and the large and divergent haplotypes occur at high frequency in the population, resulting in high within‐individual MHC‐I diversity. This functional immune gene diversity will be of critical importance for this highly threatened species’ adaptive potential.     

Finding needles in a genomic haystack: targeted capture identifies clear signatures of selection in a nonmodel plant species

Teasing apart neutral and adaptive genomic processes and identifying loci that are targets of selection can be difficult, particularly for nonmodel species that lack a reference genome. However, identifying such loci and the factors driving selection have the potential to greatly assist conservation and restoration practices, especially for the management of species in the face of contemporary and future climate change. Here, we focus on assessing adaptive genomic variation within a nonmodel plant species, the narrow‐leaf hopbush (Dodonaea viscosa ssp. angustissima), commonly used for restoration in Australia. We used a hybrid‐capture target enrichment approach to selectively sequence 970 genes across 17 populations along a latitudinal gradient from 30°S to 36°S. We analysed 8462 single‐nucleotide polymorphisms (SNPs) for F_ST outliers as well as associations with environmental variables. Using three different methods, we found 55 SNPs with significant correlations to temperature and water availability, and 38 SNPs to elevation. Genes containing SNPs identified as under environmental selection were diverse, including aquaporin and abscisic acid genes, as well as genes with ontologies relating to responses to environmental stressors such as water deprivation and salt stress. Redundancy analysis demonstrated that only a small proportion of the total genetic variance was explained by environmental variables. We demonstrate that selection has led to clines in allele frequencies in a number of functional genes, including those linked to leaf shape and stomatal variation, which have been previously observed to vary along the sampled environmental cline. Using our approach, gene regions subject to environmental selection can be readily identified for nonmodel organisms.     

MHC gene copy number variation in Tasmanian devils: implications for the spread of a contagious cancer

Tasmanian devils face extinction owing to the emergence of a contagious cancer. Devil facial tumour disease (DFTD) is a clonal cancer spread owing to a lack of major histocompatibility complex (MHC) barriers in Tasmanian devil populations. We present a comprehensive screen of MHC diversity in devils and identify 25 MHC types and 53 novel sequences, but conclude that overall levels of MHC diversity at the sequence level are low. The majority of MHC Class I variation can be explained by allelic copy number variation with two to seven sequence variants identified per individual. MHC sequences are divided into two distinct groups based on sequence similarity. DFTD cells and most devils have sequences from both groups. Twenty per cent of individuals have a restricted MHC repertoire and contain only group I or only group II sequences. Counterintuitively, we postulate that the immune system of individuals with a restricted MHC repertoire may recognize foreign MHC antigens on the surface of the DFTD cell. The implication of these results for management of DFTD and this endangered species are discussed.     

Polymorphic MHC loci in an asexual fish,the amazon molly (Poecilia formosa; Poeciliidae)

Genes of the major histocompatibility complex (MHC) encode molecules that control immune recognition and are highly polymorphic in most vertebrates. The remarkable polymorphisms at MHC loci may be maintained by selection from parasites, sexual selection, or both. If asexual species show equal (or higher) levels of polymorphisms at MHC loci as sexual ones, this would mean that sexual selection is not necessary to explain the high levels of diversity at MHC loci. In this study, we surveyed the MHC diversity of the asexual amazon molly (Poecilia formosa) and one of its sexual ancestors, the sailfin molly (P. latipinna), which lives in the same habitat. We found that the asexual molly has polymorphic MHC loci despite its clonal reproduction, yet not as polymorphic as the sexual species. Although the nucleotide diversity was similar between the asexual and sexual species, the sexual species exhibited a greater genotypic diversity compared to the asexual one from the same habitats. Within‐genome diversity was similar for MHC class I loci, but for class IIB, the sexual species had higher diversity compared to the asexual — despite the hybrid origins and higher levels of heterozygosity at microsatellite loci in the asexual species. The level of positive selection appears to be similar between the two species, which suggests that these polymorphisms are maintained by selection. Thus, our findings do not allow us to rule out the sexual selection hypothesis for the evolution of MHC diversity, and although the sexual fish has higher levels of MHC‐diversity compared to the asexual species, this may be due to differences in demography, parasites, or other factors, rather than sexual selection.     

Specific alleles at immune genes,rather than genome‐wide heterozygosity,are related to immunity and survival in the critically endangered Attwater's prairie‐chicken

The negative effects of inbreeding on fitness are serious concerns for populations of endangered species. Reduced fitness has been associated with lower genome‐wide heterozygosity and immune gene diversity in the wild; however, it is rare that both types of genetic measures are included in the same study. Thus, it is often unclear whether the variation in fitness is due to the general effects of inbreeding, immunity‐related genes or both. Here, we tested whether genome‐wide heterozygosity (20 990 SNPs) and diversity at nine immune genes were better predictors of two measures of fitness (immune response and survival) in the endangered Attwater's prairie‐chicken (Tympanuchus cupido attwateri). We found that postrelease survival of captive‐bred birds was related to alleles of the innate (Toll‐like receptors, TLRs) and adaptive (major histocompatibility complex, MHC) immune systems, but not to genome‐wide heterozygosity. Likewise, we found that the immune response at the time of release was related to TLR and MHC alleles, and not to genome‐wide heterozygosity. Overall, this study demonstrates that immune genes may serve as important genetic markers when monitoring fitness in inbred populations and that in some populations specific functional genes may be better predictors of fitness than genome‐wide heterozygosity.     

Low major histocompatibility complex diversity in the Tasmanian devil predates European settlement and may explain susceptibility to disease epidemics

The Tasmanian devil (Sarcophilus harrisii) is at risk of extinction owing to the emergence of a contagious cancer known as devil facial tumour disease (DFTD). The emergence and spread of DFTD has been linked to low genetic diversity in the major histocompatibility complex (MHC). We examined MHC diversity in historical and ancient devils to determine whether loss of diversity is recent or predates European settlement in Australia. Our results reveal no additional diversity in historical Tasmanian samples. Mainland devils had common modern variants plus six new variants that are highly similar to existing alleles. We conclude that low MHC diversity has been a feature of devil populations since at least the Mid-Holocene and could explain their tumultuous history of population crashes.     

Disease swamps molecular signatures of genetic-environmental associations to abiotic factors in Tasmanian devil (Sarcophilus harrisii) populations

Landscape genomics studies focus on identifying candidate genes under selection via spatial variation in abiotic environmental variables, but rarely by biotic factors (i.e., disease). The Tasmanian devil (Sarcophilus harrisii) is found only on the environmentally heterogeneous island of Tasmania and is threatened with extinction by a transmissible cancer, devil facial tumor disease (DFTD). Devils persist in regions of long-term infection despite epidemiological model predictions of species’ extinction, suggesting possible adaptation to DFTD. Here, we test the extent to which spatial variation and genetic diversity are associated with the abiotic environment (i.e., climatic variables, elevation, vegetation cover) and/or DFTD. We employ genetic-environment association analyses using 6886 SNPs from 3287 individuals sampled pre- and post-disease arrival across the devil's geographic range. Pre-disease, we find significant correlations of allele frequencies with environmental variables, including 365 unique loci linked to 71 genes, suggesting local adaptation to abiotic environment. The majority of candidate loci detected pre-DFTD are not detected post-DFTD arrival. Several post-DFTD candidate loci are associated with disease prevalence and were in linkage disequilibrium with genes involved in tumor suppression and immune response. Loss of apparent signal of abiotic local adaptation post-disease suggests swamping by strong selection resulting from the rapid onset of DFTD.     

Genomewide discovery of DNA polymorphisms in rice cultivars with contrasting drought and salinity stress response and their functional relevance

Next‐generation sequencing technologies provide opportunities to understand the genetic basis of phenotypic differences, such as abiotic stress response, even in the closely related cultivars via identification of large number of DNA polymorphisms. We performed whole‐genome resequencing of three rice cultivars with contrasting responses to drought and salinity stress (sensitive IR64, drought‐tolerant Nagina 22 and salinity‐tolerant Pokkali). More than 356 million 90‐bp paired‐end reads were generated, which provided about 85% coverage of the rice genome. Applying stringent parameters, we identified a total of 1 784 583 nonredundant single‐nucleotide polymorphisms (SNPs) and 154 275 InDels between reference (Nipponbare) and the three resequenced cultivars. We detected 401 683 and 662 509 SNPs between IR64 and Pokkali, and IR64 and N22 cultivars, respectively. The distribution of DNA polymorphisms was found to be uneven across and within the rice chromosomes. One‐fourth of the SNPs and InDels were detected in genic regions, and about 3.5% of the total SNPs resulted in nonsynonymous changes. Large‐effect SNPs and InDels, which affect the integrity of the encoded protein, were also identified. Further, we identified DNA polymorphisms present in the differentially expressed genes within the known quantitative trait loci. Among these, a total of 548 SNPs in 232 genes, located in the conserved functional domains, were identified. The data presented in this study provide functional markers and promising target genes for salinity and drought tolerance and present a valuable resource for high‐throughput genotyping and molecular breeding for abiotic stress traits in rice.     

Low MHC class II diversity in the Tasmanian devil (Sarcophilus harrisii)

The largest remaining carnivorous marsupial, the Tasmanian devil (Sarcophilus harrisii), is currently under threat of extinction due to a fatal contagious cancer-devil facial tumour disease. Low major histocompatibility complex (MHC) class I diversity is believed to have contributed to the transmission of the tumour allograft through devil populations. Here, we report low MHC class II variability in this species, with DA β chain genes (Saha-DAB1, 2 and 3) exhibiting very limited diversity and the sole α chain gene (Saha-DAA) monomorphic. Three, six and three alleles were found at Saha-DAB1, 2 and 3, respectively, with a predominant allele found at each locus. Heterozygosity at these three loci is low in the eastern population and modestly higher in northwestern individuals. The results are indicative of a selective sweep likely due to an infectious disease resulting in the fixation of selectively favoured alleles and depletion of genetic diversity at devil class II loci. Several attempts were made to isolate the other marsupial classical class II gene family, namely, DB, resulting in only one DBB pseudogene being found. These findings further support the view that this species has a compromised capacity to respond to pathogen evolution, emerging infectious diseases and environmental changes.     

Comparative genomics analysis in Prunoideae to identify biologically relevant polymorphisms

Prunus is an economically important genus with a wide range of physiological and biological variability. Using the peach genome as a reference, sequencing reads from four almond accessions and one sweet cherry cultivar were used for comparative analysis of these three Prunus species. Reference mapping enabled the identification of many biological relevant polymorphisms within the individuals. Examining the depth of the polymorphisms and the overall scaffold coverage, we identified many potentially interesting regions including hundreds of small scaffolds with no coverage from any individual. Non‐sense mutations account for about 70 000 of the 13 million identified single nucleotide polymorphisms (SNPs). Blast2GO analyses on these non‐sense SNPs revealed several interesting results. First, non‐sense SNPs were not evenly distributed across all gene ontology terms. Specifically, in comparison with peach, sweet cherry is found to have non‐sense SNPs in two 1‐aminocyclopropane‐1‐carboxylate synthase (ACS) genes and two 1‐aminocyclopropane‐1‐carboxylate oxidase (ACO) genes. These polymorphisms may be at the root of the nonclimacteric ripening of sweet cherry. A set of candidate genes associated with bitterness in almond were identified by comparing sweet and bitter almond sequences. To the best of our knowledge, this is the first report in plants of non‐sense SNP abundance in a genus being linked to specific GO terms.     

Genome‐wide variation in nucleotides and retrotransposons in alpine populations of Arabis alpina (Brassicaceae)

Advances in high‐throughput sequencing have promoted the collection of reference genomes and genome‐wide diversity. However, the assessment of genomic variation among populations has hitherto mainly been surveyed through single‐nucleotide polymorphisms (SNPs) and largely ignored the often major fraction of genomes represented by transposable elements (TEs). Despite accumulating evidence supporting the evolutionary significance of TEs, comprehensive surveys remain scarce. Here, we sequenced the full genomes of 304 individuals of Arabis alpina sampled from four nearby natural populations to genotype SNPs as well as polymorphic long terminal repeat retrotransposons (polymorphic TEs; i.e., presence/absence of TE insertions at specific loci). We identified 291,396 SNPs and 20,548 polymorphic TEs, comparing their contributions to genomic diversity and divergence across populations. Few SNPs were shared among populations and overall showed high population‐specific variation, whereas most polymorphic TEs segregated among populations. The genomic context of these two classes of variants further highlighted candidate adaptive loci having a putative impact on functional genes. In particular, 4.96% of the SNPs were identified as nonsynonymous or affecting start/stop codons. In contrast, 43% of the polymorphic TEs were present next to Arabis genes enriched in functional categories related to the regulation of reproduction and responses to biotic as well as abiotic stresses. This unprecedented data set, mapping variation gained from SNPs and complementary polymorphic TEs within and among populations, will serve as a rich resource for addressing microevolutionary processes shaping genome variation.     

Contrasting patterns of selection between MHC I and II across populations of Humboldt and Magellanic penguins

The evolutionary and adaptive potential of populations or species facing an emerging infectious disease depends on their genetic diversity in genes, such as the major histocompatibility complex (MHC). In birds, MHC class I deals predominantly with intracellular infections (e.g., viruses) and MHC class II with extracellular infections (e.g., bacteria). Therefore, patterns of MHC I and II diversity may differ between species and across populations of species depending on the relative effect of local and global environmental selective pressures, genetic drift, and gene flow. We hypothesize that high gene flow among populations of Humboldt and Magellanic penguins limits local adaptation in MHC I and MHC II, and signatures of selection differ between markers, locations, and species. We evaluated the MHC I and II diversity using 454 next‐generation sequencing of 100 Humboldt and 75 Magellanic penguins from seven different breeding colonies. Higher genetic diversity was observed in MHC I than MHC II for both species, explained by more than one MHC I loci identified. Large population sizes, high gene flow, and/or similar selection pressures maintain diversity but limit local adaptation in MHC I. A pattern of isolation by distance was observed for MHC II for Humboldt penguin suggesting local adaptation, mainly on the northernmost studied locality. Furthermore, trans‐species alleles were found due to a recent speciation for the genus or convergent evolution. High MHC I and MHC II gene diversity described is extremely advantageous for the long‐term survival of the species.     

Gene duplication and divergence produce divergent MHC genotypes without disassortative mating

Genes of the major histocompatibility complex (MHC) exhibit heterozygote advantage in immune defence, which in turn can select for MHC‐disassortative mate choice. However, many species lack this expected pattern of MHC‐disassortative mating. A possible explanation lies in evolutionary processes following gene duplication: if two duplicated MHC genes become functionally diverged from each other, offspring will inherit diverse multilocus genotypes even under random mating. We used locus‐specific primers for high‐throughput sequencing of two expressed MHC Class II B genes in Leach's storm‐petrels, Oceanodroma leucorhoa, and found that exon 2 alleles fall into two gene‐specific monophyletic clades. We tested for disassortative vs. random mating at these two functionally diverged Class II B genes, using multiple metrics and different subsets of exon 2 sequence data. With good statistical power, we consistently found random assortment of mates at MHC. Despite random mating, birds had MHC genotypes with functionally diverged alleles, averaging 13 amino acid differences in pairwise comparisons of exon 2 alleles within individuals. To test whether this high MHC diversity in individuals is driven by evolutionary divergence of the two duplicated genes, we built a phylogenetic permutation model. The model showed that genotypic diversity was strongly impacted by sequence divergence between the most common allele of each gene, with a smaller additional impact of monophyly of the two genes. Divergence of allele sequences between genes may have reduced the benefits of actively seeking MHC‐dissimilar mates, in which case the evolutionary history of duplicated genes is shaping the adaptive landscape of sexual selection.     

Molecular diversity and landscape genomics of the crop wild relative Triticum urartu across the Fertile Crescent

Modern plant breeding can benefit from the allelic variation that exists in natural populations of crop wild relatives that evolved under natural selection in varying pedoclimatic conditions. In this study, next‐generation sequencing was used to generate 1.3 million genome‐wide single nucleotide polymorphisms (SNPs) on ex situ collections of Triticum urartu L., the wild donor of the A^u subgenome of modern wheat. A set of 75 511 high‐quality SNPs were retained to describe 298 T. urartu accessions collected throughout the Fertile Crescent. Triticum urartu showed a complex pattern of genetic diversity, with two main genetic groups distributed sequentially from west to east. The incorporation of geographical information on sampling points showed that genetic diversity was correlated to the geographical distance (R² = 0.19) separating samples from Jordan and Lebanon, from Syria and southern Turkey, and from eastern Turkey, Iran and Iraq. The wild emmer genome was used to derive the physical positions of SNPs on the seven chromosomes of the A^u subgenome, allowing us to describe a relatively slow decay of linkage disequilibrium in the collection. Outlier loci were described on the basis of the geographic distribution of the T. urartu accessions, identifying a hotspot of directional selection on chromosome 4A. Bioclimatic variation was derived from grid data and related to allelic variation using a genome‐wide association approach, identifying several marker–environment associations (MEAs). Fifty‐seven MEAs were associated with altitude and temperature measures while 358 were associated with rainfall measures. The most significant MEAs and outlier loci were used to identify genomic loci with adaptive potential (some already reported in wheat), including dormancy and frost resistance loci. We advocate the application of genomics and landscape genomics on ex situ collections of crop wild relatives to efficiently identify promising alleles and genetic materials for incorporation into modern crop breeding.     

Chromosome‐level assembly of wild Bactrian camel genome reveals organization of immune gene loci

Camelids are characterized by their unique adaptive immune system that exhibits the generation of homodimeric heavy‐chain immunoglobulins, somatic hypermutation of T‐cell receptors, and low genetic diversity of major histocompatibility complex (MHC) genes. However, short‐read assemblies are typically highly fragmented in these gene loci owing to their repetitive and polymorphic nature. Here, we constructed a chromosome‐level assembly of wild Bactrian camel genome based on high‐coverage long‐read sequencing and chromatin interaction mapping. The assembly with a contig N50 of 5.37 Mb and a scaffold N50 of 76.03 Mb, represents the most contiguous camelid genome to date. The genomic organization of immunoglobulin heavy‐chain locus was similar between the wild Bactrian camel and alpaca, and genes encoding for conventional and heavy‐chain antibodies were intermixed. The organizations of two immunoglobulin light‐chain loci and four T cell receptor loci were also fully deciphered using the new assembly. Additionally, the complete classical MHC region was resolved into a single contig. The high‐quality assembly presented here provides an essential reference for future investigations examining the camelid immune system.     

Pathogen richness and abundance predict patterns of adaptive major histocompatibility complex variation in insular amphibians

The identification of the factors responsible for genetic variation and differentiation at adaptive loci can provide important insights into the evolutionary process and is crucial for the effective management of threatened species. We studied the impact of environmental viral richness and abundance on functional diversity and differentiation of the MHC class Ia locus in populations of the black‐spotted pond frog (Pelophylax nigromaculatus), an IUCN‐listed species, on 24 land‐bridge islands of the Zhoushan Archipelago and three nearby mainland sites. We found a high proportion of private MHC alleles in mainland and insular populations, corresponding to 32 distinct functional supertypes, and strong positive selection on MHC antigen‐binding sites in all populations. Viral pathogen diversity and abundance were reduced at island sites relative to the mainland, and islands housed distinctive viral communities. Standardized MHC diversity at island sites exceeded that found at neutral microsatellites, and the representation of key functional supertypes was positively correlated with the abundance of specific viruses in the environment (Frog virus 3 and Ambystoma tigrinum virus). These results indicate that pathogen‐driven diversifying selection can play an important role in maintaining functionally important MHC variation following island isolation, highlighting the importance of considering functionally important genetic variation and host–pathogen associations in conservation planning and management.