Advances in Animal and Veterinary Sciences

Download PDF Download ePUB
AAVS_9_3_446-452

 

 

Research Article

 

New Epizootic of Lumpy Skin Disease in Assiut- Egypt: Molecular Identification and Characterization

 

Khaled A. S. El-Khabaz1*, Eman A. M. Shosha2, Aml M. Abdel-Ra’ouf3

1Department of Animal Medicine (Infectious Diseases), Faculty of Veterinary Medicine, Assiut University, Egypt; 2Department of Microbiology and Immunology, Faculty of Veterinary Medicine, New Valley University, Egypt; 3Department of Animal Medicine (Infectious Diseases), Faculty of Veterinary Medicine, Aswan University, Egypt.

 

Abstract | An outbreak of lumpy skin disease (LSD) was recorded during the period of 2018 in the study area. Tissue samples from skin nodules were surgically collected from diseased animals (N=23) for molecular identification of the causative virus by conventional polymerase chain reaction (PCR) using specific primers designated from G protein coupled chemokine receptor (GPCR) gene. All the skin samples were found to be positive for lumpy skin disease virus (LSDV) as they gave the expected amplicon size at 554 bp. Sequencing analysis of the obtained amplified segment disclosed that the present outbreak virus has a very high identity percentage >99-100% with virulent LSDV either that’s isolated from Egypt or different foreign countries, while the percent of identity were 95.62% with the vaccinal strains of LSDV. Multiple alignment and Phylogenetic analysis of the obtained sequence manifesting that members of capripoxviruses could be segregated in different groups and the present outbreak strain grouped within the LSDV virulent strains cluster.

 

Keywords | Lumpy skin disease, capripoxvirus, GPCR gene, phylogenetic analysis, PCR

 

Received | December 07, 2020; Accepted | December 17, 2020; Published | January 15, 2021

*Correspondence | Khaled AS El-Khabaz, Department of Animal Medicine (Infectious Diseases), Faculty of Veterinary Medicine, Assiut University, Egypt; Email: khaled.sayed@aun.edu.eg

Citation | El-Khabaz KAS, Shosha EAM, Abdel-Raouf AM (2021). New epizootic of lumpy skin disease in assiut- egypt: molecular identification and characterization. Adv. Anim. Vet. Sci. 9(3): 446-452.

DOI | http://dx.doi.org/10.17582/journal.aavs/2021/9.3.446.452

ISSN (Online) | 2307-8316; ISSN (Print) | 2309-3331

Copyright © 2021 El-Khabaz et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

 

Introduction

 

Lumpy skin disease (LSD) is a highly important insect-borne cattle disease caused by a virus belonging to family poxviridae, genus capripox. The genus capripox include lumpy skin disease virus (LSDV), sheeppox virus (SPPV), and goatpox virus (GTPV), which are genetically related and host-specific (Le Goff et al., 2005).

 

Clinically LSD characterize by rise of body temperature, anorexia, appearance of intradermal skin nodules all over the body, loss of condition and enlargement of superficial lymph nodes (Ochwo et.al., 2020). There was a variation in the intensity of clinical signs according to pathogenicity of causative virus, breed, age, immunity, and production state of host, and vector spreading (OIE, 2018).

 

Diagnosis of LSD is usually based upon the pathognomic clinical signs and different laboratory techniques such as, virus isolation and identification, several PCR techniques and serological diagnosis. Serological diagnosis is time-devouring and cannot differentiate between different members of capripoxviruses (El-Bagoury et al., 2009).

 

The first outbreak of LSD in Egypt was in the year 1988 after the importation of cattle from Somalia (House et al., 1990) after that, several epidemics were recorded in different years (El-Kholy et al., 2008; Raof et al., 2010; Sharawi and Abd El-Rahim, 2011; El-Khabaz, 2014; Tamam et al., 2018; Rouby et al., 2019; Hodhod, 2020).

Previously, the disease was a threat to Africa only, but now it is well known that it has been spread to Asia (Rashid et. al., 2017; Kasem et.al., 2018) and Europe (Tasioudi et.al., 2016; Toplak et al., 2017; Kononov et al., 2019).

 

