Multiple-trait genomic evaluation for milk yield and milk quality traits using genomic and phenotypic data in buffalo in Brazil

The use of markers distributed all long the genome may increase the accuracy of the predicted additive genetic value of young animals that are candidates to be selected as reproducers. In commercial herds, due to the cost of genotyping, only some animals are genotyped and procedures, divided in two or three steps, are done in order to include these genomic data in genetic evaluation. However, genomic evaluation may be calculated using one unified step that combines phenotypic data, pedigree and genomics. The aim of the study was to compare a multiple-trait model using only pedigree information with another using pedigree and genomic data. In this study, 9,318 lactations from 3061 buffaloes were used, 384 buffaloes were genotyped using a Illumina bovine chip (Illumina Infinium (R) bovineHD BeadChip). Seven traits were analyzed milk yield (MY), fat yield (FY), protein yield (PY), lactose yield (LY), fat percentage (F%), protein percentage (P%) and somatic cell score (SCSt). Two analyses were done: one using phenotypic and pedigree information (matrix A) and in the other using a matrix based in pedigree and genomic information (one step, matrix H). The (co) variance components were estimated using multiple-trait analysis by Bayesian inference method, applying an animal model, through Gibbs sampling. The model included the fixed effects of contemporary groups (herd-year-calving season), number of milking (2 levels), and age of buffalo at calving as (co) variable (quadratic and linear effect). The additive genetic, permanent environmental, and residual effects were included as random effects in the model. The heritability estimates using matrix A were 0.25, 0.22, 0.26, 0.17, 0.37, 0.42 and 0.26 and using matrix H were 0.25, 0.24, 0.26, 0.18, 0.38, 0.46 and 0.26 for MY, FY, PY, LY, % F, % P and SCCt, respectively. The estimates of the additive genetic effect for the traits were similar in both analyses, but the accuracy were bigger using matrix H (superior to 15% for traits studied). The heritability estimates were moderated indicating genetic gain under selection. The use of genomic information in the analyses increases the accuracy. It permits a better estimation of the additive genetic value of the animals.

Genetic parameter estimates for feed efficiency and dry matter intake and their association with growth and carcass traits in Nellore cattle

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Association between ACTA1 candidate gene and performance, organs and carcass traits in broilers

This study investigates the genetic association of the SNP present in the ACTA1 gene with performance traits, organs and carcass of broilers to help marker-assisted selection of a paternal broiler line (TT) from EMBRAPA Swine and Poultry, Brazil. Genetic and phenotypic data of 1,400 broilers for 68 traits related to body performance, organ weights, weight of carcass parts, and yields as a percentage of organs and carcass parts were used. The maximum likelihood method, considering 4 analytical models, was used to analyze the genetic association between the SNP and these important economic traits. The association analysis was performed using a mixed animal model including the random effect of the animal (polygenic), and the fixed effects of sex (2 levels), hatch (5 levels) and SNP (3 levels), besides the random error. The traits significantly associated (P < 0.05) with the SNP were analyzed, along with body weight at 42 days of age (BW42), by the restricted maximum likelihood method using the multi-trait animal model to estimate genetic parameters. The analysis included the residual and additive genetic random effects and the sex-hatch fixed effect. The additive effects of the SNP were associated with breast meat (BMY), liver yield (LIVY), body weight at 35 days of age (BW35); drumstick skin (DSW), drumstick (DW) and breast (BW) weights. The heritability estimates for these traits, in addition to BW42, ranged from 0.24 ± 0.06 to 0.45 ± 0.08 for LIVY and BW35, respectively. The genetic correlation ranged from 0.02 ± 0.18 for LIVY and BMY to 0.97 ± 0.01 for BW35 and BW42. Based on the results of this study, it can be concluded that ACTA1 gene is associated with performance traits BW35, LIV and BMY, DW, BW and DW adjusted for body weight at 42 days of age. Therefore, the ACTA1 gene is an important molecular marker that could be used together with others already described to increase the economically important traits in broilers.

Comparison of wing geometry data and genetic data for assessing the population structure of Aedes aegypti

