Skip to main content
  • Research article
  • Open access
  • Published:

Karyotypic and mtDNA based characterization of Malaysian water buffalo

Abstract

Background

In Malaysia, the domestic water buffaloes (Bubalus bubalis) are classified into the swamp and the murrah buffaloes. Identification of these buffaloes is usually made via their phenotypic appearances. This study characterizes the subspecies of water buffaloes using karyotype, molecular and phylogenetic analyses. Blood of 105 buffaloes, phenotypically identified as swamp, murrah and crossbred buffaloes were cultured, terminated and harvested using conventional karyotype protocol to determine the number of chromosomes. Then, the D-loop of mitochondrial DNA of 10 swamp, 6 crossbred and 4 murrah buffaloes which were identified earlier by karyotyping were used to construct a phylogenetic tree was constructed.

Results

Karyotypic analysis confirmed that all 93 animals phenotypically identified as swamp buffaloes with 48 chromosomes, all 7 as crossbreds with 49 chromosomes, and all 5 as murrah buffaloes with 50 chromosomes. The D-loop of mitochondrial DNA analysis showed that 10 haplotypes were observed with haplotype diversity of 0.8000 ± 0.089. Sequence characterization revealed 72 variables sites in which 67 were parsimony informative sites with sequence diversity of 0.01906. The swamp and murrah buffaloes clearly formed 2 different clades in the phylogenetic tree, indicating clear maternal divergence from each other. The crossbreds were grouped within the swamp buffalo clade, indicating the dominant maternal swamp buffalo gene in the crossbreds.

Conclusion

Thus, the karyotyping could be used to differentiate the water buffaloes while genotypic analysis could be used to characterize the water buffaloes and their crossbreds.

Background

Water buffaloes, Bubalus bubalis are domesticated animals that play essential roles in agriculture, economy and food production. Water buffaloes are distributed widely around the world including the Indian subcontinent, Southeast Asia, China and across continents in Italy and Australia [1, 2]. In Malaysia, water buffaloes have great potential for meat and milk productions, especially in Sabah and a few states in Peninsular Malaysia [3,4,5].

The domestic water buffaloes are classified into two groups, based on their morphological and ecological characteristics. They are the swamp, Bubalus bubalis carabensis and the murrah, Bubalus bubalis bubalis buffaloes. In agriculture, swamp buffaloes are reared mainly for draught power in the paddy fields, oil palm and rubber plantations [6], while the murrah buffaloes produce good quality milk [7]. Morphologically, swamp buffaloes are small in size and tend to wallow in murky water such as swamp [8]. However, the murrah buffaloes are large in size with prominent horn and black jet in colour. In addition, murrah or “river” buffaloes prefer wallowing in clean water such as river, thus giving them the name [1]. Mostly, identification of the two subspecies of buffaloes is conducted based on conventional technique such as the morphological or physiological characteristics. Nevertheless, genetic predisposition, dominance and origin could not be determined using the conventional technique, particularly on crossbreds and effects of crossbreeding. Thus, the use of more accurate techniques are important to identify the phylogeny of the animals, or to compare the crossbred with their purebred parents [9, 10]. The techniques include karyotyping and molecular identification using molecular markers.

Karyotyping has been used to identify buffalo species while identification by molecular approach strengthen and supports the findings from the karyotyping. The mitochondrial DNA (mtDNA) is inherited maternally and has been influential for identifying ancestors, localizing domestication centres, tracking colonization and trading routes [11]. Most studies using mtDNA in buffaloes target the variable parts of the control region and the cytochrome b gene that look at population and phylogenetic relationships [12]. The D-loop mitochondrial DNA segment has been greatly emphasized because it acts as a crucial tool for assessing genetic relationships of individuals within species or different species. Thus, construction of phylogenetic tree should reveal the closeness of the subspecies of water buffalo. Therefore, the aims of this study were to identify the two subspecies of water buffaloes and their crossbred through karyotype and to use molecular approaches to characterize these buffaloes.

Methods

Animals

