Endogenous avian leukosis viral loci in the Red Jungle Fowl genome assembly

Short communication: diversity of endogenous avian leukosis virus subgroup E elements in 11 chicken breeds

Avian leukosis virus subgroup E (ALVE) as a kind of endogenous retroviruses extensively exists in chicken genome. The insertion of ALVE has some effects on chicken production traits and appearance. Most of the work on ALVEs has been done with commercial breeds. We present here an investigation of ALVE elements in seven Chinese domestic breeds and four standard breeds. Firstly, we established an ALVE insertion site dataset by using the obsERVer pipeline to identify ALVEs from whole-genome sequence data of eleven chicken breeds, seven Chinese domestic breeds, including Beijing You (BY), Dongxiang (DX), Luxi Game (LX), Shouguang (SG), Silkie (SK), Tibetan (TB) and Wenchang (WC), four standard breeds, including White Leghorn (WL), White Plymouth Rock (WR), Cornish (CS), and Rhode Island Red (RIR). A total of 37 ALVE insertion sites were identified and 23 of them were novel. Most of these insertion sites were distributed in intergenic regions and introns. We then used locus-specific PCR to validate the insertion sites in an expanded population with 18~60 individuals in each breed. The results showed that all predicted integration sites in 11 breeds were verified by PCR. Some ALVE insertion sites were breeds specific, and 16 out of 23 novel ALVEs were found in only one Chinese domestic chicken breed. We randomly selected three ALVE insertions including ALVE_CAU005, ALVE_ros127, and ALVE_ros276, and obtained their insertion sequences by long-range PCR and Sanger sequencing. The insertion sequences were all 7525 bp, which were full-length ALVE insertion and all of them were highly homologous to ALVE1 with similarity of 99%. Our study identified the distribution of ALVE in 11 chicken breeds, which expands the current research on ALVE in Chinese domestic breeds.

Diversity of Avian leukosis virus subgroup J in local chickens, Jiangxi, China

Avian leukosis caused by avian leukosis virus (ALV) is one of the most severe diseases endangering the poultry industry. When the eradication measures performed in commercial broilers and layers have achieved excellent results, ALV in some local chickens has gradually attracted attention. Since late 2018, following the re-outbreak of ALV-J in white feather broilers in China, AL-like symptoms also suddenly broke out in some local flocks, leading to great economic losses. In this study, a systematic epidemiological survey was carried out in eight local chicken flocks in Jiangxi Province, China, and 71 strains were finally isolated from 560 samples, with the env sequences of them being successfully sequenced. All of those new isolates belong to subgroup J but they have different molecular features and were very different from the strains that emerged in white feature broilers recently, with some strains being highly consistent with those previously isolated from commercial broilers, layers and other flocks or even isolated from USA and Russian, suggesting these local chickens have been acted as reservoirs to accumulate various ALV-J strains for a long time. More seriously, phylogenetic analysis shows that there were also many novel strains emerging and in a separate evolutionary branch, indicating several new mutated ALVs are being bred in local chickens. Besides, ALV-J strains isolated in this study can be further divided into ten groups, while there were more or fewer groups in different chickens, revealing that ALV may cross propagate in those flocks. The above analyses explain the complex background and future evolution trend of ALV-J in Chinese local chickens, providing theoretical support for the establishment of corresponding prevention and control measures.

Endogenous Retroviruses Drive Resistance and Promotion of Exogenous Retroviral Homologs

Endogenous retroviruses (ERVs) serve as markers of ancient viral infections and provide invaluable insight into host and viral evolution. ERVs have been exapted to assist in performing basic biological functions, including placentation, immune modulation, and oncogenesis. A subset of ERVs share high nucleotide similarity to circulating horizontally transmitted exogenous retrovirus (XRV) progenitors. In these cases, ERV-XRV interactions have been documented and include (a) recombination to result in ERV-XRV chimeras, (b) ERV induction of immune self-tolerance to XRV antigens, (c) ERV antigen interference with XRV receptor binding, and (d) interactions resulting in both enhancement and restriction of XRV infections. Whereas the mechanisms governing recombination and immune self-tolerance have been partially determined, enhancement and restriction of XRV infection are virus specific and only partially understood. This review summarizes interactions between six unique ERV-XRV pairs, highlighting important ERV biological functions and potential evolutionary histories in vertebrate hosts.

