You will be shown a paragraph from a biomedical research paper, followed by a question about the paper. Please extract text from the paragraph to answer the question. Do not include any text in your repsonse other than the answer.

Paragraph: Gamma-and delta-CoVs mainly infect birds , but some can also infect mammals [4, 5] . Six human CoVs (HCoVs) are known to cause human diseases. HCoV-HKU1, HCoV-OC43, HCoV-229E, and HCoV-NL63 commonly cause mild respiratory illness or asymptomatic infection; however, severe acute respiratory syndrome coronavirus (SARS-CoV) and

All sampling procedures were performed by veterinarians, with approval from Animal Ethics Committee of the Wuhan Institute of Virology (WIVH5210201). The study was conducted in accordance with the Guide for the Care and Use of Wild Mammals in Research of the People's Republic of China. Bat fecal swab and pellet samples were collected from November 2004 to November 2014 in different seasons in Southern China, as described previously [16] . Viral RNA was extracted from 200 µL of fecal swab or pellet samples using the High Pure Viral RNA Kit (Roche Diagnostics GmbH, Mannheim, Germany) as per the manufacturer's instructions. RNA was eluted in 50 µL of elution buffer, aliquoted, and stored at -80 • C. One-step hemi-nested reverse-transcription (RT-) PCR (Invitrogen, San Diego, CA, USA) was employed to detect coronavirus, as previously described [17, 18] . To confirm the bat species of an individual sample, we PCR amplified the cytochrome b (Cytob) and/or NADH dehydrogenase subunit 1 (ND1) gene using DNA extracted from the feces or swabs [19, 20] . The gene sequences were assembled excluding the primer sequences. BLASTN was used to identify host species based on the most closely related sequences with the highest query coverage and a minimum identity of 95%. Full genomic sequences were determined by one-step PCR (Invitrogen, San Diego, CA, USA) amplification with degenerate primers (Table S1 ) designed on the basis of multiple alignments of available alpha-CoV sequences deposited in GenBank or amplified with SuperScript IV Reverse Transcriptase (Invitrogen) and Expand Long Template PCR System (Roche Diagnostics GmbH, Mannheim, Germany) with specific primers (primer sequences are available upon request). Sequences of the 5' and 3' genomic ends were obtained by 5' and 3' rapid amplification of cDNA ends (SMARTer Viruses 2019, 11, 379 3 of 19 RACE 5'/3' Kit; Clontech, Mountain View, CA, USA), respectively. PCR products were gel-purified and subjected directly to sequencing. PCR products over 5kb were subjected to deep sequencing using Hiseq2500 system. For some fragments, the PCR products were cloned into the pGEM-T Easy Vector (Promega, Madison, WI, USA) for sequencing. At least five independent clones were sequenced to obtain a consensus sequence. The Next Generation Sequencing (NGS) data were filtered and mapped to the reference sequence of BatCoV HKU10 (GenBank accession number NC_018871) using Geneious 7.1.8 [21] . Genomes were preliminarily assembled using DNAStar lasergene V7 (DNAStar, Madison, WI, USA). Putative open reading frames (ORFs) were predicted using NCBI's ORF finder (https://www.ncbi.nlm.nih.gov/ orffinder/) with a minimal ORF length of 150 nt, followed by manual inspection. The sequences of the 5' untranslated region (5'-UTR) and 3'-UTR were defined, and the leader sequence, the leader and body transcriptional regulatory sequence (TRS) were identified as previously described [22] . The cleavage of the 16 nonstructural proteins coded by ORF1ab was determined by alignment of aa sequences of other CoVs and the recognition pattern of the 3C-like proteinase and papain-like proteinase. Phylogenetic trees based on nt or aa sequences were constructed using the maximum likelihood algorithm with bootstrap values determined by 1000 replicates in the MEGA 6 software package [23] . Full-length genome sequences obtained in this study were aligned with those of previously reported alpha-CoVs using MUSCLE [24] . The aligned sequences were scanned for recombination events by using Recombination Detection Program [25] . Potential recombination events as suggested by strong p-values (<10 -20 ) were confirmed using similarity plot and bootscan analyses implemented in Simplot 3.5.1 [26] . The number of synonymous substitutions per synonymous site, Ks, and the number of nonsynonymous substitutions per nonsynonymous site, Ka, for each coding region were calculated using the Ka/Ks calculation tool of the Norwegian Bioinformatics Platform (http://services.cbu.uib.no/tools/kaks) with default parameters [27] . The protein homology detection was analyzed using HHpred (https://toolkit.tuebingen.mpg.de/#/tools/hhpred) with default parameters [28] . A set of nested RT-PCRs was employed to determine the presence of viral subgenomic mRNAs in the CoV-positive samples [29] . Forward primers were designed targeting the leader sequence at the 5'-end of the complete genome, while reverse primers were designed within the ORFs. Specific and suspected amplicons of expected sizes were purified and then cloned into the pGEM-T Easy vector for sequencing. Bat primary or immortalized cells (Rhinolophus sinicus kidney immortalized cells, RsKT; Rhinolophus sinicus Lung primary cells, RsLu4323; Rhinolophus sinicus brain immortalized cells, RsBrT; Rhinolophus affinis kidney primary cells, RaK4324; Rousettus leschenaultii Kidney immortalized cells, RlKT; Hipposideros pratti lung immortalized cells, HpLuT) generated in our laboratory were all cultured in DMEM/F12 with 15% FBS. Pteropus alecto kidney cells (Paki) was maintained in DMEM/F12 supplemented with 10% FBS. Other cells were maintained according to the recommendations of American Type Culture Collection (ATCC, www.atcc.org). The putative accessory genes of the newly detected virus were generated by RT-PCR from viral RNA extracted from fecal samples, as described previously [30] . The influenza virus NS1 plasmid was generated in our lab [31] . The human bocavirus (HBoV) VP2 plasmid was kindly provided by prof. Hanzhong Wang of the Wuhan Institute of Virology, Chinese Academy of Sciences. SARS-CoV ORF7a was synthesized by Sangon Biotech. The transfections were performed with Lipofectamine 3000 Reagent (Life Technologies). Expression of these accessory genes were analyzed by Western blotting using an mAb (Roche Diagnostics GmbH, Mannheim, Germany) against the HA tag. The virus isolation was performed as previously described [12] . Briefly, fecal supernatant was acquired via gradient centrifugation and then added to Vero E6 cells, 1:10 diluted in DMEM. After incubation at 37°C for 1 h the inoculum was replaced by fresh DMEM containing 2% FBS and the antibiotic-antimycotic (Gibco, Grand Island, NY, USA). Three blind passages were carried out. Cells were checked daily for cytopathic effect. Both culture supernatant and cell pellet were examined for CoV by RT-PCR [17] . Apoptosis was analyzed as previously described [18] . Briefly, 293T cells in 12-well plates were transfected with 3 µg of expression plasmid or empty vector, and the cells were collected 24 h post transfection. Apoptosis was detected by flow cytometry using by the Annexin V-FITC/PI Apoptosis Detection Kit (YEASEN, Shanghai, China) following the manufacturer's instructions. Annexin-V-positive and PI-negative cells were considered to be in the early apoptotic phase and those stained for both Annexin V and PI were deemed to undergo late apoptosis or necrosis. All experiments were repeated three times. Student's t-test was used to evaluate the data, with p < 0.05 considered significant. HEK 293T cells were seeded in 24-well plates and then co-transfected with reporter plasmids (pRL-TK and pIFN-βIFN-or pNF-κB-Luc) [30] , as well as plasmids expressing accessory genes, empty vector plasmid pcAGGS, influenza virus NS1 [32] , SARS-CoV ORF7a [33] , or HBoV VP2 [34] . At 24 h post transfection, cells were treated with Sendai virus (SeV) (100 hemagglutinin units [HAU]/mL) or human tumor necrosis factor alpha (TNF-α; R&D system) for 6 h to activate IFNβ or NF-κB, respectively. Cell lysates were prepared, and luciferase activity was measured using the dual-luciferase assay kit (Promega, Madison, WI, USA) according to the manufacturer's instructions. Retroviruses pseudotyped with BtCoV/Rh/YN2012 RsYN1, RsYN3, RaGD, or MERS-CoV spike, or no spike (mock) were used to infect human, bat or other mammalian cells in 96-well plates. The pseudovirus particles were confirmed with Western blotting and negative-staining electromicroscopy. The production process, measurements of infection and luciferase activity were conducted, as described previously [35, 36] . The complete genome nucleotide sequences of BtCoV/Rh/YN2012 strains RsYN1, RsYN2, RsYN3, and RaGD obtained in this study have been submitted to the GenBank under MG916901 to MG916904. The surveillance was performed between November 2004 to November 2014 in 19 provinces of China. In total, 2061 fecal samples were collected from at least 12 Rhinolophus bat species ( Figure 1A ). CoVs were detected in 209 of these samples ( Figure 1B and Table 1 ). Partial RdRp sequences suggested the presence of at least 8 different CoVs. Five of these viruses are related to known species: Mi-BatCoV 1 (>94% nt identity), Mi-BatCoV HKU8 [37] (>93% nt identity), BtRf-AlphaCoV/HuB2013 [11] (>99% nt identity), SARSr-CoV [38] (>89% nt identity), and HKU2-related CoV [39] (>85% nt identity). While the other three CoV sequences showed less than 83% nt identity to known CoV species. These three viruses should represent novel CoV species. Virus isolation was performed as previously described [12] , but was not successful. identity). While the other three CoV sequences showed less than 83% nt identity to known CoV species. These three viruses should represent novel CoV species. Virus isolation was performed as previously described [12] , but was not successful. We next characterized a novel alpha-CoV, BtCoV/Rh/YN2012. It was detected in 3 R.affinis and 6 R.sinicus, respectively. Based on the sequences, we defined three genotypes, which represented by RsYN1, RsYN3, and RaGD, respectively. Strain RsYN2 was classified into the RsYN3 genotype. Four full-length genomes were obtained. Three of them were from R.sinicus (Strain RsYN1, RsYN2, and RsYN3), while the other one was from R.affinis (Strain RaGD). The sizes of these 4 genomes are between 28,715 to 29,102, with G+C contents between 39.0% to 41.3%. The genomes exhibit similar structures and transcription regulatory sequences (TRS) that are identical to those of other alpha-CoVs ( Figure 2 and Table 2 ). Exceptions including three additional ORFs (ORF3b, ORF4a and ORF4b) were observed. All the 4 strains have ORF4a & ORF4b, while only strain RsYN1 has ORF3b. The replicase gene, ORF1ab, occupies~20.4 kb of the genome. The replicase gene, ORF1ab, occupies~20.4 kb of the genome. It encodes polyproteins 1a and 1ab, which could be cleaved into 16 non-structural proteins (Nsp1-Nsp16). The 3'-end of the cleavage sites recognized by 3C-like proteinase (Nsp4-Nsp10, Nsp12-Nsp16) and papain-like proteinase (Nsp1-Nsp3) were confirmed. The proteins including Nsp3 (papain-like 2 proteas, PL2pro), Nsp5 (chymotrypsin-like protease, 3CLpro), Nsp12 (RdRp), Nsp13 (helicase), and other proteins of unknown function ( Table 3 ). The 7 concatenated domains of polyprotein 1 shared <90% aa sequence identity with those of other known alpha-CoVs ( Table 2 ), suggesting that these viruses represent a novel CoV species within the alpha-CoV. The closest assigned CoV species to BtCoV/Rh/YN2012 are BtCoV-HKU10 and BtRf-AlphaCoV/Hub2013. The three strains from Yunnan Province were clustered into two genotypes (83% genome identity) correlated to their sampling location. The third genotype represented by strain RaGD was isolated to strains found in Yunnan (<75.4% genome identity). We then examined the individual genes ( Table 2) . All of the genes showed low aa sequence identity to known CoVs. The four strains of BtCoV/Rh/YN2012 showed genetic diversity among all different genes except ORF1ab (>83.7% aa identity). Notably, the spike proteins are highly divergent among these strains. Other structure proteins (E, M, and N) are more conserved than the spike and other accessory proteins. Comparing the accessory genes among these four strains revealed that the strains of the same genotype shared a 100% identical ORF3a. However, the proteins encoded by ORF3as were highly divergent among different genotypes (<65% aa identity). The putative accessory genes were also BLASTed against GenBank records.

Question: What animals do gamma and delta coronavirus mainly infect?
birds