Genome-wide linkage analysis of multiple measures of neuroticism of 2 large cohorts from Australia and the Netherlands

CONTEXT People meeting diagnostic criteria for anxiety or depressive disorders tend to score high on the personality scale of neuroticism. Studying this personality dimension can give insights into the etiology of these important psychiatric disorders. OBJECTIVES To undertake a comprehensive genome-wide linkage study of neuroticism using large study samples that have been measured multiple times and to compare the results between countries for replication and across time within countries for consistency. DESIGN Genome-wide linkage scan. SETTING Twin individuals and their family members from Australia and the Netherlands. PARTICIPANTS Nineteen thousand six hundred thirty-five sibling pairs completed self-report questionnaires for neuroticism up to 5 times over a period of up to 22 years. Five thousand sixty-nine sibling pairs were genotyped with microsatellite markers. METHODS Nonparametric linkage analyses were conducted in MERLIN-REGRESS for the mean neuroticism scores averaged across time. Additional analyses were conducted for the time-specific measures of neuroticism from each country to investigate consistency of linkage results. RESULTS Three chromosomal regions exceeded empirically derived thresholds for suggestive linkage using mean neuroticism scores: 10p 5 Kosambi cM (cM) (Dutch study sample), 14q 103 cM (Dutch study sample), and 18q 117 cM (combined Australian and Dutch study sample), but only 14q retained significance after correction for multiple testing. These regions all showed evidence for linkage in individual time-specific measures of neuroticism and 1 (18q) showed some evidence for replication between countries. Linkage intervals for these regions all overlap with regions identified in other studies of neuroticism or related traits and/or in studies of anxiety in mice. CONCLUSIONS Our results demonstrate the value of the availability of multiple measures over time and add to the optimism reported in recent reviews for replication of linkage regions for neuroticism. These regions are likely to harbor causal variants for neuroticism and its related psychiatric disorders and can inform prioritization of results from genome-wide association studies.

A high-density association screen of 155 ion transport genes for involvement with common migraine

The clinical overlap between monogenic Familial Hemiplegic Migraine (FHM) and common migraine subtypes, and the fact that all three FHM genes are involved in the transport of ions, suggest that ion transport genes may underlie susceptibility to common forms of migraine. To test this leading hypothesis, we examined common variation in 155 ion transport genes using 5257 single nucleotide polymorphisms (SNPs) in a Finnish sample of 841 unrelated migraine with aura cases and 884 unrelated non-migraine controls. The top signals were then tested for replication in four independent migraine case-control samples from the Netherlands, Germany and Australia, totalling 2835 unrelated migraine cases and 2740 unrelated controls. SNPs within 12 genes (KCNB2, KCNQ3, CLIC5, ATP2C2, CACNA1E, CACNB2, KCNE2, KCNK12, KCNK2, KCNS3, SCN5A and SCN9A) with promising nominal association (0.00041 < P < 0.005) in the Finnish sample were selected for replication. Although no variant remained significant after adjusting for multiple testing nor produced consistent evidence for association across all cohorts, a significant epistatic interaction between KCNB2 SNP rs1431656 (chromosome 8q13.3) and CACNB2 SNP rs7076100 (chromosome 10p12.33) (pointwise P = 0.00002; global P = 0.02) was observed in the Finnish case-control sample. We conclude that common variants of moderate effect size in ion transport genes do not play a major role in susceptibility to common migraine within these European populations, although there is some evidence for epistatic interaction between potassium and calcium channel genes, KCNB2 and CACNB2. Multiple rare variants or trans-regulatory elements of these genes are not ruled out.

Within-family outliers: Segregating alleles or environmental effects? A linkage analysis of height from 5815 sibling pairs