Aedes aegypti is the most important vector of dengue viruses in tropical and subtropical regions. Because vaccines are still under development, dengue prevention depends primarily on vector control. Population genetics is a common approach in research involving Ae. aegypti. In the context of medical entomology, wing morphometric analysis has been proposed as a strong and low-cost complementary tool for investigating population structure. Therefore, we comparatively evaluated the genetic and phenotypic variability of population samples of Ae. aegypti from four sampling sites in the metropolitan area of Sao Paulo city, Brazil. The distances between the sites ranged from 7.1 to 50 km. This area, where knowledge on the population genetics of this mosquito is incipient, was chosen due to the thousands of dengue cases registered yearly. The analysed loci were polymorphic, and they revealed population structure (global F-ST = 0.062; p < 0.05) and low levels of gene flow (Nm = 0.47) between the four locations. Principal component and discriminant analyses of wing shape variables (18 landmarks) demonstrated that wing polymorphisms were only slightly more common between populations than within populations. Whereas microsatellites allowed for geographic differentiation, wing geometry failed to distinguish the samples. These data suggest that microevolution in this species may affect genetic and morphological characters to different degrees. In this case, wing shape was not validated as a marker for assessing population structure. According to the interpretation of a previous report, the wing shape of Ae. aegypti does not vary significantly because it is stabilised by selective pressure. (C) 2011 Elsevier B.V. All rights reserved.

Molecular phylogeny and phenotypic variability of clinical and environmental strains of Aspergillus flavus

Aspergillus flavus is the second most common cause of aspergillosis infection in immunocompromised patients and is responsible for the production of aflatoxins. Little is known about the population structure of A. flavus, although recent molecular and phenotypic data seem to demonstrate that different genetic lineages exist within this species. The aim of this study was to carry out a morphological, physiological, and molecular analysis of a set of clinical and environmental isolates to determine whether this variability is due to species divergence or intraspecific diversity, and to assess whether the clinical isolates form a separate group. The amdS and omtA genes were more phylogenetically informative than the other tested genes and their combined analysis inferred three main clades, with no clear distinction between clinical and environmental isolates. No important morphological and physiological differences were found between the members of the different clades, with the exception of the assimilation of D-glucosamine, which differentiates the members of the clade II from the others. (C) 2012 The British Mycological Society. Published by Elsevier Ltd. All rights reserved.

Genetic and environmental influences on blood pressure and physical activity: a study of nuclear families from Muzambinho, Brazil

Blood pressure (BP) and physical activity (PA) levels are inversely associated. Since genetic factors account for the observed variation in each of these traits, it is possible that part of their association may be related to common genetic and/or environmental influences. Thus, this study was designed to estimate the genetic and environmental correlations of BP and PA phenotypes in nuclear families from Muzambinho, Brazil. Families including 236 offspring (6 to 24 years) and their 82 fathers and 122 mothers (24 to 65 years) were evaluated. BP was measured, and total PA (TPA) was assessed by an interview (commuting, occupational, leisure time, and school time PA). Quantitative genetic modeling was used to estimate maximal heritability (h²), and genetic and environmental correlations. Heritability was significant for all phenotypes (systolic BP: h² = 0.37 ± 0.10, P < 0.05; diastolic BP: h² = 0.39 ± 0.09, P < 0.05; TPA: h² = 0.24 ± 0.09, P < 0.05). Significant genetic (r g) and environmental (r e) correlations were detected between systolic and diastolic BP (r g = 0.67 ± 0.12 and r e = 0.48 ± 0.08, P < 0.05). Genetic correlations between BP and TPA were not significant, while a tendency to an environmental cross-trait correlation was found between diastolic BP and TPA (r e = -0.18 ± 0.09, P = 0.057). In conclusion, BP and PA are under genetic influences. Systolic and diastolic BP share common genes and environmental influences. Diastolic BP and TPA are probably under similar environmental influences.

Genetic and environmental factors associated with the risk of cognitive decline and dementia