Water buffaloes (n = 105) were randomly selected from a government buffalo breeding farm in Sabah, Borneo and a private buffalo farm in Semenyih, Peninsular Malaysia. The swamp (n = 93) and crossbred (n = 7) buffaloes were selected from Sabah while the murrah (n = 5) buffaloes were selected from Semenyih, Selangor. The selection of each species was based on the phenotypic criteria. Both farms practiced natural breeding while the crossbreeding was done only in the government farm, involving breeding between the murrah males with swamp females.

The selected buffaloes were placed in paddocks for 3 days and were allowed to graze the Brachiaria grass (moisture: 6.56%, ash: 5.37%, crude fiber: 26.13%, crude fat: 5.11%, crude protein: 17.12%). Supplemented feed (moisture: 9.64%, ash: 5.44, crude fiber: 7.49, crude fat: 5.46%, crude protein: 18.15%) was provided at the rate of 1.5 kg/animal/day and drinking water was available ad libitum. The buffaloes were monitored daily and all animal handling and techniques were approved by the Institutional for Animal Care and Use Committee (IACUC), Universiti Putra Malaysia (UPM/IACUC/AUP-U017/2018).

Karyotyping analysis

Approximately 5–8 ml blood was drawn by venepuncture of the tail vein into a vacutainer tube containing heparin. The blood was immediately centrifuged at 1800 rpm for 10 min (Hettich, Zentrifugen D-7200 Tuttlingen) to produce a layer of buffy coat. Then, the buffy coat was pipetted slowly into the culture medium containing 8.0 ml of RPMI 1640, 2.0 ml of bovine calf serum, 0.1 ml of antibiotic, and 0.1 ml of mitogen (Pokeweed, PHA). The mixture was mixed gently and incubated at 37 °C for 72 h, while gently shaken twice daily to prevent leucocytes adhering to the walls of the culture flasks [13].

One hour before harvesting, 0.1 ml of colcemid was added into the culture bottle to stop the multiplying cells at metaphase. The mixture was gently shaken and incubated for another 1 h. At the end of the incubation period, the culture was gently shaken before the suspension was transferred into a 15 ml conical centrifuge tube. The suspension was spun for 5 min at 1800 rpm before the content was taken out leaving only 2.0 ml of the medium. Then, 6.0 ml of pre-warmed (37 °C) 0.075 M potassium chloride was introduced into each tube. The suspension was incubated at 37 °C C for 20 min and then centrifuged for 8 min at 1800 rpm. The supernatant was decanted with 6.0 ml fresh acetic methanol (1, glacial acetic acid: 3 methanol carnoy’s fixative) and was added into individual tube. The suspension was mixed well using pipette. The fixation step was repeated 3 times with the final suspension in 3 ml of carnoy’s fixative. The cell button was gently mixed with pasteur pipette. Two drops of 50% acetic acid were initially placed onto pre-clean slides, followed by another two drops before a drop of the cell suspension was added and spread. The slide was air-dried and another drop was put slightly away from the first. Next, Giemsa stain was used to cover the entire slide, washed with distilled water and air dried before viewing under microscope. The counting of chromosomes was performed on metaphase cells under the light microscope (Motic, China). Ten clearly observable spread of each sample were picked out and photographed (Motic Images Plus2.0).

Molecular analysis

Following karyotypic identification, 10 swamp, 6 crossbred and 4 murrah buffaloes were selected for genomic study. The DNA was extracted from the blood samples using DNeasy Blood and Tissue Kit (Qiagen) according to manufacturer’s protocol. To amplify the mitochondrial D-loop region, four primers (Table 1) were designed based on the known Bubalus bubalis mtDNA sequences [14], positioned in the conserved tRNA-Pro and tRNA-Phe. The volume of reaction mixture used for optimization of DNA was 50 μl that consisted of 25 μl of polymerase chain reaction (PCR) Master Mix (Promega GoTaq Green Master Mix), 0.4 pmol of both forward and reverse primers and ~ 60 ng/ μl of extracted DNA. Amplification was done using Applied Bioscience Thermocycler following PCR profile; preliminary denaturation at 95 °C for 3 mins, followed by 35 cycles of 94 °C for 30 s, 58 °C for 1 min and 72 °C for 1 min. This was followed by a final extension period of 72 °C for 4 mins. PCR products were electrophoresed on 1% agarose gels and purified using the Macherey-Nagel NucleoSpin Gel and PCR Clean-up according to the manufacturer’s instructions. Purified DNA was directly sequenced using Applied Biosystems BigDye Terminator v3.1 cycle sequencing kit chemistry for both forward and reverse strands.

