Review Special Issues

CRISPR engineering cardiometabolic disease models using human iPSC

  • Since their introduction, human induced pluripotent stem cells (iPSCs) have enhanced the ways human disease processes are evaluated. Only eleven years ago, Yamanaka and colleagues showed somatic cells could be reprogrammed to induced pluripotent cells (iPSC) employing four transcription factors. Since that initial demonstration, much progress has been made in establishing human cardiac cell models that recapitulate diseases using easily obtained patient somatic cells such as fibroblasts or peripheral blood cells. Investigators now have access, through established small molecule methods to differentiate iPSCs to cardiomyocytes, hepatocytes, neurons, astrocytes and many other cell types of interest. This capability provides valuable cell models of human disease derived from patient cells or CRISPR gene-edited cells that introduce or repair patient disease mutations. In many cases, these manipulations create unique cell models and address the need for scalability and reproducibility of cell samples. Consequently, human iPSC derived cells can be used to query cellular and pathological mechanisms potentially providing bio-information that animal models may be unable to report. In this article, we review iPSC derived cardiomyocyte models of cardiometabolic disease and demonstrate the power of these technologies to probe gene function for biological and clinical significance and to capture pathological signatures from genetic variation in human cardiac disease. Examples of the iPSC derived cardiometabolic disease models including channelopathies, metabolic disease and cardiomyopathy are presented.

    Citation: Cindy E. McKinney, Katherine M. Baumgarner. CRISPR engineering cardiometabolic disease models using human iPSC[J]. AIMS Cell and Tissue Engineering, 2018, 2(3): 185-202. doi: 10.3934/celltissue.2018.3.185

    Related Papers:

    [1] Martin Gugat, Mario Sigalotti . Stars of vibrating strings: Switching boundary feedback stabilization. Networks and Heterogeneous Media, 2010, 5(2): 299-314. doi: 10.3934/nhm.2010.5.299
    [2] Zhong-Jie Han, Gen-Qi Xu . Spectrum and dynamical behavior of a kind of planar network of non-uniform strings with non-collocated feedbacks. Networks and Heterogeneous Media, 2010, 5(2): 315-334. doi: 10.3934/nhm.2010.5.315
    [3] Yacine Chitour, Guilherme Mazanti, Mario Sigalotti . Stability of non-autonomous difference equations with applications to transport and wave propagation on networks. Networks and Heterogeneous Media, 2016, 11(4): 563-601. doi: 10.3934/nhm.2016010
    [4] Yaru Xie, Genqi Xu . The exponential decay rate of generic tree of 1-d wave equations with boundary feedback controls. Networks and Heterogeneous Media, 2016, 11(3): 527-543. doi: 10.3934/nhm.2016008
    [5] Avner Friedman . PDE problems arising in mathematical biology. Networks and Heterogeneous Media, 2012, 7(4): 691-703. doi: 10.3934/nhm.2012.7.691
    [6] Kaïs Ammari, Mohamed Jellouli, Michel Mehrenberger . Feedback stabilization of a coupled string-beam system. Networks and Heterogeneous Media, 2009, 4(1): 19-34. doi: 10.3934/nhm.2009.4.19
    [7] Serge Nicaise, Julie Valein . Stabilization of the wave equation on 1-d networks with a delay term in the nodal feedbacks. Networks and Heterogeneous Media, 2007, 2(3): 425-479. doi: 10.3934/nhm.2007.2.425
    [8] Dongyi Liu, Genqi Xu . Input-output L2-well-posedness, regularity and Lyapunov stability of string equations on networks. Networks and Heterogeneous Media, 2022, 17(4): 519-545. doi: 10.3934/nhm.2022007
    [9] Tong Li . Qualitative analysis of some PDE models of traffic flow. Networks and Heterogeneous Media, 2013, 8(3): 773-781. doi: 10.3934/nhm.2013.8.773
    [10] Luca Schenato, Sandro Zampieri . On rendezvous control with randomly switching communication graphs. Networks and Heterogeneous Media, 2007, 2(4): 627-646. doi: 10.3934/nhm.2007.2.627
  • Since their introduction, human induced pluripotent stem cells (iPSCs) have enhanced the ways human disease processes are evaluated. Only eleven years ago, Yamanaka and colleagues showed somatic cells could be reprogrammed to induced pluripotent cells (iPSC) employing four transcription factors. Since that initial demonstration, much progress has been made in establishing human cardiac cell models that recapitulate diseases using easily obtained patient somatic cells such as fibroblasts or peripheral blood cells. Investigators now have access, through established small molecule methods to differentiate iPSCs to cardiomyocytes, hepatocytes, neurons, astrocytes and many other cell types of interest. This capability provides valuable cell models of human disease derived from patient cells or CRISPR gene-edited cells that introduce or repair patient disease mutations. In many cases, these manipulations create unique cell models and address the need for scalability and reproducibility of cell samples. Consequently, human iPSC derived cells can be used to query cellular and pathological mechanisms potentially providing bio-information that animal models may be unable to report. In this article, we review iPSC derived cardiomyocyte models of cardiometabolic disease and demonstrate the power of these technologies to probe gene function for biological and clinical significance and to capture pathological signatures from genetic variation in human cardiac disease. Examples of the iPSC derived cardiometabolic disease models including channelopathies, metabolic disease and cardiomyopathy are presented.


    Within the framework of a nonlocal SIR model, we introduce different vaccination policies and then seek the most effective one. Vaccines are dosed at prescribed ages or at prescribed times to prescribed portions of the susceptible population. The careful choice of the ages, times and individuals to be treated may significantly influence the efficiency of a vaccination campaign. It is thus natural to tackle the resulting optimal control problem, seeking what is in some sense the best vaccination policy.

    We denote by S=S(t,a) the density of individuals at time t of age a susceptible to the disease. As usual, I=I(t,a) is the density of infected individuals at time t and of age a. The density of individuals that can not be infected by the disease is R=R(t,a), comprising individuals that recovered as well as those that are, whatever the reason, immune. As a reference to the basic properties of SIR models, originated in [1], see [2,Chapter 6], [3,Chapter 19] or [4,§ 1.5.1].

    A possible description of the evolution of S is provided by the renewal equation [4,Chapter 3] with nonlocal source term

    tS+aS=(dS(t,a)++0λ(a,a)I(t,a)da)S (1.1)

    where dS=dS(t,a)0 is the susceptibles' death rate and λ(a,a)0 quantifies the susceptible individuals of age a that are infected by individuals of age a. The nonlocal quantity on the right hand side of (1.1), namely +0λ(a,a)I(t,a)daS(t,a), represents the total number of susceptible individuals of age a that become infected at time t per unit time.

    The evolution of the infected portion I of the population is governed by

    tI+aI=(dI(t,a)+rI(t,a))I++0λ(a,a)I(t,a)daS, (1.2)

    where dI=dI(t,a)0 is the death rate of the I portion of the population and rI=rI(t,a)0 represents the fraction of infected individuals of age a that recovers at time t per unit time. Finally, the evolution of the portion R of the population that recovers or is unaffected from the disease is

    tR+aR=rI(t,a)IdR(t,a)R, (1.3)

    where dR=dR(t,a)0 is the death rate of the R population. We assume here that members of the R population are not dosed with vaccine and that individuals may enter the R population at birth. In the present setting, dealing with choices different from these latter ones only requires very minor changes and keeps falling in the present general analytic framework.

    In a first policy, vaccinations are dosed at any time t at a time dependent portion ηj(t) of the S population of the prescribed age ˉaj, with 0<ˉa1<<ˉaN. In this way, the fraction ηj(t) of susceptible population S(t,ˉaj) becomes immune. In other words, we use as control the fractions ηj:I[0,1] (for j{1,,N} and I being the, possibly unbounded, time interval under consideration). We thus have to supplement the evolution (1.1)–(1.2)–(1.3) with the vaccination conditions

    S(t,ˉaj+)=(1ηj(t))S(t,ˉaj)[tS(t,ˉaj) decreases due to vaccination]I(t,ˉaj+)=I(t,ˉaj)[the infected population is unaltered]R(t,ˉaj+)=R(t,ˉaj)+ηj(t)S(t,ˉaj)[vaccinated individuals are immunized] (1.4)

    for a.e. t>0 and for j{1,,N}. A reasonable constraint on the vaccination policy η in (1.4) is that its total cost N(η) on the considered time interval I does not exceed a prescribed maximal cost M. Measuring the total cost by the number N of dosed individuals, we have

    N(η)M where N(η)=Ni=1Iηi(t)S(t,ˉai)dt. (1.5)

    Alternatively, we also consider vaccination campaigns immunizing an age dependent portion of the whole population at given times ˉt1,,ˉtN. This amounts to substitute (1.4) with

    S(ˉtk+,a)=(1νk(a)ηj)S(ˉtk,a)[aS(ˉtk,a) decreases due to vaccination]I(ˉtk+,a)=I(ˉtk,a)[the infected population is unaltered]R(ˉtk+,a)=R(ˉtk,a)+νk(a)S(ˉtk,a)[vaccinated individuals are immunized] (1.6)

    where νk(a) is the percentage of susceptible individuals of age a that are dosed with the vaccine at time ˉtk, for k=1,,N. Now, a reasonable constraint due to the campaign cost is given by a bound on the number N of dosed individuals

    N(ν)M whereN(ν)=Ni=1R+νi(a)S(ˉti,a)da. (1.7)

    For simplicity, in both cases above, we assume that ηj, respectively νi, represents the percentages of individuals that are successfully vaccinated.

    In both cases, our aim is to seek the vaccination campaign that minimizes the total number of infected individuals of all ages over the time interval I, that is

    J=IR+φ(a)I(t,a)dadt (1.8)

    the function φ being a suitable positive age dependent weight. Here, J is a function of η in case (1.4) and of ν in case (1.6). In general, in the two cases (1.4) and (1.6), reasonable costs are thus

    J(η)+N(η)=IR+φ(a)I(t,a)dadt+Nj=1Iηj(t)S(t,ˉaj)dt (1.9)
    J(ν)+N(ν)=IR+φ(a)I(t,a)dadt+Ni=1R+νi(a)S(ˉti,a)da. (1.10)

    Clearly, in general, suitable weights can be used to modify the relative relevance of the two costs N and J in the sums above.

    The current literature offers a variety of alternative approaches to similar modeling situations. In the context of models based on ordinary differential equations, that is, without an age structure in the population, this problem has been considered, for instance, already in [5,6,7]. In the recent [8], the vaccination control enters an equation of the type (1.1) through a term uS in the source on the right hand side of (1.1), implicitly suggesting that vaccination takes place uniformly at all ages. For a recent related investigation focused on cholera, see [9].

    The next section summarizes the key analytic properties of the present SIR model, providing the necessary basis to tackle optimal control problems. Section 3 shows the qualitative properties of various vaccination strategies through numerical integrations. Section 4 presents numerical optimizations through the descent method. The paper ends with Section 5.

    Denote by I either the time interval [0,T], for a positive T, or [0,+[. Fix throughout a positive integer N representing the number of vaccination sessions.

    Both models introduced above lead to the initial – boundary value problem

    {tS+aS=(dS(t,a)++0λ(a,a)I(t,a)da)StI+aI=(dI(t,a)+rI(t,a))I++0λ(a,a)I(t,a)daStR+aR=rI(t,a)IdR(t,a)R (2.1)

    with initial datum and boundary inflow

    {S(0,a)=So(a),I(0,a)=Io(a),R(0,a)=Ro(a),aR+;S(t,0)=Sb(t),I(t,0)=Ib(t),R(t,0)=Rb(t),tI. (2.2)

    The usual case of the birth terms in the boundary data being assigned through suitable integrals of the populations in reproductive age can be recovered by mainly technical adjustments.

    We assume below that the following assumptions are satisfied:

    (λ) λC0(R+×R+;R) is bounded, with total variation in the first argument uniformly bounded with respect to the second (i.e., supaR+TV(λ(,a);R+)<+) and Lipschitz continuous in the first argument uniformly with respect to the second (i.e., there exists a C>0 such that for all a,a1,a2R+, |λ(a1,a)λ(a2,a)|C|a1a2|).

    (d-r) The maps dS,dI,dR,rI:I×R+R are bounded Caratheodory functions, with L1 norm and total variation in a uniformly bounded in t (i.e., for φ=dS,dI,dR,rI, suptIφ(t)L1(R+;R)<+ and suptITV(φ(t,);R+)<+) and Lipschitz continuous in a uniformly in t (i.e., there exists a C>0 such that for all tI and a1,a2R+, |φ(t,a2)φ(t,a1)|C|a2a1|).

    (IB) The initial datum and the boundary inflow in (2.2) satisfy So,Io,Ro(L1BV)(R+;R+) and Sb,Ib,Rb(L1BV)(I;R+).

    For the definition of Caratheodory function, of total variation and of the various functional spaces, we refer for instance to [10]. For all other analytic details specific to the present construction, see [11].

    First, we provide the basic well posedness and stability result for the model (2.1)–(2.2) in the case of the vaccination policy (1.4), namely

    {S(t,ˉaj+)=(1ηj(t))S(t,ˉaj)I(t,ˉaj+)=I(t,ˉaj)R(t,ˉaj+)=R(t,ˉaj)+ηj(t)S(t,ˉaj)j=1,,N. (2.3)

    Theorem 2.1 ([11]). Under hypotheses (λ) and (d-r), for any initial and boundary data satisfying (IB), for any choice of the positive vaccination ages ˉa1,,ˉaNR+{0} and of the control function ηBV(I;[0,1]N), problem (2.1)–(2.2)–(2.3) admits a unique solution defined on the whole interval I, depending Lipschitz continuously on the initial datum, through the L1 norm, and on η, through the L norm.

    An entirely similar result holds for the vaccination policy (1.6), which we rewrite here as

    {S(ˉti+,a)=(1νi(a)ηj)S(ˉti,a)I(ˉti+,a)=I(ˉti,a)R(ˉti+,a)=R(ˉti,a)+νi(a)S(ˉti,a).i=1,,N. (2.4)

    Theorem 2.2 ([11]). Under hypotheses (λ) and (d-r), for any initial and boundary data satisfying (IB), for any choice of the positive vaccination times ˉt1,,ˉtNI{0} and of the control function νBV(R+;[0,1]N), problem (2.1)–(2.2)–(2.4) admits a unique solution defined on the whole interval I, depending Lipschitz continuously on the initial datum, through the L1 norm, and on ν, through the L norm.

    Below, we restrict the search of optimal controls ηj in (1.4) and νi in (1.6) to functions of the form

    ηj(t)=m=1ηjχ[t1,t[(t)t0=0,t>t1,tI;νi(a)=m=1νiχ[a1,a[(a)a0=0,a>a1,aR+.

    The stability results above ensure that the costs (1.5), (1.7) and (1.8) are Lipschitz continuous functions of the various ηj and νi. Hence, a straightforward Weierstraß argument proves the existence of optimal controls minimizing the costs (1.9) or (1.10).

    We refer to [11] for the detailed proofs as well as for further stability estimates on the dependence of the solutions on the various parameters.

    We now compare different instances of both strategies (2.3) and (2.4) by means of numerical integrations.

    In all the integrations below, for the convective part we exploit the upwind method [12,§ 4.8] with mesh size Δa=2.5103 along the age axis and Δt=1.25103 along the time axis. The integrals in the right hand side in (2.1) are computed via a rectangle rule and a fractional step method [12,§ 17.1] allows to combine the convective evolution with the source term. The time axis a=0 belongs both to the physical boundary and to the numerical one, its treatment being straightforward since all characteristics move inward the domain I×R+, so that the value of the boundary data can be assigned to the solution. On the other side, along the numerical boundary a=10, the usual free flow conditions is consistent with all characteristic speeds being positive.

    Our aim in the integrations below is to clearly show the qualitative properties of the model (2.1). Indeed, we do not claim that the chosen numerical values are realistic. In particular, the most adequate units for t and a might well be different.

    With reference to the Cauchy problem (2.1)–(2.2) and to the costs (1.9)–(1.10), we use throughout the following parameters and functions:

    ds(t,a)=14a1+a,dI(t,a)=12a1+a,dR(t,a)=14a1+a,rI(t,a)=110,λ(a,a)=2e|aa|/2,φ(a)=1. (3.1)

    The mortality rates are increasing with age and higher for the infected population. The choice of λ assumes that the infection propagates preferably among individuals of similar age.

    We consider two different initial data, first the case without infection, i.e.,

    So(a)=3(1a10),Io(a)=0,Ro(a)=(1a10), (3.2)

    and then a case with a 20% of infected individuals, with respect to the S population, at time t=0, that is

    So(a)=2.5(1a10),Io(a)=0.5(1a10),Ro(a)=(1a10). (3.3)

    Both these initial data are linearly decreasing with age. Note that the total initial population is the same in the two cases (3.2) and (3.3).

    Natality is chosen so that no jump appears between the initial values at a=0+ in (3.2) and the boundary data at t=0+:

    Sb(t)=3,Ib(t)=0,Rb(t)=1, (3.4)

    where we assume that no individuals are infected already at birth.

    We first consider two extreme situations to be later used for comparisons, namely corresponding to no illness (3.2), and to a spreading illness with no vaccination strategy whatsoever.

    In other words, we integrate (2.1)–(2.2)–(3.1) first setting Io0, so that for all times t, I(t)0 and no illness ever appears, and then setting Io as in (3.3) but dose no vaccination, so that the epidemic spreads uncontrolled. The results of these two integrations are displayed in Figure 1. The two solutions are, clearly, entirely different. When I0, a sort of dynamic equilibrium is reached, the resulting contour plot being approximately invariant with respect to translation along the vertical t axis. The presence of the I population in the initial datum leads to the spread of the disease and eventually to a dramatic population decrease.

    Figure 1.  Integration of (2.1), with parameters (3.1), boundary data (3.4), initial data on the left corresponding to (3.2) and, middle and right, to (3.3). The left and middle plots display the sum S+R as a function of t and a, while the rightmost diagram displays the ratio I/(S+I+R) as a function of t and a.

    We consider now a simple vaccination strategy consisting in a single campaign, dosing vaccination at a single age ˉa for all times.

    As it is to be expected, the earlier vaccinations are dosed, the better it is in terms of the cost (1.8). In Figure 2 the ratio of the number of infected individuals vs. the total population corresponding to the vaccinations ages ˉa=0.1,0.3,0.5,1.5 are portrayed. The leftmost contour plot shows the effects of a very early vaccination, clearly far more effective than later doses.

    Figure 2.  Ratio of the infected population over the whole population in the integration of (2.1), with parameters (3.1), boundary data (3.4) and initial data as in (3.3). The different integrations result from different choices of the age ˉa at which all S individuals are vaccinated: from left to right: ˉa=0.1,0.3,0.5,1.5.

    The striking effect of this early age vaccination is confirmed in Figure 3, where we display the sum S+R as a function of t and a in the same cases of the vaccinations ages ˉa=0.1,0.3,0.5,1.5. Clearly, due to the form of the vaccination condition (2.3) and to the present smooth data (3.3)–(3.4), the sum S+R is continuous across ˉa. In the leftmost diagram, the healthy population apparently recovers from the spread of the disease. When vaccinations are dosed at a later age, i.e., in the three diagrams on the right, the relatively large portion of individuals initially vaccinated survive the disease, but juveniles get infected resulting in the solitary waves, consisting essentially of the R population, shown in Figure 3.

    Figure 3.  Healthy, i.e., S+R, population resulting from the integration of (2.1), with parameters (3.1), boundary data (3.4) and initial data as in (3.3). The different integrations result from different choices of the age ˉa, at which all S individuals are vaccinated, where, from left to right: ˉa=0.1,0.3,0.5,1.5.

    The counterpart of these solitary waves is the propagating "hole" in Figure 2.

    In other words, the vaccination campaign at the age ˉa=0.5 looks as (un)successful as the campaign at ˉa=1.5. On the contrary, dosing at ˉa=0.1 is by far more effective.

    Indeed, this strategy keeps the density of the susceptible population in (2.1) very low due to (2.3). Hence, also the nonlocal source term, which is proportional to S, is relatively small and infection hardly spreads. When vaccination is dosed at higher ages, the nonlocal source term is able to bring infection at low ages, thus making the vaccination policy far less effective.

    Table 1 confirms that, in particular for what concerns the total number of infected individuals, by far, the best strategy is the one dosing vaccines at the earliest considered age ˉa=0.1. This choice also results in the highest number of vaccinated individuals, so that also N corresponds to the most expensive policy.

    Table 1.  Value of the functionals in (1.5), (1.8), and (1.9) corresponding to the numerical integration of (2.1) with parameters (3.1), initial datum (3.3) and boundary datum (3.4). The rightmost column refers to the case where no vaccination is dosed.
    ˉa in (2.3) 0.1 0.2 0.3 0.5 1.0 1.5 none
    J in (1.8) 33.48 41.03 44.73 46.85 47.39 47.51 48.30
    N in (1.5) 10.53 5.05 2.30 0.76 0.41 0.36 0.00
    J+N in (1.9) 44.01 46.08 47.03 47.61 47.80 47.87 48.30

     | Show Table
    DownLoad: CSV

    We now simulate another single vaccination campaign, consisting in dosing all individuals of all ages at a given time ˉt, see Figure 4.

    Figure 4.  Ratio of the infected population over the whole population in the integration of (2.1), with parameters (3.1), boundary data (3.4) and initial data as in (3.3). The different integrations result from different choices of the time ˉt, at which all S individuals are vaccinated where from left to right, ˉt=0.1,0.3,0.5,1.5.

    Similarly to the integration detailed in § 3.2, campaigns at earlier times results in being more effective. The choices ˉt=0.5 and ˉt=1.5 are quite similar to each other and, in the lower right part of the (a,t)–plane, to the situation with no vaccination campaign whatsoever, as it stems from a comparison between the two rightmost diagram in Figure 4 and the rightmost one in Figure 1, see also Table 2.

    Table 2.  Value of the functionals in (1.7), (1.8) and (1.10) corresponding to the numerical integration of (2.1) with parameters (3.1), initial datum (3.3) and boundary datum (3.4). The rightmost column refers to the case where no vaccination is dosed.
    ˉt in (2.4) 0.1 0.2 0.3 0.5 1.0 1.5 none
    J in (1.8) 31.71 38.80 43.65 46.79 47.40 47.42 48.30
    N in (1.7) 9.37 5.46 2.74 0.86 0.40 0.39 0.00
    J+N in (1.10) 41.08 44.26 46.39 47.66 47.80 47.80 48.30

     | Show Table
    DownLoad: CSV

    Note that the present strategy is far less effective: vaccinations are dosed only at the given time ˉt, after which the disease spreads completely out of control, see Figure 5.

    Figure 5.  Healthy, i.e., S+R, population resulting from the integration of (2.1), with parameters (3.1), boundary data (3.4) and initial data as in (3.3). The different integrations result from different choices of the time ˉt, at which all S individuals are vaccinated where from left to right, ˉt=0.1,0.3,0.5,1.5.

    The values of the costs (1.7), (1.8) and (1.10) presented in Table 2 show that vaccinating the whole S population at ˉt=0.5 or later, makes relatively small difference, the most convenient choice being at the earliest time ˉt=0.1. Similarly to what happens in § 3.2, vaccinating all susceptible individuals at ˉt=0.1 makes the susceptible portion of the population very small wherever a>tˉt, thus essentially canceling the nonlocal source term which represents the spreading of the disease.

    We note however an evident key difference between the best vaccination age choice (i.e., ˉa=0.1) in § 3.2 and the present best vaccination time choice (i.e., ˉt=0.1). The latter one yields lower values of the chosen costs. However, as time grows, the former strategy is prone to by far preferable evolutions. Indeed, the leftmost diagram in Figure 2 (corresponding to ˉa=0.1) clearly shows a decreasing I population, while it is increasing in the latter strategy, see Figure 4, left (corresponding to ˉt=0.1).

    We now consider in more detail the search for an optimal vaccination strategy in a slightly different setting. First, we assume that infection among individuals of different ages is somewhat more difficult, slightly modifying (3.1) to

    ds(t,a)=14a1+a,dI(t,a)=12a1+a,dR(t,a)=14a1+a,rI(t,a)=110,λ(a,a)=e2|aa|,φ(t)=1, (4.1)

    and choose the initial datum

    So(a)={3.0(1a10)a[0,3][5,10]1.5(1a10)a]3,5[Io(a)={0a[0,3][5,10]1.5(1a10)a]3,5[Ro(a)=(1a10). (4.2)

    corresponding to an initial infected population that amounts to 50% of the susceptible population in the age interval ]3,5[. As boundary data, we keep the choice (3.4), corresponding to a constant natality rate for the S and R populations, and to no one being infected at birth.

    As a first reference situation, we consider the case of no vaccination being dosed to anyone. The resulting integration is displayed in Figure 6. In this reference case, the total number of infected individuals is 22.77.

    Figure 6.  Integration of (2.1)–(2.2) with parameters (4.1), initial datum (4.2) and boundary inflow (3.4). From left to right: the S, I and R plots. Note the rise in the I population in the region t3 and a[1,3], a clear effect of the nonlocal source terms in (2.1) modeling the spread of the infection.

    We only remark that Figure 6 shows a typical effect of the nonlocality of the source terms in (2.1). There is an increase in the I population for, approximately, t3 and a[1,3] due to the transmission and growth of the infection, a feature which becomes relevant in regions close to (relatively) high values of the susceptible population S.

    Below we restrict the values of the percentages ηj and νi in (2.3) and (2.4) to the interval [0,0.8]. Indeed, it is reasonable to assume that not every individual can be vaccinated and that, in some cases, vaccination may fail to immunize a small percentage of susceptible individuals. For the sake of simplicity, we introduce only three vaccination ages ˉa1, ˉa2 and ˉa3 or times ˉt1, ˉt2 and ˉt3.

    We present both sample values of the costs (1.9)–(1.10) and the results of optimization procedures. The former correspond to dosing vaccines to 60% or 80% of the S population at each prescribed age for all times, or at each prescribed time at all ages. The latter are achieved through the usual descent method, see [13,§ 7.2.2], with suitable projections that restrict to values of η1,η2,η3 and ν1,ν2,ν3 in [0,0.8]. A formal justification of the use of this method is in the Lipschitz continuous dependence of the cost (1.9) or (1.10) on the various parameters, see [11].

    In the integrations below, we keep Δa=2.5103 and Δt=1.25103. However, due to computing time limitations, when applying the steepest descent method we passed to coarse meshes, where Δa=5102 and Δt=2.5102. Correspondingly, in case of the controls (2.3), optimal controls are sought in the class of locally constant functions of the type ηj(t)=20=1ηjχ[(1)/4,/4[(t). An entirely similar procedure is followed in the case of the age dependent controls (2.4). For completeness, we recall that in the present merely Lipschitz continuous setting, the steepest descent (or gradient) method needs not converge and, even if it does, its limit may well be different from a point of minimum. Therefore, in our application of this method, we keep checking that the values of the cost (do indeed) decrease, but we can not ensure that the values obtained are indeed points of minimum.

    Now, in the model (2.1)–(2.2) with parameters (4.1), initial datum (4.2) and boundary datum (3.4) we prescribe 3 vaccination ages and, correspondingly, 3 constant values of the percentage of vaccinated individuals, so that, with reference to the notation (2.3),

    ˉa1=2.0,ˉa2=4.0,ˉa3=6.0 and ηj{0.6,0.8} for j=1,2,3. (4.3)

    First, we perform the resulting 8 numerical integrations and record the resulting costs (1.5), (1.8) and (1.9). Then, these costs are compared among each other and with the costs resulting from the steepest descent method. The numerical values obtained are in Table 3.

    Table 3.  Costs resulting from the integrations of model (2.1)–(2.2) with parameters (4.1), initial datum (4.2), boundary datum (3.4) and control (2.3) with parameters (4.3).
    η1 η2 η3 J as in (1.8) N as in (1.5) J+N as in (1.9)
    0.6 0.6 0.6 14.54 7.56 22.11
    0.6 0.6 0.8 14.38 7.87 22.25
    0.6 0.8 0.6 14.18 8.02 22.20
    0.6 0.8 0.8 14.02 8.31 22.83
    0.8 0.6 0.6 13.12 9.22 22.34
    0.8 0.6 0.8 12.96 9.52 22.48
    0.8 0.8 0.6 12.80 9.60 22.40
    0.8 0.8 0.8 12.64 9.88 22.53
    see Figure 8 16.01 4.67 20.69

     | Show Table
    DownLoad: CSV

    As is to be expected, the lowest amount of infected individuals results from the highest (and most expensive) vaccination dosing, detailed in Figure 7. There, the discontinuities in S across the ages ˉaj resulting from the conditions (2.3) are clearly visible. Choosing ˉa1=2.0 as the lowest vaccination age allows, for about t4, that individuals younger that this age are infected and, at later times, the disease may well spread among juveniles.

    Figure 7.  Contour plots of the solution to model (2.1)–(2.2) with parameters (4.1), initial datum (4.2), boundary datum (3.4) and control (2.3) with ˉa1=2.0, ˉa2=4.0, ˉa3=6.0 and η1=η2=η3=0.8.

    The steepest descent procedure yields a value of the cost J+N in (1.9) lower than any of the ones obtained with constant ηj as in (4.3), see Table 3. The optimal controls found by means of the steepest descent method are shown in Figure 8 and the corresponding solutions are in Figure 9.

    Figure 8.  Optimal time dependent controls found by means of the descent method applied to minimize the cost (1.9), corresponding to the integration in Figure 9.
    Figure 9.  Contour plots of the solutions to to model (2.1)–(2.2) with parameters (4.1), initial datum (4.2), boundary datum (3.4), optimized with respect to the controls η1,η2,η3 in (2.3).

    We now follow a procedure, in a sense, analogous to that detailed in the preceding Paragraph 4.1. Choose the times and controls

    ˉt1=0.5,ˉt2=1.0,ˉt3=1.5 and νi{0.6,0.8} for i=1,2,3. (4.4)

    The costs of the integrations resulting from these choices and from the optimization are collected in Table 4. Again, the lowest number of infected individuals is obtained with the highest number of vaccinations and the corresponding solutions are plotted in Figure 10. Given that infected individuals can not turn back susceptible, it is clear that the first vaccination of 80% of all susceptible individuals greatly hinders the spreading of the disease. The later vaccinations, dosed also to the newly born individuals, block further infections, see Figure 10.

    Table 4.  Costs resulting from the integrations of model (2.1)–(2.2) with parameters (4.1), initial datum (4.2), boundary datum (3.4) and control (2.4) with parameters (4.4).
    ν1 ν2 ν3 J as in (1.8) N as in (1.7) J+N as in (1.10)
    0.6 0.6 0.6 6.41 13.05 19.47
    0.6 0.6 0.8 6.23 13.75 19.98
    0.6 0.8 0.6 6.17 13.60 19.77
    0.6 0.8 0.8 6.06 14.10 20.16
    0.8 0.6 0.6 5.78 13.89 19.66
    0.8 0.6 0.8 5.67 14.44 20.12
    0.8 0.8 0.6 5.66 14.22 19.88
    0.8 0.8 0.8 5.59 14.65 20.24
    see Figure 11 7.80 8.18 15.98

     | Show Table
    DownLoad: CSV
    Figure 10.  Contour plots of the solution to model (2.1)–(2.2) with parameters (4.1), initial datum (4.2), boundary datum (3.4) and control (2.4) with ˉt1=0.5, ˉt2=1.0, ˉt3=1.5 and ν1=ν2=ν3=0.8.

    The strategy found through the steepest descent method, shown in Figure 11, clearly takes advantage of the particular shape of the initial datum (4.2). Indeed, the two "humps" in ν2 and ν3 roughly correspond to the wave of infected individuals propagating from t=0 onward.

    Figure 11.  Plots of the optimal age dependent controls ν1,ν2,ν3 found by means of the descent method applied to minimize the cost (1.10).

    As a consequence, note that the optimal control ν3 vanishes for small ages. This is due to the fact that newborn are healthy and that the preceding vaccination campaigns succeeded in avoiding the transmission of the infection to the youngest part of the population, see Figure 12.

    Figure 12.  Contour plots of the solutions to model (2.1)–(2.2) with parameters (4.1), initial datum (4.2), boundary datum (3.4), optimized with respect to the controls ν1,ν2,ν3 in (2.4).

    With the data chosen in (4.1) and the parameters (4.2), the best vaccination strategies turn out to be the ones of the type (2.4), i.e., those where a percentage νi(a) of all susceptible individuals of age a are dosed at time ˉti, as it stems out comparing Table 3 with Table 4. Clearly, this outcome is a consequence of the particular choices in (4.1) and (4.2), but the framework based on (2.1) and (2.3) or (2.4), and its amenability to optimization procedure, is independent of these choices.

    We addressed the issue of optimizing vaccination strategies. Sample numerical integrations show various features of the solutions to the integro–differential model (2.1)–(2.2) and the effects of different choice of the controls ηj in (2.3) or νi in (2.4). Standard numerical optimization procedures, such as the steepest descent method, can be applied to single out optimal choices of the controls.

    The present framework is amenable to a variety of different applications and extensions, only a minor part of which were explicitly considered above. For instance, natality, which is represented by the boundary datum, can be assigned as a function of the present, or past, population densities, also allowing for newborn to be infected also at birth. The present age structured formulation allows to take into account the possibility that the disease is transmitted differently in different age groups. Introducing in the S, I or R populations sexual distinctions only amounts to consider more equations, the basic analytic framework in [11] remaining essentially unaltered.

    Part of this work was supported by the PRIN 2015 project Hyperbolic Systems of Conservation Laws and Fluid Dynamics: Analysis and Applications and by the GNAMPA 2018 project Conservation Laws: Hyperbolic Games, Vehicular Traffic and Fluid dynamics.

    The IBM Power Systems Academic Initiative substantially contributed to the numerical integrations.

    The authors declare there is no conflict of interest.

    [1] Orban JC, Walrave Y, Mongardon N, et al. (2017) Causes and characteristics of death in intensive care units. Anesthesiology 126: 882–889. doi: 10.1097/ALN.0000000000001612
    [2] Lovett M, Lee K, Edwards A, et al. (2009) Vascularization strategies for tissue engineering. Tissue Eng Part B Rev 15: 353–370. doi: 10.1089/ten.teb.2009.0085
    [3] Griffith CK, Miller C, Sainson RCA, et al. (2005) Diffusion limits of an in vitro thick prevascularized tissue. Tissue Eng 11: 257–266. doi: 10.1089/ten.2005.11.257
    [4] Jain RK, Au P, Tam J, et al. (2005) Engineering vascularized tissue. Nat Biotechnol 23: 821–823. doi: 10.1038/nbt0705-821
    [5] Phelps EA, García AJ (2010) Engineering more than a cell: Vascularization strategies in tissue engineering. Curr Opin Biotechnol 21: 704–709. doi: 10.1016/j.copbio.2010.06.005
    [6] Bertassoni LE, Cecconi M, Manoharan V, et al. (2014) Hydrogel bioprinted microchannel networks for vascularization of tissue engineering constructs. Lab Chip 14: 2202–2211. doi: 10.1039/C4LC00030G
    [7] Rouwkema J, Rivron NC, van Blitterswijk CA (2008) Vascularization in tissue engineering. Trends Biotechnol 26: 434–441. doi: 10.1016/j.tibtech.2008.04.009
    [8] Ikada Y (2006) Challenges in tissue engineering. J R Soc Interface 3: 589–601. doi: 10.1098/rsif.2006.0124
    [9] Griffith LG, Naughton G (2002) Tissue engineering--current challenges and expanding opportunities. Science 295: 1009–1014. doi: 10.1126/science.1069210
    [10] Kang HW, Lee SJ, Ko IK, et al. (2016) A 3D bioprinting system to produce human-scale tissue constructs with structural integrity. Nat Biotechnol 34: 312–319. doi: 10.1038/nbt.3413
    [11] Berthiaume F, Maguire TJ, Yarmush ML (2011) CH02CH19-Yarmush. Annu Rev Chem Biomol Eng 2: 403–430. doi: 10.1146/annurev-chembioeng-061010-114257
    [12] Hollister SJ (2005) Porous scaffold design for tissue engineering. Nat Mater 4: 518–524. doi: 10.1038/nmat1421
    [13] Chan BP, Leong KW (2008) Scaffolding in tissue engineering: general approaches and tissue-specific considerations. Eur Spine J 17: 467–479. doi: 10.1007/s00586-008-0745-3
    [14] Novosel EC, Kleinhans C, Kluger PJ (2011) Vascularization is the key challenge in tissue engineering. Adv Drug Deliv Rev 63: 300–311. doi: 10.1016/j.addr.2011.03.004
    [15] Zhao X, Liu L, Wang J, et al. (2016) In vitro vascularization of a combined system based on a 3D printing technique. J Tissue Eng Regen Med 10: 833–842. doi: 10.1002/term.1863
    [16] Pashneh-Tala S, MacNeil S, Claeyssens F (2015) The tissue-engineered vascular graft-past, present, and future. Tissue Eng Part B Rev 22.
    [17] O'Brien FJ (2011) Biomaterials & scaffolds for tissue engineering. Mater Today 14: 88–95. doi: 10.1016/S1369-7021(11)70058-X
    [18] He Y, Lu F (2016) Development of synthetic and natural materials for tissue engineering applications using adipose stem cells. Stem Cells Int 2016: 5786257.
    [19] Kim JJ, Hou L, Yang G, et al. (2017) Microfibrous scaffolds enhance endothelial differentiation and organization of induced pluripotent stem cells. Cell Mol Bioeng 10: 417–432. doi: 10.1007/s12195-017-0502-y
    [20] Plunkett N, O'Brien FJ (2010) IV.3. Bioreactors in tissue engineering. Stud Health Technol Inform 152: 214–230.
    [21] Kolesky DB, Homan KA, Skylar-Scott MA, et al. (2016) Three-dimensional bioprinting of thick vascularized tissues. Proc Natl Acad Sci 113: 3179–3184. doi: 10.1073/pnas.1521342113
    [22] Kolte D, McClung JA, Aronow WS (2016) Vasculogenesis and angiogenesis. Transl Res in Coron Artery Dis 49–65.
    [23] Risau W (1997) Mechanisms of angiogenesis. Nat 386: 671–674. doi: 10.1038/386671a0
    [24] Carmeliet P, Jain RK (2000) Angiogenesis in cancer and other diseases. Nat 407: 249–257. doi: 10.1038/35025220
    [25] Carmeliet P (2000) Mechanisms of angiogenesis and arteriogenesis. Nat Med 6: 389–395. doi: 10.1038/74651
    [26] Lanza RP, Langer RS, Vacanti J (2014) Principles of Tissue Engineering. Elsevier.
    [27] Gao Y (2017) Biology of Vascular Smooth Muscle: Vasoconstriction and Dilatation, Singapore: Springer Singapore.
    [28] Ardalani H, Assadi AH, Murphy WL (2014) Structure, function, and development of blood vessels: lessons for tissue engineering. In, Engineering in Translational Medicine, London: Springer London, 155–182.
    [29] Rhodin JAG (1980) Architecture of the vessel wall. In, Comprehensive Physiology, Hoboken: John Wiley & Sons, 1–31.
    [30] Kossmann CE, Palade GE (1961) Blood capillaries of the heart and other organs. Circulation 24: 368–384. doi: 10.1161/01.CIR.24.2.368
    [31] Burton AC (1954) Relation of structure to function of the tissues of the wall of blood vessels. Physiol Rev 34: 619–642. doi: 10.1152/physrev.1954.34.4.619
    [32] Yamamoto H, Ehling M, Kato K, et al. (2015) Integrin β1 controls VE-cadherin localization and blood vessel stability. Nat Commun 6: 6429. doi: 10.1038/ncomms7429
    [33] Okabe E, Todoki K, Ito H (1990) Microcirculation: function and regulation in microvasculature. In, Dynamic Aspects of Dental Pulp, Dordrecht: Springer Netherlands, pp 151–166.
    [34] Ji S, Guvendiren M (2017) Recent advances in bioink design for 3D bioprinting of tissues and organs. Front Bioeng Biotechnol 5: 23.
    [35] Hölzl K, Lin S, Tytgat L, et al. (2016) Bioink properties before, during and after 3D bioprinting. Biofabrication 8: 032002. doi: 10.1088/1758-5090/8/3/032002
    [36] Hospodiuk M, Dey M, Sosnoski D, et al. (2017) The bioink: A comprehensive review on bioprintable materials. Biotechnol Adv 35: 217–239. doi: 10.1016/j.biotechadv.2016.12.006
    [37] Smith CM, Stone AL, Parkhill RL, et al. (2004) Three-dimensional bioassembly tool for generating viable tissue-engineered constructs. Tissue Eng 10: 1566–1576. doi: 10.1089/ten.2004.10.1566
    [38] Chen CY, Barron JA, Ringeisen BR (2006) Cell patterning without chemical surface modification: cell–cell interactions between printed bovine aortic endothelial cells (BAEC) on a homogeneous cell-adherent hydrogel. Appl Surf Sci 252: 8641–8645. doi: 10.1016/j.apsusc.2005.11.088
    [39] Gao Q, He Y, Fu J, et al. (2015) Coaxial nozzle-assisted 3D bioprinting with built-in microchannels for nutrients delivery. Biomater 61: 203–215. doi: 10.1016/j.biomaterials.2015.05.031
    [40] Norotte C, Marga FS, Niklason LE, et al. (2009) Scaffold-free vascular tissue engineering using bioprinting. Biomaterials 30: 5910–5917. doi: 10.1016/j.biomaterials.2009.06.034
    [41] Zhang YS, Arneri A, Bersini S, et al. (2016) Bioprinting 3D microfibrous scaffolds for engineering endothelialized myocardium and heart-on-a-chip. Biomater 110: 45–59. doi: 10.1016/j.biomaterials.2016.09.003
    [42] Lee VK, Kim DY, Ngo H, et al. (2014) Creating perfused functional vascular channels using 3D bio-printing technology. Biomater 35: 8092–8102. doi: 10.1016/j.biomaterials.2014.05.083
    [43] Yin Yu, Ozbolat IT (2014) Tissue strands as "bioink" for scale-up organ printing. Conf Proc IEEE Eng Med Biol Soc 2014: 1428–1431.
    [44] Wilkens CA, Rivet CJ, Akentjew TL, et al. (2016) Layer-by-layer approach for a uniformed fabrication of a cell patterned vessel-like construct. Biofabrication 9: 015001. doi: 10.1088/1758-5090/9/1/015001
    [45] Cui X, Boland T (2009) Human microvasculature fabrication using thermal inkjet printing technology. Biomater 30: 6221–6227. doi: 10.1016/j.biomaterials.2009.07.056
    [46] Darland DC, Massingham LJ, Smith SR, et al. (2003) Pericyte production of cell-associated VEGF is differentiation-dependent and is associated with endothelial survival. Dev Biol 264: 275–288. doi: 10.1016/j.ydbio.2003.08.015
    [47] Korff T, Kimmina S, Martiny-Baron G, et al. (2001) Blood vessel maturation in a 3-dimensional spheroidal coculture model: direct contact with smooth muscle cells regulates endothelial cell quiescence and abrogates VEGF responsiveness. FASEB J 15: 447–457. doi: 10.1096/fj.00-0139com
    [48] Lee VK, Lanzi AM, Ngo H, et al. (2014) Generation of Multi-scale Vascular Network System Within 3D Hydrogel Using 3D Bio-printing Technology. Cell Mol Bioeng 7: 460–472. doi: 10.1007/s12195-014-0340-0
    [49] Jia W, Gungor-Ozkerim PS, Zhang YS, et al. (2016) Direct 3D bioprinting of perfusable vascular constructs using a blend bioink. Biomater 106: 58–68. doi: 10.1016/j.biomaterials.2016.07.038
    [50] Jang J, Park HJ, Kim SW, et al. (2017) 3D printed complex tissue construct using stem cell-laden decellularized extracellular matrix bioinks for cardiac repair. Biomater 112: 264–274. doi: 10.1016/j.biomaterials.2016.10.026
    [51] Jia J, Richards DJ, Pollard S, et al. (2014) Engineering alginate as bioink for bioprinting. Acta Biomater 10: 4323–4331. doi: 10.1016/j.actbio.2014.06.034
    [52] Das S, Pati F, Choi YJ, et al. (2015) Bioprintable, cell-laden silk fibroin–gelatin hydrogel supporting multilineage differentiation of stem cells for fabrication of three-dimensional tissue constructs. Acta Biomater 11: 233–246. doi: 10.1016/j.actbio.2014.09.023
    [53] Blaeser A, Duarte Campos DF, Puster U, et al. (2016) Controlling shear stress in 3D bioprinting is a key factor to balance printing resolution and stem cell integrity. Adv Healthc Mater 5: 326–333. doi: 10.1002/adhm.201500677
    [54] Hipp J, Atala A (2008) Sources of stem cells for regenerative medicine. Stem Cell Rev 4: 3–11. doi: 10.1007/s12015-008-9010-8
    [55] Wilson KD, Wu JC (2015) Induced pluripotent stem cells. JAMA 313: 1613. doi: 10.1001/jama.2015.1846
    [56] Wong CW, Chen YT, Chien CL, et al. (2018) A simple and efficient feeder-free culture system to up-scale iPSCs on polymeric material surface for use in 3D bioprinting. Mater Sci Eng C 82: 69–79. doi: 10.1016/j.msec.2017.08.050
    [57] Faulkner-Jones A, Fyfe C, Cornelissen DJ, et al. (2015) Bioprinting of human pluripotent stem cells and their directed differentiation into hepatocyte-like cells for the generation of mini-livers in 3D. Biofabrication 7: 044102. doi: 10.1088/1758-5090/7/4/044102
    [58] Ma X, Qu X, Zhu W, et al. (2016) Deterministically patterned biomimetic human iPSC-derived hepatic model via rapid 3D bioprinting. Proc Natl Acad Sci 113: 2206–2211. doi: 10.1073/pnas.1524510113
    [59] Miller JS, Stevens KR, Yang MT, et al. (2012) Rapid casting of patterned vascular networks for perfusable engineered three-dimensional tissues. Nat Mater 11: 768–774. doi: 10.1038/nmat3357
    [60] Kolesky DB, Truby RL, Gladman AS, et al. (2014) 3D bioprinting of vascularized, heterogeneous cell-laden tissue constructs. Adv Mater 26: 3124–3130. doi: 10.1002/adma.201305506
    [61] Sukmana I (2012) Microvascular guidance: a challenge to support the development of vascularised tissue engineering construct. The Sci World J 2012: 201352.
    [62] Kim JA, Kim HN, Im SK, et al. (2015) Collagen-based brain microvasculature model in vitro using three-dimensional printed template. Biomicrofluidics 9: 024115. doi: 10.1063/1.4917508
    [63] Muscari C, Giordano E, Bonafè F, et al. (2014) Strategies affording prevascularized cell-based constructs for myocardial tissue engineering. Stem Cells Int 2014: 434169.
    [64] Bogorad MI, DeStefano J, Karlsson J, et al. (2015) Review: in vitro microvessel models. Lab Chip 15: 4242–4255. doi: 10.1039/C5LC00832H
    [65] Moon JJ, West JL (2008) Vascularization of engineered tissues: approaches to promote angio-genesis in biomaterials. Curr Top Med Chem 8: 300–310. doi: 10.2174/156802608783790983
    [66] Yang P, Huang X, Shen J, et al. (2013) Development of a new pre-vascularized tissue-engineered construct using pre-differentiated rADSCs, arteriovenous vascular bundle and porous nano-hydroxyapatide-polyamide 66 scaffold. BMC Musculoskelet Disord 14: 318. doi: 10.1186/1471-2474-14-318
    [67] Zhu W, Qu X, Zhu J, et al. (2017) Direct 3D bioprinting of prevascularized tissue constructs with complex microarchitecture. Biomater 124: 106–115. doi: 10.1016/j.biomaterials.2017.01.042
    [68] Mirabella T, MacArthur JW, Cheng D, et al. (2017) 3D-printed vascular networks direct therapeutic angiogenesis in ischaemia. Nat Biomed Eng 1: 0083. doi: 10.1038/s41551-017-0083
    [69] Gelber MK, Hurst G, Comi TJ, et al. (2018) Model-guided design and characterization of a high-precision 3D printing process for carbohydrate glass. Addit Manuf 22: 38–50. doi: 10.1016/j.addma.2018.04.026
    [70] Xu C, Lee W, Dai G, et al. (2018) Highly elastic biodegradable single-network hydrogel for cell printing. ACS Appl Mater Interfaces 10: 9969–9979. doi: 10.1021/acsami.8b01294
    [71] Hinton TJ, Jallerat Q, Palchesko RN, et al. (2015) Three-dimensional printing of complex biological structures by freeform reversible embedding of suspended hydrogels. Sci Adv 1: e1500758–e1500758. doi: 10.1126/sciadv.1500758
    [72] Suntornnond R, Tan EYS, An J, et al. (2017) A highly printable and biocompatible hydrogel composite for direct printing of soft and perfusable vasculature-like structures. Sci Rep 7: 1–11. doi: 10.1038/s41598-016-0028-x
    [73] Lee V, Singh G, Trasatti JP, et al. (2014) Design and Fabrication of Human Skin by Three-Dimensional Bioprinting. Tissue Eng Part C Methods 20: 473–484. doi: 10.1089/ten.tec.2013.0335
    [74] Qilong Z, Juan W, Huanqing C, et al. (2018) Programmed shape‐morphing scaffolds enabling facile 3D endothelialization. Adv Funct Mater 0: 1801027.
    [75] Yu Y, Moncal KK, Li J, et al. (2016) Three-dimensional bioprinting using self-assembling scalable scaffold-free "tissue strands" as a new bioink. Sci Rep 6: 28714. doi: 10.1038/srep28714
    [76] Xu C, Chai W, Huang Y, et al. (2012) Scaffold-free inkjet printing of three-dimensional zigzag cellular tubes. Biotechnol Bioeng 109: 3152–3160. doi: 10.1002/bit.24591
    [77] Zhang Y, Yu Y, Chen H, et al. (2013) Characterization of printable cellular micro-fluidic channels for tissue engineering. Biofabrication 5: 25004. doi: 10.1088/1758-5082/5/2/025004
    [78] Visser CW, Kamperman T, Karbaat LP, et al. (2018) In-air microfluidics enables rapid fabrication of emulsions, suspensions, and 3D modular (bio)materials. Sci Adv 4: 1–9.
    [79] Daniel JT, Jessop ZM, Whitaker IS (2017) 3D Bioprinting for Reconstructive Surgery Techniques and Applications. Woodhead Publishing.
    [80] Mandrycky C, Wang Z, Kim K, et al. (2016) 3D bioprinting for engineering complex tissues. Biotechnol Adv 34: 422–434. doi: 10.1016/j.biotechadv.2015.12.011
    [81] Catros S, Fricain JC, Guillotin B, et al. (2011) Laser-assisted bioprinting for creating on-demand patterns of human osteoprogenitor cells and nano-hydroxyapatite. Biofabrication 3: 025001. doi: 10.1088/1758-5082/3/2/025001
    [82] Catros S, Guillotin B, Bačáková M, et al. (2011) Effect of laser energy, substrate film thickness and bioink viscosity on viability of endothelial cells printed by laser-assisted bioprinting. Appl Surf Sci 257: 5142–5147. doi: 10.1016/j.apsusc.2010.11.049
    [83] Grogan SP, Chung PH, Soman P, et al. (2013) Digital micromirror device projection printing system for meniscus tissue engineering. Acta Biomater 9: 7218–7226. doi: 10.1016/j.actbio.2013.03.020
    [84] Sorkio A, Koch L, Koivusalo L, et al. (2018) Human stem cell based corneal tissue mimicking structures using laser-assisted 3D bioprinting and functional bioinks. Biomater 171: 57–71. doi: 10.1016/j.biomaterials.2018.04.034
    [85] Gruene M, Pflaum M, Deiwick A, et al. (2011) Adipogenic differentiation of laser-printed 3D tissue grafts consisting of human adipose-derived stem cells. Biofabrication 3: 015005. doi: 10.1088/1758-5082/3/1/015005
    [86] Foyt DA, Norman MDA, Yu TTL, et al. (2018) Exploiting advanced hydrogel technologies to address key challenges in regenerative medicine. Adv Healthc Mater 7: 1700939. doi: 10.1002/adhm.201700939
    [87] Shanjani Y, Pan CC, Elomaa L, et al. (2015) A novel bioprinting method and system for forming hybrid tissue engineering constructs. Biofabrication 7: 045008. doi: 10.1088/1758-5090/7/4/045008
    [88] Aguilar JP, Lipka M, Primo GA, et al. (2018) 3D Electrophoresis-assisted lithography (3DEAL): 3D molecular printing to create functional patterns and anisotropic hydrogels. Adv Funct Mater 28: 1–10.
    [89] Hribar KC, Meggs K, Liu J, et al. (2015) Three-dimensional direct cell patterning in collagen hydrogels with near-infrared femtosecond laser. Sci Rep 5: 17203. doi: 10.1038/srep17203
    [90] Wang Z, Jin X, Dai R, et al. (2016) An ultrafast hydrogel photocrosslinking method for direct laser bioprinting. RSC Adv 6: 21099–21104. doi: 10.1039/C5RA24910D
    [91] Miri AK, Nieto D, Iglesias L, et al. (2018) Microfluidics-enabled multimaterial maskless stereolithographic bioprinting. Adv Mater 30: e1800242. doi: 10.1002/adma.201800242
    [92] Mosadegh B, Xiong G, Dunham S, et al. (2015) Current progress in 3D printing for cardiovascular tissue engineering. Biomed Mater 10: 034002. doi: 10.1088/1748-6041/10/3/034002
    [93] Ringeisen BR, Pirlo RK, Wu PK, et al. (2013) Cell and organ printing turns 15: Diverse research to commercial transitions. MRS Bull 38: 834–843. doi: 10.1557/mrs.2013.209
    [94] Duan B (2017) State-of-the-art review of 3D bioprinting for cardiovascular tissue engineering. Ann Biomed Eng 45: 195–209. doi: 10.1007/s10439-016-1607-5
    [95] Gao Q, Liu Z, Lin Z, et al. (2017) 3D bioprinting of vessel-like structures with multilevel fluidic channels. ACS Biomater Sci Eng 3: 399–408. doi: 10.1021/acsbiomaterials.6b00643
    [96] Khademhosseini A, Camci-Unal G (2018) 3D Bioprinting in Regenerative Engineering: Principles and Applications. CRC Press.
    [97] Wu Y, Lin ZY (William), Wenger AC, et al. (2018) 3D bioprinting of liver-mimetic construct with alginate/cellulose nanocrystal hybrid bioink. Bioprinting 9: 1–6. doi: 10.1016/j.bprint.2017.12.001
    [98] Visconti RP, Kasyanov V, Gentile C, et al. (2010) Towards organ printing: engineering an intra-organ branched vascular tree. Expert Opin Biol Ther 10: 409–420. doi: 10.1517/14712590903563352
    [99] Prendergast ME, Montoya G, Pereira T, et al. (2018) Microphysiological systems: automated fabrication via extrusion bioprinting. Microphysiological Syst 2: 1–16. doi: 10.21037/mps.2017.12.01
    [100] Charbe N, McCarron PA, Tambuwala MM (2017) Three-dimensional bio-printing: A new frontier in oncology research. World J Clin Oncol 8: 21–36. doi: 10.5306/wjco.v8.i1.21
    [101] Gopinathan J, Noh I (2018) Recent trends in bioinks for 3D printing. Biomater Res 22: 11. doi: 10.1186/s40824-018-0122-1
    [102] Guvendiren M, Molde J, Soares RMD, et al. (2016) Designing biomaterials for 3d printing. ACS Biomater Sci Eng 2: 1679–1693. doi: 10.1021/acsbiomaterials.6b00121
    [103] Gu Q, Tomaskovic-Crook E, Wallace GG, et al. (2017) 3D bioprinting human induced pluripotent stem cell constructs for in situ cell proliferation and successive multilineage differentiation. Adv Healthc Mater 6: 1700175. doi: 10.1002/adhm.201700175
    [104] Ong CS, Fukunishi T, Zhang H, et al. (2017) Biomaterial-free three-dimensional bioprinting of cardiac tissue using human induced pluripotent stem cell derived cardiomyocytes. Sci Rep 7: 4566. doi: 10.1038/s41598-017-05018-4
    [105] Medvedev SP, Shevchenko AI, Zakian SM (2010) Induced pluripotent stem cells: problems and advantages when applying them in regenerative medicine. Acta Nat 2: 18–28.
    [106] Wang W, Yang J, Liu H, et al. (2011) Rapid and efficient reprogramming of somatic cells to induced pluripotent stem cells by retinoic acid receptor gamma and liver receptor homolog 1. Proc Natl Acad Sci 108: 18283–18288. doi: 10.1073/pnas.1100893108
    [107] Tidball AM, Dang LT, Glenn TW, et al. (2017) Rapid generation of human genetic loss-of-function ipsc lines by simultaneous reprogramming and gene editing. Stem Cell Reports 9: 725–731. doi: 10.1016/j.stemcr.2017.07.003
    [108] Rutz AL, Hyland KE, Jakus AE, et al. (2015) A multi-material bioink method for 3D printing tunable, cell-compatible hydrogels. Adv Mater 27: 1607–1614. doi: 10.1002/adma.201405076
    [109] Vatani M, Choi JW (2017) Direct-print photopolymerization for 3D printing. Rapid Prototyp J 23: 337–343. doi: 10.1108/RPJ-11-2015-0172
    [110] Tappa K, Jammalamadaka U (2018) Novel biomaterials used in medical 3D printing techniques. J Funct Biomater 9: 17. doi: 10.3390/jfb9010017
    [111] Aljohani W, Ullah MW, Zhang X, et al. (2018) Bioprinting and its applications in tissue engineering and regenerative medicine. Int J Biol Macromol 107: 261–275. doi: 10.1016/j.ijbiomac.2017.08.171
  • This article has been cited by:

    1. M. Ouzahra, Exponential stability of nondissipative linear system in Banach space and application to unbounded bilinear systems, 2017, 109, 01676911, 53, 10.1016/j.sysconle.2017.08.010
    2. Martin Gugat, Jan Giesselmann, G. Buttazzo, E. Casas, L. de Teresa, R. Glowinski, G. Leugering, E. Trélat, X. Zhang, Boundary feedback stabilization of a semilinear model for the flow in star-shaped gas networks, 2021, 27, 1292-8119, 67, 10.1051/cocv/2021061
    3. Martin Gugat, Stephan Gerster, On the limits of stabilizability for networks of strings, 2019, 131, 01676911, 104494, 10.1016/j.sysconle.2019.104494
    4. Charlotte Rodriguez, Networks of geometrically exact beams: Well-posedness and stabilization, 2022, 12, 2156-8472, 49, 10.3934/mcrf.2021002
    5. Martin Gugat, Günter Leugering, Ke Wang, Neumann boundary feedback stabilization for a nonlinear wave equation: A strict H2-lyapunov function, 2017, 7, 2156-8499, 419, 10.3934/mcrf.2017015
    6. Eduardo Cerpa, Emmanuelle Crépeau, Claudia Moreno, On the boundary controllability of the Korteweg–de Vries equation on a star-shaped network, 2019, 0265-0754, 10.1093/imamci/dny047
    7. Enrique Zuazua, 2013, Chapter 9, 978-3-642-32159-7, 463, 10.1007/978-3-642-32160-3_9
    8. Martin Gugat, Boundary feedback stabilization of the telegraph equation: Decay rates for vanishing damping term, 2014, 66, 01676911, 72, 10.1016/j.sysconle.2014.01.007
    9. Zhong-Jie Han, Gen-Qi Xu, Output feedback stabilisation of a tree-shaped network of vibrating strings with non-collocated observation, 2011, 84, 0020-7179, 458, 10.1080/00207179.2011.561441
    10. Farhat Shel, Exponential stability of a network of elastic and thermoelastic materials, 2013, 36, 01704214, 869, 10.1002/mma.2644
    11. Abdessamad El Alami, Imad El Harraki, Ali Boutoulout, Regional Feedback Stabilization for Infinite Semilinear Systems, 2018, 24, 1079-2724, 343, 10.1007/s10883-017-9386-3
    12. Kaïs Ammari, Farhat Shel, 2022, Chapter 2, 978-3-030-86350-0, 15, 10.1007/978-3-030-86351-7_2
    13. Martin Gugat, Markus Dick, Günter Leugering, 2013, Chapter 26, 978-3-642-36061-9, 255, 10.1007/978-3-642-36062-6_26
    14. Pierre-Olivier Lamare, Antoine Girard, Christophe Prieur, Switching Rules for Stabilization of Linear Systems of Conservation Laws, 2015, 53, 0363-0129, 1599, 10.1137/140953952
    15. Abdessamad El Alami, Ali Boutoulout, 2020, Chapter 6, 978-3-030-26148-1, 67, 10.1007/978-3-030-26149-8_6
    16. Yacine Chitour, Guilherme Mazanti, Mario Sigalotti, Stability of non-autonomous difference equations with applications to transport and wave propagation on networks, 2016, 11, 1556-1801, 563, 10.3934/nhm.2016010
    17. Yacine Chitour, Swann Marx, Guilherme Mazanti, G. Buttazzo, E. Casas, L. de Teresa, R. Glowinski, G. Leugering, E. Trélat, X. Zhang, One-dimensional wave equation with set-valued boundary damping: well-posedness, asymptotic stability, and decay rates, 2021, 27, 1292-8119, 84, 10.1051/cocv/2021067
    18. Markus Dick, Martin Gugat, Gunter Leugering, 2012, Feedback stabilization of quasilinear hyperbolic systems with varying delays, 978-1-4673-2124-2, 125, 10.1109/MMAR.2012.6347931
    19. Alaa Hayek, Serge Nicaise, Zaynab Salloum, Ali Wehbe, Existence, Uniqueness and Stabilization of Solutions of a Generalized Telegraph Equation on Star Shaped Networks, 2020, 170, 0167-8019, 823, 10.1007/s10440-020-00360-8
    20. Yacine Chitour, Guilherme Mazanti, Mario Sigalotti, Persistently damped transport on a network of circles, 2016, 369, 0002-9947, 3841, 10.1090/tran/6778
    21. Markus Dick, Martin Gugat, Günter Leugering, A strict H1-Lyapunov function and feedback stabilization for the isothermal Euler equations with friction, 2011, 1, 2155-3297, 225, 10.3934/naco.2011.1.225
    22. Fatiha Alabau-Boussouira, Vincent Perrollaz, Lionel Rosier, Finite-time stabilization of a network of strings, 2015, 5, 2156-8499, 721, 10.3934/mcrf.2015.5.721
    23. Martin Gugat, Sonja Steffensen, Dynamic boundary control games with networks of strings, 2018, 24, 1292-8119, 1789, 10.1051/cocv/2017082
    24. Martin Gugat, Markus Dick, Günter Leugering, Gas Flow in Fan-Shaped Networks: Classical Solutions and Feedback Stabilization, 2011, 49, 0363-0129, 2101, 10.1137/100799824
    25. Mohamed Ouzahra, Global stabilization of semilinear systems using switching controls, 2012, 48, 00051098, 837, 10.1016/j.automatica.2012.02.018
    26. Pierre-Olivier Lamare, Antoine Girard, Christophe Prieur, 2013, Lyapunov techniques for stabilization of switched linear systems of conservation laws, 978-1-4673-5717-3, 448, 10.1109/CDC.2013.6759922
    27. Martin Gugat, Markus Dick, Time-delayed boundary feedback stabilization of the isothermal Euler equations with friction, 2011, 1, 2156-8499, 469, 10.3934/mcrf.2011.1.469
    28. Markus Dick, Martin Gugat, Michael Herty, Günter Leugering, Sonja Steffensen, Ke Wang, 2014, Chapter 31, 978-3-319-05082-9, 487, 10.1007/978-3-319-05083-6_31
    29. Markus Dick, Martin Gugat, Günter Leugering, Classical solutions and feedback stabilization for the gas flow in a sequence of pipes, 2010, 5, 1556-181X, 691, 10.3934/nhm.2010.5.691
    30. Jon Asier Bárcena-Petisco, Márcio Cavalcante, Giuseppe Maria Coclite, Nicola De Nitti, Enrique Zuazua, Control of hyperbolic and parabolic equations on networks and singular limits, 2024, 0, 2156-8472, 0, 10.3934/mcrf.2024015
    31. Martin Gugat, Stabilization of a cyclic network of strings by nodal control, 2025, 25, 1424-3199, 10.1007/s00028-024-01030-0
    32. Yacine Chitour, Paolo Mason, Mario Sigalotti, 2025, Chapter 7, 978-3-031-82981-9, 145, 10.1007/978-3-031-82982-6_7
  • Reader Comments
  • © 2018 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(6131) PDF downloads(2020) Cited by(0)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog