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

Genome-wide characterization and comparison of endogenous retroviruses among 3 duck reference genomes

Endogenous retroviruses (ERV) are viral genomes integrated into the host genome and can be stably inherited. Although ERV sequences have been reported in some avian species' genome, the duck endogenous retroviruses (DERV) genome has yet to be quantified. This study aimed to identify ERV sequences and characterize genes near ERVs in the duck genome by utilizing LTRhavest and LTRdigest tools to forecast the duck genome and analyze the distribution of ERV copies. The results revealed 1,607, 2,031, and 1,908 full-length ERV copies in the Pekin duck (ZJU1.0), Mallard (CAU_wild_1.0), and Shaoxing duck (CAU_laying_1.0) genomes, respectively, with average lengths of 7,046, 7,027, and 6,945 bp. ERVs are mainly distributed on the 1, 2, and sex chromosomes. Phylogenetic analysis demonstrated the presence of Betaretrovirus in 3 duck genomes, whereas Alpharetrovirus was exclusively identified in the Shaoxing duck genome. Through screening, 596, 315, and 343 genes adjacent to ERV were identified in 3 duck genomes, respectively, and their functions of ERV neighboring genes were predicted. Functional enrichment analysis of ERV-adjacent genes revealed enrichment for Focal adhesion, Calcium signaling pathway, and Adherens junction in 3 duck genomes. The overlapped genes were highly expressed in 8 tissues (brain, fat, heart, kidney, liver, lung, skin, and spleen) of 8-wk-old Mallard, revealing their important expression in different tissues. Our study provides a new perspective for understanding the quantity and function of DERVs, and may also provide important clues for regulating nearby genes and affecting the traits of organisms.

Metagenomic Analysis of Viromes of Aedes Mosquitoes across India

Metagenomic analysis of Aedes aegypti and Ae. albopictus mosquitoes from diverse geographical regions of India revealed the presence of several insect viruses of human interest. Most abundant reads found in Ae. aegypti mosquitoes were of Phasi Charoen-like virus (PCLV), Choristoneura fumiferana granulovirus (CfGV), Cell fusing agent virus (CFAV), and Wenzhou sobemo-like virus 4 (WSLV4), whereas WSLV4 and CfGV constituted the highest percentage of reads in Ae. albopictus viromes. Other reads that were of low percentage included Hubei mosquito virus 2 (HMV2), Porcine astrovirus 4 (PAstV4), and Wild Boar astrovirus (WBAstV). PCLV and CFAV, which were found to be abundant in Ae. aegypti viromes were absent in Ae. albopictus viromes. Among the viromes analyzed, Ae. aegypti sampled from Pune showed the highest percentage (79.82%) of viral reads, while Ae. aegypti mosquitoes sampled from Dibrugarh showed the lowest percentage (3.47%). Shamonda orthobunyavirus (SHAV), African swine fever virus (ASFV), Aroa virus (AROAV), and Ilheus virus (ILHV), having the potential to infect vertebrates, including humans, were also detected in both mosquito species, albeit with low read numbers. Reads of gemykibivirus, avian retrovirus, bacteriophages, herpesviruses, and viruses infecting protozoans, algae, etc., were also detected in the mosquitoes. A high percentage of reads in the Ae. albopictus mosquito samples belonged to unclassified viruses and warrant further investigation. The data generated in the present work may not only lead to studies to explain the influence of these viruses on the replication and transmission of viruses of clinical importance but also to find applications as biocontrol agents against pathogenic viruses.

Comparison of Multiple Strategies for Precision Transgene Knock-In in Gallus gallus Genome via Microhomology-Mediated End Joining