Table 1 Oligonucleotides used for DNA amplification and sequencing of the bubaline mitochondrial D-loop region

The PCR products were used for sequencing (First Base, Malaysia). The sequences were edited and aligned using Geneious 11.1.12. All insertions/deletions in the alignment were used for multiple alignments. To demonstrate the phylogenetic clusters, the MEGA version 7 [15] was applied to construct maximum likelihood (ML) trees. Parameters used to construct the ML trees included using Hasegawa-Kishino-Yano model with 2000 bootstrapping replicates, Gamma distribution (G+) with 2 rate categories.

Results

The number and shape of the chromosomes of each sub-species buffalo are shown in Fig. 1. From the 105 selected buffaloes, all 93 were identified as swamp buffaloes with 48 chromosomes (Fig. 1a), 7 were crossbreds with 49 chromosomes (Fig. 1b) and 5 were murrah buffaloes with 50 chromosomes (Fig. 1c). The 48 chromosomes of swamp buffaloes consisted of 1 metacentric and 4 submetacentric chromosomes and the remaining were acrocentric, including the sex chromosomes. The crossbreds were presented with 49 chromosomes with 1 metacentric, 4 submetacentric and 18 acrocentric chromosomes. However, murrah buffaloes had 5 metacentric chromosomes and the remaining 20 pairs were acrocentric, including the sex chromosomes.

Fig. 1
figure 1

a Metaphase field exhibits 48 chromosomes showing female swamp buffalo (2n=48) by using giemsa staining method under 1000 magnification. b Metaphase field exhibits 49 chromosomes showing female crossbreed buffalo (2n=49) using giemsa staining method under 1000 magnification. c Metaphase field exhibits 50 chromosomes showing female murrah buffalo (2n=50) by using giemsa staining method under 1000 magnification

Polymerase chain reaction had successfully amplified the mitochondrial D-loop region of the buffalo samples. Both pairs of primers produced between ~ 500 to 750 bp of PCR products (Fig. 2a,b,c). The two sequences were aligned and combined to produce a 920 bp mitochondrial D-loop DNA. The mtDNA D-loop sequences had been deposited in the NCBI GenBank where the swamp (n = 10) accession numbers were MH746466-MH746475, the crossbreds (n = 6) accession numbers were MH746476-MH746481 and the murrah (n = 4) accession numbers were MH746482-MH746485. The D-loop regions of both swamp and crossbreds had slightly lower A + T content (58.28 ± 0.17% and 58.23 ± 0.15%) than the murrah buffaloes (59.48 ± 0.83%).

Fig. 2
figure 2

a Purified PCR product of primer H15773F and L16371R for swamp, crossbreed and murrah samples. The size is estimated around 500 bp. M represent 1 kb marker (Promega). b Purified PCR product of primer H16231F and L421R for murrah and crossbreed samples. The size is estimated around 500 bp. M represent 1 kb marker (Promega). c Purified PCR product of primer H16231F and L421R for swamp samples. The size is estimated around 500 bp. M represent 1 kb marker (Promega)

Comparing the DNA sequences of all 20 buffaloes (Table 2) revealed 10 different haplotypes with haplotype diversity of 0.8000 ± 0.089. The DNA sequence of the haplotype Swamp-3951 appeared with highest frequency (45%) and was shared by individuals from both swamp (n = 7) and crossbreds (n = 2). Within the murrah buffaloes, all four individuals revealed different haplotypes with haplotype diversity of 1.000 ± 0.177. The swamp buffaloes revealed three haplotypes with haplotype diversity of 0.5111 ± 0.164. Four haplotypes were observed among the crossbreds with haplotype diversity of 0.8667 ± 0.129.

