Sequencing and validation of housekeeping genes for quantitative real-time PCR during the gonadotrophic cycle of Diploptera punctata

Identification of Potential Gene Targets for Suppressing Oviposition in Holotrichia parallela Using Comparative Transcriptome Analysis

Holotrichia parallela is an important plant pest. Comparative feeding experiments showed that the egg production, oviposition duration and survival rate of H. parallela beetles were significantly higher when they fed on elm leaves than when they fed on willow or purpus privet leaves. RNA sequencing was used to determine transcriptomic changes associated with oviposition. Comparative transcriptome analysis revealed that the beetles that fed on elm and willow had a total of 171 genes with differential expression. When the beetles fed on elm and purpus privet, 3568 genes had differential expression. The vitellogenesis, ovarian serine protease, odorant-binding proteins, acyl-CoA synthetase and follicle cell proteins were commonly upregulated genes in elm-fed beetles compared with those fed on willow/purpus privet leaves. The involvement of the follicle cell protein 3C gene in the regulation of oviposition was confirmed using RNA interference. The results provide insights into the molecular mechanisms underlying oviposition in H. parallela feeding on different host plants. This study also describes a method for identifying potentially effective genes for pest control.

A single point mutation in the Bemisia tabaci cytochrome-P450 CYP6CM1 causes enhanced resistance to neonicotinoids

The tobacco whitefly, Bemisia tabaci, is a polyphagous crop pest which causes high levels of economic damage across the globe. Insecticides are often required for the effective control of this species, among which the neonicotinoid class have been particularly widely used. Deciphering the mechanisms responsible for resistance to these chemicals is therefore critical to maintain control of B. tabaci and limit the damage it causes. An important mechanism of resistance to neonicotinoids in B. tabaci is the overexpression of the cytochrome P450 gene CYP6CM1 which leads to the enhanced detoxification of several neonicotinoids. In this study we show that qualitative changes in this P450 dramatically alter its metabolic capacity to detoxify neonicotinoids. CYP6CM1 was significantly over-expressed in two strains of B. tabaci which displayed differing levels of resistance to the neonicotinoids imidacloprid and thiamethoxam. Sequencing of the CYP6CM1 coding sequence from these strains revealed four different alleles encoding isoforms carrying several amino acid changes. Expression of these alleles in vitro and in vivo provided compelling evidence that a mutation (A387G), present in two of the CYP6CM1 alleles, results in enhanced resistance to several neonicotinoids. These data demonstrate the importance of both qualitative and quantitative changes in genes encoding detoxification enzymes in the evolution of insecticide resistance and have applied implications for resistance monitoring programs.

Reference genes for qRT-PCR normalisation in different tissues, developmental stages, and stress conditions of Hypericum perforatum

Hypericum perforatum L. is a widely known medicinal herb used mostly as a remedy for depression because it contains high levels of naphthodianthrones, phloroglucinols, alkaloids, and some other secondary metabolites. Quantitative real-time PCR (qRT-PCR) is an optimized method for the efficient and reliable quantification of gene expression studies. In general, reference genes are used in qRT-PCR analysis because of their known or suspected housekeeping roles. However, their expression level cannot be assumed to remain stable under all possible experimental conditions. Thus, the identification of high quality reference genes is essential for the interpretation of qRT-PCR data. In this study, we investigated the expression of 14 candidate genes, including nine housekeeping genes (HKGs) (ACT2, ACT3, ACT7, CYP1, EF1-α, GAPDH, TUB-α, TUB-β, and UBC2) and five potential candidate genes (GSA, PKS1, PP2A, RPL13, and SAND). Three programs-GeNorm, NormFinder, and BestKeeper-were applied to evaluate the gene expression stability across four different plant tissues, four developmental stages and a set of abiotic stress and hormonal treatments. Integrating all of the algorithms and evaluations revealed that ACT2 and TUB-β were the most stable combination in different developmental stages samples and all of the experimental samples. ACT2, TUB-β, and EF1-α were identified as the three most applicable reference genes in different tissues and stress-treated samples. The majority of the conventional HKGs performed better than the potential reference genes. The obtained results will aid in improving the credibility of the standardization and quantification of transcription levels in future expression studies on H. perforatum.