GPCR are transmembrane proteins play a great role in the transmission of signals of various stimuli such as light, odorants, nucleotide and proteins from extracellular to intracellualar (Kostenis, 2004). Capripoxviruses homologue of GPCR responsible for cell proliferation in addition to suppression of the antiviral effect elicited by the host (Cao et al., 1995; Tulman et al., 2002). By phylogenetic analysis of GPCR gene, members of genus capripoxvirus can be clearly segregated into different obvious sets (Le Goff et al., 2009).

 

This study aims to molecularly identify lumpy skin disease virus, which has been circulating in Assiut governorate during the year 2018, by partial amplification of GPCR gene. Also sequencing and phylogenetic analysis of the amplified gene sequence with a subsequent comparison between the obtained LSDV sequences and these available in GenBank.

 

Materials and methods

 

Study area

Assiut governorate which is located in Upper Egypt, characterized by hot and relatively humid weather during summer and slightly cold with rear rains during winter.

 

Ethical approval

The examination and sampling were done without causing any exertion or hurt to the animals.

 

Animals

A total of 23 diseased native and mixed breed cattle were used in this study.

 

Clinical examination

The clinical examination of all diseased animals under study was carried out according to Jackson and Cockcroft (2002).

 

Samples

Including all skin layers (epidermis, dermis, and subcutis), skin nodules biopsies (n=23) were collected aseptically from different 23 diseased cattle of different ages and from different localities of Assiut governorate. All skin nodules samples were kept at -20 C in a deep freezer till used.

 

DNA extraction

DNA was extracted from skin nodules by using G-spin DNA extraction Kit (Intron Cat. no. P53077) according to the manufacturer’s instructions. The extracted DNA was then stored at –20 °C for subsequent PCR assay.

 

Primer

The primers were designated from G-Protein Coupled Chemokine Receptor (GPCR) gene according to (Tamam et al., 2018) with the following sequences:

Forward primer 5- AGT ACA GTT AGT AGC GCA ACC -3and reverse primer: 5- GGG TGA ACT ACA GCT AGG TAT-3, this is to amplify 554 bp.

 

Polymerase chain reaction

It was performed according to kit manual instructions using the thermal cycler (Biorad gradient temperature minicycler). The parameters were: initial denaturation at 95°C for 10 minutes, 40 cycles (denaturation at 95 °C for 30 seconds, annealing at 50°C for 30 seconds, extension at 72°C for 1 minute); and a final extension phase run at 72°C for 10 min.

 

Electrophoreses of PCR products were done in a 2% agarose gel. The agarose gel was prepared by dissolving 1 gram agarose in 50 ml 1X TAE buffer and heated using microwave till becoming semisolid. 2μl ethidium bromide (0.5mg/ml) (INtRON Biotechnology, Korea) was added to each 50ml agarose solution. 10µL of the PCR products as well as the DNA ladder were loaded on the gel after mixing with loading dye, and visualized by UV transilluminator and photographed (BIORAD, USA).

 

Sequencing and phylogenetic analysis

By using QIA quick gel extraction kit (Qiagen, Germany), one of the positive samples by conventional PCR was purified and sequenced in both directions using the same primer. Sequencing was done by 3500 XL Genetic Analyzer in the Genome Research Unit- Animal Health Research Institute (AHRI- Giza, Egypt).

 

Accession number

The gained partial GPCR gene sequence was submitted to the Genbank under Accession No. MT793047.

 

Database search

The obtained nucleotide sequence of partial GPCR gene of LSDV was searched in the nucleotide collection standard databases (National Center for Biotechnology Information, NCBI) using megaBLAST optimized for highly similar sequences option.

 

Multiple alignment and phylogenetic analysis

Multiple alignments were created with clustal W integrated into Molecular Evolutionary Genetics Analysis (MEGA) version X software. The phylogenetic tree was constructed by the use of maximum-likelihood method based on Tamura-Nei model in Molecular Evolutionary Genetics Analysis (MEGA) version X software (Kumar et al., 2018) with Bootstraping confidence calculated at 1000 replicates.

 

Results

 

Clinical features

All the examined animals in this study were showing the characteristic clinical features of lumpy skin disease with variation in severity ranged from mild to severe. The commonly noticed clinical signs were fever, an apparition of skin nodules with variable distribution patterns with necrosis, and sloughing of many of them leaving deep ulcers (Figure 1), loss of condition, and anorexia. In addition to the previous signs, respiratory distress, tarry like diarrhea, edema in forelimbs and/or brisket, and corneal opacity were observed in some of the examined animals in different degrees.

 

 

 