Precision exogenous gene knock-in is an attractive field for transgenic Gallus gallus (chicken) generation. In this article, we constructed multiple Precise Integration into Target Chromosome (PITCh) plasmid systems mediated by microhomology-mediated end-joining (MMEJ) for large-fragment integration in DF-1 cells and further assess the possibility of GAPDH (glyceraldehyde-3-phosphate dehydrogenase) as a genomic safe harbor for chickens. We designed three targeted sgRNAs for the all-in-one plasmid at the 3'UTR of GAPDH near the stop codon. The donor-plasmid-carrying microhomology arms correspond to sgRNA and EGFP fragments in the forward and reverse directions. MMEJ-mediated EGFP insertion can be efficiently expressed in DF-1 cells. Moreover, the differences between the forward and reverse fragments indicated that promoter interference does affect the transfection efficiency of plasmids and cell proliferation. The comparison of the 20 bp and 40 bp microhomology arms declared that the short one has higher knock-in efficiency. Even though all three different transgene insertion sites in GAPDH could be used to integrate the foreign gene, we noticed that the G2-20R-EGFP cell reduced the expression of GAPDH, and the G3-20R-EGFP cell exhibited significant growth retardation. Taken together, G1, located at the 3'UTR of GAPDH on the outer side of the last base of the terminator, can be a candidate genomic safe harbor (GSH) loci for the chicken genome. In addition, deleted-in-azoospermia-like (DAZL) and actin beta (ACTB) site-specific gene knock-in indicated that MMEJ has broad applicability and high-precision knock-in efficiency for genetically engineered chickens.

Detection of long terminal repeat loci derived from endogenous retrovirus in junglefowl using whole-genome sequencing

Endogenous retroviruses (ERVs) are genetic elements present in the genome that retain traces of past viral infections. Characterization of ERVs can provide crucial insights into avian evolution. This study aimed to identify novel long terminal repeat (LTR) loci derived from ERVs (ERV-LTRs) absent in the reference genome using whole-genome sequencing data of red junglefowl, gray junglefowl, Ceylon junglefowl, and green junglefowl. In total, 835 ERV-LTR loci were identified across the four Gallus species. The numbers of ERV-LTRs loci detected in red junglefowl and its subspecies gray junglefowl, Ceylon junglefowl, and green junglefowl were 362, 216, 193, and 128, respectively. The phylogenetic tree was congruent with previously reported trees, suggesting the potential for inferring relationships among past junglefowl populations from the identified ERV-LTR loci. Of the detected loci, 306 ERV-LTRs were identified near or within the genes, and some were associated with cell adhesion. The detected ERV-LTR sequences were classified as endogenous avian retrovirus family, avian leukosis virus subgroup E, Ovex-1, and murine leukemia virus-related ERVs. In addition, the sequence of the EAV family was divided into four patterns by combining the U3, R, and U5 regions. These findings contribute to a more comprehensive understanding of the characteristics of junglefowl ERVs.

The association between microbial community and ileal gene expression on intestinal wall thickness alterations in chickens

The dynamic development of the animal intestine with a concurrent succession of microbiota and changes in microbial community and metabolite spectrum can exert far-reaching effects on host physiology. However, the precise mechanism of mutual response between microbiota and the gut is yet to be fully elucidated. Broilers with varying developmental degrees of intestinal wall thickness were selected, and they were divided into the thick group (H type) and the thin group (B type), using multiomics data integration analysis to reveal the fundamental regulatory mechanisms of gut–microbiota interplay. Our data showed, in broilers with similar body weight, the intestinal morphological parameters were improved in H type and the diversity of microbial communities is distinguishable from each other. The beneficial bacteria such as Bifidobacterium breve was increased whereas avian endogenous retrovirus EAV-HP was decreased in the H type compared with the B type. Furthermore, microbial metabolic potentials were more active, especially the biosynthesis of folate was improved in the H type. Similarly, the consolidation of absorption, immunity, metabolism, and development was noticed in the thick group. Correlation analysis indicated that the expression levels of material transport and immunomodulatory-related genes were positively correlated with the relative abundance of several probiotics such as B. breve, Lactobacillus saerimneri, and Butyricicoccus pullicaecorum. Our findings suggest that the chickens with well-developed ileal thickness own exclusive microbial composition and metabolic potential, which is closely related to small intestinal morphogenesis and homeostasis.

