Genetic ablation of parietal cells in transgenic mice: a new model for analyzing cell lineage relationships in the gastric mucosa.

The gastric mucosa of mammalian stomach contains several differentiated cell types specialized for the secretion of acid, digestive enzymes, mucus, and hormones. Understanding whether each of these cell lineages is derived from a common stem cell has been a challenging problem. We have used a genetic approach to analyze the ontogeny of progenitor cells within mouse stomach. Herpes simplex virus 1 thymidine kinase was targeted to parietal cells within the gastric mucosa of transgenic mice, and parietal cells were ablated by treatment of animals with the antiherpetic drug ganciclovir. Ganciclovir treatment produced complete ablation of parietal cells, dissolution of gastric glands, and loss of chief and mucus-producing cells. Termination of drug treatment led to the reemergence of all major gastric epithelial cell types and restoration of glandular architecture. Our results imply the existence of a pluripotent stem cell for the gastric mucosa. Parietal cell ablation should provide a model for analyzing cell lineage relationships within the stomach as well as mechanisms underlying gastric injury and repair.

The effects of habitat fragmentation on the demography and population genetic structure of Uromys Caudimaculatus

Habitat fragmentation can have an impact on a wide variety of biological processes including abundance, life history strategies, mating system, inbreeding and genetic diversity levels of individual species. Although fragmented populations have received much attention, ecological and genetic responses of species to fragmentation have still not been fully resolved. The current study investigated the ecological factors that may influence the demographic and genetic structure of the giant white-tailed rat (Uromys caudimaculatus) within fragmented tropical rainforests. It is the first study to examine relationships between food resources, vegetation attributes and Uromys demography in a quantitative manner. Giant white-tailed rat densities were strongly correlated with specific suites of food resources rather than forest structure or other factors linked to fragmentation (i.e. fragment size). Several demographic parameters including the density of resident adults and juvenile recruitment showed similar patterns. Although data were limited, high quality food resources appear to initiate breeding in female Uromys. Where data were sufficient, influx of juveniles was significantly related to the density of high quality food resources that had fallen in the previous three months. Thus, availability of high quality food resources appear to be more important than either vegetation structure or fragment size in influencing giant white-tailed rat demography. These results support the suggestion that a species’ response to fragmentation can be related to their specific habitat requirements and can vary in response to local ecological conditions. In contrast to demographic data, genetic data revealed a significant negative effect of habitat fragmentation on genetic diversity and effective population size in U. caudimaculatus. All three fragments showed lower levels of allelic richness, number of private alleles and expected heterozygosity compared with the unfragmented continuous rainforest site. Populations at all sites were significantly differentiated, suggesting restricted among population gene flow. The combined effects of reduced genetic diversity, lower effective population size and restricted gene flow suggest that long-term viability of small fragmented populations may be at risk, unless effective management is employed in the future. A diverse range of genetic reproductive behaviours and sex-biased dispersal patterns were evident within U. caudimaculatus populations. Genetic paternity analyses revealed that the major mating system in U. caudimaculatus appeared to be polygyny at sites P1, P3 and C1. Evidence of genetic monogamy, however, was also found in the three fragmented sites, and was the dominant mating system in the remaining low density, small fragment (P2). High variability in reproductive skew and reproductive success was also found but was less pronounced when only resident Uromys were considered. Male body condition predicted which males sired offspring, however, neither body condition nor heterozygosity levels were accurate predictors of the number of offspring assigned to individual males or females. Genetic spatial autocorrelation analyses provided evidence for increased philopatry among females at site P1, but increased philopatry among males at site P3. This suggests that male-biased dispersal occurs at site P1 and female-biased dispersal at site P3, implying that in addition to mating systems, Uromys may also be able to adjust their dispersal behaviour to suit local ecological conditions. This study highlights the importance of examining the mechanisms that underlie population-level responses to habitat fragmentation using a combined ecological and genetic approach. The ecological data suggested that habitat quality (i.e. high quality food resources) rather than habitat quantity (i.e. fragment size) was relatively more important in influencing giant white-tailed rat demographics, at least for the populations studied here . Conversely, genetic data showed strong evidence that Uromys populations were affected adversely by habitat fragmentation and that management of isolated populations may be required for long-term viability of populations within isolated rainforest fragments.

Genetic polymorphisms in miRNAs targeting the estrogen receptor and their effect on breast cancer risk

