GWAS of butyrylcholinesterase activity identifies four novel loci, independent effects within BCHE and secondary associations with metabolic risk factors

Serum butyrylcholinesterase (BCHE) activity is associated with obesity, blood pressure and biomarkers of cardiovascular and diabetes risk. We have conducted a genome-wide association scan to discover genetic variants affecting BCHE activity, and to clarify whether the associations between BCHE activity and cardiometabolic risk factors are caused by variation in BCHE or whether BCHE variation is secondary to the metabolic abnormalities. We measured serum BCHE in adolescents and adults from three cohorts of Australian twin and family studies. The genotypes from approximately 2.4 million single-nucleotide polymorphisms (SNPs) were available in 8791 participants with BCHE measurements. We detected significant associations with BCHE activity at three independent groups of SNPs at the BCHE locus (P = 5.8 x 10(-262), 7.8 x 10(-47), 2.9 x 10(-12)) and at four other loci: RNPEP (P = 9.4 x 10(-16)), RAPH1-ABI2 (P = 4.1 x 10(-18)), UGT1A1 (P = 4.0 x 10(-8)) and an intergenic region on chromosome 8 (P = 1.4 x 10(-8)). These loci affecting BCHE activity were not associated with metabolic risk factors. On the other hand, SNPs in genes previously associated with metabolic risk had effects on BCHE activity more often than can be explained by chance. In particular, SNPs within FTO and GCKR were associated with BCHE activity, but their effects were partly mediated by body mass index and triglycerides, respectively. We conclude that variation in BCHE activity is due to multiple variants across the spectrum from uncommon/large effect to common/small effect, and partly results from (rather than causes) metabolic abnormalities.

Adult-type hypolactasia: Genotype-phenotype correlation

Adult-type hypolactasia (primary lactose malabsorption, lactase non-persistence) is the most common enzyme deficiency worldwide, and manifests with symptoms of lactose intolerance such as abdominal pain, gas formation and diarrhea. In humans with adult-type hypolactasia, lactase activity is high at birth, but declines during childhood to about one-tenth of the activity at birth. In 2002, a one base polymorphism C/T-13910, located 14 kilobases from the starting codon of the lactase-phlorizin hydrolase (LPH) gene was observed to be associated with the persistence of lactase activity. The T-13910 allele (C/T-13910 and T/T-13910 genotypes) associates with persistence of lactase activity throughout life, whereas the C/C-13910 genotype associates with adult-type hypolactasia. In this thesis work, the timing and mechanism of decline of lactase enzyme activity during development was studied using the C/T-13910 polymorphism as a molecular marker. We observed an excellent correlation between low lactase activity and the C/C-13910 genotype in all subjects > 12 years of age, irrespective their ethnicity. In children of African origin, the lactase activity declined somewhat earlier than among Finnish children. Furthermore, we observed an increasing imbalance in the relative lactase mRNA expression from the C-13910 and T-13910 alleles in Finnish children beginning from five years of age. The genetic test for adult-type hypolactasia showed a sensitivity of 93% and a specificity of 100% in the Finnish children and adolescents > 12 years of age. The relation of milk consumption and the milk-related abdominal complaints to the C/T-13910 genotypes associated with lactase persistence/non-persistence was studied by a questionnaire-based approach in > 2100 Finns. Both Finnish children and adults with the C/C-13910 genotype consumed significantly less dairy products compared to those with the C/T-13910 and T/T-13910 genotypes. Flatulence was the only of the abdominal symptoms of lactose intolerance that subjects with the C/C-13910 genotype reported significantly more often than those with the C/T-13910 and T/T-13910 genotypes. A minor proportion (<10%) of subjects with the C/C-13910 genotype, nevertheless, reported drinking milk without any symptoms afterwards. There was no association between cow's milk allergy starting as a newborn and adult-type hypolactasia. In an association study an increased risk of colorectal cancer was observed among those with molecular diagnosis of adult-type hypolactasia. It warrants further studies to clarify whether the increased risk observed in the Finnish population is associated with lactose or decreased intake of dairy products in these subjects.

A simple and fast two-locus quality control test to detect false positives due to batch effects in genome-wide association studies

