{'fraction_non_alphanumeric': 0.04674537881656397, 'fraction_numerical': 0.037795600839902244, 'mean_word_length': 4.774597495527728, 'pattern_counts': {'":': 0, '<': 3, '<?xml version=': 0, '>': 7, 'https://': 0, 'lorem ipsum': 0, 'www.': 3, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 3, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Prevalence Study and Genetic Typing of Bovine Viral Diarrhea Virus (BVDV) in Four Bovine Species in China

To determine the nationwide status of persistent BVDV infection in different bovine species in China and compare different test methods, a total of 1379 serum samples from clinical healthy dairy cattle, beef cattle, yaks (Bos grunniens), and water buffalo (Bubalus bubalis) were collected in eight provinces of China from 2010 to 2013. The samples were analyzed using commercial antibody (Ab) and antigen (Ag) detection kits, and RT-PCR based on the 5’-UTR and Npro gene sequencing. Results showed that the overall positive rates for BVDV Ab, Ag and RT-PCR detection were 58.09% (801/1379), 1.39% (14/1010), and 22.64% (146/645), respectively, while the individual positive rates varied among regions, species, and farms. The average Ab-positive rates for dairy cattle, beef cattle, yaks, and water buffalo were 89.49% (298/333), 63.27% (248/392), 45.38% (236/520), and 14.18% (19/134), respectively, while the Ag-positive rates were 0.00% (0/116), 0.77% (3/392), 0.82% (3/368), and 5.97% (8/134), respectively, and the nucleic acid-positive rates detected by RT-PCR were 32.06% (42/131), 13.00% (26/200), 28.89% (52/180), and 19.40% (26/134), respectively. In addition, the RT-PCR products were sequenced and 124 5’-UTR sequences were obtained. Phylogenetic analysis of the 5’-UTR sequences indicated that all of the 124 BVDV-positive samples were BVDV-1 and subtyped into either BVDV-1b (33.06%), BVDV-1m (49.19%), or a new cluster, designated as BVDV-1u (17.74%). Phylogenetic analysis based on Npro sequences confirmed this novel subtype. In conclusion, this study revealed the prevalence of BVDV-1 in bovine species in China and the dominant subtypes. The high proportion of bovines with detectable viral nucleic acids in the sera, even in the presence of high Ab levels, revealed a serious threat to bovine health.


Introduction
Bovine viral diarrhea virus (BVDV) is a single-stranded positive-sense RNA virus that is a member of the genus Pestivirus [1] and mainly affects cattle, resulting in fever, diarrhea, leucopenia, reduction in milk yield and reproductive problems [2], or no clinical symptoms, but immunosuppression [3]. Cows infected in the early gestational period with noncytopathic BVDV may produce persistently infected (PI) calves, which are mainly responsible for the spread of BVDV throughout herds via continuous viral shedding from all mucosal surfaces [4,5]. Therefore, identification and removal of such individuals is critical to the success of eradication campaigns [6]. In addition, BVDV is an etiological agent of bovine respiratory disease [7,8].
The BVDV genome contains a single open reading frame (ORF) encoding a polyprotein that is processed co-and post-transnationally into mature viral proteins. This ORF is flanked with 5'-and 3'-untranslated regions (UTRs). Based on the phylogenetic analysis of partial sequences from the 5'-UTR, the N-terminal autoprotease (Npro) or envelope glycoprotein (E2) region of the genome of this virus is usually divided into two distinct genetic species, namely BVDV-1 and BVDV-2. However, a third genetic species, BVDV-3, was recently reported [9], which was described as atypical BVDV [10]. BVDV-1 is further classified into 18 potential genetic subtypes, 1a-1t [11][12][13][14][15][16] and BVDV-2 into three subtypes (2a-2c) [17][18][19]. Originally, BVDV-2 was reported to be related to severe hemorrhagic disease, resulting in high mortality, in Canada [20], while BVDV-1 varies in virulence, as most BVDV-1 viruses only cause asymptomatic infection [21]. The high degree of antigenic and genetic diversity of BVDV causes major diagnostic and prophylactic difficulties because common diagnostic tests and vaccine production are based on viral antigens (Ags) [22]. Therefore, recognition of the variability of BVDV field strains is crucial when designing a successful control or eradication scheme at the herd level [6].
In China, the first description of BVDV infection in cattle dates back to 1980 when strain Changchun 184 was isolated from an aborted fetus. Since then, BVDV infection has been reported in beef and dairy cattle, yaks (Bos grunniens), water buffalo (Bubalus bubalis), camels, Sika deer (Cervus nippon), and swine in more than 20 regions of China with a high seroprevalence in all these species. Furthermore, the BVDV subtypes currently circulating in China are very diverse and include eight subtypes: 1a, 1b, 1c, 1d, 1m, 1o, 1p, and 1q [14,[23][24][25][26]. However, previous studies were mainly performed in northern and western China, since these areas are historically the main regions with beef and dairy cattle production.
In the last decade, increasing demand by consumers and continuously increasing prices of beef and dairy products have promoted investments in cattle ranching in southern China. This change has led to an increasing number of cattle transported from the North and West to southern parts of China and has thereby altered the distribution of diseases, including BVDV infection. Therefore, nationwide re-evaluation of BVDV prevalence in bovines is important to control BVDV infection. Second, among bovine species, surveillance of BVDV infection in water buffalo and yaks less frequent compared to that in dairy cattle. Determination of the prevalence of BVDV in these species is critical to control infection. Third, in previous reports, various surveillance methods have been used by different investigators, thus the results are not comparable. Therefore, it is necessary to compare available methods in parallel to recommend the most appropriate to determine the prevalence of BVDV infection. Finally, different surveillance methods based on varied mechanisms were conducted in parallel, thus these results are expected to help determine the true prevalence of BVDV infection by evaluating the prevalence of antibodies (Abs), antigens (Ags), and nucleic acids. The aim of this study was to facilitate an evidence-based BVDV control strategy in China.

Sample collection
This study was performed in strict accordance with the Hubei Regulations for the Administration of Affairs Concerning Experimental Animals, 2005. The study protocol was approved by the China Hubei Province Science and Technology Department (permit no: SYX-K(ER) 2010-0029). A total of 1379 serum samples were collected from eight provinces located in different regions of China from 2010 to 2013 (Table 1 and Fig 1A) which included Guangxi province in southern China, Inner Mongolia and Liaoning provinces in northern China, Qinghai province and Tibet in western China, Jiangsu province in eastern China, and Hubei and Henan provinces in central China. The species included dairy cattle, beef cattle, yaks, and water buffalo, which varied according to the primary species of the local cattle industry. The animals ranged from 2-5 years old and were clinically healthy and exhibited no reproductive problems. The herds were not vaccinated against BVDV because no commercial BVDV vaccines are available and no control programs are currently implemented in China. However, the yaks were more likely to be vaccinated with a live vaccine against classical swine fever virus (CSFV) to prevent BVDV infection due to the high incidence of BVD. Unfortunately, details regarding this heterogenic vaccination schemes were not available because the use of live CSFV vaccines is unauthorized in China. Water buffalo in Guangxi province and yaks in Qinghai province are commonly farmed species, but the herd size is usually small and the animals are free-range. Samples were collected by the local Veterinary Service Agencies for disease surveillance. In northern China, samples from beef cattle were consecutively collected at local slaughter houses for 6-9 days to ensure a sufficient sample number and diversity. In eastern and central China, samples were collected from 10% of adult cattle from one or two representative dairy farms. Blood was collected from the jugular or caudal vein of each animal and serum was isolated, transported under cool condition to our laboratory, and stored at -20°C until assayed.

Serological and antigen detection
A total of 1379 serum samples were collected and initially screened for Abs to BVDV using a commercial BVDV Ab test kit (IDEXX Laboratories, Inc., Liebefeld, Switzerland). Then, 1010 of these samples were subjected to Ag detection using the IDEXX SNAP BVDV Antigen Test kit (IDEXX Laboratories, Inc.) according to the manufacturer's recommendations. Ag detection of the other samples was not performed because the volumes were insufficient after the first round of testing.

Nested RT-PCR analysis
Serum samples were further divided into three groups according to the results of Ab detection, namely Ab-positive, Ab-negative and suspicious groups, and samples from each group were selected for RT-PCR analysis (S1 Table).
To confirm the typing results based on 5'-UTR sequences, a 411-bp product containing the Npro region was amplified by nested RT-PCR using the outer primers B32/B31 and inner primers BD1/BD3 [13,27]. For verification of newly identified subtype isolates, the outer primers B32/B31 and inner primers BD1u/BD3u were used. The primers BD1u/BD3u were specific to BVDV-1 strain M31182 isolated from a yak in Sichuan, China (GenBank accession no. JQ799141).
All primers used in this study are listed in S2 Table. RNA isolation and cDNA synthesis Total RNA was extracted from 140 μL of serum using the TIANamp virus RNA kit (Tiangen Biotech (Beijing) Co., Ltd., Beijing, China) according to the manufacturer's instructions and stored at -70°C until assayed. Reverse transcription (RT) was performed to produce cDNA using the reagents of the PrimeScript RT reagent kit with gDNA Eraser (Takara, Otsu, Shiga, Japan) following the manufacturer's protocol.

Nested PCR of cDNA and sequencing
For BVDV-1 detection, primary PCR was performed using a total volume of 25 μL containing 12.5 μL of PCR Mix (Dongsheng Biotech, Guangzhou, China), 2 μL of cDNA, 9.5 μL of sterilized H 2 O, and 0.5 mM each of the primers BVDV-1 F and BVDV-1 R. The reaction was carried out at 94°C for 5 min, followed by 35 cycles of 94°C for 30 s, 56°C for 30 s, and 72°C for 30 s, with a final elongation step of 72°C for 10 min. A 2 μL aliquot from the primary PCR was used as a template for the second PCR and the reagents and cycling conditions for the secondary PCR were the same as those for the primary PCR except 30 cycles were used. For Npro gene amplification, the reaction mixtures of primary and secondary PCR were the same as those described above. The primary PCR was carried out at 94°C for 4 min, followed by 35 cycles of 94°C for 1 min, 50°C for 30 s, and 72°C for 40 s, with a final elongation step of 72°C for 10 min, while the secondary PCR was similar except the annealing temperature was increased to 58°C.For BVDV-2 detection, the reaction mixture was the same as described above. The primary PCR reaction was included 30 cycles of 95°C for 30 s, 55°C for 30 s, and 72°C for 30 s with a final elongation step of 72°C for 10 min. The nested amplification conditions were the same as those for the primary PCR reaction, except that the annealing temperature was 53°C. Two BVDV Ag-positive samples (NMG313-1 and NMG314-65) detected by the IDEXX SNAP BVDV Antigen Test kit in this study were used as positive controls and commercial fetal bovine serum (Gibco, Grand Island, NY, USA), which was confirmed with the IDEXX SNAP BVDV Antigen Test kit and RT-PCR, was used as a negative control. The PCRamplified amplicons were checked by electrophoresis on 1% agarose gel with EB stain, purified using the TIANgel Midi Purification Kit (Tiangen Biotech (Beijing) Co., Ltd.) and sequenced by Shanghai Sangon Biological Engineering Technology & Services Co., Ltd. (Shanghai, China) using the primers listed in S2 Table with an ABI automated 3730 sequencer (Applied Biosystems, Foster City, CA, USA). Nucleotide sequences were aligned using SeqMan II sequence assembly and analysis software (DNASTAR, Inc., Madison, WI, USA). Similarities in nucleotide sequences were evaluated using the MegAlign program (DNASTAR, Inc.) and the specified sequences against existing sequences were further identified using the Basic Local Alignment Search Tool (http://www.ncbi.nlm.nih.gov/blast/Blast.cgi).

Phylogenetic analysis
Nucleotide sequences of the BVDV-1 5'-UTR (200 bp) and Npro gene (411 bp) fragments were aligned using the Clustal W program (www.clustal.org/). Further phylogenetic analysis was performed with the neighbor-joining method using the MEGA5 program [28,29], and evolutionary distances were calculated using the Kimura 2-parameter method. The robustness of the phylogenetic analysis and the significance of branch order were determined using the bootstrapping method based on 1000 replicates.
A total of 37 reference sequences of known BVDV-1 and BVDV-2 strains were retrieved from the NCBI GenBank database (http://www.ncbi.nlm.nih.gov/genbank). Sequences of BVDV isolates in this study were deposited in the GenBank database under the accession numbers KJ578795-KJ578918 (S3 Table). The nucleotide sequence of the Npro gene was analyzed using the same methods and parameters. The Npro gene sequences were submitted to the Gen-Bank database under the accession numbers KP126233-KP126243.

Statistical analysis
Prevalence was defined as the proportion of positive animals that were tested. The chi-squared test was used to analyze differences in prevalence between two groups and a probability (p) value < 0.05 was considered statistically significant ( Ã ) and p < 0.01 as very significant ( ÃÃ ).

BVDV Ag detection
A total of 1010 serum samples, which was less than the number used for Ab detection because the volume of the other 369 samples was insufficient for Ab detection, were subjected to BVDV Ag detection, which showed that only a few samples were positive with an overall rate of 1.39% (14/1010) (95% CI: 0.8%-2.3%). Contrary to the Ab profile, the positive rate of water buffalo samples was significantly highest at 5.97% (95% CI: 2.6%-11.4%) compared to that of the other species (p < 0.01) ( Table 3), while that of dairy cattle was lowest at 0.00% for 116 samples (95% CI: 0%-3.1%). For beef cattle and yaks, the Ag-positive rates lied between these values, but were less than 1% (Table 3). Furthermore, 13 of 14 Ag-positive samples were Ab-negative (S4 Table).

Nested RT-PCR detection of BVDV infection
A total of 645 samples were tested for a specific fragment of the 5'-UTR of BVDV-1. These samples were chosen from the 1243 sera samples left after the first round of testing for BVDV-1 antibodies covering Ab-positive, Ab-negative, and suspected samples (S1 Table). As shown in Table 4, 146 (22.6%) of the 645 serum samples were positive and the RT-PCR products were of the right size, as indicated by electrophoresis on 1% agarose gels (S1 Fig) and had correct sequences, as indicated by aligning the sequences with the published reference sequences retrieved from the GenBank database (S3 Table).
A comparison between Ag-positive and RT-PCR detection showed that 9 (64.3%) of the 14 Ag-positive samples were positive for BVDV nucleic acid, while five samples, including one of three from yaks and four of eight from water buffalo, were not confirmed by RT-PCR (S4 Table). The agreement between Ab and RT-PCR detection was compared (Table 4) Table 4). While we attempted to sequence all 146 positive samples, only 124 were successful due to technical problems, such as double signals during the sequencing process, encountered while processing the other 22 positive samples (S3 Table). All sequences belonged to BVDV-1. The BVDV-1 subtypes for each region were further analyzed (Table 5), which identified two main subtypes, namely 1b and 1m, collectively accounting for 82.26% (102/124) of all strains. Among these, 33.06% (41/124) were 1b and 49.19% (61/124) were 1m, respectively. The clinical strains classified as 1b shared a sequence homology of 94.0%-99.0% with the reference strains, while those classified as 1m shared a sequence homology of 89.7%-96.0%. The other 17.74% (22/124) BVDV-1 sequences did not match the known subtypes and clustered in a new branch, with a nucleotide homology of 97.5%-100.0%, designated as BVDV-1u. The sequences of these 22 samples were submitted to a BLAST search using the blastn algorithm. The results showed that they shared a homology of 91.5%-93.5% to the Chinese strain M31182 (GenBank accession no.: JQ799141.1) isolated from a yak in Sichuan, China (S2 Fig). However, this strain had not yet been classified into a known subtype [30] A phylogenetic tree of the representative clinical strains of each province and the reference strains (S3 Table) was constructed (Fig 2A) and the geographic distribution of the BVDV-1 subtypes was plotted (Fig 1B).

Prevalence of BVDV Infection in Chinese Bovines
To confirm the typing results based on 5'-UTR sequences, Npro sequences of the 11 samples (one from group 1b, seven from group 1m and three from group 1u, the new subtype) were compared. Phylogenetic analysis was performed by comparing a 411-bp region of Npro corresponding to nt 386-796 of our 11 samples with data available for BVDV reference strains representing BVDV subtypes 1a-1r, 2, and 3. The phylogenetic analysis results confirmed the classification determined by the 5'-UTR sequences (Fig 2B). Briefly, the sample JS-05059 was classified as subtype 1b and shared a sequence homology of 87.1%-88.1% with known sequences of the same subtype; samples NMG313-1, NMG314-65, XZ-84, QHTJ-303887, XZ-103, XZ-109, and XZ-25 belonged to subtype 1m and shared a sequence homology of 93.9%-95.6% with known sequences of the same subtype; and samples JS-03198, JS-03148 and GXBH-EB34 were clustered to the 1u subtype and shared a sequence homology of 89.8% with the Chinese strain M31182 retrieved from the GenBank database. In addition, the 5'-UTR sequences of these three isolates had 93.5% identity with that of strain M31182.

Discussion
It is difficult to determine the exact extent of the BVDV epidemic among different bovine species in China due to the variety of detection tests, sampling methods, species, and locations in individual reports. Generally speaking, most areas of China have reported the detection of BVDV [31]. For large-scale farms, the average seroprevalence rate reached 92.5% among dairy cows and 29.8% among beef cattle in Fujian province in southern China [32]. In western China, the seroprevalence of yaks was reportedly 53.65% in Tibet and 72.14% in Qinghai province [33]. The positive ratio of neutralization Ab in water buffalo averaged 17.25% in some areas [34]. These results demonstrated a high and variable seroprevalence of BVDV-1 among Chinese bovine species. Our results confirmed this status by demonstrating a seroprevalence of 89.49% (298/333) for dairy cows, 63.27% (248/392) for beef cattle, 45.38% (236/520) for yaks, and 14.18% (19/134) for water buffalo. These results are within previously reported ranges, but the seroprevalence of beef cattle was greatly increased in our study, which was probably due to the frequent movement of calves because of the rapid expansion of the beef industry. In addition, the seroprevelance among yaks in Qinghai province (51.36%) and Tibet (30.92%) in this study were lower than previously reported (72.14% and 53.65%, respectively) [33]. We propose the following reasons for these differences: (i) the samples came from different regions; and (ii) in some regions, yaks vaccinated with live CSFV vaccine might have improved the false positive rate in the cited previous reports or played a role in the control of BVDV spread leading to the later decrease in the seropositive rates in this study. According to local veterinarians, unauthorized use of live CSFV vaccine is administered to about 15%-20% of the yaks in some areas to prevent BVD outbreaks [26]. However, records for this vaccination were not available. In addition, it may be possible that some false positive results occurred because the test kit for BVDV Ab detection was validated in dairy and beef cattle, not yaks and water buffalo.
As proposed by Houe [35], the prevalence of BVDV Ab in cattle herds varies greatly and corresponds to the five phases of the infection cycle. We suspect that the herds sampled in this study might have been in phase B (infected herd with PI calves younger than 3-4 months old and most acute infections occur at varying rates, depending on housing of the animals) and phase C (an infected herd with PI calves older than 3-4 months old, in which seropositivity can reach > 90% with no eradication program).

Correlations between Ab and Ag testing
The concept that Ab-positive rates are inversely correlated to Ag-positive rates has be confirmed in infections with BVDV [36] and other agents s [37]. The results of this study supported these previous findings. For instance, the absolute Ab-positive rates were high, while the Ag-positive rates were low, indicating infrequent reactivation of BVDV from the latent status under a background of high Ab levels. Furthermore, the decreasing order of Ab-positive rates (dairy cattle (84.49%) > beef cattle (63.27%) > yaks (45.38%) > water buffalo (14.18%)) was opposite to that of Ag-positive rates (water buffalo (5.97%) > yaks (0.82%) > beef cattle (0.77%) > dairy cattle (0%)). Accordingly, among 14 Ag-positive samples, 13 were Ab-negative (S4 Table).
The Ag-positive rate in this study might have been underestimated because, generally speaking, the sensitivity of the colloidal gold-labeled test strip is usually less than that of an enzymelinked immunosorbent assay (ELISA). Second, according to the product data sheets, the type of sample may affect sensitivity. For example, the sensitivity of a serum sample can be 4% lower than that using tissue from a small ear notch. However, as mentioned above, the test kit for BVDV Ag detection has not been sufficiently validated in water buffalo or yaks; therefore, a higher rate of water buffalo than other species might positively influence the results. Overall, the average Ag prevalence rate of 1.4% (14/1010) for these four species in our study was close to the rates in previous reports from China and other countries [35,38].

Correlations between Ag, Ab, and nucleic acid detection
In the present study, RT-PCR demonstrated a prevalence of 22.64% (146/645) among all bovine species, which ranged from 13.00% (26/200) for beef cattle to 32.06% (42/131) for dairy cattle. These results were in agreement with those in previous reports. For instance, in previous reports, RT-PCR revealed a prevalence of BVDV of 24% (98/407) among yaks from Qinghai province [26,33], 26.85% (105/391) among dairy cattle from Ningxia province in northwestern China [25], and 29.03% (18/62) among beef cattle from northern and Eastern China [14]. Taken altogether, the proportion of bovines with viremia in China was high, although the Abpositive rates were high.
Reportedly, RT-PCR analysis of BVDV nucleic acids is more sensitive than Ag detection. For example, Oem et al [36]detected BVDV nucleic acid in 15.5% of brain samples by RT-PCR, while only 2.9% of positive samples were confirmed by Ag detection via immunohistochemical analysis and Ag capture ELISA. Similarly, there was a large difference in positive rates between BVDV nucleic acid and Ag detection in this study.
A second possible reason for the difference between results by RT-PCR and Ag detection could be that RT-PCR can detect the presence of virus in both acutely and persistently infected animals, while Ag detection was designed for only persistently infected animals. Therefore, not all samples deemed positive by RT-PCR can be detected using the Ag test [36].
On the other hand, among the 14 Ag-positive samples in this study, 64% (9/14) were according to RT-PCR results, although the five samples not confirmed by RT-PCR may be false positives. Coincidently, these five samples included four from water buffalo and one from a yak. This fact again aroused suspicion that the Ag-capture test might be problematic in testing yaks and water buffalo. According to the product instructions, retesting of Ag-positive samples after 3 weeks is recommended to exclude false positive samples. However, resampling and testing is sometimes difficult because the cattle might be slaughtered or sold, or the travel distance is economically unviable.
As shown in Table 4 and reported elsewhere [36], of the cattle with positive results by to RT-PCR, both Ab-positive and-negative results were detected. The double positive status to both BVDV Ab and nucleic acid detection might indicate an acute infection stage, while that positive to nucleic acid but negative to Ab detection might suggest either an early phase of acute infection, because development of detectable Ab usually requires 1-2 weeks after infection, or a persistent infection in which Ab detection is negative, but positive for both nucleic acid and Ag detection. The fact that more than 30% of dairy cows, 2-5 years, in China were viremic for BVDV, particularly since all of these animals appeared "healthy," was astonishing. However, considering the fact that there is no vaccination program or eradication plan in China, the virus will continue to persistently circulate among bovine populations without interference, thus it is not difficult to understand the high prevalence. Furthermore, the dominant genotype was BVDV-1 in this study, although this genotype usually does not cause serious clinical illness and thus infected bovines may understandably appear "healthy." Although BVDV infection is known to cause reproductive and productive abnormalities, there was no indication of such problems among the sampled herds.
On the other hand, the status of negative to RT-PCR detection, but positive to Ab detection might represent a stage after clearance of acute infection without persistence. As previously reported, PI with BVDV in cattle only occurs during fetal development, thus the status can be maintained for long periods until the PI cattle become pregnant and give birth to PI calves [39]. These PI cattle then develop immune tolerance at the fetal stage. However, acute infection of bovines without PI can induce protective immunity and finally clear the virus after acute infection.
In this study, we demonstrated that 1b (33%) and 1m (49%) were dominant BVDV subtypes, accounting for 82% (102/124) of the total. These results were in agreement with those of previous reports. BVDV-1a, 1c, 1d, 1o, 1p, and 1q were not detected in this study. Theoretically, the primers used for 5'-UTR amplification in this study could detect other subtypes, including 1a, 1c, 1o, 1p, and 1q [24]. Actually, we detected 1a, 1c, and 1p once each from clinical samples of diseased cattle during routine diagnosis using these primers. In this survey, all samples were collected from clinically healthy animals. Therefore, we suspect that other factors, such as individual herds, geographical distribution of the herds, sampling size, and health status, might affect the detection results. Regarding subtype 1q, a high prevalence was found in yaks in Qinghai province [26], dairy farms in Ningxia province [25], and pigs [24]. However, this study failed to detect this subtype for several possible reasons, including: (i) the samples were collected in northeast of Qinghai province in this study, whereas samples were collected in southeast areas of Qinghai and Ningxia provinces by Gong et al. [25,26]; (ii) more samples (n = 407) were collected from yaks than in our study (n = 407 vs. 120, respectively), thus a greater number of samples should be tested to confirm the prevalence of BVDV subtype 1q in future studies; and (iii) about 15%-20% of the yaks received live CSFV vaccine, thus and subtype 1q may have been introduced by BVDV-contaminated live CSFV vaccine, as discussed by Gong et al. [26]. To date, subtype 1o has only been detected in pigs.
In addition, this study is the first to report that the new subtype BVDV-1u comprised 18% of all the detected strains. This subtype shared a high homology in the sequences of the 5'-UTR (93.5%) and Npro (89.8%) genes with the untyped strain M31182 (GenBank accession no.: JQ799141.1), which was isolated from a yak in Sichuan province, located in central China, in 2010 [30]. Other than the one subtype 1u isolate detected from a yak in Qinghai province, we also detected this subtype in three other species from other regions, including seven samples from water buffalo in Guangxi province, ten from dairy cattle in Jiangsu province and one from Hubei province, and two from beef cattle in Inner Mongolia and one from Liaoning province, respectively. This is the first report to describe the distribution of BVDV-1u in varied bovine populations from different areas in China.  Table. (TIF) S1