Breast cancer is the cancer that most commonly affects women worldwide. This type of cancer is genetically complex, but is strongly linked to steroid hormone signalling systems. Because microRNAs act as translational regulators of multiple genes, including the steroid nuclear receptors, single nucleotide polymorphisms (SNPs) in microRNAs genes can have potentially wide-ranging influences on breast cancer development. Thus, this study was conducted to investigate the relationships between six SNPs (rs6977848, rs199981120, rs185641358, rs113054794, rs66461782, and rs12940701) located in four miRNA genes predicted to target the estrogen receptor (miR-148a, miR-221, miR-186, and miR-152) and breast cancer risk in Caucasian Australian women. By using high resolution melt analysis (HRM) and polymerase chain reaction- restriction fragment length polymorphism (PCR-RFLP), 487 samples including 225 controls and 262 cases were genotyped. Analysis of their genotype and allele frequencies indicated that the differences between case and control populations was not significant for rs6977848, rs66461782, and rs12940701 because their p-values are 0.81, 0.93, 0.1 which are all above the threshold value (p=0.05). Our data thus suggests that these SNPs do not affect breast cancer risk in the tested population. In addition, rs199981120, rs185641358, and rs113054794 could not be found in this population, suggesting that these SNPs do not occur in Caucasian Australians.

Symbiotic relationships : saviour siblings, family rights and biomedicine

It is now possible to combine the use of preimplantation genetic diagnosis (PGD) and tissue matching to select an IVF embryo that will, after birth, be a compatible tissue donor for an existing individual. This article analyses the ethical issues and the regulatory frameworks that intersect around the creation of tissue compatible children.

Genetic divergence and echolocation call frequency in cryptic species of Hipposideros larvatus s.l. (Chiroptera : Hipposideridae) from the Indo-Malayan region

The intermediate leaf-nosed bat (Hipposideros larvatus) is a medium-sized bat distributed throughout the Indo-Malay region. In north-east India, bats identified as H. larvatus captured at a single cave emitted echolocation calls with a bimodal distribution of peak frequencies, around either 85 kHz or 98 kHz. Individuals echolocating at 85 kHz had larger ears and longer forearms than those echolocating at 98 kHz, although no differences were detected in either wing morphology or diet, suggesting limited resource partitioning. A comparison of mitochondrial control region haplotypes of the two phonic types with individuals sampled from across the Indo-Malay range supports the hypothesis that, in India, two cryptic species are present. The Indian 98-kHz phonic bats formed a monophyletic clade with bats from all other regional populations sampled, to the exclusion of the Indian 85-kHz bats. In India, the two forms showed 12–13% sequence divergence and we propose that the name Hipposideros khasiana for bats of the 85-kHz phonic type. Bats of the 98-kHz phonic type formed a monophyletic group with bats from Myanmar, and corresponded to Hipposideros grandis, which is suggested to be a species distinct from Hipposideros larvatus. Differences in echolocation call frequency among populations did not reflect phylogenetic relationships, indicating that call frequency is a poor indicator of evolutionary history. Instead, divergence in call frequency probably occurs in allopatry, possibly augmented by character displacement on secondary contact to facilitate intraspecific communication.

A practical guide to terminology for kimberlite facies: A systematic progression from descriptive to genetic, including a pocket guide

Kimberlite terminology remains problematic because both descriptive and genetic terms are mixed together in most existing terminology schemes. In addition, many terms used in existing kimberlite terminology schemes are not used in mainstream volcanology, even though kimberlite bodies are commonly the remains of kimberlite volcanic vents and edifices. We build on our own recently published approach to kimberlite facies terminology, involving a systematic progression from descriptive to genetic. The scheme can be used for both coherent kimberlite (i.e. kimberlite that was emplaced without undergoing any fragmentation processes and therefore preserving coherent igneous textures) and fragmental kimberlites. The approach involves documentation of components, textures and assessing the degree and effects of alteration on both components and original emplacement textures. This allows a purely descriptive composite component, textural and compositional petrological rock or deposit name to be constructed first, free of any biases about emplacement setting and processes. Then important facies features such as depositional structures, contact relationships and setting are assessed, leading to a composite descriptive and genetic name for the facies or rock unit that summarises key descriptive characteristics, emplacement processes and setting. Flow charts summarising the key steps in developing a progressive descriptive to genetic terminology are provided for both coherent and fragmental facies/deposits/rock units. These can be copied and used in the field, or in conjunction with field (e.g. drill core observations) and petrographic data. Because the approach depends heavily on field scale observations, characteristics and process interpretations, only the first descriptive part is appropriate where only petrographic observations are being made. Where field scale observations are available the progression from developing descriptive to interpretative terminology can be used, especially where some petrographic data also becomes available.