Molecular identification

All the examined skin nodules samples (n=23) from diseased animals were found to be positive for LSD virus when tested by conventional PCR using a specific primer for the GPCR gene. The partial GPCR gene of LSDV was successfully amplified as the anticipated amplicon size 554 bp was clearly recorded in all skin samples from all diseased animals (Figure 2).

 

 

Sequencing and phylogenetic analysis

The BLASTn search analysis showing that the present outbreak sequence of partial GPCR gene has a very high percentage of identity >99-100% with all virulent LSDV isolated from different outbreaks in Egypt (KJ561443.1,

 

 

 

KP071937.1, MF156211.1, MH427384.1, MK736888.1, MH427385.1, MN845071.1, and MN381843.1), and in some different countries such as Turkey (KR024777.1), Sudan (FJ869369.1), KSA (MN422452.1), and Burkina Faso (FJ869352.1) (Figure 3). This percent of identity lowered to be 95.62% with LSDV vaccinal strains (Figure 4). The generated phylogenetic tree has shown four genetically clear clusters lineages (LSDV field isolates, SPPV, GTPV, and LSDV vaccinal strains) (Figure 5), and our epidemic strain (labeled with blue triangle) is found in the LSDV field isolates group.

 

Discussion

 

In Egypt, where LSD is endemic, reliance placed on vaccination using the sheep pox vaccine. Obviously, there were some faults as many epidemics appear at different time intervals. During the summer of 2018, a new epizootic of LSD has been noticed, affecting cattle of different ages and of both sexes, in different localities of the Assiut governorate. The present work directed toward molecular identification and characterization of the causative agent of such an outbreak.

 

In the present study, all the examined animals were found exhibiting the characteristic clinical signs of lumpy skin disease such as a sharp increase in body temperature above the normal level, firm nodules covering a large area of the body surface with necrosis and tearing out some of them, decreased appetite and milk production. Further clinical abnormalities such as corneal opacity, lymphadenopathy, edematous swelling in forelegs, and tarry-like diarrhea were observed in some of the diseased animals. Previously, similar clinical signs were recorded by Tuppurainen et al. (2005), Gari et al. (2010), Tasioudi et al. (2016), Tassew et.al. (2018), Rouby et al. (2019). For confirmation of these clinical suspicious, a conventional PCR (using specific primer) was employed to amplify 554 bp region of GPCR gene LSDV. The results of PCR confirmed that the causative agent of such an outbreak is LSDV. Different studies reported that LSDV was responsible for several epidemics in different Egyptian governorates along the previous few years (El-Bagoury et al., 2009; El-Khabaz, 2014; Abdallah et al., 2018; El-kady et al., 2018; Hodhod et al., 2020). This outbreak was recorded in the study area during summer months where both the temperature and relative humidity are high with an increase of insect population. Humidity and temperature greatly increase the population of insect vector (Ayelet et. al., 2014; Zeynalova et.al., 2016; Kasem et.al., 2018).

 

For its speed and accuracy, PCR is an ideal test for the diagnosis of LSD in infected animals as it efficiently detects LSDV DNA in skin samples. PCR is a very credible test for the diagnosis of LSD for its velocity, sensitivity, and specificity (El-Kholy et al., 2008; Tamam et al., 2018). In the present study, all skin samples collected from diseased animals gave positive results by PCR. LSDV is present in higher concentrations and for long period in the skin so skin samples are more superior than blood or organs for the detection of viral DNA of LSDV (Zeynalova et al., 2016).

 

Due to the continuous importation of live cattle, especially from African countries, the molecular investigation of the origin of LSDV is very important from the epidemiological point of view. The genetic characters of the recognized LSDV, from examined animals in Assiut governorate, were evaluated by interpreting the phylogenetic tree and sequence identity. The obtained sequence and all the Egyptian field isolates were found in the same cluster, indicating that the same strain circulating in different governorates of Egypt, also the stable nature of LSDV. LSDV genome characterized by high conserved nature despite its circulation in different geographical regions (Sameea Yousefi et al., 2018; Adedeji et al., 2019).

 