Most information in linkage analysis for quantitative traits comes from pairs of relatives that are phenotypically most discordant or concordant. Confounding this, within-family outliers from non-genetic causes may create false positives and negatives. We investigated the influence of within-family outliers empirically, using one of the largest genome-wide linkage scans for height. The subjects were drawn from Australian twin cohorts consisting of 8447 individuals in 2861 families, providing a total of 5815 possible pairs of siblings in sibships. A variance component linkage analysis was performed, either including or excluding the within-family outliers. Using the entire dataset, the largest LOD scores were on chromosome 15q (LOD 2.3) and 11q (1.5). Excluding within-family outliers increased the LOD score for most regions, but the LOD score on chromosome 15 decreased from 2.3 to 1.2, suggesting that the outliers may create false negatives and false positives, although rare alleles of large effect may also be an explanation. Several regions suggestive of linkage to height were found after removing the outliers, including 1q23.1 (2.0), 3q22.1 (1.9) and 5q32 (2.3). We conclude that the investigation of the effect of within-family outliers, which is usually neglected, should be a standard quality control measure in linkage analysis for complex traits and may reduce the noise for the search of common variants of modest effect size as well as help identify rare variants of large effect and clinical significance. We suggest that the effect of within-family outliers deserves further investigation via theoretical and simulation studies.

Significant evidence of one or more susceptibility loci for endometriosis with near-Mendelian inheritance on chromosome 7p13-15

BACKGROUND: Endometriosis is a common disease with a heritable component. The collaborative International Endogene Study consists of two data sets (Oxford and Australia) comprising 1176 families with multiple affected. The aim was to investigate whether the apparent concentration of cases in a proportion of families could be explained by one or more rare variants with (near-)Mendelian autosomal inheritance. METHODS AND RESULTS: Linkage analyses (aimed at finding chromosomal regions harbouring disease-predisposing genes) were conducted in families with three or more affected (Oxford: n = 52; Australia: n = 196). In the Oxford data set, a non-parametric linkage score (Kong & Cox (K&C) Log of ODds (LOD)) of 3.52 was observed on chromosome 7p (genome-wide significance P = 0.011). A parametric MOD score (equal to maximum LOD maximized over 357 possible inheritance models) of 3.89 was found at 65.72 cM (D7S510) for a dominant model with reduced penetrance. After including the Australian data set, the non-parametric K&C LOD of the combined data set was 1.46 at 57.3 cM; the parametric analysis found an MOD score of 3.30 at D7S484 (empirical significance: P = 0.035) for a recessive model with high penetrance. Critical recombinant analysis narrowed the probable region of linkage down to overlapping 6.4 Mb and 11 Mb intervals containing 48 and 96 genes, respectively. CONCLUSIONS: This is the first report to suggest that there may be one or more high-penetrance susceptibility loci for endometriosis with (near-)Mendelian inheritance.

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.

A case of true hermaphroditism reveals an unusual mechanism of twinning

Traditionally twins are classified as dizygous or fraternal and monozygous or identical (Hall Twinning, 362, 2003 and 735-743). We report a rare case of 46,XX/46,XY twins: Twin A presented with ambiguous genitalia and Twin B was a phenotypically normal male. These twins demonstrate a third, previously unreported mechanism for twinning. The twins underwent initial investigation with 17-hydroxyprogesterone and testosterone levels, pelvic ultrasound and diagnostic laparoscopy. Cytogenetic analysis was performed on peripheral blood cells and skin fibroblasts. Histological examination and Fluorescence in situ hybridization studies on touch imprints were performed on gonadal biopsies. DNA analysis using more than 6,000 DNA markers was performed on skin fibroblast samples from the twins and on peripheral blood samples from both parents. Twin A was determined to be a true hermaphrodite and Twin B an apparently normal male. Both twins had a 46,XX/46,XY chromosome complement in peripheral lymphocytes, skin fibroblasts, and gonadal biopsies. The proportion of XX to XY cells varied between the twins and the tissues evaluated. Most significantly the twins shared 100% of maternal alleles and approximately 50% of paternal alleles in DNA analysis of skin fibroblasts. The twins are chimeric and share a single genetic contribution from their mother but have two genetic contributions from their father thus supporting the existence of a third, previously unreported type of twinning.

KRAS variation and risk of endometriosis