Table 2 Comparison of the haplotype diversity (h), variable sites, parsimony informative sites and nucleotide diversity (Pi) among the swamp and crossbreed buffalo populations in Telupid and murrah buffalo population in Semenyih Farm, Selangor

Sequence characterization revealed 72 variable sites in which 67 were parsimony informative sites with sequence diversity of 0.01881. It appeared that swamp had the lowest sequence diversity of 0.00177, followed with the crossbreds (0.00269) and the highest in murrah with sequence diversity of 0.03811. Similarly, the swamp buffaloes showed lowest values for nucleotide and haplotype diversity while murrah buffaloes showed the highest values. This indicated that the murrah buffaloes used in this study were genetically more diversified than the swamp and crossbred buffaloes.

The phylogenetic tree revealed two distinct clades with extremely high posterior probability. Figure 3 shows the ability to distinguish between the swamp and the murrah buffaloes except for one murrah sample. The phylogenetic tree also revealed the fairly diversified murrah and much conserved swamp buffaloes. The 6 crossbreds were grouped in the same clade as swamp buffaloes, indicating the maternally descended from swamp buffalo.

Fig. 3
figure 3

The maximum likelihood phylogram reconstructed by MEGA 7 from 10 swamp, 4 murrah and 6 crossbreed buffaloes of the mitochondrial D-loop region, rooted by Bos taurus

Discussion

The current study shows that the swamp and murrah buffaloes differ by two pairs of chromosomes. This is due to a tandem fusion translocation of chromosomes 4 and 9 that were observed in the murrah buffalo into chromosome 1 in the swamp buffalo [16, 17]. The crossbreds had 49 chromosomes. These were in line with findings from previous studies that reported 48 chromosomes for swamp buffaloes [9, 18] and 49 chromosomes for crossbreds [19]. Thus, karyotyping appeared to confirm the phenotypic identification of the swamp, murrah and crossbred buffaloes.

The mtDNA has been broadly used in phylogenetic study since it is inherited directly from the mother, thus no change was made in the copy that makes it a suitable tool to be used to study the relationship of individual in the same species or even in different species. Besides that, the high mutation rate in the mtDNA provides the best estimation of the evolutionary relationships between individuals [20,21,22]. This study grouped the crossbreds within the same clade as the swamp buffaloes, indicating the strong swamp maternal genetic since the farm is practicing crossbred between murrah males and swamp females [4, 5]. Based on the maximum likelihood tree constructed in this study, it was found that the murrah was the ancestor of swamp buffaloes when Mur-1 was observed to be more closely related to the swamp group. This is in line with a study by Kiersten et al., (2004) [12] who reported that the two types of buffalo descended from a single domestication event some 5000 years ago. However, it contradicts the finding of Kumar et al., (2007) [23] who concluded an independent domestication of murrah and swamp buffaloes based on the sequences of mtDNA D-loop region and the cytochrome b analyses.

The murrah buffalo in this study showed high variability when all 4 samples were placed in different subclades, even though they were sampled from the same farm. On the other hand, the swamp and crossbred buffaloes showed less diversity. The murrah buffaloes were most probably originated from various sourced following importations from Indian subcontinent throughout the twentieth Century [8]. The swamp buffalo on the other hand, is the indigenous buffalo in Malaysia and many countries of the Southeast Asia with little importations [24, 25].

The crossbreds were the product of crossbreeding between male murrah and female swamp buffaloes, producing 49 chromosomes. Therefore, the crossbreds were not reproductively isolated because they shared the conserved regions of the genes from the parents. In fact, the crossbreds performed better than the swamp buffaloes especially in the meat and milk productions [7].

Conclusions

In conclusion, identification of the two buffalo sub-species and the crossbreds could be done phenotypically while confirmation was through karyotypic method, based on the number of chromosomes. In addition, analysis of mtDNA and the phylogenetic tree revealed that, the swamp buffaloes are genetically conserved and the crossbreds are dominantly swamp according to the maternal lineage using d-loop mtDNA. Thus, the combination of cytogenetic, phenotypic and mtDNA D-loop sequence analyses could be used to characterize the two buffalo subspecies and their crossbreds.

Abbreviations

DVS:

Department of Veterinary Services

G + :

Gamma distribution