Identification and characterisation of endogenous Avian Leukosis Virus subgroup E (ALVE) insertions in chicken whole genome sequencing data

Background

Endogenous retroviruses (ERVs) are the remnants of retroviral infections which can elicit prolonged genomic and immunological stress on their host organism. In chickens, endogenous Avian Leukosis Virus subgroup E (ALVE) expression has been associated with reductions in muscle growth rate and egg production, as well as providing the potential for novel recombinant viruses. However, ALVEs can remain in commercial stock due to their incomplete identification and association with desirable traits, such as ALVE21 and slow feathering. The availability of whole genome sequencing (WGS) data facilitates high-throughput identification and characterisation of these retroviral remnants.

Results

We have developed obsERVer, a new bioinformatic ERV identification pipeline which can identify ALVEs in WGS data without further sequencing. With this pipeline, 20 ALVEs were identified across eight elite layer lines from Hy-Line International, including four novel integrations and characterisation of a fast feathered phenotypic revertant that still contained ALVE21. These bioinformatically detected sites were subsequently validated using new high-throughput KASP assays, which showed that obsERVer was highly precise and exhibited a 0% false discovery rate. A further fifty-seven diverse chicken WGS datasets were analysed for their ALVE content, identifying a total of 322 integration sites, over 80% of which were novel. Like exogenous ALV, ALVEs show site preference for proximity to protein-coding genes, but also exhibit signs of selection against deleterious integrations within genes.

Conclusions

obsERVer is a highly precise and broadly applicable pipeline for identifying retroviral integrations in WGS data. ALVE identification in commercial layers has aided development of high-throughput diagnostic assays which will aid ALVE management, with the aim to eventually eradicate ALVEs from high performance lines. Analysis of non-commercial chicken datasets with obsERVer has revealed broad ALVE diversity and facilitates the study of the biological effects of these ERVs in wild and domesticated populations.

Endogenous avian leukosis virus subgroup E elements of the chicken reference genome

The chicken reference genome contains 2 endogenous avian leukosis virus subgroup E (ALVE) insertions, but gaps and unresolved repetitive sequences in previous assemblies have hindered their precise characterization. Detailed analysis of the most recent reference genome (GRCg6a) now shows both ALVEs within contiguous chromosome assemblies for the first time. ALVE6 (ALVE-JFevA) and ALVE-JFevB are both located on chromosome 1, with ALVE6 close to the p-arm telomere. ALVE-JFevB is a structurally intact element containing the ALVE gag, pol, and env genes and is capable of forming replication competent viruses. In contrast, ALVE6 contains a 3,352 bp 5′ truncation and lacks the entire 5′ long terminal repeat and gag gene. Despite this, ALVE6 remains able to produce intact envelope protein, likely due to a mutation in the recognition site for a known inhibitory miRNA (miR-155). Whole genome resequencing data sets from layers, broilers, and 3 independent sources of wild-caught red junglefowl were surveyed for the presence of each of these reference genome ALVEs. ALVE-JFevB was found in no other chicken or red junglefowl genomes, whereas ALVE6 was identified in some layers, broilers, and native breeds but not within any other red junglefowl genome. Improved assembly contiguity has facilitated better characterization of the 2 ALVEs of the chicken reference genome. However, both the limited ALVE content and unique presence of ALVE-JFevB suggests that the reference individual is unrepresentative of ancestral Gallus gallus ALVE diversity.

GGERV20, a recently integrated, segregating endogenous retrovirus in Gallus gallus