Database search and blast analysis declare very high sequence conformity >99-100% of this outbreak LSDV sequence and the Egyptian virulent isolates deposited in Genbank. The phylogenetic analysis of the attained sequence of partial GPCR gene was capable to assort capripoxviruses into 4 distinct clusters (LSDV field isolates, SPPV, GTPV, and LSDV vaccinal strains); this is greatly supported with 1000 bootstrap reiterates confidence. Similar findings were recorded by Adedeji et al. (2019) and Ochwo et al. (2020). The field LSDV lineage, including the present outbreak strain, is more closely linked to SPPV, while the vaccinal isolates were very tightly associated with GTPV. LSDV is more correlated to SPPV than to GTPV (Beard et al., 2010; El-Kenawy and El-Tholoth, 2010; Rashid et al., 2017), while Le Goff et al. (2009) reported that LSDV was closely connected with GTPV rather than SPPV. Genetic studies proved that members of genus capripoxvirus had been arisen from one common ancestor during their evolution journey (Tulman et al., 2002; Le Goff et al., 2009).

 

In the current study the phylogenetic analysis based on the GPCR gene showing that all of the field isolates of LSDV were located in separate cluster away from the vaccinal strains, this seems to be helpful in differentiation between the field isolates from the vaccinal ones. Complete GPCR gene sequence analysis may be needed to confirm or deny this hypothesis.

 

Conclusion

 

The first step of controlling infectious diseases depends mainly upon the prompt and precise diagnosis, so it seems that GPCR based PCR is a very useful technique in the diagnosis of LSD in susceptible animals. The LSDV causing this outbreak was widely linked to different virulent strains detected in other Egyptian governorates. Sequencing of GPCR gene isn’t only useful in differentiation between different members of capripoxviruses but also the distinction between filed and vaccinal strains of LSDV.

 

Conflict of interest

 

All the authors declare that they have no conflict of interest.

 

Acknowledgment

 

Nothing to declare.

 

authors contribution

 

KASE and EAMS designed the study and edited the manuscript. KASE and AMA collected the samples. All authors shared in the lab work. KASE performed the interpretation of the results. The final version of the manuscript was read and approved by all authors.

 