AD is the most common age related neurodegenerative disease in the industrialized world. Clinically AD is defined as a progressing decline of cognitive functions. Neuropathologically, AD is characterized by the aggregation of b-amyloid (Ab) peptide in the form of extracellular senile plaques, and hyperphosphorlylated tau protein in the form of intracellular neurofibrillary tangles. These neuropathological hallmarks are often accompanied by abundant microvascular damage and pronounced inflammation of the affected brain regions. In this thesis we investigated several aspects of AD focusing on the genetic aspect. We confirmed that Alpha 1 antichymotrypsin (ACT), an acute phase protein, was associated to AD subjects, being plasma levels higher in AD cases than controls. In addition, in a GWA study we demonstrated that two different gene, Clusterin and CR1 were strongly associated to AD. A single gene association not explain such a complex disease like AD. The goal should be to created a network of genetic, phenotypic and clinical data associated to AD. We used a new algorithm, the ANNs, aimed to map variables and search for connectivity among variables. We found specific variables associated to AD like cholesterol levels, the presence of variation in HMGCR enzyme and the age. Other factors such as the BMI, the amount of HDL and blood folate levels were also associated with AD. Pathogen infections, above all viral infections, have been previously associated to AD. The hypothesis suggests that virus and in particular herpes virus could enter the brain when an individual becomes older, perhaps because of a decline in the immune system. Our new hypothesis is that the presence of SNPs in our GWA gene study results in a genetic signature that might affect individual brain susceptibility to infection by herpes virus family during aging.

Comparative DNA sequence and methylation analyses of orthologous genes in humans and non-human primates: Genetic and epigenetic footprints of evolution

Welche genetische Unterschiede machen uns verschieden von unseren nächsten Verwandten, den Schimpansen, und andererseits so ähnlich zu den Schimpansen? Was wir untersuchen und auch verstehen wollen, ist die komplexe Beziehung zwischen den multiplen genetischen und epigenetischen Unterschieden, deren Interaktion mit diversen Umwelt- und Kulturfaktoren in den beobachteten phänotypischen Unterschieden resultieren. Um aufzuklären, ob chromosomale Rearrangements zur Divergenz zwischen Mensch und Schimpanse beigetragen haben und welche selektiven Kräfte ihre Evolution geprägt haben, habe ich die kodierenden Sequenzen von 2 Mb umfassenden, die perizentrischen Inversionsbruchpunkte flankierenden Regionen auf den Chromosomen 1, 4, 5, 9, 12, 17 und 18 untersucht. Als Kontrolle dienten dabei 4 Mb umfassende kollineare Regionen auf den rearrangierten Chromosomen, welche mindestens 10 Mb von den Bruchpunktregionen entfernt lagen. Dabei konnte ich in den Bruchpunkten flankierenden Regionen im Vergleich zu den Kontrollregionen keine höhere Proteinevolutionsrate feststellen. Meine Ergebnisse unterstützen nicht die chromosomale Speziationshypothese für Mensch und Schimpanse, da der Anteil der positiv selektierten Gene (5,1% in den Bruchpunkten flankierenden Regionen und 7% in den Kontrollregionen) in beiden Regionen ähnlich war. Durch den Vergleich der Anzahl der positiv und negativ selektierten Gene per Chromosom konnte ich feststellen, dass Chromosom 9 die meisten und Chromosom 5 die wenigsten positiv selektierten Gene in den Bruchpunkt flankierenden Regionen und Kontrollregionen enthalten. Die Anzahl der negativ selektierten Gene (68) war dabei viel höher als die Anzahl der positiv selektierten Gene (17). Eine bioinformatische Analyse von publizierten Microarray-Expressionsdaten (Affymetrix Chip U95 und U133v2) ergab 31 Gene, die zwischen Mensch und Schimpanse differentiell exprimiert sind. Durch Untersuchung des dN/dS-Verhältnisses dieser 31 Gene konnte ich 7 Gene als negativ selektiert und nur 1 Gen als positiv selektiert identifizieren. Dieser Befund steht im Einklang mit dem Konzept, dass Genexpressionslevel unter stabilisierender Selektion evolvieren. Die meisten positiv selektierten Gene spielen überdies eine Rolle bei der Fortpflanzung. Viele dieser Speziesunterschiede resultieren eher aus Änderungen in der Genregulation als aus strukturellen Änderungen der Genprodukte. Man nimmt an, dass die meisten Unterschiede in der Genregulation sich auf transkriptioneller Ebene manifestieren. Im Rahmen dieser Arbeit wurden die Unterschiede in der DNA-Methylierung zwischen Mensch und Schimpanse untersucht. Dazu wurden die Methylierungsmuster der Promotor-CpG-Inseln von 12 Genen im Cortex von Menschen und Schimpansen mittels klassischer Bisulfit-Sequenzierung und Bisulfit-Pyrosequenzierung analysiert. Die Kandidatengene wurden wegen ihrer differentiellen Expressionsmuster zwischen Mensch und Schimpanse sowie wegen Ihrer Assoziation mit menschlichen Krankheiten oder dem genomischen Imprinting ausgewählt. Mit Ausnahme einiger individueller Positionen zeigte die Mehrzahl der analysierten Gene keine hohe intra- oder interspezifische Variation der DNA-Methylierung zwischen den beiden Spezies. Nur bei einem Gen, CCRK, waren deutliche intraspezifische und interspezifische Unterschiede im Grad der DNA-Methylierung festzustellen. Die differentiell methylierten CpG-Positionen lagen innerhalb eines repetitiven Alu-Sg1-Elements. Die Untersuchung des CCRK-Gens liefert eine umfassende Analyse der intra- und interspezifischen Variabilität der DNA-Methylierung einer Alu-Insertion in eine regulatorische Region. Die beobachteten Speziesunterschiede deuten darauf hin, dass die Methylierungsmuster des CCRK-Gens wahrscheinlich in Adaption an spezifische Anforderungen zur Feinabstimmung der CCRK-Regulation unter positiver Selektion evolvieren. Der Promotor des CCRK-Gens ist anfällig für epigenetische Modifikationen durch DNA-Methylierung, welche zu komplexen Transkriptionsmustern führen können. Durch ihre genomische Mobilität, ihren hohen CpG-Anteil und ihren Einfluss auf die Genexpression sind Alu-Insertionen exzellente Kandidaten für die Förderung von Veränderungen während der Entwicklungsregulation von Primatengenen. Der Vergleich der intra- und interspezifischen Methylierung von spezifischen Alu-Insertionen in anderen Genen und Geweben stellt eine erfolgversprechende Strategie dar.

