Association between ORMDL3, IL1RL1 and a deletion on chromosome 17q21 with asthma risk in Australia

Genome-wide association studies followed by replication provide a powerful approach to map genetic risk factors for asthma. We sought to search for new variants associated with asthma and attempt to replicate the association with four loci reported previously (ORMDL3, PDE4D, DENND1B and IL1RL1). Genome-wide association analyses of individual single nucleotide polymorphisms (SNPs), rare copy number variants (CNVs) and overall CNV burden were carried out in 986 asthma cases and 1846 asthma-free controls from Australia. The most-associated locus in the SNP analysis was ORMDL3 (rs6503525, P = 4.8 x 10(-)(7)). Five other loci were associated with P < 10(-)(5), most notably the chemokine CXC motif ligand 14 (CXCL14) gene (rs31263, P = 7.8 x 10(-)(6)). We found no evidence for association with the specific risk variants reported recently for PDE4D, DENND1B and ILR1L1. However, a variant in IL1RL1 that is in low linkage disequilibrium with that reported previously was associated with asthma risk after accounting for all variants tested (rs10197862, gene wide P = 0.01). This association replicated convincingly in an independent cohort (P = 2.4 x 10(-)(4)). A 300-kb deletion on chromosome 17q21 was associated with asthma risk, but this did not reach experiment-wide significance. Asthma cases and controls had comparable CNV rates, length and number of genes affected by deletions or duplications. In conclusion, we confirm the association between asthma risk and variants in ORMDL3 and identify a novel risk variant in IL1RL1. Follow-up of the 17q21 deletion in larger cohorts is warranted.

Novel Multiple Sclerosis Predisposing Genetic Variants Outside the HLA Region

Multiple sclerosis (MS) is a chronic inflammatory disease of the central nervous system (CNS). Both environmental factors and several predisposing genes are required to generate MS. Despite intensive research these risk factors are still largely unknown, the pathogenesis of MS demyelination is poorly understood, and no curative treatment exists. Both prevalence and familial occurrence of MS are exceptionally high in a Finnish population subisolate, Southern Ostrobothnia, presumably due to enrichment of predisposing genetic variants within this region. Previous linkage scan on MS pedigrees from Southern Ostrobothnia detected three main MS loci on chromosomes 5p, 6p (HLA) and 17q. Linkage studies in other populations have also provided independent evidence for the location of MS susceptibility genes in these regions. Further, these loci are syntenic to the experimental autoimmune encephalomyelitis (EAE) susceptibility loci of rodents. In this thesis work an effort was made to localize MS predisposing alleles of the linked loci outside the HLA region by studying familial MS cases from the Southern Ostrobothnia isolate. Analysis of the 5p locus revealed one region, flanking the complement component 7 (C7) gene. The identified relatively rare haplotype seems to have a fairly large effect on genetic susceptibility of MS (frequency MS 12%, controls 4%; p=0.000003, OR=2.73). Evidence for association with alleles of the region and MS was seen also in more heterogeneous populations. Convincingly, plasma C7 protein levels and complement activity correlated with the risk haplotype identified. The finding stimulated us to study other complement cascade genes in MS. No evidence for association could be observed with the complement component coding genes outside 5p. A scan of the 17q locus provided evidence for association with variants of the protein kinase C alpha (PRKCA) gene (p=0.0001). Modest evidence for association with PRKCA was observed also in Canadian MS families. Finally we used a candidate gene based approach to identify potential MS loci. Mutations of DAP12 and TREM2 cause a recessively inherited CNS white matter disease PLOSL. Interestingly, DAP12 and TREM2 are located in MS regions on 6p and 19q, and we tested them as potential candidate genes in the Finnish MS sample. No evidence for association with MS was observed. This thesis provides an example of how extended families from special populations can be utilized in fine-mapping of the linked loci. A first relatively rare MS variant was identified utilizing the strength of a Finnish population subisolate. This variant seems to have an effect on activity of the complement system, which has previously been suggested to have an important role in the pathogenesis of MS.

Insights into the molecular genetics of celiac disease : Applying family- and population-based strategies