Endogenous retroviruses (ERVs) are widespread in vertebrate genomes. The recent availability of whole eukaryotic genomes has enabled their characterization in many organisms, including Gallus gallus (red jungle fowl), the progenitor of the domesticated chicken. Our bioinformatics analysis of a G. gallus ERV previously designated GGERV20 identified 35 proviruses with complete long terminal repeats (LTRs) and gag-pol open reading frames (ORFs) in the Genome Reference Consortium Chicken Build 6a, of which 8 showed potential for translation of functional retroviral polyproteins, including the integrase and reverse transcriptase enzymes. No elements were discovered with an env gene. Fifteen loci had LTR sequences with 100 % identity, indicative of recent integration. Chicken embryo fibroblast RNA-seq datasets showed reads representing the entire length of the GGERV20 provirus, supporting their potential for expressing viral proteins. To investigate the possibility that GGERV20 elements may not be fixed in the genome, we assessed the integration status of five loci in a meat-type chicken. PCRs targeting a GGERV20 locus on G. gallus chromosome one (GGERV20_1-1) reproducibly amplified both LTRs and the preintegration state, indicating that the bird from which the DNA was sampled was hemizygous at this locus. The four other loci examined only produced the preintegration state amplicons. These results reveal that GGERV20 is not fixed in the G. gallus population, and taken together with the lack of mutations seen in several provirus LTRs and their transcriptional activity, suggest that GGERV20 retroviruses have recently been and continue to be active in the chicken genome.

Genomic and transcriptional analyses of novel parvoviruses identified from dead peafowl

To identify potential pathogens responsible for a disease outbreak of cultured peafowls in China in 2013, metagenomic sequencing was conducted. The genomes of two closely related parvoviruses, namely peafowl parvovirus 1 (PePV1) and PePV2, were identified with size of 4428 bp and 4348 bp, respectively. Phylogenetic analysis revealed that both viruses are novel parvoviruses, belonging to the proposed genus Chapparvovirus of Parvoviridae. The transcriptional profile of PePV1 was analyzed by transfecting a nearly complete PePV1 genome into HEK-293T cells. Results revealed that PePV1 employs one promoter and two polyadenylation sites to start and terminate its transcriptions, with one donor site and two acceptor sites for pre-mRNA splicing. PePV1 DNA and structural protein were detected in several tissues of a dead peafowl, which appeared to have suffered enteritis, pneumonia and viremia. These results provide novel information of chapparvoviruses, and call for attention to the potential pathogens.

Domestic chickens activate a piRNA defense against avian leukosis virus

PIWI-interacting RNAs (piRNAs) protect the germ line by targeting transposable elements (TEs) through the base-pair complementarity. We do not know how piRNAs co-evolve with TEs in chickens. Here we reported that all active TEs in the chicken germ line are targeted by piRNAs, and as TEs lose their activity, the corresponding piRNAs erode away. We observed de novo piRNA birth as host responds to a recent retroviral invasion. Avian leukosis virus (ALV) has endogenized prior to chicken domestication, remains infectious, and threatens poultry industry. Domestic fowl produce piRNAs targeting ALV from one ALV provirus that was known to render its host ALV resistant. This proviral locus does not produce piRNAs in undomesticated wild chickens. Our findings uncover rapid piRNA evolution reflecting contemporary TE activity, identify a new piRNA acquisition modality by activating a pre-existing genomic locus, and extend piRNA defense roles to include the period when endogenous retroviruses are still infectious.

DOI:http://dx.doi.org/10.7554/eLife.24695.001

Viruses called retroviruses can infect animal cells and merge their genetic information with those of the animal causing damage to the animal’s genetic blueprints. Once retroviruses are integrated into a cell they can sometimes get passed down through the generations over the centuries. Almost half of the human genetic code, for example, is made from ancient retroviruses and other foreign sequences. Over time many of these ancient viruses lost the ability to infect other cells and became trapped within cells but they can still jump out and damage the animal’s genetic code under certain circumstances. These trapped foreign sequences are called transposable elements.

Animal cells produce molecules called piRNAs to shut down transposable elements. Most piRNAs are produced from genetic information that originally came from integrated retroviruses and that has been hijacked to defend the cell, a similar strategy as Crisper system in bacteria. Domestic chickens produce piRNAs against a virus called avian leukosis virus (or ALV for short) – which commonly infects domestic fowl. The virus also infected the wild ancestors of chickens, known as red jungle fowl, but these birds do not produce piRNAs. This provides an ideal setting to study the evolution of piRNAs in an animal that is not too distantly related to humans (chickens and humans both have backbones, and are therefore both warm-blooded vertebrates).