The elderly health status and its correlation with ageing biomarkers: the European Project Mark-Age

According to the latest statistics projections formulated by Eurostat, the proportion of elderly EU-27’s population aged over 65 years old is predicted to increase from 17.5 % in 2011 to 29.5 % by 2060. This "population explosion" makes extremely important to identify the different genetic and molecular mechanisms which underpin the morbidity and mortality along with new strategies able to counteract or slow down its progress. In this scenario fits the European Project MARK-AGE whose aim was to identify a robust set of biomarkers of human ageing able to discriminate between chronological and biological ageing and to derive a model for healthy ageing through the analysis of three populations from different European countries, supposed to be characterized by different ageing rate: 1. Subjects representing the “Normal” or “Physiological” aging. 2. Subjects representing the “successful” or “decelerate” aging 3. Subjects representing the “accelerated” aging. The aim of this work was to recruit and characterize volunteers, to perform an accurate analysis of the health status of elderly recruited subjects (60-79 years) verifying any possible dissimilarity in their aging trajectories, to identify a set of robust ageing biomarkers and investigate possible correlations between ageing biomarkers and health status of recruited volunteers. The model proposed by MARK-AGE Project regarding different ageing trajectories has been confirmed and several ageing biomarkers have been identified.

Genetic and hormonal factors modulate spreading depression and transient hemiparesis in mouse models of familial hemiplegic migraine type 1.

Familial hemiplegic migraine type 1 (FHM1) is an autosomal dominant subtype of migraine with aura that is associated with hemiparesis. As with other types of migraine, it affects women more frequently than men. FHM1 is caused by mutations in the CACNA1A gene, which encodes the alpha1A subunit of Cav2.1 channels; the R192Q mutation in CACNA1A causes a mild form of FHM1, whereas the S218L mutation causes a severe, often lethal phenotype. Spreading depression (SD), a slowly propagating neuronal and glial cell depolarization that leads to depression of neuronal activity, is the most likely cause of migraine aura. Here, we have shown that transgenic mice expressing R192Q or S218L FHM1 mutations have increased SD frequency and propagation speed; enhanced corticostriatal propagation; and, similar to the human FHM1 phenotype, more severe and prolonged post-SD neurological deficits. The susceptibility to SD and neurological deficits is affected by allele dosage and is higher in S218L than R192Q mutants. Further, female S218L and R192Q mutant mice were more susceptible to SD and neurological deficits than males. This sex difference was abrogated by ovariectomy and senescence and was partially restored by estrogen replacement, implicating ovarian hormones in the observed sex differences in humans with FHM1. These findings demonstrate that genetic and hormonal factors modulate susceptibility to SD and neurological deficits in FHM1 mutant mice, providing a potential mechanism for the phenotypic diversity of human migraine and aura.