A new approach to kimberlite facies terminology using a revised general approach to the nomenclature of all volcanic rocks and deposits: Descriptive to genetic

Although kimberlite pipes/bodies are usually the remains of volcanic vents, in-vent deposits, and subvolcanic intrusions, the terminology used for kimberlite rocks has largely developed independently of that used in mainstream volcanology. Existing kimberlite terminology is not descriptive and includes terms that are rarely used, used differently, and even not used at all in mainstream volcanology. In addition, kimberlite bodies are altered to varying degrees, making application of genetic terminology difficult because original components and depositional textures are commonly masked by alteration. This paper recommends an approach to the terminology for kimberlite rocks that is consistent with usage for other volcanic successions. In modern terrains the eruption and emplacement origins of deposits can often be readily deduced, but this is often not the case for old, variably altered and deformed rock successions. A staged approach is required whereby descriptive terminology is developed first, followed by application of genetic terminology once all features, including the effects of alteration on original texture and depositional features, together with contact relationships and setting, have been evaluated. Because many volcanic successions consist of both primary volcanic deposits as well as volcanic sediments, terminology must account for both possibilities.

Genetic structure of the Korean black scraper Thamnaconus modestus inferred from microsatellite marker analysis

The Korean black scraper, Thamnaconus modestus, is one of the most economically important maricultural fish species in Korea. However, the annual catch of this fish has been continuously declining over the past several decades. In this study, the genetic diversity and relationships among four wild populations and two hatchery stocks of Korean black scraper were assessed based on 16 microsatellite (MS) markers. A total of 319 different alleles were detected over all loci with an average of 19.94 alleles per locus. The hatchery stocks [mean number of alleles (N A) = 12, allelic richness (A R) = 12, expected heterozygosity (He) = 0.834] showed a slight reduction (P > 0.05) in genetic variability in comparison with wild populations (mean N A = 13.86, A R = 12.35, He = 0.844), suggesting a sufficient level of genetic variation in the hatchery populations. Similarly low levels of inbreeding and significant Hardy–Weinberg equilibrium deviations were detected in both wild and hatchery populations. The genetic subdivision among all six populations was low but significant (overall F ST = 0.008, P < 0.01). Pairwise F ST, a phylogenetic tree, and multidimensional scaling analysis suggested the existence of three geographically structured populations based on different sea basin origins, although the isolation-by-distance model was rejected. This result was corroborated by an analysis of molecular variance. This genetic differentiation may result from the co-effects of various factors, such as historical dispersal, local environment and ocean currents. These three geographical groups can be considered as independent management units. Our results show that MS markers may be suitable not only for the genetic monitoring of hatchery stocks but also for revealing the population structure of Korean black scraper populations. These results will provide critical information for breeding programs, the management of cultured stocks and the conservation of this species.

Genetics of heifer performance in 'wet' and 'dry' seasons and their relationships with steer performance in two tropical beef genotypes.

