Higher frequency of genetic variants conferring increased risk for ADRs for commonly used drugs treating cancer, AIDS and tuberculosis in persons of African descent

There is established clinical evidence for differences in drug response, cure rates and survival outcomes between different ethnic populations, but the causes are poorly understood. Differences in frequencies of functional genetic variants in key drug response and metabolism genes may significantly influence drug response differences in different populations. To assess this, we genotyped 1330 individuals of African (n=372) and European (n=958) descent for 4535 single-nucleotide polymorphisms in 350 key drug absorption, distribution, metabolism, elimination and toxicity genes. Important and remarkable differences in the distribution of genetic variants were observed between Africans and Europeans and among the African populations. These could translate into significant differences in drug efficacy and safety profiles, and also in the required dose to achieve the desired therapeutic effect in different populations. Our data points to the need for population-specific genetic variation in personalizing medicine and care.

Inherited and noninherited sources of variation in metastasis and the growth rate of tumors: Implications for tumor evolution and therapy

I have undertaken measurements of the genetic (or inherited) and nongenetic (or noninherited) components of the variability of metastasis formation and tumor diameter doubling time in more than 100 metastatic lines from each of three murine tumors (sarcoma SANH, sarcoma SA4020, and hepatocarcinoma HCA-I) syngeneic to C3Hf/Kam mice. These lines were isolated twice from lung metastases and analysed immediately thereafter to obtain the variance to spontaneous lung metastasis and tumor diameter doubling time. Additional studies utilized cells obtained from within 4 passages of isolation. Under the assumption that no genetic differences in metastasis formation or diameter doubling time existed among the cells of a given line, the variance within a line would estimate nongenetic variation. The variability derived from differences between lines would represent genetic origin. The estimates of the genetic contribution to the variation of metastasis and tumor diameter doubling time were significantly greater than zero, but only in the metastatic lines of tumor SANH was genetic variation the major source of metastatic variability (contributing 53% of the variability). In the tumor cell lines of SA4020 and HCA-I, however, the contribution of nongenetic factors predominated over genetic factors in the variability of the number of metastasis and tumor diameter doubling time. A number of other parameters examined, such as DNA content, karyotype, and selection and variance analysis with passage in vivo, indicated that genetic differences existed within the cell lines and that these differences were probably created by genetic instability. The mean metastatic propensity of the lines may have increased somewhat during their isolation and isotransplantation, but the variance was only slightly affected, if at all. Analysis of the DNA profiles of the metastatic lines of SA4020 and HCA-I revealed differences between these lines and their primary parent tumors, but not among the SANH lines and their parent tumor. Furthermore, there was a direct correlation between the extent of genetic influence on metastasis formation and the ability of the tumor cells to develop resistance to cisplatinum. Thus although nongenetic factors might predominate in contributing to metastasis formation, it is probably genetic variation and genetic instability that cause the progression of tumor cells to a more metastatic phenotype and leads to the emergence of drug resistance. ^

CYP3A5*3 and POR*28 genetic variants influence the required dose of tacrolimus in heart transplant recipients

