New Research In
Physical Sciences
Social Sciences
Featured Portals
Articles by Topic
Biological Sciences
Featured Portals
Articles by Topic
- Agricultural Sciences
- Anthropology
- Applied Biological Sciences
- Biochemistry
- Biophysics and Computational Biology
- Cell Biology
- Developmental Biology
- Ecology
- Environmental Sciences
- Evolution
- Genetics
- Immunology and Inflammation
- Medical Sciences
- Microbiology
- Neuroscience
- Pharmacology
- Physiology
- Plant Biology
- Population Biology
- Psychological and Cognitive Sciences
- Sustainability Science
- Systems Biology
Mechanics of mouse blastocyst hatching revealed by a hydrogel-based microdeformation assay
Edited by Janet Rossant, Hospital for Sick Children, University of Toronto, Toronto, ON, Canada, and approved August 15, 2018 (received for review November 15, 2017)

Significance
Blastocyst hatching is crucial for implantation of mammalian embryos and a common failure point during in vitro fertilization (IVF). We have little knowledge of the mechanical basis whereby an embryo hatches out of the zona pellucida. We have developed a technique to measure blastocyst pressure, allowing us to quantify physiological parameters and providing additional measures of efficiency in IVF optimization. We find that mechanical stretching of the zona by the blastocyst is essential for efficient hatching. Cryopreservation and thawing of embryos is common during IVF. Our technique reveals significant differences in microphysiology between fresh and thawed embryos. Our experimental and associated mathematical techniques are also applicable to other biological systems involving cavity formation, providing an approach for measuring forces in diverse contexts.
Abstract
Mammalian embryos are surrounded by an acellular shell, the zona pellucida. Hatching out of the zona is crucial for implantation and continued development of the embryo. Clinically, problems in hatching can contribute to failure in assisted reproductive intervention. Although hatching is fundamentally a mechanical process, due to limitations in methodology most studies focus on its biochemical properties. To understand the role of mechanical forces in hatching, we developed a hydrogel deformation-based method and analytical approach for measuring pressure in cyst-like tissues. Using this approach, we found that, in cultured blastocysts, pressure increased linearly, with intermittent falls. Inhibition of Na/K-ATPase led to a dosage-dependent reduction in blastocyst cavity pressure, consistent with its requirement for cavity formation. Reducing blastocyst pressure reduced the probability of hatching, highlighting the importance of mechanical forces in hatching. These measurements allowed us to infer details of microphysiology such as osmolarity, ion and water transport kinetics across the trophectoderm, and zona stiffness, allowing us to model the embryo as a thin-shell pressure vessel. We applied this technique to test whether cryopreservation, a process commonly used in assisted reproductive technology (ART), leads to alteration of the embryo and found that thawed embryos generated significantly lower pressure than fresh embryos, a previously unknown effect of cryopreservation. We show that reduced pressure is linked to delayed hatching. Our approach can be used to optimize in vitro fertilization (IVF) using precise measurement of embryo microphysiology. It is also applicable to other biological systems involving cavity formation, providing an approach for measuring forces in diverse contexts.
- mammalian embryo hatching
- embryo cryopreservation
- pressure measurement
- hydrogel deformation
- mouse blastocyst
Blastocyst formation underpins the foundation of early mammalian development, occurring before embryo implantation into the uterus. During this process, the embryo creates a pressurized cavity within it and eventually breaks out of the protective acellular outer layer, the zona pellucida. Without hatching from the zona pellucida, embryos are unable to implant into the uterus, which therefore represents a stumbling block for in vitro fertilization (IVF) technologies (1). In addition to driving molecular changes to the trophoblast cells and the uterus (2), the timing of embryo hatching is also important, with early hatching resulting in ectopic pregnancies (3) and late hatching causing the window of receptivity to be missed (4). To date, the process has been difficult to study from the mechanobiology perspective due to the lack of suitable tools for imposing and characterizing deformations in 3D tissues. Moreover, embryo assessment is a crucial step during IVF procedures, particularly for evaluating the quality of cryopreserved embryos that have been thawed (5). However, existing methods rely on qualitative scoring, and a well-characterized quantitative methodology could provide more reliable predictive accuracy.
Previous studies have calculated tissue forces by directly observing the deformation of structures of known stiffness, for example by incorporating deformable oil droplets into developing tissues to measure the anisotropic (shear) stress (6), or by culturing a layer of cells atop a 2D array of bendable hydrogel pillars (7) to measure traction. Atomic force microscopy has also been used to measure local mechanical properties of individual blastocyst cells (8) and the hydrogel shell (9⇓⇓–12). However, integrating the information into a systematic model has proven difficult due to a lack of embryonic force measurements. It is understood how embryonic tissue would respond to forces, but it is unknown what the forces are.
Mammalian preimplantation embryos have been intensely studied from the molecular biology perspective. These studies revealed much detail about the genetic regulation of cell fate decisions (13, 14) and the molecular mechanism of blastocyst cavity development. It has been shown that cavity expansion is driven by active ion transport (15), mediated by Na/K-ATPases in the outer epithelial cell layer sealed by tight junctions (16, 17). Increasing internal osmolarity then results in water transport along the osmotic gradient, which is mediated by aquaporin water channels (18, 19). These components act together to create the system responsible for internal pressurization of the cavity and eventual embryo hatching from its outer shell. The zona pellucida is composed of several glycoproteins, which form an elastic hydrogel layer around the embryo. The width of this layer decreases during embryo expansion as a precursor to hatching, partly in mechanical response to the tangential stretching imposed by the embryo’s increasing volume, and partly due to the action of degrading enzymes (20⇓–22). Although both mechanical and enzymatic factors have been suggested to be important, making a distinction between them has proven difficult due to the lack of experimental techniques.
Here, we therefore extend the idea of measuring mechanical forces using deformable media and develop an approach focused on confining embryos in hydrogels. This allows us to measure not only the pressure but also the dynamics of other important underlying processes, such as ion and water transport, and quantify the mechanical components of the hatching process. We use this approach to characterize the microphysiology of embryos and discover previously unknown blastocyst cavity pressure differences between freshly collected embryos and embryos thawed after cryopreservation.
Results and Discussion
Relating Hydrogel Deformation and Tissue Pressure.
Developing E3.5 mouse embryos were embedded in microfabricated cylindrical cavities within hydrogels of known Young’s modulus (E). Optical bright-field microscopy (Fig. 1A and Movies S1 and S2) was then used to measure channel dilation by the embryo from a (undeformed) to a + h (maximum dilated radius; SI Appendix, Fig. S1A). A simple elastic model of the channel then allowed us to calculate the pressure exerted and hence the pressure P within the embryo as follows (Eq. 1: Calculating pressure exerted upon the hydrogel):
Summary of embryo development in hydrogel cavities. Pressures in developing embryos were measured by using deformable hydrogels (Movies S1 and S2) (A). The embryos were cultured from early cavity stage ∼3.5 d after fertilization in KSOMaa Evolve media of osmolarity of 0.252 Osm. (Scale bar: 50 µm.) The elastic hydrogel made it possible to calculate the pressure (B), which was generated during volumetric expansion (C). The time lapse was taken over a period of 2 h with a resolution of 6 min.
To test whether culturing embryos within the hydrogel impairs viability, we analyzed the number of inner cell mass and trophectoderm (TE) cells, as well as the number of fragmented nuclei. Across various sets of experimental conditions, we found no significant differences between embryos cultured within hydrogels and control unconfined embryos, suggesting that our method was not harmful (SI Appendix, Fig. S1 C–E).
The Effects of Embryo Compression.
Embryo compression itself is likely to generate some additional internal pressure based on elongating an embryo from the spherical shape. This requires the zona pellucida to stretch, which is an additional source of pressure. Therefore, internal embryo pressure is the sum of measured pressure acting on the hydrogel and the added pressure due to deformation from spherical to elongated shape. To estimate the pressure originating solely from the deformation, we consider that, on average, embryos increase pressure with increasing surface area at the rate of 0.03 Pa/µm2. This can be estimated by tracking embryo area and pressure increases with time. It can also be estimated that compression results in surface area increase of less than 3,000 µm2, which indicates that internal pressure may be higher than measured by as much as 90 Pa. Since this is below the measurement error of 150 Pa, we can assume that pressure exerted on the hydrogel represents the internal embryo pressure. Also, the measurement error of 150 Pa made it impossible to measure zona pellucida free and hatched embryos, which could no longer significantly deform the hydrogel (SI Appendix, Fig. S2). Intact outer shell was required for generating a significant, measurable pressure.
Microphysiology of Blastocyst Cavity Development.
To build a mathematical model for embryo hatching, it is crucial to look at the details of how all of the different components of the system create a pressurized cavity. By measuring blastocyst pressure and volume dynamics (Fig. 1 B and C), we could see that pressure and volume rises were linear over the course of development, but interrupted by sudden drops, which likely corresponded to TE cell divisions or epithelial cell junction cycling (SI Appendix, Fig. S3). Overall, developing mouse blastocyst-stage embryos generated pressures of around 1.4 kPa and had an average volume of 0.66 × 106 μm3 (Fig. 2 A and B). These pressure and volume measurements allowed us to derive other important physiological parameters, such as blastocyst osmolarity using the van’t Hoff (ideal gas) law (Eq. 2), where R is the ideal gas constant, V is cavity volume, and cmedia is 0.252 Osm:
Blastocyst pressure and volume measurements. (A) Pressure distribution in developing mouse embryos (n = 41). The small difference between the mean and maxima reflects the fact that embryos spend most of the time close to their maximal pressure. (B) Volume distribution in mouse embryos (n = 41). The average values were generated by measuring the maximum and mean pressures each embryo generates during development and then averaging these values over the normally developing embryos used for the experiments. (C) The generated pressure could be disrupted by using a Na/K-ATPase inhibitor ouabain. Embryos responded in a dose-dependent manner, which highlighted the importance of ion transport for blastocyst cavity expansion. The 100 and 200 µM inhibitor concentrations were achieved by diluting a 100 mg/mL DMSO stock; controls contained an equivalent amount of DMSO. At least seven embryos in each condition; P values are represented by the following: n.s., not significant; **, <0.01. (D) Increasing media osmolarity reduces blastocyst cavity pressure because of the reduced ion influx into the cavity. The figure demonstrates that the rate of ion transport across the TE layer was dependent on the osmolarity balance between the outside and inside of the blastocyst cavity. At least five embryos were used per condition.
The rate of pressure and osmolarity change also makes it possible to estimate the ion flux and the number of involved Na/K-ATPases. From previous studies (15), it is known that the ion gradient is primarily generated by the activity of Na/K-ATPase enzymes, which use energy to transport ions across an ion gradient. To verify that pressure generation was also actively driven, we used ouabain to inhibit ion transport (Fig. 2C), which indeed resulted in reduced cavity pressure in a dose-dependent manner. Additionally, ion transport was also found to be dependent on media osmolarity, which indicates that the process is actively regulated (Fig. 2D).
To estimate the number of Na/K-ATPases present on the trophoblast cells, we used published data reporting ATPase turnover rates. They have been calculated from radioactivity measurements and reported for every Na/K-ATPase to be around 3,600/min (Na+) (24), 3,000/min (K+, Rb+) (25), and 2,000/min (ATP) (26). Since every ATP molecule results in transfer of two K+ ions into the cells and three Na+ outside, the average rate for Na+ transport could be estimated to be around 4,700/min. Since the dominant counterion in the media is Cl−, sodium efflux can be estimated to be responsible for one-half of the osmolyte flux, or 3.3 × 108/min. Using confocal imaging, we have estimated experimental embryos to contain 60 cells on average, 42 of which belong to the TE layer (SI Appendix, Fig. S1 C–E), which would result in ion flux per cell of 7.8 × 106/min, estimating 1.7 × 103 active ion transporters per cell. In addition to estimating the average number of active Na/K-ATPases operating in TE cells, we also looked at how embryos responded to rapid changes in osmolarity, which revealed water transfer kinetics reaching rates of 400 µm3/s, as 252 mOsm of sucrose was added to embryos developing in channels (SI Appendix, Fig. S4 A and B). Water flux during transient response is over 20 times the normal water transport rates seen during normal development (Table 1), suggesting that cavity expansion is limited and regulated by the rate of ion rather than water transport.
Summary of embryonic microphysiology
Mechanical Model of Zona Pellucida.
Although liquid inflow and sealed epithelial layer are the prerequisites for cavity formation, pressure ultimately arises due to the constricting nature of the external polysaccharide shell zona pellucida. Since shell thickness can be measured by bright-field microscopy to be less than 1/10 of the embryo diameter, the blastocyst can be modeled as a thin-walled pressure vessel (Eqs. 4–9). The model also leads to the conclusion that there exists a critical pressure and radial expansion beyond which further expansion and shell thickness decrease occurs without increase in pressure. This is broadly analogous to a balloon bursting after it is inflated beyond the critical pressure. As such, the model demonstrates that zona pellucida rupture and embryo hatching is driven by volumetric expansion, not pressure generation, as embryos do not spontaneously hatch after reaching the peak pressure but require zona thickness to decrease further until breaking. All embryos analyzed here were similar to the initial stages of balloon inflation and did not increase their radius significantly beyond the point of maximum pressure.
The elastic balloon model can be used to calculate the shell material stiffness: if the shell increases in radius from R0 to R when pressurized to P, the shell stiffness is
Characterization of blastocyst hatching. (A) Pressure required to expand the embryo according to the thin-walled pressure vessel mathematical model. After an embryo reaches the maximum pressure at ∼50% radial expansion, its further expansion does not require additional pressure. The model was calculated with stiffness of E = 31.0 kPa, mean shell thickness of t = 4.3 µm, and initial radius of R0 = 40 µm, which are corresponding to median measured embryo values. (B) Measured zona pellucida stiffness during early development (n = 38). The stiffness was calculated from measured internal pressure, estimated equivalent sphere radius based on volume, and the mean starting radius and mean shell thickness. Eventual decrease in zona pellucida thickness leads to the hatching of an embryo, which may happen through a pinhole or rupture mechanism as depicted in C and D, respectively (Movies S3 and S4). (Scale bars: 100 µm.) (E) Measured probability of embryo hatching, as a function of embryo pressure and zona pellucida thickness. The probability map was constructed from shell thickness and pressure measurements, just before the hatching events. Hatching probability approached certainty as zona pellucida thickness decreased below 2.5 µm, as shown in a probability map. We found there to be a smaller dependence on blastocyst cavity pressure as soon as it increased over 500 Pa. However, the pressure and volumetric expansion did play a role in the hatching process as shown in the Cox survival model (F). Dashed lines indicate 95% confidence intervals; P value (from Cox regression model) for 150 mM sucrose treatment was 0.0209 and for 0.5 mg/mL collagenase treatment was 0.0014.
The model also makes it possible to estimate that, given the mean rate of volume increase of 22 µm3/s, it would take 16.4 h for an embryo of 40-µm radius to increase its radius to 72 µm at which point shell thickness would decrease to one-third of the original size, at which point the embryos would be expected to hatch. The timing coincides with developmental timescale. The decrease of zona pellucida thickness can be explained and modeled by embryo microphysiology (Table 1) and the elastic model. However, in practice we found the hatching process to be more complex, defined by two distinct modes of zona failure and being dependent also on matrix degrading enzymes, which we therefore investigated.
Embryo Hatching Probability.
After the pressure buildup and volume increase, the ultimate result is embryo hatching. We identified two distinct types, which we termed “pinhole” (Fig. 3C and Movie S3) and “rupture” (Fig. 3D and Movie S4). The first type of hatching could occur at relatively high zona pellucida thickness of over 3 µm and represented 45% of the total hatching events (n = 19 of 42). The remaining embryos (55%) hatched by rupturing the shell after its thickness decreased below 2.5 µm (Fig. 3E and SI Appendix, Fig. S4C), often after partially hatching through the pinhole mechanism. This process complicated data interpretation, as hatching events classified as pinhole would be converted into rupture depending on the developmental timeline. The overall hatching probability can be summarized as a function of embryo pressure and zona thickness measurements, which were taken and interpolated during hatching events (Fig. 3E). Despite embryo hatching probability approaching unity as zona thickness decreased, the embryos still required at least 0.5 kPa to hatch.
It has previously been shown that matrix degrading enzymes are important for shell breakage (20⇓–22). To compare the relative contribution of mechanical forces and enzymatic digestion to hatching, we used media containing 0.5 mg/mL collagenase I to erode the zona pellucida and determined its effect on hatching probability. To test the role of pressure, we used 150 mM sucrose to lower the internal pressure. The experiment was fitted by the Cox survival model (Fig. 3F), and Kaplan–Meier statistic was then used to estimate the significance between treatment and control groups. Among the control group and sucrose-treated embryos, the likelihood coefficient was −1.075, which meant that for any given point in time, sucrose-treated embryos were nearly three times less likely to hatch than controls (value of P of 0.021). Conversely, collagenase-treated embryos were 3.8 times more likely to hatch than controls (P value of 0.014).
Quantifying Microphysiology of Cryopreserved Embryos.
The cryopreservation and thawing of preimplantation embryos is common during IVF, but little is known about what effect, if any, cryopreservation has upon blastocyst microphysiology and the pressure the blastocyst is able to generate. We therefore used our quantitative approach to look at how the cryopreservation/thawing process influences the embryo. We used three groups of mouse embryos, all from a C57BL/6J background: freshly flushed embryos at 2.5 dpc as controls and two groups of thawed embryos that had been cryopreserved at the 2.5 or 3.5 dpc stage. Cryopreserved embryos had been held in liquid nitrogen for at least 4 y. All embryos were cultured and measured at the stage equivalent to 3.5 dpc. We found no statistically significant differences between embryos thawed at stages 2.5 and 3.5 dpc. However, both thawed groups generated significantly lower pressure compared with freshly flushed embryos (Fig. 4A) (t test, P values < 0.05). This indicates that cryopreservation alters embryos in such a way that either their zona is less stiff, that is, “softer” (Fig. 4 B and C), or that they are unable to develop the normal levels of volume growth required to generate pressure. To determine the consequences of reduced blastocyst pressure in cryopreserved embryos, we assessed hatching efficiency and found that hatching was delayed by ∼10 h in comparison with noncryopreserved control embryos (SI Appendix, Fig. S4D). Taken together, our quantitative measurements of pressure reveal that cryopreservation alters mouse embryos in a previously unanticipated manner that has implications for hatching efficiency.
Measuring the influence of cryopreservation on blastocyst pressure. Embryos cryopreserved at 2.5 dpc (n = 15) and 3.5 dpc (n = 12) were thawed and compared with embryos freshly flushed at 2.5 dpc (n = 33). All measurements were made at the equivalent of E3.5. (A) Cryopreserved embryos generate significantly lower pressure compared with freshly collected embryos. (B) Blastocyst volume growth is similar in the three groups. (C) Zona pellucida stiffness shows a trend toward being reduced in cryopreserved embryos. Calculations of volume growth and stiffness require time derivatives, which introduce more noise and dispersion, making statistical significance harder to establish despite large differences among the means. Outliers have been excluded for clarity in the plots in B and C. *P < 0.5 and **P < 0.01, pairwise t test.
Conclusion
Mammalian embryo hatching is a biological process that is difficult to study quantitatively and model using existing approaches. We have developed a noninvasive technique to quantify embryo microphysiological parameters, which is useful not only for studying fundamental biological processes but also for evaluating embryo quality in assisted reproductive approaches. We applied this technique to developing mouse embryos, detailing microphysiological parameters as well as zona pellucida stiffness. Measuring these aspects of blastocyst development revealed the magnitude of hydraulic pressures inside blastocysts and allowed us to develop a mathematical model based on inflating a thin-shell pressure vessel. Using this framework, we showed that hatching is limited by the rate of ion transport, which ultimately draws in water, increases embryo volume, and thereby reduces zona thickness. By experimentally altering blastocyst pressure, we were able to investigate how hatching dynamics are defined by both physiology and mechanical properties of the zona.
We used our approach to study the effect of cryopreservation on blastocysts and found that thawed embryos generated a significantly lower pressure compared with fresh embryos, accompanied by a considerable delay in hatching. Hatching is a prerequisite to implantation and continued development of the conceptus. Assisted reproductive approaches increasingly utilize cryopreservation, highlighting the importance of evaluating in quantitative detail its effect not only on hatching, but other aspects of embryo microphysiology. We focused on the mouse model due to its long-standing importance in in vitro studies. Other mammalian species may have different embryo microphysiology and zona pellucida properties, and it will be interesting to uncover what differences, if any, exist. Finally, our method is broadly applicable for characterizing pressurized cavities in other contexts, such as stem cell-derived organoids or cysts.
Materials and Methods
Preparing Cylindrical Hydrogels.
Hydrogels with cylindrical cavities were prepared by casting polyacrylamide gels around small-diameter wire (40 µm). Acrylamide/bisacrylamide (37.5:1) mixture was used (5.0% by mass), which corresponded to 4.2-kPa stiffness. Ammonium persulfate (0.1%) and tetramethylethylenediamine (1%) were added to polymerize polyacrylamide. The wire was removed, forming cylindrical cavities, and the hydrogels equilibrated in embryo culture medium by changing the solution three times in 1 h. Embryos were then inserted into the channels using a glass capillary having a diameter slightly larger than the embryo, which was used to stretch open the hydrogel channel before injecting the embryos and letting the channels relax and deform the embryos.
Imaging Setup and Image Analysis.
Imaging was performed on a Zeiss Axiovert spinning-disk confocal microscope equipped with a Hamamatsu emCCD camera and a 37 °C and 5% CO2 incubator. Embryos were cultured for periods of 2 h between 3.0 and 4.0 dpc stages, recording images every 6 min. Embryos were analyzed using a combination of manual and automatic image analysis approaches. Manual analysis was used to ensure no algorithmic bias, while the high-throughput analysis relied on computational methods. The main measured parameters were channel diameter and embryo width and length, which made it possible to calculate microphysiological embryo properties. Zona thickness was measured only manually, at either of the unconfined embryo ends. Blastocyst cavity volume was determined from the total embryo volume by subtracting a constant volume of the inner cell mass cells assumed to be 2.1 × 10−5 um3. Average volumes and pressures were calculated by taking mean values of an embryo’s time lapse, summarizing them, and taking means of the means (data seen in histograms in Fig. 2).
Thin-Walled Vessel Pressure Derivation.
As is familiar to anyone who has inflated a party balloon, the inflation is difficult at first, but once the balloon is modestly inflated, it becomes much easier. This occurs because the pressure in the balloon first increases with inflation, and then, past a critical volume, the pressure in the balloon falls. This “balloon instability” ultimately has a geometric origin: As the balloon inflates, its walls become thinner and flatter and thus able to contain less pressure. We can understand this phenomenon (28) by considering a spherical rubber balloon, with natural radius R0 and thickness t. If the balloon is inflated such that its radius dilates to R = λR0, then circumferences of the balloon dilate by λ, while to conserve volume of the rubber, the thickness of the balloon must fall to tλ2. Therefore, the deformation gradient in the rubber will be as follows:
and assuming the rubber is well described by a neo-Hookean energy, the total elastic energy will be as follows:
Since in this inflation the volume of the balloon changes from its undeformed volume V0 to V = λ3V0, we can express this elastic energy as follows:
and we can calculate the pressure in the balloon as follows:
This pressure can be recast in terms of R and R0 as follows:
This form for the pressure does indeed rise from zero for modest inflations (R a little above R0), but then reaches a maximum pressure when
which occurs when R = 71/6R0 ∼ 1.38 R0, and the pressure falls with subsequent additional inflation. The form of pressure–inflation curve is shown in Fig. 3A, the essential point being that if the embryo can exert enough pressure to inflate the zona (balloon) radius by 38%, then it does not require additional pressure to expand the shell further until it breaks. The exact value of 38% is particular to the neo-Hookean rubber model, but calculations for other material models, and experiments with real balloons, indicate that the critical dilation reliably lies between 71/6 and 1.5 (29). Also, for thin-walled cylindrical pressure vessels, the longitudinal stress can be expressed as a combination of pressure, wall thickness, and the radius of the spherical part:
Acknowledgments
This work was supported by a studentship (to K.L.) from the Biotechnology and Biological Sciences Research Council (BB/J014427/1), Lithuanian Science Council Postdoctoral Award (code 09.3.3-LMT-K-712-02-0067), the Welcome Trust (203141/Z/16/Z), and a Wellcome Senior Investigator Award (103788/Z/14/Z) (to S.S.).
Footnotes
↵1Present address: Institute of Biotechnology, Vilnius University, Vilnius 10257, Lithuania.
- ↵2To whom correspondence should be addressed. Email: shankar.srinivas{at}dpag.ox.ac.uk.
Author contributions: K.L., C.R., and S.S. designed research; K.L., C.R., and C.P. performed research; C.P., B.D., and J.S.B. contributed new reagents/analytic tools; K.L., C.R., J.S.B., and S.S. analyzed data; and K.L., C.R., J.S.B., and S.S. wrote the paper.
The authors declare no conflict of interest.
This article is a PNAS Direct Submission.
This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1719930115/-/DCSupplemental.
- Copyright © 2018 the Author(s). Published by PNAS.
This open access article is distributed under Creative Commons Attribution-NonCommercial-NoDerivatives License 4.0 (CC BY-NC-ND).
References
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵
- ↵.
- du Roure O, et al.
- ↵.
- Li M, et al.
- ↵
- ↵.
- Khalilian M,
- Navidbakhsh M,
- Valojerdi MR,
- Chizari M,
- Yazdi PE
- ↵
- ↵
- ↵
- ↵.
- Plusa B, et al.
- ↵.
- Watson AJ,
- Natale DR,
- Barcroft LC
- ↵
- ↵
- ↵
- ↵
- ↵.
- Thomas M,
- Jain S,
- Kumar GP,
- Laloraya M
- ↵
- ↵
- ↵.
- Biggers JD,
- Bell JE,
- Benos DJ
- ↵
- ↵.
- Pollack LR,
- Tate EH,
- Cook JS
- ↵.
- El Mernissi G,
- Doucet A
- ↵
- ↵.
- Müller I,
- Strehlow P
- ↵
Citation Manager Formats
Sign up for Article Alerts
Jump to section
You May Also be Interested in
More Articles of This Classification
Biological Sciences
Developmental Biology
Physical Sciences
Related Content
Cited by...
- No citing articles found.