Endometriosis is a common gynaecological disease with symptoms of pelvic pain and infertility which affects 7-10% of women in their reproductive years. Activation of an oncogenic allele of Kirsten rat sarcoma viral oncogene homologue (KRAS) in the reproductive tract of mice resulted in the development of endometriosis. We hypothesized that variation in KRAS may influence risk of endometriosis in humans. Thirty tagSNPs spanning a region of 60.7 kb across the KRAS locus were genotyped using iPLEX chemistry on a MALDI-TOF MassARRAY platform in 959 endometriosis cases and 959 unrelated controls, and data were analysed for association with endometriosis. Genotypes were obtained for most individuals with a mean completion rate of 99.1%. We identified six haplotype blocks across the KRAS locus in our sample. There were no significant differences between cases and controls in the frequencies of individual single-nucleotide polymorphisms (SNPs) or haplotypes. We also developed a rapid method to screen for 11 common KRAS and BRAF mutations on the Sequenom MassARRAY system. The assay detected all mutations previously identified by direct sequencing in a panel of positive controls. No germline variants for KRAS or BRAF were detected. Our results demonstrate that any risk of endometriosis in women because of common variation in KRAS must be very small.

Quantifying northward movement rates of eastern king prawns along eastern Australia

Movement rates of eastern king prawns, Melicertus plebejus (Hess), were estimated from historical and recent conventional tag-recapture information collected across eastern Australia. Data from three studies and 2,656 tag recaptures were used. Recaptured males and females both moved east–north-east in central Queensland and north–north-east in southern Queensland and New South Wales. Over a period of one year, the estimated transition matrix reflected the species strong northerly movement and the more complex longitudinal movement, showing a very high probability of eastern movement in central Queensland and almost negligible eastern or western movement in northern New South Wales. The high exchange probability between New South Wales and Queensland waters indicated that spatial assessment models with movement rates between state jurisdictions would improve the management of this single-unit stock.

Latitudinal and seasonal effects on growth of the Australian eastern king prawn (Melicertus plebejus)

The growth of the Australian eastern king prawn (Melicertus plebejus) is understood in greater detail by quantifying the latitudinal effect. The latitudinal effect is the change in the species’ growth rate during migration. Mark–recapture data (N = 1635, latitude 22.21°S–34.00°S) presents northerly movement of the eastern king prawn, with New South Wales prawns showing substantial average movement of 140 km (standard deviation: 176 km) north. A generalized von Bertalanffy growth model framework is used to incorporate the latitudinal effect together with the canonical seasonal effect. Applying this method to eastern king prawn mark–recapture data guarantees consistent estimates for the latitudinal and seasonal effects. For M. plebejus, it was found that growth rate peaks on 25 and 29 January for males and females, respectively; is at a minimum on 27 and 31 July, respectively; and that the shape parameter, k (per year), changes by –0.0236 and –0.0556 every 1 degree of latitude south increase for males and females, respectively.

The conservation status of the world’s reptiles

Effective and targeted conservation action requires detailed information about species, their distribution, systematics and ecology as well as the distribution of threat processes which affect them. Knowledge of reptilian diversity remains surprisingly disparate, and innovative means of gaining rapid insight into the status of reptiles are needed in order to highlight urgent conservation cases and inform environmental policy with appropriate biodiversity information in a timely manner. We present the first ever global analysis of extinction risk in reptiles, based on a random representative sample of 1500 species (16% of all currently known species). To our knowledge, our results provide the first analysis of the global conservation status and distribution patterns of reptiles and the threats affecting them, highlighting conservation priorities and knowledge gaps which need to be addressed urgently to ensure the continued survival of the world’s reptiles. Nearly one in five reptilian species are threatened with extinction, with another one in five species classed as Data Deficient. The proportion of threatened reptile species is highest in freshwater environments, tropical regions and on oceanic islands, while data deficiency was highest in tropical areas, such as Central Africa and Southeast Asia, and among fossorial reptiles. Our results emphasise the need for research attention to be focussed on tropical areas which are experiencing the most dramatic rates of habitat loss, on fossorial reptiles for which there is a chronic lack of data, and on certain taxa such as snakes for which extinction risk may currently be underestimated due to lack of population information. Conservation actions specifically need to mitigate the effects of human-induced habitat loss and harvesting, which are the predominant threats to reptiles.