IACUC:

Institutional for Animal Care and Use Committee

ML:

Maximum likelihood

mtDNA:

Mitochondrial DNA

Mur-1:

Murrah 1

PCR:

Polymerase Chain Reaction

References

  1. Borghese A. Buffalo production and research. FAO Ed. REU Technical Series, vol. 67; 2005. p. 1–315.

    Google Scholar 

  2. Kandeepan G, Biswas S, Rajkumar RS. Buffalo as potential food animal. Int J Livest Prod. 2009;1(1):1–5.

    Google Scholar 

  3. Ariff OM, Sharifah NY, Hafidz AW. Status of beef industry of Malaysia. Malays J Anim Sci. 2015;18(2):1–21.

    Google Scholar 

  4. Department of Veterinary Services, DVS. Livestock statistic. Malaysia: Department of Veterinary Services; 2016.

    Google Scholar 

  5. Department of Veterinary Services, DVS. Farm Policy. Malaysia: Department of Veterinary Services; 2016.

    Google Scholar 

  6. Riedel S, Meyer M, Schlecht E, Hülsebusch C, Schiborra A. Swamp buffalo keeping – an out-dated farming activity? A case study in smallholder farming systems in Xishuangbanna, Yunnan Province, PR China. J Agric Rural Dev Trop Subtrop. 2012;113(2):137–45.

    Google Scholar 

  7. Cruz LC. Recent developments in the buffalo industry of Asia. Ital J Vet. 2010;21:7–19.

    Google Scholar 

  8. Borghese A, Mazzi M. Buffalo population and strategies in the world. In: Buffalo production and research. Rome: FAO; 2005.

    Google Scholar 

  9. Fischer H, Ulbrich F. Chromosomes of the Murrah buffalo and its crossbreds with the Asiatic swamp buffalo (Bubalus bubalis). JAnim Breed Genet. 1967;84(1–4):110–4.

    Google Scholar 

  10. Udroiu I, Sgura A. Cytogenetic tests for animal production: state of the art and perspectives. Anim Genet. 2017;48(5):505–15.

    Article  CAS  Google Scholar 

  11. Bruford MW, Bradley DG, Luikart G. DNA markers reveal the complexity of livestock domestication. Nat Rev Genet. 2003;4:900–10.

    Article  CAS  Google Scholar 

  12. Kierstein G, Vallinoto M, Silva A, Schneider MP, Iannuzzi L, Brenig B. Analysis of mitochondrial D-loop region casts new light on domestic water buffalo (Bubalus bubalis) phylogeny. Mol Phylogenet Evol. 2004;30:308–24.

    Article  CAS  Google Scholar 

  13. Yimer N, Rosnina Y. Chromosomal anomalies and infertility in farm animals: a review. Pertanika J Trop Agric Sci. 2014;37(1):1–18.

    Google Scholar 

  14. Wang NS, Chen MR, Capodiferro T, Zhang H, Lancioni Y, Miao V, Chanthakhoun V, Wanapat M, Yindee Y, Zhang H, Lu L, Caporali R, Dang HY, Lan MX, Plath H, Chen JA, Lenstra A, Achilli CL. Whole mitogenomes reveal the history of swamp buffalo: initially shaped by glacial periods and eventually modelled by domestication. Sci Rep. 2017;7(4708):1–8.

    PubMed  PubMed Central  Google Scholar 

  15. Kumar S, Stecher G, Tamura K. MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets. Mol Biol Evol. 2016;33(7):1870–4.

    Article  CAS  Google Scholar 

  16. Chaudhary SC, Sharma M, Khatta V, Katoch BS. Utilization of some industrial by products in chick starter rations. Indian J Anim Sci. 1989;59(10):1313–8.

    Google Scholar 

  17. Yindee M, Vlamings BH, Wajjwalku W, Techakumphu M, Lohachit C, Sirivaidyapong S, Thitram C, Amarasinghe SS, Alexander PA, Colenbrander B, Lenstra JA. Y-chromosomal variation confirms independent domestications of swamp and buffalo. Anim Genet. 2010;41:433–5.

    CAS  PubMed  Google Scholar 

  18. Qiu H. Bovine breeds in China. Shanghai: Shanghai Scientific and Technical Publishing House; 1986.

    Google Scholar 

  19. Degrandi TM, Marques JRF, Gunski RJ, da Costa MRT, Marques LC, Vinade L, Garnero ADV. Cytogenetic identification of four generations of crossbred buffaloes maintained in a conservation program in the Marajo island/Brazil. J Biotechnol Biodiversity. 2014;5(2):162–71.

    Google Scholar 

  20. Lei C, Zhang W, Chen H, Lu F, Liu RY, Yang XY, Zhang HC, Liu ZG, Yao LB, Lu ZF, Zhao ZL. Two maternal lineages revealed by mitochondrial DNA D-loop sequences in Chinese native water buffaloes (Bubalus bubalis). Asian Australas J Anim Sci. 2007;20:471.

    Article  CAS  Google Scholar 

  21. Lei CZ, Zhang W, Chen H, et al. Independent maternal origin of Chinese swamp buffalo (Bubalus bubalis). Anim Genet. 2007;38:97–102.

    Article  CAS  Google Scholar 

  22. Burgstaller JP, Johnston IG, Jones NS, Albrechtová J, Kolbe T, Vogl C, Futschik A, Mayrhofer C, Klein D, Sabitzer S, Blattner M, Gülly C, Poulton J, Rülicke T, Piálek J, Steinborn R, Brem G. MtDNA segregation in heteroplasmic tissues is common in vivo and modulated by haplotype differences and developmental stage. Cell Rep. 2014;7(6):2031–41.

    Article  CAS  Google Scholar 

  23. Kumar S, Nagarajan M, Sandhu JS, Kumar N, Behl V. Phylogeography and domestication of Indian river buffalo. BMC Evol Biol. 2007;7:186.

    Article  Google Scholar 

  24. Lau CH, Drinkwater RD, Yusoff K, Tan SG, Hetzel DJS, Barker JS. Genetic diversity of Asian water buffalo (Bubalus bubalis): mitochondrial DNA D-loop and cytochrome b sequence variation. Anim Genet. 1998;29:253–64.

    Article  CAS  Google Scholar 

  25. Yue X-P, Li R, Xie W-M, Xu P, Chang T-C, Liu L, Cheng F, Zhang RF, Lan X-Y, Chen H, Lei C-Z. Phylogeography and domestication of chinese swamp buffalo. PLoS One. 2013;8(2):e56552.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