The impact of erroneous genotypes having passed standard quality control (QC) can be severe in genome-wide association studies, genotype imputation, and estimation of heritability and prediction of genetic risk based on single nucleotide polymorphisms (SNP). To detect such genotyping errors, a simple two-locus QC method, based on the difference in test statistic of association between single SNPs and pairs of SNPs, was developed and applied. The proposed approach could detect many problematic SNPs with statistical significance even when standard single SNP QC analyses fail to detect them in real data. Depending on the data set used, the number of erroneous SNPs that were not filtered out by standard single SNP QC but detected by the proposed approach varied from a few hundred to thousands. Using simulated data, it was shown that the proposed method was powerful and performed better than other tested existing methods. The power of the proposed approach to detect erroneous genotypes was approximately 80% for a 3% error rate per SNP. This novel QC approach is easy to implement and computationally efficient, and can lead to a better quality of genotypes for subsequent genotype-phenotype investigations.

Replication of the association of common rs9939609 variant of FTO with increased BMI in an Australian adult twin population but no evidence for gene by environment (G x E) interaction

OBJECTIVE: To further investigate a common variant (rs9939609) in the fat mass- and obesity-associated gene (FTO), which recent genome-wide association studies have shown to be associated with body mass index (BMI) and obesity. DESIGN: We examined the effect of this FTO variant on BMI in 3353 Australian adult male and female twins. RESULTS: The minor A allele of rs9939609 was associated with an increased BMI (P=0.0007). Each additional copy of the A allele was associated with a mean BMI increase of approximately 1.04 kg/m(2) (approximately 3.71 kg). Using variance components decomposition, we estimate that this single-nucleotide polymorphism accounts for approximately 3% of the genetic variance in BMI in our sample (approximately 2% of the total variance). By comparing intrapair variances of monozygotic twins of different genotypes we were able to perform a direct test of gene by environment (G x E) interaction in both sexes and gene by parity (G x P) interaction in women, but no evidence was found for either. CONCLUSIONS: In addition to supporting earlier findings that the rs9939609 variant in the FTO gene is associated with an increased BMI, our results indicate that the associated genetic effect does not interact with environment or parity.

Rain Damage to Strawberries Grown in Southeast Queensland: Evaluation and Genetic Control.

In Queensland, Australia, strawberries (Fragaria xananassa Duchesne) are grown in open fields and rainfall events can damage fruit. Cultivars that are resistant to rain damage may reduce losses and lower risk for the growers. However, little is known about the genetic control of resistance and in a subtropical climate, unpredictable rainfall events hamper evaluation. Rain damage was evaluated on seedling and clonal trials of one breeding population comprising 645 seedling genotypes and 94 clones and on a second clonal population comprising 46 clones from an earlier crossing to make preliminary estimates of heritability. The incidence of field damage from rainfall and damage after laboratory soaking was evaluated to determine if this soaking method could be used to evaluate resistance to rain damage. Narrow-sense heritability of resistance to rain damage calculated for seedlings was low (0.21 +/- 0.15) and not significantly different from zero; however, broad-sense heritability estimates were moderate in both seedlings (0.49 +/- 0.16) and clones (0.45 +/- 0.08) from the first population and similar in clones (0.56 +/- 0.21) from the second population. Immersion of fruit in deionized water produced symptoms consistent with rain damage in the field. Lengthening the duration of soaking of 'Festival' fruit in deionized water exponentially increased the proportion of damage to fruit ranging in ripeness from immature to ripe during the first 6-h period of soaking. When eight genotypes were evaluated, the proportion of sound fruit after soaking in deionized water in the laboratory for up to 5 h was linearly related (r(2) = 0.90) to the proportion of sound fruit in the field after 89 mm of rain. The proportion of sound fruit of the breeding genotype '2008-208' and 'Festival' under soaking (0.67, 0.60) and field (0.52, 0.43) evaluations, respectively, is about the same and these genotypes may be useful sources of resistance to rain damage.

Hybridisation of Australian chickpea cultivars with wild Cicer spp. increases resistance to root-lesion nematodes (Pratylenchus thornei and P. neglectus).