Evaluation of Reference Genes for Normalization of RT-qPCR Gene Expression Data for Trichoplusia ni Cells During Antheraea pernyi (Lepidoptera: Saturniidae) Multicapsid Nucleopolyhedrovirus (AnpeNPV) Infection

Baculovirus infection impacts global gene expression in the host cell, including the expression of housekeeping genes. Evaluation of candidate reference genes during a viral infection will inform the selection of appropriate reference gene(s) for the normalization of expression data generated by Reverse Transcription Quantitative Real-timePolymerase Chain Reaction (RT-qPCR). Antheraea pernyi multicapsid nucleopolyhedrovirus (AnpeNPV) is able to infect the High Five cells (Tn-Hi5). In the present study, 10 candidate reference genes were evaluated in AnpeNPV-infected Tn-Hi5 cells. Gene expression data were analyzed using geNorm, NormFinder, BestKeeper, and RefFinder. The candidate genes were further validated as reliable reference genes for RT-qPCR by analyzing the expression of three target genes. The results of data analysis using four statistical methods showed that RPS18 was the least stable among the candidate reference genes tested. 18S rRNA and 28S rRNA were not suitable as reference genes for RT-qPCR analysis in AnpeNPV-infected Tn-Hi5 cells. Comprehensive ranking of the 10 candidate reference genes by RefFinder analysis indicated that Ann B, c45128_g1, and ACT were the top three genes. Normalization of the expression of three target genes using the candidate reference genes indicated the same expression pattern when Ann B and c45128_g1 were used as reference genes, with slight differences in the relative expression at each infection time point. Overall, Ann B and c45128_g1 were recommended to be more suitable than the most commonly used reference genes, such as ACT, GAPDH, and TUB, for RT-qPCR data normalization in AnpeNPV-infected Tn-Hi5 cells up to 48 hpi.

Forkhead box transcription factor L2 activates Fcp3C to regulate insect chorion formation

Most animals are oviparous. However, the genes regulating egg shell formation remain not very clear. In this study, we found that Nilaparvata lugens Forkhead box transcription factor L2 (NlFoxL2) directly activated follicle cell protein 3C (NlFcp3C) to regulate chorion formation. NlFoxL2 and NlFcp3C had a similar expression pattern, both highly expressed in the follicular cells of female adults. Knockdown of NlFoxL2 or NlFcp3C also resulted in the same phenotypes: obesity and female infertility. RNA interference (RNAi) results suggested that NlFcp3C is a downstream gene of NlFoxL2. Furthermore, transient expression showed that NlFoxL2 could directly activate the NlFcp3C promoter. These results suggest that NlFcp3C is a direct target gene of NlFoxL2. Depletion of NlFoxL2 or NlFcp3C prevented normal chorion formation. Our results first revealed the functions of Fcp3C and FoxL2 in regulation of oocyte maturation in an oviparous animal.

Selection of Reference Genes for the Normalization of RT-qPCR Data in Gene Expression Studies in Insects: A Systematic Review

Reverse transcriptase-quantitative polymerase chain reaction (RT-qPCR) is a reliable technique for quantifying expression levels of targeted genes during various biological processes in numerous areas of clinical and biological research. Selection of appropriate reference genes for RT-qPCR normalization is an elementary prerequisite for reliable measurements of gene expression levels. Here, by analyzing datasets published between 2008 and 2017, we summarized the current trends in reference gene selection for insect gene expression studies that employed the most widely used SYBR Green method for RT-qPCR normalization. We curated 90 representative papers, mainly published in 2013–2017, in which a total of 78 insect species were investigated in 100 experiments. Furthermore, top five journals, top 10 frequently used reference genes, and top 10 experimental factors have been determined. The relationships between the numbers of the reference genes, experimental factors, analysis tools on the one hand and publication date (year) on the other hand was investigated by linear regression. We found that the more recently the paper was published, the more experimental factors it tended to explore, and more analysis tools it used. However, linear regression analysis did not reveal a significant correlation between the number of reference genes and the study publication date. Taken together, this meta-analysis will be of great help to researchers that plan gene expression studies in insects, especially the non-model ones, as it provides a summary of appropriate reference genes for expression studies, considers the optimal number of reference genes, and reviews the average number of experimental factors and analysis tools per study.

RNA interference in the Asian Longhorned Beetle:Identification of Key RNAi Genes and Reference Genes for RT-qPCR