Stochastic growth of the Eastern King Prawn (melicertus plebejus (hess, 1865)) harvested off Eastern Australia

Stochastic growth models were fitted to length-increment data of eastern king prawns, Melicertus plebejus (Hess, 1865), tagged across eastern Australia. The estimated growth parameters and growth transition matrix are for each sex representative of the species' geographical distribution. Our study explicitly displays the stochastic nature of prawn growth. Capturing length-increment growth heterogeneity for short-lived exploited species such as prawns that cannot be readily aged is essential for length-based modelling and improved management.

A comparison of DNA pools constructed following whole genome amplification for two-stage SNP genotyping designs

Genotyping in DNA pools reduces the cost and the time required to complete large genotyping projects. The aim of the present study was to evaluate pooling as part of a strategy for fine mapping in regions of significant linkage. Thirty-nine single nucleotide polymorphisms (SNPs) were analyzed in two genomic DNA pools of 384 individuals each and results compared with data after typing all individuals used in the pools. There were no significant differences using data from either 2 or 8 heterozygous individuals to correct frequency estimates for unequal allelic amplification. After correction, the mean difference between estimates from the genomic pool and individual allele frequencies was .033. A major limitation of the use of DNA pools is the time and effort required to carefully adjust the concentration of each individual DNA sample before mixing aliquots. Pools were also constructed by combining DNA after Multiple Displacement Amplification (MDA). The MDA pools gave similar results to pools constructed after careful DNA quantitation (mean difference from individual genotyping .040) and MDA provides a rapid method to generate pools suitable for some applications. Pools provide a rapid and cost-effective screen to eliminate SNPs that are not polymorphic in a test population and can detect minor allele frequencies as low as 1% in the pooled samples. With current levels of accuracy, pooling is best suited to an initial screen in the SNP validation process that can provide high-throughput comparisons between cases and controls to prioritize SNPs for subsequent individual genotyping.

Zygosity diagnosis in the absence of genotypic data: An approach using latent class analysis

For zygosity diagnosis in the absence of genotypic data, or in the recruitment phase of a twin study where only single twins from same-sex pairs are being screened, or to provide a test for sample duplication leading to the false identification of a dizygotic pair as monozygotic, the appropriate analysis of respondents' answers to questions about zygosity is critical. Using data from a young adult Australian twin cohort (N = 2094 complete pairs and 519 singleton twins from same-sex pairs with complete responses to all zygosity items), we show that application of latent class analysis (LCA), fitting a 2-class model, yields results that show good concordance with traditional methods of zygosity diagnosis, but with certain important advantages. These include the ability, in many cases, to assign zygosity with specified probability on the basis of responses of a single informant (advantageous when one zygosity type is being oversampled); and the ability to quantify the probability of misassignment of zygosity, allowing prioritization of cases for genotyping as well as identification of cases of probable laboratory error. Out of 242 twins (from 121 like-sex pairs) where genotypic data were available for zygosity confirmation, only a single case was identified of incorrect zygosity assignment by the latent class algorithm. Zygosity assignment for that single case was identified by the LCA as uncertain (probability of being a monozygotic twin only 76%), and the co-twin's responses clearly identified the pair as dizygotic (probability of being dizygotic 100%). In the absence of genotypic data, or as a safeguard against sample duplication, application of LCA for zygosity assignment or confirmation is strongly recommended.

Quantifying northward movement rates of eastern king prawns along eastern Australia

Movement rates of eastern king prawns, Melicertus plebejus (Hess), were estimated from historical and recent conventional tag-recapture information collected across eastern Australia. Data from three studies and 2,656 tag recaptures were used. Recaptured males and females both moved east-north-east in central Queensland and north-north-east in southern Queensland and New South Wales. Over a period of one year, the estimated transition matrix reflected the species strong northerly movement and the more complex longitudinal movement, showing a very high probability of eastern movement in central Queensland and almost negligible eastern or western movement in northern New South Wales. The high exchange probability between New South Wales and Queensland waters indicated that spatial assessment models with movement rates between state jurisdictions would improve the management of this single-unit stock.