Antisense knockdown of pyruvate dehydrogenase kinase promotes the neutral lipid accumulation in the diatom Phaeodactylum tricornutum
© Ma et al.; licensee BioMed Central Ltd. 2014
Received: 22 May 2014
Accepted: 3 July 2014
Published: 9 August 2014
Microalgae have been an emerging biofuel resource; however, the germplasm improvement has been slow due to the lack of molecular tools. Pyruvate dehydrogenase kinase (PDK) deactivates the pyruvate dehydrogenase complex (PDC) which catalyzes the oxidative decarboxylation of pyruvate. Acetyl-CoA production via PDC is important in plant tissues that are active in fatty acid synthesis.
A 1261-bp cDNA of a putative PDK gene (PtPDK) was cloned from a diatom Phaeodactylum tricornutum, and PtPDK antisense knockdown transgenic diatoms were generated. Both PtPDK transcript abundance and enzyme activity were reduced significantly due to antisense knockdown of PtPDK. Neutral lipid content of transgenic diatom cells increased up to 82% as determined by Nile red staining, and fatty acid composition was not altered. Transgenic cells showed slightly lower growth rate but similar cell size with the wild type, hence retaining similar biomass productivity.
This work first obtained a successful engineered diatom regulating a key gene involved in lipid metabolism. Our findings also provide powerful indications in enhancing microalgal lipid production by metabolic engineering for biofuel industry.
KeywordsMicroalga Pyruvate dehydrogenase kinase Antisense Lipid Biofuel
Diatoms are the main component of marine phytoplankton and play an important role in providing marine primary productivity. Particularly, a range of novel metabolic pathways have been found in diatoms which were presumably acquired during their evolution , so they are of major interest for the discovery of novel biological processes which are not present in other more intensively studied model organisms . Microalgae have become an emerging biofuel resource . Many diatom species are good candidates for biofuel production due to their high lipid content and rapid growth, such as Phaeodactylum tricornutum. Genetic transformation has been established in few diatoms including Thalassiosira pseudonana, P. tricornutum,. Thus, it is possible to genetically manipulate P. tricornutum for improving traits suited to biodiesel production .
The mitochondrial pyruvate dehydrogenase complex (PDC) catalyzes the irreversible oxidative decarboxylation of pyruvate to acetyl-CoA, thus directing the carbon flow into TCA cycle for yielding NADH to satisfy the energy demands of cells . Thus PDC activity indicates the rate of the entry of pyruvate into the TCA cycle and subsequently to oxidative phosphorylation or fatty acid synthesis. The mitochondrial PDC, unlike its plastid counterpart, has an associated pyruvate dehydrogenase kinase (PDK) . The PDC activity is primarily regulated by reversible phosphorylation by PDK, which phosphorylates and deactivates PDC .
In terms of the critical role of PDK in regulating PDC activity, it is important to unravel the molecular mechanisms behind the regulation. Comparative homology modeling of PDK isozymes from Xenopus tropicalis revealed that acetyl-CoA production via PDC is important in tissues that are active in fatty acid synthesis . Overexpression of Brassica napus PDK in Arabidopsis repressed the PDC activity, leading to a decrease in seed oil content . These results propose the possibility of manipulating PDK for promoting fatty acid synthesis. In the current study, we characterized a PDK gene from P. tricornutum and investigated the effect of antisense knockdown of PDK in P. tricornutum.
Prediction and sequence analysis of putative PtPDK
Molecular characterization of transgenic diatom by PCR, qPCR and enzyme activity assay
The transcript abundance of PtPDK in transgenic diatoms in the stationary phase was determined by qPCR. As shown in Figure 2B, the transgenic lines 1 and 2 exhibited a significant decrease in transcription by 65% and 42%, respectively. It implicates that the antisense knockdown had successfully repressed the mRNA expression levels of the PtPDK gene, while still retained some basal transcription.
The enzyme activity of PDK was further determined. As shown in Figure 2C, the relative enzyme activities of transgenic lines 1 and 2 declined by 55% and 23%, respectively, compared with that of the wild type. This result demonstrated that the suppressed PtPDK expression consequently resulted in a declined enzyme activity and transcript abundance showed a high correlation with enzyme activity.
Biomass and lipid productivities of transgenic diatoms
The neutral lipid content per cell indicated by the relative fluorescence intensity in Nile red-stained cultures in the stationary phase showed an increase of 82% and 33% in cells of transgenic line 1 and 2, respectively (Figure 3B). The neutral lipid content in wild type cells was initially determined to be 23.1% in dry weight, thereby the neutral lipid content in transgenic line 1 and 2 in dry weight were calculated to be approximately 42.10% and 30.76%, respectively. Hence, the antisense knockdown of PtPDK led to a significant increase in lipid productivity in genetically improved P. tricornutum.
Fatty acid profiles were further analyzed in diatom cells in the stationary phase (Figure 3C). Only a minor difference in fatty acid composition was found between transgenic and wild type algal cells which fell in the range of the standard deviation, except for C18:0 which was apparently lower in transgenic lines. These results indicate that antisense knockdown of PDK did not have much impact on the fatty acid composition of P. tricornutum.
Morphological observation of transgenic diatom cells
The current work reported a promising trial on the genetic improvement of microalgal strains for biofuel production. The primary requirement in microalgal biofuel industry is to have excellent algal species/strains. The ideal strain should have a combination of characteristics including fast growth, high oil content, strong resistance, and be suitable for large-scale cultivation. Metabolic engineering is a promising approach to achieve such an ideal microalgal strain. It is only more recently that we have the molecular tools for metabolic engineering and investigation on the complexity of changes in microalgae at molecular level. Recently, Trentacoste et al. have tried genetic engineering of the diatom T. pseudonana to improve lipid content by knockdown of a lipase gene . The antisense strains showed an uncompromised growth comparably to the wild type, but RNAi strains showed a decreased growth. In our study, we have improved the neutral lipid contents in another diatom P. tricornutum by antisense knockdown of PDK, while the transgenic algae displayed an inconspicuously decreased growth comparably to the wild type.
Moreover, the distribution of cellular resources during the growth of microalgae is affected by the availability of light, temperature, nutrients and the amount of energy required for their use. Based on this point, for biofuel industry, the simultaneous production of biomass and oil can be addressed more rapidly by using a combination of molecular biology and environmental control in future studies.
Since de novo fatty acid synthesis occurs in the chloroplast, Roessler and Ohlrogge tried to express and target acetyl-CoA carboxylase (ACCase) to the chloroplast of the diatom Cyclotella cryptica. Although expression of ACCase was elevated, lipid yields did not exhibit significant changes . The role of mitochondrion is usually neglected in the study of lipid biosynthesis, however, it especially plays a role in affecting the lipid synthesis through energy conversion and carbon metabolism in the cell. In seeds of Arabidopsis thaliana, the plastidial pyruvate kinases were important for fatty acid synthesis . Although activity of PDC is directly inhibited by acetyl-CoA and NADH, which are products of aerobic oxidation of fatty acids, the major regulation on PDC is the covalent modification of PDC leading to the activated or inhibited activity . PDK, a member of serine/threonine kinase family, plays a role in regulating carbohydrate metabolism by deactivation of PDC via the reversible phosphorylation of pyruvate dehydrogenase (E1) α subunit to achieve ATP-dependent inactivation of pyruvate dehydrogenase . PDK has an effect on short-term regulation of pyruvate generated by glycolysis and its activity could be regulated by the molar ratios of acetyl-CoA/CoA, NADH/NAD .
In PDHK4 knockout mice, the promoted fat accumulation and reduced fatty acid oxidation resulted from activation of PDC by lower phosphorylation of PDC through repressing the expression of PDHK (pyruvate dehydrogenase kinase). The findings suggest PDHK4 deficiency could result in the altered upstream signaling components associated with the regulation of lipid metabolism . Results of Marillia et al. provided evidence that PDK may participate in the regulation of lipid biosynthesis in Arabidopsis thaliana. They hypothesized that enhanced mitochondrial PDC activity through antisense knockdown of PDK expression may lead to an increased availability of acetyl moieties converted from pyruvate which are used in the synthesis of storage lipids in the developing seed. This hypothesis is consistent with the metabolic fluxes determined in developing embryos of Brassica napus, in which mitochondrial flux largely contributed to the provision of precursors for cytosolic fatty acid elongation .
However, the study on PDC in fatty acid biosynthesis in microalgae has not been reported. Considering that complete blocking of gene expression might be harmful to the cell and antisense knockdown usually achieves partial blocking, thus, in this work, antisense knockdown of PDK was adopted to develop genetically engineered microalgae suitable for industrial application. In this study, we demonstrated that effective knockdown of PtPDK expression resulted in increased lipid content without compromising biomass much. Transgenic cells showed similar growth as well as cell size with the wild type, therefore ensuring the biomass productivity. Results also implicated that PDK plays an important role in storage lipid accumulation through the regulation of mitochondrial carbon flux. Though the mere step of antisense knockdown of PtPDK has not promoted overall lipid yield enough to achieve an “ideal” industrial strain, the trial here has paved the way for microalgal biofuel development through the approach of metabolic engineering.
Materials and methods
Microalga and culture conditions
Microalga strain Phaeodactylum tricornutum was obtained from the Freshwater Algae Culture Collection of the Institute of Hydrobiology, China (No. FACHB-863). Microalgae were grown as batch cultures in flasks containing f/2-Si medium (omitting Na2SiO3 · 9H2O). Seawater taken from the Gulf of Dayawan was filter-sterilized through 0.22 μm filters and amended with f/2 nutrients and used as medium. Cultures in liquid medium or on the plate were grown at 21 ± 1°C in an artificial climate incubator. Cool-white fluorescent tubes provided an irradiance of 200 μmol photons m−2 s−1 under a 12/12 h light/dark regime. Growth curve was determined by counting cells with Brightline hemocytometer under a microscope every day.
Cloning and analysis of putative PtPDKgene
The 1,261-bp full-length coding region of putative PtPDK was amplified from total RNA of P. tricornutum by reverse transcription PCR with primers: 5′-ACCATGGAATTCGGAAATGTGAAACC-3′ and 5′-GTAATTGGCACGGGTTCACTTGTA-3′ flanking the two ends of PtPDK cDNA. The amplified cDNA sequence was confirmed by sequencing analysis at both orientations. Amino acid sequence alignment and similarity among species were determined using ClustalX version 1.83 and BLAST on NCBI (http://blast.ncbi.nlm.nih.gov/Blast/). Phylogenetic tree of protein clusters from various species was constructed by neighbor-joining (NJ) method using software MEGA 5 . Subcellular localization of PtPDK was predicted by softwares SignalP 4.1 (http://www.cbs.dtu.dk/services/SignalP/) and iPSORT (http://ipsort.hgc.jp/predict.cgi) online.
PtPDK was cloned into a P. tricornutum transformation vector modified from a previously constructed plasmid pHY11 . In the final vector designated pHY18-antiPDK, PtPDK was cloned between fcpC promoter and fcpA terminator of fucoxanthin chlorophyll a/c binding protein gene of P. tricornutum. Chloramphenicol at the final concentration of 250 μg/ml was used as selection marker in transformed Escherichia coli and P. tricornutum cells. P. tricornutum was transformed with plasmid pHY18-antiPDK by electroporation using Bio-Rad GenePulser Xcell apparatus (Bio-Rad, USA) following the protocol of Niu et al. .
Molecular analysis of transgenic microalgae by genomic PCR and qPCR
The incorporation of constructed gene expression cassettes into diatom genomes was demonstrated by genomic DNA PCR. Genomic DNA was extracted using Universal Genomic DNA Extraction Kit Ver.3.0 (Takara). PCR was performed with primers of CAT gene including forward primer (5′-ATGGAGAAAAAAATCACTG-3′) and reverse primer (5′-TAAGCATTCTGCCGACAT-3′).
Transcriptional abundance of PtPDK was determined by quantitative real-time PCR (qPCR) performed in Boxin Co., Guangzhou, China. Total microalgae RNA was extracted from the microalage cells using an E.Z.N.A.™ Plant RNA Kit (Omega Bio-Tek) and reverse-transcribed using an Omniscript reverse transcription kit (QIAGEN) with random hexamer primers . Reactions were performed in 96-well optical reaction plates with 20 μl mixture per well, using a SYBR Green Kit (Takara) and a 7300 Sequence Detection System (Applied Biosystems) following the manufacturers’ instructions. Primers used for PtPDK were (5′-ACCATGGAATTCGGAAATGTGAAACC-3′) and (5′-GTAATTGGCACGGGTTCACTTGTA-3′) which are unique to PtPDK sequence. β-actin was used as a housekeeping marker control and primers were ACT1f (5′-AGGCAAAGCGTGGTGTTCTTA-3′) and ACT1r (5′-TCTGGGGAGCCTCAGTCAATA-3′). The Ct (threshold cycle) for each well was measured. PtPDK mRNA accumulation in transgenic and wild type cells were quantified after normalization to β-actin.
Neutral lipid content analysis
The dye Nile red has been used to detect the intracellular lipid droplets by fluorescence microscopy and quantified with a photometer and photomultiplier connected to the microscope . Nile Red (Sigma) staining was performed following Yang et al.  to detect the cellular neutral lipid contents of P. tricornutum. Thirty microliters of Nile red (0.1 mg/ml acetone solution) was added into a 3 ml aliquot of cell culture in triplicates, then mixed by rapid inversion and incubated in darkness for 20 min at room temperature. The stained cell cultures were transferred to a 96-well plate and cellular fluorescence intensity was quantitatively determined with a microplate reader Hitach F4600 (Hitach, Japan). The wavelengths used were 530 nm for excitation and 592 nm for emission. The relative fluorescence intensity values could provide quantitative comparison of neutral lipid contents between the cell cultures. For dry weight determination, neutral lipid content in P. tricornutum cells was initially measured by gravimetric method. 20 mg lyophilized microalgae were mixed with 2 ml methanol, 2 ml chloroform, and 1 ml 5% NaCl by vortex for 2 min. The mixture was centrifuged at 8000 g for 4 min and the chloroform layer was collected. The same extraction procedure was repeated three times and the combined extracts were dried by N2 flow. Then the lipid residue was dried in oven at 60°C and weighed to yield dry weight.
Fatty acid composition analysis
Total lipids were extracted according to the protocol from Yang et al.  which was modified from Lepage et al. . In particular, transgenic and wild type cells in the stationary phase (250 ml each) were harvested for lipid extraction. Fatty acid composition was determined based on fatty acid methyl ester (FAME) by gas chromatography–mass spectrometry (GC-MS) at Guangdong Institute of Microbiology, China. A 30 m × 0.25 mm × 0.25 μm DB-5 quartz capillary column was used. Fatty acids were identified with the equipped NBS spectrum library. Integrated peak areas were calculated by normalization to acquire the relative contents (percentage of weight).
Morphological observation of transgenic diatom cells
To visualize the morphological changes of transgenic cells and particularly oil bodies, Nile red (0.1 mg/ml in acetone) was added into cultured cells at a 1:100 ratio and placed in darkness for 10 min. Cells were observed under a laser-scanning confocal microscope LSM510META (Zeiss), with excitation wavelength of 525 nm and emission wavelength of 550–570 nm. Images were captured randomly from at least 20 visual fields per sample, and typical cells are presented here.
Measurement of enzyme activity of PDK
Total proteins were extracted from transgenic and wild type algal cells for enzyme activity analysis. Enzyme activity of PDK was measured by using a Pyruvate dehydrogenase kinase kit (Jiemei, China). Relative activities of PDK were represented by the production of NADH (μmol/min per mg protein from algal cells) measured by absorbance .
Nile red stained diatom cells assayed by flow cytometry
Flow cytometry is used as a rapid approach for qualitative and quantitative analysis of individual cells which can detect a variety of biochemical reactions accurately in cells , especially analysis of large number of individual cells by fluorescence staining . It has been used to measure the fluorescent properties of individual algae cells and the fluorescence intensity reflects the lipid content of stained cells . In this work, 3 ml aliquot of cell culture in the stationary phase in triplicate were stained with Nile red (0.1 mg/ml acetone solution) as above, and then analyzed using a FACS-Aria microflow cytometer (Becton Dickinson, USA). The wavelengths used were 530 nm for excitation and 592 nm for emission.
This work was supported by the National Science and Technology Support Program (2011BAD14B03), The Innovation Project of Department of Education of Guangdong (2013KJCX0024), Science and Technology Project of Guangdong (2010B030600005).
- Armbrust EV, Berges JA, Bowler C, Green BR, Martinez D, Putnam NH, Zhou S, Allen AE, Apt KE, Bechner M: The genome of the diatom Thalassiosira pseudonana: ecology, evolution, and metabolism. Science. 2004, 306: 79-86. 10.1126/science.1101156.View ArticleGoogle Scholar
- Siaut M, Heijde M, Mangogna M, Montsant A, Coesel S, Allen A, Manfredonia A, Falciatore A, Bowler C: Molecular toolbox for studying diatom biology in Phaeodactylum tricornutum. Gene. 2007, 406: 23-35. 10.1016/j.gene.2007.05.022.View ArticleGoogle Scholar
- Hu Q, Sommerfeld M, Jarvis E, Ghirardi M, Posewitz M, Seibert M, Darzins A: Microalgal triacylglycerols as feedstocks for biofuel production: perspectives and advances. Plant J. 2008, 54: 621-639. 10.1111/j.1365-313X.2008.03492.x.View ArticleGoogle Scholar
- Poulsen N, Chesley PM, Kröger N: Molecular genetic manipulation of the diatom Thalassiosira pseudonana (Bacillariophyceae). J Phycol. 2006, 42: 1059-1065. 10.1111/j.1529-8817.2006.00269.x.View ArticleGoogle Scholar
- Miyagawa A, Okami T, Kira N, Yamaguchi H, Ohnishi K, Adachi M: Research note: high efficiency transformation of the diatom Phaeodactylum tricornutum with a promoter from the diatom Cylindrotheca fusiformis. Phycol Res. 2009, 57: 142-146. 10.1111/j.1440-1835.2009.00531.x. 9.12,View ArticleGoogle Scholar
- Igamberdiev AU, Lernmark U, Gardeström P: Activity of the mitochondrial pyruvate dehydrogenase complex in plants is stimulated in the presence of malate.Mitochondrion 2014, ᅟ:ᅟ. doi:10.1016/j.mito.2014.04.006.,Google Scholar
- Thelen J, Miernyk J, Randall D: Pyruvate dehydrogenase kinase from Arabidopsis thaliana: a protein histidine kinase that phosphorylates serine residues. Biochem J. 2000, 349: 195-201. 10.1042/0264-6021:3490195. 14.25View ArticleGoogle Scholar
- Sugden MC, Holness MJ: Recent advances in mechanisms regulating glucose oxidation at the level of the pyruvate dehydrogenase complex by PDKs. Amer J Physiol-Endocrinol Metab. 2003, 284: E855-E862.View ArticleGoogle Scholar
- Tokmakov AA: Comparative homology modeling of pyruvate dehydrogenase kinase isozymes from Xenopus tropicalis reveals structural basis for their subfunctionalization. J Mol Model. 2012, 18: 2567-2576. 10.1007/s00894-011-1281-3.View ArticleGoogle Scholar
- Li R-J, Hu Z-Y, Zhang H-S, Zhan G-M, Wang H-Z, Hua W: Cloning and functions analysis of a pyruvate dehydrogenase kinase in Brassica napus. Plant Cell Rep. 2011, 30: 1533-1540. 10.1007/s00299-011-1066-2.View ArticleGoogle Scholar
- Hyka P, Lickova S, Pribyl P, Melzoch K, Kovar K: Flow cytometry for the development of biotechnological processes with microalgae. Biotechnol Adv. 2013, 31: 2-16. 10.1016/j.biotechadv.2012.04.007.View ArticleGoogle Scholar
- Trentacoste EM, Shrestha RP, Smith SR, Glé C, Hartmann AC, Hildebrand M, Gerwick WH: Metabolic engineering of lipid catabolism increases microalgal lipid accumulation without compromising growth. Proc Natl Acad Sci. 2013, 110: 19748-19753. 10.1073/pnas.1309299110.View ArticleGoogle Scholar
- Roessler PG, Bleibaum JL, Thompson GA, Ohlrogge JB: Characteristics of the gene that encodes acetyl-CoA carboxylase in the diatom Cyclotella cryptica. Ann N Y Acad Sci. 1994, 721: 250-256. 10.1111/j.1749-6632.1994.tb47398.x.View ArticleGoogle Scholar
- Baud S, Wuillème S, Dubreucq B, De Almeida A, Vuagnat C, Lepiniec L, Miquel M, Rochat C: Function of plastidial pyruvate kinases in seeds of Arabidopsis thaliana. Plant J. 2007, 52: 405-419. 10.1111/j.1365-313X.2007.03232.x.View ArticleGoogle Scholar
- Kwon H-S, Harris RA: Mechanisms responsible for regulation of pyruvate dehydrogenase kinase 4 gene expression. Adv Enzyme Regul. 2004, 44: 109-122. 10.1016/j.advenzreg.2003.11.020.View ArticleGoogle Scholar
- Connaughton S, Chowdhury F, Attia RR, Song S, Zhang Y, Elam MB, Cook GA, Park EA: Regulation of pyruvate dehydrogenase kinase isoform 4 (PDK4) gene expression by glucocorticoids and insulin. Mol Cell Endocrinol. 2010, 315: 159-167. 10.1016/j.mce.2009.08.011.View ArticleGoogle Scholar
- Pettit FH, Pelley JW, Reed LJ: Regulation of pyruvate dehydrogenase kinase and phosphatase by acetyl-CoA/CoA and NADH/NAD ratios. Biochem Biophys Res Commun. 1975, 65: 575-582. 10.1016/S0006-291X(75)80185-9.View ArticleGoogle Scholar
- Hwang B, Jeoung N, Harris R: Pyruvate dehydrogenase kinase isoenzyme 4 (PDHK4) deficiency attenuates the long-term negative effects of a high-saturated fat diet. Biochem J. 2009, 423: 243-252. 10.1042/BJ20090390.View ArticleGoogle Scholar
- Schwender J, Shachar-Hill Y, Ohlrogge JB: Mitochondrial metabolism in developing embryos of Brassica napus. J Biol Chem. 2006, 281: 34040-34047. 10.1074/jbc.M606266200.View ArticleGoogle Scholar
- Tamura K, Dudley J, Nei M, Kumar S: MEGA4: molecular evolutionary genetics analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24: 1596-1599. 10.1093/molbev/msm092.View ArticleGoogle Scholar
- Niu YF, Yang ZK, Zhang MH, Zhu CC, Yang WD, Liu JS, Li HY: Transformation of diatom Phaeodactylum tricornutum by electroporation and establishment of inducible selection marker.Biotechniques 2012, doi:10.2144/000113881.Google Scholar
- Chomczynski P, Sacchi N: The single-step method of RNA isolation by acid guanidinium thiocyanate–phenol–chloroform extraction: twenty-something years on. Nat Protoc. 2006, 1: 581-585. 10.1038/nprot.2006.83.View ArticleGoogle Scholar
- Genicot G, Leroy J, Soom AV, Donnay I: The use of a fluorescent dye, Nile red, to evaluate the lipid content of single mammalian oocytes. Theriogenology. 2005, 63: 1181-1194. 10.1016/j.theriogenology.2004.06.006.View ArticleGoogle Scholar
- Yang ZK, Niu YF, Ma YH, Xue J, Zhang MH, Yang WD, Liu JS, Lu SH, Guan Y, Li HY: Molecular and cellular mechanisms of neutral lipid accumulation in diatom following nitrogen deprivation. Biotechnol Biofuels. 2013, 6: 67-10.1186/1754-6834-6-67.View ArticleGoogle Scholar
- Lepage G, Roy CC: Improved recovery of fatty acid through direct transesterification without prior extraction or purification. J Lipid Res. 1984, 25: 1391-1396.Google Scholar
- Zou J, Qi Q, Katavic V, Marillia E-F, Taylor DC: Effects of antisense repression of an Arabidopsis thaliana pyruvate dehydrogenase kinase cDNA on plant development. Plant Mol Biol. 1999, 41: 837-849. 10.1023/A:1006393726018. 17View ArticleGoogle Scholar
- Muirhead K, Horan P, Poste G: Flow cytometry: present and future. Nat Biotechnol. 1985, 3: 337-356. 10.1038/nbt0485-337.View ArticleGoogle Scholar
- Benoist C, Hacohen N: Flow cytometry, amped up. Science. 2011, 332: 677-678. 10.1126/science.1206351.View ArticleGoogle Scholar
- Franklin NM, Stauber JL, Lim RP: Development of multispecies algal bioassays using flow cytometry. Environ Toxicol Chem. 2004, 23: 1452-1462. 10.1897/03-250.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. 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.