References

 

  • Abdallah FM, El Damaty HM, Kotb GF (2018). Sporadic cases of lumpy skin disease among cattle in Sharkia province, Egypt: Genetic characterization of lumpy skin disease virus isolates and pathological findings. Vet. World. 11:1150-1158. https://doi.org/10.14202/vetworld.2018.1150-1158
  • Adedeji AJ, Möller J, Meseko CA, Adole JA, Tekki IS, Shamaki D, Hoffmann B (2019). Molecular characterization of Capripox viruses obtained from field outbreaks in Nigeria between 2000 and 2016. Transbound. Emerg. Dis. 66:1631–1641. https://doi.org/10.1111/tbed.13197
  • Ayelet G, Haftu R, Jemberie S, Belay A, Gelaye E, Sibhat B, Skjerve E, Asmare K (2014). Lumpy skin disease in cattle in central Ethiopia: outbreak investigation and isolation and molecular detection of the virus. Rev. Sci. Tech. Off. Int. Epiz.. 33 (3): 877-887. https://doi.org/10.20506/rst.33.3.2325
  • Beard PM, Sugar S, Bazarragchaa E, Gerelmaa U, Tserendorj Sh, Tuppurainen E, Sodnomdajaa R (2010). A description of two outbreaks of capripoxvirus disease in Mongolia. Vet. Microbiol. 142:427-431. https://doi.org/10.1016/j.vetmic.2009.10.018
  • Cao JX, Gershon PD, Black DN (1995). Sequence Analysis of HindIII Q2 Fragment of Capripoxvirus Reveals a Putative Gene Encoding a G-Protein-Coupled Chemokine Receptor Homologue. Virol. J. 209: 207-212. https://doi.org/10.1006/viro.1995.1244
  • El-Bagoury GF, El-Nahas EM, El-Habbaa AS, Abdel-Nour ET (2009). Isolation and Molecular Characterization of Lumpy Skin Disease Virus from Kalubia Governorate. Benha Vet. Med. J. 20:109 -119.
  • El-kady GHEM, El-Bagoury GF, El-Nahas EM, El-Habbaa AS (2018). A Phylogeny of Envelope protein like gene of Lumpy skin disease virus, Egypt, 2014. BVMJ. 34: 338-354. https://doi.org/10.21608/bvmj.2018.47868
  • El-Kenawy AA, El-Tholoth MS (2010). Sequence Analysis of Attachment Gene of Lumpy Skin Disease and Sheep Poxviruses. Virologica Sinica. 25 :409-416. https://doi.org/10.1007/s12250-010-3150-0
  • El-Khabaz KAS (2014). Rapid laboratory diagnosis of lumpy skin disease by using PCR technique. Assiut Vet. Med. J. 60: 37-41
  • El-Kholy AA, Soliman HMT, Abdelrahman KA (2008). Polymerase chain reaction for rapid diagnosis of a recent lumpy skin disease virus incursion to Egypt. Arab J. Biotech. 11:293-302
  • Gari G, Waret-Szkuta A, Grosbois V, Jacquiet P, Roger F (2010). Risk factors associated with observed clinical lumpy skin disease in Ethiopia. Epidemiol. Infect. 138: 1657-1666. https://doi.org/10.1017/S0950268810000506
  • Hodhod A, Elgendy E, Abd El-Moniem MI, Ibrahim MS (2020). Isolation and molecular characterization of lumpy skin disease virus in Egypt during 2017-2018. EJPMR 1:96-103.
  • House JA, Wilson TM, El Nakashly S, Karim IA, Ismail I, El Danaf N, Moussa AM, Ayoub NN (1990). The isolation of lumpy skin disease virus and bovine herpesvirus from cattle in Egypt. J. Vet. Diagn. Invest., 2:111-115. https://doi.org/10.1177/104063879000200205
  • Jackson PGG, Cockcroft PD (2002). clinical examination of farm animals. A textbook 1st edi. Chapter 2. Blackwell Science Ltd, Oxford OX20EL, UK. https://doi.org/10.1002/9780470752425
  • Kasem S, Saleh M, Qasim I, Hashim O, Alkarar A, Abu-Obeida A, Gaafer A, Hussien R, AL-Sahaf A, Al-Doweriej A, Bayoumi F, Hodhood A, Abdelatif M (2018). Outbreak investigation and molecular diagnosis of Lumpy skin disease among livestock in Saudi Arabia 2016. Transbound. Emerg. Dis. 65:e494–e500. https://doi.org/10.1111/tbed.12769
  • Kononov A, Byadovskaya O, Kononova S, Yashin R, Zinyakov N, Mischenko V, Perevozchikova N, Sprygin A (2019). Detection of vaccine-like strains of lumpy skin disease virus in outbreaks in Russia in 2017. Arch. Virol. 164:1575–1585. https://doi.org/10.1007/s00705-019-04229-6
  • Kostenis E (2004). A glance at G-protein-coupled receptors for lipid mediators: a growing receptor family with remarkably diverse ligands. Pharmacol. Ther. 102: 243–257. https://doi.org/10.1016/j.pharmthera.2004.04.005
  • Kumar S, Stecher G, Li M, Knyaz C, Tamura K (2018). Molecular Evolutionary Genetics Analysis across computing platforms. Mol. Biol. Evol. 35:1547-1549. https://doi.org/10.1093/molbev/msy096
  • Le Goff C, Fakhfakh E, Chadeyras A, Aba Adulugba E, Libeau G, Hammami S, Diallo A, Albina E (2005). host-range phylogenetic grouping of capripoxviruses. Genetic typing of CaPVs. In: Applications of Gene-Based Technologies for Improving Animal Production and Health in Developing Countries, pp 727-733. Edited by HPS Makkar, GJ Viljoen (eds.), © 2005 IAEA. Printed in the Netherlands. https://doi.org/10.1007/1-4020-3312-5_58
  • Le Goff C, Lamien CE, Fakhfakh E, Chadeyras A, Aba-Adulugba E, Libeau G, Tuppurainen E, Wallace DB, Adam T, Silber R, Gulyaz V, Madani H, Caufour P, Hammami S, Diallo A, Albina E (2009). Capripoxvirus G-protein-coupled chemokine receptor: a host-range gene suitable for virus animal origin discrimination. J. Gen. Virol. 90 :1967-1977. https://doi.org/10.1099/vir.0.010686-0
  • Ochwo S, VanderWaal K, Ndekezi C, Nkamwesiga J, Munsey A, Witto SG, Nantima N, Mayanja F, Okurut ARA, Atuhaire DK, Mwiine FN (2020). Molecular detection and phylogenetic analysis of lumpy skin disease virus from outbreaks in Uganda 2017–2018. BMC Vet. Res. 16:66. https://doi.org/10.1186/s12917-020-02288-5
  • OIE (2018). lumpy skin disease. In OIE terrestrial manual. Paris: World Organisation for Animal Health. Chapter 3.4.12: 1158-1171.
  • Raof AMA, Saad AA, Ahmed LA (2010). Use of a polymerase chain reaction assay to detect lumpy skin disease virus (LSDV) in skin lesions of cattle. Epypt. J. Comp. Path. 23(1): 116-126.
  • Rashid PMA, Baba Sheikh MO, Raheem ZH, Marouf AS (2017). Molecular characterisation of lumpy skin disease virus and sheeppox virus based on P32 gene. Bulg.J. Vet. Med.20(2):131-140. https://doi.org/10.15547/bjvm.984
  • Rouby SR, Bazid A, Wasfy M, El-Sayed M (2019). Capripoxviruses: Exploring the genetic relatedness between field and vaccine strains from Egypt. Vet. World. 12(12): 1924-1930. https://doi.org/10.14202/vetworld.2019.1924-1930
  • Sameea Yousefi P, Dalir-Naghadeh B, Mardani K, Jalilzadeh-Amin G (2018). Phylogenetic analysis of the lumpy skin disease viruses in northwest of Iran. Trop. Anim. Health Prod. 50:1851-1858. https://doi.org/10.1007/s11250-018-1634-3
  • Sharawi SSA, Abd El-Rahim IHA (2011). The utility of polymerase chain reaction for diagnosis of lumpy skin disease in cattle and water buffaloes in Egypt. Rev. Sci. Tech. Off. Int. Epiz. 30(3): 821-830. https://doi.org/10.20506/rst.30.3.2075
  • Tamam SM, El-shereif NM, Shokier KA (2018). Serological and molecular characterization of recent lumpy skin disease virus isolates from naturally infected previously vaccinated cattle in Egypt. J. Vet. Med. Res. 25:80-91. https://doi.org/10.21608/jvmr.2018.43304
  • Tasioudi KE, Antoniou SE, Iliadou P, Sachpatzidis A, Plevraki E, Agianniotaki EI, Fouki C, Mangana-Vougiouka O, Chondrokouki E, Dile C (2016). Emergence of lumpy skin disease in Greece, 2015. Transbound. Emerg. Dis. 63: 260-265. https://doi.org/10.1111/tbed.12497
  • Tassew A, Assefa A, Gelaye E, Bayisa B, Ftiwi M (2018). Identification and molecular characterization of lumpy skin disease virus in east Hararghe and east Shoa zone, Oromia regional state. ARC J. Anim. Vet. Sci. 4: 1-16
  • Toplak I, Petrović T, Vidanović D, Lazić S, Šekler M, Manić M, Petrović M, Kuhar U (2017). Complete genome sequence of lumpy skin disease virus isolate SERBIA/Bujanovac/2016, detected during an outbreak in the Balkan area. Genome Announc. 5(35): e00882-17. https://doi.org/10.1128/genomeA.00882-17
  • Tulman ER, Afonso CL, Lu Z, Zsak L, Sur JH, Sandybaev NT, Kerembekova UZ, Zaitsev VL, Kutish GF, Rock DL (2002). The genomes of sheeppox and goatpox viruses. J. Virol. 76: 6054–6061. https://doi.org/10.1128/JVI.76.12.6054-6061.2002
  • Tuppurainen ESM, Venter EH, Coetzer JAW (2005). The detection of lumpy skin disease virus in samples of experimentally infected cattle using different diagnostic techniques. Onderstepoort J. Vet. Res. 72:153-164. https://doi.org/10.4102/ojvr.v72i2.213
  • Zeynalova S, Asadov K, Guliyev F, Vatani M, Aliyev V (2016). Epizootology and molecular diagnosis of lumpy skin disease among livestock in Azerbaijan. Front. Microbiol. 7: 1022. https://doi.org/10.3389/fmicb.2016.01022
  •