Ecological speciation and phenotypic plasticity affect ecosystems

Phenotypic differences among closely related populations and species can cause contrasting effects on ecosystems; however, it is unknown whether such effects result from genetic divergence, phenotypic plasticity, or both. To test this, we reared sympatric limnetic and benthic species of whitefish from a young adaptive radiation in a common garden, where the benthic species was raised on two distinct food types. We then used these fish in a mesocosm experiment to test for contrasting ecosystem effects of closely related species and of plastically induced differences within a species. We found that strong contrasting ecosystem effects resulted more frequently from genetic divergence, although they were not stronger overall than those resulting from phenotypic plasticity. Overall, our results provide evidence that genetically based differences among closely related species that evolved during a young adaptive radiation can affect ecosystems, and that phenotypic plasticity can modify the ecosystem effects of such species.

Genetic and environmental sources of covariance between reading tests used in neuropsychological assessment and IQ subtests

In this study, we examined genetic and environmental influences on covariation among two reading tests used in neuropsychological assessment (Cambridge Contextual Reading Test [CCRT], [Beardsall, L., and Huppert, F. A. ( 1994). J. Clin. Exp. Neuropsychol. 16: 232 - 242], Schonell Graded Word Reading Test [SGWRT], [ Schonell, F. J., and Schonell, P. E. ( 1960). Diagnostic and attainment testing. Edinburgh: Oliver and Boyd.]) and among a selection of IQ subtests from the Multidimensional Aptitude Battery (MAB), [Jackson, D. N. (1984). Multidimensional aptitude battery, Ontario: Research Psychologists Press.] and the Wechsler Adult Intelligence Scale-Revised (WAIS-R) [Wechsler, D. (1981). Manual for the Wechsler Adult Intelligence Scale-Revised (WAIS-R). San Antonio: The Psychological Corporation]. Participants were 225 monozygotic and 275 dizygotic twin pairs aged from 15 years to 18 years ( mean, 16 years). For Verbal IQ subtests, phenotypic correlations with the reading tests ranged from 0.44 to 0.65. For Performance IQ subtests, phenotypic correlations with the reading tests ranged from 0.23 to 0.34. Results of Structural Equation Modeling (SEM) supported a model with one genetic General factor and three genetic group factors ( Verbal, Performance, Reading). Reading performance was influenced by the genetic General factor ( accounting for 13% and 20% of the variance for the CCRT and SGWRT, respectively), the genetic Verbal factor ( explaining 17% and 19% of variance for the CCRT and SGWRT), and the genetic Reading factor ( explaining 21% of the variance for both the CCRT and SGWRT). A common environment factor accounted for 25% and 14% of the CCRT and SGWRT variance, respectively. Genetic influences accounted for more than half of the phenotypic covariance between the reading tests and each of the IQ subtests. The heritabilities of the CCRT and SGWRT were 0.54 and 0.65, respectively. Observable covariance between reading assessments used by neuropsychologists to estimate IQ and IQ subtests appears to be largely due to genetic effects.

Genetic and economic analysis of a targeted marker-assisted wheat breeding strategy

The advent of molecular markers as a tool to aid selection has provided plant breeders with the opportunity to rapidly deliver superior genetic solutions to problems in agricultural production systems. However, a major constraint to the implementation of marker-assisted selection (MAS) in pragmatic breeding programs in the past has been the perceived high relative cost of MAS compared to conventional phenotypic selection. In this paper, computer simulation was used to design a genetically effective and economically efficient marker-assisted breeding strategy aimed at a specific outcome. Under investigation was a strategy involving the integration of both restricted backcrossing and doubled haploid (DH) technology. The point at which molecular markers are applied in a selection strategy can be critical to the effectiveness and cost efficiency of that strategy. The application of molecular markers was considered at three phases in the strategy: allele enrichment in the BC1F1 population, gene selection at the haploid stage and the selection for recurrent parent background of DHs prior to field testing. Overall, incorporating MAS at all three stages was the most effective, in terms of delivering a high frequency of desired outcomes and at combining the selected favourable rust resistance, end use quality and grain yield alleles. However, when costs were included in the model the combination of MAS at the BC1F1 and haploid stage was identified as the optimal strategy. A detailed economic analysis showed that incorporation of marker selection at these two stages not only increased genetic gain over the phenotypic alternative but actually reduced the over all cost by 40%.