Asian Longhorned Beetle (ALB) Anoplophora glabripennis is a serious invasive forest pest in several countries including the United States, Canada, and Europe. RNA interference (RNAi) technology is being developed as a novel method for pest management. Here, we identified the ALB core RNAi genes including those coding for Dicer, Argonaute, and double-stranded RNA-binding proteins (dsRBP) as well as for proteins involved in dsRNA transport and the systemic RNAi. We also compared expression of six potential reference genes that could be used to normalize gene expression and selected gapdh and rpl32 as the most reliable genes among different tissues and stages of ALB. Injection of double-stranded RNA (dsRNA) targeting gene coding for inhibitor of apoptosis (IAP) into larvae and adults resulted in a significant knockdown of this gene and caused the death of 90% of the larvae and 100% of adults. No mortality of both larvae and adults injected with dsRNA targeting gene coding for green fluorescence protein (GFP, as a negative control) was observed. These data suggest that functional RNAi machinery exists in ALB and a potential RNAi-based method could be developed for controlling this insect.

Molecular chaperone Jiv promotes the RNA replication of classical swine fever virus

The nonstructural protein 2 (NS2) of classical swine fever virus (CSFV) is a self-splicing ribozyme wherein the precursor protein NS2-3 is cleaved, and the cleavage efficiency of NS2-3 is crucial to the replication of viral RNA. However, the proteolytic activity of NS2 autoprotease may be achieved through a cellular chaperone called J-domain protein interacting with viral protein (Jiv) or its fragment Jiv90, as evidence suggests that Jiv is required for the proper functioning of the NS2 protein of bovine viral diarrhea virus. Hence, the expression of Jiv may be correlated with the replication efficiency of CSFV RNA. We investigated the expression levels of Jiv and viral RNA in CSFV-infected cells and tissues using Real-time RT-PCR or Western blot analysis. The obtained results show that Jiv90 possibly plays an important role in the lifecycle of CSFV because the distribution of Jiv90 protein shows a positive correlation with the viral load of CSFV. Furthermore, the overexpression or knockdown of Jiv90 in swine cells can also significantly promote or decrease the viral load, respectively. The detection of Flow cytometry shows that the overexpression of Jiv90 prolongs the G1 phase of cell cycles but has no effect on apoptosis. These findings are likely to be of benefit in clarifying the pathogenesis of the CSFV.

Identification of Reference Genes for Quantitative Gene Expression Studies in a Non-Model Tree Pistachio (Pistacia vera L.)

The tree species, Pistacia vera (P. vera) is an important commercial product that is salt-tolerant and long-lived, with a possible lifespan of over one thousand years. Gene expression analysis is an efficient method to explore the possible regulatory mechanisms underlying these characteristics. Therefore, having the most suitable set of reference genes is required for transcript level normalization under different conditions in P. vera. In the present study, we selected eight widely used reference genes, ACT, EF1α, α-TUB, β-TUB, GAPDH, CYP2, UBQ10, and 18S rRNA. Using qRT-PCR their expression was assessed in 54 different samples of three cultivars of P. vera. The samples were collected from different organs under various abiotic treatments (cold, drought, and salt) across three time points. Several statistical programs (geNorm, NormFinder, and BestKeeper) were applied to estimate the expression stability of candidate reference genes. Results obtained from the statistical analysis were then exposed to Rank aggregation package to generate a consensus gene rank. Based on our results, EF1α was found to be the superior reference gene in all samples under all abiotic treatments. In addition to EF1α, ACT and β-TUB were the second best reference genes for gene expression analysis in leaf and root. We recommended β-TUB as the second most stable gene for samples under the cold and drought treatments, while ACT holds the same position in samples analyzed under salt treatment. This report will benefit future research on the expression profiling of P. vera and other members of the Anacardiaceae family.

RXR/USP and EcR are critical for the regulation of reproduction and the control of JH biosynthesis in Diploptera punctata

