[nodec1 facec1 elemc1]=meshacylinder([80 80 100], [110 140 100],2);
[nodec2 facec2 elemc2]=meshacylinder([80 150 100], [110 150 100],2);
elemc1(:,5)=5;
elemc2(:,5)=5;
% merge the mesh of cylinder and brain
[tempNode tempElem]=mergemesh(node,elem,nodec1,elemc1);
[nNode nElem]=mergemesh(tempNode,tempElem,nodec2,elemc2);
[tempNode tempFace]=mergemesh(node,face,nodec1,facec1);
[nNode nFace]=mergemesh(tempNode,tempFace,nodec2,facec2);
The above approach is the one given in the example of mergemesh code of iso2mesh.