The genetics of heifer performance in tropical 'wet' and 'dry' seasons, and relationships with steer performance, were studied in Brahman (BRAH) and Tropical Composite (TCOMP) (50% Bos indicus, African Sanga or other tropically adapted Bos taurus; 50% non-tropically adapted Bos taurus) cattle of northern Australia. Data were from 2159 heifers (1027 BRAH, 1132 TCOMP), representing 54 BRAH and 51 TCOMP sires. Heifers were assessed after post-weaning 'wet' (ENDWET) and 'dry' (ENDDRY) seasons. Steers were assessed post-weaning, at feedlot entry, over a 70-day feed test, and after similar to 120-day finishing. Measures studied in both heifers and steers were liveweight (LWT), scanned rump fat, rib fat and M. longissimus area (SEMA), body condition score (CS), hip height (HH), serum insulin-like growth factor-I concentration (IGF-I), and average daily gains (ADG). Additional steer measures were scanned intra-muscular fat%, flight time, and daily (DFI) and residual feed intake (RFI). Uni- and bivariate analyses were conducted for combined genotypes and for individual genotypes. Genotype means were predicted for a subset of data involving 34 BRAH and 26 TCOMP sires. A meta-analysis of genetic correlation estimates examined how these were related to the difference between measurement environments for specific traits. There were genotype differences at the level of means, variances and genetic correlations. BRAH heifers were significantly (P < 0.05) faster-growing in the 'wet' season, slower-growing in the 'dry' season, lighter at ENDDRY, and taller and fatter with greater CS and IGF-I at both ENDWET and ENDDRY. Heritabilities were generally in the 20 to 60% range for both genotypes. Phenotypic and genetic variances, and genetic correlations, were commonly lower for BRAH. Differences were often explained by the long period of tropical adaptation of B. indicus. Genetic correlations were high between corresponding measures at ENDWET and ENDDRY, positive between fat and muscle measures in TCOMP but negative in BRAH (mean of 13 estimates 0.50 and -0.19, respectively), and approximately zero between steer feedlot ADG and heifer ADG in BRAH. Numerous genetic correlations between heifers and steers differed substantially from unity, especially in BRAH, suggesting there may be scope to select differently in the sexes where that would aid the differing roles of heifers and steers in production. Genetic correlations declined as measurement environments became more different, the rates of decline (environment sensitivity) sometimes differing with genotype. Similar measures (LWT, HH and ADG; IGF-I at ENDWET in TCOMP) were genetically correlated with steer DFI in heifers as in steers. Heifer SEMA was genetically correlated with steer feedlot RFI in BRAH (0.75 +/- 0.27 at ENDWET, 0.66 +/- 0.24 at ENDDRY). Selection to reduce steer RFI would reduce SEMA in BRAH heifers but otherwise have only small effects on heifers before their first joining.

Phylogenetic relationships of unclassified, satellitic Pasteurellaceae obtained from different species of birds as demonstrated by 16S rRNA gene sequence comparison

Avian haemophili demonstrating in vitro satellitic growth, also referred to as the V-factor or NAD requirement, have mainly been classified with Avibacterium paragallinarum (Haemophilus paragallinarum), Avibacterium avium (Pasteurella avium), Avibacterium volantium (Pasteurella volantium) and Avibacterium sp. A (Pasteurella species A). The aim of the present study was to assess the taxonomic position of 18 V-factor-requiring isolates of unclassified Haemophilus-like organisms isolated from galliforme, anseriforme, columbiforme and gruiforme birds as well as kestrels and psittacine birds including budgerigars by conventional phenotypic tests and 16S rRNA gene sequencing. All isolates shared phenotypical characteristics which allowed classification with Pasteurellaceae. Haemolysis of bovine red blood cells was negative. Haemin (X-factor) was not required for growth. Maximum-likelihood phylogenetic analysis including bootstrap analysis showed that six isolates were related to the avian 16S rRNA group and were classified as Avibacterium according to 16S rRNA sequence analysis. Surprisingly, the other 12 isolates were unrelated to Avibacterium. Two isolates were unrelated to any of the known 16S rRNA groups of Pasteurellaceae. Two isolates were related to Volucribacter of the avian 16S rRNA group. Seven isolates belonged to the Testudinis 16S rRNA group and out of these, two isolates were closely related to taxa 14 and 32 of Bisgaard, whereas four other isolates were found to form a genus-like group distantly related to taxon 40 and one isolated remained distantly related to other members of the Testudinis group. One isolate was closely related to taxon 26 (a member of Actinobacillus sensu stricto). The study documented major genetic diversity among V-factor-requiring avian isolates beyond the traditional interpretation that they only belong to Avibacterium, underlining the limited value of satellitic growth for identification of avian members of Pasteurellaceae. Our study also emphasized that these organisms will never be isolated without the use of special media satisfying the V-factor requirement.

A matter of taste : genetic and environmental influences on responses to sweetness

