Unexplained length-dependence of flexural rigidity and Young's modulus of microtubules is studied using an orthotropic elastic shell model. It is showed that vibration frequencies and buckling load predicted by the accurate orthotropic shell model are much lower than that given by the approximate isotropic beam model for shorter microtubules, although the two models give almost identical results for sufficiently long microtubules. It is this inaccuracy of the isotropic beam model used by all previous researchers that leads to reported lower flexural rigidity and Young's modulus for shorter microtubules. In particular, much lower shear modulus and circumferential Young's modulus, which only weaken flexural rigidity of shorter microtubules, are responsible for the observed length-dependence of the flexural rigidity. These results confirm that longitudinal Young's modulus of microtubules is length-in dependent, and the observed length-dependence of the flexural rigidity and Young's modulus is a result of strongly anisotropic elastic properties of microtubules which have a length-dependent weakening effect on flexural rigidity of shorter microtubules. (c) 2006 Elsevier Inc. All rights reserved.