During development and reproduction the response to ecdysteroids is mediated by a heterodimeric receptor complex comprising the retinoid X receptor/ultraspiracle (RXR/USP) and the ecdysone receptor (EcR). Here, the role of these receptors in the endocrine control of reproduction is examined in the cockroach Diploptera punctata. We report the sequence of four DpRXR and three DpEcR splice variants, including the first description of a Drosophila EcRB2-like isoform in a hemimetabolous insect. DpRXR and DpEcR are broadly expressed in the tissues of adult females, with relatively high transcript levels in the corpora allata (CA), nervous tissue and ovary. Developmental profiling revealed an inverse correlation between DpRXR and DpEcR expression and the activity of the CA. RNAi-mediated depletion of DpRXR and DpEcR did not affect oocyte growth, but inhibited oviposition and impaired chorion formation. Retained oocytes exhibited a degenerating follicular epithelium and were slowly resorbed. Treated animals showed significantly higher rates of JH biosynthesis and a decrease in ecdysteroid titers at the end of vitellogenesis. Reduction of DpRXR and DpEcR expression resulted in an upregulation of genes involved in JH production and a downregulation of allatostatin receptor mRNA in the CA. Treatment with dsRNA also affected the expression of genes downstream of JH in target tissues including vitellogenin and Krüppel-homolog 1 as well as Broad-Complex, an early ecdysone response gene. Overall, results suggest that DpRXR and DpEcR are not required early in the reproductive cycle when events are JH-dependent, but do mediate critical ecdysteroid feedback to the CA late in the gonadotropic cycle.

Reference gene selection for quantitative real-time PCR normalization in larvae of three species of Grapholitini (Lepidoptera: Tortricidae)

Despite the agricultural importance of species in the Grapholitini (Lepidoptera: Tortricidae), and the value of gene expression analysis for improved population management, few gene expression studies based on quantitative real-time PCR (qPCR) have been conducted for this tribe. Part of the reason for this lack of information is that suitable reference genes, which are fundamental for accurate normalization of qPCR studies, have not been identified for the tribe. Thus, the expression stability of six potential reference genes (ACT, AK, COI, EF1, ENO and TUB) was assessed in three different tissues (whole body, midgut and cuticle) of Cryptophlebia peltastica (Meyrick), Cydia pomonella (L.) and Thaumatotibia leucotreta (Meyrick). Additionally, these reference genes were tested using T. leucotreta at different temperatures (15°C, 25°C and 35°C) with and without baculovirus infection. Suitable reference genes were identified for the whole body and midgut tissue of all three species, and for cuticle tissue of Cy. pomonella and T. leucotreta. When T. leucotreta was infected with the virus at all temperature conditions ACT, AK and EF1 were found to be the most suitable reference genes for experimental normalization. In general, for all tissue types, species and stress conditions, AK and EF1 were the best-performing reference genes. However, even though the three species analysed were closely related and within the same tribe, each species required varying gene combinations for suitable normalization. This study provides the first reference gene evaluation for the Tortricidae, and paves the way for future qPCR analysis in Tortricidae.

Identification and de novo sequencing of housekeeping genes appropriate for gene expression analyses in farmed maraena whitefish (Coregonus maraena) during crowding stress

Maraena whitefish (Coregonus maraena; synonym Coregonus lavaretus f. balticus) is a high-quality food fish in the Southern Baltic Sea belonging to the group of salmonid fishes. Coregonus sp. is successfully kept in aquaculture throughout northern Europe (e.g. in Finland, Germany, Russia) and North America. In this regard, the molecular and immunological characterisation of stress response in maraena whitefish contributes to the development of robust and fast-growing maraena whitefish breeding strains for aquaculture. Thus, in the present study, the potential housekeeping genes beta actin (ACTB), elongation factor 1 alpha (EEF1A1), glyceraldehydes-3-phosphate dehydrogenase (GAPDH), ribosomal protein 9 (RPL9), ribosomal protein 32 (RPL32) and ribosomal protein S20 (RPS20) were de novo sequenced and tested concerning their applicability as reference genes in quantitative real-time PCR (qPCR) in maraena whitefish under different stocking densities. For this purpose, tissue samples of liver, kidney, gills, head kidney, skin, adipose tissue, heart and dorsal fin were investigated. qPCR data were analysed with Normfinder tool to determine gene expression stability. DNA sequencing exposed transcribed paralogous EEF1A1A and EEF1A1B genes differing in their putative protein structure. Normfinder analysis revealed RPL9 and RPL32 as most stable, GAPDH and ACTB as least stable genes for qPCR analyses, respectively. This is the first study that provides a subset of seven de novo sequenced housekeeping genes usable as reference genes in studies of stress response in maraena whitefish.