Australian and international chickpea (Cicer arietinum) cultivars and germplasm accessions, and wild annual Cicer spp. in the primary and secondary gene pools, were assessed in glasshouse experiments for levels of resistance to the root-lesion nematodes Pratylenchus thornei and P. neglectus. Lines were grown in replicated experiments in pasteurised soil inoculated with a pure culture of either P. thornei or P. neglectus and the population density of the nematodes in the soil plus roots after 16 weeks growth was used as a measure of resistance. Combined statistical analyses of experiments (nine for P. thornei and four for P. neglectus) were conducted and genotypes were assessed using best linear unbiased predictions. Australian and international chickpea cultivars possessed a similar range of susceptibilities through to partial resistance. Wild relatives from both the primary (C. reticulatum and C. echinospermum) and secondary (C. bijugum) gene pools of chickpea were generally more resistant than commercial chickpea cultivars to either P. thornei or P. neglectus or both. Wild relatives of chickpea have probably evolved to have resistance to endemic root-lesion nematodes whereas modern chickpea cultivars constitute a narrower gene pool with respect to nematode resistance. Resistant accessions of C. reticulatum and C. echinospermum were crossed and topcrossed with desi chickpea cultivars and resistant F(4) lines were obtained. Development of commercial cultivars with the high levels of resistance to P. thornei and P. neglectus in these hybrids will be most valuable for areas of the Australian grain region and other parts of the world where alternating chickpea and wheat crops are the preferred rotation.

Human Torque teno virus: epidemiology, cell biology and immunology

Torque teno virus (TTV) was discovered in 1997 in the serum of a Japanese patient who had a post-transfusion hepatitis of unknown etiology. It is a small virus containing a circular single-stranded DNA genome which is unique among human viruses. Within a few years after its discovery, the TTVs were noted to form a large family of viruses with numerous genotypes. TTV is highly prevalent among the general population throughout the world, and persistent infections and co-infections with several genotypes occur frequently. However, the pathogenicity and the mechanism for the sustained occurrence of the virus in blood are at present unclear. To determine the prevalence of TTV in Finland, we set up PCR methods and examined the sera of asymptomatic subjects for the presence of TTV DNA and for genotype-6 DNA. TTV was found to be highly prevalent also in Finland; 85% of adults harbored TTV in their blood, and 4% were infected with genotype-6. In addition, TTV DNA was detected in a number of different tissues, with no tissue-type or symptom specificity. Most cell-biological events during TTV infections are at the moment unknown. Replicating TTV DNA has, however, been detected in liver and the hematopoietic compartment, and three mRNAs are known to be generated. To characterize TTV cell biology in more detail, we cloned in full length the genome of TTV genotype 6. We showed that in human kidney-derived cells TTV produces altogether six proteins with distinct subcellular localizations. TTV mRNA transcription was detected in all cell lines transfected with the full-length clone, and TTV DNA replicated in several of them, including those of erythroid, kidney, and hepatic origin. Furthermore, the viral DNA replication was shown to utilize the cellular DNA polymerases. Diagnoses of TTV infections have been based almost solely on PCR, whereas serological tests, measuring antibody responses, would give more information on many aspects of these infections. To investigate the TTV immunology in more detail, we produced all six TTV proteins for use as antigens in serological tests. We detected in human sera IgM and IgG antibodies to occur simultaneously with TTV DNA, and observed appearance of TTV DNA regardless of pre-existing antibodies, and disappearance of TTV DNA after antibody appearance. The genotype-6 nucleotide sequence remained stable for years within the infected subjects, suggesting that some mechanism other than mutations is used by this minute virus to evade our immune system and to establish chronic infections in immunocompetent subjects.

Genetics of Cardiovascular Disease: A Candidate Gene Study of USF1