Whole-Genome Resequencing of Red Junglefowl and Indigenous Village Chicken Reveal New Insights on the Genome Dynamics of the Species

The red junglefowl Gallus gallus is the main progenitor of domestic chicken, the commonest livestock species, outnumbering humans by an approximate ratio of six to one. The genetic control for production traits have been well studied in commercial chicken, but the selection pressures underlying unique adaptation and production to local environments remain largely unknown in indigenous village chicken. Likewise, the genome regions under positive selection in the wild red junglefowl remain untapped. Here, using the pool heterozygosity approach, we analyzed indigenous village chicken populations from Ethiopia, Saudi Arabia, and Sri Lanka, alongside six red junglefowl, for signatures of positive selection across the autosomes. Two red junglefowl candidate selected regions were shared with all domestic chicken populations. Four candidates sweep regions, unique to and shared among all indigenous domestic chicken, were detected. Only one region includes annotated genes (TSHR and GTF2A1). Candidate regions that were unique to each domestic chicken population with functions relating to adaptation to temperature gradient, production, reproduction and immunity were identified. Our results provide new insights on the consequence of the selection pressures that followed domestication on the genome landscape of the domestic village chicken.

Enhancing Targeted Genomic DNA Editing in Chicken Cells Using the CRISPR/Cas9 System

The CRISPR/Cas9 system has enabled highly efficient genome targeted editing for various organisms. However, few studies have focused on CRISPR/Cas9 nuclease-mediated chicken genome editing compared with mammalian genomes. The current study combined CRISPR with yeast Rad52 (yRad52) to enhance targeted genomic DNA editing in chicken DF-1 cells. The efficiency of CRISPR/Cas9 nuclease-induced targeted mutations in the chicken genome was increased to 41.9% via the enrichment of the dual-reporter surrogate system. In addition, the combined effect of CRISPR nuclease and yRad52 dramatically increased the efficiency of the targeted substitution in the myostatin gene using 50-mer oligodeoxynucleotides (ssODN) as the donor DNA, resulting in a 36.7% editing efficiency after puromycin selection. Furthermore, based on the effect of yRad52, the frequency of exogenous gene integration in the chicken genome was more than 3-fold higher than that without yRad52. Collectively, these results suggest that ssODN is an ideal donor DNA for targeted substitution and that CRISPR/Cas9 combined with yRad52 significantly enhances chicken genome editing. These findings could be extensively applied in other organisms.

Comprehensive profiling of retroviral integration sites using target enrichment methods from historical koala samples without an assembled reference genome

Background. Retroviral integration into the host germline results in permanent viral colonization of vertebrate genomes. The koala retrovirus (KoRV) is currently invading the germline of the koala (Phascolarctos cinereus) and provides a unique opportunity for studying retroviral endogenization. Previous analysis of KoRV integration patterns in modern koalas demonstrate that they share integration sites primarily if they are related, indicating that the process is currently driven by vertical transmission rather than infection. However, due to methodological challenges, KoRV integrations have not been comprehensively characterized.

Results. To overcome these challenges, we applied and compared three target enrichment techniques coupled with next generation sequencing (NGS) and a newly customized sequence-clustering based computational pipeline to determine the integration sites for 10 museum Queensland and New South Wales (NSW) koala samples collected between the 1870s and late 1980s. A secondary aim of this study sought to identify common integration sites across modern and historical specimens by comparing our dataset to previously published studies. Several million sequences were processed, and the KoRV integration sites in each koala were characterized.

Conclusions. Although the three enrichment methods each exhibited bias in integration site retrieval, a combination of two methods, Primer Extension Capture and hybridization capture is recommended for future studies on historical samples. Moreover, identification of integration sites shows that the proportion of integration sites shared between any two koalas is quite small.