We are extremely grateful to all the authors and all the staff from Buffaloes Breeding Centre, Telupid Sabah who took part in this study.

Funding

This experiment was sponsored by the Flagship grant, Ministry of Science, Technology and Innovation, Malaysia (Project No: 6300858–10301). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. This publication is the work of the authors and Dr. Hafandi Ahmad will serve as guarantor for the contents of this paper.

Availability of data and materials

The datasets analysed during the current study are not publicly available due to study participant confidentiality but are available from the corresponding author on reasonable request.

Author information

Authors and Affiliations

Authors

Contributions

HA, MSS and MZS conceived and designed the experiments; HA, NALS, RY and MSS contributed to sampling. NALS, RY, MJE and LSS contributed to karyotyping and molecular analysis; NALS, MJE, and NIAG contributed to data analysis; NALS, MZS and HA drafted the manuscript; All authors read and approved the final manuscript.

Corresponding author

Correspondence to Hafandi Ahmad.

Ethics declarations

Ethics approval and consent to participate

All testing related to animal care and handling was approved by the Institutional Animal Care and Use Committee (IACUC) of the Faculty of Veterinary Medicine, Universiti Putra Malaysia, Malaysia (UPM/IACUC/AUP-U017/2018, 8 January 2018). The blood samples were collected from the farm by qualified veterinarians and with the permission of the owners of the animals.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Shaari, N.‘.A.L., Jaoi-Edward, M., Loo, S. et al. Karyotypic and mtDNA based characterization of Malaysian water buffalo. BMC Genet 20, 37 (2019). https://0-doi-org.brum.beds.ac.uk/10.1186/s12863-019-0741-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12863-019-0741-0

Keywords