Cardiovascular diseases (CVD) are major contributors to morbidity and mortality worldwide. Several interacting environmental, biochemical, and genetic risk factors can increase disease susceptibility. While some of the genes involved in the etiology of CVD are known, many are yet to be discovered. During the last few decades, scientists have searched for these genes with genome-wide linkage and association methods, and with more targeted candidate gene studies. This thesis investigates variation within the upstream transcription factor 1 (USF1) gene locus in relation to CVD risk factors, atherosclerosis, and incidence and prevalence of CVD. This candidate gene was first identified in Finnish families ascertained for familial combined hyperlipidemia, a common dyslipidemia predisposing to coronary heart disease. The gene is a ubiquitously expressed transcription factor regulating expression of several genes from lipid and glucose metabolism, inflammation, and endothelial function. First, we examined association between USF1 variants and several CVD risk factors, such as lipid phenotypes, body composition measures, and metabolic syndrome, in two prospective population cohorts. Our data suggested that USF1 contributes to these CVD risk factors at the population level. Notably, the associations with quantitative measurements were mostly detected among study subjects with CVD or metabolic syndrome, suggesting complex interactions between USF1 effects and the pathophysiological state of an individual. Second, we investigated how variation at the USF1 locus contributes to atherosclerotic lesions of the coronary arteries and abdominal aorta. For this, we used two study samples of middle-aged men with detailed measurements of atherosclerosis obtained in autopsy. USF1 variation significantly associated with areas of several types of lesions, especially with calcification of the arteries. Next, we tested what effect the USF1 risk variants have on sudden cardiac death and incidence of CVD. The atherosclerosis-associated risk variant increased the risk of sudden cardiac death of the same study subjects. Furthermore, USF1 alleles associated with incidence of CVD in the Finnish population follow-up cohorts. These associations were especially prominent among women, suggesting a sex specific effect, which has also been detected in subsequent studies. Finally, as some of the low-yield DNA samples of the Finnish follow-up study cohort needed to be whole-genome amplified (WGA) prior to genotyping, we evaluated whether the produced WGA genotypes were of good quality. Although the samples giving genotype discrepancies could not be detected before genotyping with standard laboratory quality control methods, our results suggested that enhanced quality control at the time of the genotyping could identify such samples. In addition, combining two WGA reactions into one pooled DNA sample for genotyping markedly reduced the number of discrepancies and samples showing them. In conclusion, USF1 seems to have a role in the etiology of CVD. Additional studies are warranted to identify functional variants and to study interactions between USF1 and other genetic or environmental factors. This USF1 study, and other studies with low DNA yield of some samples, can benefit from whole genome amplification of the low-yield samples prior to genotyping. Careful quality control procedures are, however, needed in WGA genotyping.

Influence of plant traits on production costs and profitability of strawberry in southeast Queensland.

In Queensland the subtropical strawberry ( Fragaria * ananassa) breeding program aims to combine traits into novel genotypes that increase production efficiency. The contribution of individual plant traits to cost and income under subtropical Queensland conditions was investigated, with the overall goal of improving the profitability of the industry through the release of new strawberry cultivars. The study involved specifying the production and marketing system using three cultivars of strawberry that are currently widely grown annually in southeast Queensland, developing methods to assess the economic impact of changes to the system, and identifying plant traits that influence outcomes from the system. From May through September P (price; $ punnet -1), V (monthly mass; tonne of fruit on the market) and M (calendar month; i.e. May=5) were found to be related ( r2=0.92) by the function (SE) P=4.741(0.469)-0.001630(0.0005) V-0.226(0.102) M using data from 2006 to 2010 for the Brisbane central market. Both income and cost elements in the gross margin were subject to sensitivity analysis. 'Harvesting' and 'Handling/Packing' 'Groups' of 'Activities' were the major contributors to variable costs (each >20%) in the gross margin analysis. Within the 'Harvesting Group', the 'Picking Activity' contributed most (>80%) with the trait 'display of fruit' having the greatest (33%) influence on the cost of the 'Picking Activity'. Within the 'Handling/Packing Group', the 'Packing Activity' contributed 50% of costs with the traits 'fruit shape', 'fruit size variation' and 'resistance to bruising' having the greatest (12-62%) influence on the cost of the 'Packing Activity'. Non-plant items (e.g. carton purchases) made up the other 50% of the costs within the 'Handling/Packing Group'. When any of the individual traits in the 'Harvesting' and 'Handling/Packing' groups were changed by one unit (on a 1-9 scale) the gross margin changed by up to 1%. Increasing yield increased the gross margin to a maximum (15% above present) at 1320 g plant -1 (94% above present). A 10% redistribution of total yield from September to May increased the gross margin by 23%. Increasing fruit size increased gross margin: a 75% increase in fruit size (to ~30 g) produced a 22% increase in the gross margin. The modified gross margin analysis developed in this study allowed simultaneous estimation of the gross margin for the producer and gross value of the industry. These parameters sometimes move in opposite directions.

Multidrug-resistant extraintestinal pathogenic Escherichia coli of sequence type ST131 in animals and foods.