Sun et al. examined cells from the testicles of domestic chickens and red jungle fowl as an example of the role of piRNAs in protecting genetic information in vertebrates. The investigation revealed that piRNAs against all previously trapped viruses in the chicken’s genetic code are produced in chickens to stop them from causing more damage. Sun et al. also observed the creation of piRNAs in chickens in response to ALV that had not yet become trapped in the chicken’s genetic code. Importantly, the piRNAs could control these retroviruses while they were still infectious.

The experiments also revealed that piRNAs against ALV are produced from a single copy of ALV that is found in both domestic and wild chickens. The results showed that cells can produce new piRNAs using these pre-existing viral copies within their own genetics. This illustrates that production of piRNA from existing genetic material can be activated in response to certain cues.

Further work will seek to discover how existing genetic information becomes a source of piRNAs. In the United States, 8 billion domestic chickens are consumed each year, and a better understanding of how these birds defend themselves against viral infections could increase the productivity of the poultry industry around the world. Moreover, because other viruses trapped in the chicken’s genetic code are related to similar viruses in humans, future discoveries made in this area could help to guide research that will benefit human health as well.

DOI:http://dx.doi.org/10.7554/eLife.24695.002

Discovery of an expanded set of avian leukosis subgroup E proviruses in chickens using Vermillion, a novel sequence capture and analysis pipeline [corrected]

Transposable elements (TEs), such as endogenous retroviruses (ERVs), are common in the genomes of vertebrates. ERVs result from retroviral infections of germ-line cells, and once integrated into host DNA they become part of the host's heritable genetic material. ERVs have been ascribed positive effects on host physiology such as the generation of novel, adaptive genetic variation and resistance to infection, as well as negative effects as agents of tumorigenesis and disease. The avian leukosis virus subgroup E family (ALVE) of endogenous viruses of chickens has been used as a model system for studying the effects of ERVs on host physiology, and approximately 30 distinct ALVE proviruses have been described in the Gallus gallus genome. In this report we describe the development of a software tool, which we call Vermillion, and the use of this tool in combination with targeted next-generation sequencing (NGS) to increase the number of known proviruses belonging to the ALVE family of ERVs in the chicken genome by 4-fold, including expanding the number of known ALVE elements on chromosome 1 (Gga1) from the current 9 to a total of 40. Although we focused on the discovery of ALVE elements in chickens, with appropriate selection of target sequences Vermillion can be used to develop profiles of other families of ERVs and TEs in chickens as well as in species other than the chicken.

Genomic Landscape of Long Terminal Repeat Retrotransposons (LTR-RTs) and Solo LTRs as Shaped by Ectopic Recombination in Chicken and Zebra Finch

Transposable elements (TEs) are nearly ubiquitous among eukaryotic genomes, but TE contents vary dramatically among phylogenetic lineages. Several mechanisms have been proposed as drivers of TE dynamics in genomes, including the fixation/loss of a particular TE insertion by selection or drift as well as structural changes in the genome due to mutation (e.g., recombination). In particular, recombination can have a significant and directional effect on the genomic TE landscape. For example, ectopic recombination removes internal regions of long terminal repeat retrotransposons (LTR-RTs) as well as one long terminal repeat (LTR), resulting in a solo LTR. In this study, we focus on the intra-species dynamics of LTR-RTs and solo LTRs in bird genomes. The distribution of LTR-RTs and solo LTRs in birds is intriguing because avian recombination rates vary widely within a given genome. We used published linkage maps and whole genome assemblies to study the relationship between recombination rates and LTR-removal events in the chicken and zebra finch. We hypothesized that regions with low recombination rates would harbor more full-length LTR-RTs (and fewer solo LTRs) than regions with high recombination rates. We tested this hypothesis by comparing the ratio of full-length LTR-RTs and solo LTRs across chromosomes, across non-overlapping megabase windows, and across physical features (i.e., centromeres and telomeres). The chicken data statistically supported the hypothesis that recombination rates are inversely correlated with the ratio of full-length to solo LTRs at both the chromosome level and in 1-Mb non-overlapping windows. We also found that the ratio of full-length to solo LTRs near chicken telomeres was significantly lower than those ratios near centromeres. Our results suggest a potential role of ectopic recombination in shaping the chicken LTR-RT genomic landscape.