Celiac disease, or gluten intolerance, is triggered by dietary glutens in genetically susceptible individuals and it affects approximately 1% of the Caucasian population. The best known genetic risk factors for celiac disease are HLA DQ2 and DQ8 heterodimers, which are necessary for the development of the disease. However, they alone are not sufficient for disease induction, other risk factors are required. This thesis investigated genetic factors for celiac disease, concentrating on susceptibility loci on chromosomes 5q31-q33, 19p13 and 2q12 previously reported in genome-wide linkage and association studies. In addition, a novel genotyping method for the detection of HLA DQ2 and DQ8 coding haplotypes was validated. This study was conducted using Finnish and Hungarian family materials, and Finnish, Hungarian and Italian case-control materials. Genetic linkage and association were analysed in these materials using candidate gene and fine-mapping approaches. The results confirmed linkage to celiac disease on the chromosomal regions 5q31-q33 and 19p13. Fine-mapping on chromosome 5q31-q33 revealed several modest associations in the region, and highlighted the need for further investigations to locate the causal risk variants. The MYO9B gene on chromosome 19p13 showed evidence for linkage and association particularly with dermatitis herpetiformis, the skin manifestation of celiac disease. This implies a potential difference in the genetic background of the intestinal and skin forms of the disease, although studies on larger samplesets are required. The IL18RAP locus on chromosome 2q12, shown to be associated with celiac disease in a previous genome-wide association study and a subsequent follow-up, showed association in the Hungarian population in this study. The expression of IL18RAP was further investigated in small intestinal tissue and in peripheral blood mononuclear cells. The results showed that IL18RAP is expressed in the relevant tissues. Two putative isoforms of IL18RAP were detected by Western blot analysis, and the results suggested that the ratios and total levels of these isoforms may contribute to the aetiology of celiac disease. A novel genotyping method for celiac disease-associated HLA haplotypes was also validated in this thesis. The method utilises single-nucleotide polymorphisms tagging these HLA haplotypes with high sensitivity and specificity. Our results suggest that this method is transferable between populations, and it is suitable for large-scale analysis. In conclusion, this doctorate study provides an insight into the roles of the 5q31-q33, MYO9B, IL18RAP and HLA loci in the susceptibility to celiac disease in the Finnish, Hungarian and Italian populations, highlighting the need for further studies at these genetic loci and examination of the function of the candidate genes.

A genome wide linkage scan for dizygotic twinning in 525 families of mothers of dizygotic twins

BACKGROUND: The tendency to conceive dizygotic (DZ) twins is a complex trait influenced by genetic and environmental factors. To search for new candidate loci for twinning, we conducted a genome-wide linkage scan in 525 families using microsatellite and single nucleotide polymorphism marker panels. METHODS AND RESULTS: Non-parametric linkage analyses, including 523 families containing a total of 1115 mothers of DZ twins (MODZT) from Australia and New Zealand (ANZ) and The Netherlands (NL), produced four linkage peaks above the threshold for suggestive linkage, including a highly suggestive peak at the extreme telomeric end of chromosome 6 with an exponential logarithm of odds \[(exp)LOD] score of 2.813 (P = 0.0002). Since the DZ twinning rate increases steeply with maternal age independent of genetic effects, we also investigated linkage including only families where at least one MODZT gave birth to her first set of twins before the age of 30. These analyses produced a maximum expLOD score of 2.718 (P = 0.0002), largely due to linkage signal from the ANZ cohort, however, ordered subset analyses indicated this result is most likely a chance finding in the combined dataset. Linkage analyses were also performed for two large DZ twinning families from the USA, one of which produced a peak on chromosome 2 in the region of two potential candidate genes. Sequencing of FSHR and FIGLA, along with INHBB in MODZTs from two large NL families with family specific linkage peaks directly over this gene, revealed a potentially functional variant in the 5' untranslated region of FSHR that segregated with the DZ twinning phenotype in the Utah family. CONCLUSION: Our data provide further evidence for complex inheritance of familial DZ twinning.

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.

Dosage Compensation And Sex Determination In Drosophila - Mechanism Of Measurement Of The X/A Ratio