Identification and characterization of the NMDA receptor and its role in regulating reproduction in the cockroach Diploptera punctata

The NMDA receptor (NMDAR) plays important roles in excitatory neurotransmission and in the regulation of reproduction in mammals. NMDAR in insects comprises two subunits, NR1 and NR2. In this study, we identified two NR1 paralogs and eleven NR2 alternatively spliced variants in the cockroach Diploptera punctata. This is the first report of NR1 paralogs in insects. The tissue distributions and expression profiles of DpNR1A, DpNR1B and DpNR2 in different tissues were also investigated. Previous studies have demonstrated NMDA-stimulated biosynthesis of juvenile hormone (JH) in the corpora allata through the influx of extracellular Ca2+ in Diploptera punctata. However, our data show that the transcript levels of DpNR1A, DpNR1B and DpNR2 were low in the corpora allata. MK-801, a high-affinity antagonist of NMDAR, did not show any effect on JH biosynthesis in vitro. In addition, neither partial knockdown of DpNR2 nor in vivo treatment with a physiologically relevant dose of MK-801 resulted in any significant change in JH biosynthesis or basal oocyte growth. Injection of animals with a high dose of MK-801 (30 µg per animal per injection), which paralyzed the animals for 4–5 h, resulted in a significant decrease in JH biosynthesis on days 4 and 5. However, the reproductive events during the first gonadotrophic cycle in female D. punctata were unaffected. Thus, NMDAR does not appear to play important roles in the regulation of JH biosynthesis or mediate reproduction of female D. punctata.

Characterization of the juvenile hormone pathway in the viviparous cockroach, Diploptera punctata

Juvenile hormones (JHs) are key regulators of insect development and reproduction. The JH biosynthetic pathway is known to involve 13 discrete enzymatic steps. In the present study, we have characterized the JH biosynthetic pathway in the cockroach Diploptera punctata. The effect of exogenous JH precursors on JH biosynthesis was also determined. Based on sequence similarity, orthologs for the genes directly involved in the pathway were cloned, and their spatial and temporal transcript profiles were determined. The effect of shutting down the JH pathway in adult female cockroaches was studied by knocking down genes encoding HMG-CoA reductase (HMGR) and Juvenile hormone acid methyltransferase (JHAMT). As a result, oocyte development slowed as a consequence of reduction in JH biosynthesis. Oocyte length, fat body transcription of Vg and ovarian vitellin content significantly decreased. In addition, silencing HMGR and JHAMT resulted in a decrease in the transcript levels of other genes in the pathway.

Mode of action of allatostatins in the regulation of juvenile hormone biosynthesis in the cockroach, Diploptera punctata

The FGLamide allatostatins (FGL/ASTs) are a family of neuropeptides with pleiotropic functions, including the inhibition of juvenile hormone (JH) biosynthesis, vitellogenesis and muscle contraction. In the cockroach, Diploptera punctata, thirteen FGLa/ASTs and one allatostatin receptor (AstR) have been identified. However, the mode of action of ASTs in regulation of JH biosynthesis remains unclear. Here, we determined the tissue distribution of Dippu-AstR. And we expressed Dippu-AstR in vertebrate cell lines, and activated the receptor with the Dippu-ASTs. Our results show that all thirteen ASTs activated Dippu-AstR in a dose dependent manner, albeit with different potencies. Functional analysis of AstR in multiple cell lines demonstrated that activation of the AstR receptor resulted in elevated levels of Ca(2+) and cAMP, which suggests that Dippu-AstR can act through the Gαq and Gαs protein pathways. The study on the target of AST action reveals that FGL/AST affects JH biosynthesis prior to the entry of acetyl-CoA into the JH biosynthetic pathway.

Methoprene-tolerant (Met) knockdown in the adult female cockroach, Diploptera punctata completely inhibits ovarian development