Maternal alcohol use disorder and offspring ADHD: disentangling genetic and environmental effects using a children-of-twins design

Background. Children of alcoholics are significantly more likely to experience high-risk environmental exposures, including prenatal substance exposure, and are more likely to exhibit externalizing problems [e.g. attention deficit hyperactivity disorder (ADHD)]. While there is evidence that genetic influences and prenatal nicotine and/or alcohol exposure play separate roles in determining risk of ADHD, little has been done on determining the joint roles that genetic risk associated with maternal alcohol use disorder (AUD) and prenatal risk factors play in determining risk of ADHD. Method. Using a children-of-twins design, diagnostic telephone interview data from high-risk families (female monozygotic and dizygotic twins concordant or discordant for AUD as parents) and control families targeted from a large Australian twin cohort were analyzed using logistic regression models. Results. Offspring of twins with a history of AUD, as well as offspring of non-AUD monozygotic twins whose co-twin had AUD, were significantly more likely to exhibit ADHD than offspring of controls. This pattern is consistent with a genetic explanation for the association between maternal AUD and increased offspring risk of ADHD. Adjustment for prenatal smoking, which remained significantly predictive, did not remove the significant genetic association between maternal AUD and offspring ADHD. Conclusions. While maternal smoking during pregnancy probably contributes to the association between maternal AUD and offspring ADHD risk, the evidence for a significant genetic correlation suggests: (i) pleiotropic genetic effects, with some genes that influence risk of AUD also influencing vulnerability to ADHD; or (ii) ADHD is a direct risk-factor for AUD.

Meta-analysis of the population and phenotypic expression of CYP2C9/2C19 polymorphism on drug metabolism in different ethnicities

The term "pharmacogenetics" has been defined as the scientific study of inherited factors that affect the human drug response. Many pharmacogenetie studies have been published since 1995 and have focussed on the principal enzyme family involved in drug metabolism, the cytochrome P450 family, particularly cytochrome P4502C9 and 2C19. In order to investigate the pharmacogenetic aspect of pharmacotherapy, the relevant studies describing the association of pharmacogenetic factor(s) in drug responses must be retrieved from existing literature using a systematic review approach. In addition, the estimation of variant allele prevalence for the gene under study between different ethnic populations is important for pharmacogenetic studies. In this thesis, the prevalence of CYP2C9/2C19 alleles between different ethnicities has been estimated through meta-analysis and the population genetic principle. The clinical outcome of CYP2C9/2C19 allelic variation on the pharmacotherapy of epilepsy has been investigated; although many new antiepileptic drugs have been launched into the market, carbamazepine, phenobarbital and phenytoin are still the major agents in the pharmacotherapy of epilepsy. Therefore, phenytoin was chosen as a model AED and the effect of CYP2C9/2C19 genetic polymorphism on phenytoin metabolism was further examined.An estimation of the allele prevalence was undertaken for three CYP2C9/2C19 alleles respectively using a meta-analysis of studies that fit the Hardy-Weinberg equilibrium. The prevalence of CYP2C9*1 is approximately 81%, 96%, 97% and 94% in Caucasian, Chinese, Japanese, African populations respectively; the pooled prevalence of CYP2C19*1 is about 86%, 57%, 58% and 85% in these ethnic populations respectively. However, the studies of association between CYP2C9/2C19 polymorphism and phenytoin metabolism failed to achieve any qualitative or quantitative conclusion. Therefore, mephenytoin metabolism was examined as a probe drug for association between CYP2C19 polymorphism and mephenytoin metabolic ratio. Similarly, analysis of association between CYP2C9 polymorphism and warfarin dose requirement was undertaken.It was confirmed that subjects carrying two mutated CYP2C19 alleles have higher S/R mephenytoin ratio due to deficient CYP2C19 enzyme activity. The studies of warfarin and CYP2C9 polymorphism did not provide a conclusive result due to poor comparability between studies.The genetic polymorphism of drug metabolism enzymes has been studied extensively, however other genetic factors, such as multiple drug resistance genes (MDR) and genes encoding ion channels, which may contribute to variability in function of drug transporters and targets, require more attention in future pharmacogenetic studies of antiepileptic drugs.