BACKGROUND After heart transplantation (HTx), the interindividual pharmacokinetic variability of immunosuppressive drugs represents a major therapeutic challenge due to the narrow therapeutic window between over-immunosuppression causing toxicity and under-immunosuppression leading to graft rejection. Although genetic polymorphisms have been shown to influence pharmacokinetics of immunosuppressants, data in the context of HTx are scarce. We thus assessed the role of genetic variation in CYP3A4, CYP3A5, POR, NR1I2, and ABCB1 acting jointly in immunosuppressive drug pathways in tacrolimus (TAC) and ciclosporin (CSA) dose requirement in HTx recipients. METHODS Associations between 7 functional genetic variants and blood dose-adjusted trough (C0) concentrations of TAC and CSA at 1, 3, 6, and 12 months after HTx were evaluated in cohorts of 52 and 45 patients, respectively. RESULTS Compared with CYP3A5 nonexpressors (*3/*3 genotype), CYP3A5 expressors (*1/*3 or *1/*1 genotype) required around 2.2- to 2.6-fold higher daily TAC doses to reach the targeted C0 concentration at all studied time points (P ≤ 0.003). Additionally, the POR*28 variant carriers showed higher dose-adjusted TAC-C0 concentrations at all time points resulting in significant differences at 3 (P = 0.025) and 6 months (P = 0.047) after HTx. No significant associations were observed between the genetic variants and the CSA dose requirement. CONCLUSIONS The CYP3A5*3 variant has a major influence on the required TAC dose in HTx recipients, whereas the POR*28 may additionally contribute to the observed variability. These results support the importance of genetic markers in TAC dose optimization after HTx.

Disentangling the role of phenotypic plasticity and genetic divergence in contemporary ecotype formation during a biological invasion

The occurrence of contemporary ecotype formation through adaptive divergence of populations within the range of an invasive species typically requires standing genetic variation but can be facilitated by phenotypic plasticity. The relative contributions of both of these to adaptive trait differentiation have rarely been simultaneously quantified in recently diverging vertebrate populations. Here we study a case of intraspecific divergence into distinct lake and stream ecotypes of threespine stickleback that evolved in the past 140 years within the invasive range in Switzerland. Using a controlled laboratory experiment with full-sib crosses and treatments mimicking a key feature of ecotypic niche divergence, we test if the phenotypic divergence that we observe in the wild results from phenotypic plasticity or divergent genetic predisposition. Our experimental groups show qualitatively similar phenotypic divergence as those observed among wild adults. The relative contribution of plasticity and divergent genetic predisposition differs among the traits studied, with traits related to the biomechanics of feeding showing a stronger genetic predisposition, whereas traits related to locomotion are mainly plastic. These results implicate that phenotypic plasticity and standing genetic variation interacted during contemporary ecotype formation in this case.

Variation in life-history traits and their plasticities to elevational transplantation among seed families suggests potential for adaptative evolution of 15 tropical plant species to climate change

• Premise of the study: Because not all plant species will be able to move in response to global warming, adaptive evolution matters largely for plant persistence. As prerequisites for adaptive evolution, genetic variation in and selection on phenotypic traits are needed, but these aspects have not been studied in tropical species. We studied how plants respond to transplantation to different elevations on Mt. Kilimanjaro, Tanzania, and whether there is quantitative genetic (among-seed family) variation in and selection on life-history traits and their phenotypic plasticity to the different environments. • Methods: We reciprocally transplanted seed families of 15 common tropical, herbaceous species of the montane and savanna vegetation zone at Mt. Kilimanjaro to a watered experimental garden in the montane (1450 m) and in the savanna (880 m) zone at the mountain’s slope and measured performance, reproductive, and phenological traits. • Results: Plants generally performed worse in the savanna garden, indicating that the savanna climate was more stressful and thus that plants may suffer from future climate warming. We found significant quantitative genetic variation in all measured performance and reproductive traits in both gardens and for several measures of phenotypic plasticity in response to elevational transplantation. Moreover, we found positive selection on traits at low and intermediate trait values levelling to neutral or negative selection at high values. • Conclusions: We conclude that common plants at Mt. Kilimanjaro express quantitative genetic variation in fitness-relevant traits and in their plasticities, suggesting potential to adapt evolutionarily to future climate warming and increased temperature variability.

Late Quaternary history of North Eurasian Norway spruce ( Picea abies ) and Siberian spruce ( Picea obovata ) inferred from macrofossils, pollen and cytoplasmic DNA variation

Aim We used combined palaeobotanical and genetic data to assess whether Norway spruce (Picea abies) and Siberian spruce (Picea obovata), two major components of the Eurasian boreal forests, occupied separate glacial refugia, and to test previous hypotheses on their distinction, geographical delimitation and introgression. Location The range of Norway spruce in northern Europe and Siberian spruce in northern Asia. Methods Pollen data and recently compiled macrofossil records were summarized for the Last Glacial Maximum (LGM), late glacial and Holocene. Genetic variation was assessed in 50 populations using one maternally (mitochondrial nad1) and one paternally (chloroplast trnT–trnL) inherited marker and analysed using spatial analyses of molecular variance (SAMOVA). Results Macrofossils showed that spruce was present in both northern Europe and Siberia at the LGM. Congruent macrofossil and pollen data from the late glacial suggested widespread expansions of spruce in the East European Plain, West Siberian Plain, southern Siberian mountains and the Baikal region. Colonization was largely completed during the early Holocene, except in the formerly glaciated area of northern Europe. Both DNA markers distinguished two highly differentiated groups that correspond to Norway spruce and Siberian spruce and coincide spatially with separate LGM spruce occurrences. The division of the mtDNA variation was geographically well defined and occurred to the east of the Ural Mountains along the Ob River, whereas the cpDNA variation showed widespread admixture. Genetic diversity of both DNA markers was higher in western than in eastern populations. Main conclusions North Eurasian Norway spruce and Siberian spruce are genetically distinct and occupied separate LGM refugia, Norway spruce on the East European Plain and Siberian spruce in southern Siberia, where they were already widespread during the late glacial. They came into contact in the basin of the Ob River and probably hybridized. The lower genetic diversity in the eastern populations may indicate that Siberian spruce suffered more from past climatic fluctuations than Norway spruce.

Prehistoric genomes reveal the genetic foundation and cost of horse domestication

The domestication of the horse revolutionized warfare, trade, and the exchange of people and ideas. This at least 5,500-y-long process, which ultimately transformed wild horses into the hundreds of breeds living today, is difficult to reconstruct from archeological data and modern genetics alone. We therefore sequenced two complete horse genomes, predating domestication by thousands of years, to characterize the genetic footprint of domestication. These ancient genomes reveal predomestic population structure and a significant fraction of genetic variation shared with the domestic breeds but absent from Przewalski’s horses. We find positive selection on genes involved in various aspects of locomotion, physiology, and cognition. Finally, we show that modern horse genomes contain an excess of deleterious mutations, likely representing the genetic cost of domestication.

Dihydropyrimidinase and β-ureidopropionase gene variation and severe fluoropyrimidine-related toxicity

AIMS To assess the association of DPYS and UPB1 genetic variation, encoding the catabolic enzymes downstream of dihydropyrimidine dehydrogenase, with early-onset toxicity from fluoropyrimidine-based chemotherapy. PATIENTS & METHODS The coding and exon-flanking regions of both genes were sequenced in a discovery subset (164 patients). Candidate variants were genotyped in the full cohort of 514 patients. RESULTS & CONCLUSIONS Novel rare deleterious variants in DPYS (c.253C > T and c.1217G > A) were detected once each in toxicity cases and may explain the occurrence of severe toxicity in individual patients, and associations of common variants in DPYS (c.1-1T > C: padjusted = 0.003; OR = 2.53; 95% CI: 1.39-4.62, and c.265-58T > C: padjusted = 0.039; OR = 0.61; 95% CI: 0.38-0.97) with 5-fluorouracil toxicity were replicated.

Implications of Longterm Diameter-Limit Harvesting: Effects on Radial Growth of Red Spruce (Picea rubens) and Genetic Diversity of White Pine (Pinus strobus)

For over 3 centuries, diameter-limit harvesting has been a predominant logging method in the northeastern United States. Silvicultural theory asserts that such intensively selective harvesting can lead to genetic degradation. A decrease in softwood productivity has recently been reported in Maine - has a long history of dysgenic selection degraded the genetic resources of Maine softwoods, contributing to a decrease in growth and productivity? This study examines two aspects of potential implications of diameter-limit harvesting: effects on residual phenotypes of red spruce and impacts on genetic diversity of white pine. Radial growth of residual red spruce trees in stands experiencing 50 years of fixed diameter-limit harvesting was measured using annual increment rings and compared with residual red spruce trees in positive selection stands. Trees remaiaing after several rounds of diameter-limit harvesting exhibited sigdicantl y smaller radial sizes throughout their lives, and displayed significantly slower growth rates for the first 80 years of measured growth. These results strongly suggest that the largest and fastest-growing genotypes and their respective gene complexes determining good radial growth have been removed from the diameter-limit stand. Dysgenic selection can be observed in fixed diarneter-limit stands, resulting in a diminished genetic resource and decreased residual stand value. To examine more direct genetic implications of long-term diameter-limit harvesting, microsatellite DNA markers were implemented to study genetic diversity of eastern white pine in Maine. Three age groups of trees were studied: mature trees older than 200 years, juvenile trees 5-30 years old, and embryos. Trees were genotyped at 10 microsatellite loci. Overall genetic diversity levels of eastern white pine in Maine were extremely high, with an average observed heterozygosity of 0.762. Genetic differentiation was minimal among and between all three age groups, although an excess of heterozygotes was shown in the mature and juvenile groups that was not reflected in the embryo group, which actually had a slight heterozygote deficiency. Allele frequencies did not differ significantly between age groups, but did reveal more rare and low frequency alleles in the embryo groups than in the mature group. Overall, low frequency alleles comprise the largest portion of alleles in the sample population, with no common alleles evident overall. These results suggest that significant genetic degradation has either not occurred for white pine, or that the results of dysgenic selection have not yet emerged. It is clear, however, that selective harvesting could result in a loss of low frequency alleles, which are a primary reserve of evolutionary potential in a species. Implications of these studies affect industrial forestry, regional economics, and ecological concerns for the northeast. Long-term diameter-limit harvesting can lead to a degradation of residual phenotypes, and an overall decrease in stand quality. Potentially, a loss of low frequency, locally adapted alleles could result in a decrease of allelic richness and degradation of the regidnal genetic resource. Decreased genetic variation can lead to seriously limited evolutionary potential of species and ecosystems, particularly in rapidly changing environments. Based on these findings, I recommend a reassessment of any harvesting prescription that includes fixed diameter-limit removals, particularly for species that have low natural genetic diversity levels or a limited natural range, such as red spruce. Maintenance of a healthy genetic reserve can avoid effects of dysgenic harvesting.

Adiposity and the perilipin gene: The role of single locus and multilocus variation, obesity-related quantitative traits, and disease-related phenotypes in the Atherosclerosis Risk in Communities (ARIC) study

Obesity is a complex multifactorial disease and is a public health priority. Perilipin coats the surface of lipid droplets in adipocytes and is believed to stabilize these lipid bodies by protecting triglyceride from early lipolysis. This research project evaluated the association between genetic variation within the human perilipin (PLIN) gene and obesity-related quantitative traits and disease-related phenotypes in Non-Hispanic White (NHW) and African American (AA) participants from the Atherosclerosis Risk in Communities (ARIC) Study. ^ Multivariate linear regression, multivariate logistic regression, and Cox proportional hazards models evaluated the association between single gene variants (rs2304794, rs894160, rs8179071, and rs2304795) and multilocus variation (rs894160 and rs2304795) within the PLIN gene and both obesity-related quantitative traits (body weight, body mass index [BMI], waist girth, waist-to-hip ratio [WHR], estimated percent body fat, and plasma total triglycerides) and disease-related phenotypes (prevalent obesity, metabolic syndrome [MetS], prevalent coronary heart disease [CHD], and incident CHD). Single variant analyses were stratified by race and gender within race while multilocus analyses were stratified by race. ^ Single variant analyses revealed that rs2304794 and rs894160 were significantly related to plasma triglyceride levels in all NHWs and NHW women. Among AA women, variant rs8179071 was associated with triglyceride levels and rs2304794 was associated with risk-raising waist circumference (>0.8 in women). The multilocus effects of variants rs894160 and rs2304795 were significantly associated with body weight, waist girth, WHR, estimated percent body fat, class II obesity (BMI ≥ 35 kg/m2), class III obesity (BMI ≥ 35 kg/m2), and risk-raising WHR (>0.9 in men and >0.8 in women) in AAs. Variant rs2304795 was significantly related to prevalent MetS among AA males and prevalent CHD in NHW women; multilocus effects of the PLIN gene were associated with prevalent CHD among NHWs. Rs2304794 was associated with incident CHD in the absence of the MetS among AAs. These findings support the hypothesis that variation within the PLIN gene influences obesity-related traits and disease-related phenotypes. ^ Understanding these effects of the PLIN genotype on the development of obesity can potentially lead to tailored health promotion interventions that are more effective. ^

Genetic variability and the hypothalamic control of energy balance

Background. Obesity is a major health problem throughout the industrialized world. Despite numerous attempts to curtail the rapid growth of obesity, its incidence continues to rise. Therefore, it is crucial to better understand the etiology of obesity beyond the concept of energy balance.^ Aims. The first aim of this study was to first investigate the relationship between eating behaviors and body size. The second goal was to identify genetic variation associated with eating behaviors. Thirdly, this study aimed to examine the joint relationships between eating behavior, body size and genetic variation.^ Methods. This study utilized baseline data ascertained in young adults from the Training Interventions and Genetics of Exercise (TIGER) Study. Variables assessed included eating behavior (Emotional Eating Scale, Eating Attitudes Test-26, and the Block98 Food Frequency Questionnaire), body size (body mass index, waist and hip circumference, waist/hip ratio, and percent body fat), genetic variation in genes implicated related to the hypothalamic control of energy balance, and appropriate covariates (age, gender, race/ethnicity, smoking status, and physical activity. For the genetic association analyses, genotypes were collapsed by minor allele frequency, and haplotypes were estimated for each gene. Additionally, Bayesian networks were constructed in order to determine the relationships between genetic variation, eating behavior and body size.^ Results. We report that the EAT-26 score, Caloric intake, percent fat, fiber intake, HEAT index, and daily servings of vegetables, meats, grains, and fats were significantly associated with at least one body size measure. Multiple SNPs in 17 genes and haplotypes from 12 genes were tested for their association with body size. Variation within both DRD4 and HTR2A was found to be associated with EAT-26 score. In addition, variation in the ghrelin gene (GHRL) was significantly associated with daily Caloric intake. A significant interaction between daily servings of grains and the HEAT index and variation within the leptin receptor gene (LEPR) was shown to influence body size.^ Conclusion. This study has shown that there is a substantial genetic component to eating behavior and that genetic variation interacts with eating behavior to influence body size.^

Molecular epidemiological studies on genetic predispositions to squamous cell carcinoma of the head and neck

Two molecular epidemiological studies were conducted to examine associations between genetic variation and risk of squamous cell carcinoma of the head and neck (SCCHN). In the first study, we hypothesized that genetic variation in p53 response elements (REs) may play roles in the etiology of SCCHN. We selected and genotyped five polymorphic p53 REs as well as a most frequently studied p53 codon 72 (Arg72Pro, rs1042522) polymorphism in 1,100 non-Hispanic White SCCHN patients and 1,122 age-and sex-matched cancer-free controls recruited at The University of Texas M. D. Anderson Cancer Center. In multivariate logistic regression analysis with adjustment for age, sex, smoking and drinking status, marital status and education level, we observed that the EOMES rs3806624 CC genotype had a significant effect of protection against SCCHN risk (adjusted odds ratio= 0.79, 95% confidence interval =0.64–0.98), compared with the -838TT+CT genotypes. Moreover, a significantly increased risk associated with the combined genotypes of p53 codon 72CC and EOMES -838TT+CT was observed, especially in the subgroup of non-oropharyneal cancer patients. The values of false-positive report probability were also calculated for significant findings. In the second study, we assessed the association between SCCHN risk and four potential regulatory single nucleotide polymorphisms (SNPs) of DEC1 (deleted in esophageal cancer 1) gene, a candidate tumor suppressor gene for esophageal cancer. After adjustment for age, sex, and smoking and drinking status, the variant -606CC (i.e., -249CC) homozygotes had a significantly reduced SCCHN risk (adjusted odds ratio = 0.71, 95% confidence interval = 0.52–0.99), compared with the -606TT homozygotes. Stratification analyses showed that a reduced risk associated with the -606CC genotype was more pronounced in subgroups of non-smokers, non-drinkers, younger subjects (defined as ≤ 57 years), carriers of TP53 Arg/Arg (rs1042522) genotype, patients with oropharyngeal cancer or late-stage SCCHN. Further in silico analysis revealed that the -249 T-to-C change led to a gain of a transcription factor binding site. Additional functional analysis showed that the -249T-to-C change significantly enhanced transcriptional activity of the DEC1 promoter and the DNA-protein binding activity. We conclude that the DEC1 promoter -249 T>C (rs2012775) polymorphism is functional, modulating susceptibility to SCCHN among non-Hispanic Whites. Additional large-scale, preferably population-based studies are needed to validate our findings.^

GENETIC PREDICTORS OF HYPERGLYCEMIA DUE TO HYDROCHLOROTHIAZIDE THERAPY

Response to pharmacological treatment is variable among individuals. Some patients respond favorably to a drug while others develop adverse reactions. Early investigations showed evidence of variation in genes that code for drug receptors, drug transporters, and drug metabolizing enzymes; and pharmacogenetics appeared as the science that studies the relationship between drug response and genetic variation. Thiazide diuretics are the recommended first-line monotherapy for hypertension (i.e. SBP>140 or DBP>90). Even so, diuretics are associated with adverse metabolic side effects, such as hyperglycemia, which increase the risk of developing type 2 diabetes. Published approaches testing variation in candidate genes (e.g. the renin-angiotensin-aldosteron system (RAAS) and salt–sensitivity genes) have met with only limited success. We conducted the first genome wide association study to identify genes influencing hyperglycemia as an adverse effect of thiazide diuretics in non-Hispanic White hypertensive patients participating in the Genetic Epidemiology of Responses to Antihypertensives (GERA) and Pharmacogenomic Evaluation of Antihypertensive Responses (PEAR) clinical trials. No SNP reached the a priori defined threshold of statistical significance (p<5x10-8). We detected 50 SNPs in 9 genomic regions with suggestive p-values (p<1x10-5). Two of them, rs6870564 (p-value=3.28 X 10-6) and rs7702121 (p-value=5.09 X 10-6), were located close to biologic candidate genes, MYO and MGAT1, and one SNP in a genomic region in chromosome 6, rs7762018 (p-value=4.59 X 10-6) has been previously related to Insulin-Dependent Diabetes Mellitus (IDDM8). I conclude that 1) there are unlikely to be common SNPs with large effects on the adverse metabolic effects to hydrochlorothiazide treatment and 2) larger sample sizes are needed for pharmacogenetic studies of inter-individual variation in response to commonly prescribed medication.

Identification of genetic risk factors for cerebellar mutism in pediatric brain tumor patients

Following posterior fossa surgery for resection of childhood medulloblastoma and primitive neuroectodermal tumor (M/PNET), cerebellar mutism (CM) may develop. This is a condition of absent or diminished speech in a conscious patient with no evidence of oral apraxia, which can be accompanied by other symptoms of the posterior fossa syndrome complex, which includes ataxia and hypotonia. Little is known about the etiology. Therefore, we conducted a SNP, gene, and pathway-level analysis to assess the role of host genetic variation on the risk of CM in M/PNET subjects following treatment. Cases (n= 20) and controls (n= 53) were recruited from the Childhood Cancer Epidemiology and Prevention Center, in Houston, TX. DNA samples were genotyped using the Illumina Human 1M Quad SNP chip. Ten pathways were identified from logistic regression used to identify the marginal effect of each SNP on CM risk. The minP test was conducted to identify associations between SNPs categorized to genes and CM risk. Pathways were assessed to determine if there was a significant enrichment of genes in the pathway compared to all other pathways. There were 78 genes that reached the threshold of min P ≤0.05 in 948 genes. The Neurotoxicity pathway was the most significant pathway after adjusting for multiple comparisons (q=0.040 and q=0.005, using Fisher's exact test and a test of proportions, respectively). Most genes within the Neurotoxicity pathway that reached a threshold of minP ≤0.05 were known to have an apoptosis function, possibly inducing neuronal apoptosis in the dentatothalamocortical pathway, and may be important in CM etiology in this population. This is the first study to assess the potential role of genetic risk factors on CM. As an exploratory study, these results should be replicated in a larger sample. ^

Competition may explain the fine-scale spatial patterns and genetic structure of two co-occurring plant congeners

1. The spatial distribution of individual plants within a population and the population’s genetic structure are determined by several factors, like dispersal, reproduction mode or biotic interactions. The role of interspecific interactions in shaping the spatial genetic structure of plant populations remains largely unknown. 2. Species with a common evolutionary history are known to interact more closely with each other than unrelated species due to the greater number of traits they share. We hypothesize that plant interactions may shape the fine genetic structure of closely related congeners. 3. We used spatial statistics (georeferenced design) and molecular techniques (ISSR markers) to understand how two closely related congeners, Thymus vulgaris (widespread species) and T. loscosii (narrow endemic) interact at the local scale. Specific cover, number of individuals of both study species and several community attributes were measured in a 10 × 10 m plot. 4. Both species showed similar levels of genetic variation, but differed in their spatial genetic structure. Thymus vulgaris showed spatial aggregation but no spatial genetic structure, while T. loscosii showed spatial genetic structure (positive genetic autocorrelation) at short distances. The spatial pattern of T. vulgaris’ cover showed significant dissociation with that of T. loscosii. The same was true between the spatial patterns of the cover of T. vulgaris and the abundance of T. loscosii and between the abundance of each species. Most importantly, we found a correlation between the genetic structure of T. loscosii and the abundance of T. vulgaris: T. loscosii plants were genetically more similar when they were surrounded by a similar number of T. vulgaris plants. 5. Synthesis. Our results reveal spatially complex genetic structures of both congeners at small spatial scales. The negative association among the spatial patterns of the two species and the genetic structure found for T. loscosii in relation to the abundance of T. vulgaris indicate that competition between the two species may account for the presence of adapted ecotypes of T. loscosii to the abundance of a competing congeneric species. This suggests that the presence and abundance of close congeners can influence the genetic spatial structure of plant species at fine scales.