Multidrug-resistant Escherichia colt sequence type 131 (51131) has recently emerged as a globally distributed cause of extraintestinal infections in humans. Diverse factors have been investigated as explanations for ST131's rapid and successful dissemination, including transmission through animal contact and consumption of food, as suggested by the detection of ST131 in a number of nonhuman species. For example, ST131 has recently been identified as a cause of clinical infection in companion animals and poultry, and both host groups have been confirmed as faecal carriers of ST131. Moreover, a high degree of similarity has been shown among certain ST131 isolates from humans, companion animals, and poultry based on resistance characteristics and genomic background and human and companion animal ST131 isolates tend to exhibit similar virulence genotypes. However, most ST131 isolates from poultry appear to possess specific virulence genes that are typically absent from human and companion animal isolates, including genes associated with avian pathogenic E. coli. Since the number of reported animal and food-associated ST131 isolates is quite small, the role of nonhuman host species in the emergence, dissemination, and transmission of ST131 to humans remains unclear. Nevertheless, given the profound public health importance of the emergent ST131 clonal group, even the limited available evidence indicates a pressing need for further careful study of this significant question.

Rapid microsatellite marker development for African mahogany ( Khaya senegalensis, Meliaceae) using next-generation sequencing and assessment of its intra-specific genetic diversity.

Khaya senegalensis (African mahogany or dry-zone mahogany) is a high-value hardwood timber species with great potential for forest plantations in northern Australia. The species is distributed across the sub-Saharan belt from Senegal to Sudan and Uganda. Because of heavy exploitation and constraints on natural regeneration and sustainable planting, it is now classified as a vulnerable species. Here, we describe the development of microsatellite markers for K. senegalensis using next-generation sequencing to assess its intra-specific diversity across its natural range, which is a key for successful breeding programs and effective conservation management of the species. Next-generation sequencing yielded 93943 sequences with an average read length of 234bp. The assembled sequences contained 1030 simple sequence repeats, with primers designed for 522 microsatellite loci. Twenty-one microsatellite loci were tested with 11 showing reliable amplification and polymorphism in K. senegalensis. The 11 novel microsatellites, together with one previously published, were used to assess 73 accessions belonging to the Australian K. senegalensis domestication program, sampled from across the natural range of the species. STRUCTURE analysis shows two major clusters, one comprising mainly accessions from west Africa (Senegal to Benin) and the second based in the far eastern limits of the range in Sudan and Uganda. Higher levels of genetic diversity were found in material from western Africa. This suggests that new seed collections from this region may yield more diverse genotypes than those originating from Sudan and Uganda in eastern Africa.

Water and thermal regimes for field pea in Australia and their implications for breeding.

There is a large gap between the refined approaches to characterise genotypes and the common use of location and season as a coarse surrogate for environmental characterisation of breeding trials. As a framework for breeding, the aim of this paper is quantifying the spatial and temporal patterns of thermal and water stress for field pea in Australia. We compiled a dataset for yield of the cv. Kaspa measured in 185 environments, and investigated the associations between yield and seasonal patterns of actual temperature and modelled water stress. Correlations between yield and temperature indicated two distinct stages. In the first stage, during crop establishment and canopy expansion before flowering, yield was positively associated with minimum temperature. Mean minimum temperature below similar to 7 degrees C suggests that crops were under suboptimal temperature for both canopy expansion and radiation-use efficiency during a significant part of this early growth period. In the second stage, during critical reproductive phases, grain yield was negatively associated with maximum temperature over 25 degrees C. Correlations between yield and modelled water supply/demand ratio showed a consistent pattern with three phases: no correlation at early stages of the growth cycle, a progressive increase in the association that peaked as the crop approached the flowering window, and a progressive decline at later reproductive stages. Using long-term weather records (1957-2010) and modelled water stress for 104 locations, we identified three major patterns of water deficit nation wide. Environment type 1 (ET1) represents the most favourable condition, with no stress during most of the pre-flowering phase and gradual development of mild stress after flowering. Type 2 is characterised by increasing water deficit between 400 degree-days before flowering and 200 degree-days after flowering and rainfall that relieves stress late in the season. Type 3 represents the more stressful condition with increasing water deficit between 400 degree-days before flowering and maturity. Across Australia, the frequency of occurrence was 24% for ET1, 32% for ET2 and 43% for ET3, highlighting the dominance of the most stressful condition. Actual yield averaged 2.2 t/ha for ET1, 1.9 t/ha for ET2 and 1.4 t/ha for ET3, and the frequency of each pattern varied substantially among locations. Shifting from a nominal (i.e. location and season) to a quantitative (i.e. stress type) characterisation of environments could help improving breeding efficiency of field pea in Australia.