Diet is a major player in the maintenance of health and onset of many diseases of public health importance. The food choice is known to be largely influenced by sensory preferences. However, in many cases it is unclear whether these preferences and dietary behaviors are innate or acquired. The aim of this thesis work was to study the extent to which the individual differences in dietary responses, especially in liking for sweet taste, are influenced by genetic factors. Several traits measuring the responses to sweetness and other dietary variables were applied in four studies: in British (TwinsUK) and Finnish (FinnTwin12 and FinnTwin16) twin studies and in a Finnish migraine family study. All the subjects were adults and they participated in chemosensory measurements (taste and smell tests) and filled in food behavior questionnaires. Further, it was studied, whether the correlations among the variables are mediated by genetic or environmental factors and where in the genome the genes influencing the heritable traits are located. A study of young adult Finnish twins (FinnTwin16, n=4388) revealed that around 40% of the food use is attributable to genetic factors and that the common, childhood environment does not affect the food use even shortly after moving from the parents home. Both the family study (n=146) and the twin studies (British twins, n=663) showed that around half of the variation in the liking for sweetness is inherited. The same result was obtained both by the chemosensory measurements (heritability 41-49%) and the questionnaire variables (heritability 31-54%). By contrast, the intensity perception of sweetness or the responses to saltiness were not influenced by genetic factors. Further, a locus influencing the use-frequency of sweet foods was identified on chromosome 16p. A closer examination of the relationships among the variables based on 663 British twins revealed that several genetic and environmental correlations exist among the different measures of liking for sweetness. However, these correlations were not very strong (range 0.06-0.55) implying that the instruments used measure slightly different aspects of the phenomenon. In addition, the assessment of the associations among responses to fatty foods, dieting behaviors, and body mass index in twin populations (TwinsUK n=1027 and FinnTwin12 n=299) showed that the dieting behaviors (cognitive restraint, uncontrolled eating, and emotional eating) mediate the relationship between obesity and diet. In conclusion, the work increased the understanding of the background variables of human eating behavior. Genetic effects were shown to underlie the variation of many dietary traits, such as liking for sweet taste, use of sweet foods, and dieting behaviors. However, the responses to salty taste were shown to be mainly determined by environmental factors and thus should more easily be modifiable by dietary education, exposure, and learning than sweet taste preferences. Although additional studies are needed to characterize the genetic element located on chromosome 16 that influences the use-frequency of sweet foods, the results underline the importance of inherited factors on human eating behavior.

Genetic correlations of young bull reproductive traits and heifer puberty traits with female reproductive performance in two tropical beef genotypes in northern Australia

Genetic correlations of young bull and heifer puberty traits with measures of early and lifetime female reproductive performance were estimated in two tropical beef cattle genotypes. Heifer age at puberty was highly (r(g) = -0.71 +/- 0.11) and moderately (r(g) = -0.40 +/- 0.20) genetically correlated with pregnancy rate at first annual mating (mating 1) and lifetime annual calving rate, respectively in Brahman (BRAH). In Tropical Composite (TCOMP), heifer age at puberty was highly correlated with reproductive outcomes from the first re-breed (mating 2), mainly due to its association with lactation anoestrous interval (r(g) = 0.72 +/- 0.17). Scrotal circumference were correlated with heifer age at puberty (r(g) = -0.41 +/- 0.11 at 12 months in BRAH; -0.30 +/- 0.13 at 6 months in TCOMP) but correlations were lower with later female reproduction traits. Bull insulin-like growth factor-I was correlated with heifer age at puberty (r(g) = -0.56 +/- 0.11 in BRAH; -0.43 +/- 0.11 in TCOMP) and blood luteinising hormone concentration was moderately correlated with lactation anoestrous interval (r(g) = 0.59 +/- 0.23) in TCOMP. Semen quality traits, including mass activity, motility and percent normal sperm were genetically correlated with lactation anoestrus and female lifetime female reproductive traits in both genotypes, but the magnitudes of the relationships differed with bull age at measurement. Preputial eversion and sheath scores were genetically associated with lifetime calving and weaning rates in both genotypes. Several of the early-in-life male and female measures examined were moderately to highly genetically correlated with early and lifetime female reproduction traits and may be useful as indirect selection criteria for improving female reproduction in tropical breeds in northern Australia.

The study of genetic diversity and population structure of the Korean fleshy shrimp, Fenneropenaeus chinensis, using newly developed microsatellite markers