We propose a molecular mechanism for the intra-cellular measurement of the ratio of the number of X chromosomes to the number of sets of autosomes, a process central to both sex determination and dosage compensation in Drosophila melanogaster. In addition to the two loci, da and Sxl, which have been shown by Cline (Genetics, 90, 683, 1978)and others to be involved in these processes, we postulate two other loci, one autosomal (ω) and the other, X-linked (π). The product of the autosomal locus da stimulates ω and initiates synthesis of a limited quantity of repressor. Sxl and π ,both of which are X-linked, compete for this repressor as well as for RNA polymerase. It is assumed that Sxl has lower affinity than π for repressor as well as polymerase and that the binding of polymerase to one of these sites modulates the binding affinity of the other site for the enzyme. It can be shown that as a result of these postulated interactions transcription from the Sxl site is proportional to the X/A ratio such that the levels of Sxl+ product are low in males, high in females and intermediate in the intersexes. If, as proposed by Cline, the Sxl- product is an inhibitor of X chromosome activity, this would result in dosage compensation. The model leads to the conclusion that high levels of Sxl+ product promote a female phenotype and low levels, a male phenotype. One interesting consequence of the assumptions on which the model is based is that the level of Sxl+ product in the cell, when examined as a function of increasing repressor concentration, first goes up and then decreases, yielding a bell-shaped curve. This feature of the model provides an explanation for some of the remarkable interactions among mutants at the Sxl, da and mle loci and leads to several predictions. The proposed mechanism may also have relevance to certain other problems, such as size regulation during development, which seem to involve measurement of ratios at the cellular level.

Genome-wide association study of migraine implicates a common susceptibility variant on 8q22.1