Genome-wide association studies of quantitative traits with related individuals: Little (power) lost but much to be gained

For complex disease genetics research in human populations, remarkable progress has been made in recent times with the publication of a number of genome-wide association scans (GWAS) and subsequent statistical replications. These studies have identified new genes and pathways implicated in disease, many of which were not known before. Given these early successes, more GWAS are being conducted and planned, both for disease and quantitative phenotypes. Many researchers and clinicians have DNA samples available on collections of families, including both cases and controls. Twin registries around the world have facilitated the collection of large numbers of families, with DNA and multiple quantitative phenotypes collected on twin pairs and their relatives. In the design of a new GWAS with a fixed budget for the number of chips, the question arises whether to include or exclude related individuals. It is commonly believed to be preferable to use unrelated individuals in the first stage of a GWAS because relatives are 'over-matched' for genotypes. In this study, we quantify that for GWAS of a quantitative phenotype, relative to a sample of unrelated individuals surprisingly little power is lost when using relatives. The advantages of using relatives are manifold, including the ability to perform more quality control, the choice to perform within-family tests of association that are robust to population stratification, and the ability to perform joint linkage and association analysis. Therefore, the advantages of using relatives in GWAS for quantitative traits may well outweigh the small disadvantage in terms of statistical power.

Raw data for "Linkage disequilibrium estimation of effective population size with immigrants from divergent populations: a case study on Spanish mackerel (Scomberomorus commerson)."

Abstract of Macbeth, G. M., Broderick, D., Buckworth, R. & Ovenden, J. R. (In press, Feb 2013). Linkage disequilibrium estimation of effective population size with immigrants from divergent populations: a case study on Spanish mackerel (Scomberomorus commerson). G3: Genes, Genomes and Genetics. Estimates of genetic effective population size (Ne) using molecular markers are a potentially useful tool for the management of endangered through to commercial species. But, pitfalls are predicted when the effective size is large, as estimates require large numbers of samples from wild populations for statistical validity. Our simulations showed that linkage disequilibrium estimates of Ne up to 10,000 with finite confidence limits can be achieved with sample sizes around 5000. This was deduced from empirical allele frequencies of seven polymorphic microsatellite loci in a commercially harvested fisheries species, the narrow barred Spanish mackerel (Scomberomorus commerson). As expected, the smallest standard deviation of Ne estimates occurred when low frequency alleles were excluded. Additional simulations indicated that the linkage disequilibrium method was sensitive to small numbers of genotypes from cryptic species or conspecific immigrants. A correspondence analysis algorithm was developed to detect and remove outlier genotypes that could possibly be inadvertently sampled from cryptic species or non-breeding immigrants from genetically separate populations. Simulations demonstrated the value of this approach in Spanish mackerel data. When putative immigrants were removed from the empirical data, 95% of the Ne estimates from jacknife resampling were above 24,000.

Variants in EMX2 and PTEN do not contribute to risk of endometriosis

Endometriosis has a genetic component, and significant linkage has been found to a region on chromosome 10q. Two candidate genes, EMX2 and PTEN, implicated in both endometriosis and endometrial cancer, lie on chromosome 10q. We hypothesized that variation in EMX2 and/or PTEN could contribute to the risk of endometriosis and may account for some of the linkage signal on 10q. We genotyped single nucleotide polymorphisms (SNPs) in a case-control design to evaluate association between endometriosis and common variations in these two genes. The genotyping and statistical analysis were based on samples collected from Australian volunteers. The cases were 768 unrelated women with surgically confirmed endometriosis selected from affected sister pair (ASP) families participating in the Australian Genes behind Endometriosis Study. The controls were 768 female participants in twin studies who, based on screening questions, did not have a diagnosis of endometriosis. Genotypes of 22 SNPs in the EMX2 gene and 15 SNPs in the PTEN gene were the main outcome measures. Statistical analysis provided measures of linkage disequilibrium and association. Permutation testing showed no globally significant association between any SNPs or haplotypes and endometriosis for either gene. It is unlikely that the EMX2 or PTEN gene variants investigated contribute to risk for initiation and/or development of endometriosis.