The fleshy shrimp, Fenneropenaeus chinensis, is the family of Penaeidae and one of the most economically important marine culture species in Korea. However, its genetic characteristics have never been studied. In this study, a total of 240 wild F. chinensis individuals were collected from four locations as follows: Narodo (NRD, n = 60), Beopseongpo (BSP, n = 60), Chaesukpo (CSP, n = 60), and Cheonsuman (CSM, n = 60). Genetic variability and the relationships among four wild F. chinensis populations were analyzed using 13 newly developed microsatellite loci. Relatively high levels of genetic variability (mean allelic richness = 16.87; mean heterozygosity = 0.845) were found among localities. Among the 52 population loci, 13 showed significant deviation from the Hardy–Weinberg equilibrium. Neighbor-joining, principal coordinate, and molecular variance analyses revealed the presence of three subpopulations (NRD, CSM, BSP and CSP), which was consistent with clustering based on genetic distance. The mean observed heterozygosity values of the NRD, CSM, BSP, and CSP populations were 0.724, 0.821, 0.814, and 0.785 over all loci, respectively. These genetic variability and differentiation results of the four wild populations can be applied for future genetic improvement using selective breeding and to design suitable management guidelines for Korean F. chinensis culture.

Genotype-Phenotype Relationships in Long QT Syndrome: Role of Mental Stress, Adrenergic Activity and a Common KCNH2 Polymorphism

Long QT syndrome is a congenital or acquired arrhythmic disorder which manifests as a prolonged QT-interval on the electrocardiogram and as a tendency to develop ventricular arrhythmias which can lead to sudden death. Arrhythmias often occur during intense exercise and/or emotional stress. The two most common subtypes of LQTS are LQT1, caused by mutations in the KCNQ1 gene and LQT2, caused by mutations in the KCNH2 gene. LQT1 and LQT2 patients exhibit arrhythmias in different types of situations: in LQT1 the trigger is usually vigorous exercise whereas in LQT2 arrhythmia results from the patient being startled from rest. It is not clear why trigger factors and clinical outcome differ from each other in the different LQTS subtypes. It is possible that stress hormones such as catecholamines may show different effects depending on the exact nature of the genetic defect, or sensitivity to catecholamines varies from subject to subject. Furthermore, it is possible that subtle genetic variants of putative modifier genes, including those coding for ion channels and hormone receptors, play a role as determinants of individual sensitivity to life-threatening arrhythmias. The present study was designed to identify some of these risk modifiers. It was found that LQT1 and LQT2 patients show an abnormal QT-adaptation to both mental and physical stress. Furthermore, as studied with epinephrine infusion experiments while the heart was paced and action potentials were measured from the right ventricular septum, LQT1 patients showed repolarization abnormalities which were related to their propensity to develop arrhythmia during intense, prolonged sympathetic tone, such as exercise. In LQT2 patients, this repolarization abnormality was noted already at rest corresponding to their arrhythmic episodes as a result of intense, sudden surges in adrenergic tone, such as fright or rage. A common KCNH2 polymorphism was found to affect KCNH2 channel function as demonstrated by in vitro experiments utilizing mammalian cells transfected with the KCNH2 potassium channel as well as QT-dynamics in vivo. Finally, the present study identified a common β-1-adrenergic receptor genotype that is related a shorter QT-interval in LQT1 patients. Also, it was discovered that compound homozygosity for two common β-adrenergic polymorphisms was related to the occurrence of symptoms in the LQT1 type of long QT syndrome. The studies demonstrate important genotype-phenotype differences between different LQTS subtypes and suggest that common modifier gene polymorphisms may affect cardiac repolarization in LQTS. It will be important in the future to prospectively study whether variant gene polymorphisms will assist in clinical risk profiling of LQTS patients.

Genetic programming metamodel for rotating beams

This paper investigates the use of Genetic Programming (GP) to create an approximate model for the non-linear relationship between flexural stiffness, length, mass per unit length and rotation speed associated with rotating beams and their natural frequencies. GP, a relatively new form of artificial intelligence, is derived from the Darwinian concept of evolution and genetics and it creates computer programs to solve problems by manipulating their tree structures. GP predicts the size and structural complexity of the empirical model by minimizing the mean square error at the specified points of input-output relationship dataset. This dataset is generated using a finite element model. The validity of the GP-generated model is tested by comparing the natural frequencies at training and at additional input data points. It is found that by using a non-dimensional stiffness, it is possible to get simple and accurate function approximation for the natural frequency. This function approximation model is then used to study the relationships between natural frequency and various influencing parameters for uniform and tapered beams. The relations obtained with GP model agree well with FEM results and can be used for preliminary design and structural optimization studies.