Migraine is a common episodic neurological disorder, typically presenting with recurrent attacks of severe headache and autonomic dysfunction. Apart from rare monogenic subtypes, no genetic or molecular markers for migraine have been convincingly established. We identified the minor allele of rs1835740 on chromosome 8q22.1 to be associated with migraine (P = 5.38 x 10(-)(9), odds ratio = 1.23, 95% CI 1.150-1.324) in a genome-wide association study of 2,731 migraine cases ascertained from three European headache clinics and 10,747 population-matched controls. The association was replicated in 3,202 cases and 40,062 controls for an overall meta-analysis P value of 1.69 x 10(-)(1)(1) (odds ratio = 1.18, 95% CI 1.127-1.244). rs1835740 is located between MTDH (astrocyte elevated gene 1, also known as AEG-1) and PGCP (encoding plasma glutamate carboxypeptidase). In an expression quantitative trait study in lymphoblastoid cell lines, transcript levels of the MTDH were found to have a significant correlation to rs1835740 (P = 3.96 x 10(-)(5), permuted threshold for genome-wide significance 7.7 x 10(-)(5). To our knowledge, our data establish rs1835740 as the first genetic risk factor for migraine.

Common genetic influences underlie comorbidity of migraine and endometriosis

We examined the co-occurrence of migraine and endometriosis within the largest known collection of families containing multiple women with surgically confirmed endometriosis and in an independent sample of 815 monozygotic and 457 dizygotic female twin pairs. Within the endometriosis families, a significantly increased risk of migrainous headache was observed in women with endometriosis compared to women without endometriosis (odds ratio [OR] 1.57, 95% confidence interval [CI]: 1.12-2.21, P=0.009). Bivariate heritability analyses indicated no evidence for common environmental factors influencing either migraine or endometriosis but significant genetic components for both traits, with heritability estimates of 69 and 49%, respectively. Importantly, a significant additive genetic correlation (r(G) = 0.27, 95% CI: 0.06-0.47) and bivariate heritability (h(2)=0.17, 95% CI: 0.08-0.27) was observed between migraine and endometriosis. Controlling for the personality trait neuroticism made little impact on this association. These results confirm the previously reported comorbidity between migraine and endometriosis and indicate common genetic influences completely explain their co-occurrence within individuals. Given pharmacological treatments for endometriosis typically target hormonal pathways and a number of findings provide support for a relationship between hormonal variations and migraine, hormone-related genes and pathways are highly plausible candidates for both migraine and endometriosis. Therefore, taking into account the status of both migraine and endometriosis may provide a novel opportunity to identify the genes underlying them. Finally, we propose that the analysis of such genetically correlated comorbid traits can increase power to detect genetic risk loci through the use of more specific, homogenous and heritable phenotypes.

Strong evidence for a novel schizophrenia risk locus on chromosome 1p31.1 in homogeneous pedigrees from Tamil Nadu, India

OBJECTIVE: The study of ethnically homogeneous populations may help to identify schizophrenia risk loci. The authors conducted a genomewide linkage scan for schizophrenia in an Indian population. METHOD: Participants were 441 individuals (262 affected probands and siblings) who were recruited primarily from one ethnically homogeneous group, the Tamil Brahmin caste, although individuals from other geographically proximal castes also participated. Genotyping of 124 affected sibling pair pedigrees was performed with 402 short tandem repeat polymorphisms. Linkage analyses were conducted using nonparametric exponential LOD (logarithm of the odds ratio for linkage) scores and parametric heterogeneity LOD scores. Parametric heterogeneity scores were calculated using simple dominant and recessive models, correcting for multiple statistics. The data were examined for evidence of consanguinity. Genomewide significance levels were determined using 10,000 gene dropping simulations. RESULTS: These findings revealed genomewide significant linkage to chromosome 1p31.1, through the use of both exponential and heterogeneity LOD scores, incorporating correction for multiple statistics and mild consanguinity. The estimated sibling recurrence risk associated with this putative locus was 1.95. Analysis for heterogeneity LOD scores also detected suggestive linkage to chromosomes 13q22.1 and 16q12.2. Using 117 tag single nucleotide polymorphisms (SNPs), family-based association analyses of phosphodiesterase 4B (PDE4B), the closest schizophrenia candidate gene, detected no convincing evidence of association, suggesting that the chromosome 1 peak represents a novel risk locus. CONCLUSIONS: This is the first study-to the authors' knowledge-to report significant linkage of schizophrenia to chromosome 1p31.1. Further investigation of this chromosome region in diverse populations is warranted to identify underlying sequence variants.

Susceptibility locus on chromosome 1q23-25 for a schizophrenia subtype resembling deficit schizophrenia identified by latent class analysis

Context: Identifying susceptibility genes for schizophrenia may be complicated by phenotypic heterogeneity, with some evidence suggesting that phenotypic heterogeneity reflects genetic heterogeneity. Objective: To evaluate the heritability and conduct genetic linkage analyses of empirically derived, clinically homogeneous schizophrenia subtypes. Design: Latent class and linkage analysis. Setting: Taiwanese field research centers. Participants: The latent class analysis included 1236 Han Chinese individuals with DSM-IV schizophrenia. These individuals were members of a large affected-sibling-pair sample of schizophrenia (606 ascertained families), original linkage analyses of which detected a maximum logarithm of odds (LOD) of 1.8 (z = 2.88) on chromosome 10q22.3. Main Outcome Measures: Multipoint exponential LOD scores by latent class assignment and parametric heterogeneity LOD scores. Results: Latent class analyses identified 4 classes, with 2 demonstrating familial aggregation. The first (LC2) described a group with severe negative symptoms, disorganization, and pronounced functional impairment, resembling “deficit schizophrenia.” The second (LC3) described a group with minimal functional impairment, mild or absent negative symptoms, and low disorganization. Using the negative/deficit subtype, we detected genome-wide significant linkage to 1q23-25 (LOD = 3.78, empiric genome-wide P = .01). This region was not detected using the DSM-IV schizophrenia diagnosis, but has been strongly implicated in schizophrenia pathogenesis by previous linkage and association studies.Variants in the 1q region may specifically increase risk for a negative/deficit schizophrenia subtype. Alternatively, these results may reflect increased familiality/heritability of the negative class, the presence of multiple 1q schizophrenia risk genes, or a pleiotropic 1q risk locus or loci, with stronger genotype-phenotype correlation with negative/deficit symptoms. Using the second familial latent class, we identified nominally significant linkage to the original 10q peak region. Conclusion: Genetic analyses of heritable, homogeneous phenotypes may improve the power of linkage and association studies of schizophrenia and thus have relevance to the design and analysis of genome-wide association studies.

Molecular genetics of primary open angle glaucoma and exfoliation syndrome

Glaucoma is the second leading cause of blindness worldwide. It is a group of optic neuropathies, characterized by progressive optic nerve degeneration, excavation of the optic disc due to apoptosis of retinal ganglion cells and corresponding visual field defects. Open angle glaucoma (OAG) is a subtype of glaucoma, classified according to the age of onset into juvenile and adult- forms with a cut-off point of 40 years of age. The prevalence of OAG is 1-2% of the population over 40 years and increases with age. During the last decade several candidate loci and three candidate genes, myocilin (MYOC), optineurin (OPTN) and WD40-repeat 36 (WDR36), for OAG have been identified. Exfoliation syndrome (XFS), age, elevated intraocular pressure and genetic predisposition are known risk factors for OAG. XFS is characterized by accumulation of grayish scales of fibrillogranular extracellular material in the anterior segment of the eye. XFS is overall the most common identifiable cause of glaucoma (exfoliation glaucoma, XFG). In the past year, three single nucleotide polymorphisms (SNPs) on the lysyl oxidase like 1 (LOXL1) gene have been associated with XFS and XFG in several populations. This thesis describes the first molecular genetic studies of OAG and XFS/XFG in the Finnish population. The role of the MYOC and OPTN genes and fourteen candidate loci was investigated in eight Finnish glaucoma families. Both candidate genes and loci were excluded in families, further confirming the heterogeneous nature of OAG. To investigate the genetic basis of glaucoma in a large Finnish family with juvenile and adult onset OAG, we analysed the MYOC gene in family members. Glaucoma associated mutation (Thr377Met) was identified in the MYOC gene segregating with the disease in the family. This finding has great significance for the family and encourages investigating the MYOC gene also in other Finnish OAG families. In order to identify the genetic susceptibility loci for XFS, we carried out a genome-wide scan in the extended Finnish XFS family. This scan produced promising candidate locus on chromosomal region 18q12.1-21.33 and several additional putative susceptibility loci for XFS. This locus on chromosome 18 provides a solid starting point for the fine-scale mapping studies, which are needed to identify variants conferring susceptibility to XFS in the region. A case-control and family-based association study and family-based linkage study was performed to evaluate whether SNPs in the LOXL1 gene contain a risk for XFS, XFG or POAG in the Finnish patients. A significant association between the LOXL1 gene SNPs and XFS and XFG was confirmed in the Finnish population. However, no association was detected with POAG. Probably also other genetic and environmental factors are involved in the pathogenesis of XFS and XFG.

Search For Susceptibility Genes In Schizophrenia

Schizophrenia is a severe mental disorder affecting 0.4-1% of the population worldwide. It is characterized by impairments in the perception of reality and by significant social or occupational dysfunction. The disorder is one of the major contributors to the global burden of diseases. Studies of twins, families, and adopted children point to strong genetic components for schizophrenia, but environmental factors also play a role in the pathogenesis of disease. Molecular genetic studies have identified several potential positional candidate genes. The strongest evidence for putative schizophrenia susceptibility loci relates to the genes encoding dysbindin (DTNBP1) and neuregulin (NRG1), but studies lack impressive consistency in the precise genetic regions and alleles implicated. We have studied the role of three potential candidate genes by genotyping 28 single nucleotide polymorphisms in the DNTBP1, NRG1, and AKT1 genes in a large schizophrenia family sample consisting of 441 families with 865 affected individuals from Finland. Our results do not support a major role for these genes in the pathogenesis of schizophrenia in Finland. We have previously identified a region on chromosome 5q21-34 as a susceptibility locus for schizophrenia in a Finnish family sample. Recently, two studies reported association between the γ-aminobutyric acid type A receptor cluster of genes in this region and one study showed suggestive evidence for association with another regional gene encoding clathrin interactor 1 (CLINT1, also called Epsin 4 and ENTH). To further address the significance of these genes under the linkage peak in the Finnish families, we genotyped SNPs of these genes, and observed statistically significant association of variants between GABRG2 and schizophrenia. Furthermore, these variants also seem to affect the functioning of the working memory. Fetal events and obstetric complications are associated with schizophrenia. Rh incompatibility has been implicated as a risk factor for schizophrenia in several epidemiological studies. We conducted a family-based candidate-gene study that assessed the role of maternal-fetal genotype incompatibility at the RhD locus in schizophrenia. There was significant evidence for an RhD maternal-fetal genotype incompatibility, and the risk ratio was estimated at 2.3. This is the first candidate-gene study to explicitly test for and provide evidence of a maternal-fetal genotype incompatibility mechanism in schizophrenia. In conclusion, in this thesis we found evidence that one GABA receptor subunit, GABRG2, is significantly associated with schizophrenia. Furthermore, it also seems to affect to the functioning of the working memory. In addition, an RhD maternal-fetal genotype incompatibility increases the risk of schizophrenia by two-fold.

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.

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.