Independent of the design of the life cycle of any insect, their growth and reproduction are highly choreographed through the action of two versatile hormones: ecdysteroids and juvenile hormones (JH). However, the means by which JH can target tissues and exert its pleiotropic physiological effects is currently still not completely elucidated. Although the identity of the one JH receptor is currently still elusive, recent evidence seems to point to the product of the Methoprene-tolerant gene (Met) as the most likely contender in transducing the action of JH. Studies on the role of this transcription factor have mostly been focused on immature insect stages. In this study we used the viviparous cockroach Diploptera punctata, a favorite model in studying JH endocrinology, to examine the role of Met during reproduction. A tissue distribution and developmental profile of transcript levels was determined for Met and its downstream partners during the first gonadotropic cycle of this cockroach. Using RNA interference, our study shows that silencing Met results in an arrest of basal oocyte development; vitellogenin is no longer transcribed in the fat body and no longer taken up by the ovary. Patency is not induced in these animals which fail to produce the characteristic profile of JH biosynthesis typical of the first gonadotropic cycle. Moreover, the ultrastructure of the follicle cells showed conspicuous whorls of rough endoplasmic reticulum and a failure to form chorion. Our study describes the role of Met on a cellular and physiological level during insect reproduction, and confirms the role of Met as a key factor in the JH signaling pathway.

Selection and assessment of reference genes for quantitative PCR normalization in migratory locust Locusta migratoria (Orthoptera: Acrididae)

Locusta migratoria is a classic hemimetamorphosis insect and has caused widespread economic damage to crops as a migratory pest. Researches on the expression pattern of functional genes in L. migratoria have drawn focus in recent years, especially with the release of genome information. Real-time quantitative PCR is the most reproducible and sensitive approach for detecting transcript expression levels of target genes, but optimal internal standards are key factors for its accuracy and reliability. Therefore, it's necessary to provide a systematic stability assessment of internal control for well-performed tests of target gene expression profile. In this study, twelve candidate genes (Ach, Act, Cht2, EF1α, RPL32, Hsp70, Tub, RP49, SDH, GAPDH, 18S, and His) were analyzed with four statistical methods: the delta Ct approach, geNorm, Bestkeeper and NormFinder. The results from these analyses aimed to choose the best suitable reference gene across different experimental situations for gene profile study in L. migratoria. The result demonstrated that for different developmental stages, EF1α, Hsp70 and RPL32 exhibited the most stable expression status for all samples; EF1α and RPL32 were selected as the best reference genes for studies involving embryo and larvae stages, while SDH and RP49 were identified for adult stage. The best-ranked reference genes across different tissues are RPL32, Hsp70 and RP49. For abiotic treatments, the most appropriate genes we identified were as follows: Act and SDH for larvae subjected to different insecticides; RPL32 and Ach for larvae exposed to different temperature treatments; and Act and Ach for larvae suffering from starvation. The present report should facilitate future researches on gene expression in L. migratoria with accessibly optimal reference genes under different experimental contexts.

Inhibitional effects of metal Zn²⁺ on the reproduction of Aphis medicaginis and its predation by Harmonia axyridis

BACKGROUND

Contamination, including metals, can disturb the reproductive processes of many organisms, including both prey and predatory insects. However, there is virtually no information on the effects of high level Zinc (Zn) pollution on aphids and ladybirds. The high concentrations of Zn²⁺ or Zn pollution inhibit reproduction in the phytophagous aphid, Aphis medicaginis, and the predatory ladybird Harmonia axyridis could provide important information.

RESULTS

It was observed in this study that Zn concentrations in Vicia faba (broad bean) seeds and seedlings in all Zn²⁺ treatments were significantly higher than that in the control group, and increased with increasing Zn²⁺ concentrations in the solution. The rate of reproduction in A. medicaginis declined significantly (p<0.05) over time in the five groups fed on broad bean seedlings treated with different concentrations of Zn²⁺ solution compared with the control group. These results showed that higher concentrations of Zn²⁺ significantly inhibited the reproductive capacity of A. medicaginis. We also cloned and identified a gene encoding vitellogenin (Vg) from A. medicaginis, which has an important role in vitellogenesis, and therefore, reproduction was affected by exposure to Zn²⁺. Expression of AmVg was reduced with increasing exposure to Zn²⁺ and also in the F1-F3 generations of aphids exposed to different Zn²⁺ concentrations. Predation by H. axyridis was also reduced in aphids exposed to high-levels of Zn²⁺. Similarly, ovipositioning by H. axyridis was also reduced.

CONCLUSIONS

Our results suggest that Zn²⁺ can significantly affect the reproductive capacity of both A. medicaginis and its predator H. axyridis, the former through effects on the expression of AmVg and the latter through avoidance of aphids containing high levels of Zn²⁺.